=5 reads on both strands) for both E14-d0 and
105
E14-d6. Among these cytosines, 7,583,856 are located in CpG context, 35,228,969 in CHG
106
context, and 115,745,872 in CHH context. They represent 17.3%, 16.0% and 13.8% of the total
107
number of CpG, CHG and CHH sites in the mouse genome, respectively.
108
We first analyzed the methylation levels (ML) of CpG and non-CpG sites in E14-d0 and
109
E14-d6. Similar to previous reports (Ramsahoye et al. 2000; Lister et al. 2009; Xie et al. 2012), 5 / 30
Downloaded from genome.cshlp.org on July 14, 2014 - Published by Cold Spring Harbor Laboratory Press
110
we detected abundant DNA methylation calls (8.0%) in non-CpG context in E14-d0 and the
111
number decreased to 4.3% in E14-d6 (Figure S3A&B). Most cytosines in CpG context are
112
highly methylated in both E14-d0 and E14-d6, and their methylation levels are characterized by
113
a bimodal distribution (Figure S3C&D). The average methylation level of CpG dyads is 86.3%
114
in E14-d0, and increases to 91.5% in E14-d6 (Figure S3C&D). In contrast, most cytosines in
115
non-CpG context are either unmethylated or lowly methylated (Figure S3E-H).
116
We further determined the methylation fidelity (MF), as defined in the method section,
117
which represents the percentage of symmetrically methylated or unmethylated CpG dyads for a
118
given position. A symmetrical methylation status of CpG dyad (either methylated or
119
unmethylated) indicates successful methylation inheritance, while an asymmetrical methylation
120
status (hemi-methylated CpG dyad) indicates gain or loss of the methylation pattern. Most CpG
121
dyads have high methylation fidelity and the average methylation fidelity is as high as 88.5% and
122
91.9% for E14-d0 and E14-d6, respectively (Figure S4). While most CpG dyads maintain
123
similar methylation levels and fidelities during mES cell differentiation, our results indicate that
124
both the methylation level and fidelity increase for a considerable fraction of CpG dyads in E14-
125
d6 (Figure S5). More specifically, 24.7% of CpG dyads (17.0% increased, and 7.7% decreased)
126
have methylation level changes over 20% (Figure S5C), and 27.5% of CpG dyads (17.3%
127
increased, and 8.3% decreased) have methylation fidelity changes over 20% (Figure S5D).
128
Methylation level and fidelity vary at distinct genomic regions
129
Numerous studies indicate that DNA methylation levels tend to vary across distinct
130
genomic regions (Doi et al. 2009; Lister et al. 2009; Laurent et al. 2010). Consistently, we
131
observed low methylation levels in promoter regions, and the methylation levels are highly
132
correlated with the distance to transcription starting sites (Figure 2). Low methylation levels
133
were also observed for CGIs and CGI shores (p