Genome-wide identification and comparative analysis of diacylglycerol kinase (DGK) gene family and their expression profiling in Brassica napus under abiotic stress.

BACKGROUND: Diacylglycerol kinases (DGKs) are signaling enzymes that play pivotal roles in response to abiotic and biotic stresses by phosphorylating diacylglycerol (DAG) to form phosphatidic acid (PA). However, no comprehensive analysis of the DGK gene family had previously been reported in B. napus and its diploid progenitors (B. rapa and B. oleracea). RESULTS: In present study, we identified 21, 10, and 11 DGK genes from B. napus, B. rapa, and B. oleracea, respectively, which all contained conserved catalytic domain and were further divided into three clusters. Molecular evolutionary analysis showed that speciation and whole-genome triplication (WGT) was critical for the divergence of duplicated DGK genes. RNA-seq transcriptome data revealed that, with the exception of BnaDGK4 and BnaDGK6, BnaDGK genes have divergent expression patterns in most tissues. Furthermore, some DGK genes were upregulated or downregulated in response to hormone treatment and metal ion (arsenic and cadmium) stress. Quantitative real-time PCR analysis revealed that different BnaDGK genes contribute to seed oil content. CONCLUSIONS: Together, our results indicate that DGK genes have diverse roles in plant growth and development, hormone response, and metal ion stress, and in determining seed oil content, and lay a foundation for further elucidating the roles of DGKs in Brassica species.

Genome-Wide Identification and Comparative Expression Profile Analysis of the Long-Chain Acyl-CoA synthetase (LACS) Gene Family in Two Different Oil Content Cultivars of Brassica napus.

Long-chain acyl-CoA synthetase (LACS) is one of the key enzymes involved in fatty acid metabolism, including phospholipid biosynthesis, triacylglycerol (TAG) biosynthesis, and fatty acid ß-oxidation in plants. However, the characterization of LACSs family in seed oil biosynthesis of Brassica napus (B. napus) remains unknown. In the present study, we performed a comprehensive genome-wide analysis of this gene family in B. napus, and 34 B. napus LACS genes (BnaLACSs) were identified. Phylogenetic analysis classified the BnaLACS proteins into four groups (A, B, C, and D), which were supported by highly conserved gene structures and consensus motifs. RNA-Sequencing (RNA-Seq) and qRT-PCR combined analysis revealed that 18 BnaLACSs (BnaLACS1-2, 1-3, 1-4, 1-9, 1-10, 2-1, 2-2, 4-1, 4-2, 6-1, 6-2, 6-4, 7-1, 7-2, 8-1, 8-2, 9-3, and 9-4) were highly expressed in developmental seeds. Comparative expression analysis between extremely high oil content (P1-HO) and low oil content (P2-LO) B. napus cultivars revealed that BnaLACS6-4, BnaLACS9-3, and BnaLACS9-4 may be involved in fatty acid synthesis in chloroplast, and BnaLACS1-10 and 4-1 may play a vital role in lipid biosynthesis in B. napus, which is important for further seed oil accumulation in oilseed rape. The present study provides important information for functional characterization of BnaLACSs in seed oil metabolism in B. napus.

Transcriptome and Small RNA Sequencing Reveal the Mechanisms Regulating Harvest Index in Brassica napus.

Harvest index (HI), the ratio of harvested seed weight to total aboveground biomass weight, is an economically critical value reflecting the convergence of complex agronomic traits. HI values in rapeseed (Brassica napus) remain much lower than in other major crops, and the underlying regulatory network is largely unknown. In this study, we performed mRNA and small RNA sequencing to reveal the mechanisms shaping HI in B. napus during the seed-filling stage. A total of 8,410 differentially expressed genes (DEGs) between high-HI and low-HI accessions in four tissues (silique pericarp, seed, leaves, and stem) were identified. Combining with co-expression network, 72 gene modules were identified, and a key gene BnaSTY46 was found to participate in retarded establishment of photosynthetic capacity to influence HI. Further research found that the genes involved in circadian rhythms and response to stimulus may play important roles in HI and that their transcript levels were modulated by differentially expressed microRNAs (DEMs), and we identified 903 microRNAs (miRNAs), including 46 known miRNAs and 857 novel miRNAs. Furthermore, transporter activity-related genes were critical to enhancing HI in good cultivation environments. Of 903 miRNAs, we found that the bna-miR396-Bna.A06SRp34a/Bna.A01EMB3119 pair may control the seed development and the accumulation of storage compounds, thus contributing to higher HI. Our findings uncovered the underlying complex regulatory network behind HI and offer potential approaches to rapeseed improvement.

Identification of candidate genes regulating seed oil content by QTL mapping and transcriptome sequencing in Brassica napus.

Increasing oil production is a major goal in rapeseed (Brassica napus) molecular breeding programs. Identifying seed oil content (SOC)-related candidate genes is an important step towards achieving this goal. We performed quantitative trait locus (QTL) mapping of SOC in B. napus using a high-density SNP genetic map constructed from recombinant inbred lines and the Illumina InfiniumTM 60K SNP array. A total of 26 QTLs were detected in three years on A01, A03, A05, A06, A09, C01, C03 and C05, which accounted for 3.69%~18.47% of the phenotypic variation in SOC. Of these, 13 QTLs are reported here for the first time. 1713 candidate genes in the 26 QTLs confidence interval were obtained. We then identified differentially expressed genes (DEGs) between the high- and low-SOC accessions, to narrow down our focus to 21 candidate genes (Y1-Y21) related to SOC, and we will focus on 11 (Y1-Y11) candidate genes that contribute to the formation of high-SOC. In addition to providing insight into the genetic basis of SOC in B. napus, the loci identified and candidate genes in this study can be used in molecular breeding strategies to increase SOC in this important seed crop.

The Brassica napus fatty acid exporter FAX1-1 contributes to biological yield, seed oil content, and oil quality.

BACKGROUND: In the oilseed crop Brassica napus (rapeseed), various metabolic processes influence seed oil content, oil quality, and biological yield. However, the role of plastid membrane proteins in these traits has not been explored. RESULTS: Our genome-wide association study (GWAS) of 520 B. napus accessions identified the chloroplast membrane protein-localized FATTY ACID EXPORTER 1-1 (FAX1-1) as a candidate associated with biological yield. Seed transcript levels of BnaFAX1-1 were higher in a cultivar with high seed oil content relative to a low-oil cultivar. BnaFAX1-1 was localized to the plastid envelope. When expressed in Arabidopsis thaliana, BnaFAX1-1 enhanced biological yield (total plant dry matter), seed yield and seed oil content per plant. Likewise, in the field, B. napus BnaFAX1-1 overexpression lines (BnaFAX1-1-OE) displayed significantly enhanced biological yield, seed yield, and seed oil content compared with the wild type. BnaFAX1-1 overexpression also up-regulated gibberellic acid 4 (GA4) biosynthesis, which may contribute to biological yield improvement. Furthermore, oleic acid (C18:1) significantly increased in BnaFAX1-1 overexpression seeds. CONCLUSION: Our results indicated that the putative fatty acid exporter BnaFAX1-1 may simultaneously improve seed oil content, oil quality and biological yield in B. napus, providing new approaches for future molecular breeding.

Genome-Wide Identification and Expression Profiling of Monosaccharide Transporter Genes Associated with High Harvest Index Values in Rapeseed (Brassica napus L.).

Sugars are important throughout a plant's lifecycle. Monosaccharide transporters (MST) are essential sugar transporters that have been identified in many plants, but little is known about the evolution or functions of MST genes in rapeseed (Brassica napus). In this study, we identified 175 MST genes in B. napus, 87 in Brassica oleracea, and 83 in Brassica rapa. These genes were separated into the sugar transport protein (STP), polyol transporter (PLT), vacuolar glucose transporter (VGT), tonoplast monosaccharide transporter (TMT), inositol transporter (INT), plastidic glucose transporter (pGlcT), and ERD6-like subfamilies, respectively. Phylogenetic and syntenic analysis indicated that gene redundancy and gene elimination have commonly occurred in Brassica species during polyploidization. Changes in exon-intron structures during evolution likely resulted in the differences in coding regions, expression patterns, and functions seen among BnMST genes. In total, 31 differentially expressed genes (DEGs) were identified through RNA-seq among materials with high and low harvest index (HI) values, which were divided into two categories based on the qRT-PCR results, expressed more highly in source or sink organs. We finally identified four genes, including BnSTP5, BnSTP13, BnPLT5, and BnERD6-like14, which might be involved in monosaccharide uptake or unloading and further affect the HI of rapeseed. These findings provide fundamental information about MST genes in Brassica and reveal the importance of BnMST genes to high HI in B. napus.

Identification of candidate genes controlling oil content by combination of genome-wide association and transcriptome analysis in the oilseed crop Brassica napus.

BACKGROUND: Increasing seed oil content is one of the most important targets for rapeseed (Brassica napus) breeding. However, genetic mechanisms of mature seed oil content in Brassica napus (B. napus) remain little known. To identify oil content-related genes, a genome-wide association study (GWAS) was performed using 588 accessions. RESULTS: High-throughput genome resequencing resulted in 385,692 high-quality single nucleotide polymorphism (SNPs) with a minor allele frequency (MAF) > 0.05. We identified 17 loci that were significantly associated with seed oil content, among which 12 SNPs were distributed on the A3 (11 loci) and A1 (one loci) chromosomes, and five novel significant SNPs on the C5 (one loci) and C7 (four loci) chromosomes, respectively. Subsequently, we characterized differentially expressed genes (DEGs) between the seeds and silique pericarps on main florescences and primary branches of extremely high- and low-oil content accessions (HO and LO). A total of 64 lipid metabolism-related DEGs were identified, 14 of which are involved in triacylglycerols (TAGs) biosynthesis and assembly. Additionally, we analyzed differences in transcription levels of key genes involved in de novo fatty acid biosynthesis in the plastid, TAGs assembly and lipid droplet packaging in the endoplasmic reticulum (ER) between high- and low-oil content B. napus accessions. CONCLUSIONS: The combination of GWAS and transcriptome analyses revealed seven candidate genes located within the confidence intervals of significant SNPs. Current findings provide valuable information for facilitating marker-based breeding for higher seed oil content in B. napus.

Characterization of Fatty Acid Exporters involved in fatty acid transport for oil accumulation in the green alga Chlamydomonas reinhardtii.

BACKGROUND: In the past few decades, microalgae biofuel has become one of the most interesting sources of renewable energy. However, the higher cost of microalgae biofuel compared to that of petroleum prevented microalgae biofuel production. Therefore, the research on increasing lipid productivity from microalgae becomes more important. The lipid production source, triacylglycerol biosynthesis in microalgae requires short chain fatty acids as substrates, which are synthesized in chloroplasts. However, the transport mechanism of fatty acids from microalgae chloroplasts to cytosol remains unknown. RESULTS: cDNAs from two homologs of the Arabidopsis fatty acid exporter 1 (FAX1) were cloned from Chlamydomonas reinhardtii and were named crfax1 and crfax2. Both CrFAXs were involved in fatty acid transport, and their substrates were mainly C16 and C18 fatty acids. Overexpression of both CrFAXs increased the accumulation of the total lipid content in algae cells, and the fatty acid compositions were changed under normal TAP or nitrogen deprivation conditions. Overexpression of both CrFAXs also increased the chlorophyll content. The MGDG content was decreased but the TAG, DAG, DGDG and other lipid contents were increased in CrFAXs overexpression strains. CONCLUSION: These results reveal that CrFAX1 and CrFAX2 were involved in mediating fatty acid export for lipids biosynthesis in C. reinhardtii. In addition, overexpression of both CrFAXs obviously increased the intracellular lipid content, especially the triacylglycerol content in microalgae, which provides a potential technology for the production of more biofuels using microalgae.

Whole-genome resequencing reveals Brassica napus origin and genetic loci involved in its improvement.

Brassica napus (2n = 4x = 38, AACC) is an important allopolyploid crop derived from interspecific crosses between Brassica rapa (2n = 2x = 20, AA) and Brassica oleracea (2n = 2x = 18, CC). However, no truly wild B. napus populations are known; its origin and improvement processes remain unclear. Here, we resequence 588 B. napus accessions. We uncover that the A subgenome may evolve from the ancestor of European turnip and the C subgenome may evolve from the common ancestor of kohlrabi, cauliflower, broccoli, and Chinese kale. Additionally, winter oilseed may be the original form of B. napus. Subgenome-specific selection of defense-response genes has contributed to environmental adaptation after formation of the species, whereas asymmetrical subgenomic selection has led to ecotype change. By integrating genome-wide association studies, selection signals, and transcriptome analyses, we identify genes associated with improved stress tolerance, oil content, seed quality, and ecotype improvement. They are candidates for further functional characterization and genetic improvement of B. napus.

Association Mapping Analysis of Fatty Acid Content in Different Ecotypic Rapeseed Using mrMLM.

Brassica napus L. is a widely cultivated oil crop and provides important resources of edible vegetable oil, and its quality is determined by fatty acid composition and content. To explain the genetic basis and identify more minor loci for fatty acid content, the multi-locus random-SNP-effect mixed linear model (mrMLM) was used to identify genomic regions associated with fatty acid content in a genetically diverse population of 435 rapeseed accessions, including 77 winter-type, 55 spring-type, and 303 semi-winter-type accessions grown in different environments. A total of 149 quantitative trait nucleotides (QTNs) were found to be associated with fatty acid content and composition, including 34 QTNs that overlapped with the previously reported loci, and 115 novel QTNs. Of these, 35 novel QTNs, located on chromosome A01, A02, A03, A05, A06, A09, A10, and C02, respectively, were repeatedly detected across different environments. Subsequently, we annotated 95 putative candidate genes by BlastP analysis using sequences from Arabidopsis thaliana homologs of the identified regions. The candidate genes included 34 environmentally-insensitive genes (e.g., CER4, DGK2, KCS17, KCS18, MYB4, and TT16) and 61 environment-sensitive genes (e.g., FAB1, FAD6, FAD7, KCR1, KCS9, KCS12, and TT1) as well as genes invloved in the fatty acid biosynthesis. Among these, BnaA08g08280D and BnaC03g60080D differed in genomic sequence between the high- and low-oleic acid lines, and might thus be the novel alleles regulating oleic acid content. Furthermore, RT-qPCR analysis of these genes showed differential expression levels during seed development. Our results highlight the practical and scientific value of mrMLM or QTN detection and the accuracy of linking specific QTNs to fatty acid content, and suggest a useful strategy to improve the fatty acid content of B. napus seeds by molecular marker-assisted breeding.

Genome-Wide Association and Transcriptome Analyses Reveal Candidate Genes Underlying Yield-determining Traits in Brassica napus.

Yield is one of the most important yet complex crop traits. To improve our understanding of the genetic basis of yield establishment, and to identify candidate genes responsible for yield improvement in Brassica napus, we performed genome-wide association studies (GWAS) for seven yield-determining traits [main inflorescence pod number (MIPN), branch pod number (BPN), pod number per plant (PNP), seed number per pod (SPP), thousand seed weight, main inflorescence yield (MIY), and branch yield], using data from 520 diverse B. napus accessions from two different yield environments. In total, we detected 128 significant single nucleotide polymorphisms (SNPs), 93 of which were revealed as novel by integrative analysis. A combination of GWAS and transcriptome sequencing on 21 haplotype blocks from samples pooled by four extremely high-yielding or low-yielding accessions revealed the differential expression of 14 crucial candiate genes (such as Bna.MYB83, Bna.SPL5, and Bna.ROP3) associated with multiple traits or containing multiple SNPs associated with the same trait. Functional annotation and expression pattern analyses further demonstrated that these 14 candiate genes might be important in developmental processes and biomass accumulation, thus affecting the yield establishment of B. napus. These results provide valuable information for understanding the genetic mechanisms underlying the establishment of high yield in B. napus, and lay the foundation for developing high-yielding B. napus varieties.

A combination of genome-wide association and transcriptome analysis reveals candidate genes controlling harvest index-related traits in Brassica napus.

Harvest index (HI), the ratio of seed mass to total biomass of the aboveground plant parts, is an important trait for harvestable yield of crops. Unfortunately, HI of Brassica napus is lower than that of other economically important crops. To identify candidate genes associated with high HI, a genome-wide association study of HI and four HI-related traits was conducted with 520 B. napus accessions cultivated in both Yunnan and Chongqing. We detected 294 single nucleotide polymorphisms significantly associated with the abovementioned traits, including 79 SNPs that affected two or more traits. Differentially expressed genes between extremely high- and low-HI accessions were identified in 8 tissues at two cultivated regions. Combination of linkage disequilibrium and transcriptome analyses revealed 33 functional candidate genes located within the confidence intervals of significant SNPs associated with more than one trait, such as SHOOT GRAVITROPISM 5 (Bna.SGR5), ATP-CITRATE LYASE A-3 (Bna.ACLA-3) and CAROTENOID CLEAVAGE DIOXYGENASE 1 (Bna.CCD1), their orthologs in the Arabidopsis thaliana have been shown to play key roles in photosynthesis, inflorescence, and silique development. Our results provide insight into the molecular mechanisms underlying establishment of high-HI B. napus and lay a foundation for characterization of candidate genes aimed at developing high-HI B. napus varieties.