Multiplex CRISPR/Cas9-based genome editing for correction of dystrophin mutations that cause Duchenne muscular dystrophy

Abstract

The CRISPR/Cas9 genome-editing platform is a promising technology to correct the genetic basis of hereditary diseases. The versatility, efficiency and multiplexing capabilities of the CRISPR/Cas9 system enable a variety of otherwise challenging gene correction strategies. Here, we use the CRISPR/Cas9 system to restore the expression of the dystrophin gene in cells carrying dystrophin mutations that cause Duchenne muscular dystrophy (DMD). We design single or multiplexed sgRNAs to restore the dystrophin reading frame by targeting the mutational hotspot at exons 45-55 and introducing shifts within exons or deleting one or more exons. Following gene editing in DMD patient myoblasts, dystrophin expression is restored in vitro. Human dystrophin is also detected in vivo after transplantation of genetically corrected patient cells into immunodeficient mice. Importantly, the unique multiplex gene-editing capabilities of the CRISPR/Cas9 system facilitate the generation of a single large deletion that can correct up to 62% of DMD mutations.

ARTICLE

Received 16 Nov 2014 | Accepted 8 Jan 2015 | Published 18 Feb 2015

The CRISPR/Cas9 genome-editing platform is a promising technology to correct the genetic basis of hereditary diseases. The versatility, efciency and multiplexing capabilities of the CRISPR/Cas9 system enable a variety of otherwise challenging gene correction strategies. Here, we use the CRISPR/Cas9 system to restore the expression of the dystrophin gene in cells carrying dystrophin mutations that cause Duchenne muscular dystrophy (DMD). We design single or multiplexed sgRNAs to restore the dystrophin reading frame by targeting the mutational hotspot at exons 4555 and introducing shifts within exons or deleting one or more exons. Following gene editing in DMD patient myoblasts, dystrophin expression is restored in vitro. Human dystrophin is also detected in vivo after transplantation of genetically corrected patient cells into immunodecient mice. Importantly, the unique multiplex gene-editing capabilities of the CRISPR/Cas9 system facilitate the generation of a single large deletion that can correct up to 62% of DMD mutations.

DOI: 10.1038/ncomms7244

Multiplex CRISPR/Cas9-based genome editing for correction of dystrophin mutations that cause Duchenne muscular dystrophy

David G. Ousterout1, Ami M. Kabadi1, Pratiksha I. Thakore1, William H. Majoros2, Timothy E. Reddy3,4

& Charles A. Gersbach1,3,5

1 Department of Biomedical Engineering, Duke University, Room 136 Hudson Hall, Box 90281, Durham, North Carolina 27708, USA. 2 Program in Computational Biology and Bioinformatics, Duke University, Durham, North Carolina 27708, USA. 3 Center for Genomic and Computational Biology, Duke University, Durham, North Carolina 27708, USA. 4 Department of Biostatistics and Bioinformatics, Duke University Medical Center, Durham,North Carolina 27710, USA. 5 Department of Orthopaedic Surgery, Duke University Medical Center, Durham, North Carolina 27710, USA. Correspondence and requests for materials should be addressed to C.A.G. (email: mailto:[email protected]

Web End [email protected] ).

NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 1

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244

Genome-editing technologies use synthetic nucleases to induce cellular DNA repair mechanisms and introduce site-specic, predened genetic modications in complex

genomes1. These engineered enzymes are commonly based on zinc nger nucleases (ZFNs)1, transcription activator-like effector nucleases (TALENs)2, meganucleases3 and most recently the RNA-guided CRISPR/Cas9 system410. The nucleases create site-specic double-strand breaks at predened genomic sites that stimulate either non-homologous end joining (NHEJ) for targeted gene disruption or homologous recombination for highly efcient gene targeting. The simplicity and versatility of the CRISPR/Cas9 genome-editing system has led to rapid adoption and expansion of this technology that has proven to be remarkably robust for manipulating gene sequences in human cells. This has enabled new possibilities such as efcient multiplex gene editing for simultaneously inactivating multiple genes5,6,11. In this study, we apply the CRISPR system to repair genes mutated in hereditary disease, including capitalizing on the unique multiplex capacity of this technology to create large genomic deletions that restore gene expression.

CRISPR/Cas9 systems have been adapted from multiple bacterial species, including Streptococcus pyogenes, S. thermophilus and Neisseria meningitidis, to efciently generate targeted gene modications in human cells510,12. These systems consist of a Cas9 nuclease that is co-expressed with a single guide RNA (sgRNA) molecule. Cas9 forms a complex with the 30 end of the sgRNA, and the proteinRNA pair recognizes its genomic target by complementary base pairing between the 50 end of the sgRNA sequence and a predened 20 bp DNA sequence, known as the protospacer. By simply exchanging the 20-bp recognition sequence of the expressed sgRNA, the Cas9 nuclease can be directed to new genomic targets. The only restriction for protospacer targeting is that the sequence must be immediately followed by the protospacer-adjacent motif (PAM), a short sequence recognized by the Cas9 nuclease that is required for DNA cleavage. Several studies of the S. pyogenes CRISPR system have dened the PAM sequence for this Cas9 (SpCas9) as 50-NGG-30 and characterized the specicity of this system in human cells1319. A unique capability of the CRISPR/Cas9 system is the straightforward ability to simultaneously target multiple distinct genomic loci by co-expressing a single Cas9 protein with two or more sgRNAs5,6,11,20.

One of the most promising applications of genome editing is the correction of genetic mutations associated with hereditary disease14. Duchenne muscular dystrophy (DMD) is the most common hereditary disease; no effective treatments exist for this disorder. DMD is a severe X-linked disease that presents with progressive muscle wasting that typically leads to loss of ambulation in the second decade and death within the third decade of life because of respiratory complications or heart failure. The molecular basis of DMD is a mutation in the dystrophin gene21 that leads to the complete lack of function of this essential skeletal muscle protein. These mutations are most commonly frameshifts generated by large intragenic deletions of one or more exons. DMD is the prototypical example of a group of monogenetic hereditary diseases that can be corrected by removing internal, but unessential, regions of the mutated gene to restore the proper reading frame22,23. For example, there is a class of common deletions in the exon 4555 mutation hotspot region of the dystrophin gene that maintain the correct reading frame and lead to the expression of a truncated, but functional, dystrophin protein. Patients with this class of mutations are often asymptomatic or display mild symptoms associated with Becker muscular dystrophy, a substantially less severe disease than DMD. This has led to signicant interest in developing an oligonucleotide-mediated exon skipping strategy that will restore

the dystrophin reading frame during mRNA processing and convert DMD to a Becker-like phenotype22. Whereas early clinical trials in this area have focused on skipping exon 51 (refs 24,25), which is applicable to 13% of DMD patients, other preclinical efforts have demonstrated multi-exon skipping of the complete exon 4555 coding region with a combination treatment of up to ten oligonucleotides26,27 that could potentially address greater than 60% of known DMD patient mutations23. However, there are signicant technical and practical hurdles to designing and developing this type of complex combination therapy, in addition to the general challenges of developing any oligonucleotide-based therapy that must be continuously administered for the lifetime of the patient. In contrast to these transient mRNA-targeted oligonucleotide-mediated exon skipping strategies, genome editing has the ability to make precise and permanent changes to gene sequences that will persist after cell division. In addition, only two nucleases are necessary to delete a genomic region of any length, in contrast to exon skipping in which a distinct oligonucleotide must be designed for each exon to be removed from the mRNA transcript.

Genome editing using various designer nucleases has been proposed as a promising method to restore the native dystrophin gene in DMD patient cells2830. However, an obstacle to implementing this approach has been successfully engineering the multiple nucleases targeted to the exons and introns necessary to address a large fraction of the DMD patient population. In this study, we take advantage of the versatility of the S. pyogenes CRISPR/Cas9 system to show proof-of-concept that this approach can rapidly and efciently generate targeted frameshifts and large deletions to address commonly occurring mutations in the dystrophin gene across exons 4555. Skeletal myoblasts from DMD patients were treated with sgRNAs and SpCas9 to correct patient-specic mutations and edited cells were enriched by uorescence-activated cell sorting. Gene editing by CRISPR/Cas9 resulted in restored dystrophin mRNA transcripts and protein expression. Signicantly, we generated a large deletion of 336 kb across a mutational hotspot containing exons 4555 that is applicable to correction of greater than 60% of DMD patient mutations. This genomic deletion resulted in the loss of exons 4555 in the corresponding dystrophin transcript and restored dystrophin expression in human DMD cells. In addition, an enriched pool of gene-corrected cells demonstrated expression of human dystrophin in vivo following engraftment into immunodecient mice. CRISPR/Cas9 gene editing did not have signicant toxic effects in human myoblasts as observed by stable gene editing frequencies and minimal cytotoxicity of several sgRNAs. However, gene-editing activity was conrmed at 3 out of 50 predicted off-target sites across ve sgRNAs and CRISPR/Cas9-induced chromosomal translocations between on-target and off-target sites were detectable in immortalized cell lines. These data indicate a need to increase the specicity of this technology, although only a fraction of these events were detectable in the treated myoblasts from DMD patients. Collectively, this study provides proof-of-principle that the CRISPR/Cas9 technology is a versatile method for correcting a signicant fraction of dystrophin mutations and with continued development may serve as a general platform for treating genetic disease.

ResultsTargeting hotspot mutations in the human dystrophin gene. To utilize the CRISPR/Cas9 gene-editing platform for correcting a wide range of dystrophin mutations, we created dozens of sgRNAs targeted to the hotspot mutation region between exons 45 and 55 (Fig. 1). We selected the previously described

2 NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244 ARTICLE

a

Exon 4555 deletion (336,380 bp)

Exon 51 deletion (8001,050 bp)

44 45 46 47 48 49 50

Panel of sgRNAs to create targeted frameshifts at 5 or 3 end of each exon

Out-of-frame (4850)

b

PAM

*

Exon 51

Exon 47 Exon 51

Exon 51

Exon 47

PAM

Frame restored by small targeted deletion

Frame restored by genetic deletion of a specific exon

Frame restored by genetic deletion of an entire region of exons

c

sgRNA 1 sgRNA 2

sgRNA 1 sgRNA 2

Exon 52

Exon 47

Out-of-frame (4850)

Out-of-frame (various deletions)

Exon 47

Exon 52

51 52 53 54 55 56

d

Exon 44

Exons 4555

Exon 56

Exon 44 Exon 56

Figure 1 | CRISPR/Cas9 targeting of the dystrophin gene. (a) sgRNA sequences were designed to bind sequences in the exon 4555 mutational hotspot region of the dystrophin gene, such that gene editing could restore dystrophin expression from a wide variety of patient-specic mutations. Arrows within introns indicate sgRNA targets designed to delete entire exons from the genome. Arrows within exons indicate sgRNA targets designed to create targeted frameshifts in the dystrophin gene. (b) Example of frame correction following introduction of small insertions or deletions by NHEJ DNA repair in exon 51 using the CR3 sgRNA. (c) Schematic of multiplex sgRNA targets designed to delete exon 51 and restore the dystrophin reading frame in a patient mutation with the deletion of exons 4850. (d) Schematic of multiplex sgRNA targets designed to delete the entire exon 4555 region to address a variety of DMD patient mutations.

S. pyogenes system that utilizes a human codon-optimized SpCas9 nuclease5,6 and a chimeric single guide RNA (sgRNA) expression vector to direct efcient site-specic gene editing. Similar to our previous study targeting exon 51 with TALENs28, we selected protospacers to target the 50 and 30 ends of exons 4555, which meet the 50-NGG-30 PAM requirement of SpCas9. Small insertions or deletions created by NHEJ-based DNA repair within these exons can generate targeted frameshift mutations that address various dystrophin mutations surrounding each exon (Fig. 1a,b). For example, CR3 was designed to correct dystrophin mutations or deletions surrounding exon 51 by introducing small insertions or deletions in the 50 end of exon 51 to restore the downstream dystrophin reading frame (Fig. 1b). In addition, we designed sgRNAs to employ the multiplex capability of the CRISPR/Cas9 system and specically delete individual exons or a series of exons to restore the dystrophin reading frame, similar to the methods of oligonucleotide-based exon skipping. For this purpose, sgRNAs were targeted to the intronic regions surrounding exon 51 (Fig. 1c) or exons 4555 (Fig. 1d). These

sgRNAs were intentionally targeted to sites nearest to the downstream or upstream exon intended to be included in the resulting transcript to minimize the likelihood that the background patient deletion would include the intronic sgRNA target sites.

Screening sgRNAs targeted to the human dystrophin gene. We initially assessed gene-editing frequency in the human HEK293T cell line to rapidly determine different sgRNA targeting efciencies. As quantied by the Surveyor assay 3 days post transfection, we found that 29 out of 32 (B90%) of sgRNAs tested were able to mediate efcient gene modication at the intended locus (Table 1 and Supplementary Fig. 1). The gene-editing frequencies were stable for almost all of the sgRNAs (o25% signal change from day 3 to day 10, Table 1 and Supplementary Fig. 2), indicating that gene editing mediated by each individual sgRNA was well tolerated. An exception is CR33, which had no detectable activity at day 10, although activity may be below the sensitivity of the Surveyor assay (estimated B1%).

NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 3

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244

Table 1 | Measured activity of sgRNAs in human cells.

Target sgRNA # % Modied alleles at day 3 % Modied alleles at day 10 % Change day 10/day 3 Multiplex deletion of exon 51

Int 50 CR1 6.6 9.3 41.8 Int 50 CR2 10.3 14.0 36.2 Ex 51 CR4 11.9 14.4 21.3 Int 51 CR5 12.4 13.3 7.8

Multiplex deletion of exons 4555Int 44 CR6 16.1 16.9 4.3 Int 44 CR33 1.3 o1 n.d.

Int 44 CR34 13.2 11.0 16.6

Int 55 CR7 6.8 7.1 5.3 Int 55 CR35 22.5 20.9 7.1

Int 55 CR36 26.4 24.7 6.4

Targeted frameshiftsEx 45 CR10 14.9 16.3 9.3 Ex 45 CR11 o1 o1 n.d.

Ex 46 CR12 o1 o1 n.d. Ex 46 CR13 16.9 18.4 9.2

Ex 47 CR14 17.2 17.6 2.9 Ex 47 CR15 15.4 15.3 0.9

Ex 48 CR16 11.5 10.9 5.0

Ex 48 CR17 o1 o1 n.d. Ex 49 CR18 1.8 2.2 20.1

Ex 49 CR19 33.7 38.4 13.9 Ex 50 CR20 14.9 13.7 7.6

Ex 50 CR21 24.1 20.8 13.5

Ex 51 CR3 13.0 16.7 28.0 Ex 51 CR31 18.9 16.9 10.2

Ex 52 CR22 25.9 20.3 21.6

Ex 52 CR23 25.2 24.0 4.8

Ex 53 CR24 24.8 23.6 4.6

Ex 53 CR25 2.6 2.9 9.5 Ex 54 CR26 24.5 22.0 10.1

Ex 54 CR27 13.4 12.6 5.9

Ex 55 CR28 21.6 19.8 8.4

Ex 55 CR29 19.2 19.6 2.2

sgRNA, single guide RNA.

HEK293Ts were transfected with constructs encoding human codon-optimized SpCas9 and the indicated sgRNA. Each sgRNA was designed to modify the dystrophin gene as indicated. The frequency of gene modication at day 3 or day 10 post transfection was determined by the Surveyor assay. The ratio of measured Surveyor signal at day 3 and day 10 was calculated to quantify the stability of gene editing frequencies for each sgRNA in human cells.

Enrichment of corrected cells. We next sought to use selected sgRNAs to correct specic mutations in DMD patient myoblast cell lines. After transfection into DMD myoblasts, we observed unexpectedly low or undetectable gene modication activity as measured by the Surveyor assay (Fig. 2c, bulk population). Therefore, we used ow cytometry to select for transfected cells co-expressing green uorescent protein (GFP) through a 2A ribosomal skipping peptide linked to the SpCas9 protein (Fig. 2a), similar to previously described methods to enrich gene-modied cell populations using uorescent reporters31,32. Importantly, the addition of this uorescent reporter to the SpCas9 expression vector did not signicantly impact gene-editing activity in HEK293T cells (Fig. 2b). A low percentage of transfected myoblasts (B0.52%) expressed the uorescent reporter at 3 days after electroporation, despite high transfection efciencies of a control GFP expression plasmid (typically 470%, Fig. 2d, pmaxGFP). Given the high levels of

CRISPR/Cas9 activity in the easily transfected HEK293T cell line and previous studies demonstrating robust expression from the cytomegalovirus promoter in skeletal myoblasts, inefcient plasmid delivery by electroporation or vector-specic effects were likely responsible for the low observed gene editing efciencies in unsorted cells. After sorting the GFP-positive DMD myoblasts, we observed a substantial increase in detectable activity at most sgRNA target loci

(Fig. 2c). Therefore, all subsequent experiments used cells sorted for SpCas9 expression by expression of this uorescent reporter.

Restoration of dystrophin expression by targeted frameshifts. We have shown previously that small insertions and deletions created by NHEJ DNA repair can be used to create targeted frameshifts to correct aberrant reading frames28. Similar to this approach, we designed a sgRNA, CR3, to restore the dystrophin reading frame by introducing small insertions and deletions within exon 51 (Figs 1b and 3a). The types of insertions and deletions generated by CRISPR/Cas9 at this locus were assessed by Sanger sequencing of alleles from the genomic DNA of HEK293T cells co-transfected with expression plasmids for SpCas9 and the CR3 sgRNA (Fig. 3b). Notably, the insertions and deletions resulted in conversion to all three reading frames (Fig. 3b,c), consistent with our previous results using TALENs28. To demonstrate genetic correction in a relevant patient cell line, expression plasmids for SpCas9 and the CR3 sgRNA were electroporated into DMD myoblasts with a deletion of exons 4850 that is correctable by creating frameshifts in exon 51. The treated cells were sorted, veried to have gene modication activity by the Surveyor assay (CR3, Fig. 2c sorted population),

4 NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244 ARTICLE

a

b

CR1 +

CR2 CR3 CR4

+ + +

CMV T2A

SpCas9

eGFP

SpCas9

SpCas9

**

*

*

*

*

*

sgRNA-1

hU6

hU6

% Indels: 6.6 10.3 13.0 11.9

14.5

sgRNA-2

**

*

*

*

*

*

T2A-GFP

% Indels:

5.7

13.7

14.1

c

Bulk population

Sorted population

CR1 + + + + + + + + + + + +

CR2 CR3 CR5 CR6 CR36 CR1 CR2 CR3 CR5 CR6 CR36

** ** * **

*

*

*

% Indels: ND ND ND ND ND ND 0.8 5.0 6.5 2.1 1.0

ND

d

Empty vector

0.0%

104

pmaxGFP

SpCas9 + emptyvector SpCas9 + sgRNA

Autofluorescence

104

103

102

101

100

100 104

103

102

101

104

103

102

101

100

100 104

103

102

101

104

103

102

101

100

100 104

103

102

101

104

103

102

101

100

100

77.5% 1.2% 1.9%

101

102

103

GFP

Figure 2 | Fluorescence-activated ow sorting to enrich genetically modied DMD myoblasts. (a) A plasmid expressing a human codon-optimized SpCas9 protein linked to a GFP marker using a T2A ribosomal skipping peptide sequence was co-electroporated into human DMD myoblasts with one or two plasmids carrying sgRNA expression cassettes. (b) The indicated sgRNA expression cassettes were independently co-transfected into HEK293T cells with a separate plasmid expressing SpCas9 with (bottom) or without (top) a GFP marker linked to SpCas9 by a T2A ribosomal skipping peptide sequence. Gene modication frequencies were assessed at 3 days post transfection by the Surveyor assay. Asterisks indicate cleaved PCR amplicon fragments at the expected sizes indicative of indels generated by gene editing. (c) DMD myoblasts with deletions of exons 4850 in the dystrophin gene were treated with sgRNAs that correct the dystrophin reading frame in these patient cells. Gene modication was assessed at 20 days post-electroporation in unsorted (bulk) or GFP sorted cells. (d) GFP expression in DMD myoblasts 3 days after electroporation with indicated expression plasmids. Transfection efciencies and sorted cell populations are indicated by the gated region. ND, not determined.

and differentiated into myotubes to test for restored dystrophin expression. Expression of dystrophin protein was observed concomitant with the detectable nuclease activity (Fig. 3d and Supplementary Fig. 4). Taken together with the data from Table 1, the S. pyogenes CRISPR/Cas9 system presents a powerful method to quickly generate targeted frameshifts to address a variety of patient mutations and restore expression of the human dystrophin gene.

Multiplex CRISPR/Cas9 gene editing to delete exon 51. The multiplexing capability of the CRISPR/Cas9 system presents a novel method to efciently generate genomic deletions of specic exons for targeted gene correction. DMD patient myoblasts with background deletions correctable by exon 51 skipping were treated with two combinations of sgRNAs anking exon 51 (CR1/CR5 or CR2/CR5) and sorted to enrich for gene-edited cells as in Fig. 2. As detected by end-point PCR of the genomic DNA from these treated cells, the expected genomic deletions were only present when both sgRNAs were electroporated into the cells with SpCas9 (Fig. 4a). Sanger sequencing conrmed the expected junction of the distal chromosomal segments (Fig. 4b) for both

deletions. After differentiating the sorted myoblasts, a deletion of exon 51 from the mRNA transcript was detected only in the cells treated with both sgRNAs (Fig. 4c). Finally, restored dystrophin protein expression was detected in the treated cells concomitant with observed genome- and mRNA-level deletions of exon 51 (Fig. 4d and Supplementary Fig. 4).

Dystrophin rescue by a multi-exon large genomic deletion. Although addressing patient-specic mutations is a powerful use of the CRISPR/Cas9 system, it would be advantageous to develop a single method that can address a myriad of common patient deletions. For example, a promising strategy is to exclude the entire exon 4555 region as a method to correct up to 62% of known patient deletions22,26,27. Therefore, we tested if multiplex CRISPR/Cas9-based gene editing may be able to generate efcient deletion of the exon 4555 locus in human cells. After transfection into HEK293T cells, we detected the expected deletion of B336,000 bp by PCR of the genomic DNA (Fig. 5a). Similarly, we were able to detect this deletion by PCR of the genomic DNA from SpCas9/sgRNA-treated DMD patient cells harbouring a background deletion of exons 4850 of

NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 5

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244

a

Intron 50

Deletions

Insertions

Exon 51 PAM

sgRNA

5

b

Frame

(16) Length8

8

9 6 4 5 4 3 8

7 3

+9 +3

+3

+2

+3

(2)

(2)

(16)

12

21

+2

+3

+3

+3

+1

+2

+1

+3

+2

56

+2

+1

Length Frame

- sgRNA

c d

CR3

Total events: 17/33 (52%) +1 Frame: 3/17 (18%) +2 Frame: 7/17 (41%) +3 Frame: 7/17 (41%)

Dystrophin

GAPDH

Figure 3 | Targeted frameshifts to restore the dystrophin reading frame using CRISPR/Cas9. (a) The 50 region of exon 51 was targeted using a sgRNA, CR3, that binds immediately upstream of the rst out-of-frame stop codon. PAM: protospacer-adjacent motif. (b) The exon 51 locus was PCR amplied from HEK293T cells treated with SpCas9 and CR3 expression cassettes. Sequences of individual clones were determined by Sanger sequencing. The top sequence (bolded, exon in red) is the native, unmodied sequence. The number of clones for each sequence is indicated in parentheses. (c) Summary of total gene editing efciency and reading frame conversions resulting from gene modication shown in b.(d) Western blot for dystrophin expression in human DMD myoblasts treated with SpCas9 and the CR3 sgRNA expression cassette (Fig. 2c) to create targeted frameshifts to restore the dystrophin reading frame. Dystrophin expression was probed using an antibody against the rod-domain of the dystrophin protein after 6 days of differentiation.

unknown length (Fig. 5a). Sanger sequencing of this deletion band from the genomic DNA of treated DMD cells revealed the expected junctions of intron 44 and intron 55 immediately adjacent to the sgRNA target sites (Fig. 5b). After differentiation of treated DMD cells, the expected deletion of exons 4555 was detected in the dystrophin mRNA transcript and veried to be a fusion of exons 44 and 56 by Sanger sequencing (Fig. 5c). Based on the intensity of the modied mRNA band compared with the unmodied band, the efciency of this 336 kb deletion was less than the efciency of exon 51 deletion (Fig. 4a,c), which is expected based on previous studies showing a size-dependent decrease in the efciency of nuclease-mediated genomic deletions20,33. Nevertheless, restored protein expression was observed by western blot in the sorted cell populations containing the CRISPR/Cas9-induced deletion of exons 4555 from the genome and resulting mRNA transcripts (Fig. 5d and Supplementary Fig. 4). These data demonstrate that multiplex CRISPR/Cas9 editing presents a single universal method to

restore the dystrophin reading frame in more than 60% of DMD patient mutations.

Transplantation of corrected myoblasts into mice. A potential strategy for DMD therapy is to correct a population of autologous patient muscle progenitor cells that can be engrafted into the patients skeletal muscle tissue to rescue dystrophin expression3440. To demonstrate the ability of the corrected cells to express correctly localized human dystrophin in vivo, we transplanted a population of DMD myoblasts that were treated with sgRNAs CR1 and CR5, which ank exon 51, and sorted for expression of GFP as before (Figs 2 and 4 and Supplementary Fig. 3). After 4 weeks, muscle bres positive for human spectrin, which is expressed by both corrected and uncorrected cells, were detected in cryosections of injected muscle tissue (Fig. 6). Several of these bres were also positive for human dystrophin with expression properly localized to the sarcolemma, indicating correct function of the restored protein (Fig. 6 and Supplementary Fig. 5). The frequency of dystrophin-positive bres was low, potentially due to the effects of electroporation, extensive culture and sorting, the efciency of editing (Fig. 4a and Supplementary Fig. 3), overall dystrophin expression levels and/or the quality of available antibodies. Nevertheless, bres positive for both spectrin and dystrophin were found across several sections in multiple injected mice (Fig. 6 and Supplementary Fig. 5), validating consistent in vivo expression of the restored protein in a fraction of implanted cells. Importantly, no bres positive for human dystrophin were observed in any sections from mice injected with the untreated DMD myoblasts, conrming that the CRISPR/Cas9-modied cells were the source of human dystrophin expression.

Off-target and cytotoxicity analysis. We assessed the relative cytotoxicity of the CRISPR/Cas9 system in human cells to a subset of sgRNAs by adapting a ow cytometry-based GFP retention assay as previously described28. Minimal cytotoxicity was observed for SpCas9 co-expressed with or without sgRNAs after transfection into human cells (Fig. 7a). Several recent studies have documented activity of the CRISPR/Cas9 system at off-target loci in human cells1319. Public tools are available to assess and prioritize potential CRISPR/Cas9 activity at off-target loci based on predicted positional bias of a given mismatch in the sgRNA protospacer sequence and the total number of mismatches to the intended target site14. We used this public webserver to predict the most likely off-target sites for the sgRNAs used to correct the dystrophin gene in this study (Supplementary Table 2). The top ten potential off-target sites (OT1 to OT10) were rst assessed by the Surveyor assay in treated HEK293T cells because off-target activity in this cell line is well-known to be readily detectable due to high levels of transfection and expression of Cas9 and sgRNAs1319 and thus it would serve as a more rigorous test of specicity. We tested HE293T cells treated with SpCas9 and individual sgRNA expression cassettes for CR1, CR3, CR5, CR6 or CR36, for a total of 5 on-target and 50 off-target sites. For CR1, CR3 and CR36, one of the ten predicted off-target loci had signicant levels of gene modication (Supplementary Table 2 and Supplementary Fig. 6), consistent with other studies investigating the specicity of the S. pyogenes CRISPR/Cas9 system1318. Interestingly, the CR3 off-target sequence had substantial homology and similar modication frequencies to the intended on-target (9.3% at OT1 versus 13.3% at intended site; Supplementary Table 2 and Supplementary Fig. 6). Notably, CR3-OT1 was the only one of these three off-target sites to show signicant levels of activity in the sorted DMD myoblasts by the Surveyor assay (Fig. 7b).

6 NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244 ARTICLE

a

-sgRNA

CR1

CR2

CR5

CR1 + CR5

CR2 + CR5

Unmodified

Deletions

% Deleted:

13.6 10.5

b

CR1

CR5

CR5

Intron 50

Intron 50

PAM PAM

Intron 51

Intron 51

Intron 51 CR1/5-treated

CR2/5-treated

genomic DNA

genomic DNA

Intron 50

Intron 50

(2)

(2)

CR2PAM PAM

Intron 51

c d

-sgRNA

CR1

CR2

CR5

CR1 + CR5

CR2 + CR5

-sgRNA

CR1

CR2

CR5

CR1 + CR5

CR2 + CR5

Exon 47

Exon 52

4850 4851

* *

Dystrophin

GAPDH

Figure 4 | Deletion of exon 51 from the human genome using multiplex CRISPR/Cas9 gene editing. (a) End-point genomic PCR across the exon51 locus in human DMD myoblasts with a deletion of exons 4850. The top arrow indicates the expected position of full-length PCR amplicons and the two lower arrows indicate the expected position of PCR amplicons with deletions caused by the indicated sgRNA combinations. (b) PCR products from a were cloned and individual clones were sequenced to determine insertions and deletions present at the targeted locus. The top row shows the wild-type unmodied sequence and the triangles indicate SpCas9 cleavage sites. At the right are representative chromatograms showing the sequences of the expected deletion junctions. (c) End-point reverse transcriptionPCR analysis of dystrophin mRNA transcripts in CRISPR/Cas9-modied human D4850

DMD myoblasts treated with the indicated sgRNAs. A representative chromatogram of the expected deletion PCR product is shown at the right. Asterisk: band resulting from hybridization of the deletion product strand to the unmodied strand. (d) Rescue of dystrophin protein expression by CRISPR/Cas9 genome editing was assessed by western blot for the dystrophin protein with GAPDH as a loading control. The arrow indicates the expected restored dystrophin protein band.

Deep sequencing of on-target and off-target sites was also used to quantify gene-editing activity with greater sensitivity. Similar to the strategy with the Surveyor assay, we rst screened the 5 on-target and all 50 off-target sites in HEK293T cells (Table 2 and Supplementary Table 3). Notably, the indel frequencies measured by deep sequencing were signicantly greater than those determined by the Surveyor assay (Table 1), suggesting the Surveyor assay may underestimate overall gene editing activity. The detection limit for activity by deep sequencing was B0.010.1% indel-modied alleles, in contrast to B1% by the

Surveyor assay. Off-target sites were classied by greater than0.2% indels or a tenfold increase in indels compared with control. By this independent method, the same three off-target sites were identied by deep sequencing as by the Surveyor assay. These three off-target sites and their corresponding on-target sites were also assayed in the sorted, treated DMD myoblasts and only one of the sites showed detectable levels of activity, at a frequency approximately one-third of on-target activity (Table 2 and Supplementary Table 4). Thus, our selected sgRNAs have relatively favourable specicities, however, we cannot rule out activity at predicted off-target loci that may exist below the

sensitivity of these assays or off-target activity at other loci that were not assessed here. Notably, these sequencing results also showed a strong preference for single base-pair insertions, consistent with previous studies41, as well as enrichment for specic indels that are likely mediated by microhomology at the DNA break42 and could be used to favour conversion to specic reading frames (Supplementary Figs 7 and 8).

Nuclease activity at off-target sites may cause unintended chromosomal rearrangements by distal re-ligation between cleaved target and off-target loci on distinct chromosomes19,43. This presents a signicant concern for deletion-based gene correction strategies because of the increased potential for off-target activity by using two or more nucleases, such as in multiplex CRISPR/Cas9 gene editing. We probed for potential translocations using a highly sensitive nested genomic PCR assay to detect translocations at the validated off-target loci (Table 2 and Supplementary Table 2) during both single and multiplex CRISPR/Cas9-editing strategies. Using this assay, we readily detected translocations between on-target and off-target sites in the model HEK293T cell line that also shows higher levels of off-target activity (Fig. 7c and Supplementary Fig. 9a,b). Sanger

NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 7

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244

a

b

CR6 CR36 Intron 55

PAM

Intron 44 PAM

-sgRNA

-sgRNA

CR36

CR6 + CR36

CR6 + CR36

CR6

(5)

HEK293T

hDMD Mbs

Intron 44

c d

CR6

CR36

-sgRNA

CR6

CR36

CR6 + CR36

Exon 44

Exon 56

4850

4555

Dystrophin

GAPDH

Intron 55

Figure 5 | Deletion of the exon 4555 region in human DMD myoblasts by multiplex CRISPR/Cas9 gene editing. (a) End-point PCR of genomic DNA to detect deletion of the region between intron 44 and intron 55 after treating HEK293T cells or DMD myoblasts with the indicated sgRNAs.(b) Individual clones of PCR products of the expected size for the deletions from DMD myoblasts in a were analysed by Sanger sequencing to determine the sequences of genomic deletions present at the targeted locus. Below is a representative chromatogram showing the sequence of the expected deletion junctions. (c) End-point reverse transcriptionPCR analysis of dystrophin mRNA transcripts in CRISPR/Cas9-modied human D4850 DMD myoblasts treated with the indicated sgRNAs. A representative chromatogram of the expected deletion PCR product is shown at the right. (d) Analysis of restored dystrophin protein expression by western blot following electroporation of DMD myoblasts with sgRNAs targeted to intron 44 and/or intron 55.

Dystrophin

Spectrin

DMD untreated

100 m 100 m

100 m

Sorted population

CR1/5 treated

100 m

Figure 6 | Expression of restored human dystrophin in vivo following cell transplantation. Human D4850 DMD myoblasts were treated with

SpCas9, CR1 and CR5 to delete exon 51 and sorted for GFP expression as shown in Fig. 2. These sorted cells and untreated control cells were injected into the hind limbs of immunodecient mice and assessed for human-specic protein expression in muscle bres after 4 weeks post transplantation. Cryosections were stained with anti-human spectrin, which is expressed by both uncorrected and corrected myoblasts that have fused into mouse myobers, or anti-human dystrophin antibodies as indicated. White arrows indicate muscle bres positive for human dystrophin.

Scale bars indicate 100 mm. ND, not determined.

sequencing of the PCR amplicons conrmed the identity of the predicted translocation event for each primer pair (Supplementary Figs 10 and 11). A subset of the translocations detected in the HEK293T cells were also detectable by nested PCR in the sorted hDMD myoblasts, although the signal was considerably weaker and we were unable to conrm the sequence identity due to low yield of product (Fig. 7d and Supplementary Fig. 9a,c). Notably, we did not detect translocations using this assay in HEK293T cells or sorted DMD myoblasts treated with CR36 or CR6/CR36, respectively (Supplementary Fig. 9), that had

low levels of off-target activity at CR36-OT3 only in HEK293T cells (Table 2 and Supplementary Table 2). These results underscore the importance of selecting highly specic sgRNAs, particularly for multiplex editing applications, and show that this approach can benet from ongoing efforts to improve the specicity of the CRISPR/Cas9 system4,15,18,19,4446. However, taken together with the cytotoxicity and stable gene-editing activities observed in Table 1, these data suggest that the selected sgRNAs are able to correct the dystrophin gene without signicant toxicity and with only a single strongly predicted off-target site with detectable levels of activity.

DiscussionGenome editing is a powerful tool for correcting genetic disease and the recent development of the CRISPR/Cas9 system is dramatically accelerating progress in this area. A recent study showed the correction of the dystrophic phenotype in a mouse model of DMD by gene editing with CRISPR in single-cell mouse embryos47. Here we demonstrate the use of CRISPR to correct mutations to the human dystrophin gene in DMD patient cells. Many gene- and cell-based therapies for DMD are in preclinical development and clinical trials, and genome editing methods are compatible with many of these approaches. For example, genome editing may be combined with patient-specic cell-based therapies for DMD. Although transplantation of skeletal myoblasts, as used in this study, has been generally unsuccessful in clinical trials because of poor cell survival, engraftment, distribution and reconstitution of the progenitor cell population, several other muscle progenitor cells and pluripotent cell sources are currently under investigation and represent a potential path for this approach moving forward3440. In fact, the CRISPR/Cas9 system has been previously demonstrated to function in human stem cells9,10,48 and other human cell lines58, as well as human skeletal myoblasts in this study. Gene editing with CRISPR/Cas9 has also recently been achieved in vivo following plasmid delivery to the mouse liver49 and AAV delivery to the brain50. Furthermore, genome editing of autologous primary T cells with ZFNs is currently in clinical trials (NCT01252641, NCT00842634 and NCT01044654)51.

8 NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244 ARTICLE

a

b

100 CR3TGT sgRNA:

% Indels:

CR1 CR36

TGT TGT

OT1 OT1 OT3

Survival (%)

80

60

40

20

0

Mock

I-Scel

+ + + + + +

** *

*

* **

*

GZF3

-sgRNA

CR1

CR2

CR3

CR5

CR6

CR36

4.9 1.8 3.4 ND 1.0 ND

c

CR1

CR1

P1

P5

Primers: P1/P5 P1/P6 P2/P5 P2/P6

PCR primers: P2/P5 Sequencing primer: P2

489

432 401 339

Exp. size (bp): sgRNA:

+ + + +

chrX/intron 50

chr16

CR1-OT1

500 bp

* * * *

Primers:

sgRNA: + + + + + + + +

d

P6

P3

P1/P5

489 432 401 339 389 311 468 381

P1/P6 P2/P6 P3/P5 P3/P6 P4/P5 P4/P6

P2/P5

Exp. size (bp):

P1 CR5

Chr. X

P2

P2

Chr. X

Chr. 16

P4

500 bp * * * *

Figure 7 | Evaluation of CRISPR/Cas9 toxicity and off-target effects for deletion of human exon 51. (a) Results of a cytotoxicity assay in HEK293T cells treated with human-optimized SpCas9 and the indicated sgRNA constructs. Cytotoxicity is based on survival of GFP-positive cells that are co-transfected with the indicated nuclease. I-SceI is a well-characterized non-toxic meganuclease and GZF3 is a known toxic zinc nger nuclease. n 3 independent

transfections (mean s.e.m.). (b) Surveyor analysis at off-target sites in sorted hDMD cells treated with expression cassettes encoding SpCas9 and the

indicated sgRNAs. These three off-target sites tested in hDMD cells were identied from a panel of 50 predicted sites tested in HEK293T cells (Supplementary Fig. 6 and Supplementary Table 2). ND, not determined; OT, off-target locus; TGT, on-target locus for indicated sgRNA. (c,d) End-point nested PCR to detect chromosomal translocations in (c) HEK293T cells treated with Cas9 and CR1 or (d) sorted hDMD cells treated with Cas9, CR1 and CR5. The schematic depicts the relative location of nested primer pairs customized for each translocation event. The expected size of each band was estimated based on the primer size and the location of the predicted sgRNA cut site at each locus. Asterisks indicate bands detected at the expected size. The identities of the bands in c were veried by Sanger sequencing from each end (Supplementary Fig. 11). A representative chromatogram for the P2/P5 translocation in HEK293T cells is shown.

Importantly, gene editing with CRISPR/Cas9 did not abolish the myogenic capacity of the transfected myoblasts, as demonstrated by robust dystrophin expression in vitro and in vivo after transplantation into immunodecient mice. Thus, this strategy should be compatible with cell-based therapies for DMD that are under development3440, including induced pluripotent stem cells3739 that will likely provide a means to enhance the low engraftment efciencies observed here with extensively cultured myoblasts (Fig. 6). In addition, direct transfection of the sgRNA and Cas9 mRNA or protein to cultured cells, in contrast to the plasmid-based delivery method used here, may increase specicity and safety by reducing the duration of Cas9 expression and eliminating the possibility of random plasmid integration52,53. Delivery of Cas9 and the sgRNAs to cultured cells with non-integrating viral vectors, rather than electroporation, may also improve cell health, engraftment efciency and in vivo dystrophin expression54,55. Alternatively, delivery of the CRISPR/Cas9 system directly to skeletal and/or cardiac muscle by viral, plasmid or RNA delivery vectors is a promising strategy for in vivo genome editing and translation of this approach35,49,5659. The large size of S. pyogenes Cas9 gene (B4.2 kb) presents a challenge to its use in size-restricted adeno-associated viral vectors with typical constitutive or muscle-specic promoters. However, Cas9 genes from other species12,60, such as N. meningitidis and S. thermophilus, are small enough to

efciently package both Cas9 and sgRNA expression cassettes into single AAV vectors for in vivo gene-editing applications, as has been done with the smaller ZFNs61,62.

The S. pyogenes CRISPR/Cas9 system enabled efcient modication of nearly 90% of tested targets, consistent with other reports of robust activity of this system at diverse loci. The robustness and versatility of this technology is a signicant advancement towards at-will creation of patient-specic gene editing. However, further enhancing gene editing frequencies and minimizing off-target activity remain major challenges to customized clinical applications. Notably, recent studies have shown that low levels of dystrophin, including as little as 4% of wild-type expression, are sufcient to improve survival, motor function and cardiac function in a mouse model6365. Therefore, the levels of gene editing reported in this study and others for primary human cells, often in the range of 510% edited alleles, may be sufcient for therapeutic benet for this disease, particularly if using efcient in vivo delivery methods57,58 or a robust progenitor cell source3440 for an ex vivo strategy. Regardless, we suspect that optimized delivery and expression vehicles, as well as ongoing improvements to the CRISPR/Cas9 system4, will lead to signicant enhancements of gene correction levels.

The use of multiplexing with CRISPR/Cas9 to delete exons also presents a unique set of opportunities and challenges. This is the

NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 9

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244

Table 2 | On-target and off-target gene-editing activity by deep sequencing.

HEK293T HEK293T HEK293T HEK293T HEK293T CR1 CR3 CR5 CR6 CR36 %

Indels

Treated/ untreated

%

Indels

Treated/ untreated

%

Indels

Treated/ untreated

ON 43.347 833.596 47.12 889.13 49.82 503.20 41.74 1391.40 46.65 1227.66 OT1 2.714 100.519 14.21 1092.92 0.056 1.24 0.009 1.29 0.047 1.12 OT2 0.203 9.227 0.018 1.20 0.007 0.70 0.022 1.38 NA NA OT3 0.020 1.250 0.026 0.84 0.007 1.75 0.024 1.04 1.179 30.23 OT4 0.051 0.927 0.013 0.81 0.016 0.89 0.014 0.70 0.174 1.09 OT5 0.175 6.250 0.012 0.44 0.01 0.71 NA NA 0.138 1.12 OT6 0.057 1.295 0.038 1.81 0.064 3.05 0.141 1.10 0.062 1.19 OT7 0.117 2.294 0.032 1.07 0.019 1.00 0.116 1.20 0.02 0.87 OT8 0.068 1.478 0.011 0.58 0.127 1.40 0.008 1.60 0.215 0.67 OT9 0.089 0.967 0.014 0.70 0.018 1.06 0.042 1.24 0.014 0.78 OT10 0.030 0.612 0.01 0.53 0.023 1.10 0.013 0.68 0.21 6.77

DMD myoblasts% Indels Treated/untreated CR1-ON 11.26 261.84

CR1-OT1 0.14 3.00 CR3-ON 6.91 135.43 CR3-OT1 2.66 221.67 CR36-ON 3.23 95.12 CR36-OT3 0.04 0.86

DMD, Duchenne muscular dystrophy; NA, not applicable.

HEK293Ts were transfected with constructs encoding human codon-optimized SpCas9 and the indicated sgRNA. The frequency of indel formation at each target site (ON, bold) and the top ten predicted off-target sites (OT110) was determined by deep sequencing. Two off-target sites (CR6-OT5 and CR36-OT2) did not produce reads that met ltering criteria. Off-target sites with activity greater than0.2% or greater than tenfold increase compared with untreated controls are highlighted in red. These off-target sites and their corresponding on-target site were also assessed in the treated, sorted DMD myoblasts (bottom).

Treated/ untreated

%

Indels

Treated/ untreated

%

Indels

rst study to report the deletion of complete exons from the genome to restore dystrophin expression, in contrast to previous genome editing efforts that restored the reading frame of the dystrophin gene with small indels generated by NHEJ-based DNA repair following the action of a single nuclease28,30, inserted missing exons into the gene by homologous recombination29 or inserted a microdystrophin expression cassette into a safe harbour site by homologous recombination55. A primary benet of this new approach is that NHEJ is an active DNA repair mechanism in all cell types, whereas homologous recombination is typically only active in mitotic cells. As a result, approaches dependent on homologous recombination may not be applicable to directly modifying skeletal muscle bres, although recent evidence of genome editing in the adult liver provides some support that this may be possible49,61. Another advantage of the current method, relative to the indel-based approach, is that the protein product of the edited gene is predictable and already characterized in Becker muscular dystrophy patients with the naturally occurring deletion. This is in contrast to the random indels created by intraexonic action of a single nuclease that will lead to the creation of novel epitopes from each DNA repair event. Furthermore, the product resulting from the exon deletions will lead to restored dystrophin for every successful gene editing event, whereas modifying the gene with random indels within exons will only restore the reading frame in the one-third of editing events that leads to the correct reading frame. Additional studies on enhancing and/or biasing gene repair towards distal chromosomal re-ligation will be helpful, particularly for long deletions such as the excision of the 336 kb containing exons 4555 demonstrated in this study. Moreover, the introduction of multiple double-strand breaks increases the requirement for stringent specicity of the gene editing reagents in order to minimize the opportunity for the unintended chromosomal rearrangements that were readily detected in the treated

HEK293T cells (Fig. 7c,d and Supplementary Figs 911), but were less evident in the more relevant DMD myoblasts.

We observed that all of the sgRNAs tested were not associated with signicant cytotoxic effects in human cells. We identied three potential off-target sites out of 50 total tested sites for the ve sgRNAs used here to restore dystrophin expression (Table 2 and Supplementary Table 2). Furthermore, chromosomal trans-locations between the intended on-target sites and these off-target sites were detectable by highly sensitive nested PCR assays in HEK293T cells expressing high levels of Cas9 and sgRNAs. This is consistent with previous work showing translocations with two nucleases intentionally targeted to distinct chromosomes43. These results are also consistent with other studies that have characterized off-target activity and translocations generated by SpCas9 (refs 1319) and compare favourably with other gene-editing systems, including ZFNs and TALENs, which are sometimes cytotoxic66 and are known to also act at off-target sites67,68. Notably, the off-target activity and translocations identied in HEK293T cells, which is an immortalized and aneuploid cell line that expresses very high levels of Cas9 and sgRNA, did not occur at as high a level and in some cases were undetectable in the DMD myoblasts even though they were sorted for Cas9 expression (Table 2). This corroborates previous studies showing that lower levels of Cas9 and sgRNA can reduce off-target effects14,16. Importantly, this level of specicity may be tolerable given the severity of DMD, the lack of an apparent cytotoxic effect in human cells, and the absence of adverse events in ongoing clinical trials for genome editing with ZFNs (NCT01252641, NCT00842634 and NCT01044654)51. Furthermore, improvements to the specicity of the CRISPR/ Cas9 system, including the use of dual nickases, truncated sgRNAs, fusions to a FokI endonuclease domain, and careful sgRNA selection methods15,18,19,4446, are also under development. CRISPR/Cas9 specicity may be further enhanced

10 NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244 ARTICLE

by rational design or directed evolution of the S. pyogenes sgRNA and/or Cas9 nuclease, or characterization of novel Cas9 nucleases from other species with more stringent sgRNA or PAM requirements12,60.

This study demonstrates that CRISPR/Cas9 genome editing is a robust, easily programmable method to rapidly generate targeted frameshifts or genomic deletions that can address a variety of dystrophin mutations. Importantly, this method can reproducibly modify the dystrophin gene by deleting exons 4555, thereby addressing more than half of DMD patient deletions with a single genome-editing strategy. Further advancements in the delivery, specicity and efciency of these reagents will enhance the utility of this method for correcting the dystrophin gene and creating other custom genetic modications. Thus, CRISPR/Cas9 genome editing offers an exciting new avenue for gene therapies to treat Duchenne muscular dystrophy and other hereditary disorders.

Methods

Plasmid constructs. The expression cassettes for the S. pyogenes sgRNA and human codon-optimized Cas9 (hCas9) nuclease were used as previously described69. In order to create a uorescent reporter system to enrich CRISPR/ Cas9-modied cells, a GeneBlock (IDT) was synthesized containing a portion of the 30 end of the Cas9-coding sequence fused to a T2A skipping peptide immediately upstream of a multiple cloning site and subsequently cloned into the hCas9 expression vector. An eGFP reporter gene was then cloned into the T2A vector to allow co-translation of Cas9 and eGFP proteins from the same expression vector (hCas9-T2A-GFP, Supplementary Fig. 12).

Cell culture and transfection. HEK293T cells were obtained from the American Tissue Collection Center through the Duke Cell Culture Facility and were maintained in DMEM supplemented with 10% fetal bovine calf serum and1% penicillin/streptomycin. Immortalized myoblasts70 from a DMD patient harbouring a deletion of exons 4850 (D4850) in the dystrophin gene were maintained in skeletal muscle media (PromoCell) supplemented with 20% fetal bovine calf serum (Sigma), 50 mg ml 1 fetuin, 10 ng ml 1 human epidermal growth factor (Sigma), 1 ng ml 1 human basic broblast growth factor (Sigma), 10 mg ml 1 human insulin (Sigma), 1% GlutaMAX (Invitrogen) and 1% penicillin/

streptomycin (Invitrogen). All cell lines were maintained at 37 C and 5% CO2.

HEK293T cells were transfected with Lipofectamine 2000 (Invitrogen) with 400 ng of each expression vector according to the manufacturers protocol in 24-well plates. Immortalized myoblasts were transfected with 5 mg of each expression vector by electroporation using the Gene Pulser XCell (Bio-Rad) with PBS as an electroporation buffer using optimized conditions28. Transfection efciencies were measured by delivering an eGFP expression plasmid (pmaxGFP, Clontech) and using ow cytometry. These efciencies were routinely Z95% for HEK293T and Z70% for the immortalized myoblasts.

Surveyor assay for endogenous gene modication. CRISPR/Cas9-induced lesions at the endogenous target site were quantied using the Surveyor nuclease assay, which detects mutations characteristic of nuclease-mediated NHEJ. After transfection, cells were incubated for 3 or 10 days at 37 C and genomic DNA was extracted using the DNeasy Blood and Tissue kit (QIAGEN). The target locus was amplied by 35 cycles of PCR with the AccuPrime High Fidelity PCR kit (Invitrogen) using primers specic to each locus (Supplementary Table 5). The resulting PCR products were randomly melted and reannealed in a thermal cycler with the programme: 95 C for 240 s, followed by 85 C for 60 s, 75 C for 60 s, 65 C for 60 s, 55 C for 60 s, 45 C for 60 s, 35 C for 60 s, and 25 C for 60 s with a

0.3 C s 1 rate between steps. Following reannealing, 8 ml of PCR product was mixed with 1 ml of Surveyor Nuclease S and 1 ml of Enhancer S (Transgenomic) and incubated at 42 C for 1 h. After incubation, 6 ml of digestion product was loaded onto a 10% TBE polyacrylamide gel and run at 200 V for 30 min. The gels were stained with ethidium bromide and quantied by densitometry using the ImageLab software suite (Bio-Rad) as previously described.

Fluorescence-activated cell sorting of myoblasts. DMD myoblasts were electroporated with 5 mm each of hCas9-T2A-GFP and sgRNA expression vectors and incubated at 37 C and 5% CO2. Three days after electroporation, cells were trypsinized and collected for FACS sorting using a FACSvantage II sorting machine. GFP-positive cells were collected and expanded for analysis.

PCR-based assay to detect genomic deletions. The exon 51 or exon 4555 loci were amplied from genomic DNA by PCR (Invitrogen AccuPrime High Fidelity PCR kit) using primers anking each locus. The anking primers were CelI-CR1/ 2-F and CelI-CR5-R for exon 51 or CelI-CR6-F and CelI-CR36-R for exon 4555

analysis (Supplementary Table 5). PCR products were separated on Tris-Acetate-EDTA (TAE)-agarose gels and stained with ethidium bromide for analysis.

PCR-based detection of translocations. Loci with predicted possible translocations were amplied by a two-step nested PCR (Invitrogen AccuPrime High Fidelity PCR kit for each step) of genomic DNA from cells transfected with Cas9 alone (control) or Cas9 with sgRNA. In the rst step, translocations that may occur at each on-target and off-target sgRNA target site were amplied by 35 cycles of PCR using combinations of Surveyor primers for each locus that were modied to include restriction sites to facilitate cloning and sequencing analysis (Supplementary Table 5). One microlitre of each PCR reaction was subjected to a second round of amplication by 35 rounds of PCR using nested primer sets custom designed for each individual predicted translocation (Supplementary Table 5). Each second nested PCR primer binds within the same approximate region within the primary amplicon; however, each pair was optimized using Primer3 online bioinformatics software to ensure specic detection of each translocation. PCR amplicons corresponding to the expected length of predicted translocations and only present in cells treated with sgRNA were puried (QIAGEN Gel Extraction kit) and analysed by Sanger sequencing.

mRNA analysis. Immortalized myoblasts were differentiated into myobers by replacing the growth medium with DMEM supplemented with 1% insulin-transferrin-selenium (Invitrogen #51500056) and 1% penicillin/streptomycin (Invitrogen #15140) for 5 days before the cells were trypsinized and collected. Total RNA was isolated from these cells using the RNeasy Plus Mini Kit (QIAGEN) according to the manufacturers instructions. RNA was reverse transcribed to cDNA using the VILO cDNA synthesis kit (Life Technologies #11754) and 1.5 mg of RNA for 2 h at 42 C according to the manufacturers instructions. The target loci were amplied by 35 cycles of PCR with the AccuPrime High Fidelity PCR kit (Invitrogen) using primers annealing to exons 44 and 52 to detect exon 51 deletion by CR1/5 or CR2/5 or primers annealing to exons 44 and 60 to detect exon 4555 deletion by CR6/36 (Supplementary Table 5). PCR products were run on TAE-agarose gels and stained with ethidium bromide for analysis. The resolved PCR bands were cloned and analysed by Sanger sequencing to verify the expected exon junctions.

Western blot analysis. To assess dystrophin protein expression, immortalized myoblasts were differentiated into myobers as above for 67 days. Cells were trypsinized, collected and lysed in RIPA buffer (Sigma) supplemented with a protease inhibitor cocktail (Sigma) and the total protein amount was quantied using the bicinchoninic acid assay according to the manufacturers instructions (Pierce). Samples were then mixed with NuPAGE loading buffer (Invitrogen) and 5% b-mercaptoethanol and heated to 85 C for 10 min. Twenty-ve micrograms of protein were separated on 412% NuPAGE Bis-Tris gels (Invitrogen) with MES buffer (Invitrogen). Proteins were transferred to nitrocellulose membranes for 12 h in 1 tris-glycine transfer buffer containing 10% methanol and 0.01% SDS.

The blot was then blocked for 1 h with 5% milk in Tris-buffered saline with tween-20 (TBS-T) at room temperature. Blots were probed with the following primary antibodies: MANDYS8 to detect dystrophin (1:1,000, Sigma D8168) or rabbit anti-GAPDH (1:5,000, Cell Signaling 2118S). Blots were then incubated with mouse or rabbit horseradish peroxidase-conjugated secondary antibodies (Santa Cruz) and visualized using the ChemiDoc chemilumescent system (Bio-Rad) and Western-C ECL substrate (Bio-Rad). All uncropped gel images are provided in Supplementary Fig. 13.

Transplantation into immunodecient mice. All animal experiments were conducted under protocols approved by the Duke Institutional Animal Care and Use Committee. Cells were trypsinized, collected and washed in 1 Hanks Balanced

Salt Solution (HBSS, Sigma). Two million cells were pelleted and resuspended in 5 ml 1 HBSS (Sigma) supplemented with cardiotoxin (Sigma #C9759)

immediately before injection. These cells were transplanted into the hindlimb tibialis anterior (TA) muscle of 8-week-old male NOD.SCID.gamma mice(Duke CCIF Breeding Core) by intramuscular injection. Four weeks after injection, mice were euthanized and the TA muscles were harvested.

Immunouorescence staining. Harvested TA muscles were incubated in 30% glycerol overnight at 4 C before mounting and freezing in Optimal Cutting Temperature compound. Serial 10-mm sections were obtained by cryosectioning of the embedded muscle tissue at 20 C. Cryosections were then washed in PBS to

remove the OCT compound and subsequently blocked for 3060 min at room temperature in PBS containing 10% heat-inactivated fetal bovine serum for spectrin detection or 5% heat-inactivated fetal bovine serum for dystrophin detection. Cryosections were incubated overnight at 4 C with the following primary antibodies that are specic to human epitopes only: anti-spectrin (1:20, Leica NCL-SPEC1) or anti-dystrophin (1:2, Leica NCL-DYS3). After primary staining, spectrin or dystrophin expression was detected using a tyramide-based immunouorescence signal amplication detection kit (Life Technologies, TSA Kit #22, catalogue #T-20932,). Briey, cryosections were incubated with 1:200 goat

NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 11

& 2015 Macmillan Publishers Limited. All rights reserved.

ARTICLE NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244

anti-mouse biotin-XX secondary (Life Technologies #B2763) in blocking buffer for 1 h at room temperature. The signal was then amplied using streptavidin-HRP conjugates (1:100, from TSA Kit) in blocking buffer for 1 h at room temperature. Finally, cryosections were incubated with tyramide-AlexaFluor488 conjugates (1:100, TSA kit) in manufacturer-provided amplication buffer for 10 min at room temperature. Stained cryosections were then mounted in ProLong AntiFade (Life Technologies #P36934) and visualized with conventional uorescence microscopy.

Deep sequencing. For each of the ve gRNAs, PCR was used to amplify both the target region as well as the 50 candidate off-target sites. PCR primers included Illumina TruSeq sequencing primer sequences on the 50 ends (Supplementary

Table 6). A second round of PCR was then used to add Illumina owcell binding sequences and experiment-specic indexes 50 of the primer sequence. The resulting products were pooled and multiplex sequenced with 250 bp paired-end reads on an Illumina MiSeq instrument.

Because PCR products averaged 106 bp in length, the 30 ends of paired reads overlapped, and that overlap was used to infer complete amplicon fragments via single ungapped alignment, parameterized to score each match as 5 and each mismatch as 4. Fragments were trimmed to remove Illumina adapter and primer

sequences by performing ungapped alignment (parameterized as above to allow mismatches) of paired primer sequences to fragments and retaining only sequences between primers. Trimmed fragments were aligned to the human reference genome using BLAT. Any fragment with a top-scoring alignment that was different than the expected PCR product was discarded from downstream analysis. Each of the remaining fragments was then aligned to the reference genome sequence ofthe expected PCR product using a global afne alignment with the following parameterization: match 5, mismatch 4, gap open 5, gap extend 2.

Alignments were then trimmed to the 20-bp anking the predicted site of Cas9 nuclease activity three base pairs 50 of the PAM; for some targets, the 20-bp window extended beyond the end of the alignment and was therefore truncated at the end of the alignment. Indel (insertion and/or deletion) statistics were gathered from these windows separately for each treatment/control following demultiplexing via sequenced barcodes, by counting gaps in the query and subject sequences of the resulting truncated alignments and tabulating numbers of fragments having any indels in these windows.

Cytotoxicity assay. To quantitatively assess potential sgRNA or SpCas9 nuclease-associated cytotoxicity, HEK293T cells were transfected with 10 ng of a GFP reporter and 100 ng SpCas9 expression vector and 100 ng sgRNA expression vector using Lipofectamine 2000 according to the manufacturers instructions (Invitrogen). The percentage of GFP-positive cells was assessed at 2 and 5 days by ow cytometry. The survival rate was calculated as the decrease in GFP-positive cells from days 2 to 5 and normalized to cells transfected with an empty nuclease expression vector as described.

References

1. Gaj, T., Gersbach, C. A. & Barbas, III C. F. ZFN, TALEN, and CRISPR/Cas-based methods for genome engineering. Trends Biotechnol. 31, 397405 (2013).

2. Joung, J. K. & Sander, J. D. TALENs: a widely applicable technology for targeted genome editing. Nat. Rev. Mol. Cell. Biol. 14, 4955 (2013).

3. Silva, G. et al. Meganucleases and other tools for targeted genome engineering: perspectives and challenges for gene therapy. Curr. Gene Ther. 11, 1127 (2011).

4. Sander, J. D. & Joung, J. K. CRISPR-Cas systems for editing, regulating and targeting genomes. Nat. Biotechnol. 32, 347355 (2014).

5. Mali, P. et al. RNA-guided human genome engineering via Cas9. Science 339, 823826 (2013).

6. Cong, L. et al. Multiplex genome engineering using CRISPR/Cas systems. Science 339, 819823 (2013).

7. Jinek, M. et al. RNA-programmed genome editing in human cells. Elife 2, e00471 (2013).

8. Cho, S. W., Kim, S., Kim, J. M. & Kim, J. S. Targeted genome engineering in human cells with the Cas9 RNA-guided endonuclease. Nat. Biotechnol. 31, 230232 (2013).

9. Hou, Z. et al. Efcient genome engineering in human pluripotent stem cells using Cas9 from Neisseria meningitidis. Proc. Natl Acad. Sci. USA 110, 1564415649 (2013).

10. Yang, L. et al. Optimization of scarless human stem cell genome editing. Nucleic Acids Res. 41, 90499061 (2013).

11. Wang, H. et al. One-step generation of mice carrying mutations in multiple genes by CRISPR/Cas-mediated genome engineering. Cell 153, 910918 (2013).

12. Esvelt, K. M. et al. Orthogonal Cas9 proteins for RNA-guided gene regulation and editing. Nat. Methods 10, 11161121 (2013).

13. Fu, Y. et al. High-frequency off-target mutagenesis induced by CRISPR-Cas nucleases in human cells. Nat. Biotechnol. 31, 822826 (2013).

14. Hsu, P. D. et al. DNA targeting specicity of RNA-guided Cas9 nucleases. Nat. Biotechnol. 31, 827832 (2013).

15. Mali, P. et al. CAS9 transcriptional activators for target specicity screening and paired nickases for cooperative genome engineering. Nat. Methods 31, 833838 (2013).

16. Pattanayak, V. et al. High-throughput proling of off-target DNA cleavage reveals RNA-programmed Cas9 nuclease specicity. Nat. Biotechnol. 31, 839843 (2013).

17. Cradick, T. J., Fine, E. J., Antico, C. J. & Bao, G. CRISPR/Cas9 systems targeting beta-globin and CCR5 genes have substantial off-target activity. Nucleic Acids Res. 41, 95849592 (2013).

18. Ran, F. A. et al. Double nicking by RNA-guided CRISPR Cas9 for enhanced genome editing specicity. Cell 154, 13801389 (2013).

19. Cho, S. W. et al. Analysis of off-target effects of CRISPR/Cas-derived RNA-guided endonucleases and nickases. Genome Res. 24, 132141 (2014).

20. Canver, M. C. et al. Characterization of genomic deletion efciency mediated by clustered regularly interspaced palindromic repeats (CRISPR)/Cas9 nuclease system in mammalian cells. J. Biol. Chem. 289, 2131221324 (2014).

21. Hoffman, E. P., Brown, Jr R. H. & Kunkel, L. M. Dystrophin: the protein product of the Duchenne muscular dystrophy locus. Cell 51, 919 (1987).

22. Lu, Q. L. et al. The status of exon skipping as a therapeutic approach to duchenne muscular dystrophy. Mol. Ther. 19, 915 (2011).

23. Aartsma-Rus, A. et al. Theoretic applicability of antisense-mediated exon skipping for Duchenne muscular dystrophy mutations. Hum. Mutat. 30, 293299 (2009).

24. Cirak, S. et al. Exon skipping and dystrophin restoration in patients with Duchenne muscular dystrophy after systemic phosphorodiamidate morpholino oligomer treatment: an open-label, phase 2, dose-escalation study. Lancet 378, 595605 (2011).

25. Goemans, N. M. et al. Systemic administration of PRO051 in Duchennes muscular dystrophy. New Engl. J. Med. 364, 15131522 (2011).

26. Aoki, Y. et al. Bodywide skipping of exons 45-55 in dystrophic mdx52 mice by systemic antisense delivery. Proc. Natl Acad. Sci. USA 109, 1376313768 (2012).

27. Aartsma-Rus, A. et al. Exploring the frontiers of therapeutic exon skipping for Duchenne muscular dystrophy by double targeting within one or multiple exons. Mol. Ther. 14, 401407 (2006).

28. Ousterout, D. G. et al. Reading frame correction by targeted genome editing restores dystrophin expression in cells from Duchenne muscular dystrophy patients. Mol. Ther. 21, 17181726 (2013).

29. Popplewell, L. et al. Gene correction of a duchenne muscular dystrophy mutation by meganuclease-enhanced exon knock-in. Hum. Gene Ther. 24, 692701 (2013).

30. Rousseau, J. et al. Endonucleases: tools to correct the dystrophin gene. J. Gene Med. 13, 522537 (2011).

31. Ding, Q. et al. Enhanced efciency of human pluripotent stem cell genome editing through replacing TALENs with CRISPRs. Cell Stem Cell 12, 393394 (2013).

32. Kim, H. et al. Surrogate reporters for enrichment of cells with nuclease-induced mutations. Nat. Methods 8, 941943 (2011).

33. Lee, H. J., Kim, E. & Kim, J. S. Targeted chromosomal deletions in human cells using zinc nger nucleases. Genome Res. 20, 8189 (2010).

34. Tedesco, F. S., Dellavalle, A., Diaz-Manera, J., Messina, G. & Cossu, G. Repairing skeletal muscle: regenerative potential of skeletal muscle stem cells.J. Clin. Invest. 120, 1119 (2010).35. Konieczny, P., Swiderski, K. & Chamberlain, J. S. Gene and cell-mediated therapies for muscular dystrophy. Muscle Nerve 47, 649663 (2013).

36. Kimura, E. et al. Cell-lineage regulated myogenesis for dystrophin replacement: a novel therapeutic approach for treatment of muscular dystrophy. Hum. Mol. Genet. 17, 25072517 (2008).

37. Darabi, R. et al. Human ES- and iPS-derived myogenic progenitors restore DYSTROPHIN and improve contractility upon transplantation in dystrophic mice. Cell Stem Cell 10, 610619 (2012).

38. Filareto, A. et al. An ex vivo gene therapy approach to treat muscular dystrophy using inducible pluripotent stem cells. Nat. Commun. 4, 1549 (2013).

39. Zhao, C. et al. Recombinase-mediated reprogramming and dystrophin gene addition in mdx mouse induced pluripotent stem cells. PLoS ONE 9, e96279 (2014).

40. Tedesco, F. S. et al. Stem cell-mediated transfer of a human articial chromosome ameliorates muscular dystrophy. Sci. Transl. Med. 3, 96ra78 (2011).

41. Lin, Y. et al. CRISPR/Cas9 systems have off-target activity with insertions or deletions between target DNA and guide RNA sequences. Nucleic Acids Res. 42, 74737485 (2014).

42. Bae, S., Kweon, J., Kim, H. S. & Kim, J. S. Microhomology-based choice of Cas9 nuclease target sites. Nat. Methods 11, 705706 (2014).

43. Brunet, E. et al. Chromosomal translocations induced at specied loci in human stem cells. Proc. Natl Acad. Sci. USA 106, 1062010625 (2009).

44. Fu, Y., Sander, J. D., Reyon, D., Cascio, V. M. & Joung, J. K. Improving CRISPR-Cas nuclease specicity using truncated guide RNAs. Nat. Biotechnol. 32, 279284 (2014).

12 NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications

& 2015 Macmillan Publishers Limited. All rights reserved.

NATURE COMMUNICATIONS | DOI: 10.1038/ncomms7244 ARTICLE

45. Tsai, S. Q. et al. Dimeric CRISPR RNA-guided FokI nucleases for highly specic genome editing. Nat. Biotechnol. 32, 569576 (2014).

46. Guilinger, J. P., Thompson, D. B. & Liu, D. R. Fusion of catalytically inactive Cas9 to FokI nuclease improves the specicity of genome modication. Nat. Biotechnol. 32, 577582 (2014).

47. Long, C. et al. Prevention of muscular dystrophy in mice by CRISPR/Cas9-mediated editing of germline DNA. Science 345, 11841188 (2014).

48. Schwank, G. et al. Functional repair of CFTR by CRISPR/Cas9 in intestinal stem cell organoids of cystic brosis patients. Cell Stem Cell 13, 653658 (2013).

49. Yin, H. et al. Genome editing with Cas9 in adult mice corrects a disease mutation and phenotype. Nat. Biotechnol. 32, 551553 (2014).

50. Swiech, L. et al. In vivo interrogation of gene function in the mammalian brain using CRISPR-Cas9. Nat. Biotechnol. 33, 102106 (2014).

51. Tebas, P. et al. Gene editing of CCR5 in autologous CD4 T cells of persons infected with HIV. N. Engl. J. Med. 370, 901910 (2014).

52. Kim, S., Kim, D., Cho, S. W., Kim, J. & Kim, J. S. Highly efcient RNA-guided genome editing in human cells via delivery of puried Cas9 ribonucleoproteins. Genome Res. 24, 10121019 (2014).

53. Ramakrishna, S. et al. Gene disruption by cell-penetrating peptide-mediated delivery of Cas9 protein and guide RNA. Genome Res. 24, 10201027 (2014).

54. Lombardo, A. et al. Gene editing in human stem cells using zinc nger nucleases and integrase-defective lentiviral vector delivery. Nat. Biotechnol. 25, 12981306 (2007).

55. Benabdallah, B. F. et al. Targeted Gene addition of microdystrophin in mice skeletal muscle via human myoblast transplantation. Mol. Ther. Nucleic Acids 2, e68 (2013).

56. Pichavant, C. et al. Current status of pharmaceutical and genetic therapeutic approaches to treat DMD. Mol. Ther. 19, 830840 (2011).

57. Gregorevic, P. et al. Systemic delivery of genes to striated muscles using adeno-associated viral vectors. Nat. Med. 10, 828834 (2004).

58. Wang, Z. et al. Adeno-associated virus serotype 8 efciently delivers genes to muscle and heart. Nat. Biotechnol. 23, 321328 (2005).

59. Zangi, L. et al. Modied mRNA directs the fate of heart progenitor cells and induces vascular regeneration after myocardial infarction. Nat. Biotechnol. 31, 898907 (2013).

60. Chylinski, K., Makarova, K. S., Charpentier, E. & Koonin, E. V. Classication and evolution of type II CRISPR-Cas systems. Nucleic Acids Res. 42, 60916105 (2014).

61. Anguela, X. M. et al. Robust ZFN-mediated genome editing in adult hemophilic mice. Blood 122, 32833287 (2013).

62. Li, H. et al. In vivo genome editing restores haemostasis in a mouse model of haemophilia. Nature 475, 217221 (2011).

63. van Putten, M. et al. Low dystrophin levels increase survival and improve muscle pathology and function in dystrophin/utrophin double-knockout mice. FASEB J. 27, 24842495 (2013).

64. Li, D., Yue, Y. & Duan, D. Marginal level dystrophin expression improves clinical outcome in a strain of dystrophin/utrophin double knockout mice. PLoS ONE 5, e15286 (2010).

65. van Putten, M. et al. Low dystrophin levels in heart can delay heart failure in mdx mice. J. Mol. Cell. Cardiol. 69C, 1723 (2014).

66. Cornu, T. I. et al. DNA-binding specicity is a major determinant of the activity and toxicity of zinc-nger nucleases. Mol. Ther. 16, 352358 (2008).

67. Perez, E. E. et al. Establishment of HIV-1 resistance in CD4 T cells by

genome editing using zinc-nger nucleases. Nat. Biotechnol. 26, 808816 (2008).68. Guilinger, J. P. et al. Broad specicity proling of TALENs results in engineered nucleases with improved DNA-cleavage specicity. Nat. Methods 11, 429435 (2014).

69. Perez-Pinera, P. et al. RNA-guided gene activation by CRISPR-Cas9-based transcription factors. Nat. Methods 10, 973976 (2013).

70. Mamchaoui, K. et al. Immortalized pathological human myoblasts: towards a universal tool for the study of neuromuscular disorders. Skelet. Muscle 1, 111 (2011).

Acknowledgements

We thank Pablo Perez-Pinera for contribution of reagents used in this study, Vincent Mouly for the immortalized DMD myoblasts and Farshid Guilak for access to histology equipment. This work was supported by the Muscular Dystrophy Association (MDA277360), a Duke-Coulter Translational Partnership Grant, a Hartwell Foundation Individual Biomedical Research Award, a Basil OConnor Starter Scholar Award from the March of Dimes Foundation, NIH UH3TR000505, NIH P30AR066527 and an NIH Directors New Innovator Award (DP2-OD008586) to C.A.G. D.G.O. was supported by an American Heart Association Mid-atlantic Afliate Predoctoral Fellowship. P.I.T. was supported by a National Science Foundation Graduate Research Fellowship and an American Heart Association Mid-atlantic Afliate Predoctoral Fellowship.

Author contributions

D.G.O. conducted all experiments. A.M.K. and P.I.T. assisted with animal surgeries. D.G.O., A.M.K., P.I.T., W.H.M., T.E.R. and C.A.G. contributed to experimental design and data analysis. D.G.O. and C.A.G. wrote the manuscript with editing and approval by all authors.

Additional information

Accession codes: Sequencing data have been deposited in the NCBI Sequence Read Archive (SRA) under accession code SRP051210.

Supplementary Information accompanies this paper at http://www.nature.com/naturecommunications

Web End =http://www.nature.com/ http://www.nature.com/naturecommunications

Web End =naturecommunications

Competing nancial interests: C.A.G. and D.G.O. have led patent applications related to genome editing for Duchenne muscular dystrophy. C.A.G. is a scientic advisor to Editas Medicine, a company engaged in therapeutic development of genome editing technologies.

Reprints and permission information is available online at http://npg.nature.com/reprintsandpermissions/

Web End =http://npg.nature.com/ http://npg.nature.com/reprintsandpermissions/

Web End =reprintsandpermissions/

How to cite this article: Ousterout, D. G. et al. Multiplex CRISPR/Cas9-based genome editing for correction of dystrophin mutations that cause Duchenne muscular dystrophy. Nat. Commun. 6:6244 doi: 10.1038/ncomms7244 (2015).

NATURE COMMUNICATIONS | 6:6244 | DOI: 10.1038/ncomms7244 | http://www.nature.com/naturecommunications

Web End =www.nature.com/naturecommunications 13

& 2015 Macmillan Publishers Limited. All rights reserved.