| Literature DB >> 33899013 |
Celia Alda-Catalinas1, Melanie A Eckersley-Maslin1, Wolf Reik1,2,3.
Abstract
CRISPR/Cas9 screens are a powerful approach to identify key regulators of biological processes. By combining pooled CRISPR/Cas9 screening with single-cell RNA-sequencing readout, individual perturbations can be assessed in parallel both comprehensively and at scale. Importantly, this allows gene function and regulation to be interrogated at a cellular level in an unbiased manner. Here, we present a protocol to perform pooled CRISPR-activation screens in mouse embryonic stem cells using 10× Genomics scRNA-seq as a readout. For complete information on the generation and use of this protocol, please refer to Alda-Catalinas et al. (2020).Entities:
Keywords: CRISPR; Cell culture; Gene Expression; High Throughput Screening; RNAseq; Sequencing; Single Cell; Stem Cells
Mesh:
Year: 2021 PMID: 33899013 PMCID: PMC8055612 DOI: 10.1016/j.xpro.2021.100426
Source DB: PubMed Journal: STAR Protoc ISSN: 2666-1667
Figure 1Screening for genomic integrations of dCas9-VP64_Blast and MS2-p65-HSF1 by PCR and agarose gel electrophoresis (step 12d)
Agarose gel electrophoresis (1% agarose) showing the products resulting from PCR on genomic DNA from selected clones using the primers and conditions described in step 12, designed to amplify a region within dCas9-VP64_Blast integration (top panel) and a region within MS2-p65-HSF1_Hygro integration (bottom panel). The ladder used is HyperLadder IV (Biolane, BIO-33029); the numbers reflect different clones; the PCR positive controls (+) were plasmid DNA of dCas9-VP64_Blast (top panel) or plasmid DNA of MS2-p65-HSF1_Hygro (bottom panel); the PCR negative controls (-) were genomic DNA from the parental line (non-transduced) E14 mESCs; non-clonal (NC) shows the PCR product from the respective polyclonal line. The expected band size for a clone containing dCas9-VP64_Blast integrations is 275 bp and for MS2-p65-HSF1_Hygro is 672 bp. Asterisks above numbers indicate that those clones show clear bands of the right size to be selected for both dCas9-VP64_Blast and MS2-p65-HSF1_Hygro integrations. See also Table S1 for primer sequences.
Figure 2dCas9-VP64 and MS2-p65-HSF1 expression of two SAM mESC clones compared to the parental line E14, analyzed by q-RT-PCR (step 13e)
Analysis of dCas9-VP64 (green) and MS2-p65-HSF1 (pink) relative mRNA levels by q-RT-PCR in two different SAM mESC clones against its parental (non-transduced) cell line E14. Each clone was analyzed in three experimental replicates and the q-RT-PCR performed in two technical replicates. Data are shown as relative expression to the parental line E14 mESCs, after normalization of each transgene’s CT values to average CT values of the house-keeping genes Gapdh and CycloB1. Each dot represents the value for each experimental replicate, and bars represent mean plus and minus standard deviation. The dotted horizontal line at y=1000 indicates the minimum relative expression of both dCas9-VP64 and MS2-p65-HSF1 that a clone should have to be selected for further characterization and/or use in a screen. In this example, clone 22 was selected over clone 7. See also Table S1 for primer sequences.
Figure 3Transcriptome of two SAM mESC clones compared to the parental line E14 (step 14)
Scatterplots showing normalized gene expression in log2 reads per million (RPM) of two different SAM mESC clones against its parental (non-transduced) cell line E14 (top panels), as well as scatterplot between two replicates of non-transduced E14 (bottom panel), analyzed by RNA-sequencing, highlighting dCas9-VP64 and MS2-p65-HSF1 transcripts in black and differentially expressed endogenous genes (DEGs) in blue. The transcriptome analysis was done with two replicates of each SAM mESC clone and three replicates of E14. DEGs were determined using DESeq2 (FDR<5%), EdgeR (FDR<5%) and intensity difference filter (FDR<5%), with the high-confidence DEGs defined as the intersection between the three statistical tests and a log2 fold change >2. In this example, clone 22 was selected over clone 7 due to having less DEGs and, by q-RT-PCR, showing higher relative expression levels of dCas9-VP64 and MS2-p65-HSF1 (see also Figure 2).
Figure 4Flow cytometry gating strategy (steps 25, 27, 31 g, and 32 b)
Flow cytometry scatterplots and histograms showing a representative gating strategy for steps 25, 27, 31 g, and 32 b. After gating on single cells, DAPI- cells are selected (viable cells) and, within those, mCherry fluorescence is analyzed. mCherry+ gate should be set based on a control of non-transduced cells.
Figure 5mCherry histograms for lentiviral custom sgRNA library titration, analyzed by flow cytometry (step 25)
Flow cytometry histograms showing expression of mCherry in SAM mESCs transduced with different volumes of a custom sgRNA library cloned into CROP-sgRNA-MS2 backbone, analyzed three days after transduction and prior to puromycin selection; percentages of cells transduced (expressing mCherry) are shown.
Figure 6mCherry histograms for lentiviral custom sgRNA library titration and puromycin selection, analyzed by flow cytometry (step 27)
Flow cytometry histograms showing expression of mCherry in SAM mESCs transduced with different volumes of a custom sgRNA lentiviral library cloned into CROP-sgRNA-MS2 backbone, and selected for 4 days with 1 μL/mL puromycin; percentages of cells transduced and selected (expressing mCherry) are shown.
Figure 8Bioanalyzer traces of final sgRNA amplicon libraries (step 49l)
Bioanalyzer electropherogram showing fluorescence units (FU, y axis) vs fragment size in base pairs (bp, x axis) of a representative sgRNA amplicon library of SAM mESCs transduced with a custom CROP-sgRNA-MS2 lentiviral library. Numbers above electropherogram peaks indicate the size in bp of detected fragments, including a lower molecular weight marker of 35 bp (green) and a higher molecular weight marker of 10,380 bp (purple). sgRNA amplicon libraries derived from cells transduced with CROP-sgRNA-MS2 constructs generate amplicons of approximately 520 bp.
| REAGENT or RESOURCE | SOURCE | IDENTIFIER |
|---|---|---|
| Murine LIF | Wellcome – MRC Cambridge Stem Cell Institute | |
| TransIT transfection reagent | Mirus Bio | Cat# MIR2700 |
| Polybrene | Millipore | Cat# TR-1003-G |
| Lenti-X Concentrator | Takara | Cat# 631231 |
| DNaseI | Thermo Fisher Scientific | Cat# EN0521 |
| Blasticidin | InvivoGen | Cat# ant-bl-1 |
| Hygromycin | InvivoGen | ant-hg-1 |
| Puromycin | InvivoGen | ant-pr-1 |
| 2× KAPA HiFi Master Mix | Kapa Biosystems | Cat# KR0389 |
| Brilliant III SYBR master mix | Agilent Technologies | Cat# 600882 |
| AllPrep DNA/RNA Mini Kit | Qiagen | Cat# 80204 |
| RevertAid First-Strand cDNA synthesis kit | Thermo Fisher Scientific | Cat# K1622 |
| Chromium Next GEM Single Cell 3ʹ GEM, Library & Gel Bead Kit v3.1 | 10× Genomics | Cat# PN-1000121 or Cat#PN-1000128 for 16 or 4 reactions, respectively |
| Chromium Next GEM Chip G Single Cell Kit | 10× Genomics | Cat# PN-1000120 |
| Chromium Next GEM Chip G Single Cell Kit | 10× Genomics | Cat# PN-1000127 |
| Single Index Kit T Set A | 10× Genomics | Cat# PN-1000213 |
| AMPure XP beads | Beckman Coulter | Cat# A63881 |
| SPRIselect Reagent Kit | Beckman Coulter | Cat# B23318 |
| Nextera XT Index Kit | Illumina | Cat# FC-131-2001, FC-131-2002, FC-131-2003, FC-131-2004 or FC-131-1001 |
| Agilent Bioanalyzer High Sensitivity chip | Agilent | Cat# 5067-4626 |
| Raw and quantified sequencing data: bulk RNA-seq data of E14 and SAM mouse ESCs | GEO: | |
| Raw and quantified sequencing data: 10× Genomics CRISPRa screen dataset | GEO: | |
| Mouse: E14 embryonic stem cells | RRID: CVCL_C320 | |
| Human: HEK293T | ATCC | ATCC CRL-3216 |
| Genomic PCR primers for dCas9-VP64 and MS2-p65-HSF1 | ||
| q-RT-PCR primers | This paper | |
| Primers for sgRNA amplicon libraries | ||
| Primers for PCR amplification of a pooled sgRNA library to be cloned in CROP-sgRNA-MS2 lentiviral backbone | This paper | |
| pMD2.G | Didier Trono | Addgene plasmid #12259 |
| psPAX2 | Didier Trono | Addgene plasmid #12260 |
| Lenti dCas9-VP64_Blast | Addgene plasmid #61425 | |
| Lenti MS2-p65-HSF1_Hygro | Addgene plasmid #61426 | |
| Lenti CROP-sgRNA-MS2 | Addgene plasmid #153457 | |
| CellRanger v2.1 | ( | Github: |
| Scripts to analyse CROP-seq screens | Github: | |
| Scripts for quality control and analysis of CRISPRa CROP-seq | Github: | |
| 0.45 μm PES filter | Sartorius | Cat# 16533 |
| PCR for dCas9-VP64_Blast: | ||
|---|---|---|
| Reagent | Final concentration / amount | Amount |
| 10× DreamTaq Buffer | 1× | 2.5 μl |
| dNTP Mix (10 mM each) | 0.2 mM | 0.5 μl |
| dCas9-VP64_F primer (10 μM) | 0.4 μM | 1 μl |
| dCas9-VP64_R primer (10 μM) | 0.4 μM | 1 μl |
| Template DNA (genomic DNA from clones, or genomic DNA from non-transduced mESCs, or lenti dCas9-VP64_Blast plasmid DNA) | 0.2–1 μg for genomic DNA samples; | Variable |
| DreamTaq DNA Polymerase | 0.625 U | 0.125 μl |
| Nuclease-free water | n/a | Up to 25 μl |
| n/a | ||
sequence of forward primer (Table S1): 5′-CATCGAGCAGATCAGCGAGT-3′
sequence of reverse primer (Table S1): 5′-CGATCCGTGTCTCGTACAGG-3′
| PCR for MS2-p65-HSF1_Hygro | ||
|---|---|---|
| Reagent | Final concentration / amount | Amount |
| 10× DreamTaq Buffer | 1× | 2.5 μl |
| dNTP Mix (10 mM each) | 0.2 mM | 0.5 μl |
| MS2-p65-HSF1_F primer (10 μM) | 0.4 μM | 1 μl |
| MS2-p65-HSF1_R primer (10 μM) | 0.4 μM | 1 μl |
| Template DNA (genomic DNA from clones, or genomic DNA from non-transduced mESCs, or lenti MS2-p65-HSF1_Hygro plasmid DNA) | 0.2–1 μg for genomic DNA samples; | Variable |
| Nuclease-free water | n/a | Up to 25 μl |
| n/a | ||
sequence of forward primer (Table S1): 5′-AAGCCTGAACTCACCGCTAC-3′
sequence of reverse primer (Table S1): 5′-TTGGGAATCCCCGAACATGG-3′
| PCR cycling conditions | |||
|---|---|---|---|
| Steps | Temperature | Time | Cycles |
| Initial Denaturation | 95°C | 1 min | 1 |
| Denaturation | 95°C | 30 s | 35 |
| Annealing | 58°C | 1 min | |
| Extension | 72°C | 1 min | |
| Final Extension | 72°C | 5 min | 1 |
| Hold | 4°C | Forever | |
| Target gene | sgRNA protospacer | Mean log2 fold change expression to non-targeting controls in 10× Genomics scRNA-seq screen data described in |
|---|---|---|
| GAACTGGCTGGGATTGCGCA | 1.83 | |
| ATTTGGTAAGAGATGGCTTT | 1.24 | |
| TCCAGCCGCTGCCGAGTCCC | 0.94 | |
| GGCAGTCGGGTTCCCTCAGG | 0.76 | |
| GTCAGAGACCAGAGCAGTGG | 0.64 | |
| TAAGGAATGGCCAGCTGGCT | 0.63 | |
| GGAGGTGAAGAGGACGTTCT | 0.60 | |
| GCGTAGGTAGCTGTCCTGGT | 0.50 | |
| ACACAGGAGGACCCTCCCTC | 0.50 | |
| TTCGGAGCGAGCTCCAGCCC | 0.46 | |
| Non-targeting control | GCTTTCACGGAGGTTCGACG | N.A |
| Non-targeting control | ATGTTGCAGTTCGGCTCGAT | N.A |
| sgRNA amplicon PCR1 | ||
|---|---|---|
| Reagent | Final concentration / amount | Amount |
| Amp_sgRNA_1F primer (10 μM) | 1.5 μM | 1.5 μl |
| Amp_sgRNA_1R primer (10 μM) | 1.5 μM | 1.5 μl |
| 2× KAPA HiFi Master Mix (Kapa Biosystems, KR0389) | 1× | 5 μl |
| Full-length cDNA from 10× Genomics protocol step 2.3m | 10 ng | Variable |
| Nuclease-free water | n/a | Up to 10 μl |
| n/a | ||
sequence of forward primer (Table S1): TTTCCCATGATTCCTTCATATTTGC (primes an outer part of the U6 promoter).
sequence of reverse primer (Table S1): CTACACGACGCTCTTCCGATCT (partial TruSeq Read 1 primer).
| PCR1 Cycling Conditions | |||
|---|---|---|---|
| Steps | Temperature | Time | Cycles |
| Initial Denaturation/ Enzyme activation | 95°C | 5 min | 1 |
| Denaturation | 95°C | 30 s | Variable (monitor amplification in the screen of the q-RT-PCR instrument; stop reaction before reaching amplification plateau to avoid PCR overcycling; in our experience, at 10 |
| Annealing/Extension/Data acquisition | 65°C | 45 s | |
| sgRNA amplicon PCR2 | ||
|---|---|---|
| Reagent | Final concentration / amount | Amount |
| Amp_sgRNA_2F primer (10 μM) | 1.5 μM | 1.5 μl |
| Amp_sgRNA_2R primer (10 μM) | 1.5 μM | 1.5 μl |
| 2× KAPA HiFi Master Mix (Kapa Biosystems, KR0389) | 1× | 5 μl |
| 1:25 dilution of purified PCR1 product from step 49d | n/a | 1 μl |
| Nuclease-free water | n/a | 1 μl |
| n/a | ||
sequence of forward primer (Table S1): GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAGCTTGTGGAAAGGACGAAACAC (primes on the U6 promoter adjacent to the sgRNA protospacer sequence and adds a standard Nextera Read 2 primer).
sequence of reverse primer (Table S1): AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTC (adds the remainder of TruSeq Read 1 and the adapter P5).
| PCR2 cycling conditions | |||
|---|---|---|---|
| Steps | Temperature | Time | Cycles |
| Initial Denaturation/ Enzyme activation | 95°C | 5 min | 1 |
| Denaturation | 95°C | 30 s | Variable (monitor amplification in the screen of the q-RT-PCR instrument; stop reaction before reaching amplification plateau to avoid PCR overcycling; in our experience, at 12–15 cycles) ( |
| Annealing/Extension/Data acquisition | 65°C | 45 s | |
| sgRNA amplicon PCR3 | ||
|---|---|---|
| Reagent | Final concentration / amount | Amount |
| Nextera XT P7 index (10 μM) | 1.5 μM | 1.5 μl |
| Amp_sgRNA_2R primer (10 μM) | 1.5 μM | 1.5 μl |
| 2× KAPA HiFi Master Mix (Kapa Biosystems, KR0389) | 1× | 5 μl |
| 1:25 dilution of purified PCR2 product from step 49h | n/a | 1 μl |
| Nuclease-free water | n/a | 1 μl |
| n/a | ||
sequence of forward primer: Nextera XT P7 index in the form of: CAAGCAGAAGACGGCATACGAGAT-[8 bp Index]-GTCTCGTGGGCTCGG. These adaptors can be obtained from the Illumina kits “Nextera XT Index Kit” (Illumina, FC-131-2001, FC-131-2002, FC-131-2003, FC-131-2004 or FC-131-1001).
sequence of reverse primer (Table S1): AATGATACGGCGACCACCGAGATCTACACTCTTTCCCTACACGACGCTC (same reverse primer from PCR2).
| PCR3 cycling conditions | |||
|---|---|---|---|
| Steps | Temperature | Time | Cycles |
| Initial Denaturation/ Enzyme activation | 95°C | 5 min | 1 |
| Denaturation | 95°C | 30 s | Variable (monitor amplification in the screen of the q-RT-PCR instrument; stop reaction before reaching amplification plateau to avoid PCR overcycling; in our experience, at 10 |
| Annealing/Extension/Data acquisition | 72°C | 45 s | |