Dec 07, 2023
Perturb-seq characterizing regulators of T cell function
- Sean R. McCutcheon1,
- Adam M. Swartz1,
- Michael C. Brown1,
- Alejandro Barrera1,
- Christian McRoberts Amador1,
- Keith Siklenka1,
- Lucas Humayun1,
- Maria A. ter Weele1,
- James M. Isaacs1,
- Andrea R Daniel1,
- Timothy E. Reddy1,
- Andrew S. Allen1,
- Smita K. Nair1,
- Scott J. Antonia1,
- Charles A. Gersbach1
- 1Duke University
- Andrea R Daniel: This protocol was adapted from the work of Sean McCutcheon and colleagues in the Gersbach lab at Duke University.
- Gersbach Lab
Protocol Citation: Sean R. McCutcheon, Adam M. Swartz, Michael C. Brown, Alejandro Barrera, Christian McRoberts Amador, Keith Siklenka, Lucas Humayun, Maria A. ter Weele, James M. Isaacs, Andrea R Daniel, Timothy E. Reddy, Andrew S. Allen, Smita K. Nair, Scott J. Antonia, Charles A. Gersbach 2023. Perturb-seq characterizing regulators of T cell function. protocols.io https://dx.doi.org/10.17504/protocols.io.dm6gp3qjdvzp/v1
Manuscript citation:
McCutcheon, S.R., Swartz, A.M., Brown, M.C. et al. Transcriptional and epigenetic regulators of human CD8+ T cell function identified through orthogonal CRISPR screens.Nat Genet (2023). https://doi.org/10.1038/s41588-023-01554-0
License: This is an open access protocol distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited
Protocol status: Working
We use this protocol and it's working
Created: December 01, 2023
Last Modified: December 07, 2023
Protocol Integer ID: 91713
Keywords: Perturb-seq, CRISPRi, CRISPRa, T cells, scRNA-seq
Funders Acknowledgements:
NIH
Grant ID: HG012053
Abstract
This protocol describes methods for a Perturb-seq assay characterizing transcriptional regulators of T cell function.
Materials
Addgene: pLV hU6-gRNA hUbC-dSaCas9-KRAB-T2A-Thy1.1 (Addgene 194278) and pLV hU6-gRNA hUbC-VP64-dSaCas9-VP64-T2A-Thy1.1 (Addgene 194279).
gRNA library cloning
gRNA library cloning
Oligonucleotide pools containing 40-gRNA sequences and constant regions for polymerase chain reaction (PCR) amplification were synthesized by Twist Bioscience.
gRNA amplicons were gel extracted, PCR purified and input into 20 μl Gibson reactions (5:1 molar ratio of insert to backbone) with 200 ng of Esp3I digested and 1 × solid-phase reversible immobilization (SPRI)-selected (Beckman Coulter) plasmid backbone.
Addgene: pLV hU6-gRNA hUbC-dSaCas9-KRAB-T2A-Thy1.1 (Addgene 194278) and pLV hU6-gRNA hUbC-VP64-dSaCas9-VP64-T2A-Thy1.1 (Addgene 194279).
Gibson reactions were purified using ethanol precipitation and transformed into Lucigen’s Endura ElectroCompetent Cells.
Transformed cells were cultured overnight and plasmids were isolated using Qiagen Midi Kits.
Transfections for high-titer lentiviral production
Transfections for high-titer lentiviral production
Plate 1.2 x 106 or 7 x 106 HEK293T cells in a 6 well plate or 10 cm dish in the afternoon with 2 mL or 12 mL of complete opti-MEM (Opti-MEM‱ I Reduced Serum Medium supplemented with 1x Glutamax, 5% FBS, 1 mM Sodium Pyruvate, and 1x MEM Non-Essential Amino Acids).
The next morning, transfect HEK293T cells with 0.5 μg pMD2.G, 1.5 μg psPAX2, and 0.5 μg transgene for 6 well plates or 3.25 μg pMD2.G, 9.75 μg psPAX2, and 4.3 μg transgene for 10 cm dishes using Lipofectamine 3000.
Exchanged media 6 hours after transfection and collect and pool lentiviral supernatant at 24 hours and 48 hours after transfection.
Transduction of primary human T cells
Transduction of primary human T cells
Centrifuged lentiviral supernatant at 600g for 10 min to remove cellular debris.
Concentrate lentivirus to 50–100× the initial concentration using Lenti-X Concentrator (Takara Bio).
Transduce T cells at 5–10% v/v of concentrated lentivirus at 24 h post-activation. For dual transduction experiments, T cells were serially transduced at 24 h and 48 h.
scRNAseq
scRNAseq
CD8+CCR7+ T cells from three donors were transduced with CRISPRi and CRISPRa mini-TF gRNA libraries.
T cells were expanded for 10 days and then stained and sorted for Thy1.1+ cells (a marker to identify transduced cells).
Antibody: Thy1.1-PE, clone OX-7, Mouse IgG1, κ, 1:300, Flow cytometry, StemCell Technologies, 60024PE
An SH800 FACS Cell Sorter (Sony Biotechnology) was used for cell sorting and analysis.
For antibody staining of Thy1.1 cells were collected, spun down at 300g for 5 min, resuspended in flow buffer (1× phosphate-buffered saline (PBS), 2 mM ethylenediaminetetraacetic acid and 0.5% bovine serum albumin) with the appropriate antibody dilutions and incubated for 30 min at 4 °C on a rocker.
Cells were then washed with flow buffer, spun down at 300g for 5 min and resuspended in flow buffer for cell sorting.
Sorted cells were loaded into the Chromium X for a targeted recovery of 2 × 104 cells per donor and treatment according to the Single Cell 5′-High-Throughput Reagent Kit v2 protocol (10x Genomics).
SaCas9 gRNA sequences were captured by spiking in 2 μM of a custom primer into the reverse transcription master mix, as previously done for SpCas9 gRNA capture36. The custom primer was designed to bind to the constant region of SaCas9’s gRNA scaffold.
5′-Gene Expression (GEX) and gRNA libraries were separated using double-sided SPRI bead selection in the initial cDNA clean-up step. 5′-GEX libraries were constructed according to manufacturer’s protocol.
gRNA libraries were constructed using two sequential PCRs (PCR 1: 10 cycles, PCR 2: 25 cycles). The PCR 1 product was purified using double-sided SPRI bead selection at 0.6 × and 2 ×. Twenty percent of the purified PCR 1 product was input into PCR 2. The PCR2 product was purified using double-sided SPRI bead selection at 0.6 × and 1 ×.
All libraries were run on a High Sensitivity D1000 tape to measure the average amplicon size and quantified using Qubit’s dsDNA High Sensitivity assay.
Libraries were individually diluted to 20 nM, pooled together at desired ratios and sequenced on an Illumina NovaSeq S4 Full Flow Cell (200 cycles) with the following read allocation: Read 1, 26; i7 index, 10; Read 2, 90.
All oligos are available in Supplementary Table 5 of McCutcheon et al. Nature Genetics, 2023. https://doi.org/10.1038/s41588-023-01554-0
Processing and analyzing scRNA-seq
Processing and analyzing scRNA-seq
CellRanger v6.0.1 was used to process, demultiplex and generate UMI counts for each transcript and gRNA per cell barcode.
UMI counts tables were extracted and used for subsequent analyses in R using the Seurat69 v4.1.0 package.
Low-quality cells with <200 detected genes, >20% mitochondrial reads or <5% ribosomal reads were discarded.
DoubletFinder70 was used to identify and remove predicted doublets. All remaining high-quality cells across donors for each treatment (CRISPRi or CRISPRa) were aggregated for further analyses.
gRNAs were assigned to cells if they met the threshold (gRNA UMI >4). Cells were then grouped on the basis of gRNA identity.
For differential gene expression analysis, we compared the transcriptomic profiles of cells sharing a gRNA to cells with only NT gRNAs using Seurat’s FindMarkers function to test for DEGs with the hurdle model implemented in model-based analysis of single-cell transcriptomics (MAST).
All significant gRNA–gene links are available in Supplementary Table 3 of McCutcheon et al. Nature Genetics, 2023. https://doi.org/10.1038/s41588-023-01554-0
Upregulated DEGs were input into EnrichR’s GO Biological Process 2021 database71 for functional annotation.
Protocol references
36. Mimitou, E. P. et al. Multiplexed detection of proteins, transcriptomes, clonotypes and CRISPR perturbations in single cells. Nat. Methods 16, 409–412 (2019).
69. Butler, A., Hoffman, P., Smibert, P., Papalexi, E. & Satija, R. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Nat. Biotechnol. 36, 411–420 (2018).
70. McGinnis, C. S., Murrow, L. M. & Gartner, Z. J. DoubletFinder-doublet detection in single-cell RNA sequencing data using artificial nearest neighbors. Cell Syst. 8, 329–337 (2019).
71. Kuleshov, M. V. et al. Enrichr—a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 44, 90–97 (2016).