Why Epigenomic Tables Need Careful Reading
Epigenomic tables often place several kinds of measurements side by side, but their entries do not all mean the same thing. A DNA methylation entry may estimate the fraction of observed molecules that are methylated at a locus. A chromatin-accessibility entry may be a fragment count or a processed signal for a genomic region. The first is bounded between zero and one; the second depends on the assay and processing and may have no such bound.
That distinction matters before comparing samples or relating the measurements to other molecular layers. A methylation value of 0.7 has a proportion-based interpretation, subject to how it was measured. An accessibility value of 0.7 usually does not mean that 70% of a region is accessible. Its meaning depends on the signal’s definition and units.
The earlier material on RNA sequencing tables introduced sample-by-feature measurements for sequencing assays. Here the focus is narrower: how to read methylation fractions and accessibility measurements without treating them as interchangeable quantities.
Methylation Fractions at a Locus
For example, if 42 of 60 informative molecules are classified as methylated, the observed fraction is \(42/60=0.70\). The value describes the observed molecules at that locus in that sample; it does not establish that every cell has the same methylation state. In a bulk tissue sample, the fraction can reflect a mixture of cell types, alleles, and molecular states.
Coverage is part of the interpretation. Two fractions can be numerically equal but based on very different amounts of information. A fraction of \(8/10=0.80\) is based on ten observations, whereas \(800/1000=0.80\) is based on one thousand. Neither value alone says how representative the sampled molecules are of the underlying biological material.
A common simple model makes this dependence explicit. Suppose, as an idealization, that the \(c_{ij}\) observations are independent Bernoulli trials with the same probability \(p_{ij}\) of being classified as methylated. Then \(m_{ij}\) has a binomial distribution with size \(c_{ij}\) and probability \(p_{ij}\), and \(\hat{\beta}_{ij}\) estimates \(p_{ij}\). Under those assumptions, its standard error is estimated by \(\sqrt{\hat{\beta}_{ij}(1-\hat{\beta}_{ij})/c_{ij}}\). This is a model-based summary of sampling variability, not a guarantee that read observations are truly independent.
Worked Example: Calculating a Methylation Fraction
Suppose a hypothetical sample has 52 methylated observations among 80 informative observations at a particular CpG locus. The total coverage is 80, so the observed methylation fraction is
The estimated fraction is 0.65, or 65% of the informative observations at this locus. Under the simple independent-binomial model, an estimated standard error is
Thus, the model-based standard error is about 0.053, on the fraction scale. It describes approximate read-sampling variability under the stated model. It does not capture biological variation among independently sampled individuals, and it may understate uncertainty if observations are dependent or the assay has additional error.
Chromatin-Accessibility Signals
For example, ATAC-seq records DNA fragments associated with transposase access, while DNase-based assays record evidence related to nuclease cleavage. A table may report counts assigned to a region, or a processed signal derived from those counts. Raw counts are nonnegative integers; after normalization or transformation, values may be decimal and need not be integers. Such processing does not turn the measurement into a proportion.
An accessibility signal is evidence generated by a particular assay and analysis pipeline, not a direct, universal measurement of how “open” a region is. A higher signal can indicate more observed assay evidence under compatible measurement conditions. It does not, by itself, give the percentage of cells with accessible chromatin, the probability that a base is accessible, or a causal effect on gene activity.
Worked Example: Why Raw Accessibility Counts Can Mislead
Consider a hypothetical region measured in two samples using the same assay. Sample A has 240 assigned fragments and 12 million total mapped fragments. Sample B has 360 assigned fragments and 30 million total mapped fragments. To illustrate a simple library-size adjustment, calculate fragments per million mapped fragments:
The raw count is larger in Sample B, but the adjusted signal is larger in Sample A: 20 rather than 12 fragments per million mapped fragments. The raw-count comparison alone is therefore affected by the different total numbers of mapped fragments. This adjustment is only a simple illustration; a real analysis must use a suitable, consistently applied processing method and consider whether the samples are otherwise comparable.
The adjusted values are signals in a stated unit, not percentages. It would be incorrect to say that the region is “20% accessible” in Sample A.
Reading the Table Together
A sample-by-locus methylation table has one entry for each sample–locus combination, and a corresponding accessibility table may summarize regions rather than individual methylation sites. The locus definitions matter: a single CpG site and a wider genomic region are not necessarily the same feature. A table can therefore place values from related genomic locations together without making them identical measurements.
Methylation and accessibility may be biologically related in some contexts, but their relationship is not a universal rule. A high methylation fraction at one site does not automatically imply low accessibility across a nearby region. The relationship can depend on genomic context, cell mixture, assay design, and how regions are defined. Neither measurement alone establishes a change in gene expression or a causal mechanism.
Worked Example: Interpreting Two Epigenomic Measurements
Suppose a hypothetical table reports measurements for the same two samples at a promoter-associated region. Methylation is measured at a specified CpG site within that region. Accessibility is reported as normalized fragments per million mapped fragments using the same pipeline for both samples.
| Sample | Methylation fraction at CpG | Accessibility signal (fragments per million) |
|---|---|---|
| Sample 1 | 0.82 | 14 |
| Sample 2 | 0.24 | 48 |
At the measured CpG, Sample 1 has an observed methylation fraction of 0.82, whereas Sample 2 has 0.24. Under the stated processing, Sample 2 also has the larger accessibility signal: 48 compared with 14 fragments per million. The arithmetic difference in accessibility signal is \(48-14=34\) fragments per million.
A careful description is that these two samples show different measured methylation fractions at the specified CpG and different normalized accessibility signals for the specified region. The table alone does not show that methylation caused the accessibility difference, that the measurements represent every cell, or that the pattern applies to other loci or samples.
Assumptions, Diagnostics, and Failure Modes
Before interpreting values, check how each measurement was produced. The binomial standard-error calculation for methylation assumes independent observations with a common methylation probability at the locus. In practice, dependence among reads, amplification, mapping uncertainty, conversion or classification error, and biological mixtures can undermine that idealization. The fraction remains a useful summary of observed calls, but its uncertainty should not be treated as fully described by the simple formula when those complications are substantial.
For accessibility comparisons, a meaningful numerical contrast requires compatible feature definitions, assay units, and processing. Raw counts from samples with different sequencing depths are not directly comparable as though the total measurement effort were equal. Even processed values can be misleading if pipelines or region boundaries differ, or if technical effects are aligned with the sample groups being compared.
- Check methylation coverage. Confirm that the denominator is positive and inspect coverage alongside the fraction. Low coverage means the estimate rests on few observations.
- Check the range and encoding. A sequencing-based fraction should be between 0 and 1. A value outside that range suggests a different encoding, a processing issue, or a data error; do not silently interpret it as a fraction.
- Check signal units. Determine whether accessibility values are raw counts, normalized counts, or transformed scores. The same number can mean different things under different definitions.
- Check comparability. Confirm that samples use compatible assays, processing, and region definitions before interpreting between-sample differences.
- Check the biological unit. Read-level observations are not the same as independent biological samples. Many reads from one sample do not substitute for measurements on multiple individuals.
An absent accessibility peak or a zero count also needs context: it may mean no fragments were observed under the assay and processing, not that a region is physically inaccessible in every cell. Similarly, a methylation fraction of zero means no methylated observations were recorded among the informative observations used for that entry; it does not prove that methylation is impossible in the underlying population.
Common Mistakes and Alternatives
| Incorrect reasoning | Correction |
|---|---|
| “A methylation value of 0.7 means 70% of cells are methylated.” | It is the fraction of informative observations classified as methylated. A cell-level interpretation needs additional assumptions and evidence. |
| “An accessibility score of 0.7 is a 70% probability of open chromatin.” | Accessibility signals are assay- and pipeline-specific measurements unless explicitly defined as probabilities. |
| “The sample with more accessibility fragments must have higher accessibility.” | Raw counts may differ because of sequencing depth or other technical factors. Compare signals with compatible processing and units. |
| “Higher methylation must mean lower accessibility.” | The relationship is context-dependent. Report the measurements separately unless an appropriate analysis supports a relationship. |
Different summaries answer different questions. A methylation fraction uses methylated observations divided by total informative observations; it preserves the proportion-based interpretation. A methylated count alone gives the numerator but omits the denominator. Accessibility may be represented by raw counts when describing observed fragments, or by a suitable normalized signal when comparing compatible samples. Those are alternatives for summarizing assay evidence, not interchangeable scales. The right choice depends on the measurement process and the intended comparison.
The central limitation is that each entry is a summary of assay evidence at a defined feature. Methylation fractions are bounded but can be imprecise at low coverage; accessibility signals can be useful for comparisons but lack a universal unit. Neither table alone resolves cell-type mixture, technical bias, or causal direction. Interpret values alongside their denominators, units, feature definitions, and study design.
Check Your Understanding
Answer each question, then compare your reasoning with the solutions.
- A locus has 18 methylated observations among 30 informative observations. What is its observed methylation fraction?
- Two samples have methylation fractions of 0.50, based on coverage of 12 and 600, respectively. Are the estimates based on the same amount of information? Explain.
- One sample has 100 accessibility fragments from 5 million mapped fragments; another has 150 from 15 million. Using fragments per million, which has the larger adjusted signal?
- Why does a methylation fraction of 0.90 at a CpG and a low accessibility signal nearby not, by itself, establish that methylation caused reduced accessibility?
Solutions. (1) \(18/30=0.60\), or 60% of the informative observations. (2) No. Both observed fractions are 0.50, but the first is based on 12 observations and the second on 600; coverage affects the precision of the observed fraction. (3) The first is \(100/5=20\) fragments per million; the second is \(150/15=10\) fragments per million. The first has the larger adjusted signal. (4) The measurements are observational summaries at specified features. Their association does not establish causal direction, and the values can also reflect cell mixtures, assay properties, or other factors.