Transcriptome analysis of the oil-rich tea plant, Camellia oleifera, reveals candidate genes related to lipid metabolism

Chromosome-scale genome assembly of oil-tea tree Camellia crapnelliana

Camellia crapnelliana Tutch., belonging to the Theaceae family, is an excellent landscape tree species with high ornamental values. It is particularly an important woody oil-bearing plant species with high ecological, economic, and medicinal values. Here, we first report the chromosome-scale reference genome of C. crapnelliana with integrated technologies of SMRT, Hi-C and Illumina sequencing platforms. The genome assembly had a total length of ~2.94 Gb with contig N50 of ~67.5 Mb, and ~96.34% of contigs were assigned to 15 chromosomes. In total, we predicted 37,390 protein-coding genes, ~99.00% of which could be functionally annotated. The chromosome-scale genome of C. crapnelliana will become valuable resources for understanding the genetic basis of the fatty acid biosynthesis, and greatly facilitate the exploration and conservation of C. crapnelliana.

Genome-wide identification and comparative analysis of YABBY transcription factors in oil tea and tea tree

The plant-specific transcription factor gene family, YABBY, plays an important role in plant development and stress response. Although YABBY genes have been identified in numerous species, a comprehensive characterization of YABBYs in tea tree and oil tea has been lacking. In this study, ten and three YABBY genes were identified in Camellia sinensis and C. oleifera, respectively. YABBY proteins could be divided into five subfamilies. Most YABBY genes in the same clade had similar structures and conserved motifs. Protein evolutionary analysis revealed that FIL/YAB3 displayed high conservation in all positions, followed by INO, YAB2, YAB5, and CRC. Specific site analysis suggested that the YABBY family was polyphyletic during the evolution. Compared to C. oleifera, two segmentally duplicated gene pairs were formed in C. sinensis during recent WGD events generated 30.69 and 45.08 Mya, respectively. Cis-acting element indicated that most YABBY genes contain box4, ARE, and MYB elements. A total of 120 SSR loci were found within CsYABBYs, consisting of six types, while 48 SSR loci were identified within CoYABBY, consisting of three types. Transcriptome results revealed that CRC and INO clades were specifically expressed in floral organs. The expression of CsYABBY10 and CsYABBY5 was significantly up-regulated under drought and salt treatments, respectively, as confirmed by qRT-PCR. CoYABBY genes were more susceptible to salt stress, as CoYABBY3 increased by about 15-fold. Furthermore, functional differentiation may have occurred in duplicated genes. These discoveries provide important information for further research on YABBYs in tea tree and oil tea.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-024-03940-9.

Identification of QTLs associated with very-long chain fatty acid (VLCFA) content via linkage mapping and BSA-seq in peanut

KEY MESSAGE

Three major QTLs qA01, qB04.1 and qB05 for VLCFA content and their corresponding allele-specific markers will benefit peanut low VLCFA breeding, and a candidate gene Arahy.IF1JV3 was predicted. Peanut is a globally significant oilseed crop worldwide, and contains a high content (20%) of saturated fatty acid (SFA) in its seeds. As high level SFA intake in human dietary may increase the cardiovascular disease risk, reducing the SFA content in peanut is crucial for improving its nutritional quality. Half of the SFAs in peanut are very long-chain fatty acids (VLCFA), so reducing the VLCFA content is a feasible strategy to decrease the total SFA content. Luoaowan with extremely low VLCFA (4.80%) was crossed with Jihua16 (8.00%) to construct an F2:4 population. Three major QTLs including qA01, qB04.1 and qB05 for VLCFA content were detected with 4.43 ~ 14.32% phenotypic variation explained through linkage mapping. Meanwhile, three genomic regions on chromosomes B03, B04 and B05 were identified via BSA-seq approach. Two co-localized intervals on chromosomes B04 (100.10 ~ 103.97 Mb) and B05 (6.39 ~ 10.90 Mb) were identified. With markers developed based on SNP/InDel variations in qA01 between the two parents, the remaining interval was refined to 103.58 ~ 111.14 Mb. A candidate gene Arahy.IF1JV3 encoding a β-ketoacyl-CoA synthase was found in qA01, and its expression level in Luoaowan was significantly lower than that in Jihua16. Allele-specific markers targeting qA01, qB04.1 and qB05 were developed and validated in F4 population, and an elite line with high oleic, low VLCFA (5.05%) and low SFA (11.48%) contents was selected. This study initially revealed the genetic mechanism of VLCFA content, built a marker-assisted selection system for low VLCFA breeding, and provided an effective method to decrease the SFA content in peanut.

Integrative Metabolic and Transcriptomic Profiling in Camellia oleifera and Camellia meiocarpa Uncover Potential Mechanisms That Govern Triacylglycerol Degradation during Seed Desiccation

Camellia seed oil is a top-end quality of cooking oil in China. The oil quality and quantity are formed during seed maturation and desiccation. So far, it remains largely unresolved whether lipid degradation occurs and contributes to Camellia oil traits. In this study, three different Camellia germplasms, C. oleifera cv. Min 43 (M43), C. meiocarpa var. Qingguo (QG), and C. meiocarpa cv Hongguo (HG) were selected, their seed oil contents and compositions were quantified across different stages of seed desiccation. We found that at the late stage of desiccation, M43 and QG lost a significant portion of seed oil, while such an event was not observed in HG. To explore the molecular bases for the oil loss In M43, the transcriptomic profiling of M43 and HG was performed at the early and the late seed desiccation, respectively, and differentially expressed genes (DEGs) from the lipid metabolic pathway were identified and analyzed. Our data demonstrated that different Camellia species have diverse mechanisms to regulate seed oil accumulation and degradation, and that triacylglycerol-to-terpenoid conversion could account for the oil loss in M43 during late seed desiccation.

Genomic and genetic advances of oiltea-camellia (Camellia oleifera)

Oiltea-camellia (C. oleifera) is a widely cultivated woody oil crop in Southern China and Southeast Asia. The genome of oiltea-camellia was very complex and not well explored. Recently, genomes of three oiltea-camellia species were sequenced and assembled, multi-omic studies of oiltea-camellia were carried out and provided a better understanding of this important woody oil crop. In this review, we summarized the recent assembly of the reference genomes of oiltea-camellia, genes related to economic traits (flowering, photosynthesis, yield and oil component), disease resistance (anthracnose) and environmental stress tolerances (drought, cold, heat and nutrient deficiency). We also discussed future directions of integrating multiple omics for evaluating genetic resources and mining key genes of important traits, and the application of new molecular breeding and gene editing technologies to accelerate the breeding process of oiltea-camellia.

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.

Camellia oil Enhances Plasma Antioxidant Metabolism and Improves Plasma Lipid Metabolism in High-fat Diet-fed Rats

Living on a high-fat, high-calorie, and high-protein diet for a long period may compromise human immunity due to the long-term accumulation of free radicals and plasma lipids. The antioxidant and lipid-lowering compounds (ie polyphenols and vitamin E) in Camellia oil help to decrease the risk of numerous ailments, including cardiovascular disease (CVD), and obesity.

The aims of this study were to study the hypolipidemic and antioxidant effects of Camellia oil in high-fat-fed rats and to promote the high-value use of camellia resources. The high-fat-fed rats were administrated with 2.5, 7.5, and 15 mL/kg BW Camellia oil (Camellia oil group), and 10 mg/kg BW atorvastatin (atorvastatin group), respectively, and compared with a model group (only fed with high fat) and a control group (fed with basal diet). Several parameters were measured, including (1) body weight (BW), liver-to-BW ratio; (2) plasma total cholesterol (TC), triglyceride (TG), high-density lipoprotein cholesterol (HDL-C), and low-density lipoprotein cholesterol (LDL-C); and (3) alanine aminotransferase (ALT), alanine aminotransferase (AST), superoxide dismutase (SOD), glutathione peroxidase (GSH-Px) activity, model driven architecture (MDA) content, lipid metabolism-related genes, and antioxidant-related genes in liver tissue.

Compared with the model group, the high-fat-fed rats in the Camellia oil and atorvastatin group had significantly lower BW and liver-to-BW ratio ( P < .01), plasma TC, TG, and LDL-C levels and ALT and AST activities, but higher HDL-C levels. The relative expressions of ACAT1, DGAT2, FAS, and SREBP genes were significantly reduced in the Camellia oil and atorvastatin groups, while the relative expressions of LCAT, UCP2, MCD, and CPT-1 genes were significantly increased. The rats in the Camellia oil group showed significantly higher SOD and GSH-Px activities, significantly lower MDA content, and significantly higher relative expression of antioxidant genes (eg SOD1, GPx1, CAT, and Gclm). Thus, atorvastatin and Camellia oil exhibited significant hypolipidemic and antioxidant effects, which were better at a dose of 7.5 mL/kg (BW) of Camellia oil. Therefore, Camellia oil becomes a potential new natural resource for future research and development of antioxidant and hypolipidemic drugs, nutraceuticals, and additives.

Comparative transcriptomic analysis reveals genes related to the rapid accumulation of oleic acid in Camellia chekiangoleosa, an oil tea plant with early maturity and large fruit

Camellia chekiangoleosa has a higher oleic acid content and a shorter reproductive cycle than typical oil tea plants. It was intensively sampled over six C. chekiangoleosa seed development stages. The content of fatty acids determined by GC showed that the accumulation of fatty acids gradually increased from the S1 to S5 stages, and the maximum concentration was reached in S5. Then, fatty acids declined slightly in S6. The main fatty acid component showed the same accumulation trend as the total fatty acids, except linolenic acid, which remained at a low level throughout seed developmental stages. Changes in the expression of fatty acid accumulation-related genes were monitored using second-generation and SMRT full-length transcriptome sequencing. Finally, 18.92 G accurate and reliable data were obtained. Differential expression analysis and weighted coexpression analysis revealed two "gene modules" significantly associated with oleic acid and linoleic acid contents, and the high expression of ENR, KAS I, and KAS II, which accumulate substrates for oleic acid synthesis, was thought to be responsible for the rapid accumulation of fatty acids in the early stage. The rapid increase in fatty acids in the second stage may be closely related to the synergy between the high expression of SAD and low expression of FAD2. In addition, many transcription factors, such as ERF, GRAS, GRF, MADS, MYB and WRKY, may be involved in the fatty acid synthesis. Our data provide a rich resource for further studies on the regulation of fatty acid synthesis in C. chekiangoleosa.

The complete chloroplast genome of Camellia osmantha, an edible oil Camellia

Camellia osmantha is a new species of the Camellia genus discovered in Nanning, Guangxi, China, in 2012. It can be used as an excellent woody oil crop. There is little related research on this species in China and abroad, and its genome information is still lacking. In this study, the complete chloroplast genome sequence of C. osmantha was first reported (GenBank number: MZ128138). The whole chloroplast genome is 156,981 bp in length with a GC content of 37.28%, and it is composed of a large single copy (LSC) region of 86,647 bp, a small single copy (SSC) region of 18,284 bp, and a pair of inverted repeat (IR) regions of 26,025 bp each. The genome contains a total of 135 functional genes, including 37 transfer RNA genes, 90 protein-coding genes, and 8 ribosomal RNA genes. The maximum likelihood analysis based on 21 chloroplast genomes showed that C. osmantha and C. oleifera (MF541730.2) were the most closely related.

Integrative iTRAQ-based proteomic and transcriptomic analysis reveals the accumulation patterns of key metabolites associated with oil quality during seed ripening of Camellia oleifera

Camellia oleifera (C. oleifera) is one of the four major woody oil-bearing crops in the world and has relatively high ecological, economic, and medicinal value. Its seeds undergo a series of complex physiological and biochemical changes during ripening, which is mainly manifested as the accumulation and transformation of certain metabolites closely related to oil quality, especially flavonoids and fatty acids. To obtain new insights into the underlying molecular mechanisms, a parallel analysis of the transcriptome and proteome profiles of C. oleifera seeds at different maturity levels was conducted using RNA sequencing (RNA-seq) and isobaric tags for relative and absolute quantification (iTRAQ) complemented with gas chromatography-mass spectrometry (GC-MS) data. A total of 16,530 transcripts and 1228 proteins were recognized with significant differential abundances in pairwise comparisons of samples at various developmental stages. Among these, 317 were coexpressed with a poor correlation, and most were involved in metabolic processes, including fatty acid metabolism, α-linolenic acid metabolism, and glutathione metabolism. In addition, the content of total flavonoids decreased gradually with seed maturity, and the levels of fatty acids generally peaked at the fat accumulation stage; these results basically agreed with the regulation patterns of genes or proteins in the corresponding pathways. The expression levels of proteins annotated as upstream candidates of phenylalanine ammonia-lyase (PAL) and chalcone synthase (CHS) as well as their cognate transcripts were positively correlated with the variation in the flavonoid content, while shikimate O-hydroxycinnamoyltransferase (HCT)-encoding genes had the opposite pattern. The increase in the abundance of proteins and mRNAs corresponding to alcohol dehydrogenase (ADH) was associated with a reduction in linoleic acid synthesis. Using weighted gene coexpression network analysis (WGCNA), we further identified six unique modules related to flavonoid, oil, and fatty acid anabolism that contained hub genes or proteins similar to transcription factors (TFs), such as MADS intervening keratin-like and C-terminal (MIKC_MADS), type-B authentic response regulator (ARR-B), and basic helix-loop-helix (bHLH). Finally, based on the known metabolic pathways and WGCNA combined with the correlation analysis, five coexpressed transcripts and proteins composed of cinnamyl-alcohol dehydrogenases (CADs), caffeic acid 3-O-methyltransferase (COMT), flavonol synthase (FLS), and 4-coumarate: CoA ligase (4CL) were screened out. With this exploratory multiomics dataset, our results presented a dynamic picture regarding the maturation process of C. oleifera seeds on Hainan Island, not only revealing the temporal specific expression of key candidate genes and proteins but also providing a scientific basis for the genetic improvement of this tree species.

Identification of genes associated with the biosynthesis of unsaturated fatty acid and oil accumulation in herbaceous peony 'Hangshao' (Paeonia lactiflora 'Hangshao') seeds based on transcriptome analysis

Background

Paeonia lactiflora ‘Hangshao’ is widely cultivated in China as a traditional Chinese medicine ‘Radix Paeoniae Alba’. Due to the abundant unsaturated fatty acids in its seed, it can also be regarded as a new oilseed plant. However, the process of the biosynthesis of unsaturated fatty acids in it has remained unknown. Therefore, transcriptome analysis is helpful to better understand the underlying molecular mechanisms.

Results

Five main fatty acids were detected, including stearic acid, palmitic acid, oleic acid, linoleic acid and α-linolenic acid, and their absolute contents first increased and then decreased during seed development. A total of 150,156 unigenes were obtained by transcriptome sequencing. There were 15,005 unigenes annotated in the seven functional databases, including NR, NT, GO, KOG, KEGG, Swiss-Prot and InterPro. Based on the KEGG database, 1766 unigenes were annotated in the lipid metabolism. There were 4635, 12,304, and 18,291 DEGs in Group I (60 vs 30 DAF), Group II (90 vs 60 DAF) and Group III (90 vs 30 DAF), respectively. A total of 1480 DEGs were detected in the intersection of the three groups. In 14 KEGG pathways of lipid metabolism, 503 DEGs were found, belonging to 111 enzymes. We screened out 123 DEGs involved in fatty acid biosynthesis (39 DEGs), fatty acid elongation (33 DEGs), biosynthesis of unsaturated fatty acid (24 DEGs), TAG assembly (17 DEGs) and lipid storage (10 DEGs). Furthermore, qRT-PCR was used to analyze the expression patterns of 16 genes, including BBCP, BC, MCAT, KASIII, KASII, FATA, FATB, KCR, SAD, FAD2, FAD3, FAD7, GPAT, DGAT, OLE and CLO, most of which showed the highest expression at 45 DAF, except for DGAT, OLE and CLO, which showed the highest expression at 75 DAF.

Conclusions

We predicted that MCAT, KASIII, FATA, SAD, FAD2, FAD3, DGAT and OLE were the key genes in the unsaturated fatty acid biosynthesis and oil accumulation in herbaceous peony seed. This study provides the first comprehensive genomic resources characterizing herbaceous peony seed gene expression at the transcriptional level. These data lay the foundation for elucidating the molecular mechanisms of fatty acid biosynthesis and oil accumulation for herbaceous peony.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07339-7.

Multiple-genotypes transcriptional analysis revealed candidates genes and nucleotide variants for improvement of quality characteristics in tea (Camellia sinensis (L.) O. Kuntze)

Tea quality is a polygenic trait that exhibits tremendous genetic variability due to accumulation of array of secondary metabolites. To elucidate global molecular insights controlling quality attributes, metabolite profiling and transcriptome sequencing of twelve diverse tea cultivars was performed in tea shoots harvested during quality season. RP-HPLC-DAD analysis of quality parameters revealed significant difference in catechins, theanine and caffeine contents. Transcriptome sequencing resulted into 50,107 non-redundant transcripts with functional annotations of 81.6% (40,847) of the transcripts. Interestingly, 2872 differentially expressed transcripts exhibited significant enrichment in 38 pathways (FDR ≤ 0.05) including secondary metabolism, amino acid and carbon metabolism. Thirty-eight key candidates reportedly involved in biosynthesis of fatty acid derived volatiles, volatile terpenes, glycoside hydrolysis and key quality related pathways (flavonoid, caffeine and theanine-biosynthesis) were highly expressed in catechins-rich tea cultivars. Furthermore, enrichment of candidates involved in flavonoid biosynthesis, transcriptional regulation, volatile terpene and biosynthesis of fatty acid derived volatile in Protein-Protein Interactome network revealed well-coordinated regulation of quality characteristics in tea. Additionally, ascertainment of 23,649 non-synonymous SNPs and validation of candidate SNPs present in quality related genes suggests their potential utility in genome-wide mapping and marker development for expediting breeding of elite compound-rich tea cultivars.

Full-Length Transcriptome from Camellia oleifera Seed Provides Insight into the Transcript Variants Involved in Oil Biosynthesis

Camellia oleifera Abel., belonging to the genus Camellia of Theaceae, has been widely used as a cooking oil, lubricant, and in cosmetics. Because of complicated polyploidization and large genomes, reference genome information is still lacking. Systematic characterization of gene models based on transcriptome data is a fast and economical approach for C. oleifera. Pacific Biosciences single-molecule long-read isoform sequencing (Iso-Seq) and Illumina RNA-Seq combined with gas chromatography were performed for exploration of oil biosynthesis, accumulation, and comprehensive transcriptome analysis in C. oleifera seeds at five different developmental stages. We report the first full-length transcriptome data set of C. oleifera seeds comprising 40,143 deredundant high-quality isoforms. Among these isoforms, 37,982 were functionally annotated, and 271 (2.43%) belonged to fatty acid metabolism. A total of 8,344 full-length unique transcript models were obtained, and 8,151 (97.69%) of them produced more than two isoforms, suggesting a high degree of transcriptome complexity in C. oleifera seeds. A total of 783 alternative splicing (AS) events were identified, among which the retained intron was the most abundant. We also obtained 1,910 long noncoding RNAs (lncRNAs) and found that AS events occurred in these lncRNAs. Potential transcript variants of genes involved in oil biosynthesis were also investigated. After performing weighted correlation network analysis, we found seven "gene modules" and hub genes for each module showing a significant association with oil content. The series test of clusters classified these modules into four significant profiles based on gene expression patterns. Protein-protein interaction network analysis showed that upregulated WRI1 interacted with 17 genes encoding the enzymes playing key roles in oil synthesis. MYB and ZIP transcriptional factors also showed significant interactions with key genes involved in oil synthesis. Collectively, our data advance the knowledge of RNA isoform diversity in seeds at different developmental stages and provide a rich resource for functional studies on oil synthesis in C. oleifera.

Camellia Plant Resistance and Susceptibility to Petal Blight Disease Are Defined by the Timing of Defense Responses

The family Sclerotiniaceae includes important phytopathogens, such as Botrytis cinerea and Sclerotinia sclerotiorum, that activate plant immune responses to facilitate infection propagation. The mechanisms of plant resistance to these necrotrophic pathogens are still poorly understood. To discover mechanisms of resistance, we used the Ciborinia camelliae (Sclerotiniaceae)-Camellia spp. pathosystem. This fungus induces rapid infection of the blooms of susceptible cultivar Nicky Crisp (Camellia japonica × Camellia pitardii var. pitardii), while Camellia lutchuensis is highly resistant. Genome-wide analysis of gene expression in resistant plants revealed fast modulation of host transcriptional activity 6 h after ascospore inoculation. Ascospores induced the same defense pathways in the susceptible Camellia cultivar but much delayed and coinciding with disease development. We next tested the hypothesis that differences in defense timing influences disease outcome. We induced early defense in the susceptible cultivar using methyl jasmonate and this strongly reduced disease development. Conversely, delaying the response in the resistant species, by infecting it with actively growing fungal mycelium, increased susceptibility. The same plant defense pathways, therefore, contribute to both resistance and susceptibility, suggesting that defense timing is a critical factor in plant health, and resistance against necrotrophic pathogens may occur during the initial biotrophy-like stages.

Complementary transcriptome and proteome profiling in the mature seeds of Camellia oleifera from Hainan Island

Camellia oleifera Abel. (C. oleifera), as an important woody tree species producing edible oils in China, has attracted enormous attention due to its abundant unsaturated fatty acids and their associated benefits to human health. To reveal novel insights into the characters during the maturation period of this plant as well as the molecular basis of fatty acid biosynthesis and degradation, we conducted a conjoint analysis of the transcriptome and proteome of C. oleifera seeds from Hainan Island. Using RNA sequencing (RNA-seq) technology and shotgun proteomic method, 59,391 transcripts and 40,500 unigenes were obtained by TIGR Gene Indices Clustering Tools (TGICL), while 1691 protein species were identified from Mass Spectrometry (MS). Subsequently, all genes and proteins were employed in euKaryotic Orthologous Groups (KOG) classification, Gene Ontology (GO) annotation, and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis to investigate their essential functions. The results indicated that the most abundant pathways were biological metabolic processes. There were 946 unigenes associated with lipid metabolism at the transcriptome level, with 116 proteins at the proteome level; among these, 38 specific proteins were involved in protein-protein interactions, with the majority being related to fatty acid catabolic process. The expression levels of 21 candidate unigenes encoding target proteins were further detected by quantitative real-time polymerase chain reaction (qRT-PCR). Finally, Gas Chromatography Mass Spectrometry (GC-MS) was carried out to determine the fatty acid composition of C. oleifera oil. These findings not only deepened our understanding about the molecular mechanisms of fatty acid metabolism but also offered new evidence concerning the roles of relevant proteins in oil-bearing crops. Furthermore, the lipid-associated proteins recognized in this research might be helpful in providing a reference for the synthetic regulation of C. oleifera oil quality by genetic engineering techniques, thus resulting in potential application in agriculture.

Transcriptomic Analyses of Camellia oleifera 'Huaxin' Leaf Reveal Candidate Genes Related to Long-Term Cold Stress

‘Huaxin’ is a new high-yielding timber cultivar of Camellia oleifera of high economic value, and has been widely cultivated in the red soil hilly region of Hunan Province of the People´s Republic of China in recent years. However, its quality and production are severely affected by low temperatures during flowering. To find genes related to cold tolerance and further explore new candidategenes for chilling-tolerance, Illumina NGS (Next Generation Sequencing) technology was used to perform transcriptomic analyses of C. oleifera ‘Huaxin’ leaves under long-term cold stress. Nine cDNA libraries were sequenced, and 58.31 Gb high-quality clean reads were obtained with an average of 5.92 Gb reads for each sample. A total of 191,150 transcripts were obtained after assembly. Among them, 100,703 unigenes were generated, and 44,610 unigenes were annotated. In total, 1564 differentially expressed genes (DEGs) were identified both in the A_B and A_C gene sets. In the current study, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analyses were performed, andrevealed a group of cold-responsive genes related to hormone regulation, photosynthesis, membrane systems, and osmoregulation; these genes encoded many key proteins in plant biological processes, such as serine/threonine-protein kinase (STPK), transcription factors (TFs), fatty acid desaturase (FAD), lipid-transfer proteins (LTPs), soluble sugars synthetases, and flavonoid biosynthetic enzymes. Some physiological indicators of C. oleifera ‘Huaxin’ were determined under three temperature conditions, and the results were consistent with the molecular sequencing. In addition, the expression levels of 12 DEGs were verified using quantitative real-time polymerase chain reaction (qRT-PCR). In summary, the results of DEGs analysis together with qRT-PCR tests contribute to the understanding of cold tolerance and further exploring new candidate genes for chilling-tolerance in molecular breeding programs of C. oleifera ‘Huaxin’.

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.

Tea Plant Information Archive: a comprehensive genomics and bioinformatics platform for tea plant

Tea is the world's widely consumed nonalcohol beverage with essential economic and health benefits. Confronted with the increasing large-scale omics-data set particularly the genome sequence released in tea plant, the construction of a comprehensive knowledgebase is urgently needed to facilitate the utilization of these data sets towards molecular breeding. We hereby present the first integrative and specially designed web-accessible database, Tea Plant Information Archive (TPIA; http://tpia.teaplant.org). The current release of TPIA employs the comprehensively annotated tea plant genome as framework and incorporates with abundant well-organized transcriptomes, gene expressions (across species, tissues and stresses), orthologs and characteristic metabolites determining tea quality. It also hosts massive transcription factors, polymorphic simple sequence repeats, single nucleotide polymorphisms, correlations, manually curated functional genes and globally collected germplasm information. A variety of versatile analytic tools (e.g. JBrowse, blast, enrichment analysis, etc.) are established helping users to perform further comparative, evolutionary and functional analysis. We show a case application of TPIA that provides novel and interesting insights into the phytochemical content variation of section Thea of genus Camellia under a well-resolved phylogenetic framework. The constructed knowledgebase of tea plant will serve as a central gateway for global tea community to better understand the tea plant biology that largely benefits the whole tea industry.

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

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

Effects of Tree Shape on the Microclimate and Fruit Quality Parameters of Camellia oleifera Abel

Tree shape can affect microclimate, which may directly influence fruit growth and development. Open center and round head are two tree canopy shapes that are generally used in Camellia oleifera Abel. cultivation. To study the effects of the two canopy shapes on the microclimate and fruit quality parameters inside the canopy and investigate the correlations between microclimate and oil yield characteristics, microclimates in different developmental periods and fruit quality parameters at maturity at different positions in the two canopy shapes were determined. In this study, the open-center shape increased the light intensity and temperature, and reduced the relative humidity within the canopy, and the microclimate was more uniformly distributed within the open-center canopy. Fruit quality parameters of open-center canopied trees were also significantly increased and had more uniform distributions, especially the oil yield parameters. The content ratio of the dry seed, oil content ratio of the fresh fruit, and oil yield increased by 43.55%, 59.30% and 79.01%, in the lower layer, respectively, and 27.33%, 41.75%, and 55.43%, in the inner canopy, respectively, compared to those of round-head canopied trees. For both canopy shapes, oil yield parameters had positive correlations with light intensity and temperature, but negative correlations with relative humidity, and the correlation coefficients between light intensity and temperature at different developmental stages and oil yield parameters in round-head canopies were greater than in the open-center canopies. Additionally, regression equations between the single microclimatic factor in different developmental periods and oil yield parameters could be used to estimate the latter in the two canopy shapes. The optimum light conditions for open-center and round-head canopy shapes were 916–1893 and 1108–1259 μmol·m2·s−1, respectively. Overall, the open-center canopy had a more desirable microclimate and a more uniform microclimate distribution, especially the light intensity, which contributes to high uniform fruit quality parameters. Thus, an open-center canopy is the desirable shape in C. oleifera production.

Comparative transcriptomic analysis of high- and low-oil Camellia oleifera reveals a coordinated mechanism for the regulation of upstream and downstream multigenes for high oleic acid accumulation

Tea oil camellia (Camellia oleifera) is an important woody oil tree in southern China. However, little is known regarding the molecular mechanisms that contribute to high oleic acid accumulation in tea oil camellia. Herein, we measured the oil content and fatty acid compositions of high- and low-oil tea oil camellia seeds and investigated the global gene expression profiles by RNA-seq. The results showed that at the early, second and third seed developmental stages, a total of 64, 253, and 124 genes, respectively, were significantly differentially expressed between the high- and low-oil cultivars. Gene ontology (GO) enrichment analysis of the identified differentially expressed transcription factors (TFs; ABI3, FUS3, LEC1, WRI1, TTG2 and DOF4.6) revealed some critical GO terms associated with oil biosynthesis and fatty acid accumulation, including glycolysis, zinc ion binding, positive regulation of fatty acid biosynthetic process, triglyceride biosynthetic process, seed coat development, abscisic acid-mediated signaling pathway and embryo development. Comprehensive comparisons of transcriptomic profiles and expression analysis of multigenes based on qRT-PCR showed that coordinated high expression of the upstream genes HAD, EAR and KASI directly increased the relative levels of C16:0-ACP, which provided enough precursor resources for oleic acid biosynthesis. Continuous high expression of the SAD gene accelerated oleic acid synthesis and accumulation, and coordinated low expression of the downstream genes FAD2, FAD3, FAD7, FAD8 and FAE1 decreased the consumption of oleic acid for conversion. The coordinated regulation of these multigenes ensures the high accumulation of oleic acid in the seeds of tea oil camellia. Our data represent a comprehensive transcriptomic study of high- and low-oil tea oil camellia, not only increasing the number of sequences associated with lipid biosynthesis and fatty acid accumulation in public resource databases but also providing a scientific basis for genetic improvement of the oleic acid content in woody oil trees.

Genome-wide sequencing and metabolic annotation of Pythium irregulare CBS 494.86: understanding Eicosapentaenoic acid production

BACKGROUND

Pythium irregulare is an oleaginous Oomycete able to accumulate large amounts of lipids, including Eicosapentaenoic acid (EPA). EPA is an important and expensive dietary supplement with a promising and very competitive market, which is dependent on fish-oil extraction. This has prompted several research groups to study biotechnological routes to obtain specific fatty acids rather than a mixture of various lipids. Moreover, microorganisms can use low cost carbon sources for lipid production, thus reducing production costs. Previous studies have highlighted the production of EPA by P. irregulare, exploiting diverse low cost carbon sources that are produced in large amounts, such as vinasse, glycerol, and food wastewater. However, there is still a lack of knowledge about its biosynthetic pathways, because no functional annotation of any Pythium sp. exists yet. The goal of this work was to identify key genes and pathways related to EPA biosynthesis, in P. irregulare CBS 494.86, by sequencing and performing an unprecedented annotation of its genome, considering the possibility of using wastewater as a carbon source.

RESULTS

Genome sequencing provided 17,727 candidate genes, with 3809 of them associated with enzyme code and 945 with membrane transporter proteins. The functional annotation was compared with curated information of oleaginous organisms, understanding amino acids and fatty acids production, and consumption of carbon and nitrogen sources, present in the wastewater. The main features include the presence of genes related to the consumption of several sugars and candidate genes of unsaturated fatty acids production.

CONCLUSIONS

The whole metabolic genome presented, which is an unprecedented reconstruction of P. irregulare CBS 494.86, shows its potential to produce value-added products, in special EPA, for food and pharmaceutical industries, moreover it infers metabolic capabilities of the microorganism by incorporating information obtained from literature and genomic data, supplying information of great importance to future work.

RNA-seq data reveals a coordinated regulation mechanism of multigenes involved in the high accumulation of palmitoleic acid and oil in sea buckthorn berry pulp

BACKGROUND

Sea buckthorn is a woody oil crop in which palmitoleic acid (C16:1n7, an omega-7 fatty acid (FA)) contributes approximately 40% of the total FA content in berry pulp (non-seed tissue). However, the molecular mechanisms contributing to the high accumulation of C16:1n7 in developing sea buckthorn berry pulp (SBP) remain poorly understood.

RESULTS

We identified 1737 unigenes associated with lipid metabolism through RNA-sequencing analysis of the four developmental stages of berry pulp in two sea buckthorn lines, 'Za56' and 'TF2-36'; 139 differentially expressed genes were detected between the different berry pulp developmental stages in the two lines. Analyses of the FA composition showed that the C16:1n7 contents were significantly higher in line 'Za56' than in line 'TF2-36' in the mid-late developmental stages of SBP. Additionally, qRT-PCR analyses of 15 genes involved in FA and triacylglycerol (TAG) biosynthesis in both lines revealed that delta9-ACP-desaturase (ACP-Δ9D) competed with 3-ketoacyl-ACP-synthase II (KASII) for the substrate C16:0-ACP and that ACP-Δ9D and delta9-CoA-desaturase (CoA-Δ9D) gene expression positively correlated with C16:1n7 content; KASII and fatty acid elongation 1 (FAE1) gene expression positively correlated with C18:0 content in developing SBP. Specifically, the abundance of ACP-Δ9D and CoA-Δ9D transcripts in line 'Za56', which had a higher C16:1n7 content than line 'TF2-36', suggests that these two genes play an important role in C16:1n7 biosynthesis. Furthermore, the high expressions of the glycerol-3-phosphate dehydrogenase (GPD1) gene and the WRINKLED1 (WRI1) transcription factor contributed to increased biosynthesis of TAG precursor and FAs, respectively, in the early developmental stages of SBP, and the high expression of the diacylglycerol O-acyltransferase 1 (DGAT1) gene increased TAG assembly in the later developmental stages of SBP. Overall, we concluded that increased ACP-Δ9D and CoA-Δ9D levels coupled with decreased KASII and FAE1 activity is a critical event for high C16:1n7 accumulation and that the coordinated high expression of WRI1, GPD1, and DGAT1 genes resulted in high oil accumulation in SBP.

CONCLUSION

Our results provide a scientific basis for understanding the mechanism of high C16:1n7 accumulation in berry pulp (non-seed tissue) and are valuable to the genetic breeding programme for achieving a high quality and yield of SBP oil.

Tandem Mass Tag Based Quantitative Proteomics of Developing Sea Buckthorn Berries Reveals Candidate Proteins Related to Lipid Metabolism

Sea buckthorn ( Hippophae L.) is an economically important shrub or small tree distributed in Eurasia. Most of its well-recognized medicinal and nutraceutical products are derived from its berry oil, which is rich in monounsaturated omega-7 (C16:1) fatty acid and polyunsaturated omega-6 (C18:2) and omega-3 (C18:3) fatty acids. In this study, tandem mass tags (TMT)-based quantitative analysis was used to investigate protein profiles of lipid metabolism in sea buckthorn berries harvested 30, 50, and 70 days after flowering. In total, 8626 proteins were identified, 6170 of which were quantified. Deep analysis results for the proteins identified and related pathways revealed initial fatty acid accumulation during whole-berry development. The abundance of most key enzymes involved in fatty acid and triacylglycerol (TAG) biosynthesis peaked at 50 days after flowering, but TAG synthesis through the PDAT (phospholipid: diacylglycerol acyltransferase) pathway mostly occurred early in berry development. In addition, the patterns of proteins involved in lipid metabolism were confirmed by combined quantitative real-time polymerase chain reaction, enzyme-linked immunosorbent assay, and parallel reaction monitoring analyses. Our data on the proteomic spectrum of sea buckthorn berries provide a scientific basic for understanding lipid metabolism and related pathways in the developing berries.

Association Genetics Identifies Single Nucleotide Polymorphisms Related to Kernel Oil Content and Quality in Camellia oleifera

Camellia oleifera, as an important nonwood tree species for seed oil in China, has received enormous attention owing to its high unsaturated fatty acid contents benefited to human health. It is necessary to examine allelic diversity of key genes that are associated with oil production in C. oleifera cultivars with a large variation of fatty acid compositions. In this study, we performed the association analysis between four key genes (two CoSAD and two Cofad2) coding fatty acid desaturases and traits including oil content and fatty acid composition. We identified two single nucleotide insertion-deletion (InDel) and 362 single-nucleotide polymorphisms (SNPs) within the four candidate genes by sequencing an association population (216 accessions). Single-marker (or haplotype) and traits association tests were conducted by linkage disequilibrium (LD) approaches to detect significant marker-trait associations. Validation population (279 hybrid individuals from six full-sibs families) studies were performed to validate the function of allelic variations significantly associated. In all, 90 single marker-trait and one haplotype-trait associations were significant in association population, and these loci explained 1.87-17.93% proportion of the corresponding phenotypic variance. Further, six SNP marker-trait associations ( Q < 0.10) from Cofad2-A, CoSAD1, and CoSAD2 were successfully validated in the validation population. The SNP markers identified in this study can potentially be applied for future marker-assisted selection to improve oil content and quality in C. oleifera.

Small RNA sequencing reveals dynamic microRNA expression of important nutrient metabolism during development of Camellia oleifera fruit

To obtain insight into the function of miRNAs in the synthesis and storage of important nutrients during the development of Camellia oleifera fruit, Illumina sequencing of flower and fruit small-RNA was conducted. The results revealed that 797 miRNAs were significantly differentially expressed between flower and fruit samples of Camellia oleifera. Through integrated GO and KEGG function annotations, it was determined that the miRNA target genes were mainly involved in metabolic pathways, plant hormone signal transduction, fruit development, mitosis and regulation of biosynthetic processes. Carbohydrate accumulation genes were differentially regulated by miR156, miR390 and miR395 in the fruit growth and development process. MiR477 is the key miRNA functioning in regulation of genes and involved in fatty acid synthesis. Additionally, miR156 also has the function of regulating glycolysis and nutrient transformation genes.

Transcriptome analysis reveals a composite molecular map linked to unique seed oil profile of Neocinnamomum caudatum (Nees) Merr

BACKGROUND

Neocinnamomum caudatum (Nees) Merr., a biodiesel tree species in the subtropical areas of South China, India and Burma, is distinctive from other species in Lauraceae family and its seed oil is rich in linoleic acid (18:2) and stearic acid (18:0). However, there is little genetic information about this species so far. In this study, a transcriptomic analysis on developing seeds of N. caudatum was conducted in an attempt to discern the molecular mechanisms involving the control of the fatty acid (FA) and triacylglycerol (TAG) biosynthesis.

RESULTS

Transcriptome analysis revealed 239,703 unigenes with an average length of 436 bp and 137 putative biomarkers that are related to FA formation and TAG biosynthesis. The expression patterns of genes encoding β-ketoacyl-acyl carrier protein synthase I (KASI), β- ketoacyl-acyl carrier protein synthase II (KASII), stearoyl-ACP desaturase (SAD), fatty acid desaturase 2 (FAD2), fatty acid desaturase 8 (FAD8) and acyl-ACP thioesterase A/B (FATA/B) were further validated by qRT-PCR. These genes displayed a very similar expression pattern in two distinct assays. Moreover, sequence analysis of different FATBs from diverse plant species revealed that NcFATB is structurally different from its counterpart in other species in producing medium-chain saturated FAs. Concertedly, heterologous expression of NcFATB in E. coli BL21 (DE3) strain showed that this corresponding expressed protein, NcFATB, prefers long-chain saturated fatty acyl-ACP over medium-chain fatty acyl-ACP as substrate.

CONCLUSIONS

Transcriptome analysis of developing N. caudatum seeds revealed a composite molecular map linked to the FA formation and oil biosynthesis in this biodiesel tree species. The substrate preference of NcFATB for long-chain saturated FAs is likely to contribute to its unique seed oil profile rich in stearic acid. Our findings demonstrate that in the tree species of Lauraceae family, the FATB enzymes producing long-chain FAs are structurally distinct from those producing medium-chain FAs, thereby suggesting that the FATB genes may serve as a biomarker for the classification of tree species of Lauraceae family.

Comparative transcriptome analysis two genotypes of Acer truncatum Bunge seeds reveals candidate genes that influences seed VLCFAs accumulation

The Acer truncatum Bunge is a particular widespread forest tree species in northern China. VLCFAs are important to eukaryotes survival and play diverse roles throughout the development. So far, there are reports that the Acer truncatum seeds fatty acid (FA) rich in VLCFAs, but little is known about the physiological mechanism responsible for the biosynthesis. A total of approximately 37.07 Gbp was generated, it was comprehensive enough to determine the majority of the regulation VLCFAs biosynthesis genes. The 97,053 different unigenes were assembled and identified, and large numbers of EST-SSRs were determined. The expression profiles of crucial genes (KCS, KCR, HCD and ECR) involved in VLCFAs elongation of fatty acids were also studied. To our knowledge, the present study provides the first comprehensive of the transcriptome of Acer truncatum seeds. This transcriptome dataset have been made publicly available NCBI, we believe that it may provide new resource for future high-throughput gene expression of Acer truncatum seeds growth and development and will provide theoretical basic information for improving the yield of VLCFAs, especially nervonic acid.

Transcriptome analysis reveals crucial genes involved in the biosynthesis of nervonic acid in woody Malania oleifera oilseeds

BACKGROUND

Malania oleifera Chun et Lee (Olacaceae), an evergreen broad-leaved woody tree native to southwest China, is an important oilseed tree. Its seed oil has a high level of nervonic acid (cis-tetracos-15-enoic acid, over 60%), which is essential for human health. M. oleifera seed oil is a promising source of nervonic acid, but little is known about the physiological and molecular mechanisms underlying its biosynthesis.

RESULTS

In this study, we recorded oil accumulation at four stages of seed development. Using a high-throughput RNA-sequencing technique, we obtained 55,843 unigenes, of which 29,176 unigenes were functionally annotated. By comparison, 22,833 unigenes had a two-fold or greater expression at the fast oil accumulation stage than at the initial stage. Of these, 198 unigenes were identified as being functionally involved in diverse lipid metabolism processes (including de novo fatty acid synthesis, carbon chain elongation and modification, and triacylglycerol assembly). Key genes (encoding KCS, KCR, HCD and ECR), putatively responsible for nervonic acid biosynthesis, were isolated and their expression profiles during seed development were confirmed by quantitative real-time PCR analysis. Also, we isolated regulatory factors (such as WRI1, ABI3 and FUS3) that are putatively involved in the regulation of oil biosynthesis and seed development.

CONCLUSION

Our results provide novel data on the physiological and molecular mechanisms of nervonic acid biosynthesis and oil accumulation in M. oleifera seeds, and will also serve as a starting point for biotechnological genetic engineering for the production of nervonic acid resources.

Unraveling the Roles of Regulatory Genes during Domestication of Cultivated Camellia: Evidence and Insights from Comparative and Evolutionary Genomics

With the increasing power of DNA sequencing, the genomics-based approach is becoming a promising resolution to dissect the molecular mechanism of domestication of complex traits in trees. Genus Camellia possesses rich resources with a substantial value for producing beverage, ornaments, edible oil and more. Currently, a vast number of genetic and genomic research studies in Camellia plants have emerged and provided an unprecedented opportunity to expedite the molecular breeding program. In this paper, we summarize the recent advances of gene expression and genomic resources in Camellia species and focus on identifying genes related to key economic traits such as flower and fruit development and stress tolerances. We investigate the genetic alterations and genomic impacts under different selection programs in closely related species. We discuss future directions of integrating large-scale population and quantitative genetics and multiple omics to identify key candidates to accelerate the breeding process. We propose that future work of exploiting the genomic data can provide insights related to the targets of domestication during breeding and the evolution of natural trait adaptations in genus Camellia.

Predicting Potential Distribution and Evaluating Suitable Soil Condition of Oil Tea Camellia in China

Oil tea Camellia, as a major cash and oil crop, has a high status in the forestry cultivation systems in China. To meet the current market demand for oil tea Camellia, its potential distribution and suitable soil condition was researched, to instruct its cultivation and popularization. The potential distribution of oil tea Camellia in China was predicted by the maximum entropy model, using global environmental and soil databases. Then, we collected 10-year literature data about oil tea Camellia soil and applied multiple imputation and factor modeling for an in-depth analysis of soil suitability for growing of oil tea Camellia. The prediction indicated that oil tea Camellia was mainly distributed in Hunan, Jiangxi, Zhejiang, Hainan, East Hubei, Southwest Anhui and most of Guangdong. Climatic factors were more influential than soil factors. The minimum temperature of the coldest month, mean temperature of the coldest quarter and annual precipitation were the most significant contributors to the habitat suitability distribution. In the cultivated area of oil tea Camellia, soil fertility was poor, organic matter was the most significant factor for the soil conditions. Based on climatic and soil factor analyses, our data suggest there is a great potential to spread the oil tea Camellia cultivation industry.

Comparative transcriptomic analysis identifies genes responsible for fruit count and oil yield in the oil tea plant Camellia chekiangoleosa

Fruit yield is an important trait for horticultural crops. However, the limited fruit yield of Camellia chekiangoleosa, a novel promising oil tree, restricts the production of oil. The breeding improvement is little due to its long generation time and lack of available genomic sequence. We identified distinct fruit count phenotypes, high-yield (HY) and low-yield (LY). To understand the underlying genetic basis, we focused on global gene discovery and expression levels in floral buds, which affect fruit count. A total of 140,299 de novo assembled unigenes were obtained using RNA-seq technology, and more genes were expressed in HY than in LY. 2395 differentially expressed genes (DEGs) were identified and enriched in membrane, energy metabolism, secondary metabolism, fatty acid biosynthesis and metabolism, and 18 other metabolic pathways. Of the DEGs, twelve identified transcription factors, including AP2, mostly involve in inflorescence and flower development and in hormone networks. Key DEGs in fatty acid biosynthesis (Fab) FabB, FabF, FabZ, and AccD were highly expressed in floral buds and associated with high oil yield in fruits. We hypothesized that a potential link exists between fruit count and its oil yield. These findings help to elucidate the molecular cues affecting fruit count and oil yield.

Seed Transcriptomics Analysis in Camellia oleifera Uncovers Genes Associated with Oil Content and Fatty Acid Composition

Camellia oleifera is a major tree species for producing edible oil. Its seed oil is well known for the high level of oleic acids; however, little is known regarding the molecular mechanism of lipid biosynthesis in C. oleifera. Here, we measured the oil contents and fatty acid (FA) compositions at four developmental stages and investigated the global gene expression profiles through transcriptomics sequencing. We identified differentially-expressed genes (DEGs) among the developmental stages and found that the distribution of numbers of DEGs was associated with the accumulation pattern of seed oil. Gene Ontology (GO) enrichment analysis revealed some critical biological processes related to oil accumulation, including lipid metabolism and phosphatidylcholine metabolism. Furthermore, we investigated the expression patterns of lipid biosynthesis genes. We showed that most of the genes were identified with single or multiple copies, and some had correlated profiles along oil accumulation. We proposed that the higher levels of stearoyl-ACP desaturases (SADs) coupled with lower activities of fatty acid desaturase 2 (FAD2) might be responsive to the boost of oleic acid at the late stage of C. oleifera seeds’ development. This work presents a comprehensive transcriptomics study of C. oleifera seeds and uncovers valuable DEGs that are associated with the seed oil accumulation.

De novo transcriptome sequencing of Camellia sasanqua and the analysis of major candidate genes related to floral traits

Camellia sasanqua is one of the most famous horticultural plants in Camellia (Theaceae) due to its aesthetic appeal as landscape plant. Knowledge regarding the genetic basis of flowering time, floral aroma and color in C. sasanqua is limited, but is essential to breed new varieties with desired floral traits. Here, we described the de novo transcriptome of young leaves, flower buds and flowers of C. sasanqua. A total of 60,127 unigenes were functionally annotated based on the sequence similarity. After analysis, we found that two floral integrator genes, SOC1 and AP1, in flowering time pathway showed evidence of gene family expansion. Compared with 1-deoxy-D-xylulose-5-phosphate pathway, some genes in the mevalonate pathway were most highly expressed, suggesting that this might represent the major pathway for terpenoid biosynthesis related to floral aroma in C. sasanqua. In flavonoid biosynthesis pathway, PAL, CHI, DFR and ANS showing significantly higher expression levels in flowers and flower buds might have important role in regulation of floral color. The top five most transcription factors (TFs) families in C. sasanqua transcriptome were MYB, MIKC, C3H, FAR1 and HD-ZIP, many of which have a direct relationship with floral traits. In addition, we also identified 33,540 simple sequence repeats (SSRs) in the C. sasanqua transcriptome. Collectively, the C. sasanqua transcriptome dataset generated from this study along with the SSR markers provide a new resource for the identification of novel regulatory transcripts and will accelerate the genetic improvement of C. sasanqua breeding programs.

Transcriptome analysis of the tea oil camellia (Camellia oleifera) reveals candidate drought stress genes

BACKGROUND

The tea-oil camellia (Camellia oleifera) is the most important oil plant in southern China, and has a strong resistance to drought and barren soil. Understanding the molecular mechanisms of drought tolerance would greatly promote its cultivation and molecular breeding.

RESULTS

In total, we obtained 76,585 unigenes with an average length of 810 bp and an N50 of 1,092 bp. We mapped all the unigenes to the NCBI 'nr' (non-redundant), SwissProt, KEGG, and clusters of orthologous groups (COG) databases, where 52,531 (68.6%) unigenes were functionally annotated. According to the annotation, 46,171 (60.8%) unigenes belong to 338 KEGG pathways. We identified a series of unigenes that are related to the synthesis and regulation of abscisic acid (ABA), the activity of protective enzymes, vitamin B6 metabolism, the metabolism of osmolytes, and pathways related to the biosynthesis of secondary metabolites. After exposed to drought for 12 hours, the number of differentially-expressed genes (DEGs) between treated plants and control plants increased in the G4 cultivar, while there was no significant increase in the drought-tolerant C3 cultivar. DEGs associated with drought stress responsive pathways were identified by KEGG pathway enrichment analysis. Moreover, we found 789 DEGs related to transcription factors. Finally, according to the results of qRT-PCR, the expression levels of the 20 unigenes tested were consistent with the results of next-generation sequencing.

CONCLUSIONS

In the present study, we identified a large set of cDNA unigenes from C. oleifera annotated using public databases. Further studies of DEGs involved in metabolic pathways related to drought stress and transcription will facilitate the discovery of novel genes involved in resistance to drought stress in this commercially important plant.

Leaf transcriptome analysis of a subtropical evergreen broadleaf plant, wild oil-tea camellia (Camellia oleifera), revealing candidate genes for cold acclimation

Background

Cold tolerance is a key determinant of the geographical distribution range of a plant species and crop production. Cold acclimation can enhance freezing-tolerance of plant species through a period of exposure to low nonfreezing temperatures. As a subtropical evergreen broadleaf plant, oil-tea camellia demonstrates a relatively strong tolerance to freezing temperatures. Moreover, wild oil-tea camellia is an essential genetic resource for the breeding of cultivated oil-tea camellia, one of the four major woody oil crops in the world. The aims of our study are to identify variations in transcriptomes of wild oil-tea camellia from different latitudes and elevations, and discover candidate genes for cold acclimation.

Results

Leaf transcriptomes were obtained of wild oil-tea camellia from different elevations in Lu and Jinggang Mountains, China. Huge amounts of simple sequence repeats (SSRs), single-nucleotide polymorphisms (SNPs) and insertion/deletions (InDels) were identified. Based on SNPs, phylogenetic analysis was performed to detect genetic structure. Wild oil-tea camellia samples were genetically differentiated mainly between latitudes (between Lu and Jinggang Mountains) and then among elevations (within Lu or Jinggang Mountain). Gene expression patterns of wild oil-tea camellia samples were compared among different air temperatures, and differentially expressed genes (DEGs) were discovered. When air temperatures were below 10 °C, gene expression patterns changed dramatically and majority of the DEGs were up-regulated at low temperatures. More DEGs concerned with cold acclimation were detected at 2 °C than at 5 °C, and a putative C-repeat binding factor (CBF) gene was significantly up-regulated only at 2 °C, suggesting a stronger cold stress at 2 °C. We developed a new method for identifying significant functional groups of DEGs. Among the DEGs, transmembrane transporter genes were found to be predominant and many of them encoded transmembrane sugar transporters.

Conclusions

Our study provides one of the largest transcriptome dataset in the genus Camellia. Wild oil-tea camellia populations were genetically differentiated between latitudes. It may undergo cold acclimation when air temperatures are below 10 °C. Candidate genes for cold acclimation may be predominantly involved in transmembrane transporter activities.

Electronic supplementary material

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

Transcriptome analysis revealed the dynamic oil accumulation in Symplocos paniculata fruit

BACKGROUND

Symplocos paniculata, asiatic sweetleaf or sapphire berry, is a widespread shrub or small tree from Symplocaceae with high oil content and excellent fatty acid composition in fruit. It has been used as feedstocks for biodiesel and cooking oil production in China. Little transcriptome information is available on the regulatory molecular mechanism of oil accumulation at different fruit development stages.

RESULTS

The transcriptome at four different stages of fruit development (10, 80,140, and 170 days after flowering) of S. paniculata were analyzed. Approximately 28 million high quality clean reads were generated. These reads were trimmed and assembled into 182,904 non-redundant putative transcripts with a mean length of 592.91 bp and N50 length of 785 bp, respectively. Based on the functional annotation through Basic Local Alignment Search Tool (BLAST) with public protein database, the key enzymes involved in lipid metabolism were identified, and a schematic diagram of the pathway and temporal expression patterns of lipid metabolism was established. About 13,939 differentially expressed unigenes (DEGs) were screened out using differentially expressed sequencing (DESeq) method. The transcriptional regulatory patterns of the identified enzymes were highly related to the dynamic oil accumulation along with the fruit development of S. paniculata. In addition, quantitative real-time PCR (qRT-PCR) of six vital genes was significantly correlated with DESeq data.

CONCLUSIONS

The transcriptome sequences obtained and deposited in NCBI would enrich the public database and provide an unprecedented resource for the discovery of the genes associated with lipid metabolism pathway in S. paniculata. Results in this study will lay the foundation for exploring transcriptional regulatory profiles, elucidating molecular regulatory mechanisms, and accelerating genetic engineering process to improve the yield and quality of seed oil of S. paniculata.

Transcriptome Analysis Identifies Candidate Genes Related to Triacylglycerol and Pigment Biosynthesis and Photoperiodic Flowering in the Ornamental and Oil-Producing Plant, Camellia reticulata (Theaceae)

Camellia reticulata, which is native to Southwest China, is famous for its ornamental flowers and high-quality seed oil. However, the lack of genomic information for this species has largely hampered our understanding of its key pathways related to oil production, photoperiodic flowering process and pigment biosynthesis. Here, we first sequenced and characterized the transcriptome of a diploid C. reticulata in an attempt to identify genes potentially involved in triacylglycerol biosynthesis (TAGBS), photoperiodic flowering, flavonoid biosynthesis (FlaBS), carotenoid biosynthesis (CrtBS) pathways. De novo assembly of the transcriptome provided a catalog of 141,460 unigenes with a total length of ~96.1 million nucleotides (Mnt) and an N50 of 1080 nt. Of them, 22,229 unigenes were defined as differentially expressed genes (DEGs) across five sequenced tissues. A large number of annotated genes in C. reticulata were found to have been duplicated, and differential expression patterns of these duplicated genes were commonly observed across tissues, such as the differential expression of SOC1_a, SOC1_b, and SOC1_c in the photoperiodic flowering pathway. Up-regulation of SAD_a and FATA genes and down-regulation of FAD2_a gene in the TAGBS pathway in seeds may be relevant to the ratio of monounsaturated fatty acid (MUFAs) to polyunsaturated fatty acid (PUFAs) in seed oil. MYBF1, a transcription regulator gene of the FlaBS pathway, was found with great sequence variation and alteration of expression patterns, probably resulting in functionally evolutionary differentiation in C. reticulata. MYBA1_a and some anthocyanin-specific biosynthetic genes in the FlaBS pathway were highly expressed in both flower buds and flowers, suggesting important roles of anthocyanin biosynthesis in flower development. Besides, a total of 40,823 expressed sequence tag simple sequence repeats (EST-SSRs) were identified in the C. reticulata transcriptome, providing valuable marker resources for further basic and applied researches on this economically important Camellia plant.

Genome-Wide Identification of Genes Probably Relevant to the Uniqueness of Tea Plant (Camellia sinensis) and Its Cultivars

Tea (Camellia sinensis) is a popular beverage all over the world and a number of studies have focused on the genetic uniqueness of tea and its cultivars. However, molecular mechanisms underlying these phenomena are largely undefined. In this report, based on expression data available from public databases, we performed a series of analyses to identify genes probably relevant to the uniqueness of C. sinensis and two of its cultivars (LJ43 and ZH2). Evolutionary analyses showed that the evolutionary rates of genes involved in the pathways were not significantly different among C. sinensis, C. oleifera, and C. azalea. Interestingly, a number of gene families, including genes involved in the pathways synthesizing iconic secondary metabolites of tea plant, were significantly upregulated, expressed in C. sinensis (LJ43) when compared to C. azalea, and this may partially explain its higher content of flavonoid, theanine, and caffeine. Further investigation showed that nonsynonymous mutations may partially contribute to the differences between the two cultivars of C. sinensis, such as the chlorina and higher contents of amino acids in ZH2. Genes identified as candidates are probably relevant to the uniqueness of C. sinensis and its cultivars should be good candidates for subsequent functional analyses and marker-assisted 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 Sequencing of Lima Bean (Phaseolus lunatus) to Identify Putative Positive Selection in Phaseolus and Legumes

The identification of genes under positive selection is a central goal of evolutionary biology. Many legume species, including Phaseolus vulgaris (common bean) and Phaseolus lunatus (lima bean), have important ecological and economic value. In this study, we sequenced and assembled the transcriptome of one Phaseolus species, lima bean. A comparison with the genomes of six other legume species, including the common bean, Medicago, lotus, soybean, chickpea, and pigeonpea, revealed 15 and 4 orthologous groups with signatures of positive selection among the two Phaseolus species and among the seven legume species, respectively. Characterization of these positively selected genes using Non redundant (nr) annotation, gene ontology (GO) classification, GO term enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analyses revealed that these genes are mostly involved in thylakoids, photosynthesis and metabolism. This study identified genes that may be related to the divergence of the Phaseolus and legume species. These detected genes are particularly good candidates for subsequent functional studies.

De novo transcriptome assembly of the wild relative of tea tree (Camellia taliensis) and comparative analysis with tea transcriptome identified putative genes associated with tea quality and stress response

Background

Camellia taliensis is one of the most important wild relatives of cultivated tea tree, C. sinensis. The species extensively occupies mountainous habitats representing a wide-range abiotic tolerance and biotic resistance and thus harbors valuable gene resources that may greatly benefit genetic improvement of cultivated tea tree. However, owning to a large genome size of ~3 Gb and structurally complex genome, there are fairly limited genetic information and particularly few genomic resources publicly available for this species. To better understand the key pathways determining tea flavor and enhance tea tree breeding programs, we performed a high-throughput transcriptome sequencing for C. taliensis.

Results

In this study, approximate 241.5 million high-quality paired-end reads, accounting for ~24 Gb of sequence data, were generated from tender shoots, young leaves, flower buds and flowers using Illumina HiSeq 2000 platform. De novo assembly with further processing and filtering yielded a set of 67,923 transcripts with an average length of 685 bp and an N50 of 995 bp. Based on sequence similarity searches against public databases, a total of 39,475 transcripts were annotated with gene descriptions, conserved protein domains or gene ontology (GO) terms. Candidate genes for major metabolic pathways involved in tea quality were identified and experimentally validated using RT-qPCR. Further gene expression profiles showed that they are differentially regulated at different developmental stages. To gain insights into the evolution of these genes, we aligned them to the previously cloned orthologous genes in C. sinensis, and found that considerable nucleotide variation within several genes involved in important secondary metabolic biosynthesis pathways, of which flavone synthase II gene (FNSII) is the most variable between these two species. Moreover, comparative analyses revealed that C. taliensis shows a remarkable expansion of LEA genes, compared to C. sinensis, which might contribute to the observed stronger stress resistance of C. taliensis.

Conclusion

We reported the first large-coverage transcriptome datasets for C. taliensis using the next-generation sequencing technology. Such comprehensive EST datasets provide an unprecedented opportunity for identifying genes involved in several major metabolic pathways and will accelerate functional genomic studies and genetic improvement efforts of tea trees in the future.

Electronic supplementary material

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

CandiSSR: An Efficient Pipeline used for Identifying Candidate Polymorphic SSRs Based on Multiple Assembled Sequences

Simple sequence repeats (SSRs), also known as microsatellites, are ubiquitous short tandem duplications commonly found in genomes and/or transcriptomes of diverse organisms. They represent one of the most powerful molecular markers for genetic analysis and breeding programs because of their high mutation rate and neutral evolution. However, traditionally experimental screening of the SSR polymorphic status and their subsequent applicability to genetic studies are extremely labor-intensive and time-consuming. Thankfully, the recently decreased costs of next generation sequencing and increasing availability of large genome and/or transcriptome sequences have provided an excellent opportunity and sources for large-scale mining this type of molecular markers. However, current tools are limited. Thus we here developed a new pipeline, CandiSSR, to identify candidate polymorphic SSRs (PolySSRs) based on the multiple assembled sequences. The pipeline allows users to identify putative PolySSRs not only from the transcriptome datasets but also from multiple assembled genome sequences. In addition, two confidence metrics including standard deviation and missing rate of the SSR repetitions are provided to systematically assess the feasibility of the detected PolySSRs for subsequent application to genetic characterization. Meanwhile, primer pairs for each identified PolySSR are also automatically designed and further evaluated by the global sequence similarities of the primer-binding region, ensuring the successful rate of the marker development. Screening rice genomes with CandiSSR and subsequent experimental validation showed an accuracy rate of over 90%. Besides, the application of CandiSSR has successfully identified a large number of PolySSRs in the Arabidopsis genomes and Camellia transcriptomes. CandiSSR and the PolySSR marker sources are publicly available at: http://www.plantkingdomgdb.com/CandiSSR/index.html.