Chromatin architecture reorganization during stem cell differentiation

Chromatin structure-based prediction of recurrent noncoding mutations in cancer

Recurrence is a hallmark of cancer-driving mutations. Recurrent mutations can arise at the same site or affect the same gene at different sites. Here we identified a set of mutations arising in individual samples and altering different cis-regulatory elements that converge on a common gene via chromatin interactions. The mutations and genes identified in this fashion showed strong relevance to cancer, in contrast to noncoding mutations with site-specific recurrence only. We developed a prediction method that identifies potentially recurrent mutations on the basis of the features shared by mutations whose recurrence is observed in a given cohort. Our method was capable of accurately predicting recurrent mutations at the level of target genes but not mutations recurring at the same site. We experimentally validated predicted mutations in distal regulatory regions of the TERT gene. In conclusion, we propose a novel approach to discovering potential cancer-driving mutations in noncoding regions.

Statistical challenges in analyzing methylation and long-range chromosomal interaction data

With the rapid development of high throughput technologies such as array and next generation sequencing (NGS), genome-wide, nucleotide-resolution epigenomic data are increasingly available. In recent years, there has been particular interest in data on DNA methylation and 3-dimensional (3D) chromosomal organization, which are believed to hold keys to understand biological mechanisms, such as transcription regulation, that are closely linked to human health and diseases. However, small sample size, complicated correlation structure, substantial noise, biases, and uncertainties, all present difficulties for performing statistical inference. In this review, we present an overview of the new technologies that are frequently utilized in studying DNA methylation and 3D chromosomal organization. We focus on reviewing recent developments in statistical methodologies designed for better interrogating epigenomic data, pointing out statistical challenges facing the field whenever appropriate.

Systematic mapping of functional enhancer-promoter connections with CRISPR interference

Gene expression in mammals is regulated by noncoding elements that can affect physiology and disease, yet the functions and target genes of most noncoding elements remain unknown. We present a high-throughput approach that uses clustered regularly interspaced short palindromic repeats (CRISPR) interference (CRISPRi) to discover regulatory elements and identify their target genes. We assess >1 megabase of sequence in the vicinity of two essential transcription factors, MYC and GATA1, and identify nine distal enhancers that control gene expression and cellular proliferation. Quantitative features of chromatin state and chromosome conformation distinguish the seven enhancers that regulate MYC from other elements that do not, suggesting a strategy for predicting enhancer-promoter connectivity. This CRISPRi-based approach can be applied to dissect transcriptional networks and interpret the contributions of noncoding genetic variation to human disease.

Genome-wide association study for acute otitis media in children identifies FNDC1 as disease contributing gene

Acute otitis media (AOM) is among the most common pediatric diseases, and the most frequent reason for antibiotic treatment in children. Risk of AOM is dependent on environmental and host factors, as well as a significant genetic component. We identify genome-wide significance at a locus on 6q25.3 (rs2932989, P_meta=2.15 × 10⁻⁰⁹), and show that the associated variants are correlated with the methylation status of the FNDC1 gene (cg05678571, P=1.43 × 10⁻⁰⁶), and further show it is an eQTL for FNDC1 (P=9.3 × 10⁻⁰⁵). The mouse homologue, Fndc1, is expressed in middle ear tissue and its expression is upregulated upon lipopolysaccharide treatment. In this first GWAS of AOM and the largest OM genetic study to date, we identify the first genome-wide significant locus associated with AOM.

Acute otitis media (AOM) is an acute infection of middle ear mucosa and among the most common pediatric diseases. Here, the authors performed a genome-wide association study to link a variant in the FNDC1 locus on 6q25.3 and differential methylation status of the FNDC1 gene with predisposition to AOM.

The Chromatin Landscape of Cellular Senescence

Cellular senescence, an irreversible growth arrest triggered by a variety of stressors, plays important roles in normal physiology and tumor suppression, but accumulation of senescent cells with age contributes to the functional decline of tissues. Senescent cells undergo dramatic alterations to their chromatin landscape that affect genome accessibility and their transcriptional program. These include the loss of DNA-nuclear lamina interactions, the distension of centromeres, and changes in chromatin composition that can lead to the activation of retrotransposons. Here we discuss these findings, as well as recent advances in microscopy and genomics that have revealed the importance of the higher-order spatial organization of the genome in defining and maintaining the senescent state.

Multi-tiered Reorganization of the Genome during B Cell Affinity Maturation Anchored by a Germinal Center-Specific Locus Control Region

During the humoral immune response, B cells undergo a dramatic change in phenotype to enable antibody affinity maturation in germinal centers (GCs). Using genome-wide chromosomal conformation capture (Hi-C), we found that GC B cells undergo massive reorganization of the genomic architecture that encodes the GC B cell transcriptome. Coordinate expression of genes that specify the GC B cell phenotype-most prominently BCL6-was achieved through a multilayered chromatin reorganization process involving (1) increased promoter connectivity, (2) formation of enhancer networks, (3) 5' to 3' gene looping, and (4) merging of gene neighborhoods that share active epigenetic marks. BCL6 was an anchor point for the formation of GC-specific gene and enhancer loops on chromosome 3. Deletion of a GC-specific, highly interactive locus control region upstream of Bcl6 abrogated GC formation in mice. Thus, large-scale and multi-tiered genomic three-dimensional reorganization is required for coordinate expression of phenotype-driving gene sets that determine the unique characteristics of GC B cells.

Spatio-temporal re-organization of replication foci accompanies replication domain consolidation during human pluripotent stem cell lineage specification

Lineage specification of both mouse and human pluripotent stem cells (PSCs) is accompanied by spatial consolidation of chromosome domains and temporal consolidation of their replication timing. Replication timing and chromatin organization are both established during G1 phase at the timing decision point (TDP). Here, we have developed live cell imaging tools to track spatio-temporal replication domain consolidation during differentiation. First, we demonstrate that the fluorescence ubiquitination cell cycle indicator (Fucci) system is incapable of demarcating G1/S or G2/M cell cycle transitions. Instead, we employ a combination of fluorescent PCNA to monitor S phase progression, cytokinesis to demarcate mitosis, and fluorescent nucleotides to label early and late replication foci and track their 3D organization into sub-nuclear chromatin compartments throughout all cell cycle transitions. We find that, as human PSCs differentiate, the length of S phase devoted to replication of spatially clustered replication foci increases, coincident with global compartmentalization of domains into temporally clustered blocks of chromatin. Importantly, re-localization and anchorage of domains was completed prior to the onset of S phase, even in the context of an abbreviated PSC G1 phase. This approach can also be employed to investigate cell fate transitions in single PSCs, which could be seen to differentiate preferentially from G1 phase. Together, our results establish real-time, live-cell imaging methods for tracking cell cycle transitions during human PSC differentiation that can be applied to study chromosome domain consolidation and other aspects of lineage specification.

Capture of associated targets on chromatin links long-distance chromatin looping to transcriptional coordination

Here we describe a sensitive and novel method of identifying endogenous DNA-DNA interactions. Capture of Associated Targets on CHromatin (CATCH) uses efficient capture and enrichment of specific genomic loci of interest through hybridization and subsequent purification via complementary biotinylated oligonucleotide. The CATCH assay requires no enzymatic digestion or ligation, requires little starting material, provides high-quality data, has excellent reproducibility and is completed in less than 24 h. Efficacy is demonstrated through capture of three disparate loci, which demonstrate unique subsets of long-distance chromatin interactions enriched for both enhancer marks and oestrogen receptor-binding sites. In each experiment, CATCH-seq peaks representing long-distance chromatin interactions were centred near the TSS of genes, and, critically, the genes identified as physically interacting are shown to be transcriptionally coexpressed. These interactions could potentially create transcriptional hubs for the regulation of gene expression programmes.

Binding of nuclear factor κB to noncanonical consensus sites reveals its multimodal role during the early inflammatory response

Mammalian cells have developed intricate mechanisms to interpret, integrate, and respond to extracellular stimuli. For example, tumor necrosis factor (TNF) rapidly activates proinflammatory genes, but our understanding of how this occurs against the ongoing transcriptional program of the cell is far from complete. Here, we monitor the early phase of this cascade at high spatiotemporal resolution in TNF-stimulated human endothelial cells. NF-κB, the transcription factor complex driving the response, interferes with the regulatory machinery by binding active enhancers already in interaction with gene promoters. Notably, >50% of these enhancers do not encode canonical NF-κB binding motifs. Using a combination of genomics tools, we find that binding site selection plays a key role in NF-κΒ–mediated transcriptional activation and repression. We demonstrate the latter by describing the synergy between NF-κΒ and the corepressor JDP2. Finally, detailed analysis of a 2.8-Mbp locus using sub-kbp-resolution targeted chromatin conformation capture and genome editing uncovers how NF-κΒ that has just entered the nucleus exploits pre-existing chromatin looping to exert its multimodal role. This work highlights the involvement of topology in cis-regulatory element function during acute transcriptional responses, where primary DNA sequence and its higher-order structure constitute a regulatory context leading to either gene activation or repression.

Aberrant Epigenetic Regulation in the Pathogenesis of Systemic Lupus Erythematosus and Its Implication in Precision Medicine

Great progress has been made in the last decades in understanding the complex immune dysregulation in systemic lupus erythematosus (SLE), yet the efforts to pursue an effective treatment of SLE proved to be futile. The pathoetiology of SLE involves extremely complicated and multifactorial interaction among various genetic and epigenetic factors. Multiple gene loci predispose to disease susceptibility, and the interaction with epigenetic modifications mediated through sex, hormones, and the hypothalamo-pituitary-adrenal axis complicates susceptibility and manifestations of this disease. Finally, certain environmental and psychological factors probably trigger the disease via epigenetic mechanisms. In this review, we summarize and discuss recent epigenetic studies of SLE and suggest a personalized approach to the dissection of disease onset and therapy or precision medicine. We speculate that in the future, precision medicine based on epigenetic and genetic information could help guide more effective targeted therapeutic intervention.

Nuclear organization and 3D chromatin architecture in cognition and neuropsychiatric disorders

The current view of neuroplasticity depicts the changes in the strength and number of synaptic connections as the main physical substrate for behavioral adaptation to new experiences in a changing environment. Although transcriptional regulation is known to play a role in these synaptic changes, the specific contribution of activity-induced changes to both the structure of the nucleus and the organization of the genome remains insufficiently characterized. Increasing evidence indicates that plasticity-related genes may work in coordination and share architectural and transcriptional machinery within discrete genomic foci. Here we review the molecular and cellular mechanisms through which neuronal nuclei structurally adapt to stimuli and discuss how the perturbation of these mechanisms can trigger behavioral malfunction.

In silico identification of enhancers on the basis of a combination of transcription factor binding motif occurrences

Enhancers interact with gene promoters and form chromatin looping structures that serve important functions in various biological processes, such as the regulation of gene transcription and cell differentiation. However, enhancers are difficult to identify because they generally do not have fixed positions or consensus sequence features, and biological experiments for enhancer identification are costly in terms of labor and expense. In this work, several models were built by using various sequence-based feature sets and their combinations for enhancer prediction. The selected features derived from a recursive feature elimination method showed that the model using a combination of 141 transcription factor binding motif occurrences from 1,422 transcription factor position weight matrices achieved a favorably high prediction accuracy superior to that of other reported methods. The models demonstrated good prediction accuracy for different enhancer datasets obtained from different cell lines/tissues. In addition, prediction accuracy was further improved by integration of chromatin state features. Our method is complementary to wet-lab experimental methods and provides an additional method to identify enhancers.

CTCF and Cohesin in Genome Folding and Transcriptional Gene Regulation

Genome function, replication, integrity, and propagation rely on the dynamic structural organization of chromosomes during the cell cycle. Genome folding in interphase provides regulatory segmentation for appropriate transcriptional control, facilitates ordered genome replication, and contributes to genome integrity by limiting illegitimate recombination. Here, we review recent high-resolution chromosome conformation capture and functional studies that have informed models of the spatial and regulatory compartmentalization of mammalian genomes, and discuss mechanistic models for how CTCF and cohesin control the functional architecture of mammalian chromosomes.

Identification of epigenetic factors regulating the mesenchyme to epithelium transition by RNA interference screening in breast cancer cells

Background

In breast cancer, the epithelial to mesenchyme transition (EMT) is associated to tumour dissemination, drug resistance and high relapse risks. It is partly controlled by epigenetic modifications such as histone acetylation and methylation. The identification of genes involved in these reversible modifications represents an interesting therapeutic strategy to fight metastatic disease by inducing mesenchymal cell differentiation to an epithelial phenotype.

Methods

We designed a siRNA library based on chromatin modification-related to functional domains and screened it in the mesenchymal breast cancer cell line MDA-MB-231. The mesenchyme to epithelium transition (MET) activation was studied by following human E-CADHERIN (E-CAD) induction, a specific MET marker, and cell morphology. Candidate genes were validated by studying the expression of several differential marker genes and their impact on cell migration.

Results

The screen led to the identification of 70 gene candidates among which some are described to be, directly or indirectly, involved in EMT like ZEB1, G9a, SMAD5 and SMARCD3. We also identified the DOT1L as involved in EMT regulation in MDA-MB-231. Moreover, for the first time, KAT5 gene was linked to the maintenance of the mesenchymal phenotype.

Conclusions

A multi-parametric RNAi screening approach was developed to identify new EMT regulators such as KAT5 in the triple negative breast cancer cell line MDA-MB-231.

Electronic supplementary material

The online version of this article (doi:10.1186/s12885-016-2683-5) contains supplementary material, which is available to authorized users.

Single-cell pluripotency regulatory networks

Pluripotent stem cells (PSCs) are a popular model system for investigating development, tissue regeneration, and repair. Although much is known about the molecular mechanisms that regulate the balance between self-renewal and lineage commitment in PSCs, the spatiotemporal integration of responsive signaling pathways with core transcriptional regulatory networks are complex and only partially understood. Moreover, measurements made on populations of cells reveal only average properties of the underlying regulatory networks, obscuring their fine detail. Here, we discuss the reconstruction of regulatory networks in individual cells using novel single-cell transcriptomics and proteomics, in order to expand our understanding of the molecular basis of pluripotency, including the role of cell-cell variability within PSC populations, and ways in which networks may be controlled in order to reliably manipulate cell behavior.

From single genes to entire genomes: the search for a function of nuclear organization

The existence of different domains within the nucleus has been clear from the time, in the late 1920s, that heterochromatin and euchromatin were discovered. The observation that heterochromatin is less transcribed than euchromatin suggested that microscopically identifiable structures might correspond to functionally different domains of the nucleus. Until 15 years ago, studies linking gene expression and subnuclear localization were limited to a few genes. As we discuss in this Review, new genome-wide techniques have now radically changed the way nuclear organization is analyzed. These have provided a much more detailed view of functional nuclear architecture, leading to the emergence of a number of new paradigms of chromatin folding and how this folding evolves during development.

SMARCA4 regulates gene expression and higher-order chromatin structure in proliferating mammary epithelial cells

The packaging of DNA into chromatin plays an important role in transcriptional regulation and nuclear processes. Brahma-related gene-1 SMARCA4 (also known as BRG1), the essential ATPase subunit of the mammalian SWI/SNF chromatin remodeling complex, uses the energy from ATP hydrolysis to disrupt nucleosomes at target regions. Although the transcriptional role of SMARCA4 at gene promoters is well-studied, less is known about its role in higher-order genome organization. SMARCA4 knockdown in human mammary epithelial MCF-10A cells resulted in 176 up-regulated genes, including many related to lipid and calcium metabolism, and 1292 down-regulated genes, some of which encode extracellular matrix (ECM) components that can exert mechanical forces and affect nuclear structure. ChIP-seq analysis of SMARCA4 localization and SMARCA4-bound super-enhancers demonstrated extensive binding at intergenic regions. Furthermore, Hi-C analysis showed extensive SMARCA4-mediated alterations in higher-order genome organization at multiple resolutions. First, SMARCA4 knockdown resulted in clustering of intra- and inter-subtelomeric regions, demonstrating a novel role for SMARCA4 in telomere organization. SMARCA4 binding was enriched at topologically associating domain (TAD) boundaries, and SMARCA4 knockdown resulted in weakening of TAD boundary strength. Taken together, these findings provide a dynamic view of SMARCA4-dependent changes in higher-order chromatin organization and gene expression, identifying SMARCA4 as a novel component of chromatin organization.

Combinatorial control of plant gene expression

Combinatorial gene regulation provides a mechanism by which relatively small numbers of transcription factors can control the expression of a much larger number of genes with finely tuned temporal and spatial patterns. This is achieved by transcription factors assembling into complexes in a combinatorial fashion, exponentially increasing the number of genes that they can target. Such an arrangement also increases the specificity and affinity for the cis-regulatory sequences required for accurate target gene expression. Superimposed on this transcription factor combinatorial arrangement is the increasing realization that histone modification marks expand the regulatory information, which is interpreted by histone readers and writers that are part of the regulatory apparatus. Here, we review the progress in these areas from the perspective of plant combinatorial gene regulation, providing examples of different regulatory solutions and comparing them to other metazoans. This article is part of a Special Issue entitled: Plant Gene Regulatory Mechanisms and Networks, edited by Dr. Erich Grotewold and Dr. Nathan Springer.

Integrating epigenomic data and 3D genomic structure with a new measure of chromatin assortativity

BACKGROUND

Network analysis is a powerful way of modeling chromatin interactions. Assortativity is a network property used in social sciences to identify factors affecting how people establish social ties. We propose a new approach, using chromatin assortativity, to integrate the epigenomic landscape of a specific cell type with its chromatin interaction network and thus investigate which proteins or chromatin marks mediate genomic contacts.

RESULTS

We use high-resolution promoter capture Hi-C and Hi-Cap data as well as ChIA-PET data from mouse embryonic stem cells to investigate promoter-centered chromatin interaction networks and calculate the presence of specific epigenomic features in the chromatin fragments constituting the nodes of the network. We estimate the association of these features with the topology of four chromatin interaction networks and identify features localized in connected areas of the network. Polycomb group proteins and associated histone marks are the features with the highest chromatin assortativity in promoter-centered networks. We then ask which features distinguish contacts amongst promoters from contacts between promoters and other genomic elements. We observe higher chromatin assortativity of the actively elongating form of RNA polymerase 2 (RNAPII) compared with inactive forms only in interactions between promoters and other elements.

CONCLUSIONS

Contacts among promoters and between promoters and other elements have different characteristic epigenomic features. We identify a possible role for the elongating form of RNAPII in mediating interactions among promoters, enhancers, and transcribed gene bodies. Our approach facilitates the study of multiple genome-wide epigenomic profiles, considering network topology and allowing the comparison of chromatin interaction networks.

Analysis of chromosomal aberrations and recombination by allelic bias in RNA-Seq

Genomic instability has profound effects on cellular phenotypes. Studies have shown that pluripotent cells with abnormal karyotypes may grow faster, differentiate less and become more resistance to apoptosis. Previously, we showed that microarray gene expression profiles can be utilized for the analysis of chromosomal aberrations by comparing gene expression levels between normal and aneuploid samples. Here we adopted this method for RNA-Seq data and present eSNP-Karyotyping for the detection of chromosomal aberrations, based on measuring the ratio of expression between the two alleles. We demonstrate its ability to detect chromosomal gains and losses in pluripotent cells and their derivatives, as well as meiotic recombination patterns. This method is advantageous since it does not require matched diploid samples for comparison, is less sensitive to global expression changes caused by the aberration and utilizes already available gene expression profiles to determine chromosomal aberrations.

Chromosomal aberrations can be detected by global gene expression analysis. Here, the authors report eSNP-Karyotyping, a new method that can detect chromosomal aberrations by measuring the ratio of expression between the two alleles without comparison to a matched diploid sample.

UMI-4C for quantitative and targeted chromosomal contact profiling

We developed a targeted chromosome conformation capture (4C) approach that uses unique molecular identifiers (UMIs) to derive high-complexity quantitative chromosome contact profiles with controlled signal-to-noise ratios. UMI-4C detects chromosomal interactions with improved sensitivity and specificity, and it can easily be multiplexed to allow robust comparison of contact distributions between loci and conditions. This approach may open the way to the incorporation of contact distributions into quantitative models of gene regulation.

The Chromatin Signature of Pluripotency: Establishment and Maintenance

The revolutionary discovery that somatic cells can be reprogrammed by a defined set transcription factors to induced pluripotent stem cells (iPSCs) changed dramatically the way we perceive cell fate determination. Importantly, iPSCs, similar to embryo-derived stem cells (ESCs), are characterized by a remarkable developmental plasticity and the capacity to self-renew "indefinitely" under appropriate culture conditions, opening new avenues for personalized therapy and disease modeling. Elucidating the molecular mechanisms that maintain, induce, or alter stem cell identity is crucial for a deeper understanding of cell fate determination and potential translational applications. Intense research over the last 10 years exploiting technological advances in epigenomics and genome editing has unraveled many of the mysteries of pluripotent identity enabling novel and efficient ways to manipulate it for biomedical purposes. In this review, we focus on the chromatin and epigenetic characteristics that distinguish stem cells from somatic cells and their dynamic changes during differentiation and reprogramming.

Gene regulation at a distance: From remote enhancers to 3D regulatory ensembles

Large-scale identification of elements associated with gene expression revealed that many of them are located extremely far from gene transcriptional start sites. We review here the growing evidence that show that distal cis-acting elements provide key instructions to genes, as genetic variation affecting them is growingly identified as an importance source of phenotypic diversity and disease. We discuss the different mechanisms that allow these elements to exert their regulatory functions, in a robust and specific manner, despite the large genomic distances separating them from their target genes. We particularly focus on the role of the structural organization of the genome in guiding such regulatory interactions.

SNPsplit: Allele-specific splitting of alignments between genomes with known SNP genotypes

Sequencing reads overlapping polymorphic sites in diploid mammalian genomes may be assigned to one allele or the other. This holds the potential to detect gene expression, chromatin modifications, DNA methylation or nuclear interactions in an allele-specific fashion. SNPsplit is an allele-specific alignment sorter designed to read files in SAM/BAM format and determine the allelic origin of reads or read-pairs that cover known single nucleotide polymorphic (SNP) positions. For this to work libraries must have been aligned to a genome in which all known SNP positions were masked with the ambiguity base 'N' and aligned using a suitable mapping program such as Bowtie2, TopHat, STAR, HISAT2, HiCUP or Bismark. SNPsplit also provides an automated solution to generate N-masked reference genomes for hybrid mouse strains based on the variant call information provided by the Mouse Genomes Project. The unique ability of SNPsplit to work with various different kinds of sequencing data including RNA-Seq, ChIP-Seq, Bisulfite-Seq or Hi-C opens new avenues for the integrative exploration of allele-specific data.

SNPsplit: Allele-specific splitting of alignments between genomes with known SNP genotypes

Sequencing reads overlapping polymorphic sites in diploid mammalian genomes may be assigned to one allele or the other. This holds the potential to detect gene expression, chromatin modifications, DNA methylation or nuclear interactions in an allele-specific fashion. SNPsplit is an allele-specific alignment sorter designed to read files in SAM/BAM format and determine the allelic origin of reads or read-pairs that cover known single nucleotide polymorphic (SNP) positions. For this to work libraries must have been aligned to a genome in which all known SNP positions were masked with the ambiguity base 'N' and aligned using a suitable mapping program such as Bowtie2, TopHat, STAR, HISAT2, HiCUP or Bismark. SNPsplit also provides an automated solution to generate N-masked reference genomes for hybrid mouse strains based on the variant call information provided by the Mouse Genomes Project. The unique ability of SNPsplit to work with various different kinds of sequencing data including RNA-Seq, ChIP-Seq, Bisulfite-Seq or Hi-C opens new avenues for the integrative exploration of allele-specific data.

Ten years of induced pluripotency: from basic mechanisms to therapeutic applications

Ten years ago, the discovery that mature somatic cells could be reprogrammed into induced pluripotent stem cells (iPSCs) redefined the stem cell field and brought about a wealth of opportunities for both basic research and clinical applications. To celebrate the tenth anniversary of the discovery, the International Society for Stem Cell Research (ISSCR) and Center for iPS Cell Research and Application (CiRA), Kyoto University, together held the symposium ‘Pluripotency: From Basic Science to Therapeutic Applications’ in Kyoto, Japan. The three days of lectures examined both the mechanisms and therapeutic applications of iPSC reprogramming. Here we summarize the main findings reported, which are testament to how far the field has come in only a decade, as well as the enormous potential that iPSCs hold for the future.

On TADs and LADs: Spatial Control Over Gene Expression

The combinatorial action of transcription factors drives cell-type-specific gene expression patterns. However, transcription factor binding and gene regulation occur in the context of chromatin, which modulates DNA accessibility. High-resolution chromatin interaction maps have defined units of chromatin that are in spatial proximity, called topologically associated domains (TADs). TADs can be further classified based on expression activity, replication timing, or the histone marks or non-histone proteins associated with them. Independently, other chromatin domains have been defined by their likelihood to interact with non-DNA structures, such as the nuclear lamina. Lamina-associated domains (LADs) correlate with low gene expression and late replication timing. TADs and LADs have recently been evaluated with respect to cell-type-specific gene expression. The results shed light on the relevance of these forms of chromatin organization for transcriptional regulation, and address specifically how chromatin sequestration influences cell fate decisions during organismal development.

CRISPR/Cas9 genome editing throws descriptive 3-D genome folding studies for a loop

CRISPR/Cas9 genome editing studies have recently shed new light into the causal link between the linear DNA sequence and 3-D chromatin architecture. Here we describe current models for the folding of genomes into a nested hierarchy of higher-order structures and discuss new insights into the organizing principles governing genome folding at each length scale. WIREs Syst Biol Med 2016, 8:286-299. doi: 10.1002/wsbm.1338 For further resources related to this article, please visit the WIREs website.

CRISPR/Cas9 in Genome Editing and Beyond

The Cas9 protein (CRISPR-associated protein 9), derived from type II CRISPR (clustered regularly interspaced short palindromic repeats) bacterial immune systems, is emerging as a powerful tool for engineering the genome in diverse organisms. As an RNA-guided DNA endonuclease, Cas9 can be easily programmed to target new sites by altering its guide RNA sequence, and its development as a tool has made sequence-specific gene editing several magnitudes easier. The nuclease-deactivated form of Cas9 further provides a versatile RNA-guided DNA-targeting platform for regulating and imaging the genome, as well as for rewriting the epigenetic status, all in a sequence-specific manner. With all of these advances, we have just begun to explore the possible applications of Cas9 in biomedical research and therapeutics. In this review, we describe the current models of Cas9 function and the structural and biochemical studies that support it. We focus on the applications of Cas9 for genome editing, regulation, and imaging, discuss other possible applications and some technical considerations, and highlight the many advantages that CRISPR/Cas9 technology offers.

Nucleosome positioning changes during human embryonic stem cell differentiation

Nucleosomes are the basic unit of chromatin. Nucleosome positioning (NP) plays a key role in transcriptional regulation and other biological processes. To better understand NP we used MNase-seq to investigate changes that occur as human embryonic stem cells (hESCs) transition to nascent mesoderm and then to smooth muscle cells (SMCs). Compared to differentiated cell derivatives, nucleosome occupancy at promoters and other notable genic sites, such as exon/intron junctions and adjacent regions, in hESCs shows a stronger correlation with transcript abundance and is less influenced by sequence content. Upon hESC differentiation, genes being silenced, but not genes being activated, display a substantial change in nucleosome occupancy at their promoters. Genome-wide, we detected a shift of NP to regions of higher G+C content as hESCs differentiate to SMCs. Notably, genomic regions with higher nucleosome occupancy harbor twice as many G↔C changes but fewer than half A↔T changes, compared to regions with lower nucleosome occupancy. Finally, our analysis indicates that the hESC genome is not rearranged and has a sequence mutation rate resembling normal human genomes. Our study reveals another unique feature of hESC chromatin, and sheds light on the relationship between nucleosome occupancy and sequence G+C content.

Alterations to chromatin in intestinal macrophages link IL-10 deficiency to inappropriate inflammatory responses

Intestinal macrophages (IMs) are uniquely programmed to tolerate exposure to bacteria without mounting potent inflammatory responses. The cytokine IL-10 maintains the macrophage anti-inflammatory response such that loss of IL-10 results in chronic intestinal inflammation. To investigate how IL-10-deficiency alters IM programming and bacterial tolerance, we studied changes in chromatin accessibility in response to bacteria in macrophages from two distinct niches, the intestine and bone-marrow, from both wild-type and IL-10-deficient (Il10(-/-) ) mice. We identified chromatin accessibility changes associated with bacterial exposure and IL-10 deficiency in both bone marrow derived macrophages and IMs. Surprisingly, Il10(-/-) IMs adopted chromatin and gene expression patterns characteristic of an inflammatory response, even in the absence of bacteria. Further, when recombinant IL-10 was added to Il10(-/-) cells, it could not revert the chromatin landscape to a normal state. Our results demonstrate that IL-10 deficiency results in stable chromatin alterations in macrophages, even in the absence of bacteria. This supports a model in which IL-10-deficiency leads to chromatin alterations that contribute to a loss of IM tolerance to bacteria, which is a primary initiating event in chronic intestinal inflammation.

Mechanism of chromatin segregation to the nuclear periphery in C. elegans embryos

In eukaryotic organisms, gene regulation occurs in the context of chromatin. In the interphase nucleus, euchromatin and heterochromatin occupy distinct space during cell differentiation, with heterochromatin becoming enriched at the nuclear and nucleolar peripheries. This organization is thought to fine-tune gene expression. To elucidate the mechanisms that govern this level of genome organization, screens were carried out in C. elegans which monitored the loss of heterochromatin sequestration at the nuclear periphery. This led to the identification of a novel chromodomain protein, CEC-4 (Caenorhabditis elegans chromodomain protein 4) that mediates the anchoring of H3K9 methylation-bearing chromatin at the nuclear periphery in early to mid-stage embryos. Surprisingly, the loss of CEC-4 does not derepress genes found in heterochromatic domains, nor does it affect differentiation under standard laboratory conditions. On the other hand, CEC-4 contributes to the efficiency with which muscle differentiation is induced following ectopic expression of the master regulator, HLH-1. This is one of the first phenotypes specifically attributed to the ablation of heterochromatin anchoring.

A multi-task graph-clustering approach for chromosome conformation capture data sets identifies conserved modules of chromosomal interactions

Chromosome conformation capture methods are being increasingly used to study three-dimensional genome architecture in multiple cell types and species. An important challenge is to examine changes in three-dimensional architecture across cell types and species. We present Arboretum-Hi-C, a multi-task spectral clustering method, to identify common and context-specific aspects of genome architecture. Compared to standard clustering, Arboretum-Hi-C produced more biologically consistent patterns of conservation. Most clusters are conserved and enriched for either high- or low-activity genomic signals. Most genomic regions diverge between clusters with similar chromatin state except for a few that are associated with lamina-associated domains and open chromatin.

Nuclease Footprints in Sperm Project Past and Future Chromatin Regulatory Events

Nuclear remodeling to a condensed state is a hallmark of spermatogenesis. This is achieved by replacement of histones with protamines. Regions retaining nucleosomes may be of functional significance. To determine their potential roles, sperm from wild type and transgenic mice harboring a single copy insert of the human protamine cluster were subjected to Micrococcal Nuclease-seq. CENTIPEDE, a hierarchical Bayesian model, was used to identify multiple spatial patterns, "footprints", of MNase-seq reads along the sperm genome. Regions predicted by CENTIPEDE analysis to be bound by a regulatory factor in sperm were correlated with genomic landmarks and higher order chromatin structure datasets to identify potential roles for these factors in regulating either prior or post spermatogenic, i.e., early embryonic events. This approach linked robust endogenous protamine transcription and transgene suppression to its chromatin environment within topologically associated domains. Of the candidate enhancer-bound regulatory proteins, Ctcf, was associated with chromatin domain boundaries in testes and embryonic stem cells. The continuity of Ctcf binding through the murine germline may permit rapid reconstitution of chromatin organization following fertilization. This likely reflects its preparation for early zygotic genome activation and comparatively accelerated preimplantation embryonic development program observed in mouse as compared to human and bull.

Dynamic chromatin organization: Role in development and disease

The spatial organization of chromatin in the nucleus is important for proper regulation of gene expression. The cell-type specific transcription program is mainly controlled by distal regulatory elements, which can dynamically engage in long-range interactions with their target genes. These long-range interactions mostly occur within insulated genomic domains and are constrained by global organization of the chromatin, providing an extra layer of regulation. Genetic alterations can lead to disruption of spatial organization and consequently aberrant gene expression. In this review we will discuss the multiple layers of chromatin organization, how this organization changes during development and how its disruption can lead do aberrant gene expression and disease.

FastHiC: a fast and accurate algorithm to detect long-range chromosomal interactions from Hi-C data

MOTIVATION

How chromatin folds in three-dimensional (3D) space is closely related to transcription regulation. As powerful tools to study such 3D chromatin conformation, the recently developed Hi-C technologies enable a genome-wide measurement of pair-wise chromatin interaction. However, methods for the detection of biologically meaningful chromatin interactions, i.e. peak calling, from Hi-C data, are still under development. In our previous work, we have developed a novel hidden Markov random field (HMRF) based Bayesian method, which through explicitly modeling the non-negligible spatial dependency among adjacent pairs of loci manifesting in high resolution Hi-C data, achieves substantially improved robustness and enhanced statistical power in peak calling. Superior to peak callers that ignore spatial dependency both methodologically and in performance, our previous Bayesian framework suffers from heavy computational costs due to intensive computation incurred by modeling the correlated peak status of neighboring loci pairs and the inference of hidden dependency structure.

RESULTS

In this work, we have developed FastHiC, a novel approach based on simulated field approximation, which approximates the joint distribution of the hidden peak status by a set of independent random variables, leading to more tractable computation. Performance comparisons in real data analysis showed that FastHiC not only speeds up our original Bayesian method by more than five times, bus also achieves higher peak calling accuracy.

AVAILABILITY AND IMPLEMENTATION

FastHiC is freely accessible at:http://www.unc.edu/∼yunmli/FastHiC/ CONTACTS: : yunli@med.unc.edu or ming.hu@nyumc.org

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Mapping of Variable DNA Methylation Across Multiple Cell Types Defines a Dynamic Regulatory Landscape of the Human Genome

DNA methylation is an important epigenetic modification involved in many biological processes and diseases. Many studies have mapped DNA methylation changes associated with embryogenesis, cell differentiation, and cancer at a genome-wide scale. Our understanding of genome-wide DNA methylation changes in a developmental or disease-related context has been steadily growing. However, the investigation of which CpGs are variably methylated in different normal cell or tissue types is still limited. Here, we present an in-depth analysis of 54 single-CpG-resolution DNA methylomes of normal human cell types by integrating high-throughput sequencing-based methylation data. We found that the ratio of methylated to unmethylated CpGs is relatively constant regardless of cell type. However, which CpGs made up the unmethylated complement was cell-type specific. We categorized the 26,000,000 human autosomal CpGs based on their methylation levels across multiple cell types to identify variably methylated CpGs and found that 22.6% exhibited variable DNA methylation. These variably methylated CpGs formed 660,000 variably methylated regions (VMRs), encompassing 11% of the genome. By integrating a multitude of genomic data, we found that VMRs enrich for histone modifications indicative of enhancers, suggesting their role as regulatory elements marking cell type specificity. VMRs enriched for transcription factor binding sites in a tissue-dependent manner. Importantly, they enriched for GWAS variants, suggesting that VMRs could potentially be implicated in disease and complex traits. Taken together, our results highlight the link between CpG methylation variation, genetic variation, and disease risk for many human cell types.

4D nucleomes in single cells: what can computational modeling reveal about spatial chromatin conformation?

Genome-wide sequencing technologies enable investigations of the structural properties of the genome in various spatial dimensions. Here, we review computational techniques developed to model the three-dimensional genome in single cells versus ensembles of cells and assess their underlying assumptions. We further address approaches to study the spatio-temporal aspects of genome organization from single-cell data.

De novo deciphering three-dimensional chromatin interaction and topological domains by wavelet transformation of epigenetic profiles

Defining chromatin interaction frequencies and topological domains is a great challenge for the annotations of genome structures. Although the chromosome conformation capture (3C) and its derivative methods have been developed for exploring the global interactome, they are limited by high experimental complexity and costs. Here we describe a novel computational method, called CITD, for de novo prediction of the chromatin interaction map by integrating histone modification data. We used the public epigenomic data from human fibroblast IMR90 cell and embryonic stem cell (H1) to develop and test CITD, which can not only successfully reconstruct the chromatin interaction frequencies discovered by the Hi-C technology, but also provide additional novel details of chromosomal organizations. We predicted the chromatin interaction frequencies, topological domains and their states (e.g. active or repressive) for 98 additional cell types from Roadmap Epigenomics and ENCODE projects. A total of 131 protein-coding genes located near 78 preserved boundaries among 100 cell types are found to be significantly enriched in functional categories of the nucleosome organization and chromatin assembly. CITD and its predicted results can be used for complementing the topological domains derived from limited Hi-C data and facilitating the understanding of spatial principles underlying the chromosomal organization.

Three-dimensional disorganization of the cancer genome occurs coincident with long-range genetic and epigenetic alterations

A three-dimensional chromatin state underpins the structural and functional basis of the genome by bringing regulatory elements and genes into close spatial proximity to ensure proper, cell-type–specific gene expression profiles. Here, we performed Hi-C chromosome conformation capture sequencing to investigate how three-dimensional chromatin organization is disrupted in the context of copy-number variation, long-range epigenetic remodeling, and atypical gene expression programs in prostate cancer. We find that cancer cells retain the ability to segment their genomes into megabase-sized topologically associated domains (TADs); however, these domains are generally smaller due to establishment of additional domain boundaries. Interestingly, a large proportion of the new cancer-specific domain boundaries occur at regions that display copy-number variation. Notably, a common deletion on 17p13.1 in prostate cancer spanning the TP53 tumor suppressor locus results in bifurcation of a single TAD into two distinct smaller TADs. Change in domain structure is also accompanied by novel cancer-specific chromatin interactions within the TADs that are enriched at regulatory elements such as enhancers, promoters, and insulators, and associated with alterations in gene expression. We also show that differential chromatin interactions across regulatory regions occur within long-range epigenetically activated or silenced regions of concordant gene activation or repression in prostate cancer. Finally, we present a novel visualization tool that enables integrated exploration of Hi-C interaction data, the transcriptome, and epigenome. This study provides new insights into the relationship between long-range epigenetic and genomic dysregulation and changes in higher-order chromatin interactions in cancer.

The folding landscape of the epigenome

The role of the spatial organization of chromatin in gene regulation is a long-standing but still open question. Experimentally it has been shown that the genome is segmented into epigenomic chromatin domains that are organized into hierarchical sub-nuclear spatial compartments. However, whether this non-random spatial organization only reflects or indeed contributes-and how-to the regulation of genome function remains to be elucidated. To address this question, we recently proposed a quantitative description of the folding properties of the fly genome as a function of its epigenomic landscape using a polymer model with epigenomic-driven attractions. We propose in this article, to characterize more deeply the physical properties of the 3D epigenome folding. Using an efficient lattice version of the original block copolymer model, we study the structural and dynamical properties of chromatin and show that the size of epigenomic domains and asymmetries in sizes and in interaction strengths play a critical role in the chromatin organization. Finally, we discuss the biological implications of our findings. In particular, our predictions are quantitatively compatible with experimental data and suggest a different mean of self-interaction in euchromatin versus heterochromatin domains.

Enhancer-promoter interactions are encoded by complex genomic signatures on looping chromatin

Discriminating the gene target of a distal regulatory element from other nearby transcribed genes is a challenging problem with the potential to illuminate the causal underpinnings of complex diseases. We present TargetFinder, a computational method that reconstructs regulatory landscapes from diverse features along the genome. The resulting models accurately predict individual enhancer-promoter interactions across multiple cell lines with a false discovery rate up to 15 times smaller than that obtained using the closest gene. By evaluating the genomic features driving this accuracy, we uncover interactions between structural proteins, transcription factors, epigenetic modifications, and transcription that together distinguish interacting from non-interacting enhancer-promoter pairs. Most of this signature is not proximal to the enhancers and promoters but instead decorates the looping DNA. We conclude that complex but consistent combinations of marks on the one-dimensional genome encode the three-dimensional structure of fine-scale regulatory interactions.

The mechanisms shaping the single-cell transcriptional landscape

Recent technological and computational advances in understanding the transcriptional and chromatin features of single cells have begun answering longstanding questions in the extent and impact of biological heterogeneity. Here, we outline the intrinsic and extrinsic mechanisms that underlie the transcriptional and functional diversity within superficially homogeneous populations, and we discuss how fascinating new studies have afforded novel insight into each mechanism. The studies are chosen in part to include initial reports of novel functional genomics tools where the eventual applications will clearly have profound impact on our understanding the dynamics of cell-to-cell transcriptional variation-from individual cells to whole organisms.

CRISPR Double Cutting through the Labyrinthine Architecture of 3D Genomes

The genomes are organized into ordered and hierarchical topological structures in interphase nuclei. Within discrete territories of each chromosome, topologically associated domains (TADs) play important roles in various nuclear processes such as gene regulation. Inside TADs separated by relatively constitutive boundaries, distal elements regulate their gene targets through specific chromatin-looping contacts such as long-distance enhancer-promoter interactions. High-throughput sequencing studies have revealed millions of potential regulatory DNA elements, which are much more abundant than the mere ∼20,000 genes they control. The recently emerged CRISPR-Cas9 genome editing technologies have enabled efficient and precise genetic and epigenetic manipulations of genomes. The multiplexed and high-throughput CRISPR capabilities facilitate the discovery and dissection of gene regulatory elements. Here, we describe the applications of CRISPR for genome, epigenome, and 3D genome editing, focusing on CRISPR DNA-fragment editing with Cas9 and a pair of sgRNAs to investigate topological folding of chromatin TADs and developmental gene regulation.