Extensive population decline in the Tasmanian devil predates European settlement and devil facial tumour disease

Population genomics of a predatory mammal reveals patterns of decline and impacts of exposure to toxic toads

Mammal declines across northern Australia are one of the major biodiversity loss events occurring globally. There has been no regional assessment of the implications of these species declines for genomic diversity. To address this, we conducted a species-wide assessment of genomic diversity in the northern quoll (Dasyurus hallucatus), an Endangered marsupial carnivore. We used next generation sequencing methods to genotype 10,191 single nucleotide polymorphisms (SNPs) in 352 individuals from across a 3220-km length of the continent, investigating patterns of population genomic structure and diversity, and identifying loci showing signals of putative selection. We found strong heterogeneity in the distribution of genomic diversity across the continent, characterized by (i) biogeographical barriers driving hierarchical population structure through long-term isolation, and (ii) severe reductions in diversity resulting from population declines, exacerbated by the spread of introduced toxic cane toads (Rhinella marina). These results warn of a large ongoing loss of genomic diversity and associated adaptive capacity as mammals decline across northern Australia. Encouragingly, populations of the northern quoll established on toad-free islands by translocations appear to have maintained most of the initial genomic diversity after 16 years. By mapping patterns of genomic diversity within and among populations, and investigating these patterns in the context of population declines, we can provide conservation managers with data critical to informed decision-making. This includes the identification of populations that are candidates for genetic management, the importance of remnant island and insurance/translocated populations for the conservation of genetic diversity, and the characterization of putative evolutionarily significant units.

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

Museum genomics reveals the rapid decline and extinction of Australian rodents since European settlement

Australia has the highest historically recorded rate of mammalian extinction in the world, with 34 terrestrial species declared extinct since European colonization in 1788. Among Australian mammals, rodents have been the most severely affected by these recent extinctions; however, given a sparse historical record, the scale and timing of their decline remain unresolved. Using museum specimens up to 184 y old, we generate genomic-scale data from across the entire assemblage of Australian hydromyine rodents (i.e., eight extinct species and their 42 living relatives). We reconstruct a phylogenomic tree for these species spanning ∼5.2 million years, revealing a cumulative total of 10 million years (>10%) of unique evolutionary history lost to extinction within the past ∼150 y. We find no evidence for reduced genetic diversity in extinct species just prior to or during decline, indicating that their extinction was extremely rapid. This suggests that populations of extinct Australian rodents were large prior to European colonization, and that genetic diversity does not necessarily protect species from catastrophic extinction. In addition, comparative analyses suggest that body size and biome interact to predict extinction and decline, with larger species more likely to go extinct. Finally, we taxonomically resurrect a species from extinction, Gould's mouse (Pseudomys gouldii Waterhouse, 1839), which survives as an island population in Shark Bay, Western Australia (currently classified as Pseudomys fieldi Waite, 1896). With unprecedented sampling across a radiation of extinct and living species, we unlock a previously inaccessible historical perspective on extinction in Australia. Our results highlight the capacity of collections-based research to inform conservation and management of persisting species.

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.

Contemporary and historical selection in Tasmanian devils (Sarcophilus harrisii) support novel, polygenic response to transmissible cancer

Tasmanian devils (Sarcophilus harrisii) are evolving in response to a unique transmissible cancer, devil facial tumour disease (DFTD), first described in 1996. Persistence of wild populations and the recent emergence of a second independently evolved transmissible cancer suggest that transmissible cancers may be a recurrent feature in devils. Here, we compared signatures of selection across temporal scales to determine whether genes or gene pathways under contemporary selection (six to eight generations) have also been subject to historical selection (65-85 Myr). First, we used targeted sequencing, RAD-capture, in approximately 2500 devils in six populations to identify genomic regions subject to rapid evolution. We documented genome-wide contemporary evolution, including 186 candidate genes related to cell cycling and immune response. Then we used a molecular evolution approach to identify historical positive selection in devils compared to other marsupials and found evidence of selection in 1773 genes. However, we found limited overlap across time scales, with only 16 shared candidate genes, and no overlap in enriched functional gene sets. Our results are consistent with a novel, multi-locus evolutionary response of devils to DFTD. Our results can inform conservation by identifying high priority targets for genetic monitoring and guiding maintenance of adaptive potential in managed populations.

Offspring survival changes over generations of captive breeding

Conservation breeding programs such as zoos play a major role in preventing extinction, but their sustainability may be impeded by neutral and adaptive population genetic change. These changes are difficult to detect for a single species or context, and impact global conservation efforts. We analyse pedigree data from 15 vertebrate species – over 30,000 individuals – to examine offspring survival over generations of captive breeding. Even accounting for inbreeding, we find that the impacts of increasing generations in captivity are highly variable across species, with some showing substantial increases or decreases in offspring survival over generations. We find further differences between dam and sire effects in first- versus multi-generational analysis. Crucially, our multispecies analysis reveals that responses to captivity could not be predicted from species’ evolutionary (phylogenetic) relationships. Even under best-practice captive management, generational fitness changes that cannot be explained by known processes (such as inbreeding depression), are occurring.

Captive breeding could prevent species extinctions, but selection for captivity may decrease fitness. Here the authors analyse pedigree data on 15 long-running vertebrate breeding programs and find generational fitness changes that processes such as inbreeding depression cannot explain.

Evolution and lineage dynamics of a transmissible cancer in Tasmanian devils

Devil facial tumour 1 (DFT1) is a transmissible cancer clone endangering the Tasmanian devil. The expansion of DFT1 across Tasmania has been documented, but little is known of its evolutionary history. We analysed genomes of 648 DFT1 tumours collected throughout the disease range between 2003 and 2018. DFT1 diverged early into five clades, three spreading widely and two failing to persist. One clade has replaced others at several sites, and rates of DFT1 coinfection are high. DFT1 gradually accumulates copy number variants (CNVs), and its telomere lengths are short but constant. Recurrent CNVs reveal genes under positive selection, sites of genome instability, and repeated loss of a small derived chromosome. Cultured DFT1 cell lines have increased CNV frequency and undergo highly reproducible convergent evolution. Overall, DFT1 is a remarkably stable lineage whose genome illustrates how cancer cells adapt to diverse environments and persist in a parasitic niche.

Comparative landscape genetics reveals differential effects of environment on host and pathogen genetic structure in Tasmanian devils (Sarcophilus harrisii) and their transmissible tumour

Genetic structure in host species is often used to predict disease spread. However, host and pathogen genetic variation may be incongruent. Understanding landscape factors that have either concordant or divergent influence on host and pathogen genetic structure is crucial for wildlife disease management. Devil facial tumour disease (DFTD) was first observed in 1996 and has spread throughout almost the entire Tasmanian devil geographic range, causing dramatic population declines. Whereas DFTD is predominantly spread via biting among adults, devils typically disperse as juveniles, which experience low DFTD prevalence. Thus, we predicted little association between devil and tumour population structure and that environmental factors influencing gene flow differ between devils and tumours. We employed a comparative landscape genetics framework to test the influence of environmental factors on patterns of isolation by resistance (IBR) and isolation by environment (IBE) in devils and DFTD. Although we found evidence for broad-scale costructuring between devils and tumours, we found no relationship between host and tumour individual genetic distances. Further, the factors driving the spatial distribution of genetic variation differed for each. Devils exhibited a strong IBR pattern driven by major roads, with no evidence of IBE. By contrast, tumours showed little evidence for IBR and a weak IBE pattern with respect to elevation in one of two tumour clusters we identify herein. Our results warrant caution when inferring pathogen spread using host population genetic structure and suggest that reliance on environmental barriers to host connectivity may be ineffective for managing the spread of wildlife diseases. Our findings demonstrate the utility of comparative landscape genetics for identifying differential factors driving host dispersal and pathogen transmission.

The ecology and evolution of wildlife cancers: Applications for management and conservation

Ecological and evolutionary concepts have been widely adopted to understand host-pathogen dynamics, and more recently, integrated into wildlife disease management. Cancer is a ubiquitous disease that affects most metazoan species; however, the role of oncogenic phenomena in eco-evolutionary processes and its implications for wildlife management and conservation remains undeveloped. Despite the pervasive nature of cancer across taxa, our ability to detect its occurrence, progression and prevalence in wildlife populations is constrained due to logistic and diagnostic limitations, which suggests that most cancers in the wild are unreported and understudied. Nevertheless, an increasing number of virus-associated and directly transmissible cancers in terrestrial and aquatic environments have been detected. Furthermore, anthropogenic activities and sudden environmental changes are increasingly associated with cancer incidence in wildlife. This highlights the need to upscale surveillance efforts, collection of critical data and developing novel approaches for studying the emergence and evolution of cancers in the wild. Here, we discuss the relevance of malignant cells as important agents of selection and offer a holistic framework to understand the interplay of ecological, epidemiological and evolutionary dynamics of cancer in wildlife. We use a directly transmissible cancer (devil facial tumour disease) as a model system to reveal the potential evolutionary dynamics and broader ecological effects of cancer epidemics in wildlife. We provide further examples of tumour-host interactions and trade-offs that may lead to changes in life histories, and epidemiological and population dynamics. Within this framework, we explore immunological strategies at the individual level as well as transgenerational adaptations at the population level. Then, we highlight the need to integrate multiple disciplines to undertake comparative cancer research at the human-domestic-wildlife interface and their environments. Finally, we suggest strategies for screening cancer incidence in wildlife and discuss how to integrate ecological and evolutionary concepts in the management of current and future cancer epizootics.

Disease swamps molecular signatures of genetic-environmental associations to abiotic factors in Tasmanian devil (Sarcophilus harrisii) populations

Landscape genomics studies focus on identifying candidate genes under selection via spatial variation in abiotic environmental variables, but rarely by biotic factors (i.e., disease). The Tasmanian devil (Sarcophilus harrisii) is found only on the environmentally heterogeneous island of Tasmania and is threatened with extinction by a transmissible cancer, devil facial tumor disease (DFTD). Devils persist in regions of long-term infection despite epidemiological model predictions of species' extinction, suggesting possible adaptation to DFTD. Here, we test the extent to which spatial variation and genetic diversity are associated with the abiotic environment (i.e., climatic variables, elevation, vegetation cover) and/or DFTD. We employ genetic-environment association analyses using 6886 SNPs from 3287 individuals sampled pre- and post-disease arrival across the devil's geographic range. Pre-disease, we find significant correlations of allele frequencies with environmental variables, including 365 unique loci linked to 71 genes, suggesting local adaptation to abiotic environment. The majority of candidate loci detected pre-DFTD are not detected post-DFTD arrival. Several post-DFTD candidate loci are associated with disease prevalence and were in linkage disequilibrium with genes involved in tumor suppression and immune response. Loss of apparent signal of abiotic local adaptation post-disease suggests swamping by strong selection resulting from the rapid onset of DFTD.

Two of a kind: transmissible Schwann cell cancers in the endangered Tasmanian devil (Sarcophilus harrisii)

Devil facial tumour disease (DFTD) comprises two genetically distinct transmissible cancers (DFT1 and DFT2) endangering the survival of the Tasmanian devil (Sarcophilus harrisii) in the wild. DFT1 first arose from a cell of the Schwann cell lineage; however, the tissue-of-origin of the recently discovered DFT2 cancer is unknown. In this study, we compared the transcriptome and proteome of DFT2 tumours to DFT1 and normal Tasmanian devil tissues to determine the tissue-of-origin of the DFT2 cancer. Our findings demonstrate that DFT2 expresses a range of Schwann cell markers and exhibits expression patterns consistent with a similar origin to the DFT1 cancer. Furthermore, DFT2 cells express genes associated with the repair response to peripheral nerve damage. These findings suggest that devils may be predisposed to transmissible cancers of Schwann cell origin. The combined effect of factors such as frequent nerve damage from biting, Schwann cell plasticity and low genetic diversity may allow these cancers to develop on rare occasions. The emergence of two independent transmissible cancers from the same tissue in the Tasmanian devil presents an unprecedented opportunity to gain insight into cancer development, evolution and immune evasion in mammalian species.

Population differentiation and historical demography of the threatened snowy plover Charadrius nivosus (Cassin, 1858)

Delineating conservation units is a complex and often controversial process that is particularly challenging for highly vagile species. Here, we reassess population genetic structure and identify those populations of highest conservation value in the threatened snowy plover (Charadrius nivosus, Cassin, 1858), a partial migrant shorebird endemic to the Americas. We use four categories of genetic data—mitochondrial DNA (mtDNA), microsatellites, Z-linked and autosomal single nucleotide polymorphisms (SNPs)—to: (1) assess subspecies delineation and examine population structure (2) compare the sensitivity of the different types of genetic data to detect spatial genetic patterns, and (3) reconstruct demographic history of the populations analysed. Delineation of two traditionally recognised subspecies was broadly supported by all data. In addition, microsatellite and SNPs but not mtDNA supported the recognition of Caribbean snowy plovers (C. n. tenuirostris) and Floridian populations (eastern C. n. nivosus) as distinct genetic lineage and deme, respectively. Low migration rates estimated from autosomal SNPs (m < 0.03) reflect a general paucity of exchange between genetic lineages. In contrast, we detected strong unidirectional migration (m = 0.26) from the western into the eastern nivosus deme. Within western nivosus, we found no genetic differentiation between coastal Pacific and inland populations. The correlation between geographic and genetic distances was weak but significant for all genetic data sets. All demes showed signatures of bottlenecks occurring during the past 1000 years. We conclude that at least four snowy plover conservation units are warranted: in addition to subspecies nivosus and occidentalis, a third unit comprises the Caribbean tenuirostris lineage and a fourth unit the distinct eastern nivosus deme.

Contemporary Demographic Reconstruction Methods Are Robust to Genome Assembly Quality: A Case Study in Tasmanian Devils

Reconstructing species’ demographic histories is a central focus of molecular ecology and evolution. Recently, an expanding suite of methods leveraging either the sequentially Markovian coalescent (SMC) or the site-frequency spectrum has been developed to reconstruct population size histories from genomic sequence data. However, few studies have investigated the robustness of these methods to genome assemblies of varying quality. In this study, we first present an improved genome assembly for the Tasmanian devil using the Chicago library method. Compared with the original reference genome, our new assembly reduces the number of scaffolds (from 35,975 to 10,010) and increases the scaffold N90 (from 0.101 to 2.164 Mb). Second, we assess the performance of four contemporary genomic methods for inferring population size history (PSMC, MSMC, SMC++, Stairway Plot), using the two devil genome assemblies as well as simulated, artificially fragmented genomes that approximate the hypothesized demographic history of Tasmanian devils. We demonstrate that each method is robust to assembly quality, producing similar estimates of Ne when simulated genomes were fragmented into up to 5,000 scaffolds. Overall, methods reliant on the SMC are most reliable between ∼300 generations before present (gbp) and 100 kgbp, whereas methods exclusively reliant on the site-frequency spectrum are most reliable between the present and 30 gbp. Our results suggest that when used in concert, genomic methods for reconstructing species’ effective population size histories 1) can be applied to nonmodel organisms without highly contiguous reference genomes, and 2) are capable of detecting independently documented effects of historical geological events.

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.

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.

Evolution of Marsupial Genomes

Marsupial genomes, which are packaged into large chromosomes, provide a powerful resource for studying the mechanisms of genome evolution. The extensive and valuable body of work on marsupial cytogenetics, combined more recently with genome sequence data, has enabled prediction of the 2n = 14 karyotype ancestral to all marsupial families. The application of both chromosome biology and genome sequencing, or chromosomics, has been a necessary approach for various aspects of mammalian genome evolution, such as understanding sex chromosome evolution and the origin and evolution of transmissible tumors in Tasmanian devils. The next phase of marsupial genome evolution research will employ chromosomics approaches to begin addressing fundamental questions in marsupial genome evolution and chromosome evolution more generally. The answers to these complex questions will impact our understanding across a broad range of fields, including the genetics of speciation, genome adaptation to environmental stressors, and species management.

Transcriptomics of Tasmanian Devil (Sarcophilus Harrisii) Ear Tissue Reveals Homogeneous Gene Expression Patterns across a Heterogeneous Landscape

In an era of unprecedented global change, exploring patterns of gene expression among wild populations across their geographic range is crucial for characterizing adaptive potential. RNA-sequencing studies have successfully characterized gene expression differences among populations experiencing divergent environmental conditions in a wide variety of taxa. However, few of these studies have identified transcriptomic signatures to multivariate, environmental stimuli among populations in their natural environments. Herein, we aim to identify environmental and sex-driven patterns of gene expression in the Tasmanian devil (Sarcophilus harrisii), a critically endangered species that occupies a heterogeneous environment. We performed RNA-sequencing on ear tissue biopsies from adult male and female devils from three populations at the extremes of their geographic range. There were no transcriptome-wide patterns of differential gene expression that would be suggestive of significant, environmentally-driven transcriptomic responses. The general lack of transcriptome-wide variation in gene expression levels across the devil's geographic range is consistent with previous studies that documented low levels of genetic variation in the species. However, genes previously implicated in local adaptation to abiotic environment in devils were enriched for differentially expressed genes. Additionally, three modules of co-expressed genes were significantly associated with either population of origin or sex.

Conserving adaptive potential: lessons from Tasmanian devils and their transmissible cancer

Maintenance of adaptive genetic variation has long been a goal of management of natural populations, but only recently have genomic tools allowed identification of specific loci associated with fitness-related traits in species of conservation concern. This raises the possibility of managing for genetic variation directly relevant to specific threats, such as those due to climate change or emerging infectious disease. Tasmanian devils (Sarcophilus harrisii) face the threat of a transmissible cancer, devil facial tumor disease (DFTD), that has decimated wild populations and led to intensive management efforts. Recent discoveries from genomic and modeling studies reveal how natural devil populations are responding to DFTD, and can inform management of both captive and wild devil populations. Notably, recent studies have documented genetic variation for disease-related traits and rapid evolution in response to DFTD, as well as potential mechanisms for disease resistance such as immune response and tumor regression in wild devils. Recent models predict dynamic persistence of devils with or without DFTD under a variety of modeling scenarios, although at much lower population densities than before DFTD emerged, contrary to previous predictions of extinction. As a result, current management that focuses on captive breeding and release for maintaining genome-wide genetic diversity or demographic supplementation of populations could have negative consequences. Translocations of captive devils into wild populations evolving with DFTD can cause outbreeding depression and/or increases in the force of infection and thereby the severity of the epidemic, and we argue that these risks outweigh any benefits of demographic supplementation in wild populations. We also argue that genetic variation at loci associated with DFTD should be monitored in both captive and wild populations, and that as our understanding of DFTD-related genetic variation improves, considering genetic management approaches to target this variation is warranted in developing conservation strategies for Tasmanian devils.

Large-effect loci affect survival in Tasmanian devils (Sarcophilus harrisii) infected with a transmissible cancer

Identifying the genetic architecture of complex phenotypes is a central goal of modern biology, particularly for disease-related traits. Genome-wide association methods are a classical approach for identifying the genomic basis of variation in disease phenotypes, but such analyses are particularly challenging in natural populations due to sample size difficulties. Extensive mark-recapture data, strong linkage disequilibrium and a lethal transmissible cancer make the Tasmanian devil (Sarcophilus harrisii) an ideal model for such an association study. We used a RAD-capture approach to genotype 624 devils at ~16,000 loci and then used association analyses to assess the heritability of three cancer-related phenotypes: infection case-control (where cases were infected devils and controls were devils that were never infected), age of first infection and survival following infection. The SNP array explained much of the phenotypic variance for female survival (>80%) and female case-control (>61%). We found that a few large-effect SNPs explained much of the variance for female survival (~5 SNPs explained >61% of the total variance), whereas more SNPs (~56) of smaller effect explained less of the variance for female case-control (~23% of the total variance). By contrast, these same SNPs did not account for a significant proportion of phenotypic variance in males, suggesting that the genetic bases of these traits and/or selection differ across sexes. Loci involved with cell adhesion and cell-cycle regulation underlay trait variation, suggesting that the devil immune system is rapidly evolving to recognize and potentially suppress cancer growth through these pathways. Overall, our study provided necessary data for genomics-based conservation and management in Tasmanian devils.

Lessons learnt from the Tasmanian devil facial tumour regarding immune function in cancer

Genetic and genomic technologies have facilitated a greater understanding of the Tasmanian devil immune system and the origins, evolution and spread of devil facial tumour disease (DFTD). DFTD is a contagious cancer that has caused significant declines in devil populations across Tasmania. Immune responses to DFTD are rarely detected, allowing the cancer to pass between individuals and proliferate unimpeded. Early immunosenscence in devils appears to decrease anti-tumour immunity in older animals compared to younger animals, which may increase susceptibility to DFTD and explain high DFTD prevalence in this age group. Devils also have extremely low major histocompatibility complex (MHC) diversity, and multiple alleles are shared with the tumour, lowering histocompatibility barriers which may have contributed to DFTD evolution. DFTD actively evades immune attack by down-regulating cell-surface MHC I molecules, making it effectively invisible to the immune system. Altered MHC I profiles should activate natural killer (NK) cell anti-tumour responses, but these are absent in DFTD infection. Recent immunisation and immunotherapy using modified DFTD cells has induced an anti-DFTD immune response and regression of DFTD in some devils. Knowledge gained from immune responses to a transmissible cancer in devils will ultimately reveal useful insights into immunity to cancer in humans and other species.

The Origins and Vulnerabilities of Two Transmissible Cancers in Tasmanian Devils

Transmissible cancers are clonal lineages that spread through populations via contagious cancer cells. Although rare in nature, two facial tumor clones affect Tasmanian devils. Here we perform comparative genetic and functional characterization of these lineages. The two cancers have similar patterns of mutation and show no evidence of exposure to exogenous mutagens or viruses. Genes encoding PDGF receptors have copy number gains and are present on extrachromosomal double minutes. Drug screening indicates causative roles for receptor tyrosine kinases and sensitivity to inhibitors of DNA repair. Y chromosome loss from a male clone infecting a female host suggests immunoediting. These results imply that Tasmanian devils may have inherent susceptibility to transmissible cancers and present a suite of therapeutic compounds for use in conservation.

MHC diversity and female age underpin reproductive success in an Australian icon; the Tasmanian Devil

Devil Facial Tumour Disease (DFTD), a highly contagious cancer, has decimated Tasmanian devil (Sarcophilus harrisii) numbers in the wild. To ensure its long-term survival, a captive breeding program was implemented but has not been as successful as envisaged at its launch in 2005. We therefore investigated the reproductive success of 65 captive devil pair combinations, of which 35 produced offspring (successful pairs) whereas the remaining 30 pairs, despite being observed mating, produced no offspring (unsuccessful pairs). The devils were screened at six MHC Class I-linked microsatellite loci. Our analyses revealed that younger females had a higher probability of being successful than older females. In the successful pairs we also observed a higher difference in total number of heterozygous loci, i.e. when one devil had a high total number of heterozygous loci, its partner had low numbers. Our results therefore suggest that devil reproductive success is subject to disruptive MHC selection, which to our knowledge has never been recorded in any vertebrate. In order to enhance the success of the captive breeding program the results from the present study show the importance of using young (2-year old) females as well as subjecting the devils to MHC genotyping.

Conservation implications of limited genetic diversity and population structure in Tasmanian devils (Sarcophilus harrisii)

Tasmanian devils face a combination of threats to persistence, including Devil Facial Tumor Disease (DFTD), an epidemic transmissible cancer. We used RAD sequencing to investigate genome-wide patterns of genetic diversity and geographic population structure. Consistent with previous results, we found very low genetic diversity in the species as a whole, and we detected two broad genetic clusters occupying the northwestern portion of the range, and the central and eastern portions. However, these two groups overlap across a broad geographic area, and differentiation between them is modest (F_ST = 0.1081). Our results refine the geographic extent of the zone of mixed ancestry and substructure within it, potentially informing management of genetic variation that existed in pre-diseased populations of the species. DFTD has spread across both genetic clusters, but recent evidence points to a genomic response to selection imposed by DFTD. Any allelic variation for resistance to DFTD may be able to spread across the devil population under selection by DFTD, and/or be present as standing variation in both genetic regions.

European Colonization, Not Polynesian Arrival, Impacted Population Size and Genetic Diversity in the Critically Endangered New Zealand Kākāpō

Island endemic species are often vulnerable to decline and extinction following human settlement, and the genetic study of historical museum specimens can be useful in understanding these processes. The kākāpō (Strigops habroptilus) is a critically endangered New Zealand parrot that was formerly widespread and abundant. It is well established that both Polynesian and European colonization of New Zealand impacted the native avifauna, but the timeframe and severity of impacts have differed depending on species. Here, we investigated the relative importance of the 2 waves of human settlement on kākāpō decline, using microsatellites and mitochondrial DNA (mtDNA) to characterize recent kākāpō genetic and demographic history. We analyzed samples from 49 contemporary individuals and 54 museum specimens dating from 1884 to 1985. Genetic diversity decreased significantly between historical and contemporary kākāpō, with a decline in mean number of microsatellite alleles from 6.15 to 3.08 and in number of mtDNA haplotypes from 17 to 3. Modeling of demographic history indicated a recent population bottleneck linked to the period of European colonization (approximately 5 generations ago) but did not support a major decline linked to Polynesian settlement. Effective population size estimates were also larger for historical than contemporary kākāpō. Our findings inform contemporary kākāpō management by indicating the timeframe and possible cause of the bottleneck, which has implications for the management of extant genetic diversity. We demonstrate the broader utility of a historical perspective in understanding causes of decline and managing extinction risk in contemporary endangered species.

Detecting Selection on Temporal and Spatial Scales: A Genomic Time-Series Assessment of Selective Responses to Devil Facial Tumor Disease

Detecting loci under selection is an important task in evolutionary biology. In conservation genetics detecting selection is key to investigating adaptation to the spread of infectious disease. Loci under selection can be detected on a spatial scale, accounting for differences in demographic history among populations, or on a temporal scale, tracing changes in allele frequencies over time. Here we use these two approaches to investigate selective responses to the spread of an infectious cancer—devil facial tumor disease (DFTD)—that since 1996 has ravaged the Tasmanian devil (Sarcophilus harrisii). Using time-series ‘restriction site associated DNA’ (RAD) markers from populations pre- and post DFTD arrival, and DFTD free populations, we infer loci under selection due to DFTD and investigate signatures of selection that are incongruent among methods, populations, and times. The lack of congruence among populations influenced by DFTD with respect to inferred loci under selection, and the direction of that selection, fail to implicate a consistent selective role for DFTD. Instead genetic drift is more likely driving the observed allele frequency changes over time. Our study illustrates the importance of applying methods with different performance optima e.g. accounting for population structure and background selection, and assessing congruence of the results.

The identification of immune genes in the milk transcriptome of the Tasmanian devil (Sarcophilus harrisii)

Tasmanian devil (Sarcophilus harrisii) pouch young, like other marsupials, are born underdeveloped and immunologically naïve, and are unable to mount an adaptive immune response. The mother’s milk provides nutrients for growth and development as well as providing passive immunity. To better understand immune response in this endangered species, we set out to characterise the genes involved in passive immunity by sequencing and annotating the transcriptome of a devil milk sample collected during mid-lactation. At mid-lactation we expect the young to have heightened immune responses, as they have emerged from the pouch, encountering new pathogens. A total of 233,660 transcripts were identified, including approximately 17,827 unique protein-coding genes and 846 immune genes. The most highly expressed transcripts were dominated by milk protein genes such as those encoding early lactation protein, late lactation proteins, α-lactalbumin, α-casein and β-casein. There were numerous highly expressed immune genes including lysozyme, whey acidic protein, ferritin and major histocompatibility complex I and II. Genes encoding immunoglobulins, antimicrobial peptides, chemokines and immune cell receptors were also identified. The array of immune genes identified in this study reflects the importance of the milk in providing immune protection to Tasmanian devil young and provides the first insight into Tasmanian devil milk.

A second transmissible cancer in Tasmanian devils

Clonally transmissible cancers are somatic cell lineages that are spread between individuals via the transfer of living cancer cells. There are only three known naturally occurring transmissible cancers, and these affect dogs, soft-shell clams, and Tasmanian devils, respectively. The Tasmanian devil transmissible facial cancer was first observed in 1996, and is threatening its host species with extinction. Until now, this disease has been consistently associated with a single aneuploid cancer cell lineage that we refer to as DFT1. Here we describe a second transmissible cancer, DFT2, in five devils located in southern Tasmania in 2014 and 2015. DFT2 causes facial tumors that are grossly indistinguishable but histologically distinct from those caused by DFT1. DFT2 bears no detectable cytogenetic similarity to DFT1 and carries a Y chromosome, which contrasts with the female origin of DFT1. DFT2 shows different alleles to both its hosts and DFT1 at microsatellite, structural variant, and major histocompatibility complex (MHC) loci, confirming that it is a second cancer that can be transmitted between devils as an allogeneic, MHC-discordant graft. These findings indicate that Tasmanian devils have spawned at least two distinct transmissible cancer lineages and suggest that transmissible cancers may arise more frequently in nature than previously considered. The discovery of DFT2 presents important challenges for the conservation of Tasmanian devils and raises the possibility that this species is particularly prone to the emergence of transmissible cancers. More generally, our findings highlight the potential for cancer cells to depart from their hosts and become dangerous transmissible pathogens.

Immunology of a Transmissible Cancer Spreading among Tasmanian Devils

Devil facial tumor disease (DFTD) is a transmissible cancer that has killed most of the Tasmanian devil (Sarcophilus harrissii) population. Since the first case appeared in the mid-1990s, it has spread relentlessly across the Tasmanian devil's geographic range. As Tasmanian devils only exist in Tasmania, Australia, DFTD has the potential to cause extinction of this species. The origin of DFTD was a Schwann cell from a female devil. The disease is transmitted when devils bite each other around the facial areas, a behavior synonymous with this species. Every devil that is 'infected' with DFTD dies from the cancer. Once the DFTD cells have been transmitted, they appear to develop into a cancer without inducing an immune response. The DFTD cancer cells avoid allogeneic recognition because they do not express MHC class I molecules on the cell surface. A reduced genetic diversity and the production of immunosuppressive cytokines may also contribute.

Lack of genetic diversity across diverse immune genes in an endangered mammal, the Tasmanian devil (Sarcophilus harrisii)

The Tasmanian devil (Sarcophilus harrisii) is threatened with extinction due to the spread of devil facial tumour disease. Polymorphisms in immune genes can provide adaptive potential to resist diseases. Previous studies in diversity at immune loci in wild species have almost exclusively focused on genes of the major histocompatibility complex (MHC); however, these genes only account for a fraction of immune gene diversity. Devils lack diversity at functionally important immunity loci, including MHC and Toll-like receptor genes. Whether there are polymorphisms at devil immune genes outside these two families is unknown. Here, we identify polymorphisms in a wide range of key immune genes, and develop assays to type single nucleotide polymorphisms (SNPs) within a subset of these genes. A total of 167 immune genes were examined, including cytokines, chemokines and natural killer cell receptors. Using genome-level data from ten devils, SNPs within coding regions, introns and 10 kb flanking genes of interest were identified. We found low polymorphism across 167 immune genes examined bioinformatically using whole-genome data. From this data, we developed long amplicon assays to target nine genes. These amplicons were sequenced in 29-220 devils and found to contain 78 SNPs, including eight SNPS within exons. Despite the extreme paucity of genetic diversity within these genes, signatures of balancing selection were exhibited by one chemokine gene, suggesting that remaining diversity may hold adaptive potential. The low functional diversity may leave devils highly vulnerable to infectious disease, and therefore, monitoring and preserving remaining diversity will be critical for the long-term management of this species. Examining genetic variation in diverse immune genes should be a priority for threatened wildlife species. This study can act as a model for broad-scale immunogenetic diversity analysis in threatened species.

Emerging infectious diseases of wildlife: a critical perspective

We review the literature to distinguish reports of vertebrate wildlife disease emergence with sufficient evidence, enabling a robust assessment of emergence drivers. For potentially emerging agents that cannot be confirmed, sufficient data on prior absence (or a prior difference in disease dynamics) are frequently lacking. Improved surveillance, particularly for neglected host taxa, geographical regions and infectious agents, would enable more effective management should emergence occur. Exposure to domestic sources of infection and human-assisted exposure to wild sources were identified as the two main drivers of emergence across host taxa; the domestic source was primary for fish while the wild source was primary for other taxa. There was generally insufficient evidence for major roles of other hypothesized drivers of emergence.

Diversity in the Toll-like receptor genes of the Tasmanian devil (Sarcophilus harrisii)

The Tasmanian devil is an endangered marsupial species that has survived several historical bottlenecks and now has low genetic diversity. Here we characterize the Toll-like receptor (TLR) genes and their diversity in the Tasmanian devil. TLRs are a key innate immune gene family found in all animals. Ten TLR genes were identified in the Tasmanian devil genome. Unusually low levels of diversity were found in 25 devils from across Tasmania. We found two alleles at TLR2, TLR3 and TLR6. The other seven genes were monomorphic. The insurance population, which safeguards the species from extinction, has successfully managed to capture all of these TLR alleles, but concerns remain for the long-term survival of this species.