Identification, classification and differential expression of oleosin genes in tung tree (Vernicia fordii)

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.

Characterization of Oil Body and Starch Granule Dynamics in Developing Seeds of Brassica napus

Brassica napus is the most important oilseed crop in the world, and the lipid was stored in the oil body (OB) in the form of triacylglycerol. At present, most of studies on the relationship between oil body morphology and seed oil content in B. napus was focused on mature seeds. In the present study, the OBs in different developing seeds of B. napus with relatively high oil content (HOC) of about 50% and low oil content (LOC) of about 39% were analyzed. It was revealed that the size of OBs was first increased and then decreased in both materials. And in late seed developmental stages, the average OB size of rapeseed with HOC was higher than that of LOC, while it was reversed in the early seed developmental stages. No significant difference was observed on starch granule (SG) size in HOC and LOC rapeseed. Further results indicated that the expression of genes that involved in malonyl-CoA metabolism, fatty acid carbon chain extension, lipid metabolism, and starch synthesis in the rapeseed with HOC was higher than that of rapeseed with LOC. These results give some new insight for understanding the dynamics of OBs and SGs in embryos of B. napus.

Genomic insights into lineage-specific evolution of the oleosin family in Euphorbiaceae

BACKGROUND

Lipid droplets (LDs) present in land plants serve as an essential energy and carbon reserve for seed germination and seedling development. Oleosins, the most abundant structural proteins of LDs, comprise a small family involved in LD formation, stabilization and degradation. Despite their importance, our knowledge on oleosins is still poor in Euphorbiaceae, a large plant family that contains several important oil-bearing species.

RESULTS

To uncover lineage-specific evolution of oleosin genes in Euphorbiaceae, in this study, we performed a genome-wide identification and comprehensive comparison of the oleosin family in Euphorbiaceae species with available genome sequences, i.e. castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), Mercurialis annua, cassava (Manihot esculenta) and rubber tree (Hevea brasiliensis), and a number of five, five, five, five, eight and eight members were found, respectively. Synteny analysis revealed one-to-one collinear relationship of oleosin genes between the former four (i.e. castor bean, physic nut, tung tree and M. annua) as well as latter two species (i.e. cassava and rubber tree), whereas one-to-one and one-to-two collinear relationships were observed between physic nut and cassava, reflecting the occurrence of one recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber tree. The presence of five ortholog groups representing three previously defined clades (i.e. U, SL and SH) dates back at least to the Malpighiales ancestor, because they are also conserved in poplar (Populus trichocarpa), a tree having experienced one Salicaceae-specific recent WGD. As observed in poplar, WGD was shown to be the main driver for the family expansion in both cassava and rubber tree. Nevertheless, same retention patterns of WGD-derived duplicates observed in cassava and rubber tree are somewhat different from that of poplar, though certain homologous fragments are still present in rubber tree. Further transcriptional profiling revealed an apparent seed-predominant expression pattern of oleosin genes in physic nut, castor bean and rubber tree. Moreover, structure and expression divergence of paralogous pairs were also observed in both cassava and rubber tree.

CONCLUSION

Comparative genomics analysis of oleosin genes reported in this study improved our knowledge on lineage-specific family evolution in Euphorbiaceae, which also provides valuable information for further functional analysis and utilization of key members and their promoters.

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.

Tung Tree (Vernicia fordii) Genome Provides A Resource for Understanding Genome Evolution and Improved Oil Production

Tung tree (Vernicia fordii) is an economically important woody oil plant that produces tung oil rich in eleostearic acid. Here, we report a high-quality chromosome-scale genome sequence of tung tree. The genome sequence was assembled by combining Illumina short reads, Pacific Biosciences single-molecule real-time long reads, and Hi-C sequencing data. The size of tung tree genome is 1.12 Gb, with 28,422 predicted genes and over 73% repeat sequences. The V. fordii underwent an ancient genome triplication event shared by core eudicots but no further whole-genome duplication in the subsequent ca. 34.55 million years of evolutionary history of the tung tree lineage. Insertion time analysis revealed that repeat-driven genome expansion might have arisen as a result of long-standing long terminal repeat retrotransposon bursts and lack of efficient DNA deletion mechanisms. The genome harbors 88 resistance genes encoding nucleotide-binding sites; 17 of these genes may be involved in early-infection stage of Fusarium wilt resistance. Further, 651 oil-related genes were identified, 88 of which are predicted to be directly involved in tung oil biosynthesis. Relatively few phosphoenolpyruvate carboxykinase genes, and synergistic effects between transcription factors and oil biosynthesis-related genes might contribute to the high oil content of tung seed. The tung tree genome constitutes a valuable resource for understanding genome evolution, as well as for molecular breeding and genetic improvements for oil production.

Plant Lipid Droplets and Their Associated Proteins: Potential for Rapid Advances

Cytoplasmic lipid droplets (LDs) of neutral lipids (triacylglycerols [TAGs], sterylesters, etc.) are reserves of high-energy metabolites and other constituents for future needs. They are present in diverse cells of eukaryotes and prokaryotes. An LD has a core of neutral lipids enclosed with a monolayer of phospholipids and proteins, which play structural and/or metabolic roles. During the past 3 decades, studies of LDs in diverse organisms have blossomed after they were found to be involved in prevalent human diseases and industrial uses. LDs in plant seeds were studied before those in mammals and microbes, and the latter studies have since moved forward. Plant LDs carry a hallmark protein called oleosin, which has a long hydrophobic hairpin penetrating the TAG core and stabilizing the LD. The oleosin gene first appeared in green algae and has evolved in enhancing promoter strength, tandem repeats, and/or expression specificity, leading to the appearance of new LD organelles, such as tapetosomes in Brassicaceae. The synthesis of LDs occurs with TAG-synthesizing enzymes on the endoplasmic reticulum (ER), and nascent TAGs are sequestered in the acyl moiety region between the bilayers of phospholipids, which results in ER-LD swelling. Oleosin is synthesized on the cytosol side of the ER and extracts the LD from the ER-LD to cytosol. This extraction of LD to the cytosol is controlled solely by the innate properties of oleosin, and modified oleosin can redirect the LD to the ER lumen and then vacuoles. The breakdown of LDs requires lipase associating with core retromer and binding to peroxisomes, which then send the enzyme to LDs via tubular extensions. Two groups of LD-associated proteins, caleosin/dioxygenase/steroleosin and LD/oil body-associated proteins, participate in cellular stress defenses via enzymic activities and binding, respectively. The surface of LDs in all plant cells may be an inert refuge for these and other proteins, which exert functions on diverse cell components. Oleosin-LDs have been explored for commercial applications; successes in their uses will rely on overcoming conceptual and technical difficulties.

Relevance of Lipophilic Allergens in Food Allergy Diagnosis

PURPOSE OF REVIEW

The purpose of this review is to provide available data on a new class of allergens, the oleosins, and their diagnostic value. There is evidence that allergen extracts used for in vivo as well as in vitro diagnostic tests do not contain oleosins because these proteins are lipophilic and nearly insoluble in saline or aqueous solutions. So far, only oleosins of peanut, sesame and hazelnut have been registered as allergens. Reports on IgE-binding tests performed with oleosins of different species with sera from allergic patients show that IgE specific for oleosins are associated with severe allergic reactions which is why they should be part of the diagnostic tests in the future.

RECENT FINDINGS

Recent findings showed that oleosins purified from in shell-roasted peanuts revealed a higher IgE-binding capacity when compared to raw ones. Naturally purified as well as recombinantly produced peanut oleosins can be used in basophil activation test. The synopsis of all reports on different thermal processing of several oleosin sources and the respective data obtained with patients sera investigated via immunoblot and basophil activation test points to the recommendation that-if naturally purified oleosins are used, they should mostly be obtained from roasted food allergen sources. For immunoblot and basophil activation test, both, naturally purified oleosins as well as recombinant modified oleosins are valuable diagnostic tools.

Engineering the production of conjugated fatty acids in Arabidopsis thaliana leaves

The seeds of many nondomesticated plant species synthesize oils containing high amounts of a single unusual fatty acid, many of which have potential usage in industry. Despite the identification of enzymes for unusual oxidized fatty acid synthesis, the production of these fatty acids in engineered seeds remains low and is often hampered by their inefficient exclusion from phospholipids. Recent studies have established the feasibility of increasing triacylglycerol content in plant leaves, which provides a novel approach for increasing energy density of biomass crops. Here, we determined whether the fatty acid composition of leaf oil could be engineered to accumulate unusual fatty acids. Eleostearic acid (ESA) is a conjugated fatty acid produced in seeds of the tung tree (Vernicia fordii) and has both industrial and nutritional end-uses. Arabidopsis thaliana lines with elevated leaf oil were first generated by transforming wild-type, cgi-58 or pxa1 mutants (the latter two of which contain mutations disrupting fatty acid breakdown) with the diacylglycerol acyltransferases (DGAT1 or DGAT2) and/or oleosin genes from tung. High-leaf-oil plant lines were then transformed with tung FADX, which encodes the fatty acid desaturase/conjugase responsible for ESA synthesis. Analysis of lipids in leaves revealed that ESA was efficiently excluded from phospholipids, and co-expression of tung FADX and DGAT2 promoted a synergistic increase in leaf oil content and ESA accumulation. Taken together, these results provide a new approach for increasing leaf oil content that is coupled with accumulation of unusual fatty acids. Implications for production of biofuels, bioproducts, and plant-pest interactions are discussed.

Genome-Wide Analysis of Oleosin Gene Family in 22 Tree Species: An Accelerator for Metabolic Engineering of BioFuel Crops and Agrigenomics Industrial Applications?

Trees contribute to enormous plant oil reserves because many trees contain 50%-80% of oil (triacylglycerols, TAGs) in the fruits and kernels. TAGs accumulate in subcellular structures called oil bodies/droplets, in which TAGs are covered by low-molecular-mass hydrophobic proteins called oleosins (OLEs). The OLEs/TAGs ratio determines the size and shape of intracellular oil bodies. There is a lack of comprehensive sequence analysis and structural information of OLEs among diverse trees. The objectives of this study were to identify OLEs from 22 tree species (e.g., tung tree, tea-oil tree, castor bean), perform genome-wide analysis of OLEs, classify OLEs, identify conserved sequence motifs and amino acid residues, and predict secondary and three-dimensional structures in tree OLEs and OLE subfamilies. Data mining identified 65 OLEs with perfect conservation of the "proline knot" motif (PX5SPX3P) from 19 trees. These OLEs contained >40% hydrophobic amino acid residues. They displayed similar properties and amino acid composition. Genome-wide phylogenetic analysis and multiple sequence alignment demonstrated that these proteins could be classified into five OLE subfamilies. There were distinct patterns of sequence conservation among the OLE subfamilies and within individual tree species. Computational modeling indicated that OLEs were composed of at least three α-helixes connected with short coils without any β-strand and that they exhibited distinct 3D structures and ligand binding sites. These analyses provide fundamental information in the similarity and specificity of diverse OLE isoforms within the same subfamily and among the different species, which should facilitate studying the structure-function relationship and identify critical amino acid residues in OLEs for metabolic engineering of tree TAGs.

Characterization of the cork oak transcriptome dynamics during acorn development

BACKGROUND

Cork oak (Quercus suber L.) has a natural distribution across western Mediterranean regions and is a keystone forest tree species in these ecosystems. The fruiting phase is especially critical for its regeneration but the molecular mechanisms underlying the biochemical and physiological changes during cork oak acorn development are poorly understood. In this study, the transcriptome of the cork oak acorn, including the seed, was characterized in five stages of development, from early development to acorn maturation, to identify the dominant processes in each stage and reveal transcripts with important functions in gene expression regulation and response to water.

RESULTS

A total of 80,357 expressed sequence tags (ESTs) were de novo assembled from RNA-Seq libraries representative of the several acorn developmental stages. Approximately 7.6 % of the total number of transcripts present in Q. suber transcriptome was identified as acorn specific. The analysis of expression profiles during development returned 2,285 differentially expressed (DE) transcripts, which were clustered into six groups. The stage of development corresponding to the mature acorn exhibited an expression profile markedly different from other stages. Approximately 22 % of the DE transcripts putatively code for transcription factors (TF) or transcriptional regulators, and were found almost equally distributed among the several expression profile clusters, highlighting their major roles in controlling the whole developmental process. On the other hand, carbohydrate metabolism, the biological pathway most represented during acorn development, was especially prevalent in mid to late stages as evidenced by enrichment analysis. We further show that genes related to response to water, water deprivation and transport were mostly represented during the early (S2) and the last stage (S8) of acorn development, when tolerance to water desiccation is possibly critical for acorn viability.

CONCLUSIONS

To our knowledge this work represents the first report of acorn development transcriptomics in oaks. The obtained results provide novel insights into the developmental biology of cork oak acorns, highlighting transcripts putatively involved in the regulation of the gene expression program and in specific processes likely essential for adaptation. It is expected that this knowledge can be transferred to other oak species of great ecological value.

Identification and expression of fructose-1,6-bisphosphate aldolase genes and their relations to oil content in developing seeds of tea oil tree (Camellia oleifera)

Tea oil tree (Camellia oleifera, Co) provides a fine edible oil source in China. Tea oil from the seeds is very beneficial to human health. Fructose-1,6-bisphosphate aldolase (FBA) hydrolyzes fructose-1,6-bisphosphate into dihydroxyacetone phosphate and glyceraldehyde 3-phosphate, two critical metabolites for oil biosynthesis. The objectives of this study were to identify FBA genes and investigate the relationship between FBA gene expression and oil content in developing seeds of tea oil tree. In this paper, four developmentally up-regulated CoFBA genes were identified in Camellia oleifera seeds based on the transcriptome from two seed developmental stages corresponding to the initiation and peak stages of lipid biosynthesis. The expression of CoFBA genes, along with three key oil biosynthesis genes CoACP, CoFAD2 and CoSAD were analyzed in seeds from eight developmental stages by real-time quantitative PCR. The oil content and fatty acid composition were also analyzed. The results showed that CoFBA and CoSAD mRNA levels were well-correlated with oil content whereas CoFAD2 gene expression levels were correlated with fatty acid composition in Camellia seeds. We propose that CoFBA and CoSAD are two important factors for determining tea oil yield because CoFBA gene controls the flux of key intermediates for oil biosynthesis and CoSAD gene controls the synthesis of oleic acid, which accounts for 80% of fatty acids in tea oil. These findings suggest that tea oil yield could be improved by enhanced expression of CoFBA and CoSAD genes in transgenic plants.

Fatty acid profile and unigene-derived simple sequence repeat markers in tung tree (Vernicia fordii)

Tung tree (Vernicia fordii) provides the sole source of tung oil widely used in industry. Lack of fatty acid composition and molecular markers hinders biochemical, genetic and breeding research. The objectives of this study were to determine fatty acid profiles and develop unigene-derived simple sequence repeat (SSR) markers in tung tree. Fatty acid profiles of 41 accessions showed that the ratio of α-eleostearic acid was increasing continuously with a parallel trend to the amount of tung oil accumulation while the ratios of other fatty acids were decreasing in different stages of the seeds and that α-eleostearic acid (18∶3) consisted of 77% of the total fatty acids in tung oil. Transcriptome sequencing identified 81,805 unigenes from tung cDNA library constructed using seed mRNA and discovered 6,366 SSRs in 5,404 unigenes. The di- and tri-nucleotide microsatellites accounted for 92% of the SSRs with AG/CT and AAG/CTT being the most abundant SSR motifs. Fifteen polymorphic genic-SSR markers were developed from 98 unigene loci tested in 41 cultivated tung accessions by agarose gel and capillary electrophoresis. Genbank database search identified 10 of them putatively coding for functional proteins. Quantitative PCR demonstrated that all 15 polymorphic SSR-associated unigenes were expressed in tung seeds and some of them were highly correlated with oil composition in the seeds. Dendrogram revealed that most of the 41 accessions were clustered according to the geographic region. These new polymorphic genic-SSR markers will facilitate future studies on genetic diversity, molecular fingerprinting, comparative genomics and genetic mapping in tung tree. The lipid profiles in the seeds of 41 tung accessions will be valuable for biochemical and breeding studies.