CRISPR/Cas9-mediated genome editing efficiently creates specific mutations at multiple loci using one sgRNA in Brassica napus

CRISPR/Cas9 is a valuable tool for both basic and applied research that has been widely applied to different plant species. Nonetheless, a systematical assessment of the efficiency of this method is not available for the allotetraploid Brassica napus—an important oilseed crop. In this study, we examined the mutation efficiency of the CRISPR/Cas9 method for 12 genes and also determined the pattern, specificity and heritability of these gene modifications in B. napus. The average mutation frequency for a single-gene targeted sgRNA in the T0 generation is 65.3%. For paralogous genes located in conserved regions that were targeted by sgRNAs, we observed mutation frequencies that ranged from 27.6% to 96.6%. Homozygotes were readily found in T0 plants. A total of 48.2% of the gene mutations, including homozygotes, bi-alleles, and heterozygotes were stably inherited as classic Mendelian alleles in the next generation (T1) without any new mutations or reversions. Moreover, no mutation was found in the putative off-target sites among the examined T0 plants. Collectively, our results demonstrate that CRISPR/Cas9 is an efficient tool for creating targeted genome modifications at multiple loci that are stable and inheritable in B. napus. These findings open many doors for biotechnological applications in oilseed crops.


Results
CRISPR/Cas9 sgRNA design and vector construction. To apply CRISPR/Cas9 in B. napus, the sgRNA-Cas9 vectors from a previous report were used 34 . The sgRNAs were designed using the online tool CRISPR-P (http://cbi.hzau.edu.cn/cgi-bin/CRISPR) 35 . In one construct, two sgRNAs for each gene with a high score were selected and their expression was driven by the Arabidopsis U6-26 and U6-29 promoter, respectively (Supplemental Fig. S1A). The expression of Cas9 was driven by the CaMV 35S promoter. We planned to test a total of 12 B. napus genes involved in the regulation of plant development with diverse functions that belong to 4 gene families (Supplemental Table S1, and Supplemental Fig. S2 to S5). BnaA9.RGA, BnaC9.RGA, BnaA6.RGA, and BnaC7.RGA are paralogous genes of the BnaRGA family, orthologs of Arabidopsis REPRESSOR OF GA1-3 (RGA) gene, which acts as a master repressor in gibberellic signaling. Thus, loss-of-function mutations will cause phenotypes that mimick high level of GA 36,37 . The BnaFUL family contains three paralogs: BnaA9.FUL, BnaC2.FUL, and BnaC7.FUL, which are the orthologs of Arabidopsis FRUITFULL (FUL) regulating silique dehiscence during flower development [38][39][40] . BnaA2.DA2.1, BnaA2.DA2.2, BnaC6.DA2, BnaC5.DA1 and BnaA6.DA1, are paralogous genes of the BnaDA2 and BnaDA1 family, respectively, which are orthologs of Arabidopsis DA2 and DA1 (DA: LARGE IN CHINESE). The da2 and da1 mutants exhibit increases in organ size, consistent with these two genes serving as negative regulators of organ size 41,42 . Because BnaA2.DA2.1 (BnaA02g18880D) and BnaA2.DA2.2 (BnaA02g18890D) are adjacent genes with extremely high sequence identity (98.22%), no sgRNAs could be designed to distinguish between them. Therefore we considered BnaA2.DA2.1 and BnaA2.DA2.2 as one gene in our study, named BnaA2.DA2. To target one gene, two sgRNAs at gene-specific regions were designed. To target the paralogous genes of one gene family, two sgRNAs at the conserved regions were designed. We made a total of 10 constructs. Seven target one gene (Table 1) and 3 target the paralogs in one gene family ( Table 2). All of these constructs were used to transform oilseed callus following standard procedures (Supplemental Fig. S1B).
Mutation efficiency of CRISPR/Cas9 for single-gene targeted sgRNAs in B. napus. First, we evaluated the mutation efficiency of CRISPR/Cas9 for single-gene targeted sgRNAs in B. napus. The two sgRNAs used for one gene usually have different mutation rates 28,30 . Therefore, using two sgRNAs for each gene assured a high mutation rate. Because the paralogous genes have high sequence similarity in the allotetraploid B. napus, we confirmed the specificity of the primers for genotyping by direct sequencing (data not shown). A total of 233 Cas9-positive T0 transgenic lines for the 7 constructs were identified (Table 1). To accurately calculate the mutation rate, each target site was directly sequenced from the transgenic plants. We manually checked the sequencing chromatograms for each line. We concluded that mutations had been successfully introduced when the sequencing chromatograms indicated a nucleotide change (insertion, deletion or substitution) or multiple traces (overlapping peaks) at the sgRNA target sites (Supplemental Fig. S1C). Among the sgRNAs, the mutation rate varied from 5.3% to 100.0% (Table 1 and Supplemental Table S2). The average mutation rate was 65.3% (Table 1) is comparable to other plant species 9 . Among the plants transformed with 7 different constructs, no homologous mutations were found for 5 constructs. The percentage of homologous mutations was 2.4% and 6.25%, respectively, for only two sgRNA target sites (BnaA2.DA2-sgRNA2 and BnaA6.DA1-sgRNA2) ( Table 1).

Mutation efficiency of CRISPR/Cas9 for multiple-gene targeted sgRNAs in B. napus.
Simultaneously mutating paralogous genes in one gene family is important in allotetraploids, such as B. napus. To test the efficiency of simultaneous gene mutations in B. napus, two sgRNAs targeting the conserved sequences among gene family members were designed for three gene families (BnaRGA, BnaDA2 and BnaFUL) ( Table 2  and Supplemental Table 2), and a total of 67 Cas9-positive transgenic lines were created. The mutations were first detected using the T7 endonuclease I (T7E1) assay from 8 independent transgenic lines of BnaRGA-sgRNA (Supplemental Fig. S6). In the T7E1 assay, DNA fragments with mutations were digested by the T7E1 enzyme, whereas DNA fragments without mutations were not digested. A high mutation rate of 87.5% (7/8) occurred at the target sites of BnaA9.RGA, BnaC7.RGA and BnaC9.RGA (Supplemental Fig. S6). The DNA fragment from one line (L43) was not digested by T7E1, which is consistent with an intact target site. The results were further verified by Sanger sequencing, indicating that the T7E1 assay works well for Brassica (Supplemental Table S3).  We used Sanger sequencing of PCR amplicons to analyze transgenic plants when the T7E1 assay couldn't distinguish a homozygous mutation from a non-mutation 28 . The mutation frequency at each target site for sgRNA-mediate multiple-gene targeted mutagenesis ranged from 27.6% to 96.6% (Table 2, and Supplemental  Table S2), which was similar to sgRNA-mediated single-gene targeted mutagenesis (Table 2 vs. Table 1). Double, triple or quadruple mutations were readily detected, accounting for 3.4% to 86.2% of the mutants (Table 2). These data indicate that these sgRNAs efficiently target more than one locus in B. napus. A proportion of the homozygous mutations (5.9% to 20.7%) could also be identified in the T0 plants ( Table 2). The homozygous mutation rate was a bit higher than we observed for the single-gene targeted sgRNAs, which may be related to the targeting efficiency of different sgRNAs 9 .
Previous studies indicate that GC content may influence the efficiency of sgRNAs and that higher GC content is usually associated with higher mutation frequencies [43][44][45] . In this study, the GC content was indeed positively correlated with mutation frequencies for the two sgRNAs targeting the same gene (Supplemental Table S2). However, we observed that this rule was occasionally not followed, such as in the case of BnaC7.RGA-sgRNA and BnaDA2-sgRNA (Supplemental Table S2).
Variety and frequency of mutations caused by CRISPR/Cas9. Next, the sequencing data from all of the target sites was combined and analyzed to determine the mutation types and frequencies induced by CRISPR/ Cas9 in B. napus. In total, the sequencing results of 422 PCR amplicons were analyzed by the decoding website DSDecode (http://dsdecode.scgene.com/) 46 . Part of the results were further confirmed by TA cloning and sequencing. The results are summarized in Fig. 1, and more details are provided in Supplemental Data S1. Several types of mutations were observed: deletions, insertions, substitutions, and combined mutations (i.e., more than one mutation type in one allele). Among all of the types of mutations, 53.5% were deletions, 42.3% were insertions, 2.9% were combined mutations, and 1.3% were substitutions ( Fig. 1A and B). The mutations were predominantly short nucleotide changes (≤3-bp) (62.2%) (Fig. 1C), a majority of which (41.7%) were one nucleotide insertions (Fig. 1B). The length of deletions ranged from one bp to hundreds of bp, and 4.2% of mutations exhibited a >100-bp deletion. The longest deletion was 270 bp (Supplemental Data S1). These long fragment deletions are caused by the simultaneous repairing of two DSBs generated by two sgRNAs in one construct in our system 33 .
The cleavage site of Cas9 is usually 3-bp upstream of the PAM sequence 47 . Our results demonstrated that 91.3% of the mutations indeed occurred at this position. Additionally, 0.7%, 1.4%, 4.3% and 2.2% of the mutations occurred at the 1 st , 2 nd , 4 th and 5 th base from the PAM site, respectively (Supplemental Fig. S7A). In more detail, all of the 1-bp insertions (100%) and most of the 1-bp deletions (87.5%) were located 3-bp upstream of the PAM site (Supplemental Fig. S7B). When the base composition of the 1-bp insertions was examined, most of them were A (44.5%) or T (32.8%) insertions. Unexpectedly, the percentage of C insertions was 21.1%, which is much higher than was observed in rice (7.6%) 30 and tomato (9.3%) 28 (Fig. 1D).
We further compared the mutation types for each target gene. For example, all BnaC9.RGA-sgRNA1 mutations were deletions. Mutations in BnaA9.RGA-sgRNA1 were predominantly insertions (91.7%). The BnaA6.RGA-sgRNA1 induced substitutions at a frequency of 21.4%, which is much higher than other sgRNAs (Fig. 1E). In general, short indels (≤3-bp) were most abundant. However, the deletion length was more than 4-bp in the entire BnaC9.RGA-sgRNA1 target sites (Fig. 1E). We also observed different types of mutations for different loci targeted by the same sgRNA. For example, BnaFUL-sgRNA2 induced 81.25% deletions in BnaA9.FUL and 77.8% insertions in BnaC7.FUL (Fig. 1E). Although there were a limited number of mutations analyzed for some targets, these results provide strong evidence that the types of mutations vary at distinct target sites.
Genotypes of CRISPR/Cas9 mutants in the T0 generation. There are two alleles for each gene, both of which might be mutated by CRISPR/Cas9. Thus, CRISPR/Cas9 could produce five genotypes: homozygote (the two alleles have the same mutation), bi-allele (the two alleles have different mutations), heterozygote (only one allele is mutated), chimera (more than two different mutations exist), and WT-type (no mutation). To estimate the proportion of each genotype among the T0 mutants, we randomly chose 117 single-gene targeted and 66 multiple-gene targeted sgRNA lines containing Cas9 insertions to analyze the mutations of each targeted site by Sanger sequencing. A total of 425 PCR amplicons were decoded by the DSDecode website. Additionally, 14 PCR amplicons were analyzed by inserting them into a TA vector and sequencing 10 individual clones for each amplicon. We obtained results from 439 amplicons. The genotype data are summarized in Table 3.
Based on the genotyping results, 7.7% (34/439) sites were homozygous and 20.5% (90/439) sites were bi-allelic. Thus, a total of 28.2% had defects in both alleles (Table 3). Among all the T0 homozygotes, 70.6% (24/34) of them carried a 1-bp insertion (i1) or a 1-bp deletion (d1), indicating that these two mutation types happened at a high frequency in T0 homozygous plants. For the bi-allelic mutations, the predominant type carried a combination of i1ai1b mutations (a and b indicate two different nucleotides) (19.4%), followed by i1d1 (i: insertion, d: deletion) (12.0%) (Supplemental Fig. S8). The frequency of heterozygotes and chimeras was 16.1% (71/439) and 31.9% (140/439), respectively. No mutations were found in 23.7% (104/439) of the sites. Interestingly, we identified one line called BnaC5.DA1-sgRNA1-L16 that has an A/T insertion at the expected target site without the Cas9 insertion ( Fig. 2A Table S3). Based on these data, we suggest that it is possible to get T-DNA-free plants in the first generation by regenerating oilseed callus.
Previous studies showed that different CRISPR/Cas9-induced mutations could arise from different tissues 27,29 . To test this possibility in B. napus, we genotyped 6 single-gene targeted and 9 multiple-gene targeted sgRNA lines with mixed tissues of leaves, shoots, and flower buds. A total of 21 PCR amplicons were used for the sequencing analysis. For homozygous and bi-allelic lines, the same mutations were detected in leaves and mixed tissues (Supplemental Tables S3 and S4). However, we found that 36.4% (4/11) of heterozygous and chimeric plants had new mutations in shoots and flower buds (Supplemental Tables S3 and S4), indicating that the wild type alleles from these plants were further mutated by CRISPR/Cas9. When 3 bp or multiples of 3 bp are deleted from an exon without affecting the reading frame, a few amino acids may be deleted from the middle of a protein, which may alter its biochemical properties. BnaA6.RGA encodes a DELLA protein that negatively regulates GA signaling in many plant species. It is reported that deletion of the conserved TVHYNP domain resulted in a dwarf phenotype and insensitivity to exogenous GA 36,37 . Among the transgenic lines of BnaA6.RGA-sgRNA, two T0 lines (T0-L4 and T0-L6) were dwarf (Fig. 3A), and the phenotype was inherited in the T1 generation (Supplemental Fig. S10). According to the sequencing results, the T0-L4 and T0-L6 plants have a 6-and 12-nt deletion at the target sites, respectively (Fig. 3B, Supplemental Table S3), that causes a 2-and 4-amino acid deletion in the TVHYNP domain (Fig. 3C).      (Table 3) were examined for the genotypes at the target sites (Table 4; Supplemental Table S3). Due to the large number of samples, sequencing results from the entire set of DNA amplicons were directly decoded by DSDecode. Therefore, only homozygotes, wild type, bi-alleles and heterozygotes were clearly identifiable. We could not accurately decode chimeras using this approach. Thus, chimeras are indicated with an 'h' for heterogeneous ( Table 4). The T1 progeny of 3 T0 homozygotes were still homozygous with the same mutations, indicating that the mutations in these transgenic lines were stable (Table 4). If the bi-allele and heterozygous genotypes are inherited normally, the segregation ratio is expected to be 1xx:2xy:1yy. However, unexpected segregation ratios of 0:2:13, 3:12:0, 0:20:0 and 11:1:2 were observed among the T1 progeny of BnaDA2-sgRNA2-L12 (targeted to BnaA2.DA2), BnaDA2-sgRNA1-L20 (targeted to BnaC6.DA2), BnaRGA-sgRNA1-L27 (targeted to BnaA9.RGA) and BnaC5.DA1-sgRNA-L1, respectively ( Table 4), which might be due to unequal inheritance frequencies of the two alleles. New mutation types were identified in the T1 generation for heterozygotes. One possible reason is that Cas9 could be still functional at the non-mutated allele at the targeted region. As expected, the T1 segregation patterns of T0 chimeras were more diverse and less predictable ( Table 4). Irrespective of bi-allelic, heterozygous or chimeric mutants, homozygotes were found in T1 generation (Table 4). Unexpectedly, no T1 descendants of WT showed any novel mutations even in the presence of the Cas9 transgene (Table 4), indicating that the CRISPR/Cas9 was not functional in these transgenic lines, possibly due to lower expression of Cas9 and/or the guide RNA. Because a majority of genes have multiple paralogs with functional redundancy in B. napus (allotetraploid), single gene knockouts probably will not show obvious phenotypes 48 . In this study, the T1 quadruple mutant of BnaRGA grew longer stems than control plants ( Fig. 3D and Supplemental Fig. S11). Genotyping results indicated that the four paralogs of BnaRGA could be efficiently knocked out with homozygous or bi-allelic mutations in the T0 generation (Fig. 3E), and inherited in the T1 generation (Supplemental Table S3). Together, these data make the case that CRISPR/Cas9 has great advantages for gene function studies in B. napus.
Genetically manipulated materials without T-DNA insertions are largely favored for crop improvement and should be more public acceptable. Indeed, the T-DNA insertions in 76.2% of the T0 lines did not cosegregate with the CRISPR/Cas9-induced mutations and were therefore removed in the next generation. The average value for  T1 progeny lacking the Cas9 transgene was 10.9% when analyzed using Cas9-specific primers (Supplemental Fig. S12, Supplemental Table S5). The homozygote and bi-allele genotypes were stably passed to subsequent generations regardless of whether the T-DNA was present. These data indicate that CRISPR/Cas9 is an effective tool for the improvement of B. napus.

No off-targets were discovered in B. napus. Low-frequency examples of off-target cleavage have been
reported for CRISPR/Cas9 in plants 28,30 . To detect the off-target events in oilseed plants, potential off-target loci following PAM sequences that are highly homologous to the sgRNAs of BnaA9.RGA, BnaC9.RGA, BnaA6.RGA and BnaC7.RGA were predicted using the online tool CRISPR-P (http://cbi.hzau.edu.cn/cgi-bin/CRISPR) 35 . At least three of the most likely off-target sites for each sgRNA were examined in a total of 50 randomly selected T0 and T1 plants using gene-specific primers (Supplemental Table S6). Previous reports have indicated that the 12 nucleotides of "seed sequence" located in the target site and adjoining the PAM are critical for recognition specificity and cleavage efficiency of Cas9 30,49,50 . In the off-target sequences, there were 1 to 3 mismatches in the 'seed sequence' . No mutations were found in the putative off-target sites ( Table 5), indicating that CRISPR/ Cas9-mediated mutagenesis is specific in oilseed plants.

Discussion
Highly efficient target gene mutagenesis by CRISPR/Cas9 in B. napus. Compared to traditional mutagenesis strategies, CRISPR/Cas9 targeted genome editing is precise and efficient. Therefore, CRISPR/Cas9 is extremely useful for gene function studies and crop improvement 25,30,51,52 . In this work, to systematically assess the application of CRISPR/Cas9 in B. napus, 12 genes were selected for targeted mutagenesis. The mutation frequency ranged from 5.3% to 100%. Mutation frequencies are similar in most of plants. Based on these data, we suggest that variations in genome size do not significantly influence the efficiency of targeted genome editing mediated by the CRISPR/Cas9 system 9 . Indeed, we found that CRISPR/Cas9 efficiently created homozygous and bi-allelic mutations that were stably maintained during plant regeneration in B. napus, which has a high ploidy level. CRISPR/Cas9 induced these two genotypes not only in one target gene, but also in several paralogous genes without any reduction of efficiency (Tables 3 and 4). Due to a high ploidy level, paralogs in one gene family with functional redundancy usually exist in B. napus. Therefore, one-gene knockouts do not lead to phenotypes in B. napus. For example, we observed no significant differences in the single-gene knockout mutant of BnaA6.RGA-sgRNA and BnaA9.RGA-sgRNA relative to WT plants (Supplemental Fig. S13). In contrast, we observed significantly longer stems in the quadruple mutant of BnaRGA-sgRNA relative to WT plants ( Fig. 3D and Supplemental  Fig. S11). Our results demonstrate that using sgRNAs derived from conserved regions, CRISPR/Cas9 can simultaneously knockout a group of paralogous genes. Thus, CRISPR/Cas9 is able to create ideal materials for functional studies in oilseed rape. Based on a comparison of T0 and T1 generation plants (Table 4), it is clear that the mutations of homozygotes and most of the bi-alleles are stably inherited regardless of whether T-DNA is present. However, heterozygotes retain the wild type alleles ( Table 4) that have the potential to be mutated. Thus, new mutation types could arise from different tissues in both T0 and T1 plants probably due to cell-autonomous mutagenesis in those tissues (Supplemental Table S3), consistent with other reports 28,30 . In the meantime, the wild type alleles in the heterozygotes and chimeras could also be mutated as plants are propagated (Table 4). It is noticeable that no new mutations were generated for the WT plants when the T0 plants were propagated to the T1 generation (Table 4), which is inconsistent with the high mutagenesis efficiency of CRISPR/Cas9. We assume that the expression level of CRISPR/Cas9 is low because of position effects, post-transcriptional gene silencing, or the Cas9 protein is somehow inactivated 9,30 .  Cas9 and sgRNAs affect the genome editing efficiency. One construct expressing two sgRNAs for each targeted gene assured a high mutation rate per gene 9 . A comparison of mutation frequencies obtained using this strategy indicated that sgRNAs with higher GC contents are usually associated with higher mutation frequencies, with the exception of BnaC7.RGA-sgRNA and BnaDA2-sgRNA (Supplemental Table S2). Our findings indicate that GC content may influence the efficiency of sgRNA, which is consistent with previous work [43][44][45] . Unexpectedly, the mutation rate of each gene is usually determined by the sgRNA with the higher mutation rate. In other words, sgRNAs that induce higher mutation rates appear to compensate for sgRNAs that induce lower mutations rates. The expression level and pattern of Cas9 and sgRNAs driven by different promoters have a major influence on genome editing efficiency 9 . In Arabidopsis and maize, tissue specific or plant endogenous promoters that drive Cas9 expression greatly increased the mutation frequency compared to constitutive promoters, such as the CaMV 35S promoter 27,[53][54][55][56] . The results from soybean and liverwort demonstrated that mutation frequencies could be increased 2 to 7-fold when the intrinsic U6 promoter is used to drive sgRNA expression, compared to the Arabidopsis U6 promoter 57,58 . In contrast, high mutation frequencies were observed in B. napus when using either constitutive promoters or the Arabidopsis U6 promoter. Although these promoters work well for generating mutations in B. napus, particular sites were edited at a low frequency (e.g., BnaC7.RGA-sgRNA1, 5.3%). Tissue-specific promoters or more active promoters may increase the frequency of mutagenesis by increasing the expression level of Cas9 and sgRNAs during plant regeneration.
Potential T-DNA free mutants were generated by CRISPR/Cas9 in B. napus. T-DNA-free mutants could be generated by self-crossing or backcrossing in the T1 generation, which provides reliable material for crop improvement 28,30,58 . Indeed, 11.3% of the T1 mutant plants lacked the Cas9 transgene (Supplemental Table S5). Additionally, we found one Cas9-negative T0 line-BnaC5.DA1-sgRNA-L16-with a mutation at the target site in both the T0 and T1 generations ( Fig. 2 and Table 4). One possible explanation for this finding is that the CRISPR/ Cas9 cassette was lost during cell division or that transient expression of Cas9 is responsible for this mutation. Three pairs of primers were used to confirm the loss of the T-DNA insert (Supplemental Fig. S12). However, whole genome resequencing would be a perfect approach to rule out the presence of any T-DNA fragments 59 . For some lines, very few or no plants lost the CRISPR/Cas9 cassette, which may be due to multiple CRISPR/Cas9 cassette insertions in the genome. In this case, increasing the population size or a cross with WT is needed to remove the transgenes. Recently, two approaches have been developed for generating Cas9-free mutants in one generation. One strategy involves delivering the Cas9-sgRNA ribonucleoprotein complex or the CRISPR-Cas9 DNA/RNA into plant cells using particle bombardment or protoplast transfection 60,61 . The other strategy uses fluorescent proteins as markers to facilitate the selection of transgene-free plants 62 . These reports provide new strategies screen and isolate Cas9 free material and speed up molecular breeding in the future.
In summary, we demonstrated that the CRISPR/Cas9 is a highly efficient tool for genome editing in B. napus. We found that sgRNAs derived from conserved sequences could simultaneously induce homozygous or bi-allelic mutations at multiple loci. Therefore, CRISPR/Cas9 provides a powerful tool for studying gene function in oilseed and the fastest method for the breeding of polyploid crops. Moreover, the targeted gene modification mediated by CRISPR/Cas9 is safer for both human health and the environment because it is possible to remove the foreign DNA following the mutagenesis of the target DNA. Thus, CRISPR/Cas9 is useful for both basic and applied research in B. napus.

Material and Methods
Plant materials and growth conditions. The B. napus variety Westar was transformed with Agrobaterium. The transgenic lines and wild-type plants were grown in a greenhouse at 22 °C, 70% relative humidity, and in a photoperiod containing16 h of light/8 h of dark. Mature seeds were collected from T0 plants and germinated for 7 days at 22 °C, on petri dishes, in a photoperiod containing 16 h of light/8 h of dark. The seedlings were then transferred to soil and grown in the greenhouse.
Vector construction. The sgRNA-Cas9 plant expression vectors were constructed as previously described with minor modifications 34 . The target sgRNA sequences were designed using the web server CRISPR-P (http:// cbi.hzau.edu.cn/cgi-bin/CRISPR) 35 , and then the sequences were further analyzed using the CRISPR Primer Designer software 50 . Using pCBC-DT1T2 as the template, two AtU6 promoter-sgRNA-AtU6 terminator cassettes were amplified by PCR using the primers listed in Table S8. Then the PCR fragments were inserted into pKSE401 by Golden Gate Assembly 64 , and confirmed by Sanger sequencing. The pKSE401-sgRNA vectors were used for plant transformation.
ScieNtific RepoRts | 7: 7489 | DOI:10.1038/s41598-017-07871-9 Genotyping and T7E1 assay. To analyze the mutations caused by CRISPR/Cas9, genomic DNA was extracted from each transgenic plant using the CTAB method (Molecular Cloning, 3rd edition). The flanking sequence around the CRISPR target sites was amplified by PCR using gene-specific primers (Supplemental Table S6). First, each PCR amplicon generated from WT genomic DNA was sub-cloned into the pGEM-T easy vector (A3600, Promega, USA) to confirm the primer specificity by Sanger sequencing. Then, most of the amplicons were directly sequenced to analyze the mutations. For the complex mutations, the amplicons were first sub-cloned into the pGEM-T easy vector, and about 10 clones of each amplicon were individually sequenced by Sanger sequencing.
The T7E1 assay was carried out as previously described with minor modifications 27 . 300 ng of purified PCR products were denatured and annealed in a thermocycler using the following program: 5 min at 95 °C, 95 to 85 °C using a ramp rate of −2 °C/sec, 85 to 25 °C using a ramp rate of −0.2 °C/sec. The annealed PCR products were digested with T7E1 nuclease (M0302S, New England Biolabs) which specifically cleaves DNA with mismatches at 37 °C for 1 h. Digested PCR products were analyzed using agarose gel electrophoresis containing ethidium bromide.
Off-target analysis. The potential off-target sites were predicted using CRISPR-P (http://cbi.hzau.edu.cn/ cgi-bin/CRISPR). The top-ranking potential off-target sites containing fewer than 3-bp mismatches in the 12-bp seed sequence were selected for validation. The genomic DNA sequences surrounding the potential off-target sites were amplified by PCR using gene-specific primers (Supplemental Table S6). PCR products were analyzed by direct sequencing.
Sequencing chromatogram decoding. The online tool DSDecode (http://dsdecode.scgene.com/) 46 was used for chromatogram decoding. Sequence files (xxx.abi) and the reference gene sequences were uploaded to the server where they were decoded. Subsequently, the results were compared to the reference sequence to ensure that the cleavage site is in the target region of sgRNA.