<p>You are right that a single molecule in a single position is either methylated or not methylated. However:</p>
<p>First, assuming your organism of interest is diploid (or of higher ploidy) one of the chromosomes could be methylated, the other not. That would give you a level of 0.5 and can be found in imprinted regions (where the paternally inherited chromosome is methylated and maternal not, or vice versa).</p>
<p>Second, except if you use a pure cell population of a single cell type it is likely that there will be cellular heterogeneity in your sample, e.g. from different types of blood cells or neurons vs glial cells in the brain. Some of these may have methylated CpGs in this locus (depending on cell-type-specific gene activity) while others can be unmethylated. Depending on the relative proportions of these cell types the methylation level can be anything between 0 and 1 for a mixture of cell types.</p>
<p>Thirdly, sometimes CpGs get grouped together. For example, 5 nearby CpGs are combined to give a single methylation level for this locus, and then the level can be 0; 0.2; 0.4; 0.6; 0.8 and 1.</p>
<p>Lastly, there could be technical artifacts or technical reasons why a position was not called as methylated (insufficient quality for a single read?), leading to deviation from 0 or 1.</p>

CC BY-SA 4.0
2020-08-05T13:21:36.370Z
Thank you very much! I noticed some databases supply a single methylation level for each promoter region, which I could not completely understand before. Now I think that's the third case in your answer. Your answer is really helpful.
Jie
CC BY-SA 4.0
2020-09-16T02:39:25.893Z
Thanks for the explaination. One further question though, why would the methylation distribution are biomodal when 1) methylation could be heterozygous and 2) we have mixture of cell type where each type represent different methylation pattern?
unicorn
CC BY-SA 4.0
2020-08-05T10:12:40.117Z
2020-08-05T10:12:40.117Z
Wouter De Coster
<p>For Bisulfite Sequencing data, we can calculate <strong>methylation level</strong> for a specific locus,
<span class="math-container">$$
\mbox{methylation level} = \frac{\mbox{# methelyated reads}}{\mbox{# total reads}}
$$</span>
Here I think, for this locus, it can only have one status, methelyated or not. Or to say, <strong>methylation level</strong> here can only be 0 or 1. Why does it have other values in (0,1)?</p>

<phylogenetics><sequence-analysis><methylation>
How to understand methylation level?
# How to understand methylation level?

(Asked by: username_0 on 2020-08-05)

For Bisulfite Sequencing data, we can calculate **methylation level** for a specific locus,
$$
\mbox{methylation level} = \frac{\mbox{# methelyated reads}}{\mbox{# total reads}}
$$
Here I think, for this locus, it can only have one status, methelyated or not. Or to say, **methylation level** here can only be 0 or 1. Why does it have other values in (0,1)?

**Tags:** phylogenetics, sequence-analysis, methylation

---

## 1 answers

### Answer by username_1 on 2020-08-05. Score: 6 | Accepted answer

You are right that a single molecule in a single position is either methylated or not methylated. However:

First, assuming your organism of interest is diploid (or of higher ploidy) one of the chromosomes could be methylated, the other not. That would give you a level of 0.5 and can be found in imprinted regions (where the paternally inherited chromosome is methylated and maternal not, or vice versa).

Second, except if you use a pure cell population of a single cell type it is likely that there will be cellular heterogeneity in your sample, e.g. from different types of blood cells or neurons vs glial cells in the brain. Some of these may have methylated CpGs in this locus (depending on cell-type-specific gene activity) while others can be unmethylated. Depending on the relative proportions of these cell types the methylation level can be anything between 0 and 1 for a mixture of cell types.

Thirdly, sometimes CpGs get grouped together. For example, 5 nearby CpGs are combined to give a single methylation level for this locus, and then the level can be 0; 0.2; 0.4; 0.6; 0.8 and 1.

Lastly, there could be technical artifacts or technical reasons why a position was not called as methylated (insufficient quality for a single read?), leading to deviation from 0 or 1.

> **Comments:**
> **username_0**: Thank you very much! I noticed some databases supply a single methylation level for each promoter region, which I could not completely understand before. Now I think that's the third case in your answer. Your answer is really helpful.  
> **username_2**: Thanks for the explaination. One further question though, why would the methylation distribution are biomodal when 1) methylation could be heterozygous and 2) we have mixture of cell type where each type represent different methylation pattern?  

---

