Cis-regulatory evolution of the recently expanded Ly49 gene family.

Comparative genomics has revealed the rapid expansion of multiple gene families involved in immunity. Members within each gene family often evolved distinct roles in immunity. However, less is known about the evolution of their epigenome and cis-regulation. Here we systematically profile the epigenome of the recently expanded murine Ly49 gene family that mainly encode either inhibitory or activating surface receptors on natural killer cells. We identify a set of cis-regulatory elements (CREs) for activating Ly49 genes. In addition, we show that in mice, inhibitory and activating Ly49 genes are regulated by two separate sets of proximal CREs, likely resulting from lineage-specific losses of CRE activity. Furthermore, we find that some Ly49 genes are cross-regulated by the CREs of other Ly49 genes, suggesting that the Ly49 family has begun to evolve a concerted cis-regulatory mechanism. Collectively, we demonstrate the different modes of cis-regulatory evolution for a rapidly expanding gene family.

Methamphetamine-induced region-specific transcriptomic and epigenetic changes in the brain of male rats.

Psychostimulant methamphetamine (METH) is neurotoxic to the brain and, therefore, its misuse leads to neurological and psychiatric disorders. The gene regulatory network (GRN) response to neurotoxic METH binge remains unclear in most brain regions. Here we examined the effects of binge METH on the GRN in the nucleus accumbens, dentate gyrus, Ammon's horn, and subventricular zone in male rats. At 24 h after METH, ~16% of genes displayed altered expression and over a quarter of previously open chromatin regions - parts of the genome where genes are typically active - showed shifts in their accessibility. Intriguingly, most changes were unique to each area studied, and independent regulation between transcriptome and chromatin accessibility was observed. Unexpectedly, METH differentially impacted gene activity and chromatin accessibility within the dentate gyrus and Ammon's horn. Around 70% of the affected chromatin-accessible regions in the rat brain have conserved DNA sequences in the human genome. These regions frequently act as enhancers, ramping up the activity of nearby genes, and contain mutations linked to various neurological conditions. By sketching out the gene regulatory networks associated with binge METH in specific brain regions, our study offers fresh insights into how METH can trigger profound, region-specific molecular shifts.

Non-small Cell Lung Cancer Epigenomes Exhibit Altered DNA Methylation in Smokers and Never-smokers.

Epigenetic alterations are widespread in cancer and can complement genetic alterations to influence cancer progression and treatment outcome. To determine the potential contribution of DNAmethylation alterations to tumor phenotype in non-small cell lung cancer (NSCLC) in both smoker and never-smoker patients, we performed genome-wide profiling of DNA methylation in 17 primary NSCLC tumors and 10 matched normal lung samples using the complementary assays, methylated DNA immunoprecipitation sequencing (MeDIP-seq) and methylation sensitive restriction enzyme sequencing (MRE-seq). We reported recurrent methylation changes in the promoters of several genes, many previously implicated in cancer, including FAM83A and SEPT9 (hypomethylation), as well as PCDH7, NKX2-1, and SOX17 (hypermethylation). Although many methylation changes between tumors and their paired normal samples were shared across patients, several were specific to a particular smoking status. For example, never-smokers displayed a greater proportion of hypomethylated differentially methylated regions (hypoDMRs) and a greater number of recurrently hypomethylated promoters, including those of ASPSCR1, TOP2A, DPP9, and USP39, all previously linked to cancer. Changes outside of promoters were also widespread and often recurrent, particularly methylation loss over repetitive elements, highly enriched for ERV1 subfamilies. Recurrent hypoDMRs were enriched for several transcription factor binding motifs, often for genes involved in signaling and cell proliferation. For example, 71% of recurrent promoter hypoDMRs contained a motif for NKX2-1. Finally, the majority of DMRs were located within an active chromatin state in tissues profiled using the Roadmap Epigenomics data, suggesting that methylation changes may contribute to altered regulatory programs through the adaptation of cell type-specific expression programs.

Embryonic vitamin D deficiency programs hematopoietic stem cells to induce type 2 diabetes.

Environmental factors may alter the fetal genome to cause metabolic diseases. It is unknown whether embryonic immune cell programming impacts the risk of type 2 diabetes in later life. We demonstrate that transplantation of fetal hematopoietic stem cells (HSCs) made vitamin D deficient in utero induce diabetes in vitamin D-sufficient mice. Vitamin D deficiency epigenetically suppresses Jarid2 expression and activates the Mef2/PGC1a pathway in HSCs, which persists in recipient bone marrow, resulting in adipose macrophage infiltration. These macrophages secrete miR106-5p, which promotes adipose insulin resistance by repressing PIK3 catalytic and regulatory subunits and down-regulating AKT signaling. Vitamin D-deficient monocytes from human cord blood have comparable Jarid2/Mef2/PGC1a expression changes and secrete miR-106b-5p, causing adipocyte insulin resistance. These findings suggest that vitamin D deficiency during development has epigenetic consequences impacting the systemic metabolic milieu.

Uropathogenic Escherichia coli infection-induced epithelial trained immunity impacts urinary tract disease outcome.

Previous urinary tract infections (UTIs) can predispose one to future infections; however, the underlying mechanisms affecting recurrence are poorly understood. We previously found that UTIs in mice cause differential bladder epithelial (urothelial) remodelling, depending on disease outcome, that impacts susceptibility to recurrent UTI. Here we compared urothelial stem cell (USC) lines isolated from mice with a history of either resolved or chronic uropathogenic Escherichia coli (UPEC) infection, elucidating evidence of molecular imprinting that involved epigenetic changes, including differences in chromatin accessibility, DNA methylation and histone modification. Epigenetic marks in USCs from chronically infected mice enhanced caspase-1-mediated cell death upon UPEC infection, promoting bacterial clearance. Increased Ptgs2os2 expression also occurred, potentially contributing to sustained cyclooxygenase-2 expression, bladder inflammation and mucosal wounding-responses associated with severe recurrent cystitis. Thus, UPEC infection acts as an epi-mutagen reprogramming the urothelial epigenome, leading to urothelial-intrinsic remodelling and training of the innate response to subsequent infection.

Capture Methylation-Sensitive Restriction Enzyme Sequencing (Capture MRE-Seq) for Methylation Analysis of Highly Degraded DNA Samples.

Understanding the impact of DNA methylation within different disease contexts often requires accurate assessment of these modifications in a genome-wide fashion. Frequently, patient-derived tissues stored in long-term hospital tissue banks have been preserved using formalin-fixation paraffin-embedding (FFPE). While these samples can comprise valuable resources for studying disease, the fixation process ultimately compromises the DNA's integrity and leads to degradation. Degraded DNA can complicate CpG methylome profiling using traditional techniques, particularly when performing methylation-sensitive restriction enzyme sequencing (MRE-seq), yielding high backgrounds and resulting in lowered library complexity. Here, we describe Capture MRE-seq, a new MRE-seq protocol tailored to preserving unmethylated CpG information when using samples with highly degraded DNA. The results using Capture MRE-seq correlate well (0.92) with traditional MRE-seq calls when profiling non-degraded samples, and can recover unmethylated regions in highly degraded samples when traditional MRE-seq fails, which we validate using bisulfite sequencing-based data (WGBS) as well as methylated DNA immunoprecipitation followed by sequencing (MeDIP-seq).

Pan-cancer analysis identifies tumor-specific antigens derived from transposable elements.

Cryptic promoters within transposable elements (TEs) can be transcriptionally reactivated in tumors to create new TE-chimeric transcripts, which can produce immunogenic antigens. We performed a comprehensive screen for these TE exaptation events in 33 TCGA tumor types, 30 GTEx adult tissues and 675 cancer cell lines, and identified 1,068 TE-exapted candidates with the potential to generate shared tumor-specific TE-chimeric antigens (TS-TEAs). Whole-lysate and HLA-pulldown mass spectrometry data confirmed that TS-TEAs are presented on the surface of cancer cells. In addition, we highlight tumor-specific membrane proteins transcribed from TE promoters that constitute aberrant epitopes on the extracellular surface of cancer cells. Altogether, we showcase the high pan-cancer prevalence of TS-TEAs and atypical membrane proteins that could potentially be therapeutically exploited and targeted.

Widespread contribution of transposable elements to the rewiring of mammalian 3D genomes.

Transposable elements (TEs) are major contributors of genetic material in mammalian genomes. These often include binding sites for architectural proteins, including the multifarious master protein, CTCF, which shapes the 3D genome by creating loops, domains, compartment borders, and RNA-DNA interactions. These play a role in the compact packaging of DNA and have the potential to facilitate regulatory function. In this study, we explore the widespread contribution of TEs to mammalian 3D genomes by quantifying the extent to which they give rise to loops and domain border differences across various cell types and species using several 3D genome mapping technologies. We show that specific families and subfamilies of TEs have contributed to lineage-specific 3D chromatin structures across mammalian species. In many cases, these loops may facilitate sustained interaction between distant cis-regulatory elements and target genes, and domains may segregate chromatin state to impact gene expression in a lineage-specific manner. An experimental validation of our analytical findings using CRISPR-Cas9 to delete a candidate TE resulted in disruption of species-specific 3D chromatin structure. Taken together, we comprehensively quantify and selectively validate our finding that TEs contribute to shaping 3D genome organization and may, in some cases, impact gene regulation during the course of mammalian evolution.

The transcription factor Bach2 negatively regulates murine natural killer cell maturation and function.

BTB domain And CNC Homolog 2 (Bach2) is a transcription repressor that actively participates in T and B lymphocyte development, but it is unknown if Bach2 is also involved in the development of innate immune cells, such as natural killer (NK) cells. Here, we followed the expression of Bach2 during murine NK cell development, finding that it peaked in immature CD27+CD11b+ cells and decreased upon further maturation. Bach2 showed an organ and tissue-specific expression pattern in NK cells. Bach2 expression positively correlated with the expression of transcription factor TCF1 and negatively correlated with genes encoding NK effector molecules and those involved in the cell cycle. Lack of Bach2 expression caused changes in chromatin accessibility of corresponding genes. In the end, Bach2 deficiency resulted in increased proportions of terminally differentiated NK cells with increased production of granzymes and cytokines. NK cell-mediated control of tumor metastasis was also augmented in the absence of Bach2. Therefore, Bach2 is a key checkpoint protein regulating NK terminal maturation.

Developmental Pathways Are Epigenetically Reprogrammed during Lung Cancer Brain Metastasis.

Non-small cell lung cancer (NSCLC) is one of the most commonly diagnosed and deadliest cancers worldwide, with roughly half of all patients initially presenting with both primary and metastatic disease. While the major events in the metastatic cascade have been identified, a mechanistic understanding of how NSCLC routinely and successfully colonizes the brain is largely unknown. Recent studies have begun demonstrating the role of epigenetic misregulation during tumorigenesis and metastasis, including widespread changes in DNA methylation and histone modifications. To better understand the role of altered DNA methylation in NSCLC metastasis to the brain, we measured DNA methylation during disease progression for 12 patients, globally profiling the methylation status of normal lung, primary lung tumor, and brain metastasis samples. The variation in methylation was similar during metastatic spread and primary tumorigenesis but less coordinated across genomic features during metastasis. The greatest recurrent changes during metastatic progression were methylation gains in DNA methylation valleys (DMV) harboring the constitutive heterochromatin mark H3K9me3 as well as bivalent marks H3K27me3 and H3K4me1. In a lymph node-derived cancer cell line, EZH2 binding within DMVs was lost, accompanied by an increase in DNA methylation, exemplifying epigenetic switching. The vast majority of the differentially methylated region-associated DMVs harbored developmental genes, suggesting that altered epigenetic regulation of developmentally important genes may confer a selective advantage during metastatic progression. The characterization of epigenetic changes during NSCLC brain metastasis identified recurrent methylation patterns that may be prognostic biomarkers and contributors to disease progression. SIGNIFICANCE: Altered DNA methylation in lung cancer brain metastases corresponds with loss of EZH2 occupancy at developmental genes, which could promote stem-like phenotypes permissive of dissemination and survival in different microenvironments.

Notch1/Hes1­PTEN/AKT/IL­17A feedback loop regulates Th17 cell differentiation in mouse psoriasis­like skin inflammation.

IL­17A, the effector cytokine of T helper (Th) 17 cells, plays a crucial role in the pathogenesis of psoriasis. The Notch1 and PI3K/AKT signaling pathways are implicated in Th17 cell differentiation and IL­17A production. The present study aimed to evaluate the regulatory effect of the Notch1/hairy and enhancer of split 1 (Hes1)­PTEN/AKT/IL­17A feedback loop on Th17 cell differentiation via the PI3K/AKT inhibitor LY294002 in a mouse model of psoriasis. Mice were randomly divided into 3 groups: a control group, a model group [5% imiquimod (IMQ)­induced group] and an intervention group (5% IMQ­induced plus LY294002­treated group). Skin structural characteristics were recorded and evaluated by hematoxylin and eosin staining. The weights of the spleens and inguinal lymph nodes were measured. Th17 cell percentage, as well as the mRNA and protein expression levels of Notch1, Notch1 intracellular domain (NICD1), Hes1, PTEN, AKT, phosphorylated (p)­AKT, mTOR complex 1 (mTORC1), p­mTORC1, S6 kinase (S6K)1, S6K2 and IL­17A were detected in skin samples of the three experimental groups. Additionally, splenic mononuclear cells from model mice were treated by 10 and 50 µM LY294002 to further evaluate its regulatory effect on Notch1/Hes1­PTEN/AKT/IL­17A feedback loop. Increased Th17 cell percentage, increased expression of Notch1, NICD1, Hes1, AKT, p­AKT, mTORC1, p­mTORC1, S6K1, S6K2 and IL­17A, and decreased PTEN levels were observed in model mice alongside marked psoriasis­like skin inflammation, splenomegaly and lymphadenopathy. LY294002 treatment significantly alleviated the severity of psoriasis­like skin inflammation in the intervention mice, attenuated the degree of epidermal hyperplasia and dermal inflammatory cell infiltration, and mitigated splenomegaly and lymphadenopathy. In addition, LY294002 treatment reversed the increased Th17 cell percentage, as well as the increased expression of Notch1, NICD1, Hes1, AKT, p­AKT, mTORC1, p­mTORC1, S6K1, S6K2 and IL­17A, and the decreased expression of PTEN. In vitro study from 5% IMQ­induced mouse splenic mononuclear cells presented that high dose of LY294002 exerted more obviously regulatory effect on Notch1/Hes1­PTEN/AKT/IL­17A feedback loop. The current findings suggested that the Notch1/Hes1­PTEN/AKT/IL­17A feedback loop regulates Th17 cell differentiation within the disease environment of psoriasis. Blocking the Notch1/Hes1­PTEN/AKT/IL­17A feedback loop may thus be a potential therapeutic method for management of psoriatic inflammation.

A genome-wide CRISPR-Cas9 knockout screen identifies essential and growth-restricting genes in human trophoblast stem cells.

The recent derivation of human trophoblast stem cells (hTSCs) provides a scalable in vitro model system of human placental development, but the molecular regulators of hTSC identity have not been systematically explored thus far. Here, we utilize a genome-wide CRISPR-Cas9 knockout screen to comprehensively identify essential and growth-restricting genes in hTSCs. By cross-referencing our data to those from similar genetic screens performed in other cell types, as well as gene expression data from early human embryos, we define hTSC-specific and -enriched regulators. These include both well-established and previously uncharacterized trophoblast regulators, such as ARID3A, GATA2, and TEAD1 (essential), and GCM1, PTPN14, and TET2 (growth-restricting). Integrated analysis of chromatin accessibility, gene expression, and genome-wide location data reveals that the transcription factor TEAD1 regulates the expression of many trophoblast regulators in hTSCs. In the absence of TEAD1, hTSCs fail to complete faithful differentiation into extravillous trophoblast (EVT) cells and instead show a bias towards syncytiotrophoblast (STB) differentiation, thus indicating that this transcription factor safeguards the bipotent lineage potential of hTSCs. Overall, our study provides a valuable resource for dissecting the molecular regulation of human placental development and diseases.

Whole-genome profiling of DNA methylation and hydroxymethylation identifies distinct regulatory programs among innate lymphocytes.

Innate lymphocytes encompass a diverse array of phenotypic identities with specialized functions. DNA methylation and hydroxymethylation are essential for epigenetic fidelity and fate commitment. The landscapes of these modifications are unknown in innate lymphocytes. Here, we characterized the whole-genome distribution of methyl-CpG and 5-hydroxymethylcytosine (5hmC) in mouse innate lymphoid cell 3 (ILC3), ILC2 and natural killer (NK) cells. We identified differentially methylated regions (DMRs) and differentially hydroxymethylated regions (DHMRs) between ILC and NK cell subsets and correlated them with transcriptional signatures. We associated lineage-determining transcription factors (LDTFs) with demethylation and demonstrated unique patterns of DNA methylation/hydroxymethylation in relationship to open chromatin regions (OCRs), histone modifications and TF-binding sites. We further identified an association between hydroxymethylation and NK cell superenhancers (SEs). Using mice lacking the DNA hydroxymethylase TET2, we showed the requirement for TET2 in optimal production of hallmark cytokines by ILC3s and interleukin-17A (IL-17A) by inflammatory ILC2s. These findings provide a powerful resource for studying innate lymphocyte epigenetic regulation and decode the regulatory logic governing their identity.

Epigenetic dynamics shaping melanophore and iridophore cell fate in zebrafish.

BACKGROUND: Zebrafish pigment cell differentiation provides an attractive model for studying cell fate progression as a neural crest progenitor engenders diverse cell types, including two morphologically distinct pigment cells: black melanophores and reflective iridophores. Nontrivial classical genetic and transcriptomic approaches have revealed essential molecular mechanisms and gene regulatory circuits that drive neural crest-derived cell fate decisions. However, how the epigenetic landscape contributes to pigment cell differentiation, especially in the context of iridophore cell fate, is poorly understood. RESULTS: We chart the global changes in the epigenetic landscape, including DNA methylation and chromatin accessibility, during neural crest differentiation into melanophores and iridophores to identify epigenetic determinants shaping cell type-specific gene expression. Motif enrichment in the epigenetically dynamic regions reveals putative transcription factors that might be responsible for driving pigment cell identity. Through this effort, in the relatively uncharacterized iridophores, we validate alx4a as a necessary and sufficient transcription factor for iridophore differentiation and present evidence on alx4a's potential regulatory role in guanine synthesis pathway. CONCLUSIONS: Pigment cell fate is marked by substantial DNA demethylation events coupled with dynamic chromatin accessibility to potentiate gene regulation through cis-regulatory control. Here, we provide a multi-omic resource for neural crest differentiation into melanophores and iridophores. This work led to the discovery and validation of iridophore-specific alx4a transcription factor.

OCT4 cooperates with distinct ATP-dependent chromatin remodelers in naïve and primed pluripotent states in human.

Understanding the molecular underpinnings of pluripotency is a prerequisite for optimal maintenance and application of embryonic stem cells (ESCs). While the protein-protein interactions of core pluripotency factors have been identified in mouse ESCs, their interactome in human ESCs (hESCs) has not to date been explored. Here we mapped the OCT4 interactomes in naïve and primed hESCs, revealing extensive connections to mammalian ATP-dependent nucleosome remodeling complexes. In naïve hESCs, OCT4 is associated with both BRG1 and BRM, the two paralog ATPases of the BAF complex. Genome-wide location analyses and genetic studies reveal that these two enzymes cooperate in a functionally redundant manner in the transcriptional regulation of blastocyst-specific genes. In contrast, in primed hESCs, OCT4 cooperates with BRG1 and SOX2 to promote chromatin accessibility at ectodermal genes. This work reveals how a common transcription factor utilizes differential BAF complexes to control distinct transcriptional programs in naïve and primed hESCs.

Epigenetic Therapies in Ovarian Cancer Alter Repetitive Element Expression in a TP53-Dependent Manner.

Epithelial ovarian carcinomas are particularly deadly due to intratumoral heterogeneity, resistance to standard-of-care therapies, and poor response to alternative treatments such as immunotherapy. Targeting the ovarian carcinoma epigenome with DNA methyltransferase inhibitors (DNMTi) or histone deacetylase inhibitors (HDACi) increases immune signaling and recruits CD8+ T cells and natural killer cells to fight ovarian carcinoma in murine models. This increased immune activity is caused by increased transcription of repetitive elements (RE) that form double-stranded RNA (dsRNA) and trigger an IFN response. To understand which REs are affected by epigenetic therapies in ovarian carcinoma, we assessed the effect of DNMTi and HDACi on ovarian carcinoma cell lines and patient samples. Subfamily-level (TEtranscripts) and individual locus-level (Telescope) analysis of REs showed that DNMTi treatment upregulated more REs than HDACi treatment. Upregulated REs were predominantly LTR and SINE subfamilies, and SINEs exhibited the greatest loss of DNA methylation upon DNMTi treatment. Cell lines with TP53 mutations exhibited significantly fewer upregulated REs with epigenetic therapy than wild-type TP53 cell lines. This observation was validated using isogenic cell lines; the TP53-mutant cell line had significantly higher baseline expression of REs but upregulated fewer upon epigenetic treatment. In addition, p53 activation increased expression of REs in wild-type but not mutant cell lines. These data give a comprehensive, genome-wide picture of RE chromatin and transcription-related changes in ovarian carcinoma after epigenetic treatment and implicate p53 in RE transcriptional regulation. SIGNIFICANCE: This study identifies the repetitive element targets of epigenetic therapies in ovarian carcinoma and indicates a role for p53 in this process.

Common DNA methylation dynamics in endometriod adenocarcinoma and glioblastoma suggest universal epigenomic alterations in tumorigenesis.

Trends in altered DNA methylation have been defined across human cancers, revealing global loss of methylation (hypomethylation) and focal gain of methylation (hypermethylation) as frequent cancer hallmarks. Although many cancers share these trends, little is known about the specific differences in DNA methylation changes across cancer types, particularly outside of promoters. Here, we present a comprehensive comparison of DNA methylation changes between two distinct cancers, endometrioid adenocarcinoma (EAC) and glioblastoma multiforme (GBM), to elucidate common rules of methylation dysregulation and changes unique to cancers derived from specific cells. Both cancers exhibit significant changes in methylation over regulatory elements. Notably, hypermethylated enhancers within EAC samples contain several transcription factor binding site clusters with enriched disease ontology terms highlighting uterine function, while hypermethylated enhancers in GBM are found to overlap active enhancer marks in adult brain. These findings suggest that loss of original cellular identity may be a shared step in tumorigenesis.

Requisite Chromatin Remodeling for Myeloid and Erythroid Lineage Differentiation from Erythromyeloid Progenitors.

The mammalian SWitch/Sucrose Non-Fermentable (SWI/SNF) chromatin-remodeling BAF (BRG1/BRM-associated factor) complex plays an essential role in developmental and pathological processes. We show that the deletion of Baf155, which encodes a subunit of the BAF complex, in the Tie2(+) lineage (Baf155 (CKO) leads to defects in yolk sac myeloid and definitive erythroid (EryD) lineage differentiation from erythromyeloid progenitors (EMPs). The chromatin of myeloid gene loci in Baf155 CKO EMPs is mostly inaccessible and enriched mainly by the ETS binding motif. BAF155 interacts with PU.1 and is recruited to PU.1 target gene loci together with p300 and KDM6a. Treatment of Baf155 CKO embryos with GSK126, an H3K27me2/3 methyltransferase EZH2 inhibitor, rescues myeloid lineage gene expression. This study uncovers indispensable BAF-mediated chromatin remodeling of myeloid gene loci at the EMP stage. Future studies exploiting epigenetics in the generation and application of EMP derivatives for tissue repair, regeneration, and disease are warranted.

A map of cis-regulatory elements and 3D genome structures in zebrafish.

The zebrafish (Danio rerio) has been widely used in the study of human disease and development, and about 70% of the protein-coding genes are conserved between the two species1. However, studies in zebrafish remain constrained by the sparse annotation of functional control elements in the zebrafish genome. Here we performed RNA sequencing, assay for transposase-accessible chromatin using sequencing (ATAC-seq), chromatin immunoprecipitation with sequencing, whole-genome bisulfite sequencing, and chromosome conformation capture (Hi-C) experiments in up to eleven adult and two embryonic tissues to generate a comprehensive map of transcriptomes, cis-regulatory elements, heterochromatin, methylomes and 3D genome organization in the zebrafish Tübingen reference strain. A comparison of zebrafish, human and mouse regulatory elements enabled the identification of both evolutionarily conserved and species-specific regulatory sequences and networks. We observed enrichment of evolutionary breakpoints at topologically associating domain boundaries, which were correlated with strong histone H3 lysine 4 trimethylation (H3K4me3) and CCCTC-binding factor (CTCF) signals. We performed single-cell ATAC-seq in zebrafish brain, which delineated 25 different clusters of cell types. By combining long-read DNA sequencing and Hi-C, we assembled the sex-determining chromosome 4 de novo. Overall, our work provides an additional epigenomic anchor for the functional annotation of vertebrate genomes and the study of evolutionarily conserved elements of 3D genome organization.

DeepH&M: Estimating single-CpG hydroxymethylation and methylation levels from enrichment and restriction enzyme sequencing methods.

Increased appreciation of 5-hydroxymethylcytosine (5hmC) as a stable epigenetic mark, which defines cell identity and disease progress, has engendered a need for cost-effective, but high-resolution, 5hmC mapping technology. Current enrichment-based technologies provide cheap but low-resolution and relative enrichment of 5hmC levels, while single-base resolution methods can be prohibitively expensive to scale up to large experiments. To address this problem, we developed a deep learning-based method, "DeepH&M," which integrates enrichment and restriction enzyme sequencing methods to simultaneously estimate absolute hydroxymethylation and methylation levels at single-CpG resolution. Using 7-week-old mouse cerebellum data for training the DeepH&M model, we demonstrated that the 5hmC and 5mC levels predicted by DeepH&M were in high concordance with whole-genome bisulfite-based approaches. The DeepH&M model can be applied to 7-week-old frontal cortex and 79-week-old cerebellum, revealing the robust generalizability of this method to other tissues from various biological time points.