Adaptive Evolution and Functional Differentiation of Testis-Expressed Genes in Theria

Gene expression patterns differ in different tissues, and the expression pattern of genes in the mammalian testis is known to be extremely variable in different species. To clarify how the testis transcriptomic pattern has evolved in particular species, we examined the evolution of the adult testis transcriptome in Theria using 10 species: two marsupials (opossum and Tasmanian devil), six eutherian (placental) mammals (human, chimpanzee, bonobo, gorilla, rhesus macaque, and mouse), and two outgroup species (platypus and chicken). We show that 22 testis-expressed genes are marsupial-specific, suggesting their acquisition in the stem lineage of marsupials after the divergence from eutherians. Despite the time length of the eutherian stem lineage being similar to that of the marsupial lineage, acquisition of testis-expressed genes was not found in the stem lineage of eutherians; rather, their expression patterns differed by species, suggesting rapid gene evolution in the eutherian ancestors. Fifteen testis-expressed genes are therian-specific, and for three of these genes, the evolutionary tempo is markedly faster in eutherians than in marsupials. Our phylogenetic analysis of Rho GTPase-activating protein 28 (ARHGAP28) suggests the adaptive evolution of this gene in the eutherians, probably together with the expression pattern differentiation.

The admixed brushtail possum genome reveals invasion history in New Zealand and novel imprinted genes

Combining genome assembly with population and functional genomics can provide valuable insights to development and evolution, as well as tools for species management. Here, we present a chromosome-level genome assembly of the common brushtail possum (Trichosurus vulpecula), a model marsupial threatened in parts of their native range in Australia, but also a major introduced pest in New Zealand. Functional genomics reveals post-natal activation of chemosensory and metabolic genes, reflecting unique adaptations to altricial birth and delayed weaning, a hallmark of marsupial development. Nuclear and mitochondrial analyses trace New Zealand possums to distinct Australian subspecies, which have subsequently hybridised. This admixture allowed phasing of parental alleles genome-wide, ultimately revealing at least four genes with imprinted, parent-specific expression not yet detected in other species (MLH1, EPM2AIP1, UBP1 and GPX7). We find that reprogramming of possum germline imprints, and the wider epigenome, is similar to eutherian mammals except onset occurs after birth. Together, this work is useful for genetic-based control and conservation of possums, and contributes to understanding of the evolution of novel mammalian epigenetic traits.

Divergent patterns of meiotic double strand breaks and synapsis initiation dynamics suggest an evolutionary shift in the meiosis program between American and Australian marsupials

In eutherian mammals, hundreds of programmed DNA double-strand breaks (DSBs) are generated at the onset of meiosis. The DNA damage response is then triggered. Although the dynamics of this response is well studied in eutherian mammals, recent findings have revealed different patterns of DNA damage signaling and repair in marsupial mammals. To better characterize these differences, here we analyzed synapsis and the chromosomal distribution of meiotic DSBs markers in three different marsupial species (Thylamys elegans, Dromiciops gliorides, and Macropus eugenii) that represent South American and Australian Orders. Our results revealed inter-specific differences in the chromosomal distribution of DNA damage and repair proteins, which were associated with differing synapsis patterns. In the American species T. elegans and D. gliroides, chromosomal ends were conspicuously polarized in a bouquet configuration and synapsis progressed exclusively from the telomeres towards interstitial regions. This was accompanied by sparse H2AX phosphorylation, mainly accumulating at chromosomal ends. Accordingly, RAD51 and RPA were mainly localized at chromosomal ends throughout prophase I in both American marsupials, likely resulting in reduced recombination rates at interstitial positions. In sharp contrast, synapsis initiated at both interstitial and distal chromosomal regions in the Australian representative M. eugenii, the bouquet polarization was incomplete and ephemeral, γH2AX had a broad nuclear distribution, and RAD51 and RPA foci displayed an even chromosomal distribution. Given the basal evolutionary position of T. elegans, it is likely that the meiotic features reported in this species represent an ancestral pattern in marsupials and that a shift in the meiotic program occurred after the split of D. gliroides and the Australian marsupial clade. Our results open intriguing questions about the regulation and homeostasis of meiotic DSBs in marsupials. The low recombination rates observed at the interstitial chromosomal regions in American marsupials can result in the formation of large linkage groups, thus having an impact in the evolution of their genomes.

Monodelphis domestica Induced Pluripotent Stem Cells Reveal Metatherian Pluripotency Architecture

Marsupials have been a powerful comparative model to understand mammalian biology. However, because of the unique characteristics of their embryology, marsupial pluripotency architecture remains to be fully understood, and nobody has succeeded in developing embryonic stem cells (ESCs) from any marsupial species. We have developed an integration-free iPSC reprogramming method and established validated iPSCs from two inbred strains of a marsupial, Monodelphis domestica. The monoiPSCs showed a significant (6181 DE-genes) and highly uniform (r2 [95% CI] = 0.973 ± 0.007) resetting of the cellular transcriptome and were similar to eutherian ESCs and iPSCs in their overall transcriptomic profiles. However, monoiPSCs showed unique regulatory architecture of the core pluripotency transcription factors and were more like marsupial epiblasts. Our results suggest that POU5F1 and the splice-variant-specific expression of POU5F3 synergistically regulate the opossum pluripotency gene network. It is plausible that POU5F1, POU5F3 splice variant XM_016427856.1, and SOX2 form a self-regulatory network. NANOG expression, however, was specific to monoiPSCs and epiblasts. Furthermore, POU5F1 was highly expressed in trophectoderm cells, whereas all other pluripotency transcription factors were significantly downregulated, suggesting that the regulatory architecture of core pluripotency genes of marsupials may be distinct from that of eutherians.

Genetic and genomic architecture in eight strains of the laboratory opossum Monodelphis domestica

The gray short-tailed opossum (Monodelphis domestica) is an established laboratory-bred marsupial model for biomedical research. It is a critical species for comparative genomics research, providing the pivotal phylogenetic outgroup for studies of derived vs ancestral states of genomic/epigenomic characteristics for eutherian mammal lineages. To characterize the current genetic profile of this laboratory marsupial, we examined 79 individuals from eight established laboratory strains. Double digest restriction site-associated DNA sequencing and whole-genome resequencing experiments were performed to investigate the genetic architecture in these strains. A total of 66,640 high-quality single nucleotide polymorphisms (SNPs) were identified. We analyzed SNP density, average heterozygosity, nucleotide diversity, and population differentiation parameter Fst within and between the eight strains. Principal component and population structure analysis clearly resolve the strains at the level of their ancestral founder populations, and the genetic architecture of these strains correctly reflects their breeding history. We confirmed the successful establishment of the first inbred laboratory opossum strain LSD (inbreeding coefficient F > 0.99) and a nearly inbred strain FD2M1 (0.98 < F < 0.99), each derived from a different ancestral background. These strains are suitable for various experimental protocols requiring controlled genetic backgrounds and for intercrosses and backcrosses that can generate offspring with informative SNPs for studying a variety of genetic and epigenetic processes. Together with recent advances in reproductive manipulation and CRISPR/Cas9 techniques for Monodelphis domestica, the existence of distinctive inbred strains will enable genome editing on different genetic backgrounds, greatly expanding the utility of this marsupial model for biomedical research.

Identification and characterization of repetitive DNA in the genus Didelphis Linnaeus, 1758 (Didelphimorphia, Didelphidae) and the use of satellite DNAs as phylogenetic markers

Didelphis species have been shown to exhibit very conservative karyotypes, which mainly differ in their constitutive heterochromatin, known to be mostly composed by repetitive DNAs. In this study, we used genome skimming data combined with computational pipelines to identify the most abundant repetitive DNA families of Lutreolina crassicaudata and all six Didelphis species. We found that transposable elements (TEs), particularly LINE-1, endogenous retroviruses, and SINEs, are the most abundant mobile elements in the studied species. Despite overall similar TE proportions, we report that species of the D. albiventris group consistently present a less diverse TE composition and smaller proportions of LINEs and LTRs in their genomes than other studied species. We also identified four new putative satDNAs (sat206, sat907, sat1430 and sat2324) in the genomes of Didelphis species, which show differences in abundance and nucleotide composition. Phylogenies based on satDNA sequences showed well supported relationships at the species (sat1430) and groups of species (sat206) level, recovering topologies congruent with previous studies. Our study is one of the first attempts to present a characterization of the most abundant families of repetitive DNAs of Lutreolina and Didelphis species providing insights into the repetitive DNA composition in the genome landscape of American marsupials.

A single-cell transcriptome atlas of marsupial embryogenesis and X inactivation

Single-cell RNA sequencing of embryos can resolve the transcriptional landscape of development at unprecedented resolution. To date, single-cell RNA-sequencing studies of mammalian embryos have focused exclusively on eutherian species. Analysis of mammalian outgroups has the potential to identify deeply conserved lineage specification and pluripotency factors, and can extend our understanding of X dosage compensation. Metatherian (marsupial) mammals diverged from eutherians around 160 million years ago. They exhibit distinctive developmental features, including late implantation¹ and imprinted X chromosome inactivation², which is associated with expression of the XIST-like noncoding RNA RSX³. Here we perform a single-cell RNA-sequencing analysis of embryogenesis and X chromosome inactivation in a marsupial, the grey short-tailed opossum (Monodelphis domestica). We resolve the developmental trajectory and transcriptional signatures of the epiblast, primitive endoderm and trophectoderm, and identify deeply conserved lineage-specific markers that pre-date the eutherian-marsupial divergence. RSX coating and inactivation of the X chromosome occurs early and rapidly. This observation supports the hypothesis that-in organisms with early X chromosome inactivation-imprinted X chromosome inactivation prevents biallelic X silencing. We identify XSR, an RSX antisense transcript expressed from the active X chromosome, as a candidate for the regulator of imprinted X chromosome inactivation. Our datasets provide insights into the evolution of mammalian embryogenesis and X dosage compensation.

Marsupial genomics meet marsupial reproduction

We came from very different backgrounds, with different skills and interests. Marilyn Renfree was recognised as 'a giant of marsupial embryology'; I had spent my working life studying genes and chromosomes. We teamed up out of mutual respect (awe on my side) to form, with Des Cooper, the ARC Centre of Excellence in Kangaroo Genomics. This is the story of how our collaboration came to be, and what it has produced for our knowledge of some of the world's most remarkable animals.

Induced Pluripotent Stem Cells from a Marsupial, the Tasmanian Devil (Sarcophilus harrisii): Insight into the Evolution of Mammalian Pluripotency

We demonstrate the generation of Tasmanian devil (Sarcophilus harrisii) induced pluripotent stem cells (DeviPSCs) from dermal fibroblasts by lentiviral delivery of human transcription factors. DeviPSCs display characteristic pluripotent stem cell colony morphology, with individual cells having a high nuclear-to-cytoplasmic ratio and alkaline phosphatase activity. DeviPSCs are leukemia inhibitory factor dependent and have reactivated endogenous octamer-binding transcription factor 4 [OCT4, POU domain, class 5, transcription factor 1 (POU5F1)], POU2 [POU domain, class 5, transcription factor 3 (POU5F3)], sex determining region Y-box 2 (SOX2), Nanog homeobox (NANOG) and dosage-sensitive sex reversal, adrenal hypoplasia congenita critical region on the X chromosome, gene 1 (DAX1) genes, retained a normal karyotype, and concurrently silenced exogenous human transgenes. Notably, co-expression of both OCT4 and POU2 suggests that they are representative of cells of the epiblast, the marsupial equivalent of the inner cell mass. DeviPSCs readily form embryoid bodies and in vitro teratomas containing derivatives of all three embryonic germ layers. To date, DeviPSCs have been stably maintained for more than 45 passages. Our DeviPSCs provide an invaluable resource for studies into marsupial pluripotency and development, and they may also serve as an important tool in efforts to combat the threat of devil facial tumor disease.

Mammalian X Chromosome Dosage Compensation: Perspectives From the Germ Line

Sex chromosomes are advantageous to mammals, allowing them to adopt a genetic rather than environmental sex determination system. However, sex chromosome evolution also carries a burden, because it results in an imbalance in gene dosage between females (XX) and males (XY). This imbalance is resolved by X dosage compensation, which comprises both X chromosome inactivation and X chromosome upregulation. X dosage compensation has been well characterized in the soma, but not in the germ line. Germ cells face a special challenge, because genome wide reprogramming erases epigenetic marks responsible for maintaining the X dosage compensated state. Here we explain how evolution has influenced the gene content and germ line specialization of the mammalian sex chromosomes. We discuss new research uncovering unusual X dosage compensation states in germ cells, which we postulate influence sexual dimorphisms in germ line development and cause infertility in individuals with sex chromosome aneuploidy.

Genomic imprinting, growth and maternal–fetal interactions

In the 1980s, mouse nuclear transplantation experiments revealed that both male and female parental genomes are required for successful development to term (McGrath and Solter, 1983; Surani and Barton, 1983). This non-equivalence of parental genomes is because imprinted genes are predominantly expressed from only one parental chromosome. Uniparental inheritance of these genomic regions causes paediatric growth disorders such as Beckwith–Wiedemann and Silver–Russell syndromes (reviewed in Peters, 2014). More than 100 imprinted genes have now been discovered and the functions of many of these genes have been assessed in murine models. The first such genes described were the fetal growth factor insulin-like growth factor 2 (Igf2) and its inhibitor Igf2 receptor (Igf2r) (DeChiara et al., 1991; Lau et al., 1994; Wang et al., 1994). Since then, it has emerged that most imprinted genes modulate fetal growth and resource acquisition in a variety of ways. First, imprinted genes are required for the development of a functional placenta, the organ that mediates the exchange of nutrients between mother and fetus. Second, these genes act in an embryo-autonomous manner to affect the growth rate and organogenesis. Finally, imprinted genes can signal the nutritional status between mother and fetus, and can modulate levels of maternal care. Importantly, many imprinted genes have been shown to affect postnatal growth and energy homeostasis. Given that abnormal birthweight correlates with adverse adult metabolic health, including obesity and cardiovascular disease, it is crucial to understand how the modulation of this dosage-sensitive, epigenetically regulated class of genes can contribute to fetal and postnatal growth, with implications for lifelong health and disease.

Weird Animals, Sex, and Genome Evolution

Making my career in Australia exposed me to the tyranny of distance, but it gave me opportunities to study our unique native fauna. Distantly related animal species present genetic variation that we can use to explore the most fundamental biological structures and processes. I have compared chromosomes and genomes of kangaroos and platypus, tiger snakes and emus, devils (Tasmanian) and dragons (lizards). I particularly love the challenges posed by sex chromosomes, which, apart from determining sex, provide stunning examples of epigenetic control and break all the evolutionary rules that we currently understand. Here I describe some of those amazing animals and the insights on genome structure, function, and evolution they have afforded us. I also describe my sometimes-random walk in science and the factors and people who influenced my direction. Being a woman in science is still not easy, and I hope others will find encouragement and empathy in my story.

The evolutionary process of mammalian sex determination genes focusing on marsupial SRYs

BACKGROUND

Maleness in mammals is genetically determined by the Y chromosome. On the Y chromosome SRY is known as the mammalian male-determining gene. Both placental mammals (Eutheria) and marsupial mammals (Metatheria) have SRY genes. However, only eutherian SRY genes have been empirically examined by functional analyses, and the involvement of marsupial SRY in male gonad development remains speculative.

RESULTS

In order to demonstrate that the marsupial SRY gene is similar to the eutherian SRY gene in function, we first examined the sequence differences between marsupial and eutherian SRY genes. Then, using a parsimony method, we identify 7 marsupial-specific ancestral substitutions, 13 eutherian-specific ancestral substitutions, and 4 substitutions that occurred at the stem lineage of therian SRY genes. A literature search and molecular dynamics computational simulations support that the lineage-specific ancestral substitutions might be involved with the functional differentiation between marsupial and eutherian SRY genes. To address the function of the marsupial SRY gene in male determination, we performed luciferase assays on the testis enhancer of Sox9 core (TESCO) using the marsupial SRY. The functional assay shows that marsupial SRY gene can weakly up-regulate the luciferase expression via TESCO.

CONCLUSIONS

Despite the sequence differences between the marsupial and eutherian SRY genes, our functional assay indicates that the marsupial SRY gene regulates SOX9 as a transcription factor in a similar way to the eutherian SRY gene. Our results suggest that SRY genes obtained the function of male determination in the common ancestor of Theria (placental mammals and marsupials). This suggests that the marsupial SRY gene has a function in male determination, but additional experiments are needed to be conclusive.

Extensive genetic differentiation detected within a model marsupial, the tammar wallaby (Notamacropus eugenii)

The tammar wallaby (Notamacropus eugenii) is one of the most intensively studied of all macropodids and was the first Australasian marsupial to have its genome sequenced. However, comparatively little is known about genetic diversity and differentiation amongst the morphologically distinct allopatric populations of tammar wallabies found in Western (WA) and South Australia (SA). Here we compare autosomal and Y-linked microsatellite genotypes, as well as sequence data (~600 bp) from the mitochondrial DNA (mtDNA) control region (CR) in tammar wallabies from across its distribution. Levels of diversity at autosomal microsatellite loci were typically high in the WA mainland and Kangaroo Island (SA) populations (A = 8.9-10.6; He = 0.77-0.78) but significantly reduced in other endemic island populations (A = 3.8-4.1; He = 0.41-0.48). Autosomal and Y-linked microsatellite loci revealed a pattern of significant differentiation amongst populations, especially between SA and WA. The Kangaroo Island and introduced New Zealand population showed limited differentiation. Multiple divergent mtDNA CR haplotypes were identified within both SA and WA populations. The CR haplotypes of tammar wallabies from SA and WA show reciprocal monophyly and are highly divergent (14.5%), with levels of sequence divergence more typical of different species. Within WA tammar wallabies, island populations each have unique clusters of highly related CR haplotypes and each is most closely related to different WA mainland haplotypes. Y-linked microsatellite haplotypes show a similar pattern of divergence although levels of diversity are lower. In light of these differences, we suggest that two subspecies of tammar wallaby be recognized; Notamacropus eugenii eugenii in SA and N. eugenii derbianus in WA. The extensive neutral genetic diversity and inter-population differentiation identified within tammar wallabies should further increase the species value and usefulness as a model organism.

Chromosomal Speciation in the Genomics Era: Disentangling Phylogenetic Evolution of Rock-wallabies

The association of chromosome rearrangements (CRs) with speciation is well established, and there is a long history of theory and evidence relating to "chromosomal speciation." Genomic sequencing has the potential to provide new insights into how reorganization of genome structure promotes divergence, and in model systems has demonstrated reduced gene flow in rearranged segments. However, there are limits to what we can understand from a small number of model systems, which each only tell us about one episode of chromosomal speciation. Progressing from patterns of association between chromosome (and genic) change, to understanding processes of speciation requires both comparative studies across diverse systems and integration of genome-scale sequence comparisons with other lines of evidence. Here, we showcase a promising example of chromosomal speciation in a non-model organism, the endemic Australian marsupial genus Petrogale. We present initial phylogenetic results from exon-capture that resolve a history of divergence associated with extensive and repeated CRs. Yet it remains challenging to disentangle gene tree heterogeneity caused by recent divergence and gene flow in this and other such recent radiations. We outline a way forward for better integration of comparative genomic sequence data with evidence from molecular cytogenetics, and analyses of shifts in the recombination landscape and potential disruption of meiotic segregation and epigenetic programming. In all likelihood, CRs impact multiple cellular processes and these effects need to be considered together, along with effects of genic divergence. Understanding the effects of CRs together with genic divergence will require development of more integrative theory and inference methods. Together, new data and analysis tools will combine to shed light on long standing questions of how chromosome and genic divergence promote speciation.

Early Developmental and Evolutionary Origins of Gene Body DNA Methylation Patterns in Mammalian Placentas

Over the last 20-80 million years the mammalian placenta has taken on a variety of morphologies through both divergent and convergent evolution. Recently we have shown that the human placenta genome has a unique epigenetic pattern of large partially methylated domains (PMDs) and highly methylated domains (HMDs) with gene body DNA methylation positively correlating with level of gene expression. In order to determine the evolutionary conservation of DNA methylation patterns and transcriptional regulatory programs in the placenta, we performed a genome-wide methylome (MethylC-seq) analysis of human, rhesus macaque, squirrel monkey, mouse, dog, horse, and cow placentas as well as opossum extraembryonic membrane. We found that, similar to human placenta, mammalian placentas and opossum extraembryonic membrane have globally lower levels of methylation compared to somatic tissues. Higher relative gene body methylation was the conserved feature across all mammalian placentas, despite differences in PMD/HMDs and absolute methylation levels. Specifically, higher methylation over the bodies of genes involved in mitosis, vesicle-mediated transport, protein phosphorylation, and chromatin modification was observed compared with the rest of the genome. As in human placenta, higher methylation is associated with higher gene expression and is predictive of genic location across species. Analysis of DNA methylation in oocytes and preimplantation embryos shows a conserved pattern of gene body methylation similar to the placenta. Intriguingly, mouse and cow oocytes and mouse early embryos have PMD/HMDs but their placentas do not, suggesting that PMD/HMDs are a feature of early preimplantation methylation patterns that become lost during placental development in some species and following implantation of the embryo.

Post-natal imprinting: evidence from marsupials

Genomic imprinting has been identified in therian (eutherian and marsupial) mammals but not in prototherian (monotreme) mammals. Imprinting has an important role in optimising pre-natal nutrition and growth, and most imprinted genes are expressed and imprinted in the placenta and developing fetus. In marsupials, however, the placental attachment is short-lived, and most growth and development occurs post-natally, supported by a changing milk composition tailor-made for each stage of development. Therefore there is a much greater demand on marsupial females during post-natal lactation than during pre-natal placentation, so there may be greater selection for genomic imprinting in the mammary gland than in the short-lived placenta. Recent studies in the tammar wallaby confirm the presence of genomic imprinting in nutrient-regulatory genes in the adult mammary gland. This suggests that imprinting may influence infant post-natal growth via the mammary gland as it does pre-natally via the placenta. Similarly, an increasing number of imprinted genes have been implicated in regulating feeding and nurturing behaviour in both the adult and the developing neonate/offspring in mice. Together these studies provide evidence that genomic imprinting is critical for regulating growth and subsequently the survival of offspring not only pre-natally but also post-natally.