Molecular cytogenetic identification of B genome chromosomes linked to blackleg disease resistance in Brassica napus × B. carinata interspecific hybrids

Brassica napus Haploid and Double Haploid Production and Its Latest Applications

Rapeseed is one of the most important oil crops in the world. Increasing demand for oil and limited agronomic capabilities of present-day rapeseed result in the need for rapid development of new, superior cultivars. Double haploid (DH) technology is a fast and convenient approach in plant breeding as well as genetic research. Brassica napus is considered a model species for DH production based on microspore embryogenesis; however, the molecular mechanisms underlying microspore reprogramming are still vague. It is known that morphological changes are accompanied by gene and protein expression patterns, alongside carbohydrate and lipid metabolism. Novel, more efficient methods for DH rapeseed production have been reported. This review covers new findings and advances in Brassica napus DH production as well as the latest reports related to agronomically important traits in molecular studies employing the double haploid rapeseed lines.

In silico integration of disease resistance QTL, genes and markers with the Brassica juncea physical map

Brassica juncea (AABB), Indian mustard, is a source of disease resistance genes for a wide range of pathogens. The availability of reference genome sequences for B. juncea has made it possible to characterise the genomic structure and distribution of these disease resistance genes. Potentially functional disease resistance genes can be identified by co-localization with genetically mapped disease resistance quantitative trait loci (QTL). Here we identify and characterise disease resistance gene analogs (RGAs), including nucleotide-binding site-leucine-rich repeat (NLR), receptor-like kinase (RLK) and receptor-like protein (RLP) classes, and investigate their association with disease resistance QTL intervals. The molecular genetic marker sequences for four white rust (Albugo candida) disease resistance QTL, six blackleg (Leptosphaeria maculans) disease resistance QTL and BjCHI1, a gene cloned from B. juncea for hypocotyl rot disease, were extracted from previously published studies and used to compare with candidate RGAs. Our results highlight the complications for the identification of functional resistance genes, including the duplicated appearance of genetic markers for several resistance loci, including Ac2(t), AcB1-A4.1, AcB1-A5.1, Rlm6 and PhR2 in both the A and B genomes, due to the presence of homoeologous regions. Furthermore, the white rust loci, Ac2(t) and AcB1-A4.1, mapped to the same position on chromosome A04 and may be different alleles of the same gene. Despite these challenges, a total of nine candidate genomic regions hosting 14 RLPs, 28 NLRs and 115 RLKs were identified. This study facilitates the mapping and cloning of functional resistance genes for applications in crop improvement programs.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01309-5.

Mining of Cloned Disease Resistance Gene Homologs (CDRHs) in Brassica Species and Arabidopsis thaliana

Simple Summary

Developing cultivars with resistance genes (R genes) is an effective strategy to support high yield and quality in Brassica crops. The availability of clone R gene and genomic sequences in Brassica species and Arabidopsis thaliana provide the opportunity to compare genomic regions and survey R genes across genomic databases. In this paper, we aim to identify genes related to cloned genes through sequence identity, providing a repertoire of species-wide related R genes in Brassica crops. The comprehensive list of candidate R genes can be used as a reference for functional analysis.

Abstract

Various diseases severely affect Brassica crops, leading to significant global yield losses and a reduction in crop quality. In this study, we used the complete protein sequences of 49 cloned resistance genes (R genes) that confer resistance to fungal and bacterial diseases known to impact species in the Brassicaceae family. Homology searches were carried out across Brassica napus, B. rapa, B. oleracea, B. nigra, B. juncea, B. carinata and Arabidopsis thaliana genomes. In total, 660 cloned disease R gene homologs (CDRHs) were identified across the seven species, including 431 resistance gene analogs (RGAs) (248 nucleotide binding site-leucine rich repeats (NLRs), 150 receptor-like protein kinases (RLKs) and 33 receptor-like proteins (RLPs)) and 229 non-RGAs. Based on the position and distribution of specific homologs in each of the species, we observed a total of 87 CDRH clusters composed of 36 NLR, 16 RLK and 3 RLP homogeneous clusters and 32 heterogeneous clusters. The CDRHs detected consistently across the seven species are candidates that can be investigated for broad-spectrum resistance, potentially providing resistance to multiple pathogens. The R genes identified in this study provide a novel resource for the future functional analysis and gene cloning of Brassicaceae R genes towards crop improvement.

Development of B. carinata with super-high erucic acid content through interspecific hybridization

KEY MESSAGE

Disomic alien chromosome addition Brassica carinata lines with super-high erucic acid content were developed through interspecific hybridization with B. juncea and characterized using molecular, cytological and biochemical techniques. Brassica carinata [A.] Braun (BBCC, 2n = 34) is a climate-resilient oilseed. Its seed oil is high in erucic acid (> 40%), rendering it well suited for the production of biofuel and other bio-based applications. To enhance the competitiveness of B. carinata with high erucic B. napus (HEAR), lines with super-high erucic acid content were developed through interspecific hybridization. To this end, a fad2B null allele from Brassica juncea (AABB, 2n = 36) was introgressed into B. carinata, resulting in a B. carinata fad2B mutant with erucic acid levels of over 50%. Subsequently, the FAE allele from B. rapa spp. yellow sarson (AA, 2n = 20) was transferred to the fad2B B. carinata line, yielding lines with erucic acid contents of up to 57.9%. Molecular analysis using the Brassica 90 K Illumina Infinium™ SNP genotyping array identified these lines as disomic alien chromosome addition lines, with two extra A08 chromosomes containing the BrFAE gene. The alien chromosomes from B. rapa were clearly distinguished by molecular cytogenetics in one of the addition lines. Analysis of microspore-derived offspring and hybrids from crosses with a CMS B. carinata line showed that the transfer rate of the A08 chromosome into male gametes was over 98%, resulting in almost completely stable transmission of an A08 chromosome copy into the progeny. The increase in erucic acid levels was accompanied by changes in the proportions of other fatty acids depending on the genetic changes that were introduced in the interspecific hybrids, providing valuable insights into erucic acid metabolism in Brassica.

Exploring the gene pool of Brassica napus by genomics-based approaches

De novo allopolyploidization in Brassica provides a very successful model for reconstructing polyploid genomes using progenitor species and relatives to broaden crop gene pools and understand genome evolution after polyploidy, interspecific hybridization and exotic introgression. B. napus (AACC), the major cultivated rapeseed species and the third largest oilseed crop in the world, is a young Brassica species with a limited genetic base resulting from its short history of domestication, cultivation, and intensive selection during breeding for target economic traits. However, the gene pool of B. napus has been significantly enriched in recent decades that has been benefit from worldwide effects by the successful introduction of abundant subgenomic variation and novel genomic variation via intraspecific, interspecific and intergeneric crosses. An important question in this respect is how to utilize such variation to breed crops adapted to the changing global climate. Here, we review the genetic diversity, genome structure, and population-level differentiation of the B. napus gene pool in relation to known exotic introgressions from various species of the Brassicaceae, especially those elucidated by recent genome-sequencing projects. We also summarize progress in gene cloning, trait-marker associations, gene editing, molecular marker-assisted selection and genome-wide prediction, and describe the challenges and opportunities of these techniques as molecular platforms to exploit novel genomic variation and their value in the rapeseed gene pool. Future progress will accelerate the creation and manipulation of genetic diversity with genomic-based improvement, as well as provide novel insights into the neo-domestication of polyploid crops with novel genetic diversity from reconstructed genomes.

Brassica carinata genome characterization clarifies U's triangle model of evolution and polyploidy in Brassica

Ethiopian mustard (Brassica carinata) in the Brassicaceae family possesses many excellent agronomic traits. Here, the high-quality genome sequence of B. carinata is reported. Characterization revealed a genome anchored to 17 chromosomes with a total length of 1.087 Gb and an N50 scaffold length of 60 Mb. Repetitive sequences account for approximately 634 Mb or 58.34% of the B. carinata genome. Notably, 51.91% of 97,149 genes are confined to the terminal 20% of chromosomes as a result of the expansion of repeats in pericentromeric regions. Brassica carinata shares one whole-genome triplication event with the five other species in U's triangle, a classic model of evolution and polyploidy in Brassica. Brassica carinata was deduced to have formed ∼0.047 Mya, which is slightly earlier than B. napus but later than B. juncea. Our analysis indicated that the relationship between the two subgenomes (BcaB and BcaC) is greater than that between other two tetraploid subgenomes (BjuB and BnaC) and their respective diploid parents. RNA-seq datasets and comparative genomic analysis were used to identify several key genes in pathways regulating disease resistance and glucosinolate metabolism. Further analyses revealed that genome triplication and tandem duplication played important roles in the expansion of those genes in Brassica species. With the genome sequencing of B. carinata completed, the genomes of all six Brassica species in U's triangle are now resolved. The data obtained from genome sequencing, transcriptome analysis, and comparative genomic efforts in this study provide valuable insights into the genome evolution of the six Brassica species in U's triangle.

Current Status of the Disease-Resistant Gene(s)/QTLs, and Strategies for Improvement in Brassica juncea

Brassica juncea is a major oilseed crop in tropical and subtropical countries, especially in south-east Asia like India, China, Bangladesh, and Pakistan. The widespread cultivation of genetically similar varieties tends to attract fungal pathogens which cause heavy yield losses in the absence of resistant sources. The conventional disease management techniques are often expensive, have limited efficacy, and cause additional harm to the environment. A substantial approach is to identify and use of resistance sources within the Brassica hosts and other non-hosts to ensure sustainable oilseed crop production. In the present review, we discuss six major fungal pathogens of B. juncea: Sclerotinia stem rot (Sclerotinia sclerotiorum), Alternaria blight (Alternaria brassicae), White rust (Albugo candida), Downy mildew (Hyaloperonospora parasitica), Powdery mildew (Erysiphe cruciferarum), and Blackleg (Leptoshaeria maculans). From discussing studies on pathogen prevalence in B. juncea, the review then focuses on highlighting the resistance sources and quantitative trait loci/gene identified so far from Brassicaceae and non-filial sources against these fungal pathogens. The problems in the identification of resistance sources for B. juncea concerning genome complexity in host subpopulation and pathotypes were addressed. Emphasis has been laid on more elaborate and coordinated research to identify and deploy R genes, robust techniques, and research materials. Examples of fully characterized genes conferring resistance have been discussed that can be transformed into B. juncea using advanced genomics tools. Lastly, effective strategies for B. juncea improvement through introgression of novel R genes, development of pre-breeding resistant lines, characterization of pathotypes, and defense-related secondary metabolites have been provided suggesting the plan for the development of resistant B. juncea.

Recent Findings Unravel Genes and Genetic Factors Underlying Leptosphaeria maculans Resistance in Brassica napus and Its Relatives

Among the Brassica oilseeds, canola (Brassica napus) is the most economically significant globally. However, its production can be limited by blackleg disease, caused by the fungal pathogen Lepstosphaeria maculans. The deployment of resistance genes has been implemented as one of the key strategies to manage the disease. Genetic resistance against blackleg comes in two forms: qualitative resistance, controlled by a single, major resistance gene (R gene), and quantitative resistance (QR), controlled by numerous, small effect loci. R-gene-mediated blackleg resistance has been extensively studied, wherein several genomic regions harbouring R genes against L. maculans have been identified and three of these genes were cloned. These studies advance our understanding of the mechanism of R gene and pathogen avirulence (Avr) gene interaction. Notably, these studies revealed a more complex interaction than originally thought. Advances in genomics help unravel these complexities, providing insights into the genes and genetic factors towards improving blackleg resistance. Here, we aim to discuss the existing R-gene-mediated resistance, make a summary of candidate R genes against the disease, and emphasise the role of players involved in the pathogenicity and resistance. The comprehensive result will allow breeders to improve resistance to L. maculans, thereby increasing yield.

Creation of fertility-restored materials for Ogura CMS in Brassica oleracea by introducing Rfo gene from Brassica napus via an allotriploid strategy

Ogura CMS fertility-restored materials, with 18 chromosomes, normal seed setting, stable fertility and closer genetic background to the parent Chinese kale, were successfully developed in B. oleracea via a triploid strategy for the first time. Ogura cytoplasmic male sterility (CMS) is the most widely used sterile type in seed production for commercial hybrids of Brassica oleracea vegetables. However, the natural Ogura CMS restorer line has not been found in B. oleracea crops. In this study, the triploid strategy was used with the aim to create euploid B. oleracea progenies with the Rfo gene. The allotriploid AAC hybrid YL2 was used as a male parent to backcross with Ogura CMS Chinese kale. After successive backcrosses, the BC₂ Rfo-positive individual 16CMSF2-11 and its BC₃ progenies, with 18 chromosomes, were developed, which were morphologically identical to the parent Chinese kale. Compared with F₁ and BC₁ plants, it showed stable fertility performance, and regular meiosis behavior and could produce seeds normally under natural pollination. The genomic composition analysis of Rfo-positive progenies by using molecular markers showed that more than 87% of the C-genome components of BC₃ Rfo-progenies recovered to the parent Chinese kale, while most or all of the A_n-genome segments were lost in 16CMSF2-11 and its progenies. The results suggested that the genetic background of Rfo-positive individuals was closer to that of the parent Chinese kale along with backcrossing. Hereof, the Ogura CMS fertility-restored materials of Chinese kale were successfully created via triploid strategy for the first time, providing a bridge for utilizing the Ogura CMS B. oleracea germplasm in the future. Moreover, our study indicates that the triploid strategy is effective for transferring genes from B. napus into B. oleracea.

Genome-Wide Association Mapping Identifies Novel Loci for Quantitative Resistance to Blackleg Disease in Canola

Blackleg disease, caused by the fungal pathogen Leptosphaeria maculans, continues to be a major concern for sustainable production of canola (Brassica napus L.) in many parts of the world. The deployment of effective quantitative resistance (QR) is recognized as a durable strategy in providing natural defense to pathogens. Herein, we uncover loci for resistance to blackleg in a genetically diverse panel of canola accessions by exploiting historic recombination events which occurred during domestication and selective breeding by genome-wide association analysis (GWAS). We found extensive variation in resistance to blackleg at the adult plant stage, including for upper canopy infection. Using the linkage disequilibrium and genetic relationship estimates from 12,414 high quality SNPs, GWAS identified 59 statistically significant and "suggestive" SNPs on 17 chromosomes of B. napus genome that underlie variation in resistance to blackleg, evaluated under field and shade-house conditions. Each of the SNP association accounted for up to 25.1% of additive genetic variance in resistance among diverse panel of accessions. To understand the homology of QR genomic regions with Arabidopsis thaliana genome, we searched the synteny between QR regions with 22 ancestral blocks of Brassicaceae. Comparative analyses revealed that 25 SNP associations for QR were localized in nine ancestral blocks, as a result of genomic rearrangements. We further showed that phenological traits such as flowering time, plant height, and maturity confound the genetic variation in resistance. Altogether, these findings provided new insights on the complex genetic control of the blackleg resistance and further expanded our understanding of its genetic architecture.

Genome-Wide Mapping of Loci Associated With Resistance to Clubroot in Brassica napus ssp. napobrassica (Rutabaga) Accessions From Nordic Countries

Rutabaga [Brassica napus ssp. napobrassica (L.) Hanelt] is reported to be an excellent source of clubroot (Plasmodiophora brassicae) resistance genes. In this study, 124 rutabaga accessions from the Nordic countries (Norway, Sweden, Finland, Denmark, and Iceland) were evaluated for their reaction to five single-spore isolates representing P. brassicae pathotypes 2F, 3H, 5I, 6M, and 8N and 12 field isolates representing pathotypes 2B, 3A, 3O, 5C, 5G, 5K, 5L, 5X (two isolates, L-G2 and L-G3), 8E, 8J, and 8P. The accessions were also genotyped using a 15K Brassica SNP array and 60 PCR-based primers linked to previously identified clubroot resistance genes. Six thousand eight hundred sixty-one SNP markers were retained after filtering with TASSEL 5.0, and used to evaluate four general linear models (GLM) and four mixed linear models (MLM). The PCA + K and Q + K MLM models gave the minimal deviance of the observed from the expected distribution in quantile-quantile plots, and hence were used for SNP-clubroot association analyses. In addition, 108 alleles derived from the PCR-based markers and the phenotypic data were analyzed with the PCA + K model. Forty-five SNPs and four PCR-based markers were identified to be associated strongly with resistance to isolates representing 13 pathotypes (2F, 3H, 5I, 6M, 8N, 2B, 3A, 3O, 5C, 5G, 5K, 5L, and 8P). These markers revealed the top and bottom segments of rutabaga chromosome A03 and the middle segment of chromosome A08 as genomic hotspots associated with resistance to the different P. brassicae pathotypes.

Multicolored Fluorescent In Situ Hybridization to Assess Pairing Configurations at Metaphase I in Brassica Hybrids

Genetic diversity can be introduced into polyploid crop species through meiotic recombination by exchanges between homologous or homoeologous chromosomes. Fluorescent in situ hybridization (FISH) enables the characterization of these homologous and homoeologous chromosome pairs during meiosis and identification of structural rearrangements during mitosis in metaphase I. In this chapter, we describe a protocol for the multicolored fluorescent labeling of chromosome spreads. This protocol allows the characterization of each A and C homoeologous subgenomes in a polyploid species using a genome-specific BAC combined with specific chromosome labeling BAC sequences.

Current Status and Challenges in Identifying Disease Resistance Genes in Brassica napus

Brassica napus is an economically important crop across different continents including temperate and subtropical regions in Europe, Canada, South Asia, China and Australia. Its widespread cultivation also brings setbacks as it plays host to fungal, oomycete and chytrid pathogens that can lead to serious yield loss. For sustainable crop production, identification of resistance (R) genes in B. napus has become of critical importance. In this review, we discuss four key pathogens affecting Brassica crops: Clubroot (Plasmodiophora brassicae), Blackleg (Leptosphaeria maculans and L. biglobosa), Sclerotinia Stem Rot (Sclerotinia sclerotiorum), and Downy Mildew (Hyaloperonospora parasitica). We first review current studies covering prevalence of these pathogens on Brassica crops and highlight the R genes and QTL that have been identified from Brassica species against these pathogens. Insights into the relationships between the pathogen and its Brassica host, the unique host resistance mechanisms and how these affect resistance outcomes is also presented. We discuss challenges in identification and deployment of R genes in B. napus in relation to highly specific genetic interactions between host subpopulations and pathogen pathotypes and emphasize the need for common or shared techniques and research materials or tighter collaboration between researchers to reconcile the inconsistencies in the research outcomes. Using current genomics tools, we provide examples of how characterization and cloning of R genes in B. napus can be carried out more effectively. Lastly, we put forward strategies to breed resistant cultivars through introgressions supported by genomic approaches and suggest prospects that can be implemented in the future for a better, pathogen-resistant B. napus.

Cytogenetic and Molecular Characterization of B-Genome Introgression Lines of Brassica napus L

Brassica napus introgression lines (ILs), having B-genome segments from B. carinata, were assessed genetically for extent of introgression and phenotypically for siliqua shatter resistance. Introgression lines had 7-9% higher DNA content, were meiotically stable, and had almost normal pollen fertility/seed set. Segment introgressions were confirmed by fluorescent genomic in situ hybridization (fl-GISH), SSR analyses, and SNP studies. Genotyping with 48 B-genome specific SSRs detected substitutions from B3, B4, B6, and B7 chromosomes on 39 of the 69 ILs whereas SNP genotyping detected a total of 23 B-segments (≥3 Mb) from B4, B6, and B7 introgressed into 10 of the 19 (C1, C2, C3, C5, C6, C8, C9, A3, A9, A10) chromosomes in 17 ILs. The size of substitutions varied from 3.0 Mb on chromosome A9 (IL59) to 42.44 Mb on chromosome C2 (IL54), ranging from 7 to 83% of the recipient chromosome. Average siliqua strength in ILs was observed to be higher than that of B. napus parents (2.2-6.0 vs. 1.9-4.0 mJ) while siliqua strength in some of the lines was almost equal to that of the donor parent B. carinata (6.0 vs.7.2 mJ). These ILs, with large chunks of substituted B-genome, can prove to be a useful prebreeding resource for germplasm enhancement in B. napus, especially for siliqua shatter resistance.

Introgression of Black Rot Resistance from Brassica carinata to Cauliflower (Brassica oleracea botrytis Group) through Embryo Rescue

Black rot caused by Xanthomonas campestris pv. campestris (Xcc) is a very important disease of cauliflower (Brassica oleracea botrytis group) resulting into 10-50% yield losses every year. Since there is a dearth of availability of resistance to black rot disease in B. oleracea (C genome), therefore exploration of A and B genomes was inevitable as they have been reported to be potential reservoirs of gene(s) for resistance to black rot. To utilize these sources, interspecific hybrid and backcross progeny (B₁) were generated between cauliflower "Pusa Sharad" and Ethiopian mustard "NPC-9" employing in vitro embryo rescue technique. Direct ovule culture method was better than siliqua culture under different temperature regime periods. Hybridity testing of F₁ inter-specific plants was carried out using co-dominant SSR marker and Brassica B and C genome-specific (DB and DC) primers. Meiosis in the di-genomic (BCC) interspecific hybrid of B. oleracea botrytis group (2n = 18, CC) × B. carinata (2n = 4x = 34, BBCC) was higly disorganized and cytological analysis of pollen mother cells revealed chromosomes 2n = 26 at metaphase-I. Fertile giant pollen grain formation was observed frequently in interspecific F₁ hybrid and BC₁ plants. The F₁ inter-specific plants were found to be resistant to Xcc race 1. Segregation distortion was observed in BC₁ generation for black rot resistance and different morphological traits. The At1g70610 marker analysis confirmed successful introgression of black rot resistance in interspecific BC₁ population. This effort will go a long way in pyramiding gene(s) for resistance against black rot in Cole crops, especially cauliflower and cabbage for developing durable resistance, thus minimize dependency on bactericides.

Resynthesis of Brassica napus through hybridization between B. juncea and B. carinata

First report for the resynthesis of Brassica napus by recombining A and C genome from B. juncea and B. carinata , respectively. Also documents B genome introgressions in resynthesized B. napus. Resynthesis of Brassica napus (AACC) was achieved by hybridizing Brassica juncea (AABB) with Brassica carinata (BBCC). This was facilitated by spontaneous chromosome doubling in the F1 hybrid (ABBC) to yield octaploid (AABBBBCC), elimination of extra B genome chromosomes in the resulting octaploid and in subsequent selfed generations, aided with directed selection for fertile plants having B. napus morphology. Twenty-five plants with varying degrees of resemblance to natural B. napus were identified from 17 A5 progenies and assayed for cytogenetic stability and genetic diversity. Majority of these plants, except six (2n = 38) were hyperploids (2n = 40-56). The six plants with 2n = 38 were designated as derived B. napus types. These showed an expected meiotic configuration of 19II at metaphase-I, with 19-19 distribution at anaphase-I. Genotyping based on A and C genome specific primers confirmed genetic identity of six derived (2n = 38) B. napus plants with natural types whereas genotyping with B genome specific primers indicated introgression of B genome segments. This was also confirmed by genomic in situ hybridization (GISH). Strong signals of B genome probe were detected, proving hitherto unreported genetic exchanges between B and A/C chromosomes. These introgressions possibly occurred en route five generations of selfing. Derived plants yielded fertile hybrids in crosses with natural B. napus var. GSC 6. The selfed derived plants as evaluated in A6 plant to progeny rows were morphologically similar to natural B. napus, and meiotically stable. Agronomic assessment of these progenies revealed variation for key morpho-physiological traits. Of special interest were the progenies with plants having oil content exceeding 47% as against about 39-41% in existing cultivars.

An RNA-seq transcriptome analysis of floral buds of an interspecific Brassica hybrid between B. carinata and B. napus

Interspecific hybridizations promote gene transfer between species and play an important role in plant speciation and crop improvement. However, hybrid sterility that commonly found in the first generation of hybrids hinders the utilization of interspecific hybridization. The combination of divergent parental genomes can create extensive transcriptome variations, and to determine these gene expression alterations and their effects on hybrids, an interspecific Brassica hybrid of B. carinata × B. napus was generated. Scanning electron microscopy analysis indicated that some of the hybrid pollen grains were irregular in shape and exhibited abnormal exine patterns compared with those from the parents. Using the Illumina HiSeq 2000 platform, 39,598, 32,403 and 42,208 genes were identified in flower buds of B. carinata cv. W29, B. napus cv. Zhongshuang 11 and their hybrids, respectively. The differentially expressed genes were significantly enriched in pollen wall assembly, pollen exine formation, pollen development, pollen tube growth, pollination, gene transcription, macromolecule methylation and translation, which might be associated with impaired fertility in the F1 hybrid. These results will shed light on the mechanisms underlying the low fertility of the interspecific hybrids and expand our knowledge of interspecific hybridization.