A transmissible cancer shifts from emergence to endemism in Tasmanian devils

Emergence to endemism The emergence of a devastating transmissible facial cancer among Tasmanian devils over the past few decades has caused substantial concern for their future because these animals are already threatened by a regional distribution and other stressors. Little is known about the overall history and trajectory of this disease. Patton et al. used an epidemiological phylodynamic approach to reveal the pattern of disease emergence and spread. They found that low Tasmanian devil densities appear to be contributing to slower disease growth and spread, which is good news for Tasmanian devil persistence and suggests that care should be taken when considering options for increasing devil populations. Science, this issue p. eabb9772 Tasmanian devil facial tumor disease is now endemic in the wild, which may be good news for devil survival overall. INTRODUCTION Emerging infectious diseases pose one of the greatest threats to human health and biodiversity. Phylodynamics is an effective tool for inferring epidemiological parameters to guide intervention strategies, particularly for human viruses such as severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2). However, phylodynamic analysis has historically been limited to the study of rapidly evolving viruses and, in rare cases, bacteria. Nonetheless, application of phylodynamics to nonviral pathogens has immense potential, such as for predicting disease spread and informing the management of wildlife diseases. We conducted a phylodynamics analysis of devil facial tumor disease (DFTD), a transmissible cancer that has spread across nearly the entire geographic range of Tasmanian devils and threatens the species with extinction. DFTD is transmitted as an allograft through biting during common social interactions, susceptibility is nearly universal, and case fatality rates approach 100%. The goals of our study were to (i) characterize the geographic spread of DFTD, (ii) identify whether there are different circulating tumor lineages, and (iii) quantify rates of transmission among lineages. RATIONALE In principle, phylodynamics should be readily extended to the study of slowly evolving pathogens with large genomes through careful interrogation of genes to identify those that are measurably evolving. By testing individual genes for a clocklike signal, these genes may then be used for phylodynamic analysis. We demonstrate this proof of concept in DFTD. RESULTS We screened >11,000 genes across the DFTD genome, identifying 28 that exhibited a strong, clocklike signal, and performed the first phylodynamic analysis of a genome larger than a bacterium. We demonstrate here, contrary to field observations, that DFTD spread omnidirectionally throughout the epizootic, leaving little signal of geographic structuring of tumor lineages across Tasmania. Despite predictions of devil extinction, we found that the effective reproduction number (RE), a summary of the rate at which disease spreads, has declined precipitously after the initial epidemic spread of DFTD. Specifically, RE peaked at a high of ~3.5 shortly after the discovery of DFTD in 1996 and is now ~1 in both extant tumor lineages. This is consistent with a shift from emergence to endemism. Except for a single gene, we found little evidence for convergent molecular evolution among tumor lineages. CONCLUSION We have demonstrated that phylodynamics can be applied to virtually any pathogen. In doing so, we show that through careful interrogation of the pathogen genome, a measurably evolving set of genes can be identified to characterize epidemiological dynamics of nonviral pathogens with large genomes. By applying this approach to DFTD, we have shown that the disease appears to be transitioning from emergence to endemism. Consistent with recent models, our inference that RE ~1 predicts that coexistence between devils and DFTD is a more likely outcome than devil extinction. Therefore, our findings present cautious optimism for the continued survival of the iconic Tasmanian devil but emphasize the need for evolutionarily informed conservation management to ensure their persistence. Tasmanian devils and their transmissible cancer. Healthy (top) and DFTD-infected (bottom) Tasmanian devils. Photos: David G. Hamilton (top), Alexandra K. Fraik (bottom). 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.

[1]  Austin H. Patton,et al.  Historical and contemporary signatures of selection in response to transmissible cancer in the Tasmanian Devil (Sarcophilus harrisii) , 2020, bioRxiv.

[2]  Austin H. Patton,et al.  Comparative landscape genetics reveals differential effects of environment on host and pathogen genetic structure in Tasmanian devils (Sarcophilus harrisii) and their transmissible tumour , 2020, Molecular ecology.

[3]  Patrick J. Paddison,et al.  Spontaneous Tumor Regression in Tasmanian Devils Associated with RASL11A Activation , 2020, Genetics.

[4]  Menna E. Jones,et al.  Changes in spatial organization following an acute epizootic: Tasmanian devils and their transmissible cancer. , 2020, Global ecology and conservation.

[5]  Prasenjit Maity,et al.  COVID-19 outbreak: Migration, effects on society, global environment and prevention , 2020, Science of The Total Environment.

[6]  Yonatan H. Grad,et al.  Projecting the transmission dynamics of SARS-CoV-2 through the postpandemic period , 2020, Science.

[7]  P. Niu,et al.  Genome Composition and Divergence of the Novel Coronavirus (2019-nCoV) Originating in China , 2020, Cell Host & Microbe.

[8]  A. B. Lyons,et al.  An oral bait vaccination approach for the Tasmanian devil facial tumor diseases , 2020, Expert review of vaccines.

[9]  Austin H. Patton,et al.  Contemporary Demographic Reconstruction Methods Are Robust to Genome Assembly Quality: A Case Study in Tasmanian Devils , 2019, Molecular biology and evolution.

[10]  P. Maes,et al.  Complete Genome Sequence of a New Ebola Virus Strain Isolated during the 2017 Likati Outbreak in the Democratic Republic of the Congo , 2019, Microbiology Resource Announcements.

[11]  Menna E. Jones,et al.  Individual and temporal variation in pathogen load predicts long-term impacts of an emerging infectious disease , 2018, bioRxiv.

[12]  Menna E. Jones,et al.  Conserving adaptive potential: lessons from Tasmanian devils and their transmissible cancer , 2019, Conservation Genetics.

[13]  P. Ruiz‐Benito,et al.  Long‐term recovery of multifunctionality in Mediterranean forests depends on restoration strategy and forest type , 2019, Journal of Applied Ecology.

[14]  Austin H. Patton,et al.  The Genomic Basis of Tumor Regression in Tasmanian Devils (Sarcophilus harrisii) , 2018, Genome biology and evolution.

[15]  Menna E. Jones,et al.  Large‐effect loci affect survival in Tasmanian devils (Sarcophilus harrisii) infected with a transmissible cancer , 2018, Molecular ecology.

[16]  A. Eriksson,et al.  A metapopulation model of the spread of the Devil Facial Tumour Disease predicts the long term collapse of its host but not its extinction , 2018, bioRxiv.

[17]  Menna E. Jones,et al.  Density trends and demographic signals uncover the long-term impact of transmissible cancer in Tasmanian devils. , 2018, The Journal of applied ecology.

[18]  K. Bennett,et al.  The ERBB-STAT3 Axis Drives Tasmanian Devil Facial Tumor Disease , 2018, bioRxiv.

[19]  Daniel L. Ayres,et al.  Bayesian phylogenetic and phylodynamic data integration using BEAST 1.10 , 2018, Virus evolution.

[20]  Tanja Stadler,et al.  Detection of HIV transmission clusters from phylogenetic trees using a multi-state birth–death model , 2018, Journal of The Royal Society Interface.

[21]  Richard A Neher,et al.  TreeTime: Maximum-likelihood phylodynamic analysis , 2017, bioRxiv.

[22]  Menna E. Jones,et al.  Infection of the fittest: devil facial tumour disease has greatest effect on individuals with highest reproductive output. , 2017, Ecology letters.

[23]  M. Suchard,et al.  SpreaD3: Interactive Visualization of Spatiotemporal History and Trait Evolutionary Processes. , 2016, Molecular biology and evolution.

[24]  F. Cunningham,et al.  The Ensembl Variant Effect Predictor , 2016, Genome Biology.

[25]  Marek L Borowiec,et al.  AMAS: a fast tool for alignment manipulation and computing of summary statistics , 2016, PeerJ.

[26]  Andrew Rambaut,et al.  Exploring the temporal structure of heterochronous sequences using TempEst (formerly Path-O-Gen) , 2016, Virus evolution.

[27]  John J. Welch,et al.  The effect of genetic structure on molecular dating and tests for temporal signal , 2015, Methods in ecology and evolution.

[28]  A. B. Lyons,et al.  A second transmissible cancer in Tasmanian devils , 2015, Proceedings of the National Academy of Sciences.

[29]  J. Trent,et al.  Comparison against 186 canid whole-genome sequences reveals survival strategies of an ancient clonally transmissible canine tumor , 2015, Genome research.

[30]  Menna E. Jones,et al.  Transmissible cancer in Tasmanian devils: localized lineage replacement and host population response , 2015, Proceedings of the Royal Society B: Biological Sciences.

[31]  P. Ewald,et al.  Infection and cancer in multicellular organisms , 2015, Philosophical Transactions of the Royal Society B: Biological Sciences.

[32]  O. Pybus,et al.  Measurably evolving pathogens in the genomic era. , 2015, Trends in ecology & evolution.

[33]  A. von Haeseler,et al.  IQ-TREE: A Fast and Effective Stochastic Algorithm for Estimating Maximum-Likelihood Phylogenies , 2014, Molecular biology and evolution.

[34]  S. Vowler,et al.  N-acetyl-L-aspartyl-L-glutamate peptidase-like 2 is overexpressed in cancer and promotes a pro-migratory and pro-metastatic phenotype , 2014, Oncogene.

[35]  Dong Xie,et al.  BEAST 2: A Software Platform for Bayesian Evolutionary Analysis , 2014, PLoS Comput. Biol..

[36]  R Core Team,et al.  R: A language and environment for statistical computing. , 2014 .

[37]  Fabian Sievers,et al.  Clustal Omega, accurate alignment of very large numbers of sequences. , 2014, Methods in molecular biology.

[38]  杨凌春,et al.  Broad Institute , 2014 .

[39]  Guy Baele,et al.  Bayesian evolutionary model testing in the phylogenomics era: matching model complexity with computational efficiency , 2013, Bioinform..

[40]  Heng Li Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM , 2013, 1303.3997.

[41]  Trevor Bedford,et al.  Viral Phylodynamics , 2013, PLoS Comput. Biol..

[42]  K. Katoh,et al.  MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability , 2013, Molecular biology and evolution.

[43]  Menna E. Jones,et al.  Biting injuries and transmission of Tasmanian devil facial tumour disease. , 2013, The Journal of animal ecology.

[44]  S. Bonhoeffer,et al.  Birth–death skyline plot reveals temporal changes of epidemic spread in HIV and hepatitis C virus (HCV) , 2012, Proceedings of the National Academy of Sciences.

[45]  Carolyn Tregidgo,et al.  Genome Sequencing and Analysis of the Tasmanian Devil and Its Transmissible Cancer , 2012, Cell.

[46]  Steven Salzberg,et al.  BIOINFORMATICS ORIGINAL PAPER , 2004 .

[47]  Gonçalo R. Abecasis,et al.  The variant call format and VCFtools , 2011, Bioinform..

[48]  M. Suchard,et al.  Phylogeography takes a relaxed random walk in continuous space and time. , 2010, Molecular biology and evolution.

[49]  Christof Niehrs,et al.  Emerging links between CDK cell cycle regulators and Wnt signaling. , 2010, Trends in cell biology.

[50]  Cheryl J Briggs,et al.  Dynamics of an emerging disease drive large-scale amphibian population extinctions , 2010, Proceedings of the National Academy of Sciences.

[51]  Shelly Lachish,et al.  Transmission dynamics of Tasmanian devil facial tumor disease may lead to disease-induced extinction. , 2009, Ecology.

[52]  Gonçalo R. Abecasis,et al.  The Sequence Alignment/Map format and SAMtools , 2009, Bioinform..

[53]  M. Suchard,et al.  Smooth skyride through a rough skyline: Bayesian coalescent-based inference of population dynamics. , 2008, Molecular biology and evolution.

[54]  Clare E. Hawkins,et al.  Distribution and Impacts of Tasmanian Devil Facial Tumor Disease , 2007, EcoHealth.

[55]  K. Lafferty,et al.  Evidence for the Role of Infectious Disease in Species Extinction and Endangerment , 2006, Conservation biology : the journal of the Society for Conservation Biology.

[56]  Clare E. Hawkins,et al.  Emerging disease and population decline of an island endemic, the Tasmanian devil Sarcophilus harrisii , 2006 .

[57]  A. Pearse,et al.  Allograft theory: Transmission of devil facial-tumour disease , 2006, Nature.

[58]  O. Pybus,et al.  Bayesian coalescent inference of past population dynamics from molecular sequences. , 2005, Molecular biology and evolution.

[59]  O. Pybus,et al.  Unifying the Epidemiological and Evolutionary Dynamics of Pathogens , 2004, Science.

[60]  K. Strimmer,et al.  Exploring the demographic history of DNA sequences using the generalized skyline plot. , 2001, Molecular biology and evolution.

[61]  O. Pybus,et al.  An integrated framework for the inference of viral population history from reconstructed genealogies. , 2000, Genetics.

[62]  P. Daszak,et al.  Emerging infectious diseases of wildlife--threats to biodiversity and human health. , 2000, Science.

[63]  Peter Daszak,et al.  Extinction of a Species of Land Snail Due to Infection with a Microsporidian Parasite , 1998 .

[64]  A. Di Rienzo,et al.  Branching pattern in the evolutionary tree for human mitochondrial DNA. , 1991, Proceedings of the National Academy of Sciences of the United States of America.

[65]  S. Fields,et al.  Nucleotide sequence of human influenza A/PR/8/34 segment 2. , 1982, Nucleic acids research.