Refugia revisited: individualistic responses of species in space and time

Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

Museomics help resolving the phylogeny of snowfinches (Aves, Passeridae, Montifringilla and allies)

Historical specimens from museum collections provide a valuable source of material also from remote areas or regions of conflict that are not easily accessible to scientists today. With this study, we are providing a taxon-complete phylogeny of snowfinches using historical DNA from whole skins of an endemic species from Afghanistan, the Afghan snowfinch, Pyrgilauda theresae. To resolve the strong conflict between previous phylogenetic hypotheses, we generated novel mitogenome sequences for selected taxa and genome-wide SNP data using ddRAD sequencing for all extant snowfinch species endemic to the Qinghai-Tibet Plateau (QTP) and for an extended intraspecific sampling of the sole Central and Western Palearctic snowfinch species (Montifringilla nivalis). Our phylogenetic reconstructions unanimously refuted the previously suggested paraphyly of genus Pyrgilauda. Misplacement of one species-level taxon (Onychostruthus tazcanowskii) in previous snowfinch phylogenies was undoubtedly inferred from chimeric mitogenomes that included heterospecific sequence information. Furthermore, comparison of novel and previously generated sequence data showed that the presumed sister-group relationship between M. nivalis and the QTP endemic M. henrici was suggested based on flawed taxonomy. Our phylogenetic reconstructions based on genome-wide SNP data and on mitogenomes were largely congruent and supported reciprocal monophyly of genera Montifringilla and Pyrgilauda with monotypic Onychostruthus being sister to the latter. The Afghan endemic P. theresae likely originated from a rather ancient Pliocene out-of-Tibet dispersal probably from a common ancestor with P. ruficollis. Our extended trans-Palearctic sampling for the white-winged snowfinch, M. nivalis, confirmed strong lineage divergence between an Asian and a European clade dated to 1.5 - 2.7 million years ago (mya). Genome-wide SNP data suggested subtle divergence among European samples from the Alps and from the Cantabrian mountains.

Sympatry in a nightingale contact zone has no effect on host-specific blood parasite prevalence and lineage diversity

Parasites are a key driving force behind many ecological and evolutionary processes. Prevalence and diversity of parasites, as well as their effects on hosts, are not uniform across host species. As such, the potential parasite spillover between species can significantly influence outcomes of interspecific interactions. We screened two species of Luscinia nightingales for haemosporidian blood parasites (Plasmodium, Leucocytozoon and Haemoproteus) along an approximately 3000 km transect in Europe, incorporating areas of host distant allopatry, close allopatry and sympatry. We found significant differences in infection rates between the two host species, with common nightingales having much lower parasite prevalence than thrush nightingales (36.7% versus 83.8%). This disparity was mostly driven by Haemoproteus prevalence, which was significantly higher in thrush nightingales while common nightingales had a small, but significantly higher, Plasmodium prevalence. Furthermore, we found no effect of proximity to the contact zone on infection rate in either host species. Despite having lower infection prevalence, common nightingales were infected with a significantly higher diversity of parasite lineages than thrush nightingales, and lineage assemblages differed considerably between the two species, even in sympatry. This pattern was mostly driven by the large diversity of comparatively rare lineages, while the most abundant lineages were shared between the two host species. This suggests that, despite the close evolutionary relationships between the two nightingales, there are significant differences in parasite prevalence and diversity, regardless of the distance from the contact zone. This suggests that spillover of haemosporidian blood parasites is unlikely to contribute towards interspecific interactions in this system.

The Phylogeography of Deciduous Tree Ulmus macrocarpa (Ulmaceae) in Northern China

Disentangling how climate oscillations and geographical events significantly influence plants' genetic architecture and demographic history is a central topic in phytogeography. The deciduous ancient tree species Ulmus macrocarpa is primarily distributed throughout Northern China and has timber and horticultural value. In the current study, we studied the phylogenic architecture and demographical history of U. macrocarpa using chloroplast DNA with ecological niche modeling. The results indicated that the populations' genetic differentiation coefficient (NST) value was significantly greater than the haplotype frequency (GST) (p < 0.05), suggesting that U. macrocarpa had a clear phylogeographical structure. Phylogenetic inference showed that the putative chloroplast haplotypes could be divided into three groups, in which the group Ⅰ was considered to be ancestral. Despite significant genetic differentiation among these groups, gene flow was detected. The common ancestor of all haplotypes was inferred to originate in the middle-late Miocene, followed by the haplotype overwhelming diversification that occurred in the Quaternary. Combined with demography pattern and ecological niche modeling, we speculated that the surrounding areas of Shanxi and Inner Mongolia were potential refugia for U. macrocarpa during the glacial period in Northern China. Our results illuminated the demography pattern of U. macrocarpa and provided clues and references for further population genetics investigations of precious tree species distributed in Northern China.

Environmental and Biogeographic Drivers behind Alpine Plant Thermal Tolerance and Genetic Variation

In alpine ecosystems, elevation broadly functions as a steep thermal gradient, with plant communities exposed to regular fluctuations in hot and cold temperatures. These conditions lead to selective filtering, potentially contributing to species-level variation in thermal tolerance and population-level genetic divergence. Few studies have explored the breadth of alpine plant thermal tolerances across a thermal gradient or the underlying genetic variation thereof. We measured photosystem heat (Tcrit-hot) and cold (Tcrit-cold) thresholds of ten Australian alpine species across elevation gradients and characterised their neutral genetic variation. To reveal the biogeographical drivers of present-day genetic signatures, we also reconstructed temporal changes in habitat suitability across potential distributional ranges. We found intraspecific variation in thermal thresholds, but this was not associated with elevation, nor underpinned by genetic differentiation on a local scale. Instead, regional population differentiation and considerable homozygosity within populations may, in part, be driven by distributional contractions, long-term persistence, and migrations following habitat suitability. Our habitat suitability models suggest that cool-climate-distributed alpine plants may be threatened by a warming climate. Yet, the observed wide thermal tolerances did not reflect this vulnerability. Conservation efforts should seek to understand variations in species-level thermal tolerance across alpine microclimates.

An assessment of methods to combine evolutionary history and conservation: A case study in the Brazilian campo rupestre

Premise

Conservation policies typically focus on biodiversity hotspots. An alternative approach involves analyzing the evolutionary history of lineages in geographic areas along with their threat levels to guide conservation efforts. Mountains exhibit high levels of plant species richness and micro-endemism, and biogeographic studies commonly point to recent and rapid evolutionary radiations in these areas. Using a nearly endemic clade of legumes, our study evaluates conservation prioritization approaches in the campo rupestre, a Neotropical ecosystem associated with mountaintops that is located between two biodiversity hotspots.

Methods

We compared the EDGE and EDGE2 metrics, which combine the evolutionary distinctiveness and the extinction risk of a species in a single value. These metrics are compared with traditional metrics used to assess conservation priority, such as phylogenetic diversity.

Results

The EDGE values reported are lower than those of other studies using this metric, mostly due to the prevalence of threatened species with short phylogenetic branch lengths (low values of evolutionary distinctiveness). Certain areas of campo rupestre with relatively high phylogenetic diversity and EDGE values do not correspond to areas with high species richness, agreeing with previous studies on biodiversity hotspots.

Discussion

Our study highlights the necessity of conservation of the campo rupestres as well as advantages and disadvantages of using EDGE, EDGE2, and phylogenetic diversity for appropriate selection of conservation areas with rapid evolutionary radiations. The selection of the metrics will depend primarily on the life history of the focus group and the data availability, as well as the conservation approach.

Where east meets west: Phylogeography of the high Arctic North American brant goose

Genetic variation in Arctic species is often influenced by vicariance during the Pleistocene, as ice sheets fragmented the landscape and displaced populations to low- and high-latitude refugia. The formation of secondary contact or suture zones during periods of ice sheet retraction has important consequences on genetic diversity by facilitating genetic connectivity between formerly isolated populations. Brant geese (Branta bernicla) are a maritime migratory waterfowl (Anseriformes) species that almost exclusively uses coastal habitats. Within North America, brant geese are characterized by two phenotypically distinct subspecies that utilize disjunct breeding and wintering areas in the northern Pacific and Atlantic. In the Western High Arctic of Canada, brant geese consist of individuals with an intermediate phenotype that are rarely observed nesting outside this region. We examined the genetic structure of brant geese populations from each subspecies and areas consisting of intermediate phenotypes using mitochondrial DNA (mtDNA) control region sequence data and microsatellite loci. We found a strong east-west partition in both marker types consistent with refugial populations. Within subspecies, structure was also observed at mtDNA while microsatellite data suggested the presence of only two distinct genetic clusters. The Western High Arctic (WHA) appears to be a secondary contact zone for both Atlantic and Pacific lineages as mtDNA and nuclear genotypes were assigned to both subspecies, and admixed individuals were observed in this region. The mtDNA sequence data outside WHA suggests no or very restricted intermixing between Atlantic and Pacific wintering populations which is consistent with published banding and telemetry data. Our study indicates that, although brant geese in the WHA are not a genetically distinct lineage, this region may act as a reservoir of genetic diversity and may be an area of high conservation value given the potential of low reproductive output in this species.

The Cabrières Biota (France) provides insights into Ordovician polar ecosystems

Early Palaeozoic sites with soft-tissue preservation are predominantly found in Cambrian rocks and tend to capture past tropical and temperate ecosystems. In this study, we describe the diversity and preservation of the Cabrières Biota, a newly discovered Early Ordovician Lagerstätte from Montagne Noire, southern France. The Cabrières Biota showcases a diverse polar assemblage of both biomineralized and soft-bodied organisms predominantly preserved in iron oxides. Echinoderms are extremely scarce, while sponges and algae are abundantly represented. Non-biomineralized arthropod fragments are also preserved, along with faunal elements reminiscent of Cambrian Burgess Shale-type ecosystems, such as armoured lobopodians. The taxonomic diversity observed in the Cabrières Biota mixes Early Ordovician Lagerstätten taxa with Cambrian forms. By potentially being the closest Lagerstätte to the South Pole, the Cabrières Biota probably served as a biotic refuge amid the high-water temperatures of the Early Ordovician, and shows comparable ecological structuring to modern polar communities.

The opposed forces of differentiation and admixture across glacial cycles in the butterfly Aglais urticae

Glacial cycles lead to periodic population interbreeding and isolation in warm-adapted species, which impact genetic structure and evolution. However, the effects of these processes on highly mobile and more cold-tolerant species are not well understood. This study aims to shed light on the phylogeographic history of Aglais urticae, a butterfly species with considerable dispersal ability, and a wide Palearctic distribution reaching the Arctic. Through the analysis of genomic data, four main genetic lineages are identified: European, Sierra Nevada, Sicily/Calabria/Peloponnese, and Eastern. The results indicate that the Sardo-Corsican endemic taxon ichnusa is a distinct species. The split between the relict lineages in southern Europe and the main European lineage is estimated to have happened 400-450 thousand years ago, with admixture observed during the Quaternary glacial cycles, and still ongoing, albeit to a much smaller extent. These results suggest that these lineages may be better treated as subspecific parapatric taxa. Ecological niche modelling supported the existence of both Mediterranean and extra-Mediterranean refugia during the glacial periods, with the main one located on the Atlantic coast. Nevertheless, gene flow between populations was possible, indicating that both differentiation and admixture have acted continuously across glacial cycles in this cold-tolerant butterfly, generally balancing each other but producing differentiated lineages in the southern peninsulas. We conclude that the population dynamics and the processes shaping the population genetic structure of cold-adapted species during the Quaternary ice ages may be different than those classically accepted for warm-adapted species.

Demographic patterns of two related desert shrubs with overlapping distributions in response to past climate changes

Numerous studies have revealed that past geological events and climatic fluctuations had profoundly affected the genetic structure and demographic patterns of species. However, related species with overlapping ranges may have responded to such environmental changes in different ways. In this study, we compared the genetic structure and population dynamics of two typical desert shrubs with overlapping distributions in northern China, Nitraria tangutorum and Nitraria sphaerocarpa, based on chloroplast DNA (cpDNA) variations and species distribution models. We sequenced two cpDNA fragments (trnH-trnA and atpH-atpI) in 633 individuals sampled from 52 natural populations. Twenty-four chlorotypes, including eight rare chlorotypes, were identified, and a single dominant haplotype (H4) widely occurred in the entire geographical ranges of the two species. There were also a few distinctive chlorotypes fixed in different geographical regions. Population structure analyses suggested that the two species had significantly different levels of total genetic diversity and interpopulation differentiation, which was highly likely correlated with the special habitat preferences of the two species. A clear phylogeographic structure was identified to exist among populations of N. sphaerocarpa, but not exist for N. tangutorum. The neutral tests, together with the distribution of pairwise differences revealed that N. tangutorum experienced a sudden demographic expansion, and its expansion approximately occurred between 21 and 7 Kya before present, while a rapid range expansion was not identified for N. sphaerocarpa. The ecological niche modeling (ENM) analysis indicated that the potential ranges of two species apparently fluctuated during the past and present periods, with obvious contraction in the Last Glacial Maximum (LGM) and recolonization in the present, respectively, comparing to the Last Interglacial (LIG). These findings suggest that the two species extensively occurred in the Northwest of China before the Quaternary, and the current populations of them originated from a few separated glacial refugia following their habitat fragmentation in the Quarternary. Our results provide new insights on the impact of past geological and climatic fluctuations on the population dynamics of desert plants in northwestern China, and further enforce the hypothesis that there were several independent glacial refugia for these species during the Quaternary glaciations.

Ancient reindeer mitogenomes reveal island-hopping colonisation of the Arctic archipelagos

Climate warming at the end of the last glacial period had profound effects on the distribution of cold-adapted species. As their range shifted towards northern latitudes, they were able to colonise previously glaciated areas, including remote Arctic islands. However, there is still uncertainty about the routes and timing of colonisation. At the end of the last ice age, reindeer/caribou (Rangifer tarandus) expanded to the Holarctic region and colonised the archipelagos of Svalbard and Franz Josef Land. Earlier studies have proposed two possible colonisation routes, either from the Eurasian mainland or from Canada via Greenland. Here, we used 174 ancient, historical and modern mitogenomes to reconstruct the phylogeny of reindeer across its whole range and to infer the colonisation route of the Arctic islands. Our data shows a close affinity among Svalbard, Franz Josef Land and Novaya Zemlya reindeer. We also found tentative evidence for positive selection in the mitochondrial gene ND4, which is possibly associated with increased heat production. Our results thus support a colonisation of the Eurasian Arctic archipelagos from the Eurasian mainland and provide some insights into the evolutionary history and adaptation of the species to its High Arctic habitat.

Cliff-edge forests: Xerothermic hotspots of local biodiversity and models for future climate change

Cliffs are remarkable environments that enable the existence of microclimates. These small, isolated sites, decoupled from the regional macroclimate, play a significant role in maintaining species biodiversity, particularly in topographically homogeneous landscapes. Our study investigated the microclimate of south-exposed forests situated at the edge of sandstone cliffs in the western part of the North Alpine Foreland Basin in Switzerland and its role in local forest community composition. Using direct measurements from data loggers, as well as vegetation analyses, it was possible to quantify the microclimate of the cliff-edge forests and compare it with that of the surrounding forests. Our results highlighted the significant xerothermic and more variable nature of the cliff-edge forest microclimate, with a mean soil temperature up to 3.72°C warmer in the summer, higher annual (+28%) and daily (+250%) amplitudes of soil temperature, which frequently expose vegetation to extreme temperatures, and an 83% higher soil drying rate. These differences have a distinct influence on forest communities: cliff-edge forests are significantly different from surrounding forests. The site particularities of cliff edges support the presence of locally rare species and forest types, particularly of Scots pine. Cliff edges must therefore be considered microrefugia with a high conservation value for both xerothermic species and flora adapted to more continental climates. Moreover, the microclimate of cliff-edge forests could resemble the future climate in many ways. We argue that these small areas, which are already experiencing the future climate, can be seen as natural laboratories to better answer the following question: what will our forests look like in a few decades with accelerated climate change?

In Search of the Elusive North: Evolutionary History of the Arctic Fox (Vulpes lagopus) in the Palearctic from the Late Pleistocene to the Recent Inferred from Mitogenomic Data

Despite the high level of interest, the population history of arctic foxes during the Late Pleistocene and Holocene remains poorly understood. Here we aimed to fill gaps in the demographic and colonization history of the arctic fox by analyzing new ancient DNA data from fossil specimens aged from 50 to 1 thousand years from the Northern and Polar Urals, historic DNA from museum specimens from the Novaya Zemlya Archipelago and the Taymyr Peninsula and supplementing these data by previously published sequences of recent and extinct arctic foxes from other regions. This dataset was used for reconstruction of a time-calibrated phylogeny and a temporal haplotype network covering four time intervals: Late Pleistocene (ranging from 30 to 13 thousand years bp), Holocene (ranging from 4 to 1 thousand years bp), historical (approximately 150 years), and modern. Our results revealed that Late Pleistocene specimens showed no genetic similarity to either modern or historical specimens, thus supporting the earlier hypothesis on local extinction rather than habitat tracking.

Lake Superior's summer cooling of shorelines and adjacent inland forests: Implications for refugia of boreal forests and disjunct arctic-alpine plants

Climate refugia can serve as remnant habitat for cold-adapted species and delay forest transitions. The world's largest freshwater lake by surface area, Lake Superior, serves as a model system for understanding cooling-mediated refugia effects, as its cool summer water temperatures have maintained disjunct populations of arctic-alpine plants on its shoreline since deglaciation. It is known to affect local inland climates by providing a summer cooling effect; however, the inland temperature gradient and spatial patterns of cooling have not been well quantified. Here, we describe the extent, degree, and patterns of temperature buffering and examine drivers of buffering and disjunct plant occurrence for Lake Superior's north shore over a 3-year period at distances of 10, 100 m, 1, 10, and 100 km inland. We analyzed temperature data by year, month, summer maximum (July), and growing degree days (GDD0) for each site. Average summertime cooling at shore sites (10 m) was ~5°C cooler than reference sites (100 km inland), with a maximum difference of -19.2°C. The magnitude of cooling varied geographically, with sites further west and southeast showing little to no cooling effect, while the exposed north-central shore showed the highest degree of buffering (5.8°C cooler) and had a shorter growing season than reference sites. Finally, north-central shorelines had fewer days above 16°C, a threshold above which disjunct plants are unlikely to grow. These sites also showed the highest proportion of disjunct arctic-alpine species, reflecting the highest buffering from inland sites. On north-central shores, sites up to 10 km inland had less than 10 days per year warmer than 20°C, a threshold identified for boreal forest transition. An understanding of the extent of lake-mediated cooling on adjacent forests can better inform the risk to disjunct species, inland forests, and vegetation transition models on Lake Superior's north shore.

Is there a massive glacial-Holocene flora continuity in Central Europe?

The prevailing paradigm about the Quaternary ecological and evolutionary history of Central European ecosystems is that they were repeatedly impoverished by regional extinctions of most species during the glacial periods, followed by massive recolonizations from southern and eastern refugia during interglacial periods. Recent literature partially contradicts this view and provides evidence to re-evaluate this Postglacial Recolonization Hypothesis and develop an alternative one. We examined the long-term history of the flora of the Carpathian (Pannonian) Basin by synthesising recent advances in ecological, phylogeographical, palaeoecological and palaeoclimatological research, and analysing the cold tolerance of the native flora of a test area (Hungary, the central part of the Carpathian Basin). We found that (1) many species have likely occurred there continuously since before the Last Glacial Maximum (LGM); (2) most of the present-day native flora (1404 species, about 80%) can occur in climates as cold as or colder than the LGM (mean annual temperature ≤+3.5°C); and (3) grasslands and forests can be species-rich under an LGM-like cold climate. These arguments support an alternative hypothesis, which we call the Flora Continuity Hypothesis. It states that long-term continuity of much of the flora in the Carpathian Basin is more plausible than regional extinctions during the LGM followed by massive postglacial recolonizations. The long-term continuity of the region's flora may have fundamental implications not only for understanding local biogeography and ecology (e.g. the temporal scale of processes), but also for conservation strategies focusing on protecting ancient species-rich ecosystems and local gene pools.

Reinforcing the idea of an early dispersal of Hippopotamus amphibius in Europe: Restoration and multidisciplinary study of the skull from the Middle Pleistocene of Cava Montanari (Rome, central Italy)

A skull of Hippopotamus recovered from the area of Tor di Quinto, within the urban area of Rome (central Italy) is here redescribed. Despite being one of the most complete specimens of hippopotamuses of the European Pleistocene, the Tor di Quinto skull did not attract much research interest, due to long-standing uncertainties on its provenance. This work begun in 2021, when the skull was restored, within a large renovation project on the vertebrate exposed at the Earth Science University Museum of Sapienza University of Rome. Original sediments were found inside the cranial and mandible cavities during the restoration work, which were sampled for petrographic analyses. By combining a review of the old paleontological, archeological and geological literature published during the 19th and 20th century on the Rome basin and the correlation of these new sedimentological and petrographic information with the lithostratigraphic and synthemic units of the national geological cartography, we clarify that the Hippopotamus skull was most likely to have been collected from a quarry called Cava Montanari, from a formation dated between 560 and 460 ka. Morphological and biometric analyses clearly support an attribution of the Cava Montanari specimen to the extant species Hippopotamus amphibius. The reassessment of the stratigraphic and geological data on Cava Montanari implies that the studied specimen is the earliest confirmed occurrence of Hippopotamus amphibius in the European fossil record.

The molecular population structure of Swertia perennis (Gentianaceae) in Central Europe

Phylogeographic analysis of Swertia perennis, a typical European subalpine springtime species, revealed the existence of at least five major phylogenetic lineages. A large phylogeographic separation exists among these geographical regions, which confirms the existence of glacial refugia in the Pyrenees, but also in the Eastern and Central Alps. The results obtained from the analyzes indicate the existence of three major differences between the populations of the Alps and the Pyrenees, but also between the populations of the Alps and those of other geographical regions (Carpathians, southern Czech Republic, Sudetes and northern Poland). The studied populations from the Black Forest and from north-eastern and southern Poland are a relic of the former wider distribution of these (periglacial) genetic groups within Swertia perennis. Our results also confirm the existence of biogeographical links between the Carpathians and the Hercynian Range and the Alps. Certainly, there was an exchange of genes between populations located in the eastern Alps, the Carpathians and the Hercynian ranges (Czech Republic, Jeseníky, Sudetes, Ore Mountains). This confirms previous results of comparative studies on the genetic diversity of populations of other vascular plant species.

Deep-time paleogenomics and the limits of DNA survival

Although most ancient DNA studies have focused on the last 50,000 years, paleogenomic approaches can now reach into the early Pleistocene, an epoch of repeated environmental changes that shaped present-day biodiversity. Emerging deep-time genomic transects, including from DNA preserved in sediments, will enable inference of adaptive evolution, discovery of unrecognized species, and exploration of how glaciations, volcanism, and paleomagnetic reversals shaped demography and community composition. In this Review, we explore the state-of-the-art in paleogenomics and discuss key challenges, including technical limitations, evolutionary divergence and associated biases, and the need for more precise dating of remains and sediments. We conclude that with improvements in laboratory and computational methods, the emerging field of deep-time paleogenomics will expand the range of questions addressable using ancient DNA.

Hybridization of Atlantic puffins in the Arctic coincides with 20th-century climate change

The Arctic is experiencing the fastest rates of global warming, leading to shifts in the distribution of its biota and increasing the potential for hybridization. However, genomic evidence of recent hybridization events in the Arctic remains unexpectedly rare. Here, we use whole-genome sequencing of contemporary and 122-year-old historical specimens to investigate the origin of an Arctic hybrid population of Atlantic puffins (Fratercula arctica) on Bjørnøya, Norway. We show that the hybridization between the High Arctic, large-bodied subspecies F. a. naumanni and the temperate, smaller-sized subspecies F. a. arctica began as recently as six generations ago due to an unexpected southward range expansion of F. a. naumanni. Moreover, we find a significant temporal loss of genetic diversity across Arctic and temperate puffin populations. Our observations provide compelling genomic evidence of the impacts of recent distributional shifts and loss of diversity in Arctic communities during the 20th century.

Climate has contributed to population diversification of Daphnia galeata across Eurasia

Climate is a fundamental abiotic factor that plays a key role in driving the evolution, distribution and population diversification of species. However, there have been few investigations of genomic signatures of adaptation to local climatic conditions in cladocerans. Here, we have provided the first high-quality chromosome-level genome assembly (~143 Mb, scaffold N50 12.6 Mb) of the waterflea, Daphnia galeata, and investigated genomic variation in 22 populations from Central Europe and Eastern China. Our ecological-niche models suggested that the historic distribution of D. galeata in Eurasia was significantly affected by Quaternary climate fluctuations. We detected pronounced genomic and morphometric divergences between European and Chinese D. galeata populations. Such divergences could be partly explained by genomic signatures of thermal adaptation to distinct climate regimes: a set of candidate single-nucleotide polymorphisms (SNPs) potentially associated with climate were detected. These SNPs were in genes significantly enriched in the Gene ontology terms "determination of adult lifespan" and "translation repressor activity", and especially, mthl5 and SOD1 involved in the IIS pathway, and EIF4EBP2 involved in the target of the rapamycin signalling pathway. Our study indicates that certain alleles might be associated with particular temperature regimes, playing a functional role in shaping the population structure of D. galeata at a large geographical scale. These results highlight the potential role of molecular variation in the response to climate variation, in the context of global climate change.

Ice holes microrefugia harbor genetically and functionally distinct populations of Vaccinium vitis-idaea (Ericaceae)

In the mountain terrain, ice holes are little depressions between rock boulders that are characterized by the exit of cold air able to cool down the rock surface even in summer. This cold air creates cold microrefugia in warmer surroundings that preserve plant species probably over thousands of years under extra-zonal climatic conditions. We hypothesized that ice hole populations of the model species Vaccinium vitis-idaea (Ericaceae) show genetic differentiation from nearby zonal subalpine populations, and high functional trait distinctiveness, in agreement with genetic patterns. We genotyped almost 30,000 single nucleotide polymorphisms using restriction site-associated DNA sequencing and measured eight functional traits indicative of individual performance and ecological strategies. Genetic results showed high differentiation among the six populations suggesting isolation. On siliceous bedrock, ice hole individuals exhibited higher levels of admixture than those from subalpine populations which could have experienced more bottlenecks during demographic fluctuations related to glacial cycles. Ice hole and subalpine calcareous populations clearly separated from siliceous populations, indicating a possible effect of bedrock in shaping genetic patterns. Trait analysis reflected the bedrock effect on populations' differentiation. The significant correlation between trait and genetic distances suggests the genetic contribution in shaping intraspecific functional differentiation. In conclusion, extra-zonal populations reveal a prominent genetic and phenotypic differentiation determined by history and ecological contingency. Therefore, microrefugia populations can contribute to the overall variability of the species and lead to intraspecific-driven responses to upcoming environmental changes.

Beringia as a high-latitude engine of avian speciation

Beringia is a biogeographically dynamic region that extends from northeastern Asia into northwestern North America. This region has affected avian divergence and speciation in three important ways: (i) by serving as a route for intercontinental colonisation between Asia and the Americas; (ii) by cyclically splitting (and often reuniting) populations, subspecies, and species between these continents; and (iii) by providing isolated refugia through glacial cycles. The effects of these processes can be seen in taxonomic splits of shallow to increasing depths and in the presence of regional endemics. We review the taxa involved in the latter two processes (splitting-reuniting and isolation), with a focus on three research topics: avian diversity, time estimates of the generation of that diversity, and the regions within Beringia that might have been especially important. We find that these processes have generated substantial amounts of avian diversity, including 49 pairs of avian subspecies or species whose breeding distributions largely replace one another across the divide between the Old World and the New World in Beringia, and 103 avian species and subspecies endemic to this region. Among endemics, about one in three is recognised as a full biological species. Endemic taxa in the orders Charadriiformes (shorebirds, alcids, gulls, and terns) and Passeriformes (perching birds) are particularly well represented, although they show very different levels of diversity through evolutionary time. Endemic Beringian Charadriiformes have a 1.31:1 ratio of species to subspecies. In Passeriformes, endemic taxa have a 0.09:1 species-to-subspecies ratio, suggesting that passerine (and thus terrestrial) endemism might be more prone to long-term extinction in this region, although such 'losses' could occur through their being reconnected with wider continental populations during favourable climatic cycles (e.g. subspecies reintegration with other populations). Genetic evidence suggests that most Beringian avian taxa originated over the past 3 million years, confirming the importance of Quaternary processes. There seems to be no obvious clustering in their formation through time, although there might be temporal gaps with lower rates of diversity generation. For at least 62 species, taxonomically undifferentiated populations occupy this region, providing ample potential for future evolutionary diversification.

Effectiveness of phase synchronization in chaotic food chain model with refugia and Allee effects during seasonal fluctuations

Seasonal effects powerfully shape the population dynamics with periodic climate changes because species naturally adjust their dynamics with seasonal variations. In response to these effects, sometimes population dynamics exhibit synchrony or generate chaos. However, synchronized dynamics enhance species' persistence in naturally unstable environments; thus, it is imperative to identify parameters that alter the dynamics of an ecosystem and bring it into synchrony. This study examines how ecological parameters enable species to adapt their dynamics to seasonal changes and achieve phase synchrony within ecosystems. For this, we incorporate seasonal effects as a periodic sinusoidal function into a tri-trophic food chain system where two crucial bio-controlling parameters, Allee and refugia effects, are already present. First, it is shown that the seasonal effects disrupt the limit cycle and bring chaos to the system. Further, we perform rigorous mathematical analysis to perform the dynamical and analytical properties of the nonautonomous version of the system. These properties include sensitive dependence on initial condition (SDIC), sensitivity analysis, bifurcation results, the positivity and boundedness of the solution, permanence, ultimate boundedness, and extinction scenarios of species. The SDIC characterizes the presence of chaotic oscillations in the system. Sensitivity analysis determines the parameters that significantly affect the outcome of numerical simulations. The bifurcation study concerning seasonal parameters shows a higher dependency of species on the frequency of seasonal changes than the severity of the season. The bifurcation study also examines the bio-controlling parameters and reveals various dynamic states within the system, such as fold, transcritical branch points, and Hopf points. Moreover, the mathematical analysis of our seasonally perturbed system reveals the periodic coexistence of all species and a globally attractive solution under certain parametric constraints. Finally, we examine the role of essential parameters that contribute to phase synchrony. For this, we numerically investigate the defining role of the coupling dimension coefficient, bio-controlling parameters, and other parameters associated with seasonality. This study infers that species can tune their dynamics to seasonal effects with low seasonal frequency, whereas the species' tolerance for the severity of seasonal effects is relatively high. The research also sheds light on the correlation between the degree of phase synchrony, prey biomass levels, and the severity of seasonal forcing. This study offers valuable insights into the dynamics of ecosystems affected by seasonal perturbations, with implications for conservation and management strategies.

Tropical coral reefs in Sri Lanka are threatened due to the fluctuation of seasonal and interannual sea surface temperature

The thermal stress of oceans causes coral bleaching, which induces loss of life in coral reefs and makes them exposed to other threats which directly and indirectly affect millions of other species that inhabit the reef. However, studies focusing on how those thermal stresses affect Sri Lankan fringing reef ecosystems are scarce. Hence, the patterns of long-term and short-term fluctuations of sea surface temperature (SST) over the shallow reefs around the country were studied by separating them into different zones as the eastern coast (Passikudha, Kayankerni, Adukkuparu, Parrot Rock, and Pigeon Island), the southern coast (Beruwala Barbarian, Hikkaduwa, Unawatuna, Ahangama, Mirissa, Madiha, Polhena, and Devundara), and northern-northwestern coasts (Valiththoondal, Palk Bay, Mannar, Kalpitiya, Thalwila, and Uswatakeiyawa). The 1 km Multiscale Ultrahigh Resolution (MUR) Level 4 SST dataset was used to analyze seasonal and interannual SST variability from 2005 to 2021. The data were correlated with the Indian Ocean Dipole (IOD), Ekman velocity, and wind stress curl. The annual, seasonal, and monthly variability of SST on different coasts is significantly different. Higher increasing trends of SST from 0.0324 to 0.0411 ℃/year are observed on different coasts, and most of the time higher positive anomalies are recorded after 2014. The First Inter Monsoon (IM-1) and the month of April are more critical since they reach the maximum SST, and the minimum is in the North West Monsoon (NWM) and January. Significant positive relationships are recorded between the Indian Ocean Dipole (IOD) index and the monthly average SST on different coasts, which was robust on the southern coast. Therefore, tropical coral reefs in Sri Lanka are severely threatened due elevation of SST by global warming and climate variability.

Environmental heterogeneity increases the conservation value of small natural features in karst landscapes

Local biodiversity hotspots are often located within regions where extreme and variable environmental - e.g., climatic and soil - conditions occur. These areas are conservation priorities. Although environmental heterogeneity is recognised as an important determinant of biodiversity, studies focusing on the effects of multiple environmental heterogeneity components in the same ecosystem are scarce. Here we investigate how topography and related microclimatic variables and soil properties may influence the biodiversity and conservation value of karst landscapes. Karst landscapes of the world contain millions of dolines (i.e. bowl- or funnel-shaped depressions) that may function as 'small natural features' with a disproportionately large role in maintaining biodiversity relative to their size. We assessed the diversity of microclimates, soils and vegetation and their relationships in six microhabitats (south-facing slopes, east-facing slopes, west-facing slopes, north-facing slopes and bottoms of dolines, and the adjacent plateau) for nine large dolines in a grassland ecosystem. Although there were remarkable differences among the conservation value of these microhabitats (e.g., representation of different species groups, presence of 'climate relicts'), each microhabitat had an important role in maintaining species that are rare or absent in other microhabitats in the landscape. We found that the studied dolines exhibited highly variable environmental conditions and promoted a high diversity of vegetation types with unique species composition, contributing to the topographic, climatic, soil, vegetation and land cover heterogeneity of karst landscapes. Therefore, our findings highlight that dolines may function as local biodiversity hotspots and have a crucial conservation importance. As dolines are widespread topographic features in many karst landscapes throughout the world, our results could be directly applied to other regions as well. An integrated approach is urgently needed to provide guidelines for landscape management, promoting the retention of the microhabitat diversity of small natural features for species vulnerable to climate change and/or various disturbances.

Unrecognized diversity, genetic structuring, and phylogeography of the genus Triplophysa (Cypriniformes: Nemacheilidae) sheds light on two opposite colonization routes during Quaternary glaciation that occurred in the Qilian Mountains

In recent years, the influence of historical geological and climatic events on the evolution of flora and fauna in the Tibetan Plateau has been a hot research topic. The Qilian Mountain region is one of the most important sources of biodiversity on the Qinghai-Tibet Plateau. Many species existed in the region during the Pleistocene glacial oscillation, and the complex geographical environment provided suitable conditions for the survival of local species. The shrinkage, expansion, and transfer of the distribution range and population size of species have significant effects on genetic diversity and intraspecific differentiation. To reveal the effects of geological uplift and climate oscillation on the evolution of fish populations in the Qilian Mountains, we investigated the genetic structure, phylogenetic relationship, and phylogeographical characteristics of genus Triplophysa species in the Qilian Mountains using the mitochondrial DNA gene (COI), three nuclear genes (RAG1, sRH, and Myh6) and 11 pairs of nuclear microsatellite markers. We collected 11 species of genus Triplophysa living in the Qilian Mountains, among which Triplophysa hsutschouensis and Triplophysa papillosolabiata are widely distributed in the rivers on the northern slope of the Qilian Mountains. There was a high degree of lineage differentiation among species, and the genetic diversity of endemic species was low. The different geographical groups of T. papillosolabiata presented some allogeneic adaptation and differentiation, which was closely related to the changes in the river system. Except for the population expansion event of T. hsutschouensis during the last glacial period of the Qinghai-Tibet Plateau (0.025 MYA), the population sizes of other plateau loach species remained stable without significant population expansion. Starting from the east and west sides of the Qilian Mountains, T. hsutschouensis, and T. papillosolabiata showed two species colonization routes in opposite directions. The geological events of the uplift of the Qinghai-Tibet Plateau and the climatic oscillation of the Quaternary glaciation had a great influence on the genetic structure of the plateau loach in the Qilian Mountains, which promoted the genetic differentiation of the plateau loach and formed some unique new species. The results of this study have important guiding significance for fish habitat protection in the Qilian Mountains.

Phylogenetic divergence associated with climate oscillations and topology illustrates the dispersal history of Formosan sambar deer (Rusa unicolor swinhoii) in Taiwan

The island of Taiwan represents an ideal context for studying the effects of climatic oscillations and topographic variation on large herbivores due to its varied tropical to sub-tropical climate zones at different elevational ranges. We explored the phylogenetics of Formosan sambar deer (Rusa unicolor swinhoii) using the control region of the mitochondrial genome. We detected 18 haplotypes among 454 sequences across the island and grouped them into six regions based on SAMOVA, with 68.78% variance among regions. A Bayesian phylogenetic dendrogram revealed two spatially segregated genetic clades. Neutrality tests and Bayesian skyline plots uncovered different demographic expansion histories for the two clades. We further tested divergence times and chronology to propose potential phylogenetic scenarios, which were examined using approximate Bayesian computation. Finally, we present a credible hypothesis for a glacial refugium in the northern part of the Central Mountain Range. Subsequent secondary contact between the two clades during interglacial periods has led to the extant genetic structure of Formosan sambar deer.

Geographic Genetic Structure of Alectoris chukar in Türkiye: Post-LGM-Induced Hybridization and Human-Mediated Contaminations

Türkiye is considered an important evolutionary area for Chukar partridge (Alectoris chukar), since it is both a potential ancestral area and a diversification center for the species. Using 2 mitochondrial (Cty-b and D-loop) and 13 polymorphic microsatellite markers, we investigated the geographic genetic structure of A. chukar populations to determine how past climatic fluctuations and human activities have shaped the gene pool of this species in Türkiye. Our results indicate, firstly, that only A. chukar of the genus Alectoris is present in Türkiye (Anatolia and Thrace), with no natural or artificial gene flow from congenerics. Secondly, the geographic genetic structure of the species in Türkiye has been shaped by topographic heterogeneity, Pleistocene climatic fluctuations, and artificial transport by humans. Third, there appears to be three genetic clusters: Thracian, Eastern, and Western. Fourth, the post-LGM demographic expansion of the Eastern and Western populations has formed a hybrid zone in Central Anatolia (~8 kyBP). Fifth, the rate of China clade-B contamination in Türkiye is about 8% in mtDNA and about 12% in nuDNA, with the Southeastern Anatolian population having the highest contamination. Sixth, the Thracian population was the most genetically distinct, with the lowest genetic diversity and highest level of inbreeding and no China clad-B contamination. These results can contribute to the conservation regarding A. chukar populations, especially the Thracian population.

Phylogeography within the Peromyscus maniculatus species group: Understanding past distribution of genetic diversity and areas of refugia in western North America

The effects of anthropogenic climate change on biodiversity have been recognized on every continent, ocean, and across different taxonomic groups. Here, we study the range dynamics and demography of a cosmopolitan species: the deer mouse, Peromyscus maniculatus. We generated a multilocus SNP dataset using the ddRADseq protocol for 218 individuals across the geographic range within three western North American lineages of this species group. We evaluated population structure using several methods and explored the correlation between geographic and genetic distances. We modeled the demographic history using a site frequency spectrum approach and used a machine learning algorithm to infer current and past (Last Glacial Maximum; LGM) environmental suitability. Lastly, we explored the origin of population expansion for the identified lineages. The genome-wide SNP dataset was able to identify-three regionally distinct groups- 1) P. m. gambelii (southern California); 2) P. keeni (Pacific Northwest); 3) P. m. sonoriensis (a broad population spanning the Pacific Northwest through central California and across the Rocky Mountains into the Great Plains). Demographic analysis indicated the splits between the three populations occurred within the last 500 thousand years, with one very recent (late Holocene) split. Ecological niche models for each of these lineages predicted suitable environment present throughout their known ranges for current conditions, and a severe reduction of northern habitat in the past. The deer mouse has responded to past climate changes by expanding its range during interglacial periods and contracting its range during glacial periods leading to strong population differentiation. But lower magnitude climate change or other processes within the Holocene interglacial period led to population differentiation as well, which is likely still ongoing today given the substantial anthropogenic climate change and other landscape transformations caused by humans during the Anthropocene. By understanding the historical processes that led to the contemporary geographic distribution of biodiversity, we can determine the relative importance of different factors that shape biodiversity, now and into the future.

Ancient DNA of narrow-headed vole reveal common features of the Late Pleistocene population dynamics in cold-adapted small mammals

The narrow-headed vole, collared lemming and common vole were the most abundant small mammal species across the Eurasian Late Pleistocene steppe-tundra environment. Previous ancient DNA studies of the collared lemming and common vole have revealed dynamic population histories shaped by climatic fluctuations. To investigate the extent to which species with similar adaptations share common evolutionary histories, we generated a dataset comprised the mitochondrial genomes of 139 ancient and 6 modern narrow-headed voles from several sites across Europe and northwestern Asia covering approximately the last 100 thousand years (kyr). We inferred Bayesian time-aware phylogenies using 11 radiocarbon-dated samples to calibrate the molecular clock. Divergence of the main mtDNA lineages across the three species occurred during marine isotope stages (MIS) 7 and MIS 5, suggesting a common response of species adapted to open habitat during interglacials. We identified several time-structured mtDNA lineages in European narrow-headed vole, suggesting lineage turnover. The timing of some of these turnovers was synchronous across the three species, allowing us to identify the main drivers of the Late Pleistocene dynamics of steppe- and cold-adapted species.

A resilient and connected network of sites to sustain biodiversity under a changing climate

Motivated by declines in biodiversity exacerbated by climate change, we identified a network of conservation sites designed to provide resilient habitat for species, while supporting dynamic shifts in ranges and changes in ecosystem composition. Our 12-y study involved 289 scientists in 14 study regions across the conterminous United States (CONUS), and our intent was to support local-, regional-, and national-scale conservation decisions. To ensure that the network represented all species and ecosystems, we stratified CONUS into 68 ecoregions, and, within each, we comprehensively mapped the geophysical settings associated with current ecosystem and species distributions. To identify sites most resilient to climate change, we identified the portion of each geophysical setting with the most topoclimate variability (high landscape diversity) likely to be accessible to dispersers (high local connectedness). These "resilient sites" were overlaid with conservation priority maps from 104 independent assessments to indicate current value in supporting recognized biodiversity. To identify key connectivity areas for sustaining species movement in response to climate change, we codeveloped a fine-scale representation of human modification and ran a circuit-theory-based analysis that emphasized movement potential along geographic climate gradients. Integrating areas with high values for two or more factors, we identified a representative, resilient, and connected network of biodiverse lands covering 35% of CONUS. Because the network connects climatic gradients across 250,000 biodiversity elements and multiple resilient examples of all geophysical settings in every ecoregion, it could form the spatial foundation for targeted land protection and other conservation strategies to sustain a diverse, dynamic, and adaptive world.

Identifying the Past, Present, and Future Distribution Patterns of the Balkan Wall Lizard (Sauria: Lacertidae: Podarcis tauricus) by Ecological Niche Modelling

Collapse

Evolutionary Potential of Parthenogenesis-Bisexual Lineages within Triploid Apomictic Thelytoky in Cacopsylla ledi (Flor, 1861) (Hemiptera, Psylloidea) in Fennoscandia

A widely accepted hypothesis is that parthenogenesis is an evolutionary dead end since it is selectively advantageous in the short term only but results in lowered diversification rates. Triploid apomictic parthenogenesis might represent an exception, as in favorable environments, triploid females are able to produce rare males and diploid females. The aim of the present study was to analyze the modes of reproduction and their evolutionary implications in the parthenogenetic psyllid Cacopsylla ledi (Flor, 1861) from Fennoscandia. The cytogenetic assessment of ploidy levels and the analysis of the COI haplotype revealed two geographically separated bisexual lineages implying genuine bisexual populations. The southern lineage occurring south of latitude 65° N in Finland showed a COI haplotype different from that of parthenogenetic triploids in the same population but identical to the haplotype of specimens in a genuine bisexual population in the Czech Republic. This allows us to suggest that bisexuals in southern Fennoscandia represent the original bisexual C. ledi. By contrast, in the northern bisexual lineage north of latitude 65° N, rare males and diploid females carried the same haplotype as triploids in the same population, having been produced by the triploids. In the Kola Peninsula, a genuine bisexual population of presumably rare male/diploid female origin was discovered. As this population is geographically isolated from populations of the ancestral bisexual C. ledi, it can develop into a new bisexual species through peripatric speciation during evolution. Our findings demonstrate that apomictic triploid parthenogenesis is not necessarily an evolutionary dead end but is able to lead to the emergence of a new bisexual species of parthenogenetic origin.

From DNA barcodes to ecology: Meta-analysis of central European beetles reveal link with species ecology but also to data pattern and gaps

DNA barcoding has been used worldwide to identify biological specimens and to delimit species. It represents a cost-effective, fast, and efficient way to assess biodiversity with help of the public Barcode of Life Database (BOLD) accounting for more than 236,000 animal species and more than 10 million barcode sequences. Here, we performed a meta-analysis of available barcode data of central European Coleoptera to detect intraspecific genetic patterns among ecological groups in relation to geographic distance with the aim to investigate a possible link between infraspecific variation and species ecology. We collected information regarding feeding style, body size, as well as habitat and biotope preferences. Mantel tests and two variants of Procrustes analysis, both involving the Principal Coordinates Neighborhood Matrices (PCNM) approach, were applied on genetic and geographic distance matrices. However, significance levels were too low to further use the outcome for further trait investigation: these were in mean for all ecological guilds only 7.5, 9.4, or 15.6% for PCNM + PCA, NMDS + PCA, and Mantel test, respectively, or at best 28% for a single guild. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. Results suggest that increased numbers of species, sampling localities, and specimens for a chosen area of interest may give new insights to explore barcode data and species ecology for the scope of conservation on a larger scale. We performed a meta-analysis of available barcode data of central European beetles to detect intraspecific genetic patterns among ecological groups in relation to geographic distance, regarding feeding style, body size, as well as habitat and biotope preferences. Our study confirmed that certain ecological traits were associated with higher species diversity and foster stronger genetic differentiation. However, significance levels were too low to further use the outcome for further trait investigation.

Population dynamics and demographic history of Eurasian collared lemmings

BACKGROUND

Ancient DNA studies suggest that Late Pleistocene climatic changes had a significant effect on population dynamics in Arctic species. The Eurasian collared lemming (Dicrostonyx torquatus) is a keystone species in the Arctic ecosystem. Earlier studies have indicated that past climatic fluctuations were important drivers of past population dynamics in this species.

RESULTS

Here, we analysed 59 ancient and 54 modern mitogenomes from across Eurasia, along with one modern nuclear genome. Our results suggest population growth and genetic diversification during the early Late Pleistocene, implying that collared lemmings may have experienced a genetic bottleneck during the warm Eemian interglacial. Furthermore, we find multiple temporally structured mitogenome clades during the Late Pleistocene, consistent with earlier results suggesting a dynamic late glacial population history. Finally, we identify a population in northeastern Siberia that maintained genetic diversity and a constant population size at the end of the Pleistocene, suggesting suitable conditions for collared lemmings in this region during the increasing temperatures associated with the onset of the Holocene.

CONCLUSIONS

This study highlights an influence of past warming, in particular the Eemian interglacial, on the evolutionary history of the collared lemming, along with spatiotemporal population structuring throughout the Late Pleistocene.

Barriers to gene flow in the deepest ocean ecosystems: Evidence from global population genomics of a cosmopolitan amphipod

The deepest marine ecosystem, the hadal zone, hosts endemic biodiversity resulting from geographic isolation and environmental selection pressures. However, the pan-ocean distribution of some fauna challenges the concept that the hadal zone is a series of isolated island-like habitats. Whether this remains true at the population genomic level is untested. We investigated phylogeographic patterns of the amphipod, Bathycallisoma schellenbergi, from 12 hadal features across the Pacific, Atlantic, Indian, and Southern oceans and analyzed genome-wide single-nucleotide polymorphism markers and two mitochondrial regions. Despite a cosmopolitan distribution, populations were highly restricted to individual features with only limited gene flow between topographically connected features. This lack of connectivity suggests that populations are on separate evolutionary trajectories, with evidence of potential cryptic speciation at the Atacama Trench. Together, this global study demonstrates that the shallower ocean floor separating hadal features poses strong barriers to dispersal, driving genetic isolation and creating pockets of diversity to conserve.

Phylogeography reveals the origin of the two phenological forms of large blue, Phengaris arion (Lepidoptera: Lycaenidae)

Our main goal was to investigate the phylogeography of the butterfly Phengaris arion to reveal the evolutionary origin of its ‘spring’ and ‘summer’ forms. Molecular analyses based on highly variable microsatellites, together with Wolbachia screening, were carried out on 34 populations in Europe. We found three well-defined genetic lineages of different origins: the Apennine, the central and the eastern. The highly distinct Apennine lineage is limited by the Alps and evaluated as an Evolutionary Significant Unit (ESU). Therefore, the taxon name ligurica, described from the Ligurian coast (Italy), should not be applied to denote the ‘summer form’ of the Pannonian region. The central lineage is limited by the Carpathians and the most eastern ranges of the Alps, and lacks major range fluctuations related to glaciations, although there is evidence for extra-Mediterranean refugia in the Carpathian Basin. The eastern clade could have had refugia in central Asia. Our results exclude the potential allopatric origin of the ‘spring’ and ‘summer’ arion, and support the hypothesis that the existence of the two forms could be a result of local adaptation to the distinctive phenology of host plant flowering which is manifested in the genetic differences between them. Wolbachia infection has been ruled out as a driver of sympatric speciation in P. arion.

Phylogeography and population genetics of a widespread cold-adapted ant, Prenolepis imparis

As species arise, evolve, and diverge, they are shaped by forces that unfold across short and long time scales and at both local and vast geographic scales. It is rare, however, to be able document this history across broad sweeps of time and space in a single species. Here, we report the results of a continental-scale phylogenomic analysis across the entire range of a widespread species. We analyzed sequences of 1,402 orthologous Ultraconserved Element (UCE) loci from 75 individuals to identify population genetic structure and historical demographic patterns across the continent-wide range of a cold-adapted ant, the winter ant, Prenolepis imparis. We recovered five well-supported, genetically isolated clades representing lineages that diverged from 8.2-2.2 million years ago. These include: 1) an early diverging lineage located in Florida, 2) a lineage that spans the southern United States, 3) populations that extend across the midwestern and northeastern United States, 4) populations from the western United States, and 5) populations in southwestern Arizona and Mexico. Population genetic analyses revealed little or no gene flow among these lineages, but patterns consistent with more recent gene flow among populations within lineages, and localized structure with migration in the western United States. High support for five major geographic lineages and lack of evidence of contemporary gene flow indicate in situ diversification across the species' range, producing relatively ancient lineages that persisted through subsequent climate change and glaciation during the Quaternary.

Linking large-scale genetic structure of three Argynnini butterfly species to geography and environment

Understanding which factors and processes are associated with genetic differentiation within and among species remains a major goal in evolutionary biology. To explore differences and similarities in genetic structure and its association with geographical and climatic factors in sympatric sister species, we conducted a large‐scale (>32° latitude and >36° longitude) comparative phylogeographical study on three Argynnini butterfly species (Speyeria aglaja, Fabriciana adippe and F. niobe) that have similar life histories, but differ in ecological generalism and dispersal abilities. Analyses of nuclear (ddRAD‐sequencing derived SNP markers) and mitochondrial (COI sequences) data revealed differences between species in genetic structure and how genetic differentiation was associated with climatic factors (temperature, solar radiation, precipitation, wind speed). Geographical proximity accounted for much of the variation in nuclear and mitochondrial structure and evolutionary relationships in F. adippe and F. niobe, but only explained the pattern observed in the nuclear data in S. aglaja, for which mitonuclear discordance was documented. In all species, Iberian and Balkan individuals formed genetic clusters, suggesting isolation in glacial refugia and limited postglacial expansion. Solar radiation and precipitation were associated with the genetic structure on a regional scale in all species, but the specific combinations of environmental and geographical factors linked to variation within species were unique, pointing to species‐specific responses to common environments. Our findings show that the species share similar colonization histories, and that the same ecological factors, such as niche breadth and dispersal capacity, covary with genetic differentiation within these species to some extent, thereby highlighting the importance of comparative phylogeographical studies in sympatric sister species.

Late quaternary biotic homogenization of North American mammalian faunas

Biotic homogenization—increasing similarity of species composition among ecological communities—has been linked to anthropogenic processes operating over the last century. Fossil evidence, however, suggests that humans have had impacts on ecosystems for millennia. We quantify biotic homogenization of North American mammalian assemblages during the late Pleistocene through Holocene (~30,000 ybp to recent), a timespan encompassing increased evidence of humans on the landscape (~20,000–14,000 ybp). From ~10,000 ybp to recent, assemblages became significantly more homogenous (>100% increase in Jaccard similarity), a pattern that cannot be explained by changes in fossil record sampling. Homogenization was most pronounced among mammals larger than 1 kg and occurred in two phases. The first followed the megafaunal extinction at ~10,000 ybp. The second, more rapid phase began during human population growth and early agricultural intensification (~2,000–1,000 ybp). We show that North American ecosystems were homogenizing for millennia, extending human impacts back ~10,000 years.

Biotic homogenization, which is increased similarity in the composition of species among communities, is rising due to human activities. Using North American mammal fossil records from the past 30,000 years, this study shows that this phenomenon is ancient, beginning between 12,000 and 10,000 years ago with the extinction of the mammal megafauna.

Phylogeography and Wolbachia Infections Reveal Postglacial Recolonization Routes of the Parthenogenetic Plant Louse Cacopsylla myrtilli (W. Wagner 1947), (Hemiptera, Psylloidea)

To reveal the phylogeographic pattern of the parthenogenetic psyllid Cacopsylla myrtilli (W. Wagner 1947) (Hemiptera, Psylloidea), we sequenced a 638 bp fragment of the mitochondrial COI gene from 962 individuals. These insects originated from 46 sampling sites, which cover a significant part of the northern Palearctic distribution range of the species. The sequence data revealed 40 haplotypes, with three main (H1, H2, and H3) and 37 derived ones. The main haplotypes H1 or H2 or both were present at all sampling sites. The star-like shape of the haplotype networks indicated recent population expansion. In most cases, the derived haplotypes were specific for each country, suggesting that the main haplotypes H1 and H2 are of refugial origin, and the derived haplotypes have emerged after the postglacial recolonization process. Based on the haplotype sequences, we suggest H3 to represent the ancestral haplotype from which H1 and H2 have evolved. We suggest that the main haplotype H3 together with its derived haplotypes represents bisexual C. myrtilli, which shows a limited distribution on both sides of the border between Finland and Russia in northern Fennoscandia. The genetic diversity was the highest in Sjoa in southern Norway and also high in the White Sea region in northwest Russia. Higher diversity in Sjoa was attributed to both earlier recolonizations compared to that of the White Sea region and the absence of Wolbachia infection. We suggest that these sites were colonized from different Pleistocene refugia, i.e., from western and eastern refugia, respectively. From the White Sea region, recolonization continued eastwards to Ural Mountains and westwards to Finland and further north to Kola Peninsula. From northern Finland, recolonization continued to Finnmark, Norway, and further to Sweden and finally reached a secondary contact zone with colonizers from Norway in Central Sweden. The Caucasus and Siberian/Manchurian refugial regions have played an important role in the origin of C. myrtilli populations in Siberia and the Russian Far East.

Late Holocene anthropogenic landscape change in northwestern Europe impacted insect biodiversity as much as climate change did after the last Ice Age

Since the last Ice Age (ca 115 000-11 700 years ago), the geographical ranges of most plants and animals have shifted, expanded or contracted. Understanding the timing, geographical patterns and drivers of past changes in insect communities is essential for evaluating the biodiversity implications of future climate changes, yet our knowledge of long-term patterns is limited. We applied a network modelling approach to the recent fossil record of northwestern European beetles to investigate how their taxonomic and trait composition changed during the past 16 000 years. We found two major changes in beetle faunas 4000-3500 and 10 000-9500 years ago, coinciding with periods of human population growth in the Late Holocene and climate warming in the Early Holocene. Our results demonstrate that humans have affected insect biodiversity since at least the introduction of agropastoralism, with landscape-scale effects that can be observed at sites away from areas of direct human impact.

Population structure in Quercus suber L. revealed by nuclear microsatellite markers

Quercus suber L. is a sclerophyllous tree species native to the western Mediterranean, a region that is considered highly vulnerable to increased temperatures and severe dry conditions due to environmental changes. Understanding the population structure and demographics of Q. suber is essential in order to anticipate whether populations at greater risk and the species as a whole have the genetic background and reproductive dynamics to enable rapid adaptation. The genetic diversity of Q. suber has been subject to different studies using both chloroplast and nuclear data, but population structure patterns remain unclear. Here, we perform genetic analyses on Q. suber using 13 nuclear microsatellite markers, and analysed 17 distinct locations across the entire range of the species. Structure analyses revealed that Q. suber may contain three major genetic clusters that likely result from isolation in refugia combined with posterior admixture and putative introgression from other Quercus species. Our results show a more complex structure scenario than previously inferred for Q. suber using nuclear markers and suggest that different southern populations contain high levels of genetic variation that may contribute to the resilience of Q. suber in a context of environmental change and adaptive pressure.

Variations by degrees: Western European paleoenvironmental fluctuations across MIS 13-11

Marine Isotope Stages (MIS) 13-11 saw a major transformation in the hominin occupation of Europe, with an expansion in the scale and geographical distribution of sites and artifact assemblages. That expansion is explored here in the context of paleoenvironmental variability, focusing on geographical and chronological trends in climatic and habitat conditions at and between key Lower Paleolithic sites in Western Europe. Climatic conditions at British sites are compared across MIS 13-11, and used to test predicted values from the Oscillayers data set. Conditions at hominin and nonhominin sites are compared to explore possible limitations in hominin tolerances during MIS 13-11. Trends in conditions are explored with reference to long-term global patterns, short-term substage events, and seasonal variations. The apparent increase in the scale of hominin activity in north-western Europe during MIS 13 is surprising in light of the relatively harsh conditions of late MIS 13, and is likely to reflect significant physiological and/or behavioral adaptations, a mild south-north temperature gradient in western Europe during MIS 13, and the relatively mild, sustained conditions spanning MIS 15-13. The expanded occupation of north-western Europe during MIS 11 probably reflects the extended mild conditions of MIS 11c, since marked seasonal temperature differences and substantial behavioral changes between hominin sites in MIS 13 and 11 are not clearly evident. Site-specific conditions in south-western Europe during MIS 11 suggest milder winters, warmer summers, and reduced seasonal variability compared to north-western Europe. Some or all of these conditions may have supported larger, core populations, as may the relatively mild conditions associated with south-western European sites during MIS 12. Finally, comparisons between north-western and north-central European sites indicate relatively small differences in seasonal temperatures, suggesting that climate may only be a partial factor behind the smaller-scale occupations of north-central Europe during MIS 13-11.