The Promise and Practicality of Population Genomics Research with Endangered Species

Modern, archaeological, and paleontological DNA analysis of a human-harvested marine gastropod (Strombus pugilis) from Caribbean Panama

Although protocols exist for the recovery of ancient DNA from land snail and marine bivalve shells, marine conch shells have yet to be studied from a paleogenomic perspective. We first present reference assemblies for both a 623.7 Mbp nuclear genome and a 15.4 kbp mitochondrial genome for Strombus pugilis, the West Indian fighting conch. We next detail a method to extract and sequence DNA from conch shells and apply it to conch from Bocas del Toro, Panama across three time periods: recently-eaten and discarded (n = 3), Late Holocene (984-1258 before present [BP]) archaeological midden (n = 5), and mid-Holocene (5711-7187 BP) paleontological fossil coral reef (n = 5). These results are compared to control DNA extracted from live-caught tissue and fresh shells (n = 5). Using high-throughput sequencing, we were able to obtain S. pugilis nuclear sequence reads from shells across all age periods: up to 92.5 thousand filtered reads per sample in live-caught shell material, 4.57 thousand for modern discarded shells, 12.1 thousand reads for archaeological shells, and 114 reads in paleontological shells. We confirmed authenticity of the ancient DNA recovered from the archaeological and paleontological shells based on 5.7× higher average frequency of deamination-driven misincorporations and 15% shorter average read lengths compared to the modern shells. Reads also mapped to the S. pugilis mitochondrial genome for all but the paleontological shells, with consistent ratios of mitochondrial to nuclear mapped reads across sample types. Our methods can be applied to diverse archaeological sites to facilitate reconstructions of the long-term impacts of human behaviour on mollusc evolutionary biology.

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.

Non-invasive genotyping with a massively parallel sequencing panel for the detection of SNPs in HPA-axis genes

We designed a genotyping panel for the investigation of the genetic underpinnings of inter-individual differences in aggression and the physiological stress response. The panel builds on single nucleotide polymorphisms (SNPs) in genes involved in the three subsystems of the hypothalamic-pituitary-adrenal (HPA)-axis: the catecholamine, serotonin and corticoid metabolism. To promote the pipeline for use with wild animal populations, we used non-invasively collected faecal samples from a wild population of Assamese macaques (Macaca assamensis). We targeted loci of 46 previously reported SNPs in 21 candidate genes coding for elements of the HPA-axis and amplified and sequenced them using next-generation Illumina sequencing technology. We compared multiple bioinformatics pipelines for variant calling and variant effect prediction. Based on this strategy and the application of different quality thresholds, we identified up to 159 SNPs with different types of predicted functional effects among our natural study population. This study provides a massively parallel sequencing panel that will facilitate integrating large-scale SNP data into behavioural and physiological studies. Such a multi-faceted approach will promote understanding of flexibility and constraints of animal behaviour and hormone physiology.

Emerging and Enduring Issues in Primate Conservation Genetics

Conservation genetics is a branch of conservation biology that uses molecular data to assist in the conservation and management of imperiled populations, subspecies, and species. In this review, I examine conservation action plans (CAPs)—instrumental documents designed to influence conservation policy—for selected primate species. I use the information contained in CAPs as a means to guide this review. The primary genetics-based topics that are mentioned in CAPs are genetic connectivity, inbreeding, and subspecies/species delimitation. I discuss these topics as well as historical demographic inference and hybridization using examples from wild primate species to illustrate the myriad ways in which genetics can assist in conservation efforts. I also discuss some recent technological advances such as genomic capture techniques and the potential to do molecular work in remote locations.

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.

Is genetic drift to blame for testicular dysgenesis syndrome in Semliki chimpanzees (Pan troglodytes schweinfurthii)?

BACKGROUND

We present 3 likely cases of testicular dysgenesis syndrome (TDS) within a community of chimpanzees (Pan troglodytes schweinfurthii). We tested whether genetic drift may be the culprit, as a genetic cause has been suspected to account for TDS among other wildlife.

METHODS

We successfully sequenced a 367-bp segment spanning the first hypervariable region within the D-loop of the mitochondrial genome for 78 DNA samples.

RESULTS

We found 24 polymorphic sequence sites consisting of 7 singletons and 17 parsimony informative sites. This sample contained 9 haplotypes with a diversity index of 0.78 (SD = 0.03). All tests against the null hypothesis of neutral polymorphisms were non-significant (P > .10). The mismatch distribution of pairwise differences does not fit a Poisson's curve (raggedness index = 0.166; SSD = 0.12; P = 1).

CONCLUSIONS

Thus, we found no significant signs of genetic isolation, population expansion, or genetic bottleneck. Alternative causes of TDS and how they might pertain to this population are discussed.

A massively parallel strategy for STR marker development, capture, and genotyping

Short tandem repeat (STR) variants are highly polymorphic markers that facilitate powerful population genetic analyses. STRs are especially valuable in conservation and ecological genetic research, yielding detailed information on population structure and short-term demographic fluctuations. Massively parallel sequencing has not previously been leveraged for scalable, efficient STR recovery. Here, we present a pipeline for developing STR markers directly from high-throughput shotgun sequencing data without a reference genome, and an approach for highly parallel target STR recovery. We employed our approach to capture a panel of 5000 STRs from a test group of diademed sifakas (Propithecus diadema, n = 3), endangered Malagasy rainforest lemurs, and we report extremely efficient recovery of targeted loci—97.3–99.6% of STRs characterized with ≥10x non-redundant sequence coverage. We then tested our STR capture strategy on P. diadema fecal DNA, and report robust initial results and suggestions for future implementations. In addition to STR targets, this approach also generates large, genome-wide single nucleotide polymorphism (SNP) panels from flanking regions. Our method provides a cost-effective and scalable solution for rapid recovery of large STR and SNP datasets in any species without needing a reference genome, and can be used even with suboptimal DNA more easily acquired in conservation and ecological studies.

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.