Hidden diversity in an ecologically specialized genus of Australian marsupials, the feather-tailed gliders, Acrobates (Diprotodontia, Acrobatidae)

The marsupial Family Acrobatidae includes the smallest gliding marsupial species in the monotypic Acrobates, found only in eastern Australia, and an equally small non-gliding species in another monotypic genus Distoechurus, found only in New Guinea. We applied molecular genetic analysis to Acrobates to assess the systematic significance of variation in superficial external characters of the tail and hindfoot (pes). Deep divergence in mitochondrial and nuclear genes demonstrated the broad sympatry of two species consistent with prior morphological diagnoses. Morphological assessment of museum vouchers showed that their distributions overlap extensively in New South Wales and Victoria and include locations where a range of biological research was conducted on the assumption of the presence of a single species. Many of these studies cannot be reassessed because neither vouchers nor tissue suitable for molecular genetic identification were collected. Intriguingly, acrobatids are the only marsupial group with internal ear discs, and the two species of Acrobates show demonstrable differences in the morphology of this structure, the biological significance of which needs to be established. Both species of Acrobates occur widely in the eucalypt forests of south-eastern Australia, which appear to be subject to a growing threat from bushfires likely aggravated by anthropogenic climate change.

Unearthing the secrets of Australia's most enigmatic and cryptic mammal, the marsupial mole

The marsupial moles are arguably Australia's most enigmatic marsupials. Almost indistinguishable from placental (eutherian) moles, they provide a striking example of convergent evolution. Exploring the genome of the southern marsupial mole, we provide insights into its unusual biology. We show definitively by retrophylogenomic analysis that marsupial moles are most closely related to bandicoots and bilbies (order Peramelemorphia). We find evidence of a marked decline in marsupial mole effective population size, most likely preceding the arrival of humans in regions near its range, and potentially corresponding to periods of climatic change. Our analysis of loss of eye function-an adaptation to subterranean life-reveals a structured order of loss of gene function associated first with the lens, then cone, and finally rod cells. Last, we identify genetic changes suggestive of adaptation to an oxygen-poor environment and of its evolution of partially descended testes.

Ancient Hybridisation Fuelled Diversification in Acropora Corals

Introgression is the infiltration or flow of genes from one species to another through hybridisation followed by backcrossing. This may lead to incorrect phylogenetic reconstruction or divergence-time estimation. Acropora is a dominant genus of reef-building corals; however, whether this group has an introgression history before their diversification remains unclear, and previous divergence-time estimates of Acropora have not considered the impact of introgression. In this study, we broke through the limitation of a few genes and a few species and proved the existence of ancient introgression in the evolution of Acropora from whole-genome protein-coding sequences. We inferred 21.9% of all triplet loci (homologous loci from three different species) with a history of introgression and a series of introgression events with a genetic material contribution of up to 30.9% before diversification. Furthermore, 7756 nuclear loci were clustered into three groups using a multidimensional scaling algorithm, the heterogeneity of which resulted in different phylogenetic relationships. The diversification time of Acropora was estimated to be middle to late Miocene when we retained only the gene group with the lowest degree of introgression. The collision of Australia with the Pacific arcs and the Southeast Asian margin in the early Miocene, and a series of cooling events in the middle to late Miocene, may provide geographical and climatic conditions for the diversification of Acropora, respectively. Therefore, our results indicate that at the genome-wide level, ancient introgressive hybridisation may have promoted the radiation evolution of Acropora. Based on our results, the influence of introgression should be taken into account when reconstructing phylogenetic relationships and evaluating divergence time.

Spermatogonial stem cell technologies: applications from human medicine to wildlife conservation†

Spermatogonial stem cell (SSC) technologies that are currently under clinical development to reverse human infertility hold the potential to be adapted and applied for the conservation of endangered and vulnerable wildlife species. The biobanking of testis tissue containing SSCs from wildlife species, aligned with that occurring in pediatric human patients, could facilitate strategies to improve the genetic diversity and fitness of endangered populations. Approaches to utilize these SSCs could include spermatogonial transplantation or testis tissue grafting into a donor animal of the same or a closely related species, or in vitro spermatogenesis paired with assisted reproduction approaches. The primary roadblock to progress in this field is a lack of fundamental knowledge of SSC biology in non-model species. Herein, we review the current understanding of molecular mechanisms controlling SSC function in laboratory rodents and humans, and given our particular interest in the conservation of Australian marsupials, use a subset of these species as a case-study to demonstrate gaps-in-knowledge that are common to wildlife. Additionally, we review progress in the development and application of SSC technologies in fertility clinics and consider the translation potential of these techniques for species conservation pipelines.

The mammary hair of Monodelphis domestica and homology of the mammary pilosebacous unit

The unitary mammary gland is a synapomorphy of therian mammals and is thought to have evolved from the pilosebaceous organ in the mammalian stem lineage from which the lactogenic patch of monotremes is also derived. One of the key lines of evidence for the homology of the nipple and the lactogenic patch is that marsupials have retained a transient hair associated with developing mammary glands. However, these structures have not been documented since the early 20th-century drawings of Ernst Bresslau. In this study, we examine the developing mammary organs of Monodelphis domestica and document the presence of mammary hairs in 12-week-old females, as well as their absence after 18 weeks of age. Histochemical staining for cystine confirms the structures as keratinized hairs. Milk ducts of both juvenile and adult nipples show a division between KRT18+ luminal epithelium and KRT14+ ACTA2+ myoepithelium. These patterns match those in eutherians and suggest a conserved ductal morphology and mechanism of milk expulsion. Finally, PTHLH, a peptide hormone which promotes homeotic transformation of hairy skin into hairless nipples in the mouse, was detected in the Monodelphis milk duct during the mammary hair stage, suggesting that the mutual exclusivity of "hairless nipple" and "hair" organ identity is derived in eutherian mammals. These results reveal shared characteristics of the M. domestica nipple with both the eutherian nipple and the pilosebaceous organ, consistent with the evolutionary derivation of the mammary gland from an ancestral hair organ via developmental individualization of pilosebaceous and mammary identities.

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.

Evolution, divergence, and convergence in the mandibles of opossums (Didelphidae, Didelphimorphia)

Didelphid marsupials are considered a morphologically unspecialized group with a generalist diet that includes vertebrates, invertebrates, and plant matter. While cranium and scapula variation has already been examined within Didelphidae, variation in mandible shape, usually associated with diet or phylogeny in other mammalian groups, has not yet been properly assessed in the family. We evaluated the variation in mandible shape and size of didelphids (2470 specimens belonging to 94 species) using 2D geometric morphometrics. We classified the diet of the didelphids into four broad categories to assess whether morphospace ordination relates to dietary habits. We also provided the most comprehensive phylogeny for the family (123 out of the 126 living species) using 10 nuclear and mitochondrial genes. We then mapped mandible size and shape onto that phylogeny for 93 selected taxa and ancestral size and shapes were reconstructed by parsimony. We found phylogenetically structured variation in mandible morphology between didelphid groups, and our results indicate that they have a significant phylogenetic signal. The main axis of shape variation is poorly related to size, but the second is strongly allometric, indicating that allometry is not the main factor in shaping morphological diversity on their mandibles. Our results indicate that the shape and size of the ancestral mandible of didelphids would be similar to that of the current species of the genus Marmosa.

Vicariant speciation resulting from biogeographic barriers in the Australian tropics: The case of the red-cheeked dunnart (Sminthopsis virginiae)

Global biodiversity loss continues unabated, and in Australia, the rate of recent mammal extinctions is among the worst in the world. Meanwhile, the diversity among and within many endemic mammal species remains undescribed. This information is crucial to delineate species boundaries and thus inform decision-making for conservation. Sminthopsis virginiae (the red-cheeked dunnart) is a small, dasyurid marsupial found in four disjunct populations around the northern coast of Australia and New Guinea. There are three currently recognized subspecies, each occupying a distinct geographic location. Sminthopsis v. virginiae occurs in Queensland, S. v. rufigenis is distributed across New Guinea and the Aru Islands, and S. v. nitela has populations in the Top End of the Northern Territory and the Kimberley region of Western Australia. Previous molecular work has suggested the current subspecies definitions are not aligned with DNA sequence data, though the sampling was limited. We undertook a comprehensive genetic and morphological review of S. virginiae to clarify relationships within the species. This included mitochondrial (CR, 12S, and cytb) and nuclear (omega-globin, IRBP, and bfib7) loci, and morphometric analysis of skulls and whole wet-preserved specimens held in museums. Maximum Likelihood and Bayesian phylogenetic analyses resolved samples into two distinct clades, demarcated by the Gulf of Carpentaria in Australia's north. Sminthopsis. v. nitela was consistently separated from S. v. virginiae and S. v. rufigenis, based on the overall body and skull size and craniodental features, while S. v. virginiae and S. v. rufigenis were more difficult to distinguish from each other. Thus, we redescribed S. virginiae, recognizing two species, S. nitela (raised from subspecies) and S. virginiae (now comprising the subspecies S. v. virginiae and S. v. rufigenis). This study highlights the importance of recognizing cryptic mammal fauna to help address the gap in our knowledge about diagnosing diversity during a time of conservation crisis.

Populating a Continent: Phylogenomics Reveal the Timing of Australian Frog Diversification

The Australian continent's size and isolation make it an ideal place for studying the accumulation and evolution of biodiversity. Long separated from the ancient supercontinent Gondwana, most of Australia's plants and animals are unique and endemic, including the continent's frogs. Australian frogs comprise a remarkable ecological and morphological diversity categorized into a small number of distantly related radiations. We present a phylogenomic hypothesis based on an exon-capture dataset that spans the main clades of Australian myobatrachoid, pelodryadid hyloid, and microhylid frogs. Our time-calibrated phylogenomic-scale phylogeny identifies great disparity in the relative ages of these groups that vary from Gondwanan relics to recent immigrants from Asia and include arguably the continent's oldest living vertebrate radiation. This age stratification provides insight into the colonization of, and diversification on, the Australian continent through deep time, during periods of dramatic climatic and community changes. Contemporary Australian frog diversity highlights the adaptive capacity of anurans, particularly in response to heat and aridity, and explains why they are one of the continent's most visible faunas. [Anuran; adaptive radiation; Gondwana; phylogenetics].

Speciation across biomes: Rapid diversification with reproductive isolation in the Australian delicate mice

Phylogeographic studies of continental clades, especially when combined with palaeoclimate modelling, provide powerful insight into how environment drives speciation across climatic contexts. Australia, a continent characterized by disparate modern biomes and dynamic climate change, provides diverse opportunity to reconstruct the impact of past and present environments on diversification. Here, we use genomic-scale data (1310 exons and whole mitogenomes from 111 samples) to investigate Pleistocene diversification, cryptic diversity, and secondary contact in the Australian delicate mice (Hydromyini: Pseudomys), a recent radiation spanning almost all Australian environments. Across northern Australia, we find no evidence for introgression between cryptic lineages within Pseudomys delicatulus sensu lato, with palaeoclimate models supporting contraction and expansion of suitable habitat since the last glacial maximum. Despite multiple contact zones, we also find little evidence of introgression at a continental scale, with the exception of a potential hybrid zone in the mesic biome. In the arid zone, combined insights from genetic data and palaeomodels support a recent expansion in the arid specialist P. hermannsburgensis and contraction in the semi-arid P. bolami. In the face of repeated secondary contact, differences in sperm morphology and chromosomal rearrangements are potential mechanisms that maintain species boundaries in these recently diverged species. Additionally, we describe the western delicate mouse as a new species and recommend taxonomic reinstatement of the eastern delicate mouse. Overall, we show that speciation in an evolutionarily young and widespread clade has been driven by environmental change, and potentially maintained by divergence in reproductive morphology and chromosome rearrangements.

Facing the facts: adaptive trade-offs along body size ranges determine mammalian craniofacial scaling

The mammalian cranium (skull without lower jaw) is representative of mammalian diversity and is thus of particular interest to mammalian biologists across disciplines. One widely retrieved pattern accompanying mammalian cranial diversification is referred to as 'craniofacial evolutionary allometry' (CREA). This posits that adults of larger species, in a group of closely related mammals, tend to have relatively longer faces and smaller braincases. However, no process has been officially suggested to explain this pattern, there are many apparent exceptions, and its predictions potentially conflict with well-established biomechanical principles. Understanding the mechanisms behind CREA and causes for deviations from the pattern therefore has tremendous potential to explain allometry and diversification of the mammalian cranium. Here, we propose an amended framework to characterise the CREA pattern more clearly, in that 'longer faces' can arise through several kinds of evolutionary change, including elongation of the rostrum, retraction of the jaw muscles, or a more narrow or shallow skull, which all result in a generalised gracilisation of the facial skeleton with increased size. We define a standardised workflow to test for the presence of the pattern, using allometric shape predictions derived from geometric morphometrics analysis, and apply this to 22 mammalian families including marsupials, rabbits, rodents, bats, carnivores, antelopes, and whales. Our results show that increasing facial gracility with size is common, but not necessarily as ubiquitous as previously suggested. To address the mechanistic basis for this variation, we then review cranial adaptations for harder biting. These dictate that a more gracile cranium in larger species must represent a structural sacrifice in the ability to produce or withstand harder bites, relative to size. This leads us to propose that facial gracilisation in larger species is often a product of bite force allometry and phylogenetic niche conservatism, where more closely related species tend to exhibit more similar feeding ecology and biting behaviours and, therefore, absolute (size-independent) bite force requirements. Since larger species can produce the same absolute bite forces as smaller species with less effort, we propose that relaxed bite force demands can permit facial gracility in response to bone optimisation and alternative selection pressures. Thus, mammalian facial scaling represents an adaptive by-product of the shifting importance of selective pressures occurring with increased size. A reverse pattern of facial 'shortening' can accordingly also be found, and is retrieved in several cases here, where larger species incorporate novel feeding behaviours involving greater bite forces. We discuss multiple exceptions to a bite force-mediated influence on facial proportions across mammals which lead us to argue that ecomorphological specialisation of the cranium is likely to be the primary driver of facial scaling patterns, with some developmental constraints as possible secondary factors. A potential for larger species to have a wider range of cranial functions when less constrained by bite force demands might also explain why selection for larger sizes seems to be prevalent in some mammalian clades. The interplay between adaptation and constraint across size ranges thus presents an interesting consideration for a mechanistically grounded investigation of mammalian cranial allometry.

The Distinctive Biology and Characteristics of the Bare-Nosed Wombat (Vombatus ursinus)

The bare-nosed wombat is an iconic Australian fauna with remarkable biological characteristics and mythology. This solitary, muscular, fossorial, herbivorous marsupial from southeast Australia has continent and continental island subspeciation. Vombatiformes also contains hairy-nosed wombats (Lasiorhinus spp.); koala (Phascolarctos cinereus); and extinct megafauna, Phascolonus gigas (giant wombat), Diprotodon, and Thylacoleo (marsupial lion). Culturally important to Aboriginal people, bare-nosed wombats engineer ecosystems through digging, grazing, and defecation. Olfaction and cubic fecal aggregations appear critical for communication, including identity, courtship, and mating. Though among the largest fossorial herbivores, they have a nutrient-poor diet, a home range up to an order of magnitude smaller than expected, and a metabolism among the lowest extreme for mammals >10 kg. Metabolic depression may confer advantages over resource competitors and fossorial lifestyle protection from predators, fires, and climatic extremes. Bare-nosed wombats are loved and persecuted by European colonists. Recent population increases may reflect softening attitudes toward, and greater protections of, bare-nosed wombats.

Anaerobic gut fungal communities in marsupial hosts

The anaerobic gut fungi (AGF) inhabit the alimentary tracts of herbivores. In contrast to placental mammals, information regarding the identity, diversity, and community structure of AGF in marsupials is extremely sparse. Here, we characterized AGF communities in 61 fecal samples from 10 marsupial species belonging to four families in the order Diprotodontia: Vombatidae (wombats), Phascolarctidae (koalas), Phalangeridae (possums), and Macropodidae (kangaroos, wallabies, and pademelons). An amplicon-based diversity survey using the D2 region of the large ribosomal subunit as a phylogenetic marker indicated that marsupial AGF communities were dominated by eight genera commonly encountered in placental herbivores (Neocallimastix, Caecomyces, Cyllamyces, Anaeromyces, Orpinomyces, Piromyces, Pecoramyces, and Khoyollomyces). Community structure analysis revealed a high level of stochasticity, and ordination approaches did not reveal a significant role for the animal host, gut type, dietary preferences, or lifestyle in structuring marsupial AGF communities. Marsupial foregut and hindgut communities displayed diversity and community structure patterns comparable to AGF communities typically encountered in placental foregut hosts while exhibiting a higher level of diversity and a distinct community structure compared to placental hindgut communities. Quantification of AGF load using quantitative PCR indicated a significantly smaller load in marsupial hosts compared to their placental counterparts. Isolation efforts were only successful from a single red kangaroo fecal sample and yielded a Khoyollomyces ramosus isolate closely related to strains previously isolated from placental hosts. Our results suggest that AGF communities in marsupials are in low abundance and show little signs of selection based on ecological and evolutionary factors.IMPORTANCEThe AGF are integral part of the microbiome of herbivores. They play a crucial role in breaking down plant biomass in hindgut and foregut fermenters. The majority of research has been conducted on the AGF community in placental mammalian hosts. However, it is important to note that many marsupial mammals are also herbivores and employ a hindgut or foregut fermentation strategy for breaking down plant biomass. So far, very little is known regarding the AGF diversity and community structure in marsupial mammals. To fill this knowledge gap, we conducted an amplicon-based diversity survey targeting AGF in 61 fecal samples from 10 marsupial species. We hypothesize that, given the distinct evolutionary history and alimentary tract architecture, novel and unique AGF communities would be encountered in marsupials. Our results indicate that marsupial AGF communities are highly stochastic, present in relatively low loads, and display community structure patterns comparable to AGF communities typically encountered in placental foregut hosts. Our results indicate that marsupial hosts harbor AGF communities; however, in contrast to the strong pattern of phylosymbiosis typically observed between AGF and placental herbivores, the identity and gut architecture appear to play a minor role in structuring AGF communities in marsupials.

A supermatrix phylogeny of the world's bees (Hymenoptera: Anthophila)

The increasing availability of large molecular phylogenies has provided new opportunities to study the evolution of species traits, their origins and diversification, and biogeography; yet there are limited attempts to synthesise existing phylogenetic information for major insect groups. Bees (Hymenoptera: Anthophila) are a large group of insect pollinators that have a worldwide distribution, and a wide variation in ecology, morphology, and life-history traits, including sociality. For these reasons, as well as their major economic importance as pollinators, numerous molecular phylogenetic studies of family and genus-level relationships have been published, providing an opportunity to assemble a bee 'tree-of-life'. We used publicly available genetic sequence data, including phylogenomic data, reconciled to a taxonomic database, to produce a concatenated supermatrix phylogeny for the Anthophila comprising 4,586 bee species, representing 23% of species and 82% of genera. At family, subfamily, and tribe levels, support for expected relationships was robust, but between and within some genera relationships remain uncertain. Within families, sampling of genera ranged from 67 to 100% but species coverage was lower (17-41%). Our phylogeny mostly reproduces the relationships found in recent phylogenomic studies with a few exceptions. We provide a summary of these differences and the current state of molecular data available and its gaps. We discuss the advantages and limitations of this bee supermatrix phylogeny (available online at beetreeoflife.org), which may enable new insights into long standing questions about evolutionary drivers in bees, and potentially insects more generally.

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.

A probable koala from the Oligocene of central Australia provides insights into early diprotodontian evolution

Diprotodontians are the morphologically and ecologically most diverse order of marsupials. However, an approximately 30-million-year gap in the Australian terrestrial vertebrate fossil record means that the first half of diprotodontian evolution is unknown. Fossil taxa from immediately either side of this gap are therefore critical for reconstructing the early evolution of the order. Here we report the likely oldest-known koala relatives (Phascolarctidae), from the late Oligocene Pwerte Marnte Marnte Local Fauna (central Australia). These include coeval species of Madakoala and Nimiokoala, as well as a new probable koala (?Phascolarctidae). The new taxon, Lumakoala blackae gen. et sp. nov., was comparable in size to the smallest-known phascolarctids, with body-mass estimates of 2.2-2.6 kg. Its bunoselenodont upper molars retain the primitive metatherian condition of a continuous centrocrista, and distinct stylar cusps B and D which lacked occlusion with the hypoconid. This structural arrangement: (1) suggests a morphocline within Phascolarctidae from bunoselenodonty to selenodonty; and (2) better clarifies the evolutionary transitions between molar morphologies within Vombatomorphia. We hypothesize that the molar form of Lumakoala blackae approximates the ancestral condition of the suborder Vombatiformes. Furthermore, it provides a plausible link between diprotodontians and the putative polydolopimorphians Chulpasia jimthorselli and Thylacotinga bartholomaii from the early Eocene Tingamarra Local Fauna (eastern Australia), which we infer as having molar morphologies consistent with stem diprotodontians.

Hiding in plain sight: two new species of diminutive marsupial (Dasyuridae: Planigale) from the Pilbara, Australia

Many of Australias smaller marsupial species have been taxonomically described in just the past 50 years, and the Dasyuridae, a speciose family of carnivores, is known to harbour many cryptic taxa. Evidence from molecular studies is being increasingly utilised to help revise species boundaries and focus taxonomic efforts, and research over the past two decades has identified several undescribed genetic lineages within the dasyurid genus Planigale. Here, we describe two new species, Planigale kendricki sp. nov. (formerly known as Planigale 1) and P. tealei sp. nov. (formerly known as Planigale sp. Mt Tom Price). The two new species have broadly overlapping distributions in the Pilbara region of Western Australia. The new species are genetically distinct from each other and from all other members of the genus, at both mitochondrial and nuclear loci, and morphologically, in both external and craniodental characters. The new species are found in regional sympatry within the Pilbara but occupy different habitat types at local scales. This work makes a start at resolving the cryptic diversity within Planigale at a time when small mammals are continuing to decline throughout Australia.

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.

A review of the late Cenozoic genus Bohra (Diprotodontia: Macropodidae) and the evolution of tree-kangaroos

Tree-kangaroos of the genus Dendrolagus occupy forest habitats of New Guinea and extreme northeastern Australia, but their evolutionary history is poorly known. Descriptions in the 2000s of near-complete Pleistocene skeletons belonging to larger-bodied species in the now-extinct genus Bohra broadened our understanding of morphological variation in the group and have since helped us to identify unassigned fossils in museum collections, as well as to reassign species previously placed in other genera. Here we describe these fossils and analyse tree-kangaroo systematics via comparative osteology. Including B. planei sp. nov., B. bandharr comb. nov. and B. bila comb. nov., we recognise the existence of at least seven late Cenozoic species of Bohra, with a maximum of three in any one assemblage. All tree-kangaroos (Dendrolagina subtribe nov.) exhibit skeletal adaptations reflective of greater joint flexibility and manoeuvrability, particularly in the hindlimb, compared with other macropodids. The Pliocene species of Bohra retained the stepped calcaneocuboid articulation characteristic of ground-dwelling macropodids, but this became smoothed to allow greater hindfoot rotation in the later species of Bohra and in Dendrolagus. Tree-kangaroo diversification may have been tied to the expansion of forest habitats in the early Pliocene. Following the onset of late Pliocene aridity, some tree-kangaroo species took advantage of the consequent spread of more open habitats, becoming among the largest late Cenozoic tree-dwellers on the continent. Arboreal Old World primates and late Quaternary lemurs may be the closest ecological analogues to the species of Bohra.

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.

Evolution of visual guanylyl cyclases and their activating proteins with respect to clade and species-specific visual system adaptation

Membrane guanylyl cyclase receptors are important regulators of local cGMP production, critically influencing cell growth and differentiation as well as ion transport, blood pressure and calcium feedback of vertebrate phototransduction. Currently, seven different subtypes of membrane guanylyl cyclase receptors have been characterized. These receptors have tissue specific expression and are activated either by small extracellular ligands, changing CO2 concentrations or, in the case of visual guanylyl cyclases, intracellularly interacting Ca2+-dependent activating proteins. In this report, we focus on the visual guanylyl cyclase receptors (GCs) GC-E (gucy2d/e) and GC-F (gucy2f) and their activating proteins (GCAP1/2/3; guca1a/b/c). While gucy2d/e has been detected in all analyzed vertebrates, GC-F receptors are missing in several clades (reptiles, birds, and marsupials) and/or individual species. Interestingly, the absence of GC-F in highly visual sauropsida species with up to 4 different cone-opsins is compensated by an increased number of guanylyl cyclase activating proteins, whereas in nocturnal or visually impaired species with reduced spectral sensitivity it is consolidated by the parallel inactivation of these activators. In mammals, the presence of GC-E and GC-F is accompanied by the expression of one to three GCAPs, whereas in lizards and birds, up to five different GCAPs are regulating the activity of the single GC-E visual membrane receptor. In several nearly blind species, a single GC-E enzyme is often accompanied by a single variant of GCAP, suggesting that one cyclase and one activating protein are both sufficient and required for conferring the basic detection of light.

Total evidence phylogeny of platyrrhine primates and a comparison of undated and tip-dating approaches

There have been multiple published phylogenetic analyses of platyrrhine primates (New World monkeys) using both morphological and molecular data, but relatively few that have integrated both types of data into a total evidence approach. Here, we present phylogenetic analyses of recent and fossil platyrrhines, based on a total evidence data set of 418 morphological characters and 10.2 kilobases of DNA sequence data from 17 nuclear genes taken from previous studies, using undated and tip-dating approaches in a Bayesian framework. We compare the results of these analyses with molecular scaffold analyses using maximum parsimony and Bayesian approaches, and we use a formal information theoretic approach to identify unstable taxa. After a posteriori pruning of unstable taxa, the undated and tip-dating topologies appear congruent with recent molecular analyses and support largely similar relationships, with strong support for Stirtonia as a stem alouattine, Neosaimiri as a stem saimirine, Cebupithecia as a stem pitheciine, and Lagonimico as a stem callitrichid. Both analyses find three Greater Antillean subfossil platyrrhines (Xenothrix, Antillothrix, and Paralouatta) to form a clade that is related to Callicebus, congruent with a single dispersal event by the ancestor of this clade to the Greater Antilles. They also suggest that the fossil Proteropithecia may not be closely related to pitheciines, and that all known platyrrhines older than the Middle Miocene are stem taxa. Notably, the undated analysis found the Early Miocene Panamacebus (currently recognized as the oldest known cebid) to be unstable, and the tip-dating analysis placed it outside crown Platyrrhini. Our tip-dating analysis supports a late Oligocene or earliest Miocene (20.8-27.0 Ma) age for crown Platyrrhini, congruent with recent molecular clock analyses.

Communal nesting is the optimal strategy for heat conservation in a social marsupial: lessons from biophysical models

Endothermy, understood as the maintenance of continuous and high body temperatures owing to the combination of metabolic heat production and an insulative cover, is severely challenged in small endotherms inhabiting cold environments. As a response, social clustering combined with nest use (=communal nesting) is a common strategy for heat conservation. To quantify the actual amount of energy that is saved by this strategy, we studied the social marsupial Dromiciops gliroides (monito del monte), an endemic species of the cold forests of southern South America. It is hypothesized that sociability in this marsupial was driven by cold conditions, but evidence supporting this hypothesis is unclear. Here, we used taxidermic models ('mannequins') to experimentally test the energetic benefits of clustering combined with nest use. To do this, we fitted and compared cooling curves of solitary and grouped mannequins, within and outside of a nest, at the typical winter ambient temperatures of their habitat (5°C). We found that the strategy that minimized euthermic cost of maintenance was the combination of nest use and clustering, thus supporting communal nesting as a social adaptation to cope with the cold. Considering the basal metabolic rate of monitos, our estimates suggest that the savings represents almost half of energy consumption per day (in resting conditions). This study shows how simple biophysical models could help to evaluate bioenergetic hypotheses for social behavior in cold-adapted endotherms.

Reunion of Australasian Possums by Shared SINE Insertions

Although first posited to be of a single origin, the two superfamilies of phalangeriform marsupial possums (Phalangeroidea: brushtail possums and cuscuses and Petauroidea: possums and gliders) have long been considered, based on multiple sequencing studies, to have evolved from two separate origins. However, previous data from these sequence analyses suggested a variety of conflicting trees. Therefore, we reinvestigated these relationships by screening $\sim$200,000 orthologous short interspersed element (SINE) loci across the newly available whole-genome sequences of phalangeriform species and their relatives. Compared to sequence data, SINE presence/absence patterns are evolutionarily almost neutral molecular markers of the phylogenetic history of species. Their random and highly complex genomic insertion ensures their virtually homoplasy-free nature and enables one to compare hundreds of shared unique orthologous events to determine the true species tree. Here, we identify 106 highly reliable phylogenetic SINE markers whose presence/absence patterns within multiple Australasian possum genomes unexpectedly provide the first significant evidence for the reunification of Australasian possums into one monophyletic group. Together, our findings indicate that nucleotide homoplasy and ancestral incomplete lineage sorting have most likely driven the conflicting signal distributions seen in previous sequence-based studies. [Ancestral incomplete lineage sorting; possum genomes; possum monophyly; retrophylogenomics; SINE presence/absence.].

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.

Incomplete lineage sorting and phenotypic evolution in marsupials

Incomplete lineage sorting (ILS) makes ancestral genetic polymorphisms persist during rapid speciation events, inducing incongruences between gene trees and species trees. ILS has complicated phylogenetic inference in many lineages, including hominids. However, we lack empirical evidence that ILS leads to incongruent phenotypic variation. Here, we performed phylogenomic analyses to show that the South American monito del monte is the sister lineage of all Australian marsupials, although over 31% of its genome is closer to the Diprotodontia than to other Australian groups due to ILS during ancient radiation. Pervasive conflicting phylogenetic signals across the whole genome are consistent with some of the morphological variation among extant marsupials. We detected hundreds of genes that experienced stochastic fixation during ILS, encoding the same amino acids in non-sister species. Using functional experiments, we confirm how ILS may have directly contributed to hemiplasy in morphological traits that were established during rapid marsupial speciation ca. 60 mya.

A Chromosome-Scale Hybrid Genome Assembly of the Extinct Tasmanian Tiger (Thylacinus cynocephalus)

The extinct Tasmanian tiger or thylacine (Thylacinus cynocephalus) was a large marsupial carnivore native to Australia. Once ranging across parts of the mainland, the species remained only on the island of Tasmania by the time of European colonization. It was driven to extinction in the early 20th century and is an emblem of native species loss in Australia. The thylacine was a striking example of convergent evolution with placental canids, with which it shared a similar skull morphology. Consequently, it has been the subject of extensive study. While the original thylacine assemblies published in 2018 enabled the first exploration of the species' genome biology, further progress is hindered by the lack of high-quality genomic resources. Here, we present a new chromosome-scale hybrid genome assembly for the thylacine, which compares favorably with many recent de novo marsupial genomes. In addition, we provide homology-based gene annotations, characterize the repeat content of the thylacine genome, and show that consistent with demographic decline, the species possessed a low rate of heterozygosity even compared to extant, threatened marsupials.

Mitogenome of the extinct Desert 'rat-kangaroo' times the adaptation to aridity in macropodoids

The evolution of Australia's distinctive marsupial fauna has long been linked to the onset of continent-wide aridity. However, how this profound climate change event affected the diversification of extant lineages is still hotly debated. Here, we assemble a DNA sequence dataset of Macropodoidea-the clade comprising kangaroos and their relatives-that incorporates a complete mitogenome for the Desert 'rat-kangaroo', Caloprymnus campestris. This enigmatic species went extinct nearly 90 years ago and is known from a handful of museum specimens. Caloprymnus is significant because it was the only macropodoid restricted to extreme desert environments, and therefore calibrates the group's specialisation for increasingly arid conditions. Our robustly supported phylogenies nest Caloprymnus amongst the bettongs Aepyprymnus and Bettongia. Dated ancestral range estimations further reveal that the Caloprymnus-Bettongia lineage originated in nascent xeric settings during the middle to late Miocene, ~ 12 million years ago (Ma), but subsequently radiated into fragmenting mesic habitats after the Pliocene to mid-Pleistocene. This timeframe parallels the ancestral divergences of kangaroos in woodlands and forests, but predates their adaptive dispersal into proliferating dry shrublands and grasslands from the late Miocene to mid-Pleistocene, after ~ 7 Ma. We thus demonstrate that protracted changes in both climate and vegetation likely staged the emergence of modern arid zone macropodoids.

Microbial biogeography of the wombat gastrointestinal tract

Most herbivorous mammals have symbiotic microbes living in their gastrointestinal tracts that help with harvesting energy from recalcitrant plant fibre. The bulk of research into these microorganisms has focused on samples collected from faeces, representing the distal region of the gastrointestinal (GI) tract. However, the GI tract in herbivorous mammals is typically long and complex, containing different regions with distinct physico-chemical properties that can structure resident microbial communities. Little work has been done to document GI microbial communities of herbivorous animals at these sites. In this study, we use 16S rRNA gene sequencing to characterize the microbial biogeography along the GI tract in two species of wombats. Specifically, we survey the microbes along four major gut regions (stomach, small intestine, proximal colon, distal colon) in a single bare-nosed wombat (Vombatus ursinus) and a single southern hairy-nosed wombat (Lasiorhinus latifrons). Our preliminary results show that GI microbial communities of wombats are structured by GI region. For both wombat individuals, we observed a trend of increasing microbial diversity from stomach to distal colon. The microbial composition in the first proximal colon region was more similar between wombat species than the corresponding distal colon region in the same species. We found several microbial genera that were differentially abundant between the first proximal colon (putative site for primary plant fermentation) and distal colon regions (which resemble faecal samples). Surprisingly, only 10.6% (98) and 18.8% (206) of amplicon sequence variants (ASVs) were shared between the first proximal colon region and the distal colon region for the bare-nosed and southern hairy-nosed wombat, respectively. These results suggest that microbial communities in the first proximal colon region-the putative site of primary plant fermentation in wombats-are distinct from the distal colon, and that faecal samples may have limitations in capturing the diversity of these communities. While faeces are still a valuable and effective means of characterising the distal colon microbiota, future work seeking to better understand how GI microbiota impact the energy economy of wombats (and potentially other hindgut-fermenting mammals) may need to take gut biogeography into account.

A Burst of Numt Insertion in the Dasyuridae Family During Marsupial Evolution

Nuclear pseudogenes of mitochondrial origin (numts) are common in all eukaryotes. Our previous scan of numts in sequenced nuclear genomes suggested that the highest numt content currently known in animals is that in the gray short-tailed opossum. The present work sought to determine numt content in marsupials and to compare it to those in placental and monothematic mammals as well as in non-mammalian vertebrates. To achieve this, 70 vertebrate species with available nuclear and mitochondrial genomes were scanned for numt content. An extreme numt content was found in the Dasyuridae, with 3,450 in Sarcophilus harrisii (1,955 kb) and 2,813 in Antechinus flavipes (847 kb). The evolutionarily closest species analyzed, the extinct Thylacinus cynocephalus belonging to the Thylacindae family, had only 435 numts (238 kb). These two Dasyuridae genomes featured the highest numt content identified in animals to date. A phylogenetic analysis of numts longer than 300 bp, using a Diprotodonita mitochondrial tree, indicated a burst of numt insertion that began before the divergence of the Dasyurini and Phascogalini, reaching a peak in the early evolution of the two tribes. No comparable increase was found in the early divergent species T. cynocephalus. Divergence of the Dasyuridae tribes has been previously dated to shortly after the Miocene climate transition, characterized by a rapid temperature decline. Interestingly, deviation from optimal growth temperature is one of the environmental factors reported to increase numt insertions in a laboratory setting.

A chromosome-level genome of Antechinus flavipes provides a reference for an Australian marsupial genus with male death after mating

The 15 species of small carnivorous marsupials that comprise the genus Antechinus exhibit semelparity, a rare life-history strategy in mammals where synchronized death occurs after one breeding season. Antechinus males, but not females, age rapidly (demonstrate organismal senescence) during the breeding season and show promise as new animal models of ageing. Some antechinus species are also threatened or endangered. Here, we report a chromosome-level genome of a male yellow-footed antechinus Antechinus flavipes. The genome assembly has a total length of 3.2 Gb with a contig N50 of 51.8 Mb and a scaffold N50 of 636.7 Mb. We anchored and oriented 99.7% of the assembly on seven pseudochromosomes and found that repetitive DNA sequences occupy 51.8% of the genome. Draft genome assemblies of three related species in the subfamily Phascogalinae, two additional antechinus species (Antechinus argentus and A. arktos) and the iteroparous sister species Murexia melanurus, were also generated. Preliminary demographic analysis supports the hypothesis that climate change during the Pleistocene isolated species in Phascogalinae and shaped their population size. A transcriptomic profile across the A. flavipes breeding season allowed us to identify genes associated with aspects of the male die-off. The chromosome-level A. flavipes genome provides a steppingstone to understanding an enigmatic life-history strategy and a resource to assist the conservation of antechinuses.

The ecology and evolution of the monito del monte, a relict species from the southern South America temperate forests

The arboreal marsupial monito del monte (genus Dromiciops, with two recognized species) is a paradigmatic mammal. It is the sole living representative of the order Microbiotheria, the ancestor lineage of Australian marsupials. Also, this marsupial is the unique frugivorous mammal in the temperate rainforest, being the main seed disperser of several endemic plants of this ecosystem, thus acting as keystone species. Dromiciops is also one of the few hibernating mammals in South America, spending half of the year in a physiological dormancy where metabolism is reduced to 10% of normal levels. This capacity to reduce energy expenditure in winter contrasts with the enormous energy turnover rate they experience in spring and summer. The unique life history strategies of this living Microbiotheria, characterized by an alternation of life in the slow and fast lanes, putatively represent ancestral traits that permitted these cold-adapted mammals to survive in this environment. Here, we describe the ecological role of this emblematic marsupial, summarizing the ecophysiology of hibernation and sociality, updated phylogeographic relationships, reproductive cycle, trophic relationships, mutualisms, conservation, and threats. This marsupial shows high densities, despite presenting slow reproductive rates, a paradox explained by the unique characteristics of its three-dimensional habitat. We finally suggest immediate actions to protect these species that may be threatened in the near future due to habitat destruction and climate change.

A species-level timeline of mammal evolution integrating phylogenomic data

High-throughput sequencing projects generate genome-scale sequence data for species-level phylogenies1-3. However, state-of-the-art Bayesian methods for inferring timetrees are computationally limited to small datasets and cannot exploit the growing number of available genomes4. In the case of mammals, molecular-clock analyses of limited datasets have produced conflicting estimates of clade ages with large uncertainties5,6, and thus the timescale of placental mammal evolution remains contentious7-10. Here we develop a Bayesian molecular-clock dating approach to estimate a timetree of 4,705 mammal species integrating information from 72 mammal genomes. We show that increasingly larger phylogenomic datasets produce diversification time estimates with progressively smaller uncertainties, facilitating precise tests of macroevolutionary hypotheses. For example, we confidently reject an explosive model of placental mammal origination in the Palaeogene8 and show that crown Placentalia originated in the Late Cretaceous with unambiguous ordinal diversification in the Palaeocene/Eocene. Our Bayesian methodology facilitates analysis of complete genomes and thousands of species within an integrated framework, making it possible to address hitherto intractable research questions on species diversifications. This approach can be used to address other contentious cases of animal and plant diversifications that require analysis of species-level phylogenomic datasets.

Divergent locomotor evolution in "giant" kangaroos: Evidence from foot bone bending resistances and microanatomy

The extinct sthenurine (giant, short-faced) kangaroos have been proposed to have a different type of locomotor behavior to extant (macropodine) kangaroos, based both on physical limitations (the size of many exceeds the proposed limit for hopping) and anatomical features (features of the hind limb anatomy suggestive of weight-bearing on one leg at a time). Here, we use micro computerised tomography (micro-CT) scans of the pedal bones of six kangaroos, three sthenurine, and three macropodine, ranging from ~50 to 150 kg, to investigate possible differences in bone resistances to bending and cortical bone distribution that might relate to differences in locomotion. Using second moment of area analysis, we show differences in resistance to bending between the two subfamilies. Distribution of cortical bone shows that sthenurines had less resistant calcaneal tubers, implying a different foot posture during locomotion, and the long foot bones were more resistant to the medial bending stresses. These differences were the most pronounced between Pleistocene monodactyl sthenurines (Sthenurus stirlingi and Procoptodon browneorum) and the two species of Macropus (the extant M. giganteus and the extinct M. cf. M. titan) and support the hypothesis that these derived sthenurines employed bipedal striding. The Miocene sthenurine Hadronomas retains some more macropodine-like features of bone resistance to bending, perhaps reflecting its retention of the fifth pedal digit. The Pleistocene macropodine Protemnodon has a number of unique features, possibly indicative of a type of locomotion unlike the other kangaroos.

Molecular relationships of the red-bellied dasyure (Phascolosorex doriae) – a rare marsupial from western New Guinea

The mitochondrial genome of the rare endemic New Guinean dasyurid Phascolosorex doriae (Thomas 1886) has been used to clarify relationships within ‘phascolosoricinae’. The mitochondrial genome has the typical gene arrangement seen in other marsupials. Molecular analyses using complete mitogenomes of other dasyurids resolve the red-bellied dasyure as sister to the narrow-striped dasyure Phascolosorex dorsalis and show that these two species diverged in the early Pliocene. The invasion of emergent New Guinean rainforest habitats (in the late Miocene) by the common ancestor of Ph. doriae, Ph. dorsalis and Neophascogale lorentzii represents one of three separate such invasions by dasyurid lineages.

The value of updating GenBank accessions for supermatrix phylogeny: The case of the New Guinean marsupial carnivore genus Myoictis

Erroneous taxonomic attributions in GenBank accessions can mislead phylogenetic inference and appear to be widespread within genera. We investigate the influence of taxonomic misattributions for reconstructing the phylogeny of three-striped dasyures, which include four recognized Myoictis species (Marsupialia: Dasyuridae) that are distributed across New Guinea and nearby islands. Molecular phylogenetic studies that have focused on dasyurids consistently resolve the interrelationships of these small carnivores, grouping M. leucura with M. wavicus, and placing M. wallacei and M. melas as successively deeper divergences from these. Two recent marsupial and mammalian supermatrix phylogenies instead favour an alternative Myoictis topology that is discordant with each of these relationships. We add new nuclear and mitochondrial sequences and employ randomized accession resampling that shows the supermatrix topologies are an artefact of several outdated taxonomic attributions in GenBank. Updating these accessions brings agreement across Myoictis phylogenies with randomly resampled accessions. We encourage authors to update GenBank taxonomic attributions and we argue that an option is needed for flagging accessions that are not demonstrably incorrect, but that provide anomalous results. This would serve both as a caution for future supermatrix construction and to highlight accessions of potentially significant biological interest for further study.

Repeated dietary shifts in elapid snakes (Squamata: Elapidae) revealed by ancestral state reconstruction

Identifying the traits of ancestral organisms can reveal patterns and drivers of organismal diversification. Unfortunately, reconstructing complex multistate traits (such as diet) remains challenging. Adopting a ‘reconstruct, then aggregate’ approach in a maximum likelihood framework, we reconstructed ancestral diets for 298 species of elapid snakes. We tested whether different prey types were correlated with one another, tested for one-way contingency between prey type pairs, and examined the relationship between snake body size and dietary composition. We demonstrate that the evolution of diet was characterized by niche conservation punctuated by repeated dietary shifts. The ancestor of elapids most likely fed on reptiles and possibly amphibians, with deviations from this ancestral diet occurring repeatedly due to shifts into marine environments and changes in body size. Moreover, we demonstrate important patterns of prey use, including one-way dependency—most obviously the inclusion of eggs being dependent on a diet that already included the producers of those eggs. Despite imperfect dietary data, our approach produced a robust overview of dietary evolution. Given the paucity of natural history information for many organisms, our approach has the potential to increase the number of lineages to which ancestral state reconstructions of multistate traits can be robustly applied.

Heterothermy as the Norm, Homeothermy as the Exception: Variable Torpor Patterns in the South American Marsupial Monito del Monte (Dromiciops gliroides)

Hibernation (i.e., multiday torpor) is considered an adaptive strategy of mammals to face seasonal environmental challenges such as food, cold, and/or water shortage. It has been considered functionally different from daily torpor, a physiological strategy to cope with unpredictable environments. However, recent studies have shown large variability in patterns of hibernation and daily torpor ("heterothermic responses"), especially in species from tropical and subtropical regions. The arboreal marsupial "monito del monte" (Dromiciops gliroides) is the last living representative of the order Microbiotheria and is known to express both short torpor episodes and also multiday torpor depending on environmental conditions. However, only limited laboratory experiments have documented these patterns in D. gliroides. Here, we combined laboratory and field experiments to characterize the heterothermic responses in this marsupial at extreme temperatures. We used intraperitoneal data loggers and simultaneous measurement of ambient and body temperatures (T A and T B, respectively) for analyzing variations in the thermal differential, in active and torpid animals. We also explored how this differential was affected by environmental variables (T A, natural photoperiod changes, food availability, and body mass changes), using mixed-effects generalized linear models. Our results suggest that: (1) individuals express short bouts of torpor, independently of T A and even during the reproductive period; (2) seasonal torpor also occurs in D. gliroides, with a maximum bout duration of 5 days and a mean defended T B of 3.6 ± 0.9°C (one individual controlled T B at 0.09°C, at sub-freezing T A); (3) the best model explaining torpor occurrence (Akaike information criteria weight = 0.59) discarded all predictor variables except for photoperiod and a photoperiod by food interaction. Altogether, these results confirm that this marsupial expresses a dynamic form of torpor that progresses from short torpor to hibernation as daylength shortens. These data add to a growing body of evidence characterizing tropical and sub-tropical heterothermy as a form of opportunistic torpor, expressed as daily or seasonal torpor depending on environmental conditions.

Amazonia as the Origin and Diversification Area of Didelphidae (Mammalia: Metatheria), and a Review of the Fossil Record of the Clade

Didelphidae is the largest New World radiation of marsupials, and is mostly represented by arboreal, small- to medium-sized taxa that inhabit tropical and/or subtropical forests. The group originated and remained isolated in South America for millions of years, until the formation of the Isthmus of Panama. In this study, we present the first reconstruction of the biogeographic history of Didelphidae including all major clades, based on parametric models and stratified analyses over time. We also compiled all the pre-Quaternary fossil records of the group, and contrasted these data to our biogeographic inferences, as well as to major environmental events that occurred in the South American Cenozoic. Our results indicate the relevance of Amazonia in the early diversification of Didelphidae, including the divergence of the major clades traditionally ranked as subfamilies and tribes. Cladogeneses in other areas started in the late Miocene, an interval of intense shifts, especially in the northern portion of Andes and Amazon Basin. Occupation of other areas continued through the Pliocene, but few were only colonized in Quaternary times. The comparison between the biogeographic inference and the fossil records highlights some further steps towards better understanding the spatiotemporal evolution of the clade. Finally, our results stress that the early history of didelphids is obscured by the lack of Paleogene fossils, which are still to be unearthed from low-latitude deposits of South America.

The biogeography of Dromiciops in southern South America: Middle Miocene transgressions, speciation and associations with Nothofagus

The current distribution of the flora and fauna of southern South America is the result of drastic geological events that occurred during the last 20 million years, including marine transgressions, glaciations and active vulcanism. All these have been associated with fragmentation, isolation and subsequent expansion of the biota, south of 35°S, such as the temperate rainforest. This forest is mostly dominated by Nothofagus trees and is the habitat of the relict marsupial monito del monte, genus Dromiciops, sole survivor of the order Microbiotheria. Preliminary analyses using mtDNA proposed the existence of three main Dromiciops lineages, distributed latitudinally, whose divergence was initially attributed to recent Pleistocene glaciations. Using fossil-calibrated dating on nuclear and mitochondrial genes, here we reevaluate this hypothesis and report an older (Miocene) biogeographic history for the genus. We performed phylogenetic reconstructions using sequences from two mitochondrial DNA and four nuclear DNA genes in 159 specimens from 31 sites across Chile and Argentina. Our phylogenetic analysis resolved three main clades with discrete geographic distributions. The oldest and most differentiated clade corresponds to that of the northern distribution (35.2°S to 39.3°S), which should be considered a distinct species (D. bozinovici, sensu D'Elía et al. 2016). According to our estimations, this species shared a common ancestor with D. gliroides (southern clades) about ~13 million years ago. Divergence time estimates for the southern clades (39.6°S to 42.0°S) ranged from 9.57 to 6.5 Mya. A strong genetic structure was also detected within and between clades. Demographic analyses suggest population size stability for the northern clade (D. bozinovici), and recent demographic expansions for the central and southern clades. All together, our results suggest that the diversification of Dromiciops were initiated by the Middle Miocene transgression (MMT), the massive marine flooding that covered several lowlands of the western face of Los Andes between 37 and 48°S. The MMT resulted from an increase in global sea levels at the Miocene climatic optimum, which shaped the biogeographic origin of several species, including Nothofagus forests, the habitat of Dromiciops.

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

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

Telomeres, species differences, and unusual telomeres in vertebrates: presenting challenges and opportunities to understanding telomere dynamics

There has been increasing interest in the use of telomeres as biomarkers of stress, cellular ageing and life-histories. However, the telomere landscape is a diverse feature, with noticeable differences between species, a fact which is highlighted by the unusual telomeres of various vertebrate organisms. We broadly review differences in telomere dynamics among vertebrates, and emphasize the need to understand more about telomere processes and trends across species. As part of these species differences, we review unusual telomeres in vertebrates. This includes mega-telomeres, which are present across a diverse set of organisms, but also focusing on the unusual telomeres traits of marsupials and monotremes, which have seen little to no prior discussion, yet uniquely stand out from other unusual telomere features discovered thus far. Due to the presence of at least two unique telomere features in the marsupial family Dasyuridae, as well as to the presence of physiological strategies semelparity and torpor, which have implications for telomere life-histories in these species, we suggest that this family has a very large potential to uncover novel information on telomere evolution and dynamics.

Testing hypotheses of marsupial brain size variation using phylogenetic multiple imputations and a Bayesian comparative framework

Considerable controversy exists about which hypotheses and variables best explain mammalian brain size variation. We use a new, high-coverage dataset of marsupial brain and body sizes, and the first phylogenetically imputed full datasets of 16 predictor variables, to model the prevalent hypotheses explaining brain size evolution using phylogenetically corrected Bayesian generalized linear mixed-effects modelling. Despite this comprehensive analysis, litter size emerges as the only significant predictor. Marsupials differ from the more frequently studied placentals in displaying a much lower diversity of reproductive traits, which are known to interact extensively with many behavioural and ecological predictors of brain size. Our results therefore suggest that studies of relative brain size evolution in placental mammals may require targeted co-analysis or adjustment of reproductive parameters like litter size, weaning age or gestation length. This supports suggestions that significant associations between behavioural or ecological variables with relative brain size may be due to a confounding influence of the extensive reproductive diversity of placental mammals.

Dispersal syndromes drive the formation of biogeographical regions, illustrated by the case of Wallace's Line

AIM

Biogeographical regions (realms) reflect patterns of co-distributed species (biotas) across space. Their boundaries are set by dispersal barriers and difficulties of establishment in new locations. We extend new methods to assess these two contributions by quantifying the degree to which realms intergrade across geographical space and the contributions of individual species to the delineation of those realms. As our example, we focus on Wallace's Line, the most enigmatic partitioning of the world's faunas, where climate is thought to have little effect and the majority of dispersal barriers are short water gaps.

LOCATION

Indo-Pacific.

TIME PERIOD

Present day.

MAJOR TAXA STUDIED

Birds and mammals.

METHODS

Terrestrial bird and mammal assemblages were established in 1-degree map cells using range maps. Assemblage structure was modelled using latent Dirichlet allocation, a continuous clustering method that simultaneously establishes the likely partitioning of species into biotas and the contribution of biotas to each map cell. Phylogenetic trees were used to assess the contribution of deep historical processes. Spatial segregation between biotas was evaluated across time and space in comparison with numerous hard realm boundaries drawn by various workers.

RESULTS

We demonstrate that the strong turnover between biotas coincides with the north-western extent of the region not connected to the mainland during the Pleistocene, although the Philippines contains mixed contributions. At deeper taxonomic levels, Sulawesi and the Philippines shift to primarily Asian affinities, resulting from transgressions of a few Asian-derived lineages across the line. The partitioning of biotas sometimes produces fragmented regions that reflect habitat. Differences in partitions between birds and mammals reflect differences in dispersal ability.

MAIN CONCLUSIONS

Permanent water barriers have selected for a dispersive archipelago fauna, excluded by an incumbent continental fauna on the Sunda shelf. Deep history, such as plate movements, is relatively unimportant in setting boundaries. The analysis implies a temporally dynamic interaction between a species' intrinsic dispersal ability, physiographic barriers, and recent climate change in the genesis of Earth's biotas.

Host reproductive cycle influences the pouch microbiota of wild southern hairy-nosed wombats (Lasiorhinus latifrons)

Background

Marsupials are born much earlier than placental mammals, with most crawling from the birth canal to the protective marsupium (pouch) to further their development. However, little is known about the microbiology of the pouch and how it changes throughout a marsupial’s reproductive cycle. Here, using stringent controls, we characterized the microbial composition of multiple body sites from 26 wild Southern Hairy-nosed Wombats (SHNWs), including pouch samples from animals at different reproductive stages.

Results

Using qPCR of the 16S rRNA gene we detected a microbial community in the SHNW pouch. We observed significant differences in microbial composition and diversity between the body sites tested, as well as between pouch samples from different reproductive stages. The pouches of reproductively active females had drastically lower microbial diversity (mean ASV richness 19 ± 8) compared to reproductively inactive females (mean ASV richness 941 ± 393) and were dominated by gram positive bacteria from the Actinobacteriota phylum (81.7–90.6%), with the dominant families classified as Brevibacteriaceae, Corynebacteriaceae, Microbacteriaceae, and Dietziaceae. Three of the five most abundant sequences identified in reproductively active pouches had closest matches to microbes previously isolated from tammar wallaby pouches.

Conclusions

This study represents the first contamination-controlled investigation into the marsupial pouch microbiota, and sets a rigorous framework for future pouch microbiota studies. Our results indicate that SHNW pouches contain communities of microorganisms that are substantially altered by the host reproductive cycle. We recommend further investigation into the roles that pouch microorganisms may play in marsupial reproductive health and joey survival.

Supplementary Information

The online version contains supplementary material available at 10.1186/s42523-021-00074-8.