Fig. 4

Evolutionary conservation of the Dynamic GbM state. A Representative Dynamic GbM gene homologs from six different species. Each red bar represents a single CG, and the bar height is the percent methylation heterogeneity across cells. Dashed lines represent the cutoffs used to call heterogeneously methylated cytosines (10–85% of cells). Rightmost panel: the estimated divergence time between species and Arabidopsis thaliana, obtained from TimeTree [42]. Points represent the median divergence time estimate, and vertical dashed lines represent the 95% confidence interval. B Violin plots showing the number (left panel) and proportion (right panel) of heterogeneously methylated CGs across all homologs of Dynamic GbM genes, and ten identically sized groups of randomly selected genes. Solid lines represent the median and dashed lines represent the interquartile range. *p =  < 0.05, ***p =  < 0.0005 (two-tailed t-test)