Using genetics to understand the dynamics of wild primate populations

Ethnoprimatology reveals new extended distribution of critically endangered banded langur Presbytis femoralis (Martin, 1838) in Pahang, Malaysia: Insights from indigenous traditional knowledge and molecular analysis

The banded langur (Presbytis femoralis) is a critically endangered primate, which within Malaysia has not been known to extend significantly outside the state of Johor. Traditionally, distribution studies on this highly threatened primate have relied on conventional methods such as DNA identification, live counting, and camera trapping. However, ethnoprimatology offers an alternative approach to data collection, involving the active participation of indigenous and local communities possessing valuable knowledge and experience with local primate species. This study employed an integrated approach incorporating ethnoprimatology by utilizing pooled local expert opinion, local surveys, interviews, and fecal DNA analysis, resulting in a novel distribution range for the banded langur. The combination of expert opinions revealed this species' most optimistic distribution scenario across Johor and Pahang, inhabiting various ecosystems, including lowland forests, peat swamps, and human-modified landscapes. Further interviews and surveys conducted within the Orang Asli community in Tasik Chini and Tasek Bera have provided additional support for the revised distribution, documenting occurrences of banded langur utilization in indigenous practices, such as food consumption, cultural beliefs, medicinal applications, and craftsmanship. Phylogenetic analysis demonstrated genetic differentiation between populations in Johor and Pahang, with the populations in the southern part of Peninsular Malaysia likely serving as ancestral sources for other populations. Consequently, this study not only elucidated the updated distribution of banded langur through DNA records and direct observations but also established the efficacy of ethnoprimatology as a precursory tool for uncovering the present distribution patterns of other primate species in Malaysia.

Sequencing the orthologs of human autosomal forensic short tandem repeats provides individual- and species-level identification in African great apes

BACKGROUND

Great apes are a global conservation concern, with anthropogenic pressures threatening their survival. Genetic analysis can be used to assess the effects of reduced population sizes and the effectiveness of conservation measures. In humans, autosomal short tandem repeats (aSTRs) are widely used in population genetics and for forensic individual identification and kinship testing. Traditionally, genotyping is length-based via capillary electrophoresis (CE), but there is an increasing move to direct analysis by massively parallel sequencing (MPS). An example is the ForenSeq DNA Signature Prep Kit, which amplifies multiple loci including 27 aSTRs, prior to sequencing via Illumina technology. Here we assess the applicability of this human-based kit in African great apes. We ask whether cross-species genotyping of the orthologs of these loci can provide both individual and (sub)species identification.

RESULTS

The ForenSeq kit was used to amplify and sequence aSTRs in 52 individuals (14 chimpanzees; 4 bonobos; 16 western lowland, 6 eastern lowland, and 12 mountain gorillas). The orthologs of 24/27 human aSTRs amplified across species, and a core set of thirteen loci could be genotyped in all individuals. Genotypes were individually and (sub)species identifying. Both allelic diversity and the power to discriminate (sub)species were greater when considering STR sequences rather than allele lengths. Comparing human and African great-ape STR sequences with an orangutan outgroup showed general conservation of repeat types and allele size ranges. Variation in repeat array structures and a weak relationship with the known phylogeny suggests stochastic origins of mutations giving rise to diverse imperfect repeat arrays. Interruptions within long repeat arrays in African great apes do not appear to reduce allelic diversity.

CONCLUSIONS

Orthologs of most human aSTRs in the ForenSeq DNA Signature Prep Kit can be analysed in African great apes. Primer redesign would reduce observed variability in amplification across some loci. MPS of the orthologs of human loci provides better resolution for both individual and (sub)species identification in great apes than standard CE-based approaches, and has the further advantage that there is no need to limit the number and size ranges of analysed loci.

Population dynamics and genetic connectivity in recent chimpanzee history

Knowledge on the population history of endangered species is critical for conservation, but whole-genome data on chimpanzees (Pan troglodytes) is geographically sparse. Here, we produced the first non-invasive geolocalized catalog of genomic diversity by capturing chromosome 21 from 828 non-invasive samples collected at 48 sampling sites across Africa. The four recognized subspecies show clear genetic differentiation correlating with known barriers, while previously undescribed genetic exchange suggests that these have been permeable on a local scale. We obtained a detailed reconstruction of population stratification and fine-scale patterns of isolation, migration, and connectivity, including a comprehensive picture of admixture with bonobos (Pan paniscus). Unlike humans, chimpanzees did not experience extended episodes of long-distance migrations, which might have limited cultural transmission. Finally, based on local rare variation, we implement a fine-grained geolocalization approach demonstrating improved precision in determining the origin of confiscated chimpanzees.

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.

Maximizing the acquisition of unique reads in noninvasive capture sequencing experiments

Noninvasive samples as a source of DNA are gaining interest in genomic studies of endangered species. However, their complex nature and low endogenous DNA content hamper the recovery of good quality data. Target capture has become a productive method to enrich the endogenous fraction of noninvasive samples, such as faeces, but its sensitivity has not yet been extensively studied. Coping with faecal samples with an endogenous DNA content below 1% is a common problem when prior selection of samples from a large collection is not possible. However, samples classified as unfavourable for target capture sequencing might be the only representatives of unique specific geographical locations, or to answer the question of interest. To explore how library complexity may be increased without repeating DNA extractions and generating new libraries, in this study we captured the exome of 60 chimpanzees (Pan troglodytes) using faecal samples with very low proportions of endogenous content (<1%). Our results indicate that by performing additional hybridizations of the same libraries, the molecular complexity can be maintained to achieve higher coverage. Also, whenever possible, the starting DNA material for capture should be increased. Finally, we specifically calculated the sequencing effort needed to avoid exhausting the library complexity of enriched faecal samples with low endogenous DNA content. This study provides guidelines, schemes and tools for laboratories facing the challenges of working with noninvasive samples containing extremely low amounts of endogenous DNA.

Influence of habitat conditions on group size, social organization, and birth pattern of golden langur (Trachypithecus geei)

We studied endangered golden langurs in fragmented and altered habitats to understand the consequences of habitat conditions on group size, social organization, and birth seasonality. We selected 12 groups inhabiting forest edge and forest core of Chakrashila Wildlife Sanctuary (henceforth Chakrashila WLS) and adjoining the Abhaya rubber plantation. Each group was monitored every month from May 2013 to September 2016 and recorded the age-sex of individuals in the group. The births were recorded with the individual identity of females in five focal groups. The overall group size of golden langur was 11.3 ± 3.5_SD, and ranged between 5 and 18. The mean group size in forest core, forest edge, and rubber plantation differed significantly. We recorded a total of 46 births in 12 groups across the three different habitats. The number of infants correlates positively with adult females and group size across all the 12 groups for all the years. The number of births that occurred in all the months varied significantly across the months. Births occurred in all the months but peaked between May and September (82.6%). The mean number of births positively correlated with mean monthly rainfall. Mean inter-birth interval was 24.5 ± 1.6_SD months that did not vary between the females. It therefore appears that group size is sensitive to forest type, and births are positively related to social and environmental factors. The behavioral parameters may influence life-history traits if continuous habitat alteration persists.

Spatial and temporal distribution of Yunnan snub-nosed monkey, Rhinopithecus bieti, indices

Studying elusive species of conservation concern might be difficult for technical and ethical reasons. However, censuses can be based on the observation of activity indices. When coupled to non-invasive genetic methods this approach can provide extremely precise information about population size, individual movements and diseases. However, the design of optimal sampling is dependent on a knowledge on group distribution and possible variations of detectability of index targets. The aim of this study was to document the distribution of Yunnan snub-nosed monkey indices in space and time in that perspective. Based on transects carried out across the range of a fed population and on counts along the trail across the range of a wild group, we show that 2–3 day stays of a group in a place of some hectares were sufficient to get an homogeneous distribution of indices. Furthermore, the number of indices found were dependent on both pig presence and season. On the other hand, on a large scale of 100 km2 indices were spatially distributed as nested clusters. Indices distribution indicated a strong preference towards southern slopes and altitudes ranging between 2900 and 3400 m. Those observations pinpoint the importance of considering spatial scale to organise sampling designed to estimate population distribution.

A RAD-sequencing approach to genome-wide marker discovery, genotyping, and phylogenetic inference in a diverse radiation of primates

Until recently, most phylogenetic and population genetics studies of nonhuman primates have relied on mitochondrial DNA and/or a small number of nuclear DNA markers, which can limit our understanding of primate evolutionary and population history. Here, we describe a cost-effective reduced representation method (ddRAD-seq) for identifying and genotyping large numbers of SNP loci for taxa from across the New World monkeys, a diverse radiation of primates that shared a common ancestor ~20-26 mya. We also estimate, for the first time, the phylogenetic relationships among 15 of the 22 currently-recognized genera of New World monkeys using ddRAD-seq SNP data using both maximum likelihood and quartet-based coalescent methods. Our phylogenetic analyses robustly reconstructed three monophyletic clades corresponding to the three families of extant platyrrhines (Atelidae, Pitheciidae and Cebidae), with Pitheciidae as basal within the radiation. At the genus level, our results conformed well with previous phylogenetic studies and provide additional information relevant to the problematic position of the owl monkey (Aotus) within the family Cebidae, suggesting a need for further exploration of incomplete lineage sorting and other explanations for phylogenetic discordance, including introgression. Our study additionally provides one of the first applications of next-generation sequencing methods to the inference of phylogenetic history across an old, diverse radiation of mammals and highlights the broad promise and utility of ddRAD-seq data for molecular primatology.

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.

Genetic structure in the southernmost populations of black-and-gold howler monkeys (Alouatta caraya) and its conservation implications

Black-and-gold howler monkeys Alouatta caraya, are arboreal primates, inhabitants of Neotropical forests, highly susceptible to the yellow fever virus, considered early 'sentinels' of outbreaks, and thus, of major epidemiological importance. Currently, anthropogenic habitat loss and modifications threatens their survival. Habitat modification can prevent, reduce or change dispersal behavior, which, in turn, may influence patterns of gene flow. We explored past and contemporary levels of genetic diversity, elucidated genetic structure and identified its possible drivers, in ten populations (n = 138) located in the southernmost distribution range of the species in South America, in Argentina and Paraguay. Overall, genetic variability was moderate (ten microsatellites: 3.16 ± 0.18 alleles per locus, allelic richness of 2.93 ± 0.81, 0.443±0.025 unbiased expected heterozygosity; 22 haplotypes of 491-bp mitochondrial Control Region, haplotypic diversity of 0.930 ± 0.11, and nucleotide diversity of0.01± 0.007). Significant evidence of inbreeding was found in a population that was, later, decimated by yellow fever. Population-based gene flow measures (FST = 0.13; θST = 018), hierarchical analysis of molecular variance and Bayesian clustering methods revealed significant genetic structure, grouping individuals into four clusters. Shared haplotypes and lack of mitochondrial differentiation (non-significant θST) among some populations seem to support the hypothesis of historical dispersal via riparian forests. Current resistance analyses revealed a significant role of landscape features in modeling contemporary gene flow: continuous forest and riparian forests could promote genetic exchange, whereas disturbed forests or crop/grassland fields may restrict it. Estimates of effective population size allow anticipating that the studied populations will lose 75% of heterozygosity in less than 50 generations. Our findings suggest that anthropogenic modifications on native forests, increasingly ongoing in Northeastern Argentina, Southern Paraguay and Southeastern Brazil, might prevent the dispersal of howlers, leading to population isolation. To ensure long-term viability and maintain genetic connectivity of A. caraya remnant populations, we recommend preserving and restoring habitat continuity. To conserve the species genetic pool, as well, the four genetic clusters identified here should be considered separate Management Units and given high conservation priority. In light of our findings and considering complementary non-genetic information, we suggest upgrading the international conservation status of A. caraya to "Vulnerable".

Temporal and geographic patterns of kinship structure in common dolphins (Delphinus delphis) suggest site fidelity and female-biased long-distance dispersal

ABSTRACT

Social structure plays a crucial role in determining a species' dispersal patterns and genetic structure. Cetaceans show a diversity of social and mating systems, but their effects on dispersal and genetic structure are not well known, in part because of technical difficulties in obtaining robust observational data. Here, we combine genetic profiling and GIS analysis to identify patterns of kin distribution over time and space, to infer mating structure and dispersal patterns in short-beaked common dolphins (Delphinus delphis). This species is highly social, and exhibits weak spatial genetic structure in the Northeast Atlantic and Mediterranean Sea, thought to result from fluid social structure and low levels of site fidelity. We found that although sampled groups were not composed of closely related individuals, close kin were frequently found in the same geographic location over several years. Our results suggest that common dolphin exhibits some level of site fidelity, which could be explained by foraging for temporally varying prey resource in areas familiar to individuals. Dispersal from natal area likely involves long-distance movements of females, as males are found more frequently than females in the same locations as their close kin. Long-distance dispersal may explain the near panmixia observed in this species. By analysing individuals sampled in the same geographic location over multiple years, we avoid caveats associated with divergence-based methods of inferring sex-biased dispersal. We thus provide a unique perspective on this species' social structure and dispersal behaviour, and how it relates to the observed low levels of population genetic structure in European waters.

SIGNIFICANCE STATEMENT

Movement patterns and social interactions are aspects of wild animal's behaviour important for understanding their ecology. However, tracking these behaviours directly can be very challenging in wide-ranging species such as whales and dolphins. In this study, we used genetic information to detect how patterns of kin associations change in space and time, to infer aspects of movement and social structure. We identified previously unknown site fidelity, and suggested that dispersal usually involves females, travelling long distances from the natal area. Our data analysis strategy overcomes known limitations of previously used genetic inference methods, and provides a new approach to identify differences in dispersal between the sexes, which contribute to better understanding of the species' behaviour and ecology. In this case, we suggest that females are more likely to disperse than males, a pattern unusual amongst mammals.

Cost-effective scat-detection dogs: unleashing a powerful new tool for international mammalian conservation biology

Recently, detection dogs have been utilized to collect fecal samples from cryptic and rare mammals. Despite the great promise of this technique for conservation biology, its broader application has been limited by the high cost (tens to hundreds of thousands of dollars) and logistical challenges of employing a scat-detection dog team while conducting international, collaborative research. Through an international collaboration of primatologists and the Chinese Ministry of Public Security, we trained and used a detection dog to find scat from three species of unhabituated, free-ranging primates, for less than $3,000. We collected 137 non-human primate fecal samples that we confirmed by sequencing taxonomically informative genetic markers. Our detection dog team had a 92% accuracy rate, significantly outperforming our human-only team. Our results demonstrate that detection dogs can locate fecal samples from unhabituated primates with variable diets, locomotion, and grouping patterns, despite challenging field conditions. We provide a model for in-country training, while also building local capacity for conservation and genetic monitoring. Unlike previous efforts, our approach will allow for the wide adoption of scat-detection dogs in international conservation biology.

High Rate of Simian Immunodeficiency Virus (SIV) Infections in Wild Chimpanzees in Northeastern Gabon

The emergence of HIV-1 groups M, N, O, and P is the result of four independent cross-species transmissions between chimpanzees (cpz) and gorillas (gor) from central/south Cameroon and humans respectively. Although the first two SIVcpz were identified in wild-born captive chimpanzees in Gabon in 1989, no study has been conducted so far in wild chimpanzees in Gabon. To document the SIVcpz infection rate, genetic diversity, and routes of virus transmission, we analyzed 1458 faecal samples collected in 16 different locations across the country, and we conducted follow-up missions in two of them. We found 380 SIV antibody positive samples in 6 different locations in the north and northeast. We determined the number of individuals collected by microsatellite analysis and obtained an adjusted SIV prevalence of 39.45%. We performed parental analysis to investigate viral spread between and within communities and found that SIVs were epidemiologically linked and were transmitted by both horizontal and vertical routes. We amplified pol and gp41 fragments and obtained 57 new SIVcpzPtt strains from three sites. All strains, but one, clustered together within a specific phylogeographic clade. Given that these SIV positive samples have been collected nearby villages and that humans continue to encroach in ape's territories, the emergence of a new HIV in this area needs to be considered.

An invertebrate stomach's view on vertebrate ecology: certain invertebrates could be used as "vertebrate samplers" and deliver DNA-based information on many aspects of vertebrate ecology

Recent studies suggest that vertebrate genetic material ingested by invertebrates (iDNA) can be used to investigate vertebrate ecology. Given the ubiquity of invertebrates that feed on vertebrates across the globe, iDNA might qualify as a very powerful tool for 21st century population and conservation biologists. Here, we identify some invertebrate characteristics that will likely influence iDNA retrieval and elaborate on the potential uses of invertebrate-derived information. We hypothesize that beyond inventorying local faunal diversity, iDNA should allow for more profound insights into wildlife population density, size, mortality, and infectious agents. Based on the similarities of iDNA with other low-quality sources of DNA, a general technical framework for iDNA analyses is proposed. As it is likely that no such thing as a single ideal iDNA sampler exists, forthcoming research efforts should aim at cataloguing invertebrate properties relevant to iDNA retrieval so as to guide future usage of the invertebrate tool box.

When bugs reveal biodiversity

One of the fundamental challenges of conservation biology is gathering data on species distribution and abundance. And unless conservationists know where a species is found and in which numbers, it is very difficult to apply effective conservation efforts. In today's age of increasingly powerful monitoring tools, instant communication and online databases, one might be forgiven for thinking that such knowledge is easy to come by. However, of the approximately 5,400 terrestrial mammals on the IUCN Red List, no fewer than 789 (ca. 14%) are listed as 'Data Deficient' (IUCN 2012) – IUCN’s term for 'haven't got a clue'. Until recently, the only way to gather information of numbers and distribution of terrestrial mammals (and many other vertebrates) was through observational-based approaches such as visual records, the presence of tracks or spoor or even identification from bushmeat or hunters' trophies pinned to the walls in local villages. While recent technological developments have considerably improved the efficacy of such approaches, for example, using remote-sensing devices such as audio- or camera-traps or even remote drones (Koh & Wich 2012), there has been a growing realization of the power of molecular methods that identify mammals based on trace evidence. Suitable substrates include the obvious, such as faecal and hair samples (e.g. Vigilant et al. 2009), to the less obvious, including environmental DNA extracted from sediments, soil or water samples (e.g. Taberlet et al. 2012), and as recently demonstrated, the dietary content of blood-sucking invertebrates (Gariepy et al. 2012; Schnell et al. 2012). In this issue of Molecular Ecology, Calvignac-Spencer et al. (2013) present a potentially powerful development in this regard; diet analysis of carrion flies. With their near global distribution, and as most field biologists know, irritatingly high frequency in most terrestrial areas of conservation concern (which directly translates into ease of sampling them), the authors present extremely encouraging results that indicate how carnivorous flies may soon represent a strong weapon in the conservation arsenal.

Carrion fly-derived DNA as a tool for comprehensive and cost-effective assessment of mammalian biodiversity

Large-scale monitoring schemes are essential in assessing global mammalian biodiversity, and in this framework, leeches have recently been promoted as an indirect source of DNA from terrestrial mammal species. Carrion feeding flies are ubiquitous and can be expected to feed on many vertebrate carcasses. Hence, we tested whether fly-derived DNA analysis may also serve as a novel tool for mammalian diversity surveys. We screened DNA extracted from 201 carrion flies collected in tropical habitats of Côte d'Ivoire and Madagascar for mammal DNA using multiple PCR systems and retrieved DNA sequences from a diverse set of species (22 in Côte d'Ivoire, four in Madagascar) exploiting distinct forest strata and displaying a broad range of body sizes. Deep sequencing of amplicons generated from pools of flies performed equally well as individual sequencing approaches. We conclude that the analysis of fly-derived DNA can be implemented in a very rapid and cost-effective manner and will give a relatively unbiased picture of local mammal diversity. Carrion flies therefore represent an extraordinary and thus far unexploited resource of mammal DNA, which will probably prove useful for future inventories of wild mammal communities.

Genetic diversity and population history of a critically endangered primate, the northern muriqui (Brachyteles hypoxanthus)

Social, ecological, and historical processes affect the genetic structure of primate populations, and therefore have key implications for the conservation of endangered species. The northern muriqui (Brachyteles hypoxanthus) is a critically endangered New World monkey and a flagship species for the conservation of the Atlantic Forest hotspot. Yet, like other neotropical primates, little is known about its population history and the genetic structure of remnant populations. We analyzed the mitochondrial DNA control region of 152 northern muriquis, or 17.6% of the 864 northern muriquis from 8 of the 12 known extant populations and found no evidence of phylogeographic partitions or past population shrinkage/expansion. Bayesian and classic analyses show that this finding may be attributed to the joint contribution of female-biased dispersal, demographic stability, and a relatively large historic population size. Past population stability is consistent with a central Atlantic Forest Pleistocene refuge. In addition, the best scenario supported by an Approximate Bayesian Computation analysis, significant fixation indices (Φ(ST) = 0.49, Φ(CT) = 0.24), and population-specific haplotypes, coupled with the extirpation of intermediate populations, are indicative of a recent geographic structuring of genetic diversity during the Holocene. Genetic diversity is higher in populations living in larger areas (>2,000 hectares), but it is remarkably low in the species overall (θ = 0.018). Three populations occurring in protected reserves and one fragmented population inhabiting private lands harbor 22 out of 23 haplotypes, most of which are population-exclusive, and therefore represent patchy repositories of the species' genetic diversity. We suggest that these populations be treated as discrete units for conservation management purposes.

An ecological and conservation perspective on advances in the applied virology of zoonoses

The aim of this manuscript is to describe how modern advances in our knowledge of viruses and viral evolution can be applied to the fields of disease ecology and conservation. We review recent progress in virology and provide examples of how it is informing both empirical research in field ecology and applied conservation. We include a discussion of needed breakthroughs and ways to bridge communication gaps between the field and the lab. In an effort to foster this interdisciplinary effort, we have also included a table that lists the definitions of key terms. The importance of understanding the dynamics of zoonotic pathogens in their reservoir hosts is emphasized as a tool to both assess risk factors for spillover and to test hypotheses related to treatment and/or intervention strategies. In conclusion, we highlight the need for smart surveillance, viral discovery efforts and predictive modeling. A shift towards a predictive approach is necessary in today's globalized society because, as the 2009 H1N1 pandemic demonstrated, identification post-emergence is often too late to prevent global spread. Integrating molecular virology and ecological techniques will allow for earlier recognition of potentially dangerous pathogens, ideally before they jump from wildlife reservoirs into human or livestock populations and cause serious public health or conservation issues.

Noninvasive genome sampling in chimpanzees

The inevitable has happened: genomic technologies have been added to our noninvasive genetic sampling repertoire. In this issue of Molecular Ecology, Perry et al. (2010) demonstrate how DNA extraction from chimpanzee faeces, followed by a series of steps to enrich for target loci, can be coupled with next-generation sequencing. These authors collected sequence and single-nucleotide polymorphism (SNP) data at more than 600 genomic loci (chromosome 21 and the X) and the complete mitochondrial DNA. By design, each locus was 'deep sequenced' to enable SNP identification. To demonstrate the reliability of their data, the work included samples from six captive chimps, which allowed for a comparison between presumably genuine SNPs obtained from blood and potentially flawed SNPs deduced from faeces. Thus, with this method, anyone with the resources, skills and ambition to do genome sequencing of wild, elusive, or protected mammals can enjoy all of the benefits of noninvasive sampling.

Genomic-scale capture and sequencing of endogenous DNA from feces

Genomic-level analyses of DNA from non-invasive sources would facilitate powerful conservation and evolutionary studies in natural populations of endangered and otherwise elusive species. However, the typical low quantity and poor quality of DNA that is extracted from non-invasive samples have generally precluded such work. Here we apply a modified DNA capture protocol that, when used in combination with massively-parallel sequencing technology, facilitates efficient and highly-accurate resequencing of megabases of specified nuclear genomic regions from fecal DNA samples. We validated our approach by comparing genetic variants identified from corresponding fecal and blood DNA samples of six western chimpanzees (Pan troglodytes verus) across more than 1.5 megabases of chromosome 21, chromosome X, and the complete mitochondrial genome. Our results suggest that it is now feasible to conduct genomic studies in natural populations for which constraints on invasive sampling have otherwise long been a barrier. The data we collected also provided an opportunity to examine western chimpanzee genetic diversity at unprecedented scale. Despite high mitochondrial genome diversity (π = 0.585%), western chimpanzees have a low ratio (0.42) of X chromosomal (π = 0.034%) to autosomal (chromosome 21 π = 0.081%) sequence diversity, a pattern that may reflect an unusual demographic history of this subspecies.