OSCR

Protocol for single-cell epigenetic profiling in human organoids and tumoroids with Epi-CyTOF.

Code ↔ Paper

2 matches between paragraphs of the paper and lines of its authors' code, computed by the harvester (lexical-v1). Click a colored paragraph or line to see its counterpart.

The 2 matches · 1 of them tie a paragraph to a whole file, not to given lines: a weak match, whose lines are not tinted
  1. [1] § Step-by-step method details › Quantification and statistical analysis with R ↔ scripts/compute_all_gates.R, the whole file · a weak match · score 0.83 · quantile_val, modality_df, gates_df, df_long, unimodal, Quantification
  2. [2] § Before you begin › Innovation ↔ notebooks/Epi-CyTOF_data_analysis.Rmd, lines 9–59 · score 0.51 · epigenetic profiling, human organoid, CyTOF, single cell, technical, mass

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 · 135 lines · 4.5 KB · MIT · 1 match

  1. #' Compute marker-specific gating thresholds and signal indices
  2. #'
  3. #' This function computes gating thresholds for each marker in a long-format
  4. #' single-cell expression dataset. The gating strategy depends on whether a
  5. #' marker is classified as unimodal or multimodal based on prior modality
  6. #' assessment.
  7. #'
  8. #' For unimodal markers, gates are computed using either a quantile-based
  9. #' approach or a median/standard deviation-based rule. For multimodal markers,
  10. #' density estimation combined with peak detection is used to identify
  11. #' separation points between expression modes.
  12. #'
  13. #' Additionally, a signal index (SI) is calculated to quantify the separation
  14. #' between negative and positive populations.
  15. #'
  16. #' @param df A data.frame in long format containing at least:
  17. #' \describe{
  18. #' \item{marker}{Marker name}
  19. #' \item{expression}{Numeric expression values}
  20. #' }
  21. #'
  22. #' @param modality_df A data.frame containing modality information per marker:
  23. #' \describe{
  24. #' \item{marker}{Marker name}
  25. #' \item{unimodal}{Logical indicating whether the marker is unimodal}
  26. #' }
  27. #'
  28. #' @param lower Numeric lower bound for expression filtering (default: -Inf)
  29. #' @param upper Numeric upper bound for expression filtering (default: Inf)
  30. #' @param bw Bandwidth method for density estimation (default: "nrd0")
  31. #' @param method_unimodal Method for unimodal gating:
  32. #' \code{"quantile"} or \code{"median_sd"}
  33. #' @param quantile_val Quantile threshold used for gating when
  34. #' \code{method_unimodal = "quantile"} (default: 0.05)
  35. #' @param k_sd Number of standard deviations used when
  36. #' \code{method_unimodal = "median_sd"} (default: 2)
  37. #'
  38. #' @return A data.frame with one row per marker containing:
  39. #' \describe{
  40. #' \item{marker}{Marker name}
  41. #' \item{gate}{Estimated gating threshold}
  42. #' \item{SI}{Signal index measuring separation between negative and positive populations}
  43. #' }
  44. #'
  45. #' @details
  46. #' The signal index (SI) is defined as:
  47. #' \deqn{SI = (median(pos) - median(neg)) / (sd(pos) + sd(neg))}
  48. #'
  49. #' Density-based peak detection is performed using the \pkg{pracma} package.
  50. #'
  51. #' @importFrom stats density median sd quantile
  52. #' @importFrom pracma findpeaks
  53. #'
  54. #' @examples
  55. #' \dontrun{
  56. #' gates_df <- compute_all_gates_2(df_long, modality_df)
  57. #' }
  58. #'
  59. #' @export
  60. compute_all_gates <- function(df, modality_df,
  61. lower = -Inf, upper = Inf,
  62. bw = "nrd0",
  63. method_unimodal = c("quantile", "median_sd"),
  64. quantile_val = 0.05,
  65. k_sd = 2) {
  66. method_unimodal <- match.arg(method_unimodal)
  67. markers <- unique(df$marker)
  68. gates_SI <- list()
  69. for (m in markers) {
  70. expr <- df$expression[df$marker == m]
  71. is_unimodal <- modality_df$unimodal[modality_df$marker == m]
  72. gate_val <- NA
  73. SI <- NA
  74. if (length(expr) > 0 && !is.na(is_unimodal) && is_unimodal) {
  75. expr_in_bounds <- expr[expr >= lower & expr <= upper]
  76. if (method_unimodal == "quantile") {
  77. gate_val <- unname(quantile(expr_in_bounds, quantile_val, na.rm = TRUE))
  78. } else if (method_unimodal == "median_sd") {
  79. gate_val <- median(expr_in_bounds, na.rm = TRUE) -
  80. k_sd * sd(expr_in_bounds, na.rm = TRUE)
  81. }
  82. neg <- expr[expr < gate_val]
  83. pos <- expr[expr >= gate_val]
  84. if (length(neg) > 0 && length(pos) > 0) {
  85. SI <- (median(pos) - median(neg)) / (sd(pos) + sd(neg))
  86. }
  87. } else if (length(expr) > 0) {
  88. dens <- density(expr, bw = bw)
  89. inv_y <- max(dens$y) - dens$y
  90. minima <- pracma::findpeaks(inv_y)
  91. if (!is.null(minima)) {
  92. min_idx <- minima[, 2]
  93. peak_idx <- pracma::findpeaks(dens$y)[, 2]
  94. gate_idx <- min_idx[min_idx > min(peak_idx) & min_idx < max(peak_idx)]
  95. gate_val <- if (length(gate_idx) == 0) NA else dens$x[gate_idx]
  96. gate_val <- gate_val[gate_val >= lower & gate_val <= upper]
  97. gate_val <- sort(gate_val)[1]
  98. if (!is.na(gate_val)) {
  99. neg <- expr[expr < gate_val]
  100. pos <- expr[expr >= gate_val]
  101. if (length(neg) > 0 && length(pos) > 0) {
  102. SI <- (median(pos) - median(neg)) / (sd(pos) + sd(neg))
  103. }
  104. }
  105. }
  106. }
  107. gates_SI[[m]] <- list(gate = gate_val, SI = SI)
  108. }
  109. gate_vals <- sapply(gates_SI, function(x) x$gate)
  110. SI_vals <- sapply(gates_SI, function(x) x$SI)
  111. data.frame(marker = markers, gate = gate_vals, SI = SI_vals)
  112. }

compute_all_gates.R at commit 139f966, under MIT · at the source

Overview

Authors: Ezgi Senoglu1, Vivian Mittné2, Sevina Dietz1, Anne Eugster1, Sebastian Thieme3, Franziska Baenke2,4,5, Mareike Albert1, Claudia Peitzsch1
  1. Center for Regenerative Therapies Dresden (CRTD), TUD Dresden University of Technology, 01307 Dresden, Germany
  2. Department of Visceral, Thoracic and Vascular Surgery, Medical Faculty and University Hospital Carl Gustav Carus, Dresden University of Technology, 01307 Dresden, Germany
  3. Department of Pediatrics, Faculty of Medicine and University Hospital Carl Gustav Carus, Technische Universität Dresden, 01307 Dresden, Germany
  4. National Center for Tumor Diseases (NCT/UCC), 01037 Dresden, a partnership between DKFZ, Faculty of Medicine and University Hospital Carl Gustav Carus, TUD Dresden University of Technology, and Helmholtz-Zentrum Dresden-Rossendorf (HZDR), 01328 Dresden, Germany
  5. German Cancer Consortium (DKTK), Partner Site Dresden, 01307 Dresden, and German Cancer Research Center (DKFZ), 69120 Heidelberg, Germany
Journal: STAR protocols, volume 7, issue 3, article 104695
Dates: published online 31 July 2026; in print July 2026
Type: Research article · Language: English
License: CC BY
Identifiers: DOI 10.1016/j.xpro.2026.104695 · PMID 42541720 · PMCID PMC13453457 · OpenAlex W7171968201
Open access: gold, a free copy (OpenAlex)
Status: code verified
Categories: human (organism), other condition (population)
Methods: Statistics, Smoothing, state filtering, decompositions, Machine learning
Keywords: Bioinformatics, cell biology, single cell, mass cytometry, developmental biology, cancer, Genomics, Immunology, molecular biology, Neuroscience, organoids
Topic: Cancer Cells and Metastasis (Oncology, Medicine), according to OpenAlex
Funding: CyTOF; CRTD; TU Dresden; RRID (SCR_027145); Ezio Bonifacio; German Center for Diabetes Research; Deutsche Forschungsgemeinschaft (DFG) (FZT 111); Federal Ministry for Research; Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) (01EW2208); Sächsisches Staatsministerium für Wissenschaft und Kunst; Stanford University; Center for Regenerative Therapies TU Dresden; Emmy Noether Programme (2231/1-1, SPP EPIADAPT 563285578, AL 2231/4-1, 2231/3-1); German Centers for Health Research; DZG
Citations: not cited yet (Europe PMC); 55 references in the paper

Abstract

In this protocol, we present Epi-CyTOF, a cytometry by time of flight (CyTOF) approach for single-cell epigenome profiling in cortical organoids as a model system for human brain development. We describe steps for in silico panel design, antibody conjugation and testing, sample preparation and staining, and high-dimensional data analysis. Additionally, we provide recommendations for applying Epi-CyTOF to other model systems based on our experience with gastroesophageal tumoroids.

For complete details on the use and execution of this protocol, please refer to Ditzer et al .1

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 2 matches between paragraphs and lines of code.

Stacheltier/EPI-CyTOF

License: MIT
State: the link answers, verified on 27 September 2026
Evidence: files inventoried
Commit: 139f966968412ee930ad19077a72b98d51bf29e0, 30 April 2026
Languages: C/C++ (340), R (216), Shell (2), C (2), JavaScript (1)
Size: 2,742 files, 561 scripts
Software Heritage: not archived
Found in: the text, “Quantification and statistical analysis with R”
Holds: README, environment (renv.lock, renv/library/windows/R-4.5/x86_64-w64-mingw32/askpass/DESCRIPTION, renv/library/windows/R-4.5/x86_64-w64-mingw32/assorthead/DESCRIPTION, renv/library/windows/R-4.5/x86_64-w64-mingw32/backports/DESCRIPTION, renv/library/windows/R-4.5/x86_64-w64-mingw32/Biobase/DESCRIPTION, renv/library/windows/R-4.5/x86_64-w64-mingw32/Cairo/DESCRIPTION, renv/library/windows/R-4.5/x86_64-w64-mingw32/cli/DESCRIPTION, renv/library/windows/R-4.5/x86_64-w64-mingw32/colorRamps/DESCRIPTION, renv/library/windows/R-4.5/x86_64-w64-mingw32/diptest/DESCRIPTION, renv/library/windows/R-4.5/x86_64-w64-mingw32/doParallel/DESCRIPTION, renv/library/windows/R-4.5/x86_64-w64-mingw32/flowCore/DESCRIPTION, renv/library/windows/R-4.5/x86_64-w64-mingw32/FNN/DESCRIPTION), tests, documentation, 57 notebooks
Not found: license file, CITATION.cff, continuous integration
Tools: ggplot2 (18 files), tidyverse (8 files), broom (2 files), lme4 (2 files), patchwork (2 files), survival (2 files), UMAP (2 files), caret (1 file), SingleCellExperiment (1 file)
Availability: 1 check, the latest on 27 September 2026: the link answers
  • 27 September 2026: the link answers
562 files

ParkerICI/premessa

License: GPL-3.0
State: the link answers, verified on 27 September 2026
Evidence: files inventoried
Commit: 68b42bb984637d0f3ad6a0ecc83e9278994afc85, 6 September 2022
Languages: R (15), JavaScript (3)
Size: 63 files, 18 scripts
Software Heritage: not archived
Found in: the text, “De-barcoding of pooled samples after measurement”
Holds: README, license file, environment (DESCRIPTION), documentation
Not found: CITATION.cff, tests, continuous integration
Tools: ggplot2 (4 files), reshape2 (2 files), data.table (1 file)
Availability: 1 check, the latest on 27 September 2026: the link answers
  • 27 September 2026: the link answers
20 files

The paper's code and data availability statement is in the Data section.

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;
  • 579 scripts, each with its path and the digest of its content;
  • 2 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

No dataset and no data link were found in the paper.

Data and code availability

The Epi-CyTOF data for the cortical organoids are available in the FlowRepository: FR-FCM-Z8DC. The data analysis workflow is available within the paper-specific GitHub repository (https://github.com/Stacheltier/EPI-CyTOF).55

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, 27 September 2026: the first record

Recorded: type, language, journal, volume, issue, pages, dates, 8 authors, 11 keywords, 15 funders, 45 references, 35 RRIDs.

Cite

This paper

Senoglu, E., Mittné, V., Dietz, S., Eugster, A., Thieme, S., Baenke, F., Albert, M., & Peitzsch, C. (2026). Protocol for single-cell epigenetic profiling in human organoids and tumoroids with Epi-CyTOF. STAR protocols, 7(3), 104695. https://doi.org/10.1016/j.xpro.2026.104695

BibTeX

@article{senoglu2026protocol,
author = {Senoglu, Ezgi and Mittné, Vivian and Dietz, Sevina and Eugster, Anne and Thieme, Sebastian and Baenke, Franziska and Albert, Mareike and Peitzsch, Claudia},
title = {{Protocol for single-cell epigenetic profiling in human organoids and tumoroids with Epi-CyTOF}},
journal = {STAR protocols},
year = {2026},
month = jul,
volume = {7},
number = {3},
pages = {104695},
publisher = {Elsevier},
issn = {2666-1667},
doi = {10.1016/j.xpro.2026.104695},
url = {https://doi.org/10.1016/j.xpro.2026.104695},
pmid = {42541720},
pmcid = {PMC13453457}
}

RIS

TY - JOUR
AU - Senoglu, Ezgi
AU - Mittné, Vivian
AU - Dietz, Sevina
AU - Eugster, Anne
AU - Thieme, Sebastian
AU - Baenke, Franziska
AU - Albert, Mareike
AU - Peitzsch, Claudia
TI - Protocol for single-cell epigenetic profiling in human organoids and tumoroids with Epi-CyTOF
T2 - STAR protocols
J2 - STAR Protoc
PY - 2026
DA - 2026/07/31
VL - 7
IS - 3
SP - 104695
SN - 2666-1667
PB - Elsevier
DO - 10.1016/j.xpro.2026.104695
UR - https://doi.org/10.1016/j.xpro.2026.104695
LA - en
ER -

CSL-JSON

{
"id": "10.1016/j.xpro.2026.104695",
"type": "article-journal",
"title": "Protocol for single-cell epigenetic profiling in human organoids and tumoroids with Epi-CyTOF",
"container-title": "STAR protocols",
"author": [
{
"family": "Senoglu",
"given": "Ezgi"
},
{
"family": "Mittné",
"given": "Vivian"
},
{
"family": "Dietz",
"given": "Sevina"
},
{
"family": "Eugster",
"given": "Anne"
},
{
"family": "Thieme",
"given": "Sebastian"
},
{
"family": "Baenke",
"given": "Franziska"
},
{
"family": "Albert",
"given": "Mareike"
},
{
"family": "Peitzsch",
"given": "Claudia"
}
],
"container-title-short": "STAR Protoc",
"volume": "7",
"issue": "3",
"page": "104695",
"DOI": "10.1016/j.xpro.2026.104695",
"PMID": "42541720",
"PMCID": "PMC13453457",
"ISSN": "2666-1667",
"publisher": "Elsevier",
"URL": "https://doi.org/10.1016/j.xpro.2026.104695",
"language": "en",
"issued": {
"date-parts": [
[
2026,
7,
31
]
]
}
}

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.1038/s42003-026-10957-8 [code]
Brain defence by the extracellular matrix protein Cochlin.
Journal: Communications biology
In common: survival, caret, UMAP, 7 other tools
[2] doi:10.1093/neuonc/noag128 [code]
Spatially-resolved single-cell imaging of melanoma brain metastases identifies localized immune patterns predictive of immune checkpoint blockade response.
Journal: Neuro-oncology
In common: survival, caret, lme4, 5 other tools, other condition, 1 reference
[3] 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: survival, UMAP, broom, 6 other tools, other condition
[4] doi:10.1016/j.xcrm.2026.102682 [code]
TET CpG sequence-context-specific DNA demethylation shapes progression of IDH-mutant gliomas.
Journal: Cell reports. Medicine
In common: survival, UMAP, broom, 6 other tools, other condition
[5] doi:10.1016/j.xcrm.2026.102766 [code]
A longitudinal single-cell and spatial multiomic atlas of pediatric high-grade glioma.
Journal: Cell reports. Medicine
In common: SingleCellExperiment, UMAP, lme4, 5 other tools, other condition, 1 reference
[6] doi:10.3390/ijms27156925 [code]
XGBoost-SHAP Interpretable Modeling Identifies and Validates an Eight-Gene Biomarker for Hepatic Encephalopathy Risk Prediction in Cirrhosis.
Journal: International journal of molecular sciences
In common: survival, caret, SingleCellExperiment, 5 other tools, other condition
[7] doi:10.1038/s41586-026-10612-6 [code]
Acquired genetic and cell-state changes in IDH-mutant glioma progression.
Journal: Nature
In common: survival, SingleCellExperiment, broom, 5 other tools, other condition
[8] doi:10.1016/j.xcrm.2026.102651 [code]
Integrative CSF profiling identifies disease-specific immune responses in leptomeningeal disease.
Journal: Cell reports. Medicine
In common: SingleCellExperiment, UMAP, broom, 5 other tools, other condition
[9] doi:10.3390/ijms27104466 [code]
Uncovering the Key Circuit FOSL2/FOS/EGR3/EGR1, Contributing to the Hyperexcitability of Excitatory Neurons in the Epileptic Temporal Cortex and Hippocampus.
Journal: International journal of molecular sciences
In common: caret, UMAP, lme4, 5 other tools
[10] doi:10.1126/sciadv.aec9291 [code]
Computational mechanisms of perception in autism revealed using games inspired by rodent operant tasks.
Journal: Science advances
In common: caret, broom, lme4, 5 other tools

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.

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.