Establishment and validation of a callus tissue transformation system for German chamomile (Matricaria chamomilla L.).

BACKGROUND: German chamomile (Matricaria chamomilla L.) is an important medicinal plant, and the essential oils in the flowers have various biological activities. Genetic transformation systems are important for plant quality improvement and molecular research. To the best of our knowledge, a genetic transformation system has not yet been reported for German chamomile. RESULTS: In this study, we developed Agrobacterium-mediated transformation protocols for German chamomile callus tissues. This involved optimizing key parameters, such as hygromycin and cefotaxime concentrations, bacterial density, and infection and co-culture durations. We also performed gas chromatography-mass spectrometry analysis to identify volatile compounds in non-transgenic and transgenic callus and hairy root tissues. Furthermore, to compare and verify the callus transformation system of German chamomile, we transferred McFPS to the hairy roots of German chamomile. The results showed that the optimal conditions for Agrobacterium-mediated callus tissue transformation were as follows: explant, petiole; cefotaxime concentration, 300 mg/L; hygromycin concentration, 10 mg/L; and bacterial solution concentration, OD600 = 0.6; callus transformation efficiency was the highest when the co-culture time was 3 days. CONCLUSIONS: Establishment of a high-efficiency callus transformation system will lay the foundation for gene function identification in German chamomile.

Biosynthesis of α-Bisabolol by Farnesyl Diphosphate Synthase and α-Bisabolol Synthase and Their Related Transcription Factors in Matricaria recutita L.

The essential oil of German chamomile (Matricaria recutita L.) is widely used in food, cosmetics, and the pharmaceutical industry. α-Bisabolol is the main active substance in German chamomile. Farnesyl diphosphate synthase (FPS) and α-bisabolol synthase (BBS) are key enzymes related to the α-bisabolol biosynthesis pathway. However, little is known about the α-bisabolol biosynthesis pathway in German chamomile, especially the transcription factors (TFs) related to the regulation of α-bisabolol synthesis. In this study, we identified MrFPS and MrBBS and investigated their functions by prokaryotic expression and expression in hairy root cells of German chamomile. The results suggest that MrFPS is the key enzyme in the production of sesquiterpenoids, and MrBBS catalyzes the reaction that produces α-bisabolol. Subcellular localization analysis showed that both MrFPS and MrBBS proteins were located in the cytosol. The expression levels of both MrFPS and MrBBS were highest in the extension period of ray florets. Furthermore, we cloned and analyzed the promoters of MrFPS and MrBBS. A large number of cis-acting elements related to light responsiveness, hormone response elements, and cis-regulatory elements that serve as putative binding sites for specific TFs in response to various biotic and abiotic stresses were identified. We identified and studied TFs related to MrFPS and MrBBS, including WRKY, AP2, and MYB. Our findings reveal the biosynthesis and regulation of α-bisabolol in German chamomile and provide novel insights for the production of α-bisabolol using synthetic biology methods.

Clinical Benefits of Golden-Antrodia Camphorata Containing Antroquinonol in Liver Protection and Liver Fat Reduction After Alcoholic Hepatitis.

Objective: It has been reported that antroquinonol extracted from Golden-Antrodia camphorate exerts protective effects on liver function both in vitro and in vivo. However, the protective effects of Golden-Antrodia camphorata on liver function have not been fully investigated in human clinical studies. Therefore, the present study aimed to evaluate the beneficial effects of Golden-Antrodia camphorata on hepatic function after alcohol consumption in human subjects. Methods: A total of 80 participants with increased γ-glutamyl transferase levels (60-180 U/L) were enrolled in the current study and were randomly divided into two groups. Participants in the first group were orally administrated with 300 mg/day Golden-Antrodia camphorata (tablets), while those in the second group received placebo tablets for 12 weeks. Biochemical routine blood tests were performed at 6 and 12 weeks following the first administration. Results: At 12 weeks post the first Golden-Antrodia camphorata administration, the serum levels of aspartate aminotransferase (AST; p < 0.0001), alanine aminotransferase (ALT; p = 0.0002) and triglyceride (p = 0.0158) were notably declined in the Golden-Antrodia camphorata treatment group compared with the placebo group. No clinically significant differences were observed between the Golden-Antrodia camphorata treatment and placebo groups in terms of general safety parameters. Conclusion: A statistically significant difference was obtained in the serum levels of AST, ALT and triglycerides between the Golden-Antrodia camphorata and placebo groups. However, no clinical significance was observed in any of the safety parameters examined. Overall, these findings indicated that treatment with Golden-Antrodia camphorata exerted protective effects on liver function.

Integrated multi-omics analysis provides insights into genome evolution and phosphorus deficiency adaptation in pigeonpea (Cajanus cajan).

Pigeonpea (Cajanus cajan) is an important legume food crop and plays a crucial role in a secure food supply in many developing countries. Several previous studies have suggested that pigeonpea has great potential for phosphorus (P) deficiency tolerance, but little is known about the underlying mechanism. In this study, the physiological and molecular responses of pigeonpea roots to phosphate (Pi) starvation were investigated through integrating phenotypic, genomic, transcriptomic, metabolomic, and lipidomic analyses. The results showed that low-Pi treatment increased total root length, root surface area, and root acid phosphatase activity, and promoted the secretion of organic acids (e.g. citric acids, piscidic acids, and protocatechuic acids) and the degradation of phospholipids and other P-containing metabolites in the roots of pigeonpea. Consistent with the morphological, physiological, and biochemical changes, a large number of genes involved in these Pi-starvation responses were significantly upregulated in Pi-deficient pigeonpea roots. Among these Pi-starvation response genes upregulated by low-Pi treatment, four gene families were expanded through recent tandem duplication in the pigeonpea genome, namely phosphate transporter 1 (PHT1), phosphoethanolamine/phosphocholine phosphatase (PECP), fasciclin-like arabinogalactan protein (FLA), and glutamate decarboxylase (GAD). These gene families may be associated with Pi uptake from the soil, phospholipid recycling, root morphological remodeling, and regulation of organic acid exudation. Taken together, our results suggest that pigeonpea employs complex Pi-starvation responses to strengthen P acquisition and utilization during low-Pi stress. This study provides new insights into the genome evolution and P deficiency adaptation mechanism of pigeonpea.

Insights into acylation mechanisms: co-expression of serine carboxypeptidase-like acyltransferases and their non-catalytic companion paralogs.

Serine carboxypeptidase-like acyltransferases (SCPL-ATs) play a vital role in the diversification of plant metabolites. Galloylated flavan-3-ols highly accumulate in tea (Camellia sinensis), grape (Vitis vinifera), and persimmon (Diospyros kaki). To date, the biosynthetic mechanism of these compounds remains unknown. Herein, we report that two SCPL-AT paralogs are involved in galloylation of flavan-3-ols: CsSCPL4, which contains the conserved catalytic triad S-D-H, and CsSCPL5, which has the alternative triad T-D-Y. Integrated data from transgenic plants, recombinant enzymes, and gene mutations showed that CsSCPL4 is a catalytic acyltransferase, while CsSCPL5 is a non-catalytic companion paralog (NCCP). Co-expression of CsSCPL4 and CsSCPL5 is likely responsible for the galloylation. Furthermore, pull-down and co-immunoprecipitation assays showed that CsSCPL4 and CsSCPL5 interact, increasing protein stability and promoting post-translational processing. Moreover, phylogenetic analyses revealed that their homologs co-exist in galloylated flavan-3-ol- or hydrolyzable tannin-rich plant species. Enzymatic assays further revealed the necessity of co-expression of those homologs for acyltransferase activity. Evolution analysis revealed that the mutations of the CsSCPL5 catalytic residues may have taken place about 10 million years ago. These findings show that the co-expression of SCPL-ATs and their NCCPs contributes to the acylation of flavan-3-ols in the plant kingdom.

Cardiac Rehabilitation of Patients with Acute ST-Elevation Myocardial Infarction Undergoing Primary Percutaneous Coronary Intervention in a Han Population in Northern China: A Prospective Cohort Study.

BACKGROUND: Cardiac rehabilitation (CR) has been associated with improved cardiac function in cardiovascular diseases. Our aim was to explore the factors associated with cardiac function and CR. METHODS: This prospective cohort study had 473 STEMI patients admitted for primary percutaneous coronary intervention (PCI) who were divided into a CR group (group A, n = 104) and a non-CR group (group B, n = 369) based on whether they could complete CR. Patients' clinical features, such as age, hyperlipidemia, family history of premature coronary heart disease (FHPCHD), smoking history, body mass index (BMI, kg/m2), number of diseased vessels, arrhythmia during PCI, N-terminal pro-B-type natriuretic peptide (NT-proBNP), left ventricular ejection fraction (LVEF) and major adverse cardiac events (MACEs) at 6 months after PCI were compared. Then, the clinical characteristics of patients were further analyzed according to those with MACEs (n = 78) and those without MACEs (n = 395). RESULTS: After CR of 6 months, NT-proBNP levels (p = 0.027), 6-MWD (meter, P = 0.000), LVEF (P = 0.000) were significantly improved in group A compared to group B, but not for SMWA (P = 0.875). Multivariate analysis indicated that even though patients in group A (OR 3.06, 95% CI 1.132-8.274, p = 0.03) have a higher incidence of hyperlipidemia, their MACEs (OR 0.191, 95% CI 0.038-0.961, p = 0.05) at 6 months were significantly lower than in group B, mainly because the average patient low age (<65 years, OR 0.917, 95% CI 0.859-0.979, p = 0.01) and significant improvement of 6-MWD (OR 7.999, 95% CI 4.342-14.737, P = 0.00) and the LVEF at 6 months (OR 1.112, 95% CI 1.072-1.154, p = 0.00). Further analysis based on the MACES outcomes showed that there were 6 factors associated with the occurrence of MACEs, they were age >65 years (OR 1.032, 95% CI 1.009-1.009, p = 0.007), smoking history (OR 0.485, 95% CI 0.238-0.989, p = 0.046), education level (OR 2.646, 95% CI 1.370-5.108, p = 0.004), 6-MWD (OR 1.688, 95% CI 1.104-2.811, p = 0.044), LVEF (OR 0.958, 95% CI 0.926-0.991, p = 0.013) and CR (OR 6.271, 95% CI 2.236-17.590, p = 0.000). CONCLUSION: CR, including exercise rehabilitation, is a beneficial option to reduce MACEs in STEMI patients treated with primary PCI.

Transcriptional regulation of amino acid metabolism in response to nitrogen deficiency and nitrogen forms in tea plant root (Camellia sinensis L.).

Free amino acids, including theanine, glutamine and glutamate, contribute greatly to the pleasant taste and multiple health benefits of tea. Amino acids in tea plants are mainly synthesized in roots and transported to new shoots, which are significantly affected by nitrogen (N) level and forms. However, the regulatory amino acid metabolism genes have not been systemically identified in tea plants. Here, we investigated the dynamic changes of free amino acid contents in response to N deficiency and forms in tea plant roots, and systemically identified the genes associated amino acid contents in individual metabolism pathways. Our results showed that glutamate-derived amino acids are the most dynamic in response to various forms of N and N deficiency. We then performed transcriptomic analyses of roots treated with N deficiency and various forms of N, and differentially expressed amino acid metabolic genes in each pathway were identified. The analyses on expression patterns and transcriptional responses of metabolic genes to N treatments provided novel insights for the molecular basis of high accumulation of theanine in tea plant root. These analyses also identified potential regulatory genes in dynamic amino acid metabolism in tea plant root. Furthermore, our findings indicated that the dynamic expression levels of CsGDH, CsAlaDC, CsAspAT, CsSDH, CsPAL, CsSHMT were highly correlated with changes of amino acid contents in their corresponding pathways. Herein, this study provides comprehensive insights into transcriptional regulation of amino acid metabolism in response to nitrogen deficiency and nitrogen forms in tea plant root.

Cloning and functional analysis of three aphid alarm pheromone genes from German chamomile (Matricaria chamomilla L.).

German chamomile (Matricaria chamomilla L.) is one of the most ancient medicinal species in the world and terpenoids from their flowers have important medicinal value. We cloned three sesquiterpene synthase genes, McGDS1, McGDS2 and McGDS3, and performed sequence alignment and phylogenetic analysis. The encoded proteins possess three conserved structural features: an RRxxxxxxxxW motif, an RxR motif, and a DDxxD motif. McGDS1, McGDS2 and McGDS3 were confirmed to be (E)-farnesene synthase, germacrene D synthase, and germacrene A synthase, respectively. Subcellular localization revealed diffuse GFP reporter-gene signals in the cytoplasm and nucleus. qPCR indicated that McGDS1, McGDS2 and McGDS3, were more highly expressed in young flowers than in old flowers and the expression was highly correlated with amounts of the end-product essential oils ((E)-ß-farnesene, germacrene D and ß-elemene), with coefficients of 0.76, 0.83 and 0.68, respectively. We also established a transformation system for chamomile hairy roots. The overexpression of McGDS1, McGDS2 and McGDS3 resulted in γ-muurolene accumulation in hairy roots. The activity of three aphid alarm pheromones here forms the molecular basis for the study of the biosynthesis and regulation of volatile terpenes. Transformation of chamomile hairy roots provides a simple system in which to study terpene biosynthesis in chamomile.

The role of extracellular vesicles in prostate cancer with clinical applications.

In mammalian cells, extracellular vesicles (EVs) derived from the endosomal system carry many different kinds of bioactive molecule to deliver to recipient cells in a paracrine or endocrine manner. EVs can mediate local and systemic intercellular communications, including reeducating stromal cells, remodeling the architecture of the tumor microenvironment, modulating cancer metabolism and metastases, or even conferring drug resistance. Because the molecular and functional characteristics of prostate cancer (PCa) evolve over time, the bioactive molecule profiles/signatures of tumor-derived EVs (TDEs) reflect the real-time status of cancer cells. TDEs appear to be valuable diagnostic and prognostic biomarkers as well as potential therapeutic vehicles, suggesting their essential role in precision medicine of disease management. We summarized critical aspects of TDEs in PCa and discussed their potential clinical applications.

Phytochemical and comparative transcriptome analyses reveal different regulatory mechanisms in the terpenoid biosynthesis pathways between Matricaria recutita L. and Chamaemelum nobile L.

BACKGROUND: Matricaria recutita (German chamomile) and Chamaemelum nobile (Roman chamomile) belong to the botanical family Asteraceae. These two herbs are not only morphologically distinguishable, but their secondary metabolites - especially the essential oils present in flowers are also different, especially the terpenoids. The aim of this project was to preliminarily identify regulatory mechanisms in the terpenoid biosynthetic pathways that differ between German and Roman chamomile by performing comparative transcriptomic and metabolomic analyses. RESULTS: We determined the content of essential oils in disk florets and ray florets in these two chamomile species, and found that the terpenoid content in flowers of German chamomile is greater than that of Roman chamomile. In addition, a comparative RNA-seq analysis of German and Roman chamomile showed that 54% of genes shared > 75% sequence identity between the two species. In particular, more highly expressed DEGs (differentially expressed genes) and TF (transcription factor) genes, different regulation of CYPs (cytochrome P450 enzymes), and rapid evolution of downstream genes in the terpenoid biosynthetic pathway of German chamomile could be the main reasons to explain the differences in the types and levels of terpenoid compounds in these two species. In addition, a phylogenetic tree constructed from single copy genes showed that German chamomile and Roman chamomile are closely related to Chrysanthemum nankingense. CONCLUSION: This work provides the first insights into terpenoid biosynthesis in two species of chamomile. The candidate unigenes related to terpenoid biosynthesis will be important in molecular breeding approaches to modulate the essential oil composition of Matricaria recutita and Chamaemelum nobile.

Wnt5b integrates Fak1a to mediate gastrulation cell movements via Rac1 and Cdc42.

Focal adhesion kinase (FAK) mediates vital cellular pathways during development. Despite its necessity, how FAK regulates and integrates with other signals during early embryogenesis remains poorly understood. We found that the loss of Fak1a impaired epiboly, convergent extension and hypoblast cell migration in zebrafish embryos. We also observed a clear disturbance in cortical actin at the blastoderm margin and distribution of yolk syncytial nuclei. In addition, we investigated a possible link between Fak1a and a well-known gastrulation regulator, Wnt5b, and revealed that the overexpression of fak1a or wnt5b could cross-rescue convergence defects induced by a wnt5b or fak1a antisense morpholino (MO), respectively. Wnt5b and Fak1a were shown to converge in regulating Rac1 and Cdc42, which could synergistically rescue wnt5b and fak1a morphant phenotypes. Furthermore, we generated several alleles of fak1a mutants using CRISPR/Cas9, but those mutants only revealed mild gastrulation defects. However, injection of a subthreshold level of the wnt5b MO induced severe gastrulation defects in fak1a mutants, which suggested that the upregulated expression of wnt5b might complement the loss of Fak1a. Collectively, we demonstrated that a functional interaction between Wnt and FAK signalling mediates gastrulation cell movements via the possible regulation of Rac1 and Cdc42 and subsequent actin dynamics.

Analysis of terpenoid biosynthesis pathways in German chamomile (Matricaria recutita) and Roman chamomile (Chamaemelum nobile) based on co-expression networks.

German chamomile and Roman chamomile are the two most widely known chamomile species due to the medicinal properties of volatile compounds from their flowers. We determined the volatile compound content of different organs of these two chamomiles, and found that main volatile compounds in German chamomile were terpenoids and those in Roman chamomile were esters. Furthermore, 24 tissues from two chamomiles were sequenced and analyzed by gene co-expression network. The results showed higher terpene synthase expression levels and more modules correlated with sesquiterpenoids in German chamomile, which may explain its high sesquiterpenoid content. In both chamomiles, unigenes in volatile compound-correlated modules were significantly enriched in pathways related to plant-pathogen interactions and circadian rhythm, demonstrating that volatile compounds of chamomiles are influenced by these factors. There were ten times more unigenes related to plant-pathogen interactions in German chamomile than in Roman chamomile, which indicates German chamomile has higher resistance to pathogens.

The paracrine induction of prostate cancer progression by caveolin-1.

A subpopulation of cancer stem cells (CSCs) plays a critical role of cancer progression, recurrence, and therapeutic resistance. Many studies have indicated that castration-resistant prostate cancer (CRPC) is associated with stem cell phenotypes, which could further promote neuroendocrine transdifferentiation. Although only a small subset of genetically pre-programmed cells in each organ has stem cell capability, CSCs appear to be inducible among a heterogeneous cancer cell population. However, the inductive mechanism(s) leading to the emergence of these CSCs are not fully understood in CRPC. Tumor cells actively produce, release, and utilize exosomes to promote cancer development and metastasis, cancer immune evasion as well as chemotherapeutic resistance; the impact of tumor-derived exosomes (TDE) and its cargo on prostate cancer (PCa) development is still unclear. In this study, we demonstrate that the presence of Cav-1 in TDE acts as a potent driver to induce CSC phenotypes and epithelial-mesenchymal transition in PCa undergoing neuroendocrine differentiation through NFκB signaling pathway. Furthermore, Cav-1 in mCRPC-derived exosomes is capable of inducing radio- and chemo-resistance in recipient cells. Collectively, these data support Cav-1 as a critical driver for mCRPC progression.

Comparative Transcriptomic Analysis Reveals Regulatory Mechanisms of Theanine Synthesis in Tea (Camellia sinensis) and Oil Tea (Camellia oleifera) Plants.

Tea provides a rich taste and has healthy properties due to its variety of bioactive compounds, such as theanine, catechins, and caffeine. Theanine is the most abundant free amino acid (40%-70%) in tea leaves. Key genes related to theanine biosynthesis have been studied, but relatively little is known about the regulatory mechanisms of theanine accumulation in tea leaves. Herein, we analyzed theanine content in tea (Camellia sinensis) and oil tea (Camellia oleifera) and found it to be higher in the roots than in other tissues in both species. The theanine content was significantly higher in tea than oil tea. To explore the regulatory mechanisms of theanine accumulation, we identified genes involved in theanine biosynthesis by RNA-Seq analysis and compared theanine-related modules. Moreover, we cloned theanine synthase (TS) promoters from tea and oil tea plants and found that a difference in TS expression and cis-acting elements may explain the difference in theanine accumulation between the two species. These data provide an important resource for regulatory mechanisms of theanine accumulation in tea plants.

Activation of Aryl Hydrocarbon Receptor by Kynurenine Impairs Progression and Metastasis of Neuroblastoma.

Neuroblastoma is the most common malignant disease of infancy, and amplification of the MYCN oncogene is closely associated with poor prognosis. Recently, expression of MYCN was shown to be inversely correlated with aryl hydrocarbon receptor (AHR) expression in neuroblastoma, and overexpression of AHR downregulated MYCN expression, promoting cell differentiation. Therefore, we further investigated the potential of AHR to serve as a prognostic indicator or a therapeutic target in neuroblastoma. First, the clinical significance of AHR in neuroblastoma was examined. Positive AHR immunostaining strongly correlated with differentiated histology of neuroblastoma and predicted better survival for patients. The mouse xenograft model showed that overexpression of AHR significantly suppressed neuroblastoma tumor growth. In addition, activation of AHR by the endogenous ligand kynurenine inhibited cell proliferation and promoted cell differentiation in vitro and in vivo. kynurenine treatment also upregulated the expression of KISS1, a tumor metastasis suppressor, and attenuated metastasis in the xenograft model. Finally, analysis of KISS1 levels in neuroblastoma patient tumors using the R2: Genomics Analysis and Visualization Platform revealed that KISS1 expression positively correlated with AHR, and high KISS1 expression predicted better survival for patients. In conclusion, our results indicate that AHR is a novel prognostic biomarker for neuroblastoma, and that overexpression or activation of AHR offers a new therapeutic possibility for patients with neuroblastoma. SIGNIFICANCE: These findings show that AHR may function as a tumor suppressor in childhood neuroblastoma, potentially influencing the aetiologic and therapeutic targeting of the disease.

The tea plant reference genome and improved gene annotation using long-read and paired-end sequencing data.

Tea is a globally consumed non-alcohol beverage with great economic importance. However, lack of the reference genome has largely hampered the utilization of precious tea plant genetic resources towards breeding. To address this issue, we previously generated a high-quality reference genome of tea plant using Illumina and PacBio sequencing technology, which produced a total of 2,124 Gb short and 125 Gb long read data, respectively. A hybrid strategy was employed to assemble the tea genome that has been publicly released. We here described the data framework used to generate, annotate and validate the genome assembly. Besides, we re-predicted the protein-coding genes and annotated their putative functions using more comprehensive omics datasets with improved training models. We reassessed the assembly and annotation quality using the latest version of BUSCO. These data can be utilized to develop new methodologies/tools for better assembly of complex genomes, aid in finding of novel genes, variations and evolutionary clues associated with tea quality, thus help to breed new varieties with high yield and better quality in the future.

Basics and applications of tumor-derived extracellular vesicles.

Extracellular vesicle (EV)-mediated intercellular communication acts as a critical culprit in cancer development. The selective packaging of oncogenic molecules renders tumor-derived EVs capable of altering the tumor microenvironment and thereby modulating cancer developments that may contribute to drug resistance and cancer recurrence. Moreover, the molecular and functional characteristics of cancer through its development and posttreatment evolve over time. Tumor-derived EVs are profoundly involved in this process and can, therefore, provide valuable real-time information to reflect dynamic changes occurring within the body. Because they bear unique molecular profiles or signatures, tumor-derived EVs have been highlighted as valuable diagnostic and predictive biomarkers as well as novel therapeutic targets. In addition, the use of an advanced EV-based drug delivery system for cancer therapeutics has recently been emphasized in both basic and clinical studies. In this review, we highlight comprehensive aspects of tumor-derived EVs in oncogenic processes and their potential clinical applications.

Grb7, a Critical Mediator of EGFR/ErbB Signaling, in Cancer Development and as a Potential Therapeutic Target.

The partner of activated epidermal growth factor receptor (EGFR), growth factor receptor bound protein-7 (Grb7), a functionally multidomain adaptor protein, has been demonstrated to be a pivotal regulator for varied physiological and pathological processes by interacting with phospho-tyrosine-related signaling molecules to affect the transmission through a number of signaling pathways. In particular, critical roles of Grb7 in erythroblastic leukemia viral oncogene homolog (ERBB) family-mediated cancer development and malignancy have been intensively evaluated. The overexpression of Grb7 or the coamplification/cooverexpression of Grb7 and members of the ERBB family play essential roles in advanced human cancers and are associated with decreased survival and recurrence of cancers, emphasizing Grb7's value as a prognostic marker and a therapeutic target. Peptide inhibitors of Grb7 are being tested in preclinical trials for their possible therapeutic effects. Here, we review the molecular, functional, and clinical aspects of Grb7 in ERBB family-mediated cancer development and malignancy with the aim to reveal alternative and effective therapeutic strategies.

Global dissection of alternative splicing uncovers transcriptional diversity in tissues and associates with the flavonoid pathway in tea plant (Camellia sinensis).

BACKGROUND: Alternative splicing (AS) regulates mRNA at the post-transcriptional level to change gene function in organisms. However, little is known about the AS and its roles in tea plant (Camellia sinensis), widely cultivated for making a popular beverage tea. RESULTS: In our study, the AS landscape and dynamics were characterized in eight tissues (bud, young leaf, summer mature leaf, winter old leaf, stem, root, flower, fruit) of tea plant by Illumina RNA-Seq and confirmed by Iso-Seq. The most abundant AS (~ 20%) was intron retention and involved in RNA processes. The some alternative splicings were found to be tissue specific in stem and root etc. Thirteen co-expressed modules of AS transcripts were identified, which revealed a similar pattern between the bud and young leaves as well as a distinct pattern between seasons. AS events of structural genes including anthocyanidin reductase and MYB transcription factors were involved in biosynthesis of flavonoid, especially in vegetative tissues. The AS isoforms rather than the full-length ones were the major transcripts involved in flavonoid synthesis pathway, and is positively correlated with the catechins content conferring the tea taste. We propose that the AS is an important functional mechanism in regulating flavonoid metabolites. CONCLUSION: Our study provides the insight into the AS events underlying tea plant's uniquely different developmental process and highlights the important contribution and efficacy of alternative splicing regulatory function to biosynthesis of flavonoids.