- Open Access
Guide RNAs with embedded barcodes boost CRISPR-pooled screens
- Shiyou Zhu†1,
- Zhongzheng Cao†1, 2,
- Zhiheng Liu†1, 2,
- Yuan He†1,
- Yinan Wang1, 2,
- Pengfei Yuan3,
- Wei Li4, 5,
- Feng Tian1,
- Ying Bao1 and
- Wensheng Wei1Email authorView ORCID ID profile
© The Author(s). 2019
- Received: 10 December 2018
- Accepted: 11 January 2019
- Published: 24 January 2019
We report a new method using re-designed guide RNAs with internal barcodes (iBARs) embedded in their loop regions. Our iBAR approach outperforms the conventional method by producing screening results with much lower false-positive and false-negative rates especially with a high multiplicity of infection (MOI). Importantly, the iBAR approach reduces the starting cells at high MOI significantly with greatly improved efficiency and accuracy compared with the canonical CRISPR screens at a low MOI. This new system is particularly useful when the source of cells is limited or when it is difficult to control viral infection for in vivo screening.
The CRISPR/Cas9 system enables editing at targeted sites in the genome with high efficiency and specificity [1–3]. One of its extensive applications is to identify the functions of coding genes, non-coding RNAs, and regulatory elements through high-throughput pooled screening in combination with next-generation sequencing (NGS) analysis. By introducing a pooled single-gRNA (sgRNA) or paired-gRNA (pgRNA) library into cells expressing Cas9 or catalytically inactive Cas9 (dCas9) fused with effector domains, investigators can perform multifarious genetic screens by generating diverse mutations, large genomic deletions, transcriptional activation, or transcriptional repression [4–11]. To generate a high-quality cell library of gRNAs for any given CRISPR screen, one must use a low multiplicity of infection (MOI) during cell library construction to ensure that the majority of cells harbors one sgRNA/pgRNA to minimize the false-positive discovery rate [7, 12, 13]. To further reduce the FDR (false discovery rate) and increase data reproducibility, in-depth coverage of gRNAs and multiple biological replicates are often necessary to obtain the hit genes with high statistical significance , resulting in increased library size. Additional difficulties may arise when one performs a large number of genome-wide screens, when the cell materials for library construction are limited, or when one conducts more challenging screens (i.e., in vivo) for which it is difficult to arrange the experimental replications or control the MOI. To address these technical difficulties, we aimed to design a novel method that enabled us to use a high MOI to generate CRISPR libraries in target cells and to obtain improved statistics for the screening results. We hypothesized that both the false-positive and false-negative rates of high MOI screens could be drastically reduced if we had more experimental replicates for each of the gRNAs. If we assigned the gRNAs different barcodes that were somehow embedded within the gRNA sequences, we could then trace the performance of each gRNA multiple times within the same experiment by counting both the gRNAs and their internal barcode (iBAR) sequences.
Two papers have recently reported methods to generate random barcodes outside the sgRNA body for pooled CRISPR screening [14, 15]. Assuming each sgRNA would create both desired loss-of-function (LOF) and non-LOF alleles, calculating all reads of any given sgRNA is unable to accurately assess the importance of its targeting gene in negative screening. Much improved statistical results could be achieved by linking one UMI (unique molecular identifier) with one editing outcome of each sgRNA to enable single-cell lineage tracing so as to lower the false-negative rate , or by counting the decreased number of RSLs (random sequence labels) affiliated with sgRNAs to improve screening quality . Different from these two methods, we aimed to develop a method through embedded iBARs to enable pooled screening with CRISPR library made of viral infection at a high MOI, so as to reduce library size and improve data quality.
A design-based CRISPR sgRNAiBAR library screening method
Comparison of screenings at MOI of 0.3, 3, and 10 for essential genes involved in TcdB toxicity
We then constructed an sgRNAiBAR library covering every annotated human gene. For each of the 19,210 human genes, three unique sgRNAs were designed using our newly developed DeepRank method (“Methods” section), each of which was randomly assigned four iBAR6s. In addition, 1000 non-targeting sgRNAs, each with four iBAR6s, were included as negative controls (Additional file 2: Table S2). For the ease of statistical comparison, every set of three unique non-targeting sgRNAs was artificially named a negative control gene. The 85-nt sgRNAiBAR oligos were designed in silico (Additional file 1: Figure S3), synthesized using array synthesis, and cloned as a pooled library into a lentiviral backbone. Cas9-expressing HeLa cells were transduced with the sgRNAiBAR library lentivirus at three different MOIs (0.3, 3, and 10) with 400-fold coverage for sgRNAs to generate cell libraries, in which each sgRNAiBAR was covered 100-fold. To evaluate the effect of iBAR design for CRISPR screening at different MOIs, we performed a positive screening to identify genes that mediate the cytotoxicity of Clostridium difficile toxin B (TcdB), one of the key virulence factors of this anaerobic bacillus . We have previously reported the first identification of the functional receptor of TcdB, CSPG4 , whose coding gene was also identified and ranked at the very top from a genome-scale CRISPR library screening . In this reported CRISPR screening, UGP2 gene was also top-ranked hit, and FZD2 was identified and confirmed to encode the secondary receptor that mediates the TcdB’s killing effect on host cells. Of note, the role of FZD2 was significantly dwarfed by CSPG4 so that the FZD2 gene could only be identified using the truncated TcdB that had CSPG4-interacting region deleted . In our screens on TcdB, we used MAGeCKiBAR (Additional file 2: Table S3) and MAGeCK (Additional file 3: Table S4) to analyze data from iBAR and the conventional CRISPR screens, respectively. We consequently obtained top-ranked genes (FDR < 0.15) from both.
Identification of genes important for 6-TG-mediated cytotoxicity using the CRISPRiBAR and conventional CRISPR-pooled screens
To further evaluate the power of iBAR method, we went on conducting a screening to identify genes that modulate cellular susceptibility to 6-TG , a cancer drug that could be processed to inhibit DNA synthesis. We decided to construct the genome-scale sgRNAiBAR library at a MOI of 3 to generate a cell library with high coverage (2000-fold) for each sgRNA, in which each sgRNAiBAR was covered 500-fold. The overall read distribution of both experimental replicates was shown (Additional file 1: Figure S5a), and the reference cell libraries of both replicates reached 97% coverage of all originally designed sgRNAs (Additional file 1: Figure S5b). Over 95% of the sgRNAs in the original libraries retained three to four iBARs, indicating the good quality of libraries in which most sgRNAs had sufficient barcode variants for screening and data analysis (Additional file 1: Figure S5c). The fold change of all genes correlated well between the two biological replicates (Additional file 1: Figure S6). For the same 6-TG screening of two sgRNA library replicates, we also employed MAGeCK and MAGeCKiBAR analyses. For MAGeCKiBAR, we consequently obtained adjusted variance and mean distributions for all the sgRNAsiBAR that heightened the variance of sgRNAs with enrichment inconsistent among different iBAR repeats (Additional file 1: Figure S7).
To validate the screening results, we de novo designed and combined two sgRNAs to make a mini-pool to target each candidate gene, and each pool was introduced into HeLa cells through lentiviral infection (Additional file 3: Table S8). The effects of the sgRNA pools on cell viability against 6-TG treatment were quantified by a 3-(4,5-dimethyl-2-thiazolyl)-2,5-diphenyl-2-H-tetrazolium bromide (MTT) assay. Top 10 genes from CRISPRiBAR as well as CRISPR screens were chosen for validation. Noticeably, two non-targeting control genes were identified and ranked in the top 10 candidate list from the conventional CRISPR screen (Additional file 3: Table S6). These evident false-positive results are predictable because of the high MOI we used to generate the cell library. We successfully confirmed that the top 10 candidate genes from CRISPRiBAR of both replicates were all true-positive results; in contrast, only five genes from the top 10 candidate list from the conventional method turned out to be true positives (Fig. 4e). Among them, four genes (HPRT1, ITGB1, SRGAP2, and AKTIP) were obtained using both methods, whereas six genes (ACTR3C, PPP1R17, ACSBG1, CALM2, TCF21, and KIFAP3) were only identified and ranked at the top from CRISPRiBAR. In summary, iBAR improved accuracy with lower false-positive and false-negative rates for high MOI screens compared with conventional method.
We further assessed the performance of each sgRNAiBAR targeting the top four candidate genes (HPRT1, ITGB1, SRGAP2, and AKTIP). All the different iBARs of the enriched sgRNAs appeared to have little effect on the enrichment levels of their affiliated sgRNAs, and the order of iBARs associated with any particular sgRNA appeared to be random (Additional file 1: Figure S9), further supporting our prior notion that the iBARs did not affect the efficiency of their affiliated sgRNAs. All four HPRT1-targeting sgRNAsiBAR were significantly enriched after 6-TG treatment in both replicates (Fig. 4f). Most sgRNAsiBAR of other CRISPRiBAR identified genes were enriched after 6-TG selection (Additional file 1: Figure S10). In contrast, only a very few of sgRNAsiBAR of some top-ranked genes from conventional CRISPR screening were enriched, including FGF13 (Fig. 4g), GALR1, and two negative control genes (Additional file 1: Figure S11), leading to false-positive hits in the MAGeCK but not MAGeCKiBAR analysis (Additional file 1: Figure S12).
The library coverage was significantly increased with a high MOI in the transduction with a fixed number of cells for library construction, so the starting cells for library construction could be decreased more than 10-fold (MOI = 3) and 35-fold (MOI = 10) to match and even top the results from conventional screening at an MOI of 0.3 (Additional file 3: Table S10). Four barcodes for each sgRNA, as we designed, appeared to provide sufficient internal repeats to enable the high level of consistency between the two biological replicates using the iBAR method (Fig. 3, Fig. 4d, Additional file 3: Table S7). Therefore, in addition to the significant reduction in cells for library construction, the internal replicates offered by iBARs within the same experiment would lead to more uniform conditions and fairer comparisons versus separate biological replicates, consequently improving statistical scores. The advantage of the iBAR method would become greater when large-scale CRISPR screens in multiple cell lines are in demand or when the cell samples for screening are scarce (i.e., samples from patients or those of primary origin). Especially for in vivo screening in which the lentiviral transduction rate is hard to predict and variable conditions in different animals might greatly impact the screening outcomes, the iBAR method could be an ideal solution to resolve these technical limitations.
Because multiple cuttings decrease cell viability, CRISPR library constructed at a high MOI might have abnormal false discovery rate for negative screening [24, 25]. We therefore performed a genome-scale negative screening at an MOI of 0.3 to assess the iBAR method in calling essential genes. For positive screening using iBAR, we modified the model-estimated variance of sgRNAs with different fold change directions among barcodes to enlarge variance so that the mis-associated sgRNAs were subjected to adequate penalty. For negative screening, however, sgRNA depletion through mis-association had little effect on its consistency of fold change directions as non-functional sgRNAs remained unchanged. Therefore, we treated barcodes only as internal replicates without the penalty procedure. We indeed achieved improved statistics with higher true-positive and lower false-positive rates for negative screening using iBAR method at a low MOI than the conventional approach using gold standard essential genes  (Additional file 1: Figure S13).
Notwithstanding the technical advancement of the iBAR method to offer the same benefit of internal replications, we must be cautious with the MOI during viral transduction to generate the original cell library in negative screens based on measuring cell viability. Although massive integrates have been reported not to affect cell fitness , multiple cuttings on DNA caused by higher MOI in cells with active Cas9 have been shown to reduce cell viability [24, 25]. Strategies without cuttings, such as CRISPRi/a  or iSTOP systems , could be better choices to combine with the iBAR system for negative screening at a high MOI.
Although we had data to support that iBAR6 had little effect on the activities of sgRNAs, we would not recommend to use barcodes with consecutive T (> 4) so as to avoid any minor effects. Ultimately, 4096 types of iBAR6 provided sufficient varieties to make CRISPR libraries. In addition, the length of the iBAR is not limited to 6 nt. We have tested different lengths of iBARs and found that their lengths could be up to 50 nt without affecting functions of their affiliated sgRNAs.
(Additional file 1: Figure S14). In addition, it is not necessary to design different barcode sets for different sgRNAs. A fixed set of iBARs assigned to all sgRNAs should work as well as random assignment in library screening. Our iBAR strategy with a streamlined analytic tool MAGeCKiBAR would facilitate large-scale CRISPR screens for broad biomedical discoveries in various settings.
Cells and reagents
HeLa and HEK293T cell lines were maintained in Dulbecco’s modified Eagle’s medium (DMEM; Gibco C11995500BT) supplemented with 1% penicillin/streptomycin and 10% fetal bovine serum (FBS; CellMax BL102-02) and cultured with 5% CO2 at 37 °C. All cells were checked for the absence of mycoplasma contamination.
The lentiviral sgRNAiBAR-expressing backbone was constructed by changing the position of the BsmBI (Thermo Scientific™, ER0451) site using BstBI (NEB, R0519) and XhoI (NEB, R0146) from Plenti-sgRNA-Lib (Addgene, #53121). sgRNA- and sgRNAiBAR-expressing sequences were cloned into the backbone using the BsmBI-mediated Golden Gate cloning strategy .
Design of the genome-scale CRISPR sgRNAiBAR library
Gene annotations were retrieved from the UCSC hg38 genome, which contains 19,210 genes. For each gene, three different sgRNAs that had at least one mismatch in the 16-bp seed region in the genome with a high level of predicted targeting efficiency were designed using our newly developed DeepRank algorithm. We then randomly assigned four 6-bp iBARs (iBAR6s) to each sgRNA. We designed an additional 1000 non-targeting sgRNAs, each with four iBAR6s, to serve as negative controls.
Construction of the CRISPR sgRNAiBAR plasmid library
The 85-nt DNA oligonucleotides were designed and array synthesized. Primers (oligo-F and oligo-R) targeting the flanking sequences of oligos were used for PCR amplification. The PCR products were cloned into the lentiviral vector constructed above using the Golden Gate method . The ligation mixtures were transformed into Trans1-T1 competent cells (Transgene, CD501-03) to obtain library plasmids. Transformed clones were counted to ensure at least 100-fold coverage for the scale of the sgRNAiBAR library. The library plasmids were extracted following the standard protocol (QIAGEN 12362) and transfected into HEK293T cells with the two lentivirus package plasmids pVSVG and pR8.74 (Addgene, Inc.) to obtain the library virus. The iBAR library containing all 4096 iBAR6s for one ANTXR1-targeting sgRNA was constructed using the same protocol.
Screening of the sgRNAiBAR-ANTXR1 library containing all 4096 types of iBAR6
A total of 2 × 107 cells were plated on 150-mm Petri dishes and infected with the library lentivirus at an MOI of 0.3. After 72 h of infection, cells were re-seeded and treated with 1 μg/ml of puromycin (Solarbio P8230) for 48 h. For each replicate, 5 × 106 cells were collected for genome extraction. Screening of the sgRNAiBAR-ANTXR1 library was performed using PA/LFnDTA toxin [30, 31] after library-infected cells were cultured for 15 days . Then, sgRNA with the iBAR coding region in genomic DNA was amplified (TransGen, AP131-13) using Primer-F and Primer-R and then subjected to high-throughput sequencing analysis (Illumina HiSeq2500) using an NEBNext Ultra DNA Library Prep Kit for Illumina (NEB E7370L).
Screening of the genome-scale CRISPR/Cas9 sgRNAiBAR library for genes important for TcdB cytotoxicity and for genes essential for cell viability
A total of 1.6 × 108 cells (MOI = 0.3), 1.53 × 107 cells (MOI = 3), and 4.6 × 106 cells (MOI = 10) were plated on 150-mm Petri dishes respectively for sgRNA library construction for two replicates. Cells were infected with the library lentivirus of different MOIs and treated with 1 μg/ml of puromycin for 72 h post infection. sgRNAiBAR-integrated cells were cultured for an additional 15 days to maximize gene knockout. Cells were re-seeded onto 150-mm Petri dishes, treated by TcdB (100 pg/ml) for 10 h, and followed by the removal of the loosely attached round cells through repeated pipetting . For each round of screening, the cells were cultured in fresh medium without TcdB to reach ~ 50–60% confluence. All resistant cells in one replicate were pooled and subjected to another round of TcdB screening. For the subsequent three rounds of screening, the TcdB concentration was 125 pg/ml, 150 pg/ml, and 175 pg/ml, respectively. After four rounds of treatment, the resistant cells and untreated cells were collected for genomic DNA extraction, amplification of sgRNA, and NGS analysis. Seven pairs of primers were used for PCR amplification (Additional file 3: Table S9), and PCR products were mixed for NGS. For negative screening at an MOI of 0.3, a total of 4.6 × 107 (two replicates) sgRNAiBAR-integrated cells were cultured for 28 days before NGS decoding.
Screening of the genome-scale CRISPR/Cas9 sgRNAiBAR library for genes important for 6-TG cytotoxicity
A total of 5 × 107 cells were plated on 150-mm Petri dishes, and two replicates were obtained. Cells were infected with the library lentivirus at an MOI of 3 and treated with 1 μg/ml puromycin 72 h after infection. sgRNAiBAR-integrated cells were cultured for an additional 15 days, re-seeded at a total number of 5 × 107, and then treated with 200 ng/ml 6-TG (Selleck). For the following two rounds of screening, the 6-TG concentration was 250 ng/ml and 300 ng/ml. For each round of selection, the drug was maintained for 7 days, and the cells were cultured in fresh medium without 6-TG for another 3 days. Then, all the resistant cells in one replicate were grouped together and subjected to another round of 6-TG screening. After three rounds of treatment, the resistant cells and untreated cells were collected for genomic DNA extraction, amplification of sgRNA with iBAR regions, and deep-sequencing analysis.
Positive screening data analysis
MAGeCKiBAR is the analysis strategy developed for screens using an sgRNAiBAR library based on the MAGeCK algorithm . MAGeCKiBAR takes great advantage of Python, Pandas, NumPy, and SciPy. The analysis algorithm contains three main parts: analysis preparation, statistical tests, and rank aggregation. In the analysis preparation stage, the inputted raw counts of sgRNAsiBAR are normalized, and the coefficients of the population mean and variance are then modeled. In the statistical test stage, we use tests to determine the significance of the difference between the treatment and control normalized reads. In the rank aggregation stage, we aggregate the ranks of all the sgRNAsiBAR targeting each gene to obtain the final gene ranking.
Normalization and preparation
The means of the treatment and control counts were calculated directly, and the corresponding variance could be calculated from the mean and coefficients. For the CRISPR-iBAR analysis, we evaluated the enrichment of sgRNAs through the performances of different iBARs. We designed four iBARs for each sgRNA to serve as internal replicates. Due to the high MOI during library construction, there must be free riders of false-positive sgRNAs associated with true-positive hits. The free rider here was used to describe the sgRNAs targeting irrelevant genes that were mis-associated with functional sgRNAs to enter the same cells. We modified the variance of sgRNAsiBAR based on the enrichment directions of different iBARs for each sgRNA. If all the iBARs of one sgRNA presented the same direction of fold change, i.e., all greater or less than that of the control group, then the variance would be unchanged. However, if one sgRNA with different iBARs revealed inconsistent directions of fold change, then this kind of sgRNA would be penalized by increasing its variance. The final adjusted variance for inconsistent sgRNAsiBAR would be the model-estimated variance plus the experimental variance calculated from the Ctrl and Exp samples.
where ti is the mean of the treatment counts of the ith sgRNA and ci and vi are the mean and variance of control counts of the ith sgRNA. Because the variance is used as the denominator to calculate the score, the enlarged variance for the inconsistent sgRNAsiBAR results in lower score.
Statistical test and rank aggregation
The normal distribution was used to test the scorei of the treatment counts. The two sides of scores in a standard normal distribution provided the greater-tail and lesser-tail P value separately. To obtain the gene ranks, we used RRA, which is an appropriate method for aggregating rankings . MAGeCK adopted a modified RRA method by limiting the enriched sgRNAs . Suppose for one gene there are n sgRNAs with different iBARs in the library of M sgRNAsiBAR in total; every sgRNAiBAR has a rank in the library of R = (R1, R2, … , Rn). First, the ranks of sgRNAsiBAR should be normalized by the total number of sgRNAsiBAR in the library. We obtained the normalized rank r = (r, r2, … , rn) for each ri = Ri/M, in which 1 ≤ i ≤ n. Then, we calculated the sorted normalized ranking sr, making sr1 ≤ sr2 ≤ … ≤ srn. The sorted normalized rank follows a uniform distribution between 0 and 1. The probability βk, n(sr) in which sri ≤ ri follows a β distribution β(k, n + 1 − k), making ρ = min(β1, n, β2, n, … , βn, n). For every gene, the ρ score can be obtained by RRA and further adjusted by Bonferroni correction . We adopted MAGeCK, which developed α-RRA, to select the top α% sgRNAs from the ranking list. The P values of sgRNAs lower than a threshold (0.25 for instance) were selected. Only the top sgRNAs of one gene were considered in the RRA calculation, thus making ρ = min(β1, n, β2, n, … , βj, n), in which 1 ≤ j ≤ n.
Negative screening data analysis
During the analyzing process of positive screening at high MOI based on the iBAR strategy, we modified the model-estimated variance of sgRNAs with different fold change directions among corresponding barcodes. But for negative screening, most of the non-functional sgRNAs would be unchanged. So the variance modification algorithm based on fold change directions of corresponding barcodes becomes not sufficient to justify weather certain sgRNA is a false-positive result. Therefore, we treated barcodes as internal replicates directly. When taking iBAR into consideration, we performed two times robust rank aggregation for the negative screening rather than variance adjustment for the inconsistent sgRNAsiBAR. The first round of robust rank aggregation aggregates the sgRNAiBAR level to sgRNA level, and the second round aggregates the sgRNA level to gene level.
Validation of candidate genes
To validate each gene, we chose two sgRNAs designed in the library and cloned into a lentiviral vector with a puromycin selection marker. We mixed two sgRNA plasmids and co-transfected them into HEK293T cells with two lentiviral package plasmids (pVSVG and pR8.74) using the X-tremeGENE HP DNA transfection reagent (Roche). The HeLa cells stably expressing Cas9 were infected with the lentivirus for 3 days and treated with 1 μg/ml puromycin for 2 days. Then, 5000 cells were added into each well, and five replicates were obtained for each group. After 24 h, the experimental groups were treated with 150 ng/ml 6-TG, and the control groups were treated with normal medium for 7 days. Then, MTT (Amresco) staining and detection were performed following the standard protocol. The experimental wells treated with 6-TG were normalized to the wells without 6-TG treatment.
We acknowledge the staff of the BIOPIC sequencing facility (Peking University) for their assistance. The process of data analysis was supported by High-performance Computing Platform of Peking University. We thank the support of National Post-doctoral Innovative Talent Supporting Program.
This project was supported by funds from the National Science Foundation of China (NSFC31430025), the Beijing Advanced Innovation Center for Genomics at Peking University, and the Peking-Tsinghua Center for Life Sciences (to W.W.); the National Major Science & Technology Project for Control and Prevention of Major Infectious Diseases in China (2018ZX10301401, to S.Z.).
Availability of data and materials
The accession number for the raw sequencing data reported in this paper is NCBI Sequence Read Archive (SRA): SRP164294 (PRJNA494953) . The MAGeCK-based analysis program named as MAGeCKiBAR written by python3 for iBAR based CRISPR screening is available at https://bitbucket.org/WeiLab/mageck-ibar/. 
WW conceived and supervised the project. WW, SZ, ZC, and PY designed the experiments. SZ, ZC, and YH performed the experiments with the help from FT and YB. YW designed the oligos used to construct the iBAR library. ZL performed the data analysis with the help of WL. SZ, ZC, ZL, YH, and WW wrote the manuscript with the help of the other authors. All authors read and approved the final manuscript.
Ethics approval and consent to participate
Not applicable for this study.
Consent for publication
Not applicable for this study.
A patent has been filed relating to the data presented, which would not limit the academic use of the data and the method. W.W. is a founder and scientific advisor for EdiGene. All other authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Jinek M, Chylinski K, Fonfara I, Hauer M, Doudna JA, Charpentier E. A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. Science. 2012;337:816–21.View ArticleGoogle Scholar
- Cong L, Ran FA, Cox D, Lin S, Barretto R, Habib N, Hsu PD, Wu X, Jiang W, Marraffini LA, Zhang F. Multiplex genome engineering using CRISPR/Cas systems. Science. 2013;339:819–23.View ArticleGoogle Scholar
- Mali P, Yang L, Esvelt KM, Aach J, Guell M, DiCarlo JE, Norville JE, Church GM. RNA-guided human genome engineering via Cas9. Science. 2013;339:823–6.View ArticleGoogle Scholar
- Shalem O, Sanjana NE, Hartenian E, Shi X, Scott DA, Mikkelson T, Heckl D, Ebert BL, Root DE, Doench JG, Zhang F. Genome-scale CRISPR-Cas9 knockout screening in human cells. Science. 2014;343:84–7.View ArticleGoogle Scholar
- Wang T, Wei JJ, Sabatini DM, Lander ES. Genetic screens in human cells using the CRISPR-Cas9 system. Science. 2014;343:80–4.View ArticleGoogle Scholar
- Koike-Yusa H, Li Y, Tan EP, Velasco-Herrera Mdel C, Yusa K. Genome-wide recessive genetic screening in mammalian cells with a lentiviral CRISPR-guide RNA library. Nat Biotechnol. 2014;32:267–73.View ArticleGoogle Scholar
- Zhou Y, Zhu S, Cai C, Yuan P, Li C, Huang Y, Wei W. High-throughput screening of a CRISPR/Cas9 library for functional genomics in human cells. Nature. 2014;509:487–91.View ArticleGoogle Scholar
- Zhu S, Li W, Liu J, Chen CH, Liao Q, Xu P, Xu H, Xiao T, Cao Z, Peng J, et al. Genome-scale deletion screening of human long non-coding RNAs using a paired-guide RNA CRISPR-Cas9 library. Nat Biotechnol. 2016;34:1279–86.View ArticleGoogle Scholar
- Liu Y, Cao Z, Wang Y, Guo Y, Xu P, Yuan P, Liu Z, He Y, Wei W. Genome-wide screening for functional long noncoding RNAs in human cells by Cas9 targeting of splice sites. Nat Biotechnol. 2018;36:1203–10.View ArticleGoogle Scholar
- Gilbert LA, Horlbeck MA, Adamson B, Villalta JE, Chen Y, Whitehead EH, Guimaraes C, Panning B, Ploegh HL, Bassik MC, et al. Genome-scale CRISPR-mediated control of gene repression and activation. Cell. 2014;159:647–61.View ArticleGoogle Scholar
- Konermann S, Brigham MD, Trevino AE, Joung J, Abudayyeh OO, Barcena C, Hsu PD, Habib N, Gootenberg JS, Nishimasu H, et al. Genome-scale transcriptional activation by an engineered CRISPR-Cas9 complex. Nature. 2015;517:583–8.View ArticleGoogle Scholar
- Peng J, Zhou Y, Zhu S, Wei W. High-throughput screens in mammalian cells using the CRISPR-Cas9 system. FEBS J. 2015;282:2089–96.View ArticleGoogle Scholar
- Zhu S, Zhou Y, Wei W. Genome-wide CRISPR/Cas9 screening for high-throughput functional genomics in human cells. Methods Mol Biol. 2017;1656:175–81.View ArticleGoogle Scholar
- Michlits G, Hubmann M, Wu SH, Vainorius G, Budusan E, Zhuk S, Burkard TR, Novatchkova M, Aichinger M, Lu Y, et al. CRISPR-UMI: single-cell lineage tracing of pooled CRISPR-Cas9 screens. Nat Methods. 2017;14:1191–7.View ArticleGoogle Scholar
- Schmierer B, Botla SK, Zhang J, Turunen M, Kivioja T, Taipale J. CRISPR/Cas9 screening using unique molecular identifiers. Mol Syst Biol. 2017;13:945.View ArticleGoogle Scholar
- Shechner DM, Hacisuleyman E, Younger ST, Rinn JL. Multiplexable, locus-specific targeting of long RNAs with CRISPR-Display. Nat Methods. 2015;12:664–70.View ArticleGoogle Scholar
- Bradley KA, Mogridge J, Mourez M, Collier RJ, Young JA. Identification of the cellular receptor for anthrax toxin. Nature. 2001;414:225–9.View ArticleGoogle Scholar
- Li W, Xu H, Xiao T, Cong L, Love MI, Zhang F, Irizarry RA, Liu JS, Brown M, Liu XS. MAGeCK enables robust identification of essential genes from genome-scale CRISPR/Cas9 knockout screens. Genome Biol. 2014;15:554.View ArticleGoogle Scholar
- Lyras D, O'Connor JR, Howarth PM, Sambol SP, Carter GP, Phumoonna T, Poon R, Adams V, Vedantam G, Johnson S, et al. Toxin B is essential for virulence of Clostridium difficile. Nature. 2009;458:1176–9.View ArticleGoogle Scholar
- Yuan P, Zhang H, Cai C, Zhu S, Zhou Y, Yang X, He R, Li C, Guo S, Li S, et al. Chondroitin sulfate proteoglycan 4 functions as the cellular receptor for Clostridium difficile toxin B. Cell Res. 2015;25:157–68.View ArticleGoogle Scholar
- Tao L, Zhang J, Meraner P, Tovaglieri A, Wu X, Gerhard R, Zhang X, Stallcup WB, Miao J, He X, et al. Frizzled proteins are colonic epithelial receptors for C. difficile toxin B. Nature. 2016;538:350–5.View ArticleGoogle Scholar
- Tan YY, Epstein LB, Armstrong RD. In vitro evaluation of 6-thioguanine and alpha-interferon as a therapeutic combination in HL-60 and natural killer cells. Cancer Res. 1989;49:4431–4.PubMedGoogle Scholar
- Duan J, Nilsson L, Lambert B. Structural and functional analysis of mutations at the human hypoxanthine phosphoribosyl transferase (HPRT1) locus. Hum Mutat. 2004;23:599–611.View ArticleGoogle Scholar
- Jackson SP. Sensing and repairing DNA double-strand breaks. Carcinogenesis. 2002;23:687–96.View ArticleGoogle Scholar
- Meyers RM, Bryan JG, McFarland JM, Weir BA, Sizemore AE, Xu H, Dharia NV, Montgomery PG, Cowley GS, Pantel S, et al. Computational correction of copy number effect improves specificity of CRISPR-Cas9 essentiality screens in cancer cells. Nat Genet. 2017;49:1779–84.View ArticleGoogle Scholar
- Hart T, Brown KR, Sircoulomb F, Rottapel R, Moffat J. Measuring error rates in genomic perturbation screens: gold standards for human functional genomics. Mol Syst Biol. 2014;10:733.View ArticleGoogle Scholar
- Zhou Y, Wang P, Tian F, Gao G, Huang L, Wei W, Xie XS. Painting a specific chromosome with CRISPR/Cas9 for live-cell imaging. Cell Res. 2017;27:298–301.View ArticleGoogle Scholar
- Billon P, Bryant EE, Joseph SA, Nambiar TS, Hayward SB, Rothstein R, Ciccia A. CRISPR-mediated base editing enables efficient disruption of eukaryotic genes through induction of STOP codons. Mol Cell. 2017;67:1068–79 e1064.View ArticleGoogle Scholar
- Engler C, Gruetzner R, Kandzia R, Marillonnet S. Golden gate shuffling: a one-pot DNA shuffling method based on type IIs restriction enzymes. PLoS One. 2009;4:e5553.View ArticleGoogle Scholar
- Wei W, Lu Q, Chaudry GJ, Leppla SH, Cohen SN. The LDL receptor-related protein LRP6 mediates internalization and lethality of anthrax toxin. Cell. 2006;124:1141–54.View ArticleGoogle Scholar
- Qian L, Cai C, Yuan P, Jeong SY, Yang X, Dealmeida V, Ernst J, Costa M, Cohen SN, Wei W. Bidirectional effect of Wnt signaling antagonist DKK1 on the modulation of anthrax toxin uptake. Sci China Life Sci. 2014;57:469–81.View ArticleGoogle Scholar
- Anders S, Huber W. Differential expression analysis for sequence count data. Genome Biol. 2010;11:R106.View ArticleGoogle Scholar
- Robinson MD, Smyth GK. Small-sample estimation of negative binomial dispersion, with applications to SAGE data. Biostatistics. 2008;9:321–32.View ArticleGoogle Scholar
- Kolde R, Laur S, Adler P, Vilo J. Robust rank aggregation for gene list integration and meta-analysis. Bioinformatics. 2012;28:573–80.View ArticleGoogle Scholar
- Zhu S, Cao Z, Liu Z, He Y, Wang Y, Yuan P, Li W, Tian F, Bao Y, Wei W. Guide RNAs with embedded barcodes boost CRISPR-pooled screens. Sequence Read Archive (SRA) 2019. Dataset. https://www.ncbi.nlm.nih.gov/sra/SRP164294
- Zhu S, Cao Z, Liu Z, He Y, Wang Y, Yuan P, Li W, Tian F, Bao Y, Wei W. Guide RNAs with embedded barcodes boost CRISPR-pooled screens. Bitbucket Software. https://bitbucket.org/WeiLab/mageck-ibar/.