Negative density-dependent dispersal in the American black bear (Ursus americanus) revealed by noninvasive sampling and genotyping

Large-scale genotypic identification reveals density-dependent natal dispersal patterns in an elusive bird of prey

BACKGROUND

Natal dispersal, the distance between site of birth and site of first breeding, has a fundamental role in population dynamics and species' responses to environmental changes. Population density is considered a key driver of natal dispersal. However, few studies have been able to examine densities at both the natal and the settlement site, which is critical for understanding the role of density in dispersal. Additionally, the role of density on natal dispersal remains poorly understood in long-lived and slowly reproducing species, due to their prolonged dispersal periods and often elusive nature. We studied the natal dispersal of the white-tailed eagle (Haliaeetus albicilla) in response to local breeder densities. We investigated the effects of the number of active territories around the natal site on (a) natal dispersal distance and (b) the difference between natal and settlement site breeder density. We were interested in whether eagles showed tendencies of conspecific attraction (positive density-dependence) or intraspecific competition (negative density-dependence) and how this related to settlement site breeder density.

METHODS

We used a combination of long-term visual and genotypic identification to match individuals from their breeding site to their natal nest. We identified natal dispersal events for 355 individuals hatched between 1984 and 2015 in the Baltic Sea coast and Arctic areas of Finland. Of those, 251 were identified by their genotype.

RESULTS

Individuals born in high-density areas dispersed shorter distances than those born in low-density areas, but settled at lower density breeding sites in comparison to their natal site. Eagles born in low natal area densities dispersed farther but settled in higher density breeding sites compared to their natal site.

CONCLUSIONS

We show that eagles might be attracted by conspecifics (positive density-dependence) to identify high-quality habitats or find mates, but do not settle in the most densely populated areas. This indicates that natal dispersal is affected by an interplay of conspecific attraction and intraspecific competition, which has implications for population dynamics of white-tailed eagles, but also other top predators. Furthermore, our study demonstrates the value of long-term collection of both nestling and (non-invasive) adult DNA samples, and thereafter using genotype matching to identify individuals in long-lived and elusive species.

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.

Timing and distance of natal dispersal in Asian black bears

Dispersal has important implications for population ecology and genetics of a species through redistribution of individuals. In most mammals, males leave their natal area before they reach sexual maturity, whereas females are commonly philopatric. Here, we investigate the patterns of natal dispersal in the Asian black bear (Ursus thibetanus) based on data from 550 bears (378 males, 172 females) captured or removed in Gunma and Tochigi prefectures on central Honshu Island, Japan in 2003–2018. We used genetic data and parentage analysis to investigate sex-biased differences in the distance of natal dispersal. We further investigated the age of dispersal using spatial autocorrelation analysis, that is, the change in the correlation between genetic and geographic distances in each sex and age group. Our results revealed that male dispersal distances (mean ± SE = 17.4 ± 3.5 km) were significantly farther than female distances (4.8 ± 1.7 km), and the results were not affected by years of mast failures, a prominent forage source for this population. Based on an average adult female home range radius of 1.8 km, 96% of the males and 50% of the females dispersed. In the spatial autocorrelation analysis, the changes in the relationship between genetic and geographic distances were more pronounced in males compared to females. Males seem to mostly disperse at age 3 regardless of mast productivity, and they gradually disperse far from their home range, but young and inexperienced males may return to their natal home range in years with poor food conditions. The results suggest that factors driving the dispersal process seem to be population structure-based instead of forage availability-based. In females, a significant genetic relationship was observed among all individuals in the group with a minimum age of 6 years within a distance of 2 km, which resulted in the formation of matrilineal assemblages.

Estimation of breeding population size using DNA-based pedigree reconstruction in brown bears

Robust estimates of demographic parameters are critical for effective wildlife conservation and management but are difficult to obtain for elusive species. We estimated the breeding and adult population sizes, as well as the minimum population size, in a high-density brown bear population on the Shiretoko Peninsula, in Hokkaido, Japan, using DNA-based pedigree reconstruction. A total of 1288 individuals, collected in and around the Shiretoko Peninsula between 1998 and 2020, were genotyped at 21 microsatellite loci. Among them, 499 individuals were identified by intensive genetic sampling conducted in two consecutive years (2019 and 2020) mainly by noninvasive methods (e.g., hair and fecal DNA). Among them, both parents were assigned for 330 bears, and either maternity or paternity was assigned to 47 and 76 individuals, respectively. The subsequent pedigree reconstruction indicated a range of breeding and adult (≥4 years old) population sizes: 128-173 for female breeders and 66-91 male breeders, and 155-200 for female adults and 84-109 male adults. The minimum population size was estimated to be 449 (252 females and 197 males) in 2019. Long-term continuous genetic sampling prior to a short-term intensive survey would enable parentage to be identified in a population with a high probability, thus enabling reliable estimates of breeding population size for elusive species.

Explaining detection heterogeneity with finite mixture and non-Euclidean movement in spatially explicit capture-recapture models

Landscape structure affects animal movement. Differences between landscapes may induce heterogeneity in home range size and movement rates among individuals within a population. These types of heterogeneity can cause bias when estimating population size or density and are seldom considered during analyses. Individual heterogeneity, attributable to unknown or unobserved covariates, is often modelled using latent mixture distributions, but these are demanding of data, and abundance estimates are sensitive to the parameters of the mixture distribution. A recent extension of spatially explicit capture-recapture models allows landscape structure to be modelled explicitly by incorporating landscape connectivity using non-Euclidean least-cost paths, improving inference, especially in highly structured (riparian & mountainous) landscapes. Our objective was to investigate whether these novel models could improve inference about black bear (Ursus americanus) density. We fit spatially explicit capture-recapture models with standard and complex structures to black bear data from 51 separate study areas. We found that non-Euclidean models were supported in over half of our study areas. Associated density estimates were higher and less precise than those from simple models and only slightly more precise than those from finite mixture models. Estimates were sensitive to the scale (pixel resolution) at which least-cost paths were calculated, but there was no consistent pattern across covariates or resolutions. Our results indicate that negative bias associated with ignoring heterogeneity is potentially severe. However, the most popular method for dealing with this heterogeneity (finite mixtures) yielded potentially unreliable point estimates of abundance that may not be comparable across surveys, even in data sets with 136-350 total detections, 3-5 detections per individual, 97-283 recaptures, and 80-254 spatial recaptures. In these same study areas with high sample sizes, we expected that landscape features would not severely constrain animal movements and modelling non-Euclidian distance would not consistently improve inference. Our results suggest caution in applying non-Euclidean SCR models when there is no clear landscape covariate that is known to strongly influence the movement of the focal species, and in applying finite mixture models except when abundant data are available.

Testing the precision and sensitivity of density estimates obtained with a camera-trap method revealed limitations and opportunities

The use of camera traps in ecology helps affordably address questions about the distribution and density of cryptic and mobile species. The random encounter model (REM) is a camera-trap method that has been developed to estimate population densities using unmarked individuals. However, few studies have evaluated its reliability in the field, especially considering that this method relies on parameters obtained from collared animals (i.e., average speed, in km/h), which can be difficult to acquire at low cost and effort. Our objectives were to (1) assess the reliability of this camera-trap method and (2) evaluate the influence of parameters coming from different populations on density estimates. We estimated a reference density of black bears (Ursus americanus) in Forillon National Park (Québec, Canada) using a spatial capture-recapture estimator based on hair-snag stations. We calculated average speed using telemetry data acquired from four different bear populations located outside our study area and estimated densities using the REM. The reference density, determined with a Bayesian spatial capture-recapture model, was 2.87 individuals/10km2 [95% CI: 2.41-3.45], which was slightly lower (although not significatively different) than the different densities estimated using REM (ranging from 4.06-5.38 bears/10km2 depending on the average speed value used). Average speed values obtained from different populations had minor impacts on REM estimates when the difference in average speed between populations was low. Bias in speed values for slow-moving species had more influence on REM density estimates than for fast-moving species. We pointed out that a potential overestimation of density occurs when average speed is underestimated, that is, using GPS telemetry locations with large fix-rate intervals. Our study suggests that REM could be an affordable alternative to conventional spatial capture-recapture, but highlights the need for further research to control for potential bias associated with speed values determined using GPS telemetry data.

Harnessing paleo-environmental modeling and genetic data to predict intraspecific genetic structure

Spatially explicit simulations of gene flow within complex landscapes could help forecast the responses of populations to global and anthropological changes. Simulating how past climate change shaped intraspecific genetic variation can provide a validation of models in anticipation of their use to predict future changes. We review simulation models that provide inferences on population genetic structure. Existing simulation models generally integrate complex demographic and genetic processes but are less focused on the landscape dynamics. In contrast to previous approaches integrating detailed demographic and genetic processes and only secondarily landscape dynamics, we present a model based on parsimonious biological mechanisms combining habitat suitability and cellular processes, applicable to complex landscapes. The simulation model takes as input (a) the species dispersal capacities as the main biological parameter, (b) the species habitat suitability, and (c) the landscape structure, modulating dispersal. Our model emphasizes the role of landscape features and their temporal dynamics in generating genetic differentiation among populations within species. We illustrate our model on caribou/reindeer populations sampled across the entire species distribution range in the Northern Hemisphere. We show that simulations over the past 21 kyr predict a population genetic structure that matches empirical data. This approach looking at the impact of historical landscape dynamics on intraspecific structure can be used to forecast population structure under climate change scenarios and evaluate how species range shifts might induce erosion of genetic variation within species.

Past, present and future contributions of evolutionary biology to wildlife forensics, management and conservation

Successfully implementing fundamental concepts into concrete applications is challenging in any given field. It requires communication, collaboration and shared will between researchers and practitioners. We argue that evolutionary biology, through research work linked to conservation, management and forensics, had a significant impact on wildlife agencies and department practices, where new frameworks and applications have been implemented over the last decades. The Quebec government's Wildlife Department (MFFP: Ministère des Forêts, de la Faune et des Parcs) has been proactive in reducing the “research–implementation” gap, thanks to prolific collaborations with many academic researchers. Among these associations, our department's outstanding partnership with Dr. Louis Bernatchez yielded significant contributions to harvest management, stocking programmes, definition of conservation units, recovery of threatened species, management of invasive species and forensic applications. We discuss key evolutionary biology concepts and resulting concrete examples of their successful implementation that derives directly or indirectly from this successful partnership. While old and new threats to wildlife are bringing new challenges, we expect recent developments in eDNA and genomics to provide innovative solutions as long as the research–implementation bridge remains open.

Sex-biased dispersal and inbreeding avoidance in Hokkaido brown bears

Natal dispersal likely plays an important role in avoiding inbreeding among large carnivores. We tested the hypothesis that male-biased dispersal reduces close inbreeding by limiting the spatial overlap of opposite-sex pairs of close relatives in brown bears (Ursus arctos) in the Shiretoko Peninsula, Hokkaido, Japan. We genotyped 837 individuals collected in 1998–2017 at 21 microsatellite loci and performed parentage analysis. To calculate natal dispersal distance, we considered the site where the mother was identified as the birthplace of her offspring, and the site where the offspring were identified as their dispersed place. As predicted, we found that dispersal distances were significantly greater for males (12.4 km ± 1.0) than for females (7.7 km ± 0.9), and those for males increased from 3 years old, indicating that males begin to disperse around the time sexual maturation begins. Relatedness decreased with distance among pairs of females, and the mean relatedness was significantly higher between pairs of females than between pairs of males or between female–male pairs within 3 km. Closely related female–male pairs rarely (5–6%) resided in close proximity (< 3 km), compared with pairs of closely related females. Our study revealed that the potential for close inbreeding was low in Hokkaido brown bears because males are effective dispersers.

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.

Mother-offspring distances reflect sex differences in fine-scale genetic structure of eastern grey kangaroos

Natal dispersal affects life history and population biology and causes gene flow. In mammals, dispersal is usually male-biased so that females tend to be philopatric and surrounded by matrilineal kin, which may lead to preferential associations among female kin. Here we combine genetic analyses and behavioral observations to investigate spatial genetic structure and sex-biased dispersal patterns in a high-density population of mammals showing fission–fusion group dynamics. We studied eastern grey kangaroos (Macropus giganteus) over 2 years at Wilsons Promontory National Park, Australia, and found weak fine-scale genetic structure among adult females in both years but no structure among adult males. Immature male kangaroos moved away from their mothers at 18–25 months of age, while immature females remained near their mothers until older. A higher proportion of male (34%) than female (6%) subadults and young adults were observed to disperse, although median distances of detected dispersals were similar for both sexes. Adult females had overlapping ranges that were far wider than the maximum extent of spatial genetic structure found. Female kangaroos, although weakly philopatric, mostly encounter nonrelatives in fission–fusion groups at high density, and therefore kinship is unlikely to strongly affect sociality.

Fine-scale genetic structure analyses suggest further male than female dispersal in mountain gorillas

BACKGROUND

Molecular studies in social mammals rarely compare the inferences gained from genetic analyses with field information, especially in the context of dispersal. In this study, we used genetic data to elucidate sex-specific dispersal dynamics in the Virunga Massif mountain gorilla population (Gorilla beringei beringei), a primate species characterized by routine male and female dispersal from stable mixed-sex social groups. Specifically, we conducted spatial genetic structure analyses for each sex and linked our genetically-based observations with some key demographic and behavioural data from this population.

RESULTS

To investigate the spatial genetic structure of mountain gorillas, we analysed the genotypes of 193 mature individuals at 11 microsatellite loci by means of isolation-by-distance and spatial autocorrelation analyses. Although not all males and females disperse, female gorillas displayed an isolation-by-distance pattern among groups and a signal of dispersal at short distances from their natal group based on spatial autocorrelation analyses. In contrast, male genotypes were not correlated with spatial distance, thus suggesting a larger mean dispersal distance for males as compared to females. Both within sex and mixed-sex pairs were on average genetically more related within groups than among groups.

CONCLUSIONS

Our study provides evidence for an intersexual difference in dispersal distance in the mountain gorilla. Overall, it stresses the importance of investigating spatial genetic structure patterns on a sex-specific basis to better understand the dispersal dynamics of the species under investigation. It is currently poorly understood why some male and female gorillas disperse while others remain in the natal group. Our results on average relatedness within and across groups confirm that groups often contain close relatives. While inbreeding avoidance may play a role in driving female dispersal, we note that more detailed dyadic genetic analyses are needed to shed light on the role of inbreeding avoidance as an ultimate cause of female dispersal in mountain gorillas.

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

Understanding the factors that affect dispersal is a fundamental question in ecology and conservation biology, particularly as populations are faced with increasing anthropogenic impacts. Here we collected georeferenced genetic samples (n = 2,540) from three generations of black bears (Ursus americanus) harvested in a large (47,739 km²), geographically isolated population and used parentage analysis to identify mother-offspring dyads (n = 337). We quantified the effects of sex, age, habitat type and suitability, and local harvest density at the natal and settlement sites on the probability of natal dispersal, and on dispersal distances. Dispersal was male-biased (76% of males dispersed) but a small proportion (21%) of females also dispersed, and female dispersal distances (mean ± SE  =  48.9±7.7 km) were comparable to male dispersal distances (59.0±3.2 km). Dispersal probabilities and dispersal distances were greatest for bears in areas with high habitat suitability and low harvest density. The inverse relationship between dispersal and harvest density in black bears suggests that 1) intensive harvest promotes restricted dispersal, or 2) high black bear population density decreases the propensity to disperse. Multigenerational genetic data collected over large landscape scales can be a powerful means of characterizing dispersal patterns and causal associations with demographic and landscape features in wild populations of elusive and wide-ranging species.

Local density and group size interacts with age and sex to determine direction and rate of social dispersal in a polygynous mammal

Movement away from an area or social group in response to increasing density (density-dependent dispersal) is known for most species; why it evolves is fundamental to our understanding of ecology and evolution. However, we have yet to fully appreciate how individuals of varying conditions (e.g., age and sex) might differently consider effects of density (quorum) when deciding to disperse or not, and scale dependence in their sense of quorum. We tracked movements of all individuals of a naturalized population of feral horses (Equus ferus caballus; Sable Island National Park Reserve, Nova Scotia, Canada) during a period of rapid population growth (N increased from 375 to 484 horses from 2008 to 2010). Permanent dispersal from breeding groups (bands) was positively density dependent for all age and sex categories with respect to local density (horses/km², bounded by the 99th percentile of individual movements [8000 m]), but was negatively and positively density dependent for males and females, respectively, in relation to group (band) size. Dispersal was generally female biased, with the exception of foals which moved with their mothers (no sex effect), and for yearlings and subadults when band sizes were smaller than average, in which case males dispersed at higher rates than females. Dispersal distance was positively related to local density. We conclude that dispersal rate can be both positively and negatively density dependent for feral horses, contingent on the state of individuals and the scale at which quorum with respect to choosing to disperse or not is assessed. Scale effects and interactions of density-dependent and sex- and age-biased dispersal may have both ecological and evolutionary consequences through effects on resource and mate competition.