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

Integrative proteome and metabolome analyses reveal molecular basis of the tail resorption during the metamorphic climax of Nanorana pleskei

During the metamorphosis of anuran amphibians, the tail resorption process is a necessary and crucial change. One subject that has received relatively little or no attention is the expression patterns of proteins and metabolites in the different tail portions during metamorphosis, especially in highland amphibians. The mechanisms of tail resorption in three portions (the tip, middle and root) of the tail were investigated in N. pleskei G43 tadpole based on two omics (proteomic and metabolomic). Integrin αVβ3 was found to be high expressed in the distal portion of the tail, which could improve the sensitiveness to thyroid hormones in the distal portion of the tail. Muscle regression displayed a spatial pattern with stronger regression in distal and weaker one in proximal portion. Probably, this stronger regression was mainly performed by the proteases of proteasome from the active translation by ribosomes. The suicide model and murder model coexisted in the tail resorption. Meanwhile, fatty acids, amino acids, pyrimidine, and purine which derived from the breakdown of tissues can be used as building blocks or energy source for successful metamorphosis. Our data improved a better comprehension of the tail resorption mechanisms underlying the metamorphism of N. pleskei tadpole through identifying important participating proteins and metabolites.

The last of their kind: Is the genus Scutiger (Anura: Megophryidae) a relict element of the paleo-Transhimalaya biota?

The orographic evolution of the Himalaya-Tibet Mountain system continues to be a subject of controversy, leading to considerable uncertainty regarding the environment and surface elevation of the Tibetan Plateau during the Cenozoic era. As many geoscientific (but not paleontological) studies suggest, elevations close to modern heights exist in vast areas of Tibet since at least the late Paleogene, implicating the presence of large-scale alpine environments for more than 30 million years. To explore a recently proposed alternative model that assumes a warm temperate environment across paleo-Tibet, we carried out a phylogeographic survey using genomic analyses of samples covering the range of endemic lazy toads (Scutiger) across the Himalaya-Tibet orogen. We identified two main clades, with several, geographically distinct subclades. The long temporal gap between the stem and crown age of Scutiger may suggest high extinction rates. Diversification within the crown group, depending on the calibration, occurred either from the Mid-Miocene or Late-Miocene and continued until the Holocene. The present-day Himalayan Scutiger fauna could have evolved from lineages that existed on the southern edges of the paleo-Tibetan area (the Transhimalaya = Gangdese Shan), while extant species living on the eastern edge of the Plateau originated probably from the eastern edges of northern parts of the ancestral Tibetan area (Hoh Xil, Tanggula Shan). Based on the Mid-Miocene divergence time estimation and ancestral area reconstruction, we propose that uplift-associated aridification of a warm temperate Miocene-Tibet, coupled with high extirpation rates of ancestral populations, and species range shifts along drainage systems and epigenetic transverse valleys of the rising mountains, is a plausible scenario explaining the phylogenetic structure of Scutiger. This hypothesis aligns with the fossil record but conflicts with geoscientific concepts of high elevated Tibetan Plateau since the late Paleogene. Considering a Late-Miocene/Pliocene divergence time, an alternative scenario of dispersal from SE Asia into the East, Central, and West Himalaya cannot be excluded, although essential evolutionary and biogeographic aspects remain unresolved within this model.

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.

Genomic insights into local adaptation in the Asiatic toad Bufo gargarizans, and its genomic offset to climate warming

Genomic signatures of local adaptation have been identified in many species but remain sparsely studied in amphibians. Here, we explored genome-wide divergence within the Asiatic toad, Bufo gargarizans, to study local adaptation and genomic offset (i.e., the mismatch between current and future genotype-environment relationships) under climate warming scenarios. We obtained high-quality SNP data for 94 Asiatic toads from 21 populations in China to study spatial patterns of genomic variation, local adaptation, and genomic offset to warming in this wide-ranging species. Population structure and genetic diversity analysis based on high-quality SNPs revealed three clusters of B. gargarizans in the western, central-eastern, and northeastern portions of the species' range in China. Populations generally dispersed along two migration routes, one from the west to the central-east and one from the central-east to the northeast. Both genetic diversity and pairwise F ST were climatically correlated, and pairwise F ST was also correlated with geographic distance. Spatial genomic patterns in B. gargarizans were determined by the local environment and geographic distance. Global warming will increase the extirpation risk of B. gargarizans.

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.

The effect of long-term cold acclimation on redox state and antioxidant defense in the high-altitude frog, Nanorana pleskei

Cold hardiness is a key determinant of the distribution and abundance of ectothermic animals, and thermal acclimation can strongly influence stress tolerance phenotypes. However, the effect of cold acclimation on oxidative stress and antioxidant defenses is still not well understood. Here, we investigated the effects of long-term cold exposure (30 days at 4 °C in darkness versus 30 days at 20 °C in natural light) on the redox state and antioxidant defenses of the high-altitude frog, Nanorana pleskei, indigenous to the Tibetan plateau. We found that cold acclimation, under conditions mimicking winter, led to a significant increase in the ratio of oxidized glutathione (GSSG) to its reduced form (GSH) in liver and skeletal muscle tissues, suggesting that cold exposure induced oxidative stress in this species. Furthermore, malondialdehyde (MDA) contents were significantly augmented in heart, liver and muscle, indicating cold-related oxidative damage in these tissues. In the brain, GST activity, total antioxidant capacity (T-AOC), and vitamin C content showed a significant reduction after cold acclimation. In liver, an apparent decrease was also observed in the activities of SOD and GST, as well as T-AOC, whereas CAT and GPX activities showed a prominent increase in cold-acclimated groups. In kidney, there was a significant decrease in most antioxidant enzyme activities except for SOD and GST activity. In skeletal muscle, the activity of SOD, CAT, GR as well as T-AOC significantly decreased but GPX activity showed a significant increase in cold-acclimated frogs. These findings indicate that, in general, cold acclimation induces a suppression of the antioxidant defense system. Overall, our present study systematically describes the responses of antioxidant defenses to long-term cold acclimation and these findings contribute to extending the current understanding of the mechanisms of cold tolerance in high-altitude frogs.

Demographic history and genetic differentiation of an endemic and endangered Ulmus lamellosa (Ulmus)

BACKGROUND

Ulmus lamellosa (one of the ancient species of Ulmus) is an endemic and endangered plant that has undergone climatic oscillations and geographical changes. The elucidation of its demographic history and genetic differentiation is critical for understanding the evolutionary process and ecological adaption to forests in Northern China.

RESULTS

Polymorphic haplotypes were detected in most populations of U. lamellosa via DNA sequencing. All haplotypes were divided into three phylogeographic clades fundamentally corresponding to their geographical distribution, namely THM (Taihang Mountains), YM (Yinshan Mountains), and YSM (Yanshan Mountains) groups. The YSM group, which is regarded as ancestral, possessed higher genetic diversity and significant genetic variability in contrast to the YSM and YM groups. Meanwhile, the divergence time of intraspecies haplotypes occurred during the Miocene-Pliocene, which was associated with major Tertiary geological and/or climatic events. Different degrees of gene exchanges were identified between the three groups. During glaciation, the YSM and THM regions might have served as refugia for U. lamellosa. Based on ITS data, range expansion was not expected through evolutionary processes, except for the THM group. A series of mountain uplifts (e.g., Yanshan Mountains and Taihang Mountains) following the Miocene-Pliocene, and subsequently quaternary climatic oscillations in Northern China, further promoted divergence between U. lamellosa populations.

CONCLUSIONS

Geographical topology and climate change in Northern China played a critical role in establishing the current phylogeographic structural patterns of U. lamellosa. These results provide important data and clues that facilitate the demographic study of tree species in Northern China.

Phylogeography of Schizopygopsis malacanthus Herzenstein (Cypriniformes, Cyprinidae) in relation to the tectonic events and Quaternary climatic oscillations in the Shaluli Mountains Region

The tectonic events and Quaternary climatic oscillations in the Shaluli Mountains (the margin of the southeastern Tibetan Plateau) had an extensive effect on the genetic structure and distribution patterns of this region's terrestrial fauna and flora. It is not yet clear whether similar mechanisms influence this region's fish fauna. Schizopygopsis malacanthus is limited to high-elevation rivers and lakes, and exhibits distinct adaptations to the mountains of, and near to, the Qinghai-Tibet Plateau (QTP). Therefore, this species is a good candidate for investigating patterns of genetic variation resulting from palaeoenvironmental fluctuations in the Shaluli Mountains (China). Here, we used microsatellite and mitochondrial DNA control region sequences to analyze six populations of S. malacanthus collected from the Jinsha and Yalong River drainages. Genealogical analyses identified four maternal lineages and perhaps even four putative species, of which the Ouqu River lineage played a pivotal role during the course of the species' evolution. Two lineages from the Yalong River drainage did not cluster together, whereas those from different drainages grouped together, suggesting tectonic event impacts that possibly altered regional river drainages have been highly influential in population connectivity and gene flow. Population genetic analysis indicated that the geographic barriers and this species preference for higher elevations both played key roles in the divergence of S. malacanthus populations. Demographic tests suggested large-scale spatial synchrony in population fluctuations of S. malacanthus, accompanying dramatic Pleistocene climatic oscillations. It appears that palaeoenvironmental changes in the Shaluli Mountains influenced the distribution and evolution of studied S. malacanthus populations, which provide important information about the factors that influenced the phylogeographic history of this region's fish fauna. Additionally, our study also has implications for ongoing conservation of this vulnerable fish.

Phylogeny of spiny frogs Nanorana (Anura: Dicroglossidae) supports a Tibetan origin of a Himalayan species group

Recent advances in the understanding of the evolution of the Asian continent challenge the long-held belief of a faunal immigration into the Himalaya. Spiny frogs of the genus Nanorana are a characteristic faunal group of the Himalaya-Tibet orogen (HTO). We examine the phylogeny of these frogs to explore alternative biogeographic scenarios for their origin in the Greater Himalaya, namely, immigration, South Tibetan origin, strict vicariance. We sequenced 150 Nanorana samples from 62 localities for three mitochondrial (1,524 bp) and three nuclear markers (2,043 bp) and complemented the data with sequence data available from GenBank. We reconstructed a gene tree, phylogenetic networks, and ancestral areas. Based on the nuDNA, we also generated a time-calibrated species tree. The results revealed two major clades (Nanorana and Quasipaa), which originated in the Lower Miocene from eastern China and subsequently spread into the HTO (Nanorana). Five well-supported subclades are found within Nanorana: from the East, Central, and Northwest Himalaya, the Tibetan Plateau, and the southeastern Plateau margin. The latter subclade represents the most basal group (subgenus Chaparana), the Plateau group (Nanorana) represents the sister clade to all species of the Greater Himalaya (Paa). We found no evidence for an east-west range expansion of Paa along the Himalaya, nor clear support for a strict vicariance model. Diversification in each of the three Himalayan subclades has probably occurred in distinct areas. Specimens from the NW Himalaya are placed basally relative to the highly diverse Central Himalayan group, while the lineage from the Tibetan Plateau is placed within a more terminal clade. Our data indicate a Tibetan origin of Himalayan Nanorana and support a previous hypothesis, which implies that a significant part of the Himalayan biodiversity results from primary diversification of the species groups in South Tibet before this part of the HTO was uplifted to its recent heights.

Habitat Selection and Genetic Structure of the Endangered Frog Species Odorrana wuchuanensis (Anura: Ranidae)

Understanding the habitat selection and population genetic structure of an endangered species can play important roles in its protection. The Wuchuan odorous frog (Odorrana wuchuanensis) is endemic to the karst regions of southwest China. This frog is currently listed as "Critically Endangered" by the IUCN, but little is known about its habitat selection and population genetics. In this study, we conducted analyses of habitat selection with occurrence/absence sites and environmental data, and assessed the genetic structure between north and south populations in Guizhou provinces in China using three mitochondrial markers. The results revealed that the probability of this frog occupying cave habitats increased with higher average humidity in July and higher lowest temperature in January, but was negatively related to precipitation in January. Analyses of F statistics combined with analyses of median-joining haplotype networks and the phylogenetic tree showed low genetic differentiation between the two populations of O. wuchuanensis. Considering the small population size and geographic isolation because of the complex karst terrains, we suggest careful management practices are needed to protect this species.

Phylogeography of the Tibetan hamster Cricetulus kamensis in response to uplift and environmental change in the Qinghai-Tibet Plateau

AIM

The evolutionary process of an organism provides valuable data toward an understanding of the Earth evolution history. To investigate the relationship between the uplift of the Qinghai-Tibet Plateau (QTP) and mammalian evolution since the late Cenozoic, the geographic distribution of genetic variations in the Tibetan hamster Cricetulus kamensis was investigated using phylogeographical methods. In particular, population divergence, demographic history, genetic variation, and the prediction of species distribution area were investigated.

LOCATION

The Qinghai-Tibet Plateau.

METHODS

A total of 53 specimens, representing 13 geographic populations, were collected from the QTP. The phylogeographical pattern and demographic history of C. kamensis were analyzed, and the probable factors in the QTP uplift and the Quaternary glacial periods were inferred from one nuclear and four mitochondrial genes. Furthermore, the species distribution model (SDM) was used to predict changes in potentially suitable habitats since the last Interglacial.

RESULTS

Phylogenetic analysis demonstrated that two major genetic differentiations of the C. kamensis population occurred during the Early Pleistocene that were influenced by the Qing-Zang tectonic movement from the Middle Pliocene to the Early Pleistocene. Genetic distance between two major clades indicated low genetic divergence. Demographic history analysis showed that the C. kamensis population was affected by the Quaternary glacial period. SDM analysis indicated that C. kamensis was endemic to the QTP and the suitable habitat was affected by climate change, especially during the Last Glacial Maximum (LGM).

MAIN CONCLUSION

Our results indicated that the QTP uplift led to the population divergence of C. kamensis, and vicariance well accounted for the geographic distribution of genetic variation in C. kamensis as a result of genetic divergence and lack of gene flow. The genetic distance shows that C. alticola may be a subspecies of C. kamensis. Demographic history analysis suggests that the QTP was affected by the last glacial period. SDM analysis supports that almost the entire QTP is covered by a huge ice sheet during the LGM.

Landscape genetics reveal broad and fine-scale population structure due to landscape features and climate history in the northern leopard frog (Rana pipiens) in North Dakota

Prehistoric climate and landscape features play large roles structuring wildlife populations. The amphibians of the northern Great Plains of North America present an opportunity to investigate how these factors affect colonization, migration, and current population genetic structure. This study used 11 microsatellite loci to genotype 1,230 northern leopard frogs (Rana pipiens) from 41 wetlands (30 samples/wetland) across North Dakota. Genetic structure of the sampled frogs was evaluated using Bayesian and multivariate clustering methods. All analyses produced concordant results, identifying a major east-west split between two R. pipiens population clusters separated by the Missouri River. Substructuring within the two major identified population clusters was also found. Spatial principal component analysis (sPCA) and variance partitioning analysis identified distance, river basins, and the Missouri River as the most important landscape factors differentiating R. pipiens populations across the state. Bayesian reconstruction of coalescence times suggested the major east-west split occurred ~13-18 kya during a period of glacial retreat in the northern Great Plains and substructuring largely occurred ~5-11 kya during a period of extreme drought cycles. A range-wide species distribution model (SDM) for R. pipiens was developed and applied to prehistoric climate conditions during the Last Glacial Maximum (21 kya) and the mid-Holocene (6 kya) from the CCSM4 climate model to identify potential refugia. The SDM indicated potential refugia existed in South Dakota or further south in Nebraska. The ancestral populations of R. pipiens in North Dakota may have inhabited these refugia, but more sampling outside the state is needed to reconstruct the route of colonization. Using microsatellite genotype data, this study determined that colonization from glacial refugia, drought dynamics in the northern Great Plains, and major rivers acting as barriers to gene flow were the defining forces shaping the regional population structure of R. pipiens in North Dakota.

AutoSeqMan: batch assembly of contigs for Sanger sequences

With the wide application of DNA sequencing technology, DNA sequences are still increasingly generated through the Sanger sequencing platform. SeqMan (in the LaserGene package) is an excellent program with an easy-to-use graphical user interface (GUI) employed to assemble Sanger sequences into contigs. However, with increasing data size, larger sample sets and more sequenced loci make contig assemble complicated due to the considerable number of manual operations required to run SeqMan. Here, we present the 'autoSeqMan' software program, which can automatedly assemble contigs using SeqMan scripting language. There are two main modules available, namely, 'Classification' and 'Assembly'. Classification first undertakes preprocessing work, whereas Assembly generates a SeqMan script to consecutively assemble contigs for the classified files. Through comparison with manual operation, we showed that autoSeqMan saved substantial time in the preprocessing and assembly of Sanger sequences. We hope this tool will be useful for those with large sample sets to analyze, but with little programming experience. It is freely available at https://github.com/ Sun-Yanbo/autoSeqMan.

Understanding population structure and historical demography of Litsea auriculata (Lauraceae), an endangered species in east China

Detecting how historical and contemporary factors contribute to genetic divergence and genetic structure is a central question in ecology and evolution. We examine this question by intergrating population genetics with ecological niche modelling of Litsea auriculata (Lauraceae), which is endangered and native to east China. Geographical and environmental factors including climatic fluctuations since the last glacial maximum (LGM) have also contribute to population demography and patterns of genetic structure. L. auriculata populations underwent expansion after divergence and dramatically decreased to the current small size with relative population bottlenecks due to climate changes. Populations separated by physical geographical barrier including geographic distance and Yangtze River, as a result contemporary gene flow among L. auriculata populations showed drastic declines in comparison with historical gene flow, resulting in a high level of population divergence. Thus, patterns of genetic structure of L. auriculata can result from both geographic and environmental factors including climate changes. This information is helpful in forming conservation strategies for L. auriculata in China.

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.