Generation of Devil Facial Tumour Cells Co-Expressing MHC With CD80, CD86 or 41BBL to Enhance Tumour Immunogenicity.

The major histocompatibility complex (MHC) molecules play an integral role in the adaptive immune response to transmissible cancers through tumour antigen presentation and recognition of allogeneic MHC molecules. The transmissible devil facial tumours 1 and 2 (DFT1 and DFT2) modulate MHC-I antigen presentation to evade host immune responses and facilitate transmission of tumours cells to new Tasmanian devil (Sarcophilus harrisii) hosts. To enhance T-cell-driven tumour immunogenicity for vaccination and immunotherapy, DFT1 and DFT2 cells were co-transfected with (i) NLRC5 for MHC-I expression or CIITA for MHC-I and MHC-II expression, and (ii) a co-stimulatory molecule, either CD80, CD86 or 41BBL. The co-transfected DFT cells presented enhanced expression of MHC-I and/or MHC-II. As few devil-specific monoclonal antibodies exist, we used recombinant CTLA4 and 41BB fused to a fluorescent protein to confirm expression of cell surface CD80, CD86 and 41BBL. The capacity for these cells to induce T-cell responses including PD1 and IFNG expression was evaluated in in vitro co-culture assays with captive devil peripheral blood mononuclear cells (PBMCs). Although PBMC viability had increased, there was no evidence of enhanced T-cell activation. This system can be used to identify additional factors required to promote activation of naïve devil T-cells in vitro.

PRKACB is a novel imprinted gene in marsupials.

BACKGROUND: Genomic imprinting results in parent-of-origin-specific gene expression and, among vertebrates, is found only in therian mammals: marsupials and eutherians. A differentially methylated region (DMR), in which the methylation status of CpG dinucleotides differs between the two alleles, can mark the parental identity of imprinted genes. We developed a computational pipeline that detected CpG islands (CGIs) marked by both methylated and unmethylated signals in whole genome bisulfite sequencing data. This approach identified candidate marsupial DMRs in a publicly available koala methylome. One of these candidate DMRs was associated with PRKACB, a gene encoding the protein kinase A catalytic subunit beta. Nothing is known about the imprinting status of PRKACB in eutherian mammals although mutations of this gene are associated with endocrine neoplasia and other developmental disorders. RESULTS: In the tammar wallaby and brushtail possum there was parent-of-origin-specific DNA methylation in the PRKACB DMR in which the maternal allele was methylated and the paternal allele was unmethylated. There were multiple RNAs transcribed from this locus. Allele-specific expression analysis identified paternal expression of a PRKACB lncRNA and an mRNA isoform. Comparison of the PRKACB gene start site between marsupials and eutherians demonstrated that the CGI is longer in marsupials. The PRKACB gene product functions in the same signalling pathway as the guanine nucleotide-binding protein alpha subunit encoded at the GNAS locus, a known eutherian imprinted gene. In a mouse methylome Gnas had three differentially methylated CGIs, while in the koala methylome the GNAS locus had two unmethylated CGIs. CONCLUSIONS: We conclude that PRKACB is a novel, DMR-associated marsupial imprinted gene. Imprinting of PRKACB in marsupials and GNAS in eutherians may indicate a conserved selection pressure for imprinting of the protein kinase A signalling pathway in therians with the two lineages adapting by imprinting different genes.

Adaptive potential in the face of a transmissible cancer in Tasmanian devils.

Emerging infectious diseases (EIDs) not only cause catastrophic declines in wildlife populations but also generate selective pressures that may result in rapid evolutionary responses. One such EID is devil facial tumour disease (DFTD) in the Tasmanian devil. DFTD is almost always fatal and has reduced the average lifespan of individuals by around 2 years, likely causing strong selection for traits that reduce susceptibility to the disease, but population decline has also left Tasmanian devils vulnerable to inbreeding depression. We analysed 22 years of data from an ongoing study of a population of Tasmanian devils on Freycinet Peninsula, Tasmania, to (1) identify whether DFTD may be causing selection on body size, by estimating phenotypic and genetic correlations between DFTD and size traits, (2) estimate the additive genetic variance of susceptibility to DFTD, and (3) investigate whether size traits or susceptibility to DFTD were under inbreeding depression. We found a positive phenotypic relationship between head width and susceptibility to DFTD, but this was not underpinned by a genetic correlation. Conversely, we found a negative phenotypic relationship between body weight and susceptibility to DFTD, and there was evidence for a negative genetic correlation between susceptibility to DFTD and body weight. There was additive genetic variance in susceptibility to DFTD, head width and body weight, but there was no evidence for inbreeding depression in any of these traits. These results suggest that Tasmanian devils have the potential to respond adaptively to DFTD, although the realised evolutionary response will critically further depend on the evolution of DFTD itself.

Imprinted X chromosome inactivation in marsupials: The paternal X arrives at the egg with a silent DNA methylation profile.

X chromosome inactivation (XCI) is an epigenetic process that results in the transcriptional silencing of one X chromosome in the somatic cells of females. This phenomenon is common to both eutherian and marsupial mammals, but there are fundamental differences. In eutherians, the X chosen for silencing is random. DNA methylation on the eutherian inactive X is high at transcription start sites (TSSs) and their flanking regions, resulting in universally high DNA methylation. This contrasts XCI in marsupials where the paternally derived X is always silenced, and in which DNA methylation is low at TSSs and flanking regions. Here, we examined the DNA methylation status of the tammar wallaby X chromosome during spermatogenesis to determine the DNA methylation profile of the paternal X prior to and at fertilization. Whole genome enzymatic methylation sequencing was carried out on enriched flow-sorted populations of premeiotic, meiotic, and postmeiotic cells. We observed that the X displayed a pattern of DNA methylation from spermatogonia to mature sperm that reflected the inactive X in female somatic tissue. Therefore, the paternal X chromosome arrives at the egg with a DNA methylation profile that reflects the transcriptionally silent X in adult female somatic tissue. We present this epigenetic signature as a candidate for the long sought-after imprint for paternal XCI in marsupials.

Conversion of Mouse-Derived Hybridomas to Tasmanian Devil Recombinant IgG Antibodies.

The hybridoma method for production of monoclonal antibodies has been a cornerstone of biomedical research for several decades. Here we convert the monoclonal antibody sequence from mouse-derived hybridomas into a "devilized" recombinant antibody with devil IgG heavy chain and IgK light chain. The chimeric recombinant antibody can be used in functional assays, immunotherapy, and to improve understanding of antibodies and Fc receptors in Tasmanian devils. The process can be readily modified for other species.

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.

Two-stage evolution of mammalian adipose tissue thermogenesis.

Brown adipose tissue (BAT) is a heater organ that expresses thermogenic uncoupling protein 1 (UCP1) to maintain high body temperatures during cold stress. BAT thermogenesis is considered an overarching mammalian trait, but its evolutionary origin is unknown. We show that adipose tissue of marsupials, which diverged from eutherian mammals ~150 million years ago, expresses a nonthermogenic UCP1 variant governed by a partial transcriptomic BAT signature similar to that found in eutherian beige adipose tissue. We found that the reconstructed UCP1 sequence of the common eutherian ancestor displayed typical thermogenic activity, whereas therian ancestor UCP1 is nonthermogenic. Thus, mammalian adipose tissue thermogenesis may have evolved in two distinct stages, with a prethermogenic stage in the common therian ancestor linking UCP1 expression to adipose tissue and thermal stress. We propose that in a second stage, UCP1 acquired its thermogenic function specifically in eutherians, such that the onset of mammalian BAT thermogenesis occurred only after the divergence from marsupials.

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

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

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

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

Characterisation of defensins across the marsupial family tree.

Defensins are antimicrobial peptides involved in innate immunity, and gene number differs amongst eutherian mammals. Few studies have investigated defensins in marsupials, despite their potential involvement in immunological protection of altricial young. Here we use recently sequenced marsupial genomes and transcriptomes to annotate defensins in nine species across the marsupial family tree. We characterised 35 alpha and 286 beta defensins; gene number differed between species, although Dasyuromorphs had the largest repertoire. Defensins were encoded in three gene clusters within the genome, syntenic to eutherians, and were expressed in the pouch and mammary gland. Marsupial beta defensins were closely related to eutherians, however marsupial alpha defensins were more divergent. We identified marsupial orthologs of human DEFB3 and 6, and several marsupial-specific beta defensin lineages which may have novel functions. Marsupial predicted mature peptides were highly variable in length and sequence composition. We propose candidate peptides for future testing to elucidate the function of marsupial defensins.

Emx2 underlies the development and evolution of marsupial gliding membranes.

Phenotypic variation among species is a product of evolutionary changes to developmental programs1,2. However, how these changes generate novel morphological traits remains largely unclear. Here we studied the genomic and developmental basis of the mammalian gliding membrane, or patagium-an adaptative trait that has repeatedly evolved in different lineages, including in closely related marsupial species. Through comparative genomic analysis of 15 marsupial genomes, both from gliding and non-gliding species, we find that the Emx2 locus experienced lineage-specific patterns of accelerated cis-regulatory evolution in gliding species. By combining epigenomics, transcriptomics and in-pouch marsupial transgenics, we show that Emx2 is a critical upstream regulator of patagium development. Moreover, we identify different cis-regulatory elements that may be responsible for driving increased Emx2 expression levels in gliding species. Lastly, using mouse functional experiments, we find evidence that Emx2 expression patterns in gliders may have been modified from a pre-existing program found in all mammals. Together, our results suggest that patagia repeatedly originated through a process of convergent genomic evolution, whereby regulation of Emx2 was altered by distinct cis-regulatory elements in independently evolved species. Thus, different regulatory elements targeting the same key developmental gene may constitute an effective strategy by which natural selection has harnessed regulatory evolution in marsupial genomes to generate phenotypic novelty.

Transmissible cancers, the genomes that do not melt down.

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

Complex associations between cancer progression and immune gene expression reveals early influence of transmissible cancer on Tasmanian devils.

The world's largest extant carnivorous marsupial, the Tasmanian devil, is challenged by Devil Facial Tumor Disease (DFTD), a fatal, clonally transmitted cancer. In two decades, DFTD has spread across 95% of the species distributional range. A previous study has shown that factors such as season, geographic location, and infection with DFTD can impact the expression of immune genes in Tasmanian devils. To date, no study has investigated within-individual immune gene expression changes prior to and throughout the course of DFTD infection. To explore possible changes in immune response, we investigated four locations across Tasmania that differed in DFTD exposure history, ranging between 2 and >30 years. Our study demonstrated considerable complexity in the immune responses to DFTD. The same factors (sex, age, season, location and DFTD infection) affected immune gene expression both across and within devils, although seasonal and location specific variations were diminished in DFTD affected devils. We also found that expression of both adaptive and innate immune genes starts to alter early in DFTD infection and continues to change as DFTD progresses. A novel finding was that the lower expression of immune genes MHC-II, NKG2D and CD8 may predict susceptibility to earlier DFTD infection. A case study of a single devil with regressed tumor showed opposite/contrasting immune gene expression patterns compared to the general trends observed across devils with DFTD infection. Our study highlights the complexity of DFTD's interactions with the host immune system and the need for long-term studies to fully understand how DFTD alters the evolutionary trajectory of devil immunity.

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.

A family of olfactory receptors uniquely expanded in marsupial and monotreme genomes are expressed by a T cell subset also unique to marsupials and monotremes.

Olfactory receptors (OR), expressed on olfactory neurons, mediate the sense of smell. Recently, OR have also been shown to be expressed in non-olfactory tissues, including cells of the immune system. An analysis of single-cell transcriptomes of splenocytes of the grey short-tailed opossum (Monodelphis domestica) found OR are expressed on a subset of T cells, the γµ T cells, that are unique to marsupials and monotremes. A majority of opossum γµ T cells transcriptomes contain OR family 14 transcripts, specifically, from the OR14C subfamily. Amongst the mammals, the OR14 gene family is expanded in the genomes of marsupials and monotremes, and rarer or absent in placental mammals. In summary, here we demonstrate the intriguing correlation that a family of OR genes, abundant in the genomes of marsupials and monotremes, are ectopically expressed in a particular subset of T cells unique to the marsupials and monotremes.

A recent gibbon ape leukemia virus germline integration in a rodent from New Guinea.

Germline colonization by retroviruses results in the formation of endogenous retroviruses (ERVs). Most colonization's occurred millions of years ago. However, in the Australo-Papuan region (Australia and New Guinea), several recent germline colonization events have been discovered. The Wallace Line separates much of Southeast Asia from the Australo-Papuan region restricting faunal and pathogen dispersion. West of the Wallace Line, gibbon ape leukemia viruses (GALVs) have been isolated from captive gibbons. Two microbat species from China appear to have been infected naturally. East of Wallace's Line, the woolly monkey virus (a GALV) and the closely related koala retrovirus (KoRV) have been detected in eutherians and marsupials in the Australo-Papuan region, often vertically transmitted. The detected vertically transmitted GALV-like viruses in Australo-Papuan fauna compared to sporadic horizontal transmission in Southeast Asia and China suggest the GALV-KoRV clade originates in the former region and further models of early-stage genome colonization may be found. We screened 278 samples, seven bat and one rodent family endemic to the Australo-Papuan region and bat and rodent species found on both sides of the Wallace Line. We identified two rodents (Melomys) from Australia and Papua New Guinea and no bat species harboring GALV-like retroviruses. Melomys leucogaster from New Guinea harbored a genomically complete replication-competent retrovirus with a shared integration site among individuals. The integration was only present in some individuals of the species indicating this retrovirus is at the earliest stages of germline colonization of the Melomys genome, providing a new small wild mammal model of early-stage genome colonization.

Disease-driven top predator decline affects mesopredator population genomic structure.

Top predator declines are pervasive and often have dramatic effects on ecological communities via changes in food web dynamics, but their evolutionary consequences are virtually unknown. Tasmania's top terrestrial predator, the Tasmanian devil, is declining due to a lethal transmissible cancer. Spotted-tailed quolls benefit via mesopredator release, and they alter their behaviour and resource use concomitant with devil declines and increased disease duration. Here, using a landscape community genomics framework to identify environmental drivers of population genomic structure and signatures of selection, we show that these biotic factors are consistently among the top variables explaining genomic structure of the quoll. Landscape resistance negatively correlates with devil density, suggesting that devil declines will increase quoll genetic subdivision over time, despite no change in quoll densities detected by camera trap studies. Devil density also contributes to signatures of selection in the quoll genome, including genes associated with muscle development and locomotion. Our results provide some of the first evidence of the evolutionary impacts of competition between a top predator and a mesopredator species in the context of a trophic cascade. As top predator declines are increasing globally, our framework can serve as a model for future studies of evolutionary impacts of altered ecological interactions.

Expanding the evo-devo toolkit: generation of 3D mammary tissue from diverse mammals.

The varying pathways of mammary gland development across species and evolutionary history are underexplored, largely due to a lack of model systems. Recent progress in organoid technology holds the promise of enabling in-depth studies of the developmental adaptations that have occurred throughout the evolution of different species, fostering beneficial phenotypes. The practical application of this technology for mammary glands has been mostly confined to rodents and humans. In the current study, we have successfully created next-generation 3D mammary gland organoids from eight eutherian mammals and the first branched organoid of a marsupial mammary gland. Using mammary organoids, we identified a role for ROCK protein in regulating branching morphogenesis, a role that manifests differently in organoids from different mammals. This finding demonstrates the utility of the 3D organoid model for understanding the evolution and adaptations of signaling pathways. These achievements highlight the potential for organoid models to expand our understanding of mammary gland biology and evolution, and their potential utility in studies of lactation or breast cancer.