Multiplexed Cas9 targeting reveals genomic location effects and gRNA-based staggered breaks influencing mutation efficiency

Nature Communications, Apr 2019

Understanding the impact of guide RNA (gRNA) and genomic locus on CRISPR-Cas9 activity is crucial to design effective gene editing assays. However, it is challenging to profile Cas9 activity in the endogenous cellular environment. Here we leverage our TRIP technology to integrate ~ 1k barcoded reporter genes in the genomes of mouse embryonic stem cells. We target the integrated reporters (IRs) using RNA-guided Cas9 and characterize induced mutations by sequencing. We report that gRNA-sequence and IR locus explain most variation in mutation efficiency. Predominant insertions of a gRNA-specific nucleotide are consistent with template-dependent repair of staggered DNA ends with 1-bp 5′ overhangs. We confirm that such staggered ends are induced by Cas9 in mouse pre-B cells. To explain observed insertions, we propose a model generating primarily blunt and occasionally staggered DNA ends. Mutation patterns indicate that gRNA-sequence controls the fraction of staggered ends, which could be used to optimize Cas9-based insertion efficiency.

A PDF file should load here. If you do not see its contents the file may be temporarily unavailable at the journal website or you do not have a PDF plug-in installed and enabled in your browser.

Alternatively, you can download the file locally and open with any standalone PDF reader:

https://www.nature.com/articles/s41467-019-09551-w.pdf

Multiplexed Cas9 targeting reveals genomic location effects and gRNA-based staggered breaks influencing mutation efficiency

Multiplexed Cas9 targeting reveals genomic location effects and gRNA-based staggered breaks influencing mutation efficiency Santiago Gisler Joana P. Gon?alves 0 Waseem Akhtar 0 1 2 Johann de Jong 2 3 Alexey V. Pindyurin Lodewyk F.A. Wessels 1 2 Maarten van Lohuizen 0 0 Division of Molecular Genetics, Oncode and The Netherlands Cancer Institute , Plesmanlaan 121, Amsterdam 1066 CX , The Netherlands 1 Department of Intelligent Systems, Delft University of Technology , Van Mourik Broekmanweg 6, Delft 2628 XE , The Netherlands 2 Division of Molecular Carcinogenesis, Oncode and The Netherlands Cancer Institute , Plesmanlaan 121, Amsterdam 1066 CX , The Netherlands 3 Data & Translational Sciences Group, UCB Biosciences GmbH , Alfred-Nobel-Stra?e 10, Monheim am Rhein 40789 , Germany 4 Institute of Molecular and Cellular Biology, Siberian Branch of Russian Academy of Sciences , Acad. Lavrentiev Ave. 8, Novosibirsk 630090 , Russia 5 Division of Gene Regulation, Oncode and The Netherlands Cancer Institute - 5,6, :,; () 0 9 8 7 6 5 4 3 2 1 Understanding the impact of guide RNA (gRNA) and genomic locus on CRISPR-Cas9 activity is crucial to design effective gene editing assays. However, it is challenging to profile Cas9 activity in the endogenous cellular environment. Here we leverage our TRIP technology to integrate ~ 1k barcoded reporter genes in the genomes of mouse embryonic stem cells. We target the integrated reporters (IRs) using RNA-guided Cas9 and characterize induced mutations by sequencing. We report that gRNA-sequence and IR locus explain most variation in mutation efficiency. Predominant insertions of a gRNA-specific nucleotide are consistent with template-dependent repair of staggered DNA ends with 1-bp 5? overhangs. We confirm that such staggered ends are induced by Cas9 in mouse pre-B cells. To explain observed insertions, we propose a model generating primarily blunt and occasionally staggered DNA ends. Mutation patterns indicate that gRNA-sequence controls the fraction of staggered ends, which could be used to optimize Cas9-based insertion efficiency. G recent years, nurtured by the emergence of precision enome engineering has seen considerable progress in editing tools based on the bacterial clustered regularly interspaced short palindromic repeats (CRISPR)-associated protein 9. The CRISPR-Cas9 system complexes the endonuclease enzyme Cas9 with a guide RNA to induce double-strand breaks (DSBs) at a specific DNA locus1?5. For target DNA recognition and binding, CRISPR-Cas9 requires the presence of a short, conserved sequence known as protospacer-adjacent motif (PAM). The PAM consists of nucleotides NGG and is located downstream of the target sequence6,7. Cas9-induced DSBs activate the cellular DNA damage response, mainly through non-homologous end-joining (NHEJ) or homology-directed repair (HDR)8. NHEJ is the most common DNA repair pathway9. In NHEJ, DNA ends are processed independently without a template prior to ligation, often producing mutations at the break site. HDR relies on sequence homology for repair and therefore depends on the availability of a donor DNA template, which can be acquired from the sister chromatid in Sphase8. In genome engineering applications, error-prone repair of Cas9-induced DSBs can be exploited to disrupt the target sequence and generate gene knockouts. Exogenous genetic material can also be integrated into host DNA by providing repair templates with custom oligonucleotides flanked by homology arms. The CRISPR-Cas9 technology is used extensively for gene editing in vitro and in vivo, yet most factors controlling its nuclease activity are poorly understood. While effects of guide RNA on Cas9 nuclease efficiency and target specificity have been extensively characterized10?14, the influence of target sequence on induced mutation patterns remains unclear. Little is also known on the impact of genomic and epigenomic context at the target locus15. Early studies found that chromatin accessibility or DNA methylation affect the binding of catalytically inactive Cas9 (dCas9)16?18. Others showed that Cas9 binding and cleavage are sensitive to chromatin changes induced by nucleosome occupancy19,20 or administration of doxycycline21. Most literature suggests that genomic context influences Cas9 binding and cleavage, but effects on editing efficiency are less well understood. Several studies have observed weak correlations between epigenomic context and Cas9-induced mutation frequency at endogenous targets17,22?24. In particular, two of these studies showed that the significant effect of epigenomic context on Cas9 binding did not necessarily result in a detectable effect on Cas9-induced mutation frequency17,22. Here, we characterize Cas9-induced mutations in the genomes of mouse embryonic stem (mES) cells. We aim to survey many loci for sufficient statistical power with minimal disruption of the native environment. However, it is not trivial to scale up the number of endogenous Cas9 targets. Compromising on guide RNA specificity increases off-target effects. Targeting repetitive sequences creates challenges for alignment and mutation calling, and results might not generalize to other kinds of sequences. Both approaches can generate a large number of cleavage events per cell, eventually leading to genomic instability and unreliable Cas9 activity profiling as a result. While multiple guides could instead be used to expand Cas9 targeting, this would also introduce target heterogeneity. Alternatively, we integrate thousands of barcoded target sequences throughout the genomes of a population of mES cells using our TRIP technology25,26. In this way, we multiplex Cas9 cleavage while keeping the number of targets per cell under control. The use of TRIP reporters further enables sequence-independent analysis of effects across the targeted loci. We investigate the usefulness of these hybrid exogenousendogenous loci to profile Cas9 activity, and assess the impact of guide RNA sequence and targeted locus on induced mutation frequency and patterns. Results RNA-guided Cas9 targeting of integrated reporters. We profiled CRISPR-Cas9-induced mutations across the genomes of mES cells. First, we used TRIP25,26 to embed barcoded reporter genes randomly throughout the host DNA (Fig. 1a). We established a clonal TRIP cell line containing 36 PGK-driven integrated reporters (IRs) per cell, and a multi-promoter TRIP pool with ~1k IRs distributed heterogeneously across cells (Fig. 1b). We designed three single-guide RNAs (sgRNAs) targeting sites near the 3?-end of the IR gene body, cloned them into Cas9sgRNA plasmids and used them in independent assays (Fig. 1c). After selecting Cas9-sgRNA-carrying cells, we amplified and sequenced IR target regions to characterize induced lesions. As proof-of-concept, we performed Cas9 disruption assays using sgRNA1-3 in the TRIP cell line (Fig. 1c). By disruption, we refer to cleavage without the use of exogenous DNA. We also did disruption assays in TRIP pools to study mutations at a large number of loci. In addition, we performed editing involving the knock-in of a 21-nt single-stranded oligodeoxynucleotide (ssODN) to characterize template-dependent insertions. We used both sgRNA2/3 for disruption in TRIP pools (-ssODN), and only sgRNA2 for editing ( + ssODN) since the proximity between the sgRNA3 target and the IR barcode prevented the design of proper homology arms. Finally, we analyzed all 36 IRs in the cell line, and also the 1359 IRs with at least 30 reads in all pool assays. Variation in Cas9-induced mutation frequency across IR loci. We first analyzed Cas9-induced mutation frequencies at the targeted loci, and the effects of different factors on those frequencies. Mutation frequency was determined per IR as the fraction of reads carrying a mutation amongst all reads mapped to that specific IR. Overall, Cas9-targeted IR sequences showed high mutation frequencies genome-wide in both cell lines and pools. Cell line averages were ~50% for sgRNA1/2 and ~65% for sgRNA3, while pools reached ~30% for sgRNA1/2 and ~60% for sgRNA3 (Figs. 2a, b). We examined how Cas9-induced mutation frequencies varied with sgRNA, IR locus, ssODN, and promoter. Guide RNA sequence-driven variation in mutation frequency. In line with previous reports, different guide RNAs led to systematic variation in IR mutation frequencies1,13,27. Guide sgRNA3 was most efficient in the cell line with average 1.30-fold and 1.26-fold increases in mutation frequency relative to sgRNA1/2 (effect sizes 15% ? 1% for sgRNA3 vs. sgRNA1 and 13% ? 1% for sgRNA3 vs. sgRNA2, both with p ? 2.91 ? 10?11, two-tailed Wilcoxon rank sum tests, Fig. 2a). The difference was largely due to insertions, showing 12.5-fold and 5.26-fold increases with sgRNA3 compared to sgRNA1/2. Deletion frequency was more comparable across guides, varying by 1.23-fold and 1.10-fold between sgRNA3 and sgRNA1/2. Guide sgRNA3 was also the most efficient in TRIP pools, promoting an average 1.96-fold increase in mutation frequency (effect size 30% ? 1% sgRNA3-ssODN vs. sgRNA2-ssODN, p < 2.20 ? 10?16, two-tailed Wilcoxon signed rank test; Figs. 2b, c). Again, insertions increased by 6.76-fold while deletions varied by 1.45-fold. Overall, we found that sgRNA1-3 resulted in different levels of mutation efficiency, mostly contributed by insertions. Locus-associated variation in mutation frequency. We questioned whether certain IR loci would be more prone to mutations than others. We saw that IR mutation frequency correlated strongly across cell line assays using the three sgRNAs, with R2 > 0.85 and Promoter EGFP reporter Barcode mPGK-driven TRIP library F-test p < 3.13 ? 10?16 (Fig. 2a), and across TRIP pool assays with R2 > 0.72 and F-test p < 2.16 ? 10?16 (Fig. 2d). Correlations were lower for insertions (R2 0:30) than deletions (R2 0:72), likely due to the scarcity of insertion events. The highly reproducible mutation frequencies revealed consistent locus-specific susceptibility to Cas9-induced mutations (Figs. 2a, d). Knock-in and error-based insertion frequency. We examined the frequencies of error-based and knock-in insertions. Most insertions induced by Cas9-sgRNA2 in the TRIP pool resulted from errors in endogenous repair of Cas9-induced DSBs, with average frequencies per IR of 3.78% and 3.24% in disruption and editing experiments respectively (Fig. 2b). Knock-ins occurred only in editing assays with the integration of the designed 21-nt ssODN at the break site by HDR. Knock-in efficiency was low, as expected15, with an average of 1.74% per IR (Fig. 2b, left plot). However, knock-ins were more frequent than error-based insertions at IRs with total mutation frequency larger than 70% (Fig. 2e). Effects on mutation frequency. To quantify the effect of the above factors on IR mutation frequency, we modelled mutation sgRNA1 sgRNA2 100% frequencies in the TRIP pool as a linear function of IR locus, guide RNA, ssODN inclusion, promoter, and an interaction term for the joint (non-additive) contribution of locus and guide (Figs. 2f, g). The linear regression model yielded a goodness-of-fit of R2? 0.98 (Fig. 2f). Using multi-way ANOVA tests, we determined the effect size of each factor in the model (Fig. 2g). IR locus explained ~79.5% of the variation in mutation frequency, while guide RNA sequence was responsible for ~12.7%, and locus-guide interaction accounted for ~4.9%. ssODN and promoter had negligible effect, with less than 1% together. These results confirmed that IR locus and guide RNA are major determinants of mutation frequency. Note that IR locus encapsulates a variety of factors that make a locus unique, including genomic context or the interaction between an IR and the host DNA. Association between genomic context and mutation frequency. We sought to analyze the contribution of genomic context at IR loci to Cas9-induced mutation frequencies. Specifically, we examined the relation between the mutation frequency for IRs in TRIP pools (Fig. 3a) and transcriptional, genomic and epigenomic (TGE) features (Figs. 3b?d). Transcriptional features included IR expression in our TRIP cells, and gene expression and transcription factor binding in wild-type mES cells28. Genomic features comprised metrics such as GC content and a d 200 150 ) p b n o lili m ( ion 100 it s o p c i m o n e G 50 0 IR expression, Mutation frequency (per type, assay) IR expression Oligo insertion +ssO sgRNA2 Insertion +ssO sgRNA2 Insertion ?ssO sgRNA2 Insertion ?ssO sgRNA3 Deletion +ssO sgRNA2 Deletion ?ssO sgRNA2 Deletion ?ssO sgRNA3 +ssODN sgRNA2 ?ssODN sgRNA2 Epigenome Expression Genome Regulation (RNA polll) Regulation (epigenetic) Not significant p > = 0.05 Pearson?s correlation ?0.5 0.0 0.5 1.0 gene proximity in mES cells28. Epigenomic features included chromatin density from Hi-C assays, and chromatin immunoprecipitation (ChIP) data for numerous histone modifications, DNaseI hypersensitivity, and Lamin-B1 in mES cells28. We quantified TGE features within a region of 2 kb surrounding each IR locus and calculated their correlation with IR mutation frequency. Mutation frequency correlated weakly with TGE features (Pearson?s r = [?0.123, 0.151], (Fig. 3b?d)), but insertions showed stronger association than deletions (effect size 0.021 ? 0.007, p ? 1.91 ? 10?9, two-tailed Wilcoxon signed rank test sgRNA2 ? ssODN and sgRNA3-ssODN insertions vs. sgRNA2 ? ssODN and sgRNA3-ssODN deletions). Correlations were also higher for sgRNA3 compared with sgRNA2, more so for insertions than deletions (effect size 0.022 ? 0.008, p ? 2.67 ? 10?7 for insertions, 0.016 ? 0.009, p? 0.0001 for deletions, twotailed Wilcoxon signed rank test sgRNA3-ssODN vs. sgRNA2ssODN). We investigated whether changes induced by TRIP IRs could have disrupted the correlation between IR mutation frequency and expression-based features. For this we analyzed IR expression, quantified as IR barcode counts in cDNA normalized by IR barcode counts in genomic DNA obtained by highthroughput sequencing of a 168 bp region spanning the Cas9 target sites in TRIP mES cells. The correlation with mutation frequency was similarly weak for IR expression than for expression in wild-type mES cells, and we therefore excluded a potential impact of TRIP on these results (Fig. 3d, Supplementary Figure 3). Also reassuring was the fact that IR expression correlated strongly with wild-type TGE features (Pearson?s r = [?0.603, 0.648], Fig. 3d). Additionally, we looked into translocations upon Cas9 targeting as an eventual source of variation in mutation frequency, and found no evidence that they played a role in this regard (Supplementary Tables 1 and 2). Expression of IRs and endogenous genes showed the highest association with mutation frequency amongst six categories of TGE features (Fig. 3c). Moreover, features related to the RNA PolII complex involved in genome-wide transcription initiation produced larger effects than transcription factors. Individually, phosphorylated RNA PolII (Ser2P, Ser5P, Ser7P) and influencers of transcription initiation such as Ctr9, Taf3 and CpG islands correlated positively with insertion frequency (Fig. 3d). We saw similar effects for H3K4 methylation and other histone modifications (H3K9ac, H3K27ac, H3K36me3, H3K79me2) associated with actively transcribed regions29. Correlating negatively with insertion frequency were known indicators of transcription inhibition, namely: lamina-associated domains (Lamin-B1), heterochromatin mark H3K9me2, and the Polycomb Repressive Complex 1 (PRC1) member Cbx7 which promotes trimethylation of H3 at Lys-9 (H3K9me3)30. Characterization of Cas9-induced mutation sizes and patterns. In addition to mutation frequencies, we analyzed mutation sizes and patterns. We found that Cas9 targeting yielded mostly small mutations. On average, 73.4% of all deletions in the cell line were smaller than 10-bp (Fig. 4a). The most common sizes were {1,2}bp with sgRNA2 and {2,4,5}-bp with sgRNA3, each accounting for approximately 12 to 14% of all deletions. In contrast, sgRNA1 led to a large number of 3-bp deletions (41.8%, Fig. 4a), the majority of which denoted the loss of triplet CGG, likely at positions 1?3 or 4?6 upstream of the PAM (40.2%, Fig. 4b). Single-nucleotide insertions accounted for > 78% of all insertions with sgRNA1/2, and > 95% with sgRNA3 (Fig. 4a). We found similar size and pattern distributions in experiments using varying Cas9 concentrations (Supplementary Figure 6). We observed that deleted regions neighbored the expected Cas9 target site between nucleotides 3|4 upstream of the PAM with sgRNA1/2 (Fig. 4c). This indicates that resection might preferentially occur on one rather than both DNA ends at the break site. With sgRNA3, deletions often neighbored nucleotides 4|5 instead. In addition, each sgRNA led to predictable 1-bp insertions, namely G (63.7%) for sgRNA1, G (73.9%) for sgRNA2, and T (97.9%) for sgRNA3 (Fig. 4d). Deletion borders and 1-bp insertions were therefore highly consistent and guide-specific. Insight into Cas9 cleavage based on insertion patterns. We sought to understand the observed preference for specific nucleotide insertions under current Cas9 cleavage models (Figs. 4d and 5a). Cas9 is thought to primarily induce bluntended DSBs3,31. Besides direct re-ligation, blunt ends may be processed leading to a deletion or, crucially, templateindependent addition of nucleotides9. However, evidence from the seminal work on CRISPR-Cas93,31 and recent simulations on Cas9-domain conformation32 suggests that Cas9 may also generate staggered DSBs. Specifically, Cas9 domains RuvC and HNH could cleave between nucleotides 3|4 on the target DNA (tDNA) and 4|5 on the non-target DNA (ntDNA) upstream of the PAM32. The resulting 5? overhangs could trigger polymerasebased fill-in at position 4, producing sgRNA-specific insertions consistent with our data (Figs. 4d and 5a). Plausible Cas9 DNA cleavage and repair models. We assessed three different models of Cas9 cleavage and DNA repair (Fig. 5a). First, a blunt model producing blunt-ended DSBs primarily at 3| 43,31, and leading to insertions of a random nucleotide9. Second, a staggered model cleaving at 3|4 (tDNA) and 4|5 (ntDNA)32 and inducing the replication of nucleotide 4. Third, a combination of both blunt and staggered models. Each model determined an expected distribution of insertion counts across sites in the target sequence. In order to compare model distributions against observed insertions, we also addressed uncertainties caused by addition of nucleotides identical to neighbors in the target sequence (e.g. ambiguity in insertions AA and AA). We achieved this by redistributing the counts of ambiguous insertions based on the counts of unambiguous insertions and rules derived from the models (see Methods, Fig. 5b). Figure 5b shows observed counts of unambiguous insertions (filled bars) and model-corrected counts of ambiguous insertions (unfilled bars) obtained for a range of possible break sites on the non-target DNA in our cell line Cas9 assays using sgRNA1-3. Each row depicts the redistribution of ambiguous counts for each assay according to a different model, and the vertical shaded areas highlight the main break sites of interest located 3|4 (right) and 4|5 (left) nucleotides upstream of the PAM. For each plot, we also include an insert labelled ?expected?, with an illustration of the expected modelbased count distribution for the two sites 3|4 and 4|5. Blunt DNA cleavage-repair model. Based on the blunt model we expected to see template-independent insertions, and therefore similar insertion frequencies for every nucleotide at the break site (equal-height bars for sites 3|4 and 4|5 in ?expected? inserts, top row plots of Fig. 5b). When we redistributed the observed ambiguous counts in a nucleotide-unbiased manner, following the frequency of unambiguous insertions at each site, we did not obtain the expected nucleotide-unbiased insertion profiles. Specifically, redistributed counts showed significantly more insertions of G (sgRNA1), G (sgRNA2), and T (sgRNA3) at sites 3|4 and 4|5 (top row in Fig. 5b, note the log-transformed scale). We further note that the nucleotide imbalance would be unavoidable for the most targeted sites, 3|4 and 4|5, given the substantial frequency gap observed between the most inserted and the remaining nucleotides. Staggered DNA cleavage-repair model. Alternatively, the staggered model determined template-dependent insertions following a 1-nt 5? overhang on the opposite strand. According to this preference, we expected to see significantly more insertions of the nucleotide identical to the DNA base immediately downstream of the break site (?expected? inserts in middle row plots of Fig. 5b). For the staggered model, we redistributed ambiguous insertions in our observed counts based on both the frequency of unambiguous insertions, and the 1-nt 5? overhang template. The results obtained for this model showed the expected template-based insertion pattern at break site 4|5 for all sgRNAs (?expected? inserts vs. main plots in middle row of Fig. 5b). However, the results did not follow the expected pattern at break site 3|4 for sgRNA1 and sgRNA3. Specifically, results showed similar frequencies of the four nucleotides, typical of template-independent insertions in a blunt model, while nucleotides C (sgRNA1) and A (sgRNA3) should have been significantly more frequent than AsgRNA1 N R esgRNA2 d i u GsgRNA3 a b Deleted regions/nucleotides sgRNA1 sgRNA2 sgRNA3 sgRNA1 sgRNA2 sgRNA3 all 4I5 3I4 all 4I5 3I4 all 4I5 3I4 Sites neighboring top 10 deleted regions A C G T A C G T A C G T Inserted nucleotide in 1 bp insertions others based on the staggered model (main plots vs. ?expected? inserts in middle row of Fig. 5b). Combined blunt and staggered DNA cleavage-repair model. Finally, the combined blunt and staggered model provided the best fit: (i) most insertions occurred at sites 3|4 and 4|5, the expected primary targets of blunt and staggered cleavage on the non-target DNA (bottom row of Fig. 5b); (ii) the four nucleotides were similarly likely at site 3|4, consistent with templateindependent insertions at blunt DNA ends; and (iii) the most inserted nucleotide at site 4|5 matched the downstream neighbor, as expected upon fill-in of 1-nt 5? overhangs. A combination of the blunt model with an alternative staggered model inducing 1nt 3? overhangs at break site 3|4 (ntDNA) could possibly fit as well, although we found no reference to such a model in the literature. We note that our data could include re-cleavage events leading to a biased selection for mutations. Specifically, accurately ligated DNA ends can be re-cleaved, while mutations become nearly fixed in the population, since Cas9 is less likely to recognize and cleave at a mutated site. Over time, this could skew the ratio between mutated and wild-type sequences, and influence blunt vs. staggered patterns. For this reason, we limited the scope of our analysis to showing that both patterns occur, without quantifying how frequent each one is. Staggered DNA ends in Cas9-targeted mouse pre-B cells. We investigated whether 1-nt 5? overhangs would be generated by Cas9. For this purpose, we collected independent data generated by hairpin capture and sequencing of DNA end structures at Cas9 DSB sites (HCoDES)33. Specifically, we re-analyzed DNA ends sgRNA1 sgRNA2 sgRNA3 induced by Cas9 targeting to a region on chromosome 6 in mouse pre-B cells arrested in G1-phase and LigaseIV-deficient33. Many DNA end structures could not be uniquely mapped (Fig. 5c). We saw prevalent blunt ends at 3|4, although these could not be discerned from certain 2-nt 5?/3? overhangs. Specifically, blunt ends at 3|4 accounted for ~20% when aligning sequences to the reference as-is, as performed in the original study33. We also determined a proportion of ~41.6% upon masking the noise caused by incomplete bisulfite conversions. Shorter 1-nt overhangs were favored. Additionally, we saw a preference for 5? overhangs, particularly 1-nt 5? at 2|3 (tDNA) 3|4 (ntDNA) accounting for ~11.0%, and at 3|4 (tDNA) 4|5 (ntDNA) confounded with a 3-nt 5? overhang (~6.2%). Discussion As the CRISPR-Cas9 system is widely used for gene editing, understanding Cas9 activity across the genome is crucial to identify cleavage and mutation patterns enabling new applications or improvements. Using a combination of the CRISPR-Cas9 and TRIP technologies with high-throughput DNA sequencing, we characterized mutations at ~1k loci throughout the genomes of mES cells. We showed high reproducibility in a singlepromoter TRIP cell line with 36 IRs and pools of cells with thousands of heterogeneous multi-promoter TRIP IRs. In line with previous studies13, we saw sgRNA-dependent variation in Cas9-induced mutation frequency. Specifically, sgRNA3 led to larger mutation frequencies than sgRNA1/2. We reasoned that the high GC-content of sgRNA1/2 (75%, 70%) relative to sgRNA3 (45%) could explain the lower efficiency of sgRNA1/2, based on reports associating low and high GC-content with reduced Cas9 activity12. Overall, sgRNA1/2 produced comparable mutation frequencies, possibly given the high similarity including a 19-bp reverse complement overlap. Genomic location was the largest contributor to variation in mutation frequency. However, mutation frequency correlated weakly with TGE features of wildtype cells and reporter expression in TRIP cells. Importantly, reporter expression and mutation frequency were quantified in TRIP mES cells bearing identical integrations and epigenetic landscape. This result suggests that the effects of genomic context on mutation frequency may be influenced by other factors than those surveyed here. Although studies using doxycycline-controlled chromatin states on engineered targets previously showed effects on mutation efficiency19,21, an approach targeting endogenous sites in human cells revealed only modest effects with effect sizes that were largely dependent on guide RNA and possibly other factors22. We saw similar results in our study, with generally higher correlations between TGE features and mutation frequency using sgRNA3. We note that, due to the exogenous manipulation of chromatin states, those studies could produce larger contrasts between permissive and non-permissive states than those seen in the context of our work. Our findings were also consistent with two other studies assessing the association between expression or epigenome and mutation frequency at endogenous target sites17,24. One study reported weak correlations overall across developmental stages in zebrafish, which were slightly higher for expression than epigenetics24. The other study reported that native chromatin accessibility and DNA methylation were predictive of Cas9 binding, but Cas9 binding was uncorrelated with mutation frequency in mES cells17. Together, these findings suggest that genomic context likely influences Cas9 binding and cleavage, but that the association with mutation frequency can be further modulated by additional factors such as guide RNA sequence. It is possible that the correlation could be disrupted by stochasticity in the outcome of DNA repair, especially as a result of potential re-cleavage events. We also reason that dynamic reorganization of the regulatory landscape during the S-phase of the cell cycle could influence estimates of Cas9-induced mutation frequency but not IR expression, which is mainly produced in G1-phase. In particular, mES cells are known to exhibit unusually short G1 and long S-phases34. In S-phase, the DNA is unpacked to allow for replication, enabling unperturbed Cas9 binding to otherwise inaccessible loci. This could mean that most Cas9 cleavage in our assays occurred in S-phase under widespread DNA accessibility, masking differences between permissive and non-permissive domains, and hence lowering the association of mutation frequency with TGE features. Nevertheless, insertions were more dependent on the regulatory landscape than deletions, particularly concerning influencers of transcriptional activity genome-wide such as PolII complex, histone marks or Lamin-B1. We noted that deletions varied in size and pattern, while insertions comprised mostly sgRNA-specific single-nucleotides. Specifically, deletions could be generated by a variety of DNA processing events on blunt or staggered DNA ends, either aimed at ligation by NHEJ throughout the cell cycle35 or generation of 3? overhangs in preparation for HDR in S and G2 phases36. In contrast, insertions likely arose by gap-filling of specific overhangs by a DNA polymerase. We hypothesize that the deterministic nature of insertions could preserve the association with TGE features better than the range of processes involved in deletions, possibly yielding varied dependencies on the regulatory landscape. We observed small Cas9-induced lesions with all sgRNAs. The most common deletion was loss of trinucleotide CGG, accounting for ~40% of all deletions with sgRNA1. We attributed this event to microhomology-mediated end-joining (MMEJ), which could rely on the repetition of CGG in sgRNA1 to align the DNA ends resulting in the loss of triplet CGG37. Regions deleted with sgRNA3 neighbored nucleotides 4|5 (ntDNA) rather than the expected target site 3|4, seen with sgRNA1/2, suggesting alternative Cas9 targeting between 4|5 (ntDNA). In addition, singlenucleotide insertions showed high prevalence of a specific nucleotide per sgRNA, consistent with independent studies targeting the HPRT, AAVS1 and TREX1 genes in HCT116 cells38, the GFP, NDC1 and LBR genes in K562 cells39, as well as numerous other regions in HCT116, HEK293, and K562 human cells38?40. We determined that the frequently inserted DNA base matched nucleotide 4 upstream of the PAM, which hinted at eventual template-dependent repair of staggered DNA ends with termini at 3|4 and 4|5 on opposite strands. Different Cas9 cleavage and DNA repair models were analyzed, seeking to explain observed insertion patterns3,31,32. We propose a Cas9 cleavage model inducing primarily blunt and occasionally staggered DNA ends (Fig. 6). Most blunt ends are ligated, possibly upon resection, resulting in either wild-type or deletion. On occasion, blunt ends may lead to an insertion via template-independent addition of a random nucleotide (e.g. by Pol ?)9. We showed that insertions mainly derived from DNA ends with 1-nt 5? overhangs, specifically with termini at positions 3|4 (tDNA) and 4|5 (ntDNA)32. We noted that 1-nt 3? overhangs could also explain the observed insertions. However, re-analysis of independent data on DNA end structures generated by Cas9 targeting in mouse DNA Ligase IV-deficient pre-B cells33, revealed prevalence of blunt ends and preference for 1-nt 5? over 3? overhangs, in accordance with our hypothesis. We cannot exclude the possibility that the overhangs would arise by minimal resection of blunt ends rather than directly by Cas9 cleavage. Regardless of the process, we presented evidence that staggered DNA ends are generated and likely responsible for most insertions. The insertion patterns we observed in mES cells mostly in S-phase were corroborated by independent Cas9 targeting in human cells HEK293, HCT116 and K56238?40, as well as in mouse pre-B cells arrested in G1-phase33, using different plasmid systems. The consistency of Blunt DNA ends Template-independent insertion (Pol ??) Liga tion Wild-type sequence Exonuclease activity tion Liga Staggered DNA ends 3? 5? 5? 3? Template-dependent insertion (Pol ?, Pol ? ?) Ligation Ligation Ligation Random nucleotide insertion Deletion Sequence-specific insertion these results indicates that our findings might generalize beyond the scope of our assays. Here, for the first time, Cas9-induced mutation patterns were characterized at thousands of target sequences embedded throughout the genomes of mES cells. We laid out the likely combination of mechanisms of Cas9 cleavage and DNA repair underlying predictable 1-nt insertion patterns reported in the literature38?40. In particular, we revealed that Cas9 occasionally generates breaks with staggered DNA ends. These could be leveraged to increase knock-in efficiency and control the orientation of inserts into host DNA by homology dependent or independent insertion41. Furthermore, our data indicated that guide RNA sequence determines the frequency of staggered ends, ultimately influencing insertion and overall mutation efficiency. This result suggests that guide RNAs could be designed to maximize insertion frequency. Extensive testing of a wider range of guides will be needed to build models that predict insertion efficiency based on relevant guide RNA features. Our findings could have important implications for the optimization of Cas9mediated knock-in, which remains a major challenge in genome editing. In addition, we introduced the TRIP technology as a tool to multiplex RNA-guided Cas9 targeting to regions in reporter genes integrated genome-wide. This CRISPR-on-TRIP approach allowed us to seamlessly isolate target sequence from genomic location, and perform sequence-independent analysis of variation in Cas9 mutation efficiency and patterns in different genomic contexts. We found that genomic location is a key determinant of mutation frequency, which however correlated modestly with genomic and epigenomic context. We reasoned that guide RNA and stochasticity in the outcome of DNA repair, together with potential re-cleavage events and the reprogramming of the epigenetic landscape during the cell cycle could contribute to this result. Further investigation will be needed to clarify the impact of genomic context. Finally, we demonstrated that CRISPR-on-TRIP is a promising tool to profile Cas9 activity at a large number of target sequences scattered throughout the genome, and can be combined with other assays to study the influence of a variety of processes on Cas9 activity and induced mutation patterns (Fig. 7). Methods Construction of TRIP plasmid libraries. The monoclonal TRIP cell line used in this work was established in our previous TRIP study25. The piggyBac-based template vector pPTK-Gal4-mPGK-Puro-IRES-eGFP-sNRP-pA carrying the reporter unit in the TRIP cell line library comprised the following elements: piggyBac 5?-TR, 14 Gal4 binding sites, mPGK promoter, puromycin resistance (PuroR) coding sequence, encephalomyocarditis virus internal ribosome entry site (IRES), EGFP coding sequence, PstI site (used to clone barcodes) + DpnII site (used to map IRs), human soluble neuropilin-1 (sNRP-1) polyA signal, and piggyBac 3?-TR (67-bp). Barcoded inserts were generated through amplification of 5 ng template vector pPTK-Gal4-mPGK-Puro-IRES-eGFP-sNRP-pA (GenBank KC710227), using primers PB-barcode-long-7 (5?-GTGACACCTGCAGGATCA (N)16CTCGAGTTGTGGCCGGCCCTTGTGACTG-3?, where (N)16 denotes a random 16-nt long reporter barcode) and PB-barcode-short-7 (5?-GACATA ACGCGTATACTAGATTAACCCT-3?). After PCR purification, the PCR product was digested with restriction enzymes PstI and MluI (underlined). In parallel, the pPTK-Gal4-mPGK-Puro-IRES-eGFP-sNRP-pA vector was digested with the same restriction enzymes and then dephosphorylated. The digested PCR product was next ligated with the prepared vector using 10U of T4 DNA ligase (Promega). The resulting ligation product was transformed into electrocompetent E.coli cells and the plasmid DNA (TRIP plasmid library) was isolated using Genopure plasmid maxi kit (Roche). For the multi-promoter TRIP pool established in this study we first generated seven libraries, each containing reporter constructs with one of seven different promoters: CMV, cMyc, Hoxb1, Nanog, Oct4, p53 and PGK. These were included to assess effects of promoters with different characteristics, such as strength42, exogeneity, housekeeping status, retinoic acid-inducibility43?46, and TATA-less status45,47?50. The piggyBac-based plasmid vectors carrying the reporter unit comprised the following elements: piggyBac 5?-TR (314 bp long), promoter of interest, EGFP coding sequence, DpnII site (used to map IRs), 5-bp promoter index unique to the promoter, KpnI site (used to clone barcodes), sNRP-1 polyA signal and piggyBac 3?-TR (242-bp). Barcoded inserts were generated by amplification of 5 ng PB template vector using primers Kpn-RandomBC-1 (5?-AAAAGGTACC (N)18GAGTTGTGGCCGGCCCTTGTGACTG-3?, with (N)18 denoting a random 18-nt long reporter barcode) and BssH2-A (5?-AAAAGCGCGCATACTAGAT TAACCCTAGAAAGATAATCATATTG-3?). After PCR purification, the barcoded inserts were digested with restriction enzymes KpnI and BssHII (underlined). In parallel, the plasmid vectors were digested with restriction enzymes KpnI and MluI (the latter generates sticky ends compatible with those made by BssHII) and subsequently dephosphorylated. Ligation of the digested barcode insert into the digested plasmid vectors and electrotransformation of bacterial cells were performed as described above. The seven promoter-specific plasmid libraries were mixed together in the following molar ratios: 2:2:4:2:2:2:1 to obtain the multipromoter TRIP library. Cell culture and TRIP library transfection. Mouse embryonic stem (mES) cells EBRTcH3 expressing the tetracycline-controlled transactivator (tTA) from the endogenous ROSA26 promoter (EStTA)51 were cultured in 60% BRL cellconditioned medium in the presence of 10% fetal calf serum (FCS), leukemia inhibitory factor, MEK inhibitor PD0325901, and GSK-3 inhibitor CHIR9902152. The EBRTcH3 ES cells were provided by the te Riele lab (Netherlands Cancer Institute), which had received them from Dr. Masui (International Research Center of Japan). The EBRTcH3 ES cells were originally derived from E14tg2a ES cells53 by Masui and colleagues51. Culture dishes were coated with 0.15% gelatin and incubated at 37?C for one hour before plating. For the cell line, six million mES cells were plated and incubated for 4 h at 37?C. The cells were then transfected with 22.5 ?g of the mPGK TRIP plasmid library Cas9-induced mutations (frequency and patterns) CRISPR-on-TRIP assay gRNA-Cas9 targeting of TRIP IRs + DNA-seq Transcription Cas9-induced mutations Cas9-induced breaks (DNA end structures) DNA RNA DNA Assay 1: TRIP + RNA-seq Assay 2: CRISPR-on-TRIP H3K27ac H3K9me2 Assay 2: CRISPR-on-TRIP Hairpin DNA DNA RNA DNA Hairpin DNA Wild-type Deletion Insertion DNA RNA DNA DNA Blunt Staggered Hairpin DNA Insight RNA-guided Cas9 targeting of TRIP integrations genome-wide Epigenome Cas9-induced mutations Assay 2: TRIP + epigenome Chip-seq DNA DNA Assay 1: CRISPR-on-TRIP in ligase-inhibited cells Assay 2: Hairpin-based capture of DNA end structures together with 2.5 ?g of mouse codon-optimized version of PB transposase (mPB) plasmid54 using Lipofectamine 2000 (Invitrogen) and incubated at 37?C. After 48 h of incubation, the transfected cells were sorted through fluorescence-activated cell sorting (FACS), selecting single cells with ?medium? levels of EGFP expression, which were used to establish stable TRIP cell lines. The cell line EStTA-PB-B-18 used in this study carries 36 IRs. For the TRIP pool, nine million mES cells were plated and incubated for four hours at 37?C. The cells were then transfected with 15 ?g of the multi-promoter TRIP plasmid library mixture and 5 ?g of PB transposase plasmid mPB-L3-ERT2. TatRRR-mCherry plasmid using Lipofectamine 2000 (Invitrogen) and incubated at 37?C. After 24 h, the cells were FACS-sorted to select for a desired mCherry signal from the PB transposase plasmid26. The sorted cells were resuspended in mES medium with 1 ?M of tamoxifen to activate the PB transposase. After 24 h, the cells were washed and resuspended in mES medium without tamoxifen, in which they were cultured for 5 days, refreshing the medium every other day. Biological replicate pools were established by subculturing several aliquots of cells, which were grown for another week. Technical replicates were established by splitting each pool into two and growing each half separately for another week. The multipromoter TRIP pool used in this study was grown from approximately 500 cells containing an average of ~ 25 IRs per cell. Determination of genomic location and expression of IRs. Barcoded IR regions were extracted by inverse PCR followed by high-throughput DNA sequencing. The resulting reads were subject to quality control, including the filtering of aberrant barcodes arising from errors induced by PCR or sequencing. After pre-processing, the regions were aligned against the mouse genome assembly mm9 in order to map IR loci. Expression of IRs was determined by sequencing of reverse transcribed RNA (cDNA). Barcode abundance was also estimated by sequencing genomic DNA (gDNA), and used to normalize IR expression. These procedures were performed according to the TRIP protocol26. Association of TRIP integrations with genome-wide features. In addition to IR expression, we assessed endogenous gene expression25 and a range of regulatory elements measured genome-wide in mES cells. These data were previously collected from public sources and preprocessed aiming to maximize comparability28. Most features were obtained from ChIP-seq data on epigenetic modifications such as DNase I hypersensitivity, histone marks, and binding of transcriptional and epigenetic regulators. For ChIP-based features, a score was defined as the log2-transformed ratio between signal and control. Genomic features including GC content and gene proximity, among others, were derived from the mouse reference genome mm9 and Ensembl genes (release 66). Specifically, proximity measures were defined as the negative log2-transformed distance ( +1) to the nearest concerned genomic feature (e.g. gene, transcription start site). Chromatin compaction was estimated from Hi-C data as the rate of decay in contact probability ? between two loci with increasing genomic distance, locally approximated in 400 kb windows by a power-law function with scaling exponent ?. Finally, the association between every IR and feature pair was computed as the mean normalized score of the genomic or regulatory feature over a region of 2 kb surrounding the integration site. sgRNA and ssODN design. All oligonucleotides used in this study were purchased from Integrated DNA Technologies. We designed three sgRNAs targeting different regions within the EGFP reporter gene, using the CRISPR design tool http://crispr. mit.edu/ (Supplementary Table 3). We considered three primary criteria. First, purposed sgRNA-targeting of a sequence within the EGFP gene body in close proximity to the barcode, enabling reliable amplification of both barcode and target site. Second, reliability according to the CRISPR design tool, with no or minimal reported off-target sites. Third, EGFP recognition in both sense and anti-sense orientation, to exclude orientation-dependent effects. In addition, we designed a 141-bp single-stranded oligodeoxyribonucleotide (ssODN) template for HDRbased knock-in in with sgRNA2. Our ssODN comprised a 21-bp sequence for knock-in with ~ 60-bp homology arms at each side. Following the recommendations in https://www.addgene.org/crispr/zhang/faq/ (Supplementary Table 3), we avoided overlap between the ssODN homology arms and the barcode located downstream of the EGFP reporter gene. We also designed the junction of the homology arms < 10 bp from the Cas9 target site. Lastly, we designed the ssODN template with 58 and 62 nucleotide-long homology arms, in agreement with the recommended 50?80 range. Cloning and transfection of sgRNA-guided CRISPR/Cas9. We used human codon-optimized SpCas9 and chimeric guide RNA expression plasmid pX330-U6Chimeric_BB-CBh-hSpCas9 to complex Cas9 with different sgRNAs ( 1?3 ). In addition, we used mCherry as a fluorescent marker for visualization and sorting of Cas9-sgRNA transfected mES cells. All Cas9 targeting assays experiments were performed in triplicate. For cell line experiments, we co-transfected the Cas9sgRNAs and mPB-L3-ERT2.TatRRR-mCherry plasmids. Ten million mES cells were first seeded on a 10-cm dish. Four hours later, the cells were transfected with 13.5 ?g of Cas9-sgRNA and 1.5 ?g of mPB-L3-ERT2.TatRRR-mCherry plasmids using 45 ?l of Lipofectamin 2000 (Invitrogen). For TRIP pool experiments, we cloned mCherry-expressing Cas9-sgRNA plasmids. We first digested the Cas9sgRNA plasmids with 10 units NotI (Roche) and 10 units SbfI (NEB). We further PCR-amplified CMV-driven tatRRR-mCherry from mPB-L3-ERT2.TatRRRmCherry using primers Fragment.FOR (10 ?M) and Fragment.REV (10 ?M) (Supplementary Table 3). The digested Cas9-sgRNA plasmid and PCR-amplified tatRRR-mCherry fragment were then cloned using Gibson Assembly Master Mix (NEB). To avoid reduced transfection efficiency of mCherry-expressing Cas9sgRNA plasmid and ssODN with Lipofectamin 2000, we followed the protocol from Nucleofector? Kit for Mouse Embryonic Stem Cells (Lonza) instead. Five million mES cells per condition were trypsinized, spun down and resuspended in 90 ?l of Mouse ES Cell Nucleofector? solution. They were thereafter transfected with a total of 6.25 ?g of DNA (3.75 ?g mCherry-expressing Cas9-sgRNA and 2.5 ?g ssODN) using the program A-024 of the Nucleofector? Kit for Mouse Embryonic Stem Cells (Lonza). Transfected cells were resuspended in 500 ?l of prewarmed culture medium and plated in gelatin-coated 10 cm dishes. Cells were sorted 24 h after transfection. For disruption assays, we transfected five million mES cells with 8 ?g of mCherry-expressing Cas9-sgRNA. For editing, five million cells were transfected with 4 ?g of mCherry-expressing Cas9-sgRNA together with 1 ?g of designed ssODN containing the 21-bp sequence for knock-in with ~ 60-bp homology arms. We decided to use sgRNA2 in editing assays, given that the short distance between the region targeted by sgRNA and the IR barcode prevented the design of proper ssODN homology arms (Supplementary Figure 8). Additionally, sgRNA2 was preferred over sgRNA1 due to its larger mutation efficiency. In all experiments, transfected mCherry-positive mES cells were sorted by flow cytometry (MoFlo), collected in conditioned media containing 20% FCS, and subsequently spun down plated in conditioned media containing 10% FCS. Finally, Cas9-sgRNA-transfected mES cells were expanded for five days before isolating genomic DNA. DNA isolation and preparation of samples for sequencing. After incubation of Cas9-sgRNA-transfected mES cells, DNA was isolated using DNeasy Blood & Tissue Kit (Qiagen) and prepared for sequencing. We PCR amplified cell line DNA IR regions of 414 bp surrounding the sgRNA target sites using 10 ?M of PB-cDNAforward-1-BC primers, containing different index sequences for multiple reactions, and 10 ?M of PB-cDNA-Reverse-5 primer (Supplementary Table 4). The PCR product was sent for Illumina MiSeq. In TRIP pool experiments, we tagged amplicons using 16-nucleotide Unique IDentifiers (UIDs) to be able to detect multiple readings of the same DNA molecules (Supplementary Table 5)55. We performed all Cas9 targeting assays in triplicate. Barcoded IR regions were PCR amplified using Phusion? High-Fidelity DNA Polymerase (NEB) following the protocol in Supplementary Table 6. We used Exonuclease-I (20 units) (Enzymatics) for cleavage of single stranded DNA. Purified sequences were sent for Illumina HiSeq. Replicates were assessed for consistency of IR mutation frequencies, and merged for subsequent analyses. Validation of promoter-barcode index association. For identification purposes, each promoter was originally associated with a unique 5-nucleotide index located 44 base-pairs downstream of the IR (Supplementary Table 2). We performed Sanger sequencing on the DNA extracted for five randomly selected IRs in the TRIP pool to confirm the absence of recombination and the correct association between barcodes and corresponding promoters. We designed primers for nested PCR, including two IR specific primers, PB-Valid.3-Out primer-1 and Inner-1, and locus-specific genome spanning primers, GEMP_mvla_outer and GEMP_mval_inner (Supplementary Table 7). PCR amplification was performed in two steps using the primer combinations PB-Valid.3-Out primer-1/GEMP_mval_outer for PCR1 and Inner-1/GEMP_mval_inner for PCR2. The resulting products were diluted in water (1:5) and prepared for Sanger sequencing using 10 ?M of primers together with Big Dye terminator version 3.1. Identification and characterization of Cas9-induced lesions. We parsed the DNA reads obtained from Cas9-targeted TRIP cells in order to map barcodes and UIDs25, and extract the sequences of interest containing the 20-nucleotide region targeted by Cas9. These were aligned against the wild-type sequence to identify and characterize Cas9-induced lesions. We used semi-global alignment with the following weights: match + 2, mismatch ?2, gap opening penalty ?5, gap extension penalty ?0.5, and initial score 30. For TRIP cell line experiments with sgRNA1, sgRNA2 and sgRNA3, respectively 96%, 95% and 91% of the reads could be parsed, whereas TRIP pools yielded between 86% and 89% of successfully parsed reads. Quality control and filtering of aberrant barcodes were further performed as described in previous work25. Analysis of mutation frequencies, sizes and patterns. We analyzed sequenced IR regions with a read coverage of at least 30 in all Cas9 assays on TRIP cells. The average read coverage per IR was significantly higher in the TRIP cell line than the pool (Supplementary Figure 9). However, the cell line contained only 36 barcoded regions whereas the pool offered ~ 1k IRs. As a result, we used TRIP pool data primarily to assess overall trends of variation in IR mutation frequency throughout the genome, and high-resolution TRIP cell line data to identify mutation sizes and patterns. Relative frequencies of each lesion type were calculated per IR as the ratio between the number of reads exhibiting such lesion and all reads for the given IR (Figs. 2a?c). We plotted frequencies per lesion type for each of 36 IRs in the TRIP cell line (Fig. 2a), as well as distributions (Fig. 2b) and stacked medians stratified by promoter (Fig. 2c) considering all 1359 IRs shared by TRIP pool assays. To assess variation in IR mutation frequency with genomic locus in the TRIP pool, we fitted linear regression models describing the relationship between mutation frequencies in each pair of assays and determined corresponding R2 and F-statistic p-values (Fig. 2d). Additionally, we fitted a global linear regression model to IR mutation frequencies in all TRIP pool assays taking into account genomic locus, guide RNA, knock-in, promoter and a significant interaction term between genomic locus and guide RNA (Fig. 2f). Based on this model, we determined ?-squared values denoting the effect size or variance explained by each independent variable (Fig. 2g). For this, we relied on multi-way ANOVA tests with partial sum of squares preserving the principle of marginality (type II sum of squares). The ratio between knock-in and other insertions was calculated per IR and plotted against binned total IR mutation frequency together with the mean and confidence interval of 0.95 obtained by nonparametric bootstrapping (Fig. 2d). Association between IR total mutation frequency and regulatory element scores was determined using Pearson?s correlation, along with the corresponding t-test p-values corrected for multiple testing using the Benjamini-Hochberg procedure (Figs. 3b?d). We calculated and plotted overall deletion and insertion frequencies per size in the TRIP cell line (Fig. 4a). In addition, we determined the frequencies of individual deletion patterns (Fig. 4b) and the distributions of insertions of the four DNA bases in predominant one-nucleotide insertions (Fig. 4d). Deletions with ambiguous positions could not be distinguished. We took this ambiguity into account when calculating the frequency of deletions neighboring the expected or alternative break site. Specifically, we weighed down the frequency according to the ratio of ambiguous positions neighboring the corresponding sites (Fig. 4c). All differences between distributions reported in the manuscript were determined using nonparametric two-tailed Wilcoxon tests. All given confidence intervals are calculated for confidence level 0.95. Analysis of one-nucleotide insertion patterns. In order to assess different models3,31,32 of Cas9 cleavage, we had to first resolve ambiguity in insertion sites. Ambiguity arises due to insertions of DNA bases next to matching nucleotides in the target sequence. For instance, insertions of an A immediately up or downstream of a wild-type A cannot be distinguished, as they result in identical mutated sequences (AA and AA). Since these results are confounded, both positions represent ambiguous insertion sites for the added base A. These observations are also true for longer and more complex sequences. Additionally, the deterministic nature of variant callers means that an insertion is called where the change is identified, thus once the repetition is detected. When processing sequences 5? to 3?, insertions above would always be identified as occurring downstream of the wildtype A (AA). We redistributed accumulated insertion counts throughout ambiguous sites based on unambiguous insertions and the Cas9 cleavage and DNA repair models as follows (Fig. 5b). The wild-type non-target sequence was processed to identify ambiguous DNA base(s) at each site, essentially those matching the 5? and 3? neighbors in the wild-type sequence. Individual ambiguous sites were extended to maximal windows of consecutive ambiguous sites for the same DNA base. We performed the redistribution of ambiguous DNA base insertions within these windows. For blunt ends with template-independent insertions, without bias towards specific DNA bases, ambiguous counts were distributed to follow the relative proportions of unambiguous insertions at the sites within the window. For staggered ends repaired by template-based insertions, we classified sites within a window as likely or unlikely based on the expected insertion pattern. Considering a 1-bp 5? overhang model, ambiguous sites with downstream nucleotide identical to the inserted DNA base would be considered likely whereas those with a different downstream nucleotide would be unlikely given the conformation of the DNA ends. We re-calculated the counts for ambiguous unlikely sites as the average of the counts observed for unambiguous DNA bases at those same sites. Likely sites received the remaining counts, distributed according to the relative proportions of unambiguous DNA bases. For the combined blunt and staggered models targeting primarily between nucleotides 3?4 and 4?5 upstream of the PAM (non-target DNA), we applied the above strategies for blunt and staggered models respectively to these sites. Since the ratios between blunt and staggered were not known for the remaining sites, we used the simpler estimate based on relative proportions of unambiguous DNA base insertions. Finally, we compared the estimates against expected distributions (Figs. 5a, b). Re-analysis of Cas9-induced DNA end structure data. We downloaded targeted sequencing data of DNA ends induced by Cas9 targeting using the SRA-toolkit, after selection using the NCBI SRA Run Selector, specifically runs SRR1617071 to SRR1617082 in project PRJNA264361. The assays involved Cas9 cleavage guided to a region on chromosome 6 in mouse Ligase IV-deficient pre-B cells arrested in G1 phase33. Upon Cas9 targeting, top and bottom DNA ends were ligated into a hairpin and subject to bisulfite treatment to facilitate amplification for highthroughput DNA sequencing. We aligned the sequences in each file against the wild-type hairpin sequence assuming blunt ends between nucleotides 3?4 upstream of the PAM. For this purpose, we used BLASTn with the following options -max_target_seqs 500000 -max_hsps 1 -num_threads 30 -outfmt ?6 qacc sacc qstart qend sstart send evalue bitscore length pident nident mismatch positive gaps ppos btop?. We filtered low quality sequences and primers by keeping alignments at least 50-bp long containing at most 5 gaps. We also filtered sequences whose alignment did not span a region of at least 50-bp around the expected target site. We then used both the starting position of each alignment in the reference sequence and the BLASTn trace-back operations (BTOP) to isolate the pattern of changes relative to wild-type within a region of 30-bp around the expected target site. In parallel, we generated mutation patterns for all possible DNA end conformations within a target region of 30-bp around the expected target site. Finally, we matched the generated mutation patterns against the patterns extracted from the HCoDES data to identify the corresponding DNA end structures. Note that many of the HCoDES sequences contained substitutions at random positions, primarily TC, possibly also due to failed bisulfite conversion. Since substitutions did not show consistent patterns associated with particular DNA end conformations and had no impact in the ranking of top observed patterns, we decided to ignore them so that they would not confound otherwise valid matching patterns. Reporting Summary. Further information on experimental design is available in the Nature Research Reporting Summary linked to this article. Data availability High-throughput sequence data and processed data that support the findings of this study have been deposited at the NCBI GEO with accession GSE127752. The mapping and expression reads for the cell line containing 36 TRIP IRs correspond to samples GSM1207441, GSM1207452, GSM1207463, GSM1207474 in GEO series GSE48606. These are directly linked in our record GSE127752. The HCoDES sequence data were retrieved from the NCBI SRA repository as samples SRR1617071 to SRR1617082 from SRA study SRP049086. Source data for the figures in this manuscript and in the Supplementary Information file are available in the Source Data file. All other relevant data are available from the authors upon reasonable request. Acknowledgements We thank the NKI Genomics Core Facility and the NKI Flow Cytometry Facility for their technical support. We also thank Anton Berns, Heinz Jacobs, Hein te Riele, Tim Harmsen, Eva Brinkman and Bas van Steensel for helpful discussions and/or critical reading of the manuscript. This research was supported by a Netherlands Organization for Scientific Research (NWO) grant 823.02.007 to W.A. and M.v.L. A.V.P. was supported by an European Research Council (ERC) Advanced Grant 293662 to Bas van Steensel. Author contributions S.G., J.P.G., W.A., A.V.P., L.F.A.W. and M.v.L. designed the research. S.G., W.A. and A. V.P. performed the experimental work. W.A. and J.d.J. preprocessed the TRIP and epigenomic data. J.P.G. designed and performed the subsequent data analysis, and made most figures. S.G., J.P.G. and W.A. interpreted the results. S.G. and J.P.G. wrote the manuscript. L.F.A.W. and M.v.L. revised the manuscript, and supervised the analytical and experimental research. Additional information Supplementary Information accompanies this paper at https://doi.org/10.1038/s41467019-09551-w. Competing interests: The authors declare no competing interests. Reprints and permission information is available online at http://npg.nature.com/ reprintsandpermissions/ Journal peer review information: Nature Communications thanks the anonymous reviewers for their contribution to the peer review of this work. Peer reviewer reports are available. Publisher?s note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations. Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as 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 images or other third party material in this article are included in the article?s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article?s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/ licenses/by/4.0/. 1. Mali , P. et al. RNA-guided human genome engineering via Cas9 . Science 339 , 823 - 826 ( 2013 ). 2. Cong , L. et al. Multiplex genome engineering using CRISPR/Cas systems . Science 339 , 819 - 823 ( 2013 ). 3. Jinek , M. et al. A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity . Science 337 , 816 - 821 ( 2012 ). 4. Jiang , F. , Zhou , K. , Ma , L. , Gressel , S. & Doudna , J. A. A Cas9-guide RNA complex preorganized for target DNA recognition . Science 348 , 1477 - 1481 ( 2015 ). 5. Jinek , M. et al. Structures of Cas9 endonucleases reveal RNA-mediated conformational activation . Science 343 , 1247997 ( 2014 ). 6. Sternberg , S. H. , Redding , S. , Jinek , M. , Greene , E. C. & Doudna , J. A. DNA interrogation by the CRISPR RNA-guided endonuclease Cas9 . Nature 507 , 62 - 67 ( 2014 ). 7. Anders , C. , Niewoehner , O. , Duerst , A. & Jinek , M. Structural basis of PAMdependent target DNA recognition by the Cas9 endonuclease . Nature 513 , 569 - 573 ( 2014 ). 8. Jinek , M. et al. RNA-programmed genome editing in human cells . eLife 2 , e00471 ( 2013 ). 9. Lieber , M. R. The mechanism of double-strand DNA break repair by the nonhomologous DNA end-joining pathway . Annu. Rev. Biochem . 79 , 181 - 211 ( 2010 ). 10. Doench , J. G. et al. Rational design of highly active sgRNAs for CRISPR-Cas9- mediated gene inactivation . Nat. Biotechnol . 32 , 1262 - 1267 ( 2014 ). 11. Doench , J. G. et al. Optimized sgRNA design to maximize activity and minimize off-target effects of CRISPR-Cas9 . Nat. Biotechnol . 34 , 184 - 191 ( 2016 ). 12. Wang , T. , Wei , J. J. , Sabatini , D. M. & Lander , E. S. Genetic screens in human cells using the CRISPR-Cas9 system . Science 343 , 80 - 84 ( 2014 ). 13. Xu , H. et al. Sequence determinants of improved CRISPR sgRNA design . Genome Res . 25 , 1147 - 1157 ( 2015 ). 14. Chari , R. , Mali , P. , Moosburner , M. & Church , G. M. Unraveling CRISPRCas9 genome engineering parameters via a library-on-library approach . Nat. Methods 12 , 823 - 826 ( 2015 ). 15. Hsu , P. D. , Lander , E. S. & Zhang , F. Development and applications of CRISPR-Cas9 for genome engineering . Cell 157 , 1262 - 1278 ( 2014 ). 16. Kuscu , C. , Arslan , S. , Singh , R. , Thorpe , J. & Adli , M. Genome-wide analysis reveals characteristics of off-target sites bound by the Cas9 endonuclease . Nat. Biotechnol . 32 , 677 - 683 ( 2014 ). 17. Wu , X. et al. Genome-wide binding of the CRISPR endonuclease Cas9 in mammalian cells . Nat. Biotechnol . 32 , 670 - 676 ( 2014 ). 18. Knight , S. C. et al. Dynamics of CRISPR-Cas9 genome interrogation in living cells . Science 350 , 823 - 826 ( 2015 ). 19. Horlbeck , M. A . et al. Nucleosomes impede Cas9 access to DNA in vivo and in vitro. Elife 5 , e12677 ( 2016 ). 20. Yarrington , R. M. , Verma , S. , Schwartz , S. , Trautman , J. K. & Carroll , D. Nucleosomes inhibit target cleavage by CRISPR-Cas9 in vivo . Proc. Natl Acad. Sci. USA 115 , 9351 - 9358 ( 2018 ). 21. Chen , X. et al. Probing the impact of chromatin conformation on genome editing tools . Nucleic Acids Res . 44 , 6482 - 6492 ( 2016 ). 22. Daer , R. M. , Cutts , J. P. , Brafman , D. A. & Haynes , K. A. The impact of chromatin dynamics on Cas9-mediated genome editing in human cells . ACS Synth. Biol . 6 , 428 - 438 ( 2017 ). 23. Jensen , K. T. et al. Chromatin accessibility and guide sequence secondary structure affect CRISPR-Cas9 gene editing efficiency . FEBS Lett . 591 , 1892 - 1901 ( 2017 ). 24. Uusi-M?kel? , M. I. E. et al. Chromatin accessibility is associated with CRISPRCas9 efficiency in the zebrafish (Danio rerio) . PLoS ONE 13 , e0196238 ( 2018 ). 25. Akhtar , W. et al. Chromatin position effects assayed by thousands of reporters integrated in parallel . Cell 154 , 914 - 927 ( 2013 ). 26. Akhtar , W. et al. Using TRIP for genome-wide position effect analysis in cultured cells . Nat. Protoc . 9 , 1255 - 1281 ( 2014 ). 27. Hsu , P. D. et al. DNA targeting specificity of RNA-guided Cas9 nucleases . Nat. Biotechnol . 31 , 827 - 832 ( 2013 ). 28. de Jong, J. et al. Chromatin landscapes of retroviral and transposon integration profiles . PLoS Genet . 10 , e1004250 ( 2014 ). 29. Liang , G. et al. Distinct localization of histone H3 acetylation and H3-K4 methylation to the transcription start sites in the human genome . Proc. Natl Acad. Sci. USA 101 , 7357 - 7362 ( 2004 ). 30. Li , Q. et al. Polycomb CBX7 directly controls trimethylation of histone H3 at lysine 9 at the p16 locus . PLoS. One . 5 , e13732 ( 2010 ). 31. Gasiunas , G. , Barrangou , R. , Horvath , P. & Siksnys , V. Cas9-crRNA ribonucleoprotein complex mediates specific DNA cleavage for adaptive immunity in bacteria . Proc. Natl Acad. Sci. USA 109 , E2579 - E2586 ( 2012 ). 32. Zuo , Z. & Liu , J. Cas9-catalyzed DNA cleavage generates staggered ends: evidence from molecular dynamics simulations . Sci. Rep . 5 , 37584 ( 2016 ). 33. Dorsett , Y. et al. HCoDES reveals chromosomal DNA end structures with single-nucleotide resolution . Mol. Cell 56 , 808 - 818 ( 2014 ). 34. Stead , E. et al. Pluripotent cell division cycles are driven by ectopic Cdk2, cyclin A/E and E2F activities . Oncogene 21 , 8320 - 8333 ( 2002 ). 35. Biehs , R. et al. DNA Double-Strand Break Resection Occurs during Nonhomologous End Joining in G1 but Is Distinct from Resection during Homologous Recombination . Mol. Cell 65 , 671 - 684 e675 ( 2017 ). 36. Symington , L. S. End resection at double-strand breaks: mechanism and regulation . Cold Spring Harb. Perspect Biol . 6 , a016436 ( 2014 ). 37. McVey , M. & Lee , S. E. MMEJ repair of double-strand breaks (director's cut): deleted sequences and alternative endings . Trends Genet . 24 , 529 - 538 ( 2008 ). 38. Liao , S. , Tammaro , M. & Yan , H. Enriching CRISPR-Cas9 targeted cells by cotargeting the HPRT gene . Nucleic Acids Res . 43 , e134 ( 2015 ). 39. Brinkman , E. K. , Chen , T. , Amendola , M. & van Steensel , B. Easy quantitative assessment of genome editing by sequence trace decomposition . Nucleic Acids Res . 42 , e168 ( 2014 ). 40. van Overbeek, M. et al. DNA repair profiling reveals nonrandom outcomes at Cas9-mediated breaks . Mol. Cell 63 , 633 - 646 ( 2016 ). 41. He , X. et al. Knock-in of large reporter genes in human cells via CRISPRCas9-induced homology-dependent and independent DNA repair . Nucleic Acids Res . 44 , e85 ( 2016 ). 42. Boshart , M. et al. A very strong enhancer is located upstream of an immediate early gene of human cytomegalovirus . Cell 41 , 521 - 530 ( 1985 ). 43. Studer , M. , Popperl , H. , Marshall , H. , Kuroiwa , A. & Krumlauf , R. Role of a conserved retinoic acid response element in rhombomere restriction of Hoxb-1 . Science 265 , 1728 - 1732 ( 1994 ). 44. Ogura , T. & Evans , R. M. A retinoic acid-triggered cascade of HOXB1 gene activation . Proc. Natl Acad. Sci. USA 92 , 387 - 391 ( 1995 ). 45. Okazawa , H. et al. The oct3 gene, a gene for an embryonic transcription factor, is controlled by a retinoic acid repressible enhancer . EMBO J . 10 , 2997 - 3005 ( 1991 ). 46. Kuroda , T. et al. Octamer and Sox elements are required for transcriptional cis regulation of Nanog gene expression . Mol. Cell. Biol . 25 , 2475 - 2485 ( 2005 ). 47. Tamaru , M. et al. Selective activation of testis-specific genes in cultured rat spermatogenic cells . Biochim. Biophys. Acta 1049 , 331 - 338 ( 1990 ). 48. McBurney , M. W. et al. The mouse Pgk-1 gene promoter contains an upstream activator sequence . Nucleic Acids Res . 19 , 5755 - 5761 ( 1991 ) 49. Wu , D. Y. & Yao , Z. Isolation and characterization of the murine Nanog gene promoter . Cell Res . 15 , 317 - 324 ( 2005 ). 50. Bernard , O. , Cory , S. , Gerondakis , S. , Webb , E. & Adams , J. M. Sequence of the murine and human cellular myc oncogenes and two modes of myc transcription resulting from chromosome translocation in B lymphoid tumours . EMBO J . 2 , 2375 - 2383 ( 1983 ). 51. Masui , S. et al. An efficient system to establish multiple embryonic stem cell lines carrying an inducible expression unit . Nucleic Acids Res . 33 , e43 ( 2005 ). 52. Ying , Q. L. et al. The ground state of embryonic stem cell selfrenewal . Nature 453 , 519 - 523 ( 2008 ). 53. Hooper , M. , Hardy , K. , Handyside , A. , Hunter , S. & Monk , M. HPRTdeficient (Lesch-Nyhan) mouse embryos derived from germline colonization by cultured cells . Nature 326 , 292 - 295 ( 1987 ). 54. Cadi?anos , J. & Bradley , A. Generation of an inducible and optimized piggyBac transposon system . Nucleic Acids Res . 35 , e87 ( 2007 ). 55. Kinde , I. , Wu , J. , Papadopoulos , N. , Kinzler , K. W. & Vogelstein , B. Detection and quantification of rare mutations with massively parallel sequencing . Proc. Natl Acad. Sci. USA 108 , 9530 - 9535 ( 2011 ).


This is a preview of a remote PDF: https://www.nature.com/articles/s41467-019-09551-w.pdf

Santiago Gisler, Joana P. Gonçalves, Waseem Akhtar, Johann de Jong, Alexey V. Pindyurin, Lodewyk F. A. Wessels, Maarten van Lohuizen. Multiplexed Cas9 targeting reveals genomic location effects and gRNA-based staggered breaks influencing mutation efficiency, Nature Communications, 2019, DOI: 10.1038/s41467-019-09551-w