Phylogeographic and demographic effects of Pleistocene climatic fluctuations in a montane salamander, Plethodon fourchensis

Comparative Study of Potential Habitats for Simulium qinghaiense (Diptera: Simuliidae) in the Huangshui River Basin, Qinghai-Tibet Plateau: An Analysis Using Four Ecological Niche Models and Optimized Approaches

The Huangshui River, a vital tributary in the upper reaches of the Yellow River within the eastern Qinghai-Tibet Plateau, is home to the endemic black fly species S. qinghaiense. In this study, we conducted a systematic survey of the distribution of the species in the Huangshui River basin, revealing its predominant presence along the river's main stem. Based on four ecological niche models-MaxEnt with parameter optimization; GARP; BIOCLIM; and DOMAIN-we conduct a comparative analysis; evaluating the accuracy of AUC and Kappa values. Our findings indicate that optimizing parameters significantly improves the MaxEnt model's predictive accuracy by reducing complexity and overfitting. Furthermore, all four models exhibit higher accuracy compared to a random model, with MaxEnt demonstrating the highest AUC and Kappa values (0.9756 and 0.8118, respectively), showcasing significant superiority over the other models (p < 0.05). Evaluation of predictions from the four models elucidates that potential areas of S. qinghaiense in the Huangshui River basin are primarily concentrated in the central and southern areas, with precipitation exerting a predominant influence. Building upon these results, we utilized the MaxEnt model to forecast changes in suitable areas and distribution centers during the Last Interglacial (LIG), Mid-Holocene (MH), and future periods under three climate scenarios. The results indicate significantly smaller suitable areas during LIG and MH compared to the present, with the center of distribution shifting southeastward from the Qilian Mountains to the central part of the basin. In the future, suitable areas under different climate scenarios are expected to contract, with the center of distribution shifting southeastward. These findings provide important theoretical references for monitoring, early warning, and control measures for S. qinghaiense in the region, contributing to ecological health assessment.

Genomic analyses reveal three phylogenetic species and their evolutionary histories in the big-headed turtle

The critically endangered big-headed turtle (Platysternon megacephalum) is currently classified into three subspecies. However, the classification is still controversial and their evolutionary histories are still unclear. Here, multiple genetic analyses consistently revealed three phylogenetic groups with substantial genetic divergences and distinct demographic histories, suggesting three phylogenetic species (P. megacephalum, P. peguense, and Baise clade). Phylogeographical analyses revealed that the Red River plains and Guangxi basins are largely coincident with the boundaries between the three phylogenetic species, highlighting the key role of lowland areas in driving speciation in the big-headed turtle. The Baise clade is characterized by high-linkage disequilibrium but the lowest effective population size, indicating that the cryptic phylogenetic species is more vulnerable to human activities and environmental disturbance, and urgently needs more protection. Our findings provide fundamental insights into the taxonomy and scientific conservation of the family Platysternidae.

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.

Biodiversity of Parasitic Hymenoptera Across Sky Islands of Arkansas, United States

Sky islands are often areas of endemism and high species diversity around the world. In central North America, the Ozark and Ouachita Highlands contain numerous modestly high elevational islands that are slightly cooler and wetter than surrounding lowlands and these areas are known to contain many unique species. We studied parasitic Hymenoptera biodiversity on three sky islands in Arkansas, both in canopy and understory strata of mature post oak forests. We found that morphospecies diversity was high, and community structure varied both within (canopy versus understory) and between sky islands. Each mountaintop showed a high level of distinctiveness with few morphospecies found in multiple locations. The high intersky island variation indicates that individual sky islands in the region harbor unique Hymenoptera assemblages, a pattern that likely extends to the broader insect community. Our study suggests that these highland areas are important regions of North American biodiversity and that they should be evaluated individually for conservation efforts in order to preserve their distinctive community structure.

Hybridization, Behavioral Patterns, and Pre- and Postzygotic Isolation Between Two Recently Derived Species of Walnut-Infesting Rhagoletis Fruit Flies in the Highlands of Mexico

Episodes of isolation and secondary contact among populations of insects of Nearctic origin during Pleistocene glacial/postglacial climatic cycles had a strong evolutionary influence on the diversity of flies in the genus Rhagoletis in mountainous areas of Mexico. As a series of experiments undertaken to gather support for phylogenetic hypotheses on the origin of three walnut-infesting species in the suavis group, we examined pre- and postzygotic isolation between Rhagoletis completa Cresson, 1929 and R. ramosae Hernández-Ortiz, 1985. Despite morphological, biological, and behavioral differences, these two species were found to be capable of hybridization. Mating experiments in large enclosures revealed asymmetric sexual isolation. There were notable differences in male sexual behavior. While R. ramosae males mated exclusively on host fruit, R. completa males used fruit and alternative mating locations. During fruit-guarding and male-male contests, R. completa and R. ramosae males adopted markedly different wing postures. R. completa females were more reluctant to copulate with heterospecific males than R. ramosae females. During no choice crosses in small enclosures, there was a reduction of egg hatch for the hybrid cross of R. completa males × R. ramosae females. Our results and previous studies on reproductive isolation between other species pairs in the suavis group support a clade in which R. ramosae, R. zoqui Bush, 1966, and R. completa are close relatives all still capable of hybridizing.

Revealing the biodiversity of Chilean birds through the COI barcode approach

The mitochondrial cytochrome c oxidase subunit I (COI) gene is an effective molecular tool for the estimation of genetic variation and the identification of bird species. This molecular marker is used to differentiate among Chilean bird species by analyzing barcodes for 76 species (197 individuals), comprising 28 species with no previous barcode data and 48 species with sequences retrieved from the BOLD and GenBank databases. The DNA barcodes correctly identified 94.7% of the species analyzed (72 of 76 species). Mean intraspecific K2P distance was 0.3% (range 0-8.7%). Within the intraspecific divergence range, three species, Phrygilus gayi, Sephanoides sephanoides and Curaeus curaeus, showed relatively high intraspecific divergence (1.5-8.7%), possibly due to the presence of a species complex or geographic isolation of sub-populations. Mean interspecific K2P distance was 24.7% (range 1.3-43.5%). Consequently, the intraspecific K2P distance showed limited overlap with interspecific K2P distance. The mean intraspecific divergence in our study was similar to that found in temperate regions of South America (0.24%). However, it was approximately one order of magnitude lower than values reported for bird species in tropical regions of northern South America (1.8-2.13%). This result suggests that bird species from Chile show low levels of genetic structure and divergence. The small overlap between intra- and inter-specific distances implies that COI barcodes could be used as an effective tool to identify nearly all the Chilean bird species analyzed.

Last Glacial Maximum led to community-wide population expansion in a montane songbird radiation in highland Papua New Guinea

BACKGROUND

Quaternary climate fluctuations are an engine of biotic diversification. Global cooling cycles, such as the Last Glacial Maximum (LGM), are known to have fragmented the ranges of higher-latitude fauna and flora into smaller refugia, dramatically reducing species ranges. However, relatively less is known about the effects of cooling cycles on tropical biota.

RESULTS

We analyzed thousands of genome-wide DNA markers across an assemblage of three closely related understorey-inhabiting scrubwrens (Sericornis and Aethomyias; Aves) from montane forest along an elevational gradient on Mt. Wilhelm, the highest mountain of Papua New Guinea. Despite species-specific differences in elevational preference, we found limited differentiation within each scrubwren species, but detected a strong genomic signature of simultaneous population expansions at 27-29 ka, coinciding with the onset of the LGM.

CONCLUSION

The remarkable synchronous timing of population expansions of all three species demonstrates the importance of global cooling cycles in expanding highland habitat. Global cooling cycles have likely had strongly different impacts on tropical montane areas versus boreal and temperate latitudes, leading to population expansions in the former and serious fragmentation in the latter.

Hybridizing salamanders experience accelerated diversification

Whether hybridization generates or erodes species diversity has long been debated, but to date most studies have been conducted at small taxonomic scales. Salamanders (order Caudata) represent a taxonomic order in which hybridization plays a prevalent ecological and evolutionary role. We employed a recently developed model of trait-dependent diversification to test the hypothesis that hybridization impacts the diversification dynamics of species that are currently hybridizing. We find strong evidence supporting this hypothesis, showing that hybridizing salamander lineages have significantly greater net-diversification rates than non-hybridizing lineages. This pattern is driven by concurrently increased speciation rates and decreased extinction rates in hybridizing lineages. Our results support the hypothesis that hybridization can act as a generative force in macroevolutionary diversification.

Phylogeography of the cicada Platypleura hilpa in subtropical and tropical East Asia based on mitochondrial and nuclear genes and microsatellite markers

Climate change and geographical events play key roles in driving population genetic structure of organisms, but different scenarios were suggested for species occurring in tropical and subtropical areas. We investigated the population genetic structure, diversity and demographic history of the cicada Platypleura hilpa which occurs in coastal areas of southern China and northeastern Indo-Burma, and analysed the potential impact of climate change and geological events on its evolutionary history. Our data imply that P. hilpa comprises five main lineages with nearly unique sets of haplotypes and distinct geographic distributions. A major split of the lineages occurred in the Pleistocene. Geographic distance and geomorphic barriers serve as the primary factors shaping the genetic population structure, and several climatic factors are associated with the divergence. The potential range during the Last Glacial Maximum has apparently increased in south China and the exposed South China Sea Shelf. The Pleistocene sea-level fluctuations had profound effects on the regional genetic structure. The Beibu Gulf represents a more important geographic barrier than the Qiongzhou Strait in blocking gene flow among populations. These results contribute to a better understanding of the pressure climatic change and geographical events impose on insects in coastal areas of East Asia.

Multilocus phylogeny and cryptic diversity of white-toothed shrews (Mammalia, Eulipotyphla, Crocidura) in China

BACKGROUND

Crocidura, the most speciose mammalian genus, occurs across much of Asia, Europe and Africa. The taxonomy of Chinese representatives has been studied primarily based on cursory morphological comparisons and their molecular phylogenetic analyses remain unexplored. In order to understand the phylogeny of this group in China, we estimated the first multilocus phylogeny and conducted species delimitation, including taxon sampling throughout their distribution range.

RESULTS

We obtained one mitochondrial gene (cytb) (~ 1, 134 bp) and three nuclear genes (ApoB, BRCA1, RAG1) (~ 2, 170 bp) for 132 samples from 57 localities. Molecular analyses identified at least 14 putative species that occur within two major well-supported groups in China. Polyphyletic C. wuchihensis appears to be composed of two putative species. Two subspecies, C. rapax rapax and C. rapax kurodai should be elevated to full species status. A phylogenetic tree based on mitochondrial gene from Asian Crocidura species showed that the C. rapax rapax is embedded within C. attenuata, making the latter a paraphyletic group. Three strongly supported undescribed species (C. sp.1, C. sp.2 and C. sp.3) are revealed from Zada County of Tibet (Western China), Hongjiang County of Hunan Province (Central China) and Dongyang County of Zhejiang Province (Eastern China), Motuo County of Tibet, respectively. The divergence time estimation suggested that China's Crocidura species began to diversify during the late Pliocene (3.66 Ma) and the Early Pleistocene (2.29 Ma), followed by a series of diversifications through the Pleistocene.

CONCLUSIONS

The cryptic diversity found in this study indicated that the number of species is strongly underestimated under the current taxonomy. We propose that the three undescribed species should be evaluated using extensive taxon sampling and comprehensive morphological and morphometric approaches. Climate change since the late Pliocene and the uplift of the Qinghai-Tibet Plateau may result in the diversification and speciation of China's Crocidura species. In short, the underestimated diversity underlines the need for a taxonomic revision of Chinese Crocidura species.

Hidden species diversity in Pachyhynobius: A multiple approaches species delimitation with mitogenomes

The lack of distinct morphological features of cryptic species is a hard problem for taxonomy, especially when the taxa are closely related with considerable amounts of ancestral polymorphism. Lately, intensive coalescent-based analyses involving multiple loci have become the preferred method to assess the extent of genetic distinctiveness in otherwise phenotypically similar populations. Previously, phylogenetic studies on Pachyhynobius shangchengensis uncovered five extremely deeply divergent clades, which suggested that this species may be a cryptic species complex. In this study, we used the complete mitochondrial genome data and samples from the entire range of stout salamander (Pachyhynobius), as well as publicly available mitochondrial genomes to assess species boundaries within this genus using a suite of diverse methodologies (e.g. general mixed Yule coalescent model, Automatic Barcode Gap Discovery). The phylogenetic relationships recovered two major groups within P. shangchengensis, with one group formed by four of the six extant populations and corresponding to the central and eastern range of the Dabie mountains, while the other group encompassed two other lineages in the north west of the Dabie mountain range. The species delimitation comparison within Pachyhynobius supported the presence of recognized species within the genus, and consensus was observed across methods for the existence of up to five cryptic species within what has been traditionally considered to be P. shangchengensis. While this implies the existence of four taxa in addition to the described P. shangchengensis species, morphological data and life history information are further required to contribute to the species definition. The observed pattern of genetic variation is likely the outcome of a discontinuous habitat combined with niche conservatism, which produced the sky-island effect observed in Pachyhynobius, and which led to formation of a hidden species diversity in this genus.

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.

Past vicariance promoting deep genetic divergence in an endemic frog species of the Espinhaço Range in Brazil: The historical biogeography of Bokermannohyla saxicola (Hylidae)

The highland endemic frog Bokermannohyla saxicola occurs within the Espinhaço Range, the most extensive and continuous orogenic belt of the Brazilian territory, located in southeastern Brazil. We used mitochondrial DNA markers to test for spatial structure, to investigate the likely influence of past vicariant events, to evaluate demographic dynamics along the species range, and to understand the role of habitat discontinuities in promoting connectivity and diversity along the range. We found four major monophyletic lineages, each one associated with distinct mountain tops. The divergence time found between the four main clades clearly pre-dated the Pleistocene, except for the most recent separation. We observed no signs of population expansion for most of the sampling sites along the range, and a higher genetic diversity in the most continuous and central highland plateau, compared to smaller marginal regions. The Espinhaço Range harbors four deeply divergent lineages of B. saxicola within areas restricted by barriers for millions of years. These relatively isolated populations were kept apart by discontinuities represented by lowland habitats between mountain tops. Most of the lineage divergences occurred earlier than the Pleistocene, thus they cannot be solely explained by climatic oscillations of this epoch. However, within-lineage divergence times were all dated from the Pleistocene, suggesting an important effect in population dynamics. We also suggest that some marginal populations like those from Serra Negra and Serra de Itacambira can be the result of recent colonization events. Finally, in the southern Espinhaço region, the most continuous central highland area shows greater genetic diversity than the marginal discontinuous areas, where we have also observed a significant correlation between genetic and geographic distances. Bokermannohyla saxicola is a good model to study the biogeography of the Espinhaço Range because its high genetic structure reflects ancient as well as recent geological/climatic events, with important implications for conservation.

Phylogenetics of Kingsnakes, Lampropeltis getula Complex (Serpentes: Colubridae), in Eastern North America

Kingsnakes of the Lampropeltis getula complex range throughout much of temperate and subtropical North America. Studies over the last century have used morphology and color pattern to describe numerous subspecies. More recently, DNA analyses have made invaluable contributions to our understanding of their evolution and taxonomy. We use genetic and ecological methods to test previous hypotheses of distinct evolutionary lineages by examining 66 total snakes and 1) analyzing phylogeographic structure using 2 mtDNA loci and 1 nuclear locus, 2) estimating divergence dates and historical demography among lineages in a Bayesian coalescent framework, and 3) applying ecological niche modeling (ENM). Our molecular data and ENMs illustrate that 3 previously recognized subspecies in the eastern United States comprise well-supported monophyletic lineages that diverged during the Pleistocene. The geographic boundaries of these 3 lineages correspond closely to known biogeographic barriers (Florida peninsula, Appalachian Mountains, and Apalachicola River) previously identified for other plants and animals, indicating shared geographic influences on evolutionary history. We conclude that genetic, ecological, and morphological data support recognition of these 3 lineages as distinct species (Lampropeltis floridana, Lampropeltis getula, and Lampropeltis meansi).

Pliocene-Pleistocene lineage diversifications in the Eastern Indigo Snake (Drymarchon couperi) in the Southeastern United States

Indigo Snakes (Drymarchon; with five currently recognized species) occur from northern Argentina, northward to the United States in southern Texas and eastward in disjunct populations in Florida and Georgia. Based on this known allopatry and a difference in supralabial morphology the two United States taxa previously considered as subspecies within D. corais (Boie 1827), the Western Indigo Snake, D. melanurus erebennus (Cope 1860), and Eastern Indigo Snake, D. couperi (Holbrook 1842), are currently recognized as separate species. Drymarchon couperi is a Federally-designated Threatened species by the United States Fish and Wildlife Service under the Endangered Species Act, and currently being incorporated into a translocation program. This, combined with its disjunct distribution makes it a prime candidate for studying speciation and genetic divergence. In this study, we (1) test the hypothesis that D. m. erebennus and D. couperi are distinct lineages by analyzing 2411 base pairs (bp) of two mitochondrial (mtDNA) loci and one single copy nuclear (scnDNA) locus; (2) estimate the timing of speciation using a relaxed phylogenetics method to determine if Milankovitch cycles during the Pleistocene might have had an influence on lineage diversifications; (3) examine historical population demography to determine if identified lineages have undergone population declines, expansions, or remained stable during the most recent Milankovitch cycles; and (4) use this information to assist in an effective and scientifically sound translocation program. Our molecular data support the initial hypothesis that D. melanurus and D. couperi should be recognized as distinct species, but further illustrate that D. couperi is split into two distinct genetic lineages that correspond to historical biogeography and sea level changes in peninsular Florida. These two well-supported genetic lineages (herein termed Atlantic and Gulf lineages) illustrate a common biogeographic distributional break previously identified for other plants and animals, suggesting that these organisms might have shared a common evolutionary history related to historic sea level changes caused by Milankovitch cycles. Our estimated divergence times suggest that the most recent common ancestor (MRCA) between D. melanurus and southeastern United States Drymarchon occurred ca. 5.9Ma (95% HPD=2.5-9.8Ma; during the late Blancan of the Pleistocene through the Hemphillian of the Miocene), whereas the MRCA between the Atlantic and Gulf lineages in the southeastern United States occurred ca. 2.0Ma (95% HPD=0.7-3.7Ma; during the Irvingtonian of the Pleistocene through the Blancan of the Pliocene). During one or more glacial intervals within these times, these two lineages must have become separated and evolved independently. Despite numerous Milankovitch cycles along with associated forming of physical barriers (i.e., sea level fluctuations, high elevation sand ridges, clayey soils, and/or insufficient habitats) since their initial lineage diversification, these two lineages have likely come in and out of contact with each other many times, yet today they still illustrate near discrete geographic distributions. Although the Atlantic and Gulf lineages appear to be cryptic, a thorough study examining morphological characters should be conducted. We believe that our molecular data is crucial and should be incorporated in making conscious decisions in the management of a translocation program. We suggest that source populations for translocations include maintaining the integrity of the known genetic lineages found herein, as well as those coming from the closest areas that currently support sizable Drymarchon populations.

Thriving in the Cold: Glacial Expansion and Post-Glacial Contraction of a Temperate Terrestrial Salamander (Plethodon serratus)

The dynamic geologic history of the southeastern United States has played a major role in shaping the geographic distributions of amphibians in the region. In the phylogeographic literature, the predominant pattern of distribution shifts through time of temperate species is one of contraction during glacial maxima and persistence in refugia. However, the diverse biology and ecology of amphibian species suggest that a "one-size-fits-all" model may be inappropriate. Nearly 10% of amphibian species in the region have a current distribution comprised of multiple disjunct, restricted areas that resemble the shape of Pleistocene refugia identified for other temperate taxa in the literature. Here, we apply genetics and spatially explicit climate analyses to test the hypothesis that the disjunct regions of these species ranges are climatic refugia for species that were more broadly distributed during glacial maxima. We use the salamander Plethodon serratus as a model, as its range consists of four disjunct regions in the Southeast. Phylogenetic results show that P. serratus is comprised of multiple genetic lineages, and the four regions are not reciprocally monophyletic. The Appalachian salamanders form a clade sister to all other P. serratus. Niche and paleodistribution modeling results suggest that P. serratus expanded from the Appalachians during the cooler Last Glacial Maximum and has since been restricted to its current disjunct distribution by a warming climate. These data reject the universal applicability of the glacial contraction model to temperate taxa and reiterate the importance of considering the natural history of individual species.

Molecular phylogenetics and phylogeographic structure of Sorex bedfordiae based on mitochondrial and nuclear DNA sequences

The southeastern margin of the Tibetan Plateau is characterized by complex topography and a discontinuous landscape, creating a sky island situation. However, the way topography shapes genetic structures and demographic histories of endemic species has not been well studied. We examined the phylogeographic pattern and demographic histories of Sorex bedfordiae, a dispersal-limited small mammal, using three nuclear genes [1977bp] and two mitochondrial genes [1794bp] with comprehensive molecular approaches. We recovered five well-supported clades whose distributions are along mountain ridges and roughly subdivided by large rivers. Demographic expansions in the middle Pleistocene were strongly supported by both nuclear and mitochondrial genes. Our results support the hypothesis that sky island topography and river systems strongly affect the genetic structure of non-aquatic terrestrial species. We further clarify that S. bedfordiae and S. cylindricauda are valid sibling species, whereas S. excelsus is most likely a geographic subspecies of S. bedfordiae.

Late Pleistocene fossils and the future distribution of Rana temporaria (Amphibia, Anura) along the Apennine Peninsula (Italy)

Background

The effect of past climatic changes on the distribution of organisms is a fertile field of research that has been tackled in many different ways. Because the fossil record provides direct access to the chronological and geographic dimensions of biological events occurred in the past, it can be a useful tool for assessing range contractions and expansions related to climatic changes.

Results

Here, we provide support for the ‘recent’ shrinkage of the range of a frigophilous anuran, the common frog, Rana temporaria Linnaeus, 1758, by analyzing the amphibian fossil assemblage coming from the Grotta di Equi, a Late Pleistocene site (about 45 ka) located in the Apuan Alps (northern Apennine chain, Massa-Carrara Province, Italy). Besides, R. temporaria, the assemblage includes cf. Salamandra salamandra Linnaeus, 1758 and Bufo bufo Linnaeus, 1758.

Conclusions

The presence of R. temporaria in the Apuan Alps during a cold interglacial phase at an elevation that is much lower than the minimum current elevation in the region (about 300 m lower) supports the previously reported hypothesis that the Pleistocene coolings favored the dispersal of this frigophilous species along the Apennine chain and that its current, locally disjointed distribution is the effect of the post Pleistocene raise in temperature. Our results concur in predicting that the current global warming will further affect, negatively, the range of this frigophilous species by further reducing its distribution in the southern sectors of its range.

Niche divergence promotes rapid diversification of East African sky island white-eyes (Aves: Zosteropidae)

The Eastern Afromontane biodiversity hotspot composed of highly fragmented forested highlands (sky islands) harbours exceptional diversity and endemicity, particularly within birds. To explain their elevated diversity within this region, models founded on niche conservatism have been offered, although detailed phylogeographic studies are limited to a few avian lineages. Here, we focus on the recent songbird genus Zosterops, represented by montane and lowland members, to test the roles of niche conservatism versus niche divergence in the diversification and colonization of East Africa's sky islands. The species-rich white-eyes are a typically homogeneous family with an exceptional colonizing ability, but in contrast to their diversity on oceanic islands, continental diversity is considered depauperate and has been largely neglected. Molecular phylogenetic analysis of ~140 taxa reveals extensive polyphyly among different montane populations of Z. poliogastrus. These larger endemic birds are shown to be more closely related to taxa with divergent habitat types, altitudinal distributions and dispersal abilities than they are to populations of restricted endemics that occur in neighbouring montane forest fragments. This repeated transition between lowland and highland habitats over time demonstrate that diversification of the focal group is explained by niche divergence. Our results also highlight an underestimation of diversity compared to morphological studies that has implications for their taxonomy and conservation. Molecular dating suggests that the spatially extensive African radiation arose exceptionally rapidly (1-2.5 Ma) during the fluctuating Plio-Pleistocene climate, which may have provided the primary driver for lineage diversification.

Population structure, historical biogeography and demographic history of the alpine toad Scutiger ningshanensis in the Tsinling Mountains of Central China

Population genetic structure, historical biogeography and historical demography of the alpine toad Scutiger ningshanensis were studied using the combined data mtDNA cytochrome b (cyt b) and the mtDNA cytochrome c oxidase subunit I (COI) as the molecular markers. This species has high genetic variation. There was a significant genetic differentiation among most populations. Three lineages were detected. The phylogenetic relationship analyses and the SAMOVA (spatial analysis of molecular variance) results showed significant phylogeographic structure. 82.15% genetic variation occurred among populations whereas differentiation within populations only contributed 17.85% to the total. Mantel test results showed a significant correlation between the pairwise calculated genetic distance and pairwise calculated geographical distance of the populations (regression coefficient  = 0.001286, correlation coefficient  = 0.77051, p (rrand≥robs)  = 0.0185<0.05), indicating the existence of isolation-by-distance pattern of genetic divergence for cyt b + COI sequence, which suggests that the distribution of genetic variation is due to geographical separation rather than natural selection. The population expansion or contraction and genetic differentiation between populations or lineages could be explained by topography and the repetitive uplifts of the Tsinling Mountains and the climatic cycles during the late Pliocene and Pleistocene. S. ningshanensis experienced a rapid population expansion about 40,000 years before present. The current decline in population size was probably caused by anthropogenic disturbance. Current populations of S. ningshanensis are from different refugia though the location of these refugia could not be determined in our study. Topography, climatic changes and repetitive population expansion/contraction together led to the high level of genetic variation in S. ningshanensis. A total of three management units (MUs) was determined, which must be considered when conservation policy is made in the future.

Integrating multiple lines of evidence into historical biogeography hypothesis testing: a Bison bison case study

One of the grand goals of historical biogeography is to understand how and why species' population sizes and distributions change over time. Multiple types of data drawn from disparate fields, combined into a single modelling framework, are necessary to document changes in a species's demography and distribution, and to determine the drivers responsible for change. Yet truly integrated approaches are challenging and rarely performed. Here, we discuss a modelling framework that integrates spatio-temporal fossil data, ancient DNA, palaeoclimatological reconstructions, bioclimatic envelope modelling and coalescence models in order to statistically test alternative hypotheses of demographic and potential distributional changes for the iconic American bison (Bison bison). Using different assumptions about the evolution of the bioclimatic niche, we generate hypothetical distributional and demographic histories of the species. We then test these demographic models by comparing the genetic signature predicted by serial coalescence against sequence data derived from subfossils and modern populations. Our results supported demographic models that include both climate and human-associated drivers of population declines. This synthetic approach, integrating palaeoclimatology, bioclimatic envelopes, serial coalescence, spatio-temporal fossil data and heterochronous DNA sequences, improves understanding of species' historical biogeography by allowing consideration of both abiotic and biotic interactions at the population level.

Paleoclimatic modeling and phylogeography of least killifish, Heterandria formosa: insights into Pleistocene expansion-contraction dynamics and evolutionary history of North American Coastal Plain freshwater biota

BACKGROUND

Climatic and sea-level fluctuations throughout the last Pleistocene glacial cycle (~130-0 ka) profoundly influenced present-day distributions and genetic diversity of Northern Hemisphere biotas by forcing range contractions in many species during the glacial advance and allowing expansion following glacial retreat ('expansion-contraction' model). Evidence for such range dynamics and refugia in the unglaciated Gulf-Atlantic Coastal Plain stems largely from terrestrial species, and aquatic species Pleistocene responses remain relatively uninvestigated. Heterandria formosa, a wide-ranging regional endemic, presents an ideal system to test the expansion-contraction model within this biota. By integrating ecological niche modeling and phylogeography, we infer the Pleistocene history of this livebearing fish (Poeciliidae) and test for several predicted distributional and genetic effects of the last glaciation.

RESULTS

Paleoclimatic models predicted range contraction to a single southwest Florida peninsula refugium during the Last Glacial Maximum, followed by northward expansion. We inferred spatial-population subdivision into four groups that reflect genetic barriers outside this refuge. Several other features of the genetic data were consistent with predictions derived from an expansion-contraction model: limited intraspecific divergence (e.g. mean mtDNA p-distance = 0.66%); a pattern of mtDNA diversity (mean Hd = 0.934; mean π = 0.007) consistent with rapid, recent population expansion; a lack of mtDNA isolation-by-distance; and clinal variation in allozyme diversity with higher diversity at lower latitudes near the predicted refugium. Statistical tests of mismatch distributions and coalescent simulations of the gene tree lent greater support to a scenario of post-glacial expansion and diversification from a single refugium than to any other model examined (e.g. multiple-refugia scenarios).

CONCLUSIONS

Congruent results from diverse data indicate H. formosa fits the classic Pleistocene expansion-contraction model, even as the genetic data suggest additional ecological influences on population structure. While evidence for Plio-Pleistocene Gulf Coast vicariance is well described for many freshwater species presently codistributed with H. formosa, this species demography and diversification departs notably from this pattern. Species-specific expansion-contraction dynamics may therefore have figured more prominently in shaping Coastal Plain evolutionary history than previously thought. Our findings bolster growing appreciation for the complexity of phylogeographical structuring within North America's southern refugia, including responses of Coastal Plain freshwater biota to Pleistocene climatic fluctuations.

Integrative phylogeography of Calotriton newts (Amphibia, Salamandridae), with special remarks on the conservation of the endangered Montseny brook newt (Calotriton arnoldi)

The genus Calotriton includes two species of newts highly adapted to live in cold and fast-flowing mountain springs. The Pyrenean brook newt (Calotriton asper), restricted to the Pyrenean region, and the Montseny brook newt (Calotriton arnoldi), endemic to the Montseny massif and one of the most endangered amphibian species in Europe. In the present manuscript, we use an integrative approach including species distribution modeling (SDM), molecular analyses of mitochondrial and nuclear DNA sequence data and morphology to unravel the historical processes that have contributed to shaping the biogeography and genetic structure of the genus Calotriton, with special emphasis on the conservation of C. arnoldi. The results of the molecular analyses confirm that, despite having originated recently, being ecologically similar and geographically very close, there is no signal of hybridization between C. asper and C. arnoldi. SDM results suggest that tough environmental conditions on mountains tops during glacial periods, together with subsequent warmer periods could have prevented the contact between the two species. Within the critically endangered C. arnoldi, a high genetic structure is revealed despite its extremely small distribution range compared to C. asper. Haplotype networks, AMOVA and SAMOVA analyses suggest that two distinct groups of populations can be clearly differentiated with absence of gene flow. This is in concordance with morphological differentiation and correlates with its geographical distribution, as the two groups are situated on the eastern and western sides of a river valley that acts as a barrier. The genetic and morphological results are highly important for the ongoing conservation program of C. arnoldi and strongly justify the management of this species into at least two independent evolutionary significant units (eastern and western sectors) to guarantee the long-term population viability.

Population expanding with the phalanx model and lineages split by environmental heterogeneity: a case study of Primula obconica in subtropical China

BACKGROUND

Current and historical events have both affected the current distribution patterns and intraspecific divergence of plants. While numerous studies have focused on the Qinghai-Tibetan Plateau (QTP), the impacts of such events on the flora of subtropical China remain poorly understood. Subtropical China is famous for its highly complex topography and the limited impact from glaciation during the Pleistocene; this may have resulted in a different genetic legacy for species in this region compared to fully glaciated areas.

METHODOLOGY/PRINCIPAL FINDINGS

We used plastid and nuclear DNA sequence data and distribution modeling to analyze the divergence patterns and demographic history of Primula obconica Hance, a widespread herbaceous montane species in subtropical China. The phylogenetic analysis revealed two major lineages (lineage A and lineage B), representing a west-east split into the Yunnan and Eastern groups, and the Sichuan and Central groups, respectively. The Eastern and Central groups comprised relatively new derived haplotypes. Nested Clade Analysis and Bayesian Skyline Plot analyses both indicated that P. obconica mainly experienced a gradual expansion of populations. In addition, the simulated distribution of P. obconica during the Last Glacial Maximum was slightly larger than its present-day distribution.

CONCLUSION/SIGNIFICANCE

Our results are the first to identify a west-east migration of P. obconica. The gradual expansion pattern and a larger potential distribution range in cold periods detected for P. obconica indicate that the population expansion of this species is consistent with the phalanx model. In addition, the current patterns of genetic differentiation have persisted as a result of the extensive environmental heterogeneity that exists in subtropical China.

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.

Phylogeography of Sardinian cave salamanders (genus Hydromantes) is mainly determined by geomorphology

Detecting the factors that determine the interruption of gene flow between populations is key to understanding how speciation occurs. In this context, caves are an excellent system for studying processes of colonization, differentiation and speciation, since they represent discrete geographical units often with known geological histories. Here, we asked whether discontinuous calcareous areas and cave systems represent major barriers to gene flow within and among the five species of Sardinian cave salamanders (genus Hydromantes) and whether intraspecific genetic structure parallels geographic distance within and among caves. We generated mitochondrial cytochrome b gene sequences from 184 individuals representing 48 populations, and used a Bayesian phylogeographic approach to infer possible areas of cladogenesis for these species and reconstruct historical and current dispersal routes among distinct populations. Our results show deep genetic divergence within and among all Sardinian cave salamander species, which can mostly be attributed to the effects of mountains and discontinuities in major calcareous areas and cave systems acting as barriers to gene flow. While these salamander species can also occur outside caves, our results indicate that there is a very poor dispersal of these species between separate cave systems.