Evidence for Adaptation to the Tibetan Plateau Inferred from Tibetan Loach Transcriptomes

Whole-genome relaxed selection and molecular constraints in Triplophysa under adapted Qinghai-Tibetan Plateau

High-altitude environments are inhospitable, but Triplophysa, the largest taxon among the three major fish groups in the Qinghai-Tibetan Plateau (QTP), is an exception. However, the evolutionary profiling of the common ancestor and its contribution to the adaptation of existing QTP native species is unclear.We researched the comparative genomics of Triplophysa species and found that the genome-wide genes of Triplophysa and its ancestry have the characteristics of rapid evolution.Moreover, the rapid evolution of the ancestral genes was caused by relaxed selection. Natural selection analysis showed that more ancestral relaxed selection genes were under strongly purifying selection and showed higher expression in QTP endemic Triplophysa species.The change in natural selection might be associated with the adaptation to QTP. It should be noted that SPT5 homolog, DSIF elongation factor subunit (supt5h) experienced relaxed selection in common ancestral populations of Triplophysa but under purifying selection in extant species, which might be related to hypoxia adaptation of QTP. In summary, the extant species in different environments were used to infer the evolutionary profile of the common ancestor and to identify candidate genes based on changes in natural selection. Our work might provide new clues for understanding adaptation to extreme environments.

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

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

Selection and Gene Duplication Associated With High-Elevation Diversification in Pristimantis, the Largest Terrestrial Vertebrate Genus

The genus Pristimantis diversified in the tropical Andes mountains and is the most speciose genus of terrestrial vertebrates. Pristimantis are notable among frogs in that they thrive at high elevations (>2,000 m) and are direct developers without a tadpole stage. Despite their ecological significance, little is known about the genetic and physiological traits enabling their success. We conducted transcriptomic analysis on seven Pristimantis species sampled across elevations in the Ecuadorean Andes to explore three hypotheses for their success: (i) unique genes are under selection relative to all other frogs, (ii) common selection occurs across all direct developers, or (iii) common selection occurs across all high-elevation frog clades. Comparative analysis with 34 frog species revealed unique positive selection in Pristimantis genes related to aerobic respiration, hemostasis, signaling, cellular transportation of proteins and ions, and immunity. Additionally, we detected positive selection across all direct developers for genes associated with oxygenase activity and metal ion binding. While many genes under selection in Pristimantis were not positively selected in other high-elevation frog species, we identified some shared genes and pathways linked to lipid metabolism, innate immunity, and cellular redox processes. We observed more positive selection in duplicated- versus single-copy genes, while relaxed purifying selection was prevalent in single-copy genes. Notably, copy number of an innate immunity complement gene was positively correlated with Pristimantis species elevation. Our findings contribute novel insights into the genetic basis of adaptation in Pristimantis and provide a foundation for future studies on the evolutionary mechanisms leading to direct development and coping with high elevations.

Rapid altitude displacement induce zebrafish appearing acute high altitude illness symptoms

Rapid ascent to high-altitude areas above 2500 m often leads to acute high altitude illness (AHAI), posing significant health risks. Current models for AHAI research are limited in their ability to accurately simulate the high-altitude environment for drug screening. Addressing this gap, a novel static self-assembled water vacuum transparent chamber was developed to induce AHAI in zebrafish. This study identified 6000 m for 2 h as the optimal condition for AHAI induction in zebrafish. Under these conditions, notable behavioral changes including slow movement, abnormal exploration behavior and static behavior in the Novel tank test. Furthermore, this model demonstrated changes in oxidative stress-related markers included increased levels of malondialdehyde, decreased levels of glutathione, decreased activities of superoxide dismutase and catalase, and increased levels of inflammatory markers IL-6, IL-1β and TNF-α, and inflammatory cell infiltration and mild edema in the gill tissue, mirroring the clinical pathophysiology observed in AHAI patients. This innovative zebrafish model not only offers a more accurate representation of the high-altitude environment but also provides a high-throughput platform for AHAI drug discovery and pathogenesis research.

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.

Integrated metabolomic and transcriptomic responses to heat stress in a high-altitude fish, Triplophysa siluroides

Species in Triplophysa display strong adaptability to the extreme environment of the plateau, thus offering an ideal model to study the molecular mechanism of fish adaptation to environmental stress. In the present study, we conducted integrated analysis of the transcriptome and metabolism of liver tissue in Triplophysa siluroides under heat stress (28 °C) and control (10 °C) conditions to identify heat stress-induced genes, metabolites and pathways. RNA-Seq identified 2373 differentially expressed genes, which consisted of 1360 upregulated genes and 1013 downregulated genes, in the heat stress group vs. the control group. Genes in the heat shock protein (Hsp) family, including Hsp40, Hsp70, Hsp90 and other Hsps, were strongly upregulated by heat stress. Pathway enrichment analysis revealed that the PI3K/AKT/mTOR and protein processing in the endoplasmic reticulum (ER) pathways were significantly affected by heat stress. Metabolism sequencing identified a total of 155 differentially abundant metabolites, including 118 significantly upregulated metabolites and 37 downregulated metabolites. Combined analysis of the transcriptome and metabolism results showed that ubiquitin-dependent proteolysis and purine metabolism pathways were enhanced in response to acute heat stress to protect cells from damage under stress conditions. The results of this study may contribute to our understanding of the underlying molecular mechanism of the heat stress response in cold-water fish.

High-Altitude Andean H194R HIF2A Allele Is a Hypomorphic Allele

For over 10,000 years, Andeans have resided at high altitude where the partial pressure of oxygen challenges human survival. Recent studies have provided evidence for positive selection acting in Andeans on the HIF2A (also known as EPAS1) locus, which encodes for a central transcription factor of the hypoxia-inducible factor pathway. However, the precise mechanism by which this allele might lead to altitude-adaptive phenotypes, if any, is unknown. By analyzing whole genome sequencing data from 46 high-coverage Peruvian Andean genomes, we confirm evidence for positive selection acting on HIF2A and a unique pattern of variation surrounding the Andean-specific single nucleotide variant (SNV), rs570553380, which encodes for an H194R amino acid substitution in HIF-2α. Genotyping the Andean-associated SNV rs570553380 in a group of 299 Peruvian Andeans from Cerro de Pasco, Peru (4,338 m), reveals a positive association with increased fraction of exhaled nitric oxide, a marker of nitric oxide biosynthesis. In vitro assays show that the H194R mutation impairs binding of HIF-2α to its heterodimeric partner, aryl hydrocarbon receptor nuclear translocator. A knockin mouse model bearing the H194R mutation in the Hif2a gene displays decreased levels of hypoxia-induced pulmonary Endothelin-1 transcripts and protection against hypoxia-induced pulmonary hypertension. We conclude the Andean H194R HIF2A allele is a hypomorphic (partial loss of function) allele.

Unrecognized diversity, genetic structuring, and phylogeography of the genus Triplophysa (Cypriniformes: Nemacheilidae) sheds light on two opposite colonization routes during Quaternary glaciation that occurred in the Qilian Mountains

In recent years, the influence of historical geological and climatic events on the evolution of flora and fauna in the Tibetan Plateau has been a hot research topic. The Qilian Mountain region is one of the most important sources of biodiversity on the Qinghai-Tibet Plateau. Many species existed in the region during the Pleistocene glacial oscillation, and the complex geographical environment provided suitable conditions for the survival of local species. The shrinkage, expansion, and transfer of the distribution range and population size of species have significant effects on genetic diversity and intraspecific differentiation. To reveal the effects of geological uplift and climate oscillation on the evolution of fish populations in the Qilian Mountains, we investigated the genetic structure, phylogenetic relationship, and phylogeographical characteristics of genus Triplophysa species in the Qilian Mountains using the mitochondrial DNA gene (COI), three nuclear genes (RAG1, sRH, and Myh6) and 11 pairs of nuclear microsatellite markers. We collected 11 species of genus Triplophysa living in the Qilian Mountains, among which Triplophysa hsutschouensis and Triplophysa papillosolabiata are widely distributed in the rivers on the northern slope of the Qilian Mountains. There was a high degree of lineage differentiation among species, and the genetic diversity of endemic species was low. The different geographical groups of T. papillosolabiata presented some allogeneic adaptation and differentiation, which was closely related to the changes in the river system. Except for the population expansion event of T. hsutschouensis during the last glacial period of the Qinghai-Tibet Plateau (0.025 MYA), the population sizes of other plateau loach species remained stable without significant population expansion. Starting from the east and west sides of the Qilian Mountains, T. hsutschouensis, and T. papillosolabiata showed two species colonization routes in opposite directions. The geological events of the uplift of the Qinghai-Tibet Plateau and the climatic oscillation of the Quaternary glaciation had a great influence on the genetic structure of the plateau loach in the Qilian Mountains, which promoted the genetic differentiation of the plateau loach and formed some unique new species. The results of this study have important guiding significance for fish habitat protection in the Qilian Mountains.

Complete genome sequence of the Pogostemon cablin bacterial wilt pathogen Ralstonia solanacearum strain SY1

BACKGROUND

Ralstonia solanacearum causes bacterial wilt of Pogostemon cablin which is an important aromatic herb and also the main materials of COVID-19 therapeutic traditional drugs. However, we are lacking the information on the genomic sequences of R. solanacearum isolated from P. cablin.

OBJECTIVE

The acquisition and analysis of this whole-genome sequence of the P. cablin bacterial wilt pathogen.

METHODS

An R. solanacearum strain, named SY1, was isolated from infected P. cablin plants, and the complete genome sequence was sequenced and analyzed.

RESULTS

The SY1 strain contains a 3.70-Mb chromosome and a 2.18-Mb megaplasmid, with GC contents of 67.57% and 67.41%, respectively. A total of 3308 predicted genes were located on the chromosome and 1657 genes were located in the megaplasmid. SY1 strain has 273 unique genes compared with five representative R. solanacearum strains, and these genes were enriched in the plant-pathogen interaction pathway. SY1 possessed a higher syntenic relationship with phylotype I strains, and the arsenal of type III effectors predicted in SY1 were also more closely related to those of phylotype I strains. SY1 contained 14 and 5 genomic islands in its chromosome and megaplasmid, respectively, and two prophage sequences in its chromosome. In addition, 215 and 130 genes were annotated as carbohydrate-active enzymes and antibiotic resistance genes, respectively.

CONCLUSION

This is the first genome-scale assembly and annotation for R. solanacearum which isolated from infected P. cablin plants. The arsenal of virulence and antibiotic resistance may as the determinants in SY1 for infection of P. cablin plants.

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

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

Contrasting population differentiation in two sympatric Triplophysa loaches on the Qinghai-Tibet Plateau

Genetic differentiation in aquatic organisms is usually shaped by drainage connectivity. Sympatric aquatic species are thus expected to show similar population differentiation patterns and similar genetic responses to their habitats. Water bodies on the Qinghai-Tibet Plateau (QTP) have recently experienced dramatic physicochemical changes, threatening the biodiversity of aquatic organisms on the "roof of the world." To uncover ecological genetics in Tibetan loaches (Triplophysa)-the largest component of the QTP ichthyofauna-we characterized population differentiation patterns and adaptive mechanisms to salinity change in two sympatric and phylogenetically closely related Tibetan loaches, T. stewarti and T. stenura, by integrating population genomic, transcriptomic, and electron probe microanalysis approaches. Based on millions of genome-wide SNPs, the two Tibetan loach species show contrasting population differentiation patterns, with highly geographically structured and clear genetic differentiation among T. stewarti populations, whereas there is no such observation in T. stenura, which is also supported by otolith microchemistry mapping. While limited genetic signals of parallel adaption to salinity changes between the two species are found from either genetic or gene expression variation perspective, a catalog of genes involved in ion transport, energy metabolism, structural reorganization, immune response, detoxification, and signal transduction is identified to be related to adaptation to salinity change in Triplophysa loaches. Together, our findings broaden our understanding of the population characteristics and adaptive mechanisms in sympatric Tibetan loach species and would contribute to biodiversity conservation and management of aquatic organisms on the QTP.

Comparative transcriptome and adaptive evolution analysis on the main liver and attaching liver of Pareuchiloglanis macrotrema

Pareuchiloglanis macrotrema is a glyptosternoid fish belonging to the Siluriform family and is endemic to the Qinghai-Tibet Plateau tributaries. P. macrotrema is an ideal model for studying the adaptive evolution of fish at high altitudes. P. macrotrema has two attaching livers connected to the main liver, a common feature in most Sisoridae fishes but is a special phenomenon relative to other vertebrates. Using RNA-Seq, 42 differentially expressed genes were found between the main liver and attaching liver, of which 31 were upregulated and 11 were downregulated in the main liver. The major differentially expressed genes between the main liver and attaching liver of P. macrotrema are related to metabolism, immunity, and digestive processes. Meanwhile, a comparative transcriptome analysis was carried out on P. macrotrema fish and six non-plateau Siluriformes fishes. We found 268 positively selected genes in P. macrotrema that are related to energy metabolism, immunity, and hypoxic responses. The findings of this study highlight the gene expression differences between the main liver and attaching livers of Sisoridae fishes and provide greater insight into the evolution of Tibetan fishes.

Complete genome sequence of the kiwifruit bacterial canker pathogen Pseudomonas savastanoi strain MHT1

BACKGROUND

Pseudomonas savastanoi is an important plant pathogen that infects and causes symptoms in a variety of economically important crops, causing considerable loss of yield and quality. Because there has been no research reported to date on bacterial canker of kiwifruit (Actinidia chinensis) plants caused by P. savastanoi and, in particular, no in-depth studies of the complete genome sequence or pathogenic mechanism, long-lasting and environmentally friendly control measures against this pathogen in kiwifruit are lacking. This study therefore has both theoretical value and practical significance.

RESULTS

We report the complete genome sequence of P. savastanoi strain MHT1, which was first reported as the pathogen causing bacterial canker in kiwifruit plants. The genome consists of a 6.00-Mb chromosome with 58.5% GC content and 5008 predicted genes. Comparative genome analysis of four sequenced genomes of representative P. savastanoi strains revealed that 230 genes are unique to the MHT1 strain and that these genes are enriched in antibiotic metabolic processes and metabolic pathways, which may be associated with the drug resistance and host range observed in this strain. MHT1 showed high syntenic relationships with different P. savastanoi strains. Furthermore, MHT1 has eight conserved effectors that are highly homologous to effectors from P. syringae, Pseudomonas amygdali, and Ralstonia solanacearum strains. The MHT1 genome contains six genomic islands and two prophage sequences. In addition, 380 genes were annotated as antibiotic resistance genes and another 734 as encoding carbohydrate-active enzymes.

CONCLUSION

The whole-genome sequence of this kiwifruit bacterial canker pathogen extends our knowledge of the P. savastanoi genome, sets the stage for further studies of the interaction between kiwifruit and P. savastanoi, and provides an important theoretical foundation for the prevention and control of bacterial canker.

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

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

Genomic 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.

Adaptation of Fig Wasps (Agaodinae) to Their Host Revealed by Large-Scale Transcriptomic Data

Simple Summary

Research on fig wasps has made a considerable contribution to the understanding of insect–plant interactions. However, the molecular mechanisms underlying fig wasp host specificity are poorly understood. This study reports on a relatively large-scale transcriptomic dataset of 25 fig wasp species. We outline potential genetic mechanisms underlying the specific host adaptation by investigating changes in a gene family, in evolutionary rates, and in genes under positive selection. The transcriptome datasets reported here (1) provide new insights into the evolutionary diversification and host specificity of fig wasps and (2) contribute to a growing dataset on fig wasp genomics.

Abstract

Figs and fig wasps are highly species-specific and comprise a model system for studying co-evolution and co-speciation. The evolutionary relationships and molecular adaptations of fig wasps to their fig hosts are poorly understood, and this is in part due to limited sequence data. Here, we present large-scale transcriptomic datasets of 25 fig wasp species with the aim of uncovering the genetic basis for host specificity. Our phylogenetic results support the monophyly of all genera associated with dioecious figs, and two genera associated with monoecious figs, Eupristina and Platyscapa, were revealed to be close relatives. We identified gene loss and gain, potentially rapidly evolving genes, and genes under positive selection. Potentially functional changes were documented and we hypothesize as to how these may determine host specificity. Overall, our study provides new insights into the evolutionary diversification of fig wasps and contributes to our understanding of adaptation in this group.

Recent progress in research on the gut microbiota and highland adaptation on the Qinghai-Tibet Plateau

Microbial communities that inhabit the host's intestine influence many aspects of the host's health and bear the adaptive potential to alterations in harsh environments and diets. The Qinghai-Tibet Plateau represents one of the harshest environments in the world. Preliminary progress has been made in identifying the communities of gut microbes in Indigenous Tibetans and non-human animals. However, due to the complexity of microbial communities, the effects of gut microbes on the host's health and high-plateau adaptation remain unexplained. Herein, we review the latest progress in identifying factors affecting the gut microbiota of native Tibetans and non-human animals and highlight the complex interactions between the gut microbiota, health and highland adaptation, which provides a basis for exploring the correlations between the gut microbiota and clinical indexes in native highland residents and travellers, as well as developing microbiota-based strategies to mitigate health risks for tourists and treatments for mountain sickness during high-altitude travel in the future.

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.

Genomic Evidence for Convergent Molecular Adaptation in Electric Fishes

Fishes have independently evolved electric organs (EOs) at least six times, and the electric fields are used for communication, defense, and predation. However, the genetic basis of convergent evolution of EOs remains unclear. In this study, we conducted comparative genomic analyses to detect genes showing signatures of positive selection and convergent substitutions in electric fishes from three independent lineages (Mormyroidea, Siluriformes, and Gymnotiformes). Analysis of 4,657 orthologs between electric fishes and their corresponding control groups identified consistent evidence for accelerated evolution in electric fish lineages. A total of 702 positively selected genes (PSGs) were identified in electric fishes, and many of these genes corresponded to cell membrane structure, ion channels, and transmembrane transporter activity. Comparative genomic analyses revealed that widespread convergent amino acid substitutions occurred along the electric fish lineages. The overlap of convergent genes and PSGs was identified as adaptive convergence, and a subset of genes was putatively associated with electrical and muscular activities, especially scn4aa (a voltage-gated sodium channel gene). Our results provide hints to the genetic basis for the independent evolution of EOs during millions of years of evolution.

Comparative genomics analysis of two banana Fusarium wilt biocontrol endophytes Bacillus subtilis R31 and TR21 provides insights into their differences on phytobeneficial trait

Fusarium wilt of banana is considered one of the most destructive plant diseases. Bacillus subtilis R31 and TR21, isolated from Dendrobium sp. leaves, exhibit different phytobeneficial effects on banana Fusarium wilt bio-controlling. Here, we performed genome sequencing and comparative genomics analysis of R31 and TR21 to enhance our understanding of the different phytobeneficial traits. These results revealed that the strain-specific genes of R31 involved in sporulation, quorum sensing, and antibiotic synthesis allow R31 to present a better capacity of sporulation, rhizosphere adaptation, and quorum sensing than TR21. Selective pressure analysis indicated that the glycosylase and endo-alpha-(1- > 5)-L-arabinanase genes were strong positive selected, which may contribute to the TR21 to colonize well in banana's vascular bundles. Altogether, our findings presented here should advance further agricultural application of R31 and TR21 as two promising resources of plant growth promotion and biological control via genetic engineering.

Genomics and lipidomics analysis of the biotechnologically important oleaginous red yeast Rhodotorula glutinis ZHK provides new insights into its lipid and carotenoid metabolism

Background

Rhodotorula glutinis is recognized as a biotechnologically important oleaginous red yeast, which synthesizes numerous meritorious compounds with wide industrial usages. One of the most notable properties of R. glutinis is the formation of intracellular lipid droplets full of carotenoids. However, the basic genomic features that underlie the biosynthesis of these valuable compounds in R. glutinis have not been fully documented. To reveal the biotechnological potential of R. glutinis, the genomics and lipidomics analysis was performed through the Next-Generation Sequencing and HPLC-MS-based metabolomics technologies.

Results

Here, we firstly assemble the genome of R. glutinis ZHK into 21.8 Mb, containing 30 scaffolds and 6774 predicted genes with a N50 length of 14, 66,672 bp and GC content of 67.8%. Genome completeness assessment (BUSCO alignment: 95.3%) indicated the genome assembly with a high-quality features. According to the functional annotation of the genome, we predicted several key genes involved in lipids and carotenoids metabolism as well as certain industrial enzymes biosynthesis. Comparative genomics results suggested that most of orthologous genes have underwent the strong purifying selection within the five Rhodotorula species, especially genes responsible for carotenoids biosynthesis. Furthermore, a total of 982 lipids were identified using the lipidomics approaches, mainly including triacylglycerols, diacylglyceryltrimethylhomo-ser and phosphatidylethanolamine.

Conclusion

Using whole genome shotgun sequencing, we comprehensively analyzed the genome of R. glutinis and predicted several key genes involved in lipids and carotenoids metabolism. By performing comparative genomic analysis, we show that most of the ortholog genes have undergone strong purifying selection within the five Rhodotorula species. Furthermore, we identified 982 lipid species using lipidomic approaches. These results provided valuable resources to further advance biotechnological applications of R .glutinis.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07244-z.

Miocene Diversification and High-Altitude Adaptation of Parnassius Butterflies (Lepidoptera: Papilionidae) in Qinghai-Tibet Plateau Revealed by Large-Scale Transcriptomic Data

Simple Summary

Parnassius butterflies have contributed to fundamental studies in biogeography, insect–plant interactions, and other fields of conservation biology and ecology. However, the early evolutionary pattern and molecular adaptation mechanism of this alpine butterfly group to high altitudes in Qinghai–Tibet Plateau are poorly understood up to now. In this study, we report for the first time, a relatively large-scale transcriptomic dataset of eight Parnassius species and their two closely related papilionid species, a dated phylogeny based on hundreds of gene sequences, and potential genetic mechanisms underlying the high-altitude adaptation by investigating changes in evolutionary rates and positively selected genes. Overall, our findings indicate that the transcriptome data sets reported here can provide some new insights into the spatiotemporally evolutionary pattern and high altitude adaptation of Parnassius butterflies from the extrinsic and intrinsic view, and will support further expressional and functional studies that will help interested researchers to address evolution, biodiversity and conservation questions concerning Parnassius and other butterfly species.

Abstract

The early evolutionary pattern and molecular adaptation mechanism of alpine Parnassius butterflies to high altitudes in Qinghai–Tibet Plateau are poorly understood up to now, due to difficulties in sampling, limited sequence data, and time calibration issues. Here, we present large-scale transcriptomic datasets of eight representative Parnassius species to reveal the phylogenetic timescale and potential genetic basis for high-altitude adaptation with multiple analytic strategies using 476 orthologous genes. Our phylogenetic results strongly supported that the subgenus Parnassius formed a well-resolved basal clade, and the subgenera Tadumia and Kailasius were closely related in the phylogenetic trees. In addition, molecular dating analyses showed that the Parnassius began to diverge at about 13.0 to 14.3 million years ago (middle Miocene), correlated with their hostplant’s spatiotemporal distributions, as well as geological and palaeoenvironmental changes of the Qinghai–Tibet Plateau. Moreover, the accelerated evolutionary rate, candidate positively selected genes and their potentially functional changes were detected, probably contributed to the high-altitude adaptation of Parnassius species. Overall, our study provided some new insights into the spatiotemporally evolutionary pattern and high altitude adaptation of Parnassius butterflies from the extrinsic and intrinsic view, which will help to address evolution, biodiversity, and conservation questions concerning Parnassius and other butterfly species.

Comparative transcriptomic analysis of fireflies (Coleoptera: Lampyridae) to explore the molecular adaptations to fresh water

Aquatic insects are well adapted to freshwater environments, but the molecular basis of these adaptations remains largely unknown. Most firefly species (Coleoptera: Lampyridae) are terrestrial, but the larvae of several species are aquatic. Here, larval and adult transcriptomes from Aquatica leii (freshwater) and Lychnuris praetexta (terrestrial) were generated to test whether the genes associated with metabolic efficiency and morphology have undergone adaptive evolution to fresh water. The aquatic fireflies had a significantly lower ratio of nonsynonymous to synonymous substitutions than the terrestrial insects, indicating a genomewide evolutionary constraint in the aquatic fireflies. We identified 341 fast-evolving genes and 116 positively selected genes in the aquatic fireflies. Of these, 76 genes exhibiting both fast evolution and positive selection were primarily involved in ATP production, energy metabolism and the hypoxia response. We identified 7,271 differentially expressed genes (DEGs) in A. leii (adults versus larvae) and 8,309 DEGs in L. praetexta (adults versus larvae). DEGs specific to the aquatic firefly (n = 1,445) were screened via interspecific comparisons (A. leii versus L. praetexta) and were significantly enriched for genes involved in metabolic efficiency (e.g., ATP production, hypoxia, and immune responses) and certain aspects of morphology (e.g., cuticle chitin, tracheal and compound eye morphology). These results indicate that sequence and expression-level changes in genes associated with both metabolic efficiency and morphological attributes related to the freshwater lifestyle contributed to freshwater adaptation in fireflies. This study provides new insights into the molecular mechanisms of aquatic adaptation in insects.

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.

Whole genome sequencing and comparative genomic analysis of oleaginous red yeast Sporobolomyces pararoseus NGR identifies candidate genes for biotechnological potential and ballistospores-shooting

BACKGROUND

Sporobolomyces pararoseus is regarded as an oleaginous red yeast, which synthesizes numerous valuable compounds with wide industrial usages. This species hold biotechnological interests in biodiesel, food and cosmetics industries. Moreover, the ballistospores-shooting promotes the colonizing of S. pararoseus in most terrestrial and marine ecosystems. However, very little is known about the basic genomic features of S. pararoseus. To assess the biotechnological potential and ballistospores-shooting mechanism of S. pararoseus on genome-scale, the whole genome sequencing was performed by next-generation sequencing technology.

RESULTS

Here, we used Illumina Hiseq platform to firstly assemble S. pararoseus genome into 20.9 Mb containing 54 scaffolds and 5963 predicted genes with a N50 length of 2,038,020 bp and GC content of 47.59%. Genome completeness (BUSCO alignment: 95.4%) and RNA-seq analysis (expressed genes: 98.68%) indicated the high-quality features of the current genome. Through the annotation information of the genome, we screened many key genes involved in carotenoids, lipids, carbohydrate metabolism and signal transduction pathways. A phylogenetic assessment suggested that the evolutionary trajectory of the order Sporidiobolales species was evolved from genus Sporobolomyces to Rhodotorula through the mediator Rhodosporidiobolus. Compared to the lacking ballistospores Rhodotorula toruloides and Saccharomyces cerevisiae, we found genes enriched for spore germination and sugar metabolism. These genes might be responsible for the ballistospores-shooting in S. pararoseus NGR.

CONCLUSION

These results greatly advance our understanding of S. pararoseus NGR in biotechnological potential and ballistospores-shooting, which help further research of genetic manipulation, metabolic engineering as well as its evolutionary direction.

A Chromosome-Scale Reference Assembly of a Tibetan Loach, Triplophysa siluroides

Cobitoidea is one of the two superfamilies in Cypriniformes; however, few genomes have been sequenced for Cobitoidea fishes. Here, we obtained a total of 252.90 Gb of short Illumina reads and 31.60 Gb of long PacBio Sequel reads, representing approximate genome coverage of 256× and 50×, respectively. The final assembled genome is about 583.47 Mb with contig N50 sizes of 2.87 Mb, which accounts for 91.44% of the estimated genome size of 638.07 Mb. Using Hi-C-based chromatin contact maps, 99.31% of the genome assembly was placed into 25 chromosomes, and the N50 is 22.3 Mb. The gene annotation completeness was evaluated by BUSCO, and 2,470 of the 2,586 conserved genes (95.5%) could be found in our assembly. Repetitive elements were calculated to reach 33.08% of the whole genome. Moreover, we identified 25,406 protein-coding genes, of which 92.59% have been functionally annotated. This genome assembly will be a valuable genomic resource to understand the biology of the Tibetan loaches and will also set a stage for comparative analysis of the classification, diversification, and adaptation of fishes in Cobitoidea.

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 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.

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

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

The Antarctic sea ice alga Chlamydomonas sp. ICE-L provides insights into adaptive patterns of chloroplast evolution

BACKGROUND

The ice alga Chlamydomonas sp. ICE-L is the main contributor to primary productivity in Antarctic sea ice ecosystems and is well adapted to the extremely harsh environment. However, the adaptive mechanism of Chlamydomonas sp. ICE-L to sea-ice environment remains unclear. To study the adaptive strategies in Chlamydomonas sp. ICE-L, we investigated the molecular evolution of chloroplast photosynthetic genes that are essential for the accumulation of carbohydrate and energy living in Antarctic sea ice.

RESULTS

The 60 chloroplast protein-coding genes of Chlamydomonas sp. ICE-L were obtained, and the branch-site test detected significant signatures of positive selection on atpB, psaB, and rbcL genes in Chlamydomonas sp. ICE-L associated with the photosynthetic machinery. These positively selected genes were further identified as being under convergent evolution between Chlamydomonas sp. ICE-L and the halotolerant alga Dunaliella salina.

CONCLUSIONS

Our study provides evidence that the phototrophic component of Chlamydomonas sp. ICE-L exhibits adaptive evolution under extreme environment. The positive Darwinian selection operates on the chloroplast protein-coding genes of Antarctic ice algae adapted to extreme environment following functional-specific and lineages-specific patterns. In addition, three positively selected genes with convergent substitutions in Chlamydomonas sp. ICE-L were identified, and the adaptive modifications in these genes were in functionally important regions of the proteins. Our study provides a foundation for future experiments on the biochemical and physiological impacts of photosynthetic genes in green algae.

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

Background

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

Results

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

Conclusions

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

Electronic supplementary material

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

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.

Selection of reference genes for qRT-PCR and expression analysis of high-altitude-related genes in grassland caterpillars (Lepidoptera: Erebidae: Gynaephora) along an altitude gradient

Changes in gene expression patterns can reflect the adaptation of organisms to divergent environments. Quantitative real‐time PCR (qRT‐PCR) is an important tool for ecological adaptation studies at the gene expression level. The quality of the results of qRT‐PCR analysis largely depends on the availability of reliable reference genes (RGs). To date, reliable RGs have not been determined for adaptive evolution studies in insects using a standard approach. Here, we evaluated the reliability of 17 candidate RGs for five Gynaephora populations inhabiting various altitudes of the Tibetan Plateau (TP) using four independent (geNorm, NormFinder, BestKeeper, and the deltaCt method) and one comprehensive (RefFinder) algorithms. Our results showed that EF1‐α, RPS15, and RPS13 were the top three most suitable RGs, and a combination of these three RGs was the most optimal for normalization. Conversely, RPS2,ACT, and RPL27 were the most unstable RGs. The expression profiles of two target genes (HSP70 and HSP90) were used to confirm the reliability of the chosen RGs. Additionally, the expression patterns of four other genes (GPI,HIF1A,HSP20, and USP) associated with adaptation to extreme environments were assessed to explore the adaptive mechanisms of TP Gynaephora species to divergent environments. Each of these six target genes showed discrepant expression patterns among the five populations, suggesting that the observed expression differences may be associated with the local adaptation of Gynaephora to divergent altitudinal environments. This study is a useful resource for studying the adaptive evolution of TP Gynaephora to divergent environments using qRT‐PCR, and it also acts as a guide for selecting suitable RGs for ecological and evolutionary studies in insects.

Gene expression variations in high-altitude adaptation: a case study of the Asiatic toad (Bufo gargarizans)

Background

Genome-wide investigation of molecular mechanisms for high-altitude adaptation has attracted great attention in the last few years. In order to understand the contribution of gene expression level variations to high-altitude adaptation in Asiatic toads (Bufo gargarizans), we implemented a reciprocal transplant experiment between low- and high-altitude sites and sequenced 12 transcriptomes from brain, heart, and liver tissues.

Results

A large number of genes with expression differences (DEGs) between high- and low-altitude individuals (193 fixed and 844 plastic) were identified, and the majority of them were tissue specific. Heart displayed the largest number of DEGs, both plastic and fixed. Fixed DEGs were particularly concentrated in functions associated with muscle contraction, and the majority of them were down-regulated in high-altitude individuals. Plastic DEGs were highly concentrated in several energy metabolism related functional categories, and the majority of them were also down-regulated at high-altitude environments. In liver samples, genes associated with nutrient metabolism experienced a broad-scale expression down-regulation in high-altitude toads.

Conclusions

These broadly suppressed expression patterns at high altitudes are in strong contrast to those of endothermic homeotherms, suggesting poikilothermic vertebrates may have adopted different strategies at high altitudes. Our results strongly support that both genotypic specialization and phenotypic plasticity play crucial role in adaptation to high altitude for Asiatic toads. Poikilothermic vertebrates are among the most hypoxia-tolerant animals known, and many molecular mechanisms remain elusive. We hope that our results will provide useful directions for future research.

Electronic supplementary material

The online version of this article (doi:10.1186/s12863-017-0529-z) contains supplementary material, which is available to authorized users.

Global analysis of transcriptome sequences highlights accelerated evolution of immune genes in Danio choprae and Danio albolineatus

Danio fishes, a small type animal with short sexual cycles, are model vertebrate species. To investigate the genic evolution of this genus, the transcriptomes from Danio choprae and Danio albolineatus were sequenced by Illumina HiSeq 4000 platform. A total of 128,427,304 sequence reads from two Danio fishes were generated by Next Generation Sequencing. The resulting in two assemblies contained 88,682 and 88,029 unigenes in the Danio choprae and Danio albolineatus. Analysis of the orthologs from the Danio choprae and Danio albolineatus provided consistent evidence for the accelerated genic evolution in the Danio fishes. Several genes referring to immune functions under positive selection were identified by branch site model analysis, such as REL, GTF2E1, STAT6, MPG in Danio choprae and CYP17A1, ADORA2A, MYCN in Danio albolineatus. Our data provide novel insights into the adaptation in Danio fishes and is useful for understanding the genetic basis of adaptation in zebrafish.

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 transcriptomic analysis of Tibetan Schizothoracinae fish Gymnocypris przewalskii reveals how it adapts to a high altitude aquatic life

Background

Understanding the genetic basis of adaptation to high altitude life is of paramount importance for preserving and managing genetic diversity in highland animals. This objective has been addressed mainly in terrestrial fauna but rarely in aquatic animals. Tibetan Schizothoracinae fish is the ideal model system in evolutionary biology, carrying key insights into evolutionary genetics of speciation and adaptation at high altitude. Gymnocypris przewalskii is the newly formed Schizothoracinae fish species in the Tibetan Plateau, inhabits chronic cold, extreme saline and alkaline aquatic environment in Lake Qinghai, thus evolving the unique genomic signatures to adapt extremely severe environments.

Results

To characterize its genomic features, we assembled de novo transcriptome of G. przewalskii from Lake Qinghai. Intriguingly, by comparative genomic analyses of G. przewalskii and 8 other fish species, we identified potential expansions in gene families related to energy metabolism, transport and developmental functions, possibly underlying the adaptation to these environmental stresses. Through comprehensive molecular evolution analyses, we found that sets of genes controlling mitochondrion, ion homoeostasis, acid-base balance and innate immunity show significant signals of positive selection. Compared to previous studies on highland fishes, we failed to identify any positively selected genes related to hypoxia response.

Conclusions

Our findings provide comprehensive insights into the genetic basis of teleost fish that underlie their adaptation to extreme high altitude aquatic life on the Tibetan Plateau.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-017-0925-z) contains supplementary material, which is available to authorized users.

Characterization and functional analysis of hypoxia-inducible factor HIF1α and its inhibitor HIF1αn in tilapia

Hypoxia is a major cause of fish morbidity and mortality in the aquatic environment. Hypoxia-inducible factors are very important modulators in the transcriptional response to hypoxic stress. In this study, we characterized and conducted functional analysis of hypoxia-inducible factor HIF1α and its inhibitor HIF1αn in Nile tilapia (Oreochromis niloticus). By cloning and Sanger sequencing, we obtained the full length cDNA sequences for HIF1α (2686bp) and HIF1αn (1308bp), respectively. The CDS of HIF1α includes 15 exons encoding 768 amino acid residues and the CDS of HIF1αn contains 8 exons encoding 354 amino acid residues. The complete CDS sequences of HIF1α and HIF1αn cloned from tilapia shared very high homology with known genes from other fishes. HIF1α show differentiated expression in different tissues (brain, heart, gill, spleen, liver) and at different hypoxia exposure times (6h, 12h, 24h). HIF1αn expression level under hypoxia is generally increased (6h, 12h, 24h) and shows extremely highly upregulation in brain tissue under hypoxia. A functional determination site analysis in the protein sequences between fish and land animals identified 21 amino acid sites in HIF1α and 2 sites in HIF1αn as significantly associated sites (α = 0.05). Phylogenetic tree-based positive selection analysis suggested 22 sites in HIF1α as positively selected sites with a p-value of at least 95% for fish lineages compared to the land animals. Our study could be important for clarifying the mechanism of fish adaptation to aquatic hypoxia environment.

Regulatory Architecture of Gene Expression Variation in the Threespine Stickleback Gasterosteus aculeatus

Much adaptive evolutionary change is underlain by mutational variation in regions of the genome that regulate gene expression rather than in the coding regions of the genes themselves. An understanding of the role of gene expression variation in facilitating local adaptation will be aided by an understanding of underlying regulatory networks. Here, we characterize the genetic architecture of gene expression variation in the threespine stickleback (Gasterosteus aculeatus), an important model in the study of adaptive evolution. We collected transcriptomic and genomic data from 60 half-sib families using an expression microarray and genotyping-by-sequencing, and located expression quantitative trait loci (eQTL) underlying the variation in gene expression in liver tissue using an interval mapping approach. We identified eQTL for several thousand expression traits. Expression was influenced by polymorphism in both cis- and trans-regulatory regions. Trans-eQTL clustered into hotspots. We did not identify master transcriptional regulators in hotspot locations: rather, the presence of hotspots may be driven by complex interactions between multiple transcription factors. One observed hotspot colocated with a QTL recently found to underlie salinity tolerance in the threespine stickleback. However, most other observed hotspots did not colocate with regions of the genome known to be involved in adaptive divergence between marine and freshwater habitats.

Identification of Triplophysa species from the Qinghai-Tibetan Plateau (QTP) and its adjacent regions through DNA barcodes

The genus Triplophysa is the largest and most difficult to identity morphologically fish group of superfamily Cobitoidea with 140 currently valid species, and is mainly distributed in the Qinghai-Tibetan Plateau (QTP) and adjacent regions. Most species within this genus possess highly similar morphological characteristics for adaption to the highland environment and are very difficult to be identified only based on morphology. The traditional species identification, mainly based on external morphological diagnostic characters, leads to inconsistent results in many cases. Herein, we provided a molecular method based on mitochondrial cytochrome c subunit I (COI) for the identification of Triplophysa fishes. Thirty-three Triplophysa species, 244 individuals, were used to determine whether barcoding was effective in discriminating species for this genus. The mean intraspecific and interspecific K2P distances ranged from 0 to 14.9% (mean, 2.9%) and 0 to 23.4% (mean, 9.7%), respectively. The tree-based analysis displayed most of species formed discrete clusters with strong bootstrap support values (>90%). The results showed that most of Triplophysa species could be identified by DNA barcode and indicated DNA barcode could be used as a molecular marker for these species.

Hypoxic stress -responsive genes in air breathing catfish, Clarias magur (Hamilton 1822) and their possible physiological adaptive function

The Indian catfish, Clarias magur (previous name C. batrachus) is an air breathing fish, inhabitant of aquatic bodies characterized by low dissolved oxygen levels. It is exposed to hypoxic conditions in its natural habitat. Thus, it can be useful model to study the mechanism of hypoxia stress tolerance. In C. magur, molecular processes facilitating its adaptation to hypoxia stress remain largely unexplored, in part due to unavailability of genomic resources. The suppression subtractive hybridization technique (SSH) was employed to compare the differential expression of transcripts under experimental hypoxic conditions, to that of normoxic conditions. Twelve subtracted cDNA libraries (six each forward and reverse) were constructed from brain, heart, liver, muscle, spleen and head kidney tissues. A total of 2020 clones were screened and sequenced, resulting into 1805 high quality expressed sequence tags (ESTs). Annotation of these differentially expressed ESTs resulted into the identification of genes involved in vast majority of pathways/processes affecting metabolism, cellular processes, signal transduction and/or immune functions. Additionally, 18 potential novel genes expressed in hypoxia stress exposed fish were also identified. The study had catalogued the differentially expressed genes from hypoxia stress induced C. magur, where most of them are reported for the first time in a hypoxia-tolerant fish species. The results not only provided insights for the hypoxia stress altered cellular functions in C. magur, but also generated a valuable functional genomics resource to assist targeted studies on functional genomics and future genome projects.

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.

Transcriptome analyses provide insights into the phylogeny and adaptive evolution of the mangrove fern genus Acrostichum

The mangrove fern genus Acrostichum grows in the extremely unstable marine intertidal zone under harsh conditions, such as high salt concentrations, tidal rhythms and long-term climate changes. To explore the phylogenetic relationships and molecular mechanisms underlying adaptations in this genus, we sequenced the transcriptomes of two species of Acrostichum, A. aureum and A. speciosum, as well as a species in the sister genus, Ceratopteris thalictroides. We obtained 47,517, 36,420 and 60,823 unigenes for the three ferns, of which 24.39-45.63% were annotated using public databases. The estimated divergence time revealed that Acrostichum adapted to the coastal region during the late Cretaceous, whereas the two mangrove ferns from the Indo West-Pacific (IWP) area diverged more recently. Two methods (the modified branch-site model and the K_h method) were used to identify several positively selected genes, which may contribute to differential adaptation of the two Acrostichum species to different light and salt conditions. Our study provides abundant transcriptome data and new insights into the evolution and adaptations of mangrove ferns in the inhospitable intertidal zone.

Genetic signals of high-altitude adaptation in amphibians: a comparative transcriptome analysis

BACKGROUND

High-altitude adaptation provides an excellent system for studying how organisms cope with multiple environmental stressors and interacting genetic modifications. To explore the genetic basis of high-altitude adaptation in poikilothermic animals, we acquired transcriptome sequences from a high-altitude population and a low-altitude population of the Asiatic toad (Bufo gargarizans). Transcriptome data from another high-altitude amphibian, Rana kukunoris and its low-altitude relative R. chensiensis, which are from a previous study, were also incorporated into our comparative analysis.

RESULTS

More than 40,000 transcripts were obtained from each transcriptome, and 5107 one-to-one orthologs were identified among the four taxa for comparative analysis. A total of 29 (Bufo) and 33 (Rana) putative positively selected genes were identified for the two high-altitude species, which were mainly concentrated in nutrient metabolism related functions. Using SNP-tagging and F_ST outlier analysis, we further tested 89 other nutrient metabolism related genes for signatures of natural selection, and found that two genes, CAPN2 and ITPR1, were likely under balancing selection. We did not detect any positively selected genes associated with response to hypoxia.

CONCLUSIONS

Amphibians clearly employ different genetic mechanisms for high-altitude adaptation compared to endotherms. Modifications of genes associated with nutrient metabolism feature prominently while genes related to hypoxia tolerance appear to be insignificant. Poikilotherms represent the majority of animal diversity, and we hope that our results will provide useful directions for future studies of amphibians as well as other poikilotherms.

Mitogenomic perspectives on the origin of Tibetan loaches and their adaptation to high altitude

Tibetan loaches are the largest group of Tibetan fishes and are well adapted to the Tibetan Plateau. To investigate the origin of Tibetan loaches and their adaptations to the Tibetan Plateau, we determined 32 complete mitochondrial genomes that included 29 Tibetan loach species, two Barbatula species and Schistura longus. By combining these newly determined sequences with other previously published mitochondrial genomes, we assembled a large mitogenomic data set (11,433 bp) of 96 species in the superfamily Cobitoidea, to investigate the phylogenetic status of the genus Triplophysa. The resulting phylogeny strongly supported that the genus Triplophysa forms a monophyletic group within Nemacheilidae. Our molecular dating time suggests that the lineage leading to the Tibetan loaches and other loaches diverged approximately 23.5 Ma, which falls within the period of recent major uplifts of the Tibetan Plateau in the Early Miocene. Selection analyses revealed that the mitochondrial protein-coding genes of Tibetan loaches have larger ratios of nonsynonymous to synonymous substitutions than do those of non-Tibetan loaches, indicating that Tibetan loaches accumulated more nonsynonymous mutations than non-Tibetan loaches and exhibited rapid evolution. Two positively selected sites were identified in the ATP8 and ND1 genes.