Extract nucleosome positions from ATAC-seq fragment size patterns using ATACseqQC or NucleoATAC to understand chromatin structure at promoters and regulatory regions.
BiocManager::install(c('ATACseqQC', 'TxDb.Hsapiens.UCSC.hg38.knownGene'))pip install nucleoatacTell your AI agent what you want to do:
- "Analyze nucleosome positioning from my ATAC-seq data"
- "Generate fragment size distribution to check nucleosome periodicity"
"Plot the fragment size distribution from my ATAC-seq BAM to visualize nucleosome patterns"
"Call nucleosome positions from my ATAC-seq data using NucleoATAC"
"Analyze nucleosome positioning around TSS to visualize NFR and +1/-1 nucleosomes"
"Extract nucleosome-free region reads and analyze coverage at promoters"
- Calculate fragment size distribution to verify nucleosome periodicity
- Separate reads by fragment size (NFR, mono-, di-, tri-nucleosomal)
- Call nucleosome positions using NucleoATAC or ATACseqQC
- Generate nucleosome occupancy signal tracks
- Visualize nucleosome patterns around TSS
| Fragment Size | Origin | Meaning |
|---|---|---|
| <100 bp | Nucleosome-free | Active promoter/enhancer |
| 180-247 bp | Mono-nucleosome | Well-positioned nucleosome |
| 315-473 bp | Di-nucleosome | Regular spacing |
| 558-615 bp | Tri-nucleosome | Chromatin structure |
| Pattern | Meaning |
|---|---|
| Strong NFR peak | Active promoter/enhancer |
| Regular spacing | Well-positioned nucleosomes |
| Fuzzy positioning | Dynamic chromatin |
- Paired-end data is required for accurate fragment sizing
- Use properly deduplicated BAM files
- Filter for high-quality fragments (MAPQ > 30)
- High mitochondrial contamination can skew fragment distributions
- NucleoATAC provides base-pair resolution nucleosome positions
- ATACseqQC is useful for V-plots around TSS showing nucleosome phasing