The Coalescent Process in Models with Selection and Recombination

The statistical properties of the process describing the genealogical history of a random sample of genes at a selectively neutral locus which is linked to a locus at which natural selection operates are investigated. It is found that the equations describing this process are simple modifications of the equations describing the process assuming that the two loci are completely linked. Thus, the statistical properties of the genealogical process for a random sample at a neutral locus linked to a locus with selection follow from the results obtained for the selected locus. Sequence data from the alcohol dehydrogenase (Adh) region of Drosophila melanogaster are examined and compared to predictions based on the theory. It is found that the spatial distribution of nucleotide differences between Fast and Slow alleles of Adh is very similar to the spatial distribution predicted if balancing selection operates to maintain the allozyme variation at the Adh locus. The spatial distribution of nucleotide differences between different Slow alleles of Adh do not match the predictions of this simple model very well.     

The hitchhiking effect on linkage disequilibrium between linked neutral loci

We analyzed a three-locus model of genetic hitchhiking with one locus experiencing positive directional selection and two partially linked neutral loci. Following the original hitchhiking approach by Maynard Smith and Haigh, our analysis is purely deterministic. In the first half of the selected phase after a favored mutation has entered the population, hitchhiking may lead to a strong increase of linkage disequilibrium (LD) between the two neutral sites if both are <0.1 s away from the selected site (where s is the selection coefficient). In the second half of the selected phase, the main effect of hitchhiking is to destroy LD. This occurs very quickly (before the end of the selected phase) when the selected site is between both neutral loci. This pattern cannot be attributed to the well-known variation-reducing effect of hitchhiking but is a consequence of secondary hitchhiking effects on the recombinants created in the selected phase. When the selected site is outside the neutral loci (which are, say, <0.1s apart), however, a fast decay of LD is observed only if the selected site is in the immediate neighborhood of one of the neutral sites (i.e., if the recombination rate r between the selected site and one of the neutral sites satisfies r<0.1 s). If the selected site is far away from the neutral sites (say, r > 0.3 s), the decay rate of LD approaches that of neutrality. Averaging over a uniform distribution of initial gamete frequencies shows that the expected LD at the end of the hitchhiking phase is driven toward zero, while the variance is increased when the selected site is well outside the two neutral sites. When the direction of LD is polarized with respect to the more common allele at each neutral site, hitchhiking creates more positive than negative linkage disequilibrium. Thus, hitchhiking may have a distinctively patterned LD-reducing effect, in particular near the target of selection.     

Genetic drift in an infinite population. The pseudohitchhiking model

Selected substitutions at one locus can induce stochastic dynamics that resemble genetic drift at a closely linked neutral locus. The pseudohitchhiking model is a one-locus model that approximates these effects and can be used to describe the major consequences of linked selection. As the changes in neutral allele frequencies when hitchhiking are rapid, diffusion theory is not appropriate for studying neutral dynamics. A stationary distribution and some results on substitution processes are presented that use the theory of continuous-time Markov processes with discontinuous sample paths. The coalescent of the pseudohitchhiking model is shown to have a random number of branches at each node, which leads to a frequency spectrum that is different from that of the equilibrium neutral model. If genetic draft, the name given to these induced stochastic effects, is a more important stochastic force than genetic drift, then a number of paradoxes that have plagued population genetics disappear.     

The Effect of a Selected Locus on Linked Neutral Loci

The effects produced on linked neutral loci as a selected locus evolves towards its equilibrium value are considered. Significant effects on the neutral loci arise if the recombination fraction between the neutral and selected loci is smaller than the order of magnitude of the selective differences at the selected locus. The effect on gene frequencies at the neutral loci, that is, the hitchhiking effect, is determined, as well as the linkage disequilibrium generated by this hitchhiking effect. One of the more important findings is that significant disequilibrium can be generated between two neutral loci by the evolution of a linked selected locus. Consideration is given to the problem of determining how the effect of selection operating in natural populations can be detected, the question of the establishment of inversions in populations, and also to the nonequilibrium properties of populations.     

Evidence for genetic hitchhiking effect associated with insecticide resistance in Aedes aegypti.

Information on genetic variation within and between populations is critical for understanding the evolutionary history of mosquito populations and disease epidemiology. Previous studies with Drosophila suggest that genetic variation of selectively neutral loci in a large fraction of genome may be constrained by fixation of advantageous mutations associated with hitchhiking effect. This study examined restriction fragment length polymorphisms of four natural Aedes aegypti mosquito populations from Trinidad and Tobago, at 16 loci. These populations have been subjected to organophosphate (OP) insecticide treatments for more than two decades, while dichlor-diphenyltrichlor (DDT) was the insecticide of choice prior to this period. We predicted that genes closely linked to the OP target loci would exhibit reduced genetic variation as a result of the hitchhiking effect associated with intensive OP insecticide selection. We also predicted that genetic variability of the genes conferring resistance to DDT and loci near the target site would be similar to other unlinked loci. As predicted, reduced genetic variation was found for loci in the general chromosomal region of a putative OP target site, and these loci generally exhibited larger F(ST) values than other random loci. In contrast, the gene conferring resistance to DDT and its linked loci show polymorphisms and genetic differentiation similar to other random loci. The reduced genetic variability and apparent gene deletion in some regions of chromosome 1 likely reflect the hitchhiking effect associated with OP insecticide selection.     

Linkage disequilibrium under genetic hitchhiking in finite populations

The model of genetic hitchhiking predicts a reduction in sequence diversity at a neutral locus closely linked to a beneficial allele. In addition, it has been shown that the same process results in a specific pattern of correlations (linkage disequilibrium) between neutral polymorphisms along the chromosome at the time of fixation of the beneficial allele. During the hitchhiking event, linkage disequilibrium on either side of the beneficial allele is built up whereas it is destroyed across the selected site. We derive explicit formulas for the expectation of the covariance measure D and standardized linkage disequilibrium sigma 2D between a pair of polymorphic sites. For our analysis we use the approximation of a star-like genealogy at the selected site. The resulting expressions are approximately correct in the limit of large selection coefficients. Using simulations we show that the resulting pattern of linkage disequilibrium is quickly-i.e., in <0.1N generations-destroyed after the fixation of the beneficial allele for moderately distant neutral loci, where N is the diploid population size.     

THE EFFICACY OF DIVERGENCE HITCHHIKING IN GENERATING GENOMIC ISLANDS DURING ECOLOGICAL SPECIATION

Genes under divergent selection flow less readily between populations than other loci. This observation has led to verbal “divergence hitchhiking” models of speciation in which decreased interpopulation gene flow surrounding loci under divergent selection can generate large regions of differentiation within the genome (genomic islands). The efficacy of this model in promoting speciation depends on the size of the region affected by divergence hitchhiking. Empirical evidence is mixed, with examples of both large and small genomic islands. To address these empirical discrepancies and to formalize the theory, we present mathematical models of divergence hitchhiking, which examine neutral differentiation around selected sites. For a single locus under selection, regions of differentiation do not extend far along a chromosome away from a selected site unless both effective population sizes and migration rates are low. When multiple loci are considered, regions of differentiation can be larger. However, with many loci under selection, genome‐wide divergence occurs and genomic islands are erased. The results show that divergence hitchhiking can generate large regions of differentiation, but that the conditions under which this occurs are limited. Thus, speciation may often require multifarious selection acting on many, isolated and physically unlinked genes. How hitchhiking promotes further adaptive divergence warrants consideration.     

The Hitchhiking Effect on the Site Frequency Spectrum of DNA Polymorphisms

The level of DNA sequence variation is reduced in regions of the Drosophila melanogaster genome where the rate of crossing over per physical distance is also reduced. This observation has been interpreted as support for the simple model of genetic hitchhiking, in which directional selection on rare variants, e.g., newly arising advantageous mutants, sweeps linked neutral alleles to fixation, thus eliminating polymorphisms near the selected site. However, the frequency spectra of segregating sites of several loci from some populations exhibiting reduced levels of nucleotide diversity and reduced numbers of segregating sites did not appear different from what would be expected under a neutral equilibrium model. Specifically, a skew toward an excess of rare sites was not observed in these samples, as measured by Tajima's D. Because this skew was predicted by a simple hitchhiking model, yet it had never been expressed quantitatively and compared directly to DNA polymorphism data, this paper investigates the hitchhiking effect on the site frequency spectrum, as measured by Tajima's D and several other statistics, using a computer simulation model based on the coalescent process and recurrent hitchhiking events. The results presented here demonstrate that under the simple hitchhiking model (1) the expected value of Tajima's D is large and negative (indicating a skew toward rare variants), (2) that Tajima's test has reasonable power to detect a skew in the frequency spectrum for parameters comparable to those from actual data sets, and (3) that the Tajima's Ds observed in several data sets are very unlikely to have been the result of simple hitchhiking. Consequently, the simple hitchhiking model is not a sufficient explanation for the DNA polymorphism at those loci exhibiting a decreased number of segregating sites yet not exhibiting a skew in the frequency spectrum.     

A regression-based approach to selection mapping

Selection mapping applies the population genetics theory of hitchhiking to the localization of genomic regions containing genes under selection. This approach predicts that neutral loci linked to genes under positive selection will have reduced diversity due to their shared history with a selected locus, and thus, genome scans of diversity levels can be used to identify regions containing selected loci. Most previous approaches to this problem ignore the spatial genomic pattern of diversity expected under selection. The regression-based approach advocated in this paper takes into account the expected pattern of decreasing genetic diversity with increased proximity to a selected locus. Simulated data are used to examine the patterns of diversity under different scenarios, in order to assess the power of a regression-based approach to the identification of regions under selection. Application of this method to both simulated and empirical data demonstrates its potential to detect selection. In contrast to some other methods, the regression approach described in this paper can be applied to any marker type. Results also suggest that this approach may give more precise estimates of the location of the selected locus than alternative methods, although the power is slightly lower in some cases.     

The Coalescent Process in Models with Selection

Statistical properties of the process describing the genealogical history of a random sample of genes are obtained for a class of population genetics models with selection. For models with selection, in contrast to models without selection, the distribution of this process, the coalescent process, depends on the distribution of the frequencies of alleles in the ancestral generations. If the ancestral frequency process can be approximated by a diffusion, then the mean and the variance of the number of segregating sites due to selectively neutral mutations in random samples can be numerically calculated. The calculations are greatly simplified if the frequencies of the alleles are tightly regulated. If the mutation rates between alleles maintained by balancing selection are low, then the number of selectively neutral segregating sites in a random sample of genes is expected to substantially exceed the number predicted under a neutral model.     

Linkage disequilibria and the site frequency spectra in the su(s) and su(w(a)) regions of the Drosophila melanogaster X chromosome

Over the last decade, surveys of DNA sequence variation in natural populations of several Drosophila species and other taxa have established that polymorphism is reduced in genomic regions characterized by low rates of crossing over per physical length. Parallel studies have also established that divergence between species is not reduced in these same genomic regions, thus eliminating explanations that rely on a correlation between the rates of mutation and crossing over. Several theoretical models (directional hitchhiking, background selection, and random environment) have been proposed as population genetic explanations. In this study samples from an African population (n = 50) and a European population (n = 51) were surveyed at the su(s) (1955 bp) and su(w(a)) (3213 bp) loci for DNA sequence polymorphism, utilizing a stratified SSCP/DNA sequencing protocol. These loci are located near the telomere of the X chromosome, in a region of reduced crossing over per physical length, and exhibit a significant reduction in DNA sequence polymorphism. Unlike most previously surveyed, these loci reveal substantial skews toward rare site frequencies, consistent with the predictions of directional hitchhiking and random environment models and inconsistent with the general predictions of the background selection model (or neutral theory). No evidence for excess geographic differentiation at these loci is observed. Although linkage disequilibrium is observed between closely linked sites within these loci, many recombination events in the genealogy of the sampled alleles can be inferred and the genomic scale of linkage disequilibrium, measured in base pairs between sites, is the same as that observed for loci in regions of normal crossing over. We conclude that gene conversion must be high in these regions of low crossing over.     

Effect of strong directional selection on weakly selected mutations at linked sites: implication for synonymous codon usage

The fixation of weakly selected mutations can be greatly influenced by strong directional selection at linked loci. Here, I investigate a two-locus model in which weakly selected, reversible mutations occur at one locus and recurrent strong directional selection occurs at the other locus. This model is analogous to selection on codon usage at synonymous sites linked to nonsynonymous sites under strong directional selection. Two approximations obtained here describe the expected frequency of the weakly selected preferred alleles at equilibrium. These approximations, as well as simulation results, show that the level of codon bias declines with an increasing rate of substitution at the strongly selected locus, as expected from the well-understood theory that selection at one locus reduces the efficacy of selection at linked loci. These solutions are used to examine whether the negative correlation between codon bias and nonsynonymous substitution rates recently observed in Drosophila can be explained by this hitchhiking effect. It is shown that this observation can be reasonably well accounted for if a large fraction of the nonsynonymous substitutions on genes in the data set are driven by strong directional selection.     

Joint effects of genetic hitchhiking and background selection on neutral variation

Due to relatively high rates of strongly selected deleterious mutations, directional selection on favorable alleles (causing hitchhiking effects on linked neutral polymorphisms) is expected to occur while a deleterious mutation-selection balance is present in a population. We analyze this interaction of directional selection and background selection and study their combined effects on neutral variation, using a three-locus model in which each locus is subjected to either deleterious, favorable, or neutral mutations. Average heterozygosity is measured by simulations (1) at the stationary state under the assumption of recurrent hitchhiking events and (2) as a transient level after a single hitchhiking event. The simulation results are compared to theoretical predictions. It is shown that known analytical solutions describing the hitchhiking effect without background selection can be modified such that they accurately predict the joint effects of hitchhiking and background on linked, neutral variation. Generalization of these results to a more appropriate multilocus model (such that background selection can occur at multiple sites) suggests that, in regions of very low recombination rates, stationary levels of nucleotide diversity are primarily determined by hitchhiking, whereas in regions of high recombination, background selection is the dominant force. The implications of these results on the identification and estimation of the relevant parameters of the model are discussed.     

DNA sequence polymorphism and divergence at the erect wing and suppressor of sable loci of Drosophila melanogaster and D. simulans

Several evolutionary models of linked selection (e.g., genetic hitchhiking, background selection, and random environment) predict a reduction in polymorphism relative to divergence in genomic regions where the rate of crossing over per physical distance is restricted. We tested this prediction near the telomere of the Drosophila melanogaster and D. simulans X chromosome at two loci, erect wing (ewg) and suppressor of sable [su(s)]. Consistent with this prediction, polymorphism is reduced at both loci, while divergence is normal. The reduction is greater at ewg, the more distal of the two regions. Two models can be discriminated by comparing the observed site frequency spectra with those predicted by the models. The hitchhiking model predicts a skew toward rare variants in a sample, while the spectra under the background-selection model are similar to those of the neutral model of molecular evolution. Statistical tests of the fit to the predictions of these models require many sampled alleles and segregating sites. Thus we used SSCP and stratified DNA sequencing to cover a large number of randomly sampled alleles (approximately 50) from each of three populations. The result is a clear trend toward negative values of Tajima's D, indicating an excess of rare variants at ewg, the more distal of the two loci. One fixed difference among the populations and high FST values indicate strong population subdivision among the three populations at ewg. These results indicate genetic hitchhiking at ewg, in particular, geographically localized hitchhiking events within Africa. The reduction of polymorphism at su(s) combined with the excess of high-frequency variants in D. simulans is inconsistent with the hitchhiking and background-selection models.     

Natural selection at linked sites in humans

Theoretical and empirical work indicates that patterns of neutral polymorphism can be affected by linked, selected mutations. Under background selection, deleterious mutations removed from a population by purifying selection cause a reduction in linked neutral diversity. Under genetic hitchhiking, the rise in frequency and fixation of beneficial mutations also reduces the level of linked neutral polymorphism. Here we review the evidence that levels of neutral polymorphism in humans are affected by selection at linked sites. We then discuss four approaches for distinguishing between background selection and genetic hitchhiking based on (i) the relationship between polymorphism level and recombination rate for neutral loci with high mutation rates, (ii) relative levels of variation on the X chromosome and the autosomes, (iii) the frequency distribution of neutral polymorphisms, and (iv) population-specific patterns of genetic variation. Although the evidence for selection at linked sites in humans is clear, current methods and data do not allow us to clearly assess the relative importance of background selection and genetic hitchhiking in humans. These results contrast with those obtained for Drosophila, where the signals of positive selection are stronger.     

Analysis of a genetic hitchhiking model, and its application to DNA polymorphism data from Drosophila melanogaster

Begun and Aquadro have demonstrated that levels of nucleotide variation correlate with recombination rate among 20 gene regions from across the genome of Drosophila melanogaster. It has been suggested that this correlation results from genetic hitchhiking associated with the fixation of strongly selected mutants. The hitchhiking process can be described as a series of two-step events. The first step consists of a strongly selected substitution wiping out linked variation in a population; this is followed by a recovery period in which polymorphism can build up via neutral mutations and random genetic drift. Genetic hitchhiking has previously been modeled as a steady-state process driven by recurring selected substitutions. We show here that the characteristic parameter of this steady-state model is alpha v, the product of selection intensity (alpha = 2Ns) and the frequency of beneficial mutations v (where N is population size and s is the selective advantage of the favored allele). We also demonstrate that the steady-state model describes the hitchhiking process adequately, unless the recombination rate is very low. To estimate alpha v, we use the data of DNA sequence variation from 17 D. melanogaster loci from regions of intermediate to high recombination rates. We find that alpha v is likely to be > 1.3 x 10(-8). Additional data are needed to estimate this parameter more precisely. The estimation of alpha v is important, as this parameter determines the shape of the frequency distribution of strongly selected substitutions.      

Regions of lower crossing over harbor more rare variants in African populations of Drosophila melanogaster

A correlation between diversity levels and rates of recombination is predicted both by models of positive selection, such as hitchhiking associated with the rapid fixation of advantageous mutations, and by models of purifying selection against strongly deleterious mutations (commonly referred to as "background selection"). With parameter values appropriate for Drosophila populations, only the first class of models predicts a marked skew in the frequency spectrum of linked neutral variants, relative to a neutral model. Here, we consider 29 loci scattered throughout the Drosophila melanogaster genome. We show that, in African populations, a summary of the frequency spectrum of polymorphic mutations is positively correlated with the meiotic rate of crossing over. This pattern is demonstrated to be unlikely under a model of background selection. Models of weakly deleterious selection are not expected to produce both the observed correlation and the extent to which nucleotide diversity is reduced in regions of low (but nonzero) recombination. Thus, of existing models, hitchhiking due to the recurrent fixation of advantageous variants is the most plausible explanation for the data.     

Molecular evolution near a two-locus balanced polymorphism

Balancing selection at one locus can increase the amount of selectively neutral variation within neighboring genomic regions. Discrete phenotypic polymorphisms studied in natural populations are frequently determined by sets of interacting genes instead of alternative alleles at single loci. We extend coalescent theory to investigate balancing selection on combinations of linked genes. We find that variation at neutral sites is increased across a much larger genomic region relative to the single-locus models: the entire region lying between the two loci in balanced combination is affected to some degree. Epistatic selection maintains these high levels of neutral variation because it directly opposes the homogenizing effect of recombination. The results of the theory are discussed in relation to published gene sequence data, primarily from Drosophila.     

Background selection and population differentiation

A general analytical formula is derived, which predicts the effects of background selection on population differentiation at a neutral locus as a result of its linkage with selected loci of deleterious mutations. The theory is based on the assumptions of random mating, multiplicative fitness, and weak selection in hermaphrodite plants in the island model of population structure. The analytical results show that Fst at the neutral locus increases as a result of the effects of background selection, regardless of the dependence or independence among linked background selective loci. The increment in Fst is closely related to the magnitude of linkage disequilibria between the neutral locus and selected loci, and can be estimated by the ratio of Fst with background selection to Fst without background selection minus one. The steady-state linkage disequilibrium between a neutral locus and a selected locus in subpopulations, primarily attained by gene flow, decreases with the recombination rate, and can be enhanced when there are dependence among linked selected loci. Monte Carlo computer simulations with two- and three-locus models show that the analytical formulae perform well under general conditions. Application of the present theory may aid in analyzing the genome-wide mapping of the effect of background selection in terms of Fst.