Characterization of disease resistance genes in the Brassica napus pangenome reveals significant structural variation

Methods based on single nucleotide polymorphism (SNP), copy number variation (CNV) and presence/absence variation (PAV) discovery provide a valuable resource to study gene structure and evolution. However, as a result of these structural variations, a single reference genome is unable to cover the entire gene content of a species. Therefore, pangenomics analysis is needed to ensure that the genomic diversity within a species is fully represented. Brassica napus is one of the most important oilseed crops in the world and exhibits variability in its resistance genes across different cultivars. Here, we characterized resistance gene distribution across 50 B. napus lines. We identified a total of 1749 resistance gene analogs (RGAs), of which 996 are core and 753 are variable, 368 of which are not present in the reference genome (cv. Darmor‐bzh). In addition, a total of 15 318 SNPs were predicted within 1030 of the RGAs. The results showed that core R‐genes harbour more SNPs than variable genes. More nucleotide binding site‐leucine‐rich repeat (NBS‐LRR) genes were located in clusters than as singletons, with variable genes more likely to be found in clusters. We identified 106 RGA candidates linked to blackleg resistance quantitative trait locus (QTL). This study provides a better understanding of resistance genes to target for genomics‐based improvement and improved disease resistance.     

Variation in abundance of predicted resistance genes in the Brassica oleracea pangenome

Brassica oleracea is an important agricultural species encompassing many vegetable crops including cabbage, cauliflower, broccoli and kale; however, it can be susceptible to a variety of fungal diseases such as clubroot, blackleg, leaf spot and downy mildew. Resistance to these diseases is meditated by specific disease resistance genes analogs (RGAs) which are differently distributed across B. oleracea lines. The sequenced reference cultivar does not contain all B. oleracea genes due to gene presence/absence variation between individuals, which makes it necessary to search for RGA candidates in the B. oleracea pangenome. Here we present a comparative analysis of RGA candidates in the pangenome of B. oleracea. We show that the presence of RGA candidates differs between lines and suggests that in B. oleracea, SNPs and presence/absence variation drive RGA diversity using separate mechanisms. We identified 59 RGA candidates linked to Sclerotinia, clubroot, and Fusarium wilt resistance QTL, and these findings have implications for crop breeding in B. oleracea, which may also be applicable in other crops species.     

PCR-based approach for mining of resistant gene analogues in taro (Colocasia esculenta)

Primers based on the conserved motifs were used to isolate nucleotide-binding sites (NBS) type sequences in taro (Colocasia esculenta). Cloning and sequencing identified three taro NBS-type sequences called resistance gene analogues (RGAs) that depicted similarity to other cloned RGA sequences. The deduced amino acid sequences of the RGAs detected the presence of conserved domains, viz. P-loop, categorising them with the NBS–leucine-rich repeat class gene family. Phylogenetic characterisation of the taro RGAs along with RGAs of other plant species grouped them with the non-toll interleukin receptor subclasses of the NBS sequences. The isolation and characterisation of taro RGAs have been reported for the first time in this study. This will provide a starting point towards characterisation of candidate resistance genes in taro and can act as a reference guide for future studies.     

Selection of stable Brassica napus-Brassica juncea recombinant lines resistant to blackleg (Leptosphaeria maculans). 2. A ‘to and fro’ strategy to localise and characterise interspecific introgressions on the B. napus genome

 A new strategy to localise and characterise interspecific introgressions in the genus Brassica is presented. It consists of the localisation of RAPD specific markers from the donor species (B. juncea) by RFLP on a genetic map of the recipient (B. napus) and on the observation of the disappearance of rapeseed markers in recombinant lines. With this method, we localised an interspecific introgression of B. juncea, which confers blackleg resistance at the cotyledon stage in B. napus, on the linkage group DY17 of the previously determined B. napus genetic map. The estimated size of the substituted B. napus fragment was 39 cM, and the resistance gene was introgressed into the rapeseed genome by homologous recombination. The significance of the different strategies used and the implication of these results in breeding programs are discussed. Received: 23 August 1997 / Accepted: 13 October 1997     

A new avirulence gene of Leptosphaeria maculans,AvrLm14, identifies a resistance source in American broccoli (Brassica oleracea) genotypes

In many cultivated crops, sources of resistance to diseases are sparse and rely on introgression from wild relatives. Agricultural crops often are allopolyploids resulting from interspecific crosses between related species, which are sources of diversity for resistance genes. This is the case for Brassica napus (oilseed rape, canola), an interspecific hybrid between Brassica rapa (turnip) and Brassica oleracea (cabbage). B. napus has a narrow genetic basis and few effective resistance genes against stem canker (blackleg) disease, caused by the fungus Leptosphaeria maculans, are currently available. B. rapa diversity has proven to be a valuable source of resistance (Rlm, LepR) genes, while B. oleracea genotypes were mostly considered susceptible. Here we identified a new resistance source in B. oleracea genotypes from America, potentially effective against French L. maculans isolates under both controlled and field conditions. Genetic analysis of fungal avirulence and subsequent cloning and validation identified a new avirulence gene termed AvrLm14 and suggested a typical gene-for-gene interaction between AvrLm14 and the postulated Rlm14 gene. AvrLm14 shares all the usual characteristics of L. maculans avirulence genes: it is hosted in a genomic region enriched in transposable elements and heterochromatin marks H3K9me3, its expression is repressed during vegetative growth but shows a strong overexpression 5–9 days following cotyledon infection, and it encodes a small secreted protein enriched in cysteine residues with few matches in databases. Similar to the previously cloned AvrLm10-A, AvrLm14 contributes to reduce lesion size on susceptible cotyledons, pointing to a complex interplay between effectors promoting or reducing lesion development.     

Characterisation and genetic mapping of resistance and defence gene analogs in cocoa (Theobroma cacao L.)

Disease resistance and defence gene analog (RGA/DGA) sequences were isolated in cocoa using a PCR approach with degenerate primers designed from conserved domains of plant resistance and defence genes: the NBS (nucleotide binding site) motif present in a number of resistance genes such as the tobacco N, sub-domains of plant serine/threonine kinases such as the Pto tomato gene, and conserved domains of two defence gene families: pathogenesis-related proteins (PR) of classes 2 and 5. Nucleotide identity between thirty six sequences isolated from cocoa and known resistance or defence genes varied from 58 to 80%. Amino acid sequences translated from corresponding coding sequences produced sequences without stop codons, except for one NBS –like sequence. Most of the RGAs could be mapped on the cocoa genome and three clusters of genes could be observed : NBS-like sequences clustered in two regions located on chromosomes 7 and 10, Pto-like sequences mapped in five genome regions of which one, located on chromosome 4, corresponded to a cluster of five different sequences. PR2-like sequences mapped in two regions located on chromosome 5 and 9 respectively. An enrichment of the genetic map with microsatellite markers allowed us to identify several co-localisations of RGAs, DGAs and QTL for resistance to Phytophthora detected in several progenies, particularly on chromosome 4 where a cluster of Pto-like sequences and 4 QTL for resistance to Phytophthora were observed. Many other serious diseases affect cocoa and the candidate genes, isolated in this study, could be of broader interest in cocoa disease management.     

Sequence‐level comparative analysis of the Brassica napus genome around two stearoyl‐ACP desaturase loci

We conducted a sequence‐level comparative analyses, at the scale of complete bacterial artificial chromosome (BAC) clones, between the genome of the most economically important Brassica species, Brassica napus (oilseed rape), and those of Brassica rapa, the genome of which is currently being sequenced, and Arabidopsis thaliana. We constructed a new B. napus BAC library and identified and sequenced clones that contain homoeologous regions of the genome including stearoyl‐ACP desaturase‐encoding genes. We sequenced the orthologous region of the genome of B. rapa and conducted comparative analyses between the Brassica sequences and those of the orthologous region of the genome of A. thaliana. The proportion of genes conserved (～56%) is lower than has been reported previously between A. thaliana and Brassica (～66%). The gene models for sets of conserved genes were used to determine the extent of nucleotide conservation of coding regions. This was found to be 84.2 ± 3.9% and 85.8 ± 3.7% between the B. napus A and C genomes, respectively, and that of A. thaliana, which is consistent with previous results for other Brassica species, and 97.5 ± 3.1% between the B. napus A genome and B. rapa, and 93.1 ± 4.9% between the B. napus C genome and B. rapa. The divergence of the B. napus genes from the A genome and the B. rapa genes was greater than anticipated and indicates that the A genome ancestor of the B. napus cultivar studied was relatively distantly related to the cultivar of B. rapa selected for genome sequencing.     

A genome-wide comparison of NB-LRR type of resistance gene analogs (RGA) in the plant kingdom

Plants express resistance (R) genes to recognize invaders and prevent the spread of pathogens. To analyze nucleotide binding site, leucine-rich repeat (NB-LRR) genes, we constructed a fast pipeline to predict and classify the R gene analogs (RGAs) by applying in-house matrices. With predicted ∼37,000 RGAs, we can directly compare RGA contents across entire plant lineages, from green algae to flowering plants. We focused on the highly divergent NBLRRs in land plants following the emergence of mosses. We identified entire loss of Toll/Interleukin-1 receptor, NBLRR (TNL) in Poaceae family of monocots and interestingly from Mimulus guttatus (a dicot), which leads to the possibility of species-specific TNL loss in other sequenced flowering plants. Using RGA maps, we have elucidated a positive correlation between the cluster sizes of NB-LRRs and their numbers. The cluster members were observed to consist of the same class of NB-LRRs or their variants, which were probably generated from a single locus for an R gene. Our website (), called plant resistance gene analog (PRGA), provides useful information, such as RGA annotations, tools for predicting RGAs, and analyzing domain profiles. Therefore, PRGA provides new insights into R-gene evolution and is useful in applying RGA as markers in breeding and or systematic studies.     

Isolation,genetic variation and expression of TIR-NBS-LRR resistance gene analogs from western white pine (<Emphasis Type="Italic">Pinus monticola </Emphasis>Dougl. ex. D. Don.)

Western white pine (Pinus monticola Dougl. ex. D. Don., WWP) shows genetic variation in disease resistance to white pine blister rust (Cronartium ribicola). Most plant disease resistance (R) genes encode proteins that belong to a superfamily with nucleotide-binding site domains (NBS) and C-terminal leucine-rich repeats (LRR). In this work a PCR strategy was used to clone R gene analogs (RGAs) from WWP using oligonucleotide primers based on the conserved sequence motifs in the NBS domain of angiosperm NBS-LRR genes. Sixty-seven NBS sequences were cloned from disease-resistant trees. BLAST searches in GenBank revealed that they shared significant identity to well-characterized R genes from angiosperms, including L and M genes from flax, the tobacco N gene and the soybean gene LM6. Sequence alignments revealed that the RGAs from WWP contained the conserved motifs identified in angiosperm NBS domains, especially those motifs specific for TIR-NBS-LRR proteins. Phylogenic analysis of plant R genes and RGAs indicated that all cloned WWP RGAs can be grouped into one major branch together with well-known R proteins carrying a TIR domain, suggesting they belong to the subfamily of TIR-NBS-LRR genes. In one phylogenic tree, WWP RGAs were further subdivided into fourteen clusters with an amino acid sequence identity threshold of 75%. cDNA cloning and RT-PCR analysis with gene-specific primers demonstrated that members of 10 of the 14 RGA classes were expressed in foliage tissues, suggesting that a large and diverse NBS-LRR gene family may be functional in conifers. These results provide evidence for the hypothesis that conifer RGAs share a common origin with R genes from angiosperms, and some of them may play important roles in defense mechanisms that confer disease resistance in western white pine. Ratios of non-synonymous to synonymous nucleotide substitutions (Ka/Ks) in the WWP NBS domains were greater than 1 or close to 1, indicating that diversifying selection and/or neutral selection operate on the NBS domains of the WWP RGA family.Communicated by R. Hagemann     

Mapping of genome-wide resistance gene analogs (RGAs) in maize (<Emphasis Type="Italic">Zea mays</Emphasis> L.)

Isolation and mapping of genome-wide resistance (R) gene analogs (RGAs) is of importance in identifying candidate(s) for a particular resistance gene/QTL. Here we reported our result in mapping totally 228 genome-wide RGAs in maize. By developing RGA-tagged markers and subsequent genotyping a population consisting of 294 recombinant inbred lines (RILs), 67 RGAs were genetically mapped on maize genome. Meanwhile, in silico mapping was conducted to anchor 113 RGAs by comparing all 228 RGAs to those anchored EST and BAC/BAC-end sequences via tblastx search (E-value < 10⁻²⁰). All RGAs from different mapping efforts were integrated into the existing SSR linkage map. After accounting for redundancy, the resultant RGA linkage map was composed of 153 RGAs that were mapped onto 172 loci on maize genome, and the mapped RGAs accounted for approximate three quarters of the genome-wide RGAs in maize. The extensive co-localizations were observed between mapped RGAs and resistance gene/QTL loci, implying the usefulness of this RGA linkage map in R gene cloning via candidate gene approach. Electronic supplementary material The online version of this article (doi:) contains supplementary material, which is available to authorized users. Wenkai Xiao, Jing Zhao and Shengci Fan have contributed equally to this research.     

The Brassica napus receptor‐like protein RLM2 is encoded by a second allele of the LepR3/Rlm2 blackleg resistance locus

Leucine‐rich repeat receptor‐like proteins (LRR‐RLPs) are highly adaptable parts of the signalling apparatus for extracellular detection of plant pathogens. Resistance to blackleg disease of Brassica spp. caused by Leptosphaeria maculans is largely governed by host race‐specific R‐genes, including the LRR‐RLP gene LepR3. The blackleg resistance gene Rlm2 was previously mapped to the same genetic interval as LepR3. In this study, the LepR3 locus of the Rlm2 Brassica napus line ‘Glacier DH24287’ was cloned, and B. napus transformants were analysed for recovery of the Rlm2 phenotype. Multiple B. napus, B. rapa and B. juncea lines were assessed for sequence variation at the locus. Rlm2 was found to be an allelic variant of the LepR3 LRR‐RLP locus, conveying race‐specific resistance to L. maculans isolates harbouring AvrLm2. Several defence‐related LRR‐RLPs have previously been shown to associate with the RLK SOBIR1 to facilitate defence signalling. Bimolecular fluorescence complementation (BiFC) and co‐immunoprecipitation of RLM2‐SOBIR1 studies revealed that RLM2 interacts with SOBIR1 of Arabidopsis thaliana when co‐expressed in Nicotiana benthamiana. The interaction of RLM2 with AtSOBIR1 is suggestive of a conserved defence signalling pathway between B. napus and its close relative A. thaliana.     

Cloning and characterisation of a family of disease resistance gene analogs from wheat and barley

 The most common class of plant disease resistance (R) genes cloned so far belong to the NBS-LRR group which contain nucleotide-binding sites (NBS) and a leucine-rich repeat (LRR). Specific primer sequences derived from a previously isolated NBS-LRR sequence at the Cre3 locus, which confers resistance to cereal cyst nematode (CCN) in wheat (Triticum aestivum L.) were used in isolating a family of resistance gene analogs (RGA) through a polymerase chain reaction (PCR) cloning approach. The cloning, analysis and genetic mapping of a family of RGAs from wheat (cv ‘Chinese Spring’) and barley (Hordeum vulgare L. cvs ‘Chebec’ and ‘Harrington’) are presented. The wheat and barley RGAs contain other conserved motifs present in known R genes from other plants and share between 55–99% amino acid sequence identity to the NBS-LRR sequence at the Cre3 locus. Phylogenetic analysis of the RGAs with other cloned R genes and RGAs from various plant species indicate that they belong to a superfamily of NBS-containing genes. Two of the barley derived RGAs were mapped onto loci on chromosomes 2H (2), 5H (7) and 7H (1) using barley doubled haploid (DH) mapping populations. Some of these loci identified are associated with regions carrying resistance to CCN and corn leaf aphid. Received: 6 January 1998 / Accepted: 1 April 1998     

Cloning of fatty acid elongase1 gene and molecular identification of A and C genome in <Emphasis Type="Italic">Brassica</Emphasis> species

The fatty acid elongase 1 (FAE1) genes of Brassic napus were cloned from two cultivars, i.e. Zhongshuan No. 9 with low erucic acid content, and Zhongyou 821 with high erucic acid content, using the degenerate PCR primers. The sequence analysis showed that there was no intron within the FAE1 genes. The FAE1 genes from Zhongyou 821 contained a coding sequence of 1521 nucleotides, and those cloned from Zhongshuan No. 9 contained a 1517 bp coding sequence. Alignment of the FAE1 sequences from Brassica rapa, B. oleracea and B. napus detected 31 single nucleotide polymorphic sites (2.03%), which resulted in 7 amino-acid substitutions. Further analysis indicated that 19 SNPs were genome-specific, of which, 95% were synonymous mutations. The nucleotide substitution at position 1217 in the FAE1 genes led to a specific site of restricted cleavage. An AvrII cleavage site was present only in the C genome genes and absent in the A genome FAE1 genes. Digestion profile of the FAE1 sequences from B. rapa, B. oleracea and B. napus produced with AvrII confirmed that the FAE1 genes of B. oleracea origin was recognized and digested, while that of B. rapa origin could not. The results indicated that by AvrII cleavage it was possible to distinguish B. rapa from B. oleracea and between the A and C genome of B. napus. In addition, the FAE1 genes could be used as marker genes to detect the pollen flow of B. napus, thus providing an alternative method for risk assessment of gene flow.     

Cloning and in silico Mapping of Resistance Gene Analogues Isolated from Rice Lines Containing Known Genes for Blast Resistance

We amplified resistance gene analogues (RGAs) from the genomic DNA of 10 rice lines having varying degree of resistance to Magnaporthe grisea by using degenerate primers and various RGAs were mapped in silico on different rice chromosomes. The amplified products were grouped into 3–8 restriction fragment length polymorphic classes by using Mbo1 and Alu1 restriction enzymes. Of 98 RGAs obtained in this study, 65 RGA clones showed more than 95% homology with various RGAs sequences present in the GenBank. Phylogenetic analysis of these RGAs formed 11 groups. Using sequence homology approach, RGAs isolated in this study were physically mapped on 23 loci on chromosomes 1, 2, 3, 4, 5, 6, 7, 8, 10, 11 and 12. Twenty RGAs were mapped near to the chromosomal regions containing known genes/QTLs for rice blast, bacterial leaf blight and sheath blight resistance. Thirty‐nine RGA sequences also contained open reading frame representing signature of potential disease resistance genes.     

Selection of stable Brassica napus-B. juncea recombinant lines resistant to blackleg (Leptosphaeria maculans). 1. Identification of molecular markers, chromosomal and genomic origin of the introgression

 A scheme of selection combining selfing and backcross was applied to a B. napus line with the blackleg resistance from B. juncea in order to transfer this resistance to a winter oilseed rape variety. Cytogenetic analyses combined with cotyledon blackleg resistance tests at each generation allowed us to obtain a recombinant line showing regular meiotic behavior. The resistance is monogenic and is highly efficient under field conditions. Four-hundred RAPD primers were tested on two segregating populations by bulk segregant analysis. Three markers totally linked to the introgression were identified. The analysis of these markers on both sets of B. napus-B. nigra and B. oleracea-B. nigra addition lines revealed that they are not located on the B4 chromosome of B. nigra, which has already been shown to carry a blackleg resistance gene, but rather on the B8 chromosome. We confirmed that the resistance gene is carried by the B genome of B. juncea. Based on these data, two hypotheses, one involving chromosome rearrangements between the two B genomes of B. nigra and B. juncea, and the other based on a more probable digenic control of the resistance within B. juncea, are discussed. Received: 30 May 1997 / Accepted: 23 June 1997     

Targeted isolation,sequence analysis,and physical mapping of nonTIR NBS-LRR genes in soybean

Most cloned plant disease resistance genes (R-genes) code for proteins belonging to the nucleotide binding site (NBS) leucine-rich repeat (LRR) superfamily. NBS-LRRs can be divided into two classes based on the presence of a TIR domain (Toll and interleukin receptor-like sequence) or a coiled coil motif (nonTIR) in their N-terminus. We used conserved motifs specific to nonTIR-NBS-LRR sequences in a targeted PCR approach to generate nearly 50 genomic soybean sequences with strong homology to known resistance gene analogs (RGAs) of the nonTIR class. Phylogenetic analysis classified these sequences into four main subclasses. A representative clone from each subclass was used for genetic mapping, bacterial artificial chromosome (BAC) library screening, and construction of RGA-containing BAC contigs. Of the 14 RGAs that could be mapped genetically, 12 localized to a 25-cM region of soybean linkage group F already known to contain several classical disease resistance loci. A majority of the genomic region encompassing the RGAs was physically isolated in eight BAC contigs, together spanning more than 1 Mb of genomic sequence with at least 12 RGA copies. Phylogenetic and sequence analysis, together with genetic and physical mapping, provided insights into the genome organization and evolution of this large cluster of soybean RGAs. Received: 8 May 2001 / Accepted: 30 June 2001     

Resistance gene analogues from rice: cloning, sequencing and mapping

 Degenerate oligonucleotide primers were designed on the basis of nucleotide-binding-site (NBS) motifs conserved between resistance genes of Arabidopsis, flax and tobacco and subsequently used as PCR primers to amplify resistance gene analogues (RGA) in rice. Primers amplified a major band of approximately 500 bp. Restriction analysis of the amplified product revealed that the band was made up of several different fragments. Many of these fragments were cloned. Sixty different cloned fragments were analysed and assigned to 14 categories based on Southern blot analysis. Fourteen clones, each representing one of the 14 categories of RGAs were mapped onto the rice genetic map using a Nipponbare ( japonica)בKasalath’ (indica) mapping population consisting of 186 F₂ lines. Of the 14 clones representing each class 12 could be mapped onto five different chromosomes of rice with a major cluster of 8 RGAs on chromosome 11. Our results indicate that it is possible to use sequence homology from conserved motifs of known resistance genes to amplify candidate resistance genes from diverse plant taxa. Received: 23 September 1998 / Accepted: 28 November 1998     

Cloning of fatty acid elongase1 gene and molecular identification of A and C genome in Brassica species

The fatty acid elongase 1 (FAE1) genes of Brassic napus were cloned from two cultivars, i.e. Zhongshuan No. 9 with low erucic acid content, and Zhongyou 821 with high erucic acid content, using the degenerate PCR primers. The sequence analysis showed that there was no intron within the FAE1 genes. The FAE1 genes from Zhongyou 821 contained a coding sequence of 1521 nucleotides, and those cloned from Zhongshuan No. 9 contained a 1517 bp coding sequence. Alignment of the FAE1 sequences from Brassica rapa, B. oleracea and B. napus detected 31 single nucleotide polymorphic sites (2.03%), which resulted in 7 amino-acid substitutions. Further analysis indicated that 19 SNPs were genome-specific, of which, 95% were synonymous mutations. The nucleotide substitution at position 1217 in the FAE1 genes led to a specific site of restricted cleavage. An AvrII cleavage site was present only in the C genome genes and absent in the A genome FAE1 genes. Digestion profile of the FAE1 sequences from B. rapa, B. oleracea and B. napus produced with AvrII confirmed that the FAE1 genes of B. oleracea origin was recognized and digested, while that of B. rapa origin could not. The results indicated that by AvrII cleavage it was possible to distinguish B. rapa from B. oleracea and between the A and C genome of B. napus. In addition, the FAE1 genes could be used as marker genes to detect the pollen flow of B. napus, thus providing an alternative method for risk assessment of gene flow. Supported by the National Natural Science Foundation of China (Grant No. 30471099), Development Plan of the State Key Fundamental Research of China (Grant No. 2006CB101600), and the National High Technology and Development Program of China (Grant No. 2006AA10A113)