Abstract
Recent years have seen the development of computational tools to assist researchers in performing CRISPR-Cas9 experiment optimally. More specifically, these tools aim to maximize on-target activity (guide efficiency) while also minimizing potential off-target effects (guide specificity) by analyzing the features of the target site. Nonetheless, currently available tools cannot robustly predict experimental success as prediction accuracy depends on the approximations of the underlying model and how closely the experimental setup matches the data the model was trained on. Here, we present an overview of the available computational tools, their current limitations and future considerations. We discuss new trends around personalized health by taking genomic variants into account when predicting target sites as well as discussing other governing factors that can improve prediction accuracy.
Introduction
The CRISPR-Cas9 system allows for targeted editing of DNA in vitro. The system is targeted to the DNA via association with a guide RNA (gRNA) molecule, which binds to the targeted DNA through base complementarity and enables precise DNA cleavage (). This cleavage is then repaired via various pathways, which can be exploited for different outcomes (). Knockouts can be achieved through error prone repair via the Non-homologous End Joining pathway, which can introduce mutations and disrupt gene function. Targeted integration of a sequence (called a knock-in) can be achieved via the Homology Directed Repair pathway, which uses a provided DNA template to repair the cleavage. Activation or repression of a gene can be achieved by targeting catalytically inert Cas9 fused to a transcription activator or repressor to the promoter (). All of these approaches require the accurate and efficient targeting of the CRISPR-Cas9 system to the desired location. The success of an experiment using the CRISPR-Cas9 system therefore hinges on the correct identification of the optimal target-site and subsequent design of the complimentary gRNA (; ). While databases of validated gRNAs exist for various genomes [e.g., Cas-Database (Park et al., 2016) for knockout applications and () for gene activation/repression], these libraries are generic and may not be well-suited for specific research purposes. The design of custom gRNAs is hence frequently required.
A successful gRNA must maximize on-target activity (guide efficiency) while also minimizing potential off-target effects (guide specificity). Balancing these two requirements can be a combinatorial challenging task and as a result, significant effort in the recent years has been focused on developing computational tools to assist in the design of gRNAs. These tools are designed to assist researchers in the selection of best target sites by helping them exclude undesirable targets based on predicted low efficiency or a high potential for off-target effects. Here, we present an overview of the development of tools for the design of CRISPR-Cas9 gRNAs, their current limitations and future considerations.
Predicting On-Target Activity
Initially, CRISPR-Cas9 was thought to be able to target any 20 base-pair sequence that was flanked by a protospacer adjacent motif (PAM). Different Cas and related enzymes target different PAMs, and there is ongoing researching into designing enzymes with specific PAM recognition ability (). However, the most commonly used SpCas9, and the focus of this review, targets an NGG motif. As such, early tools for target site selection were simple pattern recognition programs that identified instances of this motif (Upadhyay and Sharma, 2014; Xie et al., 2014; Zhu et al., 2014). In some cases, information about where in a gene the target site fell (e.g., within an intron or exon) was also incorporated, allowing researchers to draw some conclusions on the likelihood of a functional effect. However, subsequent studies showed that CRISPR-Cas9 displayed a wide variety of activities across different target sites, leading to the conclusion that some target sites are inherently more effective (, ; ; , ; ; Yang et al., 2013; ; ; Shalem et al., 2014; Wang et al., 2014; ; ).
This discovery led to a series of large-scale screens of CRISPR-Cas9 activity across a variety of target sites and organisms, aimed at identifying what features contributed to targeting efficiency (; , ; ; ; ). These studies helped identify some key rules for optimizing gRNA design. This include avoiding poly-T sequences, limiting the GC content and a G immediately upstream of the PAM (i.e., an GNGG motif) (Ren et al., 2014; Shalem et al., 2014; Wong et al., 2015). Building on this research, computational methods were created for predicting on-target activity. The initial studies focused on the contribution of the target site sequence, by measuring the activity of 1000s of target sites. These studies differed in how they defined the target sites, with some considering only the 20 bp target sequence () while others included the PAM and flanking sequence (; ; Wong et al., 2015). They also differed in how they represented the target site to the mathematical model, i.e., the feature space. The studies used different combinations of position specific nucleotides and dinucleotides, global nucleotide counts, GC content, etc. More recent studies have also begun to include non-sequence information, such as thermodynamic stability of the gRNA and position of the cut site relative to the transcription start site (TSS) (; Wong et al., 2015; ).
The differences in experimental design means that each study produced a unique predictive model, with different rules for CRISPR-Cas9 activity. Supplementary Table 1 presents a selection of tools that demonstrate the variety of data types, features, and model implementations used. Despite the differences in the model, however, certain key features were repeatedly found to be important. These include position-specific nucleotides, such as a G preceding the PAM being a strong indicator of CRISPR-Cas9 activity, or global variables such as GC content and gRNA melting temperature were consistently reported as being important (Wong et al., 2015). Comparing the distribution of important features along the target site, the majority are found within the ∼10–12 bp adjacent to the PAM, a region that has become known as the seed region (). This region is typically thought to be critical for CRISPR-Cas9 activity, as this region binds the DNA first following recognition of the PAM (; Shibata et al., 2017).
The models also differed in what machine learning technique was used in their construction. While predicting activity using linear regression showed some success (), the more successful models used more complex approaches such as Random Forest (Wilson et al., 2018) and Support Vector Machines (; Wong et al., 2015; ), which consider interactions between the individual features (). The success of these more complex models suggests that there is no single feature that governs activity, but rather a combination of interactions.
Despite the extensive training of the models, the accuracy of their predictions varies widely. A recent review of different on-target efficiency models found that no model was consistently accurate across a number of independent dataset, recording high accuracy only when tested on the original training dataset (). This discrepancy is likely due to the differences in how the various studies conducted their experiments. Consistent with this, a recent review found that predictive models performed best when the CRISPR-Cas9 expression system matched the one used in the training dataset (). This would suggest that experimental conditions do affect the final model.
This same study also found that the method used to transcribe the gRNAs may also influence activity prediction. Typically, gRNAs are transcribed in cells from a U6 promoter or in vitro from a T7 promoter (Zhang et al., 2017). These promoters have differing transcription requirements, such as different polymerases and a G (for the U6 promoter) or GG (for the T7 promoter) at the TSS. These differences appear to influence any predictive model, with models performing better when applied to datasets that use the same promoter as what was used in the model’s training set (). Currently, no one predictive model is able to account for gRNA transcription method. Some pipelines (such as CRISPOR), use multiple predictive models allowing researchers to select the most appropriate score ().
It is also highly likely that the manner in which CRISPR-Cas9 activity is measured impacts the final model. There is currently no consensus in the literature in how CRISPR-Cas9 activity should be measured. Some studies measure activity by the rate at which mutations are introduced through sequencing the target site (; ), while others measure activity by the size of the phenotypic change (such as drug resistance, cell viability, or protein expression) (, ; ). While measuring activity via sequencing may prove a more direct measurement, it is also a costlier approach and does not provide information about whether the induced mutations are functional. Conversely, while phenotypic screens are easier to perform at scale they rely on the CRISPR-Cas9 introducing functional mutations, which may in turn lead to increase in false-negatives (i.e., mutations that do not cause a functional effect). These differences in experimental design likely translate to differences in the model.
reported that two of the most important variables for predicting CRISPR-Cas9 activity are the position of the target site relative to the TSS and position within the protein. However, this study was performed using a dataset that reported CRISPR-Cas9 activity based on a combination of changes in drug resistance and expression of cell-surface proteins. Given mutations near the TSS of a gene are more likely to induce a functional change, it is highly likely that the importance of target position is inflated in a phenotypic screen. In fact, a recent study comparing the impact of different training sets found that training a predictive model using sequencing-based measurements of CRISPR-Cas9 activity yields more generalizable predictions (Wilson et al., 2018). Phenotypic-trained models are governed by features such as position of the target site relative to TSS and do not generalize to other datasets.
Because the training dataset has such a strong influence on the final predictive model, it is therefore critical to know on what a model was trained before use. As a rule of thumb, phenotypic-trained models will be better suited to identifying target sites that induce functional changes but are limited to experiments with the same condition as the training set. In contrast, sequencing-based models are more universally applicable, but are only capable of predicting genotype changes not their functional result.
Predicting Off-Target Activity
Identifying potential off-target sites is typically done by repurposing computational tools used for high-throughput sequencing read alignment. Here, the target site is treated as a read and realigned back to the reference genome in order to identify similar locations that may be inadvertently targeted by the CRISPR-Cas9. Alignment of the short target sequences is typically achieved using tools such as Bowtie and BWA, which are better suited for handling short sequences compared to other traditional tools such as BLAST.
These repurposed tools, however, are not the optimal solution for this problem. Searching for potential off-target sites requires the identification of small sequence motifs (20 bp + PAM) with often many mismatches. Traditional alignment tools are not equipped to identify such small, divergent sequences. Typically, Bowtie alignments allow only up to three mismatches while BWA allows up to 5, resulting in more divergent off-targets being missed. In fact comparison of off-target identification pipelines with experimentally validated CRISPR-Cas9 off-targets shows that these traditional alignment methods not only miss the high-mismatch off-target but even some with only one mismatch (Tsai et al., 2015; ), suggesting these tools are generally poorly suited for this problem.
Implementation of new alignment methods, such as bi-directional alignments (), will be required to accurately identify all potential off-targets. Typically, aligners work by first matching a small portion of the query sequence (known as the seed) and then extending the seed out in a direction and testing the match. Bi-directional aligners work by extending the initial seed region in both directions. Using these more powerful alignment tools will be important for correctly identifying all potential off-targets.
Further complicating the matter is that not every putative off-target is actually functional (i.e., off-targets that are actually cleaved by CRISPR-Cas9). As such, naive alignment methods hence return a large number of false-positives, potentially leading to the erroneous disqualification of the optimal target site.
A recent study comparing experimentally validated off-targets and those predicted by alignment tools, showed that the prediction tools over-estimate the number of potential off-targets by up to 10-fold (). In order to reduce the number of false-positive predictions, off-target predictors often limit potential off-targets to a maximum number of mismatches and only very specific PAMs (). However, experimental studies, have shown that off-targets can differ significantly from the original target site, meaning this approach often results in false-negatives (Tsai et al., 2015, 2017; ). Predictive programs therefore need to balance the false-positives and false-negatives. To compensate for this, several studies have developed scoring algorithms, which attempt to predict the activity of a potential off-target so that false-positives can be filtered out.
The two most popular scoring methods are the MIT-Broad score () and the CFD score (). Both of these scoring algorithms are based on “synthetic” datasets, whereby a series of gRNAs targeting a specific dataset were mutated such that every one, two, and three base mismatch combination was represented. The ability of the gRNAs to cleave the target site were then measured, and the results used to construct a Linear Regression algorithm to score the off-target sites. Despite the theory behind both methods being similar, they differ in how the final model is constructed. While the MIT-Broad algorithm considers only the 20 bp target sequence (i.e., does not included the PAM), the CFD score takes the PAM sequence into account, scoring target sites as less active if they possess non-canonical PAMs. A recent comparison of the method tested their ability to accurately predict the off-target activity of different experimental datasets and concluded that the CFD score performed the best (). These methods however, are limited in the features they consider, focusing predominantly on the number and position of mismatches. Two recently developed off-target methods Elevation () and CRISTA () expand the feature set, including features such as gRNA secondary structure, genomic location and overlap with other features of interest such as DNase 1 Hypersensitive sites. These models are also capable of distinguishing between mismatches that occur through wobble pairing, and those caused by DNA/RNA bulges which may have structural implications. Inclusion of these additional features allows the models to better predict off-target activity and they outperform the CFD and MIT-Broad methods on independent datasets (; ). Supplementary Table 2 catalogs some of the more common off-target detection tools and summarizes their key differences.
Future Perspective
A key goal of future research will be to improve the accuracy of predictive models by incorporating additional features. Current methods for predicting target efficiency and specificity are based solely on the sequence of the target site. However, it is now accepted that chromatin environment (; ; ; ; ) can influence CRISPR-Cas9 activity. Early studies mapping the genome wide binding of inert Cas9 enzymes using ChIP-seq showed a preference for DNAse sensitive regions (; Wu et al., 2014; O’Geen et al., 2015), which are typically more accessible environments. This was supported by later studies which showed that high-activity target sites were often enriched for histone modifications associated with open-chromatin environments ().
A direct link between chromatin and CRISPR-Cas9 activity was shown in 2016, where a pair of studies demonstrated that the presence of nucleosomes at the target site physically blocked CRISRP-Cas9’s access and reduced overall activity (; ). The differences in chromatin environment likely explain why the same CRISPR-Cas9 target site can display different activities across cell-lines (). There is also evidence that off-target activity is influenced by chromatin accessibility, with the CROP-IT pipeline including this information into their off-target model (Singh et al., 2015). Incorporating environmental information in future predictive models will help improve accuracy and will be critical if the technology is to be applied in the clinic. Such modeling may also allow for the selective targeting of individual tissues by leveraging the differences in chromatin environments.
Incorporation of chromatin environments would likely also improve off-target predictions, which is thought to be more susceptible to chromatin accessibility. Besides chromatin information, future off-target pipelines should also focus on including variant information. A recent study demonstrated that the variance between individuals has a dramatic effect on the off-target landscape, with point mutations creating and destroying potential off-target sites (). Such information is critical for the application of CRISPR technology in almost all fields, as not taking an individual’s unique genome into account could have deleterious side-effects (, ; Scott and Zhang, 2017).
Future models may also not only be able to predict the success of CRISPR-Cas9 editing, but also the outcome. By targeting sites with microhomology and exploiting the microhomology-mediated repair pathway, researchers may be able to delete specific DNA segments and thereby control the outcome of CRISPR-Cas9 editing (; Yao et al., 2017). Additionally, a recent study found that the mutations induced by repair of CRISPR-Cas9 cleavage were non-random and determined by the target sequence (van Overbeek et al., 2016). Such a finding suggests that it would be possible to predict the mutational outcome of CRISPR-Cas9 editing, allowing for researchers to make precise edits without the need of using knock-ins.
The optimal future pipeline will incorporate all of these factors into both on- and off-target activity predictions (Figure 1). Such a pipeline could also provide a method by which experimentally validated predictions could be reintegrated into the training data for the models, to continue to improve accuracy. Future models may also predict success of other CRISPR-Cas9 applications such as knock-ins (), which involve the repair of the double strand break using a supplied template, and base-editing, where a Cas9 fusion protein converts one base into another without cleavage ().
FIGURE 1
Conclusion
Computational tools for the prediction of CRISPR-Cas9 activity are necessary for the efficient design of experiments. However, current tools are hampered by a range of issues, such as disparate training data sources, which results in models not generalizing, as well as limitations in our current understand of factors that drive CRISPR-Cas9 activity. As our understanding improves, we will be able to incorporate new features into predictive models to increase their accuracy. This will be vital for applying CRISPR-Cas9 in clinical applications, where an individual’s genomic variations may alter activity patters of CRISPR-Cas9. Until then, it is important that the data used to train a predictive model is understood before it is used, to ensure models are only applied in appropriate circumstances.
Statements
Author contributions
LW researched the software highlighted in the paper. LW and DB designed and wrote the review paper. AO contributed to Section “Future Perspective.”
Funding
The paper was supported by Commonwealth Scientific and Industrial Research Organisation and CSIRO Synthetic Biology Future Science Platform.
Conflict of interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Supplementary material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fphar.2018.00749/full#supplementary-material
References
1
AbadiS.YanW. X.AmarD.MayroseI. (2017). A machine learning approach for predicting CRISPR-Cas9 cleavage efficiencies and patterns underlying its mechanism of action.PLoS Comput. Biol.13:e1005807. 10.1371/journal.pcbi.1005807
2
BaeS.KweonJ.KimH. S.KimJ. S. (2014a). Microhomology-based choice of Cas9 nuclease target sites.Nat. Methods11705–706. 10.1038/nmeth.3015
3
BaeS.ParkJ.KimJ. S. (2014b). Cas-OFFinder: a fast and versatile algorithm that searches for potential off-target sites of Cas9 RNA-guided endonucleases.Bioinformatics301473–1475. 10.1093/bioinformatics/btu048
4
CameronP.FullerC. K.DonohoueP. D.JonesB. N.ThompsonM. S.CarterM. M.et al (2017). Mapping the genomic landscape of CRISPR-Cas9 cleavage.Nat. Methods14600–606. 10.1038/nmeth.4284
5
CanverM. C.JoungJ. K.PinelloL. (2018). Impact of genetic variation on CRISPR-Cas targeting.CRISPR J.1159–170. 10.1089/crispr.2017.0016
6
CanverM. C.LessardS.PinelloL.WuY.IlboudoY.SternE. N.et al (2017). Variant-aware saturating mutagenesis using multiple Cas9 nucleases identifies regulatory elements at trait-associated loci.Nat. Genet.49625–634. 10.1038/ng.3793
7
CanzarS.SalzbergS. L. (2017). Short read mapping: an algorithmic tour.Proc. IEEE Inst. Electr. Electron. Eng.105436–458. 10.1109/JPROC.2015.2455551
8
Cebrian-SerranoA.DaviesB. (2017). CRISPR-Cas orthologues and variants: optimizing the repertoire, specificity and delivery of genome engineering tools.Mamm. Genome28247–261. 10.1007/s00335-017-9697-9694
9
ChariR.MaliP.MoosburnerM.ChurchG. M. (2015). Unraveling CRISPR-Cas9 genome engineering parameters via a library-on-library approach.Nat. Methods12823–826. 10.1038/nmeth.3473
10
ChenF.DingX.FengY.SeebeckT.JiangY.DavisG. D. (2017). Targeted activation of diverse CRISPR-Cas systems for mammalian genome editing via proximal CRISPR targeting.Nat. Commun.8:14958. 10.1038/ncomms14958
11
CongL.RanF. A.CoxD.LinS.BarrettoR.HabibN.et al (2013). Multiplex genome engineering using CRISPR/Cas systems.Science339819–823. 10.1126/science.1231143
12
DoenchJ. G.FusiN.SullenderM.HegdeM.VaimbergE. W.DonovanK. F.et al (2016). Optimized sgRNA design to maximize activity and minimize off-target effects of CRISPR-Cas9.Nat. Biotechnol.34184–191. 10.1038/nbt.3437
13
DoenchJ. G.HartenianE.GrahamD. B.TothovaZ.HegdeM.SmithI.et al (2014). Rational design of highly active sgRNAs for CRISPR-Cas9-mediated gene inactivation.Nat. Biotechnol.321262–1267. 10.1038/nbt.3026
14
FarasatI.SalisH. M. (2016). A biophysical model of crispr/cas9 activity for rational design of genome editing and gene regulation.PLoS Comput. Biol.12:e1004724. 10.1371/journal.pcbi.1004724
15
FuY.FodenJ. A.KhayterC.MaederM. L.ReyonD.JoungJ. K.et al (2013). High-frequency off-target mutagenesis induced by CRISPR-Cas nucleases in human cells.Nat. Biotechnol.31822–826. 10.1038/nbt.2623
16
FuY.SanderJ. D.ReyonD.CascioV. M.JoungJ. K. (2014). Improving CRISPR-Cas nuclease specificity using truncated guide RNAs.Nat. Biotechnol.32279–284. 10.1038/nbt.2808
17
GaudelliN. M.KomorA. C.ReesH. A.PackerM. S.BadranA. H.BrysonD. I.et al (2017). Programmable base editing of A●T to G●C in genomic DNA without DNA cleavage.Nature551464–471. 10.1038/nature24644
18
HaeusslerM.SchönigK.EckertH.EschstruthA.MiannéJ.RenaudJ. B.et al (2016). Evaluation of off-target and on-target scoring algorithms and integration into the guide RNA selection tool CRISPOR.Genome Biol.17:148. 10.1186/s13059-016-1012-1012
19
HorlbeckM. A.GilbertL. A.VillaltaJ. E.AdamsonB.PakR. A.ChenY.et al (2016a). Compact and highly active next-generation libraries for CRISPR-mediated gene repression and activation.Elife5:e19760. 10.7554/eLife.19760
20
HorlbeckM. A.WitkowskyL. B.GuglielmiB.ReplogleJ. M.GilbertL. A.VillaltaJ. E.et al (2016b). Nucleosomes impede Cas9 access to DNA in vivo and in vitro.Elife5:e12677. 10.7554/eLife.12677
21
HsuP. D.ScottD. A.WeinsteinJ. A.RanF. A.KonermannS.AgarwalaV.et al (2013). DNA targeting specificity of RNA-guided Cas9 nucleases.Nat. Biotechnol.31827–832. 10.1038/nbt.2647
22
IsaacR. S.JiangF.DoudnaJ. A.LimW. A.NarlikarG. J.AlmeidaR. (2016). Nucleosome breathing and remodeling constrain CRISPR-Cas9 function.Elife5:e13450. 10.7554/eLife.13450
23
JinekM.ChylinskiK.FonfaraI.HauerM.DoudnaJ. A.CharpentierE. (2012). A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity.Science337816–821. 10.1126/science.1225829
24
JinekM.EastA.ChengA.LinS.MaE.DoudnaJ. (2013). RNA-programmed genome editing in human cells.Elife2:e00471. 10.7554/eLife.00471
25
KimH.KimJ. S. (2014). A guide to genome engineering with programmable nucleases.Nat. Rev. Genet.15321–334. 10.1038/nrg3686
26
KnightS. C.XieL.DengW.GuglielmiB.WitkowskyL. B.BosanacL.et al (2015). Dynamics of CRISPR-Cas9 genome interrogation in living cells.Science350823–826. 10.1126/science.aac6572
27
Koike-YusaH.LiY.TanE. P.Velasco-HerreraM. D. C.YusaK. (2014). Genome-wide recessive genetic screening in mammalian cells with a lentiviral CRISPR-guide RNA library.Nat. Biotechnol.32267–273. 10.1038/nbt.2800
28
KuscuC.ArslanS.SinghR.ThorpeJ.AdliM. (2014). Genome-wide analysis reveals characteristics of off-target sites bound by the Cas9 endonuclease.Nat. Biotechnol.32677–683. 10.1038/nbt.2916
29
La RussaM. F.QiL. S. (2015). The new state of the art: cas9 for gene activation and repression.Mol. Cell. Biol.353800–3809. 10.1128/MCB.00512-515
30
LessardS.FrancioliL.AlfoldiJ.TardifJ. C.EllinorP. T.MacArthurD. G.et al (2017). Human genetic variation alters CRISPR-Cas9 on- and off-targeting specificity at therapeutically implicated loci.Proc. Natl. Acad. Sci. U.S.A.114E11257–E11266. 10.1073/pnas.1714640114
31
ListgartenJ.WeinsteinM.KleinstiverB. P.SousaA. A.JoungJ. K.CrawfordJ.et al (2018). Prediction of off-target activities for the end-to-end design of CRISPR guide RNAs.Nat. Biomed. Eng.238–47. 10.1038/s41551-017-0178-6
32
LiuX.HommaA.SayadiJ.YangS.OhashiJ.TakumiT. (2016). Sequence features associated with the cleavage efficiency of CRISPR/Cas9 system.Sci. Rep.6:19675. 10.1038/srep19675
33
MaliP.YangL.EsveltK. M.AachJ.GuellM.DiCarloJ. E.et al (2013). RNA-guided human genome engineering via Cas9.Science339823–826. 10.1126/science.1232033
34
McKinneyB. A.ReifD. M.RitchieM. D.MooreJ. H. (2006). Machine learning for detecting gene-gene interactions: a review.Appl. Bioinform.577–88.
35
MerkleF. T.NeuhausserW. M.SantosD.ValenE.GagnonJ. A.MaasK.et al (2015). Efficient CRISPR-Cas9-mediated generation of knockin human pluripotent stem cells lacking undesired mutations at the targeted locus.Cell Rep.11875–883. 10.1016/j.celrep.2015.04.007
36
Moreno-MateosM. A.VejnarC. E.BeaudoinJ. D.FernandezJ. P.MisE. K.KhokhaM. K.et al (2015). CRISPRscan: designing highly efficient sgRNAs for CRISPR-cas9 targeting in vivo.Nat. Methods12982–988. 10.1038/nmeth.3543
37
O’GeenH.HenryI. M.BhaktaM. S.MecklerJ. F.SegalD. J. (2015). A genome-wide analysis of Cas9 binding specificity using ChIP-seq and targeted sequence capture.Nucleic Acids Res.433389–3404. 10.1093/nar/gkv137
38
ParkJ.KimJ. S.BaeS. (2016). Cas-Database: web-based genome-wide guide RNA library design for gene knockout screens using CRISPR-Cas9.Bioinformatics322017–2023. 10.1093/bioinformatics/btw103
39
RenX.YangZ.XuJ.SunJ.MaoD.HuY.et al (2014). Enhanced specificity and efficiency of the CRISPR/Cas9 system with optimized sgRNA parameters in Drosophila.Cell Rep.91151–1162. 10.1016/j.celrep.2014.09.044
40
ScottD. A.ZhangF. (2017). Implications of human genetic variation in CRISPR-based therapeutic genome editing.Nat. Med.231095–1101. 10.1038/nm.4377
41
ShalemO.SanjanaN. E.HartenianE.ShiX.ScottD. A.MikkelsonT.et al (2014). Genome-scale CRISPR-Cas9 knockout screening in human cells.Science34384–87. 10.1126/science.1247005
42
ShibataM.NishimasuH.KoderaN.HiranoS.AndoT.UchihashiT.et al (2017). Real-space and real-time dynamics of CRISPR-Cas9 visualized by high-speed atomic force microscopy.Nat. Commun.8:1430. 10.1038/s41467-017-01466-1468
43
SinghR.KuscuC.QuinlanA.QiY.AdliM. (2015). Cas9-chromatin binding information enables more accurate CRISPR off-target prediction.Nucleic Acids Res.43:e118. 10.1093/nar/gkv575
44
TsaiS. Q.NguyenN. T.Malagon-LopezJ.TopkarV. V.AryeeM. J.JoungJ. K. (2017). Circle-seq: a highly sensitive in vitro screen for genome-wide CRISPR-Cas9 nuclease off-targets.Nat. Methods14607–614. 10.1038/nmeth.4278
45
TsaiS. Q.ZhengZ.NguyenN. T.LiebersM.TopkarV. V.ThaparV.et al (2015). Guide-seq enables genome-wide profiling of off-target cleavage by CRISPR-Cas nucleases.Nat. Biotechnol.33187–197. 10.1038/nbt.3117
46
UpadhyayS. K.SharmaS. (2014). SSFinder: high throughput CRISPR-cas target sites prediction tool.Biomed. Res. Int.2014:742482. 10.1155/2014/742482
47
van OverbeekM.CapursoD.CarterM. M.ThompsonM. S.FriasE.RussC.et al (2016). DNA repair profiling reveals nonrandom outcomes at cas9-mediated breaks.Mol. Cell63633–646. 10.1016/j.molcel.2016.06.037
48
WangT.WeiJ. J.SabatiniD. M.LanderE. S. (2014). Genetic screens in human cells using the CRISPR-Cas9 system.Science34380–84. 10.1126/science.1246981
49
WilsonL. O. W.RetiD.O’BrienA. R.DunneR. A.BauerD. C. (2018). High activity target-site identification using phenotypic independent CRISPR-Cas9 core functionality.CRISPR J.1182–190. 10.1089/crispr.2017.0021
50
WongN.LiuW.WangX. (2015). WU-CRISPR: characteristics of functional guide RNAs for the CRISPR/Cas9 system.Genome Biol.16:218. 10.1186/s13059-015-0784-780
51
WuX.ScottD. A.KrizA. J.ChiuA. C.HsuP. D.DadonD. B.et al (2014). Genome-wide binding of the CRISPR endonuclease Cas9 in mammalian cells.Nat. Biotechnol.32670–676. 10.1038/nbt.2889
52
XieS.ShenB.ZhangC.HuangX.ZhangY. (2014). sgRNAcas9: a software package for designing CRISPR sgRNA and evaluating potential off-target cleavage sites.PLoS One9:e100448. 10.1371/journal.pone.0100448
53
YangH.WangH.ShivalilaC. S.ChengA. W.ShiL.JaenischR. (2013). One-step generation of mice carrying reporter and conditional alleles by CRISPR/Cas-mediated genome engineering.Cell1541370–1379. 10.1016/j.cell.2013.08.022
54
YaoX.WangX.HuX.LiuZ.LiuJ.ZhouH.et al (2017). Homology-mediated end joining-based targeted integration using CRISPR/Cas9.Cell Res.27801–814. 10.1038/cr.2017.76
55
ZhangT.GaoY.WangR.ZhaoY. (2017). Production of guide RNAs in vitro and in vivo for CRISPR using ribozymes and RNA polymerase II promoters.Bio. Protoc.7:e2148. 10.21769/BioProtoc.2148
56
ZhuL. J.HolmesB. R.AroninN.BrodskyM. H. (2014). CRISPRseek: a bioconductor package to identify target-specific guide RNAs for CRISPR-cas9 genome-editing systems.PLoS One9:e108424. 10.1371/journal.pone.0108424
Summary
Keywords
CRISPR-Cas9, bioinformatics, off-target finder, activity prediction, chromatin, machine learning
Citation
Wilson LOW, O’Brien AR and Bauer DC (2018) The Current State and Future of CRISPR-Cas9 gRNA Design Tools. Front. Pharmacol. 9:749. doi: 10.3389/fphar.2018.00749
Received
19 January 2018
Accepted
20 June 2018
Published
12 July 2018
Volume
9 - 2018
Edited by
Danilo Maddalo, Novartis (Switzerland), Switzerland
Reviewed by
Leyuan Ma, Koch Institute for Integrative Cancer Research at MIT, United States; Amir Taheri-Ghahfarokhi, AstraZeneca (Sweden), Sweden
Updates
Copyright
© 2018 Wilson, O’Brien and Bauer.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Denis C. Bauer, Denis.Bauer@CSIRO.au
This article was submitted to Cancer Molecular Targets and Therapeutics, a section of the journal Frontiers in Pharmacology
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.