RGA- and RAPD-derived SCAR markers for a Brassica B-genome introgression conferring resistance to blackleg in oilseed rape

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.

The Rlm13 Gene, a New Player of Brassica napus-Leptosphaeria maculans Interaction Maps on Chromosome C03 in Canola

Canola exhibits an extensive genetic variation for resistance to blackleg disease, caused by the fungal pathogen Leptosphaeria maculans. Despite the identification of several Avr effectors and R (race-specific) genes, specific interactions between Avr-R genes are not yet fully understood in the Brassica napus-L. maculans pathosystem. In this study, we investigated the genetic basis of resistance in an F_{2 : 3} population derived from Australian canola varieties CB-Telfer (Rlm4)/ATR-Cobbler (Rlm4) using a single-spore isolate of L. maculans, PHW1223. A genetic linkage map of the CB-Telfer/ATR-Cobbler population was constructed using 7,932 genotyping-by-sequencing-based DArTseq markers and subsequently utilized for linkage and haplotype analyses. Genetic linkage between DArTseq markers and resistance to PHW1223 isolate was also validated using the B. napus 60K Illumina Infinium array. Our results revealed that a major locus for resistance, designated as Rlm13, maps on chromosome C03. To date, no R gene for resistance to blackleg has been reported on the C subgenome in B. napus. Twenty-four candidate R genes were predicted to reside within the quantitative trait locus (QTL) region. We further resequenced both the parental lines of the mapping population (CB-Telfer and ATR-Cobbler, > 80 × coverage) and identified several structural sequence variants in the form of single-nucleotide polymorphisms (SNPs), insertions/deletions (InDels), and presence/absence variations (PAVs) near Rlm13. Comparative mapping revealed that Rlm13 is located within the homoeologous A03/C03 region in ancestral karyotype block "R" of Brassicaceae. Our results provide a "target" for further understanding the Avr-Rlm13 gene interaction as well as a valuable tool for increasing resistance to blackleg in canola germplasm.

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.

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.

A user guide to the Brassica 60K Illumina Infinium™ SNP genotyping array

The Brassica napus 60K Illumina Infinium™ SNP array has had huge international uptake in the rapeseed community due to the revolutionary speed of acquisition and ease of analysis of this high-throughput genotyping data, particularly when coupled with the newly available reference genome sequence. However, further utilization of this valuable resource can be optimized by better understanding the promises and pitfalls of SNP arrays. We outline how best to analyze Brassica SNP marker array data for diverse applications, including linkage and association mapping, genetic diversity and genomic introgression studies. We present data on which SNPs are locus-specific in winter, semi-winter and spring B. napus germplasm pools, rather than amplifying both an A-genome and a C-genome locus or multiple loci. Common issues that arise when analyzing array data will be discussed, particularly those unique to SNP markers and how to deal with these for practical applications in Brassica breeding applications.

Cloning of TTG1 gene and PCR identification of genomes A, B and C in Brassica species

Arabidopsis Transparent Testa Glabra 1 (TTG1) genes were cloned from three diploid Brassica species (B. rapa, B. nigra and B. oleracea) and two amphidiploids species (B. juncea and B. carinata) by homology cloning. TTG1 homologues identified in all the accessions of the investigated species had a coding sequence of 1,014 bp. One copy was obtained from each diploid species and two copies from each amphidiploid species. Combined analysis of the TTG1 sequences cloned in this study with those obtained from public databases demonstrated that three, forty-five and seven nucleotides were specific variations in TTG1 genes from genomes A, B and C, respectively. Primers designed with genome-specific nucleotide variations were able to distinguish among TTG1 genes originating from genomes A, B and C in Brassica. Therefore, the TTG1 gene could serve as a candidate marker gene to detect the pollen flow of Brassica and provide an alternative method for the detection of pollen drift and risk assessment of gene flow in Brassica species.

Conserved microstructure of the Brassica B Genome of Brassica nigra in relation to homologous regions of Arabidopsis thaliana, B. rapa and B. oleracea

BACKGROUND

The Brassica B genome is known to carry several important traits, yet there has been limited analyses of its underlying genome structure, especially in comparison to the closely related A and C genomes. A bacterial artificial chromosome (BAC) library of Brassica nigra was developed and screened with 17 genes from a 222 kb region of A. thaliana that had been well characterised in both the Brassica A and C genomes.

RESULTS

Fingerprinting of 483 apparently non-redundant clones defined physical contigs for the corresponding regions in B. nigra. The target region is duplicated in A. thaliana and six homologous contigs were found in B. nigra resulting from the whole genome triplication event shared by the Brassiceae tribe. BACs representative of each region were sequenced to elucidate the level of microscale rearrangements across the Brassica species divide.

CONCLUSIONS

Although the B genome species separated from the A/C lineage some 6 Mya, comparisons between the three paleopolyploid Brassica genomes revealed extensive conservation of gene content and sequence identity. The level of fractionation or gene loss varied across genomes and genomic regions; however, the greatest loss of genes was observed to be common to all three genomes. One large-scale chromosomal rearrangement differentiated the B genome suggesting such events could contribute to the lack of recombination observed between B genome species and those of the closely related A/C lineage.

Genetic linkage map construction and QTL mapping of cadmium accumulation in radish (Raphanus sativus L.)

Cadmium (Cd) is a widespread soil pollutant and poses a significant threat to human health via the food chain. Large phenotypic variations in Cd concentration of radish roots and shoots have been observed. However, the genetic and molecular mechanisms of Cd accumulation in radish remain to be elucidated. In this study, a genetic linkage map was constructed using an F(2) mapping population derived from a cross between a high Cd-accumulating cultivar NAU-Dysx and a low Cd-accumulating cultivar NAU-Yh. The linkage map consisted of 523 SRAP, RAPD, SSR, ISSR, RAMP, and RGA markers and had a total length of 1,678.2 cM with a mean distance of 3.4 cM between two markers. All mapped markers distributed on nine linkage groups (LGs) having sizes between 134.7 and 236.8 cM. Four quantitative trait loci (QTLs) for root Cd accumulation were mapped on LGs 1, 4, 6, and 9, which accounted for 9.86 to 48.64 % of all phenotypic variance. Two QTLs associated with shoot Cd accumulation were detected on LG1 and 3, which accounted for 17.08 and 29.53 % of phenotypic variance, respectively. A major-effect QTL, qRCd9 (QTL for root Cd accumulation on LG9), was identified on LG 9 flanked by NAUrp011_754 and EM5me6_286 markers with a high LOD value of 23.6, which accounted for 48.64 % of the total phenotypic variance in Cd accumulation of F(2) lines. The results indicated that qRCd9 is a novel QTL responsible for controlling root Cd accumulation in radish, and the identification of specific molecular markers tightly linked to the major QTL could be further applied for marker-assisted selection (MAS) in low-Cd content radish breeding program.

Genomic advances will herald new insights into the Brassica: Leptosphaeria maculans pathosystem

The study of the relationship between plants and phytopathogenic fungi is one of the most rapidly moving fields in the plant sciences, the findings of which have contributed to the development of new strategies and technologies to protect crops. Plants employ sophisticated mechanisms to perceive and appropriately defend themselves against pathogens. A good example of plant and pathogen evolution is the gene-for-gene interaction between the fungal pathogen Leptosphaeria maculans, the causal agent of blackleg disease, and Brassica crops. This interaction has been studied at the genetic and physiological level due to its agro-economic importance. The newly available genome sequence for Brassica spp. and L. maculans will provide the resources to study the co-evolution of this plant and pathogen. Particularly, an understanding of the co-evolution of genes responsible for virulence and resistance will lead to improved plant protection strategies for Brassica canola and provide a model to understand plant-pathogen interactions in other major crops. This review summarises the research-to-date in the study of the Brassica-L. maculans gene-for-gene interaction, with a focus on the genetics of resistance in Brassica and the wealth of information to be gained from genome sequencing efforts.

Production and characterization of interspecific somatic hybrids between Brassica oleracea var. botrytis and B. nigra and their progenies for the selection of advanced pre-breeding materials

Somatic hybridization is a potential method for gene transfer from wild relatives to cultivated crops that can overcome sexual incompatibilities of two distantly related species. In this study, interspecific asymmetric somatic hybrids of Brassica oleracea var. botrytis (cauliflower) and Brassica nigra (black mustard) were obtained by protoplast fusion and their backcrossed (BC(3)) and selfed (S(3)) offspring were analyzed. Cytological analysis showed that the B. nigra chromosomes were successively eliminated in the backcrosses with cauliflower. The fertility of the hybrid progenies was quite different due to the asynchronous and abnormal chromosome behavior of pollen mother cells (PMC) during meiosis. Analysis of sequence-related amplified polymorphism (SRAP) showed that all of these hybrids mainly had the DNA banding pattern from the two parents with some alterations. Genetically, the selfed generations were closer to B. nigra, while the backcrossed generations were closer to the cauliflower parent. Analysis of cleaved amplified polymorphic sequences (CAPS) and restriction fragment length polymorphisms (RFLP) showed that all somatic hybrids in this study contained chloroplast (cp) DNA of the donor parent black mustard, while mitochondrial (mt) DNA showed evidence of recombination and variations in the regions analyzed. Furthermore, three BC(3) plants (originated from somatic hybrids 3, 4, 10) with 2-8 B. nigra-derived chromosomes shown by genomic in situ hybridization (GISH) displayed a more cauliflower-like morphology and high resistance to black-rot. These plants were obtained as bridge materials for further analysis and breeding.

RLM3, a potential adaptor between specific TIR-NB-LRR receptors and DZC proteins

In our recent paper, we identified a TIR encoding gene, which is required for resistance against a broad range of necrotrophic fungi. Here we present this finding in a broader perspective and discuss the unique features of this gene which might explain its role as a general regulator of resistance responses against a class of pathogens that have previously not been associated to the classical resistance (R) gene type of defense.

Segregation distortion detected in six rice F2 populations generated from reciprocal hybrids at three altitudes

This paper presents investigations of segregation distortion of six rice F2 populations generated from reciprocal F1 hybrids grown at three locations varied at altitudes from 400 to 2200 m. The F1s were derived from reciprocal crosses between cv. XMG, which is a japonica landrace traditionally grown at 2650 m altitude, and cv. N34, which is a japonica restorer possessing a fertility restoring (Rf) gene and cytoplasm of male sterility (CMS) donated by an indica cultivar. Among nine morphological traits of the F2 populations, only one was in normal distribution, eight were distorted in all or at least one population. Out of 16 polymorphic PCR markers, 10 markers distributed on 7 chromosomes were significantly distorted. Among these markers, RMAN7 and RM257 were distorted in both of the reciprocal populations, which suggested that nuclear genes had strong effects on segregation distortion. The other makers were distorted only in the populations with cytoplasm donated by XMG or N34. The results indicated that segregation of DNA markers was affected by cytoplasm background. Segregation distribution was also affected by altitude, since segregation distortions of most of the markers were detected not in all the three populations generated from F1 grown at the three altitudes, but only in one population from F1 grown at one altitude. Marker M45461, which is located within Rf-1 locus, was severely distorted towards N34 in all the populations with cytoplasm donated by N34, but not in the populations with cytoplasm provided by XMG. The results indicated that interaction between CMS and Rf gene had strong effects on distortion. Results of this study indicated that japonica cytoplasm did not cause distortion favouring a special parent, but indica cytoplasm made distortion favouring a maternal parent. The results suggested that indica cytoplasm was not well compatible with japonica nuclear background, while japonica cytoplasm did not have such trouble with indica nuclei. This study also found that the six F2 populations were divergent into two groups due to difference of cytoplasm background.

Layers of defense responses to Leptosphaeria maculans below the RLM1- and camalexin-dependent resistances

Plants have evolved different defense components to counteract pathogen attacks. The resistance locus resistance to Leptosphaeria maculans 1 (RLM1) is a key factor for Arabidopsis thaliana resistance to L. maculans. The present work aimed to reveal downstream defense responses regulated by RLM1. Quantitative assessment of fungal colonization in the host was carried out using quantitative polymerase chain reaction (qPCR) and GUS expression analyses, to further characterize RLM1 resistance and the role of salicylic acid (SA), jasmonic acid (JA), and ethylene (ET) in disease development. Additional assessments of A. thaliana mutants were performed to expand our understanding of this pathosystem. Resistance responses such as lignification and the formation of vascular plugs were found to occur in an RLM1-dependent manner, in contrast to the RLM1-independent increase in reactive oxygen species at the stomata and hydathodes. Analyses of mutants defective in hormone signaling in the camalexin-free rlm1(Ler)pad3 background revealed a significant influence of JA and ET on symptom development and pathogen colonization. The overall results indicate that the defense responses of primary importance induced by RLM1 are all associated with physical barriers, and that responses of secondary importance involve complex cross-talk among SA, JA and ET. Our observations further suggest that ET positively affects fungal colonization.

Stabilization of Resistance to Leptosphaeria maculans in Brassica napus-B. juncea Recombinant Lines and Its Introgression into Spring-Type Brassica napus

The value of Katanning Early Maturing (KEM) breeding lines from Western Australia, derived from Brassica napus × B. juncea crosses, was assessed as a source of germplasm for resistance to blackleg disease (caused by Leptosphaeria maculans) in spring-type oilseed rape cultivars. The stability of blackleg resistance in these KEM lines was related to key cytological characteristics to determine why there are poor levels of introgression of this resistance into progeny. Promising recombinant KEM lines were crossed with the spring-type B. napus cv. Dunkeld, which has useful polygenic resistance to blackleg, and screened for resistance. The lines were analyzed cytologically for pairing of bivalents in each generation to aid in the selection of stable recombinant lines. KEM recombinant lines showing regular meiotic behavior and a high level of blackleg resistance were obtained for the first time. We also showed that the stable introgression of the B. juncea resistance from the KEM lines into a 'Dunkeld' background was possible. Inoculation of selfing and backcross populations with isolates of L. maculans having different AvrLm genes indicated that the B. juncea resistance gene, Rlm6, had been introgressed into a B. napus spring-type cultivar carrying polygenic resistance. The combination of both resistances would enhance the overall effectiveness of resistance against L. maculans. This is clearly needed in Australia and France where cultivars relying upon single dominant gene-based resistance for their effectiveness have proved not durable.

RLM3, a TIR domain encoding gene involved in broad-range immunity of Arabidopsis to necrotrophic fungal pathogens

Here, we describe the rapid cloning of a plant gene, Leptosphaeria maculans 3 (RLM3(Col)), which encodes a putative Toll interleukin-1 receptor-nucleotide binding (TIR-NB) class protein, which is involved in defence against the fungal pathogen L. maculans and against three other necrotrophic fungi. We have, through microarray-based case control bulk segregant comparisons of transcriptomes in pools of Col-0 x An-1 progeny, identified the absence of a locus that causes susceptibility in An-1. The significance of this locus on chromosome 4 for L. maculans resistance was supported by PCR-based mapping, and denoted resistance to RLM3(Col). Differential susceptible phenotypes in four independent T-DNA insertion lines support the hypothesis that At4g16990 is required for RLM3(Col) function. The mutants in RLM3(Col) also exhibited an enhanced susceptibility to Botrytis cinerea, Alternaria brassicicola and Alternaria brassicae. Complementations of An-1 and T-DNA mutants using overexpression of a short transcript lacking the NB-ARC domain, or a genomic clone, restored resistance to all necrotrophic fungi. The elevated expression of RLM3(Col) on B. cinerea-susceptible mutants further suggested convergence in signalling and gene regulation between defence against B. cinerea and L. maculans. In the case of L. maculans, RLM3(Col) is required for efficient callose deposition downstream of RLM1(Col).

A new method for RAPD primers selection based on primer bias in nucleotide sequence data

Sequence analysis has proved that decamer nucleotides, used as primers of RAPD (random amplified polymorphic DNA), differ with each other greatly in number of annealing sites in the Arabidopsis thaliana genome. It is called the 'primer bias' by the authors. The biased primers produce a highly variable number of amplicons by polymerase chain reaction (PCR). The number of amplicons is proved to correlate with the number of annealing sites. Therefore, a statistical method is proposed for selecting efficient primers based on the primer bias in the genomic sequence. The method was tested by experiment in A. thaliana genome, and the results demonstrate that the method outperforms routine methods and can substantially increase the efficiency of RAPD methodologies. We also proved that the expressed sequence tags (ESTs) show a highly coincident bias pattern with that of the whole genomic sequence, and can therefore be used to assess efficiencies of primers for species whose genomic sequence data are currently unknown.