Non-invasive genetic censusing and monitoring of primate populations

Southern marsh deer (Blastocerus dichotomus) populations assessed using Amplicon Sequencing on fecal samples

Populations in isolated and small fragments lose genetic variability very fast and are usually of conservation concern because they are at greater risk of local extinction. The largest native deer in South America, Blastocerus dichotomus (Illiger, 1815), is a Vulnerable species according to the IUCN categorization, which inhabits tropical and subtropical swampy areas. In Argentina, its presence has been restricted to four isolated fragments. Here we examine the genetic diversity and differentiation among three of them, including the three different patches that form the southernmost population, using 18 microsatellite markers genotyped by Amplicon Sequencing of DNA extracted from fecal samples. Genetic diversity was low (HE < 0.45) in all three populations studied. We found three genetic clusters compatible with the geographic location of the samples. We also found a metapopulation dynamics that involves the patches that make up the southernmost population, with evidence of a barrier to gene flow between two of them. Our results point to the creation of a corridor as a necessary and urgent management action. This is the first study, at the population level, employing microsatellite genotyping by Amplicon Sequencing with non-invasive samples in an endangered species.

Future coexistence with great apes will require major changes to policy and practice

The great apes-bonobos, chimpanzees, gorillas and orangutans-are critically threatened by human activities. We have destroyed their habitats, hunted them and transmitted fatal diseases to them. Yet we also conduct research on them, try to protect them and live alongside them. They are endangered, and time is running out. Here we outline what must be done to ensure that future generations continue to share this planet with great apes. We urge dialogue with those who live with great apes and interact with them often. We advocate conservation plans that acknowledge the realities of climate change, economic drivers and population growth. We encourage researchers to use technology to minimize risks to great apes. Our proposals will require substantial investment, and we identify ways to generate these funds. We conclude with a discussion of how field researchers might alter their work to protect our closest living relatives more effectively.

Camera traps unveil demography, social structure, and home range of six unhabituated Western chimpanzee groups in the Moyen Bafing National Park, Guinea

Precise estimates of population dynamics and social grouping patterns are required for effective conservation of wild animal populations. It is difficult to obtain such information on non-human great apes as they have slow reproductive rates. To gain a better understanding of demography in these populations, previous research has typically involved habituation\, a process that requires years. Here, we collected data continuously over year-long periods to monitor an unhabituated population of critically endangered Western chimpanzees (Pan troglodytes verus) in the Moyen Bafing National Park, Guinea. We used two arrays of 100 camera traps that were placed opportunistically in two distinct 100 km2 sites, named Bakoun and Koukoutamba. We identified 227 individuals in Bakoun and 207 in Koukoutamba through their unique facial features. Our camera trap data make clear that these individuals belong to six and seven closed groups, respectively. Six of those groups were near-completely sampled with an average minimum size of 46.8 individuals (range: 37-58), and a mean adult sex ratio of 1.32 (range: 0.93-2.10). We described the demographic composition of these groups and use Bayesian social network analysis to understand population structure. The network analyses suggested that the social bonds within the two populations were structured by sex homophily, with male chimpanzees being more or equally likely to be observed together than other adult associations. Through estimation of minimum convex polygons, we described the minimum home range for those groups. Compared to other chimpanzee groups living in a similar environment (mosaic savanna-forest), the Moyen Bafing region seems to host a high-density of chimpanzees with small home ranges for their group size. Our research highlights the potential of camera traps for studying the demographic composition of chimpanzee populations with high resolution and obtaining crucial information on several groups in a time-efficient and cost-effective way.

Advances in wildlife abundance estimation using pedigree reconstruction

The conservation and management of wildlife populations, particularly for threatened and endangered species are greatly aided with abundance, growth rate, and density measures. Traditional methods of estimating abundance and related metrics represent trade-offs in effort and precision of estimates. Pedigree reconstruction is an emerging, attractive alternate approach because its use of one-time, noninvasive sampling of individuals to infer the existence of unsampled individuals. However, advances in pedigree reconstruction could improve its utility, including forming a measure of precision for the method, establishing required spatial sampling effort for accurate estimates, ascertaining the spatial extent of abundance estimates derived from pedigree reconstruction, and assessing how population density affects the estimator's performance. Using established relationships for a stochastic, spatially explicit simulated moose (Alces americanus) population, pedigree reconstruction provided accurate estimates of the adult moose population size and trend. Novel bootstrapped confidence intervals performed as expected with intensive sampling but underperformed with moderate sampling efforts that could produce abundance estimates with low bias. Adult population estimates more closely reflected the total number of adults in the extant population, rather than number of adults inhabiting the area where sampling occurred. Increasing sampling effort, measured as the proportion of individuals sampled and as the proportion of a hypothetical study area, yielded similar asymptotic patterns over time. Simulations indicated a positive relationship between animal density and sampling effort required for unbiased estimates. These results indicate that pedigree reconstruction can produce accurate abundance estimates and may be particularly valuable for surveying smaller areas and low-density populations.

Estimating animal abundance at multiple scales by spatially explicit capture-recapture

Information about how animal abundance varies across landscapes is needed to inform management action but is costly and time-consuming to obtain; surveys of a single population distributed over a large area can take years to complete. Surveys employing small, spatially replicated sampling units improve efficiency, but statistical estimators rely on assumptions that constrain survey design or become less reasonable as larger areas are sampled. Efficient methods that avoid assumptions about similarity of detectability or density among replicates are therefore appealing. Using simulations and data from >3500 black bears sampled on 73 independent study areas in Ontario, Canada, we (1) quantified bias induced by unmodeled spatial heterogeneity in detectability and density; (2) evaluated novel, design-based estimators of average density across replicate study areas; and (3) evaluated two estimators of the variance of average density across study areas: an analytic estimator that assumed an underlying homogeneous spatial Poisson point process for the distribution of animals' activity centers, and an empirical estimator of variance across study areas. In simulations where detectability varied in space, assuming spatially constant detectability yielded density estimates that were negatively biased by 20% to 30%; estimating local detectability and density from local data and treating study areas as independent, equal replicates when estimating average density across study areas using the design-based estimator yielded unbiased estimates at local and landscape scales. Similarly, detectability of black bears varied among study areas and estimates of bear density at landscape scales were higher when no information was shared across study areas when estimating detectability. This approach also maximized precision (relative SEs of estimates of average black bear density ranged from 7% to 18%) and computational efficiency. In simulations, the analytic variance estimator was robust to threefold variation in local densities but the empirical estimator performed poorly. Conducting multiple, similar SECR surveys and treating them as independent replicates during analyses allowed us to efficiently estimate density at multiple scales and extents while avoiding biases caused by pooling spatially heterogeneous data. This approach enables researchers to address a wide range of ecological or management-related questions and is applicable with most types of SECR data.

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.

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.

205 newly assembled mitogenomes provide mixed evidence for rivers as drivers of speciation for Amazonian primates

Mitochondrial DNA remains a cornerstone for molecular ecology, especially for study species from which high‐quality tissue samples cannot be easily obtained. Methods using mitochondrial markers are usually reliant on reference databases, but these are often incomplete. Furthermore, available mitochondrial genomes often lack crucial metadata, such as sampling location, limiting their utility for many analyses. Here, we assembled 205 new mitochondrial genomes for platyrrhine primates, most from the Amazon and with known sampling locations. We present a dated mitogenomic phylogeny based on these samples along with additional published platyrrhine mitogenomes, and use this to assess support for the long‐standing riverine barrier hypothesis (RBH), which proposes that river formation was a major driver of speciation in Amazonian primates. Along the Amazon, Negro, and Madeira rivers, we found mixed support for the RBH. While we identified divergences that coincide with a river barrier, only some occur synchronously and also overlap with the proposed dates of river formation. The most compelling evidence is for the Amazon river potentially driving speciation within bearded saki monkeys (Chiropotes spp.) and within the smallest extant platyrrhines, the marmosets and tamarins. However, we also found that even large rivers do not appear to be barriers for some primates, including howler monkeys (Alouatta spp.), uakaris (Cacajao spp.), sakis (Pithecia spp.), and robust capuchins (Sapajus spp.). Our results support a more nuanced, clade‐specific effect of riverine barriers and suggest that other evolutionary mechanisms, besides the RBH and allopatric speciation, may have played an important role in the diversification of platyrrhines.

The complex Y-chromosomal history of gorillas

Studies of the evolutionary relationships among gorilla populations using autosomal and mitochondrial sequences suggest that male-mediated gene flow may have been important in the past, but data on the Y-chromosomal relationships among the gorilla subspecies are limited. Here, we genotyped blood and noninvasively collected fecal samples from 12 captives and 257 wild male gorillas of known origin representing all four subspecies (Gorilla gorilla gorilla, G. g. diehli, G. beringei beringei, and G. b. graueri) at 10 Y-linked microsatellite loci resulting in 102 unique Y-haplotypes for 224 individuals. We found that western lowland gorilla (G. g. gorilla) haplotypes were consistently more diverse than any other subspecies for all measures of diversity and comprised several genetically distinct groups. However, these did not correspond to geographical proximity and some closely related haplotypes were found several hundred kilometers apart. Similarly, our broad sampling of eastern gorillas revealed that mountain (G. b. beringei) and Grauer's (G. b. graueri) gorilla Y-chromosomal haplotypes did not form distinct clusters. These observations suggest structure in the ancestral population with subsequent mixing of differentiated haplotypes by male dispersal for western lowland gorillas, and postisolation migration or incomplete lineage sorting due to short divergence times for eastern gorillas.

A beacon of dung: using lemming (Lemmus lemmus) winter nests and DNA analysis of faeces to further understand predator–prey dynamics in Northern Sweden

The hypothesis that predation is the cause of the regular small rodent population oscillations observed in boreal and Arctic regions has long been debated. Within this hypothesis, it is proposed that the most likely predators to cause these destabilizing effects are sedentary specialists, with small mustelids being possible candidates. One such case would be the highly specialized least weasel (Mustela nivalis) driving the Norwegian lemming (Lemmus lemmus) cycle in Fennoscandia. These predators are often elusive and therefore distribution data can only be based on field signs, which is problematic when various mustelid species are sympatric, such as weasels and stoats (Mustela erminea). Here we present the results of using mustelid faeces in predated winter lemming nests to correctly identify the predator and thus discern which species exerts the strongest predation pressure on lemming winter populations. Samples were obtained during different phases in the lemming cycle, spanning 6 years, to account for different prey densities. Faecal mitochondrial DNA extraction and amplification of a 400-bp fragment was successful in 92/114 samples (81%); the sequencing of these samples proved that most predation occurrences (83%) could be attributed to the least weasel. These findings support the hypothesis that weasels in particular show high specificity in predation and could therefore be candidates to driving the lemming cycle in this area. We conclude that DNA analysis of faecal remains around predated nests can be a useful tool for further investigations concerning predator–prey interactions in the tundra.

Noninvasive Technologies for Primate Conservation in the 21st Century

Observing and quantifying primate behavior in the wild is challenging. Human presence affects primate behavior and habituation of new, especially terrestrial, individuals is a time-intensive process that carries with it ethical and health concerns, especially during the recent pandemic when primates are at even greater risk than usual. As a result, wildlife researchers, including primatologists, have increasingly turned to new technologies to answer questions and provide important data related to primate conservation. Tools and methods should be chosen carefully to maximize and improve the data that will be used to answer the research questions. We review here the role of four indirect methods—camera traps, acoustic monitoring, drones, and portable field labs—and improvements in machine learning that offer rapid, reliable means of combing through large datasets that these methods generate. We describe key applications and limitations of each tool in primate conservation, and where we anticipate primate conservation technology moving forward in the coming years.

FISH and Chimps: Insights into Frequency and Distribution of Sperm Aneuploidy in Chimpanzees (Pan troglodytes)

Numerical chromosomal aberrations in sperm are considered to be a major factor in infertility, early pregnancy loss and syndromes with developmental and cognitive disabilities in mammals, including primates. Despite numerous studies in human and farm animals, the incidence and importance of sperm aneuploidies in non-human primate remains mostly undetermined. Here we investigated the incidence and distribution of sperm aneuploidy in chimpanzees (Pan troglodytes), the species closest to human. We identify evolutionary conserved DNA sequences in human and chimpanzee and selected homologous sub-telomeric regions for all chromosomes to build custom probes and perform sperm-FISH analysis on more than 10,000 sperm nuclei per chromosome. Chimpanzee mean autosomal disomy rate was 0.057 ± 0.02%, gonosomes disomy rate was 0.198% and the total disomy rate was 1.497%. The proportion of X or Y gametes was respectively 49.94% and 50.06% for a ratio of 1.002 and diploidy rate was 0.053%. Our data provide for the first time an overview of aneuploidy in non-human primate sperm and shed new insights into the issues of aneuploidy origins and mechanisms.

Home range size in central chimpanzees (Pan troglodytes troglodytes) from Loango National Park, Gabon

Ranging behavior has been studied extensively in eastern (Pan troglodytes schweinfurthii) and western (P. t. verus) chimpanzees, but relatively little is known regarding home ranges of the other two subspecies (P. t. ellioti; P. t. troglodytes). In this study, we determined the home range size and space use of a habituated community (Rekambo) of central chimpanzees living in a habitat mosaic in Loango National Park, Gabon. Data on travel routes were collected during follows between January 2017 and April 2019 (N = 670,616 relocations, collected over 640 days and 5690 h of observation). We used three methods for calculating home range size (minimum convex polygon, kernel density estimation, and biased random bridges). We compare our estimates to those obtained from prior genetic and camera trap studies of the Rekambo community and contrast them with estimates from other chimpanzee communities of the four chimpanzee subspecies. Depending on the methodology used, the home range size of the Rekambo community ranged between 27.64 and 59.03 km². The location of the center of the home range remained relatively stable over the last decade, while the overall size decreased. The Rekambo home range is, therefore, one of the largest documented so far for chimpanzees outside savannah-woodland habitats. We discuss several explanations, including the presence of savannah, interspecies competition, and intercommunity interactions.

Current and Forthcoming Approaches for Benchmarking Genetic and Genomic Diversity

The current attrition of biodiversity extends beyond loss of species and unique populations to steady loss of a vast genomic diversity that remains largely undescribed. Yet the accelerating development of new techniques allows us to survey entire genomes ever faster and cheaper, to obtain robust samples from a diversity of sources including degraded DNA and residual DNA in the environment, and to address conservation efforts in new and innovative ways. Here we review recent studies that highlight the importance of carefully considering where to prioritize collection of genetic samples (e.g., organisms in rapidly changing landscapes or along edges of geographic ranges) and what samples to collect and archive (e.g., from individuals of little-known subspecies or populations, even of species not currently considered endangered). Those decisions will provide the sample infrastructure to detect the disappearance of certain genotypes or gene complexes, increases in inbreeding levels, and loss of genomic diversity as environmental conditions change. Obtaining samples from currently endangered, protected, and rare species can be particularly difficult, thus we also focus on studies that use new, non-invasive ways of obtaining genomic samples and analyzing them in these cases where other sampling options are highly constrained. Finally, biological collections archiving such samples face an inherent contradiction: their main goal is to preserve biological material in good shape so it can be used for scientific research for centuries to come, yet the technologies that can make use of such materials are advancing faster than collections can change their standardized practices. Thus, we also discuss current and potential new practices in biological collections that might bolster their usefulness for future biodiversity conservation research.

Characterization of 29 polymorphic microsatellite markers developed by genomic screening of Sumatran rhinoceros (Dicerorhinus sumatrensis)

Objective

The Sumatran rhinoceros is critically endangered, with fewer than 100 individuals surviving across its current range. Accurate census estimates of the remaining populations are essential for development and implementation of conservation plans. In order to enable molecular censusing, we here develop microsatellite markers with amplicon sizes of short length, appropriate for non-invasive fecal sampling.

Results

Due to limited sample quantity and potential lack of genome-wide diversity, Illumina sequence reads were generated from two Sumatran rhinoceros samples. Genomic screening identified reads with short tandem repeats and loci that were polymorphic within the dataset. Twenty-nine novel polymorphic microsatellite markers were characterized (A = 2.4; H_O = 0.30). These were sufficient to distinguish among individuals (P_ID < 0.0001), and to distinguish among siblings (P_ID(sib) < 0.0001). Among rhinos in Indonesia, almost all markers were established as polymorphic and effective for genotyping DNA from fecal samples. Notably, the markers amplified and displayed microsatellite polymorphisms using DNA extracted from 11 fecal samples collected non-invasively from wild Sumatran rhinoceros. These microsatellite markers provide an important resource for a census and genetic studies of wild Sumatran rhinos.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13104-021-05522-x.

Recent genetic connectivity and clinal variation in chimpanzees

Much like humans, chimpanzees occupy diverse habitats and exhibit extensive behavioural variability. However, chimpanzees are recognized as a discontinuous species, with four subspecies separated by historical geographic barriers. Nevertheless, their range-wide degree of genetic connectivity remains poorly resolved, mainly due to sampling limitations. By analyzing a geographically comprehensive sample set amplified at microsatellite markers that inform recent population history, we found that isolation by distance explains most of the range-wide genetic structure of chimpanzees. Furthermore, we did not identify spatial discontinuities corresponding with the recognized subspecies, suggesting that some of the subspecies-delineating geographic barriers were recently permeable to gene flow. Substantial range-wide genetic connectivity is consistent with the hypothesis that behavioural flexibility is a salient driver of chimpanzee responses to changing environmental conditions. Finally, our observation of strong local differentiation associated with recent anthropogenic pressures portends future loss of critical genetic diversity if habitat fragmentation and population isolation continue unabated.

A refined panel of 42 microsatellite loci to universally genotype catarrhine primates

Microsatellite genotyping is an important genetic method for a number of research questions in biology. Given that the traditional fragment length analysis using polyacrylamide gel or capillary electrophoresis has several drawbacks, microsatellite genotyping-by-sequencing (GBS) has arisen as a promising alternative. Although GBS mitigates many of the problems of fragment length analysis, issues with allelic dropout and null alleles often remain due to mismatches in primer binding sites and unnecessarily long PCR products. This is also true for GBS in catarrhine primates where cross-species amplification of loci (often human derived) is common.We therefore redesigned primers for 45 microsatellite loci based on 17 available catarrhine reference genomes. Next, we tested them in singleplex and different multiplex settings in a panel of species representing all major lineages of Catarrhini and further validated them in wild Guinea baboons (Papio papio) using fecal samples.The final panel of 42 microsatellite loci can efficiently be amplified with primers distributed into three amplification pools.With our microsatellite panel, we provide a tool to universally genotype catarrhine primates via GBS from different sample sources in a cost- and time-efficient way, with higher resolution, and comparability among laboratories and species.

Hare's affairs: Lessons learnt from a noninvasive genetic monitoring for tracking mountain hare individuals

Systematic monitoring of individuals and their abundance over time has become an important tool to provide information for conservation. For genetic monitoring studies, noninvasive sampling has emerged as a valuable approach, particularly so for elusive or rare animals. Here, we present the 5-year results of an ongoing noninvasive genetic monitoring of mountain hares (Lepus timidus) in a protected area in the Swiss Alps. We used nuclear microsatellites and a sex marker to identify individuals and assign species to noninvasively collected feces samples. Through including a marker for sex identification, we were able to assess sex ratio changes and sex-specific demographic parameters over time. Male abundance in the area showed high fluctuations and apparent survival for males was lower than for females. Generally, males and females showed only little temporary migration into and out of the study area. Additionally, using genotyped tissue samples from mountain hares, European hares (Lepus europaeus) and their hybrids, we were able to provide evidence for the first occurrence of a European hare in the study area at an elevation of 2,300 m a.s.l. in spring 2016. For future monitoring studies, we suggest to include complementary analysis methods to reliably infer species identities of the individuals analyzed and, thus, not only monitor mountain hare individual abundance, but also assess the potential threats given through competitive exclusion by and hybridization with the European hare.

Three decades of Philippine nonhuman primate studies: research gaps and opportunities for Philippine primatology

The Philippines is considered a megadiverse country. However, there are few published studies on its nonhuman primate (NHP) species (Carlito syrichta, Macaca fascicularis, and Nycticebus menagensis). Understanding the nature and extent of existing NHP research in the Philippines would help us to determine research gaps and opportunities. Thus, studies on NHPs of Philippine origin were retrieved from online databases including Google Scholar, ResearchGate, Primate Lit, the Integrated Taxonomic Information System (ITIS), and relevant authors. Of a total of 107 studies on Philippine NHP from 1989 to 2019, most focused on Philippine long-tailed macaque (52.78%), and Philippine tarsier (42.59%), with far fewer on Philippine slow loris (4.63%). Fewer studies were based on fieldwork (26.17%); more were based on research on captive animals or that undertaken in the laboratory (34.58%), or used only existing specimens or data (33.64%). Studies on wild NHPs were mostly conducted in Bohol Island. In terms of the type of research, studies on diseases (38.60%) were the most prevalent for macaques; ecological studies (23.91%) for tarsiers; evolutionary, taxonomic, and systematic studies (40%), and ecological and natural history studies (40%) for lorises. Philippine researchers and collaborations with foreign researchers contributed fewer studies (43.93%) than foreign-only researchers (56.07%). In conclusion, although research on Philippine NHPs is increasing, there is a bias with regards to the species studied, the type of research, and study location. Conservation-driven studies are also lacking. These gaps offer opportunities for further research to identify threats to the survival of NHPs in the Philippines, and for the development of plans and effective strategies for their conservation.

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.

Improving cost-efficiency of faecal genotyping: New tools for elephant species

Despite the critical need for non-invasive tools to improve monitoring of wildlife populations, especially for endangered and elusive species, faecal genetic sampling has not been adopted as regular practice, largely because of the associated technical challenges and cost. Substantial work needs to be undertaken to refine sample collection and preparation methods in order to improve sample set quality and provide cost-efficient tools that can effectively support wildlife management. In this study, we collected an extensive set of forest elephant (Loxodonta cyclotis) faecal samples throughout Gabon, Central Africa, and prepared them for genotyping using 107 single-nucleotide polymorphism assays. We developed a new quantitative polymerase chain reaction (PCR) assay targeting a 130-bp nuclear DNA fragment and demonstrated its suitability for degraded samples in all three elephant species. Using this assay to compare the efficacy of two sampling methods for faecal DNA recovery, we found that sampling the whole surface of a dung pile with a swab stored in a small tube of lysis buffer was a convenient method producing high extraction success and DNA yield. We modelled the influence of faecal quality and storage time on DNA concentration in order to provide recommendations for optimized collection and storage. The maximum storage time to ensure 75% success was two months for samples collected within 24 hours after defecation and extended to four months for samples collected within one hour. Lastly, the real-time quantitative PCR assay allowed us to predict genotyping success and pre-screen DNA samples, thus further increasing the cost-efficiency of our approach. We recommend combining the validation of an efficient sampling method, the build of in-country DNA extraction capacity for reduced storage time and the development of species-specific quantitative PCR assays in order to increase the cost-efficiency of routine non-invasive DNA analyses and expand the use of next-generation markers to non-invasive samples.

CHIIMP: An automated high-throughput microsatellite genotyping platform reveals greater allelic diversity in wild chimpanzees

Short tandem repeats (STRs), also known as microsatellites, are commonly used to noninvasively genotype wild-living endangered species, including African apes. Until recently, capillary electrophoresis has been the method of choice to determine the length of polymorphic STR loci. However, this technique is labor intensive, difficult to compare across platforms, and notoriously imprecise. Here we developed a MiSeq-based approach and tested its performance using previously genotyped fecal samples from long-term studied chimpanzees in Gombe National Park, Tanzania. Using data from eight microsatellite loci as a reference, we designed a bioinformatics platform that converts raw MiSeq reads into locus-specific files and automatically calls alleles after filtering stutter sequences and other PCR artifacts. Applying this method to the entire Gombe population, we confirmed previously reported genotypes, but also identified 31 new alleles that had been missed due to sequence differences and size homoplasy. The new genotypes, which increased the allelic diversity and heterozygosity in Gombe by 61% and 8%, respectively, were validated by replicate amplification and pedigree analyses. This demonstrated inheritance and resolved one case of an ambiguous paternity. Using both singleplex and multiplex locus amplification, we also genotyped fecal samples from chimpanzees in the Greater Mahale Ecosystem in Tanzania, demonstrating the utility of the MiSeq-based approach for genotyping nonhabituated populations and performing comparative analyses across field sites. The new automated high-throughput analysis platform (available at https://github.com/ShawHahnLab/chiimp) will allow biologists to more accurately and effectively determine wildlife population size and structure, and thus obtain information critical for conservation efforts.