Computational prediction of CTCF/cohesin-based intra-TAD loops that insulate chromatin contacts and gene expression in mouse liver

  1. Bryan J Matthews  Is a corresponding author
  2. David J Waxman  Is a corresponding author
  1. Boston University, United States
6 figures and 6 additional files

Figures

Figure 1 with 2 supplements
Features of TAD boundaries and TAD insulator function.

Profiles in A-E represent a normalized aggregate count of peaks or features along the length of all TADs, sub-divided into 100 equally-sized bins per TAD, where bin #1 is the 5’ start of the TAD and …

https://doi.org/10.7554/eLife.34077.003
Figure 1—figure supplement 1
Additional Features of TAD boundaries.

(A) CTCF binding orientations are divergent at TAD boundaries. The top plot indicates bin enrichment relative to the TAD center, as in Figure 1A. The X-axis begins at the midpoint of one TAD, …

https://doi.org/10.7554/eLife.34077.004
Figure 1—figure supplement 2
Additional Features of TAD boundaries, continued.

(A) Intra-CTCF-motif CpGs are highly hypomethylated compared to the genome-wide average and to neighboring CpGs within 10 kb (‘CpG <10 kb Away’). This is most pronounced for TAD and intra-TAD …

https://doi.org/10.7554/eLife.34077.005
Figure 2 with 7 supplements
Predicted intra-TAD loop anchors share many properties of TAD anchors.

(A) Diagram illustrating intra-TAD loop prediction based on CTCF motif orientation and CTCF and cohesin (Rad21) ChIP-seq binding strength data. Iteration was conducted until 20,000 loops were …

https://doi.org/10.7554/eLife.34077.006
Figure 2—figure supplement 1
Comparison of CTCF Features within TADs and Loop Prediction Improvements.

(A) Fewer than 25% of CTCF, CAC, and tissue ubiquitous CTCF sites are within 25 kb of a TAD boundary. While these features are strongly enriched at TAD boundaries (Figure 1A), the vast majority of …

https://doi.org/10.7554/eLife.34077.007
Figure 2—figure supplement 2
Example Screenshots of Predicted intra-TAD Loops with Observable Interactions in CH12-LX (Mouse B-Cells).

(A-C) Many of the intra-TAD loop structures that we predicted for mouse liver can be seen in the high resolution Hi-C data from the mouse B-cell lymphoma cell line CH12-LX (12). TADs are marked in …

https://doi.org/10.7554/eLife.34077.008
Figure 2—figure supplement 3
Subclasses of CTCF binding events in relation to predicted loops.

(A) Summary of CTCF sites, CAC sites, TAD loop anchors, and intra-TAD loop anchors in mouse liver based on lists in Supplementary file 1C. (1) CTCF peaks found in at least two biological replicate …

https://doi.org/10.7554/eLife.34077.009
Figure 2—figure supplement 4
Intra-TAD loop prediction in two other mouse cell types: mESC and NPC.

(A) Comparison of CAC-mediated intra-TAD loop predicted for mouse embryonic stem cells (mESCs) and neural progenitor cells (NPCs) with those predicted in mouse liver. The number of loops present …

https://doi.org/10.7554/eLife.34077.010
Figure 2—figure supplement 5
Example Screenshots for intra-TAD Loops in mESC and NPC (Mouse).

(A-C) Screenshots of intra-TAD loops predicted in mESC and NPC cells are shown below high-resolution Hi-C data for each cell type. These data provide support for both tissue-specific and …

https://doi.org/10.7554/eLife.34077.011
Figure 2—figure supplement 6
Intra-TAD loop predictions in human cell lines GM12878 and K562.

(A) Comparison of CAC-mediated intra-TAD loop predicted in human lymphoblastoid-derived cells (GM12878 cells) and in human chronic myelogenous leukemic cells (K562 cells). The number of loops …

https://doi.org/10.7554/eLife.34077.012
Figure 2—figure supplement 7
Example Screenshots for intra-TAD Loops in K562 and GM12878 (Human).

(A-C) Screenshots of intra-TAD loops predicted in K562 and GM12878 cells are shown below high-resolution Hi-C data for each cell line. The same genomic region is shown on the left and on the right …

https://doi.org/10.7554/eLife.34077.013
Figure 3 with 2 supplements
Intra-TAD loops show directional interactions and insulate chromatin marks.

(A) TADs and intra-TAD loops both show a stronger orientation of interactions downstream of the motif than other CTCF-bound regions. TAD anchors also show higher inward bias than intra-TAD loops (p<0…

https://doi.org/10.7554/eLife.34077.014
Figure 3—figure supplement 1
Direct evidence of insulation from the asymmetric read distributions for virtual 4C viewpoints anchored at intra-TAD loop anchors or adjacent upstream regions.

(A) Shown is the distribution of Hi-C sequence reads for three intra-TADs referred to in the main text, which provides evidence of directional interactions and insulation in mESCs. For all plots, …

https://doi.org/10.7554/eLife.34077.015
Figure 3—figure supplement 2
Additional features of intra-TADs and their insulation.

(A) TAD anchors are more proximal to divergently-oriented upstream loop anchors, which may explain their elevated local insulation of chromatin contacts and the more diffuse bidirectional contact …

https://doi.org/10.7554/eLife.34077.016
Figure 4 with 2 supplements
Categorization of DHS-based regulatory elements in mouse liver.

(A) Classification of set of ~70,000 open chromatin regions (DHS) identified in adult male mouse liver, based on relative intensities for a combination of H3K4me1 and H3K4me3 marks, and CTCF …

https://doi.org/10.7554/eLife.34077.017
Figure 4—figure supplement 1
Characteristics of five classes of DHS in mouse liver.

(A) Schematic for the classification of open chromatin regions (DHS) in mouse liver based on their H3K4me1, H3K4me3, and CTCF ChIP-seq signals. DHS with >4 rpm for either H3K4me3 or H3K4me1 were …

https://doi.org/10.7554/eLife.34077.018
Figure 4—figure supplement 2
Features of super-enhancers and single-TSS intra-TAD loops.

(A) Shown is the overlap of super-enhancers (SE) for all 19 liver H3K27ac ChIP-seq replicates (see Materials and methods). 503 of the super-enhancers were identified in all 19 liver samples, and are …

https://doi.org/10.7554/eLife.34077.019
Impact of intra-TAD loops on gene expression.

(A) Two possible gene targets were assigned for each super-enhancer within an intra-TAD loop, one target gene for which the TSS is within the intra-TAD loop and another target gene for which the TSS …

https://doi.org/10.7554/eLife.34077.020
Figure 6 with 1 supplement
Alb 4C-seq exemplifies intra-TAD insulation and super-enhancer interaction.

(A) The Alb promoter makes multiple directional contacts with the adjacent super-enhancer region in both male and female mouse liver, as determined by 4C-seq with a viewpoint at the Alb promoter. …

https://doi.org/10.7554/eLife.34077.021
Figure 6—figure supplement 1
Alb 4C-seq replicates and cis/trans 4C-seq signal distribution.

(A) UCSC genome browser screenshots for all biological replicates. The proximal view shows a 150 kb window around the Alb TSS, while the distal view shows a 400 kb view around the Alb TSS. For each …

https://doi.org/10.7554/eLife.34077.022

Additional files

Source Code 1

All code used for the prediction of intra-TAD loops and additional custom scripts.

https://doi.org/10.7554/eLife.34077.023
Supplementary file 1

Intra-TAD loop and CTCF coordinates relevant to Figures 13.

https://doi.org/10.7554/eLife.34077.024
Supplementary file 2

Mouse liver DHS and Super-enhancer information relevant to Figure 4.

https://doi.org/10.7554/eLife.34077.025
Supplementary file 3

Genes and GO term enrichments for genes at refined TAD boundaries and in intra-TADs.

https://doi.org/10.7554/eLife.34077.026
Supplementary file 4

Publicly available datasets used and qPCR primers used to validate CTCF/Rad21 ChIPs.

https://doi.org/10.7554/eLife.34077.027
Transparent reporting form
https://doi.org/10.7554/eLife.34077.028

Download links