Chromosome-level assembly of the gray fox (Urocyon cinereoargenteus) confirms the basal loss of PRDM9 in Canidae

Reference genome assemblies have been created from multiple lineages within the Canidae family; however, despite its phylogenetic relevance as a basal genus within the clade, there is currently no reference genome for the gray fox (Urocyon cinereoargenteus). Here, we present a chromosome-level assembly for the gray fox (U. cinereoargenteus), which represents the most contiguous, non-domestic canid reference genome available to date, with 90% of the genome contained in just 34 scaffolds and a contig N50 and scaffold N50 of 59.4 and 72.9 Megabases, respectively. Repeat analyses identified an increased number of simple repeats relative to other canids. Based on mitochondrial DNA, our Vermont sample clusters with other gray fox samples from the northeastern United States and contains slightly lower levels of heterozygosity than gray foxes on the west coast of California. This new assembly lays the groundwork for future studies to describe past and present population dynamics, including the delineation of evolutionarily significant units of management relevance. Importantly, the phylogenetic position of Urocyon allows us to verify the loss of PRDM9 functionality in the basal canid lineage, confirming that pseudogenization occurred at least 10 million years ago.

Can demographic histories explain long-term isolation and recent pulses of asymmetric gene flow between highly divergent grey fox lineages?

Secondary contact zones between deeply divergent, yet interfertile, lineages provide windows into the speciation process. North American grey foxes (Urocyon cinereoargenteus) are divided into western and eastern lineages that diverged approximately 1 million years ago. These ancient lineages currently hybridize in a relatively narrow zone of contact in the southern Great Plains, a pattern more commonly observed in smaller-bodied taxa, which suggests relatively recent contact after a long period of allopatry. Based on local ancestry inference with whole-genome sequencing (n = 43), we identified two distinct Holocene pulses of admixture. The older pulse (500-3500 YBP) reflected unidirectional gene flow from east to west, whereas the more recent pulse (70-200 YBP) of admixture was bi-directional. Augmented with genotyping-by-sequencing data from 216 additional foxes, demographic analyses indicated that the eastern lineage declined precipitously after divergence, remaining small throughout most of the late Pleistocene, and expanding only during the Holocene. Genetic diversity in the eastern lineage was highest in the southeast and lowest near the contact zone, consistent with a westward expansion. Concordantly, distribution modelling indicated that during their isolation, the most suitable habitat occurred far east of today's contact zone or west of the Great Plains. Thus, long-term isolation was likely caused by the small, distant location of the eastern refugium, with recent contact reflecting a large increase in suitable habitat and corresponding demographic expansion from the eastern refugium. Ultimately, long-term isolation in grey foxes may reflect their specialized bio-climatic niche. This system presents an opportunity for future investigation of potential pre- and post-zygotic isolating mechanisms.

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.

Genetic architecture and evolution of color variation in American black bears

Color variation is a frequent evolutionary substrate for camouflage in small mammals, but the underlying genetics and evolutionary forces that drive color variation in natural populations of large mammals are mostly unexplained. The American black bear, Ursus americanus (U. americanus), exhibits a range of colors including the cinnamon morph, which has a similar color to the brown bear, U. arctos, and is found at high frequency in the American southwest. Reflectance and chemical melanin measurements showed little distinction between U. arctos and cinnamon U. americanus individuals. We used a genome-wide association for hair color as a quantitative trait in 151 U. americanus individuals and identified a single major locus (p < 10-13). Additional genomic and functional studies identified a missense alteration (R153C) in Tyrosinase-related protein 1 (TYRP1) that likely affects binding of the zinc cofactor, impairs protein localization, and results in decreased pigment production. Population genetic analyses and demographic modeling indicated that the R153C variant arose 9.36 kya in a southwestern population where it likely provided a selective advantage, spreading both northwards and eastwards by gene flow. A different TYRP1 allele, R114C, contributes to the characteristic brown color of U. arctos but is not fixed across the range.

Pleistocene-Holocene vicariance, not Anthropocene landscape change, explains the genetic structure of American black bear (Ursus americanus) populations in the American Southwest and northern Mexico

The phylogeography of the American black bear (Ursus americanus) is characterized by isolation into glacial refugia, followed by population expansion and genetic admixture. Anthropogenic activities, including overharvest, habitat loss, and transportation infrastructure, have also influenced their landscape genetic structure. We describe the genetic structure of the American black bear in the American Southwest and northern Mexico and investigate how prehistoric and contemporary forces shaped genetic structure and influenced gene flow. Using a suite of microsatellites and a sample of 550 bears, we identified 14 subpopulations organized hierarchically following the distribution of ecoregions and mountain ranges containing black bear habitat. The pattern of subdivision we observed is more likely a product of postglacial habitat fragmentation during the Pleistocene and Holocene, rather than a consequence of contemporary anthropogenic barriers to movement during the Anthropocene. We used linear mixed‐effects models to quantify the relationship between landscape resistance and genetic distance among individuals, which indicated that both isolation by resistance and geographic distance govern gene flow. Gene flow was highest among subpopulations occupying large tracts of contiguous habitat, was reduced among subpopulations in the Madrean Sky Island Archipelago, where montane habitat exists within a lowland matrix of arid lands, and was essentially nonexistent between two isolated subpopulations. We found significant asymmetric gene flow supporting the hypothesis that bears expanded northward from a Pleistocene refugium located in the American Southwest and northern Mexico and that major highways were not yet affecting gene flow. The potential vulnerability of the species to climate change, transportation infrastructure, and the US–Mexico border wall highlights conservation challenges and opportunities for binational collaboration.

Genomic divergence, local adaptation, and complex demographic history may inform management of a popular sportfish species complex

The Neosho Bass (Micropterus velox), a former subspecies of the keystone top-predator and globally popular Smallmouth Bass (M. dolomieu), is endemic and narrowly restricted to small, clear streams of the Arkansas River Basin in the Central Interior Highlands (CIH) ecoregion, USA. Previous studies have detected some morphological, genetic, and genomic differentiation between the Neosho and Smallmouth Basses; however, the extent of neutral and adaptive divergence and patterns of intraspecific diversity are poorly understood. Furthermore, lineage diversification has likely been impacted by gene flow in some Neosho populations, which may be due to a combination of natural biogeographic processes and anthropogenic introductions. We assessed: (1) lineage divergence, (2) local directional selection (adaptive divergence), and (3) demographic history among Smallmouth Bass populations in the CIH using population genomic analyses of 50,828 single-nucleotide polymorphisms (SNPs) obtained through ddRAD-seq. Neosho and Smallmouth Bass formed monophyletic clades with 100% bootstrap support. We identified two major lineages within each species. We discovered six Neosho Bass populations (two nonadmixed and four admixed) and three nonadmixed Smallmouth Bass populations. We detected 29 SNPs putatively under directional selection in the Neosho range, suggesting populations may be locally adapted. Two populations were admixed via recent asymmetric secondary contact, perhaps after anthropogenic introduction. Two other populations were likely admixed via combinations of ancient and recent processes. These species comprise independently evolving lineages, some having experienced historical and natural admixture. These results may be critical for management of Neosho Bass as a distinct species and may aid in the conservation of other species with complex biogeographic histories.

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.

Phylogeographic analysis delimits three evolutionary significant units of least chipmunks in North America and identifies unique genetic diversity within the imperiled Peñasco population

Although least chipmunks (Neotamias minimus) are a widely distributed North American species of least concern, the southernmost population, N. m. atristriatus (Peñasco least chipmunk), is imperiled and a candidate for federal listing as a subspecies. We conducted a phylogeographic analysis across the N. minimus range to assess genomic differentiation and distinctiveness of the N. m. atristriatus population. Additionally, we leveraged the historical component of sampling to conduct a temporal analysis of N. minimus genetic diversity and also considered climate change effects on range persistence probability by projecting a species distribution model into the IPCC5 RCP 2.6 and 8.5 scenarios. We identified three geographically structured groups (West, North, and South) that were supported by both mitochondrial and nuclear data. N. m. atristriatus grouped within a unique South subclade but were not reciprocally monophyletic from N. m. operarius, and nuclear genome analyses did not separate N. m. atristriatus, N. m. caryi, and N. m. operarius. Thus, while least chipmunks in the Southwest represent an evolutionary significant unit, subspecies distinctions were not supported and listing of the Peñasco population as a Distinct Population Segment of N. m. operarius may be warranted. Our results also support consideration of populations with North and West mitogenomes as two additional evolutionary significant units. We found that N. minimus genetic diversity declined by ~87% over the last century, and our models predicted substantial future habitat contraction, including the loss of the full contemporary ranges of N. m. atristriatus, N. m. arizonensis, and N. m. chuskaensis.

Environmental genomics of Late Pleistocene black bears and giant short-faced bears

Analysis of ancient environmental DNA (eDNA) has revolutionized our ability to describe biological communities in space and time,1-3 by allowing for parallel sequencing of DNA from all trophic levels.4-8 However, because environmental samples contain sparse and fragmented data from multiple individuals, and often contain closely related species,9 the field of ancient eDNA has so far been limited to organellar genomes in its contribution to population and phylogenetic studies.5,6,10,11 This is in contrast to data from fossils12,13 where full-genome studies are routine, despite these being rare and their destruction for sequencing undesirable.14-16 Here, we report the retrieval of three low-coverage (0.03×) environmental genomes from American black bear (Ursus americanus) and a 0.04× environmental genome of the extinct giant short-faced bear (Arctodus simus) from cave sediment samples from northern Mexico dated to 16-14 thousand calibrated years before present (cal kyr BP), which we contextualize with a new high-coverage (26×) and two lower-coverage giant short-faced bear genomes obtained from fossils recovered from Yukon Territory, Canada, which date to ∼22-50 cal kyr BP. We show that the Late Pleistocene black bear population in Mexico is ancestrally related to the present-day Eastern American black bear population, and that the extinct giant short-faced bears present in Mexico were deeply divergent from the earlier Beringian population. Our findings demonstrate the ability to separately analyze genomic-scale DNA sequences of closely related species co-preserved in environmental samples, which brings the use of ancient eDNA into the era of population genomics and phylogenetics.

Mitochondrial Genomes of the United States Distribution of Gray Fox (Urocyon cinereoargenteus) Reveal a Major Phylogeographic Break at the Great Plains Suture Zone

We examined phylogeographic structure in gray fox (Urocyon cinereoargenteus) across the United States to identify the location of secondary contact zone(s) between eastern and western lineages and investigate the possibility of additional cryptic intraspecific divergences. We generated and analyzed complete mitochondrial genome sequence data from 75 samples and partial control region mitochondrial DNA sequences from 378 samples to investigate levels of genetic diversity and structure through population- and individual-based analyses including estimates of divergence (FST and SAMOVA), median joining networks, and phylogenies. We used complete mitochondrial genomes to infer phylogenetic relationships and date divergence times of major lineages of Urocyon in the United States. Despite broad-scale sampling, we did not recover additional major lineages of Urocyon within the United States, but identified a deep east-west split (∼0.8 million years) with secondary contact at the Great Plains Suture Zone and confirmed the Channel Island fox (Urocyon littoralis) is nested within U. cinereoargenteus. Genetic diversity declined at northern latitudes in the eastern United States, a pattern concordant with post-glacial recolonization and range expansion. Beyond the east-west divergence, morphologically-based subspecies did not form monophyletic groups, though unique haplotypes were often geographically limited. Gray foxes in the United States displayed a deep, cryptic divergence suggesting taxonomic revision is needed. Secondary contact at a common phylogeographic break, the Great Plains Suture Zone, where environmental variables show a sharp cline, suggests ongoing evolutionary processes may reinforce this divergence. Follow-up study with nuclear markers should investigate whether hybridization is occurring along the suture zone and characterize contemporary population structure to help identify conservation units. Comparative work on other wide-ranging carnivores in the region should test whether similar evolutionary patterns and processes are occurring.

Unraveling the mystery of the glacier bear: Genetic population structure of black bears (Ursus americanus) within the range of a rare pelage type

Glacier bears are a rare grey color morph of American black bear (Ursus americanus) found only in northern Southeast Alaska and a small portion of western Canada. We examine contemporary genetic population structure of black bears within the geographic extent of glacier bears and explore how this structure relates to pelage color and landscape features of a recently glaciated and highly fragmented landscape. We used existing radiocollar data to quantify black bear home-range size within the geographic range of glacier bears. The mean home-range size of female black bears in the study area was 13 km2 (n = 11), whereas the home range of a single male was 86.9 km2. We genotyped 284 bears using 21 microsatellites extracted from noninvasively collected hair as well as tissue samples from harvested bears. We found ten populations of black bears in the study area, including several new populations not previously identified, divided largely by geographic features such as glaciers and marine fjords. Glacier bears were assigned to four populations found on the north and east side of Lynn Canal and the north and west side of Glacier Bay with a curious absence in the nonglaciated peninsula between. Lack of genetic relatedness and geographic continuity between black bear populations containing glacier bears suggest a possible unsampled population or an association with ice fields. Further investigation is needed to determine the genetic basis and the adaptive and evolutionary significance of the glacier bear color morph to help focus black bear conservation management to maximize and preserve genetic diversity.

Non-Pessimistic Predictions of the Distributions and Suitability of Metasequoia glyptostroboides under Climate Change Using a Random Forest Model

Metasequoia glyptostroboides Hu & W. C. Cheng, which is a remarkable rare relict plant, has gradually been reduced to its current narrow range due to climate change. Understanding the comprehensive distribution of M. glyptostroboides under climate change on a large spatio-temporal scale is of great significance for determining its forest adaptation. In this study, based on 394 occurrence data and 10 bioclimatic variables, the global potential distribution of M. glyptostroboides under eight different climate scenarios (i.e., the past three, the current one, and the next four) from the Quaternary glacial to the future was simulated by a random forest model built with the biomod2 package. The key bioclimatic variables affecting the distribution of M. glyptostroboides are BIO2 (mean diurnal range), BIO1 (annual mean temperature), BIO9 (mean temperature of driest quarter), BIO6 (min temperature of coldest month), and BIO18 (precipitation of warmest quarter). The result indicates that the temperature affects the potential distribution of M. glyptostroboides more than the precipitation. A visualization of the results revealed that the current relatively suitable habitats of M. glyptostroboides are mainly distributed in East Asia and Western Europe, with a total area of approximately 6.857 × 106 km2. With the intensification of global warming in the future, the potential distribution and the suitability of M. glyptostroboides have a relatively non-pessimistic trend. Whether under the mild (RCP4.5) and higher (RCP8.5) emission scenarios, the total area of suitable habitats will be wider than it is now by the 2070s, and the habitat suitability will increase to varying degrees within a wide spatial range. After speculating on the potential distribution of M. glyptostroboides in the past, the glacial refugia of M. glyptostroboides were inferred, and projections regarding the future conditions of these places are expected to be optimistic. In order to better protect the species, the locations of its priority protected areas and key protected areas, mainly in Western Europe and East Asia, were further identified. Our results will provide theoretical reference for the long-term management of M. glyptostroboides, and can be used as background information for the restoration of other endangered species in the future.

A precarious future for distinctive peripheral populations of meadow voles (Microtus pennsylvanicus)

Conservation efforts rely on robust taxonomic assessments that should be based on critical assessment of interspecific boundaries, infraspecific variation, and potentially distinctive peripheral populations. The meadow vole (Microtus pennsylvanicus) is widely distributed across North America, including 28 morphologically defined subspecies and numerous isolated populations. Because some subspecies are of high conservation concern, we examined geographic variation across the range of the species to test existing infraspecific taxonomy in terms of local and regional diversification. We sequenced mitochondrial DNA (mtDNA) from 20 subspecies of M. pennsylvanicus and contextualized infraspecific variation through comparison of pairwise genetic distances derived from an extended data set of 63 species of Microtus. We found strong support for at least three divergent clades within M. pennsylvanicus, with observed intraspecific clade divergence exceeding that between several pairwise comparisons of sister species within Microtus. Six nuclear genes were then sequenced to test the validity of mtDNA structure and to further evaluate the possibility of cryptic, species-level diversity using Bayes factor species delimitation (BFD) analyses. BFD consistently and decisively supported multiple species based on the multilocus approach. We propose that taxonomic revision of the meadow vole is required, with the eastern clade now identified as M. pennsylvanicus (Ord 1815), the western clade as M. drummondii (Audubon and Bachman 1853), and the coastal Florida clade as M. dukecampbelli (Woods, Post, and Kilpatrick 1982). We suggest that such an arrangement would more closely reflect evolutionary history and provide critical context for further examination of distinctive southern peripheral populations that harbor novel evolutionary legacies and adaptive potential.

Phylogeography of moose in western North America

Subspecies designations within temperate species’ ranges often reflect populations that were isolated by past continental glaciation, and glacial vicariance is believed to be a primary mechanism behind the diversification of several subspecies of North American cervids. We used genetics and the fossil record to study the phylogeography of three moose subspecies (Alces alces andersoni, A. a. gigas, and A. a. shirasi) in western North America. We sequenced the complete mitochondrial genome (16,341 base pairs; n = 60 moose) and genotyped 13 nuclear microsatellites (n = 253) to evaluate genetic variation among moose samples. We also reviewed the fossil record for detections of all North American cervids to comparatively assess the evidence for the existence of a southern refugial population of moose corresponding to A. a. shirasi during the last glacial maximum of the Pleistocene. Analysis of mtDNA molecular variance did not support distinct clades of moose corresponding to currently recognized subspecies, and mitogenomic haplotype phylogenies did not consistently distinguish individuals according to subspecies groupings. Analysis of population structure using microsatellite loci showed support for two to five clusters of moose, including the consistent distinction of a southern group of moose within the range of A. a. shirasi. We hypothesize that these microsatellite results reflect recent, not deep, divergence and may be confounded by a significant effect of geographic distance on gene flow across the region. Review of the fossil record showed no evidence of moose south of the Wisconsin ice age glaciers ≥ 15,000 years ago. We encourage the integration of our results with complementary analyses of phenotype data, such as morphometrics, originally used to delineate moose subspecies, for further evaluation of subspecies designations for North American moose.

Post-Pleistocene differentiation in a Central Interior Highlands endemic salamander

AIM

For many endemic species with limited dispersal capacities, the relationship between landscape changes and species distributions is still unclear. We characterized the population structure of the endemic ringed salamander (Ambystoma annulatum) across its distribution in the Central Interior Highlands (CIH) of North America, an area of high species endemism, to infer the ecological and evolutionary history of the species.

METHODS

We sampled 498 individuals across the species distribution and characterized the population genetic structure using nuclear microsatellite and mitochondrial DNA (mtDNA) markers.

RESULTS

Ambystoma annulatum exist in two strongly supported nuclear genetic clusters across the CIH that correspond to a northern cluster that includes the Missouri Ozark populations and a southern cluster that includes the Arkansas and Oklahoma Ozarks and the Ouachita Mountains. Our demographic models estimated that these populations diverged approximately 2,700 years ago. Pairwise estimates of genetic differentiation at microsatellite and mtDNA markers indicated limited contemporary gene flow and suggest that genetic differentiation was primarily influenced by changes in the post-Pleistocene landscape of the CIH.

MAIN CONCLUSIONS

Both the geologic history and post-European settlement history of the CIH have influenced the population genetic structure of A. annulatum. The low mtDNA diversity suggests a retraction into and expansion out of refugial areas in the south-central Ozarks, during temperature fluctuations of the Pleistocene and Holocene epochs. Similarly, the estimated divergence time for the two nuclear clusters corresponds to changes in the post-Pleistocene landscape. More recently, decreased A. annulatum gene flow may be a result of increased habitat fragmentation and alteration post-European settlement.

Identifying Archaeological Bone via Non-Destructive ZooMS and the Materiality of Symbolic Expression: Examples from Iroquoian Bone Points

Today, practical, functional and symbolic choices inform the selection of raw materials for worked objects. In cases where we can discern the origin of worked bone, tooth, ivory and antler objects in the past, we assume that similar choices are being made. However, morphological species identification of worked objects is often impossible due to the loss of identifying characteristics during manufacture. Here, we describe a novel non-destructive ZooMS (Zooarchaeology by Mass Spectrometry) method which was applied to bone points from Pre-Contact St. Lawrence Iroquoian village sites in southern Quebec, Canada. The traditional ZooMS technique requires destructive analysis of a sample, which can be problematic when dealing with artefacts. Here we instead extracted proteins from the plastic bags in which the points had been stored. ZooMS analysis revealed hitherto unexpected species, notably black bear (Ursus americanus) and human (Homo sapiens sapiens), used in point manufacture. These surprising results (confirmed through genomic sequencing) highlight the importance of advancing biomolecular research in artefact studies. Furthermore, they unexpectedly and exceptionally allow us to identify and explore the tangible, material traces of the symbolic relationship between bears and humans, central to past and present Iroquoian cosmology and mythology.

Landscape genetics reveal broad and fine-scale population structure due to landscape features and climate history in the northern leopard frog (Rana pipiens) in North Dakota

Prehistoric climate and landscape features play large roles structuring wildlife populations. The amphibians of the northern Great Plains of North America present an opportunity to investigate how these factors affect colonization, migration, and current population genetic structure. This study used 11 microsatellite loci to genotype 1,230 northern leopard frogs (Rana pipiens) from 41 wetlands (30 samples/wetland) across North Dakota. Genetic structure of the sampled frogs was evaluated using Bayesian and multivariate clustering methods. All analyses produced concordant results, identifying a major east-west split between two R. pipiens population clusters separated by the Missouri River. Substructuring within the two major identified population clusters was also found. Spatial principal component analysis (sPCA) and variance partitioning analysis identified distance, river basins, and the Missouri River as the most important landscape factors differentiating R. pipiens populations across the state. Bayesian reconstruction of coalescence times suggested the major east-west split occurred ~13-18 kya during a period of glacial retreat in the northern Great Plains and substructuring largely occurred ~5-11 kya during a period of extreme drought cycles. A range-wide species distribution model (SDM) for R. pipiens was developed and applied to prehistoric climate conditions during the Last Glacial Maximum (21 kya) and the mid-Holocene (6 kya) from the CCSM4 climate model to identify potential refugia. The SDM indicated potential refugia existed in South Dakota or further south in Nebraska. The ancestral populations of R. pipiens in North Dakota may have inhabited these refugia, but more sampling outside the state is needed to reconstruct the route of colonization. Using microsatellite genotype data, this study determined that colonization from glacial refugia, drought dynamics in the northern Great Plains, and major rivers acting as barriers to gene flow were the defining forces shaping the regional population structure of R. pipiens in North Dakota.

Inferring Continuous and Discrete Population Genetic Structure Across Space

A classic problem in population genetics is the characterization of discrete population structure in the presence of continuous patterns of genetic differentiation. Especially when sampling is discontinuous, the use of clustering or assignment methods may incorrectly ascribe differentiation due to continuous processes (e.g., geographic isolation by distance) to discrete processes, such as geographic, ecological, or reproductive barriers between populations. This reflects a shortcoming of current methods for inferring and visualizing population structure when applied to genetic data deriving from geographically distributed populations. Here, we present a statistical framework for the simultaneous inference of continuous and discrete patterns of population structure. The method estimates ancestry proportions for each sample from a set of two-dimensional population layers, and, within each layer, estimates a rate at which relatedness decays with distance. This thereby explicitly addresses the "clines versus clusters" problem in modeling population genetic variation, and remedies some of the overfitting to which nonspatial models are prone. The method produces useful descriptions of structure in genetic relatedness in situations where separated, geographically distributed populations interact, as after a range expansion or secondary contact. We demonstrate the utility of this approach using simulations and by applying it to empirical datasets of poplars and black bears in North America.

Whole-genome analysis of Mustela erminea finds that pulsed hybridization impacts evolution at high latitudes

At high latitudes, climatic shifts hypothetically initiate recurrent episodes of divergence by isolating populations in glacial refugia—ice-free regions that enable terrestrial species persistence. Upon glacial recession, populations subsequently expand and often come into contact with other independently diverging populations, resulting in gene flow. To understand how recurrent periods of isolation and contact may have impacted evolution at high latitudes, we investigated introgression dynamics in the stoat (Mustela erminea), a Holarctic mammalian carnivore, using whole-genome sequences. We identify two spatio-temporally distinct episodes of introgression coincident with large-scale climatic shifts: contemporary introgression in a mainland contact zone and ancient contact ~200 km south of the contemporary zone, in the archipelagos along North America’s North Pacific Coast. Repeated episodes of gene flow highlight the central role of cyclic climates in structuring high-latitude diversity, through refugial divergence and introgressive hybridization. When introgression is followed by allopatric isolation (e.g., insularization) it may ultimately expedite divergence.

Jocelyn Colella et al. report whole-genome sequences of 10 stoats (Mustela erminea) from four regions of glacial refugia. They find evidence for two past introgressive events between lineages that coincide with interglacial periods, a pattern that may extend to other high–latitude species.

Spatial genetic structure in American black bears (Ursus americanus): female philopatry is variable and related to population history

Previously, American black bears (Ursus americanus) were thought to follow the pattern of female philopatry and male-biased dispersal. However, recent studies have identified deviations from this pattern. Such flexibility in dispersal patterns can allow individuals greater ability to acclimate to changing environments. We explored dispersal and spatial genetic relatedness patterns across ten black bear populations-including long established (historic), with known reproduction >50 years ago, and newly established (recent) populations, with reproduction recorded <50 years ago-in the Interior Highlands and Southern Appalachian Mountains, United States. We used spatially explicit, individual-based genetic simulations to model gene flow under scenarios with varying levels of population density, genetic diversity, and female philopatry. Using measures of genetic distance and spatial autocorrelation, we compared metrics between sexes, between population types (historic and recent), and among simulated scenarios which varied in density, genetic diversity, and sex-biased philopatry. In empirical populations, females in recent populations exhibited stronger patterns of isolation-by-distance (IBD) than females and males in historic populations. In simulated populations, low-density populations had a stronger indication of IBD than medium- to high-density populations; however, this effect varied in empirical populations. Condition-dependent dispersal strategies may permit species to cope with novel conditions and rapidly expand populations. Pattern-process modeling can provide qualitative and quantitative means to explore variable dispersal patterns, and could be employed in other species, particularly to anticipate range shifts in response to changing climate and habitat conditions.

A southern African origin and cryptic structure in the highly mobile plains zebra

The plains zebra (Equus quagga) is an ecologically important species of the African savannah. It is also one of the most numerous and widely distributed ungulates, and six subspecies have been described based on morphological variation. However, the within-species evolutionary processes have been difficult to resolve due to its high mobility and a lack of consensus regarding the population structure. We obtained genome-wide DNA polymorphism data from more than 167,000 loci for 59 plains zebras from across the species range, encompassing all recognized extant subspecies, as well as three mountain zebras (Equus zebra) and three Grevy's zebras (Equus grevyi). Surprisingly, the population genetic structure does not mirror the morphology-based subspecies delineation, underlining the dangers of basing management units exclusively on morphological variation. We use demographic modelling to provide insights into the past phylogeography of the species. The results identify a southern African location as the most likely source region from which all extant populations expanded around 370,000 years ago. We show evidence for inclusion of the extinct and phenotypically divergent quagga (Equus quagga quagga) in the plains zebra variation and reveal that it was less divergent from the other subspecies than the northernmost (Ugandan) extant population.

Genetic health and population monitoring of two small black bear (Ursus americanus) populations in Alabama, with a regional perspective of genetic diversity and exchange

One of the major concerns in conservation today is the loss of genetic diversity which is a frequent consequence of population isolation and small population sizes. Fragmentation of populations and persecution of carnivores has posed a substantial threat to the persistence of free ranging carnivores in North America since the arrival of European settlers. Black bears have seen significant reductions in range size from their historic extent, which is most pronounced in the southeastern United States and even more starkly in Alabama where until recently bears were reduced to a single geographically isolated population in the Mobile River Basin. Recently a second population has naturally re-established itself in northeastern Alabama. We sought to determine size, genetic diversity and genetic connectivity for these two populations in relation to other regional populations. Both populations of black bears in Alabama had small population sizes and had moderate to low genetic diversity, but showed different levels of connectivity to surrounding populations of bears. The Mobile River Basin population had a small population size at only 86 individuals (76–124, 95% C.I.), the lowest genetic diversity of compared populations (richness = 2.33, Ho and He = 0.33), and showed near complete genetic isolation from surrounding populations across multiple tests. The newly recolonizing population in northeastern Alabama had a small but growing population doubling in 3 years (34 individuals 26–43, 95% C.I.), relatively moderate genetic diversity compared to surrounding populations (richness = 3.32, Ho = 0.53, He = 0.65), and showed a high level of genetic connectivity with surrounding populations.

Consequences of severe habitat fragmentation on density, genetics, and spatial capture-recapture analysis of a small bear population

Loss and fragmentation of natural habitats caused by human land uses have subdivided several formerly contiguous large carnivore populations into multiple small and often isolated subpopulations, which can reduce genetic variation and lead to precipitous population declines. Substantial habitat loss and fragmentation from urban development and agriculture expansion relegated the Highlands-Glades subpopulation (HGS) of Florida, USA, black bears (Ursus americanus floridanus) to prolonged isolation; increasing human land development is projected to cause ≥ 50% loss of remaining natural habitats occupied by the HGS in coming decades. We conducted a noninvasive genetic spatial capture-recapture study to quantitatively describe the degree of contemporary habitat fragmentation and investigate the consequences of habitat fragmentation on population density and genetics of the HGS. Remaining natural habitats sustaining the HGS were significantly more fragmented and patchier than those supporting Florida’s largest black bear subpopulation. Genetic diversity was low (A_R = 3.57; H_E = 0.49) and effective population size was small (N_E = 25 bears), both of which remained unchanged over a period spanning one bear generation despite evidence of some immigration. Subpopulation density (0.054 bear/km²) was among the lowest reported for black bears, was significantly female-biased, and corresponded to a subpopulation size of 98 bears in available habitat. Conserving remaining natural habitats in the area occupied by the small, genetically depauperate HGS, possibly through conservation easements and government land acquisition, is likely the most important immediate step to ensuring continued persistence of bears in this area. Our study also provides evidence that preferentially placing detectors (e.g., hair traps or cameras) primarily in quality habitat across fragmented landscapes poses a challenge to estimating density-habitat covariate relationships using spatial capture-recapture models. Because habitat fragmentation and loss are likely to increase in severity globally, further investigation of the influence of habitat fragmentation and detector placement on estimation of this relationship is warranted.

Phylogeography of a widespread small carnivore, the western spotted skunk (Spilogale gracilis) reveals temporally variable signatures of isolation across western North America

We analyzed phylogeographic patterns in the western spotted skunk, Spilogale gracilis Merriam, 1890 (Carnivora: Mephitidae) in relation to historical events associated with Pre-Pleistocene Divergence (PPD) and Quaternary climate change (QCC) using mitochondrial DNA from 97 individuals distributed across Western North America. Divergence times were generated using BEAST to estimate when isolation in putative refugia occurred. Patterns and timing of demographic expansion was performed using Bayesian skyline plot. Putative climatic refugia resulting from Quaternary climate change were identified using paleoecological niche modeling and divergence dates compared to major vicariant events associated with Pre-Pleistocene conditions. We recovered three major mitochondrial clades corresponding to western North America (California, Baja, and across the Great Basin), east-central North America (Texas, central Mexico, New Mexico), and southwestern Arizona/northwestern Mexico. Time to most recent common ancestor for S. gracilis occurred ~1.36 Ma. Divergence times for each major clade occurred between 0.25 and 0.12 Ma, with signature of population expansion occurring 0.15 and 0.10 Ma. Ecological niche models identified three potential climatic refugia during the Last Interglacial, (1) west coast of California and Oregon, (2) northwestern Mexico, and (3) southern Texas/northeastern Mexico as well as two refugia during the Last Glacial Maximum, (1) western USA and (2) southern Texas/northeastern Mexico. This study supports PPD in shaping species-level diversity compared to QCC-driven changes at the intraspecific level for Spilogale, similar to the patterns reported for other small mammals (e.g., rodents and bats). Phylogeographic patterns also appear to have been shaped by both habitat and river vicariance, especially across the desert southwest. Further, continuing climate change during the Holocene coupled with anthropogenic modifications during the Anthropocene appears to be removing both of these barriers to current dispersal of western spotted skunks.

The evolutionary history of bears is characterized by gene flow across species

Bears are iconic mammals with a complex evolutionary history. Natural bear hybrids and studies of few nuclear genes indicate that gene flow among bears may be more common than expected and not limited to polar and brown bears. Here we present a genome analysis of the bear family with representatives of all living species. Phylogenomic analyses of 869 mega base pairs divided into 18,621 genome fragments yielded a well-resolved coalescent species tree despite signals for extensive gene flow across species. However, genome analyses using different statistical methods show that gene flow is not limited to closely related species pairs. Strong ancestral gene flow between the Asiatic black bear and the ancestor to polar, brown and American black bear explains uncertainties in reconstructing the bear phylogeny. Gene flow across the bear clade may be mediated by intermediate species such as the geographically wide-spread brown bears leading to large amounts of phylogenetic conflict. Genome-scale analyses lead to a more complete understanding of complex evolutionary processes. Evidence for extensive inter-specific gene flow, found also in other animal species, necessitates shifting the attention from speciation processes achieving genome-wide reproductive isolation to the selective processes that maintain species divergence in the face of gene flow.