Low MHC class II diversity in the Tasmanian devil (Sarcophilus harrisii)

Plethora of New Marsupial Genomes Informs Our Knowledge of Marsupial MHC Class II

The major histocompatibility complex (MHC) plays a vital role in the vertebrate immune system due to its role in infection, disease and autoimmunity, or recognition of "self". The marsupial MHC class II genes show divergence from eutherian MHC class II genes and are a unique taxon of therian mammals that give birth to altricial and immunologically naive young providing an opportune study system for investigating evolution of the immune system. Additionally, the MHC in marsupials has been implicated in disease associations, including susceptibility to Chlamydia pecorum infection in koalas. Due to the complexity of the gene family, automated annotation is not possible so here we manually annotate 384 class II MHC genes in 29 marsupial species. We find losses of key components of the marsupial MHC repertoire in the Dasyuromorphia order and the Pseudochiridae family. We perform PGLS analysis to show the gene losses we find are true gene losses and not artifacts of unresolved genome assembly. We investigate the associations between the number of loci and life history traits, including lifespan and reproductive output in lineages of marsupials and hypothesize that gene loss may be linked to the energetic cost and tradeoffs associated with pregnancy and reproduction. We found support for litter size being a significant predictor of the number of DBA and DBB loci, indicating a tradeoff between the energetic requirements of immunity and reproduction. Additionally, we highlight the increased susceptibility of Dasyuridae species to neoplasia and a potential link to MHC gene loss. Finally, these annotations provide a valuable resource to the immunogenetics research community to move forward and further investigate diversity in MHC genes in marsupials.

The complex effects of modern oncogenic environments on the fitness, evolution and conservation of wildlife species

Growing evidence indicates that human activities are causing cancer rates to rise in both human and wildlife populations. This is due to the inability of ancestral anti-cancer defences to cope with modern environmental risks. The evolutionary mismatch between modern oncogenic risks and evolved cancer defences has far-reaching effects on various biological aspects at different timeframes, demanding a comprehensive study of the biology and evolutionary ecology of the affected species. Firstly, the increased activation of anti-cancer defences leads to excessive energy expenditure, affecting other biological functions and potentially causing health issues like autoimmune diseases. Secondly, tumorigenesis itself can impact important fitness-related parameters such as competitiveness, predator evasion, resistance to parasites, and dispersal capacity. Thirdly, rising cancer risks can influence the species' life-history traits, often favoring early reproduction to offset fitness costs associated with cancer. However, this strategy has its limits, and it may not ensure the sustainability of the species if cancer risks continue to rise. Lastly, some species may evolve additional anti-cancer defences, with uncertain consequences for their biology and future evolutionary path. In summary, we argue that the effects of increased exposure to cancer-causing substances on wildlife are complex, ranging from immediate responses to long-term evolutionary changes. Understanding these processes, especially in the context of conservation biology, is urgently needed.

Temporal dynamics of woolly mammoth genome erosion prior to extinction

A number of species have recently recovered from near-extinction. Although these species have avoided the immediate extinction threat, their long-term viability remains precarious due to the potential genetic consequences of population declines, which are poorly understood on a timescale beyond a few generations. Woolly mammoths (Mammuthus primigenius) became isolated on Wrangel Island around 10,000 years ago and persisted for over 200 generations before becoming extinct around 4,000 years ago. To study the evolutionary processes leading up to the mammoths' extinction, we analyzed 21 Siberian woolly mammoth genomes. Our results show that the population recovered quickly from a severe bottleneck and remained demographically stable during the ensuing six millennia. We find that mildly deleterious mutations gradually accumulated, whereas highly deleterious mutations were purged, suggesting ongoing inbreeding depression that lasted for hundreds of generations. The time-lag between demographic and genetic recovery has wide-ranging implications for conservation management of recently bottlenecked populations.

Best genome sequencing strategies for annotation of complex immune gene families in wildlife

BACKGROUND

The biodiversity crisis and increasing impact of wildlife disease on animal and human health provides impetus for studying immune genes in wildlife. Despite the recent boom in genomes for wildlife species, immune genes are poorly annotated in nonmodel species owing to their high level of polymorphism and complex genomic organisation. Our research over the past decade and a half on Tasmanian devils and koalas highlights the importance of genomics and accurate immune annotations to investigate disease in wildlife. Given this, we have increasingly been asked the minimum levels of genome quality required to effectively annotate immune genes in order to study immunogenetic diversity. Here we set out to answer this question by manually annotating immune genes in 5 marsupial genomes and 1 monotreme genome to determine the impact of sequencing data type, assembly quality, and automated annotation on accurate immune annotation.

RESULTS

Genome quality is directly linked to our ability to annotate complex immune gene families, with long reads and scaffolding technologies required to reassemble immune gene clusters and elucidate evolution, organisation, and true gene content of the immune repertoire. Draft-quality genomes generated from short reads with HiC or 10× Chromium linked reads were unable to achieve this. Despite mammalian BUSCOv5 scores of up to 94.1% amongst the 6 genomes, automated annotation pipelines incorrectly annotated up to 59% of manually annotated immune genes regardless of assembly quality or method of automated annotation.

CONCLUSIONS

Our results demonstrate that long reads and scaffolding technologies, alongside manual annotation, are required to accurately study the immune gene repertoire of wildlife species.

Restoring faith in conservation action: Maintaining wild genetic diversity through the Tasmanian devil insurance program

Conservation breeding programs aim to maintain 90% wild genetic diversity, but rarely assess functional diversity. Here, we compare both genome-wide and functional diversity (in over 500 genes) of Tasmanian devils (Sarcophilus harrisii) within the insurance metapopulation and across the species’ range (64,519 km²). Populations have declined by 80% since 1996 due to a contagious cancer, devil facial tumor disease (DFTD). However, predicted local extinctions have not occurred. Recent suggestions of selection for “resistance” alleles in the wild precipitated concerns that insurance population devils may be unsuitable for translocations. Using 830 wild samples collected at 31 locations between 2012 and 2021, and 553 insurance metapopulation devils, we show that the insurance metapopulation is representative of current wild genetic diversity. Allele frequencies at DFTD-associated loci were not substantially different between captive and wild devils. Methods presented here are valuable for others investigating evolutionary potential in threatened species, particularly ones under significant selective pressures.

•

Developed target capture to assess functional diversity at over 500 genes

•

Fine-scale structure exists in the genetically depauperate Tasmanian devil

•

Insurance metapopulation is representative of wild genetic diversity

•

Allele frequencies at disease-associated loci were similar in captivity to the wild

Differences in constitutive innate immunity between divergent Australian marsupials

Understanding immunity in wildlife populations is important from both One Health and conservation perspectives. The constitutive innate immune system is the first line of defence against pathogens, and comparisons among taxa can test the impact of evolution and life history on immune function. We investigated serum bacterial killing ability (BKA) of five marsupial species that employ varying life history strategies, demonstrated to influence immunity in other vertebrates. The brushtail possum and eastern grey kangaroo had the greatest BKA, while ringtail possums and koalas had the least. These differences were independent of social structure, captivity status and phylogeny, but were associated with diet and body size. Sex and disease status had no effect on BKA in koalas, however potential for differences between wild and captive koalas warrants further investigation. The current study has provided a foundation for future investigations into how adaptive and innate immunity interact in marsupials from an eco-evolutionary perspective.

Genetic diversity and disease: The past, present, and future of an old idea

Why do infectious diseases erupt in some host populations and not others? This question has spawned independent fields of research in evolution, ecology, public health, agriculture, and conservation. In the search for environmental and genetic factors that predict variation in parasitism, one hypothesis stands out for its generality and longevity: genetically homogeneous host populations are more likely to experience severe parasitism than genetically diverse populations. In this perspective piece, I draw on overlapping ideas from evolutionary biology, agriculture, and conservation to capture the far-reaching implications of the link between genetic diversity and disease. I first summarize the development of this hypothesis and the results of experimental tests. Given the convincing support for the protective effect of genetic diversity, I then address the following questions: (1) Where has this idea been put to use, in a basic and applied sense, and how can we better use genetic diversity to limit disease spread? (2) What new hypotheses does the established disease-diversity relationship compel us to test? I conclude that monitoring, preserving, and augmenting genetic diversity is one of our most promising evolutionarily informed strategies for buffering wild, domesticated, and human populations against future outbreaks.

Improved high-throughput MHC typing for non-model species using long-read sequencing

The major histocompatibility complex (MHC) plays a critical role in the vertebrate immune system. Accurate MHC typing is critical to understanding not only host fitness and disease susceptibility, but also the mechanisms underlying host‐pathogen co‐evolution. However, due to the high degree of gene duplication and diversification of MHC genes, it is often technically challenging to accurately characterise MHC genetic diversity in non‐model species. Here we conducted a systematic review to identify common issues associated with current widely used MHC typing approaches. Then to overcome these challenges, we developed a long‐read based MHC typing method along with a new analysis pipeline. Our approach enables the sequencing of fully phased MHC alleles spanning all key functional domains and the separation of highly similar alleles as well as the removal of technical artefacts such as PCR heteroduplexes and chimeras. Using this approach, we performed population‐scale MHC typing in the Tasmanian devil (Sarcophilus harrisii), revealing previously undiscovered MHC functional diversity in this endangered species. Our new method provides a better solution for addressing research questions that require high MHC typing accuracy. Since the method is not limited by species or the number of genes analysed, it will be applicable for studying not only the MHC but also other complex gene families.

Cancer risk landscapes: A framework to study cancer in ecosystems

Cancer is a family of diseases that has been documented in most metazoan species and ecosystems. Human induced environmental changes are increasingly exposing wildlife to carcinogenic risk factors, and negative repercussions on ecosystems and on the conservation of endangered species are already been observed. It is therefore of key importance to understand the spatiotemporal variability of those risk factors and how they interact with the biosphere to mitigate their effects. Here we introduce the concept of cancer risk landscape that can be applied to understand how species are exposed to, interact with, and modify cancer risk factors. With this publication we aim to provide a framework in order to stimulate a discussion on how to mitigate cancer-causing risk factors.

Annotation of immune genes in the extinct thylacine (Thylacinus cynocephalus)

Advances in genome sequencing technology have enabled genomes of extinct species to be sequenced. However, given the fragmented nature of these genome assemblies, it is not clear whether it is possible to comprehensively annotate highly variable and repetitive genes such as those involved in immunity. As such, immune genes have only been investigated in a handful of extinct genomes, mainly in human lineages. In 2018 the genome of the thylacine (Thylacinus cynocephalus), a carnivorous marsupial from Tasmania that went extinct in 1936, was sequenced. Here we attempt to characterise the immune repertoire of the thylacine and determine similarity to its closest relative with a genome available, the Tasmanian devil (Sarcophilus harrisii), as well as other marsupials. Members from all major immune gene families were identified. However, variable regions could not be characterised, and complex families such as the major histocompatibility complex (MHC) were highly fragmented and located across multiple small scaffolds. As such, at a gene level we were unable to reconstruct full-length coding sequences for the majority of thylacine immune genes. Despite this, we identified genes encoding functionally important receptors and immune effector molecules, which suggests the functional capacity of the thylacine immune system was similar to other mammals. However, the high number of partial immune gene sequences identified limits our ability to reconstruct an accurate picture of the thylacine immune repertoire.

Characterization of reproductive gene diversity in the endangered Tasmanian devil

Interindividual variation at genes known to play a role in reproduction may impact reproductive fitness. The Tasmanian devil is an endangered Australian marsupial with low genetic diversity. Recent work has shown concerning declines in productivity in both wild and captive populations over time. Understanding whether functional diversity exists at reproductive genes in the Tasmanian devil is a key first step in identifying genes that may influence productivity. We characterized single nucleotide polymorphisms (SNPs) at 214 genes involved in reproduction in 37 Tasmanian devils. Twenty genes contained nonsynonymous substitutions, with genes involved in embryogenesis, fertilization and hormonal regulation of reproduction displaying greater numbers of nonsynonymous SNPs than synonymous SNPs. Two genes, ADAMTS9 and NANOG, showed putative signatures of balancing selection indicating that natural selection is maintaining diversity at these genes despite the species exhibiting low overall levels of genetic diversity. We will use this information in future to examine the interplay between reproductive gene variation and reproductive fitness in Tasmanian devil populations.

Transmissible cancers in mammals and bivalves: How many examples are there?: Predictions indicate widespread occurrence

Transmissible cancers are elusive and understudied parasitic life forms caused by malignant clonal cells (nine lineages are known so far). They emerge by completing sequential steps that include breaking cell cooperation, evade anti-cancer defences and shedding cells to infect new hosts. Transmissible cancers impair host fitness, and their importance as selective force is likely largely underestimated. It is, therefore, crucial to determine how common they might be in the wild. Here, we draw a parallel between the steps required for a transmissible cancer to emerge and the steps required for an intelligent civilisation to emerge in the Milky Way using a modified Drake equation. Using numerical analyses, we estimate the potential number of extant marine and bivalve species in which transmissible cancers might exist. Our results suggest that transmissible cancers are more common than expected, and that new lineages can be found by screening a large number of species.

A Tasmanian devil breeding program to support wild recovery

Tasmanian devils are threatened in the wild by devil facial tumour disease: a transmissible cancer with a high fatality rate. In response, the Save the Tasmanian Devil Program (STDP) established an 'insurance population' to enable the preservation of genetic diversity and natural behaviours of devils. This breeding program includes a range of institutions and facilities, from zoo-based intensive enclosures to larger, more natural environments, and a strategic approach has been required to capture and maintain genetic diversity, natural behaviours and to ensure reproductive success. Laboratory-based research, particularly genetics, in tandem with adaptive management has helped the STDP reach its goals, and has directly contributed to the conservation of the species in the wild. Here we review this work and show that the Tasmanian devil breeding program is a powerful example of how genetic research can be used to understand and improve reproductive success in a threatened species.

MHC Class II variability in bare-nosed wombats (Vombatus ursinus)

Studies of gene diversity are used to investigate population dynamics, including immunological fitness. Aside from the selection of an appropriate gene target, an important factor that underpins these studies is the ability to recover viable DNA samples from native animals that are protected, threatened or difficult to sample or locate such as the bare-nosed wombat (Vombatus ursinus). In this study, we used genomic DNA extracted from muscle tissue samples and also identified the optimal method to extract DNA from fresh wombat scat samples to enable further analyses to be performed using non-invasive techniques. The DNA was probed via the polymerase chain reaction using previously targeted marsupial Major Histocompatibility Complex (MHC) gene primers. These genes are highly variable and associated with binding and presentation of pathogens in the immune system. Twenty-three unique MHC Class II DAB V. ursinus gene sequences were translated to 21 unique predicted peptide sequences from 34 individual tissue or scat samples. Vombatus ursinus MHC Class II DAB gene and peptide sequences were most similar to other marsupial DNA and peptide sequences. Further analysis also indicated the likelihood of MHC Class II DAB family membership through motif identification. Additional sampling is required to assess the full level of diversity of MHC Class II DAB genes among V. ursinus populations; however, this study is the first to identify MHC genes in a wombat and will advance immunological and disease studies of the species.

MHC-associated mate choice under competitive conditions in captive versus wild Tasmanian devils

Mate choice contributes to driving evolutionary processes when animals choose breeding partners that confer genetic advantages to offspring, such as increased immunocompetence. The major histocompatibility complex (MHC) is an important group of immunological molecules, as MHC antigens bind and present foreign peptides to T-cells. Recent studies suggest that mates may be selected based on their MHC profile, leading to an association between an individual’s MHC diversity and their breeding success. In conservation, it may be important to consider mate choice in captive breeding programs, as this mechanism may improve reproductive rates. We investigated the reproductive success of Tasmanian devils in a group housing facility to determine whether increased MHC-based heterozygosity led individuals to secure more mating partners and produce more offspring. We also compared the breeding success of captive females to a wild devil population. MHC diversity was quantified using 12 MHC-linked microsatellite markers, including 11 previously characterized markers and one newly identified marker. Our analyses revealed that there was no relationship between MHC-linked heterozygosity and reproductive success either in captivity or the wild. The results of this study suggest that, for Tasmanian devils, MHC-based heterozygosity does not produce greater breeding success and that no specific changes to current captive management strategies are required with respect to preserving MHC diversity.

Fecal Viral Diversity of Captive and Wild Tasmanian Devils Characterized Using Virion-Enriched Metagenomics and Metatranscriptomics

The Tasmanian devil is an endangered carnivorous marsupial threatened by devil facial tumor disease (DFTD). While research on DFTD has been extensive, little is known about viruses in devils and whether any are of potential conservation relevance for this endangered species. Using both metagenomics based on virion enrichment and sequence-independent amplification (virion-enriched metagenomics) and metatranscriptomics based on bulk RNA sequencing, we characterized and compared the fecal viromes of captive and wild devils. A total of 54 fecal samples collected from two captive and four wild populations were processed for virome characterization using both approaches. In total, 24 novel marsupial-related viruses, comprising a sapelovirus, astroviruses, rotaviruses, picobirnaviruses, parvoviruses, papillomaviruses, polyomaviruses, and a gammaherpesvirus, were identified, as well as known mammalian pathogens such as rabbit hemorrhagic disease virus 2. Captive devils showed significantly lower viral diversity than wild devils. Comparison of the two virus discovery approaches revealed substantial differences in the number and types of viruses detected, with metatranscriptomics better suited for RNA viruses and virion-enriched metagenomics largely identifying more DNA viruses. Thus, the viral communities revealed by virion-enriched metagenomics and metatranscriptomics were not interchangeable and neither approach was able to detect all viruses present. An integrated approach using both virion-enriched metagenomics and metatranscriptomics constitutes a powerful tool for obtaining a complete overview of both the taxonomic and functional profiles of viral communities within a sample.IMPORTANCE The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

Fecal Viral Diversity of Captive and Wild Tasmanian Devils Characterized Using Virion-Enriched Metagenomics and Metatranscriptomics

The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

The Tasmanian devil is an endangered carnivorous marsupial threatened by devil facial tumor disease (DFTD). While research on DFTD has been extensive, little is known about viruses in devils and whether any are of potential conservation relevance for this endangered species. Using both metagenomics based on virion enrichment and sequence-independent amplification (virion-enriched metagenomics) and metatranscriptomics based on bulk RNA sequencing, we characterized and compared the fecal viromes of captive and wild devils. A total of 54 fecal samples collected from two captive and four wild populations were processed for virome characterization using both approaches. In total, 24 novel marsupial-related viruses, comprising a sapelovirus, astroviruses, rotaviruses, picobirnaviruses, parvoviruses, papillomaviruses, polyomaviruses, and a gammaherpesvirus, were identified, as well as known mammalian pathogens such as rabbit hemorrhagic disease virus 2. Captive devils showed significantly lower viral diversity than wild devils. Comparison of the two virus discovery approaches revealed substantial differences in the number and types of viruses detected, with metatranscriptomics better suited for RNA viruses and virion-enriched metagenomics largely identifying more DNA viruses. Thus, the viral communities revealed by virion-enriched metagenomics and metatranscriptomics were not interchangeable and neither approach was able to detect all viruses present. An integrated approach using both virion-enriched metagenomics and metatranscriptomics constitutes a powerful tool for obtaining a complete overview of both the taxonomic and functional profiles of viral communities within a sample.

IMPORTANCE The Tasmanian devil is an iconic Australian marsupial that has suffered an 80% population decline due to a contagious cancer, devil facial tumor disease, along with other threats. Until now, viral discovery in this species has been confined to one gammaherpesvirus (dasyurid herpesvirus 2 [DaHV-2]), for which captivity was identified as a significant risk factor. Our discovery of 24 novel marsupial-associated RNA and DNA viruses, and that viral diversity is lower in captive than in wild devils, has greatly expanded our knowledge of gut-associated viruses in devils and provides important baseline information that will contribute to the conservation and captive management of this endangered species. Our results also revealed that a combination of virion-enriched metagenomics and metatranscriptomics may be a more comprehensive approach for virome characterization than either method alone. Our results thus provide a springboard for continuous improvements in the way we study complex viral communities.

Complex problems need detailed solutions: Harnessing multiple data types to inform genetic management in the wild

For bottlenecked populations of threatened species, supplementation often leads to improved population metrics (genetic rescue), provided that guidelines can be followed to avoid negative outcomes. In cases where no "ideal" source populations exist, or there are other complicating factors such as prevailing disease, the benefit of supplementation becomes uncertain. Bringing multiple data and analysis types together to plan genetic management activities can help. Here, we consider three populations of Tasmanian devil, Sarcophilus harrisii, as candidates for genetic rescue. Since 1996, devil populations have been severely impacted by devil facial tumour disease (DFTD), causing significant population decline and fragmentation. Like many threatened species, the key threatening process for devils cannot currently be fully mitigated, so species management requires a multifaceted approach. We examined diversity of 31 putatively neutral and 11 MHC-linked microsatellite loci of three remnant wild devil populations (one sampled at two time-points), alongside computational diversity projections, parameterized by field data from DFTD-present and DFTD-absent sites. Results showed that populations had low diversity, connectivity was poor, and diversity has likely decreased over the last decade. Stochastic simulations projected further diversity losses. For a given population size, the effects of DFTD on population demography (including earlier age at death and increased female productivity) did not impact diversity retention, which was largely driven by final population size. Population sizes ≥500 (depending on the number of founders) were necessary for maintaining diversity in otherwise unmanaged populations, even if DFTD is present. Models indicated that smaller populations could maintain diversity with ongoing immigration. Taken together, our results illustrate how multiple analysis types can be combined to address complex population genetic challenges.

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.

Lessons learnt from the Tasmanian devil facial tumour regarding immune function in cancer

Genetic and genomic technologies have facilitated a greater understanding of the Tasmanian devil immune system and the origins, evolution and spread of devil facial tumour disease (DFTD). DFTD is a contagious cancer that has caused significant declines in devil populations across Tasmania. Immune responses to DFTD are rarely detected, allowing the cancer to pass between individuals and proliferate unimpeded. Early immunosenscence in devils appears to decrease anti-tumour immunity in older animals compared to younger animals, which may increase susceptibility to DFTD and explain high DFTD prevalence in this age group. Devils also have extremely low major histocompatibility complex (MHC) diversity, and multiple alleles are shared with the tumour, lowering histocompatibility barriers which may have contributed to DFTD evolution. DFTD actively evades immune attack by down-regulating cell-surface MHC I molecules, making it effectively invisible to the immune system. Altered MHC I profiles should activate natural killer (NK) cell anti-tumour responses, but these are absent in DFTD infection. Recent immunisation and immunotherapy using modified DFTD cells has induced an anti-DFTD immune response and regression of DFTD in some devils. Knowledge gained from immune responses to a transmissible cancer in devils will ultimately reveal useful insights into immunity to cancer in humans and other species.

Transcriptome sequencing of the long-nosed bandicoot (Perameles nasuta) reveals conservation and innovation of immune genes in the marsupial order Peramelemorphia

Bandicoots are omnivorous marsupials of the order Peramelemorphia. Conservation concerns and their unique biological characteristics suggest peramelomorphs are worthy research subjects, but knowledge of their genetics and immunology has lagged behind that of other high-profile marsupials. Here, we characterise the transcriptome of the long-nose bandicoot (Perameles nasuta), the first high-throughput data set from any peramelomorph. We investigate the immune gene repertoire of the bandicoot, with a focus on key immune gene families, and compare to previously characterised marsupial and mammalian species. We find that the immune gene complement in bandicoot is often conserved with respect to other marsupials; however, the diversity of expressed transcripts in several key families, such as major histocompatibility complex, T cell receptor μ and natural killer cell receptors, appears greater in the bandicoot than other Australian marsupials, including devil and koala. This transcriptome is an important first step for future studies of bandicoots and the bilby, allowing for population level analysis and construction of bandicoot-specific immunological reagents and assays. Such studies will be critical to understanding the immunology and physiology of Peramelemorphia and to inform the conservation of these unique marsupials.

The role of MHC genes in contagious cancer: the story of Tasmanian devils

The Tasmanian devil, a marsupial species endemic to the island of Tasmania, harbours two contagious cancers, Devil Facial Tumour 1 (DFT1) and Devil Facial Tumour 2 (DFT2). These cancers pass between individuals in the population via the direct transfer of tumour cells, resulting in the growth of large tumours around the face and neck of affected animals. While these cancers are rare, a contagious cancer also exists in dogs and five contagious cancers circulate in bivalves. The ability of tumour cells to emerge and transmit in mammals is surprising as these cells are an allograft and should be rejected due to incompatibility between Major Histocompatibility Complex (MHC) genes. As such, considerable research has focused on understanding how DFT1 cells evade the host immune system with particular reference to MHC molecules. This review evaluates the role that MHC class I expression and genotype plays in allowing DFT1 to circumvent histocompatibility barriers in Tasmanian devils. We also examine recent research that suggests that Tasmanian devils can mount an immune response to DFT1 and may form the basis of a protective vaccine against the tumour.

Low genetic variation in the MHC class II DRB gene and MHC-linked microsatellites in endangered island populations of the leopard cat (Prionailurus bengalensis) in Japan

Isolated populations of the leopard cat (Prionailurus bengalensis) on Tsushima and Iriomote islands in Japan are classified as subspecies P. b. euptilurus and P. b. iriomotensis, respectively. Because both populations have decreased to roughly 100, an understanding of their genetic diversity is essential for conservation. We genotyped MHC class II DRB exon 2 and MHC-linked microsatellite loci to evaluate the diversity of MHC genes in the Tsushima and Iriomote cat populations. We detected ten and four DRB alleles in these populations, respectively. A phylogenetic analysis showed DRB alleles from both populations to be closely related to those in other felid DRB lineages, indicating trans-species polymorphism. The MHC-linked microsatellites were more polymorphic in the Tsushima than in the Iriomote population. The MHC diversity of both leopard cat populations is much lower than in the domestic cat populations on these islands, probably due to inbreeding associated with founder effects, geographical isolation, or genetic drift. Our results predict low resistance of the two endangered populations to new pathogens introduced to the islands.

Cancer as a contagious disease

Contagious cancers are malignant cells that are physically transferred between individuals as a natural allograft, forming new clonal tumours. These cancers are highly unusual, but have emerged in 2 mammalian species, the dog and the Tasmanian devil, as well as 4 species of bivalve. The transfer of malignant cells in mammals should initiate a robust immune response and although invertebrates have a less complex immune system, these species still have mechanisms that should prevent engraftment and protect against cellular parasitism. Here the naturally occurring contagious cancers are reviewed to determine what features are important and necessary for the emergence and spread of these types of cancer, with a focus on the mammalian contagious cancers and how they successfully cross histocompatibility barriers.

Major Histocompatibility Complex Class II in the red-tailed phascogale (Phascogale calura)

Diversity in major histocompatibility complex (MHC) genes can be correlated with the level of immunological fitness of an individual or group of individuals. This study tested published primer sets designed to amplify fragments of the MHC Class II DAB and DBB genes to amplify the equivalent gene fragments in red-tailed phascogales (Phascogale calura). Seventeen genomic DNA samples extracted from phascogale muscle tissue were used to amplify the initial DAB and DBB fragments; however, only DAB PCR proved successful. The fragments were 172 bp in length between the primers and had a high level of identity to other known marsupial MHC Class II DAB gene sequences (89–98%), including those of the koala (Phascolarctos cinereus), Tasmanian devil (Sarcophilus harrisii), common brushtail possum (Trichosurus vulpecula) and several wallaby species. Multiple sequence alignment revealed limited variability of MHC Class II genes between the individuals, but eight individual sequences in total. Genomic DNA was subsequently extracted from three fresh red-tailed phascogale scat samples and DAB fragments successfully amplified. The technique will allow for red-tailed phascogales to be sampled non-invasively in the wild and to determine the level of MHC diversity among individuals in the population.

The identification of immune genes in the milk transcriptome of the Tasmanian devil (Sarcophilus harrisii)

Tasmanian devil (Sarcophilus harrisii) pouch young, like other marsupials, are born underdeveloped and immunologically naïve, and are unable to mount an adaptive immune response. The mother’s milk provides nutrients for growth and development as well as providing passive immunity. To better understand immune response in this endangered species, we set out to characterise the genes involved in passive immunity by sequencing and annotating the transcriptome of a devil milk sample collected during mid-lactation. At mid-lactation we expect the young to have heightened immune responses, as they have emerged from the pouch, encountering new pathogens. A total of 233,660 transcripts were identified, including approximately 17,827 unique protein-coding genes and 846 immune genes. The most highly expressed transcripts were dominated by milk protein genes such as those encoding early lactation protein, late lactation proteins, α-lactalbumin, α-casein and β-casein. There were numerous highly expressed immune genes including lysozyme, whey acidic protein, ferritin and major histocompatibility complex I and II. Genes encoding immunoglobulins, antimicrobial peptides, chemokines and immune cell receptors were also identified. The array of immune genes identified in this study reflects the importance of the milk in providing immune protection to Tasmanian devil young and provides the first insight into Tasmanian devil milk.

Identification and analysis of divergent immune gene families within the Tasmanian devil genome

Background

The Tasmanian devil (Sarcophilus harrisii) is being threatened with extinction in the wild by a disease known as devil facial tumour disease (DFTD). In order to prevent the spread of this disease a thorough understanding of the Tasmanian devil immune system and its response to the disease is required. In 2011 and 2012 two genome sequencing projects of the Tasmania devil were released. This has provided us with the raw data required to begin to investigate the Tasmanian devil immunome in depth. In this study we characterise immune gene families of the Tasmanian devil. We focus on immunoglobulins, T cell receptors and cytokine families.

Results

We identify and describe 119 cytokines including 40 interleukins, 39 chemokines, 8 interferons, 18 tumour necrosis family cytokines and 14 additional cytokines. Constant regions for immunoglobulins and T cell receptors were also identified. The repertoire of genes in these families was similar to the opossum, however devil specific duplications were seen and orthologs to eutherian genes not previously identified in any marsupial were also identified.

Conclusions

By using multiple data sources as well as targeted search methods, highly divergent genes across the Tasmanian devil immune system were identified and characterised. This understanding will allow for the development of devil specific assays and reagents and allow for future studies into the immune response of the Tasmanian devil immune system to DFTD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2206-9) contains supplementary material, which is available to authorized users.

Lack of genetic diversity across diverse immune genes in an endangered mammal, the Tasmanian devil (Sarcophilus harrisii)

The Tasmanian devil (Sarcophilus harrisii) is threatened with extinction due to the spread of devil facial tumour disease. Polymorphisms in immune genes can provide adaptive potential to resist diseases. Previous studies in diversity at immune loci in wild species have almost exclusively focused on genes of the major histocompatibility complex (MHC); however, these genes only account for a fraction of immune gene diversity. Devils lack diversity at functionally important immunity loci, including MHC and Toll-like receptor genes. Whether there are polymorphisms at devil immune genes outside these two families is unknown. Here, we identify polymorphisms in a wide range of key immune genes, and develop assays to type single nucleotide polymorphisms (SNPs) within a subset of these genes. A total of 167 immune genes were examined, including cytokines, chemokines and natural killer cell receptors. Using genome-level data from ten devils, SNPs within coding regions, introns and 10 kb flanking genes of interest were identified. We found low polymorphism across 167 immune genes examined bioinformatically using whole-genome data. From this data, we developed long amplicon assays to target nine genes. These amplicons were sequenced in 29-220 devils and found to contain 78 SNPs, including eight SNPS within exons. Despite the extreme paucity of genetic diversity within these genes, signatures of balancing selection were exhibited by one chemokine gene, suggesting that remaining diversity may hold adaptive potential. The low functional diversity may leave devils highly vulnerable to infectious disease, and therefore, monitoring and preserving remaining diversity will be critical for the long-term management of this species. Examining genetic variation in diverse immune genes should be a priority for threatened wildlife species. This study can act as a model for broad-scale immunogenetic diversity analysis in threatened species.

Positive selection drives the evolution of a major histocompatibility complex gene in an endangered Mexican salamander species complex

Immune gene evolution can be critical to species survival in the face of infectious disease. In particular, polymorphism in the genes of the major histocompatibility complex (MHC) helps vertebrates combat novel and diverse pathogens by increasing the number of pathogen-derived proteins that can initiate the host's acquired immune response. In this study, we used a combination of presumably adaptive and neutral markers to investigate MHC evolution in populations of five salamander species within the Ambystoma velasci complex, a group consisting of 15 recently diverged species, several of which are endangered. We isolated 31 unique MHC class II β alleles from 75 total individuals from five species in this complex. MHC heterozygosity was significantly lower than expected for all five species, and we found no clear relationship between number of MHC alleles and species range, life history, or level of heterozygosity. We inferred a phylogeny representing the evolutionary history of Ambystoma MHC, with which we found signatures of positive selection on the overall gene, putative peptide-binding residues, and allelic lineages. We identified several instances of trans-species polymorphism, a hallmark of balancing selection observed in other groups of closely related species. In contrast, we did not detect comparable allelic diversity or signatures of selection on neutral loci. Additionally, we identified 17 supertypes among the 44 unique Ambystoma alleles, indicating that these sequences may encode functionally distinct MHC variants. We therefore have strong evidence that positive selection is a major evolutionary force driving patterns of MHC polymorphism in this recently radiated species complex.

Diversity in the Toll-like receptor genes of the Tasmanian devil (Sarcophilus harrisii)

The Tasmanian devil is an endangered marsupial species that has survived several historical bottlenecks and now has low genetic diversity. Here we characterize the Toll-like receptor (TLR) genes and their diversity in the Tasmanian devil. TLRs are a key innate immune gene family found in all animals. Ten TLR genes were identified in the Tasmanian devil genome. Unusually low levels of diversity were found in 25 devils from across Tasmania. We found two alleles at TLR2, TLR3 and TLR6. The other seven genes were monomorphic. The insurance population, which safeguards the species from extinction, has successfully managed to capture all of these TLR alleles, but concerns remain for the long-term survival of this species.

Characterisation of non-classical MHC class I genes in the Tasmanian devil (Sarcophilus harrisii)

The Tasmanian devil (Sarcophilus harrisii) is a carnivorous marsupial that is under threat of extinction due to an unusual transmissible disease called Devil Facial Tumour Disease (DFTD). Previous studies on the classical MHC genes have provided important insights into immune responses in this endangered species; however, so far, very little is known about the non-classical MHC genes of this species, which can also play significant roles in the immune system. Here, we report characterisation of five non-classical class I genes in the Tasmanian devil, including Saha-UD, -UK, -UM, -MR1 and -CD1. Saha-UD has been isolated previously and is known to have low genetic polymorphism, though its categorisation as classical or non-classical gene has remained undetermined. In this study, we observed tissue-specific expression of Saha-UD, suggesting that it is more characteristic of a non-classical gene. Restricted tissue expression patterns were also observed for other genes, with an exception of Saha-MR1 being ubiquitously expressed in all examined tissues. Saha-UK, -UM and -MR1 were found to be genetically monomorphic, while four alleles were found at Saha-CD1 with signs of positive selection detected within the α1 domain. Among the four Saha-CD1 alleles, one predominant allele (Saha-CD1*01) showed a high allele frequency of 0.906 in the Tasmanian devil population, resulting in a low heterozygosity (0.188) at this locus. Alternative splicing takes place in Saha-CD1, giving rise to a full-length transcript and a splice variant lacking intact antigen-binding, β2m-binding, transmembrane and cytoplasmic domains.

Exosomes from CIITA-transfected CT26 cells enhance anti- tumor effects

AIM

To study anti-tumor effects of exosomes from class II transactivator (CIITA) gene transfected CT26 cells.

METHODS

In this study, we established an MHC class II molecule-expressing murine colon cancer cell line (CT26-CIITA) by transduction of the CIITA gene. Immune effects in vitro and tumor protective results in vivo were tested and monitored.

RESULTS

Exosomes from CT26-CIITA cells were found to contain a high level of MHC class II protein. When loaded on dendritic cells (DCs), exosomes from CT26-CIITA cells significantly increased expression of MHC class II molecules, CD86 and CD80, as compared to exosomes from CT26 cells. In vitro assays using co-culture of immunized splenocytes and exosome-loaded DCs demonstrated that CIITA- Exo enhanced splenocyte proliferation and IFN-γ production of CD4+T cells, while inhibiting IL-10 secretion. In addition, compared to exosomes from CT26 cells, CT26-CIITA-derived exosomes induced higher TNF-α and IL-12 mRNA levels. A mouse tumour preventive model showed that CT26-CIITA derived exosomes significantly inhibited tumour growth in a dose-dependent manner and significantly prolonged the survival time of tumour- bearing mice.

CONCLUSION

Our findings indicate that CT26-CIITA-released exosomes are more efficient to induce anti-tumour immune responses, suggesting a potential role of MHC class II-containing tumour exosomes as cancer vaccine candidates.

Drift rather than selection dominates MHC class II allelic diversity patterns at the biogeographical range scale in natterjack toads Bufo calamita

Study of major histocompatibility complex (MHC) loci has gained great popularity in recent years, partly due to their function in protecting vertebrates from infections. This is of particular interest in amphibians on account of major threats many species face from emergent diseases such as chytridiomycosis. In this study we compare levels of diversity in an expressed MHC class II locus with neutral genetic diversity at microsatellite loci in natterjack toad (Bufo (Epidalea) calamita) populations across the whole of the species' biogeographical range. Variation at both classes of loci was high in the glacial refugium areas (REF) and much lower in postglacial expansion areas (PGE), especially in range edge populations. Although there was clear evidence that the MHC locus was influenced by positive selection in the past, congruence with the neutral markers suggested that historical demographic events were the main force shaping MHC variation in the PGE area. Both neutral and adaptive genetic variation declined with distance from glacial refugia. Nevertheless, there were also some indications from differential isolation by distance and allele abundance patterns that weak effects of selection have been superimposed on the main drift effect in the PGE zone.

The Vertebrate Genome Annotation browser 10 years on

The Vertebrate Genome Annotation (VEGA) database (http://vega.sanger.ac.uk), initially designed as a community resource for browsing manual annotation of the human genome project, now contains five reference genomes (human, mouse, zebrafish, pig and rat). Its introduction pages have been redesigned to enable the user to easily navigate between whole genomes and smaller multi-species haplotypic regions of interest such as the major histocompatibility complex. The VEGA browser is unique in that annotation is updated via the Human And Vertebrate Analysis aNd Annotation (HAVANA) update track every 2 weeks, allowing single gene updates to be made publicly available to the research community quickly. The user can now access different haplotypic subregions more easily, such as those from the non-obese diabetic mouse, and display them in a more intuitive way using the comparative tools. We also highlight how the user can browse manually annotated updated patches from the Genome Reference Consortium (GRC).

Marsupial immunology bounding ahead

Marsupial immune responses were previously touted as ‘primitive’ but we now know that the marsupial immune system is complex and on par with that of eutherian mammals. In this manuscript we review the field of marsupial immunology, focusing on basic anatomy, developmental immunology, immunogenetics and evolution. We concentrate on advances to our understanding of marsupial immune gene architecture, made possible by the recent sequencing of the opossum, tammar wallaby and Tasmanian devil genomes. Characterisation of immune gene sequences now paves the way for the development of immunological assays that will allow us to more accurately study health and disease in marsupials.

Low major histocompatibility complex diversity in the Tasmanian devil predates European settlement and may explain susceptibility to disease epidemics

The Tasmanian devil (Sarcophilus harrisii) is at risk of extinction owing to the emergence of a contagious cancer known as devil facial tumour disease (DFTD). The emergence and spread of DFTD has been linked to low genetic diversity in the major histocompatibility complex (MHC). We examined MHC diversity in historical and ancient devils to determine whether loss of diversity is recent or predates European settlement in Australia. Our results reveal no additional diversity in historical Tasmanian samples. Mainland devils had common modern variants plus six new variants that are highly similar to existing alleles. We conclude that low MHC diversity has been a feature of devil populations since at least the Mid-Holocene and could explain their tumultuous history of population crashes.

New insights into the role of MHC diversity in devil facial tumour disease

BACKGROUND

Devil facial tumour disease (DFTD) is a fatal contagious cancer that has decimated Tasmanian devil populations. The tumour has spread without invoking immune responses, possibly due to low levels of Major Histocompatibility Complex (MHC) diversity in Tasmanian devils. Animals from a region in north-western Tasmania have lower infection rates than those in the east of the state. This area is a genetic transition zone between sub-populations, with individuals from north-western Tasmania displaying greater diversity than eastern devils at MHC genes, primarily through MHC class I gene copy number variation. Here we test the hypothesis that animals that remain healthy and tumour free show predictable differences at MHC loci compared to animals that develop the disease.

METHODOLOGY/PRINCIPAL FINDINGS

We compared MHC class I sequences in 29 healthy and 22 diseased Tasmanian devils from West Pencil Pine, a population in north-western Tasmania exhibiting reduced disease impacts of DFTD. Amplified alleles were assigned to four loci, Saha-UA, Saha-UB, Saha-UC and Saha-UD based on recently obtained genomic sequence data. Copy number variation (caused by a deletion) at Saha-UA was confirmed using a PCR assay. No association between the frequency of this deletion and disease status was identified. All individuals had alleles at Saha-UD, disproving theories of disease susceptibility relating to copy number variation at this locus. Genetic variation between the two sub-groups (healthy and diseased) was also compared using eight MHC-linked microsatellite markers. No significant differences were identified in allele frequency, however differences were noted in the genotype frequencies of two microsatellites located near non-antigen presenting genes within the MHC.

CONCLUSIONS/SIGNIFICANCE

We did not find predictable differences in MHC class I copy number variation to account for differences in susceptibility to DFTD. Genotypic data was equivocal but indentified genomic areas for further study.