The Last Glacial Maximum

Range-Wide Phylogeography and Ecological Niche Modeling Provide Insights into the Evolutionary History of the Mongolian Racerunner (Eremias argus) in Northeast Asia

The Mongolian racerunner, Eremias argus, is a small lizard endemic to Northeast Asia that can serve as an excellent model for investigating how geography and past climate change have jointly influenced the evolution of biodiversity in this region. To elucidate the processes underlying its diversification and demography, we reconstructed the range-wide phylogeographic pattern and evolutionary trajectory, using phylogenetic, population genetic, landscape genetic, Bayesian phylogeographic reconstruction and ecological niche modeling approaches. Phylogenetic analyses of the mtDNA cyt b gene revealed eight lineages that were unbounded by geographic region. The genetic structure of E. argus was mainly determined by geographic distance. Divergence dating indicated that E. argus and E. brenchleyi diverged during the Mid-Pliocene Warm Period. E. argus was estimated to have coalesced at~0.4351 Ma (Marine Isotope Stage 19). Bayesian phylogeographic diffusion analysis revealed out-of-Inner Mongolia and rapid colonization events from the end of the Last Interglacial to the Last Glacial Maximum, which is consistent with the expanded suitable range of the Last Glacial Maximum. Pre-Last Glacial Maximum growth of population is presented for most lineages of E. argus. The Glacial Maximum contraction model and the previous multiple glacial refugia hypotheses are rejected. This may be due to an increase in the amount of climatically favorable habitats in Northeast Asia. Furthermore, E. argus barbouri most likely represents an invalid taxon. The present study is the first to report a range-wide phylogeography of reptiles over such a large region in Northeast Asia. Our results make a significant contribution towards understanding the biogeography of the entire Northeast Asia.

A Late Pleistocene coastal ecosystem in French Guiana was hyperdiverse relative to today

Warmer temperatures and higher sea level than today characterized the Last Interglacial interval [Pleistocene, 128 to 116 thousand years ago (ka)]. This period is a remarkable deep-time analog for temperature and sea-level conditions as projected for 2100 AD, yet there has been no evidence of fossil assemblages in the equatorial Atlantic. Here, we report foraminifer, metazoan (mollusks, bony fish, bryozoans, decapods, and sharks among others), and plant communities of coastal tropical marine and mangrove affinities, dating precisely from a ca. 130 to 115 ka time interval near the Equator, at Kourou, in French Guiana. These communities include ca. 230 recent species, some being endangered today and/or first recorded as fossils. The hyperdiverse Kourou mollusk assemblage suggests stronger affinities between Guianese and Caribbean coastal waters by the Last Interglacial than today, questioning the structuring role of the Amazon Plume on tropical Western Atlantic communities at the time. Grassland-dominated pollen, phytoliths, and charcoals from younger deposits in the same sections attest to a marine retreat and dryer conditions during the onset of the last glacial (ca. 110 to 50 ka), with a savanna-dominated landscape and episodes of fire. Charcoals from the last millennia suggest human presence in a mosaic of modern-like continental habitats. Our results provide key information about the ecology and biogeography of pristine Pleistocene tropical coastal ecosystems, especially relevant regarding the-widely anthropogenic-ongoing global warming.

Deep-pelagic fishes: Demographic instability in a stable environment

Demographic histories are frequently a product of the environment, as populations expand or contract in response to major environmental changes, often driven by changes in climate. Meso- and bathy-pelagic fishes inhabit some of the most temporally and spatially stable habitats on the planet. The stability of the deep-pelagic could make deep-pelagic fishes resistant to the demographic instability commonly reported in fish species inhabiting other marine habitats, however the demographic histories of deep-pelagic fishes are unknown. We reconstructed the historical demography of 11 species of deep-pelagic fishes using mitochondrial and nuclear DNA sequence data. We uncovered widespread evidence of population expansions in our study species, a counterintuitive result based on the nature of deep-pelagic ecosystems. Frequency-based methods detected potential demographic changes in nine species of fishes, while extended Bayesian skyline plots identified population expansions in four species. These results suggest that despite the relatively stable nature of the deep-pelagic environment, the fishes that reside here have likely been impacted by past changes in climate. Further investigation is necessary to better understand how deep-pelagic fishes, by far Earth's most abundant vertebrates, will respond to future climatic changes.

Whole genome assembly and annotation of the King Angelfish (Holacanthus passer) gives insight into the evolution of marine fishes of the Tropical Eastern Pacific

Holacanthus angelfishes are some of the most iconic marine fishes of the Tropical Eastern Pacific (TEP). However, very limited genomic resources currently exist for the genus. In this study we: (i) assembled and annotated the nuclear genome of the King Angelfish (Holacanthus passer), and (ii) examined the demographic history of H. passer in the TEP. We generated 43.8 Gb of ONT and 97.3 Gb Illumina reads representing 75× and 167× coverage, respectively. The final genome assembly size was 583 Mb with a contig N50 of 5.7 Mb, which captured 97.5% of the complete Actinoterygii Benchmarking Universal Single-Copy Orthologs (BUSCOs). Repetitive elements accounted for 5.09% of the genome, and 33,889 protein-coding genes were predicted, of which 22,984 were functionally annotated. Our demographic analysis suggests that population expansions of H. passer occurred prior to the last glacial maximum (LGM) and were more likely shaped by events associated with the closure of the Isthmus of Panama. This result is surprising, given that most rapid population expansions in both freshwater and marine organisms have been reported to occur globally after the LGM. Overall, this annotated genome assembly provides a novel molecular resource to study the evolution of Holacanthus angelfishes, while facilitating research into local adaptation, speciation, and introgression in marine fishes.

A marine record of Patagonian ice sheet changes over the past 140,000 years

Terrestrial glacial records from the Patagonian Andes and New Zealand Alps document quasi-synchronous Southern Hemisphere-wide glacier advances during the late Quaternary. However, these records are inherently incomplete. Here, we provide a continuous marine record of western-central Patagonian ice sheet (PIS) extent over a complete glacial-interglacial cycle back into the penultimate glacial (~140 ka). Sediment core MR16-09 PC03, located at 46°S and ~150 km offshore Chile, received high terrestrial sediment and meltwater input when the central PIS extended westward. We use biomarkers, foraminiferal oxygen isotopes, and major elemental data to reconstruct terrestrial sediment and freshwater input related to PIS variations. Our sediment record documents three intervals of general PIS marginal fluctuations, during Marine Isotope Stage (MIS) 6 (140 to 135 ka), MIS 4 (~70 to 60 ka), and late MIS 3 to MIS 2 (~40 to 18 ka). These higher terrigenous input intervals occurred during sea-level low stands, when the western PIS covered most of the Chilean fjords, which today retain glaciofluvial sediments. During these intervals, high-amplitude phases of enhanced sediment supply occur at millennial timescales, reflecting increased ice discharge most likely due to a growing PIS. We assign the late MIS 3 to MIS 2 phases and, by inference, older advances to Antarctic cold stages. We conclude that the increased sediment/meltwater release during Southern Hemisphere millennial-scale cold phases was likely related to higher precipitation caused by enhanced westerly winds at the northwestern margin of the PIS. Our records complement terrestrial archives and provide evidence for PIS climate sensitivity.

Hydrogeochemical and environmental isotope study of Topusko thermal waters, Croatia

Thermal waters in Topusko (Croatia), with temperatures of up to 65 °C, have been used for heating, health, and recreational tourism for the past fifty years. Hydrogeochemical monitoring can provide insights into deeper geological processes and indicate system changes from baseline levels. It helps to identify potential anthropogenic impacts, as well as natural changes. Hydrogeochemical, geothermometrical, and environmental isotope studies of thermal waters in Topusko were conducted to improve the existing conceptual model of the Topusko hydrothermal system (THS), providing a baseline for continuous monitoring of the thermal resource. 2-year thermal springs and precipitation monitoring took place from March 2021 until March 2023. Major anions and cations, stable and radioactive isotopes (i.e. 18O, 2H, SO42-, 3H and 14C) and geothermometers were used to assess the origin of thermal waters in Topusko and their interaction with thermal aquifer. The results indicate the meteoric origin of thermal water, which was recharged in colder climatic conditions around the late Pleistocene-Early Holocene. Thermal water was last in contact with the atmosphere before approximately 9.5 kyr. Ca-HCO3 hydrochemical facies suggests carbonate dissolution as the dominant process driving the solute content. Geothermometrical results indicate an equilibrium temperature in the reservoir of 90 °C.

A renewed glance at the Palearctic golden eagle: Genetic variation in space and time

Anthropogenic pressures on nature have been causing population declines for centuries. Intensified persecution of apex predators, like the golden eagle, resulted in population bottlenecks during the 19th and 20th centuries. To study population genetics and demographic history of the golden eagle throughout its distribution, we collected museum samples from previously underrepresented regions, such as Russia and Central Asia. We used 12 microsatellite loci and a fragment of the mitochondrial DNA control region to re-evaluate phylogeography of Eurasian golden eagles and study the impacts of the population bottleneck. Our results revealed a north-south genetic gradient, expressed by the difference between Mediterranean and Holarctic lineages, as well as genetically distinct Northern Europe and Central Asia and Caucasus regions. Furthermore, Northern Europe exhibited the lowest, whereas Central Asia and Caucasus had the highest genetic diversity. Although golden eagles maintained relatively high genetic diversity, we detected genetic signatures of the recent bottleneck, including reduced genetic diversity and a decline in the effective female population size around the year 1975. Our study improves the knowledge of the genetic composition of Eurasian golden eagles and highlights the importance of understanding their historical population dynamics in the face of ongoing and future conservation efforts.

Riverscape genetics of the orangethroat darter complex

Freshwater darters belonging to the orangethroat darter species complex, or Ceasia, are widely distributed in the Central and Southern United States, with ranges that span both glaciated and unglaciated regions. Up to 15 species have been recognized in the complex, with one, Etheostoma spectabile, having a widespread northern distribution and another, Etheostoma pulchellum, having a sizeable southern distribution. The other species in the complex have much more restricted distributions in unglaciated regions of the Central Highlands. We sampled 384 darters from 52 sites covering much of the range of Ceasia and evaluated patterns of genetic diversity, genetic structure, and pre- and post-glacial patterns of range contraction and expansion. We anticipated finding much stronger signals of genetic differentiation and diversification in unglaciated regions, given the higher species diversity and levels of endemism reported there. Surprisingly, microsatellite genotyping revealed two well-differentiated genetic clusters of E. spectabile in samples from glaciated regions, one confined to the Illinois River basin and another found in the Wabash drainage and Great Lakes tributaries. This suggests that there was expansion from two isolated glacial refugia, with little subsequent post-glacial gene flow. Fish collected from throughout the unglaciated region were less genetically differentiated. Fish assigned to Etheostoma burri and Etheostoma uniporum based on collection sites and morphological characters were not genetically differentiated from E. spectabile samples from the region. Hybridization and introgression occurring in the Central Highlands may confound genetic delineation of species in this region of high endemism and diversity.

Genomic Analyses Capture the Human-Induced Demographic Collapse and Recovery in a Wide-Ranging Cervid

The glacial cycles of the Quaternary heavily impacted species through successions of population contractions and expansions. Similarly, populations have been intensely shaped by human pressures such as unregulated hunting and land use changes. White-tailed and mule deer survived in different refugia through the Last Glacial Maximum, and their populations were severely reduced after the European colonization. Here, we analyzed 73 resequenced deer genomes from across their North American range to understand the consequences of climatic and anthropogenic pressures on deer demographic and adaptive history. We found strong signals of climate-induced vicariance and demographic decline; notably, multiple sequentially Markovian coalescent recovers a severe decline in mainland white-tailed deer effective population size (Ne) at the end of the Last Glacial Maximum. We found robust evidence for colonial overharvest in the form of a recent and dramatic drop in Ne in all analyzed populations. Historical census size and restocking data show a clear parallel to historical Ne estimates, and temporal Ne/Nc ratio shows patterns of conservation concern for mule deer. Signatures of selection highlight genes related to temperature, including a cold receptor previously highlighted in woolly mammoth. We also detected immune genes that we surmise reflect the changing land use patterns in North America. Our study provides a detailed picture of anthropogenic and climatic-induced decline in deer diversity and clues to understanding the conservation concerns of mule deer and the successful demographic recovery of white-tailed deer.

New age constraints for human entry into the Americas on the north Pacific coast

The timing of the initial peopling of the Americas is unresolved. Because the archaeological record necessitates discussion of human entry from Beringia into southern North America during the last glaciation, addressing this problem routinely involves evaluating environmental parameters then targeting areas suitable for human settlement. Vertebrate remains indicate landscape quality and are a key dataset for assessing coastal migration theories and the viability of coastal routes. Here, radiocarbon dates on vertebrate specimens and archaeological sites are calibrated to document species occurrences and the ages of human settlements across the western expansion and decay of the Cordilleran Ice Sheet (CIS) during the Late Wisconsin Fraser Glaciation in four subregions of the north Pacific coast of North America. The results show archaeological sites occur after glacial maxima and are generally consistent with the age of other securely dated earliest sites in southern North America. They also highlight gaps in the vertebrate chronologies around CIS maxima in each of the subregions that point to species redistributions and extirpations and signal times of low potential for human settlement and subsistence in a key portion of the proposed coastal migration route. This study, therefore, defines new age constraints for human coastal migration theories in the peopling of the Americas debate.

The genomes of Hercules beetles reveal putative adaptive loci and distinct demographic histories in pristine North American forests

Beetles, despite their remarkable biodiversity and a long history of research, remain lacking in reference genomes annotated with structural variations in loci of adaptive significance. We sequenced and assembled high-quality chromosome-level genomes of four Hercules beetles which exhibit divergence in male horn size and shape and body colouration. The four Hercules beetle genomes were assembled to 11 pseudo-chromosomes, where the three genomes assembled using Nanopore data (Dynastes grantii, D. hyllus and D. tityus) were mapped to the genome assembled using PacBio + Hi-C data (D. maya). We demonstrated a striking similarity in genome structure among the four species. This conservative genome structure may be attributed to our use of the D. maya assembly as the reference; however, it is worth noting that such a conservative genome structure is a recurring phenomenon among scarab beetles. We further identified homologues of nine and three candidate-gene families that may be associated with the evolution of horn structure and body colouration respectively. Structural variations in Scr and Ebony2 were detected and discussed for their putative impacts on generating morphological diversity in beetles. We also reconstructed the demographic histories of the four Hercules beetles using heterozygosity information from the diploid genomes. We found that the demographic histories of the beetles closely recapitulated historical changes in suitable forest habitats driven by climate shifts.

A chromosome-level genome of electric catfish (Malapterurus electricus) provided new insights into order Siluriformes evolution

The electric catfish (Malapterurus electricus), belonging to the family Malapteruridae, order Siluriformes (Actinopterygii: Ostariophysi), is one of the six branches that has independently evolved electrical organs. We assembled a 796.75 Mb M. electricus genome and anchored 88.72% sequences into 28 chromosomes. Gene family analysis revealed 295 expanded gene families that were enriched on functions related to glutamate receptors. Convergent evolutionary analyses of electric organs among different lineage of electric fishes further revealed that the coding gene of rho guanine nucleotide exchange factor 4-like (arhgef4), which is associated with G-protein coupled receptor (GPCR) signaling pathway, underwent adaptive parallel evolution. Gene identification suggests visual degradation in catfishes, and an important role for taste in environmental adaptation. Our findings fill in the genomic data for a branch of electric fish and provide a relevant genetic basis for the adaptive evolution of Siluriformes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s42995-023-00197-8.

Millennial-scale variability of Greenland dust provenance during the last glacial maximum as determined by single particle analysis

Greenland ice core records exhibited 100-fold higher dust concentrations during the Last Glacial Maximum (LGM) than during the Holocene, and dust input temporal variability corresponded to different climate states in the LGM. While East Asian deserts, the Sahara, and European loess have been suggested as the potential source areas (PSAs) for Greenland LGM dust, millennial-scale variability in their relative contributions within the LGM remains poorly constrained. Here, we present the morphological, mineralogical, and geochemical characteristics of insoluble microparticles to constrain the provenance of dust in Greenland NEEM ice core samples covering cold Greenland Stadials (GS)-2.1a to GS-3 (~ 14.7 to 27.1 kyr ago) in the LGM. The analysis was conducted on individual particles in microdroplet samples by scanning electron microscopy with energy dispersive X-ray spectroscopy and Raman microspectroscopy. We found that the kaolinite-to-chlorite (K/C) ratios and chemical index of alteration (CIA) values were substantially higher (K/C: 1.4 ± 0.7, CIA: 74.7 ± 2.9) during GS-2.1a to 2.1c than during GS-3 (K/C: 0.5 ± 0.1, CIA: 65.8 ± 2.8). Our records revealed a significant increase in Saharan dust contributions from GS-2.1a to GS-2.1c and that the Gobi Desert and/or European loess were potential source(s) during GS-3. This conclusion is further supported by distinctly different carbon contents in particles corresponding to GS-2.1 and GS-3. These results are consistent with previous estimates of proportional dust source contributions obtained using a mixing model based on Pb and Sr isotopic compositions in NEEM LGM ice and indicate millennial-scale changes in Greenland dust provenance that are probably linked to large-scale atmospheric circulation variabilities during the LGM.

Genomic evidence of genuine wild versus admixed olive populations evolving in the same natural environments in western Mediterranean Basin

Crop-to-wild gene flow is a mechanism process widely documented, both in plants and animals. This can have positive or negative impacts on the evolution of admixed populations in natural environments, yet the phenomenon is still misunderstood in long-lived woody species, contrary to short-lived crops. Wild olive Olea europaea L. occurs in the same eco-geographical range as domesticated olive, i.e. the Mediterranean Basin (MB). Moreover, it is an allogamous and anemophilous species whose seeds are disseminated by birds, i.e. factors that drive gene flow between crops and their wild relatives. Here we investigated the genetic structure of western MB wild olive populations in natural environments assuming a homogenous gene pool with limited impact of cultivated alleles, as previously suggested. We used a target sequencing method based on annotated genes from the Farga reference genome to analyze 27 western MB olive tree populations sampled in natural environments in France, Spain and Morocco. We also target sequenced cultivated olive tree accessions from the Worldwide Olive Germplasm Bank of Marrakech and Porquerolles and from an eastern MB wild olive tree population. We combined PCA, sNMF, pairwise FST and TreeMix and clearly identified genuine wild olive trees throughout their natural distribution range along a north-south gradient including, for the first time, in southern France. However, contrary to our assumption, we highlighted more admixed than genuine wild olive trees. Our results raise questions regarding the admixed population evolution pattern in this environment, which might be facilitated by crop-to-wild gene flow.

Plankton community changes during the last 124 000 years in the subarctic Bering Sea derived from sedimentary ancient DNA

Current global warming results in rising sea-water temperatures, and the loss of sea ice in Arctic and subarctic oceans impacts the community composition of primary producers with cascading effects on the food web and potentially on carbon export rates. This study analyzes metagenomic shotgun and diatom rbcL amplicon sequencing data from sedimentary ancient DNA of the subarctic western Bering Sea that records phyto- and zooplankton community changes over the last glacial-interglacial cycles, including the last interglacial period (Eemian). Our data show that interglacial and glacial plankton communities differ, with distinct Eemian and Holocene plankton communities. The generally warm Holocene period is dominated by picosized cyanobacteria and bacteria-feeding heterotrophic protists, while the Eemian period is dominated by eukaryotic picosized chlorophytes and Triparmaceae. By contrast, the glacial period is characterized by microsized phototrophic protists, including sea ice-associated diatoms in the family Bacillariaceae and co-occurring diatom-feeding crustaceous zooplankton. Our deep-time record of plankton community changes reveals a long-term decrease in phytoplankton cell size coeval with increasing temperatures, resembling community changes in the currently warming Bering Sea. The phytoplankton community in the warmer-than-present Eemian period is distinct from modern communities and limits the use of the Eemian as an analog for future climate scenarios. However, under enhanced future warming, the expected shift toward the dominance of small-sized phytoplankton and heterotrophic protists might result in an increased productivity, whereas the community's potential of carbon export will be decreased, thereby weakening the subarctic Bering Sea's function as an effective carbon sink.

Simulating long-term wildfire impacts on boreal forest structure in Central Yakutia, Siberia, since the Last Glacial Maximum

Background

Wildfires are recognized as an important ecological component of larch-dominated boreal forests in eastern Siberia. However, long-term fire-vegetation dynamics in this unique environment are poorly understood. Recent paleoecological research suggests that intensifying fire regimes may induce millennial-scale shifts in forest structure and composition. This may, in turn, result in positive feedback on intensifying wildfires and permafrost degradation, apart from threatening human livelihoods. Most common fire-vegetation models do not explicitly include detailed individual-based tree population dynamics, but a focus on patterns of forest structure emerging from interactions among individual trees may provide a beneficial perspective on the impacts of changing fire regimes in eastern Siberia. To simulate these impacts on forest structure at millennial timescales, we apply the individual-based, spatially explicit vegetation model LAVESI-FIRE, expanded with a new fire module. Satellite-based fire observations along with fieldwork data were used to inform the implementation of wildfire occurrence and adjust model parameters.

Results

Simulations of annual forest development and wildfire activity at a study site in the Republic of Sakha (Yakutia) since the Last Glacial Maximum (c. 20,000 years BP) highlight the variable impacts of fire regimes on forest structure throughout time. Modeled annual fire probability and subsequent burned area in the Holocene compare well with a local reconstruction of charcoal influx in lake sediments. Wildfires can be followed by different forest regeneration pathways, depending on fire frequency and intensity and the pre-fire forest conditions. We find that medium-intensity wildfires at fire return intervals of 50 years or more benefit the dominance of fire-resisting Dahurian larch (Larix gmelinii (Rupr.) Rupr.), while stand-replacing fires tend to enable the establishment of evergreen conifers. Apart from post-fire mortality, wildfires modulate forest development mainly through competition effects and a reduction of the model's litter layer.

Conclusion

With its fine-scale population dynamics, LAVESI-FIRE can serve as a highly localized, spatially explicit tool to understand the long-term impacts of boreal wildfires on forest structure and to better constrain interpretations of paleoecological reconstructions of fire activity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42408-023-00238-8.

The chromosome-scale genome of Magnolia sinica (Magnoliaceae) provides insights into the conservation of plant species with extremely small populations (PSESP)

Magnolia sinica (Magnoliaceae) is a highly threatened tree endemic to southeast Yunnan, China. In this study, we generated for the first time a high-quality chromosome-scale genome sequence from M. sinica, by combining Illumina and ONT data with Hi-C scaffolding methods. The final assembled genome size of M. sinica was 1.84 Gb, with a contig N50 of ca. 45 Mb and scaffold N50 of 92 Mb. Identified repeats constituted approximately 57% of the genome, and 43,473 protein-coding genes were predicted. Phylogenetic analysis shows that the magnolias form a sister clade with the eudicots and the order Ceratophyllales, while the monocots are sister to the other core angiosperms. In our study, a total of 21 individuals from the 5 remnant populations of M. sinica, as well as 22 specimens belonging to 8 related Magnoliaceae species, were resequenced. The results showed that M. sinica had higher genetic diversity (θw = 0.01126 and θπ = 0.01158) than other related species in the Magnoliaceae. However, population structure analysis suggested that the genetic differentiation among the 5 M. sinica populations was very low. Analyses of the demographic history of the species using different models consistently revealed that 2 bottleneck events occurred. The contemporary effective population size of M. sinica was estimated to be 10.9. The different patterns of genetic loads (inbreeding and numbers of deleterious mutations) suggested constructive strategies for the conservation of these 5 different populations of M. sinica. Overall, this high-quality genome will be a valuable genomic resource for conservation of M. sinica.

Sedimentary DNA reveals the link between microbial community dynamics and climate during the late last glaciation in the offshore region of the Zambezi River, Southwest Indian Ocean

Reconstructing the relationship between microbial communities and past abrupt climate change is of great importance for understanding current biodiversity patterns and predicting changes under future climate scenarios. However, little is currently known about how microbial communities respond to changes in key environmental stages due to a lack of research in this area. Here, we examine the variability in the communities of bacteria, archaea, and fungi from sediments deposited offshore region of the Zambezi River between 21.7 and 9.6 thousand years ago (ka) (covering the last glacial maximum, or LGM, and the early Holocene) using DNA metabarcoding approach via high-throughput sequencing. The results showed that (1) microbial assemblages differed across three key time intervals, with the last deglaciation having the most homogeneous prokaryotic assemblages, while for fungal communities in the LGM, and the early Holocene and LGM differing the most; (2) the warm early Holocene showed the highest diversity, whereas the lowest diversity was found in the LGM; and (3) the selected indicator species better reflected the climatic characteristics of different environmental stages. These results highlight the power of ancient sedimentary DNA to refine our understanding of microbial dynamics in marine sedimentary systems near large rivers, thus providing a basis for better modeling ecological processes in further research.

Population Structure and Genetic Diversity Analyses Provide New Insight into the Endemic Species Aster spathulifolius Maxim. and Its Evolutionary History

Aster spathulifolius, an ecologically significant plant species native to the coastal regions of Korea and Japan, remains understudied in terms of its genetic structure and evolutionary history. In this study, we employed four chloroplast markers and the nuclear ITS region from 15 populations of A. spathulifolius from both Korea and Japan, including their islands, to unravel the spatial genetic structure, differentiation, gene flow, phylogenetic, and biogeographical relationships. Analysis based on multiple methods identified a low level of genetic diversity, genetic differentiation and gene flow among A. spathulifolius populations. Network analysis and principal coordinates analysis showed that 15 populations could be divided into two groups: mainland and island. Furthermore, UPGMA, neighbor-net, maximum-likelihood and Bayesian inference-based phylogenetic tree confirmed that these populations formed two distinct clades. Therefore, the island populations might be treated as A. spathulifolius populations rather than A. oharai populations. Divergence time analysis estimated the divergence of A. spathulifolius lineages approximately 23.09 million years ago, while ancestral area reconstruction analysis suggested Korea as the potential origin, conflicting with alternative scenarios. These findings contribute to a comprehensive understanding of the evolutionary history, genetic structure, and adaptive strategies of A. spathulifolius in coastal environments. Our study challenges previous assumptions and underscores the necessity for further population studies to elucidate the intricate dynamics of this distinctive plant species.

Millennial processes of population decline, range contraction and near extinction of the European bison

European bison (Bison bonasus) were widespread throughout Europe during the late Pleistocene. However, the contributions of environmental change and humans to their near extinction have never been resolved. Using process-explicit models, fossils and ancient DNA, we disentangle the combinations of threatening processes that drove population declines and regional extinctions of European bison through space and across time. We show that the population size of European bison declined abruptly at the termination of the Pleistocene in response to rapid environmental change, hunting by humans and their interaction. Human activities prevented populations of European bison from rebounding in the Holocene, despite improved environmental conditions. Hunting caused range loss in the north and east of its distribution, while land use change was responsible for losses in the west and south. Advances in hunting technologies from 1500 CE were needed to simulate low abundances observed in 1870 CE. While our findings show that humans were an important driver of the extinction of the European bison in the wild, vast areas of its range vanished during the Pleistocene-Holocene transition because of post-glacial environmental change. These areas of its former range have been climatically unsuitable for millennia and should not be considered in reintroduction efforts.

Evidence for a selective link between cooperation and individual recognition

The ability to recognize others is a frequent assumption of models of the evolution of cooperation. At the same time, cooperative behavior has been proposed as a selective agent favoring the evolution of individual recognition abilities. Although theory predicts that recognition and cooperation may co-evolve, data linking recognition abilities and cooperative behavior with evidence of selection are elusive. Here, we provide evidence of a selective link between individual recognition and cooperation in the paper wasp Polistes fuscatus through a combination of clinal, common garden, and population genomics analyses. We identified latitudinal clines in both rates of cooperative nesting and color pattern diversity, consistent with a selective link between recognition and cooperation. In behavioral experiments, we replicated previous results demonstrating individual recognition in cooperative and phenotypically diverse P. fuscatus from New York. In contrast, wasps from a less cooperative and phenotypically uniform Louisiana population showed no evidence of individual recognition. In a common garden experiment, groups of wasps from northern populations formed more stable and individually biased associations, indicating that recognition facilitates group stability. The strength of recent positive selection on cognition-associated loci likely to mediate individual recognition is substantially greater in northern compared with southern P. fuscatus populations. Collectively, these data suggest that individual recognition and cooperative nesting behavior have co-evolved in P. fuscatus because recognition helps stabilize social groups. This work provides evidence of a specific cognitive phenotype under selection because of social interactions, supporting the idea that social behavior can be a key driver of cognitive evolution.

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.

Dispersal from the Qinghai-Tibet plateau by a high-altitude butterfly is associated with rapid expansion and reorganization of its genome

Parnassius glacialis is a typical "Out of the QTP" alpine butterfly that originated on the Qinghai-Tibet Plateau (QTP) and dispersed into relatively low-altitude mountainous. Here we assemble a chromosome-level genome of P. glacialis and resequence 9 populations in order to explore the genome evolution and local adaptation of this species. These results indicated that the rapid accumulation and slow unequal recombination of transposable elements (TEs) contributed to the formation of its large genome. Several ribosomal gene families showed extensive expansion and selective evolution through transposon-mediated processed pseudogenes. Additionally, massive structural variations (SVs) of TEs affected the genetic differentiation of low-altitude populations. These low-altitude populations might have experienced a genetic bottleneck in the past and harbor genes with selective signatures which may be responsible for the potential adaptation to low-altitude environments. These results provide a foundation for understanding genome evolution and local adaptation for "Out of the QTP" of P. glacialis.

Unraveling the complex phylogeographic history of freshwater fishes in Lower Central America: A study of the electric fish Brachyhypopomus occidentalis

Lower Central America (LCA) has a complex biogeographic history shaped by the rise of the Isthmus of Panama and the global climatic oscillations of the Pleistocene. These events have been crucial in structuring biodiversity in LCA, but their consequences for the distribution and partitions of genetic diversity across the region remain to be elucidated. We combined complete mitochondrial genomes and nuclear ultraconserved elements (UCEs) to study the phylogeographic history and population genetic structure of the electric fish Brachyhypopomus occidentalis in LCA. Our results are consistent with the known phylogeographic history of B. occidentalis in LCA, but we update this history in several important ways that help illuminate the phylogeographic history of freshwater fishes in the region. We provide: i) support for three waves of colonization, two of which occurred prior to the final closure of the Panama Isthmus; ii) a more precise understanding of each colonization event, with evidence for a larger footprint of the first event, as well as genetic exchange across the continental divide in subsequent events; and iii) evidence for high levels of previously unrecognized population genetic structure across LCA. This updated model of colonization and diversification of B. occidentalis consists of three waves of dispersal and colonization, which triggered the evolution of geographic breaks in both nuclear and mitochondrial genomes across LCA. These processes are tightly linked to the dynamic uplift of the Isthmus, recent volcanic activity in the region, and the sea-level oscillations of the Pleistocene. These results improve previous phylogeographic inferences regarding the distribution and diversification of freshwater fishes in LCA, and generate testable hypotheses to guide future research exploring the factors shaping biodiversity in the region.

Approximate Bayesian computation and ecological niche models elucidate the demographic history and current fragmented population distribution of a Korean endemic shrub

Climatic fluctuations and geological events since the LGM are believed to have significantly impacted the population size, distribution, and mobility of many species that we observe today. In this paper, we determined the processes driving the phylogeographic structure of the Korean endemic white forsythia by combining the use of genome-wide SNPs and predicting paleoclimatic habitats during the LGM (21 kya), Early Holocene (10 kya), Mid-Holocene (6 kya), and Late Holocene (3 kya). Using a maximum of 1897 SNPs retrieved from 124 samples across nine wild populations, five environmental predictors, and the species' natural occurrence records, we aimed to infer the species' demographic history and reconstruct its possible paleodistributions with the use of approximate Bayesian computation and ecological niche models, respectively. Under this integrated framework, we found strong evidence for patterns of range shift and expansion, and population divergence events from the onset of the Holocene, resulting in the formation of its five distinct genetic units. The most highly supported model inferred that after the split of an ancestral population into the southern group and a larger central metapopulation lineage, the latter gave rise to the eastern and northern clusters, before finally dividing into two sub-central groups. While the use of molecular data allowed us to identify and refine the (phylo)genetic relationships of the species' lineages and populations, the use of ecological data helped us infer a past LGM refugium and the directions of post-glacial range dynamics. The time frames of these demographic events were shown to be congruent with climatic and geological events that affected the central Korean Peninsula during these periods. These findings gave us a better understanding of the consequences of past spatiotemporal factors that may have resulted in the current fragmented population distribution of this endangered plant.

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.

Unraveling the mitochondrial phylogenetic landscape of Thailand reveals complex admixture and demographic dynamics

The evolutionary dynamics of mitochondrial DNA within the Thai population were comprehensively explored with a specific focus on the influence of South Asian admixture. A total of 166 samples were collected through randomized sampling, ensuring a diverse representation. Our findings unveil substantial genetic and haplogroup diversity within the Thai population. We have identified 164 haplotypes categorized into 97 haplogroups, with a notable inclusion of 20 novel haplogroups. The distribution of haplogroups exhibited variations across different populations and countries. The central Thai population displayed a high diversity of haplogroups from both the M and N clades. Maternal lineage affinities were discerned between several Mainland Southeast Asia (MSEA) and South Asian populations, implying ancestral genetic connections and a substantial influence of South Asian women in establishing these relationships. f4-statistics indicates the presence of a Tibeto-Burman genetic component within the Mon population from Thailand. New findings demonstrate two phases of population expansion occurring 22,000-26,000 and 2500-3800 years ago, coinciding with the Last Glacial Maximum, and Neolithic demographic transition, respectively. This research significantly enhances our understanding of the maternal genetic history of Thailand and MSEA, emphasizing the influence of South Asian admixture. Moreover, it underscores the critical role of prior information, such as mutation rates, within the Bayesian framework for accurate estimation of coalescence times and inferring demographic history.

The maternal ancestry of the Kavaratti islanders and the last glacial maximum aftermath

The prehistoric human settlement of the Lakshadweep islands remains a mystery for various reasons. Uncertainty about the existence of indigenous tribes in these islands and the lack of folklore records present major obstacles to the reconstruction of Lakshadweep ancestry. However, with extant population data, we seek to understand the maternal ancestry of the Kavaratti islanders. Mitochondrial control region variation analysis of 80 individuals from this island shows maternal links with the populations in the northwestern region of the South Asian mainland. The founder clade R30b2, observed in the Kavaratti islanders, is so far present only in the Scheduled Castes from the Punjab region, Jat Sikhs and Nairs. All other mainland populations carry basal R30 or R30a subclades. The presence of a specific Uralic U4 lineage in our samples, in addition to the Indo-European affinity observed in the phylogeny tree, substantiates a northwestern maternal ancestry of the Kavaratti islanders and implies an ancestral admixture with early humans in the Near East at the time of the last glacial maximum (LGM). Based on our Bayesian analysis, we furthermore propose that a group bearing mostly R30b2 during the LGM recovery, moved eastward and southward, where they received Indian-specific M haplogroups. Hence, the maternal ancestry of the Kavaratti islanders is evidently a consequence of the demographic changes in the northwestern region of the Indian subcontinent caused by the Last Glacial Maximum. The haplogroup distribution pattern and nucleotide sequence data produced in this study will enrich the forensic database of the Lakshadweep islands.

1-km resolution rebound surfaces and paleotopography of glaciated North America since the Last Glacial Maximum

We present a series of 1-km spatial resolution rebound (isobase) surfaces based on publicly distributed predictions obtained from the glacio-isostatic adjustment models known as ICE-5G (VM2 L90), ICE-6G_C (VM5a) and ICE-7G_NA (VM7). Our objective is to provide readily accessible tools for a broad range of geological and paleoenvironmental studies, and to facilitate direct comparison between models' predictions and field-based observations. Rebound surfaces were interpolated at the scale of North American ice sheets (35.5°-89.5°N; 45°-165°W) and for each time increment of the models (1,000-500 yrs, between 26,000-21,000 yrs BP and present-day). The assessment of the interpolations indicates that the rebound surfaces have an overall vertical accuracy of ∼0.4 m compared to original ICE-xG outputs. These rebound surfaces were combined with the GEBCO 2021 present-day elevation grid to reconstruct the paleotopography for each time increment of the models and are all presented as raster files that can be easily integrated into geographical information systems. The resulting datasets therefore provide a unique support for geological, paleoenvironmental and archeological studies.

Phylogenomic insights into the polyphyletic nature of Altai falcons within eastern sakers (Falco cherrug) and the origins of gyrfalcons (Falco rusticolus)

The Altai falcon from Central Asia always attracted the attention of humans. Long considered a totemic bird in its native area, modern falconers still much appreciated this large-bodied and mighty bird of prey due to its rarity and unique look. The peculiar body characteristics halfway between the saker falcon (Falco cherrug) and the gyrfalcon (F. rusticolus) triggered debates about its contentious taxonomy. The weak phylogenetic signal associated with traditional genetic methods could not resolve this uncertainty. Here, we address the controversial evolutionary origin of Altai falcons by means of a genome-wide approach, Restriction-site Associated DNA sequencing, using sympatric eastern sakers falcons, allopatric western saker falcons and gyrfalcons as outgroup. This approach provided an unprecedented insight into the phylogenetic relationships of the studied populations by delivering 17,095 unlinked SNPs shedding light on the polyphyletic nature of Altai falcons within eastern sakers. Thus we concluded that the former must correspond to a low taxonomic rank, probably an ecotype or form of the latter. Also, we found that eastern sakers are paraphyletic without gyrfalcons, thus, these latter birds are best regarded as the direct sister lineage of the eastern sakers. This evolutionary relationship, corroborated also by re-analyzing the dataset with the inclusion of outgroup samples (F. biarmicus and F. peregrinus), put eastern sakers into a new light as the potential ancestral genetic source of high latitude and altitude adaptation in descendent populations. Finally, conservation genomic values hint at the stable genetic background of the studied saker populations.

Early sea ice decline off East Antarctica at the last glacial-interglacial climate transition

Antarctic climate warming and atmospheric CO2 rise during the last deglaciation may be attributed in part to sea ice reduction in the Southern Ocean. Yet, glacial-interglacial Antarctic sea ice dynamics and underlying mechanisms are poorly constrained, as robust sea ice proxy evidence is sparse. Here, we present a molecular biomarker-based sea ice record that resolves the spring/summer sea ice variability off East Antarctica during the past 40 thousand years (ka). Our results indicate that substantial sea ice reduction culminated rapidly and contemporaneously with upwelling of carbon-enriched waters in the Southern Ocean at the onset of the last deglaciation but began at least ~2 ka earlier probably driven by an increasing local integrated summer insolation. Our findings suggest that sea ice reduction and associated feedbacks facilitated stratification breakup and outgassing of CO2 in the Southern Ocean and warming in Antarctica but may also have played a leading role in initializing these deglacial processes in the Southern Hemisphere.

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.

Independent age estimates resolve the controversy of ancient human footprints at White Sands

Human footprints at White Sands National Park, New Mexico, USA, reportedly date to between ~23,000 and 21,000 years ago according to radiocarbon dating of seeds from the aquatic plant Ruppia cirrhosa. These ages remain controversial because of potential old carbon reservoir effects that could compromise their accuracy. We present new calibrated 14C ages of terrestrial pollen collected from the same stratigraphic horizons as those of the Ruppia seeds, along with optically stimulated luminescence ages of sediments from within the human footprint-bearing sequence, to evaluate the veracity of the seed ages. The results show that the chronologic framework originally established for the White Sands footprints is robust and reaffirm that humans were present in North America during the Last Glacial Maximum.

Genetic Origins and Adaptive Evolution of the Deng People on the Tibetan Plateau

The Tibetan Plateau is populated by diverse ethnic groups, but most of them are underrepresented in genomics studies compared with the Tibetans (TIB). Here, to gain further insight into the genetic diversity and evolutionary history of the people living in the Tibetan Plateau, we sequenced 54 whole genomes of the Deng people with high coverage (30-60×) and analyzed the data together with that of TIB and Sherpas, as well as 968 ancient Asian genomes and available archaic and modern human data. We identified 17.74 million novel single-nucleotide variants from the newly sequenced genomes, although the Deng people showed reduced genomic diversity and a relatively small effective population size. Compared with the other Tibetan highlander groups which are highly admixed, the Deng people are dominated by a sole ancestry that could be traced to some ancient northern East Asian populations. The divergence between Deng and Tibetan people (∼4,700-7,200 years) was more recent than that between highlanders and the Han Chinese (Deng-HAN, ∼9,000-14,000 years; TIB-HAN, 7,200-10,000 years). Adaptive genetic variants (AGVs) identified in the Deng are only partially shared with those previously reported in the TIB like HLA-DQB1, whereas others like KLHL12 were not reported in TIB. In contrast, the top candidate genes harboring AGVs as previously identified in TIB, like EPAS1 and EGLN1, do not show strong positive selection signals in Deng. Interestingly, Deng also showed a different archaic introgression scenario from that observed in the TIB. Our results suggest that convergent adaptation might be prevalent on the Tibetan Plateau.

Human occupations of upland and cold environments in inland Spain during the Last Glacial Maximum and Heinrich Stadial 1: The new Magdalenian sequence of Charco Verde II

The settlement of cold and arid environments by Pleistocene hunter-gatherers has been a heated topic in Paleolithic Archaeology and the Quaternary Sciences for years. In the Iberian Peninsula, a key area for studying human adaptations to such environments is composed by the large interior and upland regions of the northern and southern plateaus (Mesetas) and bordering areas. As, traditionally, these regions have been relatively under-investigated compared to the ecologically more favored coastal areas of the peninsula, our knowledge of the human settlement of the whole Iberian hinterland remains scarce for the Last Glacial. In this paper we present the discovery and first geoarcheological, paleoenvironmental and chronometric evidence obtained at Charco Verde II, a new site close to the southwestern foothills of the Iberian system range (Guadalajara province, Spain), bearing a sequence of Magdalenian human occupations starting at least at 20.8-21.4 ka cal BP during the Last Glacial Maximum, and covering Greenland Stadial 2 until ∼15.1-16.6 ka cal BP, including Heinrich stadial 1. As this site is located in an upland region which today faces one of the harshest climates in Iberia, such occupation sequence, occurred during some of the coldest and most arid phases of the Last Glacial, has relevant implications for our understanding of human-environment-climate interactions and population dynamics in Iberia and Western Europe. These findings support the hypothesis that the Iberian hinterland was not avoided by Upper Paleolithic hunter-gatherers due to ecological constraints, but it hosted a complex and relatively dense settlement at least in some areas, even during cold periods. This suggest, one more time, that the historical scarcity of Upper Paleolithic sites in inland Iberia is, to a significant extent, an artifact of research bias.

Deciphering complex reticulate evolution of Asian Buddleja (Scrophulariaceae): insights into the taxonomy and speciation of polyploid taxa in the Sino-Himalayan region

BACKGROUND AND AIMS

Species of the genus Buddleja in Asia are mainly distributed in the Sino-Himalayan region and form a challenging taxonomic group, with extensive hybridization and polyploidization. A phylogenetic approach to unravelling the history of reticulation in this lineage will deepen our understanding of the speciation in biodiversity hotspots.

METHODS

For this study, we obtained 80 accessions representing all the species in the Asian Buddleja clade, and the ploidy level of each taxon was determined by flow cytometry analyses. Whole plastid genomes, nuclear ribosomal DNA, single nucleotide polymorphisms and a large number of low-copy nuclear genes assembled from genome skimming data were used to investigate the reticulate evolutionary history of Asian Buddleja. Complex cytonuclear conflicts were detected through a comparison of plastid and species trees. Gene tree incongruence was also analysed to detect any reticulate events in the history of this lineage.

KEY RESULTS

Six hybridization events were detected, which are able to explain the cytonuclear conflict in Asian Buddleja. Furthermore, PhyloNet analysis combining species ploidy data indicated several allopolyploid speciation events. A strongly supported species tree inferred from a large number of low-copy nuclear genes not only corrected some earlier misinterpretations, but also indicated that there are many Asian Buddleja species that have been lumped mistakenly. Divergent time estimation shows two periods of rapid diversification (8-10 and 0-3 Mya) in the Asian Buddleja clade, which might coincide with the final uplift of the Hengduan Mountains and Quaternary climate fluctuations, respectively.

CONCLUSIONS

This study presents a well-supported phylogenetic backbone for the Asian Buddleja species, elucidates their complex and reticulate evolutionary history and suggests that tectonic activity, climate fluctuations, polyploidization and hybridization together promoted the diversification of this lineage.

Holocene deglaciation drove rapid genetic diversification of Atlantic walrus

Rapid global warming is severely impacting Arctic ecosystems and is predicted to transform the abundance, distribution and genetic diversity of Arctic species, though these linkages are poorly understood. We address this gap in knowledge using palaeogenomics to examine how earlier periods of global warming influenced the genetic diversity of Atlantic walrus (Odobenus rosmarus rosmarus), a species closely associated with sea ice and shallow-water habitats. We analysed 82 ancient and historical Atlantic walrus mitochondrial genomes (mitogenomes), including now-extinct populations in Iceland and the Canadian Maritimes, to reconstruct the Atlantic walrus' response to Arctic deglaciation. Our results demonstrate that the phylogeography and genetic diversity of Atlantic walrus populations was initially shaped by the last glacial maximum (LGM), surviving in distinct glacial refugia, and subsequently expanding rapidly in multiple migration waves during the late Pleistocene and early Holocene. The timing of diversification and establishment of distinct populations corresponds closely with the chronology of the glacial retreat, pointing to a strong link between walrus phylogeography and sea ice. Our results indicate that accelerated ice loss in the modern Arctic may trigger further dispersal events, likely increasing the connectivity of northern stocks while isolating more southerly stocks putatively caught in small pockets of suitable habitat.

The genomic footprint of whaling and isolation in fin whale populations

Twentieth century industrial whaling pushed several species to the brink of extinction, with fin whales being the most impacted. However, a small, resident population in the Gulf of California was not targeted by whaling. Here, we analyzed 50 whole-genomes from the Eastern North Pacific (ENP) and Gulf of California (GOC) fin whale populations to investigate their demographic history and the genomic effects of natural and human-induced bottlenecks. We show that the two populations diverged ~16,000 years ago, after which the ENP population expanded and then suffered a 99% reduction in effective size during the whaling period. In contrast, the GOC population remained small and isolated, receiving less than one migrant per generation. However, this low level of migration has been crucial for maintaining its viability. Our study exposes the severity of whaling, emphasizes the importance of migration, and demonstrates the use of genome-based analyses and simulations to inform conservation strategies.

Temporal matches between monarch butterfly and milkweed population changes over the past 25,000 years

In intimate ecological interactions, the interdependency of species may result in correlated demographic histories. For species of conservation concern, understanding the long-term dynamics of such interactions may shed light on the drivers of population decline. Here, we address the demographic history of the monarch butterfly, Danaus plexippus, and its dominant host plant, the common milkweed Asclepias syriaca (A. syriaca), using broad-scale sampling and genomic inference. Because genetic resources for milkweed have lagged behind those for monarchs, we first release a chromosome-level genome assembly and annotation for common milkweed. Next, we show that despite its enormous geographic range across eastern North America, A. syriaca is best characterized as a single, roughly panmictic population. Using approximate Bayesian computation with random forests (ABC-RF), a machine learning method for reconstructing demographic histories, we show that both monarchs and milkweed experienced population expansion during the most recent recession of North American glaciers 10,000-20,000 years ago. Our data also identify concurrent population expansions in both species during the large-scale clearing of eastern forests (∼200 years ago). Finally, we find no evidence that either species experienced a reduction in effective population size over the past 75 years. Thus, the well-documented decline of monarch abundance over the past 40 years is not visible in our genomic dataset, reflecting a possible mismatch of the overwintering census population to effective population size in this species.

Genomes of cultivated and wild Capsicum species provide insights into pepper domestication and population differentiation

Pepper (Capsicum spp.) is one of the earliest cultivated crops and includes five domesticated species, C. annuum var. annuum, C. chinense, C. frutescens, C. baccatum var. pendulum and C. pubescens. Here, we report a pepper graph pan-genome and a genome variation map of 500 accessions from the five domesticated Capsicum species and close wild relatives. We identify highly differentiated genomic regions among the domesticated peppers that underlie their natural variations in flowering time, characteristic flavors, and unique resistances to biotic and abiotic stresses. Domestication sweeps detected in C. annuum var. annuum and C. baccatum var. pendulum are mostly different, and the common domestication traits, including fruit size, shape and pungency, are achieved mainly through the selection of distinct genomic regions between these two cultivated species. Introgressions from C. baccatum into C. chinense and C. frutescens are detected, including those providing genetic sources for various biotic and abiotic stress tolerances.

Distributional responses to climate change of two maple species in southern China

Climate change is a major factor affecting biodiversity and species distribution, particularly of montane species. Species may respond to climate change by shifting their range to higher elevations. The southeastern Qinghai-Tibetan Plateau (QTP) and the Hengduan Mountains are considered as global biodiversity hotspots. However, information on the response of maple species to climate change in these regions was limited. Therefore, we selected two maple species that occur there and assessed changes in their habitat suitability under past, present and future climate scenarios in Biomod2. The results showed that temperature seasonality (bio4) was the most critical factor influencing their potential distributions. The distribution of potentially suitable habitat for Acer caesium and Acer stachyophyllum was predicted to be larger during the LGM compared to the present. Under the current climate scenario, the largest areas of potentially suitable habitat for these species were mainly located in southeastern Tibet, the Hengduan Mountains in northwestern Yunnan and western Sichuan, the Qinling-Daba Mountains in southern Gansu and the Wumeng-Daliang Mountains in northeastern Yunnan, western Guizhou and southeastern Sichuan. Under future climate change scenarios, the predicted loss of suitable habitat areas for these two species ranged from 13.78% to 45.71% and the increase ranged from 18.88% to 57.98%, with an overall increasing trend. The suitable habitat areas were predicted to shift towards the eastern parts of the QTP under both the pessimistic and optimistic future climate change scenarios in the 2050s and the 2070s, which became evident as global warming intensified, particularly in the eastern QTP and the Hengduan Mountains. Our results highlight the possibility that the diverse topography along altitudinal gradients in the QTP and the Hengduan Mountains may potentially mitigate the range contraction of mountain plants in response to climate warming. These findings provide a basis for planning conservation areas, planting and species conservation in the mountainous areas of southern China under the anticipated global warming.

Genetic diversity and demographic history of the leopard seal: A Southern Ocean top predator

Leopard seals (Hydrurga leptonyx) are top predators that can exert substantial top-down control of their Antarctic prey species. However, population trends and genetic diversity of leopard seals remain understudied, limiting our understanding of their ecological role. We investigated the genetic diversity, effective population size and demographic history of leopard seals to provide fundamental data that contextualizes their predatory influence on Antarctic ecosystems. Ninety leopard seals were sampled from the northern Antarctic Peninsula during the austral summers of 2008-2019 and a 405bp segment of the mitochondrial control region was sequenced for each individual. We uncovered moderate levels of nucleotide (π = 0.013) and haplotype (Hd = 0.96) diversity, and the effective population size was estimated at around 24,000 individuals (NE = 24,376; 95% CI: 16,876-33,126). Consistent with findings from other ice-breeding pinnipeds, Bayesian skyline analysis also revealed evidence for population expansion during the last glacial maximum, suggesting that historical population growth may have been boosted by an increase in the abundance of sea ice. Although leopard seals can be found in warmer, sub-Antarctic locations, the species' core habitat is centered on the Antarctic, making it inherently vulnerable to the loss of sea ice habitat due to climate change. Therefore, detailed assessments of past and present leopard seal population trends are needed to inform policies for Antarctic ecosystems.

Phylogeography of the desert scorpion illuminates a route out of Central Asia

A comprehensive understanding of phylogeography requires the integration of knowledge across different organisms, ecosystems, and geographic regions. However, a critical knowledge gap exists in the arid biota of the vast Asian drylands. To narrow this gap, here we test an "out-of-Central Asia" hypothesis for the desert scorpion Mesobuthus mongolicus by combining Bayesian phylogeographic reconstruction and ecological niche modeling. Phylogenetic analyses of one mitochondrial and three nuclear loci and molecular dating revealed that M. mongolicus represents a coherent lineage that diverged from its most closely related lineage in Central Asia about 1.36 Ma and underwent radiation ever since. Bayesian phylogeographic reconstruction indicated that the ancestral population dispersed from Central Asia gradually eastward to the Gobi region via the Junggar Basin, suggesting that the Junggar Basin has severed as a corridor for Quaternary faunal exchange between Central Asia and East Asia. Two major dispersal events occurred probably during interglacial periods (around 0.8 and 0.4 Ma, respectively) when climatic conditions were analogous to present-day status, under which the scorpion achieved its maximum distributional range. M. mongolicus underwent demographic expansion during the Last Glacial Maximum, although the predicted distributional areas were smaller than those at present and during the Last Interglacial. Development of desert ecosystems in northwest China incurred by intensified aridification might have opened up empty habitats that sustained population expansion. Our results extend the spatiotemporal dimensions of trans-Eurasia faunal exchange and suggest that species' adaptation is an important determinant of their phylogeographic and demographic responses to climate changes.

Genomic Resequencing Unravels the Genetic Basis of Domestication, Expansion, and Trait Improvement in Morus Atropurpurea

Mulberry is an economically important plant in the sericulture industry and traditional medicine. However, the genetic and evolutionary history of mulberry remains largely unknown. Here, this work presents the chromosome-level genome assembly of Morus atropurpurea (M. atropurpurea), originating from south China. Population genomic analysis using 425 mulberry accessions reveal that cultivated mulberry is classified into two species, M. atropurpurea and M. alba, which may have originated from two different mulberry progenitors and have independent and parallel domestication in north and south China, respectively. Extensive gene flow is revealed between different mulberry populations, contributing to genetic diversity in modern hybrid cultivars. This work also identifies the genetic architecture of the flowering time and leaf size. In addition, the genomic structure and evolution of sex-determining regions are identified. This study significantly advances the understanding of the genetic basis and domestication history of mulberry in the north and south, and provides valuable molecular markers of desirable traits for mulberry breeding.

Ocean current patterns drive the worldwide colonization of eelgrass (Zostera marina)

Currents are unique drivers of oceanic phylogeography and thus determine the distribution of marine coastal species, along with past glaciations and sea-level changes. Here we reconstruct the worldwide colonization history of eelgrass (Zostera marina L.), the most widely distributed marine flowering plant or seagrass from its origin in the Northwest Pacific, based on nuclear and chloroplast genomes. We identified two divergent Pacific clades with evidence for admixture along the East Pacific coast. Two west-to-east (trans-Pacific) colonization events support the key role of the North Pacific Current. Time-calibrated nuclear and chloroplast phylogenies yielded concordant estimates of the arrival of Z. marina in the Atlantic through the Canadian Arctic, suggesting that eelgrass-based ecosystems, hotspots of biodiversity and carbon sequestration, have only been present there for ~243 ky (thousand years). Mediterranean populations were founded ~44 kya, while extant distributions along western and eastern Atlantic shores were founded at the end of the Last Glacial Maximum (~19 kya), with at least one major refuge being the North Carolina region. The recent colonization and five- to sevenfold lower genomic diversity of the Atlantic compared to the Pacific populations raises concern and opportunity about how Atlantic eelgrass might respond to rapidly warming coastal oceans.

Mitogenomes revealed the history of bison colonization of Northern Plains after the Last Glacial Maximum

American bison demonstrated differential patterns of extinction, survival, and expansion since the terminal Pleistocene. We determined population dynamics of the Northern Great Plains bison using 40 mitochondrial genomes from radiocarbon dated remains with the age ranging from 12,226 to 167 calibrated years before present. Population dynamics correlated with environmental and anthropogenic factors and was characterized by three primary periods: terminal Pleistocene population growth starting 14,000 years ago, mid Holocene demographic stability between 6700 and 2700 years ago, and late Holocene population decline in the last 2700 years. Most diversification of mtDNA haplotypes occurred in the early Holocene when bison colonized new territories opened by retreating ice sheets. Holocene mtDNA lineages were not found in modern bison and lacked association with archaeological sites and morphological forms.

Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

Adaptations to a cold climate promoted social evolution in Asian colobine primates

The biological mechanisms that underpin primate social evolution remain poorly understood. Asian colobines display a range of social organizations, which makes them good models for investigating social evolution. By integrating ecological, geological, fossil, behavioral, and genomic analyses, we found that colobine primates that inhabit colder environments tend to live in larger, more complex groups. Specifically, glacial periods during the past 6 million years promoted the selection of genes involved in cold-related energy metabolism and neurohormonal regulation. More-efficient dopamine and oxytocin pathways developed in odd-nosed monkeys, which may have favored the prolongation of maternal care and lactation, increasing infant survival in cold environments. These adaptive changes appear to have strengthened interindividual affiliation, increased male-male tolerance, and facilitated the stepwise aggregation from independent one-male groups to large multilevel societies.

Different waves of postglacial recolonisation and genomic structure of bank vole populations in NE Poland

Previous studies indicated that in some species phylogeographic patterns obtained in the analysis of nuclear and mitochondrial DNA (mtDNA) markers can be different. Such mitonuclear discordance can have important evolutionary and ecological consequences. In the present study, we aimed to check whether there was any discordance between mtDNA and nuclear DNA in the bank vole population in the contact zone of its two mtDNA lineages. We analysed the population genetic structure of bank voles using genome-wide genetic data (SNPs) and diversity of sequenced heart transcriptomes obtained from selected individuals from three populations inhabiting areas outside the contact zone. The SNP genetic structure of the populations confirmed the presence of at least two genetic clusters, and such division was concordant with the patterns obtained in the analysis of other genetic markers and functional genes. However, genome-wide SNP analyses revealed the more detailed structure of the studied population, consistent with more than two bank vole recolonisation waves, as recognised previously in the study area. We did not find any significant differences between individuals representing two separate mtDNA lineages of the species in functional genes coding for protein-forming complexes, which are involved in the process of cell respiration in mitochondria. We concluded that the contemporary genetic structure of the populations and the width of the contact zone were shaped by climatic and environmental factors rather than by genetic barriers. The studied populations were likely isolated in separate Last Glacial Maximum refugia for insufficient amount of time to develop significant genetic differentiation.

Ancient bears provide insights into Pleistocene ice age refugia in Southeast Alaska

During the Late Pleistocene, major parts of North America were periodically covered by ice sheets. However, there are still questions about whether ice-free refugia were present in the Alexander Archipelago along the Southeast (SE) Alaska coast during the last glacial maximum (LGM). Numerous subfossils have been recovered from caves in SE Alaska, including American black (Ursus americanus) and brown (U. arctos) bears, which today are found in the Alexander Archipelago but are genetically distinct from mainland bear populations. Hence, these bear species offer an ideal system to investigate long-term occupation, potential refugial survival and lineage turnover. Here, we present genetic analyses based on 99 new complete mitochondrial genomes from ancient and modern brown and black bears spanning the last ~45,000 years. Black bears form two SE Alaskan subclades, one preglacial and another postglacial, that diverged >100,000 years ago. All postglacial ancient brown bears are closely related to modern brown bears in the archipelago, while a single preglacial brown bear is found in a distantly related clade. A hiatus in the bear subfossil record around the LGM and the deep split of their pre- and postglacial subclades fail to support a hypothesis of continuous occupancy in SE Alaska throughout the LGM for either species. Our results are consistent with an absence of refugia along the SE Alaska coast, but indicate that vegetation quickly expanded after deglaciation, allowing bears to recolonize the area after a short-lived LGM peak.