Complexities of gene expression patterns in natural populations of an extremophile fish (Poecilia mexicana, Poeciliidae)

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.

Contrasting transcriptomic patterns reveal a genomic basis for drought resilience in the relict fir Abies pinsapo Boiss

Climate change challenges the adaptive capacity of several forest tree species in the face of increasing drought and rising temperatures. Therefore, understanding the mechanistic connections between genetic diversity and drought resilience is highly valuable for conserving drought-sensitive forests. Nonetheless, the post-drought recovery in trees from a transcriptomic perspective has not yet been studied by comparing contrasting phenotypes. Here, experimental drought treatments, gas-exchange dynamics and transcriptomic analysis (RNA-seq) were performed in the relict and drought-sensitive fir Abies pinsapo Boiss. to identify gene expression differences over immediate (24 h) and extended drought (20 days). Post-drought responses were investigated to define resilient and sensitive phenotypes. Single nucleotide polymorphisms (SNPs) were also studied to characterize the genomic basis of A. pinsapo drought resilience. Weighted gene co-expression network analysis showed an activation of stomatal closing and an inhibition of plant growth-related genes during the immediate drought, consistent with an isohydric dynamic. During the extended drought, transcription factors, as well as cellular damage and homeostasis protection-related genes prevailed. Resilient individuals activate photosynthesis-related genes and inhibit aerial growth-related genes, suggesting a shifting shoot/root biomass allocation to improve water uptake and whole-plant carbon balance. About, 152 fixed SNPs were found between resilient and sensitive seedlings, which were mostly located in RNA-activity-related genes, including epigenetic regulation. Contrasting gene expression and SNPs were found between different post-drought resilience phenotypes for the first time in a forest tree, suggesting a transcriptomic and genomic basis for drought resilience. The obtained drought-related transcriptomic profile and drought-resilience candidate genes may guide conservation programs for this threatened tree species.

Adaption to hydrogen sulfide-rich environments: Strategies for active detoxification in deep-sea symbiotic mussels, Gigantidas platifrons

The deep-sea mussel Gigantidas platifrons is a representative species that relies on nutrition provided by chemoautotrophic endosymbiotic bacteria to survive in both hydrothermal vent and methane seep environments. However, vent and seep habitats have distinct geochemical features, with vents being more harsh than seeps because of abundant toxic chemical substances, particularly hydrogen sulfide (H₂S). Until now, the adaptive strategies of G. platifrons in a heterogeneous environment and their sulfide detoxification mechanisms are still unclear. Herein, we conducted 16S rDNA sequencing and metatranscriptome sequencing of G. platifrons collected from a methane seep at Formosa Ridge in the South China Sea and a hydrothermal vent at Iheya North Knoll in the Mid-Okinawa Trough to provide a model for understanding environmental adaption and sulfide detoxification mechanisms, and a three-day laboratory controlled Na₂S stress experiment to test the transcriptomic responses under sulfide stress. The results revealed the active detoxification of sulfide in G. platifrons gills. First, epibiotic Campylobacterota bacteria were more abundant in vent mussels and contributed to environmental adaptation by active oxidation of extracellular H₂S. Notably, a key sulfide-oxidizing gene, sulfide:quinone oxidoreductase (sqr), derived from the methanotrophic endosymbiont, was significantly upregulated in vent mussels, indicating the oxidization of intracellular sulfide by the endosymbiont. In addition, transcriptomic comparison further suggested that genes involved in oxidative phosphorylation and mitochondrial sulfide oxidization pathway played important roles in the sulfide tolerance of the host mussels. Moreover, transcriptomic analysis of Na₂S stressed mussels confirmed the upregulation of oxidative phosphorylation and sulfide oxidization genes in response to sulfide exposure. Overall, this study provided a systematic transcriptional analysis of both the active bacterial community members and the host mussels, suggesting that the epibionts, endosymbionts, and mussel host collaborated on sulfide detoxification from extracellular to intracellular space to adapt to harsh H₂S-rich environments.

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 (H₂ 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 H₂ 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 H₂ 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 H₂ S adaptation by promoting functions needed for survival and reducing functions affected by H₂ S. This study lays the groundwork for further research to directly demonstrate the role of miRNAs in adaptation to H₂ 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.

Developmental Transcriptomic Analysis of the Cave-Dwelling Crustacean, Asellus aquaticus

Cave animals are a fascinating group of species often demonstrating characteristics including reduced eyes and pigmentation, metabolic efficiency, and enhanced sensory systems. Asellus aquaticus, an isopod crustacean, is an emerging model for cave biology. Cave and surface forms of this species differ in many characteristics, including eye size, pigmentation, and antennal length. Existing resources for this species include a linkage map, mapped regions responsible for eye and pigmentation traits, sequenced adult transcriptomes, and comparative embryological descriptions of the surface and cave forms. Our ultimate goal is to identify genes and mutations responsible for the differences between the cave and surface forms. To advance this goal, we decided to use a transcriptomic approach. Because many of these changes first appear during embryonic development, we sequenced embryonic transcriptomes of cave, surface, and hybrid individuals at the stage when eyes and pigment become evident in the surface form. We generated a cave, a surface, a hybrid, and an integrated transcriptome to identify differentially expressed genes in the cave and surface forms. Additionally, we identified genes with allele-specific expression in hybrid individuals. These embryonic transcriptomes are an important resource to assist in our ultimate goal of determining the genetic underpinnings of the divergence between the cave and surface forms.

Seedling response to water stress in valley oak (Quercus lobata) is shaped by different gene networks across populations

Drought is a major stress for plants, creating a strong selection pressure for traits that enable plant growth and survival in dry environments. Many drought responses are conserved species-wide responses, while others vary among populations distributed across heterogeneous environments. We tested how six populations of the widely distributed California valley oak (Quercus lobata) sampled from contrasting climates would differ in their response to soil drying relative to well-watered controls in a common environment by measuring ecophysiological traits in 93 individuals and gene expression (RNA-seq) in 42 individuals. Populations did not differ in their adjustment of turgor loss point during soil drying, suggesting a generalized species-wide response. Differential expression analysis identified 689 genes with a common response to treatment across populations and 470 genes with population-specific responses. Weighted gene co-expression network analysis (WGCNA) identified groups of genes with similar expression patterns that may be regulated together (gene modules). Several gene modules responded differently to water stress among populations, suggesting regional differences in gene network regulation. Populations from sites with a high mean annual temperature responded to the imposed water stress with significantly greater changes in gene module expression, indicating that these populations may be locally adapted to respond to drought. We propose that this variation among valley oak populations provides a mechanism for differential tolerance to the increasingly frequent and severe droughts in California.

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.

Transcriptomic Responses of Two Ecologically Divergent Populations of Japanese Mantis Shrimp (Oratosquilla oratoria) under Thermal Stress

Simple Summary

Rising ocean temperature would change the seawater chemistry and affect the external and internal physiology of crustaceans due to their lack of certain efficient temperature regulators. In addition, the infraspecific populations of crustaceans might also have different response strategies to the rising of temperature. Therefore, we identified the transcriptomic variations to the same thermal stress between ecologically divergent populations of Oratosquilla oratoria. The aim of this study was to investigate the population-specific function genes and relevant pathways in response to thermal stress in O. oratoria. The results showed that gene-expressed variation was in a population-specific pattern, which indicated that the local environment could lead to the evolvement of changes in gene regulation, ultimately leading to adaptive divergences. Additionally, we found several genes with large pleiotropic effects in the Zhoushan population, which might indicate that the regulation mechanisms of the Zhoushan population were more efficient than those of the Qingdao population under same thermal stress. The results provided some novel insights into the local adaptive differences of the infraspecific populations of O. oratoria and other crustaceans.

Abstract

Crustaceans are generally considered more sensitive to ocean warming due to their lack of certain efficient regulators. However, the alterations in the physiology and behavior of crustaceans in response to thermal stress differ vastly even among the infraspecific populations of heterogeneous landscapes. Consequently, understanding the impact of temperature fluctuation on crustacean infraspecific populations might be essential for maintaining a sustainable persistence of populations at existing locations. In the present study, we chose the Japanese mantis shrimp (Oratosquilla oratoria) as the representative crustacean population, and conducted transcriptome analyses in two divergent O. oratoria populations (the Zhoushan and Qingdao populations) under same thermal stress (20–28 °C) to identify the population-specific expression response to thermal stress. The results showed significant differences in gene expressions, GO terms and metabolic pathways between the two populations. We hypothesized that intraspecific mutations in the same or different genes might lead to thermal adaptive divergences. Temperature increases from 20–28 °C produced significant enrichment in GO terms and altered the metabolic pathways in the Zhoushan population despite the lack of differentially expressed unigenes. Therefore, several functional genes with large pleiotropic effects may underlie the response to thermal stress in the Zhoushan population. Furthermore, the most significantly enriched biological processes of the Qingdao population were associated with the state or activity of cells and its significant enriched pathways with genetic information processing as well as immune and environmental information processing. In contrast, the differentially regulated unigenes of the Zhoushan population were primarily involved in the regulatory cellular and transcription processes and the most significant pathways found were metabolic and digestive. Consequently, the regulatory mechanisms of the Zhoushan population are probably more efficient than those of the Qingdao population under the same thermal stress.

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

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 (H₂S) 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) H₂S 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 H₂S may differentially affect biodiversity. The presence of elevated H₂S 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 H₂S-rich stream microbiome in the region. Many H₂S-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.

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.

Intrasexual competition enhances reproductive isolation between locally adapted populations

During adaptation to different habitat types, both morphological and behavioral traits can undergo divergent selection. Males often fight for status in dominance hierarchies and rank positions predict reproductive success. Ecotypes with reduced fighting abilities should have low reproductive success when migrating into habitats that harbor ecotypes with superior fighting abilities. Livebearing fishes in the Poecilia mexicana-species complex inhabit not only regular freshwater environments, but also independently colonized sulfidic (H₂S-containing) habitats in three river drainages. In the current study, we found fighting intensities in staged contests to be considerably lower in some but not all sulfidic surface ecotypes and the sulfidic cave ecotype compared with populations from non-sulfidic surface sites. This is perhaps due to selection imposed by H₂S, which hampers oxygen uptake and transport, as well as cellular respiration. Furthermore, migrants from sulfidic habitats may lose fights even if they do not show overall reduced aggressiveness, as physiological performance is likely to be challenged in the non-sulfidic environment to which they are not adapted. To test this hypothesis, we simulated migration of H₂S-adapted males into H₂S-free waters, as well as H₂S-adapted cave-dwelling males into sulfidic surface waters. We found that intruders established dominance less often than resident males, independent of whether or not they showed reduced aggressiveness overall. Our study shows that divergent evolution of male aggressive behavior may also contribute to the maintenance of genetic differentiation in this system and we call for more careful evaluation of male fighting abilities in studies on ecological speciation.

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.