The recent demographic and adaptive history of Drosophila melanogaster

Population genetic analyses of the past two decades confirmed an earlier hypothesis by L Tsacas and D Lachaise that the cosmopolitan species Drosophila melanogaster has an Afrotropical origin, and that it colonized the rest of the world only very recently. Maximum likelihood analyses based on multilocus data suggest that the putative ancestral African population expanded its size about 60,000 years ago (ya). These demographic changes were accompanied by the fixation of numerous beneficial mutations, as revealed by signatures of positive directional selection in the genome (selective sweeps). The estimated rate of adaptive substitution on the X chromosome is in the order of 10(-11) per nucleotide site per generation. Comparable (but not significantly higher) substitution rates are found in derived populations that colonized new habitats outside Africa, such as in a European population that branched off from the African lineage about 16,000 ya.     

Clonal Interference in the Evolution of Influenza

The seasonal influenza A virus undergoes rapid evolution to escape human immune response. Adaptive changes occur primarily in antigenic epitopes, the antibody-binding domains of the viral hemagglutinin. This process involves recurrent selective sweeps, in which clusters of simultaneous nucleotide fixations in the hemagglutinin coding sequence are observed about every 4 years. Here, we show that influenza A (H3N2) evolves by strong clonal interference. This mode of evolution is a red queen race between viral strains with different beneficial mutations. Clonal interference explains and quantifies the observed sweep pattern: we find an average of at least one strongly beneficial amino acid substitution per year, and a given selective sweep has three to four driving mutations on average. The inference of selection and clonal interference is based on frequency time series of single-nucleotide polymorphisms, which are obtained from a sample of influenza genome sequences over 39 years. Our results imply that mode and speed of influenza evolution are governed not only by positive selection within, but also by background selection outside antigenic epitopes: immune adaptation and conservation of other viral functions interfere with each other. Hence, adapting viral proteins are predicted to be particularly brittle. We conclude that a quantitative understanding of influenza’s evolutionary and epidemiological dynamics must be based on all genomic domains and functions coupled by clonal interference.     

Evidence for a selective sweep in the wapl region of Drosophila melanogaster

A scan of the X chromosome of a European Drosophila melanogaster population revealed evidence for the recent action of positive directional selection at individual loci. In this study we analyze one such region that showed no polymorphism in the genome scan (located in cytological division 2C10-2E1). We detect a 60.5-kb stretch of DNA encompassing the genes ph-d, ph-p, CG3835, bcn92, Pgd, wapl, and Cyp4d1, which almost completely lacks variation in the European sample. Loci flanking this region show a skewed frequency spectrum at segregating sites, strong haplotype structure, and high levels of linkage disequilibrium. Neutrality tests reveal that these data are unlikely under both the neutral equilibrium model and the simple bottleneck scenarios. In contrast, newly developed maximum-likelihood ratio tests suggest that strong selection has acted recently on the region under investigation, causing a selective sweep. Evidence that this sweep may have originated in an ancestral population in Africa is presented.     

Nonadaptive explanations for signatures of partial selective sweeps in Drosophila

A beneficial mutation that has nearly but not yet fixed in a population produces a characteristic haplotype configuration, called a partial selective sweep. Whether nonadaptive processes might generate similar haplotype configurations has not been extensively explored. Here, we consider 5 population genetic data sets taken from regions flanking high-frequency transposable elements in North American strains of Drosophila melanogaster, each of which appears to be consistent with the expectations of a partial selective sweep. We use coalescent simulations to explore whether incorporation of the species' demographic history, purifying selection against the element, or suppression of recombination caused by the element could generate putatively adaptive haplotype configurations. Whereas most of the data sets would be rejected as nonneutral under the standard neutral null model, only the data set for which there is strong external evidence in support of an adaptive transposition appears to be nonneutral under the more complex null model and in particular when demography is taken into account. High-frequency, derived mutations from a recently bottlenecked population, such as we study here, are of great interest to evolutionary genetics in the context of scans for adaptive events; we discuss the broader implications of our findings in this context.     

The effect of deleterious alleles on adaptation in asexual populations

We calculate the fixation probability of a beneficial allele that arises as the result of a unique mutation in an asexual population that is subject to recurrent deleterious mutation at rate U. Our analysis is an extension of previous works, which make a biologically restrictive assumption that selection against deleterious alleles is stronger than that on the beneficial allele of interest. We show that when selection against deleterious alleles is weak, beneficial alleles that confer a selective advantage that is small relative to U have greatly reduced probabilities of fixation. We discuss the consequences of this effect for the distribution of effects of alleles fixed during adaptation. We show that a selective sweep will increase the fixation probabilities of other beneficial mutations arising during some short interval afterward. We use the calculated fixation probabilities to estimate the expected rate of fitness improvement in an asexual population when beneficial alleles arise continually at some low rate proportional to U. We estimate the rate of mutation that is optimal in the sense that it maximizes this rate of fitness improvement. Again, this analysis relaxes the assumption made previously that selection against deleterious alleles is stronger than on beneficial alleles.     

Evidence of gene conversion associated with a selective sweep in Drosophila melanogaster

Since Drosophila melanogaster colonized Europe from tropical Africa 10 to 15 thousand years ago, it is expected that adaptation has played a major role in this species in recent times. A previously conducted multilocus scan of noncoding DNA sequences on the X chromosome in an ancestral and a derived population of D. melanogaster revealed that some loci have been affected by directional selection in the European population. We investigated if the pattern of DNA sequence polymorphism in a region surrounding one of these loci can be explained by a hitchhiking event. We found strong evidence that the studied region around the gene unc-119 was shaped by a recent selective sweep, including a valley of reduced heterozygosity of 83.4 kb, a skew in the frequency spectrum, and significant linkage disequilibrium on one side of the valley. This region, however, was interrupted by gene conversion events leading to a strong haplotype structure in the center of the valley of reduced variation.     

Hitchhiking effect of a beneficial mutation spreading in a subdivided population

A central problem in population genetics is to detect and analyze positive natural selection by which beneficial mutations are driven to fixation. The hitchhiking effect of a rapidly spreading beneficial mutation, which results in local removal of standing genetic variation, allows such an analysis using DNA sequence polymorphism. However, the current mathematical theory that predicts the pattern of genetic hitchhiking relies on the assumption that a beneficial mutation increases to a high frequency in a single random-mating population, which is certainly violated in reality. Individuals in natural populations are distributed over a geographic space. The spread of a beneficial allele can be delayed by limited migration of individuals over the space and its hitchhiking effect can also be affected. To study this effect of geographic structure on genetic hitchhiking, we analyze a simple model of directional selection in a subdivided population. In contrast to previous studies on hitchhiking in subdivided populations, we mainly investigate the range of sufficiently high migration rates that would homogenize genetic variation at neutral loci. We provide a heuristic mathematical analysis that describes how the genealogical structure at a neutral locus linked to the locus under selection is expected to change in a population divided into two demes. Our results indicate that the overall strength of genetic hitchhiking--the degree to which expected heterozygosity decreases--is diminished by population subdivision, mainly because opportunity for the breakdown of hitchhiking by recombination increases as the spread of the beneficial mutation across demes is delayed when migration rate is much smaller than the strength of selection. Furthermore, the amount of genetic variation after a selective sweep is expected to be unequal over demes: a greater reduction in expected heterozygosity occurs in the subpopulation from which the beneficial mutation originates than in its neighboring subpopulations. This raises a possibility of detecting a "hidden" geographic structure of population by carefully analyzing the pattern of a selective sweep.     

A scan of molecular variation leads to the narrow localization of a selective sweep affecting both Afrotropical and cosmopolitan populations of Drosophila melanogaster

Drosophila melanogaster originated in tropical Africa but has achieved a cosmopolitan distribution in association with human habitation. Cosmopolitan populations of D. melanogaster are known to have reduced genetic variation, particularly on the X chromosome. However, the relative importance of population bottlenecks and selective sweeps in explaining this reduction is uncertain. We surveyed variation at 31 microsatellites across a 330-kb section of the X chromosome located between the white and kirre genes. Two linked clusters of loci were observed with reduced variation and a skew toward rare alleles in both an Ecuador and a Zimbabwe population sample. Examining Zimbabwe DNA sequence polymorphism within one of these regions allowed us to localize a selective sweep to a 361-bp window within the 5' regulatory region of the roughest gene, with one nucleotide substitution representing the best candidate for the target of selection. Estimates of sweep age suggested that this fixation event occurred prior to the expansion of D. melanogaster from sub-Saharan Africa. For both putative sweep regions in our data set, cosmopolitan populations showed wider footprints of selection compared to those in Zimbabwe. This pattern appears consistent with the demographic amplification of preexisting sweep signals due to one or more population bottlenecks.     

Reduced variation around drug-resistant dhfr alleles in African Plasmodium falciparum

We have measured microsatellite diversity at 26 markers around the dhfr gene in pyrimethamine-sensitive and -resistant parasites collected in southeast Africa. Through direct comparison with diversity on sensitive chromosomes we have found significant loss of diversity across a region of 70 kb around the most highly resistant allele which is evidence of a selective sweep attributable to selection through widespread use of pyrimethamine (in combination with sulfadoxine) as treatment for malaria. Retrospective analysis through four years of direct and continuous selection from use of sulfadoxine-pyrimethamine as first-line malaria treatment on a Plasmodium falciparum population in KwaZulu Natal, South Africa, has revealed how recombination significantly narrowed the margins of the selective sweep over time. A deterministic model incorporating selection coefficients measured during the same interval indicates that the transition was toward a state of recombination-selection equilibrium. We compared loss of diversity around the same resistance allele in two populations at either extreme of the range of entomological inoculation rates (EIRs), namely, under one infective bite per year in Mpumalanga, South Africa, and more than one per day in southern Tanzania. EIRs determine effective recombination rates and are expected to profoundly influence the dimensions of the selective sweep. Surprisingly, the dimensions were broadly consistent across both populations. We conclude that despite different recombination rates and contrasting drug selection histories in neighboring countries, the region-wide movement of resistant parasites has played a key role in the establishment of resistance in these populations and the dimensions of the selective sweep are dominated by the influence of high initial starting frequencies.     

Natural selection and population history in the human angiotensinogen gene (AGT): 736 complete AGT sequences in chromosomes from around the world

Several lines of evidence suggest that patterns of genetic variability in the human angiotensinogen gene (AGT) contribute to phenotypic variability in human hypertension. The A(-6) promoter variant of AGT is associated with higher plasma angiotensinogen levels and increased risk of essential hypertension. The geographic distribution of the A(-6) variant leads to the intriguing hypothesis that the G(-6) promoter variant has been selectively advantageous outside Africa. To test these hypotheses, we investigated the roles of population history and natural selection in shaping patterns of genetic diversity in AGT, by sequencing the entire AGT gene (14400 bp) in 736 chromosomes from Africa, Asia, and Europe. We found that the A(-6) variant is present at higher frequency in African populations than in non-African populations. Neutrality tests found no evidence of a departure from selective neutrality, when whole AGT sequences were compared. However, tests restricted to sites in the vicinity of the A(-6)G polymorphism found evidence of a selective sweep. Sliding-window analyses showed that evidence of the sweep is restricted to sites in tight linkage disequilibrium (LD) with the A(-6)G polymorphism. Further, haplotypes carrying the G(-6) variant showed elevated levels of LD, suggesting that they have risen recently to high frequency. Departures from neutral expectation in some but not all regions of AGT indicate that patterns of diversity in the gene cannot be accounted for solely by population history, which would affect all regions equally. Taken together, patterns of genetic diversity in AGT suggest that natural selection has generally favored the G(-6) variant over the A(-6) variant in non-African populations. However, important localized effects may also be present.     

Soft Selective Sweeps in Complex Demographic Scenarios

Adaptation from de novo mutation can produce so-called soft selective sweeps, where adaptive alleles of independent mutational origin sweep through the population at the same time. Population genetic theory predicts that such soft sweeps should be likely if the product of the population size and the mutation rate toward the adaptive allele is sufficiently large, such that multiple adaptive mutations can establish before one has reached fixation; however, it remains unclear how demographic processes affect the probability of observing soft sweeps. Here we extend the theory of soft selective sweeps to realistic demographic scenarios that allow for changes in population size over time. We first show that population bottlenecks can lead to the removal of all but one adaptive lineage from an initially soft selective sweep. The parameter regime under which such “hardening” of soft selective sweeps is likely is determined by a simple heuristic condition. We further develop a generalized analytical framework, based on an extension of the coalescent process, for calculating the probability of soft sweeps under arbitrary demographic scenarios. Two important limits emerge within this analytical framework: In the limit where population-size fluctuations are fast compared to the duration of the sweep, the likelihood of soft sweeps is determined by the harmonic mean of the variance effective population size estimated over the duration of the sweep; in the opposing slow fluctuation limit, the likelihood of soft sweeps is determined by the instantaneous variance effective population size at the onset of the sweep. We show that as a consequence of this finding the probability of observing soft sweeps becomes a function of the strength of selection. Specifically, in species with sharply fluctuating population size, strong selection is more likely to produce soft sweeps than weak selection. Our results highlight the importance of accurate demographic estimates over short evolutionary timescales for understanding the population genetics of adaptation from de novo mutation.     

Selective sweeps of mitochondrial DNA can drive the evolution of uniparental inheritance

Although the uniparental (or maternal) inheritance of mitochondrial DNA (mtDNA) is widespread, the reasons for its evolution remain unclear. Two main hypotheses have been proposed: selection against individuals containing different mtDNAs (heteroplasmy) and selection against “selfish” mtDNA mutations. Recently, uniparental inheritance was shown to promote adaptive evolution in mtDNA, potentially providing a third hypothesis for its evolution. Here, we explore this hypothesis theoretically and ask if the accumulation of beneficial mutations provides a sufficient fitness advantage for uniparental inheritance to invade a population in which mtDNA is inherited biparentally. In a deterministic model, uniparental inheritance increases in frequency but cannot replace biparental inheritance if only a single beneficial mtDNA mutation sweeps through the population. When we allow successive selective sweeps of mtDNA, however, uniparental inheritance can replace biparental inheritance. Using a stochastic model, we show that a combination of selection and drift facilitates the fixation of uniparental inheritance (compared to a neutral trait) when there is only a single selective mtDNA sweep. When we consider multiple mtDNA sweeps in a stochastic model, uniparental inheritance becomes even more likely to replace biparental inheritance. Our findings thus suggest that selective sweeps of beneficial mtDNA haplotypes can drive the evolution of uniparental inheritance.     

The epidemiology of Huntington's disease

Summary The available information on the world distribution of Huntington's disease (HD) from population surveys and death rate analysis is summarised and discussed in the light of genetic studies. It is concluded that most European populations, both Northern and Southern, show a relatively high prevalence (4–8 per 100,000), and that the disorder may also be frequent in India and parts of central Asia. HD is notably rare in Finland and in Japan, but data for Eastern Asia and Africa are inadequate. The disorder may have been underestimated in the American black population. Populations derived from recent European imigration show frequencies and origins of HD comparable to those expected from their own origins and expansion; there is no evidence to suggest that the HD gene has spread disproportionally and its selective effect may be close to neutral. Multiple separate introductions of the gene have been the rule in large populations. Several major foci of HD exist as the result of rapid population expansion. It is likely that a number of separate mutations for HD will be shown to be responsible for the disease, but that the high frequency of HD in European populations will prove to be the result of one or a very small number of mutations, probably of great antiquity.     

The fixation probability of two competing beneficial mutations

Suppose that a beneficial mutation is undergoing a selective sweep when another beneficial mutation arises at a linked locus. We study the fixation probability of the double mutant, i.e., one (produced by recombination) that carries both mutations. Previous analysis works well for the case where the earlier beneficial mutation confers a greater selective advantage than the later mutation, but not so well in the opposite case. We present an approach to approximating the fixation probability in the case where the later mutation confers a greater selective advantage.     

Soft sweeps: molecular population genetics of adaptation from standing genetic variation

A population can adapt to a rapid environmental change or habitat expansion in two ways. It may adapt either through new beneficial mutations that subsequently sweep through the population or by using alleles from the standing genetic variation. We use diffusion theory to calculate the probabilities for selective adaptations and find a large increase in the fixation probability for weak substitutions, if alleles originate from the standing genetic variation. We then determine the parameter regions where each scenario-standing variation vs. new mutations-is more likely. Adaptations from the standing genetic variation are favored if either the selective advantage is weak or the selection coefficient and the mutation rate are both high. Finally, we analyze the probability of "soft sweeps," where multiple copies of the selected allele contribute to a substitution, and discuss the consequences for the footprint of selection on linked neutral variation. We find that soft sweeps with weaker selective footprints are likely under both scenarios if the mutation rate and/or the selection coefficient is high.     

When selection deceives phylogeographic interpretation: the case of the Mediterranean house gecko, Hemidactylus turcicus (Linnaeus, 1758)

A previous study on Hemidactylus turcicus based on mtDNA makers indicated that this gecko has a Middle-East origin, and that the current phylogeographic pattern is the result of a very rapid spread from the east to the west of the species’ range. The same study identified two distinct mitochondrial lineages with low differentiation and genetic diversity. Since H. turcicus is known to be closely associated to humanized environments, its present distribution range and phylogeography is frequently interpreted to be the result of recurrent human-mediated introductions. These conclusions used to be the same as those used to interpret the results obtained for the European populations of another gecko, Tarentola mauritanica. However, a recent study has revealed that the phylogeographic pattern of T. mauritanica is not solely the result of a recent colonization, but also of a mitochondrial selective sweep. Could the same be occurring in H. turcicus? To answer this question, two mitochondrial (12S rRNA and cytochrome b) and two nuclear genes (ACM4 and Rag2) were used in this study. From the mtDNA data we confirmed the existence of two distinct phylogeographic lineages; one occurring exclusively in the northern Mediterranean (Clade A), and another one more widespread that is the only lineage present in North Africa (Clade B). In light of these results, we could hypothesize that H. turcicus had its origin in Turkey, and from there Clade A moved to Europe and Clade B to North Africa spreading latter into Europe. However, Clade A presents significantly higher nucleotide diversity for the nuclear DNA compared to the mtDNA, and neutrality tests gave significant results for the mitochondrial data. These results suggest that the lack of mtDNA genetic diversity and structure in the European population of H. turcicus could also be due to a selective sweep, and not only because of a recent colonization. Together with the situation reported in T. mauritanica, the identification of a hitch-hiking process occurring in H. turcicus, represents two unprecedented cases of a selective sweep taking place in the same geographic area shaping the phylogeographic patterns of two unrelated genera of geckos.     

Estimating the strength of selective sweeps from deep population diversity data

Selective sweeps are typically associated with a local reduction of genetic diversity around the adaptive site. However, selective sweeps can also quickly carry neutral mutations to observable population frequencies if they arise early in a sweep and hitchhike with the adaptive allele. We show that the interplay between mutation and exponential amplification through hitchhiking results in a characteristic frequency spectrum of the resulting novel haplotype variation that depends only on the ratio of the mutation rate and the selection coefficient of the sweep. On the basis of this result, we develop an estimator for the selection coefficient driving a sweep. Since this estimator utilizes the novel variation arising from mutations during a sweep, it does not rely on preexisting variation and can also be applied to loci that lack recombination. Compared with standard approaches that infer selection coefficients from the size of dips in genetic diversity around the adaptive site, our estimator requires much shorter sequences but sampled at high population depth to capture low-frequency variants; given such data, it consistently outperforms standard approaches. We investigate analytically and numerically how the accuracy of our estimator is affected by the decay of the sweep pattern over time as a consequence of random genetic drift and discuss potential effects of recombination, soft sweeps, and demography. As an example for its use, we apply our estimator to deep sequencing data from human immunodeficiency virus populations.