Cited by 79CRISPR and Genetic EngineeringSingle-cell and spatial transcriptomicsChromatin Remodeling and Cancer

Efficient combinatorial targeting of RNA transcripts in single cells with Cas13 RNA Perturb-seq

Open access
Wessels, Hans-Hermann, Méndez-Mancilla, Alejandro, Hao, Yuhan, Papalexi, Efthymia, Mauck, William M., Lu, Lu, Morris, John A., Mimitou, Eleni P. +3 more·2023·Nature Methods

Published version of a preprint

Abstract

This work introduces CaRPool-seq that utilizes the RNA-targeting CRISPR–Cas13d system to perform combinatorial perturbations in single-cell screens. Pooled CRISPR screens coupled with single-cell RNA-sequencing have enabled systematic interrogation of gene function and regulatory networks. Here, we introduce Cas13 RNA Perturb-seq (CaRPool-seq), which leverages the RNA-targeting CRISPR–Cas13d system and enables efficient combinatorial perturbations alongside multimodal single-cell profiling. CaRPool-seq encodes multiple perturbations on a cleavable CRISPR array that is associated with a detectable barcode sequence, allowing for the simultaneous targeting of multiple genes. We compared CaRPool-seq to existing Cas9-based methods, highlighting its unique strength to efficiently profile combinatorially perturbed cells. Finally, we apply CaRPool-seq to perform multiplexed combinatorial perturbations of myeloid differentiation regulators in an acute myeloid leukemia (AML) model system and identify extensive interactions between different chromatin regulators that can enhance or suppress AML differentiation phenotypes.

Integrity signals

Registry

Crossref · 10.1038/s41592-022-01705-x

Cites

57 of 57 references verified0 citation errors · 0 not checkable

Statistics

Nothing to recompute.

Claims the paper makes

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

  • Findingown data

    In singlet cells, RNA and bcgRNA labels showed 99.2% concordance.

    “Moreover, we observed extremely high concordance between RNA and bcgRNA labels in singlet cells (99.2%) (Fig. 1d,e).”
    Quote in the paperSubstantial claim
  • Findingown data

    Transcriptomic pseudobulk profiles for all 26 targeting gRNA groups showed decreased mRNA expression for each targeted transcript, including when three genes were perturbed simultaneously.

    “When examining transcriptomic pseudobulk profiles for all 26 targeting gRNA groups, we observed decreased messenger RNA expression for each targeted transcript, even when perturbing transcripts of three genes simultaneously (15 examples in Extended Data Fig. 3f,g).”
    Quote in the paperSubstantial claim
  • Findingown data

    Indirect guide capture with an optimized configuration gave the strongest crRNA transcript detection among the tested methods, for both catalytically active and inactive Cas13 proteins.

    “While all methods successfully induced robust knockdown (Fig. 1b), we found that indirect guide capture with an optimized configuration resulted in the strongest crRNA transcript detection ability (configuration X, Fig. 1c), for both catalytically active and inactive Cas13 proteins.”
    Quote in the paperSubstantial claim
  • Findingown data

    Knockdown strength was similar for multi-gRNA crRNA arrays compared with single-gRNA perturbations.

    “Moreover, the strength of knockdown was similar for multi-gRNA crRNA arrays relative to single gRNA perturbations (Fig. 2b and Extended Data Fig. 3c–e).”
    Quote in the paperSubstantial claim
  • Findingown data

    Both synergistic and dampening effects were observed among dual perturbations.

    “We also observed both synergistic and dampening effects.”
    Quote in the paperSubstantial claim
    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
  • May combine two claims. 39
  • The quote only partly states this claim. 1
  • Classified as probably not a substantial claim. 25
  • May combine two claims.

    Findingown data

    In an AML model system, CaRPool-seq multiplexed combinatorial perturbations of myeloid differentiation regulators and found extensive interactions among chromatin regulators that can enhance or suppress AML differentiation phenotypes.

    “Finally, we apply CaRPool-seq to perform multiplexed combinatorial perturbations of myeloid differentiation regulators in an acute myeloid leukemia (AML) model system and identify extensive interactions between different chromatin regulators that can enhance or suppress AML differentiation phenotypes.”
    Quote in the paperSubstantial claim
    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
  • The quote only partly states this claim.

    Findingown data

    Targeting CD46, CD55, or CD71 with Cas13 produced an average 76.5% (±5.7%) reduction in protein levels for each targeted gene.

    “We observed, on average, a 76.5% (±5.7%) mean reduction in protein levels for each targeted gene after perturbation with Cas13 demonstrating clear evidence of robust molecular perturbation (Fig. 2b–d).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    In a species-mixing experiment of 2,387 cells, 78.5% expressed a single bcgRNA, 1.1% expressed more than one bcgRNA, and 20.4% had no detected bcgRNA.

    “Of 2,387 cells, we found that 78.5% expressed a single bcgRNA (1.1% >1 bcgRNAs; 20.4% no detected bcgRNA).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    In the CITE-seq experiment, 9,355 single-cell profiles were obtained; 74.7% expressed a single bcgRNA and 80.8% expressed at least one bcgRNA.

    “We obtained 9,355 single-cell profiles and demultiplexed them into groups based on the detected bcgRNA (Extended Data Fig. 3b; 74.7% expressed a single bcgRNA, 80.8% expressed at least one bcgRNA).”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    Dually perturbed EP300/KDM1A cells had higher expression of progenitor genes (AZU1) and reduced expression of differentiated marker genes (S100A4) compared with KDM1A perturbation alone.

    “Dually perturbed cells exhibited higher expression of progenitor genes (that is the progenitor marker AZU1), and reduced expression of differentiated marker genes (that is, myeloid marker S100A4) compared to individual KDM1A perturbation.”
    Quote in the paperSubstantial claim
    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
  • May combine two claims.

    Findingown data

    CaRPool-seq and Perturb-seq cells clustered together and grouped by gRNA identity based on ADT levels, indicating similar strength of phenotypic protein perturbation across methods.

    “CaRPool-seq and Perturb-seq cells clustered together (Fig. 3d), and grouped by gRNA identity (Fig. 3e and Extended Data Fig. 5e), again demonstrating that the strength of phenotypic protein perturbation was similar across all methods.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    EP300 knockdown (log2FC 1.57) upregulated CD11b, but dual perturbation with KDM1A (log2FC 2.05) was weaker than KDM1A knockdown alone.

    “In contrast, while individual knockdown of EP300 also leads to CD11b upregulation (log2FC 1.57), dual perturbation with KDM1A (log2FC 2.05) was weaker than the individual KDM1A knockdown.”
    Quote in the paperNot a substantial claim
    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
    • Rests on an assumption: CD14 and CD11b are immunophenotypic markers of monocyte differentiation.
  • May combine two claims.

    Findingown data

    For 33 gene pairs in cluster 1, each individual gene contributed equally to the dual perturbation response and the linear model fit strongly.

    “For 33 gene pairs in cluster 1, we saw that each individual gene’s profile contributed equally to the dual perturbation response and the linear model exhibited a strong fit.”
    Quote in the paperNot a substantial claim
    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Findingnone given

    Alternative polymerases KAPA and Q5 can lead to stronger bias in relative array abundance during library cloning.

    “During library cloning, we noticed two critical details: alternative polymerase KAPA and Q5 can lead to a stronger bias in relative array abundance.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    Cas13d-mediated gene knockdown did not introduce notable off-target effects or deleteriously alter cellular fitness.

    “We performed several analyses to demonstrate that Cas13d-mediated gene knockdown does not introduce notable off-target effects or deleteriously alter cellular fitness.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    Two stabilizing RNA elements (Zika virus-derived xrRNA1 dumbbell and evopreQ1 pseudoknot) led to a robust increase in bcgRNA UMI counts, while protein knockdown was indistinguishable among six structures.

    “While protein knockdown was indistinguishable for all six structures (Fig. 5b,c and Extended Data Fig. 9b), we found that two elements (Zika virus-derived xrRNA1 dumbbell33; evopreQ1 pseudoknot34) led to a robust increase in bcgRNA unique molecular identifier (UMI) counts (Fig. 5d).”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [33] Akiyama, B. M. et al. Zika virus produces noncoding RNAs using a multi-pseudoknot structure that confounds a cellular exonuclease. Science 354, 1148–1152 (2016). source

      Not checked in the free run.

    • cites [34] Anzalone, A. V., Lin, A. J., Zairis, S., Rabadan, R. & Cornish, V. W. Reprogramming eukaryotic translation with ligand-responsive synthetic RNA switches. Nat. Methods 13, 453–458 (2016). source

      Not checked in the free run.

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

    Findingnone given

    Reducing PCR cycles while increasing oligo pool input amounts can decrease bcgRNA reassortment.

    “Further, reducing the number of PCR cycles with increased oligo pool input amounts can decrease bcgRNA reassortment.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    KDM1A knockdown (log2FC 2.41) and HDAC3 knockdown (log2FC 0.53) individually upregulated CD11b, and dual perturbation produced a synergistic effect (log2FC 2.85).

    “For example, individual knockdown of the histone demethylase KDM1A (log2fold change (FC) 2.41) and the histone deacetylase HDAC3 (log2FC 0.53) lead to strong and weak CD11b upregulation, respectively, but dual perturbation led to a synergistic effect (log2FC 2.85).”
    Quote in the paperSubstantial claim
    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
    • Rests on an assumption: CD14 and CD11b are immunophenotypic markers of monocyte differentiation.
  • May combine two claims.

    Findingown data

    Dual perturbation of HDAC3 enhanced the KDM1A transcriptional response signature, consistent with immunophenotypic results.

    “Dual perturbation of HDAC3 enhanced the KDM1A transcriptional response signature (Fig. 4h), consistent with our previously described immunophenotypic results for these cells (Fig. 4c).”
    Quote in the paperSubstantial claim
    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
  • May combine two claims.

    Findingown data

    Each single-gene perturbation affected CD11b expression, and observed log2 fold changes were strongly correlated with gRNA enrichment from bulk CRISPR screens (R = 0.86).

    “As expected, we found that each single-gene perturbation affected CD11b expression, with observed log2 fold changes that were in strong agreement with the level of gRNA enrichment from bulk CRISPR screens (R = 0.86; Fig. 4b, Extended Data Fig. 7c–e and Supplementary Table 5).”
    Quote in the paperSubstantial claim
    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
    • Rests on an assumption: CD14 and CD11b are immunophenotypic markers of monocyte differentiation.
  • May combine two claims.

    Findingreasoning

    The average transcriptomic knockdown was 65% (s.d. 8.7%), consistently lower than the observed protein reduction.

    “The average strength of transcriptomic knockdown (mean 65%, s.d. 8.7%) was consistently reduced compared to the observed protein reduction.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Findingown data

    Additional synergistic combinations HDAC3/GFI and HDAC3/GSE1 enhanced immunophenotypic marker expression and transcriptional differentiation state.

    “Moreover, we identified additional synergistic combinations between HDAC3/GFI and HDAC3/GSE1, both of which enhanced the expression of immunophenotypic markers (Fig. 4c), as well as transcriptional differentiation state (Fig. 4j).”
    Quote in the paperSubstantial claim
    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
    • Rests on an assumption: CD14 and CD11b are immunophenotypic markers of monocyte differentiation.
  • May combine two claims.

    Findingown data

    For the 158 dual-gene perturbations, effects were strongly correlated with the corresponding single perturbations; dual perturbation was typically stronger than the average of individual knockdowns but weaker than the product.

    “We next compared the observed effects of the 158 dual gene perturbations to the effects resulting from the two corresponding single perturbations. We observed a strong correlation and found that the dual perturbation was typically stronger than the average of individual knockdowns, but weaker than the product (Fig. 4c).”
    Quote in the paperSubstantial claim
    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
  • May combine two claims.

    Findingown data

    In a bulk Cas13d CRISPR screen of 439 genes, 26 genes were selected that consistently influenced differentiation across multiple independent gRNAs and had been identified in orthogonal pooled Cas9 screens.

    “By comparing gRNA representation between low and high-expressing bins, we selected 26 genes that influenced differentiation consistently across multiple independent gRNAs and that have also been identified previously in orthogonal pooled Cas9 screens23 (Extended Data Fig. 6a–h and Supplementary Table 4).”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [23] Wang, E. et al. Surface antigen-guided CRISPR screens identify regulators of myeloid leukemia differentiation. Cell Stem Cell 28, 718–731 (2021). source

      Not checked in the free run.

    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
  • May combine two claims.

    Findingown data

    In an integrated trajectory with hematopoietic reference data, NT control cells localized to early points, while single perturbations (KDM1A, GFI1, GSE1) pushed cells further down the differentiation trajectory.

    “We found that NT control cells localized to early points, but that single perturbations (that is, KDM1A, GFI1, GSE1) pushed cells further down the differentiation trajectory (Fig. 4f), consistent with their role in enhancing leukemic differentiation.”
    Quote in the paperSubstantial claim
    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Findingown data

    The THP1 CaRPool-seq experiment detected 52,496 single cells after filtering.

    “For the THP1 experiment we detected 52,496 single cells (nFeature_RNA > 1,000, nFeature_RNA < 8,000, percent.mt < 20) after HTO demultiplexing using HTOdemux as described above.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Findingown data

    The direct capture Perturb-seq analysis detected 16,842 cells.

    “We detected 16,842 cells.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    Direct capture constructs with CS1 or CS2 enable strong knockdown for CD46, but reduced knockdown for CD55 and CD71.

    “Direct capture constructs with CS1 or CS2 enable strong knockdown for CD46, but reduced knockdown for CD55 and CD71.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Findingown data

    The final bcgRNA amplicon is 203 bp and can be sequenced with standard Illumina sequencing primers.

    “The final bcgRNA amplicon (203 bp) can be sequenced with standard Illumina sequencing primers (≥28 cycles read1 and ≥15 cycles read2) (Extended Data Figs. 2 and 3a).”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingown data

    Target gene enrichments for CD14 upregulation correlated best with CD14 target gene enrichments from Wang et al., and correlations were similarly high for the top-ranked 4-5 gRNAs.

    “This analysis indicated that, as expected, target gene enrichments for CD14 upregulation correlated best with CD14 target gene enrichments in Wang et al., and that correlations were similarly high for the top-ranked 4-5 gRNA.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites Wang et al.

      Not checked in the free run.

  • May combine two claims.

    Findingown data

    The heatmaps show a range of interactions between gene pairs, including equal contributions from both genes and dominance by one gene's perturbation signature.

    “These heatmaps visualize a range of observed interactions between gene pairs, including cases where genes contribute equally to the dual perturbation response (d-f), and where one gene’s perturbation signature dominates over the other (g-h).”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Findingown data

    31,308 cells were assigned a single bcgRNA.

    “We assigned 31,308 with a single bcgRNA.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Findingreasoning

    Model-based bcgRNA assignments (HTODemux or MultiSeqDemux) did not yield satisfying results supported by observed phenotypic changes, likely due to model limitations from high bcgRNA feature number.

    “Model-based bcgRNA assignments (HTODemux or MultiSeqDemux) did not yield satisfying results supported by the observed phenotypic changes, likely due to model limitations imposed by the high number of bcgRNA features.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Findingcitation

    Log2FC enrichments were generally more pronounced in samples using the enhanced direct repeat.

    “As we have previously described11, we noticed that log2FC enrichments were generally more pronounced in samples using the enhanced direct repeat.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [11] Wessels, H. H. et al. Massively parallel Cas13 screens reveal principles for guide RNA design. Nat. Biotechnol. 38, 722–727 (2020). source

      Not checked in the free run.

  • May combine two claims.

    Interpretationreasoning

    CaRPool-seq uses RNA-targeting Cas13d and permits efficient combinatorial perturbations together with multimodal single-cell profiling.

    “Here, we introduce Cas13 RNA Perturb-seq (CaRPool-seq), which leverages the RNA-targeting CRISPR–Cas13d system and enables efficient combinatorial perturbations alongside multimodal single-cell profiling.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Premisecitation

    The authors previously confirmed that RfxCas13d can cause strong target-RNA knockdown and derived optimal targeting rules from thousands of guide RNAs tiling multiple transcripts.

    “Recently, we confirmed that RfxCas13d can lead to striking target-RNA knockdown, and learned a set of optimal targeting rules from thousands of guide RNAs tiling multiple transcripts11.”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [11] Wessels, H. H. et al. Massively parallel Cas13 screens reveal principles for guide RNA design. Nat. Biotechnol. 38, 722–727 (2020). source

      Not checked in the free run.

  • May combine two claims.

    Premisecitation

    Cas9 perturbations are not uniformly efficient, and a considerable fraction of targeted cells may show no phenotypic effect.

    “Moreover, perturbations introduced by Cas9 are not uniformly efficient, and a considerable fraction of targeted cells may exhibit no phenotypic effects of perturbation8,9.”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [8] Michlits, G. et al. Multilayered VBC score predicts sgRNAs that efficiently generate loss-of-function alleles. Nat. Methods 17, 708–716 (2020). 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.

  • Classified as probably not a substantial claim.

    Premisecitation

    Combinatorial perturbations can reveal epistasis and other genetic interactions.

    “In particular, the ability to perform combinatorial perturbations represents an opportunity to decode complex regulatory networks, with pioneering work demonstrating the ability to identify epistasis and other genetic interactions5–7.”
    Quote in the paperSubstantial claim0 of 3 sources checked
    • cites [5] 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.

    • cites [6] Replogle, J. M. et al. Combinatorial single-cell CRISPR screens by direct guide RNA capture and targeted sequencing. Nat. Biotechnol. 38, 954–961 (2020). source

      Not checked in the free run.

    • cites [7] Norman, T. M. et al. Exploring genetic interaction manifolds constructed from rich single-cell phenotypes. Science 365, 786–793 (2019). source

      Not checked in the free run.

  • May combine two claims.

    Interpretationown data

    These findings support and provide a molecular explanation for observations that combination therapies of KDM1A antagonist and HDAC inhibitors show enhanced response.

    “These findings support and provide a molecular explanation for recent observations that combination therapies of KDM1A antagonist and HDAC inhibitors exhibit an enhanced response29.”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [29] Fiskus, W. et al. Highly effective combination of LSD1 (KDM1A) antagonist and pan-histone deacetylase inhibitor against human AML cells. Leukemia 28, 2155–2164 (2014). source

      Not checked in the free run.

    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
  • May combine two claims.

    Interpretationcitation

    Dual perturbation of KDM1A and GSE1 fell in cluster 1, consistent with a cooperative interaction via colocalization at repressed promoters to inhibit myeloid differentiation.

    “Dual perturbation of KDM1A and the transcriptional repressor GSE1 also fell in this cluster (Extended Data Fig. 8f), consistent with previous work that suggests a cooperative interaction via colocalization at repressed promoters to inhibit myeloid differentiation28.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [28] Nicosia, L. et al. Pharmacological inhibition of LSD1 triggers myeloid differentiation by targeting GSE1 oncogenic functions in AML. Oncogene https://doi.org/10.1038/s41388-021-02123-7 (2021). source

      Not checked in the free run.

    • Rests on an assumption: The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Study design and conductown data

    The single-gRNA pooled screen included four independent transductions, with CD14 sorts for all four and CD11b sorts for three.

    “In total, we prepared four independent transductions (two multiplicities of infection (MOI) and two alternative direct repeats), performed CD14 sorts for all four transduction replicates, and CD11b sorts for three transduction replicates collecting 1 × 106 to 1.5 × 106 cells per bin.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Study design and conductnone given

    During pooled screens, guide representation was maintained above 1,000×.

    “cells were split as needed always maintaining a guide representation of >1,000×.”
    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 conductreasoning

    The CaRPool-seq library contained 385 arrays with 186 single or double perturbations, each represented by two independent technical replicate gRNA combinations.

    “In total, we design 385 arrays with 186 single or double perturbations, each represented by two independent technical replicate gRNA combinations.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing). May combine two claims.

    Premisecitation

    Nelson et al. reported that structured RNA motifs on pegRNAs protect from exonuclease degradation, increase stability, and enhance efficiency.

    “While our paper was in review, Nelson et al.30 reported that the inclusion of structured RNA motifs on prime editing gRNAs (pegRNAs) led to protection from exonuclease degradation, increased stability, and enhanced efficiency.”
    Quote in the paperNot a substantial claim0 of 1 sources checked
    • cites [30] Nelson, J. W. et al. Engineered pegRNAs improve prime editing efficiency. Nat. Biotechnol. https://doi.org/10.1038/s41587-021-01039-7 (2021). source

      Not checked in the free run.

  • May combine two claims.

    Study design and conductown data

    The CaRPool-seq library achieved complete representation with minimal bias, with 90th/10th percentile crRNA read ratios of 2.6/4.8.

    “Complete library representation with minimal bias (90th percentile/10th percentile crRNA read ratio 2.6/4.8), and correct gene pair to bcgRNA linkage (>94%) was verified by Illumina sequencing (MiSeq).”
    Quote in the paperSubstantial claim
  • Classified as probably not a substantial claim.

    Premisecitation

    In prior Cas9 screens, 240 genes led to CD11b or CD14 upregulation.

    “We selected 240 genes that led to CD11b or CD14 upregulation in Cas9 screens23, in addition to 199 control genes in TLR4-signaling.”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [23] Wang, E. et al. Surface antigen-guided CRISPR screens identify regulators of myeloid leukemia differentiation. Cell Stem Cell 28, 718–731 (2021). source

      Not checked in the free run.

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

    Study design and conductown data

    The combinatorial perturbation pooled screen was performed in three transduction replicates at MOI 0.13-0.20.

    “For the combinatorial perturbation pooled screen, we performed three transduction replicates (MOI 0.13–0.20).”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Study design and conductown data

    The single-gRNA library achieved complete representation with minimal bias, with a 90th/10th percentile crRNA read ratio of 1.8 for both libraries.

    “Complete library representation with minimal bias (90th percentile/10th percentile crRNA read ratio of 1.8 for both libraries) was verified by Illumina sequencing (MiSeq).”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationown data

    The observed drop-off in combinatorial perturbation detection is consistent with the theoretical expectation for multiple independently detected gRNAs and highlights the challenge of profiling multiple perturbations with existing approaches.

    “This observed drop-off is fully consistent with the theoretical expectation of recovery for multiple independently detected gRNAs and highlights the challenge of efficiently profiling multiple perturbations with existing approaches.”
    Quote in the paperSubstantial claim
  • May combine two claims.

    Interpretationown data

    CaRPool-seq and Perturb-seq can both effectively introduce combinatorial perturbations into single cells, but CaRPool-seq has clear advantages in successfully identifying and detecting these perturbations.

    “We conclude that CaRPool-seq and Perturb-seq can both effectively introduce combinatorial perturbations into single cells. However, CaRPool-seq exhibits clear advantages in the ability to successfully identify and detect these perturbations and therefore represents an attractive approach for performing combinatorial single-cell CRISPR screens.”
    Quote in the paperSubstantial claim
  • Classified as probably not a substantial claim. May combine two claims.

    Premisecitation

    crRNAs are typically protected from degradation while complexed with Cas13, but the ends of longer bcgRNA molecules may still be accessible and susceptible to degradation.

    “While we expect that crRNA are typically protected from degradation while complexed with Cas13 (refs. 31,32), the ends of longer bcgRNA molecules may still be accessible and susceptible to degradation.”
    Quote in the paperSubstantial claim0 of 2 sources checked
    • cites [31] Méndez-Mancilla, A. et al. Chemically modified guide RNAs enhance CRISPR-Cas13 knockdown in human cells. Cell Chem. Biol. https://doi.org/10.1016/j.chembiol.2021.07.011 (2021). source

      Not checked in the free run.

    • cites [32] Zhang, C. et al. Structural basis for the RNA-guided ribonuclease activity of CRISPR-Cas13d. Cell 175, 212–223 (2018). source

      Not checked in the free run.

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

    Study design and conductnone given

    The bcgRNAs were designed as random 15mer sequences with hamming distance greater than four to one another.

    “For the bcgRNAs, we designed random 15mer sequences with hamming distance greater than four to one another.”
    Quote in the paperNot a substantial claim
  • May combine two claims.

    Study design and conductown data

    Cell survival after selection ranged between 1.7 and 5.5% at MOI <0.1, assuring high single integration probability.

    “Cell survival after selection ranged between 1.7 and 5.5% (MOI < 0.1) assuring a high single integration probability.”
    Quote in the paperSubstantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Study design and conductown data

    More than 98% of cells were BFP-positive before the 10X experiment.

    “We confirmed that >98% of cells were BFP-positive before the 10X experiment.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Study design and conductown data

    Cell viability was ≥96% before the direct capture Perturb-seq run.

    “Before the run, cell viability was determined (≥96%).”
    Quote in the paperNot a substantial claim
  • Classified as probably not a substantial claim.

    Premisecitation

    The linear regression model δab = c1δa + c2δb + ε was used to decompose double perturbation profiles, as previously introduced.

    “As previously introduced7, we model the average z-scale profiles using: δab=c1δa+c2δb+ϵ with δa and δb being the pseudobulk z-scaled profile for cells assigned to single perturbations a and b, repectively, while δab is the pseudobulk z-scaled profile for cells assigned to double perturbation ab.”
    Quote in the paperSubstantial claim0 of 1 sources checked
    • cites [7] Norman, T. M. et al. Exploring genetic interaction manifolds constructed from rich single-cell phenotypes. Science 365, 786–793 (2019). source

      Not checked in the free run.

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

    Study design and conductown data

    Cell viability was ≥95% for each CaRPool-seq experiment before the runs.

    “Before the runs, cell viability was determined ≥95% for each experiment.”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Interpretationnone given

    The Illumina 5′ PCR handle allows efficient sequencing of the bcgRNA amplicon with the first base of read2 being the first barcode base.

    “Moreover, usage of the Illumina 5′ PCR handle allows for efficient sequencing of the bcgRNA amplicon with the first base of read2 being the first barcode base.”
    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 single-gRNA pooled screen, the average mapping rate was 82.3% and the median gRNA count was 167.

    “string-match intersection with the reference to retain only perfectly matching alignments (average mapping rate 82.3%, median gRNA count 167).”
    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 combinatorial pooled screen, the average mapping rate was 97%, the median barcode read count was 635, and one barcode was not detected in input samples.

    “preprocessed reads were aligned to the barcode reference using bowtie43 (v.1.1.2) with parameters -v 1 -m 1–best –strata (average mapping rate 97%; median barcode read count 635; one barcode was not detected in input samples).”
    Quote in the paperNot a substantial claim
  • Probably not a substantial claim: may not pass S3 (load-bearing).

    Interpretationreasoning

    The bcgRNA design allowed specific amplification of a bcgRNA amplicon with a unique combination of forward and reverse primers.

    “This composition allowed the specific amplification of a bcgRNA amplicon with a unique combination of forward and reverse primers.”
    Quote in the paperNot a substantial claim
  • Classified as probably not a substantial claim. May combine two claims.

    Study design and conductcitation

    Raw counts were normalized using median of ratios as in DESeq2 and batch corrected using ComBat.

    “For each dataset, raw counts were normalized using a median of ratios method as in DESeq2 (ref. 44) and batch corrected using combat implemented in SVA (v.3.34.0)45.”
    Quote in the paperSubstantial claim
57 of 57
  • The paper relies on this without support.

    Cas13 crRNAs lack a common priming site for reverse transcription.

    “Type VI-A, C and D Cas13 crRNAs consist of a short 5′ direct repeat and a variable spacer (also called gRNA) at the 3′ end, and therefore lack a common priming site for reverse transcription (RT).”
  • The paper relies on this without support.

    Target RNAs are continuously produced and degraded before translation.

    “Target RNAs are continuously being produced and degraded before the target can be translated into protein.”
  • The paper relies on this without support.

    The THP1 MLL-AF9 NRASG12D cell line is a valid model system for AML differentiation.

    “We generated a human MLL-AF9 NRASG12D AML cell line (THP1 cells), with a stably integrated doxycycline-inducible Cas13d cassette, as a model system.”
  • The paper relies on this without support.

    CD14 and CD11b are immunophenotypic markers of monocyte differentiation.

    “cells were sorted into bins based on their surface expression of CD14 and CD11b, immunophenotypic markers of monocyte differentiation.”
  • The paper relies on this without support.

    gRNA matches shorter than 17 bases do not cause target knockdown.

    “Second, candidates were further filtered to match with at least 17 bases, as shorter matches do not lead to target knockdown and show a blastn e.value of <100.”
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?

  • High multiplicity-of-infection CRISPR screens with multiple targets per cell can complicate the calculation of distances.

    “For example, high multiplicity-of-infection CRISPR screens with multiple targets per cell can complicate the calculation of distances5,53,54.”

    From scPerturb: harmonized single-cell perturbation data (2024) · cites this paper as [5]

Cites

Flagged first, then by citations

Cited by

Sorted by citations

Authors
Wessels, H.-H., Méndez-Mancilla, A., Hao, Y., Papalexi, E., Mauck, W. M., Lu, L., Morris, J. A., Mimitou, E. P., Smibert, P., Sanjana, N. E., & Satija, R.
Affiliations
New York Genome Center (US)
Corresponding author
  • Sanjana, Neville E.
  • Satija, Rahul
Venue
Nature Methods
Published
2022-12-22
Volume
20
Issue
1
Pages
86-94
ISSN
  • 1548-7091 (ISSN-L)
  • 1548-7105
OpenAlex
W4312114042
Cited by
79
References
52
Citations by year
2026: 14 · 2025: 30 · 2024: 22 · 2023: 9 · 2022: 4
FWCI
5.31
Type
journal article
Language
en
Open access
Locations
  • Nature Methods (publishedVersion)
  • bioRxiv (Cold Spring Harbor Laboratory) (acceptedVersion)
  • PubMed (publishedVersion)
  • PubMed Central (acceptedVersion)
Funders
  • Cancer Research Institute
  • Chan Zuckerberg Initiative (HCA-A1704-01895, EOSS-0000000082)
  • York University
  • National Institutes of Health (1dp2hg009623-01)
  • Defense Advanced Research Projects Agency (D18AP00053)
  • National Human Genome Research Institute (DP2-HG-009623, DP2HG010099, RM1HG011014)
  • National Cancer Institute (DP2HG009623-01, R01CA218668)
Topics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › CRISPR and Genetic Engineering
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Single-cell and spatial transcriptomics
  • Life Sciences › Biochemistry, Genetics and Molecular Biology › Molecular Biology › Chromatin Remodeling and Cancer
Keywords
CRISPR-Cas13d, single-cell RNA sequencing, RNA-targeting CRISPR, barcode encoding, acute myeloid leukemia, myeloid differentiation, chromatin regulators, gene regulatory networks
MeSH terms
Chromatin, RNA, CRISPR-Cas Systems

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