UHRF2 commissions the completion of DNA demethylation through allosteric activation by 5hmC and K33-linked ubiquitination of XRCC1.

The transition of oxidized 5-methylcytosine (5mC) intermediates into the base excision repair (BER) pipeline to complete DNA demethylation remains enigmatic. We report here that UHRF2, the only paralog of UHRF1 in mammals that fails to rescue Uhrf1-/- phenotype, is physically and functionally associated with BER complex. We show that UHRF2 is allosterically activated by 5-hydroxymethylcytosine (5hmC) and acts as a ubiquitin E3 ligase to catalyze K33-linked polyubiquitination of XRCC1. This nonproteolytic action stimulates XRCC1's interaction with the ubiquitin binding domain-bearing RAD23B, leading to the incorporation of TDG into BER complex. Integrative epigenomic analysis in mouse embryonic stem cells reveals that Uhrf2-fostered TDG-RAD23B-BER complex is functionally linked to the completion of DNA demethylation at active promoters and that Uhrf2 ablation impedes DNA demethylation on latent enhancers that undergo poised-to-active transition during neuronal commitment. Together, these observations highlight an essentiality of 5hmC-switched UHRF2 E3 ligase activity in commissioning the accomplishment of active DNA demethylation.

Functionally distinct roles for TET-oxidized 5-methylcytosine bases in somatic reprogramming to pluripotency.

Active DNA demethylation via ten-eleven translocation (TET) family enzymes is essential for epigenetic reprogramming in cell state transitions. TET enzymes catalyze up to three successive oxidations of 5-methylcytosine (5mC), generating 5-hydroxymethylcytosine (5hmC), 5-formylcytosine (5fC), or 5-carboxycytosine (5caC). Although these bases are known to contribute to distinct demethylation pathways, the lack of tools to uncouple these sequential oxidative events has constrained our mechanistic understanding of the role of TETs in chromatin reprogramming. Here, we describe the first application of biochemically engineered TET mutants that unlink 5mC oxidation steps, examining their effects on somatic cell reprogramming. We show that only TET enzymes proficient for oxidation to 5fC/5caC can rescue the reprogramming potential of Tet2-deficient mouse embryonic fibroblasts. This effect correlated with rapid DNA demethylation at reprogramming enhancers and increased chromatin accessibility later in reprogramming. These experiments demonstrate that DNA demethylation through 5fC/5caC has roles distinct from 5hmC in somatic reprogramming to pluripotency.

The DNA Methylcytosine Dioxygenase Tet2 Sustains Immunosuppressive Function of Tumor-Infiltrating Myeloid Cells to Promote Melanoma Progression.

Ten-Eleven-Translocation-2 (Tet2) is a DNA methylcytosine dioxygenase that functions as a tumor suppressor in hematopoietic malignancies. We examined the role of Tet2 in tumor-tissue myeloid cells and found that Tet2 sustains the immunosuppressive function of these cells. We found that Tet2 expression is increased in intratumoral myeloid cells both in mouse models of melanoma and in melanoma patients and that this increased expression is dependent on an IL-1R-MyD88 pathway. Ablation of Tet2 in myeloid cells suppressed melanoma growth in vivo and shifted the immunosuppressive gene expression program in tumor-associated macrophages to a proinflammatory one, with a concomitant reduction of the immunosuppressive function. This resulted in increased numbers of effector T cells in the tumor, and T cell depletion abolished the reduced tumor growth observed upon myeloid-specific deletion of Tet2. Our findings reveal a non-cell-intrinsic, tumor-promoting function for Tet2 and suggest that Tet2 may present a therapeutic target for the treatment of non-hematologic malignancies.

Single-cell bisulfite-free 5mC and 5hmC sequencing with high sensitivity and scalability.

Existing single-cell bisulfite-based DNA methylation analysis is limited by low DNA recovery, and the measurement of 5hmC at single-base resolution remains challenging. Here, we present a bisulfite-free single-cell whole-genome 5mC and 5hmC profiling technique, named Cabernet, which can characterize 5mC and 5hmC at single-base resolution with high genomic coverage. Cabernet utilizes Tn5 transposome for DNA fragmentation, which enables the discrimination between different alleles for measuring hemi-methylation status. Using Cabernet, we revealed the 5mC, hemi-5mC and 5hmC dynamics during early mouse embryo development, uncovering genomic regions exclusively governed by active or passive demethylation. We show that hemi-methylation status can be used to distinguish between pre- and post-replication cells, enabling more efficient cell grouping when integrated with 5mC profiles. The property of Tn5 naturally enables Cabernet to achieve high-throughput single-cell methylome profiling, where we probed mouse cortical neurons and embryonic day 7.5 (E7.5) embryos, and constructed the library for thousands of single cells at high efficiency, demonstrating its potential for analyzing complex tissues at substantially low cost. Together, we present a way of high-throughput methylome and hydroxymethylome detection at single-cell resolution, enabling efficient analysis of the epigenetic status of biological systems with complicated nature such as neurons and cancer cells.

Epigenetic marks or not? The discovery of novel DNA modifications in eukaryotes.

DNA modifications add another layer of complexity to the eukaryotic genome to regulate gene expression, playing critical roles as epigenetic marks. In eukaryotes, the study of DNA epigenetic modifications has been confined to 5mC and its derivatives for decades. However, rapid developing approaches have witnessed the expansion of DNA modification reservoirs during the past several years, including the identification of 6mA, 5gmC, 4mC, and 4acC in diverse organisms. However, whether these DNA modifications function as epigenetic marks requires careful consideration. In this review, we try to present a panorama of all the DNA epigenetic modifications in eukaryotes, emphasizing recent breakthroughs in the identification of novel DNA modifications. The characterization of their roles in transcriptional regulation as potential epigenetic marks is summarized. More importantly, the pathways for generating or eliminating these DNA modifications, as well as the proteins involved are comprehensively dissected. Furthermore, we briefly discuss the potential challenges and perspectives, which should be taken into account while investigating novel DNA modifications.

UHRF2 regulates cell cycle, epigenetics and gene expression to control the timing of retinal progenitor and ganglion cell differentiation.

Ubiquitin-like, containing PHD and RING finger domains 2 (UHRF2) regulates cell cycle and binds 5-hydroxymethylcytosine (5hmC) to promote completion of DNA demethylation. Uhrf2-/- mice are without gross phenotypic defects; however, the cell cycle and epigenetic regulatory functions of Uhrf2 during retinal tissue development are unclear. Retinal progenitor cells (RPCs) produce all retinal neurons and Müller glia in a predictable sequence controlled by the complex interplay between extrinsic signaling, cell cycle, epigenetic changes and cell-specific transcription factor activation. In this study, we find that UHRF2 accumulates in RPCs, and its conditional deletion from mouse RPCs reduced 5hmC, altered gene expressions and disrupted retinal cell proliferation and differentiation. Retinal ganglion cells were overproduced in Uhrf2-deficient retinae at the expense of VSX2+ RPCs. Most other cell types were transiently delayed in differentiation. Expression of each member of the Tet3/Uhrf2/Tdg active demethylation pathway was reduced in Uhrf2-deficient retinae, consistent with locally reduced 5hmC in their gene bodies. This study highlights a novel role of UHRF2 in controlling the transition from RPCs to differentiated cell by regulating cell cycle, epigenetic and gene expression decisions.

DNA-MP: a generalized DNA modifications predictor for multiple species based on powerful sequence encoding method.

Accurate prediction of deoxyribonucleic acid (DNA) modifications is essential to explore and discern the process of cell differentiation, gene expression and epigenetic regulation. Several computational approaches have been proposed for particular type-specific DNA modification prediction. Two recent generalized computational predictors are capable of detecting three different types of DNA modifications; however, type-specific and generalized modifications predictors produce limited performance across multiple species mainly due to the use of ineffective sequence encoding methods. The paper in hand presents a generalized computational approach "DNA-MP" that is competent to more precisely predict three different DNA modifications across multiple species. Proposed DNA-MP approach makes use of a powerful encoding method "position specific nucleotides occurrence based 117 on modification and non-modification class densities normalized difference" (POCD-ND) to generate the statistical representations of DNA sequences and a deep forest classifier for modifications prediction. POCD-ND encoder generates statistical representations by extracting position specific distributional information of nucleotides in the DNA sequences. We perform a comprehensive intrinsic and extrinsic evaluation of the proposed encoder and compare its performance with 32 most widely used encoding methods on $17$ benchmark DNA modifications prediction datasets of $12$ different species using $10$ different machine learning classifiers. Overall, with all classifiers, the proposed POCD-ND encoder outperforms existing $32$ different encoders. Furthermore, combinedly over 5-fold cross validation benchmark datasets and independent test sets, proposed DNA-MP predictor outperforms state-of-the-art type-specific and generalized modifications predictors by an average accuracy of 7% across 4mc datasets, 1.35% across 5hmc datasets and 10% for 6ma datasets. To facilitate the scientific community, the DNA-MP web application is available at https://sds_genetic_analysis.opendfki.de/DNA_Modifications/.

Acetylation Enhances TET2 Function in Protecting against Abnormal DNA Methylation during Oxidative Stress.

TET proteins, by converting 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), are hypothesized, but not directly shown, to protect promoter CpG islands (CGIs) against abnormal DNA methylation (DNAm) in cancer. We define such a protective role linked to DNA damage from oxidative stress (OS) known to induce this abnormality. TET2 removes aberrant DNAm during OS through interacting with DNA methyltransferases (DNMTs) in a "Yin-Yang" complex targeted to chromatin and enhanced by p300 mediated TET2 acetylation. Abnormal gains of DNAm and 5hmC occur simultaneously in OS, and knocking down TET2 dynamically alters this balance by enhancing 5mC and reducing 5hmC. TET2 reduction results in hypermethylation of promoter CGIs and enhancers in loci largely overlapping with those induced by OS. Thus, TET2 indeed may protect against abnormal, cancer DNAm in a manner linked to DNA damage.

Dynamic changes in RNA m6A and 5 hmC influence gene expression programs during macrophage differentiation and polarisation.

RNA modifications are essential for the establishment of cellular identity. Although increasing evidence indicates that RNA modifications regulate the innate immune response, their role in monocyte-to-macrophage differentiation and polarisation is unclear. While m6A has been widely studied, other RNA modifications, including 5 hmC, remain poorly characterised. We profiled m6A and 5 hmC epitranscriptomes, transcriptomes, translatomes and proteomes of monocytes and macrophages at rest and pro- and anti-inflammatory states. Transcriptome-wide mapping of m6A and 5 hmC reveals enrichment of m6A and/or 5 hmC on specific categories of transcripts essential for macrophage differentiation. Our analyses indicate that m6A and 5 hmC modifications are present in transcripts with critical functions in pro- and anti-inflammatory macrophages. Notably, we also discover the co-occurrence of m6A and 5 hmC on alternatively-spliced isoforms and/or opposing ends of the untranslated regions (UTR) of mRNAs with key roles in macrophage biology. In specific examples, RNA 5 hmC controls the decay of transcripts independently of m6A. This study provides (i) a comprehensive dataset to interrogate the role of RNA modifications in a plastic system (ii) a resource for exploring different layers of gene expression regulation in the context of human monocyte-to-macrophage differentiation and polarisation, (iii) new insights into RNA modifications as central regulators of effector cells in innate immunity.

Deciphering the multifaceted roles of TET proteins in T-cell lineage specification and malignant transformation.

TET proteins are DNA demethylases that can oxidize 5-methylcytosine (5mC) to generate 5-hydroxymethylcytosine (5hmC) and other oxidized mC bases (oxi-mCs). Importantly, TET proteins govern cell fate decisions during development of various cell types by activating a cell-specific gene expression program. In this review, we focus on the role of TET proteins in T-cell lineage specification. We explore the multifaceted roles of TET proteins in regulating gene expression in the contexts of T-cell development, lineage specification, function, and disease. Finally, we discuss the future directions and experimental strategies required to decipher the precise mechanisms employed by TET proteins to fine-tune gene expression and safeguard cell identity.

Targeting SIRT4/TET2 Signaling Alleviates Human Keratinocyte Senescence by Reducing 5-hydroxymethylcytosine Loss.

Skin aging is characterized by wrinkle formation and increased frailty and laxity, leading to the risk of age-related skin diseases. Keratinocyte is an important component of the epidermis in skin structure, and keratinocyte senescence has been identified as a pivotal factor in skin aging development. Because epigenetic pathways play a vital role in the regulation of skin aging, we evaluated human skin samples for DNA hydroxymethylation (5-hydroxymethylcytosine; 5-hmC) and SIRT4 expressions. Results found that both 5-hmC and SIRT4 showed a significant decrease in aged human skin samples. To test the results in vitro, human keratinocytes were cultured in H2O2, which modulates skin aging in vivo. However, H2O2-induced keratinocytes showed senescence-associated protein expression and significant downregulation of 5-hmC and SIRT4 expressions. Moreover, 5-hmC-converting enzymes ten eleven translocation 2 (TET2) showed a decrease and enhanced TET2 acetylation level in H2O2-induced keratinocytes. However, the overexpression of SIRT4 in keratinocytes alleviates the senescence phenotype, such as senescence-associated protein expression, decreases the TET2 acetylation, but increases TET2 and 5-hmC expressions. Our results provide a novel relevant mechanism whereby the epigenetic regulation of keratinocytes in skin aging may be correlated with SIRT4 expression and TET2 acetylation in 5-hmC alteration. Our study may provide a potential strategy for antiskin aging, which targets the SIRT4/TET2 axis involving epigenetic modification in keratinocyte senescence.

Region-specific DNA hydroxymethylation along the malignant progression of IDH-mutant gliomas.

The majority of low-grade isocitrate dehydrogenase-mutant (IDHmt) gliomas undergo malignant progression (MP), but their underlying mechanism remains unclear. IDHmt gliomas exhibit global DNA methylation, and our previous report suggested that MP could be partly attributed to passive demethylation caused by accelerated cell cycles. However, during MP, there is also active demethylation mediated by ten-eleven translocation, such as DNA hydroxymethylation. Hydroxymethylation is reported to potentially contribute to gene expression regulation, but its role in MP remains under investigation. Therefore, we conducted a comprehensive analysis of hydroxymethylation during MP of IDHmt astrocytoma. Five primary/malignantly progressed IDHmt astrocytoma pairs were analyzed with oxidative bisulfite and the Infinium EPIC methylation array, detecting 5-hydroxymethyl cytosine at over 850,000 locations for region-specific hydroxymethylation assessment. Notably, we observed significant sharing of hydroxymethylated genomic regions during MP across the samples. Hydroxymethylated CpGs were enriched in open sea and intergenic regions (p < 0.001), and genes undergoing hydroxymethylation were significantly associated with cancer-related signaling pathways. RNA sequencing data integration identified 91 genes with significant positive/negative hydroxymethylation-expression correlations. Functional analysis suggested that positively correlated genes are involved in cell-cycle promotion, while negatively correlated ones are associated with antineoplastic functions. Analyses of The Cancer Genome Atlas clinical data on glioma were in line with these findings. Motif-enrichment analysis suggested the potential involvement of the transcription factor KLF4 in hydroxymethylation-based gene regulation. Our findings shed light on the significance of region-specific DNA hydroxymethylation in glioma MP and suggest its potential role in cancer-related gene expression and IDHmt glioma malignancy.

TET2 regulates extranodal NK/T cell lymphoma progression through regulation of DNA methylation.

The biological role of Ten-11 translocation 2 (TET2) and the conversion of 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC) in the development of extra-nodal natural killer/T-cell lymphoma (ENKTL) remains unclear. The level of 5mC and 5hmC was detected in 112 cases of ENKTL tissue specimens by immunohistochemical (IHC) staining. Subsequently, TET2 knockdown and the overexpression cell models were constructed in ENKTL cell lines. Biochemical analyses were used to assess proliferation, apoptosis, cell cycle and monoclonal formation in cells treated or untreated with L-Ascorbic acid sodium salt (LAASS). Dot-Blots were used to detect levels of genome 5mC and 5hmC. Additionally, the ILLUMINA 850k methylation chip was used to analyze the changes of TET2 regulatory genes. RNA-Seq was used to profile differentially expressed genes regulated by TET2. The global level of 5hmC was significantly decreased, while 5mC was highly expressed in ENKTL tissue. TET2 protein expression was negatively correlated with the ratio of 5mC/5hmC (p < 0.0001). The 5mC/5hmC status were related to the site of disease, clinical stage, PINK score and Ki-67 index, as well as the 5-year OS. TET2 knockdown prolonged the DNA synthesis period, increased the cloning ability of tumor cells, increased the level of 5mC and decreased the level of 5hmC in ENKTL cells. While overexpression of TET2 presented the opposite effect. Furthermore, treatment of ENKTL cells with LAASS significantly induced ENKTL cell apoptosis. These results suggest that TET2 plays an important role in ENKTL development via regulation of 5mC and 5hmC and may serve as a novel therapeutic target for ENKTL.

Evaluation of tubulin ß-3 and 5 hydroxy-methyl cytosine as diagnostic and prognostic markers in malignant melanoma.

Tubulin ß-3 staining pattern and staining intensity of 5-hydroxymethyl cytosine (5-hmC) are potential diagnostic and prognostic markers in melanocytic lesions that need further evaluation. Melanocytic nevi and primary cutaneous melanomas were immunohistochemically stained for tubulin-ß-3 and 5-hmC. Immunoreactivity and staining patterns were correlated with Breslow-thickness, clinical and pathological characteristics, and progression-free survival. Melanocytes showed positive tubulin ß-3 staining. However, in most nevi, tubulin ß-3 staining appeared as a gradient with intense cytoplasmic staining in cells of the superficial part of the lesion that faded to weak staining in the deep dermal part, while no gradient was found in deep penetrating nevi and melanomas. In 53 % of the melanomas, areas with loss of tubulin ß-3 staining were found. 5-hmC staining intensity was significantly higher in melanocytic nevi compared to melanomas. Breslow thickness in combination with low 5-hmC score and loss of tubulin-ß-3 staining was predictive for poor prognosis. As single markers, tubulin-ß-3 and 5-hmC can be useful to distinguish between melanocytic nevi and melanoma, but staining variability limits the use of 5-hmC. In melanomas measuring >1.5 mm, combination of low 5-hmC score and loss of tubulin-ß-3 staining may have prognostic value.

Loss of Uhrf1 in neural stem cells leads to activation of retroviral elements and delayed neurodegeneration.

In order to understand whether early epigenetic mechanisms instruct the long-term behavior of neural stem cells (NSCs) and their progeny, we examined Uhrf1 (ubiquitin-like PHD ring finger-1; also known as Np95), as it is highly expressed in NSCs of the developing brain and rapidly down-regulated upon differentiation. Conditional deletion of Uhrf1 in the developing cerebral cortex resulted in rather normal proliferation and neurogenesis but severe postnatal neurodegeneration. During development, deletion of Uhrf1 lead to global DNA hypomethylation with a strong activation of the intracisternal A particle (IAP) family of endogenous retroviral elements, accompanied by an increase in 5-hydroxymethylcytosine. Down-regulation of Tet enzymes rescued the IAP activation in Uhrf1 conditional knockout (cKO) cells, suggesting an antagonistic interplay between Uhrf1 and Tet on IAP regulation. As IAP up-regulation persists into postnatal stages in the Uhrf1 cKO mice, our data show the lack of means to repress IAPs in differentiating neurons that normally never express Uhrf1 The high load of viral proteins and other transcriptional deregulation ultimately led to postnatal neurodegeneration. Taken together, these data show that early developmental NSC factors can have long-term effects in neuronal differentiation and survival. Moreover, they highlight how specific the consequences of widespread changes in DNA methylation are for certain classes of retroviral elements.

Distinctive Patterns of 5-Methylcytosine and 5-Hydroxymethylcytosine in Schizophrenia.

Schizophrenia is a highly heritable neuropsychiatric disorder characterized by cognitive and social dysfunction. Genetic, epigenetic, and environmental factors are together implicated in the pathogenesis and development of schizophrenia. DNA methylation, 5-methycytosine (5mC) and 5-hydroxylcytosine (5hmC) have been recognized as key epigenetic elements in neurodevelopment, ageing, and neurodegenerative diseases. Recently, distinctive 5mC and 5hmC pattern and expression changes of related genes have been discovered in schizophrenia. Antipsychotic drugs that affect 5mC status can alleviate symptoms in patients with schizophrenia, suggesting a critical role for DNA methylation in the pathogenesis of schizophrenia. Further exploring the signatures of 5mC and 5hmC in schizophrenia and developing precision-targeted epigenetic drugs based on this will provide new insights into the diagnosis and treatment of schizophrenia.

A Bayesian hierarchical model for improving measurement of 5mC and 5hmC levels: Toward revealing associations between phenotypes and methylation states.

5-hydroxymethylcytosine (5hmC) is a methylation state linked with gene regulation, commonly found in cells of the central nervous system. 5hmC is associated with demethylation of cytosines from 5-methylcytosine (5mC) to the unmethylated state. The presence of 5hmC can be inferred by a paired experiment involving bisulfite and oxidation-bisulfite treatments on the same sample, followed by a methylation assay using a platform such as the Illumina Infinium MethylationEPIC BeadChip (EPIC). Existing methods for analysis of the resulting EPIC data are not ideal. Most approaches ignore the correlation between the two experiments and any imprecision associated with DNA damage from the additional treatment. Estimates of 5mC/5hmC levels free from these limitations are desirable to reveal associations between methylation states and phenotypes. We propose a hierarchical Bayesian method called Constrained HYdroxy Methylation Estimation (CHYME) to simultaneously estimate 5mC/5hmC signals as well as any associations between these signals and covariates or phenotypes, while accounting for the potential impact of DNA damage and dependencies induced by the experimental design. Simulations show that CHYME has valid type 1 error and better power than a range of alternative methods, including the popular OxyBS method and linear models on transformed proportions. Other methods we examined suffer from hugely inflated type 1 error for inference on 5hmC proportions. We use CHYME to explore genome-wide associations between 5mC/5hmC levels and cause of death in postmortem prefrontal cortex brain tissue samples. These analyses indicate that CHYME is a useful tool to reveal phenotypic associations with 5mC/5hmC levels.

TET2-mediated upregulation of 5-hydroxymethylcytosine in LRRC39 promoter promotes Th1 response in association with downregulated Treg response in Vogt-Koyanagi-Harada disease.

DNA 5-Hydroxymethylcytosine (5-hmC), an oxidative reaction mediated by the ten-eleven translocation (TET) family, has been reported to play an essential role in the progression of auto-inflammatory and autoimmune diseases. By far, little is known about the effect of DNA 5-hmC and the TET family on the development of Vogt-Koyanagi-Harada (VKH) disease. In this study, we discovered that the global DNA 5-hmC level and the TET activity were elevated in association with the up-regulated expression of TET2 at both mRNA and protein levels in CD4+T cells from active VKH patients compared to healthy controls. Integrated analysis of DNA 5-hmC pattern and transcription profile of CD4+ T cells revealed that 6 candidate target genes were involved in the development of VKH disease. The promoter 5-hmC and mRNA levels of leucine rich repeat containing 39 (LRRC39) were verified to be elevated in active VKH patients. Functional experiments showed that TET2 could up-regulate LRRC39 mRNA expression by increasing the promoter 5-hmC level of LRRC39 in CD4+ T cells from active VKH patients. Up-regulated LRRC39 expression could increase the frequencies of IFN-γ+ and IL-17+ CD4+ T cells as well as the secretions of IFN-γ and IL-17 in association with the decreased frequency of CD4+CD25+FOXP3+ regulatory T (Treg) cells and the reduced production of IL-10. Additionally, restoration of LRRC39 rescued TET2-silencing-mediated reduced frequency of IFN-γ+ CD4+ T cells and increased frequency of CD4+CD25+FOXP3+ Treg cells. Collectively, our study reveals a novel axis, the TET2-5-hmC-LRRC39-Th1/Treg responses axis, in the pathogenesis of VKH and provides a potential target for further investigation into the epigenetic therapy of this disease.

Quantitative classification of chromatin dynamics reveals regulators of intestinal stem cell differentiation.

Intestinal stem cell (ISC) plasticity is thought to be regulated by broadly permissive chromatin shared between ISCs and their progeny. Here, we have used a Sox9EGFP reporter to examine chromatin across ISC differentiation. We find that open chromatin regions (OCRs) can be defined as broadly permissive or dynamic in a locus-specific manner, with dynamic OCRs found primarily in loci consistent with distal enhancers. By integrating gene expression with chromatin accessibility at transcription factor (TF) motifs in the context of Sox9EGFP populations, we classify broadly permissive and dynamic chromatin relative to TF usage. These analyses identify known and potential regulators of ISC differentiation via association with dynamic changes in chromatin. Consistent with computational predictions, Id3-null mice exhibit increased numbers of cells expressing the ISC-specific biomarker OLFM4. Finally, we examine the relationship between gene expression and 5-hydroxymethylcytosine (5hmC) in Sox9EGFP populations, which reveals 5hmC enrichment in absorptive lineage-specific genes. Our data demonstrate that intestinal chromatin dynamics can be quantitatively defined in a locus-specific manner, identify novel potential regulators of ISC differentiation and provide a chromatin roadmap for further dissecting cis regulation of cell fate in the intestine.

Transition of allele-specific DNA hydroxymethylation at regulatory loci is associated with phenotypic variation in monozygotic twins discordant for psychiatric disorders.

BACKGROUND: Major psychiatric disorders such as schizophrenia (SCZ) and bipolar disorder (BPD) are complex genetic mental illnesses. Their non-Mendelian features, such as those observed in monozygotic twins discordant for SCZ or BPD, are likely complicated by environmental modifiers of genetic effects. 5-Hydroxymethylcytosine (5hmC) is an important epigenetic mark in gene regulation, and whether it is linked to genetic variants that contribute to non-Mendelian features remains largely unexplored. METHODS: We combined the 5hmC-selective chemical labeling method (5hmC-seq) and whole-genome sequencing (WGS) analysis of peripheral blood DNA obtained from monozygotic (MZ) twins discordant for SCZ or BPD to identify allelic imbalances in hydroxymethylome maps, and examined association of allele-specific hydroxymethylation (AShM) transition with disease susceptibility based on Bayes factors (BF) derived from the Bayesian generalized additive linear mixed model. We then performed multi-omics integrative analysis to determine the molecular pathogenic basis of those AShM sites. We finally employed luciferase reporter, CRISPR/Cas9 technology, electrophoretic mobility shift assay (EMSA), chromatin immunoprecipitation (ChIP), PCR, FM4-64 imaging analysis, and RNA sequencing to validate the function of interested AShM sites in the human neuroblastoma SK-N-SH cells and human embryonic kidney 293T (HEK293T) cells. RESULTS: We identified thousands of genetic variants associated with AShM imbalances that exhibited phenotypic variation-associated AShM changes at regulatory loci. These AShM marks showed plausible associations with SCZ or BPD based on their effects on interactions among transcription factors (TFs), DNA methylation levels, or other epigenomic marks and thus contributed to dysregulated gene expression, which ultimately increased disease susceptibility. We then validated that competitive binding of POU3F2 on the alternative allele at the AShM site rs4558409 (G/T) in PLLP-enhanced PLLP expression, while the hydroxymethylated alternative allele, which alleviated the POU3F2 binding activity at the rs4558409 site, might be associated with the downregulated PLLP expression observed in BPD or SCZ. Moreover, disruption of rs4558409 promoted neural development and vesicle trafficking. CONCLUSION: Our study provides a powerful strategy for prioritizing regulatory risk variants and contributes to our understanding of the interplay between genetic and epigenetic factors in mediating SCZ or BPD susceptibility.