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

Coevolution promotes the coexistence of Tasmanian devils and a fatal, transmissible cancer

Emerging infectious diseases threaten natural populations, and data-driven modeling is critical for predicting population dynamics. Despite the importance of integrating ecology and evolution in models of host-pathogen dynamics, there are few wild populations for which long-term ecological datasets have been coupled with genome-scale data. Tasmanian devil (Sarcophilus harrisii) populations have declined range wide due to devil facial tumor disease (DFTD), a fatal transmissible cancer. Although early ecological models predicted imminent devil extinction, diseased devil populations persist at low densities, and recent ecological models predict long-term devil persistence. Substantial evidence supports the evolution of both devils and DFTD, suggesting coevolution may also influence continued devil persistence. Thus, we developed an individual-based, eco-evolutionary model of devil-DFTD coevolution parameterized with nearly 2 decades of devil demography, DFTD epidemiology, and genome-wide association studies. We characterized potential devil-DFTD coevolutionary outcomes and predicted the effects of coevolution on devil persistence and devil-DFTD coexistence. We found a high probability of devil persistence over 50 devil generations (100 years) and a higher likelihood of devil-DFTD coexistence, with greater devil recovery than predicted by previous ecological models. These novel results add to growing evidence for long-term devil persistence and highlight the importance of eco-evolutionary modeling for emerging infectious diseases.

Restoring faith in conservation action: Maintaining wild genetic diversity through the Tasmanian devil insurance program

Conservation breeding programs aim to maintain 90% wild genetic diversity, but rarely assess functional diversity. Here, we compare both genome-wide and functional diversity (in over 500 genes) of Tasmanian devils (Sarcophilus harrisii) within the insurance metapopulation and across the species’ range (64,519 km²). Populations have declined by 80% since 1996 due to a contagious cancer, devil facial tumor disease (DFTD). However, predicted local extinctions have not occurred. Recent suggestions of selection for “resistance” alleles in the wild precipitated concerns that insurance population devils may be unsuitable for translocations. Using 830 wild samples collected at 31 locations between 2012 and 2021, and 553 insurance metapopulation devils, we show that the insurance metapopulation is representative of current wild genetic diversity. Allele frequencies at DFTD-associated loci were not substantially different between captive and wild devils. Methods presented here are valuable for others investigating evolutionary potential in threatened species, particularly ones under significant selective pressures.

•

Developed target capture to assess functional diversity at over 500 genes

•

Fine-scale structure exists in the genetically depauperate Tasmanian devil

•

Insurance metapopulation is representative of wild genetic diversity

•

Allele frequencies at disease-associated loci were similar in captivity to the wild

How much is enough? Sampling intensity influences estimates of reproductive variance in an introduced population

Conservation introductions to islands and fenced enclosures are increasing as in situ mitigations fail to keep pace with population declines. Few studies consider the potential loss of genetic diversity and increased inbreeding if released individuals breed disproportionately. As funding is limited and post-release monitoring expensive for conservation programs, understanding how sampling effort influences estimates of reproductive variance is useful. To investigate this relationship, we used a well-studied population of Tasmanian devils (Sarcophilus harrisii) introduced to Maria Island, Tasmania, Australia. Pedigree reconstruction based on molecular data revealed high variance in number of offspring per breeder and high proportions of unsuccessful individuals. Computational subsampling of 20%, 40%, 60%, and 80% of observed offspring resulted in inaccurate estimates of reproductive variance compared to the pedigree reconstructed with all sampled individuals. With decreased sampling effort, the proportion of inferred unsuccessful individuals was overestimated and the variance in number of offspring per breeder was underestimated. To accurately estimate reproductive variance, we recommend sampling as many individuals as logistically possible during the early stages of population establishment. Further, we recommend careful selection of colonizing individuals as they may be disproportionately represented in subsequent generations. Within the conservation management context, our results highlight important considerations for sample collection and post-release monitoring during population establishment.

Global Commitments to Conserving and Monitoring Genetic Diversity Are Now Necessary and Feasible

Global conservation policy and action have largely neglected protecting and monitoring genetic diversity—one of the three main pillars of biodiversity. Genetic diversity (diversity within species) underlies species’ adaptation and survival, ecosystem resilience, and societal innovation. The low priority given to genetic diversity has largely been due to knowledge gaps in key areas, including the importance of genetic diversity and the trends in genetic diversity change; the perceived high expense and low availability and the scattered nature of genetic data; and complicated concepts and information that are inaccessible to policymakers. However, numerous recent advances in knowledge, technology, databases, practice, and capacity have now set the stage for better integration of genetic diversity in policy instruments and conservation efforts. We review these developments and explore how they can support improved consideration of genetic diversity in global conservation policy commitments and enable countries to monitor, report on, and take action to maintain or restore genetic diversity.

First evidence of deviation from Mendelian proportions in a conservation programme

Classic Mendelian inheritance is the bedrock of population genetics and underpins pedigree-based management of animal populations. However, assumptions of Mendelian inheritance might not be upheld in conservation breeding programmes if early viability selection occurs, even when efforts are made to equalise genetic contributions of breeders. To test this possibility, we investigated deviations from Mendelian proportions in a captive metapopulation of the endangered Tasmanian devil. This marsupial population is ideal for addressing evolutionary questions in conservation due to its large size, range of enclosure types (varying in environmental conditions), good genomic resources (which aid interpretation), and the species' biology. Devil mothers give birth to more offspring than they can nurse in the pouch, providing the potential for intense viability selection amongst embryos. We used data from 140 known sire-dam-offspring triads to isolate within-family selection from population-level mechanisms (such as mate choice or inbreeding), and compared observed offspring genotypes at 123 targeted SNPs to neutral (i.e., Mendelian) expectations. We found lower offspring heterozygosity than expected, and subtle patterns that varied across a gradient of management intensity from zoo-like enclosures to semi-wild environments for some loci. Meiotic drive or maternal-foetal incompatibilities are consistent with our results, although we cannot statistically confirm these mechanisms. We found some evidence that maternal genotype affects annual litter size, suggesting that family-level patterns are driven by differential offspring mortality before birth or during early development. Our results show that deviations from Mendelian inheritance can occur in conservation programmes, despite best-practice management to prevent selection.

Big Data in Conservation Genomics: Boosting Skills, Hedging Bets, and Staying Current in the Field

A current challenge in the fields of evolutionary, ecological, and conservation genomics is balancing production of large-scale datasets with additional training often required to handle such datasets. Thus, there is an increasing need for conservation geneticists to continually learn and train to stay up-to-date through avenues such as symposia, meetings, and workshops. The ConGen meeting is a near-annual workshop that strives to guide participants in understanding population genetics principles, study design, data processing, analysis, interpretation, and applications to real-world conservation issues. Each year of ConGen gathers a diverse set of instructors, students, and resulting lectures, hands-on sessions, and discussions. Here, we summarize key lessons learned from the 2019 meeting and more recent updates to the field with a focus on big data in conservation genomics. First, we highlight classical and contemporary issues in study design that are especially relevant to working with big datasets, including the intricacies of data filtering. We next emphasize the importance of building analytical skills and simulating data, and how these skills have applications within and outside of conservation genetics careers. We also highlight recent technological advances and novel applications to conservation of wild populations. Finally, we provide data and recommendations to support ongoing efforts by ConGen organizers and instructors-and beyond-to increase participation of underrepresented minorities in conservation and eco-evolutionary sciences. The future success of conservation genetics requires both continual training in handling big data and a diverse group of people and approaches to tackle key issues, including the global biodiversity-loss crisis.

Reducing the Extinction Risk of Populations Threatened by Infectious Diseases

Extinction risk is increasing for a range of species due to a variety of threats, including disease. Emerging infectious diseases can cause severe declines in wild animal populations, increasing population fragmentation and reducing gene flow. Small, isolated, host populations may lose adaptive potential and become more susceptible to extinction due to other threats. Management of the genetic consequences of disease-induced population decline is often necessary. Whilst disease threats need to be addressed, they can be difficult to mitigate. Actions implemented to conserve the Tasmanian devil (Sarcophilus harrisii), which has suffered decline to the deadly devil facial tumour disease (DFTD), exemplify how genetic management can be used to reduce extinction risk in populations threatened by disease. Supplementation is an emerging conservation technique that may benefit populations threatened by disease by enabling gene flow and conserving their adaptive potential through genetic restoration. Other candidate species may benefit from genetic management via supplementation but concerns regarding outbreeding depression may prevent widespread incorporation of this technique into wildlife disease management. However, existing knowledge can be used to identify populations that would benefit from supplementation where risk of outbreeding depression is low. For populations threatened by disease and, in situations where disease eradication is not an option, wildlife managers should consider genetic management to buffer the host species against inbreeding and loss of genetic diversity.

Genetic Consequences of Multiple Translocations of the Banded Hare-Wallaby in Western Australia

Many Australian mammal species now only occur on islands and fenced mainland havens free from invasive predators. The range of one species, the banded hare-wallaby (Lagostrophus fasciatus), had contracted to two offshore islands in Western Australia. To improve survival, four conservation translocations have been attempted with mixed success, and all occurred in the absence of genetic information. Here, we genotyped seven polymorphic microsatellite markers in two source (Bernier Island and Dorre Island), two historic captive, and two translocated L. fasciatus populations to determine the impact of multiple translocations on genetic diversity. Subsequently, we used population viability analysis (PVA) and gene retention modelling to determine scenarios that will maximise demographic resilience and genetic richness of two new populations that are currently being established. One translocated population (Wadderin) has undergone a genetic bottleneck and lost 8.1% of its source population’s allelic diversity, while the other (Faure Island) may be inbred. We show that founder number is a key parameter when establishing new L. fasciatus populations and 100 founders should lead to high survival probabilities. Our modelling predicts that during periodic droughts, the recovery of source populations will be slower post-harvest, while 75% more animals—about 60 individuals—are required to retain adequate allelic diversity in the translocated population. Our approach demonstrates how genetic data coupled with simulations of stochastic environmental events can address central questions in translocation programmes.

Curse of the devil: molecular insights into the emergence of transmissible cancers in the Tasmanian devil (Sarcophilus harrisii)

The Tasmanian devil (Sarcophilus harrisii) is the only mammalian species known to be affected by multiple transmissible cancers. Devil facial tumours 1 and 2 (DFT1 and DFT2) are independent neoplastic cell lineages that produce large, disfiguring cancers known as devil facial tumour disease (DFTD). The long-term persistence of wild Tasmanian devils is threatened due to the ability of DFTD cells to propagate as contagious allografts and the high mortality rate of DFTD. Recent studies have demonstrated that both DFT1 and DFT2 cancers originated from founder cells of the Schwann cell lineage, an uncommon origin of malignant cancer in humans. This unprecedented finding has revealed a potential predisposition of Tasmanian devils to transmissible cancers of the Schwann cell lineage. In this review, we compare the molecular nature of human Schwann cells and nerve sheath tumours with DFT1 and DFT2 to gain insights into the emergence of transmissible cancers in the Tasmanian devil. We discuss a potential mechanism, whereby Schwann cell plasticity and frequent wounding in Tasmanian devils combine with an inherent cancer predisposition and low genetic diversity to give rise to transmissible Schwann cell cancers in devils on rare occasions.

Inbreeding depression in one of the last DFTD-free wild populations of Tasmanian devils

BACKGROUND

Vulnerable species experiencing inbreeding depression are prone to localised extinctions because of their reduced fitness. For Tasmanian devils, the rapid spread of devil facial tumour disease (DFTD) has led to population declines and fragmentation across the species' range. Here we show that one of the few remaining DFTD-free populations of Tasmanian devils is experiencing inbreeding depression. Moreover, this population has experienced a significant reduction in reproductive success over recent years.

METHODS

We used 32 microsatellite loci to examine changes in genetic diversity and inbreeding in the wild population at Woolnorth, alongside field data on breeding success from females to test for inbreeding depression.

RESULTS

Wefound that maternal internal relatedness has a negative impact on litter sizes. The results of this study imply that this population may be entering an extinction vortex and that to protect the population genetic rescue should be considered. This study provides conservation managers with useful information for managing wild devils and provides support for the "Wild Devil Recovery Program", which is currently augmenting small, isolated populations.

A comparison of pedigree, genetic and genomic estimates of relatedness for informing pairing decisions in two critically endangered birds: Implications for conservation breeding programmes worldwide

Conservation management strategies for many highly threatened species include conservation breeding to prevent extinction and enhance recovery. Pairing decisions for these conservation breeding programmes can be informed by pedigree data to minimize relatedness between individuals in an effort to avoid inbreeding, maximize diversity and maintain evolutionary potential. However, conservation breeding programmes struggle to use this approach when pedigrees are shallow or incomplete. While genetic data (i.e., microsatellites) can be used to estimate relatedness to inform pairing decisions, emerging evidence indicates this approach may lack precision in genetically depauperate species, and more effective estimates will likely be obtained from genomic data (i.e., thousands of genome-wide single nucleotide polymorphisms, or SNPs). Here, we compare relatedness estimates and subsequent pairing decisions using pedigrees, microsatellites and SNPs from whole-genome resequencing approaches in two critically endangered birds endemic to New Zealand: kakī/black stilt (Himantopus novaezelandiae) and kākāriki karaka/orange-fronted parakeet (Cyanoramphus malherbi). Our findings indicate that SNPs provide more precise estimates of relatedness than microsatellites when assessing empirical parent-offspring and full sibling relationships. Further, our results show that relatedness estimates and subsequent pairing recommendations using PMx are most similar between pedigree- and SNP-based approaches. These combined results indicate that in lieu of robust pedigrees, SNPs are an effective tool for informing pairing decisions, which has important implications for many poorly pedigreed conservation breeding programmes worldwide.

Targeted conservation genetics of the endangered chimpanzee

Populations of the common chimpanzee (Pan troglodytes) are in an impending risk of going extinct in the wild as a consequence of damaging anthropogenic impact on their natural habitat and illegal pet and bushmeat trade. Conservation management programmes for the chimpanzee have been established outside their natural range (ex situ), and chimpanzees from these programmes could potentially be used to supplement future conservation initiatives in the wild (in situ). However, these programmes have often suffered from inadequate information about the geographical origin and subspecies ancestry of the founders. Here, we present a newly designed capture array with ~60,000 ancestry informative markers used to infer ancestry of individual chimpanzees in ex situ populations and determine geographical origin of confiscated sanctuary individuals. From a test panel of 167 chimpanzees with unknown origins or subspecies labels, we identify 90 suitable non-admixed individuals in the European Association of Zoos and Aquaria (EAZA) Ex situ Programme (EEP). Equally important, another 46 individuals have been identified with admixed subspecies ancestries, which therefore over time, should be naturally phased out of the breeding populations. With potential for future re-introduction to the wild, we determine the geographical origin of 31 individuals that were confiscated from the illegal trade and demonstrate the promises of using non-invasive sampling in future conservation action plans. Collectively, our genomic approach provides an exemplar for ex situ management of endangered species and offers an efficient tool in future in situ efforts to combat the illegal wildlife trade.

An oral bait vaccination approach for the Tasmanian devil facial tumor diseases

Introduction: The Tasmanian devil (Sarcophilus harrisii) is the largest extant carnivorous marsupial. Since 1996, its population has declined by 77% primarily due to a clonal transmissible tumor, known as devil facial tumor (DFT1) disease. In 2014, a second transmissible devil facial tumor (DFT2) was discovered. DFT1 and DFT2 are nearly 100% fatal.Areas covered: We review DFT control approaches and propose a rabies-style oral bait vaccine (OBV) platform for DFTs. This approach has an extensive safety record and was a primary tool in large-scale rabies virus elimination from wild carnivores across diverse landscapes. Like rabies virus, DFTs are transmitted by oral contact, so immunizing the oral cavity and stimulating resident memory cells could be advantageous. Additionally, exposing infected devils that already have tumors to OBVs could serve as an oncolytic virus immunotherapy. The primary challenges may be identifying appropriate DFT-specific antigens and optimization of field delivery methods.Expert opinion: DFT2 is currently found on a peninsula in southern Tasmania, so an OBV that could eliminate DFT2 should be the priority for this vaccine approach. Translation of an OBV approach to control DFTs will be challenging, but the approach is feasible for combatting ongoing and future disease threats.

The Value of Reference Genomes in the Conservation of Threatened Species

Conservation initiatives are now more crucial than ever-over a million plant and animal species are at risk of extinction over the coming decades. The genetic management of threatened species held in insurance programs is recommended; however, few are taking advantage of the full range of genomic technologies available today. Less than 1% of the 13505 species currently listed as threated by the International Union for Conservation of Nature (IUCN) have a published genome. While there has been much discussion in the literature about the importance of genomics for conservation, there are limited examples of how having a reference genome has changed conservation management practice. The Tasmanian devil (Sarcophilus harrisii), is an endangered Australian marsupial, threatened by an infectious clonal cancer devil facial tumor disease (DFTD). Populations have declined by 80% since the disease was first recorded in 1996. A reference genome for this species was published in 2012 and has been crucial for understanding DFTD and the management of the species in the wild. Here we use the Tasmanian devil as an example of how a reference genome has influenced management actions in the conservation of a species.

From reference genomes to population genomics: comparing three reference-aligned reduced-representation sequencing pipelines in two wildlife species

Background

Recent advances in genomics have greatly increased research opportunities for non-model species. For wildlife, a growing availability of reference genomes means that population genetics is no longer restricted to a small set of anonymous loci. When used in conjunction with a reference genome, reduced-representation sequencing (RRS) provides a cost-effective method for obtaining reliable diversity information for population genetics. Many software tools have been developed to process RRS data, though few studies of non-model species incorporate genome alignment in calling loci. A commonly-used RRS analysis pipeline, Stacks, has this capacity and so it is timely to compare its utility with existing software originally designed for alignment and analysis of whole genome sequencing data. Here we examine population genetic inferences from two species for which reference-aligned reduced-representation data have been collected. Our two study species are a threatened Australian marsupial (Tasmanian devil Sarcophilus harrisii; declining population) and an Arctic-circle migrant bird (pink-footed goose Anser brachyrhynchus; expanding population). Analyses of these data are compared using Stacks versus two widely-used genomics packages, SAMtools and GATK. We also introduce a custom R script to improve the reliability of single nucleotide polymorphism (SNP) calls in all pipelines and conduct population genetic inferences for non-model species with reference genomes.

Results

Although we identified orders of magnitude fewer SNPs in our devil dataset than for goose, we found remarkable symmetry between the two species in our assessment of software performance. For both datasets, all three methods were able to delineate population structure, even with varying numbers of loci. For both species, population structure inferences were influenced by the percent of missing data.

Conclusions

For studies of non-model species with a reference genome, we recommend combining Stacks output with further filtering (as included in our R pipeline) for population genetic studies, paying particular attention to potential impact of missing data thresholds. We recognise SAMtools as a viable alternative for researchers more familiar with this software. We caution against the use of GATK in studies with limited computational resources or time.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5806-y) contains supplementary material, which is available to authorized users.