Unusual signatures of highly adaptable R-loci in closely-related Arabidopsis species

Plant resistance genes (R-genes) evolve rapidly in response to changing environments. What are the most remarkable signatures of fast adaptive genes, besides the commonly revealed rapid divergence and high non-synonymous substitution rate? Here we investigated these changes in five R-loci following recent differentiation between Arabidopsis thaliana and Arabidopsis lyrata. Extreme differences in evolutionary rates were observed: e.g., an overall 5.46-9.83-fold different nucleotide diversity at two R-loci between species, ten-fold higher non-synonymous substitution rates within one species versus the other, significantly different Ka/Ks ratios between species for the same R-gene, and high interspecific divergence at one R-locus. Particularly, we observed an elevated level of trans-specific polymorphism at one R-locus and a differentially maintained presence/absence polymorphism at another. The high frequency of ancestral polymorphisms amongst R-genes suggests that the persistence of some functional variation is an important evolutionary mechanism shaping genetic variation in R-genes, while the variation of presence/absence polymorphisms provides a potential mechanism for malleable activation of adaptive resistance response pathways. The distinct patterns among R-genes suggest that the same R-gene ortholog can be quickly shaped by different evolutionary processes, e.g., purifying selection in one species but positive selection in a closely-related species.     

Genetic isolation by distance in Arabidopsis thaliana: biogeography and postglacial colonization of Europe

Arabidopsis thaliana provides a useful model system for functional, evolutionary and ecological studies in plant biology. We have analysed natural genetic variation in A. thaliana in order to infer its biogeographical and historical distribution across Eurasia. We analysed 79 amplified fragment length polymorphism (AFLP) markers in 142 accessions from the species' native range, and found highly significant genetic isolation by distance among A. thaliana accessions from Eurasia and southern Europe. These spatial patterns of genetic variation suggest that A. thaliana colonized central and northern Europe from Asia and from Mediterranean Pleistocene refugia, a trend which has been identified in other species. Statistically significant levels of multilocus linkage disequilibrium suggest intermediate levels of disequilibrium among subsets of loci, and analysis of genetic relationships among accessions reveal a star or bush-like dendrogram with low bootstrap support. Taken together, it appears that there has been sufficient historical recombination in the A. thaliana genome such that accessions do not conform to a tree-like, bifurcating pattern of evolution - there is no 'ecotype phylogeny.' Nonetheless, significant isolation by distance provides a framework upon which studies of natural variation in A. thaliana may be designed and interpreted.     

Patterns of divergence among conifer ESTs and polymorphism in Pinus sylvestris identify putative selective sweeps

Finding genes that are under positive selection is a difficult task, especially in non-model organisms. Here, we have analyzed expressed sequence tag (EST) data from 4 species (Pinus pinaster, Pinus taeda, Picea glauca, and Pseudotsuga menziesii) to investigate selection patterns during their evolution and to identify genes likely to be under positive selection. To confirm selection, population samples of these genes have been sequenced in Pinus sylvestris, a species that was not included in the EST data set. The estimates of branch-specific Ka/Ks (nonsynonymous/synonymous substitution rates) across all genes in the EST data set were similar or smaller than estimates from other higher plant species. There was no evidence for the traditional indication of positive selection, Ka/Ks above 1. However, several lines of evidence based on polymorphism patterns suggest that genes with high Ka/Ks (0.20-0.52) in the EST data set are in fact more affected by positive selection in P. sylvestris than genes with low Ka/Ks (0.01-0.04). The high Ka/Ks genes have a lower level of polymorphism and more negative Tajima's D than the low Ka/Ks genes. Further, in the high Ka/Ks group, the Hudson-Kreitman-Aguade test is significant. This suggests that the EST data set is a good starting point for finding genes under positive selection in conifers and that even moderate Ka/Ks values could be indicative of selection. A group of 5 genes with high Ka/Ks collectively show evidence for positive selection within P. sylvestris.     

Rates of molecular evolution in nuclear genes of east Mediterranean scorpions

Scorpions of the genus Mesobuthus represent a useful terrestrial model system for studying molecular evolution. They are distributed on several Aegean islands and the adjacent mainland, they are believed to have low rates of dispersal, and evolutionary divergence dates of taxa are available based on biogeographic events that separated islands from each other and the mainland. Here, we present data on polymorphism and synonymous (Ks) and non-synonymous (Ka) substitution rates for nine nuclear protein-coding genes of two east Mediterranean scorpion species, Mesobuthus gibbosus and M. cyprius (Buthidae). Levels of polymorphism tend to be lower in populations from islands (mean nucleotide diversity pi = 0.0071 +/- 0.0028) than in mainland populations (mean pi = 0.0201 +/- 0.0085). By using linear regression of genetic divergence versus isolation time, we estimate Ks to be 3.17 +/- 1.54 per (site x 10(9) years), and Ka to be 0.39 +/- 0.94 per (site x 10(9) years). These estimates for both Ks and Ka are considerably lower than for many other invertebrates, such as Drosophila, and may be attributed to scorpions' mammal-like generation times (approximately 2 years) and low metabolic rates. Phylogenetic analysis using maximum likelihood revealed a phylogeny that is congruent with that expected based on biogeographic events and in which divergences at synonymous sites are proportional to the dates that the taxa are believed to have split. Tests of equality of branch lengths for the Cyprus and Crete lineages revealed that Ks-estimates are about the same in both lineages, as expected from the biogeographic events that separated the islands, but Ka was increased in the Cyprus lineage compared to the Cretan lineage.     

Evidence for a large-scale population structure among accessions of Arabidopsis thaliana: possible causes and consequences for the distribution of linkage disequilibrium

The existence of a large-scale population structure was investigated in Arabidopsis thaliana by studying patterns of polymorphism in a set of 71 European accessions. We used sequence polymorphism surveyed in 10 fragments of approximately 600 nucleotides and a set of nine microsatellite markers. Population structure was investigated using a model-based inference framework. Among the accessions studied, the presence of four groups was inferred using genetic data, without using prior information on the geographical origin of the accessions. Significant genetic isolation by geographical distance was detected at the group level, together with a geographical gradient in allelic richness across groups. These results are discussed with respect to the previously proposed scenario of postglacial colonization of Europe from putative glacial refugia. Finally, the contribution of the inferred structure to linkage disequilibrium among 171 pairs of essentially unlinked markers was also investigated. Linkage disequilibrium analysis revealed that significant associations detected in the whole sample were mainly due to genetic differentiation among the inferred groups. We discuss the implication of this finding for future association studies in A. thaliana.     

The pattern of insertion/deletion polymorphism in Arabidopsis thaliana

Little is known about variation of nucleotide insertion/deletions (indels) within species. In Arabidopsis thaliana, we investigated indel polymorphism patterns between two genome sequences and among 96 accessions at 1215 loci. Our study identified patterns in the variation of indel density, size, GC content and distribution, and a correlation between indels and substitutions. We found that the GC content in indel sequences was lower than that in non-indel sequences and that indels typically occur in regions with lower GC content. Patterns of indel frequency distribution among populations were more consistent with neutral expectation than substitution patterns. We also found that the local level of substitutions is positively correlated with indel density and negatively correlated with their distance to the closed indel, suggesting that indels play an important role in nucleotide variation.     

A multilocus sequence survey in Arabidopsis thaliana reveals a genome-wide departure from a neutral model of DNA sequence polymorphism

The simultaneous analysis of multiple genomic loci is a powerful approach to studying the effects of population history and natural selection on patterns of genetic variation of a species. By surveying nucleotide sequence polymorphism at 334 randomly distributed genomic regions in 12 accessions of Arabidopsis thaliana, we examined whether a standard neutral model of nucleotide sequence polymorphism is consistent with observed data. The average nucleotide diversity was 0.0071 for total sites and 0.0083 for silent sites. Although levels of diversity are variable among loci, no correlation with local recombination rate was observed, but polymorphism levels were correlated for physically linked loci (<250 kb). We found that observed distributions of Tajima's D- and D/D(min)- and of Fu and Li's D-, D*- and F-, F*-statistics differed significantly from the expected distributions under a standard neutral model due to an excess of rare polymorphisms and high variances. Observed and expected distributions of Fay and Wu's H were not different, suggesting that demographic processes and not selection at multiple loci are responsible for the deviation from a neutral model. Maximum-likelihood comparisons of alternative demographic models like logistic population growth, glacial refugia, or past bottlenecks did not produce parameter estimates that were more consistent with observed patterns. However, exclusion of highly polymorphic "outlier loci" resulted in a fit to the logistic growth model. Various tests of neutrality revealed a set of candidate loci that may evolve under selection.     

Molecular population genetics of accessory gland protein genes and testis-expressed genes in Drosophila mojavensis and D. arizonae

Molecular population genetic investigation of Drosophila male reproductive genes has focused primarily on melanogaster subgroup accessory gland protein genes (Acp's). Consistent with observations from male reproductive genes of numerous taxa, Acp's evolve more rapidly than nonreproductive genes. However, within the Drosophila genus, large data sets from additional types of male reproductive genes and from different species groups are lacking. Here we report findings from a molecular population genetics analysis of male reproductive genes of the repleta group species, Drosophila arizonae and D. mojavensis. We find that Acp's have dramatically higher average pairwise Ka/Ks (0.93) than testis-enriched genes (0.19) and previously reported melanogaster subgroup Acp's (0.42). Overall, 10 of 19 Acp's have Ka/Ks > 1 either in nonpolarized analyses or in at least one lineage of polarized analyses. Of the nine Acp's for which outgroup data were available, average Ka/Ks was considerably higher in D. mojavensis (2.08) than in D. arizonae (0.87). Contrasts of polymorphism and divergence suggest that adaptive protein evolution at Acp's is more common in D. mojavensis than in D. arizonae.     

Natural selection of protein structural and functional properties: a single nucleotide polymorphism perspective

Background

The rates of molecular evolution for protein-coding genes depend on the stringency of functional or structural constraints. The Ka/Ks ratio has been commonly used as an indicator of selective constraints and is typically calculated from interspecies alignments. Recent accumulation of single nucleotide polymorphism (SNP) data has enabled the derivation of Ka/Ks ratios for polymorphism (SNP A/S ratios).

Results

Using data from the dbSNP database, we conducted the first large-scale survey of SNP A/S ratios for different structural and functional properties. We confirmed that the SNP A/S ratio is largely correlated with Ka/Ks for divergence. We observed stronger selective constraints for proteins that have high mRNA expression levels or broad expression patterns, have no paralogs, arose earlier in evolution, have natively disordered regions, are located in cytoplasm and nucleus, or are related to human diseases. On the residue level, we found higher degrees of variation for residues that are exposed to solvent, are in a loop conformation, natively disordered regions or low complexity regions, or are in the signal peptides of secreted proteins. Our analysis also revealed that histones and protein kinases are among the protein families that are under the strongest selective constraints, whereas olfactory and taste receptors are among the most variable groups.

Conclusion

Our study suggests that the SNP A/S ratio is a robust measure for selective constraints. The correlations between SNP A/S ratios and other variables provide valuable insights into the natural selection of various structural or functional properties, particularly for human-specific genes and constraints within the human lineage.     

Comparative analysis of NBS domain sequences of NBS-LRR disease resistance genes from sunflower, lettuce, and chicory

Plant resistance to many types of pathogens and pests can be achieved by the presence of disease resistance (R) genes. The nucleotide binding site-leucine rich repeat (NBS-LRR) class of R-genes is the most commonly isolated class of R-genes and makes up a super-family, which is often arranged in the genome as large multi-gene clusters. The NBS domain of these genes can be targeted by polymerase chain reaction (PCR) amplification using degenerate primers. Previous studies have used PCR derived NBS sequences to investigate both ancient R-gene evolution and recent evolution within specific plant families. However, comparative studies with the Asteraceae family have largely been ignored. In this study, we address recent evolution of NBS sequences within the Asteraceae and extend the comparison to the Arabidopsis thaliana genome. Using multiple sets of primers, NBS fragments were amplified from genomic DNA of three species from the family Asteraceae: Helianthus annuus (sunflower), Lactuca sativa (lettuce), and Cichorium intybus (chicory). Analysis suggests that Asteraceae species share distinct families of R-genes, composed of genes related to both coiled-coil (CC) and toll-interleukin-receptor homology (TIR) domain containing NBS-LRR R-genes. Between the most closely related species, (lettuce and chicory) a striking similarity of CC subfamily composition was identified, while sunflower showed less similarity in structure. These sequences were also compared to the A. thaliana genome. Asteraceae NBS gene subfamilies appear to be distinct from Arabidopsis gene clades. These data suggest that NBS families in the Asteraceae family are ancient, but also that gene duplication and gene loss events occur and change the composition of these gene subfamilies over time.     

Multilocus analysis of variation using a large empirical data set: phenylpropanoid pathway genes in Arabidopsis thaliana

Detecting the signature of adaptation on nucleotide variation is often difficult in species that like Arabidopsis thaliana might have a complex demographic history. Recent re-sequencing surveys in this species provided genome-wide information that would mainly reflect its demographic history. We have used a large empirical data set (LED) as well as multilocus coalescent simulations to analyse sequence variation at loci involved in the phenylpropanoid pathway of this species. We surveyed and examined DNA sequence variation at nine of these loci (about 19.7 kb) in 23 accessions of A. thaliana and one accession of its closely related species Arabidopsis lyrata . Nucleotide variation was lower at nonsynonymous sites than at silent sites in all loci, indicating generalized functional constraint at the protein level. No association between variation and position in the metabolic pathway was detected. When the data were contrasted against the standard neutral model, significant deviations for silent variation were detected with Tajima's D , Fu's F_S and Fay and Wu's H multilocus test statistics. These deviations were in the same direction than in previous large-scale multilocus analyses, suggesting a genome-wide effect. When the nine-locus data set was contrasted against the large empirical data set, the level (Watterson's θ) and pattern of variation (Tajima's D ) detected in these loci did not deviate either at the single-locus or multilocus level from the corresponding empirical distributions. These results would support an important role of the demographic history of A. thaliana in shaping nucleotide variation at the nine studied phenylpropanoid loci. The potential and limitations of the empirical distribution approach are discussed.     

Receptor-like genes in the major resistance locus of lettuce are subject to divergent selection.

Disease resistance genes in plants are often found in complex multigene families. The largest known cluster of disease resistance specificities in lettuce contains the RGC2 family of genes. We compared the sequences of nine full-length genomic copies of RGC2 representing the diversity in the cluster to determine the structure of genes within this family and to examine the evolution of its members. The transcribed regions range from at least 7.0 to 13.1 kb, and the cDNAs contain deduced open reading frames of approximately 5. 5 kb. The predicted RGC2 proteins contain a nucleotide binding site and irregular leucine-rich repeats (LRRs) that are characteristic of resistance genes cloned from other species. Unique features of the RGC2 gene products include a bipartite LRR region with >40 repeats. At least eight members of this family are transcribed. The level of sequence diversity between family members varied in different regions of the gene. The ratio of nonsynonymous (Ka) to synonymous (Ks) nucleotide substitutions was lowest in the region encoding the nucleotide binding site, which is the presumed effector domain of the protein. The LRR-encoding region showed an alternating pattern of conservation and hypervariability. This alternating pattern of variation was also found in all comparisons within families of resistance genes cloned from other species. The Ka /Ks ratios indicate that diversifying selection has resulted in increased variation at these codons. The patterns of variation support the predicted structure of LRR regions with solvent-exposed hypervariable residues that are potentially involved in binding pathogen-derived ligands.     

An analysis of microsatellite loci in Arabidopsis thaliana: mutational dynamics and application

Microsatellite loci are among the most commonly used molecular markers. These loci typically exhibit variation for allele frequency distribution within a species. However, the factors contributing to this variation are not well understood. To expand on the current knowledge of microsatellite evolution, 20 microsatellite loci were examined for 126 accessions of the flowering plant, Arabidopsis thaliana. Substantial variability in mutation pattern among loci was found, most of which cannot be explained by the assumptions of the traditional stepwise mutation model or infinite alleles model. Here it is shown that the degree of locus diversity is strongly correlated with the number of contiguous repeats, more so than with the total number of repeats. These findings support a strong role for repeat disruptions in stabilizing microsatellite loci by reducing the substrate for polymerase slippage and recombination. Results of cluster analyses are also presented, demonstrating the potential of microsatellite loci for resolving relationships among accessions of A. thaliana.     

Length polymorphism and allele structure of trinucleotide microsatellites in natural accessions of Arabidopsis thaliana

 The objective of this work was to assess the degree of trinucleotide microsatellite length polymorphism in the selfing species Arabidopsis thaliana. PCR amplifications of 12 microsatellite loci among 49 natural populations revealed between one to eight length variants (alleles) for each locus. The average number of alleles per locus was four and the average genetic diversity index was 0.43. Divergence between length variants was investigated at the nucleotide level. Several observations emerge from the sequence data: (1) for most loci, length polymorphism results only from variations in the number of trinucleotide repeats; (2) for a few others, some variability was noted in the flanking sequences; (3) for compound and interrupted loci containing two arrays of trinucleotide repeats, length variations preferentially affect the longest one. Five of the Arabidopsis thaliana accessions were clearly composed of two sublines. In 2 other accessions, some heterozygous individual plants, probably resulting from recent outcrosses, were found. A phylogenetic tree constructed on the basis of trinucleotide microsatellite allelic diversity shows that genetic relationships among the accessions are not correlated with their geographic origin. Received: 4 November 1997 / Accepted: 3 March 1998     

Epistatic interactions among herbicide resistances in Arabidopsis thaliana: the fitness cost of multiresistance

The type of interactions among deleterious mutations is considered to be crucial in numerous areas of evolutionary biology, including the evolution of sex and recombination, the evolution of ploidy, the evolution of selfing, and the conservation of small populations. Because the herbicide resistance genes could be viewed as slightly deleterious mutations in the absence of the pesticide selection pressure, the epistatic interactions among three herbicide resistance genes (acetolactate synthase CSR, cellulose synthase IXR1, and auxin-induced AXR1 target genes) were estimated in both the homozygous and the heterozygous states, giving 27 genotype combinations in the model plant Arabidopsis thaliana. By analyzing eight quantitative traits in a segregating population for the three herbicide resistances in the absence of herbicide, we found that most interactions in both the homozygous and the heterozygous states were best explained by multiplicative effects (each additional resistance gene causes a comparable reduction in fitness) rather than by synergistic effects (each additional resistance gene causes a disproportionate fitness reduction). Dominance coefficients of the herbicide resistance cost ranged from partial dominance to underdominance, with a mean dominance coefficient of 0.07. It was suggested that the csr1-1, ixr1-2, and axr1-3 resistance alleles are nearly fully recessive for the fitness cost. More interestingly, the dominance of a specific resistance gene in the absence of herbicide varied according to, first, the presence of the other resistance genes and, second, the quantitative trait analyzed. These results and their implications for multiresistance evolution are discussed in relation to the maintenance of polymorphism at resistance loci in a heterogeneous environment.     

DNA variation in the acidic chitinase locus (ChiA) region in Arabis gemmifera and its related species

We analyzed DNA variation at the acidic chitinase (ChiA) locus of Arabis gemmifera and among its eight related species. Nucleotide diversity (pi) of the entire locus in A. gemmifera was 0.0032, which was one third that of A. thaliana. In A. gemmifera, an excess of unique polymorphisms yielded significantly negative results with the tests of Tajima and Fu and Li. The McDonald and Kreitman test revealed that the ratio of nucleotide replacement to synonymous changes in A. gemmifera was significantly greater than those between A. gemmifera and A. glabra, A. gemmifera and A. griffithiana, A. gemmifera and A. korshinskyi, A. gemmifera and A. wallichii, and A. gemmifera and A. himalaica. These results indicated that the neutrality assumption, the equilibrium population assumption, or both, could not be applied to the ChiA locus of A. gemmifera. The small size and relative isolation of local subpopulations of A. gemmifera could explain the excess of unique polymorphisms and the high proportion of replacement changes. The specific sampling scheme of this study, where one strain each was sampled from each local subpopulation might also have led to an excess of singletons. Interspecific comparison among Arabidopsis, Arabis and Cardaminopsis species showed that Ka was always lower than Ks, providing evidence against the adaptive evolution of ChiA. However, Ka/Ks was greater between closely related species than between more distant related species. ChiA had a higher level of replacement divergence and a lower level of synonymous divergence compared than did Adh. We suggest that both the mutation rate at the nucleotide level and the selective constraints at the protein level are lower in ChiA than in Adh.     

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.