Ubiquitous heterogeneity and asymmetry of the chromatin environment at regulatory elements

Intrinsic DNA binding properties demonstrated for lineage-specifying basic helix-loop-helix transcription factors

During development, transcription factors select distinct gene programs, providing the necessary regulatory complexity for temporal and tissue-specific gene expression. How related factors retain specificity, especially when they recognize the same DNA motifs, is not understood. We address this paradox using basic helix-loop-helix (bHLH) transcription factors ASCL1, ASCL2, and MYOD1, crucial mediators of lineage specification. In vivo, these factors recognize the same DNA motifs, yet bind largely different genomic sites and regulate distinct transcriptional programs. This suggests that their ability to identify regulatory targets is defined either by the cellular environment of the partially defined lineages in which they are endogenously expressed, or by intrinsic properties of the factors themselves. To distinguish between these mechanisms, we directly compared the chromatin binding properties of this subset of bHLH factors when ectopically expressed in embryonic stem cells, presenting them with a common chromatin landscape and cellular components. We find that these factors retain distinct binding sites; thus, specificity of binding is an intrinsic property not requiring a restricted landscape or lineage-specific cofactors. Although the ASCL factors and MYOD1 have some distinct DNA motif preference, it is not sufficient to explain the extent of the differential binding. All three factors can bind inaccessible chromatin and induce changes in chromatin accessibility and H3K27ac. A reiterated pattern of DNA binding motifs is uniquely enriched in inaccessible chromatin at sites bound by these bHLH factors. These combined properties define a subclass of lineage-specific bHLH factors and provide context for their central roles in development and disease.

High-Resolution Genome-Wide Mapping of Nucleosome Positioning and Occupancy Level Using Paired-End Sequencing Technology

Because of its profound influence on DNA accessibility for protein binding and thus on the regulation of diverse biological processes, nucleosome positioning has been studied for many years. In the past decade, high-throughput sequencing technologies have opened new perspectives in this research field by allowing the study of nucleosome positioning and occupancy on a genome-wide scale, therefore providing understanding on important aspects of chromatin packaging, as well as on various chromatin-template processes like transcription. In this chapter, we provide the protocol of MNase sequencing for the genome-wide mapping of nucleosomes using MNase to generate mononucleosomal DNA fragments and next-generation sequencing technology to identify their individual location.

Characteristic arrangement of nucleosomes is predictive of chromatin interactions at kilobase resolution

High-throughput chromosome conformation capture (3C) technologies, such as Hi-C, have made it possible to survey 3D genome structure. However, obtaining 3D profiles at kilobase resolution at low cost remains a major challenge. Therefore, we herein present an algorithm for precise identification of chromatin interaction sites at kilobase resolution from MNase-seq data, termed chromatin interaction site detector (CISD), and a CISD-based chromatin loop predictor (CISD_loop) that predicts chromatin–chromatin interactions (CCIs) from low-resolution Hi-C data. We show that the predictions of CISD and CISD_loop overlap closely with chromatin interaction analysis by paired-end tag sequencing (ChIA-PET) anchors and loops, respectively. The validity of CISD/CISD_loop was further supported by a 3C assay at about 5 kb resolution. Finally, we demonstrate that only modest amounts of MNase-seq and Hi-C data are sufficient to achieve ultrahigh resolution CCI maps. Our results suggest that CCIs may result in characteristic nucleosomes arrangement patterns flanking the interaction sites, and our algorithms may facilitate precise and systematic investigations of CCIs on a larger scale than hitherto have been possible.

Global mapping of transcription factor motifs in human aging

Biological aging is a complex process dependent on the interplay of cell autonomous and tissue contextual changes which occur in response to cumulative molecular stress and manifest through adaptive transcriptional reprogramming. Here we describe a transcription factor (TF) meta-analysis of gene expression datasets accrued from 18 tissue sites collected at different biological ages and from 7 different in-vitro aging models. In-vitro aging platforms included replicative senescence and an energy restriction model in quiescence (ERiQ), in which ATP was transiently reduced. TF motifs in promoter regions of trimmed sets of target genes were scanned using JASPAR and TRANSFAC. TF signatures established a global mapping of agglomerating motifs with distinct clusters when ranked hierarchically. Remarkably, the ERiQ profile was shared with the majority of in-vivo aged tissues. Fitting motifs in a minimalistic protein-protein network allowed to probe for connectivity to distinct stress sensors. The DNA damage sensors ATM and ATR linked to the subnetwork associated with senescence. By contrast, the energy sensors PTEN and AMPK connected to the nodes in the ERiQ subnetwork. These data suggest that metabolic dysfunction may be linked to transcriptional patterns characteristic of many aged tissues and distinct from cumulative DNA damage associated with senescence.

Comprehensive nucleosome mapping of the human genome in cancer progression

Altered chromatin structure is a hallmark of cancer, and inappropriate regulation of chromatin structure may represent the origin of transformation. Important studies have mapped human nucleosome distributions genome wide, but the role of chromatin structure in cancer progression has not been addressed. We developed a MNase-Transcription Start Site Sequence Capture method (mTSS-seq) to map the nucleosome distribution at human transcription start sites genome-wide in primary human lung and colon adenocarcinoma tissue. Here, we confirm that nucleosome redistribution is an early, widespread event in lung (LAC) and colon (CRC) adenocarcinoma. These altered nucleosome architectures are consistent between LAC and CRC patient samples indicating that they may serve as important early adenocarcinoma markers. We demonstrate that the nucleosome alterations are driven by the underlying DNA sequence and potentiate transcription factor binding. We conclude that DNA-directed nucleosome redistributions are widespread early in cancer progression. We have proposed an entirely new hierarchical model for chromatin-mediated genome regulation.

Genome-wide profiling of nardilysin target genes reveals its role in epigenetic regulation and cell cycle progression

Post-translational histone modifications, such as acetylation and methylation, are prerequisites for transcriptional regulation. The metalloendopeptidase nardilysin (Nrdc) is a H3K4me2-binding protein that controls thermoregulation and β-cell functions through its transcriptional coregulator function. We herein combined high-throughput ChIP-seq and RNA-seq to achieve the first genome-wide identification of Nrdc target genes. A ChIP-seq analysis of immortalized mouse embryo fibroblasts (iMEF) identified 4053 Nrdc-binding sites, most of which were located in proximal promoter sites (2587 Nrdc-binding genes). Global H3K4me2 levels at Nrdc-binding promoters slightly increased, while H3K9ac levels decreased in the absence of Nrdc. Among Nrdc-binding genes, a comparative RNA-seq analysis identified 448 candidates for Nrdc target genes, among which cell cycle-related genes were significantly enriched. We confirmed decreased mRNA and H3K9ac levels at the promoters of individual genes in Nrdc-deficient iMEF, which were restored by the ectopic introduction of Nrdc. Reduced mRNA levels, but not H3K9ac levels were fully restored by the reintroduction of the peptidase-dead mutant of Nrdc. Furthermore, Nrdc promoted cell cycle progression at multiple stages, which enhanced cell proliferation in vivo. Collectively, our integrative studies emphasize the importance of Nrdc for maintaining a proper epigenetic status and cell growth.

p53 binding sites in normal and cancer cells are characterized by distinct chromatin context

The tumor suppressor protein p53 interacts with DNA in a sequence-dependent manner. Thousands of p53 binding sites have been mapped genome-wide in normal and cancer cells. However, the way p53 selectively binds its cognate sites in different types of cells is not fully understood. Here, we performed a comprehensive analysis of 25 published p53 cistromes and identified 3,551 and 6,039 ‘high-confidence’ binding sites in normal and cancer cells, respectively. Our analysis revealed 2 distinct epigenetic features underlying p53-DNA interactions in vivo. First, p53 binding sites are associated with transcriptionally active histone marks (H3K4me3 and H3K36me3) in normal-cell chromatin, but with repressive histone marks (H3K27me3) in cancer-cell chromatin. Second, p53 binding sites in cancer cells are characterized by a lower level of DNA methylation than their counterparts in normal cells, probably related to global hypomethylation in cancers. Intriguingly, regardless of the cell type, p53 sites are highly enriched in the endogenous retroviral elements of the ERV1 family, highlighting the importance of this repeat family in shaping the transcriptional network of p53. Moreover, the p53 sites exhibit an unusual combination of chromatin patterns: high nucleosome occupancy and, at the same time, high sensitivity to DNase I. Our results suggest that p53 can access its target sites in a chromatin environment that is non-permissive to most DNA-binding transcription factors, which may allow p53 to act as a pioneer transcription factor in the context of chromatin.

Bisulfite-independent analysis of CpG island methylation enables genome-scale stratification of single cells

Conventional DNA bisulfite sequencing has been extended to single cell level, but the coverage consistency is insufficient for parallel comparison. Here we report a novel method for genome-wide CpG island (CGI) methylation sequencing for single cells (scCGI-seq), combining methylation-sensitive restriction enzyme digestion and multiple displacement amplification for selective detection of methylated CGIs. We applied this method to analyzing single cells from two types of hematopoietic cells, K562 and GM12878 and small populations of fibroblasts and induced pluripotent stem cells. The method detected 21 798 CGIs (76% of all CGIs) per cell, and the number of CGIs consistently detected from all 16 profiled single cells was 20 864 (72.7%), with 12 961 promoters covered. This coverage represents a substantial improvement over results obtained using single cell reduced representation bisulfite sequencing, with a 66-fold increase in the fraction of consistently profiled CGIs across individual cells. Single cells of the same type were more similar to each other than to other types, but also displayed epigenetic heterogeneity. The method was further validated by comparing the CpG methylation pattern, methylation profile of CGIs/promoters and repeat regions and 41 classes of known regulatory markers to the ENCODE data. Although not every minor methylation differences between cells are detectable, scCGI-seq provides a solid tool for unsupervised stratification of a heterogeneous cell population.

Genome-wide search for Zelda-like chromatin signatures identifies GAF as a pioneer factor in early fly development

BACKGROUND

The protein Zelda was shown to play a key role in early Drosophila development, binding thousands of promoters and enhancers prior to maternal-to-zygotic transition (MZT), and marking them for transcriptional activation. Recently, we showed that Zelda acts through specific chromatin patterns of histone modifications to mark developmental enhancers and active promoters. Intriguingly, some Zelda sites still maintain these chromatin patterns in Drosophila embryos lacking maternal Zelda protein. This suggests that additional Zelda-like pioneer factors may act in early fly embryos.

RESULTS

We developed a computational method to analyze and refine the chromatin landscape surrounding early Zelda peaks, using a multichannel spectral clustering. This allowed us to characterize their chromatin patterns through MZT (mitotic cycles 8-14). Specifically, we focused on H3K4me1, H3K4me3, H3K18ac, H3K27ac, and H3K27me3 and identified three different classes of chromatin signatures, matching "promoters," "enhancers" and "transiently bound" Zelda peaks. We then further scanned the genome using these chromatin patterns and identified additional loci-with no Zelda binding-that show similar chromatin patterns, resulting with hundreds of Zelda-independent putative enhancers. These regions were found to be enriched with GAGA factor (GAF, Trl) and are typically located near early developmental zygotic genes. Overall our analysis suggests that GAF, together with Zelda, plays an important role in activating the zygotic genome.

CONCLUSIONS

As we show, our computational approach offers an efficient algorithm for characterizing chromatin signatures around some loci of interest and allows a genome-wide identification of additional loci with similar chromatin patterns.

NucTools: analysis of chromatin feature occupancy profiles from high-throughput sequencing data

Background

Biomedical applications of high-throughput sequencing methods generate a vast amount of data in which numerous chromatin features are mapped along the genome. The results are frequently analysed by creating binary data sets that link the presence/absence of a given feature to specific genomic loci. However, the nucleosome occupancy or chromatin accessibility landscape is essentially continuous. It is currently a challenge in the field to cope with continuous distributions of deep sequencing chromatin readouts and to integrate the different types of discrete chromatin features to reveal linkages between them.

Results

Here we introduce the NucTools suite of Perl scripts as well as MATLAB- and R-based visualization programs for a nucleosome-centred downstream analysis of deep sequencing data. NucTools accounts for the continuous distribution of nucleosome occupancy. It allows calculations of nucleosome occupancy profiles averaged over several replicates, comparisons of nucleosome occupancy landscapes between different experimental conditions, and the estimation of the changes of integral chromatin properties such as the nucleosome repeat length. Furthermore, NucTools facilitates the annotation of nucleosome occupancy with other chromatin features like binding of transcription factors or architectural proteins, and epigenetic marks like histone modifications or DNA methylation. The applications of NucTools are demonstrated for the comparison of several datasets for nucleosome occupancy in mouse embryonic stem cells (ESCs) and mouse embryonic fibroblasts (MEFs).

Conclusions

The typical workflows of data processing and integrative analysis with NucTools reveal information on the interplay of nucleosome positioning with other features such as for example binding of a transcription factor CTCF, regions with stable and unstable nucleosomes, and domains of large organized chromatin K9me2 modifications (LOCKs). As potential limitations and problems we discuss how inter-replicate variability of MNase-seq experiments can be addressed.

Electronic supplementary material

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

GC Content Heterogeneity Transition of Conserved Noncoding Sequences Occurred at the Emergence of Vertebrates

Conserved non-coding sequences (CNSs) of Eukaryotes are known to be significantly enriched in regulatory sequences. CNSs of diverse lineages follow different patterns in abundance, sequence composition, and location. Here, we report a thorough analysis of CNSs in diverse groups of Eukaryotes with respect to GC content heterogeneity. We examined 24 fungi, 19 invertebrates, and 12 non-mammalian vertebrates so as to find lineage specific features of CNSs. We found that fungi and invertebrate CNSs are predominantly GC rich as in plants we previously observed, whereas vertebrate CNSs are GC poor. This result suggests that the CNS GC content transition occurred from the ancestral GC rich state of Eukaryotes to GC poor in the vertebrate lineage due to the enrollment of GC poor transcription factor binding sites that are lineage specific. CNS GC content is closely linked with the nucleosome occupancy that determines the location and structural architecture of DNAs.

DGW: an exploratory data analysis tool for clustering and visualisation of epigenomic marks

Background

Functional genomic and epigenomic research relies fundamentally on sequencing based methods like ChIP-seq for the detection of DNA-protein interactions. These techniques return large, high dimensional data sets with visually complex structures, such as multi-modal peaks extended over large genomic regions. Current tools for visualisation and data exploration represent and leverage these complex features only to a limited extent.

Results

We present DGW, an open source software package for simultaneous alignment and clustering of multiple epigenomic marks. DGW uses Dynamic Time Warping to adaptively rescale and align genomic distances which allows to group regions of interest with similar shapes, thereby capturing the structure of epigenomic marks. We demonstrate the effectiveness of the approach in a simulation study and on a real epigenomic data set from the ENCODE project.

Conclusions

Our results show that DGW automatically recognises and aligns important genomic features such as transcription start sites and splicing sites from histone marks. DGW is available as an open source Python package.

Genome-wide uniformity of human 'open' pre-initiation complexes

Transcription of protein-coding and noncoding DNA occurs pervasively throughout the mammalian genome. Their sites of initiation are generally inferred from transcript 5' ends and are thought to be either locally dispersed or focused. How these two modes of initiation relate is unclear. Here, we apply permanganate treatment and chromatin immunoprecipitation (PIP-seq) of initiation factors to identify the precise location of melted DNA separately associated with the preinitiation complex (PIC) and the adjacent paused complex (PC). This approach revealed the two known modes of transcription initiation. However, in contrast to prevailing views, they co-occurred within the same promoter region: initiation originating from a focused PIC, and broad nucleosome-linked initiation. PIP-seq allowed transcriptional orientation of Pol II to be determined, which may be useful near promoters where sufficient sense/anti-sense transcript mapping information is lacking. PIP-seq detected divergently oriented Pol II at both coding and noncoding promoters, as well as at enhancers. Their occupancy levels were not necessarily coupled in the two orientations. DNA sequence and shape analysis of initiation complex sites suggest that both sequence and shape contribute to specificity, but in a context-restricted manner. That is, initiation sites have the locally "best" initiator (INR) sequence and/or shape. These findings reveal a common core to pervasive Pol II initiation throughout the human genome.

metagene Profiles Analyses Reveal Regulatory Element's Factor-Specific Recruitment Patterns

ChIP-Sequencing (ChIP-Seq) provides a vast amount of information regarding the localization of proteins across the genome. The aggregation of ChIP-Seq enrichment signal in a metagene plot is an approach commonly used to summarize data complexity and to obtain a high level visual representation of the general occupancy pattern of a protein. Here we present the R package metagene, the graphical interface Imetagene and the companion package similaRpeak. Together, they provide a framework to integrate, summarize and compare the ChIP-Seq enrichment signal from complex experimental designs. Those packages identify and quantify similarities or dissimilarities in patterns between large numbers of ChIP-Seq profiles. We used metagene to investigate the differential occupancy of regulatory factors at noncoding regulatory regions (promoters and enhancers) in relation to transcriptional activity in GM12878 B-lymphocytes. The relationships between occupancy patterns and transcriptional activity suggest two different mechanisms of action for transcriptional control: i) a “gradient effect” where the regulatory factor occupancy levels follow transcription and ii) a “threshold effect” where the regulatory factor occupancy levels max out prior to reaching maximal transcription. metagene, Imetagene and similaRpeak are implemented in R under the Artistic license 2.0 and are available on Bioconductor.

MNase titration reveals differences between nucleosome occupancy and chromatin accessibility

Chromatin accessibility plays a fundamental role in gene regulation. Nucleosome placement, usually measured by quantifying protection of DNA from enzymatic digestion, can regulate accessibility. We introduce a metric that uses micrococcal nuclease (MNase) digestion in a novel manner to measure chromatin accessibility by combining information from several digests of increasing depths. This metric, MACC (MNase accessibility), quantifies the inherent heterogeneity of nucleosome accessibility in which some nucleosomes are seen preferentially at high MNase and some at low MNase. MACC interrogates each genomic locus, measuring both nucleosome location and accessibility in the same assay. MACC can be performed either with or without a histone immunoprecipitation step, and thereby compares histone and non-histone protection. We find that changes in accessibility at enhancers, promoters and other regulatory regions do not correlate with changes in nucleosome occupancy. Moreover, high nucleosome occupancy does not necessarily preclude high accessibility, which reveals novel principles of chromatin regulation.

The topography of mutational processes in breast cancer genomes

Somatic mutations in human cancers show unevenness in genomic distribution that correlate with aspects of genome structure and function. These mutations are, however, generated by multiple mutational processes operating through the cellular lineage between the fertilized egg and the cancer cell, each composed of specific DNA damage and repair components and leaving its own characteristic mutational signature on the genome. Using somatic mutation catalogues from 560 breast cancer whole-genome sequences, here we show that each of 12 base substitution, 2 insertion/deletion (indel) and 6 rearrangement mutational signatures present in breast tissue, exhibit distinct relationships with genomic features relating to transcription, DNA replication and chromatin organization. This signature-based approach permits visualization of the genomic distribution of mutational processes associated with APOBEC enzymes, mismatch repair deficiency and homologous recombinational repair deficiency, as well as mutational processes of unknown aetiology. Furthermore, it highlights mechanistic insights including a putative replication-dependent mechanism of APOBEC-related mutagenesis.

Exploration of nucleosome positioning patterns in transcription factor function

The binding of transcription factors (TFs) triggers activation of specific chromatin regions through the recruitment and activation of RNA polymerase. Unique nucleosome positioning (NP) occurs during gene expression and has been suggested to be involved in various other chromatin functions. However, the diversity of NP that can occur for each function has not been clarified. Here we used MNase-Seq data to evaluate NP around 258 cis-regulatory elements in the mouse genome. Principal component analysis of the 258 elements revealed that NP consisted of five major patterns. Furthermore, the five NP patterns had predictive power for the level of gene expression. We also demonstrated that selective NP patterns appeared around TF binding sites. These results suggest that the NP patterns are correlated to specific functions on chromatin.

Stress and glucocorticoid receptor transcriptional programming in time and space: Implications for the brain-gut axis

BACKGROUND

Chronic psychological stress is associated with enhanced abdominal pain and altered intestinal barrier function that may result from a perturbation in the hypothalamic-pituitary-adrenal (HPA) axis. The glucocorticoid receptor (GR) exploits diverse mechanisms to activate or suppress congeneric gene expression, with regulatory variation associated with stress-related disorders in psychiatry and gastroenterology.

PURPOSE

During acute and chronic stress, corticotropin-releasing hormone drives secretion of adrenocorticotropic hormone from the pituitary, ultimately leading to the release of cortisol (human) and corticosterone (rodent) from the adrenal glands. Cortisol binds with the GR in the cytosol, translocates to the nucleus, and activates the NR3C1 (nuclear receptor subfamily 3, group C, member 1 [GR]) gene. This review focuses on the rapidly developing observations that cortisol is responsible for driving circadian and ultradian bursts of transcriptional activity in the CLOCK (clock circadian regulator) and PER (period circadian clock 1) gene families, and this rhythm is disrupted in major depressive disorder, bipolar disorder, and stress-related gastrointestinal and immune disorders. Glucocorticoid receptor regulates different sets of transcripts in a tissue-specific manner, through pulsatile waves of gene expression that includes occupancy of glucocorticoid response elements located within constitutively open spatial domains in chromatin. Emerging evidence supports a potentially pivotal role for epigenetic regulation of how GR interacts with other chromatin regulators to control the expression of its target genes. Dysregulation of the central and peripheral GR regulome has potentially significant consequences for stress-related disorders affecting the brain-gut axis.

Dynamic recruitment of Ets1 to both nucleosome-occupied and -depleted enhancer regions mediates a transcriptional program switch during early T-cell differentiation

Ets1 is a sequence-specific transcription factor that plays an important role during hematopoiesis, and is essential for the transition of CD4⁻/CD8⁻ double negative (DN) to CD4⁺/CD8⁺ double positive (DP) thymocytes. Using genome-wide and functional approaches, we investigated the binding properties, transcriptional role and chromatin environment of Ets1 during this transition. We found that while Ets1 binding at distal sites was associated with active genes at both DN and DP stages, its enhancer activity was attained at the DP stage, as reflected by levels of the core transcriptional hallmarks H3K4me1/3, RNA Polymerase II and eRNA. This dual, stage-specific ability reflected a switch from non-T hematopoietic toward T-cell specific gene expression programs during the DN-to-DP transition, as indicated by transcriptome analyses of Ets1^−/− thymic cells. Coincidentally, Ets1 associates more specifically with Runx1 in DN and with TCF1 in DP cells. We also provide evidence that Ets1 predominantly binds distal nucleosome-occupied regions in DN and nucleosome-depleted regions in DP. Finally and importantly, we demonstrate that Ets1 induces chromatin remodeling by displacing H3K4me1-marked nucleosomes. Our results thus provide an original model whereby the ability of a transcription factor to bind nucleosomal DNA changes during differentiation with consequences on its cognate enhancer activity.

Functionally distinct patterns of nucleosome remodeling at enhancers in glucocorticoid-treated acute lymphoblastic leukemia

Background

Precise nucleosome positioning is an increasingly recognized feature of promoters and enhancers, reflecting complex contributions of DNA sequence, nucleosome positioning, histone modification and transcription factor binding to enhancer activity and regulation of gene expression. Changes in nucleosome position and occupancy, histone variants and modifications, and chromatin remodeling are also critical elements of dynamic transcriptional regulation, but poorly understood at enhancers. We investigated glucocorticoid receptor-associated (GR) nucleosome dynamics at enhancers in acute lymphoblastic leukemia.

Results

For the first time, we demonstrate functionally distinct modes of nucleosome remodeling upon chromatin binding by GR, which we term central, non-central, phased, and minimal. Central and non-central remodeling reflect nucleosome eviction by GR and cofactors, respectively. Phased remodeling involves nucleosome repositioning and is associated with rapidly activated enhancers and induction of gene expression. Minimal remodeling sites initially have low levels of enhancer-associated histone modification, but the majority of these regions gain H3K4me2 or H3K27Ac to become de novo enhancers. Minimal remodeling regions are associated with gene ontologies specific to decreased B cell number and mTOR inhibition and may make unique contributions to glucocorticoid-induced leukemia cell death.

Conclusions

Our findings form a novel framework for understanding the dynamic interplay between transcription factor binding, nucleosome remodeling, enhancer function, and gene expression in the leukemia response to glucocorticoids.

Electronic supplementary material

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

Genome-wide predictors of NF-κB recruitment and transcriptional activity

BACKGROUND

Inducible transcription factors (TFs) mediate transcriptional responses to environmental cues. In response to multiple inflammatory signals active NF-κB dimers enter the nucleus and trigger cell-type-, and stimulus-specific transcriptional programs. Although much is known about NF-κB inducing pathways and about locus-specific mechanisms of transcriptional control, it is poorly understood how the pre-existing chromatin landscape determines NF-κB target selection and activation. Specifically, it is not known which epigenetic marks and pre-bound TFs serve genome-wide as positive (negative) cues for active NF-κB.

RESULTS

We applied multivariate and combinatorial data mining techniques on a comprehensive dataset of DNA methylation, DNase I hypersensitivity, eight epigenetic marks, and 34 TFs to arrive at genome-wide patterns that predict NF-κB binding. Strikingly, we observed NF-κB recruitment to accessible and nucleosome-bound sites. Within nucleosomal DNA NF-κB binding was primed by H3K4me1 and H2A.Z, but also hyper-methylated DNA outside of promoters and CpG-islands. Many of these predictors showed combinatorial cooperativity and statistically significant interactions. Recruitment to pre-accessible sites was more frequent and influenced by chromatin-associated TFs. We observed that specific TF-combinations are greatly enriched for (or depleted of) NF-κB binding events.

CONCLUSIONS

We provide evidence of NF-κB binding within genomic regions that lack classical marks of activity. These pioneer binding events are relatively often associated with transcriptional regulation. Further, our predictive models indicate that specific combinations of epigenetic marks and transcription factors predetermine the NF-κB cistrome, supporting the feasibility of using statistical approaches to identify "histone codes".

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.

The histone variant H2A.Z is an important regulator of enhancer activity

Gene regulatory programs in different cell types are largely defined through cell-specific enhancers activity. The histone variant H2A.Z has been shown to play important roles in transcription mainly by controlling proximal promoters, but its effect on enhancer functions remains unclear. Here, we demonstrate by genome-wide approaches that H2A.Z is present at a subset of active enhancers bound by the estrogen receptor alpha (ERα). We also determine that H2A.Z does not influence the local nucleosome positioning around ERα enhancers using ChIP sequencing at nucleosomal resolution and unsupervised pattern discovery. We further highlight that H2A.Z-enriched enhancers are associated with chromatin accessibility, H3K122ac enrichment and hypomethylated DNA. Moreover, upon estrogen stimulation, the enhancers occupied by H2A.Z produce enhancer RNAs (eRNAs), and recruit RNA polymerase II as well as RAD21, a member of the cohesin complex involved in chromatin interactions between enhancers and promoters. Importantly, their recruitment and eRNAs production are abolished by H2A.Z depletion, thereby revealing a novel functional link between H2A.Z occupancy and enhancer activity. Taken together, our findings suggest that H2A.Z acts as an important player for enhancer functions by establishing and maintaining a chromatin environment required for RNA polymerase II recruitment, eRNAs transcription and enhancer-promoters interactions, all essential attributes of enhancer activity.

An Efficient Method for Electroporation of Small Interfering RNAs into ENCODE Project Tier 1 GM12878 and K562 Cell Lines

The Encyclopedia of DNA Elements (ENCODE) Project aims to identify all functional sequence elements in the human genome sequence by use of high-throughput DNA/cDNA sequencing approaches. To aid the standardization, comparison, and integration of data sets produced from different technologies and platforms, the ENCODE Consortium selected several standard human cell lines to be used by the ENCODE Projects. The Tier 1 ENCODE cell lines include GM12878, K562, and H1 human embryonic stem cell lines. GM12878 is a lymphoblastoid cell line, transformed with the Epstein-Barr virus, that was selected by the International HapMap Project for whole genome and transcriptome sequencing by use of the Illumina platform. K562 is an immortalized myelogenous leukemia cell line. The GM12878 cell line is attractive for the ENCODE Projects, as it offers potential synergy with the International HapMap Project. Despite the vast amount of sequencing data available on the GM12878 cell line through the ENCODE Project, including transcriptome, chromatin immunoprecipitation-sequencing for histone marks, and transcription factors, no small interfering siRNA-mediated knockdown studies have been performed in the GM12878 cell line, as cationic lipid-mediated transfection methods are inefficient for lymphoid cell lines. Here, we present an efficient and reproducible method for transfection of a variety of siRNAs into the GM12878 and K562 cell lines, which subsequently results in targeted protein depletion.

Promoter-like epigenetic signatures in exons displaying cell type-specific splicing

BACKGROUND

Pre-mRNA splicing occurs mainly co-transcriptionally, and both nucleosome density and histone modifications have been proposed to play a role in splice site recognition and regulation. However, the extent and mechanisms behind this interplay remain poorly understood.

RESULTS

We use transcriptomic and epigenomic data generated by the ENCODE project to investigate the association between chromatin structure and alternative splicing. We find a strong and significant positive association between H3K9ac, H3K27ac, H3K4me3, epigenetic marks characteristic of active promoters, and exon inclusion in a small but well-defined class of exons, representing approximately 4 % of all regulated exons. These exons are systematically maintained at comparatively low levels of inclusion across cell types, but their inclusion is significantly enhanced in particular cell types when in physical proximity to active promoters.

CONCLUSION

Histone modifications and other chromatin features that activate transcription can be co-opted to participate in the regulation of the splicing of exons that are in physical proximity to promoter regions.

CRISPR-Cas9-mediated genome editing and guide RNA design

CRISPR and CRISPR-associated (Cas) proteins, which in nature comprise the RNA-based adaptive immune system in bacteria and archaea, have emerged as particularly powerful genome editing tools owing to their unrivaled ease of use and ability to modify genomes across mammalian model systems. As such, the CRISPR-Cas9 system holds promise as a "system of choice" for functional mammalian genetic studies across biological disciplines. Here we briefly review this fast moving field, introduce the CRISPR-Cas9 system and its application to genome editing, with a focus on the basic considerations in designing the targeting guide RNA sequence.

Nucleosome positioning: resources and tools online

Nucleosome positioning is an important process required for proper genome packing and its accessibility to execute the genetic program in a cell-specific, timely manner. In the recent years hundreds of papers have been devoted to the bioinformatics, physics and biology of nucleosome positioning. The purpose of this review is to cover a practical aspect of this field, namely, to provide a guide to the multitude of nucleosome positioning resources available online. These include almost 300 experimental datasets of genome-wide nucleosome occupancy profiles determined in different cell types and more than 40 computational tools for the analysis of experimental nucleosome positioning data and prediction of intrinsic nucleosome formation probabilities from the DNA sequence. A manually curated, up to date list of these resources will be maintained at http://generegulation.info.

Structured nucleosome fingerprints enable high-resolution mapping of chromatin architecture within regulatory regions

Transcription factors canonically bind nucleosome-free DNA, making the positioning of nucleosomes within regulatory regions crucial to the regulation of gene expression. Using the assay of transposase accessible chromatin (ATAC-seq), we observe a highly structured pattern of DNA fragment lengths and positions around nucleosomes in Saccharomyces cerevisiae, and use this distinctive two-dimensional nucleosomal “fingerprint” as the basis for a new nucleosome-positioning algorithm called NucleoATAC. We show that NucleoATAC can identify the rotational and translational positions of nucleosomes with up to base-pair resolution and provide quantitative measures of nucleosome occupancy in S. cerevisiae, Schizosaccharomyces pombe, and human cells. We demonstrate the application of NucleoATAC to a number of outstanding problems in chromatin biology, including analysis of sequence features underlying nucleosome positioning, promoter chromatin architecture across species, identification of transient changes in nucleosome occupancy and positioning during a dynamic cellular response, and integrated analysis of nucleosome occupancy and transcription factor binding.

Genetic Control of Chromatin States in Humans Involves Local and Distal Chromosomal Interactions

Deciphering the impact of genetic variants on gene regulation is fundamental to understanding human disease. Although gene regulation often involves long-range interactions, it is unknown to what extent non-coding genetic variants influence distal molecular phenotypes. Here, we integrate chromatin profiling for three histone marks in lymphoblastoid cell lines (LCLs) from 75 sequenced individuals with LCL-specific Hi-C and ChIA-PET-based chromatin contact maps to uncover one of the largest collections of local and distal histone quantitative trait loci (hQTLs). Distal QTLs are enriched within topologically associated domains and exhibit largely concordant variation of chromatin state coordinated by proximal and distal non-coding genetic variants. Histone QTLs are enriched for common variants associated with autoimmune diseases and enable identification of putative target genes of disease-associated variants from genome-wide association studies. These analyses provide insights into how genetic variation can affect human disease phenotypes by coordinated changes in chromatin at interacting regulatory elements.

Nucleosome Organization in Human Embryonic Stem Cells

The fundamental repeating unit of eukaryotic chromatin is the nucleosome. Besides being involved in packaging DNA, nucleosome organization plays an important role in transcriptional regulation and cellular identity. Currently, there is much debate about the major determinants of the nucleosome architecture of a genome and its significance with little being known about its role in stem cells. To address these questions, we performed ultra-deep sequencing of nucleosomal DNA in two human embryonic stem cell lines and integrated our data with numerous epigenomic maps. Our analyses have revealed that the genome is a determinant of nucleosome organization with transcriptionally inactive regions characterized by a “ground state” of nucleosome profiles driven by underlying DNA sequences. DNA sequence preferences are associated with heterogeneous chromatin organization around transcription start sites. Transcription, histone modifications, and DNA methylation alter this “ground state” by having distinct effects on both nucleosome positioning and occupancy. As the transcriptional rate increases, nucleosomes become better positioned. Exons transcribed and included in the final spliced mRNA have distinct nucleosome profiles in comparison to exons not included at exon-exon junctions. Genes marked by the active modification H3K4m3 are characterized by lower nucleosome occupancy before the transcription start site compared to genes marked by the inactive modification H3K27m3, while bivalent domains, genes associated with both marks, lie exactly in the middle. Combinatorial patterns of epigenetic marks (chromatin states) are associated with unique nucleosome profiles. Nucleosome organization varies around transcription factor binding in enhancers versus promoters. DNA methylation is associated with increasing nucleosome occupancy and different types of methylations have distinct location preferences within the nucleosome core particle. Finally, computational analysis of nucleosome organization alone is sufficient to elucidate much of the circuitry of pluripotency. Our results, suggest that nucleosome organization is associated with numerous genomic and epigenomic processes and can be used to elucidate cellular identity.

Efficient inversions and duplications of mammalian regulatory DNA elements and gene clusters by CRISPR/Cas9

The human genome contains millions of DNA regulatory elements and a large number of gene clusters, most of which have not been tested experimentally. The clustered regularly interspaced short palindromic repeats (CRISPR)/CRISPR-associated nuclease 9 (Cas9) programed with a synthetic single-guide RNA (sgRNA) emerges as a method for genome editing in virtually any organisms. Here we report that targeted DNA fragment inversions and duplications could easily be achieved in human and mouse genomes by CRISPR with two sgRNAs. Specifically, we found that, in cultured human cells and mice, efficient precise inversions of DNA fragments ranging in size from a few tens of bp to hundreds of kb could be generated. In addition, DNA fragment duplications and deletions could also be generated by CRISPR through trans-allelic recombination between the Cas9-induced double-strand breaks (DSBs) on two homologous chromosomes (chromatids). Moreover, junctions of combinatorial inversions and duplications of the protocadherin (Pcdh) gene clusters induced by Cas9 with four sgRNAs could be detected. In mice, we obtained founders with alleles of precise inversions, duplications, and deletions of DNA fragments of variable sizes by CRISPR. Interestingly, we found that very efficient inversions were mediated by microhomology-mediated end joining (MMEJ) through short inverted repeats. We showed for the first time that DNA fragment inversions could be transmitted through germlines in mice. Finally, we applied this CRISPR method to a regulatory element of the Pcdhα cluster and found a new role in the regulation of members of the Pcdhγ cluster. This simple and efficient method should be useful in manipulating mammalian genomes to study millions of regulatory DNA elements as well as vast numbers of gene clusters.

The oxidative demethylase ALKBH3 marks hyperactive gene promoters in human cancer cells

Background

The oxidative DNA demethylase ALKBH3 targets single-stranded DNA (ssDNA) in order to perform DNA alkylation damage repair. ALKBH3 becomes upregulated during tumorigenesis and is necessary for proliferation. However, the underlying molecular mechanism remains to be understood.

Methods

To further elucidate the function of ALKBH3 in cancer, we performed ChIP-seq to investigate the genomic binding pattern of endogenous ALKBH3 in PC3 prostate cancer cells coupled with microarray experiments to examine the expression effects of ALKBH3 depletion.

Results

We demonstrate that ALKBH3 binds to transcription associated locations, such as places of promoter-proximal paused RNA polymerase II and enhancers. Strikingly, ALKBH3 strongly binds to the transcription initiation sites of a small number of highly active gene promoters. These promoters are characterized by high levels of transcriptional regulators, including transcription factors, the Mediator complex, cohesin, histone modifiers, and active histone marks. Gene expression analysis showed that ALKBH3 does not directly influence the transcription of its target genes, but its depletion induces an upregulation of ALKBH3 non-bound inflammatory genes.

Conclusions

The genomic binding pattern of ALKBH3 revealed a putative novel hyperactive promoter type. Further, we propose that ALKBH3 is an intrinsic DNA repair protein that suppresses transcription associated DNA damage at highly expressed genes and thereby plays a role to maintain genomic integrity in ALKBH3-overexpressing cancer cells. These results raise the possibility that ALKBH3 may be a potential target for inhibiting cancer progression.

Electronic supplementary material

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

Emerging applications of read profiles towards the functional annotation of the genome

Functional annotation of the genome is important to understand the phenotypic complexity of various species. The road toward functional annotation involves several challenges ranging from experiments on individual molecules to large-scale analysis of high-throughput sequencing (HTS) data. HTS data is typically a result of the protocol designed to address specific research questions. The sequencing results in reads, which when mapped to a reference genome often leads to the formation of distinct patterns (read profiles). Interpretation of these read profiles is essential for their analysis in relation to the research question addressed. Several strategies have been employed at varying levels of abstraction ranging from a somewhat ad hoc to a more systematic analysis of read profiles. These include methods which can compare read profiles, e.g., from direct (non-sequence based) alignments to classification of patterns into functional groups. In this review, we highlight the emerging applications of read profiles for the annotation of non-coding RNA and cis-regulatory elements (CREs) such as enhancers and promoters. We also discuss the biological rationale behind their formation.

Chemical profiling of the genome with anti-cancer drugs defines target specificities

Many anticancer drugs induce DNA breaks to eliminate tumor cells. The anthracycline topoisomerase II inhibitors additionally cause histone eviction. Here, we performed genome-wide high-resolution mapping of chemotherapeutic effects of various topoisomerase I and II (TopoI and II) inhibitors and integrated this mapping with established maps of genomic or epigenomic features to show their activities in different genomic regions. The TopoI inhibitor topotecan and the TopoII inhibitor etoposide are similar in inducing DNA damage at transcriptionally active genomic regions. The anthracycline daunorubicin induces DNA breaks and evicts histones from active chromatin, thus quenching local DNA damage responses. Another anthracycline, aclarubicin, has a different genomic specificity and evicts histones from H3K27me3-marked heterochromatin, with consequences for diffuse large B-cell lymphoma cells with elevated levels of H3K27me3. Modifying anthracycline structures may yield compounds with selectivity for different genomic regions and activity for different tumor types.

Occupancy by key transcription factors is a more accurate predictor of enhancer activity than histone modifications or chromatin accessibility

Background

Regulated gene expression controls organismal development, and variation in regulatory patterns has been implicated in complex traits. Thus accurate prediction of enhancers is important for further understanding of these processes. Genome-wide measurement of epigenetic features, such as histone modifications and occupancy by transcription factors, is improving enhancer predictions, but the contribution of these features to prediction accuracy is not known. Given the importance of the hematopoietic transcription factor TAL1 for erythroid gene activation, we predicted candidate enhancers based on genomic occupancy by TAL1 and measured their activity. Contributions of multiple features to enhancer prediction were evaluated based on the results of these and other studies.

Results

TAL1-bound DNA segments were active enhancers at a high rate both in transient transfections of cultured cells (39 of 79, or 56%) and transgenic mice (43 of 66, or 65%). The level of binding signal for TAL1 or GATA1 did not help distinguish TAL1-bound DNA segments as active versus inactive enhancers, nor did the density of regulation-related histone modifications. A meta-analysis of results from this and other studies (273 tested predicted enhancers) showed that the presence of TAL1, GATA1, EP300, SMAD1, H3K4 methylation, H3K27ac, and CAGE tags at DNase hypersensitive sites gave the most accurate predictors of enhancer activity, with a success rate over 80% and a median threefold increase in activity. Chromatin accessibility assays and the histone modifications H3K4me1 and H3K27ac were sensitive for finding enhancers, but they have high false positive rates unless transcription factor occupancy is also included.

Conclusions

Occupancy by key transcription factors such as TAL1, GATA1, SMAD1, and EP300, along with evidence of transcription, improves the accuracy of enhancer predictions based on epigenetic features.

Electronic supplementary material

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

ClusTrack: feature extraction and similarity measures for clustering of genome-wide data sets

Clustering is a popular technique for explorative analysis of data, as it can reveal subgroupings and similarities between data in an unsupervised manner. While clustering is routinely applied to gene expression data, there is a lack of appropriate general methodology for clustering of sequence-level genomic and epigenomic data, e.g. ChIP-based data. We here introduce a general methodology for clustering data sets of coordinates relative to a genome assembly, i.e. genomic tracks. By defining appropriate feature extraction approaches and similarity measures, we allow biologically meaningful clustering to be performed for genomic tracks using standard clustering algorithms. An implementation of the methodology is provided through a tool, ClusTrack, which allows fine-tuned clustering analyses to be specified through a web-based interface. We apply our methods to the clustering of occupancy of the H3K4me1 histone modification in samples from a range of different cell types. The majority of samples form meaningful subclusters, confirming that the definitions of features and similarity capture biological, rather than technical, variation between the genomic tracks. Input data and results are available, and can be reproduced, through a Galaxy Pages document at http://hyperbrowser.uio.no/hb/u/hb-superuser/p/clustrack. The clustering functionality is available as a Galaxy tool, under the menu option "Specialized analyzis of tracks", and the submenu option "Cluster tracks based on genome level similarity", at the Genomic HyperBrowser server: http://hyperbrowser.uio.no/hb/.

Proteogenomics analysis reveals specific genomic orientations of distal regulatory regions composed by non-canonical histone variants

Background

Histone variants play further important roles in DNA packaging and controlling gene expression. However, our understanding about their composition and their functions is limited.

Results

Integrating proteomic and genomic approaches, we performed a comprehensive analysis of the epigenetic landscapes containing the four histone variants H3.1, H3.3, H2A.Z, and macroH2A. These histones were FLAG-tagged in HeLa cells and purified using chromatin immunoprecipitation (ChIP). By adopting ChIP followed by mass spectrometry (ChIP-MS), we quantified histone post-translational modifications (PTMs) and histone variant nucleosomal ratios in highly purified mononucleosomes. Subsequent ChIP followed by next-generation sequencing (ChIP-seq) was used to map the genome-wide localization of the analyzed histone variants and define their chromatin domains. Finally, we included in our study large datasets contained in the ENCODE database. We newly identified a group of regulatory regions enriched in H3.1 and the histone variant associated with repressive marks macroH2A. Systematic analysis identified both symmetric and asymmetric patterns of histone variant occupancies at intergenic regulatory regions. Strikingly, these directional patterns were associated with RNA polymerase II (PolII). These asymmetric patterns correlated with the enhancer activities measured using global run-on sequencing (GRO-seq) data.

Conclusions

Our studies show that H2A.Z and H3.3 delineate the orientation of transcription at enhancers as observed at promoters. We also showed that enhancers with skewed histone variant patterns well facilitate enhancer activity. Collectively, our study indicates that histone variants are deposited at regulatory regions to assist gene regulation.

Electronic supplementary material

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

Integrative epigenomic analysis of differential DNA methylation in urothelial carcinoma

BACKGROUND

Urothelial carcinoma of the bladder (UC) is a common malignancy. Although extensive transcriptome analysis has provided insights into the gene expression patterns of this tumor type, the mechanistic underpinnings of differential methylation remain poorly understood. Multi-level genomic data may be used to profile the regulatory potential and landscape of differential methylation in cancer and gain understanding of the processes underlying epigenetic and phenotypic characteristics of tumors.

METHODS

We perform genome-wide DNA methylation profiling of 98 gene-expression subtyped tumors to identify between-tumor differentially methylated regions (DMRs). We integrate multi-level publically available genomic data generated by the ENCODE consortium to characterize the regulatory potential of UC DMRs.

RESULTS

We identify 5,453 between-tumor DMRs and derive four DNA methylation subgroups of UC with distinct associations to clinicopathological features and gene expression subtypes. We characterize three distinct patterns of differential methylation and use ENCODE data to show that tumor subgroup-defining DMRs display differential chromatin state, and regulatory factor binding preferences. Finally, we characterize an epigenetic switch involving the HOXA-genes with associations to tumor differentiation states and patient prognosis.

CONCLUSIONS

Genome-wide DMR methylation patterns are reflected in the gene expression subtypes of UC. UC DMRs display three distinct methylation patterns, each associated with intrinsic features of the genome and differential regulatory factor binding preferences. Epigenetic inactivation of HOX-genes correlates with tumor differentiation states and may present an actionable epigenetic alteration in UC.

Characterizing and controlling intrinsic biases of lambda exonuclease in nascent strand sequencing reveals phasing between nucleosomes and G-quadruplex motifs around a subset of human replication origins

Nascent strand sequencing (NS-seq) is used to discover DNA replication origins genome-wide, allowing identification of features for their specification. NS-seq depends on the ability of lambda exonuclease (λ-exo) to efficiently digest parental DNA while leaving RNA-primer protected nascent strands intact. We used genomics and biochemical approaches to determine if λ-exo digests all parental DNA sequences equally. We report that λ-exo does not efficiently digest G-quadruplex (G4) structures in a plasmid. Moreover, λ-exo digestion of nonreplicating genomic DNA (LexoG0) enriches GC-rich DNA and G4 motifs genome-wide. We used LexoG0 data to control for nascent strand–independent λ-exo biases in NS-seq and validated this approach at the rDNA locus. The λ-exo–controlled NS-seq peaks are not GC-rich, and only 35.5% overlap with 6.8% of all G4s, suggesting that G4s are not general determinants for origin specification but may play a role for a subset. Interestingly, we observed a periodic spacing of G4 motifs and nucleosomes around the peak summits, suggesting that G4s may position nucleosomes at this subset of origins. Finally, we demonstrate that use of Na⁺ instead of K⁺ in the λ-exo digestion buffer reduced the effect of G4s on λ-exo digestion and discuss ways to increase both the sensitivity and specificity of NS-seq.

ZNF143 provides sequence specificity to secure chromatin interactions at gene promoters

Chromatin interactions connect distal regulatory elements to target gene promoters guiding stimulus- and lineage-specific transcription. Few factors securing chromatin interactions have so far been identified. Here, by integrating chromatin interaction maps with the large collection of transcription factor-binding profiles provided by the ENCODE project, we demonstrate that the zinc-finger protein ZNF143 preferentially occupies anchors of chromatin interactions connecting promoters with distal regulatory elements. It binds directly to promoters and associates with lineage-specific chromatin interactions and gene expression. Silencing ZNF143 or modulating its DNA-binding affinity using single-nucleotide polymorphisms (SNPs) as a surrogate of site-directed mutagenesis reveals the sequence dependency of chromatin interactions at gene promoters. We also find that chromatin interactions alone do not regulate gene expression. Together, our results identify ZNF143 as a novel chromatin-looping factor that contributes to the architectural foundation of the genome by providing sequence specificity at promoters connected with distal regulatory elements.

Chromatin interactions can connect distal regulatory elements to promoters via protein factors, but few such factors have been identified. Here, the authors show that zinc-finger protein ZNF143 is a sequence-specific chromatin-looping factor that connects promoters with distal regulatory elements.

Multiple novel promoter-architectures revealed by decoding the hidden heterogeneity within the genome

An important question in biology is how different promoter-architectures contribute to the diversity in regulation of transcription initiation. A step forward has been the production of genome-wide maps of transcription start sites (TSSs) using high-throughput sequencing. However, the subsequent step of characterizing promoters and their functions is still largely done on the basis of previously established promoter-elements like the TATA-box in eukaryotes or the -10 box in bacteria. Unfortunately, a majority of promoters and their activities cannot be explained by these few elements. Traditional motif discovery methods that identify novel elements also fail here, because TSS neighborhoods are often highly heterogeneous containing no overrepresented motif. We present a new, organism-independent method that explicitly models this heterogeneity while unraveling different promoter-architectures. For example, in five bacteria, we detect the presence of a pyrimidine preceding the TSS under very specific circumstances. In tuberculosis, we show for the first time that the spacing between the bacterial 10-motif and TSS is utilized by the pathogen for dynamic gene-regulation. In eukaryotes, we identify several new elements that are important for development. Identified promoter-architectures show differential patterns of evolution, chromatin structure and TSS spread, suggesting distinct regulatory functions. This work highlights the importance of characterizing heterogeneity within high-throughput genomic data rather than analyzing average patterns of nucleotide composition.

Dynamic shifts in occupancy by TAL1 are guided by GATA factors and drive large-scale reprogramming of gene expression during hematopoiesis

We used mouse ENCODE data along with complementary data from other laboratories to study the dynamics of occupancy and the role in gene regulation of the transcription factor TAL1, a critical regulator of hematopoiesis, at multiple stages of hematopoietic differentiation. We combined ChIP-seq and RNA-seq data in six mouse cell types representing a progression from multilineage precursors to differentiated erythroblasts and megakaryocytes. We found that sites of occupancy shift dramatically during commitment to the erythroid lineage, vary further during terminal maturation, and are strongly associated with changes in gene expression. In multilineage progenitors, the likely target genes are enriched for hematopoietic growth and functions associated with the mature cells of specific daughter lineages (such as megakaryocytes). In contrast, target genes in erythroblasts are specifically enriched for red cell functions. Furthermore, shifts in TAL1 occupancy during erythroid differentiation are associated with gene repression (dissociation) and induction (co-occupancy with GATA1). Based on both enrichment for transcription factor binding site motifs and co-occupancy determined by ChIP-seq, recruitment by GATA transcription factors appears to be a stronger determinant of TAL1 binding to chromatin than the canonical E-box binding site motif. Studies of additional proteins lead to the model that TAL1 regulates expression after being directed to a distinct subset of genomic binding sites in each cell type via its association with different complexes containing master regulators such as GATA2, ERG, and RUNX1 in multilineage cells and the lineage-specific master regulator GATA1 in erythroblasts.

Chromatin structure is distinct between coding and non-coding single nucleotide polymorphisms

Background

Previous studies suggested that nucleosomes are enriched with single nucleotide polymorphisms (SNPs) in humans and that the occurrence of mutations is closely associated with CpG dinucleotides. We aimed to determine if the chromatin organization is genomic locus specific around SNPs, and if newly occurring mutations are associated with SNPs.

Results

Here, we classified SNPs according their loci and investigated chromatin organization in both CD4⁺ T cell and lymphoblastoid cell in humans. We calculated the SNP frequency around somatic mutations. The results indicated that nucleosome occupancy is different around SNPs sites in different genomic loci. Coding SNPs are mainly enriched at nucleosomes and associated with repressed histone modifications (HMs) and DNA methylation. Contrastingly, intron SNPs occur in nucleosome-depleted regions and lack HMs. Interestingly, risk-associated non-coding SNPs are also enriched at nucleosomes with HMs but associated with low GC-content and low DNA methylation level. The base-transversion allele frequency is significantly low in coding-synonymous SNPs (P < 10^-11). Another finding is that at the -1 and +1 positions relative to the somatic mutation sites, the SNP frequency was significantly higher (P < 3.2 × 10^-5).

Conclusions

The results suggested chromatin structure is different around coding SNPs and non-coding SNPs. New mutations tend to occur at the -1 and +1 position immediately near the SNPs.

A Simple but Powerful Heuristic Method for Accelerating k-Means Clustering of Large-Scale Data in Life Science

K-means clustering has been widely used to gain insight into biological systems from large-scale life science data. To quantify the similarities among biological data sets, Pearson correlation distance and standardized Euclidean distance are used most frequently; however, optimization methods have been largely unexplored. These two distance measurements are equivalent in the sense that they yield the same k-means clustering result for identical sets of k initial centroids. Thus, an efficient algorithm used for one is applicable to the other. Several optimization methods are available for the Euclidean distance and can be used for processing the standardized Euclidean distance; however, they are not customized for this context. We instead approached the problem by studying the properties of the Pearson correlation distance, and we invented a simple but powerful heuristic method for markedly pruning unnecessary computation while retaining the final solution. Tests using real biological data sets with 50-60K vectors of dimensions 10-2001 (~400 MB in size) demonstrated marked reduction in computation time for k = 10-500 in comparison with other state-of-the-art pruning methods such as Elkan's and Hamerly's algorithms. The BoostKCP software is available at http://mlab.cb.k.u-tokyo.ac.jp/~ichikawa/boostKCP/.

Nucleosome organization in the vicinity of transcription factor binding sites in the human genome

Background

The binding of transcription factors (TFs) to specific DNA sequences is an initial and crucial step of transcription. In eukaryotes, this process is highly dependent on the local chromatin state, which can be modified by recruiting chromatin remodelers. However, previous studies have focused mainly on nucleosome occupancy around the TF binding sites (TFBSs) of a few specific TFs. Here, we investigated the nucleosome occupancy profiles around computationally inferred binding sites, based on 519 TF binding motifs, in human GM12878 and K562 cells.

Results

Although high nucleosome occupancy is intrinsically encoded at TFBSs in vitro, nucleosomes are generally depleted at TFBSs in vivo, and approximately a quarter of TFBSs showed well-positioned in vivo nucleosomes on both sides. RNA polymerase near the transcription start site (TSS) has a large effect on the nucleosome occupancy distribution around the binding sites located within one kilobase to the nearest TSS; fuzzier nucleosome positioning was thus observed around these sites. In addition, in contrast to yeast, repressors, rather than activators, were more likely to bind to nucleosomal DNA in the human cells, and nucleosomes around repressor sites were better positioned in vivo. Genes with repressor sites exhibiting well-positioned nucleosomes on both sides, and genes with activator sites occupied by nucleosomes had significantly lower expression, suggesting that actions of activators and repressors are associated with the nucleosome occupancy around their binding sites. It was also interesting to note that most of the binding sites, which were not in the DNase I-hypersensitive regions, were cell-type specific, and higher in vivo nucleosome occupancy were observed at these binding sites.

Conclusions

This study demonstrated that RNA polymerase and the functions of bound TFs affected the local nucleosome occupancy around TFBSs, and nucleosome occupancy patterns around TFBSs were associated with the expression levels of target genes.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-493) contains supplementary material, which is available to authorized users.

Comparison of REST cistromes across human cell types reveals common and context-specific functions

Recent studies have shown that the transcriptional functions of REST are much broader than repressing neuronal genes in non-neuronal systems. Whether REST occupies similar chromatin regions in different cell types and how it interacts with other transcriptional regulators to execute its functions in a context-dependent manner has not been adequately investigated. We have applied ChIP-seq analysis to identify the REST cistrome in human CD4+ T cells and compared it with published data from 15 other cell types. We found that REST cistromes were distinct among cell types, with REST binding to several tumor suppressors specifically in cancer cells, whereas 7% of the REST peaks in non-neuronal cells were ubiquitously called and <25% were identified for ≥ 5 cell types. Nevertheless, using a quantitative metric directly comparing raw ChIP-seq signals, we found the majority (∼80%) was shared by ≥ 2 cell types. Integration with RNA-seq data showed that REST binding was generally correlated with low gene expression. Close examination revealed that multiple contexts were correlated with reduced expression of REST targets, e.g., the presence of a cognate RE1 motif and cellular specificity of REST binding. These contexts were shown to play a role in differential corepressor recruitment. Furthermore, transcriptional outcome was highly influenced by REST cofactors, e.g., SIN3 and EZH2 co-occupancy marked higher and lower expression of REST targets, respectively. Unexpectedly, the REST cistrome in differentiated neurons exhibited unique features not observed in non-neuronal cells, e.g., the lack of RE1 motifs and an association with active gene expression. Finally, our analysis demonstrated how REST could differentially regulate a transcription network constituted of miRNAs, REST complex and neuronal factors. Overall, our findings of contexts playing critical roles in REST occupancy and regulatory outcome provide insights into the molecular interactions underlying REST's diverse functions, and point to novel roles of REST in differentiated neurons.

Exploring the function of genetic variants in the non-coding genomic regions: approaches for identifying human regulatory variants affecting gene expression

Understanding the genetic basis of human traits/diseases and the underlying mechanisms of how these traits/diseases are affected by genetic variations is critical for public health. Current genome-wide functional genomics data uncovered a large number of functional elements in the noncoding regions of human genome, providing new opportunities to study regulatory variants (RVs). RVs play important roles in transcription factor bindings, chromatin states and epigenetic modifications. Here, we systematically review an array of methods currently used to map RVs as well as the computational approaches in annotating and interpreting their regulatory effects, with emphasis on regulatory single-nucleotide polymorphism. We also briefly introduce experimental methods to validate these functional RVs.

Chromatin features and the epigenetic regulation of pluripotency states in ESCs

In pluripotent stem cells, the interplay between signaling cues, epigenetic regulators and transcription factors orchestrates developmental potency. Flexibility in gene expression control is imparted by molecular changes to the nucleosomes, the building block of chromatin. Here, we review the current understanding of the role of chromatin as a plastic and integrative platform to direct gene expression changes in pluripotent stem cells, giving rise to distinct pluripotent states. We will further explore the concept of epigenetic asymmetry, focusing primarily on histone stoichiometry and their associated modifications, that is apparent at both the nucleosome and chromosome-wide levels, and discuss the emerging importance of these asymmetric chromatin configurations in diversifying epigenetic states and their implications for cell fate control.

Probabilistic partitioning methods to find significant patterns in ChIP-Seq data

MOTIVATION

We have witnessed an enormous increase in ChIP-Seq data for histone modifications in the past few years. Discovering significant patterns in these data is an important problem for understanding biological mechanisms.

RESULTS

We propose probabilistic partitioning methods to discover significant patterns in ChIP-Seq data. Our methods take into account signal magnitude, shape, strand orientation and shifts. We compare our methods with some current methods and demonstrate significant improvements, especially with sparse data. Besides pattern discovery and classification, probabilistic partitioning can serve other purposes in ChIP-Seq data analysis. Specifically, we exemplify its merits in the context of peak finding and partitioning of nucleosome positioning patterns in human promoters.

AVAILABILITY AND IMPLEMENTATION

The software and code are available in the supplementary material.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.