Recombination rate variation and speciation: theoretical predictions and empirical results from rabbits and mice

Recently diverged taxa may continue to exchange genes. A number of models of speciation with gene flow propose that the frequency of gene exchange will be lower in genomic regions of low recombination and that these regions will therefore be more differentiated. However, several population-genetic models that focus on selection at linked sites also predict greater differentiation in regions of low recombination simply as a result of faster sorting of ancestral alleles even in the absence of gene flow. Moreover, identifying the actual amount of gene flow from patterns of genetic variation is tricky, because both ancestral polymorphism and migration lead to shared variation between recently diverged taxa. New analytic methods have been developed to help distinguish ancestral polymorphism from migration. Along with a growing number of datasets of multi-locus DNA sequence variation, these methods have spawned a renewed interest in speciation models with gene flow. Here, we review both speciation and population-genetic models that make explicit predictions about how the rate of recombination influences patterns of genetic variation within and between species. We then compare those predictions with empirical data of DNA sequence variation in rabbits and mice. We find strong support for the prediction that genomic regions experiencing low levels of recombination are more differentiated. In most cases, reduced gene flow appears to contribute to the pattern, although disentangling the relative contribution of reduced gene flow and selection at linked sites remains a challenge. We suggest fruitful areas of research that might help distinguish between different models.     

The Effects of Background and Interference Selection on Patterns of Genetic Variation in Subdivided Populations

It is well known that most new mutations that affect fitness exert deleterious effects and that natural populations are often composed of subpopulations (demes) connected by gene flow. To gain a better understanding of the joint effects of purifying selection and population structure, we focus on a scenario where an ancestral population splits into multiple demes and study neutral diversity patterns in regions linked to selected sites. In the background selection regime of strong selection, we first derive analytic equations for pairwise coalescent times and F_ST as a function of time after the ancestral population splits into two demes and then construct a flexible coalescent simulator that can generate samples under complex models such as those involving multiple demes or nonconservative migration. We have carried out extensive forward simulations to show that the new methods can accurately predict diversity patterns both in the nonequilibrium phase following the split of the ancestral population and in the equilibrium between mutation, migration, drift, and selection. In the interference selection regime of many tightly linked selected sites, forward simulations provide evidence that neutral diversity patterns obtained from both the nonequilibrium and equilibrium phases may be virtually indistinguishable for models that have identical variance in fitness, but are nonetheless different with respect to the number of selected sites and the strength of purifying selection. This equivalence in neutral diversity patterns suggests that data collected from subdivided populations may have limited power for differentiating among the selective pressures to which closely linked selected sites are subject.     

Genomewide patterns of variation in genetic diversity are shared among populations,species and higher‐order taxa

Genomewide screens of genetic variation within and between populations can reveal signatures of selection implicated in adaptation and speciation. Genomic regions with low genetic diversity and elevated differentiation reflective of locally reduced effective population sizes (N_e) are candidates for barrier loci contributing to population divergence. Yet, such candidate genomic regions need not arise as a result of selection promoting adaptation or advancing reproductive isolation. Linked selection unrelated to lineage‐specific adaptation or population divergence can generate comparable signatures. It is challenging to distinguish between these processes, particularly when diverging populations share ancestral genetic variation. In this study, we took a comparative approach using population assemblages from distant clades assessing genomic parallelism of variation in N_e. Utilizing population‐level polymorphism data from 444 resequenced genomes of three avian clades spanning 50 million years of evolution, we tested whether population genetic summary statistics reflecting genomewide variation in N_e would covary among populations within clades, and importantly, also among clades where lineage sorting has been completed. All statistics including population‐scaled recombination rate (ρ), nucleotide diversity (π) and measures of genetic differentiation between populations (F_ST, PBS, d_xy) were significantly correlated across all phylogenetic distances. Moreover, genomic regions with elevated levels of genetic differentiation were associated with inferred pericentromeric and subtelomeric regions. The phylogenetic stability of diversity landscapes and stable association with genomic features support a role of linked selection not necessarily associated with adaptation and speciation in shaping patterns of genomewide heterogeneity in genetic diversity.     

Genetic variation in mouse salivary amylase rate of synthesis

Heterozygotes from matings of the mouse strains YBR/Cv and C3H/As have about 3 times more YBR-amylase than C3H-amylase in the saliva. The determinant for this quantitative effect is located on linkage group XVI close to or within the structural gene for salivary amylase. The quantitative effect is the result of an increase in the rate of synthesis of YBR-amylase, and the determinant is cis acting. Studies of other mouse strains suggest that regulatory genetic elements may modulate salivary amylase production.This work was supported by the Danish Natural Science Research Council and a grant from the United States Public Health Service (Grant GM-19521).     

Genetic variation and population subdivision of the endangered Iberian cyprinid Chondrostoma lusitanicum

Products of 24 presumptive enzyme loci were used to analyse the consequences on genetic structure and variation of the Chondrostoma lusitanicum population decline within the Tejo basin. This rare cyprinid is endemic to the Iberian Peninsula and has a very restricted distribution. Five samples from the Tejo catchment were compared with a sample from a small basin, the Samarra, that has not suffered obvious anthropogenic pressures and where the fish is still abundant. Heterozygosity and polymorphism were higher overall in the Samarra. Several lines of evidence indicate a high degree of population subdivision within the Tejo basin. In fact, about half of the total gene diversity detected in the Tejo population was due to differences among samples. This differentiation appeared to be caused by genetic drift and possibly differential local selection, coupled with reduced gene flow among localities. The accelerated process of habitat degradation occurring in the lowland streams of the Tejo basin will lead to the inevitable reduction of intraspecies genetic diversity.     

A coalescent model of background selection with recombination,demography and variation in selection coefficients

There is increasing evidence that background selection, the effects of the elimination of recurring deleterious mutations by natural selection on variability at linked sites, may be a major factor shaping genome-wide patterns of genetic diversity. To accurately quantify the importance of background selection, it is vital to have computationally efficient models that include essential biological features. To this end, a structured coalescent procedure is used to construct a model of background selection that takes into account the effects of recombination, recent changes in population size and variation in selection coefficients against deleterious mutations across sites. Furthermore, this model allows a flexible organization of selected and neutral sites in the region concerned, and has the ability to generate sequence variability at both selected and neutral sites, allowing the correlation between these two types of sites to be studied. The accuracy of the model is verified by checking against the results of forward simulations. These simulations also reveal several patterns of diversity that are in qualitative agreement with observations reported in recent studies of DNA sequence polymorphisms. These results suggest that the model should be useful for data analysis.     

Genomics of adaptive divergence with chromosome‐scale heterogeneity in crossover rate

Genetic differentiation between divergent populations is often greater in chromosome centres than peripheries. Commonly overlooked, this broadscale differentiation pattern is sometimes ascribed to heterogeneity in crossover rate and hence linked selection within chromosomes, but the underlying mechanisms remain incompletely understood. A literature survey across 46 organisms reveals that most eukaryotes indeed exhibit a reduced crossover rate in chromosome centres relative to the peripheries. Using simulations of populations diverging into ecologically different habitats through sorting of standing genetic variation, we demonstrate that such chromosome‐scale heterogeneity in crossover rate, combined with polygenic divergent selection, causes stronger hitchhiking and especially barriers to gene flow across chromosome centres. Without requiring selection on new mutations, this rapidly leads to elevated population differentiation in the low‐crossover centres relative to the high‐crossover peripheries of chromosomes (“Chromosome Centre‐Biased Differentiation”, CCBD). Using simulated and empirical data, we then show that strong CCBD between populations can provide evidence of polygenic adaptive divergence with a phase of gene flow. We further demonstrate that chromosome‐scale heterogeneity in crossover rate impacts analyses beyond that of population differentiation, including the inference of phylogenies and parallel adaptive evolution among populations, the detection of genetic loci under selection, and the interpretation of the strength of selection on genomic regions. Overall, our results call for a greater appreciation of chromosome‐scale heterogeneity in crossover rate in evolutionary genomics.     

Genetic population structure and call variation in a passerine bird, the satin bowerbird, Ptilonorhynchus violaceus

Geographic variation in vocalizations is widespread in passerine birds, but its origins and maintenance remain unclear. One hypothesis to explain this variation is that it is associated with geographic isolation among populations and therefore should follow a vicariant pattern similar to that typically found in neutral genetic markers. Alternatively, if environmental selection strongly influences vocalizations, then genetic divergence and vocal divergence may be disassociated. This study compared genetic divergence derived from 11 microsatellite markers with a metric of phenotypic divergence derived from male bower advertisement calls. Data were obtained from 16 populations throughout the entire distribution of the satin bowerbird, an Australian wet-forest-restricted passerine. There was no relationship between call divergence and genetic divergence, similar to most other studies on birds with learned vocalizations. Genetic divergence followed a vicariant model of evolution, with the differentiation of isolated populations and isolation-by-distance among continuous populations. Previous work on Ptilonorhynchus violaceus has shown that advertisement call structure is strongly influenced by the acoustic environment of different habitats. Divergence in vocalizations among genetically related populations in different habitats indicates that satin bowerbirds match their vocalizations to the environment in which they live, despite the homogenizing influence of gene flow. In combination with convergence of vocalizations among genetically divergent populations occurring in the same habitat, this shows the overriding importance that habitat-related selection can have on the establishment and maintenance of variation in vocalizations.     

The effects of mice on stoats in southern beech forests

Introduced stoats (Mustela erminea) are important invasive predators in southern beech (Nothofagus sp.) forests in New Zealand. In these forests, one of their primary prey species – introduced house mice (Mus musculus), fluctuate dramatically between years, driven by the irregular heavy seed‐fall (masting) of the beech trees. We examined the effects of mice on stoats in this system by comparing the weights, age structure and population densities of stoats caught on two large islands in Fiordland, New Zealand – one that has mice (Resolution Island) and one that does not (Secretary Island). On Resolution Island, the stoat population showed a history of recruitment spikes and troughs linked to beech masting, whereas the Secretary Island population had more constant recruitment, indicating that rodents are probably the primary cause for the ‘boom and bust’ population cycle of stoats in beech forests. Resolutions Island stoats were 10% heavier on average than Secretary Island stoats, supporting the hypothesis that the availability of larger prey (mice verses wētā) leads to larger stoats. Beech masting years on this island were also correlated with a higher weight for stoats born in the year of the masting event. The detailed demographic information on the stoat populations of these two islands supports previously suggested interactions among mice, stoats and beech masting. These interactions may have important consequences for the endemic species that interact with fluctuating populations of mice and stoats.     

Frequent recombination shapes the epidemic population structure of Planktothrix (Cyanoprokaryota) in Italian subalpine lakes

The planktonic genus Planktothrix, as other cyanobacteria, shows signals of both homologous and nonhomologous recombination. However, the frequency of recombination and its effect on Planktothrix population structuring is unknown. We isolated 290 Planktothrix strains from seven neighboring lakes in the subalpine Italian region and analyzed these using multilocus sequence typing. Four of six loci analyzed were polymorphic, resulting in 20 distinct multilocus genotypes. Association indices among alleles at different loci were suggestive of an “epidemic population structure,” resulting from an explosive (and temporary) dominance of one genotype against a panmictic background. ClonalFrame analyses supported this view by detecting: (i) three major clades affected by three distinct recombination events, (ii) a recombination rate about equal to the mutation rate, and (iii) the fact that recombination had an impact on introducing molecular diversity more than double the mutation rate. Furthermore, analysis of molecular variance over an annual cycle in three of seven lakes revealed that both local clonal expansion and recombination processes affected among‐lake diversity. Our observations suggest that recombination affects microevolution of Planktothrix and that an epidemic structure can emerge in populations of this genus.     

Genetic variation and population structure ofCarex breviculmis (Cyperaceae) in Korea

We determined the genetic diversity and population structures ofCarex breviculmis (Cyperaceae) populations in Korea, using genetic variations at 23 allozyme loci.C. breviculmis is a long-lived herbaceous species that is widely distributed in eastern Asia. A high level of genetic variation was found in 15 populations. Twelve enzymes revealed 23 loci, of which 11 were polymorphic (47.8%). Genetic diversity at the speciesand population levels were 0.174 and 0.146, respectively. Total genetic diversity (H_T = 0.363) and within-population genetic diversity (H_s = 0.346) were high, whereas the extent of the population divergence was relatively low (G_ST = 0.063). Deviation from random mating (Fis) within the 15 populations was 0.206. An indirect estimate of the number of migrants per generation(Nm = 3.69) indicated that gene flow was extensive among Korean populations of this species. Analysis of fixation indices revealed a substantial heterozygote deficiency in some populations and at some loci. Genetic identity between popu-lations was high, exceeding 0.956.     

Selection at linked sites in the partial selfer Caenorhabditis elegans

Natural selection can produce a correlation between local recombination rates and levels of neutral DNA polymorphism as a consequence of genetic hitchhiking and background selection. Theory suggests that selection at linked sites should affect patterns of neutral variation in partially selfing populations more dramatically than in outcrossing populations. However, empirical investigations of selection at linked sites have focused primarily on outcrossing species. To assess the potential role of selection as a determinant of neutral polymorphism in the context of partial self-fertilization, we conducted a multivariate analysis of single-nucleotide polymorphism (SNP) density throughout the genome of the nematode Caenorhabditis elegans. We based the analysis on a published SNP data set and partitioned the genome into windows to calculate SNP densities, recombination rates, and gene densities across all six chromosomes. Our analyses identify a strong, positive correlation between recombination rate and neutral polymorphism (as estimated by noncoding SNP density) across the genome of C. elegans. Furthermore, we find that levels of neutral polymorphism are lower in gene-dense regions than in gene-poor regions in some analyses. Analyses incorporating local estimates of divergence between C. elegans and C. briggsae indicate that a mutational explanation alone is unlikely to explain the observed patterns. Consequently, we interpret these findings as evidence that natural selection shapes genome-wide patterns of neutral polymorphism in C. elegans. Our study provides the first demonstration of such an effect in a partially selfing animal. Explicit models of genetic hitchhiking and background selection can each adequately describe the relationship between recombination rate and SNP density, but only when they incorporate selfing rate. Clarification of the relative roles of genetic hitchhiking and background selection in C. elegans awaits the development of specific theoretical predictions that account for partial self-fertilization and biased sex ratios.     

Low genetic differentiation in a sedentary bird: house sparrow population genetics in a contiguous landscape

The house sparrow Passer domesticus has been declining in abundance in many localities, including Finland. We studied the genetic diversity and differentiation of the house sparrow populations across Finland in the 1980s, at the onset of the species'' decline in abundance. We genotyped 472 adult males (the less dispersive sex) from 13 locations in Finland (covering a range of 400 × 800 km) and one in Sweden (Stockholm) for 13 polymorphic microsatellite markers. Our analysis of Finnish ringing records showed that natal dispersal distances are limited (90% <16 km), which confirmed earlier finding from other countries. The Finnish populations were panmictic, and genetically very homogeneous and the limited dispersal was sufficiently large to maintain their connectivity. However, all Finnish populations differed significantly from the Stockholm population, even though direct geographical distance to it was often smaller than among Finnish populations. Hence, the open sea between Finland and Sweden appears to form a dispersal barrier for this species, whereas dispersal is much less constrained across the Finnish mainland (which lacks geographical barriers). Our findings provide a benchmark for conservation biologists and emphasize the influence of landscape structure on gene flow.     

Genetic variation and population structure in the endangered greater horseshoe bat Rhinolophus ferrumequinum

Following a dramatic decline last century, the British population of the endangered greater horseshoe bat Rhinolophus ferrumequinum is highly fragmented. To examine the consequences of fragmentation and limited dispersal on patterns of genetic structure and variation, we used microsatellite markers to screen bats from around 50% of the known maternity colonies in Britain, and two areas from continental Europe. Analyses revealed that Welsh and English colonies were genetically isolated. This, and lower variability in Britain than north France, may result from either genetic drift, or the species' colonization history. Gene flow among most neighbouring colonies was not generally restricted, with one exception. These findings have important implications for the ongoing conservation management of this species.     

Genetic variation and population structure of the mosquito Anopheles jeyporiensis in southern China

Genetic differentiation among populations of Anopheles jeyporiensis was examined using 76 mtDNA COII sequences from 16 sites throughout southern China and northern Vietnam. The COII sequences are AT-rich (74.58%) and reveal high levels of diversity with 39 of 685 sites polymorphic and 50 different haplotypes present. Genetic variation is high within populations and significant geographical structure was detected at both population and regional levels. In the larger samples, the distributions of haplotypes suggest recent population expansion.     

The effect of population bottlenecks on mutation rate evolution in asexual populations

In the absence of recombination, a mutator allele can spread through a population by hitchhiking with beneficial mutations that appear in its genetic background. Theoretical studies over the past decade have shown that the survival and fixation probability of beneficial mutations can be severely reduced by population size bottlenecks. Here, we use computational modelling and evolution experiments with the yeast S. cerevisiae to examine whether population bottlenecks can affect mutator dynamics in adapting asexual populations. In simulation, we show that population bottlenecks can inhibit mutator hitchhiking with beneficial mutations and are most effective at lower beneficial mutation supply rates. We then subjected experimental populations of yeast propagated at the same effective population size to three different bottleneck regimes and observed that the speed of mutator hitchhiking was significantly slower at smaller bottlenecks, consistent with our theoretical expectations. Our results, thus, suggest that bottlenecks can be an important factor in mutation rate evolution and can in certain circumstances act to stabilize or, at least, delay the progressive elevation of mutation rates in asexual populations. Additionally, our findings provide the first experimental support for the theoretically postulated effect of population bottlenecks on beneficial mutations and demonstrate the usefulness of studying mutator frequency dynamics for understanding the underlying dynamics of fitness‐affecting mutations.     

Relative performance of Bayesian clustering software for inferring population substructure and individual assignment at low levels of population differentiation

Traditional methods for characterizing genetic differentiation among populations rely on a priori grouping of individuals. Bayesian clustering methods avoid this limitation by using linkage and Hardy–Weinberg disequilibrium to decompose a sample of individuals into genetically distinct groups. There are several software programs available for Bayesian clustering analyses, all of which describe a decrease in the ability to detect distinct clusters as levels of genetic differentiation among populations decrease. However, no study has yet compared the performance of such methods at low levels of population differentiation, which may be common in species where populations have experienced recent separation or high levels of gene flow. We used simulated data to evaluate the performance of three Bayesian clustering software programs, PARTITION, STRUCTURE, and BAPS, at levels of population differentiation below F _ST=0.1. PARTITION was unable to correctly identify the number of subpopulations until levels of F _ST reached around 0.09. Both STRUCTURE and BAPS performed very well at low levels of population differentiation, and were able to correctly identify the number of subpopulations at F _ST around 0.03. The average proportion of an individual’s genome assigned to its true population of origin increased with increasing F _ST for both programs, reaching over 92% at an F _ST of 0.05. The average number of misassignments (assignments to the incorrect subpopulation) continued to decrease as F _ST increased, and when F _ST was 0.05, fewer than 3% of individuals were misassigned using either program. Both STRUCTURE and BAPS worked extremely well for inferring the number of clusters when clusters were not well-differentiated (F _ST=0.02–0.03), but our results suggest that F _ST must be at least 0.05 to reach an assignment accuracy of greater than 97%.     

Long-term experimental evolution in Escherichia coli. V. Effects of recombination with immigrant genotypes on the rate of bacterial evolution

This study builds upon an earlier experiment that examined the dynamics of mean fitness in evolving populations of Escherichia coli in which mutations were the sole source of genetic variation. During thousands of generations in a constant environment, the rate of improvement in mean fitness of these asexual populations slowed considerably from an initially rapid pace. In this study, we sought to determine whether sexual recombination with novel genotypes would reaccelerate the rate of adaption in these populations. To that end, treatment populations were propagated for an additional 1000 generations in the same environment as their ancestors, but they were periodically allowed to mate with an immigrant pool of genetically distinct Hfr (high frequency recombination) donors. These donors could transfer genes to the resident populations by conjugation, but the donors themselves could not grow in the experimental environment. Control populations were propagated under identical conditions, but in the absence of sexual recombination with the donors. All twelve control populations retained the ancestral alleles at every locus that was scored. In contrast, the sexual recombination treatment yielded dramatic increases in genetic variation. Thus, there was a profound effect of recombination on the rate of genetic change. However, the increased genetic variation in the treatment populations had no significant effect on the rate of adaptive evolution, as measured by changes in mean fitness relative to a common competitor. We then considered three hypotheses that might reconcile these two outcomes: recombination pressure, hitchhiking of recombinant genotypes in association with beneficial mutations, and complex selection dynamics whereby certain genotypes may have a selective advantage only within a particular milieu of competitors. The estimated recombination rate was too low to explain the observed rate of genetic change, either alone or in combination with hitchhiking effects. However, we documented comple x ecological interactions among some recombinant genotypes, suggesting that our method for estimating fitness relative to a common competitor might have underestimated the rate of adaptive evolution in the treatment populations.     

Disentangling the effects of recombination,selection, and demography on the genetic variation at a major histocompatibility complex class II gene in the alpine chamois

The major histocompatibility complex (MHC) harbours some of the most polymorphic loci in vertebrate genomes. MHC genes are thought to be subject to some form of balancing selection, most likely pathogen‐mediated selection. Hence, MHC genes are excellent candidates for exploring adaptive processes. In this study, we investigated the genetic variation at exon 2 of the DRB class II MHC locus in 191 alpine chamois (Rupicapra rupicapra) from 10 populations in the eastern Alps of Italy. In particular, we were interested in distinguishing and estimating the relative impact of selective and demographic factors, while taking into account the confounding effect of recombination. The extremely high d_n/d_s ratio and the presence of trans‐species polymorphisms suggest that a strong long‐term balancing selection effect has been operating at this locus throughout the evolutionary history of this species. We analysed patterns of genetic variation within and between populations, and the mitochondrial D‐loop polymorphism patterns were analysed to provide a baseline indicator of the effects of demographic processes. These analyses showed that (i) the chamois experienced a demographic decline in the last 5000–30 000 years, most likely related to the postglacial elevation in temperature; (ii) this demographic process can explain the results of neutrality tests applied to MHC variation within populations, but cannot justify the much weaker divergence between populations implied by MHC as opposed to mitochondrial DNA; (iii) similar sets of divergent alleles are probably maintained with similar frequencies by balancing selection in different populations, and this mechanism is also operating in small isolated populations, which are strongly affected by drift.     

Genetic population structure, fitness variation and the importance of population history in remnant populations of the endangered plant Silene chlorantha (Willd.) Ehrh. (Caryophyllaceae)

Habitat fragmentation can lead to a decline of genetic diversity, a potential risk for the survival of natural populations. Fragmented populations can become highly differentiated due to reduced gene flow and genetic drift. A decline in number of individuals can result in lower reproductive fitness due to inbreeding effects. We investigated genetic variation within and between 11 populations of the rare and endangered plant Silene chlorantha in northeastern Germany to support conservation strategies. Genetic diversity was evaluated using AFLP techniques and the results were correlated to fitness traits. Fitness evaluation in nature and in a common garden approach was conducted. Our analysis revealed population differentiation was high and within population genetic diversity was intermediate. A clear population structure was supported by a Bayesian approach, AMOVA and neighbour-joining analysis. No correlation between genetic and geographic distance was found. Our results indicate that patterns of population differentiation were mainly caused by temporal and/or spatial isolation and genetic drift. The fitness evaluation revealed that pollinator limitation and habitat quality seem, at present, to be more important to reproductive fitness than genetic diversity by itself. Populations of S. chlorantha with low genetic diversity have the potential to increase in individual number if habitat conditions improve. This was detected in a single large population in the investigation area, which was formerly affected by bottleneck effects.