Nucleosome Organization in Human Embryonic Stem Cells

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.

Epigenetic regulation of early human embryo development

Studies of mammalian development have advanced our understanding of the genetic, epigenetic, and cellular processes that orchestrate embryogenesis and have uncovered new insights into the unique aspects of human embryogenesis. Recent studies have now produced the first epigenetic maps of early human embryogenesis, stimulating new ideas about epigenetic reprogramming, cell fate control, and the potential mechanisms underpinning developmental plasticity in human embryos. In this review, we discuss these new insights into the epigenetic regulation of early human development and the importance of these processes for safeguarding development. We also highlight unanswered questions and key challenges that remain to be addressed.

DNA methylation: Precise modulation of chromatin structure and dynamics

DNA methylation plays a vital role in epigenetic regulation in both plants and animals, and typically occurs at the 5-carbon position of the cytosine pyrimidine ring within the CpG dinucleotide steps. Cytosine methylation can alter DNA's geometry, mechanical and physico-chemical properties - thus influencing the molecular signaling events vital for transcription, replication and chromatin remodeling. Despite the profound effect cytosine methylation can have on DNA, the underlying atomistic mechanisms remain enigmatic. Many studies so far have produced controversial findings on how cytosine methylation dictates DNA flexibility and accessibility, nucleosome stability and dynamics. Here, we review the most recent experimental and computational studies that provide precise characterization of structure and function of cytosine methylation and its versatile roles in modulating DNA mechanics, nucleosome and chromatin structure, stability and dynamics. Moreover, the review briefly discusses the relationship between DNA methylation and nucleosome positioning, and the crosstalk between DNA methylation and histone tail modifications.

iNucs: Inter-Nucleosome Interactions

MOTIVATION

Deciphering nucleosome-nucleosome interactions is an important step towards mesoscale description of chromatin organization but computational tools to perform such analyses are not publicly available.

RESULTS

We developed iNucs, a user-friendly and efficient Python-based bioinformatics tool to compute and visualize nucleosome resolved interactions using standard pairs format input generated from pairtools.

AVAILABILITY

https://github.com/Karimi-Lab/inucs/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

NucHMM: a method for quantitative modeling of nucleosome organization identifying functional nucleosome states distinctly associated with splicing potentiality

We develop a novel computational method, NucHMM, to identify functional nucleosome states associated with cell type-specific combinatorial histone marks and nucleosome organization features such as phasing, spacing and positioning. We test it on publicly available MNase-seq and ChIP-seq data in MCF7, H1, and IMR90 cells and identify 11 distinct functional nucleosome states. We demonstrate these nucleosome states are distinctly associated with the splicing potentiality of skipping exons. This advances our understanding of the chromatin function at the nucleosome level and offers insights into the interplay between nucleosome organization and splicing processes.

Chromatin Profiling Techniques: Exploring the Chromatin Environment and Its Contributions to Complex Traits

The genetic architecture of complex traits is multifactorial. Genome-wide association studies (GWASs) have identified risk loci for complex traits and diseases that are disproportionately located at the non-coding regions of the genome. On the other hand, we have just begun to understand the regulatory roles of the non-coding genome, making it challenging to precisely interpret the functions of non-coding variants associated with complex diseases. Additionally, the epigenome plays an active role in mediating cellular responses to fluctuations of sensory or environmental stimuli. However, it remains unclear how exactly non-coding elements associate with epigenetic modifications to regulate gene expression changes and mediate phenotypic outcomes. Therefore, finer interrogations of the human epigenomic landscape in associating with non-coding variants are warranted. Recently, chromatin-profiling techniques have vastly improved our understanding of the numerous functions mediated by the epigenome and DNA structure. Here, we review various chromatin-profiling techniques, such as assays of chromatin accessibility, nucleosome distribution, histone modifications, and chromatin topology, and discuss their applications in unraveling the brain epigenome and etiology of complex traits at tissue homogenate and single-cell resolution. These techniques have elucidated compositional and structural organizing principles of the chromatin environment. Taken together, we believe that high-resolution epigenomic and DNA structure profiling will be one of the best ways to elucidate how non-coding genetic variations impact complex diseases, ultimately allowing us to pinpoint cell-type targets with therapeutic potential.

Enhancers predominantly regulate gene expression during differentiation via transcription initiation

Gene transcription occurs via a cycle of linked events, including initiation, promoter-proximal pausing, and elongation of RNA polymerase II (Pol II). A key question is how transcriptional enhancers influence these events to control gene expression. Here, we present an approach that evaluates the level and change in promoter-proximal transcription (initiation and pausing) in the context of differential gene expression, genome-wide. This combinatorial approach shows that in primary cells, control of gene expression during differentiation is achieved predominantly via changes in transcription initiation rather than via release of Pol II pausing. Using genetically engineered mouse models, deleted for functionally validated enhancers of the α- and β-globin loci, we confirm that these elements regulate Pol II recruitment and/or initiation to modulate gene expression. Together, our data show that gene expression during differentiation is regulated predominantly at the level of initiation and that enhancers are key effectors of this process.

A High-Resolution Map of Human Enhancer RNA Loci Characterizes Super-enhancer Activities in Cancer

Although enhancers play critical roles in cancer, quantifying enhancer activities in clinical samples remains challenging, especially for super-enhancers. Enhancer activities can be inferred from enhancer RNA (eRNA) signals, which requires enhancer transcription loci definition. Only a small proportion of human eRNA loci has been precisely identified, limiting investigations of enhancer-mediated oncogenic mechanisms. Here, we characterize super-enhancer regions using aggregated RNA sequencing (RNA-seq) data from large cohorts. Super-enhancers usually contain discrete loci featuring sharp eRNA expression peaks. We identify >300,000 eRNA loci in ∼377 Mb super-enhancer regions that are regulated by evolutionarily conserved, well-positioned nucleosomes and are frequently dysregulated in cancer. The eRNAs provide explanatory power for cancer phenotypes beyond that provided by mRNA expression through resolving intratumoral heterogeneity with enhancer cell-type specificity. Our study provides a high-resolution map of eRNA loci through which super-enhancer activities can be quantified by RNA-seq and a user-friendly data portal, enabling a broad range of biomedical investigations.

Nonreciprocal and Conditional Cooperativity Directs the Pioneer Activity of Pluripotency Transcription Factors

Cooperative binding of transcription factors (TFs) to chromatin orchestrates gene expression programming and cell fate specification. However, the biophysical principles of TF cooperativity remain incompletely understood. Here we use single-molecule fluorescence microscopy to study the partnership between Sox2 and Oct4, two core members of the pluripotency gene regulatory network. We find that the ability of Sox2 to target DNA inside nucleosomes is strongly affected by the translational and rotational positioning of its binding motif. In contrast, Oct4 can access nucleosomal sites with equal capacities. Furthermore, the Sox2-Oct4 pair displays nonreciprocal cooperativity, with Oct4 modulating interaction of Sox2 with the nucleosome but not vice versa. Such cooperativity is conditional upon the composite motif's residing at specific nucleosomal locations. These results reveal that pioneer factors possess distinct chromatin-binding properties and suggest that the same set of TFs can differentially regulate gene activities on the basis of their motif positions in the nucleosomal context.

DNMT3B shapes the mCA landscape and regulates mCG for promoter bivalency in human embryonic stem cells

DNMT3B is known as a de novo DNA methyltransferase. However, its preferential target sites for DNA methylation are largely unknown. Our analysis on ChIP-seq experiment in human embryonic stem cells (hESC) revealed that DNMT3B, mCA and H3K36me3 share the same genomic distribution profile. Deletion of DNMT3B or its histone-interacting domain (PWWP) demolished mCA in hESCs, suggesting that PWWP domain of DNMT3B directs the formation of mCA landscape. In contrast to the common presumption that PWWP guides DNMT3B-mediated mCG deposition, we found that deleting PWWP does not affect the mCG landscape. Nonetheless, DNMT3B knockout led to the formation of 2985 de novo hypomethylated regions at annotated promoter sites. Upon knockout, most of these promoters gain the bivalent marks, H3K4me3 and H3K27me3. We call them spurious bivalent promoters. Gene ontology analysis associated spurious bivalent promoters with development and cell differentiation. Overall, we found the importance of DNMT3B for shaping the mCA landscape and for maintaining the fidelity of the bivalent promoters in hESCs.

Human SETMAR is a DNA sequence-specific histone-methylase with a broad effect on the transcriptome

Transposons impart dynamism to the genomes they inhabit and their movements frequently rewire the control of nearby genes. Occasionally, their proteins are domesticated when they evolve a new function. SETMAR is a protein methylase with a sequence-specific DNA binding domain. It began to evolve about 50 million years ago when an Hsmar1 transposon integrated downstream of a SET-domain methylase gene. Here we show that the DNA-binding domain of the transposase targets the enzyme to transposon-end remnants and that this is capable of regulating gene expression, dependent on the methylase activity. When SETMAR was modestly overexpressed in human cells, almost 1500 genes changed expression by more than 2-fold (65% up- and 35% down-regulated). These genes were enriched for the KEGG Pathways in Cancer and include several transcription factors important for development and differentiation. Expression of a similar level of a methylase-deficient SETMAR changed the expression of many fewer genes, 77% of which were down-regulated with no significant enrichment of KEGG Pathways. Our data is consistent with a model in which SETMAR is part of an anthropoid primate-specific regulatory network centered on the subset of genes containing a transposon end.

Mitochondrial Akt Signaling Modulated Reprogramming of Somatic Cells

The signaling mechanisms controlling somatic cell reprogramming are not fully understood. In this study, we report a novel role for mitochondrial Akt1 signaling that enhanced somatic cell reprogramming efficiency. The role of mitochondrial Akt1 in somatic cell reprogramming was investigated by transducing fibroblasts with the four reprogramming factors (Oct4, Sox2, Klf4, c-Myc) in conjunction with Mito-Akt1, Mito-dnAkt1, or control virus. Mito-Akt1 enhanced reprogramming efficiency whereas Mito-dnAkt1 inhibited reprogramming. The resulting iPSCs formed embryoid bodies in vitro and teratomas in vivo. Moreover, Oct4 and Nanog promoter methylation was reduced in the iPSCs generated in the presence of Mito-Akt1. Akt1 was activated and translocated into mitochondria after growth factor stimulation in embryonic stem cells (ESCs). To study the effect of mitochondrial Akt in ESCs, a mitochondria-targeting constitutively active Akt1 (Mito-Akt1) was expressed in ESCs. Gene expression profiling showed upregulation of genes that promote stem cell proliferation and survival and down-regulation of genes that promote differentiation. Analysis of cellular respiration indicated similar metabolic profile in the resulting iPSCs and ESCs, suggesting comparable bioenergetics. These findings showed that activation of mitochondrial Akt1 signaling was required during somatic cell reprogramming.

A Scaled Framework for CRISPR Editing of Human Pluripotent Stem Cells to Study Psychiatric Disease

Scaling of CRISPR-Cas9 technology in human pluripotent stem cells (hPSCs) represents an important step for modeling complex disease and developing drug screens in human cells. However, variables affecting the scaling efficiency of gene editing in hPSCs remain poorly understood. Here, we report a standardized CRISPR-Cas9 approach, with robust benchmarking at each step, to successfully target and genotype a set of psychiatric disease-implicated genes in hPSCs and provide a resource of edited hPSC lines for six of these genes. We found that transcriptional state and nucleosome positioning around targeted loci was not correlated with editing efficiency. However, editing frequencies varied between different hPSC lines and correlated with genomic stability, underscoring the need for careful cell line selection and unbiased assessments of genomic integrity. Together, our step-by-step quantification and in-depth analyses provide an experimental roadmap for scaling Cas9-mediated editing in hPSCs to study psychiatric disease, with broader applicability for other polygenic diseases.

Exploring the Link between Nucleosome Occupancy and DNA Methylation

Near promoters, both nucleosomes and CpG sites form characteristic spatial patterns. Previously, nucleosome depleted regions were observed upstream of transcription start sites and nucleosome occupancy was reported to correlate both with CpG density and the level of CpG methylation. Several studies imply a causal link where CpG methylation might induce nucleosome formation, whereas others argue the opposite, i.e., that nucleosome occupancy might influence CpG methylation. Correlations are indeed evident between nucleosomes, CpG density and CpG methylation—at least near promoter sites. It is however less established whether there is an immediate causal relation between nucleosome occupancy and the presence of CpG sites—or if nucleosome occupancy could be influenced by other factors. In this work, we test for such causality in human genomes by analyzing the three quantities both near and away from promoter sites. For data from the human genome we compare promoter regions with given CpG densities with genomic regions without promoters but of similar CpG densities. We find the observed correlation between nucleosome occupancy and CpG density, respectively CpG methylation, to be specific to promoter regions. In other regions along the genome nucleosome occupancy is statistically independent of the positioning of CpGs or their methylation levels. Anti-correlation between CpG density and methylation level is however similarly strong in both regions. On promoters, nucleosome occupancy is more strongly affected by the level of gene expression than CpG density or CpG methylation—calling into question any direct causal relation between nucleosome occupancy and CpG organization. Rather, our results suggest that for organisms with cytosine methylation nucleosome occupancy might be primarily linked to gene expression, with no strong impact on methylation.

Biophysical characterization of the association of histones with single-stranded DNA

BACKGROUND

Despite the profound current knowledge of the architecture and dynamics of nucleosomes, little is known about the structures generated by the interaction of histones with single-stranded DNA (ssDNA), which is widely present during replication and transcription.

METHODS

Non-denaturing gel electrophoresis, transmission electron microscopy, atomic force microscopy, magnetic tweezers.

RESULTS

Histones have a high affinity for ssDNA in 0.15M NaCl ionic strength, with an apparent binding constant similar to that calculated for their association with double-stranded DNA (dsDNA). The length of DNA (number of nucleotides in ssDNA or base pairs in dsDNA) associated with a fixed core histone mass is the same for both ssDNA and dsDNA. Although histone-ssDNA complexes show a high tendency to aggregate, nucleosome-like structures are formed at physiological salt concentrations. Core histones are able to protect ssDNA from digestion by micrococcal nuclease, and a shortening of ssDNA occurs upon its interaction with histones. The purified (+) strand of a cloned DNA fragment of nucleosomal origin has a higher affinity for histones than the purified complementary (-) strand.

CONCLUSIONS

At physiological ionic strength histones have high affinity for ssDNA, possibly associating with it into nucleosome-like structures.

GENERAL SIGNIFICANCE

In the cell nucleus histones may spontaneously interact with ssDNA to facilitate their participation in the replication and transcription of chromatin.

A comparison of nucleosome organization in Drosophila cell lines

Changes in the distribution of nucleosomes along the genome influence chromatin structure and impact gene expression by modulating the accessibility of DNA to transcriptional machinery. However, the role of genome-wide nucleosome positioning in gene expression and in maintaining differentiated cell states remains poorly understood. Drosophila melanogaster cell lines represent distinct tissue types and exhibit cell-type specific gene expression profiles. They thus could provide a useful tool for investigating cell-type specific nucleosome organization of an organism's genome. To evaluate this possibility, we compared genome-wide nucleosome positioning and occupancy in five different Drosophila tissue-specific cell lines, and in reconstituted chromatin, and then tested for correlations between nucleosome positioning, transcription factor binding motifs, and gene expression. Nucleosomes in all cell lines were positioned in accordance with previously known DNA-nucleosome interactions, with helically repeating A/T di-nucleotide pairs arranged within nucleosomal DNAs and AT-rich pentamers generally excluded from nucleosomal DNA. Nucleosome organization in all cell lines differed markedly from in vitro reconstituted chromatin, with highly expressed genes showing strong nucleosome organization around transcriptional start sites. Importantly, comparative analysis identified genomic regions that exhibited cell line-specific nucleosome enrichment or depletion. Further analysis of these regions identified 91 out of 16,384 possible heptamer sequences that showed differential nucleosomal occupation between cell lines, and 49 of the heptamers matched one or more known transcription factor binding sites. These results demonstrate that there is differential nucleosome positioning between these Drosophila cell lines and therefore identify a system that could be used to investigate the functional significance of differential nucleosomal positioning in cell type specification.

Structural basis for recognition of histone H3K36me3 nucleosome by human de novo DNA methyltransferases 3A and 3B

DNA methylation is an important epigenetic modification involved in chromatin organization and gene expression. The function of DNA methylation depends on cell context and is correlated with histone modification patterns. In particular, trimethylation of Lys36 on histone H3 tail (H3K36me3) is associated with DNA methylation and elongation phase of transcription. PWWP domains of the de novo DNA methyltransferases DNMT3A and DNMT3B read this epigenetic mark to guide DNA methylation. Here we report the first crystal structure of the DNMT3B PWWP domain-H3K36me3 complex. Based on this structure, we propose a model of the DNMT3A PWWP domain-H3K36me3 complex and build a model of DNMT3A (PWWP-ADD-CD) in a nucleosomal context. The trimethylated side chain of Lys36 (H3K36me3) is inserted into an aromatic cage similar to the "Royal" superfamily domains known to bind methylated histones. A key interaction between trimethylated Lys36 and a conserved water molecule stabilized by Ser270 explains the lack of affinity of mutated DNMT3B (S270P) for the H3K36me3 epigenetic mark in the ICF (Immunodeficiency, Centromeric instability and Facial abnormalities) syndrome. The model of the DNMT3A-DNMT3L heterotetramer in complex with a dinucleosome highlights the mechanism for recognition of nucleosome by DNMT3s and explains the periodicity of de novo DNA methylation.

Non-random fragmentation patterns in circulating cell-free DNA reflect epigenetic regulation

BACKGROUND

The assessment of cell-free circulating DNA fragments, also known as a "liquid biopsy" of the patient's plasma, is an important source for the discovery and subsequent non-invasive monitoring of cancer and other pathological conditions. Although the nucleosome-guided fragmentation patterns of cell-free DNA (cfDNA) have not yet been studied in detail, non-random representation of cfDNA sequencies may reflect chromatin features in the tissue of origin at gene-regulation level.

RESULTS

In this study, we investigated the association between epigenetic landscapes of human tissues evident in the patterns of cfDNA in plasma by deep sequencing of human cfDNA samples. We have demonstrated that baseline characteristics of cfDNA fragmentation pattern are in concordance with the ones corresponding to cell lines-derived. To identify the loci differentially represented in cfDNA fragment, we mapped the transcription start sites within the sequenced cfDNA fragments and tested for association of these genomic coordinates with the relative strength and the patterns of gene expressions. Preselected sets of house-keeping and tissue specific genes were used as models for actively expressed and silenced genes. Developed measure of gene regulation was able to differentiate these two sets based on sequencing coverage near gene transcription start site.

CONCLUSION

Experimental outcomes suggest that cfDNA retains characteristics previously noted in genome-wide analysis of chromatin structure, in particular, in MNase-seq assays. Thus far the analysis of the DNA fragmentation pattern may aid further developing of cfDNA based biomarkers for a variety of human conditions.

Associations between nucleosome phasing, sequence asymmetry, and tissue-specific expression in a set of inbred Medaka species

BACKGROUND

Transcription start sites (TSSs) with pronounced and phased nucleosome arrays downstream and nucleosome-depleted regions upstream of TSSs are observed in various species.

RESULTS

We have characterized sequence variation and expression properties of this set of TSSs (which we call "Nucleocyclic TSSs") using germline and somatic cells of three medaka (Oryzias latipes) inbred isolates from different locations. We found nucleocyclic TSSs in medaka to be associated with higher gene expression and characterized by a clear boundary in sequence composition with potentially-nucleosome-destabilizing A/T-enrichment upstream (p < 10(-60)) and nucleosome- accommodating C/G-enrichment downstream (p < 10(-40)) that was highly conserved from an ancestor. A substantial genetic distance between the strains facilitated the in-depth analysis of patterns of fixed mutations, revealing a localization-specific equilibrium between the rates of distinct mutation categories that would serve to maintain the conserved sequence anisotropy around TSSs. Downstream of nucleocyclic TSSs, C to T, T to C, and other mutation rates on the sense strand increased around first nucleosome dyads and decreased around first linkers, which contrasted with genomewide mutational patterns around nucleosomes (p < 5 %). C to T rates are higher than G to A rates around nucleosome associated with germline nucleocyclic TSS sites (p < 5 %), potentially due to the asymmetric effect of transcription-coupled repair.

CONCLUSIONS

Our results demonstrate an atypical evolutionary process surrounding nucleocyclic TSSs.