Comparative analysis of nuclei isolation methods for brain single-nucleus RNA sequencing.
The 7 matches · 3 of them tie a paragraph to a whole file, not to given lines: weak matches, whose lines are not tinted
- [1] § Results › Method-specific variations in glial gene expression ↔ scripts/analysis/Fig4.UCell.R, the whole file · a weak match · score 0.80 · disease associated microglia, C1qa, UCell, homeostatic microglia, Hexb, Mbp
- [2] § Results › Isolation technique influences the cell types captured in snRNA-seq ↔ scripts/analysis/Fig2.Dotplot.R, the whole file · a weak match · score 0.78 · Slc17a7, Sv2b, Slc1a3, Gad1, Gad2, Gja1
- [3] § Results › Method-specific variations in glial gene expression ↔ scripts/analysis/Fig4.UCell.R, the whole file · a weak match · score 0.75 · C1qa, UCell, disease associated, Glul, homeostatic, Tmem119
- [4] § STAR★Methods › Quantification and statistical analysis › snRNA-seq data analysis ↔ vignettes/SAHA_vignette.Rmd, lines 62–102 · score 0.69 · Allen Brain Atlas, Cell Atlas, FindAllMarkers, RStudio, pre, Mouse
- [5] § STAR★Methods › Quantification and statistical analysis › snRNA-seq data analysis ↔ scripts/preprocess/Step1_IsolationPreProcess_DFandSX.R, lines 485–574 · score 0.66 · scMCA, elbow, Seurat, prefix, scoring, mt
- [6] § Results › Isolation technique influences the cell types captured in snRNA-seq ↔ scripts/analysis/Fig2.PropsandROGUE.R, lines 56–107 · score 0.59 · excitatory neuron, inhibitory neuron, ROGUE, oligodendrocytes, assay, astrocytes
- [7] § Results › Differences in quality-control metrics among isolation protocols ↔ scripts/preprocess/Step1_IsolationPreProcess_DFandSX.R, lines 367–421 · score 0.55 · SoupX, ambient RNA, transformed, rho, doublets, subset
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 · 80 lines · 3.4 KB · no license · 2 matches
- set.seed(1221)
- setwd("/N/project/kim_lab/Kersey_Holly/IsolationCollab")
- library(UCell)
- library(Seurat)
- library(ggplot2)
- library(gridExtra)
- #load seurat object
- isolation <- readRDS("20231206_snIsolation_collab.rds")
- #subset object based on cell type
- Idents(isolation) <- "simple_anno"
- isolation.microglia <- subset(isolation, idents="Microglia")
- isolation.astrocyte <- subset(isolation, idents="Astrocyte")
- isolation.oligo <- subset(isolation, idents="Oligo")
- #define marker genes for astrocytes
- markers.as <- list()
- markers.as$HM_AS = c("Kcnj10","Glul","Slc1a2","Slc1a3","Slc6a11","Slc16a1","Ldha","Srebf1")
- markers.as$DA_AS = c("Hmgb1","Hmgb3","Hmgb2","Id3","Gfap","S100b","Cd81","Cebpa","Ptprg")
- #define marker genes for microglia
- markers.mg <- list()
- markers.mg$HM_MG = c("Tmem119","P2ry12","P2ry13","Cx3cr1", "C1qa", "Csf1r", "Hexb")
- markers.mg$DA_MG= c("Apoe", "Cst7", "Trem2", "Itgax", "B2m", "Cst7")
- #define marker genes for oligos
- markers.og <- list()
- markers.og$HM_OL = c("Olig1","Olig2","Mog","Mbp","Mobp","Plp1","Sox10","Gpr37","Mag","Cnp","Myrf")
- #set identity to isolation method
- Idents(isolation.microglia) <- "tech"
- #Compute UCell scores
- isolation.microglia <- AddModuleScore_UCell(isolation.microglia, features = markers.mg)
- signature.names.mg <- paste0(names(markers.mg), "_UCell")
- P <- VlnPlot(isolation.microglia, features = signature.names.mg, pt.size = FALSE)
- write.csv(P[[1]][['data']], "ucell.HM_MG.csv") #save homeostatic microglia scores
- write.csv(P[[2]][['data']], "ucell.DA_MG.csv") #save disease-associated scores
- Idents(isolation.astrocyte) <- "tech"
- isolation.astrocyte <- AddModuleScore_UCell(isolation.astrocyte, features = markers.as)
- signature.names.as <- paste0(names(markers.as), "_UCell")
- Q <- VlnPlot(isolation.astrocyte, features = signature.names.as, pt.size = FALSE)
- write.csv(Q[[1]][['data']], "ucell.HM_AS.csv") #save homeostatic microglia scores
- write.csv(Q[[2]][['data']], "ucell.DA_AS.csv") #save disease-associated microglia scores
- Idents(isolation.oligo) <- "tech"
- isolation.oligo <- AddModuleScore_UCell(isolation.oligo, features = markers.og)
- signature.names.og <- paste0(names(markers.og), "_UCell")
- R <- VlnPlot(isolation.oligo, features = signature.names.og, pt.size = FALSE)
- write.csv(R[[1]][['data']], "ucell.HM_OG.csv") #save homeostatic oligo scores
- #For visualization purposes
- p1 <- VlnPlot(isolation.oligo, features = "HM_OL_UCell", pt.size = FALSE)+
- stat_summary(fun.y = median, geom='point', size = 2, colour = "black")+
- ylim(0,0.7)+
- theme_classic()+NoLegend()
- p2 <- VlnPlot(isolation.microglia, features = "DA_MG_UCell", pt.size = FALSE)+
- stat_summary(fun.y = median, geom='point', size = 2, colour = "black")+
- ylim(0,0.7)+
- theme_classic()+NoLegend()
- p3 <- VlnPlot(isolation.microglia, features = "HM_MG_UCell", pt.size = FALSE)+
- stat_summary(fun.y = median, geom='point', size = 2, colour = "black")+
- ylim(0,0.7)+
- theme_classic()+NoLegend()
- p4 <- VlnPlot(isolation.astrocyte, features = "DA_AS_UCell", pt.size = FALSE)+
- stat_summary(fun.y = median, geom='point', size = 2, colour = "black")+
- ylim(0,0.7)+
- theme_classic()+NoLegend()
- p5 <- VlnPlot(isolation.astrocyte, features = "HM_AS_UCell", pt.size = FALSE)+
- stat_summary(fun.y = median, geom='point', size = 2, colour = "black")+
- ylim(0,0.7)+
- theme_classic()+NoLegend()
- pdf(file = "Ucell.pdf")
- grid.arrange(p5,p4,p3,p2,p1, nrow=1,ncol=5)
- dev.off()
- #Plots were edited in Adobe Illustrator for stylistic components, i.e., titles, colors...
Fig4.UCell.R at commit b39d9d3, no license · at the source
Overview
- Medical Neurosciences Graduate Program, Indiana University School of Medicine, Indianapolis, IN, USA
- Stark Neurosciences Research Institute, Indiana University School of Medicine, Indianapolis, IN, USA
- Department of Medical and Molecular Genetics, Indiana University School of Medicine, Indianapolis, IN, USA
- Department of Psychological and Brain Sciences, Indiana University, Bloomington, IN, USA
Abstract
Single-nucleus RNA sequencing (snRNA-seq) enables resolving cellular heterogeneity in complex tissues by using nuclei instead of cells, overcoming limitations of single-cell RNA sequencing and enabling analysis of frozen and hard-to-isolate tissues. Despite advances in isolation techniques, systematic evaluations of their effects on nuclear integrity and subsequent data quality remain lacking, a critical gap with profound implications for rigor and reproducibility. To address this, we compared three mechanistically distinct nuclei isolation strategies with brain tissue: a sucrose gradient centrifugation-based method, a spin column-based method, and a machine-assisted platform. All methods captured diverse cell types but revealed considerable protocol-dependent differences in cell type proportions, transcriptional homogeneity, and the preservation of cell-state-specific markers. Moreover, workflows differentially influenced contamination levels from ambient, mitochondrial, and ribosomal RNAs, with the machine-assisted method exhibiting the highest overall data quality. Our findings establish nuclei isolation methodology as a critical experimental variable shaping snRNA-seq data quality and biological interpretation.
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 7 matches between paragraphs and lines of code.
jungsukimlab/NucleiIsolation
b39d9d3c75f9413a41d8855f480d05633f4d1852, 5 April 2026Availability: 1 check, the latest on 30 September 2026: the link answers
- 30 September 2026: the link answers
10 files
- scripts/
analysis/ , R, 78 linesFig2.AnnotationandUMAP.R - scripts/
analysis/ , R, 52 lines, 1 matchFig2.Dotplot.R - scripts/
analysis/ , R, 107 lines, 1 matchFig2.PropsandROGUE.R - scripts/
analysis/ , R, 152 linesFig3.MitoandRibo.R - scripts/
analysis/ , R, 241 linesFig4.MarkerComparison.R - scripts/
analysis/ , R, 80 lines, 2 matchesFig4.UCell.R - scripts/
preprocess/ , R, 575 lines, 2 matchesStep1_IsolationPreProces s_DFandSX.R - scripts/
preprocess/ , R, 33 linesStep2_job_IntegrationPos tSXandDFv2.R - scripts/
preprocess/ , R, 193 linesStep3_Clustering2Determi neResolution.R - README.md, Text, 20 lines
Zenodo 14040699
Availability: 1 check, the latest on 28 September 2026: the link answers (HTTP 200)
- 28 September 2026: the link answers (HTTP 200)
31 files
- R/
AutoAnnotate.R , R, 95 lines - R/
CorrelateDS.R , R, 34 lines - R/
Create_MarkerBased_Viz.R , R, 163 lines - R/
Create_MarkerFree_Viz.R , R, 116 lines - R/
Create_SAHA_object.R , R, 98 lines - R/
Create_SelfSimilarity_Vi , R, 67 linesz.R - R/
Downsample.R , R, 32 lines - R/
Initialize_MarkerFree.R , R, 51 lines - R/
Initialize_Markers.R , R, 60 lines - R/
Initialize_Self_Similari , R, 46 linesty.R - R/
Investigate_MarkerFree.R , R, 46 lines - R/
Investigate_Marker_Based , R, 106 lines.R - R/
Investigate_SelfSimilari , R, 54 linesty.R - R/
Marker_Diversity.R , R, 67 lines - R/
Marker_Richness.R , R, 40 lines - R/
NormalizeDS.R , R, 104 lines - R/
Run_Marker_Based.R , R, 81 lines - R/
SAHA.R , R, 37 lines - R/
SAHA_lookup_cluster.R , R, 18 lines - R/
SemiAutoAnnotate.R , R, 199 lines - R/
SimilarityDend.R , R, 43 lines - R/
Tune_Markers.R , R, 43 lines - R/
call_SAHA_plots.R , R, 37 lines - R/
data.R , R, 39 lines - R/
hello.R , R, 18 lines - R/
zzz.R , R, 142 lines - tests/
testthat.R , R, 12 lines - tests/
testthat/ , R, 6 linestest-dummy_SAHA_qs.R - vignettes/
SAHA_vignette.Rmd , R, 380 lines - LICENSE, License, 21 lines
- README.md, Text, 51 lines
neurogenetics/saha
0ced41731fb3200fd16d1e415b33d088ad42b695, 13 May 2026Availability: 1 check, the latest on 28 September 2026: the link answers
- 28 September 2026: the link answers
34 files
- R/
AutoAnnotate.R , R, 114 lines - R/
CorrelateDS.R , R, 51 lines - R/
Create_MarkerBased_Viz.R , R, 236 lines - R/
Create_MarkerFree_Viz.R , R, 147 lines - R/
Create_SAHA_object.R , R, 110 lines - R/
Create_SelfSimilarity_Vi , R, 68 linesz.R - R/
Downsample.R , R, 36 lines - R/
Generate_SAHA_Report.R , R, 40 lines - R/
Initialize_MarkerFree.R , R, 49 lines - R/
Initialize_Markers.R , R, 66 lines - R/
Initialize_Self_Similari , R, 56 linesty.R - R/
Investigate_MarkerFree.R , R, 46 lines - R/
Investigate_Marker_Based , R, 111 lines.R - R/
Investigate_Self_Similar , R, 54 linesity.R - R/
Marker_Diversity.R , R, 70 lines - R/
Marker_Richness.R , R, 53 lines - R/
NormalizeDS.R , R, 125 lines - R/
Run_Marker_Based.R , R, 77 lines - R/
SAHA.R , R, 77 lines - R/
SAHA2Seurat.R , R, 50 lines - R/
SAHA_lookup_cluster.R , R, 17 lines - R/
SemiAutoAnnotate.R , R, 386 lines - R/
Seurat2SAHA.R , R, 55 lines - R/
SimilarityDend.R , R, 43 lines - R/
Summarize_Params.R , R, 125 lines - R/
Tune_Markers.R , R, 68 lines - R/
call_SAHA_plots.R , R, 45 lines - R/
data.R , R, 39 lines - R/
zzz.R , R, 143 lines - tests/
testthat.R , R, 12 lines - tests/
testthat/ , R, 6 linestest-dummy_SAHA_qs.R - vignettes/
SAHA_vignette.Rmd , R, 380 lines, 1 match - LICENSE, License, 21 lines
- README.md, Text, 70 lines
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:
- 3 repositories of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
- 70 scripts, each with its path and the digest of its content;
- 7 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:GSE290858, at NCBI GEO; found in “Data and code availability”
Data and code availability
The snRNA-seq data generated in this study have been deposited in the GEO database: GSE290858 (https://
The code for the snRNA-seq analysis in this study has been deposited at https://
Any additional information required to reanalyze the data reported in this study is available from the lead contact upon request.
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, 30 September 2026: the first record
Recorded: type, language, journal, volume, issue, pages, dates, 7 authors, 5 keywords, 7 MeSH terms, 7 funders, 66 references, 1 RRID.
Cite
This paper
Kersey, H. N., Acri, D. J., Dabin, L. C., Hartigan, K. A., Mustaklem, R., Park, J. H., & Kim, J. (2026). Comparative analysis of nuclei isolation methods for brain single-nucleus RNA sequencing. Cell reports methods, 6(3), 101337. https://
BibTeX
@article{kersey2026compa
author = {Kersey, Holly N and Acri, Dominic J and Dabin, Luke C and Hartigan, Kelly A and Mustaklem, Richard and Park, Jung Hyun and Kim, Jungsu},
title = {{Comparative analysis of nuclei isolation methods for brain single-nucleus RNA sequencing}},
journal = {Cell reports methods},
year = {2026},
month = mar,
volume = {6},
number = {3},
pages = {101337},
publisher = {Elsevier},
issn = {2667-2375},
doi = {10.1016/
url = {https://
pmid = {41875869},
pmcid = {PMC13030979}
}
RIS
TY - JOUR
AU - Kersey, Holly N
AU - Acri, Dominic J
AU - Dabin, Luke C
AU - Hartigan, Kelly A
AU - Mustaklem, Richard
AU - Park, Jung Hyun
AU - Kim, Jungsu
TI - Comparative analysis of nuclei isolation methods for brain single-nucleus RNA sequencing
T2 - Cell reports methods
J2 - Cell Rep Methods
PY - 2026
DA - 2026/
VL - 6
IS - 3
SP - 101337
SN - 2667-2375
PB - Elsevier
DO - 10.1016/
UR - https://
LA - en
ER -
CSL-JSON
{
"id": "10.1016/
"type": "article-journal",
"title": "Comparative analysis of nuclei isolation methods for brain single-nucleus RNA sequencing",
"container-title": "Cell reports methods",
"author": [
{
"family": "Kersey",
"given": "Holly N"
},
{
"family": "Acri",
"given": "Dominic J"
},
{
"family": "Dabin",
"given": "Luke C"
},
{
"family": "Hartigan",
"given": "Kelly A"
},
{
"family": "Mustaklem",
"given": "Richard"
},
{
"family": "Park",
"given": "Jung Hyun"
},
{
"family": "Kim",
"given": "Jungsu"
}
],
"container-title-short":
"volume": "6",
"issue": "3",
"page": "101337",
"DOI": "10.1016/
"PMID": "41875869",
"PMCID": "PMC13030979",
"ISSN": "2667-2375",
"publisher": "Elsevier",
"URL": "https://
"language": "en",
"issued": {
"date-parts": [
[
2026,
3,
1
]
]
}
}
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-10382-x [code]
- ScQCenrich enables multi-metric quality control for single-cell RNA sequencing.Journal: Communications biologyIn common: circlize, ComplexHeatmap, Seurat, 3 other tools, genetics / omics, mouse, cellular / molecular, 7 references
- [2] doi:10.1038/s41514-026-00391-9 [code]
- Region-specific transcriptional signatures of brain aging in the absence of neuropathology at the single-cell level.Journal: npj agingIn common: circlize, ComplexHeatmap, Seurat, 6 other tools, genetics / omics, cellular / molecular, 5 references
- [3] doi:10.1101/gr.281113.125 [code]
- Single-nucleus multiomic profiling of the aging mouse substantia nigra reveals conserved gene alterations linked to Parkinson's disease.Journal: Genome researchIn common: circlize, ComplexHeatmap, Seurat, 6 other tools, genetics / omics, mouse, cellular / molecular, 4 references
- [4] doi:10.1038/s41593-026-02367-0 [code]
- A reproducible three-dimensional model of human brain tissue to investigate physiological and disease-associated microglia phenotypes.Journal: Nature neuroscienceIn common: circlize, ComplexHeatmap, Seurat, 6 other tools, cellular / molecular, 4 references
- [5] doi:10.1038/s42003-026-10034-0 [code]
- Region- and cell type-specific changes in gene expression in the cerebellum after classical fear conditioning.Journal: Communications biologyIn common: circlize, ComplexHeatmap, Seurat, 5 other tools, genetics / omics, mouse, cellular / molecular, 4 references
- [6] 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: circlize, ComplexHeatmap, Seurat, 5 other tools, genetics / omics, 4 references
- [7] 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: circlize, ComplexHeatmap, Seurat, 6 other tools, genetics / omics, cellular / molecular, 3 references
- [8] doi:10.1038/s41586-026-10214-2 [code]
- Multidimensional profiling of heterogeneity in supratentorial ependymomas.Journal: NatureIn common: circlize, ComplexHeatmap, Seurat, 6 other tools, genetics / omics, mouse, 2 references
- [9] doi:10.1038/s41467-026-76232-w [code]
- Th17 effector cytokines induce shared and distinct microglial and endothelial cell responses in a mouse model for post-streptococcal encephalitis.Journal: Nature communicationsIn common: circlize, ComplexHeatmap, Seurat, 5 other tools, genetics / omics, mouse, cellular / molecular, 2 references
- [10] doi:10.1038/s41467-026-74171-0 [code]
- Cluster replicability in single-cell and single-nucleus atlases of the mouse brain.Journal: Nature communicationsIn common: circlize, ComplexHeatmap, Seurat, 3 other tools, genetics / omics, mouse, cellular / molecular, 4 references
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: 3 repositories of the authors' code, each at its verified commit and with its license, 70 scripts, and 7 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:00b8cd7e3f1badca…
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.
