Estrogen-induced chromatin decondensation and nuclear re-organization linked to regional epigenetic regulation in breast cancer

BackgroundEpigenetic changes are being increasingly recognized as a prominent feature of cancer. This occurs not only at individual genes, but also over larger chromosomal domains. To investigate this, we set out to identify large chromosomal domains of epigenetic dysregulation in breast cancers.ResultsWe identify large regions of coordinate down-regulation of gene expression, and other regions of coordinate activation, in breast cancers and show that these regions are linked to tumor subtype. In particular we show that a group of coordinately regulated regions are expressed in luminal, estrogen-receptor positive breast tumors and cell lines. For one of these regions of coordinate gene activation, we show that regional epigenetic regulation is accompanied by visible unfolding of large-scale chromatin structure and a repositioning of the region within the nucleus. In MCF7 cells, we show that this depends on the presence of estrogen.ConclusionsOur data suggest that the liganded estrogen receptor is linked to long-range changes in higher-order chromatin organization and epigenetic dysregulation in cancer. This may suggest that as well as drugs targeting histone modifications, it will be valuable to investigate the inhibition of protein complexes involved in chromatin folding in cancer cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0719-9) contains supplementary material, which is available to authorized users.     

Comparison of the three-dimensional organization of sperm and fibroblast genomes using the Hi-C approach

BackgroundThe three-dimensional organization of the genome is tightly connected to its biological function. The Hi-C approach was recently introduced as a method that can be used to identify higher-order chromatin interactions genome-wide. The aim of this study was to determine genome-wide chromatin interaction frequencies using the Hi-C approach in mouse sperm cells and embryonic fibroblasts.ResultsThe obtained data demonstrate that the three-dimensional genome organizations of sperm and fibroblast cells show a high degree of similarity both with each other and with the previously described mouse embryonic stem cells. Both A- and B-compartments and topologically associated domains are present in spermatozoa and fibroblasts. Nevertheless, sperm cells and fibroblasts exhibit statistically significant differences between each other in the contact probabilities of defined loci. Tight packaging of the sperm genome results in an enrichment of long-range contacts compared with the fibroblasts. However, only 30% of the differences in the number of contacts are based on differences in the densities of their genome packages; the main source of the differences is the gain or loss of contacts that are specific for defined genome regions. We find that the dependence of the contact probability on genomic distance for sperm is close to the dependence predicted for the fractal globular folding of chromatin.ConclusionsOverall, we can conclude that the three-dimensional structure of the genome is passed through generations without being dramatically changed in sperm cells.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0642-0) contains supplementary material, which is available to authorized users.     

Large chromatin domains in pluripotent and differentiated cells

Pluripotent stem cells are able to proliferate unlimitedly and to generate all somatic cell types, thus holding a great promise in medical applications. Epigenetic modifications are believed to play crucial roles in regulating pluripotency and differentiation. Recent genome-wide studies on mammalian systems have revealed several types of large chromatin domains which are associated with higher-order organization of the genome. The elucidation of genomic distribution and dynamics of these domains have shed light on the mechanisms underling pluripotency and lineage commitment.     

Archaea: The Final Frontier of Chromatin

The three domains of life employ various strategies to organize their genomes. Archaea utilize features similar to those found in both eukaryotic and bacterial chromatin to organize their DNA. In this review, we discuss the current state of research regarding the structure–function relationships of several archaeal chromatin proteins (histones, Alba, Cren7, and Sul7d). We address individual structures as well as inferred models for higher-order chromatin formation. Each protein introduces a unique phenotype to chromatin organization, and these structures are put into the context of in vivo and in vitro data. We close by discussing the present gaps in knowledge that are preventing further studies of the organization of archaeal chromatin, on both the organismal and domain level.     

Predicting chromatin organization using histone marks

Genome-wide mapping of three dimensional chromatin organization is an important yet technically challenging task. To aid experimental effort and to understand the determinants of long-range chromatin interactions, we have developed a computational model integrating Hi-C and histone mark ChIP-seq data to predict two important features of chromatin organization: chromatin interaction hubs and topologically associated domain (TAD) boundaries. Our model accurately and robustly predicts these features across datasets and cell types. Cell-type specific histone mark information is required for prediction of chromatin interaction hubs but not for TAD boundaries. Our predictions provide a useful guide for the exploration of chromatin organization.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0740-z) contains supplementary material, which is available to authorized users.     

Whole genome bisulfite sequencing of cell-free DNA and its cellular contributors uncovers placenta hypomethylated domains

BackgroundCirculating cell-free fetal DNA has enabled non-invasive prenatal fetal aneuploidy testing without direct discrimination of the maternal and fetal DNA. Testing may be improved by specifically enriching the sample material for fetal DNA. DNA methylation may allow for such a separation of DNA; however, this depends on knowledge of the methylomes of circulating cell-free DNA and its cellular contributors.ResultsWe perform whole genome bisulfite sequencing on a set of unmatched samples including circulating cell-free DNA from non-pregnant and pregnant female donors and genomic DNA from maternal buffy coat and placenta samples. We find CpG cytosines within longer fragments are more likely to be methylated. Comparison of the methylomes of placenta and non-pregnant circulating cell-free DNA reveal many of the 51,259 identified differentially methylated regions are located in domains exhibiting consistent placenta hypomethylation across millions of consecutive bases. We find these placenta hypomethylated domains are consistently located within regions exhibiting low CpG and gene density. Differentially methylated regions identified when comparing placenta to non-pregnant circulating cell-free DNA are recapitulated in pregnant circulating cell-free DNA, confirming the ability to detect differential methylation in circulating cell-free DNA mixtures.ConclusionsWe generate methylome maps for four sample types at single-base resolution, identify a link between DNA methylation and fragment length in circulating cell-free DNA, identify differentially methylated regions between sample groups, and uncover the presence of megabase-size placenta hypomethylated domains.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0645-x) contains supplementary material, which is available to authorized users.     

Distinct polymer physics principles govern chromatin dynamics in mouse and Drosophila topological domains

Background

In higher eukaryotes, the genome is partitioned into large "Topologically Associating Domains" (TADs) in which the chromatin displays favoured long-range contacts. While a crumpled/fractal globule organization has received experimental supports at higher-order levels, the organization principles that govern chromatin dynamics within these TADs remain unclear. Using simple polymer models, we previously showed that, in mouse liver cells, gene-rich domains tend to adopt a statistical helix shape when no significant locus-specific interaction takes place.

Results

Here, we use data from diverse 3C-derived methods to explore chromatin dynamics within mouse and Drosophila TADs. In mouse Embryonic Stem Cells (mESC), that possess large TADs (median size of 840 kb), we show that the statistical helix model, but not globule models, is relevant not only in gene-rich TADs, but also in gene-poor and gene-desert TADs. Interestingly, this statistical helix organization is considerably relaxed in mESC compared to liver cells, indicating that the impact of the constraints responsible for this organization is weaker in pluripotent cells. Finally, depletion of histone H1 in mESC alters local chromatin flexibility but not the statistical helix organization. In Drosophila, which possesses TADs of smaller sizes (median size of 70 kb), we show that, while chromatin compaction and flexibility are finely tuned according to the epigenetic landscape, chromatin dynamics within TADs is generally compatible with an unconstrained polymer configuration.

Conclusions

Models issued from polymer physics can accurately describe the organization principles governing chromatin dynamics in both mouse and Drosophila TADs. However, constraints applied on this dynamics within mammalian TADs have a peculiar impact resulting in a statistical helix organization.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1786-8) contains supplementary material, which is available to authorized users.     

Relationship Between Chromatin High-order Folding and Nucleosomal Linker Twist in Nuclei of Human HeLa S3 Cells

Abstract

We have used the intercalative agent ethidium bromide to examine the association between chromatin high-order folding and the twist of internucleosomal DNA regions. The analysis was carried out on intact nuclei isolated from human HeLa S3 cells. Our data shows that alterations in the nucleosomal linker twist significantly influence the way in which a chain of nucleosomes folds to form different higher-order structures. The assay used allowed us to identify the existence of two chromatin fractions differing in their extent of high-order folding. We have also found that active gene sequences are preferentially associated with the chromatin fraction corresponding to the more extended conformation. A model is proposed to account for the effect of variations in the nucleosome linker twist on the state of chromatin folding.     

Identifying multi-locus chromatin contacts in human cells using tethered multiple 3C

Background

Several recently developed experimental methods, each an extension of the chromatin conformation capture (3C) assay, have enabled the genome-wide profiling of chromatin contacts between pairs of genomic loci in 3D. Especially in complex eukaryotes, data generated by these methods, coupled with other genome-wide datasets, demonstrated that non-random chromatin folding correlates strongly with cellular processes such as gene expression and DNA replication.

Results

We describe a genome architecture assay, tethered multiple 3C (TM3C), that maps genome-wide chromatin contacts via a simple protocol of restriction enzyme digestion and religation of fragments upon agarose gel beads followed by paired-end sequencing. In addition to identifying contacts between pairs of loci, TM3C enables identification of contacts among more than two loci simultaneously. We use TM3C to assay the genome architectures of two human cell lines: KBM7, a near-haploid chronic leukemia cell line, and NHEK, a normal diploid human epidermal keratinocyte cell line. We confirm that the contact frequency maps produced by TM3C exhibit features characteristic of existing genome architecture datasets, including the expected scaling of contact probabilities with genomic distance, megabase scale chromosomal compartments and sub-megabase scale topological domains. We also confirm that TM3C captures several known cell type-specific contacts, ploidy shifts and translocations, such as Philadelphia chromosome formation (Ph+) in KBM7. We confirm a subset of the triple contacts involving the IGF2-H19 imprinting control region (ICR) using PCR analysis for KBM7 cells. Our genome-wide analysis of pairwise and triple contacts demonstrates their preference for linking open chromatin regions to each other and for linking regions with higher numbers of DNase hypersensitive sites (DHSs) to each other. For near-haploid KBM7 cells, we infer whole genome 3D models that exhibit clustering of small chromosomes with each other and large chromosomes with each other, consistent with previous studies of the genome architectures of other human cell lines.

Conclusion

TM3C is a simple protocol for ascertaining genome architecture and can be used to identify simultaneous contacts among three or four loci. Application of TM3C to a near-haploid human cell line revealed large-scale features of chromosomal organization and multi-way chromatin contacts that preferentially link regions of open chromatin.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1236-7) contains supplementary material, which is available to authorized users.     

Comparison of Hi-C results using in-solution versus in-nucleus ligation

BackgroundChromosome conformation capture and various derivative methods such as 4C, 5C and Hi-C have emerged as standard tools to analyze the three-dimensional organization of the genome in the nucleus. These methods employ ligation of diluted cross-linked chromatin complexes, intended to favor proximity-dependent, intra-complex ligation. During development of single-cell Hi-C, we devised an alternative Hi-C protocol with ligation in preserved nuclei rather than in solution. Here we directly compare Hi-C methods employing in-nucleus ligation with the standard in-solution ligation.ResultsWe show in-nucleus ligation results in consistently lower levels of inter-chromosomal contacts. Through chromatin mixing experiments we show that a significantly large fraction of inter-chromosomal contacts are the result of spurious ligation events formed during in-solution ligation. In-nucleus ligation significantly reduces this source of experimental noise, and results in improved reproducibility between replicates. We also find that in-nucleus ligation eliminates restriction fragment length bias found with in-solution ligation. These improvements result in greater reproducibility of long-range intra-chromosomal and inter-chromosomal contacts, as well as enhanced detection of structural features such as topologically associated domain boundaries.ConclusionsWe conclude that in-nucleus ligation captures chromatin interactions more consistently over a wider range of distances, and significantly reduces both experimental noise and bias. In-nucleus ligation creates higher quality Hi-C libraries while simplifying the experimental procedure. We suggest that the entire range of 3C applications are likely to show similar benefits from in-nucleus ligation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13059-015-0753-7) contains supplementary material, which is available to authorized users.     

Epigenome overlap measure (EPOM) for comparing tissue/cell types based on chromatin states

Background

The dynamics of epigenomic marks in their relevant chromatin states regulate distinct gene expression patterns, biological functions and phenotypic variations in biological processes. The availability of high-throughput epigenomic data generated by next-generation sequencing technologies allows a data-driven approach to evaluate the similarities and differences of diverse tissue and cell types in terms of epigenomic features. While ChromImpute has allowed for the imputation of large-scale epigenomic information to yield more robust data to capture meaningful relationships between biological samples, widely used methods such as hierarchical clustering and correlation analysis cannot adequately utilize epigenomic data to accurately reveal the distinction and grouping of different tissue and cell types.

Methods

We utilize a three-step testing procedure–ANOVA, t test and overlap test to identify tissue/cell-type- associated enhancers and promoters and to calculate a newly defined Epigenomic Overlap Measure (EPOM). EPOM results in a clear correspondence map of biological samples from different tissue and cell types through comparison of epigenomic marks evaluated in their relevant chromatin states.

Results

Correspondence maps by EPOM show strong capability in distinguishing and grouping different tissue and cell types and reveal biologically meaningful similarities between Heart and Muscle, Blood & T-cell and HSC & B-cell, Brain and Neurosphere, etc. The gene ontology enrichment analysis both supports and explains the discoveries made by EPOM and suggests that the associated enhancers and promoters demonstrate distinguishable functions across tissue and cell types. Moreover, the tissue/cell-type-associated enhancers and promoters show enrichment in the disease-related SNPs that are also associated with the corresponding tissue or cell types. This agreement suggests the potential of identifying causal genetic variants relevant to cell-type-specific diseases from our identified associated enhancers and promoters.

Conclusions

The proposed EPOM measure demonstrates superior capability in grouping and finding a clear correspondence map of biological samples from different tissue and cell types. The identified associated enhancers and promoters provide a comprehensive catalog to study distinct biological processes and disease variants in different tissue and cell types. Our results also find that the associated promoters exhibit more cell-type-specific functions than the associated enhancers do, suggesting that the non-associated promoters have more housekeeping functions than the non-associated enhancers.