Search

Additional file 4: Figure S4. TAD structure and CTCF binding is preserved during the 24 hours. A. TADs overlap between time points. B. TADs size distribution at all time points. C. 50kb resolution Median Observed/Expected Hi-C signal around 1000 randomly selected TADs from ZT6, 12 and 18 plotted on the indicated time points. TADs were scaled to fit the five central bins. D. CTCF ChIP-seq motif analysis and peak overlap between ZT0 and 12. Genomic tracks of CTCF ChIP-seq signal at example regions harbouring circadian genes. E. Above, the same metaplots as in C but for 1000 randomly CTCF peaks found at ZT12 plotted using ZT12 and ZT0 Hi-C contacts. Below, metaplot using CTCF peaks from mESC and plotted using liver ZT0 and 12 Hi-C contacts. CTCF peaks are at the central bin of the metaplot. F. TAD and cTAD size comparison (p-value < 2.2e-16, Wilcoxon rank sum test).

Additional file 7: Figure S7. Core clock gene promoter vs output circadian genes interaction landscapes.A-D. Virtual 4C for Rorc, Npas2, Nr1d2 and Per2, core clock circadian gene promoters from P-CHi-C data at all time points during the day. The Rorc circadian gene promoter contacts Cgn circadian gene promoter and both peak in transcription at ZT18. E-F Virtual 4C for Dhrs3 and Ppp1r3c output circadian genes in the liver. Acrophases are written next to the gene names. Genomic tracks show significant contacts as arcs and chromatin features including liver H3K4me3, H3K4me1, H3K27ac, DNAseI, eRNAs and TADs. Expression profiles for all genes can be found in Figure S3. The interaction profiles of core clock genes display less contacts that dynamically change over time. The two output gene contact profiles show more saturated contacts that are mostly constant during the 24 hours.

Additional file 3: Figure S3. Expression profiles of the circadian genes displayed along the paper. RNA-seq signal tracks for circadian genes shown across the paper for the four timepoints during the circadian cycle. Arntl, Npas2, Cgn, Nr1d1, Dhrs3, Nr1d2, Per2, Ppp1r3c, Tef, Rnf125, Rorc and Aco2 circadian gene expression profiles are presented.

Additional file 6: Figure S6. Circadian gene promoter-promoter interactions. A. Promoter-promoter contact network at ZT0 with the circadian genes marked in blue. B. Above, number of edges (interactions) between circadian gene promoters at ZT0, 6, 12 and 18 hours of the day (blue line) compared to a random set of non circadian promoters. Below, z-score for circadian gene promoter contacts compared to the random sampling. C. Quantification of the read counts supporting all circadian gene promoter-promoter interactions compared to a random set of non circadian promoter-promoter interactions (*** p-value < 0.001, Mann Whitney test) D. Transcriptional phase distributions of circadian promoters in contact with diurnal or nocturnal circadian promoters (all p-values < 0.0001, Wilcoxon signed rank test) for our circadian intronic gene set. E. The same as in D for the circadian gene set detected through GRO-seq (Fang et al., 2014) (all p-values < 0.0001, Wilcoxon signed rank test) F. Virtual 4C landscape for the Tef circadian gene promoter from P-CHi-C data at all time points during the day. The acrophase of Tef is written next to the gene name. Expression profiles for both genes can be found in Figure S3. Genomic tracks show significant contacts as arcs and chromatin features including liver H3K4me3, H4K4me1, H3K27ac, DNAseI, eRNAs and TADs. Tef gene promoter contacts Aco2 gene promoter both with acrophase at ZT12.

Additional file 1: Figure S1. OCCs present fluctuations in histone modifications. A. Heatmap of PC1 values for stable compartments from the three independent Hi-C biological replicates and B, the merged Hi-C biological replicates. C. PC1 eigenvector correlations between every pair of timepoints (Kruskal-Wallis rank sum test p-value < 2e-16). D. Proportion of OCCs and constant compartments in the mouse genome. E. Barplot of the frequency of the different OCCs categories with compartments switching at different times of the 24 hours. For example AABA refers to OCCs with compartment assignment A at ZT0, A at ZT6, B at ZT12 and A at ZT18. F. Chromatin features (H3Kme4 and H3Kme1) at OCCs across time points (AABB, ABBB and ABAB OCCs are shown) (*** all p-values < 0.001 except for the H3K4me1 AABB, one way ANOVA and Tukey post hoc test).

Additional file 9. Review history.

Additional file 2: Figure S2. Total RNA-seq analysis at four timepoints during the circadian cycle. A. Spearman correlation analysis of the four biological replicates from ZT0 and ZT12 (r ≥ 0.85). B. Heat map of the relative transcription (Z scores) of 1,257 oscillating genes sorted by oscillation phase. C. Individual expression profiles and genome browser tracks showing examples of the RNA-seq signal for Arntl, Per1, Nr1d1, Cry1, Ppp1r3c and Gsk3a oscillating genes (Fragments Per Kilobase of transcript per Million mapped reads, FPKM) (n=4, q-value