A phylogeny of Diprotodontia (Marsupialia) based on sequences for five nuclear genes

Emx2 underlies the development and evolution of marsupial gliding membranes

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

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.

A macroevolutionary common-garden experiment reveals differentially evolvable bone organization levels in slow arboreal mammals

Eco-morphological convergence, i.e., similar phenotypes evolved in ecologically convergent taxa, naturally reproduces a common-garden experiment since it allows researchers to keep ecological factors constant, studying intrinsic evolutionary drivers. The latter may result in differential evolvability that, among individual anatomical parts, causes mosaic evolution. Reconstructing the evolutionary morphology of the humerus and femur of slow arboreal mammals, we addressed mosaicism at different bone anatomical spatial scales. We compared convergence strength, using it as indicator of evolvability, between bone external shape and inner structure, with the former expected to be less evolvable and less involved in convergent evolution, due to anatomical constraints. We identify several convergent inner structural traits, while external shape only loosely follows this trend, and we find confirmation for our assumption in measures of convergence magnitude. We suggest that future macroevolutionary reconstructions based on bone morphology should include structural traits to better detect ecological effects on vertebrate diversification.

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.

Convergent deployment of ancestral functions during the evolution of mammalian flight membranes

Lateral flight membranes, or patagia, have evolved repeatedly in diverse mammalian lineages. While little is known about patagium development, its recurrent evolution may suggest a shared molecular basis. By combining transcriptomics, developmental experiments, and mouse transgenics, we demonstrate that lateral Wnt5a expression in the marsupial sugar glider (Petaurus breviceps) promotes the differentiation of its patagium primordium. We further show that this function of Wnt5a reprises ancestral roles in skin morphogenesis predating mammalian flight and has been convergently used during patagium evolution in eutherian bats. Moreover, we find that many genes involved in limb development have been redeployed during patagium outgrowth in both the sugar glider and bat. Together, our findings reveal that deeply conserved genetic toolkits contribute to the evolutionary transition to flight in mammals.

Time-calibrated phylogeny and ecological niche models indicate Pliocene aridification drove intraspecific diversification of brushtail possums in Australia

Major aridification events in Australia during the Pliocene may have had significant impact on the distribution and structure of widespread species. To explore the potential impact of Pliocene and Pleistocene climate oscillations, we estimated the timing of population fragmentation and past connectivity of the currently isolated but morphologically similar subspecies of the widespread brushtail possum (Trichosurus vulpecula). We use ecological niche modeling (ENM) with the current fragmented distribution of brushtail possums to estimate the environmental envelope of this marsupial. We projected the ENM on models of past climatic conditions in Australia to infer the potential distribution of brushtail possums over 6 million years. D-loop haplotypes were used to describe population structure. From shotgun sequencing, we assembled whole mitochondrial DNA genomes and estimated the timing of intraspecific divergence. Our projections of ENMs suggest current possum populations were unlikely to have been in contact during the Pleistocene. Although lowered sea level during glacial periods enabled connection with habitat in Tasmania, climate fluctuation during this time would not have facilitated gene flow over much of Australia. The most recent common ancestor of sampled intraspecific diversity dates to the early Pliocene when continental aridification caused significant changes to Australian ecology and Trichosurus vulpecula distribution was likely fragmented. Phylogenetic analysis revealed that the subspecies T. v. hypoleucus (koomal; southwest), T. v. arnhemensis (langkurr; north), and T. v. vulpecula (bilda; southeast) correspond to distinct mitochondrial lineages. Despite little phenotypic differentiation, Trichosurus vulpecula populations probably experienced little gene flow with one another since the Pliocene, supporting the recognition of several subspecies and explaining their adaptations to the regional plant assemblages on which they feed.

Metatarsal fusion resisted bending as jerboas (Dipodidae) transitioned from quadrupedal to bipedal

Hind limbs undergo dramatic changes in loading conditions during the transition from quadrupedal to bipedal locomotion. For example, the most early diverging bipedal jerboas (Rodentia: Dipodidae) are some of the smallest mammals in the world, with body masses that range between 2–4 g. The larger jerboa species exhibit developmental and evolutionary fusion of the central three metatarsals into a single cannon bone. We hypothesize that small body size and metatarsal fusion are mechanisms to maintain the safety factor of the hind limb bones despite the higher ground reaction forces associated with bipedal locomotion. Using finite-element analysis to model collisions between the substrate and the metatarsals, we found that body size reduction was insufficient to reduce bone stress on unfused metatarsals, based on the scaled dynamics of larger jerboas, and that fused bones developed lower stresses than unfused bones when all metatarsals are scaled to the same size and loading conditions. Based on these results, we conclude that fusion reinforces larger jerboa metatarsals against high ground reaction forces. Because smaller jerboas with unfused metatarsals develop higher peak stresses in response to loading conditions scaled from larger jerboas, we hypothesize that smaller jerboas use alternative dynamics of bipedal locomotion to reduce the impact of collisions between the foot and substrate.

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

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.

A new family of diprotodontian marsupials from the latest Oligocene of Australia and the evolution of wombats, koalas, and their relatives (Vombatiformes)

We describe the partial cranium and skeleton of a new diprotodontian marsupial from the late Oligocene (~26–25 Ma) Namba Formation of South Australia. This is one of the oldest Australian marsupial fossils known from an associated skeleton and it reveals previously unsuspected morphological diversity within Vombatiformes, the clade that includes wombats (Vombatidae), koalas (Phascolarctidae) and several extinct families. Several aspects of the skull and teeth of the new taxon, which we refer to a new family, are intermediate between members of the fossil family Wynyardiidae and wombats. Its postcranial skeleton exhibits features associated with scratch-digging, but it is unlikely to have been a true burrower. Body mass estimates based on postcranial dimensions range between 143 and 171 kg, suggesting that it was ~5 times larger than living wombats. Phylogenetic analysis based on 79 craniodental and 20 postcranial characters places the new taxon as sister to vombatids, with which it forms the superfamily Vombatoidea as defined here. It suggests that the highly derived vombatids evolved from wynyardiid-like ancestors, and that scratch-digging adaptations evolved in vombatoids prior to the appearance of the ever-growing (hypselodont) molars that are a characteristic feature of all post-Miocene vombatids. Ancestral state reconstructions on our preferred phylogeny suggest that bunolophodont molars are plesiomorphic for vombatiforms, with full lophodonty (characteristic of diprotodontoids) evolving from a selenodont morphology that was retained by phascolarctids and ilariids, and wynyardiids and vombatoids retaining an intermediate selenolophodont condition. There appear to have been at least six independent acquisitions of very large (>100 kg) body size within Vombatiformes, several having already occurred by the late Oligocene.

Phylogenetic relationships of the cuscuses (Diprotodontia : Phalangeridae) of island Southeast Asia and Melanesia based on the mitochondrial ND2 gene

The species-level systematics of the marsupial family Phalangeridae, particularly Phalanger, are poorly understood, due partly to the family’s wide distribution across Australia, New Guinea, eastern Indonesia, and surrounding islands. In order to refine the species-level systematics of Phalangeridae, and improve our understanding of their evolution, we generated 36 mitochondrial ND2 DNA sequences from multiple species and sample localities. We combined our new data with available sequences and produced the most comprehensive molecular phylogeny for Phalangeridae to date. Our analyses (1) strongly support the monophyly of the three phalangerid subfamilies (Trichosurinae, Ailuropinae, Phalangerinae); (2) reveal the need to re-examine all specimens currently identified as ‘Phalanger orientalis’; and (3) suggest the elevation of the Solomon Island P. orientalis subspecies to species level (P. breviceps Thomas, 1888). In addition, samples of P. orientalis from Timor formed a clade, consistent with an introduction by humans from a single source population. However, further research on east Indonesian P. orientalis populations will be required to test this hypothesis, resolve inconsistencies in divergence time estimates, and locate the source population and taxonomic status of the Timor P. orientalis.

Evolution of Marsupial Genomes

Marsupial genomes, which are packaged into large chromosomes, provide a powerful resource for studying the mechanisms of genome evolution. The extensive and valuable body of work on marsupial cytogenetics, combined more recently with genome sequence data, has enabled prediction of the 2n = 14 karyotype ancestral to all marsupial families. The application of both chromosome biology and genome sequencing, or chromosomics, has been a necessary approach for various aspects of mammalian genome evolution, such as understanding sex chromosome evolution and the origin and evolution of transmissible tumors in Tasmanian devils. The next phase of marsupial genome evolution research will employ chromosomics approaches to begin addressing fundamental questions in marsupial genome evolution and chromosome evolution more generally. The answers to these complex questions will impact our understanding across a broad range of fields, including the genetics of speciation, genome adaptation to environmental stressors, and species management.

Koala and Wombat Gammaherpesviruses Encode the First Known Viral NTPDase Homologs and Are Phylogenetically Divergent from All Known Gammaherpesviruses

There is a large taxonomic gap in our understanding of mammalian herpesvirus genetics and evolution corresponding to those herpesviruses that infect marsupials, which diverged from eutherian mammals approximately 150 million years ago (mya). We compare the genomes of two marsupial gammaherpesviruses, Phascolarctid gammaherpesvirus 1 (PhaHV1) and Vombatid gammaherpesvirus 1 (VoHV1), which infect koalas (Phascolarc tos cinereus) and wombats (Vombatus ursinus), respectively. The core viral genomes were approximately 117 kbp and 110 kbp in length, respectively, sharing 69% pairwise nucleotide sequence identity. Phylogenetic analyses showed that PhaHV1 and VoHV1 formed a separate branch, which may indicate a new gammaherpesvirus genus. The genomes contained 60 predicted open reading frames (ORFs) homologous to those in eutherian herpesviruses and 20 ORFs not yet found in any other herpesvirus. Seven of these ORFs were shared by the two viruses, indicating that they were probably acquired prespeciation, approximately 30 to 40 mya. One of these shared genes encodes a putative nucleoside triphosphate diphosphohydrolase (NTPDase). NTPDases are usually found in mammals and higher-order eukaryotes, with a very small number being found in bacteria. This is the first time that an NTPDase has been identified in any viral genome. Interrogation of public transcriptomic data sets from two koalas identified PhaHV1-specific transcripts in multiple host tissues, including transcripts for the novel NTPDase. PhaHV1 ATPase activity was also demonstrated in vitro, suggesting that the encoded NTPDase is functional during viral infection. In mammals, NTPDases are important in downregulation of the inflammatory and immune responses, but the role of the PhaHV1 NTPDase during viral infection remains to be determined.IMPORTANCE The genome sequences of the koala and wombat gammaherpesviruses show that the viruses form a distinct branch, indicative of a novel genus within the Gammaherpesvirinae Their genomes contain several new ORFs, including ORFs encoding a β-galactoside α-2,6-sialyltransferase that is phylogenetically closest to poxvirus and insect homologs and the first reported viral NTPDase. NTPDases are ubiquitously expressed in mammals and are also present in several parasitic, fungal, and bacterial pathogens. In mammals, these cell surface-localized NTPDases play essential roles in thromboregulation, inflammation, and immune suppression. In this study, we demonstrate that the virus-encoded NTPDase is enzymatically active and is transcribed during natural infection of the host. Understanding how these enzymes benefit viruses can help to inform how they may cause disease or evade host immune defenses.

Evolutionary Associations of Endosymbiotic Ciliates Shed Light on the Timing of the Marsupial-Placental Split

Trichostome ciliates are among the most conspicuous protists in the gastrointestinal tract of a large variety of vertebrates. However, little is still known about phylogeny of the trichostome/vertebrate symbiotic systems, evolutionary correlations between trichostome extrinsic traits, and character-dependent diversification of trichostomes. These issues were investigated here, using the relaxed molecular clock technique along with stochastic mapping of character evolution, and binary-state speciation and extinction models. Clock analyses revealed that trichostomes colonized the vertebrate gastrointestinal tract ∼135 Ma, that is, near the paleontological minimum for the split of therian mammals into marsupials and placentals. According to stochastic mapping, the last common ancestor of trichostomes most likely invaded the hindgut of a mammal. Although multiple shifts to fish/amphibian or avian hosts and to the foregut compartments took place during the trichostome phylogeny, only transition to the foregut was recognized as a key innovation responsible for the explosive radiation of ophryoscolecid trichostomes after the Cretaceous/Tertiary boundary, when ungulates began their diversification. Since crown radiations of main trichostome lineages follow those of their mammalian hosts and are in agreement with their historic dispersal routes, the present time-calibrated phylogeny might help to elucidate controversies in the geological and molecular timing of the split between marsupials and placental mammals.

Phylogenetic analysis of the tree-kangaroos (Dendrolagus) reveals multiple divergent lineages within New Guinea

Amongst the Australasian kangaroos and wallabies (Macropodidae) one anomalous genus, the tree-kangaroos, Dendrolagus, has secondarily returned to arboreality. Modern tree-kangaroos are confined to the wet tropical forests of north Queensland, Australia (2 species) and New Guinea (8 species). Due to their behavior, distribution and habitat most species are poorly known and our understanding of the evolutionary history and systematics of the genus is limited and controversial. We obtained tissue samples from 36 individual Dendrolagus including representatives from 14 of the 17 currently recognised or proposed subspecies and generated DNA sequence data from three mitochondrial (3116 bp) and five nuclear (4097 bp) loci. Phylogenetic analysis of these multi-locus data resolved long-standing questions regarding inter-relationships within Dendrolagus. The presence of a paraphyletic ancestral long-footed and derived monophyletic short-footed group was confirmed. Six major lineages were identified: one in Australia (D. lumholtzi, D. bennettianus) and five in New Guinea (D. inustus, D. ursinus, a Goodfellow's group, D. mbaiso and a Doria's group). Two major episodes of diversification within Dendrolagus were identified: the first during the late Miocene/early Pliocene associated with orogenic processes in New Guinea and the second mostly during the early Pleistocene associated with the intensification of climatic cycling. All sampled subspecies showed high levels of genetic divergence and currently recognized species within both the Doria's and Goodfellow's groups were paraphyletic indicating that adjustments to current taxonomy are warranted.

Molecular evolution of DNMT1 in vertebrates: Duplications in marsupials followed by positive selection

DNA methylation is mediated by a conserved family of DNA methyltransferases (Dnmts). The human genome encodes three active Dnmts (Dnmt1, Dnmt3a and Dnmt3b), the tRNA methyltransferase Dnmt2, and the regulatory protein Dnmt3L. Despite their high degree of conservation among different species, genes encoding Dnmts have been duplicated and/or lost in multiple lineages throughout evolution, indicating that the DNA methylation machinery has some potential to undergo evolutionary change. However, little is known about the extent to which this machinery, or the methylome, varies among vertebrates. Here, we study the molecular evolution of Dnmt1, the enzyme responsible for maintenance of DNA methylation patterns after replication, in 79 vertebrate species. Our analyses show that all studied species exhibit a single copy of the DNMT1 gene, with the exception of tilapia and marsupials (tammar wallaby, koala, Tasmanian devil and opossum), each of which displays two apparently functional DNMT1 copies. Our phylogenetic analyses indicate that DNMT1 duplicated before the radiation of major marsupial groups (i.e., at least ~75 million years ago), thus giving rise to two DNMT1 copies in marsupials (copy 1 and copy 2). In the opossum lineage, copy 2 was lost, and copy 1 recently duplicated again, generating three DNMT1 copies: two putatively functional genes (copy 1a and 1b) and one pseudogene (copy 1ψ). Both marsupial copies (DNMT1 copies 1 and 2) are under purifying selection, and copy 2 exhibits elevated rates of evolution and signatures of positive selection, suggesting a scenario of neofunctionalization. This gene duplication might have resulted in modifications in marsupial methylomes and their dynamics.

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.

Functional morphological adaptations of the bony labyrinth in marsupials (Mammalia, Theria)

Diprotodontia represents the largest and ecologically most distinct order of marsupials occurring in Australasian being highly divers in size, locomotion, habitat preferences, feeding, and activity pattern. The spatial orientation in the habitat and therefore the three-dimensional space is detected by the vestibular system of the inner ear, more precisely by the three semicircular canals. In this study, we investigated the bony labyrinth of diprotodontian and selected non-diprotodontian marsupial mammals of almost all genera with noninvasive micro-CT scanning and 3D-reconstructions. In principal component analyses, the subterranean taxon can be separated from gliding and saltatorial taxa, whereas arboreal species can be separated from saltatorial specimens. The highest PCA loadings of this functional distinction are clearly found in the diameter of the semicircular canals, whereas the overall shape (height, width, length) of the semicircular canals is less important. Additionally, the investigated arboreal and fossorial species of South America are nested in the morphospace of the Australasian taxa. Even if a phylogenetic signal in the anatomy of the bony labyrinth cannot be excluded entirely, the main functional morphological signal of the vestibular system is found in the diameter of the semicircular canals. With the large dataset of extant marsupial mammals analysed here, the locomotion mode of extinct taxa can be inferred in future studies independent of any evidence of postcranial material.

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.

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

Germline viral "fossils" guide in silico reconstruction of a mid-Cenozoic era marsupial adeno-associated virus

Germline endogenous viral elements (EVEs) genetically preserve viral nucleotide sequences useful to the study of viral evolution, gene mutation, and the phylogenetic relationships among host organisms. Here, we describe a lineage-specific, adeno-associated virus (AAV)-derived endogenous viral element (mAAV-EVE1) found within the germline of numerous closely related marsupial species. Molecular screening of a marsupial DNA panel indicated that mAAV-EVE1 occurs specifically within the marsupial suborder Macropodiformes (present-day kangaroos, wallabies, and related macropodoids), to the exclusion of other Diprotodontian lineages. Orthologous mAAV-EVE1 locus sequences from sixteen macropodoid species, representing a speciation history spanning an estimated 30 million years, facilitated compilation of an inferred ancestral sequence that recapitulates the genome of an ancient marsupial AAV that circulated among Australian metatherian fauna sometime during the late Eocene to early Oligocene. In silico gene reconstruction and molecular modelling indicate remarkable conservation of viral structure over a geologic timescale. Characterisation of AAV-EVE loci among disparate species affords insight into AAV evolution and, in the case of macropodoid species, may offer an additional genetic basis for assignment of phylogenetic relationships among the Macropodoidea. From an applied perspective, the identified AAV “fossils” provide novel capsid sequences for use in translational research and clinical applications.

Comprehensive profiling of retroviral integration sites using target enrichment methods from historical koala samples without an assembled reference genome

Background. Retroviral integration into the host germline results in permanent viral colonization of vertebrate genomes. The koala retrovirus (KoRV) is currently invading the germline of the koala (Phascolarctos cinereus) and provides a unique opportunity for studying retroviral endogenization. Previous analysis of KoRV integration patterns in modern koalas demonstrate that they share integration sites primarily if they are related, indicating that the process is currently driven by vertical transmission rather than infection. However, due to methodological challenges, KoRV integrations have not been comprehensively characterized.

Results. To overcome these challenges, we applied and compared three target enrichment techniques coupled with next generation sequencing (NGS) and a newly customized sequence-clustering based computational pipeline to determine the integration sites for 10 museum Queensland and New South Wales (NSW) koala samples collected between the 1870s and late 1980s. A secondary aim of this study sought to identify common integration sites across modern and historical specimens by comparing our dataset to previously published studies. Several million sequences were processed, and the KoRV integration sites in each koala were characterized.

Conclusions. Although the three enrichment methods each exhibited bias in integration site retrieval, a combination of two methods, Primer Extension Capture and hybridization capture is recommended for future studies on historical samples. Moreover, identification of integration sites shows that the proportion of integration sites shared between any two koalas is quite small.

Mammals from 'down under': a multi-gene species-level phylogeny of marsupial mammals (Mammalia, Metatheria)

Marsupials or metatherians are a group of mammals that are distinct in giving birth to young at early stages of development and in having a prolonged investment in lactation. The group consists of nearly 350 extant species, including kangaroos, koala, possums, and their relatives. Marsupials are an old lineage thought to have diverged from early therian mammals some 160 million years ago in the Jurassic, and have a remarkable evolutionary and biogeographical history, with extant species restricted to the Americas, mostly South America, and to Australasia. Although the group has been the subject of decades of phylogenetic research, the marsupial tree of life remains controversial, with most studies focusing on only a fraction of the species diversity within the infraclass. Here we present the first Methaterian species-level phylogeny to include 80% of the extant marsupial species and five nuclear and five mitochondrial markers obtained from Genbank and a recently published retroposon matrix. Our primary goal is to provide a summary phylogeny that will serve as a tool for comparative research. We evaluate the extent to which the phylogeny recovers current phylogenetic knowledge based on the recovery of “benchmark clades” from prior studies—unambiguously supported key clades and undisputed traditional taxonomic groups. The Bayesian phylogenetic analyses recovered nearly all benchmark clades but failed to find support for the suborder Phalagiformes. The most significant difference with previous published topologies is the support for Australidelphia as a group containing Microbiotheriidae, nested within American marsupials. However, a likelihood ratio test shows that alternative topologies with monophyletic Australidelphia and Ameridelphia are not significantly different than the preferred tree. Although further data are needed to solidify understanding of Methateria phylogeny, the new phylogenetic hypothesis provided here offers a well resolved and detailed tool for comparative analyses, covering the majority of the known species richness of the group.

Late pleistocene Australian marsupial DNA clarifies the affinities of extinct megafaunal kangaroos and wallabies

Understanding the evolution of Australia's extinct marsupial megafauna has been hindered by a relatively incomplete fossil record and convergent or highly specialized morphology, which confound phylogenetic analyses. Further, the harsh Australian climate and early date of most megafaunal extinctions (39-52 ka) means that the vast majority of fossil remains are unsuitable for ancient DNA analyses. Here, we apply cross-species DNA capture to fossils from relatively high latitude, high altitude caves in Tasmania. Using low-stringency hybridization and high-throughput sequencing, we were able to retrieve mitochondrial sequences from two extinct megafaunal macropodid species. The two specimens, Simosthenurus occidentalis (giant short-faced kangaroo) and Protemnodon anak (giant wallaby), have been radiocarbon dated to 46-50 and 40-45 ka, respectively. This is significantly older than any Australian fossil that has previously yielded DNA sequence information. Processing the raw sequence data from these samples posed a bioinformatic challenge due to the poor preservation of DNA. We explored several approaches in order to maximize the signal-to-noise ratio in retained sequencing reads. Our findings demonstrate the critical importance of adopting stringent processing criteria when distant outgroups are used as references for mapping highly fragmented DNA. Based on the most stringent nucleotide data sets (879 bp for S. occidentalis and 2,383 bp for P. anak), total-evidence phylogenetic analyses confirm that macropodids consist of three primary lineages: Sthenurines such as Simosthenurus (extinct short-faced kangaroos), the macropodines (all other wallabies and kangaroos), and the enigmatic living banded hare-wallaby Lagostrophus fasciatus (Lagostrophinae). Protemnodon emerges as a close relative of Macropus (large living kangaroos), a position not supported by recent morphological phylogenetic analyses.

A transcriptome resource for the koala (Phascolarctos cinereus): insights into koala retrovirus transcription and sequence diversity

Background

The koala, Phascolarctos cinereus, is a biologically unique and evolutionarily distinct Australian arboreal marsupial. The goal of this study was to sequence the transcriptome from several tissues of two geographically separate koalas, and to create the first comprehensive catalog of annotated transcripts for this species, enabling detailed analysis of the unique attributes of this threatened native marsupial, including infection by the koala retrovirus.

Results

RNA-Seq data was generated from a range of tissues from one male and one female koala and assembled de novo into transcripts using Velvet-Oases. Transcript abundance in each tissue was estimated. Transcripts were searched for likely protein-coding regions and a non-redundant set of 117,563 putative protein sequences was produced. In similarity searches there were 84,907 (72%) sequences that aligned to at least one sequence in the NCBI nr protein database. The best alignments were to sequences from other marsupials. After applying a reciprocal best hit requirement of koala sequences to those from tammar wallaby, Tasmanian devil and the gray short-tailed opossum, we estimate that our transcriptome dataset represents approximately 15,000 koala genes. The marsupial alignment information was used to look for potential gene duplications and we report evidence for copy number expansion of the alpha amylase gene, and of an aldehyde reductase gene.

Koala retrovirus (KoRV) transcripts were detected in the transcriptomes. These were analysed in detail and the structure of the spliced envelope gene transcript was determined. There was appreciable sequence diversity within KoRV, with 233 sites in the KoRV genome showing small insertions/deletions or single nucleotide polymorphisms. Both koalas had sequences from the KoRV-A subtype, but the male koala transcriptome has, in addition, sequences more closely related to the KoRV-B subtype. This is the first report of a KoRV-B-like sequence in a wild population.

Conclusions

This transcriptomic dataset is a useful resource for molecular genetic studies of the koala, for evolutionary genetic studies of marsupials, for validation and annotation of the koala genome sequence, and for investigation of koala retrovirus. Annotated transcripts can be browsed and queried at http://koalagenome.org.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-786) contains supplementary material, which is available to authorized users.

The aerodynamic performance of the feathertail glider, Acrobates pygmaeus (Marsupialia: Acrobatidae)

Photographic and videographic investigations of the aerodynamic performance of four captive adult Acrobates pygmeus are described. During short (1–4 m) glides, steep angles of descent (>45°) and large angles of attack (up to 50°) were used mid-flight, and the possums failed to achieve steady velocities. During longer (>10 m) glides steady velocities of 5.3–7.5 m s–1 were achieved, and glide angles of 21–42° and angles of attack between 36° and 45° were used. The best (lowest) glide angles used at steady velocity were similar to those documented in sugar gliders, but involved lower air speeds. During steady gliding the forelimbs of Acrobates assumed an ‘elbows-out’ disposition similar to that of Petauroides volans, and unlike that of petaurid gliders. A glide polar for Acrobates shows that it uses lower air speeds and higher sinking speeds than other gliding animals for which polars are available (several birds and one bat). Few data on the radius of turn and sinking speed were obtained; nevertheless, in accord with theory, sinking speeds were higher during tighter turns. During flight the feather-like tail was used to generate pitching movements and adjust the angle of attack and hence glide angle. However, the control of turns seems to be achieved primarily through limb adjustments.

Reconstruction of the ancestral marsupial karyotype from comparative gene maps

BACKGROUND

The increasing number of assembled mammalian genomes makes it possible to compare genome organisation across mammalian lineages and reconstruct chromosomes of the ancestral marsupial and therian (marsupial and eutherian) mammals. However, the reconstruction of ancestral genomes requires genome assemblies to be anchored to chromosomes. The recently sequenced tammar wallaby (Macropus eugenii) genome was assembled into over 300,000 contigs. We previously devised an efficient strategy for mapping large evolutionarily conserved blocks in non-model mammals, and applied this to determine the arrangement of conserved blocks on all wallaby chromosomes, thereby permitting comparative maps to be constructed and resolve the long debated issue between a 2n = 14 and 2n = 22 ancestral marsupial karyotype.

RESULTS

We identified large blocks of genes conserved between human and opossum, and mapped genes corresponding to the ends of these blocks by fluorescence in situ hybridization (FISH). A total of 242 genes was assigned to wallaby chromosomes in the present study, bringing the total number of genes mapped to 554 and making it the most densely cytogenetically mapped marsupial genome. We used these gene assignments to construct comparative maps between wallaby and opossum, which uncovered many intrachromosomal rearrangements, particularly for genes found on wallaby chromosomes X and 3. Expanding comparisons to include chicken and human permitted the putative ancestral marsupial (2n = 14) and therian mammal (2n = 19) karyotypes to be reconstructed.

CONCLUSIONS

Our physical mapping data for the tammar wallaby has uncovered the events shaping marsupial genomes and enabled us to predict the ancestral marsupial karyotype, supporting a 2n = 14 ancestor. Futhermore, our predicted therian ancestral karyotype has helped to understand the evolution of the ancestral eutherian genome.

Mammalian collection on Noah's Ark: the effects of beauty, brain and body size

The importance of today's zoological gardens as the so-called "Noah's Ark" grows as the natural habitat of many species quickly diminishes. Their potential to shelter a large amount of individuals from many species gives us the opportunity to reintroduce a species that disappeared in nature. However, the selection of animals to be kept in zoos worldwide is highly selective and depends on human decisions driven by both ecological criteria such as population size or vulnerability and audience-driven criteria such as aesthetic preferences. Thus we focused our study on the most commonly kept and bred animal class, the mammals, and we asked which factors affect various aspects of the mammalian collection of zoos. We analyzed the presence/absence, population size, and frequency per species of each of the 123 mammalian families kept in the worldwide zoo collection. Our aim was to explain these data using the human-perceived attractiveness of mammalian families, their body weight, relative brain size and species richness of the family. In agreement with various previous studies, we found that the body size and the attractiveness of mammals significantly affect all studied components of the mammalian collection of zoos. There is a higher probability of the large and attractive families to be kept. Once kept, these animals are presented in larger numbers in more zoos. On the contrary, the relative mean brain size only affects the primary selection whether to keep the family or not. It does not affect the zoo population size or the number of zoos that keep the family.

Joint loads in marsupial ankles reflect habitual bipedalism versus quadrupedalism

Joint surfaces of limb bones are loaded in compression by reaction forces generated from body weight and musculotendon complexes bridging them. In general, joints of eutherian mammals have regions of high radiodensity subchondral bone that are better at resisting compressive forces than low radiodensity subchondral bone. Identifying similar form-function relationships between subchondral radiodensity distribution and joint load distribution within the marsupial postcranium, in addition to providing a richer understanding of marsupial functional morphology, can serve as a phylogenetic control in evaluating analogous relationships within eutherian mammals. Where commonalities are established across phylogenetic borders, unifying principles in mammalian physiology, morphology, and behavior can be identified. Here, we assess subchondral radiodensity patterns in distal tibiae of several marsupial taxa characterized by different habitual activities (e.g., locomotion). Computed tomography scanning, maximum intensity projection maps, and pixel counting were used to quantify radiodensity in 41 distal tibiae of bipedal (5 species), arboreal quadrupedal (4 species), and terrestrial quadrupedal (5 species) marsupials. Bipeds (Macropus and Wallabia) exhibit more expansive areas of high radiodensity in the distal tibia than arboreal (Dendrolagus, Phascolarctos, and Trichosurus) or terrestrial quadrupeds (Sarcophilus, Thylacinus, Lasiorhinus, and Vombatus), which may reflect the former carrying body weight only through the hind limbs. Arboreal quadrupeds exhibit smallest areas of high radiodensity, though they differ non-significantly from terrestrial quadrupeds. This could indicate slightly more compliant gaits by arboreal quadrupeds compared to terrestrial quadrupeds. The observed radiodensity patterns in marsupial tibiae, though their statistical differences disappear when controlling for phylogeny, corroborate previously documented patterns in primates and xenarthrans, potentially reflecting inferred limb use during habitual activities such as locomotion. Despite the complex nature of factors contributing to joint loads, broad observance of these patterns across joints and across a variety of taxa suggests that subchondral radiodensity can be used as a unifying form-function principle within Mammalia.

Horizontal transfer of OC1 transposons in the Tasmanian devil

BACKGROUND

There is growing recognition that horizontal DNA transfer, a process known to be common in prokaryotes, is also a significant source of genomic variation in eukaryotes. Horizontal transfer of transposable elements (HTT) may be especially prevalent in eukaryotes given the inherent mobility, widespread occurrence, and prolific abundance of these elements in many eukaryotic genomes.

RESULTS

Here, we provide evidence for a new case of HTT of the transposon family OposCharlie1 (OC1) in the Tasmanian devil, Sarcophilus harrisii. Bioinformatic analyses of OC1 sequences in the Tasmanian devil genome suggest that this transposon infiltrated the common ancestor of the Dasyuridae family ~17 million years ago. This estimate is corroborated by a PCR-based screen for the presence/absence of this family in Tasmanian devils and closely-related species.

CONCLUSIONS

This case of HTT is the first to be reported in dasyurids. It brings the number of animal lineages independently invaded by OC1 to 12, and adds a fourth continent to the pandemic-like pattern of invasion of this transposon. In the context of these data, we discuss the evolutionary history of this transposon family and its potential impact on the diversification of marsupials.

Small mammal investigation in spotted fever focus with DNA-barcoding and taxonomic implications on rodents species from Hainan of China

Although mammals are a well-studied group of animals, making accurate field identification of small mammals is still complex because of morphological variation across developmental stages, color variation of pelages, and often damaged osteological and dental characteristics. In 2008, small mammals were collected for an epidemiological study of a spotted fever outbreak in Hainan, China. Ten species of small mammals were identified by morphological characters in the field, most using pelage color characters only. The study is extended here, in order to assess whether DNA barcoding would be suitable as an identification tool in these small mammals. Barcode clusters showed some incongruence with morphospecies, especially for some species of Rattus and Niviventer, so molecular delineation was carried out with an expanded dataset of combined cytochrome b (Cyt-b) and cytochrome c oxidase subunit I (COI) sequences. COI sequences were successfully amplified from 83% of collected mammals, but failed in all specimens of Suncus murinus, which were thus excluded in DNA barcoding analysis. Of note, ten molecular taxonomic units were found from samples of nine morphologically identified species. Accordingly, 11 species of small mammals were present in the investigated areas, including four Rattus species, three Niviventer species, Callosciurus erythraeus, Neohylomys hainanensis, Tupaia belangeri, and Suncus murinus. Based on the results of the phylogenetic and molecular delineation analyses, the systematic status of some rodent species should be redefined. R. rattus hainanicus and R. rattus sladeni are synonyms of R. andamanensis. R. losea from China and Southeast Asia comprises two independent species: R. losea and R. sakeratensis. Finally, the taxonomic status of three putative species of Niviventer should be further confirmed according to morphological, molecular and ecological characters.

Small mammal investigation in spotted fever focus with DNA-barcoding and taxonomic implications on rodents species from Hainan of China

Although mammals are a well-studied group of animals, making accurate field identification of small mammals is still complex because of morphological variation across developmental stages, color variation of pelages, and often damaged osteological and dental characteristics. In 2008, small mammals were collected for an epidemiological study of a spotted fever outbreak in Hainan, China. Ten species of small mammals were identified by morphological characters in the field, most using pelage color characters only. The study is extended here, in order to assess whether DNA barcoding would be suitable as an identification tool in these small mammals. Barcode clusters showed some incongruence with morphospecies, especially for some species of Rattus and Niviventer, so molecular delineation was carried out with an expanded dataset of combined cytochrome b (Cyt-b) and cytochrome c oxidase subunit I (COI) sequences. COI sequences were successfully amplified from 83% of collected mammals, but failed in all specimens of Suncus murinus, which were thus excluded in DNA barcoding analysis. Of note, ten molecular taxonomic units were found from samples of nine morphologically identified species. Accordingly, 11 species of small mammals were present in the investigated areas, including four Rattus species, three Niviventer species, Callosciurus erythraeus, Neohylomys hainanensis, Tupaia belangeri, and Suncus murinus. Based on the results of the phylogenetic and molecular delineation analyses, the systematic status of some rodent species should be redefined. R. rattus hainanicus and R. rattus sladeni are synonyms of R. andamanensis. R. losea from China and Southeast Asia comprises two independent species: R. losea and R. sakeratensis. Finally, the taxonomic status of three putative species of Niviventer should be further confirmed according to morphological, molecular and ecological characters.

Expansion of CORE-SINEs in the genome of the Tasmanian devil

Background

The genome of the carnivorous marsupial, the Tasmanian devil (Sarcophilus harrisii, Order: Dasyuromorphia), was sequenced in the hopes of finding a cure for or gaining a better understanding of the contagious devil facial tumor disease that is threatening the species’ survival. To better understand the Tasmanian devil genome, we screened it for transposable elements and investigated the dynamics of short interspersed element (SINE) retroposons.

Results

The temporal history of Tasmanian devil SINEs, elucidated using a transposition in transposition analysis, indicates that WSINE1, a CORE-SINE present in around 200,000 copies, is the most recently active element. Moreover, we discovered a new subtype of WSINE1 (WSINE1b) that comprises at least 90% of all Tasmanian devil WSINE1s. The frequencies of WSINE1 subtypes differ in the genomes of two of the other Australian marsupial orders. A co-segregation analysis indicated that at least 66 subfamilies of WSINE1 evolved during the evolution of Dasyuromorphia. Using a substitution rate derived from WSINE1 insertions, the ages of the subfamilies were estimated and correlated with a newly established phylogeny of Dasyuromorphia. Phylogenetic analyses and divergence time estimates of mitochondrial genome data indicate a rapid radiation of the Tasmanian devil and the closest relative the quolls (Dasyurus) around 14 million years ago.

Conclusions

The radiation and abundance of CORE-SINEs in marsupial genomes indicates that they may be a major player in the evolution of marsupials. It is evident that the early phases of evolution of the carnivorous marsupial order Dasyuromorphia was characterized by a burst of SINE activity. A correlation between a speciation event and a major burst of retroposon activity is for the first time shown in a marsupial genome.

A comparative study of mammalian diversification pattern

Although mammals have long been regarded as a successful radiation, the diversification pattern among the clades is still poorly known. Higher-level phylogenies are conflicting and comprehensive comparative analyses are still lacking. Using a recently published supermatrix encompassing nearly all extant mammalian families and a novel comparative likelihood approach (MEDUSA), the diversification pattern of mammalian groups was examined. Both order- and family-level phylogenetic analyses revealed the rapid radiation of Boreoeutheria and Euaustralidelphia in the early mammalian history. The observation of a diversification burst within Boreoeutheria at approximately 100 My supports the Long Fuse model in elucidating placental diversification progress, and the rapid radiation of Euaustralidelphia suggests an important role of biogeographic dispersal events in triggering early Australian marsupial rapid radiation. Diversification analyses based on family-level diversity tree revealed seven additional clades with exceptional diversification rate shifts, six of which represent accelerations in net diversification rate as compared to the background pattern. The shifts gave origin to the clades Muridae+Cricetidae, Bovidae+Moschidae+Cervidae, Simiiformes, Echimyidae, Odontoceti (excluding Physeteridae+Kogiidae+Platanistidae), Macropodidae, and Vespertilionidae. Moderate to high extinction rates from background and boreoeutherian diversification patterns indicate the important role of turnovers in shaping the heterogeneous taxonomic richness observed among extant mammalian groups. Furthermore, the present results emphasize the key role of extinction on erasing unusual diversification signals, and suggest that further studies are needed to clarify the historical radiation of some mammalian groups for which MEDUSA did not detect exceptional diversification rates.