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

Genomic signatures associated with recurrent scale loss in cyprinid fish

Scale morphology represents a fundamental feature of fish and a key evolutionary trait underlying fish diversification. Despite frequent and recurrent scale loss throughout fish diversification, comprehensive genome-wide analyses of the genomic signatures associated with scale loss in divergent fish lineages remain scarce. In the current study, we investigated genome-wide signatures, specifically convergent protein-coding gene loss, amino acid substitutions, and cis-regulatory sequence changes, associated with recurrent scale loss in two divergent Cypriniformes lineages based on large-scale genomic, transcriptomic, and epigenetic data. Results demonstrated convergent changes in many genes related to scale formation in divergent scaleless fish lineages, including loss of P/Q-rich scpp genes (e.g. scpp6 and scpp7), accelerated evolution of non-coding elements adjacent to the fgf and fgfr genes, and convergent amino acid changes in genes (e.g. snap29) under relaxed selection. Collectively, these findings highlight the existence of a shared genetic architecture underlying recurrent scale loss in divergent fish lineages, suggesting that evolutionary outcomes may be genetically repeatable and predictable in the convergence of scale loss in fish.

Morphological and Histological Analysis of the Gastrointestinal Systems in Triplophysa strauchii and Triplophysa tenuis: Insights into Digestive Adaptations

Fish are vital for material cycling and energy flow in aquatic ecosystems. The genus Triplophysa, with over 100 known species, is significant in the Central Asian highlands' freshwater ecosystems. T. strauchii and T. tenuis, as representatives, occupy distinct ecological niches and face challenges from climate change and human activities. There is a lack of research on Triplophysa fishes' digestive systems, especially comparative studies, so this research aims to fill this gap. In September 2024, 40 samples of T. strauchii were collected from Sayram Lake and 40 samples of T. tenuis were collected from the Muzat River in Xinjiang. After acclimation, morphological observations (measuring fish and digestive tract parameters) and histological analyses (paraffin sectioning, HE staining, and microscopy) were carried out. The data were sorted in Excel and analyzed with an independent samples t-test in SPSS 27.0. Morphologically, T. strauchii has an obtuse snout, terminal mouth, specific upper lip papillae, and an S-shaped intestine about (1.45 ± 0.11) times its body length, while T. tenuis has an arc-shaped subterminal mouth, fringed papillae, and a spiral-shaped intestine around (0.82 ± 0.09) times its body length. Both possess a digestive tract, glands, and a hepatopancreas attached to the mesentery. Histologically, a large number of club cells were found in the oropharyngeal cavities of both species; their secretions have an adhesive effect on food, aiding food selection. Their digestive systems vary in structure and cell composition: the oropharyngeal cavity has three layers; the esophagus has four layers with more goblet cells in T. strauchii; the stomach has three regions without goblet cells and a thicker muscular layer in T. strauchii; the intestinal wall has four layers with different villi and goblet cell distributions; the hepatopancreas has lobules; and T. strauchii has a typical portal area. In conclusion, this study systematically compared the gastrointestinal systems of T. strauchii and T. tenuis for the first time, revealing significant structural differences related to their niches and feeding patterns as adaptations to specific environments. It fills the research gap, provides a basis for exploring fish ecological adaptation and environmental impacts on digestion, offers new ideas for Triplophysa protection strategies, and guides fish evolutionary biology research and Triplophysa resource protection and utilization.

Mitochondrial Genomes of Six Discogobio Species (Teleostei, Cyprinidae) and Their Phylogenetic Analysis

Discogobio is an important small freshwater economic fish in Southwestern China. In this study, we determined the complete mitochondrial genome (mitogenome) sequences of six Discogobio species by conventional overlapping PCR and Sanger sequencing. The mitogenomes were 16,591-16,605 bp in length and contained 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and a noncoding region, which was consistent with previously studied Labeoninae mitogenomes. The phylogenetic analysis revealed that the subfamily Labeoninae was divided into four major clades. Notably, Discogobio was observed to be a non-monophyletic group and mixed with Discocheilus. The divergence time of Discogobio and Discocheilus can be placed during the time of the most recent common ancestor in the late Miocene (7.80 Mya). The results of the selection pressure analyses indicated that all Discogobio fishes exhibited Ka/Ks ratios < 1, suggesting that mitochondrial function in this genus was subjected to strong purifying selection and adapted to different environments. These mitogenomes will facilitate further studies in phylogeny, taxonomy, and evolutionary biology related to the subfamily Labeoninae.

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.

Terrestrial Adaptation in Chelonoidis vicina as Revealed Based on Analysis of the Complete Mitochondrial Genome

BACKGROUND/OBJECTIVES

Mitochondrial genomes are widely used in phylogenetics and evolutionary and ecological research.

METHODS

In this study, the newest mitochondrial genome of Chelonoidis vicina was assembled and annotated. The comparative mitochondrial genome and selection pressure analyses were used to examine the terrestrial adaptive evolution characteristics of C. vicina and other terrestrial reptiles.

RESULTS

The results reveal that the mitochondrial genome of the tortoise C. vicina is consistent with that of other tortoise species, comprising 13 protein-coding genes (PCGs), 2 rRNAs, 22 tRNAs, and 1 noncoding control region (CR). The analysis of selection pressure reveals the presence of positive selection sites in the COX2, COX3, Cytb, ND3, ND4, ND4L, ND5, and ND6 genes of terrestrial reptiles. Of these, the COX2 and ND3 genes exhibited faster evolutionary rates. The mitochondrial genome structure of C. vicina is consistent with that of different terrestrial reptiles. The positive selection sites of COX2 and ND3 in terrestrial reptiles are closely related to a change in mitochondrial energy metabolism, which is possibly related to terrestrial adaptability.

CONCLUSIONS

The results of this study provide new insights into the adaptive evolution of C. vicina to terrestrial niches from a mitogenomic perspective, as well as genetic resources for the protection of C. vicina.

Revision of the loach genus Claea Kottelat, 2010 (Teleostei: Nemacheilidae) in China, with a description of a new species from the Chang-Jiang basin

A taxonomic revision of the nemacheilid genus Claea from China is provided based on molecular and morphological data. Three species are identified in this genus: C. dabryi (type species), C. wulongensis (a species originally recognized in Triplophysa), and C. minibarba sp. nov., here described from the mid-upper Chang-Jiang basin in northwestern Hubei Province, South China. This new species differs from the two congeneric species in having a deeper body, further from C. dabryi in having shorter outer rostral and maxillary barbels, and further from C. wulongensis in possessing a higher count of vertebrae, a deeper caudal-peduncle, and a narrower interorbital space. Its validity is corroborated by its monophyly recovered in a mitochondrial cytochrome b gene-based phylogenetic analysis and its significant genetic divergence of this gene with the two congeneric species.

Characterization and phylogenetic implications of newly sequenced mitochondrial genomes of cobitid fish Acantopsis Rungthipae (Boyd, Nithirojpakdee & Page, 2017)

BACKGROUND

Acantopsis rungthipae has significant ornamental and ecological value. This study aimed at structurally characterizing the A. rungthipae mitochondrial genome and elucidate its phylogenetic position in Cobitidae.

METHODS AND RESULTS

High-throughput sequencing technology was used to obtain the complete sequence of the mitochondrial genome of A. rungthipae and reconstruct a Cobitidae phylogenetic tree based on the sequence of 13 protein-coding genes. The entire mitochondrial genome of A. rungthipae was 16,600 bp, containing 22 tRNA genes, 13 protein-coding genes, 2 rRNAs, and 2 non-coding regions (D-loop and OL). The base composition showed a significant AT preference, with the highest A + T content (67.1%) in the D-loop region. Among the protein-coding genes, 12 had ATG as a typical starting codon, while only COXI had GTG as a special starting codon. Twenty-one of the tRNA genes exhibited clover structure, and only tRNA-Ser (GCT) could not fold into a clover structure because of the absence of DHU arms. The phylogenetic tree was reconstructed using the Bayesian and maximum likelihood methods and showed that A. rungthipae and Acantopsis choirorhynchos converged into one branch, and their phylogenetic relationships were relatively close.

CONCLUSIONS

Our findings supplement basic data on the A. rungthipae mitochondrial genome and deepen the understanding of the evolutionary relationships of the genus Acantopsis. Clarifying the evolutionary relationships between different species in Acantopsis lays a solid foundation for subsequent research on fish adaptation and selection pressure.

Description of three new species of Gyrodactylus von Nordmann, 1832 (Monopisthocotylea: Gyrodactylidae) on Triplophysa species (Nemacheilidae) from Qinghai-Tibet Plateau

Three new species of Gyrodactylus were identified from the body surface of the Triplophysa species from the Qinghai-Tibet Plateau, Gyrodactylus triplorienchili n. sp. on Triplophysa orientalis in northern Tibet, G. yellochili n. sp. on T. sellaefer and T. scleroptera and G. triplsellachili n. sp. on T. sellaefer and T. robusta in Lanzhou Reach of the Yellow River. The three newly identified species share the nemachili group species' characteristic of having inturning hamulus roots. Gyrodactylus triplorienchili n. sp. shared a quadrate sickle heel and a thin marginal hook sickle, two morphological traits that set them apart from G. yellochili n. sp. However, they may be identified by the distinct shapes of the sickle base and marginal hook sickle point. Gyrodactylus triplsellachili n. sp. had much larger opisthaptoral hard part size than the other two species. The three new species show relatively low interspecific differences of 2.9-5.3% p-distance for ITS1-5.85-ITS2 rDNA sequences. Phylogenetic analysis indicated that the three new species formed a well-supported monophyletic group (bp = 99) with the other nemachili group species.

The complete mitochondrial genome of Triplophysa erythraea (Huang et al. 2019) (Cypriniformes, Nemacheilidae)

In this study, the complete mitochondrial genome of Triplophysa erythraea was determined for the first time. Results showed the mitogenome was 16 565 bp, including 2 ribosomal RNA genes, 22 transfer RNA genes, 13 protein-coding genes and 2 non-coding regions. The overall base composition of its mitochondrial genome was 28.9% A, 25.2% T, 28.7% C and 17.2% G. Phylogenetic tree revealed that Triplophysa erythraea had the closest relationship with Triplophysa xiangxiensis, which was also found in Xiangxi Tujia and Miao Nationality Autonomous Prefecture, Hunan province. In general, this study provided valuable information for conservation genetics analyses of Triplophysa erythraea and further displayed the evolution of species within the genus Triplophysa.

The complete mitochondrial genome of Triplophysa grahami Regan 1906 (Cypriniformes: Nemacheilidae) and phylogenetic analysis

Triplophysa grahami Regan 1906 is a member of the family Nemacheilidae, Cypriniformes, and native loach in Yunnan. In this study, the complete mitochondrial genome (mitogenome) of T. grahami Regan 1906 was firstly reported and analyzed. The mitogenome of T. grahami Regan 1906 is 16,566 bp in length, including 13 protein-coding genes (PCGs), 22 transfer RNAs (tRNAs), two ribosomal RNAs (rRNAs), and one control region (D-loop). The arrangement and orientation of protein coding genes and RNAs in T. grahami Regan 1906 are identical to other species of Nemacheilidae. The base composition of T. grahami Regan 1906 mitogenome was 29.25% A, 28.55% T, 25.03% C, and 17.17% G. The phylogenetic analysis based on the mitogenome showed that T. grahami Regan 1906 belongs to the clade of genus Triplophysa and the monophyly of Triplophysa is identified. This study contributed valuable genetic data for T. grahami Regan 1906 and explored the phylogenetic relationships in Nemacheilidae.

Analysis of genetic diversity among three Triplophysa tenuis populations by RAD-seq

To investigate the genetic diversity of Triplophysa tenuis in the Shule River Basin of Gansu province, three populations were sequenced via RAD-seq technology. Twenty-nine microsatellite (SSR) markers with polymorphisms were finally screened to access the genetic diversity among the populations, of which 15 had high polymorphisms. The quantity of the alleles detected in the three populations of T. tenuis varied from 2 to 24. The locus with the most alleles was SSRC1, which had 24 alleles. Among the 29 SSRs, the range of effective allele number, observed heterozygosity, expected heterozygosity, and polymorphic information content were 1.246-16.615, 0.222-1, 0.198-0.940, and 0.178-0.937, respectively. Most of the identified loci were in the Hardy-Weinberg equilibrium. Analysis of the population structure revealed that the Yumen and Changma populations shared the same origin, while the Qiaowan population was different from them. The developed SSR markers discovered in this study will contribute to the conservation research on T. tenuis and the conservation of the fishery resources of the Shule River, providing scientific guidance for the development and utilization of T. tenuis resources and environmental protection.

Effect of locomotor preference on the evolution of mitochondrial genes in Bovidae

Locomotor preferences and habitat types may drive animal evolution. In this study, we speculated that locomotor preference and habitat type may have diverse influences on Bovidae mitochondrial genes. We used selection pressure and statistical analysis to explore the evolution of mitochondrial DNA (mtDNA) protein-coding genes (PCGs) from diverse locomotor preferences and habitat types. Our study demonstrates that locomotor preference (energy demand) drives the evolution of Bovidae in mtDNA PCGs. The habitat types had no significant effect on the rate of evolution in Bovidae mitochondrial genes. Our study provides deep insight into the adaptation of Bovidae.

Mitochondrial Genomic Evidence of Selective Constraints in Small-Bodied Terrestrial Cetartiodactyla

Body size may drive the molecular evolution of mitochondrial genes in response to changes in energy requirements across species of different sizes. In this study, we perform selection pressure analysis and phylogenetic independent contrasts (PIC) to investigate the association between molecular evolution of mitochondrial genome protein-coding genes (mtDNA PCGs) and body size in terrestrial Cetartiodactyla. Employing selection pressure analysis, we observe that the average non-synonymous/synonymous substitution rate ratio (ω) of mtDNA PCGs is significantly reduced in small-bodied species relative to their medium and large counterparts. PIC analysis further confirms that ω values are positively correlated with body size (R2 = 0.162, p = 0.0016). Our results suggest that mtDNA PCGs of small-bodied species experience much stronger purifying selection as they need to maintain a heightened metabolic rate. On the other hand, larger-bodied species may face less stringent selective pressures on their mtDNA PCGs, potentially due to reduced relative energy expenditure per unit mass. Furthermore, we identify several genes that undergo positive selection, possibly linked to species adaptation to specific environments. Therefore, despite purifying selection being the predominant force in the evolution of mtDNA PCGs, positive selection can also occur during the process of adaptive evolution.

Mitochondrial genome data provide insights into the phylogenetic relationships within Triplophysadalaica (Kessler, 1876) (Cypriniformes, Nemacheilidae)

Due to the detrimental effect of formaldehyde on DNA, ethanol has replaced formalin as the primary preservative for animal specimens. However, short-term formalin fixation of specimens might be applied during field collection. In an increasing number of studies, DNA extraction and sequencing have been successfully conducted from formalin-fixed specimens. Here the DNA from five specimens of Triplophysadalaica (Kessler, 1876) were extracted and performed high-throughput sequencing. Four of the specimens underwent short-term fixation with formalin and were subsequently transferred to ethanol. One was continuously stored in ethanol. No significant difference of DNA quality and amount were observed among these samples. Followed by assembly and annotation, five mitochondrial genomes ranging in length from 16,569 to 16,572 bp were obtained. Additionally, previously published data of other individuals or species were included to perform phylogenetic analyses. In the reconstructed trees, all eight individuals of T.dalaica form a monophyletic group within the Triplophysa branch. The group is divided into three clades: (1) samples from the Yellow River, (2) those from the Yangtze River, and (3) those from the Haihe River, and the Lake Dali Nur. This study sheds initial light on the phylogeographic relationships among different populations of T.dalaica, and will support the research about its evolutionary history in the future.

Adaptive evolution of mitochondrial genomes in Triplophysa cavefishes

Extreme conditions in caves pose survival challenges for cave dwellers, who gradually develop adaptive survival features. Cavefishes are one of the most successful animals among cave dwellers. Triplophysa cavefishes are an important group of cavefishes, and they show remarkable adaptability to the extreme environments of caves. However, there is a limited understanding of their adaptation mechanisms. In this study, eight complete mitochondrial genomes of Triplophysa cavefishes were newly obtained, and their genomic characteristics, including the base composition, base bias, and codon usage, were analyzed. Phylogenetic analysis was carried out based on 13 mitochondrial protein-coding genes from 44 Nemacheilidae species. This showed that Triplophysa cavefishes and non-cavefishes separate into two reciprocally monophyletic clades, suggesting a single origin of the cave phenotype. Positive selection analysis strongly suggested that the selection pressure in cavefishes is higher than that in non-cavefishes. Furthermore, the ND5 gene in cavefishes showed evidence of positive selection, which suggests that the gene may play an important role in the adaptation of cavefishes to the cave environment. Protein structure analysis of the ND5 subunit implied that the sites of positive selection in cavefishes might allow them to acquire lower ND5 protein stability, compared to that in non-cavefishes, which might help the accumulation of nonsynonymous (mildly deleterious) mutations. Together, our study revealed the genetic signatures of cave adaptation in Triplophysa cavefishes from the perspective of energy metabolism.

The complete mitochondrial genome of Triplophysa nanpanjiangensis Zhu and Cao 1988 (Cypriniformes: Nemacheilidae)

The genus Triplophysa is an ideal taxon for understanding geological evolution. In this study, for the first time, we report the complete mitochondrial genome of T. nanpanjiangensis Zhu and Cao 1988 using the Nanopore sequencing. It is a circular genome with a length of 16558 bp, comprising 22 tRNAs, 13 protein-coding genes (PCGs), two rRNAs, and one non-coding control region. The phylogenetic tree demonstrates that T. nanpanjiangensis is sister to Triplophysa zhenfengensis and placed within the genus Triplophysa. Our mitogenomic studies provide a new pathway for understanding the molecular phylogeny of the genus Triplophysa.

Age, growth, and energy storage of the subterranean fish Triplophysa rosa (Cypriniformes: Nemacheilidae) from Chongqing, China

BACKGROUND

This study explores the age, growth, and energy storage of Triplophysa rosa, a troglobitic cavefish. A total of 102 wild T. rosa specimens were collected in Wulong County, Chongqing, China, between 2018 and 2022, with otoliths used for age determination.

RESULTS

The earliest mature individuals were determined to be 4.8 years old, while the maximum ages for females and males were estimated at 15.8 years and 12.2 years, respectively. The length (L, cm)-weight (W, g) relationship was found to be the same for both sexes, following the eq. W = 0.0046 L3.03. Von Bertalanffy growth models were applied to the total length-at-age data, resulting in an asymptotic length of 23.4 cm and a K-parameter of 0.060 year-1. The body content of protein, ash, and glycogen did not show a significant correlation with the total length of T. rosa. However, both lipid and energy content exhibited a significant increase with total length. The lipid content ranged from 40.5 to 167.1 mg g-1, while the energy content ranged from 4.50 to 11.39 kJ g-1, indicating high storage features of T. rosa.

CONCLUSIONS

The results affirm that T. rosa exhibits life traits conducive to its population dynamics in cave conditions, characterized by slow growth, small size, and high lipid energy storage.

Analysis of the Complete Mitochondrial Genome of Pteronura brasiliensis and Lontra canadensis

P. brasiliensis and L. canadensis are two otter species, which successfully occupied semi-aquatic habitats and diverged from other Mustelidae. Herein, the full-length mitochondrial genome sequences were constructed for these two otter species for the first time. Comparative mitochondrial genome, selection pressure, and phylogenetic independent contrasts (PICs) analyses were conducted to determine the structure and evolutionary characteristics of their mitochondrial genomes. Phylogenetic analyses were also conducted to confirm these two otter species' phylogenetic position. The results demonstrated that the mitochondrial genome structure of P. brasiliensis and L. canadensis were consistent across Mustelidae. However, selection pressure analyses demonstrated that the evolutionary rates of mitochondrial genome protein-coding genes (PCGs) ND1, ND4, and ND4L were higher in otters than in terrestrial Mustelidae, whereas the evolutionary rates of ND2, ND6, and COX1 were lower in otters. Additionally, PIC analysis demonstrated that the evolutionary rates of ND2, ND4, and ND4L markedly correlated with a niche type. Phylogenetic analysis showed that P. brasiliensis is situated at the base of the evolutionary tree of otters, and then L. canadensis diverged from it. This study suggests a divergent evolutionary pattern of Mustelidae mitochondrial genome PCGs, prompting the otters' adaptation to semi-aquatic habitats.

Mitochondrial DNA sequencing of Kehilan and Hamdani horses from Saudi Arabia

The Arabian horse breed is well known for its purity and played a key role in the genetic improvement of other horses worldwide. The mitochondrial genome plays a vital role in maternal inheritance and it's helpful to evaluate its genetic diversity and conservation. It has higher mutation rates than nuclear DNA in vertebrates and therefore reveals phylogenetic relationships and haplotypes. In this study, the mitochondrial genome mutations in two Saudi horse strains, Kehilan and Hamdani demonstrated various changes in the gene and amino acid levels and included two other Saudi horses (Hadban and Seglawi) from the previous study for phylogenetic comparison. The whole mitochondrial genome sequencing resulted in intra and inter mtDNA variations between the studied horses. Interestingly, the Hamdani horse has nucleotide substitutions similar to those of the Hadban horse, which is reflected in the phylogenetic tree as a significantly close relationship. This type of study provides a better understanding of mitogenome structure and conservation of livestock species genetic data.

Adaptive evolution characteristics of mitochondrial genomes in genus Aparapotamon (Brachyura, Potamidae) of freshwater crabs

BACKGROUND

Aparapotamon, a freshwater crab genus endemic to China, includes 13 species. The distribution of Aparapotamon spans the first and second tiers of China's terrain ladder, showing great altitudinal differences. To study the molecular mechanisms of adaptive evolution in Aparapotamon, we performed evolutionary analyses, including morphological, geographical, and phylogenetic analyses and divergence time estimation. We sequenced the mitogenomes of Aparapotamon binchuanense and Aparapotamon huizeense for the first time and resequenced three other mitogenomes of Aparapotamon grahami and Aparapotamon gracilipedum. These sequences were combined with NCBI sequences to perform comparative mitogenome analysis of all 13 Aparapotamon species, revealing mitogenome arrangement and the characteristics of protein-coding and tRNA genes.

RESULTS

A new species classification scheme of the genus Aparapotamon has been detected and verified by different aspects, including geographical, morphological, phylogenetics and comparative mitogenome analyses. Imprints from adaptive evolution were discovered in the mitochondrial genomes of group A, including the same codon loss at position 416 of the ND6 gene and the unique arrangement pattern of the tRNA-Ile gene. Multiple tRNA genes conserved or involved in adaptive evolution were detected. Two genes associated with altitudinal adaptation, ATP8 and ND6, which experienced positive selection, were identified for the first time in freshwater crabs.

CONCLUSIONS

Geological movements of the Qinghai-Tibet Plateau and Hengduan Mountains likely strongly impacted the speciation and differentiation of the four Aparapotamon groups. After some group A species dispersed from the Hengduan Mountain Range, new evolutionary characteristics emerged in their mitochondrial genomes, facilitating adaptation to the low-altitude environment of China's second terrain tier. Ultimately, group A species spread to high latitudes along the upper reaches of the Yangtze River, showing faster evolutionary rates, higher species diversity and the widest distribution.

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.

How voles adapt to subterranean lifestyle: Insights from RNA-seq

Life under the earth surface is highly challenging and associated with a number of morphological, physiological and behavioral modifications. Subterranean niche protects animals from predators, fluctuations in environmental parameters, but is characterized by high levels of carbon dioxide and low levels of oxygen and implies high energy requirements associated with burrowing. Moreover, it lacks most of the sensory inputs available above ground. The current study describes results from RNA-seq analysis of four subterranean voles from subfamily Arvicolinae: Prometheomys schaposchnikowi, Ellobius lutescens, Terricola subterraneus, and Lasiopodomys mandarinus. Original RNA-seq data were obtained for eight species, for nine species, SRA data were downloaded from the NCBI SRA database. Additionally assembled transcriptomes of Mynomes ochrogaster and Cricetulus griseus were included in the analysis. We searched for the selection signatures and parallel amino acid substitutions in a total of 19 species. Even within this limited data set, we found significant changes of dN/dS ratio by free-ratio model analysis for subterranean Arvicolinae. Parallel substitutions were detected in genes RAD23B and PYCR2. These genes are associated with DNA repair processes and response to oxidative stress. Similar substitutions were discovered in the RAD23 genes for highly specialized subterranean Heterocephalus glaber and Fukomys damarensis. The most pronounced signatures of adaptive evolution related to subterranean niche within species of Arvicolinae subfamily were detected for Ellobius lutescens. Our results suggest that genomic adaptations can occur very quickly so far as the amount of selection signatures was found to be compliant with the degree of specialization to the subterranean niche and independent from the evolutionary age of the taxon. We found that the number of genomic signatures of selection does not depend on the age of the taxon, but is positively correlated with the degree of specialization to the subterranean niche.

Enlarged fins of Tibetan catfish provide new evidence of adaptation to high plateau

The uplift of the Tibetan Plateau significantly altered the geomorphology and climate of the Euroasia by creating large mountains and rivers. Fishes are more likely to be affected relative to other organisms, as they are largely restricted to river systems. Faced with the rapidly flowing water in the Tibetan Plateau, a group of catfish has evolved greatly enlarged pectoral fins with more numbers of fin-rays to form an adhesive apparatus. However, the genetic basis of these adaptations in Tibetan catfishes remains elusive. In this study, we performed comparative genomic analyses based on the chromosome-level genome of Glyptosternum maculatum in family Sisoridae and detected some proteins with conspicuously high evolutionary rates in particular in genes involved in skeleton development, energy metabolism, and hypoxia response. We found that the hoxd12a gene evolved faster and a loss-of-function assay of hoxd12a supports a potential role for this gene in shaping the enlarged fins of these Tibetan catfishes. Other genes with amino acid replacements and signatures of positive selection included proteins involved in low temperature (TRMU) and hypoxia (VHL) responses. Functional assays reveal that the G. maculatumTRMU allele generates more mitochondrial ATP than the ancestral allele found in low-altitude fishes. Functional assays of VHL alleles suggest that the G. maculatum allele has lower transactivation activity than the low-altitude forms. These findings provide a window into the genomic underpinnings of physiological adaptations that permit G. maculatum to survive in the harsh environment of the Tibetan Himalayas that mirror those that are convergently found in other vertebrates such as humans.

Comparative mitogenomic and evolutionary analysis of Lycaenidae (Insecta: Lepidoptera): Potential association with high-altitude adaptation

Harsh environments (e.g., hypoxia and cold temperatures) of the Qinghai-Tibetan Plateau have a substantial influence on adaptive evolution in various species. Some species in Lycaenidae, a large and widely distributed family of butterflies, are adapted to the Qinghai-Tibetan Plateau. Here, we sequenced four mitogenomes of two lycaenid species in the Qinghai-Tibetan Plateau and performed a detailed comparative mitogenomic analysis including nine other lycaenid mitogenomes (nine species) to explore the molecular basis of high-altitude adaptation. Based on mitogenomic data, Bayesian inference, and maximum likelihood methods, we recovered a lycaenid phylogeny of [Curetinae + (Aphnaeinae + (Lycaeninae + (Theclinae + Polyommatinae)))]. The gene content, gene arrangement, base composition, codon usage, and transfer RNA genes (sequence and structure) were highly conserved within Lycaenidae. TrnS1 not only lacked the dihydrouridine arm but also showed anticodon and copy number diversity. The ratios of non-synonymous substitutions to synonymous substitutions of 13 protein-coding genes (PCGs) were less than 1.0, indicating that all PCGs evolved under purifying selection. However, signals of positive selection were detected in cox1 in the two Qinghai-Tibetan Plateau lycaenid species, indicating that this gene may be associated with high-altitude adaptation. Three large non-coding regions, i.e., rrnS-trnM (control region), trnQ-nad2, and trnS2-nad1, were found in the mitogenomes of all lycaenid species. Conserved motifs in three non-coding regions (trnE-trnF, trnS1-trnE, and trnP-nad6) and long sequences in two non-coding regions (nad6-cob and cob-trnS2) were detected in the Qinghai-Tibetan Plateau lycaenid species, suggesting that these non-coding regions were involved in high-altitude adaptation. In addition to the characterization of Lycaenidae mitogenomes, this study highlights the importance of both PCGs and non-coding regions in high-altitude adaptation.

Triplophysadaryoae, a new nemacheilid loach species (Teleostei, Nemacheilidae) from the Syr Darya River basin, Central Asia

Triplophysadaryoae, new species, is described from the Sokh River, a former tributary of Syr Darya that today fails to reach the river, in the Sokh District, an exclave of Uzbekistan, surrounded by Kyrgyzstan. Triplophysadaryoae is distinguished from other species of Triplophysa in Central Asia by a truncate caudal fin with 13 or 14 branched rays, body without obvious mottling, dorsal-fin origin opposite to pelvic-fin insertion, and absence of the posterior chamber of the air bladder. Molecular data suggest that Triplophysadaryoae is closely related to T.ferganaensis from the Shakhimardan stream, a small tributary of Syr Darya in the Yordon village, another exclave of Uzbekistan in Kyrgyzstan. The two species were separated by a Kimura 2-parameter genetic distance of 2.8% in the mitochondrial DNA cytochrome c oxidase subunit I barcode region; they are also distinguished morphologically. A key to the species of Triplophysa in the Syr Darya basin and adjacent regions is provided.

Comparative Mitogenomics of Two Sympatric Catfishes of Exostoma (Siluriformes: Sisoridae) from the Lower Yarlung Tsangpo River and Its Application for Phylogenetic Consideration

The genus Exostoma is a group of stenotopic and rheophilic glyptosternine catfishes distributed in South and Southeast Asia. So far, comprehensive studies on mitogenomics referring to this genus are very scarce. In this study, we first sequenced and annotated the complete mitochondrial genomes of Exostoma tibetanum and Exostoma tenuicaudatum—two sympatric congeners from the lower Yarlung Tsangpo River, Tibet, China. The mitogenomes of both species contained 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNA genes, one light-strand origin of replication, and one control region, with lengths of 16,528 bp and 16,533 bp, respectively. The mitogenome architecture, nucleotide composition, and codon usage of protein-coding genes were almost identical between the two Exostoma species, although some estimated parameters varied. Phylogenetic analysis strongly supported the monophyly of Exostoma in the subfamily Glyptosternae, and Exostoma tibetanum had the closest relationship to Exostoma tenuicaudatum. The divergence time estimation demonstrated that these two species diverged approximately 1.51 Ma during the early Pleistocene, which was speculated to be triggered by the river system changes caused by the uplift of the southeastern Tibetan Plateau. Selection pressure analyses indicated that all protein-coding genes of Exostoma species underwent a strong purifying selection, while minority positive sites from NADH dehydrogenase complex genes were detected. These findings are expected to promote our understanding of the molecular phylogeny of the genus Exostoma and provide valuable mitogenomic resources for the subfamily Glyptosternae

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.

Adaptability and Evolution of Gobiidae: A Genetic Exploration

The Gobiidae family occupy one of the most diverse habitat ranges of all fishes. One key reason for their successful colonization of different habitats is their ability to adapt to different energy demands. This energy requirement is related to the ability of mitochondria in cells to generate energy via oxidative phosphorylation (OXPHOS). Here, we assembled three complete mitochondrial genomes of Rhinogobius shennongensis, Rhinogobius wuyanlingensis, and Chaenogobius annularis. These mitogenomes are circular and include 13 protein-coding genes (PCGs), two rRNAs, 22 tRNAs, and one non-coding control region (CR). We used comparative mitochondrial DNA (mtDNA) genome and selection pressure analyses to explore the structure and evolutionary rates of Gobiidae mitogenomics in different environments. The CmC model showed that the ω ratios of all mtDNA PCGs were <1, and that the evolutionary rate of adenosine triphosphate 8 (atp8) was faster in Gobiidae than in other mitochondrial DNA PCGs. We also found evidence of positive selection for several sites of NADH dehydrogenase (nd) 6 and atp8 genes. Thus, divergent mechanisms appear to underlie the evolution of mtDNA PCGs, which might explain the ability of Gobiidae to adapt to diverse environments. Our study provides new insights on the adaptive evolution of Gobiidae mtDNA genome and molecular mechanisms of OXPHOS.

Comparative Analysis of the Mitochondrial Genome of Galatheanthemum sp. MT-2020 (Actiniaria Galatheanthemidae) From a Depth of 9,462 m at the Mariana Trench

The hadal zone, which represents the deepest marine habitat on Earth (6,000–11,000 m), is a harsh environment mainly characterized by extremely high hydrostatic pressure, and this habitat is believed to have a high degree of endemism. The deep-sea anemone family Galatheanthemidae comprises two valid species exclusively from the hadal; however, no other information about this family is currently available. In the present study, a sea anemone was collected from a depth of 9,462 m at the Mariana Trench and was defined as Galatheanthemum sp. MT-2020 (Actiniaria Galatheanthemidae). The mitochondrial genome of Galatheanthemum sp. MT-2020 was circular, was 16,633 bp in length, and contained two ribosomal RNA genes, 13 protein-coding genes and two transfer RNA genes. The order of the genes of Galatheanthemum sp. MT-2020 was identical to that of the majority of the species of the order Actiniaria. The value of the AT-skew was the lowest in the whole mitochondrial genome, with a positive GC skew value for the atp8 gene, while other species, except Antholoba achates, had the negative values of the GC skew. Galatheanthemum sp. MT-2020 was clustered with another abyssal species, Paraphelliactis xishaensis, in the phylogenetic tree, and these species diverged in the early Jurassic approximately 200 Mya from the shallow-sea species. The usage ratio of valine, which is one of the five amino acids with the strongest barophilic properties, in the mitochondrial genomes of the two abyssal species was significantly higher than that in other species with habitats above the depth of 3,000 m. The ω (dN/dS) ratio of the genomes was 2.45-fold higher than that of the shallow-sea species, indicating a slower evolutionary rate. Overall, the present study is the first to provide a complete mitogenome of sea anemones from the hadal and reveal some characteristics that may be associated with adaptation to an extreme environment.

Effect of secondary attack by Aeromonas hydrophila on the expression level of hif genes in common carp (Cyprinus carpio)

Hypoxia-inducible factors (hifs) are involved in infectious diseases inflammatory reactions, and immune regulation. Common carp, a representative allotetraploid species that has undergone genome-wide replication events, has important research value. In this study, common carp were infected twice with Aeromonas hydrophila. Liver tissues of common carp were collected at 4 h, 12 h, 24 h, 48 h, 3 days, 7 days post-first infection and 4 h, 12 h, 24 h post-second infection. The mRNA levels of hif genes were determined at different time points. The hif2a-2, hif3a-2, hif3b-1 and hif3b-2 expression levels in the infected group were upregulated when compared with those in the control group, whereas the expression levels of other genes were downregulated after the second infection. This indicates that the effect of A. hydrophila infection on gene expression pattern is dependent on the host, pathogen, infected tissue and gene. Pressure analysis of the hif gene family revealed that the non-synonymous substitution to synonymous substitution ratio of 12 hif genes was <1, which indicated that they were in a state of purification and selection. Combined with the differences between copy genes, the polyclonal antibodies against Hif1b-1 and Hif1b-2 were successfully prepared in this study. Western blot analysis showed that the protein expression of Hif1b-1 and Hif1b-2 reached to the highest level 48 h after the first infection. After the second A. hydrophila infection, the protein expression levels of Hif1b-1 and Hif1b-2 reached the highest levels at 4 and 48 h, respectively. This may indicate that the Hif1b-1 and Hif1b-2 genes in common carp play an important role in the immune mechanism at the protein level. The findings of this study will lay the foundation for future studies on the immune regulatory function of common carp hif genes, which may aid in devising novel therapeutic strategies for common carp diseases, such as A. hydrophila infection.

Phylogeographic relationships and the evolutionary history of the Carassius auratus complex with a newly born homodiploid raw fish (2nNCRC)

Background

An important aspect of studying evolution is to understand how new species are formed and their uniqueness is maintained. Hybridization can lead to the formation of new species through reorganization of the adaptive system and significant changes in phenotype. Interestingly, eight stable strains of 2nNCRC derived from interspecies hybridization have been established in our laboratory. To examine the phylogeographical pattern of the widely distributed genus Carassius across Eurasia and investigate the possible homoploid hybrid origin of the Carassius auratus complex lineage in light of past climatic events, the mitochondrial genome (mtDNA) and one nuclear DNA were used to reconstruct the phylogenetic relationship between the C. auratus complex and 2nNCRC and to assess how demographic history, dispersal and barriers to gene flow have led to the current distribution of the C. auratus complex.

Results

As expected, 2nNCRC had a very close relationship with the C. auratus complex and similar morphological characteristics to those of the C. auratus complex, which is genetically distinct from the other three species of Carassius. The estimation of divergence time and ancestral state demonstrated that the C. auratus complex possibly originated from the Yangtze River basin in China. There were seven sublineages of the C. auratus complex across Eurasia and at least four mtDNA lineages endemic to particular geographical regions in China. The primary colonization route from China to Mongolia and the Far East (Russia) occurred during the Late Pliocene, and the diversification of other sublineages of the C. auratus complex specifically coincided with the interglacial stage during the Early and Mid-Pleistocene in China.

Conclusion

Our results support the origin of the C. auratus complex in China, and its wide distribution across Eurasia was mainly due to natural Pleistocene dispersal and recent anthropogenic translocation. The sympatric distribution of the ancestral area for both parents of 2nNCRC and the C. auratus complex, as well as the significant changes in the structure of pharyngeal teeth and morphological characteristics between 2nNCRC and its parents, imply that homoploid hybrid speciation (HHS) for C. auratus could likely have occurred in nature. The diversification pattern indicated an independent evolutionary history of the C. auratus complex, which was not separated from the most recent common ancestor of C. carassius or C. cuvieri. Considering that the paleoclimate oscillation and the development of an eastward-flowing drainage system during the Pliocene and Pleistocene in China provided an opportunity for hybridization between divergent lineages, the formation of 2nNCRC in our laboratory could be a good candidate for explaining the HHS of C. auratus in nature.

Mitogenome of Knodus borki (Cypriniformes: Characidae): genomic characterization and phylogenetic analysis

BACKGROUND

The taxonomic status of Knodu in the family Characidae is not yet clear. This study aimed to address this by sequencing and annotating Knodu borki Zarske, 2008.

MATERIALS AND RESULTS

K. borki Zarske, 2008 was sequenced using a Hiseq platform and the complete mitogenome was assembled in SPAdes v3.15.2 and SOAPdenovo2 v.2.01. The mitogenome of K. borki from Guangzhou, the first sequenced species of the genus Knodu, is 16,837 bp in length and contains 13 protein-coding genes (PCGs), two ribosomal (r) RNAs, 22 transfer (t) RNAs, and one D-loop. Among these 37 genes, 28 are encoded by the heavy strand, while nine are encoded by the light strand. Twenty-one of the tRNAs can form typical cloverleaf secondary structures, except tRNA-Ser1, which lacks dihydrouridine arms. All PCGs have the same start codon (ATG), with the exception of COI (GTG). Four PCGs (ND1, ATP8, ND4L, and ND5) have TAA as the stop codon, ND6 has TAG as the stop codon, COI has AGG as the stop codon, and the remaining seven genes have incomplete stop codons of TA-/T-(ND2, COII, COIII, ND3, ND4, and Cyt b as T-, ATP6 as TA-). Phylogenetic analysis showed that K. borki belongs to the family Characidae.

CONCLUSIONS

Our findings demonstrate that K. borki belongs to the family Characidae, due to consistency with the morphological identification. This study provides molecular information for further research on the phylogeny of the genus Knodus and for analyses of the taxonomic status of Characidae.

Signatures of Adaptation in Mitochondrial Genomes of Palearctic Subterranean Voles (Arvicolinae, Rodentia)

This study evaluates signatures of selection in the evolution of the mitochondrial DNA of voles, subfamily Arvicolinae, during the colonization of subterranean environments. The comparative sequence analysis of mitochondrial protein-coding genes of eight subterranean vole species (Prometheomys schaposchnikowi, three species of the genus Ellobius: Ellobius talpinus, Ellobius fuscocapillus and Ellobius lutescens, two species of the genus Terricola: Terricola subterraneus and Terricola daghestanicus, Lasiopodomys mandarinus, and Hyperacrius fertilis) and their closest aboveground relatives was applied using codon-substitution models. The highest number of selection signatures was detected in genes ATP8 and CYTB. The relaxation of selection was observed in most mitochondrial DNA protein-coding genes for subterranean species. The largest amount of relaxed genes is discovered in mole voles (genus Ellobius). The number of selection signatures was found to be independent of the evolutionary age of the lineage but fits the degree of specialization to the subterranean niche. The common trends of selective pressures were observed among the evolutionary ancient and highly specialized subterranean rodent families and phylogenetically young lineages of voles. It suggests that the signatures of adaptation in individual mitochondrial protein-coding genes associated with the colonization of the subterranean niche may appear within a rather short evolutionary timespan.

Complete mitochondrial genomes reveal robust phylogenetic signals and evidence of positive selection in horseshoe bats

BACKGROUND

In genus Rhinolophus, species in the Rhinolophus philippinensis and R. macrotis groups are unique because the horseshoe bats in these group have relatively low echolocation frequencies and flight speeds compared with other horseshoe bats with similar body size. The different characteristics among bat species suggest particular evolutionary processes may have occurred in this genus. To study the adaptive evidence in the mitochondrial genomes (mitogenomes) of rhinolophids, especially the mitogenomes of the species with low echolocation frequencies, we sequenced eight mitogenomes and used them for comparative studies of molecular phylogeny and adaptive evolution.

RESULTS

Phylogenetic analysis using whole mitogenome sequences produced robust results and provided phylogenetic signals that were better than those obtained using single genes. The results supported the recent establishment of the separate macrotis group. The signals of adaptive evolution discovered in the Rhinolophus species were tested for some of the codons in two genes (ND2 and ND6) that encode NADH dehydrogenases in oxidative phosphorylation system complex I. These genes have a background of widespread purifying selection. Signals of relaxed purifying selection and positive selection were found in ND2 and ND6, respectively, based on codon models and physicochemical profiles of amino acid replacements. However, no pronounced overlap was found for non-synonymous sites in the mitogenomes of all the species with low echolocation frequencies. A signal of positive selection for ND5 was found in the branch-site model when R. philippinensis was set as the foreground branch.

CONCLUSIONS

The mitogenomes provided robust phylogenetic signals that were much more informative than the signals obtained using single mitochondrial genes. Two mitochondrial genes that encoding proteins in the oxidative phosphorylation system showed some evidence of adaptive evolution in genus Rhinolophus and the positive selection signals were tested for ND5 in R. philippinensis. These results indicate that mitochondrial protein-coding genes were targets of adaptive evolution during the evolution of Rhinolophus species, which might have contributed to a diverse range of acoustic adaptations in this genus.

High-altitude adaptation in vertebrates as revealed by mitochondrial genome analyses

The high-altitude environment may drive vertebrate evolution in a certain way, and vertebrates living in different altitude environments might have different energy requirements. We hypothesized that the high-altitude environment might impose different influences on vertebrate mitochondrial genomes (mtDNA). We used selection pressure analyses and PIC (phylogenetic independent contrasts) analysis to detect the evolutionary rate of vertebrate mtDNA protein-coding genes (PCGs) from different altitudes. The results showed that the ratio of nonsynonymous/synonymous substitutions (dN/dS) in the mtDNA PCGs was significantly higher in high-altitude vertebrates than in low-altitude vertebrates. The seven rapidly evolving genes were shared by the high-altitude vertebrates, and only one positive selection gene (ND5 gene) was detected in the high-altitude vertebrates. Our results suggest the mtDNA evolutionary rate in high-altitude vertebrates was higher than in low-altitude vertebrates as their evolution requires more energy in a high-altitude environment. Our study demonstrates the high-altitude environment (low atmospheric O2 levels) drives vertebrate evolution in mtDNA PCGs.

Migration effects on the intestinal microbiota of Tibetans

BACKGROUND

Diet, environment, and genomic context have a significant impact on humans' intestinal microbiota. Moreover, migration may be accompanied by changes in human eating habits and living environment, which could, in turn, affect the intestinal microbiota. Located in southwestern China, Tibet has an average altitude of 4,000 meters and is known as the world's roof. Xianyang is situated in the plains of central China, with an average altitude of about 400 meters.

METHODS

To understand the association between intestinal microbiota and population migration, we collected the fecal samples from 30 Tibetan women on the first day (as TI1st), six months (as TI2nd), and ten months (as TI3rd) following migration from Tibet to Xianyang. Fecal samples were collected from 29 individuals (belonging to the Han women) as a control. The dietary information of the Tibetan women and the Han women was gathered. We performed a 16S rRNA gene survey of the collected fecal samples using Illumina MiSeq sequencing.

RESULTS

Following the migration, the alpha and beta diversity of Tibetan women's intestinal microbiota appeared unaffected. Linear discriminant analysis effect size (LEfSe) analysis showed that Klebsiella, Blautia, and Veillonella are potential biomarkers at TI1st, while Proteobacteria and Enterobacteriaceae were common in TI3rd. Finally, functional prediction by phylogenetic investigation of communities by reconstruction of unobserved states (PICRUSt) found no significant up-regulation or down-regulation gene pathway in the intestinal microbiota of Tibetan women after migration. The present study reveals that the higher stability in Tibetan women's intestinal microbiota was less affected by the environment and diet, indicating that Tibetan women's intestinal microbiota is relatively stable. The main limitations of the study were the small sample size and all volunteers were women.

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.

Triplophysa ferganaensis, a new loach species from Fergana Valley in Central Asia (Teleostei: Nemacheilidae)

Triplophysa ferganaensis sp. nov. is described from the Shakhimardan stream, a small tributary of the Syr Darya, which does not reach the river in Fergana Valley. It can be distinguished from other valid Triplophysa loaches based on the following combination of characters: body smooth and scaleless, lateral line complete, posterior chamber of air bladder degenerated, inner gill rakers 10-11 on the first-gill arch, outer gill rakers absent, vertebrae 4 + 35-36, 8 + 8 branched caudal-fin rays, caudal peduncle depth 2.1-2.7 times its length, two supratemporal pores, dorsal-fin origin closer to the caudal-fin base than to the snout tip, caudal fin emarginated and pelvic-fin tip reaching the anus. The new species can also be distinguished from its congeners based on the molecular analyses of mitochondrial cytochrome oxidase subunit I (coI) gene sequences. The phylogenetic position of this new species indicates that it is a sister taxon of Triplophysa tenuis.

Identification and expression analysis of lncRNA in seven organs of Rhinopithecus roxellana

Long non-coding RNA (lncRNA) represents a new direction to identify expression profiles and regulatory mechanisms in various organisms. Here, we report the first dataset of lncRNAs of the golden snub-nosed monkey (GSM), including 12,557 putative lncRNAs identified from seven organs. Compared with mRNA, GSM lncRNA had fewer exons and isoforms, and longer length. LncRNA showed more obvious tissue-specific expression than mRNA. However, for the top ten most abundant genes in each organ, mRNAs expression was more tissue-specific than lncRNAs. By identification of specifically expressed lncRNAs and mRNAs in each organ, it indicates that the expression of SEG-lncRNA (specifically expressed lncRNA) and SEG-mRNA (specifically expressed mRNA) had high correlation. In particular, combined our lncRNA and mRNA data, we identified 92 heart SEG-lncRNAs targeted ten mRNA genes in the oxidative phosphorylation pathway and upregulated the expression of these target genes such as ND4, ATP6, and ATP8. These may contribute to GSM adaption to its high-elevation environment. We also identified 171 liver SEG-lncRNAs, which targeted 27 genes associated with the metabolism of xenobiotics and leaded to high expression of these target genes in liver. These lncRNAs may play important roles in GSM adaptation to a folivory diet.

The complete mitogenome of Phymorhynchus sp. (Neogastropoda, Conoidea, Raphitomidae) provides insights into the deep-sea adaptive evolution of Conoidea

The deep-sea environment is characterized by darkness, hypoxia, and high hydrostatic pressure. Mitochondria play a vital role in energy metabolism; thus, they may endure the selection process during the adaptive evolution of deep-sea organisms. In the present study, the mitogenome of Phymorhynchus sp. from the Haima methane seep was completely assembled and characterized. This mitogenome is 16,681 bp in length and contains 13 protein-coding genes, 2 rRNAs, and 22 tRNAs. The gene order and orientation were identical to those of most sequenced conoidean gastropods. Some special elements, such as tandem repeat sequences and AT-rich sequences, which are involved in the regulation of the replication and transcription of the mitogenome, were observed in the control region. Phylogenetic analysis revealed that Conoidea is divided into two separate clades with high nodal support. Positive selection analysis revealed evidence of adaptive changes in the mitogenomes of deep-sea conoidean gastropods. Eight residues located in atp6, cox1, cytb, nad1, nad4, and nad5 were determined to have undergone positive selection. This study explores the adaptive evolution of deep-sea conoidean gastropods and provides valuable clues at the mitochondrial level regarding the exceptional adaptive ability of organisms in deep-sea environments.

Triplophysa wulongensis, a new species of cave-dwelling loach (Teleostei, Nemacheilidae) from Chongqing, Southwest China

We describe a new species of cave-dwelling loach, Triplophysawulongensissp. nov., based on specimens collected in a subterranean pool in a cave in Wulong County, Chongqing, Southwest China. The pool is connected to the Wujiang River drainage. Triplophysawulongensis differs from its congeners by the following combination of characters: eyes present, caudal fin with 18 branched rays; posterior chamber of the air bladder degenerate; stomach U-shaped; intestine without bends or loops immediately posterior to stomach; body smooth and scaleless, and lateral line complete. The mitochondrial cytochrome b sequence differs from those of other published sequences of species of Triplophysa by 14.9–24.9% in K2P distance. Phylogenetic analysis based on cytochrome b gene sequences recovered T.wulongensis as sister taxon to all other cave-dwelling species of Triplophysa.

The complete mitochondrial genome of a cave-dwelling loach Triplophysa baotianensis (Teleostei: Nemacheilidae)

Triplophysa baotianensis belong to the genus Triplophysa (Teleostei, Nemacheilidae), endemic to Guizhou Province, Southwestern China. In this study, the complete mitochondrial genome of T. baotianensis was sequenced and reported for the first time. The circular mitogenome was 16,576 bp in length and consisted of 13 protein-coding genes, 2 ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and 1 non-coding control region. The overall base composition was 30.79% A, 27.62% T, 25.46% C, and 16.13% G with 41.59% GC content. Phylogenetic analysis using mitochondrial genomes of 40 species showed that all Triplophysa species clustered as one monophyletic clade, and T. baotianensis was the closest to (T. nasobarbatula + (T. rosa + T. xiangxiensis)).

Island colonization by a ‘rheophilic’ fish: the phylogeography of Garra ceylonensis (Teleostei: Cyprinidae) in Sri Lanka

Despite exhibiting multiple morphological adaptations to living in swiftly flowing water (rheophily), Garra ceylonensis is one of the most widely distributed freshwater fish in Sri Lanka. It is thus an ideal organism to reconstruct the evolutionary history of a widespread, yet morphologically specialized, freshwater fish in a tropical-island setting. We analysed the phylogenetic and phylogeographic relationships of G. ceylonensis based on two mitochondrial and one nuclear genes. G. ceylonensis is shown to be monophyletic, with a sister-group relationship to the Indian species Garra mullya. Our results suggest a single colonization of Sri Lanka by ancestral Garra, in the late Pliocene. This suggests that the Palk Isthmus, which was exposed for most of the Pleistocene, had a hydroclimate unsuited to the dispersal of fishes such as Garra. G. ceylonensis exhibits strong phylogeographic structure: six subclades are distributed as genetically distinct populations in clusters of contiguous river basins, albeit with two exceptions. Our data reveal one or more Pleistocene extirpation events, evidently driven by aridification, with relict populations subsequently re-colonizing the island. The phylogeographic structure of G. ceylonensis suggests inter-basin dispersal largely through headwater capture, likely facilitated by free-swimming post-larvae. The Peninsular-Indian species G. mullya comprises two genetically distinct parapatric clades, which may represent distinct species.

Variations in the Mitochondrial Genome of a Goldfish-Like Hybrid [Koi Carp (♀) × Blunt Snout Bream (♂)] Indicate Paternal Leakage

Previously, a homodiploid goldfish-like fish (2n = 100; GF-L) was spontaneously generated by self-crossing a homodiploid red crucian carp-like fish (2n = 100; RCC-L), which was in turn produced via the distant hybridization of female koi carp (Cyprinus carpio haematopterus, KOC, 2n = 100) and male blunt snout bream (Megalobrama amblycephala, BSB, 2n = 48). The phenotypes and genotypes of RCC-L and GF-L differed from those of the parental species but were similar to diploid red crucian carp (2n = 100; RCC) and goldfish (2n = 100; GF), respectively. We sequenced the complete mitochondrial DNAs (mtDNAs) of the KOC, BSB, RCC-L, GF-L, and subsequent generations produced by self-crossing [the self-mating offspring of RCC-L (RCC-L-F₂) to the self-mating offspring of RCC-L-F₂ (RCC-L-F₃) and the self-mating offspring of GF-L (GF-L-F₂)]. Paternal mtDNA fragments were stably embedded in the mtDNAs of both lineages, forming chimeric DNA fragments. In addition to these chimeras, several nucleotide positions in the RCC-L and GF-L lineages differed from the parental bases, and were instead identical with RCC and GF, respectively. Moreover, RCC-L and GF-L mtDNA organization and nucleotide composition were more similar to those of RCC and GF, respectively, compared to parental mtDNA. Finally, phylogenetic analyses indicated that RCC-L and GF-L clustered with RCC and GF, not with the parental species. The molecular dating time shows that the divergence time of KOC and GF was about 21.26 Mya [95% highest posterior density (HPD): 24.41-16.67 Mya], which fell within the period of recent. The heritable chimeric DNA fragments and mutant loci identified in the mtDNA of the RCC-L and GF-L lineages provided important evidence that hybridizations might lead to changes in the mtDNA and the subsequent generation of new lineages. Our findings also demonstrated for the first time that the paternal mtDNA was transmitted into the mtDNA of homodiploid lineages (RCC-L and GF-L), which provided evidence that paternal DNA plays a role in inherited mtDNA. These evolutionary analyses in mtDNA suggest that GF might have diverged from RCC after RCC diverged from koi carp.