Signature of balancing selection in Arabidopsis

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

In common bean, the B4 disease resistance gene cluster is a complex cluster localized at the end of linkage group (LG) B4, containing at least three R specificities to the fungus Colletotrichum lindemuthianum. To investigate the evolution of this R cluster since the divergence of Andean and Mesoamerican gene pools, DNA sequences were characterized from two representative genotypes of the two major gene pools of common bean (BAT93: Mesoamerican; JaloEEP558: Andean). Sequences encoding 29 B4-CC nucleotide-binding-site-leucine-rich-repeat (B4-CNL) genes were determined-12 from JaloEEP558 and 17 from BAT93. Although sequence exchange events were identified, phylogenetic analyses revealed that they were not frequent enough to lead to homogenization of B4-CNL sequences within a haplotype. Genetic mapping based on pulsed-field gel electrophoresis separation confirmed that the B4-CNL family is a large family specific to one end of LG B4 and is present at two distinct blocks separated by 26 cM. Fluorescent in situ hybridization on meiotic pachytene chromosomes revealed that two B4-CNL blocks are located in the subtelomeric region of the short arm of chromosome 4 on both sides of a heterochromatic block (knob), suggesting that this peculiar genomic environment may favor the proliferation of a large R gene cluster.

Arabidopsis thaliana genes encoding defense signaling and recognition proteins exhibit contrasting evolutionary dynamics

The interplay between pathogen effectors, their host targets, and cognate recognition proteins provides various opportunities for antagonistic cycles of selection acting on plant and pathogen to achieve or abrogate resistance, respectively. Selection has previously been shown to maintain diversity in plant proteins involved in pathogen recognition and some of their cognate pathogen effectors. We analyzed the signatures of selection on 10 Arabidopsis thaliana genes encoding defense signal transduction proteins in plants, which are potential targets of pathogen effectors. There was insufficient evidence to reject neutral evolution for 6 genes encoding signaling components consistent with these proteins not being targets of effectors and/or indicative of constraints on their ability to coevolve with pathogen effectors. Functional constraints on effector targets may have provided the driving selective force for the evolution of guard proteins. PBS1, a known target of an effector, showed little variation but is known to be monitored by a variable guard protein. Evidence of selection maintaining diversity was present at NPR1, PAD4, and EDS1. Differences in the signatures of selection observed may reflect the numbers of effectors that target a particular protein, the presence or absence of a cognate guard protein, as well as functional constraints imposed by biochemical activities or interactions with plant proteins.

Haplotype diversity at the Pi-ta locus in cultivated rice and its wild relatives

The Pi-ta gene in rice confers resistance to races of Magnaporthe oryzae that contain AVR-Pita. Pi-ta encodes a predicted cytoplasmic receptor protein with a nucleotide-binding site and a leucine-rich domain. A panel of 51 Oryza accessions of AA genome species Oryza sativa, O. glaberrima, O. rufipogon, O. nivara, and O. barthii, and CC genome species O. officinalis were sequenced to investigate the diversity present in the exon and intron regions of the Pi-ta gene. Two major clades were identified, consisting of 16 different sequences with numerous insertion and deletions. Only one Pi-ta resistance allele was identified despite DNA sequences revealing 16 Pi-ta variants. Most differences were identified in the intron region, and obvious selection of any motif was not observed in the coding region of Pi-ta variants. Reverse-transcription polymerase chain reaction analysis of seedlings revealed that all Pi-ta variants were expressed with or without pathogen inoculation. The 15 Pi-ta variants can be translated into nine proteins highly similar to the Pi-ta protein. Resistance to M. oryzae expressing AVR-Pita correlates with alanine and susceptibility correlates with serine at position 918 of Pi-ta in most accessions examined. These data confirm that a single amino acid controlling resistance specificity underlies the evolution of resistance of Pi-ta genes in rice.

The trouble with sliding windows and the selective pressure in BRCA1

Sliding-window analysis has widely been used to uncover synonymous (silent, d_S) and nonsynonymous (replacement, d_N) rate variation along the protein sequence and to detect regions of a protein under selective constraint (indicated by d_N<d_S) or positive selection (indicated by d_N>d_S). The approach compares two or more protein-coding genes and plots estimates dˆ_S and dˆ_N from each sliding window along the sequence. Here we demonstrate that the approach produces artifactual trends of synonymous and nonsynonymous rate variation, with greater variation in dˆ_S than in dˆ_N. Such trends are generated even if the true d_S and d_N are constant along the whole protein and different codons are evolving independently. Many published tests of negative and positive selection using sliding windows that we have examined appear to be invalid because they fail to correct for multiple testing. Instead, likelihood ratio tests provide a more rigorous framework for detecting signals of natural selection affecting protein evolution. We demonstrate that a previous finding that a particular region of the BRCA1 gene experienced a synonymous rate reduction driven by purifying selection is likely an artifact of the sliding window analysis. We evaluate various sliding-window analyses in molecular evolution, population genetics, and comparative genomics, and argue that the approach is not generally valid if it is not known a priori that a trend exists and if no correction for multiple testing is applied.

Conditional neutrality at two adjacent NBS-LRR disease resistance loci in natural populations of Arabidopsis lyrata

We examined patterns of nucleotide diversity at a genomic region containing two linked candidate disease resistance (NBS-LRR) genes in seven populations of the outcrossing plant Arabidopsis lyrata. In comparison with two adjacent control genes and neutral reference genes across the genome, the NBS-LRR genes exhibited elevated nonsynonymous variation and a large number of major-effect polymorphisms causing early stop codons and/or frameshift mutations. In contrast, analysis of synonymous diversity provided no evidence that the region was subject to long-term balancing selection or recent selective sweeps in any of the seven populations surveyed. Also in contrast with earlier surveys of one of these R genes, there was no evidence that the resistance genes or the major-effect mutations were subject to elevated differentiation between populations. We suggest that conditional neutrality in the absence of the corresponding pathogen, rather than long-term balancing selection or local adaptation, may in some circumstances be a significant cause of elevated functional polymorphism at R genes. In contrast with the R genes, analysis of diversity and differentiation at the flanking FERONIA locus showed high population divergence, suggesting local adaptation on this locus controlling male-female signalling during fertilization.

Grouped nucleotide polymorphism: a major contributor to genetic variation in Arabidopsis

While information for single nucleotide polymorphism is accumulating in many organisms, little is known about the magnitude and the occurrence of nucleotide dimorphism or higher order structural polymorphisms (grouped nucleotide polymorphisms or GNPs). To address these questions, we systematically investigated the genetic variations of 996 loci in 96 Arabidopsis accessions. Our data suggest that GNP loci are highly frequent in the genomes (40.4% of 996 loci), and that 66.7% of genetic variation in 996 loci is attributed to GNPs. The frequency distribution and the linkage analysis of GNP loci demonstrate that GNPs occurred randomly, locally, ceaselessly and independently in general, although they are particularly abundant in the region near centromeres. The age distribution of GNP loci shows that the majority of these loci are under a transient phase of neutral evolution. The characteristics of GNPs imply that a molecular isolation exists in GNP loci, and most likely the indels caused isolation, which could explain how GNPs are generated and maintained.

Geographic variation in adaptation at the molecular level: a case study of plant immunity genes

Natural selection imposed by interacting species frequently varies among geographic locations and can lead to local adaptation, where alternative phenotypes are found in different populations. Little is known, however, about whether geographically variable selection acting on traits that mediate species interactions is consistent or strong enough to influence patterns of nucleotide variation at individual loci. To investigate this question, we examined patterns of nucleotide diversity and population structure at 16 plant innate immunity genes, with putative functions in defending plants against pathogens or herbivores, from six populations of teosinte (Zea mays ssp. parviglumis). Specifically, we tested whether patterns of population structure and within-population diversity at immunity genes differed from patterns found at nonimmunity (reference) loci and from neutral expectations derived from coalescent simulations of structured populations. For the majority of genes, we detected no strong evidence of geographically variable selection. However, in the wound-induced serine protease inhibitor (wip1), which inhibits the hydrolysis of dietary proteins in insect herbivores, one population showed unusually high levels of genetic differentiation, very low levels of nucleotide polymorphism, and was fixed for a novel replacement substitution in the active site of the protein. Taken together, these data suggest that wip1 experienced a recent selective sweep in one geographic region; this pattern may reflect local adaptation or an ongoing species-wide sweep. Overall, our results indicate that a signature of local adaptation at the molecular level may be uncommon-particularly for traits that are under complex genetic control.

Evidence on the molecular basis of the Ac/ac adaptive cyanogenesis polymorphism in white clover (Trifolium repens L)

White clover is polymorphic for cyanogenesis, with both cyanogenic and acyanogenic plants occurring in nature. This chemical defense polymorphism is one of the longest-studied and best-documented examples of an adaptive polymorphism in plants. It is controlled by two independently segregating genes: Ac/ac controls the presence/absence of cyanogenic glucosides; and Li/li controls the presence/absence of their hydrolyzing enzyme, linamarase. Whereas Li is well characterized at the molecular level, Ac has remained unidentified. Here we report evidence that Ac corresponds to a gene encoding a cytochrome P450 of the CYP79D protein subfamily (CYP79D15), and we describe the apparent molecular basis of the Ac/ac polymorphism. CYP79D orthologs catalyze the first step in cyanogenic glucoside biosynthesis in other cyanogenic plant species. In white clover, Southern hybridizations indicate that CYP79D15 occurs as a single-copy gene in cyanogenic plants but is absent from the genomes of ac plants. Gene-expression analyses by RT-PCR corroborate this finding. This apparent molecular basis of the Ac/ac polymorphism parallels our previous findings for the Li/li polymorphism, which also arises through the presence/absence of a single-copy gene. The nature of these polymorphisms may reflect white clover's evolutionary origin as an allotetraploid derived from cyanogenic and acyanogenic diploid progenitors.

Genomic organization, rapid evolution and meiotic instability of nucleotide-binding-site-encoding genes in a new fruit crop, "chestnut rose"

From chestnut rose, a promising fruit crop of the Rosa genus, powdery mildew disease-resistant and susceptible genotypes and their F(1) progeny were used to isolate nucleotide-binding-site (NBS)-encoding genes using 19 degenerate primer pairs and an additional cloning method called overlapping extension amplification. A total of 126 genes were harvested; of these, 38 were from a resistant parent, 37 from a susceptible parent, and 51 from F(1) progeny. A phylogenetic tree was constructed, which revealed that NBS sequences from parents and F(1) progeny tend to form a mixture and are well distributed among the branches of the tree. Mapping of these NBS genes suggested that their organization in the genome is a "tandem duplicated cluster" and, to a lesser extent, a "heterogeneous cluster." Intraspecific polymorphisms and interspecific divergence were detected by Southern blotting with NBS-encoding genes as probes. Sequencing on the nucleotide level revealed even more intraspecific variation: for the R4 gene, 9.81% of the nucleotides are polymorphic. Amino acid sites under positive selection were detected in the NBS region. Some NBS-encoding genes were meiotically unstable, which may due to recombination and deletion events. Moreover, a transposon-like element was isolated in the flanking region of NBS genes, implying a possible role for transposon in the evolutionary history of resistance genes.

Evolution and genetic population structure of prickly lettuce (Lactuca serriola) and its RGC2 resistance gene cluster

Genetic structure and diversity of natural populations of prickly lettuce (Lactuca serriola) were studied using AFLP markers and then compared with the diversity of the RGC2 disease resistance gene cluster. Screening of 696 accessions from 41 populations using 319 AFLP markers showed that eastern Turkish and Armenian populations were the most diverse populations and might be located in the origin and center of diversity of L. serriola. Screening 709 accessions using the microsatellite MSATE6 that is located in the coding region of most RGC2 homologs detected 366 different haplotypes. Again, the eastern Turkish and Armenian populations had the highest diversities at the RGC2 cluster. The diversities at the RGC2 cluster in different populations were significantly correlated with their genomewide diversities. There was significant variation of copy number of RGC2 homologs in different populations, ranging from 12 to 22 copies per genome. The nucleotide diversities of two conserved lineages (type II) of RGC2 genes (K and L) were not correlated with diversities calculated using the MSATE6 or AFLP data. We hypothesize that the high genomewide diversity and diversity of the RGC2 cluster in eastern Turkish and Armenian populations resulted from high abiotic and biotic stresses in the regions of origin of L. serriola.

Low levels of polymorphism in genes that control the activation of defense response in Arabidopsis thaliana

Plants use signaling pathways involving salicylic acid, jasmonic acid, and ethylene to defend against pathogen and herbivore attack. Many defense response genes involved in these signaling pathways have been characterized, but little is known about the selective pressures they experience. A representative set of 27 defense response genes were resequenced in a worldwide set of 96 Arabidopsis thaliana accessions, and patterns of single nucleotide polymorphisms (SNPs) were evaluated in relation to an empirical distribution of SNPs generated from either 876 fragments or 236 fragments with >400 bp coding sequence (this latter set was selected for comparisons with coding sequences) distributed across the genomes of the same set of accessions. Defense response genes have significantly fewer protein variants, display lower levels of nonsynonymous nucleotide diversity, and have fewer nonsynonymous segregating sites. The majority of defense response genes appear to be experiencing purifying selection, given the dearth of protein variation in this set of genes. Eight genes exhibit some evidence of partial selective sweeps or transient balancing selection. These results therefore provide a strong contrast to the high levels of balancing selection exhibited by genes at the upstream positions in these signaling pathways.

Widespread genetic incompatibility in C. elegans maintained by balancing selection

Natural selection is expected to eliminate genetic incompatibilities from interbreeding populations. We have discovered a globally distributed incompatibility in the primarily selfing species Caenorhabditis elegans that has been maintained despite its negative consequences for fitness. Embryos homozygous for a naturally occurring deletion of the zygotically acting gene zeel-1 arrest if their sperm parent carries an incompatible allele of a second, paternal-effect locus, peel-1. The two interacting loci are tightly linked, with incompatible alleles occurring in linkage disequilibrium in two common haplotypes. These haplotypes exhibit elevated sequence divergence, and population genetic analyses of this region indicate that natural selection is preserving both haplotypes in the population. Our data suggest that long-term maintenance of a balanced polymorphism has permitted the incompatibility to persist despite gene flow across the rest of the genome.

The relationship of host-mediated induced resistance to polymorphism in gene-for-gene relationships

Gene-for-gene relationships are a common feature of plant-parasite interactions. Polymorphism at host resistance and parasite avirulence loci is maintained if there is negative, direct frequency-dependent selection on alleles of either gene. More specifically, selection of this kind is generated when the disease is polycyclic with frequent auto-infection. When an incompatible interaction occurs between a resistant host and an avirulent parasite, systemic defenses are triggered, rendering the plant more resistant to a later attack by another parasite. However, induced resistance (IR) incurs a fitness cost to the plant. Here, the effect of IR on polymorphism in gene-for-gene interactions is investigated. First, in an infinite population model in which parasites have two generations per host generation, increasing the fitness cost of IR increases selection for susceptible plants at low disease severity, while increasing the effectiveness of IR against further parasite attacks enhances selection for resistant plants at high disease severity. This reduces the possibility of polymorphism being maintained in host and parasite populations. In finite population models, the number of plants varies over time as a function of the disease burden of the population. Polymorphism in gene-for-gene relationships is then more stable at high disease prevalence and severity if IR reactions are more costly when there is competition for resources between plants.

Characterization of natural and induced variation in the LOV1 gene, a CC-NB-LRR gene conferring victorin sensitivity and disease susceptibility in Arabidopsis

The fungus Cochliobolus victoriae, the causal agent of Victoria blight, produces a compound called victorin that is required for pathogenicity of the fungus. Victorin alone reproduces disease symptoms on sensitive plants. Victorin sensitivity and susceptibility to C. victoriae were originally described on oats but have since been identified on Arabidopsis thaliana. Victorin sensitivity and disease susceptibility in Arabidopsis are conferred by LOV1, a coiled-coil-nucleotide-binding-leucine-rich repeat (CC-NB-LRR) protein. We sequenced the LOV1 gene from 59 victorin-insensitive mutants and found that the spectrum of mutations causing LOV1 loss of function was similar to that found to cause loss of function of RPM1, a CC-NB-LRR resistance protein. Also, many of the mutated residues in LOV1 are in conserved motifs required for resistance protein function. These data indicate that LOV1 may have a mechanism of action similar to resistance proteins. Victorin sensitivity was found to be the prevalent phenotype in a survey of 30 Arabidopsis ecotypes, and we found very little genetic variation among LOV1 alleles. As selection would not be expected to preserve a functional LOV1 gene to confer victorin sensitivity and disease susceptibility, we propose that LOV1 may function as a resistance gene to a naturally-occurring pathogen of Arabidopsis.

The blast resistance gene Pi37 encodes a nucleotide binding site leucine-rich repeat protein and is a member of a resistance gene cluster on rice chromosome 1

The resistance (R) gene Pi37, present in the rice cultivar St. No. 1, was isolated by an in silico map-based cloning procedure. The equivalent genetic region in Nipponbare contains four nucleotide binding site-leucine-rich repeat (NBS-LRR) type loci. These four candidates for Pi37 (Pi37-1, -2, -3, and -4) were amplified separately from St. No. 1 via long-range PCR, and cloned into a binary vector. Each construct was individually transformed into the highly blast susceptible cultivar Q1063. The subsequent complementation analysis revealed Pi37-3 to be the functional gene, while -1, -2, and -4 are probably pseudogenes. Pi37 encodes a 1290 peptide NBS-LRR product, and the presence of substitutions at two sites in the NBS region (V239A and I247M) is associated with the resistance phenotype. Semiquantitative expression analysis showed that in St. No. 1, Pi37 was constitutively expressed and only slightly induced by blast infection. Transient expression experiments indicated that the Pi37 product is restricted to the cytoplasm. Pi37-3 is thought to have evolved recently from -2, which in turn was derived from an ancestral -1 sequence. Pi37-4 is likely the most recently evolved member of the cluster and probably represents a duplication of -3. The four Pi37 paralogs are more closely related to maize rp1 than to any of the currently isolated rice blast R genes Pita, Pib, Pi9, Pi2, Piz-t, and Pi36.

Unique nucleotide polymorphism of ankyrin gene cluster in Arabidopsis

The ankyrin (ANK) gene cluster is a part of a multigene family encoding ANK transmembrane proteins in Arabidopsis thaliana, and plays an important role in protein-protein interactions and in signal pathways. In contrast to other regions of a genome, the ANK gene cluster exhibits an extremely high level of DNA polymorphism in an approximately 5-kb region, without apparent decay. Phylogenetic analysis detects two clear, deeply differentiated haplotypes (dimorphism). The divergence between haplotypes of accession Col-0 and Ler-0 (Hap-C and Hap-L) is estimated to be 10.7%, approximately equal to the 10.5% average divergence between A. thaliana and A. lyrata. Sequence comparisons for the ANK gene cluster homologues in Col-0 indicate that the members evolve independently, and that the similarity among paralogues is lower than between alleles. Very little intralocus recombination or gene conversion is detected in ANK regions. All these characteristics of the ANK gene cluster are consistent with a tandem gene duplication and birth-and-death process. The possible mechanisms for and implications of this elevated nucleotide variation are also discussed, including the suggestion of balancing selection.

Molecular evolution of the Li/li chemical defence polymorphism in white clover (Trifolium repens L.)

White clover (Trifolium repens) is naturally polymorphic for cyanogenesis (hydrogen cyanide release following tissue damage). The ecological factors favouring cyanogenic and acyanogenic plants have been examined in numerous studies over the last half century, making this one of the best-documented examples of an adaptive polymorphism in plants. White clover cyanogenesis is controlled by two, independently segregating Mendelian genes: Ac/ac controls the presence/absence of cyanogenic glucosides; and Li/li controls the presence/absence of their hydrolysing enzyme, linamarase. In this study, we examine the molecular evolution and population genetics of Li as it relates to the cyanogenesis polymorphism. We report here that Li exists as a single-copy gene in plants possessing linamarase activity, and that the absence of enzyme activity in li/li plants is correlated with the absence of much or all of the gene from the white clover genome. Consistent with this finding, we confirm by reverse transcription-polymerase chain reaction that Li gene expression is absent in plants lacking enzyme activity. In a molecular population genetic analysis of Li and three unlinked genes using a worldwide sample of clover plants, we find an absence of nucleotide variation and statistically significant deviations from neutrality at Li; these findings are consistent with recent positive directional selection at this cyanogenesis locus.

Molecular evolution of pathogenicity-island genes in Pseudomonas viridiflava

The bacterial pathogen Pseudomonas viridiflava possesses two pathogenicity islands (PAIs) that share many gene homologs, but are structurally and phenotypically differentiated (T-PAI and S-PAI). These PAIs are paralogous, but only one is present in each isolate. While this dual presence/absence polymorphism has been shown to be maintained by balancing selection, little is known about the molecular evolution of individual genes on the PAIs. Here we investigate genetic variation of 12 PAI gene loci (7 on T-PAI and 5 on S-PAI) in 96 worldwide isolates of P. viridiflava. These genes include avirulence genes (hopPsyA and avrE), their putative chaperones (shcA and avrF), and genes encoding the type III outer proteins (hrpA, hrpZ, and hrpW). Average nucleotide diversities in these genes (pi = 0.004-0.020) were close to those in the genetic background. Large numbers of recombination events were found within PAIs and a sign of positive selection was detected in avrE. These results suggest that the PAI genes are evolving relatively freely from each other on the PAIs, rather than as a single unit under balancing selection. Evolutionarily stable PAIs may be preferable in this species because preexisting genetic variation enables P. viridiflava to respond rapidly to natural selection.

Nucleotide diversity and linkage disequilibrium in 11 expressed resistance candidate genes in Lolium perenne

BACKGROUND

Association analysis is an alternative way for QTL mapping in ryegrass. So far, knowledge on nucleotide diversity and linkage disequilibrium in ryegrass is lacking, which is essential for the efficiency of association analyses.

RESULTS

11 expressed disease resistance candidate (R) genes including 6 nucleotide binding site and leucine rich repeat (NBS-LRR) like genes and 5 non-NBS-LRR genes were analyzed for nucleotide diversity. For each of the genes about 1 kb genomic fragments were isolated from 20 heterozygous genotypes in ryegrass. The number of haplotypes per gene ranged from 9 to 27. On average, one single nucleotide polymorphism (SNP) was present per 33 bp between two randomly sampled sequences for the 11 genes. NBS-LRR like gene fragments showed a high degree of nucleotide diversity, with one SNP every 22 bp between two randomly sampled sequences. NBS-LRR like gene fragments showed very high non-synonymous mutation rates, leading to altered amino acid sequences. Particularly LRR regions showed very high diversity with on average one SNP every 10 bp between two sequences. In contrast, non-NBS LRR resistance candidate genes showed a lower degree of nucleotide diversity, with one SNP every 112 bp. 78% of haplotypes occurred at low frequency (<5%) within the collection of 20 genotypes. Low intragenic LD was detected for most R genes, and rapid LD decay within 500 bp was detected.

CONCLUSION

Substantial LD decay was found within a distance of 500 bp for most resistance candidate genes in this study. Hence, LD based association analysis is feasible and promising for QTL fine mapping of resistance traits in ryegrass.

Common sequence polymorphisms shaping genetic diversity in Arabidopsis thaliana

The genomes of individuals from the same species vary in sequence as a result of different evolutionary processes. To examine the patterns of, and the forces shaping, sequence variation in Arabidopsis thaliana, we performed high-density array resequencing of 20 diverse strains (accessions). More than 1 million nonredundant single-nucleotide polymorphisms (SNPs) were identified at moderate false discovery rates (FDRs), and approximately 4% of the genome was identified as being highly dissimilar or deleted relative to the reference genome sequence. Patterns of polymorphism are highly nonrandom among gene families, with genes mediating interaction with the biotic environment having exceptional polymorphism levels. At the chromosomal scale, regional variation in polymorphism was readily apparent. A scan for recent selective sweeps revealed several candidate regions, including a notable example in which almost all variation was removed in a 500-kilobase window. Analyzing the polymorphisms we describe in larger sets of accessions will enable a detailed understanding of forces shaping population-wide sequence variation in A. thaliana.

Genome-wide patterns of single-feature polymorphism in Arabidopsis thaliana

We used hybridization to the ATH1 gene expression array to interrogate genomic DNA diversity in 23 wild strains (accessions) of Arabidopsis thaliana (arabidopsis), in comparison with the reference strain Columbia (Col). At <1% false discovery rate, we detected 77,420 single-feature polymorphisms (SFPs) with distinct patterns of variation across the genome. Total and pair-wise diversity was higher near the centromeres and the heterochromatic knob region, but overall diversity was positively correlated with recombination rate (R(2) = 3.1%). The difference between total and pair-wise SFP diversity is a relative measure contrasting diversifying or frequency-dependent selection, similar to Tajima's D, and can be calibrated by the empirical genome-wide distribution. Each unique locus, centered on a gene, has a diversity and selection score that suggest a relative role in past evolutionary processes. Homologs of disease resistance (R) genes include members with especially high levels of diversity often showing frequency-dependent selection and occasionally evidence of a past selective sweep. Receptor-like and S-locus proteins also contained members with elevated levels of diversity and signatures of selection, whereas other gene families, bHLH, F-box, and RING finger proteins, showed more typical levels of diversity. SFPs identified with the gene expression array also provide an empirical hybridization polymorphism background for studies of gene expression polymorphism and are available through the genome browser http://signal.salk.edu/cgi-bin/AtSFP.

Intraspecific genetic variations, fitness cost and benefit of RPW8, a disease resistance locus in Arabidopsis thaliana

The RPW8 locus of Arabidopsis thaliana confers broad-spectrum resistance to powdery mildew pathogens. In many A. thaliana accessions, this locus contains two homologous genes, RPW8.1 and RPW8.2. In some susceptible accessions, however, these two genes are replaced by HR4, a homolog of RPW8.1. Here, we show that RPW8.2 from A. lyrata conferred powdery mildew resistance in A. thaliana, suggesting that RPW8.2 might have gained the resistance function before the speciation of A. thaliana and A. lyrata. To investigate how RPW8 has been maintained in A. thaliana, we examined the nucleotide sequence polymorphisms in RPW8 from 51 A. thaliana accessions, related disease reaction phenotypes to the evolutionary history of RPW8.1 and RPW8.2, and identified mutations that confer phenotypic variations. The average nucleotide diversities were high at RPW8.1 and RPW8.2, showing no sign of selective sweep. Moreover, we found that expression of RPW8 incurs fitness benefits and costs on A. thaliana in the presence and absence of the pathogens, respectively. Our results suggest that polymorphisms at the RPW8 locus in A. thaliana may have been maintained by complex selective forces, including those from the fitness benefits and costs both associated with RPW8.

Natural genetic variation in Arabidopsis: tools, traits and prospects for evolutionary ecology

BACKGROUND

The model plant Arabidopsis thaliana (Arabidopsis) shows a wide range of genetic and trait variation among wild accessions. Because of its unparalleled biological and genomic resources, the potential of Arabidopsis for molecular genetic analysis of this natural variation has increased dramatically in recent years.

SCOPE

Advanced genomics has accelerated molecular phylogenetic analysis and gene identification by quantitative trait loci (QTL) mapping and/or association mapping in Arabidopsis. In particular, QTL mapping utilizing natural accessions is now becoming a major strategy of gene isolation, offering an alternative to artificial mutant lines. Furthermore, the genomic information is used by researchers to uncover the signature of natural selection acting on the genes that contribute to phenotypic variation. The evolutionary significance of such genes has been evaluated in traits such as disease resistance and flowering time. However, although molecular hallmarks of selection have been found for the genes in question, a corresponding ecological scenario of adaptive evolution has been difficult to prove. Ecological strategies, including reciprocal transplant experiments and competition experiments, and utilizing near-isogenic lines of alleles of interest will be a powerful tool to measure the relative fitness of phenotypic and/or allelic variants.

CONCLUSIONS

As the plant model organism, Arabidopsis provides a wealth of molecular background information for evolutionary genetics. Because genetic diversity between and within Arabidopsis populations is much higher than anticipated, combining this background information with ecological approaches might well establish Arabidopsis as a model organism for plant evolutionary ecology.

Fitness Effects Associated with the Major Flowering Time GeneFRIGIDAinArabidopsis thalianain the Field

To date, the effect of natural selection on candidate genes underlying complex traits has rarely been studied experimentally, especially under ecologically realistic conditions. Here we report that the effect of selection on the flowering time gene FRIGIDA (FRI) reverses depending on the season of germination and allelic variation at the interacting gene FLOWERING LOCUS C (FLC). In field studies of 136 European accessions of Arabidopsis thaliana, accessions with putatively functional FRI alleles had higher winter survival in one FLC background in a fall-germinating cohort, but accessions with deletion null FRI alleles had greater seed production in the other FLC background in a spring-germinating cohort. Consistent with FRI's role in flowering, selection analyses suggest that the difference in winter survival can be attributed to time to bolting. However, in the spring cohort, the fitness difference was associated with rosette size. Our analyses also reveal that controlling for population structure with estimates of inferred ancestry and a geographical restriction was essential for detecting fitness associations. Overall, our results suggest that the combined effects of seasonally varying selection and epistasis could explain the maintenance of variation at FRI and, more generally, may be important in the evolution of genes underlying complex traits.

Stability of genetic polymorphism in host-parasite interactions

Allelic diversity is common at host loci involved in parasite recognition, such as the major histocompatibility complex in vertebrates or gene-for-gene relationships in plants, and in corresponding loci encoding antigenic molecules in parasites. Diverse factors have been proposed in models to account for genetic polymorphism in host-parasite recognition. Here, a simple but general theory of host-parasite coevolution is developed. Coevolution implies the existence of indirect frequency-dependent selection (FDS), because natural selection on the host depends on the frequency of a parasite gene, and vice versa. It is shown that polymorphism can be maintained in both organisms only if there is negative, direct FDS, such that the strength of natural selection for the host resistance allele, the parasite virulence allele or both declines with increasing frequency of that allele itself. This condition may be fulfilled if the parasite has more than one generation in the same host individual, a feature which is common to most diseases. It is argued that the general theory encompasses almost all factors previously proposed to account for polymorphism at corresponding host and parasite loci, including those controlling gene-for-gene interactions.

Maize Sh2 gene is constrained by natural selection but escaped domestication

In Zea mays L., we studied the molecular evolution of Shrunken2 (Sh2), a gene that encodes the large subunits of a major enzyme in endosperm starch biosynthesis, ADP-glucose pyrophosphorylase. We compared 4669 bp of the Sh2 coding region on 50 accessions of maize and teosinte. Very few nucleotide polymorphisms were found when compared with other genes in Z. mays, revealing an effect of purifying selection in the whole species that predates domestication. Additionally, the comparison of Sh2 sequences in all Z. mays subspecies and outgroups Z. diploperennis and Tripsacum dactyloides suggests the occurrence of an ancient selective sweep in the Sh2 3' region. The amount and nature of nucleotide diversity are similar in both maize and teosinte, confirming previous results that suggested that Sh2 has not been involved in maize domestication. The very low level of nucleotide diversity as well as the highly conserved protein sequence suggest that natural selection retained effective Sh2 allele(s) long before agriculture started, making human selection inefficient on this gene.

Natural variation in the Pto disease resistance gene within species of wild tomato (Lycopersicon). II. Population genetics of Pto

Disease resistance to the bacterial pathogen Pseudomonas syringae pv. tomato (Pst) in the host species Lycopersicon esculentum, the cultivated tomato, and the closely related L. pimpinellifolium is triggered by the physical interaction between the protein products of the host resistance (R) gene Pto and the pathogen avirulence genes AvrPto and AvrPtoB. Sequence variation at the Pto locus was surveyed in natural populations of seven species of Lycopersicon to test hypotheses of host-parasite coevolution and functional adaptation of the Pto gene. Pto shows significantly higher nonsynonymous polymorphism than 14 other non-R-gene loci in the same samples of Lycopersicon species, while showing no difference in synonymous polymorphism, suggesting that the maintenance of amino acid polymorphism at this locus is mediated by pathogen selection. Also, a larger proportion of ancestral variation is maintained at Pto as compared to these non-R-gene loci. The frequency spectrum of amino acid polymorphisms known to negatively affect Pto function is skewed toward low frequency compared to amino acid polymorphisms that do not affect function or silent polymorphisms. Therefore, the evolution of Pto appears to be influenced by a mixture of both purifying and balancing selection.

A molecular isolation mechanism associated with high intra-specific diversity in rice

High levels of inter-specific diversity are expected due to genetic isolation, the reproductive or geographical barriers, which lead to the accumulation of nucleotide variation. However, high levels of genetic variation are repeatedly observed even within species, notably at loci of the human major histocompatability complex and of plant resistance genes. Are molecular isolations responsible for the high intra-specific variation? To address this issue, we performed a genome-wide survey of the relationship between the possible factors that could cause genetic isolation, and the level of polymorphism, based on two rice genome comparisons. Here, we show that the levels of polymorphism in rice genes are positively correlated with the proportions of non-alignable flanking sequences, and that the correlation is observed even in single-copy genes. The physical locations of the genes were also investigated, and a strong association between the asymmetric architecture of genomes and the levels of polymorphism was revealed. These results suggest that the flank heterogeneity and the asymmetric architecture between genomes serve as isolation mechanisms at the molecular level that result in accumulation of higher genetic variation. This mechanism is of fundamental importance to understand natural genetic variation within species.

Combining population genomics and quantitative genetics: finding the genes underlying ecologically important traits

A central challenge in evolutionary biology is to identify genes underlying ecologically important traits and describe the fitness consequences of naturally occurring variation at these loci. To address this goal, several novel approaches have been developed, including 'population genomics,' where a large number of molecular markers are scored in individuals from different environments with the goal of identifying markers showing unusual patterns of variation, potentially due to selection at linked sites. Such approaches are appealing because of (1) the increasing ease of generating large numbers of genetic markers, (2) the ability to scan the genome without measuring phenotypes and (3) the simplicity of sampling individuals without knowledge of their breeding history. Although such approaches are inherently applicable to non-model systems, to date these studies have been limited in their ability to uncover functionally relevant genes. By contrast, quantitative genetics has a rich history, and more recently, quantitative trait locus (QTL) mapping has had some success in identifying genes underlying ecologically relevant variation even in novel systems. QTL mapping, however, requires (1) genetic markers that specifically differentiate parental forms, (2) a focus on a particular measurable phenotype and (3) controlled breeding and maintenance of large numbers of progeny. Here we present current advances and suggest future directions that take advantage of population genomics and quantitative genetic approaches - in both model and non-model systems. Specifically, we discuss advantages and limitations of each method and argue that a combination of the two provides a powerful approach to uncovering the molecular mechanisms responsible for adaptation.

Transgenerational changes in the genome stability and methylation in pathogen-infected plants: (virus-induced plant genome instability)

Previously, we reported the generation of a virus-induced systemic signal that increased the somatic and meiotic recombination rates in tobacco mosaic virus (TMV)-infected tobacco plants. Here, we analyzed the progeny of plants that received the signal and found that these plants also have a higher frequency of rearrangements in the loci carrying the homology to LRR region of the gene of resistance to TMV (N-gene). Analysis of the stability of repetitive elements from Nicotiana tabacum loci and 5.8S ribosomal RNA loci did not show any changes. Further analysis of the changes in the progeny of infected plants revealed that they had substantially hypermethylated genomes. At the same time, loci-specific methylation analysis showed: (1) profound hypomethylation in several LRR-containing loci; (2) substantial hypermethylation of actin loci and (3) no change in methylation in the loci of repetitive elements from N. tabacum or 5.8S ribosomal RNA. Global genome hypermethylation of the progeny is believed to be part of a general protection mechanism against stress, whereas locus-specific hypomethylation is associated with a higher frequency of rearrangements. Increased recombination events combined with the specific methylation pattern induced by pathogen attack could be a sign of an adaptive response by plants.

NATURAL SELECTION ON A POLYMORPHIC DISEASE-RESISTANCE LOCUS INIPOMOEA PURPUREA

Although disease-resistance polymorphisms are common in natural plant populations, the mechanisms responsible for this variation are not well understood. Theoretical models predict that balancing selection can maintain polymorphism within a population if the fitness effects of a resistance allele vary from a net cost to a net benefit, depending upon the extent of pathogen damage. However, there have been a few attempts to determine how commonly this mechanism operates in natural plant-pathogen interactions. Ipomoea purpurea populations are often polymorphic for resistance and susceptibility alleles at a locus that influences resistance to the fungal pathogen, Coleosporium ipomoeae. We measured the fitness effects of resistance over three consecutive years at natural and manipulated levels of damage to characterize the type of selection acting on this locus. Costs of resistance varied in magnitude from undetectable to 15.5%, whereas benefits of resistance sometimes equaled, but never exceeded, these costs. In the absence of net benefits of resistance at natural or elevated levels of disease, we conclude that selection within individual populations of I. purpurea probably does not account completely for maintenance of this polymorphism. Rather, the persistence of this polymorphism is probably best explained by a combination of variable selection and meta-population processes.

Effects of introgression and recombination on haplotype structure and linkage disequilibrium surrounding a locus encoding Bymovirus resistance in barley

We present a detailed analysis of linkage disequilibrium (LD) in the physical and genetic context of the barley gene Hv-eIF4E, which confers resistance to the barley yellow mosaic virus (BYMV) complex. Eighty-three SNPs distributed over 132 kb of Hv-eIF4E and six additional fragments genetically mapped to its flanking region were used to derive haplotypes from 131 accessions. Three haplogroups were recognized, discriminating between the alleles rym4 and rym5, which each encode for a spectrum of resistance to BYMV. With increasing map distance, haplotypes of susceptible genotypes displayed diverse patterns driven mainly by recombination, whereas haplotype diversity within the subgroups of resistant genotypes was limited. We conclude that the breakdown of LD within 1 cM of the resistance gene was generated mainly by susceptible genotypes. Despite the LD decay, a significant association between haplotype and resistance to BYMV was detected up to a distance of 5.5 cM from the resistance gene. The LD pattern and the haplotype structure of the target chromosomal region are the result of interplay between low recombination and recent breeding history.

The role of pectate lyase and the jasmonic acid defense response in Pseudomonas viridiflava virulence

Pseudomonas viridiflava is a common pathogen of Arabidopsis thaliana in wild populations, yet very little is known about mechanisms of resistance and virulence in this interaction. We examined the induced defense response of A. thaliana to several strains of P. viridiflava collected from this host by quantifying the expression of PR-1 and LOX2/PDF1.2, which serve as markers for induction of the salicylic and jasmonic acid (JA) pathways, respectively. Growth of these strains then was assessed on Col-0, the fad3/7/8 and coil-1 mutants deficient in JA- and ethylene (ET)-induced defense responses, and the sid2-1 mutant deficient in salicylic acid-induced defense responses. All strains of P. viridiflava induced high expression of LOX2 and PDF1.2 on Col-0. In contrast, PR-1 expression was delayed and reduced relative to PDF1.2 expression. Additionally, three of four P. viridiflava strains were more virulent on fad3/7/8 relative to Col-0, whereas all strains were more virulent on coil-1 relative to Col-0, indicating that P. viridiflava generally may be suppressed by JA/ET-mediated defense responses. In contrast, no increase in the growth of P. viridiflava strains was observed in the sid2-1 mutant relative to Col-0. Parallel experiments were performed with the closely related P. syringae pv. tomato for comparative purposes. In addition, we assessed the role of pectate lyase and the alternative sigma factor HrpL in P. viridiflava virulence on A. thaliana and found that pectate lyase activity is correlated with virulence, whereas the removal of pectate lyase or HrpL significantly reduced virulence.

Methods for linkage disequilibrium mapping in crops

Linkage disequilibrium (LD) mapping in plants detects and locates quantitative trait loci (QTL) by the strength of the correlation between a trait and a marker. It offers greater precision in QTL location than family-based linkage analysis and should therefore lead to more efficient marker-assisted selection, facilitate gene discovery and help to meet the challenge of connecting sequence diversity with heritable phenotypic differences. Unlike family-based linkage analysis, LD mapping does not require family or pedigree information and can be applied to a range of experimental and non-experimental populations. However, care must be taken during analysis to control for the increased rate of false positive results arising from population structure and variety interrelationships. In this review, we discuss how suitable the recently developed alternative methods of LD mapping are for crops.

Transcriptome analysis of the salivary glands of Dermacentor andersoni Stiles (Acari: Ixodidae)

Amongst blood-feeding arthropods, ticks of the family Ixodidae (hard ticks) are vectors and reservoirs of a greater variety of infectious agents than any other ectoparasite. Salivary glands of ixodid ticks secrete a large number of pharmacologically active molecules that not only facilitate feeding but also promote establishment of infectious agents. Genomic, proteomic and immunologic characterization of bioactive salivary gland molecules are, therefore, important as they offer new insights into molecular events occurring at the tick-host interface and they have implications for development of novel control strategies. The present work uses complementary DNA (cDNA) sequence analysis to identify salivary gland transcripts expressed by the Rocky Mountain wood tick, Dermacentor andersoni, a vector of the human pathogens causing Rocky Mountain spotted fever, Colorado tick fever, tularemia, and Powassan encephalitis as well as the veterinary pathogen Anaplasma marginale. Dermacentor andersoni is also capable of inducing tick paralysis. Automated single-pass DNA sequencing was conducted on 1440 randomly selected cDNA clones from the salivary glands of adult female D. andersoni collected during the early stages of feeding (18-24h). Analysis of the expressed sequence tags (ESTs) resulted in 544 singletons and 218 clusters with more than one quality read and attempts were made to assign putative functions to tick genes based on amino acid identity to published protein databases. Approximately 25.6% (195) of the sequences showed limited or no homology to previously identified gene products. A number of novel sequences were identified which presented significant sequence similarity to mammalian genes normally associated with extracellular matrix (ECM), regulation of immune responses, tumor suppression, and wound healing. Several coding sequences possessed various degrees of homology to previously described proteins from other tick species. Preliminary nucleotide variation analysis of these and other tick sequences suggests extensive nucleotide diversity, which has implications for evolution of tick feeding. Intra-species diversity studies can be a promising tool for identifying sequence variations potentially associated with phenotypic traits affecting vector-host-pathogen interactions.

Elicitors, effectors, and R genes: the new paradigm and a lifetime supply of questions

The plant basal immune system can detect broadly present microbe-associated molecular patterns (MAMPs, also called PAMPs) and induce defenses, but adapted microbes express a suite of effector proteins that often act to suppress these defenses. Plants have evolved other receptors (R proteins) that detect these pathogen effectors and activate strong defenses. Pathogens can subsequently alter or delete their recognized effectors to avoid defense elicitation, at risk of a fitness cost associated with loss of those effectors. Significant research progress is revealing, among other things, mechanisms of MAMP perception, the host defense processes and specific host proteins that pathogen effectors target, the mechanisms of R protein activation, and the ways in which pathogen effector suites and R genes evolve. These findings carry practical ramifications for resistance durability and for future resistance engineering. The present review uses numerous questions to help clarify what we know and to identify areas that are ripe for further investigation.

Resources for Genetic Variation Studies

The rapid growth of genome-wide diversity databases, as well as ongoing large-scale resequencing projects targeting genes and other functional components of our genome, provide valuable resources of natural variation at the DNA sequence level. In this review, we briefly summarize the wealth of data on DNA polymorphisms in humans, the distribution of this diversity in the genome as well as among individuals, and the consequence of recombination on its organization. These data provide a set of powerful tools that can be used to better understand inherited phenotypic variation in humans. We discuss the implications for the design of studies investigating correlations between genotypes and phenotypes, both at the fundamental level of genome function and regulation, and for the mapping of disease genes.

MAMPs and MIMPs: proposed classifications for inducers of innate immunity

Plants encode a sophisticated innate immune system. Resistance against potential pathogens often relies on active responses. Prerequisite to the induction of defences is recognition of the pathogenic threat. Significant advances have been made in our understanding of the non-self molecules that are recognized by plants and the means by which plants perceive them. Established terms describing these recognition events, including microbe-associated molecular pattern (MAMP), MAMP-receptor, effector, and resistance (R) protein, need clarification to represent our current knowledge adequately. In this review we propose criteria to classify inducers of plant defence as either MAMPs or microbe-induced molecular patterns (MIMPs). We refine the definition of MAMP to mean a molecular sequence or structure in ANY pathogen-derived molecule that is perceived via direct interaction with a host defence receptor. MIMPs are modifications of host-derived molecules that are induced by an intrinsic activity of a pathogen-derived effector and are perceived by a host defence receptor. MAMP-receptors have previously been classified separately from R-proteins as a discrete class of surveillance molecules. However, MAMP-receptors and R-proteins cannot be distinguished on the basis of their protein structures or their induced responses. We propose that MAMP-receptors and MIMP-receptors are each a subset of R-proteins. Although our review is based on examples from plant pathogens and plants, the principles discussed might prove applicable to other organisms.

A genome-wide survey of R gene polymorphisms in Arabidopsis

We used polymorphism analysis to study the evolutionary dynamics of 27 disease resistance (R) genes by resequencing the leucine-rich repeat (LRR) region in 96 Arabidopsis thaliana accessions. We compared single nucleotide polymorphisms (SNPs) in these R genes to an empirical distribution of SNP in the same sample based on 876 fragments selected to sample the entire genome. LRR regions are highly polymorphic for protein variants but not for synonymous changes, suggesting that they generate many alleles maintained for short time periods. Recombination is also relatively common and important for generating protein variants. Although none of the genes is nearly as polymorphic as RPP13, a locus previously shown to have strong signatures of balancing selection, seven genes show weaker indications of balancing selection. Five R genes are relatively invariant, indicating young alleles, but all contain segregating protein variants. Polymorphism analysis in neighboring fragments yielded inconclusive evidence for recent selective sweeps at these loci. In addition, few alleles are candidates for rapid increases in frequency expected under directional selection. Haplotype sharing analysis revealed significant underrepresentation of R gene alleles with extended haplotypes compared with 1102 random genomic fragments. Lack of convincing evidence for directional selection or selective sweeps argues against an arms race driving R gene evolution. Instead, the data support transient or frequency-dependent selection maintaining protein variants at a locus for variable time periods.

The disease resistance gene Dm3 is infrequent in natural populations of Lactuca serriola due to deletions and frequent gene conversions at the RGC2 locus

Resistance genes can exhibit heterogeneous patterns of variation. However, there are few data on their frequency and variation in natural populations. We analysed the frequency and variation of the resistance gene Dm3, which confers resistance to Bremia lactucae (downy mildew) in 1033 accessions of Lactuca serriola (prickly lettuce) from 49 natural populations. Inoculations with an isolate of Bremia lactucae carrying avirulence gene Avr3 indicated that the frequency of Dm3 in natural populations of L. serriola was very low. Molecular analysis demonstrated that Dm3 was present in only one of the 1033 wild accessions analysed. The sequence of the 5' region of Dm3 was either highly conserved among accessions, or absent. In contrast, frequent chimeras were detected in the 3' leucine-rich repeat-encoding region. Therefore low frequency of the Dm3 specificity in natural populations was due to either the recent evolution of Dm3 specificity, or deletions of the whole gene as well as variation in 3' region caused by frequent gene conversions. This is the most extensive analysis of the prevalence of a known disease resistance gene to date, and indicates that the total number of resistance genes in a species may be very high. This has implications for the scales of germplasm conservation and exploitation of sources of resistance.

The molecular genetic basis of plant adaptation

How natural selection on adaptive traits is filtered to the genetic level remains largely unknown. Theory and quantitative trait locus (QTL) mapping have provided insights into the number and effect of genes underlying adaptations, but these results have been hampered by questions of applicability to real biological systems and poor resolution, respectively. Advances in molecular technologies have expedited the cloning of adaptive genes through both forward and reverse genetic approaches. Forward approaches start with adaptive traits and attempt to characterize their underlying genetic architectures through linkage disequilibrium mapping, QTL mapping, and other methods. Reverse screens search large sequence data sets for genes that possess the signature of selection. Though both approaches have been successful in identifying adaptive genes in plants, very few, if any, of these adaptations' molecular bases have been fully resolved. The continued isolation of plant adaptive genes will lead to a more comprehensive understanding of natural selection's effect on genes and genomes.

Distribution of genetic variation within and among local populations of Arabidopsis thaliana over its species range

A detailed description of local population structure in Arabidopsis thaliana is presented, including an assessment of the genetic relatedness of individuals collected from the same field. A hierarchical sample of four individuals from 37 local populations, including North America, England, Eastern and Western Europe, and Asia, and a selection of ecotypes, were analysed for variation in Adh, ChiA, FAH1, F3H, Rpm1, Rps5, and five microsatellite loci. Twenty-eight of the 37 population samples contained individuals with identical multilocus haplotypes, 12 of which were fixed for a single haplotype. These monomorphic populations were evenly distributed over the species range. Only in North America did we find a single multilocus haplotype shared among different populations, perhaps indicating a continental founder event. Despite the occurrence of local inbreeding, a considerable amount of genetic variation was found segregating within and among local populations. A novel analysis of haplotype differences reveals that genetic differentiation occurs at every geographic scale in A. thaliana, where we find a surprising under-representation of recent migrants between local populations. This leads us to hypothesize that most dispersal between A. thaliana populations is by pollen rather than seed. Based on the structure of A. thaliana populations, it appears that regional groups of local populations may provide the most appropriate genetic material for linkage disequilibrium mapping of adaptive traits.

Genetic mechanisms and evolutionary significance of natural variation in Arabidopsis

Genomic studies of natural variation in model organisms provide a bridge between molecular analyses of gene function and evolutionary investigations of adaptation and natural selection. In the model plant species Arabidopsis thaliana, recent studies of natural variation have led to the identification of genes underlying ecologically important complex traits, and provided new insights about the processes of genome evolution, geographic population structure, and the selective mechanisms shaping complex trait variation in natural populations. These advances illustrate the potential for a new synthesis to elucidate mechanisms for the adaptive evolution of complex traits from nucleotide sequences to real-world environments.

Presence/absence polymorphism for alternative pathogenicity islands in Pseudomonas viridiflava, a pathogen of Arabidopsis

The contribution of arms race dynamics to plant-pathogen coevolution has been called into question by the presence of balanced polymorphisms in resistance genes of Arabidopsis thaliana, but less is known about the pathogen side of the interaction. Here we investigate structural polymorphism in pathogenicity islands (PAIs) in Pseudomonas viridiflava, a prevalent bacterial pathogen of A. thaliana. PAIs encode the type III secretion system along with its effectors and are essential for pathogen recognition in plants. P. viridiflava harbors two structurally distinct and highly diverged PAI paralogs (T- and S-PAI) that are integrated in different chromosome locations in the P. viridiflava genome. Both PAIs are segregating as presence/absence polymorphisms such that only one PAI ([T-PAI, nablaS-PAI] and [nablaT-PAI, S-PAI]) is present in any individual cell. A worldwide population survey identified no isolate with neither or both PAI. T-PAI and S-PAI genotypes exhibit virulence differences and a host-specificity tradeoff. Orthologs of each PAI can be found in conserved syntenic locations in other Pseudomonas species, indicating vertical phylogenetic transmission in this genus. Molecular evolutionary analysis of PAI sequences also argues against "recent" horizontal transfer. Spikes in nucleotide divergence in flanking regions of PAI and nabla-PAI alleles suggest that the dual PAI polymorphism has been maintained in this species under some form of balancing selection. Virulence differences and host specificities are hypothesized to be responsible for the maintenance of the dual PAI system in this bacterial pathogen.

Geographical variation in selection, from phenotypes to molecules

Molecular technologies now allow researchers to isolate quantitative trait loci (QTLs) and measure patterns of gene sequence variation within chromosomal regions containing important polymorphisms. I develop a simulation model to investigate gene sequence evolution within genomic regions that harbor QTLs. The QTLs influence a trait experiencing geographical variation in selection, which is common in nature and produces obvious differentiation at the phenotypic level. Counter to expectations, the simulations suggest that selection can substantially affect quantitative genetic variation without altering the amount and pattern of molecular variation at sites closely linked to the QTLs. Even with large samples of gene sequences, the likelihood of rejecting neutrality is often low. The exception is situations where strong selection is combined with low migration among demes, conditions that may be common in many plant species. The results have implications for gene sequence surveys and, perhaps more generally, for interpreting the apparently weak connection between levels of molecular and quantitative trait variation within species.

Pervasive purifying selection characterizes the evolution of I2 homologs

We sampled 384 sequences related to the Solanum pimpinellifolium (=Lycopersicon pimpinellifolium) disease resistance (R) gene 12 from six species, potato, S. demissum, tomato, eggplant, pepper, and tobacco. These species represent increasing phylogenetic distance from potato to tobacco, within the family Solanaceae. Using sequence data from the nucleotide binding site (NBS) region of this gene, we tested models of gene family evolution and inferred patterns of selection acting on the NBS gene region and I2 gene family. We find that the I2 family has diversified within the family Solanaceae for at least 14 million years and evolves through a slow birth-and-death process requiring approximately 12 million years to homogenize gene copies within a species. Analyses of selection resolved a general pattern of strong purifying selection acting on individual codon positions within the NBS and on NBS lineages through time. Surprisingly, we find nine codon positions strongly affected by positive selection and six pairs of codon positions demonstrating correlated amino acid substitutions. Evolutionary analyses serve as bioinformatic tools with which to sort through the vast R gene diversity in plants and find candidates for new resistance specificities or to identify specific amino acid positions important for biochemical function. The slow birth-and-death evolution of I2 genes suggests that some NBS-leucine rich repeat-mediated resistances may not be overcome rapidly by virulence evolution and that the natural diversity of R genes is a potentially valuable source for durable resistance.

From DNA to fitness differences: sequences and structures of adaptive variants of Colias phosphoglucose isomerase (PGI)

Colias eurytheme butterflies display extensive allozyme polymorphism in the enzyme phosphoglucose isomerase (PGI). Earlier studies on biochemical and fitness effects of these genotypes found evidence of strong natural selection maintaining this polymorphism in the wild. Here we analyze the molecular features of this polymorphism by sequencing multiple alleles and modeling their structures. PGI is a dimer with rotational symmetry. Each monomer provides a critical residue to the other monomer's catalytic center. Sequenced alleles differ at multiple amino acid positions, including cryptic charge-neutral variation, but most consistent differences among the electromorph alleles are at the charge-changing amino acid sites. Principal candidate sites of selection, identified by structural and functional analyses and by their variants' population frequencies, occur in interpenetrating loops across the interface between monomers, where they may alter subunit interactions and catalytic center geometry. Comparison to a second (and basal) species, Colias meadii, also polymorphic for PGI under natural selection, reveals one fixed amino acid difference between their PGIs, which is located in the interpenetrating loop and accompanies functional differences among their variants. We also study nucleotide variability among the PGI alleles, comparing these data to similar data from another glycolytic enzyme gene, glyceraldehyde-3-phosphate dehydrogenase. Despite extensive nonsynonymous and synonymous polymorphism at PGI in each species, the only base changes fixed between species are the two causing the amino acid replacement; this absence of synonymous fixation yields a significant McDonald-Kreitman test. Analyses of these data suggest historical population expansion. Positive peaks of Tajima's D statistic, representing regions of neutral "hitchhiking," are found around the principal candidate sites of selection. This study provides novel views of molecular-structural mechanisms, and beginnings of historical evidence, for a long-persistent balanced enzyme polymorphism at PGI in these and perhaps other species.