Cited by 194Single-cell and spatial transcriptomicsGene Regulatory Network AnalysisGene expression and cancer classification

scPerturb: harmonized single-cell perturbation data

Open access

Published version of a preprint

Abstract

Recent biotechnological advances led to growing numbers of single-cell perturbation studies, which reveal molecular and phenotypic responses to large numbers of perturbations. However, analysis across diverse datasets is typically hampered by differences in format, naming conventions, and data filtering. In order to facilitate development and benchmarking of computational methods in systems biology, we collect a set of 44 publicly available single-cell perturbation-response datasets with molecular readouts, including transcriptomics, proteomics and epigenomics. We apply uniform pre-processing and quality control pipelines and harmonize feature annotations. The resulting information resource enables efficient development and testing of computational analysis methods, and facilitates direct comparison and integration across datasets. In addition, we introduce E-statistics for perturbation effect quantification and significance testing, and demonstrate E-distance as a general distance measure for single cell data. Using these datasets, we illustrate the application of E-statistics for quantifying perturbation similarity and efficacy. The data and a package for computing E-statistics is publicly available at scperturb.org. This work provides an information resource and guide for researchers working with single-cell perturbation data, highlights conceptual considerations for new experiments, and makes concrete recommendations for optimal cell counts and read depth.

Integrity signals

Registry

Crossref · 10.1038/s41592-023-02144-y

Cites

78 of 78 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Claims the paper makes

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

  • Findingown data

    The authors assembled 44 publicly available single-cell perturbation-response datasets that include transcriptomic, proteomic, and epigenomic readouts.

    “we collect a set of 44 publicly available single-cell perturbation–response datasets with molecular readouts, including transcriptomics, proteomics and epigenomics.”
    Quote in the paperSubstantial claim
  • Findingown data

    For each scRNA-seq dataset, the authors supply count matrices in which each cell has a perturbation annotation and quality control metrics.

    “For each scRNA-seq dataset we supply count matrices, in which each cell has a perturbation annotation as well as quality control metrics.”
    Quote in the paperSubstantial claim
  • Findingown data

    Molecular readouts for the 44 datasets include transcript, protein, and epigenetic profiles.

    “Molecular readouts for the 44 single-cell perturbation–response datasets include transcript, protein and epigenetic profiles (Table 1 and Fig. 2a).”
    Quote in the paperNot a substantial claim
  • Findingown data

    The authors collected 44 publicly available perturbation-response datasets from 25 publications.

    “We collected 44 publicly available perturbation–response datasets from 25 publications (Table 1 and Supplementary Fig. 1b).”
    Quote in the paperSubstantial claim
  • Findingown data

    Measures of sample quality vary significantly across datasets.

    “The measures of sample quality vary significantly across datasets (Fig. 2b).”
    Quote in the paperSubstantial claim
  • Classified as probably not a substantial claim. 35
  • May combine two claims. 41
  • Probably not a substantial claim: may not pass S3 (load-bearing), S4 (based). May combine two claims.

    Findingnone given

    The total number of cells per dataset is usually restricted by experimental limitations although it has increased over time.

    “The total number of cells per dataset is usually restricted by experimental limitations although it has increased over time.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Findingcitation

    Three CITE-seq datasets have protein and RNA counts that are downloadable separately.

    “Three CITE-seq datasets have protein and RNA counts that are downloadable separately6,9.”
    Quote in the paperNot a substantial claim0 of 2 sources checked
    • cites [6] Frangieh, C. J. et al. Multimodal pooled Perturb-CITE-seq screens in patient models define mechanisms of cancer immune evasion. Nat. Genet. 53, 332–341 (2021). source

      Not checked in the free run.

    • cites [9] Papalexi, E. et al. Characterizing the molecular regulation of inhibitory immune checkpoints with multimodal single-cell screens. Nat. Genet. 53, 322–331 (2021). source

      Not checked in the free run.

  • May combine two claims.

    Findingown data

    32 datasets report scRNA-seq exclusively; the resource also includes scATAC-seq from three publications, one with simultaneous protein measurements.

    “While 32 datasets report scRNA-seq exclusively, we also include scATAC-seq from three publications, one with simultaneous protein measurements39.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [39] Mimitou, E. P. et al. Multiplexed detection of proteins, transcriptomes, clonotypes and CRISPR perturbations in single cells. Nat. Methods 16, 409–412 (2019). source

      Not checked in the free run.

  • May combine two claims.

    Findingown data

    Cells in 32 datasets underwent perturbation using CRISPR, and in 9 datasets the perturbation was done using drugs.

    “Cells in 32 datasets underwent perturbation using CRISPR, and in 9 datasets the perturbation was done using drugs.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingnone given

    The authors describe E-distance and E-test as tools for statistical comparisons of sets of cells and benchmark their robustness and applicability for distinguishing perturbations across datasets and modalities.

    “We also describe E-distance and E-test as tools for statistical comparisons of sets of cells and benchmark their robustness and applicability for distinguishing perturbations across datasets and modalities.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Findingown data

    Five drugs in the sci-Plex 3 dataset also appear in other drug perturbation datasets.

    “Five drugs in this dataset also appear in other drug perturbation datasets (Supplementary Table 4).”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingcitation

    The most extensive drug dataset is sci-Plex 3, which includes 188 drugs tested across three cell lines; 107 of those perturbations had significant effects on cell states according to E-test analysis.

    “The most extensive drug dataset is sci-Plex 3, which includes 188 drugs tested across three cell lines52; 107 of those perturbations had significant effects on cell states according to E-test analysis (Supplementary Table 3).”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [52] Srivatsan, S. R. et al. Massively multiplex chemical transcriptomics at single-cell resolution. Science 367, 45–51 (2020). source

      Not checked in the free run.

  • May combine two claims.

    Findingown data

    Despite the increase in E-distance with falling cell numbers, the number of significant perturbations correctly decreases with fewer cells, and only some datasets have saturated significance at the full number of cells in that dataset.

    “Despite the increase in E-distance with falling cell numbers, the number of significant perturbations correctly decreases with fewer cells, and only some datasets have saturated significance at the full number of cells in that dataset (Fig. 5b).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The number of significant perturbations under the E-test saturates at around 500 UMIs per cell, with most perturbations that were significant at the full measured read depth maintaining that significance even with far fewer counts per cell.

    “The number of significant perturbations under the E-test, however, saturates at around 500 UMIs per cell, with most perturbations that were significant at the full measured read depth maintaining that significance even with far fewer counts per cell (Fig. 5d).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The perturbations in the group more dissimilar to the unperturbed condition have a low E-distance to each other, suggesting a similar phenotype induced by perturbation of these genes (IFNGR1, IFNGR2, JAK2, STAT1).

    “The perturbations in the group more dissimilar to the unperturbed condition have a low E-distance to each other, suggesting a similar phenotype induced by perturbation of these genes (IFNGR1, IFNGR2, JAK2, STAT1).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    Cells in the top three perturbations with the largest E-distance to unperturbed cells are easily distinguishable from the gray unperturbed cells, while cells from the bottom three weakest perturbations are part of a single, uniform cloud visually indistinguishable from the unperturbed cells.

    “Cells in the top three perturbations with the largest E-distance to unperturbed cells are easily distinguishable from the gray unperturbed cells, while cells from the bottom three weakest perturbations are part of a single, uniform cloud visually indistinguishable from the unperturbed cells.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    The three perturbations with highest E-distance are double perturbations while the three with the lowest E-distance are not.

    “Indeed, the three perturbations with highest E-distance are double perturbations while the three with the lowest E-distance are not.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingnone given

    Many datasets supplied only processed data, or raw data were available only after institutional clearance.

    “For all data types, many datasets supplied only processed data, or raw data were available only after institutional clearance.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    This dataset contains >2.5 million cells, is the largest in the database, and has an average number of cells per gene detected significantly higher than in other datasets.

    “Containing >2.5 million cells, this dataset is the largest in our database, with the average number of cells, which each gene is detected in, being significantly higher than in other datasets.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Findingnone given

    The most frequently perturbed gene, MYC, is perturbed in nine datasets from three publications.

    “The most frequently perturbed gene, MYC, is perturbed in nine datasets from three publications.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    In the NormanWeissman2019 dataset, pairwise E-distances between perturbations reveal clusters of similar perturbations, including one of strongly acting CEBPA-targeting perturbations.

    “E-distances between all pairs of perturbations in the dataset NormanWeissman2019. The color scale is clipped at 5% highest and lowest percentiles. Clusters of similar perturbations are visible, for example a cluster of strongly acting perturbations targeting CEBPA at the top.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Study design and conductnone given

    The authors processed the datasets with uniform quality-control pipelines and harmonized feature annotations.

    “We apply uniform quality control pipelines and harmonize feature annotations.”
    Quote in the paperSubstantial claim
  • Classified as probably not a substantial claim. May combine two claims.

    Premisecitation

    Targeted perturbations can be modeled as affecting individual nodes of a regulatory network, which allows investigation of mechanisms and inference of regulatory interactions and their direction.

    “Experiments with targeted perturbations can be modeled as affecting individual nodes of a regulatory network model, enabling investigation of mechanistic processes and inference of regulatory interactions and their directionality11.”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [11] Pratapa, A., Jalihal, A. P., Law, J. N., Bharadwaj, A. & Murali, T. M. Benchmarking algorithms for gene regulatory network inference from single-cell transcriptomic data. Nat. Methods 17, 147–154 (2020). source

      Not checked in the free run.

  • May combine two claims.

    Interpretationreasoning

    The scPerturb resource supports development and testing of computational methods and makes comparison and integration across datasets easier.

    “The resulting information resource, scPerturb, enables development and testing of computational methods, and facilitates comparison and integration across datasets.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationown data

    E-statistics can be used to quantify how similar perturbations are and how effective they are.

    “We illustrate the application of E-statistics for quantifying similarity and efficacy of perturbations.”
    Quote in the paperSubstantial claim
    • Rests on an assumption: There is presently no convention for statistical comparison of response profiles in perturbation studies.
  • May combine two claims.

    Interpretationnone given

    E-statistics can quantify perturbation effects and test their significance.

    “We describe energy statistics (E-statistics) for quantification of perturbation effects and significance testing”
    Quote in the paperSubstantial claim
    • Rests on an assumption: There is presently no convention for statistical comparison of response profiles in perturbation studies.
  • Probably not a substantial claim: may not pass S3 (load-bearing), S4 (based).

    Interpretationnone given

    The work offers recommendations for experimental design, including optimal cell counts and read depth.

    “This work provides an information resource for researchers working with single-cell perturbation data and recommendations for experimental design, including optimal cell counts and read depth.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S2 (checkable).

    Interpretationreasoning

    Statistical comparisons between perturbations and quantification of perturbation strength should ideally use a multivariate distance measure between sets of cells.

    “Ideally, statistical comparisons between perturbations and quantification of perturbation strength should be based on a multivariate distance measure between sets of cells.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Premisecitation

    Cell-based statistical measures can identify successfully perturbed cells but currently do not quantify similarity of perturbation effects.

    “Cell-based statistical measures can be used to identify successfully perturbed cells but do not currently quantify similarity of perturbation effects2,9.”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [2] Dixit, A., Parnas, O., Li, B. & Chen, J. Perturb-seq: dissecting molecular circuits with scalable single-cell RNA profiling of pooled genetic screens. Cell 167, 1853–1866 (2016). source

      Not checked in the free run.

    • cites [9] Papalexi, E. et al. Characterizing the molecular regulation of inhibitory immune checkpoints with multimodal single-cell screens. Nat. Genet. 53, 322–331 (2021). source

      Not checked in the free run.

  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    Several distance measures for scRNA-seq have been explored, including Wasserstein distance, maximum mean discrepancy, neighborhood-based measures, and energy distance.

    “Multiple distance measures for scRNA-seq have been explored by the single-cell community in recent years, including Wasserstein distance computed in an optimal transport framework22,23, maximum mean discrepancy24, neighborhood-based measures17,18, and energy distance (E-distance)20.”
    Quote in the paperNot a substantial claim0 of 5 sources checked
    • cites [22] Chen, W. S. et al. Uncovering axes of variation among single-cell cancer specimens. Nat. Methods 17, 302–310 (2020). source

      Not checked in the free run.

    • cites [23] Schiebinger, G. et al. Optimal-transport analysis of single-cell gene expression identifies developmental trajectories in reprogramming. Cell 176, 928–943 (2019). source

      Not checked in the free run.

    • cites [24] Lotfollahi, M., Naghipourfar, M., Theis, F. J. & Wolf, F. A. Conditional out-of-distribution generation for unpaired data using transfer VAE. Bioinformatics 36, i610–i617 (2020). source

      Not checked in the free run.

    • cites [17] Burkhardt, D. B. et al. Quantifying the effect of experimental perturbations at single-cell resolution. Nat. Biotechnol. 39, 619–629 (2021). source

      Not checked in the free run.

    • cites [18] Dann, E., Henderson, N. C., Teichmann, S. A., Morgan, M. D. & Marioni, J. C. Differential abundance testing on single-cell data using k-nearest neighbor graphs. Nat. Biotechnol. 40, 245–253 (2022). source

      Not checked in the free run.

  • May combine two claims.

    Premisecitation

    Perturbation datasets have generally been too small to capture the complexity of cellular systems, so predictive models of regulatory interactions remain hard to infer.

    “Typically, however, perturbation datasets have been too small to elucidate the complexity of cellular systems; thus, accurately predictive models of regulatory interactions remain difficult to infer12.”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [12] Gross, T. & Blüthgen, N. Identifiability and experimental design in perturbation studies. Bioinformatics 36, i482–i489 (2020). source

      Not checked in the free run.

  • May combine two claims.

    Premisecitation

    A distance measure between sets of cells describes difference or similarity between expression profiles of cells treated with distinct perturbations, which can infer unique or shared mechanisms or identify perturbation targets.

    “Such a distance measure describes the difference or similarity between the expression profiles of cells treated with distinct perturbations, thus inferring unique or shared mechanisms or identifying perturbation targets, which tend to produce similar shifts in molecular profiles20,21.”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [20] Replogle, J. M. et al. Mapping information-rich genotype–phenotype landscapes with genome-scale Perturb-seq. Cell 185, 2559–2575 (2022). source

      Not checked in the free run.

    • cites [21] Tian, R. et al. Genome-wide CRISPRi/a screens in human neurons link lysosomal failure to ferroptosis. Nat. Neurosci. 24, 1020–1034 (2021). source

      Not checked in the free run.

  • May combine two claims.

    Interpretationnone given

    The energy test (E-test) is a statistically reliable tool for computational diagnostics of information content for a specific perturbation and can inform experimental design and data selection for training statistical models of perturbation effects.

    “The associated energy test (E-test) is a statistically reliable tool for computational diagnostics of information content for a specific perturbation, and can inform design of experiments and data selection for training statistical models of perturbation effects.”
    Quote in the paperSubstantial claim
    • Rests on an assumption: There is presently no convention for statistical comparison of response profiles in perturbation studies.
  • Probably not a substantial claim: may not pass S2 (checkable).

    Premisenone given

    Analyzing increasingly large perturbation datasets reliably requires efficient statistical tools that can handle many cells and perturbations.

    “Reliable analysis of increasingly large perturbation datasets requires efficient statistical tools to harness large numbers of cells and perturbations.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    Studies with mixtures of cell types have developed complex methods for quantifying similarity between heterogeneous cell populations.

    “Studies with mixtures of cell types have developed complex methods for quantifying similarity between heterogeneous cell populations17-19.”
    Quote in the paperNot a substantial claim0 of 3 sources checked
    • cites [17] Burkhardt, D. B. et al. Quantifying the effect of experimental perturbations at single-cell resolution. Nat. Biotechnol. 39, 619–629 (2021). source

      Not checked in the free run.

    • cites [18] Dann, E., Henderson, N. C., Teichmann, S. A., Morgan, M. D. & Marioni, J. C. Differential abundance testing on single-cell data using k-nearest neighbor graphs. Nat. Biotechnol. 40, 245–253 (2022). source

      Not checked in the free run.

    • cites [19] Gehring, J., Park, J. H., Chen, S., Thomson, M. & Pachter, L. Highly multiplexed single-cell RNA-seq by DNA oligonucleotide tagging of cellular proteins. Nat. Biotechnol. 38, 35–38 (2020). source

      Not checked in the free run.

  • May combine two claims.

    Premisereasoning

    Some studies combine all cells in a perturbation setting into pseudo-bulk calculations, which loses information about cell-to-cell variation in response.

    “Some studies carry out pseudo-bulk calculations by combining all cells in a given perturbation setting1,7. This means losing information about the cell-to-cell variation in response.”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [1] Datlinger, P. et al. Pooled CRISPR screening with single-cell transcriptome readout. Nat. Methods 14, 297–301 (2017). source

      Not checked in the free run.

    • cites [7] Adamson, B. et al. A multiplexed single-cell CRISPR screening platform enables systematic dissection of the unfolded protein response. Cell 167, 1867–1882 (2016). source

      Not checked in the free run.

  • May combine two claims.

    Premisecitation

    E-distance is a fundamental statistical measure of distances between point clouds that can be used in a statistical test to identify strong or weak perturbations and to distinguish perturbations affecting distinct cellular sub-processes.

    “Here, we exclusively use E-distance, a fundamental statistical measure of distances between point clouds that can be used in a statistical test to identify strong or weak perturbations as well as to distinguish between perturbations affecting distinct cellular sub-processes.”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [20] Replogle, J. M. et al. Mapping information-rich genotype–phenotype landscapes with genome-scale Perturb-seq. Cell 185, 2559–2575 (2022). source

      Not checked in the free run.

    • Rests on an assumption: There is presently no convention for statistical comparison of response profiles in perturbation studies.
  • Probably not a substantial claim: may not pass S2 (checkable).

    Premisereasoning

    Moving from single-dataset to multi-dataset analysis requires developing quantitative approaches to perturbation biology.

    “Moving from single-dataset to multi-dataset analysis requires development of quantitative approaches to perturbation biology.”
    Quote in the paperNot a substantial claim
  • Classified as probably not a substantial claim.

    Interpretationown data

    The quantification of perturbation strength and comparison of experiment-specific variables may serve as a reference for optimal design of future single-cell experiments.

    “Our quantification of perturbation strength and comparison of experiment-specific variables, such as the number of perturbations and the number of cells per perturbation, may serve as a reference for optimal design of future single-cell experiments.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Premisecitation

    Increasing sequencing depth increases UMI counts measured even for low-expression genes, reducing uncertainty associated with zero counts.

    “Increasing sequencing depth increases the UMI counts measured even for low-expression genes, reducing the uncertainty associated with zero counts34,47.”
    Quote in the paperNot a substantial claim0 of 2 sources checked
    • cites [34] Svensson, V., da Veiga Beltrame, E. & Pachter, L. A curated database reveals trends in single-cell transcriptomics. Database 2020, baaa073 (2020). source

      Not checked in the free run.

    • cites [47] Haque, A., Engel, J., Teichmann, S. A. & Lönnberg, T. A practical guide to single-cell RNA-sequencing for biomedical research and clinical applications. Genome Med. 9, 75 (2017). source

      Not checked in the free run.

  • May combine two claims.

    Premisecitation

    Unified frameworks for accessing single-cell data are in active development but do not currently support perturbation datasets or standardize perturbation annotations.

    “Unified frameworks for accessing single-cell data are in active development, but do not currently support perturbation datasets or standardize perturbation annotations37,38.”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [37] Fischer, D. S. et al. Sfaira accelerates data and model reuse in single cell genomics. Genome Biol. 22, 248 (2021). source

      Not checked in the free run.

    • cites [38] Chan Zuckerberg CELLxGENE Discover. Cellxgene Data Portal. https://cellxgene.cziscience.com/

      Not checked in the free run.

  • Classified as probably not a substantial claim.

    Premisecitation

    E-distance is a statistical measure of the distance between two distributions that provides intuition about the signal-to-noise ratio in a dataset.

    “To compare and evaluate the effect and strength of perturbations in each dataset we used the E-distance, a statistical measure of the distance between two distributions, which provides intuition about the signal-to-noise ratio in a dataset48.”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [48] Székely, G. J. & Rizzo, M. L. Energy statistics: a class of statistics based on distances. J. Stat. Plan. Inference 143, 1249–1272 (2013). source

      Not checked in the free run.

    • Rests on an assumption: There is presently no convention for statistical comparison of response profiles in perturbation studies.
  • May combine two claims.

    Premisecitation

    Existing collections of datasets are mainly for filtering and do not provide a unified format for perturbations.

    “Existing collections of datasets are primarily a means for filtering and do not supply a unified format for perturbations33-35.”
    Quote in the paperSubstantial claim0 of 3 sources checked
    • cites [33] Lance, C. et al. Multimodal single cell data integration challenge: results and lessons learned. Preprint at https://doi.org/10.1101/2022.04.11.487796 (2022). source

      Not checked in the free run.

    • cites [34] Svensson, V., da Veiga Beltrame, E. & Pachter, L. A curated database reveals trends in single-cell transcriptomics. Database 2020, baaa073 (2020). source

      Not checked in the free run.

    • cites [35] Broad Institute. Single Cell Portal. https://singlecell.broadinstitute.org/single_cell (2022).

      Not checked in the free run.

  • May combine two claims.

    Premisecitation

    Although several large databases of perturbations with bulk readouts exist, single-cell perturbation technologies are newer and their data are not unified.

    “Although several large databases of perturbations with bulk readouts exist, single-cell perturbation technologies are newer and the data are not unified31,32.”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [31] Stathias, V. et al. LINCS Data Portal 2.0: next generation access point for perturbation–response signatures. Nucleic Acids Res. 48, D431–D439 (2020). source

      Not checked in the free run.

    • cites [32] Tsherniak, A. et al. Defining a cancer dependency map. Cell 170, 564–576 (2017). source

      Not checked in the free run.

  • May combine two claims.

    Premisereasoning

    The exact value of E-distance depends on dataset-specific parameters such as sequencing depth and how cells are distributed; the E-test accounts for these differences by creating a null distribution using permutations of the data.

    “The exact value of the E-distance depends on dataset-specific parameters such as sequencing depth (Supplemental Note Section 4) and how exactly the cells are distributed; the E-test accounts for these differences by creating a null distribution using permutations of the data.”
    Quote in the paperSubstantial claim
  • Classified as probably not a substantial claim.

    Interpretationown data

    As a baseline for significant perturbations as defined by the E-test, the authors suggest at least 200-500 cells per perturbation and an average of 1,000 UMIs per cell as an experimental guideline for distinguishable perturbations.

    “As a baseline for significant perturbations, as defined by the E-test, we suggest at least 200–500 cells per perturbation (Fig. 5b) and an average of 1,000 UMIs per cell (Fig. 5d) as an experimental guideline for distinguishable perturbations.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationreasoning

    A low E-distance indicates that a perturbation did not induce a large shift in expression profiles, reflecting technical problems, weak perturbation effect, or resistance to perturbation.

    “A low E-distance indicates that a perturbation did not induce a large shift in expression profiles, reflecting technical problems in the experiment, weak effect of the perturbation, or resistance to perturbation.”
    Quote in the paperSubstantial claim
  • Classified as probably not a substantial claim.

    Premisereasoning

    For two groups of cells, E-distance relates the distribution of distances between groups (signal) to the distribution of distances within each group (noise).

    “For two groups of cells, it relates the distribution of distances of cells between the groups (‘signal’), to the distribution of distances within each group ('noise').”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationreasoning

    Any disruption of a gene in this cascade leads to a similar outcome and therefore a similar transcriptome profile.

    “In agreement with this, we observe that any disruption of a gene in this cascade leads to a similar outcome and, therefore, a similar transcriptome profile.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationreasoning

    The saturation point depends on perturbation strength and dataset heterogeneity: if all cells are similar, a small set of cells will sufficiently describe every possible response to a perturbation.

    “This saturation point depends on perturbation strength and on dataset heterogeneity: if all cells are similar, a small set of cells will sufficiently describe every possible response to a perturbation.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Interpretationcitation

    The optimal lower bounds on UMI and cell counts for a given experiment depend on downstream specific modeling tasks.

    “The optimal lower bounds on UMI and cell counts for a given experiment depend on downstream specific modeling tasks, as discussed in more detail elsewhere12.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [12] Gross, T. & Blüthgen, N. Identifiability and experimental design in perturbation studies. Bioinformatics 36, i482–i489 (2020). source

      Not checked in the free run.

  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    Another drug dataset applies combinations of three drug perturbations at varying concentrations across samples.

    “Another drug dataset applies combinations of three drug perturbations at varying concentrations across samples19.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [19] Gehring, J., Park, J. H., Chen, S., Thomson, M. & Pachter, L. Highly multiplexed single-cell RNA-seq by DNA oligonucleotide tagging of cellular proteins. Nat. Biotechnol. 38, 35–38 (2020). source

      Not checked in the free run.

  • Classified as probably not a substantial claim. May combine two claims.

    Interpretationreasoning

    Starting the analysis from sequencing reads may have improved the interoperability of datasets in this resource, but guide assignment procedures and demultiplexing algorithms are experimental set-up specific.

    “Starting the analysis from sequencing reads may have improved the interoperability of datasets in this resource, but guide assignment procedures and demultiplexing algorithms are experimental set-up specific.”
    Quote in the paperSubstantial claim
  • Classified as probably not a substantial claim.

    Interpretationreasoning

    Different experimental designs result in various covariates that may need to be accounted for in statistical testing.

    “However, different experimental designs result in various covariates that may need to be accounted for in statistical testing.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Premisecitation

    The sharing of spatial data is hampered by large file sizes and non-standardized analysis pipelines.

    “The sharing of spatial data is hampered by large file sizes and non-standardized analysis pipelines58.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [58] Bredikhin, D., Kats, I. & Stegle, O. MUON: multimodal omics analysis framework. Genome Biol. 23, 42 (2022). source

      Not checked in the free run.

  • Probably not a substantial claim: may not pass S2 (checkable), S3 (load-bearing).

    Interpretationreasoning

    Adding more datasets to this resource, or creating similar resources in the future, would be easier if there were standard formats for sharing perturbation data and single-cell annotations.

    “Adding more datasets to this resource, or the creation of similar resources in the future, would be easier if there were standard formats for sharing perturbation data, and, more generally, standard formats for sharing single-cell annotations.”
    Quote in the paperNot a substantial claim
    • Rests on an assumption: Lack of standardization in data sharing and processing across different research groups and studies hampered the creation of this resource.
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Interpretationreasoning

    Covariates arising from different experimental designs can affect the E-distance, and the authors recommend balancing cell numbers across experimental conditions if in doubt.

    “In general, covariates arising from different experimental designs can affect the E-distance, and we recommend balancing cell numbers across experimental conditions if in doubt (Supplemental Note Section 8).”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Premisenone given

    Although many processed datasets were available on GEO, there is no standard format for sharing CRISPR barcode assignments and other metadata.

    “Although many processed datasets were available on the National Center for Biotechnology Information (NCBI) Gene Expression Omnibus (GEO)55, there is no standard format for sharing CRISPR barcode assignments and other metadata.”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [55] Barrett, T. et al. NCBI GEO: archive for functional genomics data sets – update. Nucleic Acids Res. 41, D991–D995 (2013). source

      Not checked in the free run.

  • Probably not a substantial claim: may not pass S2 (checkable), S3 (load-bearing).

    Interpretationcitation

    For scATAC feature assignments specific to CRISPR perturbations, known loci of action could be used to improve feature calls.

    “In particular, for scATAC feature assignments specific to CRISPR perturbations, known loci of action could be used to improve feature calls40.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [40] Chen, H. et al. Assessment of computational methods for the analysis of single-cell ATAC-seq data. Genome Biol. 20, 241 (2019). source

      Not checked in the free run.

  • May combine two claims.

    Interpretationreasoning

    The uniform annotations in this resource enable data integration and benchmarking as well as exploration of the effects of shared perturbations across datasets.

    “The uniform annotations in this resource enable data integration and benchmarking as well as exploration of the effects of shared perturbations across datasets.”
    Quote in the paperSubstantial claim
    • Rests on an assumption: Unified datasets are key to developing generalizable machine learning methods and establishing multimodal data integration.
  • Probably not a substantial claim: may not pass S2 (checkable), S3 (load-bearing), S4 (based).

    Interpretationnone given

    Such a standardization effort would enable simpler addition of datasets to the resource.

    “Such a standardization effort would enable simpler addition of datasets to the resource.”
    Quote in the paperNot a substantial claim
    • Rests on an assumption: Lack of standardization in data sharing and processing across different research groups and studies hampered the creation of this resource.
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    There is increasing interest in spatial techniques, including single-cell resolution spatial transcriptomics, epigenomics and proteomics.

    “Moreover, there is increasing interest in spatial techniques, including single-cell resolution spatial transcriptomics, epigenomics and proteomics57.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [57] Tian, L., Chen, F. & Macosko, E. Z. The expanding vistas of spatial transcriptomics. Nat. Biotechnol. 41, 773–782 (2023). source

      Not checked in the free run.

  • Probably not a substantial claim: may not pass S2 (checkable), S3 (load-bearing), S4 (based).

    Interpretationnone given

    The unified annotations and perturbation significance testing should prove especially useful to the machine learning community for training computational models on this data.

    “The unified annotations and perturbation significance testing should prove especially useful to the machine learning community for training computational models on this data.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Study design and conductown data

    Within each dataset, p-values from the E-test were adjusted for multiple comparisons with the Holm-Sidak procedure.

    “We corrected for multiple testing using the Holm–Sidak method per dataset.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S4 (based).

    Interpretationnone given

    It is difficult to ascertain to what extent low E-distances between perturbed and unperturbed cells are caused by technical noise or simply weak perturbation effects compared with biological variation in groups of cells before and after perturbation.

    “Unfortunately, it is difficult to ascertain to what extent low E-distances between perturbed and unperturbed cells are caused by technical noise or simply weak perturbation effects compared with biological variation in groups of cells before and after perturbation.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Study design and conductown data

    For the E-distance analysis, cells with at least 1,000 UMIs and genes detected in at least 50 cells were kept, 2,000 highly variable genes were chosen, and counts were normalized and log-transformed without z-scaling.

    “Cells were retained if they had a minimum of 1,000 UMIs, and genes, with a minimum of 50 cells. A total of 2,000 highly variable genes were selected using scanpy.pp.find_variable_genes with flavor seurat_v3. We normalized the count matrix using scanpy.pp.normalize_total and log-transformed the data using scanpy.pp.log1p; we did not z-scale the data.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Study design and conductown data

    For the Spear-ATAC and CRISPR_sciATAC datasets, a barcode was assigned the most abundant sgRNA only when that guide's count was above 20 and it made up at least 80% of all guide counts.

    “We assigned the sgRNA with the highest counts to a cell barcode if the sgRNA count exceeded 20 and if that sgRNA comprised at least 80% of all sgRNA counts.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Interpretationnone given

    Experiment design choices, such as the optimal number of cells per perturbation and the sequencing depth for each cell, depend on the questions that a particular perturbation study is intended to answer, and on the strength and uniqueness of the gene expression changes caused by the perturbations.

    “Experiment design choices, such as the optimal number of cells per perturbation and the sequencing depth for each cell, depend on the questions that a particular perturbation study is intended to answer, and on the strength and uniqueness of the gene expression changes caused by the perturbations.”
    Quote in the paperNot a substantial claim
  • Classified as probably not a substantial claim. May combine two claims.

    Interpretationown data

    For perturbation distinguishability as defined by the E-test, regardless of experiment parameters, the authors recommend measuring at least 200-500 cells per perturbation and an average of 1,000 UMIs per cell.

    “For perturbation distinguishability as defined by the E-test, regardless of experiment parameters, we recommend measuring at least 200–500 cells per perturbation and an average of 1,000 UMIs per cell.”
    Quote in the paperSubstantial claim
  • Classified as probably not a substantial claim. May combine two claims.

    Study design and conductreasoning

    Cells with under 1,000 or over 100,000 mapped fragments were removed, and a minimum transcription start site enrichment score of 4 was required, to keep quality consistent and signal-to-noise sufficient across datasets.

    “To ensure a consistent and homogenous quality throughout the different datasets, we filtered out cells with fewer than 1,000 and more than 100,000 mapped fragments. Furthermore, we required a minimum transcription start site enrichment score of 4 to ensure a sufficient signal-to-noise ratio.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Premisecitation

    The scATAC-seq features used here were originally developed for non-perturbed cells, for purposes such as identifying cell types, discovering cell-type-specific regulatory elements, and reconstructing differentiation trajectories.

    “Note that these features were originally developed for scATAC-seq data on non-perturbed cells, with goals such as the identification of cell type, discovery of cell type-specific regulatory elements, or reconstruction of cellular differentiation trajectories43,61.”
    Quote in the paperNot a substantial claim0 of 2 sources checked
    • cites [43] Buenrostro, J. D. et al. Single-cell chromatin accessibility reveals principles of regulatory variation. Nature 523, 486–490 (2015). source

      Not checked in the free run.

    • cites [61] Satpathy, A. T. et al. Massively parallel single-cell chromatin landscapes of human immune cell development and intratumoral T cell exhaustion. Nat. Biotechnol. 37, 925–936 (2019). source

      Not checked in the free run.

  • Classified as probably not a substantial claim. May combine two claims.

    Study design and conductreasoning

    To avoid bias from unequal cell numbers per perturbation, each dataset was subsampled so that all perturbations had equal cell counts, after perturbations with fewer than 50 cells were removed.

    “To avoid problems due to different numbers of cells per perturbation, we subsampled each dataset such that all perturbations had the same number of cells. We removed all perturbations with fewer than 50 cells and then subsampled to the number of cells in the smallest perturbation left after filtering.”
    Quote in the paperSubstantial claim
67 of 67
  • The paper relies on this without support.

    There is presently no convention for statistical comparison of response profiles in perturbation studies.

    “There is presently no convention for statistical comparison of response profiles in perturbation studies.”
  • The paper relies on this without support.

    Unified datasets are key to developing generalizable machine learning methods and establishing multimodal data integration.

    “However, unified datasets are key to developing generalizable machine learning methods and establishing multimodal data integration.”
  • The paper relies on this without support.

    There is no consensus representation of scATAC-seq data.

    “In contrast to scRNA-seq data, which can be represented naturally as counts per gene, there is no consensus representation of scATAC-seq data.”
  • The paper relies on this without support.

    Lack of standardization in data sharing and processing across different research groups and studies hampered the creation of this resource.

    “Lack of standardization in data sharing and processing across different research groups and studies hampered the creation of this resource.”
  • The paper relies on this without support.

    No established best practices yet exist for the analysis of single-cell perturbation data.

    “To our knowledge, best practices have not yet been established for analysis of single-cell perturbation data.”
5 of 5

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?

  • Combinatorial multiplexing experiments remain prohibitively expensive and technically challenging to scale.

    “Furthermore, combinatorial multiplexing experiments remain prohibitively expensive and technically challenging to scale (7, 8).”

    From A systematic comparison of single-cell perturbation response prediction models (2026) · cites this paper as [8]

    The Abstract says perturbation datasets have typically been too small, and the Discussion describes technical complications in high-multiplicity screens. Those passages support a general scaling challenge, but the paper does not say that combinatorial multiplexing experiments are prohibitively expensive or that their experimental scaling is technically prohibitive. The Discussion instead anticipates that technological advances will decrease experimental costs. (LLM, conf=0.77)

    Experimental — accuracy not yet measured
  • Variability from experimental conditions, including transduction efficiency, RNA sequencing depth, reagent chemistry, and timing of collection, further dilutes the true perturbation-derived signal.

    “Variability from experimental conditions, including transduction efficiency, RNA sequencing depth, reagent chemistry, and timing of collection, further dilutes the true perturbation-derived signal.”

    From Predicting cellular responses to perturbation across diverse contexts with State (2026) · cites this paper as [33]

    The Results (page not provided) say sequencing depth and experiment quality affect the distinguishability of perturbations. The Discussion (page not provided) says changing the time between perturbation and cell collection may improve the signal-to-noise ratio. The cited passages do not establish transduction efficiency or reagent chemistry as contributors to the claimed signal dilution. (LLM, conf=0.86)

    Experimental — accuracy not yet measured
  • Existing methods for capturing aspects of the DE effect distribution have conceptual limitations, including lack of clear estimand, ambiguous units, or finite-sample bias.

    “However, while these approaches are useful, they have conceptual limitations. Some lack a clear estimand (such as the running-sum GSEA statistic); some have ambiguous units (such as the energy distance); some are biased at finite sample size (such as the proportion of DE genes π in iDEA, see O’Connor et al11).”
  • In Papalexi2021 and NormanWeissman2019 Perturb-seq datasets, permutations in molecular targets that are second messages to the same pathway elicit paradoxically similar phenotypes (small E-distance).

    “This type of analysis was conducted on the Papalexi2021 and the NormanWeissman2019 Perturb-seq datasets, revealing that permutations in molecular targets that are second messages to the same pathway elicit paradoxically similar phenotypes (i.e. small E-distance).”
  • scPerturb uses the E-distance metric to reveal similarities between two or more perturbations, facilitating studies for potential synergistic treatments.

    “Interestingly, the scPerturb [42] framework (referenced in the manuscript's perturbational data section) also uses the E-distance metric to reveal similarities between two or more perturbations, facilitating studies for potential synergistic treatments.”

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.

  • Conflicting distance metrics for E-distance calculations

    Major

    The data-analysis section says E-distances were calculated with Euclidean distance, while the E-distance definition says cell-wise distances were squared Euclidean. These produce different values for the paper's central statistic.

1 of 1

Cites

Flagged first, then by citations

Cited by

Sorted by citations

Authors
Peidli, S., Green, T. D., Shen, C., Gross, T., Min, J., Garda, S., Yuan, B., Schumacher, L. J., Taylor-King, J. P., Marks, D. S., Luna, A., Blüthgen, N., & Sander, C.
Affiliations
  • Humboldt-Universität zu Berlin (DE)
  • Harvard University (US)
  • Broad Institute (US)
  • MRC Centre for Regenerative Medicine (GB)
Corresponding author
  • Peidli, Stefan
  • Luna, Augustin
  • Blüthgen, Nils
  • Sander, Chris
Venue
Nature Methods
Published
2024-01-26
Volume
21
Issue
3
Pages
531-540
ISSN
  • 1548-7091 (ISSN-L)
  • 1548-7105
OpenAlex
W4391262585
Cited by
194
References
78
Citations by year
2026: 84 · 2025: 63 · 2024: 38 · 2023: 8
FWCI
24.92
Type
journal article
Language
en
Open access
Locations
  • Nature Methods (publishedVersion)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion, other-oa)
  • PubMed (publishedVersion)
  • PubMed Central (acceptedVersion, other-oa)
Funders
  • Einstein Stiftung Berlin
  • Wellcome Leap ∆Tissue Program National Resource for Network Biology
  • Deutsche Forschungsgemeinschaft (RTG2424)
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 Regulatory Network Analysis
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Gene expression and cancer classification
Keywords
single-cell perturbation, data harmonization, data interoperability, transcriptomics, proteomics, epigenomics, quality control, experimental design, cell number, sequencing depth
MeSH terms
Software, Gene Expression Profiling, Proteomics, Epigenomics, Single-Cell Analysis

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.

Claims read by deepseek-v4.1-flash under definition version 3 (extract-v3) on 2026-10-05.

Classified by jev-1.13-20260917 (gold-guideline-v0.2) on 2026-10-05.

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