Fig. 4

Reduced representation bisulfite sequencing (RRBS) results following developmental exposure to mono(2-ethylhexyl) phthalate (MEHP) and 5-azacytidine (5AC). a General statistics of RRBS analyses, showing mutual analyzed Cs in CpG context at a read depth of 10, number of tiles analyzed with regional analysis and the number of tiles overlapping different features. b–d Methylation profile spanning transcriptional start sites (TSSs), CpG islands and zebrafish conserved non-genic elements (zfCNEs). e Methylation levels relative to control overlapping different features. Error bars represent SEM; significance was calculated with two-way ANOVA, using Bonferroni multiple comparisons correction (P < 0.05). f Volcano plot of the methylation difference between MEHP and control. g Number of MEHP hypo and hyper differentially methylated regions (DMRs). h Enrichment plot of all tiles and DMRs overlapping different features. Significant enrichment calculated using a hypergeometric test. i Volcano plot of the methylation difference between 5AC and control. j Number of 5AC hypo and hyper DMRs. k Enrichment plot of all probes and DMRs overlapping different features. Significant enrichment calculated using a hypergeometric test