River islands, refugia and genetic structuring in the endemic brown frog Rana kukunoris (Anura, Ranidae) of the Qinghai-Tibetan Plateau

Genetic diversity and haplotype distribution patterns analysis of cytb and RAG2 sequences in Rana hanluica from southern China

Rana hanluica: an endemic amphibian of China, is found in the hills and mountains south of the Yangtze River. In this comprehensive study, we collected 162 samples from 14 different localities to delve into the genetic diversity of Rana hanluica using mitochondrial Cytb and nuclear RAG2 as genetic markers. Our findings reveal that the Nanling Mountains, specifically regions like Jiuyi Shan, Jinggang Shan, Mang Shan, and Qiyun Shan, are genetic hotspots harboring remarkable diversity. The research results also indicate that there is gene flow among the various populations of the species, and no distinct population structure has formed, which may be due to migration. Moreover, populations in some regions, as well as the overall population, show signs of a possible genetic bottleneck, which we speculate may have been caused by climate change. However, given the exploratory nature of our study, further investigations are warranted to confirm these observations. Through phylogenetic analyses, we uncovered indications that R. hanluica might have originated within the Nanling region, dispersing along the east-west mountain ranges, with a significant contribution originating from Jiuyi Shan. The genetic distributions uncovered through our research reflect historical migratory patterns, evident in the distinct haplotypes of the RAG2 gene between the western and eastern parts of the studied area. Moreover, Heng Shan and Yangming Shan exhibited unique genetic signatures, possibly influenced by geographic isolation, which has shaped their distinct genotypes. The insights gained from this study hold profound implications for conservation efforts. By identifying regions rich in genetic diversity and crucial gene flow corridors, we can develop more effective conservation strategies. Preserving these genetically diverse areas, especially within the Nanling Mountains, is vital for maintaining the evolutionary potential of R. hanluica. In conclusion, our research has laid a solid foundation for understanding the genetic landscape of R. hanluica, shedding light on its origins, population structures, and evolutionary trajectories. This knowledge will undoubtedly guide future research endeavors and inform conservation strategies for this endemic amphibian.

Genomic analysis reveals deep population divergence in the water snake Trimerodytes percarinatus (Serpentes, Natricidae)

Although several phylogeographic studies of Asian snakes have been conducted, most have focused on pitvipers, with non-venomous snakes, such as colubrids or natricids, remaining poorly studied. The Chinese keelback water snake (Trimerodytes percarinatus Boulenger) is a widespread, semiaquatic, non-venomous species occurring in China and southeastern Asia. Based on mitochondrial DNA (mtDNA) and single nucleotide polymorphism (SNP) data, we explored the population genetic structure, genetic diversity, and evolutionary history of this species. MtDNA-based phylogenetic analysis showed that T. percarinatus was composed of five highly supported and geographically structured lineages. SNP-based phylogenetic analysis, principal component analysis, and population structure analysis consistently revealed four distinct, geographically non-overlapping lineages, which was different from the mtDNA-based analysis in topology. Estimation of divergence dates and ancestral area of origin suggest that T. percarinatus originated ~12.68 million years ago (95% highest posterior density: 10.36-15.96 Mya) in a region covering southwestern China and Vietnam. Intraspecific divergence may have been triggered by the Qinghai-Xizang Plateau uplift. Population demographics and ecological niche modeling indicated that the effective population size fluctuated during 0.5 Mya and 0.002 Mya. Based on the data collected here, we also comment on the intraspecific taxonomy of T. percarinatus and question the validity of the subspecies T. p. suriki.

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.

Impacts of Quaternary glaciation, geological history and geography on animal species history in continental East Asia: A phylogeographic review

Continental East Asia has a mild Pleistocene climate and a complex recent geological history. Phylogeographic studies of animals over the last 30 years have produced several distinctive patterns. Glaciation refugia are numerous and are not restricted to any particular regions. Most of them are localized and species-specific, although several large refugia, for example the mountains of SW China, are shared by multiple species and have refugia-within-refugia. Furthermore, postglaciation range expansion events vary greatly in time, scale and direction. Large-scale south-to-north post-LGM expansions are few and mostly occurred in the northern regions. Additionally, several unique geographic features, including the three-step terrain of China and the northern arid belt, have significant impacts on many species histories. Overall, the impacts of Pleistocene glaciations, particularly the LGM, on species history vary drastically from nondetectable to significant. The impacts are the least for species from the southwestern region and are most dominant for species from the north. Geological events play a more significant role in shaping species history than Pleistocene climatic changes. Phylogeographic patterns among animals species are highly consistent with those of plants. Future phylogeographic endeavour in East Asia should be hypothesis-driven and seek processes that underlie common patterns. The wide use of genomic data allow accurate estimates of historical population processes and exploration of older history beyond the Pleistocene.

Uncovering the determinants of biodiversity hotspots in China: Evidence from the drivers of multiple diversity metrics on insect assemblages and implications for conservation

Understanding large-scale biodiversity patterns and underlying mechanisms during the formation process is essential for guiding conservation efforts. However, previous studies on the identification and formation mechanism of diversity hotspots in China were often limited to a single (alpha) diversity metric, while multiple (beta or zeta) diversity has rarely been used for exploring drivers and conservation actions. Here, a comprehensive species distribution dataset consisting of representative families of three insect orders was compiled to explore biodiversity hotspots based on different algorithms. Furthermore, to assess the effects of environmental factors on hotspots, we fitted generalized additive mixed-effects models (GAMMs) for species richness, generalized dissimilarity models (GDMs) and multi-site generalized dissimilarity modeling (MS-GDM) for the total beta and zeta diversity. Our results showed that biodiversity hotspots were mainly concentrated in central and southern China, especially in mountainous areas with complex topography, which indicated the insects' affinity to montane environments. Further analyses based on multiple models showed that water-energy factors exerted the strongest explanatory power for the insect assemblage diversity in hotspots of both alpha and beta (or zeta) levels. Additionally, anthropogenic factors also exerted a significant effect on hotspots, and this effect was higher for beta diversity than for alpha diversity. Overall, our study elucidates a comprehensive analysis of the identification and underlying mechanism of biodiversity hotspots in China. Despite several limitations, we still believe that our findings can provide some new insights for conservation efforts in Chinese hotspots.

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.

Wind-dispersed seeds blur phylogeographic breaks: The complex evolutionary history of Populus lasiocarpa around the Sichuan Basin

The strength of phylogeographic breaks can vary among species in the same area despite being subject to the same geological and climate history due to differences in biological traits. Several important phylogeographic breaks exist around the Sichuan Basin in Southwest China but few studies have focused on wind-dispersed plants. Here, we investigated the phylogeographic patterns and the evolutionary history of Populus lasiocarpa, a wind-pollinated and wind-dispersed tree species with a circum-Sichuan Basin distribution in southwest China. We sequenced and analyzed three plastid DNA fragments (ptDNA) and eight nuclear microsatellites (nSSRs) of 265 individuals of P. lasiocarpa from 21 populations spanning the entire distribution range. Distribution patterns based on nSSR data revealed that there are three genetic groups in P. lasiocarpa. This is consistent with the three phylogeographic breaks (Sichuan Basin, the Kaiyong Line and the 105°E line), where the Sichuan basin acts as the main barrier to gene flow between western and eastern groups. However, the distribution pattern based on ptDNA haplotypes poorly matched the phylogeographic breaks, and wind-dispersed seeds may be one of the main contributing factors. Species distribution modelling suggested a larger potential distribution in the last glacial maximum with a severe bottleneck during the last interglacial. A DIYABC model also suggested a population contraction and expansion for both western and eastern lineages. These results indicate that biological traits are likely to affect the evolutionary history of plants, and that nuclear molecular markers, which experience higher levels of gene flow, might be better indicators of phylogeographic breaks.

Chromosome-level genome assembly of a high-altitude-adapted frog (Rana kukunoris) from the Tibetan plateau provides insight into amphibian genome evolution and adaptation

BACKGROUND

The high-altitude-adapted frog Rana kukunoris, occurring on the Tibetan plateau, is an excellent model to study life history evolution and adaptation to harsh high-altitude environments. However, genomic resources for this species are still underdeveloped constraining attempts to investigate the underpinnings of adaptation.

RESULTS

The R. kukunoris genome was assembled to a size of 4.83 Gb and the contig N50 was 1.80 Mb. The 6555 contigs were clustered and ordered into 12 pseudo-chromosomes covering ~ 93.07% of the assembled genome. In total, 32,304 genes were functionally annotated. Synteny analysis between the genomes of R. kukunoris and a low latitude species Rana temporaria showed a high degree of chromosome level synteny with one fusion event between chr11 and chr13 forming pseudo-chromosome 11 in R. kukunoris. Characterization of features of the R. kukunoris genome identified that 61.5% consisted of transposable elements and expansions of gene families related to cell nucleus structure and taste sense were identified. Ninety-five single-copy orthologous genes were identified as being under positive selection and had functions associated with the positive regulation of proteins in the catabolic process and negative regulation of developmental growth. These gene family expansions and positively selected genes indicate regions for further interrogation to understand adaptation to high altitude.

CONCLUSIONS

Here, we reported a high-quality chromosome-level genome assembly of a high-altitude amphibian species using a combination of Illumina, PacBio and Hi-C sequencing technologies. This genome assembly provides a valuable resource for subsequent research on R. kukunoris genomics and amphibian genome evolution in general.

Dispersal into the Qinghai-Tibet plateau: evidence from the genetic structure and demography of the alpine plant Triosteum pinnatifidum

Triosteum pinnatifidum Maxim., an alpine plant, is traditionally used for several medicinal purposes. Here, both chloroplast DNA sequences and nuclear low copy sequence markers were used to investigate the genetic diversity and population structure of T. pinnatifidum. Materials were collected from thirteen localities in the northeast Qinghai-Tibet Plateau (QTP) and adjacent highlands and advanced analytical toolkits were used to access their origin and range shifts. The results revealed a higher level of population differentiation based on chloroplast DNA (cpDNA) concatenated sequences compared with the nuclear DNA sequences (F _ST = 0.654 for cpDNA, F _ST = 0.398 for AT103), indicating that pollen flow was still extensive in T. pinnatifidum. A decline in haplotype variation was observed from the plateau edge and adjoining highlands toward the platform of the QTP. The hypothesis "dispersal into the QTP," proposing that T. pinnatifidum experienced migration from the plateau edge and adjacent highlands to the platform, was supported. These results were in line with the hypothesis that multiple refugia exist on the plateau edge and adjacent highlands rather than on the plateau platform. Our unimodal mismatch distribution, star-like network supported a recent expansion in T. pinnatifidum.

Cryptic diversity and phylogeography of the Rhabdophis nuchalis group (Squamata: Colubridae)

Previous studies, have found that the rapid uplift of the Tibetan plateau accelerated the diversification of species. However, there are few relevant biogeographic data for the Colubridae in this region. We conducted a comprehensive study of the Rhabdophis nuchalis Group, which presently contains four nominal species, R. nuchalis, R. pentasupralabialis, R. leonardi, and R. chiwen. Building upon previous studies with specimens we have recently examined, greater interspecific and intraspecific diversity has been revealed. Here we address three questions: (1) Do the intraspecific differences represent only geographic variation within lineages, or are there cryptic species? (2) What are the interspecific relationships among members of the R. nuchalis Group? (3) What has been the biogeographic history of this species group? To resolve these questions we used four mitochondrial gene sequences and one nuclear sequence to investigate the molecular phylogenetic and geographic relationships among populations. Our molecular analysis reveals cryptic species diversity within the R. nuchalis Group, and seven clades were identified in the analysis. Ancestral area estimation suggests that the R. nuchalis Group originated in the Hengduan Mountains approximately 6.24 Mya and expanded its range northward to the Qinling-Daba Mountains. The Sichuan Basin appears to have been a barrier to migration. Species divergence seems to have been related to the rapid uplift of the Qinghai-Tibet Plateau.

Out of the Hengduan Mountains: Molecular phylogeny and historical biogeography of the Asian water snake genus Trimerodytes (Squamata: Colubridae)

The Asian water snake genus Trimerodytes (formerly Sinonatrix) is endemic to East and Southeast Asia. Although several species have been included in various phylogenetic studies previously, the evolution and relationships among members of this genus as a whole remain unexplored. In this study, we report the sequencing two protein-coding mitochondrial gene fragments (MTCYB and ND2) and three nuclear genes (c-mos, NT3, and Rag1), reconstruct interspecific phylogeny, and explore biogeography for the genus Trimerodytes. Both Bayesian inference and maximum likelihood analyses consistently recover the monophyly of Trimerodytes with strong support, with T. yapingi the sister-group to the remaining species. The divergence date and ancestral area estimation suggest that Trimerodytes likely originated in Hengduan Mountains (eastern Tibetan Plateau) in western China at 23.93 Ma (95% HPD: 17.09-31.30), and intraspecific divergence began at about 4.23 Ma (95% HPD: 2.74-6.10). Analyses support the validity of T. yunnanensis.

Multilocus phylogeography of the brown-spotted pitviper Protobothrops mucrosquamatus (Reptilia: Serpentes: Viperidae) sheds a new light on the diversification pattern in Asia

Understanding the influence of geographical events and climate changes on genetic diversity is essential in explaining current patterns of genetic structure and geographic distribution of organisms. We inferred phylogenetic relationships, investigated historical demography, explored the evolutionary history, and clarified intraspecific taxonomy of Protobothrops mucrosquamatus, which is one of the commonest and most wide-ranging Asian pitvipers. A total of 184 samples from 54 localities were sequenced and analyzed for two mitochondrial gene fragments and two nuclear genes. Phylogenetic reconstruction based on mtDNA sequences revealed the existence of a minimum of five geographically structured and well-supported lineages within P. mucrosquamatus. Based on the mtDNA gene tree, and the geographic relationship between populations allied by matrilineal lineages, a complex longitudinal and latitudinal diversification pattern was uncovered in P. mucrosquamatus. The estimated date of the origin of the species (about 5.3 Ma) and divergence of the intraspecific lineages match the rapid uplifting of Qinghai-Xizang Plateau, and is also consistent with those of some other co-distributed Asian pitvipers. Formation of the two island lineages (Taiwan and Hainan) was generally congruent with the first isolation of the islands, but the two lineages showed different relationships with the continental Asian populations in comparison with some other pitvipers. Population historical demographic analyses, based on three methods, showed that all lineages have experienced slight population expansion in and around the Dali Glacial. Tests of intraspecific taxonomy indicated that no cryptic taxon is present within this widely distributed snake.

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 structure, cryptic speciation and demographic history of the sharpbelly (Hemiculter leucisculus), a freshwater habitat generalist from southern China

BACKGROUND

Species with broad distributions frequently divide into multiple genetic forms and may therefore be viewed as "cryptic species". Here, we used the mitochondrial cytochrome b (Cytb) and 12 nuclear DNA loci to investigate phylogeographic structures of the sharpbelly (Hemiculter leucisculus) in rivers in southern China and explored how the geological and climatic factors have shaped the genetic diversity and evolutionary history of this species.

RESULTS

Our mitochondrial phylogenetic analysis identified three major lineages (lineages A, B, and C). Lineages B and C showed a relatively narrower geographic distribution, whereas lineage A was widely distributed in numerous drainages. Divergence dates suggested that H. leucisculus populations diverged between 1.61-2.38 Ma. Bayesian species delimitation analysis using 12 nuclear DNA loci indicated the three lineages probably represented three valid taxa. Isolation-with-migration (IM) analysis found substantial gene flow has occurred among the three lineages. Demographic analyses showed that lineages B and C have experienced rapid demographic expansion at 0.03 Ma and 0.08 Ma, respectively.

CONCLUSIONS

Hemiculter leucisculus populations in drainages in southern China comprise three mtDNA lineages, and each of which may represent a separate species. Intense uplift of the Qinghai-Tibetan Plateau, evolution of Asian monsoons, changes in paleo-drainages, and poor dispersal ability may have driven the divergence of the three putative species. However, gene flow occurs among the three lineages. Climatic fluctuations have a prominent impact on the populations from the lineages B and C, but exerted little influence on the lineage A.

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.

Phylogeographic studies of schizothoracine fishes on the central Qinghai-Tibet Plateau reveal the highest known glacial microrefugia

Pleistocene climatic oscillations have greatly influenced the evolutionary history and distribution pattern of most extant species. However, their effects on species on the Qinghai-Tibet Plateau (QTP) are not well understood. To investigate the effects of past climatic shifts, particularly the Last Glacial Maximum (LGM), on plateau fish, we analysed the phylogeographic structure and demographic history of five closely related taxa of the subfamily Schizothoracinae, a representative endemic taxon of the QTP, from nine endorheic lakes on the central QTP and three peripheral exorheic rivers using the mitochondrial control region (D-loop) sequence and 12 microsatellite (SSR) markers. Phylogram from D-loop haplotypes revealed two well-supported lineages (North and South) separated by the Tanggula Mountains. The results from the D-loop and SSR revealed that endorheic populations possess high genetic diversity and a unique genetic structure. The most recent demographic expansion occurred post-LGM for most endorheic populations and in the last interglacial period for Siling Co and all exorheic populations. Phylogeographic structure, together with species distribution modelling, supports the scenario of multiple glacial refugia on the QTP during the LGM and suggests that Siling Co (4540 m asl) is a cryptic glacial microrefugia for plateau fish, which would be the highest glacial microrefugia known.

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.

A Complex System of Glacial Sub-Refugia Drives Endemic Freshwater Biodiversity on the Tibetan Plateau

Although only relatively few freshwater invertebrate families are reported from the Tibetan Plateau, the degree of endemism may be high. Many endemic lineages occur within permafrost areas, raising questions about the existence of isolated intra-plateau glacial refugia. Moreover, if such refugia existed, it might be instructive to learn whether they were associated with lakes or with more dynamic ecosystems such as ponds, wetlands, or springs. To study these hypotheses, we used pulmonate snails of the plateau-wide distributed genus Radix as model group and the Lake Donggi Cona drainage system, located in the north-eastern part of the plateau, as model site. First, we performed plateau-wide phylogenetic analyses using mtDNA data to assess the overall relationships of Radix populations inhabiting the Lake Donggi Cona system for revealing refugial lineages. We then conducted regional phylogeographical analyses applying a combination of mtDNA and nuclear AFLP markers to infer the local structure and demographic history of the most abundant endemic Radix clade for identifying location and type of (sub-)refugia within the drainage system. Our phylogenetic analysis showed a high diversity of Radix lineages in the Lake Donggi Cona system. Subsequent phylogeographical analyses of the most abundant endemic clade indicated a habitat-related clustering of genotypes and several Late Pleistocene spatial/demographic expansion events. The most parsimonious explanation for these patterns would be a scenario of an intra-plateau glacial refugium in the Lake Donggi Cona drainage system, which might have consisted of isolated sub-refugia. Though the underlying processes remain unknown, an initial separation of lake and watershed populations could have been triggered by lake-level fluctuations before and during the Last Glacial Maximum. This study inferred the first intra-plateau refugium for freshwater animals on the Tibetan Plateau. It thus sheds new light on the evolutionary history of its endemic taxa and provides important insights into the complex refugial history of a high-altitude ecosystem.

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.

Vicariance and Its Impact on the Molecular Ecology of a Chinese Ranid Frog Species-Complex (Odorrana schmackeri, Ranidae)

Paleogeological events and Pleistocene climatic fluctuations have had profound influences on the genetic patterns and phylogeographic structure of species in southern China. In this study, we investigated the population genetic structure and Phylogeography of the Odorrana schmackeri species complex, mountain stream-dwelling odorous frogs, endemic to southern China. We obtained mitochondrial sequences (1,151bp) of the complete ND2 gene and two flanking tRNAs of 511 individuals from 25 sites for phylogeographic analyses. Phylogenetic reconstruction revealed seven divergent evolutionary lineages, with mean pairwise (K2P) sequence distances from 7.8% to 21.1%, except for a closer ND2 distance (3.4%). The complex geological history of southern China drove matrilineal divergence in the O. schmackeri species complex into highly structured geographical units. The first divergence between lineage A+B and other lineages (C-G) had likely been influenced by the uplift of coastal mountains of Southeast China during the Mio-Pliocene period. The subsequent divergences between the lineages C-G may have followed the formation of the Three Gorges and the intensification of the East Asian summer monsoon during the late Pliocene and early Pleistocene. Demographic analyses indicated that major lineages A and C have been experienced recent population expansion (c. 0.045–0.245 Ma) from multiple refugia prior to the Last Glacial Maximum (LGM). Molecular analysis suggest that these seven lineages may represent seven different species, three described species and four cryptic species and should at least be separated into seven management units corresponding to these seven geographic lineages for conservation.

Phylogeography of Nanorana parkeri (Anura: Ranidae) and multiple refugia on the Tibetan Plateau revealed by mitochondrial and nuclear DNA

Quaternary climatic changes have been recognized to influence the distribution patterns and evolutionary histories of extant organisms, but their effects on alpine species are not well understood. To investigate the Pleistocene climatic oscillations on the genetic structure of amphibians, we sequenced one mitochondrial and three nuclear DNA fragments in Nanorana parkeri, a frog endemic to the Tibetan Plateau, across its distribution range in the southern plateau. Mitochondrial cytochrome b (Cytb) and three nuclear genes (c-Myc2, Rhod, and Tyr) revealed two distinct lineages (i.e. the lineages East and West), which were strongly geographically structured. The split of the two divergent lineages was dated back earlier than the Middle Pleistocene, probably being associated with climatic and ecological factors. Species distribution modeling, together with the phylogeographic structuring, supported the hypothesis of multiple refugia for N. parkeri on the Tibetan Plateau during the Pleistocene glaciations, and suggested the Yarlung Zangbo valley and the Kyichu catchment to be the potential refugia. Our findings indicate that Pleistocene climatic changes have had a great impact on the evolution and demographic history of N. parkeri. Our study has important implications for conservation of this and other frog species in the Tibetan Plateau.

Multiple genetic lineages challenge the monospecific status of the West African endemic frog family Odontobatrachidae

Background

Correct species identification is crucial in different fields of biology, and in conservation. The endemic West African frog family Odontobatrachidae currently contains a single described species, Odontobatrachus natator. From western Guinea to western Côte d'Ivoire it inhabits forests around waterfalls or cascades. Based on more than 130 specimens from 78 localities, covering the entire distribution, we investigated the molecular diversity of these frogs.

Results

Our analyses included mitochondrial and nuclear genes, with a concatenated alignment of 3527 base pairs. We detected high level of genetic differentiation with five distinct lineages or operational taxonomic units (OTUs). These OTUs were also identified by two different species delimitation approaches, Generalized Mixed Yule Coalescent (GMYC) and cluster algorithm. All OTUs occur in parapatry in the Upper Guinean forests. One OTU, assigned to the “true” Odontobatrachus natator, covers the largest distribution, ranging from the border region of western Sierra Leone-Guinea to eastern Liberia. Two OTUs are restricted to western Guinea (Fouta Djallon and foothills), while two others occur in eastern Guinea and the border region of Guinea-Liberia-Côte d'Ivoire. The OTU representing O. natator consists of two divergent subclades: one restricted to the Freetown Peninsula (Sierra Leone) and the other covering all populations further inland. Environmental niche models indicated that the restricted Freetown Peninsula population is separated by unsuitable habitat from remaining populations.

Conclusion

Geographic isolation of OTUs and molecular differences comparable to species level differentiation in other frog families indicate that O. natator contains cryptic species diversity. Respective distribution patterns most probably resulted from repeated changes of forest cover (contraction and expansion) over evolutionary timescales. The survival within forest refugia that have persisted through multiple drier periods and subsequent dispersal during wetter times may best explain the observed geographic distributions of OTUs. According to the IUCN Red List range criteria each OTU should be classified as “Endangered.” If the Freetown Peninsula “natator” population is recognized as a distinct species it would warrant recognition as “Critically Endangered.” The identification of cryptic lineages highlights the urgent need to protect these frogs, all of which are endemic to small areas within the Upper Guinean biodiversity hotspot.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0346-9) contains supplementary material, which is available to authorized users.

Phylogeography of Cephalotaxus oliveri (Cephalotaxaceae) in relation to habitat heterogeneity, physical barriers and the uplift of the Yungui Plateau

Habitat heterogeneity, physical barriers, and the uplift of the Yungui Plateau were found to deeply affect the phylogeographic pattern and evolutionary history of Cephalotaxus oliveri, a perennial conifer endemic to China. In this study, we explored the phylogeography using three chloroplast sequences (trnL-trnF, trnT-trnD and atpB-rbcL) in 22 natural populations of C. oliveri distributed throughout its range. The Yungui Plateau populations of C. oliveri were revealed to origin ca. 9.15Ma by molecular clock estimation, which is consistent with rapid uplift of the Qinghai-Tibetan Plateau (QTP) ca. 8-10Ma. Additionally, geological effects of the Yungui Plateau were suggested to promote the rapid intra-specific differentiation of C. oliveri in the Pliocene and Early Pleistocene. The relatively low level of genetic diversity (h=0.719, θ=1.17×10(-3)) and high population differentiation (NST=0.771 and GST=0.642) implied restricted gene flow among populations, which was confirmed by the Nested Clade Analysis (NCA). Mismatch distribution and haplotypes network provided evidences of recent demographic population expansion. Furthermore, the statistical dispersal-vicariance analysis indicated that the center of origin was in Central China. The comparison of haplotype distribution patterns indicated that the regions of HNHPS and HBLD were the potential refugia during the Pleistocene ice ages. Our results highlighted that habitat heterogeneity and physical barriers presenting in a species range can predict genetic patterns.

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

Anthropogenic global climate changes are one of the greatest threats to biodiversity. Distribution modeling can predict the effects of climate changes and potentially their effects on genetic diversity. DNA barcoding quickly identifies patterns of genetic diversity. As a case study, we use DNA barcodes and distribution models to predict threats under climate changes in the frog Nanorana parkeri, which is endemic to the Qinghai-Tibetan Plateau. Barcoding identifies major lineages W and E. Lineage W has a single origin in a refugium and Lineage E derives from three refugia. All refugia locate in river valleys and each greatly contributes to the current level of intraspecific genetic diversity. Species distribution models suggest that global climate changes will greatly influence N. parkeri, especially in the level of genetic diversity, because two former refugia will fail to provide suitable habitat. Our pipeline provides a novel application of DNA barcoding and has important implications for the conservation of biodiversity in southern areas of the Qinghai-Tibetan Plateau.

Evolutionary diversifications of plants on the Qinghai-Tibetan Plateau

The Qinghai-Tibetan Plateau (QTP) is the highest and one of the most extensive plateaus in the world. Phylogenetic, phylogeographic, and ecological studies support plant diversifications on the QTP through multiple mechanisms such as allopatric speciation via geographic isolation, climatic oscillations and divergences, pollinator-mediated isolation, diploid hybridization and introgression, and allopolyploidy. These mechanisms have driven spectacular radiations and/or species diversifications in various groups of plants such as Pedicularis L., Saussurea DC., Rhododendron L., Primula L., Meconopsis Vig., Rhodiola L., and many lineages of gymnosperms. Nevertheless, much work is needed toward understanding the evolutionary mechanisms of plant diversifications on the QTP. Well-sampled biogeographic analyses of the QTP plants in the broad framework of the Northern Hemisphere as well as the Southern Hemisphere are still relatively few and should be encouraged in the next decade. This paper reviews recent evidence from phylogenetic and biogeographic studies in plants, in the context of rapid radiations, mechanisms of species diversifications on the QTP, and the biogeographic significance of the QTP in the broader context of both the Northern and Southern Hemisphere biogeography. Integrative multidimensional analyses of phylogeny, morphological innovations, geography, ecology, development, species interactions and diversifications, and geology are needed and should shed insights into the patterns of evolutionary assembly and radiations in this fascinating region.

A phylogeographic evaluation of the Amolops mantzorum species group: cryptic species and plateau uplift

Phylogeographic analysis brings organism phylogeny, regional geological/climatic history, and population demography together, and provides insights into species history and speciation. Using DNA sequence data from a mitochondrial marker (cytochrome b) and a nuclear marker (cmyc intron 2), we examined phylogeography of the Amolops mantzorum species group. We first constructed Bayesian and maximum-likelihood gene trees and medium-joint networks for the recovered haplotypes, and estimated divergence time for each major lineage. Species boundaries were then delineated using the general mixed Yule-coalescent model (GMYC) and a Bayesian species-delimitation method (BP&P). Potential gene flow between putative species was also estimated using the isolation-with-migration model. Furthermore, species-distribution modeling was used to probe linkage between climatic changes and speciation. Lastly, extended Bayesian skyline plotting was employed to reveal historical demography for each putative species. Our analyses clearly delimit nine species in the group, including four well-recognized species and five putative species, of which three are often listed as synonyms of A. mantzorum. The nominal species A. mantzorum may in fact include two cryptic species. Divergence-time estimates align the speciation events with the recent intense uplift of the Tibetan Plateau in the last 3.6 million years. The species-distribution modeling also reveals different habitat preferences among species that are potentially linked to climatic changes associated with the uplift. Furthermore, three species sustained continuous population growth through the last glaciation, while others expanded only after the last glaciation. The eastern escarpment of Tibet is a biodiversity hotspot; its rich species diversity, unique topography, and well-established geological history provide an excellent system for evolutionary studies.