ADAPTIVE POPULATION DIVERGENCE IN CRYPTIC COLOR-PATTERN FOLLOWING A REDUCTION IN GENE FLOW

Adaptive population divergence is often driven by divergent natural selection, but can be constrained by the homogenizing effect of gene flow between populations. Indeed, a common pattern in nature is an inverse correlation between the degree of adaptive phenotypic divergence between populations and levels of gene flow between populations. However, there is essentially no experimental data on whether this correlation arises because gene flow constrains adaptation or, conversely, because adaptive divergence causes barriers to gene flow (ecological speciation). Here, I report increased adaptive divergence in cryptic color pattern between a pair of Timema insect populations following an experimental reduction in between-population gene flow. The reduction in gene flow arose due to a natural experiment, and thus was not replicated at a second site. However, temporal replication of the trends among six generations of data, coupled with a lack of increased adaptive divergence for two other population pairs where gene flow was not manipulated (i.e., control sites), argues that the results did not arise by chance. Estimates of dispersal ability and population size further support reduced gene flow, rather than increased genetic drift, as the cause of divergence. Thus, the findings provide experimental evidence that gene flow constrains adaptation in nature.     

PHENOLOGICAL ISOLATION,GENE FLOW AND DEVELOPMENTAL DIFFERENCES AMONG LOW- AND HIGH-ELEVATION POPULATIONS OF EUPHILOTES ENOPTES (LEPIDOPTERA: LYCAENIDAE)

Populations of the specialist herbivore, Euphilotes enoptes (Lepidoptera: Lycaenidae), along three elevational transects in the mountains of central Washington state, differed markedly in the phenology of adult flight. In spite of this apparent limitation to gene flow, six allozyme loci revealed substantial gene exchange among populations along these gradients. The elevational difference, and thus the phenological difference, between populations has not influenced the extent of gene flow between them. Because the direct exchange of genes between low- and high-elevation populations is very unlikely, gene flow between them has probably occurred in a stepwise fashion via intermediate populations. It is hypothesized that such gene flow has been biased in an uphill direction due to the combined effects of source size and oviposition behavior. Adult emergence times of populations in the same region are positively correlated with elevation in a nonlinear fashion, consistent with the hypothesis that gene flow from low-elevation populations has been swamping selection at higher altitudes.     

Ecological speciation along an elevational gradient in a tropical passerine bird?

Local adaptation of populations along elevational gradients is well known, but conclusive evidence that such divergence has resulted in the origin of distinct species in parapatry remains lacking. We integrated morphological, vocal, genetic and behavioural data to test predictions pertaining to the hypothesis of parapatric ecological speciation associated with elevation in populations of a tropical montane songbird, the Grey‐breasted Wood‐wren (Henicorhina leucophrys: Troglodytidae), from the Sierra Nevada de Santa Marta, Colombia. We confirmed that two distinct populations exist along the elevational gradient. Phylogenetic analyses tentatively indicate that the two populations are not sister taxa, suggesting they did not differentiate from a single ancestor along the gradient, but rather resulted from separate colonization events. The populations showed marked divergence in morphometrics, vocalizations and genetic variation in mitochondrial and nuclear loci, and little to no evidence of hybridization. Individuals of both populations responded more strongly to their own local songs than to songs from another elevation. Although the two forms do not appear to have differentiated locally in parapatry, morphological and vocal divergence along the elevational gradient is consistent with adaptation, suggesting a possible link between adaptive evolution in morphology and songs and the origin of reproductive isolation via a behavioural barrier to gene flow. The adaptive value of phenotypic differences between populations requires additional study.     

Adaptive population divergence and directional gene flow across steep elevational gradients in a climate‐sensitive mammal

The ecological effects of climate change have been shown in most major taxonomic groups; however, the evolutionary consequences are less well‐documented. Adaptation to new climatic conditions offers a potential long‐term mechanism for species to maintain viability in rapidly changing environments, but mammalian examples remain scarce. The American pika (Ochotona princeps) has been impacted by recent climate‐associated extirpations and range‐wide reductions in population sizes, establishing it as a sentinel mammalian species for climate change. To investigate evidence for local adaptation and reconstruct patterns of genomic diversity and gene flow across rapidly changing environments, we used a space‐for‐time design and restriction site‐associated DNA sequencing to genotype American pikas along two steep elevational gradients at 30,966 SNPs and employed independent outlier detection methods that scanned for genotype‐environment associations. We identified 338 outlier SNPs detected by two separate analyses and/or replicated in both transects, several of which were annotated to genes involved in metabolic function and oxygen transport. Additionally, we found evidence of directional gene flow primarily downslope from high‐elevation populations, along with reduced gene flow at outlier loci. If this trend continues, elevational range contractions in American pikas will likely be from local extirpation rather than upward movement of low‐elevation individuals; this, in turn, could limit the potential for adaptation within this landscape. These findings are of particular relevance for future conservation and management of American pikas and other elevationally restricted, thermally sensitive species.     

No effect of elevation and fragmentation on genetic diversity and structure in Polylepis australis trees from central Argentina

Phenological differences in flowering arising along elevational gradients may be caused by either local adaptation or phenotypic plasticity. Local adaptation can lead to reproductive isolation of populations at different elevational zones and thus produce elevational genetic structuring, while phenotypic plasticity does not produce elevational genetic structuring. In this study, we examined the effects of elevation and fragmentation on genetic diversity and structure of Polylepis australis populations, where individuals exhibit phenological differences in flowering along an elevational gradient. We assessed the polymorphism of amplified fragment length polymorphism markers in adults and saplings from one conserved and one fragmented forest covering elevations from 1600 to 2600 m asl. Over 98% of variation was found within populations, and we found very low and similar genetic differentiation along elevational gradients for adults and saplings in both continuous and fragmented forests. In addition, there was no significant relationship between genetic diversity and elevation. Results indicated that phenological differences along elevational gradients are more likely caused by phenotypic plasticity than local adaptation, and fragmentation does not appear to have affected genetic diversity and differentiation in the studied populations. Results therefore imply that if necessary, seeds for reforestation purposes may be collected from different elevations to the seeding or planting sites.     

Modifying effects of phenotypic plasticity on interactions among natural selection, adaptation and gene flow

Divergent natural selection, adaptive divergence and gene flow may interact in a number of ways. Recent studies have focused on the balance between selection and gene flow in natural populations, and empirical work has shown that gene flow can constrain adaptive divergence, and that divergent selection can constrain gene flow. A caveat is that phenotypic diversification may be under the direct influence of environmental factors (i.e. it may be due to phenotypic plasticity), in addition to partial genetic influence. In this case, phenotypic divergence may occur between populations despite high gene flow that imposes a constraint on genetic divergence. Plasticity may dampen the effects of natural selection by allowing individuals to rapidly adapt phenotypically to new conditions, thus slowing adaptive genetic divergence. On the other hand, plasticity may promote future adaptive divergence by allowing populations to persist in novel environments. Plasticity may promote gene flow between selective regimes by allowing dispersers to adapt to alternate conditions, or high gene flow may result in the selection for increased plasticity. Here I expand frameworks for understanding relationships among selection, adaptation and gene flow to include the effects of phenotypic plasticity in natural populations, and highlight its importance in evolutionary diversification.     

The roles of allopatric divergence and natural selection in quantitative trait variation across a secondary contact zone in the lizard Anolis roquet

Populations of the Caribbean lizard, Anolis roquet, are thought to have experienced long periods of allopatry before recent secondary contact. To elucidate the effects of past allopatry on population divergence in A. roquet, we surveyed parallel transects across a secondary contact zone in northeastern Martinique. We used diagnostic molecular mitochondrial DNA markers to test fine‐scale association of mitochondrial DNA lineage and geological region, multivariate statistical techniques to explore quantitative trait pattern, and cline fitting techniques to model trait variation across the zone of secondary contact. We found that lineages were strongly associated with geological regions along both transects, but quantitative trait patterns were remarkably different. Patterns of morphological and mitochondrial DNA variation were consistent with a strong barrier to gene flow on the coast, whereas there were no indications of barriers to gene flow in the transitional forest. Hence, the coastal populations behaved as would be predicted by an allopatric model of divergence in this complex, while those in the transitional forest did not, despite the close proximity of the transects and their shared geological history. Patterns of geographical variation in this species complex, together with environmental data, suggest that on balance, selection regimes on either side of the secondary contact zone in the transitional forest may be more convergent, while those either side of the secondary contact zone on the coast are more divergent. Hence, the evolutionary consequences of allopatry may be strongly influenced by local natural selection regimes.     

Fine‐grained adaptive divergence in an amphibian: genetic basis of phenotypic divergence and the role of nonrandom gene flow in restricting effective migration among wetlands

Adaptive ecological differentiation among sympatric populations is promoted by environmental heterogeneity, strong local selection and restricted gene flow. High gene flow, on the other hand, is expected to homogenize genetic variation among populations and therefore prevent local adaptation. Understanding how local adaptation can persist at the spatial scale at which gene flow occurs has remained an elusive goal, especially for wild vertebrate populations. Here, we explore the roles of natural selection and nonrandom gene flow (isolation by breeding time and habitat choice) in restricting effective migration among local populations and promoting generalized genetic barriers to neutral gene flow. We examined these processes in a network of 17 breeding ponds of the moor frog Rana arvalis, by combining environmental field data, a common garden experiment and data on variation in neutral microsatellite loci and in a thyroid hormone receptor (TRβ) gene putatively under selection. We illustrate the connection between genotype, phenotype and habitat variation and demonstrate that the strong differences in larval life history traits observed in the common garden experiment can result from adaptation to local pond characteristics. Remarkably, we found that haplotype variation in the TRβ gene contributes to variation in larval development time and growth rate, indicating that polymorphism in the TRβ gene is linked with the phenotypic variation among the environments. Genetic distance in neutral markers was correlated with differences in breeding time and environmental differences among the ponds, but not with geographical distance. These results demonstrate that while our study area did not exceed the scale of gene flow, ecological barriers constrained gene flow among contrasting habitats. Our results highlight the roles of strong selection and nonrandom gene flow created by phenological variation and, possibly, habitat preferences, which together maintain genetic and phenotypic divergence at a fine‐grained spatial scale.     

Adaptive responses for seed and leaf phenology in natural populations of sessile oak along an altitudinal gradient

We assessed the adaptive potential of seed and leaf phenology in 10 natural populations of sessile oak (Quercus petraea) sampled along two altitudinal transects using common garden experiments. Population differentiation for both phenological traits was observed with high-altitude populations germinating and flushing later than low altitude ones. However, high genetic variation and heritability values were also maintained within populations, despite slightly decreasing for dates of leaf unfolding with increasing altitude. We suggest that biotic and abiotic fluctuating selection pressures within populations and high gene flow are the main mechanisms maintaining high genetic variation for these fitness related traits. Moreover, changes in selection intensity and/or selection pressures along the altitudinal gradient can explain the reduction in genetic variation observed for leaf phenology. We anticipate that the maintenance of high genetic variation will be a valuable resource for future adaptation of sessile oak populations undergoing an upslope shift caused by climate change.     

Landscape Genetic Structure of a Streamside Tree Species Euptelea pleiospermum (Eupteleaceae): Contrasting Roles of River Valley and Mountain Ridge

We used landscape genetics and statistical models to test how landscape features influence connectivity or create barriers to dispersal for a mountain riparian tree species, Euptelea pleiospermum. Young leaves from 1078 individuals belonging to 36 populations at elevations of 900–2000 m along upper reaches of four rivers were genotyped using eight nuclear microsatellite markers. We found no evidence for the unidirectional dispersal hypothesis in E. pleiospermum within each river. The linear dispersal pattern along each river valley is mostly consistent with the “classical metapopulaton” model. Mountain ridges separating rivers were genetic barriers for this wind-pollinated tree species with anemochorous seeds, whereas river valleys provided important corridors for dispersal. Gene flow among populations along elevational gradients within each river prevails over gene flow among populations at similar elevations but from different rivers. This pattern of gene flow is likely to promote elevational range shifts of plant populations and to hinder local adaptation along elevational gradients. This study provides a paradigm to determine which of the two strategies (migration or adaptation) will be adopted by mountain riparian plants under climate warming.     

Local adaptation at range edges: comparing elevation and latitudinal gradients

Local adaptation at range edges influences species’ distributions and how they respond to environmental change. However, the factors that affect adaptation, including gene flow and local selection pressures, are likely to vary across different types of range edge. We performed a reciprocal transplant experiment to investigate local adaptation in populations of Plantago lanceolata and P. major from central locations in their European range and from their latitudinal and elevation range edges (in northern Scandinavia and Swiss Alps, respectively). We also characterized patterns of genetic diversity and differentiation in populations using molecular markers. Range‐centre plants of P. major were adapted to conditions at the range centre, but performed similarly to range‐edge plants when grown at the range edges. There was no evidence for local adaptation when comparing central and edge populations of P. lanceolata. However, plants of both species from high elevation were locally adapted when compared with plants from high latitude, although the reverse was not true. This asymmetry was associated with greater genetic diversity and less genetic differentiation over the elevation gradient than over the latitudinal gradient. Our results suggest that adaptation in some range‐edge populations could increase their performance following climate change. However, responses are likely to differ along elevation and latitudinal gradients, with adaptation more likely at high‐elevation. Furthermore, based upon these results, we suggest that gene flow is unlikely to constrain adaptation in range‐edge populations of these species.     

Integrating patterns of thermal tolerance and phenotypic plasticity with population genetics to improve understanding of vulnerability to warming in a widespread copepod

Differences in population vulnerability to warming are defined by spatial patterns in thermal adaptation. These patterns may be driven by natural selection over spatial environmental gradients, but can also be shaped by gene flow, especially in marine taxa with high dispersal potential. Understanding and predicting organismal responses to warming requires disentangling the opposing effects of selection and gene flow. We begin by documenting genetic divergence of thermal tolerance and developmental phenotypic plasticity. Ten populations of the widespread copepod Acartia tonsa were collected from sites across a large thermal gradient, ranging from the Florida Keys to Northern New Brunswick, Canada (spanning over 20° latitude). Thermal performance curves (TPCs) from common garden experiments revealed local adaptation at the sampling range extremes, with thermal tolerance increasing at low latitudes and decreasing at high latitudes. The opposite pattern was observed in phenotypic plasticity, which was strongest at high latitudes. No relationship was observed between phenotypic plasticity and environmental variables. Instead, the results are consistent with the hypothesis of a trade‐off between thermal tolerance and the strength of phenotypic plasticity. Over a large portion of the sampled range, however, we observed a remarkable lack of differentiation of TPCs. To examine whether this lack of divergence is the result of selection for a generalist performance curve or constraint by gene flow, we analyzed cytochrome oxidase I mtDNA sequences, which revealed four distinct genetic clades, abundant genetic diversity, and widely distributed haplotypes. Strong divergence in thermal performance within genetic clades, however, suggests that the pace of thermal adaptation can be relatively rapid. The combined insight from the laboratory physiological experiments and genetic data indicate that gene flow constrains differentiation of TPCs. This balance between gene flow and selection has implications for patterns of vulnerability to warming. Taking both genetic differentiation and phenotypic plasticity into account, our results suggest that local adaptation does not increase vulnerability to warming, and that low‐latitude populations in general may be more vulnerable to predicted temperature change over the next century.     

GENE FLOW AND INEFFECTIVE ANTIPREDATOR BEHAVIOR IN A STREAM-BREEDING SALAMANDER

Predators often feed on prey that show ineffective antipredator behavior. Gene flow among populations may constrain evolution of effective antipredator ability in larvae of the streamside salamander, Ambystoma barbouri, a species that occupies distinctly different habitats with conflicting selection pressures. Some streams are ephemeral, where larvae should be active to feed and reach metamorphosis before stream drying. In contrast, other streams are more permanent and contain pools with predatory fish, where larvae should remain inactive to avoid fish predation. Feeding rates and predator escape behavior were assayed for laboratory-reared larvae from 15 populations. Larval survival was also compared among populations in artificial streams with natural predators. Five populations represented streams subjected to fish predation along a gradient of genetic and geographic isolation from populations without fish; the remaining 10 populations were ephemeral and without fish. Individuals from populations with fish had significantly stronger behavioral responses to fish (i.e., decreased feeding rate associated with the presence of fish and increased escape response) than individuals from fishless populations. Larvae from populations containing fish that were more isolated from fishless populations showed stronger antipredator responses than less isolated populations. Further, larvae from more isolated populations survived longer in the predation experiment, indicating that the behaviors measured were related with survival. These results suggest that gene flow between populations with conflicting selection pressures limits local adaptation in some salamander populations with fish. While previous studies have typically focused on the role of gene flow in pairs of populations, the results of this study suggest that gene flow is acting to swamp local adaptation across several populations.     

Local selection modifies phenotypic divergence among Rana temporaria populations in the presence of gene flow

In ectotherms, variation in life history traits among populations is common and suggests local adaptation. However, geographic variation itself is not a proof for local adaptation, as genetic drift and gene flow may also shape patterns of quantitative variation. We studied local and regional variation in means and phenotypic plasticity of larval life history traits in the common frog Rana temporaria using six populations from central Sweden, breeding in either open‐canopy or partially closed‐canopy ponds. To separate local adaptation from genetic drift, we compared differentiation in quantitative genetic traits (Q_ST) obtained from a common garden experiment with differentiation in presumably neutral microsatellite markers (F_ST). We found that R. temporaria populations differ in means and plasticities of life history traits in different temperatures at local, and in F_ST at regional scale. Comparisons of differentiation in quantitative traits and in molecular markers suggested that natural selection was responsible for the divergence in growth and development rates as well as in temperature‐induced plasticity, indicating local adaptation. However, at low temperature, the role of genetic drift could not be separated from selection. Phenotypes were correlated with forest canopy closure, but not with geographical or genetic distance. These results indicate that local adaptation can evolve in the presence of ongoing gene flow among the populations, and that natural selection is strong in this system.     

Climate change shifts natural selection and the adaptive potential of the perennial forb Boechera stricta in the Rocky Mountains

Heritable genetic variation is necessary for populations to evolve in response to anthropogenic climate change. However, antagonistic genetic correlations among traits may constrain the rate of adaptation, even if substantial genetic variation exists. We examine potential genetic responses to selection by comparing multivariate genetic variance–covariances of traits and fitness (multivariate Robertson–Price identities) across different environments in a reciprocal transplant experiment of the forb Boechera stricta in the Rocky Mountains. By transplanting populations into four common gardens arrayed along an elevational gradient, and exposing populations to control and snow removal treatments, we simulated future and current climates and snowmelt regimes. Genetic variation in flowering and germination phenology declined in plants moved downslope to warmer, drier sites, suggesting that these traits may have a limited ability to evolve under future climates. Simulated climate change via snow removal altered the strength of selection on flowering traits, but we found little evidence that genetic correlations among traits are likely to affect the rate of adaptation to climate change. Overall, our results suggest that climate change may alter the evolutionary potential of B. stricta, but reduced expression of genetic variation may be a larger impediment to adaptation than constraints imposed by antagonistic genetic correlations.     

Local adaptation, patterns of selection, and gene flow in the Californian serpentine sunflower (Helianthus exilis)

The traditional view of the species as the fundamental unit of evolution has been challenged by observations that in heterogeneous environments, gene flow may be too restricted to overcome the effects of local selection. Whether a species evolves as a cohesive unit depends critically on the dynamic balance between homogenizing gene flow among populations and potentially disruptive local adaptation. To examine this evolutionary balance between "global" gene flow and local selection, we studied northern Californian populations of Helianthus exilis, the serpentine sunflower, within a mosaic of contrasting serpentine and nonserpentine areas that differ considerably in soil chemistry and water availability. Local adaptation to riparian and serpentine habitats was studied in Helianthus exilis along with an analysis of gene flow patterns among populations within these habitats. Local adaptation was assessed in H. exilis during 2002 and 2003 using reciprocal transplant experiments at multiple locations within serpentine and riparian habitats. Effects of competition and germination date on the expression of local adaptation were also examined within the reciprocal transplant experiments. Local adaptation was detected in both years at the local site level and at the level of habitat. The analysis of the transplanted populations indicated that the patterns of selection differed considerably between riparian and serpentine sites. Differential survivorship occurred in serpentine habitats, whereas selection on reproductive output predominated in riparian habitats. Local adaptation was expressed only in the absence of competition. Local adaptation in terms of survivorship was most strongly expressed in treatments with delayed seed germination. Microsatellite markers were used to quantify population genetic parameters and examine the patterns of gene flow among sampled populations. Analysis of molecular markers revealed a system of population patches that freely exchange genes with each other. Strong selection seems to maintain ecotypic variation within this endemic sunflower species, while extensive gene flow among populations prevents local speciation between serpentine and riparian ecotypes.     

Evolutionary aspects of population structure for molecular and quantitative traits in the freshwater snail Radix balthica

Detecting the action of selection in natural populations can be achieved using the QST-FST comparison that relies on the estimation of FST with neutral markers, and QST using quantitative traits potentially under selection. QST higher than FST suggests the action of directional selection and thus potential local adaptation. In this article, we apply the QST-FST comparison to four populations of the hermaphroditic freshwater snail Radix balthica located in a floodplain habitat. In contrast to most studies published so far, we did not detect evidence of directional selection for local optima for any of the traits we measured: QST calculated using three different methods was never higher than FST. A strong inbreeding depression was also detected, indicating that outcrossing is probably predominant over selfing in the studied populations. Our results suggest that in this floodplain habitat, local adaptation of R. balthica populations may be hindered by genetic drift, and possibly altered by uneven gene flow linked to flood frequency.     

Speciation in a keystone plant genus is driven by elevation: a case study in New Guinean Ficus

Much of the world's insect and plant biodiversity is found in tropical and subtropical ‘hotspots’, which often include long elevational gradients. These gradients may function as ‘diversity pumps’ and contribute to both regional and local species richness. Climactic conditions on such gradients often change rapidly along short vertical distances and may result in local adaptation and high levels of population genetic structure in plants and insects. We investigated the population genetic structure of two species of Ficus (Moraceae) along a continuously forested elevational gradient in Papua New Guinea. This speciose plant genus is pollinated by tiny, species‐specific and highly coevolved chalcid wasps (Agaonidae) and represented by at least 73 species at our study gradient. We present results from two species of Ficus sampled from six elevations between 200 m and 2700 m a.s.l. (almost the entire elevational range of the genus) and 10 polymorphic microsatellite loci. These results show that strong barriers to gene flow exist between 1200 m and 1700 m a.s.l. Whereas lowland populations are panmictic across distances over 70 km, montane populations can be disjunct over 4 km, despite continuous forest cover. We suggest that the limited gene flow between populations of these two species of montane Ficus may be driven by environmental limitations on pollinator or seed dispersal in combination with local adaptation of Ficus populations. Such a mechanism may have wider implications for plant and pollinator speciation across long and continuously forested elevational gradients if generalist insect pollinators and vertebrate seed dispersers also form populations based on elevation.     

Gene flow in Dubautia arborea and D. ciliolata: the roles of ecology and isolation by distance in maintaining species boundaries despite ongoing hybridization

The relative roles of gene flow and natural selection in maintaining species differentiation have been a subject of debate for some time. The traditional view is that gene flow constrains adaptive divergence and maintains species cohesiveness. Alternatively, ecological speciation posits that the reverse is true: that adaptive ecological differentiation constrains gene flow. In this study, we examine gene flow and population differentiation among populations of two species of the Hawaiian silversword alliance, Dubautia arborea and D. ciliolata. We compare divergence in putatively neutral microsatellite markers with divergence in leaf morphometric traits, which may be selectively important or physiologically linked to selectively important traits. Gene flow between populations was found to be significant in only one of the two species, D. arborea. Leaf morphometric differentiation between species was significant, though not among populations within species. No evidence of effective genetic introgression was observed between apparently 'pure' populations of these species. Gene flow as measured by microsatellites was not correlated with geographic distance between populations, but was correlated with the linear placement of the widest part of the leaf. Because these two species are interfertile, as demonstrated by the presence of active hybrid zone, the lack of genetic introgression and the maintenance of species boundaries may be associated with natural selection on differential habitat.     

GENETIC DIFFERENTIATION AND SELECTION AGAINST MIGRANTS IN EVOLUTIONARILY REPLICATED EXTREME ENVIRONMENTS

We investigated mechanisms of reproductive isolation in livebearing fishes (genus Poecilia) inhabiting sulfidic and nonsulfidic habitats in three replicate river drainages. Although sulfide spring fish convergently evolved divergent phenotypes, it was unclear if mechanisms of reproductive isolation also evolved convergently. Using microsatellites, we found strongly reduced gene flow between adjacent populations from different habitat types, suggesting that local adaptation to sulfidic habitats repeatedly caused the emergence of reproductive isolation. Reciprocal translocation experiments indicate strong selection against immigrants into sulfidic waters, but also variation among drainages in the strength of selection against immigrants into nonsulfidic waters. Mate choice experiments revealed the evolution of assortative mating preferences in females from nonsulfidic but not from sulfidic habitats. The inferred strength of sexual selection against immigrants (RI_s) was negatively correlated with the strength of natural selection (RI_m), a pattern that could be attributed to reinforcement, whereby natural selection strengthens behavioral isolation due to reduced hybrid fitness. Overall, reproductive isolation and genetic differentiation appear to be replicated and direct consequences of local adaptation to sulfide spring environments, but the relative contributions of different mechanisms of reproductive isolation vary across these evolutionarily independent replicates, highlighting both convergent and nonconvergent evolutionary trajectories of populations in each drainage.