Abstract
Multimodal single-cell assays provide high-resolution snapshots of complex cell populations, but are mostly limited to transcriptome plus an additional modality. Here, we describe expanded CRISPR-compatible cellular indexing of transcriptomes and epitopes by sequencing (ECCITE-seq) for the high-throughput characterization of at least five modalities of information from each single cell. We demonstrate application of ECCITE-seq to multimodal CRISPR screens with robust direct single-guide RNA capture and to clonotype-aware multimodal phenotyping of cancer samples.
This is a preview of subscription content, access via your institution
Access options
Access Nature and 54 other Nature Portfolio journals
Get Nature+, our best-value online-access subscription
$29.99 / 30 days
cancel any time
Subscribe to this journal
Receive 12 print issues and online access
$259.00 per year
only $21.58 per issue
Buy this article
- Purchase on Springer Link
- Instant access to full article PDF
Prices may be subject to local taxes which are calculated during checkout
Similar content being viewed by others
Data availability
Data generated in this project have been deposited to the Gene Expression Omnibus with the accession code GSE126310.
References
Hu, Y. et al. Front. Cell. Dev. Biol. 6, 28 (2018).
Kester, L. & van Oudenaarden, A. Cell Stem Cell 23, 166–179 (2018).
Jaitin, D. A. et al. Cell 167, 1883–1896.e15 (2016).
Adamson, B. et al. Cell 167, 1867–1882.e21 (2016).
Dixit, A. et al. Cell 167, 1853–1866.e17 (2016).
Datlinger, P. et al. Nat. Methods 14, 297–301 (2017).
Hill, A. J. et al. Nat. Methods 15, 271–274 (2018).
Stoeckius, M. et al. Nat. Methods 14, 865–868 (2017).
Peterson, V. M. et al. Nat. Biotechnol. 35, 936–939 (2017).
Stoeckius, M. et al. Genome Biol. 19, 224 (2018).
Jinek, M. et al. Science 337, 816–821 (2012).
Cong, L. et al. Science 339, 819–823 (2013).
Butler, A., Hoffman, P., Smibert, P., Papalexi, E. & Satija, R. Nat. Biotechnol. 36, 411–420 (2018).
Fanok, M. H. et al. J. Investig. Dermatol. 138, 1116–1125 (2018).
Nam, A. S. et al. High throughput droplet single-cell Genotyping of Transcriptomes (GoT) reveals the cell identity dependency of the impact of somatic mutations. Preprint at https://doi.org/10.1101/444687 (2018).
Wroblewska, A. et al. Cell 175, 1141–1155.e16 (2018).
Replogle, J. M. et al. Direct capture of CRISPR guides enables scalable, multiplexed, and multi-omic Perturb-seq. Preprint at https://doi.org/10.1101/503367 (2018).
Shen, J. P. et al. Nat. Methods 14, 573–576 (2017).
Najm, F. J. et al. Nat. Biotechnol. 36, 179–189 (2018).
Kalhor, R. et al. Science 361, 893 (2018).
Van Buggenum, J. A. G. L. et al. Sci. Rep. 6, 22675 (2016).
Abrams, J. T. et al. J. Investig. Dermatol. 96, 31–37 (1991).
Kaltoft, K. et al. In Vitro Cell. Dev. Biol. 28A, 161–167 (1992).
Sanjana, N. E., Shalem, O. & Zhang, F. Nat. Methods 11, 783–784 (2014).
Patel, S. J. et al. Nature 548, 537–542 (2017).
Acknowledgements
We thank N. Ødum (University of Copenhagen) for the kind gift of cell lines. Samples from patients were obtained with the help of M.B. Natan Zommer and J.-A. Latkowski. Work in S.B.K.’s laboratory was supported by the NIH R01 grant no. HL-125816, the Colton Center for Autoimmunity, funding from the Hirschl/Weill-Coulier Trust and a grant from the Spatz Foundation. Work in the NYGC Technology Innovation lab was supported by the NIH R21 grant no. HG-009748 to P.S. and the Chan Zuckerberg Initiative grant no. HCA-A-1704-01895 to P.S. and R.S. N.E.S. is supported by NYU and NYGC startup funds, NIH/NHGRI (R00HG008171 and DP2HG010099), NIH/NCI (R01CA218668), DARPA (D18AP00053), the Sidney Kimmel Foundation, the Melanoma Research Alliance, and the Brain and Behavior Foundation. M.L. is supported by a Hope Funds for Cancer Research postdoctoral fellowship. Work in Z.O.’s laboratory was supported by NIH R35 grant GM124998. We thank L. Yang, W. Stephenson, S. Jaini and K. Pandit for helpful discussions. We thank B. Fritz from 10x Genomics for providing kits for development of 5P compatible CITE-seq reagents and B. Yeung and K. Nazor from BioLegend for providing some of the unconjugated antibodies used in this study.
Author information
Authors and Affiliations
Contributions
E.P.M. and P.S. conceived and designed the study with input from A.M., S.H., M.S., E.P., R.S., N.E.S. and S.B.K. E.P.M. performed all experiments, aided by S.H. M.S. A.C., A.H., Z.O. and S.B.K. provided samples and performed analysis on the CTCL study. A.M., M.L., T.R. and N.E.S. worked with E.P.M. to design and generate CRISPR libraries and provide experimental and analytical guidance. E.P.M. and P.S. wrote the paper.
Corresponding author
Ethics declarations
Competing interests
M.S. and P.S. are listed as co-inventors on a patent application related to this work (US provisional patent application 62/515–180).
Additional information
Publisher’s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Integrated supplementary information
Supplementary Fig. 1 ECCITE-seq enables CRISPR screens with single-cell multimodal readout.
a, Species-mixing proof-of-principle experiment: protein tag reads associated with each cell barcode. Points are colored based on species classification using transcripts as shown in Fig. 1b. About 1.2% of either human or mouse cells show cross-reactivity with mouse or human antibodies respectively. b, Mixed-species ECCITE-seq experiment demonstrating measurement of six cellular modalities. Cells of different origin were stained with ECCITE-seq and hashing antibodies, washed and combined to a single 10x run. All cells were stained with a mix of anti-human CD29 and anti-mouse CD29 antibodies. NIH-3T3 cells were split into 7 tubes and stained with 7 barcoded hashing antibodies (Hashtag-A to Hashtag-G), followed by washing and pooling. MyLa, Sez4 and PBMCs were stained with Hashtag_1, Hashtag_2 and Hashtag_3 respectively. Eventually, a cell mix comprising 65% NIH-3T3 cells, 25% PBMCs, and 5% each MyLa and Sez4 cells was loaded into the 10x Chromium for droplet formation and reverse transcription. After emulsion breakage and cDNA amplification, the distinctly sized products were separated with size selection capture beads and amplified separately. A pool consisting 88% cDNA, 7% guide-tag, 3% hashtag and 2% protein-tag library was submitted for NGS sequencing. c,d, sgRNA representation as measured by direct sgRNA capture or genomic DNA amplification of the guide variable region in two different cell lines. To assess the direct capture, we assigned a guide to each cell with a single sgRNA and quantified their proportion in the total cell population (orange). This ratio is compared to the ratio of genomic reads matching each sgRNA to the total genomic sgRNA reads. e, scRNA-seq saturation curves for the ECCITE-seq versus standard 10x V(D)J run using the K562 cells expressing the targeting sgRNA library. The same cell suspension was used into two parallel experiments; one underwent a standard 10x V(D)J run while the other was preceded by staining with CITE-seq and hashing antibodies, as well as spike-in of guide RT primer in the reverse transcription reaction. For this analysis, data were downsampled to the same total reads per cell after correcting for the different fraction of reads in cell. f. Levels of target mRNA or protein (moving median smoothing window = 151) in 1,000 cells ranked on decreasing sgRNA counts for CD46 sgRNA 2 or non-targeting sgRNA 4. g, Adjusted p-values of detecting the indicated gene expression or protein abundance change as a function of cell number. Different-size, randomized cell samples assigned to unique targeting sgRNAs were tested against equal size, randomized cell samples from the non-targeting sgRNAs groups, using the Wilcoxon rank sum test.
Supplementary Fig. 2 Surface protein and clonotype detection on PBMCs from a healthy donor and a CTCL patient.
a, Transcriptome-based clustering of PBMCs from healthy donor after removing cell doublets. Projected is the protein signal for 36 out of 49 antibodies used to stain the cells, as well as the 2nd and 3rd most abundant CD4+ or CD8+ TCR α/β or the 3rd, 4th, 5th and 6th most abundant TCR γ/δ clonotype (red). In total, we detected 1,606 TCR α/β clonotypes from 2,796 barcodes and 183 TCR γ/δ clonotypes from 269 barcodes. The top CD4+ TCR α/β clonotype (defined by TRB CDR3 sequence: CASSTLQGKETQYF, shown in Fig. 2a) accounts for ~1% of recovered clonotype-associated barcodes. b, Same analysis as panel a, using data from CTCL patient PBMCs. In total, we detected 1,738 TCR α/β clonotypes from 3,857 barcodes, and 87 TCR γ/δ clonotypes from 243 barcodes. Clonal expansion is readily apparent by the top CD4+ TCR α/β clonotype (defined by TRB CDR3 sequence: CSARFLRGGYNEQFF, shown in Fig. 2a) present in 36% of cells for which we recovered clonotype information.
Supplementary Fig. 3 Gene expression and surface marker-based clustering of the combined dataset.
a, Unsupervised clustering of PBMCs from both healthy and CTCL donors (n = 9,816) after removing cell duplicates, merging, depth normalization and cell alignment13. Each cell was colored and labelled based on unsupervised clustering information on gene expression (left) or surface marker (right). b, Heatmap of genes differentially expressed across gene expression-based (left) or surface marker-based (right) cluster assignments.
Supplementary information
Supplementary Text and Figures
Supplementary Figures 1–3 and Supplementary Tables 1–5
Supplementary Protocol
ECCITE-seq protocol
Rights and permissions
About this article
Cite this article
Mimitou, E.P., Cheng, A., Montalbano, A. et al. Multiplexed detection of proteins, transcriptomes, clonotypes and CRISPR perturbations in single cells. Nat Methods 16, 409–412 (2019). https://doi.org/10.1038/s41592-019-0392-0
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1038/s41592-019-0392-0
This article is cited by
-
deMULTIplex2: robust sample demultiplexing for scRNA-seq
Genome Biology (2024)
-
Integration of multi-modal single-cell data
Nature Biotechnology (2024)
-
Tumor reactive γδ T cells contribute to a complete response to PD-1 blockade in a Merkel cell carcinoma patient
Nature Communications (2024)
-
scPerturb: harmonized single-cell perturbation data
Nature Methods (2024)
-
scPROTEIN: a versatile deep graph contrastive learning framework for single-cell proteomics embedding
Nature Methods (2024)