Time-Series Analyses of Transcriptomes and Proteomes Reveal Molecular Networks Underlying Oil Accumulation in Canola

Analysis of WAK Genes in Nine Cruciferous Species with a Focus on Brassica napus L

The wall-associated kinase family contributes to plant cell elongation and pathogen recognition. Nine Cruciferous species were studied for identification and molecular evolution of the WAK gene family. Firstly, 178 WAK genes were identified. A phylogenetic tree was constructed of the Cruciferous WAK proteins into four categories, of which the Brassica rapa, Brassica oleracea and Brassica napus genes in the U's triangle were more closely related. The WAK gene family was unevenly distributed in B. napus chromosomal imaging, with the largest number of BnWAK genes located on chromosome C08. In the expression analysis, the expression patterns of the WAK gene family varied under different stress treatments, and some members of BnWAKs were significantly different under stress treatments. This study lays a foundation for further revealing the functional mechanisms of the WAK gene family in Brassica napus.

Transcriptome profiling, physiological, and biochemical analyses provide new insights towards drought stress response in sugar maple (Acer saccharum Marshall) saplings

Sugar maple (Acer saccharum Marshall) is a temperate tree species in the northeastern parts of the United States and is economically important for its hardwood and syrup production. Sugar maple trees are highly vulnerable to changing climatic conditions, especially drought, so understanding the physiological, biochemical, and molecular responses is critical. The sugar maple saplings were subjected to drought stress for 7, 14, and 21 days and physiological data collected at 7, 14, and 21 days after stress (DAS) showed significantly reduced chlorophyll and Normalized Difference Vegetation Index with increasing drought stress time. The drought stress-induced biochemical changes revealed a higher accumulation of malondialdehyde, proline, and peroxidase activity in response to drought stress. Transcriptome analysis identified a total of 14,099 differentially expressed genes (DEGs); 328 were common among all stress periods. Among the DEGs, transcription factors (including NAC, HSF, ZFPs, GRFs, and ERF), chloroplast-related and stress-responsive genes such as peroxidases, membrane transporters, kinases, and protein detoxifiers were predominant. GO enrichment and KEGG pathway analysis revealed significantly enriched processes related to protein phosphorylation, transmembrane transport, nucleic acids, and metabolic, secondary metabolite biosynthesis pathways, circadian rhythm-plant, and carotenoid biosynthesis in response to drought stress. Time-series transcriptomic analysis revealed changes in gene regulation patterns in eight different clusters, and pathway analysis by individual clusters revealed a hub of stress-responsive pathways. In addition, qRT-PCR validation of selected DEGs revealed that the expression patterns were consistent with transcriptome analysis. The results from this study provide insights into the dynamics of physiological, biochemical, and gene responses to progressive drought stress and reveal the important stress-adaptive mechanisms of sugar maple saplings in response to drought stress.

Genome-Wide Identification and Characterization of Oil-Body-Membrane Proteins in Polyploid Crop Brassica napus

Oil-body-membrane proteins (OBMPs) are essential structural molecules of oil bodies and also versatile metabolic enzymes involved in multiple cellular processes such as lipid metabolism, hormone signaling and stress responses. However, the global landscape for OBMP genes in oil crops is still lacking. Here, we performed genome-wide identification and characterization of OBMP genes in polyploid crop Brassica napus. B. napus contains up to 88 BnaOBMP genes including 53 oleosins, 20 caleosins and 15 steroleosins. Both whole-genome and tandem duplications have contributed to the expansion of the BnaOBMP gene family. These BnaOBMP genes have extensive sequence polymorphisms, and some harbor strong selection signatures. Various cis-acting regulatory elements involved in plant growth, phytohormones and abiotic and biotic stress responses are detected in their promoters. BnaOBMPs exhibit differential expression at various developmental stages from diverse tissues. Importantly, some BnaOBMP genes display spatiotemporal patterns of seed-specific expression, which could be orchestrated by transcriptional factors such as EEL, GATA3, HAT2, SMZ, DOF5.6 and APL. Altogether, our data lay the foundations for studying the regulatory mechanism of the seed oil storage process and provide candidate genes and alleles for the genetic improvement and breeding of rapeseed with high seed oil content.

Hairy root transformation system as a tool for CRISPR/Cas9-directed genome editing in oilseed rape (Brassica napus)

Our study examined the mutation efficiency of the CRISPR/Cas9 method for tryptophan aminotransferase BnaTAA1 genes involved in the auxin biosynthesis pathway. We made nine CRISPR/Cas9 constructs with various promoters driving the expression of a Cas9 from Staphylococcus aureus (SaCas9) or a plant-codon-optimized Streptococcus pyogenes Cas9 (pcoCas9). We developed a fast and efficient system for evaluating the variety and frequency of mutations caused by each construct using Brassica napus hairy roots. We showed that pcoCas9 is more efficient in mutating the targeted loci than SaCas9 and the presence of the NLS signal enhanced the chance of mutagenesis by 25%. The mutations were studied further in regenerated lines, and we determined the BnaTAA1 gene expression and heritability of the gene modifications in transgenic plants. Hairy root transformation combined with CRISPR/Cas9-mediated gene editing represents a fast and straightforward system for studying target gene function in the important oilseed crop B. napus.

Three strategies of transgenic manipulation for crop improvement

Heterologous expression of exogenous genes, overexpression of endogenous genes, and suppressed expression of undesirable genes are the three strategies of transgenic manipulation for crop improvement. Up to 2020, most (227) of the singular transgenic events (265) of crops approved for commercial release worldwide have been developed by the first strategy. Thirty-eight of them have been transformed by synthetic sequences transcribing antisense or double-stranded RNAs and three by mutated copies for suppressed expression of undesirable genes (the third strategy). By the first and the third strategies, hundreds of transgenic events and thousands of varieties with significant improvement of resistance to herbicides and pesticides, as well as nutritional quality, have been developed and approved for commercial release. Their application has significantly decreased the use of synthetic pesticides and the cost of crop production and increased the yield of crops and the benefits to farmers. However, almost all the events overexpressing endogenous genes remain at the testing stage, except one for fertility restoration and another for pyramiding herbicide tolerance. The novel functions conferred by the heterologously expressing exogenous genes under the control of constitutive promoters are usually absent in the recipient crops themselves or perform in different pathways. However, the endogenous proteins encoded by the overexpressing endogenous genes are regulated in complex networks with functionally redundant and replaceable pathways and are difficult to confer the desirable phenotypes significantly. It is concluded that heterologous expression of exogenous genes and suppressed expression by RNA interference and clustered regularly interspaced short palindromic repeats-cas (CRISPR/Cas) of undesirable genes are superior to the overexpression of endogenous genes for transgenic improvement of crops.

Genomic selection and genetic architecture of agronomic traits during modern rapeseed breeding

Rapeseed (Brassica napus L.) is an important oil-producing crop for the world. Its adaptation, yield and quality have been considerably improved in recent decades, but the genomic basis underlying successful breeding selection remains unclear. Hence, we conducted a comprehensive genomic assessment of rapeseed in the breeding process based on the whole-genome resequencing of 418 diverse rapeseed accessions. We unraveled the genomic basis for the selection of adaptation and agronomic traits. Genome-wide association studies identified 628 associated loci-related causative candidate genes for 56 agronomically important traits, including plant architecture and yield traits. Furthermore, we uncovered nonsynonymous mutations in plausible candidate genes for agronomic traits with significant differences in allele frequency distributions across the improvement process, including the ribosome recycling factor (BnRRF) gene for seed weight. This study provides insights into the genomic basis for improving rapeseed varieties and a valuable genomic resource for genome-assisted rapeseed breeding.

Myxospermy Evolution in Brassicaceae: A Highly Complex and Diverse Trait with Arabidopsis as an Uncommon Model

The ability to extrude mucilage upon seed imbibition (myxospermy) occurs in several Angiosperm taxonomic groups, but its ancestral nature or evolutionary convergence origin remains misunderstood. We investigated seed mucilage evolution in the Brassicaceae family with comparison to the knowledge accumulated in Arabidopsis thaliana. The myxospermy occurrence was evaluated in 27 Brassicaceae species. Phenotyping included mucilage secretory cell morphology and topochemistry to highlight subtle myxospermy traits. In parallel, computational biology was driven on the one hundred genes constituting the so-called A. thaliana mucilage secretory cell toolbox to confront their sequence conservation to the observed phenotypes. Mucilage secretory cells show high morphology diversity; the three studied Arabidopsis species had a specific extrusion modality compared to the other studied Brassicaceae species. Orthologous genes from the A. thaliana mucilage secretory cell toolbox were mostly found in all studied species without correlation with the occurrence of myxospermy or even more sub-cellular traits. Seed mucilage may be an ancestral feature of the Brassicaceae family. It consists of highly diverse subtle traits, probably underlined by several genes not yet characterized in A. thaliana or by species-specific genes. Therefore, A. thaliana is probably not a sufficient reference for future myxospermy evo-devo studies.

Transcriptomic Analysis Reveals Key Genes Involved in Oil and Linoleic Acid Biosynthesis during Artemisia sphaerocephala Seed Development

Artemisia sphaerocephala seeds are rich in polysaccharides and linoleic acid (C18:2), which have been widely used as traditional medicine and to improve food quality. The accumulation patterns and molecular regulatory mechanisms of polysaccharides during A. sphaerocephala seed development have been studied. However, the related research on seed oil and C18:2 remain unclear. For this study, A. sphaerocephala seeds at seven different development stages at 10, 20, 30, 40, 50, 60, and 70 days after flowering (designated as S1~S7), respectively, were employed as experimental samples, the accumulation patterns of oil and fatty acids (FA) and the underlying molecular regulatory mechanisms were analyzed. The results revealed that oil content increased from 10.1% to 20.0% in the early stages of seed development (S1~S2), and up to 32.0% in mature seeds, of which C18:2 accounted for 80.6% of the total FA. FA and triacylglycerol biosynthesis-related genes jointly involved in the rapid accumulation of oil in S1~S2. Weighted gene co-expression network analysis showed that transcription factors FUS3 and bHLH played a critical role in the seed oil biosynthesis. The perfect harmonization of the high expression of FAD2 with the extremely low expression of FAD3 regulated the accumulation of C18:2. This study uncovered the gene involved in oil biosynthesis and molecular regulatory mechanisms of high C18:2 accumulation in A. sphaerocephala seeds; thus, advancing research into unsaturated fatty acid metabolism in plants while generating valuable genetic resources for optimal C18:2 breeding.

Disparate Dynamics of Gene Body and cis-Regulatory Element Evolution Illustrated for the Senescence-Associated Cysteine Protease Gene SAG12 of Plants

Gene regulation networks precisely orchestrate the expression of genes that are closely associated with defined physiological and developmental processes such as leaf senescence in plants. The Arabidopsis thaliana senescence-associated gene 12 (AtSAG12) encodes a cysteine protease that is (i) involved in the degradation of chloroplast proteins and (ii) almost exclusively expressed during senescence. Transcription factors, such as WRKY53 and WRKY45, bind to W-boxes in the promoter region of AtSAG12 and play key roles in its activation. Other transcription factors, such as bZIPs, might have accessory functions in their gene regulation, as several A-boxes have been identified and appear to be highly overrepresented in the promoter region compared to the whole genome distribution but are not localized within the regulatory regions driving senescence-associated expression. To address whether these two regulatory elements exhibiting these different properties are conserved in other closely related species, we constructed phylogenetic trees of the coding sequences of orthologs of AtSAG12 and screened their respective 2000 bp promoter regions for the presence of conserved cis-regulatory elements, such as bZIP and WRKY binding sites. Interestingly, the functional relevant upstream located W-boxes were absent in plant species as closely related as Arabidopsis lyrata, whereas an A-box cluster appeared to be conserved in the Arabidopsis species but disappeared in Brassica napus. Several orthologs were present in other species, possibly because of local or whole genome duplication events, but with distinct cis-regulatory sites in different locations. However, at least one gene copy in each family analyzed carried one W-box and one A-box in its promoter. These gene differences in SAG12 orthologs are discussed in the framework of cis- and trans-regulatory factors, of promoter and gene evolution, of genetic variation, and of the enhancement of the adaptability of plants to changing environmental conditions.

A global survey of the gene network and key genes for oil accumulation in cultivated tetraploid cottons

To enrich our knowledge about gene network of fatty acid biosynthesis in cottonseed, we conducted comparative transcriptome to reveal the differences in gene expression between Gossypium hirsutum and Gossypium barbadense during cottonseed development. The prolonged expression period and increased expression abundance of oil-related genes are the main reasons for producing high seed oil content (SOC) in G. barbadense, which manifested as the bias of homeologous gene expression in Dt-subgenome after 25 day postanthesis (DPA). The dynamic expression profile showed that SAD6 and FATA are more important for oil biosynthesis in G. barbadense than that in G. hirsutum. Three key transcription factors, WRI1, NF-YB6 and DPBF2, showed their elite roles in regulating seed oil in cotton. We observed that sequence variations in the promoter region of BCCP2 genes might contribute to its divergence in expression level between the two species. Based on the quantitative trait loci (QTL) information of the seed oil content and utilizing additional G. barbadense introgression lines (ILs), we propose 21 candidate genes on the basis of their differential expression level, of which the GbSWEET and the GbACBP6 showed the potential functional to improve the oil content. Taken together, studying the different expression of oil-related genes and their genetic regulation mechanisms between G. hirsutum and G. barbadense provide new insights to understanding the mechanism of fatty acid biosynthesis network and fatty acid genetic improving breeding in cotton.

Transcriptome Analysis of Seed Weight Plasticity in Brassica napus

A critical barrier to improving crop yield is the trade-off between seed weight (SW) and seed number (SN), which has been commonly reported in several crops, including Brassica napus. Despite the agronomic relevance of this issue, the molecular factors involved in the interaction between SW and SN are largely unknown in crops. In this work, we performed a detailed transcriptomic analysis of 48 seed samples obtained from two rapeseed spring genotypes subjected to different source-sink (S-S) ratios in order to examine the relationship between SW and SN under different field conditions. A multifactorial analysis of the RNA-seq data was used to identify a group of 1014 genes exclusively regulated by the S-S ratio. We found that a reduction in the S-S ratio during seed filling induces the expression of genes involved in sucrose transport, seed weight, and stress responses. Moreover, we identified five co-expression modules that are positively correlated with SW and negatively correlated with SN. Interestingly, one of these modules was significantly enriched in transcription factors (TFs). Furthermore, our network analysis predicted several NAC TFs as major hubs underlying SW and SN compensation. Taken together, our study provides novel insights into the molecular factors associated with the SW-SN relationship in rapeseed and identifies TFs as potential targets when improving crop yield.

Temporal transcriptome profiling of developing seeds reveals candidate genes involved in oil accumulation in safflower (Carthamus tinctorius L.)

BACKGROUND

The investigation of molecular mechanisms involved in lipid metabolism plays a critical role for the genetic engineering of safflower (Carthamus tinctorius L.) to increase the oil accumulation level or to change the oil composition. Although transcript sequences are currently available for the leaves and flowers of safflower, a wide range scan of temporal transcripts at different stages of seed development has not been conducted for safflower.

RESULTS

In this study, temporal transcriptome sequencing was executed at 10, 14, 18, and 22 days after flowering (DAF) to uncover the molecular networks concerned in the biosynthesis of unsaturated fatty acids (USFAs). The results revealed that the biosynthesis of fatty acids is a dominant cellular process from 10 to 14 DAF, while degradation mainly happens after 18 DAF. Significant expression changes of two genes, stearoyl-[acyl-carrier-protein] 9-desaturase gene (SAD) from 10 to 14 DAF and oleate desaturase (FAD2-1) from 14 to 18 DAF, were detected at the transcriptomic levels, and the temporal expression patterns revealed by the transcriptomic analysis were confirmed using quantitative real-time PCR experiments. In addition, 13 candidate transcription factors (TFs) involved in regulating the expression level of the FAD2-1 gene were identified.

CONCLUSIONS

These results create a link between fatty acid biosynthesis and gene expression at different developmental stages of the seeds, provide insight into the underlying lipid metabolism, and meanwhile lay an important foundation for the genetic engineering of safflower varieties. We have identified novel candidate genes, including TFs, that are worthy of further exploration.

Comparative transcriptome and flavonoids components analysis reveal the structural genes responsible for the yellow seed coat color of Brassica rapa L

Background

Seed coat color is an important horticultural trait in Brassica crops, which is divided into two categories: brown/black and yellow. Seeds with yellow seed coat color have higher oil quality, higher protein content and lower fiber content. Yellow seed coat color is therefore considered a desirable trait in hybrid breeding of Brassica rapa, Brassica juncea and Brassica napus.

Methods

Comprehensive analysis of the abundance transcripts for seed coat color at three development stages by RNA-sequencing (RNA-seq) and corresponding flavonoids compounds by liquid chromatography-tandem mass spectrometry (LC-MS/MS) were carried out in B. rapa.

Results

We identified 41,286 unigenes with 4,989 differentially expressed genes between brown seeds (B147) and yellow seeds (B80) at the same development stage. Kyoto Encyclopedia of Genes and Genomes enrichment analysis identified 19 unigenes associated with the phenylpropanoid, flavonoid, flavone and flavonol biosynthetic pathways as involved in seed coat color formation. Interestingly, expression levels of early biosynthetic genes (BrCHS, BrCHI, BrF3H, BrF3’H and BrFLS) in the flavonoid biosynthetic pathway were down-regulated while late biosynthetic genes (BrDFR, BrLDOX and BrBAN) were hardly or not expressed in seeds of B80. At the same time, BrTT8 and BrMYB5 were down-regulated in B80. Results of LC-MS also showed that epicatechin was not detected in seeds of B80. We validated the accuracy of our RNA-seq data by RT-qPCR of nine critical genes. Epicatechin was not detected in seeds of B80 by LC-MS/MS.

Conclusions

The expression levels of flavonoid biosynthetic pathway genes and the relative content of flavonoid biosynthetic pathway metabolites clearly explained yellow seed color formation in B. rapa. This study provides a foundation for further research on the molecular mechanism of seed coat color formation.

Combining quantitative trait locus and co-expression analysis allowed identification of new candidates for oil accumulation in rapeseed

In crops there are quantitative trait loci (QTLs) in which some of the causal quantitative trait genes (QTGs) have not been functionally characterized even in the model plant Arabidopsis. We propose an approach to delineate QTGs in rapeseed by coordinating expression of genes located within QTLs and known orthologs related to traits from Arabidopsis. Using this method in developing siliques 15 d after pollination in 71 lines of rapeseed, we established an acyl-lipid metabolism co-expression network with 21 modules composed of 270 known acyl-lipid genes and 3503 new genes. The core module harbored 76 known genes involved in fatty acid and triacylglycerol biosynthesis and 671 new genes involved in sucrose transport, carbon metabolism, amino acid metabolism, seed storage protein processes, seed maturation, and phytohormone metabolism. Moreover, the core module closely associated with the modules of photosynthesis and carbon metabolism. From the co-expression network, we selected 12 hub genes to identify their putative Arabidopsis orthologs. These putative orthologs were functionally analysed using Arabidopsis knockout and overexpression lines. Four knockout mutants exhibited lower seed oil content, while the seed oil content in 10 overexpression lines was significantly increased. Therefore, combining gene co-expression network analysis and QTL mapping, this study provides new insights into the detection of QTGs and into acyl-lipid metabolism in rapeseed.

Elevating seed oil content in a polyploid crop by induced mutations in SEED FATTY ACID REDUCER genes

Plant-based oils are valuable agricultural products, and seed oil content (SOC) is the major yield component in oil crops. Increasing SOC has been successfully targeted through the selection and genetic modification of oil biosynthesis. The SOC in rapeseed declined during the seed maturation and eventually caused the final accumulated seed oil quantity. However, genes involved in oil degradation during seed maturity are not deeply studied so far. We performed a candidate gene association study using a worldwide collection of rapeseed germplasm. We identified SEED FATTY ACID REDUCER (SFAR) genes, which had a significant effect on SOC and fatty acid (FA) composition. SFAR genes belong to the GDSL lipases, and GDSL lipases have a broad range of functions in plants. After quantification of gene expression using RNA-seq and quantitative PCR, we used targeted (CRISPR-Cas mediated) and random (chemical) mutagenesis to modify turnover rates of seed oil in winter rapeseed. For the first time, we demonstrate significant increase of SOC in a crop after knocking out members of the BnSFAR4 and BnSFAR5 gene families without pleiotropic effects on seed germination, vigour and oil mobilization. Our results offer new perspectives for improving oil yield by targeted mutagenesis.

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.

Genome-Wide Identification and Characterization of FBA Gene Family in Polyploid Crop Brassica napus

Fructose-1,6-bisphosphate aldolase (FBA) is a versatile metabolic enzyme involved in multiple important processes of glycolysis, gluconeogenesis, and Calvin cycle. Despite its significance in plant biology, the identity of this gene family in oil crops is lacking. Here, we performed genome-wide identification and characterization of FBAs in an allotetraploid species, oilseed rape Brassica napus. Twenty-two BnaFBA genes were identified and divided into two groups based on integrative analyses of functional domains, phylogenetic relationships, and gene structures. Twelve and ten B. napus FBAs (BnaFBAs) were predicted to be localized in the chloroplast and cytoplasm, respectively. Notably, synteny analysis revealed that Brassica-specific triplication contributed to the expansion of the BnaFBA gene family during the evolution of B. napus. Various cis-acting regulatory elements pertinent to abiotic and biotic stresses, as well as phytohormone responses, were detected. Intriguingly, each of the BnaFBA genes exhibited distinct sequence polymorphisms. Among them, six contained signatures of selection, likely having experienced breeding selection during adaptation and domestication. Importantly, BnaFBAs showed diverse expression patterns at different developmental stages and were preferentially highly expressed in photosynthetic tissues. Our data thus provided the foundation for further elucidating the functional roles of individual BnaFBA and also potential targets for engineering to improve photosynthetic productivity in B. napus.

Overexpression of Acyl-ACP Thioesterases, CpFatB4 and CpFatB5, Induce Distinct Gene Expression Reprogramming in Developing Seeds of Brassica napus

We examined the substrate preference of Cuphea paucipetala acyl-ACP thioesterases, CpFatB4 and CpFatB5, and gene expression changes associated with the modification of lipid composition in the seed, using Brassica napus transgenic plants overexpressing CpFatB4 or CpFatB5 under the control of a seed-specific promoter. CpFatB4 seeds contained a higher level of total saturated fatty acid (FA) content, with 4.3 times increase in 16:0 palmitic acid, whereas CpFatB5 seeds showed approximately 3% accumulation of 10:0 and 12:0 medium-chain FAs, and a small increase in other saturated FAs, resulting in higher levels of total saturated FAs. RNA-Seq analysis using entire developing pods at 8, 25, and 45 days after flowering (DAF) showed up-regulation of genes for β-ketoacyl-acyl carrier protein synthase I/II, stearoyl-ACP desaturase, oleate desaturase, and linoleate desaturase, which could increase unsaturated FAs and possibly compensate for the increase in 16:0 palmitic acid at 45 DAF in CpFatB4 transgenic plants. In CpFatB5 transgenic plants, many putative chloroplast- or mitochondria-encoded genes were identified as differentially expressed. Our results report comprehensive gene expression changes induced by alterations of seed FA composition and reveal potential targets for further genetic modifications.

Comparative Transcriptome Analysis of Developing Seeds and Silique Wall Reveals Dynamic Transcription Networks for Effective Oil Production in Brassica napus L

Vegetable oil is an essential constituent of the human diet and renewable raw material for industrial applications. Enhancing oil production by increasing seed oil content in oil crops is the most viable, environmentally friendly, and sustainable approach to meet the continuous demand for the supply of vegetable oil globally. An in-depth understanding of the gene networks involved in oil biosynthesis during seed development is a prerequisite for breeding high-oil-content varieties. Rapeseed (Brassica napus) is one of the most important oil crops cultivated on multiple continents, contributing more than 15% of the world’s edible oil supply. To understand the phasic nature of oil biosynthesis and the dynamic regulation of key pathways for effective oil accumulation in B. napus, comparative transcriptomic profiling was performed with developing seeds and silique wall (SW) tissues of two contrasting inbred lines with ~13% difference in seed oil content. Differentially expressed genes (DEGs) between high- and low-oil content lines were identified across six key developmental stages, and gene enrichment analysis revealed that genes related to photosynthesis, metabolism, carbohydrates, lipids, phytohormones, transporters, and triacylglycerol and fatty acid synthesis tended to be upregulated in the high-oil-content line. Differentially regulated DEG patterns were revealed for the control of metabolite and photosynthate production in SW and oil biosynthesis and accumulation in seeds. Quantitative assays of carbohydrates and hormones during seed development together with gene expression profiling of relevant pathways revealed their fundamental effects on effective oil accumulation. Our results thus provide insights into the molecular basis of high seed oil content (SOC) and a new direction for developing high-SOC rapeseed and other oil crops.

Integrated analysis of transcriptomic and proteomic data from tree peony (P. ostii) seeds reveals key developmental stages and candidate genes related to oil biosynthesis and fatty acid metabolism

Tree peony (Paeonia section Moutan DC.) seeds are an excellent source of beneficial natural compounds that promote health, and they contain high levels of alpha-linolenic acid (ALA). In recent years, tree peony has been emerging as an oil crop. Therefore, combined analysis of the transcriptome and proteome of tree peony (P. ostii) seeds at 25, 32, 39, 53, 67, 81, 88, 95, and 109 days after pollination (DAP) was conducted to better understand the transcriptional and translational regulation of seed development and oil biosynthesis. A total of 38,482 unigenes and 2841 proteins were identified. A total of 26,912 differentially expressed genes (DEGs) and 592 differentially expressed proteins (DEPs) were clustered into three groups corresponding to the rapid growth, seed inclusion enrichment and conversion, and late dehydration and mature stages of seed development. Fifteen lipid metabolism pathways were identified at both the transcriptome and proteome levels. Pathway enrichment analysis revealed that a period of rapid fatty acid biosynthesis occurred at 53-88 DAP. Furthermore, 211 genes and 35 proteins associated with the fatty acid metabolism pathway, 63 genes and 11 proteins associated with the biosynthesis of unsaturated fatty acids (UFAs), and 115 genes and 24 proteins associated with ALA metabolism were identified. Phylogenetic analysis revealed that 16 putative fatty acid desaturase (FAD)-encoding genes clustered into four FAD groups, eight of which exhibited the highest expression at 53 DAP, suggesting that they play an important role in ALA accumulation. RT-qPCR analysis indicated that the temporal expression patterns of oil biosynthesis genes were largely similar to the RNA-seq results. The expression patterns of fatty acid metabolism- and seed development-related proteins determined by MRM were also highly consistent with the results obtained in the proteomic analysis. Correlation analysis indicated significant differences in the number and abundance of DEGs and DEPs but a high level of consistency in expression patterns and metabolic pathways. The results of the present study represent the first combined transcriptomic and proteomic analysis of tree peony seeds and provide insight into tree peony seed development and oil accumulation.

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.

Legume Cytosolic and Plastid Acetyl-Coenzyme-A Carboxylase Genes Differ by Evolutionary Patterns and Selection Pressure Schemes Acting before and after Whole-Genome Duplications

Acetyl-coenzyme A carboxylase (ACCase, E.C.6.4.1.2) catalyzes acetyl-coenzyme A carboxylation to malonyl coenzyme A. Plants possess two distinct ACCases differing by cellular compartment and function. Plastid ACCase contributes to de novo fatty acid synthesis, whereas cytosolic enzyme to the synthesis of very long chain fatty acids, phytoalexins, flavonoids, and anthocyanins. The narrow leafed lupin (Lupinus angustifolius L.) represents legumes, a plant family which evolved by whole-genome duplications (WGDs). The study aimed on the contribution of these WGDs to the multiplication of ACCase genes and their further evolutionary patterns. The molecular approach involved bacterial artificial chromosome (BAC) library screening, fluorescent in situ hybridization, linkage mapping, and BAC sequencing. In silico analysis encompassed sequence annotation, comparative mapping, selection pressure calculation, phylogenetic inference, and gene expression profiling. Among sequenced legumes, the highest number of ACCase genes was identified in lupin and soybean. The most abundant plastid ACCase subunit genes were accB. ACCase genes in legumes evolved by WGDs, evidenced by shared synteny and Bayesian phylogenetic inference. Transcriptional activity of almost all copies was confirmed. Gene duplicates were conserved by strong purifying selection, however, positive selection occurred in Arachis (accB2) and Lupinus (accC) lineages, putatively predating the WGD event(s). Early duplicated accA and accB genes underwent transcriptional sub-functionalization.

Gene network of oil accumulation reveals expression profiles in developing embryos and fatty acid composition in Upland cotton

Cottonseed oil accumulated dramatically from 20 days post-anthesis (DPA) to 30 DPA in Upland cotton (Gossypium hirsutum L.). To reveal the gene network of oil accumulation and fatty acid composition in developing embryos, embryos at 10, 20 and 30 DPA in cottonseed were sampled and used for transcriptome sequencing (RNA-Seq). In total, 8629, 7891, and 12,555 differentially expressed genes (DEGs) were identified in the comparison sets of '20 DPA vs 10 DPA', '30 DPA vs 20 DPA', and '30 DPA vs 10 DPA', respectively. The gene network highlighted the dynamic expression profiles of oil accumulation in fatty acid (FA) synthesis, FA desaturation, and triacylglycerol (TAG) biosynthesis. WRI1 and NF-YB6 were suggested elite transcription factors in regulating lipid metabolism. Compared with the gene expression levels in developing seeds, GhPDAT was highly expressed and might play a more important role than GhDGAT in transforming diacylglycerol to TAG in cotton. Expression patterns of 12 FA-biosynthesis-related genes were validated by quantitative real-time PCR (qRT-PCR) method. To reveal the reason for the high content of linoleic acid (C18:2) in cottonseed oil, we carried out a comparative analysis of gene expression levels in Upland cotton, rapeseed (Brassica napus), and oleaster (Olea europaea). Compared with in rapeseed and oleaster, GhFAD2 genes were up-regulated and GhFAD3 genes down-regulated in cottonseed, taking into account the relative high amount of C18:2 but low content of linolenic acid (C18:3) in Upland cotton. The present study offers new information to interpret the mechanism of the FA biosynthesis network and to alter FA composition in cotton breeding projects.

Systems biology of seeds: deciphering the molecular mechanisms of seed storage, dormancy and onset of germination

Seeds are heterogeneous storage reserves with wide array of storage compounds that include various soluble carbohydrates, starch polymer, storage proteins and lipids. These stored reserves comprise 70% of the world's caloric intake in the form of food and animal feed produced through sustainable agriculture, which contributes to food and nutritional security. Seed systems biology remains an enigmatic subject in understanding seed storage processes, maturation and pre-germinative metabolism. The reviews and research articles covered in this special issue of Plant Cell Reports highlight recent advances made in the area of seed biology that cover various systems biology applications such as gene regulatory networks, metabolomics, epigenetics and the role of micro-RNA in seed development.