DNA methylation dynamics and dysregulation delineated by high-throughput profiling in the mouse

Cost-effective solutions for high-throughput enzymatic DNA methylation sequencing

Characterizing DNA methylation patterns is important for addressing key questions in evolutionary biology, development, geroscience, and medical genomics. While costs are decreasing, whole-genome DNA methylation profiling remains prohibitively expensive for most population-scale studies, creating a need for cost-effective, reduced representation approaches (i.e., assays that rely on microarrays, enzyme digests, or sequence capture to target a subset of the genome). Most common whole genome and reduced representation techniques rely on bisulfite conversion, which can damage DNA resulting in DNA loss and sequencing biases. Enzymatic methyl sequencing (EM-seq) was recently proposed to overcome these issues, but thorough benchmarking of EM-seq combined with cost-effective, reduced representation strategies is currently lacking. To address this gap, we optimized the Targeted Methylation Sequencing protocol (TMS)-which profiles ~4 million CpG sites-for miniaturization, flexibility, and multispecies use at a cost of ~USD 80. First, we tested modifications to increase throughput and reduce cost, including increasing multiplexing, decreasing DNA input, and using enzymatic rather than mechanical fragmentation to prepare DNA. Second, we compared our optimized TMS protocol to commonly used techniques, specifically the Infinium MethylationEPIC BeadChip (n = 55 paired samples) and whole genome bisulfite sequencing (n = 6 paired samples). In both cases, we found strong agreement between technologies (R2 = 0.97 and 0.99, respectively). Third, we tested the optimized TMS protocol in three non-human primate species (rhesus macaques, geladas, and capuchins). We captured a high percentage (mean = 77.1%) of targeted CpG sites and produced methylation level estimates that agreed with those generated from reduced representation bisulfite sequencing (R2 = 0.98). Finally, we confirmed that estimates of 1) epigenetic age and 2) tissue-specific DNA methylation patterns are strongly recapitulated using data generated from TMS versus other technologies. Altogether, our optimized TMS protocol will enable cost-effective, population-scale studies of genome-wide DNA methylation levels across human and non-human primate species.

Ecological Realism Accelerates Epigenetic Aging in Mice

The aging of mammalian epigenomes fundamentally alters cellular functions, and such changes are the focus of many healthspan and lifespan studies. However, studies of this process typically use mouse models living under standardized laboratory conditions and neglect the impact of variation in social, physical, microbial, and other aspects of the living environment on age-related changes. We examined differences in age-associated methylation changes between traditionally laboratory-reared mice from Jackson Laboratory and "rewilded" C57BL/6J mice, which lived in an outdoor field environment at Cornell University with enhanced ecological realism. Systematic analysis of age-associated methylation dynamics in the liver indicates a genomic region-conditioned, faster epigenetic aging rate in mice living in the field than those living in the lab, implicating perturbed 3D genome conformation and liver function. Altered epigenetic aging rates were more pronounced in sites that gain methylation with age, including sites enriched for transcription factor binding related to DNA repair. These observations underscore the overlooked role of the social and physical environment in epigenetic aging with implications for both basic and applied aging research.

Epigenetic ageing clocks: statistical methods and emerging computational challenges

Over the past decade, epigenetic clocks have emerged as powerful machine learning tools, not only to estimate chronological and biological age but also to assess the efficacy of anti-ageing, cellular rejuvenation and disease-preventive interventions. However, many computational and statistical challenges remain that limit our understanding, interpretation and application of epigenetic clocks. Here, we review these computational challenges, focusing on interpretation, cell-type heterogeneity and emerging single-cell methods, aiming to provide guidelines for the rigorous construction of interpretable epigenetic clocks at cell-type and single-cell resolution.

Tracking single-cell evolution using clock-like chromatin accessibility loci

Single-cell chromatin accessibility sequencing (scATAC-seq) reconstructs developmental trajectory by phenotypic similarity. However, inferring the exact developmental trajectory is challenging. Previous studies showed age-associated DNA methylation (DNAm) changes in specific genomic regions, termed clock-like differential methylation loci (ClockDML). Age-associated DNAm could either result from or result in chromatin accessibility changes at ClockDML. As cells undergo mitosis, the heterogeneity of chromatin accessibility on clock-like loci is reduced, providing a measure of mitotic age. In this study, we developed a method, called EpiTrace, that counts the fraction of opened clock-like loci from scATAC-seq data to determine cell age and perform lineage tracing in various cell lineages and animal species. It shows concordance with known developmental hierarchies, correlates well with DNAm-based clocks and is complementary with mutation-based lineage tracing, RNA velocity and stemness predictions. Applying EpiTrace to scATAC-seq data reveals biological insights with clinically relevant implications, ranging from hematopoiesis, organ development, tumor biology and immunity to cortical gyrification.

In vitro fertilization induces reproductive changes in male mouse offspring and has multigenerational effects

In vitro fertilization (IVF) is a noncoital method of conception used to treat human infertility. Although IVF is viewed as largely safe, it is associated with adverse outcomes in the fetus, placenta, and adult offspring. Because studies focusing on the effect of IVF on the male reproductive system are limited, we used a mouse model to assess the morphological and molecular effects of IVF on male offspring. We evaluated 3 developmental stages: 18.5-day fetuses and 12- and 39-week-old adults. Regardless of age, we observed changes in testicular-to-body weight ratios, serum testosterone levels, testicular morphology, gene expression, and DNA methylation. Also, sperm showed changes in morphology and DNA methylation. To assess multigenerational phenotypes, we mated IVF-conceived and naturally conceived males with wild-type females. Offspring from IVF males exhibited decreased fetal-to-placental weight ratios and changes in placenta gene expression and morphology regardless of sex. At 12 weeks of age, offspring showed higher body weights and differences in glucose, triglyceride, insulin, total cholesterol, HDL-C, and LDL/VLDL-C levels. Both sexes showed changes in gene expression in liver, testes, and ovaries and decreased global DNA methylation. Collectively, our findings demonstrate that male IVF offspring exhibit abnormal testicular and sperm morphology and molecular alterations with a multigenerational impact.

5-aza-2-deoxycytidine improves skeletal muscle function in a mouse model for recessive RYR1-related congenital myopathy

RYR1-related congenital myopathies are rare disorders that severely impair muscle function and the quality of life of patients and their families. To date no pharmacological therapies are available to treat the severe muscle weakness of affected patients. The most severe forms of RYR1-related congenital myopathies are caused by compound heterozygous mutations (nonsense/frameshift in one allele and a missense mutation in the other), leading to reduced RyR1 protein levels and altered biochemical composition of muscles. In this pre-clinical study, we treated a mouse model carrying the RyR1 p.Q1970fsX16 + p.A4329D compound heterozygous pathogenic variants (dHT mice) for 15 weeks with 0.05 mg/kg 5-aza-2'-deoxycytidine, an FDA-approved drug targeting DNA methyltransferases. We evaluated muscle strength, calcium homeostasis and muscle proteome and report that drug treatment improves all investigated parameters in dHT mice. Importantly, the beneficial effects were particularly significant in fast twitch muscles which are the first muscles to be impaired in patients. In conclusion, this study provides proof of concept for the pharmacological treatment of patients with recessive RYR1-related congenital myopathies with the FDA approved 5-aza-2'-deoxycytidine, supporting its use in a phase 1/2 clinical trial.

Epigenetic Changes in Cerebrospinal Fluid and Blood of People With Neurosyphilis

Epigenetic changes within immune cells may contribute to neuroinflammation during bacterial infection, but their role in neurosyphilis (NS) pathogenesis and response has not yet been established. We longitudinally analyzed DNA methylation and RNA expression in cerebrospinal fluid (CSF) cells and peripheral blood mononuclear cells (PBMCs) from 11 participants with laboratory-confirmed NS (CSF Venereal Disease Research Laboratory test positive) and 11 matched controls with syphilis without NS (non-NS). DNA methylation profiles from CSF and PBMCs of participants with NS significantly differed from those of participants with non-NS. Some genes associated with these differentially methylated sites had corresponding RNA expression changes in the CSF (111/1097 [10.1%]), and included genes involved in B cell activation and insulin-response pathways. Despite antibiotic treatment, approximately 80% of CSF methylation changes associated with NS persisted, suggesting that epigenetic scars accompanying NS may persistently affect immunity following infection. Future studies must examine whether these sequelae are clinically meaningful.

Epigenetic and Mitochondrial Metabolic Dysfunction in Multiple Sclerosis: A Review of Herbal Drug Approaches and Current Clinical Trials

Multiple sclerosis (MS) is a complex autoimmune disease characterised by inflammation, demyelination, and neurodegeneration within the central nervous system (CNS). While the exact causes remain unclear, recent research highlights the significant role of epigenetic modifications and mitochondrial dysfunction in the disease's onset and progression. Epigenetic alterations, such as DNA methylation, histone modification, and microRNA regulation, influence gene expression without altering the DNA sequence, leading to immune dysregulation and inflammation. Similarly, mitochondrial dysfunction, marked by impaired oxidative phosphorylation, reduced adenosine triphosphate (ATP) production, and increased reactive oxygen species (ROS), contributes to neurodegeneration and impaired remyelination in MS. The growing interest in targeting these two interconnected mechanisms has opened new avenues for MS treatment. Herbal drugs, known for their multi-targeted effects, have shown potential in modulating epigenetic markers and enhancing mitochondrial function. Compounds such as resveratrol, curcumin, epigallocatechin-3-gallate (EGCG), quercetin, and omega-3 fatty acids demonstrate potential in regulating DNA methylation, histone deacetylation, and mitochondrial biogenesis. These natural agents offer dual-action therapies by reducing oxidative stress and inflammation while promoting neuronal survival and remyelination. This review explores the therapeutic potential of herbal drugs targeting epigenetic and mitochondrial pathways in MS, evaluating their mechanisms of action and highlighting their promise as novel therapeutic agents. While initial findings are encouraging, further research and clinical trials are required to validate the efficacy of these herbal treatments and fully understand their potential in slowing disease progression and improving patient outcomes in MS. Such exploration could pave the way for safer, multi-targeted therapies, offering new hope in the management of MS and other neurodegenerative diseases.

Epigenetic signatures in surrogate tissues are able to assess cancer risk and indicate the efficacy of preventive measures

BACKGROUND

In order to advance personalized primary cancer prevention, surrogate endpoint biomarkers in distant, easy to access tissues (i.e., field defect indicators) reflecting field cancerization in the organ at risk are essential.

METHODS

Here we utilized medroxyprogesterone acetate and 7,12-dimethylbenzanthracene to induce mammary gland cancers in mice. We assessed epigenetic signatures reflective of carcinogen exposure, cell-type composition, mitotic age, and methylation at progesterone receptor binding sites in both, the tissue at risk (normal mammary gland; field cancerization) and distant non-at-risk organs (cervix, oviduct, and blood; field defect indicators), in mice that did and did not develop mammary gland cancers.

RESULTS

We demonstrate that the anti-progestine mifepristone reduces the cancer risk by more than 50%. Importantly, the reduction in cancer risk is accompanied by a decline in both field cancerization and field defect indicators; specifically, epigenetic signatures in the cervix are predictive of mammary cancer formation but show tissue-specific directionality.

CONCLUSIONS

These data encourage further exploration of epigenetic biomarkers in certain field defect-indicating tissues with a view to monitor the efficacy of cancer prevention strategies in humans.

DNA methylation patterns are influenced by Pax3::Foxo1 expression and developmental lineage in rhabdomyosarcoma tumours forming in genetically engineered mouse models

Rhabdomyosarcoma (RMS) is a family of phenotypically myogenic paediatric cancers consisting of two major subtypes: fusion-positive (FP) RMS, most commonly involving the PAX3::FOXO1 fusion gene, formed by the fusion of paired box 3 (PAX3) and forkhead box O1 (FOXO1) genes, and fusion-negative (FN) RMS, lacking these gene fusions. In humans, DNA methylation patterns distinguish these two subtypes as well as mutation-associated subsets within these subtypes. To investigate the biological factors responsible for these methylation differences, we profiled DNA methylation in RMS tumours derived from genetically engineered mouse models (GEMMs) in which various driver mutations were introduced into different myogenic lineages. Our unsupervised analyses of DNA methylation patterns in these GEMM tumours yielded two major clusters, corresponding to high and no/low expression of Pax3::Foxo1, which mirrored the results for human FP and FN RMS tumours. Two distinct methylation-defined subsets were found for GEMM RMS tumours with no/low Pax3::Foxo1 expression: one subset enriched in Pax7 lineage tumours and a second subset enriched in myogenic factor 5 (Myf5) lineage tumours. Integrative analysis of DNA methylation and transcriptomic data in mouse and human RMS revealed a common group of differentially methylated and differentially expressed genes, highlighting a conserved set of genes functioning in both human RMS models and GEMMs of RMS. In conclusion, these studies provide insight into the roles of oncogenic fusion proteins and developmental lineages in establishing DNA methylation patterns in FP and FN RMS respectively. © 2025 The Author(s). The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland. This article has been contributed to by U.S. Government employees and their work is in the public domain in the USA.

Scalable Screening of Ternary-Code DNA Methylation Dynamics Associated with Human Traits

Epigenome-wide association studies (EWAS) are transforming our understanding of the interplay between epigenetics and complex human traits and phenotypes. We introduce the Methylation Screening Array (MSA), a new iteration of the Infinium technology for scalable and quantitative screening of trait associations of nuanced ternary-code cytosine modifications in larger, more inclusive, and stratified human populations. MSA integrates EWAS, single-cell, and cell-type-resolved methylome profiles, covering diverse human traits and diseases. Our first MSA applications yield multiple biological insights: we revealed a previously unappreciated role of 5-hydroxymethylcytosine (5hmC) in trait associations and epigenetic clocks. We demonstrated that 5hmCs complement 5-methylcytosines (5mCs) in defining tissues and cells' epigenetic identities. In-depth analyses highlighted the cell type context of EWAS and GWAS hits. Using this platform, we conducted a comprehensive human 5hmC aging EWAS, discovering tissue-invariant and tissue-specific aging dynamics, including distinct tissue-specific rates of mitotic hyper- and hypomethylation rates. These findings chart a landscape of the complex interplay of the two forms of cytosine modifications in diverse human tissues and their roles in health and disease.

Parkinson's disease-associated alterations in DNA methylation and hydroxymethylation in human brain

Epigenetic mechanisms are mediators of interactions between aging, genetics, and environmental factors in sporadic Parkinson's disease (PD). Multiple studies have explored the DNA modifications in PD, but few focus on 5-hydroxymethylcytosine (5-hmC), which is important in the central nervous system and sensitive to environmental exposures. To date, studies have not differentiated between 5-methylcytosine (5-mC) and 5-hmC or have analyzed them separately. In this study, we modeled paired 5-mC and 5-hmC data simultaneously. We identified 108 cytosines with significant PD-associated shifts between these marks in an enriched neuronal population from PD postmortem parietal cortex, within 83 genes and 34 enhancers associated with 67 genes. These data potentially link epigenetic regulation of genes related to LRRK2 and endolysosomal sort (RAB32 and AGAP1), and genes involved in neuroinflammation, the inflammasome, and neurodevelopment with early changes in PD and suggest that there are significant shifts between 5mC and 5hmC associated with PD in genes not captured by standard methods.

TRIM28-dependent developmental heterogeneity determines cancer susceptibility through distinct epigenetic states

Mutations in cancer risk genes increase susceptibility, but not all carriers develop cancer. Indeed, while DNA mutations are necessary drivers of cancer, only a small subset of mutated cells go on to cause the disease. To date, the mechanisms underlying individual cancer susceptibility remain unclear. Here, we took advantage of a unique mouse model of intrinsic developmental heterogeneity (Trim28+/D9) to investigate whether early-life epigenetic variation influences cancer susceptibility later in life. We found that heterozygosity of Trim28 is sufficient to generate two distinct early-life epigenetic states associated with differing cancer susceptibility. These developmentally primed states exhibit differential methylation patterns at typically silenced heterochromatin, detectable as early as 10 days of age. The differentially methylated loci are enriched for genes with known oncogenic potential, frequently mutated in human cancers and correlated with poor prognosis. This study provides genetic evidence that intrinsic developmental heterogeneity can prime individual, lifelong cancer susceptibility.

On correlative and causal links of replicative epimutations

The mitotic inheritability of DNA methylation as an epigenetic marker in higher-order eukaryotes has been established for >40 years. The DNA methylome and mitotic division interplay is now considered bidirectional and highly intertwined. Various epigenetic writers, erasers, and modulators shape the perceived replicative methylation dynamics. This Review surveys the principles and complexity of mitotic transmission of DNA methylation, emphasizing the awareness of mitotic aging in analyzing DNA methylation dynamics in development and disease. We reviewed how DNA methylation changes alter mitotic proliferation capacity, implicating age-related diseases like cancer. We link replicative epimutation to stem cell dysfunction, inflammatory response, cancer risks, and epigenetic clocks, discussing the causative role of DNA methylation in health and disease.

TET1 displays catalytic and non-catalytic functions in the adult mouse cortex

TET1/2/3 dioxygenases iteratively demethylate 5-methylcytosine, beginning with the formation of 5-hydroxymethylcytosine (5hmC). The post-mitotic brain maintains higher levels of 5hmC than most peripheral tissues, and TET1 ablation studies have underscored the critical role of TET1 in brain physiology. However, deletion of Tet1 precludes the disentangling of the catalytic and non-catalytic functions of TET1. Here, we dissect these functions of TET1 by comparing adult cortex of Tet1 wildtype (Tet1 WT), a novel Tet1 catalytically dead mutant (Tet1 HxD), and Tet1 knockout (Tet1 KO) mice. Using DNA methylation array, we uncover that Tet1 HxD and KO mutations perturb the methylation status of distinct subsets of CpG sites. Gene ontology (GO) analysis on specific differential 5hmC regions indicates that TET1's catalytic activity is linked to neuronal-specific functions. RNA-Seq further shows that Tet1 mutations predominantly impact the genes that are associated with alternative splicing. Lastly, we performed High-performance Liquid Chromatography Mass-Spectrometry lipidomics on WT and mutant cortices and uncover accumulation of lysophospholipids lysophosphatidylethanolamine and lysophosphatidylcholine in Tet1 HxD cortex. In summary, we show that Tet1 HxD does not completely phenocopy Tet1 KO, providing evidence that TET1 modulates distinct cortical functions through its catalytic and non-catalytic roles.

Novel DNA methylation changes in mouse lungs associated with chronic smoking

Smoking is a potent cause of asthma exacerbations, chronic obstructive pulmonary disease (COPD) and many other health defects, and changes in DNA methylation (DNAm) have been identified as a potential link between smoking and these health outcomes. However, most studies of smoking and DNAm have been done using blood and other easily accessible tissues in humans, while evidence from more directly affected tissues such as the lungs is lacking. Here, we identified DNAm patterns in the lungs that are altered by smoking. We used an established mouse model to measure the effects of chronic smoke exposure first on lung phenotype immediately after smoking and then after a period of smoking cessation. Next, we determined whether our mouse model recapitulates previous DNAm patterns observed in smoking humans, specifically measuring DNAm at a candidate gene responsive to cigarette smoke, Cyp1a1. Finally, we carried out epigenome-wide DNAm analyses using the newly released Illumina mouse methylation microarrays. Our results recapitulate some of the phenotypes and DNAm patterns observed in human studies but reveal 32 differentially methylated genes specific to the lungs which have not been previously associated with smoking. The affected genes are associated with nicotine dependency, tumorigenesis and metastasis, immune cell dysfunction, lung function decline, and COPD. This research emphasizes the need to study CS-mediated DNAm signatures in directly affected tissues like the lungs, to fully understand mechanisms underlying CS-mediated health outcomes.

Implementation of the Methyl-Seq platform to identify tissue- and sex-specific DNA methylation differences in the rat epigenome

Epigenomic annotations for the rat lag far behind those of human and mouse, despite the rat's immense utility in pharmacological and behavioral studies and the need to understand their epigenetic mechanisms. We have designed a targeted-enrichment method followed by next-generation sequencing (Methyl-Seq) to identify DNA methylation (DNAm) signatures across the rat genome. The design reflected an attempt to create a more comprehensive investigation of the rat epigenome, as it included promoters, CpG islands, and island shores of all RefSeq genes. In this study, we implemented the rat Methyl-Seq platform and tested its ability to distinguish differentially methylated regions (DMRs) among three different tissue types, three distinct brain regions, and, in the hippocampus, between males and females. These comparisons yielded DNAm differences of differing magnitudes, many of which were independently validated by bisulfite pyrosequencing, including autosomal regions that were predicted to show the least degree of difference in DNAm between males and females. Quantitative reverse transcription PCR revealed that most genes associated with the DMRs showed tissue-, brain region-, and sex-specific differences in expression. In particular, we found evidence for sex-specific DNAm and expression differences at Tubb6, Lrrn2, Tex26, and Sox5l1, all of which play important roles in neurodevelopment and have been implicated in studies examining sex differences. Our results demonstrate the utility of the rat Methyl-Seq platform and suggest the presence of DNAm differences between the male and female hippocampus. The rat Methyl-Seq has the potential to provide epigenomic insights into pharmacological and behavioral studies performed in the rat.

Genetic confounds of transgenerational epigenetic inheritance in mice

Transgenerational epigenetic inheritance in mammals remains a controversial phenomenon. A recent study by Takahashi et al. provides evidence for this mode of inheritance in mice by using a CRISPR/Cas9-based epigenetic editing technique to modify DNA methylation levels at specific promoters and then demonstrating the inheritance of the gain in methylation in offspring. In this technical commentary, we argue that the method used in the original study inherently amplifies the likelihood of genetic changes that thereafter lead to the heritability of epigenetic changes. We provide evidence that genetic changes from multiple sources do indeed occur in these experiments and explore several avenues by which these changes could be causal to the apparent inheritance of epigenetic changes. We conclude a genetic basis of inheritance cannot be ruled out and thus transgenerational epigenetic inheritance has not been adequately established by the original study.

Multi-strain probiotic administration decreases necrotizing enterocolitis severity and alters the epigenetic profile in mice

BACKGROUND

Probiotic administration may decrease the incidence of necrotizing enterocolitis (NEC) through mechanisms that are largely unknown. We investigated the effects of probiotics on intestinal epigenetics and assessed their effects on intestinal inflammation and motility using both ileum-predominant and combined ileo-colitis mouse NEC models.

METHODS

C57BL/6 J mice were gavage-fed a multi-strain probiotic from postnatal days 3-11, consisting of B. infantis, B. lactis, and S. thermophilus. From p8, mice were exposed to ileo-colitis NEC involving formula containing NEC bacteria and 0.5% DSS. DNA methylation was measured using the Infinium Methylation Assay. Gastrointestinal motility was assessed by 70 Kd FITC-dextran transit time. Probiotic colonization was measured in probiotic-fed mice by qPCR.

RESULTS

Probiotic administration caused significant changes in the small intestine's epigenetic signature, a reduction in NEC severity, and improved intestinal motility. The effects of probiotics were more pronounced in the ileo-colitis NEC model.

CONCLUSIONS

These findings shed light on the role of probiotics in two clinically relevant models of NEC, add additional insights into their underlying mechanism of action, and reveal unanticipated epigenetic modifications to the intestinal mucosa after their use.

IMPACT

These findings shed light on the role of multi-strain probiotics in two clinically relevant animal models of NEC, and add additional insights into their underlying mechanism of action This study provides a new, clinically relevant model for the study of NEC including administration of 0.5% DSS, to include ileal dominant and ileo-colonic dominant phenotypes of the disease. These results reveal that clinically relevant strains of probiotic bacteria can exert epigenetic effects on the small intestine in mice, and can attenuate the epigenetic changes induced by NEC.

Fetal programming by the parental microbiome of offspring behavior, and DNA methylation and gene expression within the hippocampus

Background

The microorganisms colonizing the gastrointestinal tract of animals, collectively referred to as the gut microbiome, affect numerous host behaviors dependent on the central nervous system (CNS). Studies comparing germ-free mice to normally colonized mice have demonstrated influences of the microbiome on anxiety-related behaviors, voluntary activity, and gene expression in the CNS. Additionally, there is epidemiologic evidence supporting an intergenerational influence of the maternal microbiome on neurodevelopment of offspring and behavior later in life. There is limited experimental evidence however directly linking the maternal microbiome to long-term neurodevelopmental outcomes, or knowledge regarding mechanisms responsible for such effects.

Results

Here we show that that the maternal microbiome has a dominant influence on several offspring phenotypes including anxiety-related behavior, voluntary activity, and body weight. Adverse outcomes in offspring were associated with features of the maternal microbiome including bile salt hydrolase activity gene expression (bsh), abundance of certain bile acids, and hepatic expression of Slc10a1. In cross-foster experiments, offspring resembled their birth dam phenotypically, despite faithful colonization in the postnatal period with the surrogate dam microbiome. Genome-wide methylation analysis of hippocampal DNA identified microbiome-associated differences in methylation of 196 loci in total, 176 of which show conserved profiles between mother and offspring. Further, single-cell transcriptional analysis revealed accompanying differences in expression of several differentially methylated genes within certain hippocampal cell clusters, and vascular expression of genes associated with bile acid transport. Inferred cell-to-cell communication in the hippocampus based on coordinated ligand-receptor expression revealed differences in expression of neuropeptides associated with satiety.

Conclusions

Collectively, these data provide proof-of-principle that the maternal gut microbiome has a dominant influence on the neurodevelopment underlying certain offspring behaviors and activities, and selectively affects genome methylation and gene expression in the offspring CNS in conjunction with that neurodevelopment.

An in vitro approach reveals molecular mechanisms underlying endocrine disruptor-induced epimutagenesis

Endocrine disrupting chemicals (EDCs) such as bisphenol S (BPS) are xenobiotic compounds that can disrupt endocrine signaling due to steric similarities to endogenous hormones. EDCs have been shown to induce disruptions in normal epigenetic programming (epimutations) and differentially expressed genes (DEGs) that predispose disease states. Most interestingly, the prevalence of epimutations following exposure to many EDCs persists over multiple generations. Many studies have described direct and prolonged effects of EDC exposure in animal models, but many questions remain about molecular mechanisms by which EDC-induced epimutations are introduced or subsequently propagated, whether there are cell type-specific susceptibilities to the same EDC, and whether this correlates with differential expression of relevant hormone receptors. We exposed cultured pluripotent (iPS), somatic (Sertoli and granulosa), and primordial germ cell-like (PGCLC) cells to BPS and found that differential incidences of BPS-induced epimutations and DEGs correlated with differential expression of relevant hormone receptors inducing epimutations near relevant hormone response elements in somatic and pluripotent, but not germ cell types. Most interestingly, we found that when iPS cells were exposed to BPS and then induced to differentiate into PGCLCs, the prevalence of epimutations and DEGs was largely retained, however, >90% of the specific epimutations and DEGs were replaced by novel epimutations and DEGs. These results suggest a unique mechanism by which an EDC-induced epimutated state may be propagated transgenerationally.

The scaffolding function of LSD1 controls DNA methylation in mouse ESCs

Lysine-specific histone demethylase 1 (LSD1), which demethylates mono- or di- methylated histone H3 on lysine 4 (H3K4me1/2), is essential for early embryogenesis and development. Here we show that LSD1 is dispensable for mouse embryonic stem cell (ESC) self-renewal but is required for mouse ESC growth and differentiation. Reintroduction of a catalytically-impaired LSD1 (LSD1MUT) recovers the proliferation capability of mouse ESCs, yet the enzymatic activity of LSD1 is essential to ensure proper differentiation. Indeed, increased H3K4me1 in Lsd1 knockout (KO) mouse ESCs does not lead to major changes in global gene expression programs related to stemness. However, ablation of LSD1 but not LSD1MUT results in decreased DNMT1 and UHRF1 proteins coupled to global hypomethylation. We show that both LSD1 and LSD1MUT control protein stability of UHRF1 and DNMT1 through interaction with HDAC1 and the ubiquitin-specific peptidase 7 (USP7), consequently, facilitating the deacetylation and deubiquitination of DNMT1 and UHRF1. Our studies elucidate a mechanism by which LSD1 controls DNA methylation in mouse ESCs, independently of its lysine demethylase activity.

The interferon γ pathway enhances pluripotency and X-chromosome reactivation in iPSC reprogramming

Reprogramming somatic cells into induced pluripotent stem cells (iPSCs) requires activation of the pluripotency network and resetting of the epigenome by erasing the epigenetic memory of the somatic state. In female mouse cells, a critical epigenetic reprogramming step is the reactivation of the inactive X chromosome. Despite its importance, a systematic understanding of the regulatory networks linking pluripotency and X-reactivation is missing. Here, we reveal important pathways for pluripotency acquisition and X-reactivation using a genome-wide CRISPR screen during neural precursor to iPSC reprogramming. In particular, we discover that activation of the interferon γ (IFNγ) pathway early during reprogramming accelerates pluripotency acquisition and X-reactivation. IFNγ stimulates STAT3 signaling and the pluripotency network and leads to enhanced TET-mediated DNA demethylation, which consequently boosts X-reactivation. We therefore gain a mechanistic understanding of the role of IFNγ in reprogramming and X-reactivation and provide a comprehensive resource of the molecular networks involved in these processes.

PRC2-AgeIndex as a universal biomarker of aging and rejuvenation

DNA methylation (DNAm) is one of the most reliable biomarkers of aging across mammalian tissues. While the age-dependent global loss of DNAm has been well characterized, DNAm gain is less characterized. Studies have demonstrated that CpGs which gain methylation with age are enriched in Polycomb Repressive Complex 2 (PRC2) targets. However, whole-genome examination of all PRC2 targets as well as determination of the pan-tissue or tissue-specific nature of these associations is lacking. Here, we show that low-methylated regions (LMRs) which are highly bound by PRC2 in embryonic stem cells (PRC2 LMRs) gain methylation with age in all examined somatic mitotic cells. We estimated that this epigenetic change represents around 90% of the age-dependent DNAm gain genome-wide. Therefore, we propose the "PRC2-AgeIndex," defined as the average DNAm in PRC2 LMRs, as a universal biomarker of cellular aging in somatic cells which can distinguish the effect of different anti-aging interventions.

Parkinson's disease-associated shifts between DNA methylation and DNA hydroxymethylation in human brain in PD-related genes, including PARK19 (DNAJC6) and PTPRN2 (IA-2β)

Background

The majority of Parkinson's disease (PD) cases are due to a complex interaction between aging, genetics, and environmental factors; epigenetic mechanisms are thought to act as important mediators of these risk factors. While multiple studies to date have explored the role of DNA modifications in PD, few focus on 5-hydroxymethylcytosine (5hmC). Because 5hmC occurs at its highest levels in the brain and is thought to be particularly important in the central nervous system, particularly in the response to neurotoxicants, it is important to explore the potential role of 5hmC in PD. This study expands on our previously published epigenome-wide association study (EWAS) performed on DNA isolated from neuron-enriched nuclei from human postmortem parietal cortex from the Banner Sun Health Research Institute Brain Bank. The study aimed to identify paired changes in 5hmC and 5mC in PD in enriched neuronal nuclei isolated from PD post-mortem parietal cortex and age- and sex-matched controls. We performed oxidative bisulfite (oxBS) conversion and paired it with our previously published bisulfite (BS)-based EWAS on the same samples to identify cytosines with significant shifts between these two related epigenetic marks. Interaction differentially modified cytosines (iDMCs) were identified using our recently published mixed-effects model for co-analyzing βmC and βhmC data.

Results

We identified 1,030 iDMCs with paired changes in 5mC and 5hmC (FDR < 0.05) that map to 695 genes, including PARK19 (DNAJC6), a familial PD gene, and PTPRN2 (IA-2), which has been previously implicated in PD in both epigenetic and mechanistic studies. The majority of iDMC-containing genes have not previously been implicated in PD and were not identified in our previous BS-based EWAS.

Conclusions

These data potentially link epigenetic regulation of the PARK19 and PTPRN2 loci in the pathogenesis of idiopathic PD. In addition, iDMC-containing genes have known functions in synaptic formation and function, cell cycle and senescence, neuroinflammation, and epigenetic regulation. These data suggest that there are significant shifts between 5mC and 5hmC associated with PD in genes relevant to PD pathogenesis that are not captured by analyzing BS-based data alone or by analyzing each mark as a distinct dataset.

Homeobox and Polycomb target gene methylation in human solid tumors

DNA methylation is an epigenetic mark that plays an important role in defining cancer phenotypes, with global hypomethylation and focal hypermethylation at CpG islands observed in tumors. These methylation marks can also be used to define tumor types and provide an avenue for biomarker identification. The homeobox gene class is one that has potential for this use, as well as other genes that are Polycomb Repressive Complex 2 targets. To begin to unravel this relationship, we performed a pan-cancer DNA methylation analysis using sixteen Illumina HM450k array datasets from TCGA, delving into cancer-specific qualities and commonalities between tumor types with a focus on homeobox genes. Our comparisons of tumor to normal samples suggest that homeobox genes commonly harbor significant hypermethylated differentially methylated regions. We identified two homeobox genes, HOXA3 and HOXD10, that are hypermethylated in all 16 cancer types. Furthermore, we identified several potential homeobox gene biomarkers from our analysis that are uniquely methylated in only one tumor type and that could be used as screening tools in the future. Overall, our study demonstrates unique patterns of DNA methylation in multiple tumor types and expands on the interplay between the homeobox gene class and oncogenesis.

Endocrine disruptor-induced epimutagenesis in vitro : Insight into molecular mechanisms

Endocrine disrupting chemicals (EDCs) such as bisphenol S (BPS) are xenobiotic compounds that can disrupt endocrine signaling following exposure due to steric similarities to endogenous hormones within the body. EDCs have been shown to induce disruptions in normal epigenetic programming (epimutations) that accompany dysregulation of normal gene expression patterns that appear to predispose disease states. Most interestingly, the prevalence of epimutations following exposure to many different EDCs often persists over multiple subsequent generations, even with no further exposure to the causative EDC. Many previous studies have described both the direct and prolonged effects of EDC exposure in animal models, but many questions remain about molecular mechanisms by which EDCs initially induce epimutations or contribute to the propagation of EDC-induced epimutations either within the exposed generation or to subsequent generations. Additional questions remain regarding the extent to which there may be differences in cell-type specific susceptibilities to various EDCs, and whether this susceptibility is correlative with expression of relevant hormone receptors and/or the location of relevant hormone response elements (HREs) in the genome. To address these questions, we exposed cultured mouse pluripotent (induced pluripotent stem [iPS]), somatic (Sertoli and granulosa), and germ (primordial germ cell like [PGCLC]) cells to BPS and measured changes in DNA methylation levels at the epigenomic level and gene expression at the transcriptomic level. We found that there was indeed a difference in cell-type specific susceptibility to EDC-induced epimutagenesis and that this susceptibility correlated with differential expression of relevant hormone receptors and, in many cases, tended to generate epimutations near relevant HREs within the genome. Additionally, however, we also found that BPS can induce epimutations in a cell type that does not express relevant receptors and in genomic regions that do not contain relevant HREs, suggesting that both canonical and non-canonical signaling mechanisms can be disrupted by BPS exposure. Most interestingly, we found that when iPS cells were exposed to BPS and then induced to differentiate into PGCLCs, the prevalence of epimutations and differentially expressed genes (DEGs) initially induced in the iPSCs was largely retained in the resulting PGCLCs, however, >90% of the specific epimutations and DEGs were not conserved but were rather replaced by novel epimutations and DEGs following the iPSC to PGCLC transition. These results are consistent with a unique concept that many EDC-induced epimutations may normally be corrected by germline and/or embryonic epigenetic reprogramming but that due to disruption of the underlying chromatin architecture induced by the EDC exposure, many novel epimutations may emerge during the reprogramming process as well. Thus, it appears that following exposure to a disruptive agent such as an EDC, a prevalence of epimutations may transcend epigenetic reprogramming even though most individual epimutations are not conserved during this process.

Changes in mouse epidermal DNA methylation during development of squamous cell carcinoma in response to UVR

Squamous cell carcinoma (SCC) is a common skin cancer, often caused by exposure to ultraviolet radiation (UVR). Recent studies have shown that changes in DNA methylation play a crucial role in the development of cancers. However, methylation patterns of SCC are not well characterised. Identifying biomarkers for the risk of developing SCC could be helpful for early detection and diagnosis and can potentially improve treatment and prevention strategies. This study aimed to investigate methylation changes in the epidermis of mice exposed to UVR for 24 weeks. We examined the DNA methylation levels of 260 199 CpGs using the Illumina Infinium Mouse Methylation BeadChip and studied the epidermis of UVR-exposed and unexposed mice every 4 weeks for 24 weeks (n = 39). We identified CpGs with large differences in methylation levels (β-values) between UVR-exposed and unexposed mice. We also observed differences in the epigenetic age of these mice. We identified CpGs in Rev, Ipmk, Rad51b, Fgfr2, Fgfr3 and Ctnnb1 that may serve as potential biomarkers for SCC risk and could be helpful for the early detection and prevention of SCC. Further investigations are necessary to determine the biological functions and clinical significance of these CpGs.

Improve-RRBS: a novel tool to correct the 3' trimming of reduced representation sequencing reads

Motivation

Reduced Representation Bisulfite Sequencing (RRBS) is a popular approach to determine DNA methylation of the CpG-rich regions of the genome. However, we observed that false positive differentially methylated sites (DMS) are also identified using the standard computational analysis.

Results

During RRBS library preparation the MspI digested DNA undergo end-repair by a cytosine at the 3' end of the fragments. After sequencing, Trim Galore cuts these end-repaired nucleotides. However, Trim Galore fails to detect end-repair when it overlaps with the 3' end of the sequencing reads. We found that these non-trimmed cytosines bias methylation calling, thus, can identify DMS erroneously. To circumvent this problem, we developed improve-RRBS, which efficiently identifies and hides these cytosines from methylation calling with a false positive rate of maximum 0.5%. To test improve-RRBS, we investigated four datasets from four laboratories and two different species. We found non-trimmed 3' cytosines in all datasets analyzed and as much as >50% of false positive DMS under certain conditions. By applying improve-RRBS, these DMS completely disappeared from all comparisons.

Availability and implementation

Improve-RRBS is a freely available python package https://pypi.org/project/iRRBS/ or https://github.com/fothia/improve-RRBS to be implemented in RRBS pipelines.

Increased hippocampal epigenetic age in the Ts65Dn mouse model of Down Syndrome

Down syndrome (DS) is a segmental progeroid genetic disorder associated with multi-systemic precocious aging phenotypes, which are particularly evident in the immune and nervous systems. Accordingly, people with DS show an increased biological age as measured by epigenetic clocks. The Ts65Dn trisomic mouse, which harbors extra-numerary copies of chromosome 21 (Hsa21)-syntenic regions, was shown to recapitulate several progeroid features of DS, but no biomarkers of age have been applied to it so far. In this pilot study, we used a mouse-specific epigenetic clock to measure the epigenetic age of hippocampi from Ts65Dn and euploid mice at 20 weeks. Ts65Dn mice showed an increased epigenetic age in comparison with controls, and the observed changes in DNA methylation partially recapitulated those observed in hippocampi from people with DS. Collectively, our results support the use of the Ts65Dn model to decipher the molecular mechanisms underlying the progeroid DS phenotypes.

Low-input and single-cell methods for Infinium DNA methylation BeadChips

The Infinium BeadChip is the most widely used DNA methylome assay technology for population-scale epigenome profiling. However, the standard workflow requires over 200 ng of input DNA, hindering its application to small cell-number samples, such as primordial germ cells. We developed experimental and analysis workflows to extend this technology to suboptimal input DNA conditions, including ultra-low input down to single cells. DNA preamplification significantly enhanced detection rates to over 50% in five-cell samples and ∼25% in single cells. Enzymatic conversion also substantially improved data quality. Computationally, we developed a method to model the background signal's influence on the DNA methylation level readings. The modified detection P-value calculation achieved higher sensitivities for low-input datasets and was validated in over 100 000 public diverse methylome profiles. We employed the optimized workflow to query the demethylation dynamics in mouse primordial germ cells available at low cell numbers. Our data revealed nuanced chromatin states, sex disparities, and the role of DNA methylation in transposable element regulation during germ cell development. Collectively, we present comprehensive experimental and computational solutions to extend this widely used methylation assay technology to applications with limited DNA.

TIME-seq reduces time and cost of DNA methylation measurement for epigenetic clock construction

Epigenetic 'clocks' based on DNA methylation have emerged as the most robust and widely used aging biomarkers, but conventional methods for applying them are expensive and laborious. Here we develop tagmentation-based indexing for methylation sequencing (TIME-seq), a highly multiplexed and scalable method for low-cost epigenetic clocks. Using TIME-seq, we applied multi-tissue and tissue-specific epigenetic clocks in over 1,800 mouse DNA samples from eight tissue and cell types. We show that TIME-seq clocks are accurate and robust, enriched for polycomb repressive complex 2-regulated loci, and benchmark favorably against conventional methods despite being up to 100-fold less expensive. Using dietary treatments and gene therapy, we find that TIME-seq clocks reflect diverse interventions in multiple tissues. Finally, we develop an economical human blood clock (R > 0.96, median error = 3.39 years) in 1,056 demographically representative individuals. These methods will enable more efficient epigenetic clock measurement in larger-scale human and animal studies.

Trophectoderm biopsy of blastocysts following IVF and embryo culture increases epigenetic dysregulation in a mouse model

STUDY QUESTION

Does trophectoderm biopsy (TEBx) of blastocysts for preimplantation genetic testing in the clinic affect normal placental and embryo development and offspring metabolic outcomes in a mouse model?

SUMMARY ANSWER

TEBx impacts placental and embryonic health during early development, with some alterations resolving and others worsening later in development and triggering metabolic changes in adult offspring.

WHAT IS KNOWN ALREADY

Previous studies have not assessed the epigenetic and morphological impacts of TEBx either in human populations or in animal models.

STUDY DESIGN, SIZE, DURATION

We employed a mouse model to identify the effects of TEBx during IVF. Three groups were assessed: naturally conceived (Naturals), IVF, and IVF + TEBx, at two developmental timepoints: embryonic day (E)12.5 (n = 40/Naturals, n = 36/IVF, and n = 36/IVF + TEBx) and E18.5 (n = 42/Naturals, n = 30/IVF, and n = 35/IVF + TEBx). Additionally, to mimic clinical practice, we assessed a fourth group: IVF + TEBx + Vitrification (Vit) at E12.5 (n = 29) that combines TEBx and vitrification. To assess the effect of TEBx in offspring health, we characterized a 12-week-old cohort (n = 24/Naturals, n = 25/IVF and n = 25/IVF + TEBx).

PARTICIPANTS/MATERIALS, SETTING, METHODS

Our mouse model used CF-1 females as egg donors and SJL/B6 males as sperm donors. IVF, TEBx, and vitrification were performed using standardized methods. Placenta morphology was evaluated by hematoxylin-eosin staining, in situ hybridization using Tpbpa as a junctional zone marker and immunohistochemistry using CD34 fetal endothelial cell markers. For molecular analysis of placentas and embryos, DNA methylation was analyzed using pyrosequencing, luminometric methylation assay, and chip array technology. Expression patterns were ascertained by RNA sequencing. Triglycerides, total cholesterol, high-, low-, and very low-density lipoprotein, insulin, and glucose were determined in the 12-week-old cohort using commercially available kits.

MAIN RESULTS AND THE ROLE OF CHANCE

We observed that at E12.5, IVF + TEBx had a worse outcome in terms of changes in DNA methylation and differential gene expression in placentas and whole embryos compared with IVF alone and compared with Naturals. These changes were reflected in alterations in placental morphology and blood vessel density. At E18.5, early molecular changes in fetuses were maintained or exacerbated. With respect to placentas, the molecular and morphological changes, although different compared to Naturals, were equivalent to the IVF group, except for changes in blood vessel density, which persisted. Of note is that most differences were sex specific. We conclude that TEBx has more detrimental effects in mid-gestation placental and embryonic tissues, with alterations in embryonic tissues persisting or worsening in later developmental stages compared to IVF alone, and the addition of vitrification after TEBx results in more pronounced and potentially detrimental epigenetic effects: these changes are significantly different compared to Naturals. Finally, we observed that 12-week IVF + TEBx offspring, regardless of sex, showed higher glucose, insulin, triglycerides, lower total cholesterol, and lower high-density lipoprotein compared to IVF and Naturals, with only males having higher body weight compared to IVF and Naturals. Our findings in a mouse model additionally support the need for more studies to assess the impact of new procedures in ART to ensure healthy pregnancies and offspring outcomes.

LARGE SCALE DATA

Data reported in this work have been deposited in the NCBI Gene Expression Omnibus under accession number GSE225318.

LIMITATIONS, REASONS FOR CAUTION

This study was performed using a mouse model that mimics many clinical IVF procedures and outcomes observed in humans, where studies on early embryos are not possible.

WIDER IMPLICATIONS OF THE FINDINGS

This study highlights the importance of assaying new procedures used in ART to assess their impact on placenta and embryo development, and offspring metabolic outcomes.

STUDY FUNDING/COMPETING INTEREST(S)

This work was funded by a National Centers for Translational Research in Reproduction and Infertility grant P50 HD068157-06A1 (M.S.B., C.C., M.M.), Ruth L. Kirschstein National Service Award Individual Postdoctoral Fellowship F32 HD107914 (E.A.R.-C.) and F32 HD089623 (L.A.V.), and National Institutes of Health Training program in Cell and Molecular Biology T32 GM007229 (C.N.H.). No conflict of interest.

Epigenetic regulation of drug metabolism in aging: utilizing epigenetics to optimize geriatric pharmacotherapy

We explore the relationship between epigenetic aging and drug metabolism. We review current evidence for changes in drug metabolism in normal aging, followed by a description of how epigenetic modifications associated with age can regulate the expression and functionality of genes. In particular, we focus on the role of epigenome-wide studies of human and mouse liver in understanding these age-related processes with respect to xenobiotic processing. We highlight genes encoding drug metabolizing enzymes and transporters revealed to be affected by epigenetic aging in these studies. We conclude that substantial evidence exists for epigenetic aging impacting drug metabolism and transport genes, but more work is needed. We further highlight the promise of pharmacoepigenetics applied to enhancing drug safety in older adults.

Single-cell methylation analysis of brain tissue prioritizes mutations that alter transcription

Relating genetic variants to behavior remains a fundamental challenge. To assess the utility of DNA methylation marks in discovering causative variants, we examined their relationship to genetic variation by generating single-nucleus methylomes from the hippocampus of eight inbred mouse strains. At CpG sequence densities under 40 CpG/Kb, cells compensate for loss of methylated sites by methylating additional sites to maintain methylation levels. At higher CpG sequence densities, the exact location of a methylated site becomes more important, suggesting that variants affecting methylation will have a greater effect when occurring in higher CpG densities than in lower. We found this to be true for a variant's effect on transcript abundance, indicating that candidate variants can be prioritized based on CpG sequence density. Our findings imply that DNA methylation influences the likelihood that mutations occur at specific sites in the genome, supporting the view that the distribution of mutations is not random.

Comparison of global DNA methylation analysis by whole genome bisulfite sequencing and the Infinium Mouse Methylation BeadChip using fresh and fresh-frozen mouse epidermis

Until recently, studying the murine methylome was restricted to sequencing-based methods. In this study we compared the global DNA methylation levels of hairless mouse epidermis using the recently released Infinium Mouse Methylation BeadChip from Illumina and whole genome bisulphite sequencing (WGBS). We also studied the effect of sample storage conditions by using fresh and fresh-frozen epidermis. The DNA methylation levels of 123,851 CpG sites covered by both the BeadChip and WGBS were compared. DNA methylation levels obtained with WGBS and the BeadChip were strongly correlated (Pearson correlation r = 0.984). We applied a threshold of 15 reads for the WGBS methylation analysis. Even at a threshold of 10 reads, we observed no substantial difference in DNA methylation levels compared with that obtained with the BeadChip. The DNA methylation levels from the fresh and the fresh-frozen samples were strongly correlated when analysed with both the BeadChip (r = 0.999) and WGBS (r = 0.994). We conclude that the two methods of analysis generally work equally well for studies of DNA methylation of mouse epidermis and find that fresh and fresh-frozen epidermis can generally be used equally well. The choice of method will depend on the specific study's aims and the available resources in the laboratory.

Comprehensive Evaluation of The Infinium Human MethylationEPIC v2 BeadChip

Infinium Methylation BeadChips are widely used to profile DNA cytosine modifications in large cohort studies for reasons of cost-effectiveness, accurate quantification, and user-friendly data analysis in characterizing these canonical epigenetic marks. In this work, we conducted a comprehensive evaluation of the updated Infinium MethylationEPIC v2 BeadChip (EPICv2). Our evaluation revealed that EPICv2 offers significant improvements over its predecessors, including expanded enhancer coverage, applicability to diverse ancestry groups, support for low-input DNA down to one nanogram, coverage of existing epigenetic clocks, cell type deconvolution panels, and human trait associations, while maintaining accuracy and reproducibility. Using EPICv2, we were able to identify epigenome and sequence signatures in cell line models of DNMT and SETD2 loss and/or hypomorphism. Furthermore, we provided probe-wise evaluation and annotation to facilitate the use of new features on this array for studying the interplay between somatic mutations and epigenetic landscape in cancer genomics. In conclusion, EPICv2 provides researchers with a valuable tool for studying epigenetic modifications and their role in development and disease.

Developmental priming of cancer susceptibility

DNA mutations are necessary drivers of cancer, yet only a small subset of mutated cells go on to cause the disease. To date, the mechanisms that determine which rare subset of cells transform and initiate tumorigenesis remain unclear. Here, we take advantage of a unique model of intrinsic developmental heterogeneity (Trim28+/D9) and demonstrate that stochastic early life epigenetic variation can trigger distinct cancer-susceptibility 'states' in adulthood. We show that these developmentally primed states are characterized by differential methylation patterns at typically silenced heterochromatin, and that these epigenetic signatures are detectable as early as 10 days of age. The differentially methylated loci are enriched for genes with known oncogenic potential. These same genes are frequently mutated in human cancers, and their dysregulation correlates with poor prognosis. These results provide proof-of-concept that intrinsic developmental heterogeneity can prime individual, life-long cancer risk.

Region-based epigenetic clock design improves RRBS-based age prediction

Recent studies suggest that epigenetic rejuvenation can be achieved using drugs that mimic calorie restriction and techniques such as reprogramming-induced rejuvenation. To effectively test rejuvenation in vivo, mouse models are the safest alternative. However, we have found that the recent epigenetic clocks developed for mouse reduced-representation bisulphite sequencing (RRBS) data have significantly poor performance when applied to external datasets. We show that the sites captured and the coverage of key CpGs required for age prediction vary greatly between datasets, which likely contributes to the lack of transferability in RRBS clocks. To mitigate these coverage issues in RRBS-based age prediction, we present two novel design strategies that use average methylation over large regions rather than individual CpGs, whereby regions are defined by sliding windows (e.g. 5 kb), or density-based clustering of CpGs. We observe improved correlation and error in our regional blood clocks (RegBCs) compared to published individual-CpG-based techniques when applied to external datasets. The RegBCs are also more robust when applied to low coverage data and detect a negative age acceleration in mice undergoing calorie restriction. Our RegBCs offer a proof of principle that age prediction of RRBS datasets can be improved by accounting for multiple CpGs over a region, which negates the lack of read depth currently hindering individual-CpG-based approaches.

Genotoxic aldehyde stress prematurely ages hematopoietic stem cells in a p53-driven manner

Aged hematopoietic stem cells (HSCs) display diminished self-renewal and a myeloid differentiation bias. However, the drivers and mechanisms that underpin this fundamental switch are not understood. HSCs produce genotoxic formaldehyde that requires protection by the detoxification enzymes ALDH2 and ADH5 and the Fanconi anemia (FA) DNA repair pathway. We find that the HSCs in young Aldh2-/-Fancd2-/- mice harbor a transcriptomic signature equivalent to aged wild-type HSCs, along with increased epigenetic age, telomere attrition, and myeloid-biased differentiation quantified by single HSC transplantation. In addition, the p53 response is vigorously activated in Aldh2-/-Fancd2-/- HSCs, while p53 deletion rescued this aged HSC phenotype. To further define the origins of the myeloid differentiation bias, we use a GFP genetic reporter to find a striking enrichment of Vwf+ myeloid and megakaryocyte-lineage-biased HSCs. These results indicate that metabolism-derived formaldehyde-DNA damage stimulates the p53 response in HSCs to drive accelerated aging.

Variation in histone configurations correlates with gene expression across nine inbred strains of mice

The Diversity Outbred (DO) mice and their inbred founders are widely used models of human disease. However, although the genetic diversity of these mice has been well documented, their epigenetic diversity has not. Epigenetic modifications, such as histone modifications and DNA methylation, are important regulators of gene expression and, as such, are a critical mechanistic link between genotype and phenotype. Therefore, creating a map of epigenetic modifications in the DO mice and their founders is an important step toward understanding mechanisms of gene regulation and the link to disease in this widely used resource. To this end, we performed a strain survey of epigenetic modifications in hepatocytes of the DO founders. We surveyed four histone modifications (H3K4me1, H3K4me3, H3K27me3, and H3K27ac), as well as DNA methylation. We used ChromHMM to identify 14 chromatin states, each of which represents a distinct combination of the four histone modifications. We found that the epigenetic landscape is highly variable across the DO founders and is associated with variation in gene expression across strains. We found that epigenetic state imputed into a population of DO mice recapitulated the association with gene expression seen in the founders, suggesting that both histone modifications and DNA methylation are highly heritable mechanisms of gene expression regulation. We illustrate how DO gene expression can be aligned with inbred epigenetic states to identify putative cis-regulatory regions. Finally, we provide a data resource that documents strain-specific variation in the chromatin state and DNA methylation in hepatocytes across nine widely used strains of laboratory mice.

Epigenetic dosage identifies two major and functionally distinct β cell subtypes

The mechanisms that specify and stabilize cell subtypes remain poorly understood. Here, we identify two major subtypes of pancreatic β cells based on histone mark heterogeneity (βHI and βLO). βHI cells exhibit ∼4-fold higher levels of H3K27me3, distinct chromatin organization and compaction, and a specific transcriptional pattern. βHI and βLO cells also differ in size, morphology, cytosolic and nuclear ultrastructure, epigenomes, cell surface marker expression, and function, and can be FACS separated into CD24+ and CD24- fractions. Functionally, βHI cells have increased mitochondrial mass, activity, and insulin secretion in vivo and ex vivo. Partial loss of function indicates that H3K27me3 dosage regulates βHI/βLO ratio in vivo, suggesting that control of β cell subtype identity and ratio is at least partially uncoupled. Both subtypes are conserved in humans, with βHI cells enriched in humans with type 2 diabetes. Thus, epigenetic dosage is a novel regulator of cell subtype specification and identifies two functionally distinct β cell subtypes.

Contextual modifiers of healthspan, lifespan, and epigenome in mice under chronic social stress

Sustained life stress and low socioeconomic status are among the major causes of aging-related diseases and decreased life expectancy. Experimental rodent models can help to identify the underlying mechanisms, yet very few studies address the long-term consequences of social stress on aging. We conducted a randomized study involving more than 300 male mice of commonly used laboratory strains (C57BL/6J, CD1, and Sv129Ev) chosen for the spontaneous aggression gradient and stress-vulnerability. Mice were exposed to a lifelong chronic psychosocial stress protocol to model social gradients in aging and disease vulnerability. Low social rank, inferred based on a discretized aggression index, was found to negatively impact lifespan in our study population. However, social rank interacted with genetic background in that low-ranking C57BL/6J, high-ranking Sv129Ev, and middle-ranking CD1 mice had lower survival, respectively, implying a cost of maintaining a given social rank that varies across strains. Machine learning linear discriminant analysis identified baseline fat-free mass as the most important predictor of mouse genetic background and social rank in the present dataset. Finally, strain and social rank differences were significantly associated with epigenetic changes, most significantly in Sv129Ev mice and in high-ranking compared to lower ranking subjects. Overall, we identified genetic background and social rank as critical contextual modifiers of aging and lifespan in an ethologically relevant rodent model of social stress, thereby providing a preclinical experimental paradigm to study the impact of social determinants of health disparities and accelerated aging.

Epigenetic state determines the in vivo efficacy of STING agonist therapy

While STING-activating agents have shown limited efficacy in early-phase clinical trials, multiple lines of evidence suggest the importance of tumor cell-intrinsic STING function in mediating antitumor immune responses. Although STING signaling is impaired in human melanoma, its restoration through epigenetic reprogramming can augment its antigenicity and T cell recognition. In this study, we show that reversal of methylation silencing of STING in murine melanoma cell lines using a clinically available DNA methylation inhibitor can improve agonist-induced STING activation and type-I IFN induction, which, in tumor-bearing mice, can induce tumor regression through a CD8⁺ T cell-dependent immune response. These findings not only provide mechanistic insight into how STING signaling dysfunction in tumor cells can contribute to impaired responses to STING agonist therapy, but also suggest that pharmacological restoration of STING signaling through epigenetic reprogramming might improve the therapeutic efficacy of STING agonists.

ARF suppression by MYC but not MYCN confers increased malignancy of aggressive pediatric brain tumors

Medulloblastoma, the most common malignant pediatric brain tumor, often harbors MYC amplifications. Compared to high-grade gliomas, MYC-amplified medulloblastomas often show increased photoreceptor activity and arise in the presence of a functional ARF/p53 suppressor pathway. Here, we generate an immunocompetent transgenic mouse model with regulatable MYC that develop clonal tumors that molecularly resemble photoreceptor-positive Group 3 medulloblastoma. Compared to MYCN-expressing brain tumors driven from the same promoter, pronounced ARF silencing is present in our MYC-expressing model and in human medulloblastoma. While partial Arf suppression causes increased malignancy in MYCN-expressing tumors, complete Arf depletion promotes photoreceptor-negative high-grade glioma formation. Computational models and clinical data further identify drugs targeting MYC-driven tumors with a suppressed but functional ARF pathway. We show that the HSP90 inhibitor, Onalespib, significantly targets MYC-driven but not MYCN-driven tumors in an ARF-dependent manner. The treatment increases cell death in synergy with cisplatin and demonstrates potential for targeting MYC-driven medulloblastoma.

TET1 Catalytic Activity is Required for Reprogramming of Imprinting Control Regions and Patterning of Sperm-Specific Hypomethylated Regions

DNA methylation erasure is required for mammalian primordial germ cell reprogramming. TET enzymes iteratively oxidize 5-methylcytosine to generate 5-hyroxymethylcytosine (5hmC), 5-formylcytosine, and 5-carboxycytosine to facilitate active genome demethylation. Whether these bases are required to promote replication-coupled dilution or activate base excision repair during germline reprogramming remains unresolved due to the lack of genetic models that decouple TET activities. Here, we generated two mouse lines expressing catalytically inactive TET1 ( Tet1-HxD ) and TET1 that stalls oxidation at 5hmC ( Tet1-V ). Tet1 -/- , Tet1 V/V , and Tet1 HxD/HxD sperm methylomes show that TET1 V and TET1 HxD rescue most Tet1 -/- hypermethylated regions, demonstrating the importance of TET1’s extra-catalytic functions. Imprinted regions, in contrast, require iterative oxidation. We further reveal a broader class of hypermethylated regions in sperm of Tet1 mutant mice that are excluded from de novo methylation during male germline development and depend on TET oxidation for reprogramming. Our study underscores the link between TET1-mediated demethylation during reprogramming and sperm methylome patterning.

Beadchip technology to detect DNA methylation in mouse faithfully recapitulates whole-genome bisulfite sequencing

Aim: To facilitate wide-scale implementation of Illumina Mouse Methylation BeadChip (MMB) technology, array-based measurement of cytosine methylation was compared with the gold-standard assessment of DNA methylation by whole-genome bisulfite sequencing (WGBS). Methods: DNA methylation across two mouse strains (C57B6 and C3H) and both sexes was assessed using the MMB and compared with previously existing deep-coverage WGBS of mice of the same strain and sex. Results & conclusion: The findings demonstrated that 93.3-99.2% of sites had similar measurements of methylation across technologies and that differentially methylated cytosines and regions identified by each technology overlap and enrich for similar biological functions, suggesting that the MMB faithfully recapitulates the findings of WGBS.

Comparative epigenome analysis using Infinium DNA methylation BeadChips

The arrival of the Infinium DNA methylation BeadChips for mice and other nonhuman mammalian species has outpaced the development of the informatics that supports their use for epigenetics study in model organisms. Here, we present informatics infrastructure and methods to allow easy DNA methylation analysis on multiple species, including domesticated animals and inbred laboratory mice (in SeSAMe version 1.16.0+). First, we developed a data-driven analysis pipeline covering species inference, genome-specific data preprocessing and regression modeling. We targeted genomes of 310 species and 37 inbred mouse strains and showed that genome-specific preprocessing prevents artifacts and yields more accurate measurements than generic pipelines. Second, we uncovered the dynamics of the epigenome evolution in different genomic territories and tissue types through comparative analysis. We identified a catalog of inbred mouse strain-specific methylation differences, some of which are linked to the strains' immune, metabolic and neurological phenotypes. By streamlining DNA methylation array analysis for undesigned genomes, our methods extend epigenome research to broad species contexts.

Dynamic DNA methylation reveals novel cis-regulatory elements in mouse hematopoiesis

Differentiation of hematopoietic stem and progenitor cells to terminally differentiated immune cells is accompanied by large-scale remodeling of the DNA methylation landscape. Although significant insights into the molecular mechanisms of hematopoietic tissue regeneration were derived from mouse models, profiling of DNA methylation has been hampered by high cost or low resolution using available methods. The recent development of the Infinium Mouse Methylation BeadChip (MMBC) array facilitates methylation profiling of the mouse genome at a single CpG resolution at affordable cost. We extended the RnBeads package to provide a computational framework for the analysis of MMBC data. This framework was applied to a newly generated reference map of mouse hematopoiesis encompassing nine different cell types. Analysis of dynamically regulated CpG sites showed progressive and unidirectional DNA methylation changes from hematopoietic stem and progenitor cells to differentiated hematopoietic cells and allowed the identification of lineage- and cell type-specific DNA methylation programs. Comparison with previously published catalogs of cis-regulatory elements (CREs) revealed 12,856 novel putative CREs that were dynamically regulated by DNA methylation (mdCREs). These mdCREs were predominantly associated with patterns of cell type-specific DNA hypomethylation and could be identified as epigenetic control regions regulating the expression of key hematopoietic genes during differentiation. In summary, we established an analysis pipeline for MMBC data sets and provide a DNA methylation atlas of mouse hematopoiesis. This resource allowed us to identify novel putative CREs involved in hematopoiesis and will serve as a platform to study epigenetic regulation of normal and malignant hematopoiesis.