Dynamic diversification history with rate upshifts in Holarctic bell‐flowers (Campanula and allies)

Campanula s.l. is one of the most speciose flowering plant lineages of the Holarctic (ca. 600 species). In the present study we sequenced three regions of the plastid genome (petD, rpl16 and trnK/matK) across a broad sample of Campanula s.l., which markedly improved phylogenetic resolution and statistical support compared to previous studies. Based on this robust phylogenetic hypothesis we estimated divergence times using BEAST, diversification rate shifts using Bayesian Analysis of Macroevolutionary Mixture (BAMM) and TreePar, and ancestral ranges using Biogeography with Bayesian (and likelihood) Evolutionary Analyses in R. Campanula s.l. is estimated to have originated during the Early Eocene but the major diversification events occurred between the Late Oligocene and Middle Miocene. Two upward diversification rate shifts were revealed by BAMM, specific to the crown nodes of two Campanula clades: CAM17, a mostly South European‐SW Asian lineage originating during the Middle Miocene and containing nearly half of all known Campanula species; and CAM15B, a SW Asian–Sino‐Himalayan lineage of nine species originating in the early Pleistocene. The dynamic diversification history of Campanula and the inferred rate shifts are discussed in a geo‐historical context.     

Plastid genomes reveal recurrent formation of allopolyploid Fragaria

Premise of the Study

Recurrent formation of polyploid taxa is a common observation in many plant groups. Haploid, cytoplasmic genomes like the plastid genome can be used to overcome the problem of homeologous genes and recombination in polyploid taxa. Fragaria (Rosaceae) contains several octo‐ and decaploid species. We use plastome sequences to infer the plastid ancestry of these taxa with special focus on the decaploid Fragaria cascadensis.

Methods

We used genome skimming of 96 polyploid Fragaria samples on a single Illumina HiSeq 3000 lane to obtain whole plastome sequences. These sequences were used for phylogenetic reconstructions and dating analyses. Ploidy of all samples was inferred with flow cytometry, and plastid inheritance was examined in a controlled cross of F. cascadensis.

Key Results

The plastid genome phylogeny shows that only the octoploid F. chiloensis is monophyletic, all other polyploid taxa were supported to be para‐ or polyphyletic. The decaploid Fragaria cascadensis has biparental plastid inheritance and four different plastid donors. Diversification of the F. cascadensis clades occurred in the last 230,000 years. The southern part of its distribution range harbors considerably higher genetic diversity, suggestive of a potential refugium.

Conclusions

Fragaria cascadensis had at least four independent origins from parents with different plastomes. In contrast, para‐ and polyphyletic taxa of the octoploid Fragaria species are best explained by incomplete lineage sorting and/or hybridization. Biogeographic patterns in F. cascadensis are probably a result of range shift during the last glacial maximum.     

Climate,physiological tolerance and sex‐biased dispersal shape genetic structure of Neotropical orchid bees

Understanding the impact of past climatic events on the demographic history of extant species is critical for predicting species' responses to future climate change. Palaeoclimatic instability is a major mechanism of lineage diversification in taxa with low dispersal and small geographical ranges in tropical ecosystems. However, the impact of these climatic events remains questionable for the diversification of species with high levels of gene flow and large geographical distributions. In this study, we investigate the impact of Pleistocene climate change on three Neotropical orchid bee species (Eulaema bombiformis, E. meriana and E. cingulata) with transcontinental distributions and different physiological tolerances. We first generated ecological niche models to identify species‐specific climatically stable areas during Pleistocene climatic oscillations. Using a combination of mitochondrial and nuclear markers, we inferred calibrated phylogenies and estimated historical demographic parameters to reconstruct the phylogeographical history of each species. Our results indicate species with narrower physiological tolerance experienced less suitable habitat during glaciations and currently exhibit strong population structure in the mitochondrial genome. However, nuclear markers with low and high mutation rates show lack of association with geography. These results combined with lower migration rate estimates from the mitochondrial than the nuclear genome suggest male‐biased dispersal. We conclude that despite large effective population sizes and capacity for long‐distance dispersal, climatic instability is an important mechanism of maternal lineage diversification in orchid bees. Thus, these Neotropical pollinators are susceptible to disruption of genetic connectivity in the event of large‐scale climatic changes.     

Phylogeographic patterns, genetic affinities and morphological differentiation between Epipactis helleborine and related lineages in a Mediterranean glacial refugium

Background and Aims

In the Mediterranean basin, the Italian peninsula has been suggested to be one of the most important glacial refugia for temperate tree species. The orchid genus Epipactis is widely represented in the Italian peninsula by widespread species and several endemic, localized taxa, including selfing and outcrossing taxa. Here the phylogenetic and phylogeographic relationships in a group of closely related taxa in Epipactis are investigated with the aim of understanding the role of this refugial area for cladogenesis and speciation in herbaceous species, such as terrestrial orchids.

Methods

Ribosomal DNA (rDNA) was employed to assess phylogenetic relationships, and plastid sequence variation in the rbcL–accD spacer was used to reveal phylogeographic patterns among plastid haplotypes using a parsimony network.

Key Results

Low genetic variation and shared ribotypes were detected in rDNA, whereas high levels of sequence variation and a strong phylogeographic structure were found in the examined plastid region. The parsimony plastid haplotype network identified two main haplotype groups, one including E. atrorubens/microphylla/muelleri/leptochila and the other including all accessions of E. helleborine and several localized and endemic taxa, with a combination of widespread and rare haplotypes detected across the Italian peninsula. A greater genetic divergence separated the Italian and other European accessions of E. helleborine.

Conclusions

Phylogenetic and phylogeographic patterns support a working hypothesis in which the Italian peninsula has only recently been colonized by Epipactis, probably during the most recent phase of the Quaternary age and, nevertheless, it acted as a remarkable centre of diversification for this orchid lineage. Changes in pollination strategy and recurrent shifts in mating system (from allogamy to autogamy) could have represented the mechanism promoting this rapid diversification and the observed high taxonomic complexity detected in the E. helleborine species complex.     

Bayesian Poisson tree processes and multispecies coalescent models shed new light on the diversification of Nawab butterflies in the Solomon Islands (Nymphalidae,Charaxinae, Polyura)

Butterflies of the genus Polyura form a widespread tropical group distributed from Pakistan to Fiji. The rare endemic Polyura epigenes Godman & Salvin, 1888 from the Solomon Islands archipelago represents a case of marked island polymorphism. We sequenced museum specimens of this species across its geographic range to study the phylogeography and genetic differentiation of populations in the archipelago. We used the Bayesian Poisson tree processes and multispecies coalescent models, to study species boundaries. We also estimated divergence times to investigate the biogeographic history of populations. Our molecular species delimitation and nuclear DNA network analyses unambiguously indicate that Malaita populations form an independent metapopulation lineage, as defined in the generalized lineage concept. This lineage, previously ranked as a subspecies, is raised to species rank under the name Polyura bicolor Turlin & Sato, 1995  stat. nov. Divergence time estimates suggest that this lineage split from its sister taxon in the late Pleistocene. At this time, the bathymetric isolation of Malaita from the rest of the archipelago probably prevented gene flow during periods of lower sea level, thereby fostering allopatric speciation. The combination of molecular species delimitation methods, morphological comparisons, and divergence time estimation is useful to study lineage diversification across intricate geographic regions.     

Lineage-specific late pleistocene expansion of an endemic subtropical gossamer-wing damselfly, Euphaea formosa, in Taiwan

Background  

Pleistocene glacial oscillations have significantly affected the historical population dynamics of temperate taxa. However, the general effects of recent climatic changes on the evolutionary history and genetic structure of extant subtropical species remain poorly understood. In the present study, phylogeographic and historical demographic analyses based on mitochondrial and nuclear DNA sequences were used. The aim was to investigate whether Pleistocene climatic cycles, paleo-drainages or mountain vicariance of Taiwan shaped the evolutionary diversification of a subtropical gossamer-wing damselfly, Euphaea formosa.     

Systematics and evolution of the Pan‐Alcidae (Aves,Charadriiformes)

Puffins, auks and their allies in the wing‐propelled diving seabird clade Pan‐Alcidae (Charadriiformes) have been proposed to be key pelagic indicators of faunal shifts in Northern Hemisphere oceans. However, most previous phylogenetic analyses of the clade have focused only on the 23 extant alcid species. Here we undertake a combined phylogenetic analysis of all previously published molecular sequence data (～ 12 kb) and morphological data (n = 353 characters) with dense species level sampling that also includes 28 extinct taxa. We present a new estimate of the patterns of diversification in the clade based on divergence time estimates that include a previously vetted set of twelve fossil calibrations. The resultant time trees are also used in the evaluation of previously hypothesized paleoclimatic drivers of pan‐alcid evolution. Our divergence dating results estimate the split of Alcidae from its sister taxon Stercorariidae during the late Eocene (～ 35 Ma), an evolutionary hypothesis for clade origination that agrees with the fossil record and that does not require the inference of extensive ghost lineages. The extant dovekie Alle alle is identified as the sole extant member of a clade including four extinct Miocene species. Furthermore, whereas an Uria + Alle clade has been previously recovered from molecular analyses, the extinct diversity of closely related Miocepphus species yields morphological support for this clade. Our results suggest that extant alcid diversity is a function of Miocene diversification and differential extinction at the Pliocene–Pleistocene boundary. The relative timing of the Middle Miocene climatic optimum and the Pliocene–Pleistocene climatic transition and major diversification and extinction events in Pan‐Alcidae, respectively, are consistent with a potential link between major paleoclimatic events and pan‐alcid cladogenesis.     

Contrasting phylogeography of two Western Palaearctic fish parasites despite similar life cycles

Aim

We used comparative phylogeography of two intestinal parasites of freshwater fish to test whether similarity in life cycle translates into concordant phylogeographical history. The thorny‐headed worms Pomphorhynchus laevis and P. tereticollis (Acanthocephala) were formerly considered as a single species with a broad geographical and host range within the Western Palaearctic.

Location

Central and eastern parts of Northern Mediterranean area, Western and Central Europe, Ponto‐Caspian Europe.

Methods

A mitochondrial marker (COI) was sequenced for 111 P. laevis and 50 P. tereticollis individuals and nuclear ITS1 and ITS2 sequences were obtained for 37 P. laevis and 21 P. tereticollis. Genetic divergence, phylogenetic relationships and divergence time were estimated for various lineages within each species, and their phylogeographical patterns were compared to known palaeogeographical events in Western Palaearctic. Biogeographical histories of each species were inferred.

Results

The two species show very different phylogeographical patterns. Five lineages were identified in P. laevis, partially matching several major biogeographical regions defined in the European riverine fish fauna. The early stages of P. laevis diversification occurred in the peri‐Mediterranean area, during the Late Miocene. Subsequent expansion across Western Europe and Russia was shaped by dispersal and vicariant events, from Middle Pliocene to Middle Pleistocene. By contrast, P. tereticollis has differentiated more recently within the Western and Central parts of Europe, and shows weak geographical and genetic structuring.

Conclusion

Our study highlights weak to moderate similarity in the phylogeographical pattern of these acanthocephalan parasites compared to their amphipod and fish hosts. The observed differences in the timing of dispersion and migration routes taken may reflect the use of a range of final hosts with different ecologies and dispersal capabilities. By using a group underrepresented in phylogeographical studies, our study is a valuable contribution to revealing the biogeography of host–parasite interactions in continental freshwaters.     

Cryptic species and co‐diversification in sand scorpions from the Karakum and Kyzylkum deserts of Central Asia

Little is known about species diversification within the deserts of Central Asia. For example, the degree of lineage divergence and timing of population differentiation, as well as potential biogeographic barriers driving diversification, are nearly unknown. Here, we analysed a multi‐locus data set for a widespread sand scorpion (Mesobuthus gorelovi) to evaluate cryptic species diversity and phylogeographic patterns across the Karakum and Kyzylkum deserts. We also combined these data with previously published sequence data to test for a signal of co‐diversification. A consensus species delimitation approach indicated that the widespread M. gorelovi is likely composed of up to five distinct species that began to diversify at the Miocene–Pliocene boundary. We observed shared patterns of lineage divergence across the Amu Darya River region in three scorpion taxa and found support for a shared history of assemblage diversification across this biogeographic barrier. Thus, major river systems appear to facilitate diversification among desert scorpions.     

Evolutionary drivers of phylogeographical diversity in the highlands of Mexico: a case study of the Crotalus triseriatus species group of montane rattlesnakes

Aim To assess the genealogical relationships of widespread montane rattlesnakes in the Crotalus triseriatus species group and to clarify the role of Late Neogene mountain building and Pleistocene pine–oak forest fragmentation in driving the diversification of Mexican highland taxa. Location Highlands of mainland Mexico and the south‐western United States (Texas, New Mexico, and Arizona). Methods A synthesis of inferences was used to address several associated questions about the biogeography of the Mexican highlands and the evolutionary drivers of phylogeographical diversity in co‐distributed taxa. We combined extensive range‐wide sampling (130 individuals representing five putative species) and mixed‐model phylogenetic analyses of 2408 base pairs of mitochondrial DNA to estimate genealogical relationships and divergence times within the C. triseriatus species group. We then assessed the tempo of diversification using a maximum likelihood framework based on the birth–death process. Estimated times of divergences provided a probabilistic temporal component and questioned whether diversification rates have remained constant or varied over time. Finally, we looked for phylogeographical patterns in other co‐distributed taxa. Results We identified eight major lineages within the C. triseriatus group, and inferred strong correspondence between maternal and geographic history within most lineages. At least one cryptic species was detected. Relationships among lineages were generally congruent with previous molecular studies, with differences largely attributable to our expanded taxonomic and geographic sampling. Estimated divergences between most major lineages occurred in the Late Miocene and Pliocene. Phylogeographical structure within each lineage appeared to have been generated primarily during the Pleistocene. Although the scale of genetic diversity recognized affected estimated rates of diversification, rates appeared to have been constant through time. Main conclusions The biogeographical history of the C. triseriatus group implies a dynamic history for the highlands of Mexico. The Neogene formation of the Transvolcanic Belt appears responsible for structuring geographic diversity among major lineages. Pleistocene glacial–interglacial climatic cycles and resultant expansions and contractions of the Mexican pine–oak forest appear to have driven widespread divergences within lineages. Climatic change, paired with the complex topography of Mexico, probably produced a myriad of species‐specific responses in co‐distributed Mexican highland taxa. The high degree of genetic differentiation recovered in our study and others suggests that the Mexican highlands may contain considerably more diversity than currently recognized.     

Phylogeography of Middle American gophersnakes: mixed responses to biogeographical barriers across the Mexican Transition Zone

Aim We used inferences of phylogeographical structure and estimates of divergence times for three species of gophersnakes (Colubridae: Pituophis) distributed across the Mexican Transition Zone (MTZ) to evaluate the postulated association of three Neogene geological events (marine seaway inundation of the Isthmus of Tehuantepec, formation of the Transvolcanic Belt across central Mexico, and secondary uplifting of the Sierra Madre Occidental) and of Pleistocene climate change with inter‐ and intraspecific diversification. Location Mexico, Guatemala, and the western United States. Methods We combined range‐wide sampling (67 individuals representing three putative species distributed across northern Middle America and western North America) and phylogenetic analyses of 1637 base pairs of mitochondrial DNA to estimate genealogical relationships and divergence times. The hypothesized concordance of inferred gene trees with geological histories was assessed using topology tests. Results We identified three major lineages of Middle American gophersnakes, and strong phylogeographical structure within each lineage. Gene trees were statistically congruent with hypothesized geological histories for two of the three postulated geological events. Estimated divergence dates and the geographical distribution of genetic variation further support mixed responses to these geological events. Considerable phylogeographical structure appears to have been generated during the Pleistocene. Main conclusions Phylogenetic and phylogeographical structure in gophersnakes distributed across northern Middle America and western North America highlights the influence of both Neogene vicariance events and Pleistocene climate change in shaping genetic diversity in this region. Despite the presence of two major geographical barriers in southern Mexico, extreme geological and environmental heterogeneity in this area may have differentially structured genetic diversity in highland taxa. To the north, co‐distributed taxa may display a more predictable pattern of diversification across the warm desert regions. Future studies should incorporate nuclear data to disentangle inferred lineage boundaries and further elucidate patterns of mitochondrial introgression.     

Speciation with gene flow in whiptail lizards from a Neotropical xeric biome

Two main hypotheses have been proposed to explain the diversification of the Caatinga biota. The riverine barrier hypothesis (RBH) claims that the São Francisco River (SFR) is a major biogeographic barrier to gene flow. The Pleistocene climatic fluctuation hypothesis (PCH) states that gene flow, geographic genetic structure and demographic signatures on endemic Caatinga taxa were influenced by Quaternary climate fluctuation cycles. Herein, we analyse genetic diversity and structure, phylogeographic history, and diversification of a widespread Caatinga lizard (Cnemidophorus ocellifer) based on large geographical sampling for multiple loci to test the predictions derived from the RBH and PCH. We inferred two well‐delimited lineages (Northeast and Southwest) that have diverged along the Cerrado–Caatinga border during the Mid‐Late Miocene (6–14 Ma) despite the presence of gene flow. We reject both major hypotheses proposed to explain diversification in the Caatinga. Surprisingly, our results revealed a striking complex diversification pattern where the Northeast lineage originated as a founder effect from a few individuals located along the edge of the Southwest lineage that eventually expanded throughout the Caatinga. The Southwest lineage is more diverse, older and associated with the Cerrado–Caatinga boundaries. Finally, we suggest that C. ocellifer from the Caatinga is composed of two distinct species. Our data support speciation in the presence of gene flow and highlight the role of environmental gradients in the diversification process.     

The timing of diversification within the most divergent parrot clade

The Strigopidae are an ancient parrot (Psittaciformes) family consisting of three extant species placed in two genera (Nestor, Strigops) and restricted to New Zealand. Their evolutionary history is clouded because the timing of divergence events within this family has variously been attributed to Pleistocene climate change or much earlier earth‐historic events. Here we examine new psittaciform DNA sequence data, and combine them with previously published sequences, to shed light on the poorly understood timing of diversification within the Strigopidae. Using calibrations indirectly derived from both psittaciform and non‐psittaciform fossils, our data indicate a Late Pliocene or Early Pleistocene (ca 1.2–3.6 mya) differentiation between the two Nestor species (kea and kaka), possibly in response to shifts in habitat distribution associated with sea level fluctuations. The unique, monotypic, nocturnal and flightless genus Strigops (kakapo) is shown to have diverged from the Nestor lineage probably ca 28–29 mya, coinciding with the potential Oligocene submergence of Zealandia when much of its landmass may have been fragmented into smaller islands, providing a setting for allopatric diversification.     

Complete plastome sequences of Equisetum arvense and Isoetes flaccida: implications for phylogeny and plastid genome evolution of early land plant lineages

Background

Understanding the forces that shaped Neotropical diversity is central issue to explain tropical biodiversity and inform conservation action; yet few studies have examined large, widespread species. Lowland tapir (Tapirus terrrestris, Perissodactyla, Tapiridae) is the largest Neotropical herbivore whose ancestors arrived in South America during the Great American Biotic Interchange. A Pleistocene diversification is inferred for the genus Tapirus from the fossil record, but only two species survived the Pleistocene megafauna extinction. Here, we investigate the history of lowland tapir as revealed by variation at the mitochondrial gene Cytochrome b, compare it to the fossil data, and explore mechanisms that could have shaped the observed structure of current populations.

Results

Separate methodological approaches found mutually exclusive divergence times for lowland tapir, either in the late or in the early Pleistocene, although a late Pleistocene divergence is more in tune with the fossil record. Bayesian analysis favored mountain tapir (T. pinchaque) paraphyly in relation to lowland tapir over reciprocal monophyly, corroborating the inferences from the fossil data these species are sister taxa. A coalescent-based analysis rejected a null hypothesis of allopatric divergence, suggesting a complex history. Based on the geographic distribution of haplotypes we propose (i) a central role for western Amazonia in tapir diversification, with a key role of the ecological gradient along the transition between Andean subcloud forests and Amazon lowland forest, and (ii) that the Amazon river acted as an barrier to gene flow. Finally, the branching patterns and estimates based on nucleotide diversity indicate a population expansion after the Last Glacial Maximum.

Conclusions

This study is the first examining lowland tapir phylogeography. Climatic events at the end of the Pleistocene, parapatric speciation, divergence along the Andean foothill, and role of the Amazon river, have similarly shaped the history of other taxa. Nevertheless further work with additional samples and loci is needed to improve our initial assessment. From a conservation perspective, we did not find a correspondence between genetic structure in lowland tapir and ecogeographic regions proposed to define conservation priorities in the Neotropics. This discrepancy sheds doubt into this scheme's ability to generate effective conservation planning for vagile species.     

Cyto‐nuclear discordance suggests complex evolutionary history in the cave‐dwelling salamander,Eurycea lucifuga

Our understanding of the evolutionary history and ecology of cave‐associated species has been driven historically by studies of morphologically adapted cave‐restricted species. Our understanding of the evolutionary history and ecology of nonrestricted cave species, troglophiles, is limited to a few studies, which present differing accounts of troglophiles’ relationship with the cave habitat, and its impact on population dynamics. Here, we used phylogenetics, demographic statistics, and population genetic methods to study lineage divergence, dates of divergence, and population structure in the Cave Salamander, Eurycea lucifuga, across its range. In order to perform these analyses, we sampled 233 individuals from 49 populations, using sequence data from three gene loci as well as genotyping data from 19 newly designed microsatellite markers. We find, as in many other species studied in a phylogeographic context, discordance between patterns inferred from mitochondrial relationships and those inferred by nuclear markers indicating a complicated evolutionary history in this species. Our results suggest Pleistocene‐based divergence among three main lineages within E. lucifuga corresponding to the western, central, and eastern regions of the range, similar to patterns seen in species separated in multiple refugia during climatic shifts. The conflict between mitochondrial and nuclear patterns is consistent with what we would expect from secondary contact between regional populations following expansion from multiple refugia.     

Geographical isolation versus dispersal: Relictual alpine grasshoppers support a model of interglacial diversification with limited hybridization

Alpine biotas are paradigmatic of the countervailing roles of geographical isolation and dispersal during diversification. In temperate regions, repeated distributional shifts driven by Pleistocene climatic oscillations produced both recurrent pulses of population fragmentation and opportunities for gene flow during range expansions. Here, we test whether a model of divergence in isolation vs. with gene flow is more likely in the diversification of flightless alpine grasshoppers of the genus Podisma from the Iberian Peninsula. The answer to this question can also provide key insights about the pace of evolution. Specifically, if the data fit a divergence in isolation model, this suggests rapid evolution of reproductive isolation. Genomic data confirm a Pleistocene origin of the species complex, and multiple analytical approaches revealed limited asymmetric historical hybridization between two taxa. Genomic-based demographic reconstructions, spatial patterns of genetic structure and range shifts inferred from palaeodistribution modelling suggest severe range contraction accompanied by declines in effective population sizes during interglacials (i.e., contemporary populations confined to sky islands are relicts) and expansions during the coldest stages of the Pleistocene in each taxon. Although limited hybridization during secondary contact leads to phylogenetic uncertainty if gene flow is not accommodated when estimating evolutionary relationships, all species exhibit strong genetic cohesiveness. Our study lends support to the notion that the accumulation of incipient differences during periods of isolation were sufficient to lead to lineage persistence, but also that the demographic changes, dispersal constraints and spatial distribution of the sky islands themselves mediated species diversification in temperate alpine biotas.     

From pampa to puna: Biogeography and diversification of a group of Neotropical obligate grassland birds (Anthus: Motacillidae)

The evolution of Neotropical birds of open landscapes remains largely unstudied. We investigate the diversification and biogeography of a group of Neotropical obligate grassland birds (Anthus: Motacillidae). We use a multilocus phylogeny of 22 taxa of Anthus to test the hypothesis that these birds radiated contemporaneously with the development of grasslands in South America. We employ the R package DDD to analyze the dynamics of Anthus diversification across time in Neotropical grasslands, explicitly testing for shifts in dynamics associated with the Miocene development of grasslands, the putative Pleistocene expansion of arid lowland biomes, and Pleistocene sundering of Andean highland grasslands. A lineage‐through‐time plot revealed increases in the number of lineages, and DDD detected shifts to a higher clade‐level carrying capacity during the late Miocene, indicating an early burst of diversification associated with grassland colonization. However, we could not corroborate the shift using power analysis, probably reflecting the small number of tips in our tree. We found evidence of a divergence at ~1 Mya between northern and southern Amazonian populations of Anthus lutescens, countering Haffer's idea of Pleistocene expansion of open biomes in the Amazon Basin. We used BioGeoBears to investigate ancestral areas and directionality of colonization of Neotropical grasslands. Members of the genus diversified into, out of, and within the Andes, within‐Andean diversification being mostly Pleistocene in origin.     

Temporal,spatial and ecological dynamics of speciation among amphi‐Beringian small mammals

Aim

Quaternary climate cycles played an important role in promoting diversification across the Northern Hemisphere, although details of the mechanisms driving evolutionary change are still poorly resolved. In a comparative phylogeographical framework, we investigate temporal, spatial and ecological components of evolution within a suite of Holarctic small mammals. We test a hypothesis of simultaneous divergence among multiple taxon pairs, investigating time to coalescence and demographic change for each taxon in response to a combination of climate and geography.

Location

Beringia, the nexus of the northern continents.

Methods

We used approximate Bayesian computation methods to test for simultaneous divergence among eight pairs of taxa, using cytochrome b gene sequences. We calculated coalescence times for eastern and western components of each pair and for the combined pairs, and relate dates to Quaternary climatic periodicity and combinations of environmental events and physical barriers. Population growth and expansion statistics were used to test evolutionary responses among taxa, including range shifts, persistence or periodic extirpation. Species distribution models (SDMs) for each taxon were used to predict their geographical ranges during the present interglacial, Last Glacial Maximum and previous interglacial.

Results

Multiple divergence events across Beringia were primarily coincident with extreme glacial cycles of the late Quaternary. Structure within Beringia is spatially consistent with at least three environmental barriers arising at different times: the Kolyma Uplands, Bering Strait and portions of the Bering Isthmus. Levels of divergence varied substantially, indicating evolutionary processes spanning deep and shallow time‐scales. The different demographics among taxa reflect their distinct ecological responses. SDMs predicted regional distributional changes through time and different spatial responses among taxa.

Main conclusions

Beringia predominantly constituted a dispersal corridor during the early Quaternary and a major centre of endemism in the late Quaternary. Coincident with severe glacial cycles, small mammal species were ‘caught’ in Beringia and diversified over multiple climatic phases. Relative genetic differentiation across Beringia appears to be related to ecological differences reflecting a gradual adaptation to Beringian environments through time. Some methodological constraints associated with resolving recent (late Quaternary) isolation events or drawing inferences from a single locus are discussed.