DNA-binding factors shape the mouse methylome at distal regulatory regions

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.

Skeletal muscle feature of different populations in large yellow croaker (Larimichthys crocea): from an epigenetic point of view

Fish skeletal muscle is composed of well-defined fiber types. In order to identify potential candidate genes affecting muscle growth and development under epigenetic regulation. Bisulfite sequencing was utilized to analyze and compare the muscle DNA methylation profiles of Larimichthys crocea inhabiting different environments. The results revealed that DNA methylation in L. crocea was predominantly CG methylation, with 2,396 differentially methylated regions (DMRs) identified through comparisons among different populations. The largest difference in methylation was observed between the ZhouShan and JinMen wild populations, suggesting that L. crocea may have undergone selection and domestication. Additionally, GO and KEGG enrichment analysis of differentially methylated genes (DMGs) revealed 626 enriched GO functional categories, including various muscle-related genes such as myh10, myf5, myf6, ndufv1, klhl31, map3k4, syn2b, sostdc1a, bag4, and hsp90ab. However, significant enrichment in KEGG pathways was observed only in the JinMen and XiangShan populations of L. crocea. Therefore, this study provides a theoretical foundation for a better understanding of the epigenetic regulation of skeletal muscle growth and development in L. crocea under different environmental conditions.

Oncoprotein LAMTOR5-mediated CHOP silence via DNA hypermethylation and miR-182/miR-769 in promotion of liver cancer growth

C/EBP homologous protein (CHOP) triggers the death of multiple cancers via endoplasmic reticulum (ER) stress. However, the function and regulatory mechanism of CHOP in liver cancer remain elusive. We have reported that late endosomal/lysosomal adapter, mitogen-activated protein kinase and mTOR activator 5 (LAMTOR5) suppresses apoptosis in various cancers. Here, we show that the transcriptional and posttranscriptional inactivation of CHOP mediated by LAMTOR5 accelerates liver cancer growth. Clinical bioinformatic analysis revealed that the expression of CHOP was low in liver cancer tissues and that its increased expression predicted a good prognosis. Elevated CHOP contributed to destruction of LAMTOR5-induced apoptotic suppression and proliferation. Mechanistically, LAMTOR5-recruited DNA methyltransferase 1 (DNMT1) to the CpG3 region (-559/-429) of the CHOP promoter and potentiated its hypermethylation to block its interaction with general transcription factor IIi (TFII-I), resulting in its inactivation. Moreover, LAMTOR5-enhanced miR-182/miR-769 reduced CHOP expression by targeting its 3'UTR. Notably, lenvatinib, a first-line targeted therapy for liver cancer, could target the LAMTOR5/CHOP axis to prevent liver cancer progression. Accordingly, LAMTOR5-mediated silencing of CHOP via the regulation of ER stress-related apoptosis promotes liver cancer growth, providing a theoretical basis for the use of lenvatinib for the treatment of liver cancer.

C2H2 proteins: Evolutionary aspects of domain architecture and diversification

The largest group of transcription factors in higher eukaryotes are C2H2 proteins, which contain C2H2-type zinc finger domains that specifically bind to DNA. Few well-studied C2H2 proteins, however, demonstrate their key role in the control of gene expression and chromosome architecture. Here we review the features of the domain architecture of C2H2 proteins and the likely origin of C2H2 zinc fingers. A comprehensive investigation of proteomes for the presence of proteins with multiple clustered C2H2 domains has revealed a key difference between groups of organisms. Unlike plants, transcription factors in metazoans contain clusters of C2H2 domains typically separated by a linker with the TGEKP consensus sequence. The average size of C2H2 clusters varies substantially, even between genomes of higher metazoans, and with a tendency to increase in combination with SCAN, and especially KRAB domains, reflecting the increasing complexity of gene regulatory networks.

Maternal Hypermethylated Genes Contribute to Intrauterine Growth Retardation of Piglets in Rongchang Pigs

The placenta is a crucial determinant of fetal survival, growth, and development. Deficiency in placental development directly causes intrauterine growth retardation (IUGR). IUGR can lead to fetal growth restriction and an increase in the mortality rate. The genetic mechanisms underlying IUGR development, however, remain unclear. In the present study, we integrated whole-genome DNA methylation and transcriptomic analyses to determine distinct gene expression patterns in various placental tissues to identify pivotal genes that are implicated with IUGR development. By performing RNA-sequencing analysis, 1487 differentially expressed genes (DEGs), with 737 upregulated and 750 downregulated genes, were identified in IUGR pigs (H_IUGR) compared with that in normal birth weight pigs (N_IUGR) (p < 0.05); furthermore, 77 miRNAs, 1331 lncRNAs, and 61 circRNAs were differentially expressed. The protein-protein interaction network analysis revealed that among these DEGs, the genes GNGT1, ANXA1, and CDC20 related to cellular developmental processes and blood vessel development were the key genes associated with the development of IUGR. A total of 495,870 differentially methylated regions were identified between the N_IUGR and H_IUGR groups, which included 25,053 differentially methylated genes (DMEs); moreover, the overall methylation level was higher in the H_IUGR group than in the N_IUGR group. Combined analysis showed an inverse correlation between methylation levels and gene expression. A total of 1375 genes involved in developmental processes, tissue development, and immune system regulation exhibited methylation differences in gene expression levels in the promoter regions and gene ontology regions. Five genes, namely, ANXA1, ADM, NRP2, SHH, and SMAD1, with high methylation levels were identified as potential contributors to IUGR development. These findings provide valuable insights that DNA methylation plays a crucial role in the epigenetic regulation of gene expression and mammalian development and that DNA-hypermethylated genes contribute to IUGR development in Rongchang pigs.

Genome-wide DNA methylation changes in human spermatogenesis

Sperm production and function require the correct establishment of DNA methylation patterns in the germline. Here, we examined the genome-wide DNA methylation changes during human spermatogenesis and its alterations in disturbed spermatogenesis. We found that spermatogenesis is associated with remodeling of the methylome, comprising a global decline in DNA methylation in primary spermatocytes followed by selective remethylation, resulting in a spermatids/sperm-specific methylome. Hypomethylated regions in spermatids/sperm were enriched in specific transcription factor binding sites for DMRT and SOX family members and spermatid-specific genes. Intriguingly, while SINEs displayed differential methylation throughout spermatogenesis, LINEs appeared to be protected from changes in DNA methylation. In disturbed spermatogenesis, germ cells exhibited considerable DNA methylation changes, which were significantly enriched at transposable elements and genes involved in spermatogenesis. We detected hypomethylation in SVA and L1HS in disturbed spermatogenesis, suggesting an association between the abnormal programming of these regions and failure of germ cells progressing beyond meiosis.

DNA methylation patterns of transcription factor binding regions characterize their functional and evolutionary contexts

BACKGROUND

DNA methylation is an important epigenetic modification which has numerous roles in modulating genome function. Its levels are spatially correlated across the genome, typically high in repressed regions but low in transcription factor (TF) binding sites and active regulatory regions. However, the mechanisms establishing genome-wide and TF binding site methylation patterns are still unclear.

RESULTS

Here we use a comparative approach to investigate the association of DNA methylation to TF binding evolution in mammals. Specifically, we experimentally profile DNA methylation and combine this with published occupancy profiles of five distinct TFs (CTCF, CEBPA, HNF4A, ONECUT1, FOXA1) in the liver of five mammalian species (human, macaque, mouse, rat, dog). TF binding sites are lowly methylated, but they often also have intermediate methylation levels. Furthermore, biding sites are influenced by the methylation status of CpGs in their wider binding regions even when CpGs are absent from the core binding motif. Employing a classification and clustering approach, we extract distinct and species-conserved patterns of DNA methylation levels at TF binding regions. CEBPA, HNF4A, ONECUT1, and FOXA1 share the same methylation patterns, while CTCF's differ. These patterns characterize alternative functions and chromatin landscapes of TF-bound regions. Leveraging our phylogenetic framework, we find DNA methylation gain upon evolutionary loss of TF occupancy, indicating coordinated evolution. Furthermore, each methylation pattern has its own evolutionary trajectory reflecting its genomic contexts.

CONCLUSIONS

Our epigenomic analyses indicate a role for DNA methylation in TF binding changes across species including that specific DNA methylation profiles characterize TF binding and are associated with their regulatory activity, chromatin contexts, and evolutionary trajectories.

Whole-Genome DNA Methylation Profiling of Intrahepatic Cholangiocarcinoma Reveals Prognostic Subtypes with Distinct Biological Drivers

Intrahepatic cholangiocarcinoma (iCCA) is the second most prevalent primary liver cancer. Although the genetic characterization of iCCA has led to targeted therapies for treating tumors with FGFR2 alterations and IDH1/2 mutations, only a limited number of patients can benefit from these strategies. Epigenomic profiles have emerged as potential diagnostic and prognostic biomarkers for improving the treatment of cancers. In this study, we conducted whole-genome bisulfite sequencing on 331 iCCAs integrated with genetic, transcriptomic, and proteomic analyses, demonstrating the existence of four DNA methylation subtypes of iCCAs (S1-S4) that exhibited unique postoperative clinical outcomes. The S1 group was an IDH1/2 mutation-specific subtype with moderate survival. The S2 subtype was characterized by the lowest methylation level and the highest mutational burden among the four subtypes and displayed upregulation of a gene-expression pattern associated with cell cycle/DNA replication. The S3 group was distinguished by high interpatient heterogeneity of tumor immunity, a gene-expression pattern associated with carbohydrate metabolism, and an enrichment of KRAS alterations. Patients with the S2 and S3 subtypes had the shortest survival among the four subtypes. Tumors in the S4 subtype, which had the best prognosis, showed global methylation levels comparable to normal controls, increased FGFR2 fusions/BAP1 mutations, and the highest copy-number variant burdens. Further integrative and functional analyses identified GBP4 demethylation, which is highly prevalent in the S2 and S3 groups, as an epigenetic oncogenic factor that regulates iCCA proliferation, migration, and invasion. Together, this study identifies prognostic methylome alterations and epigenetic drivers in iCCA.

SIGNIFICANCE

Characterization of the DNA methylome of intrahepatic cholangiocarcinoma integrated with genomic, transcriptomic, and proteomic analyses uncovers molecular mechanisms affected by genome-wide DNA methylation alterations, providing a resource for identifying potential therapeutic targets.

Systematic assessment of ISWI subunits shows that NURF creates local accessibility for CTCF

Catalytic activity of the imitation switch (ISWI) family of remodelers is critical for nucleosomal organization and DNA binding of certain transcription factors, including the insulator protein CTCF. Here we define the contribution of individual subcomplexes by deriving a panel of isogenic mouse stem cell lines, each lacking one of six ISWI accessory subunits. Individual deletions of subunits of either CERF, RSF, ACF, WICH or NoRC subcomplexes only moderately affect the chromatin landscape, while removal of the NURF-specific subunit BPTF leads to a strong reduction in chromatin accessibility and SNF2H ATPase localization around CTCF sites. This affects adjacent nucleosome occupancy and CTCF binding. At a group of sites with reduced chromatin accessibility, CTCF binding persists but cohesin occupancy is reduced, resulting in decreased insulation. These results suggest that CTCF binding can be separated from its function as an insulator in nuclear organization and identify a specific role for NURF in mediating SNF2H localization and chromatin opening at bound CTCF sites.

Spatial orchestration of the genome: topological reorganisation during X-chromosome inactivation

Genomes are organised through hierarchical structures, ranging from local kilobase-scale cis-regulatory contacts to large chromosome territories. Most notably, (sub)-compartments partition chromosomes according to transcriptional activity, while topologically associating domains (TADs) define cis-regulatory landscapes. The inactive X chromosome in mammals has provided unique insights into the regulation and function of the three-dimensional (3D) genome. Concurrent with silencing of the majority of genes and major alterations of its chromatin state, the X chromosome undergoes profound spatial rearrangements at multiple scales. These include the emergence of megadomains, alterations of the compartment structure and loss of the majority of TADs. Moreover, the Xist locus, which orchestrates X-chromosome inactivation, has provided key insights into regulation and function of regulatory domains. This review provides an overview of recent insights into the control of these structural rearrangements and contextualises them within a broader understanding of 3D genome organisation.

The role of epigenetics in rare diseases

Epigenetic control systems are based on chromatin modifications (DNA methylation, histone modifications and nucleosome positioning), which affect the local kinetics of gene expression. They play an important role in maintaining cell fate decisions, X inactivation and genomic imprinting. Aberrant chromatin states that are associated with a deleterious change in gene expression are called epimutations. An epimutation can be a primary epimutation that has occurred in the absence of any genetic change or a secondary epimutation that results from a mutation of a cis-acting regulatory element or trans-acting factor. Epimutations may play a causative role in disease, for example in imprinting disorders, or may be part of the pathogenetic mechanism as in the fragile X syndrome and in syndromes caused by a mutation affecting a chromatin modifier. For several diseases, DNA methylation testing is an important tool in the diagnostic work-up of patients.

Mutation rate heterogeneity at the sub-gene scale due to local DNA hypomethylation

Local mutation rates in human are highly heterogeneous, with known variability at the scale of megabase-sized chromosomal domains, and, on the other extreme, at the scale of oligonucleotides. The intermediate, kilobase-scale heterogeneity in mutation risk is less well characterized. Here, by analyzing thousands of somatic genomes, we studied mutation risk gradients along gene bodies, representing a genomic scale spanning roughly 1-10 kb, hypothesizing that different mutational mechanisms are differently distributed across gene segments. The main heterogeneity concerns several kilobases at the transcription start site and further downstream into 5' ends of gene bodies; these are commonly hypomutated with several mutational signatures, most prominently the ubiquitous C > T changes at CpG dinucleotides. The width and shape of this mutational coldspot at 5' gene ends is variable across genes, and corresponds to variable interval of lowered DNA methylation depending on gene activity level and regulation. Such hypomutated loci, at 5' gene ends or elsewhere, correspond to DNA hypomethylation that can associate with various landmarks, including intragenic enhancers, Polycomb-marked regions, or chromatin loop anchor points. Tissue-specific DNA hypomethylation begets tissue-specific local hypomutation. Of note, direction of mutation risk is inverted for AID/APOBEC3 cytosine deaminase activity, whose signatures are enriched in hypomethylated regions.

Comprehensive analysis of the REST transcription factor regulatory networks in IDH mutant and IDH wild-type glioma cell lines and tumors

The RE1-silencing transcription factor (REST) acts either as a repressor or activator of transcription depending on the genomic and cellular context. REST is a key player in brain cell differentiation by inducing chromatin modifications, including DNA methylation, in a proximity of its binding sites. Its dysfunction may contribute to oncogenesis. Mutations in IDH1/2 significantly change the epigenome contributing to blockade of cell differentiation and glioma development. We aimed at defining how REST modulates gene activation and repression in the context of the IDH mutation-related phenotype in gliomas. We studied the effects of REST knockdown, genome wide occurrence of REST binding sites, and DNA methylation of REST motifs in IDH wild type and IDH mutant gliomas. We found that REST target genes, REST binding patterns, and TF motif occurrence proximal to REST binding sites differed in IDH wild-type and mutant gliomas. Among differentially expressed REST targets were genes involved in glial cell differentiation and extracellular matrix organization, some of which were differentially methylated at promoters or gene bodies. REST knockdown differently impacted invasion of the parental or IDH1 mutant glioma cells. The canonical REST-repressed gene targets showed significant correlation with the GBM NPC-like cellular state. Interestingly, results of REST or KAISO silencing suggested the interplay between these TFs in regulation of REST-activated and repressed targets. The identified gene regulatory networks and putative REST cooperativity with other TFs, such as KAISO, show distinct REST target regulatory networks in IDH-WT and IDH-MUT gliomas, without concomitant DNA methylation changes. We conclude that REST could be an important therapeutic target in gliomas.

Combinatorial quantification of 5mC and 5hmC at individual CpG dyads and the transcriptome in single cells reveals modulators of DNA methylation maintenance fidelity

Inheritance of 5-methylcytosine from one cell generation to the next by DNA methyltransferase 1 (DNMT1) plays a key role in regulating cellular identity. While recent work has shown that the activity of DNMT1 is imprecise, it remains unclear how the fidelity of DNMT1 is tuned in different genomic and cell state contexts. Here we describe Dyad-seq, a method to quantify the genome-wide methylation status of cytosines at the resolution of individual CpG dinucleotides to find that the fidelity of DNMT1-mediated maintenance methylation is related to the local density of DNA methylation and the landscape of histone modifications. To gain deeper insights into methylation/demethylation turnover dynamics, we first extended Dyad-seq to quantify all combinations of 5-methylcytosine and 5-hydroxymethylcytosine at individual CpG dyads. Next, to understand how cell state transitions impact maintenance methylation, we scaled the method down to jointly profile genome-wide methylation levels, maintenance methylation fidelity and the transcriptome from single cells (scDyad&T-seq). Using scDyad&T-seq, we demonstrate that, while distinct cell states can substantially impact the activity of the maintenance methylation machinery, locally there exists an intrinsic relationship between DNA methylation density, histone modifications and DNMT1-mediated maintenance methylation fidelity that is independent of cell state.

Tet-mediated DNA methylation dynamics affect chromosome organization

DNA Methylation is a significant epigenetic modification that can modulate chromosome states, but its role in orchestrating chromosome organization has not been well elucidated. Here we systematically assessed the effects of DNA Methylation on chromosome organization with a multi-omics strategy to capture DNA Methylation and high-order chromosome interaction simultaneously on mouse embryonic stem cells with DNA methylation dioxygenase Tet triple knock-out (Tet-TKO). Globally, upon Tet-TKO, we observed weakened compartmentalization, corresponding to decreased methylation differences between CpG island (CGI) rich and poor domains. Tet-TKO could also induce hypermethylation for the CTCF binding peaks in TAD boundaries and chromatin loop anchors. Accordingly, CTCF peak generally weakened upon Tet-TKO, which results in weakened TAD structure and depletion of long-range chromatin loops. Genes that lost enhancer-promoter looping upon Tet-TKO showed DNA hypermethylation in their gene bodies, which may compensate for the disruption of gene expression. We also observed distinct effects of Tet1 and Tet2 on chromatin organization and increased DNA methylation correlation on spatially interacted fragments upon Tet inactivation. Our work showed the broad effects of Tet inactivation and DNA methylation dynamics on chromosome organization.

Biallelic non-productive enhancer-promoter interactions precede imprinted expression of Kcnk9 during mouse neural commitment

It is only partially understood how constitutive allelic methylation at imprinting control regions (ICRs) interacts with other regulation levels to drive timely parental allele-specific expression along large imprinted domains. The Peg13-Kcnk9 domain is an imprinted domain with important brain functions. To gain insights into its regulation during neural commitment, we performed an integrative analysis of its allele-specific epigenetic, transcriptomic, and cis-spatial organization using a mouse stem cell-based corticogenesis model that recapitulates the control of imprinted gene expression during neurodevelopment. We found that, despite an allelic higher-order chromatin structure associated with the paternally CTCF-bound Peg13 ICR, enhancer-Kcnk9 promoter contacts occurred on both alleles, although they were productive only on the maternal allele. This observation challenges the canonical model in which CTCF binding isolates the enhancer and its target gene on either side and suggests a more nuanced role for allelic CTCF binding at some ICRs.

Epigenomic insights into common human disease pathology

The epigenome-the chemical modifications and chromatin-related packaging of the genome-enables the same genetic template to be activated or repressed in different cellular settings. This multi-layered mechanism facilitates cell-type specific function by setting the local sequence and 3D interactive activity level. Gene transcription is further modulated through the interplay with transcription factors and co-regulators. The human body requires this epigenomic apparatus to be precisely installed throughout development and then adequately maintained during the lifespan. The causal role of the epigenome in human pathology, beyond imprinting disorders and specific tumour suppressor genes, was further brought into the spotlight by large-scale sequencing projects identifying that mutations in epigenomic machinery genes could be critical drivers in both cancer and developmental disorders. Abrogation of this cellular mechanism is providing new molecular insights into pathogenesis. However, deciphering the full breadth and implications of these epigenomic changes remains challenging. Knowledge is accruing regarding disease mechanisms and clinical biomarkers, through pathogenically relevant and surrogate tissue analyses, respectively. Advances include consortia generated cell-type specific reference epigenomes, high-throughput DNA methylome association studies, as well as insights into ageing-related diseases from biological 'clocks' constructed by machine learning algorithms. Also, 3rd-generation sequencing is beginning to disentangle the complexity of genetic and DNA modification haplotypes. Cell-free DNA methylation as a cancer biomarker has clear clinical utility and further potential to assess organ damage across many disorders. Finally, molecular understanding of disease aetiology brings with it the opportunity for exact therapeutic alteration of the epigenome through CRISPR-activation or inhibition.

One-pot trimodal mapping of unmethylated, hydroxymethylated, and open chromatin sites unveils distinctive 5hmC roles at dynamic chromatin loci

We present a method, named Mx-TOP, for profiling of three epigenetic regulatory layers-chromatin accessibility, general DNA modification, and DNA hydroxymethylation-from a single library. The approach is based on chemo-enzymatic covalent tagging of unmodified CG sites and hydroxymethylated cytosine (5hmC) along with GC sites in chromatin, which are then mapped using tag-selective base-resolution TOP-seq sequencing. Our in-depth validation of the approach revealed its sensitivity and informativity in evaluating chromatin accessibility and DNA modification interactions that drive transcriptional regulation. We employed the technology in a study of chromatin and DNA demethylation dynamics during in vitro neuronal differentiation. The study highlighted the involvement of gene body 5hmC in modulating an extensive decoupling between promoter accessibility and transcription. The importance of 5hmC in chromatin remodeling was further demonstrated by the observed resistance of the developmentally acquired open loci to the global 5hmC erasure in neuronal progenitors.

DNA methylation-dependent and -independent binding of CDX2 directs activation of distinct developmental and homeostatic genes

Precise spatiotemporal and cell type-specific gene expression is essential for proper tissue development and function. Transcription factors (TFs) guide this process by binding to developmental stage-specific targets and establishing an appropriate enhancer landscape. In turn, DNA and chromatin modifications direct the genomic binding of TFs. However, how TFs navigate various chromatin features and selectively bind a small portion of the millions of possible genomic target loci is still not well understood. Here we show that Cdx2 - a pioneer TF that binds distinct targets in developing versus adult intestinal epithelial cells - has a preferential affinity for a non-canonical CpG-containing motif in vivo. A higher frequency of this motif at embryonic and fetal Cdx2 target loci and the specifically methylated state of the CpG during development allows selective Cdx2 binding and activation of developmental enhancers and linked genes. Conversely, demethylation at these enhancers prohibits ectopic Cdx2 binding in adult cells, where Cdx2 binds its canonical motif without a CpG. This differential Cdx2 binding allows for corecruitment of Ctcf and Hnf4, facilitating the establishment of intestinal superenhancers during development and enhancers mediating adult homeostatic functions, respectively. Induced gain of DNA methylation in the adult mouse epithelium or cultured cells causes ectopic recruitment of Cdx2 to the developmental target loci and facilitates cobinding of the partner TFs. Together, our results demonstrate that the differential CpG motif requirements for Cdx2 binding to developmental versus adult target sites allow it to navigate different DNA methylation profiles and activate cell type-specific genes at appropriate times.

Characterisation and reproducibility of the HumanMethylationEPIC v2.0 BeadChip for DNA methylation profiling

BACKGROUND

The Illumina family of Infinium Methylation BeadChip microarrays has been widely used over the last 15 years for genome-wide DNA methylation profiling, including large-scale and population-based studies, due to their ease of use and cost effectiveness. Succeeding the popular HumanMethylationEPIC BeadChip (EPICv1), the recently released Infinium MethylationEPIC v2.0 BeadChip (EPICv2) claims to extend genomic coverage to more than 935,000 CpG sites. Here, we comprehensively characterise the reproducibility, reliability and annotation of the EPICv2 array, based on bioinformatic analysis of both manifest data and new EPICv2 data from diverse biological samples.

RESULTS

We find a high degree of reproducibility with EPICv1, evidenced by comparable sensitivity and precision from empirical cross-platform comparison incorporating whole genome bisulphite sequencing (WGBS), and high correlation between technical sample replicates, including between samples with DNA input levels below the manufacturer's recommendation. We provide a full assessment of probe content, evaluating genomic distribution and changes from previous array versions. We characterise EPICv2's new feature of replicated probes and provide recommendations as to the superior probes. In silico analysis of probe sequences demonstrates that probe cross-hybridisation remains a significant problem in EPICv2. By mapping the off-target sites at single nucleotide resolution and comparing with WGBS we show empirical evidence for preferential off-target binding.

CONCLUSIONS

Overall, we find EPICv2 a worthy successor to the previous Infinium methylation microarrays, however some technical issues remain. To support optimal EPICv2 data analysis we provide an expanded version of the EPICv2 manifest to aid researchers in understanding probe design, data processing, choosing appropriate probes for analysis and for integration with methylation datasets from previous versions of the Infinium Methylation BeadChip.

The impact of DNA methylation on CTCF-mediated 3D genome organization

Cytosine DNA methylation is a highly conserved epigenetic mark in eukaryotes. Although the role of DNA methylation at gene promoters and repetitive elements has been extensively studied, the function of DNA methylation in other genomic contexts remains less clear. In the nucleus of mammalian cells, the genome is spatially organized at different levels, and strongly influences myriad genomic processes. There are a number of factors that regulate the three-dimensional (3D) organization of the genome, with the CTCF insulator protein being among the most well-characterized. Pertinently, CTCF binding has been reported as being DNA methylation-sensitive in certain contexts, perhaps most notably in the process of genomic imprinting. Therefore, it stands to reason that DNA methylation may play a broader role in the regulation of chromatin architecture. Here we summarize the current understanding that is relevant to both the mammalian DNA methylation and chromatin architecture fields and attempt to assess the extent to which DNA methylation impacts the folding of the genome. The focus is in early embryonic development and cellular transitions when the epigenome is in flux, but we also describe insights from pathological contexts, such as cancer, in which the epigenome and 3D genome organization are misregulated.

Epigenetic disruption of the RARγ complex impairs its function to bookmark AR enhancer interactions required for enzalutamide sensitivity in prostate cancer

The current study in prostate cancer (PCa) focused on the genomic mechanisms at the cross-roads of pro-differentiation signals and the emergence of lineage plasticity. We explored an understudied cistromic mechanism involving RARγ's ability to govern AR cistrome-transcriptome relationships, including those associated with more aggressive PCa features. The RARγ complex in PCa cell models was enriched for canonical cofactors, as well as proteins involved in RNA processing and bookmarking. Identifying the repertoire of miR-96 bound and regulated gene targets, including those recognition elements marked by m6A, revealed their significant enrichment in the RARγ complex. RARγ significantly enhanced the AR cistrome, particularly in active enhancers and super-enhancers, and overlapped with the binding of bookmarking factors. Furthermore, RARγ expression led to nucleosome-free chromatin enriched with H3K27ac, and significantly enhanced the AR cistrome in G2/M cells. RARγ functions also antagonized the transcriptional actions of the lineage master regulator ONECUT2. Similarly, gene programs regulated by either miR-96 or antagonized by RARγ were enriched in alternative lineages and more aggressive PCa phenotypes. Together these findings reveal an under-investigated role for RARγ, modulated by miR-96, to bookmark enhancer sites during mitosis. These sites are required by the AR to promote transcriptional competence, and emphasize luminal differentiation, while antagonizing ONECUT2.

scMeFormer: a transformer-based deep learning model for imputing DNA methylation states in single cells enhances the detection of epigenetic alterations in schizophrenia

DNA methylation (DNAm), a crucial epigenetic mark, plays a key role in gene regulation, mammalian development, and various human diseases. Single-cell technologies enable the profiling of DNAm states at cytosines within the DNA sequence of individual cells, but they often suffer from limited coverage of CpG sites. In this study, we introduce scMeFormer, a transformer-based deep learning model designed to impute DNAm states for each CpG site in single cells. Through comprehensive evaluations, we demonstrate the superior performance of scMeFormer compared to alternative models across four single-nucleus DNAm datasets generated by distinct technologies. Remarkably, scMeFormer exhibits high-fidelity imputation, even when dealing with significantly reduced coverage, as low as 10% of the original CpG sites. Furthermore, we applied scMeFormer to a single-nucleus DNAm dataset generated from the prefrontal cortex of four schizophrenia patients and four neurotypical controls. This enabled the identification of thousands of differentially methylated regions associated with schizophrenia that would have remained undetectable without imputation and added granularity to our understanding of epigenetic alterations in schizophrenia within specific cell types. Our study highlights the power of deep learning in imputing DNAm states in single cells, and we expect scMeFormer to be a valuable tool for single-cell DNAm studies.

Epigenomic Characterization of Lymphoid Neoplasms

Lymphoid neoplasms represent a heterogeneous group of disease entities and subtypes with markedly different molecular and clinical features. Beyond genetic alterations, lymphoid tumors also show widespread epigenomic changes. These severely affect the levels and distribution of DNA methylation, histone modifications, chromatin accessibility, and three-dimensional genome interactions. DNA methylation stands out as a tracer of cell identity and memory, as B cell neoplasms show epigenetic imprints of their cellular origin and proliferative history, which can be quantified by an epigenetic mitotic clock. Chromatin-associated marks are informative to uncover altered regulatory regions and transcription factor networks contributing to the development of distinct lymphoid tumors. Tumor-intrinsic epigenetic and genetic aberrations cooperate and interact with microenvironmental cells to shape the transcriptome at different phases of lymphoma evolution, and intraclonal heterogeneity can now be characterized by single-cell profiling. Finally, epigenetics offers multiple clinical applications, including powerful diagnostic and prognostic biomarkers as well as therapeutic targets.

Antiproliferative, apoptosis-inducing, and GSTP1 demethylation activities of Ellagitannins isolated from Cornus alba L

This study aimed to prove the prostate cancer chemopreventive activity of compounds isolated from CA. We evaluated these compounds using the 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) assay and evaluated their NF-κB inhibitory activity and apoptosis-inducing activity using western blot analysis and flow cytometry, respectively. Their DNA methylation activity was also evaluated via a methylation-specific polymerase chain reaction in androgen-dependent (LNCaP) and androgen-independent (PC-3) prostate cancer cell lines. Camptothin B (1), cornusiin B (2), and cornusiin A (3), which were isolated in our previous work, relatively reduced the protein expression levels in PCa cells. Among them, cornusiin B (2) exhibited excellent NF-κB inhibitory activity. Also, concentration-dependently increased the unmethylated DNA content and decreased the methylated DNA content in both PC-3 and LNCaP cells. Therefore, cornusiin B (2), which was isolated from CA, has the potential to act as a chemopreventive agent for prostate cancer.

A DNA methylation haplotype block landscape in human tissues and preimplantation embryos reveals regulatory elements defined by comethylation patterns

DNA methylation and associated regulatory elements play a crucial role in gene expression regulation. Previous studies have focused primarily on the distribution of mean methylation levels. Advances in whole-genome bisulfite sequencing (WGBS) have enabled the characterization of DNA methylation haplotypes (MHAPs), representing CpG sites from the same read fragment on a single chromosome, and the subsequent identification of methylation haplotype blocks (MHBs), in which adjacent CpGs on the same fragment are comethylated. Using our expert-curated WGBS data sets, we report comprehensive landscapes of MHBs in 17 representative normal somatic human tissues and during early human embryonic development. Integrative analysis reveals MHBs as a distinctive type of regulatory element characterized by comethylation patterns rather than mean methylation levels. We show the enrichment of MHBs in open chromatin regions, tissue-specific histone marks, and enhancers, including super-enhancers. Moreover, we find that MHBs tend to localize near tissue-specific genes and show an association with differential gene expression that is independent of mean methylation. Similar findings are observed in the context of human embryonic development, highlighting the dynamic nature of MHBs during early development. Collectively, our comprehensive MHB landscapes provide valuable insights into the tissue specificity and developmental dynamics of DNA methylation.

The chromatin signatures of enhancers and their dynamic regulation

Enhancers are cis-regulatory elements that can stimulate gene expression from distance, and drive precise spatiotemporal gene expression profiles during development. Functional enhancers display specific features including an open chromatin conformation, Histone H3 lysine 27 acetylation, Histone H3 lysine 4 mono-methylation enrichment, and enhancer RNAs production. These features are modified upon developmental cues which impacts their activity. In this review, we describe the current state of knowledge about enhancer functions and the diverse chromatin signatures found on enhancers. We also discuss the dynamic changes of enhancer chromatin signatures, and their impact on lineage specific gene expression profiles, during development or cellular differentiation.

Pleiotropic influence of DNA methylation QTLs on physiological and ageing traits

DNA methylation is influenced by genetic and non-genetic factors. Here, we chart quantitative trait loci (QTLs) that modulate levels of methylation at highly conserved CpGs using liver methylome data from mouse strains belonging to the BXD family. A regulatory hotspot on chromosome 5 had the highest density of trans-acting methylation QTLs (trans-meQTLs) associated with multiple distant CpGs. We refer to this locus as meQTL.5a. Trans-modulated CpGs showed age-dependent changes and were enriched in developmental genes, including several members of the MODY pathway (maturity onset diabetes of the young). The joint modulation by genotype and ageing resulted in a more 'aged methylome' for BXD strains that inherited the DBA/2J parental allele at meQTL.5a. Further, several gene expression traits, body weight, and lipid levels mapped to meQTL.5a, and there was a modest linkage with lifespan. DNA binding motif and protein-protein interaction enrichment analyses identified the hepatic nuclear factor, Hnf1a (MODY3 gene in humans), as a strong candidate. The pleiotropic effects of meQTL.5a could contribute to variations in body size and metabolic traits, and influence CpG methylation and epigenetic ageing that could have an impact on lifespan.

MSL2 ensures biallelic gene expression in mammals

In diploid organisms, biallelic gene expression enables the production of adequate levels of mRNA1,2. This is essential for haploinsufficient genes, which require biallelic expression for optimal function to prevent the onset of developmental disorders1,3. Whether and how a biallelic or monoallelic state is determined in a cell-type-specific manner at individual loci remains unclear. MSL2 is known for dosage compensation of the male X chromosome in flies. Here we identify a role of MSL2 in regulating allelic expression in mammals. Allele-specific bulk and single-cell analyses in mouse neural progenitor cells revealed that, in addition to the targets showing biallelic downregulation, a class of genes transitions from biallelic to monoallelic expression after MSL2 loss. Many of these genes are haploinsufficient. In the absence of MSL2, one allele remains active, retaining active histone modifications and transcription factor binding, whereas the other allele is silenced, exhibiting loss of promoter-enhancer contacts and the acquisition of DNA methylation. Msl2-knockout mice show perinatal lethality and heterogeneous phenotypes during embryonic development, supporting a role for MSL2 in regulating gene dosage. The role of MSL2 in preserving biallelic expression of specific dosage-sensitive genes sets the stage for further investigation of other factors that are involved in allelic dosage compensation in mammalian cells, with considerable implications for human disease.

Predicting active enhancers with DNA methylation and histone modification

BACKGROUND

Enhancers play a crucial role in gene regulation, and some active enhancers produce noncoding RNAs known as enhancer RNAs (eRNAs) bi-directionally. The most commonly used method for detecting eRNAs is CAGE-seq, but the instability of eRNAs in vivo leads to data noise in sequencing results. Unfortunately, there is currently a lack of research focused on the noise inherent in CAGE-seq data, and few approaches have been developed for predicting eRNAs. Bridging this gap and developing widely applicable eRNA prediction models is of utmost importance.

RESULTS

In this study, we proposed a method to reduce false positives in the identification of eRNAs by adjusting the statistical distribution of expression levels. We also developed eRNA prediction models using joint gene expressions, DNA methylation, and histone modification. These models achieved impressive performance with an AUC value of approximately 0.95 for intra-cell prediction and 0.9 for cross-cell prediction.

CONCLUSIONS

Our method effectively attenuates the noise generated by stochastic RNA production, resulting in more accurate detection of eRNAs. Furthermore, our eRNA prediction model exhibited significant accuracy in both intra-cell and cross-cell validation, highlighting its robustness and potential application in various cellular contexts.

Domains Rearranged Methylase 2 maintains DNA methylation at large DNA hypomethylated shores and long-range chromatin interactions in rice

DNA methylation plays an important role in gene regulation and genomic stability. However, large DNA hypomethylated regions known as DNA methylation valleys (DMVs) or canyons have also been suggested to serve unique regulatory functions, largely unknown in rice (Oryza sativa). Here, we describe the DMVs in rice seedlings, which were highly enriched with developmental and transcription regulatory genes. Further detailed analysis indicated that grand DMVs (gDMVs) might be derived from nuclear integrants of organelle DNA (NORGs). Furthermore, Domains Rearranged Methylase 2 (OsDRM2) maintained DNA methylation at short DMV (sDMV) shores. Epigenetic maps indicated that sDMVs were marked with H3K4me3 and/or H3K27me3, although the loss of DNA methylation had a negligible effect on histone modification within these regions. In addition, we constructed H3K27me3-associated interaction maps for homozygous T-DNA insertion mutant of the gene (osdrm2) and wild type (WT). From a global perspective, most (90%) compartments were stable between osdrm2 and WT plants. At a high resolution, we observed a dramatic loss of long-range chromatin loops in osdrm2, which suffered an extensive loss of non-CG (CHG and CHH, H = A, T, or C) methylation. From another viewpoint, the loss of non-CG methylation at sDMV shores in osdrm2 could disrupt H3K27me3-mediated chromatin interaction networks. Overall, our results demonstrated that DMVs are a key genomic feature in rice and are precisely regulated by epigenetic modifications, including DNA methylation and histone modifications. OsDRM2 maintained DNA methylation at sDMV shores, while OsDRM2 deficiency strongly affected three-dimensional (3D) genome architectures.

From Genotype to Phenotype: How Enhancers Control Gene Expression and Cell Identity in Hematopoiesis

Blood comprises a wide array of specialized cells, all of which share the same genetic information and ultimately derive from the same precursor, the hematopoietic stem cell (HSC). This diversity of phenotypes is underpinned by unique transcriptional programs gradually acquired in the process known as hematopoiesis. Spatiotemporal regulation of gene expression depends on many factors, but critical among them are enhancers-sequences of DNA that bind transcription factors and increase transcription of genes under their control. Thus, hematopoiesis involves the activation of specific enhancer repertoires in HSCs and their progeny, driving the expression of sets of genes that collectively determine morphology and function. Disruption of this tightly regulated process can have catastrophic consequences: in hematopoietic malignancies, dysregulation of transcriptional control by enhancers leads to misexpression of oncogenes that ultimately drive transformation. This review attempts to provide a basic understanding of enhancers and their role in transcriptional regulation, with a focus on normal and malignant hematopoiesis. We present examples of enhancers controlling master regulators of hematopoiesis and discuss the main mechanisms leading to enhancer dysregulation in leukemia and lymphoma.

[Epigenetics, principles and examples of applications]

Since the discovery of DNA as the support of genetic information, the challenge for generations of life scientists was to understand the mechanisms underlying the process that translate the sequence of a gene to a phenotype. In the 1950s, the concept of epigenetics was defined by the British biologist Conrad H. Waddington as the study of "epigenesis" that governs the biological processes involved in the development of any organism. The term epigenetics, now best defined as "above the DNA sequence" reflects the gene-environment interactions by which genes determine traits. Since, its first description, studies underlying the mechanisms involved in these processes has led to an increasing understanding of the regulation all genome transactions such as transcription, replication, repair and the biological pathways coordinated by these mechanisms. We will discuss here the main principles regulating epigenetic processes, their roles in physiology, their evolution over the life time and their implications in medicine.

Widespread effects of DNA methylation and intra-motif dependencies revealed by novel transcription factor binding models

Several studies suggested that transcription factor (TF) binding to DNA may be impaired or enhanced by DNA methylation. We present MeDeMo, a toolbox for TF motif analysis that combines information about DNA methylation with models capturing intra-motif dependencies. In a large-scale study using ChIP-seq data for 335 TFs, we identify novel TFs that show a binding behaviour associated with DNA methylation. Overall, we find that the presence of CpG methylation decreases the likelihood of binding for the majority of methylation-associated TFs. For a considerable subset of TFs, we show that intra-motif dependencies are pivotal for accurately modelling the impact of DNA methylation on TF binding. We illustrate that the novel methylation-aware TF binding models allow to predict differential ChIP-seq peaks and improve the genome-wide analysis of TF binding. Our work indicates that simplistic models that neglect the effect of DNA methylation on DNA binding may lead to systematic underperformance for methylation-associated TFs.

Predicting regional somatic mutation rates using DNA motifs

How the locus-specificity of epigenetic modifications is regulated remains an unanswered question. A contributing mechanism is that epigenetic enzymes are recruited to specific loci by DNA binding factors recognizing particular sequence motifs (referred to as epi-motifs). Using these motifs to predict biological outputs depending on local epigenetic state such as somatic mutation rates would confirm their functionality. Here, we used DNA motifs including known TF motifs and epi-motifs as a surrogate of epigenetic signals to predict somatic mutation rates in 13 cancers at an average 23kbp resolution. We implemented an interpretable neural network model, called contextual regression, to successfully learn the universal relationship between mutations and DNA motifs, and uncovered motifs that are most impactful on the regional mutation rates such as TP53 and epi-motifs associated with H3K9me3. Furthermore, we identified genomic regions with significantly higher mutation rates than the expected values in each individual tumor and demonstrated that such cancer-related regions can accurately predict cancer types. Interestingly, we found that the same mutation signatures often have different contributions to cancer-related and cancer-independent regions, and we also identified the motifs with the most contribution to each mutation signature.

Exploring the Functional Basis of Epigenetic Aging in Relation to Body Fat Phenotypes in the Norfolk Island Cohort

DNA methylation is an epigenetic factor that is modifiable and can change over a lifespan. While many studies have identified methylation sites (CpGs) related to aging, the relationship of these to gene function and age-related disease phenotypes remains unclear. This research explores this question by testing for the conjoint association of age-related CpGs with gene expression and the relation of these to body fat phenotypes. The study included blood-based gene transcripts and intragenic CpG methylation data from Illumina 450 K arrays in 74 healthy adults from the Norfolk Island population. First, a series of regression analyses were performed to detect associations between gene transcript level and intragenic CpGs and their conjoint relationship with age. Second, we explored how these age-related expression CpGs (eCpGs) correlated with obesity-related phenotypes, including body fat percentage, body mass index, and waist-to-hip ratio. We identified 35 age-related eCpGs associated with age. Of these, ten eCpGs were associated with at least one body fat phenotype. Collagen Type XI Alpha 2 Chain (COL11A2), Complement C1s (C1s), and four and a half LIM domains 2 (FHL2) genes were among the most significant genes with multiple eCpGs associated with both age and multiple body fat phenotypes. The COL11A2 gene contributes to the correct assembly of the extracellular matrix in maintaining the healthy structural arrangement of various components, with the C1s gene part of complement systems functioning in inflammation. Moreover, FHL2 expression was upregulated under hypermethylation in both blood and adipose tissue with aging. These results suggest new targets for future studies and require further validation to confirm the specific function of these genes on body fat regulation.

Regulation of the three-dimensional chromatin organization by transposable elements in pig spleen

Like other mammalian species, the pig genome is abundant with transposable elements (TEs). The importance of TEs for three-dimensional (3D) chromatin organization has been observed in species like human and mouse, yet current understanding about pig TEs is absent. Here, we investigated the contribution of TEs for the 3D chromatin organization in three pig tissues, focusing on spleen which is crucial for both adaptive and innate immunity. We identified dozens of TE families overrepresented with CTCF binding sites, including LTR22_SS, LTR15_SS and LTR16_SSc which are pig-specific families of endogenous retroviruses (ERVs). Interestingly, LTR22_SS elements harbor a CTCF motif and create hundreds of CTCF binding sites that are associated with adaptive immunity. We further applied Hi-C to profile the 3D chromatin structure in spleen and found that TE-derived CTCF binding sites correlate with chromatin insulation and frequently overlap TAD borders and loop anchors. Notably, one LTR22_SS-derived CTCF binding site demarcate a TAD boundary upstream of XCL1, which is a spleen-enriched chemokine gene important for lymphocyte trafficking and inflammation. Overall, this study represents a first step toward understanding the function of TEs on 3D chromatin organization regulation in pigs and expands our understanding about the functional importance of TEs in mammals.

Genic constraint against nonsynonymous variation across the mouse genome

BACKGROUND

Selective constraint, the depletion of variation due to negative selection, provides insights into the functional impact of variants and disease mechanisms. However, its characterization in mice, the most commonly used mammalian model, remains limited. This study aims to quantify mouse gene constraint using a new metric called the nonsynonymous observed expected ratio (NOER) and investigate its relationship with gene function.

RESULTS

NOER was calculated using whole-genome sequencing data from wild mouse populations (Mus musculus sp and Mus spretus). Positive correlations were observed between mouse gene constraint and the number of associated knockout phenotypes, indicating stronger constraint on pleiotropic genes. Furthermore, mouse gene constraint showed a positive correlation with the number of pathogenic variant sites in their human orthologues, supporting the relevance of mouse models in studying human disease variants.

CONCLUSIONS

NOER provides a resource for assessing the fitness consequences of genetic variants in mouse genes and understanding the relationship between gene constraint and function. The study's findings highlight the importance of pleiotropy in selective constraint and support the utility of mouse models in investigating human disease variants. Further research with larger sample sizes can refine constraint estimates in mice and enable more comprehensive comparisons of constraint between mouse and human orthologues.

The Change in Whole-Genome Methylation and Transcriptome Profile under Autophagy Defect and Nitrogen Starvation

Through whole-genome bisulfite sequencing and RNA-seq, we determined the potential impact of autophagy in regulating DNA methylation in Arabidopsis, providing a solid foundation for further understanding the molecular mechanism of autophagy and how plants cope with nitrogen deficiency. A total of 335 notable differentially expressed genes (DEGs) were discovered in wild-type Arabidopsis (Col-0-N) and an autophagic mutant cultivated under nitrogen starvation (atg5-1-N). Among these, 142 DEGs were associated with hypomethylated regions (hypo-DMRs) and were upregulated. This suggests a correlation between DNA demethylation and the ability of Arabidopsis to cope with nitrogen deficiency. Examination of the hypo-DMR-linked upregulated DEGs indicated that the expression of MYB101, an ABA pathway regulator, may be regulated by DNA demethylation and the recruitment of transcription factors (TFs; ERF57, ERF105, ERF48, and ERF111), which may contribute to the growth arrest induced by abscisic acid (ABA). Additionally, we found that DNA methylation might impact the biosynthesis of salicylic acid (SA). The promoter region of ATGH3.12 (PBS3), a key enzyme in SA synthesis, was hypomethylated, combined with overexpression of PBS3 and its potential TF AT3G46070, suggesting that autophagy defects may lead to SA-activated senescence, depending on DNA demethylation. These findings suggest that DNA hypomethylation may impact the mechanism by which Arabidopsis autophagy mutants (atg5-1) respond to nitrogen deficiency, specifically in relation to ABA and SA regulation. Our evaluation of hormone levels verified that these two hormones are significantly enriched under nitrogen deficiency in atg5-1-N compared to Col-0-N.

Identification of mammalian transcription factors that bind to inaccessible chromatin

Transcription factors (TFs) are proteins that affect gene expression by binding to regulatory regions of DNA in a sequence specific manner. The binding of TFs to DNA is controlled by many factors, including the DNA sequence, concentration of TF, chromatin accessibility and co-factors. Here, we systematically investigated the binding mechanism of hundreds of TFs by analysing ChIP-seq data with our explainable statistical model, ChIPanalyser. This tool uses as inputs the DNA sequence binding motif; the capacity to distinguish between strong and weak binding sites; the concentration of TF; and chromatin accessibility. We found that approximately one third of TFs are predicted to bind the genome in a DNA accessibility independent fashion, which includes TFs that can open the chromatin, their co-factors and TFs with similar motifs. Our model predicted this to be the case when the TF binds to its strongest binding regions in the genome, and only a small number of TFs have the capacity to bind dense chromatin at their weakest binding regions, such as CTCF, USF2 and CEBPB. Our study demonstrated that the binding of hundreds of human and mouse TFs is predicted by ChIPanalyser with high accuracy and showed that many TFs can bind dense chromatin.

High temperature influences DNA methylation and transcriptional profiles in sea urchins (Strongylocentrotus intermedius)

BACKGROUND

DNA methylation plays an important role in life processes by affecting gene expression, but it is still unclear how DNA methylation is controlled and how it regulates gene transcription under high temperature stress conditions in Strongylocentrotus intermedius. The potential link between DNA methylation variation and gene expression changes in response to heat stress in S. intermedius was investigated by MethylRAD-seq and RNA-seq analysis. We screened DNA methylation driver genes in order to comprehensively elucidate the regulatory mechanism of its high temperature adaptation at the DNA/RNA level.

RESULTS

The results revealed that high temperature stress significantly affected not only the DNA methylation and transcriptome levels of S. intermedius (P < 0.05), but also growth. MethylRAD-seq analysis revealed 12,129 CG differential methylation sites and 966 CWG differential methylation sites, and identified a total of 189 differentially CG methylated genes and 148 differentially CWG methylated genes. Based on KEGG enrichment analysis, differentially expressed genes (DEGs) are mostly enriched in energy and cell division, immune, and neurological damage pathways. Further RNA-seq analysis identified a total of 1968 DEGs, of which 813 genes were upregulated and 1155 genes were downregulated. Based on the joint MethylRAD-seq and RNA-seq analysis, metabolic processes such as glycosaminoglycan degradation, oxidative phosphorylation, apoptosis, glutathione metabolism, thermogenesis, and lysosomes are regulated by DNA methylation.

CONCLUSIONS

High temperature affected the DNA methylation and expression levels of genes such as MOAP-1, GGT1 and RDH8, which in turn affects the metabolism of HPSE, Cox, glutathione, and retinol, thereby suppressing the immune, energy metabolism, and antioxidant functions of the organism and finally manifesting as stunted growth. In summary, the observations in the present study improve our understanding of the molecular mechanism of the response to high temperature stress in sea urchin.

Promoter hypermethylation of neural-related genes is compatible with stemness in solid cancers

BACKGROUND

DNA hypermethylation is an epigenetic feature that modulates gene expression, and its deregulation is observed in cancer. Previously, we identified a neural-related DNA hypermethylation fingerprint in colon cancer, where most of the top hypermethylated and downregulated genes have known functions in the nervous system. To evaluate the presence of this signature and its relevance to carcinogenesis in general, we considered 16 solid cancer types available in The Cancer Genome Atlas (TCGA).

RESULTS

All tested cancers showed significant enrichment for neural-related genes amongst hypermethylated genes. This signature was already present in two premalignant tissue types and could not be explained by potential confounders such as bivalency status or tumor purity. Further characterization of the neural-related DNA hypermethylation signature in colon cancer showed particular enrichment for genes that are overexpressed during neural differentiation. Lastly, an analysis of upstream regulators identified RE1-Silencing Transcription factor (REST) as a potential mediator of this DNA methylation signature.

CONCLUSION

Our study confirms the presence of a neural-related DNA hypermethylation fingerprint in various cancers, of genes linked to neural differentiation, and points to REST as a possible regulator of this mechanism. We propose that this fingerprint indicates an involvement of DNA hypermethylation in the preservation of neural stemness in cancer cells.

Transcription factor expression repertoire basis for epigenetic and transcriptional subtypes of colorectal cancers

Colorectal cancers (CRCs) form a heterogenous group classified into epigenetic and transcriptional subtypes. The basis for the epigenetic subtypes, exemplified by varying degrees of promoter DNA hypermethylation, and its relation to the transcriptional subtypes is not well understood. We link cancer-specific transcription factor (TF) expression alterations to methylation alterations near TF-binding sites at promoter and enhancer regions in CRCs and their premalignant precursor lesions to provide mechanistic insights into the origins and evolution of the CRC molecular subtypes. A gradient of TF expression changes forms a basis for the subtypes of abnormal DNA methylation, termed CpG-island promoter DNA methylation phenotypes (CIMPs), in CRCs and other cancers. CIMP is tightly correlated with cancer-specific hypermethylation at enhancers, which we term CpG-enhancer methylation phenotype (CEMP). Coordinated promoter and enhancer methylation appears to be driven by downregulation of TFs with common binding sites at the hypermethylated enhancers and promoters. The altered expression of TFs related to hypermethylator subtypes occurs early during CRC development, detectable in premalignant adenomas. TF-based profiling further identifies patients with worse overall survival. Importantly, altered expression of these TFs discriminates the transcriptome-based consensus molecular subtypes (CMS), thus providing a common basis for CIMP and CMS subtypes.

Activation, decommissioning, and dememorization: enhancers in a life cycle

Spatiotemporal regulation of cell type-specific gene expression is essential to convert a zygote into a complex organism that contains hundreds of distinct cell types. A class of cis-regulatory elements called enhancers, which have the potential to enhance target gene transcription, are crucial for precise gene expression programs during development. Following decades of research, many enhancers have been discovered and how enhancers become activated has been extensively studied. However, the mechanisms underlying enhancer silencing are less well understood. We review current understanding of enhancer decommissioning and dememorization, both of which enable enhancer silencing. We highlight recent progress from genome-wide perspectives that have revealed the life cycle of enhancers and how its dynamic regulation underlies cell fate transition, development, cell regeneration, and epigenetic reprogramming.

Toward a comprehensive catalog of regulatory elements

Regulatory elements are the genomic regions that interact with transcription factors to control cell-type-specific gene expression in different cellular environments. A precise and complete catalog of functional elements encoded by the human genome is key to understanding mammalian gene regulation. Here, we review the current state of regulatory element annotation. We first provide an overview of assays for characterizing functional elements, including genome, epigenome, transcriptome, three-dimensional chromatin interaction, and functional validation assays. We then discuss computational methods for defining regulatory elements, including peak-calling and other statistical modeling methods. Finally, we introduce several high-quality lists of regulatory element annotations and suggest potential future directions.

Alterations in the hepatocyte epigenetic landscape in steatosis

Fatty liver disease or the accumulation of fat in the liver, has been reported to affect the global population. This comes with an increased risk for the development of fibrosis, cirrhosis, and hepatocellular carcinoma. Yet, little is known about the effects of a diet containing high fat and alcohol towards epigenetic aging, with respect to changes in transcriptional and epigenomic profiles. In this study, we took up a multi-omics approach and integrated gene expression, methylation signals, and chromatin signals to study the epigenomic effects of a high-fat and alcohol-containing diet on mouse hepatocytes. We identified four relevant gene network clusters that were associated with relevant pathways that promote steatosis. Using a machine learning approach, we predict specific transcription factors that might be responsible to modulate the functionally relevant clusters. Finally, we discover four additional CpG loci and validate aging-related differential CpG methylation. Differential CpG methylation linked to aging showed minimal overlap with altered methylation in steatosis.

Readout of histone methylation by Trim24 locally restricts chromatin opening by p53

The genomic binding sites of the transcription factor (TF) and tumor suppressor p53 are unusually diverse with regard to their chromatin features, including histone modifications, raising the possibility that the local chromatin environment can contextualize p53 regulation. Here, we show that epigenetic characteristics of closed chromatin, such as DNA methylation, do not influence the binding of p53 across the genome. Instead, the ability of p53 to open chromatin and activate its target genes is locally restricted by its cofactor Trim24. Trim24 binds to both p53 and unmethylated histone 3 lysine 4 (H3K4), thereby preferentially localizing to those p53 sites that reside in closed chromatin, whereas it is deterred from accessible chromatin by H3K4 methylation. The presence of Trim24 increases cell viability upon stress and enables p53 to affect gene expression as a function of the local chromatin state. These findings link H3K4 methylation to p53 function and illustrate how specificity in chromatin can be achieved, not by TF-intrinsic sensitivity to histone modifications, but by employing chromatin-sensitive cofactors that locally modulate TF function.

Dissecting the roles of MBD2 isoforms and domains in regulating NuRD complex function during cellular differentiation

The Nucleosome Remodeling and Deacetylation (NuRD) complex is a crucial regulator of cellular differentiation. Two members of the Methyl-CpG-binding domain (MBD) protein family, MBD2 and MBD3, are known to be integral, but mutually exclusive subunits of the NuRD complex. Several MBD2 and MBD3 isoforms are present in mammalian cells, resulting in distinct MBD-NuRD complexes. Whether these different complexes serve distinct functional activities during differentiation is not fully explored. Based on the essential role of MBD3 in lineage commitment, we systematically investigated a diverse set of MBD2 and MBD3 variants for their potential to rescue the differentiation block observed for mouse embryonic stem cells (ESCs) lacking MBD3. While MBD3 is indeed crucial for ESC differentiation to neuronal cells, it functions independently of its MBD domain. We further identify that MBD2 isoforms can replace MBD3 during lineage commitment, however with different potential. Full-length MBD2a only partially rescues the differentiation block, while MBD2b, an isoform lacking an N-terminal GR-rich repeat, fully rescues the Mbd3 KO phenotype. In case of MBD2a, we further show that removing the methylated DNA binding capacity or the GR-rich repeat enables full redundancy to MBD3, highlighting the synergistic requirements for these domains in diversifying NuRD complex function.

Histone exchange sensors reveal variant specific dynamics in mouse embryonic stem cells

Eviction of histones from nucleosomes and their exchange with newly synthesized or alternative variants is a central epigenetic determinant. Here, we define the genome-wide occupancy and exchange pattern of canonical and non-canonical histone variants in mouse embryonic stem cells by genetically encoded exchange sensors. While exchange of all measured variants scales with transcription, we describe variant-specific associations with transcription elongation and Polycomb binding. We found considerable exchange of H3.1 and H2B variants in heterochromatin and repeat elements, contrasting the occupancy and little exchange of H3.3 in these regions. This unexpected association between H3.3 occupancy and exchange of canonical variants is also evident in active promoters and enhancers, and further validated by reduced H3.1 dynamics following depletion of H3.3-specific chaperone, HIRA. Finally, analyzing transgenic mice harboring H3.1 or H3.3 sensors demonstrates the vast potential of this system for studying histone exchange and its impact on gene expression regulation in vivo.