ORIGINAL RESEARCH article
NBS-Encoding Genes in Brassica napus Evolved Rapidly After Allopolyploidization and Co-localize With Known Disease Resistance Loci
- 1Institute of Crop and Nuclear Technology Utilization, Zhejiang Academy of Agricultural Sciences, Hangzhou, China
- 2Department of Plant Breeding, IFZ Research Centre for Biosystems, Land Use and Nutrition, Justus Liebig University Giessen, Giessen, Germany
- 3Key Laboratory of Crop Physiology, Ecology and Genetic Breeding, Ministry of Education, Agronomy College, Jiangxi Agricultural University, Nanchang, China
Genes containing nucleotide-binding sites (NBS) play an important role in pathogen resistance in plants. However, the evolutionary fate of NBS-encoding genes after formation of allotetraploid Brassica napus (AnAnCnCn, 2n = 38) is still unknown. We performed a genome-wide comparison of putatively functional NBS-encoding genes in B. napus and its progenitor species Brassica rapa (ArAr, 2n = 20) and Brassica oleracea (CoCo, 2n = 18), identifying 464, 202, and 146 putatively functional NBS-encoding genes respectively, with genes unevenly distributed in several clusters. The An-subgenome of B. napus possessed similar numbers of NBS-encoding genes (191 genes) to the Ar genome of B. rapa (202 genes) and similar clustering patterns. However, the Cn genome of B. napus had many more genes (273) than the B. oleracea Co genome (146), with different clustering trends. Only 97 NBS-encoding genes (66.4%) in B. oleracea were homologous with NBS-encoding genes in B. napus, while 176 NBS-encoding genes (87.1%) were homologous between B. rapa and B. napus. These results suggest a greater diversification of NBS-encoding genes in the C genome may have occurred after formation of B. napus. Although most NBS-encoding genes in B. napus appeared to derive from the progenitors, the birth and death of several NBS-encoding genes was also putatively mediated by non-homologous recombination. The Ka/Ks values of most homologous pairs between B. napus and the progenitor species were less than 1, suggesting purifying selection during B. napus evolution. The majority of NBS-encoding genes (60% in all species) showed higher expression levels in root tissue (out of root, leaf, stem, seed and flower tissue types). Comparative analysis of NBS-encoding genes with mapped resistance QTL against three major diseases of B. napus (blackleg, clubroot and Sclerotinia stem rot) found 204 NBS-encoding genes in B. napus located within 71 resistance QTL intervals. The majority of NBS-encoding genes were co-located with resistance QTLs against a single disease, while 47 genes were co-located with QTLs against two diseases and 3 genes were co-located with QTLs against all three. Our results revealed significant variation as well as interesting evolutionary trajectories of NBS-encoding genes in the different Brassica subgenomes, while co-localization of NBS-encoding genes and resistance QTL may facilitate resistance breeding in oilseed rape.
Plants are surrounded by a vast array of invaders, such as viruses, bacteria, fungi, nematodes and pests, many of which can cause diseases (Dangl and Jones, 2001; Meyers et al., 2003; McHale et al., 2006). In order to cope with invader attacks, plants have evolved sophisticated immune mechanisms to protect themselves against their natural enemies. The best known are the vast numbers of resistance genes in plants, which play a central role in recognizing effectors from pathogens and in triggering downstream signaling during plant responses to pathogen invasions (Liu et al., 2007; Yang et al., 2013).
Genes containing a nucleotide-binding site (NBS), namely NBS-encoding genes, constitute one of the largest plant resistance gene families (∼80%) (Meyers et al., 2003; Liu et al., 2007; Yang et al., 2013). The NBS domain was found to bind and hydrolyze ATP or GTP, and primarily functions as a signal transduction switch following pathogen recognition (Dangl and Jones, 2001; Martin et al., 2003). NBS-encoding genes typically comprise three principal domains, with the NBS domain in the middle region, flanked by a leucine-rich repeat (LRR) domain at the C-terminus, and by a toll/interleukin-1 receptor (TIR) or coiled-coil (CC) at the N-terminus (Cannon et al., 2002; Meyers et al., 2003; Shao et al., 2014; Zhang et al., 2016). The central NBS domain encodes several motifs consisting of 10–30 amino acids (aa), and is typically highly conserved (Meyers et al., 1999; Yue et al., 2012), whereas the C-terminal LRR domain exhibits high diversity and has been associated with pathogen recognition (Kobe and Deisenhofer, 1995; Leister and Katagiri, 2000; Dangl and Jones, 2001). According to the presence or absence of the N-terminal TIR domain, NBS-encoding genes are further classified into TIR-NBS-LRR (TNL) or TIR-NBS (TN) genes and non-TIR-NBS-LRR (non-TNL) or (non-TN) genes (Meyers et al., 1999). Based on the presence of the CC or other domains at the N-terminus, non-TNL and non-TN genes can be further divided into CC-NBS-LRR (CNL) or CC-NBS (CN) genes and X-NBS-LRR (XNL) or X-NBS (XN) genes (Dangl et al., 2001).
With the availability of genomic data for an increasing number of species, a set of NBS-encoding genes has been identified at the genome level in more than 30 angiosperms (Bai et al., 2002; Meyers et al., 2003; Monosi et al., 2004; Zhou et al., 2004; Yang et al., 2006, 2008a,b; Ameline-Torregrosa et al., 2008; Mun et al., 2009; Porter et al., 2009; Chen et al., 2010; Li J. et al., 2010; Li X. et al., 2010; Guo et al., 2011; Jupe et al., 2012; Lozano et al., 2012; Luo et al., 2012; Tan and Wu, 2012; Jupe et al., 2013; Lin et al., 2013; Wan et al., 2013; Andolfo et al., 2014; Arya et al., 2014; Kim et al., 2014; Shao et al., 2014; Wu et al., 2014;Yu et al., 2014; Wei C. et al., 2016; Qian et al., 2017; Xiang et al., 2017). Comparative studies of the evolutionary history of NBS-encoding genes were further performed in a number of clades over recent years to illuminate the evolutionary characteristics of NBS-encoding genes. For instance, frequent gene losses and a limited number of gene duplications were pointed out in studies targeting Cucumis sativus, C. melo, and Citrullus lanatus of the Cucurbitaceae family (Lin et al., 2013). A comparative genomic analysis of four Poaceae species revealed that the number of NBS-encoding genes in Zea mays was only half that in Sorghum bicolor and Brachypodium distachyon and a fourth of that in Oryza sativa (Li et al., 2016). Gene loss and retention patterns and a pattern of gene expansion followed by contraction were identified in various species in the Brassicaceae family (Peele et al., 2014; Zhang et al., 2016). In the Fabaceae and Rosaceae, gene expansions were also the most frequently observed evolutionary pattern (Shao et al., 2014; Jia et al., 2015; Zhang et al., 2016).
Oilseed rape (Brassica napus, AnAnCnCn, 2n = 38) is an allopolyploid that originated from spontaneous hybridization events between the two diploid Brassica species B. rapa (ArAr, 2n = 20) and B. oleracea (CoCo, 2n = 18) in the last 10,000 years (Nagaharu, 1935; Chalhoub et al., 2014). Despite the short domestication history of rapeseed (∼300–400 years; Gómez-Campo and Prakash, 1999), it is today one of the most important oil crops worldwide. Hybridization and polyploidisation often produce species with superior resistances or environmental tolerances relative to their progenitor species (Leitch and Leitch, 2008). The exact mechanisms responsible for this effect are unknown, but in rapeseed could result from hybridization, polyploidisation or domestication processes shaping genome evolution. The evolutionary impact of these processes on NBS-encoding genes, which are major players in plant disease resistance, is hence of particular interest. With the availability of genomic data for B. napus (Chalhoub et al., 2014) and its two progenitor species B. rapa (Wang et al., 2011) and B. oleracea (Liu et al., 2014; Parkin et al., 2014), NBS-encoding genes can be systematically investigated to elucidate their role in contributing to differences in disease resistance between the three species, and help to decipher the mechanisms underlying disease resistance in B. napus. Previous studies have identified and compared numbers, types and locations of NBS-encoding genes in B. rapa, B. oleracea, and B. napus. Chalhoub et al. (2014) found 425 NBS-LRRs in B. napus (245 in the C genome and 180 in the A genome), similar numbers were reported for the B. rapa and B. oleracea genomes but with only 75% conservation of synteny. Alamery et al. (2017) identified 641 NBS-LRRs in B. napus in total, of which only 365 were intact and hence putatively functional.Yu et al. (2014) identified 239 NBS-LRRs in B. oleracea, which was updated to 556 NBS-LRRs in the B. oleracea pan-genome (Bayer et al., 2018).
In the present work, a genome-wide characterization of NBS-encoding genes was performed in B. rapa, B. oleracea, and B. napus. Multiple approaches were utilized to assess the genome architecture and evolutionary characteristics of NBS genes, including genomic distribution, homologous genes, sequence similarities, selection signals, phylogenetic relationships and expression patterns. Resistance genes were also co-localized with previously identified quantitative trait loci (QTL) conferring resistance against major oilseed rape pathogens blackleg (Leptosphaeria maculans), clubroot (Plasmodiophora brassicae) and Sclerotinia stem rot (Sclerotinia sclerotiorum), highlighting potential disease resistance candidate genes for future work. This analysis provided genome level insights into the evolution of disease resistance genes in B. napus, shedding light on possible mechanisms of disease resistance for future disease resistance breeding in rapeseed.
Materials and Methods
Identification of NBS-Encoding Genes in B. napus, B. rapa, and B. oleracea
The entire genome sequences and annotation data for B. rapa (Wang et al., 2011), B. oleracea (Liu et al., 2014; Parkin et al., 2014) and B. napus (Chalhoub et al., 2014) were downloaded from the BRAD database1. NBS-encoding genes in the three species were identified by using the amino acid sequence of the Pfam NB-ARC domain (PF00931) as a “blastp” query against all known protein sequences, via the HMMER V3.0 program with “trusted cutoff” as the threshold (Finn et al., 2011). The obtained hits were further submitted to the Pfam website2 to verify the presence of the NB-ARC domain. Those proteins with verified NB-ARC domains were further subdivided into groups based on the structure of the N-terminal and C-terminal domains of the protein. By using the Pfam database3, SMART protein motif analyses4, and the COILS program5, the presence or absence of TIR, CC, LRR domains of NBS-encoding genes was identified, which was used to classify the NBS-encoding genes into different groups.
Location Anchoring and Gene Cluster Analysis of NBS-Encoding Genes in B. rapa, B. oleracea, and B. napus
The chromosome size and physical position of NBS-encoding genes were downloaded from the BRAD database6 for location anchoring of NBS-encoding genes. The visualization of NBS-encoding genes on 10, 9, and 19 chromosomes (assembled pseudomolecules) of B. rapa, B. oleracea, and B. napus was drawn by “PhenoGram Plot”.7 According to the gene cluster definition proposed by Richly et al. (2002) and Meyers et al. (2003), two or more genes located within eight ORFs of each other were treated as a gene cluster.
Detection of Duplicated Genes Within Each Species and Identification of Homologous Gene Pairs Between B. rapa, B. oleracea, and B. napus
To detect the duplicated NBS-encoding genes within individual species and to identify homologous gene pairs between species, the BLASTP program of Blast2GO (Conesa et al., 2005) was employed using protein sequences with stringent parameters, i.e., E-value cutoff (e = 0.0) (Ariyarathna and Francki, 2016). The gene pairs across genomes were visualized using the Circos software (Krzywinski et al., 2009). If no homologous gene was found for a particular gene using this stringent criterion, the top hit was selected and designated as the most similar gene.
Multiple Alignments and Phylogenetic Analysis of NBS-Encoding Genes
Multiple alignments of amino acid sequences were performed using the program “Muscle” with default options (Thompson et al., 1994). The software RAxML was then used to construct phylogenetic trees based on the Maximum Likelihood (ML) method (Stamatakis, 2015). Phylogenetic trees were drawn by the software Figtree8.
Non-synonymous/Synonymous Substitution (Ka/Ks) Ratios of Homologous Gene Pairs
Protein sequence variation between each homologous pair was firstly calculated via sequence alignment using the program Muscle with default options (Thompson et al., 1994). Based on the alignment results, the number of nucleotide variants was then calculated by reverse translation of the protein sequence to nucleotides. The ratio of nucleotide variants was calculated by the formula: ratio of nucleotide variants = number of variant nucleotides/the number of aligned nucleotides (gaps were not considered). To investigate selective pressure on NBS-encoding genes, the ratio of non-synonymous substitutions to synonymous substitutions (Ka/Ks) was calculated. Firstly, the protein sequences of NBS-encoding genes in each homologous gene family were aligned using the program Muscle with default options (Thompson et al., 1994). Then, non-synonymous substitutions (Ka) and synonymous substitutions (Ks), and the ratio between them (Ka/Ks) were calculated in each homologous gene family using the tool yn00 in the paml package (Yang, 1997). Positive and negative selective pressure were justified as a Ka/Ks ratio > 1 and Ka/Ks ratio < 1, respectively, while a ratio of 1 indicated neutral evolution (Yang, 2007). Ka/Ks ratio cutoffs of > 1.2 and < 0.8 were chosen to identify positively and negatively selected genes, respectively, with slight modification of the method of Mayer et al. (2011).
RNA-seq Data Analysis of NBS-Encoding Genes
In order to analyze the expression of the NBS-encoding genes, RNA-seq data of NBS-encoding genes of B. rapa, B. oleracea and B. napus that was generated earlier and submitted to the GEO database was downloaded. RNA-seq data of five tissues (stem transcripts, root transcripts, leaf transcripts, flower transcripts and seed transcripts) of B. napus, B. rapa, and B. oleracea were obtained from Tong et al. (2013); Liu et al. (2014), and Miao et al. (2016), respectively. Fragments per kilobase of exon model per million mapped reads (FPKM) were calculated for NBS-encoding genes. The FPKM values were log2 transformed, and a hierarchical cluster was created using the Genesis software using the hierarchical clustering algorithm (Eisen et al., 1998).
Co-localization of NBS-Encoding Genes and Mapped Quantitative Resistance Loci Against Three Major Diseases of B. napus
Resistance quantitative trait loci for blackleg (Leptosphaeria maculans and L. biglobosa), clubroot (Plasmodiophora brassicae), and Sclerotinia stem rot (Sclerotinia sclerotiorum), were collected from the literature (Hirai et al., 2004; Suwabe et al., 2006; Zhao et al., 2006; Delourme et al., 2008; Sakamoto et al., 2008; Kaur et al., 2009; Nagaoka et al., 2010; Raman et al., 2012; Jian et al., 2013; Kato et al., 2013; Mei et al., 2013; Wei et al., 2014; Wei L. et al., 2016; Larkan et al., 2016; Li et al., 2016; Sanjaya et al., 2016; Wu et al., 2016; Kumar et al., 2018; Pang et al., 2018) (Supplementary Table S1). The QTL intervals were aligned to the reference genomes of B. napus9 by BLAST analysis of the sequences of SSR primers or probe sequences linked to the QTLs, using the method of Cai et al. (2012) with slight modifications. The alignment criteria were set to allow three mismatches and one gap for a given primer pair. When a marker had multiple amplification loci on the same chromosome, an accurate position for a particular locus was determined manually by referring to the physical positions of its upstream and downstream markers.
Identification of NBS-Encoding Genes in B. napus, B. rapa, and B. oleracea
A total of 464, 202, and 146 NBS-encoding genes were identified in B. napus, B. rapa, and B. oleracea, respectively (Table 1 and Figure 1). The genome of B. napus contained 116 additional NBS-encoding genes when compared to the total number of NBS-encoding genes present in the two progenitor species, a considerable increase (33.33% more than the total number of NBS-encoding genes in the two progenitor species). According to gene structure and protein motifs, we categorized these NBS-encoding genes into six different classes: TNL (134, 76, and 19 for B. napus, B. rapa, and B. oleracea, respectively), CNL (16, 3, and 19), XNL (18, 13, and 61), TN (37, 45, and 85), CN (18, 3, and 45) and XN (37, 63, and 120) (Table 1). Except for the CNL group, which showed the same gene number in B. napus compared to the total in the two progenitor species, all the other five classes of NBS-encoding genes contained more genes in B. napus than in the two progenitor species. The greatest increase in B. napus relative to the progenitors was for the TNL gene class (39 additional genes).
Figure 1. The physical location of NBS-encoding genes in Ar genome of B. rapa, Co genome of B. oleracea and An- and Cn- subgenomes of B. napus. The circles in different colors represented types of NBS-encoding genes.
Sequence Comparison of NBS-Encoding Genes Between B. napus and Its Progenitor Species
The size range of NBS-encoding genes in B. napus was much greater than NBS-encoding genes in the progenitors (83–244393 nt; as compared to 518–15805 nt and 164–13901 nt in B. rapa and B. oleracea, respectively, Table 1). B. napus also had more exons on average than the progenitor species: on average 5.31 exons compared to 4.52 and 4.25 in B. rapa and B. oleracea respectively, and a maximum of 40 exons compared to 22 and 32 in B. rapa and B. oleracea respectively.
In order to compare variation in NBS-encoding genes between B. napus and its progenitors, BLAST was used for the 464 NBS-encoding genes in B. napus against the NBS-encoding genes in the two progenitor species, and the best matching gene pairs were compared for total gene length and number of exons. On average, a 219 bp decrease in sequence length was observed for NBS-encoding genes in B. napus when compared with the corresponding best match for NBS-encoding genes in the progenitor species (a total of 10,1719 bp decrease). By contrast, the average number of exons in B. napus relative to its progenitor species increased by 0.91 (an additional 421 exons genome-wide). These results support a trend of overall decrease in sequence length but an increase in number of exons for NBS-encoding genes in B. napus after its formation from the progenitor species.
A total of 1456 fragments (>10 bp) from 412 NBS-encoding genes in B. napus did not align to the corresponding best matching NBS-encoding genes in the progenitor species. These fragments were subsequently aligned to the whole genome sequences of B. rapa and B. oleracea using BLAST. Of these, 1083 fragments (74%) from 372 NBS-encoding genes of B. napus could be completely or partially aligned with the genomic sequences of the two progenitors (e = 1.0). These results suggest abundant non-homologous recombination between the NBS-encoding genes and other genomic sequences from progenitor species during the processes of hybridization and polyploidization, or the subsequent domestication. The remaining 373 fragments not only did not align to the progenitor genomes, but also seldom hit with other species on NCBI (e = 10). However, 8 of the 373 fragments matched known TE motifs from A. thaliana, suggesting that some of these fragments may have formed by TE-mediated mutation processes (Wicker et al., 2016). Additionally, 376 fragments with no best hits and low homology (<35%) with the progenitor genomes were further extracted for alignments using BLAST (e = 10) against all species on NCBI. 114 sequences showed no hit with any species, but three fragments matched with known TE motifs from A. thaliana, suggesting these sequences may also have arisen via mutational processes and TE insertion within B. napus. Of the 262 sequences which could be anchored, 92 sequences matched best with the genomic sequences of the progenitors (after use of a less stringent e value), suggesting drastic variation within these sequences, and 84 sequences had a highest match with the genomes of other Brassicaceae species (46.3% from Raphanus sativus and Arabis alpine). A further 31 sequences had top matches outside the Brassicaceae family: 23 sequences were best matched with some microorganisms; 26 sequences were best matched with some worms, e.g., nematodes and trematodes; and 6 sequences were best matched with some other animals.
Genomic Distribution of NBS-Encoding Genes Along Chromosomes in the Three Species
The genomic distribution of NBS-encoding genes between recent allopolyploid B. napus (AnAnCnCn) and its progenitors B. rapa (ArAr) and B. oleracea (CoCo) was determined (Figure 2). Hundred and ninety one NBS-encoding genes were identified in the An-subgenome of B. napus, which was quite similar to the number of NBS-encoding genes (202) in the Ar genome of B. rapa. The distribution of these genes on different chromosomes in the A genome was uneven, although the An-subgenome of B. napus and the Ar genome of B. rapa showed similar genomic distributions. Chromosome A09 contained the most NBS-encoding genes in both A genomes, with 38 (19.9%) and 47 (23.27%) NBS-encoding genes located on An09 and Ar09, respectively. 273 NBS-encoding genes were identified in the Cn-subgenome of B. napus, which was a considerably greater number than were identified in the Co genome of B. oleracea (146). NBS-encoding genes were also unevenly distributed between the Cn-subgenome and Co genomes. The chromosome with the most NBS-encoding genes differed between the two C genomes: Cn9 and Co7 contained the highest numbers of NBS-encoding genes in each genome, with 52 (19.26%) and 25 (17.12%) NBS-encoding genes, respectively.
Figure 2. Comparison of NBS-encoding genes between the subgenomes of B. napus (AnAnCnCn) and the genomes of B. rapa (ArAr) and B. oleracea (CoCo). (A) The comparison of gene numbers between the genomes of Ar and An, and Co and Cn. The y-axis indicates the number of NBS-encoding genes located on each individual chromosome, and the x-axis represents the chromosomes. (B) The comparison of NBS-encoding gene clusters between the Ar and An, and Co and Cn genomes. The y-axis represents the chromosomes, while the x-axis represents the number of gene clusters and the proportion of genes located in each gene cluster, respectively.
Resistance genes residing in clusters may facilitate the evolutionary process of novel resistance gene production via increased tandem duplication and gene recombination. Using the definition of a gene “cluster” defined by Richly et al. (2002) and Meyers et al. (2003) as two or more genes falling within eight ORFs, we identified 44, 27, and 90 NBS-encoding gene clusters in B. rapa, B. oleracea, and B. napus, which included 117, 81, and 256 NBS-encoding genes respectively in the three species (57.92% in B. rapa, 55.48% in B. oleracea and 55.53% in B. napus). Both the number of gene clusters and the number of genes within clusters in B. napus were much more than the sum of the two progenitor species. Chromosomes An09 and Cn9 in B. napus possessed the most abundant gene clusters genome-wide. In the two progenitor species, chromosomes Ar09 and Co7 possessed the most abundant gene clusters in B. rapa and B. oleracea, respectively.
Comparative Homology of NBS-Encoding Genes Between Genomes
NBS-encoding genes were assessed for homology between the three species’ genomes (Table 2). Three hundred and forty six NBS-encoding genes of B. napus, forming 3454 homologous gene pairs, were identified between B. napus and its two progenitor species, with uneven distribution between chromosomes (Figure 3). Chromosome Co9 had the highest number of homologous gene pairs (432 pairs) and Co4 possessed the fewest homologous gene pairs (52 pairs). Of these homologous gene pairs, 70.9 and 50.2% of NBS-encoding genes in B. napus were homologous with those from B. rapa and B. oleracea, respectively (Table 2). For the 202 NBS-encoding genes in B. rapa, 176 (87.1%) were homologous with NBS-encoding genes in B. napus, while only 97 (66.4%) of NBS-encoding genes in B. oleracea were homologous with NBS-encoding genes in B. napus (Table 2). Although the remaining NBS-encoding genes of B. rapa and B. oleracea lacked clear homologous NBS-encoding gene partners in B. napus, 24 (11.9%) and 45 (30.8%) of genes still showed best hits to NBS-genes in B. napus, and only 2 (1.0%) and 4 (2.7%) of NBS-genes in B. rapa and B. oleracea had altered NBS-domains in B. napus (Table 2).
Table 2. Homology analysis of NBS-encoding genes between B. napus (2n = AnAnCnCn), B. rapa (2n = ArAr) and B. oleracea (2n = CoCo).
Figure 3. NBS-encoding gene homologs between the B. napus, B. Rapa, and B. oleracea genomes. Colored bars represent chromosomes of the three species. Green lines indicate the relationships of homologous gene pairs between genomes.
A total of 811 homologs of NBS-encoding genes were identified between the two progenitor species (Figure 3). Of these, 128 NBS-encoding genes (63.4%) in B. rapa were homologous to NBS-encoding genes in B. oleracea, and 79 NBS-encoding genes (54.1%) in B. oleracea were homologous to NBS-encoding genes in B. rapa. As expected, NBS-encoding gene sequences varied more between the two progenitor species than between the progenitors and B. napus, due to the longer evolutionary divergence time between the two progenitor species.
Overall, 1509 NBS-encoding homologous gene pairs were identified between the An- and Cn-subgenomes of B. napus (Figure 3), which was almost twice the number observed between the Ar and Co genomes of B. rapa and B. oleracea. This result indicated the genetic assimilation of NBS-encoding genes might occur between An and Cn sub-genomes during the formation and domestication of B. napus. The proportions of genes with homologs in the An subgenome (67.5%, 129 genes) and Cn subgenome (54.8%; 149 genes) of B. napus were similar to the proportions of genes with homologs between the two progenitor species.
Analysis of Selective Pressure on NBS-Encoding Genes
Using NBS-encoding gene homologs identified between genomes, we calculated Ks (synonymous substitution rate) and Ka (non-synonymous substitution rate) values for homologous gene pairs. In total, 3454 homologous gene pairs were identified between B. napus and two progenitor species, including 2774 homologous gene pairs between B. napus and B. rapa and 680 homologous gene pairs between B. napus and B. oleracea. Overall, 98.5 and 97.34% of the homologous gene pairs between B. napus and the two progenitor species showed Ka/Ks ratios less than 1 (99.4 and 94.9% for B. rapa and B. oleracea, respectively) and 0.8 (98.27 and 93.53% for B. rapa and B. oleracea, respectively) (Table 3), and the majority of homologous genes showed Ka/Ks ratios ranging from 0.3 to 0.4. These results indicate the presence of purifying selection (selection against change) on most of the NBS-encoding genes following the divergence of B. napus from the progenitor species. The mean Ka/Ks value for all homologous gene pairs between B. napus and B. rapa was 0.38, slightly but not significantly lower than the Ka/Ks value of 0.40 between B. napus and B. oleracea. A total of 35 homologous gene pairs between B. napus and its two progenitor species showed Ka/Ks ratios of more than 1, of which 15 homologous gene pairs showed Ka/Ks ratios of more than 1.2, suggesting positive selection for beneficial mutations. Gene ontology (GO) analysis of these genes revealed that the majority of these genes had molecular functionality in ADP binding properties and were involved in signal transduction and defense response processes.
Table 3. Ka/Ks ratio of homologous gene pairs for NBS-encoding genes between the genomes of B. napus and its progenitor species B. rapa and B. Oleracea.
Analyses of Phylogenetic Patterns of NBS-Encoding Genes Between B. napus and Its Progenitor Species
Whole genome duplication is thought to be a source of complexity and diversity for plant species that may allow them to adapt to changing environmental conditions. To investigate the phylogenetic relationships of NBS-encoding genes among genomes, the coding sequences of all NBS-encoding genes in B. napus, B. rapa, and B. oleracea were analyzed and compared, and a phylogenetic tree was constructed (Supplementary Figure S1). In the phylogenetic tree, at least six evolutionary patterns were observed.
The first evolutionary pattern observed was that of straightforward inheritance of genes from progenitor species B. rapa and B. oleracea to their allopolyploid progeny species B. napus. This evolutionary pattern was commonly observed for NBS-encoding genes in B. napus. For instance in one cluster (shown in Supplementary Figure S2a), seven NBS-encoding genes of B. napus individually grouped with eight homologous NBS-encoding genes from both B. rapa and B. oleracea. In the second evolutionary pattern, NBS-encoding genes in B. napus originated from a single progenitor species (Supplementary Figure S2b). In one example (Supplementary Figure S2b), 11 NBS-encoding genes in B. napus grouped with 13 homologous NBS-encoding genes from B. rapa. The third evolutionary pattern found was a large scale expansion of NBS-encoding genes in B. napus relative to the progenitor species. A typical gene cluster with this pattern is shown in Supplementary Figure S2c, where 13 NBS-encoding genes of B. napus were grouped with only one NBS-encoding gene from B. rapa. The fourth evolutionary pattern was the loss of NBS-encoding R-genes from progenitors in B. napus. For instance in one cluster shown in Supplementary Figure S2d, nine NBS-encoding genes from the progenitor species were grouped with only three NBS-encoding R-genes from B. napus. The fifth evolutionary pattern observed is that of different sequence divergence rates between NBS-encoding genes. Some gene clusters in the phylogenetic tree had longer branch lengths than others, which suggested different sequence divergence rates between gene clusters. For instance, in a cluster shown in Supplementary Figure S2e, one NBS-encoding gene of B. napus showed an extremely long branch length relative to the other genes in the cluster, suggesting a more rapid sequence divergence rate for this gene. A sixth (rare) evolutionary pattern of gene loss from the progenitor species to B. napus was also observed for two NBS-encoding genes (Supplementary Figure S1). As individual gene clusters commonly exhibited more than one evolutionary pattern, relative prevalence of each type was difficult to assess. However, we observed at least six patterns for the evolution of NBS-encoding genes in B. napus, suggesting a complicated evolutionary trajectory for NBS-encoding genes despite the relatively short history of this young allopolyploid species.
Comparative Expression of NBS-Encoding Genes in B. rapa, B. oleracea, and B. napus
As well as variation in gene copy number and in gene sequences, variation in gene expression is also an important aspect of gene evolution. Relative expression data for the NBS-encoding genes was collected from the root, stem, leaf, flower and seed tissues of B. rapa, B. oleracea, and B. napus (Supplementary Table S1). According to the gene expression level, 28 gene expression models were determined for the NBS-encoding genes in the three species (Table 4). The first five gene expression models encompassed 62.5% of NBS-encoding genes in B. napus (290 genes). The most common model was for genes to be expressed more highly in the root only, followed by higher gene expression in root and leaf, then higher expression in leaf tissue only, then higher expression in root and flower, then higher expression in stem and leaf. In B. rapa, these first five expression models covered 121 NBS-encoding genes (59.9%). The proportion of genes fitting each of the five expression models in B. rapa were sorted from larger to smaller as follows: higher expression in root only, higher expression in stem and leaf, higher expression in root and stem, higher expression in root, stem and leaf, and higher expression in root and seed. In B. oleracea, the first five expression models covered 105 NBS-encoding genes (71.9%). The proportion of genes fitting each of the five expression models were sorted from larger to smaller as follows: higher expression in root only, no higher expression in any of the five tissues, higher expression in root and leaf, higher expression in root and stem, and higher expression in root and flower.
Table 4. Expression patterns of NBS-encoding genes in B. napus (2n = AACC), B. rapa (2n = AA), and B. oleracea (2n = CC) root, stem, leaf, flower, and seed tissues.
The majority of the NBS-encoding genes showed higher expression levels in the root in all of the three species, accounting for 59.9, 59.6, and 59.7% of the NBS-encoding genes in B. rapa, B. oleracea, and B. napus (Figure 4). In contrast, the proportion of NBS-encoding genes with higher expression in the seed was low in all the three species, accounting for 15.4, 16.4, and 6.7% of the NBS-encoding genes in B. rapa, B. oleracea, and B. napus (Figure 4). These results revealed that the NBS-encoding genes were more likely to be highly expressed in root tissue, which reflects that the roots might show stronger stress responses (or be under greater disease stress) relative to the other tissues.
Figure 4. Expression analysis of NBS-encoding genes in tissues of B. rapa, B. Oleracea, and B. napus. The proportions represent higher expression of NBS-encoding genes in each tissue.
Comparative Analysis of the Identified NBS-Encoding Genes and Mapped Quantitative Resistance Loci of Three Major Diseases of B. napus
Blackleg (Leptosphaeria maculans and L. biglobosa), clubroot (Plasmodiophora brassicae), and Sclerotinia stem rot (Sclerotinia sclerotiorum) are three key pathogens affecting B. napus. Here, we collected mapped resistance QTLs against these three major disease-causing pathogens, to identify colocalization of the identified NBS-encoding genes and these QTLs. A total of 170 chromosome intervals distributed on the 19 chromosomes of B. napus, were obtained by BLAST analysis of the flanking markers of known QTLs against the reference genome of B. napus, and included 81 loci associated with blackleg resistance, 19 loci associated with clubroot resistance, and 70 loci associated with Sclerotinia stem rot resistance (Supplementary Table S2). Of these QTL mapped to chromosome intervals, 71 were co-located with NBS-encoding genes identified in the present study (204 genes; Supplementary Table S3). Each of these 71 intervals contained from 1 to 18 NBS-encoding genes (Supplementary Table S3).
Of the 204 NBS-encoding genes located within the mapped QTLs, 165 NBS-encoding genes could be mapped within QTLs for blackleg resistance, 51 NBS-encoding genes could be mapped within the QTLs for clubroot resistance, and 41 NBS-encoding genes could be mapped within the QTLs for Sclerotinia stem rot resistance (Supplementary Tables S3, S4). Hundred and fifty four genes were co-located with resistance QTLs against a single disease. Of which, 139 genes were co-located with a resistance QTL identified for a single disease, including 105, 18, and 16 genes co-located with resistance QTLs against blackleg, clubroot and Sclerotinia stem rot, respectively (Supplementary Table S4), and 15 genes were co-located with at least two overlapping QTLs for the same disease from different studies, including 13 genes within overlapping resistance QTLs against blackleg, and two genes within overlapping resistance QTLs against clubroot (Supplementary Table S4). A further 47 NBS-encoding genes were co-located with two overlapping resistance QTLs against different diseases, including 25 genes within resistance QTLs against both blackleg and clubroot, 19 genes within resistance QTLs against both blackleg and Sclerotinia stem rot, and three genes within resistance QTLs against both clubroot and Sclerotinia stem rot (Supplementary Table S4). Three genes were also found within overlapping QTLs against all three diseases (Supplementary Table S4).
Massive Copy Number Increase in Putative Resistance Genes in B. napus Relative to Its Progenitor Species B. rapa and B. oleracea, Particularly in the C Genome
Despite the fact that B. napus is a very recent allopolyploid, <10,000 years old (Chalhoub et al., 2014), we found extremely dynamic evolution of putatively functional resistance genes in this species relative to its progenitor species B. rapa and B. oleracea. In particular, an additional 116 new NBS-encoding genes were identified (an increase of 33%), the vast majority of which were observed in the B. napus C genome relative to the progenitor B. oleracea genome. Significant sequence divergence was also observed in 1/3 of the NBS-encoding genes, and clustering patterns were also changed in the B. napus C genome relative to in B. oleracea. This pattern may be related to the much higher fraction of repetitive sequences and transposable elements in the C genome relative to the A genome (Wang et al., 2011; Chalhoub et al., 2014; Liu et al., 2014; Parkin et al., 2014): resistance genes are often located within more highly repetitive regions and hence undergo more dynamic evolutionary processes as a result, as this facilitates generation of new, potentially useful phenotypic variation for disease resistance. Our results differ significantly from those of Alamery et al. (2017), who found similar numbers of NBS-LRR genes in B. napus relative to its progenitors as well as many more genes. However, as we used functional protein coding domains as search queries and ignored pseudogenes, it is likely that the differences between our studies can be explained by this difference in approach. In support of this, 43% of NBS-LRR genes identified by Alamery et al. (2017) were thought to be non-functional in this study. Hence, our results are more likely to reflect functional evolutionary patterns for disease resistance genes in the B. napus genome.
NBS-Encoding Genes Were Preferentially Expressed in Root Tissues in All Three Species
Roots are an essential organ for a myriad of physiological processes, including water and nutrient uptake, anchoring of plants, storage of assimilates and mechanical support of the above-ground organs. Roots are constantly in contact with the soil microbial community, which consists of mostly beneficial but also pathogenic organisms. Thus, plants have evolved root-specific defense responses against soil-borne pathogens, viruses, pests and oomycetes, including pre-formed physical and chemical barriers and inducible defenses (Dixon, 2001;De Coninck et al., 2015). NBS-encoding genes are the largest class of known disease resistance genes for defense against various pathogens, including bacteria, fungi, oomycetes, viruses, and nematodes. Therefore, it is not surprising that abundant NBS-encoding genes may actually be involved in beneficial root–pathogen interactions in the soil (Okubara and Paulitz, 2005). For instance, the tomato I-2 gene, a member of the NBS-LRR family, confers resistance to the wilt pathogen Fusarium oxysporum f.sp. lycopersici (Mes et al., 2000), and is preferentially expressed in lateral root primordia of young roots and in vascular regions of mature roots (Okubara and Paulitz, 2005). Two NBS-LRR genes (Mi-1 and Hero A) from tomato and two NBS-LRR genes (Gpa2 and Gro1-4) from potato individually confer root-specific resistance against root knot nematodes and cyst nematodes (Williamson and Kumar, 2006). Our study found that NBS-encoding genes were predominantly expressed in root tissue in each of B. rapa, B. oleracea, and B. napus, which is congruent with the greater pressure from micro-organisms on the root than on the other above-ground organs. Future studies may elucidate whether these genes are involved in disease resistance or in beneficial interactions with soil micro-organisms.
Diverse Patterns Explain the Rapid Evolution of NBS-Encoding Genes in B. napus After Formation From Its Diploid Progenitor Species
Phylogenetic analysis of coding sequences of NBS-encoding genes in the three species revealed at least six distinct evolutionary patterns, falling into three groups. We summarize these evolutionary patterns in Figure 5. The first group represented the straightforward inheritance of NBS-encoding genes directly from the progenitor species to B. napus (the first two patterns identified in the results). Within this group, NBS-encoding genes in B. napus were inherited directly from the two progenitor species, either as one or more copies from both progenitors or as one or more copies from only one of the two progenitor species. Given that B. rapa and B. oleracea diverged from each other about 1.44 Mya (Arias et al., 2014), followed by the very recent formation of B. napus as an allopolyploid species < 10,000 years ago (Chalhoub et al., 2014), it is not surprising that these two patterns of direct inheritance from both progenitor species as well as inheritance of species-specific genes were both commonly observed.
Figure 5. The evolutionary routes of NBS-encoding genes in B. napus relative to its two progenitor species B. rapa and B. oleracea.
The second group represented the expansion and contraction of NBS-encoding genes in B. napus relative to its two progenitor species (evolutionary patterns 3 and 4 in the Results section). In many species, such as grapevine, poplar, rice, and legumes, large scale expansions and contractions of NBS-encoding genes have also been found (Zhou et al., 2004; Yang et al., 2008a,b; Zheng et al., 2016). Overall, previous research suggests that plant resistance genes have very rapid expansion and contraction rates, as they evolve in response to interactions with pathogens (Faulkner and Robatzek, 2012; Nepal and Benson, 2015). Hence, the expansion and contraction of NBS-encoding genes in B. napus is expected to facilitate rapid evolution of resistance, particularly as a response to domestication.
The third group represented the birth and death of NBS-encoding genes in B. napus relative to its two progenitor species, including evolutionary pattern 5 from the results, where genes experienced large-scale sequence variation (effectively generating a novel gene), and evolutionary pattern 6 from the results, where genes from the progenitor species disappeared in B. napus. As an example of putative gene birth, two genes in the phylogenetic tree, GSBRNA2T00096083001 and GSBRNA2T00117236001, showed extreme sequence divergence from the other genes in the corresponding cluster. Sequence comparison of these two genes with the best matching genes in progenitor species revealed that these two genes combined sequences from several genes in the progenitor species, including NBS-encoding genes providing the NBS domains and non-NBS-encoding genes providing the other sequences, to generate new NBS-encoding genes in B. napus (Figure 6A). A putative mechanism for this is non-homologous recombination, either as a result of ancestral homoeology relationships (A-C or from triplicated regions within the A and C genomes; Parkin et al., 2014) or ectopically via the gene copies themselves as a substrate. Two genes (Bra000758 and Bol030521 from progenitor species), were also found to have been separated from their NBS-encoding domains in B. napus, indicating the “death” of these two genes in B. napus. Sequence comparison of these two genes with the best match genes in B. napus revealed that nearly half of the sequence of these two genes was retained (from approximately 1700 bp to the end of the gene), but that the NBS domain had disappeared (Figure 6B). These results suggest that both gene birth and death for NBS- encoding genes in B. napus are based on processes affecting original NBS-encoding genes, which are either subjected to large scale sequence variation to generate new genes or domain deletion resulting in gene death, with both processes putatively mediated by non-homologous recombination.
Figure 6. Sequence comparison of novel NBS-encoding genes in B. napus with the best matching genes in its progenitor species B. rapa and B. oleracea, and missing genes with the best matching genes in B. napus. (A) Sequence comparison of two novel genes in B. napus (GSBRNA2T00096083001 and GSBRNA2T00117236001) with the best matching genes in B. rapa and B. oleracea. (B) Sequence comparison of two missing genes (Bra000758 and Bol030521) with the best matching genes in B. napus. The colorful bars represent homologous blocks. The digits on the bars indicate the beginning and end of each homologous block.
A caveat to our results is that single genome sequences of B. rapa, B. oleracea, and B. napus are being used to infer species and population-wide evolution patterns; a shortfall of many current genomic studies where multiple reference genomes are not yet available for individual species (Golicz et al., 2015). As the exact progenitor genotypes of B. napus are also unknown (Chalhoub et al., 2014), it is also difficult to predict with certainty which genomic variants are representative of evolutionary patterns post-allopolyploid formation, and which may have been inherited from different genotypes of the diploid progenitors than those used to create the genome references. This may also explain some of the results where the best gene match to NBS-encoding genes in B. napus was outside the Brassica genus: although horizontal gene transfer is a possibility, distinguishing this from incomplete lineage sorting with currently available genomic resources is challenging. Genome sequence quality can also play a major role: although the Brassica genome assemblies are predicted to cover most of the gene space (Wang et al., 2011;Chalhoub et al., 2014; Parkin et al., 2014), these genome references are still incomplete, which can significantly affect the numbers of genes detected. Nevertheless, our results predict some interesting patterns for resistance gene evolution in B. napus, forming a basis for further investigation pending additional genomic resources in this young allopolyploid crop.
More NBS-Encoding Genes Were Identified in B. rapa Than in Close Relative B. oleracea
Many more NBS-encoding genes were found in B. rapa than in B. oleracea. It has been proposed previously by molecular phylogenies that a clade containing B. oleracea, B. rapa and their respective wild relatives originated 6.47 Mya, late in the Miocene, and that finally B. oleracea varieties and their wild relatives diversified from this clade 1.44 Mya (Arias et al., 2014). The later diversification of B. oleracea relative to B. rapa might offer an explanation for why the Co genome of B. oleracea, despite being significantly larger (630 Mb) than the Ar genome of B. rapa (485 Mb), contained fewer functional NBS-encoding genes. Additionally, B. oleracea varieties and their wild relatives were found to diversified in the northeastern Mediterranean and later spread through the rest of Europe (Sinskaia, 1928; Arias et al., 2014), while multiple divergence regions have been proposed for B. rapa, including Europe, central China, west Asia and north Africa (Zohary et al., 2012; Annisa et al., 2013; Guo et al., 2014). These multiple centers for divergence in B. rapa might provide another explanation for the greater prevalence of putatively functional NBS-genes in the Ar-genome relative to the Co genome: these may have evolved in response to the need to adapt to variable environments and disease pressures, and later spread between subpopulations. In future, with the availability of more, high-quality genome references, it may be interesting to compare NBS-encoding gene divergence between different genotypes and cultivars of B. rapa and B. oleracea (as well as B. napus), in order to fully elucidate the different evolutionary trajectories of resistance genes in response to genetic diversification, domestication and other genomic and environmental selective pressures.
A High Proportion of Mapped QTLs Against the Three Major Diseases of B. napus Were Repeatedly Identified and Showed Close Associations With Predicted NBS-Encoding Genes
We reviewed current studies covering QTL mapping of blackleg, clubroot, and Sclerotinia stem rot, the three major diseases of B. napus. We found that QTLs of these three diseases were unevenly distributed on chromosomes of B. napus, with the number of QTLs ranging from four to 14 (Table 5). This uneven distribution of QTLs across the chromosomes of B. napus was similar to the uneven distribution of NBS-encoding genes we observed. The three chromosomes with the most QTLs detected for resistance against these pathogens were An09 with 14 QTLs, Cn3 with 14 QTLs, and Cn4 with 13 QTLs; this coincides with the finding that An09 and Cn3 were two of the three chromosomes with the most NBS-encoding genes (Figure 2). The similar distribution of disease resistance QTLs and NBS-encoding genes supports the close associations of NBS encoding genes with disease resistance against blackleg, clubroot and Sclerotinia stem rot in B. napus.
We also investigated overlapping and neighboring (within 1 Mb) QTL intervals to find repeatedly identified QTLs: approximately half (42/81) of resistance QTLs against blackleg were repeatedly identified. A total of 13 of these loci were identified for blackleg resistance (on An01, An02, An08, An09, An10, Cn1, Cn2, Cn4, Cn5, Cn6, and Cn8). One locus on chromosome An03 contained 4 of 19 resistance QTLs against clubroot, most likely corresponding to the major resistance gene CRa previously identified in this region (Ueno et al., 2012). Of the Sclerotinia stem rot resistance QTLs, 20 of 70 were repeatedly identified in 8 loci on An02, An06, An09, Cn2, Cn4, Cn6, and Cn9 (Table 5 and Supplementary Table S2). Co-localized QTL detected across multiple studies are highly credible candidates for resistance against these three diseases. Across all QTL, 41.76% of intervals were co-located with NBS-encoding genes, which verified the close associations between the NBS-encoding genes and QTLs against these pathogens. Three NBS-encoding genes, (GSBRNA2T00081287001, GSBRNA2T00021121001, and GSBRNA2T00098068001) were located within overlapping QTLs against all three diseases, while 47 NBS-encoding genes were located within overlapping QTLs against two diseases. These results suggested that several common factors might exist for B. napus that confer resistance against clubroot, blackleg and Sclerotinia stem rot, and that these key factors may comprise NBS-encoding genes.
YF and DF designed the research. YF, YZ, BL, DZ, and HY analyzed the data. YF wrote the manuscript. AM assisted with interpretation of results, manuscript writing, and revision.
This study was supported financially by Project of creative agriculture derived R&D and demonstrations of horticultural varieties (2018C02057), The Key Project of Novel Variety Breeding of Zhejiang Province (2016C02050-8), The Key Project of Variety Breeding of Seven Major Crops (2016YFD0101306), and Project of mechanized rapeseed Breeding of the downstream of Yangtze River (2018YFD0100602). AM is supported by an Emmy Noether DFG grant (MA6473/1-1).
Conflict of Interest Statement
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fpls.2019.00026/full#supplementary-material
Alamery, S., Tirnaz, S., Bayer, P., Alamery, S., Tirnaz, S., Bayer, P., et al. (2017). Genome-wide identification and comparative analysis of NBS-LRR resistance genes in Brassica napus. Crop Pasture Sci. 69, 79–94. doi: 10.1071/CP17214
Ameline-Torregrosa, C., Wang, B. B., O’Bleness, M. S., Deshpande, S., Zhu, H., Roe, B., et al. (2008). Identification and characterization of nucleotide-binding site-leucine-rich repeat genes in the model plant Medicago truncatula. Plant Physiol. 146, 5–21. doi: 10.1104/pp.107.104588
Andolfo, G., Sanseverino, W., Aversano, R., Frusciante, L., and Ercolano, M. R. (2014). Genome-wide identification and analysis of candidate genes for disease resistance in tomato. Mol. Breed. 33, 227–233. doi: 10.1007/s11032-013-9928-7
Arias, T., Beilstein, M. A., Tang, M., McKain, M. R., and Pires, J. C. (2014). Diversification times among Brassica (Brassicaceae) crops suggest hybrid formation after 20 million years of divergence. Am. J. Bot. 101, 86–91. doi: 10.3732/ajb.1300312
Ariyarathna, H. C. K., and Francki, M. G. (2016). Phylogenetic relationships and protein modelling revealed two distinct subfamilies of group II HKT genes between crop and model grasses. Genome 59, 509–517. doi: 10.1139/gen-2016-0035
Arya, P., Kumar, G., Acharya, V., and Singh, A. K. (2014). Genome-wide identification and expression analysis of NBS-encoding genes in Malus x domestica and expansion of NBS genes family in Rosaceae. PLoS One 9:e107987. doi: 10.1371/journal.pone.0107987
Bai, J., Pennill, L. A., Ning, J., Lee, S. W., Ramalingam, J., Webb, C. A., et al. (2002). Diversity in nucleotide binding site–leucine-rich repeat genes in cereals. Genome Res. 12, 1871–1884. doi: 10.1101/gr.454902
Bayer, P. E., Golicz, A. A., Tirnaz, S., Chan, C. K., Edwards, D., and Batley, J. (2018). Variation in abundance of predicted resistance genes in the Brassica oleracea pangenome. Plant Biotechnol. J. doi: 10.1111/pbi.13015 [Epub ahead of print].
Cai, G., Yang, Q., Yang, Q., Zhao, Z., Chen, H., Wu, J., et al. (2012). Identification of candidate genes of QTLs for seed weight in Brassica napus through comparative mapping among Arabidopsis and Brassica species. BMC Genet. 13:105. doi: 10.1186/1471-2156-13-105
Cannon, S. B., Zhu, H., Baumgarten, A. M., Spangler, R., May, G., Cook, D. R., et al. (2002). Diversity, distribution, and ancient taxonomic relationships within the TIR and non-TIR NBS-LRR resistance gene subfamilies. J. Mol. Evol. 54, 548–562. doi: 10.1007/s0023901-0057-2
Chalhoub, B., Denoeud, F., Liu, S., Parkin, I. A., Tang, H., Wang, X., et al. (2014). Early allopolyploid evolution in the post-Neolithic Brassica napus oilseed genome. Science 345, 950–953. doi: 10.1126/science.1253435
Chen, Q., Han, Z., Jiang, H., Tian, D., and Yang, S. (2010). Strong positive selection drives rapid diversification of R-genes in Arabidopsis relatives. J. Mol. Evol. 70, 137–148. doi: 10.1007/s00239-009-9316-4
Conesa, A., Götz, S., García-Gómez, J. M., Terol, J., Talón, M., and Robles, M. (2005). Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics 21, 3674–3676. doi: 10.1093/bioinformatics/bti610
De Coninck, B., Timmermans, P., Vos, C., Cammue, B. P., and Kazan, K. (2015). What lies beneath: belowground defense strategies in plants. Trends Plant Sci. 20, 91–101. doi: 10.1016/j.tplants.2014.09.007
Delourme, R., Piel, N., Horvais, R., Pouilly, N., Domin, C., Vallée, P., et al. (2008). Molecular and phenotypic characterization of near-isogenic lines at QTL for quantitative resistance to Leptosphaeria maculans, in oilseed rape (Brassica napus, l.). Theor. Appl. Genet. 117, 1055–1067. doi: 10.1007/s00122-008-0844-x
Eisen, M. B., Spellman, P. T., Brown, P. O., and Botstein, D. (1998). Cluster analysis and display of genome-wide expression patterns. Proc. Natl. Acad. Sci. U.S.A. 95, 14863–14868. doi: 10.1073/pnas.95.25.14863
Golicz, A. A., Martinez, P. A., Zander, M., Patel, D. A., Van De Wouw, A. P., Visendi, P., et al. (2015). Gene loss in the fungal canola pathogen Leptosphaeria maculans. Funct. Integr. Genomics 15, 189–196. doi: 10.1007/s10142-014-0412-1
Guo, Y. L., Fitz, J., Schneeberger, K., Ossowski, S., Cao, J., and Weigel, D. (2011). Genome-wide comparison of nucleotide-binding site leucine-rich repeat-encoding genes in Arabidopsis. Plant Physiol. 157, 757–769. doi: 10.1104/pp.111.181990
Hirai, M., Harada, T., Kubo, N., Tsukada, M., Suwabe, K., and Matsumoto, S. (2004). A novel locus for Clubroot resistance in Brassica rapa and its linkage markers. Theor. Appl. Genet. 108, 639–643. doi: 10.1007/s00122-003-1475-x
Jian, W., Cai, G., Tu, J., Li, L., Sheng, L., Luo, X., et al. (2013). Identification of QTLS for resistance to Sclerotinia Stem Rot and Bnac. IGMT5.a as a candidate gene of the major resistant QTL SRC6 in Brassica napus. PLoS One 8:e67740. doi: 10.1371/journal.pone.0067740
Jupe, F., Pritchard, L., Etherington, G. J., MacKenzie, K., Cock, P. J., Wright, F., et al. (2012). Identification and localisation of the NB-LRR gene family within the potato genome. BMC Genet. 13:75. doi: 10.1186/1471-2164-13-75
Jupe, F., Witek, K., Verweij, W., Śliwka, J., Pritchard, L., Etherington, G. J., et al. (2013). Resistance gene enrichment sequencing (RenSeq) enables reannotation of the NB-LRR gene family from sequenced plant genomes and rapid mapping of resistance loci in segregating populations. Plant J. 76, 530–544. doi: 10.1111/tpj.12307
Kato, T., Hatakeyama, K., Fukino, N., and Matsumoto, S. (2013). Fine mapping of the Clubroot resistance gene CRb and development of a useful selectable marker in Brassica rapa. Breed. Sci. 63, 116–124. doi: 10.1270/jsbbs.63.116
Kaur, S., Cogan, N. O. I., Ye, G., Baillie, R. C., Hand, M. L., Ling, A. E., et al. (2009). Genetic map construction and QTL mapping of resistance to blackleg (Leptosphaeria maculans) disease in Australian canola (Brassica napus L.) cultivars. Theor. Appl. Genet. 120, 71–83. doi: 10.1007/s00122-009-1160-9
Kim, S., Park, M., Yeom, S. I., Kim, Y. M., Lee, J. M., Lee, H. A., et al. (2014). Genome sequence of the hot pepper provides insights into the evolution of pungency in Capsicum species. Nat. Genet. 46, 270–278. doi: 10.1038/ng.2877
Krzywinski, M., Schein, J., Birol, I., Connors, J., Gascoyne, R., Horsman, D., et al. (2009). Circos: an information aesthetic for comparative genomics. Genome Res. 19, 1639–1645. doi: 10.1101/gr.092759.109
Kumar, V., Paillard, S., Fopa-Fomeju, B., Falentin, C., Deniot, G., Baron, C., et al. (2018). Multi-year linkage and association mapping confirm the high number of genomic regions involved in oilseed rape quantitative resistance to blackleg. Theor. Appl. Genet. 131, 1627–1643. doi: 10.1007/s00122-018-3103-9
Larkan, N. J., Raman, H., Lydiate, D. J., Robinson, S. J., Yu, F., Barbulescu, D. M., et al. (2016). Multi-environment QTL studies suggest a role for cysteine-rich protein kinase genes in quantitative resistance to blackleg disease in Brassica napus. BMC Plant Biol. 16:183. doi: 10.1186/s12870-016-0877-2
Leister, R. T., and Katagiri, F. (2000). A resistance gene product of the nucleotide binding site-leucine rich repeats class can form a complex with bacterial avirulence proteins in vivo. Plant J. 22, 345–354. doi: 10.1046/j.1365-313x.2000.00744.x
Li, J., Ding, J., Zhang, W., Zhang, Y., Tang, P., Chen, J. Q., et al. (2010). Unique evolutionary pattern of numbers of gramineous NBS–LRR genes. Mol. Genet. Genomics 283, 427–438. doi: 10.1007/s00438-010-0527-6
Li, L., Luo, Y., Chen, B., Xu, K., Zhang, F., Li, H., et al. (2016). A genome-wide association study reveals new loci for resistance to Clubroot disease in Brassica napus. Front. Plant Sci. 7:1483. doi: 10.3389/fpls.2016.01483
Li, X., Cheng, Y., Ma, W., Zhao, Y., Jiang, H., and Zhang, M. (2010). Identification and characterization of NBS-encoding disease resistance genes in Lotus japonicus. Plant Syst. Evol. 289, 101–110. doi: 10.1007/s00606-010-0331-0
Lin, X., Zhang, Y., Kuang, H., and Chen, J. (2013). Frequent loss of lineages and deficient duplications accounted for low copy number of disease resistance genes in Cucurbitaceae. BMC Genomics 14:335. doi: 10.1186/1471-2164-14-335
Liu, J., Liu, X., Dai, L., and Wang, G. (2007). Recent progress in elucidating the structure, function and evolution of disease resistance genes in plants. J. Genet. Genomics 34, 765–776. doi: 10.1016/S1673-8527(07)60087-3
Liu, S., Liu, Y., Yang, X., Tong, C., Edwards, D., Parkin, I. A., et al. (2014). The Brassica oleracea genome reveals the asymmetrical evolution of polyploid genomes. Nat. Commun. 5:3930. doi: 10.1038/ncomms4930
Lozano, R., Ponce, O., Ramirez, M., Mostajo, N., and Orjeda, G. (2012). Genome-wide identification and mapping of NBS-encoding resistance genes in Solanum tuberosum group phureja. PLoS One 7:e34775. doi: 10.1371/journal.pone.0034775
Luo, S., Zhang, Y., Hu, Q., Chen, J., Li, K., Lu, C., et al. (2012). Dynamic nucleotide-binding site and leucine-rich repeat-encoding genes in the grass family. Plant Physiol. 159, 197–210. doi: 10.1104/pp.111.192062
Martin, G. B., Bogdanove, A. J., and Sessa, G. (2003). Understanding the functions of plant disease resistance proteins. Annu. Rev. Plant Biol. 54, 23–61. doi: 10.1146/annurev.arplant.54.031902.135035
Mayer, K. F., Martis, M., Hedley, P. E., Šimková, H., Liu, H., Morris, J. A., et al. (2011). Unlocking the barley genome by chromosomal and comparative genomics. Plant Cell 23, 1249–1263. doi: 10.1105/tpc.110.082537
Mei, J., Lu, K., Wei, D., Liu, Y., Disi, J. O., Li, J., et al. (2013). Identification of genomic regions involved in resistance against Sclerotinia sclerotiorum from wild Brassica oleracea. Theor. Appl. Genet. 126, 549–556. doi: 10.1007/s00122-012-2000-x
Mes, J. J., Van Doorn, A. A., Wijbrandi, J., Simons, G., Cornelissen, B. J., and Haring, M. A. (2000). Expression of the Fusarium resistance gene I-2 colocalizes with the site of fungal containment. Plant J. 23, 183–193. doi: 10.1046/j.1365-313x.2000.00765.x
Meyers, B. C., Dickerman, A. W., Michelmore, R. W., Sivaramakrishnan, S., Sobral, B. W., and Young, N. D. (1999). Plant disease resistance genes encode members of an ancient and diverse protein family within the nucleotide binding superfamily. Plant J. 20, 317–332. doi: 10.1046/j.1365-313X.1999.t01-1-00606.x
Miao, L., Zhang, L., Raboanatahiry, N., Lu, G., Zhang, X., Xiang, J., et al. (2016). Transcriptome analysis of stem and globally comparison with other tissues in Brassica napus. Front. Plant Sci. 7:1403. doi: 10.3389/fpls.2016.01403
Mun, J. H., Yu, H. J., Park, S., and Park, B. S. (2009). Genome-wide identification of NBS-encoding resistance genes in Brassica rapa. Mol. Genet. Genomics 282, 617–631. doi: 10.1007/s00438-009-0492-0
Nagaoka, T., Doullah, M. A. U., Matsumoto, S., Kawasaki, S., Ishikawa, T., Hori, H., et al. (2010). Identification of QTLS that control clubroot resistance in Brassica oleracea and comparative analysis of Clubroot resistance genes between B. rapa and B. oleracea. Theor. Appl. Genet. 120, 1335–1346. doi: 10.1007/s00122-010-1259-z
Pang, W., Fu, P., Li, X., Zhan, Z., Yu, S., and Piao, Z. (2018). Identification and mapping of the Clubroot resistance gene CRd in Chinese cabbage (Brassica rapa ssp. pekinensis). Front. Plant Sci. 9:653. doi: 10.3389/fpls.2018.00653
Parkin, I. A., Koh, C., Tang, H., Robinson, S. J., Kagale, S., Clarke, W. E., et al. (2014). Transcriptome and methylome profiling reveals relics of genome dominance in the mesopolyploid Brassica oleracea. Genome Biol. 15:R77. doi: 10.1186/gb-2014-15-6-r77
Peele, H. M., Guan, N., Fogelqvist, J., and Dixelius, C. (2014). Loss and retention of resistance genes in five species of the Brassicaceae family. BMC Plant Biol. 14:298. doi: 10.1186/s12870-014-0298-z
Porter, B. W., Paidi, M., Ming, R., Alam, M., Nishijima, W. T., and Zhu, Y. J. (2009). Genome-wide analysis of Carica papaya reveals a small NBS resistance gene family. Mol. Genet. Genomics 281, 609–626. doi: 10.1007/s00438-009-0434-x
Qian, L. H., Zhou, G. C., Sun, X. Q., Lei, Z., Zhang, Y. M., Xue, J. Y., et al. (2017). Distinct patterns of gene gain and loss: diverse evolutionary modes of NBS-encoding genes in three Solanaceae crop species. G3 7, 1577–1585. doi: 10.1534/g3.117.040485
Raman, R., Taylor, B., Marcroft, S., Stiller, J., Eckermann, P., Coombes, N., et al. (2012). Molecular mapping of qualitative and quantitative loci for resistance to Leptosphaeria maculans causing blackleg disease in canola (Brassica napus L.). Theor. Appl. Genet. 125, 405–418 doi: 10.1007/s00122-012-1842-6.
Richly, E., Kurth, J., and Leister, D. (2002). Mode of amplification and reorganization of resistance genes during recent Arabidopsis thaliana evolution. Mol. Biol. Evol. 19, 76–84. doi: 10.1093/oxfordjournals.molbev.a003984
Sakamoto, K., Saito, A., Hayashida, N., Taguchi, G., and Matsumoto, E. (2008). Mapping of isolate-specific QTLs for Clubroot resistance in Chinese cabbage (Brassica rapa L. ssp. pekinensis). Theor. Appl. Genet. 117, 759–767. doi: 10.1007/s00122-008-0817-0
Sanjaya, G., Myrtle, H., Jonathan, D., Kyla, H., Parkin, I. A. P., Hegedus, D. D., et al. (2016). Microsatellite markers used for genome-wide association mapping of partial resistance to Sclerotinia sclerotiorum in a world collection of Brassica napus. Mol. Breed. 36, 1–13.
Shao, Z. Q., Zhang, Y. M., Hang, Y. Y., Xue, J. Y., Zhou, G. C., Wu, P., et al. (2014). Long-term evolution of nucleotide-binding site-leucine-rich repeat genes: understanding gained from and beyond the legume family. Plant Physiol. 166, 217–234. doi: 10.1104/pp.114.243626
Suwabe, K., Tsukazaki, H., Iketani, H., Hatakeyama, K., Kondo, M., Fujimura, M., et al. (2006). Simple sequence repeat-based comparative genomics between Brassica rapa and Arabidopsis thaliana: the genetic origin of Clubroot resistance. Genetics 173, 309–319. doi: 10.1534/genetics.104.038968
Thompson, J. D., Higgins, D. G., and Gibson, T. J. (1994). CLUSTAL W: improving the sensitivity of progressive multiple sequence alignment through sequence weighting, position-specific gap penalties and weight matrix choice. Nucleic Acids Res. 22, 4673–4680. doi: 10.1093/nar/22.22.4673
Tong, C., Wang, X., Yu, J., Wu, J., Li, W., Huang, J., et al. (2013). Comprehensive analysis of RNA-seq data reveals the complexity of the transcriptome in Brassica rapa. BMC Genomics 14:689. doi: 10.1186/1471-2164-14-689
Ueno, H., Matsumoto, E., Aruga, D., Kitagawa, S., Matsumura, H., and Hayashida, N. (2012). Molecular characterization of the CRa gene conferring clubroot resistance in Brassica rapa. Plant Mol. Biol. 80, 621–629. doi: 10.1007/s11103-012-9971-5
Wan, H., Yuan, W., Bo, K., Shen, J., Pang, X., and Chen, J. (2013). Genome-wide analysis of NBS-encoding disease resistance genes in Cucumis sativus and phylogenetic study of NBS-encoding genes in Cucurbitaceae crops. BMC Genomics 14:109. doi: 10.1186/1471-2164-14-109
Wei, D., Mei, J., Fu, Y., Disi, J. O., Li, J., and Qian, W. (2014). Quantitative trait loci analyses for resistance to Sclerotinia sclerotiorum, and flowering time in Brassica napus. Mol. Breed. 34, 1797–1804. doi: 10.1007/s11032-014-0139-7
Wei, L., Jian, H., Lu, K., Filardo, F., Yin, N., Liu, L., et al. (2016). Genome-wide association analysis and differential expression analysis of resistance to Sclerotinia Stem Rot in Brassica napus. Plant Biotechnol. J. 14, 1368–1380. doi: 10.1111/pbi.12501
Wicker, T., Yu, Y., Haberer, G., Mayer, K. F., Marri, P. R., Rounsley, S., et al. (2016). DNA transposon activity is associated with increased mutation rates in genes of rice and other grasses. Nat. Commun. 7:12790. doi: 10.1038/ncomms12790
Wu, J., Zhao, Q., Liu, S., Shahid, M., Lan, L., Cai, G., et al. (2016). Genome-wide association study identifies new loci for resistance to Sclerotinia Stem Rot in Brassica napus. Front. Plant Sci. 7:1418. doi: 10.3389/fpls.2016.01418
Wu, P., Shao, Z. Q., Wu, X. Z., Wang, Q., Wang, B., Chen, J. Q., et al. (2014). Loss/retention and evolution of NBS-encoding genes upon whole genome triplication of Brassica rapa. Gene 540, 54–61. doi: 10.1016/j.gene.2014.01.082
Xiang, L., Liu, J., Wu, C., Deng, Y., Cai, C., Zhang, X., et al. (2017). Genome-wide comparative analysis of NBS-encoding genes in four Gossypium species. BMC Genomics 18:292. doi: 10.1186/s12864-017-3682-x
Yang, S., Feng, Z., Zhang, X., Jiang, K., Jin, X., Hang, Y., et al. (2006). Genome-wide investigation on the genetic variations of rice disease resistance genes. Plant Mol. Biol. 62, 181–193. doi: 10.1007/s11103-006-9012-3
Yang, S., Zhang, X., Yue, J. X., Tian, D., and Chen, J. Q. (2008b). Recent duplications dominate NBS-encoding gene expansion in two woody species. Mol. Genet. Genomics 280, 187–198. doi: 10.1007/s00438-008-0355-0
Yang, S., Li, J., Zhang, X., Zhang, Q., Huang, J., Chen, J. Q., et al. (2013). Rapidly evolving R genes in diverse grass species confer resistance to rice blast disease. Proc. Natl. Acad. Sci. U.S.A. 110, 18572–18577. doi: 10.1073/pnas.1318211110
Yu, J., Tehrim, S., Zhang, F., Tong, C., Huang, J., Cheng, X., et al. (2014). Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana. BMC Genomics 15:3. doi: 10.1186/1471-2164-15-3
Yue, J. X., Meyers, B. C., Chen, J. Q., Tian, D., and Yang, S. (2012). Tracing the origin and evolutionary history of plant nucleotide-binding site–leucine-rich repeat (NBS-LRR) genes. New Phytol. 193, 1049–1063. doi: 10.1111/j.1469-8137.2011.04006.x
Zhang, Y. M., Shao, Z. Q., Wang, Q., Hang, Y. Y., Xue, J. Y., Wang, B., et al. (2016). Uncovering the dynamic evolution of nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes in Brassicaceae. J. Integr. Plant Biol. 58, 165–177. doi: 10.1111/jipb.12365
Zhao, J., Udall, J. A., Quijada, P. A., Grau, C. R., Meng, J., and Osborn, T. C. (2006). Quantitative trait loci for resistance to Sclerotinia sclerotiorum and its association with a homeologous non-reciprocal transposition in Brassica napus L. Theor. Appl. Genet. 112, 509–516. doi: 10.1007/s00122-005-0154-5
Zheng, F., Wu, H., Zhang, R., Li, S., He, W., Wong, F. L., et al. (2016). Molecular phylogeny and dynamic evolution of disease resistance genes in the legume family. BMC Genomics 17:402. doi: 10.1186/s12864-016-2736-9
Zhou, T., Wang, Y., Chen, J. Q., Araki, H., Jing, Z., Jiang, K., et al. (2004). Genome-wide identification of NBS genes in japonica rice reveals significant expansion of divergent non-TIR NBS-LRR genes. Mol. Genet. Genomics 271, 402–415. doi: 10.1007/s00438-004-0990-z
Zohary, D., Hopf, M., and Weiss, E. (2012). Domestication of Plants in the Old World: the Origin and Spread of Domesticated Plants in Southwest Asia, Europe, and the Mediterranean Basin. Oxford: Oxford University Press. doi: 10.1093/acprof:osobl/9780199549061.001.0001
Keywords: nucleotide-binding site (NBS), disease resistance, Brassica napus, Brassica oleracea, Brassica rapa
Citation: Fu Y, Zhang Y, Mason AS, Lin B, Zhang D, Yu H and Fu D (2019) NBS-Encoding Genes in Brassica napus Evolved Rapidly After Allopolyploidization and Co-localize With Known Disease Resistance Loci. Front. Plant Sci. 10:26. doi: 10.3389/fpls.2019.00026
Received: 02 July 2018; Accepted: 09 January 2019;
Published: 30 January 2019.
Edited by:Carl J. Rothfels, University of California, Berkeley, United States
Reviewed by:Tong Zhang, Pacific Northwest National Laboratory (DOE), United States
Christian Jung, Christian-Albrechts-Universität zu Kiel, Germany
Copyright © 2019 Fu, Zhang, Mason, Lin, Zhang, Yu and Fu. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.