An Integrated Bioinformatics Analysis Reveals Divergent Evolutionary Pattern of Oil Biosynthesis in High- and Low-Oil Plants

Lipid Metabolism and Improvement in Oilseed Crops: Recent Advances in Multi-Omics Studies

Oilseed crops are rich in plant lipids that not only provide essential fatty acids for the human diet but also play important roles as major sources of biofuels and indispensable raw materials for the chemical industry. The regulation of lipid metabolism genes is a major factor affecting oil production. In this review, we systematically summarize the metabolic pathways related to lipid production and storage in plants and highlight key research advances in characterizing the genes and regulatory factors influencing lipid anabolic metabolism. In addition, we integrate the latest results from multi-omics studies on lipid metabolism to provide a reference to better understand the molecular mechanisms underlying oil anabolism in oilseed crops.

Acyl-CoA-dependent and acyl-CoA-independent avocado acyltransferases positively influence oleic acid content in nonseed triacylglycerols

In higher plants, acyl-CoA:diacylglycerol acyltransferase (DGAT) and phospholipid:diacylglycerol acyltransferase (PDAT) catalyze the terminal step of triacylglycerol (TAG) synthesis in acyl-CoA-dependent and -independent pathways, respectively. Avocado (Persea americana) mesocarp, a nonseed tissue, accumulates significant amounts of TAG (~70% by dry weight) that is rich in heart-healthy oleic acid (18:1). The oil accumulation stages of avocado mesocarp development coincide with high expression levels for type-1 DGAT (DGAT1) and PDAT1, although type-2 DGAT (DGAT2) expression remains low. The strong preference for oleic acid demonstrated by the avocado mesocarp TAG biosynthetic machinery represents lucrative biotechnological opportunities, yet functional characterization of these three acyltransferases has not been explored to date. We expressed avocado PaDGAT1, PaDGAT2, and PaPDAT1 in bakers' yeast and leaves of Nicotiana benthamiana. PaDGAT1 complemented the TAG biosynthesis deficiency in the quadruple mutant yeast strain H1246, and substantially elevated total cellular lipid content. In vitro enzyme assays showed that PaDGAT1 prefers oleic acid compared to palmitic acid (16:0). Both PaDGAT1 and PaPDAT1 increased the lipid content and elevated oleic acid levels when expressed independently or together, transiently in N. benthamiana leaves. These results indicate that PaDGAT1 and PaPDAT1 prefer oleate-containing substrates, and their coordinated expression likely contributes to sustained TAG synthesis that is enriched in oleic acid. This study establishes a knowledge base for future metabolic engineering studies focused on exploitation of the biochemical properties of PaDGAT1 and PaPDAT1.

Identification of QTN-by-environment interactions and their candidate genes for soybean seed oil-related traits using 3VmrMLM

INTRODUCTION

Although seed oil content and its fatty acid compositions in soybean were affected by environment, QTN-by-environment (QEIs) and gene-by-environment interactions (GEIs) were rarely reported in genome-wide association studies.

METHODS

The 3VmrMLM method was used to associate the trait phenotypes, measured in five to seven environments, of 286 soybean accessions with 106,013 SNPs for detecting QTNs and QEIs.

RESULTS

Seven oil metabolism genes (GmSACPD-A, GmSACPD-B, GmbZIP123, GmSWEET39, GmFATB1A, GmDGAT2D, and GmDGAT1B) around 598 QTNs and one oil metabolism gene GmFATB2B around 54 QEIs were verified in previous studies; 76 candidate genes and 66 candidate GEIs were predicted to be associated with these traits, in which 5 genes around QEIs were verified in other species to participate in oil metabolism, and had differential expression across environments. These genes were found to be related to soybean seed oil content in haplotype analysis. In addition, most candidate GEIs were co-expressed with drought response genes in co-expression network, and three KEGG pathways which respond to drought were enriched under drought stress rather than control condition; six candidate genes were hub genes in the co-expression networks under drought stress.

DISCUSSION

The above results indicated that GEIs, together with drought response genes in co-expression network, may respond to drought, and play important roles in regulating seed oil-related traits together with oil metabolism genes. These results provide important information for genetic basis, molecular mechanisms, and soybean breeding for seed oil-related traits.

An expanded role for the transcription factor WRINKLED1 in the biosynthesis of triacylglycerols during seed development

The transcription factor WRINKLED1 (WRI1) is known as a master regulator of fatty acid synthesis in developing oilseeds of Arabidopsis thaliana and other species. WRI1 is known to directly stimulate the expression of many fatty acid biosynthetic enzymes and a few targets in the lower part of the glycolytic pathway. However, it remains unclear to what extent and how the conversion of sugars into fatty acid biosynthetic precursors is controlled by WRI1. To shortlist possible gene targets for future in-planta experimental validation, here we present a strategy that combines phylogenetic foot printing of cis-regulatory elements with additional layers of evidence. Upstream regions of protein-encoding genes in A. thaliana were searched for the previously described DNA-binding consensus for WRI1, the ASML1/WRI1 (AW)-box. For about 900 genes, AW-box sites were found to be conserved across orthologous upstream regions in 11 related species of the crucifer family. For 145 select potential target genes identified this way, affinity of upstream AW-box sequences to WRI1 was assayed by Microscale Thermophoresis. This allowed definition of a refined WRI1 DNA-binding consensus. We find that known WRI1 gene targets are predictable with good confidence when upstream AW-sites are phylogenetically conserved, specifically binding WRI1 in the in vitro assay, positioned in proximity to the transcriptional start site, and if the gene is co-expressed with WRI1 during seed development. When targets predicted in this way are mapped to central metabolism, a conserved regulatory blueprint emerges that infers concerted control of contiguous pathway sections in glycolysis and fatty acid biosynthesis by WRI1. Several of the newly predicted targets are in the upper glycolysis pathway and the pentose phosphate pathway. Of these, plastidic isoforms of fructokinase (FRK3) and of phosphoglucose isomerase (PGI1) are particularly corroborated by previously reported seed phenotypes of respective null mutations.

An integrated omics analysis reveals the gene expression profiles of maize, castor bean, and rapeseed for seed oil biosynthesis

BACKGROUND

Seed storage lipids are valuable for human diet and for the sustainable development of mankind. In recent decades, many lipid metabolism genes and pathways have been identified, but the molecular mechanisms that underlie differences in seed oil biosynthesis in species with developed embryo and endosperm are not fully understood.

RESULTS

We performed comparative genome and transcriptome analyses of castor bean and rapeseed, which have high seed oil contents, and maize, which has a low seed oil content. These results revealed the molecular underpinnings of the low seed oil content in maize. First of all, transcriptome analyses showed that more than 61% of the lipid- and carbohydrate-related genes were regulated in castor bean and rapeseed, but only 20.1% of the lipid-related genes and 22.5% of the carbohydrate-related genes were regulated in maize. Then, compared to castor bean and rapeseed, fewer lipid biosynthesis genes but more lipid metabolism genes were regulated in the maize embryo. More importantly, most maize genes encoding lipid-related transcription factors, triacylglycerol (TAG) biosynthetic enzymes, pentose phosphate pathway (PPP) and Calvin Cycle proteins were not regulated during seed oil synthesis, despite the presence of many homologs in the maize genome. Additionally, we observed differential regulation of vital oil biosynthetic enzymes and extremely high expression levels of oil biosynthetic genes in castor bean, which were consistent with the rapid accumulation of oil in castor bean developing seeds.

CONCLUSIONS

Compared to high-oil seeds (castor bean and rapeseed), less oil biosynthetic genes were regulated during the seed development in low-oil seed (maize). These results shed light on molecular mechanisms of lipid biosynthesis in maize, castor bean, and rapeseed. They can provide information on key target genes that may be useful for future experimental manipulation of oil production in oil plants.

Identification and analysis of oil candidate genes reveals the molecular basis of cottonseed oil accumulation in Gossypium hirsutum L

Based on the integration of QTL-mapping and regulatory network analyses, five high-confidence stable QTL regions, six candidate genes and two microRNAs that potentially affect the cottonseed oil content were discovered. Cottonseed oil is increasingly becoming a promising target for edible oil with its high content of unsaturated fatty acids. In this study, a recombinant inbred line (RIL) cotton population was constructed to detect quantitative trait loci (QTLs) for the cottonseed oil content. A total of 39 QTLs were detected across eight different environments, of which five QTLs were stable. Forty-three candidate genes potentially involved in carbon metabolism, fatty acid synthesis and triacylglycerol biosynthesis processes were further obtained in the stable QTL regions. Transcriptome analysis showed that nineteen of these candidate genes expressed during  the developing cottonseed ovules and may affect the cottonseed oil content. Besides, transcription factor (TF) and microRNA (miRNA) co-regulatory network analyses based on the nineteen candidate genes suggested that six genes, two core miRNAs (ghr-miR2949b and ghr-miR2949c), and one TF GhHSL1 were considered to be closely associated with the cottonseed oil content. Moreover, four vital genes were validated by quantitative real-time PCR (qRT-PCR). These results provide insights into the oil accumulation mechanism in developing cottonseed ovules through the construction of a detailed oil accumulation model.

Integrating omics approaches to discover and prioritize candidate genes involved in oil biosynthesis in soybean

Soybean is a major source of edible protein and oil. Oil content is a quantitative trait that is significantly determined by genetic and environmental factors. Over the past 30 years, a large volume of soybean genetic, genomic, and transcriptomic data have been accumulated. Nevertheless, integrative analyses of such data remain scarce, in spite of their importance for crop improvement. We hypothesized that the co-occurrence of genomic regions for oil-related traits in different studies may reveal more stable regions encompassing important genetic determinants of oil content and quality in soybean. We integrated publicly available data, obtained with distinct techniques, to discover and prioritize candidate genes involved in oil biosynthesis and regulation in soybean. We detected key fatty acid biosynthesis genes (e.g., BCCP2 and ACCase, FADs, KAS family proteins) and several transcription factors, which are likely regulators of oil biosynthesis. In addition, we identified new candidates for seed oil accumulation and quality, such as Glyma.03G213300 and Glyma.19G160700, which encode a translocator protein homolog and a histone acetyltransferase, respectively. Further, oil and protein genomic hotspots are strongly associated with breeding and not with domestication, suggesting that soybean domestication prioritized other traits. The genes identified here are promising targets for breeding programs and for the development of soybean lines with increased oil content and quality.

Domestication and improvement genes reveal the differences of seed size- and oil-related traits in soybean domestication and improvement

•

Due to reduced diversity, it is essential to map domesticated and improved genes.

•

13 known and 442 candidate genes were mined for seed size- and oil-related traits.

•

All the genes were used to explain trait changes in domestication and improvement.

•

56 domesticated and 15 improved genes may be valuable for future soybean breeding.

•

This study provides useful gene resources for future breeding and biology research.

To address domestication and improvement studies of soybean seed size- and oil-related traits, a series of domesticated and improved regions, loci, and candidate genes were identified in 286 soybean accessions using domestication and improvement analyses, genome-wide association studies, quantitative trait locus (QTL) mapping and bulked segregant analyses in this study. As a result, 534 candidate domestication regions (CDRs) and 458 candidate improvement regions (CIRs) were identified in this study and integrated with those in five and three previous studies, respectively, to obtain 952 CDRs and 538 CIRs; 1469 loci for soybean seed size- and oil-related traits were identified in this study and integrated with those in Soybase to obtain 433 QTL clusters. The two results were intersected to obtain 245 domestication and 221 improvement loci for the above traits. Around these trait-related domestication and improvement loci, 7 domestication and 7 improvement genes were found to be truly associated with these traits, and 372 candidate domestication and 87 candidate improvement genes were identified using gene expression, SNP variants in genome, miRNA binding, KEGG pathway, DNA methylation, and haplotype analysis. These genes were used to explain the trait changes in domestication and improvement. As a result, the trait changes can be explained by their frequencies of elite haplotypes, base mutations in coding region, and three factors affecting their expression levels. In addition, 56 domestication and 15 improvement genes may be valuable for future soybean breeding. This study can provide useful gene resources for future soybean breeding and molecular biology research.

High oil accumulation in tuber of yellow nutsedge compared to purple nutsedge is associated with more abundant expression of genes involved in fatty acid synthesis and triacylglycerol storage

BACKGROUND

Yellow nutsedge is a unique plant species that can accumulate up to 35% oil of tuber dry weight, perhaps the highest level observed in the tuber tissues of plant kingdom. To gain insight into the molecular mechanism that leads to high oil accumulation in yellow nutsedge, gene expression profiles of oil production pathways involved carbon metabolism, fatty acid synthesis, triacylglycerol synthesis, and triacylglycerol storage during tuber development were compared with purple nutsedge, the closest relative of yellow nutsedge that is poor in oil accumulation.

RESULTS

Compared with purple nutsedge, high oil accumulation in yellow nutsedge was associated with significant up-regulation of specific key enzymes of plastidial RubisCO bypass as well as malate and pyruvate metabolism, almost all fatty acid synthesis enzymes, and seed-like oil-body proteins. However, overall transcripts for carbon metabolism toward carbon precursor for fatty acid synthesis were comparable and for triacylglycerol synthesis were similar in both species. Two seed-like master transcription factors ABI3 and WRI1 were found to display similar transcript patterns but were expressed at 6.5- and 14.3-fold higher levels in yellow nutsedge than in purple nutsedge, respectively. A weighted gene co-expression network analysis revealed that ABI3 was in strong transcriptional coordination with WRI1 and other key oil-related genes.

CONCLUSIONS

These results implied that pyruvate availability and fatty acid synthesis in plastid, along with triacylglycerol storage in oil bodies, rather than triacylglycerol synthesis in endoplasmic reticulum, are the major factors responsible for high oil production in tuber of yellow nutsedge, and ABI3 most likely plays a critical role in regulating oil accumulation. This study is of significance with regard to understanding the molecular mechanism controlling carbon partitioning toward oil production in oil-rich tuber and provides a valuable reference for enhancing oil accumulation in non-seed tissues of crops through genetic breeding or metabolic engineering.

An evolutionary population structure model reveals pleiotropic effects of GmPDAT for traits related to seed size and oil content in soybean

Seed oil traits in soybean that are of benefit to human nutrition and health have been selected for during crop domestication. However, these domesticated traits have significant differences across various evolutionary types. In this study, we found that the integration of evolutionary population structure (evolutionary types) with genome-wide association studies increased the power of gene detection, and it identified one locus for traits related to seed size and oil content on chromosome 13. This domestication locus, together with another one in a 200-kb region, was confirmed by the GEMMA and EMMAX software. The candidate gene, GmPDAT, had higher expressional levels in high-oil and large-seed accessions than in low-oil and small-seed accessions. Overexpression lines had increased seed size and oil content, whereas RNAi lines had decreased seed size and oil content. The molecular mechanism of GmPDAT was deduced based on results from linkage analysis for triacylglycerols and on histocytological comparisons of transgenic soybean seeds. Our results illustrate a new approach for identifying domestication genes with pleiotropic effects.

Three-dimensional genetic networks among seed oil-related traits, metabolites and genes reveal the genetic foundations of oil synthesis in soybean

Although the biochemical and genetic basis of lipid metabolism is clear in Arabidopsis, there is limited information concerning the relevant genes in Glycine max (soybean). To address this issue, we constructed three-dimensional genetic networks using six seed oil-related traits, 52 lipid metabolism-related metabolites and 54 294 SNPs in 286 soybean accessions in total. As a result, 284 and 279 candidate genes were found to be significantly associated with seed oil-related traits and metabolites by phenotypic and metabolic genome-wide association studies and multi-omics analyses, respectively. Using minimax concave penalty (MCP) and smoothly clipped absolute deviation (SCAD) analyses, six seed oil-related traits were found to be significantly related to 31 metabolites. Among the above candidate genes, 36 genes were found to be associated with oil synthesis (27 genes), amino acid synthesis (four genes) and the tricarboxylic acid (TCA) cycle (five genes), and four genes (GmFATB1a, GmPDAT, GmPLDα1 and GmDAGAT1) are already known to be related to oil synthesis. Using this information, 133 three-dimensional genetic networks were constructed, 24 of which are known, e.g. pyruvate-GmPDAT-GmFATA2-oil content. Using these networks, GmPDAT, GmAGT and GmACP4 reveal the genetic relationships between pyruvate and the three major nutrients, and GmPDAT, GmZF351 and GmPgs1 reveal the genetic relationships between amino acids and seed oil content. In addition, GmCds1, along with average temperature in July and the rainfall from June to September, influence seed oil content across years. This study provides a new approach for the construction of three-dimensional genetic networks and reveals new information for soybean seed oil improvement and the identification of gene function.

Effect of marker segregation distortion on high density linkage map construction and QTL mapping in Soybean (Glycine max L.)

Marker segregation distortion is a natural phenomenon. Severely distorted markers are usually excluded in the construction of linkage maps. We investigated the effect of marker segregation distortion on linkage map construction and quantitative trait locus (QTL) mapping. A total of 519 recombinant inbred lines of soybean from orthogonal and reciprocal crosses between LSZZH and NN493-1 were genotyped by specific length amplified fragment markers and seed linoleic acid content was measured in three environments. As a result, twenty linkage groups were constructed with 11,846 markers, including 1513 (12.77%) significantly distorted markers, on 20 chromosomes, and the map length was 2475.86 cM with an average marker-interval of 0.21 cM. The inclusion of distorted markers in the analysis was shown to not only improve the grouping of the markers from the same chromosomes, and the consistency of linkage maps with genome, but also increase genome coverage by markers. Combining genotypic data from both orthogonal and reciprocal crosses decreased the proportion of distorted markers and then improved the quality of linkage maps. Validation of the linkage maps was confirmed by the high collinearity between positions of markers in the soybean reference genome and in linkage maps and by the high consistency of 24 QTL regions in this study compared with the previously reported QTLs and lipid metabolism related genes. Additionally, linkage maps that include distorted markers could add more information to the outputs from QTL mapping. These results provide important information for linkage mapping, gene cloning and marker-assisted selection in soybean.

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

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

Genetic control of fatty acid composition in coconut (Cocos nucifera), African oil palm (Elaeis guineensis), and date palm (Phoenix dactylifera)

Predominant gene isoforms and expression bias in lipid metabolism pathways are highly conserved between oil-producing Arecaceae crop species coconut and oil palm, but diverge in non-oil-producing species date palm. Coconut (Cocos nucifera), African oil palm (Elaeis guineensis) and date palm (Phoenix dactylifera) are three major crop species in the Arecaceae family for which genome sequences have recently become available. Coconut and African oil palm both store oil in their endosperms, while date palm fruits contain very little oil. We analyzed fatty acid composition in three coconut tissues (leaf, endosperm and embryo) and in two African oil palm tissues (leaf and mesocarp), and identified 806, 840 and 848 lipid-related genes in 22 lipid metabolism pathways from the coconut, African oil palm and date palm genomes, respectively. The majority of lipid-related genes were highly homologous and retained in homologous segments between the three species. Genes involved in the conversion of pyruvate to fatty acid had a five-to-sixfold higher expression in the coconut endosperm and oil palm mesocarp than in the leaf or embryo tissues based on Fragments Per Kilobase of transcript per Million mapped reads values. A close evolutionary relationship between predominant gene isoforms and high conservation of gene expression bias in the lipid and carbohydrate gene metabolism pathways was observed for the two oil-producing species coconut and oil palm, differing from that of date palm, a non-oil-producing species. Our results elucidate the similarities and differences in lipid metabolism between the three major Arecaceae crop species, providing important information for physiology studies as well as breeding for fatty acid composition and oil content in these crops.

An integrated omics analysis reveals molecular mechanisms that are associated with differences in seed oil content between Glycine max and Brassica napus

BACKGROUND

Rapeseed (Brassica napus L.) and soybean (Glycine max L.) seeds are rich in both protein and oil, which are major sources of biofuels and nutrition. Although the difference in seed oil content between soybean (~ 20%) and rapeseed (~ 40%) exists, little is known about its underlying molecular mechanism.

RESULTS

An integrated omics analysis was performed in soybean, rapeseed, Arabidopsis (Arabidopsis thaliana L. Heynh), and sesame (Sesamum indicum L.), based on Arabidopsis acyl-lipid metabolism- and carbon metabolism-related genes. As a result, candidate genes and their transcription factors and microRNAs, along with phylogenetic analysis and co-expression network analysis of the PEPC gene family, were found to be largely associated with the difference between the two species. First, three soybean genes (Glyma.13G148600, Glyma.13G207900 and Glyma.12G122900) co-expressed with GmPEPC1 are specifically enriched during seed storage protein accumulation stages, while the expression of BnPEPC1 is putatively inhibited by bna-miR169, and two genes BnSTKA and BnCKII are co-expressed with BnPEPC1 and are specifically associated with plant circadian rhythm, which are related to seed oil biosynthesis. Then, in de novo fatty acid synthesis there are rapeseed-specific genes encoding subunits β-CT (BnaC05g37990D) and BCCP1 (BnaA03g06000D) of heterogeneous ACCase, which could interfere with synthesis rate, and β-CT is positively regulated by four transcription factors (BnaA01g37250D, BnaA02g26190D, BnaC01g01040D and BnaC07g21470D). In triglyceride synthesis, GmLPAAT2 is putatively inhibited by three miRNAs (gma-miR171, gma-miR1516 and gma-miR5775). Finally, in rapeseed there was evidence for the expansion of gene families, CALO, OBO and STERO, related to lipid storage, and the contraction of gene families, LOX, LAH and HSI2, related to oil degradation.

CONCLUSIONS

The molecular mechanisms associated with differences in seed oil content provide the basis for future breeding efforts to improve seed oil content.

Selection for a Zinc-Finger Protein Contributes to Seed Oil Increase during Soybean Domestication

Seed oil is a momentous agronomical trait of soybean (Glycine max) targeted by domestication in breeding. Although multiple oil-related genes have been uncovered, knowledge of the regulatory mechanism of seed oil biosynthesis is currently limited. We demonstrate that the seed-preferred gene GmZF351, encoding a tandem CCCH zinc finger protein, is selected during domestication. Further analysis shows that GmZF351 facilitates oil accumulation by directly activating WRINKLED1, BIOTIN CARBOXYL CARRIER PROTEIN2, 3-KETOACYL-ACYL CARRIER PROTEIN SYNTHASE III, DIACYLGLYCEROL O-ACYLTRANSFERASE1, and OLEOSIN2 in transgenic Arabidopsis (Arabidopsis thaliana) seeds. Overexpression of GmZF351 in transgenic soybean also activates lipid biosynthesis genes, thereby accelerating seed oil accumulation. The ZF351 haplotype from the cultivated soybean group and the wild soybean (Glycine soja) subgroup III correlates well with high gene expression level, seed oil contents and promoter activity, suggesting that selection of GmZF351 expression leads to increased seed oil content in cultivated soybean. Our study provides novel insights into the regulatory mechanism for seed oil accumulation, and the manipulation of GmZF351 may have great potential in the improvement of oil production in soybean and other related crops.

Adaptive evolution of seed oil content in angiosperms: accounting for the global patterns of seed oils

Background

Studies of the biogeographic distribution of seed oil content in plants are fundamental to understanding the mechanisms of adaptive evolution in plants as seed oil is the primary energy source needed for germination and establishment of plants. However, seed oil content as an adaptive trait in plants is poorly understood. Here, we examine the adaptive nature of seed oil content in 168 angiosperm families occurring in different biomes across the world. We also explore the role of multiple seed traits like seed oil content and composition in plant adaptation in a phylogenetic and nonphylogenetic context.

Result

It was observed that the seed oil content in tropical plants (28.4 %) was significantly higher than the temperate plants (24.6 %). A significant relationship between oil content and latitude was observed in three families Papaveraceae, Sapindaceae and Sapotaceae indicating that selective forces correlated with latitude influence seed oil content. Evaluation of the response of seed oil content and composition to latitude and the correlation between seed oil content and composition showed that multiple seed traits, seed oil content and composition contribute towards plant adaptation. Investigation of the presence or absence of phylogenetic signals across 168 angiosperm families in 62 clades revealed that members of seven clades evolved to have high or low seed oil content independently as they did not share a common evolutionary path.

Conclusion

The study provides us an insight into the biogeographical distribution and the adaptive role of seed oil content in plants. The study indicates that multiple seed traits like seed oil content and the fatty acid composition of the seed oils determine the fitness of the plants and validate the adaptive hypothesis that seed oil quantity and quality are crucial to plant adaptation.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0752-7) contains supplementary material, which is available to authorized users.