A Nonparametric Test Reveals Selection for Rapid Flowering in the Arabidopsis Genome

The detection of footprints of natural selection in genetic polymorphism data is fundamental to understanding the genetic basis of adaptation, and has important implications for human health. The standard approach has been to reject neutrality in favor of selection if the pattern of variation at a candidate locus was significantly different from the predictions of the standard neutral model. The problem is that the standard neutral model assumes more than just neutrality, and it is almost always possible to explain the data using an alternative neutral model with more complex demography. Today's wealth of genomic polymorphism data, however, makes it possible to dispense with models altogether by simply comparing the pattern observed at a candidate locus to the genomic pattern, and rejecting neutrality if the pattern is extreme. Here, we utilize this approach on a truly genomic scale, comparing a candidate locus to thousands of alleles throughout the Arabidopsis thaliana genome. We demonstrate that selection has acted to increase the frequency of early-flowering alleles at the vernalization requirement locus FRIGIDA. Selection seems to have occurred during the last several thousand years, possibly in response to the spread of agriculture. We introduce a novel test statistic based on haplotype sharing that embraces the problem of population structure, and so should be widely applicable.     

The Pattern of Polymorphism in Arabidopsis thaliana

We resequenced 876 short fragments in a sample of 96 individuals of Arabidopsis thaliana that included stock center accessions as well as a hierarchical sample from natural populations. Although A. thaliana is a selfing weed, the pattern of polymorphism in general agrees with what is expected for a widely distributed, sexually reproducing species. Linkage disequilibrium decays rapidly, within 50 kb. Variation is shared worldwide, although population structure and isolation by distance are evident. The data fail to fit standard neutral models in several ways. There is a genome-wide excess of rare alleles, at least partially due to selection. There is too much variation between genomic regions in the level of polymorphism. The local level of polymorphism is negatively correlated with gene density and positively correlated with segmental duplications. Because the data do not fit theoretical null distributions, attempts to infer natural selection from polymorphism data will require genome-wide surveys of polymorphism in order to identify anomalous regions. Despite this, our data support the utility of A. thaliana as a model for evolutionary functional genomics.     

Hitchhiking mapping reveals a candidate genomic region for natural selection in three-spined stickleback chromosome VIII

Identification of genes and genomic regions under directional natural selection has become one of the major goals in evolutionary genetics, but relatively little work to this end has been done by applying hitchhiking mapping to wild populations. Hitchhiking mapping starts from a genome scan using a randomly spaced set of molecular markers followed by a fine-scale analysis in the flanking regions of the candidate regions under selection. We used the hitchhiking mapping approach to narrow down a selective sweep in the genomic region flanking a candidate locus (Stn90) in chromosome VIII in the three-spined stickleback (Gasterosteus aculeatus). Twenty-four microsatellite markers were screened in an approximately 800-kb region around the candidate locus in three marine and four freshwater populations. The patterns of genetic diversity and differentiation in the candidate region were compared to those of a putatively neutral set of markers. The Bayesian FST-test indicated an elevated genetic differentiation, deviating significantly from neutral expectations, at a continuous region of approximately 20 kb upstream from the candidate locus. Furthermore, a method developed for an array of microsatellite markers rejected neutrality in a region of approximately 90 kb flanking the candidate locus supporting the selective sweep hypothesis. Likewise, the genomewide pattern of genetic diversity differed from the candidate region in a bottleneck analysis suggesting that selection, rather than demography, explains the reduced genetic diversity at the candidate interval. The neutrality tests suggest that the selective sweep had occurred mainly in the Lake Pulmanki population, but the results from bottleneck analyses indicate that selection might have operated in other populations as well. These results suggest that the narrow interval around locus Stn90 has likely been under directional selection, but the region contains several predicted genes, each of which can be the actual targets of selection. Understanding of the functional significance of this genomic region in an ecological context will require a more detailed sequence analysis.     

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.     

Genetic analysis of differentiation among breeding ponds reveals a candidate gene for local adaptation in Rana arvalis

One of the main questions in evolutionary and conservation biology is how geographical and environmental features of the landscape shape neutral and adaptive genetic variation in natural populations. The identification of genomic polymorphisms that account for adaptive variation can aid in finding candidate loci for local adaptation. Consequently, a comparison of spatial patterns in neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection at the landscape scale. Many amphibians breed in wetlands, which differ in environmental conditions and in the degree of isolation, enhancing the potential for local adaptation. We used microsatellite markers to measure genetic differentiation among 17 local populations of Rana arvalis breeding in a network of wetlands. We found that locus RC08604 deviated from neutral expectations, suggesting that it is a good candidate for directional selection. We used a genetic network analysis to show that the allele distribution in this locus is correlated with habitat characteristics, whereas this was not the case at neutral markers that displayed a different allele distribution and population network in the study area. The graph approach illustrated the genomic heterogeneity (neutral loci vs. the candidate locus for directional selection) of gene exchange and genetic divergence among populations under directional selection. Limited gene flow between wetlands was only observed at the candidate genomic region under directional selection. RC08604 is partially located inside an up‐regulated thyroid‐hormone receptor (TRβ) gene coordinating the expression of other genes during metamorphosis and appears to be linked with variation in larval life‐history traits found among R. arvalis populations. We suggest that directional selection on genes coding larval life‐history traits is strong enough to maintain the divergence in these genomic regions, reducing the effective recombination of locally adapted alleles but not in other regions of the genome. Integrating this knowledge into conservation plans at the landscape scale will improve the design of management strategies to preserve adaptive genetic diversity in wetland networks.     

The Case for Selection at CCR5-Δ32

The C-C chemokine receptor 5, 32 base-pair deletion (CCR5-Δ32) allele confers strong resistance to infection by the AIDS virus HIV. Previous studies have suggested that CCR5-Δ32 arose within the past 1,000 y and rose to its present high frequency (5%–14%) in Europe as a result of strong positive selection, perhaps by such selective agents as the bubonic plague or smallpox during the Middle Ages. This hypothesis was based on several lines of evidence, including the absence of the allele outside of Europe and long-range linkage disequilibrium at the locus. We reevaluated this evidence with the benefit of much denser genetic maps and extensive control data. We find that the pattern of genetic variation at CCR5-Δ32 does not stand out as exceptional relative to other loci across the genome. Moreover using newer genetic maps, we estimated that the CCR5-Δ32 allele is likely to have arisen more than 5,000 y ago. While such results can not rule out the possibility that some selection may have occurred at C-C chemokine receptor 5 (CCR5), they imply that the pattern of genetic variation seen atCCR5-Δ32 is consistent with neutral evolution. More broadly, the results have general implications for the design of future studies to detect the signs of positive selection in the human genome.     

Testing the generalized neutrality hypothesis

Various tests of the hypothesis of selective neutrality based on gene frequency are now available. These tests take as null hypothesis the concept of “strict neutrality”: all new mutants are required to be selectively identical to each other. For evolutionary questions, however, (as opposed to those of genetic polymorphism), a wider null hypothesis might be of interest. Since deleterious alleles have essentially no evolutionary importance, one might wish to test the null hypothesis that only neutral or deleterious mutations occur. The principal alternative to this hypothesis is that there exists heterotic selection of some form for some alleles tending to maintain a level of genetic polymorphism higher than that under neutrality. In this paper an assessment is made of the usefulness of a test of strict neutrality first proposed by this author (Ewens, 1972) as a test of null hypothesis of “generalized neutrality,” i.e. that only neutral or deleterious alleles occur. At the same time some remarks will be made about estimation of the fundamental parameter θ defining these processes.     

Functional Divergence Caused by Ancient Positive Selection of a Drosophila Hybrid Incompatibility Locus

Interspecific hybrid lethality and sterility are a consequence of divergent evolution between species and serve to maintain the discrete identities of species. The evolution of hybrid incompatibilities has been described in widely accepted models by Dobzhansky and Muller where lineage-specific functional divergence is the essential characteristic of hybrid incompatibility genes. Experimentally tractable models are required to identify and test candidate hybrid incompatibility genes. Several Drosophila melanogaster genes involved in hybrid incompatibility have been identified but none has yet been shown to have functionally diverged in accordance with the Dobzhansky-Muller model. By introducing transgenic copies of the X-linked Hybrid male rescue (Hmr) gene into D. melanogaster from its sibling species D. simulans and D. mauritiana, we demonstrate that Hmr has functionally diverged to cause F1 hybrid incompatibility between these species. Consistent with the Dobzhansky-Muller model, we find that Hmr has diverged extensively in the D. melanogaster lineage, but we also find extensive divergence in the sibling-species lineage. Together, these findings implicate over 13% of the amino acids encoded by Hmr as candidates for causing hybrid incompatibility. The exceptional level of divergence at Hmr cannot be explained by neutral processes because we use phylogenetic methods and population genetic analyses to show that the elevated amino-acid divergence in both lineages is due to positive selection in the distant past—at least one million generations ago. Our findings suggest that multiple substitutions driven by natural selection may be a general phenomenon required to generate hybrid incompatibility alleles.     

Allelic diversity at the Mhc-DQA locus of woodmouse populations (Apodemus sylvaticus) present in the islands and mainland of the northern Mediterranean

Aim Polymorphism at neutral markers and at MHC loci in rodent populations living on islands is generally low. The main genetic factors that may contribute to a reduced level of genetic variability are genetic drift, reduced gene flow and founder events. We investigated the pattern of polymorphism at the second exon of the Mhc‐DQA gene in island populations of Apodemus sylvaticus and in their mainland counterparts to investigate the pattern of MHC polymorphism in a phylogeographical context and to assess the impact of insularity on diversity at this locus. Location Eight north Mediterranean populations of Apodemus sylvaticus were studied, including five island populations (Majorca, Minorca, Porquerolles, Port‐Cros and Sicily) and three mainland populations. Methods cDNA sequencing and nucleotide sequences analyses. Synonymous and non‐synonymous substitutions were examined at the PBR and non‐PBR sites. The DQA allelic distribution in populations was compared with the woodmouse phylogeography. Results This study presents novel DQA alleles. High polymorphism of the DQA locus is recorded in natural populations of A. sylvaticus with 13 alleles being widely distributed irrespective of the geographical origin and palaeoclimatic history of populations. The DQA locus does not show the expected pattern for non‐synonymous substitutions at the PBR sites. However, island populations show a weak loss of polymorphism in comparison with their mainland counterparts. Main conclusions The DQA locus in the woodmouse seems to be subject to weak selection and does not allow resolution of phylogeographical relationships among European woodmouse populations. The presence of at least three alleles in island populations and the maintenance of five alleles between the two European lineages over 1.5 Myr suggest that balancing selection may act within populations, and more precisely within island populations, to maintain genetic variability. This study shows that phylogeographical studies are a prerequisite for any genetic investigation of selected genes in natural populations.     

The pattern of polymorphism in Arabidopsis thaliana

We resequenced 876 short fragments in a sample of 96 individuals of Arabidopsis thaliana that included stock center accessions as well as a hierarchical sample from natural populations. Although A. thaliana is a selfing weed, the pattern of polymorphism in general agrees with what is expected for a widely distributed, sexually reproducing species. Linkage disequilibrium decays rapidly, within 50 kb. Variation is shared worldwide, although population structure and isolation by distance are evident. The data fail to fit standard neutral models in several ways. There is a genome-wide excess of rare alleles, at least partially due to selection. There is too much variation between genomic regions in the level of polymorphism. The local level of polymorphism is negatively correlated with gene density and positively correlated with segmental duplications. Because the data do not fit theoretical null distributions, attempts to infer natural selection from polymorphism data will require genome-wide surveys of polymorphism in order to identify anomalous regions. Despite this, our data support the utility of A. thaliana as a model for evolutionary functional genomics.     

Bayesian analysis suggests that most amino acid replacements in Drosophila are driven by positive selection

One of the principal goals of population genetics is to understand the processes by which genetic variation within species (polymorphism) becomes converted into genetic differences between species (divergence). In this transformation, selective neutrality, near neutrality, and positive selection may each play a role, differing from one gene to the next. Synonymous nucleotide sites are often used as a uniform standard of comparison across genes on the grounds that synonymous sites are subject to relatively weak selective constraints and so may, to a first approximation, be regarded as neutral. Synonymous sites are also interdigitated with nonsynonymous sites and so are affected equally by genomic context and demographic factors. Hence a comparison of levels of polymorphism and divergence between synonymous sites and amino acid replacement sites in a gene is potentially informative about the magnitude of selective forces associated with amino acid replacements. We have analyzed 56 genes in which polymorphism data from D. simulans are compared with divergence from a reference strain of D. melanogaster. The framework of the analysis is Bayesian and assumes that the distribution of selective effects (Malthusian fitnesses) is Gaussian with a mean that differs for each gene. In such a model, the average scaled selection intensity (gamma = N(e)s) of amino acid replacements eligible to become polymorphic or fixed is -7.31, and the standard deviation of selective effects within each locus is 6.79 (assuming homoscedasticity across loci). For newly arising mutations of this type that occur in autosomal or X-linked genes, the average proportion of beneficial mutations is 19.7%. Among the amino acid polymorphisms in the sample, the expected average proportion of beneficial mutations is 47.7%, and among amino acid replacements that become fixed the average proportion of beneficial mutations is 94.3%. The average scaled selection intensity of fixed mutations is +5.1. The presence of positive selection is pervasive with the single exception of kl-5, a Y-linked fertility gene. We find no evidence that a significant fraction of fixed amino acid replacements is neutral or nearly neutral or that positive selection drives amino acid replacements at only a subset of the loci. These results are model dependent and we discuss possible modifications of the model that might allow more neutral and nearly neutral amino acid replacements to be fixed.     

Contrasted patterns of selection on MHC-linked microsatellites in natural populations of the Malagasy plague reservoir

Plague (Yersinia pestis infection) is a highly virulent rodent disease that persists in many natural ecosystems. The black rat (Rattus rattus) is the main host involved in the plague focus of the central highlands of Madagascar. Black rat populations from this area are highly resistant to plague, whereas those from areas in which the disease is absent (low altitude zones of Madagascar) are susceptible. Various lines of evidence suggest a role for the Major Histocompatibility Complex (MHC) in plague resistance. We therefore used the MHC region as a candidate for detecting signatures of plague-mediated selection in Malagasy black rats, by comparing population genetic structures for five MHC-linked microsatellites and neutral markers in two sampling designs. We first compared four pairs of populations, each pair including one population from the plague focus and one from the disease-free zone. Plague-mediated selection was expected to result in greater genetic differentiation between the two zones than expected under neutrality and this was observed for one MHC-class I-linked locus (D20Img2). For this marker as well as for four other MHC-linked loci, a geographic pattern of genetic structure was found at local scale within the plague focus. This pattern would be expected if plague selection pressures were spatially variable. Finally, another MHC-class I-linked locus (D20Rat21) showed evidences of balancing selection, but it seems more likely that this selection would be related to unknown pathogens more widely distributed in Madagascar than plague.     

Ascertainment Bias and the Pattern of Nucleotide Diversity at the Human ALDH2 Locus in a Japanese Population

Many East Asian human populations harbor a high-frequency deficiency allele for the aldehyde dehydrogenase 2 (ALDH2) enzyme, a critical protein involved in the metabolism of ethanol. Here we use resequencing and long-range SNP haplotype data from a Japanese sample to test whether patterns of nucleotide diversity and linkage disequilibrium at this locus are compatible with a standard neutral model of evolution. Examination of the pattern of polymorphism at a locus such as this, where the frequency of a common allele is known a priori, introduces an ascertainment bias that must be corrected for in analyses of the frequency spectrum of polymorphisms. We apply a flexible and generally applicable simulation approach to correct for this bias in our ALDH2 data and, also, to explore the effect of bias on the commonly used summary statistics Tajima’s D, Fu and Li’s D, and Fay and Wu’s H. Our study finds no evidence that the pattern of genetic variation at ALDH2 differs from that expected under a standard neutral model. However, our general examination of ascertainment bias indicates that a priori knowledge of segregating alleles greatly affects the expected distributions of summary statistics. Under many parameter combinations we find that ascertainment bias introduces an elevated rate of false positives when summary statistics are used to test for deviations from a standard neutral model. However, we also show that over a wide range of conditions the power of all summary statistics can be greatly increased by incorporating prior knowledge of segregating alleles. [Reviewing Editor: Dr. Martin Kreitman]     

Effect of natural selection on the duplicated lysyl oxidase gene in Atlantic salmon

We examined the polymorphism of the lysyl oxidase (LOX) locus, involved in the initiation of muscle collagen cross-linking, in three populations of Atlantic salmon with different life histories and growth rates and compared it with a closely related species (rainbow trout). Up to four alleles were observed per individual, probably as a consequence of the tetraploid origin of the salmonid genome. We found high polymorphism in the LOX locus (16 alleles expressed in total and several low frequency private alleles) in two natural Atlantic salmon populations and extremely reduced diversity in a farmed population (3 alleles) with low density of collagen crosslinks. We also assessed the relative role of selection in maintaining LOX genetic variability in Atlantic salmon. Results from several neutrality tests suggest that selection is playing a role in shaping diversity at the LOX locus. Positive selection was inferred by three different likelihood phylogeny-based methods and one selected site, identified by all three different methods (PAML, FEL and REL) was located within the “copper-talon” characteristic of LOX proteins. We suggest that the retention of four alleles in the salmon LOX locus could be related to its multiple functions.     

Effect of balancing selection on spatial genetic structure within populations: theoretical investigations on the self-incompatibility locus and empirical studies in Arabidopsis halleri

The effect of selection on patterns of genetic structure within and between populations may be studied by contrasting observed patterns at the genes targeted by selection with those of unlinked neutral marker loci. Local directional selection on target genes will produce stronger population genetic structure than at neutral loci, whereas the reverse is expected for balancing selection. However, theoretical predictions on the intensity of this signal under precise models of balancing selection are still lacking. Using negative frequency-dependent selection acting on self-incompatibility systems in plants as a model of balancing selection, we investigated the effect of such selection on patterns of spatial genetic structure within a continuous population. Using numerical simulations, we tested the effect of the type of self-incompatibility system, the number of alleles at the self-incompatibility locus and the dominance interactions among them, the extent of gene dispersal, and the immigration rate on spatial genetic structure at the selected locus and at unlinked neutral loci. We confirm that frequency-dependent selection is expected to reduce the extent of spatial genetic structure as compared to neutral loci, particularly in situations with low number of alleles at the self-incompatibility locus, high frequency of codominant interactions among alleles, restricted gene dispersal and restricted immigration from outside populations. Hence the signature of selection on spatial genetic structure is expected to vary across species and populations, and we show that empirical data from the literature as well as data reported here on three natural populations of the herb Arabidopsis halleri confirm these theoretical results.     

AFLP genome scans suggest divergent selection on colour patterning in allopatric colour morphs of a cichlid fish

Genome scan-based tests for selection are directly applicable to natural populations to study the genetic and evolutionary mechanisms behind phenotypic differentiation. We conducted AFLP genome scans in three distinct geographic colour morphs of the cichlid fish Tropheus moorii to assess whether the extant, allopatric colour pattern differentiation can be explained by drift and to identify markers mapping to genomic regions possibly involved in colour patterning. The tested morphs occupy adjacent shore sections in southern Lake Tanganyika and are separated from each other by major habitat barriers. The genome scans revealed significant genetic structure between morphs, but a very low proportion of loci fixed for alternative AFLP alleles in different morphs. This high level of polymorphism within morphs suggested that colour pattern differentiation did not result exclusively from neutral processes. Outlier detection methods identified six loci with excess differentiation in the comparison between a bluish and a yellow-blotch morph and five different outlier loci in comparisons of each of these morphs with a red morph. As population expansions and the genetic structure of Tropheus make the outlier approach prone to false-positive signals of selection, we examined the correlation between outlier locus alleles and colour phenotypes in a genetic and phenotypic cline between two morphs. Distributions of allele frequencies at one outlier locus were indeed consistent with linkage to a colour locus. Despite the challenges posed by population structure and demography, our results encourage the cautious application of genome scans to studies of divergent selection in subdivided and recently expanded populations.     

The Pattern of Neutral Molecular Variation under the Background Selection Model

Stochastic simulations of the infinite sites model were used to study the behavior of genetic diversity at a neutral locus in a genomic region without recombination, but subject to selection against deleterious alleles maintained by recurrent mutation (background selection). In large populations, the effect of background selection on the number of segregating sites approaches the effct on nucleotide site diversity, i.e., the reduction in genetic variability caused by background selection resembles that caused by a simple reduction in effective population size. We examined, by coalescence-based methods, the power of several tests for the departure from neutral expectation of the frequency spectra of alleles in samples from randomly mating populations (TAJIMA's, FU and LI's, and WATTERSON's tests). All of the tests have low power unless the selection against mutant alleles is extremely weak. In Drosophila, significant TAJIMA's tests are usually not obtained with empirical data sets from loci in genomic regions with restricted recombination frequencies and that exhibit low genetic diversity. This is consistent with the operation of background selection as opposed to selective sweeps. It remains to be decided whether background selection is sufficient to explain the observed extent of reduction in diversity in regions of restricted recombination.     

Unusual pattern of nucleotide sequence variation at the OS-E and OS-F genomic regions of Drosophila simulans

Nucleotide variation at the genomic region encompassing the odorant-binding protein genes OS-E and OS-F (OS region) was surveyed in two populations of Drosophila simulans, one from Europe and the other from Africa. We found that the European population shows an atypical and large haplotype structure, which extends throughout the approximately 5-kb surveyed genomic region. This structure is depicted by two major haplotype groups segregating at intermediate frequency in the sample, one haplogroup with nearly no variation, and the other at levels more typical for this species. This pattern of variation was incompatible with neutral predictions for a population at a stationary equilibrium. Nevertheless, neutrality tests contrasting polymorphism and divergence data fail to detect any departure from the standard neutral model in this species, whereas they confirm the non-neutral behavior previously observed at the OS-E gene in D. melanogaster. Although positive Darwinian selection may have been responsible for the observed unusual nucleotide variation structure, coalescent simulation results do not allow rejecting the hypothesis that the pattern was generated by a recent bottleneck in the history of European populations of D. simulans.     

Population Genomics: Whole-Genome Analysis of Polymorphism and Divergence in Drosophila simulans

The population genetic perspective is that the processes shaping genomic variation can be revealed only through simultaneous investigation of sequence polymorphism and divergence within and between closely related species. Here we present a population genetic analysis of Drosophila simulans based on whole-genome shotgun sequencing of multiple inbred lines and comparison of the resulting data to genome assemblies of the closely related species, D. melanogaster and D. yakuba. We discovered previously unknown, large-scale fluctuations of polymorphism and divergence along chromosome arms, and significantly less polymorphism and faster divergence on the X chromosome. We generated a comprehensive list of functional elements in the D. simulans genome influenced by adaptive evolution. Finally, we characterized genomic patterns of base composition for coding and noncoding sequence. These results suggest several new hypotheses regarding the genetic and biological mechanisms controlling polymorphism and divergence across the Drosophila genome, and provide a rich resource for the investigation of adaptive evolution and functional variation in D. simulans.