We investigated replicated ecological speciation in the livebearing fish Poecilia mexicana and P. sulphuraria (Poeciliidae), which inhabit freshwater habitats and have also colonized multiple sulfidic springs in southern Mexico. These springs exhibit extreme hypoxia and high concentrations of hydrogen sulfide, which is lethal to most metazoans. We used phylogenetic analyses to test whether springs were independently colonized, performed phenotypic assessments of body and gill morphology variation to identify convergent patterns of trait differentiation, and conducted an eco-toxicological experiment to detect differences in sulfide tolerances among ecotypes. Our results indicate that sulfidic springs were colonized by three different lineages, two within P. mexicana and one representing P. sulphuraria. Colonization occurred earlier in P. sulphuraria, whereas invasion of sulfidic springs in P. mexicana was more recent, such that each population is more closely related to neighboring populations from adjacent nonsulfidic habitats. Sulfide spring fish also show divergence from nonsulfidic phenotypes and a phenotypic convergence toward larger heads, larger gills, and increased tolerance to H 2 S. Together with previous studies that indicated significant reproductive isolation between fish from sulfidic and nonsulfidic habitats, this study provides evidence for repeated ecological speciation in the independent sulfide spring populations of P. mexicana and P. sulphuraria. K E Y W O R D S :Ecological speciation, hydrogen sulfide, hypoxia, local adaptation, morphological differentiation, mtDNA phylogeny, natural selection, Poecilia.

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.

The process of ecological speciation drives the evolution of locally adapted and reproductively isolated populations in response to divergent natural selection. In Southern Mexico, several lineages of the freshwater fish species of the genus Poecilia have independently colonized toxic, hydrogen sulfide-rich springs. Even though ecological speciation processes are increasingly well understood in this system, aligning the taxonomy of these fish with evolutionary processes has lagged behind. While some sulfide spring populations are classified as ecotypes of Poecilia mexicana, others, like P. sulphuraria, have been described as highly endemic species. Our study particularly focused on elucidating the taxonomy of the long described sulfide spring endemic, Poecilia thermalis Steindachner 1863, and investigates if similar evolutionary patterns of phenotypic trait divergence and reproductive isolation are present as observed in other sulfidic species of Poecilia. We applied a geometric morphometric approach to assess body shape similarity to other sulfidic and non-sulfidic fish of the genus Poecilia. We also conducted phylogenetic and population genetic analyses to establish the phylogenetic relationships of P. thermalis and used a population genetic approach to determine levels of gene flow among Poecilia from sulfidic and non-sulfidic sites. Our results indicate that P. thermalis' body shape has evolved in convergence with other sulfide spring populations in the genus. Phylogenetic analyses placed P. thermalis as most closely related to one population of P. sulphuraria, and population genetic analyses demonstrated that P. thermalis is genetically isolated from both P. mexicana ecotypes and P. sulphuraria. Based on these findings, we make taxonomic recommendations for P. thermalis. Overall, our study verifies the role of hydrogen sulfide as a main factor shaping convergent, phenotypic evolution and the emergence of reproductive isolation between Poecilia populations residing in adjacent sulfidic and non-sulfidic environments.

Hydrogen sulfide (H2S) is a potent toxicant interfering with oxidative phosphorylation in mitochondria and creating extreme environmental conditions in aquatic ecosystems. The mechanistic basis of adaptation to perpetual exposure to H2S remains poorly understood. We investigated evolutionarily independent lineages of livebearing fishes that have colonized and adapted to springs rich in H2S and compared their genome-wide gene expression patterns with closely related lineages from adjacent, nonsulfidic streams. Significant differences in gene expression were uncovered between all sulfidic and nonsulfidic population pairs. Variation in the number of differentially expressed genes among population pairs corresponded to differences in divergence times and rates of gene flow, which is consistent with neutral drift driving a substantial portion of gene expression variation among populations. Accordingly, there was little evidence for convergent evolution shaping large-scale gene expression patterns among independent sulfide spring populations. Nonetheless, we identified a small number of genes that was consistently differentially expressed in the same direction in all sulfidic and nonsulfidic population pairs. Functional annotation of shared differentially expressed genes indicated upregulation of genes associated with enzymatic H2S detoxification and transport of oxidized sulfur species, oxidative phosphorylation, energy metabolism, and pathways involved in responses to oxidative stress. Overall, our results suggest that modification of processes associated with H2S detoxification and toxicity likely complement each other to mediate elevated H2S tolerance in sulfide spring fishes. Our analyses allow for the development of novel hypotheses about biochemical and physiological mechanisms of adaptation to extreme environments.

Replicated ecological gradients are prime systems to study processes of molecular evolution underlying ecological divergence. Here, we investigated the repeated adaptation of the neotropical fish Poecilia mexicana to habitats containing toxic hydrogen sulphide (H2 S) and compared two population pairs of sulphide-adapted and ancestral fish by sequencing population pools of >200 individuals (Pool-Seq). We inferred the evolutionary processes shaping divergence and tested the hypothesis of increase of parallelism from SNPs to molecular pathways. Coalescence analyses showed that the divergence occurred in the face of substantial bidirectional gene flow. Population divergence involved many short, widely dispersed regions across the genome. Analyses of allele frequency spectra suggest that differentiation at most loci was driven by divergent selection, followed by a selection-mediated reduction of gene flow. Reconstructing allelic state changes suggested that selection acted mainly upon de novo mutations in the sulphide-adapted populations. Using a corrected Jaccard index to quantify parallel evolution, we found a negligible proportion of statistically significant parallel evolution of Jcorr  = 0.0032 at the level of SNPs, divergent genome regions (Jcorr  = 0.0061) and genes therein (Jcorr  = 0.0091). At the level of metabolic pathways, the overlap was Jcorr  = 0.2545, indicating increasing parallelism with increasing level of biological integration. The majority of pathways contained positively selected genes in both sulphide populations. Hence, adaptation to sulphidic habitats necessitated adjustments throughout the genome. The largely unique evolutionary trajectories may be explained by a high proportion of de novo mutations driving the divergence. Our findings favour Gould's view that evolution is often the unrepeatable result of stochastic events with highly contingent effects.