Genomics of Adaptation to Multiple Concurrent Stresses: Insights from Comparative Transcriptomics of a Cichlid Fish from One of Earth’s Most Extreme Environments, the Hypersaline Soda Lake Magadi in Kenya, East Africa

Characterization of transcriptome changes in saline stress adaptation on Leuciscus merzbacheri using PacBio Iso-Seq and RNA-Seq

Leuciscus merzbacheri is a native fish species found exclusively in the Junggar Basin in Xinjiang. It exhibits remarkable adaptability, thriving in varying water conditions such as the saline waters, the semi-saline water, and the freshwater. Despite its significant economic and ecological value, the underlying mechanisms of its remarkable salinity tolerance remain elusive. Our study marks the first time the full-length transcriptome of L. merzbacheri has been reported, utilizing RNA-Seq and PacBio Iso-Seq technologies. We found that the average length of the full-length transcriptome is 1,780 bp, with an N50 length of 2,358 bp. We collected RNA-Seq data from gill, liver, and kidney tissues of L. merzbacheri from both saline water and freshwater environments and conducted comparative analyses across these tissues. Further analysis revealed significant enrichment in several key functional gene categories and signalling pathways related to stress response and environmental adaptation. The findings provide a valuable genetic resource for further investigation into saline-responsive candidate genes, which will deepen our understanding of teleost adaptation to extreme environmental stress. This knowledge is crucial for the future breeding and conservation of native fish species.

Chromatin accessibility in gill tissue identifies candidate genes and loci associated with aquaculture relevant traits in tilapia

The Nile tilapia (Oreochromis niloticus) accounts for ∼9% of global freshwater finfish production however, extreme cold weather and decreasing freshwater resources has created the need to develop resilient strains. By determining the genetic bases of aquaculture relevant traits, we can genotype and breed desirable traits into farmed strains. We generated ATAC-seq and gene expression data from O. niloticus gill tissues, and through the integration of SNPs from 27 tilapia species, identified 1168 highly expressed genes (4% of all Nile tilapia genes) with highly accessible promoter regions with functional variation at transcription factor binding sites (TFBSs). Regulatory variation at these TFBSs is likely driving gene expression differences associated with tilapia gill adaptations, and differentially segregate in freshwater and euryhaline tilapia species. The generation of novel integrative data revealed candidate genes e.g., prolactin receptor 1 and claudin-h, genetic relationships, and loci associated with aquaculture relevant traits like salinity and osmotic stress acclimation.

The Biological Hierarchy, Time, and Temporal 'Omics in Evolutionary Biology: A Perspective

Sequencing data-genomics, transcriptomics, epigenomics, proteomics, and metabolomics-have revolutionized biological research, enabling a more detailed study of processes, ranging from subcellular to evolutionary, that drive biological organization. These processes, collectively, are responsible for generating patterns of phenotypic variation and can operate over dramatically different timescales (milliseconds to billions of years). While researchers often study phenotypic variation at specific levels of biological organization to isolate processes operating at that particular scale, the varying types of sequence data, or 'omics, can also provide complementary inferences to link molecular and phenotypic variation to produce an integrated view of evolutionary biology, ranging from molecular pathways to speciation. We briefly describe how 'omics has been used across biological levels and then demonstrate the utility of integrating different types of sequencing data across multiple biological levels within the same study to better understand biological phenomena. However, single-time-point studies cannot evaluate the temporal dynamics of these biological processes. Therefore, we put forward temporal 'omics as a framework that can better enable researchers to study the temporal dynamics of target processes. Temporal 'omics is not infallible, as the temporal sampling regime directly impacts inferential ability. Thus, we also discuss the role the temporal sampling regime plays in deriving inferences about the environmental conditions driving biological processes and provide examples that demonstrate the impact of the sampling regime on biological inference. Finally, we forecast the future of temporal 'omics by highlighting current methodological advancements that will enable temporal 'omics to be extended across species and timescales. We extend this discussion to using temporal multi-omics to integrate across the biological hierarchy to evaluate and link the temporal dynamics of processes that generate phenotypic variation.

Integrated application of multi-omics approach and biochemical assays provides insights into physiological responses to saline-alkaline stress in the gills of crucian carp (Carassius auratus)

Given the decline of freshwater resources in recent years, the accessible space for freshwater aquaculture is rapidly shrinking, and aquaculture in saline-alkaline water has become a critical approach to meet the rising demand. However, the molecular mechanism behind the adverse effects of saline-alkaline water on fish and the regulatory mechanism in fish tolerance remains unclear. Here, adult crucian carp (Carassius auratus) were exposed to 60 mmol/L NaHCO₃ for 30 days. It was observed that long-term carbonate alkalinity (CA) exposure not only caused gill oxidative stress but also changed the levels of several physiological parameters associated with ammonia transport, including blood ammonia, urea nitrogen (BUN), glutamine (Gln), and glutamine synthetase (GS). According to the metabolomics study, differential metabolites (DMs) engaged in various metabolic pathways, such as glycerophospholipid metabolism, sphingolipid metabolism, and arachidonic acid metabolism. In addition, transcriptomics data showed that differentially expressed genes (DEGs) were closely related to ammonia transport, apoptosis, and immunological response. In general, comprehensive multi-omics and biochemical analysis revealed that crucian carp might adopt Rh glycoprotein as a carrier to mediate ammonia transport and increase glutamine and urea synthesis under long-term high saline-alkaline stress to mitigate the adverse effects of blocked ammonia excretion. Simultaneously, saline-alkaline stress caused the destruction of the antioxidant system and the disorder of lipid metabolism in the crucian carp gills, which induced apoptosis and immunological response. To our knowledge, this is the first study to investigate fish's molecular and metabolic mechanisms under saline-alkaline stress using integrated metabolomics, transcriptomics, and biochemical assays. Overall, the results of this study provided new insights into the molecular mechanism behind the adverse effects of saline-alkaline water on fish and the regulatory mechanism in fish tolerance.

Transcriptome analysis reveals the mechanism of alkalinity exposure on spleen oxidative stress, inflammation and immune function of Luciobarbus capito

Saline-alkali land is distributed all over the world, and it affects the economic development of fisheries. The alkalinity in water is related to the accumulation of carbonate, so the is generally higher. To understand how alkalinity impacts the immune response in Luciobarbus capito, we performed transcriptomic profiles for spleen, the immune organ of Luciobarbus capito which were underwent alkalinity exposure. Totally there are 47,727,954, 53,987,820 and 51,398,546 high quality clean reads obtained from the control groups, and 46,996,982, 49,650,460 and 45,964,986 clean reads from the alkalinity exposure groups. Among them, 611 genes were differently expressed, including 534 upregulated and 77 down-regulated genes. The identified genes were enriched using databases of Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG). It was found that differentially expressed genes in Luciobarbus capito spleen tissue were enriched into 14 GO pathways, and differentially expressed genes in Luciobarbus capito spleen were enriched into 25 corresponding KEGG pathways under alkalinity stress. Inflammation and immune function genes and pathways were identified and validated by quantitative real-time RT-PCR. Our results showed that alkalinity exposure leads to inflammation and immunoregulation in spleen of Luciobarbus capito. These results provide new insights for unveiling the biological effects of alkalinity in Luciobarbus capito.

Salt to conserve: a review on the ecology and preservation of hypersaline ecosystems

When it comes to the investigation of key ecosystems in the world, we often omit salt from the ecological recipe. In fact, despite occupying almost half of the volume of inland waters and providing crucial services to humanity and nature, inland saline ecosystems are often overlooked in discussions regarding the preservation of global aquatic resources of our planet. As a result, our knowledge of the biological and geochemical dynamics shaping these environments remains incomplete and we are hesitant in framing effective protective strategies against the increasing natural and anthropogenic threats faced by such habitats. Hypersaline lakes, water bodies where the concentration of salt exceeds 35 g/l, occur mainly in arid and semiarid areas resulting from hydrological imbalances triggering the accumulation of salts over time. Often considered the 'exotic siblings' within the family of inland waters, these ecosystems host some of the most extremophile communities worldwide and provide essential habitats for waterbirds and many other organisms in already water-stressed regions. These systems are often highlighted as natural laboratories, ideal for addressing central ecological questions due to their relatively low complexity and simple food web structures. However, recent studies on the biogeochemical mechanisms framing hypersaline communities have challenged this archetype, arguing that newly discovered highly diverse communities are characterised by specific trophic interactions shaped by high levels of specialisation. The main goal of this review is to explore our current understanding of the ecological dynamics of hypersaline ecosystems by addressing four main research questions: (i) why are hypersaline lakes unique from a biological and geochemical perspective; (ii) which biota inhabit these ecosystems and how have they adapted to the high salt conditions; (iii) how do we protect biodiversity from increasing natural and anthropogenic threats; and (iv) which scientific tools will help us preserve hypersaline ecosystems in the future? First, we focus on the ecological characterisation of hypersaline ecosystems, illustrate hydrogeochemical dynamics regulating such environments, and outline key ecoregions supporting hypersaline systems across the globe. Second, we depict the diversity and functional aspects of key taxa found in hypersaline lakes, from microorganisms to plants, invertebrates, waterbirds and upper trophic levels. Next, we describe ecosystem services and discuss possible conservation guidelines. Finally, we outline how cutting-edge technologies can provide new insights into the study of hypersaline ecology. Overall, this review sheds further light onto these understudied ecosystems, largely unrecognised as important sources of unique biological and functional diversity. We provide perspectives for key future research avenues, and advocate that the conservation of hypersaline lakes should not be taken with 'a grain of salt'.

Gene expression remodelling and immune response during adaptive divergence in an African cichlid fish

Variation in gene expression contributes to ecological speciation by facilitating population persistence in novel environments. Likewise, immune responses can be of relevance in speciation driven by adaptation to different environments. Previous studies examining gene expression differences between recently diverged ecotypes have often relied on only one pair of populations, targeted the expression of only a subset of genes or used wild-caught individuals. Here, we investigated the contribution of habitat-specific parasites and symbionts and the underlying immunological abilities of ecotype hosts to adaptive divergence in lake-river population pairs of the cichlid fish Astatotilapia burtoni. To shed light on the role of phenotypic plasticity in adaptive divergence, we compared parasite and microbiota communities, immune response, and gene expression patterns of fish from natural habitats and a lake-like pond set-up. In all investigated population pairs, lake fish were more heavily parasitized than river fish, in terms of both parasite taxon composition and infection abundance. The innate immune response in the wild was higher in lake than in river populations and was elevated in a river population exposed to lake parasites in the pond set-up. Environmental differences between lake and river habitat and their distinct parasite communities have shaped differential gene expression, involving genes functioning in osmoregulation and immune response. Most changes in gene expression between lake and river samples in the wild and in the pond set-up were based on a plastic response. Finally, gene expression and bacterial communities of wild-caught individuals and individuals acclimatized to lake-like pond conditions showed shifts underlying adaptive phenotypic plasticity.

Exploring the Expression of Cardiac Regulators in a Vertebrate Extremophile: The Cichlid Fish Oreochromis (Alcolapia) alcalica

Although it is widely accepted that the cellular and molecular mechanisms of vertebrate cardiac development are evolutionarily conserved, this is on the basis of data from only a few model organisms suited to laboratory studies. Here, we investigate gene expression during cardiac development in the extremophile, non-model fish species, Oreochromis (Alcolapia) alcalica. We first characterise the early development of O. alcalica and observe extensive vascularisation across the yolk prior to hatching. We further investigate heart development by identifying and cloning O. alcalica orthologues of conserved cardiac transcription factors gata4, tbx5, and mef2c for analysis by in situ hybridisation. Expression of these three key cardiac developmental regulators also reveals other aspects of O. alcalica development, as these genes are expressed in developing blood, limb, eyes, and muscle, as well as the heart. Our data support the notion that O. alcalica is a direct-developing vertebrate that shares the highly conserved molecular regulation of the vertebrate body plan. However, the expression of gata4 in O. alcalica reveals interesting differences in the development of the circulatory system distinct from that of the well-studied zebrafish. Understanding the development of O. alcalica embryos is an important step towards providing a model for future research into the adaptation to extreme conditions; this is particularly relevant given that anthropogenic-driven climate change will likely result in more freshwater organisms being exposed to less favourable conditions.

Adaptation of the carbamoyl-phosphate synthetase enzyme in an extremophile fish

Tetrapods and fish have adapted distinct carbamoyl-phosphate synthase (CPS) enzymes to initiate the ornithine urea cycle during the detoxification of nitrogenous wastes. We report evidence that in the ureotelic subgenus of extremophile fish Oreochromis Alcolapia, CPS III has undergone convergent evolution and adapted its substrate affinity to ammonia, which is typical of terrestrial vertebrate CPS I. Unusually, unlike in other vertebrates, the expression of CPS III in Alcolapia is localized to the skeletal muscle and is activated in the myogenic lineage during early embryonic development with expression remaining in mature fish. We propose that adaptation in Alcolapia included both convergent evolution of CPS function to that of terrestrial vertebrates, as well as changes in development mechanisms redirecting CPS III gene expression to the skeletal muscle.

Transcriptome changes for Nile tilapia (Oreochromis niloticus) in response to alkalinity stress

Nile tilapia is an important economic fish in the world because of its fast growth, high meat yield and strong adaptability. It is more adaptable to high alkalinity than common freshwater fish and provides valuable material for developing alkaline-tolerant strains and understanding the adaptation mechanism of fish to extreme environmental stress. In this study, we employed high throughput RNA sequencing to reveal the tissues (gill, kidney and liver) transcriptome differences of O. niloticus at different carbonate alkalinities (FW, AW40 and AW60). A total of 1,369,381,790 raw reads were obtained, including 496,441,232 reads in FW group, 437,907,696 reads in AW40 and 435,032,862 reads in AW60. In addition, 484,555,626 reads in gill, 451,618,224 reads in kidney and 433,207,940 reads in liver. A large number of stress-regulated changes were detected comprehensively. We focused on 3 significantly change pathways (steroid biosynthesis, drug metabolism and protein digestion/absorption) and 17 DEGs (HMG-CoA reductase, UDP-glucuronosyltransferase, and carbonic anhydrase etc.) which were shared among compared groups (AW40 vs FW, AW60 vs FW, AW40 vs 60 AW60) in gill, kidney and liver, respectively. These pathways/genes are sensitive to alkalinity stress and crucial to the alkalinity adaptation of tilapia. Overall, we found a large number of candidate genes, which encode important regulators of stress tolerance and ultimately contribute to future alkaline-tolerant fish breeding. Among these genes, lipid metabolism (involving signal transduction), detoxification and immune related genes are more prominent to the response and adaptability of fish to alkalinity stress.

The Molecular Basis of Freshwater Adaptation in Prawns: Insights from Comparative Transcriptomics of Three Macrobrachium Species

Elucidating the molecular basis of adaptation to different environmental conditions is important because adaptive ability of a species can shape its distribution, influence speciation, and also drive a variety of evolutionary processes. For crustaceans, colonization of freshwater habitats has significantly impacted diversity, but the molecular basis of this process is poorly understood. In the current study, we examined three prawn species from the genus Macrobrachium (M. australiense, M. tolmerum, and M. novaehollandiae) to better understand the molecular basis of freshwater adaptation using a comparative transcriptomics approach. Each of these species naturally inhabit environments with different salinity levels; here, we exposed them to the same experimental salinity conditions (0‰ and 15‰), to compare expression patterns of candidate genes that previously have been shown to influence phenotypic traits associated with freshwater adaptation (e.g., genes associated with osmoregulation). Differential gene expression analysis revealed 876, 861, and 925 differentially expressed transcripts under the two salinities for M. australiense, M. tolmerum, and M. novaehollandiae, respectively. Of these, 16 were found to be unannotated novel transcripts and may be taxonomically restricted or orphan genes. Functional enrichment and molecular pathway mapping revealed 13 functionally enriched categories and 11 enriched molecular pathways that were common to the three Macrobrachium species. Pattern of selection analysis revealed 26 genes with signatures of positive selection among pairwise species comparisons. Overall, our results indicate that the same key genes and similar molecular pathways are likely to be involved with freshwater adaptation widely across this decapod group; with nonoverlapping sets of genes showing differential expression (mainly osmoregulatory genes) and signatures of positive selection (genes involved with different life history traits).

Bidirectional adaptive introgression between two ecologically divergent sparrow species

Natural hybrid zones can be used to dissect the mechanisms driving key evolutionary processes by allowing us to identify genomic regions important for establishing reproductive isolation and that allow for transfer of adaptive variation. We leverage whole-genome data in a system where two bird species, the saltmarsh (Ammospiza caudacuta) and Nelson's (A. nelsoni) sparrow, hybridize despite their relatively high background genetic differentiation and past ecological divergence. Adaptive introgression is plausible in this system because Nelson's sparrows are recent colonists of saltwater marshes, compared to the specialized saltmarsh sparrow that has a longer history of saltmarsh adaptation. Comparisons among whole-genome sequences of 34 individuals from allopatric and sympatric populations show that ongoing gene flow is shaping the genomic landscape, with allopatric populations exhibiting genome-wide F_ST estimates close to double of that observed in sympatry. We characterized patterns of introgression across the genome and identify regions that exhibit biased introgression into hybrids from one parental species. These regions offer compelling candidates for genes related to tidal marsh adaptations suggesting that adaptive introgression may be an important consequence of hybridization. These findings highlight the value of considering the landscapes of both genome-wide introgression and divergence when characterizing the evolutionary forces that drive speciation.

Phylogenomics uncovers early hybridization and adaptive loci shaping the radiation of Lake Tanganyika cichlid fishes

Lake Tanganyika is the oldest and phenotypically most diverse of the three East African cichlid fish adaptive radiations. It is also the cradle for the younger parallel haplochromine cichlid radiations in Lakes Malawi and Victoria. Despite its evolutionary significance, the relationships among the main Lake Tanganyika lineages remained unresolved, as did the general timescale of cichlid evolution. Here, we disentangle the deep phylogenetic structure of the Lake Tanganyika radiation using anchored phylogenomics and uncover hybridization at its base, as well as early in the haplochromine radiation. This suggests that hybridization might have facilitated these speciation bursts. Time-calibrated trees support that the radiation of Tanganyika cichlids coincided with lake formation and that Gondwanan vicariance concurred with the earliest splits in the cichlid family tree. Genes linked to key innovations show signals of introgression or positive selection following colonization of lake habitats and species' dietary adaptations are revealed as major drivers of colour vision evolution. These findings shed light onto the processes shaping the evolution of adaptive radiations.

Comparative transcriptomics of cyprinid minnows and carp in a common wild setting: a resource for ecological genomics in freshwater communities

Comparative transcriptomics can now be conducted on organisms in natural settings, which has greatly enhanced understanding of genome–environment interactions. Here, we demonstrate the utility and potential pitfalls of comparative transcriptomics of wild organisms, with an example from three cyprinid fish species (Teleostei:Cypriniformes). We present extensively filtered and annotated transcriptome assemblies that provide a valuable resource for studies of genome evolution (e.g. polyploidy), ecological and morphological diversification, speciation, and shared and unique responses to environmental variation in cyprinid fishes. Our results and analyses address the following points: (i) ‘essential developmental genes’ are shown to be ubiquitously expressed in a diverse suite of tissues across later ontogenetic stages (i.e. juveniles and adults), making these genes are useful for assessing the quality of transcriptome assemblies, (ii) the influence of microbiomes and other exogenous DNA, (iii) potentially novel, species-specific genes, and (iv) genomic rearrangements (e.g. whole genome duplication). The data we present provide a resource for future comparative work in cypriniform fishes and other taxa across a variety of sub-disciplines, including stress response, morphological diversification, community ecology, ecotoxicology, and climate change.

Transcriptomic response to thermal and salinity stress in introduced and native sympatric Palaemon caridean shrimps

Organisms develop local adaptations to cope with spatially and temporally variable environments such as estuarine habitats, where abiotic parameters such as salinity and temperature fluctuate continuously. Studying the regulation of gene expression in a variable environment allows us to understand the underlying molecular mechanisms of these adaptations and the relative roles of the genetic and plastic response. The transcriptomes of the European native Palaemon longirostris (PL) and the introduced P. macrodactylus (PM) shrimps are described and compared after an experiment simulating summer conditions in the Guadalquivir Estuary, Spain. Specimens, collected in the Guadalquivir Estuary, were maintained at a temperature and salinity of 20 °C and 5 ppt for the control, and 30 °C and 15 ppt for the stress treatment. A large amount of differential gene expression was observed: 16,013 and 2,594 for PL and PM respectively. Functionally annotated unigenes revealed some differences, with PL seemingly having to face stronger physiological stress than PM. Thus, PM seems to have greater resistance than PL under conditions of high temperature and salinity. These results constitute a step forward in the understanding of the underlying molecular mechanisms of genetic adaptation of native invertebrates, and alien taxa that have successfully invaded estuaries in temperate regions around the world.

The Role of Alternative Splicing and Differential Gene Expression in Cichlid Adaptive Radiation

Species diverge eco-morphologically through the continuous action of natural selection on functionally important structures, producing alternative adaptive morphologies. In cichlid fishes, the oral and pharyngeal jaws are such key structures. Adaptive variation in jaw morphology contributes to trophic specialization, which is hypothesized to fuel their rapid speciation in the East African Great Lakes. Much is known about the genes involved in cichlid jaw and craniofacial development. However, it is still unclear what salient sources of variation gave rise to trophic-niche specialization, facilitating adaptive radiation. Here, we explore two sources of transcriptional variation that may underlie species-specific disparities in jaw morphology. Using whole transcriptome RNA-sequencing, we analyze differences in gene expression and alternative splicing, at the end of postlarval development, in fully functional jaws of six species of cichlids from the Lake Tanganyika tribe Tropheini. Our data reveal a surprisingly high degree of alternative splicing events compared with gene expression differences among species and trophic types. This suggests that differential trophic adaptation of the jaw apparatus may have been shaped by transcriptional rewiring of splicing as well as gene expression variation during the rapid radiation of the Tropheini. Specifically, genes undergoing splicing across most species were found to be enriched for pharyngeal jaw gene ontology terms. Overall, jaw transcriptional patterns at postlarval developmental stage were highly dynamic and species-specific. In conclusion, this work indicates that shifts in alternative splicing could have played a more important role in cichlid adaptive radiation, and possibly adaptive radiation in general, than currently recognized.

Transcriptomic responses to environmental change in fishes: Insights from RNA sequencing

The need to better understand how plasticity and evolution affect organismal responses to environmental variability is paramount in the face of global climate change. The potential for using RNA sequencing (RNA-seq) to study complex responses by non-model organisms to the environment is evident in a rapidly growing body of literature. This is particularly true of fishes for which research has been motivated by their ecological importance, socioeconomic value, and increased use as model species for medical and genetic research. Here, we review studies that have used RNA-seq to study transcriptomic responses to continuous abiotic variables to which fishes have likely evolved a response and that are predicted to be affected by climate change (e.g., salinity, temperature, dissolved oxygen concentration, and pH). Field and laboratory experiments demonstrate the potential for individuals to respond plastically to short- and long-term environmental stress and reveal molecular mechanisms underlying developmental and transgenerational plasticity, as well as adaptation to different environmental regimes. We discuss experimental, analytical, and conceptual issues that have arisen from this work and suggest avenues for future study.

Signals of selection in conditionally expressed genes in the diversification of three horned beetle species

Species radiations may be facilitated by phenotypic differences already present within populations, such as those arising through sex-specific development or developmental processes biased towards particular reproductive or trophic morphs. We sought to test this hypothesis by utilizing a comparative transcriptomic approach to contrast among- and within-species differentiation using three horned beetle species in the genus Onthophagus. These three species exhibit differences along three phenotypic axes reflective of much of the interspecific diversity present within the genus: horn location, polarity of sexual dimorphism and degree of nutritional sensitivity. Our approach combined de novo transcript assembly, assessment of amino acid substitutions (dN/dS) across orthologous gene pairs and integration of gene function and conditional gene expression data. We identified 17 genes across the three species pairs related to axis patterning, development and metabolism with dN/dS > 1 and detected elevated dN/dS in genes related to metabolism and biosynthesis in the most closely related species pair, which is characterized by a loss of nutritional polyphenism and a reversal of sexual dimorphism. Further, we found that genes that are conditionally expressed (i.e. as a function of sex, nutrition or body region) within one of our focal species also showed significantly stronger signals of positive or relaxed purifying selection between species divergent along the same morphological axis (i.e. polarity of sexual dimorphism, degree of nutritional sensitivity or location of horns). Our findings thus reveal a positive relationship between intraspecific differentiation due to condition-specific development and genetic divergences among species.

Comprehensive transcriptomic analysis of Tibetan Schizothoracinae fish Gymnocypris przewalskii reveals how it adapts to a high altitude aquatic life

Background

Understanding the genetic basis of adaptation to high altitude life is of paramount importance for preserving and managing genetic diversity in highland animals. This objective has been addressed mainly in terrestrial fauna but rarely in aquatic animals. Tibetan Schizothoracinae fish is the ideal model system in evolutionary biology, carrying key insights into evolutionary genetics of speciation and adaptation at high altitude. Gymnocypris przewalskii is the newly formed Schizothoracinae fish species in the Tibetan Plateau, inhabits chronic cold, extreme saline and alkaline aquatic environment in Lake Qinghai, thus evolving the unique genomic signatures to adapt extremely severe environments.

Results

To characterize its genomic features, we assembled de novo transcriptome of G. przewalskii from Lake Qinghai. Intriguingly, by comparative genomic analyses of G. przewalskii and 8 other fish species, we identified potential expansions in gene families related to energy metabolism, transport and developmental functions, possibly underlying the adaptation to these environmental stresses. Through comprehensive molecular evolution analyses, we found that sets of genes controlling mitochondrion, ion homoeostasis, acid-base balance and innate immunity show significant signals of positive selection. Compared to previous studies on highland fishes, we failed to identify any positively selected genes related to hypoxia response.

Conclusions

Our findings provide comprehensive insights into the genetic basis of teleost fish that underlie their adaptation to extreme high altitude aquatic life on the Tibetan Plateau.

Electronic supplementary material

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

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.