- Short Report
- Open Access
TIR-NBS-LRR genes are rare in monocots: evidence from diverse monocot orders
BMC Research Notes volume 2, Article number: 197 (2009)
Plant resistance (R) gene products recognize pathogen effector molecules. Many R genes code for proteins containing nucleotide binding site (NBS) and C-terminal leucine-rich repeat (LRR) domains. NBS-LRR proteins can be divided into two groups, TIR-NBS-LRR and non-TIR-NBS-LRR, based on the structure of the N-terminal domain. Although both classes are clearly present in gymnosperms and eudicots, only non-TIR sequences have been found consistently in monocots. Since most studies in monocots have been limited to agriculturally important grasses, it is difficult to draw conclusions. The purpose of our study was to look for evidence of these sequences in additional monocot orders.
Using degenerate PCR, we amplified NBS sequences from four monocot species (C. blanda, D. marginata, S. trifasciata, and Spathiphyllum sp.), a gymnosperm (C. revoluta) and a eudicot (C. canephora). We successfully amplified TIR-NBS-LRR sequences from dicot and gymnosperm DNA, but not from monocot DNA. Using databases, we obtained NBS sequences from additional monocots, magnoliids and basal angiosperms. TIR-type sequences were not present in monocot or magnoliid sequences, but were present in the basal angiosperms. Phylogenetic analysis supported a single TIR clade and multiple non-TIR clades.
We were unable to find monocot TIR-NBS-LRR sequences by PCR amplification or database searches. In contrast to previous studies, our results represent five monocot orders (Poales, Zingiberales, Arecales, Asparagales, and Alismatales). Our results establish the presence of TIR-NBS-LRR sequences in basal angiosperms and suggest that although these sequences were present in early land plants, they have been reduced significantly in monocots and magnoliids.
Plants recognize pathogens using both non-specific and specific mechanisms. Pattern recognition receptors (PRRs) mediate non-specific recognition by interacting with microbe- or pathogen-associated molecular patterns (MAMPs or PAMPs), while the products of plant resistance (R) genes recognize specific pathogen molecules [1, 2]. Disease resistance is the only known function for R genes, which appear to have a gene-for-gene relationship with pathogen avirulence (avr) genes .
Many R genes code for proteins containing nucleotide binding site (NBS) and C-terminal leucine-rich repeat (LRR) domains. The NBS domain of plant R genes (also called the NB-ARC domain) shares homology with human APAF-1 and C. elegans CED-4, proteins involved in regulating cell death . NBS-LRR proteins can be divided into two groups, TIR-NBS-LRR and non-TIR-NBS-LRR, based on the structure of the N-terminal domain (Figure 1) [5, 6].
The NBS domain from R genes is relatively conserved and contains type-specific motifs (Table 1). The final residue of the kinase-2 motif is especially useful for classifying a sequence as TIR or non-TIR . TIR-type NBS sequences are relatively homogeneous and form a single clade, while non-TIR sequences form multiple clades that likely originated before the split between angiosperms and gymnosperms [8, 9].
The TIR class is found in bryophytes , and both TIR and non-TIR sequences are found in gymnosperms [11, 12]. While both classes are present in eudicots, studies in monocots have reported only non-TIR sequences [7–9, 13]. It is thought that TIR-NBS-LRRs either never developed in monocots  or have been lost [7–9, 13]. However, four resistance gene analogs (RGAs) from the Triticum-Thinopyrum alien addition line TAi-27 have a kinase-2 motif consistent with TIR-NBS-LRR sequences .
Studies of NBS-LRR sequences in monocots have been limited to agriculturally important species in the grass family (Poaceae). Recent studies from Zingiber and Musa species (order Zingiberales) reported only non-TIR type sequences [15–18]. Since there are ten orders of monocots , we are limited in our ability to make generalizations based on information from only two orders. To further investigate the presence of TIR-NBS-LRR sequences in monocots, we combined PCR and bioinformatics to obtain data from additional monocots as well as magnoliids and basal angiosperms (Figure 2).
We amplified sequences from four monocot species representing three monocot orders (Figure 2): Draceana marginata and Sansevieria trifasciata (Asparagales), Spathiphyllum sp. (Alismatales), and Carex blanda (Poales). For comparison, we included a gymnosperm (Cycas revoluta) and a dicot (Coffea canephora). We obtained sequences from a total of 60 PCR products that resulted in 24 unique NBS sequences (Table 2). We found non-TIR type sequences in all plants tested except the cycad, but only two unique TIR-type NBS sequences, one each from C. revoluta and C. canephora.
Using Pfam  and GenBank , we retrieved 17 monocot sequences (ten from Musa acuminata, four from Elaeis guineensis, and three from Zingiber species), all of which we classified as non-TIR-NBS-LRR sequences based on the kinase-2 motif. Although we did not find any new TIR-type sequences from monocots, the search confirmed the similarity of the Triticum-Thinopyrum sequences .
In addition to monocot sequences, we retrieved two sequences from Persea americana (magnoliid) and seven sequences from basal angiosperms (five from Nuphar advena and two from Amborella trichopoda). Based on the kinase-2 motif, both P. americana sequences were non-TIR and all five N. advena sequences were TIR-type sequences. The A. trichopoda sequences have a glutamic acid in the diagnostic position, but downstream motifs similar to TIR-type NBS sequences (data not shown). For comparison, we also retrieved 37 Pinus (gymnosperm) and six Physcomitrella patens (bryophyte) sequences.
We eliminated redundant sequences within a species (>70% identity), resulting in an analysis of 53 plant sequences (Table 3 and Additional Files 1, 2 &3) and human APAF-1 as an outgroup sequence. As much of the NBS domain as was available for each sequence was used for phylogenetic analysis using parsimony criteria (Figure 3). Based on the scaffold tree (see Methods), we compared our clades with those previously reported . All fifteen sequences that we identified as TIR-type NBS sequences based on consensus motifs formed a single clade that was well-supported by bootstrap analysis (85%). The non-TIR-type NBS sequences formed several clades, but many were not well-supported. The well-supported non-TIR clades correspond to non-TIR clades 3 and 4 in Cannon's analysis , while non-TIR clades 1 and 2 are more ambiguous.
Previous studies of plant NBS-LRR sequences have suggested that only non-TIR-NBS-LRR sequences are present in monocots [7–9, 13], and the sequences from the Triticum-Thinopyrum addition line  have not been mentioned in later studies [17, 18, 22, 23]. A study in Agrostis species  reported two TIR-NBS-LRR sequences (Genbank: EE284250, EE284257). However, when these cDNA sequences are translated, they do not contain an open reading frame consistent with a NBS domain (data not shown), so it is unclear that these represent monocot TIR-NBS-LRR genes.
Our PCR strategy amplified TIR-NBS-LRR sequences from dicot and gymnosperm DNA, but not from monocot DNA
In spite of attempts to bias amplification and cloning toward TIR-type NBS sequences, we did not find TIR-NBS-LRR sequences in any of the monocots we tested, although we easily cloned and sequenced TIR-type sequences from a gymnosperm, eudicot, and Arabidopsis control reactions. Although we expected to find the TIR class in dicots, a previous study in coffee did not report any TIR-type sequences .
Our results support the hypothesis that TIR-type NBS sequences are rare in monocots. We used diverse monocot taxa, including a species closely related to grasses (C. blanda) and a species from a basal monocot order (Spathiphyllum sp.). As with any PCR study, we cannot eliminate the possibility that TIR-NBS-LRR sequences in monocots are too divergent for our primers to amplify. Species-specific amplification has been reported  and more comparative work is needed to confirm that there are definitive consensus sequences for these motifs that are well-conserved across diverse taxa.
Database searches and phylogenetic analysis show that TIR-type NBS sequences are present in basal angiosperms, but are rare in monocots and magnoliids
Using Pfam and GenBank, we obtained NBS sequences from monocots, a magnoliid, basal angiosperms, gymnosperms, and a bryophyte. The pine and moss sequences had been classified previously [10–12] and provided diverse lineages for comparison to the predominantly monocot and dicot sequences in our study. Based on the kinase-2 motif, TIR-type sequences were absent from monocots and magnoliids (with the exception of the reported Triticum-Thinopyrum sequences), but were present in basal angiosperms, gymnosperms, and bryophytes (Figure 2).
Our phylogenetic analysis (Figure 3) was consistent with previous analyses that showed a single TIR clade and multiple non-TIR clades [8, 9]. The full NBS domain was not available for some sequences used in the analysis. We expect that the phylogenetic relationships will be clarified as more sequence becomes available. Our clear non-TIR clades corresponded to N3 and N4 (Cannon), with N1 and N2 split into several poorly-supported clades. Cannon's N4 clade did not include monocots, while our analysis placed a Z. officinale sequence in this clade (Figure 3). Cannon reported that N1.2 might be monocot specific, but our corresponding clade was not well-supported. Based on the current analysis, N3.2 may be monocot specific.
We expected to find both dicots and gymnosperms represented across TIR and non-TIR clades, but gymnosperm sequences were only found in the TIR and N4 clades. Both magnoliid sequences were non-TIR type sequences, and all basal angiosperm sequences were in the TIR clade (Figure 3). As more basal angiosperm sequences become available, we expect to also find non-TIR-type NBS sequences. Our results suggest that TIR-type sequences are rare in both magnoliids and monocots.
TIR-NBS-LRR vs. the TIR domain
The rarity of TIR-NBS-LRR sequences in monocots does not necessarily reflect on the abundance of the TIR domain itself. Two similar protein families that may act as adapters (TIR-NBS and TIR-X) are found in monocots, but in low numbers compared to dicots and gymnosperms . Many sequences in the databases are fragments, and some predicted non-TIR-NBS-LRR proteins may be members of these new families .
Classification of NBS-LRR proteins into TIR and non-TIR (Figure 1) is based on consensus motifs within the NBS domain (Table 1). Although these are assumed to be diagnostic , the presence or absence of a TIR domain has not usually been confirmed. We cannot eliminate the possibility that these motifs are not diagnostic. Further sequencing of the N-terminal regions of these genes is needed to confirm that our categorization is correct.
We were unable to find monocot TIR-NBS-LRR sequences by PCR amplification or database searches. In contrast to previous studies, our results represent five monocot orders (Poales, Zingiberales, Arecales, Asparagales, and Alismatales) as well as basal angiosperms and magnoliids. Establishing the presence of TIR-type NBS sequences in basal angiosperms fills a gap in our knowledge of these important genes. Our results suggest that although TIR-type NBS-LRR sequences were present in early land plants, they have been reduced significantly in monocots. The sequences from the Triticum-Thinopyrum line  remain the only reported monocot TIR-NBS-LRR sequences. We do not know when these sequences were lost, but the P. americana sequences suggest that TIR-NBS-LRR sequences are rare in magnoliids as well. It is not clear whether these sequences were lost independently in both lineages or prior to their divergence. Further sequencing from additional taxa and confirmation that the motifs in the NBS region are diagnostic will be helpful in clarifying the evolutionary history of plant R genes.
We vegetatively propagated Carex blanda individuals collected from Kansas  and grew Arabidopsis thaliana from seeds (Ruth Shaw, University of Minnesota) in the greenhouse. Draceana marginata, Sansevieria trifasciata, Spathiphyllum sp., Cycas revoluta, and Coffea canephora came from the University of Kansas greenhouse.
Genomic DNA isolation
We isolated genomic DNA from frozen plant material by grinding with a mortar and pestle in extraction buffer (0.2 M Tris-HCl pH 7.5, 0.25 M NaCl, 25 mM EDTA pH 8.0, 0.5% SDS) and incubating at 65°C for 10 minutes. We performed three extractions with phenol:chloroform:isoamyl alcohol (25:24:1); the final extraction used a phase-lock tube (Eppendorf). We precipitated the DNA and resuspended it in 100 μL10 mM Tris pH 8.0 with 10 μg RNaseA.
We amplified 200 ng of genomic DNA with degenerate primers (Invitrogen). We used two 5' primers (P-loop/Kinase-1a region) and twelve 3' primers (GLPL and RNBS-D regions; Figure 4 and Table 4). The 50 μL reactions contained 2.5 units Platinum Taq polymerase (Invitrogen), 5 μL 10× buffer, 1.5 μL 50 mM MgCl2, 1 μL 10 mM dNTP mix, and 5 μL each degenerate primer (10 μM). The cycling parameters included initial denaturation at 94°C for 2 minutes, 35 cycles of 94°C for 30 seconds, 50°C for 30 seconds, 72°C for 1 minute, and a final extension at 72°C for 10 minutes.
Cloning and sequencing
We purified PCR products from agarose gels using the QIAquick gel extraction kit (QIAgen), ligated into the pCR4-TOPO vector (Invitrogen), and transformed into maximum efficiency DH5α competent cells. We isolated plasmid DNA from overnight cultures using a standard alkaline lysis protocol, digested with EcoRI, and sequenced representative clones with standard primers (T3/T7) at the University of Kansas sequencing facility or ACGT, Inc. We typically sequenced between one and three clones from each fragment. Sequences were submitted to GenBank with accession numbers EF687860-EF687864 and EF687876-EF687894.
We retrieved monocot (excluding Poales), gymnosperm, and bryophyte sequences from Pfam by viewing the NB-ARC domain (PF00931) species distribution. We retrieved sequences from GenBank by searching with known plant NBS sequences [see Additional File 2] and using the taxonomy reports to identify non-dicot, non-grass plant sequences with significant similarity (e < 0.05). We used EditSeq (Lasergene) for all sequence viewing and editing, and excluded sequences with internal stops or that lacked the diagnostic kinase-2 region.
We aligned sequences using the ClustalW algorithm in MegAlign (Lasergene) with manual adjustments as necessary. The alignment contained a core region of approximately 170 amino acids. We generated a scaffold tree based on the Pfam seed alignment and previous analyses [see Additional File 3] . Phylogenetic analysis with parsimony criteria was performed using PAUP* 4.0 beta10 .
Chisholm ST, Coaker G, Day B, Staskawicz BJ: Host-microbe interactions: shaping the evolution of the plant immune response. Cell. 2006, 124: 803-814. 10.1016/j.cell.2006.02.008.
Jones JDG, Dangl JL: The plant immune system. Nature. 2006, 444: 323-329. 10.1038/nature05286.
Hammond-Kosack KE, Jones JDG: Plant disease resistance genes. Annu Rev Plant Physiol Plant Mol Biol. 1997, 48: 575-607. 10.1146/annurev.arplant.48.1.575.
Biezen van der EA, Jones JDG: The NB-ARC domain: a novel signalling motif shared by plant resistance gene products and regulators of cell death in animals. Curr Biol. 1998, 8: R226-R227. 10.1016/S0960-9822(98)70145-9.
Dangl JL, Jones JDG: Plant pathogens and integrated defence responses to infection. Nature. 2001, 411: 826-833. 10.1038/35081161.
Hulbert SH, Webb CA, Smith SM, Sun Q: Resistance gene complexes: evolution and utilization. Annu Rev Phytopathol. 2001, 39: 285-312. 10.1146/annurev.phyto.39.1.285.
Meyers BC, Dickerman AW, Michelmore RW, Sivaramakrishnan S, Sobral BW, Young ND: Plant disease resistance genes encode members of an ancient and diverse protein family within the nucleotide-binding superfamily. Plant J. 1999, 20: 317-332. 10.1046/j.1365-313X.1999.t01-1-00606.x.
Cannon SB, Zhu H, Baumgarten AM, Spangler R, May G, Cook DR, Young ND: Diversity, distribution, and ancient taxonomic relationships within the TIR and non-TIR NBS-LRR resistance gene subfamilies. J Mol Evol. 2002, 54: 548-562. 10.1007/s00239-001-0057-2.
Bai J, Pennill LA, Ning J, Lee SW, Ramalingam J, Webb CA, Zhao B, Sun Q, Nelson JC, Leach JE, Hulbert SH: Diversity in nucleotide binding site--leucine-rich repeat genes in cereals. Genome Res. 2002, 12: 1871-1884. 10.1101/gr.454902.
Akita M, Valkonen JP: A novel gene family in moss (Physcomitrella patens) shows sequence homology and a phylogenetic relationship with the TIR-NBS class of plant disease resistance genes. J Mol Evol. 2002, 55 (5): 595-605. 10.1007/s00239-002-2355-8.
Liu JJ, Ekramoddoullah AK: Isolation, genetic variation and expression of TIR-NBS-LRR resistance gene analogs from western white pine (Pinus monticola Dougl. ex. D. Don.). Mol Genet Genomics. 2003, 270 (5): 432-441. 10.1007/s00438-003-0940-1.
Jermstad K, Sheppard L, Kinloch B, Delfino-Mix A, Ersoz E, Krutovsky K, Neale D: Isolation of a full-length CC-NBS-LRR resistance gene analog candidate from sugar pine showing low nucleotide diversity. Tree Genet Genomes. 2006, 2: 76-85. 10.1007/s11295-005-0029-6.
Pan Q, Wendel J, Fluhr R: Divergent evolution of plant NBS-LRR resistance gene homologues in dicot and cereal genomes. J Mol Evol. 2000, 50: 203-213.
Jiang S-M, Hu J, Yin W-B, Chen Y-H, Wang RR-C, Hu Z-M: Cloning of resistance gene analogs located on the alien chromosome in an addition line of wheat-Thinopyrum intermedium. Theor Appl Genet. 2005, 111: 923-931. 10.1007/s00122-005-0022-3.
Aswati Nair R, Thomas G: Isolation, characterization and expression studies of resistance gene candidates (RGCs) from Zingiber spp. Theor Appl Genet. 2007, 116: 123-134. 10.1007/s00122-007-0652-8.
Azhar M, Heslop-Harrison J: Genomes, diversity and resistance gene analogues in Musa species. Cytogenet Genome Res. 2008, 121: 59-66. 10.1159/000124383.
Miller RN, Bertioli DJ, Baurens FC, Santos CM, Alves PC, Martins NF, Togawa RC, Souza MT, Pappas GJ: Analysis of non-TIR NBS-LRR resistance gene analogs is Musa acuminata Colla: Isolation, RFLP marker development, and physical mapping. BMC Plant Biol. 2008, 8: 15-10.1186/1471-2229-8-15.
Pei X, Li S, Jiang Y, Zhang Y, Wang Z, Jia S: Isolation, characterization, and phylogenetic analysis of the resistance gene analogues (RGAs) in banana (Musa spp.). Plant Sci. 2007, 172: 1166-1174. 10.1016/j.plantsci.2007.02.019.
Angiosperm Phylogeny Group: An update of the Angiosperm Phylogeny Group classification for the orders and families of flowering plants: APG II. Bot J Linn Soc. 2003, 141: 399-436. 10.1046/j.1095-8339.2003.t01-1-00158.x.
Finn RD, Mistry J, Schuster-Bockler B, Griffiths-Jones S, Hollich V, Lassmann T, Moxon S, Marshall M, Khanna A, Durbin R, et al: Pfam: clans, web tools and services. Nucleic Acids Res. 2006, 34: D247-D251. 10.1093/nar/gkj149.
Benson DA, Karsch-Mizrachi I, Lipman DJ, Ostell J, Wheeler DL: GenBank. Nucleic Acids Res. 2007, 35: D21-D25. 10.1093/nar/gkl986.
Budak H, Su S, Ergen N: Revealing constitutively expressed resistance genes in Agrostis species using PCR-based motif-directed RNA fingerprinting. Genet Res. 2006, 88: 165-175. 10.1017/S0016672307008518.
Ameline-Torregrosa C, Wang B-B, O'Bleness MS, Deshpande S, Zhu H, Roe B, Young ND, Cannon SB: Identification and characterization of nucleotide-binding-site-leucine-rich repeat genes in the model plant Medicago truncatula. Plant Physiol. 2008, 146: 5-21. 10.1104/pp.107.104588.
Noir S, Combes MC, Anthony F, Lashermes P: Origin, diversity and evolution of NBS-type disease-resistance gene homologues in coffee trees (Coffea L.). Mol Genet Genomics. 2001, 265: 654-662. 10.1007/s004380100459.
Meyers BC, Morgante M, Michelmore RW: TIR-X and TIR-NBS proteins:two new families related to disease resistance TIR-NBS-LRR proteins encoded in Arabidopsis and other plant genomes. Plant J. 2002, 32: 77-92. 10.1046/j.1365-313X.2002.01404.x.
Alexander H, Price S, Houser R, Finch D, Tourtellot M: Is there a reduction in disease and pre-dispersal predation at the border of a host plant's range? Field and herbarium studies of Carex blanda. J Ecol. 2007, 95: 446-457. 10.1111/j.1365-2745.2007.01228.x.
Swofford D: PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). Version 4. 1998, Sinauer Associates, Sunderland, Massachusetts
Moffett P, Farnham G, Peart J, Baulcombe DC: Interaction between domains of a plant NBS-LRR protein in disease resistance-related cell death. EMBO J. 2002, 21: 4511-4519. 10.1093/emboj/cdf453.
Aarts N, Metz M, Holub E, Staskawicz BJ, Daniels MJ, Parker JE: Different requirements for EDS1 and NDR1 by disease resistance genes define at least two R gene-mediated signaling pathways in Arabidopsis. Proc Natl Acad Sci USA. 1998, 95: 10306-10311. 10.1073/pnas.95.17.10306.
DeYoung BJ, Innes RW: Plant NBS-LRR proteins in pathogen sensing and host defense. Nature Immunol. 2006, 7: 1243-1249. 10.1038/ni1410.
We thank T. Hildebrand for encouraging the work in diverse monocot species, P. Cartwright and C. Haufler for use of lab equipment, K. Nus for greenhouse maintenance, and H. Yu at ACGT, Inc. for special pricing of sequencing services. We thank B. Williams for technical support. L. Hileman and M. Mort provided helpful comments on the manuscript.
DEKT was supported by the KU-Haskell IRACDA Project (NIH:5K12GM063651-02).
The authors declare that they have no competing interests.
DT participated in sample collection, carried out the PCR study, performed the database searches and phylogenetic analysis, and drafted the manuscript. HA participated in sample collection, study design, and manuscript revision. Both authors read and approved the final manuscript.