Diversification of non-TIR class NB-LRR genes in relation to whole-genome duplication events in Arabidopsis

An NLR paralog Pit2 generated from tandem duplication of Pit1 fine-tunes Pit1 localization and function

NLR family proteins act as intracellular receptors. Gene duplication amplifies the number of NLR genes, and subsequent mutations occasionally provide modifications to the second gene that benefits immunity. However, evolutionary processes after gene duplication and functional relationships between duplicated NLRs remain largely unclear. Here, we report that the rice NLR protein Pit1 is associated with its paralogue Pit2. The two are required for the resistance to rice blast fungus but have different functions: Pit1 induces cell death, while Pit2 competitively suppresses Pit1-mediated cell death. During evolution, the suppression of Pit1 by Pit2 was probably generated through positive selection on two fate-determining residues in the NB-ARC domain of Pit2, which account for functional differences between Pit1 and Pit2. Consequently, Pit2 lost its plasma membrane localization but acquired a new function to interfere with Pit1 in the cytosol. These findings illuminate the evolutionary trajectory of tandemly duplicated NLR genes after gene duplication.

Genome-Wide Identification and Evolution of Receptor-Like Kinases (RLKs) and Receptor like Proteins (RLPs) in Brassica juncea

Brassica juncea, an allotetraploid species, is an important germplasm resource for canola improvement, due to its many beneficial agronomic traits, such as heat and drought tolerance and blackleg resistance. Receptor-like kinase (RLK) and receptor-like protein (RLP) genes are two types of resistance gene analogues (RGA) that play important roles in plant innate immunity, stress response and various development processes. In this study, genome wide analysis of RLKs and RLPs is performed in B. juncea. In total, 493 RLKs (LysM-RLKs and LRR-RLKs) and 228 RLPs (LysM-RLPs and LRR-RLPs) are identified in the genome of B. juncea, using RGAugury. Only 13.54% RLKs and 11.79% RLPs are observed to be grouped within gene clusters. The majority of RLKs (90.17%) and RLPs (52.83%) are identified as duplicates, indicating that gene duplications significantly contribute to the expansion of RLK and RLP families. Comparative analysis between B. juncea and its progenitor species, B. rapa and B. nigra, indicate that 83.62% RLKs and 41.98% RLPs are conserved in B. juncea, and RLPs are likely to have a faster evolution than RLKs. This study provides a valuable resource for the identification and characterisation of candidate RLK and RLP genes.

Genome-wide identification, characterization, and evolutionary analysis of NBS-encoding resistance genes in barley

In this study, a systematic analysis of Nucleotide-Binding Site (NBS) disease resistance (R) gene family in the barley, Hordeum vulgare L. cv. Bowman, genome was performed. Using multiple computational analyses, we could identify 96 regular NBS-encoding genes and characterize them on the bases of structural diversity, conserved protein signatures, genomic distribution, gene duplications, differential expression, selection pressure, codon usage, regulation by microRNAs and phylogenetic relationships. Depending on the presence or absence of CC and LRR domains; the identified NBS genes were assigned to four distinct groups; NBS-LRR (53.1%), CC-NBS-LRR (14.6%), NBS (26%), and CC-NBS (6.3%). NBS-associated domain analysis revealed the presence of signal peptides, zinc fingers, diverse kinases, and other structural features. Eighty-five of the identified NBS-encoding genes were mapped onto the seven barley chromosomes, revealing that 50% of them were located on chromosomes 7H, 2H, and 3H, with a tendency of NBS genes to be clustered in the distal telomeric regions of the barley chromosomes. Nine gene clusters, representing 22.35% of total mapped barley NBS-encoding genes, were found, suggesting that tandem duplication stands for an important mechanism in the expansion of this gene family in barley. Phylogenetic analysis determined 31 HvNBS orthologs from rice and Brachypodium. 87 out of 96 HvNBSs were supported by expression evidence, exhibiting various and quantitatively uneven expression patterns across distinct tissues, organs, and development stages. Fourteen potential miRNA-R gene target pairs were further identified, providing insight into the regulation of NBS genes expression. These findings offer candidate target genes to engineer disease-resistant barley genotypes, and promote our understanding of the evolution of NBS-encoding genes in Poaceae crops.

A unique RPW8-encoding class of genes that originated in early land plants and evolved through domain fission, fusion, and duplication

Duplication, lateral gene transfer, domain fusion/fission and de novo domain creation play a key role in formation of initial common ancestral protein. Abundant protein diversities are produced by domain rearrangements, including fusions, fissions, duplications, and terminal domain losses. In this report, we explored the origin of the RPW8 domain and examined the domain rearrangements that have driven the evolution of RPW8-encoding genes in land plants. The RPW8 domain first emerged in the early land plant, Physcomitrella patens, and it likely originated de novo from a non-coding sequence or domain divergence after duplication. It was then incorporated into the NBS-LRR protein to create a main sub-class of RPW8-encoding genes, the RPW8-NBS-encoding genes. They evolved by a series of genetic events of domain fissions, fusions, and duplications. Many species-specific duplication events and tandemly duplicated clusters clearly demonstrated that species-specific and tandem duplications played important roles in expansion of RPW8-encoding genes, especially in gymnosperms and species of the Rosaceae. RPW8 domains with greater Ka/Ks values than those of the NBS domains indicated that they evolved faster than the NBS domains in RPW8-NBSs.

Large-Scale Analyses of Angiosperm Nucleotide-Binding Site-Leucine-Rich Repeat Genes Reveal Three Anciently Diverged Classes with Distinct Evolutionary Patterns

Nucleotide-binding site-leucine-rich repeat (NBS-LRR) genes make up the largest plant disease resistance gene family (R genes), with hundreds of copies occurring in individual angiosperm genomes. However, the expansion history of NBS-LRR genes during angiosperm evolution is largely unknown. By identifying more than 6,000 NBS-LRR genes in 22 representative angiosperms and reconstructing their phylogenies, we present a potential framework of NBS-LRR gene evolution in the angiosperm. Three anciently diverged NBS-LRR classes (TNLs, CNLs, and RNLs) were distinguished with unique exon-intron structures and DNA motif sequences. A total of seven ancient TNL, 14 CNL, and two RNL lineages were discovered in the ancestral angiosperm, from which all current NBS-LRR gene repertoires were evolved. A pattern of gradual expansion during the first 100 million years of evolution of the angiosperm clade was observed for CNLs. TNL numbers remained stable during this period but were eventually deleted in three divergent angiosperm lineages. We inferred that an intense expansion of both TNL and CNL genes started from the Cretaceous-Paleogene boundary. Because dramatic environmental changes and an explosion in fungal diversity occurred during this period, the observed expansions of R genes probably reflect convergent adaptive responses of various angiosperm families. An ancient whole-genome duplication event that occurred in an angiosperm ancestor resulted in two RNL lineages, which were conservatively evolved and acted as scaffold proteins for defense signal transduction. Overall, the reconstructed framework of angiosperm NBS-LRR gene evolution in this study may serve as a fundamental reference for better understanding angiosperm NBS-LRR genes.

Genome-wide analysis of the gene families of resistance gene analogues in cotton and their response to Verticillium wilt

BACKGROUND

Gossypium raimondii is a Verticillium wilt-resistant cotton species whose genome encodes numerous disease resistance genes that play important roles in the defence against pathogens. However, the characteristics of resistance gene analogues (RGAs) and Verticillium dahliae response loci (VdRLs) have not been investigated on a global scale. In this study, the characteristics of RGA genes were systematically analysed using bioinformatics-driven methods. Moreover, the potential VdRLs involved in the defence response to Verticillium wilt were identified by RNA-seq and correlations with known resistance QTLs.

RESULTS

The G. raimondii genome encodes 1004 RGA genes, and most of these genes cluster in homology groups based on high levels of similarity. Interestingly, nearly half of the RGA genes occurred in 26 RGA-gene-rich clusters (Rgrcs). The homology analysis showed that sequence exchanges and tandem duplications frequently occurred within Rgrcs, and segmental duplications took place among the different Rgrcs. An RNA-seq analysis showed that the RGA genes play roles in cotton defence responses, forming 26 VdRLs inside in the Rgrcs after being inoculated with V. dahliae. A correlation analysis found that 12 VdRLs were adjacent to the known Verticillium wilt resistance QTLs, and that 5 were rich in NB-ARC domain-containing disease resistance genes.

CONCLUSIONS

The cotton genome contains numerous RGA genes, and nearly half of them are located in clusters, which evolved by sequence exchanges, tandem duplications and segmental duplications. In the Rgrcs, 26 loci were induced by the V. dahliae inoculation, and 12 are in the vicinity of known Verticillium wilt resistance QTLs.

Identification and characterization of rye genes not expressed in allohexaploid triticale

Background

One of the most important evolutionary processes in plants is polyploidization. The combination of two or more genomes in one organism often initially leads to changes in gene expression and extensive genomic reorganization, compared to the parental species. Hexaploid triticale (x Triticosecale) is a synthetic hybrid crop species generated by crosses between T. turgidum and Secale cereale. Because triticale is a recent synthetic polyploid it is an important model for studying genome evolution following polyploidization. Molecular studies have demonstrated that genomic sequence changes, consisting of sequence elimination or loss of expression of genes from the rye genome, are common in triticale. High-throughput DNA sequencing allows a large number of genes to be surveyed, and transcripts from the different homeologous copies of the genes that have high sequence similarity can be better distinguished than hybridization methods previously employed.

Results

The expression levels of 23,503 rye cDNA reference contigs were analyzed in 454-cDNA libraries obtained from anther, root and stem from both triticale and rye, as well as in five 454-cDNA data sets created from triticale seedling shoot, ovary, stigma, pollen and seed tissues to identify the classes of rye genes silenced or absent in the recent synthetic hexaploid triticale. Comparisons between diploid rye and hexaploid triticale detected 112 rye cDNA contigs (~0.5%) that were totally undetected by expression analysis in all triticale tissues, although their expression was relatively high in rye tissues. Non-expressed rye genes were found to be strikingly less similar to their closest BLASTN matches in the wheat genome or in the other Triticum genomes than a test set of 200 random rye genes. Genes that were not detected in the RNA-seq data were further characterized by testing for their presence in the triticale genome by PCR using genomic DNA as a template.

Conclusion

Genes with low similarity between rye sequences and their closest matches in the Triticum genome have a higher probability to be repressed or absent in the allopolyploid genome.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1480-x) contains supplementary material, which is available to authorized users.

Evolution and Conservation of Plant NLR Functions

In plants and animals, nucleotide-binding domain and leucine-rich repeats (NLR)-containing proteins play pivotal roles in innate immunity. Despite their similar biological functions and protein architecture, comparative genome-wide analyses of NLRs and genes encoding NLR-like proteins suggest that plant and animal NLRs have independently arisen in evolution. Furthermore, the demonstration of interfamily transfer of plant NLR functions from their original species to phylogenetically distant species implies evolutionary conservation of the underlying immune principle across plant taxonomy. In this review we discuss plant NLR evolution and summarize recent insights into plant NLR-signaling mechanisms, which might constitute evolutionarily conserved NLR-mediated immune mechanisms.

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.

Overview of tomato (Solanum lycopersicum) candidate pathogen recognition genes reveals important Solanum R locus dynamics

To investigate the genome-wide spatial arrangement of R loci, a complete catalogue of tomato (Solanum lycopersicum) and potato (Solanum tuberosum) nucleotide-binding site (NBS) NBS, receptor-like protein (RLP) and receptor-like kinase (RLK) gene repertories was generated. Candidate pathogen recognition genes were characterized with respect to structural diversity, phylogenetic relationships and chromosomal distribution. NBS genes frequently occur in clusters of related gene copies that also include RLP or RLK genes. This scenario is compatible with the existence of selective pressures optimizing coordinated transcription. A number of duplication events associated with lineage-specific evolution were discovered. These findings suggest that different evolutionary mechanisms shaped pathogen recognition gene cluster architecture to expand and to modulate the defence repertoire. Analysis of pathogen recognition gene clusters associated with documented resistance function allowed the identification of adaptive divergence events and the reconstruction of the evolution history of these loci. Differences in candidate pathogen recognition gene number and organization were found between tomato and potato. Most candidate pathogen recognition gene orthologues were distributed at less than perfectly matching positions, suggesting an ongoing lineage-specific rearrangement. Indeed, a local expansion of Toll/Interleukin-1 receptor (TIR)-NBS-leucine-rich repeat (LRR) (TNL) genes in the potato genome was evident. Taken together, these findings have implications for improved understanding of the mechanisms of molecular adaptive selection at Solanum R loci.

Genome-wide identification and characterization of nucleotide binding site leucine-rich repeat genes in linseed reveal distinct patterns of gene structure

Plants employ different disease-resistance genes to detect pathogens and to induce defense responses. The largest class of these genes encodes proteins with nucleotide binding site (NBS) and leucine-rich repeat (LRR) domains. To identify the putative NBS-LRR encoding genes from linseed, we analyzed the recently published linseed genome sequence and identified 147 NBS-LRR genes. The NBS domain was used for phylogeny construction and these genes were classified into two well-known families, non-TIR (CNL) and TIR related (TNL), and formed eight clades in the neighbor-joining bootstrap tree. Eight different gene structures were observed among these genes. An unusual domain arrangement was observed in the TNL family members, predominantly in the TNL-5 clade members belonging to class D. About 12% of the genes observed were linseed specific. The study indicated that the linseed genes probably have an ancient origin with few progenitor genes. Quantitative expression analysis of five genes showed inducible expression. The in silico expression evidence was obtained for a few of these genes, and the expression was not correlated with the presence of any particular regulatory element or with unusual domain arrangement in those genes. This study will help in understanding the evolution of these genes, the development of disease resistant varieties, and the mechanism of disease resistance in linseed.

Evolution of a complex disease resistance gene cluster in diploid Phaseolus and tetraploid Glycine

We used a comparative genomics approach to investigate the evolution of a complex nucleotide-binding (NB)-leucine-rich repeat (LRR) gene cluster found in soybean (Glycine max) and common bean (Phaseolus vulgaris) that is associated with several disease resistance (R) genes of known function, including Rpg1b (for Resistance to Pseudomonas glycinea1b), an R gene effective against specific races of bacterial blight. Analysis of domains revealed that the amino-terminal coiled-coil (CC) domain, central nucleotide-binding domain (NB-ARC [for APAF1, Resistance genes, and CED4]), and carboxyl-terminal LRR domain have undergone distinct evolutionary paths. Sequence exchanges within the NB-ARC domain were rare. In contrast, interparalogue exchanges involving the CC and LRR domains were common, consistent with both of these regions coevolving with pathogens. Residues under positive selection were overrepresented within the predicted solvent-exposed face of the LRR domain, although several also were detected within the CC and NB-ARC domains. Superimposition of these latter residues onto predicted tertiary structures revealed that the majority are located on the surface, suggestive of a role in interactions with other domains or proteins. Following polyploidy in the Glycine lineage, NB-LRR genes have been preferentially lost from one of the duplicated chromosomes (homeologues found in soybean), and there has been partitioning of NB-LRR clades between the two homeologues. The single orthologous region in common bean contains approximately the same number of paralogues as found in the two soybean homeologues combined. We conclude that while polyploidization in Glycine has not driven a stable increase in family size for NB-LRR genes, it has generated two recombinationally isolated clusters, one of which appears to be in the process of decay.

Relative evolutionary rates of NBS-encoding genes revealed by soybean segmental duplication

It is well known that nucleotide binding site (NBS)-encoding genes are duplicate-rich and fast-evolving genes. However, there is little information on the relative importance of tandem and segmental NBS duplicates and their exact evolutionary rates. The two rounds of large-scale duplication that have occurred in soybean provide a unique opportunity to investigate these issues. Comparison of NBS and non-NBS genes on segments of syntenic homoeologs shows that NBS-encoding genes evolve at least 1.5-fold faster (~1.5-fold higher synonymous and approximately 2.3-fold higher nonsynonymous substitution rates) and lose their genes approximately twofold faster than the flanking non-NBS genes. Compared with segmental duplicates, tandem NBS duplicates are more abundant in soybean, suggesting that tandem duplication is the major driving force in the expansion of NBS genes. Notably, significant sequence exchanges along with significantly positive selection were detected in most tandem-duplicated NBS gene families. The results suggest that the rapid evolution of NBS genes may be due to the combined effects of diversifying selection and frequent sequence exchanges. Interestingly, TIR-NBS-LRR genes (TNLs) have a higher nucleotide substitution rate than non-TNLs, indicating that these types of NBS genes may have a rather different evolutionary pattern. It is important to determine the exact relative evolutionary rates of TNL, non-TNL, and non-NBS genes in order to understand how fast the host plant can adjust its response to rapidly evolving pathogens in a coevolutionary context.

A nomadic subtelomeric disease resistance gene cluster in common bean

The B4 resistance (R) gene cluster is one of the largest clusters known in common bean (Phaseolus vulgaris [Pv]). It is located in a peculiar genomic environment in the subtelomeric region of the short arm of chromosome 4, adjacent to two heterochromatic blocks (knobs). We sequenced 650 kb spanning this locus and annotated 97 genes, 26 of which correspond to Coiled-Coil-Nucleotide-Binding-Site-Leucine-Rich-Repeat (CNL). Conserved microsynteny was observed between the Pv B4 locus and corresponding regions of Medicago truncatula and Lotus japonicus in chromosomes Mt6 and Lj2, respectively. The notable exception was the CNL sequences, which were completely absent in these regions. The origin of the Pv B4-CNL sequences was investigated through phylogenetic analysis, which reveals that, in the Pv genome, paralogous CNL genes are shared among nonhomologous chromosomes (4 and 11). Together, our results suggest that Pv B4-CNL was derived from CNL sequences from another cluster, the Co-2 cluster, through an ectopic recombination event. Integration of the soybean (Glycine max) genome data enables us to date more precisely this event and also to infer that a single CNL moved from the Co-2 to the B4 cluster. Moreover, we identified a new 528-bp satellite repeat, referred to as khipu, specific to the Phaseolus genus, present both between B4-CNL sequences and in the two knobs identified at the B4 R gene cluster. The khipu repeat is present on most chromosomal termini, indicating the existence of frequent ectopic recombination events in Pv subtelomeric regions. Our results highlight the importance of ectopic recombination in R gene evolution.

Genome-wide identification of NBS-encoding resistance genes in Brassica rapa

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.

Differential accumulation of retroelements and diversification of NB-LRR disease resistance genes in duplicated regions following polyploidy in the ancestor of soybean

The genomes of most, if not all, flowering plants have undergone whole genome duplication events during their evolution. The impact of such polyploidy events is poorly understood, as is the fate of most duplicated genes. We sequenced an approximately 1 million-bp region in soybean (Glycine max) centered on the Rpg1-b disease resistance gene and compared this region with a region duplicated 10 to 14 million years ago. These two regions were also compared with homologous regions in several related legume species (a second soybean genotype, Glycine tomentella, Phaseolus vulgaris, and Medicago truncatula), which enabled us to determine how each of the duplicated regions (homoeologues) in soybean has changed following polyploidy. The biggest change was in retroelement content, with homoeologue 2 having expanded to 3-fold the size of homoeologue 1. Despite this accumulation of retroelements, over 77% of the duplicated low-copy genes have been retained in the same order and appear to be functional. This finding contrasts with recent analyses of the maize (Zea mays) genome, in which only about one-third of duplicated genes appear to have been retained over a similar time period. Fluorescent in situ hybridization revealed that the homoeologue 2 region is located very near a centromere. Thus, pericentromeric localization, per se, does not result in a high rate of gene inactivation, despite greatly accelerated retrotransposon accumulation. In contrast to low-copy genes, nucleotide-binding-leucine-rich repeat disease resistance gene clusters have undergone dramatic species/homoeologue-specific duplications and losses, with some evidence for partitioning of subfamilies between homoeologues.

Recent duplications dominate NBS-encoding gene expansion in two woody species

Most disease resistance genes in plants encode NBS-LRR proteins. However, in woody species, little is known about the evolutionary history of these genes. Here, we identified 459 and 330 respective NBS-LRRs in grapevine and poplar genomes. We subsequently investigated protein motif composition, phylogenetic relationships and physical locations. We found significant excesses of recent duplications in perennial species, compared with those of annuals, represented by rice and Arabidopsis. Consequently, we observed higher nucleotide identity among paralogs and a higher percentage of NBS-encoding genes positioned in numerous clusters in the grapevine and poplar. These results suggested that recent tandem duplication played a major role in NBS-encoding gene expansion in perennial species. These duplication events, together with a higher probability of recombination revealed in this study, could compensate for the longer generation time in woody perennial species e.g. duplication and recombination could serve to generate novel resistance specificities. In addition, we observed extensive species-specific expansion in TIR-NBS-encoding genes. Non-TIR-NBS-encoding genes were poly- or paraphyletic, i.e. genes from three or more plant species were nested in different clades, suggesting different evolutionary patterns between these two gene types.

Recent insights into R gene evolution

SUMMARY Plants are under strong evolutionary pressure to maintain surveillance against pathogens. Resistance (R) gene-dependent recognition of pathogen avirulence (Avr) determinants plays a major role in plant defence. Here we highlight recent insights into the molecular mechanisms and selective forces that drive the evolution of NB-LRR (nucleotide binding-leucine-rich repeat) resistance genes. New implications for models of R gene evolution have been raised by demonstrations that R proteins can detect cognate Avr proteins indirectly by 'guarding' virulence targets, and by evidence that R protein signalling is regulated by intramolecular interactions between different R functional domains. Comparative genomic surveys of NB-LRR diversity in different species have revealed ancient NB-LRR lineages that are unequally represented among plant taxa, consistent with a Birth and Death Model of evolution. The physical distribution of NB-LRRs in plant genomes indicates that tandem and segmental duplication are important factors in R gene proliferation. The majority of R genes reside in clusters, and the frequency of recombination between clustered genes can vary strikingly, even within a single cluster. Biotic and abiotic factors have been shown to increase the frequency of recombination in reporter transgene-based assays, suggesting that external stressors can affect genome stability. Fitness penalties have been associated with some R genes, and population studies have provided evidence for maintenance of ancient R allelic diversity by balancing selection. The available data suggest that different R genes can follow strikingly distinct evolutionary trajectories, indicating that it will be difficult to formulate universally applicable models of R gene evolution.