Comparison of DNA quantity and quality from fecal samples of mammals transported in ethanol or lysis buffer

Using fecal microbial community profiles through sequencing approaches helps to unravel the intimate interplay between health, wellness, and diet in wild animals with their environment. Ensuring the proper preservation of fecal samples before processing is crucial to ensure reliable results. In this study, we evaluated the efficiency of two different preservation methods, considering the following criteria: DNA yield, quality and integrity, and microbial community structure based on Oxford Nanopore amplicon sequencing of the V3-V4 region of bacterial 16S rRNA and protozoa 18S rRNA genes. Eighteen matched pairs of mammalian fecal samples were collected and transported in 99.8% ethanol and lysis buffer; processing occurred between 55 and 461 days post-collection. Wilcoxon signed-rank tests were used to analyze quantitative measurements for paired samples. The A260/280 ratio, a measure of nucleic acid purity, was assessed descriptively for each media, and the Bartlett test evaluated dispersion of this ratio. A Fisher test was performed to compare the number of positive reactions for DNA extraction or PCR amplification of the 16S and 18S rRNA genes between both media. The concentration of total DNA and amplicons, as well as the number of reads obtained in sequencing, was significantly higher in the samples preserved with lysis buffer compared to ethanol, with magnitudes up to three times higher. Electrophoretic analysis of total DNA and amplicons further confirmed superior DNA integrity in lysis buffer preserved samples. The A260/280 values  obtained using the lysis buffer were of optimal purity (mean: 1.92) and with little dispersion (SD: 0.27); on the other hand, the ethanol samples also presented an excellent average quality (mean: 1.94), but they were dispersed (SD: 1.10). For molecular studies using mammalian feces, the lysis buffer reagent proved to be a reliable solution for their collection, conservation, and storage.

A simple and effective method to enrich endogenous DNA from mammalian faeces

Utilization of faeces has long been a popular approach for genetic and ecological studies of wildlife. However, the success of molecular marker genotyping and genome resequencing is often unpredictable due to insufficient enrichment of endogenous DNA in the total faecal DNA that is dominated by bacterial DNA. Here, we report a simple and cheap method named PEERS to predominantly lyse animal cells over bacteria by using sodium dodecyl sulphate so as to discharge endogenous DNA into liquid phase before bacterial DNA. By brief centrifugation, total DNA with enriched endogenous fraction can be extracted from the supernatant using routine methods. Our assessments showed that the endogenous DNA extracted by PEERS was significantly enriched for various types of faeces from different species, preservation time and conditions. It significantly improves the genotyping correctness and efficiency of genome resequencing with the total additional cost of $ 0.1 and a short incubation step to treat a faecal sample. We also provide methods to assess the enrichment efficiency of mitochondrial and nuclear DNA and models to predict the usability of faecal DNA for genotyping of short tandem repeat, single-nucleotide polymorphism and whole-genome resequencing.

An experimental design for obtaining DNA of a target species and its diet from a single non-invasive genetic protocol

Next-generation sequencing technology has enabled accurate insights into the diet of wildlife species. The protocols for faecal sample collection and DNA extraction for diet analysis have differed from those focusing on target species, even in most studies combining questions on both aspects. We designed an experiment to evaluate two protocols using 11 parameters and select a single one that will generate both target species (Asiatic wild ass, Equus hemionus, in Israel) and diet DNA, as an effective strategy to minimise time, effort, and cost without hampering efficiency. In Protocol A, we swabbed the outer surface of faecal boluses and extracted DNA using a Stool Kit, while for Protocol B, we homogenised faecal matter from inside the bolus followed by extraction using a Powersoil Kit. Protocol A performed significantly better for four parameters, which included, for the target species, microsatellite amplification success and the quantity of the GAPDH gene; and for its diet, the number of exact sequence variants (ESVs) obtained at genus level and plant genus richness. However, there was no significant difference in the amplification success of sex-linked and plant markers, total reads at genus level, number of genera obtained and plant genus composition. Although we chose Protocol A, both protocols yielded results for the target species and its diet, demonstrating that one single protocol can be used for both purposes, although a pilot study is recommended to optimise the protocol for specific systems. This strategy may also be useful for studies combining target species and their gut microbiome and parasitic load.

Intestinal Parasites of Neotropical Wild Jaguars, Pumas, Ocelots, and Jaguarundis in Colombia: Old Friends Brought Back from Oblivion and New Insights

Neotropical wild felids (NWF) are obligate carnivore species present in Central and South America, and some are considered endangered due to constantly decreasing populations. NWF can become infected by a wide range of protozoan and metazoan parasites, some of them affecting their health conditions and others having anthropozoonotic relevance. Parasitological studies on NWF are still very scarce, and most data originated from dead or captive animals. On this account, the current study aimed to characterize gastrointestinal parasites of free-ranging jaguars (Panthera onca), pumas (Puma concolor), ocelots (Leopardus pardalis), and jaguarundis (Herpailurus yagouaroundi), i.e., four out of six NWF species endemic to Colombia. Fecal samples from jaguars (n = 10) and ocelots (n = 4) were collected between 2012 and 2017 as part of the Jaguar Corridor Initiative from six geographic locations in Colombia. In addition, cestode specimens were obtained during puma and jaguarundi necropsies. Scat samples were processed by standardized sodium acetate-acetic acid-formalin (SAF), sedimentation, and flotation techniques and by carbol fuchsin-stained fecal smears. Morphological evaluation of feces showed the presence of one cestode (Spirometra sp.), a nematode (Toxocara cati), an acanthocephalan (Oncicola sp.), and one cyst-forming coccidian (Cystoisospora-like oocysts). Feces oocysts were submitted to a Toxoplasma gondii-specific PCR for species identification, but no product was amplified. The cestodes isolated from a puma and jaguarundi were molecularly characterized by sequencing cytochrome c oxidase subunit I, identifying them as Taenia omissa and as a T. omissa sister lineage, respectively. These results collectively demonstrate the potential role of NWF as natural reservoir hosts for neglected zoonotic parasites (e.g., Spirometra sp., T. cati) and highlight their possible role in parasite transmission to human communities. Due to public health concerns, the occurrence of these parasites should be monitored in the future for appropriate zoonotic management practices in conservation strategies and wild felid health management programs.

Genetic diversity, population structure, and immigration, in a partially hunted puma population of south-central Argentina

Carnivores are decreasing globally due in part to anthropogenic ecological disturbances. In Argentina, human activities have fragmented wildlife habitat, thereby intensifying puma–livestock conflict and leading to population control of the predator species by hunting. We investigated genetic variability and population structure of pumas (Puma concolor) from three south-central Argentine provinces with two different management policies for the species: full protection versus legal hunting. All genetic estimates were based on 83 individuals genotyped at 25 species-specific microsatellite loci. The overall genetic diversity was high (observed heterozygosity = 0.63), but lower than in other South American populations. Spatial analyses revealed the presence of two bottlenecked genetic clusters with very similar diversity and low gene flow (3% per generation) between them. However, analyses based on a priori separated groups showed that gene flow follows increasing values of hunting pressure, converging to the area with the greatest number of individuals harvested. Our results suggest that hunting pressure likely is contributing to the gene flow pattern, limiting pumas’ movements and creating a metapopulation dynamic among geographic subpopulations. Integrated demographic and genetic approaches are needed to better understand pumas’ movements across the landscape and adopt successful management plans to achieve long-term population viability.

Does faecal matter reflect location? An initial assessment of isotopic variability between consumed prey remains and faecal matter for wild jaguars. ISOTOPES IN ENVIRONMENTAL AND HEALTH STUDIES 2019;55:478-498. [PMID: 31402710 DOI: 10.1080/10256016.2019.1648263]

Faecal isotopic analysis may complement other non-invasive wildlife survey tools for monitoring landscape use by carnivores, such as motion-detecting cameras and non-invasive genetic sampling. We analysed carbon, nitrogen, and strontium isotopes in faecal matter produced by jaguars (Panthera onca) as well as bones from consumed prey at the Mountain Pine Ridge Forest Reserve (MPR) in Belize, Central America. The MPR is ideally suited for a spatial isotope study as vegetation and geology both vary considerably. The isotopic composition of faecal matter should reflect the habitat and geology where consumed prey lived. We used bone from consumed prey recovered from jaguar scats as a proxy for diet. Faecal matter and bone showed comparable spatial isotopic trends, suggesting that the isotopic composition of jaguar faeces can be used to detect foraging in different habitats (pine forest versus broadleaf forest) or on different geologies (Mesozoic carbonates; Palaeozoic granite, contact metamorphics, and metasediments). This result is reassuring as bones are not always present in carnivore scats. Studying landscape use by cryptic and wide-ranging carnivore species like jaguars remains challenging. Isotopic analysis of faecal matter complements the existing array of non-invasive spatial monitoring tools.

Gut microbiota and their putative metabolic functions in fragmented Bengal tiger population of Nepal

Bengal tigers (Panthera tigris tigris) serve a pivotal role as an apex predator in forest ecosystems. To increase our knowledge on factors impacting the viability and health of this endangered species, we studied the gut microbiota in 32 individual Bengal tigers from three geographically separated areas (Chitwan National Park (CNP), Bardia National Park (BNP) and Suklaphanta Wildlife Reserve (SWR)) in Nepal, using noninvasive genetic sampling methods. Gut microbiota influence the immune system, impact various physiological functions, and modulates metabolic reactions, that ultimately impact the host health, behavior and development. Across the tiger populations in Nepal, we found significant differences in the composition of microbial communities based on their geographic locations. Specifically, we detected significant differences between CNP and the other two protected areas (CNP vs BNP: pseudo t = 1.944, P = 0.006; CNP vs SWR: pseudo t = 1.9942, P = 0.0071), but no differences between BNP and SWR. This mirrors what has been found for tiger gene flow in the same populations, suggesting gut microbiota composition and host gene flow may be linked. Furthermore, predictive metagenome functional content analysis (PICRUSt) revealed a higher functional enrichment and diversity for significant gut microbiota in the Chitwan tiger population and the lowest enrichment and diversity in Suklaphanta. The CNP tiger population contained higher proportions of microbiota that are associated with predicted functions relevant for metabolism of amino acid, lipid, xenobiotics biodegradation, terpenoides and polyketides than the SWR population. We conclude the tiger population structure, gut microbiota profile and associated functional metabolic categories are correlated, with geographically most separated CNP and SWR tiger population having the most distinct and different host genotype and microbiota profiles. Our work dramatically expands the understanding of tiger microbiota in wild populations and provides a valuable case study on how to investigate genetic diversity at different hierarchical levels, including hosts as well as their microbial communities.

Stream Transport and Retention of Environmental DNA Pulse Releases in Relation to Hydrogeomorphic Scaling Factors

The DNA of aquatic organisms can be identified in water sampled from freshwater ecosystems to detect species presence. Because these DNA-based methods (termed environmental DNA, eDNA) confirm species presence by proxy of DNA in water, the processes influencing eDNA transport and removal from water are critical to the method's efficacy and interpretation of results. Previous studies of aquatic eDNA transport and fate have employed uncontrolled field experiments, controlled studies in experimental streams, and laboratory column tests. As a step toward understanding the processes controlling eDNA transport and retention, we released and tracked experimental pulses of white sturgeon eDNA (novel to the system) in five fourth-order stream reaches with varied hydrology and geomorphology. We found strong support that stream water transient storage controls eDNA areal uptake rate (or spiraling length). We calculated the median spiraling length to be ∼260 m. Down channel slope correlated with transient storage, suggesting that this slope could be used as a proximate measure of eDNA removal into the benthic zone. Our results suggest that sampling effort should be increased in reaches with longer transient storage (or lower slopes) to compensate for the increase in eDNA retention.

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.

Genetic and genomic monitoring with minimally invasive sampling methods

The decreasing cost and increasing scope and power of emerging genomic technologies are reshaping the field of molecular ecology. However, many modern genomic approaches (e.g., RAD-seq) require large amounts of high-quality template DNA. This poses a problem for an active branch of conservation biology: genetic monitoring using minimally invasive sampling (MIS) methods. Without handling or even observing an animal, MIS methods (e.g., collection of hair, skin, faeces) can provide genetic information on individuals or populations. Such samples typically yield low-quality and/or quantities of DNA, restricting the type of molecular methods that can be used. Despite this limitation, genetic monitoring using MIS is an effective tool for estimating population demographic parameters and monitoring genetic diversity in natural populations. Genetic monitoring is likely to become more important in the future as many natural populations are undergoing anthropogenically driven declines, which are unlikely to abate without intensive adaptive management efforts that often include MIS approaches. Here, we profile the expanding suite of genomic methods and platforms compatible with producing genotypes from MIS, considering factors such as development costs and error rates. We evaluate how powerful new approaches will enhance our ability to investigate questions typically answered using genetic monitoring, such as estimating abundance, genetic structure and relatedness. As the field is in a period of unusually rapid transition, we also highlight the importance of legacy data sets and recommend how to address the challenges of moving between traditional and next-generation genetic monitoring platforms. Finally, we consider how genetic monitoring could move beyond genotypes in the future. For example, assessing microbiomes or epigenetic markers could provide a greater understanding of the relationship between individuals and their environment.

Assessment of genetic diversity, population structure, and gene flow of tigers (Panthera tigris tigris) across Nepal's Terai Arc Landscape

With fewer than 200 tigers (Panthera tigris tigris) left in Nepal, that are generally confined to five protected areas across the Terai Arc Landscape, genetic studies are needed to provide crucial information on diversity and connectivity for devising an effective country-wide tiger conservation strategy. As part of the Nepal Tiger Genome Project, we studied landscape change, genetic variation, population structure, and gene flow of tigers across the Terai Arc Landscape by conducting Nepal’s first comprehensive and systematic scat-based, non-invasive genetic survey. Of the 770 scat samples collected opportunistically from five protected areas and six presumed corridors, 412 were tiger (57%). Out of ten microsatellite loci, we retain eight markers that were used in identifying 78 individual tigers. We used this dataset to examine population structure, genetic variation, contemporary gene flow, and potential population bottlenecks of tigers in Nepal. We detected three genetic clusters consistent with three demographic sub-populations and found moderate levels of genetic variation (H_e = 0.61, A_R = 3.51) and genetic differentiation (F_ST = 0.14) across the landscape. We detected 3–7 migrants, confirming the potential for dispersal-mediated gene flow across the landscape. We found evidence of a bottleneck signature likely caused by large-scale land-use change documented in the last two centuries in the Terai forest. Securing tiger habitat including functional forest corridors is essential to enhance gene flow across the landscape and ensure long-term tiger survival. This requires cooperation among multiple stakeholders and careful conservation planning to prevent detrimental effects of anthropogenic activities on tigers.

The impact of endogenous content, replicates and pooling on genome capture from faecal samples

Target-capture approach has improved over the past years, proving to be very efficient tool for selectively sequencing genetic regions of interest. These methods have also allowed the use of noninvasive samples such as faeces (characterized by their low quantity and quality of endogenous DNA) to be used in conservation genomic, evolution and population genetic studies. Here we aim to test different protocols and strategies for exome capture using the Roche SeqCap EZ Developer kit (57.5 Mb). First, we captured a complex pool of DNA libraries. Second, we assessed the influence of using more than one faecal sample, extract and/or library from the same individual, to evaluate its effect on the molecular complexity of the experiment. We validated our experiments with 18 chimpanzee faecal samples collected from two field sites as a part of the Pan African Programme: The Cultured Chimpanzee. Those two field sites are in Kibale National Park, Uganda (N = 9) and Loango National Park, Gabon (N = 9). We demonstrate that at least 16 libraries can be pooled, target enriched through hybridization, and sequenced allowing for the genotyping of 951,949 exome markers for population genetic analyses. Further, we observe that molecule richness, and thus, data acquisition, increase when using multiple libraries from the same extract or multiple extracts from the same sample. Finally, repeated captures significantly decrease the proportion of off-target reads from 34.15% after one capture round to 7.83% after two capture rounds, supporting our conclusion that two rounds of target enrichment are advisable when using complex faecal samples.

Comparing hair-morphology and molecular methods to identify fecal samples from Neotropical felids

To avoid certain problems encountered with more-traditional and invasive methods in behavioral-ecology studies of mammalian predators, such as felids, molecular approaches have been employed to identify feces found in the field. However, this method requires a complete molecular biology laboratory, and usually also requires very fresh fecal samples to avoid DNA degradation. Both conditions are normally absent in the field. To address these difficulties, identification based on morphological characters (length, color, banding, scales and medullar patterns) of hairs found in feces could be employed as an alternative. In this study we constructed a morphological identification key for guard hairs of eight Neotropical felids (jaguar, oncilla, Geoffroy's cat, margay, ocelot, Pampas cat, puma and jaguarundi) and compared its efficiency to that of a molecular identification method, using the ATP6 region as a marker. For this molecular approach, we simulated some field conditions by postponing sample-conservation procedures. A blind test of the identification key obtained a nearly 70% overall success rate, which we considered equivalent to or better than the results of some molecular methods (probably due to DNA degradation) found in other studies. The jaguar, puma and jaguarundi could be unequivocally discriminated from any other Neotropical felid. On a scale ranging from inadequate to excellent, the key proved poor only for the margay, with only 30% of its hairs successfully identified using this key; and have intermediate success rates for the remaining species, the oncilla, Geoffroy's cat, ocelot and Pampas cat, were intermediate. Complementary information about the known distributions of felid populations may be necessary to substantially improve the results obtained with the key. Our own molecular results were even better, since all blind-tested samples were correctly identified. Part of these identifications were made from samples kept in suboptimal conditions, with some samples remaining outdoors for up to seven days, simulating conditions in the field. It appears that both methods can be used, depending on the available laboratory facilities and on the expected results.

Genetic Diversity and Population Structure of Mesoamerican Jaguars (Panthera onca): Implications for Conservation and Management

Mesoamerican jaguars (Panthera onca) have been extirpated from over 77% of their historic range, inhabiting fragmented landscapes at potentially reduced population sizes. Maintaining and restoring genetic diversity and connectivity across human-altered landscapes has become a major conservation priority; nonetheless large-scale genetic monitoring of natural populations is rare. This is the first regional conservation genetic study of jaguars to primarily use fecal samples collected in the wild across five Mesoamerican countries: Belize, Costa Rica, Guatemala, Honduras, and Mexico. We genotyped 445 jaguar fecal samples and examined patterns of genetic diversity and connectivity among 115 individual jaguars using data from 12 microsatellite loci. Overall, moderate levels of genetic variation were detected (NA = 4.50 ± 1.05, AR = 3.43 ± 0.22, HE = 0.59 ± 0.04), with Mexico having the lowest genetic diversity, followed by Honduras, Guatemala, Belize, and Costa Rica. Population-based gene flow measures (FST = 0.09 to 0.15, Dest = 0.09 to 0.21), principal component analysis, and Bayesian clustering applied in a hierarchical framework revealed significant genetic structure in Mesoamerican jaguars, roughly grouping individuals into four genetic clusters with varying levels of admixture. Gene flow was highest among Selva Maya jaguars (northern Guatemala and central Belize), whereas genetic differentiation among all other sampling sites was moderate. Genetic subdivision was most pronounced between Selva Maya and Honduran jaguars, suggesting limited jaguar movement between these close geographic regions and ultimately refuting the hypothesis of contemporary panmixia. To maintain a critical linkage for jaguars dispersing through the Mesoamerican landscape and ensure long-term viability of this near threatened species, we recommend continued management and maintenance of jaguar corridors. The baseline genetic data provided by this study underscores the importance of understanding levels of genetic diversity and connectivity to making informed management and conservation decisions with the goal to maintain functional connectivity across the region.

Non-invasive genetic sampling reveals diet shifts, but little difference in endoparasite richness and faecal glucocorticoids, in Belizean felids inside and outside protected areas

Many Neotropical felids are threatened with extinction due to direct effects of habitat destruction and/or human persecution. However, indirect and synergistic effects of human-felid conflict remain under-studied and potentially include increased stress and diet shifts that may negatively impact felid health. We hypothesized that faecal glucocorticoid metabolites (FGM) and endoparasite species richness (ESR) would be higher, and diet would shift, for felids outside protected areas where conflict occurs. In north-western Belize, a scat-detector dog located 336 faecal samples, identified to species and individual using DNA analyses. DNA amplification success was substantially higher within protected areas than outside. We detected jaguar, puma, ocelot, jaguarundi and domestic cat. FGMs were higher in puma and jaguarundi than in other felids, while ESR was similar across felids with domestic cats exhibiting the highest number of genera. Diet partitioning occurred among felids, but domestic cats may compete with ocelot and jaguarundi for small prey. Outside of protected areas, large cats shifted their diet to smaller prey and livestock remains were not found. Contrary to our hypotheses, FGM and ESR did not differ inside versus outside protected areas, but sample sizes were low in human-modified areas. We provide a baseline on wild felid adrenal activity, endoparasites and diet and suggest improvements to increase sample sizes outside protected areas. Our research provides a template for expanding non-invasive sampling approaches more widely across the range of Neotropical felids.

A Comparative Analysis of Genetic Diversity and Structure in Jaguars (Panthera onca), Pumas (Puma concolor), and Ocelots (Leopardus pardalis) in Fragmented Landscapes of a Critical Mesoamerican Linkage Zone

With increasing anthropogenic impact and landscape change, terrestrial carnivore populations are becoming more fragmented. Thus, it is crucial to genetically monitor wild carnivores and quantify changes in genetic diversity and gene flow in response to these threats. This study combined the use of scat detector dogs and molecular scatology to conduct the first genetic study on wild populations of multiple Neotropical felids coexisting across a fragmented landscape in Belize, Central America. We analyzed data from 14 polymorphic microsatellite loci in 1053 scat samples collected from wild jaguars (Panthera onca), pumas (Puma concolor), and ocelots (Leopardus pardalis). We assessed levels of genetic diversity, defined potential genetic clusters, and examined gene flow for the three target species on a countrywide scale using a combination of individual- and population-based analyses. Wild felids in Belize showed moderate levels of genetic variation, with jaguars having the lowest diversity estimates (HE = 0.57 ± 0.02; AR = 3.36 ± 0.09), followed by pumas (HE = 0.57 ± 0.08; AR = 4.20 ± 0.16), and ocelots (HE = 0.63 ± 0.03; AR = 4.16 ± 0.08). We observed low to moderate levels of genetic differentiation for all three target species, with jaguars showing the lowest degree of genetic subdivision across the country, followed by ocelots and pumas. Although levels of genetic diversity and gene flow were still fairly high, we detected evidence of fine-scale genetic subdivision, indicating that levels of genetic connectivity for wild felids in Belize are likely to decrease if habitat loss and fragmentation continue at the current rate. Our study demonstrates the value of understanding fine-scale patterns of gene flow in multiple co-occurring felid species of conservation concern, which is vital for wildlife movement corridor planning and prioritizing future conservation and management efforts within human-impacted landscapes.