Integrative metabolome and transcriptome analyses reveal the coloration mechanism in Camellia oleifera petals with different color.

BACKGROUND: Camellia olelfera petals are colorful, and have high ornamental value. However, the color formation mechanism of C. olelfera petals with different color is still unclear. In our study, WGCNA method was applied to integrate metabolites and transcriptomes to investigate the coloration mechanism of four C. olelfera cultivars with different petal colors. RESULTS: Here, a total of 372 flavonoids were identified (including 27 anthocyanins), and 13 anthocyanins were significantly differentially accumulated in C. olelfera petals. Among them, cyanidin-3-O-(6''-O-p-Coumaroyl) glucoside was the main color constituent in pink petals, cyanidin-3-O-glucoside, cyanidin-3-O-galactoside, cyanidin-3-O-rutinoside, and cyanidin-3-O-(6''-O-malonyl) glucoside were the main contributors to candy pink petals, and peonidin-3-O-glucoside was the important color substance responsible for the red petals of C. oleifera. Furthermore, six structural genes (Co4CL1, CoF3H1, CoF3'H, CoANS, CoUGT75C1-4, and CoUGT75C1-5), three MYBs (CoMYB1, CoMYB4, and CoMYB44-3), three bHLHs (CobHLH30, CobHLH 77, and CobHLH 79-1), and two WRKYs (CoWRKY7 and CoWRKY22) could be identified candidate genes related to anthocyanins biosynthesis and accumulation, and lead to the pink and red phenotypes. The regulatory network of differentially accumulated anthocyanins and the anthocyanins related genes in C. olelfera petals were established. CONCLUSIONS: These findings elucidate the molecular basis of the coloration mechanisms of pink and red color in C. olelfera petals, and provided valuable target genes for future improvement of petals color in C. olelfera.

Coloration differences in three Camellia reticulata Lindl. cultivars: 'Tongzimian', 'Shizitou' and 'Damanao'.

Camellia reticulata Lindl., also known as Yunnan Camellia, is an important ornamental plant in China, especially for its large and stunning flowers. A comprehensive understanding of their coloration mechanisms can aid breeders in developing new cultivars and improving their ornamental value; however, it is still unclear in Yunnan Camellia, especially in mixed-color flowers. In this study, we conducted metabolic and transcriptomic comparison analyses to investigate the coloration differences in three Yunnan Camellia cultivars: C. reticulata 'Shizitou' (SZT), C. reticulata 'Damanao' (MN), and C. reticulata 'Tongzimian' (TZM). Our results revealed that the initial flowering stage may play a critical role in the color change of MN. Metabolome analysis demonstrated that cyanidin was the primary anthocyanin in SZT and MN's red region, while its content was low in TZM and MN's white region. According to the transcriptome analysis, the anthocyanins biosynthesis pathway was reconstructed in Yunnan Camellia, and the low expression of CHS was detected in TZM and MN's white region, while ANR maintained a high expression level, which may lead to the low content of cyanidin in them. Transcription factors MYBs, bHLH, and bZIP may play a key role in regulating anthocyanin-structural genes. The co-expression analysis showed that the meristem tissue may play a crucial role in the formation of the mixed white-red color in MN. Our study enriched the genetic basis of flower coloration differences in Yunnan Camellia which will be a valuable genomic resource to understanding the biology of coloration formation and for breeding the Camellia cultivars.

Genome-Wide Analysis of MADS-Box Gene Family Reveals CjSTK as a Key Regulator of Seed Abortion in Camellia japonica.

The plant MADS-box transcription factor family is a major regulator of plant flower development and reproduction, and the AGAMOUS-LIKE11/SEEDSTICK (AGL11/STK) subfamily plays conserved functions in the seed development of flowering plants. Camellia japonica is a world-famous ornamental flower, and its seed kernels are rich in highly valuable fatty acids. Seed abortion has been found to be common in C. japonica, but little is known about how it is regulated during seed development. In this study, we performed a genome-wide analysis of the MADS-box gene the in C. japonica genome and identified 126 MADS-box genes. Through gene expression profiling in various tissue types, we revealed the C/D-class MADS-box genes were preferentially expressed in seed-related tissues. We identified the AGL11/STK-like gene, CjSTK, and showed that it contained a typical STK motif and exclusively expressed during seed development. We found a significant increase in the CjSTK expression level in aborted seeds compared with normally developing seeds. Furthermore, overexpression of CjSTK in Arabidopsis thaliana caused shorter pods and smaller seeds. Taken together, we concluded that the fine regulation of the CjSTK expression at different stages of seed development is critical for ovule formation and seed abortion in C. japonica. The present study provides evidence revealing the regulation of seed development in 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.

Chromosome-level genome of Camellia lanceoleosa provides a valuable resource for understanding genome evolution and self-incompatibility.

The section Oleifera (Theaceae) has attracted attention for the high levels of unsaturated fatty acids found in its seeds. Here, we report the chromosome-scale genome of the sect. Oleifera using diploid wild Camellia lanceoleosa with a final size of 3.00 Gb and an N50 scaffold size of 186.43 Mb. Repetitive sequences accounted for 80.63% and were distributed unevenly across the genome. Camellia lanceoleosa underwent a whole-genome duplication event approximately 65 million years ago (65 Mya), prior to the divergence of C. lanceoleosa and Camellia sinensis (approx. 6-7 Mya). Syntenic comparisons of these two species elucidated the genomic rearrangement, appearing to be driven in part by the activity of transposable elements. The expanded and positively selected genes in C. lanceoleosa were significantly enriched in oil biosynthesis, and the expansion of homomeric acetyl-coenzyme A carboxylase (ACCase) genes and the seed-biased expression of genes encoding heteromeric ACCase, diacylglycerol acyltransferase, glyceraldehyde-3-phosphate dehydrogenase and stearoyl-ACP desaturase could be of primary importance for the high oil and oleic acid content found in C. lanceoleosa. Theanine and catechins were present in the leaves of C. lanceoleosa. However, caffeine can not be dectected in the leaves but was abundant in the seeds and roots. The functional and transcriptional divergence of genes encoding SAM-dependent N-methyltransferases may be associated with caffeine accumulation and distribution. Gene expression profiles, structural composition and chromosomal location suggest that the late-acting self-incompatibility of C. lanceoleosa is likely to have favoured a novel mechanism co-occurring with gametophytic self-incompatibility. This study provides valuable resources for quantitative and qualitative improvements and genome assembly of polyploid plants in sect. Oleifera.

Novel insight into anthocyanin metabolism and molecular characterization of its key regulators in Camellia sasanqua.

Flower color is a trait that affects the ornamental value of a plant. Camellia sasanqua is a horticultural plant with rich flower color, but little is known about the regulatory mechanism of color diversity in this plant. Here, the anthocyanin profile of 20 C. sasanqua cultivars revealed and quantified 11 anthocyanin derivatives (five delphinidin-based and six cyanidin-based anthocyanins) for the first time. Cyanidin-3-O-(6-O-(E)-p-coumaroyl)-glucoside was the main contributor to flower base color, and the accumulation of cyanidin and delphinidin derivatives differed in the petals. To further explore the molecular mechanism of color divergence, a transcriptome analysis was performed using C. sasanqua cultivars 'YingYueYe', 'WanXia', 'XueYueHua', and'XiaoMeiGui'. The co-expression network related to differences in delphinidin and cyanidin derivatives accumulation was identified. Eleven candidate genes encoding key enzymes (e.g., F3H, F3'H, and ANS) were involved in anthocyanin biosynthesis. Moreover, 27 transcription factors were screened as regulators of the two types of accumulating anthocyanins. The association was suggested by correlation analysis between the expression levels of the candidate genes and the different camellia cultivars. We concluded that cyanidin and delphinidin derivatives are the major drivers of color diversity in C. sasanqua. This finding provides valuable resources for the study of flower color in C. sasanqua and lays a foundation for genetic modification of anthocyanin biosynthesis.

Ectopic expression of Camellia oleifera Abel. gibberellin 20-oxidase gene increased plant height and promoted secondary cell walls deposition in Arabidopsis.

MAIN CONCLUSION: Ectopic expression of Camellia oleifera Abel. gibberellin 20-oxidase 1 caused a taller phenotype, promoted secondary cell wall deposition, leaf enlargement, and early flowering, and reduced chlorophyll and anthocyanin accumulation and seed enlargement phenotype in Arabidopsis. Plant height and secondary cell wall (SCW) deposition are important plant traits. Gibberellins (GAs) play important roles in regulating plant height and SCWs deposition. Gibberellin 20-oxidase (GA20ox) is an important enzyme involved in GA biosynthesis. In the present study, we identified a GA synthesis gene in Camellia oleifera. The total length of the CoGA20ox1 gene sequence was 1146 bp, encoding 381 amino acids. Transgenic plants with CoGA20ox1 had a taller phenotype; a seed enlargement phenotype; promoted SCWs deposition, leaf enlargement, and early flowering; and reduced chlorophyll and anthocyanin accumulation. Genetic analysis showed that the mutant ga20ox1-3 Arabidopsis partially rescued the phenotype of CoGA20ox1 overexpression plants. The results showed that CoGA20ox1 participates in the growth and development of C. oleifera. The morphological changes in CoGA20ox1 overexpressed plants provide a theoretical basis for further exploration of GA biosynthesis and analysis of the molecular mechanism in C. oleifera.

CamelliA-based simultaneous imaging of Ca2+ dynamics in subcellular compartments.

As a universal second messenger, calcium (Ca2+) transmits specific cellular signals via a spatiotemporal signature generated from its extracellular source and internal stores. Our knowledge of the mechanisms underlying the generation of a Ca2+ signature is hampered by limited tools for simultaneously monitoring dynamic Ca2+ levels in multiple subcellular compartments. To overcome the limitation and to further improve spatiotemporal resolutions, we have assembled a molecular toolset (CamelliA lines) in Arabidopsis (Arabidopsis thaliana) that enables simultaneous and high-resolution monitoring of Ca2+ dynamics in multiple subcellular compartments through imaging different single-colored genetically encoded calcium indicators. We uncovered several Ca2+ signatures in three types of Arabidopsis cells in response to internal and external cues, including rapid oscillations of cytosolic Ca2+ and apical plasma membrane Ca2+ influx in fast-growing Arabidopsis pollen tubes, the spatiotemporal relationship of Ca2+ dynamics in four subcellular compartments of root epidermal cells challenged with salt, and a shockwave-like Ca2+ wave propagating in laser-wounded leaf epidermis. These observations serve as a testimony to the wide applicability of the CamelliA lines for elucidating the subcellular sources contributing to the Ca2+ signatures in plants.

Small RNA profiling reveals that an ovule-specific microRNA, cja-miR5179, targets a B-class MADS-box gene in Camellia japonica.

BACKGROUND AND AIMS: The functional specialization of microRNA and its target genes is often an important factor in the establishment of spatiotemporal patterns of gene expression that are essential to plant development and growth. In different plant lineages, understanding the functional conservation and divergence of microRNAs remains to be explored. METHODS: To identify small regulatory RNAs underlying floral patterning, we performed a tissue-specific profiling of small RNAs in various floral organs from single and double flower varieties (flowers characterized by multiple layers of petals) in Camellia japonica. We identified cja-miR5179, which belongs to a deeply conserved microRNA family that is conserved between angiosperms and basal plants but frequently lost in eudicots. We characterized the molecular function of cja-miR5179 and its target - a B-function MADS-box gene - through gene expression analysis and transient expression assays. KEY RESULTS: We showed that cja-miR5179 is exclusively expressed in ovule tissues at the early stage of floral development. We found that cja-miR5179 targets the coding sequences of a DEFICIENS-like B-class gene (CjDEF) mRNA, which is located in the K motif of the MADS-box domain; and the target sites of miR5179/MADS-box were consistent in Camellia and orchids. Furthermore, through a petal transient-expression assay, we showed that the BASIC PENTACYSTEINE proteins bind to the GA-rich motifs in the cja-miR5179 promoter region and suppresses its expression. CONCLUSIONS: We propose that the regulation between miR5179 and a B-class MADS-box gene in C. japonica has a deep evolutionary origin before the separation of monocots and dicots. During floral development of C. japonica, cja-miR5179 is specifically expressed in the ovule, which may be required for the inhibition of CjDEF function. This work highlights the evolutionary conservation as well as functional divergence of small RNAs in floral development.

Integrated Analysis of Transcriptome and Metabolome Provides Insight into Camellia oleifera Oil Alleviating Fat Accumulation in High-Fat Caenorhabditis elegans.

Camellia oil (CO) is a high medicinal and nutritional value edible oil. However, its ability to alleviate fat accumulation in high-fat Caenorhabditis elegans has not been well elucidated. Therefore, this study aimed to investigate the effect of CO on fat accumulation in high-fat C. elegans via transcriptome and metabolome analysis. The results showed that CO significantly reduced fat accumulation in high-fat C. elegans by 10.34% (Oil Red O method) and 11.54% (TG content method), respectively. Furthermore, CO primarily altered the transcription levels of genes involved in longevity regulating pathway. Specifically, CO decreased lipid storage in high-fat C. elegans by inhibiting fat synthesis. In addition, CO supplementation modulated the abundance of metabolic biomarkers related to pyrimidine metabolism and riboflavin metabolism. The integrated transcriptome and metabolome analyses indicated that CO supplementation could alleviate fat accumulation in high-fat C. elegans by regulating retinol metabolism, drug metabolism-cytochrome P450, metabolism of xenobiotics by cytochrome P450, ascorbate and aldarate metabolism, and pentose and glucuronate interconversions. Overall, these findings highlight the potential health benefits of CO that could potentially be used as a functional edible oil.

Genome-Wide Detection of SPX Family and Profiling of CoSPX-MFS3 in Regulating Low-Phosphate Stress in Tea-Oil Camellia.

Camellia oleifera a member of the family Theaceae, is a phosphorus (P) tolerator native to southern China. The SPX gene family critically regulates plant growth and development and maintains phosphate (Pi) homeostasis. However, the involvement of SPX genes in Pi signaling in Tea-Oil Camellia remains unknown. In this work, 20 SPX genes were identified and categorized into four subgroups. Conserved domains, motifs, gene structure, chromosomal location and gene duplication events were also investigated in the SPX gene family. Defense and stress responsiveness cis-elements were identified in the SPX gene promoters, which participated in low-Pi stress responses. Based on transcriptome data and qRT-PCR results, nine CoSPX genes had similar expression patterns and eight genes (except CoPHO1H3) were up-regulated at 30 days after exposure to low-Pi stress. CoSPX-MFS3 was selected as a key candidate gene by WGCNA analysis. CoSPX-MFS3 was a tonoplast protein. Overexpression of CoSPX-MFS3 in Arabidopsis promoted the accumulation of total P content and decreased the anthocyanin content. Overexpression of CoSPX-MFS3 could enhance low-Pi tolerance by increased biomass and organic acid contents in transgenic Arabidopsis lines. Furthermore, the expression patterns of seven phosphate starvation genes were higher in transgenic Arabidopsis than those in the wild type. These results highlight novel physiological roles of the SPX family genes in C. oleifera under low-Pi stress, and lays the foundation for a deeper knowledge of the response mechanism of C. oleifera to low-Pi stress.

Acclimation in plants - the Green Hub consortium.

Acclimation is the capacity to adapt to environmental changes within the lifetime of an individual. This ability allows plants to cope with the continuous variation in ambient conditions to which they are exposed as sessile organisms. Because environmental changes and extremes are becoming even more pronounced due to the current period of climate change, enhancing the efficacy of plant acclimation is a promising strategy for mitigating the consequences of global warming on crop yields. At the cellular level, the chloroplast plays a central role in many acclimation responses, acting both as a sensor of environmental change and as a target of cellular acclimation responses. In this Perspective article, we outline the activities of the Green Hub consortium funded by the German Science Foundation. The main aim of this research collaboration is to understand and strategically modify the cellular networks that mediate plant acclimation to adverse environments, employing Arabidopsis, tobacco (Nicotiana tabacum) and Chlamydomonas as model organisms. These efforts will contribute to 'smart breeding' methods designed to create crop plants with improved acclimation properties. To this end, the model oilseed crop Camelina sativa is being used to test modulators of acclimation for their potential to enhance crop yield under adverse environmental conditions. Here we highlight the current state of research on the role of gene expression, metabolism and signalling in acclimation, with a focus on chloroplast-related processes. In addition, further approaches to uncovering acclimation mechanisms derived from systems and computational biology, as well as adaptive laboratory evolution with photosynthetic microbes, are highlighted.

Identification of long non-coding RNAs and microRNAs involved in anther development in the tropical Camellia oleifera.

BACKGROUND: Explored the molecular science of anther development is important for improving productivity and overall yield of crops. Although the role of regulatory RNAs, including long non-coding RNAs (lncRNAs) and microRNAs (miRNAs), in regulating anther development has been established, their identities and functions in Camellia oleifera, an important industrial crop, have yet not been clearly explored. Here, we report the identification and characterization of genes, lncRNAs and miRNAs during three stages of the tropical C. oleifera anther development by single-molecule real-time sequencing, RNA sequencing and small RNA sequencing, respectively. RESULTS: These stages, viz. the pollen mother cells stage, tetrad stage and uninucleate pollen stage, were identified by analyzing paraffin sections of floral buds during rapid expansion periods. A total of 18,393 transcripts, 414 putative lncRNAs and 372 miRNAs were identified, of which 5,324 genes, 115 lncRNAs, and 44 miRNAs were differentially accumulated across three developmental stages. Of these, 44 and 92 genes were predicted be regulated by 37 and 30 differentially accumulated lncRNAs and miRNAs, respectively. Additionally, 42 differentially accumulated lncRNAs were predicted as targets of 27 miRNAs. Gene ontology enrichment indicated that potential target genes of lncRNAs were enriched in photosystem II, regulation of autophagy and carbohydrate phosphatase activity, which are essential for anther development. Functional annotation of genes targeted by miRNAs indicated that they are relevant to transcription and metabolic processes that play important roles in microspore development. An interaction network was built with 2 lncRNAs, 6 miRNAs and 10 mRNAs. Among these, miR396 and miR156 family were up-regulated, while their targets, genes (GROWTH REGULATING FACTORS and SQUAMOSA PROMOTER BINDING PROTEIN-LIKE genes) and lncRNAs, were down-regulated. Further, the trans-regulated targets of these lncRNAs, like wall-associated kinase2 and phosphomannose isomerase1, are involved in pollen wall formation during anther development. CONCLUSIONS: This study unravels lncRNAs, miRNAs and miRNA-lncRNA-mRNA networks involved in development of anthers of the tropical C. oleifera lays a theoretical foundation for further elucidation of regulatory roles of lncRNAs and miRNAs in anther development.

De novo transcriptome assembly of the cotyledon of Camellia oleifera for discovery of genes regulating seed germination.

BACKGROUND: Camellia oleifera (C.oleifera) is one of the most important wood oil species in the world. C.oleifera was propagated by nurse seedling grafting. Since the morphology of rootstocks has a significant impact on grafting efficiency and seedling quality, it is necessary to understand the molecular mechanism of morphogenesis for cultivating high-quality and controllable rootstocks. However, the genomic resource for this species is relatively limited, which hinders us from fully understanding the molecular mechanisms of seed germination in C.oleifera. RESULTS: In this paper, using transcriptome sequencing, we measured the gene expression in the C.oleifera cotyledon in different stages of development and the global gene expression profiles. Approximately 45.4 gigabases (GB) of paired-end clean reads were assembled into 113,582 unigenes with an average length of 396 bp. Six public protein databases annotate 61.5% (68,217) of unigenes. We identified 11,391 differentially expressed genes (DEGs) throughout different stages of germination. Enrichment analysis revealed that DEGs were mainly involved in hormone signal transduction and starch sucrose metabolism pathways. The gravitropism regulator UNE10, the meristem regulators STM, KNAT1, PLT2, and root-specific transcription factor WOX11 all have higher gene expression levels in the CAM2 stage (seed soaking), which indicates that the cotyledon-regulated program for germination had initiated when the seeds were imbibition. Our data showed differentially reprogrammed to multiple hormone-related genes in cotyledons during C.oleifera seed germination. CONCLUSION: Cotyledons play vital roles, both as the main nutrient provider and as one primary instructor for seed germination and seedling growth. Together, our study will significantly enrich the genomic resources of Camellia and help us understand the molecular mechanisms of the development in the seed germination and seedling growth of C.oleifera. It is helpful to culture standard and superior quality rootstock for C.oleifera breeding.

The metabolic origins of non-photorespiratory CO2 release during photosynthesis: a metabolic flux analysis.

Respiration in the light (RL) releases CO2 in photosynthesizing leaves and is a phenomenon that occurs independently from photorespiration. Since RL lowers net carbon fixation, understanding RL could help improve plant carbon-use efficiency and models of crop photosynthesis. Although RL was identified more than 75 years ago, its biochemical mechanisms remain unclear. To identify reactions contributing to RL, we mapped metabolic fluxes in photosynthesizing source leaves of the oilseed crop and model plant camelina (Camelina sativa). We performed a flux analysis using isotopic labeling patterns of central metabolites during 13CO2 labeling time course, gas exchange, and carbohydrate production rate experiments. To quantify the contributions of multiple potential CO2 sources with statistical and biological confidence, we increased the number of metabolites measured and reduced biological and technical heterogeneity by using single mature source leaves and quickly quenching metabolism by directly injecting liquid N2; we then compared the goodness-of-fit between these data and data from models with alternative metabolic network structures and constraints. Our analysis predicted that RL releases 5.2 µmol CO2 g-1 FW h-1 of CO2, which is relatively consistent with a value of 9.3 µmol CO2 g-1 FW h-1 measured by CO2 gas exchange. The results indicated that ≤10% of RL results from TCA cycle reactions, which are widely considered to dominate RL. Further analysis of the results indicated that oxidation of glucose-6-phosphate to pentose phosphate via 6-phosphogluconate (the G6P/OPP shunt) can account for >93% of CO2 released by RL.

cDNA cloning, prokaryotic expression, and functional analysis of squalene synthase (SQS) in Camellia vietnamensis Huang.

Camellia vietnamensis Huang, which belongs to Camellia oleifera, is a traditional Chinese medicinal plant widely planted on Hainan Island. Tea saponin is an important functional component of C. vietnamensis, and squalene is the precursor substance that controls its formation. Squalene synthase (SQS: EC 2.5.1.21) synthesizes squalene from 2 molecules of farnesyl pyrophosphate (FPP). In this study, 1683 bp of the C. vietnamensis SQS gene, designated as CvSQS, was cloned and encoded 414 amino acids. Bioinformatics and phylogenetic tree analysis revealed the high homology of CvSQS with squalene synthases from other plants. For soluble proteins, the carboxy-terminal deleted CvSQS was obtained for expression in Escherichia coli Transetta (DE3), and the recombinant protein with a weight of 42.5 kDa was detected using SDS-PAGE and Western blot. After an enzymatic reaction, the presence of squalene in the product was analyzed using GC-MS detection, which indicated that CvSQS had catalytic activity. The tissue specificity of CvSQS and its presence in seeds at various ripening stages was detected by q-RT PCR. CvSQS had the highest transcriptional level in leaves, followed by seeds, roots, and flowers; the amount of CvSQS in the seeds was highest in September. The identification and functional characterization of CvSQS is essential for further studies on the regulation mechanism of tea saponin in C. vietnamensis.

Whole-Genome Identification and Analysis of Multiple Gene Families Reveal Candidate Genes for Theasaponin Biosynthesis in Camellia oleifera.

Camellia oleifera is an economically important oilseed tree. Seed meals of C. oleifera have a long history of use as biocontrol agents in shrimp farming and as cleaning agents in peoples' daily lives due to the presence of theasaponins, the triterpene saponins from the genus Camellia. To characterize the biosynthetic pathway of theasaponins in C. oleifera, members of gene families involved in triterpenoid biosynthetic pathways were identified and subjected to phylogenetic analysis with corresponding members in Arabidopsis thaliana, Camellia sinensis, Actinidia chinensis, Panax ginseng, and Medicago truncatula. In total, 143 triterpenoid backbone biosynthetic genes, 1169 CYP450s, and 1019 UGTs were identified in C. oleifera. The expression profiles of triterpenoid backbone biosynthetic genes were analyzed in different tissue and seed developmental stages of C. oleifera. The results suggested that MVA is the main pathway for triterpenoid backbone biosynthesis. Moreover, the candidate genes for theasaponin biosynthesis were identified by WGCNA and qRT-PCR analysis; these included 11 CYP450s, 14 UGTs, and eight transcription factors. Our results provide valuable information for further research investigating the biosynthetic and regulatory network of theasaponins.

Integrative Physiological and Transcriptomic Analysis Reveals the Transition Mechanism of Sugar Phloem Unloading Route in Camellia oleifera Fruit.

Sucrose phloem unloading plays a vital role in photoassimilate distribution and storage in sink organs such as fruits and seeds. In most plants, the phloem unloading route was reported to shift between an apoplasmic and a symplasmic pattern with fruit development. However, the molecular transition mechanisms of the phloem unloading pathway still remain largely unknown. In this study, we applied RNA sequencing to profile the specific gene expression patterns for sucrose unloading in C. oleifera fruits in the apo- and symplasmic pathways that were discerned by CF fluoresce labelling. Several key structural genes were identified that participate in phloem unloading, such as PDBG11, PDBG14, SUT8, CWIN4, and CALS10. In particular, the key genes controlling the process were involved in callose metabolism, which was confirmed by callose staining. Based on the co-expression network analysis with key structural genes, a number of transcription factors belonging to the MYB, C2C2, NAC, WRKY, and AP2/ERF families were identified to be candidate regulators for the operation and transition of phloem unloading. KEGG enrichment analysis showed that some important metabolism pathways such as plant hormone metabolism, starch, and sucrose metabolism altered with the change of the sugar unloading pattern. Our study provides innovative insights into the different mechanisms responsible for apo- and symplasmic phloem unloading in oil tea fruit and represents an important step towards the omics delineation of sucrose phloem unloading transition in crops.

Integrated Transcriptome and Metabolome Analysis Reveals Key Metabolites Involved in Camellia oleifera Defense against Anthracnose.

Camellia oleifera (Ca. oleifera) is a woody tree species cultivated for the production of edible oil from its seed. The growth and yield of tea-oil trees are severely affected by anthracnose (caused by Colletotrichum gloeosporioides). In this study, the transcriptomic and metabolomic analyses were performed to detect the key transcripts and metabolites associated with differences in the susceptibility between anthracnose-resistant (ChangLin150) and susceptible (ChangLin102) varieties of Ca. oleifera. In total, 5001 differentially expressed genes (DEGs) were obtained, of which 479 DEGs were common between the susceptible and resistant varieties and further analyzed. KEGG enrichment analysis showed that these DEGs were significantly enriched in tyrosine metabolism, phenylpropanoid biosynthesis, flavonoid biosynthesis and isoquinoline alkaloid biosynthesis pathways. Furthermore, 68 differentially accumulated metabolites (DAMs) were detected, including flavonoids, such as epicatechin, phenethyl caffeate and procyanidin B2. Comparison of the DEGs and DAMs revealed that epicatechin, procyanidin B2 and arachidonic acid (peroxide free) are potentially important. The expression patterns of genes involved in flavonoid biosynthesis were confirmed by qRT-PCR. These results suggested that flavonoid biosynthesis might play an important role in the fight against anthracnose. This study provides valuable molecular information about the response of Ca. oleifera to Co. gloeosporioides infection and will aid the selection of resistant varieties using marker-assisted breeding.

Transcriptome Analysis Reveals Putative Induction of Floral Initiation by Old Leaves in Tea-Oil Tree (Camellia oleifera 'changlin53').

Floral initiation is a major phase change in the spermatophyte, where developmental programs switch from vegetative growth to reproductive growth. It is a key phase of flowering in tea-oil trees that can affect flowering time and yield, but very little is known about the molecular mechanism of floral initiation in tea-oil trees. A 12-year-old Camellia oleifera (cultivar 'changlin53') was the source of experimental materials in the current study. Scanning electron microscopy was used to identify the key stage of floral initiation, and transcriptome analysis was used to reveal the transcriptional regulatory network in old leaves involved in floral initiation. We mined 5 DEGs related to energy and 55 DEGs related to plant hormone signal transduction, and we found floral initiation induction required a high level of energy metabolism, and the phytohormones signals in the old leaves regulate floral initiation, which occurred at stage I and II. Twenty-seven rhythm-related DEGs and 107 genes associated with flowering were also identified, and the circadian rhythm interacted with photoperiod pathways to induce floral initiation. Unigene0017292 (PSEUDO-RESPONSE REGULATOR), Unigene0046809 (LATE ELONGATED HYPOCOTYL), Unigene0009932 (GIGANTEA), Unigene0001842 (CONSTANS), and Unigene0084708 (FLOWER LOCUS T) were the key genes in the circadian rhythm-photoperiod regulatory network. In conjunction with morphological observations and transcriptomic analysis, we concluded that the induction of floral initiation by old leaves in C. oleifera 'changlin53' mainly occurred during stages I and II, floral initiation was completed during stage III, and rhythm-photoperiod interactions may be the source of the main signals in floral initiation induced by old leaves.