Transcriptomic and Translatomic Analyses Reveal Insights into the Developmental Regulation of Secondary Metabolism in the Young Shoots of Tea Plants (Camellia sinensis L.)

Ribosome Profiling and RNA Sequencing Reveal Translation and Transcription Regulation under Acute Heat Stress in Rainbow Trout (Oncorhynchus mykiss, Walbaum, 1792) Liver

Rainbow trout (Oncorhynchus mykiss, Walbaum, 1792) is an important economic cold-water fish that is susceptible to heat stress. To date, the heat stress response in rainbow trout is more widely understood at the transcriptional level, while little research has been conducted at the translational level. To reveal the translational regulation of heat stress in rainbow trout, in this study, we performed a ribosome profiling assay of rainbow trout liver under normal and heat stress conditions. Comparative analysis of the RNA-seq data with the ribosome profiling data showed that the folding changes in gene expression at the transcriptional level are moderately correlated with those at the translational level. In total, 1213 genes were significantly altered at the translational level. However, only 32.8% of the genes were common between both levels, demonstrating that heat stress is coordinated across both transcriptional and translational levels. Moreover, 809 genes exhibited significant differences in translational efficiency (TE), with the TE of these genes being considerably affected by factors such as the GC content, coding sequence length, and upstream open reading frame (uORF) presence. In addition, 3468 potential uORFs in 2676 genes were identified, which can potentially affect the TE of the main open reading frames. In this study, Ribo-seq and RNA-seq were used for the first time to elucidate the coordinated regulation of transcription and translation in rainbow trout under heat stress. These findings are expected to contribute novel data and theoretical insights to the international literature on the thermal stress response in fish.

Dissecting infectious bronchitis virus-induced host shutoff at the translation level

Viruses have evolved a range of strategies to utilize or manipulate the host's cellular translational machinery for efficient infection, although the mechanisms by which infectious bronchitis virus (IBV) manipulates the host translation machinery remain unclear. In this study, we firstly demonstrate that IBV infection causes host shutoff, although viral protein synthesis is not affected. We then screened 23 viral proteins, and identified that more than one viral protein is responsible for IBV-induced host shutoff, the inhibitory effects of proteins Nsp15 were particularly pronounced. Ribosome profiling was used to draw the landscape of viral mRNA and cellular genes expression model, and the results showed that IBV mRNAs gradually dominated the cellular mRNA pool, the translation efficiency of the viral mRNAs was lower than the median efficiency (about 1) of cellular mRNAs. In the analysis of viral transcription and translation, higher densities of RNA sequencing (RNA-seq) and ribosome profiling (Ribo-seq) reads were observed for structural proteins and 5' untranslated regions, which conformed to the typical transcriptional characteristics of nested viruses. Translational halt events and the number of host genes increased significantly after viral infection. The translationally paused genes were enriched in translation, unfolded-protein-related response, and activation of immune response pathways. Immune- and inflammation-related mRNAs were inefficiently translated in infected cells, and IBV infection delayed the production of IFN-β and IFN-λ. Our results describe the translational landscape of IBV-infected cells and demonstrate new strategies by which IBV induces host gene shutoff to promote its replication.

IMPORTANCE

Infectious bronchitis virus (IBV) is a γ-coronavirus that causes huge economic losses to the poultry industry. Understanding how the virus manipulates cellular biological processes to facilitate its replication is critical for controlling viral infections. Here, we used Ribo-seq to determine how IBV infection remodels the host's biological processes and identified multiple viral proteins involved in host gene shutoff. Immune- and inflammation-related mRNAs were inefficiently translated, the translation halt of unfolded proteins and immune activation-related genes increased significantly, benefitting IBV replication. These data provide new insights into how IBV modulates its host's antiviral responses.

Translatomics reveals the role of dietary calcium addition in regulating muscle fat deposition in pigs

Intramuscular fat (IMF) in pork holds significant importance for economic performance within the pig industry and dietary calcium supplementation enhances the accumulation of intramuscular fat. Additionally, calcium ions inhibit translation and reduce protein synthesis. However, the mechanism by which calcium regulates IMF deposition in muscle through translation remains largely unknown. In this study, we compared the ribosome profiles of the longissimus dorsi muscles of Duroc × Landrace × Large white pigs from the normal calcium (NC) group or calcium supplement (HC) group by Ribo-seq, and RNA-seq. By integrating multiple-omics analysis, we further discovered 437 genes that were transcriptionally unchanged but translationally altered and these genes were significantly enriched in the oxidative phosphorylation signaling pathway. Furthermore, experimental data showed that inhibiting the expression of COX10 and mtND4L increased triglyceride accumulation in C2C12 cells, providing new targets for intramuscular fat deposition. Finally, this work links dietary calcium, translation regulation and IMF deposition, providing a new strategy for both meat quality and economic performance within the pig industry.

Integrated analysis of transcriptome, translatome and proteome reveals insights into yellow catfish (Pelteobagrus fulvidraco) brain in response to hypoxia

Brain plays a central role in adapting to environmental changes and is highly sensitive to the oxygen level. Although previous studies investigated the molecular response of brain exposure to acute hypoxia in fish, the lack of studies at the translational level hinders further understanding of the regulatory mechanism response to hypoxia from multi-omics levels. Yellow catfish (Pelteobagrus fulvidraco) is an important freshwater aquaculture species; however, hypoxia severely restricts the sustainable development of its breeding industry. In the present study, the transcriptome, translatome, and proteome were integrated to study the global landscapes of yellow catfish brain response to hypoxia. The evidently increased amount of cerebral cortical cells with oedema and pyknotic nuclei has been observed in hypoxia group of yellow catfish. A total of 2750 genes were significantly changed at the translational level. Comparative transcriptional and translational analysis suggested the HIF-1 signaling pathway, autophagy and glycolysis/gluconeogenesis were up-regulated after hypoxia exposure. KEGG enrichment of translational efficiency (TE) differential genes suggested that the lysosome and autophagy were highly enriched. Our result showed that yellow catfish tends to inhibit the TE of genes by increasing the translation of uORFs to adapt to hypoxia. Correlation analysis showed that transcriptome and translatome exhibit higher correlation. In summary, this study demonstrated that hypoxia dysregulated the cerebral function of yellow catfish at the transcriptome, translatome, and proteome, which provides a better understanding of hypoxia adaptation in teleost.

An integrated metabolomic and transcriptomic analysis reveals the dynamic changes of key metabolites and flavor formation over Tieguanyin oolong tea production

To interpret the formation characteristic flavor during oolong tea manufacturing process, the dynamic changes of key flavor components in samples from various processing steps of Tieguanyin oolong tea production were investigated using widely-targeted metabolomic and the transcriptomic approaches. As a result, a total of 1078 metabolites were determined, of which 62 compounds were identified as biomarkers significantly changed over the manufacturing process. Quantitative determination of the total 50,343 transcripts showed 7480 of them were co-expressed different genes. Glutamic acid served as a critical metabolism hub and a signaling molecule for diverse stress responses. Additionally, the targeted quantification results showed that the contents of catechins and xanthine alkaloids in dried tea were dramatically decreased by 20.19% and 7.15% respectively than those in fresh leaves, which potentially contributed to the alleviation of astringent or bitter palates, promoting the characteristic mellow and rich flavor of Tieguanyin oolong tea.

Gene mining and genomics-assisted breeding empowered by the pangenome of tea plant Camellia sinensis

Tea is one of the world's oldest crops and is cultivated to produce beverages with various flavours. Despite advances in sequencing technologies, the genetic mechanisms underlying key agronomic traits of tea remain unclear. In this study, we present a high-quality pangenome of 22 elite cultivars, representing broad genetic diversity in the species. Our analysis reveals that a recent long terminal repeat burst contributed nearly 20% of gene copies, introducing functional genetic variants that affect phenotypes such as leaf colour. Our graphical pangenome improves the efficiency of genome-wide association studies and allows the identification of key genes controlling bud flush timing. We also identified strong correlations between allelic variants and flavour-related chemistries. These findings deepen our understanding of the genetic basis of tea quality and provide valuable genomic resources to facilitate its genomics-assisted breeding.

Translational Landscape of Medicago truncatula Seedlings under Salt Stress

The expression of plant genes under salt stress at the transcriptional level has been extensively studied. However, less attention has been paid to gene translation regulation under salt stress. In this study, Ribo-seq and RNA-seq analyses were conducted in Medicago truncatula seedlings grown under normal and salt stress conditions. The results showed that salt stress significantly altered the gene expression at the transcriptional and translational levels, with 2755 genes showing significant changes only at the translational level. Salt stress significantly inhibited the gene translation efficiency. Small ORFs (including uORFs in the 5'UTR, dORFs in 3'UTRs, and sORFs in lncRNAs) were identified throughout the genome of M. truncatula. The efficiency of gene translation was simultaneously regulated by the uORFs, dORFs, and miRNAs. In summary, our results provide valuable information about translatomic resources and new insights into plant responses to salt stress.

Comprehensive Metabolite Profile Uncovers the Bioactive Components, Antioxidant and Antibacterial Activities in Wild Tree Peony Leaves

Tree peonies (Paeonia Section Moutan)-including nine wild species, which belong to subsections Vaginatae and Delavayanae-are economically important plants with ornamental, nutritional, and medicinal applications. In this study, for the first time, we determined the bioactive components and antioxidant activities and antibacterial activities of the newly grown leaves of nine wild tree peony species (WTPS). A total of 276 bioactive components were identified through non-targeted metabolomics; more than 80% of the 276 metabolites identified are terpenoids and flavonoids. A total of 42 differential metabolites were quantitatively determined. The main differential metabolites were Paeoniflorin, Luteoloside, Hyperin, Apigenin-7-glucoside, Rhoifolin, and Cantharidin. Such a high terpenoid and flavonoid content of the leaf extracts renders them as species with strong antibacterial capacities, and most of the bacteria tested showed greater sensitivity derived from the members of subsection Vaginatae than those of subsection Delavayanae. All WTPS have significant antioxidant activity; this activity is attributed to high levels of the total phenolic content (TPC) and total flavonoid content (TFC), of which, among the nine WTPS, P. lutea has the strongest antioxidant capacity. Our results provided a theoretical basis for the in-deep application of tree peony leaves for food, medical, and pharmaceutical industries.

Dynamic Translational Landscape Revealed by Genome-Wide Ribosome Profiling under Drought and Heat Stress in Potato

The yield and quality of potatoes, an important staple crop, are seriously threatened by high temperature and drought stress. In order to deal with this adverse environment, plants have evolved a series of response mechanisms. However, the molecular mechanism of potato's response to environmental changes at the translational level is still unclear. In this study, we performed transcriptome- and ribosome-profiling assays with potato seedlings growing under normal, drought, and high-temperature conditions to reveal the dynamic translational landscapes for the first time. The translational efficiency was significantly affected by drought and heat stress in potato. A relatively high correlation (0.88 and 0.82 for drought and heat stress, respectively) of the fold changes of gene expression was observed between the transcriptional level and translational level globally based on the ribosome-profiling and RNA-seq data. However, only 41.58% and 27.69% of the different expressed genes were shared by transcription and translation in drought and heat stress, respectively, suggesting that the transcription or translation process can be changed independently. In total, the translational efficiency of 151 (83 and 68 for drought and heat, respectively) genes was significantly changed. In addition, sequence features, including GC content, sequence length, and normalized minimal free energy, significantly affected the translational efficiencies of genes. In addition, 28,490 upstream open reading frames (uORFs) were detected on 6463 genes, with an average of 4.4 uORFs per gene and a median length of 100 bp. These uORFs significantly affected the translational efficiency of downstream major open reading frames (mORFs). These results provide new information and directions for analyzing the molecular regulatory network of potato seedlings in response to drought and heat stress.

Ribosome footprint profiling enables elucidating the systemic regulation of fatty acid accumulation in Acer truncatum

BACKGROUND

The accumulation of fatty acids in plants covers a wide range of functions in plant physiology and thereby affects adaptations and characteristics of species. As the famous woody oilseed crop, Acer truncatum accumulates unsaturated fatty acids and could serve as the model to understand the regulation and trait formation in oil-accumulation crops. Here, we performed Ribosome footprint profiling combing with a multi-omics strategy towards vital time points during seed development, and finally constructed systematic profiling from transcription to proteomes. Additionally, we characterized the small open reading frames (ORFs) and revealed that the translational efficiencies of focused genes were highly influenced by their sequence features.

RESULTS

The comprehensive multi-omics analysis of lipid metabolism was conducted in A. truncatum. We applied the Ribo-seq and RNA-seq techniques, and the analyses of transcriptional and translational profiles of seeds collected at 85 and 115 DAF were compared. Key members of biosynthesis-related structural genes (LACS, FAD2, FAD3, and KCS) were characterized fully. More meaningfully, the regulators (MYB, ABI, bZIP, and Dof) were identified and revealed to affect lipid biosynthesis via post-translational regulations. The translational features results showed that translation efficiency tended to be lower for the genes with a translated uORF than for the genes with a non-translated uORF. They provide new insights into the global mechanisms underlying the developmental regulation of lipid metabolism.

CONCLUSIONS

We performed Ribosome footprint profiling combing with a multi-omics strategy in A. truncatum seed development, which provides an example of the use of Ribosome footprint profiling in deciphering the complex regulation network and will be useful for elucidating the metabolism of A. truncatum seed oil and the regulatory mechanisms.

The Impact of Different Withering Approaches on the Metabolism of Flavor Compounds in Oolong Tea Leaves

In this study, complementary metabolomic and proteomic analyses were conducted on the solar- and indoor-withered oolong tea leaves, and freshly plucked leaves as the control, for the purpose to reveal the mechanisms underlying the initial formation of some flavor determinants during the early stage of oolong tea processing. As a result, a total of 978 non-volatile compounds and 152 volatile compounds were identified, the flavonoids and several esters were differently accumulated in various tea samples. In total, 7048 proteins were qualitatively and quantitatively determined, the analysis on pathway enrichment showed that phenylpropanoid, flavonoid metabolisms, and protein processing in endoplasmic reticulum were the major pathways discriminating the different tea samples. The joint protein-metabolite analysis showed that the multiple stresses such as dehydration, heat, and ultra-violet irradiation occurred during the withering step induced the dynamic and distinct changes in the biochemical network in the treated leaves compared to fresh leaves. The significant decreases in flavonoids, xanthine alkaloids, and several amino acids contributed to the alleviation of bitter or astringent taste of withered leaves, although the decomposition of L-theanine resulted in the loss of umami flavor over the solar-withering step. Moreover, the fruity or floral aromas, especially volatile terpenoids and phenylpropanoids/benzenoids, were retained or accumulated in the solar withered leaves, potentially aiding the formation of a better characteristic flavor of oolong tea made by indoor withered tea leaves. Distinct effects of solar- and indoor-withering methods on the flavor determinant formation provide a novel insight into the relationship between the metabolite accumulation and flavor formation during the withering step of oolong tea production.

Multiomics analysis of the mechanisms behind flavonoid differences between purple and green tender shoots of Camellia sinensis var. assamica

Flavonoids are rich in tea plants (Camellia sinensis), and responsible for the flavor and healthful benefits of tea beverage. The anthocyanin levels in the purple tender shoots are higher than in the general green leaves of tea plant, which provide special materials to search metabolic mechanisms of flavonoid enrichment in plant. In this work, flavonoid differences between purple and green shoots from tea cultivars "Zijuan" (ZJ) and "Yunkang10" (YK-10) were investigated through metabolomic analysis, and mechanisms for their difference were surveyed by comparative transcriptomic and proteomic analysis. Levels of 34 flavonoids were different between ZJ and YK-10 shoots. Among them, 8 and 6 were marker metabolites in ZJ and YK-10, respectively. The differentially expressed genes (DEGs), differentially expressed proteins (DEPs), and different-level metabolites (DLMs) between ZJ and YK-10 were researched, respectively; and interactions including DEG-DLM, DEP-DLM, DEG-DEP, and DEG-DEP-DLM were analyzed; the contents of 18 characteristic flavonoids in tea leaves and expressions of 34 flavonoid metabolic genes were measured to verify the omics results. Integrated above analyses, a proposed model of flavonoids biosynthesis in tea shoots were established. The differential expression of the leucoanthocyanidin reductase (LAR), anthocyanidin synthase (ANS), anthocyanidin reductase (ANR), UDPG-flavonoid glucosyltransferase (UGT) 75L12 and 94P1 at gene level, and the ANS, ANR, and UGT78A15 at protein level, were closely associated with differences in flavonoids between ZJ and YK-10 shoot. Together, this study provides new information on the flavonoid accumulation mechanism in tea plant.

Integrative Transcriptomic and Phytohormonal Analyses Provide Insights into the Cold Injury Recovery Mechanisms of Tea Leaves

Tea plant is susceptible to low temperature, while the cold injury recovery mechanisms of tea leaves are still unclear. Windbreak has an effective and gradient range of protecting tea plants. Tea plants with increasing cold damage degree have varying recovery status accordingly, which are the ideal objects for investigating the cold injury recovery mechanisms of tea leaves. Here, we investigated the transcriptome and phytohormone profiles of tea leaves with different cold injury degrees in recovery (adjacent to the windbreak), and the levels of chlorophylls, malondialdehyde, major phytohormones as well as the activities of peroxidase (POD) and superoxide dismutase (SOD) were also measured. The results showed the content of total chlorophylls and the activity of POD in mature tea leaves gradually decreased with the distance to windbreak, while SOD showed the opposite. The major phytohormones were highly accumulated in the moderately cold-injured tea leaves. The biosynthesis of abscisic acid (ABA) was enhanced in the moderate cold damaged tea leaves, suggesting that ABA plays an important role in the cold response and resistance of tea plants. The transcriptomic result showed that the samples in different rows were well discriminated, and the pathways of plant-pathogen interaction and flavonoid biosynthesis were enriched based on KEGG analysis. WRKY, GRAS and NAC were the top classes of transcription factors differentially expressed in the different cold-injured tea leaves. Thus, windbreak is effective to protect adjacent tea plants from cold wave, and phytohormones importantly participate in the cold injury recovery of tea leaves.

Untargeted Metabolomics and Transcriptomics Reveal the Mechanism of Metabolite Differences in Spring Tender Shoots of Tea Plants of Different Ages

The metabolites in the tender shoots of the tea plant are the material basis for the determination of tea quality. The composition and abundance of these metabolites are affected by many key factors, and the tea plant’s age is one of them. However, the effect of plant age on the tender shoot metabolites of tea cultivars of different genotypes is poorly understood. Therefore, we used a combination of untargeted metabolomics and transcriptomics to analyze the differential mechanism behind the differences in the metabolites of the spring tender shoots of 7- and 40-year-old tea plants of two tea cultivars of different genotypes. We found that plant age could significantly change the metabolites in the spring tender shoots of tea plants and that flavonoids, and amino acids and their derivatives, were predominant among the differential metabolites. The quantities of most flavonoids in the aged tea plants of different genotypes were upregulated, which was caused by the upregulated expression of differential genes in the flavonoid biosynthesis pathway. We further discovered that 11 key structural genes play key regulatory roles in the changes in the flavonoid contents of tea plants of different plant ages. However, the influence of plant age on amino acids and their derivatives might be cultivar-specific. By characterizing and evaluating the quality-related metabolites of tea cultivars of two different genotypes at different plant ages, we found that whether an old tea plant (40 years old) can produce high-quality tea is related to the genotype of the tea plant.

Quantitative trait loci mapping for free amino acid content using an albino population and SNP markers provides insight into the genetic improvement of tea plants

Free amino acids are one of the main chemical components in tea, and they contribute to the pleasant flavor, function, and quality of tea, notably the level of theanine. Here, a high-density genetic map was constructed to characterize quantitative trait loci (QTLs) for free amino acid content. A total of 2688 polymorphic SNP markers were obtained using genotyping-by-sequencing (GBS) based on 198 individuals derived from a pseudotestcross population of "Longjing 43" × "Baijiguan", which are elite and albino tea cultivars, respectively. The 1846.32 cM high-density map with an average interval of 0.69 cM was successfully divided into 15 linkage groups (LGs) ranging from 93.41 cM to 171.28 cM. Furthermore, a total of 4 QTLs related to free amino acid content (theanine, glutamate, glutamine, aspartic acid and arginine) identified over two years were mapped to LG03, LG06, LG11 and LG14. The phenotypic variation explained by these QTLs ranged from 11.8% to 23.7%, with an LOD score from 3.56 to 7.7. Furthermore, several important amino acid metabolic pathways were enriched based on the upregulated differentially expressed genes (DEGs) among the offspring. These results will be essential for fine mapping genes involved in amino acid pathways and diversity, thereby providing a promising avenue for the genetic improvement of tea plants.

Multiomics Analysis Reveals New Insights into the Apple Fruit Quality Decline under High Nitrogen Conditions

Excessive application of nitrogen (N) fertilizer is common in Chinese apple production. High N reduced the contents of soluble sugar and total flavonoids by 16.05 and 19.01%, respectively, resulting in poor fruit quality. Moreover, high N increased the total N and decreased the total C and C/N ratio of apple fruits. On the basis of the transcriptomic, proteomic, and metabolomic analyses, the global network was revealed. High N inhibited the accumulation of carbohydrates (sucrose, glucose, and trehalose) and flavonoids (rhamnetin-3-O-rutinoside, rutin, and trihydroxyisoflavone-7-O-galactoside) in fruits, and more C skeletons were used to synthesize amino acids and their derivatives (especially low C/N ratio, e.g., arginine) to be transferred to N metabolism. This study revealed new insights into the decline in soluble sugar and flavonoids caused by high N, and hub genes (MD07G1172700, MD05G1222800, MD16G1227200, MD01G1174400, and MD02G1207200) and hub proteins (PFK, gapN, and HK) were obtained.

Chromatin accessibility and translational landscapes of tea plants under chilling stress

Plants have evolved regulatory mechanisms at multiple levels to regulate gene expression in order to improve their cold adaptability. However, limited information is available regarding the stress response at the chromatin and translational levels. Here, we characterize the chromatin accessibility, transcriptional, and translational landscapes of tea plants in vivo under chilling stress for the first time. Chilling stress significantly affected both the transcription and translation levels as well as the translation efficiency of tea plants. A total of 3010 genes that underwent rapid and independent translation under chilling stress were observed, and they were significantly enriched in the photosynthesis-antenna protein and phenylpropanoid biosynthesis pathways. A set of genes that were significantly responsive to cold at the transcription and translation levels, including four (+)-neomenthol dehydrogenases (MNDs) and two (E)-nerolidol synthases (NESs) arranged in tandem on the chromosomes, were also found. We detected potential upstream open reading frames (uORFs) on 3082 genes and found that tea plants may inhibit the overall expression of genes by enhancing the translation of uORFs under chilling stress. In addition, we identified distal transposase hypersensitive sites (THSs) and proximal THSs and constructed a transcriptional regulatory network for tea plants under chilling stress. We also identified 13 high-confidence transcription factors (TFs) that may play a crucial role in cold regulation. These results provide valuable information regarding the potential transcriptional regulatory network in plants and help to clarify how plants exhibit flexible responses to chilling stress.

Cotton transcriptome analysis reveals novel biological pathways that eliminate reactive oxygen species (ROS) under sodium bicarbonate (NaHCO3) alkaline stress

Alkaline stress is one of the abiotic stresses limiting cotton production. Though RNA-Seq analyses, have been conducted to investigate genome-wide gene expression in response to alkaline stress in plants, the response of sodium bicarbonate (NaHCO₃) stress-related genes in cotton has not been reported. To explore the mechanisms of cotton response to this alkaline stress, we used next-generation sequencing (NGS) technology to study transcriptional changes of cotton under NaHCO₃ alkaline stress. A total of 18,230 and 11,177 differentially expressed genes (DEGs) were identified in cotton roots and leaves, respectively. Gene ontology (GO) analysis indicated the enrichment of DEGs involved in various stimuli or stress responses. Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis showed that DEGs associated with plant hormone signal transduction, amino acid biosynthesis, and biosynthesis of secondary metabolites were regulated in response to the NaHCO₃ stress. We further analyzed genes enriched in secondary metabolic pathways and found that secondary metabolites were regulated to eliminate the reactive oxygen species (ROS) and improve the cotton tolerance to the NaHCO₃ stress. In this study, we learned that the toxic effect of NaHCO₃ was more profound than that of NaOH at the same pH. Thus, Na⁺, HCO₃^- and pH had a great impact on the growth of cotton plant. The novel biological pathways and candidate genes for the cotton tolerance to NaHCO₃ stress identified from the study would be useful in the genetic improvement of the alkaline tolerance in cotton.

Transcriptional analysis reveals key insights into seasonal induced anthocyanin degradation and leaf color transition in purple tea (Camellia sinensis (L.) O. Kuntze)

Purple-tea, an anthocyanin rich cultivar has recently gained popularity due to its health benefits and captivating leaf appearance. However, the sustainability of purple pigmentation and anthocyanin content during production period is hampered by seasonal variation. To understand seasonal dependent anthocyanin pigmentation in purple tea, global transcriptional and anthocyanin profiling was carried out in tea shoots with two leaves and a bud harvested during in early (reddish purple: S1_RP), main (dark gray purple: S2_GP) and backend flush (moderately olive green: S3_G) seasons. Of the three seasons, maximum accumulation of total anthocyanin content was recorded in S2_GP, while least amount was recorded during S3_G. Reference based transcriptome assembly of 412 million quality reads resulted into 71,349 non-redundant transcripts with 6081 significant differentially expressed genes. Interestingly, key DEGs involved in anthocyanin biosynthesis [PAL, 4CL, F3H, DFR and UGT/UFGT], vacuolar trafficking [ABC, MATE and GST] transcriptional regulation [MYB, NAC, bHLH, WRKY and HMG] and Abscisic acid signaling pathway [PYL and PP2C] were significantly upregulated in S2_GP. Conversely, DEGs associated with anthocyanin degradation [Prx and lac], repressor TFs and key components of auxin and ethylene signaling pathways [ARF, AUX/IAA/SAUR, ETR, ERF, EBF1/2] exhibited significant upregulation in S3_G, correlating positively with reduced anthocyanin content and purple coloration. The present study for the first-time elucidated genome-wide transcriptional insights and hypothesized the involvement of anthocyanin biosynthesis activators/repressor and anthocyanin degrading genes via peroxidases and laccases during seasonal induced leaf color transition in purple tea. Futuristically, key candidate gene(s) identified here can be used for genetic engineering and molecular breeding of seasonal independent anthocyanin-rich tea cultivars.