Genome of the marsupial Monodelphis domestica reveals innovation in non-coding sequences

A three-dimensional stereotaxic atlas of the gray short-tailed opossum (Monodelphis domestica) brain

The gray short-tailed opossum (Monodelphis domestica) is a small marsupial gaining recognition as a laboratory animal in biomedical research. Despite numerous studies on opossum neuroanatomy, a consistent and comprehensive neuroanatomical reference for this species is still missing. Here we present the first three-dimensional, multimodal atlas of the Monodelphis opossum brain. It is based on four complementary imaging modalities: high resolution ex vivo magnetic resonance images, micro-computed tomography scans of the cranium, images of the face of the cutting block, and series of sections stained with the Nissl method and for myelinated fibers. Individual imaging modalities were reconstructed into a three-dimensional form and then registered to the MR image by means of affine and deformable registration routines. Based on a superimposition of the 3D images, 113 anatomical structures were demarcated and the volumes of individual regions were measured. The stereotaxic coordinate system was defined using a set of cranial landmarks: interaural line, bregma, and lambda, which allows for easy expression of any location within the brain with respect to the skull. The atlas is released under the Creative Commons license and available through various digital atlasing web services.

Resolving kangaroo phylogeny and overcoming retrotransposon ascertainment bias

Reconstructing phylogeny from retrotransposon insertions is often limited by access to only a single reference genome, whereby support for clades that do not include the reference taxon cannot be directly observed. Here we have developed a new statistical framework that accounts for this ascertainment bias, allowing us to employ phylogenetically powerful retrotransposon markers to explore the radiation of the largest living marsupials, the kangaroos and wallabies of the genera Macropus and Wallabia. An exhaustive in silico screening of the tammar wallaby (Macropus eugenii) reference genome followed by experimental screening revealed 29 phylogenetically informative retrotransposon markers belonging to a family of endogenous retroviruses. We identified robust support for the enigmatic swamp wallaby (Wallabia bicolor) falling within a paraphyletic genus, Macropus. Our statistical approach provides a means to test for incomplete lineage sorting and introgression/hybridization in the presence of the ascertainment bias. Using retrotransposons as “molecular fossils”, we reveal one of the most complex patterns of hemiplasy yet identified, during the rapid diversification of kangaroos and wallabies. Ancestral state reconstruction incorporating the new retrotransposon phylogenetic information reveals multiple independent ecological shifts among kangaroos into more open habitats, coinciding with the Pliocene onset of increased aridification in Australia from ~3.6 million years ago.

Marsupial and monotreme cathelicidins display antimicrobial activity, including against methicillin-resistant Staphylococcus aureus

With the growing demand for new antibiotics to combat increasing multi-drug resistance, a family of antimicrobial peptides known as cathelicidins has emerged as potential candidates. Expansions in cathelicidin-encoding genes in marsupials and monotremes are of specific interest as the peptides they encode have evolved to protect immunologically naive young in the harsh conditions of the pouch and burrow. Our previous work demonstrated that some marsupial and monotreme cathelicidins have broad-spectrum antibacterial activity and kill resistant bacteria, but the activity of many cathelicidins is unknown. To investigate associations between peptide antimicrobial activity and physiochemical properties, we tested 15 cathelicidin mature peptides from tammar wallaby, grey short-tailed opossum, platypus and echidna for antimicrobial activity against a range of bacterial and fungal clinical isolates. One opossum cathelicidin ModoCath4, tammar wallaby MaeuCath7 and echidna Taac-CATH1 had broad-spectrum antibacterial activity and killed methicillin-resistant Staphylococcus aureus. However, antimicrobial activity was reduced in the presence of serum or whole blood, and non-specific toxicity was observed at high concentrations. The active peptides were highly charged, potentially increasing binding to microbial surfaces, and contained amphipathic helical structures, which may facilitate membrane permeabilisation. Peptide sequence homology, net charge, amphipathicity and alpha helical content did not correlate with antimicrobial activity. However active peptides contained a significantly higher percentage of cationic residues than inactive ones, which may be used to predict active peptides in future work. Along with previous studies, our results indicate that marsupial and monotreme cathelicidins show potential for development as novel therapeutics to combat increasing antimicrobial resistance.

Satellite DNA: An Evolving Topic

Satellite DNA represents one of the most fascinating parts of the repetitive fraction of the eukaryotic genome. Since the discovery of highly repetitive tandem DNA in the 1960s, a lot of literature has extensively covered various topics related to the structure, organization, function, and evolution of such sequences. Today, with the advent of genomic tools, the study of satellite DNA has regained a great interest. Thus, Next-Generation Sequencing (NGS), together with high-throughput in silico analysis of the information contained in NGS reads, has revolutionized the analysis of the repetitive fraction of the eukaryotic genomes. The whole of the historical and current approaches to the topic gives us a broad view of the function and evolution of satellite DNA and its role in chromosomal evolution. Currently, we have extensive information on the molecular, chromosomal, biological, and population factors that affect the evolutionary fate of satellite DNA, knowledge that gives rise to a series of hypotheses that get on well with each other about the origin, spreading, and evolution of satellite DNA. In this paper, I review these hypotheses from a methodological, conceptual, and historical perspective and frame them in the context of chromosomal organization and evolution.

The transcriptome of the Didelphis virginiana opossum kidney OK proximal tubule cell line

The OK cell line derived from the kidney of a female opossum Didelphis virginiana has proven to be a useful model in which to investigate the unique regulation of ion transport and membrane trafficking mechanisms in the proximal tubule (PT). Sequence data and comparison of the transcriptome of this cell line to eutherian mammal PTs would further broaden the utility of this culture model. However, the genomic sequence for D. virginiana is not available and although a draft genome sequence for the opossum Monodelphis domestica (sequenced in 2012 by the Broad Institute) exists, transcripts sequenced from both species show significant divergence. The M. domestica sequence is not highly annotated, and the majority of transcripts are predicted rather than experimentally validated. Using deep RNA sequencing of the D. virginiana OK cell line, we characterized its transcriptome via de novo transcriptome assembly and alignment to the M. domestica genome. The quality of the de novo assembled transcriptome was assessed by the extent of homology to sequences in nucleotide and protein databases. Gene expression levels in the OK cell line, from both the de novo transcriptome and genes aligned to the M. domestica genome, were compared with publicly available rat kidney nephron segment expression data. Our studies demonstrate the expression in OK cells of numerous PT-specific ion transporters and other key proteins relevant for rodent and human PT function. Additionally, the sequence and expression data reported here provide an important resource for genetic manipulation and other studies on PT cell function using these cells.

Ontogeny of myosin isoform expression and prehensile function in the tail of the gray short-tailed opossum (Monodelphis domestica)

Terrestrial opossums use their semiprehensile tail for grasping nesting materials as opposed to arboreal maneuvering. We relate the development of this adaptive behavior with ontogenetic changes in myosin heavy chain (MHC) isoform expression from 21 days to adulthood. Monodelphis domestica is expected to demonstrate a progressive ability to flex the distal tail up to age 7 mo, when it should exhibit routine nest construction. We hypothesize that juvenile stages (3-7 mo) will be characterized by retention of the neonatal isoform (MHC-Neo), along with predominant expression of fast MHC-2X and -2B, which will transition into greater MHC-1β and -2A isoform content as development progresses. This hypothesis was tested using Q-PCR to quantify and compare gene expression of each isoform with its protein content determined by gel electrophoresis and densitometry. These data were correlated with nesting activity in an age-matched sample of each age group studied. Shifts in regulation of MHC gene transcripts matched well with isoform expression. Notably, mRNA for MHC-Neo and -2B decrease, resulting in little-to-no isoform translation after age 7 mo, whereas mRNA for MHC-1β and -2A increase, and this corresponds with subtle increases in content for these isoforms into late adulthood. Despite the tail remaining intrinsically fast-contracting, a critical growth period for isoform transition is observed between 7 and 13 mo, correlating primarily with use of the tail during nesting activities. Functional transitions in MHC isoforms and fiber type properties may be associated with muscle "tuning" repetitive nest remodeling tasks requiring sustained contractions of the caudal flexors.NEW & NOTEWORTHY Little is understood about skeletal muscle development as it pertains to tail prehensility in mammals. This study uses an integrative approach of relating both MHC gene and protein expression with behavioral and morphometric changes to reveal a predominant fast MHC expression with subtle isoform transitions in caudal muscle across ontogeny. The functional shifts observed are most notably correlated with increased tail grasping for nesting activities.

Platypus and opossum calcitonins exhibit strong activities, even though they belong to mammals

In mammalian assay systems, calcitonin peptides of non-mammalian species exhibit stronger activity than those of mammals. Recently, comparative analyses of a wide-range of species revealed that platypus and opossum, which diverged early from other mammals, possess calcitonins that are more similar in amino acid sequence to those of non-mammals than mammals. We herein determined whether platypus and opossum calcitonins exhibit similar biological activities to those of non-mammalian calcitonins using an assay of actin ring formation in mouse osteoclasts. We also compared the dose-dependent effects of each calcitonin on cAMP production in osteoclasts. Consistent with the strong similarities in their primary amino acid sequences, platypus and opossum calcitonins disrupted actin rings with similar efficacies to that of salmon calcitonin. Human calcitonin exhibited the weakest inhibitory potency and required a 100-fold higher concentration (EC₅₀=3×10^-11M) than that of salmon calcitonin (EC₅₀=2×10^-13M). Platypus and opossum calcitonins also induced cAMP production in osteoclast cultures with the same efficacies as that of salmon calcitonin. Thus, platypus and opossum calcitonins exhibited strong biological activities, similar to those of the salmon. In addition, phylogenetic analysis revealed that platypus and opossum calcitonins clustered with the salmon-type group but not human- or porcine-type group. These results suggest that platypus and opossum calcitonins are classified into the salmon-type group, in terms of the biological activities and amino acid sequences.

Opossum APOBEC1 is a DNA mutator with retrovirus and retroelement restriction activity

APOBEC3s (A3s) are single-stranded DNA cytosine deaminases that provide innate immune defences against retroviruses and mobile elements. A3s are specific to eutherian mammals because no direct homologs exist at the syntenic genomic locus in metatherian (marsupial) or prototherian (monotreme) mammals. However, the A3s in these species have the likely evolutionary precursors, the antibody gene deaminase AID and the RNA/DNA editing enzyme APOBEC1 (A1). Here, we used cell culture-based assays to determine whether opossum A1 restricts the infectivity of retroviruses including human immunodeficiency virus type 1 (HIV-1) and the mobility of LTR/non-LTR retrotransposons. Opossum A1 partially inhibited HIV-1, as well as simian immunodeficiency virus (SIV), murine leukemia virus (MLV), and the retrotransposon MusD. The mechanism of inhibition required catalytic activity, except for human LINE1 (L1) restriction, which was deamination-independent. These results indicate that opossum A1 functions as an innate barrier to infection by retroviruses such as HIV-1, and controls LTR/non-LTR retrotransposition in marsupials.

The Unique Paired Retinal Vessels of the Gray Short-Tailed Opossum (Monodelphis domestica) and Their Relationship to Astrocytes and Microglial Cells

In marsupials that possess a retinal vasculature, the arterial and venous segments, down to the smallest calibre capillaries, have been shown to occur in pairs. This pattern is seen in the marsupial central nervous system (CNS) but not in other tissues in this group or in any tissues in eutherian mammals. The present study aimed to determine if the gray short-tailed opossum (Monodelphis domestica), a south American marsupial, possesses double retinal vessels. Secondly, we investigated the relationship between vessels and astrocytes and microglia, which are known to play pivotal roles in the blood retinal barrier and immune surveillance respectively. Eyes from M. domestica between 2 months and 33 months of age were examined by bright field and fluorescein angiography, resin histology, and wholemount immunostaining. Retinal vessels in this marsupial always occur in closely related pairs with the arteriolar limb usually on the vitread aspect. Branches penetrate the retina to form layers of paired capillaries as far as the outer nuclear layer. Dense networks of GFAP⁺ astrocytes enveloped the vitread aspect of vessels. No particularly strong association with blood vessels and ramified Iba1⁺ and Ib4⁺ microglia was noted. M. domestica possessed the unusual paired vasculature and capillary loops arrangement previously described in the marsupial CNS. These observations in a small laboratory-friendly marsupial open up new frontiers to investigate the factors that regulate paired blood vessel development and the functional significance of this arrangement when compared to the anastomotic pattern observed in the retina of eutherian mammals. Anat Rec, 300:1391-1400, 2017. © 2017 Wiley Periodicals, Inc.

Implications of monotreme and marsupial chromosome evolution on sex determination and differentiation

Studies of chromosomes from monotremes and marsupials endemic to Australasia have provided important insight into the evolution of their genomes as well as uncovering fundamental differences in their sex determination/differentiation pathways. Great advances have been made this century into solving the mystery of the complicated sex chromosome system in monotremes. Monotremes possess multiple different X and Y chromosomes and a candidate sex determining gene has been identified. Even greater advancements have been made for marsupials, with reconstruction of the ancestral karyotype enabling the evolutionary history of marsupial chromosomes to be determined. Furthermore, the study of sex chromosomes in intersex marsupials has afforded insight into differences in the sexual differentiation pathway between marsupials and eutherians, together with experiments showing the insensitivity of the mammary glands, pouch and scrotum to exogenous hormones, led to the hypothesis that there is a gene (or genes) on the X chromosome responsible for the development of either pouch or scrotum. This review highlights the major advancements made towards understanding chromosome evolution and how this has impacted on our understanding of sex determination and differentiation in these interesting mammals.

Variants in the host genome may inhibit tumour growth in devil facial tumours: evidence from genome-wide association

Devil facial tumour disease (DFTD) has decimated wild populations of Tasmanian devils (Sarcophilus harrisii) due to its ability to avoid immune detection and pass from host to host by biting. A small number of devils have been observed to spontaneously recover from the disease which is otherwise fatal. We have sequenced the genomes of these rare cases and compared them to the genomes of devils who succumbed to the disease. Genome-wide association, based on this limited sampling, highlighted two key genomic regions potentially associated with ability to survive DFTD. Following targeted genotyping in additional samples, both of these loci remain significantly different between cases and controls, with the PAX3 locus retaining significance at the 0.001 level, though genome-wide significance was not achieved. We propose that PAX3 may be involved in a regulatory pathway that influences the slowing of tumour growth and may allow more time for an immune response to be mounted in animals with regressed tumours. This provides an intriguing hypothesis for further research and could provide a novel route of treatment for this devastating disease.

FGF2, FGF3 and FGF4 expression pattern during molars odontogenesis in Didelphis albiventris

Odontogenesis is guided by a complex signaling cascade in which several molecules, including FGF2-4, ensure all dental groups development and specificity. Most of the data on odontogenesis derives from rodents, which does not have all dental groups. Didelphis albiventris is an opossum with the closest dentition to humans, and the main odontogenesis stages occur when the newborns are in the pouch. In this study, D. albiventris postnatals were used to characterize the main stages of their molars development; and also to establish FGF2, FGF3 and FGF4 expression pattern. D. albiventris postnatals were processed for histological and indirect immunoperoxidase analysis of the tooth germs. Our results revealed similar dental structures between D. albiventris and mice. However, FGF2, FGF3 and FGF4 expression patterns were observed in a larger number of dental structures, suggesting broader functions for these molecules in this opossum species. The knowledge of the signaling that determinates odontogenesis in an animal model with complete dentition may contribute to the development of therapies for the replacement of lost teeth in humans. This study may also contribute to the implementation of D. albiventris as model for Developmental Biology studies.

Diversity as Opportunity: Insights from 600 Million Years of AHR Evolution

The aryl hydrocarbon receptor (AHR) was for many years of interest only to pharmacologists and toxicologists. However, this protein has fundamental roles in biology that are being revealed through studies in diverse animal species. The AHR is an ancient protein. AHR homologs exist in most major groups of modern bilaterian animals, including deuterostomes (chordates, hemichordates, echinoderms) and the two major clades of protostome invertebrates [ecdysozoans (e.g. arthropods and nematodes) and lophotrochozoans (e.g. molluscs and annelids)]. AHR homologs also have been identified in cnidarians such as the sea anemone Nematostella and in the genome of Trichoplax, a placozoan. Bilaterians, cnidarians, and placozoans form the clade Eumetazoa, whose last common ancestor lived approximately 600 million years ago (MYA). The presence of AHR homologs in modern representatives of all these groups indicates that the original eumetazoan animal possessed an AHR homolog. Studies in invertebrates and vertebrates reveal parallel functions of AHR in the development and function of sensory neural systems, suggesting that these may be ancestral roles. Vertebrate animals are characterized by the expansion and diversification of AHRs, via gene and genome duplications, from the ancestral protoAHR into at least five classes of AHR-like proteins: AHR, AHR1, AHR2, AHR3, and AHRR. The evolution of multiple AHRs in vertebrates coincided with the acquisition of high-affinity binding of halogenated and polynuclear aromatic hydrocarbons and the emergence of adaptive functions involving regulation of xenobiotic-metabolizing enzymes and roles in adaptive immunity. The existence of multiple AHRs may have facilitated subfunction partitioning and specialization of specific AHR types in some taxa. Additional research in diverse model and non-model species will continue to enrich our understanding of AHR and its pleiotropic roles in biology and toxicology.

A new developmental mechanism for the separation of the mammalian middle ear ossicles from the jaw

Multiple mammalian lineages independently evolved a definitive mammalian middle ear (DMME) through breakdown of Meckel's cartilage (MC). However, the cellular and molecular drivers of this evolutionary transition remain unknown for most mammal groups. Here, we identify such drivers in the living marsupial opossum Monodelphis domestica, whose MC transformation during development anatomically mirrors the evolutionary transformation observed in fossils. Specifically, we link increases in cellular apoptosis and TGF-BR2 signalling to MC breakdown in opossums. We demonstrate that a simple change in TGF-β signalling is sufficient to inhibit MC breakdown during opossum development, indicating that changes in TGF-β signalling might be key during mammalian evolution. Furthermore, the apoptosis that we observe during opossum MC breakdown does not seemingly occur in mouse, consistent with homoplastic DMME evolution in the marsupial and placental lineages.

A new mammalian model system for thalidomide teratogenesis: Monodelphis domestica

From 1957 to 1962, thalidomide caused birth defects in >10,000 children. While the drug was pulled from the market, thalidomide is currently prescribed to treat conditions including leprosy. As a result, a new generation of babies with thalidomide defects is being born in the developing world. This represents a serious problem, as the mechanisms by which thalidomide disrupts development remain unresolved. This lack of resolution is due, in part, to the absence of an appropriate mammalian model for thalidomide teratogenesis. We test the hypothesis that opossum (Monodelphis domestica) is well suited to model human thalidomide defects. Results suggest that opossum embryos exposed to thalidomide display a range of phenotypes (e.g., heart, craniofacial, limb defects) and penetrance similar to humans. Furthermore, all opossums with thalidomide defects exhibit vascular disruptions. Results therefore support the hypotheses that opossums make a good mammalian model for thalidomide teratogenesis, and that thalidomide can severely disrupt angiogenesis in mammals.

Evolution and Diversity of Transposable Elements in Vertebrate Genomes

Transposable elements (TEs) are selfish genetic elements that mobilize in genomes via transposition or retrotransposition and often make up large fractions of vertebrate genomes. Here, we review the current understanding of vertebrate TE diversity and evolution in the context of recent advances in genome sequencing and assembly techniques. TEs make up 4-60% of assembled vertebrate genomes, and deeply branching lineages such as ray-finned fishes and amphibians generally exhibit a higher TE diversity than the more recent radiations of birds and mammals. Furthermore, the list of taxa with exceptional TE landscapes is growing. We emphasize that the current bottleneck in genome analyses lies in the proper annotation of TEs and provide examples where superficial analyses led to misleading conclusions about genome evolution. Finally, recent advances in long-read sequencing will soon permit access to TE-rich genomic regions that previously resisted assembly including the gigantic, TE-rich genomes of salamanders and lungfishes.

Identification of the mRNA encoding interleukin-6 and its receptor, interleukin-6 receptor α, in five marsupial species

Expressed coding sequences for interleukin-6 (IL-6) and interleukin-6 receptor α (IL-6R) were examined in five marsupial species. Full length expressed coding sequences for IL-6 and IL-6R were identified and characterized in the gray short-tailed opossum (Monodelphis domestica). For IL-6, ∼225 bp fragments of the mRNA sequence were identified in the red-tailed phascogale (Phascogale calura), kultarr (Antechinomys laniger), and stripe-faced dunnart (Sminthopsis macroura), while ∼563 bp fragments of mRNA encoding IL-6R were identified in the red-tailed phascogale, kultarr, stripe-face dunnart and fat-tailed dunnart (Sminthopsis crassicaudata). Relative expression levels of IL-6 and IL-6R were examined in the heart, muscle, lung, liver, spleen and kidney of adult red-tailed phascogales, and IL-6 gene expression was found to be significantly higher in the lung and spleen than the other tissues examined, while the expression of IL-6R was significantly higher in the liver, lung and spleen. These results now serve as a reference point for examining the role and levels of IL-6 and IL-6R in the health and disease of these marsupial species. The pro-tumorigenic nature of IL-6 is of particular interest, and the identification of these IL-6 and IL-6R coding sequences provides a platform for further work to evaluate the potential role of IL-6 in marsupial cancers.

The UT family of MHC class I loci unique to non-eutherian mammals has limited polymorphism and tissue specific patterns of expression in the opossum

BACKGROUND

The Major Histocompatibility Complex (MHC) class I family of genes encode for molecules that have well-conserved structures, but have evolved to perform diverse functions. The availability of the gray, short-tailed opossum, Monodelphis domestica whole genome sequence has allowed for analysis of MHC class I gene content in this marsupial. Utilization of a novel method to search for MHC related domain structures revealed a previously unknown family of MHC class I-related genes. These genes, named UT1-17, are clustered on chromosome 1 in the opossum, unlinked to the MHC region. UT genes are only found in marsupial and monotreme genomes, consistent with being ancient in mammals yet lost in eutherian mammals. This study investigates the expression and polymorphism of the UT loci in the opossum to gain insight into their possible function.

RESULTS

Of the 17 opossum UT genes, most have restricted tissue transcription patterns, with the thymus and skin being the most common sites. Full-length structure of 11 UT transcripts revealed genes varying between five and eight exons, typical for class I family members. There were only two alternative splice variants found. The UT genes also have limited polymorphism and little evidence of positive selection. One locus, UT8, was chosen for further analysis due to its conservation amongst marsupials and generic characteristics. UT8 transcription is limited to developing αβ thymocytes, and is absent from mature αβ T cells in peripheral lymphoid tissues.

CONCLUSION

The overall characteristics and features of UT genes including low polymorphism and restricted tissue expression make it likely that the molecules encoded by UT genes perform roles other than antigenic peptide presentation.

Weird mammals provide insights into the evolution of mammalian sex chromosomes and dosage compensation

The deep divergence of mammalian groups 166 and 190 million years ago (MYA) provide genetic variation to explore the evolution of DNA sequence, gene arrangement and regulation of gene expression in mammals. With encouragement from the founder of the field, Mary Lyon, techniques in cytogenetics and molecular biology were progressively adapted to characterize the sex chromosomes of kangaroos and other marsupials, platypus and echidna-and weird rodent species. Comparative gene mapping reveals the process of sex chromosome evolution from their inception 190 MYA (they are autosomal in platypus) to their inevitable end (the Y has disappeared in two rodent lineages). Our X and Y are relatively young, getting their start with the evolution of the sex-determining SRY gene, which triggered progressive degradation of the Y chromosome. Even more recently, sex chromosomes of placental mammals fused with an autosomal region which now makes up most of the Y. Exploration of gene activity patterns over four decades showed that dosage compensation via X-chromosome inactivation is unique to therian mammals, and that this whole chromosome control process is different in marsupials and absent in monotremes and reptiles, and birds. These differences can be exploited to deduce how mammalian sex chromosomes and epigenetic silencing evolved.

Evolution of the Genotype-to-Phenotype Map and the Cost of Pleiotropy in Mammals

Evolutionary studies have long emphasized that the genetic architecture of traits holds important microevolutionary consequences. Yet, studies comparing the genetic architecture of traits across species are rare, and discussions of the evolution of genetic systems are made on theoretical arguments rather than on empirical evidence. Here, we compared the genetic architecture of cranial traits in two different mammalian model organisms: the gray short-tailed opossum, Monodelphis domestica, and the laboratory mouse, Mus musculus We show that both organisms share a highly polygenic genetic architecture for craniofacial traits, with many loci of small effect. However, these two model species differ significantly in the overall degree of pleiotropy, N, of the genotype-to-phenotype map, with opossums presenting a higher average N They also diverge in their degree of genetic modularity, with opossums presenting less modular patterns of genetic association among traits. We argue that such differences highlight the context dependency of gene effects, with developmental systems shaping the variational properties of genetic systems. Finally, we also demonstrate based on the opossum data that current measurements for the relationship between the mutational effect size and N need to be re-evaluated in relation to the importance of the cost of pleiotropy for mammals.

Coordinately Co-opted Multiple Transposable Elements Constitute an Enhancer for wnt5a Expression in the Mammalian Secondary Palate

Acquisition of cis-regulatory elements is a major driving force of evolution, and there are several examples of developmental enhancers derived from transposable elements (TEs). However, it remains unclear whether one enhancer element could have been produced via cooperation among multiple, yet distinct, TEs during evolution. Here we show that an evolutionarily conserved genomic region named AS3_9 comprises three TEs (AmnSINE1, X6b_DNA and MER117), inserted side-by-side, and functions as a distal enhancer for wnt5a expression during morphogenesis of the mammalian secondary palate. Functional analysis of each TE revealed step-by-step retroposition/transposition and co-option together with acquisition of a binding site for Msx1 for its full enhancer function during mammalian evolution. The present study provides a new perspective suggesting that a huge variety of TEs, in combination, could have accelerated the diversity of cis-regulatory elements involved in morphological evolution.

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

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

Developmental and Regulatory Functions of Na(+) Channel Non-pore-forming β Subunits

Voltage-gated Na(+) channels (VGSCs) isolated from mammalian neurons are heterotrimeric complexes containing one pore-forming α subunit and two non-pore-forming β subunits. In excitable cells, VGSCs are responsible for the initiation of action potentials. VGSC β subunits are type I topology glycoproteins, containing an extracellular amino-terminal immunoglobulin (Ig) domain with homology to many neural cell adhesion molecules (CAMs), a single transmembrane segment, and an intracellular carboxyl-terminal domain. VGSC β subunits are encoded by a gene family that is distinct from the α subunits. While α subunits are expressed in prokaryotes, β subunit orthologs did not arise until after the emergence of vertebrates. β subunits regulate the cell surface expression, subcellular localization, and gating properties of their associated α subunits. In addition, like many other Ig-CAMs, β subunits are involved in cell migration, neurite outgrowth, and axon pathfinding and may function in these roles in the absence of associated α subunits. In sum, these multifunctional proteins are critical for both channel regulation and central nervous system development.

Emergence and Evolution of Hominidae-Specific Coding and Noncoding Genomic Sequences

Family Hominidae, which includes humans and great apes, is recognized for unique complex social behavior and intellectual abilities. Despite the increasing genome data, however, the genomic origin of its phenotypic uniqueness has remained elusive. Clade-specific genes and highly conserved noncoding sequences (HCNSs) are among the high-potential evolutionary candidates involved in driving clade-specific characters and phenotypes. On this premise, we analyzed whole genome sequences along with gene orthology data retrieved from major DNA databases to find Hominidae-specific (HS) genes and HCNSs. We discovered that Down syndrome critical region 4 (DSCR4) is the only experimentally verified gene uniquely present in Hominidae. DSCR4 has no structural homology to any known protein and was inferred to have emerged in several steps through LTR/ERV1, LTR/ERVL retrotransposition, and transversion. Using the genomic distance as neutral evolution threshold, we identified 1,658 HS HCNSs. Polymorphism coverage and derived allele frequency analysis of HS HCNSs showed that these HCNSs are under purifying selection, indicating that they may harbor important functions. They are overrepresented in promoters/untranslated regions, in close proximity of genes involved in sensory perception of sound and developmental process, and also showed a significantly lower nucleosome occupancy probability. Interestingly, many ancestral sequences of the HS HCNSs showed very high evolutionary rates. This suggests that new functions emerged through some kind of positive selection, and then purifying selection started to operate to keep these functions.

Characterisation of CD4 T cells in healthy and diseased koalas (Phascolarctos cinereus) using cell-type-specific monoclonal antibodies

The koala (Phascolarctos cinereus) is an arboreal herbivorous marsupial that is an Australian icon. Koalas in many parts of Australia are under multiple threats including habitat destruction, dog attacks, vehicular accidents, and infectious diseases such as Chlamydia spp. and the koala retrovirus (KoRV), which may contribute to the incidence of lymphoma and leukaemia in this species. Due to a lack of koala-specific immune reagents and assays there is currently no way to adequately analyse the immune response in healthy, diseased or vaccinated animals. This paper reports the production and characterisation of the first anti-koala CD4 monoclonal antibody (mAb). The koala CD4 gene was identified and used to develop recombinant proteins for mAb production. Fluorochrome-conjugated anti-CD4 mAb was used to measure the levels of CD4(+) lymphocytes collected from koala spleens (41.1%, range 20-45.1%) lymph nodes (36.3%, range 19-55.9%) and peripheral blood (23.8%, range 17.3-35%) by flow cytometry. Biotin-conjugated anti-CD4 mAb was used for western blot to determine an approximate size of 52 kDa for the koala CD4 molecule and used in immunohistochemistry to identify CD4(+) cells in the paracortical region and germinal centres of spleen and lymph nodes. Using the anti-CD4 mab we showed that CD4 cells from vaccinated, but not control, koalas proliferated following in vitro stimulation with UV-inactivated Chlamydia pecorum and recombinant chlamydial antigens. Since CD4(+) T cells have been shown to play a pivotal role in clearing chlamydial infection in both human and mouse infections, using this novel antibody will help determine the role CD4(+) T cells play in protection against chlamydial infection in koalas and also enhance our knowledge of how KoRV affects the koala immune system.

The Ensembl gene annotation system

The Ensembl gene annotation system has been used to annotate over 70 different vertebrate species across a wide range of genome projects. Furthermore, it generates the automatic alignment-based annotation for the human and mouse GENCODE gene sets. The system is based on the alignment of biological sequences, including cDNAs, proteins and RNA-seq reads, to the target genome in order to construct candidate transcript models. Careful assessment and filtering of these candidate transcripts ultimately leads to the final gene set, which is made available on the Ensembl website. Here, we describe the annotation process in detail.Database URL: http://www.ensembl.org/index.html.

Angioarchitecture and morphology of temporomandibular joint of Monodelphis domestica

The opossum Monodelphis domestica presents movement of the temporomandibular joint (TMJ) reflecting adaptation to eating habits similar to movement in humans, but the structure of the TMJ is not yet known. Thus, nine young M. domestica, of both sexes were weighed, anesthetized with xylazine (10 mg kg(-1) ), and ketamine (70 mg kg(-1) ) and processed for: 1. The analyses of the macroscopic angioarchitecture after latex injection, as well as the topography of the TMJ; 2. The analysis of microvascularization after injection of Mercox resin and corrosion of soft tissue with NaOH using scanning electron microscopy and; 3. The histological evaluation of the TMJ with an optical microscope. Macroscopic analysis of the latex injected vessels revealed the distribution of the arteries from the common carotid artery, receiving branches of the superficial temporal and maxillary arteries. The mandibular condyle has the long axis in the lateral-lateral direction, and is convex in the anterior-posterior direction. Its topography was determined in relation to the eye and external acoustic meatus. With scanning electron microscopy, microvascularization consists of arterioles of varying diameter (85-15 µm) of the meandering capillary network in the retrodiscal region, and a network of straight capillaries in the TMJ anterior region. Via light microscopy the TMJ has similar histological features to those of humans. These macroscopic, microscopic and ultrastructural data from TMJ of the M. domestica could be a suitable model for TMJ physiology and pathophysiology studies for then speculate on possible human studies. Microsc. Res. Tech. 79:806-813, 2016. © 2016 Wiley Periodicals, Inc.

Massive contribution of transposable elements to mammalian regulatory sequences

Barbara McClintock discovered the existence of transposable elements (TEs) in the late 1940s and initially proposed that they contributed to the gene regulatory program of higher organisms. This controversial idea gained acceptance only much later in the 1990s, when the first examples of TE-derived promoter sequences were uncovered. It is now known that half of the human genome is recognizably derived from TEs. It is thus important to understand the scope and nature of their contribution to gene regulation. Here, we provide a timeline of major discoveries in this area and discuss how transposons have revolutionized our understanding of mammalian genomes, with a special emphasis on the massive contribution of TEs to primate evolution. Our analysis of primate-specific functional elements supports a simple model for the rate at which new functional elements arise in unique and TE-derived DNA. Finally, we discuss some of the challenges and unresolved questions in the field, which need to be addressed in order to fully characterize the impact of TEs on gene regulation, evolution and disease processes.

Enlightening the contribution of the dark matter to the X chromosome inactivation process in mammals

X-chromosome inactivation (XCI) in mammals represents an exceptional example of transcriptional co-regulation occurring at the level of an entire chromosome. XCI is considered as a means to compensate for gene dosage imbalance between sexes, yet the largest part of the chromosome is composed of repeated elements of different nature and origins. Here we consider XCI from a repeat point of view, interrogating the mechanisms for inactivating X chromosome-derived repeated sequences and discussing the contribution of repetitive elements to the silencing process itself and to its evolution.

Accelerated Evolution of Enhancer Hotspots in the Mammal Ancestor

Mammals have evolved remarkably different sensory, reproductive, metabolic, and skeletal systems. To explore the genetic basis for these differences, we developed a comparative genomics approach to scan whole-genome multiple sequence alignments to identify regions that evolved rapidly in an ancestral lineage but are conserved within extant species. This pattern suggests that ancestral changes in function were maintained in descendants. After applying this test to therian mammals, we identified 4,797 accelerated regions, many of which are noncoding and located near developmental transcription factors. We then used mouse transgenic reporter assays to test if noncoding accelerated regions are enhancers and to determine how therian-specific substitutions affect their activity in vivo. We discovered enhancers with expression specific to the therian version in brain regions involved in the hormonal control of milk ejection, uterine contractions, blood pressure, temperature, and visual processing. This work underscores the idea that changes in developmental gene expression are important for mammalian evolution, and it pinpoints candidate genes for unique aspects of mammalian biology.

Epithelial sodium channel (ENaC) family: Phylogeny, structure-function, tissue distribution, and associated inherited diseases

The epithelial sodium channel (ENaC) is composed of three homologous subunits and allows the flow of Na(+) ions across high resistance epithelia, maintaining body salt and water homeostasis. ENaC dependent reabsorption of Na(+) in the kidney tubules regulates extracellular fluid (ECF) volume and blood pressure by modulating osmolarity. In multi-ciliated cells, ENaC is located in cilia and plays an essential role in the regulation of epithelial surface liquid volume necessary for cilial transport of mucus and gametes in the respiratory and reproductive tracts respectively. The subunits that form ENaC (named as alpha, beta, gamma and delta, encoded by genes SCNN1A, SCNN1B, SCNN1G, and SCNN1D) are members of the ENaC/Degenerin superfamily. The earliest appearance of ENaC orthologs is in the genomes of the most ancient vertebrate taxon, Cyclostomata (jawless vertebrates) including lampreys, followed by earliest representatives of Gnathostomata (jawed vertebrates) including cartilaginous sharks. Among Euteleostomi (bony vertebrates), Actinopterygii (ray finned-fishes) branch has lost ENaC genes. Yet, most animals in the Sarcopterygii (lobe-finned fish) branch including Tetrapoda, amphibians and amniotes (lizards, crocodiles, birds, and mammals), have four ENaC paralogs. We compared the sequences of ENaC orthologs from 20 species and established criteria for the identification of ENaC orthologs and paralogs, and their distinction from other members of the ENaC/Degenerin superfamily, especially ASIC family. Differences between ENaCs and ASICs are summarized in view of their physiological functions and tissue distributions. Structural motifs that are conserved throughout vertebrate ENaCs are highlighted. We also present a comparative overview of the genotype-phenotype relationships in inherited diseases associated with ENaC mutations, including multisystem pseudohypoaldosteronism (PHA1B), Liddle syndrome, cystic fibrosis-like disease and essential hypertension.

Genetic Diversity on the Human X Chromosome Does Not Support a Strict Pseudoautosomal Boundary

Unlike the autosomes, recombination between the X chromosome and the Y chromosome is often thought to be constrained to two small pseudoautosomal regions (PARs) at the tips of each sex chromosome. PAR1 spans the first 2.7 Mb of the proximal arm of the human sex chromosomes, whereas the much smaller PAR2 encompasses the distal 320 kb of the long arm of each sex chromosome. In addition to PAR1 and PAR2, there is a human-specific X-transposed region that was duplicated from the X to the Y chromosome. The X-transposed region is often not excluded from X-specific analyses, unlike the PARs, because it is not thought to routinely recombine. Genetic diversity is expected to be higher in recombining regions than in nonrecombining regions because recombination reduces the effect of linked selection. In this study, we investigated patterns of genetic diversity in noncoding regions across the entire X chromosome of a global sample of 26 unrelated genetic females. We found that genetic diversity in PAR1 is significantly greater than in the nonrecombining regions (nonPARs). However, rather than an abrupt drop in diversity at the pseudoautosomal boundary, there is a gradual reduction in diversity from the recombining through the nonrecombining regions, suggesting that recombination between the human sex chromosomes spans across the currently defined pseudoautosomal boundary. A consequence of recombination spanning this boundary potentially includes increasing the rate of sex-linked disorders (e.g., de la Chapelle) and sex chromosome aneuploidies. In contrast, diversity in PAR2 is not significantly elevated compared to the nonPARs, suggesting that recombination is not obligatory in PAR2. Finally, diversity in the X-transposed region is higher than in the surrounding nonPARs, providing evidence that recombination may occur with some frequency between the X and Y chromosomes in the X-transposed region.

Characterization of DNA methylation as a function of biological complexity via dinucleotide inter-distances

We perform a statistical study of the distances between successive occurrences of a given dinucleotide in the DNA sequence for a number of organisms of different complexity. Our analysis highlights peculiar features of the CG dinucleotide distribution in mammalian DNA, pointing towards a connection with the role of such dinucleotide in DNA methylation. While the CG distributions of mammals exhibit exponential tails with comparable parameters, the picture for the other organisms studied (e.g. fish, insects, bacteria and viruses) is more heterogeneous, possibly because in these organisms DNA methylation has different functional roles. Our analysis suggests that the distribution of the distances between CG dinucleotides provides useful insights into characterizing and classifying organisms in terms of methylation functionalities.

Transposons, Genome Size, and Evolutionary Insights in Animals

The relationship between genome size and the percentage of transposons in 161 animal species evidenced that variations in genome size are linked to the amplification or the contraction of transposable elements. The activity of transposable elements could represent a response to environmental stressors. Indeed, although with different trends in protostomes and deuterostomes, comprehensive changes in genome size were recorded in concomitance with particular periods of evolutionary history or adaptations to specific environments. During evolution, genome size and the presence of transposable elements have influenced structural and functional parameters of genomes and cells. Changes of these parameters have had an impact on morphological and functional characteristics of the organism on which natural selection directly acts. Therefore, the current situation represents a balance between insertion and amplification of transposons and the mechanisms responsible for their deletion or for decreasing their activity. Among the latter, methylation and the silencing action of small RNAs likely represent the most frequent mechanisms.

Mucosal-associated invariant T cells in autoimmunity, immune-mediated diseases and airways disease

Mucosal-associated invariant T (MAIT) cells are a novel class of innate-like T cells, expressing a semi-invariant T-cell receptor (TCR) and able to recognize small molecules presented on the non-polymorphic MHC-related protein 1. Their intrinsic effector-memory phenotype, enabling secretion of pro-inflammatory cytokines, and their relative abundance in humans imply a significant potential to contribute to autoimmune processes. However, as MAIT cells were unknown until recently and specific immunological tools were unavailable, little is known of their roles in disease. Here I review observations from clinical studies and animal models of autoimmune and immune-mediated diseases including the roles of MAIT cells in systemic lupus erythematosus, rheumatoid arthritis, multiple sclerosis, inflammatory bowel disease and airways diseases. MAIT cell deficiencies are frequently observed in peripheral blood, and at sites of disease such as the airways in asthma. However, MAIT cells have a specific sensitivity to suppression by therapeutic corticosteroids that may confound many of these observations, as may the tendency of the surface marker CD161 to activation-induced down-regulation. Nonetheless, the dependence on bacteria for the development of MAIT cells suggests a potentially important protective role linking the influences of early life microbial exposures and subsequent development of autoimmunity. Conversely, MAIT cells could contribute to chronic inflammation either through TCR-independent activation, or potentially by TCR recognition of as yet undiscovered ligands. Future research will be greatly facilitated by the immunological tools that are now available, including murine genetic models and human and murine specific tetramers.

The Interrelationships of Placental Mammals and the Limits of Phylogenetic Inference

Placental mammals comprise three principal clades: Afrotheria (e.g., elephants and tenrecs), Xenarthra (e.g., armadillos and sloths), and Boreoeutheria (all other placental mammals), the relationships among which are the subject of controversy and a touchstone for debate on the limits of phylogenetic inference. Previous analyses have found support for all three hypotheses, leading some to conclude that this phylogenetic problem might be impossible to resolve due to the compounded effects of incomplete lineage sorting (ILS) and a rapid radiation. Here we show, using a genome scale nucleotide data set, microRNAs, and the reanalysis of the three largest previously published amino acid data sets, that the root of Placentalia lies between Atlantogenata and Boreoeutheria. Although we found evidence for ILS in early placental evolution, we are able to reject previous conclusions that the placental root is a hard polytomy that cannot be resolved. Reanalyses of previous data sets recover Atlantogenata + Boreoeutheria and show that contradictory results are a consequence of poorly fitting evolutionary models; instead, when the evolutionary process is better-modeled, all data sets converge on Atlantogenata. Our Bayesian molecular clock analysis estimates that marsupials diverged from placentals 157-170 Ma, crown Placentalia diverged 86-100 Ma, and crown Atlantogenata diverged 84-97 Ma. Our results are compatible with placental diversification being driven by dispersal rather than vicariance mechanisms, postdating early phases in the protracted opening of the Atlantic Ocean.

Preliminary genomic survey and sequence analysis of the complement system in non-eutherian mammals

The complement system is a major mediator of the vertebrate immune system, which functions in both innate and specific immune responses. It comprises more than 30 proteins working to remove foreign cells by way of anaphylatoxins, opsonins or the membrane attack complex. Over the last few years, whole genome sequences of non-eutherian mammals (marsupials and a monotreme), the gray short-tailed opossum (Monodelphis domestica), tammar wallaby (Macropus eugenii), Tasmanian devil (Sarcophilus harrisii), koala (Phascolarctos cinereus) and platypus (Ornithorhynchus anatinus), have become publicly available. Using these sequences, we have identified an array of complement components in non-eutherians using online search tools and algorithms. Of 57 complement and complement-related genes investigated, we identified 46 in the gray short-tailed opossum genome, 27 in the tammar wallaby genome, 44 in the Tasmanian devil genome, 47 in the koala genome and 40 in the platypus genome. The results of this study confirm the presence of key complement components in the immune repertoire of non-eutherian mammals and provide a platform for future studies on immune protection in young marsupials.

The devil is in the details: Transposable element analysis of the Tasmanian devil genome

The third marsupial genome was sequenced from the Tasmanian devil (Sarcophilus harrisii), a species that currently is driven to extinction by a rare transmissible cancer. The transposable element (TE) landscape of the Tasmanian devil genome revealed that the main driver of retrotransposition the Long INterspersed Element 1 (LINE1) seem to have become inactivated during the past 12 million years. Strangely, the Short INterspersed Elements (SINE), that normally hijacks the LINE1 retrotransposition system, became inactive prior to LINE1 at around 30 million years ago. The SINE inactivation was in vitro verified in several species. Here I discuss that the apparent LINE1 inactivation might be caused by a genome assembly artifact. The repetitive fraction of any genome is highly complex to assemble and the observed problems are not unique to the Tasmanian devil genome.

Evolution of vertebrate sex chromosomes and dosage compensation

Differentiated sex chromosomes in mammals and other vertebrates evolved independently but in strikingly similar ways. Vertebrates with differentiated sex chromosomes share the problems of the unequal expression of the genes borne on sex chromosomes, both between the sexes and with respect to autosomes. Dosage compensation of genes on sex chromosomes is surprisingly variable - and can even be absent - in different vertebrate groups. Systems that compensate for different gene dosages include a wide range of global, regional and gene-by-gene processes that differ in their extent and their molecular mechanisms. However, many elements of these control systems are similar across distant phylogenetic divisions and show parallels to other gene silencing systems. These dosage systems cannot be identical by descent but were probably constructed from elements of ancient silencing mechanisms that are ubiquitous among vertebrates and shared throughout eukaryotes.

Identification and analysis of divergent immune gene families within the Tasmanian devil genome

Background

The Tasmanian devil (Sarcophilus harrisii) is being threatened with extinction in the wild by a disease known as devil facial tumour disease (DFTD). In order to prevent the spread of this disease a thorough understanding of the Tasmanian devil immune system and its response to the disease is required. In 2011 and 2012 two genome sequencing projects of the Tasmania devil were released. This has provided us with the raw data required to begin to investigate the Tasmanian devil immunome in depth. In this study we characterise immune gene families of the Tasmanian devil. We focus on immunoglobulins, T cell receptors and cytokine families.

Results

We identify and describe 119 cytokines including 40 interleukins, 39 chemokines, 8 interferons, 18 tumour necrosis family cytokines and 14 additional cytokines. Constant regions for immunoglobulins and T cell receptors were also identified. The repertoire of genes in these families was similar to the opossum, however devil specific duplications were seen and orthologs to eutherian genes not previously identified in any marsupial were also identified.

Conclusions

By using multiple data sources as well as targeted search methods, highly divergent genes across the Tasmanian devil immune system were identified and characterised. This understanding will allow for the development of devil specific assays and reagents and allow for future studies into the immune response of the Tasmanian devil immune system to DFTD.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-2206-9) contains supplementary material, which is available to authorized users.