Evolution and lineage dynamics of a transmissible cancer in Tasmanian devils.

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.

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.

Density trends and demographic signals uncover the long-term impact of transmissible cancer in Tasmanian devils.

1. Monitoring the response of wild mammal populations to threatening processes is fundamental to effective conservation management. This is especially true for infectious diseases, which may have dynamic and therefore unpredictable interactions with their host. 2. We investigate the long-term impact of a transmissible cancer, devil facial tumour disease (DFTD), on the endemic Tasmanian devil. We analyse trends in devil spot-light counts and density across the area impacted by the disease. We investigate the demographic parameters which might be driving these trends, and use spatial capture-recapture models to examine whether DFTD has affected home range size. 3. We found that devils have declined by an average of 77% in areas affected by DFTD, and that there is a congruent trend of ongoing small decline in spotlight counts and density estimates. Despite this, devils have persisted to date within each of nine monitoring sites. One site is showing as yet unexplained small increases in density 8-10 years after the emergence of DFTD. 4. We also found the prevalence of DFTD has not abated despite large declines in density and that diseased sites continue to be dominated by young devils. The long-term impact of the disease has been partially offset by increased fecundity in the form of precocial breeding in 1-year-old females, and more pouch young per female in diseased sites. The lower densities resulting from DFTD did not affect home range size. 5. Synthesis and applications. Transmission of devil facial tumour disease continues despite large declines in devil density over multiple generations. Plasticity in life history traits has ameliorated the impact of devil facial tumour disease, however broad-scale trends in density show ongoing decline. In light of this, devil facial tumour disease and the impact of stochastic events on the reduced densities wrought by the disease, continue to threaten devils. In the absence of methods to manage disease in wild populations, we advocate managing the low population densities resulting from disease rather than disease per se.

Monitoring, imperfect detection, and risk optimization of a Tasmanian devil insurance population.

Most species are imperfectly detected during biological surveys, which creates uncertainty around their abundance or presence at a given location. Decision makers managing threatened or pest species are regularly faced with this uncertainty. Wildlife diseases can drive species to extinction; thus, managing species with disease is an important part of conservation. Devil facial tumor disease (DFTD) is one such disease that led to the listing of the Tasmanian devil (Sarcophilus harrisii) as endangered. Managers aim to maintain devils in the wild by establishing disease-free insurance populations at isolated sites. Often a resident DFTD-affected population must first be removed. In a successful collaboration between decision scientists and wildlife managers, we used an accessible population model to inform monitoring decisions and facilitate the establishment of an insurance population of devils on Forestier Peninsula. We used a Bayesian catch-effort model to estimate population size of a diseased population from removal and camera trap data. We also analyzed the costs and benefits of declaring the area disease-free prior to reintroduction and establishment of a healthy insurance population. After the monitoring session in May-June 2015, the probability that all devils had been successfully removed was close to 1, even when we accounted for a possible introduction of a devil to the site. Given this high probability and the baseline cost of declaring population absence prematurely, we found it was not cost-effective to carry out any additional monitoring before introducing the insurance population. Considering these results within the broader context of Tasmanian devil management, managers ultimately decided to implement an additional monitoring session before the introduction. This was a conservative decision that accounted for uncertainty in model estimates and for the broader nonmonetary costs of mistakenly declaring the area disease-free.

Analysis of extended HLA haplotypes in multiple sclerosis and narcolepsy families confirms a predisposing effect for the class I region in Tasmanian MS patients.

Human leucocyte antigen (HLA)-DRB1*15 is associated with predisposition to multiple sclerosis (MS), although conjecture surrounds the possible involvement of an alternate risk locus in the class I region of the HLA complex. We have shown previously that an alternate MS risk allele(s) may be encompassed by the telomerically extended DRB1*15 haplotype, and here, we have attempted to map the putative variant. Thirteen microsatellite markers encompassing a 6.79-megabase (D6S2236-G51152) region, and the DRB1 and DQB1 genes, were genotyped in 166 MS simplex families and 104 control families from the Australian State of Tasmania and 153 narcolepsy simplex families (trios) from the USA. Complementary approaches were used to investigate residual predisposing effects of microsatellite alleles comprising the extended DRB1*15 haplotype taking into account the strong predisposing effect of DRB1*15: (1) Disease association of the extended DRB1*15 haplotype was compared for MS and narcolepsy families--predisposing effects were observed for extended class I microsatellite marker alleles in MS families, but not narcolepsy families; (2) disease association of the extended DRB1*15 haplotype was investigated after conditioning MS and control haplotypes on the absence of DRB1*15--a significant predisposing effect was observed for a 627-kb haplotype (D6S258 allele 8-MOGCA allele 4; MOG, myelin oligodendrocyte glycoprotein) spanning the extended class I region. MOGCA allele 4 displayed the strongest predisposing effect in DRB1*15-conditioned haplotypes (p = 0.0006; OR 2.83 [1.54-5.19]). Together, these data confirm that an alternate MS risk locus exists in the extended class I region in Tasmanian MS patients independent of DRB1*15.

On the utility of data from the International HapMap Project for Australian association studies.

We compare patterns of linkage disequilibrium (LD) for 633 SNPs in two regions between samples collected in two Australian states and HapMap samples collected from Utah residents of Northern and Western (NW) European ancestry (CEU). Patterns of LD in the Australian and HapMap samples are similar, and tag SNPs chosen using HapMap genotypes perform almost as well on Australian samples as tags chosen using Australian genotypes.