Phylogeographical structure and demographic expansion in the endemic alpine stream salamander (Hynobiidae: Batrachuperus) of the Qinling Mountains

Comparative phylogeography reveals dissimilar genetic differentiation patterns in two sympatric amphibian species

Global climate change is expected to have a profound effect on species distribution. Due to the temperature constraints, some narrow niche species could shift their narrow range to higher altitudes or latitudes. In this study, we explored the correlation between species traits, genetic structure, and geographical range size. More specifically, we analyzed how these variables are affected by differences in fundamental niche breadth or dispersal ability in the members of two sympatrically distributed stream-dwelling amphibian species (frog, Quasipaa yei; salamander, Pachyhynobius shangchengensis), in Dabie Mountains, East China. Both species showed relatively high genetic diversity in most geographical populations and similar genetic diversity patterns (JTX, low; BYM, high) correlation with habitat changes and population demography. Multiple clustering analyses were used to disclose differentiation among the geographical populations of these two amphibian species. Q. yei disclosed the relatively shallow genetic differentiation, while P. shangchengensis showed an opposite pattern. Under different historical climatic conditions, all ecological niche modeling disclosed a larger suitable habitat area for Q. yei than for P. shangchengensis; these results indicated a wider environment tolerance or wider niche width of Q. yei than P. shangchengensis. Our findings suggest that the synergistic effects of environmental niche variation and dispersal ability may help shape genetic structure across geographical topology, particularly for species with extremely narrow distribution.

Evidence for Admixture and Rapid Evolution During Glacial Climate Change in an Alpine Specialist

The pace of current climate change is expected to be problematic for alpine flora and fauna, as their adaptive capacity may be limited by small population size. Yet, despite substantial genetic drift following post-glacial recolonization of alpine habitats, alpine species are notable for their success surviving in highly heterogeneous environments. Population genomic analyses demonstrating how alpine species have adapted to novel environments with limited genetic diversity remain rare, yet are important in understanding the potential for species to respond to contemporary climate change. In this study, we explored the evolutionary history of alpine ground beetles in the Nebria ingens complex, including the demographic and adaptive changes that followed the last glacier retreat. We first tested alternative models of evolutionary divergence in the species complex. Using millions of genome-wide SNP markers from hundreds of beetles, we found evidence that the N. ingens complex has been formed by past admixture of lineages responding to glacial cycles. Recolonization of alpine sites involved a distributional range shift to higher elevation, which was accompanied by a reduction in suitable habitat and the emergence of complex spatial genetic structure. We tested several possible genetic pathways involved in adaptation to heterogeneous local environments using genome scan and genotype-environment association approaches. From the identified genes, we found enriched functions associated with abiotic stress responses, with strong evidence for adaptation to hypoxia-related pathways. The results demonstrate that despite rapid demographic change, alpine beetles in the N. ingens complex underwent rapid physiological evolution.

Geological events and climate change drive diversification and speciation of mute cicadas in eastern continental Asia

The poor mobility of nymphs living underground, usually for many years and the weak flying ability of adults make cicadas unique for evolutionary biology and bio-geographical study. Cicadas of the genus Karenia are unusual in Cicadidae in lacking the timbals that produce sound. Population differentiation, genetic structure, dispersal and evolutionary history of the eastern Asian mute cicada Karenia caelatata were investigated based on morphological, acoustic and molecular data. The results reveal a high level of genetic differentiation in this species. Six independent clades with nearly unique sets of haplotypes corresponding to geographically isolated populations are recognized. Genetic and geographic distances are significantly correlated among lineages. The phenotypic differentiation is generally consistent with the high levels of genetic divergence across populations. Results of ecological niche modeling suggest that the potential distribution range of this mountain-habitat specialist during the Last Glacial Maximum was broader than its current range, indicating this species had benefited from the climate change during the early Pleistocene in southern China. Geological events such as orogeny in Southwest China and Pleistocene climate oscillations have driven the differentiation and divergence of this species, and basins, plains and rivers function as natural "barriers" to block the gene flow. Besides significant genetic divergence being found among clades, the populations occurring in the Wuyi Mountains and the Hengduan Mountains are significantly different in the calling song structure from other populations. This may have resulted from significant population differentiation and subsequent adaptation of related populations. We conclude that ecological differences in habitats, coupled with geographical isolation, have driven population divergence and allopatric speciation. This study provides a plausible example of incipient speciation in Cicadidae and improves understanding of population differentiation, acoustic signal diversification and phylogeographic relationships of this unusual cicada species. It informs future studies on population differentiation, speciation and phylogeography of other mountain-habitat insects in the East Asian continent.

Pleistocene climate and geomorphology drive the evolution and phylogeographic pattern of Triplophysa robusta (Kessler, 1876)

Montane systems provide excellent opportunities to study the rapid radiation influenced by geological and climatic processes. We assessed the role of Pleistocene climatic oscillations and mountain building on the evolution history of Triplophysa robusta, a cold-adapted species restricted to high elevations in China. We found seven differentiated sublineages of T. robusta, which were established during the Mid Pleistocene 0.87–0.61 Mya. The species distribution modeling (SDM) showed an expansion of T. robusta during the Last Glacial Maximum (LGM) and a considerable retraction during the Last Interglacial (LIG). The deep divergence between Clade I distributed in Qinling Mountains and Clade II in Northeastern Qinghai-Tibet Plateau (QTP) was mainly the result of a vicariance event caused by the rapid uplifting of Qinling Mountains during the Early Pleistocene. While the middling to high level of historical gene flow among different sublineages could be attributed to the dispersal events connected to the repetition of the glacial period during the Pleistocene. Our findings suggested that frequent range expansions and regressions due to Pleistocene glaciers likely have been crucial for driving the phylogeographic pattern of T. robusta. Finally, we urge a burning question in future conservation projection on the vulnerable cold-adapted species endemic to high elevations, as they would be negatively impacted by the recent rapid climate warming.

Metabolites with antioxidant and α-glucosidase inhibitory activities produced by the endophytic fungi Aspergillus niger from Pachysandra terminalis

One new compound and 13 known compounds were isolated from Aspergillus niger, a plant endophytic fungus of Pachysandra terminalis collected from Qinling Mountains, Xi'an, China. The structure of new compound 1 was classically determined by extensive spectroscopic analysis. Compounds 5, 6, 8, and 14 were firstly reported from Aspergillus, while compound 2 was isolated from A. niger for the first time. All isolated compounds were further evaluated for their antioxidant and α-glucosidase inhibitory activities. Compounds 2 and 3 exhibited significant antioxidant activities with IC50 values of 31.64 μm and 24.32 μm, respectively, similar to the positive control ascorbic acid. Additionally, compound 1 displayed remarkable inhibitory activity against α-glucosidase with an IC50 value of 96.25 μm, which was 3.4-fold more potent than that of the positive control acarbose. Compound 1 has great potential for development as a new lead compound owing to its simple structure and remarkable biological activity.

The phylogeographical pattern of the Amur minnow Rhynchocypris lagowskii (Cypriniformes: Cyprinidae) in the Qinling Mountains

In this study, the phylogeographical pattern of the Amur minnow (Rhynchocypris lagowskii) widely distributed in the cold freshwaters of the Qinling Mountains was examined. A total of 464 specimens from 48 localities were sequenced at a 540-bp region of the mitochondrial cytochrome b (Cytb) gene, and 69 haplotypes were obtained. The mean ratio of the number of synonymous and nonsynonymous substitutions per site (dN/dS) was 0.028 and indicated purifying selection. Haplotype diversity (h) and nucleotide diversity (π) of natural populations of R. lagowskii varied widely between distinct localities. Phylogenetic trees based on Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) methods, and network analysis showed five well-differentiated lineages, but these did not completely correspond to localities and geographic distribution. Meanwhile, analysis of molecular variances (AMOVA) indicated the highest proportion of genetic variation was attributed to the differentiation between populations rather than by our defined lineages. In addition, there was no significant correlation between the pairwise Fst values and geographic distance (p > .05). Based on the molecular clock calibration, the time to the most recent common ancestor (TMRCA) was estimated to have emerged from the Late Miocene to the Early Pleistocene. Finally, the results of demographic history based on the neutrality test, mismatch distribution, and Bayesian skyline plot (BSP) analyses showed that collectively, the populations were stable during the Pleistocene while one lineage (lineage E) probably underwent a slight contraction during the Middle Pleistocene and a rapid expansion from the Middle to the Late Pleistocene. Therefore, the study suggests the current phylogeographical pattern of R. lagowskii was likely shaped by geological events that led to vicariance followed by dispersal and secondary contact, river capture, and climatic oscillation during the Late Miocene to the Early Pleistocene in the Qinling Mountains.

A Phylogeographical Analysis of the Beetle Pest Species Callosobruchus chinensis (Linnaeus, 1758) in China

Callosobruchus chinensis (Coleoptera Bruchidae), is a pest of different varieties of legumes. In this paper, a phylogeographical analysis of C. chinensis was conducted to provide knowledge for the prevention and control of C. chinensis. A total of 224 concatenated mitochondrial sequences were obtained from 273 individuals. Suitable habitat shifts were predicted by the distribution modelling (SDM). Phylogeny, genetic structure and population demographic history were analyzed using multiple software. Finally, the least-cost path (LCP) method was used to identify possible dispersal corridors and genetic connectivity. The SDM results suggested that the distribution of C. chinensis experienced expansion and contraction with changing climate. Spatial distribution of mtDNA haplotypes showed there was partial continuity among different geographical populations of C. chinensis, except for the Hohhot (Inner Mongolia) population. Bayesian skyline plots showed that the population had a recent expansion during 0.0125 Ma and 0.025 Ma. The expansion and divergent events were traced back to Quaternary glaciations. The LCP method confirmed that there were no clear dispersal routes. Our findings indicated that climatic cycles of the Pleistocene glaciations, unsuitable climate and geographic isolation played important roles in the genetic differentiation of C. chinensis. Human activities weaken the genetic differentiation between populations. With the change in climate, the suitable areas of C. chinensis will disperse greatly in the future.

UCE Phylogenomics, detection of a putative hybrid population, and one older mitogenomic node age of Batrachuperus salamanders

The prevalence of incomplete lineage sorting complicates the examination of hybridization and species-level paraphyly with gene trees of a small number of loci. In Asian mountain salamanders of the genus Batrachuperus, possible hybridization and species paraphyly had been identified by utilizing mitochondrial genealogy and fixed allozyme differences. Here we sampled 2909 UCEs in 44 local populations from all six Batrachuperus species, inferred gene and species trees, compared them with mitochondrial and allozyme results, and examined the potential hybridization and species paraphyly. The clustering pattern of single-locus trees, increased proportion of heterozygous SNPs, allele frequency-based migration edge estimation, and intrapopulation long branches (as expected from an increase of genetic lineage and nucleotide diversity) support that an eastern B. karlschmidti population has experienced admixture with B. tibetanus. On the 2909-UCE concatenated and species trees, lower nodal supports were observed when similar proportions of loci agreed with alternative topologies, i.e., a reciprocal monophyly between a Pengxian lineage and the remainder of B. pinchonii (0.379) or a paraphyly of the latter with respect to Pengxian (0.362). The UCE phylogenomics agreed with the relatively recent groupings in the allozyme dendrogram. Despite incomplete lineage sorting, the mitochondrial trees were similar to the UCE trees for deeper relationships of the genus. However, one significant branch-length level discordance was identified. The branch between the common ancestor of B. daochengensis and B. yenyuanensis and common ancestor of the genus was approximately three times shorter on the mitochondrial tree than on the UCE tree, suggesting that the split of the mitochondrial lineages was likely a few million years earlier than the split of species. This finding supports considering possible ancestral polymorphism when interpreting different divergence dates estimated from mitochondrial and genome-wide data.

Evolutionary history of the oriental fire-bellied toad (Bombina orientalis) in Northeast China

The evolutionary history of a species is generally affected by the combination of geological events and climate fluctuations. By analyzing the population features, genetic structure and the effective population historical dynamics of existing species, the population evolutionary history can be reestablished. In recent years, geological evidence shows that the Yilan-Yitong fault zone located in Northeast Asia experienced strong and frequent geological changes in the late Quaternary period. Species population history has been shaped by the combination of the complex climatic conditions of the Quaternary and Pleistocene glacial interglacial cycles and palaeogeological events in Northeast Asia and it has become a research focus for evolutionary biology researchers. In this study, mitochondrial and microsatellite molecular markers were used to reveal the population features, genetic structure, and the effective population historical dynamics of the Oriental fire-bellied toad (Bombina orientalis). The results showed that the strong seismic activity of the Yilan-Yitong fault zone in the late Quaternary period was the main reason for the population differentiation of Oriental fire-bellied toad in northeast China. The Quaternary Pleistocene glacial interglacial cycles led to the significant bottleneck effect of the western population located in the Maoer mountain area. As a result, the western population has low genetic diversity. Recent gene flow between eastern and western populations and historical evidence of population expansion proved that the dispersal behavior of the western populations was the main cause of the low genetic diversity and mitochondrial and nuclear discordance. Human economic activity may be the mainly driving factor. These evidences showed that the comprehensive influence of geology, climate, human activities and other factors should be considered in the process of exploring the evolutionary history of species.

Spatial Genetic Structure and Demographic History of the Wild Boar in the Qinling Mountains, China

Simple Summary

The wild boar is native to the temperate region of Eurasia, which is now one of the most widely distributed mammals worldwide. The recent expansion in the wild boar population has attracted a lot of attention, which may cause great damage to ecosystems. Elucidating the patterns of the population structure, genetic diversity, population origin, and colonization route of wild boar is very helpful in the conservation and management of wild populations. Phylogeographic analysis has proven to be a powerful tool. Here, 82 samples of wild boars in 16 sampling locations were collected in Qinling Mountains (QM). Genetic analysis was conducted based on the mitochondrial control region and nuclear genes. The level of genetic diversity of wild boars in QM was lower than the total population in East Asia, but higher than European population. No obvious phylogeographic pattern were found. The effective population size was under demographic equilibrium in the past.

Abstract

Species dispersal patterns and population genetic structure can be influenced by geographical features. Qinling Mountains (QM) provide an excellent area for phylogeographic study. The phylogeography of Asian-wide wild boars revealed the colonization route. However, the impact of the QM on genetic diversity, genetic structure and population origin is still poorly understood. In this study, genetic analysis of wild boar in the QM was conducted based on the mitochondrial control region (943 bp) and twelve microsatellite loci of 82 individuals in 16 sampling locations. Overall genetic haplotype diversity was 0.86, and the nucleotide diversity was 0.0079. A total of 17 new haplotypes were detected. The level of genetic diversity of wild boars in QM was lower than in East Asia, but higher than in Europe. Phylogenetic analysis showed the weak genetic divergence in QM. Mismatch analysis, neutrality tests, and Bayesian Skyline Plot (BSP) results revealed that the estimates of effective population size were under demographic equilibrium in the past. Spatial analysis of molecular variance indicated no obvious phylogeographic structure.

Drainage basins serve as multiple glacial refugia for alpine habitats in the Sierra Nevada Mountains, California

The evolutionary histories of alpine species are often directly associated with responses to glaciation. Deep divergence among populations and complex patterns of genetic variation have been inferred as consequences of persistence within glacier boundaries (i.e., on nunataks), while shallow divergence and limited genetic variation are assumed to result from expansion from large refugia at the edge of ice shields (i.e., massifs de refuge). However, for some species, dependence on specific microhabitats could profoundly influence their spatial and demographic response to glaciation, and such a simple dichotomy may obscure the localization of actual refugia. In this study, we use the Nebria ingens complex (Coleoptera: Carabidae), a water-affiliated ground beetle lineage, to test how drainage basins are linked to their observed population structure. By analysing mitochondrial COI gene sequences and genome-wide single nucleotide polymorphisms, we find that the major drainage systems of the Sierra Nevada Mountains in California best explain the population structure of the N. ingens complex. In addition, we find that an intermediate morphotype within the N. ingens complex is the product of historical hybridization of N. riversi and N. ingens in the San Joaquin basin during glaciation. This study highlights the importance of considering ecological preferences in how species respond to climate fluctuations and provides an explanation for discordances that are often observed in comparative phylogeographical studies.

Evolutionary history of endangered and relict tree species Dipteronia sinensis in response to geological and climatic events in the Qinling Mountains and adjacent areas

Geological and climatic events are considered to profoundly affect the evolution and lineage divergence of plant species. However, the evolutionary histories of tree species that have responded to past geological and climate oscillations in central China's mountainous areas remain mostly unknown. In this study, we assessed the evolutionary history of the endangered and relict tree species Dipteronia sinensis in the Qinling Mountains (QM) and adjacent areas in East Asia based on variations in the complete chloroplast genomes (cpDNA) and reduced-genomic scale single nucleotide polymorphisms (SNPs). Population structure and phylogenetic analysis based on the cpDNA variations suggested that D. sinensis could be divided into two intraspecific genetic lineages in the eastern and western sides of the QM (EQM and WQM, respectively) in East Asia. Molecular dating suggested that the intraspecific divergence of D. sinensis occurred approximately 39.2 million years ago during the later Paleogene. It was significantly correlated with the orogeny of the QM, where the formation of this significant geographic barrier in the region may have led to the divergence of independent lineages. Bayesian clustering and demographic analysis showed that intraspecific gene flow was restricted between the EQM and WQM lineages. Isolation-with-migration analysis indicated that the two genetic lineages experienced significant demographic expansions after the Pleistocene ice ages. However, the genetic admixture was determined in some populations between the two lineages by the large scale of SNP variations due to DNA incompatibility, the large significant population size, and rapid gene flow of nuclear DNA markers. Our results suggest that two different conservation and management units should be constructed for D. sinensis in the EQM and WQM areas. These findings provide novel insights into the unprecedented effects of tectonic changes and climatic oscillations on lineage divergence and plant population evolution in the QM and adjacent areas in East Asia.

Two new species of Cerapanorpa (Mecoptera, Panorpidae) from the Qinling and Minshan mountains

Two new species of Cerapanorpa Gao, Ma & Hua, 2016 are described from the Qinling and Minshan mountains, respectively. Cerapanorpa qinlingensis sp. nov. can be readily distinguished from its congeners by the elongate hypovalves and the extremely developed basal process of gonostylus in male genitalia. Cerapanorpa minshana sp. nov. is characterized by its bifurcated parameres and a cluster of long black bristles on the inner apex of the gonocoxite. The number of species of Cerapanorpa is raised to 21. An updated key to species of Cerapanorpa is presented.

Phylogeography and demographic history of Gyrodactylus konovalovi (Monogenoidea: Gyrodactylidae), an ectoparasite on the East Asia Amur minnow (Cyprinidae) in Central China

Gyrodactylus konovalovi is an ectoparasite on the Amur minnow (Rhynchocypris lagowskii) that is widely distributed in the cold fresh waters of East Asia. In the present study, the phylogeography and demographic history of G. konovalovi and the distribution of its host in the Qinling Mountains are examined. A total of 79 individual parasites was sequenced for a 528 bp region of the mitochondrial NADH dehydrogenase subunit 5 (ND5) gene, and 25 haplotypes were obtained. The substitution rate (dN/dS) was 0.068 and indicated purifying selection. Haplotype diversity (h) and nucleotide diversity (π) varied widely in the Qinling Mountains. Phylogenetic trees based on Bayesian inference (BI), maximum likelihood (ML), and maximum parsimony (MP) methods and network analysis revealed that all haplotypes were consistently well-supported in three different lineages, indicating a significant geographic distribution pattern. There was a significant positive correlation between genetic differentiation (F st) and geographic distance. The results of mismatch distribution, neutrality test and Bayesian skyline plot analyses showed that whole populations underwent population contraction during the Pleistocene. Based on the molecular clock calibration, the most common ancestor was estimated to have emerged in the middle Pleistocene. Our study suggests for the first time that a clearly phylogeography of G. konovalovi was shaped by geological events and climate fluctuations, such as orogenesis, drainage capture changes, and vicariance, during the Pleistocene in the Qinling Mountains.

Late Pleistocene speciation of three closely related tree peonies endemic to the Qinling-Daba Mountains, a major glacial refugium in Central China

Determining the factors promoting speciation is a major task in ecological and evolutionary research and can be aided by phylogeographic analysis. The Qinling-Daba Mountains (QDM) located in central China form an important geographic barrier between southern subtropical and northern temperate regions, and exhibit complex topography, climatic, and ecological diversity. Surprisingly, few phylogeographic analyses and studies of plant speciation in this region have been conducted. To address this issue, we investigated the genetic divergence and evolutionary histories of three closely related tree peony species (Paeonia qiui, P. jishanensis, and P. rockii) endemic to the QDM. Forty populations of the three tree peony species were genotyped using 22 nuclear simple sequence repeat markers (nSSRs) and three chloroplast DNA sequences to assess genetic structure and phylogenetic relationships, supplemented by morphological characterization and ecological niche modeling (ENM). Morphological and molecular genetic analyses showed the three species to be clearly differentiated from each other. In addition, coalescent analyses using DIYABC conducted on nSSR variation indicated that the species diverged from each other in the late Pleistocene, while ecological niche modeling (ENM) suggested they occupied a larger area during the Last Glacial Maximum (LGM) than at present. The combined genetic evidence from nuclear and chloroplast DNA and the results of ENM indicate that each species persisted through the late Pleistocene in multiple refugia in the Qinling, Daba, and Taihang Mountains with divergence favored by restricted gene flow caused by geographic isolation, ecological divergence, and limited pollen and seed dispersal. Our study contributes to a growing understanding of the origin and population structure of tree peonies and provides insights into the high level of plant endemism present in the Qinling-Daba Mountains of Central China.

Discovery of MicroRNAs from Batrachuperus yenyuanensis Using Deep Sequencing and Prediction of Their Targets

MicroRNAs (miRNAs), a family of ∼22-nucleotide non-coding single-stranded RNA molecules, are considered as key post-transcriptional regulators of gene expression that regulate various biological processes in living organism. Many miRNAs have been identified in animals; however, few have been reported in Hynobiidae species. The present study is aimed to identify a full repertoire of miRNAs in Batrachuperus yenyuanensis (Yenyuan stream salamander), which would significantly increase our knowledge of miRNAs in amphibians. A small RNA library was constructed from B. yenyuanensis and sequenced using deep sequencing. As a result, 1,717,751 clean reads were obtained, representing 356 known and 80 novel miRNAs. Additionally, expression levels of eight randomly selected miRNAs in B. yenyuanensis were confirmed using the stem-loop quantitative real-time reverse transcription PCR. In addition, 13,972 targets were predicted for these identified miRNAs, although the physiological functions of many of these targets remain unknown. Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis suggested that the predicted targets are involved in a variety of physiological regulatory functions in B. yenyuanensis. These results provide useful information for further research on the miRNAs involved in the growth and development of B. yenyuanensis, as well as adaptation of this species to its high-altitude habitats.

Long-term sky islands generate highly divergent lineages of a narrowly distributed stream salamander (Pachyhynobius shangchengensis) in mid-latitude mountains of East Asia

BACKGROUND

Climate oscillation may have a profound effect on species distributions, gene flow patterns and population demography. In response to environmental change, those species restricted to montane habitats experienced expansions and contractions along elevation gradients, which can drive differentiation among sky islands.

RESULTS

The Shangcheng stout salamander (Pachyhynobius shangchengensis) is a cool stream amphibian restricted to high-elevation areas in the Dabie Mountains, East China. In the present study, we used mtDNA genes (Cyt b and ND2) of 193 individuals and 12 nuclear microsatellite loci genotyped on 370 individuals, representing 6 populations (JTX, KHJ, MW, TTZ, BYM and KJY) across the taxon's distribution area, to investigate their genetic variation and evolutionary history of P. shangchengensis. Most populations showed unusually high levels of genetic diversity. Phylogenetic analyses revealed five monophyletic clades with divergence times ranging from 3.96 to 1.4 Mya. Accordingly, significant genetic differentiation was present between these populations. Bayesian skyline plot analyses provided that all populations underwent long-term population expansions since the last inter-glacial (0.13 Mya ~ 0.12 Mya). Msvar analyses found recent signals of population decline for two northern populations (JTX and KHJ) reflecting a strong bottleneck (approximately 15-fold decrease) during the mid-Holocene (about 6000 years ago). Ecological niche modelling has shown a discontinuity in suitable habitats for P. shangchengensis under different historical climatic conditions.

CONCLUSIONS

Our results suggest that the niche conservatism of P. shangchengensis and sky island effects may have led to long-term isolation between populations. In sky island refuges, the mid-latitude Dabie Mountains have provided a long-term stable environment for P. shangchengensis, which has led to the accumulation of genetic diversity and has promoted genetic divergence.