Nucleosome repositioning links DNA (de)methylation and differential CTCF binding during stem cell development

Single-stranded pre-methylated 5mC adapters uncover the methylation profile of plasma ultrashort Single-stranded cell-free DNA

Whole-genome bisulfite sequencing (BS-Seq) measures cytosine methylation changes at single-base resolution and can be used to profile cell-free DNA (cfDNA). In plasma, ultrashort single-stranded cfDNA (uscfDNA, ∼50 nt) has been identified together with 167 bp double-stranded mononucleosomal cell-free DNA (mncfDNA). However, the methylation profile of uscfDNA has not been described. Conventional BS-Seq workflows may not be helpful because bisulfite conversion degrades larger DNA into smaller fragments, leading to erroneous categorization as uscfDNA. We describe the '5mCAdpBS-Seq' workflow in which pre-methylated 5mC (5-methylcytosine) single-stranded adapters are ligated to heat-denatured cfDNA before bisulfite conversion. This method retains only DNA fragments that are unaltered by bisulfite treatment, resulting in less biased uscfDNA methylation analysis. Using 5mCAdpBS-Seq, uscfDNA had lower levels of DNA methylation (∼15%) compared to mncfDNA and was enriched in promoters and CpG islands. Hypomethylated uscfDNA fragments were enriched in upstream transcription start sites (TSSs), and the intensity of enrichment was correlated with expressed genes of hemopoietic cells. Using tissue-of-origin deconvolution, we inferred that uscfDNA is derived primarily from eosinophils, neutrophils, and monocytes. As proof-of-principle, we show that characteristics of the methylation profile of uscfDNA can distinguish non-small cell lung carcinoma from non-cancer samples. The 5mCAdpBS-Seq workflow is recommended for any cfDNA methylation-based investigations.

Nucleosome reorganisation in breast cancer tissues

BACKGROUND

Nucleosome repositioning in cancer is believed to cause many changes in genome organisation and gene expression. Understanding these changes is important to elucidate fundamental aspects of cancer. It is also important for medical diagnostics based on cell-free DNA (cfDNA), which originates from genomic DNA regions protected from digestion by nucleosomes.

RESULTS

We have generated high-resolution nucleosome maps in paired tumour and normal tissues from the same breast cancer patients using MNase-assisted histone H3 ChIP-seq and compared them with the corresponding cfDNA from blood plasma. This analysis has detected single-nucleosome repositioning at key regulatory regions in a patient-specific manner and common cancer-specific patterns across patients. The nucleosomes gained in tumour versus normal tissue were particularly informative of cancer pathways, with ~ 20-fold enrichment at CpG islands, a large fraction of which marked promoters of genes encoding DNA-binding proteins. The tumour tissues were characterised by a 5-10 bp decrease in the average distance between nucleosomes (nucleosome repeat length, NRL), which is qualitatively similar to the differences between pluripotent and differentiated cells. This effect was correlated with gene activity, differential DNA methylation and changes in local occupancy of linker histone variants H1.4 and H1X.

CONCLUSIONS

Our study offers a novel resource of high-resolution nucleosome maps in breast cancer patients and reports for the first time the effect of systematic decrease of NRL in paired tumour versus normal breast tissues from the same patient. Our findings provide a new mechanistic understanding of nucleosome repositioning in tumour tissues that can be valuable for patient diagnostics, stratification and monitoring.

Apparent RNA bridging between PRC2 and chromatin is an artifact of non-specific chromatin precipitation upon RNA degradation

Polycomb repressive complex 2 (PRC2) modifies chromatin to maintain repression of genes specific for other cell lineages. In vitro, RNA inhibits PRC2 activity, but the effect of RNA on PRC2 in cells is less clear, with studies concluding that RNA either antagonizes or promotes PRC2 chromatin association. The addition of RNase A to chromatin immunoprecipitation reactions has been reported to reduce detection of PRC2 target sites, suggesting the existence of RNA bridges connecting PRC2 to chromatin. Here, we show that the apparent loss of PRC2 chromatin association after RNase A treatment is due to non-specific chromatin precipitation. RNA degradation precipitates chromatin out of solution, thereby masking enrichment of specific DNA sequences in chromatin immunoprecipitation reactions. Maintaining chromatin solubility by the addition of poly-L-glutamic acid rescues detection of PRC2 chromatin occupancy upon RNA degradation. These findings undermine support for the model that RNA bridges PRC2 and chromatin in cells.

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.

Detection of new pioneer transcription factors as cell-type-specific nucleosome binders

Wrapping of DNA into nucleosomes restricts accessibility to DNA and may affect the recognition of binding motifs by transcription factors. A certain class of transcription factors, the pioneer transcription factors, can specifically recognize their DNA binding sites on nucleosomes, initiate local chromatin opening, and facilitate the binding of co-factors in a cell-type-specific manner. For the majority of human pioneer transcription factors, the locations of their binding sites, mechanisms of binding, and regulation remain unknown. We have developed a computational method to predict the cell-type-specific ability of transcription factors to bind nucleosomes by integrating ChIP-seq, MNase-seq, and DNase-seq data with details of nucleosome structure. We have demonstrated the ability of our approach in discriminating pioneer from canonical transcription factors and predicted new potential pioneer transcription factors in H1, K562, HepG2, and HeLa-S3 cell lines. Last, we systematically analyzed the interaction modes between various pioneer transcription factors and detected several clusters of distinctive binding sites on nucleosomal DNA.

Detection of new pioneer transcription factors as cell-type specific nucleosome binders

Wrapping of DNA into nucleosomes restricts accessibility to the DNA and may affect the recognition of binding motifs by transcription factors. A certain class of transcription factors, the pioneer transcription factors, can specifically recognize their DNA binding sites on nucleosomes, may initiate local chromatin opening and facilitate the binding of co-factors in a cell-type-specific manner. For the majority of human pioneer transcription factors, the locations of their binding sites, mechanisms of binding and regulation remain unknown. We have developed a computational method to predict the cell-type-specific ability of transcription factors to bind nucleosomes by integrating ChIP-seq, MNase-seq and DNase-seq data with details of nucleosome structure. We have demonstrated the ability of our approach in discriminating pioneer from canonical transcription factors and predicted new potential pioneer transcription factors in H1, K562, HepG2 and HeLa cell lines. Lastly, we systemically analyzed the interaction modes between various pioneer transcription factors and detected several clusters of distinctive binding sites on nucleosomal DNA.

Nucleosome repositioning in chronic lymphocytic leukemia

The location of nucleosomes in the human genome determines the primary chromatin structure and regulates access to regulatory regions. However, genome-wide information on deregulated nucleosome occupancy and its implications in primary cancer cells is scarce. Here, we conducted a genome-wide comparison of high-resolution nucleosome maps in peripheral blood B cells from patients with chronic lymphocytic leukemia (CLL) and healthy individuals at single-base-pair resolution. Our investigation uncovered significant changes of nucleosome positioning in CLL. Globally, the spacing between nucleosomes-the nucleosome repeat length (NRL)-is shortened in CLL. This effect is stronger in the more aggressive IGHV-unmutated CLL subtype than in the IGHV-mutated CLL subtype. Changes in nucleosome occupancy at specific sites are linked to active chromatin remodeling and reduced DNA methylation. Nucleosomes lost or gained in CLL marks differential binding of 3D chromatin organizers such as CTCF as well as immune response-related transcription factors and delineated mechanisms of epigenetic deregulation. The principal component analysis of nucleosome occupancy in cancer-specific regions allowed the classification of samples between cancer subtypes and normal controls. Furthermore, patients could be better assigned to CLL subtypes according to differential nucleosome occupancy than based on DNA methylation or gene expression. Thus, nucleosome positioning constitutes a novel readout to dissect molecular mechanisms of disease progression and to stratify patients. Furthermore, we anticipate that the global nucleosome repositioning detected in our study, such as changes in the NRL, can be exploited for liquid biopsy applications based on cell-free DNA to stratify patients and monitor disease progression.

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.

DNA strand asymmetry generated by CpG hemimethylation has opposing effects on CTCF binding

CpG methylation generally occurs on both DNA strands and is essential for mammalian development and differentiation. Until recently, hemimethylation, in which only one strand is methylated, was considered to be simply a transitory state generated during DNA synthesis. The discovery that a subset of CCCTC-binding factor (CTCF) binding sites is heritably hemimethylated suggests that hemimethylation might have an unknown biological function. Here we show that the binding of CTCF is profoundly altered by which DNA strand is methylated and by the specific CTCF binding motif. CpG methylation on the motif strand can inhibit CTCF binding by up to 7-fold, whereas methylation on the opposite strand can stimulate binding by up to 4-fold. Thus, hemimethylation can alter binding by up to 28-fold in a strand-specific manner. The mechanism for sensing methylation on the opposite strand requires two critical residues, V454 and S364, within CTCF zinc fingers 7 and 4. Similar to methylation, CpG hydroxymethylation on the motif strand can inhibit CTCF binding by up to 4-fold. However, hydroxymethylation on the opposite strand removes the stimulatory effect. Strand-specific methylation states may therefore provide a mechanism to explain the transient and dynamic nature of CTCF-mediated chromatin interactions.

Paternal Combined Botanicals Contribute to the Prevention of Estrogen Receptor-Negative Mammary Cancer in Transgenic Mice

BACKGROUND

Parental nutritional interventions have considerably affected gametogenesis and embryogenesis, leading to the differential susceptibility of offspring to chronic diseases such as cancer. Moreover, combinatorial bioactive diets are more efficacious in ameliorating epigenetic aberrations in tumorigenesis.

OBJECTIVES

We sought to investigate the transgenerational influence and epigenetic regulation of paternal sulforaphane (SFN)-rich broccoli sprouts (BSp) and epigallocatechin-3-gallate (EGCG)-rich green tea polyphenols (GTPs) consumption in the prevention of estrogen receptor-negative [ER(-)] mammary cancer in transgenic mice.

METHODS

Human breast cancer cells were used to detect cell viability and epigenetic-related gene expression after treatment with EGCG and/or SFN. Twenty-four C3 or HER2/neu males were randomly assigned into 4 groups and treated with control, 26% BSp (w/w) in food, 0.5% GTPs (w/v) in drinking water or combined BSp and GTPs for 7 wk before mating. Tumor growth of nontreated female pups was monitored weekly for 19 wk (C3) and 25 wk (HER2/neu). Tumor- and epigenetic-related protein expression and enzyme activities in mammary tumors were measured. Sperms were isolated from treated males for RNA sequencing and reduced-representation bisulfite sequencing analysis. Data were analyzed with a 2-factor or 3-factor analysis of variance.

RESULTS

EGCG and SFN inhibited breast cancer cell growth via epigenetic regulation. Combined BSp and GTPs synergistically (combination index < 1) suppressed tumor growth over time (P < 0.001) in 2 mouse models. Key tumor-related proteins were found differentially expressed (P < 0.05) along with epigenetic regulations in offspring mammary tumors. The transcriptome profile of sperm derived from dietary-treated males revealed differentially expressed genes correlated with spermatogenesis and breast cancer progression. DNA methylomes of the sperm and further integrated analysis with transcriptomes indicate that DNA methylation alone may not contribute to sufficient regulation in dietary-treated sperm pronucleus, leading to offspring tumor suppression.

CONCLUSIONS

Collectively, paternal consumption of combined BSp and GTPs shows potential for preventing ER(-) mammary cancer through transgenerational effects. J Nutr 2023;xx:xx-xx.

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.

CTCF and Its Multi-Partner Network for Chromatin Regulation

Architectural proteins are essential epigenetic regulators that play a critical role in organizing chromatin and controlling gene expression. CTCF (CCCTC-binding factor) is a key architectural protein responsible for maintaining the intricate 3D structure of chromatin. Because of its multivalent properties and plasticity to bind various sequences, CTCF is similar to a Swiss knife for genome organization. Despite the importance of this protein, its mechanisms of action are not fully elucidated. It has been hypothesized that its versatility is achieved through interaction with multiple partners, forming a complex network that regulates chromatin folding within the nucleus. In this review, we delve into CTCF's interactions with other molecules involved in epigenetic processes, particularly histone and DNA demethylases, as well as several long non-coding RNAs (lncRNAs) that are able to recruit CTCF. Our review highlights the importance of CTCF partners to shed light on chromatin regulation and pave the way for future exploration of the mechanisms that enable the finely-tuned role of CTCF as a master regulator of chromatin.

Epigenetic signals that direct cell type-specific interferon beta response in mouse cells

The antiviral response induced by type I interferon (IFN) via the JAK-STAT signaling cascade activates hundreds of IFN-stimulated genes (ISGs) across human and mouse tissues but varies between cell types. However, the links between the underlying epigenetic features and the ISG profile are not well understood. We mapped ISGs, binding sites of the STAT1 and STAT2 transcription factors, chromatin accessibility, and histone H3 lysine modification by acetylation (ac) and mono-/tri-methylation (me1, me3) in mouse embryonic stem cells and fibroblasts before and after IFNβ treatment. A large fraction of ISGs and STAT-binding sites was cell type specific with promoter binding of a STAT1/2 complex being a key driver of ISGs. Furthermore, STAT1/2 binding to putative enhancers induced ISGs as inferred from a chromatin co-accessibility analysis. STAT1/2 binding was dependent on the chromatin context and positively correlated with preexisting H3K4me1 and H3K27ac marks in an open chromatin state, whereas the presence of H3K27me3 had an inhibitory effect. Thus, chromatin features present before stimulation represent an additional regulatory layer for the cell type-specific antiviral response.

Diverse silent chromatin states modulate genome compartmentalization and loop extrusion barriers

The relationships between chromosomal compartmentalization, chromatin state and function are poorly understood. Here by profiling long-range contact frequencies in HCT116 colon cancer cells, we distinguish three silent chromatin states, comprising two types of heterochromatin and a state enriched for H3K9me2 and H2A.Z that exhibits neutral three-dimensional interaction preferences and which, to our knowledge, has not previously been characterized. We find that heterochromatin marked by H3K9me3, HP1α and HP1β correlates with strong compartmentalization. We demonstrate that disruption of DNA methyltransferase activity greatly remodels genome compartmentalization whereby domains lose H3K9me3-HP1α/β binding and acquire the neutrally interacting state while retaining late replication timing. Furthermore, we show that H3K9me3-HP1α/β heterochromatin is permissive to loop extrusion by cohesin but refractory to CTCF binding. Together, our work reveals a dynamic structural and organizational diversity of the silent portion of the genome and establishes connections between the regulation of chromatin state and chromosome organization, including an interplay between DNA methylation, compartmentalization and loop extrusion.

The use of mass spectrometry in a proteome-centered multiomics study of human pituitary adenomas

A pituitary adenoma (PA) is a common intracranial neoplasm, and is a complex, chronic, and whole-body disease with multicausing factors, multiprocesses, and multiconsequences. It is very difficult to clarify molecular mechanism and treat PAs from the single-factor strategy model. The rapid development of multiomics and systems biology changed the paradigms from a traditional single-factor strategy to a multiparameter systematic strategy for effective management of PAs. A series of molecular alterations at the genome, transcriptome, proteome, peptidome, metabolome, and radiome levels are involved in pituitary tumorigenesis, and mutually associate into a complex molecular network system. Also, the center of multiomics is moving from structural genomics to phenomics, including proteomics and metabolomics in the medical sciences. Mass spectrometry (MS) has been extensively used in phenomics studies of human PAs to clarify molecular mechanisms, and to discover biomarkers and therapeutic targets/drugs. MS-based proteomics and proteoform studies play central roles in the multiomics strategy of PAs. This article reviews the status of multiomics, multiomics-based molecular pathway networks, molecular pathway network-based pattern biomarkers and therapeutic targets/drugs, and future perspectives for personalized, predeictive, and preventive (3P) medicine in PAs.

Mechanisms that regulate the activities of TET proteins

The ten-eleven translocation (TET) family of dioxygenases consists of three members, TET1, TET2, and TET3. All three TET enzymes have Fe⁺² and α-ketoglutarate (α-KG)-dependent dioxygenase activities, catalyzing the 1st step of DNA demethylation by converting 5-methylcytosine (5mC) to 5-hydroxymethylcytosine (5hmC), and further oxidize 5hmC to 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC). Gene knockout studies demonstrated that all three TET proteins are involved in the regulation of fetal organ generation during embryonic development and normal tissue generation postnatally. TET proteins play such roles by regulating the expression of key differentiation and fate-determining genes via (1) enzymatic activity-dependent DNA methylation of the promoters and enhancers of target genes; and (2) enzymatic activity-independent regulation of histone modification. Interacting partner proteins and post-translational regulatory mechanisms regulate the activities of TET proteins. Mutations and dysregulation of TET proteins are involved in the pathogenesis of human diseases, specifically cancers. Here, we summarize the research on the interaction partners and post-translational modifications of TET proteins. We also discuss the molecular mechanisms by which these partner proteins and modifications regulate TET functioning and target gene expression. Such information will help in the design of medications useful for targeted therapy of TET-mutant-related diseases.

Enhanced nucleosome assembly at CpG sites containing an extended 5-methylcytosine analogue

Methylation of cytosine to 5-methylcytosine (mC) at CpG sites is a prevalent reversible epigenetic mark in vertebrates established by DNA methyltransferases (MTases); the attached methyl groups can alter local structure of DNA and chromatin as well as binding of dedicated proteins. Nucleosome assembly on methylated DNA has been studied extensively, however little is known how the chromatin structure is affected by larger chemical variations in the major groove of DNA. Here, we studied the nucleosome formation in vitro on DNA containing an extended 5mC analog, 5-(6-azidohex-2-ynyl)cytosine (ahyC) installed at biological relevant CpG sites. We found that multiple ahyC residues on 80-Widom and Hsp70 promoter DNA fragments proved compatible with nucleosome assembly. Moreover, unlike mC, ahyC increases the affinity of histones to the DNA, partially altering nucleosome positioning, stability, and the action of chromatin remodelers. Based on molecular dynamics calculations, we suggest that these new features are due to increased DNA flexibility at ahyC-modified sites. Our findings provide new insights into the biophysical behavior of modified DNA and open new ways for directed design of synthetic nucleosomes.

DNA sequence-dependent formation of heterochromatin nanodomains

The mammalian epigenome contains thousands of heterochromatin nanodomains (HNDs) marked by di- and trimethylation of histone H3 at lysine 9 (H3K9me2/3), which have a typical size of 3–10 nucleosomes. However, what governs HND location and extension is only partly understood. Here, we address this issue by introducing the chromatin hierarchical lattice framework (ChromHL) that predicts chromatin state patterns with single-nucleotide resolution. ChromHL is applied to analyse four HND types in mouse embryonic stem cells that are defined by histone methylases SUV39H1/2 or GLP, transcription factor ADNP or chromatin remodeller ATRX. We find that HND patterns can be computed from PAX3/9, ADNP and LINE1 sequence motifs as nucleation sites and boundaries that are determined by DNA sequence (e.g. CTCF binding sites), cooperative interactions between nucleosomes as well as nucleosome-HP1 interactions. Thus, ChromHL rationalizes how patterns of H3K9me2/3 are established and changed via the activity of protein factors in processes like cell differentiation.

The ability to predict epigenetic regulation is an important challenge in biology. Here the authors describe heterochromatin nanodomains (HNDs) and compare four different types of H3K9me2/3-marked HNDs in mouse embryonic stem cells. They further develop a computational framework to predict genome-wide HND maps from DNA sequence and protein concentrations, at single-nucleotide resolution.

DNA methylation cues in nucleosome geometry, stability and unwrapping

Cytosine methylation at the 5-carbon position is an essential DNA epigenetic mark in many eukaryotic organisms. Although countless structural and functional studies of cytosine methylation have been reported, our understanding of how it influences the nucleosome assembly, structure, and dynamics remains obscure. Here, we investigate the effects of cytosine methylation at CpG sites on nucleosome dynamics and stability. By applying long molecular dynamics simulations on several microsecond time scale, we generate extensive atomistic conformational ensembles of full nucleosomes. Our results reveal that methylation induces pronounced changes in geometry for both linker and nucleosomal DNA, leading to a more curved, under-twisted DNA, narrowing the adjacent minor grooves, and shifting the population equilibrium of sugar-phosphate backbone geometry. These DNA conformational changes are associated with a considerable enhancement of interactions between methylated DNA and the histone octamer, doubling the number of contacts at some key arginines. H2A and H3 tails play important roles in these interactions, especially for DNA methylated nucleosomes. This, in turn, prevents a spontaneous DNA unwrapping of 3-4 helical turns for the methylated nucleosome with truncated histone tails, otherwise observed in the unmethylated system on several microseconds time scale.

NucPosDB: a database of nucleosome positioning in vivo and nucleosomics of cell-free DNA

Nucleosome positioning is involved in many gene regulatory processes happening in the cell, and it may change as cells differentiate or respond to the changing microenvironment in a healthy or diseased organism. One important implication of nucleosome positioning in clinical epigenetics is its use in the “nucleosomics” analysis of cell-free DNA (cfDNA) for the purpose of patient diagnostics in liquid biopsies. The rationale for this is that the apoptotic nucleases that digest chromatin of the dying cells mostly cut DNA between nucleosomes. Thus, the short pieces of DNA in body fluids reflect the positions of nucleosomes in the cells of origin. Here, we report a systematic nucleosomics database — NucPosDB — curating published nucleosome positioning datasets in vivo as well as datasets of sequenced cell-free DNA (cfDNA) that reflect nucleosome positioning in situ in the cells of origin. Users can select subsets of the database by a number of criteria and then obtain raw or processed data. NucPosDB also reports the originally determined regions with stable nucleosome occupancy across several individuals with a given condition. An additional section provides a catalogue of computational tools for the analysis of nucleosome positioning or cfDNA experiments and theoretical algorithms for the prediction of nucleosome positioning preferences from DNA sequence. We provide an overview of the field, describe the structure of the database in this context, and demonstrate data variability using examples of different medical conditions. NucPosDB is useful both for the analysis of fundamental gene regulation processes and the training of computational models for patient diagnostics based on cfDNA. The database currently curates ~ 400 publications on nucleosome positioning in cell lines and in situ as well as cfDNA from > 10,000 patients and healthy volunteers. For open-access cfDNA datasets as well as key MNase-seq datasets in human cells, NucPosDB allows downloading processed mapped data in addition to the regions with stable nucleosome occupancy. NucPosDB is available at https://generegulation.org/nucposdb/.

Advances in measuring DNA methylation

DNA methylation is one of the most important components of epigenetics, which plays essential roles in maintaining genome stability and regulating gene expression. In recent years, DNA methylation measuring methods have been continuously optimized. Combined with next generation sequencing technologies, these approaches have enabled the detection of genome-wide cytosine methylation at single-base resolution. In this paper, we review the development of 5-methylcytosine and its oxidized derivatives measuring methods, and recent advancement of single-cell epigenome sequencing technologies, offering more referable information for the selection and optimization of DNA methylation sequencing technologies and related research.

DNA 5-hydroxymethylcytosine in pediatric central nervous system tumors may impact tumor classification and is a positive prognostic marker

Background

Nucleotide-specific 5-hydroxymethylcytosine (5hmC) remains understudied in pediatric central nervous system (CNS) tumors. 5hmC is abundant in the brain, and alterations to 5hmC in adult CNS tumors have been reported. However, traditional approaches to measure DNA methylation do not distinguish between 5-methylcytosine (5mC) and its oxidized counterpart 5hmC, including those used to build CNS tumor DNA methylation classification systems. We measured 5hmC and 5mC epigenome-wide at nucleotide resolution in glioma, ependymoma, and embryonal tumors from children, as well as control pediatric brain tissues using tandem bisulfite and oxidative bisulfite treatments followed by hybridization to the Illumina Methylation EPIC Array that interrogates over 860,000 CpG loci.

Results

Linear mixed effects models adjusted for age and sex tested the CpG-specific differences in 5hmC between tumor and non-tumor samples, as well as between tumor subtypes. Results from model-based clustering of tumors was used to test the relation of cluster membership with patient survival through multivariable Cox proportional hazards regression. We also assessed the robustness of multiple epigenetic CNS tumor classification methods to 5mC-specific data in both pediatric and adult CNS tumors. Compared to non-tumor samples, tumors were hypohydroxymethylated across the epigenome and tumor 5hmC localized to regulatory elements crucial to cell identity, including transcription factor binding sites and super-enhancers. Differentially hydroxymethylated loci among tumor subtypes tended to be hypermethylated and disproportionally found in CTCF binding sites and genes related to posttranscriptional RNA regulation, such as DICER1. Model-based clustering results indicated that patients with low 5hmC patterns have poorer overall survival and increased risk of recurrence. Our results suggest 5mC-specific data from OxBS-treated samples impacts methylation-based tumor classification systems giving new opportunities for further refinement of classifiers for both pediatric and adult tumors.

Conclusions

We identified that 5hmC localizes to super-enhancers, and genes commonly implicated in pediatric CNS tumors were differentially hypohydroxymethylated. We demonstrated that distinguishing methylation and hydroxymethylation is critical in identifying tumor-related epigenetic changes. These results have implications for patient prognostication, considerations of epigenetic therapy in CNS tumors, and for emerging molecular neuropathology classification approaches.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13148-021-01156-9.

Analyses of Promoter , Enhancer, and Nucleosome Organization in Mammalian Cells by MNase-Seq

MNase-Seq is a genome-wide procedure that allows mapping of DNA associated to nucleosomes following micrococcal nuclease digestion. It is a rapid and robust technology useful for the analysis of chromatin properties genome-wide at the resolution of mono-nucleosomes. Here, we describe how to produce high-resolution nucleosome maps of cells grown in suspension or adherent mammalian cells. After only three steps: nuclei or cell preparation, native MNase digestion and DNA purification, libraries for high-throughput sequencing can be prepared. Genome-wide nucleosome maps allow analyzing chromatin opening at promoters or enhancers, nucleosome displacement, or labile nucleosome occupancy depending on the digestion condition used. As presented, MNase-Seq is a versatile tool for investigating chromatin dynamics, regulation, and to define open chromatin regions of regulatory elements in mammalian genomes.

DNA G-Quadruplexes Contribute to CTCF Recruitment

G-quadruplex (G4) sites in the human genome frequently colocalize with CCCTC-binding factor (CTCF)-bound sites in CpG islands (CGIs). We aimed to clarify the role of G4s in CTCF positioning. Molecular modeling data suggested direct interactions, so we performed in vitro binding assays with quadruplex-forming sequences from CGIs in the human genome. G4s bound CTCF with Kd values similar to that of the control duplex, while respective i-motifs exhibited no affinity for CTCF. Using ChIP-qPCR assays, we showed that G4-stabilizing ligands enhance CTCF occupancy at a G4-prone site in STAT3 gene. In view of the reportedly increased CTCF affinity for hypomethylated DNA, we next questioned whether G4s also facilitate CTCF recruitment to CGIs via protecting CpG sites from methylation. Bioinformatics analysis of previously published data argued against such a possibility. Finally, we questioned whether G4s facilitate CTCF recruitment by affecting chromatin structure. We showed that three architectural chromatin proteins of the high mobility group colocalize with G4s in the genome and recognize parallel-stranded or mixed-topology G4s in vitro. One of such proteins, HMGN3, contributes to the association between G4s and CTCF according to our bioinformatics analysis. These findings support both direct and indirect roles of G4s in CTCF recruitment.

The Roles of DNA Demethylases in Triple-Negative Breast Cancer

Triple-negative breast cancers (TNBCs) are very heterogenous, molecularly diverse, and are characterized by a high propensity to relapse or metastasize. Clinically, TNBC remains a diagnosis of exclusion by the lack of hormone receptors (Estrogen Receptor (ER) and Progesterone Receptor (PR)) as well as the absence of overexpression and/or amplification of HER2. DNA methylation plays an important role in breast cancer carcinogenesis and TNBCs have a distinct DNA methylation profile characterized by marked hypomethylation and lower gains of methylations compared to all other subtypes. DNA methylation is regulated by the balance of DNA methylases (DNMTs) and DNA demethylases (TETs). Here, we review the roles of TETs as context-dependent tumor-suppressor genes and/or oncogenes in solid tumors, and we discuss the current understandings of the oncogenic role of TET1 and its therapeutic implications in TNBCs.

NUCOME: A comprehensive database of nucleosome organization referenced landscapes in mammalian genomes

Background

Nucleosome organization is involved in many regulatory activities in various organisms. However, studies integrating nucleosome organization in mammalian genomes are very limited mainly due to the lack of comprehensive data quality control (QC) assessment and uneven data quality of public data sets.

Results

The NUCOME is a database focused on filtering qualified nucleosome organization referenced landscapes covering various cell types in human and mouse based on QC metrics. The filtering strategy guarantees the quality of nucleosome organization referenced landscapes and exempts users from redundant data set selection and processing. The NUCOME database provides standardized, qualified data source and informative nucleosome organization features at a whole-genome scale and on the level of individual loci.

Conclusions

The NUCOME provides valuable data resources for integrative analyses focus on nucleosome organization. The NUCOME is freely available at http://compbio-zhanglab.org/NUCOME.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12859-021-04239-9.

RSC primes the quiescent genome for hypertranscription upon cell-cycle re-entry

Quiescence is a reversible G0 state essential for differentiation, regeneration, stem-cell renewal, and immune cell activation. Necessary for long-term survival, quiescent chromatin is compact, hypoacetylated, and transcriptionally inactive. How transcription activates upon cell-cycle re-entry is undefined. Here we report robust, widespread transcription within the first minutes of quiescence exit. During quiescence, the chromatin-remodeling enzyme RSC was already bound to the genes induced upon quiescence exit. RSC depletion caused severe quiescence exit defects: a global decrease in RNA polymerase II (Pol II) loading, Pol II accumulation at transcription start sites, initiation from ectopic upstream loci, and aberrant antisense transcription. These phenomena were due to a combination of highly robust Pol II transcription and severe chromatin defects in the promoter regions and gene bodies. Together, these results uncovered multiple mechanisms by which RSC facilitates initiation and maintenance of large-scale, rapid gene expression despite a globally repressive chromatin state.

Effects of DNA Methylation on TFs in Human Embryonic Stem Cells

DNA methylation is an important epigenetic mechanism for gene regulation. The conventional view of DNA methylation is that DNA methylation could disrupt protein-DNA interactions and repress gene expression. Several recent studies reported that DNA methylation could alter transcription factors (TFs) binding sequence specificity in vitro. Here, we took advantage of the large sets of ChIP-seq data for TFs and whole-genome bisulfite sequencing data in many cell types to perform a systematic analysis of the protein-DNA methylation in vivo. We observed that many TFs could bind methylated DNA regions, especially in H1-hESC cells. By locating binding sites, we confirmed that some TFs could bind to methylated CpGs directly. The different proportion of CpGs at TF binding specificity motifs in different methylation statuses shows that some TFs are sensitive to methylation and some could bind to the methylated DNA with different motifs, such as CEBPB and CTCF. At the same time, TF binding could interactively alter local DNA methylation. The TF hypermethylation binding sites extensively overlap with enhancers. And we also found that some DNase I hypersensitive sites were specifically hypermethylated in H1-hESC cells. At last, compared with TFs' binding regions in multiple cell types, we observed that CTCF binding to high methylated regions in H1-hESC were not conservative. These pieces of evidence indicate that TFs that bind to hypermethylation DNA in H1-hESC cells may associate with enhancers to regulate special biological functions.

Genomic properties of variably methylated retrotransposons in mouse

BACKGROUND

Transposable elements (TEs) are enriched in cytosine methylation, preventing their mobility within the genome. We previously identified a genome-wide repertoire of candidate intracisternal A particle (IAP) TEs in mice that exhibit inter-individual variability in this methylation (VM-IAPs) with implications for genome function.

RESULTS

Here we validate these metastable epialleles and discover a novel class that exhibit tissue specificity (tsVM-IAPs) in addition to those with uniform methylation in all tissues (constitutive- or cVM-IAPs); both types have the potential to regulate genes in cis. Screening for variable methylation at other TEs shows that this phenomenon is largely limited to IAPs, which are amongst the youngest and most active endogenous retroviruses. We identify sequences enriched within cVM-IAPs, but determine that these are not sufficient to confer epigenetic variability. CTCF is enriched at VM-IAPs with binding inversely correlated with DNA methylation. We uncover dynamic physical interactions between cVM-IAPs with low methylation ranges and other genomic loci, suggesting that VM-IAPs have the potential for long-range regulation.

CONCLUSION

Our findings indicate that a recently evolved interplay between genetic sequence, CTCF binding, and DNA methylation at young TEs can result in inter-individual variability in transcriptional outcomes with implications for phenotypic variation.

Hippocampal epigenetic and insulin-like growth factor alterations in noninvasive versus invasive mechanical ventilation in preterm lambs

BACKGROUND

The brain of chronically ventilated preterm human infants is vulnerable to collateral damage during invasive mechanical ventilation (IMV). Damage is manifest, in part, by learning and memory impairments, which are hippocampal functions. A molecular regulator of hippocampal development is insulin-like growth factor 1 (IGF1). A gentler ventilation strategy is noninvasive respiratory support (NRS). We tested the hypotheses that NRS leads to greater levels of IGF1 messenger RNA (mRNA) variants and distinct epigenetic profile along the IGF1 gene locus in the hippocampus compared to IMV.

METHODS

Preterm lambs were managed by NRS or IMV for 3 or 21 days. Isolated hippocampi were analyzed for IGF1 mRNA levels and splice variants for promoter 1 (P1), P2, and IGF1A and 1B, DNA methylation in P1 region, and histone covalent modifications along the gene locus.

RESULTS

NRS had significantly greater levels of IGF1 P1 (predominant transcript), and 1A and 1B mRNA variants compared to IMV at 3 or 21 days. NRS also led to more DNA methylation and greater occupancy of activating mark H3K4 trimethylation (H3K4me³), repressive mark H3K27me³, and elongation mark H3K36me³ compared to IMV.

CONCLUSIONS

NRS leads to distinct IGF1 mRNA variant levels and epigenetic profile in the hippocampus compared to IMV.

IMPACT

Our study shows that 3 or 21 days of NRS of preterm lambs leads to distinct IGF1 mRNA variant levels and epigenetic profile in the hippocampus compared to IMV. Preterm infant studies suggest that NRS leads to better neurodevelopmental outcomes later in life versus IMV. Also, duration of IMV is directly related to hippocampal damage; however, molecular players remain unknown. NRS, as a gentler mode of respiratory management of preterm neonates, may reduce damage to the immature hippocampus through an epigenetic mechanism.

Epigenetics of obesity

The pathophysiology of obesity is complex and includes changes in eating behavior, genetic, epigenetic, environmental factors, and much more. To date, ~40 genetic polymorphisms are associated with obesity and fat distribution. However, since these options do not fully explain the inheritance of obesity, other options, such as epigenetic changes, need to be considered. Epigenetic modifications affect gene expression without changing the deoxyribonucleic acid sequence. In addition, environmental exposure during critical periods of development can affect the epigenetic tags and lead to obesity. A deeper understanding of the epigenetic mechanisms underlying obesity can aid in prevention based on lifestyle changes. This review focuses on the role of epigenetic modifications in the development of obesity and related conditions.

CTCF as a boundary factor for cohesin-mediated loop extrusion: evidence for a multi-step mechanism

Mammalian genome structure is closely linked to function. At the scale of kilobases to megabases, CTCF and cohesin organize the genome into chromatin loops. Mechanistically, cohesin is proposed to extrude chromatin loops bidirectionally until it encounters occupied CTCF DNA-binding sites. Curiously, loops form predominantly between CTCF binding sites in a convergent orientation. How CTCF interacts with and blocks cohesin extrusion in an orientation-specific manner has remained a mechanistic mystery. Here, we review recent papers that have shed light on these processes and suggest a multi-step interaction between CTCF and cohesin. This interaction may first involve a pausing step, where CTCF halts cohesin extrusion, followed by a stabilization step of the CTCF-cohesin complex, resulting in a chromatin loop. Finally, we discuss our own recent studies on an internal RNA-Binding Region (RBRi) in CTCF to elucidate its role in regulating CTCF clustering, target search mechanisms and chromatin loop formation and future challenges.

Aging of preleukemic thymocytes drives CpG island hypermethylation in T-cell acute lymphoblastic leukemia

Cancer cells display DNA hypermethylation at specific CpG islands in comparison to their normal healthy counterparts, but the mechanism that drives this so-called CpG island methylator phenotype (CIMP) remains poorly understood. Here, we show that CpG island methylation in human T-cell acute lymphoblastic leukemia (T-ALL) mainly occurs at promoters of Polycomb Repressor Complex 2 (PRC2) target genes that are not expressed in normal or malignant T-cells and which display a reciprocal association with H3K27me3 binding. In addition, we revealed that this aberrant methylation profile reflects the epigenetic history of T-ALL and is established already in pre-leukemic, self-renewing thymocytes that precede T-ALL development. Finally, we unexpectedly uncover that this age-related CpG island hypermethylation signature in T-ALL is completely resistant to the FDA-approved hypomethylating agent Decitabine. Altogether, we here provide conceptual evidence for the involvement of a pre-leukemic phase characterized by self-renewing thymocytes in the pathogenesis of human T-ALL.

DNA methylation methods: Global DNA methylation and methylomic analyses

DNA methylation provides a pivotal layer of epigenetic regulation in eukaryotes that has significant involvement for numerous biological processes in health and disease. The function of methylation of cytosine bases in DNA was originally proposed as a "silencing" epigenetic marker and focused on promoter regions of genes for decades. Improved technologies and accumulating studies have been extending our understanding of the roles of DNA methylation to various genomic contexts including gene bodies, repeat sequences and transcriptional start sites. The demand for comprehensively describing DNA methylation patterns spawns a diversity of DNA methylation profiling technologies that target its genomic distribution. These approaches have enabled the measurement of cytosine methylation from specific loci at restricted regions to single-base-pair resolution on a genome-scale level. In this review, we discuss the different DNA methylation analysis technologies primarily based on the initial treatments of DNA samples: bisulfite conversion, endonuclease digestion and affinity enrichment, involving methodology evolution, principles, applications, and their relative merits. This review may offer referable information for the selection of various platforms for genome-wide analysis of DNA methylation.

Linker histone epitopes are hidden by in situ higher-order chromatin structure

BACKGROUND

Histone H1 is the most mobile histone in the cell nucleus. Defining the positions of H1 on chromatin in situ, therefore, represents a challenge. Immunoprecipitation of formaldehyde-fixed and sonicated chromatin, followed by DNA sequencing (xChIP-seq), is traditionally the method for mapping histones onto DNA elements. But since sonication fragmentation precedes ChIP, there is a consequent loss of information about chromatin higher-order structure. Here, we present a new method, xxChIP-seq, employing antibody binding to fixed intact in situ chromatin, followed by extensive washing, a second fixation, sonication and immunoprecipitation. The second fixation is intended to prevent the loss of specifically bound antibody during washing and subsequent sonication and to prevent antibody shifting to epitopes revealed by the sonication process. In many respects, xxChIP-seq is comparable to immunostaining microscopy, which also involves interaction of the primary antibody with fixed and permeabilized intact cells. The only epitopes displayed after immunostaining are the "exposed" epitopes, not "hidden" by the fixation of chromatin higher-order structure. Comparison of immunoprecipitated fragments between xChIP-seq versus xxChIP-seq should indicate which epitopes become inaccessible with fixation and identify their associated DNA elements.

RESULTS

We determined the genomic distribution of histone variants H1.2 and H1.5 in human myeloid leukemia cells HL-60/S4 and compared their epitope exposure by both xChIP-seq and xxChIP-seq, as well as high-resolution microscopy, illustrating the influences of preserved chromatin higher-order structure in situ. We found that xChIP and xxChIP H1 signals are in general negatively correlated, with differences being more pronounced near active regulatory regions. Among the intriguing observations, we find that transcription-related regions and histone PTMs (i.e., enhancers, promoters, CpG islands, H3K4me1, H3K4me3, H3K9ac, H3K27ac and H3K36me3) exhibit significant deficiencies (depletions) in H1.2 and H1.5 xxChIP-seq reads, compared to xChIP-seq. These observations suggest the existence of in situ transcription-related chromatin higher-order structures stabilized by formaldehyde.

CONCLUSION

Comparison of H1 xxChIP-seq to H1 xChIP-seq allows the development of hypotheses on the chromosomal localization of (stabilized) higher-order structure, indicated by the generation of "hidden" H1 epitopes following formaldehyde crosslinking. Changes in H1 epitope exposure surrounding averaged chromosomal binding sites or epigenetic modifications can also indicate whether these sites have chromatin higher-order structure. For example, comparison between averaged active or inactive promoter regions suggests that both regions can acquire stabilized higher-order structure with hidden H1 epitopes. However, the H1 xChIP-seq comparison cannot define their differences. Application of the xxChIP-seq versus H1 xChIP-seq method is particularly relevant to chromatin-associated proteins, such as linker histones, that play dynamic roles in establishing chromatin higher-order structure.

Structure encoding in DNA

It is proposed that transposons and related long non-coding RNA define the fine structure of body parts. Although morphogens have long been known to direct the formation of many gross structures in early embryonic development, they do not have the necessary precision to define a structure down to the individual cellular level. Using the distinction between procedural and declarative knowledge in information processing as an analogy, it is hypothesized that DNA encodes fine structure in a manner that is different from the genetic code for proteins. The hypothesis states that repeated or near-repeated sequences that are in transposons and non-coding RNA define body part structures. As the cells in a body part go through the epigenetic process of differentiation, the action of methylation serves to inactivate all but the relevant structure definitions and some associated cell type genes. The transposons left active will then physically modify the DNA sequence in the heterochromatin to establish the local context in the three-dimensional body part structure. This brings the encoded definition of the cell type to the histone. The histone code for that cell type starts the regulatory cascade that turns on the genes associated with that particular type of cell, transforming it from a multipotent cell to a fully differentiated cell. This mechanism creates structures in the musculoskeletal system, the organs of the body, the major parts of the brain, and other systems.

GsmPlot: a web server to visualize epigenome data in NCBI

Background

Epigenetic regulation is essential in regulating gene expression across a variety of biological processes. Many high-throughput sequencing technologies have been widely used to generate epigenetic data, such as histone modification, transcription factor binding sites, DNA modifications, chromatin accessibility, and etc. A large scale of epigenetic data is stored in NCBI Gene Expression Omnibus (GEO). However, it is a great challenge to reanalyze these large scale and complex data, especially for researchers who do not specialize in bioinformatics skills or do not have access to expensive computational infrastructure.

Results

GsmPlot can simply accept GSM IDs to automatically download NCBI data or can accept user’s private bigwig files as input to plot the concerned data on promoters, exons or any other user-defined genome locations and generate UCSC visualization tracks. By linking public data repository and private data, GsmPlot can spark data-driven ideas and hence promote the epigenetic research.

Conclusions

GsmPlot web server allows convenient visualization and efficient exploration of any NCBI epigenetic data in any genomic region without need of any bioinformatics skills or special computing resources. GsmPlot is freely available at https://gsmplot.deqiangsun.org/.

CTCF-dependent chromatin boundaries formed by asymmetric nucleosome arrays with decreased linker length

The CCCTC-binding factor (CTCF) organises the genome in 3D through DNA loops and in 1D by setting boundaries isolating different chromatin states, but these processes are not well understood. Here we investigate chromatin boundaries in mouse embryonic stem cells, defined by the regions with decreased Nucleosome Repeat Length (NRL) for ∼20 nucleosomes near CTCF sites, affecting up to 10% of the genome. We found that the nucleosome-depleted region (NDR) near CTCF is asymmetrically located >40 nucleotides 5'-upstream from the centre of CTCF motif. The strength of CTCF binding to DNA and the presence of cohesin is correlated with the decrease of NRL near CTCF, and anti-correlated with the level of asymmetry of the nucleosome array. Individual chromatin remodellers have different contributions, with Snf2h having the strongest effect on the NRL decrease near CTCF and Chd4 playing a major role in the symmetry breaking. Upon differentiation, a subset of preserved, common CTCF sites maintains asymmetric nucleosome pattern and small NRL. The sites which lost CTCF upon differentiation are characterized by nucleosome rearrangement 3'-downstream, with unchanged NDR 5'-upstream of CTCF motifs. Boundaries of topologically associated chromatin domains frequently contain several inward-oriented CTCF motifs whose effects, described above, add up synergistically.

CTCF confers local nucleosome resiliency after DNA replication and during mitosis

The access of Transcription Factors (TFs) to their cognate DNA binding motifs requires a precise control over nucleosome positioning. This is especially important following DNA replication and during mitosis, both resulting in profound changes in nucleosome organization over TF binding regions. Using mouse Embryonic Stem (ES) cells, we show that the TF CTCF displaces nucleosomes from its binding site and locally organizes large and phased nucleosomal arrays, not only in interphase steady-state but also immediately after replication and during mitosis. Correlative analyses suggest this is associated with fast gene reactivation following replication and mitosis. While regions bound by other TFs (Oct4/Sox2), display major rearrangement, the post-replication and mitotic nucleosome positioning activity of CTCF is not unique: Esrrb binding regions are also characterized by persistent nucleosome positioning. Therefore, selected TFs such as CTCF and Esrrb act as resilient TFs governing the inheritance of nucleosome positioning at regulatory regions throughout the cell-cycle.

TET-Catalyzed 5-Carboxylcytosine Promotes CTCF Binding to Suboptimal Sequences Genome-wide

The mechanisms supporting dynamic regulation of CTCF-binding sites remain poorly understood. Here we describe the TET-catalyzed 5-methylcytosine derivative, 5-carboxylcytosine (5caC), as a factor driving new CTCF binding within genomic DNA. Through a combination of in vivo and in vitro approaches, we reveal that 5caC generally strengthens CTCF association with DNA and facilitates binding to suboptimal sequences. Dramatically, profiling of CTCF binding in a cellular model that accumulates genomic 5caC identified ~13,000 new CTCF sites. The new sites were enriched for overlapping 5caC and were marked by an overall reduction in CTCF motif strength. As CTCF has multiple roles in gene expression, these findings have wide-reaching implications and point to induced 5caC as a potential mechanism to achieve differential CTCF binding in cells.

Quantitative MNase-seq accurately maps nucleosome occupancy levels

Micrococcal nuclease (MNase) is widely used to map nucleosomes. However, its aggressive endo-/exo-nuclease activities make MNase-seq unreliable for determining nucleosome occupancies, because cleavages within linker regions produce oligo- and mono-nucleosomes, whereas cleavages within nucleosomes destroy them. Here, we introduce a theoretical framework for predicting nucleosome occupancies and an experimental protocol with appropriate spike-in normalization that confirms our theory and provides accurate occupancy levels over an MNase digestion time course. As with human cells, we observe no overall differences in nucleosome occupancies between Drosophila euchromatin and heterochromatin, which implies that heterochromatic compaction does not reduce MNase accessibility of linker DNA.

TET1 regulates fibroblast growth factor 8 transcription in gonadotropin releasing hormone neurons

Fibroblast growth factor 8 (FGF8) is a potent morphogen that regulates the ontogenesis of gonadotropin-releasing hormone (GnRH) neurons, which control the hypothalamus-pituitary-gonadal (HPG) axis, and therefore reproductive success. Indeed, FGF8 and FGFR1 deficiency severely compromises vertebrate reproduction in mice and humans and is associated with Kallmann Syndrome (KS), a congenital disease characterized by hypogonadotropic hypogonadism associated with anosmia. Our laboratory demonstrated that FGF8 signaling through FGFR1, both of which are KS-related genes, is necessary for proper GnRH neuron development in mice and humans. Here, we investigated the possibility that non-genetic factors, such as the epigenome, may contribute to KS onset. For this purpose, we developed an embryonic explant model, utilizing the mouse olfactory placode (OP), the birthplace of GnRH neurons. We show that TET1, which converts 5-methylcytosine residues (5mC) to 5-hydroxymethylated cytosines (5hmC), controls transcription of Fgf8 during GnRH neuron ontogenesis. Through MeDIP and ChIP RT-qPCR we found that TET1 bound to specific CpG islands on the Fgf8 promoter. We found that the temporal expression of Fgf8 correlates with not only TET1 binding, but also with 5hmC enrichment. siRNA knockdown of Tet1 reduced Fgf8 and Fgfr1 mRNA expression. During this time period, Fgf8 also switched histone status, most likely via recruitment of EZH2, a major component of the polycomb repressor complex-2 (PRC2) at E13.5. Together, these studies underscore the significance of epigenetics and chromatin modifications to temporally regulated genes involved in KS.

Acute depletion of CTCF directly affects MYC regulation through loss of enhancer-promoter looping

Numerous pieces of evidence support the complex, 3D spatial organization of the genome dictates gene expression. CTCF is essential to define topologically associated domain boundaries and to facilitate the formation of insulated chromatin loop structures. To understand CTCF's direct role in global transcriptional regulation, we integrated the miniAID-mClover3 cassette to the endogenous CTCF locus in a human pediatric B-ALL cell line, SEM, and an immortal erythroid precursor cell line, HUDEP-2, to allow for acute depletion of CTCF protein by the auxin-inducible degron system. In SEM cells, CTCF loss notably disrupted intra-TAD loops and TAD integrity in concurrence with a reduction in CTCF-binding affinity, while showing no perturbation to nuclear compartment integrity. Strikingly, the overall effect of CTCF's loss on transcription was minimal. Whole transcriptome analysis showed hundreds of genes differentially expressed in CTCF-depleted cells, among which MYC and a number of MYC target genes were specifically downregulated. Mechanically, acute depletion of CTCF disrupted the direct interaction between the MYC promoter and its distal enhancer cluster residing ∼1.8 Mb downstream. Notably, MYC expression was not profoundly affected upon CTCF loss in HUDEP-2 cells suggesting that CTCF could play a B-ALL cell line specific role in maintaining MYC expression.

Chromatin accessibility and the regulatory epigenome

Physical access to DNA is a highly dynamic property of chromatin that plays an essential role in establishing and maintaining cellular identity. The organization of accessible chromatin across the genome reflects a network of permissible physical interactions through which enhancers, promoters, insulators and chromatin-binding factors cooperatively regulate gene expression. This landscape of accessibility changes dynamically in response to both external stimuli and developmental cues, and emerging evidence suggests that homeostatic maintenance of accessibility is itself dynamically regulated through a competitive interplay between chromatin-binding factors and nucleosomes. In this Review, we examine how the accessible genome is measured and explore the role of transcription factors in initiating accessibility remodelling; our goal is to illustrate how chromatin accessibility defines regulatory elements within the genome and how these epigenetic features are dynamically established to control gene expression.

BORIS Expression in Ovarian Cancer Precursor Cells Alters the CTCF Cistrome and Enhances Invasiveness through GALNT14

High-grade serous carcinoma (HGSC) is the most aggressive and predominant form of epithelial ovarian cancer and the leading cause of gynecologic cancer-related death. We have previously shown that CTCFL (also known as BORIS, Brother of the Regulator of Imprinted Sites) is expressed in most ovarian cancers, and is associated with global and promoter-specific DNA hypomethylation, advanced tumor stage, and poor prognosis. To explore its role in HGSC, we expressed BORIS in human fallopian tube secretory epithelial cells (FTSEC), the presumptive cells of origin for HGSC. BORIS-expressing cells exhibited increased motility and invasion, and BORIS expression was associated with alterations in several cancer-associated gene expression networks, including fatty acid metabolism, TNF signaling, cell migration, and ECM-receptor interactions. Importantly, GALNT14, a glycosyltransferase gene implicated in cancer cell migration and invasion, was highly induced by BORIS, and GALNT14 knockdown significantly abrogated BORIS-induced cell motility and invasion. In addition, in silico analyses provided evidence for BORIS and GALNT14 coexpression in several cancers. Finally, ChIP-seq demonstrated that expression of BORIS was associated with de novo and enhanced binding of CTCF at hundreds of loci, many of which correlated with activation of transcription at target genes, including GALNT14. Taken together, our data indicate that BORIS may promote cell motility and invasion in HGSC via upregulation of GALNT14, and suggests BORIS as a potential therapeutic target in this malignancy. IMPLICATIONS: These studies provide evidence that aberrant expression of BORIS may play a role in the progression to HGSC by enhancing the migratory and invasive properties of FTSEC.

Linking aberrant chromatin features in chronic lymphocytic leukemia to transcription factor networks

In chronic lymphocytic leukemia (CLL), a diverse set of genetic mutations is embedded in a deregulated epigenetic landscape that drives cancerogenesis. To elucidate the role of aberrant chromatin features, we mapped DNA methylation, seven histone modifications, nucleosome positions, chromatin accessibility, binding of EBF1 and CTCF, as well as the transcriptome of B cells from CLL patients and healthy donors. A globally increased histone deacetylase activity was detected and half of the genome comprised transcriptionally downregulated partially DNA methylated domains demarcated by CTCF CLL samples displayed a H3K4me3 redistribution and nucleosome gain at promoters as well as changes of enhancer activity and enhancer linkage to target genes. A DNA binding motif analysis identified transcription factors that gained or lost binding in CLL at sites with aberrant chromatin features. These findings were integrated into a gene regulatory enhancer containing network enriched for B-cell receptor signaling pathway components. Our study predicts novel molecular links to targets of CLL therapies and provides a valuable resource for further studies on the epigenetic contribution to the disease.

DNA (de)methylation in embryonic stem cells controls CTCF-dependent chromatin boundaries

Coordinated changes of DNA (de)methylation, nucleosome positioning, and chromatin binding of the architectural protein CTCF play an important role for establishing cell-type-specific chromatin states during differentiation. To elucidate molecular mechanisms that link these processes, we studied the perturbed DNA modification landscape in mouse embryonic stem cells (ESCs) carrying a double knockout (DKO) of the Tet1 and Tet2 dioxygenases. These enzymes are responsible for the conversion of 5-methylcytosine (5mC) into its hydroxymethylated (5hmC), formylated (5fC), or carboxylated (5caC) forms. We determined changes in nucleosome positioning, CTCF binding, DNA methylation, and gene expression in DKO ESCs and developed biophysical models to predict differential CTCF binding. Methylation-sensitive nucleosome repositioning accounted for a significant portion of CTCF binding loss in DKO ESCs, whereas unmethylated and nucleosome-depleted CpG islands were enriched for CTCF sites that remained occupied. A number of CTCF sites also displayed direct correlations with the CpG modification state: CTCF was preferentially lost from sites that were marked with 5hmC in wild-type (WT) cells but not from 5fC-enriched sites. In addition, we found that some CTCF sites can act as bifurcation points defining the differential methylation landscape. CTCF loss from such sites, for example, at promoters, boundaries of chromatin loops, and topologically associated domains (TADs), was correlated with DNA methylation/demethylation spreading and can be linked to down-regulation of neighboring genes. Our results reveal a hierarchical interplay between cytosine modifications, nucleosome positions, and DNA sequence that determines differential CTCF binding and regulates gene expression.

X-ray irradiation induces subtle changes in the genome-wide distribution of DNA hydroxymethylation with opposing trends in genic and intergenic regions

DNA hydroxymethylation has gained attention as an intermediate in the process of DNA demethylation. More recently, 5-hydroxymethylcytosine has been recognized as an independent epigenetic mark that can persist over time and that exerts influence on gene regulation and other biological processes. Deregulation of this DNA modification has been linked to tumorigenesis and a variety of other diseases. The impact of irradiation on DNA hydroxymethylation is poorly understood. In this study we exposed lung fibroblasts (IMR90) to 0.5 Gy and 2 Gy of X-rays, respectively. We characterized radiation induced changes of DNA hydroxymethylation 1 h, 6 h, 24 h and 120 h after exposure employing immunoprecipitation and subsequent deep sequencing of the genomic fraction enriched for hydroxymethylated DNA. Transcriptomic response to irradiation was analyzed for time points 6 h and 24 h post exposure by means of RNA sequencing. Irradiated and sham-irradiated samples shared the same overall distribution of 5-hydroxymethylcytosines with respect to genomic features such as promoters and exons. The frequency of 5-hydroxymethylcytosine peaks differentially detected in irradiated samples increased in genic regions over time, while the opposing trend was observed for intergenic regions. Onset and extent of this effect was dose dependent. Moreover, we demonstrated a biased distribution of 5-hmC alterations at CpG islands and sites occupied by the DNA binding protein CTCF. In summary, our study provides new insights into the epigenetic response to irradiation. Our data highlight genomic features more prone to irradiation induced changes of DNA hydroxymethylation, which might impact early and late onset effects of irradiation.

Transcription factor activity and nucleosome organization in mitosis

Mitotic bookmarking transcription factors (BFs) maintain the capacity to bind to their targets during mitosis, despite major rearrangements of the chromatin. While they were thought to propagate gene regulatory information through mitosis by statically occupying their DNA targets, it has recently become clear that BFs are highly dynamic in mitotic cells. This represents both a technical and a conceptual challenge to study and understand the function of BFs: First, formaldehyde has been suggested to be unable to efficiently capture these transient interactions, leading to profound contradictions in the literature; and second, if BFs are not permanently bound to their targets during mitosis, it becomes unclear how they convey regulatory information to daughter cells. Here, comparing formaldehyde to alternative fixatives we clarify the nature of the chromosomal association of previously proposed BFs in embryonic stem cells: While ESRRB can be considered as a canonical BF that binds at selected regulatory regions in mitosis, SOX2 and POU5F1 (also known as OCT4) establish DNA sequence-independent interactions with the mitotic chromosomes, either throughout the chromosomal arms (SOX2) or at pericentromeric regions (POU5F1). Moreover, we show that ordered nucleosomal arrays are retained during mitosis at ESRRB bookmarked sites, whereas regions losing transcription factor binding display a profound loss of order. By maintaining nucleosome positioning during mitosis, ESRRB might ensure the rapid post-mitotic re-establishment of functional regulatory complexes at selected enhancers and promoters. Our results provide a mechanistic framework that reconciles dynamic mitotic binding with the transmission of gene regulatory information across cell division.

5-Formylcytosine organizes nucleosomes and forms Schiff base interactions with histones in mouse embryonic stem cells

Nucleosomes are the basic unit of chromatin that help the packaging of genetic material while controlling access to the genetic information. The underlying DNA sequence, together with transcription-associated proteins and chromatin remodelling complexes, are important factors that influence the organization of nucleosomes. Here, we show that the naturally occurring DNA modification, 5-formylcytosine (5fC) is linked to tissue-specific nucleosome organization. Our study reveals that 5fC is associated with increased nucleosome occupancy in vitro and in vivo. We demonstrate that 5fC-associated nucleosomes at enhancers in the mammalian hindbrain and heart are linked to elevated gene expression. Our study also reveals the formation of a reversible-covalent Schiff base linkage between lysines of histone proteins and 5fC within nucleosomes in a cellular environment. We define their specific genomic loci in mouse embryonic stem cells and look into the biological consequences of these DNA-histone Schiff base sites. Collectively, our findings show that 5fC is a determinant of nucleosome organization and plays a role in establishing distinct regulatory regions that control transcription.