Physiological, metabolomic, and transcriptomic reveal metabolic pathway alterations in Gymnocypris przewalskii due to cold exposure.

Teleost fish have evolved various adaptations that allow them to tolerate cold water conditions. However, the underlying mechanism of this adaptation is poorly understood in Tibetan Plateau fish. RNA-seq combined with liquid chromatography‒mass spectrometry (LC‒MS/MS) metabolomics was used to investigate the physiological responses of a Tibetan Plateau-specific teleost, Gymnocypris przewalskii, under cold conditions. The 8-month G. przewalskii juvenile fish were exposed to cold (4 ℃, cold acclimation, CA) and warm (17 ℃, normal temperature, NT) temperature water for 15 days. Then, the transcript profiles of eight tissues, including the brain, gill, heart, intestine, hepatopancreas, kidney, muscle, and skin, were evaluated by transcriptome sequencing. The metabolites of the intestine, hepatopancreas, and muscle were identified by LC‒MS/MS. A total of 5,745 differentially expressed genes (DEGs) were obtained in the CA group. The key DEGs were annotated using Gene Ontology and Kyoto Encyclopedia of Genes and Genomes analysis. The DEGs from the eight tissues were significantly enriched in spliceosome pathways, indicating that activated alternative splicing is a critical biological process that occurs in the tissues to help fish cope with cold stress. Additionally, 82, 97, and 66 differentially expressed metabolites were identified in the intestine, hepatopancreas, and muscle, respectively. Glutathione metabolism was the only overlapping significant pathway between the transcriptome and metabolome analyses in these three tissues, indicating that an activated antioxidative process was triggered during cold stress. In combination with the multitissue transcriptome and metabolome, we established a physiology-gene‒metabolite interaction network related to energy metabolism during cold stress and found that gluconeogenesis and long-chain fatty acid metabolism played critical roles in glucose homeostasis and energy supply.

Fatty acid metabolism and antioxidant capacity in Gymnocypris przewalskii (Kessler, 1876) response to thermal stress.

The Qinghai-Tibet Plateau is undergoing a wet-warming transition, which could affect the survival of the native fish. However, the tolerance and physiological response to thermal stress is rarely studied in Gymnocypris przewalskii, a rare native fish in the Tibetan plateau. In this study, first, we detected the thermal tolerance of five groups of six-month G. przewalskii which acclimated at 8, 12, 16, 20, and 24 °C for two weeks, respectively, by critical thermal methodology. Then, through heat challenge, we detected the metabolites, key enzyme activities, and gene expressions involved in metabolism and antioxidant in the hepatopancreas when the temperatures increased from 16 °C to 18, 20, 22, 24, 26, and 28 °C for 12 h, respectively. The results showed that although the fish are sensitive to high temperatures, the quick acclimation at mild high temperatures could significantly improve the tolerance to acute high-temperature stress in juvenile G. przewalskii. During the heat challenge study, blood glucose significantly increased at heat stress (P < 0.05). At the same time, total cholesterol (TC), triglyceride (TG), and free fatty acid (FFA) significantly decreased when the temperature rose continuously to 20 °C. Metabolic enzyme activities of carnitine palmityl transferase I (CPT-Ⅰ), acetyl-CoA carboxylase (ACC), and fatty acid synthase (FAS) significantly decreased at 20 °C (P < 0.05). Superoxide dismutase (SOD) and antioxidant capacity (T-AOC) significantly increased at 20 °C (P < 0.05). The relative transcript levels of genes involved in antioxidant and glycolysis/gluconeogenesis were markedly higher than the control at 20-26 °C (P < 0.05). The genes involved in fatty acid biosynthesis or metabolism showed different expression patterns under heat stress. Heat shock protein 70 (Hsp70) and Hsp90 were significantly higher than the control at 18 °C and 26 °C, respectively. These results confirmed the prediction that G. przewalskii is sensitive to high temperatures, so conservation efforts should pay more attention to the warming damage.

Population genomics analysis to identify ion and water transporter genes involved in the adaptation of Tibetan naked carps to brackish water.

Understanding how evolutionary processes shape the genetic variations and influence the response of species to environmental alterations is critical for biodiversity conservation and molecular breeding. Gymnocypris przewalskii przewalskii is the only known cyprinid fish that dwells in the brackish water of Lake Qinghai on the Qinghai-Tibetan Plateau. To reveal the genetic basis of its adaptation to high salinity and alkalinity, whole-genome sequencing was performed in G. p. przewalskii and its freshwater relatives Gymnocypris eckloni and Gymnocypris przewalskii ganzihonensis. Compared with freshwater species, lower genetic diversity and higher linkage disequilibrium were observed in G. p. przewalskii. Selective sweep analysis identified 424 core-selective genes enriched in transport activities. Transfection analysis showed that genetic changes in the positively selected gene aquaporin 3 (AQP3) improved cell viability after salt treatment, suggesting its involvement in brackish water adaptation. Our analysis indicates that ion and water transporter genes experienced intensive selection, which might have contributed to the maintenance of high osmolality and ion content in G. p. przewalskii. The current study identified key molecules involved in the adaptation of fish to brackish water, providing valuable genomic resources for the molecular breeding of salt-tolerant fish.

Genome-wide identification and expression of the peroxisome proliferator-activated receptor gene family in the Tibetan highland fish Gymnocypris przewalskii.

Peroxisome proliferator-activated receptor (PPAR) plays an important role in the regulation of lipid metabolism and has been widely identified in diverse species. Gymnocypris przewalskii is a native fish of the Qinghai Tibetan Plateau that survives in a chronically cold environment. In the current study, we conducted genome-wide identification of PPAR genes, revealing the existence of seven PPARs in the G. przewalskii genome. Collinearity was observed between two copies of PPARαb and PPARγ in G. przewalskii, suggesting that the additional copy might be gained through whole genome duplication. Both phylogenetic and multiple sequence alignment analyses indicated that PPARs in G. przewalskii were conserved with teleosts. The cold treatment (10 °C and 4 °C) led to the developmental delay of G. przewalskii embryos. Continuous expression of PPARs was observed during the embryonic development of G. przewalskii under normal and cold conditions, with significantly different transcriptional patterns. These results indicated that PPARs participated in the embryonic development of G. przewalskii, and were involved in the cold response during development. The current study proposed a potential role of PPARs in the cold response in the embryonic development of G. przewalskii, which shed light on understanding cold adaptation in Tibetan highland fish.

Genome-wide characterization of the Elovl gene family in Gymnocypris przewalskii and their potential roles in adaptation to cold temperature.

The elongase of the very long-chain fatty acids (Elovls) gene family in fish has more diversity than in other vertebrates, which plays several critical roles in fatty acid synthesis and low-temperature stress adaptation. Gymnocypris przewalskii settles in plateau lakes with cold and resource-poor settings, and the evolution and function of Elovl genes in this fish are unknown. In the study, to identify the Elovl genes in G. przewalskii, the genome-wide identification and phylogenetic analysis of the gene members have been conducted with the expression profile of different tissues under cold stress. Fatty acid compositions, meanwhile, were detected in both the hepatopancreas and skeletal muscle during cold adaptation. A total of 21 Elovl members have been identified from the genome of G. przewalskii, belonging to Elovl1, Elovl2, Elovl4, Elovl5, Elovl6, Elovl7, and Elovl8 subgroups, with conserved ELO domain and four common motifs. Phylogenetic analysis revealed that subfamilies Elovl1 and Elovl7, Elov2, and Elovl5 have a closer genetic relationship, while the Elovl6 class was classed into an independent clade. Synteny analysis showed that whole-genome duplication, tandem duplicates, and gene conversion could drive the Elovls family expansion in G. przewalskii. The Ka/Ks and RELAX analysis showed distinguishing positive selection traces in ORF sequences of gpElovl2. Transcriptional data showed that different gpElovl subtypes exhibited a tissue-specific expression. Subtypes gpElovl1a, gpElovl2 and gpElovl6l were highly expressed induced by cold stress, as well as fatty acid metabolism-related genes, including Acyl-CoA synthetase long-chain gene (Ascl1a-1) and Stearyl-CoA desaturase gene (Scd1a-1). In addition, monounsaturated fatty acid (MUFA) and polyunsaturated fatty acid (PUFA) contents of the hepatopancreas and skeletal muscle were significantly increased under 15-day cold stress. These results provide a better understanding of fish Elovl genes and their roles in cold adaptation.

Physiological, morphological and transcriptomic responses of Tibetan naked carps (Gymnocypris przewalskii) to salinity variations.

Gymnocypris przewalskii is a native cyprinid fish that dwells in the Lake Qinghai with salinity of 12-13‰. It migrates annually to the freshwater rivers for spawning, experiencing the significant changes in salinity. In the present study, we performed the physiological, morphological and transcriptomic analyses to understand the osmoregulation in G. przewalskii. The physiological assay showed that the osmotic pressure of G. przewalskii was almost isosmotic to the brackish lake water. The low salinity reduced its ionic concentrations and osmotic pressure. The plasticity of gill microstructure was linked to the salinity variations, including the presence of mucus and intact tight junctions in brackish water and the development of the mitochondria-rich cells and the loosened tight junctions in freshwater. RNA-seq analysis identified 1926 differentially expressed genes, including 710 and 1216 down- and up-regulated genes in freshwater, which were enriched in ion transport, cell-cell adhesion, and mucus secretion. Genes in ion uptake were activated in low salinity, and mucus pathways and tight junction showed the higher transcription in brackish water. The isosmoticity between the body fluid and the environment suggested G. przewalskii was in the metabolic-saving condition in the brackish water. The decreased salinity disrupted this balance, which activated the ion uptake in freshwater to maintain osmotic homeostasis. The gill remodeling was involved in this process through the development of the mitochondria-rich cells to enhance ion uptake. The current finding provided insights into the potential mechanisms of G. przewalskii to cope with salinity alteration.

Effectiveness assessment of using riverine water eDNA to simultaneously monitor the riverine and riparian biodiversity information.

Both aquatic and terrestrial biodiversity information can be detected in riverine water environmental DNA (eDNA). However, the effectiveness of using riverine water eDNA to simultaneously monitor the riverine and terrestrial biodiversity information remains unidentified. Here, we proposed that the monitoring effectiveness could be approximated by the transportation effectiveness of land-to-river and upstream-to-downstream biodiversity information flows and described by three new indicators. Subsequently, we conducted a case study in a watershed on the Qinghai-Tibet Plateau. The results demonstrated that there was higher monitoring effectiveness on summer or autumn rainy days than in other seasons and weather conditions. The monitoring of the bacterial biodiversity information was more efficient than the monitoring of the eukaryotic biodiversity information. On summer rainy days, 43-76% of species information in riparian sites could be detected in adjacent riverine water eDNA samples, 92-99% of species information in riverine sites could be detected in a 1-km downstream eDNA sample, and half of dead bioinformation (the bioinformation labeling the biological material that lacked life activity and fertility) could be monitored 4-6 km downstream for eukaryotes and 13-19 km downstream for bacteria. The current study provided reference method and data for future monitoring projects design and for future monitoring results evaluation.

Riverbed Substrate Requirements for Natural Reproduction of Gymnocypris przewalskii.

Gymnocypris przewalskii (i.e., Qinghai Lake naked carp) is a migratory fish species that lives in highland brackish water. It is important to understand the abiotic environment required by this fish to reproduce naturally so that its habitat can be protected and the wild population can be conserved. Here, artificial simulation and spawning ground substrate transformation experiments were conducted to examine the riverbed substrate requirements for G. przewalskii to naturally reproduce. Using various techniques (in vitro markers, videography, and Ethovision XT behavior tracking), this study systematically investigated the riverbed substrate preferences of G. przewalskii as well as the characteristics and effectiveness of natural reproduction induced by pebble riverbed substrate. The findings can be summarized as follows: (1) the habitat preferences of G. przewalskii differed significantly between various riverbed substrate, with pebble substrate being preferred during natural reproduction, and sand substrate being preferred pre- and post-spawning, and (2) the natural reproduction of G. przewalskii was heavily reliant on pebble riverbed substrate. Specifically, pebble substrate significantly improved spawn quantity and fertilization rate. These findings provide scientific evidence for the improvement and restoration of G. przewalskii spawning grounds, and insights regarding the artificial bionic reproduction of G. przewalskii.

Characterization and complexity of transcriptome in Gymnocypris przewalskii using single-molecule long-read sequencing and RNA-seq.

The Tibetan Schizothoracinae fish Gymnocypris przewalskii has the ability to adapt to the extreme plateau environment, making it an ideal biological material for evolutionary biology research. However, the lack of well-annotated reference genomes has limited the study of the molecular genetics of G. przewalskii. To characterize its transcriptome features, we first used long-read sequencing technology in combination with RNA-seq for transcriptomic analysis. A total of 159,053 full-length (FL) transcripts were captured by Iso-Seq, having a mean length of 3,445 bp with N50 value of 4,348. Of all FL transcripts, 145,169 were well-annotated in the public database and 134,537 contained complete open reading frames. There were 4,149 pairs of alternative splicing events, of which three randomly selected were defined by RT-PCR and sequencing, and 13,293 long non-coding RNAs detected, based on all-vs.-all BLAST. A total of 118,185 perfect simple sequence repeats were identified from FL transcripts. The FL transcriptome might provide basis for further research of G. przewalskii.

Transcriptomic profiling reveals molecular regulation of seasonal reproduction in Tibetan highland fish, Gymnocypris przewalskii.

BACKGROUND: The Tibetan highland fish, Gymnocypris przewalskii, migrates from Lake Qinghai to its spawning grounds every summer. This seasonal reproduction is critically regulated by intrinsic and extrinsic signals. However, the molecular mechanisms that process environmental oscillations to initiate the seasonal mating are largely unknown. RESULTS: A transcriptomic analysis was conducted on the brain and gonad of male and female G. przewalskii in reproductive and nonreproductive seasons. We obtained 2034, 760, 1158 and 17,856 differentially expressed genes between the reproductively active and dormant female brain, male brain, ovary and testis. Among these genes, DIO2 was upregulated in the reproductively active brain and gonad of both males and females. Neuroactive ligand-receptor genes were activated in male and female brain. Functional enrichment analysis suggested that retinol metabolism was uniquely stimulated in reproductively active males. Genes involved in GnRH signaling and sex hormone synthesis exhibited higher expression levels in brain and gonad during the reproductive season. A co-expression network classified all the genes into 9 modules. The network pinpointed CDC42 as the hub gene that connected the pathways in responsible for modulating reproduction in G. przewalskii. Meanwhile, the sex pheromone receptor gene prostaglandin receptor was identified to link to multiple endocrine receptors, such as GnRHR2 in the network. CONCLUSIONS: The current study profiled transcriptomic variations between reproductively active and dormant fish, highlighting the potential regulatory mechanisms of seasonal reproduction in G. przewalskii. Our data suggested that the seasonal regulation of reproduction in G. przewalskii was controlled by the external stimulation of photoperiodic variations. The activated transcription of neuroendocrine and sex hormone synthesis genes contributed to seasonal reproduction regulation in G. przewalskii, which was presumably influenced by the increased day-length during the breeding season.

Studies of the cryopreservation condition of Gymnocypris przewalskii spermatozoa.

The endemic naked carp (Gymnocypris przewalskii Kessler, 1876) plays an important role in the maintenance of the distinctive ecological system of Lake Qinghai at 3.2km altitude on the Qinghai-Tibet Plateau in China. This study aimed to develop a cryopreservation protocol for Gymnocypris przewalskii spermatozoa. Semen was collected from mature individuals during migration and frozen using the liquid nitrogen vapor method. The influence of different cryoprotectants and three extenders on the post-thaw quality of the sperm was analyzed. The highest sperm motility rate and longest motility time after cryopreservation were achieved by combining Ringer's solution with 15% ethylene glycol (P<0.05). The fertilization rate of this cryopreserved semen was 15.26±4.54%. This study thus provides a valuable method for the cryopreservation of the sperm of this important endangered fish species.

Gymnocypris przewalskii decreases cytosolic carbonic anhydrase expression to compensate for respiratory alkalosis and osmoregulation in the saline-alkaline lake Qinghai.

Naked carp (Gymnocypris przewalskii), endemic to the saline-alkaline Lake Qinghai, have the capacity to tolerate combined high salinity and alkalinity, but migrate to spawn in freshwater rivers each year. In this study, the full-length cDNA of the cytosolic carbonic anhydrase c isoform of G. przewalskii (GpCAc) was amplified and sequenced; mRNA levels and enzyme activity of GpCAc and blood chemistry were evaluated to understand the compensatory responses as the naked carp returned to the saline-alkaline lake after spawning. We found that GpCAc had a total length of 1400 bp and encodes a peptide of 260 amino acids. Comparison of the deduced amino acid sequences and phylogenetic analysis showed that GpCAc was a member of the cytosolic carbonic anhydrase II-like c family. Cytosolic-carbonic-anhydrase-c-specific primers were used to analyze the tissue distribution of GpCAc mRNA expression. Expression of GpCAc mRNA was found in brain, gill, liver, kidney, gut, and muscle tissues, but primarily in the gill and posterior kidney; however, none was evident in red blood cells. Transferring fish from river water to lake water resulted in a respiratory alkalosis, osmolality, and ion rise in the blood, as well as significant decreases in the expression and enzyme activity of GpCAc in both the gill and kidney within 96 h. These results indicate that GpCAc may play an important role in the acclimation to both high salinity and carbonate alkalinity. Specifically, G. przewalskii decreases cytosolic carbonic anhydrase c expression to compensate for a respiratory alkalosis and to aid in osmoregulation during the transition from river to saline-alkaline lake.

Transcriptome-wide identification, molecular evolution and expression analysis of Toll-like receptor family in a Tibet fish, Gymnocypris przewalskii.

Toll-like receptors (TLR) are key components of innate immunity that play significant roles in immune defense against pathogens invasion. Recent frequent outbreaks of the "white spot disease" caused by parasitic infection in farmed Tibetan fishes had resulted in great economic losses. However, to our knowledge, the roles of TLRs in mediating immune response to parasitic infection in Tibetan fishes remain to be determined. Here, we performed data-mining on a widely-farmed Tibetan fish (Gymnocypris przewalskii or Gp) transcriptome to determine the genetic variation and expression pattern of TLRs. We totally obtained 14 GpTLRs and identified 5 with a complete coding sequence. Phylogenetic analysis verified their identities and supported the classification of TLRs into six families as in other vertebrates. The TLR family motifs, such as leucine rich repeat (LRR) and Toll/interleukin (IL)-1 receptor (TIR) domain, are conserved in GpTLR1-5. Selective pressure test demonstrated that all known GpTLRs are under purifying selection, except GpTLR4 underwent positive selection. Further, site model analysis suggested that 11 positively selected sites are found in LRR domain of GpTLR4. Three positively selected sites are located on outside surface of TLR4 3D structure, indicating that function of GpTLR4 may be affected. Tissue specific expression analysis showed all GpTLRs are present in gill, head-kidney and spleen but the relative abundance varied among tissues. In response to parasite Ichthyophthirius multifiliis infection, 5 GpTLR (GpTLR1, -2, -4, -9 and -20) expressions were induced. Intriguingly, GpTLR4 was significantly up-regulated in gills, while GpTLR19 and GpTLR21 unexpectedly showed no any change. In summary, these results revealed the first genomic resources of TLR family and several parasitic infection responsive TLRs in Tibetan fish. These findings provide key information for future studies aiming to understand the molecular mechanisms underlying the immune response to pathogen invasion in Tibetan fishes.

Characterization of two paralogous myostatin genes and evidence for positive selection in Tibet fish: Gymnocypris przewalskii.

Myostatin (mstn) is an important member of TGF-ß superfamily, a muscle growth inhibitor. Though mstn has been identified in many organisms, little is known about this gene in highland fish, Gymnocypris przewalskii endemic to the Qinghai-Tibetan Plateau. In this study, we first cloned two paralogous mstn genes (mstn1 and mstn2) from G. przewalskii through homologue cloning. The 3D structures of both Mstn proteins varied in the numbers of ß-sheets and conformations of α-helices. The branch-site model showed that mstn1 has undergone positive selection, and two positively selected sites (107M and 181T) were located on the random coils of the 3D protein structure. Expression patterns indicated that the mstn1 expressed widely, while the mstn2 only expressed in the muscle and brain. During the early stage of embryo development, the expression levels of both mstn paralogous genes showed different increasing trends. These results suggest that it is diverging in two mstn paralogues of G. przewalskii via specific differences in gene structure, protein structure, selection pressure and gene expression patterns. Taken together, this study provides novel contribution on the research topics of growth related gene function and mechanism of highland fish in extreme aquatic environment on the Qinghai-Tibetan Plateau.

Cloning and characterization of two subunits of calcineurin cDNA in naked carp (Gymnocypris przewalskii) from Lake Qinghai, China.

The naked carp (Gymnocypris przewalskii), a native teleost, plays an important role in maintenance of the ecological balance in the system of Lake Qinghai (altitude, 3.2 km) on the Qinghai-Tibet Plateau in China. Calcineurin (CN) is the only member of the serine/threonine phosphatase family that can be activated by both Ca2+ and calmodulin (CaM) and involved in many important physiological processes such as salt tolerance/adaption. In this report, cDNAs of CN catalytic subunit paralogue isoforms: GpCAα (GenBank accession no.JQ407043), GpCAγ (GenBank accession no. JQ407043), and CN regulatory subunit (GpCB) (GenBank accession no. JQ410473), were isolated from Gymnocypris przewalskii and their expression patterns in embryos developmentwere characterized. Gene expression profile demonstrated that GpCA and GpCB mRNA was distributed ubiquitously in all embryonic stages and showed decline until final stage of development. Immunohistologicalanalysis revealed CN localization in different tissues including kidney, heart, brain, spermary, and gill. Collectively, these results provide molecular basis and clues to further understand the role of CN during embryos development and its function in tissues for the adaptation mechanism of naked carp.