Gene co-expression network analysis reveals coordinated regulation of three characteristic secondary biosynthetic pathways in tea plant (Camellia sinensis)

A method for mining condition-specific co-expressed genes in Camellia sinensis based on k-means clustering

BACKGROUND

As one of the world's most important beverage crops, tea plants (Camellia sinensis) are renowned for their unique flavors and numerous beneficial secondary metabolites, attracting researchers to investigate the formation of tea quality. With the increasing availability of transcriptome data on tea plants in public databases, conducting large-scale co-expression analyses has become feasible to meet the demand for functional characterization of tea plant genes. However, as the multidimensional noise increases, larger-scale co-expression analyses are not always effective. Analyzing a subset of samples generated by effectively downsampling and reorganizing the global sample set often leads to more accurate results in co-expression analysis. Meanwhile, global-based co-expression analyses are more likely to overlook condition-specific gene interactions, which may be more important and worthy of exploration and research.

RESULTS

Here, we employed the k-means clustering method to organize and classify the global samples of tea plants, resulting in clustered samples. Metadata annotations were then performed on these clustered samples to determine the "conditions" represented by each cluster. Subsequently, we conducted gene co-expression network analysis (WGCNA) separately on the global samples and the clustered samples, resulting in global modules and cluster-specific modules. Comparative analyses of global modules and cluster-specific modules have demonstrated that cluster-specific modules exhibit higher accuracy in co-expression analysis. To measure the degree of condition specificity of genes within condition-specific clusters, we introduced the correlation difference value (CDV). By incorporating the CDV into co-expression analyses, we can assess the condition specificity of genes. This approach proved instrumental in identifying a series of high CDV transcription factor encoding genes upregulated during sustained cold treatment in Camellia sinensis leaves and buds, and pinpointing a pair of genes that participate in the antioxidant defense system of tea plants under sustained cold stress.

CONCLUSIONS

To summarize, downsampling and reorganizing the sample set improved the accuracy of co-expression analysis. Cluster-specific modules were more accurate in capturing condition-specific gene interactions. The introduction of CDV allowed for the assessment of condition specificity in gene co-expression analyses. Using this approach, we identified a series of high CDV transcription factor encoding genes related to sustained cold stress in Camellia sinensis. This study highlights the importance of considering condition specificity in co-expression analysis and provides insights into the regulation of the cold stress in Camellia sinensis.

Deciphering the regulatory networks involved in mild and severe salt stress responses in the roots of wild grapevine Vitis vinifera spp. sylvestris

Transcriptional regulatory networks are pivotal components of plant's response to salt stress. However, plant adaptation strategies varied as a function of stress intensity, which is mainly modulated by climate change. Here, we determined the gene regulatory networks based on transcription factor (TF) TF_gene co-expression, using two transcriptomic data sets generated from the salt-tolerant "Tebaba" roots either treated with 50 mM NaCl (mild stress) or 150 mM NaCl (severe stress). The analysis of these regulatory networks identified specific TFs as key regulatory hubs as evidenced by their multiple interactions with different target genes related to stress response. Indeed, under mild stress, NAC and bHLH TFs were identified as central hubs regulating nitrogen storage process. Moreover, HSF TFs were revealed as a regulatory hub regulating various aspects of cellular metabolism including flavonoid biosynthesis, protein processing, phenylpropanoid metabolism, galactose metabolism, and heat shock proteins. These processes are essentially linked to short-term acclimatization under mild salt stress. This was further consolidated by the protein-protein interaction (PPI) network analysis showing structural and plant growth adjustment. Conversely, under severe salt stress, dramatic metabolic changes were observed leading to novel TF members including MYB family as regulatory hubs controlling isoflavonoid biosynthesis, oxidative stress response, abscisic acid signaling pathway, and proteolysis. The PPI network analysis also revealed deeper stress defense changes aiming to restore plant metabolic homeostasis when facing severe salt stress. Overall, both the gene co-expression and PPI network provided valuable insights on key transcription factor hubs that can be employed as candidates for future genetic crop engineering programs.

Regulation of cardiac fibroblasts reprogramming into cardiomyocyte-like cells with a cocktail of small molecule compounds

Myocardial infarction results in extensive cardiomyocyte apoptosis, leading to the formation of noncontractile scar tissue. Given the limited regenerative capacity of adult mammalian cardiomyocytes, direct reprogramming of cardiac fibroblasts (CFs) into cardiomyocytes represents a promising therapeutic strategy for myocardial repair, and small molecule drugs might offer a more attractive alternative to gene editing approaches in terms of safety and clinical feasibility. This study aimed to reprogram rat CFs into cardiomyocytes using a small molecular chemical mixture comprising CHIR99021, Valproic acid, Dorsomorphin, SB431542, and Forskolin. Immunofluorescence analysis revealed a significant increase in the expression of cardiomyocyte-specific markers, including cardiac troponin T (cTnT), Connexin 43 (Cx43), α-actinin, and Tbx5. Changes in intracellular calcium ion levels and Ca2+ signal transfer between adjacent cells were monitored using a calcium ion fluorescence probe. mRNA sequencing analysis demonstrated the upregulation of genes associated with cardiac morphogenesis, myocardial differentiation, and muscle fiber contraction during CF differentiation induced by the small-molecule compounds. Conversely, the expression of fibroblast-related genes was downregulated. These findings suggest that chemical-induced cell fate conversion of rat CFs into cardiomyocyte-like cells is feasible, offering a potential therapeutic solution for myocardial injury.

Comparative analysis of differential gene expression reveals novel insights into the heteroblastic foliage functional traits of Pinus massoniana seedlings

Pinus massoniana needles, rich in medicinal polysaccharides and flavonoids, undergo heteroblastic foliage, transitioning from primary needles (PN) to secondary needles (SN) during growth, resulting in altered functional traits. Despite its significance, the molecular regulatory mechanisms governing these traits remain unclear. This study employs Iso-Seq and RNA-Seq analyses to explore differentially expressed genes (DEGs) associated with functional traits throughout the main growth season of heteroblastic foliage. Co-expression network analysis identified 34 hub genes and 17 key transcription factors (TFs) influencing light-harvesting antenna, photosystem I and II, crucial in photosynthesis regulation. Additionally, 14 genes involved in polysaccharide metabolism pathways, synthesizing sucrose, glucose, UDP sugars, and xylan, along with four genes in flavonoid biosynthesis pathways, regulating p-coumaroyl-CoA, quercetin, galangin, and myricetin production, exhibited differential expression between PN and SN. Further analysis unveils a highly interconnected network among these genes, forming a pivotal cascade of TFs and DEGs. Therefore, heteroblastic changes significantly impact needle functional traits, potentially affecting the pharmacological properties of PN and SN. Thus, these genomic insights into understanding the molecular-level differences of heteroblastic foliage, thereby establishing a foundation for advancements in the pharmaceutical industry related to needle-derived products.

Tea-Derived Polyphenols Enhance Drought Resistance of Tea Plants (Camellia sinensis) by Alleviating Jasmonate-Isoleucine Pathway and Flavonoid Metabolism Flow

Extreme drought weather has occurred frequently in recent years, resulting in serious yield loss in tea plantations. The study of drought in tea plantations is becoming more and more intensive, but there are fewer studies on drought-resistant measures applied in actual production. Therefore, in this study, we investigated the effect of exogenous tea polyphenols on the drought resistance of tea plant by pouring 100 mg·L-1 of exogenous tea polyphenols into the root under drought. The exogenous tea polyphenols were able to promote the closure of stomata and reduce water loss from leaves under drought stress. Drought-induced malondialdehyde (MDA) accumulation in tea leaves and roots was also significantly reduced by exogenous tea polyphenols. Combined transcriptomic and metabolomic analyses showed that exogenous tea polyphenols regulated the abnormal responses of photosynthetic and energy metabolism in leaves under drought conditions and alleviated sphingolipid metabolism, arginine metabolism, and glutathione metabolism in the root system, which enhanced the drought resistance of tea seedlings. Exogenous tea polyphenols induced jasmonic acid-isoleucine (JA-ILE) accumulation in the root system, and the jasmonic acid-isoleucine synthetase gene (TEA028623), jasmonic acid ZIM structural domain proteins (JAMs) synthesis genes (novel.22237, TEA001821), and the transcription factor MYC2 (TEA014288, TEA005840) were significantly up-regulated. Meanwhile, the flavonoid metabolic flow was significantly altered in the root; for example, the content of EGCG, ECG, and EGC was significantly increased. Thus, exogenous tea polyphenols enhance the drought resistance of tea plants through multiple pathways.

The Laccase Family Gene CsLAC37 Participates in Resistance to Colletotrichum gloeosporioides Infection in Tea Plants

Fungal attacks have become a major obstacle in tea plantations. Colletotrichum gloeosporioides is one of the most devastating fungal pathogens in tea plantations that can severely affect tea yield and quality. However, the molecular mechanism of resistance genes involved in anthracnose is still largely unknown in tea plants. Here, we found that the laccase gene CsLAC37 was involved in the response to fungal infection based on a transcriptome analysis. The full-length CDS of CsLAC37 was cloned, and its protein sequence had the closest relationship with the Arabidopsis AtLAC15 protein compared to other AtLACs. Tissue-specific expression analysis showed that CsLAC37 had higher expression levels in mature leaves and stems than in the other tissues. Subcellular localization showed that the CsLAC37 protein was predominantly localized in the cell membrane. The expression levels of CsLAC37 were upregulated at different time points under cold, salt, SA, and ABA treatments. qRT-PCR confirmed that CsLAC37 responded to both Pestalotiopsis-like species and C. gloeosporioides infections. Functional validation showed that the hydrogen peroxide (H2O2) content increased significantly, and POD activity decreased in leaves after antisense oligonucleotide (AsODN) treatment compared to the controls. The results demonstrated that CsLAC37 may play an important role in resistance to anthracnose, and the findings provide a theoretical foundation for molecular breeding of tea varieties with resistance to fungal diseases.

Identification and expression analysis of the KNOX genes during organogenesis and stress responseness in Camellia sinensis (L.) O. Kuntze

Tea plant (Camellia sinensis L.), whose leaves are the major reproductive organs, has been cultivated and consumed widely for its economic and health benefits. The Knotted1-like Homeobox (KNOX) proteins play significant roles in leaf morphology formation and development. However, the functions of KNOX proteins in tea plants are still unknown. Here, 11 CsKNOX genes from the tea plants were cloned and divided into Class I, II, and KNATM clades based on their protein sequences. These 11 CsKNOX genes were mapped on 8 out of 15 tea plant chromosomes, all localized in the nucleus. Specific spatiotemporal expression patterns of CsKNOX genes were found in various tissues and different development periods of buds, flowers, and roots of tea plants. Meanwhile, transcript levels of CsKNOX in tea leaves were strongly correlated with the accumulation of flavan-3-ols and proanthocyanidins. It was found that most of the CsKNOX genes could respond to drought, salt, cold, and exogenous MeJA and GA3 by analysis of transcriptomics data and promoter elements. The protein interaction analysis showed that CsKNOX could cooperate with CsAS1 and other critical functional proteins. In conclusion, this research provided the basic information for the functions of the CsKNOX family during organogenesis and stress response in tea plants, which was necessary for further functional characterization verification.

The R2R3-MYB transcription factor CsMYB42 regulates theanine biosynthesis in albino tea leaves

Theanine is a unique secondary metabolite in tea plants and contributes to the umami taste and health benefits of tea. However, theanine biosynthesis in tea plants is not fully understood, and its mechanism of transcriptional regulation remains poorly reported. Theanine content was significantly correlated with the expression of theanine biosynthesis-related gene CsGS1c and transcription factor CsMYB42 in different leaf positions and picking times, but there was no significant correlation in different tissues of albino tea plant 'Anjibaicha'. This suggests that CsMYB42 may regulate CsGS1c to synthesize theanine in albino tea leaves, and the regulation is tissue specific. CsMYB42 is a nuclear-localized R2R3-MYB transcription factor gene with transcriptional activation activity. Yeast one-hybrid assay and electrophoretic mobility shift assay confirmed the direct binding of CsMYB42 to the promoter of CsGS1c. Luciferase assay showed that CsMYB42 activates the CsGS1c expression. Furthermore, the inhibition of CsMYB42 using an antisense oligonucleotide in tea leaves decreased CsGS1c expression and theanine content. These results indicate that CsMYB42 plays a crucial role in activating the expression of CsGS1c and may be involved in the biosynthesis of theanine in albino tea leaves. This study provides fresh insights into the tissue-specific regulation of theanine biosynthesis, which laid a foundation for breeding high-theanine tea plants.

Comparative transcriptome analysis of sensory genes from the antenna and abdomen of Quadrastichus mendeli Kim

Quadrastichus mendeli Kim is one of the most important parasitoids of Leptocybe invasa Fisher et La Salle, which is an invasive gall-making pest in eucalyptus plantations in the world. Gall-inducing insects live within plant tissues and induce tumor-like growths that provide the insects with food, shelter, and protection from natural enemies. Empirical evidences showed that sensory genes play a key role in the host location of parasitoids. So far, what kind of sensory genes regulate parasitoids to locate gall-inducing insects has not been uncovered. In this study, sensory genes in the antenna and abdomen of Q. mendeli were studied using high-throughput sequencing. In total, 181,543 contigs was obtained from the antenna and abdomen transcriptome of Q. mendeli. The major sensory genes (chemosensory proteins, CSPs; gustatory receptors, GRs; ionotropic receptors, IRs; odorant binding proteins, OBPs; odorant receptors, ORs; and sensory neuron membrane proteins, SNMPs) were identified, and phylogenetic analyses were performed with these genes from Q. mendeli and other model insect species. The gene co-expression network constructed by WGCNA method is robust and reliable. There were 10,314 differentially expressed genes (DEGs), and among them, 99 genes were DEGs. A comprehensive sequence resource with desirable quality was built by comparative transcriptome of the antenna and abdomen of Q. mendeli, enriching the genomic platform of Q. mendeli.

Low light stress promotes new tiller regeneration by changing source-sink relationship and activating expression of expansin genes in wheat

Low light stress seriously decreased wheat grain number through the formation of aborted spike during the reproductive period and induced new tiller regeneration to offset the loss of grain number. However, the mechanism by which plants coordinate spike aborted growth and the regeneration of new tillers remains unknown. To better understand this coordinated process, morphological, physiological and transcriptomic analyses were performed under low light stress at the young microspore stage. Our findings indicated that leaves exhausted most stored carbohydrates in 1 day of darkness. However, spike and uppermost internode (UI) were converted from sink to source, due to increased abscisic acid (ABA) content and decreased cytokinin content. During this process, genes encoding amylases, Sugars Will Eventually be Exported Transporters (SWEET) and sucrose transporters or sucrose carriers (SUT/SUC) were upregulated in spike and UI, which degraded starch into soluble sugars and loaded them into the phloem. Subsequently, soluble sugars were transported to tiller node (TN) where cytokinin and auxin content increased and ABA content decreased, followed by unloading into TN cells by upregulated cell wall invertase (CWINV) genes and highly expressed H⁺ /hexose symporter genes. Finally, expansin genes integrated the sugar pathway and hormone pathway, and regulate the formation of new tillers directly.

Tandemly duplicated CYP82Ds catalyze 14-hydroxylation in triptolide biosynthesis and precursor production in Saccharomyces cerevisiae

Triptolide is a valuable multipotent antitumor diterpenoid in Tripterygium wilfordii, and its C-14 hydroxyl group is often selected for modification to enhance both the bioavailability and antitumor efficacy. However, the mechanism for 14-hydroxylation formation remains unknown. Here, we discover 133 kb of tandem duplicated CYP82Ds encoding 11 genes on chromosome 12 and characterize CYP82D274 and CYP82D263 as 14-hydroxylases that catalyze the metabolic grid in triptolide biosynthesis. The two CYP82Ds catalyze the aromatization of miltiradiene, which has been repeatedly reported to be a spontaneous process. In vivo assays and evaluations of the kinetic parameters of CYP82Ds indicate the most significant affinity to dehydroabietic acid among multiple intermediates. The precursor 14-hydroxy-dehydroabietic acid is successfully produced by engineered Saccharomyces cerevisiae. Our study provides genetic elements for further elucidation of the downstream biosynthetic pathways and heterologous production of triptolide and of the currently intractable biosynthesis of other 14-hydroxyl labdane-type secondary metabolites.

Cytochrome P450s in plant terpenoid biosynthesis: discovery, characterization and metabolic engineering

As the largest family of natural products, terpenoids play valuable roles in medicine, agriculture, cosmetics and food. However, the traditional methods that rely on direct extraction from the original plants not only produce low yields, but also result in waste of resources, and are not applicable at all to endangered species. Modern heterologous biosynthesis is considered a promising, efficient, and sustainable production method, but it relies on the premise of a complete analysis of the biosynthetic pathway of terpenoids, especially the functionalization processes involving downstream cytochrome P450s. In this review, we systematically introduce the biotech approaches used to discover and characterize plant terpenoid-related P450s in recent years. In addition, we propose corresponding metabolic engineering approaches to increase the effective expression of P450 and improve the yield of terpenoids, and also elaborate on metabolic engineering strategies and examples of heterologous biosynthesis of terpenoids in Saccharomyces cerevisiae and plant hosts. Finally, we provide perspectives for the biotech approaches to be developed for future research on terpenoid-related P450.

Weighted gene co-expression network analysis identifies genes related to HG Type 0 resistance and verification of hub gene GmHg1

INTRODUCTION

The soybean cyst nematode (SCN) is a major disease in soybean production thatseriously affects soybean yield. At present, there are no studies on weighted geneco-expression network analysis (WGCNA) related to SCN resistance.

METHODS

Here, transcriptome data from 36 soybean roots under SCN HG Type 0 (race 3) stresswere used in WGCNA to identify significant modules.

RESULTS AND DISCUSSION

A total of 10,000 differentially expressed genes and 21 modules were identified, of which the module most related to SCN was turquoise. In addition, the hub gene GmHg1 with high connectivity was selected, and its function was verified. GmHg1 encodes serine/threonine protein kinase (PK), and the expression of GmHg1 in SCN-resistant cultivars ('Dongnong L-204') and SCN-susceptible cultivars ('Heinong 37') increased significantly after HG Type 0 stress. Soybean plants transformed with GmHg1-OX had significantly increased SCN resistance. In contrast, the GmHg1-RNAi transgenic soybean plants significantly reduced SCN resistance. In transgenic materials, the expression patterns of 11 genes with the same expression trend as the GmHg1 gene in the 'turquoise module' were analyzed. Analysis showed that 11genes were co-expressed with GmHg1, which may be involved in the process of soybean resistance to SCN. Our work provides a new direction for studying the Molecular mechanism of soybean resistance to SCN.

A Practical Guide to Inferring Multi-Omics Networks in Plant Systems

The inference of gene regulatory networks can reveal molecular connections underlying biological processes and improve our understanding of complex biological phenomena in plants. Many previous network studies have inferred networks using only one type of omics data, such as transcriptomics. However, given more recent work applying multi-omics integration in plant biology, such as combining (phospho)proteomics with transcriptomics, it may be advantageous to integrate multiple omics data types into a comprehensive network prediction. Here, we describe a state-of-the-art approach for integrating multi-omics data with gene regulatory network inference to describe signaling pathways and uncover novel regulators. We detail how to download and process transcriptomics and (phospho)proteomics data for network inference, using an example dataset from the plant hormone signaling field. We provide a step-by-step protocol for inference, visualization, and analysis of an integrative multi-omics network using currently available methods. This chapter serves as an accessible guide for novice and intermediate bioinformaticians to analyze their own datasets and reanalyze published work.

Integrated multi-omics reveals the molecular mechanisms underlying efficient phosphorus use under phosphate deficiency in elephant grass (Pennisetum purpureum)

Phosphorus (P) is an essential macronutrient element for plant growth, and deficiency of inorganic phosphate (Pi) limits plant growth and yield. Elephant grass (Pennisetum purpureum) is an important fodder crop cultivated widely in tropical and subtropical areas throughout the world. However, the mechanisms underlying efficient P use in elephant grass under Pi deficiency remain poorly understood. In this study, the physiological and molecular responses of elephant grass leaves and roots to Pi deficiency were investigated. The results showed that dry weight, total P concentration, and P content decreased in Pi-deprived plants, but that acid phosphatase activity and P utilization efficiency (PUE) were higher than in Pi-sufficient plants. Regarding Pi starvation-responsive (PSR) genes, transcriptomics showed that 59 unigenes involved in Pi acquisition and transport (especially 18 purple acid phosphatase and 27 phosphate transporter 1 unigenes) and 51 phospholipase unigenes involved in phospholipids degradation or Pi-free lipids biosynthesis, as well as 47 core unigenes involved in the synthesis of phenylpropanoids and flavonoids, were significantly up-regulated by Pi deprivation in leaves or roots. Furthermore, 43 unigenes related to Pi-independent- or inorganic pyrophosphate (PPi)-dependent bypass reactions were markedly up-regulated in Pi-deficient leaves, especially five UDP-glucose pyrophosphorylase and 15 phosphoenolpyruvate carboxylase unigenes. Consistent with PSR unigene expression changes, metabolomics revealed that Pi deficiency significantly increased metabolites of Pi-free lipids, phenylpropanoids, and flavonoids in leaves and roots, but decreased phospholipid metabolites. This study reveals the mechanisms underlying the responses to Pi starvation in elephant grass leaves and roots, which provides candidate unigenes involved in efficient P use and theoretical references for the development of P-efficient elephant grass varieties.

The Current Developments in Medicinal Plant Genomics Enabled the Diversification of Secondary Metabolites' Biosynthesis

Medicinal plants produce important substrates for their adaptation and defenses against environmental factors and, at the same time, are used for traditional medicine and industrial additives. Plants have relatively little in the way of secondary metabolites via biosynthesis. Recently, the whole-genome sequencing of medicinal plants and the identification of secondary metabolite production were revolutionized by the rapid development and cheap cost of sequencing technology. Advances in functional genomics, such as transcriptomics, proteomics, and metabolomics, pave the way for discoveries in secondary metabolites and related key genes. The multi-omics approaches can offer tremendous insight into the variety, distribution, and development of biosynthetic gene clusters (BGCs). Although many reviews have reported on the plant and medicinal plant genome, chemistry, and pharmacology, there is no review giving a comprehensive report about the medicinal plant genome and multi-omics approaches to study the biosynthesis pathway of secondary metabolites. Here, we introduce the medicinal plant genome and the application of multi-omics tools for identifying genes related to the biosynthesis pathway of secondary metabolites. Moreover, we explore comparative genomics and polyploidy for gene family analysis in medicinal plants. This study promotes medicinal plant genomics, which contributes to the biosynthesis and screening of plant substrates and plant-based drugs and prompts the research efficiency of traditional medicine.

FunOrder 2.0 - a method for the fully automated curation of co-evolved genes in fungal biosynthetic gene clusters

Coevolution is an important biological process that shapes interacting proteins - may it be physically interacting proteins or consecutive enzymes in a metabolic pathway, such as the biosynthetic pathways for secondary metabolites. Previously, we developed FunOrder, a semi-automated method for the detection of co-evolved genes, and demonstrated that FunOrder can be used to identify essential genes in biosynthetic gene clusters from different ascomycetes. A major drawback of this original method was the need for a manual assessment, which may create a user bias and prevents a high-throughput application. Here we present a fully automated version of this method termed FunOrder 2.0. In the improved version, we use several mathematical indices to determine the optimal number of clusters in the FunOrder output, and a subsequent k-means clustering based on the first three principal components of a principal component analysis of the FunOrder output to automatically detect co-evolved genes. Further, we replaced the BLAST tool with the DIAMOND tool as a prerequisite for using larger proteome databases. Potentially, FunOrder 2.0 may be used for the assessment of complete genomes, which has not been attempted yet. However, the introduced changes slightly decreased the sensitivity of this method, which is outweighed by enhanced overall speed and specificity.

Weighted Gene Correlation Network Analysis (WGCNA) of Arabidopsis Somatic Embryogenesis (SE) and Identification of Key Gene Modules to Uncover SE-Associated Hub Genes

Somatic embryogenesis (SE), which occurs naturally in many plant species, serves as a model to elucidate cellular and molecular mechanisms of embryo patterning in plants. Decoding the regulatory landscape of SE is essential for its further application. Hence, the present study was aimed at employing Weighted Gene Correlation Network Analysis (WGCNA) to construct a gene coexpression network (GCN) for Arabidopsis SE and then identifying highly correlated gene modules to uncover the hub genes associated with SE that may serve as potential molecular targets. A total of 17,059 genes were filtered from a microarray dataset comprising four stages of SE, i.e., stage I (zygotic embryos), stage II (proliferating tissues at 7 days of induction), stage III (proliferating tissues at 14 days of induction), and stage IV (mature somatic embryos). This included 1,711 transcription factors and 445 EMBRYO DEFECTIVE genes. GCN analysis identified a total of 26 gene modules with the module size ranging from 35 to 3,418 genes using a dynamic cut tree algorithm. The module-trait analysis revealed that four, four, seven, and four modules were associated with stages I, II, III, and IV, respectively. Further, we identified a total of 260 hub genes based on the degree of intramodular connectivity. Validation of the hub genes using publicly available expression datasets demonstrated that at least 78 hub genes are potentially associated with embryogenesis; of these, many genes remain functionally uncharacterized thus far. In silico promoter analysis of these genes revealed the presence of cis-acting regulatory elements, “soybean embryo factor 4 (SEF4) binding site,” and “E-box” of the napA storage-protein gene of Brassica napus; this suggests that these genes may play important roles in plant embryo development. The present study successfully applied WGCNA to construct a GCN for SE in Arabidopsis and identified hub genes involved in the development of somatic embryos. These hub genes could be used as molecular targets to further elucidate the molecular mechanisms underlying SE in plants.

Comparative transcriptomic analysis unveils the deep phylogeny and secondary metabolite evolution of 116 Camellia plants

Camellia plants include more than 200 species of great diversity and immense economic, ornamental, and cultural values. We sequenced the transcriptomes of 116 Camellia plants from almost all sections of the genus Camellia. We constructed a pan-transcriptome of Camellia plants with 89 394 gene families and then resolved the phylogeny of genus Camellia based on 405 high-quality low-copy core genes. Most of the inferred relationships are well supported by multiple nuclear gene trees and morphological traits. We provide strong evidence that Camellia plants shared a recent whole genome duplication event, followed by large expansions of transcription factor families associated with stress resistance and secondary metabolism. Secondary metabolites, particularly those associated with tea quality such as catechins and caffeine, were preferentially heavily accumulated in the Camellia plants from section Thea. We thoroughly examined the expression patterns of hundreds of genes associated with tea quality, and found that some of them exhibited significantly high expression and correlations with secondary metabolite accumulations in Thea species. We also released a web-accessible database for efficient retrieval of Camellia transcriptomes. The reported transcriptome sequences and obtained novel findings will facilitate the efficient conservation and utilization of Camellia germplasm towards a breeding program for cultivated tea, camellia, and oil-tea plants.

Effects of growth years on ginsenoside biosynthesis of wild ginseng and cultivated ginseng

Background

Ginsenoside, as the main active substance in ginseng, has the function of treating various diseases. However, the ginsenosides content of cultivated ginseng is obviously affected by the growth years, but the molecular mechanism is not clear. In addition, there are significant differences in morphology and physiology between wild ginseng and cultivated ginseng, and the effect of growth years on ginsenoside synthesis not yet understood in wild ginseng.

Results

Transcriptome sequencing on the roots, stems and leaves of cultivated ginseng and wild ginseng with different growth years was performed in this study, exploring the effect of growth years on gene expression in ginseng. The number of differentially expressed genes (DEGs) from comparison groups in cultivated ginseng was higher than that in wild ginseng. The result of weighted gene co-expression network analysis (WGCNA) showed that growth years significantly affected the gene expression of Mitogen-activated protein kinases (MAPK) signaling pathway and terpenoid backbone biosynthesis pathway in cultivated ginseng, but had no effects in wild ginseng. Furthermore, the growth years had significant effects on the genes related to ginsenoside synthesis in cultivated ginseng, and the effects were different in the roots, stems and leaves. However, it had little influence on the expression of genes related to ginsenoside synthesis in wild ginseng. Growth years might affect the expression of genes for ginsenoside synthesis by influencing the expression of these transcription factors (TFs), like my elob lastosis (MYB), NAM, ATAF1 and 2, and CUC2 (NAC), APETALA2/ethylene-responsive factor (AP2/ERF), basic helix-loop-helix (bHLH) and WRKY, etc., thereby affecting the content of ginsenosides.

Conclusions

This study complemented the gaps in the genetic information of wild ginseng in different growth periods and helped to clarify the potential mechanisms of the effect of growth years on the physiological state in wild ginseng and cultivated ginseng, which also provided a new insight into the mechanism of ginsenoside regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08570-0.

Distance correlation application to gene co-expression network analysis

BACKGROUND

To construct gene co-expression networks, it is necessary to evaluate the correlation between different gene expression profiles. However, commonly used correlation metrics, including both linear (such as Pearson's correlation) and monotonic (such as Spearman's correlation) dependence metrics, are not enough to observe the nature of real biological systems. Hence, introducing a more informative correlation metric when constructing gene co-expression networks is still an interesting topic.

RESULTS

In this paper, we test distance correlation, a correlation metric integrating both linear and non-linear dependence, with other three typical metrics (Pearson's correlation, Spearman's correlation, and maximal information coefficient) on four different arrays (macrophage and liver) and RNA-seq (cervical cancer and pancreatic cancer) datasets. Among all the metrics, distance correlation is distribution free and can provide better performance on complex relationships and anti-outlier. Furthermore, distance correlation is applied to Weighted Gene Co-expression Network Analysis (WGCNA) for constructing a gene co-expression network analysis method which we named Distance Correlation-based Weighted Gene Co-expression Network Analysis (DC-WGCNA). Compared with traditional WGCNA, DC-WGCNA can enhance the result of enrichment analysis and improve the module stability.

CONCLUSIONS

Distance correlation is better at revealing complex biological relationships between gene profiles compared with other correlation metrics, which contribute to more meaningful modules when analyzing gene co-expression networks. However, due to the high time complexity of distance correlation, the implementation requires more computer memory.

Gene Co-expression Network Analysis

Gene co-expression analysis is a data analysis technique that helps identify groups of genes with similar expression patterns across several different conditions. By means of these techniques, different groups have been able to assign putative metabolic pathways and functions to understudied genes and to identify novel metabolic regulation networks for different metabolites. Some groups have even used network comparative studies to understand the evolution of these networks from green algae to land plants. In this chapter, we will go over the basic definitions required to understand network topology and gene module identification. Additionally, we offer the reader a walk-through a standard analysis pipeline as implemented in the package WGCNA that takes as input raw fastq files and obtains co-expressed gene clusters and representative gene expression patterns from each module for downstream applications.

From gene to biomolecular networks: a review of evidences for understanding complex biological function in plants

Although at the infancy stage, biomolecular network biology is a comprehensive approach to understand complex biological function in plants. Recent advancements in the accumulation of multi-omics data coupled with computational approach have accelerated our current understanding of the complexities of gene function at the system level. Biomolecular networks such as protein-protein interaction, co-expression and gene regulatory networks have extensively been used to decipher the intricacies of transcriptional reprogramming of hundreds of genes and their regulatory interaction in response to various environmental perturbations mainly in the model plant Arabidopsis. This review describes recent applications of network-based approaches to understand the biological functions in plants and focuses on the challenges and opportunities to harness the full potential of the approach.

Light control of catechin accumulation is mediated by photosynthetic capacity in tea plant (Camellia sinensis)

BACKGROUND

Catechins are crucial in determining the flavour and health benefits of tea, but it remains unclear that how the light intensity regulates catechins biosynthesis. Therefore, we cultivated tea plants in a phytotron to elucidate the response mechanism of catechins biosynthesis to light intensity changes.

RESULTS

In the 250 μmol·m^- 2·s^- 1 treatment, the contents of epigallocatechin, epigallocatechin gallate and total catechins were increased by 98.94, 14.5 and 13.0% respectively, compared with those in the 550 μmol·m^- 2·s^- 1 treatment. Meanwhile, the photosynthetic capacity was enhanced in the 250 μmol·m^- 2·s^- 1 treatment, including the electron transport rate, net photosynthetic rate, transpiration rate and expression of related genes (such as CspsbA, CspsbB, CspsbC, CspsbD, CsPsbR and CsGLK1). In contrast, the extremely low or high light intensity decreased the catechins accumulation and photosynthetic capacity of the tea plants. The comprehensive analysis revealed that the response of catechins biosynthesis to the light intensity was mediated by the photosynthetic capacity of the tea plants. Appropriately high light upregulated the expression of genes related to photosynthetic capacity to improve the net photosynthetic rate (Pn), transpiration rate (Tr), and electron transfer rate (ETR), which enhanced the contents of substrates for non-esterified catechins biosynthesis (such as EGC). Meanwhile, these photosynthetic capacity-related genes and gallic acid (GA) biosynthesis-related genes (CsaroB, CsaroDE1, CsaroDE2 and CsaroDE3) co-regulated the response of GA accumulation to light intensity. Eventually, the epigallocatechin gallate content was enhanced by the increased contents of its precursors (EGC and GA) and the upregulation of the CsSCPL gene.

CONCLUSIONS

In this study, the catechin content and photosynthetic capacity of tea plants increased under appropriately high light intensities (250 μmol·m^- 2·s^- 1 and 350 μmol·m^- 2·s^- 1) but decreased under extremely low or high light intensities (150 μmol·m^- 2·s^- 1 or 550 μmol·m^- 2·s^- 1). We found that the control of catechin accumulation by light intensity in tea plants is mediated by the plant photosynthetic capacity. The research provided useful information for improving catechins content and its light-intensity regulation mechanism in tea plant.

Research progress on the response of tea catechins to drought stress

Drought stress (DS) is the most important abiotic stress affecting yield and quality of tea worldwide. DS causes oxidative stress to cells due to the accumulation of reactive oxygen species (ROS). As non-enzymatic antioxidants, tea catechins can scavenge excess ROS in response to DS. Further, catechin accumulation contributes to the formation of oxidative polymerization products (e.g. theaflavins and thearubigins) that improve the quality of black tea. However, there are no systematic reports on the response of tea catechins to DS. First, we reviewed the available literature on the response of tea plants to DS. Second, we summarized the current knowledge of ROS production in tea leaves under DS and typical antioxidant response mechanisms. Third, we conducted a detailed review of the changes in catechin levels in tea under different drought conditions. We found that the total amounts of catechin and o-quinone increased under DS conditions. We propose that the possible mechanisms underlying tea catechin accumulation under DS conditions include (i) autotrophic formation of o-quinone, (ii) polymerization of proanthocyanidins that directly scavenge excess ROS, and (iii) formation of metal ion complexes and by influencing the antioxidant systems that indirectly eliminate excess ROS. Finally, we discuss ways of potentially improving black tea quality using drought before picking in the summer/fall dry season. In summary, we mainly discuss the antioxidant mechanisms of tea catechins under DS and the possibility of using drought to improve black tea quality. Our review provides a theoretical basis for the production of high-quality black tea under DS conditions. © 2021 Society of Chemical Industry.

Weighted gene co-expression network analysis unveils gene networks regulating folate biosynthesis in maize endosperm

Folates are essential elements for human growth and development, and their deficiency can lead to serious disorders. Waxy maize is a rich source of folates; however, the regulatory mechanism underlying folate biosynthesis in the endosperm remains unclear. Here, we examined changes in the folate content of maize endosperm collected at 15, 18, 21, 24, and 27 days after pollination (DAP) using liquid chromatograph-mass spectrometry and identified genes related to folate biosynthesis using transcriptome sequencing data. The results showed that 5-methyl-tetrahydrofolate and 5,10-methylene tetrahydrofolate were the main storage forms of folates in the endosperm, and their contents were relatively high at 21-24 days. We also identified 569, 3183, 4365, and 5513 differentially expressed genes (DEGs) in different days around milk stage. Functional annotation revealed 518 transcription factors (TFs) belonging to 33 families exhibiting specific expression in at least one sampling time. The key hub genes involved in folate biosynthesis were identified by weighted gene co-expression network analysis. In total, 24,976 genes were used to construct a co-expression network with 29 co-expression modules, among which the brown and purple modules were highly related to folate biosynthesis. Further, 187 transcription factors in the brown and purple modules were considered potential transcription factors related to endosperm folate biosynthesis. These results may improve the understanding of the molecular mechanism underlying folate biosynthesis in waxy maize and lead to the development of nutritionally fortified varieties.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s13205-021-02974-7.

FunOrder: A robust and semi-automated method for the identification of essential biosynthetic genes through computational molecular co-evolution

Secondary metabolites (SMs) are a vast group of compounds with different structures and properties that have been utilized as drugs, food additives, dyes, and as monomers for novel plastics. In many cases, the biosynthesis of SMs is catalysed by enzymes whose corresponding genes are co-localized in the genome in biosynthetic gene clusters (BGCs). Notably, BGCs may contain so-called gap genes, that are not involved in the biosynthesis of the SM. Current genome mining tools can identify BGCs, but they have problems with distinguishing essential genes from gap genes. This can and must be done by expensive, laborious, and time-consuming comparative genomic approaches or transcriptome analyses. In this study, we developed a method that allows semi-automated identification of essential genes in a BGC based on co-evolution analysis. To this end, the protein sequences of a BGC are blasted against a suitable proteome database. For each protein, a phylogenetic tree is created. The trees are compared by treeKO to detect co-evolution. The results of this comparison are visualized in different output formats, which are compared visually. Our results suggest that co-evolution is commonly occurring within BGCs, albeit not all, and that especially those genes that encode for enzymes of the biosynthetic pathway are co-evolutionary linked and can be identified with FunOrder. In light of the growing number of genomic data available, this will contribute to the studies of BGCs in native hosts and facilitate heterologous expression in other organisms with the aim of the discovery of novel SMs.

The discovery and description of novel fungal secondary metabolites promises novel antibiotics, pharmaceuticals, and other useful compounds. A way to identify novel secondary metabolites is to express the corresponding genes in a suitable expression host. Consequently, a detailed knowledge or an accurate prediction of these genes is necessary. In fungi, the genes are co-localized in so-called biosynthetic gene clusters. Notably, the clusters may also contain genes that are not necessary for the biosynthesis of the secondary metabolites, so-called gap genes. We developed a method to detect co-evolved genes within the clusters and demonstrated that essential genes are co-evolving and can thus be differentiated from the gap genes. This adds an additional layer of information, which can support researchers with their decisions on which genes to study and express for the discovery of novel secondary metabolites.

Identification of Key Transcription Factors Related to Bacterial Spot Resistance in Pepper through Regulatory Network Analyses

Bacterial spot (BS), caused by Xanthomonas campestris pv. Vesicatoria (Xcv), severely affects the quality and yield of pepper. Thus, breeding new pepper cultivars with enhanced resistance to BS can improve economic benefits for pepper production. Identification of BS resistance genes is an essential step to achieve this goal. However, very few BS resistance genes have been well characterized in pepper so far. In this study, we reanalyzed public multiple time points related to RNA-seq data sets from two pepper cultivars, the Xcv-susceptible cultivar ECW and the Xcv-resistant cultivar VI037601, post Xcv infection. We identified a total of 3568 differentially expressed genes (DEGs) between two cultivars post Xcv infection, which were mainly involved in some biological processes, such as Gene Ontology (GO) terms related to defense response to bacterium, immune system process, and regulation of defense response, etc. Through weighted gene co-expression network analysis (WGCNA), we identified 15 hub (Hub) transcription factor (TF) candidates in response to Xcv infection. We further selected 20 TFs from the gene regulatory network (GRN) potentially involved in Xcv resistance response. Finally, we predicted 4 TFs, C3H (p-coumarate 3-hydroxylase), ERF (ethylene-responsive element binding factor), TALE (three-amino-acid-loop-extension), and HSF (heat shock transcription factor), as key factors responsible for BS disease resistance in pepper. In conclusion, our study provides valuable resources for dissecting the underlying molecular mechanism responsible for Xcv resistance in pepper. Additionally, it also provides valuable references for mining transcriptomic data to identify key candidates for disease resistance in horticulture crops.

Gene-Metabolite Network Analysis Revealed Tissue-Specific Accumulation of Therapeutic Metabolites in Mallotus japonicus

Mallotus japonicus is a valuable traditional medicinal plant in East Asia for applications as a gastrointestinal drug. However, the molecular components involved in the biosynthesis of bioactive metabolites have not yet been explored, primarily due to a lack of omics resources. In this study, we established metabolome and transcriptome resources for M. japonicus to capture the diverse metabolite constituents and active transcripts involved in its biosynthesis and regulation. A combination of untargeted metabolite profiling with data-dependent metabolite fragmentation and metabolite annotation through manual curation and feature-based molecular networking established an overall metabospace of M. japonicus represented by 2129 metabolite features. M. japonicus de novo transcriptome assembly showed 96.9% transcriptome completeness, representing 226,250 active transcripts across seven tissues. We identified specialized metabolites biosynthesis in a tissue-specific manner, with a strong correlation between transcripts expression and metabolite accumulations in M. japonicus. The correlation- and network-based integration of metabolome and transcriptome datasets identified candidate genes involved in the biosynthesis of key specialized metabolites of M. japonicus. We further used phylogenetic analysis to identify 13 C-glycosyltransferases and 11 methyltransferases coding candidate genes involved in the biosynthesis of medicinally important bergenin. This study provides comprehensive, high-quality multi-omics resources to further investigate biological properties of specialized metabolites biosynthesis in M. japonicus.

Transcriptome profiling reveal key hub genes in co-expression networks involved in Iridoid glycosides biosynthetic machinery in Picrorhiza kurroa

Picrorhiza kurroa is a medicinal herb rich in hepatoprotective iridoid glycosides, picroside-I (P-I) and picroside-II (P-II). The biosynthetic machinery of picrosides is poorly understood, therefore, 'no-direction' gene co-expression networks were used to extract linked/closed and separated interactions in terpenoid glycosides-specific sub-networks. Transcriptomes generated from different organs, varying for P-I and P-II contents such as shoots grown at 15 and 25 °C and nursery-grown shoots, stolons, and roots resulted in 47,726, 44,958, 40,117, 66,979, and 55,578 annotated transcripts, respectively. Occurrence of 2810 ± 136 nodes and 15,626 ± 696 edges in these networks indicated intense, co-expressed, closed loop interactions. Either deregulation/inhibition of abscisic acid (ABA) biosynthesis/signaling or constitutive degradation of ABA resulted in organ-specific accumulation of P-I and P-II. Biosynthesis, condensation and glucosylation of isoprene units may occur in shoots, roots or stolons; but addition of phenylpropanoid moiety and further modification/s of the iridoid backbone occurs mainly inside vacuoles in roots.

Combined Analysis of Volatile Terpenoid Metabolism and Transcriptome Reveals Transcription Factors Related to Terpene Synthase in Two Cultivars of Dendrobium officinale Flowers

Dendrobium officinale is a kind of traditional Chinese herbal medicine. Its flowers could be used as health care tea for its aroma flavor and medicinal value. Most recent studies demonstrated that terpenoids are the main components of the aromatic compounds in the flowers, but the biosynthesis of terpenoids is poorly understood in D. officinale. In the experiment, the flowers from two cultivars of D. officinale with different smells were collected. The transcriptome analysis and combined volatile terpenoids determination were performed to identify the genes related to the biosynthesis of the terpenoids. The results showed that the different products of volatile terpenoids are α-thujene, linalool, α-terpineol, α-phellandrene, γ-muurolene, α-patchoulene, and δ-elemene in two cultivar flowers. The transcriptome analysis detected 25,484 genes in the flowers. And 18,650 differentially expressed genes were identified between the two cultivars. Of these genes, 253 genes were mapped to the terpenoid metabolism pathway. Among these genes, 13 terpene synthase (TPS) genes may have correlations with AP2/ERF, WRKY, MYB, bHLH, and bZIP transcription factors by weighted gene co-expression network analysis (WGCNA). The transcription factors have regulatory effects on TPS genes. These results may provide ideas for the terpenoid biosynthesis and regulatory network of D. officinale flowers.

Genome-wide association analysis identified molecular markers associated with important tea flavor-related metabolites

The characteristic secondary metabolites in tea (theanine, caffeine, and catechins) are important factors contributing to unique tea flavors. However, there has been relatively little research on molecular markers related to these metabolites. Thus, we conducted a genome-wide association analysis of the levels of these tea flavor-related metabolites in three seasons. The theanine, caffeine, and catechin levels in Population 1 comprising 191 tea plant germplasms were examined, which revealed that their heritability exceeded 0.5 in the analyzed seasons, with the following rank order (highest to lowest heritabilities): (+)-catechin > (-)-gallocatechin gallate > caffeine = (-)-epicatechin > (-)-epigallocatechin-3-gallate > theanine > (-)-epigallocatechin > (-)-epicatechin-3-gallate > catechin gallate > (+)-gallocatechin. The SNPs detected by amplified-fragment SNP and methylation sequencing divided Population 1 into three groups and seven subgroups. An association analysis yielded 307 SNP markers related to theanine, caffeine, and catechins that were common to all three seasons. Some of the markers were pleiotropic. The functional annotation of 180 key genes at the SNP loci revealed that FLS, UGT, MYB, and WD40 domain-containing proteins, as well as ATP-binding cassette transporters, may be important for catechin synthesis. KEGG and GO analyses indicated that these genes are associated with metabolic pathways and secondary metabolite biosynthesis. Moreover, in Population 2 (98 tea plant germplasm resources), 30 candidate SNPs were verified, including 17 SNPs that were significantly or extremely significantly associated with specific metabolite levels. These results will provide a foundation for future research on important flavor-related metabolites and may help accelerate the breeding of new tea varieties.

The Wild Sugarcane and Sorghum Kinomes: Insights Into Expansion, Diversification, and Expression Patterns

The protein kinase (PK) superfamily is one of the largest superfamilies in plants and the core regulator of cellular signaling. Despite this substantial importance, the kinomes of sugarcane and sorghum have not been profiled. Here, we identified and profiled the complete kinomes of the polyploid Saccharum spontaneum (Ssp) and Sorghum bicolor (Sbi), a close diploid relative. The Sbi kinome was composed of 1,210 PKs; for Ssp, we identified 2,919 PKs when disregarding duplications and allelic copies, and these were related to 1,345 representative gene models. The Ssp and Sbi PKs were grouped into 20 groups and 120 subfamilies and exhibited high compositional similarities and evolutionary divergences. By utilizing the collinearity between the species, this study offers insights into Sbi and Ssp speciation, PK differentiation and selection. We assessed the PK subfamily expression profiles via RNA-Seq and identified significant similarities between Sbi and Ssp. Moreover, coexpression networks allowed inference of a core structure of kinase interactions with specific key elements. This study provides the first categorization of the allelic specificity of a kinome and offers a wide reservoir of molecular and genetic information, thereby enhancing the understanding of Sbi and Ssp PK evolutionary history.

BEM: Mining Coregulation Patterns in Transcriptomics via Boolean Matrix Factorization

MOTIVATION

The matrix factorization is an important way to analyze coregulation patterns in transcriptomic data, which can reveal the tumor signal perturbation status and subtype classification. However, current matrix factorization methods do not provide clear bicluster structure. Furthermore, these algorithms are based on the assumption of linear combination, which may not be sufficient to capture the coregulation patterns.

RESULTS

We presented a new algorithm for Boolean matrix factorization (BMF) via expectation maximization (BEM). BEM is more aligned with the molecular mechanism of transcriptomic coregulation and can scale to matrix with over 100 million data points. Synthetic experiments showed that BEM outperformed other BMF methods in terms of reconstruction error. Real-world application demonstrated that BEM is applicable to all kinds of transcriptomic data, including bulk RNA-seq, single-cell RNA-seq and spatial transcriptomic datasets. Given appropriate binarization, BEM was able to extract coregulation patterns consistent with disease subtypes, cell types or spatial anatomy.

AVAILABILITY AND IMPLEMENTATION

Python source code of BEM is available on https://github.com/LifanLiang/EM_BMF.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Integrated Transcriptome, microRNA, and Phytochemical Analyses Reveal Roles of Phytohormone Signal Transduction and ABC Transporters in Flavor Formation of Oolong Tea (Camellia sinensis) during Solar Withering

The unique aroma and flavor of oolong tea develop during the withering stage of postharvest processing. We explored the roles of miRNA-related regulatory networks during tea withering and their effects on oolong tea quality. We conducted transcriptome and miRNA analyses to identify differentially expressed (DE) miRNAs and target genes among fresh leaves, indoor-withered leaves, and solar-withered leaves. We identified 32 DE-miRNAs and 41 target genes involved in phytohormone signal transduction and ABC transporters. Further analyses indicated that these two pathways regulated the accumulation of flavor-related metabolites during tea withering. Flavonoid accumulation was correlated with the miR167d_1-ARF-GH3, miR845-ABCC1-3/ABCC2, miR166d-5p_1-ABCC1-2, and miR319c_3-PIF-ARF modules. Terpenoid content was correlated with the miR171b-3p_2-DELLA-MYC2 and miR166d-5p_1-ABCG2-MYC2 modules. These modules inhibited flavonoid biosynthesis and enhanced terpenoid biosynthesis in solar-withered leaves. Low auxin and gibberellic acid contents and circRNA-related regulatory networks also regulated the accumulation of flavor compounds in solar-withered leaves. Our analyses reveal how solar withering produces high-quality oolong tea.

Synthetic Biology of Cannabinoids and Cannabinoid Glucosides in Nicotiana benthamiana and Saccharomyces cerevisiae

Phytocannabinoids are a group of plant-derived metabolites that display a wide range of psychoactive as well as health-promoting effects. The production of pharmaceutically relevant cannabinoids relies on extraction and purification from cannabis (Cannabis sativa) plants yielding the major constituents, Δ9-tetrahydrocannabinol and cannabidiol. Heterologous biosynthesis of cannabinoids in Nicotiana benthamiana or Saccharomyces cerevisiae may provide cost-efficient and rapid future production platforms to acquire pure and high quantities of both the major and the rare cannabinoids as well as novel derivatives. Here, we used a meta-transcriptomic analysis of cannabis to identify genes for aromatic prenyltransferases of the UbiA superfamily and chalcone isomerase-like (CHIL) proteins. Among the aromatic prenyltransferases, CsaPT4 showed CBGAS activity in both N. benthamiana and S. cerevisiae. Coexpression of selected CsaPT pairs and of CHIL proteins encoding genes with CsaPT4 did not affect CBGAS catalytic efficiency. In a screen of different plant UDP-glycosyltransferases, Stevia rebaudiana SrUGT71E1 and Oryza sativa OsUGT5 were found to glucosylate olivetolic acid, cannabigerolic acid, and Δ9-tetrahydrocannabinolic acid. Metabolic engineering of N. benthamiana for production of cannabinoids revealed intrinsic glucosylation of olivetolic acid and cannabigerolic acid. S. cerevisiae was engineered to produce olivetolic acid glucoside and cannabigerolic acid glucoside.

Genome-wide identification of the HDAC family proteins and functional characterization of CsHD2C, a HD2-type histone deacetylase gene in tea plant (Camellia sinensis L. O. Kuntze)

The histone deacetylases (HDACs) are involved in growth, development and stress responses in many plants. However, the functions of HDACs in tea plant (Camellia sinensis L. O. Kuntze) and other woody plants remain unclear. Here, 18 CsHDAC genes were identified by genome-wide analysis in tea plant. The phylogenetic analysis demonstrated that the CsHDAC proteins were divided into three subfamilies, namely, the RPD3/HDA1 subfamily (8 members), the SIR2 subfamily (4 members) and the plant specific HD2 subfamily (6 members). The expression patterns showed that most members of CsHDACs family were regulated by different abiotic stress. High correlation was found between the expression of the CsHDACs and the accumulation of theanine, catechin, EGCG and other metabolites in tea plant. Most of the CsHDAC proteins were negative regulators. We further studied a specific gene CsHD2C (NCBI-ID: KY364373) in tea plant, which is the homolog of AtHD2C, encoded a protein of 306 aa. CsHD2C was highly expressed in leaves, young buds and stems. The transcription of CsHD2C was inhibited by ABA, NaCl and low temperature. It was found localized in the nucleus when fused with a YFP reporter gene. Overexpression of CsHD2C can rescue the phenotype related to different abiotic stresses in the mutant of AtHD2C in Arabidopsis. The stress-responsive genes RD29A, RD29B, ABI1 and ABI2 were also investigated to understand the regulating role of CsHD2C under abiotic stresses. We also found that CsHD2C could renew the change of acetylation level for histone H4 and the RNAP-II occupancy accumulation in the promoter of abiotic stress responses gene in the hd2c Arabidopsis mutant. Together, our results suggested that CsHD2C may act as a positive regulator in abiotic stress responses in tea plant.

Identifying Conserved Functional Gene Modules Underlying the Dynamic Regulation of Tea Plant Development and Secondary Metabolism

Tea plants adjust development and metabolism by integrating environmental and endogenous signals in complex but poorly defined gene networks. Here, we present an integrative analysis framework for the identification of conserved modules controlling important agronomic traits using a comprehensive collection of RNA-seq datasets in Camellia plants including 189 samples. In total, 212 secondary metabolism-, 182 stress response-, and 182 tissue development-related coexpressed modules were revealed. Functional modules (e.g., drought response, theobromine biosynthesis, and new shoot development-related modules) and potential regulators that were highly conserved across diverse genetic backgrounds and/or environmental conditions were then identified by cross-experiment comparisons and consensus clustering. Moreover, we investigate the preservation of gene networks between Camellia sinensis and other Camellia species. This revealed that the coexpression patterns of several recently evolved modules related to secondary metabolism and environmental adaptation were rewired and showed higher connectivity in tea plants. These conserved modules are excellent candidates for modeling the core mechanism of tea plant development and secondary metabolism and should serve as a great resource for hypothesis generation and tea quality improvement.

Screening of Potential Key Transcripts Involved in Planarian Regeneration and Analysis of Its Regeneration Patterns by PacBio Long-Read Sequencing

Dugesia japonica is an excellent animal model for studying the regeneration mechanism due to its characteristics of rapid regeneration and easy breeding. PacBio sequencing was performed on the intact planarians (In) and regenerating planarians of 1 day (1d), 3 days (3d), and 5 days (5d) after amputation. The aim of this study is to deeply profile the transcriptome of D. japonica and to evaluate its regenerate changes. Using robust statistical analysis, we identified 5931, 5115, and 4669 transcripts differentially expressed between 1d and In, 3d and In, 5d and In, respectively. A total of 63 key transcripts were screened from these DETs. These key transcripts enhance the expression in different regenerate stages respectively to regulate specific processes including signal transduction, mitosis, protein synthesis, transport and degradation, apoptosis, neural development, and energy cycling. Finally, according to the biological processes involved in these potential key transcripts, we propose a hypothesis of head regeneration model about D. japonica. In addition, the weighted gene co-expression network analysis provides a new way to screen key transcripts from large amounts of data. Together, these analyses identify a number of potential key regulators controlling proliferation, differentiation, apoptosis, and signal transduction. What's more, this study provides a powerful data foundation for further research on planarians regeneration.

Integrative expression network analysis of microRNA and gene isoforms in sacred lotus

Background

Gene expression is complex and regulated by multiple molecular mechanisms, such as miRNA-mediated gene inhibition and alternative-splicing of pre-mRNAs. However, the coordination of interaction between miRNAs with different splicing isoforms, and the change of splicing isoform in response to different cellular environments are largely unexplored in plants. In this study, we analyzed the miRNA and mRNA transcriptome from lotus (Nelumbo nucifera), an economically important flowering plant.

Results

Through RNA-seq analyses on miRNAs and their target genes (isoforms) among six lotus tissues, expression of most miRNAs seem to be negatively correlated with their targets and tend to be tissue-specific. Further, our results showed that preferential interactions between miRNAs and hub gene isoforms in one coexpression module which is highly correlated with leaf. Intriguingly, for many genes, their corresponding isoforms were assigned to different co-expressed modules, and they exhibited more divergent mRNA structures including presence and absence of miRNA binding sites, suggesting functional divergence for many isoforms is escalated by both structural and expression divergence. Further detailed functional enrichment analysis of miRNA targets revealed that miRNAs are involved in the regulation of lotus growth and development by regulating plant hormone-related pathway genes.

Conclusions

Taken together, our comprehensive analyses of miRNA and mRNA transcriptome elucidate the coordination of interaction between miRNAs and different splicing isoforms, and highlight the functional divergence of many transcript isoforms from the same locus in lotus.

Shading Effects on Leaf Color Conversion and Biosynthesis of the Major Secondary Metabolites in the Albino Tea Cultivar "Yujinxiang"

Albino became a novel kind of tea cultivar in China recently. In this study, transcriptome and whole-genome bisulfite sequencing (WGBS) were employed to investigate the shading effects on leaf color conversion and biosynthesis of three major secondary metabolites in the albino tea cultivar "Yujinxiang". The increased leaf chlorophyll level was likely the major cause for shaded leaf greening from young pale or yellow leaf. In comparison with the control, the total catechin level of the shading group was significantly decreased and the abundance of caffeine was markedly increased, while the theanine level was nearly not influenced. Meanwhile, differentially expressed genes (DEGs) enriched in some biological processes and pathways were identified by transcriptome analysis. Furthermore, whole-genome DNA methylation analysis revealed that the global genomic DNA methylation patterns of the shading period were remarkably altered in comparison with the control. In addition, differentially methylated regions (DMRs) and the DMR-related DEG analysis indicated that the DMR-related DEGs were the critical participants in biosynthesis of the major secondary metabolites. These findings suggest that DNA methylation is probably responsible for changes in the contents of the major secondary metabolites in Yujinxiang.

Tea plant genomics: achievements, challenges and perspectives

Tea is among the world's most widely consumed non-alcoholic beverages and possesses enormous economic, health, and cultural values. It is produced from the cured leaves of tea plants, which are important evergreen crops globally cultivated in over 50 countries. Along with recent innovations and advances in biotechnologies, great progress in tea plant genomics and genetics has been achieved, which has facilitated our understanding of the molecular mechanisms of tea quality and the evolution of the tea plant genome. In this review, we briefly summarize the achievements of the past two decades, which primarily include diverse genome and transcriptome sequencing projects, gene discovery and regulation studies, investigation of the epigenetics and noncoding RNAs, origin and domestication, phylogenetics and germplasm utilization of tea plant as well as newly developed tools/platforms. We also present perspectives and possible challenges for future functional genomic studies that will contribute to the acceleration of breeding programs in tea plants.

Optimized sequencing depth and de novo assembler for deeply reconstructing the transcriptome of the tea plant, an economically important plant species

Background

Tea is the oldest and among the world’s most popular non-alcoholic beverages, which has important economic, health and cultural values. Tea is commonly produced from the leaves of tea plants (Camellia sinensis), which belong to the genus Camellia of family Theaceae. In the last decade, many studies have generated the transcriptomes of tea plants at different developmental stages or under abiotic and/or biotic stresses to investigate the genetic basis of secondary metabolites that determine tea quality. However, these results exhibited large differences, particularly in the total number of reconstructed transcripts and the quality of the assembled transcriptomes. These differences largely result from limited knowledge regarding the optimized sequencing depth and assembler for transcriptome assembly of structurally complex plant species genomes.

Results

We employed different amounts of RNA-sequencing data, ranging from 4 to 84 Gb, to assemble the tea plant transcriptome using five well-known and representative transcript assemblers. Although the total number of assembled transcripts increased with increasing sequencing data, the proportion of unassembled transcripts became saturated as revealed by plant BUSCO datasets. Among the five representative assemblers, the Bridger package shows the best performance in both assembly completeness and accuracy as evaluated by the BUSCO datasets and genome alignment. In addition, we showed that Bridger and BinPacker harbored the shortest runtimes followed by SOAPdenovo and Trans-ABySS.

Conclusions

The present study compares the performance of five representative transcript assemblers and investigates the key factors that affect the assembly quality of the transcriptome of the tea plants. This study will be of significance in helping the tea research community obtain better sequencing and assembly of tea plant transcriptomes under conditions of interest and may thus help to answer major biological questions currently facing the tea industry.