When can noninvasive samples provide sufficient information in conservation genetics studies?

Non-invasive real-time genomic monitoring of the critically endangered kākāpō

We used non-invasive real-time genomic approaches to monitor one of the last surviving populations of the critically endangered kākāpō (Strigops habroptilus). We first established an environmental DNA metabarcoding protocol to identify the distribution of kākāpō and other vertebrate species in a highly localized manner using soil samples. Harnessing real-time nanopore sequencing and the high-quality kākāpō reference genome, we then extracted species-specific DNA from soil. We combined long read-based haplotype phasing with known individual genomic variation in the kākāpō population to identify the presence of individuals, and confirmed these genomically informed predictions through detailed metadata on kākāpō distributions. This study shows that individual identification is feasible through nanopore sequencing of environmental DNA, with important implications for future efforts in the application of genomics to the conservation of rare species, potentially expanding the application of real-time environmental DNA research from monitoring species distribution to inferring fitness parameters such as genomic diversity and inbreeding.

Genetic evidence of differential dispersal pattern in the Asiatic wild dog: Comparing two populations with different pack sizes

Introduction

Dispersal is a multi-causal, crucial life-history event in shaping the genetic and behavioral structure of mammals. We assessed the dispersal pattern of dholes aka Asiatic wild dog (Cuon alpinus), a social monogamous mammal at two tiger reserves of Maharashtra with different degrees of pack size and competition with tigers i.e. Tadoba-Andhari (TATR, smaller pack size, higher tiger density) and Nawegaon-Nagzira (NNTR, larger pack size, lower tiger density).

Methods

We used the microsatellite data of 174 individual genotypes (98 males and 67 females) to assess the dispersal pattern of dholes from two populations with varying pack size, tiger density, and landscape connectivity using gene flow as a proxy. We compared the population structure, pairwise F statistics, assignment index, and relatedness across a spatial scale.

Results and discussion

Overall, the results suggested a difference in sex-bias dispersal pattern for the two sub-populations, exhibiting significant results for female-biased dispersal in the TATR population with a smaller pack size and higher tiger density. Our study highlights the variability in sex-biased dispersal patterns in two different populations which could be the consequence of different variables such as pack size, tiger density, and geographical scale. The study warrants further quantitative investigation including several factors such as individual behavior, pack composition, pack size, tiger density, etc. In the present Anthropocene era, determining the sex bias in dispersal patterns for a short-range, pack-living carnivore will help in devising an effective conservation management plan for their long-term survival.

Genetic analyses reveal demographic decline and population differentiation in an endangered social carnivore, Asiatic wild dog

Deforestation and agricultural intensification have resulted in an alarming change in the global land cover over the past 300 years, posing a threat to species conservation. Dhole is a monophyletic, social canid and, being an endangered and highly forest-dependent species, is more prone to the loss of favorable habitat in the Anthropocene. We determined the genetic differentiation and demographic history of dhole across the tiger reserves of Maharashtra using the microsatellite data of 305 individuals. Simulation-based analyses revealed a 77-85% decline in the major dhole sub-populations. Protected areas have provided refuge to the historically declining dhole population resulting in clustering with strong genetic structure in the remnant dhole population. The historical population decline coincides with the extreme events in the landscape over the past 300 years. The study highlights the pattern of genetic differentiation and diversity of a highly forest-dependent species which can be associated with the loss of forest cover outside tiger reserves. It also warrants attention to develop conservation plans for the remnant surviving population of dholes in India.

Evaluation of the Minimum Sampling Design for Population Genomic and Microsatellite Studies: An Analysis Based on Wild Maize

Massive parallel sequencing (MPS) is revolutionizing the field of molecular ecology by allowing us to understand better the evolutionary history of populations and species, and to detect genomic regions that could be under selection. However, the economic and computational resources needed generate a tradeoff between the amount of loci that can be obtained and the number of populations or individuals that can be sequenced. In this work, we analyzed and compared two simulated genomic datasets fitting a hierarchical structure, two extensive empirical genomic datasets, and a dataset comprising microsatellite information. For all datasets, we generated different subsampling designs by changing the number of loci, individuals, populations, and individuals per population to test for deviations in classic population genetics parameters (H_S, F_IS, F_ST). For the empirical datasets we also analyzed the effect of sampling design on landscape genetic tests (isolation by distance and environment, central abundance hypothesis). We also tested the effect of sampling a different number of populations in the detection of outlier SNPs. We found that the microsatellite dataset is very sensitive to the number of individuals sampled when obtaining summary statistics. F_IS was particularly sensitive to a low sampling of individuals in the simulated, genomic, and microsatellite datasets. For the empirical and simulated genomic datasets, we found that as long as many populations are sampled, few individuals and loci are needed. For the empirical datasets, we found that increasing the number of populations sampled was important in obtaining precise landscape genetic estimates. Finally, we corroborated that outlier tests are sensitive to the number of populations sampled. We conclude by proposing different sampling designs depending on the objectives.

The curse of observer experience: Error in noninvasive genetic sampling

Noninvasive genetic sampling (NGS) is commonly used to study elusive or rare species where direct observation or capture is difficult. Little attention has been paid to the potential effects of observer bias while collecting noninvasive genetic samples in the field, however. Over a period of 7 years, we examined whether different observers (n = 58) and observer experience influenced detection, amplification rates, and correct species identification of 4,836 gray wolf (Canis lupus) fecal samples collected in Idaho and Yellowstone National Park, USA and southwestern Alberta, Canada (2008-2014). We compared new observers (n = 33) to experienced observers (n = 25) and hypothesized experience level would increase the overall success of using NGS techniques in the wild. In contrast to our hypothesis, we found that new individuals were better than experienced observers at detecting and collecting wolf scats and correctly identifying wolf scats from other sympatric carnivores present in the study areas. While adequate training of new observers is crucial for the successful use of NGS techniques, attention should also be directed to experienced observers. Observer experience could be a curse because of their potential effects on NGS data quality arising from fatigue, boredom or other factors. The ultimate benefit of an observer to a project is a combination of factors (i.e., field savvy, local knowledge), but project investigators should be aware of the potential negative effects of experience on NGS sampling.

Population-level inferences from environmental DNA-Current status and future perspectives

Environmental DNA (eDNA) extracted from water samples has recently shown potential as a valuable source of population genetic information for aquatic macroorganisms. This approach offers several potential advantages compared with conventional tissue-based methods, including the fact that eDNA sampling is noninvasive and generally more cost-efficient. Currently, eDNA approaches have been limited to single-marker studies of mitochondrial DNA (mtDNA), and the relationship between eDNA haplotype composition and true haplotype composition still needs to be thoroughly verified. This will require testing of bioinformatic and statistical software to correct for erroneous sequences, as well as biases and random variation in relative sequence abundances. However, eDNA-based population genetic methods have far-reaching potential for both basic and applied research. In this paper, we present a brief overview of the achievements of eDNA-based population genetics to date, and outline the prospects for future developments in the field, including the estimation of nuclear DNA (nuDNA) variation and epigenetic information. We discuss the challenges associated with eDNA samples as opposed to those of individual tissue samples and assess whether eDNA might offer additional types of information unobtainable with tissue samples. Lastly, we provide recommendations for determining whether an eDNA approach would be a useful and suitable choice in different research settings. We limit our discussion largely to contemporary aquatic systems, but the advantages, challenges, and perspectives can to a large degree be generalized to eDNA studies with a different spatial and temporal focus.

Standardization and validation of a panel of cross-species microsatellites to individually identify the Asiatic wild dog (Cuon alpinus)

Background

The Asiatic wild dog or dhole (Cuon alpinus) is a highly elusive, monophyletic, forest dwelling, social canid distributed across south and Southeast Asia. Severe pressures from habitat loss, prey depletion, disease, human persecution and interspecific competition resulted in global population decline in dholes. Despite a declining population trend, detailed information on population size, ecology, demography and genetics is lacking. Generating reliable information at landscape level for dholes is challenging due to their secretive behaviour and monomorphic physical features. Recent advances in non-invasive DNA-based tools can be used to monitor populations and individuals across large landscapes. In this paper, we describe standardization and validation of faecal DNA-based methods for individual identification of dholes. We tested this method on 249 field-collected dhole faeces from five protected areas of the central Indian landscape in the state of Maharashtra, India.

Results

We tested a total of 18 cross-species markers and developed a panel of 12 markers for unambiguous individual identification of dholes. This marker panel identified 101 unique individuals from faecal samples collected across our pilot field study area. These loci showed varied level of amplification success (57–88%), polymorphism (3–9 alleles), heterozygosity (0.23–0.63) and produced a cumulative misidentification rate or PID_(unbiased) and PID_(sibs) value of 4.7 × 10⁻¹⁰ and 1.5 × 10⁻⁴, respectively, indicating a high statistical power in individual discrimination from poor quality samples.

Conclusion

Our results demonstrated that the selected panel of 12 microsatellite loci can conclusively identify dholes from poor quality, non-invasive biological samples and help in exploring various population parameters. This genetic approach would be useful in dhole population estimation across its range and will help in assessing population trends and other genetic parameters for this elusive, social carnivore.

Beyond Biodiversity: Can Environmental DNA (eDNA) Cut It as a Population Genetics Tool?

Population genetic data underpin many studies of behavioral, ecological, and evolutionary processes in wild populations and contribute to effective conservation management. However, collecting genetic samples can be challenging when working with endangered, invasive, or cryptic species. Environmental DNA (eDNA) offers a way to sample genetic material non-invasively without requiring visual observation. While eDNA has been trialed extensively as a biodiversity and biosecurity monitoring tool with a strong taxonomic focus, it has yet to be fully explored as a means for obtaining population genetic information. Here, we review current research that employs eDNA approaches for the study of populations. We outline challenges facing eDNA-based population genetic methodologies, and suggest avenues of research for future developments. We advocate that with further optimizations, this emergent field holds great potential as part of the population genetics toolkit.

Complex problems need detailed solutions: Harnessing multiple data types to inform genetic management in the wild

For bottlenecked populations of threatened species, supplementation often leads to improved population metrics (genetic rescue), provided that guidelines can be followed to avoid negative outcomes. In cases where no "ideal" source populations exist, or there are other complicating factors such as prevailing disease, the benefit of supplementation becomes uncertain. Bringing multiple data and analysis types together to plan genetic management activities can help. Here, we consider three populations of Tasmanian devil, Sarcophilus harrisii, as candidates for genetic rescue. Since 1996, devil populations have been severely impacted by devil facial tumour disease (DFTD), causing significant population decline and fragmentation. Like many threatened species, the key threatening process for devils cannot currently be fully mitigated, so species management requires a multifaceted approach. We examined diversity of 31 putatively neutral and 11 MHC-linked microsatellite loci of three remnant wild devil populations (one sampled at two time-points), alongside computational diversity projections, parameterized by field data from DFTD-present and DFTD-absent sites. Results showed that populations had low diversity, connectivity was poor, and diversity has likely decreased over the last decade. Stochastic simulations projected further diversity losses. For a given population size, the effects of DFTD on population demography (including earlier age at death and increased female productivity) did not impact diversity retention, which was largely driven by final population size. Population sizes ≥500 (depending on the number of founders) were necessary for maintaining diversity in otherwise unmanaged populations, even if DFTD is present. Models indicated that smaller populations could maintain diversity with ongoing immigration. Taken together, our results illustrate how multiple analysis types can be combined to address complex population genetic challenges.

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.

Biases of STRUCTURE software when exploring introduction routes of invasive species

Population genetic methods are widely used to retrace the introduction routes of invasive species. The unsupervised Bayesian clustering algorithm implemented in STRUCTURE is amongst the most frequently used of these methods, but its ability to provide reliable information about introduction routes has never been assessed. We simulated microsatellite datasets to evaluate the extent to which the results provided by STRUCTURE were misleading for the inference of introduction routes. We focused on an invasion scenario involving one native and two independently introduced populations, because it is the sole scenario that can be rejected when obtaining a particular clustering with a STRUCTURE analysis at K = 2 (two clusters). Results were classified as "misleading" or "non-misleading". We investigated the influence of effective size, bottleneck severity and number of loci on the type and frequency of misleading results. We showed that misleading STRUCTURE results were obtained for 10% of all simulated datasets. Our results highlighted two categories of misleading output. The first occurs when the native population has a low level of diversity. In this case, the two introduced populations may be very similar, despite their independent introduction histories. The second category results from convergence issues in STRUCTURE for K = 2, with strong bottleneck severity and/or large numbers of loci resulting in high levels of differentiation between the three populations. Overall, the risk of being misled by STRUCTURE in the context of introduction routes inferences is moderate, but it is important to remain cautious when low genetic diversity or genuine multimodality between runs are involved.

Non-invasive genetic monitoring involving citizen science enables reconstruction of current pack dynamics in a re-establishing wolf population

Background

Carnivores are re-establishing in many human-populated areas, where their presence is often contentious. Reaching consensus on management decisions is often hampered by a dispute over the size of the local carnivore population. Understanding the reproductive dynamics and individual movements of the carnivores can provide support for management decisions, but individual-level information can be difficult to obtain from elusive, wide-ranging species. Non-invasive genetic sampling can yield such information, but makes subsequent reconstruction of population history challenging due to incomplete population coverage and error-prone data. Here, we combine a collaborative, volunteer-based sampling scheme with Bayesian pedigree reconstruction to describe the pack dynamics of an establishing grey wolf (Canis lupus) population in south-west Finland, where wolf breeding was recorded in 2006 for the first time in over a century.

Results

Using DNA extracted mainly from faeces collected since 2008, we identified 81 individual wolves and assigned credible full parentages to 70 of these and partial parentages to a further 9, revealing 7 breeding pairs. Individuals used a range of strategies to obtain breeding opportunities, including dispersal to established or new packs, long-distance migration and inheriting breeding roles. Gene flow occurred between all packs but inbreeding events were rare.

Conclusions

These findings demonstrate that characterizing ongoing pack dynamics can provide detailed, locally-relevant insight into the ecology of contentious species such as the wolf. Involving various stakeholders in data collection makes these results more likely to be accepted as unbiased and hence reliable grounds for management decisions.

Electronic supplementary material

The online version of this article (10.1186/s12898-017-0154-8) contains supplementary material, which is available to authorized users.

Assessing SNP genotyping of noninvasively collected wildlife samples using microfluidic arrays

Noninvasively collected samples are a common source of DNA in wildlife genetic studies. Currently, single nucleotide polymorphism (SNP) genotyping using microfluidic arrays is emerging as an easy-to-use and cost-effective methodology. Here we assessed the performance of microfluidic SNP arrays in genotyping noninvasive samples from grey wolves, European wildcats and brown bears, and we compared results with traditional microsatellite genotyping. We successfully SNP-genotyped 87%, 80% and 97% of the wolf, cat and bear samples, respectively. Genotype recovery was higher based on SNPs, while both marker types identified the same individuals and provided almost identical estimates of pairwise differentiation. We found that samples for which all SNP loci were scored had no disagreements across the three replicates (except one locus in a wolf sample). Thus, we argue that call rate (amplification success) can be used as a proxy for genotype quality, allowing the reduction of replication effort when call rate is high. Furthermore, we used cycle threshold values of real-time PCR to guide the choice of protocols for SNP amplification. Finally, we provide general guidelines for successful SNP genotyping of degraded DNA using microfluidic technology.

Genetic Divergence in Domestic Japanese Quail Inferred from Mitochondrial DNA D-Loop and Microsatellite Markers

To assess the genetic diversity of domestic Japanese quail (Coturnix japonica) populations, and their genetic relationships, we examined mitochondrial DNA (mtDNA) D-loop sequences and microsatellite markers for 19 Japanese quail populations. The populations included nine laboratory lines established in Japan (LWC, Quv, RWN, WE, AWE, AMRP, rb-TKP, NIES-L, and W), six meat-type quail lines reimported from Western countries (JD, JW, Estonia, NIES-Br, NIES-Fr, and NIES-Hn), one commercial population in Japan, and three wild quail populations collected from three Asian areas. The phylogenetic tree of mtDNA D-loop sequences revealed two distinct haplotype groups, Dloop-Group1 and Dloop-Group2. Dloop-Group1 included a dominant haplotype representing most of the quail populations, including wild quail. Dloop-Group2 was composed of minor haplotypes found in several laboratory lines, two meat-type lines, and a few individuals in commercial and wild quail populations. Taking the breeding histories of domestic populations into consideration, these results suggest that domestic quail populations may have derived from two sources, i.e., domestic populations established before and after World War II in Japan. A discriminant analysis of principal components and a Bayesian clustering analysis with microsatellite markers indicated that the domestic populations are clustered into four genetic groups. The two major groups were Microsat-Group1, which contained WE, and four WE-derived laboratory lines (LWC, Quv, RWN, and AWE), and Microsat-Group2 consisting of NIES-L, JD, JW, Estonia, NIES-Br, NIES-Fr, NIES-Hn, W, and commercial and wild populations. The remaining two lines (AMRP and rb-TKP) were each clustered into a separate clade. This hierarchical genetic difference between domestic quail populations is attributed to the genetic background derived from two different genetic sources-the pre-war and post-war populations-which is well supported by their breeding histories.

COMPARISON OF THREE SHORT-TERM IMMOBILIZATION REGIMES IN WILD VERREAUX'S SIFAKAS (PROPITHECUS VERREAUXI): KETAMINE-XYLAZINE, KETAMINE-XYLAZINE-ATROPINE, AND TILETAMINE-ZOLAZEPAM

Although research on lemurid primates in Madagascar has been ongoing for several decades, reports on different drug regimes to immobilize wild lemurs are limited. This study compares the efficacy, reliability, and side effects of ketamine-xylazine, ketamine-xylazine-atropine, and tiletamine-zolazepam immobilization in wild Verreaux's sifakas (Propithecus verreauxi). In the course of a long-term study in Kirindy Forest, western Madagascar, eight animals each received a mixture of ketamine (5.32±1.71 mg/kg) and xylazine (0.56±0.19 mg/kg) (KX; 7 males, 1 female) and ketamine (6.58±1.36 mg/kg), xylazine (1.28±0.28 mg/kg), and atropine (0.013±0.003 mg/kg) (KXA; 5 males, 3 females), respectively, and 14 individuals received tiletamine-zolazepam (7.73±1.37 mg/kg) (TZ; 9 males, 5 females). Induction was smooth in all protocols, but showed considerable variation in duration when animals had received KXA. Immobilization as well as recovery lasted significantly longer with TZ than with KX (P<0.05). Occurrence of side effects was not significantly different between the protocols; however, excessive salivation, involuntary muscular contractions, and vocalization only occurred in animals immobilized with TZ. Heart rate measurement at 10 min after onset of complete immobilization yielded significantly higher values if the animals had been immobilized with TZ compared to KX (P<0.05). Heart rate decreased from the first to the second measurement for the KX- and KXA-immobilized animals, whereas immobilization with TZ resulted in an increase in heart rate. The results suggest that KX produces good, but short, immobilization in Verreaux's sifakas at approximately 5 mg/kg ketamine and 0.5 mg/kg xylazine and a smoother and shorter recovery phase than 5 to 10 mg/kg TZ, whereas adding atropine to KX did not provide any benefits.

From promise to practice: pairing non-invasive sampling with genomics in conservation

Conservation genomics has become an increasingly popular term, yet it remains unclear whether the non-invasive sampling that is essential for many conservation-related studies is compatible with the minimum requirements for harnessing next-generation sequencing technologies. Here, we evaluated the feasibility of using genotyping-by-sequencing of non-invasively collected hair samples to simultaneously identify and genotype single nucleotide polymorphisms (SNPs) in a climate-sensitive mammal, the American pika (Ochotona princeps). We identified and genotyped 3,803 high-confidence SNPs across eight sites distributed along two elevational transects using starting DNA amounts as low as 1 ng. Fifty-five outlier loci were detected as candidate gene regions under divergent selection, constituting potential targets for future validation. Genome-wide estimates of gene diversity significantly and positively correlated with elevation across both transects, with all low elevation sites exhibiting significant heterozygote deficit likely due to inbreeding. More broadly, our results highlight a range of issues that must be considered when pairing genomic data collection with non-invasive sampling, particularly related to field sampling protocols for minimizing exogenous DNA, data collection strategies and quality control steps for enhancing target organism yield, and analytical approaches for maximizing cost-effectiveness and information content of recovered genomic data.