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.

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.

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.

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.

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.

Quantifying 25 years of disease-caused declines in Tasmanian devil populations: host density drives spatial pathogen spread.

Infectious diseases are strong drivers of wildlife population dynamics, however, empirical analyses from the early stages of pathogen emergence are rare. Tasmanian devil facial tumour disease (DFTD), discovered in 1996, provides the opportunity to study an epizootic from its inception. We use a pattern-oriented diffusion simulation to model the spatial spread of DFTD across the species' range and quantify population effects by jointly modelling multiple streams of data spanning 35 years. We estimate the wild devil population peaked at 53 000 in 1996, less than half of previous estimates. DFTD spread rapidly through high-density areas, with spread velocity slowing in areas of low host densities. By 2020, DFTD occupied >90% of the species' range, causing 82% declines in local densities and reducing the total population to 16 900. Encouragingly, our model forecasts the population decline should level-off within the next decade, supporting conservation management focused on facilitating evolution of resistance and tolerance.

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.

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.

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.

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.

Spontaneous Tumor Regression in Tasmanian Devils Associated with RASL11A Activation.

Spontaneous tumor regression has been documented in a small proportion of human cancer patients, but the specific mechanisms underlying tumor regression without treatment are not well understood. Tasmanian devils are threatened with extinction from a transmissible cancer due to universal susceptibility and a near 100% case fatality rate. In over 10,000 cases, <20 instances of natural tumor regression have been detected. Previous work in this system has focused on Tasmanian devil genetic variation associated with the regression phenotype. Here, we used comparative and functional genomics to identify tumor genetic variation associated with tumor regression. We show that a single point mutation in the 5' untranslated region of the putative tumor suppressor RASL11A significantly contributes to tumor regression. RASL11A was expressed in regressed tumors but silenced in wild-type, nonregressed tumors, consistent with RASL11A downregulation in human cancers. Induced RASL11A expression significantly reduced tumor cell proliferation in vitro The RAS pathway is frequently altered in human cancers, and RASL11A activation may provide a therapeutic treatment option for Tasmanian devils as well as a general mechanism for tumor inhibition.

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.

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.

Tracing the rise of malignant cell lines: Distribution, epidemiology and evolutionary interactions of two transmissible cancers in Tasmanian devils.

Emerging infectious diseases are rising globally and understanding host-pathogen interactions during the initial stages of disease emergence is essential for assessing potential evolutionary dynamics and designing novel management strategies. Tasmanian devils (Sarcophilus harrisii) are endangered due to a transmissible cancer-devil facial tumour disease (DFTD)-that since its emergence in the 1990s, has affected most populations throughout Tasmania. Recent studies suggest that devils are adapting to the DFTD epidemic and that disease-induced extinction is unlikely. However, in 2014, a second and independently evolved transmissible cancer-devil facial tumour 2 (DFT2)-was discovered at the d'Entrecasteaux peninsula, in south-east Tasmania, suggesting that the species is prone to transmissible cancers. To date, there is little information about the distribution, epidemiology and effects of DFT2 and its interaction with DFTD. Here, we use data from monitoring surveys and roadkills found within and adjacent to the d'Entrecasteaux peninsula to determine the distribution of both cancers and to compare their epidemiological patterns. Since 2012, a total of 51 DFTD tumours have been confirmed among 26 individuals inside the peninsula and its surroundings, while 40 DFT2 tumours have been confirmed among 23 individuals, and two individuals co-infected with both tumours. All devils with DFT2 were found within the d'Entrecasteaux peninsula, suggesting that this new transmissible cancer is geographically confined to this area. We found significant differences in tumour bodily location in DFTD and DFT2, with non-facial tumours more commonly found in DFT2. There was a significant sex bias in DFT2, with most cases reported in males, suggesting that since DFT2 originated from a male host, females might be less susceptible to this cancer. We discuss the implications of our results for understanding the epidemiological and evolutionary interactions of these two contemporary transmissible cancers and evaluating the effectiveness of potential management strategies.

Individual and temporal variation in pathogen load predicts long-term impacts of an emerging infectious disease.

Emerging infectious diseases increasingly threaten wildlife populations. Most studies focus on managing short-term epidemic properties, such as controlling early outbreaks. Predicting long-term endemic characteristics with limited retrospective data is more challenging. We used individual-based modeling informed by individual variation in pathogen load and transmissibility to predict long-term impacts of a lethal, transmissible cancer on Tasmanian devil (Sarcophilus harrisii) populations. For this, we employed approximate Bayesian computation to identify model scenarios that best matched known epidemiological and demographic system properties derived from 10 yr of data after disease emergence, enabling us to forecast future system dynamics. We show that the dramatic devil population declines observed thus far are likely attributable to transient dynamics (initial dynamics after disease emergence). Only 21% of matching scenarios led to devil extinction within 100 yr following devil facial tumor disease (DFTD) introduction, whereas DFTD faded out in 57% of simulations. In the remaining 22% of simulations, disease and host coexisted for at least 100 yr, usually with long-period oscillations. Our findings show that pathogen extirpation or host-pathogen coexistence are much more likely than the DFTD-induced devil extinction, with crucial management ramifications. Accounting for individual-level disease progression and the long-term outcome of devil-DFTD interactions at the population-level, our findings suggest that immediate management interventions are unlikely to be necessary to ensure the persistence of Tasmanian devil populations. This is because strong population declines of devils after disease emergence do not necessarily translate into long-term population declines at equilibria. Our modeling approach is widely applicable to other host-pathogen systems to predict disease impact beyond transient dynamics.

Sex bias in ability to cope with cancer: Tasmanian devils and facial tumour disease.

Knowledge of the ecological dynamics between hosts and pathogens during the initial stages of disease emergence is crucial to understanding the potential for evolution of new interspecific interactions. Tasmanian devil (Sarcophilus harrisii) populations have declined precipitously owing to infection by a transmissible cancer (devil facial tumour disease, DFTD) that emerged approximately 20 years ago. Since the emergence of DFTD, and as the disease spreads across Tasmania, the number of devils has dropped up to 90% across 80% of the species's distributional range. As a result, the disease is expected to act as a strong selective force on hosts to develop mechanisms of tolerance and/or resistance to the infection. We assessed the ability of infected devils to cope with infection, which translates into host tolerance to the cancer, by using the reaction norm of the individual body condition by tumour burden. We found that body condition of infected hosts is negatively affected by cancer progression. Males and females presented significant differences in their tolerance levels to infection, with males suffering declines of up to 25% of their body condition, in contrast to less than 5% in females. Sex-related differences in tolerance to cancer progression may select for changes in life-history strategies of the host and could also alter the selective environment for the tumours.

The Genomic Basis of Tumor Regression in Tasmanian Devils (Sarcophilus harrisii).

Understanding the genetic basis of disease-related phenotypes, such as cancer susceptibility, is crucial for the advancement of personalized medicine. Although most cancers are somatic in origin, a small number of transmissible cancers have been documented. Two such cancers have emerged in the Tasmanian devil (Sarcophilus harrisii) and now threaten the species with extinction. Recently, cases of natural tumor regression in Tasmanian devils infected with the clonally contagious cancer have been detected. We used whole-genome sequencing and FST-based approaches to identify the genetic basis of tumor regression by comparing the genomes of seven individuals that underwent tumor regression with those of three infected individuals that did not. We found three highly differentiated candidate genomic regions containing several genes related to immune response and/or cancer risk, indicating that the genomic basis of tumor regression was polygenic. Within these genomic regions, we identified putative regulatory variation in candidate genes but no nonsynonymous variation, suggesting that natural tumor regression may be driven, at least in part, by differential host expression of key loci. Comparative oncology can provide insight into the genetic basis of cancer risk, tumor development, and the pathogenicity of cancer, particularly due to our limited ability to monitor natural, untreated tumor progression in human patients. Our results support the hypothesis that host immune response is necessary for triggering tumor regression, providing candidate genes that may translate to novel treatments in human and nonhuman cancers.

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.