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 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 two 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.

Adaptive potential in the face of a transmissible cancer in Tasmanian devils.

Emerging infectious diseases (EIDs) not only cause catastrophic declines in wildlife populations but also generate selective pressures that may result in rapid evolutionary responses. One such EID is devil facial tumour disease (DFTD) in the Tasmanian devil. DFTD is almost always fatal and has reduced the average lifespan of individuals by around 2 years, likely causing strong selection for traits that reduce susceptibility to the disease, but population decline has also left Tasmanian devils vulnerable to inbreeding depression. We analysed 22 years of data from an ongoing study of a population of Tasmanian devils on Freycinet Peninsula, Tasmania, to (1) identify whether DFTD may be causing selection on body size, by estimating phenotypic and genetic correlations between DFTD and size traits, (2) estimate the additive genetic variance of susceptibility to DFTD, and (3) investigate whether size traits or susceptibility to DFTD were under inbreeding depression. We found a positive phenotypic relationship between head width and susceptibility to DFTD, but this was not underpinned by a genetic correlation. Conversely, we found a negative phenotypic relationship between body weight and susceptibility to DFTD, and there was evidence for a negative genetic correlation between susceptibility to DFTD and body weight. There was additive genetic variance in susceptibility to DFTD, head width and body weight, but there was no evidence for inbreeding depression in any of these traits. These results suggest that Tasmanian devils have the potential to respond adaptively to DFTD, although the realised evolutionary response will critically further depend on the evolution of DFTD itself.

Do Tasmanian devil declines impact ecosystem function?

Tasmanian eucalypt forests are among the most carbon-dense in the world, but projected climate change could destabilize this critical carbon sink. While the impact of abiotic factors on forest ecosystem carbon dynamics have received considerable attention, biotic factors such as the input of animal scat are less understood. Tasmanian devils (Sarcophilus harrisii)-an osteophageous scavenger that can ingest and solubilize nutrients locked in bone material-may subsidize plant and microbial productivity by concentrating bioavailable nutrients (e.g., nitrogen and phosphorus) in scat latrines. However, dramatic declines in devil population densities, driven by the spread of a transmissible cancer, may have underappreciated consequences for soil organic carbon (SOC) storage and forest productivity by altering nutrient cycling. Here, we fuse experimental data and modeling to quantify and predict future changes to forest productivity and SOC under various climate and scat-quality futures. We find that devil scat significantly increases concentrations of nitrogen, ammonium, phosphorus, and phosphate in the soil and shifts soil microbial communities toward those dominated by r-selected (e.g., fast-growing) phyla. Further, under expected increases in temperature and changes in precipitation, devil scat inputs are projected to increase above- and below-ground net primary productivity and microbial biomass carbon through 2100. In contrast, when devil scat is replaced by lower-quality scat (e.g., from non-osteophageous scavengers and herbivores), forest carbon pools are likely to increase more slowly, or in some cases, decline. Together, our results suggest often overlooked biotic factors will interact with climate change to drive current and future carbon pool dynamics in Tasmanian forests.

Intergenomic signatures of coevolution between Tasmanian devils and an infectious cancer.

Coevolution is common and frequently governs host-pathogen interaction outcomes. Phenotypes underlying these interactions often manifest as the combined products of the genomes of interacting species, yet traditional quantitative trait mapping approaches ignore these intergenomic interactions. Devil facial tumor disease (DFTD), an infectious cancer afflicting Tasmanian devils (Sarcophilus harrisii), has decimated devil populations due to universal host susceptibility and a fatality rate approaching 100%. Here, we used a recently developed joint genome-wide association study (i.e., co-GWAS) approach, 15 y of mark-recapture data, and 960 genomes to identify intergenomic signatures of coevolution between devils and DFTD. Using a traditional GWA approach, we found that both devil and DFTD genomes explained a substantial proportion of variance in how quickly susceptible devils became infected, although genomic architectures differed across devils and DFTD; the devil genome had fewer loci of large effect whereas the DFTD genome had a more polygenic architecture. Using a co-GWA approach, devil-DFTD intergenomic interactions explained ~3× more variation in how quickly susceptible devils became infected than either genome alone, and the top genotype-by-genotype interactions were significantly enriched for cancer genes and signatures of selection. A devil regulatory mutation was associated with differential expression of a candidate cancer gene and showed putative allele matching effects with two DFTD coding sequence variants. Our results highlight the need to account for intergenomic interactions when investigating host-pathogen (co)evolution and emphasize the importance of such interactions when considering devil management strategies.

Complex associations between cancer progression and immune gene expression reveals early influence of transmissible cancer on Tasmanian devils.

The world's largest extant carnivorous marsupial, the Tasmanian devil, is challenged by Devil Facial Tumor Disease (DFTD), a fatal, clonally transmitted cancer. In two decades, DFTD has spread across 95% of the species distributional range. A previous study has shown that factors such as season, geographic location, and infection with DFTD can impact the expression of immune genes in Tasmanian devils. To date, no study has investigated within-individual immune gene expression changes prior to and throughout the course of DFTD infection. To explore possible changes in immune response, we investigated four locations across Tasmania that differed in DFTD exposure history, ranging between 2 and >30 years. Our study demonstrated considerable complexity in the immune responses to DFTD. The same factors (sex, age, season, location and DFTD infection) affected immune gene expression both across and within devils, although seasonal and location specific variations were diminished in DFTD affected devils. We also found that expression of both adaptive and innate immune genes starts to alter early in DFTD infection and continues to change as DFTD progresses. A novel finding was that the lower expression of immune genes MHC-II, NKG2D and CD8 may predict susceptibility to earlier DFTD infection. A case study of a single devil with regressed tumor showed opposite/contrasting immune gene expression patterns compared to the general trends observed across devils with DFTD infection. Our study highlights the complexity of DFTD's interactions with the host immune system and the need for long-term studies to fully understand how DFTD alters the evolutionary trajectory of devil immunity.

Disease-driven top predator decline affects mesopredator population genomic structure.

Top predator declines are pervasive and often have dramatic effects on ecological communities via changes in food web dynamics, but their evolutionary consequences are virtually unknown. Tasmania's top terrestrial predator, the Tasmanian devil, is declining due to a lethal transmissible cancer. Spotted-tailed quolls benefit via mesopredator release, and they alter their behaviour and resource use concomitant with devil declines and increased disease duration. Here, using a landscape community genomics framework to identify environmental drivers of population genomic structure and signatures of selection, we show that these biotic factors are consistently among the top variables explaining genomic structure of the quoll. Landscape resistance negatively correlates with devil density, suggesting that devil declines will increase quoll genetic subdivision over time, despite no change in quoll densities detected by camera trap studies. Devil density also contributes to signatures of selection in the quoll genome, including genes associated with muscle development and locomotion. Our results provide some of the first evidence of the evolutionary impacts of competition between a top predator and a mesopredator species in the context of a trophic cascade. As top predator declines are increasing globally, our framework can serve as a model for future studies of evolutionary impacts of altered ecological interactions.

Over-eruption in marsupial carnivore teeth: compensation for a constraint.

Pronounced over-eruption of the canine teeth, causing the cervical enamel margin to extend beyond the alveolar bone and exposing the root, occurs with age and growth in Australian marsupial carnivores, much more than in eco-morphologically equivalent placental carnivores. Suppression of functional tooth replacement is characteristic of marsupials, where most placentals have the primitive diphyodont pattern of two generations of incisor, canine and premolar teeth. Canine and molar tooth dimensions of four species of marsupial carnivores (thylacine Thylacinus cynocephalus, Tasmanian devil Sarcophilus harrisii and two quolls Dasyurus spp.) and canine dimensions of seven eco-morphologically equivalent placental carnivore species were measured from museum specimens. Canine dimensions were measured in a time series on live wild-living individual devils and quolls. The canine teeth and to a lesser extent the molar teeth of marsupial carnivores continue to erupt through life, resulting in a net increase in tooth height and diameter, a phenomenon not evident in placental carnivores. Potential mechanisms causing over-eruption include tooth wear and gradual release of occlusal pressure as the individual grows. Over-eruption in marsupial carnivores may be a compensatory response for tooth size limits imposed by monophyodont tooth replacement, ensuring that animal's teeth are scaled to jaw size from juvenile to adulthood.

The tumour is in the detail: Local phylogenetic, population and epidemiological dynamics of a transmissible cancer in Tasmanian devils.

Infectious diseases are a major threat for biodiversity conservation and can exert strong influence on wildlife population dynamics. Understanding the mechanisms driving infection rates and epidemic outcomes requires empirical data on the evolutionary trajectory of pathogens and host selective processes. Phylodynamics is a robust framework to understand the interaction of pathogen evolutionary processes with epidemiological dynamics, providing a powerful tool to evaluate disease control strategies. Tasmanian devils have been threatened by a fatal transmissible cancer, devil facial tumour disease (DFTD), for more than two decades. Here we employ a phylodynamic approach using tumour mitochondrial genomes to assess the role of tumour genetic diversity in epidemiological and population dynamics in a devil population subject to 12 years of intensive monitoring, since the beginning of the epidemic outbreak. DFTD molecular clock estimates of disease introduction mirrored observed estimates in the field, and DFTD genetic diversity was positively correlated with estimates of devil population size. However, prevalence and force of infection were the lowest when devil population size and tumour genetic diversity was the highest. This could be due to either differential virulence or transmissibility in tumour lineages or the development of host defence strategies against infection. Our results support the view that evolutionary processes and epidemiological trade-offs can drive host-pathogen coexistence, even when disease-induced mortality is extremely high. We highlight the importance of integrating pathogen and population evolutionary interactions to better understand long-term epidemic dynamics and evaluating disease control strategies.

The evolution of two transmissible cancers in Tasmanian devils.

Tasmanian devils have spawned two transmissible cancer lineages, named devil facial tumor 1 (DFT1) and devil facial tumor 2 (DFT2). We investigated the genetic diversity and evolution of these clones by analyzing 78 DFT1 and 41 DFT2 genomes relative to a newly assembled, chromosome-level reference. Time-resolved phylogenetic trees reveal that DFT1 first emerged in 1986 (1982 to 1989) and DFT2 in 2011 (2009 to 2012). Subclone analysis documents transmission of heterogeneous cell populations. DFT2 has faster mutation rates than DFT1 across all variant classes, including substitutions, indels, rearrangements, transposable element insertions, and copy number alterations, and we identify a hypermutated DFT1 lineage with defective DNA mismatch repair. Several loci show plausible evidence of positive selection in DFT1 or DFT2, including loss of chromosome Y and inactivation of MGA, but none are common to both cancers. This study reveals the parallel long-term evolution of two transmissible cancers inhabiting a common niche in Tasmanian devils.

Cathelicidin-3 Associated With Serum Extracellular Vesicles Enables Early Diagnosis of a Transmissible Cancer.

The identification of practical early diagnostic biomarkers is a cornerstone of improved prevention and treatment of cancers. Such a case is devil facial tumor disease (DFTD), a highly lethal transmissible cancer afflicting virtually an entire species, the Tasmanian devil (Sarcophilus harrisii). Despite a latent period that can exceed one year, to date DFTD diagnosis requires visual identification of tumor lesions. To enable earlier diagnosis, which is essential for the implementation of effective conservation strategies, we analyzed the extracellular vesicle (EV) proteome of 87 Tasmanian devil serum samples using data-independent acquisition mass spectrometry approaches. The antimicrobial peptide cathelicidin-3 (CATH3), released by innate immune cells, was enriched in serum EV samples of both devils with clinical DFTD (87.9% sensitivity and 94.1% specificity) and devils with latent infection (i.e., collected while overtly healthy, but 3-6 months before subsequent DFTD diagnosis; 93.8% sensitivity and 94.1% specificity). Although high expression of antimicrobial peptides has been mostly related to inflammatory diseases, our results suggest that they can be also used as accurate cancer biomarkers, suggesting a mechanistic role in tumorous processes. This EV-based approach to biomarker discovery is directly applicable to improving understanding and diagnosis of a broad range of diseases in other species, and these findings directly enhance the capacity of conservation strategies to ensure the viability of the imperiled Tasmanian devil population.

Extracellular vesicle proteomes of two transmissible cancers of Tasmanian devils reveal tenascin-C as a serum-based differential diagnostic biomarker.

The iconic Tasmanian devil (Sarcophilus harrisii) is endangered due to the transmissible cancer Devil Facial Tumour Disease (DFTD), of which there are two genetically independent subtypes (DFT1 and DFT2). While DFT1 and DFT2 can be differentially diagnosed using tumour biopsies, there is an urgent need to develop less-invasive biomarkers that can detect DFTD and distinguish between subtypes. Extracellular vesicles (EVs), the nano-sized membrane-enclosed vesicles present in most biofluids, represent a valuable resource for biomarker discovery. Here, we characterized the proteome of EVs from cultured DFTD cells using data-independent acquisition-mass spectrometry and an in-house spectral library of > 1500 proteins. EVs from both DFT1 and DFT2 cell lines expressed higher levels of proteins associated with focal adhesion functions. Furthermore, hallmark proteins of epithelial-mesenchymal transition were enriched in DFT2 EVs relative to DFT1 EVs. These findings were validated in EVs derived from serum samples, revealing that the mesenchymal marker tenascin-C was also enriched in EVs derived from the serum of devils infected with DFT2 relative to those infected with DFT1 and healthy controls. This first EV-based investigation of DFTD increases our understanding of the cancers' EVs and their possible involvement in DFTD progression, such as metastasis. Finally, we demonstrated the potential of EVs to differentiate between DFT1 and DFT2, highlighting their potential use as less-invasive liquid biopsies for the Tasmanian devil.

Spatial variation in gene expression of Tasmanian devil facial tumors despite minimal host transcriptomic response to infection.

BACKGROUND: Transmissible cancers lie at the intersection of oncology and infectious disease, two traditionally divergent fields for which gene expression studies are particularly useful for identifying the molecular basis of phenotypic variation. In oncology, transcriptomics studies, which characterize the expression of thousands of genes, have identified processes leading to heterogeneity in cancer phenotypes and individual prognoses. More generally, transcriptomics studies of infectious diseases characterize interactions between host, pathogen, and environment to better predict population-level outcomes. Tasmanian devils have been impacted dramatically by a transmissible cancer (devil facial tumor disease; DFTD) that has led to widespread population declines. Despite initial predictions of extinction, populations have persisted at low levels, due in part to heterogeneity in host responses, particularly between sexes. However, the processes underlying this variation remain unknown. RESULTS: We sequenced transcriptomes from healthy and DFTD-infected devils, as well as DFTD tumors, to characterize host responses to DFTD infection, identify differing host-tumor molecular interactions between sexes, and investigate the extent to which tumor gene expression varies among host populations. We found minimal variation in gene expression of devil lip tissues, either with respect to DFTD infection status or sex. However, 4088 genes were differentially expressed in tumors among our sampling localities. Pathways that were up- or downregulated in DFTD tumors relative to normal tissues exhibited the same patterns of expression with greater intensity in tumors from localities that experienced DFTD for longer. No mRNA sequence variants were associated with expression variation. CONCLUSIONS: Expression variation among localities may reflect morphological differences in tumors that alter ratios of normal-to-tumor cells within biopsies. Phenotypic variation in tumors may arise from environmental variation or differences in host immune response that were undetectable in lip biopsies, potentially reflecting variation in host-tumor coevolutionary relationships among sites that differ in the time since DFTD arrival.

Isotopic niche variation in Tasmanian devils Sarcophilus harrisii with progression of devil facial tumor disease.

Devil facial tumor disease (DFTD) is a transmissible cancer affecting Tasmanian devils Sarcophilus harrisii. The disease has caused severe population declines and is associated with demographic and behavioral changes, including earlier breeding, younger age structures, and reduced dispersal and social interactions. Devils are generally solitary, but social encounters are commonplace when feeding upon large carcasses. DFTD tumors can disfigure the jaw and mouth and so diseased individuals might alter their diets to enable ingestion of alternative foods, to avoid conspecific interactions, or to reduce competition. Using stable isotope analysis (δ13C and δ15N) of whiskers, we tested whether DFTD progression, measured as tumor volume, affected the isotope ratios and isotopic niches of 94 infected Tasmanian devils from six sites in Tasmania, comprising four eucalypt plantations, an area of smallholdings and a national park. Then, using tissue from 10 devils sampled before and after detection of tumors and 8 devils where no tumors were detected, we examined whether mean and standard deviation of δ13C and δ15N of the same individuals changed between healthy and diseased states. δ13C and δ15N values were generally not related to tumor volume in infected devils, though at one site, Freycinet National Park, δ15N values increased significantly as tumor volume increased. Infection with DFTD was not associated with significant changes in the mean or standard deviation of δ13C and δ15N values in individual devils sampled before and after detection of tumors. Our analysis suggests that devils tend to maintain their isotopic niche in the face of DFTD infection and progression, except where ecological conditions facilitate a shift in diets and feeding behaviors, demonstrating that ecological context, alongside disease severity, can modulate the behavioral responses of Tasmanian devils to DFTD.

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.

Infectious disease and sickness behaviour: tumour progression affects interaction patterns and social network structure in wild Tasmanian devils.

Infectious diseases, including transmissible cancers, can have a broad range of impacts on host behaviour, particularly in the latter stages of disease progression. However, the difficulty of early diagnoses makes the study of behavioural influences of disease in wild animals a challenging task. Tasmanian devils (Sarcophilus harrisii) are affected by a transmissible cancer, devil facial tumour disease (DFTD), in which tumours are externally visible as they progress. Using telemetry and mark-recapture datasets, we quantify the impacts of cancer progression on the behaviour of wild devils by assessing how interaction patterns within the social network of a population change with increasing tumour load. The progression of DFTD negatively influences devils' likelihood of interaction within their network. Infected devils were more active within their network late in the mating season, a pattern with repercussions for DFTD transmission. Our study provides a rare opportunity to quantify and understand the behavioural feedbacks of disease in wildlife and how they may affect transmission and population dynamics in general.

A transmissible cancer shifts from emergence to endemism in Tasmanian devils.

Emerging infectious diseases pose one of the greatest threats to human health and biodiversity. Phylodynamics is often used to infer epidemiological parameters essential for guiding intervention strategies for human viruses such as severe acute respiratory syndrome coronavirus 2 (SARS-Cov-2). Here, we applied phylodynamics to elucidate the epidemiological dynamics of Tasmanian devil facial tumor disease (DFTD), a fatal, transmissible cancer with a genome thousands of times larger than that of any virus. Despite prior predictions of devil extinction, transmission rates have declined precipitously from ~3.5 secondary infections per infected individual to ~1 at present. Thus, DFTD appears to be transitioning from emergence to endemism, lending hope for the continued survival of the endangered Tasmanian devil. More generally, our study demonstrates a new phylodynamic analytical framework that can be applied to virtually any pathogen.

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.

Blood Parasites in Endangered Wildlife-Trypanosomes Discovered During a Survey of Haemoprotozoa from the Tasmanian Devil.

The impact of emerging infectious diseases is increasingly recognised as a major threat to wildlife. Wild populations of the endangered Tasmanian devil, Sarcophilus harrisii, are experiencing devastating losses from a novel transmissible cancer, devil facial tumour disease (DFTD); however, despite the rapid decline of this species, there is currently no information on the presence of haemoprotozoan parasites. In the present study, 95 Tasmanian devil blood samples were collected from four populations in Tasmania, Australia, which underwent molecular screening to detect four major groups of haemoprotozoa: (i) trypanosomes, (ii) piroplasms, (iii) Hepatozoon, and (iv) haemosporidia. Sequence results revealed Trypanosoma infections in 32/95 individuals. Trypanosoma copemani was identified in 10 Tasmanian devils from three sites and a second Trypanosoma sp. was identified in 22 individuals that were grouped within the poorly described T. cyclops clade. A single blood sample was positive for Babesia sp., which most closely matched Babesia lohae. No other blood protozoan parasite DNA was detected. This study provides the first insight into haemoprotozoa from the Tasmanian devil and the first identification of Trypanosoma and Babesia in this carnivorous marsupial.

Telomere Length is a Susceptibility Marker for Tasmanian Devil Facial Tumor Disease.

Telomeres protect chromosomes from degradation during cellular replication. In humans, it is well-documented that excessive telomere degradation is one mechanism by which cells can become cancerous. Increasing evidence from wildlife studies suggests that telomere length is positively correlated with survival and health and negatively correlated with disease infection intensity. The recently emerged devil facial tumor disease (DFTD) has led to dramatic and rapid population declines of the Tasmanian devil throughout its geographic range. Here, we tested the hypothesis that susceptibility to DFTD is negatively correlated with telomere length in devils across three populations with different infection histories. Our findings suggest telomere length is correlated with DFTD resistance in three ways. First, devils from a population with the slowest recorded increase in DFTD prevalence (West Pencil Pine) have significantly longer telomeres than those from two populations with rapid and exponential increases in prevalence (Freycinet and Narawantapu). Second, using extensive mark-recapture data obtained from a long-term demographic study, we found that individuals with relatively long telomeres tend to be infected at a significantly later age than those with shorter telomeres. Third, a hazard model showed devils with longer telomeres tended to become infected at a lower rate than those with shorter telomeres. This research provides a rare study of telomere length variation and its association with disease in a wildlife population. Our results suggest that telomere length may be a reliable marker of susceptibility to DFTD and assist with future management of this endangered species.