Transmissible cancers are clonal lineages that spread through populations via contagious cancer cells. Although rare in nature, two facial tumor clones affect Tasmanian devils. Here we perform comparative genetic and functional characterization of these lineages. The two cancers have similar patterns of mutation and show no evidence of exposure to exogenous mutagens or viruses. Genes encoding PDGF receptors have copy number gains and are present on extrachromosomal double minutes. Drug screening indicates causative roles for receptor tyrosine kinases and sensitivity to inhibitors of DNA repair. Y chromosome loss from a male clone infecting a female host suggests immunoediting. These results imply that Tasmanian devils may have inherent susceptibility to transmissible cancers and present a suite of therapeutic compounds for use in conservation.

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.

Streptococcus agalactiae or Group B Streptococcus (GBS) are opportunistic bacteria that can cause lethal sepsis in children and immuno-compromised patients. Their genome is a reservoir of mobile genetic elements that can be horizontally transferred. Among them, integrative and conjugative elements (ICEs) and the smaller integrative and mobilizable elements (IMEs) primarily reside in the bacterial chromosome, yet have the ability to be transferred between cells by conjugation. ICEs and IMEs are therefore a source of genetic variability that participates in the spread of antibiotic resistance. Although IMEs seem to be the most prevalent class of elements transferable by conjugation, they are poorly known. Here, we have studied a GBS-IME, termed IMESag-rpsI, which is widely distributed in GBS despite not carrying any apparent virulence trait. Analyses of 240 whole genomes showed that IMESag-rpsI is present in approximately 47% of the genomes, has a roughly constant size (approx. 9 kb) and is always integrated at a single location, the 3'-end of the gene encoding the ribosomal protein S9 (rpsI). Based on their genetic variation, several IMESag-rpsI types were defined (A-J) and classified in clonal complexes (CCs). CC1 was the most populated by IMESag-rpsI (more than 95%), mostly of type-A (71%). One CC1 strain (S. agalactiae HRC) was deep-sequenced to understand the rationale underlying type-A IMESag-rpsI enrichment in GBS. Thirteen open reading frames were identified, one of them encoding a protein (MobSag) belonging to the broadly distributed family of relaxases MOBV1 Protein MobSag was purified and, by a newly developed method, shown to cleave DNA at a specific dinucleotide. The S. agalactiae HRC-IMESag-rpsI is able to excise from the chromosome, as shown by the presence of circular intermediates, and it harbours a fully functional mobilization module. Further, the mobSag gene encoded by this mobile element is able to promote plasmid transfer among pneumococcal strains, suggesting that MobSag facilitates the spread of IMESag-rpsI and that this spread would explain the presence of the same IMESag-rpsI type in GBS strains belonging to different CCs.

The rates and patterns of somatic mutation in normal tissues are largely unknown outside of humans1-7. Comparative analyses can shed light on the diversity of mutagenesis across species, and on long-standing hypotheses about the evolution of somatic mutation rates and their role in cancer and ageing. Here we performed whole-genome sequencing of 208 intestinal crypts from 56 individuals to study the landscape of somatic mutation across 16 mammalian species. We found that somatic mutagenesis was dominated by seemingly endogenous mutational processes in all species, including 5-methylcytosine deamination and oxidative damage. With some differences, mutational signatures in other species resembled those described in humans8, although the relative contribution of each signature varied across species. Notably, the somatic mutation rate per year varied greatly across species and exhibited a strong inverse relationship with species lifespan, with no other life-history trait studied showing a comparable association. Despite widely different life histories among the species we examined-including variation of around 30-fold in lifespan and around 40,000-fold in body mass-the somatic mutation burden at the end of lifespan varied only by a factor of around 3. These data unveil common mutational processes across mammals, and suggest that somatic mutation rates are evolutionarily constrained and may be a contributing factor in ageing.

About half of all cancers have somatic integrations of retrotransposons. Here, to characterize their role in oncogenesis, we analyzed the patterns and mechanisms of somatic retrotransposition in 2,954 cancer genomes from 38 histological cancer subtypes within the framework of the Pan-Cancer Analysis of Whole Genomes (PCAWG) project. We identified 19,166 somatically acquired retrotransposition events, which affected 35% of samples and spanned a range of event types. Long interspersed nuclear element (LINE-1; L1 hereafter) insertions emerged as the first most frequent type of somatic structural variation in esophageal adenocarcinoma, and the second most frequent in head-and-neck and colorectal cancers. Aberrant L1 integrations can delete megabase-scale regions of a chromosome, which sometimes leads to the removal of tumor-suppressor genes, and can induce complex translocations and large-scale duplications. Somatic retrotranspositions can also initiate breakage-fusion-bridge cycles, leading to high-level amplification of oncogenes. These observations illuminate a relevant role of L1 retrotransposition in remodeling the cancer genome, with potential implications for the development of human tumors.

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.

Autonomous replication and segregation of mitochondrial DNA (mtDNA) creates the potential for evolutionary conflict driven by emergence of haplotypes under positive selection for 'selfish' traits, such as replicative advantage. However, few cases of this phenomenon arising within natural populations have been described. Here, we survey the frequency of mtDNA horizontal transfer within the canine transmissible venereal tumour (CTVT), a contagious cancer clone that occasionally acquires mtDNA from its hosts. Remarkably, one canine mtDNA haplotype, A1d1a, has repeatedly and recently colonised CTVT cells, recurrently replacing incumbent CTVT haplotypes. An A1d1a control region polymorphism predicted to influence transcription is fixed in the products of an A1d1a recombination event and occurs somatically on other CTVT mtDNA backgrounds. We present a model whereby 'selfish' positive selection acting on a regulatory variant drives repeated fixation of A1d1a within CTVT cells.

Most cancers are characterized by the somatic acquisition of genomic rearrangements during tumour evolution that eventually drive the oncogenesis. Here, using multiplatform sequencing technologies, we identify and characterize a remarkable mutational mechanism in human hepatocellular carcinoma caused by Hepatitis B virus, by which DNA molecules from the virus are inserted into the tumour genome causing dramatic changes in its configuration, including non-homologous chromosomal fusions, dicentric chromosomes and megabase-size telomeric deletions. This aberrant mutational mechanism, present in at least 8% of all HCC tumours, can provide the driver rearrangements that a cancer clone requires to survive and grow, including loss of relevant tumour suppressor genes. Most of these events are clonal and occur early during liver cancer evolution. Real-time timing estimation reveals some HBV-mediated rearrangements occur as early as two decades before cancer diagnosis. Overall, these data underscore the importance of characterising liver cancer genomes for patterns of HBV integration.