Capwire: a R package for estimating population census size from non‐invasive genetic sampling

Assessing the role of individual foxes in environmental contamination with Echinococcus multilocularis through faecal samples

Key parasite transmission parameters are difficult to obtain from elusive wild animals. For Echinococcus multilocularis, the causative agent of alveolar echinococcosis (AE), the red fox is responsible for most of the environmental contamination in Europe. The identification of individual spreaders of E. multilocularis environmental contamination is crucial to improving our understanding of the ecology of parasite transmission in areas of high endemicity and optimising the effectiveness of prevention and control measures in the field. Genetic faecal sampling appears to be a feasible method to gain information about the faecal deposition of individual animals. We conducted a 4 year faecal sampling study in a village that is highly endemic for E. multilocularis, to assess the feasibility of individual identification and sexing of foxes to describe individual infection patterns. Individual fox identification from faecal samples was performed by obtaining reliable genotypes from 14 microsatellites and one sex locus, coupled with the detection of E. multilocularis DNA, first using captive foxes and then by environmental sampling. From a collection of 386 fox stools collected between 2017 and 2020, tested for the presence of E. multilocularis DNA, 180 were selected and 124 samples were successfully genotyped (68.9%). In total, 45 unique individual foxes were identified and 26 associated with at least one sample which tested positive for E. multilocularis (Em(+)). Estimation of the population size showed the fox population to be between 29 and 34 individuals for a given year and 67 individuals over 4 years. One-third of infected individuals (9/26 Em(+) foxes) deposited 2/3 of the faeces which tested positive for E. multilocularis (36/60 Em(+) stools). Genetic investigation showed a significantly higher average number of multiple stools for females than males, suggesting that the two sexes potentially defecated unequally in the studied area. Three partially overlapping clusters of fox faeces were found, with one cluster concentrating 2/3 of the total E. multilocularis-positive faeces. Based on these findings, we estimated that 12.5 million E. multilocularis eggs were produced during the study period, emphasizing the high contamination level of the environment and the risk of exposure faced by the parasite hosts.

Nutritional ecology of a prototypical generalist predator, the red fox (Vulpes vulpes)

Generalist species, which exploit a wide range of food resources, are expected to be able to combine available resources as to attain their specific macronutrient ratio (percentage of caloric intake of protein, lipids and carbohydrates). Among mammalian predators, the red fox Vulpes vulpes is a widespread, opportunistic forager: its diet has been largely studied, outlining wide variation according to geographic and climatic factors. We aimed to check if, throughout the species' European range, diets vary widely in macronutrient composition or foxes can combine complementary foods to gain the same nutrient intake. First, we assessed fox's intake target in the framework of nutritional geometry. Secondly, we aimed to highlight the effects of unbalanced diets on fox density, which was assumed as a proxy for Darwinian fitness, as assessed in five areas of the western Italian Alps. Unexpectedly, the target macronutrient ratio of the fox (52.4% protein-, 38.7% lipid- and 8.9% carbohydrate energy) was consistent with that of hypercarnivores, such as wolves and felids, except for carbohydrate intakes in urban and rural habitats. The inverse relation between density and the deviation of observed macronutrient ratios from the intake target suggests that fox capability of surviving in a wide range of habitats may not be exempt from fitness costs and that nutrient availability should be regarded among the biotic factors affecting animal abundance and distribution.

Genetic Analysis and Status of Brown Bear Sub-Populations in Three National Parks of Greece Functioning as Strongholds for the Species’ Conservation

In order to optimize the appropriate conservation actions for the brown bear (Ursus arctos L.) population in Greece, we estimated the census (Nc) and effective (Ne) population size as well as the genetic status of brown bear sub-populations in three National Parks (NP): Prespa (MBPNP), Pindos (PINDNP), and Rhodopi (RMNP). The Prespa and Pindos sub-populations are located in western Greece and the Rhodopi population is located in eastern Greece. We extracted DNA from 472 hair samples and amplified through PCR 10 microsatellite loci. In total, 257 of 472 samples (54.5%) were genotyped for 6–10 microsatellite loci. Genetic analysis revealed that the Ne was 35, 118, and 61 individuals in MBPNP, PINDNP, and RMNP, respectively, while high levels of inbreeding were found in Prespa and Rhodopi but not in Pindos. Moreover, analysis of genetic structure showed that the Pindos population is genetically distinct, whereas Prespa and Rhodopi show mutual overlaps. Finally, we found a notable gene flow from Prespa to Rhodopi (10.19%) and from Rhodopi to Prespa (14.96%). Therefore, targeted actions for the conservation of the bears that live in the abovementioned areas must be undertaken, in order to ensure the species’ viability and to preserve the corridors that allow connectivity between the bear sub-populations in Greece.

Faecal DNA-based genetic survey of a relict Eurasian otter (Lutra lutra) population (Sila Massif, S Italy)

Faecal DNA-based genetic analysis is a suitable tool for assessing both population size and genetic diversity of threatened and elusive species. We applied microsatellite analysis and mtDNA sequencing for investigating the southernmost Italian (Sila Massif, Calabria Region) population of the Eurasian otter (Lutra lutra). This relict population, filed as extinct in the mid-1980s, is currently expanding but still quite isolated. On the two main rivers hosting otters permanently since 2014, we collected 47 spraints, out of which 24 (51.1%) were successfully genotyped (on average 2.0 alleles per locus). Thirteen individuals were identified: seven females and three males (sex identification success of 76.9%). Population size was assessed as 16 individuals (13–22), corresponding to a density of 0.15 (0.13–0.21) ind/km. Successfully amplified mtDNA samples (N = 16) confirmed the occurrence of a haplotype—H10—which had been previously reported only for Southern Italy, bringing new evidence of the unicity of the Italian otter population. Although density values complied with those reported for the core area of otter Italian range, the small size and genetic isolation of this population require special attention. To assist the ongoing re-colonisation of the Sila Massif, habitat management should aim to enhance fish availability and connectivity with the core area.

A Noninvasive Genetic Insight into the Spatial and Social Organization of an Endangered Population of the Eurasian Otter (Lutra lutra, Mustelidae, Carnivora)

The Eurasian otter is endangered in Italy, only surviving in southern river basins. The spatial and social structure of a population living at the border of the current range was explored through a noninvasive genetic study along 174 km of the Sangro river. Sampling was conducted in 2011 and 2012, collecting spraints and anal jellies at 62 marking sites. Samples were successfully genotyped at 13 nuclear microsatellites and the ZFX/ZFY locus for molecular sexing, resulting in 14 distinct genotypes (4 females, 2 possible females, 8 males), from 35 marking sites. Mean captures/recaptures rate was 3.8 captures/individual, with males being recaptured more frequently than females. Spatial overlap among individuals was analyzed through a linear regression model fitted against sibship categories and sex pairing. Nine out of the fourteen genotyped individuals belonged to three full-sib clusters, while five individuals had no full-sibs in the population. Full-sibs overlapped more than half-sibs, while male–male pairs showed significantly higher spatial overlap than both male–female and female–female pairs. Estimated mean density was 0.152 otters/km and 2.4 individuals/10 × 10 km grid cell. Accordingly, the 3440 grid cells of otter occurrence in Italy could likely host about 8000 otters, suggesting the current population has become larger than the minimum viable population size.

An efficient method for simultaneous species, individual, and sex identification via in-solution single nucleotide polymorphism capture from low-quality scat samples

Understanding predator population dynamics is important for conservation management because of the critical roles predators play within ecosystems. Noninvasive genetic sampling methods are useful for the study of predators like canids that can be difficult to capture or directly observe. Here, we introduce the FAECES* method (Fast and Accurate Enrichment of Canid Excrement for Species* and other analyses) which expands the toolbox for canid researchers and conservationists by using in-solution hybridization sequence capture to produce single nucleotide polymorphism (SNP) genotypes for multiple canid species from scat-derived DNA using a single enrichment. We designed a set of hybridization probes to genotype both coyotes (Canis latrans) and kit foxes (Vulpes macrotis) at hundreds of polymorphic SNP loci and we tested the probes on both tissues and field-collected scat samples. We enriched and genotyped by sequencing 52 coyote and 70 kit fox scats collected in and around a conservation easement in the Nevada Mojave Desert. We demonstrate that the FAECES* method produces genotypes capable of differentiating coyotes and kit foxes, identifying individuals and their sex, and estimating genetic diversity and effective population sizes, even using highly degraded, low-quantity DNA extracted from scat. We found that the study area harbours a large and diverse population of kit foxes and a relatively smaller population of coyotes. By replicating our methods in the future, conservationists can assess the impacts of management decisions on canid populations. The method can also be adapted and applied more broadly to enrich and sequence multiple loci from any species of interest using scat or other noninvasive genetic samples.

Neonicotinoid Pesticides Cause Mass Fatalities of Native Bumble Bees: A Case Study From Wilsonville, Oregon, United States

In June of 2013 an application of dinotefuran on an ornamental planting of European linden trees (Tilia cordata Mill. [Malvales: Malvalceae]) in a shopping mall parking lot in Wilsonville, Oregon provoked the largest documented pesticide kill of bumble bees in North America. Based on geographic information systems and population genetic analysis, we estimate that between 45,830 and 107,470 bumble bees originating from between 289 and 596 colonies were killed during this event. Dinotefuran is a neonicotinoid that is highly effective in exterminating and/or harming target pest insects and non-target beneficial insects. Analysis to detect the concentration of pesticides in flowers that received foliar application revealed that the minimum reported dinotefuran concentration of a sampled T. cordata flower was 7.4 ppm, or in excess of 737% above the LC50 of the beneficial pollinator, the honey bee (Apis mellifera Linnaeus, 1758 [Hymenoptera: Apidae]). Furthermore, sampled Vosnesensky bumble bees (Bombus vosnesenskii Radoskowski, 1862 [Hymenoptera: Apidae]) were found to have an average dinotefuran concentration of 0.92 ppm at the time of death, which exceeds the maximum LC50 of A. mellifera (0.884 ppm). Our study underscores the lethal impact of the neonicotinoid pesticide dinotefuran on pollinating insect populations in a suburban environment. To our knowledge, the documentation and impact of pesticide kills on wild populations of beneficial insects has not been widely reported in the scientific literature. It is likely that the vast majority of mass pesticide kills of beneficial insects across other environments go unnoticed and unreported.

Spatial ecology of the stone marten in an Alpine area: combining camera-trapping and genetic surveys

A species’ potential distribution can be modelled adequately only if no factor other than habitat availability affects its occurrences. Space use by stone marten Martes foina is likely to be affected by interspecific competition with the strictly related pine marten Martes martes, the latter being able to outcompete the first species in forested habitats. Hence, to point out the environmental factors which determine the distribution and density of the stone marten, a relatively understudied mesocarnivore, we applied two non-invasive survey methods, camera-trapping and faecal-DNA based genetic analysis, in an Alpine area where the pine marten was deemed to be absent (Val Grande National Park N Italy). Camera trapping was conducted from October 2014 to November 2015, using up to 27 cameras. Marten scats were searched for between July and November 2015 and, to assess density, in spring 2017. Species identification was accomplished by a PCR-RFLP method, while 17 autosomal microsatellites were used for individual identification. The stone marten occurred in all available habitats (83% of trapping sites and 73.2% of scats); nonetheless, habitat suitability, as assessed using MaxEnt, depended on four major land cover variables—rocky grasslands, rocks and debris, beech forests and chestnut forests—, martens selecting forests and avoiding open rocky areas. Sixteen individuals were identified, of which 14 related to each other, possibly forming six different groups. Using capwire estimators, density was assessed as 0.95 (0.7–1.3) ind/km2. In the study area, the widespread stone marten selected forested areas, attaining density values like those reported for the pine marten in northern Europe and suggesting that patterns of habitat selection may depend on the relative abundance of the two competing martens.

Estimating red deer (Cervus elaphus) population size based on non-invasive genetic sampling

Some deer species are of conservation concern; others are officially managed as a food source or for their trophies, whereas in many regions, deer are regarded as overabundant or even as a nuisance causing damages. Regardless of local management issues, in most cases, reliable data on deer population sizes and sex ratios are lacking. Non-invasive genetic approaches are promising tools for the estimation of population size and structure. We developed and tested a non-invasive genetic approach for red deer (Cervus elaphus) population size and density estimation based on faeces collected from three free-ranging red deer populations in south-western Germany. Altogether, we genotyped 2762 faecal samples, representing 1431 different individuals. We estimated population density for both sexes separately using two different approaches: spatially explicit capture-recapture (SECR) approach and a single-session urn model (CAPWIRE). The estimated densities of both approaches were similar for all three study areas, ranging between total densities of 3.3 (2.5–4.4) and 8.5 (6.4–11.3) red deer/km2. The estimated sex ratios differed significantly between the studied populations (ranging between 1:1.1 and 1:1.7), resulting in considerable consequences for management. In further research, the issues of population closure and approximation of the effectively sampled area for density estimation should be addressed. The presented approach can serve as a valuable tool for the management of deer populations, and to our knowledge, it represents the only sex-specific approach for estimation of red deer population size and density.

Noninvasive Sampling Reveals Short-Term Genetic Rescue in an Insular Red Fox Population

Genetic factors in the decline of small populations are extremely difficult to study in nature. We leveraged a natural experiment to investigate evidence of inbreeding depression and genetic rescue in a remnant population of subalpine-specialized Sierra Nevada red foxes (Vulpes vulpes necator) using noninvasive genetic monitoring during 2010-2017. Only 7 individuals were detected in the first 2 years. These individuals assigned genetically to the historical population and exhibited genetic hallmarks of inbreeding and no evidence of reproduction. Two years into the study, we detected 2 first-generation immigrant males from a recently expanding population of red foxes in the Great Basin Desert. Through annual resampling of individuals (634 red fox DNA samples, 41 individuals) and molecular reconstruction of pedigrees, we documented 1-3 litters/year for 5 years, all descended directly or indirectly from matings involving immigrant foxes. The observed heterozygosity and allelic richness of the population nearly doubled in 2 years. Abundance increased, indicative of a rapidly expanding population. Throughout the study, adult survival was high. Restoration of gene flow apparently improved the demographic trajectory of this population in the short term. Whether these benefits continue in the longer term could depend on numerous factors, such as maintenance of any locally adapted alleles. This study highlights the value of noninvasive genetic monitoring to assess rapidly shifting conditions in small populations. Uncertainties about the longer-term trajectory of this population underscore the need to continue monitoring and to research potential for both negative and positive aspects of continued genetic infusion.

Wildfire reveals transient changes to individual traits and population responses of a native bumble bee Bombus vosnesenskii

Fire-induced changes in the abundance and distribution of organisms, especially plants, can alter resource landscapes for mobile consumers driving bottom-up effects on their population sizes, morphologies and reproductive potential. We expect these impacts to be most striking for obligate visitors of plants, like bees and other pollinators, but these impacts can be difficult to interpret due to the limited information provided by forager counts in the absence of survival or fitness proxies. Increased bumble bee worker abundance is often coincident with the pulses of flowers that follow recent fire. However, it is unknown if observed postfire activity is due to underlying population growth or a stable pool of colonies recruiting more foragers to abundant resource patches. This distinction is necessary for determining the net impact of disturbance on bumble bees: are there population-wide responses or do just a few colonies reap the rewards? We estimated colony abundance before and after fire in burned and unburned areas using a genetic mark-recapture framework. We paired colony abundance estimates with measures of body size, counts of queens, and estimates of foraging and dispersal to assess changes in worker size, reproductive output, and landscape-scale movements. Higher floral abundance following fire not only increased forager abundance but also the number of colonies from which those foragers came. Importantly, despite a larger population size, we also observed increased mean worker size. Two years following fire, queen abundance was higher in both burned and unburned sites, potentially due to the dispersal of queens from burned into unburned areas. The effects of fire were transient; within two growing seasons, worker abundance was substantially reduced across the entire sampling area and body sizes were similar between burned and unburned sites. Our results reveal how disturbance can temporarily release populations from resource limitation, boosting the genetic diversity, body size, and reproductive output of populations. Given that the effects of fire on bumble bees acted indirectly through pulsed resource availability, it is likely our results are generalizable to other situations, such as habitat restorations, where resource density is enhanced within the landscape.

Abundance, density, and social structure of African forest elephants (Loxodonta cyclotis) in a human-modified landscape in southwestern Gabon

Population monitoring is critical to effective conservation, but forest living taxa can be difficult to directly observe. This has been true of African forest elephants (Loxodonta cyclotis), for which we have limited information regarding population size and social behavior despite their threatened conservation status. In this study, we estimated demographic parameters using genetic capture-recapture of forest elephants in the southern Industrial Corridor of the Gamba Complex of Protected Areas in southwestern Gabon, which is considered a global stronghold for forest elephants. Additionally, we examined social networks, predicting that we would find matrilineal structure seen in both savanna and forest elephants. Given 95% confidence intervals, we estimate population size in the sampled area to be between 754 and 1,502 individuals and our best density estimate ranges from 0.47 to 0.80 elephants per km². When extrapolated across the entire Industrial Corridor, this estimate suggests an elephant population size of 3,033 to 6,043 based on abundance or 1,684 to 2,832 based on density, approximately 40–80% smaller than previously suggested. Our social network analysis revealed approximately half of network components included females with different mitochondrial haplotypes suggesting a wider range of variation in forest elephant sociality than previously thought. This study emphasizes the threatened status of forest elephants and demonstrates the need to further refine baseline estimates of population size and knowledge on social behavior in this taxon, both of which will aid in determining how population dynamics in this keystone species may be changing through time in relation to increasing conservation threats.

Marking behaviour and census of Eurasian otters (Lutra lutra) in riverine habitats: what can scat abundances and non-invasive genetic sampling tell us about otter numbers?

Growing human-wildlife conflicts and legal conservation obligations increased the need for precise information on Eurasian otter (Lutra lutra) population parameters for species status assessment and wildlife management measures. Scat surveys have become the method of choice to monitor species distribution range, abundance and habitat use. Although methodological concerns exist, scat abundance is often used as main indicator for otter population estimates and trends. To evaluate whether scat counts serve as good proxy for otter densities in linear stream habitats, we modelled the relationship between old/fresh scat abundance and otter numbers on two spatial scales, including also seasonal effects. Actual otter densities and marking behaviour were ascertained by genotyping of otter faeces collected at 218 marking sites along 132.6 km of four salmonid streams. Otter densities in our study areas ranged from 0.16 to 0.28 otters per stream kilometre. The age of scats and the level of spatial scale were determining factors when evaluating the suitability of scat abundance as index for otter densities. Covering large parts of an otter’s territory, fresh scat abundance increased with otter density, irrespectively of season. On the spatial scale of single marking sites, the relationship was not well supported, and no relationship could be found when using old scats in the models. Otter marking behaviour and differences in marking site use intensity, as well as seasonal differences in traceability and accumulation of old scats were discussed to bias the relationship between scat abundance and otter numbers.

Noninvasive population assessment of moose (Alces alces) by SNP genotyping of fecal pellets

Noninvasive genetic studies of wild animals enable the recovery of information infeasible to obtain using other means. However, the low quantity and quality of noninvasively collected DNA often challenge the retrieval of reliable genotypes, which may cause biases in downstream analyses. In this study, we optimized SNP (single nucleotide polymorphism) genotyping of fecal samples from moose (Alces alces) with the main purpose of exploring the potential of using noninvasively retrieved genotypes for individual- and sex identification. Fecal pellets were collected during the late winter of 2016 on the Swedish island of Öland in the Baltic Sea and DNA was extracted and genotyped using 86 autosomal, six sex-specific and five species diagnostic SNPs. The SNP error rate of the quality filtered dataset was 0.06 and the probability of identity for siblings below 0.001. Following a thorough quality filtering process, 182 reliable genotypes were obtained, corresponding to 100 unique individuals (37 males, 63 females), with an estimated male proportion of 37% (± 9%). The population size, estimated using two different capture-mark-recapture approaches, was found to be in the range of 115–156 individuals (95% CI). Furthermore, moose on Öland showed significantly lower heterozygosity levels (zHexp = −5.51, N = 69, pHexp = 3.56·10−8, zHobs = −3.58, N = 69, pHobs = 3.38·10−4) and appeared genetically differentiated from moose on the Swedish mainland. Thus, we show that quality controlled noninvasively derived SNP genotypes can be highly informative for individual and population monitoring in a large ungulate.

Population genetics reveals high connectivity of giant panda populations across human disturbance features in key nature reserve

The giant panda is an example of a species that has faced extensive historical habitat fragmentation, and anthropogenic disturbance and is assumed to be isolated in numerous subpopulations with limited gene flow between them. To investigate the population size, health, and connectivity of pandas in a key habitat area, we noninvasively collected a total of 539 fresh wild giant panda fecal samples for DNA extraction within Wolong Nature Reserve, Sichuan, China. Seven validated tetra-microsatellite markers were used to analyze each sample, and a total of 142 unique genotypes were identified. Nonspatial and spatial capture-recapture models estimated the population size of the reserve at 164 and 137 individuals (95% confidence intervals 153-175 and 115-163), respectively. Relatively high levels of genetic variation and low levels of inbreeding were estimated, indicating adequate genetic diversity. Surprisingly, no significant genetic boundaries were found within the population despite the national road G350 that bisects the reserve, which is also bordered with patches of development and agricultural land. We attribute this to high rates of migration, with four giant panda road-crossing events confirmed within a year based on repeated captures of individuals. This likely means that giant panda populations within mountain ranges are better connected than previously thought. Increased development and tourism traffic in the area and throughout the current panda distribution pose a threat of increasing population isolation, however. Maintaining and restoring adequate habitat corridors for dispersal is thus a vital step for preserving the levels of gene flow seen in our analysis and the continued conservation of the giant panda meta-population in both Wolong and throughout their current range.

Comparing performance of multiple non-invasive genetic capture–recapture methods for abundance estimation: a case study with the Sonoran pronghorn Antilocapra americana sonoriensis

Demographic monitoring is required in threatened species management, yet effective and efficient monitoring is challenging for species that are difficult to capture or susceptible to capture stress. One possible monitoring approach for such species is non-invasive genetic sampling with capture–recapture methods (genetic capture–recapture). We evaluated the performance of genetic capture–recapture in a challenging model system, monitoring the threatened Sonoran pronghorn Antilocapra americana sonoriensis. In an effort to determine the best (i.e. efficient, accurate, precise, cost-effective) method for abundance estimation, we used simulations to examine the optimal genetic capture–recapture faecal sampling design for this population. We simulated encounter histories for 100–300 individuals, with 0.33–3.33 samples/individual/session, in 1–3 sampling sessions. We explored trade-offs between sample size, number of sessions and multi-session (MARK) versus single-session (capwire) closed capture–recapture abundance estimators, and an accurate and precise estimate. We also compared the cost between the genetic capture–recapture approaches and current aerial monitoring methods. Abundance was biased positively in capwire and negatively in MARK. Bias increased and precision decreased with fewer samples/individual/session. Annual genetic capture–recapture monitoring cost was nearly twice the cost of aerial surveys, although genetic capture–recapture methods provided much higher precision. However at the current estimated abundance (c. 200), the same level of precision achieved with aerial methods can be obtained by collecting 0.75 samples/individual in a single session, for an annual cost saving of > USD 4,000. This approach of comparing estimator performance and cost can easily be applied to other systems and is a useful evaluation for managers to implement prior to designing capture–recapture studies.

Non-invasive genetic censusing and monitoring of primate populations

Knowing the density or abundance of primate populations is essential for their conservation management and contextualizing socio-demographic and behavioral observations. When direct counts of animals are not possible, genetic analysis of non-invasive samples collected from wildlife populations allows estimates of population size with higher accuracy and precision than is possible using indirect signs. Furthermore, in contrast to traditional indirect survey methods, prolonged or periodic genetic sampling across months or years enables inference of group membership, movement, dynamics, and some kin relationships. Data may also be used to estimate sex ratios, sex differences in dispersal distances, and detect gene flow among locations. Recent advances in capture-recapture models have further improved the precision of population estimates derived from non-invasive samples. Simulations using these methods have shown that the confidence interval of point estimates includes the true population size when assumptions of the models are met, and therefore this range of population size minima and maxima should be emphasized in population monitoring studies. Innovations such as the use of sniffer dogs or anti-poaching patrols for sample collection are important to ensure adequate sampling, and the expected development of efficient and cost-effective genotyping by sequencing methods for DNAs derived from non-invasive samples will automate and speed analyses.

Combining noninvasive genetics and a new mammalian sex-linked marker provides new tools to investigate population size, structure and individual behaviour: An application to bats

Monitoring wild populations is crucial for their effective management. Noninvasive genetic methods provide robust data from individual free-ranging animals, which can be used in capture-mark-recapture (CMR) models to estimate demographic parameters without capturing or disturbing them. However, sex- and status-specific behaviour, which may lead to differences in detection probabilities, is rarely considered in monitoring. Here, we investigated population size, sex ratio, sex- and status-related behaviour in 19 Rhinolophus hipposideros maternity colonies (Northern France) with a noninvasive genetic CMR approach (using faeces) combined with parentage assignments. The use of the DDX3X/Y-Mam sexual marker designed in this study, which shows inter- and intrachromosomal length polymorphism across placental mammals, together with eight polymorphic microsatellite markers, produced high-quality genetic data with limited genotyping errors and allowed us to reliably distinguish different categories of individuals (males, reproductive and nonreproductive females) and to estimate population sizes. We showed that visual counts represent well-adult female numbers and that population composition in maternity colonies changes dynamically during the summer. Before parturition, colonies mainly harbour pregnant and nonpregnant females with a few visiting males, whereas after parturition, colonies are mainly composed of mothers and their offspring with a few visiting nonmothers and males. Our approach gives deeper insight into sex- and status-specific behaviour, a prerequisite for understanding population dynamics and developing effective monitoring and management strategies. Provided sufficient samples can be obtained, this approach can be readily applied to a wide range of species.

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

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

Environmental DNA from Residual Saliva for Efficient Noninvasive Genetic Monitoring of Brown Bears (Ursus arctos)

Noninvasive genetic sampling is an important tool in wildlife ecology and management, typically relying on hair snaring or scat sampling techniques, but hair snaring is labor and cost intensive, and scats yield relatively low quality DNA. New approaches utilizing environmental DNA (eDNA) may provide supplementary, cost-effective tools for noninvasive genetic sampling. We tested whether eDNA from residual saliva on partially-consumed Pacific salmon (Oncorhynchus spp.) carcasses might yield suitable DNA quality for noninvasive monitoring of brown bears (Ursus arctos). We compared the efficiency of monitoring brown bear populations using both fecal DNA and salivary eDNA collected from partially-consumed salmon carcasses in Southeast Alaska. We swabbed a range of tissue types from 156 partially-consumed salmon carcasses from a midseason run of lakeshore-spawning sockeye (O. nerka) and a late season run of stream-spawning chum (O. keta) salmon in 2014. We also swabbed a total of 272 scats from the same locations. Saliva swabs collected from the braincases of salmon had the best amplification rate, followed by swabs taken from individual bite holes. Saliva collected from salmon carcasses identified unique individuals more quickly and required much less labor to locate than scat samples. Salmon carcass swabbing is a promising method to aid in efficient and affordable monitoring of bear populations, and suggests that the swabbing of food remains or consumed baits from other animals may be an additional cost-effective and valuable tool in the study of the ecology and population biology of many elusive and/or wide-ranging species.

Genetic censusing identifies an unexpectedly sizeable population of an endangered large mammal in a fragmented forest landscape

Background

As habitat degradation and fragmentation continue to impact wildlife populations around the world, it is critical to understand the behavioral flexibility of species in these environments. In Uganda, the mostly unprotected forest fragment landscape between the Budongo and Bugoma Forests is a potential corridor for chimpanzees, yet little is known about the status of chimpanzee populations in these fragments.

Results

From 2011 through 2013, we noninvasively collected 865 chimpanzee fecal samples across 633 km² and successfully genotyped 662 (77%) at up to 14 microsatellite loci. These genotypes corresponded to 182 chimpanzees, with a mean of 3.5 captures per individual. We obtained population size estimates of 256 (95% confidence interval 246–321) and 319 (288–357) chimpanzees using capture-with-replacement and spatially explicit capture–recapture models, respectively. The spatial clustering of associated genotypes suggests the presence of at least nine communities containing a minimum of 8–33 individuals each. Putative community distributions defined by the locations of associated genotypes correspond well with the distribution of 14 Y-chromosome haplotypes.

Conclusions

These census figures are more than three times greater than a previous estimate based on an extrapolation from small-scale nest count surveys that tend to underestimate population size. The distribution of genotype clusters and Y-chromosome haplotypes together indicate the presence of numerous male philopatric chimpanzee communities throughout the corridor habitat. Our findings demonstrate that, despite extensive habitat loss and fragmentation, chimpanzees remain widely distributed and exhibit distinct community home ranges. Our results further imply that elusive and rare species may adapt to degraded habitats more successfully than previously believed. Their long-term persistence is unlikely, however, if protection is not afforded to them and habitat loss continues unabated.

Electronic supplementary material

The online version of this article (doi:10.1186/s12898-015-0052-x) contains supplementary material, which is available to authorized users.

Evaluating noninvasive genetic sampling techniques to estimate large carnivore abundance

Monitoring large carnivores is difficult because of intrinsically low densities and can be dangerous if physical capture is required. Noninvasive genetic sampling (NGS) is a safe and cost-effective alternative to physical capture. We evaluated the utility of two NGS methods (scat detection dogs and hair sampling) to obtain genetic samples for abundance estimation of coyotes, black bears and Canada lynx in three areas of Newfoundland, Canada. We calculated abundance estimates using program capwire, compared sampling costs, and the cost/sample for each method relative to species and study site, and performed simulations to determine the sampling intensity necessary to achieve abundance estimates with coefficients of variation (CV) of <10%. Scat sampling was effective for both coyotes and bears and hair snags effectively sampled bears in two of three study sites. Rub pads were ineffective in sampling coyotes and lynx. The precision of abundance estimates was dependent upon the number of captures/individual. Our simulations suggested that ~3.4 captures/individual will result in a < 10% CV for abundance estimates when populations are small (23-39), but fewer captures/individual may be sufficient for larger populations. We found scat sampling was more cost-effective for sampling multiple species, but suggest that hair sampling may be less expensive at study sites with limited road access for bears. Given the dependence of sampling scheme on species and study site, the optimal sampling scheme is likely to be study-specific warranting pilot studies in most circumstances.

Real-time assessment of hybridization between wolves and dogs: combining noninvasive samples with ancestry informative markers

Wolves and dogs provide a paradigmatic example of the ecological and conservation implications of hybridization events between wild and domesticated forms. However, our understanding of such implications has been traditionally hampered by both high genetic similarity and the difficulties in obtaining tissue samples (TS), which limit our ability to assess ongoing hybridization events. To assess the occurrence and extension of hybridization in a pack of wolf-dog hybrids in northwestern Iberia, we compared the power of 52 nuclear markers implemented on TS with a subset of 13 ancestry informative markers (AIMs) typed in noninvasive samples (NIS). We demonstrate that the 13 AIMs are as accurate as the 52 markers that were chosen without regard to the power to differentiate between wolves and dogs, also having the advantage of being rapidly screened on NIS. The efficiency of AIMs significantly outperformed ten random sets of similar size and an additional commercial set of 18 markers. Bayesian clustering analysis implemented on AIMs and NIS identified nine hybrids, two wolves and two dogs. Four hybrids were unambiguously assigned to F1xWolf backcrosses. Our approach (AIMs + NIS) overcomes previous difficulties related to sample availability and informative power of markers, allowing a quick identification of wolf-dog hybrids in the first phases of hybridization episodes. This provides managers with a reliable tool to evaluate hybridization and estimate the success of their actions. This approach may be easily adapted for other pairs of wild/domesticated species, thus improving our understanding of the introgression of domestication genes into natural populations.