Recurrent evolution of host-specialized races in a globally distributed parasite

The outcome of coevolutionary interactions is predicted to vary across landscapes depending on local conditions and levels of gene flow, with some populations evolving more extreme specializations than others. Using a globally distributed parasite of colonial seabirds, the tick Ixodes uriae, we examined how host availability and geographic isolation influences this process. In particular, we sampled ticks from 30 populations of six different seabird host species, three in the Southern Hemisphere and three in the Northern Hemisphere. We show that parasite races have evolved independently on hosts of both hemispheres. Moreover, the degree of differentiation between tick races varied spatially within each region and suggests that the divergence of tick races is an ongoing process that has occurred multiple times across isolated areas. As I. uriae is vector to the bacterium responsible for Lyme disease Borrelia burgdorferi sensu lato, these results may have important consequence for the epidemiology of this disease. With the increased occurrence of novel interspecific interactions due to global change, these results also stress the importance of the combined effects of gene flow and selection for parasite diversification.     

Coevolution and maladaptation

Many of the most commonly cited examples of exquisite adaptationare of coevolved symbioses. As we learn more about the coevolutionaryprocess, however, it is becoming increasingly evident that coevolutionmay also keep populations moderately maladapted much of thetime. As a result, coevolving populations may only rarely occupyadaptive peaks, because the selective landscape is under continualchange through reciprocal selection on the species themselves.These shifting patterns of coadaptation are further shaped bythe geographic structure of most species. Selection mosaicsacross landscapes and coevolutionary hotspots can favor differentevolutionary trajectories in different populations. The combinedaction of gene flow, random genetic drift, and local extinctionof populations may then continually remold these local patterns,creating a geographic mosaic in the degrees of maladaptationfound within local interactions. Recent mathematical modelsof the geographic mosaic of coevolution suggest that complexmosaics of maladaptation are a likely consequence of spatiallystructured species interactions. These models indicate thatthe spatial structure of maladaptation may depend upon the typeof coevolutionary interaction, the underlying selection mosaic,and patterns of gene flow across landscapes. By maintaininglocal polymorphisms and driving the divergence of populations,coevolution may produce spatial patterns of maladaptation thatare a source of ongoing innovation and diversification in speciesinteractions.     

Coevolutionary clines across selection mosaics

Abstract. Much of the dynamics of coevolution may be driven by the interplay between geographic variation in reciprocal selection (selection mosaics) and the homogenizing action of gene flow. We develop a genetic model of geographically structured coevolution in which gene flow links coevolving communities that may differ in both the direction and magnitude of reciprocal selection. The results show that geographically structured coevolution may lead to allele-frequency clines within both interacting species when fitnesses are spatially uniform or spatially heterogeneous. Furthermore, the results show that the behavior and shape of clines differ dramatically among different types of coevolutionary interaction. Antagonistic interactions produce dynamic clines that change shape rapidly through time, producing shifting patterns of local adaptation and maladaptation. Unlike antagonistic interactions, mutualisms generate stable equilibrium patterns that lead to fixed spatial patterns of adaptation. Interactions that vary between mutualism and antagonism produce both equilibrium and dynamic clines. Furthermore, the results demonstrate that these interactions may allow mutualisms to persist throughout the geographic range of an interaction, despite pockets of locally antagonistic selection. In all cases, the coevolved spatial patterns of allele frequencies are sensitive to the relative contributions of gene flow, selection, and overall habitat size, indicating that the appropriate scale for studies of geographically structured coevolution depends on the relative contributions of each of these factors.     

Toxicity and population structure of the Rough‐Skinned Newt (Taricha granulosa) outside the range of an arms race with resistant predators

Species interactions, and their fitness consequences, vary across the geographic range of a coevolutionary relationship. This spatial heterogeneity in reciprocal selection is predicted to generate a geographic mosaic of local adaptation, wherein coevolutionary traits are phenotypically variable from one location to the next. Under this framework, allopatric populations should lack variation in coevolutionary traits due to the absence of reciprocal selection. We examine phenotypic variation in tetrodotoxin (TTX) toxicity of the Rough‐Skinned Newt (Taricha granulosa) in regions of allopatry with its TTX‐resistant predator, the Common Garter Snake (Thamnophis sirtalis). In sympatry, geographic patterns of phenotypic exaggeration in toxicity and toxin‐resistance are closely correlated in prey and predator, implying that reciprocal selection drives phenotypic variation in coevolutionary traits. Therefore, in allopatry with TTX‐resistant predators, we expect to find uniformly low levels of newt toxicity. We characterized TTX toxicity in northwestern North America, including the Alaskan panhandle where Ta. granulosa occur in allopatry with Th. sirtalis. First, we used microsatellite markers to estimate population genetic structure and determine if any phenotypic variation in toxicity might be explained by historical divergence. We found northern populations of Ta. granulosa generally lacked population structure in a pattern consistent with northern range expansion after the Pleistocene. Next, we chose a cluster of sites in Alaska, which uniformly lacked genetic divergence, to test for phenotypic divergence in toxicity. As predicted, overall levels of newt toxicity were low; however, we also detected unexpected among‐ and within‐population variation in toxicity. Most notably, a small number of individuals contained large doses of TTX that rival means of toxic populations in sympatry with Th. sirtalis. Phenotypic variation in toxicity, despite limited neutral genetic divergence, suggests that factors other than reciprocal selection with Th. sirtalis likely contribute to geographic patterns of toxicity in Ta. granulosa.     

Hot Spots, Cold Spots, and the Geographic Mosaic Theory of Coevolution

Species interactions commonly coevolve as complex geographic mosaics of populations shaped by differences in local selection and gene flow. We use a haploid matching-alleles model for coevolution to evaluate how a pair of species coevolves when fitness interactions are reciprocal in some locations ("hot spots") but not in others ("cold spots"). Our analyses consider mutualistic and antagonistic interspecific interactions and a variety of gene flow patterns between hot and cold spots. We found that hot and cold spots together with gene flow influence coevolutionary dynamics in four important ways. First, hot spots need not be ubiquitous to have a global influence on evolution, although rare hot spots will not have a disproportionate impact unless selection is relatively strong there. Second, asymmetries in gene flow can influence local adaptation, sometimes creating stable equilibria at which species experience minimal fitness in hot spots and maximal fitness in cold spots, or vice versa. Third, asymmetries in gene flow are no more important than asymmetries in population regulation for determining the maintenance of local polymorphisms through coevolution. Fourth, intraspecific allele frequency differences among hot and cold spot populations evolve under some, but not all, conditions. That is, selection mosaics are indeed capable of producing spatially variable coevolutionary outcomes across the landscapes over which species interact. Altogether, our analyses indicate that coevolutionary trajectories can be strongly shaped by the geographic distribution of coevolutionary hot and cold spots, and by the pattern of gene flow among populations.     

Cryptic lineage divergence in marine environments: genetic differentiation at multiple spatial and temporal scales in the widespread intertidal goby Gobiosoma bosc

The adaptive radiation of the seven‐spined gobies (Gobiidae: Gobiosomatini) represents a classic example of how ecological specialization and larval retention can drive speciation through local adaptation. However, geographically widespread and phenotypically uniform species also do occur within Gobiosomatini. This lack of phenotypic variation across large geographic areas could be due to recent colonization, widespread gene flow, or stabilizing selection acting across environmental gradients. We use a phylogeographic approach to test these alternative hypotheses in the naked goby Gobiosoma bosc, a widespread and phenotypically invariable intertidal fish found along the Atlantic Coast of North America. Using DNA sequence from 218 individuals sampled at 15 localities, we document marked intraspecific genetic structure in mitochondrial and nuclear genes at three main geographic scales: (i) between Gulf of Mexico and Atlantic Coast, (ii) between the west coast of the Florida peninsula and adjacent Gulf of Mexico across the Apalachicola Bay, and (iii) at local scales of a few hundred kilometers. Clades on either side of Florida diverged about 8 million years ago, whereas some populations along the East Cost show divergent phylogroups that have differentiated within the last 200,000 years. The absence of noticeable phenotypic or ecological differentiation among lineages suggests the role of stabilizing selection on ancestral phenotypes, together with isolation in allopatry due to reduced dispersal and restricted gene flow, as the most likely explanation for their divergence. Haplotype phylogenies and spatial patterns of genetic diversity reveal frequent population bottlenecks followed by rapid population growth, particularly along the Gulf of Mexico. The magnitude of the genetic divergence among intraspecific lineages suggests the existence of cryptic species within Gobiosoma and indicates that modes of speciation can vary among lineages within Gobiidae.     

Gene flow and geographically structured coevolution

Many of the dynamic properties of coevolution may occur at the level of interacting populations, with local adaptation acting as a force of diversification, as migration between populations homogenizes these isolated interactions. This interplay between local adaptation and migration may be particularly important in structuring interactions that vary from mutualism to antagonism across the range of an interacting set of species, such as those between some plants and their insect herbivores, mammals and trypanosome parasites, and bacteria and plasmids that confer antibiotic resistance. Here we present a simple geographically structured genetic model of a coevolutionary interaction that varies between mutualism and antagonism among communities linked by migration. Inclusion of geographic structure with gene flow alters the outcomes of local interactions and allows the maintenance of allelic polymorphism across all communities under a range of selection intensities and rates of migration. Furthermore, inclusion of geographic structure with gene flow allows fixed mutualisms to be evolutionarily stable within both communities, even when selection on the interaction is antagonistic within one community. Moreover, the model demonstrates that the inclusion of geographic structure with gene flow may lead to considerable local maladaptation and trait mismatching as predicted by the geographic mosaic theory of coevolution.     

Parasite local adaptation in a geographic mosaic

A central prediction of the geographic mosaic theory of coevolution is that coevolving interspecific interactions will show varying degrees of local maladaptation. According to the theory, much of this local maladaptation is driven by selection mosaics and spatially intermingled coevolutionary hot and cold spots, rather than a simple balance between gene flow and selection. Here I develop a genetic model of host-parasite coevolution that is sufficiently general to incorporate selection mosaics, coevolutionary hot and cold spots, and a diverse array of genetic systems of infection/resistance. Results from this model show that the selection mosaics experienced by the interacting species are an important determinant of the sign and magnitude of local maladaptation. In some cases, this effect may be stronger than a previously described effect of relative rates of parasite and host gene flow. These results provide the first theoretical evidence that selection mosaics and coevolutionary hot and cold spots per se determine the magnitude and sign of local maladaptation. At the same time, however, these results demonstrate that coevolution in a geographic mosaic can lead to virtually any pattern of local adaptation or local maladaptation. Consequently, empirical studies that describe only patterns of local adaptation or maladaptation do not provide evidence either for or against the theory.     

Consequences of divergence and introgression for speciation in Andean cloud forest birds

Divergence with gene flow is well documented and reveals the influence of ecological adaptation on speciation. Yet, it remains intuitive that gene exchange inhibits speciation in many scenarios, particularly among ecologically similar populations. The influence of gene flow on the divergence of populations facing similar selection pressures has received less empirical attention than scenarios where differentiation is coupled with local environmental adaptation. I used a paired study design to test the influence of genomic divergence and introgression on plumage differentiation between ecologically similar allopatric replacements of Andean cloud forest birds. Through analyses of short‐read genome‐wide sequences from over 160 individuals in 16 codistributed lineages, I found that plumage divergence is associated with deep genetic divergence, implicating a prominent role of geographic isolation in speciation. By contrast, lineages that lack plumage divergence across the same geographic barrier are more recently isolated or exhibit a signature of secondary genetic introgression, indicating a negative relationship between gene flow and divergence in phenotypic traits important to speciation. My results suggest that the evolutionary outcomes of cycles of isolation and divergence in this important theatre of biotic diversification are sensitive to time spent in the absence of gene flow.     

Variation in infectivity and aggressiveness in space and time in wild host-pathogen systems: causes and consequences

Variation in host resistance and in the ability of pathogens to infect and grow (i.e. pathogenicity) is important as it provides the raw material for antagonistic (co)evolution and therefore underlies risks of disease spread, disease evolution and host shifts. Moreover, the distribution of this variation in space and time may inform us about the mode of coevolutionary selection (arms race vs. fluctuating selection dynamics) and the relative roles of G × G interactions, gene flow, selection and genetic drift in shaping coevolutionary processes. Although variation in host resistance has recently been reviewed, little is known about overall patterns in the frequency and scale of variation in pathogenicity, particularly in natural systems. Using 48 studies from 30 distinct host–pathogen systems, this review demonstrates that variation in pathogenicity is ubiquitous across multiple spatial and temporal scales. Quantitative analysis of a subset of extensively studied plant–pathogen systems shows that the magnitude of within‐population variation in pathogenicity is large relative to among‐population variation and that the distribution of pathogenicity partly mirrors the distribution of host resistance. At least part of the variation in pathogenicity found at a given spatial scale is adaptive, as evidenced by studies that have examined local adaptation at scales ranging from single hosts through metapopulations to entire continents and – to a lesser extent – by comparisons of pathogenicity with neutral genetic variation. Together, these results support coevolutionary selection through fluctuating selection dynamics. We end by outlining several promising directions for future research.     

Local adaptation in the Linum marginale-Melampsora lini host-pathogen interaction

The potential for local adaptation between pathogens and their hosts has generated strong theoretical and empirical interest with evidence both for and against local adaptation reported for a range of systems. We use the Linum marginale-Melampsora lini plant-pathogen system and a hierarchical spatial structure to investigate patterns of local adaptation within a metapopulation characterised by epidemic dynamics and frequent extinction of pathogen populations. Based on large sample sizes and comprehensive cross-inoculation trials, our analyses demonstrate strong local adaptation by Melampsora to its host populations, with this effect being greatest at regional scales, as predicted from the broader spatial scales at which M. lini disperses relative to L. marginale. However, there was no consistent trend for more distant pathogen populations to perform more poorly. Our results further show how the coevolutionary interaction between hosts and pathogens can be influenced by local structure such that resistant hosts select for generally virulent pathogens, while susceptible hosts select for more avirulent pathogens. Empirically, local adaptation has generally been tested in two contrasting ways: (1) pathogen performance on sympatric versus allopatric hosts; and (2) sympatric versus allopatric pathogens on a given host population. In situations where no host population is more resistant or susceptible than others when averaged across pathogen populations (and likewise, no pathogen population is more virulent or avirulent than others), results from these tests should generally be congruent. We argue that this is unlikely to be the case in the metapopulation situations that predominate in natural host-pathogen interactions, thus requiring tests that control simultaneously for variation in plant and pathogen populations.     

Evolution of population structure in an estuarine‐dependent marine fish

Restriction site‐associated DNA (RAD) sequencing was used to characterize neutral and adaptive genetic variation among geographic samples of red drum, Sciaenops ocellatus, an estuarine‐dependent fish found in coastal waters along the southeastern coast of the United States (Atlantic) and the northern Gulf of Mexico (Gulf). Analyses of neutral and outlier loci revealed three genetically distinct regional clusters: one in the Atlantic and two in the northern Gulf. Divergence in neutral loci indicated gradual genetic change and followed a linear pattern of isolation by distance. Divergence in outlier loci was at least an order of magnitude greater than divergence in neutral loci, and divergence between the regions in the Gulf was twice that of divergence between other regions. Discordance in patterns of genetic divergence between outlier and neutral loci is consistent with the hypothesis that the former reflects adaptive responses to environmental factors that vary on regional scales, while the latter largely reflects drift processes. Differences in basic habitat, initiated by glacial retreat and perpetuated by contemporary oceanic and atmospheric forces interacting with the geomorphology of the northern Gulf, followed by selection, appear to have led to reduced gene flow among red drum across the northern Gulf, reinforcing differences accrued during isolation and resulting in continued divergence across the genome. This same dynamic also may pertain to other coastal or nearshore fishes (18 species in 14 families) where genetically or morphologically defined sister taxa occur in the three regions.     

Metapopulation structure of a seed-predator weevil and its host plant in arms race coevolution

Although the importance of gene flow in the geographic structuring of host-parasite interactions has been well discussed, little is known about how dispersal drives the spatial dynamics of other types of coevolutionary interactions in nature. We evaluated the roles of gene flow in the geographically structured processes of a predator-prey arms race involving a seed-predatory weevil with a long mouthpart and its host camellia plant with a thick fruit coat. Molecular genetic analyses showed that both weevil and camellia populations were structured at a spatial scale of several kilometers. Importantly, the spatial pattern of the migration of weevils, but not that of camellias, imposed significant effects on the geographic configuration of the levels of coevolutionary escalation. This result suggests that even if migration is limited in one species (camellia), local coevolution with the other species that migrates between neighboring localities (weevil) can reduce the interpopulation difference in the local adaptive optima of the former species. Thus, gene flow of a species potentially homogenizes the local biological environments provided by the species and thereby promotes the evolutionary convergence of its coevolving counterparts. Consequently, by focusing on coevolutionary interactions in natural communities, "indirect" effects of gene flow on the adaptive divergence of organisms could be identified.     

Gut microbial diversity across a contact zone for California voles: Implications for lineage divergence of hosts and mitonuclear mismatch in the assembly of the mammalian gut microbiome

Gut microbial diversity is thought to reflect the co‐evolution of microbes and their hosts as well as current host‐specific attributes such as genetic background and environmental setting. To explore interactions among these parameters, we characterized variation in gut microbiome composition of California voles (Microtus californicus) across a contact zone between two recently diverged lineages of this species. Because this contact zone contains individuals with mismatched mitochondrial‐nuclear genomes (cybrids), it provides an important opportunity to explore how different components of the genotype contribute to gut microbial diversity. Analyses of bacterial 16S rRNA sequences and joint species distribution modelling revealed that host genotypes and genetic differentiation among host populations together explained more than 50% of microbial community variation across our sampling transect. The ranked importance (most to least) of factors contributing to gut microbial diversity in our study populations were: genome‐wide population differentiation, local environmental conditions, and host genotypes. However, differences in microbial communities among vole populations (β‐diversity) did not follow patterns of lineage divergence (i.e., phylosymbiosis). Instead, among‐population variation was best explained by the spatial distribution of hosts, as expected if the environment is a primary source of gut microbial diversity (i.e., dispersal limitation hypothesis). Across the contact zone, several bacterial taxa differed in relative abundance between the two parental lineages as well as among individuals with mismatched mitochondrial and nuclear genomes. Thus, genetic divergence among host lineages and mitonuclear genomic mismatches may also contribute to microbial diversity by altering interactions between host genomes and gut microbiota (i.e., hologenome speciation hypothesis).     

Current geography masks dynamic history of gene flow during speciation in northern Australian birds

Genome divergence is greatly influenced by gene flow during early stages of speciation. As populations differentiate, geographic barriers can constrain gene flow and so affect the dynamics of divergence and speciation. Current geography, specifically disjunction and continuity of ranges, is often used to predict the historical gene flow during the divergence process. We test this prediction in eight meliphagoid bird species complexes codistributed in four regions. These regions are separated by known biogeographical barriers across northern Australia and Papua New Guinea. We find that bird populations currently separated by terrestrial habitat barriers within Australia and marine barriers between Australia and Papua New Guinea have a range of divergence levels and probability of gene flow not associated with current range connectivity. Instead, geographic distance and historical range connectivity better predict divergence and probability of gene flow. In this dynamic environmental context, we also find support for a nonlinear decrease of the probability of gene flow during the divergence process. The probability of gene flow initially decreases gradually after a certain level of divergence is reached. Its decrease then accelerates until the probability is close to zero. This implies that although geographic connectivity may have more of an effect early in speciation, other factors associated with higher divergence may play a more important role in influencing gene flow midway through and later in speciation. Current geographic connectivity may then mislead inferences regarding potential for gene flow during speciation under a complex and dynamic history of geographic and reproductive isolation.     

Geographic structure in a widespread plant-mycorrhizal interaction: pines and false truffles

Mutualistic interactions are likely to exhibit a strong geographic mosaic in their coevolutionary dynamics, but the structure of geographic variation in these interactions is much more poorly characterized than in host-parasite interactions. We used a cross-inoculation experiment to characterize the scales and patterns at which geographic structure has evolved in an interaction between three pine species and one ectomycorrhizal fungus species along the west coast of North America. We found substantial and contrasting patterns of geographic interaction structure for the plants and fungi. The fungi exhibited a clinal pattern of local adaptation to their host plants across the geographic range of three coastal pines. In contrast, plant growth parameters were unaffected by fungal variation, but varied among plant populations and species. Both plant and fungal performance measures varied strongly with latitude. This set of results indicates that in such widespread species interactions, interacting species may evolve asymmetrically in a geographic mosaic because of differing evolutionary responses to clinally varying biotic and abiotic factors.     

PERCHED AT THE MITO‐NUCLEAR CROSSROADS: DIVERGENT MITOCHONDRIAL LINEAGES CORRELATE WITH ENVIRONMENT IN THE FACE OF ONGOING NUCLEAR GENE FLOW IN AN AUSTRALIAN BIRD

Relationships among multilocus genetic variation, geography, and environment can reveal how evolutionary processes affect genomes. We examined the evolution of an Australian bird, the eastern yellow robin Eopsaltria australis, using mitochondrial (mtDNA) and nuclear (nDNA) genetic markers, and bioclimatic variables. In southeastern Australia, two divergent mtDNA lineages occur east and west of the Great Dividing Range, perpendicular to latitudinal nDNA structure. We evaluated alternative scenarios to explain this striking discordance in landscape genetic patterning. Stochastic mtDNA lineage sorting can be rejected because the mtDNA lineages are essentially distinct geographically for > 1500 km. Vicariance is unlikely: the Great Dividing Range is neither a current barrier nor was it at the Last Glacial Maximum according to species distribution modeling; nuclear gene flow inferred from coalescent analysis affirms this. Female philopatry contradicts known female‐biased dispersal. Contrasting mtDNA and nDNA demographies indicate their evolutionary histories are decoupled. Distance‐based redundancy analysis, in which environmental temperatures explain mtDNA variance above that explained by geographic position and isolation‐by‐distance, favors a nonneutral explanation for mitochondrial phylogeographic patterning. Thus, observed mito‐nuclear discordance accords with environmental selection on a female‐linked trait, such as mtDNA, mtDNA–nDNA interactions or genes on W‐chromosome, driving mitochondrial divergence in the presence of nuclear gene flow.     

Drift versus selection as drivers of phenotypic divergence at small spatial scales: The case of Belgjarskógur threespine stickleback

Divergence in phenotypic traits is facilitated by a combination of natural selection, phenotypic plasticity, gene flow, and genetic drift, whereby the role of drift is expected to be particularly important in small and isolated populations. Separating the components of phenotypic divergence is notoriously difficult, particularly for multivariate phenotypes. Here, we assessed phenotypic divergence of threespine stickleback (Gasterosteus aculeatus) across 19 semi‐interconnected ponds within a small geographic region (~7.5 km²) using comparisons of multivariate phenotypic divergence (PST), neutral genetic (FST), and environmental (EST) variation. We found phenotypic divergence across the ponds in a suite of functionally relevant phenotypic traits, including feeding, defense, and swimming traits, and body shape (geometric morphometric). Comparisons of PSTs with FSTs suggest that phenotypic divergence is predominantly driven by neutral processes or stabilizing selection, whereas phenotypic divergence in defensive traits is in accordance with divergent selection. Comparisons of population pairwise PSTs with ESTs suggest that phenotypic divergence in swimming traits is correlated with prey availability, whereas there were no clear associations between phenotypic divergence and environmental difference in the other phenotypic groups. Overall, our results suggest that phenotypic divergence of these small populations at small geographic scales is largely driven by neutral processes (gene flow, drift), although environmental determinants (natural selection or phenotypic plasticity) may play a role.     

Genomic divergence across ecological gradients in the Central African rainforest songbird (Andropadus virens)

The little greenbul, a common rainforest passerine from sub‐Saharan Africa, has been the subject of long‐term evolutionary studies to understand the mechanisms leading to rainforest speciation. Previous research found morphological and behavioural divergence across rainforest–savannah transition zones (ecotones), and a pattern of divergence with gene flow suggesting divergent natural selection has contributed to adaptive divergence and ecotones could be important areas for rainforests speciation. Recent advances in genomics and environmental modelling make it possible to examine patterns of genetic divergence in a more comprehensive fashion. To assess the extent to which natural selection may drive patterns of differentiation, here we investigate patterns of genomic differentiation among populations across environmental gradients and regions. We find compelling evidence that individuals form discrete genetic clusters corresponding to distinctive environmental characteristics and habitat types. Pairwise F_ST between populations in different habitats is significantly higher than within habitats, and this differentiation is greater than what is expected from geographic distance alone. Moreover, we identified 140 SNPs that showed extreme differentiation among populations through a genomewide selection scan. These outliers were significantly enriched in exonic and coding regions, suggesting their functional importance. Environmental association analysis of SNP variation indicates that several environmental variables, including temperature and elevation, play important roles in driving the pattern of genomic diversification. Results lend important new genomic evidence for environmental gradients being important in population differentiation.     

Genomic Adaptations to Salinity Resist Gene Flow in the Evolution of Floridian Watersnakes

The migration-selection balance often governs the evolution of lineages, and speciation with gene flow is now considered common across the tree of life. Ecological speciation is a process that can facilitate divergence despite gene flow due to strong selective pressures caused by ecological differences; however, the exact traits under selection are often unknown. The transition from freshwater to saltwater habitats provides strong selection targeting traits with osmoregulatory function. Several lineages of North American watersnakes (Nerodia spp.) are known to occur in saltwater habitat and represent a useful system for studying speciation by providing an opportunity to investigate gene flow and evaluate how species boundaries are maintained or degraded. We use double digest restriction-site associated DNA sequencing to characterize the migration-selection balance and test for evidence of ecological divergence within the Nerodia fasciata-clarkii complex in Florida. We find evidence of high intraspecific gene flow with a pattern of isolation-by-distance underlying subspecific lineages. However, we identify genetic structure indicative of reduced gene flow between inland and coastal lineages suggesting divergence due to isolation-by-environment. This pattern is consistent with observed environmental differences where the amount of admixture decreases with increased salinity. Furthermore, we identify significantly enriched terms related to osmoregulatory function among a set of candidate loci, including several genes that have been previously implicated in adaptation to salinity stress. Collectively, our results demonstrate that ecological differences, likely driven by salinity, cause strong divergent selection which promotes divergence in the N. fasciata-clarkii complex despite significant gene flow.