ARTICLE   Open Access    

Neuroprotective comparisons and bioactive profiles of green tea and black tea: in vitro cellular experiments, metabolomics, and network pharmacology analysis

  • # Authors contributed equally: Huan Wu, Juan Wan

More Information
  • This study aimed to compare the neuroprotective effects of green tea and black tea, made from the same raw materials, on an Aβ25-35-induced PC12 cell model, using Ultra-High Performance Liquid Chromatography-Tandem Mass Spectrometry (UPLC-MS/MS) and network pharmacology approaches. This addresses the gap in current research, which has extensively explored the neuroprotective properties of green tea and its components but has paid less attention to black tea. The findings indicate that both teas can alleviate Aβ-induced neurodegenerative changes by reducing inflammation, mitochondrial disruption, and other cellular stressors. Notably, black tea showed higher effectiveness, enriching more differentially expressed genes within critical pathways and exhibiting a broader spectrum of bioactive compounds. Its protein-protein interaction network also suggested that black tea acted on a wider range of potential targets. However, these results are preliminary and emphasize the importance of the complex interplay of bioactive components in tea, advocating for further comparative studies to fully understand their neuroprotective mechanisms.
  • As the world's third pole and Asia's water tower, the Qinghai-Xizang Plateau (QXP) acts as a vital ecological security barrier for the world[1]. In addition, the QXP is also one of the important biodiversity hotspots and harbours many rare resource plants[2]. In recent decades, with global climate change and human activities, the alpine meadow ecosystem on the QXP is facing great risks of grassland degradation and land desertification, which would greatly affect the safety of ecosystems around the world[3]. Therefore, there is an urgent need to repair the desertified lands on the QXP. Among the numerous prevention and control technologies/methods for desertified lands, the construction of stable vegetation adapted to the local climate environment has become the consensus on the sustainable development of desert ecosystems[4]. Due to the long-term stresses of low oxygen, strong ultraviolet radiation, drought, and the harsh conditions of alpine environments, plants on the QXP typically grow very slowly. Therefore, once vegetation degradation occurs on the QXP, restoration efforts become considerably more challenging. Compared to other plant species, desert-adapted plant species exhibit greater resilient to the harsh environmental conditions of the QXP, making them ideal for restoring degraded vegetation. Therefore, plant species naturally adapted to deserts on the QXP are the optimal choice for vegetation reconstruction in alpine desertified lands.

    Over the past decades, numerous researches have been carried out to dissect the genetic basis of plateau adaptation in species endemic to the extreme environments on the QXP[5,6]. Among these studies, most were focused on human and animal species[79], while studies on plant adaptation to the QXP have started to gain more attentions in recent years, driven by an increased recognition of the importance of plant diversity for ecosystem resilience[1012]. Based on comparative transcriptome, metabolome, and/or genome analysis, an increasing number of studies have substantiated that most of the plants thriving on the QXP possess an abundance of secondary metabolites and robust genetic resources tailored to withstand its severe natural conditions[13,14]. Most studies have focused on sister species within genera, however, due to the long history of species differentiation, population dynamics, and/or breeding systems vary significantly across species, making it rare to investigate the genetic basis of altitude gradient adaptation within the same plant species. Nevertheless, plants native to the deserts of the QXP exhibit remarkable tolerance to multiple stresses such as UV-B radiation, drought, cold, and hypoxia, rendering them optimal pioneer species for restoring desertified lands in this region. Furthermore, a clarified molecular basis of local adaptation in different ecotypes within a plant species can not only forecast the evolutionary trajectory of plant adaptation to future climate changes[15,16], but also provide crucial molecular insights and genetic resources for improving and selecting plant species capable of surviving extreme environments and severe climate fluctuations on the QXP[17]. Meanwhile, due to long-term adaptation to the harsh environment and special climatic conditions of the QXP, indigenous plants have produced many secondary metabolites with medicinal value during their adaptive evolution. Therefore, multi-omics studies among natural populations of plant species are necessary to unravel the molecular metabolic mechanisms underlying the adaptation of desert plants to altitude gradients on the QXP, which can not only reconcile the conflicts between local agricultural development and ecological conservation in these fragile ecosystems, but also provide valuable insights into the indigenous nature of high-altitude medicinal plants.

    Agriophyllum squarrosum (L.) Moq., commonly known as sandrice, thrives in the vast deserts and sandy landscapes of arid and semi-arid regions throughout the interior of Asia (www.efloras.org). Field investigations underscore its remarkable ecological versatility, thriving at altitudes from about 50 to over 4,000 m, particularly adapting well to the harsh desert conditions of the QXP, where it exhibits strong growth and reproductive success[18,19]. Additionally, sandrice plays a crucial ecological role in reducing wind velocity by at least 90% when withered, and it enriches nutrient-poor soils with carbon and nitrogen, significantly sustaining and restoring fragile desert ecosystems[18,20]. Moreover, despite growing in the infertile sandy soils, sandrice seeds are rich in essential nutrients like amino acids, crude fiber, and polyunsaturated fatty acids, making it an ideal natural functional food[21]. Furthermore, its above-ground parts are abundant in bio-actives, including flavonoids, organic acids, terpenoids, and alkaloids[22], and have been traditionally used in Mongolian medicine for treating kidney inflammation, dyspepsia, fever, and pain relief[23]. Notably, recent studies have indicated that sandrice shows great potential as both an antibiotic substitute and a functional forage crop in antibiotic-free ruminant farming, owing to the abundance of bio-active compounds found in its aerial parts[24,25]. Therefore, exploring the mechanism of alpine adaptation in sandrice would not only help combat desertification in the plateau region, but also enhance our understanding of the factors contributing to the high medicinal quality of its alpine ecotypes.

    Previous biogeographic studies have underscored notable genetic divergence among sandrice populations across heterogeneous deserts and sandy lands. However, minimal genetic divergence was observed between the alpine group and its adjacent central desert group, despite notable habitat heterogeneity between them[19,26]. Notably, based on metabolomic analysis and common garden experiments, variations were observed in the accumulation of medicinal metabolites with significant pharmacological activity, such as flavonoids, among populations from different altitudinal habitats, even under the same environmental conditions.[27,28]. Recently, cold-stress treatment among ecotypes from different altitude gradients further indicates that flavonoids are crucial for sandrice to defend against low temperatures[29]. These phenomena suggest that flavonoid biosynthesis pathways in sandrice may have been favored through long-term local adaptation to the QXP, and further contributes to its distinctive medicinal properties.

    To investigate the apparent paradox of significant differences in secondary metabolite accumulation despite similar genetic backgrounds in sandrice, and to further verify the role of flavonoid biosynthesis pathway in the alpine adaptation of sandrice, this study conducted the first in situ metabolome and transcriptome analyses comparing two ecotypes of sandrice. These ecotypes occupy habitats at different altitudes while sharing the same genetic composition: ETL from the alpine group (2,917 m altitude) and CC from the central desert group (1,530 m altitude). Then, population genetic analysis across 22 natural populations was carried out to determine whether these adaptive functional genes underwent directional evolution along the altitude gradient, compared to neutral genes. This endeavour aims to elucidate the molecular metabolic basis underlying sandrice's adpatation to harsh alpine desert environments, especially the role of metabolites with medicinal value in plant adaptive evolution, which will further provide molecular guidance and genetic resources for the restoration of desertification on the QXP and facilitate molecular marker-assisted breeding to enhance the medicinal quality of this promising plant sepcies.

    Based on the neutral genetic markers, two wild ecotypes that shared the same genetic composition were collected at the same mature growth period from different altitudinal habitats. One was located in Ertala (ETL), Gonghe county, Qinghai province (36°11'39.48'' N, 100°31'28.26'' E, 2,917 m) on the QXP, while the other one was located in Changcheng county (CC), Gansu province (37°54'10.98'' N, 102°54'4.2'' E, 1,530 m) in the southern edge of the Tengger Desert (Fig. 1a). To minimize variations caused by different growth stages, the ETL ecotype was collected in mid-September 2016, while the CC ecotype was collected in early October. Both the CC and ETL ecotypes were in the early reproductive stage at the time of collection. All collected samples were accurately identified by Prof. Xiao-Fei Ma, Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences. Fresh tissues from each ecotype, including leaves, stems and spikes, were promptly flash-frozen in liquid nitrogen and stored at −80 °C in an ultra-low temperature freezer for further extraction of metabolites and RNA.

    Figure 1.  Comparative metabolomic profilings of two altitudinal ecotypes of A. squarrosum. (a) Geographical localization of two altitudinal ecotypes of A. squarrosum provenances. (b) PCA plot of two ecotypes metabolomes revealing the first two principal components with t[1] = 30.5% and t[2] = 16.7%. The meaning of each abbreviation is as follows: L (leaf), St. (stem), Sp. (spike). Each samples have three biological replicates. (c) Heat map showing standardized contents of partial classified metabolites in different tissues of ETL and CC ecotype. These metabolites are classed into two groups according to their chemical characters.

    As the non-targeted approach provides the advantage of discovering important metabolites that might otherwise remain undetected with a targeted approach, to identify the key metabolites and metabolic pathways, particularly the major secondary metabolic pathways, involved in plateau adaptation for sandrice, samples with three biological replicates of each ecotype were prepared for UPLC-MS non-targeted metabonomics using LC-ESI-Q TRAP-MS/MS systems at Metware Biotechnology Co., Ltd (Wuhan, China). The quantification of the metabolites was carried out in MRM mode and the analytical conditions were as the study of Chen et al.[30]. Analyst 1.6.1 and MultiQuant 3.0.2 software were used for data set acquisition, peak recognition, and normalization. Metabolites were annotated by mapping them to the self-built database MWDB (Metware Biotechnology Co., Ltd. Wuhan, China) as well as public databases to identify their chemical structures. Quantitative analysis of metabolites was carried out by a multi-reaction monitoring mode (MRM) on a triple quadrupole mass spectrometer.

    To further determine the differentially enriched metabolites (DEMs) between the two ecotypes, PCA (Principal Component Analysis) and PLS-DA (PLS Discriminant Analysis) were performed with SIMCA 13.0 software. DEMs were determined based on relative content, with thresholds set at a variable importance in the projection (VIP) value of ≥ 1 and a fold change of ≥ 2 or ≤ 0.5.

    Total RNA was extracted from each tissue using RNAprep Pure Plant Kit (Tiangen Biotech Co., Ltd, Bejing, China). Double-strand cDNA was constructed following the study of Shi et al.[31]. These fragments were firstly purified with QiaQuick PCR extraction kit (QIAGEN Inc., Valencia, CA, USA) according to the construction, then were end-repaired with A added to the 3' ends, and finally ligated to sequencing adaptors. The ligated cDNA products under size-selected demands were further concentrated by PCR amplification to construct the cDNA libraries. Library preparations were sequenced on an Illumina HiseqTM3000 platform with a 150-bp paired-end mode. Raw sequenced data have been submitted to the NCBI database with the accession number PRJNA659807.

    Raw reads containing unknown sequences ('N') exceeding 5% and low-quality reads (with a base quality less than Q20) were eliminated from the dataset. The remaining filtered clean reads were then utilized for de novo assembly with Trinity version 2.4.0[32]. According to the pair-end information, contigs were clustered and assembled into sequences as long as possible after removing redundancies, and then the clustered longest contigs were subsequently amalgamated into the total unigenes. Functional annotation of each unigene was performing BLASTx searches against the public protein and/or nucleotide databases (such as the NCBI Nr, Nt databases, Swiss-Prot protein database, KOG database, GO database, InterPro, and the KEGG database) with an E-value cutoff of 1e-5. Differentially expressed genes (DEGs) between different tissues of the two ecotypes were estimated by DESeq2. A significance threshold was set with a p-value less than 0.05 and an absolute value of fold-change over 2 to determine significant differential expression[33]. Enrichment analysis of DEGs in the Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways was performed using KOBAS software (version 2.0.12) and visualized with ggplot2.

    To explore the potential regulatory relationships between genes and metabolites, the relative expression of DEGs and the relative contents of DEMs within two interested pathways of phenylpropanoid and flavonoids biosynthesis pathways were firstly standardized respectively using z-score standardization. Then, Pearson correlation coefficients between DEGs and DEMs were calculated by R version 4.2.3. Correlation pairs were selected with an absolute value threshold over 0.9 and p-value lower than 0.05. Finally, the gene-metabolite co-expression network was visualized with Cytoscape-v3.7.2 (available at www.gnu.org/licenses/lgpl-2.1.html).

    To explore orthologous genes that have potentially undergoing adaptive differentiation between two ecotypes in response to differing altitudinal environments, the transcriptome sequences of each ecotype were first assembled separately by Trinity version 2.4.0[32]. Subsequently, the resulting clean sequences were searched by BLASTp version 2.2.31 under the threshold of E-value < 1e-5, and then were predicted and translated into protein-coding sequences by TransDecoder version 3.0.0. Meanwhile, OrthoMCL-V2.0.9 software was employed to discern potential orthologs and paralogs among the protein sequences derived from each ecotype's transcriptome. Then, pairwise comparisons were further conducted on putative single-copy orthologs to estimate selective pressure, and ParaAT was used to parallelly construct protein-CDS alignments for these orthologs[34]. Finally, the synonymous substitution rates (Ks), non-synonymous rates (Ka) and Ka/Ks value were calculated for each putative single-copy homologous gene between the two ecotypes with KaKs_Calculator 2.0[35] under the YN model of approximate method[36]. Pairs with Ks > 0.1 were excluded to avoid potential paralogs. The positive selection genes (PSGs) with Ka/Ks value higher than 1 were further verified by the codeml program in PAML[37].

    To further verify whether sandrice has significantly diverged among ecotypes along with different altitude gradients, population analysis was also conducted among populations of sandrice inhabiting different altitudes. A total of 22 sandrice populations with four to six individuals for each population were collected, including the alpine group from the Qaidam Basin with an altitude of 2,000−3,500 m, the middle altitude group from Tengger desert, Ulan Buhe Desert, Kubuqi Desert, and Mu Us sandy land with an altitude of 1,000−2,000 m, and the low altitude group from Otindag sandy land, Horqin sandy land, and Hulun Buir sandy land with altitude of 0−1,000 m (Supplementary Table S1). Meanwhile, a population of A. minus was also collected from Gurbantunggut desert as the outgroup. All the samples were selected apart from > 50 m for each individual. Fresh leaves were dried and preserved in silica gel, and voucher specimens were deposited in the Northwest Institute of Eco-Environment and Resources, Chinese Academy of Sciences.

    Total genomic DNA was isolated from the tissue samples with TIANGEN Plant Genomic DNA Kit (Tiangen Biotech Co., Ltd, Bejing, China) following the manufacturer's protocol. All the DNA samples were quantified by Qubit assay HS kit (Life Technologies, Burlington, ON, Canada) with assays read on a Qubit v2.0 (Life Technologies). In total, 24 pairs of primers were designed by PRIMER version 5.0 based on the RNA-seq data (Supplementary Table S2). Among them, genes under positive selection identified by selective pressure analysis, especially those associated with the known stress-resistance pathways of phenylpropanoid and flavonoid biosynthesis pathways[13,3841], were identified as candidate genes. Conversely, neutral genes annotated with irrelevant functions to stress resistance were considered reference genes (Supplementary Table S2). All these gene fragments were amplified across these 22 populations using 2 × Taq Plus highfidelity PCR MasterMix (Tiangen, Beijing, China) in a Gene-Amp PCR system 9700 DNA Thermal Cycler (PE Applied Biosystems, Norwalk, USA) following the programs listed in Supplementary Table S2. PCR products were purified with TIAN quick Midi Purification Kits (Tiangen, Beijing, China) and then were Sanger-sequenced with both forward and reverse primers by GENEWIZ company (Tianjin, China). Multiple sequences were aligned and adjusted manually by BIOEDIT version 7.2.6.1 software[42]. The structures of all gene fragments were defined by alignment with their corresponding mRNA sequences and their best hits of BLAST on ESTs (Expressed Sequence tags) from NCBI. All these new sequences were deposited in GenBank under accession numbers OM338032-OM338057, OP846852-OP846955.

    Genetic diversity was estimated for each gene fragment in three groups with different altitudinal gradients by calculating the number of segregating sites (S), nucleotide diversity statistics (θw; π)[43,44], the number of haplotypes (Nh) and haplotype diversity (He)[45] for all sites, silent sites, and nonsynonymous sites with DnaSP version 5.10[46]. Meanwhile, the fixation index (FST) of each loci was also computed among high, middle, and low altitude groups to estimate the genetic divergence degree with AMOVA in the program Arlequin version 3.1.1 with default settings[47].

    Furthermore, to detect whether there were any signals of evolutionary adaptation to different altitudinal gradients habitats, neutrality tests were performed for each fragment with several methods, including Tajima's D test , Fu & Li's D* and F* test[43], Fay & Wu's H test[48], DH test[49], McDonald-Kreitman (MK) test[50] and the maximum frequency of derived mutations (MFDM) test[51]. Finally, to understand the evolutionary relationships and patterns of these putative alpine adaptive genes across different altitude populations, genealogical topologies were constructed using the median-joining (MJ) model in NETWORK Version 4.6.1.259[52] for their haplotypes.

    Based on the non-targeted UPLC-MS/MS metabolic profiling, a total of 506 metabolites were detected. Among these metabolites, 244 metabolites could be matched to known biochemical structures (Supplementary Fig. S1), including 39 flavonoids, 39 amino acids and derivatives, 27 nucleotide and derivatives, 19 polyphenols, 18 vitamins, 16 alkaloids, 11 lipids, seven organic acids, eight coumarins, eight terpnoids, nine carbohydrates, and 43 additional compounds. PCA analysis revealed distinct clustering of metabolites, segregating into ETL and CC groups based on two altitudinal ecotypes and different tissue types (Fig. 1b). Remarkably, metabolite profiles in leaves differed from those in stems and spikes. However, metabolites in stem samples exhibited similarities or equivalences to those in spikes.

    Among the 244 metabolites, the most prevalent secondary metabolites were identified as flavonoids, alkaloids, and polyphenols. In comparison with CC, ETL exhibited higher levels of total hesperetin, quercetin, betaine, trigonelline, caffeic acid, and ferulic acid. Conversely, CC displayed greater accumulation of total tricin, chrysoeriol, etamiphylline, theobromine, sinapic acid, and p-coumaric acid (Fig. 1c). Among three tissues of the two ecotypes, the most remarkable and largest number of DEMs were identified as flavonoid and polyphenol compounds, followed by nucleotides and derivatives, as well as lipids. In ETL, the contents of eriodictyol chalcone and ferulic acid O-hexoside significantly accumulated, whereas tricin 7-O-rutinoside showed a notable accumulation in CC. Compared to CC, the leaf of ETL exhibited significant increases in apigenin 7-O-rutinoside and quercetin-3-beta-O-galactoside, while the content of chrysoeriol O-hexoside decreased significantly. In the stem and spike of ETL, hesperetin 5-O-glucoside content was significantly higher, whereas two glycosylated tricins (tricin O-rutinoside and tricin 5-O-hexoside) were significantly reduced. Additionally, kaempferide showed a significant decrease in levels specifically in the stem of ETL (Table 1).

    Table 1.  List of metabolites in phenylpropanoid and flavonoid biosynthesis pathways in two altitudinal ecotypes of A. squarrosum.
    Class Metabolite name ETL vs CC
    Leaf Stem Spike
    FC VIP FC VIP FC VIP
    Flavonoid Naringenin 7-O-glucoside 0.98 0.01 0.61 0.56 0.52 0.56
    Hesperetin 5-O-glucoside 1.30 0.27 2.00 1.22 2.38 1.22
    Apigenin 7-O-rutinoside 105.03 1.59 2.44 0.66 2.18 0.66
    Luteolin 7-O-glucoside 1.44 0.67 1.31 0.43 0.93 0.43
    Luteolin O-hexoside 1.77 0.96 1.33 0.78 0.92 0.78
    Chrysoeriol O-hexoside 0.39↓ 1.52 0.45 0.94 0.65 0.94
    Chrysoeriol C-hexoside 0.42 0.57 0.33 0.72 0.26 0.72
    Selgin O-hexoside 1.03 0.54 1.16 0.53 1.07 0.53
    Tricin O-rutinoside 0.32 0.89 0.17↓ 1.02 0.32↓ 1.02
    Tricin 7-O-rutinoside 0.34↓ 1.09 0.20↓ 1.27 0.40↓ 1.27
    Tricin 5-O-hexoside 0.83 0.12 0.43↓ 1.19 0.38↓ 1.19
    Eriodictyol chalcone 5.93 1.39 3.43 1.42 2.62 1.42
    Quercetin-3-beta-O-galactoside 4.23 1.03 1.19 0.08 1.07 0.08
    Kaempferol 3-O-glucoside 1.55 0.96 1.24 0.38 0.95 0.38
    Kaempferide 1.23 0.69 0.48↓ 1.42 0.62 1.42
    Phenylpropanoids p-Coumaric acid 0.83 1.22 0.69 0.88 0.79 0.88
    Caffeic acid 1.74 1.00 1.00 0.22 2.28 0.22
    Ferulic acid O-hexoside 3.46 1.26 4.20 1.48 2.75 1.48
    The relative abundance of metabolites were displayed in Supplementary Table S1. FC, the fold change of metabolites when ETL compared CC; VIP, variable importance in the projection value. Bold font indicates significantly changed metabolites (FC ≥ 2 or ≤ 0.5,VIP ≥ 0.5). represents increased metabolites in ETL. ↓ represents decreased metabolites in ETL.
     | Show Table
    DownLoad: CSV

    After the removal of sequences with low quality, poly-N, and adaptors, 884,051,398 clean reads were filtered from the 128.96 Gb raw data, with Q30 over 81.13% and the average N50 of 690.38 bp (Supplementary Table S3). Subsequently, these high-quality trimmed clean reads were de novo assembled into contigs and then a joint transcript of 135,686 unigenes. Finally, 69,977 annotated unigenes were identified across all transcripts, each represented in at least one database (Supplementary Table S4). The most abundant KOG terms identified in the unigenes included those related to general function prediction only, signal transduction mechanisms, post-translational modification, protein turnover, and chaperones (Supplementary Fig. S2a). The GO enrichment analysis of the entire transcriptome revealed that the annotated unigenes predominantly participated in metabolic processes and cellular processes (Supplementary Fig. S2b).

    In comparison to the gene expression levels in CC, 16,680 unigenes were up-regulated and 12,773 unigenes were down-regulated in ETL. Additionally, Venn analysis of DEGs revealed that the number of tissue-specific DEGs was higher in the leaf and stem than in the spike (Fig. 2a). The highest number of DEGs with the strongest differential expression level was detected in leaves. KEGG analysis revealed that DEGs of three tissues between the high and middle ecotypes of sandrice were predominantly enriched in pathways such as photosynthesis, starch and sucrose metabolism, flavonoid synthesis, phenylpropanoid synthesis, ribosomal regulation, and carbon metabolism. (Fig. 2bd).

    Figure 2.  The state and KEGG pathway enrichment of differentially expressed genes (DEGs) in different tissues from high- and middle-ecotypes of A. squarrosum. (a) Venn diagram indicating the overlapping and unique up-regulated (left) and down-regulated (right). (b)−(d) KEGG pathway enrichment analysis of DEGs in leaf, stem and spike. The number of genes is indicated by the size of the circle and the color of the circle shows significant enrichment through P-value. The top 20 pathways with the minimum P-value are shown in each tissue.

    To explore candidate genes involved in high-altitude adaptation, comparative transcriptome analysis was further conducted. A total of 10,275 pairs of single-copy putative orthologous genes were identified after filtering out those with unexpected stop codons. Among them, 127 pairs of orthologous genes showed signs of positive selection with Ka/Ks values greater than 1 (Supplementary Table S5). Although no significantly over-represented KEGG categories were detected, and half of these positively selected genes (PSGs) could not be well-matched to several annotation databases, some PSGs were still putatively annotated to functions related to DNA repair, response to stress, and metabolism (Table 2).

    Table 2.  Typical differentiated genes for high elevation adaptation.
    Unigene ID ETL_Leaf
    FPKM
    CC_leaf
    FPKM
    ETL_spike
    FPKM
    CC_spike
    FPKM
    ETL_stem
    FPKM
    CC_stem
    FPKM
    La Ka Ks Ka/Ks Gene annotation Possible functions and biological process
    TR86494|c0_g1_i1 0 0 1.71 2.76 0.20 0.34 1044 0.00787 0.00365 2.157 Shikimate O-hydroxycinnamoyltransferase (HCT) Phenylpropanoid biosynthesis;
    Flavonoid biosynthesis
    TR80792|c3_g1_i2 51.14 25.12 6.13 1.94 29.94 4.86 1176 0.00780 0.00423 1.845 Flavonol synthase (FLS) Phenylpropanoid biosynthesis;
    Flavonoid biosynthesis
    TR45359|c0_g2_i1 29.58 5.35 1.73 0.27 7.19 1.72 489 0.0164 0.00899 1.828 Caffeoyl-CoA O-methyltransferase (CCoAOMT) Phenylpropanoid biosynthesis;
    Flavonoid biosynthesis
    TR934|c0_g1_i1 0.94 0.84 2.68 2.48 2.49 1.38 1143 0.00288 0.00272 1.056 A/G-specific adenine DNA glycosylase (ANG ) Base excision repair
    TR74739|c0_g2_i1 2.13 3.45 3.52 3.52 3.50 2.85 1791 0.00514 0.00318 1.615 Uracil-DNA glycosylase (UNG) Base excision repair; Immune diseases
    TR39350|c0_g4_i1 1.59 3.59 1.53 2.28 9.18 6.21 1434 0.00358 0.00330 1.087 Vegetative cell wall protein gp1(GP1) Defense response
    TR39755|c0_g1_i1 7.79 2.74 5.65 4.64 7.37 3.86 1089 0.00505 0.00345 1.463 Elongator complex protein 4 (ELP4) Response to oxidative stress, abscisic acid signaling pathway
    TR40190|c7_g5_i1 25.76 28.38 23.06 20.04 31.05 27.46 300 0.02602 0.00831 3.130 Calcium-dependent protein kinase 1 (CDPK1) Plant-pathogen interaction; Response to drought, cold stress; Salicylic acid biosynthesis
    FPKM, gene expression level; La, the length of candidate genes (bp); Ka, nonsynonymous substitution rate; Ks, synonymous substitution rate; Ka/Ks, selective strength; p-value, the value computed by Fisher exact test when calculating Ka/Ks
     | Show Table
    DownLoad: CSV

    For example, two genes (TR934|c0_g1_i1, TR74739|c0_g2_i1), which encoded A/G-specific adenine DNA glycosylase (ANG) and uracil-DNA glycosylase (UNG) participated in DNA base excision repair; TR39350|c0_g4_i1 encoding Vegetative cell wall protein gp1 (GP1) took part in defense response; TR39755|c0_g1_i1 encoding a subunit of elongator complex (ELP4), mediated ABA signaling pathway and manifested oxidative stress resistance. TR40190|c7_g5_i1, which encoded putative calcium-dependent protein kinase family protein (CDPK1), played a vital role in pathogen resistance abiotic stress and salicylic acid biosynthesis. Besides, three genes encoding shikimate O-hydroxycinnamoyltransferase (HCT), flavonol synthase (FLS), and caffeoyl-CoA O-methyltransferase (CCoAOMT), which are involved in phenylpropanoid and flavonoid pathway, were suggested to be under strong positive selection between the two altitudinal ecotypes (Table 2).

    The integrated analysis between DEMs and DEGs in the phenylpropanoid and flavonoid pathway revealed a significant relationship between the accumulation of metabolisms and the expression of key genes. For example, PAL, CHS, CHI, F3H, and FLS showed significant up-regulation in the high-altitude ecotype ETL compared to CC, consisting of significantly elevated enrichment of corresponding downstream flavonoid and phenylpropanoid contents in ETL (Fig. 3a). Subsequently, as illustrated in Fig. 3b, correlation analysis between gene expression in flavonoid and phenylpropane metabolic pathways and the enrichment of metabolites further revealed that the expression of the COMT1 is positively correlated with the accumulation of most metabolites in the flavonoid and phenylpropanoid metabolic pathways, except for luteolin-O-hexoside and quercetin 3-O-glucoside, which exhibit negative correlations. In terms of differential accumulation of metabolites between the two altitude ecotypes, the level of naringenin 5-O-glucoside directly correlated with the differential expression of CHI, CHS, and F3H. Meanwhile, the accumulation of chrysin 5-O-hexoside was positively correlated with the expression of CHS and COMT. The high expression of these genes would reduce the availability of the substrate naringenin for chrysin synthesis, thereby prompting the production of products from alternative pathways (Fig. 3a). Differences in the accumulation of ferulic acid O-hexoside correlated directly with the differential expression of PAL. Furthermore, the content of quercetin 3-O-glucoside was negatively correlated with the expression levels of CHI, CHS, 4CL, and F3H genes in the high-altitude ecotype of sandrice, leading to reduced quercetin synthesis compared to the mid-altitude ecotype.

    Figure 3.  Phenylpropanoid and flavonoid pathways in A. squarrosum. (a) The schematic representation of gene and metabolite changes in phenylpropanoid and flavonoid pathways. The dotted line represented unreported pathway. (b) The interacted network constituted from genes and metabolites co-expression in phenylpropanoid and flavonoid pathways.

    Population genetics analysis revealed that putative alpine adaptive genes exhibited higher levels of nucleotide diversity compared to each reference locus, as well as the haplotype diversity (Supplementary Table S6). Interestingly, further analysis of genetic diversity across different altitude gradient groups revealed discernible differences, despite the statistical insignificance (Supplementary Fig. S3). Specifically, gene segments at higher altitudes exhibited notably lower nucleotide diversity. Notably, certain gene segments in high-altitude populations, such as PAL, C3H, FOMT, FNS, F3H, CYP75B4, and CHS, displayed nucleotide diversity as low as zero (Supplementary Fig. S3).

    Fixation index (FST) values were also estimated to assess genetic differentiation among high, middle, and low altitude populations with population genetic data from candidate genes and reference loci, supplemented by SNPs obtained through restriction site-associated DNA sequencing (RAD-seq), which provided a representation of genome-wide variation. Significant genetic differentiation was observed among the high-, middle-, and low-altitude populations (Fig. 4). Specifically, between the high and middle altitude populations, several candidate genes showed notable levels of FST: UNG (FST = 0.410), CDPK1 (FST = 0.228), FLS (FST = 0.144), GP1 (FST = 0.104), and CCoAOMT (FST = 0.099), which were all higher than the average genome-wide FST level (FST = 0.051), indicating significant genetic divergence. Within the phenylpropanoid and flavonoids pathway, F3'H (FST = 0.525), and CHI (FST = 0.300) exhibited the highest levels of genetic differentiation between the high and middle populations.

    Figure 4.  Genetic differentiation among different altitudinal populations in A. squarrosum. (a) High- and middle-altitude populations. (b) High- and low-altitude populations. (c) Middle- and low-altitude populations. The red dashed line represents the average genome-wide FST level.

    Neutrality tests were further conducted for all candidate genes and reference loci among 22 populations of sandrice along with altitudinal clines. Among these, four PSGs (ELP4, GP1, FLS, and HCT) were fixed in the high-altitude group with only one allele, as well as nine genes involved in phenylpropanoid and flavonoids pathways, such as FLS, HCT, PAL, C3H, FOMT, FNS, F3H, CYP75B4, and CHS (Supplementary Table S6). Furthermore, as indicated by Tajima's D, Fu & Li's D*, F* values, UNG showed robust signal of positive selection in the high-altitude populations (Supplementary Table S7). Interestingly, for PSG CCoAOMT, involved in the phenylpropanoid and flavonoids pathway, Tajima's D and F* were significantly greater than zero. The MK test of CCoAOMT was also significant, suggesting that it was under balancing selection at the population level. However, the previous Ka/Ks value indicated positive selection for this gene (Table 2). Furthermore, haplotype network and topology analysis of CCoAOMT across different altitude populations showed that H1 might be the ancestral haplotype, while haplotypes H3 and H4 were specific to the high-altitude populations (Fig. 5a). Combined with distinct expression patterns of its alleles in high and middle altitude ecotypes, as revealed by transcriptome data featuring a non-synonymous mutation (Fig. 5b), CCoAOMT appears to have been selected due to ecological differentiation.

    Figure 5.  Haplotype distribution of CCoAOMT in A. squarrosum populations and the expression of two alleles of this gene based on transcriptome data. (a) Haplotype topology of CCoAOMT. (b) Expression levels and mutation sites of the two CCoAOMT alleles in high-altitude and middle-altitude ecotypes of A. squarrosum.

    The QXP, adjacent to arid Central Asia, offers diverse habitats for biomes, however, it is highly sensitive to ongoing climate changes, particularly extensive desertification of alpine meadows[3]. Vegetation colonizing these desertified areas has evolved adaptive traits to cope with extreme environmental conditions resulting from climate change scenarios, and are also a priority candidate for vegetation reconstruction in alpine desertified lands[17]. However, few studies have explored the molecular basis of adaptation in plant species native to alpine desert ecosystems, especially concerning ecotypes across altitude gradients. As a pioneering species in vegetation restoration of desertified lands, sandrice provides compelling evidence that long-term local adaptation to multiple stresses drives adaptive divergence. This adaptation could significantly impact the success of ecological restoration and development in alpine grassland ecosystems threatened by desertification. Meanwhile, it also provides a solid foundation for improving and developing the medicinal value of sandrice.

    Metabolites, particularly secondary metabolites, are pivotal for shaping species-specific traits and are integral to how plants respond to challenging environmental conditions[53]. Previous research indicates that in adapting to alpine conditions, plants have evolved to produce a rich array of secondary metabolites, including phenylpropanoids, flavonoids, ascorbic acid, etc. Besides the outstanding pharmaceutical values, these compounds are also pivotal in helping plants resist environmental stressors, providing numerous advantages such as antioxidative properties, the scavenging of reactive oxygen species (ROS), UV-B radiation absorption, enhanced cold tolerance, and the maintenance of osmotic balance[54]. Phenolics such as flavonoids and phenylpropanoids have even been demonstrated to exhibit species-specific distribution patterns that accumulate along altitudinal gradients in certain plant species[13,41,55]. In this study, compared to the mid-altitude ecotype of CC, significant enrichments of secondary metabolites were observed in the alpine ecotype of ETL, particularly in the phenylpropanoid and flavonoid biosynthesis pathways, such as hesperetin, betaine, quercetin, apigenin, caffeic acid, ferulic acid, etc (Fig. 2c). Interestingly, among these DEMs, apigenin, nobiletin, and ferulic acid were primarily found in glycoside forms (Table 1), which demonstrated potent antioxidant activity by effectively scavenging ROS to safeguard cellular functions and biotic stressor resistance[55,56]. Significant accumulations of flavonoid glycosides have also been found in qingke (Hordeum vulgare L. var. nudum), a crop that has been cultivated and exposed to long-term and intense UV-B radiation on the QXP[10]. Previous studies based on common garden experiments have demonstrated flavonoid accumulation in sandrice is positively correlated with temperature and UV-B radiation, but negatively affected by precipitation and sunshine duration[27,28]. Thus, for the alpine ecotype of sandrice, the significant accumulation of metabolisms, especially the flavonoids, and phenylpropanoids, would hypothetically be induced for adaptation to the long-term harsh alpine desert environment stressors, which further contribute to the high medicinal quality of alpine ecotype.

    Besides detecting a substantial number of DEMs within the phenylpropanoid and flavonoid biosynthesis pathways between mid-altitude and alpine ecotypes, a significant clustering of DEGs associated with the phenylpropanoid and flavonoid biosynthesis pathways was also observed between these two ecotypes (Fig. 2). Correlation analysis of DEMs and DEGs further demonstrated that the differential expression of constitutive genes in the phenylpropanoid and flavonoid biosynthesis pathways would contribute to the divergent enrichments of the downstream DEMs in sandrice. For instance, it has been observed that the transcriptional up-regulation of flavonoid biosynthesis genes (such as PAL, CHS, CHI, FLS, F3H, etc.) significantly enhances the accumulation of hesperetin, apigenin, and quercetin (Fig. 3a, b). These metabolites have been documented to boost plant resistance against UV-B radiation, drought, extreme temperatures, pathogens, and oxidative damage[55,56].

    Meanwhile, Ka/Ks analysis between the two altitudinal ecotypes identified that DEGs involved in the phenylpropanoid and flavonoid biosynthesis pathways, such as FLS, CCoAOMT, and HCT, experienced significant positive selection (Table 2). Furthermore, population genetic studies also identified alleles of nine genes (FLS, HCT, PAL, C3H, FOMT, FNS, F3H, CYP75B4, and CHS) out of the 15 tested genes from the phenylpropanoid and flavonoid biosynthesis pathways were fixed in the high-altitude populations (Supplementary Table S6). Previous studies based on RAD-seq and several neutral genetic markers have elucidated that there is no significant genetic differentiation between the high- and middle-altitude populations. This indicates that the high-altitude populations and mid-altitude populations share the same origination and dispersal patterns[19,26]. Thus, the fixation of alleles for these genes in the high-altitude populations might have occurred after populations' colonization on the QXP. Previous simulation and population genetic analyses in wild barley have revealed that advantageous alleles could be selectively preserved and tend to become fixed within the populations under strong environmental pressures[57]. To survive the harsh desert conditions of the QXP, advantageous alleles in the genes related to the phenylpropanoid and flavonoid biosynthesis pathways might be selectively retained. Then, under the environment's directional selection, these alleles could progressively replace other alleles and eventually establish themselves within the high-altitude populations.

    Interestingly, Ka/Ks analysis revealed that CCoAOMT, another critical gene in the biosynthesis of flavonoids and phenylpropanoids, is under diversifying selection between two altitudinal ecotypes (Table 2). Additionally, signals of balancing selection on this gene were also observed among the high-altitude populations, and the mid-altitude populations as well (Supplementary Table S7). CCoAOMT is a key enzyme that catalyzes the methylation of caffeoyl-CoA to feruloyl-CoA. This reaction is crucial in the production of monolignols, the fundamental building blocks of lignin. Lignin is vital for maintaining the structural integrity of plants and enhancing their resistance to pathogens, and also play a significant role in enabling plants to withstand abiotic stresses such as cold and drought[58,59]. In recent years, increasing research underscores that balancing selection, by preserving genetic diversity, is a fundamental evolutionary mechanism significantly contributing to the adaptability and survival of species in changing environments, ensuring populations can cope with new challenges and thrive across diverse ecological niches[26,60]. This result suggested that the balancing selection of genes involved in the phenylpropanoid and flavonoid biosynthesis pathways plays a pivotal role in enabling sandrice to endure and prosper under the varied extreme conditions of desert environments, including the harsher environment of alpine deserts.

    In general, to survive and thrive in the diverse and extreme desert environment for sandrice, functional genes such as those involved in phenylpropanoid and flavonoid biosynthesis pathways experienced balanced selection among populations with different ecological niches, which could preserve high genetic diversity to ensure populations that can cope with diverse challenges. However, populations of sandrice inhabiting the alpine deserts endure even harsher conditions compared to those in other northern deserts. Under such environmental stress and selective pressure, a greater number of functional genes of phenylpropanoid and flavonoid biosynthesis pathways undergo directional selection, resulting in a shift in the population's genetic makeup toward certain favorable alleles. In some cases, this process leads to the presence of only one advantageous allele for specific functional genes within the population, influencing gene expression and subsequently regulating the high accumulation of downstream flavonoids and phenylpropanoids to enable the population to thrive in the harsh environments of the QXP deserts. As a result of long-term local adaptation, the accumulation of flavonoids and phenylpropanoids in sandrice has significantly diverged among ecotypes from different altitudes, even within the same common garden[27,28].

    In addition to the phenylpropanoid and flavonoid biosynthesis pathways, pathways of photosynthesis, starch and sucrose metabolism, flavonoid synthesis, phenylpropanoid synthesis, ribosomal regulation, and carbon metabolism, and several genes related to chronic hypoxia, oxidative stress, DNA damage repair, and stress response regulation, such as ANG, UNG, PRX3, ELP4, CDPK1, and GP1, are also suggested to play crucial roles in sandrice's adaptation to the harsh desert environment of QXP (Fig. 2, Table 2). To gain a comprehensive understanding of sandrice's adaptation mechanisms to alpine deserts and the formation of medicinal components in sandrice under high-altitude environmental factors, further genetic evidence is needed to validate the molecular functions and regulatory relationships among these adaptive alleles and pathways, the expression of functional genes, and the subsequent synthesis and accumulation of those anti-stress metabolites.

  • The authors confirm contribution to the paper as follows: study conception and design: Ma XF, Yan X, Yin X; data collection: Yin X, Qian C, Fan X, Zhou S; analysis and interpretation of results: Yin X, Qian C, Fang T, Yang J, Chang Y; draft manuscript preparation: Qian C, Yin X, Yan X, Chang Y. All authors reviewed the results and approved the final version of the manuscript.

  • All the raw sequenced data was submitted to GenBank (www.ncbi.nlm.nih.gov) (accession numbers: PRJNA659807; OM338032-OM338057, OP846852-OP846955).

  • This research was supported by the National Natural Science Foundation of China (Grant No. 32271695); Lanzhou Youth Science and Technology Talent Innovation Project (Grant No. 2023-QN-140); Chinese Academy of Sciences Strategic Biological Resources Program (Grant No. KFJ-BRP-007-015); Gansu Province to Guide the Development of Scientific and Technological Innovation Special Fund Competitive Project (Grant No. Y939BD1001), and National Natural Science Foundation of China (NSFC, Grant Nos 32171608 and 32201378).

  • The authors declare that they have no conflict of interest.

  • Supplemental Tables S1 Compound information in green tea (GT).
    Supplemental Tables S2 Compound information in black tea (BT).
    Supplemental Tables S3 Analysis of the main differential active ingredients associated with neurodegenerative changes PPI in green tea (GT) compared to black tea (BT).
    Supplemental Tables S4 Analysis of the main differential active ingredients associated with neurodegenerative changes PPI in black tea (BT) compared to green tea (GT).
  • [1]

    Opare SKA, Rauk A. 2019. Pseudopeptide designed to inhibit oligomerization and redox chemistry in Alzheimer's disease. The Journal of Physical Chemistry 123:5206−15

    doi: 10.1021/acs.jpcb.9b01665

    CrossRef   Google Scholar

    [2]

    Samanta S, Rajasekhar K, Babagond V, Govindaraju T. 2019. Small molecule inhibits metal-dependent and -independent multifaceted toxicity of Alzheimer's disease. ACS Chemical Neuroscience 10:3611−21

    doi: 10.1021/acschemneuro.9b00216

    CrossRef   Google Scholar

    [3]

    Chen Z, Zhong C. 2014. Oxidative stress in Alzheimer's disease. Neurosci Bull 30:271−81

    doi: 10.1007/s12264-013-1423-y

    CrossRef   Google Scholar

    [4]

    Forloni G, Balducci C. 2018. Alzheimer's disease, oligomers, and inflammation. Journal of Alzheimer’s Disease 62:1261−76

    doi: 10.3233/JAD-170819

    CrossRef   Google Scholar

    [5]

    Shoshan-Barmatz V, Nahon-Crystal E, Shteinfer-Kuzmine A, Gupta R. 2018. VDAC1, mitochondrial dysfunction, and Alzheimer's disease. Pharmacological Research 131:87−101

    doi: 10.1016/j.phrs.2018.03.010

    CrossRef   Google Scholar

    [6]

    Naseri NN, Wang H, Guo J, Sharma M, Luo W. 2019. The complexity of tau in Alzheimer's disease. Neuroscience Letters 705:183−94

    doi: 10.1016/j.neulet.2019.04.022

    CrossRef   Google Scholar

    [7]

    Xu T, Niu C, Zhang X, Dong M. 2018. β-Ecdysterone protects SH-SY5Y cells against β-amyloid-induced apoptosis via c-Jun N-terminal kinase- and Akt-associated complementary pathways. Laboratory Investigation 98:489−99

    doi: 10.1038/s41374-017-0009-0

    CrossRef   Google Scholar

    [8]

    Neumann U, Ufer M, Jacobson LH, Rouzade-Dominguez ML, Huledal G, et al. 2018. The BACE-1 inhibitor CNP520 for prevention trials in Alzheimer's disease. EMBO Molecular Medicine 10:e9316

    doi: 10.15252/emmm.201809316

    CrossRef   Google Scholar

    [9]

    Reiss AB, Arain HA, Stecker MM, Siegart NM, Kasselman LJ. 2018. Amyloid toxicity in Alzheimer's disease. Reviews in the Neurosciences 29:613−27

    doi: 10.1515/revneuro-2017-0063

    CrossRef   Google Scholar

    [10]

    Chao AC, Chen CH, Wu MH, Hou BY, Yang DI. 2020. Roles of Id1/HIF-1 and CDK5/HIF-1 in cell cycle reentry induced by amyloid-beta peptide in post-mitotic cortical neuron. Biochimica et Biophysica Acta - Molecular Cell Research 1867:118628

    doi: 10.1016/j.bbamcr.2019.118628

    CrossRef   Google Scholar

    [11]

    Caldeira GL, Ferreira IL, Rego AC. 2013. Impaired transcription in Alzheimer's disease: key role in mitochondrial dysfunction and oxidative stress. Journal of Alzheimer’s Disease 34:115−31

    doi: 10.3233/JAD-121444

    CrossRef   Google Scholar

    [12]

    Spangenberg EE, Green KN. 2017. Inflammation in Alzheimer's disease: Lessons learned from microglia-depletion models. Brain, Behavior, and Immunity 61:1−11

    doi: 10.1016/j.bbi.2016.07.003

    CrossRef   Google Scholar

    [13]

    Zhou L, Chen D, Huang XM, Long F, Cai H, et al. 2017. Wnt5a promotes cortical neuron survival by inhibiting cell-cycle activation. Frontiers in Cellular Neuroscience 11:281

    doi: 10.3389/fncel.2017.00281

    CrossRef   Google Scholar

    [14]

    Schimidt HL, Garcia A, Martins A, Mello-Carpes PB, Carpes FP. 2017. Green tea supplementation produces better neuroprotective effects than red and black tea in Alzheimer-like rat model. Food Research International 100:442−48

    doi: 10.1016/j.foodres.2017.07.026

    CrossRef   Google Scholar

    [15]

    Hidalgo FJ, Delgado RM, Zamora R. 2017. Protective effect of phenolic compounds on carbonyl-amine reactions produced by lipid-derived reactive carbonyls. Food Chemistry 229:388−95

    doi: 10.1016/j.foodchem.2017.02.097

    CrossRef   Google Scholar

    [16]

    Pan SY, Nie Q, Tai HC, Song XL, Tong YF, et al. 2022. Tea and tea drinking: China's outstanding contributions to the mankind. Chinese Medicine 17:27

    doi: 10.1186/s13020-022-00571-1

    CrossRef   Google Scholar

    [17]

    Fei T, Fei J, Huang F, Xie T, Xu J, et al. 2017. The anti-aging and anti-oxidation effects of tea water extract in Caenorhabditis elegans. Experimental Gerontology 97:89−96

    doi: 10.1016/j.exger.2017.07.015

    CrossRef   Google Scholar

    [18]

    Pan H, Gao Y, Tu Y. 2016. Mechanisms of Body Weight Reduction by Black Tea Polyphenols. Molecules 21:1659

    doi: 10.3390/molecules21121659

    CrossRef   Google Scholar

    [19]

    Zhao T, Li C, Wang S, Song X. 2022. Green tea (Camellia sinensis): a review of its phytochemistry, pharmacology, and toxicology. Molecules 27:3909

    doi: 10.3390/molecules27123909

    CrossRef   Google Scholar

    [20]

    Wang RH, Zhu XF, Qian W, Tang HY, Jiang J, et al. 2018. Effect of tea polyphenols on copper adsorption and manganese release in two variable-charge soils. Journal of Geochemical Exploration 190:374−80

    doi: 10.1016/j.gexplo.2018.04.006

    CrossRef   Google Scholar

    [21]

    Deb S, Dutta A, Phukan BC, Manivasagam T, Justin Thenmozhi A, et al. 2019. Neuroprotective attributes of L-theanine, a bioactive amino acid of tea, and its potential role in Parkinson's disease therapeutics. Neurochemistry International 129:104478

    doi: 10.1016/j.neuint.2019.104478

    CrossRef   Google Scholar

    [22]

    Zhao T, Tang H, Xie L, Zheng Y, Ma Z, et al. 2019. Scutellaria baicalensis Georgi. (Lamiaceae): a review of its traditional uses, botany, phytochemistry, pharmacology and toxicology. Journal of Pharmacy and Pharmacology 71:1353−69

    doi: 10.1111/jphp.13129

    CrossRef   Google Scholar

    [23]

    Tanaka T, Matsuo Y. 2020. Production Mechanisms of Black Tea Polyphenols. Chemical & Pharmaceutical Bulletin 68:1131−42

    doi: 10.1248/cpb.c20-00295

    CrossRef   Google Scholar

    [24]

    Li X, Smid SD, Lin J, Gong Z, Chen S, et al. 2019. Neuroprotective and anti-amyloid β effect and main chemical profiles of white tea: comparison against green, oolong and black tea. Molecules 24:1926

    doi: 10.3390/molecules24101926

    CrossRef   Google Scholar

    [25]

    Wan J, Feng M, Pan W, Zheng X, Xie X, et al. 2021. Inhibitory effects of six types of tea on aging and high-fat diet-related amyloid formation activities. Antioxidants 10:1513

    doi: 10.3390/antiox10101513

    CrossRef   Google Scholar

    [26]

    Zheng X, Feng M, Wan J, Shi Y, Xie X, et al. 2021. Anti-damage effect of theaflavin-3'-gallate from black tea on UVB-irradiated HaCaT cells by photoprotection and maintaining cell homeostasis. Journal of Photochemistry and Photobiology B, Biology 224:112304

    doi: 10.1016/j.jphotobiol.2021.112304

    CrossRef   Google Scholar

    [27]

    Ashraf GM, Greig NH, Khan TA, Hassan I, Tabrez S, et al. 2014. Protein misfolding and aggregation in Alzheimer's disease and type 2 diabetes mellitus. CNS & Neurological Disorders Drug Targets 13:1280−93

    doi: 10.2174/1871527313666140917095514

    CrossRef   Google Scholar

    [28]

    Cai S, Yang H, Zeng K, Zhang J, Zhong N, et al. 2016. EGCG Inhibited Lipofuscin Formation Based on Intercepting Amyloidogenic β-Sheet-Rich Structure Conversion. PLoS One 11:e0152064

    doi: 10.1371/journal.pone.0152064

    CrossRef   Google Scholar

    [29]

    Skene JHP. 1989. Axonal growth-associated proteins. Annual Review of Neuroscience 12:127−56

    doi: 10.1146/annurev.ne.12.030189.001015

    CrossRef   Google Scholar

    [30]

    Jöbstl E, Fairclough JPA, Davies AP, Williamson MP. 2005. Creaming in black tea. Journal of Agricultural and Food Chemistry 53:7997−8002

    doi: 10.1021/jf0506479

    CrossRef   Google Scholar

    [31]

    Lim IK. 2006. TIS21 /BTG2/PC3 as a link between ageing and cancer: cell cycle regulator and endogenous cell death molecule. Journal of Cancer Research and Clinical Oncology 132:417−26

    doi: 10.1007/s00432-006-0080-1

    CrossRef   Google Scholar

    [32]

    Loreto A, Hill CS, Hewitt VL, Orsomando G, Angeletti C, et al. 2020. Mitochondrial impairment activates the Wallerian pathway through depletion of NMNAT2 leading to SARM1-dependent axon degeneration. Neurobiology of Disease 134:104678

    doi: 10.1016/j.nbd.2019.104678

    CrossRef   Google Scholar

    [33]

    Cambron M, D'Haeseleer M, Laureys G, Clinckers R, Debruyne J, De Keyser J. 2012. White-matter astrocytes, axonal energy metabolism, and axonal degeneration in multiple sclerosis. Journal of Cerebral Blood Flow and Metabolism 32:413−24

    doi: 10.1038/jcbfm.2011.193

    CrossRef   Google Scholar

    [34]

    Chaturvedi RK, Flint Beal M. 2013. Mitochondrial diseases of the brain. Free Radical Biology & Medicine 63:1−29

    doi: 10.1016/j.freeradbiomed.2013.03.018

    CrossRef   Google Scholar

    [35]

    Kalia LV, Gingrich JR, Salter MW. 2004. Src in synaptic transmission and plasticity. Oncogene 23:8007−16

    doi: 10.1038/sj.onc.1208158

    CrossRef   Google Scholar

    [36]

    Malm T, Loppi S, Kanninen KM. 2016. Exosomes in Alzheimer's disease. Neurochemistry International 97:193−99

    doi: 10.1016/j.neuint.2016.04.011

    CrossRef   Google Scholar

    [37]

    Lezi E, Swerdlow RH. 2012. Mitochondria in neurodegeneration. Advances in Experimental Medicine and Biology 942:269−86

    doi: 10.1007/978-94-007-2869-1_12

    CrossRef   Google Scholar

    [38]

    Zhao X, Fang J, Li S, Gaur U, Xing X, et al. 2019. Artemisinin attenuated hydrogen peroxide H2O2-induced oxidative injury in SH-SY5Y and hippocampal neurons via the activation of AMPK pathway. International Journal of Molecular Sciences 20:2680

    doi: 10.3390/ijms20112680

    CrossRef   Google Scholar

    [39]

    Crouch PJ, Harding SM, White AR, Camakaris J, Bush AI, et al. 2008. Mechanisms of A beta mediated neurodegeneration in Alzheimer's disease. The International Journal of Biochemistry & Cell Biology 40:181−98

    doi: 10.1016/j.biocel.2007.07.013

    CrossRef   Google Scholar

    [40]

    Padmanabhan S, Burke RE. 2018. Induction of axon growth in the adult brain: A new approach to restoration in Parkinson's disease. Movement Disorders 33:62−70

    doi: 10.1002/mds.27209

    CrossRef   Google Scholar

    [41]

    Yin H, Xu L, Porter NA. 2011. Free radical lipid peroxidation: mechanisms and analysis. Chemical Reviews 111:5944−72

    doi: 10.1021/cr200084z

    CrossRef   Google Scholar

    [42]

    Orlowski RZ. 1999. The role of the ubiquitin-proteasome pathway in apoptosis. Cell Death and Differentiation 6:303−13

    doi: 10.1038/sj.cdd.4400505

    CrossRef   Google Scholar

    [43]

    Kocaturk NM, Gozuacik D. 2018. Crosstalk Between Mammalian Autophagy and the Ubiquitin-Proteasome System. Frontiers in Cell and Developmental Biology 6:128

    doi: 10.3389/fcell.2018.00128

    CrossRef   Google Scholar

    [44]

    Martin-Rincon M, Pérez-López A, Morales-Alamo D, Perez-Suarez I, de Pablos-Velasco P, et al. 2019. Exercise mitigates the loss of muscle mass by attenuating the activation of autophagy during severe energy deficit. Nutrients 11:2824

    doi: 10.3390/nu11112824

    CrossRef   Google Scholar

    [45]

    Zaffagnini G, Savova A, Danieli A, Romanov J, Tremel S, et al. 2018. p62 filaments capture and present ubiquitinated cargos for autophagy. The EMBO Journal 37:e98308

    doi: 10.15252/embj.201798308

    CrossRef   Google Scholar

    [46]

    Daemen S, van Zandvoort MAMJ, Parekh SH, Hesselink MKC. 2016. Microscopy tools for the investigation of intracellular lipid storage and dynamics. Molecular Metabolism 5:153−63

    doi: 10.1016/j.molmet.2015.12.005

    CrossRef   Google Scholar

    [47]

    Loix M, Wouters E, Vanherle S, Dehairs J, McManaman JL, et al. 2022. Perilipin-2 limits remyelination by preventing lipid droplet degradation. Cellular and Molecular Life Sciences 79:515

    doi: 10.1007/s00018-022-04547-0

    CrossRef   Google Scholar

    [48]

    Veliova M, Petcherski A, Liesa M, Shirihai OS. 2020. The biology of lipid droplet-bound mitochondria. Seminars in Cell & Developmental Biology 108:55−64

    doi: 10.1016/j.semcdb.2020.04.013

    CrossRef   Google Scholar

    [49]

    Moreno JA, Radford H, Peretti D, Steinert JR, Verity N, et al. 2012. Sustained translational repression by eIF2α-P mediates prion neurodegeneration. Nature 485:507−11

    doi: 10.1038/nature11058

    CrossRef   Google Scholar

    [50]

    Jishi A, Qi X, Miranda HC. 2021. Implications of mRNA translation dysregulation for neurological disorders. Seminars in Cell & Deve lopmental Biology 114:11−19

    doi: 10.1016/j.semcdb.2020.09.005

    CrossRef   Google Scholar

    [51]

    Meier S, Bell M, Lyons DN, Rodriguez-Rivera J, Ingram A, et al. 2016. Pathological tau promotes neuronal damage by impairing ribosomal function and decreasing protein synthesis. The Journal of Neuroscience 36:1001−7

    doi: 10.1523/JNEUROSCI.3029-15.2016

    CrossRef   Google Scholar

    [52]

    Stein KC, Morales-Polanco F, van der Lienden J, Rainbolt TK, Frydman J. 2022. Ageing exacerbates ribosome pausing to disrupt cotranslational proteostasis. Nature 601:637−42

    doi: 10.1038/s41586-021-04295-4

    CrossRef   Google Scholar

  • Cite this article

    Wu H, Wan J, Yuan J, Xie M, Nie Q, et al. 2024. Neuroprotective comparisons and bioactive profiles of green tea and black tea: in vitro cellular experiments, metabolomics, and network pharmacology analysis. Beverage Plant Research 4: e017 doi: 10.48130/bpr-0024-0019
    Wu H, Wan J, Yuan J, Xie M, Nie Q, et al. 2024. Neuroprotective comparisons and bioactive profiles of green tea and black tea: in vitro cellular experiments, metabolomics, and network pharmacology analysis. Beverage Plant Research 4: e017 doi: 10.48130/bpr-0024-0019

Figures(6)

Article Metrics

Article views(2986) PDF downloads(331)

ARTICLE   Open Access    

Neuroprotective comparisons and bioactive profiles of green tea and black tea: in vitro cellular experiments, metabolomics, and network pharmacology analysis

Beverage Plant Research  4 Article number: e017  (2024)  |  Cite this article

Abstract: This study aimed to compare the neuroprotective effects of green tea and black tea, made from the same raw materials, on an Aβ25-35-induced PC12 cell model, using Ultra-High Performance Liquid Chromatography-Tandem Mass Spectrometry (UPLC-MS/MS) and network pharmacology approaches. This addresses the gap in current research, which has extensively explored the neuroprotective properties of green tea and its components but has paid less attention to black tea. The findings indicate that both teas can alleviate Aβ-induced neurodegenerative changes by reducing inflammation, mitochondrial disruption, and other cellular stressors. Notably, black tea showed higher effectiveness, enriching more differentially expressed genes within critical pathways and exhibiting a broader spectrum of bioactive compounds. Its protein-protein interaction network also suggested that black tea acted on a wider range of potential targets. However, these results are preliminary and emphasize the importance of the complex interplay of bioactive components in tea, advocating for further comparative studies to fully understand their neuroprotective mechanisms.

    • Increasing evidence suggests the imbalance between the accumulation and clearance of amyloid β (Aβ) and Aβ-related peptides plays a crucial role in the onset of Alzheimer’s disease (AD)[1,2]. Aβ is known to trigger various stress responses, including oxidative stress, inflammation, and mitochondrial dysfunction[35]. Externally, Aβ forms senile plaques outside of cells, while intracellular hyperphosphorylated tau protein forms neurofibrillary tangles and cell apoptosis, both of which are important factors leading to the development and progression of AD[68].

      Moreover, Aβ exerts multiple toxic effects, such as causing ion leakage, disruption of cellular calcium ion balance, and induction of cell apoptosis[9]. Beyond external accumulation, Aβ can also accumulate in cell membranes, organelles, and the cell nucleus, leading to mitochondrial dysfunction, endoplasmic reticulum stress, and abnormal cell cycle[10,11]. This buildup triggers microglia activation, leading to the production and release of pro-inflammatory cytokines, including IL-6 and TNF-α[12]. Additionally, Aβ can stimulate the reactivation of the mitotic cycle in neurons, hastening their apoptosis[13].

      Numerous investigations have indicated that tea and its functional components can delay brain aging[10,1418]. Key active ingredients in green tea, such as catechins and theanine, are noted for their ability to prevent the accumulation of Aβ aggregates, while also providing anti-inflammatory, antioxidant, and axonal growth stimulation[16,1922]. Black tea accounts for approximately 80% of the world's tea production and has the quality characteristics of 'red tea, red soup, red leaves, and mellow taste'[23]. Research has found that black tea can inhibit the formation of Aβ aggregates and has significant protective effects against Aβ-induced neurotoxicity[24]. In a rat model of attention deficit hyperactivity disorder (ADHD), black tea has been shown to effectively counteract Aβ-induced neurotoxicity and prevent memory loss[10,25,26]. The exact mechanisms behind black tea's protective effects in the brain remain to be fully elucidated.

      This research established a neurodegeneration model of Alzheimer's disease (AD) using PC12 cells exposed to Aβ25-35 to compare the effectiveness of green tea and black teas, derived from the same batch of leaves, in mitigating Aβ-induced cellular dysfunction. Furthermore, this study employed a comprehensive metabolomic analysis of tea components, alongside network pharmacology approaches, to explore the neuroprotective qualities of both green tea and black tea.

    • The tea raw materials used in the experiment were commercially available green tea and black tea. Fresh tea leaves of the same variety and from the same period of time were used to make the tea. 3-(4,5-Dimethyl-2-Thiazolyl)-2,5-Diphenyl Tetrazolium Bromide (MTT), and Dimethyl sulfoxide (DMSO) were purchased from Sigma-Aldrich (St Louis, MO, USA). RIPA lysate, protease inhibitor mixture, and anti-fluorescence quenching agent (containing DAPI) were purchased from Beyotime Biotechnology (Shanghai, China). In addition to anti-GAPDH, anti-IL-6, anti-NF-κB, anti-Histone H3, anti-Sirt1, anti-Cyclin B1, anti-Cyclin D1, anti-Bax, anti-Bcl-2 and anti-AMPK (Cell Signaling, Boston, MA, USA), the following primary antibodies were also used for Western blot analysis: anti-multiubiquitin (UPs) (Medical & Biological Laboratories co. Ltd, Tokyo, Japan), anti-4-HNE (Millipore, Boston, USA), anti-sequestosome-1 (p62) (Epitomics, Burlingame, CA, USA), anti-Klotho (Novusbio, Littleton, USA) and anti-RAGE (Santa Cruz Biotechnology, Dallas, TX, USA). Rabbit anti-Goat IgG-HRP Antibody, Goat anti-Rabbit IgG-HRP Antibody, Goat anti-Mouse IgG-HRP Antibody were purchased from Absin Biotechnology (Shanghai, China). Western chemiluminescent horseradish peroxidase substrate was purchased from Millipore (Boston, USA).

    • To accurately simulate the neuroprotective effects of tea consumption, this study utilized a simulated tea brewing method and a direct freeze-drying process to accurately preserve the components and particle structures of the tea infusion. Specifically, 100 g each of green tea and black tea were accurately weighed and placed into conical flasks containing 1 L of boiling water (100 °C), followed by a 45-min immersion in a 100 °C water bath, filtered, and then subjected to another liter of boiling water for an additional 45 min in the water bath. The infusions from each tea were combined, cooled to room temperature, pre-frozen at −20 °C for 24 h, and freeze-dried to produce green tea and black tea water extracts, subsequently stored at −20 °C.

    • 25-35 (1 mM), Aβ25-35 (1 mM)/green tea (1 mM) and Aβ25-35 (1 mM)/black tea (1 mM) were preincubation in sterile water for 7 d in a constant temperature incubator at 37 °C to obtain Aβ25-35, Aβ25-35/green tea, Aβ25-35/black tea proteins.

    • Thioflavin T (ThT) was used to analyze the β-sheet structure generated during amyloid aggregation[25]. Aβ25-35 (50 μM), Aβ25-35 (50 μM)/green tea (50 μg/mL) and Aβ25-35 (50 μM)/black Tea (50 μg/mL) were incubated at 37 °C for 7 d. In subsequent experiments, these groups of proteins were abbreviated as Aβ25-35, Aβ25-35/green tea and Aβ25-35/black tea, respectively. ThT was assayed as follows: 40 μL of different protein solutions and 160 μL of ThT working solution were aspirated in a 96-well plate. The fluorescence intensity was measured on a multifunctional microplate apparatus at an excitation wavelength of 440 nm and an emission wavelength of 485 nm.

    • PC12 cells were inoculated at 1 × 105/mL in a culture dish containing DMEM medium with 10% FBS and incubated in a constant temperature incubator at 37 °C containing 5% CO2. The medium was changed every 2 d, and when the cells reached 60%−70% fusion rate, the cells were treated with Aβ25-35, Aβ25-35/green tea, and Aβ25-35/black tea for 24 h. For the control group (Control) cells, equal amounts of sterile water were added.

    • PC12 cells were seeded at 1 × 104 cells/well in a 96-well plate and cultured for 24 h. Cells were incubated with the different protein samples (Aβ25-35, Aβ25-35/green tea, Aβ25-35/black tea) prepared as above and the control group was added with the same amount of sterile water. After 24 h, the cell viability of different treatment groups was detected using the MTT method[26].

    • Cells were inoculated in 6-well culture plates at 1.5 × 105 cells/well and cultured as described above. Then the supernatant was aspirated, the cells were washed twice by adding pre-cooled PBS, and 200 μL of lysate was added, and the supernatant was collected by centrifugation (4 °C, 12,000 rpm, 5 min) to determine the ATP content (Beyotime Biotechnology, Shanghai, China).

    • Cells were inoculated in 6-well culture plates at a density of 1.5 × 105 cells/well, and cell culture was performed as described above. The cells were collected and suspended in diluted DCFH-DA (10 μM) and incubated in a cell culture incubator at 37 °C for 20 min. The cells were washed three times with a serum-free cell culture medium to remove the DCFH-DA that had not entered the cells fully. Then, the fluorescence values were detected at 488 nm excitation wavelength and 525 nm emission wavelength using a fluorescence zymograph.

    • Cells were inoculated at a density of 5 × 104 cells/well onto cell crawls in 24-well culture plates, and cell culture was performed as described above. The supernatant was then aspirated and the cells were washed twice by adding pre-cooled PBS. Cell culture medium (250 μL) was added to each well, then 250 μL of JC-1 staining working solution was added to each well and incubated in a cell incubator at 37 °C for 20 min. 300 μL of JC-1 staining buffer (1×) was added to wash the cells twice, followed by blocking the slides with an anti-fluorescent bursting agent (containing DAPI). Images were acquired using a fluorescence microscope (Zeiss, Axio scope. A1). The average fluorescence intensity of the images was measured using ZEN software package. Three measurements were taken, each time measuring an area of the same size.

    • Cellular lipid droplet staining was performed using the BODIPY fluorescent probe (Invitrogen, New York, USA). Cells were treated according to the methods described above. After the experiment, the cells were collected and stained strictly following the kit's instructions. Finally, photographs were taken with a fluorescence microscope (Zeiss, Axio scope. A1), and the average fluorescence intensities of different treatment groups were calculated using ZEN software.

    • Cells were inoculated at 6 × 105 cells/well in a 10 cm diameter cell culture dish, and cell culture was performed as described above. The supernatant was then aspirated, cells were washed twice by adding pre-cooled PBS, and 500 μL of RIPA lysate containing protease inhibitors and the supernatant was collected by centrifugation (4 °C, 12,000 rpm, 20 min). The protein concentration was determined by the BCA protein analysis kit (Pierce, Grand Island, NY, USA). Equal amounts of protein from each sample were separated on 8%−10% SDS-PAGE gels. Proteins from the gels were transferred to PVDF (Millipore, USA) membranes. The membranes were closed in TBST buffer containing 5% skim milk powder for 1 h. The membranes were then incubated overnight at 4 °C in the primary antibody at dilution ratios of 1:500−1:2000. After washing in TBST buffer, the membranes were incubated with appropriate secondary antibodies for 90 min at room temperature. Signals were detected using enhanced chemiluminescence (ECL) detection kit (Millipore, USA) and analyzed by Western blot densitometry using Image-J software[26].

    • ELISA detection kits (Jiangsu Feiya Biotechnology, Yancheng, China) were used to determine the contents of GAP43 and TNF-α in the protein samples of different treatment groups prepared above.

    • Cell experiments were performed according to the methods above. After the experiment, the collected cell samples were snap-frozen in liquid nitrogen, frozen in dry ice, and sent to BGI Company (Shenzhen, China) for subsequent RNA-seq detection (BGISEQ-500 sequencer). In order to explore molecular pathways and networks more comprehensively from the transcriptome level, a sufficient number of differential genes (DEGs) were screened with |fold change|≥1.2 and Q ≤ 0.05 for heatmap, KEGG and gene interaction analysis. All analyses were performed using the online bioinformatics platform Dr. Tom (biosys.bgi.com/) provided by BGI.

    • Pre-treatment, extract analysis, metabolite identification and quantification of green tea and black tea samples were performed at BGI-Shenzhen (Shenzhen, China) according to their standard procedures. Untargeted Metabolomics analysis was performed by UPLC-MS/MS technique using a high-resolution mass spectrometer Q Exactive HF (Thermo Fisher Scientific, USA) with separate data acquisition in both positive and negative ion modes to improve metabolite coverage. The VIP values of the first two principal components were modelled using Partial Least Squares Method-Discriminant Analysis (PLS-DA), and the results of Fold change and Student's t test obtained from univariate analysis were combined to screen for differential metabolites.

    • The screened compounds were converted into the Standard Simplified Molecular Input Line Entry System (SMILES) by Pubchem (https://pubchem.ncbi.nlm.nih.gov/) and screened for drug similarity by blood-brain barrier (BBB) permeation and SwissADME calculations (www.swissadme.ch), active compounds with a BBB of 'yes' and Bioavailability Score > 0.17 are considered to have good bioavailability. Probability > 0 was used as the screening threshold, pharmacological targets were obtained using the SwissTargetPrediction database. In the GeneCard database (www.genecards.org), the disease genes of 'Neurodegenerative changes' were searched for intersection targets for subsequent analysis. The Venn diagram was used to evaluate the potential targets of the intersection genes of green tea and black tea with neurodegenerative changes. GO and KEGG analysis of cross-targeted genes using the 'custom analysis' module in Metascape (https://metascape.org/gp/index.html), with species set to Homo sapiens and omicstudio (www.omicstudio.cn) visualized KEGG analysis of target genes. PPI analysis of targets was performed by STRING database (https://string-db.org/), with the target species set to Homo sapiens and a confidence level greater than 0.09. Then, the PPI network was visualized using Cytoscape software (version 3.9.1), with images showing only the core targets with degree values greater than the mean or above.

    • Statistical analysis was performed using the GraphPad Prism 8.01 software package. Combined with Turkey's multiple comparison test, one-way ANOVA was used to test the significance of differences, and the results were expressed as the mean ± standard deviation. p < 0.05 was significant, and p < 0.01 was judged to be extremely significant.

    • UPLC-MS/MS technology was utilized for untargeted metabolomics analysis, with data acquisition in both positive and negative ionization modes to enhance the range of detected metabolites. The base peak ion chromatograms revealed that the total ion current (TIC) profiles for metabolite identification overlapped significantly, demonstrating consistent retention times and peak intensities, which suggests reliable instrument signal stability (Fig. 1a & b). Partial least squares-discriminant analysis (PLS-DA) showed that samples of green and black tea were distinctly separated within the 95% confidence interval ellipse, effectively differentiating between the two tea varieties based on their component qualities (Fig. 1c). Utilizing criteria such as a VIP score ≥ 1, a Fold-Change range of 0.83 ≤ Fold-Change ≤ 1.2, and a q-value < 0.05, a total of 753 differential metabolites were screened. Among these, black tea demonstrated an increase in 570 metabolites and a decrease in 183 metabolites when compared to green tea (Fig. 1e). A heatmap further illustrated the propertional differences in various components between green tea and black tea (Fig. 1d).

      Figure 1. 

      Untargeted metabolomics analysis of aqueous extracts of green tea and black tea. (a), (b) Base peak ion chromatogram (BPC). (c) PLSDA Score Chart. (d) Heat map of differential metabolite clustering. Classification of compounds according to 'Family'. (e) Differential metabolite volcano map, VIP ≥ 1, 0.83 ≤ Fold-change of Difference ≤ 1.2, q-value < 0.05.

    • Research has shown that the formation of β-sheet structures plays a pivotal role in the early stages of amyloidgenesis[27,28]. Fluorescence detection revealed that the fluorescence intensity of Aβ25-35 was about sevenfold higher than that of the control (p < 0.01), suggesting a significant increase in β-sheet structures within the Aβ25-35 group. Treatment with both green tea and black tea markedly reduced the formation of these β-sheet structures (p < 0.01) (Fig. 2a). For simplicity, in the subsequent sections of the article, the treatments with Aβ25-35 (50 μM), Aβ25-35 (50 μM) combined with green tea (50 μg/ml), and Aβ25-35 (50 μM) combined with black tea (50 μg/ml) are referred to as Aβ25-35, Aβ25-35/green tea and Aβ25-35/black tea, respectively.

      Figure 2. 

      Green tea and black tea inhibits Aβ25-35-induced degenerative changes in differentiated PC12 cells. (a) ThT assay for β-fold structure content of different Aβ25-35 incubated samples. (b) MTT assay for cell viability. (c) DAPI fluorescence staining (bar = 20 μm). (d) ELISA assay for GAP43 protein expression level. (e) Westen-blotting for Aβ and other protein aggregate-related pathways. (f) Westen-blotting detection of cell cycle and apoptosis-related pathways. Compared with the Control group, # p < 0.05, ## p < 0.01. Compared with Aβ25-35 group, ** p < 0.01, n = 3.

      When well-differentiated PC12 cells were exposed to the various Aβ25-35 samples for 24 h, the viability of cells treated with Aβ25-35 alone decreased significantly (p < 0.01), as evidenced by enhanced DAPI staining and observations of fragmented or enlarged nuclei, indicating DNA damage. Conversely, green tea and black tea treatments significantly mitigated the toxicity induced by Aβ25-35. Notably, cell viability in the Aβ25-35/black tea group improved (p < 0.01), with DAPI staining results comparable to those of the control group (Fig. 2b & c).

      The expression of GAP43 in neurons plays a pivotal role in axon elongation, synapse formation, and neural germination during development[29]. Bax and Bcl-2 are homologous water-soluble related proteins, and overexpression of Bax antagonizes the protective effect of Bcl-2, leading to cell death[30]. Cyclin B1 drives the G2/M phase transition, and Cyclin D1 regulates the G1/S phase transition, both key cell cycle regulators[31]. ELISA and Western blot analyses (Fig. 2df) revealed that in the Aβ25-35-treated cells, the levels of GAP43 and Bcl-2 were significantly reduced (p < 0.01), whereas the levels of Bax, p62, UPs, β-Amyloid, RAGE, Cyclin B1 and Cyclin D1 were significantly elevated (p < 0.05 or p < 0.01). Conversely, these proteins exhibited reverse expression patterns in Aβ25-35/black tea group cells. The findings suggest that Aβ25-35 contributes to axonal atrophy, apoptosis, and the accumulation of toxic aggregates, thereby leading to aberrant cell cycle activation. Black tea exhibited a notable protective effect against the toxic stress induced by Aβ25-35, performing better than green tea in this regard.

    • Axons are rich in mitochondria, and mitochondrial dysfunction leads to neuronal axonal degeneration, the root cause of neurodegenerative diseases[3234]. In cells treated with Aβ25-35, there was a significant reduction in mitochondrial membrane potential (MMP) and ATP levels (p < 0.01), accompanied by an increase in ROS (p < 0.01) (Fig. 3ad). Further, ELISA and Western blot assays showed that Aβ25-35 promoted the nuclear translocation of NF-κB (p < 0.01), elevated the levels of TNF-α (p < 0.01), and decreased AMPK and Sirt1 protein levels (p < 0.01) (Fig. 3e & f). BODIPY staining for lipid droplets demonstrated that Aβ25-35 treatment led to increased cellular lipid deposition (p < 0.01) (Fig. 3g & h). These findings imply that Aβ25-35 triggers inflammatory pathways and suppresses cellular metabolism. Black tea exhibited a notable protective effect against mitochondrial damage, reduced inflammation and lipid accumulation, and enhanced cellular metabolic functions (p < 0.01), outperforming green tea in these respects.

      Figure 3. 

      Anti-inflammatory and metabolism-promoting functions of green tea and black tea. (a), (b) JC-1 staining plot with statistics (bar = 50 μm). (c) Fluorescence enzyme marker for ROS level. (d) ATP content assay. (e) ELISA for cellular TNF-α protein expression level. (f) Western-blotting for inflammation and energy metabolism-related pathways. (g), (h) BODIPY fluorescence staining with statistics for lipid droplets (bar = 10 μm). Compared with the Control group, ## p < 0.01; Compared with the Aβ25-35 group, * p < 0.05, ** p < 0.01, n = 3.

    • In comparison to the Aβ25-35 group, the DEGs in both the Aβ25-35/green tea and the Aβ25-35/black tea groups were significantly increased, with the DEGs in the Aβ25-35/black group being approximately 2.9 times higher than those in the Aβ25-35 group (Fig. 4a). KEGG pathway analysis revealed that these DEGs were predominantly involved in processes such as ribosomes function, related neurodegeneration, oxidative phosphorylation, AGE-RAGE signaling pathway, RNA transport, and DNA replication, among others (Fig. 4bd).

      Figure 4. 

      Transcriptome analysis of different treatment groups. (a) Comparative analysis of DEGs between groups. (b)−(d) KEGG analysis of DEGs in different treatment groups. (e) Heatmap and interaction network analysis of DEGs in different treatment groups with DEGs of Aβ25-35/black tea as a reference.

      Heatmap and protein-protein interaction network (PPI) analysis of DEGs in different treatment groups with DEGs of Aβ25-35/black tea as a reference. The expression patterns of DEGs in the tea-treated groups counteracted those observed in the Aβ25-35 group. Specifically, the DEGs in the Aβ25-35/black tea group predominantly enhanced pathways related to signal transduction, cell growth, and regeneration, while suppressing pathways associated with neurodegeneration, translation, transcription, protein folding, sorting and degradation, as well as networks related to endocrine and metabolic diseases. Moreover, black tea’s impact on the transcriptome was more pronounced than that of green tea (Fig. 4e).

    • To uncover the neuroprotective mechanisms of green tea and black tea on highly differentiated neuronal cells, network pharmacology approaches were employed to identify potential targets and pathways. Utilizing ADME (Absorption, Distribution, Metabolism, and Excretion) criteria, specifically a Bioavailability Score greater than 0.17, 46 active components in green tea and 155 in black tea that were upregulated were identified. Among these, 66 compounds from black tea and 14 from green tea were predicted to cross the blood-brain barrier (detailed compound information is provided in Supplemental Tables S1 & S2). The investigation highlighted 389 potential targets for green tea and 455 for black tea in combating neurodegenerative changes. GO (Gene Ontology) enrichment analysis revealed that both green tea and black tea target processes related to phosphorylation and protein kinase activity. Notably, green tea was associated with phosphotransferase activity, and alcohol group as acceptor, whereas black tea was linked to protein phosphorylation (Fig. 5a & b). KEGG analysis showed that the potential targets of both teas significantly impacted pathways involving neuroactive ligand-receptor interaction, cancer, PI3K-Akt signaling pathway, lipid and atherosclerosis, and MAPK signaling pathway. Despite the similarity in targeted pathways, black tea was found to enrich a greater number of genes within these pathways compared to green tea (Fig. 5c& d).

      Figure 5. 

      Network pharmacological analysis of green tea and black tea against neurodegenerative changes. (a), (b) VENN and GO analysis. (c), (d) KEGG analysis.

      The analysis of active component-target interaction maps revealed that green tea combats neurodegenerative changes through key components such as benzene and its derivatives, endogenous metabolites, nucleic acids, amino acids, alkaloids, coumarins, terpenoids. Black tea, while encompassing these types of components, also includes organic acids, lipids, phenolic acids, and flavonoids, offering a wider and more varied spectrum of neuroprotective targets. Black tea stands out with 113 main differential active compounds, compared to green tea’s 77, showcasing a greater diversity through higher compound interconnectivity. This suggests that black tea exhibits more complex biological activities (Fig. 6) (details of the compounds are provided in Supplemental Tables S3 & S4). SRC, identified as a pivotal target of black tea, is integral to central nervous system development and prevalent in neurons, highlighting black tea's nuanced neuroprotective capabilities[35].

      Figure 6. 

      Network diagram of the active component-key target of green tea and black tea in the intervention of neurodegenerative changes. (a) Association analysis plot for green tea's differential active components with protein-protein interactions (PPI). (b) Similar plot for black tea. Circles in the diagram now represent target genes, with their size and color intensity indicating the extent of their network connections and interaction strength, respectively. Quadrilaterals represent active components, with their size reflecting the degree of linkage to other genes.

    • This research indicates that both green tea and black tea share similar protective mechanisms against Aβ25-35-induced neurodegenerative changes in PC12 cells. Yet, black tea outperforms green tea in safeguarding mitochondria, preventing DNA damage, reducing lipid accumulation, and hindering the formation of aggregates (Figs 2 & 3). Further, transcriptomic analysis and network pharmacology studies suggest that black tea impacts a wider and more complex network of transcriptome and cellular signaling pathways, showcasing a richer array of bioactive compounds indicating a greater variety of bioactive components (Figs 4& 6).

      Aβ’s multiple toxic effects on neurons include inducing apoptosis, which contributes to neuronal loss[36]; and mitochondrial accumulation, which disrupts normal function[37]. These findings align with our study’s outcomes. Black tea pre-incubation can attenuate Aβ25-35-induced apoptosis of PC12 cells, enhance cell viability, and downregulate the expression of apoptotic pathway-related proteins (Fig. 2f). Mitochondrial dysfunction plays a key role in aging-related neurodegenerative diseases[38]. Aβ can impair neuronal function by damaging mitochondrial electron transport chains and inhibiting ATP production[39]. Our findings show that the pre-incubation with black tea significantly restores mitochondrial membrane potential, increases ATP levels, and elevates the expression of proteins in energy-metabolism pathways (Fig. 3a, df). For highly differentiated cells to undergo normal differentiation and maintain physiological functions, they must remain in a stationary phase of their cycle. Promoting axonal differentiation has been identified as crucial in combating neurodegenerative diseases[40]. Black tea appears to prevent neuronal degeneration by enforcing cycle arrest, outperforming green tea in this regard and aligning with previous research findings.

      Oxidative stress induced by Aβ leads to the oxidation of polyunsaturated fatty acids (PUFAs) in cell membranes through free radical chain reactions, resulting in the formation of lipid hydroperoxides such as 4-HNE[41]. These lipid hydroperoxides, including 4-HNE, can bind to nucleophilic functional groups in proteins, nucleic acids, and membrane lipids, thereby contributing to the impairment of autophagy[41]. When autophagy is disrupted, protein aggregates and damaged organelles accumulate within cells[42,43], leading to the cellular buildup of p62 and ubiquitin-modified proteins[44,45]. Black tea treatment effectively reduces the levels of ubiquitin- and p62-modified proteins, and suppresses inflammatory pathways (Fig. 2e & 3f).

      Excessive lipid storage is closely associated with metabolic abnormalities[46]. The accumulation of lipids in microglial cells can induced a pro-inflammatory state, hindering the repair mechanisms of the central nervous system. The proper management of lipid storage within cells is crucial for the maintenance of cellular energy balance[47,48]. Black tea has been shown to more effectively reduce the aberrant lipid accumulation induced by Aβ in PC12 cells compared to green tea (Fig. 3h).

      Given the prolonged half-lives of proteins in the brain, changes in mRNA translation can have enduring impacts on neuronal cells and potentially contribute to disease processes[49,50]. The biogenesis of ribosomes, which is critical for both the growth and upkeep of neuronal cells, is compromised under amyloid stress, leading to neuronal atrophy and loss of synaptic connections[4952]. Hence, targeting the dysregulation in protein synthesis might offer a viable strategy for restoring neuronal functionality. Transcriptome data revealed that black tea counters the disruptions in ribosome-related pathways induced by Aβ25-35. Analysis of DEGs suggested that black tea predominantly enhances axonal growth and signaling pathways, while it diminishes pathways involved in neurodegeneration, translation, and protein folding (Fig. 4e).

      Black tea exhibited significant metabolic regulation through the synergistic action of multiple components. The neurodegenerative mechanisms of green tea and black tea were explored using network pharmacology, highlighting how both teas influenced similar signaling pathways. However, the pathways influenced by black tea were characterized by a notably higher count of DEGs (Fig. 5c & d). Analysis of the effective component-target interactions showed that black tea, compared to green tea, contains distinctive components such as organic acids, lipids, phenolic acids, and flavonoids. The metabolic pathways affected by these components are more complex and interconnected (Fig. 6). Notably, the green tea and black tea used in our experiments were derived from the same raw materials harvested in early spring. Differences in the content of common active components like amino acids and caffeine between the two types of tea were minimal, leading to their exclusion from the differential metabolite analysis.

      In our investigation, we employed in vitro cell models to assess the neuroprotective properties of green tea and black tea, with a focus on Homo sapien target genes in network pharmacology analysis. Our findings indicated that black tea impacts a wider array of genes, pointing to a more potent anti-neurodegenerative activity. Nonetheless, these results should be approached with caution. Given the complexity of neurodegenerative conditions and the myriad of bioactive compounds present in tea, further confirmation via animal studies is imperative. Future research will extend to in vivo experiments to validate the neuroprotective efficacy of black tea and to elucidate its underlying mechanisms more thoroughly. This strategy aims to provide clearer insights into our findings and bolster the evidence supporting black tea's contribution to neuroprotection.

    • The authors confirm contribution to the paper as follows: methodology, data curation: Wu H, Wan J; investigation, visualization: Yuan J, Xie M, Nie Q; draft manuscript preparation: Wu H, Wan J, Cai S; manuscript review: Cai S; resources: Liu Z; funding acquisition, supervision: Liu Z, Cai S. All authors reviewed the results and approved the final version of the manuscript.

    • Due to administrative requirements, the datasets generated and/or analyzed during this study are not publicly available, but can be obtained from the corresponding author upon reasonable request.

      • This research was funded by the Guangxi Innovation Driven Development Special Fund Project (No. AA20302018), the National Key R&D Program of China (2018YFC1604405), the Key R&D Program of Hunan Province (2020WK2017), National Tea Industry Technology System Research Project of China (CARS-19-C01), National Natural Science Foundation Project of China (31471590, 31100501), and Self-Science Foundation of Hunan Province, China (2019jj50237).

      • The authors declare that they have no conflict of interest. Zhonghua Liu is the Editorial Board member of Beverage Plant Research who was blinded from reviewing or making decisions on the manuscript. The article was subject to the journal's standard procedures, with peer-review handled independently of this Editorial Board member and the research groups.

      • # Authors contributed equally: Huan Wu, Juan Wan

      • Supplemental Tables S1 Compound information in green tea (GT).
      • Supplemental Tables S2 Compound information in black tea (BT).
      • Supplemental Tables S3 Analysis of the main differential active ingredients associated with neurodegenerative changes PPI in green tea (GT) compared to black tea (BT).
      • Supplemental Tables S4 Analysis of the main differential active ingredients associated with neurodegenerative changes PPI in black tea (BT) compared to green tea (GT).
      • Copyright: © 2024 by the author(s). Published by Maximum Academic Press, Fayetteville, GA. This article is an open access article distributed under Creative Commons Attribution License (CC BY 4.0), visit https://creativecommons.org/licenses/by/4.0/.
    Figure (6)  References (52)
  • About this article
    Cite this article
    Wu H, Wan J, Yuan J, Xie M, Nie Q, et al. 2024. Neuroprotective comparisons and bioactive profiles of green tea and black tea: in vitro cellular experiments, metabolomics, and network pharmacology analysis. Beverage Plant Research 4: e017 doi: 10.48130/bpr-0024-0019
    Wu H, Wan J, Yuan J, Xie M, Nie Q, et al. 2024. Neuroprotective comparisons and bioactive profiles of green tea and black tea: in vitro cellular experiments, metabolomics, and network pharmacology analysis. Beverage Plant Research 4: e017 doi: 10.48130/bpr-0024-0019

Catalog

  • About this article

/

DownLoad:  Full-Size Img  PowerPoint
Return
Return