NEUTRAL GENE FLOW ACROSS SINGLE LOCUS CLINES

Persistence of a Geographically-Stable Hybrid Zone in Puerto Rican Dwarf Geckos

Determining the mechanisms that create and maintain biodiversity is a central question in ecology and evolution. Speciation is the process that creates biodiversity. Speciation is mediated by incompatibilities that lead to reproductive isolation between divergent populations and these incompatibilities can be observed in hybrid zones. Gecko lizards are a speciose clade possessing an impressive diversity of behavioral and morphological traits. In geckos, however, our understanding of the speciation process is negligible. To address this gap, we used genetic sequence data (both mitochondrial and nuclear markers) to revisit a putative hybrid zone between Sphaerodactylus nicholsi and Sphaerodactylus townsendi in Puerto Rico, initially described in 1984. First, we addressed discrepancies in the literature on the validity of both species. Second, we sampled a 10-km-wide transect across the putative hybrid zone and tested explicit predictions about its dynamics using cline models. Third, we investigated potential causes for the hybrid zone using species distribution modeling and simulations; namely, whether unique climatic variables within the hybrid zone might elicit selection for intermediate phenotypes. We find strong support for the species-level status of each species and no evidence of movement, or unique climatic variables near the hybrid zone. We suggest that this narrow hybrid zone is geographically stable and is maintained by a combination of dispersal and selection. Thus, this work has identified an extant model system within geckos that that can be used for future investigations detailing genetic mechanisms of reproductive isolation in an understudied vertebrate group.

Genetic Incompatibilities Between Mitochondria and Nuclear Genes: Effect on Gene Flow and Speciation

The process of speciation is, according to the biological species concept, the reduction in gene flow between genetically diverging populations. Most of the previous theoretical studies analyzed the effect of nuclear genetic incompatibilities on gene flow. There is, however, an increasing number of empirical examples suggesting that cytoplasmically inherited genetic elements play an important role in speciation. Here, we present a theoretical analysis of mitochondrial driven speciation, in which genetic incompatibilities occur between mitochondrial haplotypes and nuclear alleles. Four population genetic models with mainland-island structure were analyzed that differ with respect to the type of incompatibility and the underlying genetics. Gene flow reduction was measured on selectively neutral alleles of an unlinked locus and quantified by the effective migration rate. Analytical formulae for the different scenarios were derived using the fitness graph method. For the models with haploid genetics, we found that mito-nuclear incompatibilities (MtNI) are as strong as nuclear-nuclear incompatibilities (NNI) in reducing gene flow at the unlinked locus, but only if males and females migrate in equal number. For models with diploid genetics, we found that MtNI reduce gene flow stronger than NNI when incompatibilities are recessive, but weaker when they are dominant. For both haploid and diploid MtNI, we found that gene flow reduction is stronger if females are the migrating sex, but weaker than NNI when males are the migrating sex. These results encourage further examination on the role of mitochondria on genetic divergence and speciation and point toward specific factors (e.g., migrating sex) that could be the focus of an empirical test.

Varied Genomic Responses to Maladaptive Gene Flow and Their Evidence

Adaptation to a local environment often occurs in the face of maladaptive gene flow. In this perspective, I discuss several ideas on how a genome may respond to maladaptive gene flow during adaptation. On the one hand, selection can build clusters of locally adaptive alleles at fortuitously co-localized loci within a genome, thereby facilitating local adaptation with gene flow ('allele-only clustering'). On the other hand, the selective pressure to link adaptive alleles may drive co-localization of the actual loci relevant for local adaptation within a genome through structural genome changes or an evolving intra-genomic crossover rate ('locus clustering'). While the expected outcome is, in both cases, a higher frequency of locally adaptive alleles in some genome regions than others, the molecular units evolving in response to gene flow differ (i.e., alleles versus loci). I argue that, although making this distinction is important, we commonly lack the critical empirical evidence to do so. This is mainly because many current approaches are biased towards detecting local adaptation in genome regions with low crossover rates. The importance of low-crossover genome regions for adaptation with gene flow, such as in co-localizing relevant loci within a genome, thus remains unclear. Future empirical investigations should address these questions by making use of comparative genomics, where multiple de novo genome assemblies from species evolved under different degrees of genetic exchange are compared. This research promises to advance our understanding of how a genome adapts to maladaptive gene flow, thereby promoting adaptive divergence and reproductive isolation.

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

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

The Strength of Selection against Neanderthal Introgression

Hybridization between humans and Neanderthals has resulted in a low level of Neanderthal ancestry scattered across the genomes of many modern-day humans. After hybridization, on average, selection appears to have removed Neanderthal alleles from the human population. Quantifying the strength and causes of this selection against Neanderthal ancestry is key to understanding our relationship to Neanderthals and, more broadly, how populations remain distinct after secondary contact. Here, we develop a novel method for estimating the genome-wide average strength of selection and the density of selected sites using estimates of Neanderthal allele frequency along the genomes of modern-day humans. We confirm that East Asians had somewhat higher initial levels of Neanderthal ancestry than Europeans even after accounting for selection. We find that the bulk of purifying selection against Neanderthal ancestry is best understood as acting on many weakly deleterious alleles. We propose that the majority of these alleles were effectively neutral—and segregating at high frequency—in Neanderthals, but became selected against after entering human populations of much larger effective size. While individually of small effect, these alleles potentially imposed a heavy genetic load on the early-generation human–Neanderthal hybrids. This work suggests that differences in effective population size may play a far more important role in shaping levels of introgression than previously thought.

A small percentage of Neanderthal DNA is present in the genomes of many contemporary human populations due to hybridization tens of thousands of years ago. Much of this Neanderthal DNA appears to be deleterious in humans, and natural selection is acting to remove it. One hypothesis is that the underlying alleles were not deleterious in Neanderthals, but rather represent genetic incompatibilities that became deleterious only once they were introduced to the human population. If so, reproductive barriers must have evolved rapidly between Neanderthals and humans after their split. Here, we show that observed patterns of Neanderthal ancestry in modern humans can be explained simply as a consequence of the difference in effective population size between Neanderthals and humans. Specifically, we find that on average, selection against individual Neanderthal alleles is very weak. This is consistent with the idea that Neanderthals over time accumulated many weakly deleterious alleles that in their small population were effectively neutral. However, after introgressing into larger human populations, those alleles became exposed to purifying selection. Thus, rather than being the result of hybrid incompatibilities, differences between human and Neanderthal effective population sizes appear to have played a key role in shaping our present-day shared ancestry.

The genomic signature of parallel adaptation from shared genetic variation

Parallel adaptation is common and may often occur from shared genetic variation, but the genomic consequences of this process remain poorly understood. We first use individual-based simulations to demonstrate that comparisons between populations adapted in parallel to similar environments from shared variation reveal a characteristic genomic signature around a selected locus: a low-divergence valley centred at the locus and flanked by twin peaks of high divergence. This signature is initiated by the hitchhiking of haplotype tracts differing between derived populations in the broader neighbourhood of the selected locus (driving the high-divergence twin peaks) and shared haplotype tracts in the tight neighbourhood of the locus (driving the low-divergence valley). This initial hitchhiking signature is reinforced over time because the selected locus acts as a barrier to gene flow from the source to the derived populations, thus promoting divergence by drift in its close neighbourhood. We next empirically confirm the peak-valley-peak signature by combining targeted and RAD sequence data at three candidate adaptation genes in multiple marine (source) and freshwater (derived) populations of threespine stickleback. Finally, we use a genome-wide screen for the peak-valley-peak signature to discover additional genome regions involved in parallel marine-freshwater divergence. Our findings offer a new explanation for heterogeneous genomic divergence and thus challenge the standard view that peaks in population divergence harbour divergently selected loci and that low-divergence regions result from balancing selection or localized introgression. We anticipate that genome scans for peak-valley-peak divergence signatures will promote the discovery of adaptation genes in other organisms.

Biased sex-ratio and sex-biased heterozygote disadvantage affect the maintenance of a genetic polymorphism and the properties of hybrid zones

The evolution of biodiversity is a major issue of modern biology, and it is becoming increasingly topical as the ongoing erosion of diversity puts serious threats on human well-being. An elementary mechanism that allows maintaining diversity is the interplay between dispersal and heterozygote selective disadvantage, which can lead to self-sustainable spatial genetic structures and is central to the stability of hybrid zones. Theoretical studies supporting the importance of this mechanism assume a balanced sex-ratio and a heterozygote disadvantage equally affecting both sexes, despite the multiplicity of empirical evidence that (i) adult sex-ratio is usually biased towards either male or female and that (ii) heterozygote disadvantage often affects a single sex. We expanded the existing theory by weighting the strength of selection against heterozygote according to the biased in sex-ratio and in heterozygote disadvantage. The range of conditions allowing for the maintenance of polymorphism can then either double or vanish. We discuss the implications of such finding for birds, mammals and insects diversity. Finally, we provide simple analytical predictions about the effect of those biased on the width and speed of hybrid zones and on the time for the spread of beneficial mutations through such zones.

Genomics of isolation in hybrids

Hybrid zones are common in nature and can offer critical insights into the dynamics and components of reproductive isolation. Hybrids between diverged lineages are particularly informative about the genetic architecture of reproductive isolation, because introgression in an admixed population is a direct measure of isolation. In this paper, we combine simulations and a new statistical model to determine the extent to which different genetic architectures of isolation leave different signatures on genome-level patterns of introgression. We found that reproductive isolation caused by one or several loci of large effect caused greater heterogeneity in patterns of introgression than architectures involving many loci with small fitness effects, particularly when isolating factors were closely linked. The same conditions that led to heterogeneous introgression often resulted in a reasonable correspondence between outlier loci and the genetic loci that contributed to isolation. However, demographic conditions affected both of these results, highlighting potential limitations to the study of the speciation genomics. Further progress in understanding the genomics of speciation will require large-scale empirical studies of introgression in hybrid zones and model-based analyses, as well as more comprehensive modelling of the expected levels of isolation with different demographies and genetic architectures of isolation.

Establishment of new mutations under divergence and genome hitchhiking

Theoretical models addressing genome-wide patterns of divergence during speciation are needed to help us understand the evolutionary processes generating empirical patterns. Here, we examine a critical issue concerning speciation-with-gene flow: to what degree does physical linkage (r < 0.5) of new mutations to already diverged genes aid the build-up of genomic islands of differentiation? We used simulation and analytical approaches to partition the probability of establishment for a new divergently selected mutation when the mutation (i) is the first to arise in an undifferentiated genome (the direct effect of selection), (ii) arises unlinked to any selected loci (r = 0.5), but within a genome that has some already diverged genes (the effect of genome-wide reductions in gene flow for facilitating divergence, which we term 'genome hitchhiking'), and (iii) arises in physical linkage to a diverged locus (divergence hitchhiking). We find that the strength of selection acting directly on a new mutation is generally the most important predictor for establishment, with divergence and genomic hitchhiking having smaller effects. We outline the specific conditions under which divergence and genome hitchhiking can aid mutation establishment. The results generate predictions about genome divergence at different points in the speciation process and avenues for further work.

The efficacy of divergence hitchhiking in generating genomic islands during ecological speciation

Genes under divergent selection flow less readily between populations than other loci. This observation has led to verbal "divergence hitchhiking" models of speciation in which decreased interpopulation gene flow surrounding loci under divergent selection can generate large regions of differentiation within the genome (genomic islands). The efficacy of this model in promoting speciation depends on the size of the region affected by divergence hitchhiking. Empirical evidence is mixed, with examples of both large and small genomic islands. To address these empirical discrepancies and to formalize the theory, we present mathematical models of divergence hitchhiking, which examine neutral differentiation around selected sites. For a single locus under selection, regions of differentiation do not extend far along a chromosome away from a selected site unless both effective population sizes and migration rates are low. When multiple loci are considered, regions of differentiation can be larger. However, with many loci under selection, genome-wide divergence occurs and genomic islands are erased. The results show that divergence hitchhiking can generate large regions of differentiation, but that the conditions under which this occurs are limited. Thus, speciation may often require multifarious selection acting on many, isolated and physically unlinked genes. How hitchhiking promotes further adaptive divergence warrants consideration.

Patterns of introgression across an expanding hybrid zone: analysing historical patterns of gene flow using nonequilibrium approaches

Previous studies suggest that the hybrid zone between two taxa in the Piriqueta caroliniana complex in central Florida (south-eastern North America) has recently expanded with hybrids replacing parental genotypes across a broad region of the Florida peninsula. Here I use patterns of genetic disequilibria and levels of differentiation among populations to infer historical patterns of introgression across this broad hybrid zone. There were strong positive associations among taxon-specific alleles at the southern extent of hybridization, but disequilibria values were close to zero across the central and northern portions of the hybrid zone. Levels of among-population differentiation in the central portion of the hybrid zone were relatively low, and increased towards the northern extent of hybridization. The high levels of disequilibria are coincident with the sharp clines at the southern end of the hybrid zone, suggesting that there is a tension zone in this region that is maintained by selection against hybrids and dispersal from parental regions. The levels of disequilibria within populations and patterns of differentiation among populations are consistent with historical introgression and northward expansion of this hybrid zone, which may have slowed or ceased in recent generations.

Clines maintained by overdominant selection in hybrid zones

Hybrid zone models often consider environment-independent selection to operate against all hybrids. However, empirical studies suggest that hybrids may be as fit or fitter than the hybridizing parental taxa in some environments. In this study we develop a novel mathematical model to explore the effects of one form of hybrid superiority on the genetic structure of hybrid zones. Our primary goals were to investigate the allele frequency clines at a locus experiencing overdominant selection and at a linked neutral or underdominant locus. Our results indicate that overdominant selection results in flat equilibrium allele frequency clines throughout the hybrid zone and an excess of heterozygosity relative to Hardy-Weinberg equilibrium. However, the genetic clines at linked neutral or underdominant loci tend not to reflect this overdominance even when the loci are tightly linked. Overall, we conclude that overdominance is unlikely to be detected in genetic surveys unless many loci are assayed.

A comparison of multilocus clines maintained by environmental adaptation or by selection against hybrids

There has recently been considerable debate over the relative importance of selection against hybrids ("endogenous" selection) vs. adaptation to different environments ("exogenous") in maintaining stable hybrid zones and hence in speciation. Single-locus models of endogenous and exogenous viability selection generate clines of similar shape, but the comparison has not been extended to multilocus systems, which are both quantitatively and qualitatively very different from the single-locus case. Here we develop an analytical multilocus model of differential adaptation across an environmental transition and compare it to previous heterozygote disadvantage models. We show that the shape of clines generated by exogenous selection is indistinguishable from that generated by endogenous selection. A stochastic simulation model is used to test the robustness of the analytical description to the effects of drift and strong selection, and confirms the prediction that pairwise linkage disequilibria are predominantly generated by migration. However, although analytical predictions for the width of clines maintained by heterozygote disadvantage fit well with the simulation results, those for environmental adaptation are consistently too narrow; reasons for the discrepancy are discussed. There is a smooth transition between a system in which a set of loci effectively act independently of each other and one in which they act as a single nonrecombining unit.

Global dispersal reduces local diversity

Metapopulation models and stepping-stone models in genetics are based on very different underlying dispersal structures, yet it can be difficult to distinguish the behaviour of the two kinds of models. We demonstrate a striking qualitative difference in the equilibrium behaviour possible with these two kinds of dispersal. If, in a local patch, there are multiple stable equilibria (and consequently an unstable equilibrium), we demonstrate that, for the spatial system with a metapopulation structure, at equilibrium every patch has to be near one of the stable equilibria. This contrasts with the clinal structure possible with a stepping-stone or continuous space model; thus the result can be used to deduce qualitative information about the form of dispersal from observations of allele frequencies.