Find research datasets worth reusing

Data from: Evolution of population genetic structure of the British roe deer by natural and anthropogenic processes (Capreolus capreolus)

Human influence typically impacts on natural populations of conservation interest. These interactions are varied and sometimes complex, and may be negative and unintended or associated with conservation and management strategy. Understanding the details of how these interactions influence and are influenced by natural evolutionary processes is essential to the development of effective conservation strategies. In this study we investigate a species in Britain that has experienced both negative impact through over-hunting in historical times and management efforts through culls and translocations. At the same time, there are regional populations that have been less affected by human influence. We use mtDNA and nuclear microsatellite DNA markers to investigate patterns of connectivity and diversity and find multiple insular populations in Britain that likely evolved within the Holocene (when the habitat was free of ice). We identify three concurrent processes. First, surviving indigenous populations show highly provincial patterns of philopatry, maintaining and generating population structure on a small geographic scale. Second, founder populations into habitat extirpated of native populations have expanded but remained largely insular. Third, introductions into established populations generate some admixture. We discuss the implications for the evolution of diversity of the integration of natural processes with anthropogenic influences on population size and distribution.

opencc-zeroDec 2013View details →

Data from: An investigation for population maintenance mechanism in a miniature garden: genetic connectivity or independence of small islet populations of the Ryukyu five-lined skink

The Ryukyu five-lined skink (Plestiodon marginatus) is an island lizard that is even found in tiny islets with less than half a hectare of habitat area. We hypothesized that the island populations are maintained under frequent gene flow among the islands or independent of each other. To test our hypotheses, we investigated genetic structure of 21 populations from 11 land-bridge islands that were connected during the latest glacial age, and four isolated islands. Analyses using mitochondrial cytochrome b gene sequence (n = 67) and 10 microsatellite loci (n = 235) revealed moderate to high levels of genetic differentiation, existence of many private alleles/haplotypes in most islands, little contemporary migration, a positive correlation between genetic variability and island area, and a negative correlation between relatedness and island area. These evidences suggest a strong effect of independent genetic drift as opposed to gene flow, favoring the isolation hypothesis even in tiny islet populations. Isolation-by-distance effect was demonstrated and it became more prominent when the four isolated islands were excluded, suggesting that the pattern is a remnant of the land-bridge age. In a few island populations, however, the possibility of occasional overwater dispersals was partially supported and therefore could not be ruled out.

opencc-zeroDec 2010View details →

Data from: The effect of variable frequency of sexual reproduction on the genetic structure of natural populations of a cyclical parthenogen

Cyclical parthenogens are a valuable system in which to empirically test theoretical predictions as to the genetic consequences of sexual reproduction in natural populations, particularly if the frequency of sexual relative to asexual reproduction can be quantified. In this study we utilized a series of lake populations of the cyclical parthenogen, Daphnia pulicaria, that vary consistently in their investment in sexual reproduction, to address the questions of whether the ecological variation in investment in sex is detectable at the genetic level, and if so, whether the genetic patterns seen are consistent with theoretical predictions. We show that there is variation in the genetic structure of these populations in a manner consistent with their investment in sexual reproduction. Populations engaging in a high frequency of sex were in Hardy-Weinberg and gametic phase equilibrium, and showed little genotypic differentiation across sampled years. In contrast, populations with a lower frequency of sex deviated widely from equilibrium, had reduced multi-locus clonal diversity, and showed significant temporal genotypic deviation.

Data from: Influence of introduction history on genetic variation in introduced populations: a case study of Oregon Chub

Population introductions and reintroductions have become a common tool for conserving threatened species, but oftentimes introduced populations have reduced the genetic diversity compared with the source population they were founded from. Population introductions played an important role in the recovery of the Oregon Chub Oregonichthys crameri, a small floodplain minnow found in western Oregon. Unlike many introduction efforts, introduced populations of Oregon Chub were founded using large numbers of individuals (hundreds in many cases) and each population had a unique introduction history (e.g., number of founders, source populations selected, duration of the introduction effort). We used microsatellite loci to examine 13 introduced populations and their respective sources to evaluate how well the introduction program captured genetic diversity present in the wild populations. Genetic variation was reduced by roughly 25% in one introduced population, and three introduced populations showed evidence of a genetic bottleneck due to heterozygote excess. Populations introduced from multiple sources had greater genetic diversity than populations from a single source. When multiple source populations were used, all source populations contributed genetic material to the introduced population, though the proportional contribution from each source population varied. Using correlation analyses and general linear models, we explored the relationship between introduction history variables and genetic diversity. Our top-ranked models included genetic diversity in the source population, and this variable had the highest variable importance weight (0.999), but the number of founders and the number of source populations were also important. Overall, the Oregon Chub introduction program was highly successful at capturing the genetic variation observed in natural populations. Results of this study will be useful for planning future population introductions for Oregon Chub and other species of conservation concern.

Data from: Parasite genetic distance and local adaptation in coevolving bacteria-bacteriophage populations

Antagonistic coevolution between hosts and parasites can lead to local adaptation (LA), such that parasite fitness is greatest in sympatric hosts (or vice versa). The magnitude of LA typically increases with geographic distance, which is assumed to be because genetic (and hence phenotypic) distance increases with geographic distance. Here we explicitly test the relationships between parasite genetic and phenotypic distance and LA using isolates of coevolved viral parasites (lytic bacteriophage ϕ2) and the host bacterium Pseudomonas fluorescens SBW25. We find positive relationships between parasite genotype and infectivity phenotype, but the strength of the relationship was greater when infectivity was defined by the identity of hosts that could be infected rather than the actual number of hosts infected (host range), and when measurements were compared within rather than among populations. Crucially, we find a monotonic relationship between LA and genetic distance across phage isolates from different populations, although in contrast to many geographic studies, parasite LA decreased with genetic distance. These results can be explained by the fact that bacteria can rapidly adapt to phage infectivity mutations, but that evolved resistance has a degree of specificity to the local phage population. Our results show that antagonistic coevolution alone can result in predictable links between genetic distance and host-parasite local adaptation.

Data from: Post-glacial northward expansion and genetic differentiation between migratory and sedentary populations of the broad-tailed hummingbird (Selasphorus platycercus)

Unlike other migratory hummingbirds in North America, the broad-tailed hummingbird (Selasphorus platycercus) exhibits both long-distance migratory behaviour in the USA and sedentary behaviour in Mexico and Guatemala. We examined the evolution of migration linked to its northward expansion using a multi-perspective approach. We analysed variation in morphology, mitochondrial and nuclear DNA, estimated migration rates between migratory and sedentary populations, compared divergence times with the occurrence of Quaternary climate events, and constructed species distribution models to predict where migratory and sedentary populations resided during the Last Glacial Maximum (LGM) and Last Interglacial (LIG) events. Our results are consistent with a recent northward population expansion driven by migration from southern sedentary populations. Phylogeographical analyses and population genetics methods revealed that migratory populations in the USA and sedentary populations in Mexico of the platycercus subspecies form one admixed population, and that sedentary populations from southern Mexico and Guatemala (guatemalae) undertook independent evolutionary trajectories. Species distribution modelling revealed that the species is a niche tracker and that the climate conditions associated with modern obligate migrants in the USA were not present during the LIG, which provides indirect evidence for recent migratory behaviour in broad-tailed hummingbirds on the temporal scale of glacial cycles. The finding that platycercus hummingbirds form one genetic population and that suitable habitat for migratory populations was observed in eastern Mexico during the LIG also suggests that the conservation of overwintering sites is crucial for obligate migratory populations currently facing climate change effects.

Data from: A comparison of single-sample estimators of effective population sizes from genetic marker data

In molecular ecology and conservation genetics studies, the important parameter of effective population size (Ne) is increasingly estimated from a single sample of individuals taken at random from a population and genotyped at a number of marker loci. Several estimators are developed, based on the information of linkage disequilibrium (LD), heterozygote excess (HE), molecular coancestry (MC) and sibship frequency (SF) in marker data. The most popular is the LD estimator, because it is more accurate than HE and MC estimators and is simpler to calculate than SF estimator. However, little is known about the accuracy of LD estimator relative to that of SF and about the robustness of all single-sample estimators when some simplifying assumptions (e.g. random mating, no linkage, no genotyping errors) are violated. This study fills the gaps and uses extensive simulations to compare the biases and accuracies of the four estimators for different population properties (e.g. bottlenecks, nonrandom mating, haplodiploid), marker properties (e.g. linkage, polymorphisms) and sample properties (e.g. numbers of individuals and markers) and to compare the robustness of the four estimators when marker data are imperfect (with allelic dropouts). Extensive simulations show that SF estimator is more accurate, has a much wider application scope (e.g. suitable to nonrandom mating such as selfing, haplodiploid species, dominant markers) and is more robust (e.g. to the presence of linkage and genotyping errors of markers) than the other estimators. An empirical data set from a Yellowstone grizzly bear population was analysed to demonstrate the use of the SF estimator in practice.

opencc-zeroDec 2011View details →

Data from: Large-scale natural disturbance alters genetic population structure of the sailfin molly, Poecilia latipinna

Many inferences about contemporary rates of gene flow are based on the assumption that the observed genetic structure among populations is stable. Recent studies have uncovered several cases in which this assumption is tenuous. Most of those studies have focused on the effects that regular environmental fluctuations can have on genetic structure and gene flow patterns. Occasional catastrophic disturbances could also alter either the distribution of habitat or the spatial distribution of organisms in a way that affects population structure. However, evidence of such effects is sparse in the literature because it is difficult to obtain. Hurricanes, in particular, have the potential to exert dramatic effects on population structure of organisms found on islands or coral reefs or in near shore and coastal habitats. Here we draw on a historic genetic data set and new data to suggest that the genetic structure of sailfin molly (Poecilia latipinna) populations in north Florida was altered dramatically by an unusually large and uncommon type of storm surge associated with Hurricane Dennis in 2005. We compare the spatial pattern of genetic variation in these populations after Hurricane Dennis to the patterns described in an earlier study in this same area. We use comparable genetic data from another region of Florida, collected in the same two periods, to estimate the amount of change expected from typical temporal variation in population structure. The comparative natural history of sailfin mollies in these two regions indicates that the change in population structure produced by the storm surge is not the result of many local extinctions with recolonization from a few refugia but emerged from a pattern of mixing and redistribution.

opencc-zeroDec 2014View details →

Data from: The red coral (Corallium rubrum) transcriptome: a new resource for population genetics and local adaptation studies

The question of species survival and evolution in heterogeneous environments has long been a subject for study. Indeed, it is often difficult to identify the molecular basis of adaptation to contrasted environments, and nongenetic effects increase the difficulty to disentangle fixed effects, such as genetic adaptation, from variable effects, such as individual phenotypic plasticity, in adaptation. Nevertheless, this question is also of great importance for understanding the evolution of species in a context of climate change. The red coral (Corallium rubrum) lives in the Mediterranean Sea, where at depths ranging from 5 to 600 m, it meets very contrasted thermal conditions. The shallowest populations of this species suffered from mortality events linked with thermal anomalies that have highlighted thermotolerance differences between individuals. We provide here a new transcriptomic resource, as well as candidate markers for the study of local adaptation. We sequenced the transcriptome of six individuals from 5 m and six individuals from 40 m depth at the same site of the Marseilles bay, after a period of common garden acclimatization. We found differential expression maintained between the two depths even after common garden acclimatization, and we analysed the polymorphism pattern of these samples. We highlighted contigs potentially implicated in the response to thermal stress, which could be good candidates for the study of thermal adaptation for the red coral. Some of these genes are also involved in the response to thermal stress in other corals. Our method enables the identification of candidate loci of local adaptation useful for other nonmodel organisms.

opencc-zeroDec 2018View details →

Data from: A population genetic assessment of taxonomic species: the case of Lake Malawi cichlid fishes

Organisms sampled for population level research are typically assigned to species by morphological criteria. But if those criteria are limited to one sex or life stage, or the organisms come from a complex of closely related forms, the species assignments may misdirect analyses. The impact of such sampling can be assessed from the correspondence of genetic clusters, identified only from patterns of genetic variation, to the species identified using only phenotypic criteria. We undertook this protocol with the rock-dwelling mbuna cichlids of Lake Malawi, for which species within genera are usually identified by investigators using adult male coloration patterns. Given high local endemism of male color patterns, and considerable allele sharing among species, there persists considerable taxonomic uncertainty in these fishes. Over 700 individuals from a single transect were photographed, genotyped, and separately assigned: (1) to morphospecies using photographs; and (2) to genetic clusters using five widely used methods. Overall, the correspondence between clustering methods was strong for larger clusters, but methods varied widely in estimated number of clusters. The correspondence between morphospecies and genetic clusters was also strong for larger clusters, as well as some smaller clusters for some methods. These analyses generally affirm (1) adult male-limited sampling and (2) the taxonomic status of Lake Malawi mbuna, as the species in our study largely appear to be well-demarcated genetic entities. More generally, our analyses highlight the challenges for clustering methods when the number of populations is unknown, especially in cases of highly uneven sample sizes.

opencc-zeroDec 2016View details →

Data from: Sex, size and timing: sampling design for reliable population genetics analyses using microsatellite data

1. Population genetics is used in a wide variety of fields such as ecology and biodiversity conservation. How estimated genetic characteristics of natural populations can be influenced by the sampling design has been a long-standing concern. Multiple simulation and empirical studies illustrated the influence of both sample size and polymorphism of markers. However, our review of studies on butterfly population genetics indicates no consensus on sample size for the estimation of genetic diversity or differentiation. Furthermore, other aspects of sampling design (sex ratio and timing of sampling) were not addressed and their potential impact on genetic parameter estimates rarely explored. 2. Using a large empirical dataset (with spatial and temporal replicates) collected on a butterfly species, Boloria aquilonaris, as well as simulated datasets reflecting (1) three scenarios of migration-genetic drift equilibrium and (2) one scenario of parameter stabilization after 100,000 generations, we quantified the impacts of three aspects of genetic sampling design (namely sample size, sex ratio, and timing of sampling) on the estimation of allele frequencies and its potential downstream impact on the estimation of genetic parameters. 3. With empirical data, we found that sample size and timing of sampling strongly affected the accuracy of allele frequencies and the downstream analyses, while sex ratio did not. Our results were consistent across spatial and temporal replicates. Also, with simulated data, we showed that the genetic sampling design had limited effect in systems where dispersal outweighs genetic drift, while it can have major consequences on our understanding of the genetic diversity and population differentiation in systems dominated by genetic drift (such as most study systems with conservation concerns). 4. We advocate for careful consideration of all aspects of the sampling design in population genetics studies, i.e. a sufficient number of samples, while ensuring similar sex ratio among sampling locations and collecting with timing appropriate to the question under study. This is particularly important when the study aims at species conservation.

opencc-zeroDec 2010View details →

Data from: Elevated genetic diversity of mitochondrial genes in asexual populations of bark lice ("Psocoptera": Echmepteryx hageni)

Asexual reproduction is commonly thought to be associated with low genetic diversity in animals. Echmepteryx hageni (Insecta: "Psocoptera") is one of several psocopteran species that are primarily parthenogenetic, but also exists in small, isolated sexual populations. We used mitochondrial DNA sequences to investigate the population history and genealogical relationships between the sexual and asexual forms of this species. The asexual population of E. hageni exhibits extremely high mitochondrial haplotype diversity (H = 0.98), whereas the sexual forms had significantly lower haplotypic diversity (H = 0.25, after correcting for sample size). This diversity in asexuals represents one the greatest genetic diversities reported for asexual animals in the literature. Nucleotide diversities were also higher in asexual compared to sexual populations (π = 0.0071 vs. 0.00027). Compared to other reported estimates of π in insects, asexual nucleotide diversity is high, but not remarkably elevated. Three hypotheses might explain the elevated genetic diversity of asexual populations: 1) larger effective population size, 2) greater mutation rate, or 3) possible recent origin of sexuals. In addition, phylogeographic analysis revealed little geographic structure among asexual E. hageni, although specimens from the upper Midwest form a single clade and are genetically differentiated. The mismatch distribution and neutrality tests indicate a historical population size increase, possibly associated with expansion from glacial refugia.

Data from: Multivariate selection and intersexual genetic constraints in a wild bird population

When traits are genetically correlated between the sexes, the response to selection in one sex can be altered by indirect selection in the other sex, a type of genetic constraint commonly referred to as intralocus sexual conflict (ISC). While potentially common, ISC has rarely been studied in wild populations. In this study, we applied a multivariate framework to quantify the microevolutionary impacts of ISC over a set of morphological traits (wing length, tarsus length, bill depth, and bill length) in a wild population of great tits (Parus major) from Wytham Woods, UK. Specifically, we quantified the impact of cross-sex genetic covariances (the B matrix) on the additive genetic variance for relative fitness expected to be generated by directional selection and additive genetic (co)variance. Together, multivariate sex-specific selection and additive genetic (co)variance were expected to generate additive genetic variance for relative fitness that was uncorrelated between the sexes (cross-sex genetic correlation = -0.003, 95% CI = -0.83, 0.83). Gender load, defined as the expected reduction in additive genetic variance for relative fitness generated by the traits under study due to sex-specific effects, was estimated at 50% (95% CI = 13%, 86%). This study provides novel insights into the evolution of sexual dimorphism in great tits and illustrates how quantitative genetics and selection analyses can be combined in a multivariate framework to quantify the expected microevolutionary impacts of ISC.

Data from: Population size and time since island isolation determine genetic diversity loss in insular frog populations

Understanding the factors that contribute to loss of genetic diversity in fragmented populations is crucial for conservation measurements. Land-bridge archipelagoes offer ideal model systems for identifying the long-term effects of these factors on genetic variations in wild populations. In this study, we used 9 microsatellite markers to quantify genetic diversity and differentiation of 810 pond frogs (Pelophylax nigromaculataus) from 24 islands of the Zhoushan Archipelago and 3 sites on nearby mainland China and estimated the effects of the island area, population size, time since island isolation, distance to the mainland and distance to the nearest larger island on reduced genetic diversity of insular populations. The mainland populations displayed higher genetic diversity than insular populations. Genetic differentiations and no obvious gene flow were detected among the frog populations on the islands. Hierarchical partitioning analysis showed that only time since island isolation (square root transformed) and population size (log transformed) significantly contributed to insular genetic diversity. These results suggest that decreased genetic diversity and genetic differentiations among insular populations may have been caused by random genetic drift following isolation by rising sea levels during the Holocene. The results provide strong evidence for a relationship between retained genetic diversity and population size and time since island isolation for pond frogs on the islands, consistent with the prediction of the neutral theory for finite populations. Our study highlights the importance of the size and estimated isolation time of populations in understanding the mechanisms of genetic diversity loss and differentiation in fragmented wild populations.

Data from: Divergent selection on, but no genetic conflict over, female and male timing and rate of reproduction in a human population

The sexes often have different phenotypic optima for important life-history traits, and because of a largely shared genome this can lead to a conflict over trait expression. In mammals, the obligate costs of reproduction are higher for females, making reproductive timing and rate especially liable to conflict between the sexes. While studies from wild vertebrates support such sexual conflict, it remains unexplored in humans. We used a pedigreed human population from preindustrial Finland to estimate sexual conflict over age at first and last reproduction, reproductive lifespan and reproductive rate. We found that the phenotypic selection gradients differed between the sexes. We next established significant heritabilities in both sexes for all traits. All traits, except reproductive rate, showed strongly positive intersexual genetic correlations and were strongly genetically correlated with fitness in both sexes. Moreover, the genetic correlations with fitness were almost identical in men and women. For reproductive rate, the intersexual correlation and the correlation with fitness were weaker but again similar between the sexes. Thus, in this population, an apparent sexual conflict at the phenotypic level did not reflect an underlying genetic conflict over the studied reproductive traits. These findings emphasize the need for incorporating genetic perspectives into studies of human life-history evolution.

Data from: Multiple processes drive genetic structure of humpback whale (Megaptera novaeangliae) populations across spatial scales

Elucidating patterns of population structure for species with complex life histories, and disentangling the processes driving such patterns, remains a significant analytical challenge. Humpback whale (Megaptera novaeangliae) populations display complex genetic structures that have not been fully resolved at all spatial scales. We generated a data set of nuclear markers for 3,575 samples spanning the seven breeding stocks and substocks found in the South Atlantic and western and northern Indian Oceans. For the total sample, and males and females separately, we assessed genetic diversity, tested for genetic differentiation between putative populations and isolation by distance, estimated the number of genetic clusters without a priori population information, and estimated rates of gene flow using maximum likelihood and Bayesian approaches. At the ocean basin scale, structure is governed by geographic distance (IBD p<0.05) and female fidelity to breeding areas, in line with current understanding of the drivers of broad-scale population structure. Consistent with previous studies, the Arabian Sea breeding stock was highly genetically differentiated (FST 0.034-0.161; p<0.01 for all comparisons). However, the breeding stock boundary between west South Africa and east Africa was more porous than expected based on genetic differentiation, cluster, and gene flow analyses. Instances of male-fidelity to breeding areas and relatively high rates of dispersal for females were also observed between the three substocks in the western Indian Ocean. This mismatch between demographic units and current management boundaries may have ramifications for assessments of the status and continued protections of populations still in recovery from commercial whaling.

opencc-zeroDec 2016View details →

Data from: Population genomic analysis suggests strong influence of river network on spatial distribution of genetic variation in invasive saltcedar across the southwestern US

Understanding the complex influences of landscape and anthropogenic elements that shape the population genetic structure of invasive species provides insight into patterns of colonization and spread. The application of landscape genomics techniques to these questions may offer detailed, previously undocumented insights into factors influencing species invasions. We investigated the spatial pattern of genetic variation and the influences of landscape factors on population similarity in the invasive riparian shrub saltcedar (Tamarix L.) by analyzing 1,997 genome-wide SNP markers for 259 individuals from 25 populations collected throughout the southwestern US. Our results revealed a broad-scale spatial genetic differentiation of saltcedar populations between the Colorado and Rio Grande river basins and identified potential barriers to population similarity along both river systems. River pathways most strongly contributed to population similarity. In contrast, low temperature and dams likely served as barriers to population similarity. We hypothesize that large-scale geographic patterns in genetic diversity resulted from a combination of early introductions from distinct populations, the subsequent influence of natural selection, dispersal barriers, and founder effects during range expansion.

opencc-zeroDec 2013View details →

Data from: Indirect genetic effects underlie oxygen-limited thermal tolerance within a coastal population of chinook salmon

With global temperatures projected to surpass the limits of thermal tolerance for many species, evaluating the heritable variation underlying thermal tolerance is critical for understanding the potential for adaptation to climate change. We examined the evolutionary potential of thermal tolerance within a population of chinook salmon (Oncorhynchus tshawytscha) by conducting a full-factorial breeding design and measuring the thermal performance of cardiac function and the critical thermal maximum (CTmax) of offspring from each family. Additive genetic variation in offspring phenotype was mostly negligible, although these direct genetic effects explained 53% of the variation in resting heart rate (fH). Conversely, maternal effects had a significant influence on resting fH, scope for fH, cardiac arrhythmia temperature and CTmax. These maternal effects were associated with egg size, as indicated by strong relationships between the mean egg diameter of mothers and offspring thermal tolerance. Because egg size can be highly heritable in chinook salmon, our finding indicates that the maternal effects of egg size constitute an indirect genetic effect contributing to thermal tolerance. Such indirect genetic effects could accelerate evolutionary responses to the selection imposed by rising temperatures and could contribute to the population-specific thermal tolerance that has recently been uncovered among Pacific salmon populations.

Data from: Conserved genetic architecture underlying recombination rate variation in a wild population of Soay sheep (Ovis aries)

Meiotic recombination breaks down linkage disequilibrium and forms new haplotypes, meaning that it is an important driver of diversity in eukaryotic genomes. Understanding the causes of variation in recombination rate is important in interpreting and predicting evolutionary phenomena and for understanding the potential of a population to respond to selection. However, despite attention in model systems, there remains little data on how recombination rate varies at the individual level in natural populations. Here, we used extensive pedigree and high-density SNP information in a wild population of Soay sheep (Ovis aries) to investigate the genetic architecture of individual autosomal recombination rate. Individual rates were high relative to other mammal systems, and were higher in males than in females (autosomal map lengths of 3748 cM and 2860 cM, respectively). The heritability of autosomal recombination rate was low but significant in both sexes (h2 = 0.16 & 0.12 in females and males, respectively). In females, 46.7% of the heritable variation was explained by a sub-telomeric region on chromosome 6; a genome-wide association study showed the strongest associations at the locus RNF212, with further associations observed at a nearby ~374kb region of complete linkage disequilibrium containing three additional candidate loci, CPLX1, GAK and PCGF3. A second region on chromosome 7 containing REC8 and RNF212B explained 26.2% of the heritable variation in recombination rate in both sexes. Comparative analyses with 40 other sheep breeds showed that haplotypes associated with recombination rates are both old and globally distributed. Both regions have been implicated in rate variation in mice, cattle and humans, suggesting a common genetic architecture of recombination rate variation in mammals.