Phylogeography of a widespread terrestrial vertebrate in a barely-studied Palearctic region: green toads (Bufo viridis subgroup) indicate glacial refugia in Eastern Central Asia

The Fast Evolution of the Stenobothrini Grasshoppers (Orthoptera, Acrididae, and Gomphocerinae) Revealed by an Analysis of the Control Region of mtDNA, with an Emphasis on the Stenobothrus eurasius Group

The two cryptic grasshopper species of the genus Stenobothrus, S. eurasius and S. hyalosuperficies, demonstrate different acoustic behavior despite a strong similarity in morphology. A hybridization between these species is possible in the contact zone; however, there are little molecular data about the relationships of these species. The analysis of the mtDNA control region (CR) reveals that haplotypes of S. hyalosuperficies have more in common with the more distant Stenobothrus species than with the closely related S. eurasius. In the contact zone, S. eurasius has mt-haplotypes shared with S. hyalosuperficies, which might indicate an introgression of mtDNA from S. hyalosuperficies to the S. eurasius gene pool. We also analyze the structure and evolutionary rate of the mtDNA CR for the Stenobothrus genus and estimate the time of divergence of the species within the genus. The phylogenetic tree of the tribe Stenobothrini reconstructed with either the CR or COI gave the same four groups. The phylogenetic tree of the Stenobothrus genus has a star-like topology with each mtDNA haplotype found in any analyzed species, except S. eurasius, which forms a separate branch. The maximum degree of incomplete lineage sorting can demonstrate either ancestral polymorphism or introgression.

The Silk roads: phylogeography of Central Asian dice snakes (Serpentes: Natricidae) shaped by rivers in deserts and mountain valleys

Influenced by rapid changes in climate and landscape features since the Miocene, widely distributed species provide suitable models to study the environmental impact on their evolution and current genetic diversity. The dice snake Natrix tessellata, widely distributed in the Western Palearctic is one such species. We aimed to resolve a detailed phylogeography of N. tessellata with a focus on the Central Asian clade with 4 and the Anatolia clade with 3 mitochondrial lineages, trace their origin, and correlate the environmental changes that affected their distribution through time. The expected time of divergence of both clades began at 3.7 Mya in the Pliocene, reaching lineage differentiation approximately 1 million years later. The genetic diversity in both clades is rich, suggesting different ancestral areas, glacial refugia, demographic changes, and colonization routes. The Caspian lineage is the most widespread lineage in Central Asia, distributed around the Caspian Sea and reaching the foothills of the Hindu Kush Mountains in Afghanistan, and Eastern European lowlands in the west. Its distribution is limited by deserts, mountains, and cold steppe environments. Similarly, Kazakhstan and Uzbekistan lineages followed the Amu Darya and the Syr Darya water systems in Central Asia, with ranges delimited by the large Kyzylkum and Karakum deserts. On the western side, there are several lineages within the Anatolia clade that converged in the central part of the peninsula with 2 being endemic to Western Asia. The distribution of both main clades was affected by expansion from their Pleistocene glacial refugia around the Caspian Sea and in the valleys of Central Asia as well as by environmental changes, mostly through aridification.

Species distribution modelling sheds light on the widespread distribution of Sphingonotus (Sphingonotus) rubescens (Orthoptera: Acrididae: Oedipodinae)

Sphingonotus (Sphingonotus) rubescens (Walker, 1870) occurs from the Atlantic Islands to central and southern Asia, but its ecological preferences and the potential drivers that shaped its currently extensive distribution remain unknown. We review the known occurrence data for S. rubescens and reconstruct its current and palaeoclimatic niche (Last Glacial Maximum and mid-Holocene) using species distribution modelling (SDM). We examine how climatic suitability might have created potential past migratory pathways shaping its current distribution. Moreover, we test the utility of SDM approaches to flag dubious records detected in the assembled dataset. The results reveal new distributional records from four countries. Climatic models indicate high levels of overlap between current and palaeoclimatic models, with stability of large suitable areas through time. Furthermore, we observe that suitability values associated with dubious records are lower than the distribution of suitability values within the known distribution of S. rubescens. Climatic stability of suitable areas through time for S. rubescens might have aided the expansion and maintenance of its current wide distribution. Furthermore, our results support previous studies indicating the usefulness of SDM tools for the detection of doubious occurrences.

Environmental and geological changes in the Tarim Basin promoted the phylogeographic formation of Phrynocephalus forsythii (Squamata: Agamidae)

Environmental factors can promote genetic divergence among populations. The Tarim Basin has experienced a series of environmental and geological changes since the late Pliocene. Phrynocephalus forsythii (Agamidae) has a continuous and circular distribution around the Tarim Basin, and the evolutionary history of this species remains unclear. In the present study, mitochondrial DNA and single nucleotide polymorphism sequences were obtained from 195 P. forsythii samples in the Tarim Basin to examine the phylogeographic structure and evolutionary history of this species. All populations of P. forsythii formed three distinct clusters: the Minfeng, low-elevation, and high-elevation groups. The Minfeng group was the first to separate from all other groups at 4.26 Ma after the original desert environment emerged in the Tarim Basin. The uplift of the Kunlun Mountains created a new high-altitude environment along the southwestern edge of the basin and promoted the divergence of low-elevation and high-elevation groups at 3.67 Ma. Subsequently, the gradual development of deserts and rivers in the Tarim Basin accelerated the dispersal of P. forsythii, eventually forming the current distribution pattern. Our results indicate that ecological separation caused by environmental and geological changes in the Tarim Basin was the main cause for the genetic divergence of P. forsythii.

Impacts of late Quaternary environmental change on the long-tailed ground squirrel (Urocitellus undulatus) in Mongolia

Impacts of Quaternary environmental changes on mammal faunas of central Asia remain poorly understood due to a lack of geographically comprehensive phylogeographic sampling for most species. To help address this knowledge gap, we conducted the most extensive molecular analysis to date of the long-tailed ground squirrel (Urocitellus undulatus Pallas 1778) in Mongolia, a country that comprises the southern core of this species' range. Drawing on material from recent collaborative field expeditions, we genotyped 128 individuals at 2 mitochondrial genes (cytochrome b and cytochrome oxidase I; 1 797 bp total). Phylogenetic inference supports the existence of two deeply divergent infraspecific lineages (corresponding to subspecies U. u. undulatus and U. u. eversmanni), a result in agreement with previous molecular investigations but discordant with patterns of range-wide craniometric and external phenotypic variation. In the widespread westerneversmanni lineage, we recovered geographically-associated clades from the: (a) Khangai, (b) Mongolian Altai, and (c) Govi Altai mountain ranges. Phylogeographic structure in U. u. eversmanni is consistent with an isolation-by-distance model; however, genetic distances are significantly lower than among subspecies, and intra-clade relationships are largely unresolved. The latter patterns, as well as the relatively higher nucleotide polymorphism of populations from the Great Lakes Depression of northwestern Mongolia, suggest a history of range shifts into these lowland areas in response to Pleistocene glaciation and environmental change, followed by upslope movements and mitochondrial lineage sorting with Holocene aridification. Our study illuminates possible historical mechanisms responsible for U. undulatus genetic structure and contributes to a framework for ongoing exploration of mammalian response to past and present climate change in central Asia.

Evolutionary melting pots: a biodiversity hotspot shaped by ring diversifications around the Black Sea in the Eastern tree frog (Hyla orientalis)

Hotspots of intraspecific genetic diversity, which are of primary importance for the conservation of species, have been associated with glacial refugia, that is areas where species survived the Quaternary climatic oscillations. However, the proximate mechanisms generating these hotspots remain an open issue. Hotspots may reflect the long-term persistence of large refugial populations; alternatively, they may result from allopatric differentiation between small and isolated populations, that later admixed. Here, we test these two scenarios in a widely distributed species of tree frog, Hyla orientalis, which inhabits Asia Minor and southeastern Europe. We apply a fine-scale phylogeographic survey, combining fast-evolving mitochondrial and nuclear markers, with a dense sampling throughout the range, as well as ecological niche modelling, to understand what shaped the genetic variation of this species. We documented an important diversity centre around the Black Sea, composed of multiple allopatric and/or parapatric diversifications, likely driven by a combination of Pleistocene climatic fluctuations and complex regional topography. Remarkably, this diversification forms a ring around the Black Sea, from the Caucasus through Anatolia and eastern Europe, with terminal forms coming into contact and partially admixing in Crimea. Our results support the view that glacial refugia generate rather than host genetic diversity and can also function as evolutionary melting pots of biodiversity. Moreover, we report a new case of ring diversification, triggered by a large, yet cohesive dispersal barrier, a very rare situation in nature. Finally, we emphasize the Black Sea region as an important centre of intraspecific diversity in the Palearctic with implications for conservation.

Phylogeographic Structure of a Tethyan Relict Capparis spinosa (Capparaceae) Traces Pleistocene Geologic and Climatic Changes in the Western Himalayas, Tianshan Mountains, and Adjacent Desert Regions

Complex geological movements more or less affected or changed floristic structures, while the alternation of glacials and interglacials is presumed to have further shaped the present discontinuous genetic pattern of temperate plants. Here we consider Capparis spinosa, a xeromorphic Tethyan relict, to discuss its divergence pattern and explore how it responded in a stepwise fashion to Pleistocene geologic and climatic changes. 267 individuals from 31 populations were sampled and 24 haplotypes were identified, based on three cpDNA fragments (trnL-trnF, rps12-rpl20, and ndhF). SAMOVA clustered the 31 populations into 5 major clades. AMOVA suggests that gene flow between them might be restricted by vicariance. Molecular clock dating indicates that intraspecific divergence began in early Pleistocene, consistent with a time of intense uplift of the Himalaya and Tianshan Mountains, and intensified in mid-Pleistocene. Species distribution modeling suggests range reduction in the high mountains during the Last Glacial Maximum (LGM) as a result of cold climates when glacier advanced, while gorges at midelevations in Tianshan appear to have served as refugia. Populations of low-altitude desert regions, on the other hand, probably experienced only marginal impacts from glaciation, according to the high levels of genetic diversity.

Low phylogeographic structure of Rhinella arunco (Anura: Bufonidae), an endemic amphibian from the Chilean Mediterranean hotspot

Background

The Mediterranean zone of central Chile (30° to 38°S) is one of the 25 diversity hotspots in the world. However, there are few phylogeographic studies which identify the factors that have influenced population diversification in the fauna of this area. In this study, we investigated the phylogeographic structure of Rhinella arunco, an anuran endemic to Mediterranean Chile, using 160 individuals from 23 localities representative of its entire distribution (32° to 38°S).

Results

The haplotype network revealed four haplogroups, three of which overlap geographically and only one of which has an exclusive geographic distribution. An analysis of molecular variance indicated that neither watershed limits nor the main rivers in the current distribution of this species have been important geographic barriers against the dispersion of individuals. Finally, the Geneland analysis showed three population units, one of which concurs with one of the haplogroups found in the haplotype network. Together, these analyses indicated a low level of phylogeographic structure for this species. On the other hand, the highest levels of intrapopulational genetic variation were concentrated in the central part of the distribution (33° to 34°S), which may indicate an effect of Pleistocene glaciations on the genetic diversity of the populations in the extreme south of its range.

Conclusions

The low phylogeographic structure observed in R. arunco is a rarely documented pattern for amphibians and contrasts with the phylogeographic studies of other vertebrates which inhabit the same zone. This result may be attributed to a series of attributes of bufonids, related mainly to water retention and their reproductive biology, which have allowed them to disperse and colonize an enormous variety of environments.

Retreating or standing: responses of forest species and steppe species to climate change in arid Eastern Central Asia

Background

The temperature in arid Eastern Central Asia is projected to increase in the future, accompanied by increased variability of precipitation. To investigate the impacts of climate change on plant species in this area, we selected two widespread species as candidates, Clematis sibirica and C. songorica, from montane coniferous forest and arid steppe habitats respectively.

Methodology/Principal Findings

We employed a combined approach of molecular phylogeography and species distribution modelling (SDM) to predict the future responses of these two species to climate change, utilizing evidence of responses from the past. Genetic data for C. sibirica shows a significant phylogeographical signal (N_ST > F_ST, P<0.05) and demographic contraction during the glacial-interglacial cycles in the Pleistocene. This forest species would likely experience range reduction, though without genetic loss, in the face of future climate change. In contrast, SDMs predict that C. songorica, a steppe species, should maintain a consistently stable potential distribution under the Last Glacial Maximum (LGM) and the future climatic conditions referring to its existing potential distribution. Molecular results indicate that the presence of significant phylogeographical signal in this steppe species is rejected and this species contains a high level of genetic differentiation among populations in cpDNA, likely benefiting from stable habitats over a lengthy time period.

Conclusions/Significance

Evidence from the molecular phylogeography of these two species, the forest species is more sensitive to past climate changes than the steppe species. SDMs predict that the forest species will face the challenge of potential range contraction in the future more than the steppe species. This provides a perspective on ecological management in arid Eastern Central Asia, indicating that increased attention should be paid to montane forest species, due to their high sensitivity to disturbance.

Biogeographical origin and speciation of the Anthocoris nemorum group

The Anthocoris nemorum group belongs to the Anthocoridae (Hemiptera), and is an important group of predators of agricultural pests. A phylogeny was constructed in conjunction with dispersal-vicariance analysis of the Anthocoris nemorum group species in order to discern the relationships between the phylogeographical structuring of A. nemorum group species, and the effects of the Qinghai-Tibet plateau uplift. The divergence times were estimated using Bayesian inference as implemented in BEAST. A portion of the mitochondrial COI gene (1406 bp) and 16S rDNA (932 bp) were chosen as molecular markers to reconstruct evolutionary relationships among 10 species in this group. The combined approach, based on phylogeny, estimation of node dates, and dispersal-vicariance analyses, indicated that the phylogeographical structuring of A. nemorum has been primarily shaped by the two main periods of the Qinghai-Tibet plateau uplift. The DIVA optimal reconstructions suggest that A. nemorum diversified from the Miocene to the Pleistocene from a center of origin in the Hengduan Mountains. The rapid uplift of Mountain ranges associated with the uplift of the entire Qinghai-Tibet plateau may have promoted rapid divergence in the A. nemorum group. Vicariance and dispersal were both essential in shaping the present distribution patterns of A. nemorum.

Comparative population genetics and phylogeography of two lacertid lizards (Eremias argus and E. brenchleyi) from China

Eremias argus and Eremias brenchleyi are lacertid lizards that are sympatric throughout the distribution of E. brenchleyi. We sequenced partial mitochondrial DNA from cytochrome (cyt) b gene for 106 individuals of E. argus from nine localities, and for 45 individuals of E. brenchleyi from five localities, in central and northern parts of North China. We determined 53 cyt b haplotypes from the E. argus samples, and 27 cyt b haplotypes from the E. brenchleyi samples. Only E. brenchleyi had followed a stepping-stone model of dispersal. Partitioned Bayesian phylogenetic analysis reveals that E. argus and E. brenchleyi are reciprocally monophyletic, and the divergence time between the two species was dated to about 4.1±1.2 million years ago. Geographical structuring of haplotypes is more significant in E. brenchleyi than in E. argus. Haplotypes of E. brenchleyi could be divided into four groups by the Yellow River and Taihang Mountains. Within-population genetic diversity indices are correlated neither with latitude nor with longitude. We calculated significant among-population structure for both species (E. argus: Φ(ST)=0.608, P<0.001; E. brenchleyi: Φ(ST)=0.925, P<0.001). Eremias brenchleyi has four independent management units, while E. argus has a more homogeneous genetic structure across its range. Our data show that: (1) the pattern seen in North American and European species that southern populations have higher genetic diversity as consequence of post-glaciation dispersal is absent in the two Chinese lizards; (2) the Yellow River and Taihang Mountains may have acted as important barriers to gene flow only in E. brenchleyi; and (3) genetic structure differs between the two lizards that differ in habitat preference and dispersal ability.