Fig. 7: Analysis of ChIP-seq data revealed that H3R17me2 is positively correlated with H3K27me3.

a, b Heatmaps showing the Log₂FC of H3R17me2 and H3K27me3 across the gene body, based on normalized ChIP-seq read coverage. Twenty-five fly heads were collected at ZT8 and used for ChIP-seq. c, d Boxplot showing the Log₂FC value of H3R17me2 and H3K27me3 across the gene body, respectively, Corresponding to (a, b). n = 13967 genes. The middle line within the box plots corresponds to the median; the lower and upper hinges correspond to the first and third quartiles; the upper whisker extends from the hinge to the largest value no further than 1.5 times the interquartile range; and the lower whisker extends from the hinge to the smallest value at most 1.5 times the interquartile range. P-values from unpaired two-tailed Student’s t-test are indicated, where ****P < 0.0001. e, f Heatmap showing H3R17me2 and H3K27me3 enrichment at dilp2, dilp3 and dilp5 locus based on normalized ChIP-seq read coverage. g Heatmap showing read counts by RNA-seq at dilp2, dilp3 and dilp5 locus. h–k Genome-wide correlation plots of H3R17me2 and H3K27me3 occupancies at 3-day-old w¹¹¹⁸, 3-day-old LKRSDH^MB01942, 40-day-old w¹¹¹⁸ and 40-day-old LKRSDH^MB01942. The correlation coefficient (R) is indicated. l, m Genome-wide correlation plots of H3R17me2 and H3K27me3 occupancies at 3-day-old LKRSDH^MB01942 / 3-day-old w¹¹¹⁸ and 3-day-old w¹¹¹⁸ / 40-day-old w¹¹¹⁸, respectively. n Genome-wide correlation plots of RNA-seq signal at 3-day-old LKRSDH^MB01942 / 3-day-old w¹¹¹⁸ and 3-day-old w¹¹¹⁸ / 40-day-old w¹¹¹⁸. o Heatmaps showing the dynamics of H3K27me3, H3R17me2 and H3K27me3/H3R17me2 co-modified genes in 3-day-old w¹¹¹⁸, 3-day-old LKRSDH^MB01942, 40-day-old w¹¹¹⁸ and 40-day-old LKRSDH^MB01942 (left panel). The ratios of H3K27me3, H3R17me2 and H3K27me3/H3R17me2 co-modified genes are also shown (middle panel). The red line indicates the fraction of domains that is marked only by H3K27me3, the green line represents the fraction of domains that is marked only by H3R17me2, and the blue line indicates the fraction of H3K27me3/H3R17me2 co-modified genes. The top five GO enrichment of different gene clusters is shown on the right panel. Source data are provided as a Source Data file.