Molecular phylogeny, biogeography, and habitat preference evolution of marsupials

Marsupials exhibit great diversity in ecology and morphology. However, compared to their sister group, the placental mammals, our understanding of many aspects of marsupial evolution remains limited. We use 101 mitochondrial genomes and data from 26 nuclear loci to reconstruct a dated phylogeny including 97% of extant genera and 58% of modern marsupial species. This tree allows us to analyze the evolution of habitat preference and geographic distributions of marsupial species through time. We found a pattern of mesic-adapted lineages evolving to use more arid and open habitats, which is broadly consistent with regional climate and environmental change. However, contrary to the general trend, several lineages subsequently appear to have reverted from drier to more mesic habitats. Biogeographic reconstructions suggest that current views on the connectivity between Australia and New Guinea/Wallacea during the Miocene and Pliocene need to be revised. The antiquity of several endemic New Guinean clades strongly suggests a substantially older period of connection stretching back to the Middle Miocene, and implies that New Guinea was colonized by multiple clades almost immediately after its principal formation.

Evidence that human chlamydia pneumoniae was zoonotically acquired

Zoonotic infections are a growing threat to global health. Chlamydia pneumoniae is a major human pathogen that is widespread in human populations, causing acute respiratory disease, and has been associated with chronic disease. C. pneumoniae was first identified solely in human populations; however, its host range now includes other mammals, marsupials, amphibians, and reptiles. Australian koalas (Phascolarctos cinereus) are widely infected with two species of Chlamydia, C. pecorum and C. pneumoniae. Transmission of C. pneumoniae between animals and humans has not been reported; however, two other chlamydial species, C. psittaci and C. abortus, are known zoonotic pathogens. We have sequenced the 1,241,024-bp chromosome and a 7.5-kb cryptic chlamydial plasmid of the koala strain of C. pneumoniae (LPCoLN) using the whole-genome shotgun method. Comparative genomic analysis, including pseudogene and single-nucleotide polymorphism (SNP) distribution, and phylogenetic analysis of conserved genes and SNPs against the human isolates of C. pneumoniae show that the LPCoLN isolate is basal to human isolates. Thus, we propose based on compelling genomic and phylogenetic evidence that humans were originally infected zoonotically by an animal isolate(s) of C. pneumoniae which adapted to humans primarily through the processes of gene decay and plasmid loss, to the point where the animal reservoir is no longer required for transmission.

Morphologic and genomic characterisation of the koala Chlamydia pneumoniae strain

Chlamydia pneumoniae is a common human and animal pathogen associated with a wide range of upper and lower respiratory tract infections. In more recent years there has been increasing evidence to suggest a link between C. pneumoniae and chronic diseases in humans, including atherosclerosis, stroke and Alzheimer’s disease. C. pneumoniae human strains show little genetic variation, indicating that the human-derived strain originated from a common ancestor in the recent past. Despite extensive information on the genetics and morphology processes of the human strain, knowledge concerning many other hosts (including marsupials, amphibians, reptiles and equines) remains virtually unexplored. The koala (Phascolarctos cinereus) is a native Australian marsupial under threat due to habitat loss, predation and disease. Koalas are very susceptible to chlamydial infections, most commonly affecting the conjunctiva, urogenital tract and/or respiratory tract. To address this gap in the literature, the present study (i) provides a detailed description of the morphologic and genomic architecture of the C. pneumoniae koala (and human) strain, and shows that the koala strain is microscopically, developmentally and genetically distinct from the C. pneumoniae human strain, and (ii) examines the genetic relationship of geographically diverse C. pneumoniae isolates from human, marsupial, amphibian, reptilian and equine hosts, and identifies two distinct lineages that have arisen from animal-to-human cross species transmissions. Chapter One of this thesis explores the scientific problem and aims of this study, while Chapter Two provides a detailed literature review of the background in this field of work. Chapter Three, the first results chapter, describes the morphology and developmental stages of C. pneumoniae koala isolate LPCoLN, as revealed by fluorescence and transmission electron microscopy. The profile of this isolate, when cultured in HEp-2 human epithelial cells, was quite different to the human AR39 isolate. Koala LPCoLN inclusions were larger; the elementary bodies did not have the characteristic pear-shaped appearance, and the developmental cycle was completed within a shorter period of time (as confirmed by quantitative real-time PCR). These in vitro findings might reflect biological differences between koala LPCoLN and human AR39 in vivo. Chapter Four describes the complete genome sequence of the koala respiratory pathogen, C. pneumoniae LPCoLN. This is the first animal isolate of C. pneumoniae to be fully-sequenced. The genome sequence provides new insights into genomic ‘plasticity’ (organisation), evolution and biology of koala LPCoLN, relative to four complete C. pneumoniae human genomes (AR39, CWL029, J138 and TW183). Koala LPCoLN contains a plasmid that is not shared with any of the human isolates, there is evidence of gene loss in nucleotide salvage pathways, and there are 10 hot spot genomic regions of variation that were previously not identified in the C. pneumoniae human genomes. Sequence (partial-length) from a second, independent, wild koala isolate (EBB) at several gene loci confirmed that the koala LPCoLN isolate was representative of a koala C. pneumoniae strain. The combined sequence data provides evidence that the C. pneumoniae animal (koala LPCoLN) genome is ancestral to the C. pneumoniae human genomes and that human infections may have originated from zoonotic infections. Chapter Five examines key genome components of the five C. pneumoniae genomes in more detail. This analysis reveals genomic features that are shared by and/or contribute to the broad ecological adaptability and evolution of C. pneumoniae. This analysis resulted in the identification of 65 gene sequences for further analysis of intraspecific variation, and revealed some interesting differences, including fragmentation, truncation and gene decay (loss of redundant ancestral traits). This study provides valuable insights into metabolic diversity, adaptation and evolution of C. pneumoniae. Chapter Six utilises a subset of 23 target genes identified from the previous genomic comparisons and makes a significant contribution to our understanding of genetic variability among C. pneumoniae human (11) and animal (6 amphibian, 5 reptilian, 1 equine and 7 marsupial hosts) isolates. It has been shown that the animal isolates are genetically diverse, unlike the human isolates that are virtually clonal. More convincing evidence that C. pneumoniae originated in animals and recently (in the last few hundred thousand years) crossed host species to infect humans is provided in this study. It is proposed that two animal-to-human cross species events have occurred in the context of the results, one evident by the nearly clonal human genotype circulating in the world today, and the other by a more animal-like genotype apparent in Indigenous Australians. Taken together, these data indicate that the C. pneumoniae koala LPCoLN isolate has morphologic and genomic characteristics that are distinct from the human isolates. These differences may affect the survival and activity of the C. pneumoniae koala pathogen in its natural host, in vivo. This study, by utilising the genetic diversity of C. pneumoniae, identified new genetic markers for distinguishing human and animal isolates. However, not all C. pneumoniae isolates were genetically diverse; in fact, several isolates were highly conserved, if not identical in sequence (i.e. Australian marsupials) emphasising that at some stage in the evolution of this pathogen, there has been an adaptation/s to a particular host, providing some stability in the genome. The outcomes of this study by experimental and bioinformatic approaches have significantly enhanced our knowledge of the biology of this pathogen and will advance opportunities for the investigation of novel vaccine targets, antimicrobial therapy, or blocking of pathogenic pathways.

Chlamydia pneumoniae is genetically diverse in animals and appears to have crossed the host barrier to humans on (at least) two occasions

Chlamydia pneumoniae is a common human and animal pathogen associated with a wide range of diseases. Since the first isolation of C. pneumoniae TWAR in 1965, all human isolates have been essentially clonal, providing little evolutionary insight. To address this gap, we investigated the genetic diversity of 30 isolates from diverse geographical locations, from both human and animal origin (amphibian, reptilian, equine and marsupial). Based on the level of variation that we observed at 23 discreet gene loci, it was clearly evident that the animal isolates were more diverse than the isolates of human origin. Furthermore, we show that C. pneumoniae isolates could be grouped into five major genotypes, A-E, with A, B, D and E genotypes linked by geographical location, whereas genotype C was found across multiple continents. Our evidence strongly supports two separate animal-to-human cross species transfer events in the evolutionary history of this pathogen. The C. pneumoniae human genotype identified in the USA, Canada, Taiwan, Iran, Japan, Korea and Australia (non- Indigenous) most likely originated from a single amphibian or reptilian lineage, which appears to have been previously geographically widespread. We identified a separate human lineage present in two Australian Indigenous isolates (independent geographical locations). This lineage is distinct and is present in Australian amphibians as well as a range of Australian marsupials.

Mammalian phylogeny

Living mammals can be divided into three subclasses (monotremes, marsupials and placentals) and within these, about 27 orders. Final resolution of the relationships between the orders is only now being achieved with the increased availability of deoxyribonucleic acid (DNA) sequences. Highlights include the deep division of placental mammals into African (Afrotheria), South American (Xenarthra) and northern hemisphere (Boreoeutheria) super-orders, and the finding that the once considered primitive ‘Insectivora’ and ‘Edentata’ clades, in fact, have members distributed widely among these super-orders. Another surprise finding from DNA studies has been that whale origins lie among the even-toed ungulates (Artiodactyla). Our order, Primates is most closely related to the flying lemurs and next, the tree shrews. With the mammal phylogeny becoming well resolved, it is increasingly being used as a framework for inferring evolutionary and ecological processes, such as adaptive radiation.

Combined mitochondrial and nuclear DNA sequences resolve the interrelations of the major Australasian marsupial radiations

Australasian marsupials include three major radiations, the insectivorous/carnivorous Dasyuromorphia, the omnivorous bandicoots (Peramelemorphia), and the largely herbivorous diprotodontians. Morphologists have generally considered the bandicoots and diprotodontians to be closely related, most prominently because they are both syndactylous (with the 2nd and 3rd pedal digits being fused). Molecular studies have been unable to confirm or reject this Syndactyla hypothesis. Here we present new mitochondrial (mt) genomes from a spiny bandicoot (Echymipera rufescens) and two dasyurids, a fat-tailed dunnart (Sminthopsis crassicaudata) and a northern quoll (Dasyurus hallucatus). By comparing trees derived from pairwise base-frequency differences between taxa with standard (absolute, uncorrected) distance trees, we infer that composition bias among mt protein-coding and RNA sequences is sufficient to mislead tree reconstruction. This can explain incongruence between trees obtained from mt and nuclear data sets. However, after excluding major sources of compositional heterogeneity, both the “reduced-bias” mt and nuclear data sets clearly favor a bandicoot plus dasyuromorphian association, as well as a grouping of kangaroos and possums (Phalangeriformes) among diprotodontians. Notably, alternatives to these groupings could only be confidently rejected by combining the mt and nuclear data. Elsewhere on the tree, Dromiciops appears to be sister to the monophyletic Australasian marsupials, whereas the placement of the marsupial mole (Notoryctes) remains problematic. More generally, we contend that it is desirable to combine mt genome and nuclear sequences for inferring vertebrate phylogeny, but as separately modeled process partitions. This strategy depends on detecting and excluding (or accounting for) major sources of nonhistorical signal, such as from compositional nonstationarity.

Molecules, morphology, and ecology indicate a recent, amphibious ancestry for echidnas

The semiaquatic platypus and terrestrial echidnas (spiny anteaters) are the only living egg-laying mammals (monotremes). The fossil record has provided few clues as to their origins and the evolution of their ecological specializations; however, recent reassignment of the Early Cretaceous Teinolophos and Steropodon to the platypus lineage implies that platypuses and echidnas diverged >112.5 million years ago, reinforcing the notion of monotremes as living fossils. This placement is based primarily on characters related to a single feature, the enlarged mandibular canal, which supplies blood vessels and dense electrosensory receptors to the platypus bill. Our reevaluation of the morphological data instead groups platypus and echidnas to the exclusion of Teinolophos and Steropodon and suggests that an enlarged mandibular canal is ancestral for monotremes (partly reversed in echidnas, in association with general mandibular reduction). A multigene evaluation of the echidna–platypus divergence using both a relaxed molecular clock and direct fossil calibrations reveals a recent split of 19–48 million years ago. Platypus-like monotremes (Monotrematum) predate this divergence, indicating that echidnas had aquatically foraging ancestors that reinvaded terrestrial ecosystems. This ecological shift and the associated radiation of echidnas represent a recent expansion of niche space despite potential competition from marsupials. Monotremes might have survived the invasion of marsupials into Australasia by exploiting ecological niches in which marsupials are restricted by their reproductive mode. Morphology, ecology, and molecular biology together indicate that Teinolophos and Steropodon are basal monotremes rather than platypus relatives, and that living monotremes are a relatively recent radiation.

Inferring kangaroo phylogeny from incongruent nuclear and mitochondrial genes

The marsupial genus Macropus includes three subgenera, the familiar large grazing kangaroos and wallaroos of M. (Macropus) and M. (Osphranter), as well as the smaller mixed grazing/browsing wallabies of M. (Notamacropus). A recent study of five concatenated nuclear genes recommended subsuming the predominantly browsing Wallabia bicolor (swamp wallaby) into Macropus. To further examine this proposal we sequenced partial mitochondrial genomes for kangaroos and wallabies. These sequences strongly favour the morphological placement of W. bicolor as sister to Macropus, although place M. irma (black-gloved wallaby) within M. (Osphranter) rather than as expected, with M. (Notamacropus). Species tree estimation from separately analysed mitochondrial and nuclear genes favours retaining Macropus and Wallabia as separate genera. A simulation study finds that incomplete lineage sorting among nuclear genes is a plausible explanation for incongruence with the mitochondrial placement of W. bicolor, while mitochondrial introgression from a wallaroo into M. irma is the deepest such event identified in marsupials. Similar such coalescent simulations for interpreting gene tree conflicts will increase in both relevance and statistical power as species-level phylogenetics enters the genomic age. Ecological considerations in turn, hint at a role for selection in accelerating the fixation of introgressed or incompletely sorted loci. More generally the inclusion of the mitochondrial sequences substantially enhanced phylogenetic resolution. However, we caution that the evolutionary dynamics that enhance mitochondria as speciation indicators in the presence of incomplete lineage sorting may also render them especially susceptible to introgression.

Taxonomy and palaeoecology of Quaternary faunas from caves in eastern tropical Queensland : a record of broad-scale environmental change

Analysis of fossils from cave deposits at Mount Etna (eastern-central Queensland) has established that a species-rich rainforest palaeoenvironment existed in that area during the middle Pleistocene. This unexpected finding has implications for several fields (e.g., biogeography/phylogeography of rainforest-adapted taxa, and the impact of climate change on rainforest communities), but it was unknown whether the Mount Etna sites represented a small refugial patch of rainforest or was more widespread. In this study numerous bone deposits in caves in north-east Queensland are analysed to reconstruct the environmental history of the area during the late Quaternary. Study sites are in the Chillagoe/Mitchell Palmer and Broken River/Christmas Creek areas. The cave fossil records in these study areas are compared with dated (middle Pleistocene-Holocene) cave sites in the Mount Etna area. Substantial taxonomic work on the Mount Etna faunas (particularly dasyurid marsupials and murine rodents) is also presented as a prerequisite for meaningful comparison with the study sites further north. Middle Pleistocene sites at Mount Etna contain species indicative of a rainforest palaeoenvironment. Small mammal assemblages in the Mount Etna rainforest sites (>500-280 ka) are unexpectedly diverse and composed almost entirely of new species. Included in the rainforest assemblages are lineages with no extant representatives in rainforest (e.g., Leggadina), one genus previously known only from New Guinea (Abeomelomys), and forms that appear to bridge gaps between related but morphologically-divergent extant taxa ('B-rat' and 'Pseudomys C'). Curiously, some taxa (e.g., Melomys spp.) are notable for their absence from the Mount Etna rainforest sites. After 280 ka the rainforest faunas are replaced by species adapted to open, dry habitats. At that time the extinct ‘rainforest’ dasyurids and rodents are replaced by species that are either extant or recently extant. By the late Pleistocene all ‘rainforest’ and several ‘dry’ taxa are locally or completely extinct, and the small mammal fauna resembles that found in the area today. The faunal/environmental changes recorded in the Mount Etna sites were interpreted by previous workers as the result of shifts in climate during the Pleistocene. Many samples from caves in the Chillagoe/Mitchell-Palmer and Broken River/Christmas Creek areas are held in the Queensland Museum’s collection. These, supplemented with additional samples collected in the field as well as samples supplied by other workers, were systematically and palaeoecologically analysed for the first time. Palaeoecological interpretation of the faunal assemblages in the sites suggests that they encompass a similar array of palaeoenvironments as the Mount Etna sites. ‘Rainforest’ sites at the Broken River are here interpreted as being of similar age to those at Mount Etna, suggesting the possibility of extensive rainforest coverage in eastern tropical Queensland during part of the Pleistocene. Likewise, faunas suggesting open, dry palaeoenvironments are found at Chillagoe, the Broken River and Mount Etna, and may be of similar age. The 'dry' faunal assemblage at Mount Etna (Elephant hole Cave) dates to 205-170 ka. Dating of one of the Chillagoe sites (QML1067) produced a maximum age for the deposit of approximately 200 ka, and the site is interpreted as being close to that age, supporting the interpretation of roughly contemporaneous deposition at Mount Etna and Chillagoe. Finally, study sites interpreted as being of late Pleistocene-Holocene age show faunal similarities to sites of that age near Mount Etna. This study has several important implications for the biogeography and phylogeography of murine rodents, and represents a major advance in the study of the Australian murine fossil record. Likewise the survey of the northern study areas is the first systematic analysis of multiple sites in those areas, and is thus a major contribution to knowledge of tropical Australian faunas during the Quaternary. This analysis suggests that climatic changes during the Pleistocene affected a large area of eastern tropical Queensland in similar ways. Further fieldwork and dating is required to properly analyse the geographical extent and timing of faunal change in eastern tropical Queensland.

Molecular characterisation and expression analysis of interferon gamma in response to natural chlamydia infection in the koala, phascolarctos cinereus

Interferon gamma (IFNγ) is a key Th1 cytokine, with a principal role in the immune response against intracellular organisms such as Chlamydia. Along with being responsible for significant morbidity in human populations, Chlamydia is also responsible for wide spread infection and disease in many animal hosts, with reports that many Australian koala subpopulations are endemically infected. An understanding of the role played by IFNγ in koala chlamydial diseases is important for the establishment of better prophylactic and therapeutic approaches against chlamydial infection in this host. A limited number of IFNγ sequences have been published from marsupials and no immune reagents to measure expression have been developed. Through preliminary analysis of the koala transcriptome, we have identified the full coding sequence of the koala IFNγ gene. Transcripts were identified in spleen and lymph node tissue samples. Phylogenetic analysis demonstrated that koala IFNγ is closely related to other marsupial IFNγ sequences and more distantly related to eutherian mammals. To begin to characterise the role of this important cytokine in the koala's response to chlamydial infection, we developed a quantitative real time PCR assay and applied it to a small cohort of koalas with and without active chlamydial disease, revealing significant differences in expression patterns between the groups. Description of the IFNγ sequence from the koala will not only assist in understanding this species' response to its most important pathogen but will also provide further insight into the evolution of the marsupial immune system

Muramidases found in the foregut microbiome of the Tammar wallaby can direct cell aggregation and biofilm formation

We describe here the role of muramidases present in clones of metagenomic DNA that result in cell aggregation and biofilm formation by Escherichia coli. The metagenomic clones were obtained from uncultured Lachnospiraceae-affiliated bacteria resident in the foregut microbiome of the Tammar wallaby. One of these fosmid clones (p49C2) was chosen for more detailed studies and a variety of genetic methods were used to delimit the region responsible for the phenotype to an open reading frame of 1425 bp. Comparative sequence analysis with other fosmid clones giving rise to the same phenotype revealed the presence of muramidase homologues with the same modular composition. Phylogenetic analysis of the fosmid sequence data assigned these fosmid inserts to recently identified, but uncultured, phylogroups of Lachnospiraceae believed to be numerically dominant in the foregut microbiome of the Tammar wallaby. The muramidase is a modular protein containing putative N-acetylmuramoyl--alanine amidase and an endo-β-N-acetylglucosaminidase catalytic module, with a similar organization and functional properties to some Staphylococcal autolysins that also confer adhesive properties and biofilm formation. We also show here that the cloned muramidases result in the production of extracellular DNA, which appears to be the key for biofilm formation and autoaggregation. Collectively, these findings suggest that biofilm formation and cell aggregation in gut microbiomes might occur via the concerted action of carbohydrate-active enzymes and the production of extracellular DNA to serve as a biofilm scaffold.

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

Practical guidelines for planning and implementing fauna translocations

A successful translocation involves many complex factors, including a genetically appropriate source population that can sustain harvest, social and governmental support, assessment of disease transmission risk and a release site with appropriately secure habitat that can support population establishment and persistance. This information is typically discussed during staturory approval processes and can take considerable time. However, following approval, for translocations of most fauna, the initial critical step involves the inherently stressful process of capture, holding, transportation and release. This process is unpredictable and novel, and is especially challenging for wild animals when they are confined in close proximity to conspecifics and humans. In contrast, captive-reared animals have to cope with the unfamiliar challenges of finding food and shelter, along with coping with competition and predation. Little has been written in the scientific literature about the translocation process. This is unsurprising because this process has usually been the realm of skilled practioners, often with animal husbandry backgrounds, rather than research scientists. Highly skilled intuition, observation and the translocation practioner's equivalent of a 'green thumb' often guides the way. However, theory and experimentation, particularly on the effects of stress, is available and this work is invaluable for a successful translocation. Here, we provide a brief description of the translocation process, and discussion of what stress is and how it can be managed. We then provide practical guidelines for the successful translocation of invertebrates, lizards, turtles, passerine birds, marsupials and bats, using examples from Australia and New Zealand.

Phylogenetic analysis of human Chlamydia pneumoniae strains reveals a distinct Australian indigenous clade that predates European exploration of the continent

Background The obligate intracellular bacterium Chlamydia pneumoniae is a common respiratory pathogen, which has been found in a range of hosts including humans, marsupials and amphibians. Whole genome comparisons of human C. pneumoniae have previously highlighted a highly conserved nucleotide sequence, with minor but key polymorphisms and additional coding capacity when human and animal strains are compared. Results In this study, we sequenced three Australian human C. pneumoniae strains, two of which were isolated from patients in remote indigenous communities, and compared them to all available C. pneumoniae genomes. Our study demonstrated a phylogenetically distinct human C. pneumoniae clade containing the two indigenous Australian strains, with estimates that the most recent common ancestor of these strains predates the arrival of European settlers to Australia. We describe several polymorphisms characteristic to these strains, some of which are similar in sequence to animal C. pneumoniae strains, as well as evidence to suggest that several recombination events have shaped these distinct strains. Conclusions Our study reveals a greater sequence diversity amongst both human and animal C. pneumoniae strains, and suggests that a wider range of strains may be circulating in the human population than current sampling indicates.