Conservation of Glutathione Transferase mRNA and Protein Sequences Similar to Human and Horse Alpha Class GST A3-3 across Dog, Goat, and Opossum Species.

The glutathione transferase A3-3 (GST A3-3) homodimeric enzyme is the most efficient enzyme that catalyzes isomerization of the precursors of testosterone, estradiol, and progesterone in the gonads of humans and horses. However, the presence of GST A3-3 orthologs with equally high ketosteroid isomerase activity has not been verified in other mammalian species, even though pig and cattle homologs have been cloned and studied. Identifying GSTA3 genes is a challenge because of multiple GSTA gene duplications (e.g., 12 in the human genome); consequently, the GSTA3 gene is not annotated in most genomes. To improve our understanding of GSTA3 gene products and their functions across diverse mammalian species, we cloned homologs of the horse and human GSTA3 mRNAs from the testes of a dog, goat, and gray short-tailed opossum, the genomes of which all currently lack GSTA3 gene annotations. The resultant novel GSTA3 mRNA and inferred protein sequences had a high level of conservation with human GSTA3 mRNA and protein sequences (≥70% and ≥64% identities, respectively). Sequence conservation was also apparent for the 12 residues of the "H-site" in the 222 amino acid GSTA3 protein that is known to interact with the steroid substrates. Modeling predicted that the dog GSTA3-3 may be a more active ketosteroid isomerase than the corresponding goat or opossum enzymes. However, expression of the GSTA3 gene was higher in liver than in other dog tissue. Our results improve understanding of the active sites of mammalian GST A3-3 enzymes, inhibitors of which might be useful for reducing steroidogenesis for medical purposes, such as fertility control or treatment of steroid-dependent diseases.

Origin and Evolution of Marsupial-specific Imprinting Clusters Through Lineage-specific Gene Duplications and Acquisition of Promoter Differential Methylation.

Genomic imprinting is a parent-of-origin-specific expression phenomenon that plays fundamental roles in many biological processes. In animals, imprinting is only observed in therian mammals, with ∼200 imprinted genes known in humans and mice. The imprinting pattern in marsupials has been minimally investigated by examining orthologs to known eutherian imprinted genes. To identify marsupial-specific imprinting in an unbiased way, we performed RNA-seq studies on samples of fetal brain and placenta from the reciprocal cross progeny of two laboratory opossum stocks. We inferred allele-specific expression for >3,000 expressed genes and discovered/validated 13 imprinted genes, including three previously known imprinted genes, Igf2r, Peg10, and H19. We estimate that marsupials imprint ∼60 autosomal genes, which is a much smaller set compared with eutherians. Among the nine novel imprinted genes, three noncoding RNAs have no known homologs in eutherian mammals, while the remaining genes have important functions in pluripotency, transcription regulation, nucleolar homeostasis, and neural differentiation. Methylation analyses at promoter CpG islands revealed differentially methylated regions in five of these marsupial-specific imprinted genes, suggesting that differential methylation is a common mechanism in the epigenetic regulation of marsupial imprinting. Clustering and co-regulation were observed at marsupial imprinting loci Pou5f3-Npdc1 and Nkrfl-Ipncr2, but eutherian-type multi-gene imprinting clusters were not detected. Also differing from eutherian mammals, the brain and placenta imprinting profiles are remarkably similar in opossums, presumably due to the shared origin of these organs from the trophectoderm. Our results contribute to a fuller understanding of the origin, evolution, and mechanisms of genomic imprinting in therian mammals.

Strategies for meiotic sex chromosome dynamics and telomeric elongation in Marsupials.

During meiotic prophase I, homologous chromosomes pair, synapse and recombine in a tightly regulated process that ensures the generation of genetically variable haploid gametes. Although the mechanisms underlying meiotic cell division have been well studied in model species, our understanding of the dynamics of meiotic prophase I in non-traditional model mammals remains in its infancy. Here, we reveal key meiotic features in previously uncharacterised marsupial species (the tammar wallaby and the fat-tailed dunnart), plus the fat-tailed mouse opossum, with a focus on sex chromosome pairing strategies, recombination and meiotic telomere homeostasis. We uncovered differences between phylogroups with important functional and evolutionary implications. First, sex chromosomes, which lack a pseudo-autosomal region in marsupials, had species specific pairing and silencing strategies, with implications for sex chromosome evolution. Second, we detected two waves of γH2AX accumulation during prophase I. The first wave was accompanied by low γH2AX levels on autosomes, which correlated with the low recombination rates that distinguish marsupials from eutherian mammals. In the second wave, γH2AX was restricted to sex chromosomes in all three species, which correlated with transcription from the X in tammar wallaby. This suggests non-canonical functions of γH2AX on meiotic sex chromosomes. Finally, we uncover evidence for telomere elongation in primary spermatocytes of the fat-tailed dunnart, a unique strategy within mammals. Our results provide new insights into meiotic progression and telomere homeostasis in marsupials, highlighting the importance of capturing the diversity of meiotic strategies within mammals.

Evolutionary Characterization of the Short Protein SPAAR.

Microproteins (<100 amino acids) are receiving increasing recognition as important participants in numerous biological processes, but their evolutionary dynamics are poorly understood. SPAAR is a recently discovered microprotein that regulates muscle regeneration and angiogenesis through interactions with conserved signaling pathways. Interestingly, SPAAR does not belong to any known protein family and has known homologs exclusively among placental mammals. This lack of distant homology could be caused by challenges in homology detection of short sequences, or it could indicate a recent de novo emergence from a noncoding sequence. By integrating syntenic alignments and homology searches, we identify SPAAR orthologs in marsupials and monotremes, establishing that SPAAR has existed at least since the emergence of mammals. SPAAR shows substantial primary sequence divergence but retains a conserved protein structure. In primates, we infer two independent evolutionary events leading to the de novo origination of 5' elongated isoforms of SPAAR from a noncoding sequence and find evidence of adaptive evolution in this extended region. Thus, SPAAR may be of ancient origin, but it appears to be experiencing continual evolutionary innovation in mammals.

Developmental and evolutionary dynamics of cis-regulatory elements in mouse cerebellar cells.

Organ development is orchestrated by cell- and time-specific gene regulatory networks. In this study, we investigated the regulatory basis of mouse cerebellum development from early neurogenesis to adulthood. By acquiring snATAC-seq (single-nucleus assay for transposase accessible chromatin using sequencing) profiles for ~90,000 cells spanning 11 stages, we mapped cerebellar cell types and identified candidate cis-regulatory elements (CREs). We detected extensive spatiotemporal heterogeneity among progenitor cells and a gradual divergence in the regulatory programs of cerebellar neurons during differentiation. Comparisons to vertebrate genomes and snATAC-seq profiles for ∼20,000 cerebellar cells from the marsupial opossum revealed a shared decrease in CRE conservation during development and differentiation as well as differences in constraint between cell types. Our work delineates the developmental and evolutionary dynamics of gene regulation in cerebellar cells and provides insights into mammalian organ development.

A Chromosome-Level Genome of the Agile Gracile Mouse Opossum (Gracilinanus agilis).

There are more than 100 species of American didelphid marsupials (opossums and mouse opossums). Limited genomic resources for didelphids exists, with only two publicly available genome assemblies compared with dozens in the case of their Australasian counterparts. This discrepancy impedes evolutionary and ecological research. To address this gap, we assembled a high-quality chromosome-level genome of the agile gracile mouse opossum (Gracilinanus agilis) using a combination of stLFR sequencing, polishing with mate-pair data, and anchoring onto pseudochromosomes using Hi-C. This species employs a rare life-history strategy, semelparity, and all G. agilis males and most females die at the end of their first breeding season after succumbing to stress and exhaustion. The 3.7-Gb chromosome-level assembly, with 92.6% anchored onto pseudochromosomes, has a scaffold N50 of 683.5 Mb and a contig N50 of 56.9 kb. The genome assembly shows high completeness, with a mammalian BUSCO score of 88.1%. Around 49.7% of the genome contains repetitive elements. Gene annotation yielded 24,425 genes, of which 83.9% were functionally annotated. The G. agilis genome is an important resource for future studies of marsupial biology, evolution, and conservation.

Molecular Phylogeny of Cryptonanus (Didelphidae: Thylamyini): Evidence for a recent and complex diversification in South American open biomes.

Systematic revisions of South American marsupials have contributed to our knowledge about genus and species diversity in the last decades, including studies of the most recently described genus Cryptonanus (Didelphidae), currently comprising four recognized species. Herein we provide the first phylogeny for these mouse opossums based on comprehensive sampling, including representatives from all nominal taxa, encompassing most of the geographic distribution of the genus while also extending its known range. The taxonomic status of Cryptonanus species was explored by analyses of multiple mitochondrial and nuclear DNA markers to assess phylogenetic relationships and to provide divergence time estimates, species delimitations and biogeographical hypotheses. Cryptonanus monophyly remained highly supported despite the inclusion of abundant new data from more than a hundred specimens, comprising 10 independent evolutionary lineages. Species-complexes within valid nominal taxa reveal higher species richness in the genus. Based on divergence estimates from a dated phylogeny, we suggest that Cryptonanus diversified along the Quaternary, with speciation events occurring well into the Pleistocene. The best supported biogeographical hypothesis endorses speciation by vicariance and subset speciation across open formations in shaping the evolutionary history of this didelphid genus, strongly associated with dry tropical landscapes of South America.

Molecular Evolution of CatSper in Mammals and Function of Sperm Hyperactivation in Gray Short-Tailed Opossum.

Males have evolved species-specifical sperm morphology and swimming patterns to adapt to different fertilization environments. In eutherians, only a small fraction of the sperm overcome the diverse obstacles in the female reproductive tract and successfully migrate to the fertilizing site. Sperm arriving at the fertilizing site show hyperactivated motility, a unique motility pattern displaying asymmetric beating of sperm flagella with increased amplitude. This motility change is triggered by Ca2+ influx through the sperm-specific ion channel, CatSper. However, the current understanding of the CatSper function and its molecular regulation is limited in eutherians. Here, we report molecular evolution and conservation of the CatSper channel in the genome throughout eutherians and marsupials. Sequence analyses reveal that CatSper proteins are slowly evolved in marsupials. Using an American marsupial, gray short-tailed opossum (Monodelphis domestica), we demonstrate the expression of CatSper in testes and its function in hyperactivation and unpairing of sperm. We demonstrate that a conserved IQ-like motif in CatSperζ is required for CatSperζ interaction with the pH-tuned Ca2+ sensor, EFCAB9, for regulating CatSper activity. Recombinant opossum EFCAB9 can interact with mouse CatSperζ despite high sequence divergence of CatSperζ among CatSper subunits in therians. Our finding suggests that molecular characteristics and functions of CatSper are evolutionarily conserved in gray short-tailed opossum, unraveling the significance of sperm hyperactivation and fertilization in marsupials for the first time.

Anatomical comparison across heads, fore- and hindlimbs in mammals using network models.

Animal body parts evolve with variable degrees of integration that nonetheless yield functional adult phenotypes: but, how? The analysis of modularity with Anatomical Network Analysis (AnNA) is used to quantitatively determine phenotypic modules based on the physical connection among anatomical elements, an approach that is valuable to understand developmental and evolutionary constraints. We created anatomical network models of the head, forelimb, and hindlimb of two taxa considered to represent a 'generalized' eutherian (placental: mouse) and metatherian (marsupial: opossum) anatomical configuration and compared them with our species, which has a derived eutherian configuration. In these models, nodes represent anatomical units and links represent their physical connection. Here, we aimed to identify: (1) the commonalities and differences in modularity between species, (2) whether modules present a potential phylogenetic character, and (3) whether modules preferentially reflect either developmental or functional aspects of anatomy, or a mix of both. We predicted differences between networks of metatherian and eutherian mammals that would best be explained by functional constraints, versus by constraints of development and/or phylogeny. The topology of contacts between bones, muscles, and bones + muscles showed that, among all three species, skeletal networks were more similar than musculoskeletal networks. There was no clear indication that humans and mice are more alike when compared to the opossum overall, even though their musculoskeletal and skeletal networks of fore- and hindlimbs are slightly more similar. Differences were greatest among musculoskeletal networks of heads and next of forelimbs, which showed more variation than hindlimbs, supporting previous anatomical studies indicating that in general the configuration of the hindlimbs changes less across evolutionary history. Most observations regarding the anatomical networks seem to be best explained by function, but an exception is the adult opossum ear ossicles. These ear bones might form an independent module because the incus and malleus are involved in forming a functional primary jaw that enables the neonate to attach to the teat, where this newborn will complete its development. Additionally, the human data show a specialized digit 1 module (thumb/big toe) in both limb types, likely the result of functional and evolutionary pressures, as our ape ancestors had highly movable big toes and thumbs.

Meiotic Cells Counteract Programmed Retrotransposon Activation via RNA-Binding Translational Repressor Assemblies.

Retrotransposon proliferation poses a threat to germline integrity. While retrotransposons must be activated in developing germ cells in order to survive and propagate, how they are selectively activated in the context of meiosis is unclear. We demonstrate that the transcriptional activation of Ty3/Gypsy retrotransposons and host defense are controlled by master meiotic regulators. We show that budding yeast Ty3/Gypsy co-opts binding sites of the essential meiotic transcription factor Ndt80 upstream of the integration site, thereby tightly linking its transcriptional activation to meiotic progression. We also elucidate how yeast cells thwart Ty3/Gypsy proliferation by blocking translation of the retrotransposon mRNA using amyloid-like assemblies of the RNA-binding protein Rim4. In mammals, several inactive Ty3/Gypsy elements are undergoing domestication. We show that mammals utilize equivalent master meiotic regulators (Stra8, Mybl1, Dazl) to regulate Ty3/Gypsy-derived genes in developing gametes. Our findings inform how genes that are evolving from retrotransposons can build upon existing regulatory networks during domestication.

Transcriptome and translatome co-evolution in mammals.

Gene-expression programs define shared and species-specific phenotypes, but their evolution remains largely uncharacterized beyond the transcriptome layer1. Here we report an analysis of the co-evolution of translatomes and transcriptomes using ribosome-profiling and matched RNA-sequencing data for three organs (brain, liver and testis) in five mammals (human, macaque, mouse, opossum and platypus) and a bird (chicken). Our within-species analyses reveal that translational regulation is widespread in the different organs, in particular across the spermatogenic cell types of the testis. The between-species divergence in gene expression is around 20% lower at the translatome layer than at the transcriptome layer owing to extensive buffering between the expression layers, which especially preserved old, essential and housekeeping genes. Translational upregulation specifically counterbalanced global dosage reductions during the evolution of sex chromosomes and the effects of meiotic sex-chromosome inactivation during spermatogenesis. Despite the overall prevalence of buffering, some genes evolved faster at the translatome layer-potentially indicating adaptive changes in expression; testis tissue shows the highest fraction of such genes. Further analyses incorporating mass spectrometry proteomics data establish that the co-evolution of transcriptomes and translatomes is reflected at the proteome layer. Together, our work uncovers co-evolutionary patterns and associated selective forces across the expression layers, and provides a resource for understanding their interplay in mammalian organs.

Opossum Cathelicidins Exhibit Antimicrobial Activity Against a Broad Spectrum of Pathogens Including West Nile Virus.

This study aimed to characterize cathelicidins from the gray short-tailed opossum in silico and experimentally validate their antimicrobial effects against various pathogenic bacteria and West Nile virus (WNV). Genome-wide in silico analysis against the current genome assembly of the gray short-tailed opossum yielded 56 classical antimicrobial peptides (AMPs) from eight different families, among which 19 cathelicidins, namely ModoCath1 - 19, were analyzed in silico to predict their antimicrobial domains and three of which, ModoCath1, -5, and -6, were further experimentally evaluated for their antimicrobial activity, and were found to exhibit a wide spectrum of antimicroial effects against a panel of gram-positive and gram-negative bacterial strains. In addition, these peptides displayed low-to-moderate cytotoxicity in mammalian cells as well as stability in serum and various salt and pH conditions. Circular dichroism analysis of the spectra resulting from interactions between ModoCaths and lipopolysaccharides (LPS) showed formation of a helical structure, while a dual-dye membrane disruption assay and scanning electron microscopy analysis revealed that ModoCaths exerted bactericidal effects by causing membrane damage. Furthermore, ModoCath5 displayed potent antiviral activity against WNV by inhibiting viral replication, suggesting that opossum cathelicidins may serve as potentially novel antimicrobial endogenous substances of mammalian origin, considering their large number. Moreover, analysis of publicly available RNA-seq data revealed the expression of eight ModoCaths from five different tissues, suggesting that gray short-tailed opossums may be an interesting source of cathelicidins with diverse characteristics.

Didelphis albiventris: an overview of unprecedented transcriptome sequencing of the white-eared opossum.

BACKGROUND: The white-eared opossum (Didelphis albiventris) is widely distributed throughout Brazil and South America. It has been used as an animal model for studying different scientific questions ranging from the restoration of degraded green areas to medical aspects of Chagas disease, leishmaniasis and resistance against snake venom. As a marsupial, D. albiventris can also contribute to the understanding of the molecular mechanisms that govern the different stages of organogenesis. Opossum joeys are born after only 13 days, and the final stages of organogenesis occur when the neonates are inside the pouch, depending on lactation. As neither the genome of this opossum species nor its transcriptome has been completely sequenced, the use of D. albiventris as an animal model is limited. In this work, we sequenced the D. albiventris transcriptome by RNA-seq to obtain the first catalogue of differentially expressed (DE) genes and gene ontology (GO) annotations during the neonatal stages of marsupial development. RESULTS: The D. albiventris transcriptome was obtained from whole neonates harvested at birth (P0), at 5 days of age (P5) and at 10 days of age (P10). The de novo assembly of these transcripts generated 85,338 transcripts. Approximately 30% of these transcripts could be mapped against the amino acid sequences of M. domestica, the evolutionarily closest relative of D. albiventris to be sequenced thus far. Among the expressed transcripts, 2077 were found to be DE between P0 and P5, 13,780 between P0 and P10, and 1453 between P5 and P10. The enriched GO terms were mainly related to the immune system, blood tissue development and differentiation, vision, hearing, digestion, the CNS and limb development. CONCLUSIONS: The elucidation of opossum transcriptomes provides an out-group for better understanding the distinct characteristics associated with the evolution of mammalian species. This study provides the first transcriptome sequences and catalogue of genes for a marsupial species at different neonatal stages, allowing the study of the mechanisms involved in organogenesis.

Developmental dynamics of lncRNAs across mammalian organs and species.

Although many long noncoding RNAs (lncRNAs) have been identified in human and other mammalian genomes, there has been limited systematic functional characterization of these elements. In particular, the contribution of lncRNAs to organ development remains largely unexplored. Here we analyse the expression patterns of lncRNAs across developmental time points in seven major organs, from early organogenesis to adulthood, in seven species (human, rhesus macaque, mouse, rat, rabbit, opossum and chicken). Our analyses identified approximately 15,000 to 35,000 candidate lncRNAs in each species, most of which show species specificity. We characterized the expression patterns of lncRNAs across developmental stages, and found many with dynamic expression patterns across time that show signatures of enrichment for functionality. During development, there is a transition from broadly expressed and conserved lncRNAs towards an increasing number of lineage- and organ-specific lncRNAs. Our study provides a resource of candidate lncRNAs and their patterns of expression and evolutionary conservation across mammalian organ development.

Gene expression across mammalian organ development.

The evolution of gene expression in mammalian organ development remains largely uncharacterized. Here we report the transcriptomes of seven organs (cerebrum, cerebellum, heart, kidney, liver, ovary and testis) across developmental time points from early organogenesis to adulthood for human, rhesus macaque, mouse, rat, rabbit, opossum and chicken. Comparisons of gene expression patterns identified correspondences of developmental stages across species, and differences in the timing of key events during the development of the gonads. We found that the breadth of gene expression and the extent of purifying selection gradually decrease during development, whereas the amount of positive selection and expression of new genes increase. We identified differences in the temporal trajectories of expression of individual genes across species, with brain tissues showing the smallest percentage of trajectory changes, and the liver and testis showing the largest. Our work provides a resource of developmental transcriptomes of seven organs across seven species, and comparative analyses that characterize the development and evolution of mammalian organs.

The SCJ Small Parsimony Problem for Weighted Gene Adjacencies.

Reconstructing ancestral gene orders in a given phylogeny is a classical problem in comparative genomics. Most existing methods compare conserved features in extant genomes in the phylogeny to define potential ancestral gene adjacencies, and either try to reconstruct all ancestral genomes under a global evolutionary parsimony criterion, or, focusing on a single ancestral genome, use a scaffolding approach to select a subset of ancestral gene adjacencies, generally aiming at reducing the fragmentation of the reconstructed ancestral genome. In this paper, we describe an exact algorithm for the Small Parsimony Problem that combines both approaches. We consider that gene adjacencies at internal nodes of the species phylogeny are weighted, and we introduce an objective function defined as a convex combination of these weights and the evolutionary cost under the Single-Cut-or-Join (SCJ) model. The weights of ancestral gene adjacencies can, e.g., be obtained through the recent availability of ancient DNA sequencing data, which provide a direct hint at the genome structure of the considered ancestor, or through probabilistic analysis of gene adjacencies evolution. We show the NP-hardness of our problem variant and propose a Fixed-Parameter Tractable algorithm based on the Sankoff-Rousseau dynamic programming algorithm that also allows to sample co-optimal solutions. We apply our approach to mammalian and bacterial data providing different degrees of complexity. We show that including adjacency weights in the objective has a significant impact in reducing the fragmentation of the reconstructed ancestral gene orders. An implementation is available at http://github.com/nluhmann/PhySca.

Functional Conservation of a Developmental Switch in Mammals since the Jurassic Age.

ThPOK is a "master regulator" of T lymphocyte lineage choice, whose presence or absence is sufficient to dictate development to the CD4 or CD8 lineages, respectively. Induction of ThPOK is transcriptionally regulated, via a lineage-specific silencer element, SilThPOK. Here, we take advantage of the available genome sequence data as well as site-specific gene targeting technology, to evaluate the functional conservation of ThPOK regulation across mammalian evolution, and assess the importance of motif grammar (order and orientation of TF binding sites) on SilThPOK function in vivo. We make three important points: First, the SilThPOK is present in marsupial and placental mammals, but is not found in available genome assemblies of nonmammalian vertebrates, indicating that it arose after divergence of mammals from other vertebrates. Secondly, by replacing the murine SilThPOK in situ with its marsupial equivalent using a knockin approach, we demonstrate that the marsupial SilThPOK supports correct CD4 T lymphocyte lineage-specification in mice. To our knowledge, this is the first in vivo demonstration of functional equivalency for a silencer element between marsupial and placental mammals using a definitive knockin approach. Finally, we show that alteration of the position/orientation of a highly conserved region within the murine SilThPOK is sufficient to destroy silencer activity in vivo, demonstrating that motif grammar of this "solid" synteny block is critical for silencer function. Dependence of SilThPOK function on motif grammar conserved since the mid-Jurassic age, 165 Ma, suggests that the SilThPOK operates as a silenceosome, by analogy with the previously proposed enhanceosome model.

Transcriptomic analysis of skin pigmentation variation in the Virginia opossum (Didelphis virginiana).

Skin pigmentation and coat pigmentation are two of the best-studied examples of traits under natural selection given their quantifiable fitness interactions with the environment (e.g., camouflage) and signalling with other organisms (e.g., warning coloration). Previous morphological studies have found that skin pigmentation variation in the Virginia opossum (Didelphis virginiana) is associated with variation in precipitation and temperatures across its distribution range following Gloger's rule (lighter pigmentation in temperate environments). To investigate the molecular mechanism associated with skin pigmentation variation, we used RNA-Seq and quantified gene expression of wild opossums from tropical and temperate populations. Using differential expression analysis and a co-expression network approach, we found that expression variation in genes with melanocytic and immune functions is significantly associated with the degree of skin pigmentation variation and may be underlying this phenotypic difference. We also found evidence suggesting that the Wnt/ß-catenin signalling pathway might be regulating the depigmentation observed in temperate populations. Based on our study results, we present several alternative hypotheses that may explain Gloger's rule pattern of skin pigmentation variation in opossum, including changes in pathogen diversity supporting a pathogen-resistant hypothesis, thermal stress associated with temperate environments, and pleiotropic and epistatic interactions between melanocytic and immune genes.

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.

Superior ab initio identification, annotation and characterisation of TEs and segmental duplications from genome assemblies.

Transposable Elements (TEs) are mobile DNA sequences that make up significant fractions of amniote genomes. However, they are difficult to detect and annotate ab initio because of their variable features, lengths and clade-specific variants. We have addressed this problem by refining and developing a Comprehensive ab initio Repeat Pipeline (CARP) to identify and cluster TEs and other repetitive sequences in genome assemblies. The pipeline begins with a pairwise alignment using krishna, a custom aligner. Single linkage clustering is then carried out to produce families of repetitive elements. Consensus sequences are then filtered for protein coding genes and then annotated using Repbase and a custom library of retrovirus and reverse transcriptase sequences. This process yields three types of family: fully annotated, partially annotated and unannotated. Fully annotated families reflect recently diverged/young known TEs present in Repbase. The remaining two types of families contain a mixture of novel TEs and segmental duplications. These can be resolved by aligning these consensus sequences back to the genome to assess copy number vs. length distribution. Our pipeline has three significant advantages compared to other methods for ab initio repeat identification: 1) we generate not only consensus sequences, but keep the genomic intervals for the original aligned sequences, allowing straightforward analysis of evolutionary dynamics, 2) consensus sequences represent low-divergence, recently/currently active TE families, 3) segmental duplications are annotated as a useful by-product. We have compared our ab initio repeat annotations for 7 genome assemblies to other methods and demonstrate that CARP compares favourably with RepeatModeler, the most widely used repeat annotation package.