10. Figure legends

Figure 1: Hepatic DNA methylation landscape in K. marmoratus . A: Frequency distribution of MSPI fragment size. Mean fragment size was 95 bp and median 88 bp; B: Overall distribution of analysed fragments in genomic elements; C: Frequency distribution of DNA methylation % for MSPI fragments. Distribution of MSPI fragments follows a bimodal distribution, which is expected given the binary character of DNA methylation in cells; D: Boxplot of DNA methylation (%) of fragments according to their location in genomic elements (including medians and quantile 25-75%).

Figure 2: Distribution of lowly methylated (<20% methylation) (A) and highly methylated (> 80% methylation) (B) fragments among the genomic regions (gene body, promoter, upstream, intergenic) of the liver reduced representative genome of the mangrove rivulus.

Figure 3: A: Volcano plot representations of significance and differential methylation in analysed fragments. Dotted lines represent thresholds of significance: -10LogP20, corresponding to P0.01 and mean methylation difference 10% A: 4 ng/L vs control; B:120 ng/L vs control; C: 4 ng/L vs 120 ng/L. B: Overall distribution of DMFs in genomic elements; C: Venn Diagram showing specific and common DMFs between group comparisons; D: Methylation levels (%) of common fragments between 4 vs Ctl and 120 vs Ctl comparisons (Da,Db), between 4 vs Ctl and 4 vs 120 comparisons (c,d,e) and between 120 vs Ctl and 4 vs 120 ng/L comparisons (f). For clearer representation of methylation differences, y axes start at 40 %.

Figure 4: Network of genes corresponding to significant DMFs (located in promoter region or in gene body) at 4 ng/L compared to control. Blue: hypermethylation at 4 ng/L compared to control, Yellow: hypomethylation at 4 ng/L compared to control.

11. Figures