Hexaploid-bridged introgression broadens genetic diversity and enriches clubroot-resistance alleles in Brassica juncea

Introduction: The narrow genetic base of current Brassica juncea (AABB) has limited the genetic improvement of this crop.

Methods: In this study, we developed novel germplasm by crossing synthetic hexaploid bridge lines (AABBBB and AAAABB) with diploid progenitors B. rapa (AA) and B. nigra (BB), followed by successive selfing and marker-assisted selection.

Results and Discussion: Although the F1 hybrids exhibited wide variation in fertility and morphology, these traits stabilized at levels comparable to current B. juncea by the F2 and F3 generations. Whole-genome resequencing confirmed that the new-type lines not only constitute a unique genetic group but also show a genomic shift toward the current B. juncea cluster in the F2 populations. Consequently, these new-type lines represent a significant new source of genetic diversity. As a practical application, we introgressed clubroot resistance from the European fodder turnip ECD04 (B. rapa) into B. juncea using the hexaploid bridge strategy coupled with MAS. Over four selection cycles (F2 to F5), the donor allele frequency surged from 31.3% to 54.4%, showing a strong negative correlation with disease severity (r = -0.98) and yielding highly resistant F5 lines.

Conclusions: This study demonstrates that the hexaploid bridge strategy enables rapid genetic base broadening and efficient trait improvement in B. juncea.

1 Introduction

Brassica juncea (B. juncea, AABB, 2n = 36), an amphidiploid species originating from natural hybridization between the diploid progenitors Brassica nigra (B. nigra, BB, 2n = 16) and Brassica rapa (B. rapa, AA, 2n = 20) and spontaneous chromosome doubling (Kang et al., 2021). It is widely cultivated as a multipurpose crop for oilseed, vegetable, and condiment production (Tian and Deng, 2020). This species is valued for its high adaptability to abiotic stresses, including drought (Chen et al., 2013b), salinity, and pod-shattering resistance (Rahman et al., 2022). These attributes make it well-suited for cultivation in marginal environments (Limbalkar et al., 2021). However, compared to its diploid progenitors, the genetic base of cultivated B. juncea remains narrow (Hittorf et al., 2020; Zhang et al., 2022), limiting its potential for further genetic improvement through conventional breeding (Singh et al., 2022).

Traditional approaches for gene introgression from diploid progenitors into B. juncea, such as direct hybridization between B. rapa and B. nigra, are hindered by strong pre- and post-zygotic reproductive barriers, including hybrid sterility and embryo abortion, often necessitating embryo rescue with low success rates (Liu et al., 2018). Moreover, backcrossing resistance genes from B. rapa into B. juncea is labor-intensive and time-consuming, typically requiring six to eight generations to achieve sufficient recovery of the recipient genome (Zhu et al., 2022). To overcome these limitations, hexaploid bridge lines with the AABBBB or AAAABB genome, developed from crosses of B. juncea with B. rapa or B. nigra, have demonstrated enhanced crossability with the other diploid progenitor and exhibited stable fertility (> 90% pollen viability and 2–4 seeds per pod) from the S1 to S3 generation, enabling efficient introgression of alleles from diploid progenitors into B. juncea under field conditions (Mei et al., 2022).

Although the synthetic hexaploid bridge strategy has demonstrated potential for introducing exotic genomic segments into B. juncea (Mei et al., 2022), the fate of these segments under continuous selfing remains poorly understood, particularly regarding fertility restoration, phenotypic stabilization, and genomic composition across generations (Yang et al., 2022a). It is well known that early-generation interspecific hybrids often suffer from meiotic abnormalities, including multivalent formation and genome loss, leading to compromised fertility (Xiong et al., 2011; Gaeta et al., 2007). Whether fertility and phenotypic stability can be restored in subsequent generations is critical to the practical utility of this strategy.

In addition to broadening genetic variation, the hexaploid bridge also offers a foundation for developing novel germplasm with key agronomic traits such as resistance to clubroot disease (caused by Plasmodiophora brassicae), which is absent in the current B. juncea gene pool but present in B. rapa. European fodder turnip ECD04 is a well-known B. rapa line that carrying multiple CR loci and exhibits strong resistance to several Plasmodiophora brassicae races (Yang et al., 2022b). To systematically introgress these CR genes from ECD04 into the B. juncea background, marker-assisted selection (MAS) offers a viable approach for tracking and enriching desirable resistance alleles in early generations (Tomita et al., 2013; Dakouri et al., 2013). However, the efficiency of this combined strategy by using MAS within a hexaploid-derived population to rapidly increase resistance allele dosage and achieve effective phenotypic resistance has not been evaluated.

To address these questions, we developed a population of new-type B. juncea by employing the hexaploid bridge strategy and advanced it through successive selfing generations. Our study aimed to (i) track fertility restoration over selfing generations, (ii) analyze phenotypic and genomic dynamics, and (iii) demonstrate practical trait improvement by introgressing CR from ECD04 into the B. juncea background via MAS. This work not only establishes the hexaploid bridge as an efficient strategy to diversify the B. juncea gene pool, but also provides a novel and diverse germplasm resource for B. juncea breeding.

2 Materials and methods

2.1 Plant materials

An effective strategy in developing B. juncea-type hybrids, i.e., the hexaploid bridging strategy was developed in our previous study (Mei et al., 2022) and employed in the present study to generate new-type B. juncea germplasm. In brief, 32 B. nigra accessions were randomly cross-pollinated with three A^jA^jA^rA^rB^jB^j hexaploid genotypes (designated MS2, MS3 and MS4) to produce A^jA^rBⁿB^j-type F₁ hybrids, and 20 B. rapa accessions were crossed with one A^jA^jB^jB^jBⁿBⁿ (designated MS6) hexaploid to produce A^jA^rBⁿB^j-type F₁ hybrids. Resultant F₁ hybrids were self-pollinated to produce F₂, F₃ and other generations. From the above B. rapa panel, the European fodder turnip ECD04 (AA) — a known donor harboring multiple resistance genes against Plasmodiophora brassicae — was specifically crossed with the hexaploid bridge line MS6 to initiate CR gene transfer. Phenotypic observations (e.g., leaf morphology) and imaging were conducted concurrently at designated growth stages (e.g., seedling stage) to ensure direct comparability across all genotypes. Fertility was assessed by counting the number of seeds per pod from 30 consecutive pods collected from the primary inflorescence of each plant. All materials, including the parental lines (B. nigra, B. rapa, hexaploids MS2-MS6 and ECD04), the derived new-type lines (F₁–F₅), and ten control accessions of current B. juncea (vegetable and oil types) from Southwest University (SWU), were sown and transplanted on the same dates in the experimental field at SWU. A uniform planting design was maintained with 30 cm between rows and 25 cm within rows. Detailed information for all accessions and crosses is provided in Supplementary Table S1.

2.2 Genome resequencing and data analysis

To understand the genetic characteristic of new-type B. juncea germplasm, 61 new-type B. juncea individuals were randomly selected from the F₁ and F₂ generations for whole-genome resequencing, together with 20 B. nigra parents, the four hexaploid parents, three B. rapa parents, and ten current B. juncea accessions. Young leaf samples from these materials were submitted to Wuhan Benagen Biotech Co., Ltd. for whole-genome resequencing using BGISEQ-500. Sequencing yielded an average depth ranging from 9.3× to 16.6× per sample. Following whole-genome resequencing, the genomic sequence data were subjected to quality assessment and filtering using fastp, which included: removal of adapter sequences, trimming of poly G and poly X tails from reads, discarding reads containing more than 5 N bases, removing reads in which the proportion of bases with quality scores below 15 exceeded 40%, eliminating reads with an average base quality score lower than 20, and discarding reads shorter than 50 bp after filtering.

Clean reads were aligned to the reference genome of B. juncea (BjuIR: Brassica juncea information resource) using BWA. Variants were called using Genome Analysis Toolkit (https://wiki.rc.usf.edu/index.php/Genome_Analysis_ToolKit_GATK), with SNPs and InDels filtered based on specific criteria. Genetic distances among materials were calculated based on the SNPs and InDels using Powermarker V 3.25 (https://brcwebportal.cos.ncsu.edu/powermarker/). Phylogenetic trees were constructed using MEGA software (https://www.megasoftware.net/), and principal component analysis (PCA) was performed using Rstudio.

2.3 DNA extraction

Genomic DNA was isolated from fresh young leaf tissues of all plant materials using a modified cetyltrimethylammonium bromide (CTAB) method. Approximately 100 mg of tissue was ground in liquid nitrogen and incubated in pre-warmed CTAB extraction buffer at 65 °C for 45 min. After chloroform: isoamyl alcohol (24:1) extraction, DNA was precipitated with isopropanol, washed with 70% ethanol, and dissolved in nuclease-free water. DNA concentration and quality were assessed using a NanoDrop spectrophotometer and agarose gel electrophoresis.

2.4 Marker-assisted selection

To track the introgression of CR genes, genotyping was performed using four PCR-based markers (Powell et al., 1996): CR-m090a (linked to PbBa8.1), CR-BSA3 (linked to PbBa1.1), and CR-S14R14 and CR-S17R17 (both linked to BraA.CR.b) (Chen et al., 2013a; Hirani et al., 2018; primer sequences in Supplementary Table S2). All F2 individuals from the cross MS6 × ECD04 were screened using these markers, and those carrying the maximum number of ECD04-specific alleles were self-pollinated Selection was carried out in subsequent generations (F3 to F5) by genotyping ten plants per line with the same markers. PCR was carried out in 10 µL reactions containing approximately 50 ng genomic DNA, 0.5 µL each of forward and reverse primer (10 µM), 5 µL of 2× PCR Master Mix (Vazyme, China), and nuclease-free water to volume. The thermal cycling conditions were as follows: initial denaturation at 94 °C for 5 min; 10 cycles of touchdown PCR (94 °C for 30 s, annealing from 62 °C to 57 °C decreasing by 0.5 °C per cycle for 30 s, 72 °C for 30 s); followed by 25 cycles of standard PCR (94 °C for 30 s, 55 °C for 30 s, 72 °C for 30 s); and a final extension at 72 °C for 5 min. Reactions were held at 16 °C post-amplification. The products were separated by electrophoresis on a 2% agarose gel. Clear bands exceeding the background intensity were scored as 1 (presence of the resistance allele), while faint or absent bands were scored as 0. The ECD04 allele frequency (AF, %) was calculated as.

$$AF=\frac{Np}{Nt}\times 100\%$$

where Np indicates the number of plants identified as positive for the molecular marker, and Nt represents the total number of plants assayed. The line graph depicting ECD04 allele frequencies for each marker was generated using GraphPad Prism.

2.5 Clubroot resistance evaluation

To evaluate clubroot resistance, the plants were inoculated with Plasmodiophora brassicae race 4 using the injection method in lab (Jubault et al., 2013) under controlled conditions (22–25 °C, 60% relative humidity). Three replicates were set up for each tested line, with 30 seeds sown per replicate. Disease incidence (DI) and severity were assessed 42 days post-inoculation. Based on root symptoms and gall size, disease severity was classified on a 0–4 scale (Chen et al., 2021), with 0 indicating no symptoms and 4 representing the most severe root deformation. The DI was calculated as follows (Ludwig-Müller et al., 2017):

$$DI=\frac{1\times N1+2\times N2+3\times N3+4\times N4}{4\times Nt}\times 100$$

where N1-N4 indicate the number of plants in each disease severity category, and Nt represents the total number of inoculated plants. The correlation between the DI and the ECD04 AF was calculated and graphically presented using GraphPad Prism.

3 Result

3.1 Fertility characteristics of new-type B. juncea

Through hybridization between hexaploids and diploids, followed by successive selfing, we developed 180 new-type B. juncea lines, comprising 76 F₁ hybrids (from direct crosses), 104 advanced selfing lines (46 F₂, 26 F₃, 23 F₄ and 9 F₅ generation). To evaluate the fertility of the new-type B. juncea lines, we quantified the seed set rate, measured as seeds per pod, for the F₁, F₂, and F₃ generations. Statistical analysis across various generations (Figure 1A) revealed that in the F₁ and F₂ generations, hybrids derived from MS2, MS3, and MS4 (A^jA^jA^rA^rB^jB^j) exhibited higher seed set rate compared to those derived from MS6 (A^jA^jB^jB^jBⁿBⁿ), suggesting that the different hybridization combinations significantly affect the fertility of hybrids, i.e., the cross of hexaploid A^jA^jA^rA^rB^jB^j with B. nigra is easier to produce hybrids with better fertility in early generations than the cross between hexaploid A^jA^jB^jB^jBⁿBⁿ and B. rapa. More importantly, the data demonstrated that the seed set rate of the hybrids significantly improved with successive generations, particularly noticeable in the hybrids derived from MS6. This indicates that successive selfing can effectively enhance the fertility of all hybrids, regardless of whether they are hybrids of A^jA^jA^rA^rB^jB^j × A^rA^r or A^jA^jB^jB^jBⁿBⁿ × BⁿBⁿ.

Figure 1

Figure 1. Variation in seeds per pod and phenotypic traits in F₁, F₂, and F₃ new-type B . juncea. Different letters above bars indicate significant difference at p < 0.05. Images show the variation in leaf shape and color of new-type B . juncea; bars =10 cm. (A) Seeds per pod in F₁, F₂, and F₃ derived from crosses with different hexaploid parents. MS2 (A^jA^jA^rA^rB^jB^j), MS3 (A^jA^jA^rA^rB^jB^j), MS4 (A^jA^jA^rA^rB^jB^j) and MS6 (A^jA^jB^jB^jBⁿBⁿ) denote the specific hexaploid parent used in each cross combination. (B-E) Phenotypic traits in F₁ plants (images are derived from Mei et al., 2022). (F, G) Additional examples of phenotypic traits in F₁ plants (newly developed in this study). The F₁ plants of the new-type B . juncea exhibited distinct morphological characteristics, including leaf trichomes and purple venation (Figures F) as well as pinnately lobed leaves with prominent venation . (G-J) Phenotypic traits in F₂ plants. In the F₂ generation, typical F₁ traits such as trichomes and serrated leaf margins became less pronounced. (K-M) Phenotypic traits in F₃ plants. By the F₃ generation, a subset of plants converged morphologically toward existing cultivated B . juncea types (K , L), while some materials retained the leaf glossiness characteristic of B . rapa (M).

3.2 Phenotypic diversity of new-type B. juncea

Consistent with fertility restoration, the new-type B. juncea lines showed dynamic morphological changes from the F₁ to F₃ generations (Figures 1B-M). The F₁ plants exhibited considerable phenotypic variation, including pinnately lobed leaves, leaf trichomes, purple venation and stem, and prominent venation. In the F₂ generation, while phenotypic variation persisted, a few plants began to converge toward current B. juncea traits, with certain F₁ features such as trichomes and serrated leaf margins becoming less distinct. By the F₃ generation, some plants clearly converged toward the morphology of current B. juncea cultivars. These morphological changes demonstrate that successive selfing generations enabled the selection of plants with a high degree of phenotypic identity to current B. juncea from the new-type lines.

3.3 Genetic diversity of new-type B. juncea

The observed convergence in both fertility and phenotype suggested underlying genetic changes. To elucidate the genetic characteristics of the new-type B. juncea, whole-genome resequencing was conducted on the 61 new-type B. juncea lines (24 F₁ and 37 F₂), alongside 23 parental B. rapa and B. nigra accessions, 10 current B. juncea accessions (including vegetative and oilseed types), and four hexaploid parents. After strict filtering, 24.6 million high-confidence SNPs and small indels were retained, spanning the A and B sub-genomes, with a mean density of 1.36 million variants per chromosome.

As revealed by phylogenomic analysis (Figure 2A), both F₁ and F₂ generations exhibited substantial genetic variation as compared to current B. juncea accessions. Most F₁ individuals occupied a long-branch lineage diverging prior to the clade of current B. juncea, whereas F₂ individuals formed a paraphyletic group partially overlapping with current B. juncea. Consistent with this phylogenetic pattern, the principal component analysis separated most of the F₁ progeny from current B. juncea along Principal Component 1 (PC1) (19.2% of variance) and Principal Component 1 (PC2) (8.6%), while more F₂ individuals were notably shifted toward the current B. juncea cloud along PC1 (Figure 2B).

Figure 2

Figure 2. Phylogenomic and principal component analysis results of F₁and F₂ new-type B . juncea, current B . juncea, B . rapa B . nigra and hexaploid. The accessions of F₁, F_2, current B . juncea, B . rapa, B . nigra and hexaploids are denoted by the suffixes F₁, F_2, Bju, Bar, Bni and Hex, respectively. F₁ represented by teal, F₂ represented by purple, current B . jucea represented by red, B . nigra represented by blue, B . rapa represented by pink, hexaploid represented by yellow. The F₁ individuals are genetically distinct and scattered, while the F₂ individuals show a shift toward the cluster of current B. juncea. (A) Phylogenomic result; branches in different colors represent different materials. (B) Principal component analysis result; different colors and distinct shapes represent different materials. The principal components PC1 and PC2 accounted for 19.2% and 8.6% of the total variance, respectively. The confidence ellipses represent the 95% confidence intervals for the F₁ and F₂ populations.

Genetic distance (GD) analysis further quantified these relationships (Table 1). The mean within-group GD was lowest in current B. juncea (GD_Cj = 0.090), intermediate in the F₂ (GD_F2 = 0.152), and highest in the F₁ (GD_F1 = 0.168). The between-group GD to current B. juncea was slightly lower for the F₂ (GD_F2-Cj = 0.136) than for the F₁ (GD_F1-Cj = 0.146), suggesting that the F₂ lines are still markedly different from current B. juncea. However, it also suggests that genetic recombination had already produced some individuals with genomic compositions more closely aligned with current B. juncea. Collectively, these results confirm that the new-type lines constitute a unique genetic group, harboring a spectrum of diversity that includes both distinct lineages and F₂ individuals with genomic compositions closer to current B. juncea.

Table 1

Table 1. CD within and among F₁and F₂ new-type B. juncea, current B. juncea, B. rapa B. nigra and hexaploid.

3.4 Development of clubroot resistant B. juncea

To validate the practical application value of the hexaploid bridge hybridization strategy for trait improvement, we adopted a MAS in the hybrids between 6MS6 (A^jA^jB^jB^jBⁿBⁿ) with the European fodder turnip ECD04 which is well known B. rapa accession harboring multiple CR loci. Four PCR-based markers tightly linked to these CR loci were assayed in F₂-F₅ progeny (Supplementary Table S3), and plants carrying the highest number of ECD04-alleles were advanced to the next generation.

According to the four CR-linked markers (Figure 3A), a progressive enrichment was presented since the mean frequency of ECD04-alleles rose from 31.3% in the F₂ to 54.4% in the F₅, representing a 73.8% rise over the selection period (Figure 3B). More than half of the F₅ families (e.g., 5G06, 5G09, 5G10, 5G11 and 5G13) retained ≥ 50% donor alleles, while a subset (5G16-5G19) still carried< 50% ECD04 alleles (Table 2). Laboratory phenotyping for CR was carried out on F₅ lines with AF ≥ 50%, alongside parental controls. A strong negative linear correlation was observed between ECD04 AF and disease index (r = -0.98, p < 0.0001; Figure 3C), indicating that the allele dosage is a reliable predictor of resistance levels.

Figure 3

Figure 3. Changes in ECD04 allele frequency (AF, %) across generations, its negative correlation with disease resistance and genotyping of the four CR-linked markers. (A) Agarose gel electrophoresis images of PCR amplification using four CR-linked markers for parental and F₅ individuals. M: DNA molecular weight standard. Arrows indicate the specific bands corresponding to the resistance alleles. The resistance alleles were assigned based on previously reported linkage: markers CR−m090a and CR−BSA3 are linked to the resistance genes PbBa8.1 and PbBa1.1, respectively (Chen et al., 2013a); markers CR−S14R14 and CR−S17R17 are linked to the resistance locus BraA.CR.b (Hirani et al., 2018). The resistance−donor parent ECD04 is a homozygous resistant accession: for markers CR−m090a, CR−S14R14 and CR−S17R17 it shows only a single band (the resistance allele); for marker CR−BSA3, although multiple bands were amplified, the arrow points to the target resistance band. Lane labels 5G06, 5G11 and 5G18 denote F₅ families. Each lane contains DNA from a single plant that was randomly selected from the corresponding family. (B) Changes in ECD04 AF of different markers across generations. (C) The negative linear correlation between ECD04 AF and DI; r = -0.98, p < 0.0001.

Table 2

Table 2. Molecular markers and disease incidence statistics of clubroot in F₅ new-type B. juncea.

Concerning the detailed resistance (Figure 4, Table 2), the hexaploid parent and current B. juncea exhibited nearly 100% incidence and DI > 60, whereas the resistant parent ECD04 showed near immunity (disease incidence = 4.3%, DI = 1.7). Two F₅ families with AF = 77.5% (5G06 and 5G11) were virtually homozygous for the target haplotype (especially markers S14R14 and m090a) and recorded the lowest DI values (11.5 and 12.5, respectively), indicating the successful introgression and concentration of ECD04-derived CR loci into the B. juncea background. Three other F₅ families with AF between 50% and 75% displayed moderate resistance, with DI of 23.1-31.2 and disease incidence of 63%-75%, indicating the need for further selfing to achieve full homozygosity.

Figure 4

Figure 4. Lab-inoculation results and field phenotype of F₅ new-type B . juncea. (A) Phenotyping of clubroot inoculation. Current B . juncea and hexaploid were susceptible controls and ECD04 served as the resistant control; inoculated with Plasmodiophora brassicae race 4 for 42 d, bar = 10 cm. (B) Disease severity grading distribution. Based on the severity of root symptoms, the disease rating scale is divided into 5 grades: grade 0 indicates no symptoms, grade 4 represents the most severe infection and different colors represent different disease grades. (C) Field phenotype of F₅ family 5G06 at seedling stage, bar = 10 cm.

4 Discussion

4.1 New-type B. juncea, a novel genetic resource for B. juncea breeding

The strong reproductive barriers that impede traditional breeding, such as hybrid sterility and embryo abortion, have severely restricted the introgression of novel alleles from diploid progenitors into the cultivated B. juncea gene pool (Liu et al., 2018; Zhu et al., 2022). In this study, synthetic hexaploids (AABBBB/AAAABB), proposed as a bridge for gene transfer (Mei et al., 2022), were employed successfully to achieve efficient introgression of genomic components from both progenitor species into B. juncea. The newly developed materials represent a distinct and unique gene pool, not merely a random genetic mixture.

Phylogenetic, principal component, and genetic distance analyses consistently demonstrated that the new-type lines, particularly the F₁ generation, formed a cluster separate from natural B. juncea. Although a genetic shift toward natural B. juncea was observed in the F₂ generation, the genetic distance between the new-type lines and natural B. juncea remained substantially larger than that within natural B. juncea itself. These lines therefore retain a unique set of alleles from the progenitor species, collectively establishing a valuable and novel genetic reservoir. This population can serve as a diverse foundation set for future gene discovery and breeding, directly addressing the issue of genetic constriction in current cultivars.

4.2 Rapid fertility restoration and genomic stabilization through selfing

The practical utilization of this novel gene pool first requires overcoming a common challenge in newly synthesized polyploids, i.e., meiotic instability, which often leads to reduced fertility (Tian et al., 2023). The significantly lower seeds per pod in the F₁ generation in our study confirmed such instability, reflecting a form of genomic shock characterized by rapid genomic reorganization in response to aberrant homologous chromosome pairing (McClintock, 1984; Xiong et al., 2011; Martin and Jouve, 1992). Nevertheless, the fertility was rapidly restored in subsequent generations. This recovery is consistent with a similar scenario (Xiong et al., 2011), wherein meiotic stability is gradually established through restricted homologous pairing, epigenetic silencing of conflict regions, and selective elimination of cytologically unstable lineages. Concurrently, phenotypic convergence from diploid-parent-like morphological variation in the F₁ to current B. juncea type traits in the F₃ was possibly also driven by recurrent selfing. The increase in homozygosity through selfing not only facilitated stable chromosome pairing but also purged deleterious genomic interactions, thereby accelerating the emergence of a cytologically and phenotypically stable allopolyploid.

4.3 Efficient trait introgression and a clear dosage effect of CR loci

Building upon this rapid stabilization, MAS was effectively applied within the complex polyploid background. Four cycles of targeted selection substantially enriched CR alleles from ECD04, demonstrating the high efficiency and precision of the hexaploid bridge strategy for trait improvement. Notably, a strong linear relationship was observed between allelic dosage of the resistance loci and disease severity, indicating that the introgressed CR genes likely act in a dominant or additive manner (Li et al., 2016). This relationship enabled reliable phenotypic prediction in early generations based on genotypic data, significantly improving selection efficiency.

However, our findings also reveal a limitation of the dosage model. Family 5G13, which was fixed for the major locus S17 and carried 75% donor alleles across the four target loci, still exhibited 63% disease incidence. We hypothesize that the following factors may be responsible: (i) linkage drag, in which undesirable genomic segments tightly linked to the target CR loci were co-introgressed, potentially carrying alleles that compromise resistance or plant fitness (Jabran et al., 2023); (ii) residual heterozygosity at loci BSA3, S14, and m090a (87.5% donor alleles, not yet fixed), or the presence of tightly linked deleterious segments near S17, may have compromised the expected resistance; (iii) background effects from the recipient B. juncea genome, where interactions between the introgressed alleles and the genetic background may modulate the final resistance level. This case illustrates a general limitation of relying solely on foreground selection in marker-assisted breeding (Dakouri et al., 2013). Therefore, incorporating routine background selection is recommended to break such unfavorable linkages and accelerate recovery of the recurrent parent genome (Nannuru et al., 2025), ultimately facilitating the development of agronomically superior lines with stable and complete clubroot resistance.

5 Conclusion and perspectives

In summary, this study establishes an efficient framework for polyploid crop improvement via synthetic hexaploid bridges. The substantial genetic diversity within this novel germplasm pool created in this study offers a promising resource for discovering alleles underlying valuable agronomic traits from B. rapa and B. nigra. This strategy not only broadens the genetic base of B. juncea but also enables precision breeding for clubroot resistance. Future efforts will focus on multi-environment field evaluations of advanced lines to assess their yield potential and adaptability, paving the way for commercial deployment.

Data availability statement

The raw whole-genome sequencing data generated in this study have been deposited, in accordance with the journal’s requirements, in the BIG Submission database of the National Genomics Data Center (NGDC) in China. The data can be accessed via the BioProject accession number PRJCA049343 at the following link: https://ngdc.cncb.ac.cn/bioproject/browse/PRJCA049343.

Author contributions

YWu: Data curation, Formal analysis, Investigation, Software, Visualization, Writing – original draft. JL: Data curation, Formal analysis, Investigation, Visualization, Writing – review & editing. YD: Data curation, Formal analysis, Investigation, Writing – review & editing. FC: Data curation, Investigation, Resources, Validation, Writing – review & editing. PB: Resources, Writing – review & editing. LP: Resources, Writing – review & editing. XZ: Resources, Writing – review & editing. YWa: Writing – review & editing. WQ: Methodology, Project administration, Supervision, Writing – review & editing. JM: Funding acquisition, Methodology, Supervision, Validation, Writing – review & editing.

Funding

The author(s) declared that financial support was received for this work and/or its publication. This work was supported by Southwest University Fuling Academy Joint Construction Project (FLYJY202405), National Key R&D Program (2023YFD1201402) and Modern Pioneer District and County Project of Agricultural Science and Technology in Fuling of Chongqing (Modern Biotechnology Breeding of Mustard).

Acknowledgments

The authors wish to acknowledge Prof. Qinfei Li and Dr. Fang Yue in Southwest University, China, to provide the hexaploid lines.

Conflict of interest

The author(s) declared that this work was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Generative AI statement

The author(s) declared that generative AI was not used in the creation of this manuscript.

Any alternative text (alt text) provided alongside figures in this article has been generated by Frontiers with the support of artificial intelligence and reasonable efforts have been made to ensure accuracy, including review by the authors wherever possible. If you identify any issues, please contact us.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2025.1729594/full#supplementary-material

References

Chen, J., Jing, J., Zhan, Z., Zhang, T., Zhang, C., and Piao, Z. (2013a). Identification of novel QTLs for isolate-specific partial resistance to Plasmodiophora brassicae in Brassica rapa. PloS One 8, e85307. doi: 10.1371/journal.pone.0085307

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, W., Li, Y., Yan, R., Ren, L., Liu, F., Zeng, L., et al. (2021). SnRK1.1-mediated resistance of Arabidopsis thaliana to clubroot disease is inhibited by the novel Plasmodiophora brassicae effector PBZF1. Mol. Plant Pathol. 22, 1057–1069. doi: 10.1111/mpp.13095

PubMed Abstract | Crossref Full Text | Google Scholar

Chen, S., Wan, Z., Nelson, M. N., Chauhan, J. S., Redden, R., Burton, W. A., et al. (2013b). Evidence from genome-wide simple sequence repeat markers for a polyphyletic origin and secondary centers of genetic diversity of Brassica juncea in China and India. J. Heredity 104, 416–427. doi: 10.1093/jhered/est015

PubMed Abstract | Crossref Full Text | Google Scholar

Dakouri, A., McCallum, B. D., Radovanovic, N., and Cloutier, S. (2013). Molecular and phenotypic characterization of seedling and adult plant leaf rust resistance in a world wheat collection. Mol. Breed. 32, 663–677. doi: 10.1007/s11032-013-9899-8

PubMed Abstract | Crossref Full Text | Google Scholar

Gaeta, R. T., Pires, J. C., Iniguez-Luy, F., Leon, E., and Osborn, T. C. (2007). Genomic changes in resynthesized brassica napus and their effect on gene expression and phenotype. Plant Cell 19, 3403–3417. doi: 10.1105/tpc.107.054346

PubMed Abstract | Crossref Full Text | Google Scholar

Hirani, A. H., Gao, F., Liu, J., Fu, G., Wu, C., McVetty, P. B., et al. (2018). Combinations of independent dominant loci conferring clubroot resistance in all four turnip accessions (Brassica rapa) from the European clubroot differential set. Front. Plant Sci. 9. doi: 10.3389/fpls.2018.01628

PubMed Abstract | Crossref Full Text | Google Scholar

Hittorf, M., Letsch Praxmarer, S., Windegger, A., Bass, D., Kirchmair, M., and Neuhauser, S. (2020). Revised taxonomy and expanded biodiversity of the Phytomyxea (Rhizaria, Endomyxa). J. Eukaryotic Microbiol. 67, 648–659. doi: 10.1111/jeu.12817

PubMed Abstract | Crossref Full Text | Google Scholar

Jabran, M., Ali, M. A., Zahoor, A., Muhae-Ud-Din, G., Liu, T., Chen, W., et al. (2023). Intelligent reprogramming of wheat for enhancement of fungal and nematode disease resistance using advanced molecular techniques. Front. Plant Sci. 14, 1132699. doi: 10.3389/fpls.2023.1132699

PubMed Abstract | Crossref Full Text | Google Scholar

Jubault, M., Lariagon, C., Taconnat, L., Renou, J., Gravot, A., Delourme, R., et al. (2013). Partial resistance to clubroot in Arabidopsis is based on changes in the host primary metabolism and targeted cell division and expansion capacity. Funct. Integr. Genomics 13, 191–205. doi: 10.1007/s10142-013-0312-9

PubMed Abstract | Crossref Full Text | Google Scholar

Kang, L., Qian, L., Zheng, M., Chen, L., Chen, H., Yang, L., et al. (2021). Genomic insights into the origin, domestication and diversification of Brassica juncea. Nat. Genet. 53, 1392–1402. doi: 10.1038/s41588-021-00922-y

PubMed Abstract | Crossref Full Text | Google Scholar

Li, L., Luo, Y., Chen, B., Xu, K., Zhang, F., Li, H., et al. (2016). A genome-wide association study reveals new loci for resistance to clubroot disease in brassica napus. Front. Plant Sci. 7. doi: 10.3389/fpls.2016.01483

PubMed Abstract | Crossref Full Text | Google Scholar

Limbalkar, O. M., Singh, R., Kumar, P., Nanjundan, J., Parihar, C. M., Vasisth, P., et al. (2021). Deployment of Brassica carinata. Braun derived Brassica juncea (L.) czern. lines for improving heterosis and water use efficiency under water deficit stress conditions. Front. Plant Science. 12. doi: 10.3389/fpls.2021.765645

PubMed Abstract | Crossref Full Text | Google Scholar

Liu, Y., Xu, A., Liang, F., Yao, X., Wang, Y., Liu, X., et al. (2018). Screening of clubroot-resistant varieties and transfer of clubroot resistance genes to Brassica napus using distant hybridization. Breed. Sci. 68, 258–267. doi: 10.1270/jsbbs.17125

PubMed Abstract | Crossref Full Text | Google Scholar

Ludwig-Müller, J., Auer, S., Jülke, S., and Marschollek, S. (2017). “Manipulation of auxin and cytokinin balance during the Plasmodiophora brassicae - Arabidopsis thaliana interaction,” in Auxins and cytokinins in plant biology: methods and protocols (New York, NY: Springer), 41–60.

PubMed Abstract | Google Scholar

Martin, A. and Jouve, N. (1992). “Cytogenetics of F1 and their progenies.“ in Distant Hybridization of Crop Plants (Berlin, Heidelberg: Springer), 82–105.

Google Scholar

McClintock, B. (1984). The significance of responses of the genome to challenge. Science 226, 792–801. doi: 10.1126/science.15739260

PubMed Abstract | Crossref Full Text | Google Scholar

Mei, J., Liu, J., Yue, F., Chen, Y., Ming, J., Xiong, Z., et al. (2022). Broadening the genetic base of Brassica juncea by introducing genomic components from B. rapa and B. nigra via digenomic allohexaploid bridging. Crop J. 10, 672–679. doi: 10.1016/j.cj.2021.09.008

Crossref Full Text | Google Scholar

Nannuru, V. K. R., Dieseth, J. A., Lillemo, M., and Meuwissen, T. H. E. (2025). Evaluating genomic selection and speed breeding for Fusarium head blight resistance in wheat using stochastic simulations. Molecular Breeding 45, 14. doi: 10.1007/s11032-024-01527-z

PubMed Abstract | Crossref Full Text | Google Scholar

Powell, W., Morgante, M., Andre, C., Hanafey, M., Vogel, J., Tingey, S., et al. (1996). The comparison of RFLP, RAPD, AFLP and SSR (microsatellite) markers for germplasm analysis. Mol. Breed. 2, 225–238. doi: 10.1007/BF00564200

Crossref Full Text | Google Scholar

Rahman, M., Shahadat, M. K., Munmun, T. S., Kobir, S., Maukeeb, A. R., Nasim, F. A., et al. (2022). Performance of mustard varieties under saline prone areas of Bangladesh. Afr. J. Agric. Res. 18, 608–616. doi: 10.5897/ajar2022.15955

Crossref Full Text | Google Scholar

Singh, K. H., Singh, L., Parmar, N., Kumar, S., Nanjundan, J., Singh, G., et al. (2022). Molecular characterization and genetic diversity analysis in Indian mustard (Brassica juncea L. Czern & Coss.) varieties using SSR markers. . PloS One 17, e0272914. doi: 10.1371/journal.pone.0272914

PubMed Abstract | Crossref Full Text | Google Scholar

Tian, Y. and Deng, F. (2020). Phytochemistry and biological activity of mustard (Brassica juncea): a review. CyTA-Journal Food 18, 704–718. doi: 10.1080/19476337.2020.1833988

Crossref Full Text | Google Scholar

Tian, Z., Ji, C., Xie, Z., Shi, X., Tian, B., Cao, G., et al. (2023). Integrated cytological and transcriptomic analysis reveals insights into pollen fertility in newly synthetic Brassica allohexaploids. Front. Plant Sci. 13. doi: 10.3389/fpls.2022.1096804

PubMed Abstract | Crossref Full Text | Google Scholar

Tomita, H., Shimizu, M., Asad-ud Doullah, M., Fujimoto, R., and Okazaki, K. (2013). Accumulation of quantitative trait loci conferring broad-spectrum clubroot resistance in Brassica oleracea. Mol. Breed. 32, 889–900. doi: 10.1007/s11032-013-9918-9

Crossref Full Text | Google Scholar

Xiong, Z., Gaeta, R. T., and Pires, J. C. (2011). Homoeologous shuffling and chromosome compensation maintain genome balance in resynthesized allopolyploid Brassica napus. Proc. Natl. Acad. Sci. United States America 108, 7908–7913. doi: 10.1073/pnas.1014138108

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, Z., Jiang, Y., Gong, J., Li, Q., Dun, B., Liu, D., et al. (2022b). R gene triplication confers European fodder turnip with improved clubroot resistance. Plant Biotechnol. J. 20, 1502–1517. doi: 10.1111/pbi.13827

PubMed Abstract | Crossref Full Text | Google Scholar

Yang, S., Zhang, K., Lu, C., Chen, G., Huang, Q., Ulhassan, Z., et al. (2022a). Genetic variation and stability analysis of an artificially synthesized allohexaploid Brassica for breeding innovations. Agronomy-Basel 12, 2843. doi: 10.3390/agronomy12112843

Crossref Full Text | Google Scholar

Zhang, L., Li, X., Chang, L., Wang, T., Liang, J., Lin, R., et al. (2022). Expanding the genetic variation of Brassica juncea by introgression of the Brassica rapa genome. Horticulture Res. 9, uhab054. doi: 10.1093/hr/uhab054

PubMed Abstract | Crossref Full Text | Google Scholar

Zhu, M., Yang, L., Zhang, Y., Zhuang, M., Ji, J., Hou, X., et al. (2022). Introgression of clubroot resistant gene into Brassica oleracea L. from Brassica rapa based on homoeologous exchange. Horticulture Res. 9, uhac195. doi: 10.1093/hr/uhac195/6679138

PubMed Abstract | Crossref Full Text | Google Scholar