Genome-wide analysis of NBS-LRR-encoding genes in Arabidopsis

The Arabidopsis genome contains approximately 200 genes that encode proteins with similarity to the nucleotide binding site and other domains characteristic of plant resistance proteins. Through a reiterative process of sequence analysis and reannotation, we identified 149 NBS-LRR-encoding genes in the Arabidopsis (ecotype Columbia) genomic sequence. Fifty-six of these genes were corrected from earlier annotations. At least 12 are predicted to be pseudogenes. As described previously, two distinct groups of sequences were identified: those that encoded an N-terminal domain with Toll/Interleukin-1 Receptor homology (TIR-NBS-LRR, or TNL), and those that encoded an N-terminal coiled-coil motif (CC-NBS-LRR, or CNL). The encoded proteins are distinct from the 58 predicted adapter proteins in the previously described TIR-X, TIR-NBS, and CC-NBS groups. Classification based on protein domains, intron positions, sequence conservation, and genome distribution defined four subgroups of CNL proteins, eight subgroups of TNL proteins, and a pair of divergent NL proteins that lack a defined N-terminal motif. CNL proteins generally were encoded in single exons, although two subclasses were identified that contained introns in unique positions. TNL proteins were encoded in modular exons, with conserved intron positions separating distinct protein domains. Conserved motifs were identified in the LRRs of both CNL and TNL proteins. In contrast to CNL proteins, TNL proteins contained large and variable C-terminal domains. The extant distribution and diversity of the NBS-LRR sequences has been generated by extensive duplication and ectopic rearrangements that involved segmental duplications as well as microscale events. The observed diversity of these NBS-LRR proteins indicates the variety of recognition molecules available in an individual genotype to detect diverse biotic challenges.     

Genome-wide identification of NBS-encoding resistance genes in <Emphasis Type="Italic">Brassica rapa</Emphasis>

Nucleotide-binding site (NBS)-encoding resistance genes are key plant disease-resistance genes and are abundant in plant genomes, comprising up to 2% of all genes. The availability of genome sequences from several plant models enables the identification and cloning of NBS-encoding genes from closely related species based on a comparative genomics approach. In this study, we used the genome sequence of Brassica rapa to identify NBS-encoding genes in the Brassica genome. We identified 92 non-redundant NBS-encoding genes [30 CC-NBS-LRR (CNL) and 62 TIR-NBS-LRR (TNL) genes] in approximately 100 Mbp of B. rapa euchromatic genome sequence. Despite the fact that B. rapa has a significantly larger genome than Arabidopsis thaliana due to a recent whole genome triplication event after speciation, B. rapa contains relatively small number of NBS-encoding genes compared to A. thaliana, presumably because of deletion of redundant genes related to genome diploidization. Phylogenetic and evolutionary analyses suggest that relatively higher relaxation of selective constraints on the TNL group after the old duplication event resulted in greater accumulation of TNLs than CNLs in both Arabidopsis and Brassica genomes. Recent tandem duplication and ectopic deletion are likely to have played a role in the generation of novel Brassica lineage-specific resistance genes.     

Genome-Wide Identification and Characterization of Nucleotide-Binding Site (NBS) Resistance Genes in Pineapple

Pineapple is a major tropical fruit and the most important crop processing CAM photosynthesis. It originated in southwest Brazil and northeast Paraguay and survived the harsh, semi-arid environment. Disease resistance genes have contributed to the survival and thriving of this species. The largest class of disease resistance (R) genes in plants consists of genes encoding nucleotide-binding site (NBS) domains. The sequenced genome of pineapple (Ananas comosus (L.) Merr.) provides a resource for analyzing the NBS-encoding genes in this species. A total of 177 NBS-encoding genes were identified using automated and manual analysis criteria, and these represent about 0.6 % of the total number of predicted pineapple genes. Five genes identified here contained the N-terminal Toll/Interleukin-l receptor (TIR) domain, and 46 genes carried the N-terminal Coiled-Coil (CC) motif. A majority of these NBS-encoding genes (84 %) contained a leucine-rich repeat (LRR) domain. A total of 130 of 177 (73 %) of these NBS-encoding genes were distributed across 20 pineapple linkage groups. The identification and characterization of NBS genes in pineapple yielded a valuable genomic resource and improved understanding of R genes in pineapple, which will facilitate the development of disease resistant pineapple cultivars.     

Two TIR:NB:LRR genes are required to specify resistance to Peronospora parasitica isolate Cala2 in Arabidopsis

Resistance responses that plants deploy in defence against pathogens are often triggered following a recognition event mediated by resistance (R) genes. The encoded R proteins usually contain a nucleotide-binding site (NB) and a leucine-rich repeat (LRR) domain. They are further classified into those that contain an N-terminal coiled coil (CC) motif or a Toll interleukin receptor (TIR) domain. Such R genes, when transferred into a susceptible plant of the same or closely related species, usually impart full resistance capability. We have used map-based cloning and mutation analysis to study the recognition of Peronospora parasitica (RPP)2 (At) locus in Arabidopsis accession Columbia (Col-0), which is a determinant of specific recognition of P. parasitica (At) isolate Cala2. Genetic mapping located RPP2 to a 200-kb interval on chromosome 4, which contained four adjacent TIR:NB:LRR genes. Mutational analysis revealed three classes of genes involved in specifying resistance to Cala2. One class, which resulted in pleiotropic effects on resistance to other P. parasitica (At) isolates, was unlinked to the RPP2 locus; this class included AtSGT1b. The other two classes were mapped within the interval and were specific to Cala2 resistance. Representatives of each of these classes were sequenced, and mutations were found in one or the other of two (RPP2A and RPP2B) of the four TIR:NB:LRR genes. RPP2A and RPP2B complemented their specific mutations, but failed to impart resistance when present alone, and it is concluded that both genes are essential determinants for isolate-specific recognition of Cala2. RPP2A has an unusual structure with a short LRR domain at the C-terminus, preceded by two potential but incomplete TIR:NB domains. In addition, the RPP2A LRR domain lacks conserved motifs found in all but three other TIR:NB:LRR class proteins. In contrast, RPP2B has a complete TIR:NB:LRR structure. It is concluded that RPP2A and RPP2B cooperate to specify Cala2 resistance by providing recognition or signalling functions lacked by either partner protein.     

A genome-wide genetic map of NB-LRR disease resistance loci in potato

Like all plants, potato has evolved a surveillance system consisting of a large array of genes encoding for immune receptors that confer resistance to pathogens and pests. The majority of these so-called resistance or R proteins belong to the super-family that harbour a nucleotide binding and a leucine-rich-repeat domain (NB-LRR). Here, sequence information of the conserved NB domain was used to investigate the genome-wide genetic distribution of the NB-LRR resistance gene loci in potato. We analysed the sequences of 288 unique BAC clones selected using filter hybridisation screening of a BAC library of the diploid potato clone RH89-039-16 (S. tuberosum ssp. tuberosum) and a physical map of this BAC library. This resulted in the identification of 738 partial and full-length NB-LRR sequences. Based on homology of these sequences with known resistance genes, 280 and 448 sequences were classified as TIR-NB-LRR (TNL) and CC-NB-LRR (CNL) sequences, respectively. Genetic mapping revealed the presence of 15 TNL and 32 CNL loci. Thirty-six are novel, while three TNL loci and eight CNL loci are syntenic with previously identified functional resistance genes. The genetic map was complemented with 68 universal CAPS markers and 82 disease resistance trait loci described in literature, providing an excellent template for genetic studies and applied research in potato.     

Full-genome analysis of resistance gene homologues in rice

The availability of the rice genome sequence enabled the global characterization of nucleotide-binding site (NBS)–leucine-rich repeat (LRR) genes, the largest class of plant disease resistance genes. The rice genome carries approximately 500 NBS–LRR genes that are very similar to the non-Toll/interleukin-1 receptor homology region (TIR) class (class 2) genes of Arabidopsis but none that are homologous to the TIR class genes. Over 100 of these genes were predicted to be pseudogenes in the rice cultivar Nipponbare, but some of these are functional in other rice lines. Over 80 other NBS-encoding genes were identified that belonged to four different classes, only two of which are present in dicotyledonous plant sequences present in databases. Map positions of the identified genes show that these genes occur in clusters, many of which included members from distantly related groups. Members of phylogenetic subgroups of the class 2 NBS–LRR genes mapped to as many as ten different chromosomes. The patterns of duplication of the NBS–LRR genes indicate that they were duplicated by many independent genetic events that have occurred continuously through the expansion of the NBS–LRR superfamily and the evolution of the modern rice genome. Genetic events, such as inversions, that inhibit the ability of recently duplicated genes to recombine promote the divergence of their sequences by inhibiting concerted evolution.Electronic Supplementary Material Supplementary material is available for this article at     

Genome-wide identification and mapping of NBS-encoding resistance genes in Solanum tuberosum group phureja

The majority of disease resistance (R) genes identified to date in plants encode a nucleotide-binding site (NBS) and leucine-rich repeat (LRR) domain containing protein. Additional domains such as coiled-coil (CC) and TOLL/interleukin-1 receptor (TIR) domains can also be present. In the recently sequenced Solanum tuberosum group phureja genome we used HMM models and manual curation to annotate 435 NBS-encoding R gene homologs and 142 NBS-derived genes that lack the NBS domain. Highly similar homologs for most previously documented Solanaceae R genes were identified. A surprising ～41% (179) of the 435 NBS-encoding genes are pseudogenes primarily caused by premature stop codons or frameshift mutations. Alignment of 81.80% of the 577 homologs to S. tuberosum group phureja pseudomolecules revealed non-random distribution of the R-genes; 362 of 470 genes were found in high density clusters on 11 chromosomes.     

Genome-wide comparative analysis of NBS-encoding genes between Brassica species and Arabidopsis thaliana

Background

Plant disease resistance (R) genes with the nucleotide binding site (NBS) play an important role in offering resistance to pathogens. The availability of complete genome sequences of Brassica oleracea and Brassica rapa provides an important opportunity for researchers to identify and characterize NBS-encoding R genes in Brassica species and to compare with analogues in Arabidopsis thaliana based on a comparative genomics approach. However, little is known about the evolutionary fate of NBS-encoding genes in the Brassica lineage after split from A. thaliana.

Results

Here we present genome-wide analysis of NBS-encoding genes in B. oleracea, B. rapa and A. thaliana. Through the employment of HMM search and manual curation, we identified 157, 206 and 167 NBS-encoding genes in B. oleracea, B. rapa and A. thaliana genomes, respectively. Phylogenetic analysis among 3 species classified NBS-encoding genes into 6 subgroups. Tandem duplication and whole genome triplication (WGT) analyses revealed that after WGT of the Brassica ancestor, NBS-encoding homologous gene pairs on triplicated regions in Brassica ancestor were deleted or lost quickly, but NBS-encoding genes in Brassica species experienced species-specific gene amplification by tandem duplication after divergence of B. rapa and B. oleracea. Expression profiling of NBS-encoding orthologous gene pairs indicated the differential expression pattern of retained orthologous gene copies in B. oleracea and B. rapa. Furthermore, evolutionary analysis of CNL type NBS-encoding orthologous gene pairs among 3 species suggested that orthologous genes in B. rapa species have undergone stronger negative selection than those in B .oleracea species. But for TNL type, there are no significant differences in the orthologous gene pairs between the two species.

Conclusion

This study is first identification and characterization of NBS-encoding genes in B. rapa and B. oleracea based on whole genome sequences. Through tandem duplication and whole genome triplication analysis in B. oleracea, B. rapa and A. thaliana genomes, our study provides insight into the evolutionary history of NBS-encoding genes after divergence of A. thaliana and the Brassica lineage. These results together with expression pattern analysis of NBS-encoding orthologous genes provide useful resource for functional characterization of these genes and genetic improvement of relevant crops.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-3) contains supplementary material, which is available to authorized users.     

Dramatic Number Variation of R Genes in Solanaceae Species Accounted for by a Few R Gene Subfamilies

Most disease resistance genes encode nucleotide-binding-site (NBS) and leucine-rich-repeat (LRR) domains, and the NBS-LRR encoding genes are often referred to as R genes. Using newly developed approach, 478, 485, 1,194, 1,665, 2,042 and 374 R genes were identified from the genomes of tomato Heinz1706, wild tomato LA716, potato DM1-3, pepper Zunla-1 and wild pepper Chiltepin and tobacco TN90, respectively. The majority of R genes from Solanaceae were grouped into 87 subfamilies, including 16 TIR-NBS-LRR (TNL) and 71 non-TNL subfamilies. Each subfamily was annotated manually, including identification of intron/exon structure and intron phase. Interestingly, TNL subfamilies have similar intron phase patterns, while the non-TNL subfamilies have diverse intron phase due to frequent gain of introns. Prevalent presence/absence polymorphic R gene loci were found among Solanaceae species, and an integrated map with 427 R loci was constructed. The pepper genome (2,042 in Chiltepin) has at least four times of R genes as in tomato (478 in Heinz1706). The high number of R genes in pepper genome is due to the amplification of R genes in a few subfamilies, such as the Rpi-blb2 and BS2 subfamilies. The mechanism underlying the variation of R gene number among different plant genomes is discussed.     

Systematic Analysis and Comparison of Nucleotide-Binding Site Disease Resistance Genes in a Diploid Cotton Gossypium raimondii

Plant disease resistance genes are a key component of defending plants from a range of pathogens. The majority of these resistance genes belong to the super-family that harbors a Nucleotide-binding site (NBS). A number of studies have focused on NBS-encoding genes in disease resistant breeding programs for diverse plants. However, little information has been reported with an emphasis on systematic analysis and comparison of NBS-encoding genes in cotton. To fill this gap of knowledge, in this study, we identified and investigated the NBS-encoding resistance genes in cotton using the whole genome sequence information of Gossypium raimondii. Totally, 355 NBS-encoding resistance genes were identified. Analyses of the conserved motifs and structural diversity showed that the most two distinct features for these genes are the high proportion of non-regular NBS genes and the high diversity of N-termini domains. Analyses of the physical locations and duplications of NBS-encoding genes showed that gene duplication of disease resistance genes could play an important role in cotton by leading to an increase in the functional diversity of the cotton NBS-encoding genes. Analyses of phylogenetic comparisons indicated that, in cotton, the NBS-encoding genes with TIR domain not only have their own evolution pattern different from those of genes without TIR domain, but also have their own species-specific pattern that differs from those of TIR genes in other plants. Analyses of the correlation between disease resistance QTL and NBS-encoding resistance genes showed that there could be more than half of the disease resistance QTL associated to the NBS-encoding genes in cotton, which agrees with previous studies establishing that more than half of plant resistance genes are NBS-encoding genes.     

Six amino acid changes confined to the leucine-rich repeat beta-strand/beta-turn motif determine the difference between the P and P2 rust resistance specificities in flax

At least six rust resistance specificities (P and P1 to P5) map to the complex P locus in flax. The P2 resistance gene was identified by transposon tagging and transgenic expression. P2 is a member of a small multigene family and encodes a protein with nucleotide binding site (NBS) and leucine-rich repeat (LRR) domains and an N-terminal Toll/interleukin-1 receptor (TIR) homology domain, as well as a C-terminal non-LRR (CNL) domain of approximately 150 amino acids. A related CNL domain was detected in almost half of the predicted Arabidopsis TIR-NBS-LRR sequences, including the RPS4 and RPP1 resistance proteins, and in the tobacco N protein, but not in the flax L and M proteins. Presence or absence of this domain defines two subclasses of TIR-NBS-LRR resistance genes. Truncations of the P2 CNL domain cause loss of function, and evidence for diversifying selection was detected in this domain, suggesting a possible role in specificity determination. A spontaneous rust-susceptible mutant of P2 contained a G-->E amino acid substitution in the GLPL motif, which is conserved in the NBS domains of plant resistance proteins and the animal cell death control proteins APAF-1 and CED4, providing direct evidence for the importance of this motif in resistance gene function. A P2 homologous gene isolated from a flax line expressing the P resistance specificity encodes a protein with only 10 amino acid differences from the P2 protein. Chimeric gene constructs indicate that just six of these amino acid changes, all located within the predicted beta-strand/beta-turn motif of four LRR units, are sufficient to alter P2 to the P specificity.     

Plant NBS-LRR proteins: adaptable guards

The majority of disease resistance genes in plants encode nucleotide-binding site leucine-rich repeat (NBS-LRR) proteins. This large family is encoded by hundreds of diverse genes per genome and can be subdivided into the functionally distinct TIR-domain-containing (TNL) and CC-domain-containing (CNL) subfamilies. Their precise role in recognition is unknown; however, they are thought to monitor the status of plant proteins that are targeted by pathogen effectors.     

Uncovering the dynamic evolution of nucleotide‐binding site‐leucine‐rich repeat (NBS‐LRR) genes in Brassicaceae

Plant genomes harbor dozens to hundreds of nucleotide-binding site-leucine-rich repeat(NBS-LRR) genes;however,the long-term evolutionary history of these resistance genes has not been fully understood. This study focuses on five Brassicaceae genomes and the Carica papaya genome to explore changes in NBS-LRR genes that have taken place in this Rosid II lineage during the past 72 million years. Various numbers of NBS-LRR genes were identified from Arabidopsis lyrata(198),A. thaliana(165),Brassica rapa(204),Capsella rubella(127),Thellungiella salsuginea(88),and C. papaya(51). In each genome,the identified NBS-LRR genes were found to be unevenly distributed among chromosomes and most of them were clustered together.Phylogenetic analysis revealed that,before and after Brassicaceae speciation events,both toll/interleukin-1receptor-NBS-LRR(TNL) genes and non-toll/interleukin-1receptor-NBS-LRR(n TNL) genes exhibited a pattern of first expansion and then contraction,suggesting that both subclasses of NBS-LRR genes were responding to pathogen pressures synchronically. Further,by examining the gain/loss of TNL and n TNL genes at different evolutionary nodes,this study revealed that both events often occurred more drastically in TNL genes. Finally,the phylogeny of n TNL genes suggested that this NBS-LRR subclass is composed o two separate ancient gene types: RPW8-NBS-LRR and Coiled-coil-NBS-LRR.     

Genome-Wide Analysis in Wild and Cultivated <Emphasis Type="Italic">Oryza</Emphasis> Species Reveals Abundance of NBS Genes in Progenitors of Cultivated Rice

NBS-encoding genes play a critical role in the plant defense system. Wild relatives of crop plants are rich reservoirs of plant defense genes. Here, we performed a stringent genome-wide identification of NBS-encoding genes in three cultivated and eight wild Oryza species, representing three different genomes (AA, BB, and FF) from four continents. A total of 2688 NBS-encoding genes were identified from 11 Oryza genomes. All the three progenitor species of cultivated rice, namely O. barthii, O. rufipogon, and O. nivara, were the richest reservoir of NBS-encoding genes (214, 313, and 307 respectively). Interestingly, the two Asian cultivated species showed a contrasting pattern in the number of NBS-encoding genes. While indica subspecies maintained nearly equal number of NBS genes as its progenitor (309 and 313), the japonica subspecies had retained only two third in the course of evolution (213 and 307). Other major sources for NBS-encoding genes could be (i) O. longistaminata since it had the highest proportion of NBS-encoding genes and (ii) O. glumaepatula as it clustered distinctly away from the rest of the AA genome species. The present study thus revealed that NBS-encoding genes can be exploited from the primary gene pool for disease resistance breeding in rice.     

The Ma gene for complete-spectrum resistance to Meloidogyne species in Prunus is a TNL with a huge repeated C-terminal post-LRR region

Root-knot nematode (RKN) Meloidogyne species are major polyphagous pests of most crops worldwide, and cultivars with durable resistance are urgently needed because of nematicide bans. The Ma gene from the Myrobalan plum (Prunus cerasifera) confers complete-spectrum, heat-stable, and high-level resistance to RKN, which is remarkable in comparison with the Mi-1 gene from tomato (Solanum lycopersicum), the sole RKN resistance gene cloned. We report here the positional cloning and the functional validation of the Ma locus present at the heterozygous state in the P.2175 accession. High-resolution mapping totaling over 3,000 segregants reduced the Ma locus interval to a 32-kb cluster of three Toll/Interleukin1 Receptor-Nucleotide Binding Site-Leucine-Rich Repeat (LRR) genes (TNL1-TNL3), including a pseudogene (TNL2) and a truncated gene (TNL3). The sole complete gene in this interval (TNL1) was validated as Ma, as it conferred the same complete-spectrum and high-level resistance (as in P.2175) using its genomic sequence and native promoter region in Agrobacterium rhizogenes-transformed hairy roots and composite plants. The full-length cDNA (2,048 amino acids) of Ma is the longest of all Resistance genes cloned to date. Its TNL structure is completed by a huge post-LRR (PL) sequence (1,088 amino acids) comprising five repeated carboxyl-terminal PL exons with two conserved motifs. The amino-terminal region (213 amino acids) of the LRR exon is conserved between alleles and contrasts with the high interallelic polymorphisms of its distal region (111 amino acids) and of PL domains. The Ma gene highlights the importance of these uncharacterized PL domains, which may be involved in pathogen recognition through the decoy hypothesis or in nuclear signaling.     

三裂叶薯NBS-LRR类抗病基因的筛选鉴定与结构分析

为发掘甘薯近缘野生种三裂叶薯(Ipomoea triloba)的NBS-LRR类抗病基因,从基因数据库中对三裂叶薯基因组序列进行了筛选、鉴定和分析。结果表明,从三裂叶薯的98 025个基因中,筛选到282个编码NBS-LRR类蛋白的基因,其中N型80个,NL型83个,CN型28个,CNL型57个,TN型10个,TNL型23个,RN型1个。三裂叶薯的16条染色体上均含有NBS-LRR家族基因,数量最多的染色体含有65个,最少的只有1个。三裂叶薯基因组共有55个基因簇,包含了63.5%的NBS-LRR家族基因。在NBS-LRR抗病基因家族中,CNL和TNL亚家族分别对应到7和11个保守结构域。这为三裂叶薯抗性资源的利用提供了科学参考。