Extremophile Poeciliidae: multivariate insights into the complexity of speciation along replicated ecological gradients

Evolutionary Rate Shifts in Coding and Regulatory Regions Underpin Repeated Adaptation to Sulfidic Streams in Poeciliid Fishes

Adaptation to extreme environments often involves the evolution of dramatic physiological changes. To better understand how organisms evolve these complex phenotypic changes, the repeatability and predictability of evolution, and possible constraints on adapting to an extreme environment, it is important to understand how adaptive variation has evolved. Poeciliid fishes represent a particularly fruitful study system for investigations of adaptation to extreme environments due to their repeated colonization of toxic hydrogen sulfide-rich springs across multiple species within the clade. Previous investigations have highlighted changes in the physiology and gene expression in specific species that are thought to facilitate adaptation to hydrogen sulfide-rich springs. However, the presence of adaptive nucleotide variation in coding and regulatory regions and the degree to which convergent evolution has shaped the genomic regions underpinning sulfide tolerance across taxa are unknown. By sampling across seven independent lineages in which nonsulfidic lineages have colonized and adapted to sulfide springs, we reveal signatures of shared evolutionary rate shifts across the genome. We found evidence of genes, promoters, and putative enhancer regions associated with both increased and decreased convergent evolutionary rate shifts in hydrogen sulfide-adapted lineages. Our analysis highlights convergent evolutionary rate shifts in sulfidic lineages associated with the modulation of endogenous hydrogen sulfide production and hydrogen sulfide detoxification. We also found that regions with shifted evolutionary rates in sulfide spring fishes more often exhibited convergent shifts in either the coding region or the regulatory sequence of a given gene, rather than both.

Abiotic factors and heavy metals defining eco-physiological niche in fish

Ecosystem dynamics undergoing alterations in structure and function highlights the need to look into the relations between ecological parameters and organismal fitness and tolerance. Ecophysiological studies are used to understand how organisms adapt to and cope up with environmental stress. Current study uses a process-based approach to model physiochemical parameters regarding seven different fish species. Species respond to climatic variations via acclimation or adaptation through physiological plasticity. Four sites are differentiated into two types based on the water quality parameters and metal contamination. Seven fish species are clustered into two groups, each group depicting separate pattern of response in similar habitat. In this manner, biomarkers from three different physiological axes- stress, reproduction, and neurology were taken to determine the organism's ecological niche. Cortisol, Testosterone, Estradiol, and AChE are the signature molecules estimated for the said physiological axes. The ordination technique, nonmetric multidimensional scaling, has been utilized to visualize the differentiated physiological response to changing environmental conditions. Then, Bayesian Model Averaging (BMA) was used to identify the factors that play a key role in refining the stress physiology and determining the niche. Current study confirms different species belonging to similar habitats respond to various environmental and physiological factors in a different manner as various biomarkers respond in a species-specific pattern that induces the choice of habitat preference controlling its ecophysiological niche. In the present study, it is quite apparent that adaptive mechanism of fish to environmental stress is achieved through modification of physiological mechanisms through a panel of biochemical markers. These markers organize a cascade of physiological event at various levels including reproduction.

Phenotypic responses to oil pollution in a poeciliid fish

Pollution damages ecosystems around the globe and some forms of pollution, like oil pollution, can be either man-made or derived from natural sources. Despite the pervasiveness of oil pollution, certain organisms are able to colonise polluted or toxic environments, yet we only have a limited understanding of how they are affected by it. Here, we analysed phenotypic responses to oil pollution in guppies (Poecilia reticulata) living in oil-polluted habitats across southern Trinidad. We analysed body-shape and life-history traits for 352 individuals from 11 independent populations, six living in oil-polluted environments (including the naturally oil-polluted Pitch Lake), and five stemming from non-polluted habitats. Based on theory of, and previous studies on, responses to environmental stressors, we predicted guppies from oil-polluted waters to have larger heads and shallower bodies, to be smaller, to invest more into reproduction, and to produce more but smaller offspring compared to guppies from non-polluted habitats. Contrary to most of our predictions, we uncovered strong population-specific variation regardless of the presence of oil pollution. Moreover, guppies from oil-polluted habitats were characterised by increased body size; rounder, deeper bodies with increased head size; and increased offspring size, when compared to their counterparts from non-polluted sites. This suggests that guppies in oil-polluted environments are not only subject to the direct negative effects of oil pollution, but might gain some (indirect) benefits from other concomitant environmental factors, such as reduced predation and reduced parasite load. Our results extend our knowledge of organismal responses to oil pollution and highlight the importance of anthropogenic pollution as a source of environmental variation. They also emphasise the understudied ecological heterogeneity of extreme environments.

Multivariate analysis of morphology, behaviour, growth and developmental timing in hybrids brings new insights into the divergence of sympatric Arctic charr morphs

BACKGROUND

Studying the development of fitness related traits in hybrids from populations diverging in sympatry is a fundamental approach to understand the processes of speciation. However, such traits are often affected by covariance structures that complicate the comprehension of these processes, especially because the interactive relationships between traits of different nature (e.g. morphology, behaviour, life-history) remain largely unknown in this context. In a common garden setup, we conducted an extensive examination of a large suit of traits putatively involved in the divergence of two morphs of Arctic charr (Salvelinus alpinus), and investigated the consequences of potential patterns of trait covariance on the phenotype of their hybrids. These traits were measured along ontogeny and involved growth, yolk sac resorption, developmental timing (hatching and the onset of exogeneous feeding), head morphology and feeding behaviour.

RESULTS

Growth trajectories provided the strongest signal of phenotypic divergence between the two charr. Strikingly, the first-generation hybrids did not show intermediate nor delayed growth but were similar to the smallest morph, suggesting parental biases in the inheritance of growth patterns. However, we did not observe extensive multivariate trait differences between the two morphs and their hybrids. Growth was linked to head morphology (suggesting that morphological variations in early juveniles relate to simple allometric effects) but this was the only strong signal of covariance observed between all the measured traits. Furthermore, we did not report evidence for differences in overall phenotypic variance between morphs, nor for enhanced phenotypic variability in their hybrids.

CONCLUSION

Our study shed light on the multivariate aspect of development in a context of adaptive divergence. The lack of evidence for the integration of most traits into a single covariance structure suggested that phenotypic constraints may not always favour nor impede divergence toward ecological niches differing in numerous physical and ecological variables, as observed in the respective habitats of the two charr. Likewise, the role of hybridization as a disruptive agent of trait covariance may not necessarily be significant in the evolution of populations undergoing resource polymorphism.

A Novel Body Plan Alters Diversification of Body Shape and Genitalia in Live-Bearing Fish

Major evolutionary innovations can greatly influence subsequent evolution. While many major transitions occurred in the deep past, male live-bearing fishes (family Poeciliidae) more recently evolved a novel body plan. This group possesses a three-region axial skeleton, with one region—the ano-urogenital region—representing a unique body region accommodating male genitalic structures (gonopodial complex). Here we evaluate several hypotheses for the evolution of diversity in this region and examine its role in the evolution of male body shape. Examining Gambusia fishes, we tested a priori predictions for (1) joint influence of gonopodial-complex traits on mating performance, (2) correlated evolution of gonopodial-complex traits at macro- and microevolutionary scales, and (3) predator-driven evolution of gonopodial-complex traits in a post-Pleistocene radiation of Bahamas mosquitofish. We found the length of the sperm-transfer organ (gonopodium) and its placement along the body (gonopodial anterior transposition) jointly influenced mating success, with correlational selection favoring particular trait combinations. Despite these two traits functionally interacting during mating, we found no evidence for their correlated evolution at macro- or microevolutionary scales. In contrast, we did uncover correlated evolution of modified vertebral hemal spines (part of the novel body region) and gonopodial anterior transposition at both evolutionary scales, matching predictions of developmental connections between these components. Developmental linkages in the ano-urogenital region apparently play key roles in evolutionary trajectories, but multiple selective agents likely act on gonopodium length and cause less predictable evolution. Within Bahamas mosquitofish, evolution of hemal-spine morphology, and gonopodial anterior transposition across predation regimes was quite predictable, with populations evolving under high predation risk showing more modified hemal spines with greater modifications and a more anteriorly positioned gonopodium. These changes in the ano-urogenital vertebral region have facilitated adaptive divergence in swimming abilities and body shape between predation regimes. Gonopodium surface area, but not length, evolved as predicted in Bahamas mosquitofish, consistent with a previously suggested tradeoff between natural and sexual selection on gonopodium size. These results provide insight into how restructured body plans offer novel evolutionary solutions. Here, a novel body region—originally evolved to aid sperm transfer—was apparently co-opted to alter whole-organism performance, facilitating phenotypic diversification.

Does evolution of echolocation calls and morphology in Molossus result from convergence or stasis?

Although many processes of diversification have been described to explain variation of morphological traits within clades that have obvious differentiation among taxa, not much is known about these patterns in complexes of cryptic species. Molossus is a genus of bats that is mainly Neotropical, occurring from the southeastern United States to southern Argentina, including the Caribbean islands. Molossus comprises some groups of species that are morphologically similar but phylogenetically divergent, and other groups of species that are genetically similar but morphologically distinct. This contrast allows investigation of unequal trait diversification and the evolution of morphological and behavioural characters. In this study, we assessed the role of phylogenetic history in a genus of bat with three cryptic species complexes, and evaluated if morphology and behavior are evolving concertedly. The Genotype by Sequence genomic approach was used to build a species-level phylogenetic tree for Molossus and to estimate the ancestral states of morphological and echolocation call characters. We measured the correlation of phylogenetic distances to morphological and echolocation distances, and tested the relationship between morphology and behavior when the effect of phylogeny is removed. Morphology evolved via a mosaic of convergence and stasis, whereas call design was influenced exclusively through local adaptation and convergent evolution. Furthermore, the frequency of echolocation calls is negatively correlated with the size of the bat, but other characters do not seem to be evolving in concert. We hypothesize that slight variation in both morphology and behaviour among species of the genus might result from niche specialization, and that traits evolve to avoid competition for resources in similar environments.

A century later: Adaptive plasticity and rapid evolution contribute to geographic variation in invasive mosquitofish

One century after their introduction to Europe, eastern mosquitofish (Gambusia holbrooki) represent a natural experiment to determine the relative contributions of adaptive plasticity and rapid evolutionary change in creating large-scale geographic variation in phenotypes. We evaluated the population-genetic structure and invasion history based on allele length polymorphisms of 15 nuclear microsatellites, which we quantified for N = 660 individuals from 23 populations sampled in 2013 across the invasive range of G. holbrooki in Europe. We analysed body-shape and life-history variation in N = 1331 individuals from 36 populations, sampled in 2013 and 2017, and tested heritability of phenotypic differences in a subset of four populations using a common-garden experiment. The genetic structure of wild-caught individuals suggested a single introduction for all European mosquitofish, which were genetically impoverished compared to their native counterparts. We found some convergent patterns of phenotypic divergence across native and invasive climatic gradients (e.g., increased body size in colder/more northern populations); however, several phenotypic responses were not consistent between sampling years, pointing towards plastic phenotypes. Our analysis of common-garden reared individuals uncovered moderate heritability estimates only for two measures of male body size (intraclass correlation coefficient, ICC = 0.628 and 0.556) and offspring fat content (ICC = 0.734), while suggesting high levels of plasticity in most other phenotypic traits (ICC ≤ 0.407). Our results highlight the importance of phenotypic plasticity in invasive species during range expansions and demonstrate that strong selective pressures-in this case towards increased body size in colder environments-simultaneously promote rapid evolutionary divergence.

Effect of salinity on zinc toxicity (ZnCl2 and ZnO nanomaterials) in the mosquitofish (Gambusia sexradiata)

Zn is an essential trace metal in living beings. However, excessive concentrations can cause toxic effects even in the aquatic biota. Zn is widely used in different industrial sectors, which has increased its presence in aquatic environments. To assess the acute toxicity of Zn, bioassays were performed with the fish Gambusia sexradiata for a 96-h exposure using ZnCl₂ (0 and 15 salinity) and ZnO nanomaterials (0 salinity). The mean lethal concentrations (LC₅₀-96 h) for ZnCl₂ were 25.36 (19.64-32.76) and 177.91 (129.39-244.63) mg Zn L^-1 to 0 and 15 salinity, respectively. The increased concentration of ZnCl₂ showed a dose-response relationship; similarly, the increase in salinity significantly reduces the toxicity of Zn. Characterisation of ZnO nanomaterials was carried out by FTIR, DRX, SEM, DLS and zeta potential. The FTIR spectra showed the characteristic band of Zn-O vibration at 364 cm^-1, while DRX presents the hexagonal wurtzite structure with an average crystallite size of 40 nm. SEM micrographs reveal rod-like shapes with lengths and diameters of 40-350 nm and 90 nm, respectively. Agglomerates of 423 nm in water suspension were obtained by DLS and zeta potential of + 14.4 mV. Under these conditions, no mortality was observed due to the rapid flocculation/precipitation of ZnO nanomaterials, which involved brief interaction periods of Zn in the water column with the fish. Gambusia sexradiata is affected by increased Zn concentrations in hard water conditions, and salinity changes modified Zn toxicity, placing it as a suitable model for toxicity tests for this type of particles.

Mitochondria and the Origin of Species: Bridging Genetic and Ecological Perspectives on Speciation Processes

Mitochondria have been known to be involved in speciation through the generation of Dobzhansky-Muller incompatibilities, where functionally neutral co-evolution between mitochondrial and nuclear genomes can cause dysfunction when alleles are recombined in hybrids. We propose that adaptive mitochondrial divergence between populations can not only produce intrinsic (Dobzhansky-Muller) incompatibilities, but could also contribute to reproductive isolation through natural and sexual selection against migrants, post-mating prezygotic isolation, as well as by causing extrinsic reductions in hybrid fitness. We describe how these reproductive isolating barriers can potentially arise through adaptive divergence of mitochondrial function in the absence of mito-nuclear coevolution, a departure from more established views. While a role for mitochondria in the speciation process appears promising, we also highlight critical gaps of knowledge: (1) many systems with a potential for mitochondrially-mediated reproductive isolation lack crucial evidence directly linking reproductive isolation and mitochondrial function; (2) it often remains to be seen if mitochondrial barriers are a driver or a consequence of reproductive isolation; (3) the presence of substantial gene flow in the presence of mito-nuclear incompatibilities raises questions whether such incompatibilities are strong enough to drive speciation to completion; and (4) it remains to be tested how mitochondrial effects on reproductive isolation compare when multiple mechanisms of reproductive isolation coincide. We hope this perspective and the proposed research plans help to inform future studies of mitochondrial adaptation in a manner that links genotypic changes to phenotypic adaptations, fitness, and reproductive isolation in natural systems, helping to clarify the importance of mitochondria in the formation and maintenance of biological diversity.

Evolution in caves: selection from darkness causes spinal deformities in teleost fishes

Only few fish species have successfully colonized subterranean habitats, but the underlying biological constraints associated with this are still poorly understood. Here, we investigated the influence of permanent darkness on spinal-column development in one species (Midas cichlid, Amphilophus citrinellus) with no known cave form, and one (Atlantic molly, Poecilia mexicana) with two phylogenetically young cave forms. Specifically, fish were reared under a normal light : dark cycle or in permanent darkness (both species). We also surveyed wild-caught cave and surface ecotypes of P. mexicana In both species, permanent darkness was associated with significantly higher rates of spinal deformities (especially in A. citrinellus). This suggests strong developmental (intrinsic) constraints on the successful colonization of subterranean environments in teleost fishes and might help explain the relative paucity of cave-adapted lineages. Our results add depth to our understanding of the aspects of selection driving trait divergence and maintaining reproductive isolation in cave faunas.

Small-scale phenotypic differentiation along complex stream gradients in a non-native amphipod

Background

Selective landscapes in rivers are made up by an array of selective forces that vary from source to downstream regions or between seasons, and local/temporal variation in fitness maxima can result in gradual spatio-temporal variation of phenotypic traits. This study aimed at establishing freshwater amphipods as future model organisms to study adaptive phenotypic diversification (evolutionary divergence and/or adaptive plasticity) along stream gradients.

Methods

We collected Gammarus roeselii from 16 sampling sites in the Rhine catchment during two consecutive seasons (summer and winter). Altogether, we dissected n = 1648 individuals and quantified key parameters related to morphological and life-history diversification, including naturally selected (e.g., gill surface areas) as well as primarily sexually selected traits (e.g., male antennae). Acknowledging the complexity of selective regimes in streams and the interrelated nature of selection factors, we assessed several abiotic (e.g., temperature, flow velocity) and biotic ecological parameters (e.g., conspecific densities, sex ratios) and condensed them into four principal components (PCs).

Results

Generalized least squares models revealed pronounced phenotypic differentiation in most of the traits investigated herein, and components of the stream gradient (PCs) explained parts of the observed differences. Depending on the trait under investigation, phenotypic differentiation could be ascribed to variation in abiotic conditions, anthropogenic disturbance (influx of thermally polluted water), or population parameters. For example, female fecundity showed altitudinal variation and decreased with increasing conspecific densities, while sexual dimorphism in the length of male antennae—used for mate finding and assessment—increased with increasing population densities and towards female-biased sex ratios.

Conclusions

We provide a comprehensive protocol for comparative analyses of intraspecific variation in life history traits in amphipods. Whether the observed phenotypic differentiation over small geographical distances reflects evolutionary divergence or plasticity (or both) remains to be investigated in future studies. Independent of the mechanisms involved, variation in several traits is likely to have consequences for ecosystem functions. For example, leaf-shredding in G. roeselii strongly depends on body size, which varied in dependence of several ecological parameters.

Electronic supplementary material

The online version of this article (10.1186/s12983-019-0327-8) contains supplementary material, which is available to authorized users.

Testes mass in the livebearing fish Brachyrhaphis rhabdophora (Poeciliidae) varies hypoallometrically with body size but not between predation environments

In this study, we considered potential causes of variation in testis size in the livebearing fish Brachyrhaphis rhabdophora. We evaluated variation in testes mass among individual males and among populations that occupy different selective environments. First, we predicted that small males should allocate more to testes mass than large males (i.e., hypoallometric pattern) based on a sperm competition argument. Second, based on life history theory and associated differences in mortality rates between populations that coexist with many fish predators and those with few predators, we predicted that males in high-predation environments should allocate more to testes mass than males in habitats with few predators. Our results showed that small males allocated proportionally more to testes mass than larger males (slope of testes mass to body mass was hypoallometric). However, there was no effect of predator environment on testes mass independent of body size differences. In this system, size-specific patterns of reproductive allocation in males (hypoallometry) differ from that seen in females (hyperallometry). Allocation to testes mass may respond to differences in mortality rate through selection on body size.

Evolution of morphological crypsis in the Tetramorium caespitum ant species complex (Hymenoptera: Formicidae)

Cryptic species are morphologically very similar to each other. To what extent stasis or convergence causes crypsis and whether ecology influences the evolution of crypsis has remained unclear. The Tetramorium caespitum complex is one of the most intricate examples of cryptic species in ants. Here, we test three hypotheses concerning the evolution of its crypsis: H1: The complex is monophyletic. H2: Morphology resulted from evolutionary stasis. H3: Ecology and morphology evolved concertedly. We confirmed (H1) monophyly of the complex; (H2) a positive relation between morphological and phylogenetic distances, which indicates a very slow loss of similarity over time and thus stasis; and (H3) a positive relation between only one morphological character and a proxy of the ecological niche, which indicates concerted evolution of these two characters, as well as a negative relation between p-values of correct species identification and altitude, which suggests that species occurring in higher altitudes are more cryptic. Our data suggest that species-specific morphological adaptations to the ecological niche are exceptions in the complex, and we consider the worker morphology in this complex as an adaptive solution for various environments.

Female Choice Undermines the Emergence of Strong Sexual Isolation between Locally Adapted Populations of Atlantic Mollies (Poecilia mexicana)

Divergent selection between ecologically dissimilar habitats promotes local adaptation, which can lead to reproductive isolation (RI). Populations in the Poecilia mexicana species complex have independently adapted to toxic hydrogen sulfide and show varying degrees of RI. Here, we examined the variation in the mate choice component of prezygotic RI. Mate choice tests across drainages (with stimulus males from another drainage) suggest that specific features of the males coupled with a general female preference for yellow color patterns explain the observed variation. Analyses of male body coloration identified the intensity of yellow fin coloration as a strong candidate to explain this pattern, and common-garden rearing suggested heritable population differences. Male sexual ornamentation apparently evolved differently across sulfide-adapted populations, for example because of differences in natural counterselection via predation. The ubiquitous preference for yellow color ornaments in poeciliid females likely undermines the emergence of strong RI, as female discrimination in favor of own males becomes weaker when yellow fin coloration in the respective sulfide ecotype increases. Our study illustrates the complexity of the (partly non-parallel) pathways to divergence among replicated ecological gradients. We suggest that future work should identify the genomic loci involved in the pattern reported here, making use of the increasing genomic and transcriptomic datasets available for our study system.

Molecular evolution and expression of oxygen transport genes in livebearing fishes (Poeciliidae) from hydrogen sulfide rich springs

Hydrogen sulfide (H2S) is a natural toxicant in some aquatic environments that has diverse molecular targets. It binds to oxygen transport proteins, rendering them non-functional by reducing oxygen-binding affinity. Hence, organisms permanently inhabiting H2S-rich environments are predicted to exhibit adaptive modifications to compensate for the reduced capacity to transport oxygen. We investigated 10 lineages of fish of the family Poeciliidae that have colonized freshwater springs rich in H2S—along with related lineages from non-sulfidic environments—to test hypotheses about the expression and evolution of oxygen transport genes in a phylogenetic context. We predicted shifts in the expression of and signatures of positive selection on oxygen transport genes upon colonization of H2S-rich habitats. Our analyses indicated significant shifts in gene expression for multiple hemoglobin genes in lineages that have colonized H2S-rich environments, and three hemoglobin genes exhibited relaxed selection in sulfidic compared to non-sulfidic lineages. However, neither changes in gene expression nor signatures of selection were consistent among all lineages in H2S-rich environments. Oxygen transport genes may consequently be predictable targets of selection during adaptation to sulfidic environments, but changes in gene expression and molecular evolution of oxygen transport genes in H2S-rich environments are not necessarily repeatable across replicated lineages.

Shared and unique patterns of phenotypic diversification along a stream gradient in two congeneric species

Stream ecosystems show gradual variation of various selection factors, which can result in a zonation of species distributions and gradient evolution of morphological and life-history traits within species. Identifying the selective agents underlying such phenotypic evolution is challenging as different species could show shared and/or unique (species-specific) responses to components of the river gradient. We studied a stream gradient inhabited by two mosquitofishes (genus Gambusia) in the Río Grijalva basin in southern Mexico and found a patchy distribution pattern of both congeners along a stretch of 100 km, whereby one species was usually dominant at a given site. We uncovered both shared and unique patterns of diversification: some components of the stream gradient, including differences in piscine predation pressure, drove shared patterns of phenotypic divergence, especially in females. Other components of the gradient, particularly abiotic factors (max. annual temperature and temperature range) resulted in unique patterns of divergence, especially in males. Our study highlights the complexity of selective regimes in stream ecosystems. It exemplifies that even closely related, congeneric species can respond in unique ways to the same components of the river gradient and shows how both sexes can exhibit quite different patterns of divergence in multivariate phenotypic character suites.

The Evolutionary Ecology of Animals Inhabiting Hydrogen Sulfide–Rich Environments

Hydrogen sulfide (H2S) is a respiratory toxicant that creates extreme environments tolerated by few organisms. H2S is also produced endogenously by metazoans and plays a role in cell signaling. The mechanisms of H2S toxicity and its physiological functions serve as a basis to discuss the multifarious strategies that allow animals to survive in H2S-rich environments. Despite their toxicity, H2S-rich environments also provide ecological opportunities, and complex selective regimes of covarying abiotic and biotic factors drive trait evolution in organisms inhabiting H2S-rich environments. Furthermore, adaptation to H2S-rich environments can drive speciation, giving rise to biodiversity hot spots with high levels of endemism in deep-sea hydrothermal vents, cold seeps, and freshwater sulfide springs. The diversity of H2S-rich environments and their inhabitants provides ideal systems for comparative studies of the effects of a clear-cut source of selection across vast geographic and phylogenetic scales, ultimately informing our understanding of how environmental stressors affect ecological and evolutionary processes.

Extremophile Poeciliidae: multivariate insights into the complexity of speciation along replicated ecological gradients

Background

Replicate population pairs that diverge in response to similar selective regimes allow for an investigation of (a) whether phenotypic traits diverge in a similar and predictable fashion, (b) whether there is gradual variation in phenotypic divergence reflecting variation in the strength of natural selection among populations, (c) whether the extent of this divergence is correlated between multiple character suites (i.e., concerted evolution), and (d) whether gradual variation in phenotypic divergence predicts the degree of reproductive isolation, pointing towards a role for adaptation as a driver of (ecological) speciation. Here, we use poeciliid fishes of the genera Gambusia and Poecilia that have repeatedly evolved extremophile lineages able to tolerate high and sustained levels of toxic hydrogen sulfide (H₂S) to answer these questions.

Results

We investigated evolutionary divergence in response to H₂S in Gambusia spp. (and to a lesser extent Poecilia spp.) using a multivariate approach considering the interplay of life history, body shape, and population genetics (nuclear miscrosatellites to infer population genetic differentiation as a proxy for reproductive isolation). We uncovered both shared and unique patterns of evolution: most extremophile Gambusia predictably evolved larger heads and offspring size, matching a priori predictions for adaptation to sulfidic waters, while variation in adult life histories was idiosyncratic. When investigating patterns for both genera (Gambusia and Poecilia), we found that divergence in offspring-related life histories and body shape were positively correlated across populations, but evidence for individual-level associations between the two character suites was limited, suggesting that genetic linkage, developmental interdependencies, or pleiotropic effects do not explain patterns of concerted evolution. We further found that phenotypic divergence was positively correlated with both environmental H₂S-concentration and neutral genetic differentiation (a proxy for gene flow).

Conclusions

Our results suggest that higher toxicity exerts stronger selection, and that divergent selection appears to constrain gene flow, supporting a scenario of ecological speciation. Nonetheless, progress toward ecological speciation was variable, partially reflecting variation in the strength of divergent selection, highlighting the complexity of selective regimes even in natural systems that are seemingly governed by a single, strong selective agent.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0705-1) contains supplementary material, which is available to authorized users.