Genetics, Vol. 163, 735-746, February 2003, Copyright © 2003

Natural Selection for Polymorphism in the Disease Resistance Gene Rps2 of Arabidopsis thaliana

Rodney Mauricio1,a,b, Eli A. Stahl1,a,c, Tonia Korvesa, Dacheng Tiana, Martin Kreitmana,c, and Joy Bergelsona,c
a Department of Ecology and Evolution, University of Chicago, Chicago, Illinois 60637,
b Department of Genetics, University of Georgia, Athens, Georgia 30602-7223
c Committee on Genetics, University of Chicago, Chicago, Illinois 60637

Corresponding author: Joy Bergelson, 1101 E. 57th St., Chicago, IL 60637., jbergels{at}midway.uchicago.edu (E-mail)

Communicating editor: M. AGUADÉ


*  ABSTRACT
*TOP
 *ABSTRACT
*MATERIALS AND METHODS
 *RESULTS
 *DISCUSSION
 *LITERATURE CITED

Pathogen resistance is an ecologically important phenotype increasingly well understood at the molecular genetic level. In this article, we examine levels of avrRpt2-dependent resistance and Rps2 locus DNA sequence variability in a worldwide sample of 27 accessions of Arabidopsis thaliana. The rooted parsimony tree of Rps2 sequences drawn from a diverse set of ecotypes includes a deep bifurcation separating major resistance and susceptibility clades of alleles. We find evidence for selection maintaining these alleles and identify the N-terminal part of the leucine-rich repeat region as a probable target of selection. Additional protein variants are found within the two major clades and correlate well with measurable differences among ecotypes in resistance to the avirulence gene avrRpt2 of the pathogen Pseudomonas syringae. Long-lived polymorphisms have been observed for other resistance genes of A. thaliana; the Rps2 data suggest that the long-term maintenance of phenotypic variation in resistance genes may be a general phenomenon and are consistent with diversifying selection acting in concert with selection to maintain variation.

PLANTS are attacked by a multitude of pathogens and can respond to a subset of these attacks by mounting an induced defense response (BURDON 1987 Down). The initial step in the induction of a defense response involves a genetic interaction between a specific allele of a disease resistance (R) gene and a complementary pathogen avirulence (avr) gene, the so-called gene-for-gene interaction (FLOR 1956 Down, FLOR 1971 Down; STASKAWICZ et al. 1995 Down). In Arabidopsis thaliana, the Rps2 resistance gene confers resistance to pathogens with the avirulence gene avrRpt2 in the pathogen Pseudomonas syringae (DONG et al. 1991 Down; WHALEN et al. 1991 Down; KUNKEL et al. 1993 Down; YU et al. 1993 Down; BENT et al. 1994 Down; MINDRINOS et al. 1994 Down). Recently, P. syringae strains have been found to infect A. thaliana in natural populations (JAKOB et al. 2002 Down).

The RPS2 protein contains a nucleotide-binding site (NBS) and a leucine-rich repeat (LRR) region, two characteristics of a large family of plant R genes (e.g., SALMERON et al. 1996 Down; THOMAS et al. 1997 Down; MCDOWELL et al. 1998 Down; MEYERS et al. 1998 Down; ELLIS et al. 1999 Down; NOEL et al. 1999 Down; BITTNER-EDDY et al. 2000 Down; LUCK et al. 2000 Down). The LRR region is thought to function in pathogen recognition and thereby determine resistance specificity (JONES and JONES 1997 Down; LEISTER and KATAGIRI 2000 Down; TAO et al. 2000 Down; AXTELL et al. 2001 Down). Within the LRR, solvent-exposed amino acid residues framed by conserved aliphatic residues are predicted to make direct contacts with the avirulence gene product or avr gene-dependent factor(s) (JONES and JONES 1997 Down). Evolutionary analyses point to the framed, solvent-exposed residues as exhibiting very fast substitution rates due to positive Darwinian selection (PARNISKE et al. 1997 Down; MEYERS et al. 1998 Down; BITTNER-EDDY et al. 2000 Down; BERGELSON et al. 2001 Down; MONDRAGON-PALOMINO et al. 2002 Down), consistent with their direct role in pathogen (i.e., avirulence gene) recognition. Other regions may also determine recognition (ELLIS et al. 1999 Down; LUCK et al. 2000 Down), however, and R gene-mediated resistance levels can also depend on other host factors (BANERJEE et al. 2001 Down).

Disease resistance genes are often polymorphic for resistance and susceptibility alleles (KUNKEL 1996 Down; STAHL et al. 1999 Down; ELLIS et al. 2000 Down; BERGELSON et al. 2001 Down; HOLUB 2001 Down). CAICEDO et al. 1999 Down examined patterns of polymorphism among eight independent alleles of Rps2 and found evidence of two divergent classes. Statistical tests of the data failed to detect evidence for natural selection, but several features of the data led the authors to suggest that selection, nonetheless, might be important at Rps2. First, the locus contained a high level of nucleotide polymorphism, with almost half of the polymorphisms resulting in amino acid changes. Second, the unrooted gene tree structure included one long branch separating a susceptibility allele (present in accessions Wu-0 and Zu-0-1) from a cluster of more closely related resistance and susceptibility alleles, a structure consistent with balancing selection maintaining Rps2 polymorphism. Finally, the tree indicated a preponderance of amino acid changes between more closely related alleles, suggesting that diversifying selection may have generated Rps2 sequence variation.

Here we extend the results of CAICEDO et al. 1999 Down by carrying out statistical tests of selective neutrality and balancing selection at the Rps2 locus with a larger sample of A. thaliana accessions and a sequence from the closely related congener, A. lyrata. We relate quantitative resistance phenotypes to the evolutionary history of the alleles and identify RPS2 mutations that may confer phenotypic variation. We also test for associations of Rps2 sequence variation and the geographic origin of alleles. The data are discussed in reference to the evolutionary processes thought to underlie plant disease resistance polymorphism.


*  MATERIALS AND METHODS
*TOP
 *ABSTRACT
 *MATERIALS AND METHODS
*RESULTS
 *DISCUSSION
 *LITERATURE CITED

Plant materials:
Twenty-seven accessions of A. thaliana were chosen to create a worldwide sample for Rps2 sequencing representing the major geographic regions in the species' distribution (Table 1). Twelve of these accessions were taken from collections of J. Bergelson and R. Mauricio. Fifteen were selected from those at the Arabidopsis Biological Resource Center (ABRC) at random, except avoiding an excess of accessions from any one country. These accessions were obtained from the ABRC, and seeds from single individuals were harvested to create single-seed stocks for producing the plant materials used in the study. Two individuals of A. lyrata from Indiana (collected by R. Mauricio and D. Jacobson) were used to determine a consensus sequence of the locus for this species.


 
View this table:
In this window
In a new window

  		Table 1. Arabidopsis thaliana accessions and their avrRpt2-dependent resistance phenotypes

Phenotype assessment:
Resistance phenotypes to the P. syringae avirulence gene, avrRpt2, were determined in all but six of the sequenced accessions, as well as the "Columbia" accession and the mutant, rps2-101C (in a Columbia background). Plants were grown from seed in Promix soil with a 12-hr day length at 20°. When the plants were 3–4 weeks old, one entire new leaf was infiltrated with P. syringae pv. tomato strain DC3000 at OD of 0.0002 using a blunt 1-cc syringe. The pathogen strain used in these infections contained a plasmid: either pLAFR3 or pLABL18. The pLABL18 plasmid is identical to the pLAFR3 plasmid, but contains an additional 3.6-kb fragment containing the avrRpt2 gene (WHALEN et al. 1991 Down). Three days postinfection, bacterial levels were measured by grinding standard hole-punch-size leaf punctures in 10 mM MgSO4 and plating dilutions on King's medium B with 40 mg/ml tetracycline. Five to eight replicate plants of each accession were infected with each bacterial strain per experiment. Phenotyping was replicated in at least two experiments for each sequenced accession, except for Po-1 and Mt-0, which were tested in only one experiment.

Plants were identified as resistant or susceptible by comparing the growth (colony-forming units per leaf punch, log-transformed) of the two pathogen strains in an analysis of variance (ANOVA) that included experimental day for accessions evaluated on multiple days. Accessions in which the pathogen strain with pLABL18 grew significantly less than the strain with pLAFR3 (Table 1, column 3) were designated "resistant." Other accessions were designated "susceptible."

For those accessions designated resistant, resistance was quantified by comparing pathogen growth in each focal accession with growth in the Columbia accession. Resistance relative to Columbia (Table 1, column 4) was calculated by dividing the difference in the growth of the strains pLABL18 and pLAFR3 in the focal accession by the difference in the growth of the strains in Columbia, as assessed on the same experimental days. Gaps in the distribution of relative resistance values were used to categorize accessions by degree of resistance. Accessions in the group with lower resistance than that of Columbia were labeled "mildly resistant" (mR), and those in the group with higher values than that of Columbia were labeled "strongly resistant" (sR). To determine the significance level for the degree of resistance relative to Columbia, we evaluated the significance of the interaction between accession (Columbia vs. the focal accession) and the pathogen strain (containing pLAFR3 or pLAB18) using ANOVA.

DNA sequence determination:
For each accession, DNA was extracted from young rosette leaves using protocols described previously (BERGELSON et al. 1998 Down). The region encompassing Rps2 was amplified in three overlapping amplicons, using primers GTTAGTTGGGTGGCGGGAGAG and GGCACAACCGAAACAACTGAGG, AACGGAGACTAAAACAGCCC and GACATGCATCTTCACC, and GTGGATCCATGCTAGTCACATTGAAGTTC and GACCTTTTTATTCCTTTTTCCG, in standard PCR protocols. Both strands were sequenced throughout the region using internal primers (sequences available from the authors) and ABI (Applied Biosystems, Foster City, CA) dye terminator sequencing chemistry. Sequences for each accession were compiled and aligned using Sequencher 3.0 (Genecodes, Ann Arbor, MI). A single consensus sequence for A. lyrata was generated from partial sequences of the two A. lyrata individuals. A small number of sites in our A. lyrata sequences were polymorphic; in each case one of the two alleles included the base found in A. thaliana, and we assigned the A. thaliana base to the consensus A. lyrata sequence for analyses. Multiple large insertions and deletions between the A. thaliana and A. lyrata sequences in the 5' and 3' noncoding regions substantially decreased the number of sites at which between-species comparisons could be made. As a consequence, some polymorphism analyses were conducted without considering the outgroup sequence.

We found differences between our sequences from accessions Wu-0 and Zu-0 and those reported by CAICEDO et al. 1999 Down for the same accessions. In particular, CAICEDO et al. 1999 Down did not detect mutations at positions 1279, 2554, and 3085 (Fig 1), and they found that variants at positions 3396 and 3502 were shared by Wu-0 and Zu-0 (Fig 1). Variation within accessions has been noted (BREYNE et al. 1999 Down; CAICEDO et al. 1999 Down) and may reflect ecotype seed collection from multiple (or heterozygous) individuals in the field.



View larger version (34K):
In this window
In a new window
Download PPT slide
  		Figure 1. Rps2 region polymorphic sites. Shown are positions in the alignment and the bases of A. lyrata (consensus) and A. thaliana sequences. Periods represent ancestral bases inferred from A. lyrata, and bases indicate derived polymorphic mutations. Amino acid replacement mutations are indicated relative to the wild-type Columbia accession (GenBank accession no. AL049483), which is identical to Bla-2, C2-1, and Gott-20. The region encoding the RPS2 leucine-rich repeat (LRR) region is indicated by the line above the amino acid replacement mutations, with the thicker line indicating its 5' half.

Population genetic analyses:
Silent (noncoding and synonymous) and amino acid replacement (nonsynonymous) polymorphism and divergence (Jukes-Cantor corrected) calculations were conducted using DnaSP (ROZAS and ROZAS 1997 Down). Genealogy estimation was conducted by parsimony using PAUP (SWOFFORD 1996 Down), with 500 bootstrap replicates. Standard tests of a panmictic population, neutral mutation model utilized coalescent simulations with a fixed number of segregating sites (HUDSON 1993 Down), with programs available from R. Hudson (http://home.uchicago.edu/~rhudson1). Analyses testing for heterogeneity of polymorphism to divergence ratios were conducted using DNASlider (MCDONALD 1998 Down); sliding window average G-statistics were analyzed for scaled population recombination rates, from RSlider = 0 to 100 (RSlider = 4NerL, where Ne is the effective population size, r is the recombination rate per base pair per generation and L is the length of the analyzed region), with the most conservative P values obtained for RSlider = 7 for the entire sequenced region and RSlider = 6 for the coding sequence. Linkage disequilibrium was tested for in 2 x 2 contingency tables by Fisher's exact tests, using shareware available from W. Engels. Differentiation among "populations" (groups of sequences defined by phenotype or geographic origin) was calculated as FST = 1 - {pi}W/{pi}T, where {pi}W is average pairwise nucleotide difference within populations and {pi}T is that in total (HUDSON et al. 1992B Down), and was tested by resampling with sequences permuted across groups (following HUDSON et al. 1992A Down; HOLSINGER and MASON-GAMER 1996 Down; BERGELSON et al. 1998 Down), using programs written by E. A. Stahl.

We analyzed a coalescent model with selection and recombination as described in TIAN et al. 2002 Down. Selection was assumed to maintain two alleles at fixed frequency 0.81 (and 0.19), acting at the beginning of the LRR (site 2654 in our alignment). The scaled per-base-pair recombination rate R = 4Ne(1 - s)r = 0.00057 uses published estimates of effective population size and selfing rate for A. thaliana (see TIAN et al. 2002 Down) and a recombination rate per meiosis estimated from regression of genetic and physical positions of markers near Rps2 (2.71 cM/Mb, r2 = 0.96; data from the Arabidopsis Genome Resource, http://ukcrop.net/agr, markers mi475, SEP2B, m600, PG11, DD1, mi123, RLK5, mi232, prha, g8300, and mi431). The scaled mutation rate between selected alleles (0.0125) was adjusted to fit the observed data near the selected site.


*  RESULTS
*TOP
 *ABSTRACT
 *MATERIALS AND METHODS
 *RESULTS
*DISCUSSION
 *LITERATURE CITED

avrRpt2-dependent resistance phenotypes:
For each of 21 accessions, we compared the growth of P. syringae strain DC3000 with avrRpt2 and DC3000 without avrRpt2. If an accession is resistant, the growth of the strain with avrRpt2 should be significantly less than the growth of the strain without avrRpt2. The log of growth of the pathogen without avrRpt2 minus that of the pathogen with avrRpt2 is listed in Table 1; this measure of resistance is unitless since it is equivalent to the log of the ratio of growth for the two pathogen strains. The results of our ANOVAs indicate that 17 of the 21 accessions tested were resistant. Accessions BG-4, Po-1, Zu-0, and Knox-2 and the Columbia rps2 mutant showed no indication of resistance. Statistical designations of resistance and susceptibility were consistent with observed disease symptoms.

We determined whether resistant accessions inhibited bacterial growth of DC3000 with avrRpt2 to different extents by comparing bacterial growth in each line relative to this same measure in a common paired control line, Columbia. Relative resistance values (Table 1) ranged from 0.285 to 1.61. Gaps in the distribution of relative resistance values, between 0.67 and 0.945 and between 1.14 and 1.39, allowed us to group alleles into three operational subclasses of resistance, mild (mR), intermediate (R), and strong (sR). We used relative resistance values rather than ANOVA P values (Table 1) to categorize accessions because the power to detect differences from Columbia varied among accessions. The mR group included accessions AB-7, GR-6, Wu-0, Yo-0, and Cvi-0, and the sR group included Pog-0, RLD-1, Co-1, and Tsu-0.

Low growth of DC3000 without avrRpt2 in Pu-8 suggested partial resistance to the DC3000 background; additional resistance in the presence of avrRpt2 indicated that Pu-8 is resistant, but we were unable to measure its relative resistance. Also, Wu-0 has been reported previously as susceptible (CAICEDO et al. 1999 Down) although it exhibits growth and symptoms consistent with an intermediate phenotype (KUNKEL et al. 1993 Down; this study). It is possible that CAICEDO et al. 1999 Down studied a different genotype within Wu-0 (see MATERIALS AND METHODS).

Molecular variation at Rps2:
We surveyed DNA sequence variability in 27 accessions from throughout the species range, including the accessions whose resistance phenotypes we determined (Table 1), and from the closely related species A. lyrata. The sequenced region spans 4248 base pairs (bp) in A. thaliana accession Columbia (GenBank accession no. AL049483), from 1003 bp upstream of the Rps2 start codon to 521 bp downstream of its stop codon. Our survey yielded a 4461-bp alignment including the A. lyrata sequence, with 3755 sites at which polymorphism and divergence were ascertained (Table 2).


 
View this table:
In this window
In a new window

  		Table 2. Levels of variability across the Rps2 locus and among RPS2 functional domains

The data including the outgroup sequence revealed a total of 197 nucleotide differences fixed between A. lyrata and all A. thaliana sequences and 58 single nucleotide polymorphisms distinguishing 18 haplotypes in the 27 A. thaliana alleles (Fig 1; Table 2). Within the Rps2 coding sequence, we detected 55 nonsynonymous (amino acid changing) differences between species and 20 nonsynonymous polymorphisms. The Rps2 coding sequence reading frame is intact in all individuals, despite two one-codon insertions in A. lyrata relative to A. thaliana at Columbia residues 741 and 771 (both in the LRR region) and a four-codon deletion at 877 (near the RPS2 C terminus). We also introduced one-codon insertion/deletions (indels; in both A. lyrata and A. thaliana) at Columbia residues 86 and 737, where the two species differ at all three nucleotide positions; these three-base differences were not included in polymorphism analysis. We found numerous indels between species in noncoding regions and five indel polymorphisms all outside of the coding sequence. A homonucleotide run at 821 varied between two A. lyrata individuals, but in A. thaliana no microsatellites were detected. No heterozygous sites were detected in A. thaliana individuals. Overall levels of polymorphism and divergence at Rps2 (Table 2) fall within the range seen at other loci in A. thaliana and A. lyrata (KAWABE and MIYASHITA 1999 Down; PURUGGANAN and SUDDITH 1999 Down; AGUADE 2001 Down).

In Table 3, levels of polymorphism and divergence in the LRR region are presented. Within this region, the ß-pleated sheet structural motif consensus sequence (JONES and JONES 1997 Down) allows framed solvent-exposed amino acid residues, specific candidates for positive selection, to be analyzed and compared with conserved structural residues and nonconserved residues between frames. Significantly greater Ka than Ks between R gene paralogs at framed exposed residues (MEYERS et al. 1998 Down; BERGELSON et al. 2001 Down) has provided strong evidence for positive selection on plant R genes. In contrast, synonymous and nonsynonymous divergence reveals no evidence for positive selection on Rps2 (framed exposed residues, Ka = 0.033, Ks = 0.12); functional constraint is evident for all categories of LRR region amino acid residues (Table 3). Contingency tables comparing synonymous and replacement polymorphism and divergence (MCDONALD and KREITMAN 1991 Down) also provide no indication of selectively driven protein evolution (Table 3).


 
View this table:
In this window
In a new window

  		Table 3. Polymorphism and divergence within the RPS2 leucine-rich repeat (LRR) region

Evidence for balancing selection at Rps2:
Fig 2 shows a parsimony tree inferred from silent and nonsynonymous polymorphism and divergence, with accession name and avrRpt2-dependent resistance phenotype shown for each allele. The Rps2 gene tree reveals the presence of two highly supported major clades. This haplotype structure is evident for synonymous as well as amino acid replacement polymorphisms, but only for polymorphisms falling in the middle of the coding sequence. Tests for nonrandom associations between all pairs of nonsingleton polymorphisms (Fig 3) reveal that linkage disequilibrium is clustered within a central segment of the Rps2 coding sequence. Indeed only polymorphisms in this segment show significant linkage disequilibrium after correction for multiple tests of association. Outside of this central segment of the Rps2 coding sequence, the data reveal little haplotype structure (Fig 1 and Fig 3).



View larger version (17K):
In this window
In a new window
Download PPT slide
  		Figure 2. Phylogeny of Rps2 sequences based on parsimony analysis of silent, synonymous, and amino acid replacement variability. Accession names are indicated for each Rps2 sequence with avrRpt2-dependent phenotype in boldface type (R, resistant; mR, mildly resistant; sR, strongly resistant; S, susceptible). The tree shown is one of three most parsimonious trees (length 265, consistency index 0.974) that differ only in the resistance (upper) clade. Numbers of mutations are shown above branches, with proportional branch length. Bootstrap values >90% are shown below branches (boldface italics).



View larger version (77K):
In this window
In a new window
Download PPT slide
  		Figure 3. Linkage disequilibrium between polymorphisms in the Rps2 region. The Rps2 region diagram shows the coding sequence (box) with RPS2 functional regions (LZ, leucine zipper; NBS, nucleotide-binding site; LRR, leucine-rich repeat). On the horizontal lines below, singleton polymorphisms (small hash marks) and nonsingleton polymorphisms (sample frequency two or greater, large hash marks) are indicated for silent/synonymous polymorphisms (top line) and amino acid replacement polymorphisms in the coding sequence (bottom line). In the triangle at bottom, Fisher's exact test P values for each pair of nonsingleton polymorphic sites are indicated by shading, P > 10-2 (white), 10-3 < P < 10-2 (stippled), 10-4 < P < 10-3 (shaded stipple), and P < 10-4 (black). Only P values <10-4 (black) remain significant after Bonferroni correction.

Sliding window analysis of nucleotide diversity between the two major clades (Fig 4) shows a peak of silent polymorphism in the center of the coding sequence—the 300 bp 5' of the region encoding the RPS2 LRR region and in the 5' half of the region encoding the RPS2 LRR region itself (hereafter referred to as the 5' LRR region)—corresponding to the region containing the cluster of polymorphic sites in linkage disequilibrium. Peak nucleotide diversity between the two major clades reaches {pi}b = 0.086 in the Rps2 5' LRR region, a value approaching estimates of silent divergence between species. Clustering of silent polymorphism within this segment of the Rps2 coding sequence results in significant heterogeneity in the ratio of polymorphism to divergence across the sequenced region (sliding window average G: entire region, P <= 0.004; coding sequence, P <= 0.0014; MCDONALD 1998 Down). Variation at Rps2, therefore, is not compatible with an equilibrium model of selective neutrality in a panmictic population.



View larger version (22K):
In this window
In a new window
Download PPT slide
  		Figure 4. Sliding window analysis of silent (noncoding and synonymous) divergence between resistance and susceptibility clades of Rps2 alleles. Average numbers of pairwise differences per site within the window are shown with a solid line. Predicted levels under a coalescent model with selection and recombination (dashed line) assume that selection acts at the beginning of the LRR region (2654) and maintains Rps2 polymorphism at frequency 22/27 = 0.81, with independently estimated recombination rate 0.00057 and fitted mutation rate between selected alleles 0.0125 (see MATERIALS AND METHODS). Expected levels under neutrality (dotted line) are calculated as divergence in the window times the ratio of averages across the region of polymorphism ({pi}) and divergence, multiplied by the expected time to the most recent common ancestor for sample size 27 relative to expected average pair-wise coalescence time, The window is 150 silent sites wide, slid by 10-site increments. Beneath the sliding window plot the corresponding functional regions of RPS2 are shown, with amino acid differences between the clades indicated by asterisks (*).

Alleles from resistant and susceptible accessions are not scattered throughout the Rps2 gene tree, but are grouped together; therefore, we refer to the two major clades as the resistance (R) clade and susceptibility (S) clade. We tested for a significant association between Rps2 sequence variation and avrRpt2-dependent resistance variation by analyzing differentiation (an FST estimator based on nucleotide diversities; HUDSON et al. 1992B Down) between phenotypes. Overall differentiation between phenotypes is highly significant (S, mR, R, and sR; FST = 0.52, P < 0.001). Pairwise comparisons between phenotypes reveal significant differentiation between S and each of R, mR, and sR (FST >= 0.47, P <= 0.019), marginally significant differentiation between R and mR phenotypes (FST = 0.12, P = 0.09), and no significant differentiation for other comparisons between resistant phenotypes (FST < 0.005, P > 0.3). Thus, sequence variation at Rps2 correlates with avrRpt2-dependent disease resistance, suggestive of causal links between the two (see DISCUSSION).

Geographic differentiation:
In contrast to avrRpt2-dependent resistance, accessions from the same geographic region are scattered throughout the Rps2 gene tree (Fig 2). We categorized accessions into five regions, (1) Eastern Europe, Asia, and Africa; (2) Central and Northern Europe; (3) Western and Southern Europe; (4) Eastern North America; and (5) Western North America, on the basis of the recent expansion of A. thaliana from Western Asia and Eastern Europe to its current worldwide distribution (PRICE et al. 1994 Down; SHARBEL et al. 2000 Down). Rps2 sequence variation reveals no differentiation among regions (overall FST = 0.043, P = 0.3; for all pairs of regions FST <= 0.14, P >= 0.15). In addition, Rps2 variation does not differentiate North America from other continents (FST = 0.043, P > 0.5), revealing no evidence for a founder effect in the colonization of the Western hemisphere by A. thaliana. These results are typical of studies of a single or few loci and moderate sample sizes in A. thaliana (INNAN et al. 1996 Down; BERGELSON et al. 1998 Down; KAWABE and MIYASHITA 1999 Down).


*  DISCUSSION
*TOP
 *ABSTRACT
 *MATERIALS AND METHODS
 *RESULTS
 *DISCUSSION
*LITERATURE CITED

Previously, CAICEDO et al. 1999 Down found a high level of polymorphism at the Rps2 locus and two highly divergent alleles, suggestive of balancing selection, but a statistical test (Tajima's D) could not reject selective neutrality. Here we find statistical evidence in support of the selection hypothesis and tentatively identify the Rps2 5' LRR region as the target of selection. An Rps2 sequence from sister species A. lyrata and a larger sampling of alleles allowed us to detect a clustering of polymorphism relative to divergence exceeding that possible under selective neutrality in a panmictic population. This result rules out the possibility that the region of high polymorphism is a mutational hotspot, since mutation rate heterogeneity would affect both polymorphism and divergence.

Our statistical confirmation of a peak of polymorphism should not be taken, in and of itself, as a strong refutation of neutral evolution. For example, INNAN et al. 1996 Down identified a short highly diverged stretch in exon 4 of Adh, as well as in three adjacent sequence "blocks." While the authors argue in favor of balancing selection acting on exon 4 (owing to amino acid replacement differences between the alleles), they raise the possibility that population structure and history produced the "dimorphism" seen throughout the locus. Biallelic variation has also been found at several other loci in A. thaliana (KAWABE et al. 1997 Down; KAWABE and MIYASHITA 1999 Down; STAHL et al. 1999 Down; PURUGGANAN and SUDDITH 1999 Down; AGUADE 2001 Down; HAUSER et al. 2001 Down; TIAN et al. 2002 Down), adding to the appeal of a population structure hypothesis.

We favor balancing selection as an explanation for Rps2 variation, on the basis of features of the data that distinguish our results from those of other studies that find biallelic variation but favor a population structure hypothesis. As indicated in the sliding window analysis (Fig 4), most of the variation is present in the coding segment of the Rps2 gene and overlaps with the functional domain of the protein implicated in pathogen recognition. Seven amino acid replacement changes separate the R and S clades, four in the LRR region, and the suggestion that differences between Rps2 allelic classes could be functional is consistent with a role of selection. Furthermore, accessions' Rps2 allelic classes correspond closely with their resistance phenotypes. Since selection can act only if functionally distinct alleles exist, a correspondence between phenotype and genotype provides additional evidence in support of balancing selection. Others have also pointed to the importance of possible functional differences distinguishing diverged alleles. For example, HAUSER et al. 2001 Down found two divergent alleles across part of the region in their analysis of polymorphism in Glabrous1 (Gl1), a candidate gene for leaf trichome density variation; they argued against selection because the divergence was not in the coding region of the gene and variation in trichome density did not correlate with Gl1 sequence variation. KAWABE et al. 2000 Down found divergent alleles of the cytosolic phosphoglucose isomerase (PgiC) and favored balancing selection because the alleles produced distinct allozymes (but note that phenotypic properties of the allozymes were not investigated).

Balancing selection is expected to lead to a signature in which neutral variation accumulates between the alleles surrounding the site(s) under selection. This signature is a simple manifestation of the genealogical correlation of tightly linked sites: as a balanced polymorphism becomes old, so too do the genealogical ages of sites tightly linked to it. In HUDSON and KAPLAN's (1988) coalescent treatment of balanced polymorphism, the physical scale of neutral polymorphism linked to the site under selection is, to a first approximation, determined by a balance between the origination of new neutral mutations (governed by the scaled neutral mutation rate, 4Neu, where Ne is the effective population size and u is the neutral mutation rate per site per generation) and the decay of the linkage disequilibrium between these mutations and the site under selection (governed by the scaled recombination rate, 4Ner, where r is the per generation recombination rate between adjacent sites). Even for a highly self-fertilizing species, balancing selection can be expected to produce a relatively sharp peak of neutral polymorphism linked to a site under selection (NORDBORG et al. 1996 Down; NORDBORG 1997 Down). Based on available genetic and population genetic estimates of mutation and recombination rates in A. thaliana (TIAN et al. 2002 Down; see MATERIALS AND METHODS and Fig 4), the peak of polymorphism seen at Rps2 is compatible with theoretical predictions for a balanced polymorphism at the 5' end of the region of the gene that encodes the RPS2 LRR region (Fig 4).

We note that the balancing selection analysis is based on a constant-size panmictic population model and does not take into account departures from this model in the demographic history of A. thaliana. Nevertheless, given that the peak of polymorphism is restricted to within the Rps2 coding sequence, that polymorphisms within the peak are not in linkage disequilibrium with polymorphisms outside it, and that significant linkage disequilibrium is rarely observed between loci in A. thaliana (INNAN et al. 1997 Down; NORDBORG et al. 2002 Down), we can identify the region of the peak, which includes the N-terminal half of the RPS2 LRR region, as the target of natural selection.

In our balancing selection analysis, the best-fit mutation rate between allelic classes was found to be equal to 0.01, i.e., on the order of one-hundredth the rate of neutral coalescence (1/Ne). Higher mutation rates between the selected alleles would lead to more recent common ancestry between them, and if large enough may not result in an observable peak of polymorphism even with balancing selection. Many kinds of mutations can cause loss of function; therefore the rate of origination of new susceptibility alleles might be expected to be quite high. An ancient balanced polymorphism between a resistance and a susceptibility allele would imply that selection favors one susceptibility allele over others and that the rate of origination of this particular susceptibility allele is low. Alternative resistance alleles, on the other hand, might be expected to have a low rate of origination. Thus, the observation of a signature of selection between the two major Rps2 clades is consistent with the hypothesis that the two major allelic classes of Rps2 contain functional resistance alleles. Indeed, BANERJEE et al. 2001 Down showed that the susceptibility allele of Po-1 is partially functional against avrRpt2 when in the Col-0 genetic background. Note that the designation of resistance or susceptibility in this study is based only on the ability to recognize one specific avirulence gene, avrRpt2. We propose that the alleles represented by the Rps2 resistance and susceptibility clades encode distinct specificities against natural pathogens in wild populations. The recent finding of infection by P. syringae in natural populations of A. thaliana (JAKOB et al. 2002 Down) makes this a realistic possibility.

Rps2 exhibits marked sequence variability in association with phenotypic variation. Seven of the nine phenotypic changes that would be inferred by simply mapping phenotypes onto the Rps2 gene tree are associated with amino acid polymorphisms, six with polymorphisms in the LRR region (Fig 1 and Fig 2). Polymorphisms that distinguish the R and S clades are found upstream of the LRR region (not shown) and in nonconserved residues between LRR frames (Fig 5); these changes could confer phenotypic variation that is maintained by natural selection (ELLIS et al. 1999 Down; LUCK et al. 2000 Down). Polymorphisms associated with other phenotypic changes on the tree include framed solvent-exposed residues and conserved residues between frames (Fig 5). While we cannot rule out the possibility that changes at other loci contribute to phenotypic variation in these accessions, we suggest that these amino acid polymorphisms should be candidates for further study of RPS2 function (AXTELL et al. 2001 Down). Moreover, besides conferring phenotypic variation that is maintained by selection, hypervariability of amino acid residues N-terminal to the LRR region and in the N-terminal half of the LRR region may be consistent with diversifying selection on Rps2.



View larger version (44K):
In this window
In a new window
Download PPT slide
  		Figure 5. The RPS2 LRR region, with polymorphic mutations. The amino acid sequence taken from JONES and JONES 1997 Down is shown, and codon number of the rightmost residue in each row is shown on the right. Residues matching the LRR consensus (at bottom) are shown in boldface type, and the vertical lines bracket the structural motif frame. Residues that differ between the resistance and susceptibility clades are indicated in red, and residues that differ in association with phenotypic changes within the clades are indicated in blue.

Previous studies have found evidence for rapid adaptive substitution rates in LRR region solvent-exposed residues among R gene paralogs (MEYERS et al. 1998 Down; BERGELSON et al. 2001 Down; HOLUB 2001 Down; MONDRAGON-PALOMINO et al. 2002 Down). In contrast, previous studies have not found evidence for positive selection at two Arabidopsis R genes that exhibit signatures of balancing selection, Rpm1 (STAHL et al. 1999 Down) and Rps5 (TIAN et al. 2002 Down). The possibility that the major alleles of Rps2 represent a functional balanced polymorphism suggests that the maintenance of variation by natural selection may be a general feature of R gene evolutionary dynamics. At Rps2, we do not find evidence for adaptive protein evolution between species, but we do observe marked amino acid variability that could be consistent with diversifying selection. It remains to be seen whether even faster evolving R genes can also support balanced polymorphisms, as exemplified by genes of mammalian major histocompatibility complex (HUGHES and NEI 1988 Down) and plant self-incompatibility loci (CLARK 1993 Down).


*  FOOTNOTES

Sequence data from this article have been deposited with the EMBL/Genbank Data Libraries under accession nos. AF487796–AF487823. Back
1 These authors contributed equally to this work. Back


*  ACKNOWLEDGMENTS

We thank A. Berry, I. Cetl, M. Nachman, G. Robellen, O. Savolainen, and J. Winterer for collecting A. thaliana seeds, as well as the Arabidopsis Biological Resource Center at Ohio State University for providing seeds of A. thaliana accessions. F. Ausubel provided seed of rps2-101C. We acknowledge the assistance of M. Aguadé and the reviewers who provided careful and helpful reviews. This work was funded by a Sloan Foundation/National Science Foundation Fellowship in Molecular Evolution and University of Georgia Faculty Research grant to R.M., a Sloan Foundation/Department of Energy Fellowship in Computational Molecular Biology to E.A.S., and a Packard Fellowship and National Institutes of Health awards GM-57994 and GM-62504 to J.B.

Manuscript received August 7, 2002; Accepted for publication November 11, 2002.


*  LITERATURE CITED
*TOP
 *ABSTRACT
 *MATERIALS AND METHODS
 *RESULTS
 *DISCUSSION
 *LITERATURE CITED

AGUADÉ, M., 2001  Nucleotide sequence variation at two genes of the phenylpropanoid pathway, the FAH1 and F3H genes, in Arabidopsis thaliana.. Mol. Biol. Evol. 18:1-9.[Abstract/Free Full Text]

AXTELL, M. J., T. W. MCNELLIS, M. B. MUDGETT, C. S. HSU, and B. J. STASKAWICZ, 2001  Mutational analysis of the Arabidopsis RPS2 disease resistance gene and the corresponding Pseudomonas syringae avrRpt2 avirulence gene. Mol. Plant-Microbe Interact. 14:181-188.[Medline]

BANERJEE, D., X. ZHANG, and A. F. BENT, 2001  The leucine-rich repeat domain can determine effective interaction between RPS2 and other host factors in Arabidopsis RPS2-mediated disease resistance. Genetics 158:439-450.[Abstract/Free Full Text]

BENT, A. F., B. N. KUNKEL, D. DAHLBECK, K. L. BROWN, and R. SCHMIDT et al., 1994  RPS2 of Arabidopsis thaliana: a leucine-rich repeat class of plant disease resistance genes. Science 265:1856-1860.[Abstract/Free Full Text]

BERGELSON, J., E. A. STAHL, S. DUDEK, and M. KREITMAN, 1998  Genetic variation within and among populations. Genetics 148:1311-1323.[Abstract/Free Full Text]

BERGELSON, J., M. KREITMAN, E. A. STAHL, and D. TIAN, 2001  Evolutionary dynamics of plant R-genes. Science 292:2281-2285.[Abstract/Free Full Text]

BITTNER-EDDY, P. D., I. R. CRUTE, E. B. HOLUB, and J. L. BEYNON, 2000  RPP13 is a simple locus in Arabidopsis thaliana for alleles that specify downy mildew resistance to different avirulence determinants in Peronospora parasitica.. Plant J. 21:177-188.[Medline]

BREYNE, P., D. ROMBAUT, A. VAN GYSEL, M. VAN MONTAGU, and T. GERATS, 1999  AFLP analysis of genetic diversity within and between Arabidopsis thaliana ecotypes. Mol. Gen. Genet. 261:627-634.[Medline]

BURDON, J. J., 1987 Diseases and Plant Population Biology. Cambridge University Press, Cambridge, UK.

CAICEDO, A. L., B. A. SCHAAL, and B. N. KUNKEL, 1999  Diversity and molecular evolution of the RPS2 resistance gene in Arabidopsis thaliana.. Proc. Natl. Acad. Sci. USA 96:302-306.[Abstract/Free Full Text]

CLARK, A. G., 1993 Evolutionary inferences from molecular characterization of self-incompatibility alleles, pp. 79–108 in Mechanisms of Molecular Evolution, edited by N. TAKAHATA and A. G. CLARK. Sinauer Associates, Sunderland, MA.

DONG, X., M. MINDRINOS, K. R. DAVIS, and F. M. AUSUBEL, 1991  Induction of Arabidopsis defense genes by virulent and avirulent Pseudomonas syringae strains and by a cloned avirulence gene. Plant Cell 3:61-72.[Abstract/Free Full Text]

ELLIS, J., P. DODDS, and T. PRYOR, 2000  Structure, function and evolution of plant disease resistance genes. Curr. Opin. Plant Biol. 3:278-284.[Medline]

ELLIS, J. G., G. J. LAWRENCE, J. E. LUCK, and P. N. DODDS, 1999  Identification of regions in alleles of the flax rust resistance gene L that determine differences in gene-for-gene specificity. Plant Cell 11:495-506.[Abstract/Free Full Text]

FLOR, H. H., 1956  The complementary genic systems in flax and flax rust. Adv. Genet. 8:29-54.

FLOR, H. H., 1971  Current status of the gene-for-gene concept. Annu. Rev. Phytopathol. 9:275-296.

HAUSER, M.-T., B. HARR, and C. SCHLOTTERER, 2001  Trichome distribution in Arabidopsis thaliana and its close relative Arabidopsis lyrata: molecular analysis of the candidate gene GLABROUS1.. Mol. Biol. Evol. 18:1754-1763.[Abstract/Free Full Text]

HOLSINGER, K. E. and R. J. MASON-GAMER, 1996  Hierarchical analysis of nucleotide diversity in geographically structured populations. Genetics 142:629-639.[Abstract]

HOLUB, E. B., 2001  The arms race is ancient history in Arabidopsis, the wildflower. Nat. Rev. Genet. 2:516-527.[Medline]

HUDSON, R. R., 1993 The how and why of generating gene genealogies, pp. 23–36 in Mechanisms of Molecular Evolution, edited by N. TAKAHATA and A. G. CLARK. Sinauer Associates, Sunderland, MA.

HUDSON, R. R. and N. L. KAPLAN, 1988  The coalescent process in models with selection and recombination. Genetics 120:831-840.[Abstract/Free Full Text]

HUDSON, R. R., D. D. BOOS, and N. L. KAPLAN, 1992a  A statistical test to detect geographic subdivision. Mol. Biol. Evol 9:138-151.[Abstract]

HUDSON, R. R., M. SLATKIN, and W. P. MADISON, 1992b  Estimation of levels of gene flow from DNA sequence data. Genetics 132:583-589.[Abstract]

HUGHES, A. L. and M. NEI, 1988  Pattern of nucleotide substitution at major histocompatibility complex class-I loci reveals overdominant selection. Nature 335:167-170.[Medline]

INNAN, H., F. TAJIMA, R. TERAUCHI, and N. T. MIYASHITA, 1996  Intragenic recombination in the Adh locus of the wild plant Arabidopsis thaliana.. Genetics 143:1761-1770.[Abstract]

INNAN, H., R. TERAUCHI, and N. T. MIYASHITA, 1997  Microsatellite polymorphism in natural populations of the wild plant Arabidopsis thaliana.. Genetics 146:1441-1452.[Abstract]

JAKOB, K., E. M. GOSS, H. ARAKI, T. VAN, and M. KREITMAN et al., 2002  Pseudomonas viridiflava and P. syringae—natural pathogens of Arabidopsis thaliana. Mol. Plant-Microbe Interact. 15:1195-1203.[Medline]

JONES, D. A. and J. D. G. JONES, 1997  The role of leucine-rich repeat proteins in plant defenses. Adv. Bot. Res. 24:89-167.

KAWABE, A. and N. T. MIYASHITA, 1999  DNA variation in the basic chitinase locus (ChiB) region of the wild plant Arabidopsis thaliana.. Genetics 153:1445-1453.[Abstract/Free Full Text]

KAWABE, A., H. INNAN, R. TERAUCHI, and N. T. MIYASHITA, 1997  Nucleotide polymorphism in the acidic chitinase locus (ChiA) region of the wild plant Arabidopsis thaliana.. Mol. Biol. Evol. 14:1303-1315.[Abstract]

KAWABE, A., K. YAMANE, and N. T. MIYASHITA, 2000  DNA polymorphism at the cytosolic phosphoglucose isomerase (PgiC) locus of the wild plant Arabidopsis thaliana. Genetics 156:1339-1347.[Abstract/Free Full Text]

KUNKEL, B. N., 1996  A useful weed put to work: genetic analysis of disease resistance in Arabidopsis thaliana.. Trends Genet. 12:63-69.[Medline]

KUNKEL, B. N., A. F. BENT, D. DAHLBECK, R. W. INNES, and B. J. STASKAWICZ, 1993  RPS2, an Arabidopsis disease resistance locus specifying recognition of Pseudomonas syringae strains expressing the avirulence gene avrRpt2.. Plant Cell 5:865-875.[Abstract/Free Full Text]

LEISTER, R. T. and F. KATAGIRI, 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.[Medline]

LUCK, J. E., G. J. LAWRENCE, P. N. DODDS, K. W. SHEPHERD, and J. G. ELLIS, 2000  Regions outside of the leucine-rich repeats of flax rust resistance proteins play a role in specificity determination. Plant Cell 12:1367-1377.[Abstract/Free Full Text]

MCDONALD, J. H., 1998  Improved tests for heterogeneity across a region of DNA sequence in the ratio of polymorphism to divergence. Mol. Biol. Evol 15:377-384.[Abstract]

MCDONALD, J. H. and M. KREITMAN, 1991  Adaptive protein evolution at the Adh locus in Drosophila. Nature 351:652-654.[Medline]

MCDOWELL, J. M., M. DHANDAYDHAM, T. A. LONG, M. G. M. AARTS, and S. GOFF et al., 1998  Intragenic recombination and diversifying selection contribute to the evolution of downy mildew resistance at the RPP8 locus of Arabidopsis. Plant Cell 10:1861-1874.[Abstract/Free Full Text]

MEYERS, B. C., K. A. SHEN, P. ROHANI, B. S. GAUT, and R. W. MICHELMORE, 1998  Receptor-like genes in the major resistance locus of lettuce are subject to divergent selection. Plant Cell 10:1833-1846.[Abstract/Free Full Text]

MINDRINOS, M., F. KATAGIRI, G.-L. YU, and F. M. AUSUBEL, 1994  The A. thaliana disease resistance gene RPS2 encodes a protein containing a nucleotide-binding site and leucine-rich repeats. Cell 78:1089-1099.[Medline]

MONDRAGÓN-PALOMINO, M., B. C. MEYERS, R. W. MICHELMORE, and B. S. GAUT, 2002  Patterns of positive selection in the complete NBS-LRR gene family of Arabidopsis thaliana.. Genome Res. 12:1305-1315.[Abstract/Free Full Text]

NOEL, L., T. L. MOORES, E. A. VAN DER BIEZEN, M. PARNISKE, and M. J. DANIELS et al., 1999  Pronounced intraspecific haplotype divergence at the RPP5 complex disease resistance locus of Arabidopsis. Plant Cell 11:2099-2112.[Abstract/Free Full Text]

NORDBORG, M., 1997  Structured coalescent processes on different time scales. Genetics 146:1501-1514.[Abstract]

NORDBORG, M., B. CHARLESWORTH, and D. CHARLESWORTH, 1996  Increased levels of polymorphism surrounding selectively maintained sites in highly selfing species. Proc. Roy. Soc. Lond. Ser. B 263:1033-1039.[Abstract/Free Full Text]

NORDBORG, M., J. O. BOREVITZ, J. BERGELSON, C. C. BERRY, and J. CHORY et al., 2002  The extent of linkage disequilibrium in the highly selfing species Arabidopsis thaliana.. Nat. Genet. 30:190-193.[Medline]

PARNISKE, M., K. E. HAMMOND-KOSACK, C. GOLSTEIN, C. M. THOMAS, and D. A. JONES et al., 1997  Novel disease resistance specificities result from sequence exchange between tandemly repeated genes at the Cf-4/9 locus of tomato. Cell 91:821-832.[Medline]

PRICE, R. A., J. D. PALMER and I. A. AL-SHEHBAZ, 1994 Systematic relationships of Arabidopsis: a molecular and morphological perspective, pp. 7–19 in Arabidopsis, edited by E. M. MEYEROWITZ and C. R. SOMERVILLE. Cold Spring Harbor Laboratory Press, Cold Spring Harbor, NY.

PURUGGANAN, M. D. and J. I. SUDDITH, 1999  Molecular population genetics of floral homeotic loci: departures from the equilibrium-neutral model at the APETALA3 and PISTILLATA genes of Arabidopsis thaliana.. Genetics 151:839-848.[Abstract/Free Full Text]

ROZAS, J. and R. ROZAS, 1997  DnaSP version 2.0: a novel software package for extensive molecular population genetics analysis. Comput. Appl. Biosci. 13:307-311.

SALMERON, J. M., G. E. D. OLDROYD, C. M. T. ROMMENS, S. R. SCOFIELD, and H. S. KIM et al., 1996  Tomato Prf is a member of the leucine-rich repeat class of plant disease resistance genes and lies embedded within the Pto kinase gene cluster. Cell 86:123-133.[Medline]

SHARBEL, T. F., B. HAUBOLD, and T. MITCHELL-OLDS, 2000  Genetic isolation by distance in Arabidopsis thaliana: biogeography and post-glacial colonization of Europe. Mol. Ecol. 9:2109-2118.[Medline]

STAHL, E. A., G. DWYER, R. MAURICIO, M. KREITMAN, and J. BERGELSON, 1999  Dynamics of disease resistance polymorphism at the Rpm1 locus of Arabidopsis. Nature 400:667-671.[Medline]

STASKAWICZ, B. J., F. M. AUSUBEL, B. J. BAKER, J. G. ELLIS, and J. D. G. JONES, 1995  Molecular genetics of plant disease resistance. Science 292:661-667.

SWOFFORD, D., 1996 PAUP: Phylogenetic Analysis Using Parsimony (and Other Methods), Version 4. Sinauer Associates, Sunderland, MA.

TAO, Y., F. YUAN, R. T. LEISTER, F. M. AUSUBEL, and F. KATAGIRI, 2000  Mutational analysis of the Arabidopsis nucleotide binding site-leucine-rich repeat resistance gene RPS2.. Plant Cell 12:2541-2554.[Abstract/Free Full Text]

THOMAS, C. M., D. A. JONES, M. PARNISKE, K. HARRISON, and P. J. BALINT-KURTI et al., 1997  Characterization of the tomato Cf-4 gene for resistance to Cladosporium fulvum identifies sequences that determine recognitional specificity in Cf-4 and Cf-9. Plant Cell 9:2209-2224.[Abstract]

TIAN, D., H. ARAKI, E. A. STAHL, J. BERGELSON, and M. KREITMAN, 2002  Signature of balancing selection in Arabidopsis. Proc. Natl. Acad. Sci. USA 99:11525-11530.[Abstract/Free Full Text]

WHALEN, M. C., R. W. INNES, A. F. BENT, and B. J. STASKAWICZ, 1991  Identification of Pseudomonas syringae pathogens of Arabidopsis and a bacterial locus determining avirulence on both Arabidopsis and soybean. Plant Cell 3:49-59.[Abstract/Free Full Text]

YU, G.-L., F. KATAGIRI, and F. M. AUSUBEL, 1993  Arabidopsis mutations at the RPS2 locus result in loss of resistance to Pseudomonas syringae strains expressing the avirulence gene avrRpt2.. Mol. Plant Microbe Interact. 6:434-443.[Medline]




This article has been cited by other articles:


Home page
 GeneticsHome page
L. Huang, S. Brooks, W. Li, J. Fellers, J. C. Nelson, and B. Gill
Evolution of New Disease Specificity at a Simple Resistance Locus in a Crop-Weed Complex: Reconstitution of the Lr21 Gene in Wheat
Genetics, June 1, 2009; 182(2): 595 - 602.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
K. S. Caldwell and R. W. Michelmore
Arabidopsis thaliana Genes Encoding Defense Signaling and Recognition Proteins Exhibit Contrasting Evolutionary Dynamics
Genetics, February 1, 2009; 181(2): 671 - 684.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
V. Geffroy, C. Macadre, P. David, A. Pedrosa-Harand, M. Sevignac, C. Dauga, and T. Langin
Molecular Analysis of a Large Subtelomeric Nucleotide-Binding-Site-Leucine-Rich-Repeat Family in Two Representative Genotypes of the Major Gene Pools of Phaseolus vulgaris
Genetics, February 1, 2009; 181(2): 405 - 419.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
I. Ashikawa, N. Hayashi, H. Yamane, H. Kanamori, J. Wu, T. Matsumoto, K. Ono, and M. Yano
Two Adjacent Nucleotide-Binding Site-Leucine-Rich Repeat Class Genes Are Required to Confer Pikm-Specific Rice Blast Resistance
Genetics, December 1, 2008; 180(4): 2267 - 2276.
[Abstract] [Full Text] [PDF]

Home page
 Am. J. Bot.Home page
A. L. Caicedo
Geographic diversity cline of R gene homologs in wild populations of Solanum pimpinellifolium (Solanaceae)
Am. J. Botany, March 1, 2008; 95(3): 393 - 398.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
H. Kuang, H. J. van Eck, D. Sicard, R. Michelmore, and E. Nevo
Evolution and Genetic Population Structure of Prickly Lettuce (Lactuca serriola) and Its RGC2 Resistance Gene Cluster
Genetics, March 1, 2008; 178(3): 1547 - 1558.
[Abstract] [Full Text] [PDF]

Home page
 Plant CellHome page
R. M.P. Van Poecke, M. Sato, L. Lenarz-Wyatt, S. Weisberg, and F. Katagiri
Natural Variation in RPS2-Mediated Resistance among Arabidopsis Accessions: Correlation between Gene Expression Profiles and Phenotypic Responses
PLANT CELL, December 1, 2007; 19(12): 4046 - 4060.
[Abstract] [Full Text] [PDF]

Home page
 Proc R Soc BHome page
C.J.E Metcalf and D.N Koons
Environmental uncertainty, autocorrelation and the evolution of survival
Proc R Soc B, September 7, 2007; 274(1622): 2153 - 2160.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
U. Orgil, H. Araki, S. Tangchaiburana, R. Berkey, and S. Xiao
Intraspecific Genetic Variations, Fitness Cost and Benefit of RPW8, A Disease Resistance Locus in Arabidopsis thaliana
Genetics, August 1, 2007; 176(4): 2317 - 2333.
[Abstract] [Full Text] [PDF]

Home page
 Proc. Natl. Acad. Sci. USAHome page
J. O. Borevitz, S. P. Hazen, T. P. Michael, G. P. Morris, I. R. Baxter, T. T. Hu, H. Chen, J. D. Werner, M. Nordborg, D. E. Salt, et al.
Genome-wide patterns of single-feature polymorphism in Arabidopsis thaliana
PNAS, July 17, 2007; 104(29): 12057 - 12062.
[Abstract] [Full Text] [PDF]

Home page
 Nucleic Acids ResHome page
A. Boyko, P. Kathiria, F. J. Zemp, Y. Yao, I. Pogribny, and I. Kovalchuk
Transgenerational changes in the genome stability and methylation in pathogen-infected plants: (Virus-induced plant genome instability)
Nucleic Acids Res., March 12, 2007; 35(5): 1714 - 1725.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
L. E. Rose, R. W. Michelmore, and C. H. Langley
Natural Variation in the Pto Disease Resistance Gene Within Species of Wild Tomato (Lycopersicon). II. Population Genetics of Pto
Genetics, March 1, 2007; 175(3): 1307 - 1319.
[Abstract] [Full Text] [PDF]

Home page
 Mol Biol EvolHome page
D. S. Guttman, S. J. Gropp, R. L. Morgan, and P. W. Wang
Diversifying Selection Drives the Evolution of the Type III Secretion System Pilus of Pseudomonas syringae
Mol. Biol. Evol., December 1, 2006; 23(12): 2342 - 2354.
[Abstract] [Full Text] [PDF]

Home page
 Plant CellHome page
E. G. Bakker, C. Toomajian, M. Kreitman, and J. Bergelson
A Genome-Wide Survey of R Gene Polymorphisms in Arabidopsis
PLANT CELL, August 1, 2006; 18(8): 1803 - 1818.
[Abstract] [Full Text] [PDF]

Home page
 Am. J. Bot.Home page
I. M. Ehrenreich and M. D. Purugganan
The molecular genetic basis of plant adaptation
Am. J. Botany, July 1, 2006; 93(7): 953 - 962.
[Abstract] [Full Text] [PDF]

Home page
 Proc. Natl. Acad. Sci. USAHome page
P. N. Dodds, G. J. Lawrence, A.-M. Catanzariti, T. Teh, C.-I. A. Wang, M. A. Ayliffe, B. Kobe, and J. G. Ellis
From the Cover: Direct protein interaction underlies gene-for-gene specificity and coevolution of the flax resistance genes and flax rust avirulence genes
PNAS, June 6, 2006; 103(23): 8888 - 8893.
[Abstract] [Full Text] [PDF]

Home page
 Proc R Soc BHome page
M. Konior, J. Radwan, M. Kolodziejczyk, and L. Keller
Strong association between a single gene and fertilization efficiency of males and fecundity of their mates in the bulb mite
Proc R Soc B, February 7, 2006; 273(1584): 309 - 314.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
D. J. Kliebenstein, M. A. L. West, H. van Leeuwen, K. Kim, R. W. Doerge, R. W. Michelmore, and D. A. St. Clair
Genomic Survey of Gene Expression Diversity in Arabidopsis thaliana
Genetics, February 1, 2006; 172(2): 1179 - 1189.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
J. Shen, H. Araki, L. Chen, J.-Q. Chen, and D. Tian
Unique Evolutionary Mechanism in R-Genes Under the Presence/Absence Polymorphism in Arabidopsis thaliana
Genetics, February 1, 2006; 172(2): 1243 - 1250.
[Abstract] [Full Text] [PDF]

Home page
 Mol Biol EvolHome page
D. A. Moeller and P. Tiffin
Genetic Diversity and the Evolutionary History of Plant Immunity Genes in Two Species of Zea
Mol. Biol. Evol., December 1, 2005; 22(12): 2480 - 2490.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
L. E. Rose, C. H. Langley, A. J. Bernal, and R. W. Michelmore
Natural Variation in the Pto Pathogen Resistance Gene Within Species of Wild Tomato (Lycopersicon). I. Functional Analysis of Pto Alleles
Genetics, September 1, 2005; 171(1): 345 - 357.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
J. M. Cork and M. D. Purugganan
High-Diversity Genes in the Arabidopsis Genome
Genetics, August 1, 2005; 170(4): 1897 - 1911.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
A. Takahashi and T. Takano-Shimizu
A High-Frequency Null Mutant of an Odorant-Binding Protein Gene, Obp57e, in Drosophila melanogaster
Genetics, June 1, 2005; 170(2): 709 - 718.
[Abstract] [Full Text] [PDF]

Home page
 Appl. Environ. Microbiol.Home page
G. Wichmann, D. Ritchie, C. S. Kousik, and J. Bergelson
Reduced Genetic Variation Occurs among Genes of the Highly Clonal Plant Pathogen Xanthomonas axonopodis pv. vesicatoria, Including the Effector Gene avrBs2
Appl. Envir. Microbiol., May 1, 2005; 71(5): 2418 - 2432.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
K. J. Schmid, S. Ramos-Onsins, H. Ringys-Beckstein, B. Weisshaar, and T. Mitchell-Olds
A Multilocus Sequence Survey in Arabidopsis thaliana Reveals a Genome-Wide Departure From a Neutral Model of DNA Sequence Polymorphism
Genetics, March 1, 2005; 169(3): 1601 - 1615.
[Abstract] [Full Text] [PDF]

Home page
 Mol Biol EvolHome page
S. I. Wright and B. S. Gaut
Molecular Population Genetics and the Search for Adaptive Evolution in Plants
Mol. Biol. Evol., March 1, 2005; 22(3): 506 - 519.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
J. G. Bishop, D. R. Ripoll, S. Bashir, C. M. B. Damasceno, J. D. Seeds, and J. K. C. Rose
Selection on Glycine {beta}-1,3-Endoglucanase Genes Differentially Inhibited by a Phytophthora Glucanase Inhibitor Protein
Genetics, February 1, 2005; 169(2): 1009 - 1019.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
E. M. Goss, M. Kreitman, and J. Bergelson
Genetic Diversity, Recombination and Cryptic Clades in Pseudomonas viridiflava Infecting Natural Populations of Arabidopsis thaliana
Genetics, January 1, 2005; 169(1): 21 - 35.
[Abstract] [Full Text] [PDF]

Home page
 Proc. Natl. Acad. Sci. USAHome page
A. L. Caicedo and B. A. Schaal
Heterogeneous evolutionary processes affect R gene diversity in natural populations of Solanum pimpinellifolium
PNAS, December 14, 2004; 101(50): 17444 - 17449.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
P. Tiffin, R. Hacker, and B. S. Gaut
Population Genetic Evidence for Rapid Changes in Intraspecific Diversity and Allelic Cycling of a Specialist Defense Gene in Zea
Genetics, September 1, 2004; 168(1): 425 - 434.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
P. Tiffin
Comparative Evolutionary Histories of Chitinase Genes in the Genus Zea and Family Poaceae
Genetics, July 1, 2004; 167(3): 1331 - 1340.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
M. J. Clauss and T. Mitchell-Olds
Functional Divergence in Tandemly Duplicated Arabidopsis thaliana Trypsin Inhibitor Genes
Genetics, March 1, 2004; 166(3): 1419 - 1436.
[Abstract] [Full Text] [PDF]

Home page
 GeneticsHome page
L. E. Rose, P. D. Bittner-Eddy, C. H. Langley, E. B. Holub, R. W. Michelmore, and J. L. Beynon
The Maintenance of Extreme Amino Acid Diversity at the Disease Resistance Gene, RPP13, in Arabidopsis thaliana
Genetics, March 1, 2004; 166(3): 1517 - 1527.
[Abstract] [Full Text] [PDF]