PreprintCited by 15,508Single-cell and spatial transcriptomicsGene expression and cancer classificationGene Regulatory Network Analysis

The DOI matches its Crossref record; 97 of the 97 references it cites were found in a registry or an index; none of the 47 sources its claims cite is checked for support yet.

Integrated analysis of multimodal single-cell data

Open access

Preprint — published as 10.1016/j.cell.2021.04.048 · Paper page

Abstract

The simultaneous measurement of multiple modalities, known as multimodal analysis, represents an exciting frontier for single-cell genomics and necessitates new computational methods that can define cellular states based on multiple data types. Here, we introduce ‘weighted-nearest neighbor’ analysis, an unsupervised framework to learn the relative utility of each data type in each cell, enabling an integrative analysis of multiple modalities. We apply our procedure to a CITE-seq dataset of hundreds of thousands of human white blood cells alongside a panel of 228 antibodies to construct a multimodal reference atlas of the circulating immune system. We demonstrate that integrative analysis substantially improves our ability to resolve cell states and validate the presence of previously unreported lymphoid subpopulations. Moreover, we demonstrate how to leverage this reference to rapidly map new datasets, and to interpret immune responses to vaccination and COVID-19. Our approach represents a broadly applicable strategy to analyze single-cell multimodal datasets, including paired measurements of RNA and chromatin state, and to look beyond the transcriptome towards a unified and multimodal definition of cellular identity. Availability: Installation instructions, documentation, tutorials, and CITE-seq datasets are available at http://www.satijalab.org/seurat

Overview

AI snapshot

Written from the paper’s full text.

TL;DR

The paper introduces weighted-nearest neighbor (WNN) analysis, an unsupervised framework that learns cell-specific weights for each modality and integrates multimodal single-cell data into a joint neighbor graph. Applied to CITE-seq and other multimodal datasets, WNN improves cell-state resolution, enables a 210,911-cell multimodal PBMC reference atlas, and supports reference-based mapping of vaccination and COVID-19 immune responses.

Key findings

  • WNN integration improved resolution of cell states across datasets: in a CITE-seq dataset of 8,617 cord blood mononuclear cells, it separated CD4+ and CD8+ T cells while retaining cDC identity; in 30,672 human bone marrow mononuclear cells, it resolved T cell and progenitor states not separable by either modality alone.
  • The multimodal PBMC atlas contained 210,911 total cells, including 161,764 10x 3′ cells profiled with 228 antibodies and 49,147 10x 5′ cells; WNN clustering identified 57 clusters, grouped into 8 level-1 and 30 level-2 annotations.
  • Benchmarking against totalVI and MOFA+ showed WNN had superior or equivalent performance for reconstructing 25 proteins and RNA transcriptome features, with speed improvements up to 15-fold on the full dataset.
  • Reference mapping of a COVID-19 scRNA-seq dataset identified MAIT cell depletion, validated by CyTOF with Pearson correlation R = 0.911; vaccination analysis found strong day-3 interferon responses in CD14+ and CD16+ monocytes and cDC2s, including a shared module of 62 interferon-responsive genes and Siglec-1 (CD169) protein induction.

Methods

  • Design: WNN constructs independent k-nearest-neighbor graphs per modality, performs within- and cross-modality prediction, converts prediction accuracy to cell-specific modality weights via a softmax, and builds a weighted nearest-neighbor graph for downstream visualization, clustering, and trajectory analysis.
  • Data: analyses included CITE-seq datasets of 8,617 cord blood mononuclear cells (10 antibodies), 30,672 bone marrow mononuclear cells (25 antibodies), and 210,911 PBMCs (228 antibodies); plus 11,351 10x Multiome ATAC+RNA PBMCs, 4,725 ASAP-seq PBMCs (227 antibodies), and 34,774 SHARE-seq mouse skin cells.
  • Main analysis: WNN was implemented in Seurat v4; clustering used the WNN graph with SLM, differential expression used Wilcoxon tests, and reference mapping used supervised PCA (sPCA) with anchor-based transfer to project query scRNA-seq datasets onto the multimodal reference.

Limitations

  • Stated: WNN requires a dimensional reduction to describe neighborhood structure and is not compatible with categorical or low-dimensional data.
  • Stated: WNN assumes modalities do not define conflicting sets of cell states, which may be problematic when integrating morphological, functional, and molecular data.
  • Stated: the circulating immune atlas was constructed from PBMCs and therefore contains few cells with no nuclei (erythrocytes) or multi-lobed nuclei (granulocytes); evident in the text, developing/differentiated neutrophils absent from the reference could not be successfully mapped, and small marker panels may not isolate truly homogeneous high-resolution subsets.

Model-written summary. Check it against the paper before relying on it.

Integrity signals

Registry

Crossref · 10.1101/2020.10.12.335331

Cites

97 of 97 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Verification

Claims the paper makes

0 of 47 cited sources checked · 0 supported · 0 flagged

  • It may state two separate claims as one. 71
  • Probably not a substantial claim. 38
102 of 102
  • The paper relies on this without support.

    The eight volunteers and three time points are representative of the circulating immune system and vaccination response.

    “We obtained PBMC samples from eight volunteers enrolled in an HIV vaccine trial (Elizaga et al., 2018; Li et al., 2017), with ages spanning from 20–49 years.”
    “For each subject, PBMCs were collected at three time points: immediately before (day 0), 3 days, and 7 days following administration of a VSV-vectored HIV vaccine (Figure 3A).”
  • The paper relies on this without support.

    The 24 human volunteers sampled are representative of the healthy immune system.

    “Because these results are consistent in both pre- and post-vaccination time points, they likely describe general characteristics of the healthy immune system.”
  • The paper relies on this without support.

    The reference dataset's cell states are sufficiently comprehensive to serve as a reference for mapping query PBMC datasets.

    “having constructed a multimodal reference of the immune system, we sought to leverage this dataset to assist in the analysis and interpretation of additional single-cell experiments profiling human PBMCs (queries)”
  • The paper relies on this without support.

    The protein data held out from mapping can serve as ground truth for assessing annotation accuracy.

    “We verified the accuracy of our predictions using the query protein data, which was held out of the reference mapping procedure”
  • The authors state this as an assumption.

    WNN assumes that modalities do not define conflicting sets of cell states.

    “Second, WNN assumes that modalities do not define conflicting sets of cell states.”
  • The paper relies on this without support.

    Cells originating from rare states should not have the same bandwidth constraint as cells originating from abundant states.

    “We considered that cells originating from rare states should not have the same bandwidth constraint as cells originating from abundant states, and therefore considered a modified approach to select kernel bandwidths.”
  • The paper relies on this without support.

    The 20 cells with the lowest non-zero Jaccard similarity to cell i represent cells that exhibit some similarity with cell i but are unlikely to reside in the same molecular state.

    “We expect that these represent cells that exhibit some similarity with cell i, but are unlikely to reside in the same molecular state.”
  • The paper relies on this without support.

    The first LSI dimension is typically correlated with technical metrics in ATAC-seq data.

    “excluding the first dimension as this is typically correlated with technical metrics in ATAC-seq data”
  • The paper relies on this without support.

    If the WNN procedure is performing well, gene expression levels within the examined neighborhoods should represent a homogeneous population with low levels of variability.

    “If the WNN procedure is performing well, gene expression levels within these neighborhoods should represent a ‘homogeneous’ population with low levels of variability.”
  • The paper relies on this without support.

    The statistical power of per-gene differential expression tests across time points within a cell state is heavily dependent on the abundance of the cell state.

    “We note that when calculating DE genes and proteins within a cell state, across experimental time points, the statistical power of these per-gene tests is heavily dependent on the abundance of the cell state.”
10 of 10

Claims other papers make about this one

Each is checked against this paper’s own text: does it say what the citing paper says it does?

Self-check

Percentages, sums and reported differences re-computed from the paper’s own numbers.

Each contradiction between two passages listed here was confirmed by an independent check.

  • Flow-cytometry validation uses a different post-vaccination day

    Major

    The results say flow cytometry validated the shift between day 0 and day 3, but the methods say whole blood for that assay was collected before vaccination and on days 1 and 7-not day 3.

  • NK-cell cluster count changes from six to five

    Minor

    The atlas results count six NK-cell clusters, but the subsequent description of NK-cell heterogeneity says NK cells were subdivided into five clusters.

  • Figure 5B names a CD103-positive, CD49-positive group despite mutually exclusive expression

    Minor

    The caption describes CD103 and CD49a expression as mutually exclusive and elsewhere names the comparison group as CD103-negative/CD49-positive. Its description of Figure 5B instead calls that group CD103-positive/CD49-positive.

  • Targeted marker panels are assigned to different supplementary tables

    Minor

    The results direct readers to Table S2 for the targeted panels, whereas the methods say the panels are reported in Table S3.

  • Conflicting percentages for the same annotation benchmark

    Minor

    For query cells receiving different annotations from Seurat and scArches, the Figure S7 caption reports stronger support for Seurat in 73.8% of cases, while the benchmarking methods report 79.4%.

5 of 5

Citations

Cites

Flagged first, then by citations

Show all 91 references

Cited by

Sorted by citations

See all 15,508

Metadata

Authors
Hao, Y., Hao, S., Andersen-Nissen, E., Mauck, W. M., Zheng, S., Butler, A., Lee, M. J., Wilk, A. J., Darby, C., Zagar, M., Hoffman, P., Stoeckius, M., Papalexi, E., Mimitou, E. P., Jain, J., Srivastava, A., Stuart, T., Fleming, L. B., Yeung, B., Rogers, A. J., McElrath, J. M., Blish, C. A., Gottardo, R., Smibert, P., & Satija, R.
Affiliations
  • New York Genome Center (US)
  • New York University (US)
  • Stanford Medicine (US)
  • BioLegend (United States) (US)
  • Chan Zuckerberg Biohub San Francisco (US)
Corresponding author
  • Smibert, Peter
  • Satija, Rahul
Venue
Cell
Publisher
openRxiv
Published
2020-10-12
Volume
184
Issue
13
Pages
3573-3587.e29
ISSN
0092-8674 (ISSN-L)
OpenAlex
W3164692211
Cited by
15,508
References
91
Citations by year
2026: 2908 · 2025: 3703 · 2024: 3990 · 2023: 2850 · 2022: 1684 · 2021: 345 · 2020: 12
FWCI
731.42
Type
journal article
Language
en
Open access
Locations
  • Cell (publishedVersion, cc-by)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, cc-by-nc-nd)
  • PubMed Central (publishedVersion, cc-by)
  • PubMed (publishedVersion)
Funders
  • Bill and Melinda Gates Foundation (OPP 1113682)
  • Chan Zuckerberg Initiative (HCA-A1704-01895, EOSS-0000000082, U19AI128914)
  • National Institutes of Health (1dp2hg009623-01, OT2OD026673, AI068618)
  • National Human Genome Research Institute (DP2-HG-009623, R21HG009748, RM1HG011014)
  • National Institute of Allergy and Infectious Diseases (U01AI068618, U19AI089986, UM1AI068618, UM1AI068614, UM1AI068635, U19AI128914)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Gene expression and cancer classification
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Gene Regulatory Network Analysis
Keywords
CITE-seq, peripheral blood mononuclear cells, cell identity, antibody profiling
MeSH terms
COVID-19, SARS-CoV-2, Animals, Cell Line, Humans, Immunity, Leukocytes, Mononuclear, Lymphocytes, Vaccination, 3T3 Cells, Sequence Analysis, RNA, Gene Expression Profiling, Mice, Single-Cell Analysis, Transcriptome

About these checks

Registry, retraction, Cites and statistics results are computed, not judged, with one exception: a printed reference that no index lists is read once by a model before it counts as a citation error. Claims are extracted by a model and kept only when their quote appears verbatim in the text; each cited source is then checked, with a model’s help, for whether it says what the claim relies on. Those verdicts can be wrong — read the claim’s quote and the cited source before relying on one.

A language model drew the claims from the paper on 2026-10-10, under version 3 of the claim definition (model deepseek-v4.1-flash, prompt extract-v3).

A second model labelled each claim on 2026-10-10 (model jev-1.13-20260917, definition revision gold-guideline-v0.2).

Statistics and arithmetic are recomputed, not judged.

Confidence bands are shown only for fields calibrated on our annotated test set. None is calibrated yet, so no band is shown.

What counts as a claim ›How we measure confidence ›How we check this ›

Are you an author, or do you see a problem with this page? Request a review