Guogang Li

Local adaptation of Gymnocypris przewalskii (Cyprinidae) on the Tibetan Plateau

Divergent selection among environments affects species distributions and can lead to speciation. In this article, we investigated the transcriptomes of two ecotypes of scaleless carp (Gymnocypris przewalskii przewalskii and G. p. ganzihonensis) from the Tibetan Plateau. We used a transcriptome sequencing approach to screen approximately 250,000 expressed sequence tags (ESTs) from the gill and kidney tissues of twelve individuals from the Ganzi River and Lake Qinghai to understand how this freshwater fish has adapted to an ecological niche shift from saline to freshwater. We identified 9,429 loci in the gill transcriptome and 12,034 loci in the kidney transcriptome with significant differences in their expression, of which 242 protein-coding genes exhibited strong positive selection (Ka/Ks > 1). Many of the genes are involved in ion channel functions (e.g., Ca2+-binding proteins), immune responses (e.g., nephrosin) or cellular water absorption functions (e.g., aquaporins). These results have potentially broad importance in understanding shifts from saline to freshwater habitats. Furthermore, this study provides the first transcriptome of G. przewalskii, which will facilitate future ecological genomics studies and aid in the identification of genes underlying adaptation and incipient ecological speciation.

Ongoing Speciation in the Tibetan Plateau Gymnocypris Species Complex

Local adaptation towards divergent ecological conditions often results in genetic differentiation and adaptive phenotypic divergence. To illuminate the ecological distinctiveness of the schizothoracine fish, we studied a Gymnocypris species complex consisting of three morphs distributed across four bodies of water (the Yellow River, Lake Qinghai, the Ganzi River and Lake Keluke) in the Northeast Tibetan Plateau. We used a combination of mitochondrial (16S rRNA and Cyt b) and nuclear (RAG-2) genetic sequences to investigate the phylogeography of these morphs based on a sample of 277 specimens. Analysis of gill rakers allowed for mapping of phenotypic trajectories along the phylogeny. The phylogenetic and morphological analyses showed that the three sparsely rakered morphs were present at two extremes of the phylogenetic tree: the Yellow River morphs were located at the basal phylogenetic split, and the Lake Keluke and Ganzi River morphs at the peak, with the densely rakered Lake Qinghai morphs located between these two extremes. Age estimation further indicated that the sparsely rakered morphs constituted the oldest and youngest lineages, whereas the densely rakered morph was assigned to an intermediate-age lineage. These results are most compatible with the process of evolutionary convergence or reversal. Disruptive natural selection due to divergent habitats and dietary preferences is likely the driving force behind the formation of new morphs, and the similarities between their phenotypes may be attributable to the similarities between their forms of niche tracking associated with food acquisition. This study provides the first genetic evidence for the occurrence of convergence or reversal in the schizothoracine fish of the Tibetan Plateau at small temporal scales.

Characterization of two paralogous myostatin genes and evidence for positive selection in Tibet fish: Gymnocypris przewalskii

Myostatin (mstn) is an important member of TGF-β superfamily, a muscle growth inhibitor. Though mstn has been identified in many organisms, little is known about this gene in highland fish, Gymnocypris przewalskii endemic to the Qinghai-Tibetan Plateau. In this study, we first cloned two paralogous mstn genes (mstn1 and mstn2) from G. przewalskii through homologue cloning. The 3D structures of both Mstn proteins varied in the numbers of β-sheets and conformations of α-helices. The branch-site model showed that mstn1 has undergone positive selection, and two positively selected sites (107M and 181T) were located on the random coils of the 3D protein structure. Expression patterns indicated that the mstn1 expressed widely, while the mstn2 only expressed in the muscle and brain. During the early stage of embryo development, the expression levels of both mstn paralogous genes showed different increasing trends. These results suggest that it is diverging in two mstn paralogues of G. przewalskii via specific differences in gene structure, protein structure, selection pressure and gene expression patterns. Taken together, this study provides novel contribution on the research topics of growth related gene function and mechanism of highland fish in extreme aquatic environment on the Qinghai-Tibetan Plateau.

Taxonomic status of a population of Gymoncypris waddelli Regan, 1905 (Cypriniformes: Schizothoracinae) distributed in Pengqu River, Tibet, China

Gymnocypris waddelli, a highland cold-water fish, is distributed among rivers and lakes of the southern Qinghai-Tibet Plateau. Although previously described as G. waddelli, specimens from the Pengqu River are morphologically distinguishable from those from Lake Yamzhog Yumco, which is the type locality of G. waddelli. We investigated morphological variations of G. waddelli specimens from Pengqu River and Lake Yamzhog Yumco using not only traditional morphological methods but also newly-developed geometric morphometric analysis; furthermore, mitochondrial Cytochrome b (Cyt b) gene analysis was employed to explore the phylogenetic position of the Pengqu River population. Our morphological analyses suggested that G. waddelli from the Pengqu River had an obtuse snout, nearly straight oral fissure, shorter lower jaw than upper jaw, broad mouth and sparse inside gill rakers of the first arch relative to the Lake Yamzhog Yumco population. The strongly supported monophyly of Pengqu population according to mitochondrial datasets was obtained, indicating genetic differences between Pengqu River and Lake Yamzhog Yumco populations. Given their substantial genetic and morphological divergence, specimens from Pengqu River warrant recognition as a distinct species, here described and named Gymnocypris pengquensis.

Phylogeography of Diptychus maculatus (Cyprinidae) endemic to the northern margin of the QTP and Tien Shan region

BackgroundPhylogeography and historical demography of the cyprinid fish Diptychus maculatus (subfamily Schizothoracinae) are evaluated across three river systems in the Northern Qinghai-Tibetan Plateau (QTP) and Tien Shan range: the Indus River, Tarim River and Ili River.ResultsResults from both mtDNA (16S rRNA, Cyt b and D-loop) and nucDNA (RAG-2) resolved four reciprocally monophyletic clades, representing populations from Indus River, South Tarim River, North Tarim River and Ili River, respectively. The divergence times was estimated to be 1.5–2.5 Mya. It is consistent with the hypothesis that the split of four clades is the consequence of vicariance resulting from both the intensive uplift of QTP and Tien Shan as well as the resultant expansion of the Taklimakan Desert. Several lines of evidences indicate dynamic demographic histories for the populations, with late Pleistocene and Holocene population bottlenecks and expansions except the Indus River.ConclusionsOur results clearly depicted the phylogenetic relationship of D. maculatus from Indus River, Tarim River and Ili River. The analyses implicated the relationship among the distribution of D. maculatus, paleo-drainages and geographic events, and implied the existence of the South Tarim River in history.

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

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

Phylogeography of Schizopygopsis stoliczkai (Cyprinidae) in Northwest Tibetan Plateau area

Abstract Schizopygopsis stoliczkai (Cyprinidae, subfamily Schizothoracinae) is one of the major freshwater fishes endemic to the northwestern margin of the Tibetan Plateau. In the current study, we used mitochondrial DNA markers cytochrome b (Cyt b) and 16S rRNA (16S), as well as the nuclear marker, the second intron of the nuclear beta‐actin gene (Act2), to uncover the phylogeography of S. stoliczkai. In total, we obtained 74 haplotypes from 403 mitochondrial concatenated sequences. The mtDNA markers depict the phylogenetic structures of S. stoliczkai, which consist of clade North and clade South. The split time of the two clades is dated back to 4.27 Mya (95% HPD = 1.96–8.20 Mya). The estimated split time is earlier than the beginning of the ice age of Pleistocene (2.60 Mya), suggesting that the northwestern area of the Tibetan Plateau probably contain at least two glacial refugia for S. stoliczkai. SAMOVA supports the formation of four groups: (i) the Karakash River group; (ii) The Lake Pangong group; (iii) the Shiquan River group; (iv) the Southern Basin group. Clade North included Karakash River, Lake Pangong, and Shiquan River groups, while seven populations of clade South share the haplotypes. Genetic diversity, star‐like network, BSP analysis, as well as negative neutrality tests indicate recent expansions events of S. stoliczkai. Conclusively, our results illustrate the phylogeography of S. stoliczkai, implying the Shiquan River is presumably the main refuge for S. stoliczkai.

Elevational diversity gradients of Tibetan loaches: The relative roles of ecological and evolutionary processes

Abstract It is widely believed that species richness patterns (SRPs) are shaped by both ecological and evolutionary processes. However, the relative roles of these processes remain unclear, especially for aquatic organisms. In this study, we integrated ecological and evolutionary measures to tease apart the relative influences of these factors on the SRP of Tibetan loaches along an extensive elevational gradient. We found that the Tibetan loaches displayed a richness pattern that peaked at midelevations. The mean annual temperature (MAT), mid‐domain effect (MDE), and summed age of colonization (SAC, complex of colonization age and colonization frequency) were the main drivers, accounting for 85%, 51%, and 88% of the variations in the SRP, respectively. The three predictors had very high combined effects (MAT‐MDE‐SAC, MAT‐SAC, and MDE‐SAC were 44%, 38%, and 6%, respectively). Our analyses suggested that energy input, time‐for‐speciation, and species dispersal may directly guide the SRP or mediate it by geometric constraints. Conclusively, the SRP of the Tibetan loaches with elevation is the outcome of interactions between biogeographical processes and regional ecological conditions.