skip to main content


The NSF Public Access Repository (NSF-PAR) system and access will be unavailable from 5:00 PM ET until 11:00 PM ET on Friday, June 21 due to maintenance. We apologize for the inconvenience.

Title: Spatial variation in gene expression of Tasmanian devil facial tumors despite minimal host transcriptomic response to infection
Abstract 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.  more » « less
Award ID(s):
Author(s) / Creator(s):
; ; ; ; ; ; ; ; ; ;
Date Published:
Journal Name:
BMC Genomics
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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.

    more » « less
  2. Abstract

    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, whereasDFTDfaded 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 theDFTD‐induced devil extinction, with crucial management ramifications. Accounting for individual‐level disease progression and the long‐term outcome of devil–DFTDinteractions 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.

    more » « less
  3. Abstract

    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.

    more » « less
  4. Abstract

    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.

    more » « less
  5. Abstract

    A plethora of intrinsic and environmental factors have been shown to influence the length of telomeres, the protector of chromosome ends. Despite the growing interest in infection–telomere interactions, there is very limited knowledge on how transmissible cancers influence telomere maintenance. An emblematic example of transmissible cancer occurs in the Tasmanian devil (Sarcophilus harrisii), whose populations have been dramatically reduced by infectious cancer cells. To investigate associations between telomere dynamics and the transmissible cancer, we used longitudinal data from a Tasmanian devil population that has been exposed to the disease for over 15 years. We detected substantial temporal variation in individual telomere length (TL), and a positive significant association between TL and age, as well as a marginally significant trend for devils with devil facial tumour disease (DFTD) having longer telomeres. A proportional hazard analysis yielded no significant effect of TL on the development of DFTD. Like previous studies, we show the complexity that TL dynamics may exhibit across the lifetime of organisms. Our work highlights the importance of long‐term longitudinal sampling for understanding the effects of wildlife diseases on TL.

    more » « less