A scalable Tn5-based method for genome-wide DNA methylation profiling in development and disease.
The 18 matches · 4 of them tie a paragraph to a whole file, not to given lines: weak matches, whose lines are not tinted
- [1] § Methods › EM-seq quantification of DNA methylation reduction upon inhibitor treatment ↔ figures/extfig5a_methylkit_bpDNAm.R, lines 45–80 · score 0.95 · bismarkCoverage, hi.perc, lo.count, lo.perc, methRead, methylKit
- [2] § Methods › EM-seq quantification of DNA methylation reduction upon inhibitor treatment ↔ figures/fig4e_methylKit_sankeyplot.R, lines 1–43 · score 0.90 · bismarkCoverage, hi.perc, lo.count, lo.perc, methRead, MethylDackel
- [3] § Methods › Chromatin state and ChIPseeker annotation of CUT&Tag peaks ↔ figures/extfig2e_chipseeker_addcpg.R, lines 21–60 · score 0.73 · gene body, CpG islands, UTRs, introns, Exons, ChIPseeker
- [4] § Methods › Preprocessing of CmeCUT&Tag-BS/EM ↔ snakePipes_methCnT/snakePipes_WGBS.sh, the whole file · a weak match · score 0.71 · WGBS pipeline, snakePipes, peak calling, aligner, BAM
- [5] § Methods › DNA methylation profiling in differential regions with CmeCUT&Tag-BS ↔ figures/fig4e_methylKit_sankeyplot.R, lines 45–104 · score 0.71 · methylKit, percMethylation, DiffBind regions, mid, 80 %, 20 %
- [6] § Methods › CG count and methylation profiling of CUT&Tag peaks ↔ figures/fig1f_enrichment_vs_WGBS.py, lines 118–120 · score 0.69 · 10–20 %, 90–100 %, bins, enrichment
- [7] § Methods › Pre-processing, alignment, and normalization of CUT&Tag data ↔ snakePipes_methCnT/snakePipes_WGBS.sh, the whole file · a weak match · score 0.69 · bamCoverage, snakePipes, aligner, deepTools, spike, genome
- [8] § Methods › Chromatin state and ChIPseeker annotation of CUT&Tag peaks ↔ utils/annotate_peaksets.R, lines 1–42 · score 0.68 · annoDb, annotatePeak, TSS, hs, gene, overlap
- [9] § Results › MBD–Tn5 fusion design for DNA methylation targeting ↔ figures/fig1e_plot_distributions.py, lines 44–67 · score 0.66 · xMeCP2, xMBD2, MBD1, NTD, IDR
- [10] § Results › MBD–Tn5 enriches CpG-rich methylated regions in iPSCs ↔ utils/intervene_constructs.sh, the whole file · a weak match · score 0.63 · xMeCP2, xMBD2, NTD, IDR, CpG, nuclei
- [11] § Methods › DNA methylation profiling in differential regions with CmeCUT&Tag-BS ↔ figures/extfig5a_methylkit_bpDNAm.R, lines 45–80 · score 0.63 · methylKit, percMethylation, DNAme, minCov, coverage
- [12] § Results › MBD–Tn5 fusion design for DNA methylation targeting ↔ utils/intervene_constructs.sh, the whole file · a weak match · score 0.62 · xMeCP2, xMBD2, NTD, IDR, DNA
- [13] § Results › MBD–Tn5 enriches CpG-rich methylated regions in iPSCs ↔ figures/fig1e_plot_distributions.py, lines 44–67 · score 0.61 · xMeCP2, xMBD2, NTD, IDR, CpG
- [14] § Results › MBD–Tn5 fusion design for DNA methylation targeting ↔ figures/extfig2e_chipseeker_addcpg.R, lines 96–138 · score 0.60 · MBDseq, MeCP2, RRBS, H9, MeDIP, histone
- [15] § Methods › crossNN prediction of tumor biopsies ↔ figures/fig4e_methylKit_sankeyplot.R, lines 1–43 · score 0.58 · Bismark coverage, MethylDackel, BS, hg38, filtering
- [16] § Results › MBD–Tn5 enriches CpG-rich methylated regions in iPSCs ↔ figures/extfig2e_chipseeker_addcpg.R, lines 21–60 · score 0.58 · gene bodies, CpG island, exons, promoters, overlap
- [17] § Methods › Single-cell CmeCUT&Tag ↔ sc_CmeCUT-Tag/sc_cme_cut-tag.Rmd, lines 194–197 · score 0.57 · FindClusters, FindNeighbors, clustering, Tag
- [18] § Methods › Comparison of MBD constructs ↔ figures/extfig2e_chipseeker_addcpg.R, lines 1–19 · score 0.53 · CpG island, UCSC, track, Browser, hg38, Genome
Paper
Loaded from Europe PMC by your browser, not stored by OSCR: doi.org · Europe PMC
The paper is loaded when this pane is shown.
The authors' code
R · 140 lines · 5.1 KB · no license · 4 matches
- # Load libraries
- library(ChIPseeker)
- library(org.Hs.eg.db)
- library(GenomicFeatures)
- library(dplyr)
- library(rtracklayer)
- library(clusterProfiler)
- library(stringr)
- library(ggplot2)
- library(purrr)
- txdb_ensembl <- makeTxDbFromEnsembl(organism="Homo sapiens", release='113')
- # Load CpG island annotation (see https://github.com/Bioconductor/AnnotationHub/issues/42)
- library(rtracklayer)
- session <- browserSession()
- genome(session) <- "hg38"
- query <- ucscTableQuery(session, table="cpgIslandExt")
- cpg_islands <- track(query) # GRanges object
- seqlevelsStyle(cpg_islands) <- seqlevelsStyle(txdb_ensembl)
- # Define a helper function to customize the annotation dataframe
- customize_anno <- function(peakAnno){
- # extract the annotation
- df <- peakAnno@anno %>% as.data.frame()
- # clean annotation
- df2 <- df %>%
- mutate(
- # remove trailing "(<=1kb)" etc
- annotation_minimal = annotation %>%
- str_remove_all("\\(.*?\\)") %>% # remove any parentheses content
- str_remove_all("\\s"),
- # extract numbers if present: "exon 1 of", "intron 2 of", etc.
- exon_no = str_extract(annotation, "(?<=exon )\\d+"),
- intron_no = str_extract(annotation, "(?<=intron )\\d+"),
- annotation_keepfirst = case_when(
- !is.na(exon_no) & exon_no == "1" ~ "1stExon",
- !is.na(exon_no) ~ "OtherExon",
- !is.na(intron_no) & intron_no == "1" ~ "1stIntron",
- !is.na(intron_no) ~ "OtherIntron",
- TRUE ~ annotation_minimal
- ),
- annotation_simple = case_when(
- str_detect(annotation, "Promoter") ~ "Promoter",
- str_detect(annotation, "exon|UTR|intron") ~ "GeneBody",
- TRUE ~ "Intergenic"
- )
- )
- # add cpg island info
- gr <- as.GRanges(peakAnno)
- hits <- findOverlaps(gr, cpg_islands)
- cpg_flag <- rep("-CpGisland", length(gr))
- cpg_flag[unique(queryHits(hits))] <- "+CpGisland"
- df2$CpG_island <- cpg_flag
- df2$annotation_minimal_CpG <- interaction(df2$annotation_minimal, df2$CpG_island, sep = '')
- df2$annotation_keepfirst_CpG <- interaction(df2$annotation_keepfirst, df2$CpG_island, sep = '')
- df2$annotation_simple_CpG <- interaction(df2$annotation_simple, df2$CpG_island, sep = '')
- return(df2)
- }
- # read anno
- peakAnnoList_all <- readRDS('all_peakAnnoList.rds')
- dfs <- lapply(peakAnnoList_all, customize_anno)
- saveRDS(dfs, "all_peakAnno_annotDFs.rds")
- # plot annotation of the list
- plot_annotation_bar <- function(dfs, col="annotation_simple_CpG"){
- if (is.null(names(dfs))) names(dfs) <- paste0("dataset_", seq_along(dfs))
- summarize_categories <- function(df, dataset, col) {
- df %>%
- count(.data[[col]], name = "n") %>%
- mutate(dataset = dataset, prop = n / sum(n)) %>%
- rename(category = !!col)
- }
- anno_bar_df <- map2_dfr(dfs, names(dfs), summarize_categories, col)
- p <- ggplot(anno_bar_df, aes(x = dataset, y = prop, fill = category)) +
- geom_bar(stat = "identity", position = "fill") +
- scale_y_continuous(labels = scales::percent_format()) +scale_fill_brewer(palette = "Set2") +
- labs(x = NULL, y = "Proportion", fill = "Annotation",
- title = "Annotation composition across datasets") +
- coord_flip() +
- theme_bw()
- return(list(
- anno_bar_df = anno_bar_df,
- plot=p
- ))
- }
- anno <- plot_annotation_bar(dfs)
- anno$plot
- # Subset for more meaningful plots
- histones <- dfs[names(dfs) %>% str_detect("H3")]
- cpg <- dfs[names(dfs) %>% str_detect("cpg")]
- bs <- dfs[names(dfs) %>% str_detect("ES")]
- cg_bs <- bs[names(bs) %>% str_detect("CG")]
- ch_bs<- bs[names(bs) %>% str_detect("CH")]
- rrbs <- bs[names(bs) %>% str_detect("RRBS")]
- high_bs <- dfs[names(dfs) %>% str_detect("CGN.h9_hESC_80-101|ES_80-101")]
- mbdseq <- dfs[names(dfs) %>% str_detect("MBDseq")]
- medip <- dfs[names(dfs) %>% str_detect("roadmap")]
- nuclei_constructs <- dfs[c('nuclei_2xMBD2.broadPeak_optimal','nuclei_2xMeCP2.broadPeak_optimal','nuclei_4xMeCP2.broadPeak_optimal','nuclei_NTD-MeCP2-IDR.broadPeak_optimal')]
- gDNA_constructs <- dfs[c('gDNA_2xMBD2.broadPeak_optimal','gDNA_2xMeCP2.broadPeak_optimal','gDNA_4xMeCP2.broadPeak_optimal','gDNA_NTD-MeCP2-IDR.broadPeak_optimal')]
- mbd2 <- dfs[c('gDNA_2xMBD2.broadPeak_optimal','nuclei_2xMBD2.broadPeak_optimal')]
- mbd2_intersect <- dfs[names(dfs) %>% str_detect("MBD2")]
- selected <- c(
- high_bs, cpg, histones, nuclei_constructs
- )
- names(selected) <- sub(".broadPeak_optimal", "", names(selected))
- pdf('cutomized_highBS_histone_nucleiConstructs.pdf',width = 8, height = 6)
- plot_annotation_bar(selected)$plot
- dev.off()
- selected <- c(
- cg_bs, rrbs, mbdseq, medip, mbd2
- )
- names(selected) <- sub(".broadPeak_optimal", "", names(selected))
- names(selected) <- sub("_chipr_optimal_filtered", "", names(selected))
- pdf('cutomized_MethCnT_MeDIP_MBDseq.pdf',width = 8, height = 6)
- plot_annotation_bar(selected)$plot
- dev.off()
- selected <- c(
- high_bs, mbdseq, medip, mbd2
- )
- names(selected) <- sub(".broadPeak_optimal", "", names(selected))
- names(selected) <- sub("_chipr_optimal_filtered", "", names(selected))
- pdf('cutomized_high_MethCnT_MeDIP_MBDseq.pdf', width = 8, height = 6)
- plot_annotation_bar(selected)$plot
- dev.off()
extfig2e_chipseeker_addcpg.R at commit f0f1063, no license · at the source
Overview
- Ecole Polytechnique Federale de Lausanne (EPFL), School of Life Sciences, Brain Mind Institute, EpiGN—NeuroNA Chair in Epigenomics of Neurodevelopmental Disorders, Station 19, Lausanne, Switzerland
- Campus Biotech, Chemin des Mines 9, Geneve, Switzerland
- UniversitätsSpital Zürich, Schmelzbergstrasse 12, Zürich, Switzerland
Abstract
DNA methylation is a key epigenetic modification involved in development and disease, including cancer, and serves as a biomarker for diagnosis. Current detection methods, such as bisulfite sequencing, provide base-pair resolution but require high sequencing depth and cost. Here, we developed CmeCUT&
Reproduced under the paper's license (CC BY), from the paper cited above.
Repositories
Its files are read in the Code ↔ Paper reader above, with 18 matches between paragraphs and lines of code.
EpiGN-EPFL/CmeCUT-Tag
f0f106319447018fc16798402e2051cbb7af3af7, 26 February 2026Availability: 1 check, the latest on 28 September 2026: the link answers
- 28 September 2026: the link answers
19 files
- figures/
extfig2e_chipseeker_addc , R, 140 lines, 4 matchespg.R - figures/
extfig5a_methylkit_bpDNA , R, 82 lines, 2 matchesm.R - figures/
extfig7c_WGBS_boxplot.R , R, 105 lines - figures/
fig1b_plot_colorbar.py , Python, 99 lines - figures/
fig1e_plot_distributions , Python, 346 lines, 2 matches.py - figures/
fig1f_enrichment_vs_WGBS , Python, 190 lines, 1 match.py - figures/
fig3c_multiBSummary_plot , Python, 200 linesPCA.py - figures/
fig3e_plot_sorted_colorb , Python, 65 linesar.py - figures/
fig4e_methylKit_sankeypl , R, 104 lines, 3 matchesot.R - sc_CmeCUT-Tag/
sc_cme_cut-tag.Rmd , R, 234 lines, 1 match - snakePipes_methCnT/
snakePipes_WGBS.sh , Shell, 67 lines, 2 matches - utils/
annotate_peaksets.R , R, 92 lines, 1 match - utils/
chipr.sh , Shell, 12 lines - utils/
chromHMM_plot.R , R, 154 lines - utils/
demuxlet_cellline.sh , Shell, 18 lines - utils/
diffbind.R , R, 58 lines - utils/
intervene_constructs.sh , Shell, 50 lines, 2 matches - utils/
subsample_bamfiles.sh , Shell, 20 lines - README.md, Text, 110 lines
Zenodo 19555724
Availability: 1 check, the latest on 28 September 2026: the link answers (HTTP 200)
- 28 September 2026: the link answers (HTTP 200)
19 files
- figures/
extfig2e_chipseeker_addc , R, 140 linespg.R - figures/
extfig5a_methylkit_bpDNA , R, 82 linesm.R - figures/
extfig7c_WGBS_boxplot.R , R, 105 lines - figures/
fig1b_plot_colorbar.py , Python, 99 lines - figures/
fig1e_plot_distributions , Python, 346 lines.py - figures/
fig1f_enrichment_vs_WGBS , Python, 190 lines.py - figures/
fig3c_multiBSummary_plot , Python, 200 linesPCA.py - figures/
fig3e_plot_sorted_colorb , Python, 65 linesar.py - figures/
fig4e_methylKit_sankeypl , R, 104 linesot.R - sc_CmeCUT-Tag/
sc_cme_cut-tag.Rmd , R, 234 lines - snakePipes_methCnT/
snakePipes_WGBS.sh , Shell, 67 lines - utils/
annotate_peaksets.R , R, 92 lines - utils/
chipr.sh , Shell, 12 lines - utils/
chromHMM_plot.R , R, 154 lines - utils/
demuxlet_cellline.sh , Shell, 18 lines - utils/
diffbind.R , R, 58 lines - utils/
intervene_constructs.sh , Shell, 50 lines - utils/
subsample_bamfiles.sh , Shell, 20 lines - README.md, Text, 110 lines
Code availability
All generated code is available on GitHub: https://
Reproduced under the paper's license (CC BY), from the paper cited above.
Tracing map
Proposed by the machine: these links were found in the paper and verified at the source, without human review. The map will receive a Zenodo DOI once one of the paper's authors has validated it with their ORCID.
What the map holds:
- 2 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 36 scripts, each with its path and the digest of its content;
- 18 matches between paragraphs of the paper and lines of the code (method lexical-v1);
- neither the text of the paper nor the code itself.
Its JSON (tracing-map.json) is deposited on Zenodo with its DOI once the map is validated.
Data
Datasets cited
- geo:GSE320203, at NCBI GEO; found in “Data availability”
Data Availability Statement
The data supporting the findings of this study are available from the corresponding authors upon request. The data generated in this study have been deposited in the GEO database under accession code GSE320203 (https://
All generated code is available on GitHub: https://
Reproduced under the paper's license (CC BY), from the paper cited above.
Versions
The history of this record: each version stored by the harvester or made by a correction of its authors or of the maintainers of its code, and what changed in its facts. The texts of the paper (its abstract, its availability statements) are not part of it; versions that changed only those are not listed.
Version 1, 28 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 7 authors, 2 keywords, 8 MeSH terms, 1 funder, 64 references.
Cite
This paper
Hu, H., Simonet, N., Bilgiç, E. N., Murray, H., Reimann, R., Rechsteiner, M., & Zenk, F. (2026). A scalable Tn5-based method for genome-wide DNA methylation profiling in development and disease. Nature communications, 17(1), 6736. https://
BibTeX
@article{hu2026scalable,
author = {Hu, Hanrong and Simonet, Nahuel and Bilgiç, Ece Naz and Murray, Heather and Reimann, Regina and Rechsteiner, Markus and Zenk, Fides},
title = {{A scalable Tn5-based method for genome-wide DNA methylation profiling in development and disease}},
journal = {Nature communications},
year = {2026},
month = may,
volume = {17},
number = {1},
pages = {6736},
publisher = {Nature Publishing Group},
issn = {2041-1723},
doi = {10.1038/
url = {https://
pmid = {42173849},
pmcid = {PMC13385352}
}
RIS
TY - JOUR
AU - Hu, Hanrong
AU - Simonet, Nahuel
AU - Bilgiç, Ece Naz
AU - Murray, Heather
AU - Reimann, Regina
AU - Rechsteiner, Markus
AU - Zenk, Fides
TI - A scalable Tn5-based method for genome-wide DNA methylation profiling in development and disease
T2 - Nature communications
J2 - Nat Commun
PY - 2026
DA - 2026/
VL - 17
IS - 1
SP - 6736
SN - 2041-1723
PB - Nature Publishing Group
DO - 10.1038/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1038/
"type": "article-journal",
"title": "A scalable Tn5-based method for genome-wide DNA methylation profiling in development and disease",
"container-title": "Nature communications",
"author": [
{
"family": "Hu",
"given": "Hanrong"
},
{
"family": "Simonet",
"given": "Nahuel"
},
{
"family": "Bilgiç",
"given": "Ece Naz"
},
{
"family": "Murray",
"given": "Heather"
},
{
"family": "Reimann",
"given": "Regina"
},
{
"family": "Rechsteiner",
"given": "Markus"
},
{
"family": "Zenk",
"given": "Fides"
}
],
"container-title-short":
"volume": "17",
"issue": "1",
"page": "6736",
"DOI": "10.1038/
"PMID": "42173849",
"PMCID": "PMC13385352",
"ISSN": "2041-1723",
"publisher": "Nature Publishing Group",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
5,
22
]
]
}
}
The tracing map gets a citation of its own once an author has validated it and it has a DOI.
Similar papers
The papers with a page that share the most with this one: the tools found in their code, their categories, datasets, cited references and authors, the rarest counting most.
- [1] doi:10.1016/j.celrep.2026.117073 [code]
- Single-cell epigenomics uncovers heterochromatin instability and transcription factor dysfunction during mouse brain aging.Journal: Cell reportsIn common: Snakemake, Subread (featureCounts), SAMtools, 11 other tools, genetics / omics, 2 references
- [2] doi:10.1038/s41467-026-71803-3 [code]
- Charting the transition from in vitro gliogenesis to the in vivo maturation of human glial progenitor cells transplanted into the hypomyelinated mouse brain.Journal: Nature communicationsIn common: SAMtools, pheatmap, Seurat, 8 other tools, genetics / omics, 6 references
- [3] doi:10.1038/s41586-026-10512-9 [code]
- Astrocyte glucocorticoid receptor signalling restricts neuronal plasticity.Journal: NatureIn common: Subread (featureCounts), SAMtools, clusterProfiler, 8 other tools, 5 references
- [4] doi:10.1126/sciadv.aed2952 [code]
- Activation of transposable elements is linked to a region- and cell type-specific interferon response in Parkinson's disease.Journal: Science advancesIn common: Subread (featureCounts), SAMtools, clusterProfiler, 10 other tools, 2 references
- [5] doi:10.1186/s13059-026-04177-w [code]
- Genomic sequence evolution underlying human neocortical interareal diversification.Journal: Genome biologyIn common: Snakemake, SAMtools, pheatmap, 9 other tools, genetics / omics, 3 references
- [6] doi:10.1038/s42003-026-10957-8 [code]
- Brain defence by the extracellular matrix protein Cochlin.Journal: Communications biologyIn common: Subread (featureCounts), SAMtools, clusterProfiler, 11 other tools, zebrafish
- [7] doi:10.1038/s41467-026-69944-6 [code]
- Multi-modal dissection of cell-type specific TDP-43 pathology in the motor cortex.Journal: Nature communicationsIn common: SAMtools, clusterProfiler, pheatmap, 8 other tools, genetics / omics, 4 references
- [8] doi:10.1093/bioinformatics/btag592 [code]
- Network-based stratification of allele-specific expression reveals patient subgroups in Huntington's disease.Journal: Bioinformatics (Oxford, England)In common: Subread (featureCounts), SAMtools, clusterProfiler, 10 other tools, genetics / omics, 1 reference
- [9] doi:10.1016/j.xcrm.2026.102766 [code]
- A longitudinal single-cell and spatial multiomic atlas of pediatric high-grade glioma.Journal: Cell reports. MedicineIn common: clusterProfiler, pheatmap, Seurat, 9 other tools, genetics / omics, 2 references
- [10] doi:10.21203/rs.3.rs-9927928/v1 [code]
- Genome-wide and allele-resolved maps of the radial architecture of the mouse genomeJournal: Research Square (preprint)In common: Subread (featureCounts), SAMtools, clusterProfiler, 8 other tools, genetics / omics, 1 reference
Contribute
The authors of this paper can claim it, correct its record and validate its tracing map, and the maintainers of its code (its owner, or a public member of its organization) correct what it says of their repository; anyone signed in can ask for its removal. Every request goes to OSCR's own machine, which answers it; your account page follows them.
Sign in with ORCID to claim this paper as one of its authors, correct its record or validate its tracing map: when the paper's metadata lists your ORCID iD, you are recognized at once. Maintainers of its code: sign in with GitHub, then claim the repository on your account page.
Claim this paper
Correct its record
Say what each link of this record is, remove the ones that are not the paper's, add the ones that are missing. The correction becomes a new version of the record, in its Versions section.
Validate its tracing map
You validate the map as this page shows it: 2 repositories of the authors' code, each at its verified commit and with its license, 36 scripts, and 18 matches between paragraphs and code (see the Code and Map sections). It then receives a DOI on Zenodo, with you (your ORCID iD) and OSCR as its creators; the code itself is not deposited.
The map's fingerprint: sha256:4d343e56dc672679…
Add the badge to its README
The badge links the code to this page. Copy one of these into the README of the paper's code: only you decide where it goes, and nothing is changed for you.
Markdown
[, paste the snippet at the top, then “Commit changes…” and, to review it first, “Create a new branch and start a pull request”. You open the pull request; OSCR asks for no permission.
Request its removal
To ask OSCR to remove this record, the copies of its authors' scripts or its tracing map, use the removal request page: signed in, you say who you are, what to remove and why, then review and confirm the request. Published rules decide every request (how).
Discussion, reproductions, activity
Discussion: questions and error reports about this paper and its code, from signed-in readers and its authors. It opens with sign-in.
Reproductions: reports from readers who ran the authors' code: what they reproduced, with which environment, commit and data. It opens with sign-in.
Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.
