Background: Vernonia galamensis native to Africa is an annual oleaginous plant of Asteraceae family. As a newly established industrial oil crop, this plant produces high level (> 70%) of vernolic acid (cis-12-epoxyoctadeca-cis-9-enoic acid), which is an unusual epoxy fatty acid (EFA) with multiple industrial applications. Here, transcriptome analysis and fatty acid profiling from developing V. galamensis seeds were integrated to uncover the critical metabolic pathways responsible for high EFA accumulation, aiming to identify the target genes that could be used in the biotechnological production of high-value oils. Results: Based on oil accumulation dynamics of V. galamensis seeds, we harvested seed samples from three stages (17, 38, and 45 days after pollination, DAP) representing the initial, fast and final EFA accumulation phases, and one mixed sample from different tissues for RNA-sequencing, with three biological replicates for each sample. Using Illumina platform, we have generated a total of 265 million raw cDNA reads. After filtering process, de novo assembly of clean reads yielded 67,114 unigenes with an N50 length of 1316 nt. Functional annotation resulted in the identification of almost all genes involved in diverse lipid-metabolic pathways, including the novel fatty acid desaturase/epoxygenase, diacylglycerol acyltransferases, and phospholipid:diacylglycerol acyltransferases. Expression profiling revealed that various genes associated with acyl editing, fatty acid β-oxidation, triacylglycerol assembly and oil-body formation had greater expression levels at middle developmental stage (38 DAP), which were consistent with the fast accumulation of EFA in V. galamensis developing seed, these genes were detected to play fundamental roles in EFA production. In addition, we isolated some transcription factors (such as WRI1, FUS3 and ABI4), which putatively regulated the production of V. galamensis seed oils. The transient expression of the selected genes resulted in a synergistic increase of EFA-enriched TAG accumulation in tobacco leaves. Transcriptome data were further confirmed by quantitative real-time PCR for twelve key genes in EFA biosynthesis. Finally, a comprehensive network for high EFA accumulation in V. galamensis seed was established. Conclusions: Our findings provide new insights into molecular mechanisms underlying the natural epoxy oil production in V. galamensis. A set of genes identified here could be used as the targets to develop other oilseeds highly accumulating valued epoxy oils for commercial production.

In spite of increasing use in the food industry, high relative levels of palmitic acid (C16:0) in cottonseed oil imposes harmful effects on human health when overconsumed in the diet. The limited understanding of the mechanism in controlling fatty acid composition has become a significant obstacle for breeding novel cotton varieties with high-quality oil. Fatty acyl–acyl carrier protein (ACP) thioesterase B (FatBs) are a group of enzymes which prefer to hydrolyze the thioester bond from saturated acyl-ACPs, thus playing key roles in controlling the accumulation of saturated fatty acids. However, FatB members and their roles in cotton are largely unknown. In this study, a genome-wide characterization of FatB members was performed in allotetraploid upland cotton, aiming to explore the GhFatBs responsible for high accumulations of C16:0 in cotton seeds. A total of 14 GhFatB genes with uneven distribution on chromosomes were identified from an upland cotton genome and grouped into seven subfamilies through phylogenetic analysis. The six key amino acid residues (Ala, Trys, Ile, Met, Arg and Try) responsible for substrate preference were identified in the N-terminal acyl binding pocket of GhFatBs. RNA-seq and qRT-PCR analysis revealed that the expression profiles of GhFatB genes varied in multiple cotton tissues, with eight GhFatBs (GhA/D-FatB3, GhA/D-FatB4, GhA/D-FatB5, and GhA/D-FatB7) having high expression levels in developing seeds. In particular, expression patterns of GhA-FatB3 and GhD-FatB4 were positively correlated with the dynamic accumulation of C16:0 during cotton seed development. Furthermore, heterologous overexpression assay of either GhA-FatB3 or GhD-FatB4 demonstrated that these two GhFatBs had a high substrate preference to 16:0-ACP, thus contributing greatly to the enrichment of palmitic acid in the tested tissues. Taken together, these findings increase our understanding on fatty acid accumulation and regulation mechanisms in plant seeds. GhFatBs, especially GhA-FatB3 and GhD-FatB4, could be molecular targets for genetic modification to reduce palmitic acid content or to optimize fatty acid profiles in cotton and other oil crops required for the sustainable production of healthy edible oil.