Genome- and transcriptome-wide association studies reveal the genetic basis and the breeding history of seed glucosinolate content in Brassica napus

Exploring silique number in Brassica napus L.: Genetic and molecular advances for improving yield

Silique number is a crucial yield-related trait for the genetic enhancement of rapeseed (Brassica napus L.). The intricate molecular process governing the regulation of silique number involves various factors. Despite advancements in understanding the mechanisms regulating silique number in Arabidopsis (Arabidopsis thaliana) and rice (Oryza sativa), the molecular processes involved in controlling silique number in rapeseed remain largely unexplored. In this review, we identify candidate genes and review the roles of genes and environmental factors in regulating rapeseed silique number. We use genetic regulatory networks for silique number in Arabidopsis and grain number in rice to uncover possible regulatory pathways and molecular mechanisms involved in regulating genes associated with rapeseed silique number. A better understanding of the genetic network regulating silique number in rapeseed will provide a theoretical basis for the genetic improvement of this trait and genetic resources for the molecular breeding of high-yielding rapeseed.

Natural variation in BnaA07.MKK9 confers resistance to Sclerotinia stem rot in oilseed rape

Sclerotinia stem rot (SSR), caused by the necrotrophic fungus Sclerotinia sclerotiorum, is one of the most devastating diseases for several major oil-producing crops. Despite its impact, the genetic basis of SSR resistance in plants remains poorly understood. Here, through a genome-wide association study, we identify a key gene, BnaA07. MKK9, that encodes a mitogen-activated protein kinase kinase that confers SSR resistance in oilseed rape. Our functional analyses reveal that BnaA07.MKK9 interacts with BnaC03.MPK3 and BnaC03.MPK6 and phosphorylates them at the TEY activation motif, triggering a signaling cascade that initiates biosynthesis of ethylene, camalexin, and indole glucosinolates, and promotes accumulation of H2O2 and the hypersensitive response, ultimately conferring resistance. Furthermore, variations in the coding sequence of BnaA07.MKK9 alter its kinase activity and improve SSR resistance by ~30% in cultivars carrying the advantageous haplotype. These findings enhance our understanding of SSR resistance and may help engineer novel diversity for future breeding of oilseed rape.

Functional genomics of Brassica napus: Progresses, challenges, and perspectives

Brassica napus, commonly known as rapeseed or canola, is a major oil crop contributing over 13% to the stable supply of edible vegetable oil worldwide. Identification and understanding the gene functions in the B. napus genome is crucial for genomic breeding. A group of genes controlling agronomic traits have been successfully cloned through functional genomics studies in B. napus. In this review, we present an overview of the progress made in the functional genomics of B. napus, including the availability of germplasm resources, omics databases and cloned functional genes. Based on the current progress, we also highlight the main challenges and perspectives in this field. The advances in the functional genomics of B. napus contribute to a better understanding of the genetic basis underlying the complex agronomic traits in B. napus and will expedite the breeding of high quality, high resistance and high yield in B. napus varieties.

Two-Component System Genes in Brassica napus: Identification, Analysis, and Expression Patterns in Response to Abiotic and Biotic Stresses

The two-component system (TCS), consisting of histidine kinases (HKs), histidine phosphotransfer proteins (HPs) and response regulators (RRs) in eukaryotes, plays pivotal roles in regulating plant growth, development, and responses to environment stimuli. However, the TCS genes were poorly characterized in rapeseed, which is an important tetraploid crop in Brassicaceae. In this work, a total of 182 BnaTCS genes were identified, including 43 HKs, 16 HPs, and 123 RRs, which was more than that in other crops due to segmental duplications during the process of polyploidization. It was significantly different in genetic diversity between the three subfamilies, and some members showed substantial genetic differentiation among the three rapeseed ecotypes. Several hormone- and stress-responsive cis-elements were identified in the putative promoter regions of BnaTCS genes. Furthermore, the expression of BnaTCS genes under abiotic stresses, exogenous phytohormone, and biotic stresses was analyzed, and numerous candidate stress-responsive genes were screened out. Meanwhile, using a natural population with 505 B. napus accessions, we explored the genetic effects of BnaTCS genes on salt tolerance by association mapping analysis and detected some significant association SNPs/genes. The result will help to further understand the functions of TCS genes in the developmental and stress tolerance improvement in B. napus.

Targeted editing of multiple homologues of GTR1 and GTR2 genes provides the ideal low-seed, high-leaf glucosinolate oilseed mustard with uncompromised defence and yield

Glucosinolate content in the two major oilseed Brassica crops-rapeseed and mustard has been reduced to the globally accepted Canola quality level (<30 μmoles/g of seed dry weight, DW), making the protein-rich seed meal useful as animal feed. However, the overall lower glucosinolate content in seeds as well as in the other parts of such plants renders them vulnerable to biotic challenges. We report CRISPR/Cas9-based editing of glucosinolate transporter (GTR) family genes in mustard (Brassica juncea) to develop ideal lines with the desired low seed glucosinolate content (SGC) while maintaining high glucosinolate levels in the other plant parts for uncompromised plant defence. Use of three gRNAs provided highly efficient and precise editing of four BjuGTR1 and six BjuGTR2 homologues leading to a reduction of SGC from 146.09 μmoles/g DW to as low as 6.21 μmoles/g DW. Detailed analysis of the GTR-edited lines showed higher accumulation and distributional changes of glucosinolates in the foliar parts. However, the changes did not affect the plant defence and yield parameters. When tested against the pathogen Sclerotinia sclerotiorum and generalist pest Spodoptera litura, the GTR-edited lines displayed a defence response at par or better than that of the wild-type line. The GTR-edited lines were equivalent to the wild-type line for various seed yield and seed quality traits. Our results demonstrate that simultaneous editing of multiple GTR1 and GTR2 homologues in mustard can provide the desired low-seed, high-leaf glucosinolate lines with an uncompromised defence and yield.

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.

The current scenario and future perspectives of transgenic oilseed mustard by CRISPR-Cas9

PURPOSE

Production of a designer crop having added attributes is the primary goal of all plant biotechnologists. Specifically, development of a crop with a simple biotechnological approach and at a rapid pace is most desirable. Genetic engineering enables us to displace genes among species. The newly incorporated foreign gene(s) in the host genome can create a new trait(s) by regulating the genotypes and/or phenotypes. The advent of the CRISPR-Cas9 tools has enabled the modification of a plant genome easily by introducing mutation or replacing genomic fragment. Oilseed mustard varieties (e.g., Brassica juncea, Brassica nigra, Brassica napus, and Brassica carinata) are one such plants, which have been transformed with different genes isolated from the wide range of species. Current reports proved that the yield and value of oilseed mustard has been tremendously improved by the introduction of stably inherited new traits such as insect and herbicide resistance. However, the genetic transformation of oilseed mustard remains incompetent due to lack of potential plant transformation systems. To solve numerous complications involved in genetically modified oilseed mustard crop varieties regeneration procedures, scientific research is being conducted to rectify the unwanted complications. Thus, this study provides a broader overview of the present status of new traits introduced in each mentioned varieties of oilseed mustard plant by different genetical engineering tools, especially CRISPR-Cas9, which will be useful to improve the transformation system of oilseed mustard crop plants.

METHODS

This review presents recent improvements made in oilseed mustard genetic engineering methodologies by using CRISPR-Cas9 tools, present status of new traits introduced in oilseed mustard plant varieties.

RESULTS

The review highlighted that the transgenic oilseed mustard production is a challenging process and the transgenic varieties of oilseed mustard provide a powerful tool for enhanced mustard yield. Over expression studies and silencing of desired genes provide functional importance of genes involved in mustard growth and development under different biotic and abiotic stress conditions. Thus, it can be expected that in near future CRISPR can contribute enormously in improving the mustard plant's architecture and develop stress resilient oilseed mustard plant species.

Genome-wide association study reveals candidate genes controlling root system architecture under low phosphorus supply at seedling stage in Brassica napus

Optimal root system architecture (RSA) is essential for vigorous growth and yield in crops. Plants have evolved adaptive mechanisms in response to low phosphorus (LP) stress, and one of those is changes in RSA. Here, more than five million single-nucleotide polymorphisms (SNPs) obtained from whole-genome re-sequencing data (WGR) of an association panel of 370 oilseed rape (Brassica napus L.) were used to conduct a genome-wide association study (GWAS) of RSA traits of the panel at LP in "pouch and wick" system. Fifty-two SNPs were forcefully associated with lateral root length (LRL), total root length (TRL), lateral root density (LRD), lateral root number (LRN), mean lateral root length (MLRL), and root dry weight (RDW) at LP. There were significant correlations between phenotypic variation and the number of favorable alleles of the associated loci on chromosomes A06 (chrA06_20030601), C03 (chrC03_3535483), and C07 (chrC07_42348561), respectively. Three candidate genes (BnaA06g29270D, BnaC03g07130D, and BnaC07g43230D) were detected by combining transcriptome, candidate gene association analysis, and haplotype analysis. Cultivar carrying "CCGC" at BnaA06g29270DHap1, "CAAT" at BnaC03g07130DHap1, and "ATC" at BnaC07g43230DHap1 had greater LRL, LRN, and RDW than lines carrying other haplotypes at LP supply. The RSA of a cultivar harboring the three favorable haplotypes was further confirmed by solution culture experiments. These findings define exquisite insights into genetic architectures underlying B. napus RSA at LP and provide valuable gene resources for root breeding.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01411-2.

New vegetable varieties of Brassica rapa and Brassica napus with modified glucosinolate content obtained by mass selection approach

Background

Glucosinolates (GSLs) constitute a characteristic group of secondary metabolites present in the Brassica genus. These compounds confer resistance to pests and diseases. Moreover, they show allelopathic and anticarcinogenic effects. All those effects are dependent on the chemical structure of the GSL. The modification of the content of specific GSLs would allow obtaining varieties with enhanced resistance and/or improved health benefits. Moreover, the attainment of varieties with the same genetic background but with divergent GSLs concentration will prompt the undertaking of studies on their biological effects.

Objective and Methods

The objective of this study was to evaluate the efficacy of two divergent mass selection programs to modify GSL content in the leaves of two Brassica species: nabicol (Brassica napus L.), selected by glucobrassicanapin (GBN), and nabiza (Brassica rapa L.), selected by gluconapin (GNA) through several selection cycles using cromatographic analysis.

Results

The response to selection fitted a linear regression model with no signs of variability depletion for GSL modification in either direction, but with higher efficiency in reducing the selected GSL than in the increasing. The selection was also effective in other parts of the plant, suggesting that there is a GSL translocation in the plant or a modification in their synthesis pathway that is not-organ specific. There was an indirect response of selection in other GSL; thus this information should be considered when designing breeding programs. Finally, populations obtained by selection have the same agronomic performance or even better than the original population.

Conclusion

Therefore, mass selection seems to be a good method to modify the content of specific GSL in Brassica crops.

Integrating genome-wide association studies with selective sweep reveals genetic loci associated with tolerance to low phosphate availability in Brassica napus

Oilseed rape (Brassica napus L.; B. napus) is an important oil crop worldwide. However, the genetic mechanisms of B. napus adaptations to low phosphate (P) stress are largely unknown. In this study, a genome-wide association study (GWAS) identified 68 SNPs significantly associated with seed yield (SY) under low P (LP) availability, and 7 SNPs significantly associated with phosphorus efficiency coefficient (PEC) in two trials. Among these SNPs, two, chrC07__39807169 and chrC09__14194798, were co-detected in two trials, and BnaC07.ARF9 and BnaC09.PHT1;2 were identified as candidate genes of them, respectively, by combining GWAS with quantitative reverse-transcription PCR (qRT-PCR). There were significant differences in the gene expression level of BnaC07.ARF9 and BnaC09.PHT1;2 between P-efficient and -inefficiency varieties at LP. SY_LP had a significant positive correlation with the gene expression level of both BnaC07.ARF9 and BnaC09.PHT1;2. BnaC07.ARF9 and BnaA01.PHR1 could directly bind the promoters of BnaA01.PHR1 and BnaC09.PHT1;2, respectively. Selective sweep analysis was conducted between ancient and derived B. napus, and detected 1280 putative selective signals. Within the selected region, a large number of genes related to P uptake, transport, and utilization were detected, such as purple acid phosphatase (PAP) family genes and phosphate transporter (PHT) family genes. These findings provide novel insights into the molecular targets for breeding P efficiency varieties in B. napus.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-023-01399-9.

Genetic variation in BnGRP1 contributes to low phosphorus tolerance in Brassica napus

Lack of phosphorus (P) is a major environmental factor affecting rapeseed (Brassica napus. L) root growth and development. For breeding purposes, it is crucial to identify the molecular mechanisms of root system architecture (RSA) traits underlying low P tolerance in rapeseed. The natural variations in the glycine-rich protein gene, BnGRP1, were analyzed in the natural population of 400 rapeseed cultivars under low P stress through genome-wide association study (GWAS) and transcriptome analyses. Based on 11 SNP mutations in BnGRP1 sequence, ten types of haplotypes (Hap) were formed. Compared with the other types, the cultivar of BnGRP1Hap1 type in the panel demonstrated the longest root length and heaviest root weight. BnGRP1Hap1 overexpression in rapeseed depicted the ability to enhance its resistance in response to low P tolerance. CRISPR/Cas9-derived BnGRP1Hap4 knockout mutations in rapeseed can lead to sensitivity to low P stress. Furthermore, BnGRP1Hap1 influences the expression of phosphate transporter 1 (PHT1) genes associated with P absorption. Overall, the findings of this study highlight new mechanisms of GRP1 genes in enhancing low P tolerance in rapeseed.

CRISPR-mediated technology for seed oil improvement in rapeseed: Challenges and future perspectives

Rapeseed not only provide considerable amount of edible oil with high nutritional properties but can also be used as a raw material for biofuel production in many industries. It is therefore in high demand to bring genetic changes in order to fulfill the need of human and of industries. Though traditional breeding techniques such as hybridization and mutagenesis remained the top methods for long time to create improved varieties in oilseed rape. Clustered regularly interspaced short palindromic repeats (CRISPR) is becoming one of the most valuable gene editing technologies that allow precise genome engineering, and open new ways for research in plant functional genomics. Though CRISPR has been used in many other crops for genetic improvement it is expected to be an effective tool for genome editing and molecular design in oilseed rape for seed oil improvement. This mini review will discuss and summarize the past and ongoing research and development in rapeseed in terms of seed oil improvement and fatty acid composition using CRISPR technology. In addition, the factors that hinder the efficiency of this tool and how to eliminate those factors will be briefly summarized. The improvement of CRISPR technology for getting better results in oilseed rape will also be considered here. This minireview will open new windows for researchers in Brassica napus oil improvement research and genetic improvement using CRISPR technology.

Developing multifunctional crops by engineering Brassicaceae glucosinolate pathways

Glucosinolates (GSLs), found mainly in species of the Brassicaceae family, are one of the most well-studied classes of secondary metabolites. Produced by the action of myrosinase on GSLs, GSL-derived hydrolysis products (GHPs) primarily defend against biotic stress in planta. They also significantly affect the quality of crop products, with a subset of GHPs contributing unique food flavors and multiple therapeutic benefits or causing disagreeable food odors and health risks. Here, we explore the potential of these bioactive functions, which could be exploited for future sustainable agriculture. We first summarize our accumulated understanding of GSL diversity and distribution across representative Brassicaceae species. We then systematically discuss and evaluate the potential of exploited and unutilized genes involved in GSL biosynthesis, transport, and hydrolysis as candidate GSL engineering targets. Benefiting from available information on GSL and GHP functions, we explore options for multifunctional Brassicaceae crop ideotypes to meet future demand for food diversification and sustainable crop production. An integrated roadmap is subsequently proposed to guide ideotype development, in which maximization of beneficial effects and minimization of detrimental effects of GHPs could be combined and associated with various end uses. Based on several use-case examples, we discuss advantages and limitations of available biotechnological approaches that may contribute to effective deployment and could provide novel insights for optimization of future GSL engineering.

Genome-Wide Association Study of Glucosinolate Metabolites (mGWAS) in Brassica napus L

Glucosinolates (GSLs) are secondary plant metabolites that are enriched in rapeseed and related Brassica species, and they play important roles in defense due to their anti-nutritive and toxic properties. Here, we conducted a genome-wide association study of six glucosinolate metabolites (mGWAS) in rapeseed, including three aliphatic glucosinolates (m145 gluconapin, m150 glucobrassicanapin and m151 progoitrin), one aromatic glucosinolate (m157 gluconasturtiin) and two indole glucosinolates (m165 indolylmethyl glucosinolate and m172 4-hydroxyglucobrassicin), respectively. We identified 113 candidate intervals significantly associated with these six glucosinolate metabolites. In the genomic regions linked to the mGWAS peaks, 187 candidate genes involved in glucosinolate biosynthesis (e.g., BnaMAM1, BnaGGP1, BnaSUR1 and BnaMYB51) and novel genes (e.g., BnaMYB44, BnaERF025, BnaE2FC, BnaNAC102 and BnaDREB1D) were predicted based on the mGWAS, combined with analysis of differentially expressed genes. Our results provide insight into the genetic basis of glucosinolate biosynthesis in rapeseed and should facilitate marker-based breeding for improved seed quality in Brassica species.

Going broad and deep: sequencing-driven insights into plant physiology, evolution, and crop domestication

Deep sequencing is a term that has become embedded in the plant genomic literature in recent years and with good reason. A torrent of (largely) high-quality genomic and transcriptomic data has been collected and most of this has been publicly released. Indeed, almost 1000 plant genomes have been reported (www.plabipd.de) and the 2000 Plant Transcriptomes Project has long been completed. The EarthBioGenome project will dwarf even these milestones. That said, massive progress in understanding plant physiology, evolution, and crop domestication has been made by sequencing broadly (across a species) as well as deeply (within a single individual). We will outline the current state of the art in genome and transcriptome sequencing before we briefly review the most visible of these broad approaches, namely genome-wide association and transcriptome-wide association studies, as well as the compilation of pangenomes. This will include both (i) the most commonly used methods reliant on single nucleotide polymorphisms and short InDels and (ii) more recent examples which consider structural variants. We will subsequently present case studies exemplifying how their application has brought insight into either plant physiology or evolution and crop domestication. Finally, we will provide conclusions and an outlook as to the perspective for the extension of such approaches to different species, tissues, and biological processes.

Transcriptomic comparison of seeds and silique walls from two rapeseed genotypes with contrasting seed oil content

Silique walls play pivotal roles in contributing photoassimilates and nutrients to fuel seed growth. However, the interaction between seeds and silique walls impacting oil biosynthesis is not clear during silique development. Changes in sugar, fatty acid and gene expression during Brassica napus silique development of L192 with high oil content and A260 with low oil content were investigated to identify key factors affecting difference of their seed oil content. During the silique development, silique walls contained more hexose and less sucrose than seeds, and glucose and fructose contents in seeds and silique walls of L192 were higher than that of A260 at 15 DAF, and sucrose content in the silique walls of L192 were lower than that of A260 at three time points. Genes related to fatty acid biosynthesis were activated over time, and differences on fatty acid content between the two genotypes occurred after 25 DAF. Genes related to photosynthesis expressed more highly in silique walls than in contemporaneous seeds, and were inhibited over time. Gene set enrichment analysis suggested photosynthesis were activated in L192 at 25 and 35 DAF in silique walls and at both 15 and 35 DAF in the seed. Expressions of sugar transporter genes in L192 was higher than that in A260, especially at 35 DAF. Expressions of genes related to fatty acid biosynthesis, such as BCCP2s, bZIP67 and LEC1s were higher in L192 than in A260, especially at 35 DAF. Meanwhile, genes related to oil body proteins were expressed at much lower levels in L192 than in A260. According to the WGCNA results, hub modules, such as ME.turquoise relative to photosynthesis, ME.green relative to embryo development and ME.yellow relative to lipid biosynthesis, were identified and synergistically regulated seed development and oil accumulation. Our results are helpful for understanding the mechanism of oil accumulation of seeds in oilseed rape for seed oil content improvement.

Natural variation and artificial selection at the BnaC2.MYB28 locus modulate Brassica napus seed glucosinolate

The degradation products of glucosinolates (GSLs) greatly lower the nutritional value of rapeseed (Brassica napus) meal; thus, reduction of seed GSL content (SGC) has become an important objective of rapeseed breeding. In our previous study, we finely mapped a major QTL (qGSL-C2) for SGC to a 49-kb collinear region on B. rapa chromosome A2. Here, we experimentally validated that BnaC2.MYB28, encoding an R2R3-MYB transcription factor, is the causal gene of qGSL-C2. BnaC2.MYB28 is a nucleus-localized protein mainly expressed in vegetative tissues. Knockout of BnaC2.MYB28 in the high-SGC parent G120 reduced SGC to a value lower than that in the low-SGC parent ZY50, while overexpression of BnaC2.MYB28 in both parental lines (G120 and ZY50) led to extremely high SGC, indicating that BnaC2.MYB28 acts as a positive regulator of SGC in both parents. Molecular characterization revealed that BnaC2.MYB28 forms a homodimer and specifically interacts with BnaMYC3. Moreover, BnaC2.MYB28 can directly activate the expression of GSL biosynthesis genes. Differential expression abundance resulting from the polymorphic promoter sequences, in combination with the different capability in activating downstream genes involved in aliphatic GSL biosynthesis, caused the functional divergence of BnaC2.MYB28 in SGC regulation between the parents. Natural variation of BnaC2.MYB28 was highly associated with SGC in natural germplasm and has undergone artificial selection in modern low-GSL breeding. This study provides important insights into the core function of BnaC2.MYB28 in regulating SGC and a promising strategy for manipulating SGC in rapeseed.

Genome-wide detection of genotype environment interactions for flowering time in Brassica napus

Flowering time is strongly related to the environment, while the genotype-by-environment interaction study for flowering time is lacking in Brassica napus. Here, a total of 11,700,689 single nucleotide polymorphisms in 490 B. napus accessions were used to associate with the flowering time and related climatic index in eight environments using a compressed variance-component mixed model, 3VmrMLM. As a result, 19 stable main-effect quantitative trait nucleotides (QTNs) and 32 QTN-by-environment interactions (QEIs) for flowering time were detected. Four windows of daily average temperature and precipitation were found to be climatic factors highly correlated with flowering time. Ten main-effect QTNs were found to be associated with these flowering-time-related climatic indexes. Using differentially expressed gene (DEG) analysis in semi-winter and spring oilseed rapes, 5,850 and 5,511 DEGs were found to be significantly expressed before and after vernalization. Twelve and 14 DEGs, including 7 and 9 known homologs in Arabidopsis, were found to be candidate genes for stable QTNs and QEIs for flowering time, respectively. Five DEGs were found to be candidate genes for main-effect QTNs for flowering-time-related climatic index. These candidate genes, such as BnaFLCs, BnaFTs, BnaA02.VIN3, and BnaC09.PRR7, were further validated by the haplotype, selective sweep, and co-expression networks analysis. The candidate genes identified in this study will be helpful to breed B. napus varieties adapted to particular environments with optimized flowering time.

Comprehensive transcriptional variability analysis reveals gene networks regulating seed oil content of Brassica napus

BACKGROUND

Regulation of gene expression plays an essential role in controlling the phenotypes of plants. Brassica napus (B. napus) is an important source for the vegetable oil in the world, and the seed oil content is an important trait of B. napus.

RESULTS

We perform a comprehensive analysis of the transcriptional variability in the seeds of B. napus at two developmental stages, 20 and 40 days after flowering (DAF). We detect 53,759 and 53,550 independent expression quantitative trait loci (eQTLs) for 79,605 and 76,713 expressed genes at 20 and 40 DAF, respectively. Among them, the local eQTLs are mapped to the adjacent genes more frequently. The adjacent gene pairs are regulated by local eQTLs with the same open chromatin state and show a stronger mode of expression piggybacking. Inter-subgenomic analysis indicates that there is a feedback regulation for the homoeologous gene pairs to maintain partial expression dosage. We also identify 141 eQTL hotspots and find that hotspot87-88 co-localizes with a QTL for the seed oil content. To further resolve the regulatory network of this eQTL hotspot, we construct the XGBoost model using 856 RNA-seq datasets and the Basenji model using 59 ATAC-seq datasets. Using these two models, we predict the mechanisms affecting the seed oil content regulated by hotspot87-88 and experimentally validate that the transcription factors, NAC13 and SCL31, positively regulate the seed oil content.

CONCLUSIONS

We comprehensively characterize the gene regulatory features in the seeds of B. napus and reveal the gene networks regulating the seed oil content of B. napus.

Genome-wide association study reveals a GLYCOGEN SYNTHASE KINASE 3 gene regulating plant height in Brassica napus

Rapeseed (Brassica napus) is an allotetraploid crop that is the main source of edible oils and feed proteins in the world. The ideal plant architecture breeding is a major objective of rapeseed breeding and determining the appropriate plant height is a key element of the ideal plant architecture. Therefore, this study aims to improve the understanding of the genetic controls underlying plant height. The plant heights of 230 rapeseed accessions collected worldwide were investigated in field experiments over two consecutive years in Wuhan, China. Whole-genome resequencing of these accessions yielded a total of 1,707,194 informative single nucleotide polymorphisms (SNPs) that were used for genome-wide association analysis (GWAS). GWAS and haplotype analysis showed that BnaA01g09530D, which encodes BRASSINOSTEROID-INSENSITIVE 2 and belongs to the GLYCOGEN SYNTHASE KINASE 3 (GSK3) family, was significantly associated with plant height in B. napus. Moreover, a total of 31 BnGSK3s with complete domains were identified from B. napus genome and clustered into four groups according to phylogenetic analysis, gene structure, and motif distribution. The expression patterns showed that BnGSK3s exhibited significant differences in 13 developmental tissues in B. napus, suggesting that BnGSK3s may be involved in tissue-specific development. Sixteen BnGSK3 genes were highly expressed the in shoot apical meristem, which may be related to plant height or architecture development. These results are important for providing new haplotypes of plant height in B. napus and for extending valuable genetic information for rapeseed genetic improvement of plant architecture.

Perspectives for integrated insect pest protection in oilseed rape breeding

In the past, breeding for incorporation of insect pest resistance or tolerance into cultivars for use in integrated pest management schemes in oilseed rape/canola (Brassica napus) production has hardly ever been approached. This has been largely due to the broad availability of insecticides and the complexity of dealing with high-throughput phenotyping of insect performance and plant damage parameters. However, recent changes in the political framework in many countries demand future sustainable crop protection which makes breeding approaches for crop protection as a measure for pest insect control attractive again. At the same time, new camera-based tracking technologies, new knowledge-based genomic technologies and new scientific insights into the ecology of insect-Brassica interactions are becoming available. Here we discuss and prioritise promising breeding strategies and direct and indirect breeding targets, and their time-perspective for future realisation in integrated insect pest protection of oilseed rape. In conclusion, researchers and oilseed rape breeders can nowadays benefit from an array of new technologies which in combination will accelerate the development of improved oilseed rape cultivars with multiple insect pest resistances/tolerances in the near future.

Genetic dissection of branch architecture in oilseed rape (Brassica napus L.) germplasm

Branch architecture is an important factor influencing rapeseed planting density, mechanized harvest, and yield. However, its related genes and regulatory mechanisms remain largely unknown. In this study, branch angle (BA) and branch dispersion degree (BD) were used to evaluate branch architecture. Branch angle exhibited a dynamic change from an increase in the early stage to a gradual decrease until reaching a stable state. Cytological analysis showed that BA variation was mainly due to xylem size differences in the vascular bundle of the branch junction. The phenotypic analysis of 327 natural accessions revealed that BA in six environments ranged from 24.3° to 67.9°, and that BD in three environments varied from 4.20 cm to 21.4 cm, respectively. A total of 115 significant loci were detected through association mapping in three models (MLM, mrMLM, and FarmCPU), which explained 0.53%-19.4% of the phenotypic variations. Of them, 10 loci were repeatedly detected in different environments and models, one of which qBAD.A03-2 was verified as a stable QTL using a secondary segregation population. Totally, 1066 differentially expressed genes (DEGs) were identified between branch adaxial- and abaxial- sides from four extremely large or small BA/BD accessions through RNA sequencing. These DEGs were significantly enriched in the pathways related to auxin biosynthesis and transport as well as cell extension such as indole alkaloid biosynthesis, other glycan degradation, and fatty acid elongation. Four known candidate genes BnaA02g16500D (PIN1), BnaA03g10430D (PIN2), BnaC03g06250D (LAZY1), and BnaC06g20640D (ARF17) were identified by both GWAS and RNA-seq, all of which were involved in regulating the asymmetric distribution of auxins. Our identified association loci and candidate genes provide a theoretical basis for further study of gene cloning and genetic improvement of branch architecture.

Integration of multi-omics technologies for crop improvement: Status and prospects

With the rapid development of next-generation sequencing (NGS), multi-omics techniques have been emerging as effective approaches for crop improvement. Here, we focus mainly on addressing the current status and future perspectives toward omics-related technologies and bioinformatic resources with potential applications in crop breeding. Using a large amount of omics-level data from the functional genome, transcriptome, proteome, epigenome, metabolome, and microbiome, clarifying the interaction between gene and phenotype formation will become possible. The integration of multi-omics datasets with pan-omics platforms and systems biology could predict the complex traits of crops and elucidate the regulatory networks for genetic improvement. Different scales of trait predictions and decision-making models will facilitate crop breeding more intelligent. Potential challenges that integrate the multi-omics data with studies of gene function and their network to efficiently select desirable agronomic traits are discussed by proposing some cutting-edge breeding strategies for crop improvement. Multi-omics-integrated approaches together with other artificial intelligence techniques will contribute to broadening and deepening our knowledge of crop precision breeding, resulting in speeding up the breeding process.

Brassica Genus Seeds: A Review on Phytochemical Screening and Pharmacological Properties

Traditionally, Brassica species are widely used in traditional medicine, human food, and animal feed. Recently, special attention has been dedicated to Brassica seeds as source of health-promoting phytochemicals. This review provides a summary of recent research on the Brassica seed phytochemistry, bioactivity, dietary importance, and toxicity by screening the major online scientific database sources and papers published in recent decades by Elsevier, Springer, and John Wiley. The search was conducted covering the period from January 1964 to July 2022. Phytochemically, polyphenols, glucosinolates, and their degradation products were the predominant secondary metabolites in seeds. Different extracts and their purified constituents from seeds of Brassica species have been found to possess a wide range of biological properties including antioxidant, anticancer, antimicrobial, anti-inflammatory, antidiabetic, and neuroprotective activities. These valuable functional properties of Brassica seeds are related to their richness in active compounds responsible for the prevention and treatment of various chronic diseases such as obesity, diabetes, cancer, and COVID-19. Currently, the potential properties of Brassica seeds and their components are the main focus of research, but their toxicity and health risks must also be accounted for.

Genomic Regions Associated With Seed Meal Quality Traits in Brassica napus Germplasm

The defatted Brassica napus (rapeseed) meal can be high-protein feed for livestock as the protein value of rapeseed meal is higher than that of the majority of other vegetable proteins. Extensive work has already been carried out on developing canola rapeseed where the focus was on reducing erucic acid and glucosinolate content, with less consideration to other antinutritional factors such as tannin, phytate, sinapine, crude fiber, etc. The presence of these antinutrients limits the use and marketing of rapeseed meals and a significant amount of it goes unused and ends up as waste. We investigated the genetic architecture of crude protein, methionine, tryptophan, total phenols, β-carotene, glucosinolates (GLSs), phytate, tannins, sinapine, and crude fiber content of defatted seed meal samples by conducting a genome-wide association study (GWAS), using a diversity panel comprising 96 B. napus genotypes. Genotyping by sequencing was used to identify 77,889 SNPs, spread over 19 chromosomes. Genetic diversity and phenotypic variations were generally high for the studied traits. A total of eleven genotypes were identified which showed high-quality protein, high antioxidants, and lower amount of antinutrients. A significant negative correlation between protein and limiting amino acids and a significant positive correlation between GLS and phytic acid were observed. General and mixed linear models were used to estimate the association between the SNP markers and the seed quality traits and quantile-quantile (QQ) plots were generated to allow the best-fit algorithm. Annotation of genomic regions around associated SNPs helped to predict various trait-related candidates such as ASP2 and EMB1027 (amino acid biosynthesis); HEMA2, GLU1, and PGM (tryptophan biosynthesis); MS3, CYSD1, and MTO1 (methionine biosynthesis); LYC (β-carotene biosynthesis); HDR and ISPF (MEP pathway); COS1 (riboflavin synthesis); UGT (phenolics biosynthesis); NAC073 (cellulose and hemicellulose biosynthesis); CYT1 (cellulose biosynthesis); BGLU45 and BGLU46 (lignin biosynthesis); SOT12 and UGT88A1 (flavonoid pathway); and CYP79A2, DIN2, and GSTT2 (GLS metabolism), etc. The functional validation of these candidate genes could confirm key seed meal quality genes for germplasm enhancement programs directed at improving protein quality and reducing the antinutritional components in B. napus.

The application of CRISPR/Cas technologies to Brassica crops: current progress and future perspectives

Brassica species are a global source of nutrients and edible vegetable oil for humans. However, all commercially important Brassica crops underwent a whole-genome triplication event, hindering the development of functional genomics and breeding programs. Fortunately, clustered regularly interspaced short palindromic repeat (CRISPR)/CRISPR-associated (Cas) technologies, by allowing multiplex and precise genome engineering, have become valuable genome-editing tools and opened up new avenues for biotechnology. Here, we review current progress in the use of CRISPR/Cas technologies with an emphasis on the latest breakthroughs in precise genome editing. We also summarize the application of CRISPR/Cas technologies to Brassica crops for trait improvements. Finally, we discuss the challenges and future directions of these technologies for comprehensive application in Brassica crops. Ongoing advancement in CRISPR/Cas technologies, in combination with other achievements, will play a significant role in the genetic improvement and molecular breeding of Brassica 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.

Enhancing canola breeding by editing a glucosinolate transporter gene lacking natural variation

A low seed glucosinolate resource was developed in polyploid B. napus using a method that identifies the functions of genes with rare or no genetic variation.

Multi-omics analysis dissects the genetic architecture of seed coat content in Brassica napus

Background

Brassica napus is an important vegetable oil source worldwide. Seed coat content is a complex quantitative trait that negatively correlates with the seed oil content in B. napus.

Results

Here we provide insights into the genetic basis of natural variation of seed coat content by transcriptome-wide association studies (TWAS) and genome-wide association studies (GWAS) using 382 B. napus accessions. By population transcriptomic analysis, we identify more than 700 genes and four gene modules that are significantly associated with seed coat content. We also characterize three reliable quantitative trait loci (QTLs) controlling seed coat content by GWAS. Combining TWAS and correlation networks of seed coat content-related gene modules, we find that BnaC07.CCR-LIKE (CCRL) and BnaTT8s play key roles in the determination of the trait by modulating lignin biosynthesis. By expression GWAS analysis, we identify a regulatory hotspot on chromosome A09, which is involved in controlling seed coat content through BnaC07.CCRL and BnaTT8s. We then predict the downstream genes regulated by BnaTT8s using multi-omics datasets. We further experimentally validate that BnaCCRL and BnaTT8 positively regulate seed coat content and lignin content. BnaCCRL represents a novel identified gene involved in seed coat development. Furthermore, we also predict the key genes regulating carbon allocation between phenylpropane compounds and oil during seed development in B. napus.

Conclusions

This study helps us to better understand the complex machinery of seed coat development and provides a genetic resource for genetic improvement of seed coat content in B. napus breeding.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02647-5.

Automatic identification and annotation of MYB gene family members in plants

BACKGROUND

MYBs are among the largest transcription factor families in plants. Consequently, members of this family are involved in a plethora of processes including development and specialized metabolism. The MYB families of many plant species were investigated in the last two decades since the first investigation looked at Arabidopsis thaliana. This body of knowledge and characterized sequences provide the basis for the identification, classification, and functional annotation of candidate sequences in new genome and transcriptome assemblies.

RESULTS

A pipeline for the automatic identification and functional annotation of MYBs in a given sequence data set was implemented in Python. MYB candidates are identified, screened for the presence of a MYB domain and other motifs, and finally placed in a phylogenetic context with well characterized sequences. In addition to technical benchmarking based on existing annotation, the transcriptome assembly of Croton tiglium and the annotated genome sequence of Castanea crenata were screened for MYBs. Results of both analyses are presented in this study to illustrate the potential of this application. The analysis of one species takes only a few minutes depending on the number of predicted sequences and the size of the MYB gene family. This pipeline, the required bait sequences, and reference sequences for a classification are freely available on github: https://github.com/bpucker/MYB_annotator .

CONCLUSIONS

This automatic annotation of the MYB gene family in novel assemblies makes genome-wide investigations consistent and paves the way for comparative studies in the future. Candidate genes for in-depth analyses are presented based on their orthology to previously characterized sequences which allows the functional annotation of the newly identified MYBs with high confidence. The identification of orthologs can also be harnessed to detect duplication and deletion events.

Multi-Functional Development and Utilization of Rapeseed: Comprehensive Analysis of the Nutritional Value of Rapeseed Sprouts

Rapeseed is the third largest oil crop in the world and the largest oil crop in China. The multi-functional development and utilization of rapeseed is an effective measure for the high-quality development of rapeseed industry in China. In this study, several basic nutrients of eight rapeseed sprouts and five bean sprouts (3–5 varieties each) were determined, including sugar, crude protein, crude fiber, vitamin E, minerals, fatty acids, amino acids, and glucosinolates. Data analysis revealed that compared with bean sprouts, rapeseed sprouts were nutritionally balanced and were richer in active nutrients such as glucose, magnesium, selenium, vitamin E, and glucosinolate. Moreover, rapeseed sprouts exhibited reasonable amino acid composition and abundant unsaturated fatty acids (accounting for 90.32% of the total fatty acids). All these results indicated the potential of rapeseed sprout as a functional vegetable. Subsequently, three dominant nutrients including vitamin E, glucosinolate, and selenium were investigated in seeds and sprouts of 44 B. napus L. varieties. The results showed that germination raised the ratio of α-tocopherol/γ-tocopherol from 0.53 in seeds to 9.65 in sprouts, greatly increasing the content of α-tocopherol with the strongest antioxidant activity among the eight isomers of vitamin E. Furthermore, germination promoted the conversion and accumulation of glucosinolate components, especially, glucoraphanin with strong anti-cancer activity with its proportion increased from 1.06% in seeds to 1.62% in sprouts. In addition, the contents of selenium, vitamin E, and glucosinolate in rapeseed sprouts were highly correlated with those in seeds. Furthermore, these three dominant nutrients varied greatly within B. napus varieties, indicating the great potential of rapeseed sprouts to be further bio-enhanced. Our findings provide reference for the multi-purpose development and utilization of rapeseed, lay a theoretical foundation for the development of rapeseed sprout into a functional vegetable, and provide a novel breeding direction.

Ancient Biosyntheses in an Oil Crop: Glucosinolate Profiles in Limnanthes alba and Its Relatives (Limnanthaceae, Brassicales)

The glucosinolate (GSL) profiles of four Limnanthaceae species, including the oil crop Limnanthes alba (meadowfoam), were investigated by an ultrahigh-performance liquid chromatography-quadrupole time-of-flight tandem mass spectrometry (UHPLC-QToF-MS/MS) analysis of desulfoGSLs after desulfation of native GSLs, supplemented by NMR of desulfated 2-hydroxy-2-methylpropylGSL and 3-methoxybenzylGSL. Leaves, roots, and seeds were investigated, providing an overview of biosynthetic capabilities in the genera Floerkea and Limnanthes. Methoxyl groups on benzylGSLs were in meta but not para positions; two 3,5-disubstituted benzylGSLs are tentatively proposed. 2-Hydroxy-2-methylpropylGSL was accompanied by an isomer that was not a previously reported GSL. The combined GSL profile of the family included GSLs derived from valine, leucine, isoleucine, phenylalanine, and tyrosine, and possibly methionine and tryptophan. Substituted indole GSLs and GSLs derived from chain-elongated amino acids or alanine were searched for but not detected. Hypothetic glycosides of GSLs were detected at low levels. Based on biochemical interpretation, we suggest biosynthetic schemes and gene families (CYP79C, GSOH) relevant for tailoring GSL profiles in Limnanthes crops.