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

Identification of genes associated with fatty acid biosynthesis based on 214 safflower core germplasm

BACKGROUND

Safflower (Carthamus tinctorius L.) is an oilseed crop with substantial medicinal and economic value. However, the methods for constructing safflower core germplasm resources are limited, and the molecular mechanisms of lipid biosynthesis in safflower seeds are not well understood.

RESULTS

In this study, 11 oil-related quantitative traits and 50 pairs of InDel markers were used to assess the diversity of a collection of 605 safflower germplasms. The original safflower germplasm exhibited rich phenotypic diversity, with high variation for most of the phenotypic traits under investigation. Similarly, high genetic diversity was evaluated in the original germplasm, in which the mean Shannon's information index (I), observed heterozygosity (H0), and expected heterozygosity (He) were 0.553, 0.182, and 0.374, respectively. Four subgroups with strong genetic structures were identified and a core germplasm of 214 cultivars was constructed, which is well represented in the original germplasm. Meanwhile, differential expression analysis of the transcriptomes of high and low linoleic acid safflower varieties at two stages of seed development identified a total of 47 genes associated with lipid biosynthesis. High expression of the genes KAS II and SAD enhanced the synthesis and accumulation of oleic acid, while FAD genes like FAD2 (Chr8G0104100), FAD3, FAD7 and FAD8 promoted the consumption of oleic acid conversion. The coordinated regulation of these multiple genes ensures the high accumulation of oleic acid in safflower seed oil.

CONCLUSIONS

Based on these findings, a core germplasm of 214 cultivars was constructed and 47 candidate genes related to unsaturated fatty acid biosynthesis and lipid accumulation were identified. These results not only provide guidance for further studies to elucidate the molecular basis of oil lipid accumulation in safflower seeds, but also contribute to safflower cultivar improvements.

Oilomics: An important branch of foodomics dealing with oil science and technology

Oil is one of three nutritious elements. The application of omics techniques in the field of oil science and technology is attracted increasing attention. Oilomics, which emerged as an important branch of foodomics, has been widely used in various aspects of oil science and technology. However, there are currently no articles systematically reviewing the application of oilomics. This paper aims to provide a critical overview of the advantages and value of oilomics technology compared to traditional techniques in various aspects of oil science and technology, including oil nutrition, oil processing, oil quality, safety, and traceability. Moreover, this article intends to review major issues in oilomics and give a comprehensive, critical overview of the current state of the art, future challenges and trends in oilomics, with a view to promoting the optimal application and development of oilomics technology in oil science and technology.

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.

Regulation of seed oil accumulation by lncRNAs in Brassica napus

BACKGROUND

Studies have indicated that long non-coding RNAs (lncRNAs) play important regulatory roles in many biological processes. However, the regulation of seed oil biosynthesis by lncRNAs remains largely unknown.

RESULTS

We comprehensively identified and characterized the lncRNAs from seeds in three developing stages in two accessions of Brassica napus (B. napus), ZS11 (high oil content) and WH5557 (low oil content). Finally, 8094 expressed lncRNAs were identified. LncRNAs MSTRG.22563 and MSTRG.86004 were predicted to be related to seed oil accumulation. Experimental results show that the seed oil content is decreased by 3.1-3.9% in MSTRG.22563 overexpression plants, while increased about 2% in MSTRG.86004, compared to WT. Further study showed that most genes related to lipid metabolism had much lower expression, and the content of some metabolites in the processes of respiration and TCA (tricarboxylic acid) cycle was reduced in MSTRG.22563 transgenic seeds. The expression of genes involved in fatty acid synthesis and seed embryonic development (e.g., LEC1) was increased, but genes related to TAG assembly was decreased in MSTRG.86004 transgenic seeds.

CONCLUSION

Our results suggest that MSTRG.22563 might impact seed oil content by affecting the respiration and TCA cycle, while MSTRG.86004 plays a role in prolonging the seed developmental time to increase seed oil accumulation.

Detection and Quantification of Anthracnose Pathogen Colletotrichum fructicola in Cultivated Tea-Oil Camellia Species from Southern China Using a DNA-Based qPCR Assay

Tea-oil Camellia species as edible-oil producing trees are widely cultivated in southern China. Camellia anthracnose that is mainly caused by Colletotrichum fructicola is a major disease of tea-oil trees. However, rapid detection and precise quantification of C. fructicola in different Camellia species that are crucial for the fundamental study of this pathosystem and effective disease management remain largely unexplored. Here, we developed a sensitive, rapid, and accurate method for quantifying C. fructicola growth in different Camellia species using a quantitative PCR assay. Amplified C. fructicola DNA using ITS-specific primers is relatively compared with the amplification of Camellia oleifera using the TUB gene. We determined that the fungal growth is tightly associated with the disease development in Ca. oleifera following C. fructicola infection in a time-course manner. This assay is highly sensitive, as fungal growth was detected in six different inoculated tea-oil Camellia species without visible disease lesion symptoms. Additionally, this method was validated by quantifying the Camellia anthracnose in orchards that did not show any disease symptoms. This assay enables the rapid, highly sensitive, and precise detection and quantification of C. fructicola growth in different tea-oil Camellia species, which will have a practical application for early diagnosis of anthracnose disease under asymptomatic conditions in Camellia breeding and field and will facilitate the development of tea-oil trees and C. fructicola interaction as a mold system to study woody plant and fungal pathogens interaction.

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.

Maturity Grading and Identification of Camellia oleifera Fruit Based on Unsupervised Image Clustering

Maturity grading and identification of Camellia oleifera are prerequisites to determining proper harvest maturity windows and safeguarding the yield and quality of Camellia oil. One problem in Camellia oleifera production and research is the worldwide confusion regarding the grading and identification of Camellia oleifera fruit maturity. To solve this problem, a Camellia oleifera fruit maturity grading and identification model based on the unsupervised image clustering model DeepCluster has been developed in the current study. The proposed model includes the following two branches: a maturity grading branch and a maturity identification branch. The proposed model jointly learns the parameters of the maturity grading branch and maturity identification branch and used the maturity clustering assigned from the maturity grading branch as pseudo-labels to update the parameters of the maturity identification branch. The maturity grading experiment was conducted using a training set consisting of 160 Camellia oleifera fruit samples and 2628 Camellia oleifera fruit digital images collected using a smartphone. The proposed model for grading Camellia oleifera fruit samples and images in training set into the following three maturity levels: unripe (47 samples and 883 images), ripe (62 samples and 1005 images), and overripe (51 samples and 740 images). Results suggest that there was a significant difference among the maturity stages graded by the proposed method with respect to seed oil content, seed soluble protein content, seed soluble sugar content, seed starch content, dry seed weight, and moisture content. The maturity identification experiment was conducted using a testing set consisting of 160 Camellia oleifera fruit digital images (50 unripe, 60 ripe, and 50 overripe) collected using a smartphone. According to the results, the overall accuracy of maturity identification for Camellia oleifera fruit was 91.25%. Moreover, a Gradient-weighted Class Activation Mapping (Grad-CAM) visualization analysis reveals that the peel regions, crack regions, and seed regions were the critical regions for Camellia oleifera fruit maturity identification. Our results corroborate a maturity grading and identification application of unsupervised image clustering techniques and are supported by additional physical and quality properties of maturity. The current findings may facilitate the harvesting process of Camellia oleifera fruits, which is especially critical for the improvement of Camellia oil production and quality.

Maturity Stage Discrimination of Camellia oleifera Fruit Using Visible and Near-Infrared Hyperspectral Imaging

The maturity of Camellia oleifera fruit is one of the most important indicators to optimize the harvest day, which, in turn, results in a high yield and good quality of the produced Camellia oil. A hyperspectral imaging (HSI) system in the range of visible and near-infrared (400–1000 nm) was employed to assess the maturity stages of Camellia oleifera fruit. Hyperspectral images of 1000 samples, which were collected at five different maturity stages, were acquired. The spectrum of each sample was extracted from the identified region of interest (ROI) in each hyperspectral image. Spectral principal component analysis (PCA) revealed that the first three PCs showed potential for discriminating samples at different maturity stages. Two classification models, including partial least-squares discriminant analysis (PLS-DA) and principal component analysis discriminant analysis (PCA-DA), based on the raw or pre-processed full spectra, were developed, and performances were compared. Using a PLS-DA model, based on second-order (2nd) derivative pre-processed spectra, achieved the highest results of correct classification rates (CCRs) of 99.2%, 98.4%, and 97.6% in the calibration, cross-validation, and prediction sets, respectively. Key wavelengths selected by PC loadings, two-dimensional correlation spectroscopy (2D-COS), and the uninformative variable elimination and successive projections algorithm (UVE+SPA) were applied as inputs of the PLS-DA model, while UVE-SPA-PLS-DA built the optimal model with the highest CCR of 81.2% in terms of the prediction set. In a confusion matrix of the optimal simplified model, satisfactory sensitivity, specificity, and precision were acquired. Misclassification was likely to occur between samples at maturity stages two, three, and four. Overall, an HSI with effective selected variables, coupled with PLS-DA, could provide an accurate method and a reference simple system by which to rapidly discriminate the maturity stages of Camellia oleifera fruit samples.

Antifungal action and induction of resistance by Bacillus sp. strain YYC 155 against Colletotrichum fructicola for control of anthracnose disease in Camellia oleifera

Anthracnose disease caused by Colletotrichum fructicola is a serious disease that can afflict Camellia oleifera. Biological control is a rapidly growing approach for the management of plant diseases. In this study, we investigated the bio-control efficiency and the defense responses of an endophytic Bacillus tequilensis strain YYC 155, which was isolated from the root nodules of the Crotalaria pallida against anthracnose disease, caused by C. fructicola in C. oleifera. B. tequilensis YYC 155 exhibited significant inhibitory activity against anthracnose disease, caused by C. fructicola in C. oleifera. YYC 155 can secrete extracellular hydrolases, such as chitinase and β-1, 3-glucanase, which produce lipopeptides that are antimicrobial and forms strong biofilms. In addition, in treatment with YYC 155, the cell membranes of C. fructicola were injured and the leakage of cell contents from the mycelia of the pathogen was increased. Spraying 1 × 10⁷ cells mL^–1 bacterial suspension of YYC 155 on C. oleifera leaves enhanced the activity of key enzymes in C. oleifera associated with the phenylpropanoid pathway and increased the content of phenolic compounds and flavonoids. Results of our study indicate that B. tequilensis YYC 155 may potentially represent an effective biocontrol agent against anthracnose disease in C. oleifera.

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.

Critical metabolic pathways and SAD/FADs, WRI1s, and DGATs cooperate for high-oleic acid oil production in developing oil tea (Camellia oleifera) seeds

Oil tea trees produce high-quality edible oils with desirably high oleic acid (18:1) and low linoleic (18:2) and linolenic (18:3) fatty acid (FA) levels, but limited understanding of tea oil biosynthesis and regulation has become a significant obstacle for the breeding of high-yield and -quality oil tea varieties. By integrating metabolite and transcriptome analyses of developing oil tea seeds, we dissected the critical metabolic pathways, including glycolysis, fatty acid, and triacylglycerol (TAG) biosynthesis, as well as genes essential for tea seed oil production. Two plastidic stearoyl-acyl carrier protein desaturases (CoSAD1 and 2) and two endoplasmic reticulum-localized FA desaturases (CoFAD2 and 3) were functionally characterized as responsible for high 18:1 and low 18:2 and 18:3 proportions in tea oils. Two diacylglycerol O-acyltransferases (CoDGAT1 and 2) that may prefer to synthesize 18:1-TAG were functionally characterized and might be also important for high 18:1-TAG production. The highly expressed CoWRI1a and b were identified and characterized as activators of glycolysis and regulators of directing source carbon flux into FA biosynthesis in developing oil tea seeds. The upregulated CoSADs with downregulated CoFAD2 and CoFAD3 at the late seed developmental stages mainly accounted for high 18:1 levels. Two CoDGATs might be responsible for assembling TAGs with oleoyl acyl chains, whilst two CoWRI1s regulated carbons from parental sources, partitioning into oil production in oil tea embryo sinks. This study provides a deep understanding of the biosynthesis of tea seed oils and information on genes that may be used as molecular markers to breed oil tea varieties with higher oil yield and quality.

Identification of miRNA-mRNA Regulatory Modules Involved in Lipid Metabolism and Seed Development in a Woody Oil Tree (Camellia oleifera)

Tea oil camellia (Camellia oleifera), an important woody oil tree, is a source of seed oil of high nutritional and medicinal value that is widely planted in southern China. However, there is no report on the identification of the miRNAs involved in lipid metabolism and seed development in the high- and low-oil cultivars of tea oil camellia. Thus, we explored the roles of miRNAs in the key periods of oil formation and accumulation in the seeds of tea oil camellia and identified miRNA–mRNA regulatory modules involved in lipid metabolism and seed development. Sixteen small RNA libraries for four development stages of seed oil biosynthesis in high- and low-oil cultivars were constructed. A total of 196 miRNAs, including 156 known miRNAs from 35 families, and 40 novel miRNAs were identified, and 55 significantly differentially expressed miRNAs were found, which included 34 upregulated miRNAs, and 21 downregulated miRNAs. An integrated analysis of the miRNA and mRNA transcriptome sequence data revealed that 10 miRNA–mRNA regulatory modules were related to lipid metabolism; for example, the regulatory modules of ath-miR858b–MYB82/MYB3/MYB44 repressed seed oil biosynthesis, and a regulation module of csi-miR166e-5p–S-ACP-DES6 was involved in the formation and accumulation of oleic acid. A total of 23 miRNA–mRNA regulatory modules were involved in the regulation of the seed size, such as the regulatory module of hpe-miR162a_L-2–ARF19, involved in early seed development. A total of 12 miRNA–mRNA regulatory modules regulating growth and development were identified, such as the regulatory modules of han-miR156a_L+1–SPL4/SBP2, promoting early seed development. The expression changes of six miRNAs and their target genes were validated using quantitative real-time PCR, and the targeting relationship of the cpa-miR393_R-1–AFB2 regulatory module was verified by luciferase assays. These data provide important theoretical values and a scientific basis for the genetic improvement of new cultivars of tea oil camellia in the future.

Comparative study on fruit development and oil synthesis in two cultivars of Camellia oleifera

BACKGROUND

The oil-tea tree (Camellia oleifera Abel.) is a woody tree species that produces edible oil in the seed. C. oleifera oil has high nutritional value and is also an important raw material for medicine and cosmetics. In China, due to the uncertainty on maturity period and oil synthesis mechanism of many C. oleifera cultivars, growers may harvest fruits prematurely, which could not maximize fruit and oil yields. In this study, our objective was to explore the mechanism and differences of oil synthesis between two Camellia oleifera cultivars for a precise definition of the fruit ripening period and the selection of appropriate cultivars.

RESULTS

The results showed that 'Huashuo' had smaller fruits and seeds, lower dry seed weight and lower expression levels of fatty acid biosynthesis genes in July. We could not detect the presence of oil and oil bodies in 'Huashuo' seeds until August, and oil and oil bodies were detected in 'Huajin' seeds in July. Moreover, 'Huashuo' seeds were not completely blackened in October with up to 60.38% of water and approximately 37.98% of oil in seed kernels whose oil content was much lower than normal mature seed kernels. The oil bodies in seed endosperm cells of 'Huajin' were always higher than those of 'Huashuo' from July to October.

CONCLUSION

Our results confirmed that C. oleifera 'Huashuo' fruits matured at a lower rate compared to 'Huajin' fruits and that 'Huajin' seeds entered the oil synthesis period earlier than 'Huashuo' seeds. Moreover, 'Huashuo' fruits did not mature during the Frost's Descent period (October 23-24 each year).

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.

Dynamic transcriptome analysis identifies genes related to fatty acid biosynthesis in the seeds of Prunus pedunculata Pall

BACKGROUND

Prunus pedunculata Pall, the deciduous shrub of Amygdalus subgenus in Rosaceae, is a new kind of desert oil-bearing tree. It has a long story of being planted in the West and North of China for sand fixation and desert control. In addition, the seeds of P. pedunculata are rich of oil, especially the monounsaturated fatty acid and polyunsaturated fatty acid. However, little is known about the molecular mechanisms of oil accumulation during the seed development of P. pedunculata.

RESULTS

The seeds of P. pedunculata from three independent plants at 10, 18, 24, 31, 39, 45, 59 and 73 days after flowering (DAF) were obtained and the oil compositions were evaluated. It showed that oleic acid was the dominant type of oil content in the mature seeds (from 32.724% at 10DAF to 72.06% at 73DAF). Next, transcriptome sequencing for the developing seeds produced 988.795 million high quality reads and TRINITY assembled 326,271 genes for the first transcriptome for P. pedunculata. After the assembled transcriptome was evaluated by BUSCO with 85.9% completeness, we identified 195,342, 109,850 and 121,897 P. pedunculata genes aligned to NR, GO and KEGG pathway databases, respectively. Then, we predicted 23,229 likely proteins from the assembled transcriptome and identified 1917 signal peptides and 5512 transmembrane related proteins. In the developing seeds we detected 91,362 genes (average FPKM > 5) and correlation analysis indicated three possible development stages - early (10 ~ 24DAF), middle (31 ~ 45DAF) and late (59 ~ 73DAF). We next analyzed the differentially expressed genes (DEGs) in the developing seeds. Interestingly, compared to 10DAF the number of DEGs was increased from 4406 in 18DAF to 27,623 in 73DAF. Based on the gene annotation, we identified 753, 33, 8 and 645 DEGs related to the fatty acid biosynthesis, lipid biosynthesis, oil body and transcription factors. Notably, GPAT, DGD1, LACS2, UBC and RINO were highly expressed at the early development stage, ω6-FAD, SAD, ACP, ACCA and AHG1 were highly expressed at the middle development stage, and LACS6, DGD1, ACAT1, AGPAT, WSD1, EGY2 and oleosin genes were highly expressed at the late development stage.

CONCLUSIONS

This is the first time to study the developing seed transcriptome of P. pedunculata and our findings will provide a valuable resource for future studies. More importantly, it will improve our understanding of molecular mechanisms of oil accumulation in P. pedunculata.

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.

Global Transcriptome and Correlation Analysis Reveal Cultivar-Specific Molecular Signatures Associated with Fruit Development and Fatty Acid Determination in Camellia oleifera Abel

Background

Oil-tea Camellia is a very important edible oil plant widely distributed in southern China. Tea oil extracted from the oil-tea Camellia seeds is beneficial to health and is considered as a health edible oil. We attempt to identify genes related to fatty acid biosynthesis in an oil-tea Camellia seed kernel, generated a comprehensive transcriptome analysis of the seed kernel at different developmental stages, and explore optimal picking time of fruit. Material and Methods. A gas chromatography-mass spectrometer was used to detect the content of various fatty acids in samples. Transcriptome analysis was performed to detect gene dynamics and corresponding functions.

Results

Multiple phenotypic data were counted in detail, including the oil content, oleic acid content, linoleic acid content, linolenic acid content, fruit weight, fruit height, fruit diameter, single seed weight, seed length, and seed width in different developmental stages, which indicate that a majority of indicators increased with the development of oil-tea Camellia. The transcriptomics was conducted to perform a comprehensive and system-level view on dynamic gene expression networks for different developmental stages. Short Time-series Expression Miner (STEM) analysis of XL106 (the 6 time points) and XL210 (8 time points) was performed to screen related fatty acid (FA) gene set, from which 1041 candidate genes related to FA were selected in XL106 and 202 related genes were screened in XL210 based on GO and KEGG enrichment. Then, candidate genes and trait dataset were combined to conduct correlation analysis, and 10 genes were found to be strongly connected with several key traits.

Conclusions

The multiple phenotypic data revealed the dynamic law of changes during the picking stage. Transcriptomic analysis identified a large number of potential key regulatory factors that can control the oil content of dried kernels, oleic acid, linoleic acid, linolenic acid, fresh seed rate, and kernel-to-seed ratio, thereby providing a new insight into the molecular networks underlying the picking stage of oil-tea Camellia, which provides a theoretical basis for the optimal fruit picking point.

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.

De novo transcriptome assembly for the five major organs of Zanthoxylum armatum and the identification of genes involved in terpenoid compound and fatty acid metabolism

Background

Zanthoxylum armatum (Z. armatum) is a highly economically important tree that presents a special numbing taste. However, the underlying regulatory mechanism of the numbing taste remains poorly understood. Thus, the elucidation of the key genes associated with numbing taste biosynthesis pathways is critical for providing genetic information on Z. armatumand the breeding of high-quality germplasms of this species.

Results

Here, de novo transcriptome assembly was performed for the five major organs of Z. armatum, including the roots, stems, leaf buds, mature leaves and fruits. A total of 111,318 unigenes were generated with an average length of 1014 bp. Additionally, a large number of SSRs were obtained to improve our understanding of the phylogeny and genetics of Z. armatum. The organ-specific unigenes of the five major samples were screened and annotated via GO and KEGG enrichment analysis. A total of 53 and 34 unigenes that were exclusively upregulated in fruit samples were identified as candidate unigenes for terpenoid biosynthesis or fatty acid biosynthesis, elongation and degradation pathways, respectively. Moreover, 40 days after fertilization (Fr4 stage) could be an important period for the accumulation of terpenoid compounds during the fruit development and maturation of Z. armatum. The Fr4 stage could be a key point at which the first few steps of the fatty acid biosynthesis process are promoted, and the catalysis of subsequent reactions could be significantly induced at 62 days after fertilization (Fr6 stage).

Conclusions

The present study realized de novo transcriptome assembly for the five major organs of Z. armatum. To the best of our knowledge, this study provides the first comprehensive analysis revealing the genes underlying the special numbing taste of Z. armatum. The assembled transcriptome profiles expand the available genetic information on this species and will contribute to gene functional studies, which will aid in the engineering of high-quality cultivars of Z. armatum.