Disentangling the relationship of the Australian marsupial orders using retrotransposon and evolutionary network analyses

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.

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.

The limits of convergence: the roles of phylogeny and dietary ecology in shaping non-avian amniote crania

Cranial morphology is remarkably varied in living amniotes and the diversity of shapes is thought to correspond with feeding ecology, a relationship repeatedly demonstrated at smaller phylogenetic scales, but one that remains untested across amniote phylogeny. Using a combination of morphometric methods, we investigate the links between phylogenetic relationships, diet and skull shape in an expansive dataset of extant toothed amniotes: mammals, lepidosaurs and crocodylians. We find that both phylogeny and dietary ecology have statistically significant effects on cranial shape. The three major clades largely partition morphospace with limited overlap. Dietary generalists often occupy clade-specific central regions of morphospace. Some parallel changes in cranial shape occur in clades with distinct evolutionary histories but similar diets. However, members of a given clade often present distinct cranial shape solutions for a given diet, and the vast majority of species retain the unique aspects of their ancestral skull plan, underscoring the limits of morphological convergence due to ecology in amniotes. These data demonstrate that certain cranial shapes may provide functional advantages suited to particular dietary ecologies, but accounting for both phylogenetic history and ecology can provide a more nuanced approach to inferring the ecology and functional morphology of cryptic or extinct amniotes.

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.

The genomic view of diversification

The process of species diversification is traditionally summarized by a single tree, the species tree, whose reconstruction from molecular data is hindered by frequent conflicts between gene genealogies. Here, we argue that instead of seeing these conflicts as nuisances, we can exploit them to inform the diversification process itself. We adopt a gene-based view of diversification to model the ubiquitous presence of gene flow between diverging lineages, one of the most important processes explaining disagreements among gene trees. We propose a new framework for modelling the joint evolution of gene and species lineages relaxing the hierarchy between the species tree and gene trees inherent to the standard view, as embodied in a popular model known as the multispecies coalescent (MSC). We implement this framework in two alternative models called the gene-based diversification models (GBD): (a) GBD-forward following all evolving genomes through time and (b) GBD-backward based on coalescent theory. They feature four parameters tuning colonization, gene flow, genetic drift and genetic differentiation. We propose an inference method based on differences between gene trees. Applied to two empirical data sets prone to gene flow, we find better support for the GBD-backward model than for the MSC model. Along with the increasing awareness of the extent of gene flow, this work shows the importance of considering the richer signal contained in genomic histories, rather than in the mere species tree, to better apprehend the complex evolutionary history of species.

Genome of the Tasmanian tiger provides insights into the evolution and demography of an extinct marsupial carnivore

The Tasmanian tiger or thylacine (Thylacinus cynocephalus) was the largest carnivorous Australian marsupial to survive into the modern era. Despite last sharing a common ancestor with the eutherian canids ~160 million years ago, their phenotypic resemblance is considered the most striking example of convergent evolution in mammals. The last known thylacine died in captivity in 1936 and many aspects of the evolutionary history of this unique marsupial apex predator remain unknown. Here we have sequenced the genome of a preserved thylacine pouch young specimen to clarify the phylogenetic position of the thylacine within the carnivorous marsupials, reconstruct its historical demography and examine the genetic basis of its convergence with canids. Retroposon insertion patterns placed the thylacine as the basal lineage in Dasyuromorphia and suggest incomplete lineage sorting in early dasyuromorphs. Demographic analysis indicated a long-term decline in genetic diversity starting well before the arrival of humans in Australia. In spite of their extraordinary phenotypic convergence, comparative genomic analyses demonstrated that amino acid homoplasies between the thylacine and canids are largely consistent with neutral evolution. Furthermore, the genes and pathways targeted by positive selection differ markedly between these species. Together, these findings support models of adaptive convergence driven primarily by cis-regulatory evolution.

The evolution of the macrophage-specific enhancer (Fms intronic regulatory element) within the CSF1R locus of vertebrates

The Csf1r locus encodes the receptor for macrophage colony-stimulating factor, which controls the proliferation, differentiation and survival of macrophages. The 300 bp Fms intronic regulatory element (FIRE), within the second intron of Csf1r, is necessary and sufficient to direct macrophage-specific transcription. We have analysed the conservation and divergence of the FIRE DNA sequence in vertebrates. FIRE is present in the same location in the Csf1r locus in reptile, avian and mammalian genomes. Nearest neighbor analysis based upon this element alone largely recapitulates phylogenies inferred from much larger genomic sequence datasets. One core element, containing binding sites for AP1 family and the macrophage-specific transcription factor, PU.1, is conserved from lizards to humans. Around this element, the FIRE sequence is conserved within clades with the most conserved elements containing motifs for known myeloid-expressed transcription factors. Conversely, there is little alignment between clades outside the AP1/PU.1 element. The analysis favours a hybrid between “enhanceosome” and “smorgasbord” models of enhancer function, in which elements cooperate to bind components of the available transcription factor milieu.

Total evidence phylogeny and evolutionary timescale for Australian faunivorous marsupials (Dasyuromorphia)

BACKGROUND

The order Dasyuromorphia is a diverse radiation of faunivorous marsupials, comprising >80 modern species in Australia and New Guinea. It includes dasyurids, the numbat (the myrmecobiid Myrmecobius fasciatus) and the recently extinct thylacine (the thylacinid Thylacinus cyncocephalus). There is also a diverse fossil record of dasyuromorphians and "dasyuromorphian-like" taxa known from Australia. We present the first total evidence phylogenetic analyses of the order, based on combined morphological and molecular data (including a novel set of 115 postcranial characters), to resolve relationships and calculate divergence dates. We use this information to analyse the diversification dynamics of modern dasyuromorphians.

RESULTS

Our morphology-only analyses are poorly resolved, but our molecular and total evidence analyses confidently resolve most relationships within the order, and are strongly congruent with recent molecular studies. Thylacinidae is the first family to diverge within the order, and there is strong support for four tribes within Dasyuridae (Dasyurini, Phascogalini, Planigalini and Sminthopsini). Among fossil taxa, Ankotarinja and Keeuna do not appear to be members of Dasyuromorphia, whilst Barinya and Mutpuracinus are of uncertain relationships within the order. Divergence dates calculated using total evidence tip-and-node dating are younger than both molecular node-dating and total evidence tip-dating, but appear more congruent with the fossil record and are relatively insensitive to calibration strategy. The tip-and-node divergence dates indicate that Dasyurini, Phascogalini and Sminthopsini began to radiate almost simultaneously during the middle-to-late Miocene (11.5-13.1 MYA; composite 95% HPD: 9.5-15.9 MYA); the median estimates for these divergences are shortly after a drop in global temperatures (the middle Miocene Climatic Transition), and coincide with a faunal turnover event in the mammalian fossil record of Australia. Planigalini radiated much later, during the latest Miocene to earliest Pliocene (6.5 MYA; composite 95% HPD: 4.4-8.9 MYA); the median estimates for these divergences coincide with an increase in grass pollen in the Australian palynological record that suggests the development of more open habitats, which are preferred by modern planigale species.

CONCLUSIONS

Our results provide a phylogenetic and temporal framework for interpreting the evolution of modern and fossil dasyuromorphians, but future progress will require a much improved fossil record.

Skeleton of an unusual, cat-sized marsupial relative (Metatheria: Marsupialiformes) from the middle Eocene (Lutetian: 44-43 million years ago) of Turkey

We describe a near-complete, three-dimensionally preserved skeleton of a metatherian (relative of modern marsupials) from the middle Eocene (Lutetian: 44–43 million years ago) Lülük member of the Uzunçarşıdere Formation, central Turkey. With an estimated body mass of 3–4 kg, about the size of a domestic cat (Felis catus) or spotted quoll (Dasyurus maculatus), it is an order of magnitude larger than the largest fossil metatherians previously known from the Cenozoic of the northern hemisphere. This new taxon is characterised by large, broad third premolars that probably represent adaptations for hard object feeding (durophagy), and its craniodental morphology suggests the capacity to generate high bite forces. Qualitative and quantitative functional analyses of its postcranial skeleton indicate that it was probably scansorial and relatively agile, perhaps broadly similar in locomotor mode to the spotted quoll, but with a greater capacity for climbing and grasping. Bayesian phylogenetic analysis of a total evidence dataset comprising 259 morphological characters and 9kb of DNA sequence data from five nuclear protein-coding genes, using both undated and “tip-and-node dating” approaches, place the new taxon outside the marsupial crown-clade, but within the clade Marsupialiformes. It demonstrates that at least one metatherian lineage evolved to occupy the small-medium, meso- or hypo-carnivore niche in the northern hemisphere during the early Cenozoic, at a time when there were numerous eutherians (placentals and their fossil relatives) filling similar niches. However, the known mammal fauna from Uzunçarşıdere Formation appears highly endemic, and geological evidence suggests that this region of Turkey was an island for at least part of the early Cenozoic, and so the new taxon may have evolved in isolation from potential eutherian competitors. Nevertheless, the new taxon reveals previously unsuspected ecomorphological disparity among northern hemisphere metatherians during the first half of the Cenozoic.

RUNX2 repeat variation does not drive craniofacial diversity in marsupials

BACKGROUND

Runt-related transcription factor 2 (RUNX2) is a transcription factor essential for skeletal development. Variation within the RUNX2 polyglutamine / polyalanine (QA) repeat is correlated with facial length within orders of placental mammals and is suggested to be a major driver of craniofacial diversity. However, it is not known if this correlation exists outside of the placental mammals.

RESULTS

Here we examined the correlation between the RUNX2 QA repeat ratio and facial length in the naturally evolving sister group to the placental mammals, the marsupials. Marsupials have a diverse range of facial lengths similar to that seen in placental mammals. Despite their diversity there was almost no variation seen in the RUNX2 QA repeat across individuals spanning the entire marsupial infraclass. The extreme conservation of the marsupial RUNX2 QA repeat indicates it is under strong purifying selection. Despite this, we observed an unexpectedly high level of repeat purity.

CONCLUSIONS

Unlike within orders of placental mammals, RUNX2 repeat variation cannot drive craniofacial diversity in marsupials. We propose conservation of the marsupial RUNX2 QA repeat is driven by the constraint of accelerated ossification of the anterior skeleton to facilitate life in the pouch. Thus, marsupials must utilize alternate pathways to placental mammals to drive craniofacial evolution.

On the development of the chondrocranium and the histological anatomy of the head in perinatal stages of marsupial mammals

An overview of the literature on the chondrocranium of marsupial mammals reveals a relative conservatism in shape and structures. We document the histological cranial anatomy of individuals representing Monodelphis domestica, Dromiciops gliroides, Perameles sp. and Macropus eugenii. The marsupial chondrocranium is generally characterized by the great breadth of the lamina basalis, absence of pila metoptica and large otic capsules. Its most anterior portion (cupula nasi anterior) is robust, and anterior to it there are well-developed tactile sensory structures, functionally important in the neonate. Investigations of ossification centers at and around the nasal septum are needed to trace the presence of certain bones (e.g., mesethmoid, parasphenoid) across marsupial taxa. In many adult marsupials, the tympanic floor is formed by at least three bones: alisphenoid (alisphenoid tympanic process), ectotympanic and petrosal (rostral and caudal tympanic processes); the squamosal also contributes in some diprotodontians. The presence of an entotympanic in marsupials has not been convincingly demonstrated. The tubal element surrounding the auditory tube in most marsupials is fibrous connective tissue rather than cartilage; the latter is the case in most placentals recorded to date. However, we detected fibrocartilage in a late juvenile of Dromiciops, and a similar tissue has been reported for Tarsipes. Contradictory reports on the presence of the tegmen tympani can be found in the literature. We describe a small tegmen tympani in Macropus. Several heterochronic shifts in the timing of development of the chondocranium and associated structures (e.g., nerves, muscles) and in the ossification sequence have been interpreted as largely being influenced by functional requirements related to the altriciality of the newborn marsupial during early postnatal life. Comparative studies of chondocranial development of mammals can benefit from a solid phylogenetic framework, research on non-classical model organisms, and integration with imaging and sectional data derived from computer-tomography.

Cellular Scaling Rules for the Brains of Marsupials: Not as "Primitive" as Expected

In the effort to understand the evolution of mammalian brains, we have found that common relationships between brain structure mass and numbers of nonneuronal (glial and vascular) cells apply across eutherian mammals, but brain structure mass scales differently with numbers of neurons across structures and across primate and nonprimate clades. This suggests that the ancestral scaling rules for mammalian brains are those shared by extant nonprimate eutherians - but do these scaling relationships apply to marsupials, a sister group to eutherians that diverged early in mammalian evolution? Here we examine the cellular composition of the brains of 10 species of marsupials. We show that brain structure mass scales with numbers of nonneuronal cells, and numbers of cerebellar neurons scale with numbers of cerebral cortical neurons, comparable to what we have found in eutherians. These shared scaling relationships are therefore indicative of mechanisms that have been conserved since the first therians. In contrast, while marsupials share with nonprimate eutherians the scaling of cerebral cortex mass with number of neurons, their cerebella have more neurons than nonprimate eutherian cerebella of a similar mass, and their rest of brain has fewer neurons than eutherian structures of a similar mass. Moreover, Australasian marsupials exhibit ratios of neurons in the cerebral cortex and cerebellum over the rest of the brain, comparable to artiodactyls and primates. Our results suggest that Australasian marsupials have diverged from the ancestral Theria neuronal scaling rules, and support the suggestion that the scaling of average neuronal cell size with increasing numbers of neurons varies in evolution independently of the allocation of neurons across structures.

Bandicoot fossils and DNA elucidate lineage antiquity amongst xeric-adapted Australasian marsupials

Bandicoots (Peramelemorphia) are a unique order of Australasian marsupials whose sparse fossil record has been used as prima facie evidence for climate change coincident faunal turnover. In particular, the hypothesized replacement of ancient rainforest-dwelling extinct lineages by antecedents of xeric-tolerant extant taxa during the late Miocene (~10 Ma) has been advocated as a broader pattern evident amongst other marsupial clades. Problematically, however, this is in persistent conflict with DNA phylogenies. We therefore determine the pattern and timing of bandicoot evolution using the first combined morphological + DNA sequence dataset of Peramelemorphia. In addition, we document a remarkably archaic new fossil peramelemorphian taxon that inhabited a latest Quaternary mosaic savannah-riparian forest ecosystem on the Aru Islands of Eastern Indonesia. Our phylogenetic analyses reveal that unsuspected dental homoplasy and the detrimental effects of missing data collectively obscure stem bandicoot relationships. Nevertheless, recalibrated molecular clocks and multiple ancestral area optimizations unanimously infer an early diversification of modern xeric-adapted forms. These probably originated during the late Palaeogene (30-40 Ma) alongside progenitors of other desert marsupials, and thus occupied seasonally dry heterogenous habitats long before the onset of late Neogene aridity.

The Skull of Epidolops ameghinoi from the Early Eocene Itaboraí Fauna, Southeastern Brazil, and the Affinities of the Extinct Marsupialiform Order Polydolopimorphia

The skull of the polydolopimorphian marsupialiform Epidolops ameghinoi is described in detail for the first time, based on a single well-preserved cranium and associated left and right dentaries plus additional craniodental fragments, all from the early Eocene (53–50 million year old) Itaboraí fauna in southeastern Brazil. Notable craniodental features of E. ameghinoi include absence of a masseteric process, very small maxillopalatine fenestrae, a prominent pterygoid fossa enclosed laterally by a prominent ectopterygoid crest, an absent or tiny transverse canal foramen, a simple, planar glenoid fossa, and a postglenoid foramen that is immediately posterior to the postglenoid process. Most strikingly, the floor of the hypotympanic sinus was apparently unossified, a feature found in several stem marsupials but absent in all known crown marsupials. “Type II” marsupialiform petrosals previously described from Itaboraí plausibly belong to E. ameghinoi; in published phylogenetic analyses, these petrosals fell outside (crown-clade) Marsupialia. “IMG VII” tarsals previously referred to E. ameghinoi do not share obvious synapomorphies with any crown marsupial clade, nor do they resemble those of the only other putative polydolopimorphians represented by tarsal remains, namely the argyrolagids. Most studies have placed Polydolopimorphia within Marsupialia, related to either Paucituberculata, or to Microbiotheria and Diprotodontia. However, diprotodonty almost certainly evolved independently in polydolopimorphians, paucituberculatans and diprotodontians, and Epidolops does not share obvious synapomorphies with any marsupial order. Epidolops is dentally specialized, but several morphological features appear to be more plesiomorphic than any crown marsupial. It seems likely Epidolops that falls outside Marsupialia, as do morphologically similar forms such as Bonapartherium and polydolopids. Argyrolagids differ markedly in their known morphology from Epidolops but share some potential apomorphies with paucituberculatans. It is proposed that Polydolopimorphia as currently recognised is polyphyletic, and that argyrolagids (and possibly other taxa currently included in Argyrolagoidea, such as groeberiids and patagoniids) are members of Paucituberculata. This hypothesis is supported by Bayesian non-clock phylogenetic analyses of a total evidence matrix comprising DNA sequence data from five nuclear protein-coding genes, indels, retroposon insertions, and morphological characters: Epidolops falls outside Marsupialia, whereas argyrolagids form a clade with the paucituberculatans Caenolestes and Palaeothentes, regardless of whether the Type II petrosals and IMG VII tarsals are used to score characters for Epidolops or not. There is no clear evidence for the presence of crown marsupials at Itaboraí, and it is possible that the origin and early evolution of Marsupialia was restricted to the “Austral Kingdom” (southern South America, Antarctica, and Australia).

Multiple Evolutionary Origins of Ubiquitous Cu2+ and Zn2+ Binding in the S100 Protein Family

The S100 proteins are a large family of signaling proteins that play critical roles in biology and disease. Many S100 proteins bind Zn2+, Cu2+, and/or Mn2+ as part of their biological functions; however, the evolutionary origins of binding remain obscure. One key question is whether divalent transition metal binding is ancestral, or instead arose independently on multiple lineages. To tackle this question, we combined phylogenetics with biophysical characterization of modern S100 proteins. We demonstrate an earlier origin for established S100 subfamilies than previously believed, and reveal that transition metal binding is widely distributed across the tree. Using isothermal titration calorimetry, we found that Cu2+ and Zn2+ binding are common features of the family: the full breadth of human S100 paralogs-as well as two early-branching S100 proteins found in the tunicate Oikopleura dioica-bind these metals with μM affinity and stoichiometries ranging from 1:1 to 3:1 (metal:protein). While binding is consistent across the tree, structural responses to binding are quite variable. Further, mutational analysis and structural modeling revealed that transition metal binding occurs at different sites in different S100 proteins. This is consistent with multiple origins of transition metal binding over the evolution of this protein family. Our work reveals an evolutionary pattern in which the overall phenotype of binding is a constant feature of S100 proteins, even while the site and mechanism of binding is evolutionarily labile.

When Genomics Is Not Enough: Experimental Evidence for a Decrease in LINE-1 Activity During the Evolution of Australian Marsupials

The autonomous transposable element LINE-1 is a highly abundant element that makes up between 15% and 20% of therian mammal genomes. Since their origin before the divergence of marsupials and placental mammals, LINE-1 elements have contributed actively to the genome landscape. A previous in silico screen of the Tasmanian devil genome revealed a lack of functional coding LINE-1 sequences. In this study we present the results of an in vitro analysis from a partial LINE-1 reverse transcriptase coding sequence in five marsupial species. Our experimental screen supports the in silico findings of the genome-wide degradation of LINE-1 sequences in the Tasmanian devil, and identifies a high frequency of degraded LINE-1 sequences in other Australian marsupials. The comparison between the experimentally obtained LINE-1 sequences and reference genome assemblies suggests that conclusions from in silico analyses of retrotransposition activity can be influenced by incomplete genome assemblies from short reads.