De novo assembly and transcriptome characterization: novel insights into catechins biosynthesis in Camellia sinensis

Transcriptome analysis of immature xylem in the Chinese fir at different developmental phases

Background.Chinese fir [Cunninghamia lanceolata (Lamb.) Hook.] is one of the most important native tree species for timber production in southern China. An understanding of overall fast growing stage, stem growth stage and senescence stage cambium transcriptome variation is lacking. We used transcriptome sequencing to identify the repertoire of genes expressed during development of xylem tissue in Chinese fir, aiming to delineate the molecular mechanisms of wood formation. Results. We carried out transcriptome sequencing at three different cultivation ages (7Y, 15Y and 21Y) generating 68.71 million reads (13.88 Gbp). A total of 140,486 unigenes with a mean size of 568.64 base pairs (bp) were obtained via de novo assembly. Of these, 27,427 unigenes (19.52%) were further annotated by comparison to public protein databases. A total of 5,331 (3.79%) unigenes were mapped into 118 pathways by searching against the Kyoto Encyclopedia of Genes and Genomes Pathway database (KEGG). Differentially expressed genes (DEG) analysis identified 3, 16 and 5,899 DEGs from the comparison of 7Y vs. 15Y, 7Y vs. 21Y and 15Y vs. 21Y, respectively, in the immature xylem tissues, including 2,638 significantly up-regulated and 3,280 significantly down-regulated genes. Besides, five NAC transcription factors, 190 MYB transcription factors, and 34 WRKY transcription factors were identified respectively from Chinese fir transcriptome. Conclusion. Our results revealed the active transcriptional pathways and identified the DEGs at different cultivation phases of Chinese fir wood formation. This transcriptome dataset will aid in understanding and carrying out future studies on the molecular basis of Chinese fir wood formation and contribute to future artificial production and applications.

The Dark-Purple Tea Cultivar 'Ziyan' Accumulates a Large Amount of Delphinidin-Related Anthocyanins

Recently, we developed a novel tea cultivar 'Ziyan' with distinct purple leaves. There was a significant correlation between leaf color and anthocyanin pigment content in the leaves. A distinct allocation of metabolic flow for B-ring trihydroxylated anthocyanins and catechins in 'Ziyan' was observed. Delphinidin, cyanidin, and pelargonidin (88.15 mg/100 g FW in total) but no other anthocyanin pigments were detected in 'Ziyan', and delphinidin (70.76 mg/100 g FW) was particularly predominant. An analysis of the catechin content in 'Ziyan' and eight other cultivars indicated that 'Ziyan' exhibits a preference for synthesizing B-ring trihydroxylated catechins (with a proportion of 74%). The full-length cDNA sequences of flavonoid pathway genes were isolated by RNA-Seq coupled with conventional TA cloning, and their expression patterns were characterized. Purple-leaved cultivars had lower amounts of total catechins, polyphenols, and water extract than ordinary non-anthocyanin cultivars but similar levels of caffeine. Because dark-purple-leaved Camellia species are rare in nature, this study provides new insights into the interplay between the accumulations of anthocyanins and other bioactive components in tea leaves.

Transcriptomic profiling of the salt stress response in excised leaves of the halophyte Beta vulgaris ssp. maritima

Beta vulgaris ssp. maritima is a halophytic relative of cultivated beets. In the present work a transcriptome response to acute salt stress imposed to excised leaves of sea beet was investigated. Salt treatments consisted of adding NaCl directly to the transpiration stream by immersing the petioles of excised leaves into the salt solutions. Sequencing libraries were generated from leaves subjected to either moderate or strong salt stress. Control libraries were constructed from untreated leaves. Sequencing was performed using the Illumina MiSeq platform. We obtained 32970 unigenes by assembling the pooled reads from all the libraries with Trinity software. Screening the nr database returned 18,362 sequences with functional annotation. Using the reference transcriptome we identified 1,246 genes that were differentially expressed after 48 h of NaCl stress. Genes related to several cellular functions such as membrane transport, osmoprotection, molecular chaperoning, redox metabolism or protein synthesis were differentially expressed in response to salt stress. The response of sea beet leaves to salt treatments was marked out by transcriptomic up-regulation of genes related to photosynthetic carbon fixation, ribosome biogenesis, cell wall-building and cell wall expansion. Furthermore, several novel and undescribed transcripts were responsive to salinity in leaves of sea beet.

Selection of suitable reference genes for qRT-PCR normalization during leaf development and hormonal stimuli in tea plant (Camellia sinensis)

Tea plant (Camellia sinensis) leaf is an important non-alcoholic beverage resource. The application of quantitative real time polymerase chain reaction (qRT-PCR) has a profound significance for the gene expression studies of tea plant, especially when applied to tea leaf development and metabolism. In this study, nine candidate reference genes (i.e., CsACT7, CsEF-1α, CseIF-4α, CsGAPDH, CsPP2A, CsSAND, CsTBP, CsTIP41, and CsTUB) of C. sinensis were cloned. The quantitative expression data of these genes were investigated in five tea leaf developmental stages (i.e., 1st, 2nd, 3rd, 4th, and older leaves) and normal growth tea leaves subjected to five hormonal stimuli (i.e., ABA, GA, IAA, MeJA, and SA), and gene expression stability was calculated using three common statistical algorithms, namely, geNorm, NormFinder, and Bestkeeper. Results indicated that CsTBP and CsTIP41 were the most stable genes in tea leaf development and CsTBP was the best gene under hormonal stimuli; by contrast, CsGAPDH and CsTUB genes showed the least stability. The gene expression profile of CsNAM gene was analyzed to confirm the validity of the reference genes in this study. Our data provide basis for the selection of reference genes for future biological research in the leaf development and hormonal stimuli of C. sinensis.

De novo Transcriptome Analysis and Molecular Marker Development of Two Hemarthria Species

Hemarthria R. Br. is an important genus of perennial forage grasses that is widely used in subtropical and tropical regions. Hemarthria grasses have made remarkable contributions to the development of animal husbandry and agro-ecosystem maintenance; however, there is currently a lack of comprehensive genomic data available for these species. In this study, we used Illumina high-throughput deep sequencing to characterize of two agriculturally important Hemarthria materials, H. compressa "Yaan" and H. altissima "1110." Sequencing runs that used each of four normalized RNA samples from the leaves or roots of the two materials yielded more than 24 million high-quality reads. After de novo assembly, 137,142 and 77,150 unigenes were obtained for "Yaan" and "1110," respectively. In addition, a total of 86,731 "Yaan" and 48,645 "1110" unigenes were successfully annotated. After consolidating the unigenes for both materials, 42,646 high-quality SNPs were identified in 10,880 unigenes and 10,888 SSRs were identified in 8330 unigenes. To validate the identified markers, high quality PCR primers were designed for both SNPs and SSRs. We randomly tested 16 of the SNP primers and 54 of the SSR primers and found that the majority of these primers successfully amplified the desired PCR product. In addition, high cross-species transferability (61.11-87.04%) of SSR markers was achieved for four other Poaceae species. The amount of RNA sequencing data that was generated for these two Hemarthria species greatly increases the amount of genomic information available for Hemarthria and the SSR and SNP markers identified in this study will facilitate further advancements in genetic and molecular studies of the Hemarthria genus.

Cs-miR156 is involved in the nitrogen form regulation of catechins accumulation in tea plant (Camellia sinensis L.)

The nitrogen source affects the growth of tea plants and regulates the accumulation of catechins in the leaves. In this report, we assessed the influences of NH4(+) and NO3(-) on plant growth, catechins accumulation and associated gene expression. Compared with the preferential nitrogen source NH4(+), when NO3(-) was supplied as the sole nitrogen source, tea plants showed similar symptoms with the nitrogen-free treatments and showed lower nitrogen, free amino acid accumulation, chlorophyll content and biomass gain, indicating NO3(-) was not efficiently used by these plants. However, the total shoot catechins content was significantly higher for NO3(-) treatments than that for NH4(+) treatment or combined NH4(+)+NO3(-) treatment, suggesting that, in addition to its influence on plant growth, the nitrogen form regulated the accumulation of catechins in tea. The expression of catechins biosynthesis-related genes was associated with the regulation of catechins accumulation and composition changes mediated by nitrogen form. PAL, CHS, CHI, and DFR genes exhibited higher expression levels in plants supplied with NO3(-), in which the transcript level of DFR in the shoots was significantly correlated with the catechins content. In the end, we identified a new function for the Cs-miR156, which was drastically induced through NH4(+). Moreover, a potential mechanism of the Cs-miR156 pathway in regulating catechins biosynthesis in tea plants has been suggested, with particular respect to nitrogen forms. Cs-miR156 might repress the expression of the target gene SPL to regulate the DFR gene, which plays a vital role in catechins biosynthesis.

Identification, classification, and expression profiles of heat shock transcription factors in tea plant (Camellia sinensis) under temperature stress

In vascular plants, heat shock transcription factors (Hsfs) regulate heat stress response by regulating the expression of heat shock proteins. This study systematically and comprehensively analyzed the Hsf family in tea plant [Camellia sinensis (L.) O. Kuntze]. A total of 16 CsHsfs were identified from the transcriptome database of tea plant and analyzed for their phylogenetic relationships, motifs, and physicochemical characteristics. On the basis of the phylogenetic comparison of tea plant with Arabidopsis thaliana, Populus trichocarpa, Theobroma cacao, and Oryza sativa, the CsHsfs were classified into three classes, namely, A (56.25%), B (37.50%), and C (6.25%). Heat mapping showed that the expression profiles of CsHsf genes under non-stress conditions varied among four tea plant cultivars, namely, 'Yunnanshilixiang', 'Chawansanhao', 'Ruchengmaoyecha', and 'Anjibaicha'. Six CsHsf genes (CsHsfA1a, CsHsfA1b, CsHsfA6, CsHsfB1, CsHsfB2b, and CsHsfC1) were selected from classes A, B, and C to analyze the expression profiles of CsHsf genes through quantitative real-time PCR in 'Yingshuang', 'Anjibaicha', and 'Yunnanshilixiang' under high (38 °C) or low (4 °C) temperature stress. Temperature stress positively or negatively regulated all of the selected CsHsf genes, and the expression levels evidently varied even among CsHsf genes belonging to the same class. This study provided a relatively detailed summary of Hsfs in tea plant and may serve as a reference for further studies on the mechanism of temperature stress regulation by CsHsfs.

Transcriptome Analysis Reveals Key Flavonoid 3'-Hydroxylase and Flavonoid 3',5'-Hydroxylase Genes in Affecting the Ratio of Dihydroxylated to Trihydroxylated Catechins in Camellia sinensis

The ratio of dihydroxylated to trihydroxylated catechins (RDTC) is an important indicator of tea quality and biochemical marker for the study of genetic diversity. It is reported to be under genetic control but the underlying mechanism is not well understood. Flavonoid 3′-hydroxylase (F3′H) and flavonoid 3′,5′-hydroxylase (F3′5′H) are key enzymes involved in the formation of dihydroxylated and trihydroxylated catechins. The transcriptome and HPLC analysis of tea samples from Longjing43 and Zhonghuang2 under control and shading treatment were performed to assess the F3′H and F3′5′H genes that might affect RDTC. A total of 74.7 million reads of mRNA seq (2×101bp) data were generated. After de novo assembly, 109,909 unigenes were obtained, and 39,982 of them were annotated using 7 public databases. Four key F3′H and F3′5′H genes (including CsF3′5′H1, CsF3′H1, CsF3′H2 and CsF3′H3) were identified to be closely correlated with RDTC. Shading treatment had little effect on RDTC, which was attributed to the stable expression of these key F3′H and F3′5′H genes. The correlation of the coexpression of four key genes and RDTC was further confirmed among 13 tea varieties by real time PCR and HPLC analysis. The coexpression of three F3′H genes and a F3′5′H gene may play a key role in affecting RDTC in Camellia sinensis. The current results may establish valuable foundation for further research about the mechanism controlling catechin composition in tea.

Transcriptome-wide identification of Camellia sinensis WRKY transcription factors in response to temperature stress

Tea plant [Camellia sinensis (L.) O. Kuntze] is a leaf-type healthy non-alcoholic beverage crop, which has been widely introduced worldwide. Tea is rich in various secondary metabolites, which are important for human health. However, varied climate and complex geography have posed challenges for tea plant survival. The WRKY gene family in plants is a large transcription factor family that is involved in biological processes related to stress defenses, development, and metabolite synthesis. Therefore, identification and analysis of WRKY family transcription factors in tea plant have a profound significance. In the present study, 50 putative C. sinensis WRKY proteins (CsWRKYs) with complete WRKY domain were identified and divided into three Groups (Group I-III) on the basis of phylogenetic analysis results. The distribution of WRKY family transcription factors among plantae, fungi, and protozoa showed that the number of WRKY genes increased in higher plant, whereas the number of these genes did not correspond to the evolutionary relationships of different species. Structural feature and annotation analysis results showed that CsWRKY proteins contained WRKYGQK/WRKYGKK domains and C2H2/C2HC-type zinc-finger structure: D-X18-R-X1-Y-X2-C-X4-7-C-X23-H motif; CsWRKY proteins may be associated with the biological processes of abiotic and biotic stresses, tissue development, and hormone and secondary metabolite biosynthesis. Temperature stresses suggested that the candidate CsWRKY genes were involved in responses to extreme temperatures. The current study established an extensive overview of the WRKY family transcription factors in tea plant. This study also provided a global survey of CsWRKY transcription factors and a foundation of future functional identification and molecular breeding.

Transcriptomic and phytochemical analysis of the biosynthesis of characteristic constituents in tea (Camellia sinensis) compared with oil tea (Camellia oleifera)

BACKGROUND

Tea plants (Camellia sinensis) are used to produce one of the most important beverages worldwide. The nutritional value and healthful properties of tea are closely related to the large amounts of three major characteristic constituents including polyphenols (mainly catechins), theanine and caffeine. Although oil tea (Camellia oleifera) belongs to the genus Camellia, this plant lacks these three characteristic constituents. Comparative analysis of tea and oil tea via RNA-Seq would help uncover the genetic components underlying the biosynthesis of characteristic metabolites in tea.

RESULTS

We found that 3,787 and 3,359 bud genes, as well as 4,042 and 3,302 leaf genes, were up-regulated in tea and oil tea, respectively. High-performance liquid chromatography (HPLC) analysis revealed high levels of all types of catechins, theanine and caffeine in tea compared to those in oil tea. Activation of the genes involved in the biosynthesis of these characteristic compounds was detected by RNA-Seq analysis. In particular, genes encoding enzymes involved in flavonoid, theanine and caffeine pathways exhibited considerably different expression levels in tea compared to oil tea, which were also confirmed by quantitative RT-PCR (qRT-PCR).

CONCLUSION

We assembled 81,826 and 78,863 unigenes for tea and oil tea, respectively, based on their differences at the transcriptomic level. A potential connection was observed between gene expression and content variation for catechins, theanine and caffeine in tea and oil tea. The results demonstrated that the metabolism was activated during the accumulation of characteristic metabolites in tea, which were present at low levels in oil tea. From the molecular biological perspective, our comparison of the transcriptomes and related metabolites revealed differential regulatory mechanisms underlying secondary metabolic pathways in tea versus oil tea.

Transcriptome-based discovery of AP2/ERF transcription factors related to temperature stress in tea plant (Camellia sinensis)

Tea plant (Camellia sinensis) is an important natural resource for the global supply of non-alcoholic beverage production. The extension of tea plant cultivation is challenged by biotic and abiotic stresses. Transcription factors (TFs) of the APETALA 2 (AP2)/ethylene-responsive factor (ERF) family are involved in growth and anti-stresses through multifaceted transcriptional regulation in plants. This study comprehensively analyzed AP2/ERF family TFs from C. sinensis on the basis of the transcriptome sequencing data of four tea plant cultivars, namely, 'Yunnanshilixiang', 'Chawansanhao', 'Ruchengmaoyecha', and 'Anjibaicha'. A total of 89 putative AP2/ERF transcription factors with full-length AP2 domain were identified from C. sinensis and classified into five subfamilies, namely, AP2, dehydration-responsive-element-binding (DREB), ERF, related to ABI3/VP (RAV), and Soloist. All identified CsAP2/ERF genes presented relatively stable expression levels in the four tea plant cultivars. Many groups also showed cultivar specificity. Five CsAP2/ERF genes from each AP2/ERF subfamily (DREB, ERF, AP2, and RAV) were related to temperature stresses; these results indicated that AP2/ERF TFs may play important roles in abnormal temperature stress response in C. sinensis.

Global transcriptome and gene regulation network for secondary metabolite biosynthesis of tea plant (Camellia sinensis)

Background

Major secondary metabolites, including flavonoids, caffeine, and theanine, are important components of tea products and are closely related to the taste, flavor, and health benefits of tea. Secondary metabolite biosynthesis in Camellia sinensis is differentially regulated in different tissues during growth and development. Until now, little was known about the expression patterns of genes involved in secondary metabolic pathways or their regulatory mechanisms. This study aimed to generate expression profiles for C. sinensis tissues and to build a gene regulation model of the secondary metabolic pathways.

Results

RNA sequencing was performed on 13 different tissue samples from various organs and developmental stages of tea plants, including buds and leaves of different ages, stems, flowers, seeds, and roots. A total of 43.7 Gbp of raw sequencing data were generated, from which 347,827 unigenes were assembled and annotated. There were 46,693, 8446, 3814, 10,206, and 4948 unigenes specifically expressed in the buds and leaves, stems, flowers, seeds, and roots, respectively. In total, 1719 unigenes were identified as being involved in the secondary metabolic pathways in C. sinensis, and the expression patterns of the genes involved in flavonoid, caffeine, and theanine biosynthesis were characterized, revealing the dynamic nature of their regulation during plant growth and development. The possible transcription factor regulation network for the biosynthesis of flavonoid, caffeine, and theanine was built, encompassing 339 transcription factors from 35 families, namely bHLH, MYB, and NAC, among others. Remarkably, not only did the data reveal the possible critical check points in the flavonoid, caffeine, and theanine biosynthesis pathways, but also implicated the key transcription factors and related mechanisms in the regulation of secondary metabolite biosynthesis.

Conclusions

Our study generated gene expression profiles for different tissues at different developmental stages in tea plants. The gene network responsible for the regulation of the secondary metabolic pathways was analyzed. Our work elucidated the possible cross talk in gene regulation between the secondary metabolite biosynthetic pathways in C. sinensis. The results increase our understanding of how secondary metabolic pathways are regulated during plant development and growth cycles, and help pave the way for genetic selection and engineering for germplasm improvement.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1773-0) contains supplementary material, which is available to authorized users.

Transcriptome Analysis Reveals Key Flavonoid 3'-Hydroxylase and Flavonoid 3',5'-Hydroxylase Genes in Affecting the Ratio of Dihydroxylated to Trihydroxylated Catechins in Camellia sinensis

The ratio of dihydroxylated to trihydroxylated catechins (RDTC) is an important indicator of tea quality and biochemical marker for the study of genetic diversity. It is reported to be under genetic control but the underlying mechanism is not well understood. Flavonoid 3'-hydroxylase (F3'H) and flavonoid 3',5'-hydroxylase (F3'5'H) are key enzymes involved in the formation of dihydroxylated and trihydroxylated catechins. The transcriptome and HPLC analysis of tea samples from Longjing43 and Zhonghuang2 under control and shading treatment were performed to assess the F3'H and F3'5'H genes that might affect RDTC. A total of 74.7 million reads of mRNA seq (2×101bp) data were generated. After de novo assembly, 109,909 unigenes were obtained, and 39,982 of them were annotated using 7 public databases. Four key F3'H and F3'5'H genes (including CsF3'5'H1, CsF3'H1, CsF3'H2 and CsF3'H3) were identified to be closely correlated with RDTC. Shading treatment had little effect on RDTC, which was attributed to the stable expression of these key F3'H and F3'5'H genes. The correlation of the coexpression of four key genes and RDTC was further confirmed among 13 tea varieties by real time PCR and HPLC analysis. The coexpression of three F3'H genes and a F3'5'H gene may play a key role in affecting RDTC in Camellia sinensis. The current results may establish valuable foundation for further research about the mechanism controlling catechin composition in tea.

Comparative transcriptome analysis to reveal genes involved in wheat hybrid necrosis

Wheat hybrid necrosis is an interesting genetic phenomenon that is found frequently and results in gradual death or loss of productivity of wheat. However, the molecular basis and mechanisms of this genetic phenomenon are still not well understood. In this study, the transcriptomes of wheat hybrid necrosis F₁ and its parents (Neimai 8 and II469) were investigated using digital gene expression (DGE). A total of 1300 differentially expressed genes were identified, indicating that the response to hybrid necrosis in wheat is complicated. The assignments of the annotated genes based on Gene Ontology (GO) revealed that most of the up-regulated genes belong to “universal stress related”, “DNA/RNA binding”, “protein degradation” functional groups, while the down-regulated genes belong to “carbohydrate metabolism” and “translation regulation” functional groups. These findings suggest that these pathways were affected by hybrid necrosis. Our results provide preliminarily new insight into the underlying molecular mechanisms of hybrid necrosis and will help to identify important candidate genes involved in wheat hybrid necrosis.