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

Population structure and adaptability analysis of Schizothorax o'connori based on whole-genome resequencing

BACKGROUND

Schizothorax o'connori is an endemic fish distributed in the upper and lower reaches of the Yarlung Zangbo River in China. It has experienced a fourth round of whole gene replication events and is a good model for exploring the genetic differentiation and environmental adaptability of fish in the Qinghai-Tibet Plateau. The uplift of the Qinghai-Tibet Plateau has led to changes in the river system, thereby affecting gene exchange and population differentiation between fish populations. With the release of fish whole genome data, whole genome resequencing has been widely used in genetic evolutionary analysis and screening of selected genes in fish, which can better elucidate the genetic basis and molecular environmental adaptation mechanisms of fish. Therefore, our purpose of this study was to understand the population structure and adaptive characteristics of S. o'connori using the whole-genome resequencing method.

RESULTS

The results showed that 23,602,746 SNPs were identified from seven populations, mostly distributed on chromosomes 2 and 23. There was no significant genetic differentiation between the populations, and the genetic diversity was relatively low. However, the Zangga population could be separated from the Bomi, Linzhi, and Milin populations in the cluster analysis. Based on historical dynamics analysis of the population, the size of the ancestral population of S. o'connori was affected by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Glacial Age. The selected sites were mostly enriched in pathways related to DNA repair and energy metabolism.

CONCLUSION

Overall, the whole-genome resequencing analysis provides valuable insights into the population structure and adaptive characteristics of S. o'connori. There was no obvious genetic differentiation at the genome level between the S. o'connori populations upstream and downstream of the Yarlung Zangbo River. The current distribution pattern and genetic diversity are influenced by the late accelerated uplift of the Qinghai Tibet Plateau and the Fourth Ice Age. The selected sites of S. o'connori are enriched in the energy metabolism and DNA repair pathways to adapt to the low temperature and strong ultraviolet radiation environment at high altitude.

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.

Fish Community Structure and Biomass Particle-Size Spectrum in the Upper Reaches of the Jinsha River (China)

To understand the characteristics of the fish community structure and biomass particle-size spectrum in the upper reaches of the Jinsha River, fish and environmental surveys were conducted in 21 segments of the upper reaches of the Jinsha River in September 2019 and June 2020. A total of 4062 fish belonging to 2 orders, 5 families, 18 genera, and 28 species were collected. Among them, Cyprinidae fish were the most abundant (14 species), accounting for 50.00%. The Shannon index and Pielou evenness index values varied from 0.402-1.770 and 0.254-0.680, respectively. The dominant species of fish were Triplophysa stenura, Schizothorax wangchiachii, and Schizopygopsis malacanthus. Redundancy analysis (RDA) was used to analyse the relationship between the fish community and environmental factors. Velocity, altitude, and dissolved oxygen were the main influencing factors of fish community structure differences in the upper reaches of the Jinsha River. The abundance/biomass curves showed that the fish communities in the upper reaches of the Jinsha River were moderately or severely disturbed. The standardized biomass particle-size spectrum of fish showed that the degree of disturbance of fish in tributaries was much lower than that in the main stream. Compared with the historical data, the fish community structure in the Jinsha River has changed significantly, with the number of exotic species increasing, and the individual fish showing miniaturization and younger ages. It is suggested that habitat conservation strategies be adopted in the upper tributaries of the Jinsha River to provide a reference for the restoration of fishery resources and the conservation of fish diversity in the Yangtze River.

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.

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.

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.

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.

Dietary Diversification and Specialization in Neotropical Bats Facilitated by Early Molecular Evolution

Dietary adaptation is a major feature of phenotypic and ecological diversification, yet the genetic basis of dietary shifts is poorly understood. Among mammals, Neotropical leaf-nosed bats (family Phyllostomidae) show unmatched diversity in diet; from a putative insectivorous ancestor, phyllostomids have radiated to specialize on diverse food sources including blood, nectar, and fruit. To assess whether dietary diversification in this group was accompanied by molecular adaptations for changing metabolic demands, we sequenced 89 transcriptomes across 58 species and combined these with published data to compare ∼13,000 protein coding genes across 66 species. We tested for positive selection on focal lineages, including those inferred to have undergone dietary shifts. Unexpectedly, we found a broad signature of positive selection in the ancestral phyllostomid branch, spanning genes implicated in the metabolism of all major macronutrients, yet few positively selected genes at the inferred switch to plantivory. Branches corresponding to blood- and nectar-based diets showed selection in loci underpinning nitrogenous waste excretion and glycolysis, respectively. Intriguingly, patterns of selection in metabolism genes were mirrored by those in loci implicated in craniofacial remodeling, a trait previously linked to phyllostomid dietary specialization. Finally, we show that the null model of the widely-used branch-site test is likely to be misspecified, with the implication that the test is too conservative and probably under-reports true cases of positive selection. Our findings point to a complex picture of adaptive radiation, in which the evolution of new dietary specializations has been facilitated by early adaptations combined with the generation of new genetic variation.

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.

Chromosomal genome of Triplophysa bleekeri provides insights into its evolution and environmental adaptation

BACKGROUND

Intense stresses caused by high-altitude environments may result in noticeable genetic adaptions in native species. Studies of genetic adaptations to high elevations have been largely limited to terrestrial animals. How fish adapt to high-elevation environments is largely unknown. Triplophysa bleekeri, an endemic fish inhabiting high-altitude regions, is an excellent model to investigate the genetic mechanisms of adaptation to the local environment. Here, we assembled a chromosomal genome sequence of T. bleekeri, with a size of ∼628 Mb (contig and scaffold N50 of 3.1 and 22.9 Mb, respectively). We investigated the origin and environmental adaptation of T. bleekeri based on 21,198 protein-coding genes in the genome.

RESULTS

Compared with fish species living at low altitudes, gene families associated with lipid metabolism and immune response were significantly expanded in the T. bleekeri genome. Genes involved in DNA repair exhibit positive selection for T. bleekeri, Triplophysa siluroides, and Triplophysa tibetana, indicating that adaptive convergence in Triplophysa species occurred at the positively selected genes. We also analyzed whole-genome variants among samples from 3 populations. The results showed that populations separated by geological and artificial barriers exhibited obvious differences in genetic structures, indicating that gene flow is restricted between populations.

CONCLUSIONS

These results will help us expand our understanding of environmental adaptation and genetic diversity of T. bleekeri and provide valuable genetic resources for future studies on the evolution and conservation of high-altitude fish species such as T. bleekeri.

Comparative Genomics Identifies Putative Signatures of Sociality in Spiders

Comparative genomics has begun to elucidate the genomic basis of social life in insects, but insight into the genomic basis of spider sociality has lagged behind. To begin, to characterize genomic signatures associated with the evolution of social life in spiders, we performed one of the first spider comparative genomics studies including five solitary species and two social species, representing two independent origins of sociality in the genus Stegodyphus. We found that the two social spider species had a large expansion of gene families associated with transport and metabolic processes and an elevated genome-wide rate of molecular evolution compared with the five solitary spider species. Genes that were rapidly evolving in the two social species relative to the five solitary species were enriched for transport, behavior, and immune functions, whereas genes that were rapidly evolving in the solitary species were enriched for energy metabolism processes. Most rapidly evolving genes in the social species Stegodyphus dumicola were broadly expressed across four tissues and enriched for transport functions, but 12 rapidly evolving genes showed brain-specific expression and were enriched for social behavioral processes. Altogether, our study identifies putative genomic signatures and potential candidate genes associated with spider sociality. These results indicate that future spider comparative genomic studies, including broader sampling and additional independent origins of sociality, can further clarify the genomic causes and consequences of social life.

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.

Comprehensive transcriptome data for endemic Schizothoracinae fish in the Tibetan Plateau

The Schizothoracinae fishes, endemic species in the Tibetan Plateau, are considered as ideal models for highland adaptation and speciation investigation. Despite several transcriptome studies for highland fishes have been reported before, the transcriptome information of Schizothoracinae is still lacking. To obtain comprehensive transcriptome data for Schizothoracinae, the transcriptome of a total of 183 samples from 14 representative Schizothoracinae species, were sequenced and de novo assembled. As a result, about 1,363 Gb transcriptome clean data was obtained. After the assembly, we obtain 76,602-154,860 unigenes for each species with sequence N50 length of 1,564-2,143 bp. More than half of the unigenes were functionally annotated by public databases. The Schizothoracinae fishes in this work exhibited diversified ecological distributions, phenotype characters and feeding habits; therefore, the comprehensive transcriptome data of those species provided valuable information for the environmental adaptation and speciation of Schizothoracinae in the Tibetan Plateau.

Divergent selection and drift shape the genomes of two avian sister species spanning a saline-freshwater ecotone

The role of species divergence due to ecologically based divergent selection-or ecological speciation-in generating and maintaining biodiversity is a central question in evolutionary biology. Comparison of the genomes of phylogenetically related taxa spanning a selective habitat gradient enables discovery of divergent signatures of selection and thereby provides valuable insight into the role of divergent ecological selection in speciation. Tidal marsh ecosystems provide tractable opportunities for studying organisms' adaptations to selective pressures that underlie ecological divergence. Sharp environmental gradients across the saline-freshwater ecotone within tidal marshes present extreme adaptive challenges to terrestrial vertebrates. Here, we sequence 20 whole genomes of two avian sister species endemic to tidal marshes-the saltmarsh sparrow (Ammospiza caudacutus) and Nelson's sparrow (A. nelsoni)-to evaluate the influence of selective and demographic processes in shaping genome-wide patterns of divergence. Genome-wide divergence between these two recently diverged sister species was notably high (genome-wide F ST = 0.32). Against a background of high genome-wide divergence, regions of elevated divergence were widespread throughout the genome, as opposed to focused within islands of differentiation. These patterns may be the result of genetic drift resulting from past tidal march colonization events in conjunction with divergent selection to different environments. We identified several candidate genes that exhibited elevated divergence between saltmarsh and Nelson's sparrows, including genes linked to osmotic regulation, circadian rhythm, and plumage melanism-all putative candidates linked to adaptation to tidal marsh environments. These findings provide new insights into the roles of divergent selection and genetic drift in generating and maintaining biodiversity.

Adaptive Evolution of the Eda Gene and Scales Loss in Schizothoracine Fishes in Response to Uplift of the Tibetan Plateau

Schizothoracine is the predominant wild fish subfamily of the Tibetan plateau (TP). Their scales, pharyngeal teeth and barbels have gradually regressed with increasing altitude. Schizothoracine have been divided into three groups: primitive, specialized and highly specialized. Ectodysplasin-A (Eda) has been considered as a major gene that contributes to the development of skin appendages. The present study cloned the Eda genes of 51 Schizothoracine fish species which represent the three groups and five Barbinae species. Phylogenetic analyses indicated that Eda may have acted as the genetic trigger for scale loss in the Schizothoracine. Furthermore, 14 single nucleotide polymorphisms (SNPs) and two deletions (18 bp and 6 bp in size), were also detected in the Eda coding sequence of the highly specialized group compared to the primitive group. The same SNPs and two indels result in four non-synonymous and two G-X-Y and 1 XY motif indels, which possibly contribute to significant structure changes in the Eda gene. The domain including (G-X-Y)_n motif in the Eda gene is relatively conserved amongst teleosts. Based on the above results, we hypothesize that the evolution of Eda gene might be associated with the scale loss in Schizothoracine fishes in response to the phased uplift of the TP.

Transcriptome Analysis Based on RNA-Seq in Understanding Pathogenic Mechanisms of Diseases and the Immune System of Fish: A Comprehensive Review

In recent years, with the advent of next-generation sequencing along with the development of various bioinformatics tools, RNA sequencing (RNA-Seq)-based transcriptome analysis has become much more affordable in the field of biological research. This technique has even opened up avenues to explore the transcriptome of non-model organisms for which a reference genome is not available. This has made fish health researchers march towards this technology to understand pathogenic processes and immune reactions in fish during the event of infection. Recent studies using this technology have altered and updated the previous understanding of many diseases in fish. RNA-Seq has been employed in the understanding of fish pathogens like bacteria, virus, parasites, and oomycetes. Also, it has been helpful in unraveling the immune mechanisms in fish. Additionally, RNA-Seq technology has made its way for future works, such as genetic linkage mapping, quantitative trait analysis, disease-resistant strain or broodstock selection, and the development of effective vaccines and therapies. Until now, there are no reviews that comprehensively summarize the studies which made use of RNA-Seq to explore the mechanisms of infection of pathogens and the defense strategies of fish hosts. This review aims to summarize the contemporary understanding and findings with regard to infectious pathogens and the immune system of fish that have been achieved through RNA-Seq technology.