Genome editing has potential for the targeted correction of germline mutations. Here we describe the correction of the heterozygous MYBPC3 mutation in human preimplantation embryos with precise CRISPR–Cas9-based targeting accuracy and high homology-directed repair efficiency by activating an endogenous, germline-specific DNA repair response. Induced double-strand breaks (DSBs) at the mutant paternal allele were predominantly repaired using the homologous wild-type maternal gene instead of a synthetic DNA template. By modulating the cell cycle stage at which the DSB was induced, we were able to avoid mosaicism in cleaving embryos and achieve a high yield of homozygous embryos carrying the wild-type MYBPC3 gene without evidence of off-target mutations. The efficiency, accuracy and safety of the approach presented suggest that it has potential to be used for the correction of heritable mutations in human embryos by complementing preimplantation genetic diagnosis. However, much remains to be considered before clinical applications, including the reproducibility of the technique with other heterozygous mutations.
More than 10,000 monogenic inherited disorders have been identified, affecting millions of people worldwide. Among these are autosomal dominant mutations, where inheritance of a single copy of a defective gene can result in clinical symptoms. Genes in which dominant mutations manifest as late-onset adult disorders include BRCA1 and BRCA2, which are associated with a high risk of breast and ovarian cancers1, and MYBPC3, mutation of which causes hypertrophic cardiomyopathy (HCM)2. Because of their delayed manifestation, these mutations escape natural selection and are often transmitted to the next generation. Consequently, the frequency of some of these founder mutations in particular human populations is very high. For example, the MYBPC3 mutation is found at frequencies ranging from 2% to 8%3 in major Indian populations, and the estimated frequency of both BRCA1 and BRCA2 mutations among Ashkenazi Jews exceeds 2%4.
HCM is a myocardial disease characterized by left ventricular hypertrophy, myofibrillar disarray and myocardial stiffness; it has an estimated prevalence of 1:500 in adults5 and manifests clinically with heart failure. HCM is the commonest cause of sudden death in otherwise healthy young athletes. HCM, while not a uniformly fatal condition, has a tremendous impact on the lives of individuals, including physiological (heart failure and arrhythmias), psychological (limited activity and fear of sudden death), and genealogical concerns. MYBPC3 mutations account for approximately 40% of all genetic defects causing HCM and are also responsible for a large fraction of other inherited cardiomyopathies, including dilated cardiomyopathy and left ventricular non-compaction6. MYBPC3 encodes the thick filament-associated cardiac myosin-binding protein C (cMyBP-C), a signalling node in cardiac myocytes that contributes to the maintenance of sarcomeric structure and regulation of both contraction and relaxation2.
Current treatment options for HCM provide mostly symptomatic relief without addressing the genetic cause of the disease. Thus, the development of novel strategies to prevent germline transmission of founder mutations is desirable. One approach for preventing second-generation transmission is preimplantation genetic diagnosis (PGD) followed by selection of non-mutant embryos for transfer in the context of an in vitro fertilization (IVF) cycle. When only one parent carries a heterozygous mutation, 50% of the embryos should be mutation-free and available for transfer, while the remaining carrier embryos are discarded. Gene correction would rescue mutant embryos, increase the number of embryos available for transfer and ultimately improve pregnancy rates.
Recent developments in precise genome-editing techniques and their successful applications in animal models have provided an option for correcting human germline mutations. In particular, CRISPR–Cas9 is a versatile tool for recognizing specific genomic sequences and inducing DSBs7,8,9,10. DSBs are then resolved by endogenous DNA repair mechanisms, preferentially using a non-homologous end-joining (NHEJ) pathway. Obviously, NHEJ is inappropriate for gene correction applications because it introduces additional mutations in the form of insertions or deletions at the DSB site, commonly referred to as indels. In some cases, however, targeted cells activate an alternative DNA repair pathway called homology-directed repair (HDR) that rebuilds the DSB site using the non-mutant homologous chromosome or a supplied exogenous DNA molecule as a template, leading to correction of the mutant allele11,12. At present, CRISPR–Cas9 is predominantly used to introduce mutations and in the generation of gene knockouts using intrinsic NHEJ. Because HDR efficiency is relatively low8, applications of genome editing for gene therapy have been limited.
In early attempts, the introduction of genome editing constructs into one-cell embryos (zygotes), resulted in multicellular embryos or offspring with mosaic outcomes in individual cells13,14. Also, off-target mutations that could be introduced into the developing embryo remained an undesirable possibility.
We sought to investigate human gamete and embryo DNA repair mechanisms activated in response to CRISPR–Cas9-induced DSBs. In an effort to demonstrate the proof-of-principle that heterozygous gene mutations can be corrected in human gametes or early embryos, we focused on the MYBPC3 mutation that has been implicated in HCM. Although homozygous mutations with no PGD alternative would have been most desirable for gene correction, generating homozygous human embryos for research purposes is practically impossible. Homozygous MYBPC3 mutations in adults are extremely rare owing to the severity of the clinical symptoms and early onset of the disease. Therefore, we specifically targeted the heterozygous four-base-pair (bp) deletion in the MYBPC3 gene in human zygotes introduced by heterozygous, carrier sperm while oocytes obtained from healthy donors provided the wild-type allele. By accurate analysis of cleaving embryos at the single-cell level, we show high targeting efficiency and specificity of preselected CRISPR–Cas9 constructs. Moreover, DSBs in the mutant paternal MYBPC3 gene were preferentially repaired using the wild-type oocyte allele as a template, suggesting an alternative, germline-specific DNA repair response. Mechanisms responsible for mosaicism in embryos were also investigated and a proposed solution to minimize their occurrence developed—namely the co-injection of sperm and CRISPR–Cas9 components into metaphase II (MII) oocytes.
Subject with a heterozygous MYBPC3∆GAGT deletion
An adult male patient with well-documented familial HCM caused by a heterozygous dominant 4-bp GAGT deletion (g.9836_9839 del., NC_000011.10) in exon 16 of MYBPC3, currently managed with an implantable cardioverter defibrillator and antiarrhythmic medications, agreed to donate skin, blood and semen samples. Skin fibroblast cultures were expanded and used to generate heterozygous patient induced pluripotent stem cells (iPSCs) as described previously15. Two single-guide RNA (sgRNA)–Cas916,17,18 constructs were designed to target this specific MYBPC3∆GAGT deletion (Extended Data Fig. 1a, b) along with two exogenous single-stranded oligodeoxynucleotide (ssODN) templates encoding homology arms to the targeted region (Extended Data Table 1). To differentiate from the wild-type allele, two synonymous single-nucleotide substitutions were introduced into each ssODN template. In addition, ssODN-2 nucleotide substitutions provided an additional restriction enzyme (BstBI) recognition site (Extended Data Fig. 1a, b).
The efficacy and specificity of each construct were tested by transfecting patient iPSCs. Cells were electroporated together with ssODN, Cas9 and sgRNA expression plasmids and subcloned, and the targeted region for each clone was analysed by sequencing (Extended Data Fig. 1c). Of 61 iPSC clones transfected with CRISPR–Cas9-1, 44 (72.1%) were not targeted, as evidenced by the presence of both intact wild-type and intact mutant alleles. Among targeted clones, 10 of 17 (58.8%) were repaired by NHEJ and contained various indels adjacent to the mutation site (Extended Data Fig. 1d, e and Supplementary Table 1). The remaining seven clones were repaired by HDR using ssODN-1 as judged by the presence of the marker nucleotide substitutions. Thus, the total targeting efficiency for CRISPR–Cas9-1 was 27.9% (17/61). Among the targeted clones, only 41.2% (7/17) were repaired by HDR (Extended Data Fig. 1e). The targeting efficiency with CRISPR–Cas9-2 was 13.1% (23/175) and the HDR was considerably lower at 13% (3/23). Of note, among the three HDR-repaired iPSC clones, two were repaired using the ssODN-2 template while the third clone contained intact wild-type sequences in both alleles (Extended Data Fig. 1d, e and Supplementary Table 1), indicating HDR using the wild-type allele.
The wild-type allele in all iPSC clones analysed remained intact, demonstrating high fidelity of sgRNAs.
We also directly compared CRISPR–Cas9-1 and CRISPR–Cas9-2 in patient iPSCs transfected with preassembled Cas9 ribonucleoproteins (RNPs). Targeted deep sequencing demonstrated that CRISPR–Cas9-1 had higher HDR efficiency (Extended Data Fig. 1f). On-target mutations were not detected in wild-type embryonic stem (ES) cells (H9) carrying both wild-type MYBPC3 alleles, demonstrating high specificity of CRISPR–Cas9-1. On the basis of these outcomes, we selected CRISPR–Cas9-1 (hereafter referred to as CRISPR–Cas9), with higher efficiency of HDR-based gene correction, for subsequent studies.
HDR efficiency in heterozygous MYBPC3∆GAGT zygotes
We next evaluated targeting outcomes in human zygotes produced by fertilizing healthy donor oocytes with sperm from the patient heterozygous for the MYBPC3 mutation. As direct introduction of Cas9 protein is more efficient than using a plasmid9,19, recombinant Cas9 protein microinjection was adopted, injecting a mixture of sgRNA, Cas9 protein and ssODN DNA into the cytoplasm of pronuclear stage zygotes 18 h after fertilization. Injected zygotes and intact controls were cultured for 3 days before each embryonic blastomere was isolated and individually analysed by sequencing (Fig. 1). Cytoplasmic microinjection of CRISPR–Cas9 was confirmed visually (Supplementary Video 1) and shown to be efficient with a 97.1% (68/70) survival rate after injection and development rates comparable to controls (Extended Data Table 2).
Sequencing of 83 individual blastomeres collected from 19 control embryos revealed that 9 (47.4%) were homozygous wild type (MYBPC3WT/WT) and 10 (52.6%) were heterozygous, carrying the wild-type maternal and mutant paternal alleles (MYBPC3WT/∆GAGT) (Fig. 2a and Supplementary Table 2). This distribution was expected and confirms that the heterozygous patient sperm sample contained equal numbers of wild-type and mutant spermatozoa with similar motility and fertilization efficiency.
Among 54 CRISPR–Cas9-injected embryos, 36 (66.7%) were uniformly homozygous for the wild-type allele with each blastomere containing MYBPC3WT/WT, while 18 (33.3%) were uniform or mosaic heterozygous (Fig. 2a and Supplementary Table 3). In this group of 18, 5 embryos were uniformly heterozygous with each blastomere containing the intact wild-type and intact mutant allele (MYBPC3WT/∆GAGT) and 13 were mosaic, each containing blastomeres carrying more than one genotype (Fig. 2a and Supplementary Table 3). Each mosaic embryo contained at least one heterozygous blastomere with wild-type and either the intact ∆GAGT deletion or the ∆GAGT deletion plus additional indels, suggesting that these embryos originated from heterozygous (MYBPC3WT/∆GAGT) zygotes after fertilization by mutant sperm (Fig. 2b and Extended Data Table 3). Remarkably, the majority of the remaining sister blastomeres in eight mosaic embryos (numbers 1, 2, 4, 6, 7, 10, 11 and 12 in Fig. 2b) were homozygous for the wild-type allele (MYBPC3WT/WT). Overall, 52.2% (35/67) of individual blastomeres within mosaic embryos were homozygous MYBPC3WT/WT (Fig. 2c and Supplementary Table 3). As these embryos originated from MYBPC3WT/∆GAGT zygotes, individual blastomeres are likely to have repaired the MYBPC3∆GAGT deletion by HDR using the maternal wild-type allele as a template instead of the injected ssODNs. This conclusion was corroborated by the observation that correction occurred in blastomeres of mosaic embryos not injected with ssODNs (Fig. 2b). Among the other genotypes, four mosaic embryos (numbers 5, 8, 10 and 13 in Fig. 2b) contained blastomeres with an intact, mutant allele (MYBPC3WT/∆GAGT), albeit with most (29.9%) also containing additional small deletions (1–20 bp length, n = 16) or insertions (1 bp, n = 3) adjacent to the DSB site (MYBPC3WT/∆GAGT-indel), characteristic of NHEJ. One blastomere carried a 10-bp deletion and a 5-bp insertion (Extended Data Table 3 and Supplementary Table 3) and mosaic embryo #9 displayed four various NHEJ genotypes in its blastomeres, suggesting that targeting and NHEJ repair had occurred independently multiple times after the first zygotic division.
The overall targeting efficiency in human embryos was 72.2% (13/18), much higher than in iPSCs exposed to the same construct at 27.9% (17/61) (Fig. 2d, e), probably owing to more efficient delivery of the CRISPR–Cas9 constructs by zygote microinjection as compared to transfection in iPSCs. Even more remarkably, the majority of targeted blastomeres (63.6%, 35/55) resolved the DSBs by HDR using the wild-type allele, also markedly different from what was seen in iPSCs (Fig. 2d and Extended Data Fig. 2a). We did not find any evidence of HDR using exogenous ssODN, suggesting that HDR is guided exclusively by the wild-type maternal allele.
As HDR calculations are based on mosaic embryos only, it is probable that some targeted heterozygous (MYBPC3WT/∆GAGT) zygotes repaired the mutant allele in all blastomeres using the wild-type template (MYBPC3WT/WT). These HDR-repaired, uniform embryos would be indistinguishable from their wild-type homozygous counterparts, thereby increasing the portion of MYBPC3WT/WT embryos in the CRISPR–Cas9 injected group. Indeed, 66.7% (36/54) of injected embryos were homozygous WT/WT, a substantial increase over the WT/WT yield (47.4%, 9/19) in control non-injected embryos (Fig. 2f). Similar to observations in iPSCs, the wild-type allele was intact in all human embryos, corroborating the specificity of the sgRNA.
In summary, these results clearly demonstrate that CRISPR–Cas9 gene targeting in human zygotes is exceptionally efficient and that DSBs in the mutant paternal allele were predominantly repaired through HDR. Furthermore, HDR was exclusively directed by the homologous wild-type allele present on the maternal chromosome. This outcome suggests that human embryos employ different DNA repair mechanisms than do somatic or pluripotent cells, probably reflecting evolutionary requirements for stringent control over genome fidelity in the germline.
Targeting gametes eliminates mosaicism
Mosaicism in gene-targeted human embryos is unacceptable in clinical applications. The presence even of a single mutant blastomere within a mosaic embryo would make detection by PGD problematic; therefore, molecular mechanisms responsible for mosaicism were investigated. Analysis of targeting outcomes in the majority of mosaic, zygote-injected human embryos revealed only two different genotypes (MYBPC3WT/HDR and MYBPC3WT/∆GAGT-indel, or MYBPC3WT/HDR and MYBPC3WT/∆GAGT; Fig. 2b and Extended Data Table 3). Embryos #5 and #9 were the exceptions, containing three or more genotypes (Extended Data Table 3). These findings suggest that CRISPR–Cas9 targeted at least two mutant sperm alleles despite injection into the zygote. There are two different possible explanations for this outcome: 1) at the time of injection, a zygote had completed the S-phase of the cell cycle with DNA replication and already produced two mutant alleles20; 2) CRISPR–Cas9 remained active, continuing to target after zygotic division.
We hypothesized that both situations could be abrogated if CRISPR–Cas9 was co-injected with sperm into the M-phase oocyte during fertilization by intracytoplasmic sperm injection (ICSI), allowing genome editing to occur when the sperm undoubtedly still contains a single mutant copy21. In addition, the extended time of exposure to MII cytoplasm could allow CRISPR–Cas9 components to degrade before DNA replication resulted in two or more mutant alleles (Fig. 3a). Therefore, CRISPR–Cas9 was mixed with the sperm suspension and co-injected into 75 MII oocytes during ICSI (Supplementary Video 1); no difference was observed in the survival, fertilization and cleavage rates between CRISPR–Cas9-injected and intact control oocytes (Extended Data Table 2). At day 3 after fertilization, embryos were disaggregated and each individual blastomere was analysed as described above for S-phase-injected zygotes. Blastomeres from 16 of 58 (27.6%) M-phase-injected embryos were uniformly heterozygous, carrying an intact wild-type maternal allele along with NHEJ-repaired mutant paternal sequences carrying various indels (MYBPC3WT/∆GAGT-indel) (Fig. 3b and Supplementary Table 4). The remaining 42 (72.4%) were MYBPC3WT/WT. Of these, the vast majority (41/42) were uniformly homozygous embryos consisting of blastomeres carrying indistinguishable MYBPC3WT/WT alleles. Interestingly, the remaining embryo (M2-WT42) contained four blastomeres with MYBPC3WT/WT but HDR-repaired with ssODN, while the other three sister blastomeres were MYBPC3WT/WT, suggesting HDR using the maternal wild-type allele (Supplementary Table 4). No heterozygous blastomeres with intact mutant alleles (MYBPC3WT/∆GAGT) were detected, indicating 100% targeting efficiency in the M-phase-injected group compared to 72.2% efficiency in the S-phase-injected zygotes (Fig. 2d and Fig. 3b). More importantly, all sister blastomeres in all but one embryo carried identical genotypes, indicating a marked reduction in mosaicism in M-phase-injected embryos. The only mosaic embryo had all blastomeres repaired by HDR (either wild-type or ssODN as a template). Thus, this embryo, with every blastomere carrying repaired MYBPC3WT/WT, would be eligible for transfer.
The yield of MYBPC3WT/WT embryos (72.4%, 42/58) in the M-phase-injected group was significantly higher than in untreated controls (47.4%, 9/19) (Fig. 3c, P Fig. 3d). To rule out the possibility that the observed increase in WT/WT embryos in CRISPR–Cas9-injected zygotes and oocytes was due to allele drop-out during PCR and Sanger sequencing, we validated genotypes by independent on-target deep sequencing (Extended Data Fig. 2b). The estimated HDR-based repair and increase in WT/WT embryos for the S-phase-injected and M-phase-injected groups were 16.7% (9/54) and 22.4% (13/58), respectively (Fig. 3e). In summary, delivery of CRISPR–Cas9 into MII oocytes provides more efficient targeting than injection into zygotes while eliminating mosaicism.
Development and cytogenetics of repaired embryos
To examine the effect of gene correction on preimplantation development, we cultured CRISPR–Cas9-injected embryos to blastocysts. Similar to intact controls, 72.7% (16/22) of M-phase-injected embryos developed to the eight-cell stage and 50.0% (11/22) progressed to blastocysts (Student’s t-test, P > 0.05) (Fig. 4a, b). In an effort to provide additional insights into the developmental competence of gene-corrected blastocysts, and to obtain sufficient cellular material for detailed cytogenetic studies, we established six ES cell lines from CRISPR–Cas9-injected blastocysts and one from controls. On-target analysis revealed that four CRISPR–Cas9-treated ES cell lines (ES-WT1, ES-WT2, ES-WT3 and ES-WT4) and one control cell line (ES-C1) were MYBPC3WT/WT whereas the remaining two CRISPR–Cas9-injected cell lines (ES-Mut1 and ES-Mut2) were MYBPC3WT/∆GAGT-indel (Extended Data Fig. 3a). These results corroborate the exceptionally high targeting efficiency of CRISPR–Cas9 in M-phase-injected human embryos.
Cytogenetic G-banding analysis revealed that ES-WT1, ES-WT4, ES-Mut1 and ES-Mut2 carried normal diploid karyotypes with no evidence of detectable numerical or structural chromosomal rearrangements (Extended Data Fig. 3a). Notably, ES-WT2, ES-WT3 and the control line ES-C1 exhibited a pericentric inversion on chromosome 10 (Extended Data Fig. 3a–d). As both treated and control ES cells showed this chromosomal rearrangement, we reasoned that it was contributed by the sperm and could be inherited. Indeed, analysis of the patient’s skin fibroblast-derived iPSCs showed the same inversion, indicating that this inversion was balanced. In summary, CRISPR–Cas9-treated human embryos displayed normal development to blastocysts and ES cells without cytogenetic abnormalities.
Off-target consequences in repaired human embryos
Apart from the overall targeting and HDR efficacy and mosaicism, one safety concern regarding clinical application of gene correction in human embryos is that CRISPR–Cas9 can induce undesirable off-target mutations at genome regions that are highly homologous to the targeted sequence22,23,24. Therefore, we conducted a comprehensive, whole-genome sequencing (WGS) analysis of the patient’s genomic DNA by digested genome sequencing (Digenome-seq)25,26. Potential off-target sequences were identified by digestion of iPSC-derived, cell-free genomic DNA with CRISPR–Cas9 followed by WGS. Sequencing reads of CRISPR–Cas9-digested genomic DNA are vertically aligned at on-/off-target sites in IGV viewer25,27. By contrast, undigested genomic sites are aligned in a staggered manner at those loci. In addition, improved Digenome-seq provides DNA cleavage scores for potential off-target sites based on alignment patterns of sequence reads26. Digested iPSC DNA produced uniform cleavage patterns at both on-target and potential off-target sites (Extended Data Fig. 4a, b). We identified 16 potential off-target sites with a DNA cleavage score higher than 2.5 (ref. 26) (Extended Data Fig. 5a and Extended Data Table 4). Sequencing analysis of these 16 sites with Web Logo (http://weblogo.berkeley.edu/) confirmed that they are indeed highly homologous to the on-target MYBPC3 mutant allele28 (Extended Data Fig. 5b). Furthermore, we identified seven additional sites with DNA cleavage scores of 0.1 or greater and with ten or fewer nucleotides mismatched in the human genome (Extended Data Table 4). Next, we sequenced and analysed all these sites in each individual blastomere from two untreated control embryos (C2 and C10 from Supplementary Table 2), two mosaic S-phase-injected embryos (Mos1 and Mos7), one uniform, non-mosaic S-phase-injected embryo (WT15 from Supplementary Table 3) and two M-phase-injected embryos (M2-WT10 and M2-Mut7 from Supplementary Table 4) (Extended Data Table 5). All on-target indels in each blastomere were corroborated and results were identical to the Sanger sequencing results. In addition, indels were not detected in any blastomeres known to be carrying either intact WT/WT or WT/Mut alleles at the target site (Extended Data Fig. 5c and Supplementary Table 5). More importantly, indels were also not detected in 23 off-target loci examined in 28 screened blastomeres (Extended Data Fig. 5d). In selected blastomeres (Extended Data Fig. 5c), we extended off-target screening by WGS. Potential off-target sites were examined by comparing genomic variants found in intact, control embryos (C2.3, C10.2) with those in CRISPR–Cas9-injected embryos (Mos1.1, W15.4, Mos7.2, M2-WT10.1, and M2-Mut7.1) (Extended Data Table 6). After filtering out annotated variants in the dbSNP database, we found 19–71 potential off-target sites with indels in each blastomere obtained from CRISPR–Cas9-injected embryos (lanes 4 and 7 in Extended Data Table 6). All of these sites contained repeated sequences such as poly-A or poly-GT repeats (lanes 5 and 8), suggesting that indels found at these sites were caused by sequencing errors rather than Cas9-catalysed, off-target DNA cleavage. These WGS results support our Digenome-seq conclusions that gene correction did not induce any detectable off-target mutations in selected blastomeres (Extended Data Table 6).
We also investigated whether CRISPR–Cas9 targeting induced global off-target genetic variations and genome instability by performing whole-exome sequencing (WES) in CRISPR–Cas9-treated ES cells and compared the results to those of control ES cells and corresponding egg and sperm donor blood DNA. WES analysis revealed a large number of variants in all samples when compared to the hg19 reference genome (Supplementary Table 6). The majority of these variants were also present in egg or sperm donors and found in the dbSNP and 1000genomes databases. Some variants detected in ES cells showed decreased fractions matching the population hotspots, indicating the potential effect of experimental procedures, including embryo culture and ES cell derivation and culture (Extended Data Fig. 6a and Supplementary Table 6). Three treated ES cell lines and a control line (ES-Mut1, ES-WT1, ES-WT2 and ES-C1) showed similar statistics in all variant categories and were comparable to gamete donor profiles (egg donors 1 and 2, sperm donor). ES-WT3 exhibited an increase in variant numbers but this sample did not have a control sibling ES cell line for comparison (Extended Data Fig. 6b, c and Supplementary Table 6). Next, we investigated potential off-target effects in ES cells and identified a total of 685 potential off-target sites using full sensitive aligner Batmis (V3.00)29 (Supplementary Table 7). Variants that were also present in the gamete donors were filtered out as inherited (Supplementary Table 8). Notably, analysis of these sites did not reveal any variants. Taken together, these Digenome-seq, WGS and WES results demonstrate high on-targeting specificity of CRISPR–Cas9 in human embryos without any off-target effects.
DSBs induced by genome editing are primarily resolved via error-prone NHEJ, and such repair approaches are predominantly used to generate gene knockouts in cells and organisms17,30,31. By contrast, HDR, although occurring at a substantially lower efficiency, is necessary for gene correction, particularly when human germline gene therapy is envisioned. We show here that DSBs in human gametes and zygotes are preferentially resolved using an endogenous HDR mechanism, exclusively directed by the wild-type allele as a repair template. By contrast, HDR efficiency in iPSCs was significantly lower and primarily achieved through an exogenous DNA template. This striking difference implies that human gametes and embryos employ a different DNA damage response system, perhaps reflecting the evolutionary importance of maintaining germline genome integrity32. If, as seems likely, gametes and zygotes endure an increased number of DSBs during meiotic recombination and segregation, an efficient genome repair capacity would be critical33 and unique zygotic DNA repair machinery might rely entirely on maternal oocyte factors deposited and stored during maturation since zygotes are transcriptionally silent. Recent studies suggest that oocytes might employ an ataxia-telangiectasia mutated (ATM)-mediated DNA damage signalling (DDS) pathway that regulates repair of DSBs via a homologous recombination mechanism34. Thus, Cas9-induced DNA breaks are likely to attract the existing native oocyte machinery reserved for repair of meiotic recombination-induced DSBs. On the basis of our results, it may not even be necessary to provide exogenous oligo templates for gene correction in heterozygous human embryos.
The efficiency of CRISPR–Cas9 was evaluated recently in a mouse study involving a heterozygous dominant mutation in the Crygc gene responsible for an inherited form of cataracts. Some HDR-repaired events used sequences from the wild-type allele from the homologous chromosome12. However, HDR via the exogenous oligo template was more frequent, with three of four pups carrying corrected Crygc genes with a DNA sequence from the exogenous oligo and only one from the wild-type allele. In a study involving human heterozygous embryos, HDR was exclusively directed by the exogenous DNA template with no evidence of wild-type allele-based repair35. As these results were derived from bulk DNA from whole embryos rather than individual blastomeres, it is possible that instances of HDR via the wild-type allele were overlooked.
Despite remarkable targeting efficiency and high HDR frequency, some CRISPR–Cas9-treated human embryos demonstrated NHEJ-induced indels and thus would not be suitable for transfer. Therefore, genome editing approaches must be further optimized before clinical application of germline correction can be considered. Modifications in genome editing by inhibiting NHEJ mechanisms while enhancing HDR pathways have been reported36,37. Other approaches have been focused on manipulating the cell cycle or modifying the donor ssDNA design11,30. While some of these developments significantly improved HDR outcomes in the context of cultured cells, their relevance to embryonic gene correction remains unknown. In addition, supplementary exposure of human gametes or embryos to small molecules and/or inhibitors may adversely affect embryonic development.
Nonhuman primate studies have demonstrated that CRISPR–Cas9 injection into monkey zygotes can disrupt wild-type genes, with the resultant full-term offspring carrying the mutations and associated phenotypes38,39. Similar to the outcomes seen in mice and other animals, genome-edited human preimplantation embryos and newborn monkeys display mosaic targeting genotypes in their cells and tissues, suggesting that DSBs and subsequent repair do not occur at the single mutant allele stage35,40,41. As discussed above, mosaicism in gene-corrected human embryos would be difficult to detect and would, therefore, pose serious safety issues for possible clinical applications. Modifications involving shortening the half-life of Cas9 activity reduced, but did not completely eliminate, the manifestations of mosaicism in monkey embryos41, while introduction of CRISPR–Cas9 into early-stage zygotes reduced mosaicism in mouse embryos42.
Importantly, we demonstrate that the delivery of CRISPR–Cas9 into M-phase oocytes abolished mosaicism in cleaving embryos, suggesting that gene targeting and editing efficiencies are strongly associated with DNA synthesis and the cell cycle phase42. The choice of DSB repair using either NHEJ or HDR may be dependent on the cell cycle phase, with HDR restricted to late S and G2 phases when DNA replication is completed and sister chromatids are available as repair templates11. In particular, HDR mechanisms were downregulated at the M and early G1 phases, thus favouring NHEJ-induced genome editing43. However, we did not observe reduced HDR efficiency even when CRISPR–Cas9 was delivered into MII oocytes at the time of ICSI. One possible explanation is that the DNA repair response is different in germ cell meiotic M phase compared to mitotic M phase in cultured cells. Alternatively, the DSBs may have occurred at the M or G1 phase, while the HDR repair followed later at the S or G2 phase of the cell cycle.
Extensive reports on potential off-target DNA damage induced by CRISPR–Cas9 beyond the intended targeting region have been published. In particular, Cas9 overexpression via plasmid transfection and subsequent high enzyme concentrations increased off-site targeting9. In our studies with human oocytes and zygotes, we used purified recombinant Cas9 protein instead of plasmid, which may enhance the specificity while shortening enzymatic exposure time, thereby diminishing off-site targeting. Screening by Digenome-seq, WGS and WES did not detect any off-target mutations in multiple individual blastomeres from CRISPR–Cas9-injected human embryos or ES cells. In addition, WES demonstrated maintenance of genome stability after gene correction. CRISPR–Cas9 targeting seen in our investigation was very accurate, providing assurance regarding safety concerns.
Other issues that need to be addressed include the applicability of this gene correction approach to other heterozygous mutations in humans. It may be challenging to repair homozygous mutations in human embryos when both alleles are mutant and wild-type allele-based HDR mechanisms cannot be employed.
As indicated, PGD may be a viable option for heterozygous couples at risk of producing affected offspring. In cases when only one parent carries a heterozygous mutation, 50% of embryos should be mutant. In contrast, targeted gene correction can potentially rescue a substantial portion of mutant human embryos, thus increasing the number of embryos available for transfer.
Regulations for research on human gametes and embryos
The regulatory framework surrounding the use of human gametes and embryos for this research was based on the guidelines set by the Oregon Health & Science University (OHSU) Stem Cell Research Oversight Committee (OSCRO). In 2008, OSCRO established policy and procedural guidelines formally defining the use of human embryos and their derivatives at OHSU, informed by the National Academy of Sciences’ Guidelines. These policies and guidelines permitted the procurement of gametes and embryos for research purposes, the creation of human embryos specifically for research, genetic manipulation of human gametes and embryos, creation of human embryonic stem cell lines and molecular analyses. Together, OSCRO and the OHSU Institutional Review Board (IRB) worked concurrently to review and monitor applications for research studies involving human embryos at OHSU.
Human embryo and embryonic stem cell research policies and principles at OHSU were vetted over the course of a decade informed by the NAS guidelines, and subsequently affirmed by new guidelines released in 2015 by the Hinxton Group, the International Society for Stem Cell Research (ISSCR), and 2017 recommendations by the NAS and National Academy of Medicine joint panel on human genome editing.
As part of the review process, OHSU convened additional ad hoc committees to evaluate the scientific merit and ethical justification of the proposed study: the OHSU Innovative Research Advisory Panel (IRAP) and a Scientific Review Committee (SRC). Members of both committees were independent and their names were kept confidential from the research team; OHSU Research Integrity supervised all committee meetings, documentation, and formal recommendations.
While international discussions were in their infancy, the OHSU Innovative Research Advisory Panel (IRAP) Committee was tasked with deliberating ethical considerations related to using gene correction technology in human embryos for basic research at OHSU. The committee was composed of eleven members from internal and external sources: a lay member, a clinical ObGyn physician, three bioethicists, an OHSU Institutional Ethics committee member, three former OSCRO members, a clinical geneticist, and a clinician. Upon completion of the review, the IRAP recommended allowing this research “with significant oversight and continued dialogue, the use of gene correction technologies in human embryos for the purpose of answering basic science questions needed to evaluate germline gene correction prior to the use in human models,” at OHSU.
The established track record of the study team to uphold strict confidentiality and regulatory requirements paved the way for full OHSU IRB study approval in 2016, contingent upon strict continuing oversight which includes: a phased scientific approach requiring evaluation of results on the safety and efficacy of germline gene correction in iPSCs before approving studies on human pre-implantation embryos; external bi-annual monitoring of all regulatory documents regarding human subjects; bi-annual Data Safety Monitoring Committee (DSMC) review; and annual continuing review by the OHSU IRB. The DSMC is required to remain active for the length of the approved IRB protocol and consists of four members: a lay member, an ethicist, a geneticist, and a reproductive endocrinologist. This committee conducts full review of all donations, the subsequent uses of these samples, and participant adverse events. The DSMC provides formal recommendations to the study team and IRB at the completion of each meeting.
The robust regulatory framework set forth by OHSU clearly specified that informed consent could be obtained only if prospective donors were made aware of the sensitive nature of the study. The consent form clearly presented the scientific rationale for the study; stating (in both the Clinical Research Consent Summary and the Purpose section of the consent form) that gene editing tools will be used on eggs, sperm, and/or embryos to evaluate the safety and efficacy of gene correction for heritable diseases. Additionally, consent form language clearly stated that genetic testing would be conducted in addition to creation of preimplantation embryos and embryonic stem cell lines for in vitro analyses and stored for future use. The incidental discovery of genetic information that might be important to the donors’ healthcare is a possible outcome when engaging in this type of research. Informed consent documents provided the donor with the option to receive this information or not. Written informed consent was obtained before all study-related procedures on current, IRB-approved, study-specific consent forms.
Healthy gamete donors were recruited locally, via print and web-based advertising. Homozygous and heterozygous adult patients with known heritable MYBPC3 mutations were sought; however, only three adult heterozygous patients were identified by OHSU Knight Cardiovascular Institute physicians and referred to the research team, one of whom agreed to participate in the study.
Controlled ovarian stimulation
Research oocyte donors were evaluated before study inclusion as previously reported; standard IVF protocols and procedures for ovarian stimulation were described previously44. Oocyte donation cycles were managed by OHSU Fertility physicians. Immediately following oocyte retrieval, recovered gametes were transferred to the research laboratory. All study-related procedures took place at the OHSU Center for Embryonic Cell and Gene Therapy. Following oocyte retrieval, cumulus-oocyte complexes (COCs) were treated with hyaluronidase to disaggregate cumulus and granulosa cells. Mature metaphase II (MII) oocytes were placed in Global Medium (LifeGlobal, IVFonline) supplemented with 10% serum substitute supplement (Global 10% medium) at 37 °C in 6% CO2 and covered with tissue culture oil (Sage IVF, Cooper Surgical).
All research donors were compensated for their time, effort, and discomfort associated with the donation process at rates similar to those used for gamete donation for fertility purposes.
Intracytoplasmic sperm injection (ICSI)
MII oocytes were placed into a 50-μl micromanipulation droplet of HTF (modified human tubal fluid) with HEPES 10% medium. The droplet was covered with tissue culture oil. The dish was then mounted on the stage of an inverted microscope (Olympus IX71) equipped with a stage warmer (http://www.tokaihit.com) and Narishige micromanipulators. Oocytes were fertilized by ICSI using frozen and thawed sperm. Fertilization was determined approximately 18 h after ICSI by noting the presence of two pronuclei and second polar body extrusion.
CRISPR–Cas9 injection into zygote or oocytes
For S-phase injections, zygotes were collected 18 h after ICSI and placed in a micromanipulation drop. The CRISPR–Cas9 mixture, containing Cas9 protein (200 ng/μl), sgRNA (100 ng/μl) and ssODN (200 ng/μl), was then injected into the cytoplasm. Injected zygotes were cultured in Global 10% medium at 37 °C in 6% CO2, 5% O2 and 89% N2 for up to 3 days to the 4–8-cell stage. For M-phase injections, CRISPR–Cas9 was co-injected with sperm during ICSI. A single sperm was first washed in a 4-μl drop of mixture containing Cas9 protein, sgRNA, and ssODN as described above.
Blastomere isolation, whole-genome amplification and Sanger sequencing
Zonae pellucidae from 4–8-cell stage embryos were removed by brief exposure to acidic Tyrode solution (NaCl 8 mg/ml, KCl 0.2 mg/ml, CaCl2.2H2O 2.4 mg/ml, MgCl2.6H2O 0.1 mg/ml, glucose 1 mg/ml, PVP 0.04 mg/ml). Zona-free embryos were briefly (30 s) exposed to a trypsin solution (0.15% in EDTA containing Ca- and Mg-free PBS) before manual disaggregation into single blastomeres with a small bore pipette. A total of 830 blastomeres were isolated from 131 embryos, including 19 from control, 54 from zygote-injected and 58 from M-phase-injected groups. Individual blastomeres were transferred into 0.2-ml PCR tubes containing 4 μl PBS and placed into a freezer at −80° C until further use. Whole-genome amplification was performed using a REPLI-g Single Cell Kit (Qiagen). Successful amplification was evaluated by PCR for MYBPC3 and Sanger sequencing. Briefly, amplified DNA was diluted 1/100 and the on-target region for MYBPC3 was amplified using a PCR Platinum SuperMix High Fidelity Kit (Life Technologies) with primer set: F 5′-CCCCCACCCAGGTACATCTT-3′ and R 5′-CTAGTGCACAGTGCATAGTG-3′. PCR products of 534 base pairs (bp) were purified, Sanger sequenced and analysed by Sequencher v5.0 (GeneCodes). Of 830 blastomeres, 730 (88%) resulted in successful libraries and produced PCR products for MYBPC3 while the remaining 100 blastomeres (12%) failed to generate PCR products and were excluded from the study.
iPSC derivation and transfection with CRISPR–Cas9
Patient iPSCs were derived from skin fibroblasts with a CytoTune-iPS Reprogramming Kit (Life Technologies), according to the manufacturer’s protocol. Cell lines were cultured in mTeSR1 medium (STEMCELL technology) at 37 °C in a humidified atmosphere containing 5% CO2. To test CRISPR–Cas9, 2 × 105 iPSCs were dissociated into single cells (using Accutase from STEMCELL technology, or TrypLe from Invitrogen). For the CRISPR–Cas9-1 construct (in the J.-S.K. laboratory), Cas9 expression plasmid (p3 s-Cas9HC, 2.4 μg), sgRNA expression plasmid (pU6-sgRNA, 1.6 μg) and ssODN-1 (100 pmol, IDT) were transfected into iPSCs using an Amaxa P3 Primary Cell 4D-Nucleofector Kit (Program CB-150) according to the manufacturer’s protocol. Three days after transfection, ~5,000 cells were plated onto a Matrigel-coated culture dish and cultured for clonal propagation and individual clone selection. For the CRISPR–Cas9-2 construct (in the J.C.I.B. laboratory), Cas9 expression plasmid (pCAG-1BPNLS-Cas9-1BPNLS, 15ug), sgRNA expression plasmids (pCAGmCherry-MYBPC3gRNA, 15ug) and 30 μg ssODN-2 (IDT) were co-transfected by electroporation using the BioRad Gene Pulser II (a single 320-V, 200-μF pulse at room temperature) with a 0.4-cm gap cuvette. Cells were plated at high density on 6-well plates coated with Matrigel. Two to three days after electroporation, iPSCs were harvested and subjected to clonal selection. All cell lines were negative for mycoplasma contamination. For direct comparisons of CRISPR–Cas9-1 and CRISPR–Cas9-2, Cas9 RNP complexes, composed of the recombinant Cas9 protein (15 μg) and sgRNA (20 μg), were co-transfected with ssODN-1 (50–200 pmol, IDT) into iPSCs (2 × 105 cells) via electroporation as described above. Three days after transfection, indel and HDR efficiencies were analysed by targeted deep sequencing.
Recombinant Cas9 protein and in vitro transcription of sgRNA
Recombinant Cas9 protein was purchased from ToolGen, Inc. The sgRNA was synthesized by in vitro transcription using T7 polymerase (New England Biolabs) as described previously45. In brief, sgRNA templates were generated by annealing and extension of two oligonucleotides (Extended Data Table 1). Then, in vitro transcription was performed by incubating sgRNA templates with T7 RNA polymerase supplemented with NTPs (Jena Bioscience) and RNase inhibitor (New England Biolabs) overnight at 37 °C. In vitro transcribed RNA was then treated with DNase I (New England Biolabs) for 30 min at 37 °C, and purified using MinElute Cleanup kit (Qiagen).
Targeted deep sequencing, genomic DNA cleavage, WGS and Digenome sequencing
To analyse HDR and NHEJ frequencies, on-target and off-target regions were amplified using Phusion polymerase (New England Biolabs). PCR amplicons were subjected to paired-end sequencing using Illumina Miniseq. Cas-analyzer was used for analysing indel and HDR frequencies46,47. The primers used for amplification are listed in Extended Data Table 5. Genomic DNA was isolated from patient iPSCs using a DNeasy Tissue Kit (Qiagen). Digenome-seq was performed as described25,26. In brief, 20 μg genomic DNA was cleaved by incubating recombinant Cas9 protein (16.7 μg) and in vitro transcribed sgRNA (12.5 μg) in 1× NEB buffer 3.1(100 mM NaCl, 50 mM Tris-HCl, 10 mM MgCl2, 100 μg/ml BSA, pH 7.9) at 37 °C for 3 h. Cas9- and sgRNA-treated genomic DNA was treated with 50 μg/ml RNase A (Sigma Aldrich) at 37 °C for 30 min, and purified with a DNeasy Tissue Kit (Qiagen). WGS and Digenome sequencing were performed as described previously25,26. In brief, 1 μg genomic DNA was fragmented and ligated with adaptors using TruSeq DNA libraries. DNA libraries were subjected to WGS using an Illumina HiSeq X Ten Sequencer at Macrogen (30× to 40×). The sequence file was aligned to the human reference genome hg19 from UCSC with the following mapping program and parameters using Isaac aligner48: base quality cutoff, 15; keep duplicate reads, yes; variable read length support, yes; realign gaps, no; and adaptor clipping, yes (adaptor: AGATCGGAAGAGC*, *GCTCTTCCGATCT). In vitro DNA cleavage sites were identified computationally using a DNA cleavage scoring system described previously26. Indel frequencies of 23 genomic loci with DNA cleavage score above the 0.1 cutoff value were individually examined in individual blastomeres by targeted deep sequencing. Primers for detecting indel frequencies of Digenome-seq captured sites are listed in Extended Data Table 5.
Analysis of off-target effects in CRISPR–Cas9-injected human embryos by WGS
WGS was performed using an Illumina HiSeq X Ten sequencer with a sequencing depth of 30× to 40× (Macrogen, South Korea). Sequences from each blastomere were processed to obtain total variants (lane 1 in Extended Data Table 6) using the Isaac variant calling program48. Annotated variants, including dbSNPs and all novel SNPs (substitution changes), were filtered out, and novel indel sites were identified (lane 2 in Extended Data Table 6). Cas-OFFinder46 was used to extract potential off-target sequences that differed from the on-target sequence by up to 7-nucleotide mismatches or up to 5-nucleotide mismatches with a DNA bulge of up to 2 nucleotides. Indel sites found in each blastomere were compared to homologous sites identified by Cas-OFFinder and potential off-target sites were identified (lanes 3 and 6 in Extended Data Table 6). Then, we excluded potential off-target sites, which were found in intact control embryos (C2.3 and C10.2, lanes 4 and 7 in Extended Data Table 6). Finally, we determined whether each of these potential off-target sites was caused by CRISPR-Cas9 by inspecting sequences with Integrative Genomics Viewer27 (IGV, lanes 5 and 8 in Extended Data Table 6).
Whole-exome sequencing and data analyses
Whole-exome sequencing (WES) was performed using genomic DNAs isolated from peripheral blood of the sperm donor and two egg donors (egg donor1 and egg donor2) and ES cells derived from individual human embryos (ES-WT1, ES-Mut1 and ES-C1 were from egg donor1; ES-WT2 and ES-WT3 were from egg donor2). ES-WT1, ES-WT2 and ES-WT3 were from treated wild-type embryos. ES-C1 was from an untreated wild-type embryo. ES-Mut1 was from a treated heterozygous mutant embryo. Sequencing libraries were prepared according to the instructions for Illumina library preparation. Exome capture was done using an Agilent V5 chip. Sequencing was done on an Illumina Hiseq 4000 platform with paired-end 101 (PE101) strategy at a depth of 100×. All sequencing data were first processed by filtering adaptor sequences and removing low quality reads or reads with a high percentage of N bases using SOAPnuke (1.5.2) software (http://soap.genomics.org.cn/) developed by BGI, and clean reads were generated for each library. Clean data were paired-end aligned using the Burrows-Wheeler Aligner (BWA) program version 0.7.12 to the human genome assembly hg19. Duplicate reads in alignment BAM files were identified using MarkDuplicates in Picard (1.54). The alignment results were processed by RealignerTargetCreator, IndelRealigner and BaseRecalibrator modules in GATK (3.3.0). Variant detection was performed using HaplotypeCaller tool in GATK. SNV and indel information was extracted and filtered by VQSR in GATK and annotated by AnnoDB (v3).
The guide sequence (GGGTGGAGTTTGTGAAGTAT) was aligned to the human genome assembly hg19 to identify potential off-target sites by full sensitive aligner Batmis (V3.00), allowing a maximum of five mismatches globally and a maximum of two mismatches in the core region (12 bp adjacent to the PAM site). Inherited variants from parents and all novel SNPs (substitution changes) were filtered out, and novel indels located within the off-target site plus flanking 20-bp region were defined as off-target variants.
Student’s t-test was used for the comparisons in Fig. 4a. One-tailed Fisher’s test was used for the comparisons in Fig. 2f, Fig. 3c, and Extended Data Fig. 1e. One-way ANOVA with Bonferroni test was used for the comparisons in Extended Data Table 2. A P value less than 0.05 was considered significant. No statistical methods were used to predetermine sample size. The experiments were randomized and the investigators were blinded to allocation during experiments and outcome assessment whenever possible.
The sequencing data sets, including WGS, WES and Digenome-seq, generated during the current study are not being made publicly available owing to concerns that the data could reveal the research participant’s genetic identity, and revealing the identity would be against the participant’s wishes and consent. However, the data will be made available to researchers from the corresponding author on reasonable request, dependent upon OHSU IRB approval.
Antoniou, A. et al. Average risks of breast and ovarian cancer associated with BRCA1 or BRCA2 mutations detected in case series unselected for family history: a combined analysis of 22 studies. Am. J. Hum. Genet. 72, 1117–1130 (2003)
Carrier, L., Mearini, G., Stathopoulou, K. & Cuello, F. Cardiac myosin-binding protein C (MYBPC3) in cardiac pathophysiology. Gene 573, 188–197 (2015)
Dhandapany, P. S. et al. A common MYBPC3 (cardiac myosin binding protein C) variant associated with cardiomyopathies in South Asia. Nat. Genet. 41, 187–191 (2009)
Struewing, J. P. et al. The risk of cancer associated with specific mutations of BRCA1 and BRCA2 among Ashkenazi Jews. N. Engl. J. Med. 336, 1401–1408 (1997)
Maron, B. J. et al. Prevalence of hypertrophic cardiomyopathy in a general population of young adults. Echocardiographic analysis of 4111 subjects in the CARDIA Study. Circulation 92, 785–789 (1995)
Schlossarek, S., Mearini, G. & Carrier, L. Cardiac myosin-binding protein C in hypertrophic cardiomyopathy: mechanisms and therapeutic opportunities. J. Mol. Cell. Cardiol. 50, 613–620 (2011)
Hsu, P. D., Lander, E. S. & Zhang, F. Development and applications of CRISPR–Cas9 for genome engineering. Cell 157, 1262–1278 (2014)
Mali, P. et al. RNA-guided human genome engineering via Cas9. Science 339, 823–826 (2013)
Kim, S., Kim, D., Cho, S. W., Kim, J. & Kim, J. S. Highly efficient RNA-guided genome editing in human cells via delivery of purified Cas9 ribonucleoproteins. Genome Res. 24, 1012–1019 (2014)
Cong, L. et al. Multiplex genome engineering using CRISPR/Cas systems. Science 339, 819–823 (2013)
Lin, S., Staahl, B. T., Alla, R. K. & Doudna, J. A. Enhanced homology-directed human genome engineering by controlled timing of CRISPR/Cas9 delivery. eLife 3, e04766 (2014)
Wu, Y. et al. Correction of a genetic disease in mouse via use of CRISPR–Cas9. Cell Stem Cell 13, 659–662 (2013)
Guo, X. & Li, X. J. Targeted genome editing in primate embryos. Cell Res. 25, 767–768 (2015)
Wang, H. et al. One-step generation of mice carrying mutations in multiple genes by CRISPR/Cas-mediated genome engineering. Cell 153, 910–918 (2013)
Kang, E. et al. Age-related accumulation of somatic mitochondrial DNA mutations in adult-derived human iPSCs. Cell Stem Cell 18, 625–636 (2016)
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, 230–232 (2013)
Kim, H. & Kim, J. S. A guide to genome engineering with programmable nucleases. Nat. Rev. Genet. 15, 321–334 (2014)
Jinek, M. et al. A programmable dual-RNA-guided DNA endonuclease in adaptive bacterial immunity. Science 337, 816–821 (2012)
Aida, T. et al. Cloning-free CRISPR/Cas system facilitates functional cassette knock-in in mice. Genome Biol. 16, 87 (2015)
Capmany, G., Taylor, A., Braude, P. R. & Bolton, V. N. The timing of pronuclear formation, DNA synthesis and cleavage in the human 1-cell embryo. Mol. Hum. Reprod. 2, 299–306 (1996)
Suzuki, T., Asami, M. & Perry, A. C. Asymmetric parental genome engineering by Cas9 during mouse meiotic exit. Sci. Rep. 4, 7621 (2014)
Fu, Y. et al. High-frequency off-target mutagenesis induced by CRISPR–Cas nucleases in human cells. Nat. Biotechnol. 31, 822–826 (2013)
Hsu, P. D. et al. DNA targeting specificity of RNA-guided Cas9 nucleases. Nat. Biotechnol. 31, 827–832 (2013)
Cho, S. W. et al. Analysis of off-target effects of CRISPR/Cas-derived RNA-guided endonucleases and nickases. Genome Res. 24, 132–141 (2014)
Kim, D. et al. Digenome-seq: genome-wide profiling of CRISPR–Cas9 off-target effects in human cells. Nat. Methods 12, 237–243 (2015)
Kim, D., Kim, S., Kim, S., Park, J. & Kim, J. S. Genome-wide target specificities of CRISPR–Cas9 nucleases revealed by multiplex Digenome-seq. Genome Res. 26, 406–415 (2016)
Robinson, J. T. et al. Integrative genomics viewer. Nat. Biotechnol. 29, 24–26 (2011)
Schneider, T. D. & Stephens, R. M. Sequence logos: a new way to display consensus sequences. Nucleic Acids Res. 18, 6097–6100 (1990)
Tennakoon, C., Purbojati, R. W. & Sung, W. K. BatMis: a fast algorithm for k-mismatch mapping. Bioinformatics 28, 2122–2128 (2012)
Richardson, C. D., Ray, G. J., DeWitt, M. A., Curie, G. L. & Corn, J. E. Enhancing homology-directed genome editing by catalytically active and inactive CRISPR–Cas9 using asymmetric donor DNA. Nat. Biotechnol. 34, 339–344 (2016)
Doudna, J. A. & Charpentier, E. The new frontier of genome engineering with CRISPR–Cas9. Science 346, 1258096 (2014)
Luo, Y. et al. Hypersensitivity of primordial germ cells to compromised replication-associated DNA repair involves ATM-p53-p21 signaling. PLoS Genet. 10, e1004471 (2014)
Lange, J. et al. The landscape of mouse meiotic double-strand break formation, processing, and repair. Cell 167, 695–708 (2016)
Titus, S. et al. Impairment of BRCA1-related DNA double-strand break repair leads to ovarian aging in mice and humans. Sci. Transl. Med. 5, 172ra21 (2013)
Tang, L. et al. CRISPR/Cas9-mediated gene editing in human zygotes using Cas9 protein. Mol. Genet. Genomics 292, 525–533 (2017)
Chu, V. T. et al. Increasing the efficiency of homology-directed repair for CRISPR–Cas9-induced precise gene editing in mammalian cells. Nat. Biotechnol. 33, 543–548 (2015)
Maruyama, T. et al. Increasing the efficiency of precise genome editing with CRISPR–Cas9 by inhibition of nonhomologous end joining. Nat. Biotechnol. 33, 538–542 (2015)
Niu, Y. et al. Generation of gene-modified cynomolgus monkey via Cas9/RNA-mediated gene targeting in one-cell embryos. Cell 156, 836–843 (2014)
Kang, Y. et al. CRISPR/Cas9-mediated Dax1 knockout in the monkey recapitulates human AHC-HH. Hum. Mol. Genet. 24, 7255–7264 (2015)
Liang, P. et al. CRISPR/Cas9-mediated gene editing in human tripronuclear zygotes. Protein Cell 6, 363–372 (2015)
Tu, Z. et al. Promoting Cas9 degradation reduces mosaic mutations in non-human primate embryos. Sci. Rep. 7, 42081 (2017)
Hashimoto, M., Yamashita, Y. & Takemoto, T. Electroporation of Cas9 protein/sgRNA into early pronuclear zygotes generates non-mosaic mutants in the mouse. Dev. Biol. 418, 1–9 (2016)
Orthwein, A. et al. Mitosis inhibits DNA double-strand break repair to guard against telomere fusions. Science 344, 189–193 (2014)
Tachibana, M. et al. Towards germline gene therapy of inherited mitochondrial diseases. Nature 493, 627–631 (2013)
Kim, J. M., Kim, D., Kim, S. & Kim, J. S. Genotyping with CRISPR–Cas-derived RNA-guided endonucleases. Nat. Commun. 5, 3157 (2014)
Bae, S., Park, J. & Kim, J. S. Cas-OFFinder: a fast and versatile algorithm that searches for potential off-target sites of Cas9 RNA-guided endonucleases. Bioinformatics 30, 1473–1475 (2014)
Park, J., Lim, K., Kim, J. S. & Bae, S. Cas-analyzer: an online tool for assessing genome editing results using NGS data. Bioinformatics 33, 286–288 (2017)
Raczy, C. et al. Isaac: ultra-fast whole-genome secondary analysis on Illumina sequencing platforms. Bioinformatics 29, 2041–2043 (2013)
We acknowledge the OHSU Institutional Review Board (IRB), Innovative Research Advisory Panel (IRAP), Scientific Review Committee (SRC) and Data Safety Monitoring Committee (DSMC) for oversight and guidance on this study. We thank all study participants for tissue donations; the Women’s Health Research Unit staff, IVF laboratory staff, University Fertility Consultants and the Reproductive Endocrinology and Infertility Division in the Department of Obstetrics and Gynecology, OHSU for support and procurement of human gametes; S. Olson and Research Cytogenetics Laboratory at OHSU for cytogenetic analysis of ES cells; S. Cooper from the Wallace Division of Smiths Medical for donating ICSI and gamete manipulation micropipettes; Y. Wang, T. Wu and Y. Shen from BGI-Shenzhen for help with sample preparation, sequencing and data analyses; M. Ku from the H. A. and Mary K. Chapman Charitable Foundations Genomic Sequencing Core of the Salk Institute for next generation sequencing; and E. Aizawa and R. Hernandez-Benitez from the laboratory of JCIB for assistance. Studies conducted at OHSU were supported by OHSU institutional funds. Work in the laboratory of J.-S.K. was supported by the Institute for Basic Science (IBS-R021-D1). Work in the laboratory of J.C.I.B. was supported by the G. Harold and Leila Y. Mathers Charitable Foundation, the Moxie Foundation and the Leona M. and Harry B. Helmsley Charitable Trust. Work at BGI was supported by the Shenzhen Municipal Government of China (DRC-SZ  884).
J.-S.K. is a co-founder and shareholder of ToolGen, Inc, and S.M. is a co-founder and shareholder of Mitogenome Therapeutics, Inc and Health and Science Center “M1” Inc and declare competing financial interests. The other authors declare no competing financial interests.
Reviewer Information Nature thanks F. Lanner and the other anonymous reviewer(s) for their contribution to the peer review of this work.
Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Extended data figures and tables
a, b, Schematic of CRISPR–Cas9-1 and CRISPR–Cas9-2 constructs. Both systems consist of a single-chain chimaeric sgRNA designed to target the MYBPC3∆GAGT deletion and Cas9 protein. Exogenous single-stranded oligodeoxynucleotide (ssODN) templates encoding homology arms to the targeted region were designed for each system to facilitate HDR. Synonymous single-nucleotide substitutions were introduced into each ssODN template as indicated by red letters. Underline in 1a represents wild-type sequence of 4-bp deletion in the ssODN-1. In addition, the ssODN-2 nucleotide substitutions provide a restriction enzyme (BstBI) recognition site as indicated with black open box. c, Patient iPSCs were transfected with CRISPR–Cas9 plasmids by electroporation and individual single iPSC clones were analysed. d, Representative chromatographs showing untargeted iPSC clone with heterozygous mutant (left), targeted iPSC clone with gene corrected via HDR using ssODN-2 as repair template (middle) and targeted iPSC clone with gene corrected via HDR using wild-type sequence as template. e, Targeting and HDR efficiency comparison between CRISPR–Cas9-1 and CRISPR–Cas9-2. The asterisk indicates a significantly higher targeting efficiency of CRISPR–Cas9-1 compared to CRISPR–Cas9-2 (P Supplementary Table 1. f, HDR and NHEJ efficiency in wild-type ES cells (H9) and patient iPSCs transfected with preassembled Cas9 ribonucleoproteins (RNPs).
Extended Data Figure 2 Differences in DNA repair responses in human embryos and iPSCs and independent validation of on-target genotyping.
a, Different HDR-based DNA repair was observed in patient embryos compared to iPSCs. In human embryos, induced DSBs in the mutant paternal gene were repaired exclusively using the wild-type maternal allele as a template. In striking contrast, HDR in iPSCs was directed by the external ssODN-1 template. b, Verification of blastomere on-target genotypes by deep sequencing. DNA samples from 86 individual blastomeres representing different treatments and controls were sequenced using the MiSeq platform. Total reads, wild-type reads and mutant reads are presented as a mean. Distribution of mutant reads in heterozygous (34.3–64%) or homozygous (0.01–0.1%) blastomeres in experimental embryos was similar to controls.
a, Origin and genotypes of ES cells derived from CRISPR–Cas9-injected embryos. b, Karyotype of ES-WT3 showing a pericentric inversion on chromosome 10. c, Magnified image showing normal (left chromosome) and inverted chromosome 10 (right chromosome, inverted region indicated by arrow). d, Chromosome 10 ideogram showing the location of the inversion.
a, Representative Integrative Genomics Viewer (IGV) images produced using CRISPR–Cas9 at the on-target site. b, Representative IGV images showing CRISPR–Cas9-induced DNA cleavage at the potential off-target sites. Red arrows indicate DNA cleavage sites at each off-target site.
a, Genome-wide Circos plots showing DNA cleavage scores. Cas9-only-treated DNA is shown in grey, and CRISPR–Cas9-treated DNA is in blue. b, Sequence logo obtained via WebLogo using Digenome-captured sites (DNA cleavage score >2.5). On-target sequence is indicated below the sequence logo. PAM sequence is shown in blue. c, On-target indels for 28 blastomeres detected by Digenome-seq. Only blastomeres carrying NHEJ were captured by Digenome-seq. Blastomeres indicated in red fonts were also analysed by WGS (Extended Data Table 6). d, Indel frequencies for potential off-target sites captured by Digenome-seq in CRISPR–Cas9-treated (n = 5) and untreated control embryos (n = 2). Mismatch nucleotides are shown in red. PAM sequence is shown in blue. OnT: on target site; OT: off-target. Percentages of mutagenic indels on the x-axis are presented on a logarithmic scale. For source data, see Supplementary Table 5.
a, SNP and indel fractions of variants identified in each sample presented in dbSNP and 1000genomes databases. b, Numbers of coding region variants (SNPs) in each sample divided into indicated categories. c, Numbers of coding region variants (indels) divided into indicated categories. Hg19 was used as a reference genome for variant calling. Data presented as mean with interquartile range (black horizontal bars). For source data, see Supplementary Tables 6–8.
Genotypes of individual iPSC clones
Genotypes of individual blastomeres in control embryos
Genotypes of individual blastomeres in zygote injected embryos
Genotypes of individual blastomeres in M-phase injected embryos
Mutagenic indel frequencies at potential off-target sites captured by Digenome-sequencing
Whole exome sequencing (WES) analyses of variants identified using Hg19 as reference genome
List of predicted off-target sites
WES analyses of non-inherited variants
CRISPR/Cas9 injection into human oocytes and zygotes. S-phase CRISPR/Cas9 injection was performed 18 hours after ICSI. M-phase injections were co-injected with sperm during ICSI.
About this article
Cite this article
Ma, H., Marti-Gutierrez, N., Park, S. et al. Correction of a pathogenic gene mutation in human embryos. Nature 548, 413–419 (2017). https://doi.org/10.1038/nature23305
Science China Life Sciences (2020)
Biomedicine & Pharmacotherapy (2020)
Current Problems in Cardiology (2020)
Suppression of unwanted CRISPR-Cas9 editing by co-administration of catalytically inactivating truncated guide RNAs
Nature Communications (2020)
Nature Genetics (2020)