Variation in mutation, recombination, and transposition rates in Drosophila melanogaster and Drosophila simulans.

The rates of mutation, recombination, and transposition are core parameters in models of evolution. They impact genetic diversity, responses to ongoing selection, and levels of genetic load. However, even for key evolutionary model species such as Drosophila melanogaster and Drosophila simulans, few estimates of these parameters are available, and we have little idea of how rates vary between individuals, sexes, or populations. Knowledge of this variation is fundamental for parameterizing models of genome evolution. Here, we provide direct estimates of mutation, recombination, and transposition rates and their variation in a West African and a European population of D. melanogaster and a European population of D. simulans Across 89 flies, we observe 58 single-nucleotide mutations, 286 crossovers, and 89 transposable element (TE) insertions. Compared to the European D. melanogaster, we find the West African population has a lower mutation rate (1.67 × 10-9 site-1 gen-1 vs. 4.86 × 10-9 site-1 gen-1) and a lower transposition rate (8.99 × 10-5 copy-1 gen-1 vs. 23.36 × 10-5 copy-1 gen-1), but a higher recombination rate (3.44 cM/Mb vs. 2.06 cM/Mb). The European D. simulans population has a similar mutation rate to European D. melanogaster, but a significantly higher recombination rate and a lower, but not significantly different, transposition rate. Overall, we find paternal-derived mutations are more frequent than maternal ones in both species. Our study quantifies the variation in rates of mutation, recombination, and transposition among different populations and sexes, and our direct estimates of these parameters in D. melanogaster and D. simulans will benefit future studies in population and evolutionary genetics.

Understanding the Genetic Basis of Variation in Meiotic Recombination: Past, Present, and Future.

Meiotic recombination is a fundamental feature of sexually reproducing species. It is often required for proper chromosome segregation and plays important role in adaptation and the maintenance of genetic diversity. The molecular mechanisms of recombination are remarkably conserved across eukaryotes, yet meiotic genes and proteins show substantial variation in their sequence and function, even between closely related species. Furthermore, the rate and distribution of recombination shows a huge diversity within and between chromosomes, individuals, sexes, populations, and species. This variation has implications for many molecular and evolutionary processes, yet how and why this diversity has evolved is not well understood. A key step in understanding trait evolution is to determine its genetic basis-that is, the number, effect sizes, and distribution of loci underpinning variation. In this perspective, I discuss past and current knowledge on the genetic basis of variation in recombination rate and distribution, explore its evolutionary implications, and present open questions for future research.

Genetic architecture of inbreeding depression may explain its persistence in a population of wild red deer.

Inbreeding depression is of major concern in declining populations, but relatively little is known about its genetic architecture in wild populations, such as the degree to which it is composed of large or small effect loci and their distribution throughout the genome. Here, we combine fitness and genomic data from a wild population of red deer to investigate the genomic distribution of inbreeding effects. Based on the runs of homozygosity (ROH)-based inbreeding coefficient, FROH, we use chromosome-specific inbreeding coefficients (FROHChr) to explore whether the effect of inbreeding varies between chromosomes. Under the assumption that within an individual the probability of being identical-by-descent is equal across all chromosomes, we used a multi-membership model to estimate the deviation of FROHChr from the average inbreeding effect. This novel approach ensures effect sizes are not overestimated whilst maximising the power of our available dataset of >3000 individuals genotyped on >35,000 autosomal SNPs. We find that most chromosomes confer a minor reduction in fitness-related traits, which when these effects are summed, results in the observed inbreeding depression in birth weight, survival and lifetime breeding success. However, no chromosomes had a significant detrimental effect compared to the overall effect of inbreeding, indicating no major effect loci. We conclude that in this population, inbreeding depression is likely the result of multiple mildly or moderately deleterious mutations spread across all chromosomes, which are difficult to detect with statistical confidence. Such mutations will be inefficiently purged, which may explain the persistence of inbreeding depression in this population.

Selection, recombination and population history effects on runs of homozygosity (ROH) in wild red deer (Cervus elaphus).

The distribution of runs of homozygosity (ROH) may be shaped by a number of interacting processes such as selection, recombination and population history, but little is known about the importance of these mechanisms in shaping ROH in wild populations. We combined an empirical dataset of >3000 red deer genotyped at >35,000 genome-wide autosomal SNPs and evolutionary simulations to investigate the influence of each of these factors on ROH. We assessed ROH in a focal and comparison population to investigate the effect of population history. We investigated the role of recombination using both a physical map and a genetic linkage map to search for ROH. We found differences in ROH distribution between both populations and map types indicating that population history and local recombination rate have an effect on ROH. Finally, we ran forward genetic simulations with varying population histories, recombination rates and levels of selection, allowing us to further interpret our empirical data. These simulations showed that population history has a greater effect on ROH distribution than either recombination or selection. We further show that selection can cause genomic regions where ROH is common, only when the effective population size (Ne) is large or selection is particularly strong. In populations having undergone a population bottleneck, genetic drift can outweigh the effect of selection. Overall, we conclude that in this population, genetic drift resulting from a historical population bottleneck is most likely to have resulted in the observed ROH distribution, with selection possibly playing a minor role.

Contemporary selection on MHC genes in a free-living ruminant population.

Genes within the major histocompatibility complex (MHC) are the most variable identified in vertebrates. Pathogen-mediated selection is believed to be the main force maintaining MHC diversity. However, relatively few studies have demonstrated contemporary selection on MHC genes. Here, we examine associations between MHC variation and several fitness measurements including total fitness and five fitness components, in 3400 wild Soay sheep (Ovis aries) monitored between 1989 and 2012. In terms of total fitness, measured as lifetime breeding success of all individuals born, we found haplotypes named C and D were associated with decreased and increased male total fitness respectively. In terms of fitness components, juvenile survival was associated with haplotype divergence while individual haplotypes (C, D and F) were associated with adult fitness components. Consistent with the increased male total fitness, the rarest haplotype D has increased in frequency throughout the study period more than expected under neutral expectations. Our results demonstrate that contemporary natural selection is acting on MHC class II genes in Soay sheep and that the mode of selection on specific fitness components can be different mode from selection on total fitness.

Taking quantitative genomics into the wild.

We organized this special issue to highlight new work and review recent advances at the cutting edge of 'wild quantitative genomics'. In this editorial, we will present some history of wild quantitative genetic and genomic studies, before discussing the main themes in the papers published in this special issue and highlighting the future outlook of this dynamic field.

Genomic analysis reveals a polygenic architecture of antler morphology in wild red deer (Cervus elaphus).

Sexually selected traits show large variation and rapid evolution across the animal kingdom, yet genetic variation often persists within populations despite apparent directional selection. A key step in solving this long-standing paradox is to determine the genetic architecture of sexually selected traits to understand evolutionary drivers and constraints at the genomic level. Antlers are a form of sexual weaponry in male red deer (Cervus elaphus). On the island of Rum, Scotland, males with larger antlers have increased breeding success, yet there has been no evidence of any response to selection at the genetic level. To try and understand the mechanisms underlying this observation, we investigate the genetic architecture of ten antler traits and their principal components using genomic data from >38,000 SNPs. We estimate the heritabilities and genetic correlations of the antler traits using a genomic relatedness approach. We then use genome-wide association and haplotype-based regional heritability to identify regions of the genome underlying antler morphology, and an empirical Bayes approach to estimate the underlying distributions of allele effect sizes. We show that antler morphology is highly repeatable over an individual's lifetime, heritable and has a polygenic architecture and that almost all antler traits are positively genetically correlated with some loci identified as having pleiotropic effects. Our findings suggest that a large mutational target and genetic covariances among antler traits, in part maintained by pleiotropy, are likely to contribute to the maintenance of genetic variation in antler morphology in this population.

Associations between MHC class II variation and phenotypic traits in a free-living sheep population.

Pathogen-mediated selection (PMS) is thought to maintain the high level of allelic diversity observed in the major histocompatibility complex (MHC) class II genes. A comprehensive way to demonstrate contemporary selection is to examine associations between MHC variation and individual fitness. As individual fitness is hard to measure, many studies examine associations between MHC variation and phenotypic traits, including direct or indirect measures of adaptive immunity thought to contribute to fitness. Here, we tested associations between MHC class II variation and five phenotypic traits measured in free-living sheep captured in August: weight, strongyle faecal egg count, and plasma IgA, IgE and IgG immunoglobulin titres against the gastrointestinal nematode parasite Teladorsagia circumcincta. We found no association between MHC class II variation and weight or strongyle faecal egg count. We did, however, find associations between MHC class II variation and immunoglobulin levels which varied with isotype, age and sex. Our results suggest associations between MHC and phenotypic traits are more likely to be found for traits more closely associated with pathogen defence than integrative traits such as bodyweight and highlight the association between MHC variation and antibodies in wild populations.

Genomic prediction in the wild: A case study in Soay sheep.

Genomic prediction, the technique whereby an individual's genetic component of their phenotype is estimated from its genome, has revolutionised animal and plant breeding and medical genetics. However, despite being first introduced nearly two decades ago, it has hardly been adopted by the evolutionary genetics community studying wild organisms. Here, genomic prediction is performed on eight traits in a wild population of Soay sheep. The population has been the focus of a >30 year evolutionary ecology study and there is already considerable understanding of the genetic architecture of the focal Mendelian and quantitative traits. We show that the accuracy of genomic prediction is high for all traits, but especially those with loci of large effect segregating. Five different methods are compared, and the two methods that can accommodate zero-effect and large-effect loci in the same model tend to perform best. If the accuracy of genomic prediction is similar in other wild populations, then there is a real opportunity for pedigree-free molecular quantitative genetics research to be enabled in many more wild populations; currently the literature is dominated by studies that have required decades of field data collection to generate sufficiently deep pedigrees. Finally, some of the potential applications of genomic prediction in wild populations are discussed.

Vitamin D status is heritable and under environment-dependent selection in the wild.

Vitamin D has a well-established role in skeletal health and is increasingly linked to chronic disease and mortality in humans and companion animals. Despite the clear significance of vitamin D for health and obvious implications for fitness under natural conditions, no longitudinal study has tested whether the circulating concentration of vitamin D is under natural selection in the wild. Here, we show that concentrations of dietary-derived vitamin D2 and endogenously produced vitamin D3  metabolites are heritable and largely polygenic in a wild population of Soay sheep (Ovis aries). Vitamin D2  status was positively associated with female adult survival, and vitamin D3  status predicted female fecundity in particular, good environment years when sheep density and competition for resources was low. Our study provides evidence that vitamin D status has the potential to respond to selection, and also provides new insights into how vitamin D metabolism is associated with fitness in the wild.

Recombination rates in pigs differ between breeds, sexes and individuals, and are associated with the RNF212, SYCP2, PRDM7, MEI1 and MSH4 loci.

BACKGROUND: Recombination is a fundamental part of mammalian meiosis that leads to the exchange of large segments of DNA between homologous chromosomes and is therefore an important driver of genetic diversity in populations. In breeding populations, understanding recombination is of particular interest because it can break up unfavourable linkage phases between alleles and produce novel combinations of alleles that could be exploited in selection. In this study, we used dense single nucleotide polymorphism (SNP) genotype data and pedigree information to analyse individual and sex-specific variation and genetic architecture of recombination rates within and between five commercially selected pig breeds. RESULTS: In agreement with previous studies, recombination rates were higher in females than in males for all breeds and for all chromosomes, except 1 and 13, for which male rates were slightly higher. Total recombination rate differed between breeds but the pattern of recombination along the chromosomes was well conserved across breeds for the same sex. The autosomal linkage maps spanned a total length of 1731 to 1887 cM for males and of 2231 to 2515 cM for females. Estimates of heritability for individual autosomal crossover count ranged from 0.04 to 0.07 for males and from 0.08 to 0.11 for females. Fourteen genomic regions were found to be associated with individual autosomal crossover count. Of these, four were close to or within candidate genes that have previously been associated with individual recombination rates in pigs and other mammals, namely RNF212, SYCP2 and MSH4. Two of the identified regions included the PRDM7 and MEI1 genes, which are known to be involved in meiosis but have not been previously associated with variation in individual recombination rates. CONCLUSIONS: This study shows that genetic variation in autosomal recombination rate persists in domesticated species under strong selection, with differences between closely-related breeds and marked differences between the sexes. Our findings support results from other studies, i.e., that individual crossover counts are associated with the RNF212, SYCP2 and MSH4 genes in pig. In addition, we have found two novel candidate genes associated with the trait, namely PRDM7 and MEI1.

The genetic architecture of helminth-specific immune responses in a wild population of Soay sheep (Ovis aries).

Much of our knowledge of the drivers of immune variation, and how these responses vary over time, comes from humans, domesticated livestock or laboratory organisms. While the genetic basis of variation in immune responses have been investigated in these systems, there is a poor understanding of how genetic variation influences immunity in natural, untreated populations living in complex environments. Here, we examine the genetic architecture of variation in immune traits in the Soay sheep of St Kilda, an unmanaged population of sheep infected with strongyle gastrointestinal nematodes. We assayed IgA, IgE and IgG antibodies against the prevalent nematode Teladorsagia circumcincta in the blood plasma of > 3,000 sheep collected over 26 years. Antibody levels were significantly heritable (h2 = 0.21 to 0.57) and highly stable over an individual's lifespan. IgA levels were strongly associated with a region on chromosome 24 explaining 21.1% and 24.5% of heritable variation in lambs and adults, respectively. This region was adjacent to two candidate loci, Class II Major Histocompatibility Complex Transactivator (CIITA) and C-Type Lectin Domain Containing 16A (CLEC16A). Lamb IgA levels were also associated with the immunoglobulin heavy constant loci (IGH) complex, and adult IgE levels and lamb IgA and IgG levels were associated with the major histocompatibility complex (MHC). This study provides evidence of high heritability of a complex immunological trait under natural conditions and provides the first evidence from a genome-wide study that large effect genes located outside the MHC region exist for immune traits in the wild.

Sex-dependent dominance at a single locus maintains variation in age at maturity in salmon.

Males and females share many traits that have a common genetic basis; however, selection on these traits often differs between the sexes, leading to sexual conflict. Under such sexual antagonism, theory predicts the evolution of genetic architectures that resolve this sexual conflict. Yet, despite intense theoretical and empirical interest, the specific loci underlying sexually antagonistic phenotypes have rarely been identified, limiting our understanding of how sexual conflict impacts genome evolution and the maintenance of genetic diversity. Here we identify a large effect locus controlling age at maturity in Atlantic salmon (Salmo salar), an important fitness trait in which selection favours earlier maturation in males than females, and show it is a clear example of sex-dependent dominance that reduces intralocus sexual conflict and maintains adaptive variation in wild populations. Using high-density single nucleotide polymorphism data across 57 wild populations and whole genome re-sequencing, we find that the vestigial-like family member 3 gene (VGLL3) exhibits sex-dependent dominance in salmon, promoting earlier and later maturation in males and females, respectively. VGLL3, an adiposity regulator associated with size and age at maturity in humans, explained 39% of phenotypic variation, an unexpectedly large proportion for what is usually considered a highly polygenic trait. Such large effects are predicted under balancing selection from either sexually antagonistic or spatially varying selection. Our results provide the first empirical example of dominance reversal allowing greater optimization of phenotypes within each sex, contributing to the resolution of sexual conflict in a major and widespread evolutionary trade-off between age and size at maturity. They also provide key empirical evidence for how variation in reproductive strategies can be maintained over large geographical scales. We anticipate these findings will have a substantial impact on population management in a range of harvested species where trends towards earlier maturation have been observed.

Captive-bred Atlantic salmon released into the wild have fewer offspring than wild-bred fish and decrease population productivity.

The release of captive-bred animals into the wild is commonly practised to restore or supplement wild populations but comes with a suite of ecological and genetic consequences. Vast numbers of hatchery-reared fish are released annually, ostensibly to restore/enhance wild populations or provide greater angling returns. While previous studies have shown that captive-bred fish perform poorly in the wild relative to wild-bred conspecifics, few have measured individual lifetime reproductive success (LRS) and how this affects population productivity. Here, we analyse data on Atlantic salmon from an intensely studied catchment into which varying numbers of captive-bred fish have escaped/been released and potentially bred over several decades. Using a molecular pedigree, we demonstrate that, on average, the LRS of captive-bred individuals was only 36% that of wild-bred individuals. A significant LRS difference remained after excluding individuals that left no surviving offspring, some of which might have simply failed to spawn, consistent with transgenerational effects on offspring survival. The annual productivity of the mixed population (wild-bred plus captive-bred) was lower in years where captive-bred fish comprised a greater fraction of potential spawners. These results bolster previous empirical and theoretical findings that intentional stocking, or non-intentional escapees, threaten, rather than enhance, recipient natural populations.

Maternally derived anti-helminth antibodies predict offspring survival in a wild mammal.

The transfer of antibodies from mother to offspring provides crucial protection against infection to offspring during early life in humans and domestic and laboratory animals. However, few studies have tested the consequences of variation in maternal antibody transfer for offspring fitness in the wild. Further, separating the immunoprotective effects of antibodies from their association with nutritional resources provided by mothers is difficult. Here, we measured plasma levels of total and parasite-specific antibodies in neonatal (less than 10 days old) wild Soay sheep over 25 years to quantify variation in maternal antibody transfer and test its association with offspring survival. Maternal antibody transfer was predicted by maternal age and previous antibody responses, and was consistent within mothers across years. Neonatal total IgG antibody levels were positively related to early growth, suggesting they reflected nutritional transfer. Neonatal parasite-specific IgG levels positively predicted first-year survival, independent of lamb weight, total IgG levels and subsequent lamb parasite-specific antibody levels. This relationship was partly mediated via an indirect negative association with parasite burden. We show that among-female variation in maternal antibody transfer can have long-term effects on offspring growth, parasite burden and fitness in the wild, and is likely to impact naturally occurring host-parasite dynamics.

Life history trade-offs at a single locus maintain sexually selected genetic variation.

Sexual selection, through intra-male competition or female choice, is assumed to be a source of strong and sustained directional selection in the wild. In the presence of such strong directional selection, alleles enhancing a particular trait are predicted to become fixed within a population, leading to a decrease in the underlying genetic variation. However, there is often considerable genetic variation underlying sexually selected traits in wild populations, and consequently, this phenomenon has become a long-discussed issue in the field of evolutionary biology. In wild Soay sheep, large horns confer an advantage in strong intra-sexual competition, yet males show an inherited polymorphism for horn type and have substantial genetic variation in their horn size. Here we show that most genetic variation in this trait is maintained by a trade-off between natural and sexual selection at a single gene, relaxin-like receptor 2 (RXFP2). We found that an allele conferring larger horns, Ho(+), is associated with higher reproductive success, whereas a smaller horn allele, Ho(P), confers increased survival, resulting in a net effect of overdominance (that is, heterozygote advantage) for fitness at RXFP2. The nature of this trade-off is simple relative to commonly proposed explanations for the maintenance of sexually selected traits, such as genic capture ('good genes') and sexually antagonistic selection. Our results demonstrate that by identifying the genetic architecture of trait variation, we can determine the principal mechanisms maintaining genetic variation in traits under strong selection and explain apparently counter-evolutionary observations.

Natural Selection on Antihelminth Antibodies in a Wild Mammal Population.

An effective immune response is expected to confer fitness benefits through improved resistance to parasites but also incur energetic costs that negatively impact fitness-related traits, such as reproduction. The fitness costs and benefits of an immune response are likely to depend on host age, sex, and levels of parasite exposure. Few studies have examined the full extent to which patterns of natural selection on immune phenotypes vary across demographic groups and environments in the wild. Here, we assessed natural selection on plasma levels of three functionally distinct isotypes (IgA, IgE, and IgG) of antibodies against a prevalent nematode parasite measured in a wild Soay sheep population over 26 years. We found little support for environment-dependent selection or reproductive costs. However, antibody levels were negatively associated with parasite egg counts and positively associated with subsequent survival, albeit in a highly age- and isotype-dependent manner. Raised levels of antiparasite IgA best predicted reduced egg counts, but this did not predict survival in lambs. In adults increased antiparasite IgG predicted reduced egg counts, and in adult females IgG levels also positively predicted overwinter survival. Our results highlight the potential importance of age- and sex-dependent selection on immune phenotypes in nature and show that patterns of selection can vary even among functionally related immune markers.

The potential of shifting recombination hotspots to increase genetic gain in livestock breeding.

BACKGROUND: This study uses simulation to explore and quantify the potential effect of shifting recombination hotspots on genetic gain in livestock breeding programs. METHODS: We simulated three scenarios that differed in the locations of quantitative trait nucleotides (QTN) and recombination hotspots in the genome. In scenario 1, QTN were randomly distributed along the chromosomes and recombination was restricted to occur within specific genomic regions (i.e. recombination hotspots). In the other two scenarios, both QTN and recombination hotspots were located in specific regions, but differed in whether the QTN occurred outside of (scenario 2) or inside (scenario 3) recombination hotspots. We split each chromosome into 250, 500 or 1000 regions per chromosome of which 10% were recombination hotspots and/or contained QTN. The breeding program was run for 21 generations of selection, after which recombination hotspot regions were kept the same or were shifted to adjacent regions for a further 80 generations of selection. We evaluated the effect of shifting recombination hotspots on genetic gain, genetic variance and genic variance. RESULTS: Our results show that shifting recombination hotspots reduced the decline of genetic and genic variance by releasing standing allelic variation in the form of new allele combinations. This in turn resulted in larger increases in genetic gain. However, the benefit of shifting recombination hotspots for increased genetic gain was only observed when QTN were initially outside recombination hotspots. If QTN were initially inside recombination hotspots then shifting them decreased genetic gain. DISCUSSION: Shifting recombination hotspots to regions of the genome where recombination had not occurred for 21 generations of selection (i.e. recombination deserts) released more of the standing allelic variation available in each generation and thus increased genetic gain. However, whether and how much increase in genetic gain was achieved by shifting recombination hotspots depended on the distribution of QTN in the genome, the number of recombination hotspots and whether QTN were initially inside or outside recombination hotspots. CONCLUSIONS: Our findings show future scope for targeted modification of recombination hotspots e.g. through changes in zinc-finger motifs of the PRDM9 protein to increase genetic gain in production species.

Effect of manipulating recombination rates on response to selection in livestock breeding programs.

BACKGROUND: In this work, we performed simulations to explore the potential of manipulating recombination rates to increase response to selection in livestock breeding programs. METHODS: We carried out ten replicates of several scenarios that followed a common overall structure but differed in the average rate of recombination along the genome (expressed as the length of a chromosome in Morgan), the genetic architecture of the trait under selection, and the selection intensity under truncation selection (expressed as the proportion of males selected). Recombination rates were defined by simulating nine different chromosome lengths: 0.10, 0.25, 0.50, 1, 2, 5, 10, 15 and 20 Morgan, respectively. One Morgan was considered to be the typical chromosome length for current livestock species. The genetic architecture was defined by the number of quantitative trait variants (QTV) that affected the trait under selection. Either a large (10,000) or a small (1000 or 500) number of QTV was simulated. Finally, the proportions of males selected under truncation selection as sires for the next generation were equal to 1.2, 2.4, 5, or 10 %. RESULTS: Increasing recombination rate increased the overall response to selection and decreased the loss of genetic variance. The difference in cumulative response between low and high recombination rates increased over generations. At low recombination rates, cumulative response to selection tended to asymptote sooner and the genetic variance was completely eroded. If the trait under selection was affected by few QTV, differences between low and high recombination rates still existed, but the selection limit was reached at all rates of recombination. CONCLUSIONS: Higher recombination rates can enhance the efficiency of breeding programs to turn genetic variation into response to selection. However, to increase response to selection significantly, the recombination rate would need to be increased 10- or 20-fold. The biological feasibility and consequences of such large increases in recombination rates are unknown.

Low but significant genetic differentiation underlies biologically meaningful phenotypic divergence in a large Atlantic salmon population.

Despite decades of research assessing the genetic structure of natural populations, the biological meaning of low yet significant genetic divergence often remains unclear due to a lack of associated phenotypic and ecological information. At the same time, structured populations with low genetic divergence and overlapping boundaries can potentially provide excellent models to study adaptation and reproductive isolation in cases where high-resolution genetic markers and relevant phenotypic and life history information are available. Here, we combined single nucleotide polymorphism (SNP)-based population inference with extensive phenotypic and life history data to identify potential biological mechanisms driving fine-scale subpopulation differentiation in Atlantic salmon (Salmo salar) from the Teno River, a major salmon river in Europe. Two sympatrically occurring subpopulations had low but significant genetic differentiation (FST  = 0.018) and displayed marked differences in the distribution of life history strategies, including variation in juvenile growth rate, age at maturity and size within age classes. Large, late-maturing individuals were virtually absent from one of the two subpopulations, and there were significant differences in juvenile growth rates and size at age after oceanic migration between individuals in the respective subpopulations. Our findings suggest that different evolutionary processes affect each subpopulation and that hybridization and subsequent selection may maintain low genetic differentiation without hindering adaptive divergence.