Violin plots of DNA methylation in promoter (Tetraodon, n = 12,896; Xenopus, n = 12,704; human liver, n = 22,680; human lung, n = 23,012), first exon (Tetraodon, n = 11,887; Xenopus, n = 10,361; human liver, n = 20,383; human lung, n = 20,704), first intron (Tetraodon, n = 11,420; Xenopus, n = 12,202; human liver, n = 20,029; human lung, n = 20,757), rest of exons (Tetraodon, n = 12,618; Xenopus, n = 12,662; human liver, n = 18,961; human lung, n = 19,331) and rest of introns (Tetraodon, n = 11,840; Xenopus, n = 11,905; human liver, n = 16,930; human lung, n = 17,007) divided into deciles based on increasing ranking of gene expression. Box plots with rotated kernel density plots at both sides indicate the interquartile range, and white central dots the median of the distribution. Correlations between DNA methylation and gene expression were measured using Spearman’s rank correlation coefficient (ρ), and the significance levels are reported as follows: ***p < 0.001