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.

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.

Microbial mechanisms of C/N/S geochemical cycling during low-water-level sediment remediation in urban rivers

Low-water-level regulation has been effectively implemented in the restoration of urban river sediments in Guangzhou City, China. Further investigation is needed to understand the microbial mechanisms involved in pollutant degradation in low-water-level environments. This study examined sediment samples from nine rivers, including low-water-level rivers (LW), tidal waterways (TW), and enclosed rivers (ER). Metagenomic high-throughput sequencing and the Diting pipeline were utilized to investigate the microbial mechanisms involved in sediment C/N/S geochemical cycling during low-water-level regulation. The results reveal that the degree of pollution in LW sediment is lower compared to TW and ER sediment. LW sediment exhibits a higher capacity for pollutant degradation and elimination of black, odorous substances due to its stronger microbial methane oxidation, nitrification, denitrification, anammox, and oxidation of sulfide, sulfite, and thiosulfate. Conversely, TW and ER sediment showcase greater microbial methanogenesis, anaerobic fermentation, and sulfide generation abilities, leading to the persistence of black, odorous substances. Factors such as grit and silt content, nitrate, and ammonia concentrations impacted microbial metabolic pathways. Low-water-level regulation improved the micro-environment for functional microbes, facilitating pollutant removal and preventing black odorous substance accumulation. These findings provide insights into the microbial mechanisms underlying low-water-level regulation technology for sediment restoration in urban rivers.

Functional convergence in gastric lysozymes of foregut-fermenting rodents, ruminants, and primates is not attributed to convergent molecular evolution

Convergent evolution is a widespread phenomenon. While there are many examples of convergent evolution at the phenotypic scale, convergence at the molecular level has been more difficult to identify. A classic example of convergent evolution across scales is that of the digestive lysozyme found in ruminants and Colobine monkeys. These herbivorous species rely on foregut fermentation, which has evolved to function more optimally under acidic conditions. Here, we explored if rodents with similar dietary strategies and digestive morphologies have convergently evolved a lysozyme with digestive functions. At the phenotypic level, we find that rodents with bilocular stomach morphologies exhibited a lysozyme that maintained higher relative activities at low pH values, similar to the lysozymes of ruminants and Colobine monkeys. Additionally, the lysozyme of Peromyscus leucopus shared a similar predicted protonation state as that observed in previously identified digestive lysozymes. However, we found limited evidence of positive selection acting on the lysozyme gene in foregut-fermenting species and did not identify patterns of convergent molecular evolution in this gene. This study emphasizes that phenotypic convergence need not be the result of convergent genetic modifications, and we encourage further exploration into the mechanisms regulating convergence across biological 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.

Cold-Driven Hemoglobin Evolution in Antarctic Notothenioid Fishes Prior to Hemoglobin Gene Loss in White-Blooded Icefishes

Expression of multiple hemoglobin isoforms with differing physiochemical properties likely helps species adapt to different environmental and physiological conditions. Antarctic notothenioid fishes inhabit the icy Southern Ocean and display fewer hemoglobin isoforms, each with less affinity for oxygen than temperate relatives. Reduced hemoglobin multiplicity was proposed to result from relaxed selective pressure in the cold, thermally stable, and highly oxygenated Antarctic waters. These conditions also permitted the survival and diversification of white-blooded icefishes, the only vertebrates living without hemoglobin. To understand hemoglobin evolution during adaptation to freezing water, we analyzed hemoglobin genes from 36 notothenioid genome assemblies. Results showed that adaptation to frigid conditions shaped hemoglobin gene evolution by episodic diversifying selection concomitant with cold adaptation and by pervasive evolution in Antarctic notothenioids compared to temperate relatives, likely a continuing adaptation to Antarctic conditions. Analysis of hemoglobin gene expression in adult hematopoietic organs in various temperate and Antarctic species further revealed a switch in hemoglobin gene expression underlying hemoglobin multiplicity reduction in Antarctic fish, leading to a single hemoglobin isoform in adult plunderfishes and dragonfishes, the sister groups to icefishes. The predicted high hemoglobin multiplicity in Antarctic fish embryos based on transcriptomic data, however, raises questions about the molecular bases and physiological implications of diverse hemoglobin isoforms in embryos compared to adults. This analysis supports the hypothesis that the last common icefish ancestor was vulnerable to detrimental mutations affecting the single ancestral expressed alpha- and beta-globin gene pair, potentially predisposing their subsequent loss.

Genomic data support the taxonomic validity of Middle American livebearers Poeciliopsis gracilis and Poeciliopsis pleurospilus (Cyprinodontiformes: Poeciliidae)

Poeciliopsis (Cyprinodontiformes: Poeciliidae) is a genus comprised of 25 species of freshwater fishes. Several well-known taxonomic uncertainties exist within the genus, especially in relation to the taxonomic status of Poeciliopsis pleurospilus and P. gracilis. However, to date, no studies have been conducted to specifically address the taxonomic status of these two species. The goal of this study was to examine the taxonomic validity of P. pleurospilus and P. gracilis using genomic data (ddRADseq) in phylogenetic, population genetic, and species delimitation frameworks. Multiple analyses support the recognition of both taxa as distinct species and also permits us to revise their respective distributions. A species delimitation analysis indicates that P. pleurospilus and P. gracilis are distinct species, each of which consists of two distinct lineages that are geographically structured. Phylogenetic and population genetic analyses provide clear evidence that individuals of P. gracilis are distributed north and west of the Isthmus of Tehuantepec in both Pacific and Atlantic river systems in Mexico, whereas individuals of P. pleurospilus are distributed in both Atlantic and Pacific river systems south and east of the Isthmus of Tehuantepec, from southern Mexico to Honduras.

The use of non-lethal sampling for transcriptomics to assess the physiological status of wild fishes

Fishes respond to different abiotic and biotic stressors through changes in gene expression as a part of an integrated physiological response. Transcriptomics approaches have been used to quantify gene expression patterns as a reductionist approach to understand responses to environmental stressors in animal physiology and have become more commonly used to study wild fishes. We argue that non-lethal sampling for transcriptomics should become the norm for assessing the physiological status of wild fishes, especially when there are conservation implications. Processes at the level of the transcriptome provide a "snapshot" of the cellular conditions at a given time; however, by using a non-lethal sampling protocol, researchers can connect the transcriptome profile with fitness-relevant ecological endpoints such as reproduction, movement patterns and survival. Furthermore, telemetry is a widely used approach in fisheries to understand movement patterns in the wild, and when combined with transcriptional profiling, provides arguably the most powerful use of non-lethal sampling for transcriptomics in wild fishes. In this review, we discuss the different tissues that can be successfully incorporated into non-lethal sampling strategies, which is particularly useful in the context of the emerging field of conservation transcriptomics. We briefly describe different methods for transcriptional profiling in fishes from high-throughput qPCR to whole transcriptome approaches. Further, we discuss strategies and the limitations of using transcriptomics for non-lethally studying fishes. Lastly, as 'omics' technology continues to advance, transcriptomics paired with different omics approaches to study wild fishes will provide insight into the factors that regulate phenotypic variation and the physiological responses to changing environmental conditions in the future.

Diurnal Changes in Hypoxia Shape Predator-Prey Interaction in a Bird-Fish System

Animals often face changing environments, and behavioral flexibility allows them to rapidly and adaptively respond to abiotic factors that vary more or less regularly. However, abiotic factors that affect prey species do not necessarily affect their predators. Still, the prey’s response might affect the predator indirectly, yet evidence from the wild for such a classical bottom-up effect of abiotic factors shaping several trophic levels remains sparse. In many aquatic environments, daily changes in oxygen concentrations occur frequently. When oxygen levels drop to hypoxic levels, many fishes respond with aquatic surface respiration (ASR), during which they obtain oxygen by skimming the upper, oxygenated surface layer. By increasing time at the surface, fish become more vulnerable to fish-eating birds. We explored these cascading effects in a sulfidic spring system that harbors the endemic sulphur molly (Poecilia sulphuraria) as prey species and several fish-eating bird species. Sulfide-rich springs pose harsh conditions as hydrogen sulfide (H2S) is lethal to most metazoans and reduces dissolved oxygen (DO). Field sampling during three daytimes indicated that water temperatures rose from morning to (after)noon, resulting in the already low DO levels to decrease further, while H2S levels showed no diurnal changes. The drop in DO levels was associated with a decrease in time spent diving in sulphur mollies, which corresponded with an increase in ASR. Interestingly, the laboratory-estimated threshold at which the majority of sulphur mollies initiate ASR (ASR50: <1.7 mg/L DO) was independent of temperature and this value was exceeded daily when hypoxic stress became more severe toward noon. As fish performed ASR, large aggregations built up at the water surface over the course of the day. As a possible consequence of fish spending more time at the surface, we found high activity levels of fish-eating birds at noon and in the afternoon. Our study reveals that daily fluctuations in water’s oxygen levels have the potential to alter predator-prey interactions profoundly and thus highlights the joined actions of abiotic and biotic factors shaping the evolution of a prey species.

Genomic signals found using RNA sequencing show signatures of selection and subtle population differentiation in walleye (Sander vitreus) in a large freshwater ecosystem

RNA sequencing is an effective approach for studying aquatic species yielding both physiological and genomic data. However, its population genetic applications are not well-characterized. We investigate this possible role for RNA sequencing for population genomics in Lake Winnipeg, Manitoba, Canada, walleye (Sander vitreus). Lake Winnipeg walleye represent the largest component of the second-largest freshwater fishery in Canada. In the present study, large female walleye were sampled via nonlethal gill biopsy over two years at three spawning sites representing a latitudinal gradient in the lake. Genetic variation from sequenced mRNA was analyzed for neutral and adaptive markers to investigate population structure and possible adaptive variation. We find low population divergence (F ST = 0.0095), possible northward gene flow, and outlier loci that vary latitudinally in transcripts associated with cell membrane proteins and cytoskeletal function. These results indicate that Lake Winnipeg walleye may be effectively managed as a single demographically connected metapopulation with contributing subpopulations and suggest genomic differences possibly underlying observed phenotypic differences. Despite its high cost relative to other genotyping methods, RNA sequencing data can yield physiological in addition to genetic information discussed here. We therefore argue that it is useful for addressing diverse molecular questions in the conservation of freshwater species.

Multiplicity and Polymorphism of Fish Hemoglobins

The diversity of fish hemoglobins and the association with oxygen availability and physiological requirements during the life cycle has attracted scientists since the first report on multiple hemoglobin in fishes (Buhler and Shanks 1959). The functional heterogeneity of the fish hemoglobins enables many species to tolerate hypoxic conditions and exhausting swimming, but also to maintain the gas pressure in the swim bladder at large depths. The hemoglobin repertoire has further increased in various species displaying polymorphic hemoglobin variants differing in oxygen binding properties. The multiplicity of fish hemoglobins as particularly found in the tetraploid salmonids strongly contrasts with the complete loss of hemoglobins in Antarctic icefishes and illustrates the adaptive radiation in the oxygen transport of this successful vertebrate group.

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'.

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.