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.

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).

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.

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.

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.

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.

Improved mammalian family phylogeny using gap-rare multiple sequence alignment: A timetree of extant placentals and marsupials.

The timing of mammalian diversification in relation to the Cretaceous-Paleogene (KPg) mass extinction continues to be a subject of substantial debate. Previous studies have either focused on limited taxonomic samples with available whole-genome data or relied on short sequence alignments coupled with extensive species samples. In the present study, we improved an existing dataset from the landmark study of Meredith et al. (2011) by filling in missing fragments and further generated another dataset containing 120 taxa and 98 exonic markers. Using these two datasets, we then constructed phylogenies for extant mammalian families, providing improved resolution of many conflicting relationships. Moreover, the timetrees generated, which were calibrated using appropriate molecular clock models and multiple fossil records, indicated that the interordinal diversification of placental mammals initiated before the Late Cretaceous period. Additionally, intraordinal diversification of both extant placental and marsupial lineages accelerated after the KPg boundary, supporting the hypothesis that the availability of numerous vacant ecological niches subsequent to the mass extinction event facilitated rapid diversification. Thus, our results support a scenario of placental radiation characterized by both basal cladogenesis and active interordinal divergences spanning from the Late Cretaceous into the Paleogene.

Three Blind Moles: Molecular Evolutionary Insights on the Tempo and Mode of Convergent Eye Degeneration in Notoryctes typhlops (Southern Marsupial Mole) and Two Chrysochlorids (Golden Moles).

Golden moles (Chrysochloridae) and marsupial moles (Notoryctidae) are textbook examples of convergent evolution. Both taxa are highly adapted to subterranean lifestyles and have powerful limbs for digging through the soil/sand, ears that are adapted for low-frequency hearing, vestigial eyes that are covered by skin and fur, and the absence of optic nerve connections between the eyes and the brain. The eyes of marsupial moles also lack a lens as well as retinal rods and cones. Two hypotheses have been proposed to account for the greater degeneracy of the eyes of marsupial moles than golden moles. First, marsupial moles may have had more time to adapt to their underground habitat than other moles. Second, the eyes of marsupial moles may have been rapidly and recently vestigialized to (1) reduce the injurious effects of sand getting into the eyes and (2) accommodate the enlargement of lacrimal glands that keep the nasal cavity moist and prevent the entry of sand into the nasal passages during burrowing. Here, we employ molecular evolutionary methods on DNA sequences for 38 eye genes, most of which are eye-specific, to investigate the timing of relaxed selection (=neutral evolution) for different groups of eye-specific genes that serve as proxies for distinct functional components of the eye (rod phototransduction, cone phototransduction, lens/cornea). Our taxon sampling included 12 afrothere species, of which two are golden moles (Amblysomus hottentotus, Chrysochloris asiatica), and 28 marsupial species including two individuals of the southern marsupial mole (Notoryctes typhlops). Most of the sequences were mined from databases, but we also provide new genome data for A. hottentotus and one of the two N. typhlops individuals. Even though the eyes of golden moles are less degenerate than the eyes of marsupial moles, there are more inactivating mutations (e.g., frameshift indels, premature stop codons) in their cone phototransduction and lens/cornea genes than in orthologous genes of the marsupial mole. We estimate that cone phototransduction recovery genes were inactivated first in each group, followed by lens/cornea genes and then cone phototransduction activation genes. All three groups of genes were inactivated earlier in golden moles than in marsupial moles. For the latter, we estimate that lens/cornea genes were inactivated ~17.8 million years ago (MYA) when stem notoryctids were burrowing in the soft soils of Australian rainforests. Selection on phototransduction activation genes was relaxed much later (5.38 MYA), during the early stages of Australia's aridification that produced coastal sand plains and eventually sand dunes. Unlike cone phototransduction activation genes, rod phototransduction activation genes are intact in both golden moles and one of the two individuals of N. typhlops. A second marsupial mole individual has just a single inactivating mutation in one of the rod phototransduction activation genes (PDE6B). One explanation for this result is that some rod phototransduction activation genes are pleiotropic and are expressed in extraocular tissues, possibly in conjunction with sperm thermotaxis.

Comparative genomics of the T cell receptor µ locus in marsupials and monotremes.

T cells are a primary component of the vertebrate adaptive immune system. There are three mammalian T cell lineages based on their T cell receptors (TCR). The αß T cells and γδ T cells are ancient and found broadly in vertebrates. The more recently discovered γµ T cells are uniquely mammalian and only found in marsupials and monotremes. In this study, we compare the TCRµ locus (TRM) across the genomes of two marsupials, the gray short-tailed opossum and Tasmanian devil, and one monotreme, the platypus. These analyses revealed lineage-specific duplications, common to all non-eutherian mammals described. There is conserved synteny in the TRM loci of both marsupials but not in the monotreme. Our results are consistent with an ancestral cluster organization which was present in the last common mammalian ancestor which underwent lineage-specific duplications and divergence among the non-eutherian mammals.

Tasmanian devil cathelicidins exhibit anticancer activity against Devil Facial Tumour Disease (DFTD) cells.

The Tasmanian devil (Sarcophilus harrisii) is endangered due to the spread of Devil Facial Tumour Disease (DFTD), a contagious cancer with no current treatment options. Here we test whether seven recently characterized Tasmanian devil cathelicidins are involved in cancer regulation. We measured DFTD cell viability in vitro following incubation with each of the seven peptides and describe the effect of each on gene expression in treated cells. Four cathelicidins (Saha-CATH3, 4, 5 and 6) were toxic to DFTD cells and caused general signs of cellular stress. The most toxic peptide (Saha-CATH5) also suppressed the ERBB and YAP1/TAZ signaling pathways, both of which have been identified as important drivers of cancer proliferation. Three cathelicidins induced inflammatory pathways in DFTD cells that may potentially recruit immune cells in vivo. This study suggests that devil cathelicidins have some anti-cancer and inflammatory functions and should be explored further to determine whether they have potential as treatment leads.

Historical RNA expression profiles from the extinct Tasmanian tiger.

Paleogenomics continues to yield valuable insights into the evolution, population dynamics, and ecology of our ancestors and other extinct species. However, DNA sequencing cannot reveal tissue-specific gene expression, cellular identity, or gene regulation, which are only attainable at the transcriptional level. Pioneering studies have shown that useful RNA can be extracted from ancient specimens preserved in permafrost and historical skins from extant canids, but no attempts have been made so far on extinct species. We extract, sequence, and analyze historical RNA from muscle and skin tissue of a ∼130-year-old Tasmanian tiger (Thylacinus cynocephalus) preserved in desiccation at room temperature in a museum collection. The transcriptional profiles closely resemble those of extant species, revealing specific anatomical features such as slow muscle fibers or blood infiltration. Metatranscriptomic analysis, RNA damage, tissue-specific RNA profiles, and expression hotspots genome-wide further confirm the thylacine origin of the sequences. RNA sequences are used to improve protein-coding and noncoding annotations, evidencing missing exonic loci and the location of ribosomal RNA genes while increasing the number of annotated thylacine microRNAs from 62 to 325. We discover a thylacine-specific microRNA isoform that could not have been confirmed without RNA evidence. Finally, we detect traces of RNA viruses, suggesting the possibility of profiling viral evolution. Our results represent the first successful attempt to obtain transcriptional profiles from an extinct animal species, providing thought-to-be-lost information on gene expression dynamics. These findings hold promising implications for the study of RNA molecules across the vast collections of natural history museums and from well-preserved permafrost remains.

Patterns of ontogenetic evolution across extant marsupials reflect different allometric pathways to ecomorphological diversity.

The relatively high level of morphological diversity in Australasian marsupials compared to that observed among American marsupials remains poorly understood. We undertake a comprehensive macroevolutionary analysis of ontogenetic allometry of American and Australasian marsupials to examine whether the contrasting levels of morphological diversity in these groups are reflected in their patterns of allometric evolution. We collate ontogenetic series for 62 species and 18 families of marsupials (n = 2091 specimens), spanning across extant marsupial diversity. Our results demonstrate significant lability of ontogenetic allometric trajectories among American and Australasian marsupials, yet a phylogenetically structured pattern of allometric evolution is preserved. Here we show that species diverging more than 65 million years ago converge in their patterns of ontogenetic allometry under animalivorous and herbivorous diets, and that Australasian marsupials do not show significantly greater variation in patterns of ontogenetic allometry than their American counterparts, despite displaying greater magnitudes of extant ecomorphological diversity.

Torpor-responsive microRNAs in the heart of the Monito del monte, Dromiciops gliroides.

The marsupial Monito del monte (Dromiciops gliroides) utilizes both daily and seasonal bouts of torpor to preserve energy and prolong survival during periods of cold and unpredictable food availability. Torpor involves changes in cellular metabolism, including specific changes to gene expression that is coordinated in part, by the posttranscriptional gene silencing activity of microRNAs (miRNA). Previously, differential miRNA expression has been identified in D. gliroides liver and skeletal muscle; however, miRNAs in the heart of Monito del monte remained unstudied. In this study, the expression of 82 miRNAs was assessed in the hearts of active and torpid D. gliroides, finding that 14 were significantly differentially expressed during torpor. These 14 miRNAs were then used in bioinformatic analyses to identify Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways that were predicted to be most affected by these differentially expressed miRNAs. Overexpressed miRNAs were predicted to primarily regulate glycosaminoglycan biosynthesis, along with various signaling pathways such as Phosphoinositide-3-kinase/protein kinase B and transforming growth factor-ß. Similarly, signaling pathways including phosphatidylinositol and Hippo were predicted to be regulated by the underexpression of miRNAs during torpor. Together, these results suggest potential molecular adaptations that protect against irreversible tissue damage and enable continued cardiac and vascular function despite hypothermia and limited organ perfusion during torpor.

Pedomorphosis in the ancestry of marsupial mammals.

Within mammals, different reproductive strategies (e.g., egg laying, live birth of extremely underdeveloped young, and live birth of well-developed young) have been linked to divergent evolutionary histories. How and when developmental variation across mammals arose is unclear. While egg laying is unquestionably considered the ancestral state for all mammals, many long-standing biases treat the extreme underdeveloped state of marsupial young as the ancestral state for therian mammals (clade including both marsupials and placentals), with the well-developed young of placentals often considered the derived mode of development. Here, we quantify mammalian cranial morphological development and estimate ancestral patterns of cranial shape development using geometric morphometric analysis of the largest comparative ontogenetic dataset of mammals to date (165 specimens, 22 species). We identify a conserved region of cranial morphospace for fetal specimens, after which cranial morphology diversified through ontogeny in a cone-shaped pattern. This cone-shaped pattern of development distinctively reflected the upper half of the developmental hourglass model. Moreover, cranial morphological variation was found to be significantly associated with the level of development (position on the altricial-precocial spectrum) exhibited at birth. Estimation of ancestral state allometry (size-related shape change) reconstructs marsupials as pedomorphic relative to the ancestral therian mammal. In contrast, the estimated allometries for the ancestral placental and ancestral therian were indistinguishable. Thus, from our results, we hypothesize that placental mammal cranial development most closely reflects that of the ancestral therian mammal, while marsupial cranial development represents a more derived mode of mammalian development, in stark contrast to many interpretations of mammalian evolution.

The evolution of two transmissible cancers in Tasmanian devils.

Tasmanian devils have spawned two transmissible cancer lineages, named devil facial tumor 1 (DFT1) and devil facial tumor 2 (DFT2). We investigated the genetic diversity and evolution of these clones by analyzing 78 DFT1 and 41 DFT2 genomes relative to a newly assembled, chromosome-level reference. Time-resolved phylogenetic trees reveal that DFT1 first emerged in 1986 (1982 to 1989) and DFT2 in 2011 (2009 to 2012). Subclone analysis documents transmission of heterogeneous cell populations. DFT2 has faster mutation rates than DFT1 across all variant classes, including substitutions, indels, rearrangements, transposable element insertions, and copy number alterations, and we identify a hypermutated DFT1 lineage with defective DNA mismatch repair. Several loci show plausible evidence of positive selection in DFT1 or DFT2, including loss of chromosome Y and inactivation of MGA, but none are common to both cancers. This study reveals the parallel long-term evolution of two transmissible cancers inhabiting a common niche in Tasmanian devils.

Origin and Evolution of Marsupial-specific Imprinting Clusters Through Lineage-specific Gene Duplications and Acquisition of Promoter Differential Methylation.

Genomic imprinting is a parent-of-origin-specific expression phenomenon that plays fundamental roles in many biological processes. In animals, imprinting is only observed in therian mammals, with ∼200 imprinted genes known in humans and mice. The imprinting pattern in marsupials has been minimally investigated by examining orthologs to known eutherian imprinted genes. To identify marsupial-specific imprinting in an unbiased way, we performed RNA-seq studies on samples of fetal brain and placenta from the reciprocal cross progeny of two laboratory opossum stocks. We inferred allele-specific expression for >3,000 expressed genes and discovered/validated 13 imprinted genes, including three previously known imprinted genes, Igf2r, Peg10, and H19. We estimate that marsupials imprint ∼60 autosomal genes, which is a much smaller set compared with eutherians. Among the nine novel imprinted genes, three noncoding RNAs have no known homologs in eutherian mammals, while the remaining genes have important functions in pluripotency, transcription regulation, nucleolar homeostasis, and neural differentiation. Methylation analyses at promoter CpG islands revealed differentially methylated regions in five of these marsupial-specific imprinted genes, suggesting that differential methylation is a common mechanism in the epigenetic regulation of marsupial imprinting. Clustering and co-regulation were observed at marsupial imprinting loci Pou5f3-Npdc1 and Nkrfl-Ipncr2, but eutherian-type multi-gene imprinting clusters were not detected. Also differing from eutherian mammals, the brain and placenta imprinting profiles are remarkably similar in opossums, presumably due to the shared origin of these organs from the trophectoderm. Our results contribute to a fuller understanding of the origin, evolution, and mechanisms of genomic imprinting in therian mammals.