Comprehensive Transcriptome Analysis of Six Catfish Species from an Altitude Gradient Reveals Adaptive Evolution in Tibetan Fishes

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.

Current Status of Omics in Biological Quality Elements for Freshwater Biomonitoring

Freshwater ecosystems have been experiencing various forms of threats, mainly since the last century. The severity of this adverse scenario presents unprecedented challenges to human health, water supply, agriculture, forestry, ecological systems, and biodiversity, among other areas. Despite the progress made in various biomonitoring techniques tailored to specific countries and biotic communities, significant constraints exist, particularly in assessing and quantifying biodiversity and its interplay with detrimental factors. Incorporating modern techniques into biomonitoring methodologies presents a challenging topic with multiple perspectives and assertions. This review aims to present a comprehensive overview of the contemporary advancements in freshwater biomonitoring, specifically by utilizing omics methodologies such as genomics, metagenomics, transcriptomics, proteomics, metabolomics, and multi-omics. The present study aims to elucidate the rationale behind the imperative need for modernization in this field. This will be achieved by presenting case studies, examining the diverse range of organisms that have been studied, and evaluating the potential benefits and drawbacks associated with the utilization of these methodologies. The utilization of advanced high-throughput bioinformatics techniques represents a sophisticated approach that necessitates a significant departure from the conventional practices of contemporary freshwater biomonitoring. The significant contributions of omics techniques in the context of biological quality elements (BQEs) and their interpretations in ecological problems are crucial for biomonitoring programs. Such contributions are primarily attributed to the previously overlooked identification of interactions between different levels of biological organization and their responses, isolated and combined, to specific critical conditions.

Molecular Mechanisms Underlying Vertebrate Adaptive Evolution: A Systematic Review

Adaptive evolution is a process in which variation that confers an evolutionary advantage in a specific environmental context arises and is propagated through a population. When investigating this process, researchers have mainly focused on describing advantageous phenotypes or putative advantageous genotypes. A recent increase in molecular data accessibility and technological advances has allowed researchers to go beyond description and to make inferences about the mechanisms underlying adaptive evolution. In this systematic review, we discuss articles from 2016 to 2022 that investigated or reviewed the molecular mechanisms underlying adaptive evolution in vertebrates in response to environmental variation. Regulatory elements within the genome and regulatory proteins involved in either gene expression or cellular pathways have been shown to play key roles in adaptive evolution in response to most of the discussed environmental factors. Gene losses were suggested to be associated with an adaptive response in some contexts. Future adaptive evolution research could benefit from more investigations focused on noncoding regions of the genome, gene regulation mechanisms, and gene losses potentially yielding advantageous phenotypes. Investigating how novel advantageous genotypes are conserved could also contribute to our knowledge of adaptive evolution.

Cytochemical and comparative transcriptome analyses elucidate the formation and ecological adaptation of three types of pollen coat in Zingiberaceae

BACKGROUND

The pollen ornate surface of flowering plants has long fascinated and puzzled evolutionary biologists for their variety. Each pollen grain is contained within a pollen wall consisting of intine and exine, over which the lipoid pollen coat lies. The cytology and molecular biology of the development of the intine and exine components of the pollen wall are relatively well characterised. However, little is known about the pollen coat, which confers species specificity. We demonstrate three types of pollen coat in Zingiberaceae, a mucilage-like pollen coat and a gum-like pollen coat, along with a pollen coat more typical of angiosperms. The morphological differences between the three types of pollen coat and the related molecular mechanisms of their formation were studied using an integrative approach of cytology, RNA-seq and positive selection analysis.

RESULTS

Contrary to the 'typical' pollen coat, in ginger species with a mucilage-like (Caulokaempferia coenobialis, Cco) or gum-like (Hornstedtia hainanensis, Hhn) pollen coat, anther locular fluid was still present at the bicellular pollen (BCP) stage of development. Nevertheless, there were marked differences between these species: there were much lower levels of anther locular fluid in Hhn at the BCP stage and it contained less polysaccharide, but more lipid, than the locular fluid of Cco. The set of specific highly-expressed (SHE) genes in Cco was enriched in the 'polysaccharide metabolic process' annotation term, while 'fatty acid degradation' and 'metabolism of terpenoids and polyketides' were significantly enriched in SHE-Hhn.

CONCLUSIONS

Our cytological and comparative transcriptome analysis showed that different types of pollen coat depend on the residual amount and composition of anther locular fluid at the BCP stage. The genes involved in 'polysaccharide metabolism' and 'transport' in the development of a mucilage-like pollen coat and in 'lipid metabolism' and 'transport' in the development of a gum-like pollen coat probably evolved under positive selection in both cases. We suggest that the shift from a typical pollen coat to a gum-like or mucilage-like pollen coat in flowering plants is an adaptation to habitats with high humidity and scarcity of pollinators.

Transcriptomic and metabolomic analyses reveal the altitude adaptability and evolution of different-colored flowers in alpine Rhododendron species

Understanding the molecular mechanisms and evolutionary process of plant adaptation to the heterogeneous environment caused by altitude gradients in plateau mountain ecosystems can provide novel insight into species' responses to global changes. Flower color is the most conspicuous and highly diverse trait in nature. Herein, the gene expression patterns, evolutionary adaptation and metabolites changes of different-colored flowers of alpine Rhododendron L. species along altitude gradients were investigated based on a combined analysis of transcriptomics and metabolomics. Differentially expressed genes were found to be related to the biosynthesis of carbohydrates, fatty acids, amino acids and flavonoids, suggesting their important roles in the altitude adaptability of Rhododendron species. The evolution rate of high-altitude species was faster than that of low-altitude species. Genes related to DNA repair, mitogen-activated protein kinase and ABA signal transduction, and lipoic acid and propanoate metabolism were positively selected in the flowers of high-altitude Rhododendron species and those associated with carotenoid biosynthesis pathway, ABA signal transduction and ethylene signal transduction were positively selected in low-altitude species. These results indicated that the genes with differentiated expressions or functions exhibit varying evolution during the adaptive divergence of heterogeneous environment caused by altitude gradients. Flower-color variation might be attributed to the significant differences in gene expression or metabolites related to sucrose, flavonoids and carotenoids at the transcription or metabolism levels of Rhododendron species. This work suggests that Rhododendron species have multiple molecular mechanisms in their adaptation to changing environments caused by altitude gradients.

Genomic and functional evidence reveals convergent evolution in fishes on the Tibetan Plateau

High-altitude environments are strong drivers of adaptive evolution in endemic organisms. However, little is known about the genetic mechanisms of convergent adaptation among different lineages, especially in fishes. There are three independent fish groups on the Tibetan Plateau: Tibetan Loaches, Schizothoracine fishes and Glyptosternoid fishes; all are well adapted to the harsh environmental conditions. They represent an excellent example of convergent evolution but with an unclear genetic basis. We used comparative genomic analyses between Tibetan fishes and fishes from low altitudes and detected genomic signatures of convergent evolution in fishes on the Tibetan Plateau. The Tibetan fishes exhibited genome-wide accelerated evolution in comparison with a control set of fishes from low altitudes. A total of 368 positively selected genes were identified in Tibetan fishes, which were enriched in functional categories related to energy metabolism and hypoxia response. Widespread parallel amino acid substitutions were detected among the Tibetan fishes and a subset of these substitutions occurred in positively selected genes associated with high-altitude adaptation. Functional assays suggested that von Hippel-Lindau (VHL) tumour suppressor genes from Tibetan fishes enhance hypoxia-inducible factor (HIF) activity convergently under hypoxia compared to low-altitude fishes. The results provide genomic and functional evidence supporting convergent genetic mechanisms for high-altitude adaptation in fishes on the Tibetan Plateau.

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.

Transcriptomic evidences of local thermal adaptation for the native fish Colossoma macropomum (Cuvier, 1818)

Brazil has five climatically distinct regions, with an annual average temperature difference up to 14 ºC between the northern and southern extremes. Environmental variation of this magnitude can lead to new genetic patterns among farmed fish populations. Genetically differentiated populations of tambaqui (Colossoma macropomum Cuvier, 1818), an important freshwater fish for Brazilian continental aquaculture, may be associated with regional adaptation. In this study, we selected tambaquis raised in two thermally distinct regions, belonging to different latitudes, to test this hypothesis. De novo transcriptome analysis was performed to compare the significant differences of genes expressed in the liver of juvenile tambaqui from a northern population (Balbina) and a southeastern population (Brumado). In total, 2,410 genes were differentially expressed (1,196 in Balbina and 1,214 in Brumado). Many of the genes are involved in a multitude of biological functions such as biosynthetic processes, homeostasis, biorhythm, immunity, cell signaling, ribosome biogenesis, modification of proteins, intracellular transport, structure/cytoskeleton, and catalytic activity. Enrichment analysis based on biological networks showed a different protein interaction profile for each population, whose encoding genes may play potential functions in local thermal adaptation of fish to their respective farming environments.

Analysis of Natural Selection of Immune Genes in Spinibarbus caldwelli by Transcriptome Sequencing

Spinibarbus caldwelli is an omnivorous cyprinid fish that is distributed widely in China. To investigate the adaptive evolution of S. caldwelli, the muscle transcriptome was sequenced by Illumina HiSeq 4000 platform. A total of 80,447,367 reads were generated by next-generation sequencing. Also, 211,386 unigenes were obtained by de novo assembly. Additionally, we calculated that the divergence time between S. caldwelli and Sinocyclocheilus grahami is 23.14 million years ago (Mya). And both of them diverged from Ctenopharyngodon idellus 46.95 Mya. Furthermore, 38 positive genes were identified by calculating Ka/Ks ratios from 9225 orthologs. Among them, several immune-related genes were identified as positively selected, such as POLR3B, PIK3C3, TOPORS, FASTKD3, CYPLP1A1, and UACA. Our results throw light on the nature of the natural selection of S. caldwelli and contribute to future immunological and transcriptome studies.

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.

De novo transcriptome assembly and sex-biased gene expression in the gonads of Amur catfish (Silurus asotus)

Amur catfish is extensively distributed and cultured in Asian countries. Despite of economic importance, the genomic information of this species remains limited. A reference transcriptome of Amur catfish was assembled and the sex-biased gene expression in the gonads was characterized using RNA-sequencing. The assembled transcriptome of Amur catfish consisted of 74,840 transcripts. The N50, mean length and max length of transcripts are 1970, 1235 and 16,748 bp. Putative sex-specific transcripts were identified and sex-specific expression of the representative genes was verified by RT-PCR. Differential expression analysis identified 5401 ovary-biased and 5618 testis-biased genes. The ovary-biased genes were mainly enriched in pathways such as RNA transport and ribosome biogenesis in eukaryotes. The testis-biased genes were enriched in calcium signaling and cytokine-cytokine receptor interaction, etc. Our data provide a valuable genomic resource for further investigating the genetic basis of sex determination, sex differentiation and sexual dimorphism of catfish.

Transcriptomic Responses of Two Ecologically Divergent Populations of Japanese Mantis Shrimp (Oratosquilla oratoria) under Thermal Stress

Simple Summary

Rising ocean temperature would change the seawater chemistry and affect the external and internal physiology of crustaceans due to their lack of certain efficient temperature regulators. In addition, the infraspecific populations of crustaceans might also have different response strategies to the rising of temperature. Therefore, we identified the transcriptomic variations to the same thermal stress between ecologically divergent populations of Oratosquilla oratoria. The aim of this study was to investigate the population-specific function genes and relevant pathways in response to thermal stress in O. oratoria. The results showed that gene-expressed variation was in a population-specific pattern, which indicated that the local environment could lead to the evolvement of changes in gene regulation, ultimately leading to adaptive divergences. Additionally, we found several genes with large pleiotropic effects in the Zhoushan population, which might indicate that the regulation mechanisms of the Zhoushan population were more efficient than those of the Qingdao population under same thermal stress. The results provided some novel insights into the local adaptive differences of the infraspecific populations of O. oratoria and other crustaceans.

Abstract

Crustaceans are generally considered more sensitive to ocean warming due to their lack of certain efficient regulators. However, the alterations in the physiology and behavior of crustaceans in response to thermal stress differ vastly even among the infraspecific populations of heterogeneous landscapes. Consequently, understanding the impact of temperature fluctuation on crustacean infraspecific populations might be essential for maintaining a sustainable persistence of populations at existing locations. In the present study, we chose the Japanese mantis shrimp (Oratosquilla oratoria) as the representative crustacean population, and conducted transcriptome analyses in two divergent O. oratoria populations (the Zhoushan and Qingdao populations) under same thermal stress (20–28 °C) to identify the population-specific expression response to thermal stress. The results showed significant differences in gene expressions, GO terms and metabolic pathways between the two populations. We hypothesized that intraspecific mutations in the same or different genes might lead to thermal adaptive divergences. Temperature increases from 20–28 °C produced significant enrichment in GO terms and altered the metabolic pathways in the Zhoushan population despite the lack of differentially expressed unigenes. Therefore, several functional genes with large pleiotropic effects may underlie the response to thermal stress in the Zhoushan population. Furthermore, the most significantly enriched biological processes of the Qingdao population were associated with the state or activity of cells and its significant enriched pathways with genetic information processing as well as immune and environmental information processing. In contrast, the differentially regulated unigenes of the Zhoushan population were primarily involved in the regulatory cellular and transcription processes and the most significant pathways found were metabolic and digestive. Consequently, the regulatory mechanisms of the Zhoushan population are probably more efficient than those of the Qingdao population under the same thermal stress.

Comparative transcriptomic analysis reveals genetic divergence and domestication genes in Diospyros

BACKGROUND

Persimmon (Diospyros kaki) is the most economically cultivated species belonging to the genus Diospyros. However, little is known about the interspecific diversity and mechanism of domestication, partly due to the lack of genomic information that is available for closely related species of D. kaki (DK). Here, we performed transcriptome sequencing on nine samples, including DK, a variety of DK and seven closely related species, to evaluate the interspecific genetic divergence and to identify candidate genes involved in persimmon domestication.

RESULTS

We obtained a total of 483,421 unigenes with N50 at 1490 bp in the nine Diospyros samples and identified 2603 orthogroups that were shared among all the samples using OrthoMCL analysis. A phylogenetic tree was established based on the tandem 2603 one-to-one single copy gene alignments, showing that DK was closely related to D. kaki var. silvestris (DKV) and that it clustered with the clade of D. deyangnsis (DD) and was farthest from the D. cathayensis (DC) species. The nonsynonymous substitutions (Ka), via synonymous substitution (Ks) ratios, was directly proportional to the genetic relationship of the different species. The higher the Ka/Ks ratios, the longer the distance was. Moreover, 31 positively selected genes (PSGs) involved in carbohydrate metabolism and phenolic metabolism were identified and isolated, and nearly all PSGs except the MATE gene had a high expression in the DK or DKV species. It was hypothesized that these genes might contribute to the domestication of the DK species. Finally, we developed the expressed sequence tag-simple sequence repeat (EST-SSR) and identified 2 unique amplicons DKSSR10 and DKSSR39: the former was absent in the DC species but was present in the other species, the latter had a long amplification product in the DJ species.

CONCLUSION

This study presents the first transcriptome resources for the closely related species of persimmon and reveals interspecific genetic divergence. It is speculated that DK is derived from the hybridization of DD and DO species. Furthermore, our analysis suggests candidate PSGs that may be crucial for the adaptation, domestication, and speciation of persimmon relatives and suggests that DKSSR10 and DKSSSR39 could potentially serve as species-specific molecular markers.

Comparative transcriptomics of 3 high-altitude passerine birds and their low-altitude relatives

High-altitude environments present strong stresses for living organisms, which have driven striking phenotypic and genetic adaptations. While previous studies have revealed multiple genetic adaptations in high-altitude species, how evolutionary history (i.e., phylogenetic background) contributes to similarity in genetic adaptations to high-altitude environments is largely unknown, in particular in a group of birds. We explored this in 3 high-altitude passerine birds from the Qinghai-Tibet Plateau and their low-altitude relatives in lowland eastern China. We generated transcriptomic data for 5 tissues across these species and compared sequence changes and expression shifts between high- and low-altitude pairs. Sequence comparison revealed that similarity in all 3 high-altitude species was high for genes under positive selection (218 genes) but low in amino acid substitutions (only 4 genes sharing identical amino acid substitutions). Expression profiles for all genes identified a tissue-specific expression pattern (i.e., all species clustered by tissue). By contrast, an altitude-related pattern was observed in genes differentially expressed between all 3 species pairs and genes associated with altitude, suggesting that the high-altitude environment may drive similar expression shifts in the 3 high-altitude species. Gene expression level, gene connectivity, and the interactions of these 2 factors with altitude were correlated with evolutionary rates. Our results provide evidence for how gene sequence changes and expression shifts work in a concerted way in a group of high-altitude birds, leading to similar evolution routes in response to high-altitude environmental stresses.

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.

Comparative transcriptomics of sympatric species of coral reef fishes (genus: Haemulon)

Background

Coral reefs are major hotspots of diversity for marine fishes, yet there is still ongoing debate on the mechanisms that promote divergence in these rich ecosystems. Our understanding of how diversity originates in this environment could be enhanced by investigating the evolutionary dynamics of closely related fishes with overlapping ranges. Here, we focus on grunts of the genus Haemulon, a group of coral reef fishes with 15 species in the Western Atlantic, 11 of which are syntopic.

Methods

Wild fish samples from three sympatric species of the Caribbean: Haemulon flavolineatum, H. carbonarium and H. macrostomum, were collected while SCUBA diving. RNA was extracted from livers, and the transcriptomes were assembled and annotated to investigate positive selection (Pairwise d _N/d _S) and patterns of gene expression between the three species.

Results

Pairwise d _N/d _S analyses showed evidence of positive selection for genes associated with immune response, cranial morphology and formation of the anterior-posterior axis. Analyses of gene expression revealed that despite their sympatric distribution, H. macrostomum showed upregulation of oxidation-reduction machinery, while there was evidence for activation of immune response in H. carbonarium.

Discussion

Overall, our analyses suggest closely related grunts show important differences in genes associated with body shape and feeding morphology, a result in-line with previous morphological studies in the group. Further, despite their overlapping distribution they interact with their environment in distinct fashions. This is the largest compendium of genomic information for grunts thus far, representing a valuable resource for future studies in this unique group of coral reef fishes.

Comparative transcriptomic analysis of Tibetan Gynaephora to explore the genetic basis of insect adaptation to divergent altitude environments

Adaptation of insects to different altitudes remain largely unknown, especially those endemic to the Tibetan Plateau (TP). Here, we generated the transcriptomes of Gynaephora menyuanensis and G. alpherakii, inhabiting different high altitudes on the TP, and used these and the previously available transcriptomic and genomic sequences from low-altitude insects to explore potential genetic basis for divergent high-altitude adaptation in Gynaephora. An analysis of 5,869 orthologous genes among Gynaephora and other three low-altitude insects uncovered that fast-evolving genes and positively selected genes (PSGs) in the two Gynaephora species were enriched in energy metabolism and hypoxia response categories (e.g. mitochondrion, oxidation-reduction process, and response to oxidative stress). Particularly, mTOR signaling pathway involving hypoxia was enriched by PSGs, indicating this well-known pathway in mammal hypoxia adaptation may be an important signaling system in Gynaephora. Furthermore, some PSGs were associated with response to hypoxia (e.g. cytochrome proteins), cold (e.g. dehydrogenase) and DNA repair (e.g. DNA repair proteins). Interestingly, several insect-specific genes that were associated with exoskeleton and cuticle development (e.g. chitinase and ecdysteroids) had experienced positive selection, suggesting the specific adaptive mechanisms in insects. This study is favourable for understanding the adaptive evolution of Gynaephora and even TP insects to divergent altitudes.

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

BACKGROUND

Organisms living at high altitudes face low oxygen and temperature conditions; thus, the genetic mechanisms underlying the adaptations in these organisms merit investigation. The glyptosternoid fish, Creteuchiloglanis macropterus mainly inhabits regions with gradual increases in altitudes along the Nujiang River and might serve as an appropriate evolutionary model for detecting adaptation processes in environments with altitude changes.

RESULTS

We constructed eleven RNA-sequencing (RNA-seq) libraries of C. macropterus collected from five locations at different altitudes to identify the genetic signatures of high-altitude adaptation. The comparative genomic analysis indicated that C. macropterus has an accelerated evolutionary rate compared with that of fishes in the lowland, and fishes at higher altitudes might evolve faster. Functional enrichment analysis of the fast-evolving and positively selected genes, differentially expressed genes and highly expressed genes, showed that these genes were involved in many functions related to energy metabolism and hypoxia.

CONCLUSIONS

Our study provides evidence of high-altitude adaptation in C. macropterus, and the detected adaptive genes might be a resource for future investigations of adaptations to high-altitude environments in other fishes.

Population genetics analysis of the Nujiang catfish Creteuchiloglanis macropterus through a genome-wide single nucleotide polymorphisms resource generated by RAD-seq

Advances in genome scanning using high-throughput sequencing technologies has led to a revolution in studies of non-model organisms. The glyptosternoid fish Creteuchiloglanis macropterus, is widely distributed in the main stem and tributaries of the Nujiang River basin. Here, we analyzed IIB restriction-site-associated DNA (2b-RAD) sequences and mitochondrial DNA sequences, to assess the genomic signature of adaptation by detecting and estimating the degree of genetic differentiation among ten Creteuchiloglanis macropterus populations from the Nujiang River. The analyses revealed significant population differentiation among the up-tributaries, main stem, mid-tributary and low-tributary. Annotation of contigs containing outlier SNPs revealed that the candidate genes showed significant enrichment in several important biological process terms between up-tributaries and low-tributary, and exhibited prominent enrichment in the term macromolecular metabolic process between all tributaries and the main stem. Population dynamics analyses indicated that the Late Pleistocene glaciations strongly influenced the demographic history of C. macropterus. Our results provide strong evidence for the utility of RAD-seq in population genetics studies, and our generated SNP resource should provide a valuable tool for population genomics studies of C. macropterus in the future.

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.

Comprehensive Transcriptome Analysis Provides Evidence of Local Thermal Adaptation in Three Loaches (Genus: Misgurnus)

The geographic distribution of three Misgurnus species, M. anguillicaudatus, M. bipartitus, and M. mohoity, displays a specific pattern in China, coincident with temperature zones. In this study, we sequenced the transcriptomes of these three species and used the sequences to investigate the lineage-specific adaptations within the genus Misgurnus. In total, 51 orphan genes (19 in M. anguillicaudatus, 18 in M. bipartitus, and 14 in M. mohoity) that may contribute to the species-specific adaptations were identified. An analysis of 1392 one-to-one orthologous genes revealed significantly higher ratios of nonsynonymous-to-synonymous substitutions in the M. mohoity lineage than in M. anguillicaudatus. The genes displaying signatures of positive selection and rapid evolution in Misgurnus were involved in four function categories, (1) energy metabolism; (2) signal transduction; (3) membrane; and (4) cell proliferation or apoptosis, implying that these candidate genes play critical roles in the thermal adaptation of the fish to their living environments. We also detected more than five positively selected sites in cldn15lb and isca1, which function as important factors in paracellular Na⁺ transport and Fe/S cluster assembly, respectively. Overall, our study provides valuable insights into the adaptive evolution of loaches from different temperature zones in China and is a foundation for future studies to clarify the genetic basis of temperature adaptation in fishes.