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

Cathelicidin antimicrobial peptides mediate immune protection of marsupial neonates

Marsupial neonates are born with immature immune systems, making them vulnerable to pathogens. While neonates receive maternal protection, they can also independently combat pathogens, though the mechanisms remain unknown. Using the sugar glider ( Petaurus breviceps ), we investigated marsupial neonatal immune defenses. Cathelicidin antimicrobial peptides are highly expressed in neonates, with their coordinated expression mediated by enhancer sharing within gene clusters and long-range interactions between clusters. Moreover, cathelicidins modulate immune responses and have potent antimicrobial effects, sufficient to provide protection in a mouse model of sepsis. Lastly, cathelicidins have a complex evolutionary history, where marsupials and monotremes are the only tetrapods that retained two cathelicidin clusters. Thus, cathelicidins are critical mediators of marsupial immunity, and their evolution reflects the life history-specific immunological needs of these animals.

Y Chromosome Loss and Implications for Oncology

The Y chromosome has recognized functions in promoting male sex determination and regulating aspects of fertility. However, recent work has demonstrated important roles for the Y chromosome and Y-encoded genes in multiple domains of male health, including cancer. It is well established that males experience shorter lifespans than females, and this sex bias on overall mortality is accentuated in populations with longer life expectancy, in part related to elevated rates of cancer. The majority of human malignancies exhibit a sex bias with elevated frequencies in males. For many of these cancer types, the disparity has not been explained by environmental risk factors such as tobacco use. Notably, loss of the Y chromosome (LOY) detected in blood cells, termed mosaic LOY, is a common event that is related to advancing age and is associated with a shortened lifespan. Mosaic LOY is linked to increased incidence and mortality across a range of malignancies. Furthermore, tumors arising in different anatomic sites exhibit different frequencies of partial or complete Y chromosome loss. Causal oncogenic or tumor-suppressive roles have been documented for several Y-encoded genes, such as lysine-specific demethylase 5 D, that exert pleiotropic effects on cellular functions by virtue of genome-wide regulation of gene activity. In this review, we discuss aspects of the Y chromosome relevant to oncology. The recent completion of the entire human Y-chromosome sequence provides a reference map of Y-encoded genes and regulatory elements to enable causal molecular studies that may explain and exploit the marked disparity in male cancer risk and mortality.

slc26a12-A novel member of the slc26 family, is located in tandem with slc26a2 in coelacanths, amphibians, reptiles, and birds

Solute carrier family 26 (Slc26) is a family of anion exchangers with 11 members in mammals (named Slc26a1-a11). Here, we identified a novel member of the slc26 family, slc26a12, located in tandem with slc26a2 in the genomes of several vertebrate lineages. BLAST and synteny analyses of various jawed vertebrate genome databases revealed that slc26a12 is present in coelacanths, amphibians, reptiles, and birds but not in cartilaginous fishes, lungfish, mammals, or ray-finned fishes. In some avian and reptilian lineages such as owls, penguins, egrets, and ducks, and most turtles examined, slc26a12 was lost or pseudogenized. Phylogenetic analysis showed that Slc26a12 formed an independent branch with the other Slc26 members and Slc26a12, Slc26a1 and Slc26a2 formed a single branch, suggesting that these three members formed a subfamily in Slc26. In jawless fish, hagfish have two genes homologous to slc26a2 and slc26a12, whereas lamprey has a single gene homologous to slc26a2. African clawed frogs express slc26a12 in larval gills, skin, and fins. These results show that slc26a12 was present at least before the separation of lobe-finned fish and tetrapods; the name slc26a12 is appropriate because the gene duplication occurred in the distant past.

A genomic hotspot of diversifying selection and structural change in the hoary bat (Lasiurus cinereus)

Background

Previous work found that numerous genes positively selected within the hoary bat (Lasiurus cinereus) lineage are physically clustered in regions of conserved synteny. Here I further validate and expand on those finding utilizing an updated L. cinereus genome assembly and additional bat species as well as other tetrapod outgroups.

Methods

A chromosome-level assembly was generated by chromatin-contact mapping and made available by DNAZoo (www.dnazoo.org). The genomic organization of orthologous genes was extracted from annotation data for multiple additional bat species as well as other tetrapod clades for which chromosome-level assemblies were available from the National Center for Biotechnology Information (NCBI). Tests of branch-specific positive selection were performed for L. cinereus using PAML as well as with the HyPhy package for comparison.

Results

Twelve genes exhibiting significant diversifying selection in the L. cinereus lineage were clustered within a 12-Mb genomic window; one of these (Trpc4) also exhibited diversifying selection in bats generally. Ten of the 12 genes are landmarks of two distinct blocks of ancient synteny that are not linked in other tetrapod clades. Bats are further distinguished by frequent structural rearrangements within these synteny blocks, which are rarely observed in other Tetrapoda. Patterns of gene order and orientation among bat taxa are incompatible with phylogeny as presently understood, implying parallel evolution or subsequent reversals. Inferences of positive selection were found to be robust to alternative phylogenetic topologies as well as a strong shift in background nucleotide composition in some taxa.

Discussion

This study confirms and further localizes a genomic hotspot of protein-coding divergence in the hoary bat, one that also exhibits an increased tempo of structural change in bats compared with other mammals. Most genes in the two synteny blocks have elevated expression in brain tissue in humans and model organisms, and genetic studies implicate the selected genes in cranial and neurological development, among other functions.

Genome of the endangered eastern quoll (Dasyurus viverrinus) reveals signatures of historical decline and pelage color evolution

The eastern quoll (Dasyurus viverrinus) is an endangered marsupial native to Australia. Since the extirpation of its mainland populations in the 20th century, wild eastern quolls have been restricted to two islands at the southern end of their historical range. Eastern quolls are the subject of captive breeding programs and attempts have been made to re-establish a population in mainland Australia. However, few resources currently exist to guide the genetic management of this species. Here, we generated a reference genome for the eastern quoll with gene annotations supported by multi-tissue transcriptomes. Our assembly is among the most complete marsupial genomes currently available. Using this assembly, we infer the species' demographic history, identifying potential evidence of a long-term decline beginning in the late Pleistocene. Finally, we identify a deletion at the ASIP locus that likely underpins pelage color differences between the eastern quoll and the closely related Tasmanian devil (Sarcophilus harrisii).

Conserved Noncoding Elements Evolve Around the Same Genes Throughout Metazoan Evolution

Conserved noncoding elements (CNEs) are DNA sequences located outside of protein-coding genes that can remain under purifying selection for up to hundreds of millions of years. Studies in vertebrate genomes have revealed that most CNEs carry out regulatory functions. Notably, many of them are enhancers that control the expression of homeodomain transcription factors and other genes that play crucial roles in embryonic development. To further our knowledge of CNEs in other parts of the animal tree, we conducted a large-scale characterization of CNEs in more than 50 genomes from three of the main branches of the metazoan tree: Cnidaria, Mollusca, and Arthropoda. We identified hundreds of thousands of CNEs and reconstructed the temporal dynamics of their appearance in each lineage, as well as determining their spatial distribution across genomes. We show that CNEs evolve repeatedly around the same genes across the Metazoa, including around homeodomain genes and other transcription factors; they also evolve repeatedly around genes involved in neural development. We also show that transposons are a major source of CNEs, confirming previous observations from vertebrates and suggesting that they have played a major role in wiring developmental gene regulatory mechanisms since the dawn of animal evolution.

Translating the Timing of Developmental Benchmarks in Short-Tailed Opossums (Monodelphisdomestica) to Facilitate Comparisons with Commonly Used Rodent Models

INTRODUCTION

The gray short-tailed opossum, Monodelhis domestica (M. domestica), is a widely used marsupial model species that presents unique advantages for neurodevelopmental studies. Notably their extremely altricial birth allows manipulation of postnatal pups at timepoints equivalent to embryonic stages of placental mammals. A robust literature exists on the development of short-tailed opossums, but many researchers working in the more conventional model species of mice and rats may find it daunting to identify the appropriate age at which to conduct experiments.

METHODS

Here, we present detailed staging diagrams taken from photographic observations of 40 individual pups, in 6 litters, over 25 timepoints across postnatal development. We also present a comparative neurodevelopmental timeline of short-tailed opossums (M. domestica), the house mouse (Mus musculus), and the laboratory rat (Rattus norvegicus) during embryonic as well as postnatal development, using timepoints taken from this study and a review of existing literature, and use this dataset to present statistical models comparing the opossum to the rat and mouse.

RESULTS

One aim of this research was to aid in testing the generalizability of results found in rodents to other mammalian brains, such as the more distantly related metatherians. However, this broad dataset also allows the identification of potential heterochronies in opossum development compared to rats and mice. In contrast to previous work, we found broad similarity between the pace of opossum neural development with that of rats and mice. We also found that development of some systems was accelerated in the opossum, such as the forelimb motor plant, oral motor control, and some aspects of the olfactory system, while the development of the cortex, some aspects of the retina, and other aspects of the olfactory system are delayed compared to the rat and mouse.

DISCUSSION

The pace of opossum development is broadly similar to that of mice and rats, which underscores the usefulness of this species as a compliment to the more commonly used rodents. Many features that differ the most between opossums and rats and mice were either clustered around the day of birth and were features that have functional importance for the pup immediately after or during birth, or were features that have reduced functional importance for the pup until later in postnatal development, given that it is initially attached to the mother.

Expanding the evo-devo toolkit: generation of 3D mammary tissue from diverse mammals

The varying pathways of mammary gland development across species and evolutionary history are underexplored, largely due to a lack of model systems. Recent progress in organoid technology holds the promise of enabling in-depth studies of the developmental adaptations that have occurred throughout the evolution of different species, fostering beneficial phenotypes. The practical application of this technology for mammary glands has been mostly confined to rodents and humans. In the current study, we have successfully created next-generation 3D mammary gland organoids from eight eutherian mammals and the first branched organoid of a marsupial mammary gland. Using mammary organoids, we identified a role for ROCK protein in regulating branching morphogenesis, a role that manifests differently in organoids from different mammals. This finding demonstrates the utility of the 3D organoid model for understanding the evolution and adaptations of signaling pathways. These achievements highlight the potential for organoid models to expand our understanding of mammary gland biology and evolution, and their potential utility in studies of lactation or breast cancer.

Genomic and demographic processes differentially influence genetic variation across the human X chromosome

Many forces influence genetic variation across the genome including mutation, recombination, selection, and demography. Increased mutation and recombination both lead to increases in genetic diversity in a region-specific manner, while complex demographic patterns shape patterns of diversity on a more global scale. While these processes act across the entire genome, the X chromosome is particularly interesting because it contains several distinct regions that are subject to different combinations and strengths of these forces: the pseudoautosomal regions (PARs) and the X-transposed region (XTR). The X chromosome thus can serve as a unique model for studying how genetic and demographic forces act in different contexts to shape patterns of observed variation. We therefore sought to explore diversity, divergence, and linkage disequilibrium in each region of the X chromosome using genomic data from 26 human populations. Across populations, we find that both diversity and substitution rate are consistently elevated in PAR1 and the XTR compared to the rest of the X chromosome. In contrast, linkage disequilibrium is lowest in PAR1, consistent with the high recombination rate in this region, and highest in the region of the X chromosome that does not recombine in males. However, linkage disequilibrium in the XTR is intermediate between PAR1 and the autosomes, and much lower than the non-recombining X. Finally, in addition to these global patterns, we also observed variation in ratios of X versus autosomal diversity consistent with population-specific evolutionary history as well. While our results were generally consistent with previous work, two unexpected observations emerged. First, our results suggest that the XTR does not behave like the rest of the recombining X and may need to be evaluated separately in future studies. Second, the different regions of the X chromosome appear to exhibit unique patterns of linked selection across different human populations. Together, our results highlight profound regional differences across the X chromosome, simultaneously making it an ideal system for exploring the action of evolutionary forces as well as necessitating its careful consideration and treatment in genomic analyses.

Conserved enhancers control notochord expression of vertebrate Brachyury

The cell type-specific expression of key transcription factors is central to development and disease. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three conserved Brachyury-controlling notochord enhancers, T3, C, and I, in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, in cis deletion of all three enhancers in mouse abolishes Brachyury/T/Tbxt expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. The three Brachyury-driving notochord enhancers are conserved beyond mammals in the brachyury/tbxtb loci of fishes, dating their origin to the last common ancestor of jawed vertebrates. Our data define the vertebrate enhancers for Brachyury/T/TBXTB notochord expression through an auto-regulatory mechanism that conveys robustness and adaptability as ancient basis for axis development.

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.

Dynamic Spatial-temporal Expression Ratio of X Chromosome to Autosomes but Stable Dosage Compensation in Mammals

In the evolutionary model of dosage compensation, per-allele expression level of the X chromosome has been proposed to have twofold up-regulation to compensate its dose reduction in males (XY) compared to females (XX). However, the expression regulation of X-linked genes is still controversial, and comprehensive evaluations are still lacking. By integrating multi-omics datasets in mammals, we investigated the expression ratios including X to autosomes (X:AA ratio) and X to orthologs (X:XX ratio) at the transcriptome, translatome, and proteome levels. We revealed a dynamic spatial-temporal X:AA ratio during development in humans and mice. Meanwhile, by tracing the evolution of orthologous gene expression in chickens, platypuses, and opossums, we found a stable expression ratio of X-linked genes in humans to their autosomal orthologs in other species (X:XX ≈ 1) across tissues and developmental stages, demonstrating stable dosage compensation in mammals. We also found that different epigenetic regulations contributed to the high tissue specificity and stage specificity of X-linked gene expression, thus affecting X:AA ratios. It could be concluded that the dynamics of X:AA ratios were attributed to the different gene contents and expression preferences of the X chromosome, rather than the stable dosage compensation.

Conserved enhancer logic controls the notochord expression of vertebrate Brachyury

The cell type-specific expression of key transcription factors is central to development. Brachyury/T/TBXT is a major transcription factor for gastrulation, tailbud patterning, and notochord formation; however, how its expression is controlled in the mammalian notochord has remained elusive. Here, we identify the complement of notochord-specific enhancers in the mammalian Brachyury/T/TBXT gene. Using transgenic assays in zebrafish, axolotl, and mouse, we discover three Brachyury-controlling notochord enhancers T3, C, and I in human, mouse, and marsupial genomes. Acting as Brachyury-responsive, auto-regulatory shadow enhancers, deletion of all three enhancers in mouse abolishes Brachyury/T expression selectively in the notochord, causing specific trunk and neural tube defects without gastrulation or tailbud defects. Sequence and functional conservation of Brachyury-driving notochord enhancers with the brachyury/tbxtb loci from diverse lineages of fishes dates their origin to the last common ancestor of jawed vertebrates. Our data define the enhancers for Brachyury/T/TBXTB notochord expression as ancient mechanism in axis development.

Blood transcriptome responses to PFOA and GenX treatment in the marsupial biomedical model Monodelphis domestica

Introduction: Perfluoroalkyl and poly-fluoroalkyl substances (PFASs) are widely used in industrial and consumer products. Due to their environmental persistence and bioaccumulation, PFASs can be found in the blood of humans and wild animals all over the world. Various fluorinated alternatives such as GenX have been developed to replace the long-chain PFASs, but there is limited information about their potential toxicity. Methods:The current study developed blood culture protocols to assess the response to toxic compounds in the marsupial, Monodelphis domestica. After whole-blood culture conditions were tested and optimized, changes in gene expression in response to PFOA and GenX treatment were assessed. Results: More than 10,000 genes were expressed in the blood transcriptomes with and without treatment. Both PFOA and GenX treatment led to significant changes in the whole blood culture transcriptomes. A total of 578 and 148 differentially expressed genes (DEGs) were detected in the PFOA and GenX treatment groups, 32 of which overlapped. Pathway enrichment analysis revealed that DEGs involved in developmental processes were upregulated after PFOA exposure, while those enriched for metabolic and immune system processes were downregulated. GenX exposure upregulated genes associated with fatty acid transport pathways and inflammatory processes, which is consistent with previous studies using rodent models. Discussion: To our knowledge, this study is the first to investigate the effect of PFASs in a marsupial model. The findings provide supportive evidence for significant transcriptomic alterations, suggesting that this mammalian model may provide a mechanism for exploring the potential toxicity of PFOA and GenX.

Combining vertebrate mitochondrial 12S rRNA gene sequencing and shotgun metagenomic sequencing to investigate the diet of the leopard cat (Prionailurus bengalensis) in Korea

The leopard cat (Prionailurus bengalensis), an endangered species in South Korea, is a small feline widely distributed in Asia. Here, we investigated the diet of leopard cats in the inland areas of Korea by examining their fecal contents using vertebrate mitochondrial 12S rRNA gene sequencing and shotgun metagenomic sequencing. Shotgun metagenomic sequencing revealed that the feces were rich in DNA not only of vertebrates but also of arthropods and plants, but care should be taken when using shotgun metagenomic sequencing to identify vertebrates at low taxonomic levels (e.g., genus level), as it was often erroneous. Meanwhile, vertebrate mitochondrial 12S rRNA gene sequencing was found to be accurate in the genus-level identification, as the genera identified were consistent with the Korean fauna. We found that small mammals such as murids were their main prey. By using these two sequencing methods in combination, this study demonstrated that accurate information about the overall dietary content and vertebrate prey of leopard cats could be obtained. We expect that the continued community efforts to expand the genome database of wildlife, including vertebrates, will alleviate the problem of erroneous identification of prey at low taxonomic levels by shotgun metagenomic sequencing in the near future.

GENESPACE tracks regions of interest and gene copy number variation across multiple genomes

The development of multiple chromosome-scale reference genome sequences in many taxonomic groups has yielded a high-resolution view of the patterns and processes of molecular evolution. Nonetheless, leveraging information across multiple genomes remains a significant challenge in nearly all eukaryotic systems. These challenges range from studying the evolution of chromosome structure, to finding candidate genes for quantitative trait loci, to testing hypotheses about speciation and adaptation. Here, we present GENESPACE, which addresses these challenges by integrating conserved gene order and orthology to define the expected physical position of all genes across multiple genomes. We demonstrate this utility by dissecting presence–absence, copy-number, and structural variation at three levels of biological organization: spanning 300 million years of vertebrate sex chromosome evolution, across the diversity of the Poaceae (grass) plant family, and among 26 maize cultivars. The methods to build and visualize syntenic orthology in the GENESPACE R package offer a significant addition to existing gene family and synteny programs, especially in polyploid, outbred, and other complex genomes.

The genome is the complete DNA sequence of an individual. It is a crucial foundation for many studies in medicine, agriculture, and conservation biology. Advances in genetics have made it possible to rapidly sequence, or read out, the genome of many organisms. For closely related species, scientists can then do detailed comparisons, revealing similar genes with a shared past or a common role, but comparing more distantly related organisms remains difficult.

One major challenge is that genes are often lost or duplicated over evolutionary time. One way to be more confident is to look at ‘synteny’, or how genes are organized or ordered within the genome. In some groups of species, synteny persists across millions of years of evolution. Combining sequence similarity with gene order could make comparisons between distantly related species more robust.

To do this, Lovell et al. developed GENESPACE, a software that links similarities between DNA sequences to the order of genes in a genome. This allows researchers to visualize and explore related DNA sequences and determine whether genes have been lost or duplicated. To demonstrate the value of GENESPACE, Lovell et al. explored evolution in vertebrates and flowering plants. The software was able to highlight the shared sequences between unique sex chromosomes in birds and mammals, and it was able to track the positions of genes important in the evolution of grass crops including maize, wheat, and rice.

Exploring the genetic code in this way could lead to a better understanding of the evolution of important sections of the genome. It might also allow scientists to find target genes for applications like crop improvement. Lovell et al. have designed the GENESPACE software to be easy for other scientists to use, allowing them to make graphics and perform analyses with few programming skills.

Molecular evolution of vision-related genes may contribute to marsupial photic niche adaptations

Vision plays an essential role in the life of many animals. While most mammals are night-active (nocturnal), many have adapted to novel light environments. This includes diurnal (day-active) and crepuscular (twilight-active) species. Here, we used integrative approaches to investigate the molecular evolution of 112 vision-related genes across 19 genomes representing most marsupial orders. We found that four genes (GUCA1B, GUCY2F, RGR, and SWS2) involved in retinal phototransduction likely became functionally redundant in the ancestor of marsupials, a group of largely obligate nocturnal mammals. We also show evidence of rapid evolution and positive selection of bright-light vision genes in the common ancestor of Macropus (kangaroos, wallaroos, and wallabies). Macropus-specific amino acid substitutions in opsin genes (LWS and SWS1), in particular, may be an adaptation for crepuscular vision in this genus via opsin spectral sensitivity tuning. Our study set the stage for functional genetics studies and provides a stepping stone to future research efforts that fully capture the visual repertoire of marsupials.

The GmSNAP11 Contributes to Resistance to Soybean Cyst Nematode Race 4 in Glycine max

Soybean cyst nematode (SCN) has devastating effects on soybean production, making it crucial to identify genes conferring SCN resistance. Here we employed next-generation sequencing-based bulked segregant analysis (BSA) to discover genomic regions, candidate genes, and diagnostic markers for resistance to SCN race 4 (SCN4) in soybean. Phenotypic analysis revealed highly significant differences among the reactions of 145 recombinant inbred lines (RILs) to SCN4. In combination with euclidean distance (ED) and Δsingle-nucleotide polymorphism (SNP)-index analyses, we identified a genomic region on Gm11 (designated as rhg1-paralog) associated with SCN4 resistance. Overexpression and RNA interference analyzes of the two candidate genes identified in this region (GmPLAC8 and GmSNAP11) revealed that only GmSNAP11 significantly contributes to SCN4 resistance. We developed a diagnostic marker for GmSNAP11. Using this marker, together with previously developed markers for SCN-resistant loci, rhg1 and Rhg4, we evaluated the relationship between genotypes and SCN4 resistance in 145 RILs and 30 soybean accessions. The results showed that all the SCN4-resistant lines harbored all the three loci, however, some lines harboring the three loci were still susceptible to SCN4. This suggests that these three loci are necessary for the resistance to SCN4, but they alone cannot confer full resistance. The GmSNAP11 and the diagnostic markers developed could be used in genomic-assisted breeding to develop soybean varieties with increased resistance to SCN4.

Incomplete lineage sorting and phenotypic evolution in marsupials

Incomplete lineage sorting (ILS) makes ancestral genetic polymorphisms persist during rapid speciation events, inducing incongruences between gene trees and species trees. ILS has complicated phylogenetic inference in many lineages, including hominids. However, we lack empirical evidence that ILS leads to incongruent phenotypic variation. Here, we performed phylogenomic analyses to show that the South American monito del monte is the sister lineage of all Australian marsupials, although over 31% of its genome is closer to the Diprotodontia than to other Australian groups due to ILS during ancient radiation. Pervasive conflicting phylogenetic signals across the whole genome are consistent with some of the morphological variation among extant marsupials. We detected hundreds of genes that experienced stochastic fixation during ILS, encoding the same amino acids in non-sister species. Using functional experiments, we confirm how ILS may have directly contributed to hemiplasy in morphological traits that were established during rapid marsupial speciation ca. 60 mya.

Epigenetic clock and methylation studies in marsupials: opossums, Tasmanian devils, kangaroos, and wallabies

The opossum (Monodelphis domestica), with its sequenced genome, ease of laboratory care and experimental manipulation, and unique biology, is the most used laboratory marsupial. Using the mammalian methylation array, we generated DNA methylation data from n = 100 opossum samples from the ear, liver, and tail. We contrasted postnatal development and later aging effects in the opossum methylome with those in mouse (Mus musculus, C57BL/6 J strain) and other marsupial species such as Tasmanian devil, kangaroos, and wallabies. While the opossum methylome is similar to that of mouse during postnatal development, it is distinct from that shared by other mammals when it comes to the age-related gain of methylation at target sites of polycomb repressive complex 2. Our immunohistochemical staining results provide additional support for the hypothesis that PRC2 activity increases with later aging in mouse tissues but remains constant in opossum tissues. We present several epigenetic clocks for opossums that are distinguished by their compatibility with tissue type (pan-tissue and blood clock) and species (opossum and human). Two dual-species human-opossum pan-tissue clocks accurately measure chronological age and relative age, respectively. The human-opossum epigenetic clocks are expected to provide a significant boost to the attractiveness of opossum as a biological model. Additional epigenetic clocks for Tasmanian devil, red kangaroos and other species of the genus Macropus may aid species conservation efforts.

Divergent evolution of developmental timing in the neocortex revealed by marsupial and eutherian transcriptomes

Only mammals evolved a neocortex, which integrates sensory-motor and cognitive functions. Significant diversifications in the cellular composition and connectivity of the neocortex occurred between the two main therian groups: marsupials and eutherians. However, the developmental mechanisms underlying these diversifications are largely unknown. Here, we compared the neocortical transcriptomes of Sminthopsis crassicaudata, a mouse-sized marsupial, with those of eutherian mice at two developmentally equivalent time points corresponding to deeper and upper layer neuron generation. Enrichment analyses revealed more mature gene networks in marsupials at the early stage, which reverted at the later stage, suggesting a more precocious but protracted neuronal maturation program relative to birth timing of cortical layers. We ranked genes expressed in different species and identified important differences in gene expression rankings between species. For example, genes known to be enriched in upper-layer cortical projection neuron subtypes, such as Cux1, Lhx2 and Satb2, likely relate to corpus callosum emergence in eutherians. These results show molecular heterochronies of neocortical development in Theria, and highlight changes in gene expression and cell type composition that may underlie neocortical evolution and diversification.

This article has an associated ‘The people behind the papers’ interview.

Why sequence all eukaryotes?

Life on Earth has evolved from initial simplicity to the astounding complexity we experience today. Bacteria and archaea have largely excelled in metabolic diversification, but eukaryotes additionally display abundant morphological innovation. How have these innovations come about and what constraints are there on the origins of novelty and the continuing maintenance of biodiversity on Earth? The history of life and the code for the working parts of cells and systems are written in the genome. The Earth BioGenome Project has proposed that the genomes of all extant, named eukaryotes-about 2 million species-should be sequenced to high quality to produce a digital library of life on Earth, beginning with strategic phylogenetic, ecological, and high-impact priorities. Here we discuss why we should sequence all eukaryotic species, not just a representative few scattered across the many branches of the tree of life. We suggest that many questions of evolutionary and ecological significance will only be addressable when whole-genome data representing divergences at all of the branchings in the tree of life or all species in natural ecosystems are available. We envisage that a genomic tree of life will foster understanding of the ongoing processes of speciation, adaptation, and organismal dependencies within entire ecosystems. These explorations will resolve long-standing problems in phylogenetics, evolution, ecology, conservation, agriculture, bioindustry, and medicine.

Mammalian tumor-like organs. 1. The role of tumor-like normal organs and atypical tumor organs in the evolution of development (carcino-evo-devo)

Background

Earlier I hypothesized that hereditary tumors might participate in the evolution of multicellular organisms. I formulated the hypothesis of evolution by tumor neofunctionalization, which suggested that the evolutionary role of hereditary tumors might consist in supplying evolving multicellular organisms with extra cell masses for the expression of evolutionarily novel genes and the origin of new cell types, tissues, and organs. A new theory—the carcino-evo-devo theory—has been developed based on this hypothesis.

Main text

My lab has confirmed several non-trivial predictions of this theory. Another non-trivial prediction is that evolutionarily new organs if they originated from hereditary tumors or tumor-like structures, should recapitulate some tumor features in their development. This paper reviews the tumor-like features of evolutionarily novel organs. It turns out that evolutionarily new organs such as the eutherian placenta, mammary gland, prostate, the infantile human brain, and hoods of goldfishes indeed have many features of tumors. I suggested calling normal organs, which have many tumor features, the tumor-like organs.

Conclusion

Tumor-like organs might originate from hereditary atypical tumor organs and represent the part of carcino-evo-devo relationships, i.e., coevolution of normal and neoplastic development. During subsequent evolution, tumor-like organs may lose the features of tumors and the high incidence of cancer and become normal organs without (or with almost no) tumor features.

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.

Rates of Protein Evolution across the Marsupial Phylogeny: Heterogeneity and Link to Life-History Traits

Despite the importance of effective population size (Ne) in evolutionary and conservation biology, it remains unclear what factors have an impact on this quantity. The Nearly Neutral Theory of Molecular Evolution predicts a faster accumulation of deleterious mutations (and thus a higher dN/dS ratio) in populations with small Ne; thus, measuring dN/dS ratios in different groups/species can provide insight into their Ne. Here, we used an exome data set of 1,550 loci from 45 species of marsupials representing 18 of the 22 extant families, to estimate dN/dS ratios across the different branches and families of the marsupial phylogeny. We found a considerable heterogeneity in dN/dS ratios among families and species, which suggests significant differences in their Ne. Furthermore, our multivariate analyses of several life-history traits showed that dN/dS ratios (and thus Ne) are affected by body weight, body length, and weaning age.

Between the Devil and the Deep Blue Sea: Non-Coding RNAs Associated with Transmissible Cancers in Tasmanian Devil, Domestic Dog and Bivalves

Currently there are nine known examples of transmissible cancers in nature. They have been observed in domestic dog, Tasmanian devil, and six bivalve species. These tumours can overcome host immune defences and spread to other members of the same species. Non-coding RNAs (ncRNAs) are known to play roles in tumorigenesis and immune system evasion. Despite their potential importance in transmissible cancers, there have been no studies on ncRNA function in this context to date. Here, we present possible applications of the CRISPR/Cas system to study the RNA biology of transmissible cancers. Specifically, we explore how ncRNAs may play a role in the immortality and immune evasion ability of these tumours.

Developmental and comparative immunology single-cell transcriptome analysis of the B-cell repertoire reveals the usage of immunoglobulins in the gray short-tailed opossum (Monodelphis domestica)

B-cells are key to humoral immunity, are found in multiple lymphoid organs, and have the unique ability to mediate the production of antigen-specific antibodies in the presence of pathogens. The marsupial immunoglobulin (Ig) heavy (H) chain locus encodes four constant region isotypes, IgA, IgG, IgM and IgE, but no IgD, and there are two light (L) chain isotypes, lambda (Igλ) and kappa (Igκ). To gain an understanding of the marsupial humoral immune system, B-cell transcriptomes generated by single-cell RNA sequencing from gray short-tailed opossum (Monodelphis domestica) splenocytes, and peripheral blood mononuclear cells were analysed. The cells used were from a single unimmunized animal and the majority of B-cells were transcribing IgM heavy chains. The ratio of Ig light chain use was roughly 2:1, Igλ:Igκ in this individual. This was not predicted due to Igκ being the more complex of the two L chain loci. The variable (V) gene segment pairs used in individual B-cells confirm greater diversity provided by the L chain V. This study is the first to report on using single cell analysis to investigate Ig repertoires in a marsupial and confirms a number of prior hypothesis, as well as revealing some surprises.

Factors Regulating the Activity of LINE1 Retrotransposons

LINE-1 (L1) is a class of autonomous mobile genetic elements that form somatic mosaicisms in various tissues of the organism. The activity of L1 retrotransposons is strictly controlled by many factors in somatic and germ cells at all stages of ontogenesis. Alteration of L1 activity was noted in a number of diseases: in neuropsychiatric and autoimmune diseases, as well as in various forms of cancer. Altered activity of L1 retrotransposons for some pathologies is associated with epigenetic changes and defects in the genes involved in their repression. This review discusses the molecular genetic mechanisms of the retrotransposition and regulation of the activity of L1 elements. The contribution of various factors controlling the expression and distribution of L1 elements in the genome occurs at all stages of the retrotransposition. The regulation of L1 elements at the transcriptional, post-transcriptional and integration into the genome stages is described in detail. Finally, this review also focuses on the evolutionary aspects of L1 accumulation and their interplay with the host regulation system.

Seq Your Destiny: Neural Crest Fate Determination in the Genomic Era

Neural crest stem/progenitor cells arise early during vertebrate embryogenesis at the border of the forming central nervous system. They subsequently migrate throughout the body, eventually differentiating into diverse cell types ranging from neurons and glia of the peripheral nervous system to bones of the face, portions of the heart, and pigmentation of the skin. Along the body axis, the neural crest is heterogeneous, with different subpopulations arising in the head, neck, trunk, and tail regions, each characterized by distinct migratory patterns and developmental potential. Modern genomic approaches like single-cell RNA- and ATAC-sequencing (seq) have greatly enhanced our understanding of cell lineage trajectories and gene regulatory circuitry underlying the developmental progression of neural crest cells. Here, we discuss how genomic approaches have provided new insights into old questions in neural crest biology by elucidating transcriptional and posttranscriptional mechanisms that govern neural crest formation and the establishment of axial level identity. Expected final online publication date for the Annual Review of Genetics, Volume 55 is November 2021. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Circular RNA repertoires are associated with evolutionarily young transposable elements

Circular RNAs (circRNAs) are found across eukaryotes and can function in post-transcriptional gene regulation. Their biogenesis through a circle-forming backsplicing reaction is facilitated by reverse-complementary repetitive sequences promoting pre-mRNA folding. Orthologous genes from which circRNAs arise, overall contain more strongly conserved splice sites and exons than other genes, yet it remains unclear to what extent this conservation reflects purifying selection acting on the circRNAs themselves. Our analyses of circRNA repertoires from five species representing three mammalian lineages (marsupials, eutherians: rodents, primates) reveal that surprisingly few circRNAs arise from orthologous exonic loci across all species. Even the circRNAs from orthologous loci are associated with young, recently active and species-specific transposable elements, rather than with common, ancient transposon integration events. These observations suggest that many circRNAs emerged convergently during evolution - as a byproduct of splicing in orthologs prone to transposon insertion. Overall, our findings argue against widespread functional circRNA conservation.

Ancient viral integrations in marsupials: a potential antiviral defence

Marsupial viruses are understudied compared to their eutherian mammal counterparts, although they may pose severe threats to vulnerable marsupial populations. Genomic viral integrations, termed 'endogenous viral elements' (EVEs), could protect the host from infection. It is widely known past viral infections and EVEs play an active role in antiviral defence in invertebrates and plants. This study aimed to characterise actively transcribed EVEs in Australian marsupial species, because they may play an integral role in cellular defence against viruses. This study screened publicly available RNA sequencing data sets (n = 35) and characterised 200 viral transcripts from thirteen Australian marsupial species. Of the 200 transcripts, 188 originated from either Bornaviridae, Filoviridae, or Parvoviridae EVEs. The other twelve transcripts were from putative active infections from members of the Herpesviridae and Anelloviridae, and Hepadnaviridae. EVE transcripts (n = 188) were mapped to marsupial genomes (where available, n = 5/13) to identify the genomic insertion sites. Of the 188 transcripts, 117 mapped to 39 EVEs within the koala, bare-nosed wombat, tammar wallaby, brushtail possum, and Tasmanian devil genomes. The remaining eight animals had no available genome (transcripts n = 71). Every marsupial has Bornaviridae, Filoviridae, and Parvoviridae EVEs, a trend widely observed in eutherian mammals. Whilst eutherian bornavirus EVEs are predominantly nucleoprotein-derived, marsupial bornavirus EVEs demonstrate a surprising replicase gene bias. We predicted these widely distributed EVEs were conserved within marsupials from ancient germline integrations, as many were over 65 million years old. One bornavirus replicase EVE, present in six marsupial genomes, was estimated to be 160 million years old, predating the American-Australian marsupial split. We considered transcription of these EVEs through small non-coding RNA as an ancient viral defence. Consistent with this, in koala small RNA sequence data sets, we detected Bornaviridae replicase and Filoviridae nucleoprotein produced small RNA. These were enriched in testis tissue, suggesting they could protect marsupials from vertically transmitted viral integrations.

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.

Targeted gene disruption in a marsupial, Monodelphis domestica, by CRISPR/Cas9 genome editing

Marsupials represent one of three extant mammalian subclasses with very unique characteristics not shared by other mammals. Most notably, much of the development of neonates immaturely born after a relatively short gestation takes place in the external environment. Among marsupials, the gray short-tailed opossum (Monodelphis domestica; hereafter "the opossum") is one of very few established laboratory models. Due to many biologically unique characteristics and experimentally advantageous features, the opossum is used as a prototype species for basic research on marsupial biology.^1,2 However, in vivo studies of gene function in the opossum, and thus marsupials in general, lag far behind those of eutherian mammals due to the lack of reliable means to manipulate their genomes. In this study, we describe the successful generation of genome edited opossums by a combination of refined methodologies in reproductive biology and embryo manipulation. We took advantage of the opossum's resemblance to popular rodent models, such as the mouse and rat, in body size and breeding characteristics. First, we established a tractable pipeline of reproductive technologies, from induction of ovulation, timed copulation, and zygote collection to embryo transfer to pseudopregnant females, that warrant an essential platform to manipulate opossum zygotes. Further, we successfully demonstrated the generation of gene knockout opossums at the Tyr locus by microinjection of pronuclear stage zygotes using CRISPR/Cas9 genome editing, along with germline transmission of the edited alleles to the F1 generation. This study provides a critical foundation for venues to expand mammalian reverse genetics into the metatherian subclass.

Evolution of mouse circadian enhancers from transposable elements

BACKGROUND

Transposable elements are increasingly recognized as a source of cis-regulatory variation. Previous studies have revealed that transposons are often bound by transcription factors and some have been co-opted into functional enhancers regulating host gene expression. However, the process by which transposons mature into complex regulatory elements, like enhancers, remains poorly understood. To investigate this process, we examined the contribution of transposons to the cis-regulatory network controlling circadian gene expression in the mouse liver, a well-characterized network serving an important physiological function.

RESULTS

ChIP-seq analyses reveal that transposons and other repeats contribute ~ 14% of the binding sites for core circadian regulators (CRs) including BMAL1, CLOCK, PER1/2, and CRY1/2, in the mouse liver. RSINE1, an abundant murine-specific SINE, is the only transposon family enriched for CR binding sites across all datasets. Sequence analyses and reporter assays reveal that the circadian regulatory activity of RSINE1 stems from the presence of imperfect CR binding motifs in the ancestral RSINE1 sequence. These motifs matured into canonical motifs through point mutations after transposition. Furthermore, maturation occurred preferentially within elements inserted in the proximity of ancestral CR binding sites. RSINE1 also acquired motifs that recruit nuclear receptors known to cooperate with CRs to regulate circadian gene expression specifically in the liver.

CONCLUSIONS

Our results suggest that the birth of enhancers from transposons is predicated both by the sequence of the transposon and by the cis-regulatory landscape surrounding their genomic integration site.

Contemporary and historical selection in Tasmanian devils (Sarcophilus harrisii) support novel, polygenic response to transmissible cancer

Tasmanian devils (Sarcophilus harrisii) are evolving in response to a unique transmissible cancer, devil facial tumour disease (DFTD), first described in 1996. Persistence of wild populations and the recent emergence of a second independently evolved transmissible cancer suggest that transmissible cancers may be a recurrent feature in devils. Here, we compared signatures of selection across temporal scales to determine whether genes or gene pathways under contemporary selection (six to eight generations) have also been subject to historical selection (65-85 Myr). First, we used targeted sequencing, RAD-capture, in approximately 2500 devils in six populations to identify genomic regions subject to rapid evolution. We documented genome-wide contemporary evolution, including 186 candidate genes related to cell cycling and immune response. Then we used a molecular evolution approach to identify historical positive selection in devils compared to other marsupials and found evidence of selection in 1773 genes. However, we found limited overlap across time scales, with only 16 shared candidate genes, and no overlap in enriched functional gene sets. Our results are consistent with a novel, multi-locus evolutionary response of devils to DFTD. Our results can inform conservation by identifying high priority targets for genetic monitoring and guiding maintenance of adaptive potential in managed populations.

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.

Large vs small genomes in Passiflora: the influence of the mobilome and the satellitome

While two lineages of retrotransposons were more abundant in larger Passiflora genomes, the satellitome was more diverse and abundant in the smallest genome analysed. Repetitive sequences are ubiquitous and fast-evolving elements responsible for size variation and large-scale organization of plant genomes. Within Passiflora genus, a tenfold variation in genome size, not attributed to polyploidy, is known. Here, we applied a combined in silico and cytological approach to study the organization and diversification of repetitive elements in three species of this genus representing its known range in genome size variation. Sequences were classified in terms of type and repetitiveness and the most abundant were mapped to chromosomes. We identified long terminal repeat (LTR) retrotransposons as the most abundant elements in the three genomes, showing a considerable variation among species. Satellite DNAs (satDNAs) were less representative, but highly diverse between subgenera. Our results clearly confirm that the largest genome species (Passiflora quadrangularis) presents a higher accumulation of repetitive DNA sequences, specially Angela and Tekay elements, making up most of its genome. Passiflora cincinnata, with intermediate genome and from the same subgenus, showed similarity with P. quadrangularis regarding the families of repetitive DNA sequences, but in different proportions. On the other hand, Passiflora organensis, the smallest genome, from a different subgenus, presented greater diversity and the highest proportion of satDNA. Altogether, our data indicates that while large genomes evolved by an accumulation of retrotransposons, the smallest genome known for the genus has evolved by diversification of different repeat types, particularly satDNAs.

New Insights into the Mammalian Egg Zona Pellucida

Mammalian oocytes are surrounded by an extracellular coat called the zona pellucida (ZP), which, from an evolutionary point of view, is the most ancient of the coats that envelope vertebrate oocytes and conceptuses. This matrix separates the oocyte from cumulus cells and is responsible for species-specific recognition between gametes, preventing polyspermy and protecting the preimplantation embryo. The ZP is a dynamic structure that shows different properties before and after fertilization. Until very recently, mammalian ZP was believed to be composed of only three glycoproteins, ZP1, ZP2 and ZP3, as first described in mouse. However, studies have revealed that this composition is not necessarily applicable to other mammals. Such differences can be explained by an analysis of the molecular evolution of the ZP gene family, during which ZP genes have suffered pseudogenization and duplication events that have resulted in differing models of ZP protein composition. The many discoveries made in recent years related to ZP composition and evolution suggest that a compilation would be useful. Moreover, this review analyses ZP biosynthesis, the role of each ZP protein in different mammalian species and how these proteins may interact among themselves and with other proteins present in the oviductal lumen.

On the importance of evolutionary constraint for regulatory sequence identification

Regulation of gene expression relies on the activity of specialized genomic elements, enhancers or silencers, distributed over sometimes large distance from their target gene promoters. A significant part of vertebrate genomes consists in such regulatory elements, but their identification and that of their target genes remains challenging, due to the lack of clear signature at the nucleotide level. For many years the main hallmark used for identifying functional elements has been their sequence conservation between genomes of distant species, indicative of purifying selection. More recently, genome-wide biochemical assays have opened new avenues for detecting regulatory regions, shifting attention away from evolutionary constraints. Here, we review the respective contributions of comparative genomics and biochemical assays for the definition of regulatory elements and their targets and advocate that both sequence conservation and preserved synteny, taken as signature of functional constraint, remain essential tools in this task.

Establishment of Long-Term Primary Cortical Neuronal Cultures From Neonatal Opossum Monodelphis domestica

Primary dissociated neuronal cultures have become a standard model for studying central nervous system (CNS) development. Such cultures are predominantly prepared from the hippocampus or cortex of rodents (mice and rats), while other mammals are less used. Here, we describe the establishment and extensive characterization of the primary dissociated neuronal cultures derived from the cortex of the gray South American short-tailed opossums, Monodelphis domestica. Opossums are unique in their ability to fully regenerate their CNS after an injury during their early postnatal development. Thus, we used cortex of postnatal day (P) 3-5 opossum to establish long-surviving and nearly pure neuronal cultures, as well as mixed cultures composed of radial glia cells (RGCs) in which their neurogenic and gliogenic potential was confirmed. Both types of cultures can survive for more than 1 month in vitro. We also prepared neuronal cultures from the P16-18 opossum cortex, which were composed of astrocytes and microglia, in addition to neurons. The long-surviving opossum primary dissociated neuronal cultures represent a novel mammalian in vitro platform particularly useful to study CNS development and regeneration.

Annotation of immune genes in the extinct thylacine (Thylacinus cynocephalus)

Advances in genome sequencing technology have enabled genomes of extinct species to be sequenced. However, given the fragmented nature of these genome assemblies, it is not clear whether it is possible to comprehensively annotate highly variable and repetitive genes such as those involved in immunity. As such, immune genes have only been investigated in a handful of extinct genomes, mainly in human lineages. In 2018 the genome of the thylacine (Thylacinus cynocephalus), a carnivorous marsupial from Tasmania that went extinct in 1936, was sequenced. Here we attempt to characterise the immune repertoire of the thylacine and determine similarity to its closest relative with a genome available, the Tasmanian devil (Sarcophilus harrisii), as well as other marsupials. Members from all major immune gene families were identified. However, variable regions could not be characterised, and complex families such as the major histocompatibility complex (MHC) were highly fragmented and located across multiple small scaffolds. As such, at a gene level we were unable to reconstruct full-length coding sequences for the majority of thylacine immune genes. Despite this, we identified genes encoding functionally important receptors and immune effector molecules, which suggests the functional capacity of the thylacine immune system was similar to other mammals. However, the high number of partial immune gene sequences identified limits our ability to reconstruct an accurate picture of the thylacine immune repertoire.

The developmental hourglass model and recapitulation: An attempt to integrate the two models

Recapitulation is a hypothetical concept that assumes embryogenesis of an animal parallels its own phylogenetic history, sequentially developing from more ancestral features to more derived ones. This concept predicts that the earliest developmental stage of various animals should represent the most evolutionarily conserved patterns. Recent transcriptome‐based studies, on the other hand, have reported that mid‐embryonic, organogenetic periods show the highest level of conservation (the developmental hourglass model). This, however, does not rule out the possibility that recapitulation would still be detected after the mid‐embryonic period. In accordance with this, recapitulation‐like morphological features are enriched in late developmental stages. Moreover, our recent chromatin accessibility‐based study provided molecular evidence for recapitulation in the mid‐to‐late embryogenesis of vertebrates, as newly evolved gene regulatory elements tended to be activated at late embryonic stages. In this review, we revisit the recapitulation hypothesis, together with recent molecular‐based studies that support the developmental hourglass model. We contend that the recapitulation hypothesis does not entirely contradict the developmental hourglass model and that these two may even coexist in later embryonic stages of vertebrates. Finally, we review possible mechanisms underlying the recapitulation pattern of chromatin accessibility together with the hourglass‐like evolutionary conservation in vertebrate embryogenesis.

Recapitulation pattern has been reported for chromatin accessibility during the mid‐to‐late embryogenesis.

The observed recapitulation pattern and the developmental hourglass model may coexist.

The possible evolutionary mechanisms underlying tendencies of embryonic evolution were discussed.

Potential APOBEC-mediated RNA editing of the genomes of SARS-CoV-2 and other coronaviruses and its impact on their longer term evolution

Members of the APOBEC family of cytidine deaminases show antiviral activities in mammalian cells through lethal editing in the genomes of small DNA viruses, herpesviruses and retroviruses, and potentially those of RNA viruses such as coronaviruses. Consistent with the latter, APOBEC-like directional C→U transitions of genomic plus-strand RNA are greatly overrepresented in SARS-CoV-2 genome sequences of variants emerging during the COVID-19 pandemic. A C→U mutational process may leave evolutionary imprints on coronavirus genomes, including extensive homoplasy from editing and reversion at targeted sites and the occurrence of driven amino acid sequence changes in viral proteins. If sustained over longer periods, this process may account for the previously reported marked global depletion of C and excess of U bases in human seasonal coronavirus genomes. This review synthesizes the current knowledge on APOBEC evolution and function and the evidence of their role in APOBEC-mediated genome editing of SARS-CoV-2 and other coronaviruses.

-

SARS-CoV-2 sequence variants contain an overabundance of C- > U transitions

-

C- > U transitions are the hallmark of the activity of APOBEC cytosine deaminases

-

Further work is needed to determine APOBEC's role in coronavirus evolution

Red-tailed phascogales: A review of their biology and importance as model marsupial species

There are many limitations when using traditional laboratory species. Limits on variation, may result in limited outcomes, at both the species and individual level, due to different individuals/species having diverse physiological processes, or differing molecular and genetic mechanisms. By using a variety of model species, we will be able to develop creative solutions to biological problems and identify differences of which we were not previously aware. The laboratory mouse has been a suitable model species for various mammalian studies, however most are bred specifically for laboratory research with limited variability due to selective breeding. Marsupial models offer unique research opportunities compared to eutherian models. We believe that there should be an expansion in marsupial model species, and the introduction of the red-tailed phascogale (Phascogale calura), a dasyurid marsupial, should be one of them. Phascogales are easily managed in captivity, and there are now multiple studies involving their development, reproduction, nutrition, behavior and immune system, which can serve as a baseline for future studies. The addition of the phascogale as a model species will improve future mammalian studies by introducing variability and offer alternate solutions to biological problems, particularly in the areas of genetics, nutrition, immunology, the neuro-endocrine system, and ageing, due to their semelparous reproductive strategy and hence, subsequent predictive physiology. In this review, we provide information based on existing research on red-tailed phascogales to support their inclusion as a model species.

Marsupial chromosomics: bridging the gap between genomes and chromosomes

Marsupials have unique features that make them particularly interesting to study, and sequencing of marsupial genomes is helping to understand their evolution. A decade ago, it was a huge feat to sequence the first marsupial genome. Now, the advances in sequencing technology have made the sequencing of many more marsupial genomes possible. However, the DNA sequence is only one component of the structures it is packaged into: chromosomes. Knowing the arrangement of the DNA sequence on each chromosome is essential for a genome assembly to be used to its full potential. The importance of combining sequence information with cytogenetics has previously been demonstrated for rapidly evolving regions of the genome, such as the sex chromosomes, as well as for reconstructing the ancestral marsupial karyotype and understanding the chromosome rearrangements involved in the Tasmanian devil facial tumour disease. Despite the recent advances in sequencing technology assisting in genome assembly, physical anchoring of the sequence to chromosomes is required to achieve a chromosome-level assembly. Once chromosome-level assemblies are achieved for more marsupials, we will be able to investigate changes in the packaging and interactions between chromosomes to gain an understanding of the role genome architecture has played during marsupial evolution.

Opossum milk IgG is from maternal circulation and timing of transfer correlates with neonatal immune development

Marsupials, with short gestation times, have more complex and changing patterns of milk composition than eutherians. Maternal immunoglobulins (Ig) that confer immunity on offspring are among the components that change during marsupial lactation. In the present study we quantified the abundance of mammary transcripts encoding Ig heavy chains and their corresponding transporters in the laboratory opossum Monodelphis domestica. IgA transcripts were the most abundant in opossum mammary and, with IgM, increased in abundance linearly from birth to weaning. Similarly, the Fc receptor for IgA, the poly-Ig receptor, also increased in abundance throughout lactation. There were few transcripts for IgG or IgE within the opossum mammaries. This is in contrast with reports for Australian marsupial species. Transcripts for the Neonatal Fc Receptor (FcRN), which transports IgG, were detected throughout lactation, and opossum milk is known to contain IgG. Therefore, milk IgG is likely to be taken from the maternal circulation, rather than resulting from local production. There is a parallel increase in FcRN in the newborn gut that declines around the time when neonates have matured to the point where they can make their own IgG. These results are consistent with a transfer of maternal Ig that is coordinated with the development of the neonatal immune system.