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

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.

Population Genomic Analysis of Two Endemic Schizothoracins Reveals Their Genetic Differences and Underlying Selection Associated with Altitude and Temperature

Simple Summary

Diptychus maculates and Gymnodiptychus dybowskii are two rare aboriginal fishes in the Xinjiang Uygur Autonomous Region. In recent years, due to overfishing and habitat fragmentation caused by construction of water conservancy and hydropower projects, the fishery resources have decreased sharply. Understanding the genetic background is of great significance for resource protection. In this study, we revealed the similar trends of population genetic diversities in these two species collected from the Tarim River and the Yili River. In addition, outlier SNPs associated with temperature and altitude were detected in both of them, indicating that Schizothoracinae fishes represented by D. maculates and G. dybowskii were still under the selection pressure of plateau environments.

Abstract

Schizothoracins are a group of cyprinid fishes distributed throughout the Qinghai–Tibet Plateau, which can be classified in three grades: primitive, specialised and highly specialised according to adaptation ability to plateau environments. As the only specialised schizothoracins in Xinjiang, China, Diptychus maculates and Gymnodiptychus dybowskii are ideal materials for adaptive evolution research. Based on single-nucleotide polymorphism (SNP) loci detected by specific-locus amplified fragment (SLAF) technology, the genome-wide genetic diversities of these two species from nine sites in Xinjiang were evaluated. D.maculates in the Muzat River (BM) and G. dybowskii in the Kaidu River (LKG) presented the lowest genetic diversity levels, whereas D. maculates in the Kumarik River (BK) and G.dybowskii in the Kashi River (LK) were just the opposite. Cluster and principal component analysis demonstrated a distant genetic affinity between D. maculates in the Tashkurgan River (BT) and other populations. Outlier SNP loci were discovered both in D. maculates and G. dybowskii. The coalescent Bayenv and latent factor mixed model (LFMM) methods showed that a total of thirteen and eighteen SNPs in D. maculates were associated with altitude and temperature gradient, respectively. No intersection was revealed in G. dybowskii. The results indicated that D. maculates was subject to much greater divergent selection pressure. A strong signal of isolation-by-distance (IBD) was detected across D. maculates (Mantel test, rs = 0.65; p = 0.05), indicating an evident geographical isolation in the Tarim River. Isolation-by-environment (IBE) analysis implied that temperature and altitude selections were more intensive in D. maculates, with greater environmental variation resulting in weak gene flow.

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.