Deciphering the transcriptional regulatory networks that control size, color, and oil content in Brassica rapa seeds

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.

Comparative transcriptome and co-expression network analysis revealed the genes associated with senescence and polygalacturonase activity involved in pod shattering of rapeseed

BACKGROUND

The pod shattering (PS) trait negatively affects the crop yield in rapeseed especially under dry conditions. To better understand the trait and cultivate higher resistance varieties, it's necessary to identify key genes and unravel the PS mechanism thoroughly.

RESULTS

In this study, we conducted a comparative transcriptome analysis between two materials significantly different in silique shatter resistance lignin deposition and polygalacturonase (PG) activity. Here, we identified 10,973 differentially expressed genes at six pod developmental stages. We found that the late pod development stages might be crucial in preparing the pods for upcoming shattering events. GO enrichment results from K-means clustering and weighed gene correlation network analysis (WGCNA) both revealed senescence-associated genes play an important role in PS. Two hub genes Bna.A05ABI5 and Bna.C03ERF/AP2-3 were selected from the MEyellow module, which possibly regulate the PS through senescence-related mechanisms. Further investigation found that senescence-associated transcription factor Bna.A05ABI5 upregulated the expression of SAG2 and ERF/AP2 to control the shattering process. In addition, the upregulation of Bna.C03ERF/AP2-3 is possibly involved in the transcription of downstream SHP1/2 and LEA proteins to trigger the shattering mechanism. We also analyzed the PS marker genes and found Bna.C07SHP1/2 and Bna.PG1/2 were significantly upregulated in susceptible accession. Furthermore, the role of auxin transport by Bna.WAG2 was also observed, which could reduce the PG activity to enhance the PS resistance through the cell wall loosening process.

CONCLUSION

Based on comparative transcriptome evaluation, this study delivers insights into the regulatory mechanism primarily underlying the variation of PS in rapeseed. Taken together, these results provide a better understanding to increase the yield of rapeseed by reducing the PS through better engineered crops.

Multi-omics revolution to promote plant breeding efficiency

Crop production is the primary goal of agricultural activities, which is always taken into consideration. However, global agricultural systems are coming under increasing pressure from the rising food demand of the rapidly growing world population and changing climate. To address these issues, improving high-yield and climate-resilient related-traits in crop breeding is an effective strategy. In recent years, advances in omics techniques, including genomics, transcriptomics, proteomics, and metabolomics, paved the way for accelerating plant/crop breeding to cope with the changing climate and enhance food production. Optimized omics and phenotypic plasticity platform integration, exploited by evolving machine learning algorithms will aid in the development of biological interpretations for complex crop traits. The precise and progressive assembly of desire alleles using precise genome editing approaches and enhanced breeding strategies would enable future crops to excel in combating the changing climates. Furthermore, plant breeding and genetic engineering ensures an exclusive approach to developing nutrient sufficient and climate-resilient crops, the productivity of which can sustainably and adequately meet the world's food, nutrition, and energy needs. This review provides an overview of how the integration of omics approaches could be exploited to select crop varieties with desired traits.

Seed coat colour of Indian mustard [Brassica juncea (L.) Czern. and Coss.] is associated with Bju.TT8 homologs identifiable by targeted functional markers

Seed coat colour is an important trait in Indian mustard. Breeding for seed coat colour needs precise knowledge of mode of inheritance and markers linked to it. The present study was focussed on genetics and development of functional markers for seed coat colour. F₁s (direct and reciprocal) and F₂ populations were developed by crossing two contrasting parents for seed coat colour (DRMRIJ-31, brown seeded and RLC-3, yellow seeded). Phenotypic results have shown that the seed coat colour trait was under the influence of maternal effect and controlled by digenic-duplicate gene action. Further, Bju.TT8 homologs of both parents (DRMRIJ-31 and RLC-3) were cloned and sequenced. Sequencing results of Bju.TT8 homologs revealed that in RLC-3, gene Bju.ATT8 had an insertion of 1279bp in the 7^th exon; whereas, gene Bju.BTT8 had an SNP (C→T) in the 7^th exon. These two mutations were found to be associated with yellow seed coat colour. Using sequence information, functional markers were developed for both Bju.TT8 homologs, validated on F₂ population and were found highly reliable with no recombination between the markers and the phenotype. Further, these markers were subjected to a germplasm assembly of Indian mustard, and their allelic combination for the seed coat colour genes has been elucidated. The comparative genomics of TT8 genes revealed high degree of similarity between and across the Brassica species, and the respective diploid progenitors in tetraploid Brassica species are the possible donors of TT8 homologs. This study will help in the marker-assisted breeding for seed coat colour, and aid in understanding seed coat colour genetics more precisely.

24-Epibrassinolide Promotes Fatty Acid Accumulation and the Expression of Related Genes in Styrax tonkinensis Seeds

Styrax tonkinensis, whose seeds are rich in unsaturated fatty acids (UFAs), is a high oil value tree species, and the seed oil has perfect biodiesel properties. Therefore, the elucidation of the effect of 24-epibrassinolide (EBL) on fatty acid (FA) concentration and the expression of FA biosynthesis-related genes is critical for deeply studying the seed oil in S. tonkinensis. In this study, we aimed to investigate the changing trend of FA concentration and composition and identify candidate genes involved in FA biosynthesis under EBL treatment using transcriptome sequencing and GC-MS. The results showed that 5 μmol/L of EBL (EBL5) boosted the accumulation of FA and had the hugest effect on FA concentration at 70 days after flowering (DAF). A total of 20 FAs were identified; among them, palmitic acid, oleic acid, linoleic acid, and linolenic acid were the main components. In total, 117,904 unigenes were detected, and the average length was 1120 bp. Among them, 1205 unigenes were assigned to ‘lipid translations and metabolism’ in COG categories, while 290 unigenes were assigned to ‘biosynthesis of unsaturated fatty acid’ in KEGG categories. Twelve important genes related to FA biosynthesis were identified, and their expression levels were confirmed by quantitative real-time PCR. KAR, KASIII, and accA, encoding FA biosynthesis-related enzymes, all expressed the highest at 70 DAF, which was coincident with a rapid rise in FA concentration during seed development. FAD2 and FATB conduced to UFA and saturated fatty acids (SFA) accumulation, respectively. EBL5 induced the expression of FA biosynthesis-related genes. The concentration of FA was increased after EBL5 application, and EBL5 also enhanced the enzyme activity by promoting the expression of genes related to FA biosynthesis. Our research could provide a reference for understanding the FA biosynthesis of S. tonkinensis seeds at physiological and molecular levels.

Spatio-temporal transcriptome profiling and subgenome analysis in Brassica napus

Brassica napus is an important oil crop and an allotetraploid species. However, the detailed analysis of gene function and homoeologous gene expression in all tissues at different developmental stages was not explored. In this study, we performed a global transcriptome analysis of 24 vegetative and reproductive tissues at six developmental stages (totally 111 tissues). These samples were clustered into eight groups. The gene functions of silique pericarp were similar to roots, stems and leaves. In particular, glucosinolate metabolic process was associated with root and silique pericarp. Genes involved in protein phosphorylation were often associated with stamen, anther and the early developmental stage of seeds. Transcription factor (TF) genes were more specific than structural genes. A total of 17 100 genes that were preferentially expressed in one tissue (tissue-preferred genes, TPGs), including 889 TFs (5.2%), were identified in the 24 tissues. Some TPGs were identified as hub genes in the co-expression network analysis, and some TPGs in different tissues were involved in different hormone pathways. About 67.0% of the homoeologs showed balanced expression, whereas biased expression of homoeologs was associated with structural divergence. In addition, the spatiotemporal expression of homoeologs was related to the presence of transposable elements (TEs) and regulatory elements (REs); more TEs and fewer REs in the promoters resulted in divergent expression in different tissues. This study provides a valuable transcriptional map for understanding the growth and development of B. napus, for identifying important genes for future crop improvement, and for exploring gene expression patterns in the B. napus.

Multi-omics analysis reveals the mechanism of seed coat color formation in Brassica rapa L

Multi-omics analysis of the transcriptome, metabolome and genome identified major and minor loci and candidate genes for seed coat color and explored the mechanism of flavonoid metabolites biosynthesis in Brassica rapa. Yellow seed trait is considered an agronomically desirable trait with great potential for improving seed quality of Brassica crops. Mechanisms of the yellow seed trait are complex and not well understood. In this study, we performed an integrated metabolome, transcriptome and genome-wide association study (GWAS) on different B. rapa varieties to explore the mechanisms underlying the seed coat color formation. A total of 2,499 differentially expressed genes and 116 differential metabolites between yellow and black seeds with strong association with the flavonoid biosynthesis pathway was identified. In addition, 330 hub genes involved in the seed coat color formation, and the most significantly differential flavonoids biosynthesis were detected based on weighted gene co-expression network analysis. Metabolite GWAS analysis using the contents of 42 flavonoids in developing seeds of 159 B. rapa lines resulted in the identification of 1,626 quantitative trait nucleotides (QTNs) and 37 chromosomal intervals, including one major locus on chromosome A09. A combination of QTNs detection, transcriptome and functional analyses led to the identification of 241 candidate genes that were associated with different flavonoid metabolites. The flavonoid biosynthesis pathway in B. rapa was assembled based on the identified flavonoid metabolites and candidate genes. Furthermore, BrMYB111 members (BraA09g004490.3C and BraA06g034790.3C) involved in the biosynthesis of taxifolin were functionally analyzed in vitro. Our findings lay a foundation and provide a reference for systematically investigating the mechanism of seed coat color in B. rapa and in the other plants.

A major yellow-seed QTL on chromosome A09 significantly increases the oil content and reduces the fiber content of seed in Brassica napus

A major yellow-seed QTL on chromosome A09 significantly increases the oil content and reduces the fiber content of seed in Brassica napus. The yellow-seed trait (YST) has always been a main breeding objective for rapeseed because yellow-seeded B. napus generally contains higher oil contents, fewer pigments and polyphenols and lower fiber content than black-seeded B. napus, although the mechanism controlling this correlation remains unclear. In this study, QTL mapping was implemented for YST based on a KN double haploid population derived from the hybridization of yellow-seeded B. napus N53-2 with a high oil content and black-seeded Ken-C8 with a relatively low oil content. Ten QTLs were identified, including four stable QTLs that could be detected in multiple environments. A major QTL, cqSC-A09, on chromosome A09 was identified by both QTL mapping and BSR-Seq technology, and explained more than 41% of the phenotypic variance. The major QTL cqSC-A09 for YST not only controls the seed color but also affects the oil and fiber contents in seeds. More importantly, the advantageous allele could increase the oil content and reduce the pigment and fiber content at the same time. This is the first QTL reported to control seed color, oil content and fiber content simultaneously with a large effect and has great application value for breeding high oil varieties with high seed quality. Important candidate genes, including BnaA09. JAZ1, BnaA09. GH3.3 and BnaA09. LOX3, were identified for cqSC-A09 by combining sequence variation annotation, expression differences and an interaction network, which lays a foundation for further cloning and breeding applications in the future.

Comparative Analysis of Seed Transcriptome and Coexpression Analysis Reveal Candidate Genes for Enhancing Seed Size/Weight in Brassica juncea

Seed size/weight is a multigenic trait that is governed by complex transcriptional regulatory pathways. An understanding of the genetic basis of seed size is of great interest in the improvement of seed yield and quality in oilseed crops. A global transcriptome analysis was performed at the initial stages of seed development in two lines of Brassica juncea, small-seeded EH-2 and large-seeded PJ. The anatomical analyses revealed significant differences in cell number and cell size in the outer layer of the seed coat between EH-2 and PJ. Pairwise comparisons at each developmental stage identified 5,974 differentially expressed genes (DEGs) between the two lines, of which 954 genes belong to different families of transcription factors. Two modules were found to be significantly correlated with an increased seed size using weighted gene coexpression network analysis. The DEG and coexpression datasets were integrated with the thousand seed weight (Tsw) quantitative trait loci (QTL) mapped earlier in the EPJ (EH-2 × PJ) doubled haploid (DH) population, which identified forty potential key components controlling seed size. The candidate genes included genes regulating the cell cycle, cell wall biogenesis/modification, solute/sugar transport, and hormone signaling. The results provide a valuable resource to widen the current understanding of regulatory mechanisms underlying seed size in B. juncea.

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.

Abscisic Acid Improves Linoleic Acid Accumulation Possibly by Promoting Expression of EgFAD2 and Other Fatty Acid Biosynthesis Genes in Oil Palm Mesocarp

Abscisic acid plays an important role in fruit development. However, the effect of ABA on fatty acid biosynthesis in oil palm is still unknown. In this study, ABA treatments (CK, A1-A4) were applied to oil palm fruit at 16 WAP (weeks after pollination), and fatty acids in the mesocarp at 24 WAP were analyzed by GC-MS. Results showed that linoleic acid content under treatment A2 (20 μM ABA) was significantly higher (slightly increased by 8.33%) than the control. Therefore, mesocarp samples of A2, and the control at 16, 20, and 24 WAP was sampled for RNA-Seq. KEGG pathway enrichment analysis showed that 43 genes were differentially expressed in the fatty acid biosynthesis pathway, of which expression of EgFAD2 (unigene 105050201) under 20 μM ABA treatment was 1.84-fold higher than in the control at 20 WAP. Further sequence analysis found that unigene 105050201 had more ABA-responsive elements (ABRE), complete conserved domains, and a C-terminal signaling motif among two FAD2 copies. Furthermore, WGCNA and correlation analysis showed co-expression of EgFAD2 (unigene 105050201) with transcription factors (TFs) (WRI1, AP2-EREBP, bZIP, bHLH, C2C2-Dof, MYB, NAC, and WRKY), ABA signaling genes (PYR, PP2C, SnRK, and ABI5), and other genes involved in fatty acid biosynthesis (FATA, FATB, LACS, SAD, Oleosins, and so on). These results indicated that ABA treatment promoted the expression of FAD2 and other genes involved in fatty acid biosynthesis, which possibly resulted in the accumulation of linoleic acid. This study will be helpful for understanding the possible mechanisms through which ABA affects fatty acid biosynthesis and their accumulation in the mesocarp of oil palm.

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.

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.

Unravelling the Complex Interplay of Transcription Factors Orchestrating Seed Oil Content in Brassica napus L

Transcription factors (TFs) and their complex interplay are essential for directing specific genetic programs, such as responses to environmental stresses, tissue development, or cell differentiation by regulating gene expression. Knowledge regarding TF-TF cooperations could be promising in gaining insight into the developmental switches between the cultivars of Brassica napus L., namely Zhongshuang11 (ZS11), a double-low accession with high-oil- content, and Zhongyou821 (ZY821), a double-high accession with low-oil-content. In this regard, we analysed a time series RNA-seq data set of seed tissue from both of the cultivars by mainly focusing on the monotonically expressed genes (MEGs). The consideration of the MEGs enables the capturing of multi-stage progression processes that are orchestrated by the cooperative TFs and, thus, facilitates the understanding of the molecular mechanisms determining seed oil content. Our findings show that TF families, such as NAC, MYB, DOF, GATA, and HD-ZIP are highly involved in the seed developmental process. Particularly, their preferential partner choices as well as changes in their gene expression profiles seem to be strongly associated with the differentiation of the oil content between the two cultivars. These findings are essential in enhancing our understanding of the genetic programs in both cultivars and developing novel hypotheses for further experimental studies.