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

The Himalayan uplift and evolution of aquatic biodiversity across Asia: Snowtrout (Cyprininae: Schizothorax) as a test case

Global biodiversity hotspots are often remote, tectonically active areas undergoing climatic fluctuations, such as the Himalaya Mountains and neighboring Qinghai-Tibetan Plateau (QTP). They provide biogeographic templates upon which endemic biodiversity can be mapped to infer diversification scenarios. Yet, this process can be somewhat opaque for the Himalaya, given substantial data gaps separating eastern and western regions. To help clarify, we evaluated phylogeographic and phylogenetic hypotheses for a widespread fish (Snowtrout: Cyprininae; Schizothorax) by sequencing 1,140 base pair of mtDNA cytochrome-b (cytb) from Central Himalaya samples (Nepal: N = 53; Bhutan: N = 19), augmented with 68 GenBank sequences (N = 60 Schizothorax/N = 8 outgroups). Genealogical relationships (N = 132) were analyzed via maximum likelihood (ML), Bayesian (BA), and haplotype network clustering, with clade divergence estimated via TimeTree. Snowtrout seemingly originated in Central Asia, dispersed across the QTP, then into Bhutan via southward-flowing tributaries of the east-flowing Yarlung-Tsangpo River (YLTR). Headwaters of five large Asian rivers provided dispersal corridors from Central into eastern/southeastern Asia. South of the Himalaya, the YLTR transitions into the Brahmaputra River, facilitating successive westward colonization of Himalayan drainages first in Bhutan, then Nepal, followed by far-western drainages subsequently captured by the (now) westward-flowing Indus River. Two distinct Bhutanese phylogenetic groups were recovered: Bhutan-1 (with three subclades) seemingly represents southward dispersal from the QTP; Bhutan-2 apparently illustrates northward colonization from the Lower Brahmaputra. The close phylogenetic/phylogeographic relationships between the Indus River (Pakistan) and western tributaries of the Upper Ganges (India/Nepal) potentially implicate an historic, now disjunct connection. Greater species-divergences occurred across rather than within-basins, suggesting vicariance as a driver. The Himalaya is a component of the Earth's largest glacial reservoir (i.e., the "third-pole") separate from the Arctic/Antarctic. Its unique aquatic biodiversity must be defined and conserved through broad, trans-national collaborations. Our study provides an initial baseline for this process.

Phylogenetic relationship and taxonomic status of Gymnocypris eckloni (Schizothoracinae) based on specific locus amplified fragments sequencing

Accurately delimiting phylogenetic relationships and taxonomic status is important for understanding species diversity and distributions and devising effective strategies for biodiversity conservation. However, species delimitation is controversial in Gymnocypris eckloni, a schizothoracine fish endemic to the Qinghai–Tibetan Plateau. The aim of this study is robustly identifying the phylogeny of G. eckloni in the Yellow River (YR) population and Qaidam basin (QB) population. The specific-locus amplified fragments sequencing (SLAF-seq) is employed with comprehensively sampling of schizothoracine fishes. In total, 350,181,802 clean reads and 5,114,096 SNPs are identified from SLAF-seq. Phylogenetic analysis recovers a non-monophyletic population of G. eckloni between YR and QB populations, representing an independent phylogenetic relationship between the two populations. Species delimitation analyses by SNAPPER and GMYC methods using the genome-wide SNP data confirm that their taxonomic statuses are separated. This study highlights the importance of further reconsidering clearer taxonomy, which would improve the genetic diversity conservation of Tibetan highland fishes.

Chromosome-level genome of Tibetan naked carp (Gymnocypris przewalskii) provides insights into Tibetan highland adaptation

Gymnocypris przewalskii, a cyprinid fish endemic to the Qinghai-Tibetan Plateau, has evolved unique morphological, physiological and genetic characteristics to adapt to the highland environment. Herein, we assembled a high-quality G. przewalskii tetraploid genome with a size of 2.03 Gb and scaffold N50 of 44.93 Mb, which was anchored onto 46 chromosomes. The comparative analysis suggested that gene families related to highland adaptation were significantly expanded in G. przewalskii. According to the G. przewalskii genome, we evaluated the phylogenetic relationship of 13 schizothoracine fishes, and inferred that the demographic history of G. przewalskii was strongly associated with geographic and eco-environmental alterations. We noticed that G. przewalskii experienced whole-genome duplication, and genes preserved post duplication were functionally associated with adaptation to high salinity and alkalinity. In conclusion, a chromosome-scale G. przewalskii genome provides an important genomic resource for teleost fish, and will particularly promote our understanding of the molecular evolution and speciation of fish in the highland environment.

Geological and climatic influences on population differentiation of the Phrynocephalus vlangalii species complex (Sauria: Agamidae) in the northern Qinghai-Tibet Plateau

Extremely heterogeneous topography and complex paleoclimate history of the Qinghai-Tibet Plateau (QTP) have a key role in promoting genetic divergence and lineage/species formation. Here, we sequenced one nuclear and three mitochondrial markers of 532 individuals from the entire range of the Phrynocephalus vlangalii species complex including two species, P. putjatai and P. vlangalii, endemic to the northern QTP. We integrated multilocus phylogeny, demographic analysis and geographic barrier detection to evaluate the population structure and dynamics. We found a new mitochondrial clade (PV-I) in the Gonghe County population of P. vlangalii, partial mitochondrial DNA replacement within P. vlangalii and complete mitochondrial DNA replacement between P. putjatai and P. vlangalii. Neutrality test, mismatch distribution analysis and Extended Bayesian Skyline Plot (EBSP) analysis all supported a significant expansion of the Qaidam Basin population of P. vlangalii (PV-II-2) from 0.091 to 0.026 Ma after Penultimate Glaciation. The uplift of the Arjin and Anyemanqen Mountains during the Kunhuang Movement (∼1.2 Ma) split populations of P. vlangalii in Akesai, Qaidam Basin and source of the Yellow River. The uplift of the Elashan Mountains during the second phase of the Qingzang Movement (∼2.5 Ma) contributed to the divergence of the Gonghe County population of P. vlangalii from other conspecific populations. The third phase of the Qingzang Movement (∼1.7 Ma) contributed to the divergence of the Xinghai population of P. vlangalii from P. putjatai and to the divergence of the northern populations of P. putjatai from the southern conspecific populations. Our data support the idea that the geological and climatic changes following the orogeny of the QTP may have promoted population differentiation and shaped the current population patterns of the P. vlangalii species complex in the northeastern QTP.

Convergent evolution misled taxonomy in schizothoracine fishes (Cypriniformes: Cyprinidae)

Highly specialized grade (HSG; genera Gymnocypris, Oxygymnocypris, Schizopygopsis, Platypharodon and Chuanchia) of the Schizothoracinae (Cypriniformes: Cyprinidae) are endemic to the Qinghai-Tibet Plateau (QTP). Previously, two distinct ecomorphs were recognized according to trophic traits. One was a limnetic omnivore with normal lower jaw morphology, terminal mouth, and moderate or dense gill rakers, mostly inhabiting in open water of lakes, including Gymnocypris and Oxygymnocypris. Another was a benthic feeder with inferior mouth, sparse gill rakers and sharp horny sheath on the lower jaw for scraping of attached prey off hard substrates, including Schizopygopsis, Platypharodon and Chuanchia. However, traditional taxonomy of HSG based on these trophic traits presented extensive conflicts with the molecular studies in recent years. The possible cause could be convergent evolution in morphology, retention of ancestral polymorphisms or mitochondrial introgression, but these hypotheses could not be assessed due to incomplete taxon sampling and only mitochondrial data employed in previous works. Here, we conducted the most comprehensive molecular analysis on HSG fishes to date, using four mitochondrial loci and 152,464 genome-wide SNPs, and including 21 of 24 putative species and one undescribed Schizopygopsis species. Both SNP and mtDNA trees confirmed extensive paraphyly of genera Gymnocypris and Schizopygopsis, where species often were clustered together by watershed instead of by genus. Basal split into the north clade B and the south clade C (ca. 3.03 Ma) approximately by the Tanggula-Tanitawen Mountains in SLAF tree coincided with a violent uplift of the QTP during the phase A of 'Qingzang movement' (ca. 3.6 Ma). Ancestral state reconstruction of the trophic ecomorph showed that the limnetic omnivore ecomorph had evolved repeatedly in clade B and C. Furthermore, we presented a striking case of convergent evolution between two 'subspecies' Gymnocypris chui chui and G. chui longimandibularis, which had diverged as early as two million years ago (ca. 2.42 Ma). Ecological analyses revealed that similar food utilization, particularly in zooplankton, was the main underlying driving force. This work showed an example of taxonomy with the most extensive errors at the genus/species levels due to convergent evolution and suggested that trophic traits could be misleading in fish taxonomy. Therefore, we propose a major generic revision for HSG species.

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

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.