- Open Access
A large-scale whole-genome sequencing analysis reveals highly specific genome editing by both Cas9 and Cpf1 (Cas12a) nucleases in rice
- Xu Tang†1,
- Guanqing Liu†2, 3,
- Jianping Zhou†1,
- Qiurong Ren1,
- Qi You2, 3,
- Li Tian1,
- Xuhui Xin1,
- Zhaohui Zhong1,
- Binglin Liu1,
- Xuelian Zheng1,
- Dengwei Zhang1,
- Aimee Malzahn4,
- Zhiyun Gong2,
- Yiping Qi4, 5Email author,
- Tao Zhang2, 3Email author and
- Yong Zhang1Email authorView ORCID ID profile
© The Author(s). 2018
- Received: 1 May 2018
- Accepted: 27 May 2018
- Published: 4 July 2018
Targeting specificity has been a barrier to applying genome editing systems in functional genomics, precise medicine and plant breeding. In plants, only limited studies have used whole-genome sequencing (WGS) to test off-target effects of Cas9. The cause of numerous discovered mutations is still controversial. Furthermore, WGS-based off-target analysis of Cpf1 (Cas12a) has not been reported in any higher organism to date.
We conduct a WGS analysis of 34 plants edited by Cas9 and 15 plants edited by Cpf1 in T0 and T1 generations along with 20 diverse control plants in rice. The sequencing depths range from 45× to 105× with read mapping rates above 96%. Our results clearly show that most mutations in edited plants are created by the tissue culture process, which causes approximately 102 to 148 single nucleotide variations (SNVs) and approximately 32 to 83 insertions/deletions (indels) per plant. Among 12 Cas9 single guide RNAs (sgRNAs) and three Cpf1 CRISPR RNAs (crRNAs) assessed by WGS, only one Cas9 sgRNA resulted in off-target mutations in T0 lines at sites predicted by computer programs. Moreover, we cannot find evidence for bona fide off-target mutations due to continued expression of Cas9 or Cpf1 with guide RNAs in T1 generation.
Our comprehensive and rigorous analysis of WGS data across multiple sample types suggests both Cas9 and Cpf1 nucleases are very specific in generating targeted DNA modifications and off-targeting can be avoided by designing guide RNAs with high specificity.
Bacterial type II CRISPR-Cas9 systems can effectively induce RNA-guided DNA double strand breaks (DSBs) , making them popular tools for genome editing in bacteria , animal cells , mammalian systems [4–7], and plants [8–11]. The most widely used Streptococcus pyogenes Cas9 (SpCas9) uses ~ 20 nucleotides (nt) of a single guide RNA (sgRNA) to recognize a complementary target DNA site along with an NGG protospacer adjacent motif (PAM) [1, 12]. More recently, type V CRISPR-Cpf1 (CRISPR-Cas12a) was shown to mediate efficient genome editing in human cells  and plants [14–16]. Cpf1(Cas12a) uses ~ 23 nt of an RNA guide to target DNA with a TTTV PAM . RNA-guided nucleases (RGNs) such as Cas9 and Cpf1 represent versatile genome editing tools that promise to advance basic science, enable personalized medicine, and accelerate crop breeding. However, Cas9 may cause undesired off-target mutations due to sgRNAs recognizing DNA sequences with one to a few nucleotide mismatches, albeit with reduced nuclease binding and cleavage activity [1, 6, 17, 18]. Although similar rules apply to Cpf1, recent studies in human cells [19, 20] have shown Cpf1 is generally more specific than Cas9.
Understanding the scope of off-target mutations in Cas9- or Cpf1-edited crops is critical for research and regulation. Previously, whole-genome sequencing (WGS) was applied for detecting off-target mutations by Cas9 in Arabidopsis , rice , and tomato . Unfortunately, these studies either only looked at potential off-target sites predicted by computer programs or fell short of full analysis of all the mutations identified by WGS in edited plants. Without inclusion of enough necessary controls, such WGS studies had limited power for isolating off-target mutations in edited plants because they were unable to fully assess the levels of preexisting mutations, spontaneous mutations, and mutations caused by tissue culture- and Agrobacterium-mediated transformation. Genome-wide identification of off-target mutations by Cas9 or Cpf1 will be empowered only if all background mutations can be isolated. Furthermore, WGS-based off-target analysis of Cpf1 has not been reported in any higher organism. In recent years, WGS studies on Cas9-edited mice have generated contrasting results; one study found few off-target mutations  while another found many . This controversy raised the urgency for comprehensive and rigorous analyses of off-target mutations using WGS in edited animals and plants. We reasoned a large-scale and well-designed study is required for comprehensive assessment of off-target effects in crops by Cas9 and Cpf1, two leading CRISPR genome editing systems. Here, we describe a large-scale WGS study to assess off-target effects of Cas9 and Cpf1 in rice, an important food crop. Our results suggest off-target mutations of Cas9 and Cpf1 are largely negligible when compared to spontaneous mutations or mutations caused by tissue culture and Agrobacterium infection in edited plants. The resulting knowledge is likely to serve as an important reference for plant researchers and regulatory agencies.
Detection of off-target, spontaneous, and background mutations
To assess mutations generated by tissue culture and Agrobacterium infection, we produced and sequenced four types of control plants: tissue culture only, tissue culture with Agrobacterium, tissue culture with Agrobacterium transformation of Cas9 without sgRNA, and tissue culture with Agrobacterium transformation of Cpf1 without crRNA (Fig. 1b and Additional file 1: Figure S4). Tissue culture is known to be mutagenic, causing somaclonal variations . Indeed, the two tissue culture-only samples contained an average of 114 SNVs and 36 indels (Fig. 2a, b), resulting in a background mutation rate of 1.86 × 10− 7, which is similar to the rates (1.7–3.3 × 10− 7) previously published . Importantly, similar numbers of SNVs were observed from Agrobacterium-infected or Cas9/Cpf1 backbone-transformed plants (Fig. 2a). These three controls generated ~ 15 to 41 more indels compared to tissue culture-only samples (Fig. 2b), suggesting Agrobacterium infection is mutagenic with a preference for introducing indels. This warrants further investigation as these three controls show large variations on indel counts. We mapped all identified mutations from these four control types to the rice genome across 12 chromosomes (Additional file 1: Figure S5). Further analysis of the genome-wide distribution of these background mutations revealed high enrichment of SNVs in transposable elements (TEs) and repeats (Fig. 2c), as well as high enrichment of indels in repeats (Fig. 2d).
SNVs and indels identified in edited T0 plants are largely background mutations
Identification of true off-target mutations in T0 plants
Cas9-E sgRNA was predicted by CRISPOR and Cas-OFFinder to contain six off-target sites when up to a 3-nt mismatch was allowed (Fig. 4a and Additional file 2: Table S3). However, no off-target mutations were found at these predicted sites. Although the two Cas9-E T0 lines shared seven SNVs and three indels (Fig. 4b), these ten shared mutations had very poor sequence homology to the target site (Fig. 4e). Only five sites contained the NGG PAM. Among them, the site sharing highest sequence homology with the target site still contained a 10-nt mismatch, making it unlikely to be a true off-target site. Unlike indels found in Cas9-J samples, these putative off-target mutations are mostly SNVs (Fig. 4b, f). Furthermore, both independent T0 lines always carried the same mutant alleles (Fig. 4f and Additional file 3: Figure S11). These observations suggest that the ten shared mutations of two Cas9-E T0 lines were not caused by Cas9, but were pre-existing mutations from a parental line.
Cas9 was previously shown to induce off-target mutations at sites with missing or extra nucleotides when compared to the target site, which forms bulges when targeted by guide RNAs . To detect such off-target mutations, we extracted all T0 mutation site flanking sequences (25 bp upstream and downstream) and aligned them to corresponding sgRNA/crRNA sequences using BLAST. Only Cas9-J1 and Cas9-J2 samples had alignments to the Cas9-J-sgRNA01 target (15 in Cas9-J1 and ten in Cas9-J2); other samples had no hit. None of the detected mutations were caused by bulge-forming DNA–sgRNA recognition. We also investigated whether DNA translocation events were induced by Cas9 or Cpf1 by searching for structural variants (SVs) and gene fusion events in the whole rice genome. We did not detect any translocation event in all T0 lines. Given the level of nuclease-induced DNA translocation can be used for assessing targeting specificity , absence of detectable translocation events in all T0 samples here indicates these Cas9 and Cpf1 reagents are indeed very specific, limiting cleavage activity almost exclusively to the target sites.
No evidence of off-target mutations in T1 plants
These new mutations were mapped to the rice genome together with new mutations that were discovered in WT plants across two generations (Additional file 3: Figure S17). The genome distribution of these new mutations in T1 lines also showed enrichment in repeats (Additional file 3: Figure S16), consistent with the spontaneous mutations discovered in the WT (Fig. 2c, d). Detailed analysis of SNVs among all sample types revealed T1 lines have higher rates of G:C > A:T transitions than T0 lines (Additional file 3: Figure S18), consistent with the observation on spontaneous mutations in Arabidopsis . Further analysis of T1 lines either with or without the Cas9 transgene did not reveal any difference in the numbers of new SNVs and indels among these two subpopulations (Fig. 5c, d). By applying similar methods from the analysis of T0 plants, we were unable to identify any off-target mutations by Cas9 or Cpf1 in T1 lines. Given most T1 lines analyzed still carry the RGN constructs, our results suggest continued expression of Cas9 or Cpf1 constructs did not cause de novo off-target mutations in T1 lines.
To further assess the new mutations found in T1 lines, we calculated and compared the allele frequency of SNVs and indels among four groups: tissue culture controls, T0 plants, T1 plants, and WT (Fig. 5e). The tissue culture controls and Cas9/Cpf1 T0 lines share strikingly similar (mostly heterozygous-like) allele frequency distribution. This reiterates our earlier conclusion that all mutations in T0 samples (except a few found in Cas9-J samples) are background mutations. By contrast, T1 plants show more homozygous-like SNVs (0.75 to 1.0 in allele frequency) and somatic-like indels (0 to 0.25 in allele frequency). This trend of rapidly fixing SNVs and the increase of somatic indels in T1 is interesting and relatively in line with the observation in WT plants.
More attention has been given to the specificity of CRISPR-Cas RGN systems in humans than in animals or plants due to medicinal applications of RGNs. Earlier WGS studies in human cells found low incidence of off-target mutations by Cas9 [37, 38]. Recently, two WGS off-target studies in mice showed conflicting results [24, 25]. However, the study that claimed unexpected large-scale off-target effects by Cas9 may be flawed due to limitations in its experimental design and WGS data analysis . Given the wide adoption of CRISPR-Cas systems in agriculture, with genome-edited crop products reaching market in record time , it becomes urgent to conduct large-scale and exhaustive WGS analysis of off-target effects by Cas9 and Cpf1 (Cas12a), two leading RGN systems, in agriculturally important crops. Such studies will help assess the safety of Cas9 and Cpf1 in precise crop breeding as well as provide valuable information to scientists, breeders, regulators, and consumers.
In this study, we conducted a large-scale WGS analysis for detecting potential off-target mutations caused by 12 Cas9 sgRNAs and three Cpf1 crRNAs in rice, an important food crop. We confirmed WGS-identified mutations by Sanger sequencing at randomly selected sites with a 100% success rate (Additional file 2: Table S4), which is consistent with the high quality of our WGS data. Our experimental design took into account background mutations caused by tissue culture and Agrobacterium-mediated transformation, pre-existing mutations in parents, and spontaneous mutations that arise from seed propagation. Through sequencing 20 control plants of different types and 49 Cas9 or Cpf1-edited T0 and T1 plants, we only found true off-target mutations in two T0 lines expressing Cas9 protein with Cas9-J-sgRNA01. Importantly, these empirically validated off-target sites can be readily predicted computationally. Our examination of T1 plants that continue to carry Cas9-sgRNA or Cpf1-crRNA did not reveal off-target mutations, suggesting continued presence of the RGN reagents with varying activity in plants does not cause off-target mutations if the guide RNAs are well-designed for specificity. This observation is also highly significant because it encourages the use of Cas9 and Cpf1 in certain breeding applications that may require expression of RGNs across several generations. For example, a RGN cassette may be introduced from a transgenic line into a transformation-recalcitrant variety of the same plant species for genome editing with simple genetic crossing.
Our study also provided insights into avoiding off-target effects of Cas9 and Cpf1 in edited crops. To minimize off-target effects, many systems have been developed, including paired Cas9 nickases , high fidelity Cas9 proteins [42–44], FokI-dCas9 fusions [45, 46], truncated sgRNAs , and ribonucleotide protein (RNP) delivery . To assess and identify off-target sites, in vivo [18, 36] and in vitro [49–51] tools have also been developed in human cells, which may be applied in plants. Our WGS analysis with WT SpCas9 and LbCpf1 proteins did not find off-target mutations for 14 out of 15 guide RNAs tested in T0 and T1 plants, suggesting utilization of a high-fidelity enzyme, which are typically of lower activity, may be unnecessary in crop applications. When a mismatch up to 3 nt of the protospacer is allowed, Cas9-OFFinder programs predicted a total of 37 off-target sites for 7 out of 11 Cas9 sgRNAs. Yet, we could not detect any mutations at these putative off-target sites. Alternatively, Cas9-OFFinder predicted all the off-target sites that we identified for Cas9-J-sgRNA01; many of the sites have just 1-nt mismatch to the protospacer of the target site. Therefore, we can deduce a simple rule to alleviate off-target effects: make sure even the highest scored potential off-target sites will have at least a 2-nt mismatch to the seed sequence of the protospacer. We note this may not always be possible if the target sequence shares many homologous sequences in the genome. For example, maize has a very repetitive genome and wheat has A, B, D sub-genomes that share high similarity. In these cases, targeted amplicon sequencing using next-generation sequencing technologies may be an appropriate and cost-effective method to look for off-target mutations .
Finally, we hope our data can be a valuable reference for regulatory agencies and other entities. It is reasonable and necessary to scrutinize any new technology for its efficacy and safety. Cas9 and Cpf1, as new crop-breeding technologies, are no exception. Although Cas9-based off-target effects have been studied by WGS in plants [21–23], our study differs from previous studies significantly with regard to scale, depth, and comprehensiveness. Our research also represents the first report of using WGS to assess off-targeting by Cpf1 in any edited higher eukaryotic organism. We could not find any off-target mutations in 47 out of 49 rice plants edited by 11 Cas9-sgRNA and three Cpf1-crRNA constructs. This precise level of genome modification is in stunning contrast to many conventional breeding technologies. For example, we found that even the safest breeding approach, harvesting seeds from parental lines, introduces ~ 30 to 50 spontaneous mutations into the next generation in rice. We also observed ~ 200 tissue culture-introduced somaclonal variations per rice plant, even though few are affecting coding sequences. In conclusion, our data support a recent call to “Regulate genome-edited products, not genome editing itself” .
Plant material and growth conditions
This study used the rice variety Nipponbare (Oryza sativa L. ssp. Japonica cv. Nipponbare). All plants were grown in growth chambers under controlled environmental conditions with a 16/8 h light/dark regime at 28 °C and 60% relative humidity.
Plasmids encoding Cas9 and a single sgRNA were generated by ligating annealed oligos with a 4-bp overhang into a BsaI digested backbone (either pZHY988 or pTX172) [54–56]. Plasmids with two sgRNAs were created by ligating pZHY988 with a 485-bp fragment, after digestion with BsaI. This 485-bp fragment contains two sgRNAs generated by overlap extension PCR . All CRISPR-Cpf1 (CRISPR-Cas12a) nuclease expression vectors were reported in our previous study . The sequences of all primers used to construct vectors are shown in Additional file 2: Table S5.
Rice stable transformation
Agrobacterium-mediated rice transformation was performed as described in published protocols  with slight modification [16, 56]. The binary vectors were introduced into Agrobacterium tumefaciens strain EHA105 by the freeze-thaw method . For rice transformation, dehusked seeds were sterilized with 70% ethanol for 1 min. Afterwards, seeds were washed five times with sterile water, then further sterilized for 15 min with a 2.5% sodium hypochlorite solution containing a drop of Tween 20. The washing and sterilization step were repeated, this time without addition of Tween. Seeds were then rinsed an additional five times before being dried on sterilized filter paper and cultured on solid medium at 28 °C in a dark growth chamber for 2–3 weeks. Actively growing calli were collected for subculture at 28 °C in the dark for 1–2 weeks. Agrobacterium cultures were collected and resuspended in liquid medium (OD600 = 0.06–0.1) containing 100 μM acetosyringone. Rice calli were immersed in the Agrobacterium suspension for 30 min, then dried on sterilized filter paper and co-cultured for 3 days on solid medium at 25 °C in a dark growth chamber. The infected calli were moved to a sterile plastic bottle and washed five times with sterile water to remove excessive Agrobacterium. After being dried on sterilized filter paper, these calli were transferred onto screening medium at 28 °C in a dark growth chamber for 5 weeks. During the screening stage, infected calli were transferred to fresh screening medium every 2 weeks. After the screening stage, actively growing calli were moved onto regenerative medium for regeneration at 28 °C with a 16 h light/8 h dark cycle. After 3–4 weeks, transgenic seedlings were transferred to sterile plastic containers containing fresh solid medium and grown for 2–3 weeks before being transferred into soil. Transgenic rice plants were grown in a growth chamber at 28 °C with a 16 h light/8 h dark cycle.
Mutagenesis analysis at target sites
Genomic DNA was extracted from transgenic plants using the CTAB method . The genomic region flanking the CRISPR target site for each gene was amplified and sequenced. Samples with heterozygous and biallelic mutations were decoded using CRISP-ID .
WGS and data analysis
For each sample, about 1 g of fresh leaves were collected from seedlings between 5 and 6 weeks old. DNA samples were extracted using the Plant Genome DNA Kit (Tiangen) as described by the manufacturer. All 69 samples were sequenced by Bionova (Beijing, China) using the Illumina X10 platform. Adapters were trimmed using SKEWER (v. 0.2.2)  and the Illumina TruSeq adapter. Cleaned reads were mapped to rice reference sequence TIGR7 (http://rice.plantbiology.msu.edu/)  with BWA (v. 0.7.15) software . The Genome Analysis Toolkit (GATK)  was used to realign reads near indels and recalibrate base quality scores by following GATK best practices . A known SNPs and indels database for GATK best practices was downloaded from Rice SNP-Seek Database (http://snp-seek.irri.org/) . Whole genome SNVs were detected by LoFreq , MuTect2 , and VarScan2 . Whole genome indels were identified using MuTect2 , VarScan2, and Pindel . Bedtools  and BCFtools  were used to process overlapping SNVs/indels. Off-target sites were predicted with CRISPOR  online and Cas-OFFinder software  by allowing up to 10-nt mismatch. A genome-wide map of mutations was plotted with Circos software . Structural variants and translocation events were analyzed using TopHat2  with the ‘—fusion-search’ parameter, DELLY  with default parameters, and manually checking with IGV software . NCBI BLAST+ with the ‘-task blastn-short’ parameter for off-target mutation site analysis, which includes mismatch, deletion, and insertion. Read-mapping screenshots were from Golden Helix GenomeBrowse ® visualization tool v2.1. Data processing and analyses were completed using R and Python. One T1 sample of Cas9-CC2 was excluded from analyses due to contamination of fungal DNA.
We thank Dr. Yusheng Qin at Soil and Fertilizer Research Institute in Sichuan Academy of Agricultural Sciences for plant care.
This work was supported by the National Natural Science Foundation of China (31330017 and 31771486) and the Sichuan Youth Science and Technology Foundation (2017JQ0005) to YZ. This work was also supported by grants, including the National Transgenic Major Project of China (2018ZX08022001–003) to YZ, and TZ, the Jiangsu Specially-Appointed Professor and the Priority Academic Program Development of Jiangsu Higher Education Institutions (PAPD) to TZ, and startup funds provided by University of Maryland to YQ.
Availability of data and materials
The WGS raw data reported in this article are listed in Additional file 2: Table S2 and have been deposited in the Sequence Read Archive in National Center for Biotechnology Information (NCBI) under the accession number PRJNA420933 (https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA420933)  and the Genome Sequence Archive in Beijing Institute of Genomics (BIG) under the accession number PRJCA000656 (http://bigd.big.ac.cn/bioproject/browse/PRJCA000656) .
YZ proposed the project. YQ, TZ, and YZ designed the experiments. XT, JZ, DZ, and AM generated CRISPR-Cas nuclease vectors. XT, JZ, LT, XX, DZ, and XZ generated stable transgenic rice and identified the rice mutants. XT, JZ, QR, ZZ, BL, and XZ prepared samples for WGS and confirmed WGS results with Sanger sequencing. GL, QY, ZG, and TZ performed WGS data analysis. YQ, TZ, and YZ analyzed the data and wrote the manuscript. All authors participated in discussion and revision of the manuscript. All authors read and approved the final manuscript.
The 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.Google Scholar
- Jiang W, Bikard D, Cox D, Zhang F, Marraffini LA. RNA-guided editing of bacterial genomes using CRISPR-Cas systems. Nat Biotechnol. 2013;31:233–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Hwang WY, Fu Y, Reyon D, Maeder ML, Tsai SQ, Sander JD, et al. Efficient genome editing in zebrafish using a CRISPR-Cas system. Nat Biotechnol. 2013;31:227–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Jinek M, East A, Cheng A, Lin S, Ma E, Doudna J. RNA-programmed genome editing in human cells. elife. 2013;2:e00471.View ArticlePubMedPubMed CentralGoogle Scholar
- Mali P, Yang L, Esvelt KM, Aach J, Guell M, DiCarlo JE, et al. RNA-guided human genome engineering via Cas9. Science. 2013;339:823–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Cong L, Ran FA, Cox D, Lin S, Barretto R, Habib N, et al. Multiplex genome engineering using CRISPR/Cas systems. Science. 2013;339:819–23.View ArticlePubMedPubMed CentralGoogle Scholar
- Cho SW, Kim S, Kim JM, Kim JS. Targeted genome engineering in human cells with the Cas9 RNA-guided endonuclease. Nat Biotechnol. 2013;31:230–2.View ArticlePubMedGoogle Scholar
- Li JF, Norville JE, Aach J, McCormack M, Zhang D, Bush J, et al. Multiplex and homologous recombination-mediated genome editing in Arabidopsis and Nicotiana benthamiana using guide RNA and Cas9. Nat Biotechnol. 2013;31:688–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Shan Q, Wang Y, Li J, Zhang Y, Chen K, Liang Z, et al. Targeted genome modification of crop plants using a CRISPR-Cas system. Nat Biotechnol. 2013;31:686–8.View ArticlePubMedGoogle Scholar
- Nekrasov V, Staskawicz B, Weigel D, Jones JD, Kamoun S. Targeted mutagenesis in the model plant Nicotiana benthamiana using Cas9 RNA-guided endonuclease. Nat Biotechnol. 2013;31:691–3.Google Scholar
- Lowder LG, Zhang D, Baltes NJ, Paul JW, Tang X, Zheng X, et al. A CRISPR/Cas9 toolbox for multiplexed plant genome editing and transcriptional regulation. Plant Physiol. 2015;169:971–85.View ArticlePubMedPubMed CentralGoogle Scholar
- Mojica FJ, Diez-Villasenor C, Garcia-Martinez J, Almendros C. Short motif sequences determine the targets of the prokaryotic CRISPR defence system. Microbiology. 2009;155:733–40.View ArticlePubMedGoogle Scholar
- Zetsche B, Gootenberg JS, Abudayyeh OO, Slaymaker IM, Makarova KS, Essletzbichler P, et al. Cpf1 is a single RNA-guided endonuclease of a class 2 CRISPR-Cas system. Cell. 2015;163:759–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Tang X, Lowder LG, Zhang T, Malzahn AA, Zheng X, Voytas DF, et al. A CRISPR-Cpf1 system for efficient genome editing and transcriptional repression in plants. Nat Plants. 2017;3:17018.View ArticlePubMedGoogle Scholar
- Endo A, Masafumi M, Kaya H, Toki S. Efficient targeted mutagenesis of rice and tobacco genomes using Cpf1 from Francisella novicida. Sci Rep. 2016;6:38169.Google Scholar
- Zhong Z, Zhang Y, You Q, Tang X, Ren Q, Liu S, et al. Plant genome editing using FnCpf1 and LbCpf1 nucleases at redefined and altered PAM sites. Mol Plant. 2018. https://doi.org/10.1016/j.molp.2018.03.008.
- Fu Y, Foden JA, Khayter C, Maeder ML, Reyon D, Joung JK, et al. High-frequency off-target mutagenesis induced by CRISPR-Cas nucleases in human cells. Nat Biotechnol. 2013;31:822–6.View ArticlePubMedPubMed CentralGoogle Scholar
- Tsai SQ, Zheng Z, Nguyen NT, Liebers M, Topkar VV, Thapar V, et al. Guide-seq enables genome-wide profiling of off-target cleavage by CRISPR-Cas nucleases. Nat Biotechnol. 2015;33:187–97.View ArticlePubMedGoogle Scholar
- Kim D, Kim J, Hur JK, Been KW, Yoon SH, Kim JS. Genome-wide analysis reveals specificities of Cpf1 endonucleases in human cells. Nat Biotechnol. 2016;34:863–8.View ArticlePubMedGoogle Scholar
- Kleinstiver BP, Tsai SQ, Prew MS, Nguyen NT, Welch MM, Lopez JM, et al. Genome-wide specificities of CRISPR-Cas Cpf1 nucleases in human cells. Nat Biotechnol. 2016;34:869–74.Google Scholar
- Feng Z, Mao Y, Xu N, Zhang B, Wei P, Yang DL, et al. Multigeneration analysis reveals the inheritance, specificity, and patterns of CRISPR/Cas-induced gene modifications in Arabidopsis. Proc Natl Acad Sci U S A. 2014;111:4632–7.View ArticlePubMedPubMed CentralGoogle Scholar
- Zhang H, Zhang J, Wei P, Zhang B, Gou F, Feng Z, et al. The CRISPR/Cas9 system produces specific and homozygous targeted gene editing in rice in one generation. Plant Biotechnol J. 2014;12:797–807.View ArticlePubMedGoogle Scholar
- Nekrasov V, Wang C, Win J, Lanz C, Weigel D, Kamoun S. Rapid generation of a transgene-free powdery mildew resistant tomato by genome deletion. Sci Rep. 2017;7:482.View ArticlePubMedPubMed CentralGoogle Scholar
- Iyer V, Shen B, Zhang W, Hodgkins A, Keane T, Huang X, et al. Off-target mutations are rare in Cas9-modified mice. Nat Methods. 2015;12:479.View ArticlePubMedGoogle Scholar
- Schaefer KA, Wu WH, Colgan DF, Tsang SH, Bassuk AG, Mahajan VB. Unexpected mutations after CRISPR-Cas9 editing in vivo. Nat Methods. 2017;14:547–8.View ArticlePubMedPubMed CentralGoogle Scholar
- Yang N, Xu XW, Wang RR, Peng WL, Cai L, Song JM, et al. Contributions of Zea mays subspecies mexicana haplotypes to modern maize. Nat Commun. 2017;8:1874.View ArticlePubMedPubMed CentralGoogle Scholar
- Ossowski S, Schneeberger K, Lucas-Lledo JI, Warthmann N, Clark RM, Shaw RG, et al. The rate and molecular spectrum of spontaneous mutations in Arabidopsis thaliana. Science. 2010;327:92–4.View ArticlePubMedGoogle Scholar
- Yang S, Wang L, Huang J, Zhang X, Yuan Y, Chen JQ, et al. Parent-progeny sequencing indicates higher mutation rates in heterozygotes. Nature. 2015;523:463–7.View ArticlePubMedGoogle Scholar
- Evans DA. Somaclonal variation--genetic basis and breeding applications. Trends Genet. 1989;5:46–50.View ArticlePubMedGoogle Scholar
- Wei FJ, Kuang LY, Oung HM, Cheng SY, Wu HP, Huang LT, et al. Somaclonal variation does not preclude the use of rice transformants for genetic screening. Plant J. 2016;85:648–59.View ArticlePubMedGoogle Scholar
- Haeussler M, Schonig K, Eckert H, Eschstruth A, Mianne J, Renaud JB, et al. Evaluation of off-target and on-target scoring algorithms and integration into the guide RNA selection tool CRISPOR. Genome Biol. 2016;17:148.View ArticlePubMedPubMed CentralGoogle Scholar
- Bae S, Park J, Kim JS. Cas-offinder: a fast and versatile algorithm that searches for potential off-target sites of Cas9 RNA-guided endonucleases. Bioinformatics. 2014;30:1473–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Hsu PD, Scott DA, Weinstein JA, Ran FA, Konermann S, Agarwala V, et al. DNA targeting specificity of RNA-guided Cas9 nucleases. Nat Biotechnol. 2013;31:827–32.Google Scholar
- Meng X, Hu X, Liu Q, Song X, Gao C, Li J, et al. Robust genome editing of CRISPR-Cas9 at NAG PAMs in rice. Sci China Life Sci. 2018;61:122–25.Google Scholar
- Lin Y, Cradick TJ, Brown MT, Deshmukh H, Ranjan P, Sarode N, et al. CRISPR/Cas9 systems have off-target activity with insertions or deletions between target DNA and guide RNA sequences. Nucleic Acids Res. 2014;42:7473–85.Google Scholar
- Frock RL, Hu J, Meyers RM, Ho YJ, Kii E, Alt FW. Genome-wide detection of DNA double-stranded breaks induced by engineered nucleases. Nat Biotechnol. 2015;33:179–86.View ArticlePubMedGoogle Scholar
- Veres A, Gosis BS, Ding Q, Collins R, Ragavendran A, Brand H, et al. Low incidence of off-target mutations in individual CRISPR-Cas9 and TALEN targeted human stem cell clones detected by whole-genome sequencing. Cell Stem Cell. 2014;15:27–30.Google Scholar
- Smith C, Gore A, Yan W, Abalde-Atristain L, Li Z, He C, et al. Whole-genome sequencing analysis reveals high specificity of CRISPR/Cas9 and TALEN-based genome editing in human iPSCs. Cell Stem Cell. 2014;15:12–3.Google Scholar
- Editorial. CRISPR off-targets: a reassessment. Nat Methods. 2018;15:229–30.Google Scholar
- Waltz E. CRISPR-edited crops free to enter market, skip regulation. Nat Biotechnol. 2016;34:582.View ArticlePubMedGoogle Scholar
- Ran FA, Hsu PD, Lin CY, Gootenberg JS, Konermann S, Trevino AE, et al. Double nicking by RNA-guided CRISPR Cas9 for enhanced genome editing specificity. Cell. 2013;154:1380–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Kleinstiver BP, Pattanayak V, Prew MS, Tsai SQ, Nguyen NT, Zheng Z, et al. High-fidelity CRISPR-Cas9 nucleases with no detectable genome-wide off-target effects. Nature. 2016;529:490–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Slaymaker IM, Gao L, Zetsche B, Scott DA, Yan WX, Zhang F. Rationally engineered Cas9 nucleases with improved specificity. Science. 2016;351:84–8.View ArticlePubMedGoogle Scholar
- Chen JS, Dagdas YS, Kleinstiver BP, Welch MM, Sousa AA, Harrington LB, et al. Enhanced proofreading governs CRISPR-Cas9 targeting accuracy. Nature. 2017;550:407–10.View ArticlePubMedPubMed CentralGoogle Scholar
- Tsai SQ, Wyvekens N, Khayter C, Foden JA, Thapar V, Reyon D, et al. Dimeric CRISPR RNA-guided FokI nucleases for highly specific genome editing. Nat Biotechnol. 2014;32:569–76.Google Scholar
- Guilinger JP, Thompson DB, Liu DR. Fusion of catalytically inactive Cas9 to FokI nuclease improves the specificity of genome modification. Nat Biotechnol. 2014;32:577–82.Google Scholar
- Fu Y, Sander JD, Reyon D, Cascio VM, Joung JK. Improving CRISPR-Cas nuclease specificity using truncated guide RNAs. Nat Biotechnol. 2014;32:279–84.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim S, Kim D, Cho SW, Kim J, Kim JS. Highly efficient RNA-guided genome editing in human cells via delivery of purified Cas9 ribonucleoproteins. Genome Res. 2014;24:1012–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Cameron P, Fuller CK, Donohoue PD, Jones BN, Thompson MS, Carter MM, et al. Mapping the genomic landscape of CRISPR-Cas9 cleavage. Nat Methods. 2017;14:600–6.View ArticlePubMedGoogle Scholar
- Kim D, Bae S, Park J, Kim E, Kim S, Yu HR, et al. Digenome-seq: genome-wide profiling of CRISPR-Cas9 off-target effects in human cells. Nat Methods. 2015;12:237–43.View ArticlePubMedGoogle Scholar
- Tsai SQ, Nguyen NT, Malagon-Lopez J, Topkar VV, Aryee MJ, Joung JK. Circle-seq: a highly sensitive in vitro screen for genome-wide CRISPR-Cas9 nuclease off-targets. Nat Methods. 2017;14:607–14.View ArticlePubMedPubMed CentralGoogle Scholar
- You Q, Zhong Z, Ren Q, Hassan F, Zhang Y, Zhang T. CRISPRMatch: an automatic calculation and visualization tool for high-throughput CRISPR genome-editing data analysis. Int J Biol Sci. 14:858–62.Google Scholar
- Carroll D, Van Eenennaam AL, Taylor JF, Seger J, Voytas DF. Regulate genome-edited products, not genome editing itself. Nat Biotechnol. 2016;34:477–9.View ArticlePubMedGoogle Scholar
- Tang X, Zheng X, Qi Y, Zhang D, Cheng Y, Tang A, et al. A single transcript CRISPR-Cas9 system for efficient genome editing in plants. Mol Plant. 2016;9:1088–91.View ArticlePubMedGoogle Scholar
- Zhou J, Deng K, Cheng Y, Zhong Z, Tian L, Tang X, et al. CRISPR-Cas9 based genome editing reveals new insights into microRNA function and regulation in rice. Front Plant Sci. 2017;8:1598.View ArticlePubMedPubMed CentralGoogle Scholar
- Zheng X, Yang S, Zhang D, Zhong Z, Tang X, Deng K, et al. Effective screen of CRISPR/Cas9-induced mutants in rice by single-strand conformation polymorphism. Plant Cell Rep. 2016;35:1545–54.View ArticlePubMedGoogle Scholar
- Hiei Y, Ohta S, Komari T, Kumashiro T. Efficient transformation of rice (Oryza sativa l.) mediated by Agrobacterium and sequence analysis of the boundaries of the T-DNA. Plant J. 1994;6:271–82.View ArticlePubMedGoogle Scholar
- Weigel D, Glazebrook J. Transformation of Agrobacterium using the freeze-thaw method. CSH Protoc. 2006. https://doi.org/10.1101/pdb.prot4666.
- Stewart CN Jr, Via LE. A rapid ctab DNA isolation technique useful for rapd fingerprinting and other PCR applications. BioTechniques. 1993;14:748–50.PubMedGoogle Scholar
- Dehairs J, Talebi A, Cherifi Y, Swinnen JV. CRISP-ID: decoding CRISPR mediated indels by Sanger sequencing. Sci Rep. 2016;6:28973.Google Scholar
- Jiang H, Lei R, Ding SW, Zhu S. Skewer: a fast and accurate adapter trimmer for next-generation sequencing paired-end reads. BMC Bioinformatics. 2014;15:182.View ArticlePubMedPubMed CentralGoogle Scholar
- Kawahara Y, de la Bastide M, Hamilton JP, Kanamori H, McCombie WR, Ouyang S, et al. Improvement of the Oryza sativa Nipponbare reference genome using next generation sequence and optical map data. Rice. 2013;6:4.Google Scholar
- Li H. Aligning sequence reads, clone sequences and assembly contigs with BWA-MEM. In: ArXiv e-prints. 2013;arXiv:1303.3997v2.Google Scholar
- McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, et al. The genome analysis toolkit: a mapreduce framework for analyzing next-generation DNA sequencing data. Genome Res. 2010;20:1297–303.View ArticlePubMedPubMed CentralGoogle Scholar
- Van der Auwera GA, Carneiro MO, Hartl C, Poplin R, Del Angel G, Levy-Moonshine A, et al. From FastQ data to high confidence variant calls: the genome analysis toolkit best practices pipeline. Curr Protoc Bioinformatics. 2013;43:11.10.1–33.Google Scholar
- Mansueto L, Fuentes RR, Borja FN, Detras J, Abriol-Santos JM, Chebotarov D, et al. Rice SNP-seek database update: new SNPs, indels, and queries. Nucleic Acids Res. 2017;45:D1075–D81.Google Scholar
- Wilm A, Aw PP, Bertrand D, Yeo GH, Ong SH, Wong CH, et al. LoFreq: a sequence-quality aware, ultra-sensitive variant caller for uncovering cell-population heterogeneity from high-throughput sequencing datasets. Nucleic Acids Res. 2012;40:11189–201.Google Scholar
- Cibulskis K, Lawrence MS, Carter SL, Sivachenko A, Jaffe D, Sougnez C, et al. Sensitive detection of somatic point mutations in impure and heterogeneous cancer samples. Nat Biotechnol. 2013;31:213–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Koboldt DC, Zhang Q, Larson DE, Shen D, McLellan MD, Lin L, et al. VarScan 2: somatic mutation and copy number alteration discovery in cancer by exome sequencing. Genome Res. 2012;22:568–76.Google Scholar
- Ye K, Schulz MH, Long Q, Apweiler R, Ning Z. Pindel: a pattern growth approach to detect break points of large deletions and medium sized insertions from paired-end short reads. Bioinformatics. 2009;25:2865–71.View ArticlePubMedPubMed CentralGoogle Scholar
- Quinlan AR. BEDTools: the swiss-army tool for genome feature analysis. Curr Protoc Bioinformatics. 2014;47:1121–34.Google Scholar
- Li H. A statistical framework for SNP calling, mutation discovery, association mapping and population genetical parameter estimation from sequencing data. Bioinformatics. 2011;27:2987–93.View ArticlePubMedPubMed CentralGoogle Scholar
- Krzywinski M, Schein J, Birol I, Connors J, Gascoyne R, Horsman D, et al. Circos: an information aesthetic for comparative genomics. Genome Res. 2009;19:1639–45.View ArticlePubMedPubMed CentralGoogle Scholar
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14:R36.Google Scholar
- Rausch T, Zichner T, Schlattl A, Stutz AM, Benes V, Korbel JO. DELLY: structural variant discovery by integrated paired-end and split-read analysis. Bioinformatics. 2012;28:i333–i9.Google Scholar
- Thorvaldsdottir H, Robinson JT, Mesirov JP. Integrative Genomics Viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14:178–92.Google Scholar
- Tang X, Liu G, Zhou J, Ren Q, You Q, Tian L, Xin X, Zhong Z, Liu B, Zheng X, Zhang D, Malzahn A, Gong Z, Qi Y, Zhang T, Zhang Y. A large-scale whole-genome sequencing analysis reveals highly specific genome editing by both Cas9 and Cpf1 (Cas12a) nucleases in rice. https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA420933.
- Tang X, Liu G, Zhou J, Ren Q, You Q, Tian L, Xin X, Zhong Z, Liu B, Zheng X, Zhang D, Malzahn A, Gong Z, Qi Y, Zhang T, Zhang Y. A large-scale whole-genome sequencing analysis reveals highly specific genome editing by both Cas9 and Cpf1 (Cas12a) nucleases in rice http://bigd.big.ac.cn/bioproject/browse/PRJCA000656.