1 Basics of MNase-seq

1.1 What MNase-seq measures

MNase-seq (Micrococcal Nuclease sequencing) is a chromatin-mapping method that profiles where DNA is protected by nucleosomes and, to a limited extent, where DNA is relatively exposed. Cells or nuclei are treated with micrococcal nuclease (MNase), an enzyme that preferentially cuts accessible DNA. Because nucleosome-bound DNA is more resistant, digestion enriches for short fragments that correspond largely to mono-nucleosomal DNA. After sequencing and mapping, the resulting fragment patterns are used to infer nucleosome positions (or nucleosome “centers” or dyads), nucleosome occupancy (how often nucleosomes are present), and chromatin-scale structural features reflected in fragment spacing.

The method is especially informative for nucleosome-scale organization, where it can detect preferred placement along the genome and compare chromatin states across conditions.

1.2 MNase enzyme properties and digestion behavior

MNase is a nonspecific endo-exonuclease that cleaves DNA in a pattern influenced by local DNA accessibility. In chromatin, DNA that is not tightly wrapped around histone octamers is typically digested more quickly than nucleosome-protected segments. The digestion endpoint therefore strongly affects the fragment distribution: partial digestion yields a mixture of subnucleosomal and nucleosomal fragments, while extensive digestion can further chew nucleosome-protected DNA and reduce recoverable mononucleosome-sized material.

Because MNase activity depends on enzyme concentration, incubation time, temperature, buffer composition, and how the chromatin is prepared, digestion is usually tuned so that mononucleosome-sized fragments are enriched without excessive overdigestion.

1.3 Mononucleosome fragment selection

After MNase digestion, the DNA mixture is processed to isolate fragments centered around the mono-nucleosome size range. These fragments are expected to include DNA wrapped around one nucleosome plus additional flanking DNA that MNase leaves intact at the chosen digestion level. In practice, “mononucleosome” fragments are selected by size filtering (often around a characteristic length window), producing a sequencing library that is biased toward nucleosome-associated DNA.

This selection step is a major determinant of what the sequencing signal reflects: the chosen size window shapes the inferred nucleosome positions and affects how much contamination from di- or subnucleosomal fragments contributes to downstream analyses.

1.4 Fragment size and nucleosome periodicity intuition

When nucleosomes occupy chromatin regularly, MNase can produce a repeating pattern of fragment sizes and DNA cleavage products along the genome. The nucleosome repeat length—an effective spacing between adjacent nucleosomes—can manifest as periodicity in coverage and fragment characteristics. Even when nucleosome arrays are not perfectly regular, the fragment distribution often contains hints of ordered chromatin packaging.

In sequencing read space, periodic signals can appear in aggregate profiles or in genome-wide analyses, serving as an intuitive readout for nucleosome organization beyond single-position calling.

2 Experimental workflow

2.1 Sample preparation

2.1.1 Cell vs. nuclei MNase treatment

MNase-seq can be performed starting from whole cells or isolated nuclei. Nuclei preparations often provide more consistent access to chromatin while limiting unwanted cytoplasmic components. Whole-cell treatment can be simpler operationally, but differences in permeabilization and enzyme access may change digestion kinetics and the resulting fragment profiles. The choice can therefore affect comparability across samples, and it can influence the calibration needed for digestion optimization.

MNase-seq is commonly done with either native chromatin (no crosslinking) or with mild crosslinking to preserve chromatin states. Crosslinking can stabilize nucleosome occupancy and reduce shifts caused by processing, but it can also complicate subsequent steps because crosslinks must be reversed and the DNA-protein chemistry can alter accessibility to MNase. Without crosslinking, digestion occurs on more dynamic chromatin and may be more sensitive to handling conditions.

Whether crosslinking is used typically depends on the biological question and on practical considerations about how stable the relevant chromatin features are expected to be during sample processing.

2.2 MNase digestion optimization

2.2.1 Titration to achieve appropriate digestion levels

Optimization usually centers on finding a digestion endpoint that enriches mono-nucleosome-sized DNA. Titration can involve varying MNase concentration and monitoring the DNA fragment distribution by gel electrophoresis. A useful target is a clear mononucleosome band with limited remaining higher-molecular-weight material and minimal degradation into very small fragments. Because different cell types and chromatin states digest at different rates, titration is typically performed per sample type and sometimes per batch.

If digestion is too mild, the signal may include more di- and oligo-nucleosomal fragments, complicating interpretation of nucleosome positions. If digestion is too aggressive, mononucleosome fragments can be lost to further cutting.

2.2.2 Time, temperature, and enzyme handling

Digestion kinetics depend on incubation time and temperature, with MNase activity increasing under conditions that support enzyme function. Enzyme handling is also important: MNase stock quality, freeze-thaw cycles, and mixing consistency can contribute to variability. Standardized protocols often specify how quickly samples are processed after adding MNase, how digestion is quenched, and how conditions are balanced across replicates.

Even small deviations in handling can shift fragment-size distributions and alter the balance between nucleosome-protected and accessible DNA.

2.3 DNA purification and size selection

2.3.1 Gel-based vs. bead/column-based approaches

After digestion, DNA is purified and then size selected to enrich mononucleosome fragments. Gel-based methods (for example, excising the appropriate band) provide direct control over fragment length but can be labor-intensive and may reduce yield. Bead- or column-based size selection can be easier to scale and automate, though it requires calibration to match the intended fragment window.

The chosen approach affects both recovery efficiency and the exact fragment-length distribution that enters sequencing.

2.3.2 Recovering mononucleosome-sized DNA

Recovering mononucleosome-sized DNA requires balancing specificity with yield. If the selected window is too narrow, fewer reads may be obtained and stochastic variation can increase. If the window is too wide, fragments derived from di-nucleosome or subnucleosomal DNA may contaminate the mononucleosome signal. Many pipelines therefore include a downstream filtration step that further refines fragment selection using the observed length and mapping characteristics.

2.4 Library preparation and sequencing

2.4.1 End repair, adapters, and PCR enrichment

Library construction converts selected DNA fragments into sequencing-ready molecules by end-repair, adapter ligation, and PCR amplification. Enzymatic steps can introduce artifacts such as bias toward particular sequence contexts or fragment ends, and PCR can preferentially amplify certain molecules. Some library designs incorporate strategies to mitigate duplication effects or to control amplification cycles, but all MNase-seq experiments require attention to amplification consistency.

Because MNase-seq already involves enrichment via digestion and size selection, careful control of PCR conditions helps ensure that read counts reflect biology rather than overamplification artifacts.

2.4.2 Read length and sequencing depth guidance

Read length should be sufficient to map fragments uniquely for the target genome and to support nucleosome positioning calling. Sequencing depth determines how well nucleosome occupancy patterns and differential comparisons can be estimated, particularly in low-signal regions. Depth requirements vary with genome size, expected nucleosome density, and the planned level of statistical comparison, such as whether the goal is robust differential occupancy across many genomic features.

A practical aspect is matching depth across replicates and conditions to avoid conflating biological differences with sampling noise.

2.5 Controls and experimental replication

2.5.1 Technical vs. biological replicates

Technical replicates measure variance introduced by library preparation, sequencing, or mapping, while biological replicates capture variation across independently prepared samples. For chromatin mapping, biological replication is usually essential because nucleosome patterns can vary due to growth conditions, cell-cycle effects, and sample handling. Technical replication can be useful for troubleshooting but cannot substitute for biological coverage.

Replication strategy should also consider how digestion optimization is performed, since tuning parameters may differ across biological batches.

2.5.2 Negative/blank controls and spike-ins (if used)

Blank controls (for example, no-enzyme or no-template controls) help identify contamination and procedural background. Spike-in controls can be used to facilitate cross-sample normalization, particularly when comparing experiments with different efficiencies or sequencing yields. When used, spike-ins provide an external reference that is independent of the endogenous genome.

Whether spike-ins are adopted depends on experimental goals and budgeting, but their inclusion can improve interpretability in large comparative studies.

3 Data processing and quality control

3.1 Read preprocessing

3.1.1 Adapter trimming and quality filtering

Sequencing reads must be cleaned by removing adapter sequences and low-quality bases. Trimming decisions affect alignment rates and the accuracy of fragment boundary inference, which in turn influences nucleosome calling. Quality thresholds are often set to retain high-confidence bases while reducing the impact of sequencing artifacts.

Preprocessing also includes discarding reads that fall below a minimum length, since short fragments can map ambiguously or distort fragment-length distributions.

3.1.2 Duplicate handling and alignment considerations

MNase-seq libraries often contain PCR duplicates due to amplification steps. Removing duplicates can reduce biases from overrepresented molecules, but overly aggressive duplicate removal may discard genuine biological signal when coverage is limited. Alignment considerations include how mapping quality is filtered and how read pairs are handled, especially for paired-end sequencing where the fragment boundaries can be inferred from concordant alignments.

A consistent policy for duplicates and alignment filters should be applied across all samples being compared.

3.2 Read alignment to the reference genome

3.2.1 Choosing aligners and parameter basics

Reads are mapped using genome aligners that accommodate short fragments and handle indels appropriately. Parameter choices such as mismatch allowances, seed length, and treatment of soft clipping can influence mapping rates and bias toward particular regions. For nucleosome-scale analyses, alignment confidence matters because fragment end positions help determine nucleosome centers and occupancy profiles.

Aligners differ in speed and behavior, but the key requirement is stable performance across the dataset and clear documentation of settings.

3.2.2 Handling multimapping and repetitive regions

Repetitive genomic segments can produce multimapping reads. Pipelines commonly apply strategies such as restricting to uniquely mapped reads or using probabilistic assignment where appropriate. Excluding multimappers reduces noise but can remove signal from repetitive regions that may be biologically relevant. Alternatively, retaining them can introduce ambiguity and inflate coverage.

The decision typically reflects both the research focus and the quality of the genome annotation and repeat masking.

3.3 Nucleosomal fragment filtering

3.3.1 Selecting mononucleosome-sized reads

After alignment, reads are filtered by inferred fragment length to enrich for mononucleosome-derived fragments. Length thresholds can be based on the known selection window from laboratory size filtering, as well as on observed fragment-length distributions in the sequencing data. Additional filters may exclude reads likely to originate from subnucleosomal or di-nucleosomal fragments.

This stage is crucial for reducing contamination that would otherwise blur nucleosome position calling.

3.3.2 Filtering strategies for bias reduction

Beyond length, filtering can include quality thresholds for mapping, restrictions on fragment ends near mappability gaps, and exclusion of reads in low-complexity regions. Some workflows also adjust fragment start-end definitions depending on library chemistry and strand orientation. The goal is to reduce systematic biases that could be misinterpreted as biological occupancy changes.

Well-documented filtering criteria are important for reproducibility, especially when comparing across studies.

3.4 QC metrics

3.4.1 Fragment length distribution checks

The fragment-length distribution is assessed both before and after filtering. A strong mononucleosome peak indicates appropriate digestion and selection. The relative abundance of subnucleosomal fragments suggests under-recovery due to over-digestion or selection issues, while an elevated di-nucleosome fraction can indicate incomplete digestion or a too-permissive selection window.

If the peak shifts across replicates, it can point to digestion inconsistencies or library preparation variability.

3.4.2 Nucleosome occupancy and coverage uniformity

Coverage uniformity is evaluated by examining how reads distribute across the genome and within mappable regions. Sparse datasets can produce noisy occupancy tracks, while overly uneven coverage may indicate technical problems such as poor library complexity, mapping issues, or sampling bias. Summaries might include fraction of reads mapping to the genome, duplication rates, and the coverage breadth across known nucleosome-rich regions.

3.4.3 Periodicity signals and signal-to-noise

Periodicity can be assessed by aggregate analyses of nucleosome-centered profiles or by methods that quantify periodic spacing. Strong signal-to-noise supports effective mononucleosome enrichment and reasonably consistent nucleosome-scale structure. Weak periodicity does not automatically mean failure—biological chromatin states can differ—but it can help diagnose problems with digestion or filtering.

3.5 Batch effects and normalization

3.5.1 Scaling strategies across libraries

Libraries differ in sequencing depth and potentially in enrichment efficiency. Normalization aims to make occupancy comparisons meaningful by accounting for differences in total read counts and, when available, external controls. Common scaling approaches use per-library normalization factors derived from total mapped reads, counts in a subset of the genome, or spike-in signals if incorporated.

Normalization should be applied consistently so that downstream differences reflect biology rather than measurement scale.

3.5.2 Comparing conditions and replicates

Comparisons across conditions require careful handling of replicate structure. Statistical models often incorporate replicate variability rather than treating libraries as independent observations without uncertainty. Exploratory plots, correlation analyses, and clustering based on normalized signals can reveal outliers. When correlations are low, re-checking digestion QC, filtering choices, and mapping consistency is typically warranted before interpreting biological differences.

4 Downstream analysis: nucleosome positioning

4.1 Building nucleosome maps

4.1.1 Calling nucleosome dyads and centers

Nucleosomes are often represented by a dyad or center coordinate derived from the aligned fragment ends. Since a mono-nucleosomal fragment spans DNA around one nucleosome, the dyad position can be inferred by offsetting from fragment start and/or end positions using empirical or protocol-specific conventions. Pipelines differ slightly in how they define the dyad for forward and reverse strands, but the overall objective is to align sequencing signal to a consistent nucleosome coordinate.

Once dyads are defined, they can be used to generate nucleosome positioning tracks across the genome.

4.1.2 Generating occupancy tracks

Occupancy tracks quantify how frequently a nucleosome is detected at each genomic position (or within a small window around dyad coordinates). Aggregating read alignments into fixed bins produces coverage-like profiles. These tracks can then be smoothed or transformed for visualization, but the underlying representation should match the statistical assumptions of any downstream tests.

Occupancy maps enable comparisons between genomic regions and between experimental conditions.

4.2 Positioning metrics

4.2.1 Signal averaging around genomic features

A common approach is to average nucleosome signals around annotated genomic landmarks, such as transcription start sites or other gene-related features. Averaging reduces noise and highlights general positioning tendencies, though it can mask variability across different loci. The choice of window size and alignment to genomic coordinates influences the interpretability of these average profiles.

Care is taken to ensure that the genomic feature definitions are consistent with the genome build used for mapping.

4.2.2 Calling phased nucleosome arrays (conceptual)

Phased arrays refer to coordinated nucleosome placement at a fixed spacing relative to a reference point, yielding periodic positioning. Detecting such arrays often relies on assessing periodicity strength and alignment to expected repeat lengths. Because chromatin is frequently heterogeneous, phased calling is usually probabilistic or uses criteria that trade sensitivity for specificity.

In conceptual terms, stronger phasing suggests more structured nucleosome organization, while weak phasing indicates irregular placement or mixed chromatin states.

4.3 Differential nucleosome occupancy

4.3.1 Defining windows and statistical comparisons

Differential analysis typically partitions the genome into windows (fixed-size bins or feature-centered regions) and tests whether occupancy differs between conditions. Statistical frameworks may model count data with negative binomial-like assumptions or apply other normalization and dispersion estimates. Window selection should consider nucleosome-scale resolution and the expected extent of nucleosome shifts.

Interpretation is sensitive to multiple testing correction and to whether replicate variability is incorporated.

4.3.2 Interpreting changes in occupancy vs. positioning

Changes in occupancy indicate a gain or loss of nucleosome presence, while shifts in positioning can reflect how nucleosomes slide along DNA without necessarily changing their overall count. Some analyses focus on occupancy differences, others on positioning shifts, and some attempt to separate these effects. Distinguishing these facets requires careful choice of metrics and visualization.

A key caution is that differential patterns can be influenced by digestion efficiency, fragment selection, or sequence-dependent cleavage behavior.

4.4 Sequence bias assessment

4.4.1 GC bias and sequence-dependent MNase sensitivity

MNase digestion can vary with DNA sequence composition, including GC content and local sequence features, which can lead to systematic differences in fragment recovery independent of nucleosome presence. Assessing GC bias involves comparing fragment or coverage patterns to sequence composition across the genome. If uncorrected, sequence-dependent bias may masquerade as chromatin changes.

Some pipelines mitigate bias by filtering or by modeling sequence effects, especially when interpreting differential occupancy.

Nucleosome positioning can correlate with DNA motifs and sequence patterns that influence nucleosome affinity and stability. However, because MNase sensitivity is also sequence-dependent, motif interpretations must consider whether observed differences reflect true nucleosome placement or differential digestability. When feasible, analyses compare nucleosome signal patterns with motif annotations and evaluate consistency across replicates and conditions.

Motif-related results are most credible when supported by consistent occupancy and positioning changes that cannot be fully explained by sequence cleavage properties.

5 Downstream analysis: chromatin features

5.1 Nucleosome occupancy around regulatory regions

5.1.1 Promoters and gene bodies (typical analyses)

MNase-seq is often used to study chromatin around promoters and within gene bodies. Typical analyses examine occupancy profiles over transcription start sites, promoter-proximal regions, and across gene-length intervals. Patterns can suggest whether nucleosomes are enriched, depleted, or repositioned in association with transcriptional regulation. Because gene annotation quality and strand conventions affect coordinate mapping, analyses rely on consistent gene models.

Results are usually interpreted alongside independent measures of transcriptional activity and chromatin state where available.

5.1.2 Enhancers and other cis-regulatory elements (general)

For enhancers and other cis-regulatory elements, nucleosome occupancy patterns can reflect regulatory chromatin architecture. Analyses may involve intersecting nucleosome maps with regulatory annotations and comparing average occupancy around these regions. Since enhancer definitions can vary by dataset and cell type, conclusions depend on the accuracy of the regulatory annotations used.

As with promoters, it is important to avoid overinterpreting occupancy changes without considering possible sequence-dependent MNase effects.

5.2 Chromatin accessibility at the nucleosome level

5.2.1 Relating MNase digestion pattern to accessibility

Although MNase-seq is not identical to dedicated accessibility assays, the digestion pattern can be used to infer differences in how DNA is protected within nucleosomal and flanking regions. Regions with altered accessibility may show shifts in mononucleosome enrichment, fragment-size distributions, or signal strength near nucleosome centers. Aggregated patterns around nucleosome edges can sometimes provide a proxy for how tightly DNA is wrapped or how nucleosome stability differs.

The interpretation is inherently indirect because MNase cleavage reflects both protection and enzyme susceptibility.

5.2.2 Distinguishing occupancy from digestability

A core analytical challenge is separating whether a change in signal arises from altered nucleosome presence (occupancy) or from changes in how susceptible the DNA is to digestion (digestability). Sequence bias and chromatin context can affect digestability, so comparisons that interpret “accessibility” should be supported by fragment-length and bias analyses. Some workflows compare patterns across different digestion levels or incorporate models of fragment distribution to better disentangle these contributions.

A cautious framing is typically used: MNase-seq can provide nucleosome-scale accessibility-related information, but it should not be treated as a direct measure of open chromatin in the same way as assays that specifically target accessible ends.

5.3 Higher-order chromatin organization (practical angles)

5.3.1 Fragment periodicity as an indicator of structure

Higher-order organization can sometimes be inferred indirectly from periodicity and organization of nucleosome spacing. Strong periodic patterns in specific genomic contexts can indicate more regular packing or structured nucleosome arrays. Aggregate analyses may reveal periodicity that is not obvious at single-locus resolution.

However, periodicity signals are influenced by chromatin state, digestion conditions, and sequence effects, so they are often treated as supportive evidence rather than definitive structural proof.

5.3.2 Limitations in interpreting large-scale organization

Large-scale chromatin folding and long-range interactions are not directly measured by MNase-seq. Fragment periodicity and local nucleosome patterns can suggest certain structural regimes, but they do not capture 3D genome contacts. Interpreting higher-order organization therefore requires careful restraint and, ideally, integration with complementary assays that probe chromatin structure more directly.

Overinterpretation is a common risk when periodicity is treated as a direct readout of long-range architecture.

6.1 MNase-seq with different digestion strategies

Variants of MNase-seq include modifications to digestion time, enzyme concentration, and quenching steps, as well as changes in selection criteria for fragment sizes. Some protocols aim for gentler digestion to better preserve subnucleosomal information, while others emphasize enrichment of mono-nucleosome fragments. Others may optimize for specific types of chromatin, such as actively transcribed regions that may display different nucleosome stability.

Different strategies can alter what is most strongly represented in sequencing data.

6.2 MNase titration strategies across organisms/cell types

Digestion calibration is not universal across organisms or across cell types. Chromatin composition, nucleosome repeat lengths, and histone variants can affect MNase sensitivity and the resulting fragment sizes. Therefore, titration strategies often involve preliminary pilot digestions, followed by selection of an endpoint that produces a robust mononucleosome signal. In comparative studies, differences in digestion efficiency must be addressed so that inferred chromatin differences are not artifacts of protocol tuning.

A careful approach includes documenting the titration endpoint and ensuring replicates share comparable settings.

6.3 Comparisons to ATAC-seq and DNase-seq (methodological contrasts)

ATAC-seq and DNase-seq measure chromatin accessibility by targeting exposed DNA and detecting fragment ends generated by transposase or DNase activity. These methods directly focus on accessibility signatures, whereas MNase-seq emphasizes nucleosome-protected fragments enriched by MNase digestion. Consequently, MNase-seq and accessibility assays can correlate but may also diverge due to differences in what each enzyme captures and how sequence and chromatin context affect cleavage.

Comparative interpretation often considers whether changes reflect nucleosome presence, nucleosome stability, or general exposure to cleavage.

6.4 Micrococcal nuclease vs. other chromatin nuclease approaches

Other nucleases used for chromatin mapping can differ in cleavage preferences, fragment size patterns, and sensitivity to nucleosome stability. Protocols that use alternative nucleases may capture different aspects of nucleosome protection or chromatin accessibility. As a result, cross-protocol comparisons require careful normalization and interpretation.

Even when the general goal is nucleosome-scale mapping, enzyme choice can shape the effective signal.

7 Common pitfalls and troubleshooting

7.1 Under- vs. over-digestion effects

Under-digestion can lead to excessive di-nucleosome or higher-order fragments, reducing the clarity of mono-nucleosome-based positioning. Over-digestion can consume nucleosome-protected DNA, shifting the fragment distribution toward smaller sizes and decreasing signal quality for nucleosome calling. Gel-based monitoring during optimization is typically the best early indicator of whether digestion level is appropriate.

If fragment profiles differ substantially between replicates, downstream results may reflect technical variation.

7.2 Library artifacts and PCR bias

Library artifacts can include uneven adapter ligation, preferential amplification, or poor complexity due to low input material. PCR bias may distort occupancy by overrepresenting particular fragments, particularly when amplification cycles are high. Evaluating duplication rates, library complexity metrics, and replicate concordance can help diagnose these issues.

If QC shows unexpectedly high duplication or poor mapping, revisiting input quality and amplification strategy is often necessary.

7.3 Mapping challenges and repeat regions

Repeated sequences and segmental duplications can complicate alignment. If a substantial fraction of reads map ambiguously, filtering choices can alter effective coverage, potentially biasing occupancy estimates. Additionally, mismatches between genome builds or annotation versions and the reference used for mapping can lead to systematic coordinate issues.

A remedy typically involves consistent reference usage, repeat-aware mapping settings, and careful reporting of mapping and filtration policies.

7.4 Over-interpretation of digestion “accessibility”

A frequent pitfall is treating MNase-derived patterns as direct accessibility measurements. Because MNase cleavage depends on both nucleosome occupancy and sequence-dependent digestability, changes in digestion profiles can reflect multiple mechanisms. Interpreting “accessibility” without assessing fragment-length shifts and bias can lead to incorrect biological conclusions.

A cautious stance frames MNase-seq as providing nucleosome-scale protection and enrichment patterns, with indirect accessibility inference.

7.5 Reproducibility and documentation checklist

Reproducibility depends on documenting key experimental and computational choices, including digestion conditions, quenching procedures, size selection boundaries, library preparation chemistry, sequencing platform, mapping reference genome, alignment parameters, and all filtering criteria. Quality control results—especially fragment-length distributions and replicate correlations—should be recorded. Sharing analysis code and specifying normalization methods supports consistent reanalysis.

A structured checklist helps reduce accidental variability between experiments and across labs.

8 Reporting and reproducibility

8.1 Methods to include in publications

Publications using MNase-seq typically include enough protocol detail for other researchers to reproduce key steps: sample preparation type (cells vs. nuclei), whether crosslinking was used, digestion optimization approach and endpoint selection, DNA purification and size selection strategy, library construction workflow, sequencing platform and read configuration, and mapping reference version. For transparency, the fragment size window used for mononucleosome enrichment should be specified.

If multiple experimental variants were used, each should be clearly described.

8.2 Reporting QC and normalization choices

QC reporting commonly includes fragment-length distribution summaries, mapping statistics (including fraction mapped and alignment quality filters), duplication or complexity metrics, replicate correlations, and any periodicity assessment if used. Normalization strategies must be explicitly stated, such as how read counts are scaled and how batch effects or spike-ins are handled.

For differential analyses, reporting the statistical approach, window definitions, multiple-testing corrections, and thresholds for significance improves interpretability.

8.3 Sharing data, metadata, and analysis code

Reproducibility is strengthened by depositing sequencing data in public repositories and providing comprehensive metadata: experimental conditions, replicate identifiers, digestion parameters, library batch identifiers, and any normalization references. Sharing analysis code, including preprocessing, alignment settings, filtering, and statistical scripts, enables independent verification. Where possible, providing processed intermediate outputs and clear pipeline versions reduces the effort needed to rerun workflows.

Good reporting also includes describing how results were visualized and how figures were generated from the underlying data.