Phylogeographic Pattern of the Assassin Bug Sycanus bifidus Inferred from Mitochondrial Genomes and Nuclear Genes

The assassin bug Sycanus bifidus has a wide distribution across southern China. This study explored its distribution and evolution by analyzing mitochondrial and nuclear ribosomal RNA genes, revealing how Pleistocene climate and geological changes shaped its phylogeography. We identified two main clades, A and B, that diverged in the Middle Pleistocene. Hainan Island's populations form a unique group within Clade A, suggesting that the Qiongzhou Strait served as a dispersal corridor during glaciation. Rising sea levels likely separated the Hainan population afterward. Ecological niche modeling showed that both populations have been viable since the last interglacial period, with demographic analyses indicating possible expansions during the Middle and Late Pleistocene, driven by favorable climates. This study highlights the significant effects of Pleistocene sea-level and climatic changes on the distribution and evolution of S. bifidus in China.

Identification of Risk Areas for Gloydius Snakebites in South Korea

Snakebites can pose a significant threat to human health as the destruction of natural habitats and increased human intrusion into ecosystems result in more frequent encounters with snakes. Mitigation measures for snakebites are particularly crucial for hiking trails where transportation of snakebite victims to medical facilities is challenging due to limited emergency resources and difficult access. This study employed a random forest-based species distribution model approach to investigate the potential habitats of Gloydius spp., specifically Gloydius saxatilis, Gloydius brevicaudus, and Gloydius ussuriensis, in South Korea and to assess the snakebite risk in national parks. Potential habitats of Gloydius spp. were identified and visualized by overlaying binary maps derived from species distribution models (SDMs) of each Gloydius spp. that corresponded to high-risk snakebite areas. In addition, hiking trails with high snakebite risk in the national parks were identified after demonstrating the statistical correlation between the potential habitat distribution of Gloydius spp. and the actual snakebite incidents in major regions of South Korea. The primary environmental variables determining Gloydius spp. habitat were the topographic position index, slope, and the annual average of the maximum and minimum temperatures. The potential habitat of G. saxatilis generally appeared in high-altitude mountainous areas, mostly in the eastern part of the study area. Favorable habitats for G. brevicaudus and G. ussuriensis were predominantly located in mountainous areas throughout the study area, with the exception of some high-altitude mountainous terrain in the east. The number of snakebite incidents per 10,000 people was significantly correlated with the area ratio of Gloydius spp. potential habitat (Spearman's rho = 0.638, p < 0.01). The proportion of snakebite risk areas among national parks in South Korea ranged from 18% to 57%. This study can support practical solutions to prevent injuries and fatalities among hikers due to snakebites by identifying areas with a high risk of snakebite accidents at the hiking-trail level.

Phylogenetic relationships between three Korean pit viper Gloydius (Serpentes: Crotalinae) species using mitochondrial DNA genes

BACKGROUND

Molecular phylogenetic studies of the Asian pit viper genus Gloydius have been widely published in Asia, but Korea population have not been conducted till date.

OBJECTIVE

This study aimed to analyze the phylogenetic relationships of three Gloydius species (G. saxatilis, G. brevicaudus, and G. ussuriensis) from Korea with other Gloydius species, based on Cytochrome b and ND4.

METHODS

We compared 160 samples representing the three species with those of 17 reference species and their phylogenetic status and genetic diversity were analyzed with concatenated sequences of two mitochondrial DNA.

RESULTS

Korean G. brevicaudus and G. saxatilis showed high haplotype diversity and relatively low and moderate nucleotide diversity, respectively. Although G. ussuriensis showed high genetic diversity, it was low in the Baengnyeong Island population. The phylogenetic tree represented two major lineages. One major lineage comprised G. ussuriensis, G. tsushimaensis, G. blomhoffii, and G. brevicaudus. The Chinese G. ussuriensis belonged to the same clade as the Korean G. ussuriensis and was closely related to the Baengnyeong Island population. Moreover, G. tsushimaensis was closely related to G. ussuriensis from southwestern Korean and Jeju Island populations. The other major lineage comprised the remaining 12 species and G. saxatilis. Korean G. saxatilis was closely related to G. saxatilis, G. shedanoensis, and G. intermedius from China.

CONCLUSION

The phylogenetic status of the Korean Gloydius species in comparison with the other Gloydius species was identified. We suggesting the conservation management unit for the Baengnyeong Island population, while the current conservation status of Korean G. saxatilis is suggested to be revised to a higher level.

Ancient volcanos as species pumps: A case study of freshwater amphipods in Northeast Asia

Volcano-tectonic processes have been viewed as primary drivers in the formation of present-day diversity. Volcanos associated with mountain uplifts drive allopatric speciation through vicariance and may impact the surrounding areas like species pump or species attractor. However, the application of these hypotheses to aquatic fauna has rarely been tested explicitly. We tested these hypotheses in the Changbai Mountains (Mts), which are one of the most typical, active volcanic ranges in Northeast (NE) Asia with a long and turbulent geological history. The Gammarus nekkensis species complex of amphipod crustaceans, widely distributed throughout NE Asia with poor dispersal abilities and a long evolutionary history, is a suitable model for testing hypotheses of species pump or species attractor. Phylogenetic and ancestral range reconstructions demonstrated that the studied amphipod originated from the Changbai Mts ~27 Ma and diverged into eastern (Clade I) and western (Clade II) clades, which corresponds well with the initial volcanic eruption of the Changbai Mts in the Late Oligocene. The subsequent diversifications of subclades CI-3, CII-1a and CII-2a were probably driven by second and third eruptions of the Changbai Mts during the Miocene. In particular, the Changbai lineages had spread to the Russian Far East multiple times since the Early Miocene, and widely colonized the region during the Pleistocene. Our discoveries suggest that the ancient volcanos of the Changbai Mts act as species pumps in NE Asia, resulted in burst of diversification around the Changbai Mts and subsequent dispersals into adjacent regions.

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.

Update on Distribution and Conservation Status of Amphibians in the Democratic People's Republic of Korea: Conclusions Based on Field Surveys, Environmental Modelling, Molecular Analyses and Call Properties

Determining the range, status, ecology and behaviour of species from areas where surveys and samplings are uncommon or difficult to conduct is a challenge, such as in the Democratic People's Republic of Korea (DPR Korea). Here, we used genetic samples, field surveys, call recordings, photographic identification and a literature review to estimate the presence, range and status of amphibians in the DPR Korea. From our combined results and based on the IUCN Red List categories and criteria, we were able to estimate the national threat levels for most species. Our results demonstrated the presence of 18 native species and the suspected presence of Karsenia koreana and two Onychodactylus species. We reported the first record for Rana uenoi in the vicinity of Pyongyang using molecular tools and similarly confirmed the presence of Dryophytes japonicus at the same location. Based on distribution and modelling, we can expect the contact zone between species within the Rana and Onychodactylus genera to be located along the Changbai Massif, a mountain range that marks a shift in ecoregions and acts as a barrier to dispersion. The species richness was higher in the lowlands and at lower latitudes, with such areas populated by up to 11 species, while more northern regions were characterised by species richness of about half of that value. The combination of ecological models and known threats resulted in the recommendation of ten species as threatened at the national level following the IUCN Red List categories and criteria. This high number of threatened species was anticipated based on the high threat level to amphibians in bordering nations and globally. While the ecology of species in the DPR Korea is still understudied, we argue that species relying on agricultural wetlands such as rice paddies are not under imminent threat due to the enduring presence of extensive agricultural landscapes with low rates of chemical use and mechanisation. The maintenance of such landscapes is a clear benefit to amphibian species, in contrast to more industrialised agricultural landscapes in neighbouring nations. In comparison, the status of species dependent on forested habitats is unclear and threat levels are likely to be higher because of deforestation, as in neighbouring nations.

Dispersal and mating patterns determine the fate of naturally dispersed populations: evidence from Bombina orientalis

Background

In contrast to the explosive increase of a population following biological invasion, natural dispersal, i.e., when a population disperses from its original range into a new range, is a passive process that is affected by resources, the environment, and other factors. Natural dispersal is also negatively impacted by genetic drift and the founder effect. Although the fates of naturally dispersed populations are unknown, they can adapt evolutionarily over time to the new environment. Can naturally dispersed populations evolve beneficial adaptive strategies to offset these negative effects to maintain their population in a stable state?

Results

The current study addressed this question by focusing on the toad Bombina orientalis, the population of which underwent natural dispersal following the Last Glacial Maximum in Northeast Asia. Population genetic approaches were used to determine the genetic structure, dispersal pattern, and mating system of the population of B. orientalis in northeast China (Northern population). The results showed that this northern population of B. orientalis is a typical naturally dispersed population, in which the stable genetic structure and high level of genetic diversity of the population have been maintained through the long-distance biased dispersal behavior of males and the pattern of promiscuity within the population.

Conclusions

Our findings suggest that naturally dispersed populations can evolve effective adaptive strategies to maintain a stable population. Different species may have different strategies. The relevance of these maintenance mechanisms for naturally dispersed populations provide a new perspective for further understanding the processes of speciation and evolution.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12862-021-01844-3.

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.

Impact of the Miocene orogenesis on Kaloula spp. radiation and implication of local refugia on genetic diversification

The phylogeography of the Kaloula genus in East Asia is still poorly understood. One of the difficulties is the absence of fossils to corroborate molecular dating estimates. Here, we examined the mitochondrial structure of Kaloula spp. in East Asia and focused on the impact of glaciations on the northernmost species: Kaloula borealis. We determined the phylogenetic relationships, molecular dating, and genetic connectivity assessments within the genus from 1211 bp of concatenated mitochondrial 12S and 16S. The relaxed clock analyses reveal the emergence of Kaloula spp. common ancestor in East and Southeast Asia between the Eocene and Oligocene, c. 38.47 Ma (24.69-53.65). The genetic diversification of lineages then increased on the East Asian Mainland during the Lower Miocene, c. 20.10 (8.73-30.65), most likely originating from the vicariance and radiation triggered by the orogeny of the Qinghai-Tibetan Plateau. Later, the dispersal towards the North East Asian Mainland during the Upper Miocene drove the population diversification of K. borealis c. 9.01 Ma (3.66-15.29). Finally, the central mainland population became isolated following orogenesis events and diverged into K. rugifera during the Pliocene, c. 3.06 Ma (0.02-10.90). The combination of population genetic and barrier analyses revealed a significant genetic isolation between populations of Kaloula spp. matching with the massive Qinling-Daba Mountain chain located in south-central China. Finally, we highlight a young divergence within the Eastern Mainland population of K. borealis, possibly attributed to refugia in south eastern China from which populations later expanded.

Phylogeography and ecological niche modeling reveal evolutionary history of Leiolepis ocellata (Squamata, Leiolepidae)

Leiolepis ocellata is a lizard species distributing in topographically diverse habitats in northern Thailand. To explore its evolutionary history, 113 samples of L. ocellata were collected from 11 localities covering its distributional range in northern Thailand, and sequenced for mtDNA fragments (Cyt b and ND2). Pairwise comparisons across sampling localities yielded significant genetic differentiation (F ST and Jost's D) but no clear pattern of isolation by distance could be demonstrated based on the Mantel test. Phylogenetic and network analyses highlighted six haplogroups. Their divergence times were estimated to occur during the Pleistocene, much more recent than major orogenic events affecting northern Thailand. Instead, the results suggested that lineage divergences, of particularly eastern and western haplogroups of the region, coincided with the major rivers in the region (Yom river and Ping river, respectively), indicating vicariance in response to riverine barriers. Furthermore, ecological niche modeling suggested an expansion of suitable habitats of L. ocellata, when LGM-liked conditions. This expansion potentially facilitated their dispersal among adjacent localities leading to lineage diversification and genetic admixture, after the riverine divergence.

The phylogeny, phylogeography, and diversification history of the westernmost Asian cobra (Serpentes: Elapidae: Naja oxiana) in the Trans-Caspian region

We conducted a comprehensive analysis of the phylogenetic, phylogeographic, and demographic relationships of Caspian cobra (Naja oxiana; Eichwald, 1831) populations based on a concatenated dataset of two mtDNA genes (cyt b and ND4) across the species' range in Iran, Afghanistan, and Turkmenistan, along with other members of Asian cobras (i.e., subgenus Naja Laurenti, 1768). Our results robustly supported that the Asiatic Naja are monophyletic, as previously suggested by other studies. Furthermore, N. kaouthia and N. sagittifera were recovered as sister taxa to each other, and in turn sister clades to N. oxiana. Our results also highlighted the existence of a single major evolutionary lineage for populations of N. oxiana in the Trans-Caspian region, suggesting a rapid expansion of this cobra from eastern to western Asia, coupled with a rapid range expansion from east of Iran toward the northeast. However, across the Iranian range of N. oxiana, subdivision of populations was not supported, and thus, a single evolutionary significant unit is proposed for inclusion in future conservation plans in this region.

Regional assemblages shaped by historical and contemporary factors: Evidence from a species-rich insect group

Understanding diversity patterns requires accounting for the roles of both historical and contemporary factors in the assembly of communities. Here, we compared diversity patterns of two moth assemblages sampled from Taihang and Yanshan mountains in Northern China and performed ancestral range reconstructions using the Multi-State Speciation and Extinction model, to track the origins of these patterns. Further, we estimated diversification rates of the two moth assemblages and explored the effects of contemporary ecological factors. From 7,788 specimens we identified 835 species belonging to 23 families, using both DNA barcode analysis and morphology. Moths in Yanshan mountains showed higher species diversity than in Taihang mountains. Ancestral range analysis indicated Yanshan as the origin, with significant historical dispersals from Yanshan to Taihang. Asymmetrical diversification, population expansion, along with frequent and considerable gene flow were detected between communities. Moreover, dispersal limitation or the joint effect of environment filtering and dispersal limitation were inferred as main driving forces shaping current diversity patterns. In summary, we demonstrate that a multiscale (community, population and species level) analysis incorporating both historical and contemporary factors can be useful in delineating factors contributing to community assembly and patterning in diversity.

Diversity of the Bosmina (Cladocera: Bosminidae) in China, revealed by analysis of two genetic markers (mtDNA 16S and a nuclear ITS)

BACKGROUND

China is an important biogeographical zone in which the genetic legacies of the Tertiary and Quaternary periods are abundant, and the contemporary geography environment plays an important role in species distribution. Therefore, many biogeographical studies have focused on the organisms of the region, especially zooplankton, which is essential in the formation of biogeographical principles. Moreover, the generality of endemism also reinforces the need for detailed regional studies of zooplankton. Bosmina, a group of cosmopolitan zooplankton, is difficult to identify by morphology, and no genetic data are available to date to assess this species complex in China. In this study, 48 waterbodies were sampled covering a large geographical and ecological range in China, the goal of this research is to explore the species distribution of Bosmina across China and to reveal the genetic information of this species complex, based on two genetic markers (a mtDNA 16S and a nuclear ITS). The diversity of taxa in the Bosmina across China was investigated using molecular tools for the first time.

RESULTS

Two main species were detected in 35 waterbodies: an endemic east Asia B. fatalis, and the B. longirostris that has a Holarctic distribution. B. fatalis had lower genetic polymorphism and population differentiation than B. longirostris. B. fatalis was preponderant in central and eastern China, whereas B. longirostris was dominated in western China. The third lineage (B. hagmanni) was only detected in a reservoir (CJR) of eastern China (Guangdong province). Bosmina had limited distribution on the Tibetan plateau.

CONCLUSIONS

This study revealed that the biogeography of Bosmina appear to be affected by historical events (Pleistocene glaciations) and contemporary environment (such as altitude, eutrophication and isolated habitat).

Phylogeography of the Assassin Bug Sphedanolestes impressicollis in East Asia Inferred From Mitochondrial and Nuclear Gene Sequences

The assassin bug, Sphedanolestes impressicollis (Hemiptera: Reduviidae), is widely distributed in East Asia. It is an ideal model for evaluating the effects of climatic fluctuation and geographical events on the distribution patterns of East Asian reduviids. Here, we used two mitochondrial genes and one nuclear gene to investigate the phylogeographic pattern of the assassin bug based on comprehensive sampling in China, Japan, South Korea, Vietnam, and Laos. High levels of genetic differentiation were detected among the geographic populations classified into the northern and southern groups. A significant correlation was detected between genetic and geographical distances. The East China Sea land bridge served as a “dispersal corridor” during Pleistocene glaciation. The estimated divergence time indicated that the northern group may have separated from the eastern Chinese populations when the sea level rapidly rose during the “Ryukyu Coral Sea Stage” and the East China Sea land bridge was completely submerged. Demographic history and ecological niche modeling suggested that appropriate climatic conditions may have accounted for the rapid spread across the Korean Peninsula and Japan during the late Pleistocene. Our study underscores the pivotal roles of the Pleistocene sea level changes and climatic fluctuations in determining the distribution patterns of East Asian reduviids.

Evolutionary history and postglacial colonization of an Asian pit viper (Gloydius halys caucasicus) into Transcaucasia revealed by phylogenetic and phylogeographic analyses

It has been generally acknowledged that glacial climates at the time of the Pleistocene altered the patterns of species distributions, prompting latitudinal and altitudinal distribution shifts in several species, including poikilothermic species commonly known for their thermal sensitivity. However, the historical phylogeographic patterns of such species have remained largely unknown. Here, we present the historical biogeographic, phylogenetic, and phylogeographic relationships of the Caucasian pit viper, G. h. caucasicus, based on two mtDNA (cyt b and ND4) and one nDNA (c-mos) genes. This pit viper represents the westernmost member of the Crotalinae subfamily in the Palearctic and occurs in a variety of habitats, from 30 m to 3,000 m above sea level. In Iran, it is distributed on the northern and southern slopes of the Alborz Mountains, rendering it a target for phylogenetic and phylogeographic studies of a terrestrial poikilothermic animal. Our study identified four Iranian lineages of G. h. caucasicus along the northeastern to northwestern slopes of the Alborz Mountains and southern Azerbaijan (Talysh Mountains). Diversification of the Iranian lineages highlights population expansion and subsequent isolation into four plausible refugial areas during the Quaternary paleo-climatic oscillations, confirmed by our molecular dating and historical biogeographic analyses. The results of coalescence-based simulations support the incursion of the species from northeastern Iran to the western end of the Alborz, and then toward Transcaucasia via two directions: northern and southern slopes of the Alborz Mountains. Furthermore, our results clearly implied that G. h. caucasicus should be elevated to species rank and further referred to as G. caucasicus (Nikolsky, 1916).

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.

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.

Multiple refugia from penultimate glaciations in East Asia demonstrated by phylogeography and ecological modelling of an insect pest

BACKGROUND

Refugial populations in Quaternary glaciations are critical to understanding the evolutionary history and climatic interactions of many extant species. Compared with the well-studied areas of Europe and Northern America, refugia of species in eastern Asia remain largely unknown. Here, we investigated the phylogeographic history of a globally important insect pest, the oriental fruit moth Grapholita molesta, in its native range of China.

RESULTS

Genetic structure analyses unveiled three distinct groups and a set of populations with admixture. Approximate Bayesian Computation (ABC) analyses support range expansion of this moth from southwest groups of Yunnan and Sichuan to northern and eastern China. A set of admixed populations was found around these two ancestral groups. This pattern of genetic structure points to two refugia located in the Yunnan region and Sichuan Basin. The split of the two refugia was dated to 329.2 thousand years ago in the penultimate glacial period. One of the lineages was exclusively found around the Sichuan Basin, indicating the formation of endemic populations in this refugium. Ecological niche model analysis suggested a shrinking distribution from the LIG period to the MID period in the Sichuan lineage but a wide and stable distribution in the other lineage.

CONCLUSIONS

Our results for the first time suggest that Yunnan and Sichuan jointly served as two large-scale refugia in eastern Asia in Quaternary glaciations, helping to maintain genetic diversity overall.

Phylogenetic surveys on the newt genus Tylototriton sensu lato (Salamandridae, Caudata) reveal cryptic diversity and novel diversification promoted by historical climatic shifts

Global climatic transitions and Tibetan Plateau uplifts are hypothesized to have profoundly impacted biodiversity in southeastern Asia. To further test the hypotheses related to the impacts of these incidents, we investigated the diversification patterns of the newt genus Tylototriton sensu lato, distributed across the mountain ranges of southeastern Asia. Gene-tree and species-tree analyses of two mitochondrial genes and two nuclear genes revealed five major clades in the genus, and suggested several cryptic species. Dating estimates suggested that the genus originated in the early-to-middle Miocene. Under different species delimitating scenarios, diversification analyses with birth-death likelihood tests indicated that the genus held a higher diversification rate in the late Miocene-to-Pliocene era than that in the Pleistocene. Ancestral area reconstructions indicated that the genus originated from the northern Indochina Peninsula. Accordingly, we hypothesized that the Miocene Climatic Transition triggered the diversification of the genus, and the reinforcement of East Asian monsoons associated with the stepwise uplifts of the Tibetan Plateau promoted the radiation of the genus in southeastern Asia during the Miocene-to-Pliocene period. Quaternary glacial cycles likely had limited effects on speciation events in the genus, but mainly had contributions on their intraspecific differentiations.

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.

Phylogeography of the Rice Spittle Bug (Callitettix versicolor) Implies Two Long-Term Mountain Barriers in South China

South China is a region of remarkable topographic complexity. However, the impact of climate fluctuations in the Pleistocene on the local fauna and especially insects has not been extensively studied. We integrated mitochondrial DNA (mtDNA) and microsatellite data of the rice spittle bug, Callitettix versicolor, to determine the genetic structure, potential biogeographic barriers, and historical demography of this species. The mtDNA data revealed two distinct lineages (Western and Eastern) congruent with the geographically separated western and eastern sub-regions of the Hengduan Mountains. The Eastern lineage was subdivided into two sub-lineages, E1 and E2, congruent with the geographically separated northern and southern sub-regions of the Dabie Mountains. E2 was further subdivided into two sub-groups, E2-1 and E2-2, with a hybrid zone (Guizhou and Hunan Provinces) in which their areas were contiguous. The genetic structures constructed using mtDNA were corroborated by four clusters (G1-G4) of microsatellite data. The populations of each cluster were nearly consistent with a sub-lineage of the mtDNA gene tree (G1-G4 corresponded to the Western, E1, E2-2 and E2-1 lineages, respectively). The divergence time estimated between the Western and Eastern lineages was 1.17 (0.50-2.37) to 0.89 (0.39-1.78) Mya, indicating that the lineages diversified on both geographic and temporal scales. The historical demography of the Eastern lineage showed continuous population growth after the Last Interglacial (LIG) and a stable population during the Last Glacial Maximum (LGM) period. However, the Western lineage remained largely unchanged during the LIG and LGM periods. This suggests that the historical demography of C. versicolor is probably related not only to the paleoclimate of South China, but also to the geological restriction and specific habitat preferences of species.

Influence of geology and human activity on the genetic structure and demography of the Oriental fire-bellied toad (Bombina orientalis)

The Oriental fire-bellied toad (Bombina orientalis) is a commonly used study organism, but knowledge of its evolutionary history is incomplete. We analyze sequence data from four genetic markers (mtDNA genes encoding cytochrome c oxidase subunit I, cytochrome b, and 12S-16S rRNA; nuDNA gene encoding recombination activating gene 2) from 188 individuals across its range in Northeast Asia to elucidate phylogeographic patterns and to identify the historic events that shaped its evolutionary history. Although morphologically similar across its range, B. orientalis exhibits phylogeographic structure, which we infer was shaped by geologic, climatic, and anthropogenic events. Phylogenetic and divergence-dating analyses recover four genetically distinct groups of B. orientalis: Lineage 1-Shandong Province and Beijing (China); Lineage 2-Bukhan Mountain (Korea); Lineage 3-Russia, Northeast China, and northern South Korea; and Lineage 4-South Korea. Lineage 2 was previously unknown. Additionally, we discover an area of secondary contact on the Korean Peninsula, and infer a single dispersal event as the origin of the insular Jeju population. Skyline plots estimate different population histories for the four lineages: Lineages 1 and 2 experienced population decreases, Lineage 3 remained stable, while Lineage 4 experienced a sharp increase during the Holocene. The timing of the population expansion of Lineage 4 coincides with the advent of rice cultivation, which may have facilitated the increase in population size by providing additional breeding habitat.

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.

Phylogeographical analysis of a cold-temperate freshwater fish, the Amur sleeper (Perccottus glenii) in the Amur and Liaohe River basins of Northeast Asia

The Amur sleeper Perccottus glenii (Perciformes, Gobioidei, Odontobutidae) is well known as an invasive fish in the river basins of Eastern and Central Europe, but its genetic background is unavailable across its native habitats in northeast Asia. In this study, we used the mitochondrial cytochrome b gene by sampling 19 populations of P. glenii across its native distributional areas of Liaohe and Amur River basins to explore its evolutionary history. Phylogenetic analyses identified three major clades within P. glenii, among which Clade A and Clade B were co-distributed in the Liaohe and Amur River basins, and Clade C was restricted to the latter. Molecular dating showed that the splits of Clades A, B and C have happened in the late Early-early Middle Pleistocene and the most recent common ancestors of these clades have been presented in the late Middle-early Late Pleistocene. The P. glenii showed very high levels of genetic structure among populations (ΦST = 0.801), probably due to the characters of its life histories with very limited dispersal ability. The admixture of different clades in some populations of P. glenii probably reflects historical secondary contact. These findings indicate that Pleistocene climatic oscillation and river capture were major determinants for genetic variations and evolutionary history of the P. glenii.

The phylogeographical pattern and conservation of the Chinese cobra (Naja atra) across its range based on mitochondrial control region sequences

The vulnerable Chinese cobra (Naja atra) ranges from southeastern China south of the Yangtze River to northern Vietnam and Laos. Large mountain ranges and water bodies may influence the pattern of genetic diversity of this species. We sequenced the mitochondrial DNA control region (1029 bp) using 285 individuals collected from 23 localities across the species' range and obtained 18 sequences unique to Taiwan from GenBank for phylogenetic and population analysis. Two distinct clades were identified, one including haplotypes from the two westernmost localities (Hekou and Miyi) and the other including haplotypes from all sampling sites except Miyi. A strong population structure was found (Φst = 0.76, P<0.0001) with high haplotype diversity (h = 1.00) and low nucleotide diversity (π = 0.0049). The Luoxiao and Nanling Mountains act as historical geographical barriers limiting gene exchange. In the haplotype network there were two “star” clusters. Haplotypes from populations east of the Luoxiao Mountains were represented within one cluster and haplotypes from populations west of the mountain range within the other, with haplotypes from populations south of the Nanling Mountains in between. Lineage sorting between mainland and island populations is incomplete. It remains unknown as to how much adaptive differentiation there is between population groups or within each group. We caution against long-distance transfers within any group, especially when environmental differences are apparent.

Proteomic and biochemical analyses of short-tailed pit viper (Gloydius brevicaudus) venom: age-related variation and composition-activity correlation

We conducted an in-depth analysis of the proteomic and biochemical profiles of the venom of neonate and adult short-tailed pit vipers (Gloydius brevicaudus). Identified proteins were assigned to a few main toxin families. Disintegrin, phospholipase A2 (PLA2), serine proteinase, cysteine-rich secretory protein, C-type lectin-like protein, l-amino acid oxidase and snake venom metalloproteinase (SVMP) were detected in both venoms, while 5'-nucleotidase was detected only in the adult venom. SVMP was the predominant protein family in both venoms (neonate: 65.7%; adult: 64.4%), followed by PLA2 (neonate: 13.4%; adult: 25.0%). Antivenomic analysis revealed that commercial G. brevicaudus antivenom almost neutralized the chromatographic peaks with medium and high molecular masses in both venoms, but did not completely recognize peaks with low molecular mass. Toxicological and enzymatic activities show remarkable age-related variation in G. brevicaudus venom, probably resulting from variation in venom composition. Our data demonstrate age-related variation across venomics, antivenomics and biochemical profiles of G. brevicaudus venom, and have implications for the management of G. brevicaudus bites, including improving antivenom preparation by combining both venoms.

BIOLOGICAL SIGNIFICANCE

This study investigates the composition and biochemical activity of neonate and adult Gloydius brevicaudus venoms. We found remarkable age-related variation in venom biological activity, likely the result of variation in venom composition. Antivenomics analysis was used to explore difference in neonate and adult G. brevicaudus venoms. Our findings have implications for the diagnosis and clinical management of G. brevicaudus bites, and the design of venom mixtures that will increase the efficacy of commercial antivenom. This article is part of a Special Issue entitled: Proteomics of non-model organisms.

Phylogeography of the Rickett's big-footed bat, Myotis pilosus (Chiroptera: Vespertilionidae): a novel pattern of genetic structure of bats in China

Background

China is characterized by complex topographic structure and dramatic palaeoclimatic changes, making species biogeography studies particularly interesting. Previous researchers have also demonstrated multiple species experienced complex population histories, meanwhile multiple shelters existed in Chinese mainland. Despite this, species phylogeography is still largely unexplored. In the present study, we used a combination of microsatellites and mitochondrial DNA (mtDNA) to investigate the phylogeography of the east Asian fish-eating bat (Myotis pilosus).

Results

Phylogenetic analyses showed that M. pilosus comprised three main lineages: A, B and C, which corresponded to distinct geographic populations of the Yangtze Plain (YTP), Sichuan Basin (SCB) and North and South of China (NSC), respectively. The most recent common ancestor of M. pilosus was dated as 0.25 million years before present (BP). Population expansion events were inferred for populations of Clade C, North China Plain region, Clade B and YunGui Plateau region at 38,700, 15,900, 4,520 and 4,520 years BP, respectively. Conflicting results were obtained from mtDNA and microsatellite analyses; strong population genetic structure was obtained from mtDNA data but not microsatellite data. The microsatellite data indicated that genetic subdivision fits an isolation-by-distance (IBD) model, but the mtDNA data failed to support this model.

Conclusions

Our results suggested that Pleistocene climatic oscillations might have had a profound influence on the demographic history of M. pilosus. Spatial genetic structures of maternal lineages that are different from those observed in other sympatric bats species may be as a result of interactions among special population history and local environmental factors. There are at least three possible refugia for M. pilosus during glacial episodes. Apparently contradictory genetic structure patterns of mtDNA and microsatellite could be explained by male-mediated gene flow among populations. This study also provides insights on the necessity of conservation of M. pilosus populations to conserve this genetic biodiversity, especially in the areas of YTP, SCB and NSC regions.

Phylogeography of Silver Pheasant (Lophura nycthemera L.) across China: aggregate effects of refugia, introgression and riverine barriers

The role of Pleistocene glacial cycles in forming the contemporary genetic structure of organisms has been well studied in China with a particular focus on the Tibetan Plateau. However, China has a complex topography and diversity of local climates, and how glacial cycles may have shaped the subtropical and tropical biota of the region remains mostly unaddressed. To investigate the factors that affected the phylogeography and population history of a widely distributed and nondeciduous forest species, we analysed morphological characters, mitochondrial DNA sequences and nuclear microsatellite loci in the Silver Pheasant (Lophura nycthemera). In a pattern generally consistent with phenotypic clusters, but not nominal subspecies, deeply divergent mitochondrial lineages restricted to different geographic regions were detected. Coalescent simulations indicated that the time of main divergence events corresponded to major glacial periods in the Pleistocene and gene flow was only partially lowered by drainage barriers between some populations. Intraspecific cytonuclear discordance was revealed in mitochondrial lineages from Hainan Island and the Sichuan Basin with evidence of nuclear gene flow from neighbouring populations into the latter. Unexpectedly, hybridization was revealed in Yingjiang between the Silver Pheasant and Kalij Pheasant (Lophura leucomelanos) with wide genetic introgression at both the mtDNA and nuclear levels. Our results highlight a novel phylogeographic pattern in a subtropical area generated from the combined effects of climate oscillation, partial drainage barriers and interspecific hybridization. Cytonuclear discordance combined with morphological differentiation implies that complex historical factors shaped the divergence process in this biodiversity hot spot area of southern China.

Phylogeographic patterns of mtDNA variation revealed multiple glacial refugia for the frog species Feirana taihangnica endemic to the Qinling Mountains

Diversification patterns and demography of montane species are affected by Pleistocene climate fluctuations. Empirical cases from the Qinling Mountains (QM) region, which is a major biogeographic divider of East Asia, are few. We used DNA sequence data of the complete mitochondrial ND2 gene to detect effects of the Pleistocene glaciations on phylogeographic profiles of a frog species, Feirana taihangnica, which is endemic to the QM. Four distinct lineages consisting of seven sublineages were revealed. The strongest signal of biogeographical structure (F(ct) = 0.971, P < 0.01) was found when populations were grouped according to these seven sublineages. One narrow secondary contact zone was detected in the middle QM between the lineage from middle QM and the lineage from eastern QM. Coalescent simulations indicated that this species colonized the QM region by a stepping-stone model. Divergences among lineages had likely been influenced by the uplift of the Tibetan Plateau during the late Miocene-to-late Pleistocene, as well as by the Pleistocene climatic cycles. Coalescent simulations also suggested that F. taihangnica populations have persisted through the Pleistocene glacial periods in multiple refugia across the QM region. Demographic analyses indicated that all lineages, except the lineage in the Funiu Mountains, have been experienced postglacial expansion of population size and distribution range. In conclusion, Pleistocene climate fluctuations and tectonic changes during the late Miocene-late Pleistocene have profoundly influenced the phylogeography and historical demography of F. taihangnica.

Genealogy and Demographic History of a Widespread Amphibian throughout Indochina

Relatively little is known about spatial patterns of cryptic diversity in tropical species and the processes that generate them. Few studies examine the geographic distribution of genetic lineages in Southeast Asia, an area hypothesized to harbor substantial cryptic diversity. We investigated the evolutionary history of Asian tree frogs of the Polypedates leucomystax complex (n = 172) based on 1800 bp of the mtDNA genes ND1 and cytochrome b and tested hypotheses pertaining to climate, geology, and dispersal patterns. Analyses revealed substantial genetic diversity and lineage divergence throughout the region with evidence for widespread sympatric lineages and a general north versus south clustering. Relaxed molecular clock analysis and tests for demographic expansion identified an initial cladogenesis during the Miocene with subsequent Plio-Pleistocene diversification, with the former corresponding to periods of increased aridity and the onset of monsoonal weather systems. Rates of diversification were relatively constant until the Early Pleistocene when rates increased exponentially. We found equivocal evidence for isolation-by-distance and a potential role of some landscape features as partial barriers to dispersal. Finally, our analyses showed that divergence between insular and mainland populations occurred before Homo sapiens colonized Southeast Asia, suggesting that historical human-mediated dispersal did not drive insular diversification. Our results suggested that demographic expansion in the Late Pleistocene resulted in widespread sympatric lineages in the P. leucomystax complex throughout southern China and Indochina and further clarified the evolutionary history of lineages within P. leucomystax.

Phylogeography of the Mekong mud snake (Enhydris subtaeniata): the biogeographic importance of dynamic river drainages and fluctuating sea levels for semiaquatic taxa in Indochina

During the Cenozoic, Southeast Asia was profoundly affected by plate tectonic events, dynamic river systems, fluctuating sea levels, shifting coastlines, and climatic variation, which have influenced the ecological and evolutionary trajectories of the Southeast Asian flora and fauna. We examined the role of these paleogeographic factors on shaping phylogeographic patterns focusing on a species of semiaquatic snake, Enhydris subtaeniata (Serpentes: Homalopsidae) using sequence data from three mitochondrial fragments (cytochrome b, ND4, and ATPase-2785 bp). We sampled E. subtaeniata from seven locations in three river drainage basins that encompassed most of this species' range. Genetic diversities were typically low within locations but high across locations. Moreover, each location had a unique suite of haplotypes not shared among locations, and pairwise φ(ST) values (0.713-0.998) were highly significant between all location pairs. Relationships among phylogroups were well resolved and analysis of molecular variance (AMOVA) revealed strong geographical partitioning of genetic variance among the three river drainage basins surveyed. The genetic differences observed among the populations of E. subtaeniata were likely shaped by the Quaternary landscapes of Indochina and the Sunda Shelf. Historically, the middle and lower Mekong consisted of strongly dissected river valleys separated by low mountain ranges and much of the Sunda Shelf consisted of lowland river valleys that served to connect faunas associated with major regional rivers. It is thus likely that the contemporary genetic patterns observed among populations of E. subtaeniata are the result of their histories in a complex terrain that created abundant opportunities for genetic isolation and divergence yet also provided lowland connections across now drowned river valleys.

Postglacial colonization of the Qinling Mountains: phylogeography of the swelled vent frog (Feirana quadranus)

BACKGROUND

The influence of Pleistocene climatic fluctuations on intraspecific diversification in the Qinling-Daba Mountains of East Asia remains poorly investigated. We tested hypotheses concerning refugia during the last glacial maximum (LGM) in this region by examining the phylogeography of the swelled vent frog (Feirana quadranus; Dicroglossidae, Anura, Amphibia).

METHODOLOGY/PRINCIPAL FINDINGS

We obtained complete mitochondrial ND2 gene sequences of 224 individuals from 34 populations of Feirana quadranus for phylogeographic analyses. Additionally, we obtained nuclear tyrosinase gene sequences of 68 F. quadranus, one F. kangxianensis and three F. taihangnica samples to test for mitochondrial introgression among them. Phylogenetic analyses based on all genes revealed no introgression among them. Phylogenetic analyses based on ND2 datasets revealed that F. quadranus was comprised of six lineages which were separated by deep valleys; the sole exception is that the Main Qinling and Micang-Western Qinling lineages overlap in distribution. Analyses of population structure indicated restricted gene flow among lineages. Coalescent simulations and divergence dating indicated that the basal diversification within F. quadranus may be associated with the dramatic uplifts of the Tibetan Plateau during the Pliocene. Coalescent simulations indicated that Wuling, Daba, and Western Qinling-Micang-Longmen Mountains were refugia for F. quadranus during the LGM. Demographic analyses indicated that the Daba lineage experienced population size increase prior to the LGM but the Main Qinling and the Micang-Western Qinling lineages expanded in population size and range after the LGM, and the other lineages almost have stable population size or slight slow population size decline.

CONCLUSIONS/SIGNIFICANCE

The Qinling-Daba Mountains hosted three refugia for F. quadranus during the LGM. Populations that originated in the Daba Mountains colonized the Main Qinling Mountains after the LGM. Recent sharp expansion of the Micang-Western Qinling and Main Qinling lineages probably contribute to their present-day secondary contact.