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.

Human habitat modification, not apex scavenger decline, drives isotopic niche variation in a carnivore community.

Top carnivores can influence the structure of ecological communities, primarily through competition and predation; however, communities are also influenced by bottom-up forces such as anthropogenic habitat disturbance. Top carnivore declines will likely alter competitive dynamics within and amongst sympatric carnivore species. Increasing intraspecific competition is generally predicted to drive niche expansion and/or individual specialisation, while interspecific competition tends to constrain niches. Using stable isotope analysis of whiskers, we studied the effects of Tasmanian devil Sarcophilus harrisii declines upon the population- and individual-level isotopic niches of Tasmanian devils and sympatric spotted-tailed quolls Dasyurus maculatus subsp. maculatus. We investigated whether time since the onset of devil decline (a proxy for severity of decline) and landscape characteristics affected the isotopic niche breadth and overlap of devil and quoll populations. We quantified individual isotopic niche breadth for a subset of Tasmanian devils and spotted-tailed quolls and assessed whether between-site population niche variation was driven by individual-level specialisation. Tasmanian devils and spotted-tailed quolls demonstrated smaller population-level isotopic niche breadths with increasing human-modified habitat, while time since the onset of devil decline had no effect on population-level niche breadth or interspecific niche overlap. Individual isotopic niche breadths of Tasmanian devils and spotted-tailed quolls were narrower in human-modified landscapes, likely driving population isotopic niche contraction, however, the degree of individuals' specialisation relative to one another remained constant. Our results suggest that across varied landscapes, mammalian carnivore niches can be more sensitive to the bottom-up forces of anthropogenic habitat disturbance than to the top-down effects of top carnivore decline.

Over-eruption in marsupial carnivore teeth: compensation for a constraint.

Pronounced over-eruption of the canine teeth, causing the cervical enamel margin to extend beyond the alveolar bone and exposing the root, occurs with age and growth in Australian marsupial carnivores, much more than in eco-morphologically equivalent placental carnivores. Suppression of functional tooth replacement is characteristic of marsupials, where most placentals have the primitive diphyodont pattern of two generations of incisor, canine and premolar teeth. Canine and molar tooth dimensions of four species of marsupial carnivores (thylacine Thylacinus cynocephalus, Tasmanian devil Sarcophilus harrisii and two quolls Dasyurus spp.) and canine dimensions of seven eco-morphologically equivalent placental carnivore species were measured from museum specimens. Canine dimensions were measured in a time series on live wild-living individual devils and quolls. The canine teeth and to a lesser extent the molar teeth of marsupial carnivores continue to erupt through life, resulting in a net increase in tooth height and diameter, a phenomenon not evident in placental carnivores. Potential mechanisms causing over-eruption include tooth wear and gradual release of occlusal pressure as the individual grows. Over-eruption in marsupial carnivores may be a compensatory response for tooth size limits imposed by monophyodont tooth replacement, ensuring that animal's teeth are scaled to jaw size from juvenile to adulthood.

Top predator restricts the niche breadth of prey: effects of assisted colonization of Tasmanian devils on a widespread omnivorous prey.

Few landscape-scale experiments test the effects of predators on the abundance and distribution of prey across habitat gradients. We use the assisted colonization of a top predator, the Tasmanian devil (Sarcophilus harrisii), to test the impacts of predation on the abundance, habitat use and temporal activity of a widespread prey species, the omnivorous common brushtail possum (Trichosurus vulpecula). Before introduction of devils to Maria Island, Tasmania, Australia, in 2012, possums were abundant in open grasslands as well as forests. Predation by devils caused high mortality of possums in grasslands, but individuals with access to trees had a higher survival probability. Possum abundance declined across the whole island from 2012-2016, as possums disappeared almost completely from grasslands and declined in drier forests with more open understorey. Abundance remained stable in wet forests, which are not preferred habitat for possums but provide better refuge from devils. Abundance and habitat use of possums remained unchanged at a control site on the adjacent Tasmanian mainland, where the devil population was low and stable. This study demonstrates how spatial variation in predator-caused mortality can limit both abundance and habitat breadth in generalist prey species, excluding them entirely from certain habitats.

Dominant carnivore loss benefits native avian and invasive mammalian scavengers.

Scavenging by large carnivores is integral for ecosystem functioning by limiting the build-up of carrion and facilitating widespread energy flows. However, top carnivores have declined across the world, triggering trophic shifts within ecosystems. Here, we compare findings from previous work on predator decline against areas with recent native mammalian carnivore loss. Specifically, we investigate top-down control on utilization of experimentally placed carcasses by two mesoscavengers-the invasive feral cat and native forest raven. Ravens profited most from carnivore loss, scavenging for five times longer in the absence of native mammalian carnivores. Cats scavenged on half of all carcasses in the region without dominant native carnivores. This was eight times more than in areas where other carnivores were at high densities. All carcasses persisted longer than the three-week monitoring period in the absence of native mammalian carnivores, while in areas with high carnivore abundance, all carcasses were fully consumed. Our results reveal that top-carnivore loss amplifies impacts associated with carnivore decline-increased carcass persistence and carrion access for smaller scavengers. This suggests that even at low densities, native mammalian carnivores can fulfil their ecological functions, demonstrating the significance of global carnivore conservation and supporting management approaches, such as trophic rewilding.

Extracellular vesicle proteomes of two transmissible cancers of Tasmanian devils reveal tenascin-C as a serum-based differential diagnostic biomarker.

The iconic Tasmanian devil (Sarcophilus harrisii) is endangered due to the transmissible cancer Devil Facial Tumour Disease (DFTD), of which there are two genetically independent subtypes (DFT1 and DFT2). While DFT1 and DFT2 can be differentially diagnosed using tumour biopsies, there is an urgent need to develop less-invasive biomarkers that can detect DFTD and distinguish between subtypes. Extracellular vesicles (EVs), the nano-sized membrane-enclosed vesicles present in most biofluids, represent a valuable resource for biomarker discovery. Here, we characterized the proteome of EVs from cultured DFTD cells using data-independent acquisition-mass spectrometry and an in-house spectral library of > 1500 proteins. EVs from both DFT1 and DFT2 cell lines expressed higher levels of proteins associated with focal adhesion functions. Furthermore, hallmark proteins of epithelial-mesenchymal transition were enriched in DFT2 EVs relative to DFT1 EVs. These findings were validated in EVs derived from serum samples, revealing that the mesenchymal marker tenascin-C was also enriched in EVs derived from the serum of devils infected with DFT2 relative to those infected with DFT1 and healthy controls. This first EV-based investigation of DFTD increases our understanding of the cancers' EVs and their possible involvement in DFTD progression, such as metastasis. Finally, we demonstrated the potential of EVs to differentiate between DFT1 and DFT2, highlighting their potential use as less-invasive liquid biopsies for the Tasmanian devil.

Spatial variation in gene expression of Tasmanian devil facial tumors despite minimal host transcriptomic response to infection.

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.

A triple threat: high population density, high foraging intensity and flexible habitat preferences explain high impact of feral cats on prey.

Alien mammalian carnivores have contributed disproportionately to global loss of biodiversity. In Australia, predation by the feral cat and red fox is one of the most significant causes of the decline of native vertebrates. To discover why cats have greater impacts on prey than native predators, we compared the ecology of the feral cat to a marsupial counterpart, the spotted-tailed quoll. Individual prey are 20-200 times more likely to encounter feral cats, because of the combined effects of cats' higher population densities, greater intensity of home-range use and broader habitat preferences. These characteristics also mean that the costs to the prey of adopting anti-predator behaviours against feral cats are likely to be much higher than adopting such behaviours in response to spotted-tailed quolls, due to the reliability and ubiquity of feral cat cues. These results help explain the devastating impacts of cats on wildlife in Australia and other parts of the world.

Contemporary and historical selection in Tasmanian devils (Sarcophilus harrisii) support novel, polygenic response to transmissible cancer.

Tasmanian devils (Sarcophilus harrisii) are evolving in response to a unique transmissible cancer, devil facial tumour disease (DFTD), first described in 1996. Persistence of wild populations and the recent emergence of a second independently evolved transmissible cancer suggest that transmissible cancers may be a recurrent feature in devils. Here, we compared signatures of selection across temporal scales to determine whether genes or gene pathways under contemporary selection (six to eight generations) have also been subject to historical selection (65-85 Myr). First, we used targeted sequencing, RAD-capture, in approximately 2500 devils in six populations to identify genomic regions subject to rapid evolution. We documented genome-wide contemporary evolution, including 186 candidate genes related to cell cycling and immune response. Then we used a molecular evolution approach to identify historical positive selection in devils compared to other marsupials and found evidence of selection in 1773 genes. However, we found limited overlap across time scales, with only 16 shared candidate genes, and no overlap in enriched functional gene sets. Our results are consistent with a novel, multi-locus evolutionary response of devils to DFTD. Our results can inform conservation by identifying high priority targets for genetic monitoring and guiding maintenance of adaptive potential in managed populations.

The effects of weather variability on patterns of genetic diversity in Tasmanian bettongs.

While the effects of climate (long-term, prevailing weather) on species abundance, range and genetic diversity have been widely studied, short-term, localized variations in atmospheric conditions (i.e., weather) can also rapidly alter species' geographical ranges and population sizes, but little is known about how they affect genetic diversity. We investigated the relationship between weather and range-wide genetic diversity in a marsupial, Bettongia gaimardi, using dynamic species distribution models (SDMs). Genetic diversity was lower in parts of the range where the weather-based SDM predicted high variability in probability of B. gaimardi occurrence during 1950-2009. This is probably an effect of lower population sizes and extinction-recolonization cycles in places with highly variable weather. Spatial variation in genetic diversity was also better predicted by mean probabilities of B. gaimardi occurrence from weather- than climate-based SDMs. Our results illustrate the importance of weather in driving population dynamics and species distributions on decadal timescales and thereby in affecting genetic diversity. Modelling the links between changing weather patterns, species distributions and genetic diversity will allow researchers to better forecast biological impacts of climate change.

A native apex predator limits an invasive mesopredator and protects native prey: Tasmanian devils protecting bandicoots from cats.

Apex predators can limit the abundance and behaviour of mesopredators, thereby reducing predation on smaller species. We know less about whether native apex predators are effective in suppressing invasive mesopredators, a major global driver of vertebrate extinctions. We use the severe disease-induced decline of an apex predator, the Tasmanian devil, as a natural experiment to test whether devils limit abundance of invasive feral cats and in turn protect smaller native prey. Cat abundance was c. 58% higher where devils had declined, which in turn negatively affected a smaller native prey species. Devils had a stronger limiting effect on cats than on a native mesopredator, suggesting apex predators may have stronger suppressive effects on evolutionarily naive species than coevolved species. Our results highlight how disease in one species can affect the broader ecosystem. We show that apex predators not only regulate native species but can also confer resistance to the impacts of invasive populations. Apex predators could therefore be a powerful but underutilised tool to prevent biodiversity loss.

Contemporary Demographic Reconstruction Methods Are Robust to Genome Assembly Quality: A Case Study in Tasmanian Devils.

Reconstructing species' demographic histories is a central focus of molecular ecology and evolution. Recently, an expanding suite of methods leveraging either the sequentially Markovian coalescent (SMC) or the site-frequency spectrum has been developed to reconstruct population size histories from genomic sequence data. However, few studies have investigated the robustness of these methods to genome assemblies of varying quality. In this study, we first present an improved genome assembly for the Tasmanian devil using the Chicago library method. Compared with the original reference genome, our new assembly reduces the number of scaffolds (from 35,975 to 10,010) and increases the scaffold N90 (from 0.101 to 2.164 Mb). Second, we assess the performance of four contemporary genomic methods for inferring population size history (PSMC, MSMC, SMC++, Stairway Plot), using the two devil genome assemblies as well as simulated, artificially fragmented genomes that approximate the hypothesized demographic history of Tasmanian devils. We demonstrate that each method is robust to assembly quality, producing similar estimates of Ne when simulated genomes were fragmented into up to 5,000 scaffolds. Overall, methods reliant on the SMC are most reliable between ∼300 generations before present (gbp) and 100 kgbp, whereas methods exclusively reliant on the site-frequency spectrum are most reliable between the present and 30 gbp. Our results suggest that when used in concert, genomic methods for reconstructing species' effective population size histories 1) can be applied to nonmodel organisms without highly contiguous reference genomes, and 2) are capable of detecting independently documented effects of historical geological events.

Infectious disease and sickness behaviour: tumour progression affects interaction patterns and social network structure in wild Tasmanian devils.

Infectious diseases, including transmissible cancers, can have a broad range of impacts on host behaviour, particularly in the latter stages of disease progression. However, the difficulty of early diagnoses makes the study of behavioural influences of disease in wild animals a challenging task. Tasmanian devils (Sarcophilus harrisii) are affected by a transmissible cancer, devil facial tumour disease (DFTD), in which tumours are externally visible as they progress. Using telemetry and mark-recapture datasets, we quantify the impacts of cancer progression on the behaviour of wild devils by assessing how interaction patterns within the social network of a population change with increasing tumour load. The progression of DFTD negatively influences devils' likelihood of interaction within their network. Infected devils were more active within their network late in the mating season, a pattern with repercussions for DFTD transmission. Our study provides a rare opportunity to quantify and understand the behavioural feedbacks of disease in wildlife and how they may affect transmission and population dynamics in general.

Comparative landscape genetics reveals differential effects of environment on host and pathogen genetic structure in Tasmanian devils (Sarcophilus harrisii) and their transmissible tumour.

Genetic structure in host species is often used to predict disease spread. However, host and pathogen genetic variation may be incongruent. Understanding landscape factors that have either concordant or divergent influence on host and pathogen genetic structure is crucial for wildlife disease management. Devil facial tumour disease (DFTD) was first observed in 1996 and has spread throughout almost the entire Tasmanian devil geographic range, causing dramatic population declines. Whereas DFTD is predominantly spread via biting among adults, devils typically disperse as juveniles, which experience low DFTD prevalence. Thus, we predicted little association between devil and tumour population structure and that environmental factors influencing gene flow differ between devils and tumours. We employed a comparative landscape genetics framework to test the influence of environmental factors on patterns of isolation by resistance (IBR) and isolation by environment (IBE) in devils and DFTD. Although we found evidence for broad-scale costructuring between devils and tumours, we found no relationship between host and tumour individual genetic distances. Further, the factors driving the spatial distribution of genetic variation differed for each. Devils exhibited a strong IBR pattern driven by major roads, with no evidence of IBE. By contrast, tumours showed little evidence for IBR and a weak IBE pattern with respect to elevation in one of two tumour clusters we identify herein. Our results warrant caution when inferring pathogen spread using host population genetic structure and suggest that reliance on environmental barriers to host connectivity may be ineffective for managing the spread of wildlife diseases. Our findings demonstrate the utility of comparative landscape genetics for identifying differential factors driving host dispersal and pathogen transmission.

Individual and temporal variation in pathogen load predicts long-term impacts of an emerging infectious disease.

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, whereas DFTD faded 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 the DFTD-induced devil extinction, with crucial management ramifications. Accounting for individual-level disease progression and the long-term outcome of devil-DFTD interactions 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.

Conserving adaptive potential: lessons from Tasmanian devils and their transmissible cancer.

Maintenance of adaptive genetic variation has long been a goal of management of natural populations, but only recently have genomic tools allowed identification of specific loci associated with fitness-related traits in species of conservation concern. This raises the possibility of managing for genetic variation directly relevant to specific threats, such as those due to climate change or emerging infectious disease. Tasmanian devils (Sarcophilus harrisii) face the threat of a transmissible cancer, devil facial tumor disease (DFTD), that has decimated wild populations and led to intensive management efforts. Recent discoveries from genomic and modeling studies reveal how natural devil populations are responding to DFTD, and can inform management of both captive and wild devil populations. Notably, recent studies have documented genetic variation for disease-related traits and rapid evolution in response to DFTD, as well as potential mechanisms for disease resistance such as immune response and tumor regression in wild devils. Recent models predict dynamic persistence of devils with or without DFTD under a variety of modeling scenarios, although at much lower population densities than before DFTD emerged, contrary to previous predictions of extinction. As a result, current management that focuses on captive breeding and release for maintaining genome-wide genetic diversity or demographic supplementation of populations could have negative consequences. Translocations of captive devils into wild populations evolving with DFTD can cause outbreeding depression and/or increases in the force of infection and thereby the severity of the epidemic, and we argue that these risks outweigh any benefits of demographic supplementation in wild populations. We also argue that genetic variation at loci associated with DFTD should be monitored in both captive and wild populations, and that as our understanding of DFTD-related genetic variation improves, considering genetic management approaches to target this variation is warranted in developing conservation strategies for Tasmanian devils.

Sex bias in ability to cope with cancer: Tasmanian devils and facial tumour disease.

Knowledge of the ecological dynamics between hosts and pathogens during the initial stages of disease emergence is crucial to understanding the potential for evolution of new interspecific interactions. Tasmanian devil (Sarcophilus harrisii) populations have declined precipitously owing to infection by a transmissible cancer (devil facial tumour disease, DFTD) that emerged approximately 20 years ago. Since the emergence of DFTD, and as the disease spreads across Tasmania, the number of devils has dropped up to 90% across 80% of the species's distributional range. As a result, the disease is expected to act as a strong selective force on hosts to develop mechanisms of tolerance and/or resistance to the infection. We assessed the ability of infected devils to cope with infection, which translates into host tolerance to the cancer, by using the reaction norm of the individual body condition by tumour burden. We found that body condition of infected hosts is negatively affected by cancer progression. Males and females presented significant differences in their tolerance levels to infection, with males suffering declines of up to 25% of their body condition, in contrast to less than 5% in females. Sex-related differences in tolerance to cancer progression may select for changes in life-history strategies of the host and could also alter the selective environment for the tumours.

Top carnivore decline has cascading effects on scavengers and carrion persistence.

Top carnivores have suffered widespread global declines, with well-documented effects on mesopredators and herbivores. We know less about how carnivores affect ecosystems through scavenging. Tasmania's top carnivore, the Tasmanian devil (Sarcophilus harrisii), has suffered severe disease-induced population declines, providing a natural experiment on the role of scavenging in structuring communities. Using remote cameras and experimentally placed carcasses, we show that mesopredators consume more carrion in areas where devils have declined. Carcass consumption by the two native mesopredators was best predicted by competition for carrion, whereas consumption by the invasive mesopredator, the feral cat (Felis catus), was better predicted by the landscape-level abundance of devils, suggesting a relaxed landscape of fear where devils are suppressed. Reduced discovery of carcasses by devils was balanced by the increased discovery by mesopredators. Nonetheless, carcasses persisted approximately 2.6-fold longer where devils have declined, highlighting their importance for rapid carrion removal. The major beneficiary of increased carrion availability was the forest raven (Corvus tasmanicus). Population trends of ravens increased 2.2-fold from 1998 to 2017, the period of devil decline, but this increase occurred Tasmania-wide, making the cause unclear. This case study provides a little-studied potential mechanism for mesopredator release, with broad relevance to the vast areas of the world that have suffered carnivore declines.

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.

Infection of the fittest: devil facial tumour disease has greatest effect on individuals with highest reproductive output.

Emerging infectious diseases rarely affect all members of a population equally and determining how individuals' susceptibility to infection is related to other components of their fitness is critical to understanding disease impacts at a population level and for predicting evolutionary trajectories. We introduce a novel state-space model framework to investigate survival and fecundity of Tasmanian devils (Sarcophilus harrisii) affected by a transmissible cancer, devil facial tumour disease. We show that those devils that become host to tumours have otherwise greater fitness, with higher survival and fecundity rates prior to disease-induced death than non-host individuals that do not become infected, although high tumour loads lead to high mortality. Our finding that individuals with the greatest reproductive value are those most affected by the cancer demonstrates the need to quantify both survival and fecundity in context of disease progression for understanding the impact of disease on wildlife populations.