Fine mapping and candidate gene analysis of a seed glucosinolate content QTL, qGSL-C2, in rapeseed (Brassica napus L.)

Exploration of the molecular mechanism behind a novel natural genic male-sterile mutation of 1205A in Brassica napus

The use of a male sterility hybrid seed production system has resulted in a significant increase in rapeseed yields by over 20%. Nevertheless, the mechanisms underlying male sterility remain largely unexamined. This study presents a spontaneous recessive genic male-sterile (RGMS) mutant of 1205A, which was employed to establish two two-line hybrid production systems: 1205AB and NT7G132AB. Cytological investigations reveal that the mutation occurs at the early microspore stage, resulting in premature degradation of pollen. Through inheritance analysis, linkage mapping, and bulked-segregant analysis sequencing (BSA-Seq), a single gene locus, designated Bna1205ams1, was identified within the QTL region on chrC03 (15.36-18.90 Mb). The development of three newly co-segregated kompetitive allele-specific PCR (KASP) markers, in conjunction with two traditional co-segregated markers, allowed for the refinement of the QTL of Bna1205ams1 to a segment of 181.47 kb. This refinement facilitated the identification of a candidate gene, BnaC03g27700D, through functional and expression analyses. Furthermore, the subcellular localization of BnaC03g27700D was examined. Metabolic fluctuations associated with the fertility gene were observed, particularly in processes related to aborted tapetal programmed cell death (PCD), which may contribute to reduced pollen fertility with abnormal pollen exine. A strong correlation was also established between BnaC03g27700D and thirteen metabolites. This study not only offers valuable insights into the research and practical application of plant male sterility but also serves as a case study on the genetic regulatory mechanisms governing male sterility.

Quantitative Trait Loci Identification and Candidate Genes Characterization for Indole-3-Carbinol Content in Seedlings of Brassica napus

Brassica napus is a member of the cruciferous family with rich glucosinolate (GSL) content, particularly glucobrassicin (3-indolylmethyl glucosinolate, I3M), that can be metabolized into indole-3-carbinol (I3C), a compound with promising anticancer properties. To unravel the genetic mechanism influencing I3C content in rapeseed seedlings, a comprehensive study was undertaken with a doubled haploid (DH) population. By quantitative trait loci (QTL) mapping, seven QTL that were located on A01, A07, and C04 were identified, with the most significant contribution to phenotypic variation observed on chromosome A07 (11.78%). The genes within the QTL confidence intervals (CIs) include transcription factors (TFs) and glycosyltransferases. After co-expression analysis, GSL-related regulatory network of TFs-targets was constructed and two TFs, BnaA07.ERF019 and BnaA07.NAC92, were identified as possible regulators in GSL biosynthesis. Three IGMT (glucosinolate methyltransferases) genes were found within the CIs that expressed higher in seedlings with more I3C, indicating their roles in I3C synthesis regulation. Molecular docking studies validated the binding capability of I3M to IGMTs, and those within the I3C QTL CIs have the strongest binding energy. These new discoveries offer critical insights into the genetic regulation of I3C content in rapeseed seedlings and establish a foundation for breeding high-I3C rapeseed varieties with potential health-promoting properties.

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.

Quantitative trait locus mapping and improved resistance to sclerotinia stem rot in a backbone parent of rapeseed (Brassica napus L.)

There are three main challenges to improving sclerotinia stem rot (SSR) resistance in rapeseed (Brassica napus L.). First, breeding materials such as the backbone parents have not been extensively investigated, making the findings of previous studies difficult to directly implement. Second, SSR resistance and flowering time (FT) loci are typically linked; thus, use of these loci requires sacrifice of the rapeseed growth period. Third, the SSR resistance loci in susceptible materials are often neglected, thereby reducing the richness of resistant resources. This study was conducted to investigate the stem resistance, disease index, and FT of a doubled haploid population consisting of 151 lines constructed from the backbone parent 19514A and conventional rapeseed cultivar ZY50 within multiple environments. Quantitative trait locus (QTL) mapping revealed 13 stem resistance QTLs, 9 disease index QTLs, and 20 FT QTLs. QTL meta-analysis showed that uqA04, uqC03.1, and uqC03.2 were repeatable SSR resistance QTLs derived from different parents but not affected by the FT. Based on these three QTLs, we proposed a strategy for improving the SSR resistance of 19514A and ZY50. This study improves the understanding of the resistance to rapeseed SSR and genetic basis of FT and demonstrates that SSR resistance QTLs can be mined from parents with a minimal resistance level difference, thereby supporting the application of backbone parents in related research and resistance improvement.

Construction of a Quantitative Genomic Map, Identification and Expression Analysis of Candidate Genes for Agronomic and Disease-Related Traits in Brassica napus

Rapeseed is the second most important oil crop in the world. Improving seed yield and seed oil content are the two main highlights of the research. Unfortunately, rapeseed development is frequently affected by different diseases. Extensive research has been made through many years to develop elite cultivars with high oil, high yield, and/or disease resistance. Quantitative trait locus (QTL) analysis has been one of the most important strategies in the genetic deciphering of agronomic characteristics. To comprehend the distribution of these QTLs and to uncover the key regions that could simultaneously control multiple traits, 4,555 QTLs that have been identified during the last 25 years were aligned in one unique map, and a quantitative genomic map which involved 128 traits from 79 populations developed in 12 countries was constructed. The present study revealed 517 regions of overlapping QTLs which harbored 2,744 candidate genes and might affect multiple traits, simultaneously. They could be selected to customize super-rapeseed cultivars. The gene ontology and the interaction network of those candidates revealed genes that highly interacted with the other genes and might have a strong influence on them. The expression and structure of these candidate genes were compared in eight rapeseed accessions and revealed genes of similar structures which were expressed differently. The present study enriches our knowledge of rapeseed genome characteristics and diversity, and it also provided indications for rapeseed molecular breeding improvement in the future.

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.

Evolution and functional diversification of R2R3-MYB transcription factors in plants

R2R3-MYB genes (R2R3-MYBs) form one of the largest transcription factor gene families in the plant kingdom, with substantial structural and functional diversity. However, the evolutionary processes leading to this amazing functional diversity have not yet been clearly established. Recently developed genomic and classical molecular technologies have provided detailed insights into the evolutionary relationships and functions of plant R2R3-MYBs. Here, we review recent genome-level and functional analyses of plant R2R3-MYBs, with an emphasis on their evolution and functional diversification. In land plants, this gene family underwent a large expansion by whole genome duplications and small-scale duplications. Along with this population explosion, a series of functionally conserved or lineage-specific subfamilies/groups arose with roles in three major plant-specific biological processes: development and cell differentiation, specialized metabolism, and biotic and abiotic stresses. The rapid expansion and functional diversification of plant R2R3-MYBs are highly consistent with the increasing complexity of angiosperms. In particular, recently derived R2R3-MYBs with three highly homologous intron patterns (a, b, and c) are disproportionately related to specialized metabolism and have become the predominant subfamilies in land plant genomes. The evolution of plant R2R3-MYBs is an active area of research, and further studies are expected to improve our understanding of the evolution and functional diversification of this gene family.

Identifying Key Metabolites Associated with Glucosinolate Biosynthesis in Response to Nitrogen Management Strategies in Two Rapeseed (Brassica napus) Varieties

A high glucosinolate (GSL) concentration, an undesirable substance, has severely restricted rapeseed (Brassica species) development. We performed widely targeted metabolomics analysis based on the ultra-high-performance liquid chromatography-electrospray ionization-tandem mass spectrometry (UPLC-ESI-MS/MS) technology to analyze the metabolic profiles and identify the differential metabolites and GSL components in response to different nitrogen (N) levels in two rapeseed varieties. A total of 341 metabolites and 38 GSL components were detected in the seeds. A total of 188 differential metabolites, including 34 GSL components, were identified in response to different treatments, which were mapped into 2-oxocarboxylic acid metabolism, tryptophan metabolism, and GSL biosynthesis. Key indicators of GSL components highly responsible for different N levels under two contrasting varieties were recognized, i.e., 1-methylpropyl GSL and 4-methylthiobutyl GSL. This study suggests that the efficient N management and variety selection are important strategies for developing rapeseed with low GSLs.

BnaC7.ROT3, the causal gene of cqSL-C7, mediates silique length by affecting cell elongation in Brassica napus

Siliques are a major carbohydrate source of energy for later seed development in rapeseed (Brassica napus). Thus, silique length has received great attention from breeders. We previously detected a novel quantitative trait locus cqSL-C7 that controls silique length in B. napus. Here, we further validated the cqSL-C7 locus and isolated its causal gene (BnaC7.ROT3) by map-based cloning. In 'Zhongshuang11' (parent line with long siliques), BnaC7.ROT3 encodes the potential cytochrome P450 monooxygenase CYP90C1, whereas in 'G120' (parent line with short siliques), a single nucleotide deletion in the fifth exon of BnaC7.ROT3 results in a loss-of-function truncated protein. Sub-cellular localization and expression pattern analysis revealed that BnaC7.ROT3 is a membrane-localized protein mainly expressed in leaves, flowers and siliques. Cytological observations showed that the cells in silique walls of BnaC7.ROT3-transformed positive plants were longer than those of transgene-negative plants in the background of 'G120', suggesting that BnaC7.ROT3 affects cell elongation. Haplotype analysis demonstrated that most alleles of BnaC7.ROT3 are favorable in B. napus germplasms, and its homologs may also be involved in silique length regulation. Our findings provide novel insights into the regulatory mechanisms of natural silique length variations and valuable genetic resources for the improvement of silique length in rapeseed.

Identification and Fine Mapping of a Locus Related to Leaf Up-Curling Trait (Bnuc3) in Brassica napus

Leaf trait is an important target trait in crop breeding programs. Moderate leaf curling may be a help for improving crop yield by minimizing the shadowing by leaves. Mining locus for leaf curling trait is of significance for plant genetics and breeding researches. The present study identified a novel rapeseed accession with up-curling leaf, analyzed the up-curling leaf trait inheritance, and fine mapped the locus for up-curling leaf property (Bnuc3) in Brassica napus. Genetic analysis revealed that the up-curling leaf trait is controlled by a single dominant locus, named BnUC3. We performed an association study of BnUC3 with single nucleotide polymorphism (SNP) markers using a backcross population derived from the homozygous up-curling leaf line NJAU-M1295 and the canola variety ‘zhongshuang11’ with typical flat leaves, and mapped the BnUC3 locus in a 1.92 Mb interval of chromosome A02 of B. napus. To further map BnUC3, 232 simple sequence repeat (SSR) primers and four pairs of Insertion/Deletion (InDel) primers were developed for the mapping interval. Among them, five SSR markers and two InDel markers were polymorphic. By these markers, the mapping interval was narrowed to 92.0 kb using another F₂ population. This fine mapping interval has 11 annotated genes among which BnaA02T0157000ZS were inferred to be candidate casual genes for up-curling leaf based on the cloned sequence analysis, gene functionality, and gene expression analysis. The current study laid a foundational basis for further elucidating the mechanism of BnUC3 and breeding of variety with up-curling leaf.

A 24,482-bp deletion is associated with increased seed weight in Brassica napus L

A major QTL for seed weight was fine-mapped in rapeseed, and a 24,482-bp deletion likely mediates the effect through multiple pathways. Exploration of the genes controlling seed weight is critical to the improvement of crop yield and elucidation of the mechanisms underlying seed formation in rapeseed (Brassica napus L.). We previously identified the quantitative trait locus (QTL) qSW.C9 for the thousand-seed weight (TSW) in a double haploid population constructed from F₁ hybrids between the parental accessions HZ396 and Y106. Here, we confirmed the phenotypic effects associated with qSW.C9 in BC₃F₂ populations and fine-mapped the candidate causal locus to a 266-kb interval. Sequence and expression analyses revealed that a 24,482-bp deletion in HZ396 containing six predicted genes most likely underlies qSW.C9. Differential gene expression analysis and cytological observations suggested that qSW.C9 affects both cell proliferation and cell expansion through multiple signaling pathways. After genotyping of a rapeseed diversity panel to define the haplotype structure, it could be concluded that the selection of germplasm with two specific markers may be effective in improving the seed weight of rapeseed. This study provides a solid foundation for the identification of the causal gene of qSW.C9 and offers a promising target for the breeding of higher-yielding rapeseed.

Fine mapping of the major QTLs for biochemical variation of sulforaphane in broccoli florets using a DH population

Glucoraphanin is a major secondary metabolite found in Brassicaceae vegetables, especially broccoli, and its degradation product sulforaphane plays an essential role in anticancer. The fine mapping of sulforaphane metabolism quantitative trait loci (QTLs) in broccoli florets is necessary for future marker-assisted selection strategies. In this study, we utilized a doubled haploid population consisting of 176 lines derived from two inbred lines (86,101 and 90,196) with significant differences in sulforaphane content, coupled with extensive genotypic and phenotypic data from two independent environments. A linkage map consisting of 438 simple sequence repeats markers was constructed, covering a length of 1168.26 cM. A total of 18 QTLs for sulforaphane metabolism in broccoli florets were detected, 10 were detected in 2017, and the other 8 were detected in 2018. The LOD values of all QTLs ranged from 3.06 to 14.47, explaining 1.74-7.03% of the biochemical variation between two years. Finally, 6 QTLs (qSF-C3-1, qSF-C3-2, qSF-C3-3, qSF-C3-5, qSF-C3-6 and qSF-C7) were stably detected in more than one environment, each accounting for 4.54-7.03% of the phenotypic variation explained (PVE) and a total of 30.88-34.86% of PVE. Our study provides new insights into sulforaphane metabolism in broccoli florets and marker-assisted selection breeding in Brassica oleracea crops.

Genetic Characterization of Russian Rapeseed Collection and Association Mapping of Novel Loci Affecting Glucosinolate Content

Rapeseed is the second most common oilseed crop worldwide. While the start of rapeseed breeding in Russia dates back to the middle of the 20th century, its widespread cultivation began only recently. In contrast to the world’s rapeseed genetic variation, the genetic composition of Russian rapeseed lines remained unexplored. We have addressed this question by performing genome-wide genotyping of 90 advanced rapeseed accessions provided by the All-Russian Research Institute of Oil Crops (VNIIMK). Genome-wide genetic analysis demonstrated a clear difference between Russian rapeseed varieties and the rapeseed varieties from the rest of the world, including the European ones, indicating that rapeseed breeding in Russia proceeded in its own independent direction. Hence, genetic determinants of agronomical traits might also be different in Russian rapeseed lines. To assess it, we collected the glucosinolate content data for the same 90 genotyped accessions obtained during three years and performed an association mapping of this trait. We indeed found that the loci significantly associated with glucosinolate content variation in the Russian rapeseed collection differ from those previously reported for the non-Russian rapeseed lines.