Analysis of the erythropoietin of a Tibetan Plateau schizothoracine fish (Gymnocypris dobula) reveals enhanced cytoprotection function in hypoxic environments

Transcriptomic analysis of immune-related genes in Pacific white snook (Centropomus viridis) gills infected with the monogenean parasite Rhabdosynochus viridisi

The parasite Rhabdosynochus viridisi (Platyhelminthes: Monogenea) infects the Pacific white snook Centropomus viridis gills and can cause adverse effects in the aquaculture industry. The immune responses of Pacific white snook to monogenean infections are poorly understood. Thus, this study aimed to identify differentially expressed genes (DEGs) in the gills of Pacific white snook juveniles experimentally infected with R. viridisi, emphasizing immune-related genes and pathways activated or suppressed during the infection. RNA sequencing was performed on the gills of uninfected (control) and infected fish. The algorithm Seq2Fun was selected without a reference transcriptome to map the reads to transcripts of fishes available from a database for gene orthologs (EcoOmics) and obtain the counting table. The ExpressAnalyst software was used for differential expression and functional analyses. A total of 20,106 transcripts were found, and 1430 (7 %) were differentially expressed genes (DEGs) between infected and control groups. We identified 860 (60 %) downregulated and 570 (40 %) upregulated genes. Thirteen canonical pathways after the Kyoto Encyclopedia of Genes and Genomes (KEGG) database were overrepresented, and most of the DEGs were downregulated, suggesting the inactivation of these pathways. The functions of most of the DEGs with higher fold change found in this study are poorly understood in fish. Even though the well-known pro-inflammatory cytokines remained unchanged in infected gills of C. viridis, and transforming growth factor β (tgfβ) was downregulated, interleukin-17 ligands il17d and il17a/f1, as well as C-X-C motif chemokine receptor 2 (cxcr2) genes were upregulated, indicating that the infection with R. viridisi promotes Th17-like immunity. Overexpression of plasma B cell activity markers such as immunoglobulin light chain-like genes and the v-set pre-B cell surrogate light chain 3 (vpreb3) was also detected in this study. The possible implications of DEGs related to calcium imbalance, hypoxia adaptation, hemostasis, and immunity are discussed. These results will support future studies to improve the prevention and treatment of monogenean infections in finfish aquaculture.

The effect of tris (1,3-dichloro-2-propyl) phosphate on the early embryonic heart development of Oryzias melastigma

The flame retardant tri (1, 3-dichloro-2-propyl) phosphate (TDCIPP) is widely present in environmental media and organisms. People have paid much attention to the growth and developmental toxicity of TDCIPP, but there is little information about its cardiotoxicity and potential mechanisms. In this study, marine medaka (Oryzias melastigma) embryos were exposed to TDCIPP solutions (0, 0.05, 0.5, 5, and 50 μg/L) for 21 days to investigate the adverse effects of TDCIPP on cardiac development. The results showed that TDCIPP exposure altered the heart rate at different stages of embryonic development. In addition, 50 μg/L TDCIPP resulted in increased sinus venosus (SV)-bulbus arteriosus (BA) distance, pericardial cysts, and cardiac linearization in newly hatched fish. During embryonic development, the expression level of key genes regulating cardiac development is disturbed. The early stage of cardiac development is the sensitive window period for the toxic effects of TDCIPP. Oxidative stress was observed in newly hatched juveniles, but no significant lipid peroxidation damage was observed. In addition, vitellogenin (VTG) levels in juvenile fish were significantly reduced. Our results show that TDCIPP exposure induces cardiotoxicity in marine medaka embryos, which is induced in the early stages and promotes heart defects by amplifying inflammatory responses at a later stage.

Comparative mitochondrial genomics analysis of selected species of Schizothoracinae sub family to explore the differences at mitochondrial DNA level

A comprehensive analysis of the whole mitochondrial genomes of the Schizothoracinae subfamily of the family Cyprinidae has been revealed for the first time. The species analyzed include Schizothorax niger, Schizothorax esocinus, Schizothorax labiatus and Schizothorax plagoistomus. The total mitochondrial DNA (mtDNA) length was determined to be 16585 bp, 16583 bp, 16582 bp and 16576 bp, respectively with 13 protein-coding genes, 2 rRNA genes, 22 tRNA genes and 2 non-coding area genes. The combined mean base compositions of the four species were as follows: A: 29.91 % T: 25.47 % G: 17.65 % C 27.01 %. The range of the GC content is 45-44 %, respectively. All protein coding genes (PCGs) commenced with the typical ATG codon, except for the cytochrome c oxidase subunit 1 (COX1) gene with GTG. The analysis of vital amino acid biosynthesis genes (COX1, ATPase 6, ATPase 8) in four different species revealed no significant differences. All 13 PCGs had Ka/Ks ratios that were all lesser than one, demonstrating purifying selection on those molecules. These tRNA genes were predicted to fold into the typical cloverleaf secondary structures with normal base pairing and ranged in size from 66 to 75 nucleotides. Additionally, the phylogenetic tree analysis revealed that S. esocinus species that was most alike to S. labiatus. This study provides critical data for phylogenetic analysis of the Schizothoracinae subfamily, which will help to resolve taxonomic difficulties and identify evolutionary links. Detailed mtDNA data are an invaluable resource for studying genetic diversity, population structure, and gene flow. Understanding genetic makeup can help inform conservation plans, identify unique populations, and track genetic variation to ensure effective preservation.

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.

Expression and Variations in EPO Associated with Oxygen Metabolism in Tibetan Sheep

After a long period of adaptive evolution, Tibetan sheep have adapted to the plateau environment in terms of genetics, physiology and biochemistry, but the mechanism of hypoxia adaptation has not been fully elucidated, and the functional genes and molecular mechanisms regulating the hypoxia adaptation of Tibetan sheep need to be further studied. In this study, Tibetan sheep were selected as the research object, and the mRNA expression levels of the hypoxa-related gene EPO in heart, lung, kidney, liver, spleen and longissimus dorsi muscle of Hu sheep (100 m) and Tibetan sheep at different altitudes (2500 m, 3500 m, 4500 m) were assessed by RT-qPCR. The SNPs loci were detected by sequencing and Kompetitive Allele-Specific PCR (KASP) technology, then the correlation between genetic polymorphism and blood gas was analyzed. The results show that the expression of the EPO gene was the highest in the kidney, indicating that the expression of EPO gene had tissue differences. The expression levels of the EPO gene in the heart, lung and liver of Tibetan sheep at a 4500 m altitude were significantly higher than those in Hu sheep (p < 0.05), and the levels in the hearts of Tibetan sheep increased with the increase in altitude. Three mutations were identified in the EPO gene, the SNPs (g.855 A > C) in exon 1 and the SNPs (g.1985 T > G and g.2115 G > C) in exon 4, which were named EPO-SNP1, EPO-SNP2 and EPO-SNP3, respectively, and all three SNPs showed three genotypes. Correlation analysis showed that g.2115 G > C sites were significantly correlated with pO2 (p < 0.05), and haplotype combinations were significantly correlated with pO2 (p < 0.05). Thesee results suggest that the expression of the EPO gene is altitude-differentiated and organ-differentiated, and the EPO gene variants have significant effects on pO2, which may be beneficial to the adaptation of Tibetan sheep to hypoxia stress.

Measurement of Fitness and Predatory Ability of Four Predatory Mite Species in Tibetan Plateau under Laboratory Conditions

Predatory mites are biological control agents used in many countries against various vegetable pests, particularly spider mites. Despite the significant presence of predatory mites in the Tibetan plateau, there is limited research on their potential against spider mites in the area. This study investigated the fitness parameters and performance against TSSM of four predatory, including Amblyseius swirskii (Athias-Henriot) and three species from the genus Neoseiulus (Neoseiulus californicus (McGregor), Neoseiulus barkeri (Hughes), and Neoseiulus cucumeris (Oudemans)), originally collected from fields in the Tibetan Plateau. Compared to the other three predatory species, A. swirskii exhibited the highest fecundity (11.60 ± 0.34) and the highest pre-adult survival rate (83.33 ± 3.33%). Since their juvenile survival rate (SR) was extremely low (13.33% ± 5.77%), most N. barkeri nymphs died before emergence. Compared to the other three predatory mites, A. swirskii showed the highest predation capacity against adult TSSMs at 15 d post-release (14.28 ± 2.24). Based on the results, A. swirskii was the most effective, and N. barkeri was the least effective in controlling two-spotted mites in the Tibetan Plateau among the four species tested in this study. Collectively, these findings imply notable advantages in employing A. swirskii for controlling two-spotted mites in the Tibetan Plateau. This study informs the development of a feasible biological control method based on suitable predatory mite species to manage TSSMs in the Tibetan Plateau.

Micro(nano)plastics in marine medaka: Entry pathways and cardiotoxicity with triphenyltin

The simultaneous presence of micro(nano)plastics (MNPs) and pollutants represents a prevalent environmental challenge that necessitates understanding their combined impact on toxicity. This study examined the distribution of 5 μm (PS-MP5) and 50 nm (PS-NP50) polystyrene plastic particles during the early developmental stages of marine medaka (Oryzias melastigma) and assessed their combined toxicity with triphenyltin (TPT). Results showed that 2 mg/L PS-MP5 and PS-NP50 could adhere to the embryo surface. PS-NP50 can passively enter the larvae and accumulate predominantly in the intestine and head, while PS-MP5 cannot. Nonetheless, both types can be actively ingested by the larvae and distributed in the intestine. 2 mg/L PS-MNPs enhance the acute toxicity of TPT. Interestingly, high concentrations of PS-NP50 (20 mg/L) diminish the acute toxicity of TPT due to their sedimentation properties and interactions with TPT. 200 μg/L PS-MNPs and 200 ng/L TPT affect complement and coagulation cascade pathways and cardiac development of medaka larvae. PS-MNPs exacerbate TPT-induced cardiotoxicity, with PS-NP50 exhibiting stronger effects than PS-MP5, which may be related to the higher adsorption capacity of NPs to TPT and their ability to enter the embryos before hatching. This study elucidates the distribution of MNPs during the early developmental stages of marine medaka and their effects on TPT toxicity, offering a theoretical foundation for the ecological risk assessment of MNPs.

The complete mitochondrial genome of a new species of the genus Schizothorax from Sichuan, China (Cypriniformes: Cyprinidae)

Schizothorax gulinensis sp. nov. is a new species of the genus Schizothorax from Sichuan, China (Cypriniformes: Cyprinidae). In this study, we have first reported the complete mitochondrial genome of S. gulinensis with Illumina sequencing. There were 16,587 nucleotide pairs in the mitochondrial genome (mitogenome) of S. gulinensis, including 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), and 22 transfer RNAs (tRNAs), as well as one non-coding control region (CR). The proportion of nucleotides in mitochondrial genome was 29.67% (A), 25.45% (T), 17.84% (G), 27.05% (C), and A + T content was 55.12%. All PCGs have the same start codon of the standard ATG, excepting for that of NADH dehydrogenase subunit 1 (nad1) which was the ATC, NADH dehydrogenase subunit 5 (nad5) which was the ATT and cytochrome c oxidase 1 (cox1) which was the ATC. Phylogenetic analysis results supported that S. gulinensis was closely related to Schizothorax grahami. The complete mitochondrial sequence of S. gulinensis will contribute to mitochondrial genome database and provide useful resources for population genetics and evolution analyses.

Comparative analyses reveal potential genetic mechanisms for high-altitude adaptation of Schizopygopsis fishes based on chromosome-level genomes

The schizothoracine fishes, widely distributed in the Qinghai-Tibetan Plateau and its adjacent areas, are considered as ideal models for investigation of high-altitude adaptation. Schizophygopsis are one group of the highly specialized schizothoracine fishes, and the genetic basis for their high-altitude adaptation is poorly understood. In this study, we performed comparative genomics analyses to investigate the potential genetic mechanisms for high-altitude adaptation of Schizopygopsis malacanthus and Schizopygopsis pylzovi based on the chromosome-level genomes. Functional enrichment analysis revealed that many expanded gene families in Schizopygopsis were associated with immune response while many contracted gene families were functionally associated with olfaction. Among the 123 positively selected genes (PSGs), angpt2a was detected in HIF-1 signaling pathway and possibly related to the hypoxia adaptation of Schizopygopsis. Furthermore, two PSGs cox15 and ndufb10 were distributed in thermogenesis, and there was a Schizopygopsis-specific missense mutation in cox15 (Gln115Glu), which possibly contributed to the cold temperature adaptation of the Schizopygopsis. Kyoto Encyclopedia of Genes and Genomes enrichment of the PSGs revealed three significant pathways including metabolic pathways, cell cycle, and homologous recombination and Gene Ontology enrichment analysis of the PSGs revealed several categories associated with DNA repair, cellular response to DNA damage stimulus, and metabolic process. Chromosome-scale characterization of olfactory receptor (OR) repertoires indicated that Schizopygopsis had the least number of OR genes, and the OR gene contraction was possibly caused by the limited food variety and the environmental factors such as lower air pressure, lower humidity, and lower temperature. Our study will help expand our understanding of the potential adaptive mechanism of Schizopygopsis to cope with the high-altitude conditions.

Exposure to manganese (II) chloride induces developmental toxicity, oxidative stress and inflammatory response in Marine medaka (Oryzias melastigma) embryos

Manganese (Mn) is an essential metal for organisms, but high levels can induce serious toxicity. To date, the toxic mechanism of Mn to marine fish is still poorly understood. In the present study, Oryzias melastigma embryos were exposed to different concentrations of MnCl₂ (0-152.00 mg/L) to investigate its effect on early development. The results showed that exposure to MnCl₂ caused developmental toxicity to embryos, including increased heart rate, delayed hatching time, decreased hatching rate and increased malformation rate. MnCl₂ exposure could induce oxidative stress in O. melastigma embryos, as indicated by increased the contents of malondialdehyde (MDA) and the activities of the antioxidant enzymes (superoxide dismutase (SOD), glutathione peroxidase (GPx) and catalase (CAT)). The heart might be an important target organ for MnCl₂ because of cardiac malformations and disruption in the expression of cardiac development-related genes (ATPase, epo, fg8g, cox1, cox2, bmp4 and gata4). In addition, the expression levels of stress- (omTERT and p53) and inflammation-related genes (TNFα and il1β) were significantly up-regulated, suggesting that MnCl₂ can trigger stress and inflammatory response in O. melastigma embryos. In conclusion, this study demonstrated that MnCl₂ exposure can induce developmental toxicity, oxidative stress and inflammatory response in O. melastigma embryos, providing insights into the toxic mechanism of Mn to the early development of marine fish.

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.

Characterization of EPO H131S as a key mutation site in the hypoxia-adaptive evolution of Gymnocypris dobula

Erythropoietin (EPO) is a glycoprotein hormone involved in proerythropoiesis, antioxidation, and antiapoptosis. It also contributes to cellular immune function in high-altitude species, such as the schizothoracine fish Gymnocypris dobula. Six mutation sites previously identified in EPO from G. dobula (GD-EPO) were injected into zebrafish embryos, and their effects were compared with EPO from the low-altitude schizothoracine Schizothorax prenanti. The key mutation site in GD-EPO was identified as H131S. Under hypoxic conditions, the levels of superoxide dismutase and malondialdehyde were decreased, whereas that of nitric oxide was increased in zebrafish injected with GD-EPO compared with those injected with S. prenanti-EPO (SP-EPO). The results suggest that EPO in high-altitude schizothoracine species is both antioxidative and antiapoptotic, driven by the H131S mutation site. Thus, this enhanced the ability of this species to adapt to the high-altitude hypoxic environment. These results provide a basis for investigating further the hypoxia adaptation mechanisms of teleosts.

Integrated Metabolomics and Transcriptomic Analysis of Hepatopancreas in Different Living Status Macrobrachium nipponense in Response to Hypoxia

As the basic element of aerobic animal life, oxygen participates in most physiological activities of animals. Hypoxia stress is often the subject of aquatic animal research. Macrobrachium nipponense, an economically important aquatic animal in southern China, has been affected by hypoxia for many years and this has resulted in a large amount of economic loss due to its sensitivity to hypoxia; Metabolism and transcriptome data were combined in the analysis of the hepatopancreas of M. nipponense in different physiological states under hypoxia; A total of 108, 86, and 48 differentially expressed metabolites (DEMs) were found in three different comparisons (survived, moribund, and dead shrimps), respectively. Thirty-two common DEMs were found by comparing the different physiological states of M. nipponense with the control group in response to hypoxia. Twelve hypoxia-related genes were identified by screening and analyzing common DEMs. GTP phosphoenolpyruvate carboxykinase (PEPCK) was the only differentially expressed gene that ranked highly in transcriptome analysis combined with metabolome analysis. PEPCK ranked highly both in transcriptome analysis and in combination with metabolism analysis; therefore, it was considered to have an important role in hypoxic response. This manuscript fills the one-sidedness of the gap in hypoxia transcriptome analysis and reversely deduces several new genes related to hypoxia from metabolites. This study contributes to the clarification of the molecular process associated with M. nipponense under hypoxic stress.

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.

Transcriptome analysis of hepatopancreas from different living states oriental river prawn (Macrobrachium nipponense) in response to hypoxia

As an important economical freshwater prawn, Macrobrachium nipponense has difficulty with adapting to hypoxia. In this study, comparative transcriptome analysis was used for the first time to explore the differences between different living states of Macrobrachium nipponense under hypoxia. A total of 94.22 Gb clean reads were obtained and assembled into 54,688 unigenes. A total of 224, 266, and 750 differently expressed genes were found in the comparison of the control and death groups, the control and moribund groups, and the control and survived groups, respectively. Three signal pathways closely related to hypoxia were found by enriching of the signal pathways in three comparison groups. In addition, much attention was focused on the differential genes in these pathways. Oxidative stress related genes, such as 70 kDa heat shock protein, phosphoenolpyruvate carboxykinase and cyclooxygenase were differentially expressed in different comparisons. After comparing with previous studies, cyclooxygenase was found to be an important hypoxia-related gene that is fully involved in the hypoxic response. Interestingly, two new genes with no Nr annotation were found in this manuscript. This manuscript will enrich our understanding of oxidative stress response to hypoxia and provide a theoretical basis for the subsequent solution of apoptosis caused by hypoxia.

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.

Genome of Tetraploid Fish Schizothorax o'connori Provides Insights into Early Re-diploidization and High-Altitude Adaptation

Whole-genome duplications (WGDs) of Schizothoracinae are believed to have played a significant role in speciation and environmental adaptation on the Qinghai-Tibet Plateau (QTP). Here, we present a genome for Schizothorax o'connori, a QTP endemic fish and showed the species as a young tetraploid with a recent WGD later than ∼1.23 mya. We exhibited that massive insertions between duplicated genomes caused by transposon bursts could induce mutagenesis in adjacent sequences and alter the expression of neighboring genes, representing an early re-diploidization process in a polyploid genome after WGD. Meanwhile, we found that many genes involved in DNA repair and folate transport/metabolism experienced natural selection and might contribute to the environmental adaptation of this species. Therefore, the S. o'connori genome could serve as a young tetraploid model for investigations of early re-diploidization in polyploid genomes and offers an invaluable genetic resource for environmental adaptation studies of the endemic fish of the QTP.

Teratogenic effects of environmentally relevant concentrations of phenanthrene on the early development of marine medaka (Oryzia melastigma)

Polycyclic aromatic hydrocarbons (PAHs) are ubiquitous pollutants in marine environments and have arouse great concern since they pose adverse effects to marine ecosystem. To determine the potential impacts of environmentally relevant PAHs on early life stages of marine fish, this study exposed embryos of marine medaka (Oryzias melastigma) to 0, 2, 10, 50, and 250 μg/L of phenanthrene (Phe), one of the most abundant PAHs. The results demonstrated that Phe exposure decreased hatching rates, delayed hatching time of embryos, and increased deformity rate of newly-hatched larvae. Exposure to 10 and 50 μg/L Phe decreased the survival rate of marine medaka larvae at 28 days post-fertilization (dpf), and no embryo successfully hatched in 250 μg/L Phe exposure group. Morphology results showed that 10, 50, and 250 μg/L Phe exposure significantly retarded the development of embryos, and 2, 10, and 50 μg/L caused yolk sac edema and pericardial edema in newly-hatched larvae, indicating that low concentrations of Phe could induce developmental cardiac toxicity. Furthermore, the changes in the expression of heart development-related genes were determined, and the results showed that Phe-induced cardiac malformation might be related with fgf8, bmp4, smyd1, ATPase and gata4 genes. Overall, environmentally relevant PAHs could disrupt heart morphogenesis and hatching process of marine medaka, which might have profound consequences for sustainability of fish population.

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.

Hypoxia-inducible factor 1α from a high-altitude fish enhances cytoprotection and elevates nitric oxide production in hypoxic environment

Hypoxia-inducible factors (HIFs) are master transcription factor regulating hypoxic responses in vertebrates. Species of Schizothoracine, a sub-family of cyprinidae, are highly endemic to the hypoxic Qinghai-Tibetan Plateau (QTP). What roles the HIFs play in hypoxic adaptation in the Schizothoracine fish is little known. In this study, the HIF-1α/B gene from Gymnocypris dobula (Gd) was characterized. The predicted protein for Gd-HIF-1α/B contains the main domains (bHLH, PAS, PAC, ODD, N-TAD, and C-TAD). Moreover, a specific mutation that the proline hydroxylation motif (LXXLAP) mutated into PxxLAP was observed in Gd-HIF-1α/B CODD domain, which may lead to changes in the function. To clarify whether HIF-1α/B of G. dobula possesses hypoxic adaptive features, Gd-HIF1α/B and Schizothorax prenanti-HIF1α/B (Sp-HIF1α/B) were cloned into an expression vector and transfected into 293T cells. Cell viability was found to be significantly higher in cells transfected with Gd-HIF-1α/B than those transfected with Sp-HIF-1α/B under hypoxic conditions. In addition, G. dobula HIF-1α/B showed stronger activity in transactivating the expression of nitric oxide (NO)-synthesizing enzyme, NOS2B under hypoxia stresses than the orthologous gene from S. prenanti, which were accompanied with upregulated expressions of NOS2B in heart of G. dobula, which may attribute to elevated NO levels detected in G. dobula than the lower land species. These results indicated that the HIF-1α plays an important role in mediating the iNOS signaling system in the process of evolutionary adaptation of the Schizothoracine to the highland environment.

The Genome Landscape of Tibetan Sheep Reveals Adaptive Introgression from Argali and the History of Early Human Settlements on the Qinghai-Tibetan Plateau

Tibetan sheep are the most common and widespread domesticated animals on the Qinghai-Tibetan Plateau (QTP) and have played an essential role in the permanent human occupation of this high-altitude region. However, the precise timing, route, and process of sheep pastoralism in the QTP region remain poorly established, and little is known about the underlying genomic changes that occurred during the process. Here, we investigate the genomic variation in Tibetan sheep using whole-genome sequences, single nucleotide polymorphism arrays, mitochondrial DNA, and Y-chromosomal variants in 986 samples throughout their distribution range. We detect strong signatures of selection in genes involved in the hypoxia and ultraviolet signaling pathways (e.g., HIF-1 pathway and HBB and MITF genes) and in genes associated with morphological traits such as horn size and shape (e.g., RXFP2). We identify clear signals of argali (Ovis ammon) introgression into sympatric Tibetan sheep, covering 5.23-5.79% of their genomes. The introgressed genomic regions are enriched in genes related to oxygen transportation system, sensory perception, and morphological phenotypes, in particular the genes HBB and RXFP2 with strong signs of adaptive introgression. The spatial distribution of genomic diversity and demographic reconstruction of the history of Tibetan sheep show a stepwise pattern of colonization with their initial spread onto the QTP from its northeastern part ∼3,100 years ago, followed by further southwest expansion to the central QTP ∼1,300 years ago. Together with archeological evidence, the date and route reveal the history of human expansions on the QTP by the Tang-Bo Ancient Road during the late Holocene. Our findings contribute to a depth understanding of early pastoralism and the local adaptation of Tibetan sheep as well as the late-Holocene human occupation of the QTP.

Function-related Drivers of Skull Morphometric Variation and Sexual Size Dimorphism in a Subterranean Rodent, Plateau Zokor (Eospalax baileyi)

Sexual dimorphism is prevalent in most living organisms. The difference in size between sexes of a given species is generally known as sexual size dimorphism (SSD). The magnitude of the SSD is determined by Rensch's rule where size dimorphism increases with increasing body size when the male is the larger sex and decreases with increasing average body size when the female is the larger sex. The unique underground environment that zokors (Eospalax baileyi) live under in the severe habitat of the Qinghai‐Tibetan Plateau (QTP) could create SSD selection pressures that may or may not be supported by Rensch's rule, making this scientific question worthy of investigation. In this study, we investigated the individual variation between sexes in body size and SSD of plateau zokors using measurements of 19 morphological traits. We also investigated the evolutionary mechanisms underlying SSD in plateau zokors. Moreover, we applied Rensch's rule to all extant zokor species. Our results showed male‐biased SSD in plateau zokors: The body‐ and head‐related measurements were greater in males than in females. Linear regression analysis between body length, body weight, and carcass weight showed significant relationships with some traits such as skull length, lower incisor length, and tympanic bulla width, which might support our prediction that males have faster growth rates than females. Further, the SSD pattern corroborated the assumption of Rensch's rule in plateau zokors but not in the other zokor species. Our findings suggest that the natural underground habitat and behavioral differences between sexes can generate selection pressures on male traits and contribute to the evolution of SSD in plateau zokors.

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.

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.