DNA barcodes and species distribution models evaluate threats of global climate changes to genetic diversity: a case study from Nanorana parkeri (Anura: Dicroglossidae)

Exploring Paleogene Tibet's warm temperate environments through target enrichment and phylogenetic niche modelling of Himalayan spiny frogs (Paini, Dicroglossidae)

The Cenozoic topographic development of the Himalaya-Tibet orogen (HTO) substantially affected the paleoenvironment and biodiversity patterns of High Asia. However, concepts on the evolution and paleoenvironmental history of the HTO differ massively in timing, elevational increase and sequence of surface uplift of the different elements of the orogen. Using target enrichment of a large set of transcriptome-derived markers, ancestral range estimation and paleoclimatic niche modelling, we assess a recently proposed concept of a warm temperate paleo-Tibet in Asian spiny frogs of the tribe Paini and reconstruct their historical biogeography. That concept was previously developed in invertebrates. Because of their early evolutionary origin, low dispersal capacity, high degree of local endemism, and strict dependence on temperature and humidity, the cladogenesis of spiny frogs may echo the evolution of the HTO paleoenvironment. We show that diversification of main lineages occurred during the early to Mid-Miocene, while the evolution of alpine taxa started during the late Miocene/early Pliocene. Our distribution and niche modelling results indicate range shifts and niche stability that may explain the modern disjunct distributions of spiny frogs. They probably maintained their (sub)tropical or (warm)temperate preferences and moved out of the ancestral paleo-Tibetan area into the Himalaya as the climate shifted, as opposed to adapting in situ. Based on ancestral range estimation, we assume the existence of low-elevation, climatically suitable corridors across paleo-Tibet during the Miocene along the Kunlun, Qiangtang and/or Gangdese Shan. Our results contribute to a deeper understanding of the mechanisms and processes of faunal evolution in the HTO.

Impacts of climate change on herpetofauna diversity in the Qinghai-Tibetan Plateau

Although numerous studies on the impacts of climate change on biodiversity have been published, only a handful are focused on the intraspecific level or consider population-level models (separate models per population). We endeavored to fill this knowledge gap relative to the Qinghai-Tibetan plateau (QTP) by combining species distribution modeling (SDMs) with population genetics (i.e., population-level models) and phylogenetic methods (i.e., phylogenetic tree reconstruction and phylogenetic diversity analyses). We applied our models to 11 endemic and widely distributed herpetofauna species inhabiting high elevations in the QTP. We aimed to determine the influence of environmental heterogeneity on species' responses to climate change, the magnitude of climate-change impacts on intraspecific diversity, and the relationship between species range loss and intraspecific diversity losses under 2 shared socioeconomic pathways (SSP245 and SSP585) and 3 future periods (2050s, 2070s, and 2090s). The effects of global climatic change were more pronounced at the intraspecific level (22% of haplotypes lost and 36% of populations lost) than the morphospecies level in the SSP585 climate change scenario. Maintenance of genetic diversity was in general determined by a combination of factors including range changes, species genetic structure, and the part of the range predicted to be lost. This is owing to the fact that the loss and survival of populations were observed in species irrespective of the predicted range changes (contraction or expansion). In the southeast (mountainous regions), climate change had less of an effect on range size (>100% in 3 species) than in central and northern QTP plateau regions (range size <100% in all species). This may be attributed to environmental heterogeneity, which provided pockets of suitable climate in the southeast, whereas ecosystems in the north and central regions were homogeneous. Generally, our results imply that mountainous regions with high environmental heterogeneity and high genetic diversity may buffer the adverse impacts of climate change on species distribution and intraspecific diversity. Therefore, genetic structure and characteristics of the ecosystem may be crucial for conservation under climate change.

A New Species of Nanorana (Anura: Dicroglossidae) from Northwestern Yunnan, China, with Comments on the Taxonomy of Nanorana arunachalensis and Allopaa

The genus Nanorana contains three subgenera, namely Nanorana, Paa, and Chaparana, and currently, there are four species known to science in Nanorana (Nanorana). In this study, we describe a new species belonging to the subgenus Nanorana from northwestern Yunnan, China. Phylogenetically, the new species, Nanorana laojunshanensissp. nov., is the sister to the clade of N. pleskei and N. ventripunctata. Morphologically, the new species can be distinguished from known congeners by the combination of following characters: present tympanum, equal fingers I and II, small body size, yellow ventral surface of limbs, distinct vomerine teeth, indistinct subarticular tubercles, head width greater than head length, slender supratympanic fold, absent dorsolateral fold, nuptial spines present on fingers I and II in adult males, absent vocal sac, and paired brown spines on the chest. Moreover, we suggest moving the genus Allopaa into Nanorana (Chaparana) and consider that N. arunachalensis is neither an Odorrana species nor a member of the subfamily Dicroglossinae (therefore Nanorana), but probably represents a distinct genus closely related to Ingerana or belongs to Ingerana, pending more data. Additionally, we consider that Nanorana minica deserves the rank of an independent subgenus, and we suggest assigning N. arnoldi, N. blanfordii, N. ercepeae, N. polunini, N. rarica, N. rostandi, N. vicina, N. xuelinensis, and N. zhaoermii into the subgenus Paa and placing N. kangxianensis, N. phrynoides, and N. sichuanensis in the subgenus Chaparana.

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.

Prediction of Potential Suitable Distribution Areas of Quasipaa spinosa in China Based on MaxEnt Optimization Model

Quasipaa spinosa is a large cold-water frog unique to China, with great ecological and economic value. In recent years, due to the impact of human activities on the climate, its habitat has been destroyed, resulting in a sharp decline in natural population resources. Based on the existing distribution records of Q. spinosa, this study uses the optimized MaxEnt model and ArcGis 10.2 software to screen out 10 factors such as climate and altitude to predict its future potential distribution area because of climate change. The results show that when the parameters are FC = LQHP and RM = 3, the MaxEnt model is optimal and AUC values are greater than 0.95. The precipitation of the driest month (bio14), temperature seasonality (bio4), elevation (ele), isothermality (bio3), and the minimum temperature of coldest month (bio6) were the main environmental factors affecting the potential range of the Q. spinosa. At present, high-suitability areas are mainly in the Hunan, Fujian, Jiangxi, Chongqing, Guizhou, Anhui, and Sichuan provinces of China. In the future, the potential distribution area of Q. spinosa may gradually extend to the northwest and north. The low-concentration emissions scenario in the future can increase the area of suitable habitat for Q. spinosa and slow down the reduction in the amount of high-suitability areas to a certain extent. In conclusion, the habitat of Q. spinosa is mainly distributed in southern China. Because of global climate change, the high-altitude mountainous areas in southern China with abundant water resources may be the main potential habitat area of Q. spinosa. Predicting the changes in the distribution patterns of Q. spinosa can better help us understand the biogeography of Q. spinosa and develop conservation strategies to minimize the impacts of climate change.

Phylogeography of the Plateau Pika (Ochotona curzoniae) in Response to the Uplift of the Qinghai-Tibet Plateau

The evolution and current distribution of species on the Qinghai-Tibet Plateau have been significantly impacted by historical occurrences, including the uplift of the plateau and the Quaternary climate upheaval. As a remnant species, the plateau pika (Ochotona curzoniae) is a great model for researching historical events. In this study, 302 samples from 42 sample sites were utilized to analyze the impact of historical events on the evolution and distribution pattern of plateau pikas. The genetic diversity, patterns of differentiation, and historical dynamics of the plateau pika were investigated using molecular markers that included four mitochondrial genes (COI, D-loop, Cytb, and 12S rRNA) and three nuclear genes (GHR, IRBP, and RAG1). The results showed that: (1) The genetic diversity of the plateau pika was high in the Tibetan Plateau (Hd = 0.9997, π = 0.01205), and the plateau pika evolved into five lineages that occupied different geographical areas, with lineage 1 (Group 1) in the south of the Yarlung Zangbo River, lineage 2 (Group 2) in the hinterland of the plateau, lineage 3 (Group 3) in the northeastern part of the plateau, lineage 4 (Group 4) in the Hengduan Mountains, and lineage 5 (Group 5) in the eastern part of the plateau. (2) The gene flow among the five lineages was low, and the differentiation level was high (Nm < 0.25; Fst > 0.25), indicating that the geographical barriers between the five lineages, such as the Yarlung Zangbo River, the Qaidam-Ghuong-Guide Basin, and the Lancang River, effectively promoted the population differentiation of the plateau pika. (3) The plateau pika first spread from the Hengduan Mountains to the entire Qinghai-Tibet Plateau and then conducted small-scale migration and dispersal in several refuges across the plateau in response to climate changes during the glacial and interglacial periods. (4) Except for Group 1 and Group 4, all the other populations exhibited a rapid expansion between 0.06 and 0.01 Mya, but the expansion was considerably delayed or halted by the effects of climate change during the last glacial maximum (0.02 Mya). Overall, the plateau pika on the Qinghai-Tibet Plateau exhibits high genetic diversity, and topographic obstacles, including mountains, valleys, and basins, created by the uplift of the plateau and climatic changes since the Quaternary period have played an important role in the differentiation and historical dynamics of the plateau pika population.

Assessment of Amphibians Vulnerability to Climate Change in China

Global climate change is considered to be one of the main threats to organisms. As poikilothermic animals, amphibians are in particular sensitive because they cannot adapt to the dramatic climate change through active physiological regulation. Using 104 representative species, the present study conducted an assessment of amphibians vulnerability to climate change in China through the combination of two approaches. Specifically, 18 vulnerability criteria belonging to five categories (i.e., thermal tolerance, individual reproductive, population diffusion and diversity, food and habitat, and climate conditions) were first selected and scored based on literatures and experts opinions. Species were then ranked into three levels of climate change vulnerability (i.e., high, moderate, and low) by calculating vulnerability scores and conducting natural breaks analyses, as well as performing a principal coordinate analysis (PCoA) and k-means cluster analyses, respectively. To integrate the two results, a matrix with the ranks from each result was developed to produce a final integrated list. Our results indicated that the 104 amphibian species were classified into three types by natural breaks, with 54 low vulnerable species, 41 moderately vulnerable species, and nine highly vulnerable species. Based on the results of PCoA and k-means cluster analyses, five species were highly vulnerable, 38 species were moderately vulnerable, and 61 species were low vulnerable. The combination of the two ranks suggested that 36 species such as Hyla tsinlingensis and Liangshantriton taliangensis were of low vulnerability, 54 species such as Echinotriton chinhaiensis and Hynobius chinensis were of moderate vulnerability, and 14 species such as Ichthyophis kohtaoensis and Zhangixalus prasinatus were of high vulnerability. Overall, our results indicated that climate change could have strong potential effects on amphibians in China. And the highly vulnerable species such as Ichthyophis kohtaoensis, Zhangixalus prasinatus, and Theloderma corticale should be the priority in future conservation activities.

Genetic diversity in frogs linked to past and future climate changes on the roof of the world

Mountains, representing storehouses of biodiversity, endemism and threatened species, are biodiversity hotspots of great conservation importance. However, increasing evidence indicates that mountain species throughout the world are responding to climate change, past or contemporary, by shifting their geographic distributions and patterns of genetic diversity, potentially affecting their adaptive capacity and increasing risk of extinction. Using the iconic high-elevation frog Nanorana parkeri as indicator, we showed how spatial analyses of climatic stability combined with genetic data allow unravelling amphibian responses to past and future climate changes on 'the roof of the world'-the Qinghai-Tibetan Plateau. We found that areas along the Yarlung Tsangpo Valley were climatically more stable relative to other regions, apparently serving as a large climatic refugium during Quaternary glaciations, but that these areas will likely be affected by future climate change. As populations closer to Quaternary refugia usually had higher genetic diversity, current genetic diversity can be explained in the largest part by distance to historically stable areas, outweighing other historical and contemporary factors. Along with the dynamics of suitable range, a fluctuating habitat fragmentation supported the pattern of historical changes in genetic diversity (N_e ) over time. Our results emphasize strong relationships between amphibian genetic diversity, past range dynamics and where to preserve suitable habitats in the face of future climate changes. More generally, our findings highlighted a central role of refugia during Quaternary climatic fluctuations, and how isolation from refugia may have modulated amphibian genetic diversity across the Qinghai-Tibetan Plateau.

The first complete mitochondrial genome sequence of Nanorana parkeri and Nanorana ventripunctata (Amphibia: Anura: Dicroglossidae), with related phylogenetic analyses

Members of the Nanorana genus (family Dicroglossidae) are often referred to as excellent model species with which to study amphibian adaptations to extreme environments and also as excellent keystone taxa for providing insights into the evolution of the Dicroglossidae. However, a complete mitochondrial genome is currently only available for Nanorana pleskei. Thus, we analyzed the complete mitochondrial genomes of Nanorana parkeri and Nanorana ventripunctata to investigate their evolutionary relationships within Nanorana and their phylogenetic position in the family Dicroglossidae. Our results showed that the genomes of N. parkeri (17,837 bp) and N. ventripunctata (18,373 bp) encode 13 protein‐coding genes (PCGs), two ribosomal RNA genes, 23 transfer RNA (tRNA) genes, and a noncoding control region. Overall sequences and genome structure of the two species showed high degree of similarity with N. pleskei, although the motif structures and repeat sequences of the putative control region showed clear differences among these three Nanorana species. In addition, a tandem repeat of the tRNA‐Met gene was found located between the tRNA‐Gln and ND2 genes. On both the 5′ and 3′‐sides, the control region possessed distinct repeat regions; however, the CSB‐2 motif was not found in N. pleskei. Based on the nucleotide sequences of 13 PCGs, our phylogenetic analyses, using Bayesian inference and maximum‐likelihood methods, illustrate the taxonomic status of Nanorana with robust support showing that N. ventripunctata and N. pleskei are more closely related than they are to N. parkeri. In conclusion, our analyses provide a more robust and reliable perspective on the evolutionary history of Dicroglossidae than earlier analyses, which used only a single species (N. pleskei).

Selection and environmental adaptation along a path to speciation in the Tibetan frog Nanorana parkeri

Tibetan frogs, Nanorana parkeri, are differentiated genetically but not morphologically along geographical and elevational gradients in a challenging environment, presenting a unique opportunity to investigate processes leading to speciation. Analyses of whole genomes of 63 frogs reveal population structuring and historical demography, characterized by highly restricted gene flow in a narrow geographic zone lying between matrilines West (W) and East (E). A population found only along a single tributary of the Yalu Zangbu River has the mitogenome only of E, whereas nuclear genes of W comprise 89-95% of the nuclear genome. Selection accounts for 579 broadly scattered, highly divergent regions (HDRs) of the genome, which involve 365 genes. These genes fall into 51 gene ontology (GO) functional classes, 14 of which are likely to be important in driving reproductive isolation. GO enrichment analyses of E reveal many overrepresented functional categories associated with adaptation to high elevations, including blood circulation, response to hypoxia, and UV radiation. Four genes, including DNAJC8 in the brain, TNNC1 and ADORA1 in the heart, and LAMB3 in the lung, differ in levels of expression between low- and high-elevation populations. High-altitude adaptation plays an important role in maintaining and driving continuing divergence and reproductive isolation. Use of total genomes enabled recognition of selection and adaptation in and between populations, as well as documentation of evolution along a stepped cline toward speciation.

Phylogeographic investigation and ecological niche modelling of the endemic frog species Nanorana pleskei revealed multiple refugia in the eastern Tibetan Plateau

The largest plateau Tibetan Plateau supplied an excellent opportunity to investigate the influence of the Pleistocene events on the high-elevation species. To test for the alternative hypotheses of Pleistocene glacial refugia, we used partial sequences of two mitochondrial genes and one nuclear gene to examine the phylogeographic patterns of the endemic frog species Nanorana pleskei across its known range in the eastern Tibetan Plateau, and conducted species distribution modelling (SDM) to explore changes of its distribution range through current and paleo periods. In all data sets, the species was divided into lineage north occupying open plateau platform and lineage south colonizing the mountainous plateau. The divergence of two major clades was estimated at the early Pleistocene. In mtDNA, lineage north contained northeastern and northwestern sublineages, and lineage south had two overlapping-distributed sublineages. Different lineages possessed distinct demographic characteristics, i.e., subdivision in the northeastern sublineage, historical bottleneck effects and recent expansions in the northwestern sublineage and the southeastern sublineage. SDMs depicted that stable suitable habitats had existed in the upper-middle streams of the Yellow River, Dadu River, Jinsha River and Yalong River. These regions were also recognized as the ancestral areas of different lineages. In conclusion, Nanorana pleskei lineages have probably experienced long-term separations. Stable suitable habitats existing in upper-middle streams of major rivers on the eastern Tibetan Plateau and distinct demographic dynamics of different lineages indicated that the lineages possessed independent evolutionary processes in multiple glacial refugia. The findings verified the profound effects of Pleistocene climatic fluctuations on the plateau endemic species.

Mountains too high and valleys too deep drive population structuring and demographics in a Qinghai-Tibetan Plateau frog Nanorana pleskei (Dicroglossidae)

Pleistocene glacial–interglacial climatic oscillations greatly shaped the current genetic structure of many species. However, geographic features may influence the impact of climatic cycling. Distinct geographic and environmental characters between northern and southern parts of the eastern Qinghai–Tibetan Plateau (EQTP) facilitate explorations into the impacts of geographic features on species. The northern parts of EQTP contain large areas of marsh, and the environment is rather homogeneous. In contrast, the southern EQTP harbors complex alpine valleys and a much more heterogeneous setting. We evaluate DNA sequence variation from both the mitochondrial and nuclear genomes in Nanorana pleskei, a species endemic to the EQTP. Hypothesis testing on the evolutionary history of N. pleskei indicates that northern populations can disperse freely, but alpine valleys isolate southern populations. Demographic histories between northern and southern populations also differ. Northern populations appear to have experienced population expansions, while southern frogs exhibit a far more stable demographic history. By combining climatic analyses and species' distribution models, our study suggests that geographic and environmental features drive the differences between the northern and southern EQTP.

Genetic diversity and evolutionary history of the Schizothorax species complex in the Lancang River (upper Mekong)

The genus Schizothorax (Cyprinidae), one of the most diverse genera of ichthyofauna of the Qinghai‐Tibetan Plateau (QTP), is a good candidate for investigating patterns of genetic variation and evolutionary mechanisms. In this study, sequences from the mitochondrial control region, the cytochrome b gene, and two nuclear genes were used to re‐examine the genetic diversity and investigate the evolutionary history of the Schizothorax species complex inhabiting the Lancang River. Three maternal clades were detected in the Schizothorax species complex, but frequent nuclear allele sharing also occurred among the three maternal clades. A discrepancy between topologies of mitochondrial and nuclear loci might result from introgression or/and incomplete lineage sorting. The divergence of the clades of the Schizothorax species complex was closely related to the Late Pliocene and Early Pleistocene orogenesis of the QTP and Southwest Mountains of China. Demographic analyses indicated that the species complex subsequently persisted in situ with stable populations during Pleistocene glacial cycling, which suggested that Pleistocene climate changes did not exert a remarkable influence on the species complex. Our study provides a comprehensive analysis of the genetic diversity and evolutionary history of the Schizothorax species complex in the Lancang River.