The role of Beringia in human adaptation to Arctic conditions based on results of genomic studies of modern and ancient populations.

The results of studies in Quaternary geology, archeology, paleoanthropology and human genetics demonstrate that the ancestors of Native Americans arrived in mid-latitude North America mainly along the Pacific Northwest Coast, but had previously inhabited the Arctic and during the last glacial maximum were in a refugium in Beringia, a land bridge connecting Eurasia and North America. The gene pool of Native Americans is represented by unique haplogroups of mitochondrial DNA and the Y chromosome, the evolutionary age of which ranges from 13 to 22 thousand years. The results of a paleogenomic analysis also show that during the last glacial maximum Beringia was populated by human groups that had arisen as a result of interaction between the most ancient Upper Paleolithic populations of Northern Eurasia and newcomer groups from East Asia. Approximately 20 thousand years ago the Beringian populations began to form, and the duration of their existence in relative isolation is estimated at about 5 thousand years. Thus, the adaptation of the Beringians to the Arctic conditions could have taken several millennia. The adaptation of Amerindian ancestors to high latitudes and cold climates is supported by genomic data showing that adaptive genetic variants in Native Americans are associated with various metabolic pathways: melanin production processes in the skin, hair and eyes, the functioning of the cardiovascular system, energy metabolism and immune response characteristics. Meanwhile, the analysis of the existing hypotheses about the selection of some genetic variants in the Beringian ancestors of the Amerindians in connection with adaptation to the Arctic conditions (for example, in the FADS, ACTN3, EDAR genes) shows the ambiguity of the testing results, which may be due to the loss of some traces of the "Beringian" adaptation in the gene pools of modern Native Americans. The most optimal strategy for further research seems to be the search for adaptive variant.

Ultraconserved elements improve resolution of marmot phylogeny and offer insights into biogeographic history.

Marmots (Marmota spp.) comprise a lineage of large-bodied ground squirrels that diversified rapidly in the Pleistocene, when the planet quickly transitioned to a drier, colder, and highly seasonal climate-particularly at high latitudes. Fossil evidence indicates the genus spread from North America, across Beringia, and into the European Alps over the course of only a few million years, beginning in the late Pliocene. Marmots are highly adapted to survive long and severely cold winters, and this likely favored their expansion and diversification over this time period. Previous phylogenetic studies have identified two major subgenera of marmots, but the timing of important speciation events and some species relationships have been difficult to resolve. Here we use ultraconserved elements and mitogenomes, with samples from all 15 extant species, to more precisely retrace how and when marmots came to inhabit a vast Holarctic range. Our results indicate marmots arose in North America in the Early Miocene (∼16.3 Ma) and dispersed across the Bering Land Bridge in the Pliocene (∼3-4 Ma); in addition, our fossil-calibrated timeline is suggestive of the rise and spread of open grasslands as being particularly important to marmot diversification. The woodchuck (M. monax) and the Alaska marmot (M. broweri) are found to be more closely related to the Eurasian species than to the other North American species. Paraphyly is evident in the bobak marmot (M. bobak) and the hoary marmot (M. caligata), and in the case of the latter the data are highly suggestive of a second, cryptic species in the Cascade Mountains of Washington.

Late Cenozoic history and the role of Beringia in assembling a Holarctic cestode species complex.

The dynamic climate history that drove sea level fluctuation during past glacial periods mediated the movement of organisms between Asia and North America via the Bering Land Bridge. Investigations of the biogeographic histories of small mammals and their parasites demonstrate facets of a complex history of episodic geographic colonization and refugial isolation that structured diversity across the Holarctic. We use a large multi-locus nuclear DNA sequence dataset to robustly resolve relationships within the cestode genus Arostrilepis (Cyclophyllidea: Hymenolepididae), a widespread parasite of predominantly arvicoline rodents (voles, lemmings). Using this phylogeny, we confirm that several Asian Arostrilepis lineages colonized North America during up to four distinct glacial periods in association with different rodent hosts, consistent with taxon-pulse dynamics. A previously inferred westward dispersal across the land bridge is rejected. We also refine interpretations of past host colonization, providing evidence for several distinct episodes of expanding host range, which probably contributed to diversification by Arostrilepis. Finally, Arostrilepis is shown to be paraphyletic with respect to Hymenandrya thomomyis, a parasite of pocket gophers, confirming that ancient Arostrilepis species colonized new host lineages upon arriving in North America.

Population genomics indicate three different modes of divergence and speciation with gene flow in the green-winged teal duck complex.

The processes leading to divergence and speciation can differ broadly among taxa with different life histories. We examine these processes in a small clade of ducks with historically uncertain relationships and species limits. The green-winged teal (Anas crecca) complex is a Holarctic species of dabbling duck currently categorized as three subspecies (Anas crecca crecca, A. c. nimia, and A. c. carolinensis) with a close relative, the yellow-billed teal (Anas flavirostris) from South America. A. c. crecca and A. c. carolinensis are seasonal migrants, while the other taxa are sedentary. We examined divergence and speciation patterns in this group, determining their phylogenetic relationships and the presence and levels of gene flow among lineages using both mitochondrial and genome-wide nuclear DNA obtained from 1,393 ultraconserved element (UCE) loci. Phylogenetic relationships using nuclear DNA among these taxa showed A. c. crecca, A. c. nimia, and A. c. carolinensis clustering together to form one polytomous clade, with A. flavirostris sister to this clade. This relationship can be summarized as (crecca, nimia, carolinensis)(flavirostris). However, whole mitogenomes revealed a different phylogeny: (crecca, nimia)(carolinensis, flavirostris). The best demographic model for key pairwise comparisons supported divergence with gene flow as the probable speciation mechanism in all three contrasts (crecca-nimia, crecca-carolinensis, and carolinensis-flavirostris). Given prior work, gene flow was expected among the Holarctic taxa, but gene flow between North American carolinensis and South American flavirostris (M âˆ¼0.1-0.4 individuals/generation), albeit low, was not expected. Three geographically oriented modes of divergence are likely involved in the diversification of this complex: heteropatric (crecca-nimia), parapatric (crecca-carolinensis), and (mostly) allopatric (carolinensis-flavirostris). Our study shows that ultraconserved elements are a powerful tool for simultaneously studying systematics and population genomics in systems with historically uncertain relationships and species limits.

Beringia and the peopling of the Western Hemisphere.

Did Beringian environments represent an ecological barrier to humans until less than 15 000 years ago or was access to the Americas controlled by the spatial-temporal distribution of North American ice sheets? Beringian environments varied with respect to climate and biota, especially in the two major areas of exposed continental shelf. The East Siberian Arctic Shelf ('Great Arctic Plain' (GAP)) supported a dry steppe-tundra biome inhabited by a diverse large-mammal community, while the southern Bering-Chukchi Platform ('Bering Land Bridge' (BLB)) supported mesic tundra and probably a lower large-mammal biomass. A human population with west Eurasian roots occupied the GAP before the Last Glacial Maximum (LGM) and may have accessed mid-latitude North America via an interior ice-free corridor. Re-opening of the corridor less than 14 000 years ago indicates that the primary ancestors of living First Peoples, who already had spread widely in the Americas at this time, probably dispersed from the NW Pacific coast. A genetic 'arctic signal' in non-arctic First Peoples suggests that their parent population inhabited the GAP during the LGM, before their split from the former. We infer a shift from GAP terrestrial to a subarctic maritime economy on the southern BLB coast before dispersal in the Americas from the NW Pacific coast.

Investigating the reliability of metapodials as taxonomic Indicators for Beringian horses.

The metapodials of extinct horses have long been regarded as one of the most useful skeletal elements to determine taxonomic identity. However, recent research on both extant and extinct horses has revealed the possibility for plasticity in metapodial morphology, leading to notable variability within taxa. This calls into question the reliability of metapodials in species identification, particularly for species identified from fragmentary remains. Here, we use ten measurements of metapodials from 203 specimens of four Pleistocene horse species from eastern Beringia to test whether there are significant differences in metapodial morphology that support the presence of multiple species. We then reconstruct the body masses for every specimen to assess the range in body size within each species and determine whether species differ significantly from one another in mean body mass. We find that that taxonomic groups are based largely on the overall size of the metapodial, and that all metapodial measurements are highly autocorrelated. We also find that mean body mass differs significantly among most, but not all, species. We suggest that metapodial measurements are unreliable taxonomic indicators for Beringian horses given evidence for plasticity in metapodial morphology and their clear reflection of differences in body mass. We recommend future studies use more reliable indicators of taxonomy to identify Beringian horse species, particularly from localities from which fossils of several species have been recovered. Supplementary Information: The online version contains supplementary material available at 10.1007/s10914-022-09626-4.

Genetic drift drives rapid speciation of an Arctic insular endemic shrew (Sorex pribilofensis).

Episodes of Quaternary environmental change shaped the genomes of extant species, influencing their response to contemporary environments, which are changing rapidly. Island endemics are among the most vulnerable to such change, accounting for a disproportionate number of recent extinctions. To prevent extinctions and conserve island biodiversity it is vital to combine knowledge of species' ecologies with their complex evolutionary histories. The Bering Sea has a history of cyclical island isolation and reconnection, coupled with modern rates of climate change that exceed global averages. The endangered Pribilof Island shrew (Sorex pribilofensis) is endemic to St. Paul Island, Alaska, which was isolated from mainland Beringia ~14,000 years ago by rising sea levels. Using ~11,000 single nucleotide polymorphisms, 17 microsatellites and mitochondrial sequence data, we test predictions about the evolutionary processes driving shrew speciation across Beringia. Our data show considerable differentiation of S. pribilofensis from mainland sibling species, relative to levels of divergence between mainland shrews. We also find a genome-wide loss of diversity and extremely low Ne for S. pribilofensis. We then show that intraspecific genetic diversity is significantly related to interspecific divergence, and that differentiation between S. pribilofensis and other Beringian shrews is highest across loci that are fixed in S. pribilofensis, indicating that strong drift has driven differentiation of this island species. Our findings show that drift as a consequence of Arctic climate cycling can rapidly reshape insular biodiversity. Arctic island species that lack genomic diversity and have evolved in response to past climate may have limited ability to respond to modern environmental changes.

Demography and evolutionary history of grey wolf populations around the Bering Strait.

Glacial and interglacial periods throughout the Pleistocene have been substantial drivers of change in species distributions. Earlier analyses suggested that modern grey wolves (Canis lupus) trace their origin to a single Late Pleistocene Beringian population that expanded east and westwards, starting c. 25,000 years ago (ya). Here, we examined the demographic and phylogeographic histories of extant populations around the Bering Strait with wolves from two inland regions of the Russian Far East (RFE) and one coastal and two inland regions of North-western North America (NNA), genotyped for 91,327 single nucleotide polymorphisms. Our results indicated that RFE and NNA wolves had a common ancestry until c. 34,400 ya, suggesting that these populations started to diverge before the previously proposed expansion out of Beringia. Coastal and inland NNA populations diverged c. 16,000 ya, concordant with the minimum proposed date for the ecological viability of the migration route along the Pacific Northwest coast. Demographic reconstructions for inland RFE and NNA populations reveal spatial and temporal synchrony, with large historical effective population sizes that declined throughout the Pleistocene, possibly reflecting the influence of broadscale climatic changes across continents. In contrast, coastal NNA wolves displayed a consistently lower effective population size than the inland populations. Differences between the demographic history of inland and coastal wolves may have been driven by multiple ecological factors, including historical gene flow patterns, natural landscape fragmentation, and more recent anthropogenic disturbance.

Estimating Movement Rates Between Eurasian and North American Birds That Are Vectors of Avian Influenza.

Avian influenza (AI) is a zoonotic disease that will likely be involved in future pandemics. Because waterbird movements are difficult to quantify, determining the host-specific risk of Eurasian-origin AI movements into North America is challenging. We estimated relative rates of movements, based on long-term evolutionary averages of gene flow, between Eurasian and North American waterbird populations to obtain bidirectional baseline rates of the intercontinental movements of these AI hosts. We used population genomics and coalescent-based demographic models to obtain these gene-flow-based movement estimates. Inferred rates of movement between these continental populations varies greatly among species. Within dabbling ducks, gene flow, relative to effective population size, varies from ∼3 to 24 individuals/generation between Eurasian and American wigeons (Mareca penelope and Mareca americana) to ∼100-300 individuals/generation between continental populations of northern pintails (Anas acuta). These are evolutionary long-term averages and provide a solid foundation for understanding the relative risks of each of these host species in potential intercontinental AI movements. We scale these values to census size for evaluation in that context. In addition to being AI hosts, many of these bird species are also important in the subsistence diets of Alaskans, increasing the risk of direct bird-to-human exposure to Eurasian-origin AI virus. We contrast species-specific rates of intercontinental movements with the importance of each species in Alaskan diets to understand the relative risk of these taxa to humans. Assuming roughly equivalent AI infection rates among ducks, greater scaup (Aythya marila), mallard (Anas platyrhynchos), and northern pintail (Anas acuta) were the top three species presenting the highest risks for intercontinental AI movement both within the natural system and through exposure to subsistence hunters. Improved data on AI infection rates in this region could further refine these relative risk assessments. These directly comparable, species-based intercontinental movement rates and relative risk rankings should help in modeling, monitoring, and mitigating the impacts of intercontinental host and AI movements.

Estimación de las tasas de movimiento entre aves euroasiáticas y norteamericanas que son vectores de la influenza aviar. La influenza aviar es una enfermedad zoonótica que probablemente estará involucrada en futuras pandemias. Debido a que los movimientos de aves acuáticas son difíciles de cuantificar, La determinación del riesgo específico de hospedador de los movimientos de influenza aviar de origen euroasiático en América del Norte es un desafío. Se estimaron las tasas relativas de movimientos, sobre la base de promedios evolutivos a largo plazo del flujo de genes, entre las poblaciones de aves acuáticas euroasiáticas y norteamericanas para obtener tasas de referencia bidireccionales de los movimientos intercontinentales de estos huéspedes de influenza aviar. Se utilizó genómica de poblaciones y modelos demográficos basados en la teoría de la coalescencia para obtener estas estimaciones de movimiento basadas en el flujo de genes. Las tasas inferidas de movimiento entre estas poblaciones continentales varían mucho entre especies. Dentro de los patos chapuceros, el flujo de genes, en relación con el tamaño efectivo de la población, varía aproximadamente de 3 a 24 individuos/generación entre los silbones europeos y americanos (Mareca penelope y Mareca americana) hasta aproximadamente entre 100 a 300 individuos/generación entre poblaciones continentales de ánades rabudos (Anas acuta). Estos son promedios evolutivos a largo plazo y proporcionan una base sólida para comprender los riesgos relativos de cada una de estas especies hospedadoras en posibles movimientos intercontinentales de la influenza aviar. Se evaluaron estos valores al tamaño del censo para evaluarlos en ese contexto. Además de ser huéspedes de influenza aviar, muchas de estas especies de aves también son importantes en las dietas de subsistencia de los habitantes de Alaska, lo que aumenta el riesgo de exposición directa de las aves al ser humano por el virus de la influenza aviar de origen euroasiático. Se contrastaron las tasas específicas de especies de movimientos intercontinentales con la importancia de cada especie en las dietas de personas en Alaska para comprender el riesgo relativo de estos taxones para los humanos. Suponiendo tasas de infección por influenza aviar aproximadamente equivalentes entre patos, el porrón bastardo o pato boludo mayor (Aythya marila), el ánade real (Anas platyrhynchos) y el ánade rabudo eran las tres especies principales que presentaban los mayores riesgos para el movimiento de influenza aviar intercontinental tanto dentro del sistema natural como a través de la exposición a cazadores de subsistencia. La mejora de los datos sobre las tasas de infección por influenza aviar en esta región podría mejorar aún más estas evaluaciones de riesgo relativo. Estas tasas de movimiento intercontinental directamente comparables, basadas en especies, y clasificaciones de riesgo relativo deberían ayudar a modelar, monitorear y mitigar los impactos de los movimientos intercontinentales de huéspedes y de la influenza aviar.

SNP-based phylogenomic inference in Holarctic ground squirrels (Urocitellus).

Resolution of rapid evolutionary radiatons requires harvesting maximal signal from phylogenomic datasets. However, studies of non-model clades often target conserved loci that are characterized by reduced information content, which can negatively affect gene tree precision and species tree accuracy. Single nucleotide polymorphism (SNP)-based methods are an underutilized but potentially valuable tool for estimating phylogeny and divergence times because they do not rely on resolved gene trees, allowing information from many or all variant loci to be leveraged in species tree reconstruction. We evaluated the utility of SNP-based methods in resolving phylogeny of Holarctic ground squirrels (Urocitellus), a radiation that has been difficult to disentangle, even in prior phylogenomic studies. We inferred phylogeny from a dataset of >3,000 ultraconserved element loci (UCEs) using two methods (SNAPP, SVDquartets) and compared our results with a new mitogenome phylogeny. We also systematically evaluated how phasing of UCEs improves per-locus information content, inference of topology, and other parameters within each of these SNP-based methods. Phasing improved topological resolution and branch length estimation at shallow levels (within species complexes), but less so at deeper levels, likely reflecting true uncertainty due to ancestral polymorphisms segregating in rapidly diverging lineages. We resolved key clades in Urocitellus and present targeted opportunities for future phylogenomic inquiry. Our results also extend the roadmap for use of SNPs to address vertebrate radiations and inform comparative analyses at multiple temporal scales.

Lions and brown bears colonized North America in multiple synchronous waves of dispersal across the Bering Land Bridge.

The Bering Land Bridge connecting North America and Eurasia was periodically exposed and inundated by oscillating sea levels during the Pleistocene glacial cycles. This land connection allowed the intermittent dispersal of animals, including humans, between Western Beringia (far northeast Asia) and Eastern Beringia (northwest North America), changing the faunal community composition of both continents. The Pleistocene glacial cycles also had profound impacts on temperature, precipitation and vegetation, impacting faunal community structure and demography. While these palaeoenvironmental impacts have been studied in many large herbivores from Beringia (e.g., bison, mammoths, horses), the Pleistocene population dynamics of the diverse guild of carnivorans present in the region are less well understood, due to their lower abundances. In this study, we analyse mitochondrial genome data from ancient brown bears (Ursus arctos; n = 103) and lions (Panthera spp.; n = 39), two megafaunal carnivorans that dispersed into North America during the Pleistocene. Our results reveal striking synchronicity in the population dynamics of Beringian lions and brown bears, with multiple waves of dispersal across the Bering Land Bridge coinciding with glacial periods of low sea levels, as well as synchronous local extinctions in Eastern Beringia during Marine Isotope Stage 3. The evolutionary histories of these two taxa underline the crucial biogeographical role of the Bering Land Bridge in the distribution, turnover and maintenance of megafaunal populations in North America.

Molecular phylogeny, historical biogeography and revised classification of andrenine bees (Hymenoptera: Andrenidae).

The mining bee subfamily Andreninae (Hymenoptera: Andrenidae) is a widely distributed and diverse group of ground-nesting solitary bees, including numerous species known to be important pollinators. Most of the species diversity of Andreninae is concentrated in the mainly Holarctic genus Andrena, comprising ca. 1550 described species. The subfamily and especially the genus have remained relatively neglected by recent molecular phylogenetic studies, with current classifications relying largely on morphological characters. We sampled ultraconserved element (UCE) sequences from 235 taxa, including all andrenine genera and 98 out of 104 currently recognized Andrena subgenera. Using 419,858 aligned nucleotide sites from 1009 UCE loci, we present a comprehensive molecular phylogenetic analysis of the subfamily. Our analysis supports the recognition of seven distinct genera in the Andreninae: Alocandrena, Ancylandrena, Andrena, Cubiandrena, Euherbstia, Megandrena, and Orphana. Within the genus Andrena, present-day subgeneric concepts revealed high degrees of paraphyly and polyphyly, due to strong homoplasy of morphological characters, necessitating a thorough, extensive revision of the higher classification of the genus. Based on our findings, we place the subgenus Calcarandrena in synonymy with Andrena (Lepidandrena); Hyperandrena, Nemandrena, Scoliandrena, Tylandrena and Zonandrena with A. (Melandrena); Distandrena, Fumandrena and Proxiandrena with A. (Micrandrena); Carandrena with A. (Notandrena); Agandrena with A. (Plastandrena); Geandrena and Xanthandrena with A. (Ptilandrena); Xiphandrena with A. (Scrapteropsis); and Platygalandrena and Poliandrena with A. (Ulandrena) (new synonymies). We additionally reestablish the groups known as Opandrena and Truncandrena as valid subgenera of Andrena. Our results also show that the MRCA of Andrena + Cubiandrena dispersed from the New World to the Palaearctic probably during the Eocene-early Oligocene, followed by 10-14 Neogene dispersal events from the Palaearctic to the Nearctic and 1-6 Neogene dispersals back into the Palaearctic, all within the genus Andrena.

Late Pleistocene shrub expansion preceded megafauna turnover and extinctions in eastern Beringia.

The collapse of the steppe-tundra biome (mammoth steppe) at the end of the Pleistocene is used as an important example of top-down ecosystem cascades, where human hunting of keystone species led to profound changes in vegetation across high latitudes in the Northern Hemisphere. Alternatively, it is argued that this biome transformation occurred through a bottom-up process, where climate-driven expansion of shrub tundra (Betula, Salix spp.) replaced the steppe-tundra vegetation that grazing megafauna taxa relied on. In eastern Beringia, these differing hypotheses remain largely untested, in part because the precise timing and spatial pattern of Late Pleistocene shrub expansion remains poorly resolved. This uncertainty is caused by chronological ambiguity in many lake sediment records, which typically rely on radiocarbon (14C) dates from bulk sediment or aquatic macrofossils-materials that are known to overestimate the age of sediment layers. Here, we reexamine Late Pleistocene pollen records for which 14C dating of terrestrial macrofossils is available and augment these data with 14C dates from arctic ground-squirrel middens and plant macrofossils. Comparing these paleovegetation data with a database of published 14C dates from megafauna remains, we find the postglacial expansion of shrub tundra preceded the regional extinctions of horse (Equus spp.) and mammoth (Mammuthus primigenius) and began during a period when the frequency of 14C dates indicates large grazers were abundant. These results are not consistent with a model of top-down ecosystem cascades and support the hypothesis that climate-driven habitat loss preceded and contributed to turnover in mammal communities.

Comparative phylogeography uncovers evolutionary past of Holarctic dragonflies.

Here, we investigate the evolutionary history of five northern dragonfly species to evaluate what role the last glaciation period may have played in their current distributions. We look at the population structure and estimate divergence times for populations of the following species: Aeshna juncea (Linnaeus), Aeshna subarctica Walker, Sympetrum danae (Sulzer), Libellula quadrimaculata Linnaeus and Somatochlora sahlbergi Trybom across their Holarctic range. Our results suggest a common phylogeographic pattern across all species except for S. sahlbergi. First, we find that North American and European populations are genetically distinct and have perhaps been separated for more than 400,000 years. Second, our data suggests that, based on genetics, populations from the Greater Beringian region (Beringia, Japan and China) have haplotypes that cluster with North America or Europe depending on the species rather than having a shared geographic affinity. This is perhaps a result of fluctuating sea levels and ice sheet coverage during the Quaternary period that influenced dispersal routes and refugia. Indeed, glacial Beringia may have been as much a transit zone as a refugia for dragonflies. Somatochlora sahlbergi shows no genetic variation across its range and therefore does not share the geographic patterns found in the other circumboreal dragonflies studied here. Lastly, we discuss the taxonomic status of Sympetrum danae, which our results indicate is a species complex comprising two species, one found in Eurasia through Beringia, and the other in North America east and south of Beringia. Through this study we present a shared history among different species from different families of dragonflies, which are influenced by the climatic fluctuations of the past.

Divergence, gene flow, and speciation in eight lineages of trans-Beringian birds.

Determining how genetic diversity is structured between populations that span the divergence continuum from populations to biological species is key to understanding the generation and maintenance of biodiversity. We investigated genetic divergence and gene flow in eight lineages of birds with a trans-Beringian distribution, where Asian and North American populations have likely been split and reunited through multiple Pleistocene glacial cycles. Our study transects the speciation process, including eight pairwise comparisons in three orders (ducks, shorebirds and passerines) at population, subspecies and species levels. Using ultraconserved elements (UCEs), we found that these lineages represent conditions from slightly differentiated populations to full biological species. Although allopatric speciation is considered the predominant mode of divergence in birds, all of our best divergence models included gene flow, supporting speciation with gene flow as the predominant mode in Beringia. In our eight lineages, three were best described by a split-migration model (divergence with gene flow), three best fit a secondary contact scenario (isolation followed by gene flow), and two showed support for both models. The lineages were not evenly distributed across a divergence space defined by gene flow (M) and differentiation (FST ), instead forming two discontinuous groups: one with relatively shallow divergence, no fixed single nucleotide polymorphisms (SNPs), and high rates of gene flow between populations; and the second with relatively deeply divergent lineages, multiple fixed SNPs, and low gene flow. Our results highlight the important role that gene flow plays in avian divergence in Beringia.

Seasonal trends in underwater ambient noise near St. Lawrence Island and the Bering Strait.

We measured spatial and temporal patterns of ambient noise in dynamic, relatively pristine Arctic marine habitats and evaluate the contributions of environmental and human noise sources. Long-term acoustic recorders were deployed around St. Lawrence Island and the Bering Strait region within key feeding and migratory corridors for protected species that are inherently important to Native Alaskan cultures. Over 3000 h of data from 14 recorders at nine sites were obtained from October 2014 to June 2017. Spatial and temporal ambient noise patterns were quantified with percentile statistics in 1/3rd-octave bands (0.02-8 kHz). Ice presence strongly influenced ambient noise by influencing the physical environment and presence of marine mammals. High variability in noise was observed within and between sites, largely as a function of ice presence and associated factors. Acute contributions of biological and anthropogenic sources to local ambient noise are compared to monthly averages, demonstrating how they influence Arctic soundscapes.

Protura in Arctic Regions, with Description of Mastodonentomon n. gen. (Acerentomidae, Nipponentominae) and a Key to Known Arctic Taxa.

Protura are widespread, but their presence in the Arctic was first noted only ca. 70 years ago and is still little acknowledged. This work compiles taxonomic information on proturans in the Arctic regions and adds unpublished data from Northern Siberia. Currently, this fauna is represented by 23 species in two orders and 14 genera. The large cosmopolitan genus Eosentomon is represented by only four species, whereas Acerentomidae is much more diverse, with 19 species in 13 genera (eight Nipponentominae, five Acerentominae). Most of the Arctic species possess a larger number of setae than species living in temperate regions. Based on several unique characters, a new genus, Mastodonentomon, is erected for Nipponentomon macleani, and the species is re-described with the original description supplemented with new characters, including head chaetotaxy, seta length, and porotaxy. Proturan occurrence in the Arctic is limited to Beringia, but the majority of species have restricted distributions and none have been found in both the American Arctic and Siberia. This implies relict origins and high levels of proturan endemism in the Arctic. This emerging view on biogeographical history is, however, hampered by the limited extent of available data, which highlights the need for considerably greater survey efforts. A key to Arctic proturans is provided to facilitate further studies.

Impact of past climate warming on genomic diversity and demographic history of collared lemmings across the Eurasian Arctic.

The Arctic climate was warmer than today at the last interglacial and the Holocene thermal optimum. To reveal the impact of past climate-warming events on the demographic history of an Arctic specialist, we examined both mitochondrial and nuclear genomic variation in the collared lemming (Dicrostonyx torquatus, Pallas), a keystone species in tundra communities, across its entire distribution in northern Eurasia. The ancestral phylogenetic position of the West Beringian group and divergence time estimates support the hypothesis of continental range contraction to a single refugial area located in West Beringia during high-magnitude warming of the last interglacial, followed by westward recolonization of northern Eurasia in the last glacial period. The West Beringian group harbors the highest mitogenome diversity and its inferred demography indicates a constantly large effective population size over the Late Pleistocene to Holocene. This suggests that northward forest expansion during recent warming of the Holocene thermal optimum did not affect the gene pool of the collared lemming in West Beringia but reduced genomic diversity and effective population size in all other regions of the Eurasian Arctic. Demographic inference from genomic diversity was corroborated by species distribution modeling showing reduction in species distribution during past climate warming. These conclusions are supported by recent paleoecological evidence suggesting smaller temperature increases and moderate northward forest advances in the extreme northeast of Eurasia during the Late Pleistocene-to-Holocene warming events. This study emphasizes the importance of West Beringia as a potential refugium for cold-adapted Arctic species under ongoing climate warming.

Phylogeography of higher Diptera in glacial and postglacial grasslands in western North America.

BACKGROUND: Pleistocene glaciations have had an important impact on the species distribution and community composition of the North American biota. Species survived these glacial cycles south of the ice sheets and/or in other refugia, such as Beringia. In this study, we assessed, using mitochondrial DNA from three Diptera species, whether flies currently found in Beringian grasslands (1) survived glaciation as disjunct populations in Beringia and in the southern refugium; (2) dispersed northward postglacially from the southern refugium; or (3) arose by a combination of the two. Samples were collected in grasslands in western Canada: Prairies in Alberta and Manitoba; the Peace River region (Alberta); and the southern Yukon Territory. We sequenced two gene regions (658 bp of cytochrome c oxidase subunit I, 510 bp of cytochrome b) from three species of higher Diptera: one with a continuous distribution across grassland regions, and two with disjunct populations between the regions. We used a Bayesian approach to determine population groupings without a priori assumptions and performed analysis of molecular variance (AMOVA) and exact tests of population differentiation (ETPD) to examine their validity. Molecular dating was used to establish divergence times. RESULTS: Two geographically structured populations were found for all species: a southern Prairie and Peace River population, and a Yukon population. Although AMOVA did not show significant differentiation between populations, ETPD did. Divergence time between Yukon and southern populations predated the Holocene for two species; the species with an ambiguous divergence time had high haplotype diversity, which could suggest survival in a Beringian refugium. CONCLUSIONS: Populations of Diptera in Yukon grasslands could have persisted in steppe habitats in Beringia through Pleistocene glaciations. Current populations in the region appear to be a mix of Beringian relict populations and, to a lesser extent, postglacial dispersal northward from southern prairie grasslands.

Y Chromosome Sequences Reveal a Short Beringian Standstill, Rapid Expansion, and early Population structure of Native American Founders.

The Americas were the last inhabitable continents to be occupied by humans, with a growing multidisciplinary consensus for entry 15-25 thousand years ago (kya) from northeast Asia via the former Beringia land bridge [1-4]. Autosomal DNA analyses have dated the separation of Native American ancestors from the Asian gene pool to 23 kya or later [5, 6] and mtDNA analyses to ∼25 kya [7], followed by isolation ("Beringian Standstill" [8, 9]) for 2.4-9 ky and then a rapid expansion throughout the Americas. Here, we present a calibrated sequence-based analysis of 222 Native American and relevant Eurasian Y chromosomes (24 new) from haplogroups Q and C [10], with four major conclusions. First, we identify three to four independent lineages as autochthonous and likely founders: the major Q-M3 and rarer Q-CTS1780 present throughout the Americas, the very rare C3-MPB373 in South America, and possibly the C3-P39/Z30536 in North America. Second, from the divergence times and Eurasian/American distribution of lineages, we estimate a Beringian Standstill duration of 2.7 ky or 4.6 ky, according to alternative models, and entry south of the ice sheet after 19.5 kya. Third, we describe the star-like expansion of Q-M848 (within Q-M3) starting at 15 kya [11] in the Americas, followed by establishment of substantial spatial structure in South America by 12 kya. Fourth, the deep branches of the Q-CTS1780 lineage present at low frequencies throughout the Americas today [12] may reflect a separate out-of-Beringia dispersal after the melting of the glaciers at the end of the Pleistocene.