A genome-wide departure from the standard neutral model in natural populations of Drosophila

We analyze nucleotide polymorphism data for a large number of loci in areas of normal to high recombination in Drosophila melanogaster and D. simulans (24 and 16 loci, respectively). We find a genome-wide, systematic departure from the neutral expectation for a panmictic population at equilibrium in natural populations of both species. The distribution of sequence-based estimates of 2Nc across loci is inconsistent with the assumptions of the standard neutral theory, given the observed levels of nucleotide diversity and accepted values for recombination and mutation rates. Under these assumptions, most estimates of 2Nc are severalfold too low; in other words, both species exhibit greater intralocus linkage disequilibrium than expected. Variation in recombination or mutation rates is not sufficient to account for the excess of linkage disequilibrium. While an equilibrium island model does not seem to account for the data, more complicated forms of population structure may. A proper test of alternative demographic models will require loci to be sampled in a more consistent fashion.     

Modelling bacterial speciation

A central problem in understanding bacterial speciation is how clusters of closely related strains emerge and persist in the face of recombination. We use a neutral Fisher-Wright model in which genotypes, defined by the alleles at 140 house-keeping loci, change in each generation by mutation or recombination, and examine conditions in which an initially uniform population gives rise to resolved clusters. Where recombination occurs at equal frequency between all members of the population, we observe a transition between clonal structure and sexual structure as the rate of recombination increases. In the clonal situation, clearly resolved clusters are regularly formed, break up or go extinct. In the sexual situation, the formation of distinct clusters is prevented by the cohesive force of recombination. Where the rate of recombination is a declining log-linear function of the genetic distance between the donor and recipient strain, distinct clusters emerge even with high rates of recombination. These clusters arise in the absence of selection, and have many of the properties of species, with high recombination rates and thus sexual cohesion within clusters and low rates between clusters. Distance-scaled recombination can thus lead to a population splitting into distinct genotypic clusters, a process that mimics sympatric speciation. However, empirical estimates of the relationship between sequence divergence and recombination rate indicate that the decline in recombination is an insufficiently steep function of genetic distance to generate species in nature under neutral drift, and thus that other mechanisms should be invoked to explain speciation in the presence of recombination.     

Estimation of hominoid ancestral population sizes under bayesian coalescent models incorporating mutation rate variation and sequencing errors

Estimation of population parameters for the common ancestors of humans and the great apes is important in understanding our evolutionary history. In particular, inference of population size for the human-chimpanzee common ancestor may shed light on the process by which the 2 species separated and on whether the human population experienced a severe size reduction in its early evolutionary history. In this study, the Bayesian method of ancestral inference of Rannala and Yang (2003. Bayes estimation of species divergence times and ancestral population sizes using DNA sequences from multiple loci. Genetics. 164:1645-1656) was extended to accommodate variable mutation rates among loci and random species-specific sequencing errors. The model was applied to analyze a genome-wide data set of approximately 15,000 neutral loci (7.4 Mb) aligned for human, chimpanzee, gorilla, orangutan, and macaque. We obtained robust and precise estimates for effective population sizes along the hominoid lineage extending back approximately 30 Myr to the cercopithecoid divergence. The results showed that ancestral populations were 5-10 times larger than modern humans along the entire hominoid lineage. The estimates were robust to the priors used and to model assumptions about recombination. The unusually low X chromosome divergence between human and chimpanzee could not be explained by variation in the male mutation bias or by current models of hybridization and introgression. Instead, our parameter estimates were consistent with a simple instantaneous process for human-chimpanzee speciation but showed a major reduction in X chromosome effective population size peculiar to the human-chimpanzee common ancestor, possibly due to selective sweeps on the X prior to separation of the 2 species.     

A neutral explanation for the correlation of diversity with recombination rates in humans

One of the most striking findings to emerge from the study of genomic patterns of variation is that regions with lower recombination rates tend to have lower levels of intraspecific diversity but not of interspecies divergence. This uncoupling of variation within and between species has been widely interpreted as evidence that natural selection shapes patterns of genetic variability genomewide. We revisited the relationship between diversity, divergence, and recombination in humans, using data from closely related species and better estimates of recombination rates than previously available. We show that regions that experience less recombination have reduced divergence to chimpanzee and to baboon, as well as lower levels of diversity. This observation suggests that mutation and recombination are associated processes in humans, so that the positive correlation between diversity and recombination may have a purely neutral explanation. Consistent with this hypothesis, diversity levels no longer increase significantly with recombination rates after correction for divergence to chimpanzee.     

Equilibrium processes cannot explain high levels of short- and medium-range linkage disequilibrium in the domesticated grass Sorghum bicolor

Patterns of linkage disequilibrium (LD) are of interest because they provide evidence of both equilibrium (e.g., mating system or long-term population structure) and nonequilibrium (e.g., demographic or selective) processes, as well as because of their importance in strategies for identifying the genetic basis of complex phenotypes. We report patterns of short and medium range (up to 100 kb) LD in six unlinked genomic regions in the partially selfing domesticated grass, Sorghum bicolor. The extent of allelic associations in S. bicolor, as assessed by pairwise measures of LD, is higher than in maize but lower than in Arabidopsis, in qualitative agreement with expectations based on mating system. Quantitative analyses of the population recombination parameter, rho, however, based on empirical estimates of rates of recombination, mutation, and self-pollination, show that LD is more extensive than expected under a neutral equilibrium model. The disparity between rho and the population mutation parameter, , is similar to that observed in other species whose population history appears to be complex. From a practical standpoint, these results suggest that S. bicolor is well suited for association studies using reasonable numbers of markers, since LD typically extends at least several kilobases but has largely decayed by 15 kb.     

Effects of Mitochondrial DNA Rate Variation on Reconstruction of Pleistocene Demographic History in a Social Avian Species,Pomatostomus superciliosus

Mitochondrial sequence data is often used to reconstruct the demographic history of Pleistocene populations in an effort to understand how species have responded to past climate change events. However, departures from neutral equilibrium conditions can confound evolutionary inference in species with structured populations or those that have experienced periods of population expansion or decline. Selection can affect patterns of mitochondrial DNA variation and variable mutation rates among mitochondrial genes can compromise inferences drawn from single markers. We investigated the contribution of these factors to patterns of mitochondrial variation and estimates of time to most recent common ancestor (T_MRCA) for two clades in a co-operatively breeding avian species, the white-browed babbler Pomatostomus superciliosus. Both the protein-coding ND3 gene and hypervariable domain I control region sequences showed departures from neutral expectations within the superciliosus clade, and a two-fold difference in T_MRCA estimates. Bayesian phylogenetic analysis provided evidence of departure from a strict clock model of molecular evolution in domain I, leading to an over-estimation of T_MRCA for the superciliosus clade at this marker. Our results suggest mitochondrial studies that attempt to reconstruct Pleistocene demographic histories should rigorously evaluate data for departures from neutral equilibrium expectations, including variation in evolutionary rates across multiple markers. Failure to do so can lead to serious errors in the estimation of evolutionary parameters and subsequent demographic inferences concerning the role of climate as a driver of evolutionary change. These effects may be especially pronounced in species with complex social structures occupying heterogeneous environments. We propose that environmentally driven differences in social structure may explain observed differences in evolutionary rate of domain I sequences, resulting from longer than expected retention times for matriarchal lineages in the superciliosus clade.     

Estimating the Relative Roles of Recombination and Point Mutation in the Generation of Single Locus Variants in Campylobacter jejuni and Campylobacter coli

Single locus variants (SLVs) are bacterial sequence types that differ at only one of the seven canonical multilocus sequence typing (MLST) loci. Estimating the relative roles of recombination and point mutation in the generation of new alleles that lead to SLVs is helpful in understanding how organisms evolve. The relative rates of recombination and mutation for Campylobacter jejuni and Campylobacter coli were estimated at seven different housekeeping loci from publically available MLST data. The probability of recombination generating a new allele that leads to an SLV is estimated to be roughly seven times more than that of mutation for C. jejuni, but for C. coli recombination and mutation were estimated to have a similar contribution to the generation of SLVs. The majority of nucleotide differences (98?% for C. jejuni and 85?% for C. coli) between strains that make up an SLV are attributable to recombination. These estimates are much larger than estimates of the relative rate of recombination to mutation calculated from more distantly related isolates using MLST data. One explanation for this is that purifying selection plays an important role in the evolution of Campylobacter. A simulation study was performed to test the performance of our method under a range of biologically realistic parameters. We found that our method performed well when the recombination tract length was longer than 3?kb. For situations in which recombination may occur with shorter tract lengths, our estimates are likely to be an underestimate of the ratio of recombination to mutation, and of the importance of recombination for creating diversity in closely related isolates. A parametric bootstrap method was applied to calculate the uncertainty of these estimates.     

Estimating ancestral population sizes and divergence times

This article presents a new method for jointly estimating species divergence times and ancestral population sizes. The method improves on previous ones by explicitly incorporating intragenic recombination, by utilizing orthologous sequence data from closely related species, and by using a maximum-likelihood framework. The latter allows for efficient use of the available information and provides a way of assessing how much confidence we should place in the estimates. I apply the method to recently collected intergenic sequence data from humans and the great apes. The results suggest that the human-chimpanzee ancestral population size was four to seven times larger than the current human effective population size and that the current human effective population size is slightly >10,000. These estimates are similar to previous ones, and they appear relatively insensitive to assumptions about the recombination rates or mutation rates across loci.     

Frequent beneficial mutations during single-colony serial transfer of Streptococcus pneumoniae

The appearance of new mutations within a population provides the raw material for evolution. The consistent decline in fitness observed in classical mutation accumulation studies has provided support for the long-held view that deleterious mutations are more common than beneficial mutations. Here we present results of a study using a mutation accumulation design with the bacterium Streptococcus pneumoniae in which the fitness of the derived populations increased. This rise in fitness was associated specifically with adaptation to survival during brief stationary phase periods between single-colony population bottlenecks. To understand better the population dynamics behind this unanticipated adaptation, we developed a maximum likelihood model describing the processes of mutation and stationary-phase selection in the context of frequent population bottlenecks. Using this model, we estimate that the rate of beneficial mutations may be as high as 4.8×10(-4) events per genome for each time interval corresponding to the pneumococcal generation time. This rate is several orders of magnitude higher than earlier estimates of beneficial mutation rates in bacteria but supports recent results obtained through the propagation of small populations of Escherichia coli. Our findings indicate that beneficial mutations may be relatively frequent in bacteria and suggest that in S. pneumoniae, which develops natural competence for transformation, a steady supply of such mutations may be available for sampling by recombination.     

The effect of intragenic recombination on the number of alleles in a finite population

A two-site infinite allele model is constructed to study the effect of intragenic recombination on the number of neutral alleles and the distribution of their frequencies in a finite population. The results of theory and Monte Carlo simulation of the two-site model demonstrate that intragenic recombination significantly increases the mean and variance of the number of alleles when the rates of mutation and recombination are as large as the reciprocal of the population size. Data from natural populations indicate that this may be a significant process in generating variation and determining its distribution.     

The Population and Evolutionary Dynamics of Homologous Gene Recombination in Bacteria

In bacteria, recombination is a rare event, not a part of the reproductive process. Nevertheless, recombination—broadly defined to include the acquisition of genes from external sources, i.e., horizontal gene transfer (HGT)—plays a central role as a source of variation for adaptive evolution in many species of bacteria. Much of niche expansion, resistance to antibiotics and other environmental stresses, virulence, and other characteristics that make bacteria interesting and problematic, is achieved through the expression of genes and genetic elements obtained from other populations of bacteria of the same and different species, as well as from eukaryotes and archaea. While recombination of homologous genes among members of the same species has played a central role in the development of the genetics and molecular biology of bacteria, the contribution of homologous gene recombination (HGR) to bacterial evolution is not at all clear. Also, not so clear are the selective pressures responsible for the evolution and maintenance of transformation, the only bacteria-encoded form of HGR. Using a semi-stochastic simulation of mutation, recombination, and selection within bacterial populations and competition between populations, we explore (1) the contribution of HGR to the rate of adaptive evolution in these populations and (2) the conditions under which HGR will provide a bacterial population a selective advantage over non-recombining or more slowly recombining populations. The results of our simulation indicate that, under broad conditions: (1) HGR occurring at rates in the range anticipated for bacteria like Streptococcus pneumoniae, Escherichia coli, Haemophilus influenzae, and Bacillus subtilis will accelerate the rate at which a population adapts to environmental conditions; (2) once established in a population, selection for this capacity to increase rates of adaptive evolution can maintain bacteria-encoded mechanisms of recombination and prevent invasion of non-recombining populations, even when recombination engenders a modest fitness cost; and (3) because of the density- and frequency-dependent nature of HGR in bacteria, this capacity to increase rates of adaptive evolution is not sufficient as a selective force to provide a recombining population a selective advantage when it is rare. Under realistic conditions, homologous gene recombination will increase the rate of adaptive evolution in bacterial populations and, once established, selection for higher rates of evolution will promote the maintenance of bacteria-encoded mechanisms for HGR. On the other hand, increasing rates of adaptive evolution by HGR is unlikely to be the sole or even a dominant selective pressure responsible for the original evolution of transformation.     

A robust measure of HIV-1 population turnover within chronically infected individuals

A simple nonparameteric test for population structure was applied to temporally spaced samples of HIV-1 sequences from the gag-pol region within two chronically infected individuals. The results show that temporal structure can be detected for samples separated by about 22 months or more. The performance of the method, which was originally proposed to detect geographic structure, was tested for temporally spaced samples using neutral coalescent simulations. Simulations showed that the method is robust to variation in samples sizes and mutation rates, to the presence/absence of recombination, and that the power to detect temporal structure is high. By comparing levels of temporal structure in simulations to the levels observed in real data, we estimate the effective intra-individual population size of HIV-1 to be between 10(3) and 10(4) viruses, which is in agreement with some previous estimates. Using this estimate and a simple measure of sequence diversity, we estimate an effective neutral mutation rate of about 5 x 10(-6) per site per generation in the gag-pol region. The definition and interpretation of estimates of such "effective" population parameters are discussed.     

Estimation of the rate and effect of new beneficial mutations in asexual populations

The rate and effect of available beneficial mutations are key parameters in determining how a population adapts to a new environment. However, these parameters are poorly known, in large part because of the difficulty of designing and interpreting experiments to examine the rare and intrinsically stochastic process of mutation occurrence. We present a new approach to estimate the rate and selective advantage of beneficial mutations that underlie the adaptation of asexual populations. We base our approach on the analysis of experiments that track the effect of newly arising beneficial mutations on the dynamics of a neutral marker in evolving bacterial populations and develop efficient estimators of mutation rate and selective advantage. Using extensive simulations, we evaluate the accuracy of our estimators and conclude that they are quite robust to the use of relatively low experimental replication. To validate the predictions of our model, we compare theoretical and experimentally determined estimates of the selective advantage of the first beneficial mutation to fix in a series of ten replicate populations. We find that our theoretical predictions are not significantly different from experimentally determined selection coefficients. Application of our method to suitably designed experiments will allow estimation of how population evolvability depends on demographic and initial fitness parameters.     

Fluctuation Tests: How Reliable Are the Estimates of Mutation Rates?

Fifty one years ago, Luria and Delbruck published in Genetics a paper that was to become a classic. In it they proved, beyond all reasonable doubt, that bacteria were mutating to phage resistance long before they could have encountered any bacteriophage. Luria and Delbruck also showed how the same experimental data could be used to estimate bacterial mutation rates. Since that time and in many different contexts the methods that they introduced have been used to estimate mutation rates. However, little seems to be known about the errors to be expected in such estimates. In what follows I examine how much uncertainty in the estimates is to be expected merely on the basis of the stochastic variability inherent in the sampling process. On the basis of this examination I question a few traditional ideas and conclude with some practical suggestions. The results were obtained by simulation. It is my hope that they may inspire others to provide a rigorous theoretical basis for such calculations.     

Microevolution of Helicobacter pylori during Prolonged Infection of Single Hosts and within Families

Our understanding of basic evolutionary processes in bacteria is still very limited. For example, multiple recent dating estimates are based on a universal inter-species molecular clock rate, but that rate was calibrated using estimates of geological dates that are no longer accepted. We therefore estimated the short-term rates of mutation and recombination in Helicobacter pylori by sequencing an average of 39,300 bp in 78 gene fragments from 97 isolates. These isolates included 34 pairs of sequential samples, which were sampled at intervals of 0.25 to 10.2 years. They also included single isolates from 29 individuals (average age: 45 years) from 10 families. The accumulation of sequence diversity increased with time of separation in a clock-like manner in the sequential isolates. We used Approximate Bayesian Computation to estimate the rates of mutation, recombination, mean length of recombination tracts, and average diversity in those tracts. The estimates indicate that the short-term mutation rate is 1.4×10⁻⁶ (serial isolates) to 4.5×10⁻⁶ (family isolates) per nucleotide per year and that three times as many substitutions are introduced by recombination as by mutation. The long-term mutation rate over millennia is 5–17-fold lower, partly due to the removal of non-synonymous mutations due to purifying selection. Comparisons with the recent literature show that short-term mutation rates vary dramatically in different bacterial species and can span a range of several orders of magnitude.     

Variation in estimated recombination rates across human populations

Recently it has been reported that recombination hotspots appear to be highly variable between humans and chimpanzees, and there is evidence for between-person variability in hotspots, and evolutionary transience. To understand the nature of variation in human recombination rates, it is important to describe patterns of variability across populations. Direct measurement of recombination rates remains infeasible on a large scale, and population-genetic approaches can be imprecise, and are affected by demographic history. Reports to date have suggested broad similarity in recombination rates at large genomic scales and across human populations. Here, we examine recombination rate estimates at a finer population and genomic scale: 28 worldwide populations and 107 SNPs in a 1 Mb stretch of chromosome 22q. We employ analysis of variance of recombination rate estimates, corrected for differences in effective population size using genome-wide microsatellite mutation rate estimates. We find substantial variation in fine-scale rates between populations, but reduced variation within continental groups. All effects examined (SNP-pair, region, population and interactions) were highly significant. Adjustment for effective population size made little difference to the conclusions. Observed hotspots tended to be conserved across populations, albeit at varying intensities. This holds particularly for populations from the same region, and also to a considerable degree across geographical regions. However, some hotspots appear to be population-specific. Several results from studies on the population history of humans are in accordance with our analysis. Our results suggest that between-population variation in DNA sequences may underly recombination rate variation.     

Periodic Selection, Infectious Gene Exchange and the Genetic Structure of E. COLI Populations

As a consequence of sequential replacements by clones of higher fitness (periodic selection), bacterial populations would be continually purged of genetic variability, and the fate of selectively neutral alleles in very large populations of bacteria would be similar to that in demes of sexually reproducing organisms with small genetically effective population sizes. The significance of periodic selection in reducing genetic variability in these clonally reproducing species is dependent on the amount of genetic exchange between clones (recombination). In an effort to determine the relationship between the rates of periodic selection, recombination and the genetically effective sizes of bacterial populations, a model for periodic selection and infectious gene exchange has been developed and its properties analyzed. It shows that, for a given periodic selection regime, genetically effective population size increases exponentially with the rate of recombination.—With the parameters of this model in the range anticipated for natural populations of E. coli, the purging effects of periodic selection on genetic variability are significant; individual populations or lineages of this bacterial species would have very small genetically effective population sizes.—Based on this result, some other a priori considerations and a review of the results of epidemiological and genetic variability studies, it is postulated that E. coli is composed of a relatively limited number of geographically widespread and genetically nearly isolated and monomorphic lineages. The implications of these considerations of the genetic structure of E. coli populations of the interpretation of protein variation and the neutral gene hypothesis are discussed.     

The relationship of nucleotide polymorphism, recombination rate and selection in wild tomato species

We analyzed the effects of mating system and recombination rate on single nucleotide polymorphisms using 14 single-copy nuclear loci from single populations of five species of wild tomatoes (Solanum section Lycopersicon). The taxa investigated comprise two self-compatible (SC) and three self-incompatible (SI) species. The observed reduction in nucleotide diversity in the SC populations compared to the SI populations is much stronger than expected under the neutral effects of the mating system on effective population size. Importantly, outgroup sequences available for 11 of the 14 loci yield strong positive correlations between silent nucleotide diversity and silent divergence, indicative of marked among-locus differences in mutation rates and/or selective constraints. Furthermore, using a physical estimate of local recombination rates, we find that silent nucleotide diversity (but not divergence) is positively correlated with recombination rate in two of the SI species. However, this correlation is not nearly as strong as in other well-characterized species (in particular, Drosophila). We propose that nucleotide diversity in Lycopersicon is dominated mainly by differences in neutral mutation rates and/or selective constraints among loci, demographic processes (such as population subdivision), and background selection. In addition, we hypothesize that the soil seed bank plays an important role in the maintenance of the large genetic diversity in the SI species (in particular L. peruvianum).     

Microsatellite variability in wild populations of the house mouse is not influenced by differences in chromosomal recombination rates

Differences in recombination rates along the chromosomes can influence the evolution of neutral loci via hitchhiking effects. Generally, these effects should be stronger in regions of low recombination than in regions of high recombination. Detailed information on physical and genetic maps in the house mouse now allows an assessment of the correlation between neutral variability and recombination rates at given chromosomal regions. We chose 29 microsatellite loci from chromosomal regions which show differences in recombination rates and tested their variability in samples from five wild populations of Mus musculus musculus and M. m. domesticus . Our results provide no evidence for a correlation between microsatellite variability and recombination rates. This suggests that the high average mutation rate of microsatellites in mammals counterbalances the effects of long-range hitchhiking in the mouse genome.  © 2005 The Linnean Society of London, Biological Journal of the Linnean Society , 2005, 84 , 629–635.     

An approximate bayesian estimator suggests strong, recurrent selective sweeps in Drosophila

The recurrent fixation of newly arising, beneficial mutations in a species reduces levels of linked neutral variability. Models positing frequent weakly beneficial substitutions or, alternatively, rare, strongly selected substitutions predict similar average effects on linked neutral variability, if the product of the rate and strength of selection is held constant. We propose an approximate Bayesian (ABC) polymorphism-based estimator that can be used to distinguish between these models, and apply it to multi-locus data from Drosophila melanogaster. We investigate the extent to which inference about the strength of selection is sensitive to assumptions about the underlying distributions of the rates of substitution and recombination, the strength of selection, heterogeneity in mutation rate, as well as the population''s demographic history. We show that assuming fixed values of selection parameters in estimation leads to overestimates of the strength of selection and underestimates of the rate. We estimate parameters for an African population of D. melanogaster (ŝ∼2E−03, ) and compare these to previous estimates. Finally, we show that surveying larger genomic regions is expected to lend much more discriminatory power to the approach. It will thus be of great interest to apply this method to emerging whole-genome polymorphism data sets in many taxa.