Evidence of high-altitude adaptation in the glyptosternoid fish, Creteuchiloglanis macropterus from the Nujiang River obtained through transcriptome analysis

PacBio Full-Length Transcriptome of a Tetraploid Sinocyclocheilus multipunctatus Provides Insights into the Evolution of Cavefish

Sinocyclocheilus multipunctatus is a second-class nationally protected wild animal in China. As one of the cavefish, S. multipunctatus has strong adaptability to harsh subterranean environments. In this study, we used PacBio SMRT sequencing technology to generate a first representative full-length transcriptome for S. multipunctatus. Sequence clustering analysis obtained 232,126 full-length transcripts. Among all transcripts, 40,487 were annotated in public databases, while 70,300 microsatellites, 2384 transcription factors, and 16,321 long non-coding RNAs were identified. The phylogenetic tree showed that S. multipunctatus shows a closer relationship to Carassius auratus and Cyprinus carpio, phylogenetically diverging from the common ancestor ~14.74 million years ago (Mya). We also found that between 15.6 and 17.5 Mya, S. multipunctatus also experienced an additional whole-genome duplication (WGD) event, which may have promoted the species evolution of S. multipunctatus. Meanwhile, the overall rates of evolutionary of polyploid S. multipunctatus were significantly higher than those of the other cyprinids, and 220 positively selected genes (PSGs) were identified in two sub-genomes of S. multipunctatus. These PSGs are likely to fulfill critical roles in the process of adapting to diverse cave environments. This study has the potential to facilitate future investigations into the genomic characteristics of S. multipunctatus and provide valuable insights into revealing the evolutionary history of polyploid S. multipunctatus.

The relevance of genetic structure in ecotype designation and conservation management

The concept of ecotypes is complex, partly because of its interdisciplinary nature, but the idea is intrinsically valuable for evolutionary biology and applied conservation. The complex nature of ecotypes has spurred some confusion and inconsistencies in the literature, thereby limiting broader theoretical development and practical application. We provide suggestions for how incorporating genetic analyses can ease confusion and help define ecotypes. We approach this by systematically reviewing 112 publications across taxa that simultaneously mention the terms ecotype, conservation and management, to examine the current use of the term in the context of conservation and management. We found that most ecotype studies involve fish, mammals and plants with a focus on habitat use, which at 60% was the most common criterion used for categorization of ecotypes. Only 53% of the studies incorporated genetic analyses, and major discrepancies in available genomic resources among taxa could have contributed to confusion about the role of genetic structure in delineating ecotypes. Our results show that the rapid advances in genetic methods, also for nonmodel organisms, can help clarify the spatiotemporal distribution of adaptive and neutral genetic variation and their relevance to ecotype designations. Genetic analyses can offer empirical support for the ecotype concept and provide a timely measure of evolutionary potential, especially in changing environmental conditions. Genetic variation that is often difficult to detect, including polygenic traits influenced by small contributions from several genes, can be vital for adaptation to rapidly changing environments. Emerging ecotypes may signal speciation in progress, and findings from genome‐enabled organisms can help clarify important selective factors driving ecotype development and persistence, and thereby improve preservation of interspecific genetic diversity. Incorporation of genetic analyses in ecotype studies will help connect evolutionary biology and applied conservation, including that of problematic groups such as natural hybrid organisms and urban or anthropogenic ecotypes.

Genome sequencing and transcriptomic analysis of the Andean killifish Orestias ascotanensis reveals adaptation to high-altitude aquatic life

Orestias ascotanensis (Cyprinodontidae) is a teleost pupfish endemic to springs feeding into the Ascotan saltpan in the Chilean Altiplano (3,700 m.a.s.l.) and represents an opportunity to study adaptations to high-altitude aquatic environments. We have de novo assembled the genome of O. ascotanensis at high coverage. Comparative analysis of the O. ascotanensis genome showed an overall process of contraction, including loss of genes related to G-protein signaling, chemotaxis and signal transduction, while there was expansion of gene families associated with microtubule-based movement and protein ubiquitination. We identified 818 genes under positive selection, many of which are involved in DNA repair. Additionally, we identified novel and conserved microRNAs expressed in O. ascotanensis and its closely-related species, Orestias gloriae. Our analysis suggests that positive selection and expansion of genes that preserve genome stability are a potential adaptive mechanism to cope with the increased solar UV radiation to which high-altitude animals are exposed to.

Genome and population evolution and environmental adaptation of Glyptosternon maculatum on the Qinghai-Tibet Plateau

Persistent uplift means the Qinghai-Tibet Plateau (QTP) is an ideal natural laboratory to investigate genome evolution and adaptation within highland environments. However, how paleogeographic and paleoclimatic events influence the genome and population of endemic fish species remains unclear. Glyptosternon maculatum is an ancient endemic fish found on the QTP and the only critically endangered species in the Sisoridae family. Here, we found that major transposons in the G. maculatum genome showed episodic bursts, consistent with contemporaneous geological and climatic events during the QTP formation. Notably, histone genes showed significant expansion in the G. maculatum genome, which may be mediated by long interspersed nuclear elements (LINE) repetitive element duplications. Population analysis showed that ancestral G. maculatum populations experienced two significant depressions 2.6 million years ago (Mya) and 10 000 years ago, exhibiting excellent synchronization with Quaternary glaciation and the Younger Dryas, respectively. Thus, we propose that paleogeography and paleoclimate were dominating driving forces for population dynamics in endemic fish on the QTP. Tectonic movements and temperature fluctuation likely destroyed the habitat and disrupted the drainage connectivity among populations. These factors may have caused severe bottlenecks and limited migration among ancestral G. maculatum populations, resulting in the low genetic diversity and endangered status of the species today.

Targeting the Mild-Hypoxia Driving Force for Metabolic and Muscle Transcriptional Reprogramming of Gilthead Sea Bream (Sparus aurata) Juveniles

Simple Summary

Reduced oxygen availability generates a number of adaptive features across all the animal kingdom, and the goal of this study was targeting the mild-hypoxia driving force for metabolic and muscle transcriptional reprogramming of gilthead sea bream juveniles. Attention was focused on blood metabolic and muscle transcriptomic landmarks before and after exhaustive exercise. Our results after mild-hypoxia conditioning highlighted an increased contribution of lipid metabolism to whole energy supply to preserve the aerobic energy production, a better swimming performance regardless of changes in feed intake, as well as reduced protein turnover and improved anaerobic fitness with the restoration of normoxia.

Abstract

On-growing juveniles of gilthead sea bream were acclimated for 45 days to mild-hypoxia (M-HYP, 40–60% O₂ saturation), whereas normoxic fish (85–90% O₂ saturation) constituted two different groups, depending on if they were fed to visual satiety (control fish) or pair-fed to M-HYP fish. Following the hypoxia conditioning period, all fish were maintained in normoxia and continued to be fed until visual satiation for 3 weeks. The time course of hypoxia-induced changes was assessed by changes in blood metabolic landmarks and muscle transcriptomics before and after exhaustive exercise in a swim tunnel respirometer. In M-HYP fish, our results highlighted a higher contribution of aerobic metabolism to whole energy supply, shifting towards a higher anaerobic fitness following normoxia restoration. Despite these changes in substrate preference, M-HYP fish shared a persistent improvement in swimming performance with a higher critical speed at exercise exhaustion. The machinery of muscle contraction and protein synthesis and breakdown was also largely altered by mild-hypoxia conditioning, contributing this metabolic re-adjustment to the positive regulation of locomotion and to the catch-up growth response during the normoxia recovery period. Altogether, these results reinforce the presence of large phenotypic plasticity in gilthead sea bream, and highlights mild-hypoxia as a promising prophylactic measure to prepare these fish for predictable stressful events.

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.

Microgeographical adaptation corresponds to elevational distributions of congeneric montane grasshoppers

Local adaptation can occur at small spatial scales relative to the dispersal capacity of species. Alpine ecosystems have sharp environmental clines that offer an opportunity to investigate the effects of fine-scale shifts in species' niche breadth on adaptive genetic processes. Here we examine two grasshopper species endemic to the Australian Alps (Kosciuscola spp.) that differ in elevational niche breadth: one broader, K. usitatus (1400-2200 m), and one narrower, K. tristis (1600-2000 m). We examine signatures of selection with respect to environmental and morphological variables in two mountain regions using F_ST outlier tests and environmental association analyses (EAAs) applied to single nucleotide polymorphism (SNP) data (K. usitatus: 9017 SNPs, n = 130; K. tristis: 7363 SNPs, n = 135). Stronger genetic structure was found in the more narrowly distributed K. tristis, which showed almost twice the number of SNPs under putative selection (10.8%) compared with K. usitatus (5.3%). When examining SNPs in common across species (n = 3058), 260 SNPs (8.5%) were outliers shared across species, and these were mostly associated with elevation, a proxy for temperature, suggesting parallel adaptive processes in response to climatic drivers. Additive polygenic scores (an estimate of the cumulative signal of selection across all candidate loci) were nonlinearly and positively correlated with elevation in both species. However, a steeper correlation in K. tristis indicated a stronger signal of spatially varying selection towards higher elevations. Our study illustrates that the niche breadth of co-occurring and related species distributed along the same environmental cline is associated with differences in patterns of microgeographical adaptation.

Transcriptome analysis of the growth performance of hybrid mandarin fish after food conversion

During recent years, China has become a hotspot for the domestication of mandarin fish, and this is of great commercial value. Although the food preference of domesticated mandarin fish has been studied, little is known about genes regulating their growth. We raised hybrid mandarin fish on artificial feed for 3 months, the results showed that the survival rate of hybrid mandarin fish was 60.00%. Their total length and body weight were 18.34 ±0.43 cm and 100.44 ±4.87 g. The absolute length and weight gain rates were 0.14 cm/d and 1.08 g/d, respectively. Finally, RNA sequencing (RNA-Seq) was performed to identify potential genes and pathways activated in response to growth performance. The transcriptome analysis generated 68, 197 transcripts and 45,871 unigenes. Among them, 1025 genes were up-regulated and 593 genes were down-regulated between the fast- and slow-growth fish. Finally, we obtained 32 differentially expressed genes, which were mainly related to fatty acid biosynthesis (e.g. FASN and ACACB), collecting duct acid secretion (e.g. ATP6E and KCC4), cell cycle (e.g. CDC20 and CCNB), and the insulin-like growth factor (IGF) system (IGFBP1). These pathways might be related to the growth of hybrid mandarin fish. In addition, more potential single nucleotide polymorphisms (SNPs) were detected in the fast-growth fish than in the slow-growth fish. The results suggest that the interaction of metabolism and abundant alleles might determine the growth of hybrid mandarin fish after food conversion.

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.

Comprehensive Transcriptome Analysis Reveals Insights into Phylogeny and Positively Selected Genes of Sillago Species

Sillago species lives in the demersal environments and face multiple stressors, such as localized oxygen depletion, sulfide accumulation, and high turbidity. In this study, we performed transcriptome analyses of seven Sillago species to provide insights into the phylogeny and positively selected genes of this species. After de novo assembly, 82,024, 58,102, 63,807, 85,990, 102,185, 69,748, and 102,903 unigenes were generated from S. japonica, S. aeolus, S. sp.1, S. sihama, S. sp.2, S. parvisquamis, and S. sinica, respectively. Furthermore, 140 shared orthologous exon markers were identified and then applied to reconstruct the phylogenetic relationships of the seven Sillago species. The reconstructed phylogenetic structure was significantly congruent with the prevailing morphological and molecular biological view of Sillago species relationships. In addition, a total of 44 genes were identified to be positively selected, and these genes were potential participants in the stress response, material (carbohydrate, amino acid and lipid) and energy metabolism, growth and differentiation, embryogenesis, visual sense, and other biological processes. We suspected that these genes possibly allowed Sillago species to increase their ecological adaptation to multiple environmental stressors.

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.