Mapping genes for resistance to Leptosphaeria maculans in Brassica juncea

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.

Candidate Rlm6 resistance genes against Leptosphaeria. maculans identified through a genome-wide association study in Brassica juncea (L.) Czern

One hundred and sixty-seven B. juncea varieties were genotyped on the 90K Brassica assay (42,914 SNPs), which led to the identification of sixteen candidate genes for Rlm6. Brassica species are at high risk of severe crop loss due to pathogens, especially Leptosphaeria maculans (the causal agent of blackleg). Brassica juncea (L.) Czern is an important germplasm resource for canola improvement, due to its good agronomic traits, such as heat and drought tolerance and high blackleg resistance. The present study is the first using genome-wide association studies to identify candidate genes for blackleg resistance in B. juncea based on genome-wide SNPs obtained from the Illumina Infinium 90 K Brassica SNP array. The verification of Rlm6 in B. juncea was performed through a cotyledon infection test. Genotyping 42,914 single nucleotide polymorphisms (SNPs) in a panel of 167 B. juncea lines revealed a total of seven SNPs significantly associated with Rlm6 on chromosomes A07 and B04 in B. juncea. Furthermore, 16 candidate Rlm6 genes were found in these regions, defined as nucleotide binding site leucine-rich-repeat (NLR), leucine-rich repeat RLK (LRR-RLK) and LRR-RLP genes. This study will give insights into the blackleg resistance in B. juncea and facilitate identification of functional blackleg resistance genes which can be used in Brassica breeding.

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.

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.

Identification of resistance loci in Chinese and Canadian canola/rapeseed varieties against Leptosphaeria maculans based on genome-wide association studies

Background

The fungal pathogen Leptosphaeria maculans (Lm). causes blackleg disease on canola/rapeseed in many parts of the world. It is important to use resistant cultivars to manage the disease and minimize yield losses. In this study, twenty-two Lm isolates were used to identify resistance genes in a collection of 243 canola/rapeseed (Brassica napus L.) accessions from Canada and China. These Lm isolates carry different compliments of avirulence genes, and the investigation was based on a genome-wide association study (GWAS) and genotype-by-sequencing (GBS).

Results

Using the CROP-SNP pipeline, a total of 81,471 variants, including 78,632 SNPs and 2839 InDels, were identified. The GWAS was performed using TASSEL 5.0 with GLM + Q model. Thirty-two and 13 SNPs were identified from the Canadian and Chinese accessions, respectively, tightly associated with blackleg resistance with P values < 1 × 10^− 4. These SNP loci were distributed on chromosomes A03, A05, A08, A09, C01, C04, C05, and C07, with the majority of them on A08 followed by A09 and A03. The significant SNPs identified on A08 were all located in a 2010-kb region and associated with resistance to 12 of the 22 Lm isolates. Furthermore, 25 resistance gene analogues (RGAs) were identified in these regions, including two nucleotide binding site (NBS) domain proteins, fourteen RLKs, three RLPs and six TM-CCs. These RGAs can be the potential candidate genes for blackleg resistance.

Conclusion

This study provides insights into potentially new genomic regions for discovery of additional blackleg resistance genes. The identified regions associated with blackleg resistance in the germplasm collection may also contribute directly to the development of canola varieties with novel resistance genes against blackleg of canola.

An update on the arsenal: mining resistance genes for disease management of Brassica crops in the genomic era

Brassica species include many economically important crops that provide nutrition and health-promoting substances to humans worldwide. However, as with all crops, their production is constantly threatened by emerging viral, bacterial, and fungal diseases, whose incidence has increased in recent years. Traditional methods of control are often costly, present limited effectiveness, and cause environmental damage; instead, the ideal approach is to mine and utilize the resistance genes of the Brassica crop hosts themselves. Fortunately, the development of genomics, molecular genetics, and biological techniques enables us to rapidly discover and apply resistance (R) genes. Herein, the R genes identified in Brassica crops are summarized, including their mapping and cloning, possible molecular mechanisms, and application in resistance breeding. Future perspectives concerning how to accurately discover additional R gene resources and efficiently utilize these genes in the genomic era are also discussed.

Inheritance Pattern and Molecular Markers for Resistance to Blackleg Disease in Cabbage

The inheritance and causal loci for resistance to blackleg, a devastating disease of Brassicaceous crops, are yet to be known in cabbage (Brassica oleracea L.). Here, we report the pattern of inheritance and linked molecular marker for this trait. A segregating BC₁ population consisting of 253 plants was raised from resistant and susceptible parents, L29 (♀) and L16 (♂), respectively. Cotyledon resistance bioassay of BC₁ population, measured based on a scale of 0-9 at 12 days after inoculation with Leptosphaeria maculans isolate 03-02 s, revealed the segregation of resistance and ratio, indicative of dominant monogenic control of the trait. Investigation of potential polymorphism in the previously identified differentially expressed genes within the collinear region of 'B. napus blackleg resistant loci Rlm1' in B. oleracea identified two insertion/deletion (InDel) mutations in the intron and numerous single nucleotide polymorphisms (SNPs) throughout the LRR-RLK gene Bol040029, of which six SNPs in the first exon caused the loss of two LRR domains in the susceptible line. An InDel marker, BLR-C-InDel based on the InDel mutations, and a high resolution melting (HRM) marker, BLR-C-2808 based on the SNP C²⁸⁰⁸T in the second exon were developed, which predicated the resistance status of the BC₁ population with 80.24%, and of 24 commercial inbred lines with 100% detection accuracy. This is the first report of inheritance and molecular markers linked with blackleg resistance in cabbage. This study will enhance our understanding of the trait, and will be helpful in marker assisted breeding aiming at developing resistant cabbage varieties.

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.

Single R Gene Introgression Lines for Accurate Dissection of the Brassica - Leptosphaeria Pathosystem

Seven blackleg resistance (R) genes (Rlm1, Rlm2, Rlm3, Rlm4, LepR1, LepR2 & LepR3) were each introgressed into a common susceptible B. napus doubled-haploid (DH) line through reciprocal back-crossing, producing single-R gene introgression lines (ILs) for use in the pathological and molecular study of Brassica-Leptosphaeria interactions. The genomic positions of the R genes were defined through molecular mapping and analysis with transgenic L. maculans isolates was used to confirm the identity of the introgressed genes where possible. Using L. maculans isolates of contrasting avirulence gene (Avr) profiles, we preformed extensive differential pathology for phenotypic comparison of the ILs to other B. napus varieties, demonstrating the ILs can provide for the accurate assessment of Avr-R gene interactions by avoiding non-Avr dependant alterations to resistance responses which can occur in some commonly used B. napus varieties. Whole-genome SNP-based assessment allowed us to define the donor parent introgressions in each IL and provide a strong basis for comparative molecular dissection of the pathosystem.

Co-linearity and divergence of the A subgenome of Brassica juncea compared with other Brassica species carrying different A subgenomes

Background

There are three basic Brassica genomes (A, B, and C) and three parallel sets of subgenomes distinguished in the diploid Brassica (i.e.: B. rapa, A^rA^r; B. nigra, BⁿⁱBⁿⁱ; B. oleracea, C^oC^o) and the derived allotetraploid species (i.e.: B. juncea, A^jA^jB^jB^j; B. napus, AⁿAⁿCⁿCⁿ; B. carinata, B^cB^cC^cC^c). To understand subgenome differentiation in B. juncea in comparison to other A genome-carrying Brassica species (B. rapa and B. napus), we constructed a dense genetic linkage map of B. juncea, and conducted population genetic analysis on diverse lines of the three A-genome carrying Brassica species using a genotyping-by-sequencing approach (DArT-seq).

Results

A dense genetic linkage map of B. juncea was constructed using an F₂ population derived from Sichuan Yellow/Purple Mustard. The map included 3329 DArT-seq markers on 18 linkage groups and covered 1579 cM with an average density of two markers per cM. Based on this map and the alignment of the marker sequences with the physical genome of Arabidopsis thaliana, we observed strong co-linearity of the ancestral blocks among the different A subgenomes but also considerable block variation. Comparative analyses at the level of genome sequences of B. rapa and B. napus, and marker sequence anchored on the genetic map of B. juncea, revealed a total of 30 potential inversion events across large segments and 20 potential translocation events among the three A subgenomes. Population genetic analysis on 26 accessions of the three A genome-carrying Brassica species showed that the highest genetic distance were estimated when comparing A^j-Aⁿ than between Aⁿ-A^r and A^j-A^r subgenome pairs.

Conclusions

The development of the dense genetic linkage map of B. juncea with informative DArT-seq marker sequences and availability of the reference sequences of the A^r, and AⁿCⁿ genomes allowed us to compare the A subgenome structure of B. juncea (A^j) . Our results suggest that strong co-linearity exists among the three A Brassica genomes (A^r, Aⁿ and A^j) but with apparent subgenomic variation. Population genetic analysis on three A-genome carrying Brassica species support the idea that B. juncea has distinct genomic diversity, and/or evolved from a different A genome progenitor of B. napus.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2343-1) contains supplementary material, which is available to authorized users.

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

Provide evidence that the Brassica B genome chromosome B3 carries blackleg resistance gene, and also the B genome chromosomes were inherited several generations along with B. napus chromosomes. Blackleg disease caused by fungus Leptosphaeria maculans causes significant yield losses in Brassica napus. Brassica carinata possesses excellent resistance to this disease. To introgress blackleg resistance, crosses between B. napus cv. Westar and B. carinata were done. The interspecific-hybrids were backcrossed twice to Westar and self-pollinated three times to produce BC2S3 families. Doubled haploid lines (DH1) were produced from one blackleg resistant family. SSR markers were used to study the association between B genome chromosome(s) and blackleg resistance. The entire B3 chromosome of B. carinata was associated with blackleg resistance in DH1. A second DH population (DH2) was produced from F1s of resistant DH1 lines crossed to blackleg susceptible B. napus cv. Polo where resistance was found to be associated with SSR markers from the middle to bottom of the B3 and top of the B8 chromosomes. The results demonstrated that the B3 chromosome carried gene(s) for blackleg resistance. Genomic in situ hybridization (GISH) and GISH-like analysis of the DH2 lines revealed that susceptible lines, in addition to B. napus chromosomes, possessed one pair of B genome chromosomes (2n = 40), while resistant lines had either one (2n = 40) or two pairs (2n = 42) of B chromosomes. The molecular and GISH data suggested that the B chromosome in the susceptible lines was B7, while it was difficult to confirm the identity of the B chromosomes in the resistant lines. Also, B chromosomes were found to be inherited over several generations along with B. napus chromosomes.

An avirulence gene, AvrLmJ1, from the blackleg fungus, Leptosphaeria maculans, confers avirulence to Brassica juncea cultivars

The fungus Leptosphaeria maculans causes blackleg of Brassica species. Here, we report the mapping and subsequent cloning of an avirulence gene from L. maculans. This gene, termed AvrLmJ1, confers avirulence towards all three Brassica juncea cultivars tested. Analysis of RNA-seq data showed that AvrLmJ1 is housed in a region of the L. maculans genome which contains only one gene that is highly expressed in planta. The closest genes are 57 and 33 kb away and, like other avirulence genes of L. maculans, AvrLmJ1 is located within an AT-rich, gene-poor region of the genome. The encoded protein is 141 amino acids, has a predicted signal peptide and is cysteine rich. Two virulent isolates contain a premature stop codon in AvrLmJ1. Complementation of an isolate that forms cotyledonary lesions on B. juncea with the wild-type allele of AvrLmJ1 confers avirulence towards all three B. juncea cultivars tested, suggesting that the gene may confer species-specific avirulence activity.

Comparative genomics of Brassicaceae crops

The family Brassicaceae is one of the major groups of the plant kingdom and comprises diverse species of great economic, agronomic and scientific importance, including the model plant Arabidopsis. The sequencing of the Arabidopsis genome has revolutionized our knowledge in the field of plant biology and provides a foundation in genomics and comparative biology. Genomic resources have been utilized in Brassica for diversity analyses, construction of genetic maps and identification of agronomic traits. In Brassicaceae, comparative sequence analysis across the species has been utilized to understand genome structure, evolution and the detection of conserved genomic segments. In this review, we focus on the progress made in genetic resource development, genome sequencing and comparative mapping in Brassica and related species. The utilization of genomic resources and next-generation sequencing approaches in improvement of Brassica crops is also discussed.

Identification and mapping of a novel blackleg resistance locus LepR4 in the progenies from Brassica napus × B. rapa subsp. sylvestris

Blackleg, caused by Leptosphaeria maculans, is one of the most economically important diseases of Brassica napus worldwide. Two blackleg-resistant lines, 16S and 61446, were developed through interspecific hybridization between B. napus and B. rapa subsp. sylvestris and backcrossing to B. napus. Classical genetic analysis demonstrated that a single recessive gene in both lines conferred resistance to L. maculans and that the resistance alleles were allelic. Using BC(1) progeny derived from each resistant plant, this locus was mapped to B. napus linkage group N6 and was flanked by microsatellite markers sN2189b and sORH72a in an interval of about 10 cM, in a region equivalent to about 6 Mb of B. rapa DNA sequence. This new resistance gene locus was designated as LepR4. The two lines were evaluated for resistance to a wide range of L. maculans isolates using cotyledon inoculation tests under controlled environment conditions, and for stem canker resistance in blackleg field nurseries. Results indicated that line 16S, carrying LepR4a, was highly resistant to all isolates tested on cotyledons and had a high level of stem canker resistance under field conditions. Line 61446, carrying LepR4b, was only resistant to some of the isolates tested on cotyledons and was weakly resistant to stem canker under field conditions.

The Brassica napus blackleg resistance gene LepR3 encodes a receptor-like protein triggered by the Leptosphaeria maculans effector AVRLM1

LepR3, found in the Brassica napus cv 'Surpass 400', provides race-specific resistance to the fungal pathogen Leptosphaeria maculans, which was overcome after great devastation in Australia in 2004. We investigated the LepR3 locus to identify the genetic basis of this resistance interaction. We employed a map-based cloning strategy, exploiting collinearity with the Arabidopsis thaliana and Brassica rapa genomes to enrich the map and locate a candidate gene. We also investigated the interaction of LepR3 with the L. maculans avirulence gene AvrLm1 using transgenics. LepR3 was found to encode a receptor-like protein (RLP). We also demonstrated that avirulence towards LepR3 is conferred by AvrLm1, which is responsible for both the Rlm1 and LepR3-dependent resistance responses in B. napus. LepR3 is the first functional B. napus disease resistance gene to be cloned. AvrLm1's interaction with two independent resistance loci, Rlm1 and LepR3, highlights the need to consider redundant phenotypes in 'gene-for-gene' interactions and offers an explanation as to why LepR3 was overcome so rapidly in parts of Australia.

Introgression of B-genome chromosomes in a doubled haploid population of Brassica napus x B. carinata

The Brassica B-genome species possess many valuable agronomic and disease resistance traits. To transfer traits from the B genome of B. carinata into B. napus, an interspecific cross between B. napus and B. carinata was performed and a doubled haploid (DH) population was generated from the BC2S3 generation. Successful production of interspecific DH lines as identified using B-genome microsatellite markers is reported. Five percent of DH lines carry either intact B-genome chromosomes or chromosomes that have deletions. All of the DH lines have linkage group J13/B7 in common. This was further confirmed using B. nigra genomic DNA in a fluorescent in situ hybridization assay where the B-genome chromosomes were visualized and distinguished from the A- and C-genome chromosomes. The 60 DH lines were also evaluated for morphological traits in the field for two seasons and were tested for resistance to blackleg, caused by Leptosphaeria maculans, under greenhouse conditions. Variation in the DH population followed a normal distribution for several agronomic traits and response to blackleg. The lines with B-genome chromosomes were significantly different (p < 0.01) from the lines without B-genome chromosomes for both morphological and seed quality traits such as days to flowering, days to maturity, and erucic acid content.