Evolution of the genomic rate of recombination in mammals

Rates of recombination vary considerably between species. Despite the significance of this observation for evolutionary biology and genetics, the evolutionary mechanisms that contribute to these interspecific differences are unclear. On fine physical scales, recombination rates appear to evolve rapidly between closely related species, but the mode and tempo of recombination rate evolution on the broader scale is poorly understood. Here, we use phylogenetic comparative methods to begin to characterize the evolutionary processes underlying average genomic recombination rates in mammals. We document a strong phylogenetic effect in recombination rates, indicating that more closely related species tend to have more similar average rates of recombination. We demonstrate that this phylogenetic signal is not an artifact of errors in recombination rate estimation and show that it is robust to uncertainty in the mammalian phylogeny. Neutral evolutionary models present good fits to the data and we find no evidence for heterogeneity in the rate of evolution in recombination across the mammalian tree. These results suggest that observed interspecific variation in average genomic rates of recombination is largely attributable to the steady accumulation of neutral mutations over evolutionary time. Although single recombination hotspots may live and die on short evolutionary time scales, the strong phylogenetic signal in genomic recombination rates indicates that the pace of evolution on this scale may be considerably slower.     

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.     

Higher differentiation among subspecies of the house mouse (Mus musculus) in genomic regions with low recombination

In the early stages of reproductive isolation, genomic regions of reduced recombination are expected to show greater levels of differentiation, either because gene flow between species is reduced in these regions or because the effects of selection at linked sites within species are enhanced in these regions. Here, we study the patterns of DNA sequence variation at 27 autosomal loci among populations of Mus musculus musculus, M. m. domesticus, and M. m. castaneus, three subspecies of house mice with collinear genomes. We found that some loci exhibit considerable shared variation among subspecies, while others exhibit fixed differences. We used an isolation-with-gene-flow model to estimate divergence times and effective population sizes (N(e) ) and to disentangle ancestral variation from gene flow. Estimates of divergence time indicate that all three subspecies diverged from one another within a very short period of time approximately 350,000 years ago. Overall, N(e) for each subspecies was associated with the degree of genetic differentiation: M. m. musculus had the smallest N(e) and the greatest proportion of monophyletic gene genealogies, while M. m. castaneus had the largest N(e) and the smallest proportion of monophyletic gene genealogies. M. m. domesticus and M. m. musculus were more differentiated from each other than either were from M. m. castaneus, consistent with greater reproductive isolation between M. m. domesticus and M. m. musculus. F(ST) was significantly greater at loci experiencing low recombination rates compared to loci experiencing high recombination rates in comparisons between M. m. castaneus and M. m. musculus or M. m. domesticus. These results provide evidence that genomic regions with less recombination show greater differentiation, even in the absence of chromosomal rearrangements.     

Predictors of genomic differentiation within a hybrid taxon

Hybridization is increasingly recognized as an important evolutionary force. Novel genetic methods now enable us to address how the genomes of parental species are combined in hybrid lineages. However, we still do not know the relative importance of admixed proportions, genome architecture and local selection in shaping hybrid genomes. Here, we take advantage of the genetically divergent island populations of Italian sparrow on Crete, Corsica and Sicily to investigate the predictors of genomic variation within a hybrid taxon. We test if differentiation is affected by recombination rate, selection, or variation in ancestry proportions. We find that the relationship between recombination rate and differentiation is less pronounced within hybrid lineages than between the parent species, as expected if purging of minor parent ancestry in low recombination regions reduces the variation available for differentiation. In addition, we find that differentiation between islands is correlated with differences in signatures of selection in two out of three comparisons. Signatures of selection within islands are correlated across all islands, suggesting that shared selection may mould genomic differentiation. The best predictor of strong differentiation within islands is the degree of differentiation from house sparrow, and hence loci with Spanish sparrow ancestry may vary more freely. Jointly, this suggests that constraints and selection interact in shaping the genomic landscape of differentiation in this hybrid species.     

Recombination rate variation in mice from an isolated island

Recombination rate is a heritable trait that varies among individuals. Despite the major impact of recombination rate on patterns of genetic diversity and the efficacy of selection, natural variation in this phenotype remains poorly characterized. We present a comparison of genetic maps, sampling 1212 meioses, from a unique population of wild house mice (Mus musculus domesticus) that recently colonized remote Gough Island. Crosses to a mainland reference strain (WSB/EiJ) reveal pervasive variation in recombination rate among Gough Island mice, including subchromosomal intervals spanning up to 28% of the genome. In spite of this high level of polymorphism, the genomewide recombination rate does not significantly vary. In general, we find that recombination rate varies more when measured in smaller genomic intervals. Using the current standard genetic map of the laboratory mouse to polarize intervals with divergent recombination rates, we infer that the majority of evolutionary change occurred in one of the two tested lines of Gough Island mice. Our results confirm that natural populations harbour a high level of recombination rate polymorphism and highlight the disparities in recombination rate evolution across genomic scales.     

Divergence between the Drosophila pseudoobscura and D. persimilis genome sequences in relation to chromosomal inversions

As whole-genome sequence assemblies accumulate, a challenge is to determine how these can be used to address fundamental evolutionary questions, such as inferring the process of speciation. Here, we use the sequence assemblies of Drosophila pseudoobscura and D. persimilis to test hypotheses regarding divergence with gene flow. We observe low differentiation between the two genome sequences in pericentromeric and peritelomeric regions. We interpret this result as primarily a remnant of the correlation between levels of variation and local recombination rate observed within populations. However, we also observe lower differentiation far from the fixed chromosomal inversions distinguishing these species and greater differentiation within and near these inversions. This finding is consistent with models suggesting that chromosomal inversions facilitate species divergence despite interspecies gene flow. We also document heterogeneity among the inverted regions in their degree of differentiation, suggesting temporal differences in the origin of each inverted region consistent with the inversions arising during a process of divergence with gene flow. While this study provides insights into the speciation process using two single-genome sequences, it was informed by lower throughput but more rigorous examinations of polymorphism and divergence. This reliance highlights the need for complementary genomic and population genetic approaches for tackling fundamental evolutionary questions such as speciation.     

Reduced nucleotide variability at an androgen-binding protein locus (Abpa) in house mice: evidence for positive natural selection.

Previous work has shown that the gene for the alpha subunit of androgen-binding protein, Abpa, may be involved in premating isolation between different subspecies of the house mouse, Mus musculus. We investigated patterns of DNA sequence variation at Abpa within and between species of mice to test several predictions of a model of neutral molecular evolution. Intraspecific variation among 10 Mus musculus domesticus alleles was compared with divergence between M. m. domesticus and M. caroli for Abpa and two X-linked genes, Glra2 and Amg. No variation was observed at Abpa within M. m. domesticus. The ratio of polymorphism to divergence was significantly lower at Abpa than at Glra2 and Amg, despite the fact that all three genes experience similar rates of recombination. Interspecific comparisons among M. m. domesticus, Mus musculus musculus, Mus musculus castaneus, Mus spretus, Mus spicilegus, and Mus caroli revealed that the ratio of nonsynonymous substitutions to synonymous substitutions on a per-site basis (Ka/Ks) was generally greater than one. The combined observations of no variation at Abpa within M. m: domesticus and uniformly high Ka/Ks values between species suggest that positive directional selection has acted recently at this locus.     

No Evidence that Infection Alters Global Recombination Rate in House Mice

Recombination rate is a complex trait, with genetic and environmental factors shaping observed patterns of variation. Although recent studies have begun to unravel the genetic basis of recombination rate differences between organisms, less attention has focused on the environmental determinants of crossover rates. Here, we test the effect of one ubiquitous environmental pressure–bacterial infection–on global recombination frequency in mammals. We applied MLH1 mapping to assay global crossover rates in male mice infected with the pathogenic bacterium Borrelia burgdorferi, the causative agent of Lyme Disease, and uninfected control animals. Despite ample statistical power to identify biologically relevant differences between infected and uninfected animals, we find no evidence for a global recombination rate response to bacterial infection. Moreover, broad-scale patterns of crossover distribution, including the number of achiasmate bivalents, are not affected by infection status. Although pathogen exposure can plastically increase recombination in some species, our findings suggest that recombination rates in house mice may be resilient to at least some forms of infection stress. This negative result motivates future experiments with alternative house mouse pathogens to evaluate the generality of this conclusion.     

Genetic analysis of genome-scale recombination rate evolution in house mice

The rate of meiotic recombination varies markedly between species and among individuals. Classical genetic experiments demonstrated a heritable component to population variation in recombination rate, and specific sequence variants that contribute to recombination rate differences between individuals have recently been identified. Despite these advances, the genetic basis of species divergence in recombination rate remains unexplored. Using a cytological assay that allows direct in situ imaging of recombination events in spermatocytes, we report a large (∼30%) difference in global recombination rate between males of two closely related house mouse subspecies (Mus musculus musculus and M. m. castaneus). To characterize the genetic basis of this recombination rate divergence, we generated an F2 panel of inter-subspecific hybrid males (n = 276) from an intercross between wild-derived inbred strains CAST/EiJ (M. m. castaneus) and PWD/PhJ (M. m. musculus). We uncover considerable heritable variation for recombination rate among males from this mapping population. Much of the F2 variance for recombination rate and a substantial portion of the difference in recombination rate between the parental strains is explained by eight moderate- to large-effect quantitative trait loci, including two transgressive loci on the X chromosome. In contrast to the rapid evolution observed in males, female CAST/EiJ and PWD/PhJ animals show minimal divergence in recombination rate (∼5%). The existence of loci on the X chromosome suggests a genetic mechanism to explain this male-biased evolution. Our results provide an initial map of the genetic changes underlying subspecies differences in genome-scale recombination rate and underscore the power of the house mouse system for understanding the evolution of this trait.     

Sex-specific variation in the genome-wide recombination rate

In most species that reproduce sexually, successful gametogenesis requires recombination during meiosis. The number and placement of crossovers (COs) vary among individuals, with females and males often presenting the most striking contrasts. Despite the recognition that the sexes recombine at different rates (heterochiasmy), existing data fail to answer the question of whether patterns of genetic variation in recombination rate are similar in the two sexes. To fill this gap, we measured the genome-wide recombination rate in both sexes from a panel of wild-derived inbred strains from multiple subspecies of house mice (Mus musculus) and from a few additional species of Mus. To directly compare recombination rates in females and males from the same genetic backgrounds, we applied established methods based on immunolocalization of recombination proteins to inbred strains. Our results reveal discordant patterns of genetic variation in the two sexes. Whereas male genome-wide recombination rates vary substantially among strains, female recombination rates measured in the same strains are more static. The direction of heterochiasmy varies within two subspecies, Mus musculus molossinus and Mus musculus musculus. The direction of sex differences in the length of the synaptonemal complex and CO positions is consistent across strains and does not track sex differences in genome-wide recombination rate. In males, contrasts between strains with high recombination rate and strains with low recombination rate suggest more recombination is associated with stronger CO interference and more double-strand breaks. The sex-specific patterns of genetic variation we report underscore the importance of incorporating sex differences into recombination research.     

A pronounced evolutionary shift of the pseudoautosomal region boundary in house mice

The pseudoautosomal region (PAR) is essential for the accurate pairing and segregation of the X and Y chromosomes during meiosis. Despite its functional significance, the PAR shows substantial evolutionary divergence in structure and sequence between mammalian species. An instructive example of PAR evolution is the house mouse Mus musculus domesticus (represented by the C57BL/6J strain), which has the smallest PAR among those that have been mapped. In C57BL/6J, the PAR boundary is located just ~700 kb from the distal end of the X chromosome, whereas the boundary is found at a more proximal position in Mus spretus, a species that diverged from house mice 2-4 million years ago. In this study we used a combination of genetic and physical mapping to document a pronounced shift in the PAR boundary in a second house mouse subspecies, Mus musculus castaneus (represented by the CAST/EiJ strain), ~430 kb proximal of the M. m. domesticus boundary. We demonstrate molecular evolutionary consequences of this shift, including a marked lineage-specific increase in sequence divergence within Mid1, a gene that resides entirely within the M. m. castaneus PAR but straddles the boundary in other subspecies. Our results extend observations of structural divergence in the PAR to closely related subspecies, pointing to major evolutionary changes in this functionally important genomic region over a short time period.     

Genetic and Evolutionary Correlates of Fine-Scale Recombination Rate Variation in Drosophila persimilis

Recombination is fundamental to meiosis in many species and generates variation on which natural selection can act, yet fine-scale linkage maps are cumbersome to construct. We generated a fine-scale map of recombination rates across two major chromosomes in Drosophila persimilis using 181 SNP markers spanning two of five major chromosome arms. Using this map, we report significant fine-scale heterogeneity of local recombination rates. However, we also observed “recombinational neighborhoods,” where adjacent intervals had similar recombination rates after excluding regions near the centromere and telomere. We further found significant positive associations of fine-scale recombination rate with repetitive element abundance and a 13-bp sequence motif known to associate with human recombination rates. We noted strong crossover interference extending 5–7 Mb from the initial crossover event. Further, we observed that fine-scale recombination rates in D. persimilis are strongly correlated with those obtained from a comparable study of its sister species, D. pseudoobscura. We documented a significant relationship between recombination rates and intron nucleotide sequence diversity within species, but no relationship between recombination rate and intron divergence between species. These results are consistent with selection models (hitchhiking and background selection) rather than mutagenic recombination models for explaining the relationship of recombination with nucleotide diversity within species. Finally, we found significant correlations between recombination rate and GC content, supporting both GC-biased gene conversion (BGC) models and selection-driven codon bias models. Overall, this genome-enabled map of fine-scale recombination rates allowed us to confirm findings of broader-scale studies and identify multiple novel features that merit further investigation.     

Clonal reproduction shapes evolution in the lizard malaria parasite Plasmodium floridense

The preponderant clonal evolution hypothesis (PCE) predicts that frequent clonal reproduction (sex between two clones) in many pathogens capable of sexual recombination results in strong linkage disequilibrium and the presence of discrete genetic subdivisions characterized by occasional gene flow. We expand on the PCE and predict that higher rates of clonal reproduction will result in: (1) morphologically cryptic species that exhibit (2) low within‐species variation and (3) recent between‐species divergence. We tested these predictions in the Caribbean lizard malaria parasite Plasmodium floridense using 63 single‐infection samples in lizards collected from across the parasite's range, and sequenced them at two mitochondrial, one apicoplast, and five nuclear genes. We identified 11 provisionally cryptic species within P. floridense, each of which exhibits low intraspecific variation and recent divergence times between species (some diverged approximately 110,000 years ago). Our results are consistent with the hypothesis that clonal reproduction can profoundly affect diversification of species capable of sexual recombination, and suggest that clonal reproduction may have led to a large number of unrecognized pathogen species. The factors that may influence the rates of clonal reproduction among pathogens are unclear, and we discuss how prevalence and virulence may relate to clonal reproduction.     

Natural selection affects multiple aspects of genetic variation at putatively neutral sites across the human genome

A major question in evolutionary biology is how natural selection has shaped patterns of genetic variation across the human genome. Previous work has documented a reduction in genetic diversity in regions of the genome with low recombination rates. However, it is unclear whether other summaries of genetic variation, like allele frequencies, are also correlated with recombination rate and whether these correlations can be explained solely by negative selection against deleterious mutations or whether positive selection acting on favorable alleles is also required. Here we attempt to address these questions by analyzing three different genome-wide resequencing datasets from European individuals. We document several significant correlations between different genomic features. In particular, we find that average minor allele frequency and diversity are reduced in regions of low recombination and that human diversity, human-chimp divergence, and average minor allele frequency are reduced near genes. Population genetic simulations show that either positive natural selection acting on favorable mutations or negative natural selection acting against deleterious mutations can explain these correlations. However, models with strong positive selection on nonsynonymous mutations and little negative selection predict a stronger negative correlation between neutral diversity and nonsynonymous divergence than observed in the actual data, supporting the importance of negative, rather than positive, selection throughout the genome. Further, we show that the widespread presence of weakly deleterious alleles, rather than a small number of strongly positively selected mutations, is responsible for the correlation between neutral genetic diversity and recombination rate. This work suggests that natural selection has affected multiple aspects of linked neutral variation throughout the human genome and that positive selection is not required to explain these observations.     

Recombination Modulates How Selection Affects Linked Sites in Drosophila

One of the most influential observations in molecular evolution has been a strong association between local recombination rate and nucleotide polymorphisms across the genome. This is interpreted as evidence for ubiquitous natural selection. The alternative explanation, that recombination is mutagenic, has been rejected by the absence of a similar association between local recombination rate and nucleotide divergence between species. However, many recent studies show that recombination rates are often very different even in closely related species, questioning whether an association between recombination rate and divergence between species has been tested satisfactorily. To circumvent this problem, we directly surveyed recombination across approximately 43% of the D. pseudoobscura physical genome in two separate recombination maps and 31% of the D. miranda physical genome, and we identified both global and local differences in recombination rate between these two closely related species. Using only regions with conserved recombination rates between and within species and accounting for multiple covariates, our data support the conclusion that recombination is positively related to diversity because recombination modulates Hill–Robertson effects in the genome and not because recombination is predominately mutagenic. Finally, we find evidence for dips in diversity around nonsynonymous substitutions. We infer that at least some of this reduction in diversity resulted from selective sweeps and examine these dips in the context of recombination rate.     

Molecular phylogeny of Rodentia in the descent of genusBandicota

The interrelationships of murids and other rodent families as well as the evolutionary descent of multiple β-globins of murines are deduced from parsimony trees of relevant globin sequences. Our results show that Caviidae arises first, followed by Sciuridae which joins Muridae. In the murid line of descent Spalacinae arises first followed by two branches, one representing Cricentinae and Arvicolinae and the other Murinae. Although the rates of evolution of globin genes in the different rodent families are different, the murid branches show more or less a uniform rate of evolution of β globins. We have used this information to show that mouse-rat divergence occurred around 20 million years ago. The evolutionary rationale for the presence and the expression of different β-globin genes in murid populations is also discussed. Based on mitochondrial DNA restriction fragment analysis, the between-species relationships ofRattus rattus rufescens, Bandicota indica andBandicota bengalensis have been assessed and the time of divergence of the two bandicoot rats estimated.     

Nucleotide variation in wild and inbred mice

The house mouse is a well-established model organism, particularly for studying the genetics of complex traits. However, most studies of mice use classical inbred strains, whose genomes derive from multiple species. Relatively little is known about the distribution of genetic variation among these species or how variation among strains relates to variation in the wild. We sequenced intronic regions of five X-linked loci in large samples of wild Mus domesticus and M. musculus, and we found low levels of nucleotide diversity in both species. We compared these data to published data from short portions of six X-linked and 18 autosomal loci in wild mice. We estimate that M. domesticus and M. musculus diverged <500,000 years ago. Consistent with this recent divergence, some gene genealogies were reciprocally monophyletic between these species, while others were paraphyletic or polyphyletic. In general, the X chromosome was more differentiated than the autosomes. We resequenced classical inbred strains for all 29 loci and found that inbred strains contain only a small amount of the genetic variation seen in wild mice. Notably, the X chromosome contains proportionately less variation among inbred strains than do the autosomes. Moreover, variation among inbred strains derives from differences between species as well as from differences within species, and these proportions differ in different genomic regions. Wild mice thus provide a reservoir of additional genetic variation that may be useful for mapping studies. Together these results suggest that wild mice will be a valuable complement to laboratory strains for studying the genetics of complex traits.     

Evidence for widespread selection in shaping the genomic landscape during speciation of Populus

Increasing our understanding of how evolutionary processes drive the genomic landscape of variation is fundamental to a better understanding of the genomic consequences of speciation. However, genome‐wide patterns of within‐ and between‐ species variation have not been fully investigated in most forest tree species despite their global ecological and economic importance. Here, we use whole‐genome resequencing data from four Populus species spanning the speciation continuum to reconstruct their demographic histories and investigate patterns of diversity and divergence within and between species. Using Populus trichocarpa as an outgroup species, we further infer the genealogical relationships and estimate the extent of ancient introgression among the three aspen species (Populus tremula, Populus davidiana and Populus tremuloides) throughout the genome. Our results show substantial variation in these patterns along the genomes with this variation being strongly predicted by local recombination rates and the density of functional elements. This implies that the interaction between recurrent selection and intrinsic genomic features has dramatically sculpted the genomic landscape over long periods of time. In addition, our findings provide evidence that, apart from background selection, recent positive selection and long‐term balancing selection have also been crucial components in shaping patterns of genome‐wide variation during the speciation process.     

Estimating divergence time and ancestral effective population size of Bornean and Sumatran orangutan subspecies using a coalescent hidden Markov model

Due to genetic variation in the ancestor of two populations or two species, the divergence time for DNA sequences from two populations is variable along the genome. Within genomic segments all bases will share the same divergence-because they share a most recent common ancestor-when no recombination event has occurred to split them apart. The size of these segments of constant divergence depends on the recombination rate, but also on the speciation time, the effective population size of the ancestral population, as well as demographic effects and selection. Thus, inference of these parameters may be possible if we can decode the divergence times along a genomic alignment. Here, we present a new hidden Markov model that infers the changing divergence (coalescence) times along the genome alignment using a coalescent framework, in order to estimate the speciation time, the recombination rate, and the ancestral effective population size. The model is efficient enough to allow inference on whole-genome data sets. We first investigate the power and consistency of the model with coalescent simulations and then apply it to the whole-genome sequences of the two orangutan sub-species, Bornean (P. p. pygmaeus) and Sumatran (P. p. abelii) orangutans from the Orangutan Genome Project. We estimate the speciation time between the two sub-species to be thousand years ago and the effective population size of the ancestral orangutan species to be , consistent with recent results based on smaller data sets. We also report a negative correlation between chromosome size and ancestral effective population size, which we interpret as a signature of recombination increasing the efficacy of selection.