OSCR

CellCover defines marker gene panels capturing developmental progression in neocortical neural stem cell identity.

Code ↔ Paper

The paper beside its authors' code: matches between them have not been computed for this paper yet.

Paper

Loaded from Europe PMC by your browser; OSCR does not store it: doi.org · Europe PMC

The paper is loaded when this pane is shown.

The authors' code

R · 94 lines · 2.3 KB · no license

  1. #' @export
  2. CellCover = function(J = 3, alpha = 0.05, genes_per_class = 6000, te=0.1,
  3. mat, CellTypeNames, CellTypeLabels,prev = NULL)
  4. {
  5. selectedtrain = Find_CTGroups(CellTypeLabels, CellTypeNames)
  6. trainfreq = CompPropMatrix(selectedtrain, mat, paral = FALSE)
  7. nct = length(selectedtrain)
  8. ng = nrow(mat)
  9. library(slam)
  10. library(gurobi)
  11. marginmatrix=matrix(NA,ng,nct)
  12. for (currenttype in 1:nct)
  13. {
  14. for (i in 1:ng)
  15. {
  16. temp=mean(trainfreq[i,-currenttype])
  17. marginmatrix[i,currenttype]=trainfreq[i,currenttype]/temp
  18. }
  19. }
  20. currentmarker=list()
  21. for (currenttype in 1:nct)
  22. {
  23. td = sort(marginmatrix[,currenttype], decreasing = TRUE)[genes_per_class]
  24. if (td < 1){
  25. td = 1
  26. }
  27. keepgeneindex=intersect(which(marginmatrix[,currenttype]>td),which(trainfreq[,currenttype]>=te))
  28. G=length(keepgeneindex)
  29. group=selectedtrain[[currenttype]]
  30. tempA=mat[keepgeneindex,group]
  31. N=ncol(tempA)
  32. A=matrix(0,G,N)
  33. for (i in 1:G)
  34. {
  35. A[i,which(tempA[i,]>0)]=1
  36. }
  37. if (is.null(prev) == FALSE)
  38. {
  39. z0 = prev[[currenttype]] #only the index
  40. holder1 = rep(0,ng)
  41. for (i in 1:length(z0))
  42. {
  43. holder1[z0[i]] = 1
  44. }
  45. holder2 = rep(0,G)
  46. for (i in 1:G)
  47. {
  48. holder2[i] = holder1[keepgeneindex[i]]
  49. }
  50. AA = matrix(0, N+1+G,G+N)
  51. AA[1:N,1:G]=matrix(t(A),N,G)
  52. AA[1:N,(G+1):(G+N)]=diag(J,N,N)
  53. AA[N+1,]=c(rep(0,G),rep(1,N))
  54. AA[(N+2):(N+1+G),1:G] = diag(1,G,G)
  55. rhs=c(rep(J,N),alpha*N,holder2)
  56. sense= c(rep('>=',N),'<=',rep('>=',G))
  57. }
  58. else{
  59. AA=matrix(NA,N+1,G+N)
  60. AA[1:N,1:G]=matrix(t(A),N,G)
  61. AA[1:N,(G+1):(G+N)]=diag(J,N,N)
  62. AA[N+1,]=c(rep(0,G),rep(1,N))
  63. rhs = c(rep(J,N),alpha*N)
  64. sense = c(rep('>=',N),'<=')
  65. }
  66. model = list()
  67. model$A = AA
  68. model$obj = c(1/(marginmatrix[keepgeneindex,currenttype]-td+0.00001),rep(0,N))
  69. model$modelsense = 'min'
  70. model$modelname = 'poolserch'
  71. model$rhs = rhs
  72. model$sense = sense
  73. model$vtype = 'B'
  74. params = list()
  75. params$MIPFocus = 0
  76. params$PoolSearchMode = 0
  77. #params$PoolSolutions = 2
  78. solutions = gurobi(model, params = params)
  79. markerlist=keepgeneindex[which(solutions$x[1:G]>0.9)]
  80. currentmarker[[currenttype]]=markerlist
  81. }
  82. return(currentmarker)}

CellCover.R at commit 9c3d331, no license · at the source

Overview

Authors: Lanlan Ji1, An Wang1, Shreyash Sonthalia2, Seungmae Seo3, Daniel Q Naiman1, Laurent Younes1,4, Carlo Colantuoni2,5, Donald Geman1,4
  1. Department of Applied Mathematics and Statistics, Johns Hopkins University Baltimore United States
  2. Departments of Neurology and Neuroscience, Johns Hopkins University Baltimore United States
  3. Department of Natural Sciences, University of Maryland Eastern Shore Princess Anne United States
  4. Center for Imaging Science, Johns Hopkins University Baltimore United States
  5. Institute for Genome Sciences, University of Maryland School of Medicine Baltimore United States
Institutions: Johns Hopkins University (United States); University of Maryland Eastern Shore (United States); University of Maryland, Baltimore (United States)
Journal: eLife, volume 14, article RP107531
Dates: published online 28 July 2026
Type: Research article · Language: English
License: CC BY
Identifiers: DOI 10.7554/elife.107531 · PMID 42517345 · PMCID PMC13412324 · OpenAlex W4415402955
Open access: gold, a free copy (OpenAlex)
Status: code verified (how the code is found and verified)
Categories: genetics / omics (modality), human (organism), mouse (organism), cellular / molecular (subfield)
Methods: Statistics
Keywords: human, primate, mouse, single-cell RNA-seq, marker genes, neocortex
MeSH: Neocortex*, Neural Stem Cells*, Neurogenesis*, Animals, Biomarkers, Gene Expression Regulation, Developmental, Humans, Mice, Neurodevelopment, Sequence Analysis, RNA, Single-Cell Analysis, Single-Cell Gene Expression Analysis (* major topic)
Journal subjects: Neuroscience
Topic: Single-cell and spatial transcriptomics (Molecular Biology, Biochemistry, Genetics and Molecular Biology), according to OpenAlex
Funding: National Cancer Institute (R01CA200859); National Science Foundation (2124230); Data sharing and visualization via NeMO Analytics (R24MH114815, R01DC019370)
Citations: not cited yet (Europe PMC); 53 references in the paper

Abstract

Defining cell classes is central to the analysis of growing single-cell RNA sequencing (scRNA-seq) atlases. Marker genes are most often identified by differential expression (DE) methods that assess genes one at a time, ignoring the redundancy and complementarity revealed when genes are considered jointly. Working with binarized expression data, we instead seek discriminating panels of genes that together are specific to a cell type, framing marker-panel selection as a variant of the minimal set-covering problem in combinatorial optimization. This formulation efficiently searches the vast space of candidate panels, exploits the large cell numbers typical of scRNA-seq, and is robust to zero-inflation. Using blood and brain data, we show that our method, CellCover, reduces gene redundancy and captures cell-class-specific signals distinct from those found by DE. Transfer-learning experiments across mouse, primate, and human data demonstrate that CellCover identifies conserved cell classes in neocortical neurogenesis and tracks developmental progression in progenitors and neurons. Examining outer radial glia markers across mammals, we find that transcriptomic elements of this key cell type likely arose in rodent gliogenic precursors before the full program emerged in the primate lineage.

Reproduced under the paper's license (CC BY), from the paper cited above.

Repository

Its files are read in the Code ↔ Paper reader above.

lanlanji/CoveringPackage

License: none: the authors keep all their rights
State: the link answers, verified on 27 September 2026
Evidence: files inventoried
Commit: 9c3d331e65f0d9fe5860d13abbec219f3dab4059, 12 April 2023
Languages: R (4)
Size: 11 files, 4 scripts
Software Heritage: not archived
Found in: “Data availability”
Holds: README, environment (DESCRIPTION)
Not found: license file, CITATION.cff, tests, continuous integration, documentation
Availability: 1 check, the latest on 27 September 2026: the link answers
  • 27 September 2026: the link answers
5 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:

  • 1 repository of the authors' code, each at its verified commit, with its license and how the link was found in the paper;
  • 4 scripts, each with its path and the digest of its content;
  • no match between paragraphs and code yet;
  • 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 availability

While we have not generated any new data in this report, all the public RNA-seq data used in the analyses is freely explorable at NeMO Analytics (https://nemoanalytics.org/expression.html?gene_symbol=SOX2&gene_symbol_exact_match=1&is_multigene=0&layout_id=CellCover). Links to the data we used and original data sources are available at NeMO Analytics. We invite researchers to interrogate all the public data resources we examine here in the NeMO Analytics multi-omics data exploration environment that is designed to provide biologists with no programming expertise the ability to perform powerful analyses across collections of data in brain development (Sonthalia et al., 2026). Users can visualize individual genes (NeMO: Individual genes in cortex (https://nemoanalytics.org/p?&l=CellCover&g=EOMES) and NeMO: Individual genes in blood (https://nemoanalytics.org/p?&l=immune&g=CSF1R)) or groups of genes simultaneously across multiple datasets, including automated transfer of the CellCover gene marker panels we have identified across cell types and developmental time (NeMO: Telley 3 CellCover Panels (https://nemoanalytics.org/p?p=p&l=CellCover&c=TelleyCellCover3sets&algo=binary), NeMO: Telley 12 CellCover Panels (https://nemoanalytics.org/p?p=p&l=CellCover&c=TelleyCellCover12sets&algo=binary), NeMO: Sorted Brain Cell CellCover Panels (https://nemoanalytics.org/p?p=p&l=CellCover&c=LiuCellCoverSorted&algo=binary), and NeMO: Blood 34 CellCover Panels (https://nemoanalytics.org/p?p=p&l=immune&c=CellCoverHao34&algo=binary)). CellCover is available in CellCover R (https://github.com/lanlanji/CoveringPackage) and CellCover Python (https://pypi.org/project/CellCover/).

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, pages, dates, 8 authors, 6 keywords, 12 MeSH terms, 3 funders, 50 references.

Cite

This paper

Ji, L., Wang, A., Sonthalia, S., Seo, S., Naiman, D. Q., Younes, L., Colantuoni, C., & Geman, D. (2026). CellCover defines marker gene panels capturing developmental progression in neocortical neural stem cell identity. eLife, 14, RP107531. https://doi.org/10.7554/elife.107531

BibTeX

@article{ji2026cellcover,
author = {Ji, Lanlan and Wang, An and Sonthalia, Shreyash and Seo, Seungmae and Naiman, Daniel Q and Younes, Laurent and Colantuoni, Carlo and Geman, Donald},
title = {{CellCover defines marker gene panels capturing developmental progression in neocortical neural stem cell identity}},
journal = {eLife},
year = {2026},
month = jul,
volume = {14},
pages = {RP107531},
publisher = {eLife Sciences Publications, Ltd},
issn = {2050-084X},
doi = {10.7554/elife.107531},
url = {https://doi.org/10.7554/elife.107531},
pmid = {42517345},
pmcid = {PMC13412324}
}

RIS

TY - JOUR
AU - Ji, Lanlan
AU - Wang, An
AU - Sonthalia, Shreyash
AU - Seo, Seungmae
AU - Naiman, Daniel Q
AU - Younes, Laurent
AU - Colantuoni, Carlo
AU - Geman, Donald
TI - CellCover defines marker gene panels capturing developmental progression in neocortical neural stem cell identity
T2 - eLife
J2 - eLife
PY - 2026
DA - 2026/07/28
VL - 14
SP - RP107531
SN - 2050-084X
PB - eLife Sciences Publications, Ltd
DO - 10.7554/elife.107531
UR - https://doi.org/10.7554/elife.107531
LA - en
ER -

CSL-JSON

{
"id": "10.7554/elife.107531",
"type": "article-journal",
"title": "CellCover defines marker gene panels capturing developmental progression in neocortical neural stem cell identity",
"container-title": "eLife",
"author": [
{
"family": "Ji",
"given": "Lanlan"
},
{
"family": "Wang",
"given": "An"
},
{
"family": "Sonthalia",
"given": "Shreyash"
},
{
"family": "Seo",
"given": "Seungmae"
},
{
"family": "Naiman",
"given": "Daniel Q"
},
{
"family": "Younes",
"given": "Laurent"
},
{
"family": "Colantuoni",
"given": "Carlo"
},
{
"family": "Geman",
"given": "Donald"
}
],
"container-title-short": "eLife",
"volume": "14",
"page": "RP107531",
"DOI": "10.7554/elife.107531",
"PMID": "42517345",
"PMCID": "PMC13412324",
"ISSN": "2050-084X",
"publisher": "eLife Sciences Publications, Ltd",
"URL": "https://doi.org/10.7554/elife.107531",
"language": "en",
"issued": {
"date-parts": [
[
2026,
7,
28
]
]
}
}

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/s41593-026-02204-4 [code]
NeMO Analytics: a compendium of transcriptomic data for the exploration of neocortical development.
Journal: Nature neuroscience
In common: genetics / omics, mouse, 17 references, author Carlo Colantuoni
[2] doi:10.1038/s41467-025-61316-w [code]
Early developmental origins of cortical disorders modeled in human neural stem cells.
Journal: Nature communications
In common: genetics / omics, 6 references, author Carlo Colantuoni
[3] doi:10.7554/elife.90325 [code]
Temporal transcriptomic dynamics in developing macaque neocortex.
Journal: eLife
In common: genetics / omics, mouse, 9 references
[4] doi:10.1126/sciadv.adn9998 [code]
Multiple parallel cell lineages in the developing mammalian cerebral cortex.
Journal: Science advances
In common: genetics / omics, mouse, 9 references
[5] doi:10.1093/nsr/nwae213 [code]
Epigenomic landscapes during prefrontal cortex development and aging in rhesus.
Journal: National science review
In common: genetics / omics, 7 references
[6] doi:10.1016/j.xgen.2024.100581 [code]
SIMS: A deep-learning label transfer tool for single-cell RNA sequencing analysis.
Journal: Cell genomics
In common: genetics / omics, mouse, 7 references
[7] doi:10.1038/s41467-025-66069-0 [code]
Syngap1 and the development of murine neocortical progenitor cells.
Journal: Nature communications
In common: mouse, 7 references
[8] doi:10.3389/fcell.2022.884748 [code]
Single-cell and single-nuclei RNA sequencing as powerful tools to decipher cellular heterogeneity and dysregulation in neurodegenerative diseases.
Journal: Frontiers in cell and developmental biology
In common: genetics / omics, cellular / molecular, 7 references
[9] doi:10.1073/pnas.2314802121
ERK signaling expands mammalian cortical radial glial cells and extends the neurogenic period.
Journal: Proceedings of the National Academy of Sciences of the United States of America
In common: mouse, cellular / molecular, 6 references
[10] doi:10.1073/pnas.2416757122
NOTCH, ERK, and SHH signaling respectively control the fate determination of cortical glia and olfactory bulb interneurons.
Journal: Proceedings of the National Academy of Sciences of the United States of America
In common: mouse, cellular / molecular, 6 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.

Add the badge to the code's README

OSCR: the badge of this paper

If you are an author or a maintainer of this paper's code, you can add this badge. It says what OSCR knows of the paper (its code, its validated tracing map) and links the code to this page. Paste one of these snippets under the README's title. Keep the two marks around it, so that the badge can be replaced in place and is never added twice. OSCR changes nothing in your repository.

Markdown
<!-- oscr-badge:start -->
[![OSCR](https://openscicode.org/badge/doi/10.7554/elife.107531.svg)](https://openscicode.org/badge/doi/10.7554/elife.107531)
<!-- oscr-badge:end -->
reStructuredText
.. oscr-badge:start

.. image:: https://openscicode.org/badge/doi/10.7554/elife.107531.svg
   :target: https://openscicode.org/badge/doi/10.7554/elife.107531
   :alt: OSCR

.. oscr-badge:end
HTML
<!-- oscr-badge:start -->
<a href="https://openscicode.org/badge/doi/10.7554/elife.107531"><img src="https://openscicode.org/badge/doi/10.7554/elife.107531.svg" alt="OSCR"></a>
<!-- oscr-badge:end -->

On GitHub, open the README, choose Edit (the pencil), paste the snippet under the title, then "Commit changes".

If you choose to, you can have the badge kept up to date for you, with oscr badge --write from the command line or with the GitHub Action in your repository's workflows. Every other form and option (a card, AsciiDoc, LaTeX, shields.io) is on the badge 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: code errors and mismatches between the code and the paper about this paper and its code, from readers and its authors, beside the paper and its tracing map.

Report a mismatch between the code and the paper · Report a code error (signed in)

Reproductions: reports from readers who ran the authors' code and did not get the paper's result, what they ran, with which environment, commit and data.

Report a reproduction failure (signed in)

Versions of its code: the releases of its code tied to a version of this paper (the preprint, the accepted manuscript, the version of record…), each with the tracing map versioned with it.

Activity: what happens around this paper: new versions of its record, its map's validation, discussions and reproductions. It opens with sign-in.