How closely does genetic diversity in finite populations conform to predictions of neutral theory? Large deficits in regions of low recombination

Levels of genetic diversity in finite populations are crucial in conservation and evolutionary biology. Genetic diversity is required for populations to evolve and its loss is related to inbreeding in random mating populations, and thus to reduced population fitness and increased extinction risk. Neutral theory is widely used to predict levels of genetic diversity. I review levels of genetic diversity in finite populations in relation to predictions of neutral theory. Positive associations between genetic diversity and population size, as predicted by neutral theory, are observed for microsatellites, allozymes, quantitative genetic variation and usually for mitochondrial DNA (mtDNA). However, there are frequently significant deviations from neutral theory owing to indirect selection at linked loci caused by balancing selection, selective sweeps and background selection. Substantially lower genetic diversity than predicted under neutrality was found for chromosomes with low recombination rates and high linkage disequilibrium (compared with 'normally' recombining chromosomes within species and adjusted for different copy numbers and mutation rates), including W (median 100% lower) and Y (89% lower) chromosomes, dot fourth chromosomes in Drosophila (94% lower) and mtDNA (67% lower). Further, microsatellite genetic and allelic diversity were lost at 12 and 33% faster rates than expected in populations adapting to captivity, owing to widespread selective sweeps. Overall, neither neutral theory nor most versions of the genetic draft hypothesis are compatible with all empirical results.     

The distribution of copia-type retrotransposons and the evolutionary history of tomato and related wild species

Retrotransposons are mobile genetic elements that amplify throughout the genome and may be important contributors of genetic diversity. Their distribution is influenced by element behaviour and host-driven controls. We analysed the distribution of three copia-type retrotransposons, ToRTL1, T135 and Tnt1 using sequence-specific amplification polymorphism in self-compatible (SC) and incompatible (SI) species of Solanum subsection Lycopersicon, and genetically mapped polymorphic insertions in S. lycopersicum (tomato). The majority of polymorphic insertions (61%) are located in centromeric regions of the tomato genome. A significant positive relationship was detected between insertion polymorphisms and mating system, independent of selection as most insertions were found to be neutral. As insertion patterns successfully inferred interspecific relationships of Solanum subsection Lycopersicon, our results suggest that the distribution of ToRTL1, T135 and Tnt1 may essentially be determined by selection removing strongly deleterious insertions, with genetic drift and mating system, but not recombination rate, playing important roles.     

Testing for effects of recombination rate on nucleotide diversity in natural populations of Arabidopsis lyrata

We investigated DNA sequence diversity for loci on chromosomes 1 and 2 in six natural populations of Arabidopsis lyrata and tested for the role of natural selection in structuring genomewide patterns of variability, specifically examining the effects of recombination rate on levels of silent polymorphism. In contrast with theoretical predictions from models of genetic hitchhiking, maximum-likelihood-based analyses of diversity and divergence do not suggest reduction of diversity in the region of suppressed recombination near the centromere of chromosome 1, except in a single population from Russia, in which the pericentromeric region may have undergone a local selective sweep or demographic process that reduced variability. We discuss various possibilities that might explain why nucleotide diversity in most A. lyrata populations is not related to recombination rate, including genic recombination hotspots, and low gene density in the low recombination rate region.     

Breakdown of self-incompatibility in the perennial Arabidopsis lyrata (Brassicaceae) and its genetic consequences

Mating systems in plants are known to be highly labile traits, with frequent transitions from outcrossing to selfing. The genetic basis for breakdown in self-incompatibility (SI) systems has been studied, but data on variation in selfing rates in species for which the molecular basis of SI is known are rare. This study surveyed such variation in Arabidopsis lyrata (Brassicaceae), which is often considered an obligately outcrossing species, to examine the causes and genetic consequences of changes in its breeding system. Based on controlled self-pollinations in the greenhouse, three populations from the Great Lakes region of North America included a minority of self-compatible (SC) individuals, while two showed larger proportions of SC individuals and all populations contained some individuals capable of setting selfed seeds. Loss of SI was not associated with particular haplotypes at the S-locus (as estimated by alleles amplified at the SRK locus, the gene controlling female specificity) and all populations contained similar numbers of SRK alleles, suggesting that some other genetic factor is responsible for modifying the SI reaction. The loss of SI has resulted in an effective shift in the mating system, as the two populations with a high frequency of SC individuals showed significantly lower microsatellite-based multilocus outcrossing rates and higher inbreeding coefficients than the other populations. Based on microsatellites, observed heterozygosities and genetic diversity were also significantly depressed in these populations. These findings provide the unique opportunity to examine in detail the consequences of mating system changes within a species with a well-characterized SI system.     

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.     

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.     

Patterns of nucleotide diversity in two species of Mimulus are affected by mating system and asymmetric introgression

The evolutionary transition from outcrossing to self-fertilization has far-reaching implications for patterns of intraspecific genetic diversity and the potential for speciation. Using DNA sequence variation at two nuclear loci, we examined the divergence history of two closely related species of Mimulus. To investigate the effects of mating system and introgressive hybridization on the outcrossing M. guttatus and the selfing M. nasutus, we inspected nucleotide diversity within and between natural populations spanning the species' geographic ranges. High sequence similarity among populations of the selfing M. nasutus points to a single evolutionary origin for the species. Consistent with their distinct mating systems, all genetic variation in M. nasutus is distributed among populations, whereas M. guttatus exhibits appreciable levels of nucleotide diversity within populations. Silent genetic diversity is extensive in M. guttatus (mean theta(sil)/site = 0.077) and greatly exceeds the predicted twofold elevation in neutral variation for outcrossers relative to selfers. The finding of several M. guttatus sequences that share complete identity with sequences from M. nasutus suggests that recent asymmetric introgression may have occurred. We argue that exceptionally high nucleotide diversity in M. guttatus is consistent with a long-term history of directional introgression from M. nasutus to M. guttatus throughout the divergence of these two species.     

Pervasive recombination and sympatric genome diversification driven by frequency-dependent selection in Borrelia burgdorferi, the Lyme disease bacterium

How genomic diversity within bacterial populations originates and is maintained in the presence of frequent recombination is a central problem in understanding bacterial evolution. Natural populations of Borrelia burgdorferi, the bacterial agent of Lyme disease, consist of diverse genomic groups co-infecting single individual vertebrate hosts and tick vectors. To understand mechanisms of sympatric genome differentiation in B. burgdorferi, we sequenced and compared 23 genomes representing major genomic groups in North America and Europe. Linkage analysis of >13,500 single-nucleotide polymorphisms revealed pervasive horizontal DNA exchanges. Although three times more frequent than point mutation, recombination is localized and weakly affects genome-wide linkage disequilibrium. We show by computer simulations that, while enhancing population fitness, recombination constrains neutral and adaptive divergence among sympatric genomes through periodic selective sweeps. In contrast, simulations of frequency-dependent selection with recombination produced the observed pattern of a large number of sympatric genomic groups associated with major sequence variations at the selected locus. We conclude that negative frequency-dependent selection targeting a small number of surface-antigen loci (ospC in particular) sufficiently explains the maintenance of sympatric genome diversity in B. burgdorferi without adaptive divergence. We suggest that pervasive recombination makes it less likely for local B. burgdorferi genomic groups to achieve host specialization. B. burgdorferi genomic groups in the northeastern United States are thus best viewed as constituting a single bacterial species, whose generalist nature is a key to its rapid spread and human virulence.     

Using multilocus sequence data to assess population structure, natural selection, and linkage disequilibrium in wild tomatoes

We employed a multilocus approach to examine the effects of population subdivision and natural selection on DNA polymorphism in 2 closely related wild tomato species (Solanum peruvianum and Solanum chilense), using sequence data for 8 nuclear loci from populations across much of the species' range. Both species exhibit substantial levels of nucleotide variation. The species-wide level of silent nucleotide diversity is 18% higher in S. peruvianum (pi(sil) approximately 2.50%) than in S. chilense (pi(sil) approximately 2.12%). One of the loci deviates from neutral expectations, showing a clinal pattern of nucleotide diversity and haplotype structure in S. chilense. This geographic pattern of variation is suggestive of an incomplete (ongoing) selective sweep, but neutral explanations cannot be entirely dismissed. Both wild tomato species exhibit moderate levels of population differentiation (average F(ST) approximately 0.20). Interestingly, the pooled samples (across different demes) exhibit more negative Tajima's D and Fu and Li's D values; this marked excess of low-frequency polymorphism can only be explained by population (or range) expansion and is unlikely to be due to population structure per se. We thus propose that population structure and population/range expansion are among the most important evolutionary forces shaping patterns of nucleotide diversity within and among demes in these wild tomatoes. Patterns of population differentiation may also be impacted by soil seed banks and historical associations mediated by climatic cycles. Intragenic linkage disequilibrium (LD) decays very rapidly with physical distance, suggesting high recombination rates and effective population sizes in both species. The rapid decline of LD seems very promising for future association studies with the purpose of mapping functional variation in wild tomatoes.     

Species and recombination effects on DNA variability in the tomato genus

Population genetics theory predicts that strong selection for rare, beneficial mutations or against frequent, deleterious mutations reduces polymorphism at linked neutral (or weakly selected) sites. The reduction of genetic variation is expected to be more severe when recombination rates are lower. In outbreeding species, low recombination rates are usually confined to certain chromosomal regions, such as centromeres and telomeres. In contrast, in predominantly selfing species, the rarity of double heterozygotes leads to a reduced effective recombination rate in the whole genome. We investigated the effects of restricted recombination on DNA polymorphism in these two cases, analyzing five Lycopersicon species with contrasting mating systems: L. chilense, L. hirsutum, L. peruvianum, L. chmielewskii, and L. pimpinellifolium, of which only the first three species have self-incompatibility alleles. In each species, we determined DNA sequence variation of five single-copy genes located in chromosomal regions with either high or low recombination rate. We found that the mating system has a highly significant effect on the level of polymorphism, whereas recombination has only a weak influence. The effect of recombination on levels of polymorphism in Lycopersicon is much weaker than in other well-studied species, including Drosophila. To explain these observations, we discuss a number of hypotheses, invoking selection, recombination, and demographic factors associated with the mating system. We also provide evidence that L. peruvianum, showing a level of polymorphism (almost 3%) that is comparable to the level of divergence in the whole genus, is the ancestral species from which the other species of the genus Lycopersicon have originated relatively recently.     

Patterns of DNA Variability at X-Linked Loci in Mus Domesticus

Introns of four X-linked genes (Hprt, Plp, Glra2, and Amg) were sequenced to provide an estimate of nucleotide diversity at nuclear genes within the house mouse and to test the neutral prediction that the ratio of intraspecific polymorphism to interspecific divergence is the same for different loci. Hprt and Plp lie in a region of the X chromosome that experiences relatively low recombination rates, while Glra2 and Amg lie near the telomere of the X chromosome, a region that experiences higher recombination rates. A total of 6022 bases were sequenced in each of 10 Mus domesticus and one M. caroli. Average nucleotide diversity (π) for introns within M. domesticus was quite low (π = 0.078%). However, there was substantial variation in the level of heterozygosity among loci. The two telomeric loci, Glra2 and Amg, had higher ratios of polymorphism to divergence than the two loci experiencing lower recombination rates. These results are consistent with the hypothesis that heterozygosity is reduced in regions with lower rates of recombination, although sampling of additional genes is needed to establish whether there is a general correlation between heterozygosity and recombination rate as in Drosophila melanogaster.     

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.     

DNA variation in a conifer,Cryptomeria japonica (Cupressaceae sensu lato)

We investigated the nucleotide variation of a conifer, Cryptomeria japonica, and the divergence between this species and its closest relative, Taxodium distichum, at seven nuclear loci (Acl5, Chi1, Ferr, GapC, HemA, Lcyb, and Pat). Samples of C. japonica were collected from three areas, Kantou-Toukai, Hokuriku, and Iwate. No apparent geographic differentiation was found among these samples. However, the frequency spectrum of the nucleotide polymorphism revealed excesses of intermediate-frequency variants, which suggests that the population was not panmictic and a constant size in the past. The average nucleotide diversity, pi, for silent sites was 0.00383. However, values of pi for silent sites vary among loci. Comparisons of polymorphism to divergence among loci (the HKA test) showed that the polymorphism at the Acl5 locus was significantly lower. We also observed a nearly significant excess of replacement polymorphisms at the Lcyb locus. These results suggested possibilities of natural selection acting at some of the loci. Intragenic recombination was detected only once at the Chi1 locus and was not detected at the other loci. The low level of population recombination rate, 4Nr, seemed to be due to both low level of recombination, r, and small population size, N.     

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.     

Genetic variation and population structure of a threatened timber tree <Emphasis Type="Italic">Dalbergia cochinchinensis</Emphasis> in Cambodia

Dalbergia cochinchinensis Pierre ex Laness. (Fabaceae) is a commercially important tree in Southeast Asia. Although this species is under legal protections, illegal logging and disorderly developments have reduced its populations, and the conservation of this species is currently of much concern. In this study, we determined nucleotide sequences at six chloroplasts and ten nuclear loci in four populations of D. cochinchinensis in Cambodia, followed by population genetic analyses. The average silent nucleotide diversity over the nuclear loci, excluding one with an exceptionally high value, was 0.0057 in the entire population, and the mean F _ST across the nuclear loci between each population pair was between 0.135 and 0.467. Thus, the nucleotide diversity in the studied populations was not low compared with that in other tree species, and the level of population differentiation was high. Neutrality test statistics indicated a recent reduction of population size and a subdivision of the population within this species. The divergence times and migration rates were estimated with a likelihood-based method assuming the isolation with migration model. Based on the results, the three populations split 68,000–138,000 years ago, possibly corresponding to the start of the last glacial period, and the level of gene flow among the populations was very low thereafter. Moreover, after the split, population sizes were reduced considerably. Notably, the nucleotide diversity in an insertion sequence in a noncoding region of nuclear C4H was much higher than the mean nucleotide diversity in silent sites across other nuclear genes, indicating that the region was affected by selection.     

Comparison of DNA and protein polymorphisms between humans and chimpanzees

To examine the nucleotide diversity at silent (synonymous + intron + untranslated) and non-silent (nonsynonymous) sites in chimpanzees and humans, genes at six nuclear loci from two chimpanzees were sequenced. The average silent diversity was 0.19%, which was significantly higher than that in humans (0.05%). This observation suggests a significantly larger effective population size and a higher extent of neutral polymorphism in chimpanzees than in humans. On the other hand, the non-silent nucleotide diversity is similar in both species, resulting in a larger fraction of neutral mutations at non-silent sites in humans than in chimpanzees. Other types of polymorphism data were collected from the literature or databases to examine whether or not they are consistent with the nuclear DNA sequence polymorphism observed here. The nucleotide diversity at both silent and non-silent sites in mitochondrial (mt) DNA genes was compatible with that of the nuclear genes. Microsatellite loci showed a similar high extent of heterozygosity in both species, perhaps due to the combined effect of a high mutation rate and a recent population expansion in humans. At protein loci, humans are more heterozygous than chimpanzees, and the estimated fraction of neutral alleles in humans (0.84) is much larger than that in chimpanzees (0.26). These data show that the neutral fraction in non-silent changes is relatively large in the human population. This difference may be due to a relaxation of the functional constraint against proteins in the human lineage. To evaluate this possibility, it will be necessary to examine nucleotide sequences in relation to the physiological or biochemical properties of proteins.     

Genetic diversity and structure in two species of Leavenworthia with self-incompatible and self-compatible populations

Self-fertilization is a common mating system in plants and is known to reduce genetic diversity, increase genetic structure and potentially put populations at greater risk of extinction. In this study, we measured the genetic diversity and structure of two cedar glade endemic species, Leavenworthia alabamica and L. crassa. These species have self-incompatible (SI) and self-compatible (SC) populations and are therefore ideal for understanding how the mating system affects genetic diversity and structure. We found that L. alabamica and L. crassa had high species-level genetic diversity (H_e=0.229 and 0.183, respectively) and high genetic structure among their populations (F_ST=0.45 and 0.36, respectively), but that mean genetic diversity was significantly lower in SC compared with SI populations (SC vs SI, H_e for L. alabamica was 0.065 vs 0.206 and for L. crassa was 0.084 vs 0.189). We also found significant genetic structure using maximum-likelihood clustering methods. These data indicate that the loss of SI leads to the loss of genetic diversity within populations. In addition, we examined genetic distance relationships between SI and SC populations to analyze possible population history and origins of self-compatibility. We find there may have been multiple origins of self-compatibility in L. alabamica and L. crassa. However, further work is required to test this hypothesis. Finally, given their high genetic structure and that individual populations harbor unique alleles, conservation strategies seeking to maximize species-level genetic diversity for these or similar species should protect multiple populations.     

Demographic history has influenced nucleotide diversity in European Pinus sylvestris populations

To infer the role of natural selection in shaping standing genetic diversity, it is necessary to assess the genomewide impact of demographic history on nucleotide diversity. In this study we analyzed sequence diversity of 16 nuclear loci in eight Pinus sylvestris populations. Populations were divided into four geographical groups on the basis of their current location and the geographical history of the region: northern Europe, central Europe, Spain, and Turkey. There were no among-group differences in the level of silent nucleotide diversity, which was approximately 0.005/bp in all groups. There was some evidence that linkage disequilibrium extended further in northern Europe than in central Europe: the estimates of the population recombination rate parameter, rho, were 0.0064 and 0.0294, respectively. The summary statistics of nucleotide diversity in central and northern European populations were compatible with an ancient bottleneck rather than the standard neutral model.