The evolution of two transmissible cancers in Tasmanian devils

Extant and extinct bilby genomes combined with Indigenous knowledge improve conservation of a unique Australian marsupial

Ninu (greater bilby, Macrotis lagotis) are desert-dwelling, culturally and ecologically important marsupials. In collaboration with Indigenous rangers and conservation managers, we generated the Ninu chromosome-level genome assembly (3.66 Gbp) and genome sequences for the extinct Yallara (lesser bilby, Macrotis leucura). We developed and tested a scat single-nucleotide polymorphism panel to inform current and future conservation actions, undertake ecological assessments and improve our understanding of Ninu genetic diversity in managed and wild populations. We also assessed the beneficial impact of translocations in the metapopulation (N = 363 Ninu). Resequenced genomes (temperate Ninu, 6; semi-arid Ninu, 6; and Yallara, 4) revealed two major population crashes during global cooling events for both species and differences in Ninu genes involved in anatomical and metabolic pathways. Despite their 45-year captive history, Ninu have fewer long runs of homozygosity than other larger mammals, which may be attributable to their boom-bust life history. Here we investigated the unique Ninu biology using 12 tissue transcriptomes revealing expression of all 115 conserved eutherian chorioallantoic placentation genes in the uterus, an XY1Y2 sex chromosome system and olfactory receptor gene expansions. Together, we demonstrate the holistic value of genomics in improving key conservation actions, understanding unique biological traits and developing tools for Indigenous rangers to monitor remote wild populations.

Transmissible cancers, the genomes that do not melt down

Evolutionary theory predicts that the accumulation of deleterious mutations in asexually reproducing organisms should lead to genomic decay. Clonally reproducing cell lines, i.e., transmissible cancers, when cells are transmitted as allografts/xenografts, break these rules and survive for centuries and millennia. The currently known 11 transmissible cancer lineages occur in dogs (canine venereal tumour disease), in Tasmanian devils (devil facial tumor diseases, DFT1 and DFT2), and in bivalves (bivalve transmissible neoplasia). Despite the mutation loads of these cell lines being much higher than observed in human cancers, they have not been eliminated in space and time. Here, we provide potential explanations for how these fascinating cell lines may have overcome the fitness decline due to the progressive accumulation of deleterious mutations and propose that the high mutation load may carry an indirect positive fitness outcome. We offer ideas on how these host-pathogen systems could be used to answer outstanding questions in evolutionary biology. The recent studies on the evolution of these clonal pathogens reveal key mechanistic insight into transmissible cancer genomes, information that is essential for future studies investigating how these contagious cancer cell lines can repeatedly evade immune recognition, evolve, and survive in the landscape of highly diverse hosts.

The Evolution of Transglutaminases Underlies the Origin and Loss of Cornified Skin Appendages in Vertebrates

Transglutaminases (TGMs) cross-link proteins by introducing covalent bonds between glutamine and lysine residues. These cross-links are essential for epithelial cornification which enables tetrapods to live on land. Here, we investigated which evolutionary adaptations of vertebrates were associated with specific changes in the family of TGM genes. We determined the catalog of TGMs in the main clades of vertebrates, performed a comprehensive phylogenetic analysis of TGMs, and localized the distribution of selected TGMs in tissues. Our data suggest that TGM1 is the phylogenetically oldest epithelial TGM, with orthologs being expressed in the cornified teeth of the lamprey, a basal vertebrate. Gene duplications led to the origin of TGM10 in stem vertebrates, the origin of TGM2 in jawed vertebrates, and an increasing number of epithelium-associated TGM genes in the lineage leading to terrestrial vertebrates. TGM9 is expressed in the epithelial egg tooth, and its evolutionary origin in stem amniotes coincided with the evolution of embryonic development in eggs that are surrounded by a protective shell. Conversely, viviparous mammals have lost both the epithelial egg tooth and TGM9. TGM3 and TGM6 evolved as regulators of cornification in hair follicles and underwent pseudogenization upon the evolutionary loss of hair in cetaceans. Taken together, this study reveals the gain and loss of vertebrate TGM genes in association with the evolution of cornified skin appendages and suggests an important role of TGM9 in the evolution of amniotes.

Genome of the endangered eastern quoll (Dasyurus viverrinus) reveals signatures of historical decline and pelage color evolution

The eastern quoll (Dasyurus viverrinus) is an endangered marsupial native to Australia. Since the extirpation of its mainland populations in the 20th century, wild eastern quolls have been restricted to two islands at the southern end of their historical range. Eastern quolls are the subject of captive breeding programs and attempts have been made to re-establish a population in mainland Australia. However, few resources currently exist to guide the genetic management of this species. Here, we generated a reference genome for the eastern quoll with gene annotations supported by multi-tissue transcriptomes. Our assembly is among the most complete marsupial genomes currently available. Using this assembly, we infer the species' demographic history, identifying potential evidence of a long-term decline beginning in the late Pleistocene. Finally, we identify a deletion at the ASIP locus that likely underpins pelage color differences between the eastern quoll and the closely related Tasmanian devil (Sarcophilus harrisii).

Translating genomic advances into biodiversity conservation

A key action of the new Global Biodiversity Framework is the maintenance of genetic diversity in all species to safeguard their adaptive potential. To achieve this goal, a translational mindset, which aims to convert results of basic research into direct practical benefits, needs to be applied to biodiversity conservation. Despite much discussion on the value of genomics to conservation, a disconnect between those generating genomic resources and those applying it to biodiversity management remains. As global efforts to generate reference genomes for non-model species increase, investment into practical biodiversity applications is critically important. Applications such as understanding population and multispecies diversity and longitudinal monitoring need support alongside education for policymakers on integrating the data into evidence-based decisions. Without such investment, the opportunity to revolutionize global biodiversity conservation using genomics will not be fully realized.

Medicinal Chemistry Gone Wild

Earth is currently experiencing a mass extinction event. The flora and fauna of our planet are experiencing mass die-offs from a multitude of factors, with wildlife disease emerging as one parameter where medicinal chemists are equipped to intervene. While contemporary medicinal chemistry focuses on human health, many traditional pharmaceutical companies have historic roots in human health, animal health, and plant health. This trifecta of health sciences perfectly maps to the current field of One Health, which recognizes that optimal health outcomes can only be achieved through the health of humans, plants, animals, and their shared environments. This Perspective imagines a world where state-of-the-art medicinal chemistry tactics are used to prevent the extinction of endangered species and points to preliminary work in the emerging area of conservation medicine.

No evidence that a transmissible cancer has shifted from emergence to endemism in Tasmanian devils

Tasmanian devils are endangered by a transmissible cancer known as Tasmanian devil facial tumour 1 (DFT1). A 2020 study by Patton et al. (Science 370, eabb9772 (doi:10.1126/science.abb9772)) used genome data from DFT1 tumours to produce a dated phylogenetic tree for this transmissible cancer lineage, and thence, using phylodynamics models, to estimate its epidemiological parameters and predict its future trajectory. It concluded that the effective reproduction number for DFT1 had declined to a value of one, and that the disease had shifted from emergence to endemism. We show that the study is based on erroneous mutation calls and flawed methodology, and that its conclusions cannot be substantiated.

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.

In vitro competition between two transmissible cancers and potential implications for their host, the Tasmanian devil

Since the emergence of a transmissible cancer, devil facial tumour disease (DFT1), in the 1980s, wild Tasmanian devil populations have been in decline. In 2016, a second, independently evolved transmissible cancer (DFT2) was discovered raising concerns for survival of the host species. Here, we applied experimental and modelling frameworks to examine competition dynamics between the two transmissible cancers in vitro. Using representative cell lines for DFT1 and DFT2, we have found that in monoculture, DFT2 grows twice as fast as DFT1 but reaches lower maximum cell densities. Using co-cultures, we demonstrate that DFT2 outcompetes DFT1: the number of DFT1 cells decreasing over time, never reaching exponential growth. This phenomenon could not be replicated when cells were grown separated by a semi-permeable membrane, consistent with exertion of mechanical stress on DFT1 cells by DFT2. A logistic model and a Lotka-Volterra competition model were used to interrogate monoculture and co-culture growth curves, respectively, suggesting DFT2 is a better competitor than DFT1, but also showing that competition outcomes might depend on the initial number of cells, at least in the laboratory. We provide theories how the in vitro results could be translated to observations in the wild and propose that these results may indicate that although DFT2 is currently in a smaller geographic area than DFT1, it could have the potential to outcompete DFT1. Furthermore, we provide a framework for improving the parameterization of epidemiological models applied to these cancer lineages, which will inform future disease management.

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.

Comparative genomics of monotremes provides insights into the early evolution of mammalian epidermal differentiation genes

The function of the skin as a barrier against the environment depends on the differentiation of epidermal keratinocytes into highly resilient corneocytes that form the outermost skin layer. Many genes encoding structural components of corneocytes are clustered in the epidermal differentiation complex (EDC), which has been described in placental and marsupial mammals as well as non-mammalian tetrapods. Here, we analyzed the genomes of the platypus (Ornithorhynchus anatinus) and the echidna (Tachyglossus aculeatus) to determine the gene composition of the EDC in the basal clade of mammals, the monotremes. We report that mammal-specific subfamilies of EDC genes encoding small proline-rich proteins (SPRRs) and late cornified envelope proteins as well as single-copy EDC genes such as involucrin are conserved in monotremes, suggesting that they have originated in stem mammals. Monotremes have at least one gene homologous to the group of filaggrin (FLG), FLG2 and hornerin (HRNR) in placental mammals, but no clear one-to-one pairwise ortholog of either FLG, FLG2 or HRNR. Caspase-14, a keratinocyte differentiation-associated protease implicated in the processing of filaggrin, is encoded by at least 3 gene copies in the echidna. Our results reveal evolutionarily conserved and clade-specific features of the genetic regulation of epidermal differentiation in monotremes.

Whole-genome doubling in tissues and tumors

The overwhelming majority of proliferating somatic human cells are diploid, and this genomic state is typically maintained across successive cell divisions. However, failures in cell division can induce a whole-genome doubling (WGD) event, in which diploid cells transition to a tetraploid state. While some WGDs are developmentally programmed to produce nonproliferative tetraploid cells with specific cellular functions, unscheduled WGDs can be catastrophic: erroneously arising tetraploid cells are ill-equipped to cope with their doubled cellular and chromosomal content and quickly become genomically unstable and tumorigenic. Deciphering the genetics that underlie the genesis, physiology, and evolution of whole-genome doubled (WGD+) cells may therefore reveal therapeutic avenues to selectively eliminate pathological WGD+ cells.

Comparative genomics of the T cell receptor μ locus in marsupials and monotremes

T cells are a primary component of the vertebrate adaptive immune system. There are three mammalian T cell lineages based on their T cell receptors (TCR). The αβ T cells and γδ T cells are ancient and found broadly in vertebrates. The more recently discovered γμ T cells are uniquely mammalian and only found in marsupials and monotremes. In this study, we compare the TCRμ locus (TRM) across the genomes of two marsupials, the gray short-tailed opossum and Tasmanian devil, and one monotreme, the platypus. These analyses revealed lineage-specific duplications, common to all non-eutherian mammals described. There is conserved synteny in the TRM loci of both marsupials but not in the monotreme. Our results are consistent with an ancestral cluster organization which was present in the last common mammalian ancestor which underwent lineage-specific duplications and divergence among the non-eutherian mammals.

Centuries of genome instability and evolution in soft-shell clam, Mya arenaria, bivalve transmissible neoplasia

Transmissible cancers are infectious parasitic clones that metastasize to new hosts, living past the death of the founder animal in which the cancer initiated. We investigated the evolutionary history of a cancer lineage that has spread though the soft-shell clam (Mya arenaria) population by assembling a chromosome-scale soft-shell clam reference genome and characterizing somatic mutations in transmissible cancer. We observe high mutation density, widespread copy-number gain, structural rearrangement, loss of heterozygosity, variable telomere lengths, mitochondrial genome expansion and transposable element activity, all indicative of an unstable cancer genome. We also discover a previously unreported mutational signature associated with overexpression of an error-prone polymerase and use this to estimate the lineage to be >200 years old. Our study reveals the ability for an invertebrate cancer lineage to survive for centuries while its genome continues to structurally mutate, likely contributing to the evolution of this lineage as a parasitic cancer.

Agential autonomy and biological individuality

What is a biological individual? How are biological individuals individuated? How can we tell how many individuals there are in a given assemblage of biological entities? The individuation and differentiation of biological individuals are central to the scientific understanding of living beings. I propose a novel criterion of biological individuality according to which biological individuals are autonomous agents. First, I articulate an ecological-dynamical account of natural agency according to which, agency is the gross dynamical capacity of a goal-directed system to bias its repertoire to respond to its conditions as affordances. Then, I argue that agents or agential dynamical systems can be agentially dependent on, or agentially autonomous from, other agents and that this agential dependence/autonomy can be symmetrical or asymmetrical, strong or weak. Biological individuals, I propose, are all and only those agential dynamical systems that are strongly agentially autonomous. So, to determine how many individuals there are in a given multiagent aggregate, such as multicellular organism, a colony, symbiosis, or a swarm, we first have to identify how many agential dynamical systems there are, and then what their relations of agential dependence/autonomy are. I argue that this criterion is adequate to the extent that it vindicates the paradigmatic cases, and explains why the paradigmatic cases are paradigmatic, and why the problematic cases are problematic. Finally, I argue for the importance of distinguishing between agential and causal dependence and show the relevance of agential autonomy for understanding the explanatory structure of evolutionary developmental biology.

Somatic evolution of marine transmissible leukemias in the common cockle, Cerastoderma edule

Transmissible cancers are malignant cell lineages that spread clonally between individuals. Several such cancers, termed bivalve transmissible neoplasia (BTN), induce leukemia-like disease in marine bivalves. This is the case of BTN lineages affecting the common cockle, Cerastoderma edule, which inhabits the Atlantic coasts of Europe and northwest Africa. To investigate the evolution of cockle BTN, we collected 6,854 cockles, diagnosed 390 BTN tumors, generated a reference genome and assessed genomic variation across 61 tumors. Our analyses confirmed the existence of two BTN lineages with hemocytic origins. Mitochondrial variation revealed mitochondrial capture and host co-infection events. Mutational analyses identified lineage-specific signatures, one of which likely reflects DNA alkylation. Cytogenetic and copy number analyses uncovered pervasive genomic instability, with whole-genome duplication, oncogene amplification and alkylation-repair suppression as likely drivers. Satellite DNA distributions suggested ancient clonal origins. Our study illuminates long-term cancer evolution under the sea and reveals tolerance of extreme instability in neoplastic genomes.

Historical RNA expression profiles from the extinct Tasmanian tiger

Paleogenomics continues to yield valuable insights into the evolution, population dynamics, and ecology of our ancestors and other extinct species. However, DNA sequencing cannot reveal tissue-specific gene expression, cellular identity, or gene regulation, which are only attainable at the transcriptional level. Pioneering studies have shown that useful RNA can be extracted from ancient specimens preserved in permafrost and historical skins from extant canids, but no attempts have been made so far on extinct species. We extract, sequence, and analyze historical RNA from muscle and skin tissue of a ∼130-year-old Tasmanian tiger (Thylacinus cynocephalus) preserved in desiccation at room temperature in a museum collection. The transcriptional profiles closely resemble those of extant species, revealing specific anatomical features such as slow muscle fibers or blood infiltration. Metatranscriptomic analysis, RNA damage, tissue-specific RNA profiles, and expression hotspots genome-wide further confirm the thylacine origin of the sequences. RNA sequences are used to improve protein-coding and noncoding annotations, evidencing missing exonic loci and the location of ribosomal RNA genes while increasing the number of annotated thylacine microRNAs from 62 to 325. We discover a thylacine-specific microRNA isoform that could not have been confirmed without RNA evidence. Finally, we detect traces of RNA viruses, suggesting the possibility of profiling viral evolution. Our results represent the first successful attempt to obtain transcriptional profiles from an extinct animal species, providing thought-to-be-lost information on gene expression dynamics. These findings hold promising implications for the study of RNA molecules across the vast collections of natural history museums and from well-preserved permafrost remains.

Cancerogenic parasites in veterinary medicine: a narrative literature review

Parasite infection is one of the many environmental factors that can significantly contribute to carcinogenesis and is already known to be associated with a variety of malignancies in both human and veterinary medicine. However, the actual number of cancerogenic parasites and their relationship to tumor development is far from being fully understood, especially in veterinary medicine. Thus, the aim of this review is to investigate parasite-related cancers in domestic and wild animals and their burden in veterinary oncology. Spontaneous neoplasia with ascertained or putative parasite etiology in domestic and wild animals will be reviewed, and the multifarious mechanisms of protozoan and metazoan cancer induction will be discussed.

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.