- Methodology article
- Open Access
UDiTaS™, a genome editing detection method for indels and genome rearrangements
BMC Genomics volume 19, Article number: 212 (2018)
Understanding the diversity of repair outcomes after introducing a genomic cut is essential for realizing the therapeutic potential of genomic editing technologies. Targeted PCR amplification combined with Next Generation Sequencing (NGS) or enzymatic digestion, while broadly used in the genome editing field, has critical limitations for detecting and quantifying structural variants such as large deletions (greater than approximately 100 base pairs), inversions, and translocations.
To overcome these limitations, we have developed a Uni-Directional Targeted Sequencing methodology, UDiTaS, that is quantitative, removes biases associated with variable-length PCR amplification, and can measure structural changes in addition to small insertion and deletion events (indels), all in a single reaction. We have applied UDiTaS to a variety of samples, including those treated with a clinically relevant pair of S. aureus Cas9 single guide RNAs (sgRNAs) targeting CEP290, and a pair of S. pyogenes Cas9 sgRNAs at T-cell relevant loci. In both cases, we have simultaneously measured small and large edits, including inversions and translocations, exemplifying UDiTaS as a valuable tool for the analysis of genome editing outcomes.
UDiTaS is a robust and streamlined sequencing method useful for measuring small indels as well as structural rearrangements, like translocations, in a single reaction. UDiTaS is especially useful for pre-clinical and clinical application of gene editing to measure on- and off-target editing, large and small.
Common assays for genome editing involve PCR amplification of a targeted genomic region and subsequent analysis, either by endonuclease cleavage at base mismatches [1, 2] or sequencing [3,4,5]. However, PCR-mediated assays are fundamentally unable to measure structural changes to the genome in conjunction with small indels. Unintended translocations and other structural changes have been specifically called out for investigation in genome editing therapies by the NIH recombinant DNA advisory Committee (RAC)  and the FDA . Measuring structural changes has recently become more feasible using a method called AMP-Seq (Anchored Multiplex PCR sequencing)  that is intended for clinical detection of oncogenic translocations and a similar method called HTGTS (High throughput, genome-wide translocation sequencing)  or LAM-HTGTS (linear amplification-mediated-HTGTS) . GUIDE-seq , a modification of AMP-seq, is a powerful tool to capture de novo off-target editing by CRISPR RNA-guided nucleases. All these methods utilize a target specific primer in addition to an adapter ligated universal priming site on sheared DNA to achieve “uni-directional” amplification, and sequencing. However, DNA shearing is a cumbersome step in the library preparation process for all these methods: it requires specialized equipment and it is not readily amenable to studies with large numbers of samples. Consequently, it has not been broadly applied in the gene editing field. In addition, shearing DNA has been shown to induce DNA damage that results in base miscalling [12, 13], something potentially problematic when trying to assess gene editing frequencies at low levels (e.g.: less than 1%).
Tagmentation has emerged as a straightforward tool that simultaneously fragments and adds adapters in a fast (~ 10 min) enzymatic reaction [14,15,16]. In UDiTaS, a custom-designed Tn5 transposon has been developed, which contains the full-length Illumina forward (i5) adapter, a sample barcode, and unique molecule identifier (UMI) (Fig. 1a). The full UDiTaS process (Fig. 1b) is comprised of the tagmentation step, followed by two PCR steps: round one PCR uses the target specific anchor primer, and round two adds the reverse (i7) Illumina adapter and an additional sample barcode. UDiTaS has been tested with multiple primers and genomic DNA samples and is amenable to all labs using genome editing technologies with access to NGS equipment.
As a case study for assessing complex gene editing events, a pair of S. aureus Cas9 (SaCas9) sgRNA were used to target a region within intron 26 of the CEP290 gene. Reduction of CEP290 expression leads to several human diseases including a blindness condition termed Leber’s Congenital Amaurosis Type 10 (LCA10) . The most frequent cause of LCA10 is a rare single nucleotide variant in intron 26 that creates a splice donor, leading to an additional exon, and prematurely truncating the protein . Removal of the deleterious intron 26 splice donor site is predicted to restore CEP290 gene function. A single guide RNA pair, CEP290-64 and CEP290-323, that cuts 1176 base pairs around the splice donor has been selected from an internal screen (not shown) for further characterization.
To illustrate the utility of UDiTaS, U-2 OS cells were transfected with linear DNA fragments expressing sgRNA CEP290-64 and CEP290-323 along with a plasmid expressing SaCas9; after three days genomic DNA was isolated, UDiTaS libraries created, sequenced on a MiSeq, and the data processed through a bioinformatics pipeline (Additional file 1 Figure S1 and Additional file 2). Using a targeted primer flanking the guide 323 cut site (Fig. 1c) a range of edits and rearrangements around the expected cut sites were observed (Additional file 3: Figure S2a-d), automatically classified, and tallied (Fig. 2a). Editing that resulted in small indel events were observed, as expected, at a rate of ~ 17%. In addition, junctions from the desired ~ 1.1 kb deletion were also present at ~ 40%. Notably, inversions of the ~ 1.1 kb fragment between the two cut sites were also observed at ~ 18%, comparable to the deletion. Other lower frequency junctions were also observed at ~ 0.75%, including translocations between homologous or sister chromosomes at the identical cut sites (Additional file 3: Figure S2 d).
To characterize the method’s linearity, accuracy, and Lower Limit of Detection (LLoD) for the desired large editing events, we constructed a stable cell line and plasmids that contained the CEP290 intron 26 wild type locus, the deletion, and the inversion. Our HEK293 cells have three copies of the CEP290 locus and the clone we created has two deletions and one inversion. Genomic DNA (gDNA) from this edited HEK293 cell line was mixed with parental, non-edited gDNA to generate a titration of the ~ 1.1 kb deletion and inversion across a five-log range. Plotting expected versus UDiTaS measured editing rates for the deletion and inversion showed excellent correlation down to approximately 0.1% (Fig. 2b). When comparing experiments on the same samples we observed high reproducibility, R2 = 0.99 (Additional file 4: Figure S3).
The UDiTaS protocol uses 50 ng of input DNA, which is equivalent to approximately 14,300 human haplomes. Assuming a binomial sampling distribution and ~ 20% process yield, one would expect, with 95% confidence, 2-3 observations at 0.1% in this UDiTaS library; consistent with our observations (Additional file 5: Figure S4). Increasing the sensitivity is theoretically possible by increasing the input DNA along with sequencing read depth; as both are needed to increase the number of unique UMIs in the analysis. On-target mapping rates to the genome for this study were > 95%, also indicative of the process being robust and productive (Additional file 6: Figure S5).
To further demonstrate the linearity of UDiTaS and to compare it to AMP-Seq, reference plasmids synthesized to contain the intron 26 wild type locus, large deletion, and inversion, were used as samples. These plasmids ranging from ~ 2200 to ~ 714,000 genome equivalents, were spiked into mouse genomic DNA and processed through the UDiTaS and Amp-Seq methods. UDiTaS showed excellent linearity down to the lowest dilution (Additional file 7: Figure S6a). As expected, UDiTaS and Amp-Seq yielded similar results. However, UDiTaS libraries were more linear with higher complexity given the same DNA input material (Additional file 7: Figure S6a). We attribute this to the more efficient tagmentation process compared to shearing and adapter ligation, and fewer, more streamlined processing steps that increased the overall yield. Similar results were obtained with UDiTaS when DNA was diluted without carrier mouse genomic DNA (Additional file 8: Figure S7a) further demonstrating the robustness of the process.
UDiTaS can also measure inter-chromosomal translocation rates. To demonstrate this, we simultaneously nucleofected two S. pyogenes CAS9 (SpCas9) ribonucleoprotein (RNP) complexes with sgRNAs targeting TRAC and B2M genes into activated human CD4+ T-cells. UDiTaS libraries from these cells were prepared with primers flanking each guide site. Of ten possible end joining outcomes, seven were theoretically possible to measure with the two targeting primers used in the study; and all seven of those events were detected and quantified in the experiment. These included typical NHEJ (non-homologous end joining) editing at the cut sites as well as balanced, acentric, and dicentric fusions between homologous or sister chromosomes as well as between distinct unrelated chromosomes (Fig. 2c). Although several studies have been published characterizing translocations in the context of gene editing [10, 19,20,21], none have holistically measured all events, relying on translocation assays that need to be pieced together and do not contextualize smaller indels rates. To our knowledge this is the first comprehensive study quantifying double gene editing translocation rates with indels. This study demonstrates inter-chromosomal translocation rates of ~ 2.5% with on-target indel editing of ~ 82% and ~ 91% (Fig. 2c) respectively for each guide RNA.
Linearity and an LLoD of ~ 0.01% of the assay at TRAC and B2M translocation loci were characterized with plasmids in a similar fashion as for the CEP290 locus (Additional file 7: Figure S6b, c and Additional file 8: Figure S7b). UDiTaS, when compared to AMP-Seq, was significantly more linear for both primers and the NGS libraries had higher complexity (unique molecules per read) further demonstrating the high sensitivity of UDiTaS to detect translocation events. Finally, because there were high single guide editing rates, we were able to cross-compare UDiTaS with methods commonly used for indel measurement and found excellent concordance with UDiTaS (Additional file 9: Figure S8).
We have developed a sequencing and analysis methodology that enables simultaneous measurement of small indels and larger structural rearrangements such as large deletions, inversions, and translocations. The UDiTaS method is robust, scalable, and available to any lab practicing genome editing with access to Next Generation Sequencing. Of note, this methodology may be useful in other gene editing settings. For example, multiplexing the anchor primers along with low (50 ng) input DNA will enable panels of candidate off-target editing sites to be monitored when samples are limiting. In addition, the custom transposon described here has potential to improve methods that utilize DNA shearing along with an anchor primer, such as GUIDE-Seq. We show excellent accuracy and linearity of the method for several primers. With that said, a potential limitation of the method is measurement inaccuracy that can emerge from biases of Tn5* transposition, primer binding, or amplification. As UDiTaS assays are developed, especially those with clinical implications, it is and will be important to validate and calibrate the assays with standards, as shown here using plasmids, to ensure accuracy.
UDiTaS is an important new sequencing methodology for genome editing detection and analysis. It enables the accurate quantification of intended or unintended large structural changes in addition to small, more typical indels and SNVs that arise from single and dual gene edits. Detection methods, like UDiTaS, are especially important in therapeutic gene editing settings, where editing needs to be carefully monitored to assess efficacy and therapeutic risk.
UDiTaS and NGS methods
Detailed protocols for AMP-seq and UDiTaS are provided as Additional file 2. Modifications to the UDiTaS protocol are described in the Plasmid Sensitivity and T-cell TRAC/B2M experiments. Plasmids and Oligos are listed in Table 1 and Table 2.
The analysis pipeline was built using python code that calls additional software for specialized steps (see Supplemental Fig. 1). Code is available at https://github.com/editasmedicine/uditas. Briefly, it consists of the following steps:
Demultiplexing. Sequencing reads were first demultiplexed into the different experiments in the run using the appropriate sequencing barcodes, allowing up to one mismatch in each barcode. UMIs for each read were extracted for further downstream analysis.
Trimming. 3′ adapters were trimmed using cutadapt , version 1.9.1
Create reference amplicons. We used the expected cut sites to build reference amplicons with the expected chromosomal rearrangements: wild type, large deletion, inversion, translocation, etc.
Alignment. Paired reads were then globally aligned to all the reference amplicons using bowtie2 , version 2.1.0. Finally, samtools  (version 1.3-5-g664cc5f) was used to create and index sorted bam files.
Alignment analysis. Reads completely covering a window around the predicted junctions (15 bp) were extracted and the total number of unique UMIs counted.
Final genome wide analysis. Finally, reads that could not be mapped to the reference amplicons were extracted and mapped globally using bowtie2 to the appropriate background reference genome.
U-2 OS bulk editing transfection experiment
U-2 OS cells (ATCC) were maintained in DMEM, high glucose with Glutamax and sodium pyruvate (ThermoFisher), 10% Fetal Bovine Serum, and supplemented with 1% penicillin/streptomycin. Cells were transfected by Lonza nucleofection using the 4D nucleofector system. Briefly, 250,000 cells were transfected with 1.5μg plasmid pAF003 expressing SaCas9 driven by CMV promoter and 500 ng of linear DNA fragment expressing gRNAs driven by U6 promoter (250 ng each guide). Cells were nucleofected using the SE kit and pulse code DN-100 and plated in 6-well plates. Cells were cultured for 3 days post-nucleofection and 3 transfection technical replicates were pooled together. Genomic DNA was isolated using the Agencourt DNAdvance kit (Beckman Coulter) according to manufacturer’s instructions.
HEK293 cell line creation at the CEP290 locus
Hek293 cells (ATCC) were maintained in DMEM, high glucose with Glutamax and sodium pyruvate (ThermoFisher), 10% Fetal Bovine Serum, and supplemented with 1% penicillin/streptomycin. Cells were transfected using the Mirus TransIT 293 kit, according to manufacturer’s instructions. Briefly, 120,000 cells were seeded in a well of a 24-well plate 24 h pre-transfection. Cells were transfected with 750 ng plasmid pAF003 expressing SaCas9 driven by CMV promoter and 250 ng of linear DNA fragment expressing gRNAs driven by U6 promoter (125 ng each guide). Following expansion, cells were trypsinized, diluted and re-plated in 96-well plates at a dilution of approximately 1 cell per every 3 wells. Cells were visually monitored to ensure single cell colonies and expanded into 24-well plates. To determine editing, genomic DNA was isolated from clones using the Agencourt DNAdvance kit (Beckman Coulter) according to manufacturer’s instructions. Clones were screened by ddPCR and verified by Sanger sequencing.
T cell – TRAC / B2M
Streptococcus pyogenes guide RNAs targeting the B2M and TRAC loci were generated by in vitro transcription of a PCR product using the T7-ScribeTM Standard RNA IVT Kit (CELLSCRIPT) following the manufacturer’s protocol. The PCR product for the in vitro transcription reaction was generated using plasmid PLA380 as a template and the indicated forward primers (OLI7076 for B2M, and OLI7077 for TRAC) with a common reverse primer (OLI7078).
The in vitro transcribed guide RNAs were complexed to wild type Cas9 protein at a molar ratio of 2:1 to generate ribonucleoprotein (RNP). The complexation integrity was evaluated by differential scanning fluorimetry (DSF). In brief, 5 μL of complexed RNP was diluted in 5 μL 2X Dye Mix. The 2X Dye Mix was generated from the 5000X stock SYPRO Orange Protein Gel Stain dye (Life Technologies, S6651) in 10× HEPES-Saline solution with MgCl2 (Boston Bio Products, C-6767) diluted to 1X in nuclease-free water. The complexed samples and uncomplexed protein controls were placed in a 384-well plate and placed in a BioRad thermocycler using the following protocol: 1 min at 20 °C, Melt Curve from 20 °C to 95 °C with increment changes of 1 °C, 1 min at 4 °C. Successful complexation is defined as a clear temperature shift between uncomplexed control samples and complexed RNP.
Human T cells were isolated from buffy coats using Miltenyi CD4 microbeads following the manufacturer’s protocol. On day 0, T-cells were activated using Dynabeads® Human T-Activator CD3/CD28 for T Cell Expansion and Activation (ThermoFisher Scientific). Beads were removed on day 2. On day 4, cells were counted using Trypan Blue (ThermoFisher Scientific) and TC20™ Automated Cell Counter (Bio-Rad) according to manufacturer’s protocol. For each condition 500,000 T cells were resuspended in 22 μL of Primary Cell Nucleofector Solution P2 (Lonza) containing 2 μM of total RNP. Samples were transferred to 16-well NucleocuvetteTM Strips (Lonza) and electroporated using program DS130 of the 4D-NucleofectorTM System (Lonza). Cells were subsequently transferred to untreated 96-well round bottom plates and cultured in 200 μL of X-Vivo 15 media (Lonza) containing 5% Human AB Serum (Gemini BioProduct), 1.6 mg/mL N-acetylcysteine (Sigma), 2 mM L-alanyl-L-glutamine (Thermo Scientific), 50 IU/ mL IL-2 (Peprotech), 5 ng/mL IL-7 (Peprotech) and 0.5 ng/mL IL-15 (Peprotech).
Four days post nucleofection, samples were pelleted and genomic DNA purified using the Agencourt DNAdvance kit (Beckman Coulter) according to the manufacturer’s protocol.
The CRISPR targeted genomic region of TRAC was PCR amplified for subsequent Sanger sequencing using primers OLI11371 and OLI11403. Amplification was performed in a 50 μl reaction volume, consisting of 10 μL of 5X Phusion HF buffer, 0.5 μM forward primer, 0.5 μM reverse primer, 200 μM dNTP, 1.5 μL DMSO, 0.5 μl of Phusion polymerase and 25 ng of gDNA template. PCR conditions were as follows: 30 s at 98 °C for initial denaturation, followed by 40 cycles of 10s at 98 °C for denaturation, 15 s at 64 °C for annealing, 30s at 72 °C for extension, and 5 min at 72 °C for the final extension. The PCR product was purified using (1.8×) Agencourt AMPure XP beads (Beckman Coulter Agencourt AMPure XP - PCR Purification #A63882) as per the manufacturer’s protocol. The amplified locus fragments were then cloned into pCR4-TOPO vectors using the ZeroBlunt TOPO Cloning Kit (Life Technologies Zero Blunt TOPO PCR Cloning Kit for Sequencing with One Shot TOP10 Chemically Competent E. coli #K287540) and transformed in One Shot Top10 chemically competent Escherichia coli cells. Cells were plated on Carbenicillin LB agar plates and incubated overnight at 37 °C. Plasmid DNA from 96 colonies per sample was sequenced by Genewiz, Inc. using an M13 reverse primer. Analysis of indel rates was done with the Geneious Software package (Biomatters, https://www.geneious.com).
For Illumina amplicon sequencing, two rounds of amplification were performed: round 1 targets the TRAC region, and round 2 adds the full-length Illumina adapter sequence. Round 1 was performed in a 12 μl reaction volume, consisting of 6 μL of NEBNext® Ultra™ II Q5® Master Mix (New England Biolabs), 0.125 μM forward primer (OLI4610), 0.125 μM reverse primer (OLI4611), and 20 ng of gDNA template. PCR conditions were as follows: 30 s at 98 °C for initial denaturation, followed by 20 cycles of 10s at 98 °C for denaturation, 15 s at 60 °C for annealing, 30s at 72 °C for extension, and 5 min at 72 °C for the final extension. The PCR product was purified using (0.9×) Agencourt AMPure XP beads (Beckman Coulter Agencourt AMPure XP - PCR Purification #A63882) as per the manufacturer’s protocol. Round 2 was performed in a 12 μl reaction volume, consisting of 6 μL of NEBNext® Ultra™ II Q5® Master Mix (New England Biolabs), 0.5 μM forward primer (Illumina forward), 0.5 μM reverse primer (Illumina reverse), and 20 ng of gDNA template. PCR conditions were as follows: 30 s at 98 °C for initial denaturation, followed by 20 cycles of 10s at 98 °C for denaturation, 15 s at 60 °C for annealing, 30s at 72 °C for extension, and 5 min at 72 °C for the final extension. The PCR product was purified using (0.9×) Agencourt AMPure XP beads (Beckman Coulter Agencourt AMPure XP - PCR Purification #A63882) as per the manufacturer’s protocol followed by size 300-1200 bp size selection on the BluePippin (Sage Science, Beverly, MA) and loaded on the Illumina MiSeq with 10% phiX. Analysis of indel rates was done as described in Bothmer et al. .
UDiTaS was performed according to the detailed protocol with the following modification. After tagmentation, the enzyme was inactivated with the addition of 1 μL of 0.2% SDS, pipette mixing and 5 min room temperature incubation. The tagmented DNA was added directly into round 1 PCR using primers OLI6259 and OLI6256.
Plasmid sensitivity experiments
For CEP290 plasmid-based sensitivity experiments PLA370, PLA367, and PLA379 were used (Additional file 8: Figure S7a). For the TRAC/B2M translocation plasmid-based sensitivity experiments PLA365, PLA366, PLA377, and PLA378 were used (Additional file 8: Figure S7b). Plasmid concentrations were determined using a NanoDrop2000 Spectrophotometer and working dilutions of 10 ng/μl were generated for all plasmids. In brief, for CEP290 sensitivity experiments the first sample consists of a 50% mix of PLA370 (Inversion) and PLA367 (Large Deletion) and contains no control plasmid (PLA379). Subsequently, 10 dilutions were generated by serially diluting the PLA370/PLA367 mix (sqrt10 dilution factor) into control plasmids (PLA379) maintaining a total plasmid concentration of 10 ng/μL throughout the different dilutions. The last sample consisted of only control plasmids (PLA379). For TRAC/B2M sensitivity experiments the first sample consists of a 50% mix of PLA365 and PLA366 and contains no control plasmids (PLA377/PLA378). Subsequently, 10 dilutions were generated by serially diluting PLA366 (sqrt10 dilution factor) into an equal mix of control plasmids PLA377/PLA378 maintaining a total plasmid concentration of 10 ng/μL throughout the different dilutions. The last sample consisted of only control plasmids (equal mix between PLA377/PLA378). All samples were subsequently subjected to UDiTaS: TRAC/B2M translocation samples were amplified with OLI6256 and OLI6259, while CEP290 plasmids were amplified with OLI6062. The UDiTaS protocol was applied with the following modifications: 6 cycles for first and second round PCR.For plasmids spiked into mouse DNA, different amounts of unique plasmids were mixed into mouse gDNA. For CEP290 plasmid spike in experiments, plasmids with the estimated copy number shown in Table 3 were spiked into 10 ng/μL mouse gDNA. For TRAC/B2M plasmid spike in experiments, plasmids were spiked into 10 ng/μL mouse gDNA as described in Table 4.
Anchored multiplex PCR sequencing
Genome-wide, unbiased identification of DSBs enabled by sequencing
High throughput, genome-wide translocation sequencing
Leber’s Congenital Amaurosis Type 10
Lower Limit of Detection
Next Generation Sequencing
Non-homologous end joining
Polymerase Chain Reaction
Recombinant DNA Advisory Committee
- RNP complex:
S. aureus Cas9
single guide RNAs
S. pyogenes CAS9
Unidirectional Targeted Sequencing
Unique molecule identifier
Mashal RD, Koontz J, Sklar J. Detection of mutations by cleavage of DNA heteroduplexes with bacteriophage resolvases. Nat Genet. 1995;9:177–83.
Vouillot L, Thelie A, Pollet N. Comparison of T7E1 and surveyor mismatch cleavage assays to detect mutations triggered by engineered nucleases. G3. 2015;5:407–15. https://doi.org/10.1534/g3.114.015834.
Hendel A, Fine EJ, Bao G, Porteus MH. Quantifying on- and off-target genome editing. Trends Biotechnol. 2015;33:132–40. Elsevier Ltd
Tycko J, Myer VE, Hsu PD. Methods for optimizing CRISPR-Cas9 genome editing specificity. Mol Cell. 2016;63:355–70. Elsevier Inc
Brinkman EK, Chen T, Amendola M, van Steensel B. Easy quantitative assessment of genome editing by sequence trace decomposition. Nucleic Acids Res. 2014;42:e168.
C D (NIHOD), Atkins MB, Paula C, T CW, Patrick H, Peter H, et al. RAC review: phase I trial of autologous T cells engineered to express NYESO-1 TCR and gene edited to eliminate endogenous TCR and PD-1 2016.
Witten CM. Update on the Office of Cellular , tissue , and gene therapies. La Jolla, California: Cell and Gene Meeting on the Mesa; 2016. https://www.youtube.com/watch?v=MpxTSCY6HjM&index=8&list=PLbTBF__p6d_KL57ZH06cC1W3SgInm4YIv. Accessed 8 Feb 2018
Zheng Z, Liebers M, Zhelyazkova B, Cao Y, Panditi D, Lynch KD, et al. Anchored multiplex PCR for targeted next-generation sequencing. Nat Med. 2014;20:1479–84. Nature Publishing Group
Chiarle R, Zhang Y, Frock RLL, Lewis SMM, Molinie B, Ho Y-J, et al. Genome-wide translocation sequencing reveals mechanisms of chromosome breaks and rearrangements in B cells. Cell. 2011;147:107–19. https://doi.org/10.1016/j.cell.2011.07.049. Elsevier Inc
Hu J, Meyers RM, Dong J, Panchakshari RA, Alt FW, Frock RL. Detecting DNA double-stranded breaks in mammalian genomes by linear amplification-mediated high-throughput genome-wide translocation sequencing. Nat Protoc. 2016;11:853–71. https://doi.org/10.1038/nprot.2016.043. Nature Publishing Group
Tsai SQ, Zheng Z, Nguyen NT, Liebers M, Topkar VV, Thapar V, et al. GUIDE-seq enables genome-wide profiling of off-target cleavage by CRISPR-Cas nucleases. Nat. Biotechnol. 2014;33:187–97.
Chen L, Liu P, Evans TC, Ettwiller LM. DNA damage is a pervasive cause of sequencing errors, directly confounding variant identification. Science. 2017;355:752–6. https://doi.org/10.1126/science.aai8690.
Costello M, Pugh TJ, Fennell TJ, Stewart C, Lichtenstein L, Meldrim JC, et al. Discovery and characterization of artifactual mutations in deep coverage targeted capture sequencing data due to oxidative DNA damage during sample preparation. Nucleic Acids Res. 2013;41:1–12.
Adey A, Morrison HG, Asan, Xun X, Kitzman JO, Turner EH, et al. Rapid, low-input, low-bias construction of shotgun fragment libraries by high-density in vitro transposition. Genome Biol. 2010;11:R119. https://doi.org/10.1186/gb-2010-11-12-r119. BioMed Central Ltd
Picelli S, Björklund ÅK, Reinius B, Sagasser S, Winberg G, Sandberg R. Tn5 transposase and tagmentation procedures for massively scaled sequencing projects. Genome Res. 2014;24:2033–40.
Wang Q, Gu L, Adey A, Radlwimmer B, Wang W, Hovestadt V, et al. Tagmentation-based whole-genome bisulfite sequencing. Nat Protoc. 2013;8:2022–32.
Drivas TG, Wojno AP, Tucker BA, Stone EM, Bennett J. Basal exon skipping and genetic pleiotropy: a predictive model of disease pathogenesis. Sci Transl Med. 2015;7:291ra97. https://doi.org/10.1126/scitranslmed.aaa5370.
den Hollander AI, Koenekoop RK, Yzer S, Lopez I, Arends ML, Voesenek KEJ, et al. Mutations in the CEP290 (NPHP6) gene are a frequent cause of Leber congenital amaurosis. Am J Hum Genet. 2006;79:556–61.
Jiang J, Zhang L, Zhou X, Chen X, Huang G, Li F, et al. Induction of site-specific chromosomal translocations in embryonic stem cells by CRISPR/Cas9. Sci Rep. 2016;6:21918. https://doi.org/10.1126/scitranslmed.aaa5370. Nature Publishing Group
Ghezraoui H, Piganeau M, Renouf B, Renaud JB, Sallmyr A, Ruis B, et al. Chromosomal translocations in human cells are generated by canonical nonhomologous end-joining. Mol Cell. 2014;55:829–42. Elsevier Inc
Frock RL, Hu J, Meyers RM, Ho Y-J, Kii E, Alt FW. Genome-wide detection of DNA double-stranded breaks induced by engineered nucleases. Nat Biotechnol. 2015;33:179–86.
Martin M. Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnetjournal. 2011;17:10–2. https://doi.org/10.14806/ej.17.1.200.
Langmead B, Salzberg SL. Fast gapped-read alignment with bowtie 2. Nat Methods. 2012;9:357–9.
Li H, Handsaker B, Wysoker A, Fennell T, Ruan J, Homer N, et al. The sequence alignment/map format and SAMtools. Bioinformatics. 2009;25:2078–9.
Bothmer A, et al. Characterization of the interplay between DNA repair and CRISPR/Cas9-induced DNA lesions at an endogenous locus. Nat Commun. 2017;8:13905.
Thorvaldsdóttir H, Robinson JT, Mesirov JP. Integrative genomics viewer (IGV): high-performance genomics data visualization and exploration. Brief Bioinform. 2013;14:178–92.
Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, et al. Integrative genomics viewer. Nat Biotechnol. 2011;29:24–6.
The authors would like to thank all the people at Editas for helpful discussions and support. We would especially would like to thank Cecilia Cotta for translocation insights and discussion about translocations, Christina Lee for method development, William Selleck for Tn5 protein purification, Tim Fennell for primer design algorithms, and Mike Dinsmore for informatics support.
This work was funded by Editas Medicine.
Availability of data and materials
All NGS data are available via the NCBI-SRA database: PRJNA433666. Analysis code is available in https://github.com/editasmedicine/uditas.
Ethics approval and consent to participate
Consent for publication
All authors were employed by Editas Medicine when the work was carried out. No other interests are declared.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Figure S1. Schematic of the bioinformatics pipeline for UDiTaS analysis. (PPTX 61 kb)
Detailed Protocols for UDiTaS and AMP-Seq methods. (DOCX 237 kb)
Figure S2. Example editing events. Examples of various editing events in the U-2 OS bulk editing experiment shown in the Integrated Genome Viewer (IGV) [26, 27]. A schematic on top of each view depicts the observed editing event. Reads colored in red/blue were aligned to the top/bottom genomic reference DNA sequence. Note that small indels are observed in addition to the junctions formed from the larger structural changes. These indels likely arose due to repair pathway activity prior to rearrangement. a. 323 site small indels. b. 323-64 large desired 1.1 kb deletion junction. c. 323-64 large desired 1.1 kb inversion junction. d. 323 homologous junction. (PPTX 147 kb)
Figure S3. UDiTaS reproducibility. Identical samples were run in UDiTaS using either SDS addition or Zymo column purification after tagmentation. Measured values for the various constructs are reproducible and highly correlated across a wide range of concentrations. (PPTX 3928 kb)
Figure S4. Binomial power calculation applied to UDiTaS. A simulated binomial distribution, plotting editing frequency (e.g.: probability of success) vs. number of unique molecular identifiers (e.g.: trials) for a given number of expected observations (1, 2, or 3). Graphs on the left are 95% confidence and right 99% confidence. (PPTX 86 kb)
Figure S5. Genome mapping rates for UDiTaS. Individual reads map to the expected genome site with high frequency indicating the robustness of the assay. Ten distinct samples for primer OLI6062 are plotted on the x-axis and the y-axis shows the percentage or reads mapping to the expected reference amplicon for each sample. (PPTX 3767 kb)
Figure S6. UDiTaS characterization and comparison to AMP-Seq with plasmid standards. Plasmids containing the CEP290 structural variants a. or the TRAC-B2M balanced translocation b. and c. were synthesized and contain engineered unique SNPs in the insert to identify the plasmid after sequencing. The plasmids were diluted at various levels into mouse genomic DNA and processed through UDiTaS and AMP-Seq using primers for CEP290 a., B2M b. and TRAC c. The number of input plasmids versus the number of plasmids detected is plotted for both UDiTaS and AMP-Seq. Linear regression models and 95% confidence model predictions are displayed on the plots. The parameter β determines the linearity of the method, with values close to 1 indicating more linearity. We used ANOVA p-values to examine differences in β for UDiTaS and AMP-Seq. Below each plot, the table displays the total number of fastq reads sequenced in the reaction, the number of reads mapped to the wild-type amplicon (the most abundant one) and the final number of UMIs counted, for both UDiTaS and AMP-Seq. At all tested loci, UDiTaS shows greater linearity and number of UMIs detected when compared to AMP-Seq. (PPTX 12722 kb)
Figure S7. UDiTaS characterization of plasmid standards without carrier DNA. To ensure that the carrier mouse genomic DNA was not influencing the UDiTaS reaction, additional sets of UDiTaS reactions were run with plasmids in the absence of any carrier DNA. a. CEP290 plasmids with the Wild Type, Large Deletion, and Large Insertion (PLA379, PLA367, and PLA370) and b. B2M-TRAC plasmids with the B2M, TRAC, and both balanced translocations (PLA377, PLA378, PLA365, and PLA366) were diluted as described in the methods. The DNA plasmids mixtures were process through UDiTaS and the analysis pipeline. Plotted is the expected frequency for a given structural variant vs. measured frequency for a structural variant (x = y is the grey line). Accuracy and linearity appear to be excellent for both loci with all four primers, with an LLOD of ~ 0.01%-0.1%. (PPTX 991 kb)
Figure S8. Comparison of Indel rates between UDiTaS and other methods. T-Cells edited with the TRAC + B2M guides were analyzed for indel editing at the TRAC locus using PCR-amplification followed by Sanger Sequencing or NGS, in addition to UDiTaS with two different anchor primers. Indel rates were very similar between the methods. (PPTX 81 kb)
About this article
Cite this article
Giannoukos, G., Ciulla, D.M., Marco, E. et al. UDiTaS™, a genome editing detection method for indels and genome rearrangements. BMC Genomics 19, 212 (2018) doi:10.1186/s12864-018-4561-9
- Gene editing
- Next generation sequencing
- Translocation detection