Influence of drift and admixture on population structure of American black bears (Ursus americanus) in the Central Interior Highlands, USA, 50 years after translocation

Gut microbiome as a key monitoring indicator for reintroductions of captive animals

Reintroduction programs seek to restore degraded populations and reverse biodiversity loss. To examine the hypothesis that gut symbionts could be used as an indicator of reintroduction success, we performed intensive metagenomic monitoring over 10 years to characterize the ecological succession and adaptive evolution of the gut symbionts of captive giant pandas reintroduced to the wild. We collected 63 fecal samples from 3 reintroduced individuals and 22 from 9 wild individuals and used 96 publicly available samples from another 3 captive individuals. By microbial composition analysis, we identified 3 community clusters of the gut microbiome (here termed enterotypes) with interenterotype succession that was closely related to the reintroduction process. Each of the 3 enterotypes was identified based on significant variation in the levels of 1 of 3 genera: Clostridium, Pseudomonas, and Escherichia. The enterotype of captive pandas was Escherichia. This enterotype was gradually replaced by the Clostridium enterotype during the wild-training process, which in turn was replaced by the Pseudomonas enterotype that resembled the enterotype of wild pandas, an indicator of conversion to wildness and a successful reintroduction. We also isolated 1 strain of Pseudomonas protegens from the wild enterotype, a previously reported free-living microbe, and found that its within-host evolution contributed to host dietary adaptation in the wild. Monitoring gut microbial structure provides a novel, noninvasive tool that can be used as an indicator of successful reintroduction of a captive individual to the wild.

Genetic mixing in conservation translocations increases diversity of a keystone threatened species, Bettongia lesueur

Translocation programmes are increasingly being informed by genetic data to monitor and enhance conservation outcomes for both natural and established populations. These data provide a window into contemporary patterns of genetic diversity, structure and relatedness that can guide managers in how to best source animals for their translocation programmes. The inclusion of historical samples, where possible, strengthens monitoring by allowing assessment of changes in genetic diversity over time and by providing a benchmark for future improvements in diversity via management practices. Here, we used reduced representation sequencing (ddRADseq) data to report on the current genetic health of three remnant and seven translocated boodie (Bettongia lesueur) populations, now extinct on the Australian mainland. In addition, we used exon capture data from seven historical mainland specimens and a subset of contemporary samples to compare pre-decline and current diversity. Both data sets showed the significant impact of population founder source (whether multiple or single) on the genetic diversity of translocated populations. Populations founded by animals from multiple sources showed significantly higher genetic diversity than the natural remnant and single-source translocation populations, and we show that by mixing the most divergent populations, exon capture heterozygosity was restored to levels close to that observed in pre-decline mainland samples. Relatedness estimates were surprisingly low across all contemporary populations and there was limited evidence of inbreeding. Our results show that a strategy of genetic mixing has led to successful conservation outcomes for the species in terms of increasing genetic diversity and provides strong rationale for mixing as a management strategy.

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.

Urbanization and Population Genetic Structure of the Panama City crayfish (Procambarus econfinae)

For species with geographically restricted distributions, the impacts of habitat loss and fragmentation on long-term persistence may be particularly pronounced. We examined the genetic structure of Panama City crayfish (PCC), Procambarus econfinae, whose historical distribution is limited to an area approximately 145 km2, largely within the limits of Panama City and eastern Bay County, FL. Currently, PCC occupy approximately 28% of its historical range, with suitable habitat composed of fragmented patches in the highly urbanized western portion of the range and managed plantations in the more contiguous eastern portion of the range. We used 1640 anonymous single-nucleotide polymorphisms to evaluate the effects of anthropogenic habitat modification on the genetic diversity and population structure of 161 PCC sampled from across its known distribution. First, we examined urban habitat patches in the west compared with less-developed habitat patches in the east. Second, we used approximate Bayesian computation to model inferences on the demographic history of eastern and western populations. We found anthropogenic habitat modifications explain the genetic structure of PCC range-wide. Clustering analyses revealed significant genetic structure between and within eastern and western regions. Estimates of divergence between east and west were consistent with urban growth in the mid-20th century. PCC have low genetic diversity and high levels of inbreeding and relatedness, indicating populations are small and isolated. Our results suggest that PCC have been strongly affected by habitat loss and fragmentation and management strategies, including legal protection, translocations, or reintroductions, may be necessary to ensure long-term persistence.

Population genetics of the wolverine in Finland: the road to recovery?

After decades, even centuries of persecution, large carnivore populations are widely recovering in Europe. Considering the recent recovery of the wolverine (Gulo gulo) in Finland, our aim was to evaluate genetic variation using 14 microsatellites and mtDNA control region (579 bp) in order (1) to determine whether the species is represented by a single genetic population within Finland, (2) to quantify the genetic diversity, and (3) to estimate the effective population size. We found two major genetic clusters divided between eastern and northern Finland based on microsatellites (FST = 0.100) but also a significant pattern of isolation by distance. Wolverines in western Finland had a genetic signature similar to the northern cluster, which can be explained by former translocations of wolverines from northern to western Finland. For both main clusters, most estimates of the effective population size Ne were below 50. Nevertheless, the genetic diversity was higher in the eastern cluster (HE = 0.57, AR = 4.0, AP = 0.3) than in the northern cluster (HE = 0.49, AR = 3.7, AP = 0.1). Migration between the clusters was low. Two mtDNA haplotypes were found: one common and identical to Scandinavian wolverines; the other rare and not previously detected. The rare haplotype was more prominent in the eastern genetic cluster. Combining all available data, we infer that the genetic population structure within Finland is shaped by a recent bottleneck, isolation by distance, human-aided translocations and postglacial recolonization routes.

Giraffe translocation population viability analysis

Most populations of giraffes have declined in recent decades, leading to the recent IUCN decision to upgrade the species to Vulnerable status, and some subspecies to Endangered. Translocations have been used as a conservation tool to re-introduce giraffes to previously occupied areas or establish new populations, but guidelines for founding populations are lacking. To provide general guidelines for translocation projects regarding feasibility, we simulated various scenarios of translocated giraffe populations to identify viable age and sex distributions of founding populations using population viability analysis (PVA) implemented in Vortex software. We explored the parameter space for demography and the genetic load, examining how variation in founding numbers and sex ratios affected 100 yr probability of population extinction and genetic diversity. We found that even very small numbers of founders (N ≤ 10 females) can appear to be successful in the first decades due to transient positive population growth, but with moderate population growth rate and moderate genetic load, long-term population viability (probability of extinction <0.01) was only achieved with ≥30 females and ≥3 males released. To maintain >95% genetic diversity of the source population in an isolated population, 50 females and 5 males are recommended to compose the founding population. Sensitivity analyses revealed first-year survival and reproductive rate were the simulation parameters with the greatest proportional influence on probability of extinction and genetic diversity. These simulations highlight important considerations for translocation success and data gaps including true genetic load in wild giraffe populations.

UrsaPlex: An STR multiplex for forensic identification of North American black bear (Ursus americanus)

UrsaPlex is a forensic quality 5-dye multiplexed tetranucleotide STR and sex identification panel for individual genetic identification of North American black bears (Ursus americanus). The panel is validated for the identification of black bears involved in human-wildlife conflict events and poaching investigations. This is the first single multiplex panel composed solely of tetranucleotide STRs derived from black bear and bear-specific sex markers. UrsaPlex produces complete genetic profiles from as little as 78 pg of DNA template and has a probability of identity of 2.63 × 10^-13. The panel has also been tested for utility in other ursids, and our results indicate with minor modifications, UrsaPlex should prove valuable in identification investigations involving these species as well.

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.

Modeling Heterogeneity in the Genetic Architecture of Ethnically Diverse Groups Using Random Effect Interaction Models

In humans, most genome-wide association studies have been conducted using data from Caucasians and many of the reported findings have not replicated in other populations. This lack of replication may be due to statistical issues (small sample sizes or confounding) or perhaps more fundamentally to differences in the genetic architecture of traits between ethnically diverse subpopulations. What aspects of the genetic architecture of traits vary between subpopulations and how can this be quantified? We consider studying effect heterogeneity using Bayesian random effect interaction models. The proposed methodology can be applied using shrinkage and variable selection methods, and produces useful information about effect heterogeneity in the form of whole-genome summaries (e.g., the proportions of variance of a complex trait explained by a set of SNPs and the average correlation of effects) as well as SNP-specific attributes. Using simulations, we show that the proposed methodology yields (nearly) unbiased estimates when the sample size is not too small relative to the number of SNPs used. Subsequently, we used the methodology for the analyses of four complex human traits (standing height, high-density lipoprotein, low-density lipoprotein, and serum urate levels) in European-Americans (EAs) and African-Americans (AAs). The estimated correlations of effects between the two subpopulations were well below unity for all the traits, ranging from 0.73 to 0.50. The extent of effect heterogeneity varied between traits and SNP sets. Height showed less differences in SNP effects between AAs and EAs whereas HDL, a trait highly influenced by lifestyle, exhibited a greater extent of effect heterogeneity. For all the traits, we observed substantial variability in effect heterogeneity across SNPs, suggesting that effect heterogeneity varies between regions of the genome.

How many came home? Evaluating ex situ conservation of green turtles in the Cayman Islands

Ex situ management is an important conservation tool that allows the preservation of biological diversity outside natural habitats while supporting survival in the wild. Captive breeding followed by re-introduction is a possible approach for endangered species conservation and preservation of genetic variability. The Cayman Turtle Centre Ltd was established in 1968 to market green turtle (Chelonia mydas) meat and other products and replenish wild populations, thought to be locally extirpated, through captive breeding. We evaluated the effects of this re-introduction programmme using molecular markers (13 microsatellites, 800-bp D-loop and simple tandem repeat mitochondrial DNA sequences) from captive breeders (N = 257) and wild nesting females (N = 57) (sampling period: 2013-2015). We divided the captive breeders into three groups: founders (from the original stock), and then two subdivisions of F₁ individuals corresponding to two different management strategies, cohort 1995 ("C1995") and multicohort F₁ ("MCF1"). Loss of genetic variability and increased relatedness was observed in the captive stock over time. We found no significant differences in diversity among captive and wild groups, and similar or higher levels of haplotype variability when compared to other natural populations. Using parentage and sibship assignment, we determined that 90% of the wild individuals were related to the captive stock. Our results suggest a strong impact of the re-introduction programmme on the present recovery of the wild green turtle population nesting in the Cayman Islands. Moreover, genetic relatedness analyses of captive populations are necessary to improve future management actions to maintain genetic diversity in the long term and avoid inbreeding depression.

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.

Turtle soup, Prohibition, and the population genetic structure of Diamondback Terrapins (Malaclemys terrapin)

Diamondback terrapins (Malaclemys terrapin) were a popular food item in early twentieth century America, and were consumed in soup with sherry. Intense market demand for terrapin meat resulted in population declines, notably along the Atlantic seaboard. Efforts to supply terrapins to markets resulted in translocation events, as individuals were moved about to stock terrapin farms. However, in 1920 the market for turtle soup buckled with the enactment of the eighteenth amendment to the United States’ Constitution—which initiated the prohibition of alcoholic drinks—and many terrapin fisheries dumped their stocks into local waters. We used microsatellite data to show that patterns of genetic diversity along the terrapin’s coastal range are consistent with historical accounts of translocation and cultivation activities. We identified possible instances of human-mediated dispersal by estimating gene flow over historical and contemporary timescales, Bayesian model testing, and bottleneck tests. We recovered six genotypic clusters along the Gulf and Atlantic coasts with varying degrees of admixture, including increased contemporary gene flow from Texas to South Carolina, from North Carolina to Maryland, and from North Carolina to New York. In addition, Bayesian models incorporating translocation events outperformed stepping-stone models. Finally, we were unable to detect population bottlenecks, possibly due to translocation reintroducing genetic diversity into bottlenecked populations. Our data suggest that current patterns of genetic diversity in the terrapin were altered by the demand for turtle soup followed by the enactment of alcohol prohibition. In addition, our study shows that population genetic tools can elucidate metapopulation dynamics in taxa with complex genetic histories impacted by anthropogenic activities.

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.

Population genetic structure and comparative diversity of smallmouth bass Micropterus dolomieu: congruent patterns from two genomes

Genetic diversity and divergence patterns of smallmouth bass Micropterus dolomieu spawning groups are analysed across its northern native range with mtDNA cytochrome b gene sequences and eight unlinked nuclear DNA microsatellite loci. Results reveal high levels of genetic variability and significant differences in allelic representation among populations (mtDNA: mean ± s.e., H_D = 0·50 ± 0·06, mean ± s.e., θ_ST = 0·41 ± 0·02 and microsatellites: mean ± s.e. H_O = 0·46 ± 0·03, mean ± s.e. θ_ST = 0·25 ± 0·01). The distributions of 28 variant mtDNA haplotypes, which differ by an average of 3·94 nucleotides (range = 1-8), denote divergent representation among geographic areas. Microsatellite data support nine primary population groups, whose high self-assignment probabilities likewise display marked divergence. Genetic patterns demonstrate: (1) high genetic diversity in both genomes, (2) significant divergence among populations, probably resulting from natal site homing and low lifetime migration, (3) support for three post-glacial refugia that variously contributed to the current northern populations, which remain evident today despite waterway connectivity and (4) a weak yet significant genetic isolation by geographic distance pattern, indicating that other processes affect the differences among populations, such as territoriality and site fidelity.

Black bear (Ursus americanus) functional resource selection relative to intraspecific competition and human risk

The spatial scales at which animals make behavioral trade-offs is assumed to relate to the scales at which factors most limiting resources and increasing mortality risk occur. We used global positioning system collar locations of 29 reproductive-age female black bears (Ursus americanus Pallas, 1780) in three states to assess resource selection relative to bear population-specific density, an index of vegetation productivity, riparian corridors, or two road classes of and within home ranges during spring–summer of 2009–2013. Female resource selection was best explained by functional responses to vegetation productivity across nearly all populations and spatial scales, which appeared to be influenced by variation in bear density (i.e., intraspecific competition). Behavioral trade-offs were greatest at the landscape scale, but except for vegetation productivity, were consistent for populations across spatial scales. Females across populations selected locations nearer to tertiary roads, but females in Michigan and Mississippi selected main roads and avoided riparian corridors, whereas females in Missouri did the opposite, suggesting population-level trade-offs between resource (e.g., food) acquisition and mortality risks (e.g., vehicle collisions). Our study emphasizes that female bear population-level resource selection can be influenced by multiple spatially dependent factors, and that scale-dependent functional behavior should be identified for management of bears across their range.

Spatiotemporal factors affecting detection of black bears during noninvasive capture–recapture surveys

Accounting for low and heterogeneous detection probabilities in large mammal capture–recapture sampling designs is a persistent challenge. Our objective was to improve understanding of ecological and biological factors driving detection using multiple data sources from an American black bear ( Ursus americanus ) DNA hair trap study in south-central Missouri. We used Global Positioning System telemetry and remote camera data to examine how a bear’s distance to traps, probability of space use, sex-specific behavior, and temporal sampling frame affect detection probability and number of hair samples collected at hair traps. Regression analysis suggested that bear distance to nearest hair trap was the best predictor of detection probability and indicated that detection probability at encounter was 0.15 and declined to &lt; 0.05 at nearest distances &gt; 330 m from hair traps. From remote camera data, number of hair samples increased with number of visits, but the proportion of hair samples from known visits declined 39% from early June to early August. Bears appeared attracted to lured hair traps from close distances and we recommend a hair trap density of 1 trap/2.6 km 2 with spatial coverage that encompasses potentially large male home ranges. We recommend sampling during the late spring and early summer molting period to increase hair deposition rates.

Trap array configuration influences estimates and precision of black bear density and abundance

Spatial capture-recapture (SCR) models have advanced our ability to estimate population density for wide ranging animals by explicitly incorporating individual movement. Though these models are more robust to various spatial sampling designs, few studies have empirically tested different large-scale trap configurations using SCR models. We investigated how extent of trap coverage and trap spacing affects precision and accuracy of SCR parameters, implementing models using the R package secr. We tested two trapping scenarios, one spatially extensive and one intensive, using black bear (Ursus americanus) DNA data from hair snare arrays in south-central Missouri, USA. We also examined the influence that adding a second, lower barbed-wire strand to snares had on quantity and spatial distribution of detections. We simulated trapping data to test bias in density estimates of each configuration under a range of density and detection parameter values. Field data showed that using multiple arrays with intensive snare coverage produced more detections of more individuals than extensive coverage. Consequently, density and detection parameters were more precise for the intensive design. Density was estimated as 1.7 bears per 100 km2 and was 5.5 times greater than that under extensive sampling. Abundance was 279 (95% CI = 193-406) bears in the 16,812 km2 study area. Excluding detections from the lower strand resulted in the loss of 35 detections, 14 unique bears, and the largest recorded movement between snares. All simulations showed low bias for density under both configurations. Results demonstrated that in low density populations with non-uniform distribution of population density, optimizing the tradeoff among snare spacing, coverage, and sample size is of critical importance to estimating parameters with high precision and accuracy. With limited resources, allocating available traps to multiple arrays with intensive trap spacing increased the amount of information needed to inform parameters with high precision.