Skip to main content

A new class of constitutively active super-enhancers is associated with fast recovery of 3D chromatin loops

Abstract

Background

Super-enhancers or stretch enhancers are clusters of active enhancers that often coordinate cell-type specific gene regulation during development and differentiation. In addition, the enrichment of disease-associated single nucleotide polymorphism in super-enhancers indicates their critical function in disease-specific gene regulation. However, little is known about the function of super-enhancers beyond gene regulation.

Results

In this study, through a comprehensive analysis of super-enhancers in 30 human cell/tissue types, we identified a new class of super-enhancers which are constitutively active across most cell/tissue types. These ‘common’ super-enhancers are associated with universally highly expressed genes in contrast to the canonical definition of super-enhancers that assert cell-type specific gene regulation. In addition, the genome sequence of these super-enhancers is highly conserved by evolution and among humans, advocating their universal function in genome regulation. Integrative analysis of 3D chromatin loops demonstrates that, in comparison to the cell-type specific super-enhancers, the cell-type common super-enhancers present a striking association with rapidly recovering loops.

Conclusions

In this study, we propose that a new class of super-enhancers may play an important role in the early establishment of 3D chromatin structure.

Background

Super-enhancers or stretch enhancers are defined by a strong enrichment of mediators and transcription-regulating proteins, appearing to play a deterministic role in cellular identity by controlling the expression of cell-type specific genes [1, 2]. Previous studies have revealed the critical function of super-enhancers during development and differentiation [3]. The enrichment of disease-associated single nucleotide polymorphism (SNP) in super-enhancers compared to that of typical enhancers proposed a substantial link between super-enhancers and many complex human diseases [2]. In addition, recent studies have proposed potential functions of super-enhancers in the extremely long-range chromatin communications and the establishment of 3D chromatin loops [4]. These results suggest a universal role of super-enhancers in genome regulation apart from cell-type specific gene regulation, but little is known about the mechanisms underlying these various functions. To extend the current knowledge of super-enhancers and their biological roles, we conducted a comprehensive analysis of super-enhancer activities across 30 human cell/tissue types. Our analysis suggests that a substantial number of super-enhancers exhibits prevalent activities across cell-types in terms of H3K27ac signals and that these non-canonical super-enhancers are involved in the formation of fast recovering chromatin loops.

Results

Genomic landscape of super-enhancer domains

To characterize different modes of super-enhancers across various human cell/tissue types, we first investigated the genomic distribution of super-enhancers across 30 human cell/tissue types, including 5 H1-derived early lineages, 8 immortalized cell lines, and 17 human postmortem tissues. Super-enhancers were defined using ROSE algorithm [2, 5] based on H3K27ac ChIP-seq signal (see Methods), which resulted in 17,916 super-enhancers with an average of 597 super-enhancers per cell/tissue type. Merging the super-enhancers in all of the cell/tissue types resulted in 6039 putative super-enhancer regions, named super-enhancer domains (SE domains) (Fig. 1a, see Methods). The super-enhancer domains presented a mean length of 32Kb, which was slightly longer than the individual super-enhancers with the mean length of 26Kb. Excluding any one of 30 samples did not significantly affect the overall distribution of length of super-enhancer domains (one-way Analysis of Variance, p-value > 0.5), which excludes the possibility of a critical bias in the definition of super-enhancer domains caused by one or more of samples. In total, around 6.32% of the human genome was marked by super-enhancer domains, and 48.3% of the super-enhancer domains consisted of multiple super-enhancers identified in at least two or more cell/tissue types analyzed, suggesting a recurrent formation of super-enhancers in specific genomic regions.

Fig. 1
figure 1

Super-enhancer domains can be grouped into distinct categories. a A definition of super-enhancer domains. b A smooth scatter plot between tau score (x-axis) and the number of cell/tissue types associated with super-enhancer domains. Vertical arrowheads indicate 20 and 80% percentile of tau score. Red, blue, and green colored box indicate unique, non-unique and common super-enhancer domains, respectively. c-e Heatmaps of background-subtracted H3K27ac RPM of super-enhancer domains across cell/tissue types. Color intensity indicates log2 transformed H3K27ac RPM of super-enhancer domains. f-h Genome browser snapshots showing H3K27ac signals of each class of super-enhancer domain examples. Unique super-enhancer domain nearby SOX1, specific to neuronal progenitor cells (f). Non-unique super-enhancer domain nearby TJP3, specific to digestive system related tissues (g). Common super-enhancer domain nearby MALAT1 and NEAT1 (h). H3K27ac signal is uniformly scaled in the range of 0–20. Cell/tissue types are grouped according to their lineage or functional similarity

Identification of common super-enhancer domains

In contrast to the previous notion that super-enhancers are highly cell/tissue-type specific, some super-enhancer domains exhibited a strong enrichment of H3K27ac signals in a surprisingly large number of cell/tissue types. To assess cell/tissue-type specificity of the super-enhancer domains, we utilized tau score [6], a measurement of tissue-specificity commonly used in gene expression studies (see Methods). Tau score of 0 indicates ubiquitous super-enhancer domain activity, while 1 represents tissue-specific activation of the domain. The super-enhancer domain activities based on H3K27ac signals showed a wide spectrum of tau scores, reflecting a varying degree of cell/tissue-type specificity for each super-enhancer domain (Fig. 1b). It should be noted that some super-enhancer domains were annotated by only one cell/tissue type, but showed low tau scores. Through manual inspection, we found that these domains presented moderately high H3K27ac signals in other cell/tissue types (See Additional file 1: Figure S1), but the signals were insufficient to pass the threshold for super-enhancer calling by a small margin. For this reason, tau score was used as the main factor to define super-enhancer specificity.

In order to categorize super-enhancer domains with different cell/tissue-type specificity, we defined cell/tissue-type common and cell/tissue-type unique super-enhancer domains. The domains located within 20% or lower percentile of tau score were defined as common super-enhancer domains (1208 domains, 20.0%, Fig. 1c). On the other hand, domains that were called only in one cell/tissue type and located at 80% or higher percentile of tau score were defined as unique super-enhancer domains (1098 domains, 18.2%, Fig. 1e). Remaining domains were defined to be non-unique super-enhancer domains (3733 domains, 61.8%, Fig. 1d). To ensure that cancer cell lines or other immortalized cell lines did not produce critical bias by producing a large portion of common super-enhancers, we conducted hierarchical clustering between samples using the H3K27ac signals in the super-enhancer domains of available cell/tissue types (See Additional file 1: Figure S2). The result demonstrated that, while tissue samples with close relatedness grouped together, cancer cell lines or other immortalized cell lines were not concentrated in the same cluster.

A notable example of unique super-enhancer domains was found in neural progenitor cells derived from human embryonic stem cells (H1-ESC). Two super-enhancer domains were uniquely defined in neural progenitor cells at downstream of SOX1 (Fig. 1f), a gene known to be highly expressed only in neural progenitor cells and to play an important role in nervous system development. Non-unique super-enhancer domains showed high activity in several cell/tissue types but low activity in the rest. For example, a non-unique super-enhancer domain was called for tissues associated with the digestive system such as stomach, small bowel, and sigmoid colon (Fig. 1g). There was also a high H3K27ac signal in the same domain in esophagus tissue, although it was not called as a super-enhancer. The putative target gene of the domain is TJP3 (tight junction protein 3), which is involved in junctional integrity of the intestinal cells. This example suggests a possibility that the formation of non-unique super-enhancer domains is a key regulator for shared functions in several tissues. Common super-enhancer domains showed high ubiquitous H3K27ac signals across the cell/tissue types, discordant with the previous notion of cell/tissue-type specificity of super-enhancers. An example of common super-enhancer domains was found near the noncoding RNA genes, NEAT1 and MALAT1, which are highly expressed across all cell/tissue types and previously suggested to function in a general biological process, such as the association with nuclear speckle formation [7] (Fig. 1h).

Distinct biological function of target genes in common super-enhancer domains

Super-enhancers are presumably associated with cell/tissue-type specific gene regulation, but the identification of distinct classes of super-enhancer domains raises the possibility that distinct modes of action exist for super-enhancers. We investigated the expression patterns of putative target genes in each class of super-enhancer domains. As expected, putative target genes of unique super-enhancer domains were highly expressed in the corresponding cell/tissue type (Fig. 2a and b, see Methods). In contrast, genes associated with common super-enhancer domains showed universally high expression (Fig. 2a and b). Although there was a significant enrichment of housekeeping genes [8] in putative target genes of common super-enhancer domains, the majority (82.4%) were non-housekeeping genes, indicating that housekeeping genes cannot fully explain the function of common super-enhancer domains. Gene set enrichment analysis of genes associated with unique super-enhancer domains showed enrichment in cellular identity functions in the corresponding cell/tissue type (Table 1, see Methods), whereas genes associated with common super-enhancer domains were enriched in basic cellular functions, such as transcriptional regulation, cell motility, and regulation of cell proliferation (Table 2). The overall pattern of expression and the significance of enriched pathways were moderately affected but largely stayed consistent, when an alternative definition of putative target genes was used including the nearest 3 genes from the midpoint of each super-enhancer domain (See Additional file 1: Figure S3a, b and Table S1, S2). Our analysis revealed distinct biological functional enrichment in putative target genes between common and unique super-enhancer domains.

Fig. 2
figure 2

Functional characterization of super-enhancer domain classes. a Heatmaps showing log2(RPKM+ 1) value of putative target gene of unique super-enhancer domains (top), non-unique super-enhancer domains (middle), and common super-enhancer domains (bottom). b Histogram of tau score for putative target gene expression in each class of super-enhancer domains. Low tau score indicates a universal expression pattern

Table 1 GO analysis of putative target genes of unique super-enhancer domains in lymphoblast and prefrontal cortex
Table 2 GO analysis of putative target genes of common super-enhancer domains

Genomic properties of common super-enhancer domains

To reason that super-enhancers play a key role in gene regulation, one would expect to see a strong correlation between the activity of a super-enhancer domain and the expression of its putative target gene. In general, putative target gene expression was significantly correlated with the activities of super-enhancer domains (one-sample t-test, p-value < 10− 15 for common, non-unique, and unique super-enhancer domains) (Fig. 3a, see Methods). However, common super-enhancer domains showed the lowest correlation coefficients, when compared to unique and non-unique super-enhancer domains (KS test, *** p-value < 10− 15). For example, a common super-enhancer domain shows considerable activities with enriched H3K27ac signals, but its putative target gene, MLL5, does not show a concordant expression level with the super-enhancer domain activities in all of the cell/tissue types (Fig. 3b). This result raised the possibility that common super-enhancer domains are responsible for an additional biological role, other than the expression regulation of genes within their close proximity.

Fig. 3
figure 3

Genomic characterization of common super-enhancer domains. a Boxplots of correlation between log2(H3K27ac RPM + 1) in super-enhancer domains and RPKM of putative target genes (n = 1208 for common SE domains, n = 3733 for non-unique SE domains, n = 1098 for unique SE domains). Pearson correlation coefficient was calculated for each super-enhancer domain across 24 cell/tissue types. Significance of difference between super-enhancer domain categories was evaluated by KS test (*** p-value < 5 × 10− 12). b An example of common super-enhancer domain showing high H3K27ac signal across cell/tissue types but showing irregular RNA expression level. Left panel presents genome browser snapshot of H3K27ac signals at super-enhancer domain. Right panel presents RNA-seq RPKM of MLL5 gene. Pearson correlation coefficient (PCC) between H3K27ac signal at a super enhancer domain and the MLL5 RPKM is shown below. c Boxplots of mean sequence conservation score for each super-enhancer domain (n = 1098 for unique SE domains, n = 3733 for non-unique SE domains, n = 1208 for common SE domains). Significance of difference between super-enhancer domain categories was evaluated with KS test (*** p-value < 1.1 × 10− 5 between random and others, p-value < 5.1 × 10− 12 for the rest). d Boxplots of mean CDTS for each super-enhancer domain (n = 1098 for unique SE domains, n = 3733 for non-unique SE domains, n = 1208 for common SE domains). Significance of difference between super-enhancer domain categories is evaluated with KS test (*** p-value < 1.7 × 10− 9 for unique and random, p-value < 3.0 × 10− 14 for the rest). e An example of adjacent non-unique super-enhancer domain and common super-enhancer domain illustrating evolutionary conservation of common super-enhancer domain. Shown together are 100-way PhastCons and CDTS. 100-way PhastCons track has the scale from 0 to 1, and CDTS track is scaled from − 20 to + 20. Common and unique super-enhancer domain are shown as right and left bar, respectively, in ‘Super-enhancer domains’ annotation

Further investigation revealed that common super-enhancer domains are often located in evolutionarily conserved genomic regions compared to other super-enhancer domains and random expectation (Fig. 3c, see Methods). However, we did not observe a strong conservation score in the unique super-enhancer domains compared to random expectation. Similarly, genomic areas of common super-enhancer domains also tend to have a low context-dependent tolerance score (CDTS) [9], which represents the frequency of sequence variation across the human population (Fig. 3d, see Methods). The lower CDTS indicates a sequence is more resistant to the expected variation. Thus, low CDTS for common super-enhancer domains suggests, again, a potential function of these domains in universal genome regulation. For example, a common super-enhancer domain nearby NR2F2 gene presents both a high evolutionary conservation score and a low CDTS compared to the adjacent non-unique super-enhancer domain (Fig. 3e).

A striking association of common super-enhancer domains during early establishment of 3D chromatin loops

A previous study has revealed an association between super-enhancers and 3D chromatin structure [4], where 3D chromatin loops have a wide range of recovery rate following cohesin degradation. Loop domains with a fast recovery rate contained a greater number of super-enhancers compared to slow-recovering loops. We surmised that common super-enhancer domains play a significant role in the loop recovery since rapid recovery of 3D chromatin structure could be more critical in the universal genome regulation compared to the cell/tissue-type specific regulation.

To test our hypothesis, we examined the enrichment of 336 super-enhancer domains called from HCT-116 human colorectal cancer cell lines in both loop anchors (two contact loci) and loop domains (between two anchors). For each time point from 20 min to 180 min after the degradation of cohesin and subsequent start of recovery, we classified chromatin loops into fast, slow, and moderate-recovery loops by the intensity percentile (top 10% defined as fast recovery loop and bottom 10% as slow recovery loop) (Fig. 4a upper panel, see Methods). The overlap between super-enhancer domain and loop anchor or loop domain was calculated as the measurement of super-enhancer domain enrichment (Fig. 4a lower panel, see Methods). At an early time point (20 min after cohesin recovery), fast recovering loop anchors and domains were strikingly enriched in common super-enhancer domains but not in unique super-enhancer domains (Fig. 4b upper panels). In contrast, slow-recovering loop anchors were significantly enriched in unique super-enhancer domains but not in common super-enhancer domains. In order to confirm the result not only for super-enhancer domains of HCT-116 but also of the 30 cell/tissue types, we further expanded the analysis and observed concordant results (Fig. 4b, lower panels). This result supports the cell/tissue-type independent role of common super-enhancers in 3D chromatin loop formation.

Fig. 4
figure 4

Common super-enhancer domains are associated with the early establishment of 3D chromatin loops. a The schematic overview of calculated loop recovery rate (upper panel) and super-enhancer domain enrichment (lower panel). b Heatmaps showing mean enrichment value of super-enhancer domains called from HCT-116 cell line (upper two panels) and all samples (lower two panels) at 20 min after cohesin recovery in loop anchors (left panels) and loop domains (right panels). Empirical p-values were calculated through random expectation values based on 10,000 iterations (* p-value < 0.05, ** p-value < 0.01, *** p-value < 0.001). c Bar plots of super-enhancer domain enrichment in the top 10% fast recovery loops for each time point in loop anchors (left) and loop domains (right) (two-sample t-test, * p-value < 0.05, ** p-value < 0.01, *** p-value < 0.001). d Model diagram of modes of super-enhancer functions. Canonical cell-type specific gene regulation (left) and 3D genome regulation of common super-enhancer domains (right)

We expanded the analysis to include later time points (see Methods) and revealed that enrichment of common super-enhancer domains at both loop anchors and domains was the highest at the earliest time point (20 min), and gradually decreased in later time points (Fig. 4c). On the other hand, enrichment of unique super-enhancer domains significantly increased in the later time points (two-sample t-test, * p-value < 0.05, ** p-value < 0.01, *** p-value < 0.001). Our result implies that common super-enhancer domains are often associated with the initial establishment of 3D chromatin loops compared to other super-enhancer domains.

Discussion and conclusion

We present a systematic approach to assess and characterize super-enhancers in a wide array of human cell/tissue types. The effective computation of background normalized super-enhancer activity and the use of simple, yet powerful mathematical expression to evaluate the tissue specificity have led us to identify a new class of super-enhancers. These super-enhancer domains exhibit a high degree of universality across many cell/tissue types – a novel aspect of super-enhancers that has never been spotlighted in previous literature. The universally active super-enhancers, which we defined as common super-enhancer domains, display a strong association with fast recovering chromatin loops after a sequential cohesin removal and restoration. Although the enrichment of super-enhancers at fast-recovering loops is previously reported [4], our analysis further reveals that common super-enhancer domains are strikingly enriched in the fast-recovering loops by more than 12-folds, compared to unique super-enhancer domains. The implications that our data bring forth to the function of super-enhancers in shaping chromatin organization may be outlined as follows. First, common super-enhancers may facilitate recruiting a substantial amount of structural proteins, such as cohesin and CTCF, to expedite loop recovery. A sequential model is also an interesting possibility, where the cell-type nonspecific formation of fast recovering loops facilitated by common super-enhancer guide the genome folding in a stepwise manner to lead following cell-type specific loop formation. A hierarchical model where cell-type specific loop formation by unique super-enhancers requires the higher-order loop formation promoted by common super-enhancers is also reasonable. Further experiments will be required to test our hypothesis for the validation.

To conclude, we propose two distinct modes of super-enhancers; one is cell-type specific gene regulation mainly mediated by unique super-enhancer domains and second is 3D genome regulation mediated by common super-enhancer domains as a non-canonical function (Fig. 4d). Although there is a limitation to generalize our conclusion since we investigated only one cell type, we shed light on common super-enhancers as a new potential mechanism underlying the early establishment of 3D chromatin loops.

Methods

Human cell/tissue types

In this study, 30 human cell/tissue types were examined. This includes H1 embryonic stem cell and its derived cell types (H1 embryonic stem cell, mesendoderm, mesenchymal stem cell, neuronal progenitor cell, and trophoblast), 8 cell lines (GM12878, HUVEC, IMR90, K562, NHEK, HMEC, HEK293, and HCT-116), and 17 tissue types (adrenal gland, aorta, bladder, dorsolateral prefrontal cortex, esophagus, gastric, lung, ovary, pancreas, psoas, right atrium, right ventricle, liver, sigmoid colon, small bowel, spleen, and thymus).

Super-enhancer call

H3K27ac and input ChIP-seq reads were downloaded from Roadmap Epigenomics, ENCODE project and GEO database [10, 11]. Sources are specified in Additional file 1: Table S3. Reads were aligned to the human reference genome (hg19 assembly) using BWA-mem [12]. Unmapped and poorly mapped (MAPQ < 10) reads were removed. PCR duplicates were also removed using Picard Markduplicates. The remaining reads were then used to call peaks using MACS2 [13]. Default parameters were used except q-value 0.01 was replaced with the p-value less than 1E-05. Super-enhancers were defined by stitching peaks using ROSE [2, 5] with default parameters, except TSS exclusion zone size was adjusted to 2500 bp. TSS exclusion was used because the H3K27ac signal is enriched for both active enhancers and promoters. Super-enhancers that span ENCODE blacklist regions were discarded.

Calling super-enhancer domains

Super-enhancer domains were defined by merging the super-enhancers called for 30 samples using BEDTools (see Fig. 1a). Super-enhancers with any overlap were merged into one super-enhancer domain.

Cell/tissue-type specificity measurement with tau value

Tau score was utilized to assess the cell/tissue-type specificity of super-enhancer domain activities and gene expression. Tau score is a measure commonly used to quantify tissue-specificity of expression as it enables a threshold-free evaluation of cell/tissue-type specificity. Tau score ranges from 0 to 1, where 0 indicates general, and 1 indicates specific. Tau score is calculated as \( \tau =\frac{\sum_{i=1}^n\left(1-{\widehat{x}}_i\right)}{n-1} \); \( {\widehat{x}}_i=\frac{x_i}{\underset{1\le \mathrm{i}\le \mathrm{n}}{\max}\left({x}_i\right)} \) [6], where x indicates cell/tissue-type specificity and n denotes the number of cell/tissue types.

Gene expression data

Uniformly processed and quantified RPKM of protein-coding genes for 10 cell lines (H1 embryonic stem cell, mesendoderm, trophoblast, mesenchymal stem cell, neuronal progenitor cell, GM12878, HMEC, HUVEC, K562, and NHEK) and 14 tissues (aorta, esophagus, gastric, liver, lung, ovary, pancreas, psoas, right atrium and ventricle, sigmoid colon, small intestine, thymus, and spleen) was obtained from Roadmap epigenomics (“RNA-seq uniform processing and quantification for consolidated epigenomes”). The data had been normalized between cell/tissue types to have the same number of aligned reads on coding exons. Super-enhancer domains called for none of the 24 cell/tissue types are excluded in gene expression analysis.

Target gene assignment of super-enhancer domains

Putative target gene of a super-enhancer domain was defined as the gene with the closest transcriptional start site to the midpoint of the super-enhancer domain. For extended analysis, alternative definition utilized the 3 genes with transcriptional start site closest to the midpoint of super-enhancer domains.

GO enrichment analysis using DAVID

For putative target genes of super-enhancer domains, GO terms were analyzed using DAVID 6.8 [14]. Significant GO_BP annotations were selected using p-value corrected by Benjamini-Hochberg procedure with threshold 0.05 and sorted by the p-value. For the analysis of the putative target genes of unique super-enhancer domains, we selected lymphoblast and prefrontal cortex as examples. Due to the small number of putative target genes of unique super-enhancer domains in each cell type, p-values tend to be insignificant. Thus, we presented up to top 10 enriched GO_BP annotations sorted with p-values without cutoff.

Housekeeping genes

The putative target genes of super-enhancer domains were compared to the list of housekeeping genes provided by Eisenberg, E. et al. [8].

Correlation between SE domain activity and expression level of its putative target gene

Pearson correlation coefficients between background-subtracted log2(RPM + 1) value of H3K27ac ChIP-seq signal on super-enhancer domains and RPKM of their putative target gene were calculated across 24 cell/tissue. Super-enhancer domains called for none of the 24 cell/tissue types were excluded in the correlation analysis.

Conservation and tolerance scores of super-enhancer domains

Evolutionary sequence conservation score calculated from multiple alignments of 100 vertebrates using phastCons method was obtained from the UCSC annotation database (phastCons100way). SumData column of 24 bp resolution conservation score was mapped onto super-enhancer domains, and the mean value of mapped conservation score was used. CDTS referring to the variability of a nucleotide calculated for 11,257 human genomes was obtained from a previous study [9]. CDTS of the 1 bp resolution was mapped onto super-enhancer domains, and the mean of mapped CDTS was used.

Hi-C matrix of HCT-116 following cohesin degradation and subsequent recovery

Hi-C matrix data of HCT-116 with a 25 kb resolution for control, cohesin degraded, and 20, 40, 60, and 180 min since subsequent recovery of cohesion, were obtained from the previous study (GEO accession number GSE104334) [4]. Each matrix was normalized by Knight-Ruiz algorithm. As the original read depth of each matrix was different, we normalized matrices to have the same total interaction frequency within 100 kb distance since short-range interactions are known to be not affected by cohesin degradation. Chromatin loops were identified using HiCCUPS [15] with default parameters for medium resolution matrix. False positive chromatin loops from karyotypic abnormalities of HCT-116 cell line were filtered. The observed/expected enrichment filter of 4.5 was applied, as suggested by the previous study.

Super-enhancer enrichment in loop anchors and domains

The recovery rate of loop anchors (interacting loci) and domains (region between two anchors) was defined as the ratio of loop strength at each time point compared to the strength in untreated cells (see Fig. 4a, upper panel). Loops with monotonically increasing strength with regard to the time were selected for the analysis. Loop anchors were padded by 5 kb in obtaining the super-enhancer domain enrichment. “Fraction” value of super-enhancer domains in loop anchors or domains was calculated by dividing the density of overlapping super-enhancer domains by the expected density (see Fig. 4a, lower panel). The expected density of each class of super-enhancer domain is the proportion of the corresponding domains in the whole hg19 genome. The expected value of fraction for randomly shuffled super-enhancer domains is 1, with an infinite number of trials. The super-enhancer domain enrichment of a loop is obtained by dividing the observed fraction value by the random expectation of the fraction. The random expectation of fraction is the average of fractions calculated for 10,000 shuffled super-enhancer domains. Empirical p-value of a fraction where the distribution is defined as the result of 10,000 trials, is utilized as the significance measure of the enrichment. When comparing two enrichments, the p-value of two-sample t-test was utilized as the significance measure.

Abbreviations

CDTS:

Context-dependent tolerance score

ChIP-seq:

Chromatin immunoprecipitation sequencing

H3K27ac:

Histone H3 lysine 27 acetylation

KS test:

Kolmogorov-Smirnov test

RPKM:

Reads per kilobase million

RPM:

Reads per million mapped reads

SE domain:

Super-enhancer domain

TSS:

Transcription start site

References

  1. Parker SC, Stitzel ML, Taylor DL, Orozco JM, Erdos MR, Akiyama JA, et al. Chromatin stretch enhancer states drive cell-specific gene regulation and harbor human disease risk variants. Proc Natl Acad Sci U S A. 2013;110(44):17921–6.

    Article  CAS  Google Scholar 

  2. Whyte WA, Orlando DA, Hnisz D, Abraham BJ, Lin CY, Kagey MH, et al. Master transcription factors and mediator establish super-enhancers at key cell identity genes. Cell. 2013;153(2):307–19.

    Article  CAS  Google Scholar 

  3. Adam RC, Yang H, Rockowitz S, Larsen SB, Nikolova M, Oristian DS, et al. Pioneer factors govern super-enhancer dynamics in stem cell plasticity and lineage choice. Nature. 2015;521(7552):366–70.

    Article  CAS  Google Scholar 

  4. Rao SSP, Huang SC, Glenn St Hilaire B, Engreitz JM, Perez EM, Kieffer-Kwon KR, et al. Cohesin Loss Eliminates All Loop Domains. Cell 2017;171(2):305–320 e24.

  5. Loven J, Hoke HA, Lin CY, Lau A, Orlando DA, Vakoc CR, et al. Selective inhibition of tumor oncogenes by disruption of super-enhancers. Cell. 2013;153(2):320–34.

    Article  CAS  Google Scholar 

  6. Kryuchkova-Mostacci N, Robinson-Rechavi M. A benchmark of gene expression tissue-specificity metrics. Brief Bioinform. 2017;18(2):205–14.

    CAS  PubMed  Google Scholar 

  7. West JA, Davis CP, Sunwoo H, Simon MD, Sadreyev RI, Wang PI, et al. The long noncoding RNAs NEAT1 and MALAT1 bind active chromatin sites. Mol Cell. 2014;55(5):791–802.

    Article  CAS  Google Scholar 

  8. Eisenberg E, Levanon EY. Human housekeeping genes, revisited. Trends Genet. 2013;29(10):569–74.

    Article  CAS  Google Scholar 

  9. di Iulio J, Bartha I, Wong EHM, Yu HC, Lavrenko V, Yang D, et al. The human noncoding genome defined by genetic diversity. Nat Genet. 2018;50(3):333–7.

    Article  Google Scholar 

  10. Roadmap Epigenomics C, Kundaje A, Meuleman W, Ernst J, Bilenky M, Yen A, et al. Integrative analysis of 111 reference human epigenomes. Nature. 2015;518(7539):317–30.

    Article  Google Scholar 

  11. Leung D, Jung I, Rajagopal N, Schmitt A, Selvaraj S, Lee AY, et al. Integrative analysis of haplotype-resolved epigenomes across human tissues. Nature. 2015;518(7539):350–4.

    Article  CAS  Google Scholar 

  12. Li H, Durbin R. Fast and accurate short read alignment with burrows-wheeler transform. Bioinformatics. 2009;25(14):1754–60.

    Article  CAS  Google Scholar 

  13. Zhang Y, Liu T, Meyer CA, Eeckhoute J, Johnson DS, Bernstein BE, et al. Model-based analysis of ChIP-Seq (MACS). Genome Biol. 2008;9(9):R137.

    Article  Google Scholar 

  14. Huang d W, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4(1):44–57.

    Article  CAS  Google Scholar 

  15. Durand NC, Shamim MS, Machol I, Rao SS, Huntley MH, Lander ES, et al. Juicer provides a one-click system for analyzing loop-resolution hi-C experiments. Cell Syst. 2016;3(1):95–8.

    Article  CAS  Google Scholar 

Download references

Acknowledgments

We thank all members of Dr. Jung’s laboratory for their comments and criticisms on this study.

Funding

Ministry of Science, ICT, and Future Planning through the National Research Foundation in Republic of Korea [2017R1C1B2008838 to I. Jung]; Korean Ministry of Health and Welfare [HI17C0328 to I. Jung]. This work was supported by the KAIST Future Systems Healthcare Project from the Ministry of Science and ICT. Publication charges come from HI17C0328.

Availability of data and materials

All data generated or analyzed during this study are included in this published article.

About this supplement

This article has been published as part of BMC Bioinformatics Volume 20 Supplement 3, 2019: Selected articles from the 17th Asia Pacific Bioinformatics Conference (APBC 2019): bioinformatics. The full contents of the supplement are available online at https://bmcbioinformatics.biomedcentral.com/articles/supplements/volume-20-supplement-3.

Author information

Authors and Affiliations

Authors

Contributions

JR, HK, and IJ conceived the study. JR led the data analysis with assistance from HK and DY. JR prepared the manuscript with assistance from AJL and IJ. All authors read and commented on the manuscript. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Inkyung Jung.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Additional file

Additional file 1:

Figure S1. Ubiquitous activity of super-enhancers defined in one cell/tissue type. 196 super-enhancer domains are called for only one cell/tissue type but showed constitutive H3K27ac signals across cell/tissue types and also located at less than 20% percentile of tau score. The heatmap of background subtracted log2 transformed H3K27ac signals of the domains across cell/tissue types is shown. Figure S2. Hierarchical clustering of samples using H3K27ac signal on super-enhancer domains. Complete-link hierarchical clustering of distance between samples with respect to log2 transformed H3K27ac RPM of SE domains. H1 human embryonic cell line and its derived cell lines are shown in brown and immortalized cell lines are shown in blue. Cancer cell lines are shown in bold text. Figure S3. Functional characterization of super-enhancer domain classes with alternative putative target genes. a, Heatmaps showing log2(RPKM+ 1) value of putative target gene of unique super-enhancer domains (top), non-unique super-enhancer domains (middle), and common super-enhancer domains (bottom). b, Histogram of tau score for putative target gene expression in each class of super-enhancer domains. Low tau score indicates universal expression pattern. Table S1. GO analysis of alternative putative target genes of unique super-enhancer domains. Table S2. GO analysis of alternative putative target genes of common super-enhancer domains. Table S3. Source of H3K27ac and input ChIP-seq reads. (PDF 574 kb)

Rights and permissions

Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Ryu, J., Kim, H., Yang, D. et al. A new class of constitutively active super-enhancers is associated with fast recovery of 3D chromatin loops. BMC Bioinformatics 20 (Suppl 3), 127 (2019). https://doi.org/10.1186/s12859-019-2646-3

Download citation

  • Published:

  • DOI: https://doi.org/10.1186/s12859-019-2646-3

Keywords