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.

Why and when should organisms elongate their telomeres? Elaborations on the 'excess resources elongation' and 'last resort elongation' hypotheses.

Telomere length and telomere shortening are thought to be critical cellular attributes and processes that are related to an individual's life span and fitness. The general pattern across most taxa is that after birth telomere length gradually decreases with age. Telomere protection and restoration mechanisms are usually assumed to reduce the rate of shortening or at most keep telomere length constant. However, here we have compiled a list of 26 articles showing that there is an increasing number of studies reporting apparent elongation of telomeres (i.e., a net increase in TL from timet to timet+1) often in a considerable proportion of the individuals studied. Moreover, the few studies which have studied telomere elongation in detail show that increases in telomere length are unlikely to be due to measurement error alone. In this article, we argue that episodes of telomere elongation deserve more attention as they could reflect individual strategies to optimise life histories and maximise fitness, which may not be reflected in the overall telomere dynamics patterns. We propose that patterns of telomere (net) elongation may be partly determined by other factors than those causing telomere shortening, and therefore deserve analyses specifically targeted to investigate the occurrence of telomere elongation. We elaborate on two ecological hypotheses that have been proposed to explain patterns of telomere elongation (the 'excess resources elongation' and the 'last resort elongation' hypothesis) and we discuss the current evidence for (or against) these hypotheses and propose ways to test them.

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.

Assessment of avian health status: suitability and constraints of the Zoetis VetScan VS2 blood analyser for ecological and evolutionary studies.

Biochemical analyses of blood can decipher physiological conditions of living animals and unravel mechanistic underpinnings of life-history strategies and trade-offs. Yet, researchers in ecology and evolution often face constraints in which methods to apply, not least due to blood volume restrictions or field settings. Here, we test the suitability of a portable biochemical analyser (Zoetis VetScan VS2) for ecological and evolutionary studies that may help solve those problems. Using as little as 80 µl of whole-bird blood from free-living Jackdaws (Corvus monedula) and captive Zebra Finches (Taeniopygia guttata), we show that eight (out of 10) blood analytes show high repeatability after short-term storage (approximately 2 h) and six after 12 h storage time. Handling stress had a clear impact on all except two analytes by 16 min after catching. Finally, six analytes showed consistency within individuals over a period of 30 days, and three even showed individual consistency over a year. Taken together, we conclude that the VetScan VS2 captures biologically relevant variation in blood analytes using just 80 µl of whole blood and, thus, provides valuable physiological measurements of (small) birds sampled in semi-field and field conditions.

Impacts of heavy metal pollution on the ionomes and transcriptomes of Western mosquitofish (Gambusia affinis).

Our understanding of the mechanisms mediating the resilience of organisms to environmental change remains lacking. Heavy metals negatively affect processes at all biological scales, yet organisms inhabiting contaminated environments must maintain homeostasis to survive. Tar Creek in Oklahoma, USA, contains high concentrations of heavy metals and an abundance of Western mosquitofish (Gambusia affinis), though several fish species persist at lower frequency. To test hypotheses about the mechanisms mediating the persistence and abundance of mosquitofish in Tar Creek, we integrated ionomic data from seven resident fish species and transcriptomic data from mosquitofish. We predicted that mosquitofish minimize uptake of heavy metals more than other Tar Creek fish inhabitants and induce transcriptional responses to detoxify metals that enter the body, allowing them to persist in Tar Creek at higher density than species that may lack these responses. Tar Creek populations of all seven fish species accumulated heavy metals, suggesting mosquitofish cannot block uptake more efficiently than other species. We found population-level gene expression changes between mosquitofish in Tar Creek and nearby unpolluted sites. Gene expression differences primarily occurred in the gill, where we found upregulation of genes involved with lowering transfer of metal ions from the blood into cells and mitigating free radicals. However, many differentially expressed genes were not in known metal response pathways, suggesting multifarious selective regimes and/or previously undocumented pathways could impact tolerance in mosquitofish. Our systems-level study identified well characterized and putatively new mechanisms that enable mosquitofish to inhabit heavy metal-contaminated environments.

Telomeres in ecology and evolution: A review and classification of hypotheses.

Research on telomeres in the fields of ecology and evolution has been rapidly expanding over the last two decades. This has resulted in the formulation of a multitude of, often name-given, hypotheses related to the associations between telomeres and life-history traits or fitness-facilitating processes (and the mechanisms underlying them). However, the differences (or similarities) between the various hypotheses, which can originate from different research fields, are often not obvious. Our aim here is therefore to give an overview of the hypotheses that are of interest in ecology and evolution and to provide two frameworks that help discriminate among them. We group the hypotheses (i) based on their association with different research questions, and (ii) using a hierarchical approach that builds on the assumptions they make, such as about causality of telomere length/shortening and/or the proposed functional consequences of telomere shortening on organism performance. Both our frameworks show that there exist parallel lines of thoughts in different research fields. Moreover, they also clearly illustrate that there are in many cases competing hypotheses within clusters, and that some of these even have contradictory assumptions and/or predictions. We also touch upon two topics in telomere research that would benefit from further conceptualization. This review should help researchers, both those familiar with and those new to the subject, to identify future avenues of research.

Parallel shifts of visual sensitivity and body coloration in replicate populations of extremophile fish.

Visual sensitivity and body pigmentation are often shaped by both natural selection from the environment and sexual selection from mate choice. One way of quantifying the impact of the environment is by measuring how traits have changed after colonization of a novel habitat. To do this, we studied Poecilia mexicana populations that have repeatedly adapted to extreme sulphidic (H2 S-containing) environments. We measured visual sensitivity using opsin gene expression, as well as body pigmentation, for populations in four independent drainages. Both visual sensitivity and body pigmentation showed significant parallel shifts towards greater medium-wavelength sensitivity and reflectance in sulphidic populations. Altogether we found that sulphidic habitats select for differences in visual sensitivity and pigmentation. Shifts between habitats may be due to both differences in the water's spectral properties and correlated ecological changes.

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.

microRNA expression variation as a potential molecular mechanism contributing to adaptation to hydrogen sulphide.

microRNAs (miRNAs) are post-transcriptional regulators of gene expression and can play an important role in modulating organismal development and physiology in response to environmental stress. However, the role of miRNAs in mediating adaptation to diverse environments in natural study systems remains largely unexplored. Here, we characterized miRNAs and their expression in Poecilia mexicana, a species of small fish that inhabits both normal streams and extreme environments in the form of springs rich in toxic hydrogen sulphide (H2 S). We found that P. mexicana has a similar number of miRNA genes as other teleosts. In addition, we identified a large population of mature miRNAs that were differentially expressed between locally adapted populations in contrasting habitats, indicating that miRNAs may contribute to P. mexicana adaptation to sulphidic environments. In silico identification of differentially expressed miRNA-mRNA pairs revealed, in the sulphidic environment, the downregulation of miRNAs predicted to target mRNAs involved in sulphide detoxification and cellular homeostasis, which are pathways essential for life in H2 S-rich springs. In addition, we found that predicted targets of upregulated miRNAs act in the mitochondria (16.6% of predicted annotated targets), which is the main site of H2 S toxicity and detoxification, possibly modulating mitochondrial function. Together, the differential regulation of miRNAs between these natural populations suggests that miRNAs may be involved in H2 S adaptation by promoting functions needed for survival and reducing functions affected by H2 S. This study lays the groundwork for further research to directly demonstrate the role of miRNAs in adaptation to H2 S. Overall, this study provides a critical stepping-stone towards a comprehensive understanding of the regulatory mechanisms underlying the adaptive variation in gene expression in a natural system.

Convergent evolution of conserved mitochondrial pathways underlies repeated adaptation to extreme environments.

Extreme environments test the limits of life; yet, some organisms thrive in harsh conditions. Extremophile lineages inspire questions about how organisms can tolerate physiochemical stressors and whether the repeated colonization of extreme environments is facilitated by predictable and repeatable evolutionary innovations. We identified the mechanistic basis underlying convergent evolution of tolerance to hydrogen sulfide (H2S)-a toxicant that impairs mitochondrial function-across evolutionarily independent lineages of a fish (Poecilia mexicana, Poeciliidae) from H2S-rich springs. Using comparative biochemical and physiological analyses, we found that mitochondrial function is maintained in the presence of H2S in sulfide spring P. mexicana but not ancestral lineages from nonsulfidic habitats due to convergent adaptations in the primary toxicity target and a major detoxification enzyme. Genome-wide local ancestry analyses indicated that convergent evolution of increased H2S tolerance in different populations is likely caused by a combination of selection on standing genetic variation and de novo mutations. On a macroevolutionary scale, H2S tolerance in 10 independent lineages of sulfide spring fishes across multiple genera of Poeciliidae is correlated with the convergent modification and expression changes in genes associated with H2S toxicity and detoxification. Our results demonstrate that the modification of highly conserved physiological pathways associated with essential mitochondrial processes mediates tolerance to physiochemical stress. In addition, the same pathways, genes, and-in some instances-codons are implicated in H2S adaptation in lineages that span 40 million years of evolution.

Functional consequences of phenotypic variation between locally adapted populations: Swimming performance and ventilation in extremophile fish.

Natural selection drives the evolution of traits to optimize organismal performance, but optimization of one aspect of performance can influence other aspects of performance. Here, we asked how phenotypic variation between locally adapted fish populations affects locomotion and ventilation, testing for functional trade-offs and trait-performance correlations. Specifically, we investigated two populations of livebearing fish (Poecilia mexicana) that inhabit distinct habitat types (hydrogen-sulphide-rich springs and adjacent nonsulphidic streams). For each individual, we quantified different metrics of burst swimming during simulated predator attacks, steady swimming and gill ventilation. Coinciding with predictions, we documented significant population differences in all aspects of performance, with fish from sulphidic habitats exhibiting higher steady swimming performance and higher ventilation capacity, but slower burst swimming. There was a significant functional trade-off between steady and burst swimming, but not between different aspects of locomotion and ventilation. Although our findings about population differences in locomotion performance largely parallel the results from previous studies, we provide novel insights about how morphological variation might impact ventilation and ultimately oxygen acquisition. Overall, our analyses provided insights into the functional consequences of previously documented phenotypic variation, which will help to disentangle the effects of different sources of selection that may coincide along complex environmental gradients.

Temperature effects on performance and physiology of two prairie stream minnows.

Earth's atmosphere has warmed by ~1°C over the past century and continues to warm at an increasing rate. Effects of atmospheric warming are already visible in most major ecosystems and are evident across all levels of biological organization. Linking functional responses of individuals to temperature is critical for predicting responses of populations and communities to global climate change. The southern redbelly dace Chrosomus erythrogaster and the central stoneroller Campostoma anomalum are two minnows (Cyprinidae) that commonly occur in the Flint Hills region of the USA but show different patterns of occurrence, with dace largely occupying headwater reaches and stonerollers persisting in both headwater and intermediate-sized streams. We tested for differences between species in critical thermal maximum, energy metabolism, sustained swimming and activity over an ecologically relevant temperature gradient of acclimation temperatures. Typically, metrics increased with acclimation temperature for both species, although stoneroller activity decreased with temperature. We observed a significant interaction between species and temperature for critical thermal maxima, where stonerollers only had higher critical thermal maxima at the coldest temperature and at warm temperatures compared to the dace. We did not find evidence suggesting differences in the energy metabolism of dace and stonerollers. We detected interspecific differences in sustained swimming performance, with dace having higher swimming speed than stonerollers regardless of acclimation temperature. Finally, there was a significant interaction between temperature and species for activity; dace activity was higher at intermediate and warm temperatures compared to stonerollers. We observed subtle interspecific differences in how performance metrics responded to temperature that did not always align with observed patterns of distribution for these species. Thus, other ecological factors likely are important drivers of distributional patterns in these species.

Expression analyses of cave mollies (Poecilia mexicana) reveal key genes involved in the early evolution of eye regression.

Eye regression occurs across cave-dwelling populations of many species and is often coupled with a decrease or loss in eye function. Teleost fishes are among the few vertebrates to undergo widespread colonization of caves and often exhibit eye regression with blindness. Cave populations of the poeciliid fish Poecilia mexicana (cave molly) exhibit reduced-albeit functional-eyes, offering the opportunity to investigate partial eye regression. We sequenced eye transcriptomes of cave and surface populations of P. mexicana to identify differentially expressed genes that potentially underlie eye regression in cave mollies. We identified 28 significantly differentially expressed genes, 20 of which were directly related to light sensitivity, eye structure and visual signaling. Twenty-six of these genes were downregulated in cave compared to surface populations. Functional enrichment analysis revealed eye-related gene ontologies that were under-represented in cave mollies. In addition, a set of co-expressed genes related to vision and circadian rhythm was correlated with habitat type (cave versus surface). Our study suggests that differential gene expression plays a key role in the beginning evolutionary stages of eye regression in P. mexicana, shedding further light on regressive evolution in cavefish.

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 (H2S)-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 H2S 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 H2S) 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'.

Detection of changes in mitochondrial hydrogen sulfide i n vivo in the fish model Poecilia mexicana (Poeciliidae).

In this paper, we outline the use of a mitochondria-targeted ratiometric mass spectrometry probe, MitoA, to detect in vivo changes in mitochondrial hydrogen sulfide (H2S) in Poecilia mexicana (family Poeciliidae). MitoA is introduced via intraperitoneal injection into the animal and is taken up by mitochondria, where it reacts with H2S to form the product MitoN. The MitoN/MitoA ratio can be used to assess relative changes in the amounts of mitochondrial H2S produced over time. We describe the use of MitoA in the fish species P. mexicana to illustrate the steps for adopting the use of MitoA in a new organism, including extraction and purification of MitoA and MitoN from tissues followed by tandem mass spectrometry. In this proof-of-concept study we exposed H2S tolerant P. mexicana to 59 µM free H2S for 5 h, which resulted in increased MitoN/MitoA in brain and gills, but not in liver or muscle, demonstrating increased mitochondrial H2S levels in select tissues following whole-animal H2S exposure. This is the first time that accumulation of H2S has been observed in vivo during whole-animal exposure to free H2S using MitoA. This article has an associated First Person interview with the first author of the paper.

Correlated divergence of female and male genitalia in replicated lineages with ongoing ecological speciation.

Divergence of genital traits among lineages has the potential to serve as a reproductive isolating barrier when copulation, insemination, and fertilization are inhibited by incompatibilities between female and male genitalia. Despite widespread evidence for genital trait diversity among closely related lineages and coevolution of female and male genitalia within lineages, few studies have investigated genital evolution during the early stages of speciation. We quantified genital variation in replicated population pairs of Poecilia mexicana with ongoing ecological speciation between sulfidic (H2 S containing) and nearby nonsulfidic habitats. These analyses revealed rapid and correlated divergence of female and male genitalia across evolutionarily independent population pairs exposed to divergent selection regimes. Both sexes exhibited convergent evolution of genital traits among populations inhabiting similar habitat types. Our results demonstrate that genital evolution can occur during the early stages of speciation-with-gene-flow, potentially as a result of variation in the intensity of sexual conflict among populations. Our results suggest genitalia may contribute to early stages of divergence and challenge the generality of previously suggested mechanisms of genital evolution in poeciliids.

Bacterial Diversity in Replicated Hydrogen Sulfide-Rich Streams.

Extreme environments typically require costly adaptations for survival, an attribute that often translates to an elevated influence of habitat conditions on biotic communities. Microbes, primarily bacteria, are successful colonizers of extreme environments worldwide, yet in many instances, the interplay between harsh conditions, dispersal, and microbial biogeography remains unclear. This lack of clarity is particularly true for habitats where extreme temperature is not the overarching stressor, highlighting a need for studies that focus on the role other primary stressors (e.g., toxicants) play in shaping biogeographic patterns. In this study, we leveraged a naturally paired stream system in southern Mexico to explore how elevated hydrogen sulfide (H2S) influences microbial diversity. We sequenced a portion of the 16S rRNA gene using bacterial primers for water sampled from three geographically proximate pairings of streams with high (> 20 µM) or low (~ 0 µM) H2S concentrations. After exploring bacterial diversity within and among sites, we compared our results to a previous study of macroinvertebrates and fish for the same sites. By spanning multiple organismal groups, we were able to illuminate how H2S may differentially affect biodiversity. The presence of elevated H2S had no effect on overall bacterial diversity (p = 0.21), a large effect on community composition (25.8% of variation explained, p < 0.0001), and variable influence depending upon the group-whether fish, macroinvertebrates, or bacteria-being considered. For bacterial diversity, we recovered nine abundant operational taxonomic units (OTUs) that comprised a core H2S-rich stream microbiome in the region. Many H2S-associated OTUs were members of the Epsilonproteobacteria and Gammaproteobacteria, which both have been implicated in endosymbiotic relationships between sulfur-oxidizing bacteria and eukaryotes, suggesting the potential for symbioses that remain to be discovered in these habitats.

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.

Publisher Correction: Sex-specific evolution during the diversification of live-bearing fishes.

In the version of this Article originally published, some production notes starting "Should we change..." were mistakenly left in at the end of the section 'Sexual selection'; these notes have now been removed.

Correlated evolution of thermal niches and functional physiology in tropical freshwater fishes.

The role of ecology in phenotypic and species diversification is widely documented. Nonetheless, numerous nonadaptive processes can shape realized niches and phenotypic variation in natural populations, complicating inferences about adaptive evolution at macroevolutionary scales. We tested for evolved differences in thermal tolerances and their association with the realized thermal niche (including metrics describing diurnal and seasonal patterns of temperature extremes and variability) across a genus of tropical freshwater fishes reared in a standardized environment. There was limited evolution along the thermal niche axis associated with variation in maximum temperature and in upper thermal limits. In contrast, there was considerable diversification along the first major axis of the thermal niche associated with minimum temperatures and in lower thermal limits. Across our adaptive landscape analyses, 70% of species exhibited evidence of divergence in thermal niches. Most importantly, the first two major axes of thermal niche variation were significantly correlated with variation in lower thermal limits. Our results indicate adaptation to divergent thermal niches and adaptive evolution of related functional traits, and highlight the importance of divergence in lower thermal limits for the evolution of tropical biodiversity.