Colonisation of toxic environments drives predictable life-history evolution in livebearing fishes (Poeciliidae)

Integrative analyses of convergent adaptation in sympatric extremophile fishes

The evolution of independent lineages along replicated environmental transitions frequently results in convergent adaptation, yet the degree to which convergence is present across multiple levels of biological organization is often unclear. Additionally, inherent biases associated with shared ancestry and variation in selective regimes across geographic replicates often pose challenges for confidently identifying patterns of convergence. We investigated a system in which three species of poeciliid fishes sympatrically occur in a toxic spring rich in hydrogen sulfide (H2S) and an adjacent nonsulfidic stream to examine patterns of adaptive evolution across levels of biological organization. We found convergence in morphological and physiological traits and genome-wide patterns of gene expression among all three species. In addition, there were shared signatures of selection on genes encoding H2S toxicity targets in the mitochondrial genomes of each species. However, analyses of nuclear genomes revealed neither evidence for substantial genomic islands of divergence around genes involved in H2S toxicity and detoxification nor substantial congruence of strongly differentiated regions across population pairs. These non-convergent, heterogeneous patterns of genomic divergence may indicate that sulfide tolerance is highly polygenic, with shared allele frequency shifts present at many loci with small effects along the genome. Alternatively, H2S tolerance may involve substantial genetic redundancy, with non-convergent, lineage-specific variation at multiple loci along the genome underpinning similar changes in phenotypes and gene expression. Overall, we demonstrate variability in the extent of convergence across organizational levels and highlight the challenges of linking patterns of convergence across scales.

Thermal tolerance in an extremophile fish from Mexico is not affected by environmental hypoxia

The thermal ecology of ectotherm animals has gained considerable attention in the face of human-induced climate change. Particularly in aquatic species, the experimental assessment of critical thermal limits (CTmin and CTmax) may help to predict possible effects of global warming on habitat suitability and ultimately species survival. Here we present data on the thermal limits of two endemic and endangered extremophile fish species, inhabiting a geothermally heated and sulfur-rich spring system in southern Mexico: The sulfur molly (Poecilia sulphuraria) and the widemouth gambusia (Gambusia eurystoma). Besides physiological challenges induced by toxic hydrogen sulfide and related severe hypoxia during the day, water temperatures have been previously reported to exceed those of nearby clearwater streams. We now present temperature data for various locations and years in the sulfur spring complex and conducted laboratory thermal tolerance tests (CTmin and CTmax) both under normoxic and severe hypoxic conditions in both species. Average CTmax limits did not differ between species when dissolved oxygen was present. However, critical temperature (CTmax=43.2°C) in P. sulphuraria did not change when tested under hypoxic conditions, while G. eurystoma on average had a lower CTmax when oxygen was absent. Based on this data we calculated both species' thermal safety margins and used a TDT (thermal death time) model framework to relate our experimental data to observed temperatures in the natural habitat. Our findings suggest that both species live near their thermal limits during the annual dry season and are locally already exposed to temperatures above their critical thermal limits. We discuss these findings in the light of possible physiological adaptions of the sulfur-adapted fish species and the anthropogenic threats for this unique system.

Selection on standing genetic variation mediates convergent evolution in extremophile fish

Hydrogen sulfide is a toxic gas that disrupts numerous biological processes, including energy production in the mitochondria, yet fish in the Poecilia mexicana species complex have independently evolved sulfide tolerance several times. Despite clear evidence for convergence at the phenotypic level in these fishes, it is unclear if the repeated evolution of hydrogen sulfide tolerance is the result of similar genomic changes. To address this gap, we used a targeted capture approach to sequence genes associated with sulfide processes and toxicity from five sulfidic and five nonsulfidic populations in the species complex. By comparing sequence variation in candidate genes to a reference set, we identified similar population structure and differentiation, suggesting that patterns of variation in most genes associated with sulfide processes and toxicity are due to demographic history and not selection. But the presence of tree discordance for a subset of genes suggests that several loci are evolving divergently between ecotypes. We identified two differentiation outlier genes that are associated with sulfide detoxification in the mitochondria that have signatures of selection in all five sulfidic populations. Further investigation into these regions identified long, shared haplotypes among sulfidic populations. Together, these results reveal that selection on standing genetic variation in putatively adaptive genes may be driving phenotypic convergence in this species complex.

Whole-genome sequencing illuminates multifaceted targets of selection to humic substances in Eurasian perch

Extreme environments are inhospitable to the majority of species, but some organisms are able to survive in such hostile conditions due to evolutionary adaptations. For example, modern bony fishes have colonized various aquatic environments, including perpetually dark, hypoxic, hypersaline and toxic habitats. Eurasian perch (Perca fluviatilis) is among the few fish species of northern latitudes that is able to live in very acidic humic lakes. Such lakes represent almost "nocturnal" environments; they contain high levels of dissolved organic matter, which in addition to creating a challenging visual environment, also affects a large number of other habitat parameters and biotic interactions. To reveal the genomic targets of humic-associated selection, we performed whole-genome sequencing of perch originating from 16 humic and 16 clear-water lakes in northern Europe. We identified over 800,000 single nucleotide polymorphisms, of which >10,000 were identified as potential candidates under selection (associated with >3000 genes) using multiple outlier approaches. Our findings suggest that adaptation to the humic environment may involve hundreds of regions scattered across the genome. Putative signals of adaptation were detected in genes and gene families with diverse functions, including organism development and ion transportation. The observed excess of variants under selection in regulatory regions highlights the importance of adaptive evolution via regulatory elements, rather than via protein sequence modification. Our study demonstrates the power of whole-genome analysis to illuminate the multifaceted nature of humic adaptation and provides the foundation for further investigation of causal mutations underlying phenotypic traits of ecological and evolutionary importance.

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.

Epigenetic inheritance of DNA methylation changes in fish living in hydrogen sulfide-rich springs

Environmental factors can promote phenotypic variation through alterations in the epigenome and facilitate adaptation of an organism to the environment. Although hydrogen sulfide is toxic to most organisms, the fish Poecilia mexicana has adapted to survive in environments with high levels that exceed toxicity thresholds by orders of magnitude. Epigenetic changes in response to this environmental stressor were examined by assessing DNA methylation alterations in red blood cells, which are nucleated in fish. Males and females were sampled from sulfidic and nonsulfidic natural environments; individuals were also propagated for two generations in a nonsulfidic laboratory environment. We compared epimutations between the sexes as well as field and laboratory populations. For both the wild-caught (F0) and the laboratory-reared (F2) fish, comparing the sulfidic and nonsulfidic populations revealed evidence for significant differential DNA methylation regions (DMRs). More importantly, there was over 80% overlap in DMRs across generations, suggesting that the DMRs have stable generational inheritance in the absence of the sulfidic environment. This is an example of epigenetic generational stability after the removal of an environmental stressor. The DMR-associated genes were related to sulfur toxicity and metabolic processes. These findings suggest that adaptation of P. mexicana to sulfidic environments in southern Mexico may, in part, be promoted through epigenetic DNA methylation alterations that become stable and are inherited by subsequent generations independent of the environment.

Predation and Resource Availability Interact to Drive Life-History Evolution in an Adaptive Radiation of Livebearing Fish

Predation risk and resource availability are two primary factors predicted by theory to drive the evolution of life histories. Yet, disentangling their roles in life-history evolution in the wild is challenging because (1) the two factors often co-vary across environments, and (2) environmental effects on phenotypes can mask patterns of genotypic evolution. Here, we use the model system of the post-Pleistocene radiation of Bahamas mosquitofish (Gambusia hubbsi) inhabiting blue holes to provide a strong test of the roles of predation and resources in life-history evolution, as the two factors do not co-vary in this system and we attempted to minimize environmental effects by raising eight populations under common laboratory conditions. We tested a priori predictions of predation- and resource-driven evolution in five life-history traits. We found that life-history evolution in Bahamas mosquitofish largely reflected complex interactions in the effects of predation and resource availability. High predation risk has driven the evolution of higher fecundity, smaller offspring size, more frequent reproduction, and slower growth rate—but this predation-driven divergence primarily occurred in environments with relatively high resource availability, and the effects of resources on life-history evolution was generally greater within environments having high predation risk. This implies that resource-driven selection on life histories overrides selection from predators when resources are particularly scarce. While several results matched a priori predictions, with the added nuance of interdependence among selective agents, some did not. For instance, only resource levels, not predation risk, explained evolutionary change in male age at maturity, with more rapid sexual maturation in higher-resource environments. We also found faster (not slower) juvenile growth rates within low-resource and low-predation environments, probably caused by selection in these high-competition scenarios favoring greater growth efficiency. Our approach, using common-garden experiments with a natural system of low- and high-predation populations that span a continuum of resource availability, provides a powerful way to deepen our understanding of life-history evolution. Overall, it appears that life-history evolution in this adaptive radiation has resulted from a complex interplay between predation and resources, underscoring the need for increased attention on more sophisticated interactions among selective agents in driving phenotypic diversification.

What Drives Life-History Variation in the Livebearing Fish Poeciliopsis prolifica? An Assessment of Multiple Putative Selective Agents

Life-history traits are directly linked to fitness, and therefore, can be highly adaptive. Livebearers have been used as models for understanding the evolution of life histories due to their wide diversity in these traits. Several different selective pressures, including population density, predation, and resource levels, can shape life-history traits. However, these selective pressures are usually considered independently in livebearers and we lack a clear understanding of how they interact in shaping life-history evolution. Furthermore, selective pressures such as interspecific competition are rarely considered as drivers of life-history evolution in poeciliids. Here we test the simultaneous effects of several potential selective pressures on life-history traits in the livebearing fish Poeciliopsis prolifica. We employ a multi-model inference approach. We focus on four known agents of selection: resource availability, stream velocity, population density, and interspecific competition, and their effect on four life-history traits: reproductive allocation, superfetation, number of embryos, and individual embryo size. We found that models with population density and interspecific competition alone were strongly supported in our data and, hence, indicated that these two factors are the most important selective agents for most life-history traits, except for embryo size. When population density and interspecific competition increase there is an increase in each of the three life-history traits (reproductive allocation, superfetation, and number of embryos). For individual embryo size, we found that all single-agent models were equivalent and it was unclear which selective agent best explained variation. We also found that models that included population density and interspecific competition as direct effects were better supported than those that included them as indirect effects through their influence on resource availability. Our study underscores the importance of interspecific competitive interactions on shaping life-history traits and suggests that these interactions should be considered in future life-history studies.

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.

Molecular and morphological convergence to sulfide-tolerant fishes in a new species of Jenynsia (Cyprinodontiformes: Anablepidae), the first extremophile member of the family

Freshwater sulfide springs have extreme environmental conditions that only few vertebrate species can tolerate. These species often develop a series of morphological and molecular adaptations to cope with the challenges of life under the toxic and hypoxic conditions of sulfide springs. In this paper, we described a new fish species of the genus Jenynsia, Anablepidae, from a sulfide spring in Northwestern Argentina, the first in the family known from such extreme environment. Jenynsia sulfurica n. sp. is diagnosable by the lack of scales on the pre-pelvic area or the presence of a single row of scales, continuous or not, from the isthmus to the bases of the pelvic fins. Additionally, it presents a series of morphological and molecular characteristics that appear convergent with those seen in other fish species (e.g., Poeciliids) inhabiting sulfide springs. Most notably, J. sulfurica has an enlarged head and postorbital area compared to other fish of the genus and a prognathous lower jaw with a hypertrophied lip, thought to facilitate respiration at the air-water interface. Analyses of cox1 sequence showed that J. sulfurica has two unique mutations resulting in amino acid substitutions convergent to those seen in Poeciliids from sulfide springs and known to provide a physiological mechanism related to living in sulfide environments. A phylogenetic analysis, including molecular and morphological characters, placed J. sulfurica as sister taxa to J. alternimaculata, a species found in nearby, non-sulfide habitats directly connected to the sulfide springs. Thus, it can be inferred that the selection imposed by the presence of H2S has resulted in the divergence between these two species and has potentially served as a barrier to gene flow.

Local ancestry analysis reveals genomic convergence in extremophile fishes

The molecular basis of convergent phenotypes is often unknown. However, convergence at a genomic level is predicted when there are large population sizes, gene flow among diverging lineages or strong genetic constraints. We used whole-genome resequencing to investigate genomic convergence in fishes ( Poecilia spp.) that have repeatedly colonized hydrogen sulfide (H₂S)-rich environments in Mexico. We identified genomic similarities in both single nucleotide polymorphisms (SNPs) and structural variants (SVs) among independently derived sulfide spring populations, with approximately 1.2% of the genome being shared among sulfidic ecotypes. We compared these convergent genomic regions to candidate genes for H₂S adaptation identified from transcriptomic analyses and found that a significant proportion of these candidate genes (8%) were also in regions where sulfidic individuals had similar SNPs, while only 1.7% were in regions where sulfidic individuals had similar SVs. Those candidate genes included genes involved in sulfide detoxification, the electron transport chain (the main toxicity target of H₂S) and other processes putatively important for adaptation to sulfidic environments. Regional genomic similarity across independent populations exposed to the same source of selection is consistent with selection on standing variation or introgression of adaptive alleles across divergent lineages. However, combined with previous analyses, our data also support that adaptive changes in mitochondrially encoded subunits arose independently via selection on de novo mutations. Pressing questions remain on what conditions ultimately facilitate the independent rise of adaptive alleles at the same loci in separate populations, and thus, the degree to which evolution is repeatable or predictable. This article is part of the theme issue 'Convergent evolution in the genomics era: new insights and directions'.

Concordant Changes in Gene Expression and Nucleotides Underlie Independent Adaptation to Hydrogen-Sulfide-Rich Environments

The colonization of novel environments often involves changes in gene expression, protein coding sequence, or both. Studies of how populations adapt to novel conditions, however, often focus on only one of these two processes, potentially missing out on the relative importance of different parts of the evolutionary process. In this study, our objectives were 1) to better understand the qualitative concordance between conclusions drawn from analyses of differential expression and changes in genic sequence and 2) to quantitatively test whether differentially expressed genes were enriched for sites putatively under positive selection within gene regions. To achieve this, we compared populations of fish (Poecilia mexicana) that have independently adapted to hydrogen-sulfide-rich environments in southern Mexico to adjacent populations residing in nonsulfidic waters. Specifically, we used RNA-sequencing data to compare both gene expression and DNA sequence differences between populations. Analyzing these two different data types led to similar conclusions about which biochemical pathways (sulfide detoxification and cellular respiration) were involved in adaptation to sulfidic environments. Additionally, we found a greater overlap between genes putatively under selection and differentially expressed genes than expected by chance. We conclude that considering both differential expression and changes in DNA sequence led to a more comprehensive understanding of how these populations adapted to extreme environmental conditions. Our results imply that changes in both gene expression and DNA sequence-sometimes at the same loci-may be involved in adaptation.

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.

Brain differences in ecologically differentiated sticklebacks

Populations that have recently diverged offer a powerful model for studying evolution. Ecological differences are expected to generate divergent selection on multiple traits, including neurobiological ones. Animals must detect, process, and act on information from their surroundings and the form of this information can be highly dependent on the environment. We might expect different environments to generate divergent selection not only on the sensory organs, but also on the brain regions responsible for processing sensory information. Here, we test this hypothesis using recently evolved reproductively isolated species pairs of threespine stickleback fish Gasterosteus aculeatus that have well-described differences in many morphological and behavioral traits correlating with ecological differences. We use a state-of-the-art method, magnetic resonance imaging, to get accurate volumetric data for 2 sensory processing regions, the olfactory bulbs and optic tecta. We found a tight correlation between ecology and the size of these brain regions relative to total brain size in 2 lakes with intact species pairs. Limnetic fish, which rely heavily on vision, had relatively larger optic tecta and smaller olfactory bulbs compared with benthic fish, which utilize olfaction to a greater extent. Benthic fish also had larger total brain volumes relative to their body size compared with limnetic fish. These differences were erased in a collapsed species pair in Enos Lake where anthropogenic disturbance has led to intense hybridization. Together these data indicate that evolution of sensory processing regions can occur rapidly and independently.

Convergent evolution of reduced energy demands in extremophile fish

Convergent evolution in organismal function can arise from nonconvergent changes in traits that contribute to that function. Theory predicts that low resource availability and high maintenance costs in extreme environments select for reductions in organismal energy demands, which could be attained through modifications of body size or metabolic rate. We tested for convergence in energy demands and underlying traits by investigating livebearing fish (genus Poecilia) that have repeatedly colonized toxic, hydrogen sulphide-rich springs. We quantified variation in body size and routine metabolism across replicated sulphidic and non-sulphidic populations in nature, modelled total organismal energy demands, and conducted a common-garden experiment to test whether population differences had a genetic basis. Sulphidic populations generally exhibited smaller body sizes and lower routine metabolic rates compared to non-sulphidic populations, which together caused significant reductions in total organismal energy demands in extremophile populations. Although both mechanisms contributed to variation in organismal energy demands, variance partitioning indicated reductions of body size overall had a greater effect than reductions of routine metabolism. Finally, population differences in routine metabolism documented in natural populations were maintained in common-garden reared individuals, indicating evolved differences. In combination with other studies, these results suggest that reductions in energy demands may represent a common theme in adaptation to physiochemical stressors. Selection for reduced energy demand may particularly affect body size, which has implications for life history evolution in extreme environments.

The roles of plasticity and evolutionary change in shaping gene expression variation in natural populations of extremophile fish

The notorious plasticity of gene expression responses and the complexity of environmental gradients complicate the identification of adaptive differences in gene regulation among populations. We combined transcriptome analyses in nature with common-garden and exposure experiments to establish cause-effect relationships between the presence of a physiochemical stressor and expression differences, as well as to test how evolutionary change and plasticity interact to shape gene expression variation in natural systems. We studied two evolutionarily independent population pairs of an extremophile fish (Poecilia mexicana) living in toxic, hydrogen sulphide (H₂ S)-rich springs and adjacent nontoxic habitats and assessed genomewide expression patterns of wild-caught and common-garden-raised individuals exposed to different concentrations of H₂ S. We found that 7.7% of genes that were differentially expressed between sulphidic and nonsulphidic ecotypes remained differentially expressed in the laboratory, indicating that sources of selection other than H₂ S-or plastic responses to other environmental factors-contribute substantially to gene expression patterns observed in the wild. Concordantly differentially expressed genes in the wild and the laboratory were primarily associated with H₂ S detoxification, sulphur processing and metabolic physiology. While shared, ancestral plasticity played a minor role in shaping gene expression variation observed in nature, we documented evidence for evolved population differences in the constitutive expression as well as the H₂ S inducibility of candidate genes. Mechanisms underlying gene expression variation also varied substantially across the two ecotype pairs. These results provide a springboard for studying evolutionary modifications of gene regulatory mechanisms that underlie expression variation in locally adapted populations.

Genome-scale data reveal that endemic Poecilia populations from small sulphidic springs display no evidence of inbreeding

Populations with limited ranges can be highly vulnerable to changes in their environment and are, thus, of high conservation concern. Populations that experience human-induced range reductions are often highly inbred and lack genetic diversity, but it is unknown whether this is also the case for populations with naturally small ranges. The fishes Poecilia sulphuraria (listed as critically endangered) and Poecilia thermalis, which are endemic to small hydrogen sulphide-rich springs in southern Mexico, are examples of such populations with inherently small habitats. We used geometric morphometrics and population genetics to quantify phenotypic and genetic variation within and among two populations of P. sulphuraria and one population of P. thermalis. Principal component analyses revealed phenotypic and genetic differences among the populations. Evidence for inbreeding was low compared to populations that have undergone habitat reduction. The genetic data were also used to infer the demographic history of these populations to obtain estimates for effective population sizes and migration rates. Effective population sizes were large given the small habitats of these populations. Our results imply that these three endemic extremophile populations should each be considered separately for conservation purposes. Additionally, this study suggests that populations in naturally small habitats may have lower rates of inbreeding and higher genetic diversity than expected, and therefore may be better equipped to handle environmental perturbations than anticipated. We caution, however, that the inferred lack of inbreeding and the large effective population sizes could potentially be a result of colonization by genetically diverse ancestors.

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.

Selection is stronger in early-versus-late stages of divergence in a Neotropical livebearing fish

How selection acts to drive trait evolution at different stages of divergence is of fundamental importance in our understanding of the origins of biodiversity. Yet, most studies have focused on a single point along an evolutionary trajectory. Here, we provide a case study evaluating the strength of divergent selection acting on life-history traits at early-versus-late stages of divergence in Brachyrhaphis fishes. We find that the difference in selection is stronger in the early-diverged population than the late-diverged population, and that trait differences acquired early are maintained over time.

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.

Does personality affect premating isolation between locally-adapted populations?

Background

One aspect of premating isolation between diverging, locally-adapted population pairs is female mate choice for resident over alien male phenotypes. Mating preferences often show considerable individual variation, and whether or not certain individuals are more likely to contribute to population interbreeding remains to be studied. In the Poecilia mexicana-species complex different ecotypes have adapted to hydrogen sulfide (H₂S)-toxic springs, and females from adjacent non-sulfidic habitats prefer resident over sulfide-adapted males. We asked if consistent individual differences in behavioral tendencies (animal personality) predict the strength and direction of the mate choice component of premating isolation in this system.

Results

We characterized focal females for their personality and found behavioral measures of ‘novel object exploration’, ‘boldness’ and ‘activity in an unknown area’ to be highly repeatable. Furthermore, the interaction term between our measures of exploration and boldness affected focal females’ strength of preference (SOP) for the resident male phenotype in dichotomous association preference tests. High exploration tendencies were coupled with stronger SOPs for resident over alien mating partners in bold, but not shy, females. Shy and/or little explorative females had an increased likelihood of preferring the non-resident phenotype and thus, are more likely to contribute to rare population hybridization. When we offered large vs. small conspecific stimulus males instead, less explorative females showed stronger preferences for large male body size. However, this effect disappeared when the size difference between the stimulus males was small.

Conclusions

Our results suggest that personality affects female mate choice in a very nuanced fashion. Hence, population differences in the distribution of personality types could be facilitating or impeding reproductive isolation between diverging populations depending on the study system and the male trait(s) upon which females base their mating decisions, respectively.

Electronic supplementary material

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

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.

Unique evolutionary trajectories in repeated adaptation to hydrogen sulphide-toxic habitats of a neotropical fish (Poecilia mexicana)

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.

Mechanisms Underlying Adaptation to Life in Hydrogen Sulfide-Rich 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.

Context dependence of maternal effects: testing assumptions of optimal egg size, differential, and sex allocation models

If offspring develop in adverse conditions, the maternal component of their phenotypic variation might increase due to the stronger dependence of offspring traits on parental investment. This should result in increased parental investment to individual offspring, as assumed by the model of optimal egg size. The opposite pattern, i.e., stronger dependence of offspring fitness on parental investment and consequently larger parental investment under good conditions is assumed by both the theory of differential allocation if attractive males provide material benefits, and reproductive compensation if they invest less into paternal care. Another influential idea is the Trivers-Willard model, which assumes sex-specific dependence of offspring fitness on parental investment. Here we tested these ideas by examining the effects of egg size on offspring fitness across many postnatal contexts in the Collared Flycatcher Ficedula albicollis. We employed a cross-fostering design that generated variation in egg size within nests and used brood means of fledgling mass as a functional measure of the quality of rearing conditions. Effects of egg size on three offspring traits, including lifetime reproductive success of recruits, were more pronounced in low-quality broods. These results support the assumption of the model of optimal egg size. Based on female preference for males providing material benefits, this pattern could support differential allocation, if attractive males invest less in paternal care, or reproductive compensation, if they invest more. By comparison, we did not find any evidence for sex specificity of fitness returns that might explain sex monomorphism of egg size in this species. The challenge for future studies will be the integration of components of parental investment and offspring fitness into their global measures and testing how the former affects the latter across gradients of postnatal conditions.

Speciation in Freshwater Fishes

The extraordinary species richness of freshwater fishes has attracted much research on mechanisms and modes of speciation. We here review research on speciation in freshwater fishes in light of speciation theory, and place this in a context of broad-scale diversity patterns in freshwater fishes. We discuss several major repeated themes in freshwater fish speciation and the speciation mechanisms they are frequently associated with. These include transitions between marine and freshwater habitats, transitions between discrete freshwater habitats, and ecological transitions within habitats, as well as speciation without distinct niche shifts. Major research directions in the years to come include understanding the transition from extrinsic environment-dependent to intrinsic reproductive isolation and its influences on species persistence and understanding the extrinsic and intrinsic constraints to speciation and how these relate to broad-scale diversification patterns through time.

Patterns of Macroinvertebrate and Fish Diversity in Freshwater Sulphide Springs

Extreme environments are characterised by the presence of physicochemical stressors and provide unique study systems to address problems in evolutionary ecology research. Sulphide springs provide an example of extreme freshwater environments; because hydrogen sulphide’s adverse physiological effects induce mortality in metazoans even at micromolar concentrations. Sulphide springs occur worldwide, but while microbial communities in sulphide springs have received broad attention, little is known about macroinvertebrates and fish inhabiting these toxic environments. We reviewed qualitative occurrence records of sulphide spring faunas on a global scale and present a quantitative case study comparing diversity patterns in sulphidic and adjacent non-sulphidic habitats across replicated river drainages in Southern Mexico. While detailed studies in most regions of the world remain scarce, available data suggests that sulphide spring faunas are characterised by low species richness. Dipterans (among macroinvertebrates) and cyprinodontiforms (among fishes) appear to dominate the communities in these habitats. At least in fish, there is evidence for the presence of highly endemic species and populations exclusively inhabiting sulphide springs. We provide a detailed discussion of traits that might predispose certain taxonomic groups to colonize sulphide springs, how colonizers subsequently adapt to cope with sulphide toxicity, and how adaptation may be linked to speciation processes.

Microhabitat use, population densities, and size distributions of sulfur cave-dwelling Poecilia mexicana

The Cueva del Azufre in Tabasco, Mexico, is a nutrient-rich cave and its inhabitants need to cope with high levels of dissolved hydrogen sulfide and extreme hypoxia. One of the successful colonizers of this cave is the poeciliid fish Poecilia mexicana, which has received considerable attention as a model organism to examine evolutionary adaptations to extreme environmental conditions. Nonetheless, basic ecological data on the endemic cave molly population are still missing; here we aim to provide data on population densities, size class compositions and use of different microhabitats. We found high overall densities in the cave and highest densities at the middle part of the cave with more than 200 individuals per square meter. These sites have lower H₂S concentrations compared to the inner parts where most large sulfide sources are located, but they are annually exposed to a religious harvesting ceremony of local Zoque people called La Pesca. We found a marked shift in size/age compositions towards an overabundance of smaller, juvenile fish at those sites. We discuss these findings in relation to several environmental gradients within the cave (i.e., differences in toxicity and lighting conditions), but we also tentatively argue that the annual fish harvest during a religious ceremony (La Pesca) locally diminishes competition (and possibly, cannibalism by large adults), which is followed by a phase of overcompensation of fish densities.

H2S exposure elicits differential expression of candidate genes in fish adapted to sulfidic and non-sulfidic environments

Disentangling the effects of plasticity, genetic variation, and their interactions on organismal responses to environmental stressors is a key objective in ecological physiology. We quantified the expression of five candidate genes in response to hydrogen sulfide (H2S) exposure in fish (Poecilia mexicana, Poeciliidae) from a naturally sulfide-rich environment as well as an ancestral, non-sulfidic population to test for constitutive and environmentally dependent population differences in gene expression patterns. Common garden raised individuals that had never encountered environmental H2S during their lifetime were subjected to short or long term H2S exposure treatments or respective non-sulfidic controls. The expression of genes involved in responses to H2S toxicity (cytochrome c oxidase, vascular endothelial growth factor, and cytochrome P450-2J6), H2S detoxification (sulfide:quinone oxidoreductase), and endogenous H2S production (cystathionine γ lyase) was determined in both gill and liver tissues by real time PCR. The results indicated complex changes in expression patterns that--depending on the gene--not only differed between organs and populations, but also on the type of H2S exposure. Populations differences, both constitutive and H2S exposure dependent (i.e., plastic), in gene expression were particularly evident for sulfide:quinone oxidoreductase, vascular endothelial growth factor, and to a lesser degree for cytochrome P450-2J6. Our study uncovered putatively adaptive modifications in gene regulation that parallel previously documented adaptive changes in phenotypic traits.