Late‐glacial recolonization and phylogeography of European red deer (Cervus elaphus L.)

The Pleistocene was an epoch of extreme climatic and environmental changes. How individual species responded to the repeated cycles of warm and cold stages is a major topic of debate. For the European fauna and flora, an expansion–contraction model has been suggested, whereby temperate species were restricted to southern refugia during glacial times and expanded northwards during interglacials, including the present interglacial (Holocene). Here, we test this model on the red deer (Cervus elaphus) a large and highly mobile herbivore, using both modern and ancient mitochondrial DNA from the entire European range of the species over the last c. 40 000 years. Our results indicate that this species was sensitive to the effects of climate change. Prior to the Last Glacial Maximum (LGM) haplogroups restricted today to South‐East Europe and Western Asia reached as far west as the UK. During the LGM, red deer was mainly restricted to southern refugia, in Iberia, the Balkans and possibly in Italy and South‐Western Asia. At the end of the LGM, red deer expanded from the Iberian refugium, to Central and Northern Europe, including the UK, Belgium, Scandinavia, Germany, Poland and Belarus. Ancient DNA data cannot rule out refugial survival of red deer in North‐West Europe through the LGM. Had such deer survived, though, they were replaced by deer migrating from Iberia at the end of the glacial. The Balkans served as a separate LGM refugium and were probably connected to Western Asia with genetic exchange between the two areas.     

Contrasting genetic patterns between two coexisting Eleutherococcus species in northern China

Climate oscillations are the key factors to understand the patterns in modern biodiversity. East Asia harbors the most diverse temperate flora, largely because an extensive terrestrial ice cap was absent during repeated Pleistocene glaciation–interglacial cycles. Comparing the demographic histories of species that are codistributed and are close relatives may provide insight into how the process of climate change influences species ranges. In this study, we compared the spatial genetic structure and demographic histories of two coexisting Eleutherococcus species, Eleutherococcus senticosus and E. sessiliflorus. Both species are distributed in northern China, regions that are generally considered to be sensitive to climatic fluctuations. These regions once hosted temperate forest, but this temperate forest was replaced by tundra and taiga forest during the Last Glacial Maximum (LGM), according to pollen records. Using three chloroplast DNA fragments, we assessed the genetic structure of 20 and 9 natural populations of E. senticosus and E. sessiliflorus, respectively. Extremely contrasting genetic patterns were found between the two species; E. sessiliflorus had little genetic variation, whereas E. senticosus had considerably higher levels of genetic variation (15 haplotypes). We speculated that a recent severe bottleneck may have resulted in the extremely low genetic diversity in E. sessiliflorus. In E. senticosus, populations in Northeast China (NEC) harbored all of the haplotypes found in this species and included private haplotypes. The populations in NEC had higher levels of genetic diversity than did those from North China (NC). Therefore, we suggest that both the NC and NEC regions can sustain LGM refugia and that lineage admixture from multiple refugia took place after the LGM elevated the local genetic diversity in NEC. In NEC, multiple genetic hot spots were found in the Changbai Mountains and the Xiaoxing'an Range, which implied that multiple locations in NEC may sustain LGM refugia, even in the Xiaoxing'an Range.     

Genetic and paleomodelling evidence of the population expansion of the cattle egret Bubulcus ibis in Africa during the climatic oscillations of the Late Pleistocene

Increasing aridity during glacial periods produced the retraction of forests and the expansion of arid and semi‐arid environments in Africa, with consequences for birds. Cattle egret Bubulcus ibis is a dispersive species that prefers semi‐arid environments and requires proximity to bodies of water. We expected that climatic oscillations led to the expansion of the range of the cattle egret during arid periods, such as the Last Maximum Glacial (LGM) and contraction of distribution during the Last Interglacial (LIG) period, resulting in contact of populations previously isolated. We investigated this hypothesis by evaluating the genetic structure and population history of 15 cattle egret breeding colonies located in west and South Africa using the mitochondrial DNA (mtDNA) control region, mtDNA ATPase 8 and 6, and an intron of nuclear gene transforming growth factor‐beta 2. Occurrence data and bioclimatic information were used to generate ecological niche models of three periods (present, LGM and LIG). We used the genetic and paleomodelling data to assess the responses of the cattle egret from Africa to the climatic oscillations during the late Pleistocene. Genetic data revealed low levels of genetic differentiation, signs of isolation‐by‐distance, as well as recent increases in effective population size that started during the LGM. The observed low genetic structure may be explained by recent colonization events due to the demographic expansion following the last glacial period and by dispersal capacity of this species. The paleomodels corroborated the expansion during the LGM, and a more restricted potential distribution during the LIG. Our findinds supports the hypothesis that the species range of the cattle egret expanded during arid periods and contracted during wet periods.     

Uncommon paleodistribution patterns of Chrysolophus pheasants in east Asia: explanations and implications

Some modeling studies indicated that the past distributions of species in east Asia during the Last Interglacial (LIG) and Last Glacial Maximum (LGM) periods differ from those of European and North American species and the deviant Asian distribution pattern is known under the term ‘pre‐LGM expansion’. It represents the unusually similar distribution patterns between the current and the LGM scenario. However, there is still no satisfying explanation for this phenomenon so far. Therefore, we took the two recently separated pheasant species of genus Chrysolophus in east Asia as an example to test the pattern by performing ecological niche models. The main findings of this study include: 1) the paleodistributions of these two pheasants also corresponded to the ‘pre‐LGM expansion’ pattern; 2) climatic similarity results from mobility‐oriented parity analysis also revealed similar pattern for both species; 3) climate regimes of east Asia showed patterns different from those in Europe and North America in a climate shift towards drier conditions and stronger seasonality and to more extreme temperatures of the coldest months particularly during the LIG; 4) the two Chrysolophus species occupied significantly different ecological niches according to current distribution. We suggest that ecological segregation established in allopatric glacial refugia should be the main determinants for the separation of two Chrysolophus species until they came into extant post‐Pleistocene contact.     

East Asian allopatry and north Eurasian sympatry in Long‐tailed Tit lineages despite similar population dynamics during the late Pleistocene

Eurasia is a large continent characterized by heterogeneous environments. Glacial cycles during the late Pleistocene have had variable impacts on the avifauna across Eurasia. Bird populations from South‐East Asia show stability through the Last Glacial Maximum (LGM), while populations from Europe exhibit evidence of post‐LGM expansion. We investigated the phylogeography of the Long‐tailed Tit (Aegithalos caudatus), which spans the longitudinal breadth of Eurasia to test how climatic history and regional topographical complexity affected populations and diversification within the species complex. Our results show that two lineages from central and southern China (lineages C and D) segregate geographically, while lineages across northern Eurasia (lineage A and B) show substantial sympatry. Bayesian estimates for the timing of diversification suggest that the four lineages diverged during the middle Pleistocene, splitting in parallel and undergoing concurrent demographic histories since divergence. A. caudatus lineages experienced similar and synchronous population size dynamics during glacial cycles before the LGM. We conclude that the difference in geo‐topologic complexity may be an important factor that led to the variation in secondary admixture between northern Eurasian and eastern Asian lineages.     

Profound climatic effects on two East Asian black-throated tits (Ave: Aegithalidae), revealed by ecological niche models and phylogeographic analysis

Although a number of studies have assessed the effects of geological and climatic changes on species distributions in East Asian, we still have limited knowledge of how these changes have impacted avian species in south-western and southern China. Here, we aim to study paleo-climatic effects on an East Asian bird, two subspecies of black-throated tit (A. c. talifuensis-concinnus) with the combined analysis of phylogeography and Ecological Niche Models (ENMs). We sequenced three mitochondrial DNA markers from 32 populations (203 individuals) and used phylogenetic inferences to reconstruct the intra-specific relationships among haplotypes. Population genetic analyses were undertaken to gain insight into the demographic history of these populations. We used ENMs to predict the distribution of target species during three periods; last inter-glacial (LIG), last glacial maximum (LGM) and present. We found three highly supported, monophyletic MtDNA lineages and different historical demography among lineages in A. c. talifuensis-concinnus. These lineages formed a narrowly circumscribed intra-specific contact zone. The estimated times of lineage divergences were about 2.4 Ma and 0.32 Ma respectively. ENMs predictions were similar between present and LGM but substantially reduced during LIG. ENMs reconstructions and molecular dating suggest that Pleistocene climate changes had triggered and shaped the genetic structure of black-throated tit. Interestingly, in contrast to profound impacts of other glacial cycles, ENMs and phylogeographic analysis suggest that LGM had limited effect on these two subspecies. ENMs also suggest that Pleistocene climatic oscillations enabled the formation of the contact zone and thus support the refuge theory.     

Multiple Pleistocene refugia and Holocene range expansion of an abundant southwestern American desert plant species (Melampodium leucanthum,Asteraceae)

Pleistocene climatic fluctuations had major impacts on desert biota in southwestern North America. During cooler and wetter periods, drought‐adapted species were isolated into refugia, in contrast to expansion of their ranges during the massive aridification in the Holocene. Here, we use Melampodium leucanthum (Asteraceae), a species of the North American desert and semi‐desert regions, to investigate the impact of major aridification in southwestern North America on phylogeography and evolution in a widespread and abundant drought‐adapted plant species. The evidence for three separate Pleistocene refugia at different time levels suggests that this species responded to the Quaternary climatic oscillations in a cyclic manner. In the Holocene, once differentiated lineages came into secondary contact and intermixed, but these range expansions did not follow the eastwardly progressing aridification, but instead occurred independently out of separate Pleistocene refugia. As found in other desert biota, the Continental Divide has acted as a major migration barrier for M. leucanthum since the Pleistocene. Despite being geographically restricted to the eastern part of the species’ distribution, autotetraploids in M. leucanthum originated multiple times and do not form a genetically cohesive group.     

Tracing the legacy of the early Hainan Islanders - a perspective from mitochondrial DNA

Background  

Hainan Island is located around the conjunction of East Asia and Southeast Asia, and during the Last Glacial Maximum (LGM) was connected with the mainland. This provided an opportunity for the colonization of Hainan Island by modern human in the Upper Pleistocene. Whether the ancient dispersal left any footprints in the contemporary gene pool of Hainan islanders is debatable.     

Spatial genetic structure and historical demography of East Asian wild boar

Pleistocene climatic fluctuations may have had a profound impact on the evolutionary history of many species. The geographical pattern of European wild boar (Sus scrofa) is clearly studied, and it was greatly influenced by ancient climatic events, especially the Last Glacial Maximum. Previous research on genetic variation has mainly focused on the origin and distribution histories of domestic pigs. However, some questions have not been answered, including those concerning the genetic diversity, geographical pattern and possible historic influence of climate on East Asian wild boar (EAWB). Employing the control region of mtDNA (511 bp), we investigated the contributions of historic climate, which possibly shaped the genetic pattern of wild boar. Given that the level of genetic diversity of wild boars is higher in East Asia than in Europe, 172 haplotypes were detected from 680 individuals. Phylogenetic analysis demonstrated the complex phylogeographic structure of EAWB. Mismatch analysis, neutrality tests and the Bayesian Skyline Plot results all retrieved signals of a rapid population expansion, which might have played an important role in driving the formation of complex spatial genetic structure. Genetic data and species distribution modelling showed that the Last Glacial Maximum had weak effect on the distribution of the EAWB. We suggest that, in shaping spatial genetic structure in East Asian, long-term gene flow and population history played more important roles than Pleistocene climate fluctuations.     

Refugia,colonization and diversification of an arid‐adapted bird: coincident patterns between genetic data and ecological niche modelling

Phylogeographical studies are common in boreal and temperate species from the Palaearctic, but scarce in arid‐adapted species. We used nuclear and mitochondrial markers to investigate phylogeography and to estimate chronology of colonization events of the trumpeter finch Bucanetes githagineus, an arid‐adapted bird. We used 271 samples from 16 populations, most of which were fresh samples but including some museum specimens. Microsatellite data showed no clear grouping according to the sampling locations. Microsatellite and mitochondrial data showed the clearest differentiation between Maghreb and Canary Islands and between Maghreb and Western Sahara. Mitochondrial data suggest differentiation between different Maghreb populations and among Maghreb and Near East populations, between Iberian Peninsula and Canary Islands, as well as between Western Sahara and Maghreb. Our coalescence analyses indicate that the trumpeter finch colonized North Africa during the humid Marine Isotope Stage 5 (MIS5) period of the Sahara region 125 000 years ago. We constructed an ecological niche model (ENM) to estimate the geographical distribution of climatically suitable habitats for the trumpeter finch. We tested whether changes in the species range in relation to glacial–interglacial cycles could be responsible for observed patterns of genetic diversity and structure. Modelling results matched with those from genetic data as the species' potential range increases in interglacial scenarios (in the present climatic scenario and during MIS5) and decreases in glacial climates (during the last glacial maximum, LGM, 21 000 years ago). Our results suggest that the trumpeter finch responded to Pleistocene climatic changes by expanding and contracting its range.     

Repeated landmass reformation limits diversification in the widespread littoral zone mosquito Anopheles sundaicus sensu lato in the Indo‐Oriental Region

Southeast Asia harbours abundant biodiversity, hypothesized to have been generated by Pliocene and Pleistocene climatic and environmental change. Vicariance between the island of Borneo, the remaining Indonesian archipelago and mainland Southeast Asia caused by elevated sea levels during interglacial periods has been proposed to lead to diversification in the littoral zone mosquito Anopheles (Cellia) sundaicus (Rodenwaldt) sensu lato. To test this biogeographical hypothesis, we inferred the population history and assessed gene flow of A. sundaicus s.l. sampled from 18 populations across its pan‐Asian species range, using sequences from mitochondrial cytochrome c oxidase subunit 1 (CO1), the internal transcribed spacer 2 (ITS2) and the mannose phosphate isomerase (Mpi) gene. A hypothesis of ecological speciation for A. sundaicus involving divergent adaptation to brackish and freshwater larval habitats was also previously proposed, based on a deficiency of heterozygotes for Mpi allozyme alleles in sympatry. This hypothesis was not supported by Mpi sequence data, which exhibited no fixed differences between brackish and freshwater larval habitats. Mpi and CO1 supported the presence of up to eight genetically distinct population groupings. Counter to the hypothesis of three allopatric species, divergence was often no greater between Borneo, Sumatra/Java and the Southeast Asian mainland than it was between genetic groupings within these landmasses. An isolation‐with‐migration (IM) model indicates recurrent gene flow between the current major landmasses. Such gene flow would have been possible during glacial periods when the current landmasses merged, presenting opportunities for dispersal along expanding and contracting coastlines. Consequently, Pleistocene climatic variation has proved a homogenizing, rather than diversifying, force for A. sundaicus diversity.     

Whole‐genome resequencing reveals the pleistocene temporal dynamics of Branchiostoma belcheri and Branchiostoma floridae populations

Global climatic fluctuations governed the ancestral demographic histories of species and contributed to place the current population status into a more extensive ecological and evolutionary context. Genetic variations will leave unambiguous signatures in the patterns of intraspecific genetic variation in extant species since the genome of each individual is an imperfect mosaic of the ancestral genomes. Here, we report the genome sequences of 20 Branchiostoma individuals by whole‐genome resequencing strategy. We detected over 140 million genomic variations for each Branchiostoma individual. In particular, we applied the pairwise sequentially Markovian coalescent (PSMC) method to estimate the trajectories of changes in the effective population size (N_e) of Branchiostoma population during the Pleistocene. We evaluated the threshold of sequencing depth for proper inference of demographic histories using PSMC was ≥25×. The PSMC results highlight the role of historical global climatic fluctuations in the long‐term population dynamics of Branchiostoma. The inferred ancestral N_e of the Branchiostoma belcheri populations from Zhanjiang and Xiamen (China) seawaters was different in amplitude before the first (mutation rate = 3 × 10^?9) or third glaciation (mutation rate = 9 × 10^?9) of the Pleistocene, indicating that the two populations most probably started to evolve in isolation in their respective seas after the first or third glaciation of the Pleistocene. A pronounced population bottleneck coinciding with the last glacial maximum was observed in all Branchiostoma individuals, followed by a population expansion occurred during the late Pleistocene. Species that have experienced long‐term declines may be especially vulnerable to recent anthropogenic activities. Recently, the industrial pollution and the exploitation of sea sand have destroyed the harmonious living environment of amphioxus species. In the future, we need to protect the habitat of Branchiostoma and make full use of these detected genetic variations to facilitate the functional study of Branchiostoma for adaptation to local environments.     

Shaped by uneven Pleistocene climate: mitochondrial phylogeographic pattern and population history of white wagtail Motacilla alba (Aves: Passeriformes)

We studied the phylogeography and population history of the white wagtail Motacilla alba, which has a vast breeding range, covering areas with different Pleistocene climatic histories. The mitochondrial NADH dehydrogenase subunit II gene (ND2) and Control Region (CR) were analyzed for 273 individuals from 45 localities. Our data comprised all nine subspecies of white wagtail. Four primary clades were inferred (M, N, SW and SE), with indications of M. grandis being nested within M. alba. The oldest split was between two haplotypes from the endemic Moroccan M. a. subpersonata (clade M) and the others, at 0.63–0.96 Mya; other divergences were at 0.31–0.38 Mya. The entire differentiation falls within the part of the Pleistocene characterized by Milankovitch cycles of large amplitudes and durations. Clade N was distributed across the northern Palearctic; clade SW in southwestern Asia plus the British Isles and was predicted by Ecological niche models (ENMs) to occur also in central and south Europe; and clade SE was distributed in central and east Asia. The deep divergence within M. a. subpersonata may reflect retention of ancestral haplotypes. Regional differences in historical climates have had different impacts on different populations: clade N expanded after the last glacial maximum (LGM), whereas milder Pleistocene climate of east Asia allowed clade SE a longer expansion time (since MIS 5); clade SW expanded over a similarly long time as clade SE, which is untypical for European species. ENMs supported these conclusions in that the northern part of the Eurasian continent was unsuitable during the LGM, whereas southern parts remained suitable. The recent divergences and poor structure in the mitochondrial tree contrasts strongly with the pronounced, well defined phenotypical differentiation, indicating extremely fast plumage divergence.     

Climatic stability in the Brazilian Cerrado: implications for biogeographical connections of South American savannas,species richness and conservation in a biodiversity hotspot

Aim To investigate the historical distribution of the Cerrado across Quaternary climatic fluctuations and to generate historical stability maps to test: (1) whether the ‘historical climate’ stability hypothesis explains squamate reptile richness in the Cerrado; and (2) the hypothesis of Pleistocene connections between savannas located north and south of Amazonia. Location The Cerrado, a savanna biome and a global biodiversity hotspot distributed mainly in central Brazil. Methods We generated occurrence datasets from 1000 presence points randomly selected from the entire distribution of the Cerrado, as determined by two spatial definitions. We modelled the potential Cerrado distribution by implementing a maximum‐entropy machine‐learning algorithm across four time projections: current, mid‐Holocene (6 ka), Last Glacial Maximum (LGM, 21 ka) and Last Interglacial (LIG, 120 ka). We generated historical stability maps (refugial areas) by overlapping presence/absence projections of all scenarios, and checked consistencies with qualitative comparisons with available fossil pollen records. We built a spatially explicit simultaneous autoregressive model to explore the relationship between current climate, climatic stability, and squamate species richness. Results Models predicted the LGM and LIG as the periods of narrowest and widest Cerrado distributions, respectively, and were largely corroborated by palynological evidence. We found evidence for two savanna corridors (eastern coastal during the LIG, and Andean during the LGM) and predicted a large refugial area in the north‐eastern Cerrado (Serra Geral de Goiás refugium). Variables related to climatic stability predicted squamate richness better than present climatic variables did. Main conclusions Our results indicate that Bolivian savannas should be included within the Cerrado range and that the Cerrado’s biogeographical counterparts are not Chaco and Caatinga but rather the disjunct savannas of the Guyana shield plateaus. Climatic stability is a good predictor of Cerrado squamate richness, and our stability maps could be used in future studies to test diversity patterns and genetic signatures of different taxonomic groups and as a higher‐order landscape biodiversity surrogate for conservation planning.     

Diverse processes shape deep phylogeographical divergence in Cobitis sinensis (Teleostei: Cobitidae) in East Asia

Little has been known about the impacts of past vicariance events on the phylogeography and population structure of freshwater fishes in East Asia. The aims of this study are to assess the genetic variability with extensive sampling throughout the range of Chinese spiny loach, Cobitis sinensis, and to infer the genetic structure and evolutionary history of populations. Cobitis sinensis in China may have initiated from two ancestral populations, namely Yangtze and Pearl Rivers, which diverged about 7.24 MYA likely due to drainage systems alteration. In the phylogroup I, a southward dispersal event occurred from East China (Yangtze River) to south ZheMin and Hainan subregions, followed by eastward dispersal from ZheMin to south Taiwan. In the phylogroup II, eastward colonization took place from Pearl River to north Taiwan in the late Pliocene, coupled with loss of genetic diversity in the island populations. This study showed that Cenozoic tectonic movements and climatic and sea‐level fluctuations may have shaped the genetic structure of C. sinensis in concert. Highly diverged mtDNA sequences suggest existence of cryptic species in morphospecies C. sinensis.     

A NEW SPECIES OF DICERORHINUS (RHINOCEROTIDAE) FROM THE PLIO‐PLEISTOCENE OF MYANMAR

Abstract: A skull and mandible of the new species Dicerorhinus gwebinensis sp. nov. of Rhinocerotidae (Mammalia, Perissodactyla) is described. The material is collected from the upper part of the Irrawaddy sediments (Plio‐Pleistocene) in central Myanmar. D. gwebinensis sp. nov. is morphologically more similar to the extant species D. sumatrensis (Sumatran rhinoceros) than to other species of the genus but differs from D. sumatrensis in having the comparatively shorter nasal, the more concave dorsal profile of the skull, the more elevated occiput and presence of molar crista in M3/. This is the first discovery of Dicerorhinus in the upper Miocene to lower Pleistocene of the Indian subcontinent and Mainland Southeast Asia, and fills the chronological and geographical gap of this lineage in Asia. The Dicerorhinus clade probably migrated into Southeast Asia from East Asia by the Pliocene or early Pleistocene. This hypothesis is supported by the scarcity or absence of this clade in the Neogene mammalian fauna of the Indian Subcontinent.     

Genetic patterns and changes in availability of suitable habitat support a colonisation history of a North American perennial plant

• Climatic fluctuations during the Pleistocene influenced the geographical distribution of plant species across the southern region of California. Following an integrative approach, we combined genetic data analysis with Environmental Niche Models (ENMs) to assess the historical range expansion of Yucca schidigera, a long‐lived desert perennial native of the Baja California Peninsula.
• We genotyped 240 individuals with seven nuclear microsatellite to investigate genetic diversity distribution across 13 populations. Indeed, we used Environmental Niche Models to examine the changes on the distribution of suitable climatic conditions for this species during the LIG (~120 ka), LGM (~22 ka) and Mid Holocene (~6 ka).
• We detected high genetic diversity across Y. schidigera populations (A_R = 9.94 ± 0.38 SE; H_exp = 0.791 ± 0.011 SE) with genetic variation decreasing significantly with latitude (allelic richness: R ² = 0.38, P = 0.023; expected heterocigosity: R² = 0.32, P = 0.042). We observed low, but significant genetic differentiation (F_ST = 0.0678; P < 0.001) which was consistent with the parapatric distribution of the three genetic groupings detected by the Bayesian clustering algorithm. The ENMs suggest that suitable habitat for this species increased since the LGM.
• Our results support a range expansion of Y. schidigera across northwestern Baja California during the late Quaternary. Genetic data suggest that colonization of the current distribution followed a southward directionality as suitable climatic conditions became widely available in this region. High genetic variation across our sample suggests large historic effective population sizes for this section of the geographical range.

     

Incomplete lineage sorting or secondary admixture: disentangling historical divergence from recent gene flow in the Vinous‐throated parrotbill (Paradoxornis webbianus)

Although Pleistocene glaciations had a major impact on the population genetic patterns of many species in North America and Europe, it remains unclear how these climatic fluctuations contributed to species diversification in East Asia. One reason for this is the difficulty of distinguishing genetic admixture following secondary contact from incomplete lineage sorting, both of which can generate similar patterns of genetic variation. Using a combination of multilocus analyses and coalescent simulation, we explore how these two processes occurred in the Pleistocene evolutionary history of a widespread East Asian bird, the Vinous‐throated parrotbill, Paradoxornis webbianus. Maximum likelihood (ML) tree identified two major mitochondrial lineages, which are geographically separated in most parts of its range, but are sympatric at a few sampling sites. NJ tree and Structure analysis of microsatellite data set revealed an extensive level of admixture and little population structure, suggesting recent admixture between two formerly separated groups. Networks from nuclear DNA data sets, however, did not indicate any geographically isolated groups but rather a panmictic population, thus support incomplete lineage sorting. By using coalescent simulation approaches, we show that both processes did occur, although at different temporal scales. During the Pleistocene glaciations, probably around 0.1–0.5 Ma (the Marine Isotope Stage 6, MIS6), P. webbianus contracted into two separate refugia, and subsequently accumulated genetic divergence. During the interglacial MIS5, the species expanded into previously glaciated areas allowing the once separated groups to come into contact and become admixed. Taken together, our results indicate the current genetic variation within P. webbianus is a combination pattern of widespread distribution in pre‐Pleistocene, then contraction and fragmentation into separated refugia during glacial advance, followed by recently postglacial expansion and admixture.     

Isolation by instability: Historical climate change shapes population structure and genomic divergence of treefrogs in the Neotropical Cerrado savanna

Although the impact of Pleistocene glacial cycles on the diversification of the tropical biota was once dismissed, increasing evidence suggests that Pleistocene climatic fluctuations greatly affected the distribution and population divergence of tropical organisms. Landscape genomic analyses coupled with paleoclimatic distribution models provide a powerful way to understand the consequences of past climate changes on the present‐day tropical biota. Using genome‐wide SNP data and mitochondrial DNA, combined with projections of the species distribution across the late Quaternary until the present, we evaluate the effect of paleoclimatic shifts on the genetic structure and population differentiation of Hypsiboas lundii, a treefrog endemic to the South American Cerrado savanna. Our results show a recent and strong genetic divergence in H. lundii across the Cerrado landscape, yielding four genetic clusters that do not seem congruent with any current physical barrier to gene flow. Isolation by distance (IBD) explains some of the population differentiation, but we also find strong support for past climate changes promoting range shifts and structuring populations even in the presence of IBD. Post‐Pleistocene population persistence in four main areas of historical stable climate in the Cerrado seems to have played a major role establishing the present genetic structure of this treefrog. This pattern is consistent with a model of reduced gene flow in areas with high climatic instability promoting isolation of populations, defined here as “isolation by instability,” highlighting the effects of Pleistocene climatic fluctuations structuring populations in tropical savannas.