PRKACB is a novel imprinted gene in marsupials.

BACKGROUND: Genomic imprinting results in parent-of-origin-specific gene expression and, among vertebrates, is found only in therian mammals: marsupials and eutherians. A differentially methylated region (DMR), in which the methylation status of CpG dinucleotides differs between the two alleles, can mark the parental identity of imprinted genes. We developed a computational pipeline that detected CpG islands (CGIs) marked by both methylated and unmethylated signals in whole genome bisulfite sequencing data. This approach identified candidate marsupial DMRs in a publicly available koala methylome. One of these candidate DMRs was associated with PRKACB, a gene encoding the protein kinase A catalytic subunit beta. Nothing is known about the imprinting status of PRKACB in eutherian mammals although mutations of this gene are associated with endocrine neoplasia and other developmental disorders. RESULTS: In the tammar wallaby and brushtail possum there was parent-of-origin-specific DNA methylation in the PRKACB DMR in which the maternal allele was methylated and the paternal allele was unmethylated. There were multiple RNAs transcribed from this locus. Allele-specific expression analysis identified paternal expression of a PRKACB lncRNA and an mRNA isoform. Comparison of the PRKACB gene start site between marsupials and eutherians demonstrated that the CGI is longer in marsupials. The PRKACB gene product functions in the same signalling pathway as the guanine nucleotide-binding protein alpha subunit encoded at the GNAS locus, a known eutherian imprinted gene. In a mouse methylome Gnas had three differentially methylated CGIs, while in the koala methylome the GNAS locus had two unmethylated CGIs. CONCLUSIONS: We conclude that PRKACB is a novel, DMR-associated marsupial imprinted gene. Imprinting of PRKACB in marsupials and GNAS in eutherians may indicate a conserved selection pressure for imprinting of the protein kinase A signalling pathway in therians with the two lineages adapting by imprinting different genes.

Imprinted X chromosome inactivation in marsupials: The paternal X arrives at the egg with a silent DNA methylation profile.

X chromosome inactivation (XCI) is an epigenetic process that results in the transcriptional silencing of one X chromosome in the somatic cells of females. This phenomenon is common to both eutherian and marsupial mammals, but there are fundamental differences. In eutherians, the X chosen for silencing is random. DNA methylation on the eutherian inactive X is high at transcription start sites (TSSs) and their flanking regions, resulting in universally high DNA methylation. This contrasts XCI in marsupials where the paternally derived X is always silenced, and in which DNA methylation is low at TSSs and flanking regions. Here, we examined the DNA methylation status of the tammar wallaby X chromosome during spermatogenesis to determine the DNA methylation profile of the paternal X prior to and at fertilization. Whole genome enzymatic methylation sequencing was carried out on enriched flow-sorted populations of premeiotic, meiotic, and postmeiotic cells. We observed that the X displayed a pattern of DNA methylation from spermatogonia to mature sperm that reflected the inactive X in female somatic tissue. Therefore, the paternal X chromosome arrives at the egg with a DNA methylation profile that reflects the transcriptionally silent X in adult female somatic tissue. We present this epigenetic signature as a candidate for the long sought-after imprint for paternal XCI in marsupials.

Coat colour in marsupials: genetic variants at the ASIP locus determine grey and black fur of the brushtail possum.

The possession of fur or hair is a defining characteristic of mammals and can occur in a variety of colours and patterns. While genetic determinants of coat colour are well described in eutherian 'placental' mammals, the other major mammalian infraclass, marsupials, is grossly understudied. The fur of the common brushtail possum (Trichosurus vulpecula), an iconic native mammal found throughout Australia and introduced into Aotearoa New Zealand, possesses two main colour morphs: grey and black. To identify genetic variants associated with coat colour, we performed a genome-wide association study (GWAS) with genotype by sequencing (GBS) data. Single nucleotide variants (SNVs) on chromosome 3, close to the agouti signalling protein (ASIP) gene that controls the temporal and spatial distribution of pigments in eutherian mammals, were identified. Fine-mapping identified a C>T variant at chr3:100483705 that results in a ASIP:p.Arg115Cys missense substitution, and animals homozygous for this variant have black fur. In addition to uncovering the first genetic determinant of coat colour in a natural marsupial population, comparative analysis of ASIP in divergent marsupial species identified the dasyurids as having accelerated evolution, reflecting their well described diversity of coat colour and pattern.

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.

Pronounced sequence specificity of the TET enzyme catalytic domain guides its cellular function.

TET (ten-eleven translocation) enzymes catalyze the oxidation of 5-methylcytosine bases in DNA, thus driving active and passive DNA demethylation. Here, we report that the catalytic domain of mammalian TET enzymes favor CGs embedded within basic helix-loop-helix and basic leucine zipper domain transcription factor-binding sites, with up to 250-fold preference in vitro. Crystal structures and molecular dynamics calculations show that sequence preference is caused by intrasubstrate interactions and CG flanking sequence indirectly affecting enzyme conformation. TET sequence preferences are physiologically relevant as they explain the rates of DNA demethylation in TET-rescue experiments in culture and in vivo within the zygote and germ line. Most and least favorable TET motifs represent DNA sites that are bound by methylation-sensitive immediate-early transcription factors and octamer-binding transcription factor 4 (OCT4), respectively, illuminating TET function in transcriptional responses and pluripotency support.

Stage-specific regulation of signalling pathways to differentiate pluripotent stem cells to cardiomyocytes with ventricular lineage.

BACKGROUND: Pluripotent stem cells (PSCs) can be an ideal source of differentiation of cardiomyocytes in vitro and during transplantation to induce cardiac regeneration. However, differentiation of PSCs into a heterogeneous population is associated with an increased incidence of arrhythmia following transplantation. We aimed to design a protocol to drive PSCs to a ventricular lineage by regulating Wnt and retinoic acid (RA) signalling pathways. METHODS: Mouse embryonic stem cells were cultured either in monolayers or three-dimensional hanging drop method to form embryonic bodies (EBs) and exposed to different treatments acting on Wnt and retinoic acid signalling. Samples were collected at different time points to analyse cardiomyocyte-specific markers by RT-PCR, flow cytometry and immunofluorescence. RESULTS: Treatment of monolayer and EBs with Wnt and RA signalling pathways and ascorbic acid, as a cardiac programming enhancer, resulted in the formation of an immature non-contractile cardiac population that expressed many of the putative markers of cardiac differentiation. The population exhibited upregulation of ventricular specific markers while suppressing the expression of pro-atrial and pro-sinoatrial markers. Differentiation of EBs resulted in early foetal like non-contractile ventricular cardiomyocytes with an inherent propensity to contract when stimulated. CONCLUSION: Our results provide the first evidence of in vitro differentiation that mimics the embryonic morphogenesis towards ventricular specific cardiomyocytes through regulation of Wnt and RA signalling pathways.

Castration delays epigenetic aging and feminizes DNA methylation at androgen-regulated loci.

In mammals, females generally live longer than males. Nevertheless, the mechanisms underpinning sex-dependent longevity are currently unclear. Epigenetic clocks are powerful biological biomarkers capable of precisely estimating chronological age and identifying novel factors influencing the aging rate using only DNA methylation data. In this study, we developed the first epigenetic clock for domesticated sheep (Ovis aries), which can predict chronological age with a median absolute error of 5.1 months. We have discovered that castrated male sheep have a decelerated aging rate compared to intact males, mediated at least in part by the removal of androgens. Furthermore, we identified several androgen-sensitive CpG dinucleotides that become progressively hypomethylated with age in intact males, but remain stable in castrated males and females. Comparable sex-specific methylation differences in MKLN1 also exist in bat skin and a range of mouse tissues that have high androgen receptor expression, indicating that it may drive androgen-dependent hypomethylation in divergent mammalian species. In characterizing these sites, we identify biologically plausible mechanisms explaining how androgens drive male-accelerated aging.

Estimating Global Methylation and Erasure Using Low-Coverage Whole-Genome Bisulfite Sequencing (WGBS ).

Whole-genome bisulfite sequencing (WGBS) is a popular method for characterizing cytosine methylation because it is fully quantitative and has base-pair resolution. While WGBS is prohibitively expensive for experiments involving many samples, low-coverage WGBS can accurately determine global methylation and erasure at similar cost to high-performance liquid chromatography (HPLC) or enzyme-linked immunosorbent assays (ELISA). Moreover, low-coverage WGBS has the capacity to distinguish between methylation in different cytosine contexts (e.g., CG, CHH, and CHG), can tolerate low-input material (<100 cells), and can detect the presence of overrepresented DNA originating from mitochondria or amplified ribosomal DNA. In addition to describing a WGBS library construction and quantitation approach, here we detail computational methods to predict the accuracy of low-coverage WGBS using empirical bootstrap samplers and theoretical estimators similar to those used in election polling. Using examples, we further demonstrate how non-independent sampling of cytosines can alter the precision of error calculation and provide methods to improve this.

Primordial germ cell expression of SSEA1 and DDX4 (VASA) in female Trichosurus vulpecula (Marsupialia) reveals conserved and unique molecular patterns during marsupial germ cell development.

Development of primordial germ cells (PGCs: precursors to adult gametes) is a key process in vertebrate sexual differentiation. Marsupials are ideal to investigate this phenomenon because much of PGC migration and development unusually occurs postnatally in pouch young. However, investigation of the molecular dynamics underpinning PGC development is restricted to one marsupial model species: the tammar wallaby (Macropus eugenii). Given the reproductive diversity among clades, marsupial PGCs likely exhibit diversity in molecular patterns that could help uncover their developmental dynamics. Here we characterise PGC marker expression (SSEA1 and DDX4) in developing ovaries of the brushtail possum, Trichosurus vulpecula. Female germ cells expressed DDX4 from 6 days postpartum (dpp) and almost all germ cells expressed DDX4 by meiosis (40 dpp), consistent with M. eugenii and eutherian mammals. In contrast, PGCs and oogonia expressed SSEA1 from 12 dpp, throughout proliferation and until entry into meiosis (40-63 dpp). SSEA1 expression was temporally distinct from that of M. eugenii, in which SSEA1 expression persists only until 14 dpp, indicating differential expression between marsupial species at equivalent stages of germ cell development. Hence, the molecular characteristics of M. eugenii germ cells cannot be assumed for all marsupials, as at least one key molecule exhibits species-specific expression.

Uncoupling Molecular Testing for SARS-CoV-2 From International Supply Chains.

The rapid global rise of COVID-19 from late 2019 caught major manufacturers of RT-qPCR reagents by surprise and threw into sharp focus the heavy reliance of molecular diagnostic providers on a handful of reagent suppliers. In addition, lockdown and transport bans, necessarily imposed to contain disease spread, put pressure on global supply lines with freight volumes severely restricted. These issues were acutely felt in New Zealand, an island nation located at the end of most supply lines. This led New Zealand scientists to pose the hypothetical question: in a doomsday scenario where access to COVID-19 RT-qPCR reagents became unavailable, would New Zealand possess the expertise and infrastructure to make its own reagents onshore? In this work we describe a review of New Zealand's COVID-19 test requirements, bring together local experts and resources to make all reagents for the RT-qPCR process, and create a COVID-19 diagnostic assay referred to as HomeBrew (HB) RT-qPCR from onshore synthesized components. This one-step RT-qPCR assay was evaluated using clinical samples and shown to be comparable to a commercial COVID-19 assay. Through this work we show New Zealand has both the expertise and, with sufficient lead time and forward planning, infrastructure capacity to meet reagent supply challenges if they were ever to emerge.

Hairpin-bisulfite sequencing of cells exposed to decitabine documents the process of DNA demethylation.

Although the mechanism of DNA demethylating drugs has been understood for many years, the direct effect of these drugs on methylation of the complementary strands of DNA has not been formally demonstrated. By using hairpin-bisulphite sequencing, we describe the kinetics and pattern of DNA methylation following treatment of cells by the DNA methyltransferase 1 (DNMT1) inhibitor, decitabine. As expected, we demonstrate complete loss of methylation on the daughter strand following S-phase in selected densely methylated genes in synchronized Jurkat cells. Thereafter, cells showed a heterogeneous pattern of methylation reflecting replication of the unmethylated strand and restoration of methylation.

DNA methylation study of Huntington's disease and motor progression in patients and in animal models.

Although Huntington's disease (HD) is a well studied Mendelian genetic disorder, less is known about its associated epigenetic changes. Here, we characterize DNA methylation levels in six different tissues from 3 species: a mouse huntingtin (Htt) gene knock-in model, a transgenic HTT sheep model, and humans. Our epigenome-wide association study (EWAS) of human blood reveals that HD mutation status is significantly (p < 10-7) associated with 33 CpG sites, including the HTT gene (p = 6.5 × 10-26). These Htt/HTT associations were replicated in the Q175 Htt knock-in mouse model (p = 6.0 × 10-8) and in the transgenic sheep model (p = 2.4 × 10-88). We define a measure of HD motor score progression among manifest HD cases based on multiple clinical assessments. EWAS of motor progression in manifest HD cases exhibits significant (p < 10-7) associations with methylation levels at three loci: near PEX14 (p = 9.3 × 10-9), GRIK4 (p = 3.0 × 10-8), and COX4I2 (p = 6.5 × 10-8). We conclude that HD is accompanied by profound changes of DNA methylation levels in three mammalian species.

The Genetics and Epigenetics of Sex Change in Fish.

Fish show extraordinary sexual plasticity, changing sex naturally as part of their life cycle or reversing sex because of environmental stressors. This plasticity shows that sexual fate is not an irreversible process but the result of an ongoing tug-of-war for supremacy between male and female signaling networks. The behavioral, gonadal, and morphological changes involved in this process are well described, yet the molecular events that underpin those changes remain poorly understood. Epigenetic modifications emerge as a critical link between environmental stimuli, the onset of sex change, and subsequent maintenance of sexual phenotype. Here we synthesize current knowledge of sex change, focusing on the genetic and epigenetic processes that are likely involved in the initiation and regulation of sex change. We anticipate that better understanding of sex change in fish will shed new light on sex determination and development in vertebrates and on how environmental perturbations affect sexual fate.

Stress, novel sex genes, and epigenetic reprogramming orchestrate socially controlled sex change.

Bluehead wrasses undergo dramatic, socially cued female-to-male sex change. We apply transcriptomic and methylome approaches in this wild coral reef fish to identify the primary trigger and subsequent molecular cascade of gonadal metamorphosis. Our data suggest that the environmental stimulus is exerted via the stress axis and that repression of the aromatase gene (encoding the enzyme converting androgens to estrogens) triggers a cascaded collapse of feminizing gene expression and identifies notable sex-specific gene neofunctionalization. Furthermore, sex change involves distinct epigenetic reprogramming and an intermediate state with altered epigenetic machinery expression akin to the early developmental cells of mammals. These findings reveal at a molecular level how a normally committed developmental process remains plastic and is reversed to completely alter organ structures.

Zebrafish preserve global germline DNA methylation while sex-linked rDNA is amplified and demethylated during feminisation.

The germline is the only cellular lineage capable of transferring genetic information from one generation to the next. Intergenerational transmission of epigenetic memory through the germline, in the form of DNA methylation, has been proposed; however, in mammals this is largely prevented by extensive epigenetic erasure during germline definition. Here we report that, unlike mammals, the continuously-defined 'preformed' germline of zebrafish does not undergo genome-wide erasure of DNA methylation during development. Our analysis also uncovers oocyte-specific germline amplification and demethylation of an 11.5-kb repeat region encoding 45S ribosomal RNA (fem-rDNA). The peak of fem-rDNA amplification coincides with the initial expansion of stage IB oocytes, the poly-nucleolar cell type responsible for zebrafish feminisation. Given that fem-rDNA overlaps with the only zebrafish locus identified thus far as sex-linked, we hypothesise fem-rDNA expansion could be intrinsic to sex determination in this species.

Natural sex change in fish.

Sexual fate can no longer be considered an irreversible deterministic process that once established during early embryonic development, plays out unchanged across an organism's life. Rather, it appears to be a dynamic process, with sexual phenotype determined through an ongoing battle for supremacy between antagonistic male and female developmental pathways. That sexual fate is not final and is actively regulated via the suppression or activation of opposing genetic networks creates the potential for flexibility in sexual phenotype in adulthood. Such flexibility is seen in many fish, where sex change is a usual and adaptive part of the life cycle. Many fish are sequential hermaphrodites, beginning life as one sex and changing sometime later to the other. Sequential hermaphrodites include species capable of female-to-male (protogynous), male-to-female (protandrous), or bidirectional (serial) sex change. These natural forms of sex change involve coordinated transformations across multiple biological systems, including behavioral, anatomical, neuroendocrine and molecular axes. Here we review the biological processes underlying this amazing transformation, focusing particularly on the molecular aspects, where new genomic technologies are beginning to help us understand how sex change is initiated and regulated at the molecular level.

Bio-On-Magnetic-Beads (BOMB): Open platform for high-throughput nucleic acid extraction and manipulation.

Current molecular biology laboratories rely heavily on the purification and manipulation of nucleic acids. Yet, commonly used centrifuge- and column-based protocols require specialised equipment, often use toxic reagents, and are not economically scalable or practical to use in a high-throughput manner. Although it has been known for some time that magnetic beads can provide an elegant answer to these issues, the development of open-source protocols based on beads has been limited. In this article, we provide step-by-step instructions for an easy synthesis of functionalised magnetic beads, and detailed protocols for their use in the high-throughput purification of plasmids, genomic DNA, RNA and total nucleic acid (TNA) from a range of bacterial, animal, plant, environmental and synthetic sources. We also provide a bead-based protocol for bisulfite conversion and size selection of DNA and RNA fragments. Comparison to other methods highlights the capability, versatility, and extreme cost-effectiveness of using magnetic beads. These open-source protocols and the associated webpage (https://bomb.bio) can serve as a platform for further protocol customisation and community engagement.

Simple Synthesis of Functionalized Paramagnetic Beads for Nucleic Acid Purification and Manipulation.

The purification of nucleic acids is one of the most common procedures employed in modern molecular biology laboratories. Typically, commercial column-based protocols are utilized to isolate DNA or RNA from various sources. However, these methods not only require specialized equipment, but are also extremely expensive for high-throughput applications. Although an elegant answer to this issue can be provided by paramagnetic beads, bead-based open-source protocols have been limited in the past. Here, we provide an easy to follow step-by-step manual for the synthesis of paramagnetic beads, as well as their functionalization with either a silica- or a carboxyl-surface that can be used to replace the commercial columns with self-made magnetic beads. Together with a variety of detailed protocols for their use in high-throughput nucleic acids extractions, this bead synthesis method forms the recently published open platform Bio-On-Magnetic-Beads (BOMB), which is available on PLOS Biology ( Oberacker et al., 2019 ). Updated protocols can be found on the associated webpage (https://bomb.bio).

The elephant shark methylome reveals conservation of epigenetic regulation across jawed vertebrates.

BACKGROUND: Methylation of CG dinucleotides constitutes a critical system of epigenetic memory in bony vertebrates, where it modulates gene expression and suppresses transposon activity. The genomes of studied vertebrates are pervasively hypermethylated, with the exception of regulatory elements such as transcription start sites (TSSs), where the presence of methylation is associated with gene silencing. This system is not found in the sparsely methylated genomes of invertebrates, and establishing how it arose during early vertebrate evolution is impeded by a paucity of epigenetic data from basal vertebrates. METHODS:  We perform whole-genome bisulfite sequencing to generate the first genome-wide methylation profiles of a cartilaginous fish, the elephant shark Callorhinchus milii. Employing these to determine the elephant shark methylome structure and its relationship with expression, we compare this with higher vertebrates and an invertebrate chordate using published methylation and transcriptome data.  Results: Like higher vertebrates, the majority of elephant shark CG sites are highly methylated, and methylation is abundant across the genome rather than patterned in the mosaic configuration of invertebrates. This global hypermethylation includes transposable elements and the bodies of genes at all expression levels. Significantly, we document an inverse relationship between TSS methylation and expression in the elephant shark, supporting the presence of the repressive regulatory architecture shared by higher vertebrates. CONCLUSIONS:  Our demonstration that methylation patterns in a cartilaginous fish are characteristic of higher vertebrates imply the conservation of this epigenetic modification system across jawed vertebrates separated by 465 million years of evolution. In addition, these findings position the elephant shark as a valuable model to explore the evolutionary history and function of vertebrate methylation.