Application of large-scale parentage analysis for investigating natal dispersal in highly vagile vertebrates: a case study of American black bears (Ursus americanus)

Anthropogenic impacts at the interface of animal spatial and social behaviour

Human disturbance is contributing to widespread, global changes in the distributions and densities of wild animals. These anthropogenic impacts on wildlife arise from multiple bottom-up and top-down pathways, including habitat loss, resource provisioning, climate change, pollution, infrastructure development, hunting and our direct presence. Animal behaviour is an important mechanism linking these disturbances to population outcomes, although these behavioural pathways are often complex and can remain obscured when different aspects of behaviour are studied in isolation from one another. The spatial-social interface provides a lens for understanding how an animal's spatial and social environments interact to determine its spatial and social phenotype (i.e. measurable characteristics of an individual), and how these phenotypes interact and feed back to reshape environments. Here, we review studies of animal behaviour at the spatial-social interface to understand and predict how human disturbance affects animal movement, distribution and intraspecific interactions, with consequences for the conservation of populations and ecosystems. By understanding the spatial-social mechanisms linking human disturbance to conservation outcomes, we can better design management interventions to mitigate undesired consequences of disturbance.This article is part of the theme issue 'The spatial-social interface: a theoretical and empirical integration'.

Evidence of population genetic structure in Ecuadorian Andean bears

Wildlife conservation in Andean countries is a global priority because of the high levels of biodiversity and endemism. Historically, these countries have had limited resources to monitor wildlife (e.g., through genetic tools) and establish conservation programs. Focusing on the study and emblematic use of a few charismatic species has been a strategic approach to direct efforts for conservation and development planning. Consequently, the Andean bear is a flagship and umbrella species for highly biodiverse Andean countries like Ecuador. The few studies exploring the population genetics of this species have concluded that it has low genetic diversity and few units for conservation as populations appear to be well connected. However, these results might be attributed to ascertainment bias as studies have been performed with heterologous molecular markers. Here, using both mtDNA sequences and species-specific microsatellite markers, we show that Andean bears in Ecuador have population structure. Additionally, we found through the study of three Ecuadorian populations that the species might have a higher genetic diversity than we previously thought. These results could support the revision of research priorities, conservation, and planning strategies to improve connectivity for this species which occurs in crucial biodiversity hotspots.

Evaluating density-weighted connectivity of black bears (Ursus americanus) in Glacier National Park with spatial capture-recapture models

BACKGROUND

Improved understanding of wildlife population connectivity among protected area networks can support effective planning for the persistence of wildlife populations in the face of land use and climate change. Common approaches to estimating connectivity often rely on small samples of individuals without considering the spatial structure of populations, leading to limited understanding of how individual movement links to demography and population connectivity. Recently developed spatial capture-recapture (SCR) models provide a framework to formally connect inference about individual movement, connectivity, and population density, but few studies have applied this approach to empirical data to support connectivity planning.

METHODS

We used mark-recapture data collected from 924 genetic detections of 598 American black bears (Ursus americanus) in 2004 with SCR ecological distance models to simultaneously estimate density, landscape resistance to movement, and population connectivity in Glacier National Park northwest Montana, USA. We estimated density and movement parameters separately for males and females and used model estimates to calculate predicted density-weighted connectivity surfaces.

RESULTS

Model results indicated that landscape structure influences black bear density and space use in Glacier. The mean density estimate was 16.08 bears/100 km2 (95% CI 12.52-20.6) for females and 9.27 bears/100 km2 (95% CI 7.70-11.14) for males. Density increased with forest cover for both sexes. For male black bears, density decreased at higher grizzly bear (Ursus arctos) densities. Drainages, valley bottoms, and riparian vegetation decreased estimates of landscape resistance to movement for male and female bears. For males, forest cover also decreased estimated resistance to movement, but a transportation corridor bisecting the study area strongly increased resistance to movement presenting a barrier to connectivity.

CONCLUSIONS

Density-weighed connectivity surfaces highlighted areas important for population connectivity that were distinct from areas with high potential connectivity. For black bears in Glacier and surrounding landscapes, consideration of both vegetation and valley topography could inform the placement of underpasses along the transportation corridor in areas characterized by both high population density and potential connectivity. Our study demonstrates that the SCR ecological distance model can provide biologically realistic, spatially explicit predictions to support movement connectivity planning across large landscapes.

Partitioning variance in reproductive success, within years and across lifetimes

Variance in reproductive success (s k 2 , with k = number of offspring) plays a large role in determining the rate of genetic drift and the scope within which selection acts. Various frameworks have been proposed to parse factors that contribute to s k 2 , but none has focused on age-specific values of ϕ = s k 2 / k ¯ , which indicate the degree to which reproductive skew is overdispersed (compared to the random Poisson expectation) among individuals of the same age and sex. Instead, within-age effects are generally lumped with residual variance and treated as "noise." Here, an ANOVA sums-of-squares framework is used to partition variance in annual and lifetime reproductive success into between-group and within-group components. For annual reproduction, the between-age effect depends on age-specific fecundity (b x), but relatively few empirical data are available on the within-age effect, which depends on ϕ x. By defining groups by age-at-death rather than age, the same ANOVA framework can be used to partition variance in lifetime reproductive success (LRS) into between-group and within-group components. Analytical methods are used to develop null-model expectations for random contributions to within-group and between-group components. For analysis of LRS, random variation in longevity appears as part of the between-group variance, and effects (if any) of skip breeding and persistent individual differences contribute to the within-group variance. Simulations are used to show that the methods for variance partitioning are asymptotically unbiased. Practical application is illustrated with empirical data for annual reproduction in American black bears and lifetime reproduction in Dutch great tits. Results show that overdispersed within-age variance (1) dominates annual s k 2 in both male and female black bears, (2) is the primary factor that reduces annual effective size to a fraction of the number of adults, and (3) represents most of the opportunity for selection. In contrast, about a quarter of the variance in LRS in great tits can be attributed to random variation in longevity, and most of the rest is due to modest differences in fecundity with age estimated for a single cohort of females. R code is provided that reads generic input files for annual and lifetime reproductive success and allows users to conduct variance partitioning with their own data.

Characterization of 35 new microsatellite markers for the blacktip reef shark (Carcharhinus melanopterus) and cross-species amplification in eight other shark species

BACKGROUND

Shark species are overfished at a global scale, as they are poached for the finning industry or are caught as bycatch. Efficient conservation measures require fine-scale spatial and temporal studies to characterize shark habitat use, infer migratory habits, analyze relatedness, and detect population genetic differentiation. Gathering these types of data is costly and time-consuming, especially when it requires collection of shark tissue samples.

METHODS AND RESULTS

Genetic tools, such as microsatellite markers, are the most economical sampling method for collecting genetic data, as they enable the estimation of genetic diversity, population structure and parentage relationships and are thus an efficient way to inform conservation strategies. Here, a set of 45 microsatellite loci was tested on three blacktip reef shark (Carcharhinus melanopterus) populations from three Polynesian islands: Moorea, Morane and Tenararo. The set was composed of 10 previously published microsatellite markers and 35 microsatellite markers that were developed specifically for C. melanopterus as part of the present study. The 35 novel and 10 existing loci were cross-amplified on eight additional shark species (Carcharhinus amblyrhynchos, C. longimanus, C. sorrah, Galeocerdo cuvier, Negaprion acutidens, Prionacea glauca, Rhincodon typus and Sphyrna lewini). These species had an average of 69% of successful amplification, considered if at least 50% of the individual samples being successfully amplified per species and per locus.

CONCLUSIONS

This novel microsatellite marker set will help address numerous knowledge gaps that remain, concerning genetic stock identification, shark behavior and reproduction via parentage analysis.

Identification of familial networks reveals sex-specific density dependence in the dispersal and reproductive success of an endangered ungulate

Density is an important demographic parameter that is commonly overlooked in studies of wild populations. Here, we examined the effects of variable spatially explicit density on a range of demographic parameters in a wild population of a cryptic ungulate, boreal woodland caribou (Rangifer tarandus caribou). Using non-invasive genetic sampling, we applied spatial capture–recapture methods with landscape covariates to estimate the density of boreal woodland caribou across a 108,806 km2 study area. We then created a familial network from the reconstructed parent–offspring relationships to determine whether spatial density influenced sex-specific individual reproductive success, female pregnancy status, and dispersal distance. We found that animal density varied greatly in response to land cover types and disturbance; animal density was most influenced by landscape composition and distance to roads varying from 0 in areas with >20% deciduous cover to 270 caribou per 1,000 km2 in areas presenting contiguous older coniferous cover. We found that both male and female reproductive success varied with density, with males showing a higher probability of having offspring in higher-density areas, and the opposite for females. No differences were found in female pregnancy rates occurring in high- and low-density areas. Dispersal distances varied with density, with offspring moving shorter distances when parents were found in higher-density areas. Familial networks showed lower-closeness centrality and lower-degree centrality for females in higher-density areas, indicating that females found in higher-density areas tend to be less broadly associated with animals across the range. Although high-density areas do reflect good-quality caribou habitat, the observed decreased closeness and degree centrality measures, dispersal rates, and lower female recruitment rates suggest that remnant habitat patches across the landscape may create population sinks.

Intact landscape promotes gene flow and low genetic structuring in the threatened Eastern Massasauga Rattlesnake

Genetic structuring of wild populations is dependent on environmental, ecological, and life-history factors. The specific role environmental context plays in genetic structuring is important to conservation practitioners working with rare species across areas with varying degrees of fragmentation. We investigated fine-scale genetic patterns of the federally threatened Eastern Massasauga Rattlesnake (Sistrurus catenatus) on a relatively undisturbed island in northern Michigan, USA. This species often persists in habitat islands throughout much of its distribution due to extensive habitat loss and distance-limited dispersal. We found that the entire island population exhibited weak genetic structuring with spatially segregated variation in effective migration and genetic diversity. The low level of genetic structuring contrasts with previous studies in the southern part of the species' range at comparable fine scales (~7 km), in which much higher levels of structuring were documented. The island population's genetic structuring more closely resembles that of populations from Ontario, Canada, that occupy similarly intact habitats. Intrapopulation variation in effective migration and genetic diversity likely corresponds to the presence of large inland lakes acting as barriers and more human activity in the southern portion of the island. The observed genetic structuring in this intact landscape suggests that the Eastern Massasauga is capable of sufficient interpatch movements to reduce overall genetic structuring and colonize new habitats. Landscape mosaics with multiple habitat patches and localized barriers (e.g., large water bodies or roads) will promote gene flow and natural colonization for this declining species.

Spatial genetic structure and landscape connectivity in black bears: Investigating the significance of using different land cover datasets and classifications in landscape genetics analyses

Landscape genetic analyses allow detection of fine-scale spatial genetic structure (SGS) and quantification of effects of landscape features on gene flow and connectivity. Typically, analyses require generation of resistance surfaces. These surfaces characteristically take the form of a grid with cells that are coded to represent the degree to which landscape or environmental features promote or inhibit animal movement. How accurately resistance surfaces predict association between the landscape and movement is determined in large part by (a) the landscape features used, (b) the resistance values assigned to features, and (c) how accurately resistance surfaces represent landscape permeability. Our objective was to evaluate the performance of resistance surfaces generated using two publicly available land cover datasets that varied in how accurately they represent the actual landscape. We genotyped 365 individuals from a large black bear population (Ursus americanus) in the Northern Lower Peninsula (NLP) of Michigan, USA at 12 microsatellite loci, and evaluated the relationship between gene flow and landscape features using two different land cover datasets. We investigated the relative importance of land cover classification and accuracy on landscape resistance model performance. We detected local spatial genetic structure in Michigan's NLP black bears and found roads and land cover were significantly correlated with genetic distance. We observed similarities in model performance when different land cover datasets were used despite 21% dissimilarity in classification between the two land cover datasets. However, we did find the performance of land cover models to predict genetic distance was dependent on the way the land cover was defined. Models in which land cover was finely defined (i.e., eight land cover classes) outperformed models where land cover was defined more coarsely (i.e., habitat/non-habitat or forest/non-forest). Our results show that landscape genetic researchers should carefully consider how land cover classification changes inference in landscape genetic studies.

Comparative genomics approach to evolutionary process connectivity

The influence of species life history traits and historical demography on contemporary connectivity is still poorly understood. However, these factors partly determine the evolutionary responses of species to anthropogenic landscape alterations. Genetic connectivity and its evolutionary outcomes depend on a variety of spatially dependent evolutionary processes, such as population structure, local adaptation, genetic admixture, and speciation. Over the last years, population genomic studies have been interrogating these processes with increasing resolution, revealing a large diversity of species responses to spatially structured landscapes. In parallel, multispecies meta-analyses usually based on low-genome coverage data have provided fundamental insights into the ecological determinants of genetic connectivity, such as the influence of key life history traits on population structure. However, comparative studies still lack a thorough integration of macro- and micro-evolutionary scales to fully realize their potential. Here, I present how a comparative genomics framework may provide a deeper understanding of evolutionary process connectivity. This framework relies on coupling the inference of long-term demographic and selective history with an assessment of the contemporary consequences of genetic connectivity. Standardizing this approach across several species occupying the same landscape should help understand how spatial environmental heterogeneity has shaped the diversity of historical and contemporary connectivity patterns in different taxa with contrasted life history traits. I will argue that a reasonable amount of genome sequence data can be sufficient to resolve and connect complex macro- and micro-evolutionary histories. Ultimately, implementing this framework in varied taxonomic groups is expected to improve scientific guidelines for conservation and management policies.

Unsustainable anthropogenic mortality disrupts natal dispersal and promotes inbreeding in leopards

Anthropogenic mortality of wildlife is typically inferred from measures of the absolute decline in population numbers. However, increasing evidence suggests that indirect demographic effects including changes to the age, sex, and social structure of populations, as well as the behavior of survivors, can profoundly impact population health and viability. Specifically, anthropogenic mortality of wildlife (especially when unsustainable) and fragmentation of the spatial distribution of individuals (home-ranges) could disrupt natal dispersal mechanisms, with long-term consequences to genetic structure, by compromising outbreeding behavior and gene flow. We investigate this threat in African leopards (Panthera pardus pardus), a polygynous felid with male-biased natal dispersal. Using a combination of spatial (home-range) and genetic (21 polymorphic microsatellites) data from 142 adult leopards, we contrast the structure of two South African populations with markedly different histories of anthropogenically linked mortality. Home-range overlap, parentage assignment, and spatio-genetic autocorrelation together show that historical exploitation of leopards in a recovering protected area has disrupted and reduced subadult male dispersal, thereby facilitating opportunistic male natal philopatry, with sons establishing territories closer to their mothers and sisters. The resultant kin-clustering in males of this historically exploited population is comparable to that of females in a well-protected reserve and has ultimately led to localized inbreeding. Our findings demonstrate novel evidence directly linking unsustainable anthropogenic mortality to inbreeding through disrupted dispersal in a large, solitary felid and expose the genetic consequences underlying this behavioral change. We therefore emphasize the importance of managing and mitigating the effects of unsustainable exploitation on local populations and increasing habitat fragmentation between contiguous protected areas by promoting in situ recovery and providing corridors of suitable habitat that maintain genetic connectivity.

Accounting for Age Structure and Spatial Structure in Eco-Evolutionary Analyses of a Large, Mobile Vertebrate

The idealized concept of a population is integral to ecology, evolutionary biology, and natural resource management. To make analyses tractable, most models adopt simplifying assumptions, which almost inevitably are violated by real species in nature. Here, we focus on both demographic and genetic estimates of effective population size per generation (Ne), the effective number of breeders per year (Nb), and Wright's neighborhood size (NS) for black bears (Ursus americanus) that are continuously distributed in the northern lower peninsula of Michigan, United States. We illustrate practical application of recently developed methods to account for violations of 2 common, simplifying assumptions about populations: 1) reproduction occurs in discrete generations and 2) mating occurs randomly among all individuals. We use a 9-year harvest dataset of >3300 individuals, together with genetic determination of 221 parent-offspring pairs, to estimate male and female vital rates, including age-specific survival, age-specific fecundity, and age-specific variance in fecundity (for which empirical data are rare). We find strong evidence for overdispersed variance in reproductive success of same-age individuals in both sexes, and we show that constraints on litter size have a strong influence on results. We also estimate that another life-history trait that is often ignored (skip breeding by females) has a relatively modest influence, reducing Nb by 9% and increasing Ne by 3%. We conclude that isolation by distance depresses genetic estimates of Nb, which implicitly assume a randomly mating population. Estimated demographic NS (100, based on parent-offspring dispersal) was similar to genetic NS (85, based on regression of genetic distance and geographic distance), indicating that the >36000 km2 study area includes about 4-5 black-bear neighborhoods. Results from this expansive data set provide important insight into effects of violating assumptions when estimating evolutionary parameters for long-lived, free-ranging species. In conjunction with recently developed analytical methodology, the ready availability of nonlethal DNA sampling methods and the ability to rapidly and cheaply survey many thousands of molecular markers should facilitate eco-evolutionary studies like this for many more species in nature.

Testing the influence of habitat experienced during the natal phase on habitat selection later in life in Scandinavian wolves

Natal habitat preference induction (NHPI) occurs when characteristics of the natal habitat influence the future habitat selection of an animal. However, the influence of NHPI after the dispersal phase has received remarkably little attention. We tested whether exposure to humans in the natal habitat helps understand why some adult wolves Canis lupus may approach human settlements more than other conspecifics, a question of both ecological and management interest. We quantified habitat selection patterns within home ranges using resource selection functions and GPS data from 21 wolf pairs in Scandinavia. We identified the natal territory of each wolf with genetic parental assignment, and we used human-related characteristics within the natal territory to estimate the degree of anthropogenic influence in the early life of each wolf. When the female of the adult wolf pair was born in an area with a high degree of anthropogenic influence, the wolf pair tended to select areas further away from humans, compared to wolf pairs from natal territories with a low degree of anthropogenic influence. Yet the pattern was statistically weak, we suggest that our methodological approach can be useful in other systems to better understand NHPI and to inform management  about human-wildlife interactions.

No place like home? A test of the natal habitat-biased dispersal hypothesis in Scandinavian wolves

Natal dispersal is an important mechanism for the viability of populations. The influence of local conditions or experience gained in the natal habitat could improve fitness if dispersing individuals settle in an area with similar habitat characteristics. This process, defined as 'natal habitat-biased dispersal' (NHBD), has been used to explain distribution patterns in large carnivores, but actual studies evaluating it are rare. We tested whether grey wolf Canis lupus territory establishment was influenced by the habitat characteristics of the natal territory using the long-term monitoring of the Scandinavian wolf population. We paired the locations of natal and established territories, accounted for available habitats along the dispersing route, and compared their habitat characteristics for 271 wolves during 1998-2012. Wolves with the shortest dispersal distances established in natal-like habitat types more than expected by chance, whereas wolves that dispersed longer distances did not show NHBD. The pattern was consistent for male and female wolves, with females showing more NHBD than males. Chances to detect NHBD increased with the size of habitat defined as available. This highlights the importance of considering the biological characteristics of the studied species when defining habitat availability. Our methodological approach can prove useful to inform conservation and management to identify habitats to be selected by reintroduced or naturally expanding populations.

Beyond the snapshot: Landscape genetic analysis of time series data reveal responses of American black bears to landscape change

Landscape genetic studies typically focus on the evolutionary processes that give rise to spatial patterns that are quantified at a single point in time. Although landscape change is widely recognized as a strong driver of microevolutionary processes, few landscape genetic studies have directly evaluated the change in spatial genetic structure (SGS) over time with concurrent changes in landscape pattern. We introduce a novel approach to analyze landscape genetic data through time. We demonstrate this approach using genotyped samples (n = 569) from a large black bear (Ursus americanus) population in Michigan (USA) that were harvested during 3 years (2002, 2006, and 2010). We identified areas that were consistently occupied over this 9-year period and quantified temporal variation in SGS. Then, we evaluated alternative hypotheses about effects of changes in landscape features (e.g., deforestation or crop conversion) on fine-scale SGS among years using spatial autoregressive modeling and model selection. Relative measures of landscape change such as magnitude of landscape change (i.e., number of patches changing from suitable to unsuitable states or vice versa), and during later periods, measures of fragmentation (i.e., patch aggregation and cohesion) were associated with change in SGS. Our results stress the importance of conducting time series studies for the conservation and management of wildlife inhabiting rapidly changing landscapes.

Spatial genetic patterns indicate mechanism and consequences of large carnivore cohabitation within development

Patterns of human development are shifting from concentrated housing toward sprawled housing intermixed with natural land cover, and wildlife species increasingly persist in close proximity to housing, roads, and other anthropogenic features. These associations can alter population dynamics and evolutionary trajectories. Large carnivores increasingly occupy urban peripheries, yet the ecological consequences for populations established entirely within urban and exurban landscapes are largely unknown. We applied a spatial and landscape genetics approach, using noninvasively collected genetic data, to identify differences in black bear spatial genetic patterns across a rural‐to‐urban gradient and quantify how development affects spatial genetic processes. We quantified differences in black bear dispersal, spatial genetic structure, and migration between differing levels of development within a population primarily occupying areas with >6 houses/km² in western Connecticut. Increased development disrupted spatial genetic structure, and we found an association between increased housing densities and longer dispersal. We also found evidence that roads limited gene flow among bears in more rural areas, yet had no effect among bears in more developed ones. These results suggest dispersal behavior is condition‐dependent and indicate the potential for landscapes intermixing development and natural land cover to facilitate shifts toward increased dispersal. These changes can affect patterns of range expansion and the phenotypic and genetic composition of surrounding populations. We found evidence that subpopulations occupying more developed landscapes may be sustained by male‐biased immigration, creating potentially detrimental demographic shifts.

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.

Detecting black bear source-sink dynamics using individual-based genetic graphs

Source-sink dynamics affects population connectivity, spatial genetic structure and population viability for many species. We introduce a novel approach that uses individual-based genetic graphs to identify source-sink areas within a continuously distributed population of black bears (Ursus americanus) in the northern lower peninsula (NLP) of Michigan, USA. Black bear harvest samples (n = 569, from 2002, 2006 and 2010) were genotyped at 12 microsatellite loci and locations were compared across years to identify areas of consistent occupancy over time. We compared graph metrics estimated for a genetic model with metrics from 10 ecological models to identify ecological factors that were associated with sources and sinks. We identified 62 source nodes, 16 of which represent important source areas (net flux > 0.7) and 79 sink nodes. Source strength was significantly correlated with bear local harvest density (a proxy for bear density) and habitat suitability. Additionally, resampling simulations showed our approach is robust to potential sampling bias from uneven sample dispersion. Findings demonstrate black bears in the NLP exhibit asymmetric gene flow, and individual-based genetic graphs can characterize source-sink dynamics in continuously distributed species in the absence of discrete habitat patches. Our findings warrant consideration of undetected source-sink dynamics and their implications on harvest management of game species.

Demographic mechanisms underpinning genetic assimilation of remnant groups of a large carnivore

Current range expansions of large terrestrial carnivores are occurring following human-induced range contraction. Contractions are often incomplete, leaving small remnant groups in refugia throughout the former range. Little is known about the underlying ecological and evolutionary processes that influence how remnant groups are affected during range expansion. We used data from a spatially explicit, long-term genetic sampling effort of grizzly bears (Ursus arctos) in the Northern Continental Divide Ecosystem (NCDE), USA, to identify the demographic processes underlying spatial and temporal patterns of genetic diversity. We conducted parentage analysis to evaluate how reproductive success and dispersal contribute to spatio-temporal patterns of genetic diversity in remnant groups of grizzly bears existing in the southwestern (SW), southeastern (SE) and east-central (EC) regions of the NCDE. A few reproductively dominant individuals and local inbreeding caused low genetic diversity in peripheral regions that may have persisted for multiple generations before eroding rapidly (approx. one generation) during population expansion. Our results highlight that individual-level genetic and reproductive dynamics play critical roles during genetic assimilation, and show that spatial patterns of genetic diversity on the leading edge of an expansion may result from historical demographic patterns that are highly ephemeral.