Genome-wide analysis in the mouse embryo reveals the importance of DNA methylation for transcription integrity

LncRNA FPASL suppresses fibroblast proliferation through its DNA methylation via DNMT3b in hypertrophic scar

Long noncoding RNAs (lncRNAs) are increasingly being implicated as key regulators of cell proliferation, apoptosis, and differentiation. However, the molecular mechanisms of specific lncRNAs in the context of hypertrophic scar remain largely unclear. Here, we find that the lncRNA FPASL (fibroblast proliferation-associated LncRNA) is downregulated in HS, and FPASL reduces fibroblast proliferation and colony formation and blocks cell cycle progression. Using GO annotation enrichment analysis along with AZC (a specific inhibitor of DNA methylation), we identify that DNA methylation is responsible for downregulating FPASL in hypertrophic scar. Subsequent studies demonstrate that high expression of DNMT3b inhibits FPASL expression in HS. Mechanistic study reveals a significant increase in fibroblast proliferation after transfection with LNA-FPASL, which is further inhibited by knockdown of DNMT3b. Thus, our study reveals that DNMT3b mediates hypermethylation of the lncRNA FPASL promoter and the downregulation of lncRNA FPASL promotes fibroblast proliferation in hypertrophic scar.

Mammalian DNA methylome dynamics: mechanisms, functions and new frontiers

DNA methylation is a highly conserved epigenetic modification that plays essential roles in mammalian gene regulation, genome stability and development. Despite being primarily considered a stable and heritable epigenetic silencing mechanism at heterochromatic and repetitive regions, whole genome methylome analysis reveals that DNA methylation can be highly cell-type specific and dynamic within proximal and distal gene regulatory elements during early embryonic development, stem cell differentiation and reprogramming, and tissue maturation. In this Review, we focus on the mechanisms and functions of regulated DNA methylation and demethylation, highlighting how these dynamics, together with crosstalk between DNA methylation and histone modifications at distinct regulatory regions, contribute to mammalian development and tissue maturation. We also discuss how recent technological advances in single-cell and long-read methylome sequencing, along with targeted epigenome-editing, are enabling unprecedented high-resolution and mechanistic dissection of DNA methylome dynamics.

Promoter-Adjacent DNA Hypermethylation Can Downmodulate Gene Expression: TBX15 in the Muscle Lineage

TBX15, which encodes a differentiation-related transcription factor, displays promoter-adjacent DNA hypermethylation in myoblasts and skeletal muscle (psoas) that is absent from non-expressing cells in other lineages. By whole-genome bisulfite sequencing (WGBS) and enzymatic methyl-seq (EM-seq), these hypermethylated regions were found to border both sides of a constitutively unmethylated promoter. To understand the functionality of this DNA hypermethylation, we cloned the differentially methylated sequences (DMRs) in CpG-free reporter vectors and tested them for promoter or enhancer activity upon transient transfection. These cloned regions exhibited strong promoter activity and, when placed upstream of a weak promoter, strong enhancer activity specifically in myoblast host cells. In vitro CpG methylation targeted to the DMR sequences in the plasmids resulted in 86−100% loss of promoter or enhancer activity, depending on the insert sequence. These results as well as chromatin epigenetic and transcription profiles for this gene in various cell types support the hypothesis that DNA hypermethylation immediately upstream and downstream of the unmethylated promoter region suppresses enhancer/extended promoter activity, thereby downmodulating, but not silencing, expression in myoblasts and certain kinds of skeletal muscle. This promoter-border hypermethylation was not found in cell types with a silent TBX15 gene, and these cells, instead, exhibit repressive chromatin in and around the promoter. TBX18, TBX2, TBX3 and TBX1 display TBX15-like hypermethylated DMRs at their promoter borders and preferential expression in myoblasts. Therefore, promoter-adjacent DNA hypermethylation for downmodulating transcription to prevent overexpression may be used more frequently for transcription regulation than currently appreciated.

Active DNA demethylation of developmental cis-regulatory regions predates vertebrate origins

DNA methylation [5-methylcytosine (5mC)] is a repressive gene-regulatory mark required for vertebrate embryogenesis. Genomic 5mC is tightly regulated through the action of DNA methyltransferases, which deposit 5mC, and ten-eleven translocation (TET) enzymes, which participate in its active removal through the formation of 5-hydroxymethylcytosine (5hmC). TET enzymes are essential for mammalian gastrulation and activation of vertebrate developmental enhancers; however, to date, a clear picture of 5hmC function, abundance, and genomic distribution in nonvertebrate lineages is lacking. By using base-resolution 5mC and 5hmC quantification during sea urchin and lancelet embryogenesis, we shed light on the roles of nonvertebrate 5hmC and TET enzymes. We find that these invertebrate deuterostomes use TET enzymes for targeted demethylation of regulatory regions associated with developmental genes and show that the complement of identified 5hmC-regulated genes is conserved to vertebrates. This work demonstrates that active 5mC removal from regulatory regions is a common feature of deuterostome embryogenesis suggestive of an unexpected deep conservation of a major gene-regulatory module.

MC profiling: a novel approach to analyze DNA methylation heterogeneity in genome-wide bisulfite sequencing data

DNA methylation is an epigenetic mark implicated in crucial biological processes. Most of the knowledge about DNA methylation is based on bulk experiments, in which DNA methylation of genomic regions is reported as average methylation. However, average methylation does not inform on how methylated cytosines are distributed in each single DNA molecule. Here, we propose Methylation Class (MC) profiling as a genome-wide approach to the study of DNA methylation heterogeneity from bulk bisulfite sequencing experiments. The proposed approach is built on the concept of MCs, groups of DNA molecules sharing the same number of methylated cytosines. The relative abundances of MCs from sequencing reads incorporates the information on the average methylation, and directly informs on the methylation level of each molecule. By applying our approach to publicly available bisulfite-sequencing datasets, we individuated cell-to-cell differences as the prevalent contributor to methylation heterogeneity. Moreover, we individuated signatures of loci undergoing imprinting and X-inactivation, and highlighted differences between the two processes. When applying MC profiling to compare different conditions, we identified methylation changes occurring in regions with almost constant average methylation. Altogether, our results indicate that MC profiling can provide useful insights on the epigenetic status and its evolution at multiple genomic regions.

DNA methyltransferases 3A and 3B target specific sequences during mouse gastrulation

In mammalian embryos, DNA methylation is initialized to maximum levels in the epiblast by the de novo DNA methyltransferases DNMT3A and DNMT3B before gastrulation diversifies it across regulatory regions. Here we show that DNMT3A and DNMT3B are differentially regulated during endoderm and mesoderm bifurcation and study the implications in vivo and in meso-endoderm embryoid bodies. Loss of both Dnmt3a and Dnmt3b impairs exit from the epiblast state. More subtly, independent loss of Dnmt3a or Dnmt3b leads to small biases in mesoderm-endoderm bifurcation and transcriptional deregulation. Epigenetically, DNMT3A and DNMT3B drive distinct methylation kinetics in the epiblast, as can be predicted from their strand-specific sequence preferences. The enzymes compensate for each other in the epiblast, but can later facilitate lineage-specific methylation kinetics as their expression diverges. Single-cell analysis shows that differential activity of DNMT3A and DNMT3B combines with replication-linked methylation turnover to increase epigenetic plasticity in gastrulation. Together, these findings outline a dynamic model for the use of DNMT3A and DNMT3B sequence specificity during gastrulation.

Evidence that direct inhibition of transcription factor binding is the prevailing mode of gene and repeat repression by DNA methylation

Cytosine methylation efficiently silences CpG-rich regulatory regions of genes and repeats in mammalian genomes. To what extent this entails direct inhibition of transcription factor (TF) binding versus indirect inhibition via recruitment of methyl-CpG-binding domain (MBD) proteins is unclear. Here we show that combinatorial genetic deletions of all four proteins with functional MBDs in mouse embryonic stem cells, derived neurons or a human cell line do not reactivate genes or repeats with methylated promoters. These do, however, become activated by methylation-restricted TFs if DNA methylation is removed. We identify several causal TFs in neurons, including ONECUT1, which is methylation sensitive only at a motif variant. Rampantly upregulated retrotransposons in methylation-free neurons feature a CRE motif, which activates them in the absence of DNA methylation via methylation-sensitive binding of CREB1. Our study reveals methylation-sensitive TFs in vivo and argues that direct inhibition, rather than indirect repression by the tested MBD proteins, is the prevailing mechanism of methylation-mediated repression at regulatory regions and repeats.

NNMT-DNMT1 Axis is Essential for Maintaining Cancer Cell Sensitivity to Oxidative Phosphorylation Inhibition

Lacking a clear understanding of the molecular mechanism determining cancer cell sensitivity to oxidative phosphorylation (OXPHOS) inhibition limits the development of OXPHOS-targeting cancer treatment. Here, cancer cell lines sensitive or resistant to OXPHOS inhibition are identified by screening. OXPHOS inhibition-sensitive cancer cells possess increased OXPHOS activity and silenced nicotinamide N-methyltransferase (NNMT) expression. NNMT expression negatively correlates with OXPHOS inhibition sensitivity and functionally downregulates the intracellular levels of S-adenosyl methionine (SAM). Expression of DNA methyltransferase 1 (DNMT1), a SAM consumer, positively correlates with OXPHOS inhibition sensitivity. NNMT overexpression and DNMT1 inhibition render OXPHOS inhibition-sensitive cancer cells resistant. Importantly, treatments of OXPHOS inhibitors (Gboxin and Berberine) hamper the growth of mouse tumor xenografts by OXPHOS inhibition sensitive but not resistant cancer cells. What's more, the retrospective study of 62 tumor samples from a clinical trial demonstrates that administration of Berberine reduces the tumor recurrence rate of NNMTlow /DNMT1high but not NNMThigh /DNMT1low colorectal adenomas (CRAs). These results thus reveal a critical role of the NNMT-DNMT1 axis in determining cancer cell reliance on mitochondrial OXPHOS and suggest that NNMT and DNMT1 are faithful biomarkers for OXPHOS-targeting cancer therapies.

Single-cell multi-omics profiling links dynamic DNA methylation to cell fate decisions during mouse early organogenesis

BACKGROUND

Perturbation of DNA methyltransferases (DNMTs) and of the active DNA demethylation pathway via ten-eleven translocation (TET) methylcytosine dioxygenases results in severe developmental defects and embryonic lethality. Dynamic control of DNA methylation is therefore vital for embryogenesis, yet the underlying mechanisms remain poorly understood.

RESULTS

Here we report a single-cell transcriptomic atlas from Dnmt and Tet mutant mouse embryos during early organogenesis. We show that both the maintenance and de novo methyltransferase enzymes are dispensable for the formation of all major cell types at E8.5. However, DNA methyltransferases are required for silencing of prior or alternative cell fates such as pluripotency and extraembryonic programmes. Deletion of all three TET enzymes produces substantial lineage biases, in particular, a failure to generate primitive erythrocytes. Single-cell multi-omics profiling moreover reveals that this is linked to a failure to demethylate distal regulatory elements in Tet triple-knockout embryos.

CONCLUSIONS

This study provides a detailed analysis of the effects of perturbing DNA methylation on mouse organogenesis at a whole organism scale and affords new insights into the regulatory mechanisms of cell fate decisions.

Early Expression of Tet1 and Tet2 in Mouse Zygotes Altered DNA Methylation Status and Affected Embryonic Development

Ten-eleven translocation (Tet) dioxygenases can induce DNA demethylation by catalyzing 5-methylcytosine(5mC) to 5-hydroxymethylcytosine(5hmC), and play important roles during mammalian development. In mouse, Tet1 and Tet2 are not expressed in pronucleus-staged embryos and are not involved in the genomic demethylation of early zygotes. Here, we investigated the influence of Tet1 and Tet2 on methylation of parental genomes by ectopically expressing Tet1 and Tet2 in zygotes. Immunofluorescence staining showed a marked 5hmC increase in the maternal pronucleus after injection of Tet1 or Tet2 mRNA into zygotes. Whole-genome bisulfite sequencing further revealed that Tet2 greatly enhanced the global demethylation of both parental genomes, while Tet1 only promoted the paternal demethylation. Tet1 and Tet2 overexpression altered the DNA methylation across genomes, including various genic elements and germline-specific differently methylated regions. Tet2 exhibited overall stronger demethylation activity than Tet1. Either Tet1 or Tet2 overexpression impaired preimplantation embryonic development. These results demonstrated that early expression of Tet1 and Tet2 could substantially alter the zygotic methylation landscape and damage embryonic development. These findings provide new insights into understanding the function of Tet dioxygenases and the mechanism of DNA methylation in relation to embryogenesis.

Genome-wide detection of imprinted differentially methylated regions using nanopore sequencing

Imprinting is a critical part of normal embryonic development in mammals, controlled by defined parent-of-origin (PofO) differentially methylated regions (DMRs) known as imprinting control regions. Direct nanopore sequencing of DNA provides a means to detect allelic methylation and to overcome the drawbacks of methylation array and short-read technologies. Here, we used publicly available nanopore sequencing data for 12 standard B-lymphocyte cell lines to acquire the genome-wide mapping of imprinted intervals in humans. Using the sequencing data, we were able to phase 95% of the human methylome and detect 94% of the previously well-characterized, imprinted DMRs. In addition, we found 42 novel imprinted DMRs (16 germline and 26 somatic), which were confirmed using whole-genome bisulfite sequencing (WGBS) data. Analysis of WGBS data in mouse (Mus musculus), rhesus monkey (Macaca mulatta), and chimpanzee (Pan troglodytes) suggested that 17 of these imprinted DMRs are conserved. Some of the novel imprinted intervals are within or close to imprinted genes without a known DMR. We also detected subtle parental methylation bias, spanning several kilobases at seven known imprinted clusters. At these blocks, hypermethylation occurs at the gene body of expressed allele(s) with mutually exclusive H3K36me3 and H3K27me3 allelic histone marks. These results expand upon our current knowledge of imprinting and the potential of nanopore sequencing to identify imprinting regions using only parent-offspring trios, as opposed to the large multi-generational pedigrees that have previously been required.

A Critical Role for p53 during the HPV16 Life Cycle

Human papillomaviruses (HPV) are causative agents in ano-genital and oral cancers; HPV16 is the most prevalent type detected in human cancers. The HPV16 E6 protein targets p53 for proteasomal degradation to facilitate proliferation of the HPV16 infected cell. However, in HPV16 immortalized cells E6 is predominantly spliced (E6*) and unable to degrade p53. Here, we demonstrate that human foreskin keratinocytes immortalized by HPV16 (HFK+HPV16), and HPV16 positive oropharyngeal cancers, retain significant expression of p53. In addition, p53 levels increase in HPV16+ head and neck cancer cell lines following treatment with cisplatin. Introduction of full-length E6 into HFK+HPV16 resulted in attenuation of cellular growth (in hTERT immortalized HFK, E6 expression promoted enhanced proliferation). An understudied interaction is that between E2 and p53 and we investigated whether this was important for the viral life cycle. We generated mutant genomes with E2 unable to interact with p53 resulting in profound phenotypes in primary HFK. The mutant induced hyper-proliferation, but an ultimate arrest of cell growth; β-galactosidase staining demonstrated increased senescence, and COMET assays showed increased DNA damage compared with HFK+HPV16 wild-type cells. There was failure of the viral life cycle in organotypic rafts with the mutant HFK resulting in premature differentiation and reduced proliferation. The results demonstrate that p53 expression is critical during the HPV16 life cycle, and that this may be due to a functional interaction between E2 and p53. Disruption of this interaction has antiviral potential. IMPORTANCE Human papillomaviruses are causative agents in around 5% of all cancers. There are currently no antivirals available to combat these infections and cancers, therefore it remains a priority to enhance our understanding of the HPV life cycle. Here, we demonstrate that an interaction between the viral replication/transcription/segregation factor E2 and the tumor suppressor p53 is critical for the HPV16 life cycle. HPV16 immortalized cells retain significant expression of p53, and the critical role for the E2-p53 interaction demonstrates why this is the case. If the E2-p53 interaction is disrupted then HPV16 immortalized cells fail to proliferate, have enhanced DNA damage and senescence, and there is premature differentiation during the viral life cycle. Results suggest that targeting the E2-p53 interaction would have therapeutic benefits, potentially attenuating the spread of HPV16.

Comprehensive analysis to identify noncoding RNAs mediated upregulation of maternal embryonic leucine zipper kinase (MELK) correlated with poor prognosis in hepatocellular

Object: Maternal embryonic leucine zipper kinase (MELK) is involved in the development and progression of various cancers. This work investigated the usefulness of MELK in the prediction of hepatocellular carcinoma (HCC) prognosis.

Methods: Information on MELK expression was obtained by pan-cancer analysis using The Cancer Genome Atlas (TCGA) database. The TCGA-liver hepatic cancer (TCGA-LIHC), Oncomine datasets, International Cancer Genome Consortium (ICGC) datasets were used to investigate MELK expression in HCC. The prognostic roles of MELK in HCC were assessed by univariate and multivariate survival analyses. The underlying mechanism for noncoding RNAs (ncRNAs) involved in MELK expression was investigated by in silico studies, correlation, methylation, and survival analyses. The relationships between MELK expression and immune cells, immune markers, and checkpoint markers were also analyzed.

Results: (1) MELK was identified as an independent predictor of overall survival (OS) in HCC patients (MELK high vs. low expression, HR 2.469; 95% CI 1.217–5.008; p = 0.012) in a multivariate Cox analysis, with a concordance index (C-index) value of 0.727 (95% CI 0.750–0.704). (2) The noncoding RNA miR3142HG and the LINC00265/has-miR-101-3p axis were found to regulate MELK expression in HCC tissue. (3) MELK levels were linked to various immune functions, including tumor infiltration and the expression of immune checkpoints and biomarkers in HCC.

Conclusion: MELK may have an oncogenic function in HCC and was found to be up-regulated by ncRNAs and associated with immune cell infiltration and unfavorable prognosis.

CDK2 Inhibition Enhances Antitumor Immunity by Increasing IFN Response to Endogenous Retroviruses

Inhibitors of cyclin-dependent kinase-2 (CDK2) are commonly used against several solid tumors, and their primary mechanisms of action were thought to include cell proliferation arrest, induction of cancer cell apoptosis and induction of differentiation. Here, we found that CDK2 inhibition by either small molecular inhibitors or genetic Cdk2 deficiency promoted antitumor immunity in murine models of fibrosarcoma and lung carcinoma. Mechanistically, CDK2 inhibition reduced phosphorylation of RB protein and transcription of E2F-mediated DNA methyltransferase 1 (DNMT1), which resulted in increased expression of endogenous retroviral RNA and type I IFN (IFN-I) response. The increased IFN-I response subsequently promoted antitumor immunity by enhancing tumor antigen presentation and CD8+ T-cell infiltration. Our studies provide evidence that inhibition of CDK2 in cancer cells suppresses tumor growth by enhancing antitumor immune responses in the tumor microenvironment, suggesting a new mechanism to enhance antitumor immunity by CDK2 inhibitors.

A comparative methylome analysis reveals conservation and divergence of DNA methylation patterns and functions in vertebrates

Background

Cytosine DNA methylation is a heritable epigenetic mark present in most eukaryotic groups. While the patterns and functions of DNA methylation have been extensively studied in mouse and human, their conservation in other vertebrates remains poorly explored. In this study, we interrogated the distribution and function of DNA methylation in primary fibroblasts of seven vertebrate species including bio-medical models and livestock species (human, mouse, rabbit, dog, cow, pig, and chicken).

Results

Our data highlight both divergence and conservation of DNA methylation patterns and functions. We show that the chicken genome is hypomethylated compared to other vertebrates. Furthermore, compared to mouse, other species show a higher frequency of methylation of CpG-rich DNA. We reveal the conservation of large unmethylated valleys and patterns of DNA methylation associated with X-chromosome inactivation through vertebrate evolution and make predictions of conserved sets of imprinted genes across mammals. Finally, using chemical inhibition of DNA methylation, we show that the silencing of germline genes and endogenous retroviruses (ERVs) are conserved functions of DNA methylation in vertebrates.

Conclusions

Our data highlight conserved properties of DNA methylation in vertebrate genomes but at the same time point to differences between mouse and other vertebrate species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12915-022-01270-x.

Common mechanisms cannot explain time- and dose-dependent DNA methylation changes in earthworms exposed to cadmium

DNA hypermethylation caused by environmental pollutants like cadmium (Cd) has already been demonstrated in many invertebrates, including earthworms. However, the exact epigenetic mechanisms that drive this hypermethylation are largely unknown and even basic DNA methylation and demethylation processes are hardly characterized. Therefore, we used an important bioindicator, the earthworm Lumbricus terrestris, as a model organism to determine time- and dose-dependent effects of Cd on global and gene-specific DNA methylation and its underlying mechanisms. We revealed Cd-induced adenine and cytosine hypermethylation using specific antibodies in dot blots and found that the methylation level of adenine compared to cytosine changed even to a bigger extent. However, the levels of hydroxymethylated cytosine did not differ between treatment groups. General methylation and demethylation components like methyltransferases (DNMT1 and 3), and ten-eleven translocation (TET) genes were confirmed in L. terrestris by quantitative RealTime PCR. However, neither gene expression, nor DNMT and TET enzyme activity showed significant differences in the Cd exposure groups. Using bisulfite conversion and sequencing, gene body methylation (gbm) of metallothionein 2 (MT2), one of the most important detoxification proteins, was characterized. Cd-dependent changes in MT2 gbm could, however, not be correlated to MT2 gene activity evaluated by quantitative RealTime PCR. Future directions as well as missing links are discussed in the present study hinting towards the importance of studying epigenetic marks and mechanistic insights in a broad variety of species to deepen our knowledge on the effects of changing environmental conditions.

Interplay between chromatin marks in development and disease

DNA methylation (DNAme) and histone post-translational modifications (PTMs) have important roles in transcriptional regulation. Although many reports have characterized the functions of such chromatin marks in isolation, recent genome-wide studies reveal surprisingly complex interactions between them. Here, we focus on the interplay between DNAme and methylation of specific lysine residues on the histone H3 tail. We describe the impact of genetic perturbation of the relevant methyltransferases in the mouse on the landscape of chromatin marks as well as the transcriptome. In addition, we discuss the specific neurodevelopmental growth syndromes and cancers resulting from pathogenic mutations in the human orthologues of these genes. Integrating these observations underscores the fundamental importance of crosstalk between DNA and histone H3 methylation in development and disease.

DNA Methylation and Detoxification in the Earthworm Lumbricus terrestris Exposed to Cadmium and the DNA Demethylation Agent 5-aza-2'-deoxycytidine

Earthworms are well-established model organisms for testing the effects of heavy metal pollution. How DNA methylation affects cadmium (Cd) detoxification processes such as the expression of metallothionein 2 (MT2), however, is largely unknown. We therefore exposed Lumbricus terrestris to 200 mg concentrations of Cd and 5-aza-2'-deoxycytidine (Aza), a demethylating agent, and sampled tissue and coelomocytes, cells of the innate immune system, for 48 h. MT2 transcription significantly increased in the Cd- and Cd-Aza-treated groups. In tissue samples, a significant decrease in MT2 in the Aza-treated group was detected, showing that Aza treatment inhibits basal MT2 gene activity but has no effect on Cd-induced MT2 levels. Although Cd repressed the gene expression of DNA-(cytosine-5)-methyltransferase-1 (DNMT1), which is responsible for maintaining DNA methylation, DNMT activity was unchanged, meaning that methylation maintenance was not affected in coelomocytes. The treatment did not influence DNMT3, which mediates de novo methylation, TET gene expression, which orchestrates demethylation, and global levels of hydroxymethylcytosine (5hmC), a product of the demethylation process. Taken together, this study indicates that Aza inhibits basal gene activity, in contrast to Cd-induced MT2 gene expression, but does not affect global DNA methylation. We therefore conclude that Cd detoxification based on the induction of MT2 does not relate to DNA methylation changes.

Genetic variation at mouse and human ribosomal DNA influences associated epigenetic states

BACKGROUND

Ribosomal DNA (rDNA) displays substantial inter-individual genetic variation in human and mouse. A systematic analysis of how this variation impacts epigenetic states and expression of the rDNA has thus far not been performed.

RESULTS

Using a combination of long- and short-read sequencing, we establish that 45S rDNA units in the C57BL/6J mouse strain exist as distinct genetic haplotypes that influence the epigenetic state and transcriptional output of any given unit. DNA methylation dynamics at these haplotypes are dichotomous and life-stage specific: at one haplotype, the DNA methylation state is sensitive to the in utero environment, but refractory to post-weaning influences, whereas other haplotypes entropically gain DNA methylation during aging only. On the other hand, individual rDNA units in human show limited evidence of genetic haplotypes, and hence little discernible correlation between genetic and epigenetic states. However, in both species, adjacent units show similar epigenetic profiles, and the overall epigenetic state at rDNA is strongly positively correlated with the total rDNA copy number. Analysis of different mouse inbred strains reveals that in some strains, such as 129S1/SvImJ, the rDNA copy number is only approximately 150 copies per diploid genome and DNA methylation levels are < 5%.

CONCLUSIONS

Our work demonstrates that rDNA-associated genetic variation has a considerable influence on rDNA epigenetic state and consequently rRNA expression outcomes. In the future, it will be important to consider the impact of inter-individual rDNA (epi)genetic variation on mammalian phenotypes and diseases.

Intragenic CpG Islands and Their Impact on Gene Regulation

The mammalian genome is depleted in CG dinucleotides, except at protected regions where they cluster as CpG islands (CGIs). CGIs are gene regulatory hubs and serve as transcription initiation sites and are as expected, associated with gene promoters. Advances in genomic annotations demonstrate that a quarter of CGIs are found within genes. Such intragenic regions are repressive environments, so it is surprising that CGIs reside here and even more surprising that some resist repression and are transcriptionally active within a gene. Hence, intragenic CGI positioning within genes is not arbitrary and is instead, selected for. As a wealth of recent studies demonstrate, intragenic CGIs are embedded within genes and consequently, influence ‘host’ gene mRNA isoform length and expand transcriptome diversity.

Analysis of epigenetic features characteristic of L1 loci expressed in human cells

Only a select few L1 loci in the human genome are expressed in any given cell line or organ, likely to minimize damage done to the genome. The epigenetic features and requirements of expressed L1 loci are currently unknown. Using human cells and comprehensive epigenetic analysis of individual expressed and unexpressed L1 loci, we determined that endogenous L1 transcription depends on a combination of epigenetic factors, including open chromatin, activating histone modifications, and hypomethylation at the L1 promoter. We demonstrate that the L1 promoter seems to require interaction with enhancer elements for optimal function. We utilize epigenetic context to predict the expression status of L1Hs loci that are poorly mappable with RNA-Seq. Our analysis identified a population of ‘transitional’ L1 loci that likely have greater potential to be activated during the epigenetic dysregulation seen in tumors and during aging because they are the most responsive to targeted CRISPR-mediated delivery of trans-activating domains. We demonstrate that an engineered increase in endogenous L1 mRNA expression increases Alu mobilization. Overall, our findings present the first global and comprehensive analysis of epigenetic status of individual L1 loci based on their expression status and demonstrate the importance of epigenetic context for L1 expression heterogeneity.

Evolution of DNA Methylation Across Ecdysozoa

DNA methylation is a crucial, abundant mechanism of gene regulation in vertebrates. It is less prevalent in many other metazoan organisms and completely absent in some key model species, such as Drosophila melanogaster and Caenorhabditis elegans. We report here a comprehensive study of the presence and absence of DNA methyltransferases (DNMTs) in 138 Ecdysozoa, covering Arthropoda, Nematoda, Priapulida, Onychophora, and Tardigrada. Three of these phyla have not been investigated for the presence of DNA methylation before. We observe that the loss of individual DNMTs independently occurred multiple times across ecdysozoan phyla. We computationally predict the presence of DNA methylation based on CpG rates in coding sequences using an implementation of Gaussian Mixture Modeling, MethMod. Integrating both analysis we predict two previously unknown losses of DNA methylation in Ecdysozoa, one within Chelicerata (Mesostigmata) and one in Tardigrada. In the early-branching Ecdysozoa Priapulus caudatus, we predict the presence of a full set of DNMTs and the presence of DNA methylation. We are therefore showing a very diverse and independent evolution of DNA methylation in different ecdysozoan phyla spanning a phylogenetic range of more than 700 million years.

Exogenous melatonin protects preimplantation embryo development from decabromodiphenyl ethane-induced circadian rhythm disorder and endogenous melatonin reduction

Decabromodiphenyl ethane (DBDPE) is a novel flame retardant that is widely used in plastics, electronic products, building materials and textiles. Our previous studies have revealed the oocyte toxicity of DBDPE, but the effect of DBDPE on preimplantation embryo development has not been reported. Here, we investigated whether and how DBDPE exposure affects preimplantation embryo development. Adult female mice were orally exposed to DBDPE (0, 5, 50, 500 μg/kg bw/day) for 14 days. First, we found that after DBDPE exposure, mice showed obvious circadian rhythm disorder. Moreover, the development of preimplantation embryos was inhibited in DBDPE-exposed mice after pregnancy. Then, we further explored and revealed that DBDPE exposure reduced the endogenous melatonin (MLT) level during pregnancy, thereby inhibiting the development of preimplantation embryos. Furthermore, we discovered that exogenous MLT supplementation (15 mg/kg bw/day) rescued the inhibition of preimplantation embryo development induced by DBDPE, and a mechanistic study demonstrated that exogenous MLT inhibited the overexpression of ROS and DNA methylation at the 5-position of cytosine (5-mC) in DBDPE-exposed preimplantation embryos. Simultaneously, MLT ameliorated the DBDPE-induced mitochondrial dysfunction by increasing the mitochondrial membrane potential (MMP), ATP, and Trp1 expression. Additionally, MLT restored DBDPE-induced changes in zona pellucida (ZP) hardness and trophectoderm (TE) cortical tension. Finally, the protective effect of MLT on embryos ameliorated the adverse reproductive outcomes (dead fetus, fetus with abnormal liver, fetal weight loss) induced by DBDPE. Collectively, DBDPE induced preimplantation embryo damage leading to adverse reproductive outcomes, and MLT has emerged as a potential tool to rescue adverse reproductive outcomes induced by DBDPE.

The DNMT1 inhibitor GSK-3484862 mediates global demethylation in murine embryonic stem cells

Background

DNA methylation plays an important role in regulating gene expression in mammals. The covalent DNMT1 inhibitors 5-azacytidine and decitabine are widely used in research to reduce DNA methylation levels, but they impart severe cytotoxicity which limits their demethylation capability and confounds interpretation of experiments. Recently, a non-covalent inhibitor of DNMT1 called GSK-3484862 was developed by GlaxoSmithKline. We sought to determine whether GSK-3484862 can induce demethylation more effectively than 5-azanucleosides. Murine embryonic stem cells (mESCs) are an ideal cell type in which to conduct such experiments, as they have a high degree of DNA methylation but tolerate dramatic methylation loss.

Results

We determined the cytotoxicity and optimal concentration of GSK-3484862 by treating wild-type (WT) or Dnmt1/3a/3b triple knockout (TKO) mESC with different concentrations of the compound, which was obtained from two commercial sources. Concentrations of 10 µM or below were readily tolerated for 14 days of culture. Known DNA methylation targets such as germline genes and GLN-family transposons were upregulated within 2 days of the start of GSK-3484862 treatment. By contrast, 5-azacytidine and decitabine induced weaker upregulation of methylated genes and extensive cell death. Whole-genome bisulfite sequencing showed that treatment with GSK-3484862 induced dramatic DNA methylation loss, with global CpG methylation levels falling from near 70% in WT mESC to less than 18% after 6 days of treatment with GSK-3484862. The treated cells showed a methylation level and pattern similar to that observed in Dnmt1-deficient mESCs.

Conclusions

GSK-3484862 mediates striking demethylation in mESCs with minimal non-specific toxicity.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13072-021-00429-0.

Differential DNA methylation across environments has no effect on gene expression in the eastern oyster

It has been hypothesized that environmentally induced changes to gene body methylation could facilitate adaptive transgenerational responses to changing environments. We compared patterns of global gene expression (Tag-seq) and gene body methylation (reduced representation bisulfite sequencing) in 80 eastern oysters Crassostrea virginica from six full-sib families, common gardened for 14 months at two sites in the northern Gulf of Mexico that differed in mean salinity. At the time of sampling, oysters from the two sites differed in mass by 60% and in parasite loads by nearly two orders of magnitude. They also differentially expressed 35% of measured transcripts. However, we observed differential methylation at only 1.4% of potentially methylated loci in comparisons between individuals from these different environments, and little correspondence between differential methylation and differential gene expression. Instead, methylation patterns were largely driven by genetic differences among families, with a PERMANOVA analysis indicating nearly a two orders of magnitude greater number of genes differentially methylated between families than between environments. An analysis of CpG observed/expected values (CpG O/E) across the C. virginica genome showed a distinct bimodal distribution, with genes from the first cluster showing the lower CpG O/E values, greater methylation and higher and more stable gene expression, while genes from the second cluster showed lower methylation, and lower and more variable gene expression. Taken together, the differential methylation results suggest that only a small portion of the C. virginica genome is affected by environmentally induced changes in methylation. At this point, there is little evidence to suggest that environmentally induced methylation states would play a leading role in regulating gene expression responses to new environments.

Repression of germline genes by PRC1.6 and SETDB1 in the early embryo precedes DNA methylation-mediated silencing

Silencing of a subset of germline genes is dependent upon DNA methylation (DNAme) post-implantation. However, these genes are generally hypomethylated in the blastocyst, implicating alternative repressive pathways before implantation. Indeed, in embryonic stem cells (ESCs), an overlapping set of genes, including germline "genome-defence" (GGD) genes, are upregulated following deletion of the H3K9 methyltransferase SETDB1 or subunits of the non-canonical PRC1 complex PRC1.6. Here, we show that in pre-implantation embryos and naïve ESCs (nESCs), hypomethylated promoters of germline genes bound by the PRC1.6 DNA-binding subunits MGA/MAX/E2F6 are enriched for RING1B-dependent H2AK119ub1 and H3K9me3. Accordingly, repression of these genes in nESCs shows a greater dependence on PRC1.6 than DNAme. In contrast, GGD genes are hypermethylated in epiblast-like cells (EpiLCs) and their silencing is dependent upon SETDB1, PRC1.6/RING1B and DNAme, with H3K9me3 and DNAme establishment dependent upon MGA binding. Thus, GGD genes are initially repressed by PRC1.6, with DNAme subsequently engaged in post-implantation embryos.

Transposable Element Dynamics and Regulation during Zygotic Genome Activation in Mammalian Embryos and Embryonic Stem Cell Model Systems

Transposable elements (TEs) are mobile genetic sequences capable of duplicating and reintegrating at new regions within the genome. A growing body of evidence has demonstrated that these elements play important roles in host genome evolution, despite being traditionally viewed as parasitic elements. To prevent ectopic activation of TE transposition and transcription, they are epigenetically silenced in most somatic tissues. Intriguingly, a specific class of TEs-retrotransposons-is transiently expressed at discrete phases during mammalian development and has been linked to the establishment of totipotency during zygotic genome activation (ZGA). While mechanisms controlling TE regulation in somatic tissues have been extensively studied, the significance underlying the unique transcriptional reactivation of retrotransposons during ZGA is only beginning to be uncovered. In this review, we summarize the expression dynamics of key retrotransposons during ZGA, focusing on findings from in vivo totipotent embryos and in vitro totipotent-like embryonic stem cells (ESCs). We then dissect the functions of retrotransposons and discuss how their transcriptional activities are finetuned during early stages of mammalian development.

SETD2-mediated epigenetic regulation of noncanonical Wnt5A during osteoclastogenesis

To define the role of SETD2 in the WNT5a signaling in the context of osteoclastogenesis, we exploited two different models: in vitro osteoclast differentiation, and K/BxN serum-induced arthritis model. We found that SETD2 and WNT5a were upregulated during osteoclast differentiation and after induction of arthritis. Using gain- and loss-of-function approaches in the myeloid cell, we confirmed that SETD2 regulated the osteoclast markers, and WNT5a via modulating active histone marks by enriching H3K36me3, and by reducing repressive H3K27me3 mark. Additionally, during osteoclastic differentiation, the transcription of Wnt5a was also associated with the active histone H3K9 and H4K8 acetylations. Mechanistically, SETD2 directed induction of NF-κβ expression facilitated the recruitment of H3K9Ac and H4K8Ac around the TSS region of the Wnt5a gene, thereby, assisting osteoclast differentiation. Together these findings for the first time revealed that SETD2 mediated epigenetic regulation of Wnt5a plays a critical role in osteoclastogenesis and induced arthritis. Model for the Role of SETD2 dependent regulation of osteoclastic differentiation. A In monocyte cells SETD2-dependent H3K36 trimethylation help to create open chromatin region along with active enhancer mark, H3K27Ac. This chromatin state facilitated the loss of a suppressive H3K27me3 mark. B Additionally, SETD2 mediated induction of NF-κβ expression leads to the recruitment of histone acetyl transferases, P300/PCAF, to the Wnt5a gene and establish H3K9Ac and H4K8Ac marks. Along with other activation marks, these acetylation marks help in Wnt5a transcription which leads to osteoclastogenesis.

Disabling de novo DNA methylation in embryonic stem cells allows an illegitimate fate trajectory

Genome remethylation is essential for mammalian development but specific reasons are unclear. Here we examined embryonic stem (ES) cell fate in the absence of de novo DNA methyltransferases. We observed that ES cells deficient for both Dnmt3a and Dnmt3b are rapidly eliminated from chimeras. On further investigation we found that in vivo and in vitro the formative pluripotency transition is derailed toward production of trophoblast. This aberrant trajectory is associated with failure to suppress activation of Ascl2Ascl2 encodes a bHLH transcription factor expressed in the placenta. Misexpression of Ascl2 in ES cells provokes transdifferentiation to trophoblast-like cells. Conversely, Ascl2 deletion rescues formative transition of Dnmt3a/b mutants and improves contribution to chimeric epiblast. Thus, de novo DNA methylation safeguards against ectopic activation of Ascl2 However, Dnmt3a/b-deficient cells remain defective in ongoing embryogenesis. We surmise that multiple developmental transitions may be secured by DNA methylation silencing potentially disruptive genes.

DNA Methylation Dynamics in the Female Germline and Maternal-Effect Mutations That Disrupt Genomic Imprinting

Genomic imprinting is an epigenetic marking process that results in the monoallelic expression of a subset of genes. Many of these ‘imprinted’ genes in mice and humans are involved in embryonic and extraembryonic growth and development, and some have life-long impacts on metabolism. During mammalian development, the genome undergoes waves of (re)programming of DNA methylation and other epigenetic marks. Disturbances in these events can cause imprinting disorders and compromise development. Multi-locus imprinting disturbance (MLID) is a condition by which imprinting defects touch more than one locus. Although most cases with MLID present with clinical features characteristic of one imprinting disorder. Imprinting defects also occur in ‘molar’ pregnancies-which are characterized by highly compromised embryonic development-and in other forms of reproductive compromise presenting clinically as infertility or early pregnancy loss. Pathogenic variants in some of the genes encoding proteins of the subcortical maternal complex (SCMC), a multi-protein complex in the mammalian oocyte, are responsible for a rare subgroup of moles, biparental complete hydatidiform mole (BiCHM), and other adverse reproductive outcomes which have been associated with altered imprinting status of the oocyte, embryo and/or placenta. The finding that defects in a cytoplasmic protein complex could have severe impacts on genomic methylation at critical times in gamete or early embryo development has wider implications beyond these relatively rare disorders. It signifies a potential for adverse maternal physiology, nutrition, or assisted reproduction to cause epigenetic defects at imprinted or other genes. Here, we review key milestones in DNA methylation patterning in the female germline and the embryo focusing on humans. We provide an overview of recent findings regarding DNA methylation deficits causing BiCHM, MLID, and early embryonic arrest. We also summarize identified SCMC mutations with regard to early embryonic arrest, BiCHM, and MLID.

Extracellular Matrix Stiffness Regulates DNA Methylation by PKCα-Dependent Nuclear Transport of DNMT3L

Extracellular matrix (ECM) stiffness has profound effects on the regulation of cell functions. DNA methylation is an important epigenetic modification governing gene expression. However, the effects of ECM stiffness on DNA methylation remain elusive. Here, it is reported that DNA methylation is sensitive to ECM stiffness, with a global hypermethylation under stiff ECM condition in mouse embryonic stem cells (mESCs) and embryonic fibroblasts compared with soft ECM. Stiff ECM enhances DNA methylation of both promoters and gene bodies, especially the 5' promoter regions of pluripotent genes. The enhanced DNA methylation is functionally required for the loss of pluripotent gene expression in mESCs grown on stiff ECM. Further experiments reveal that the nuclear transport of DNA methyltransferase 3-like (DNMT3L) is promoted by stiff ECM in a protein kinase C α (PKCα)-dependent manner and DNMT3L can be binding to Nanog promoter regions during cell-ECM interactions. These findings unveil DNA methylation as a novel target for the mechanical sensing mechanism of ECM stiffness, which provides a conserved mechanism for gene expression regulation during cell-ECM interactions.

BANP opens chromatin and activates CpG-island-regulated genes

The majority of gene transcripts generated by RNA polymerase II in mammalian genomes initiate at CpG island (CGI) promoters^1,2, yet our understanding of their regulation remains limited. This is in part due to the incomplete information that we have on transcription factors, their DNA-binding motifs and which genomic binding sites are functional in any given cell type^3-5. In addition, there are orphan motifs without known binders, such as the CGCG element, which is associated with highly expressed genes across human tissues and enriched near the transcription start site of a subset of CGI promoters^6-8. Here we combine single-molecule footprinting with interaction proteomics to identify BTG3-associated nuclear protein (BANP) as the transcription factor that binds this element in the mouse and human genome. We show that BANP is a strong CGI activator that controls essential metabolic genes in pluripotent stem and terminally differentiated neuronal cells. BANP binding is repelled by DNA methylation of its motif in vitro and in vivo, which epigenetically restricts most binding to CGIs and accounts for differential binding at aberrantly methylated CGI promoters in cancer cells. Upon binding to an unmethylated motif, BANP opens chromatin and phases nucleosomes. These findings establish BANP as a critical activator of a set of essential genes and suggest a model in which the activity of CGI promoters relies on methylation-sensitive transcription factors that are capable of chromatin opening.

Single cell RNA-seq reveals genes vital to in vitro fertilized embryos and parthenotes in pigs

Successful early embryo development requires the correct reprogramming and configuration of gene networks by the timely and faithful execution of zygotic genome activation (ZGA). However, the regulatory principle of molecular elements and circuits fundamental to embryo development remains largely obscure. Here, we profiled the transcriptomes of single zygotes and blastomeres, obtained from in vitro fertilized (IVF) or parthenogenetically activated (PA) porcine early embryos (1- to 8-cell), focusing on the gene expression dynamics and regulatory networks associated with maternal-to-zygote transition (MZT) (mainly maternal RNA clearance and ZGA). We found that minor and major ZGAs occur at 1-cell and 4-cell stages for both IVF and PA embryos, respectively. Maternal RNAs gradually decay from 1- to 8-cell embryos. Top abundantly expressed genes (CDV3, PCNA, CDR1, YWHAE, DNMT1, IGF2BP3, ARMC1, BTG4, UHRF2 and gametocyte-specific factor 1-like) in both IVF and PA early embryos identified are of vital roles for embryo development. Differentially expressed genes within IVF groups are different from that within PA groups, indicating bi-parental and maternal-only embryos have specific sets of mRNAs distinctly decayed and activated. Pathways enriched from DEGs showed that RNA associated pathways (RNA binding, processing, transport and degradation) could be important. Moreover, mitochondrial RNAs are found to be actively transcribed, showing dynamic expression patterns, and for DNA/H3K4 methylation and transcription factors as well. Taken together, our findings provide an important resource to investigate further the epigenetic and genome regulation of MZT events in early embryos of pigs.

Hydroxycarbamide effects on DNA methylation and gene expression in myeloproliferative neoplasms

Hydroxycarbamide (HC, hydroxyurea) is a cytoreductive drug inducing cell cycle blockade. However, emerging evidence suggests that HC plays a role in the modulation of transcription through the activity of transcription factors and DNA methylation. Examining the global mechanism of action of HC in the context of myeloproliferative neoplasms (MPNs), for which HC is the first-line treatment, will provide a better understanding of its molecular effects. To explore the effects of HC genome-wide, transcriptomic analyses were performed on two clinically relevant cell types at different stages of differentiation treated with HC in a murine MPN model. This study was replicated in MPN patients by profiling genome-wide gene expression and DNA methylation using patient blood samples collected longitudinally, before and following HC exposure. The effects of HC on the transcriptome were not only associated with cell cycle interruption but also with hematopoietic functions. Moreover, a group of genes were restored to normal expression levels in murine hematopoietic stem cells (HSCs) following drug treatment, including the master regulator of hematopoiesis, RUNX1. In humans, HC significantly modifies DNA methylation levels in HSCs at several distal regulatory regions, which we show to be associated with SPI1 binding sites and at the SPI1 locus itself. We have identified novel targets of HC that include pivotal transcription factors involved in hematopoiesis, and for the first time we report abnormal methylation patterns in MPN patients at the master regulator gene SPI1 and its distal binding sites, which HC is able to restore to normal levels.

E2F6 initiates stable epigenetic silencing of germline genes during embryonic development

In mouse development, long-term silencing by CpG island DNA methylation is specifically targeted to germline genes; however, the molecular mechanisms of this specificity remain unclear. Here, we demonstrate that the transcription factor E2F6, a member of the polycomb repressive complex 1.6 (PRC1.6), is critical to target and initiate epigenetic silencing at germline genes in early embryogenesis. Genome-wide, E2F6 binds preferentially to CpG islands in embryonic cells. E2F6 cooperates with MGA to silence a subgroup of germline genes in mouse embryonic stem cells and in embryos, a function that critically depends on the E2F6 marked box domain. Inactivation of E2f6 leads to a failure to deposit CpG island DNA methylation at these genes during implantation. Furthermore, E2F6 is required to initiate epigenetic silencing in early embryonic cells but becomes dispensable for the maintenance in differentiated cells. Our findings elucidate the mechanisms of epigenetic targeting of germline genes and provide a paradigm for how transient repression signals by DNA-binding factors in early embryonic cells are translated into long-term epigenetic silencing during mouse development.

DNA methylation targets CpG island promoters of germline genes to repress their expression in mouse somatic cells. Here the authors show that a transcription factor E2F6 is required to target CpG island DNA methylation and epigenetic silencing to germline genes during early mouse development.

Epigenetic clocks reveal a rejuvenation event during embryogenesis followed by aging

The notion that the germ line does not age goes back to the 19th-century ideas of August Weismann. However, being metabolically active, the germ line accumulates damage and other changes over time, i.e., it ages. For new life to begin in the same young state, the germ line must be rejuvenated in the offspring. Here, we developed a multi-tissue epigenetic clock and applied it, together with other aging clocks, to track changes in biological age during mouse and human prenatal development. This analysis revealed a significant decrease in biological age, i.e., rejuvenation, during early stages of embryogenesis, followed by an increase in later stages. We further found that pluripotent stem cells do not age even after extensive passaging and that the examined epigenetic age dynamics is conserved across species. Overall, this study uncovers a natural rejuvenation event during embryogenesis and suggests that the minimal biological age (ground zero) marks the beginning of organismal aging.

Loosening chromatin and dysregulated transcription: a perspective on cryptic transcription during mammalian aging

Cryptic transcription, the initiation of transcription from non-promoter regions within a gene body, is a type of transcriptional dysregulation that occurs throughout eukaryotes. In mammals, cryptic transcription is normally repressed at the level of chromatin, and this process is increased upon perturbation of complexes that increase intragenic histone H3 lysine 4 methylation or decrease intragenic H3 lysine 36 methylation, DNA methylation, or nucleosome occupancy. Significantly, similar changes to chromatin structure occur during aging, and, indeed, recent work indicates that cryptic transcription is elevated during aging in mammalian stem cells. Although increased cryptic transcription is known to promote aging in yeast, whether elevated cryptic transcription also contributes to mammalian aging is unclear. There is ample evidence that perturbations known to increase cryptic transcription are deleterious in embryonic and adult stem cells, and in some cases phenocopy certain aging phenotypes. Furthermore, an increase in cryptic transcription requires or impedes pathways that are known to have reduced function during aging, potentially exacerbating other aging phenotypes. Thus, we propose that increased cryptic transcription contributes to mammalian stem cell aging.

An effective method to resolve ambiguous bisulfite-treated reads

BACKGROUND

The combination of the bisulfite treatment and the next-generation sequencing is an important method for methylation analysis, and aligning the bisulfite-treated reads (BS-reads) is the critical step for the downstream applications. As bisulfite treatment reduces the complexity of the sequences, a large portion of BS-reads might be aligned to multiple locations of the reference genome ambiguously, called multireads. These multireads cannot be employed in the downstream applications since they are likely to introduce artifacts. To identify the best mapping location of each multiread, existing Bayesian-based methods calculate the probability of the read at each position by considering how does it overlap with unique mapped reads. However, [Formula: see text]% of multireads are not overlapped with any unique reads, which are unresolvable for existing method.

RESULTS

Here we propose a novel method (EM-MUL) that not only rescues multireads overlapped with unique reads, but also uses the overall coverage and accurate base-level alignment to resolve multireads that cannot be handled by current methods. We benchmark our method on both simulated datasets and real datasets. Experimental results show that it is able to align more than 80% of multireads to the best mapping position with very high accuracy.

CONCLUSIONS

EM-MUL is an effective method designed to accurately determine the best mapping position of multireads in BS-reads. For the downstream applications, it is useful to improve the methylation resolution on the repetitive regions of genome. EM-MUL is free available at https://github.com/lmylynn/EM-MUL.

Phenotypically concordant but epigenetically discordant monozygotic dichorionic diamniotic twins with Beckwith-Wiedemann syndrome

Beckwith-Wiedemann syndrome (BWS) is an imprinting disorder caused by (epi)genetic alterations. The incidence of monozygotic (MZ) twins in BWS is higher than in the general population. Most MZ twins with BWS are female and have phenotypical discordance: one twin is clinically diagnosed with BWS, while the other shows a mild or normal phenotype. The most frequent (epi)genetic alteration in MZ twins is loss of methylation of imprinting control region 2 (ICR2-LOM) at 11p15.5. Intriguingly, ICR2-LOM is usually found in the peripheral blood leukocytes (PBL) of both twins, even if they are clinically discordant. Here, we present a rare pair of MZ dichorionic diamniotic female twins with BWS and concordant phenotypes (a Beckwith-Wiedemann spectrum score of 5 in each twin). Molecular analysis of genomic DNA from PBL revealed ICR2-LOM in one twin but not the other. Our analyses suggest that ICR2-LOM occurred between days 1 and 3 after fertilization, followed by twinning. We speculate that during embryogenesis, ICR2-LOM cells were distributed to the hematopoietic stem cells in different ratios in the two fetuses, and also to commonly affected tissues, such as the tongue, in similar ratios, although we were unable to analyze any tissues other than PBL.

Staying true to yourself: mechanisms of DNA methylation maintenance in mammals

DNA methylation is essential to development and cellular physiology in mammals. Faulty DNA methylation is frequently observed in human diseases like cancer and neurological disorders. Molecularly, this epigenetic mark is linked to other chromatin modifications and it regulates key genomic processes, including transcription and splicing. Each round of DNA replication generates two hemi-methylated copies of the genome. These must be converted back to symmetrically methylated DNA before the next S-phase, or the mark will fade away; therefore the maintenance of DNA methylation is essential. Mechanistically, the maintenance of this epigenetic modification takes place during and after DNA replication, and occurs within the very dynamic context of chromatin re-assembly. Here, we review recent discoveries and unresolved questions regarding the mechanisms, dynamics and fidelity of DNA methylation maintenance in mammals. We also discuss how it could be regulated in normal development and misregulated in disease.

Interplay between Histone and DNA Methylation Seen through Comparative Methylomes in Rare Mendelian Disorders

DNA methylation (DNAme) profiling is used to establish specific biomarkers to improve the diagnosis of patients with inherited neurodevelopmental disorders and to guide mutation screening. In the specific case of mendelian disorders of the epigenetic machinery, it also provides the basis to infer mechanistic aspects with regard to DNAme determinants and interplay between histone and DNAme that apply to humans. Here, we present comparative methylomes from patients with mutations in the de novo DNA methyltransferases DNMT3A and DNMT3B, in their catalytic domain or their N-terminal parts involved in reading histone methylation, or in histone H3 lysine (K) methylases NSD1 or SETD2 (H3 K36) or KMT2D/MLL2 (H3 K4). We provide disease-specific DNAme signatures and document the distinct consequences of mutations in enzymes with very similar or intertwined functions, including at repeated sequences and imprinted loci. We found that KMT2D and SETD2 germline mutations have little impact on DNAme profiles. In contrast, the overlapping DNAme alterations downstream of NSD1 or DNMT3 mutations underlines functional links, more specifically between NSD1 and DNMT3B at heterochromatin regions or DNMT3A at regulatory elements. Together, these data indicate certain discrepancy with the mechanisms described in animal models or the existence of redundant or complementary functions unforeseen in humans.

Mitotic inheritance of DNA methylation: more than just copy and paste

Decades of investigation on DNA methylation have led to deeper insights into its metabolic mechanisms and biological functions. This understanding was fueled by the recent development of genome editing tools and our improved capacity for analyzing the global DNA methylome in mammalian cells. This review focuses on the maintenance of DNA methylation patterns during mitotic cell division. We discuss the latest discoveries of the mechanisms for the inheritance of DNA methylation as a stable epigenetic memory. We also highlight recent evidence showing the rapid turnover of DNA methylation as a dynamic gene regulatory mechanism. A body of work has shown that altered DNA methylomes are common features in aging and disease. We discuss the potential links between methylation maintenance mechanisms and disease-associated methylation changes.

CXCR4 Regulates Temporal Differentiation via PRC1 Complex in Organogenesis of Epithelial Glands

CXC-chemokine receptor type 4 (CXCR4), a 7-transmembrane receptor family member, displays multifaceted roles, participating in immune cell migration, angiogenesis, and even adipocyte metabolism. However, the activity of such a ubiquitously expressed receptor in epithelial gland organogenesis has not yet been fully explored. To investigate the relationship between CXCL12/CXCR4 signaling and embryonic glandular organogenesis, we used an ex vivo culture system with live imaging and RNA sequencing to elucidate the transcriptome and protein-level signatures of AMD3100, a potent abrogating reagent of the CXCR4-CXCL12 axis, imprinted on the developing organs. Immunostaining results showed that CXCR4 was highly expressed in embryonic submandibular gland, lung, and pancreas, especially at the periphery of end buds containing numerous embryonic stem/progenitor cells. Despite no significant increase in apoptosis, AMD3100-treated epithelial organs showed a retarded growth with significantly slower branching and expansion. Further analyses with submandibular glands revealed that such responses resulted from the AMD3100-induced precocious differentiation of embryonic epithelial cells, losing mitotic activity. RNA sequencing analysis revealed that inhibition of CXCR4 significantly down-regulated polycomb repressive complex (PRC) components, known as regulators of DNA methylation. Treatment with PRC inhibitor recapitulated the AMD3100-induced precocious differentiation. Our results indicate that the epigenetic modulation by the PRC-CXCR12/CXCR4 signaling axis is crucial for the spatiotemporal regulation of proliferation and differentiation of embryonic epithelial cells during embryonic glandular organogenesis.

Epigenetic activities in erythroid cell gene regulation

Interest in the role of epigenetic mechanisms in human biology has exponentially increased over the past several decades. The multitude of opposing and context-dependent chromatin-modifying enzymes/coregulator complexes is just beginning to be understood at a molecular level. This science has benefitted tremendously from studies of erythropoiesis, in which a series of β-globin genes are in sequence turned "on" and "off," serving as a fascinating model of coordinated gene expression. We, therefore, describe here epigenetic complexes about which we know most, using erythropoiesis as the context. The biochemical insights lay the foundation for proposing and developing novel treatments for diseases of red cells and of erythropoiesis, identifying for example epigenetic enzymes that can be drugged to manipulate β-globin locus regulation, to favor activation of unmutated fetal hemoglobin over mutated adult β-globin genes to treat sickle cell disease and β-thalassemias. Other potential translational applications are in redirecting hematopoietic commitment decisions, as treatment for bone marrow failure syndromes.

A DNA methylation state transition model reveals the programmed epigenetic heterogeneity in human pre-implantation embryos

BACKGROUND

During mammalian early embryogenesis, expression and epigenetic heterogeneity emerge before the first cell fate determination, but the programs causing such determinate heterogeneity are largely unexplored.

RESULTS

Here, we present MethylTransition, a novel DNA methylation state transition model, for characterizing methylation changes during one or a few cell cycles at single-cell resolution. MethylTransition involves the creation of a transition matrix comprising three parameters that represent the probabilities of DNA methylation-modifying activities in order to link the methylation states before and after a cell cycle. We apply MethylTransition to single-cell DNA methylome data from human pre-implantation embryogenesis and elucidate that the DNA methylation heterogeneity that emerges at promoters during this process is largely an intrinsic output of a program with unique probabilities of DNA methylation-modifying activities. Moreover, we experimentally validate the effect of the initial DNA methylation on expression heterogeneity in pre-implantation mouse embryos.

CONCLUSIONS

Our study reveals the programmed DNA methylation heterogeneity during human pre-implantation embryogenesis through a novel mathematical model and provides valuable clues for identifying the driving factors of the first cell fate determination during this process.