key: cord-0939710-xvi0asl9 authors: Zhang, Beibei; Zhou, Jiamu; Li, Miao; Wei, Yuanmeng; Wang, Jiaojiao; Wang, Yange; Shi, Pingling; Li, Xiaoli; Huang, Zixu; Tang, He; Song, Zongming title: Evaluation of CRISPR/Cas9 site-specific function and validation of sgRNA sequence by a Cas9/sgRNA-assisted reverse PCR technique date: 2021-03-04 journal: Anal Bioanal Chem DOI: 10.1007/s00216-021-03173-2 sha: c15585dcbe1c12f59f97ff48fb60bfd8665d5a72 doc_id: 939710 cord_uid: xvi0asl9 The effective application of the clustered regularly interspaced short palindromic repeat (CRISPR)/Cas9 system in biology, medicine and other fields is hindered by the off-target effects and loci-affinity of Cas9-sgRNA, especially at a genome-wide scale. In order to eliminate the occurrence of off-target effects and evaluate loci-affinity by CRISPR/Cas9 site-specific detection and screening of high-affinity sgRNA sequences, respectively, we develop a CRISPR/Cas9-assisted reverse PCR method for site-specific detection and sgRNA sequence validation. The detection method based on PCR can be used directly in the laboratory with PCR reaction conditions, without the need for an additional detection system, and the whole process of detection can be completed within 2 h. Therefore, it can be easily popularized with a PCR instrument. Finally, this method is fully verified by detecting multiple forms of site mutations and evaluating the affinity of a variety of sgRNA sequences for the CRISPR/Cas9 system. In sum, it provides an effective new analysis tool for CRISPR/Cas9 genome editing-related research. [Figure: see text] SUPPLEMENTARY INFORMATION: The online version contains supplementary material available at 10.1007/s00216-021-03173-2. CRISPR is widely found in the genomes of bacteria and archaea [1] , as was clearly defined in 2002 [2] . It has dramatically transformed our ability to edit the genomes of various organisms, especially in eukaryotes and animal models of human disease [3, 4] . With continuous development, in recent years CRISPR/Cas has become widely used in disciplines from basic biology to biotechnology and medicine [5] [6] [7] . Genome editing of mammalian cells was accomplished by the CRISPR/Cas9 system in 2013 [8, 9] . In 2014, genome-wide function screens in human cells were achieved [10, 11] . While the CRISPR/Cas system has been applied mainly to the study of basic biology, in 2017 Gootenberg et al. successfully engineered the system for innovative applications in nucleic acid detection [12] , and numerous diverse detection techniques associated with CRISPR/ Cas systems were subsequently developed for exploring different research fields [6, [13] [14] [15] [16] . Several researchers have explored treatments for various diseases through CRISPR/Cas systems [17] [18] [19] [20] [21] [22] [23] [24] . However, despite the invaluable advantages of the CRISPR/Cas system (for example, easily designed and highly effective) and its wide application in biological and medical research, there are some limitations that need to be overcome. The effective application of CRISPR/Cas systems in medicine and other fields is hindered by off-target effects, especially on a genome-wide scale [25, 26] . When an expected target DNA is modified with the CRISPR/Cas system, unpredictable genome damage is simultaneously triggered by offtarget effects. When used in clinical trials, this can cause unanticipated risk to the health of patients. In order to eliminate or reduce the occurrence of off-target effects, several researchers have extensively explored methods for optimizing the specificity of the CRISPR/Cas system. Labun et al. designed a web tool (CHOPCHOP) for CRISPR genome engineering which is able to avoid off-target effects by identifying single-guide RNA (sgRNA) targets [27] . Even though some sgRNA selection tools contain an extensive database [28] [29] [30] , the database is still insufficient for all genomic loci of various organisms. In addition, variations have been noted in the cleavage efficiency of the CRISPR/Cas system with different genomic loci [31] , and these variances cannot be estimated with current sgRNA selection tools. In order to evaluate and validate a sitespecific sgRNA sequence, an electrochemiluminescence assay method was developed by Liu et al. [32] for detecting nucleic acid with multiple site mutations in vitro, and a label-free CRISPR/Cas9 assay method based on CuNPs was established by Hu et al. [33] . sgRNA sequence validation and site-specific detection can be achieved with these two methods. However, they have not achieved the popularity of PCR-based methods. Polymerase chain reaction (PCR) was introduced in 1985 [34] , and as the most common nucleic acid amplification method, it has been widely used in almost all laboratories related to biology and medicine. Although there are many novel detection methods which have certain advantages, PCR is still the most commonly used. Detection methods improved on the basis of PCR can be used directly in the laboratory with PCR reaction conditions, without an additional detection system. Therefore, PCR-based methods are easily popularized. We have developed a CRISPR-or Cas9/ sgRNA-assisted reverse PCR (CARP) method to detect and type target DNAs such as the various genotypes of the human papillomavirus [35] . CARP, derived from Cas9-sgRNA, has high specificity, so it is ideally suited for site-specific detection of CRISPR systems. This paper demonstrates that CARP can be used for the evaluation of CRISPR/Cas9 site-specific function and screening of high-efficiency sgRNA before delivering Cas9-sgRNA vectors into cells, animal models and patients. Firstly, the single-base mutation detection ability of CARP was verified. Then we explored the effect of mutations at different sites of target DNA and sgRNA on the cleavage efficiency of Cas9-sgRNA. Finally, the effects of different sgRNAs on the cleavage efficiency of Cas9 were successfully evaluated by CARP. Thus, in this study, we developed a sensitive, cost-effective and time-saving CRISPR-assisted reverse PCR method for site-specific detection and sgRNA sequence validation. All oligonucleotides used in this study were synthesized by Sangon Biotechnology Co. Ltd. (Shanghai). sgRNA templates were generated in three steps (a complementary extension step and two fusion PCR amplification steps using primers listed in Table 1 ). Firstly, F1 and R were annealed and extended by PCR reaction for 10 cycles. A second PCR amplification step was then performed with primers F2 and sgR using the products of the previous step as template for 30 cycles. Thirdly, the last PCR reaction was performed with primers F3 and sgR using the products of the previous step as template for 30 cycles, after which the products of the third PCR were purified with a PCR purification kit. Finally, sgRNAs were synthesized by in vitro transcription. The transcription reaction was performed by incubating the purified products with RNA transcription reagents (T7 RNA polymerase, buffer and rNTPs, New England Biolabs) overnight at 37°C. The transcription products were mixed with TRIzol solution, and then successively purified with extraction of chloroform and isopropanol, and precipitated with ethanol. Purified RNA was dissolved in RNase-free ddH 2 O, quantified by spectrometry and detected by polyacrylamide gel electrophoresis. The sgRNA was synthesized by in vitro transcription according to the manufacturer's instructions. The CDS of vascular endothelial growth factor A (VEGFA, NCBI Reference Sequence: NM_001171628.1) was synthesized and integrated into pUC57 plasmid by Sangon Biotechnology Co. Ltd. All substrate DNAs (one nonmismatched sequence and 10 mutations) were amplified with the universal primers UP1 and UP2 (Table 1) , which are complementary with the pUC57 plasmid sequence. The 20 μL PCR reaction system consisted of 10 μL 2× premix Taq (Takara), 500 nM UP1, 500 nM UP2, and 10 ng recombinant plasmid DNA. The PCR procedure was as follows: 95°C for 5 min; 30 cycles of 95°C for 30 s, 60°C for 20 s, and 72°C for 60 s; and 72°C for 5 min. The PCR products were run with 1% agarose gel and recovered with a DNA Gel Extraction Kit. The purified DNA fragments were quantified with NanoDrop (Thermo) and used as the substrate DNAs for the CARP assay. All substrate DNAs were cleaved with a pair of Cas9/sgRNA complexes. The cleavage reaction, which consisted of 1 μM Cas9 Nuclease (NEB), 1× Cas9 Nuclease Reaction Buffer, 300 nM sgRNA a (a1 or a1m1; Table 1 ), and 300 nM sgRNA b (b1 or b1m1; Table 1 ), was pre-incubated for assembling Cas9/ sgRNA complexes (10 min at 25°C). The non-mismatched substrate DNAs were separately mixed with the pre-incubated product and cleaved at 37°C for 5 min. Then the cleavage product was linked with T4 Ligase. The ligation reaction reagents (15 μL) were as follows: 5 U T4 DNA Ligase (Thermo), 1× T4 Ligase Buffer, and 5 μL cleavage product. The ligation reaction was incubated at 22°C for 10 min. Finally, 1 μL ligation product was mixed with 10 μL 2× premix Taq (Takara), 500 nM RP1 and 500 nM RP2 for a traditional PCR (PCR) reaction. The PCR procedure was as follows: 95°C for 5 min; 30 cycles of 95°C for 20 s, 60°C for 20 s, and 72°C for 30s; and 72°C for 5 min. The PCR product was detected with 2% agarose gel electrophoresis. The cleavage reaction, which consisted of 1 μM Cas9 Nuclease (NEB), 1× Cas9 Nuclease Reaction Buffer, 300 nM sgRNA a (a1m1 to a1m5, a1 to a9; only one of these sgRNA was added in each reaction), and 300 nM sgRNA b (b1m1 to b1m5, b1 to b9; only one of these sgRNA was added in each reaction), was pre-incubated for assembling Cas9/ sgRNA complexes (10 min at 25°C), along with various negative control reactions. The substrate DNAs (gradient detection: 20, 2, 0.2, 0.02, 0.002, 0 ng; other related reactions were all 20 ng) were separately mixed with pre-incubated product and cleaved at 37°C for 5 min. Then the cleavage product was linked with T4 Ligase. The ligation reaction reagents (15 μL) were as follows: 5 U T4 DNA Ligase (Thermo), 1× T4 Ligase Buffer, and 5 μL cleavage product. The ligation reaction was incubated at 22°C for 10 min. Finally, 1 μL ligation product was mixed with 10 μL 2× SYBR Green Master Mix (Yeasen), 500 nM RP1 and 500 nM RP2 for quantitative PCR (qPCR) reaction. The qPCR procedure was as follows: 95°C for 10 min; 40 cycles of 95°C for 15 s; and 60°C for 1 min. The reactions were carried out on a StepOnePlus real-time PCR apparatus (ABI). The CARP experimental process is schematically shown in Fig. 1 . There are three main steps in CARP: (1) target DNA is specifically cut by a pair of Cas9/sgRNA complexes; (2) the cleaved DNA is rapidly linked by T4 DNA ligase; (3) the ligated DNA is efficiently amplified with PCR (PCR or qPCR). The quality of sgRNAs is crucial in this study. sgRNA templates were generated by a complementary extension step and two fusion PCR amplification steps. In the first step, two complementary chains were annealed and extended to generate a general sgRNA frame template. In the second step, 20 bp specific sequences complemented with target DNA were added to the general sgRNA frame template by fusion PCR. In the last step, T7 promoter sequences were added to the amplification products of the previous step by fusion PCR. The final PCR products were used as templates to synthesize sgRNAs. The products and sgRNAs for all three steps were respectively detected by agarose gel electrophoresis (see Supplementary Information [ESM] Fig. S1 ) and polyacrylamide gel electrophoresis (ESM Fig. S2) . The results showed that 20 bp specific sequences and T7 promoter sequences were successively added to the general sgRNA frame template, and complete sgRNA templates were successfully generated. Thus, a good foundation for the whole research was established. After sgRNAs were generated, we explored whether single-base mutations could be discriminated using CARP. Firstly, target DNA was cut by a pair of Cas9/sgRNA complexes. Then cleavage products were linked by T4 DNA Ligase. Finally, the ligated products were detected with PCR and qPCR. If target DNA was cleaved and linked successfully by a pair of Cas9/sgRNA complexes and T4 DNA Ligase, the PCR amplification was carried out normally, and there would be a 144 bp band in the electrophoretic lane. On the contrary, there were no right bands in lanes of various control groups (Fig. 2a) . When higher concentrations of target DNA were detected by PCR-based CARP, the bands were all in the right place (Fig. 2b) . It preliminarily proved that CARP has sufficient specificity to detect the single-base mutation. In order to further validate the specificity of CARP, the single-base mutation was detected by qPCR-based CARP (Fig. 3) . Firstly, the single-base mutation was perfectly discriminated by comparing the Ct values of different qPCR reactions (Fig. 3b and c) , and there was a negative correlation between the number of ligated products and Ct value. Then the sensitivity of CARP was determined by detecting 20, 2, 0.2, 0.02, 0.002 and 0 ng target DNA, respectively. It was shown that as few as 0.02 ng target DNA could be detected (Fig. 3d) , compared with previous methods (limit of detection of 30 nM and 0.13 nM reported by Liu et al. and Hu et al., respectively [32, 33] , whereas our substrate DNA is 1766 bp, and the total volume of the reaction is 30 μL, so 0.02 ng/30 μL can be approximately converted into 0.6 pM), and therefore the sensitivity of CARP is high enough to accomplish evaluation of CRISPR/Cas9 site-specific function and validation of sgRNA sequences. A standard curve of Ct value and a linear function with R 2 > 0.99 were thus generated, indicating excellent stability of the detection method (Fig. 3f ). If the melting curve was used to help analyze the experimental results, the sensitivity of CARP was higher; the minimum sensitivity was 0.002 ng (Fig. 3d ). In addition, the whole detection process of CARP could be completed within 2 h. The results fully demonstrated that CARP is able to distinguish single-base mutation. Therefore, CARP has the potential to quickly evaluate site-specific and validate sgRNA sequences. After determining the ability of CARP to detect single-base mutations, we investigated the site-specific of Cas9-sgRNA. First, site mutation of DNA was detected with qPCR-based CARP, and all target DNAs containing site mutations were treated with the same Cas9-sgRNA complexes. In the process of gene editing with the CRISPR/Cas system, a protospacer adjacent motif (PAM) sequence was found to be indispensable. Only DNA sequences containing PAM can be recognized and cleaved by Cas9-sgRNA. Because target DNA contained three bases of the PAM and its adjacent two bases (Fig. 4a) . The results show that the adjacent base mutation of PAM has little effect on cleavage efficiency of Cas9-sgRNA. Once the mutation site is in the PAM sequence, the cleavage efficiency of Cas9-sgRNA is significantly reduced. This proves that the complete PAM sequence is critical to excellent cleavage activity and affinity of the CRISPR/Cas9 system. Spontaneous or chemical-inducible single-base mutations are common phenomena in genomic regions. Thus it is important to determine whether the cleavage activity of the CRISPR/Cas9 system can be affected by these mutations. We chose five important mutations for investigation: methylated cytosine [36, 37] , two U-G base mutations (the mutation can be caused by intracellular folate deficiency) [38] , and two I-C base mutations (caused by nitrite stimulation of cells) [39] . The results demonstrated that these mutations have little influence on the cleavage efficiency of the CRISPR/Cas9 system (Fig. 4b) . Therefore, the effects of these mutations can be excluded in the design of Cas9-sgRNA. The effects of single-base mismatch of sgRNA on cleavage efficiency The specificity of the CRISPR/Cas9 system is largely determined by sgRNA which contains recognition sequences of target DNA. sgRNAs are then designed to verify whether off-target DNA can be successfully distinguished by CARP. In order to ensure the reliability of the study, we respectively designed two groups of sgRNAs based on two different sites, with each group containing five single-base mismatch sgRNAs and mismatch sites adjacent to PAM. The results revealed that the cleavage activity of the CRISPR/sgRNA system was greatly affected by a single-base mismatch of sgRNA. A small fraction of target DNA was still cut (Fig. 5) , indicating that nonspecific cleavage could occur at undesigned sites. Therefore, CARP can be used to effectively avoid off-target effects before CRISPR/Cas9-related genome engineering events by exploring unpredictable DNA cleavage sites with different off-target substrate designs. (Fig. 6) . The optimal and ineffective sgRNAs for cleavage efficiency can be successfully differentiated by comparing Ct values. The results suggest that the CARP method has excellent ability to screen sgRNAs with the highest affinity for CRISPR/Cas9-related gene editing. The CRISPR/Cas system is generally one of three major types: types I, II and III [40] . Compared with types I and III, the type II system only contains a single guided RNA (sgRNA) and one protease, such as Cas9, which is the hallmark protease of a type II system guided by sgRNA to complete a double-stranded DNA (dsDNA) cleavage event. Because of the simplicity and more mature commercial applications, we chose Cas9 as the cleavage enzyme for the CARP assay. But that does not mean that other CRISPR enzymes cannot be employed in this way. Various studies have reported the use of other types of CRISPR enzymes for nucleic acid detection. In 2017, Gootenberg et al. combined Cas13a with isothermal amplification to develop a specific high-sensitivity enzymatic reporter unlocking (SHERLOCK) technique, which provides rapid nucleic acid detection with ultra-high sensitivity and specificity [12] . Subsequently, a number of CRISPR-based nucleic acid detection technologies were established, for example, a multiplexed and portable nucleic acid detection platform with Cas13, Cas12a, and Csm6 [6] , a CRISPR/Cas12a-based platform for rapid and visual nucleic acid detection [41] , and an ultra-sensitive and high-throughput CRISPR-powered COVID-19 diagnostic method [42] . Both these common enzymes and other CRISPR enzymes that have the ability to cleave DNA can be used to build an inclusive CARP detection system. If off-target activity can be avoided, the CRISPR/Cas9 system has the ability to distinguish single-base mutations [25] . To overcome off-target activity, several approaches have been explored. A web selection tool for sgRNA, CHOPCHOP, was developed by Labun et al. to avoid off-target activity of the CRISPR/Cas9 system [27] . Double nicking (combining a Cas9 nickase with paired guide RNAs) used to enhance genome editing specificity can reduce off-target activity by 50to 1500-fold [43] . In this study, all sgRNA were designed with CHOPCHOP to reduce off-target activity. Double nicking was the first step of the CARP system to ensure the specificity of the detection, so that target DNA could be specifically cleaved by a pair of Cas9-sgRNA complexes. Then the cleavage products were linked by T4 DNA ligase. Finally, PCR reactions were carried out, and a pair of reverse primers was cleverly designed for amplifying target DNA. If cleavage or ligation reaction was unsuccessful, the reverse primers would elongate along the DNA template back to back, and the amplification would fail. If the target DNA was cut and the cleavage products were linked to each other into longer DNAs, the reverse primers would be changed into normal primers, and thus the target DNA was successfully amplified. On the basis of these features, the CARP system has high specificity, which was fully verified by detecting various single-base mutations. PCR is the most common amplification method. As an exponential amplification technique, its amplification efficiency for DNA is excellent. The amplification efficiency of the CARP depends on PCR, which is the third step of the overall detection procedure, and the minimum detectable quantity in this work was 0.02 ng. Therefore, the sensitivity of our method is high enough to detect low-concentration substrate DNA. High sensitivity is a key criterion in an ideal in vitro detection method; the other key criterion is time efficiency. Compared with CRISPR/Cas9-related isothermal amplification and next-generation sequencing [44, 45] , our method is a time-saving process for prescreening gene sites and sgRNA sequences. The whole CARP detection procedure can be completed within 2 h. On the basis of the advantages above, CARP will be a promising method for evaluating sitespecific functions of CRISPR/Cas9 and validating sgRNA sequences. In summary, we have developed a specific, sensitive, and time-saving Cas9/sgRNA-assisted reverse PCR method to evaluate site-specific function of CRISPR/Cas9 and validate sgRNA sequence. This method was fully verified by detecting variety site mutations and evaluating the affinity of sgRNA sequences. The detection method based on PCR can be used directly in traditional reaction conditions of PCR, without the need to establish an additional detection system. Therefore, it can be popularized easily in laboratories that have PCR instrument. According to the advantages above, it will provide a new effective analysis tool to help us explore CRISPR/Cas9 genome engineering related research. CRISPR-a widespread system that provides acquired resistance against phages in bacteria and archaea Identification of genes that are associated with DNA repeats in prokaryotes Cas9 as a versatile tool for engineering biology CRISPR-based technologies for the manipulation of eukaryotic genomes Development and applications of CRISPR-Cas9 for genome engineering Multiplexed and portable nucleic acid detection platform with Cas13, Cas12a, and Csm6 Advance genome editing technologies in the treatment of human diseases: CRISPR therapy Multiplex genome engineering using CRISPR/Cas systems RNA-guided human genome engineering via Cas9 Genetic screens in human cells using the CRISPR-Cas9 system Genome-scale CRISPR-Cas9 knockout screening in human cells Nucleic acid detection with CRISPR-Cas13a/C2c2 Field-deployable viral diagnostics using CRISPR-Cas13 Detection of unamplified target genes via CRISPR-Cas9 immobilized on a graphene field-effect transistor SHERLOCK: nucleic acid detection with CRISPR nucleases CRISPR-Cas12-based detection of SARS-CoV-2 CRISPR-LbCpf1 prevents choroidal neovascularization in a mouse model of age-related macular degeneration Development of a gene-editing approach to restore vision loss in Leber congenital amaurosis type 10 CRISPR-Cas9 corrects Duchenne muscular dystrophy exon 44 deletion mutations in mice and human cells CRISPR/ Cas9-mediated miR-29b editing as a treatment of different types of muscle atrophy in mice Correction of muscular dystrophies by CRISPR gene editing CRISPR-engineered T cells in patients with refractory cancer Modeling breast cancer using CRISPR-Cas9-mediated engineering of human breast organoids Safety and feasibility of CRISPR-edited T cells in patients with refractory nonsmall-cell lung cancer DNA targeting specificity of RNA-guided Cas9 nucleases Enhanced guide-RNA design and targeting analysis for precise CRISPR genome editing of single and consortia of industrially relevant and non-model organisms CHOPCHOP v3: expanding the CRISPR web toolbox beyond genome editing CRISPRscan: designing highly efficient sgRNAs for CRISPR-Cas9 targeting in vivo CRISPR-DO for genome-wide CRISPR design and optimization GuideScan software for improved single and paired CRISPR guide RNA design Highly efficient targeted mutagenesis of Drosophila with the CRISPR/Cas9 system In vitro evaluation of CRISPR/Cas9 function by an electrochemiluminescent assay Label-free CRISPR/Cas9 assay for site-specific nucleic acid detection Enzymatic amplification of beta-globin genomic sequences and restriction site analysis for diagnosis of sickle cell anemia Evaluation of CRISPR/Cas9 site-specific function and validation of sgRNA sequence by a Cas9/sgRNA-assisted Detection of target DNA with a novel Cas9/sgRNAs-associated reverse PCR (CARP) technique The DNA (cytosine-5) methyltransferases Bone marrow cells from vitamin B12-and folate-deficient patients misincorporate uracil into DNA Nitrite-induced deamination and hypochlorite-induced oxidation of DNA in intact human respiratory tract epithelial cells Evolution and classification of the CRISPR-Cas systems Cas12aVDet: a CRISPR/Cas12a-based platform for rapid and visual nucleic acid detection Ultrasensitive and high-throughput CRISPR-powered COVID-19 diagnosis Double nicking by RNA-guided CRISPR Cas9 for enhanced genome editing specificity Clustered regularly interspaced short palindromic repeats/Cas9 triggered isothermal amplification for site-specific nucleic acid detection A high-throughput screening strategy for detecting CRISPR-Cas9 induced mutations using next-generation sequencing Publisher's note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations Acknowledgments This work was supported by grants from the Key R & D and Promotion Program of Henan Province (192102310075) and the Joint construction program of Henan medical science and technology research plan (LHGJ20200067). Conflict of interest The authors declare that they have no conflict of interest. The online version contains supplementary material available at https://doi.org/10.1007/s00216-021-03173-2.