Analysis of hypoxia-inducible factor alpha polyploidization reveals adaptation to Tibetan Plateau in the evolution of schizothoracine fish

Evolution of Key Oxygen-Sensing Genes Is Associated with Hypoxia Tolerance in Fishes

Low dissolved oxygen (hypoxia) is recognized as a major threat to aquatic ecosystems worldwide. Because oxygen is paramount for the energy metabolism of animals, understanding the functional and genetic drivers of whole-animal hypoxia tolerance is critical to predicting the impacts of aquatic hypoxia. In this study, we investigate the molecular evolution of key genes involved in the detection of and response to hypoxia in ray-finned fishes: the prolyl hydroxylase domain (PHD)-hypoxia-inducible factor (HIF) oxygen-sensing system, also known as the EGLN (egg-laying nine)-HIF oxygen-sensing system. We searched fish genomes for HIFA and EGLN genes, discovered new paralogs from both gene families, and analyzed protein-coding sites under positive selection. The physicochemical properties of these positively selected amino acid sites were summarized using linear discriminants for each gene. We employed phylogenetic generalized least squares to assess the relationship between these linear discriminants for each HIFA and EGLN and hypoxia tolerance as reflected by the critical oxygen tension (Pcrit) of the corresponding species. Our results demonstrate that Pcrit in ray-finned fishes correlates with the physicochemical variation of positively selected sites in specific HIFA and EGLN genes. For HIF2A, two linear discriminants captured more than 90% of the physicochemical variation of these sites and explained between 20% and 39% of the variation in Pcrit. Thus, variation in HIF2A among fishes may contribute to their capacity to cope with aquatic hypoxia, similar to its proposed role in conferring tolerance to high-altitude hypoxia in certain lineages of terrestrial vertebrates.

Linking HIF oxygen-sensing system diversity to hypoxia fitness in Eleutheronema: Molecular characterization and transcriptional response to hypoxia exposure

Hypoxia is a mounting environmental problem affecting coastal waters globally, posing severe consequences for biodiversity and marine life. Metazoans respond to hypoxia stress via the hypoxia-inducible factor (HIF) pathway, but few studies have addressed the gene diversity of the functionally important HIF-pathway. Understanding whether functional diversity exists in the HIF-pathway is a key first step in identifying genes that may impact hypoxia fitness. Here, we leveraged whole-genome resequencing data and bioinformatics tools to identify the key members of the HIF-pathway (HIFα/β, EGLN, and VHL) and genetic diversity in the threatened Eleutheronema. Phylogenetic analysis revealed that teleost-specific duplicates of epas1 (epas1a/b) were followed by the loss of one of each hif1α and hif1αl in Eleutheronema species. Strong collinearity and similarity of gene characteristics suggested the functional conservation of the HIF-pathway during Eleutheronema evolution. Purifying selection was the major theme in HIF-pathway evolution, leading to a reduction in genetic diversity. Substantially low nucleotide diversity of the HIF-pathway was observed among populations, which might indicate the loss of hypoxia fitness in Eleutheronema. Additionally, the normoxic presence of the HIF-pathway differed among tissues and was species-dependent, indicating their diverse roles during development. Significant regulation of HIF-pathway expression levels was observed across tissues under hypoxic conditions, suggesting critical roles in the hypoxia stress response. Moreover, variant molecular characters suggested different roles in response to hypoxia of the HIF-pathway, which were reflected in the different expression patterns across tissues. Our present study provides novel insights into the interplay between gene diversity within the HIF-pathway and hypoxia fitness in threatened Eleutheronema. We highlighted the importance of HIF-pathway-mediated transcriptional regulation in response to hypoxia stress, which provided valuable information for the genetic mechanisms underlying hypoxia adaptation in fish. The bioinformatic methods developed here have broad applications for other species.

Post-Subfunctionalization Functions of HIF-1αA and HIF-1αB in Cyprinid Fish: Fine-Tuning Mitophagy and Apoptosis Regulation Under Hypoxic Stress

Hypoxia-inducible factor 1 (HIF-1) is a crucial transcriptional factor that can restore oxygen balance in the body by regulating multiple vital activities. Two HIF-1α copies were retained in cyprinid fish after experiencing a teleost-specific genome duplication. How the "divergent collaboration" of HIF-1αA and HIF-1αB proceeds in regulating mitophagy and apoptosis under hypoxic stress in cells of cyprinid fish remains unclear. In this study, zebrafish HIF-1αA/B expression plasmids were constructed and transfected into the epithelioma papulosum cyprini cells and were subjected to hypoxic stress. HIF-1αA induced apoptosis through promoting ROS generation and mitochondrial depolarization when cells were subjected to oxygen deficiency. Conversely, HIF-1αB was primarily responsible for mitophagy induction, prompting ATP production to mitigate apoptosis. HIF-1αA did not induce mitophagy in the mitochondria and lysosomes co-localization assay but it was involved in the regulation of different mitophagy pathways. Over-expression of HIF-1αA increased the expression of bnip3, fundc1, Beclin1, and foxo3, suggesting it has a dual role in mitochondrial autophagy and cell death. Each duplicated copy also experienced functional divergence and target shifting in the regulation of complexes in the mitochondrial electron transport chain (ETC). Our findings shed light on the post-subfunctionalization function of HIF-1αA and HIF-1αB in zebrafish to fine-tune regulation of mitophagy and apoptosis following hypoxia exposure.

Role of HIF in fish inflammation

The hypoxia-inducing factor (HIF) is a central transcription factor in cellular oxygen sensing and regulation. It is common that the inflammation always appears in many diseases, like infectious diseases in fishes, and the inflammation is often accompanied by hypoxia, as a hallmark of inflammation. Besides coordinating cellular responses to low oxygen, HIF-mediated hypoxia signaling pathway is also crucial for immune responses such as the regulations of innate immune cell phenotype and function, as well as metabolic reprogramming under the inflammation. However, the understanding of the molecular mechanisms by which HIFs regulate the inflammatory response in fish is still very limited. Here, we review the characteristics of HIF as well as its roles in innate immune cells and the infections caused by bacteria and viruses. The regulatory effects of HIF on the metabolic reprogramming of innate immune cells are also discussed and the future research directions are outlooked. This paper will serve as a reference for elucidating the molecular mechanism of HIF regulating inflammation and identifying treatment strategies to target HIF for fish disease.

Functional exploration of SNP mutations in HIF2αb gene correlated with hypoxia tolerance in blunt snout bream (Megalobrama amblycephala)

Blunt snout bream (Megalobrama amblycephala) is sensitive to hypoxia environment. Hypoxia-inducible factor (HIF) is the most critical factor in the HIF pathway, which strictly regulates the hypoxia stress process of fish. In this study, we found six hifα genes in blunt snout bream that demonstrated different expressions under hypoxia conditions. In HEK293T cells, all six hifαs were detected to activate the HRE region by luciferase reporter assay. More importantly, we identified two linkage-disequilibrium SNP sites at exon 203 and 752 of the hif2αb gene in blunt snout bream. Haplotype II (A²⁰³A⁷⁵²) and its homozygous diplotype II (A²⁰³A²⁰³A⁷⁵²A⁷⁵²) appeared frequently in a selected strain of blunt snout bream with hypoxia tolerance. Diplotype II has a lower oxygen tension threshold for loss of equilibrium (LOE_crit) over a similar range of temperatures. Moreover, its erythrocyte number increased significantly (p < 0.05) than those in diplotype I and diplotype III strains at 48 h of hypoxia. The enzymes related with hypoxia tolerant traits, i.e., reduced glutathione, superoxide dismutase, and catalase, were also significantly (p < 0.05) induced in diplotype II than in diplotype I or III. In addition, the expression of epo in the liver of diplotype II was significantly (p < 0.01) higher than that in the diplotype I or III strains at 48 h of hypoxia. Taken together, our results found that the hypoxia-tolerant-related diplotype II of hif2αb has the potential to be used as a molecular marker in future genetic breeding of hypoxia-tolerant strain.

CRISPR-Cas9-Mediated NPC1 Gene Deletion Enhances HEK 293 T Cell Adhesion by Regulating E-Cadherin

NPC1 gene encodes a transmembrane glycoprotein on the late endosome/lysosomal membrane. Its mutation leads to a rare and aggravated autosomal recessive neurovisceral condition, termed Niemann-Pick disease type C1 (NPC1), which is characterized by progressive neurodegeneration, visceral symptoms, and premature death. To investigate the influence of NPC1 gene deletion on cell morphology, adhesion, proliferation, and apoptosis, CRISPR-Cas9 technology was used to knockout the NPC1 gene in HEK 293 T cells. Sanger sequencing, western blotting, and immunofluorescence were used to confirm successful NPC1 ablation. Filipin staining results indicated that deletion of NPC1 gene led to accumulation of unesterified cholesterol in HEK 293 T cells. Phalloidin staining results revealed cell aggregation, synapse shortening, nuclear enlargement, and cytoskeleton filamentous actin thinning in HEK 293 T cells with NPC1 gene mutation. Furthermore, NPC1 gene mutated HEK 293 T cell showed enhanced cell adhesion, inhibited cell proliferation, and increased cell apoptosis. In addition, NPC1 gene mutations significantly increased the protein expression levels of E-cadherin and γ-catenin and significantly decreased the protein expression levels of Wnt 3a, c-Myc, and cyclin D1. These results suggest that NPC1 may regulate cell adhesion by affecting the cadherin-catenin complex through E-cadherin, and that the classical Wnt signaling pathway may be inhibited by restricting β-catenin from entering the nucleus to inhibit cell proliferation.

Genomic analysis of hypoxia inducible factor alpha in ray-finned fishes reveals missing Ohnologs and evidence of widespread positive selection

As aquatic hypoxia worsens on a global scale, fishes will become increasingly challenged by low oxygen, and understanding the molecular basis of their response to hypoxia may help to better define the capacity of fishes to cope with this challenge. The hypoxia inducible factor (HIF) plays a critical role in the molecular response to hypoxia by activating the transcription of genes that serve to improve oxygen delivery to the tissues or enhance the capacity of tissues to function at low oxygen. The current study examines the molecular evolution of genes encoding the oxygen-dependent HIFα subunit (HIFA) in the ray-finned fishes (Actinopterygii). Genomic analyses demonstrate that several lineages retain four paralogs of HIFA predicted from two rounds of genome duplication at the base of vertebrate evolution, broaden the known distribution of teleost-specific HIFA paralogs, and provide evidence for salmonid-specific HIFA duplicates. Evolution of the HIFA gene family is characterized by widespread episodic positive selection at amino acid sites that potentially mediate protein stability, protein-protein interactions, and transcriptional regulation. HIFA transcript abundance depends upon paralog, tissue, and fish lineage. A phylogenetically-informed gene nomenclature is proposed along with avenues for future research on this critical family of transcription factors.

Insights Into miRNA-mRNA Regulatory Mechanisms of Cold Adaptation in Gymnocypris eckloni: Ubiquitin-Mediated Proteolysis Is Pivotal for Adaptive Energy Metabolism

This study aimed to understand cold stress adaptations mechanism in fish. Thus, the transcriptional response to cold conditions in Gymnocypris eckloni was evaluated using RNA-seq and microRNA (miRNA)-seq analyses. Low-temperature (LT) group G. eckloni was cultivated outdoors in waters cooled to 2–4°C for 3 weeks, while individuals in the control temperature (CT) group were exposed to 14–16°C. Significantly different responses were observed in both mRNA and miRNA expression profiles, with more mRNAs (1,833 and 1,869 mRNAs were up- and downregulated, respectively) and fewer miRNAs (15 and 6 were up- and downregulated, respectively) observed in the LT group individuals relative to the CT group individuals. A miRNA-mRNA network involved in the regulation of G. eckloni responses to cold stress was constructed; this network included ubiquitin-mediated proteolysis, protein processing, and oxidative phosphorylation. These results provided new insights into mechanisms of cold tolerance by fish, including decreased metabolic activity in addition to proteolysis.

Parallel introgression, not recurrent emergence, explains apparent elevational ecotypes of polyploid Himalayan snowtrout

The recurrence of similar evolutionary patterns within different habitats often reflects parallel selective pressures acting upon either standing or independently occurring genetic variation to produce a convergence of phenotypes. This interpretation (i.e. parallel divergences within adjacent streams) has been hypothesized for drainage-specific morphological 'ecotypes' observed in polyploid snowtrout (Cyprinidae: Schizothorax). However, parallel patterns of differential introgression during secondary contact are a viable alternative hypothesis. Here, we used ddRADseq (N = 35 319 de novo and N = 10 884 transcriptome-aligned SNPs), as derived from Nepali/Bhutanese samples (N = 48 each), to test these competing hypotheses. We first employed genome-wide allelic depths to derive appropriate ploidy models, then a Bayesian approach to yield genotypes statistically consistent under the inferred expectations. Elevational 'ecotypes' were consistent in geometric morphometric space, but with phylogenetic relationships at the drainage level, sustaining a hypothesis of independent emergence. However, partitioned analyses of phylogeny and admixture identified subsets of loci under selection that retained genealogical concordance with morphology, suggesting instead that apparent patterns of morphological/phylogenetic discordance are driven by widespread genomic homogenization. Here, admixture occurring in secondary contact effectively 'masks' previous isolation. Our results underscore two salient factors: (i) morphological adaptations are retained despite hybridization and (ii) the degree of admixture varies across tributaries, presumably concomitant with underlying environmental or anthropogenic factors.

The evolutionary and physiological significance of the Hif pathway in teleost fishes

The hypoxia-inducible factor (HIF) pathway is a key regulator of cellular O2 homeostasis and an important orchestrator of the physiological responses to hypoxia (low O2) in vertebrates. Fish can be exposed to significant and frequent changes in environmental O2, and increases in Hif-α (the hypoxia-sensitive subunit of the transcription factor Hif) have been documented in a number of species as a result of a decrease in O2. Here, we discuss the impact of the Hif pathway on the hypoxic response and the contribution to hypoxia tolerance, particularly in fishes of the cyprinid lineage, which includes the zebrafish (Danio rerio). The cyprinids are of specific interest because, unlike in most other fishes, duplicated paralogs of the Hif-α isoforms arising from a teleost-specific genome duplication event have been retained. Positive selection has acted on the duplicated paralogs of the Hif-α isoforms in some cyprinid sub-families, pointing to adaptive evolutionary change in the paralogs. Thus, cyprinids are valuable models for exploring the evolutionary significance and physiological impact of the Hif pathway on the hypoxic response. Knockout in zebrafish of either paralog of Hif-1α greatly reduces hypoxia tolerance, indicating the importance of both paralogs to the hypoxic response. Here, with an emphasis on the cardiorespiratory system, we focus on the role of Hif-1α in the hypoxic ventilatory response and the regulation of cardiac function. We explore the effects of the duration of the hypoxic exposure (acute, sustained or intermittent) on the impact of Hif-1α on cardiorespiratory function and compare relevant data with those from mammalian systems.

Genomic Signature of Shifts in Selection and Alkaline Adaptation in Highland Fish

Understanding how organisms adapt to aquatic life at high altitude is fundamental in evolutionary biology. This objective has been addressed primarily related to hypoxia adaptation by recent comparative studies, whereas highland fish has also long suffered extreme alkaline environment, insight into the genomic basis of alkaline adaptation has rarely been provided. Here, we compared the genomes or transcriptomes of 15 fish species, including two alkaline tolerant highland fish species and their six alkaline intolerant relatives, three alkaline tolerant lowland fish species, and four alkaline intolerant species. We found putatively consistent patterns of molecular evolution in alkaline tolerant species in a large number of shared orthologs within highland and lowland fish taxa. Remarkably, we identified consistent signatures of accelerated evolution and positive selection in a set of shared genes associated with ion transport, apoptosis, immune response, and energy metabolisms in alkaline tolerant species within both highland and lowland fish taxa. This is one of the first comparative studies that began to elucidate the consistent genomic signature of alkaline adaptation shared by highland and lowland fish. This finding also highlights the adaptive molecular evolution changes that support fish adapting to extreme environments at high altitude.

Specific cyprinid HIF isoforms contribute to cellular mitochondrial regulation

Hypoxia-inducible factor 1 (HIF-1) functions as a master regulator of the cellular response to hypoxic stress. Two HIF-1α paralogs, HIF-1αA and HIF-1αB, were generated in euteleosts by the specific, third round of genome duplication, but one paralog was later lost in most families with the exception of cyprinid fish. How these duplicates function in mitochondrial regulation and whether their preservation contributes to the hypoxia tolerance demonstrated by cyprinid fish in freshwater environments is not clear. Here we demonstrated the divergent function of these two zebrafish Hif-1a paralogs through cellular approaches. The results showed that Hif-1aa played a role in tricarboxylic acid cycle by increasing the expression of Citrate synthase and the activity of mitochondrial complex II, and it also enhanced mitochondrial membrane potential and ROS production by reducing free Ca²⁺ in the cytosol. Hif-1ab promoted intracellular ATP content by up-regulating the activity of mitochondrial complexes I, III and IV and the expression of related genes. Furthermore, both the two zebrafish Hif-1a paralogs promoted mitochondrial mass and the expression level of mtDNA, contributing to mitochondrial biogenesis. Our study reveals the divergent functions of Hif-1aa and Hif-1ab in cellular mitochondrial regulation.

Experimental study of hypoxia-induced changes in gene expression in an Asian pika, Ochotona dauurica

Acclimation to environmental changes driven by alterations in gene expression will serve as an important response for some species facing rapid Anthropogenic climate change. Pikas, genus Ochotona, are particularly vulnerable to climate change and current trends suggest that only the highest, coldest elevations within their ranges may remain suitable habitat for these species. In this study we aimed to assess the role of changes in gene expression in potentially facilitating elevational movements in pikas by measuring gene expression in the only known captive pika population, Ochotona dauurica, in response to hypoxic conditions. Using a controlled experiment, we exposed four male pikas to oxygen concentrations characteristic of sea-level, 2,000 m, and 4,000 m for 5 days each. Using blood samples collected after each treatment, we used RNAseq to determine if candidate pathways were undergoing significant changes in gene expression at different levels of oxygen (~100%, ~77%, and ~61% of sea-level oxygen concentrations). Gene set enrichment analyses showed that gene sets associated with the oxidative phosphorylation pathway and electron transport chain were significantly enriched for up-regulated genes in the 4,000 m samples compared to samples from the same individuals at lower-elevation conditions. Up-regulation of these pathways is consistent with known mechanisms of oxygen compensation. Our results suggest that these pikas have the acclimation capacity to tolerate oxygen concentrations characteristic of any elevation within their species range and that gene expression can be changed in a matter of days to accommodate drastically different oxygen concentrations. Thus, rapid and radical elevational movements that may be required of some pika species to avoid warmer temperatures in the Anthropocene will likely not be limited by hypoxic stress.

Transcriptome-Wide Patterns of the Genetic and Expression Variations in Two Sympatric Schizothoracine Fishes in a Tibetan Plateau Glacier Lake

Sympatric speciation remains a central focus of evolutionary biology. Although some evidence shows speciation occurring in this way, little is known about the gene expression evolution and the characteristics of population genetics as species diverge. Two closely related Gymnocypris fish (Gymnocypris chui and Gymnocypris scleracanthus), which come from a small glacier lake in the Tibetan Plateau, Lake Langcuo, exist a possible incipient sympatric adaptive ecological speciation. We generated large amounts of RNA-Seq data from multiple individuals and tissues from each of the two species and compared gene expression patterns and genetic polymorphisms between them. Ordination analysis separated samples by organ rather than by species. The degree of expression difference between organs within and between species was different. Phylogenetic analyses indicated that the two closely related taxa formed a monophyletic complex. Population structure analysis displayed two distinctly divergent clusters of G. chui and G. scleracanthus populations. By contrast, G. scleracanthus population genetic diversity is higher than that of G. chui. Considerable sites of the two populations were differentiated with a coefficient of F_ST = 0.25–0.50, implying that a small proportion of loci nevertheless exhibited deep divergence in two comparisons. Concomitantly, putatively selected genes during speciation revealed functional categories are enriched in bone morphogenesis, cell growth, neurogenetics, enzyme activity, and binding activity in G. chui population. In contrast, nutrition and localization were highlighted in G. scleracanthus. Collectively, morphological traits and dietary preference combine with genetic variation and expression variation, probably contributed to the incipient speciation of two sympatric populations.

Analysis of Multiplicity of Hypoxia-Inducible Factors in the Evolution of Triplophysa Fish (Osteichthyes: Nemacheilinae) Reveals Hypoxic Environments Adaptation to Tibetan Plateau

HIF (Hypoxia-inducible factor) gene family members function as master regulators of cellular and systemic oxygen homeostasis during changes in oxygen availability. Qinghai-Tibet Plateau is a natural laboratory for for long-term hypoxia and cold adaptation. In this context, T. scleroptera that is restricted to >3500 m high-altitude freshwater rivers was selected as the model to compare with a representative species from the plain, P. dabryanus. We cloned different HIF-α and carried out a phylogenetic analysis from invertebrates to vertebrates for identifying HIF-α genes and analyzing their evolutionary history. Intriguingly, the HIF-α has undergone gene duplications might be due to whole-genome duplication (WGD) events during evolution. PAML analysis indicated that HIF-1αA was subjected to positive selection acted on specific sites in Triplophysa lineages. To investigate the relationship between hypoxia adaptation and the regulation of HIF-α stability by pVHL in plateau and plain fish, a series of experiments were carried out. Comparison the luciferase transcriptional activity and protein levels of HIF-αs and the differing interactions of HIF-αs with pVHL, show clear differences between plateau and plain fish. T. scleroptera pVHL could enhance HIF-α transcriptional activity under hypoxia, and functional validation through pVHL protein mutagenesis showed that these mutations increased the stability of HIF-α and its hetero dimerization affinity to ARNT. Our research shows that missense mutations of pVHL induced evolutionary molecular adaptation in Triplophysa fishes living in high altitude hypoxic environments.

Phylogenetic perspective on the relationships and evolutionary history of the Acipenseriformes

The Acipenseriformes, as one of the earliest extant vertebrates, plays an important role in the evolution of fishes and even the whole vertebrates. Here we collected and analyzed all complete mitochondrial genomes of Acipenseriformes species. Phylogenetic analyses demonstrated that the polytomous branch included Acipenseridae and Polyodontidae formed five clades. The Polyodontidae clade and the Scaphirhynchus clade both were monophyletic group, whereas the Acipenser species and the Huso species both were polyphyletic group. The Bayesian divergence times showed that the origin time for Acipenseriformes was at 318.0 Mya, which was similar to the some previous results of 312.1 Mya, 346.9 Mya and 389.7 Mya. The result was in good consistent with the paleontological data available and the split time of the Pacific and Atlantic Oceans from the Jurassic to the Cretaceous (Laurasia splits in North America and Eurasia). The dN/dS ratios showed the evolutionary rates gradually slow down in five major Acipenseriformes clades from the Clade A (the Pacific sturgeons species) to Clade C (the genus Scaphirhynchus), which was related to the process of geographical formation.

Hypoxia-inducible factor 1α from a high-altitude fish enhances cytoprotection and elevates nitric oxide production in hypoxic environment

Hypoxia-inducible factors (HIFs) are master transcription factor regulating hypoxic responses in vertebrates. Species of Schizothoracine, a sub-family of cyprinidae, are highly endemic to the hypoxic Qinghai-Tibetan Plateau (QTP). What roles the HIFs play in hypoxic adaptation in the Schizothoracine fish is little known. In this study, the HIF-1α/B gene from Gymnocypris dobula (Gd) was characterized. The predicted protein for Gd-HIF-1α/B contains the main domains (bHLH, PAS, PAC, ODD, N-TAD, and C-TAD). Moreover, a specific mutation that the proline hydroxylation motif (LXXLAP) mutated into PxxLAP was observed in Gd-HIF-1α/B CODD domain, which may lead to changes in the function. To clarify whether HIF-1α/B of G. dobula possesses hypoxic adaptive features, Gd-HIF1α/B and Schizothorax prenanti-HIF1α/B (Sp-HIF1α/B) were cloned into an expression vector and transfected into 293T cells. Cell viability was found to be significantly higher in cells transfected with Gd-HIF-1α/B than those transfected with Sp-HIF-1α/B under hypoxic conditions. In addition, G. dobula HIF-1α/B showed stronger activity in transactivating the expression of nitric oxide (NO)-synthesizing enzyme, NOS2B under hypoxia stresses than the orthologous gene from S. prenanti, which were accompanied with upregulated expressions of NOS2B in heart of G. dobula, which may attribute to elevated NO levels detected in G. dobula than the lower land species. These results indicated that the HIF-1α plays an important role in mediating the iNOS signaling system in the process of evolutionary adaptation of the Schizothoracine to the highland environment.

Transcriptomic signature of rapidly evolving immune genes in a highland fish

Recent genome-wide studies have begun to elucidate the genomic basis of hypoxia, long-term cold and high saline and alkaline adaptation in highland fish, and a number of key genes contributed to its highland adaptation were identified. An increasing number of studies indicated that immune genes of Tibetan endemic fish species underwent positive selection towards functional shift, while the insight into immune gene repertoire of Tibetan highland fishes from genome-wide studies has largely lagged behind. In this study, we performed one of the first comparative genomics study in particular focusing on the signatures of immune genes in a highland fish, Gymnocypris przewalskii based on immune-relevant tissue transcriptome assemblies. We identified seven putative rapidly evolving immune genes with elevated molecular evolutionary rate (dN/dS) relative to lowland fish species. Using tissue-transcriptome data, we found most of rapidly evolving immune genes were broadly expressed in head-kidney, spleen, gills and skin tissues, which significantly enriched for complement activation and inflammatory response processes. In addition, we found a set of complement activation related genes underwent accelerated evolution and showed consistently repressed expression patterns in response to parasite Ichthyophthirius multifiliis infection. Moreover, we detected a number of immune genes involved in adaptive immune system exhibited distinct signature of upregulated expression patterns after parasite infection. Taken together, this study provided putative transcriptomic signatures of rapidly evolving immune genes, and will gain the insight into Schizothoracine fish adaptation to high-altitude extreme aquatic environments including diversified pathogen challenge.

Disentangling the interplay of positive and negative selection forces that shaped mitochondrial genomes of Gammarus pisinnus and Gammarus lacustris

We hypothesized that the mitogenome of Gammarus lacustris (GL), native to the Qinghai-Tibet Plateau, might exhibit genetic adaptations to the extreme environmental conditions associated with high altitudes (greater than 3000 m). To test this, we also sequenced the mitogenome of Gammarus pisinnus (GP), whose native range is close to the Tibetan plateau, but at a much lower altitude (200-1500 m). The two mitogenomes exhibited conserved mitochondrial architecture, but low identity between genes (55% atp8 to 76.1% cox1). Standard (homogeneous) phylogenetic models resolved Gammaridae as paraphyletic, but 'heterogeneous' CAT-GTR model as monophyletic. In indirect support of our working hypothesis, GL, GP and Gammarus fossarum exhibit evidence of episodic diversifying selection within the studied Gammaroidea dataset. The mitogenome of GL generally evolves under a strong purifying selection, whereas GP evolves under directional (especially pronounced in atp8) and/or relaxed selection. This is surprising, as GP does not inhabit a unique ecological niche compared to other gammarids. We propose that this rapid evolution of the GP mitogenome may be a reflection of its relatively recent speciation and heightened non-adaptive (putatively metabolic rate-driven) mutational pressures. To test these hypotheses, we urge sequencing mitogenomes of remaining Gammarus species populating the same geographical range as GP.

Effects of dietary supplementation of three strains of Lactococcus lactis on HIFs genes family expression of the common carp following Aeromonas hydrophila infection

HIFs (Hypoxia inducible factors) are the main regulators of the expression change of oxygen-dependent genes, in addition, they also play important roles in immune regulation. HIFs participate in infectious diseases and inflammatory responses, providing us a new therapeutic target for the treatment of diseases. In this study, 16 HIFs were identified in common carp genome database. Comparative genomics analysis showed large expansion of HIF gene family and approved the four round whole genome duplication (WGD) event in common carp. To further understand the function of HIFs, the domain architectures were predicted. All HIF proteins had the conserved HLH-PAS domain, which were essential for them to form dimer and bind to the downstream targets. The differences in domain of HIFα and HIFβ might result in their different functions. Phylogenetic analysis revealed that all HIFs were divided into two subfamilies and the HIFs in common carp were clustered with their teleost counterparts indicating they are highly conservative during evolution. In addition, the tissue distribution was examined by RT-PCR showed that most of HIF genes had a wide range of tissue distribution but exhibited tissue-specific expression patterns. The expression divergences were observed between the copy genes, for example, HIF1A-1, HIF2A-1, ARNT-2 had wide tissue distribution while their copies had limited tissue distribution, proving the function divergence of copies post the WGD event. In order to find an effective activation of HIFs and apply to treatment of aquatic diseases, we investigate the dietary supplementation effects of different strains of Lactococcus lactis on the expression of HIFα subfamily members in kidney of common carp infected with A. hydrophila. In addition, all of the HIF genes have a high expression in the early stages of infection, and decreased in the treatment time point of 48 h in common carp. This phenomenon confirms that as a switch, the main function of HIFs is to regulate the production of immune response factors in early infection. So activation of the switch may be an effective method for infectious disease treatment. As expected, the treatment groups improved the expression of HIFs compared with the control group, and the effects of the three strains are different. The strain1 of L. lactis had a stronger induction on HIF genes than strain2 and strain3, and it might be applied as a potential activation of HIF genes for disease treatment. So, adding befitting L. lactis maybe a well method to activate the HIF genes to protect them from mycobacterial infection.

Evolution of metazoan oxygen-sensing involved a conserved divergence of VHL affinity for HIF1α and HIF2α

Duplication of ancestral hypoxia-inducible factor (HIF)α coincided with the evolution of vertebrate species. Paralogs HIF1α and HIF2α are the most well-known factors for modulating the cellular transcriptional profile following hypoxia. However, how the processes of natural selection acted upon the coding region of these two genes to optimize the cellular response to hypoxia during evolution remains unclear. A key negative regulator of HIFα is von Hippel-Lindau (VHL) tumour suppressor protein. Here we show that evolutionarily-relevant substitutions can modulate a secondary contact between HIF1α Met561 and VHL Phe91. Notably, HIF1α binds more tightly than HIF2α to VHL due to a conserved Met to Thr substitution observed in the vertebrate lineage. Similarly, substitution of VHL Phe91 with Tyr, as seen in invertebrate species, decreases VHL affinity for both HIF1α and HIF2α. We propose that vertebrate evolution involved a more complex hypoxia response with fine-tuned divergence of VHL affinity for HIF1α and HIF2α.

Paralogs HIF1α and HIF2α are important modulators regulating cellular transcriptional profile following hypoxia. Here, the authors investigate evolutionary substitutions that fine tune the interaction between HIFα and their regulator VHL in the vertebrate and invertebrate lineages.

Divergence, evolution and adaptation in ray-finned fish genomes

With the rapid development of next-generation sequencing technologies and bioinformatics, over 50 ray-finned fish genomes by far have been sequenced with high quality. The genomic work provides abundant genetic resources for deep understanding of divergence, evolution and adaptation in the fish genomes. They are also instructive for identification of candidate genes for functional verification, molecular breeding, and development of novel marine drugs. As an example of other omics data, the Fish-T1K project generated a big database of fish transcriptomes to integrate with these published fish genomes for potential applications. In this review, we highlight the above-mentioned recent investigations and core topics on the ray-finned fish genome research, with a main goal to obtain a deeper understanding of fish biology for theoretical and practical applications.

Convergent evolution misled taxonomy in schizothoracine fishes (Cypriniformes: Cyprinidae)

Highly specialized grade (HSG; genera Gymnocypris, Oxygymnocypris, Schizopygopsis, Platypharodon and Chuanchia) of the Schizothoracinae (Cypriniformes: Cyprinidae) are endemic to the Qinghai-Tibet Plateau (QTP). Previously, two distinct ecomorphs were recognized according to trophic traits. One was a limnetic omnivore with normal lower jaw morphology, terminal mouth, and moderate or dense gill rakers, mostly inhabiting in open water of lakes, including Gymnocypris and Oxygymnocypris. Another was a benthic feeder with inferior mouth, sparse gill rakers and sharp horny sheath on the lower jaw for scraping of attached prey off hard substrates, including Schizopygopsis, Platypharodon and Chuanchia. However, traditional taxonomy of HSG based on these trophic traits presented extensive conflicts with the molecular studies in recent years. The possible cause could be convergent evolution in morphology, retention of ancestral polymorphisms or mitochondrial introgression, but these hypotheses could not be assessed due to incomplete taxon sampling and only mitochondrial data employed in previous works. Here, we conducted the most comprehensive molecular analysis on HSG fishes to date, using four mitochondrial loci and 152,464 genome-wide SNPs, and including 21 of 24 putative species and one undescribed Schizopygopsis species. Both SNP and mtDNA trees confirmed extensive paraphyly of genera Gymnocypris and Schizopygopsis, where species often were clustered together by watershed instead of by genus. Basal split into the north clade B and the south clade C (ca. 3.03 Ma) approximately by the Tanggula-Tanitawen Mountains in SLAF tree coincided with a violent uplift of the QTP during the phase A of 'Qingzang movement' (ca. 3.6 Ma). Ancestral state reconstruction of the trophic ecomorph showed that the limnetic omnivore ecomorph had evolved repeatedly in clade B and C. Furthermore, we presented a striking case of convergent evolution between two 'subspecies' Gymnocypris chui chui and G. chui longimandibularis, which had diverged as early as two million years ago (ca. 2.42 Ma). Ecological analyses revealed that similar food utilization, particularly in zooplankton, was the main underlying driving force. This work showed an example of taxonomy with the most extensive errors at the genus/species levels due to convergent evolution and suggested that trophic traits could be misleading in fish taxonomy. Therefore, we propose a major generic revision for HSG species.

Evolution for extreme living: variation in mitochondrial cytochrome c oxidase genes correlated with elevation in pikas (genus Ochotona)

The genus Ochotona (pikas) is a clade of cold-tolerant lagomorphs that includes many high-elevation species. Pikas offer a unique opportunity to study adaptations and potential limitations of an ecologically important mammal to high-elevation hypoxia. We analyzed the evolution of 3 mitochondrial genes encoding the catalytic core of cytochrome c oxidase (COX) in 10 pika species occupying elevations from sea level to 5000 m. COX is an enzyme highly reliant on oxygen and essential for cell function. One amino acid property, the equilibrium constant (ionization of COOH), was found to be under selection in the overall protein complex. We observed a strong relationship between the net value change in this property and the elevation each species occupies, with higher-elevation species having potentially more efficient proteins. We also found evidence of selection in low-elevation species for potentially less efficient COX, perhaps trading efficiency for heat production in the absence of hypoxia. Our results suggest that different pika species may have evolved elevation-specific COX proteins, specialization that may indicate limitations in their ability to shift their elevational ranges in response to future climate change.

Transcriptome Analysis Provides Insights Into the Adaptive Responses to Hypoxia of a Schizothoracine Fish (Gymnocypris eckloni)

The schizothoracine fish endemic to the Qinghai-Tibetan Plateau are comparatively well adapted to aquatic environments with low oxygen partial pressures. However, few studies have used transcriptomic profiling to investigate the adaptive responses of schizothoracine fish tissues to hypoxic stress. This study compared the transcriptomes of Gymnocypris eckloni subjected to 72 h of hypoxia (Dissolved oxygen, DO = 3.0 ± 0.1 mg/L) to those of G. eckloni under normoxia (DO = 8.4 ± 0.1 mg/L). To identify the potential genes and pathways activated in response to hypoxic stress, we collected muscle, liver, brain, heart, and blood samples from normoxic and hypoxic fish for RNA-Seq analysis. We annotated 337,481 gene fragments. Of these, 462 were differentially expressed in the hypoxic fish as compared to the normoxic fish. Under hypoxia, the transcriptomic profiles of the tissues differed, with muscle the most strongly affected by hypoxia. Our data indicated that G. eckloni underwent adaptive changes in gene expression in response to hypoxia. Several strategies used by G. eckloni to cope with hypoxia were similar to those used by other fish, including a switch from aerobic oxidation to anaerobic glycolysis and the suppression of major energy-requiring processes. However, G. eckloni used an additional distinct strategy to survive hypoxic environments: a strengthening of the antioxidant system and minimization of ischemic injury. Here, we identified several pathways and related genes involved in the hypoxic response of the schizothoracine fish. This study provides insights into the mechanisms used by schizothoracine fish to adapt to hypoxic environments.

Hypoxia-inducible transcription factors in fish: expression, function and interconnection with the circadian clock

The hypoxia-inducible transcription factors are key regulators for the physiological response to low oxygen availability. In vertebrates, typically three Hif-α isoforms, Hif-1α, Hif-2α and Hif-3α, are expressed, each of which, together with Hif-1β, may form a functional heterodimer under hypoxic conditions, controlling expression of hundreds of genes. A teleost-specific whole-genome duplication complicates the analysis of isoform-specific functions in fish, but recent studies suggest that the existence of paralogues of a specific isoform opens up the possibility for a subfunctionalization. In contrast to during development inside the uterus, fish eggs are freely accessible and studies analyzing Hif expression in fish embryos during development have revealed that Hif proteins are not only controlling the hypoxic response, but are also crucial for proper development and organ differentiation. Significant advances have been made in our knowledge about tissue-specific functions of Hif proteins, especially with respect to gill or gonadal tissue. The hypoxia signalling pathway is known to be tightly and mutually intertwined with the circadian clock in zebrafish and mammals. Recently, a mechanistic explanation for the hypoxia-induced dampening of the transcriptional clock was detected in zebrafish, including also metabolically induced alterations of cellular redox signalling. In turn, MAP kinase-mediated H2O2 signalling modulates the temporal expression of Hif-1α protein, similar to the redox regulation of the circadian clock itself. Once again, the zebrafish has emerged as an excellent model organism with which to explore these specific functional aspects of basic eukaryotic cell biology.

Molecular evolution of myoglobin in the Tibetan Plateau endemic schizothoracine fish (Cyprinidae, Teleostei) and tissue-specific expression changes under hypoxia

Myoglobin (Mb) is an oxygen-binding hemoprotein that was once thought to be exclusively expressed in oxidative myocytes of skeletal and cardiac muscle where it serves in oxygen storage and facilitates intracellular oxygen diffusion. In this study, we cloned the coding sequence of the Mb gene from four species, representing three groups, of the schizothoracine fish endemic to the Qinghai-Tibetan Plateau (QTP), then conducted molecular evolution analyses. We also investigated tissue expression patterns of Mb and the expression response to moderate and severe hypoxia at the mRNA and protein levels in a representative of the highly specialized schizothoracine fish species, Schizopygopsis pylzovi. Molecular evolution analyses showed that Mb from the highly specialized schizothoracine fish have undergone positive selection and one positively selected residue (81L) was identified, which is located in the F helix, close to or in contact with the heme. We present tentative evidence that the Mb duplication event occurred in the ancestor of the schizothoracine and Cyprininae fish (common carp and goldfish), and that the Mb2 paralog was subsequently lost in the schizothoracine fish. In S. pylzovi, Mb mRNA is expressed in various tissues with the exception of the intestine and gill, but all such tissues, including the liver, muscle, kidney, brain, eye, and skin, expressed very low levels of Mb mRNA (< 8.0%) relative to that of the heart. The trace levels of Mb expression in non-muscle tissues are perhaps the major reason why non-muscle Mb remained undiscovered for so long. The expression response of the Mb gene to hypoxia at the mRNA and protein levels was strikingly different in S. pylzovi compared to that found in the common carp, medaka, zebrafish, and goldfish, suggesting that the hypoxia response of Mb in fish may be species and tissue-specific. Notably, severe hypoxia induced significant expression of Mb at the mRNA and protein levels in the S. pylzovi heart, which suggests Mb has a major role in the supply of oxygen to the heart of Tibetan Plateau fish.

Adaptive evolution of interferon regulatory factors is not correlated with body scale reduction or loss in schizothoracine fish

Body scales in teleost fish are the first line of defense in protecting the fish from invading pathogens. However, the relationship between the adaptive evolution of immune-related genes and changes in the body scale-covering of fish has not been previously studied. Schizothoracine fish, characterized by progressive reduction of body scales from the primitive to the highly specialized species, are a good group to investigate this relationship. We obtained 11 IRF genes (IRF1-11) from 14 of schizothoracine fish representing primitive, specialized, and highly specialized species, of which seven IRF genes (IRF2, IRF3, IRF5, IRF6, IRF7, IRF8 and IRF9) contained the complete CDS. Sequence analysis demonstrated the deletion or insertion of 4-7 amino acids in IRF2, IRF3, IRF6 and IRF9, which seems to be a common phenomenon in the schizothoracine fish. Selection pressure analysis supported the hypothesis that positive selection has been driving the rapid evolution of IRFs in specific lineages of the schizothoracine fish. This, however, is not correlated with body scale reduction or loss in the evolution of these IRFs. Remarkably, the deletion or insertion found in IRF protein sequences presented a regular pattern corresponding to the scale-covering changes in schizothoracine fish. Our study provides evidence for positive selection in the IRF family, contributing to a better understanding of the adaptive evolution of immune-related genes in schizothoracine fish in response to environmental changes of the Qinghai-Tibetan Plateau.

Genomic signature of highland adaptation in fish: a case study in Tibetan Schizothoracinae species

Background

Genome-wide studies on highland adaptation mechanism in terrestrial animal have been widely reported with few available for aquatic animals. Tibetan Schizothoracinae species are ideal model systems to study speciation and adaptation of fish. The Schizothoracine fish, Gymnocypris przewalskii ganzihonensis had underwent the ecological niche shift from salt water to freshwater, and also experienced a recent split from Gymnocypris przewalskii przewalskii. In addition, G. p. ganzihonensis inhabited harsh aquatic environment including low temperature and hypoxia as well as other Schizothoracinae species, its genetic mechanism of highland adaptation have yet to be determined.

Results

Our study used comparative genomic analysis based on the transcriptomic data of G. p. ganzihonensis and other four fish genome datasets to investigate the genetic basis of highland adaptation in Schizothoracine fish. We found that Schizothoracine fish lineage on the terminal branch had an elevated dN/dS ratio than its ancestral branch. A total of 202 gene ontology (GO) categories involved into transport, energy metabolism and immune response had accelerated evolutionary rates than zebrafish. Interestingly, we also identified 162 genes showing signature of positive selection (PSG) involved into energy metabolism, transport and immune response in G. p. ganzihonesis. While, we failed to find any PSG related to hypoxia response as previous studies.

Conclusions

Comparative genomic analysis based on G. p. ganzihonensis transcriptome data revealed significant genomic signature of accelerated evolution ongoing within Tibetan Schizothoracinae species lineage. Molecular evolution analysis suggested that genes involved in energy metabolism, transport and immune response functions in Schizothoracine fish underwent positive selection, especially in innate immunity including toll-like receptor signaling pathway genes. Taken together, our result as a case study in Schizothoracinae species provides novel insights in understanding the aquatic animal adaptation to extreme environment on the Tibetan Plateau, and also provides valuable genomic resource for further functional verification studies.

Electronic supplementary material

The online version of this article (10.1186/s12864-017-4352-8) contains supplementary material, which is available to authorized users.

Integrated mRNA and microRNA transcriptome analyses reveal regulation of thermal acclimation in Gymnocypris przewalskii: A case study in Tibetan Schizothoracine fish

Environmental acclimation is important episode in wildlife occupation of the high-altitude Tibetan Plateau (TP). Transcriptome-wide studies on thermal acclimation mechanism in fish species are rarely revealed in Tibetan Plateau fish at high altitude. Thus, we used mRNA and miRNA transcriptome sequencing to investigate regulation of thermal acclimation in larval Tibetan naked carp, Gymnocypris przewalskii. We first remodeled the regulation network of mRNA and miRNA in thermal acclimation, and then identified differential expression of miRNAs and target mRNAs enriched in metabolic and digestive pathways. Interestingly, we identified two candidate genes contributed to normal skeletal development. The altered expression of these gene groups could potentially be associated with the developmental issues of deformity and induced larval death. Our results have three important implications: first, these findings provide strong evidences to support our hypothesis that G. przewalskii possess ability to build heat-tolerance against the controversial issue. Second, this study shows that transcriptional and post-transcriptional regulations are extensively involved in thermal acclimation. Third, the integrated mRNA and microRNA transcriptome analyses provide a large number of valuable genetic resources for future studies on environmental stress response in G. przewalskii and as a case study in Tibetan Schizothoracine fish.

Identification, molecular evolution of toll-like receptors in a Tibetan schizothoracine fish (Gymnocypris eckloni) and their expression profiles in response to acute hypoxia

Hypoxia plays an important role in regulating a variety of physiological responses as well as in pathological situations, but to date the roles of Toll-like receptors (TLRs) in fish in response to hypoxia are still poorly understood. Here, we sequenced the transcriptome of G. eckloni and identified the members of TLR family by scanning transcriptome, and then investigated the expression profiles of a complete set of TLRs in G. eckloni in response to acute hypoxia (4 h at DO = 0.3 ± 0.1 mg/L). The de novo-assembled transcriptome consisted of a total of 162,235 transcripts, further clustered into 110,231 unigenes. Based on the transcriptome, a total of 18 TLRs were identified in G. eckloni, and of them three TLRs (TLR5, TLR8 and TLR22) possessed two distinct paralogous genes. The duplicated genes of TLR22 were discovered for the first time in cyprinid fish, but did not origin from a recent duplication event. Of them TLR22b may be specific for schizothoracine fish, at least for G. eckloni. Phylogenetic analysis supported the classification of TLRs into six families as in other vertebrates but was partly different from the previous study. The sliding window analysis showed strong signals of positive selection in TLR2, TLR 4, TLR 5a, TLR 7, TLR 19, TLR 20, TLR 21, TLR 22a and TLR 22b, but most codons under positive selection were located in the putative LRR regions. The mRNA expression of most TLRs in head kidney, spleen and gill decreased significantly or remained unchanged under acute hypoxia, whereas acute hypoxia increased expressions of TLR2 and TLR3 in head kidney, of TLR8a, TLR12 and TLR19 in spleen, and of TLR1 in gill, suggesting tissues-specific expressions of TLRs play important roles in mediating innate immune responses for host defense against tissue damages or physiological changes induced by hypoxia.

Genome-wide identification of genes probably relevant to the adaptation of schizothoracins (Teleostei: Cypriniformes) to the uplift of the Qinghai-Tibet Plateau

BACKGROUND

Molecular adaptation to the severe environments present during the uplift of the Qinghai-Tibet Plateau has attracted the attention of researchers. The divergence of the three specialization groups of schizothoracins (Primitive, Specialized and Highly Specialized) may correspond to the three phases of plateau uplift. Based on the transcripts of representative species of the three specialized groups and an outgroup, genes in schizothoracins that may have played important roles during the adaptation to new environments were investigated.

RESULTS

The contigs of Gymnodiptychus dybowskii and Schizothorax pseudaksaiensis were compared with those of Gymnocypris przewalskii ganzihonensis and the outgroup Sinocyclocheilus angustiporus, and 5,894 ortholog groups with an alignment length longer than 90 nt after deleting gaps were retained. Evolutionary analyses indicated that the average evolutionary rate of the branch leading to the Specialized group was faster than that of the branch leading to the Highly Specialized group. Moreover, the numbers of gene categories in which more than half of the genes evolved faster than the average values of the genome were 117 and 15 along the branches leading to the Specialized and Highly Specialized groups, respectively. A total of 40, 36, and 55 genes were likely subject to positive selection along the branches leading to the Primitive, Specialized and Highly Specialized groups, respectively, and many of these genes are likely relevant to adaptation to the cold temperatures, low oxygen concentrations, and strong ultraviolet radiation that result from elevation.

CONCLUSIONS

By selecting representative species of the three groups of schizothoracins and applying next-generation sequencing technology, several candidate genes corresponding to adaptation to the three phases of plateau uplift were identified. Some of the genes identified in this report that were likely subject to positive selection are good candidates for subsequent evolutionary and functional analyses of adaptation to high altitude.

Genetic Adaptation of Schizothoracine Fish to the Phased Uplifting of the Qinghai-Tibetan Plateau

Many species of Schizothoracine, a subfamily of Cyprinidae, are highly endemic to the Qinghai–Tibetan Plateau (QTP). To characterize the adaptive changes associated with the Schizothoracine expansion at high altitudes, we sequenced tissue transcriptomes of two highland and two subhighland Schizothoracines and analyzed gene evolution patterns by comparing with lowland cyprinids. Phylogenetic tree reconstruction and divergence time estimation indicated that the common ancestor of Schizothoracine fish lived ∼32.7 million years ago (MYA), coinciding with the timing of the first phase of QTP uplifting. Both high- and subhigh-Schizothoracines demonstrated elevated dN/dS ratios in the protein-coding genes compared to lowland cyprinids, from which some biological processes implicated in altitude adaptation were commonly identified. On the other hand, the highland and subhighland lineages presented drastically divergent landscapes of positively selected genes (PSGs), enriched with very different gene ontology (GO) profiles, including those in “sensory organ morphogenesis,” “regulation of protein ubiquitination,” “blood circulation,” and “blood vessel development.” These results indicated different selection pressures imposed on the highland and subhighland lineages of the Schizothoracine subfamily, with a higher number of genes in the high-altitude species involved in adaptations such as sensory perception, blood circulation, and protein metabolism. Our study indicated divergent genetic adaptations in the aquatic species facing the phased uplifting of QTP.

The microRNA repertoire of Tibetan naked carp Gymnocypris przewalskii: A case study in Schizothoracinae fish on the Tibetan Plateau

Tibetan naked carp Gymnocypris przewalskii is an ideal model system to study highland adaptation of fish, because it evolved specific genetic and phenotypic characteristics to adapt to chronic cold and alkaline environments in Lake Qinghai. MicroRNAs (miRNAs) are small noncoding RNAs that regulating gene expression post-transcriptionally in a wide range of biological processes. In this study, we focus on the role of miRNAs in adaptation of G. przewalskii to extreme conditions in Lake Qinghai. We generate the first miRNAome of G. przewalskii in Schizothoracinae fish. Using several genomic resources, we inferred 341 conserved miRNAs belonged to 152 miRNA families and 43 novel miRNAs in G. przewalskii, and also identified 15 teleost-specific miRNAs. Using a large scale of conserved miRNAs, we constructed a high-confidence phylogenetic tree between teleost and mammals than mitochondria and nuclear genes. In addition, we found that several miRNA family (e.g. miR-10 and let-7) members highly expressed in G. przewalskii, which may function in multiple biological processes. Finally, we predicted a total of 34,258 miRNA targets genes. Conserved miRNAs target genes participating in signal transduction, cell differentiation and biosynthetic process, and showed signature of functional constraint. While novel miRNAs in a species displayed species-specific targets and involved in ion binding, transport and oxidoreductase activity, may affect the expression patterns of targets with signature of gene family expansion or positive selection under extreme environment. Taken together, this study demonstrated that miRNAs may involve into roles of adaptation of G. przewalskii to highland aquatic environment, and also provide insights into miRNA regulatory network in Schizothoracinae fish as a case study.

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.

Positive Darwinian selection within interferon regulatory factor genes of Gymnocypris przewalskii (Cyprinidae) on the Tibetan Plateau

Tibetan Plateau (TP) had experienced phased uplift, resulting in inhospitable environment of low temperature, hypoxia and high ultraviolet radiation for Tibetan wildlife. Many organisms can well adapt to TP, it is of ecological and evolutionary interest to untangle how organisms adapt to extreme environment on TP through evolution. Previous studies mainly focused on hypoxia and metabolism related genes, but we know little about the evolutionary history of immune genes in Tibetan wildlife. In this study, we first identified 10 interferon regulatory factor (IRF) genes from Tibetan naked carp Gymnocypris przewalskii. Within this gene family, IRF3, IRF5, IRF7 and IRF8 contained positive selection sites. Evidences indicated that positive selection may lead to IRF genes functional alternations, presumably driving genes towards adaptation to the environmental changes. Taken together, our results suggested 4 candidate genes as interesting targets for further experimental confirmation of their functional variations and contributions to high altitude adaptation in Tibet fish.

Analysis of the erythropoietin of a Tibetan Plateau schizothoracine fish (Gymnocypris dobula) reveals enhanced cytoprotection function in hypoxic environments

BACKGROUND

Erythropoietin (EPO) is a glycoprotein hormone that plays a principal regulatory role in erythropoiesis and initiates cell homeostatic responses to environmental challenges. The Qinghai-Tibet Plateau is a natural laboratory for hypoxia adaptation. Gymnocypris dobula is a highly specialized plateau schizothoracine fish that is restricted to > 4500 m high-altitude freshwater rivers and ponds in the Qinghai-Tibet Plateau. The role of EPO in the adaptation of schizothoracine fish to hypoxia is unknown.

RESULTS

The EPO and EPO receptor genes from G. dobula and four other schizothoracine fish from various altitudinal habitats were characterized. Schizothoracine EPOs are predicted to possess 2-3 N-glycosylation (NGS) sites, 4-5 casein kinase II phosphorylation (CK2) sites, 1-2 protein kinase C (PKC) phosphorylation sites, and four conserved cysteine residues within four helical domains, with variations in the numbers of NGS and CK2 sites in G. dobula. PAML analysis indicated a d N/d S value (ω) = 1.112 in the G. dobula lineage, and a few amino acids potentially under lineage-specific positive selection were detected within the G. dobula EPO. Similarly, EPO receptors of the two high-altitude schizothoracines (G. dobula and Ptychobarbus kaznakovi), were found to be statistically on the border of positive selection using the branch-site model (P-value = 0.096), and some amino acids located in the ligand-binding domain and the fibronectin type III domain were identified as potentially positive selection sites. Tissue EPO expression profiling based on transcriptome sequencing of three schizothoracines (G. dobula, Schizothorax nukiangensis Tsao, and Schizothorax prenanti) showed significant upregulation of EPO expression in the brain and less significantly in the gill of G. dobula. The elevated expression together with the rapid evolution of the EPO gene in G. dobula suggested a possible role for EPO in adaptation to hypoxia. To test this hypothesis, Gd-EPO and Sp-EPO were cloned into an expression vector and transfected into the cultured cell line 293 T. Significantly higher cell viability was observed in cells transfected with Gd-EPO than cells harboring Sp-EPO when challenged by hypoxia.

CONCLUSION

The deduced EPO proteins of the schizothoracine fish contain characteristic structures and important domains similar to EPOs from other taxa. The presence of potentially positive selection sites in both EPO and EPOR in G. dobula suggest possible adaptive evolution in the ligand-receptor binding activity of the EPO signaling cascade in G. dobula. Functional study indicated that the EPO from high-altitude schizothoracine species demonstrated features of hypoxic adaptation by reducing toxic effects or improving cell survival when expressed in cultured cells, providing evidence of molecular adaptation to hypoxic conditions in the Qinghai-Tibet Plateau.

Hypoxia-Inducible Factor α and Hif-prolyl Hydroxylase Characterization and Gene Expression in Short-Time Air-Exposed Mytilus galloprovincialis

Aquatic organisms experience environmental hypoxia as a result of eutrophication and naturally occurring tidal cycles. Mytilus galloprovincialis, being an anoxic/hypoxic-tolerant bivalve, provides an excellent model to investigate the molecular mechanisms regulating oxygen sensing. Across the animal kingdom, inadequacy in oxygen supply is signalled predominantly by hypoxia-inducible factors (HIF) and Hif-prolyl hydroxylases (PHD). In this study, hif-α 5'-end and partial phd mRNA sequences from M. galloprovincialis were obtained. Phylogenetic and molecular characterization of both HIF-α and PHD putative proteins showed shared key features with the respective orthologues from animals strongly suggesting their crucial involvement in the highly conserved oxygen sensing pathway. Both transcripts displayed a tissue-specific distribution with prominent expression in gills. Quantitative gene expression analysis of hif-α and phd mRNAs from gills of M. galloprovincialis demonstrated that both these key sensors are transcriptionally modulated by oxygen availability during the short-time air exposure and subsequent re-oxygenation treatments proving that they are critical players of oxygen-sensing mechanisms in mussels. Remarkably, hif-α gene expression showed a prompt and transient response suggesting the precocious implication of this transcription factor in the early phase of the adaptive response to hypoxia in Mytilus. HIF-α and PHD proteins were modulated in a time-dependent manner with trends comparable to mRNA expression patterns, thus suggesting a central role of their transcriptional regulation in the hypoxia tolerance strategies in marine bivalves. These results provide molecular information about the effects of oxygen deficiency and identify hypoxia-responsive biomarker genes in mussels applicable in ecotoxicological studies of natural marine areas.

Polyphyletic origins of schizothoracine fish (Cyprinidae, Osteichthyes) and adaptive evolution in their mitochondrial genomes

The schizothoracine fish, also called snow trout, are members of the Cyprinidae, and are the most diversified teleost fish in the Qinghai-Tibetan Plateau (QTP). Clarifying the evolutionary history of the schizothoracine fish is therefore important for better understanding the biodiversity of the QTP. Although morphological and molecular phylogenetic studies have supported the monophyly of the Schizothoracinae, a recent molecular phylogenetic study based on the mitochondrial genome questioned the monophyly of this taxon. However, the phylogenetic analysis of that study was on the basis of only three schizothoracine species, and the support values were low. In this report, we inferred the phylogenetic tree on the basis of mitochondrial genome data including 21 schizothoracine species and five closely related species, and the polyphyletic origins of the Schizothoracinae were strongly supported. The tree further suggests that the Schizothoracinae consists of two clades, namely the "morphologically specialized clade" and the "morphologically primitive clade", and that these two clades migrated independently of each other to the QTP and adapted to high altitude. We also detected in their mitochondrial genomes strong signals of positive selection, which probably represent evidence of high-altitude adaptation. In the case of the morphologically specialized clade, positive selection mainly occurred during the Late Paleocene to the Early Oligocene. Its migration also seems to have occurred in the Early Eocene, and this timing is consistent with the drastic uplifting of the QTP. On the other hand, positive selection in the morphologically primitive clade has mainly occurred since the Late Miocene. Because its members are thought to have migrated to the QTP recently, it is possible that they are now undergoing high-altitude adaptation.

Analysis of hypoxia-inducible factor alpha polyploidization reveals adaptation to Tibetan Plateau in the evolution of schizothoracine fish

BACKGROUND

Hypoxia-inducible factor (HIF) is a master regulator that mediates major changes in gene expression under hypoxic conditions. Though HIF family has been identified in many organisms, little is known about this family in schizothoracine fish.

RESULTS

Duplicated hif-α (hif-1αA, hif-1αB, hif-2αA, and hif-2αB) genes were identified in schizothoracine fish. All the deduced HIF-α proteins contain the main domains (bHLH-PAS, ODDD, and TAD), also found in humans. Evidence suggests a Cyprinidae-specific deletion, specifically, a conserved proline hydroxylation motif LxxLAP, in the NODD domain of schizothoracine fish HIF-1αA. In addition, a schizothoracine-specific mutation was observed in the CODD domain of the specialized and highly specialized schizothoracine fish HIF-1αB, which is the proline hydroxylation motif mutated into PxxLAP. Standard and stochastic branch-site codon model analysis indicated that only HIF-1αB has undergone positive selection, which may have led to changes in function. To confirm this hypothesis, HIF-αs tagged with Myc were transfected into HEK 293 T cells. Each HIF-1αB was found to significantly upregulate luciferase activity under normoxic and hypoxic conditions, which indicated that the HIF-1αB protein was more stable than other HIF-αs.

CONCLUSIONS

All deduced HIF-α proteins of schizothoracine fish contain important domains, like their mammalian counterparts, and each HIF-α is shorter than that of human. Our experiments reveal that teleost-specific duplicated hif-α genes played different roles under hypoxic conditions, and HIF-1αB may be the most important regulator in the adaptation of schizothoracine fish to the environment of the Tibetan Plateau.