Fig. 6: Loci that are susceptible to spontaneous epiallele formation have a greater contact frequency with H3K9me2 region.

a Hi-C data was used for identifying significant three-dimensional contacts between genomic regions. The schematic plot depicts the interaction from regions possessing H3K9me2 and genes. b Enrichment test for genes that have contacts with H3K9me2. Yellow bars show the proportion of each type of genes that contacts with H3K9me2. Grey bars show the proportion of each type genes in comparison with all coding genes. Fisher’s exact test was used for enrichment test. p-value is based on one-sided test with alternative hypothesis that odds ratio is greater than 1. c Distribution of distance between three types of genes (gbM, UM, teM) and H3K9me2 regions. Two-sided Wilcoxon tests was used for pairwise comparison of each two groups without adjustment for multiple comparisons. d Distribution of the number of contacts between three types of genes (gbM, UM, teM) and H3K9me2 regions. gbM genes were further classified into a group with ectopic mCHG (gbM; mCHG+) and a group without ectopic mCHG (gbM; mCHG-) based on ddm1 mutant line data. Two-sided Wilcoxon tests was used for pairwise comparison of each two groups without adjustment for multiple comparisons. Box plots show a median center line, the lower and upper hinges are the first and third quartiles. Whiskers represent 1.5x the interquartile range. Source data underlying Fig. 6b–d are provided as a Source Data file.