Advertisement
Methods and Resources| Volume 294, ISSUE 4, P1142-1151, January 25, 2019

Download started.

Ok

Optimizing CRISPR/Cas9 technology for precise correction of the Fgfr3-G374R mutation in achondroplasia in mice

  • Kai Miao
    Footnotes
    Affiliations
    Cancer Center, Faculty of Health Sciences, Macau SAR

    Centre for Precision Medicine Research and Training, Faculty of Health Sciences, University of Macau, Macau SAR
    Search for articles by this author
  • Xin Zhang
    Footnotes
    Affiliations
    Cancer Center, Faculty of Health Sciences, Macau SAR

    Transgenic and Knockout Core, Faculty of Health Sciences, University of Macau, Macau SAR, China
    Search for articles by this author
  • Sek Man Su
    Affiliations
    Cancer Center, Faculty of Health Sciences, Macau SAR

    Centre for Precision Medicine Research and Training, Faculty of Health Sciences, University of Macau, Macau SAR
    Search for articles by this author
  • Jianming Zeng
    Affiliations
    Cancer Center, Faculty of Health Sciences, Macau SAR

    Centre for Precision Medicine Research and Training, Faculty of Health Sciences, University of Macau, Macau SAR
    Search for articles by this author
  • Zebin Huang
    Affiliations
    Cancer Center, Faculty of Health Sciences, Macau SAR

    Centre for Precision Medicine Research and Training, Faculty of Health Sciences, University of Macau, Macau SAR
    Search for articles by this author
  • Un In Chan
    Affiliations
    Cancer Center, Faculty of Health Sciences, Macau SAR

    Transgenic and Knockout Core, Faculty of Health Sciences, University of Macau, Macau SAR, China
    Search for articles by this author
  • Xiaoling Xu
    Correspondence
    To whom correspondence may be addressed. Tel.: 853-8822-4925; Fax: 853-8822-2911
    Affiliations
    Cancer Center, Faculty of Health Sciences, Macau SAR

    Centre for Precision Medicine Research and Training, Faculty of Health Sciences, University of Macau, Macau SAR

    Transgenic and Knockout Core, Faculty of Health Sciences, University of Macau, Macau SAR, China
    Search for articles by this author
  • Chu-Xia Deng
    Correspondence
    To whom correspondence may be addressed. Tel.: 853-8822-4997; Fax: 853-8822-2314
    Affiliations
    Cancer Center, Faculty of Health Sciences, Macau SAR

    Centre for Precision Medicine Research and Training, Faculty of Health Sciences, University of Macau, Macau SAR
    Search for articles by this author
  • Author Footnotes
    1 Both authors contributed equally to the results of this work.
Open AccessPublished:November 28, 2018DOI:https://doi.org/10.1074/jbc.RA118.006496
      CRISPR/Cas9 is a powerful technology widely used for genome editing, with the potential to be used for correcting a wide variety of deleterious disease-causing mutations. However, the technique tends to generate more indels (insertions and deletions) than precise modifications at the target sites, which might not resolve the mutation and could instead exacerbate the initial genetic disruption. We sought to develop an improved protocol for CRISPR/Cas9 that would correct mutations without unintended consequences. As a case study, we focused on achondroplasia, a common genetic form of dwarfism defined by missense mutation in the Fgfr3 gene that results in glycine to arginine substitution at position 374 in mice in fibroblast growth factor receptor 3 (Fgfr3-G374R), which corresponds to G380R in humans. First, we designed a GFP reporter system that can evaluate the cutting efficiency and specificity of single guide RNAs (sgRNAs). Using the sgRNA selected based on our GFP reporter system, we conducted targeted therapy of achondroplasia in mice. We found that we achieved higher frequency of precise correction of the Fgfr3-G374R mutation using Cas9 protein rather than Cas9 mRNA. We further demonstrated that targeting oligos of 100 and 200 nucleotides precisely corrected the mutation at equal efficiency. We showed that our strategy completely suppressed phenotypes of achondroplasia and whole genome sequencing detected no off-target effects. These data indicate that improved protocols can enable the precise CRISPR/Cas9-mediated correction of individual mutations with high fidelity.

      Introduction

      CRISPR/Cas9 is a RNA-guided genome editing tool developed from a microbial adaptive immune defense system (
      • Cong L.
      • Ran F.A.
      • Cox D.
      • Lin S.
      • Barretto R.
      • Habib N.
      • Hsu P.D.
      • Wu X.
      • Jiang W.
      • Marraffini L.A.
      • Zhang F.
      Multiplex genome engineering using CRISPR/Cas systems.
      ,
      • Mali P.
      • Yang L.
      • Esvelt K.M.
      • Aach J.
      • Guell M.
      • DiCarlo J.E.
      • Norville J.E.
      • Church G.M.
      RNA-guided human genome engineering via Cas9.
      • Knott G.J.
      • Doudna J.A.
      CRISPR-Cas guides the future of genetic engineering.
      ). Because of its high efficiency and easy-to-use nature, the CRISPR/Cas9 system is revolutionizing many research fields, including gene functional study, medical research, and genomic editing. Since its introduction into mammalian cells in 2013, many important progresses have been made in various scientific fields (
      • Chen S.
      • Sun H.
      • Miao K.
      • Deng C.X.
      CRISPR-Cas9: from genome editing to cancer res.
      ,
      • Luther D.C.
      • Lee Y.W.
      • Nagaraj H.
      • Scaletti F.
      • Rotello V.M.
      Delivery approaches for CRISPR/Cas9 therapeutics in vivo: advances and challenges.
      • Chen S.
      • Sanjana N.E.
      • Zheng K.
      • Shalem O.
      • Lee K.
      • Shi X.
      • Scott D.A.
      • Song J.
      • Pan J.Q.
      • Weissleder R.
      • Lee H.
      • Zhang F.
      • Sharp P.A.
      Genome-wide CRISPR screen in a mouse model of tumor growth and metastasis.
      ). In particular for gene therapy, several studies have demonstrated the power of CRISPR/Cas9 technology in mouse models for the correction of human hereditary genetic diseases, such as, cataract disorder (
      • Wu Y.
      • Liang D.
      • Wang Y.
      • Bai M.
      • Tang W.
      • Bao S.
      • Yan Z.
      • Li D.
      • Li J.
      Correction of a genetic disease in mouse via use of CRISPR-Cas9.
      ), tyrosinemia (
      • Yin H.
      • Xue W.
      • Chen S.
      • Bogorad R.L.
      • Benedetti E.
      • Grompe M.
      • Koteliansky V.
      • Sharp P.A.
      • Jacks T.
      • Anderson D.G.
      Genome editing with Cas9 in adult mice corrects a disease mutation and phenotype.
      ), thalassemia (
      • Ou Z.
      • Niu X.
      • He W.
      • Chen Y.
      • Song B.
      • Xian Y.
      • Fan D.
      • Tang D.
      • Sun X.
      The combination of CRISPR/Cas9 and iPSC technologies in the gene therapy of human β-thalassemia in mice.
      ), Duchenne muscular dystrophy (
      • Tabebordbar M.
      • Zhu K.
      • Cheng J.K.W.
      • Chew W.L.
      • Widrick J.J.
      • Yan W.X.
      • Maesner C.
      • Wu E.Y.
      • Xiao R.
      • Ran F.A.
      • Cong L.
      • Zhang F.
      • Vandenberghe L.H.
      • Church G.M.
      • Wagers A.J.
      In vivo gene editing in dystrophic mouse muscle and muscle stem cells.
      ,
      • Nelson C.E.
      • Hakim C.H.
      • Ousterout D.G.
      • Thakore P.I.
      • Moreb E.A.
      • Castellanos Rivera R.M.
      • Madhavan S.
      • Pan X.
      • Ran F.A.
      • Yan W.X.
      • Asokan A.
      • Zhang F.
      • Duan D.
      • Gersbach C.A.
      In vivo genome editing improves muscle function in a mouse model of Duchenne muscular dystrophy.
      • Long C.
      • Amoasii L.
      • Mireault A.A.
      • McAnally J.R.
      • Li H.
      • Sanchez-Ortiz E.
      • Bhattacharyya S.
      • Shelton J.M.
      • Bassel-Duby R.
      • Olson E.N.
      Postnatal genome editing partially restores dystrophin expression in a mouse model of muscular dystrophy.
      ), liver diseases (
      • Yang Y.
      • Wang L.
      • Bell P.
      • McMenamin D.
      • He Z.
      • White J.
      • Yu H.
      • Xu C.
      • Morizono H.
      • Musunuru K.
      • Batshaw M.L.
      • Wilson J.M.
      A dual AAV system enables the Cas9-mediated correction of a metabolic liver disease in newborn mice.
      ,
      • Yin H.
      • Song C.Q.
      • Dorkin J.R.
      • Zhu L.J.
      • Li Y.
      • Wu Q.
      • Park A.
      • Yang J.
      • Suresh S.
      • Bizhanova A.
      • Gupta A.
      • Bolukbasi M.F.
      • Walsh S.
      • Bogorad R.L.
      • Gao G.
      • Weng Z.
      • Dong Y.
      • Koteliansky V.
      • Wolfe S.A.
      • Langer R.
      • Xue W.
      • Anderson D.G.
      Therapeutic genome editing by combined viral and non-viral delivery of CRISPR system components in vivo.
      ) and so on. These approaches are potentially translatable for human clinical therapy. Meanwhile studies also revealed some technical difficulties that CRISPR/Cas9 system encounters, such as it preferably generates more small insertions and deletions (Indels) than precise modifications, and off-target effect, i.e. binding and cutting at sites with shared homology (
      • Knott G.J.
      • Doudna J.A.
      CRISPR-Cas guides the future of genetic engineering.
      ,
      • Martin F.
      • Sánchez-Hernandez S.
      • Gutiérrez-Guerrero A.
      • Pinedo-Gomez J.
      • Benabdellah K.
      Biased and unbiased methods for the detection of off-target cleavage by CRISPR/Cas9: an overview.
      ,
      • Rastogi A.
      • Murik O.
      • Bowler C.
      • Tirichine L.
      PhytoCRISP-Ex: a web-based and stand-alone application to find specific target sequences for CRISPR/CAS editing.
      • Haeussler M.
      • Schönig K.
      • Eckert H.
      • Eschstruth A.
      • Mianné J.
      • Renaud J.B.
      • Schneider-Maunoury S.
      • Shkumatava A.
      • Teboul L.
      • Kent J.
      • Joly J.S.
      • Concordet J.P.
      Evaluation of off-target and on-target scoring algorithms and integration into the guide RNA selection tool CRISPOR.
      ).
      Achondroplasia is the most common genetic form of dwarfism inherited as an autosomal dominant disorder (
      • Ornitz D.M.
      • Legeai-Mallet L.
      Achondroplasia: development, pathogenesis, and therapy.
      ). People with achondroplasia have impaired ability in the longitudinal growth of long bone from endochondral ossification at the epiphyseal growth plate, leading to the short status. Homozygous achondroplasia predisposes its carriers to neonatal lethal condition, whereas heterozygous achondroplastic adults can be as short as 62.8 cm (24.7 in). More than 96% of patients with achondroplasia have a Gly to Arg transition at position 1138 (G1138A) of the fibroblast growth factor receptor 3 (FGFR3) gene, resulting in the Gly to Arg substitution at position 380 of the FGFR3 protein (
      • Bellus G.A.
      • Hefferon T.W.
      • Ortiz de Luna R.I.
      • Hecht J.T.
      • Horton W.A.
      • Machado M.
      • Kaitila I.
      • McIntosh I.
      • Francomano C.A.
      Achondroplasia is defined by recurrent G380R mutations of FGFR3.
      ,
      • Stoilov I.
      • Kilpatrick M.W.
      • Tsipouras P.
      A common FGFR3 gene mutation is present in achondroplasia but not in hypochondroplasia.
      ), which corresponds to the Gly to Arg at 374 (G1120A for DNA) in mouse Fgfr3 (
      • Wang J.M.
      • Du. X.L.
      • Li C.L.
      • Yin L.J.
      • Chen B.
      • Sun J.
      • Su N.
      • Zhao L.
      • Song R.H.
      • Song W.W.
      • Chen L.
      • Deng C.X.
      Gly374Arg mutation in Fgfr3 causes achondroplasia in mice.
      ). Although human growth hormone has been used to aid growth in other dwarfism, it does not help people with achondroplasia. The only way is the controversial surgery of limb-lengthening, which will lead to patients suffering huge pain. Therefore, developing a radical treatment is desired.
      Using our previously generated mouse Fgfr3-G374R (corresponds to Fgfr3G1120A) achondroplasia model (
      • Wang J.M.
      • Du. X.L.
      • Li C.L.
      • Yin L.J.
      • Chen B.
      • Sun J.
      • Su N.
      • Zhao L.
      • Song R.H.
      • Song W.W.
      • Chen L.
      • Deng C.X.
      Gly374Arg mutation in Fgfr3 causes achondroplasia in mice.
      ), here we show that high frequency of targeted correction of the G374R can be achieved by co-injection into zygotes of Cas9 protein and a single guide RNA (sgRNA)
      The abbreviations used are: sgRNA
      single guide RNA
      HDR
      homology-directed repair
      nt
      nucleotide(s)
      PAM
      protospacer adjacent motif
      lnl
      Loxp-neo-Loxp
      oligo
      oligonucleotide
      ssDO
      single strand donor DNA oligo
      EGFP
      enhanced green fluorescent protein.
      optimized for specifically targeting the mutant allele. The correction occurs via homology-directed repair (HDR) based on an exogenously supplied single strand donor DNA without inducing nonspecific off targeting events. The resulting mice showed normal body size, were fertile, and able to transmit the corrected allele to their progeny. Thus, our study provides proof for using the CRISPR/Cas9 system to correct this dominate genetic disease.

      Results

      Evaluation of genome editing efficiency and specificity of sgRNAs against Fgfr3-G374R by sgRNA testers

      Three overlapping sgRNAs of 20 nt were designed to spanning the Gly to Arg mutation in Fgfr3 gene (Fig. S1A), the top 10 predicted off-target locus with NGG PAM (protospacer adjacent motif) sequence were listed for further analysis (Fig. S1, B–D). The highest homology for these sgRNAs is WT Fgfr3 that contains 1 mismatched base at the “Gly to Arg” mutation site. This is followed by 1 mismatches of the Cnnm2 gene to sgRNA1; 2 mismatches of noncoding region locus Chr2: 10643681 to sgRNA2, and 2 mismatches of Extl3 genes to sgRNA3, respectively. Increasing mismatches to these sgRNAs up to 4 bases in some other genes of varying homology are shown in Fig. S1, B–D.
      We then tested the cutting efficiency of the sgRNAs using a GFP reporter plasmid (pCAG-EGx-xFP-BbsI), which carries tandemly duplicated defective GFPs with a mutation in their 3′ and 5′, respectively. A dual BbsI cutting site was also introduced between two mutated GFPs for cloning of sgDNA+NGG to be tested (Fig. 1A). For each sgRNA, we cloned the mutant form and the corresponding WT form, which carries one base mismatch with the mutant form, into the tester construct (Fig. 1B) to generate 3 pairs of tester reporters. Then, each of the reporters alone or together with their corresponding CRISPR/Cas9 construct were transfected into 293 cells to test their cutting efficiency and specificity.
      Figure thumbnail gr1
      Figure 1Cutting efficiency and off-target assay for sgRNAs. A, map for pCAG-EGx-xFP-BbsI tester, a dual BbsI cutting site was flanked by a 3′-mutated GFP and a 5′-mutated GFP. sgRNA oligos plus NGG PAM sequence are inserted into the pCAG-EGxxFP-BbsI plasmid through the dual BbsI cutting site to generate the tester reporter for each single target sgRNA. B, sequence of sgRNA oligos used for targeting vector (PX330) and tester vector cloning. C, tester assay for the Fgfr3 mutant sgRNA targeting constructs. PX330-sgRNA2 and its corresponding mutant tester or WT tester are co-transfected into 293FT cells, 24 h later, the cells were examined for GFP intensity to estimate the cutting efficiency. D, quantification of GFP positive cells and its intensity related to testers construct generated by sgRNA2.
      Our data indicated that transfection of these testers alone did not generate obvious GFP signaling, which is consistent with the fact that tandem duplicated GFPs are defective in the pCAG-EGx-xFP-BbsI construct. On the other hand, co-transfection of all three pairs of testers with their corresponding sgRNA-CRISPR/Cas9 (PX330-sgRNA) constructs generated GFP positive cells. This is because the mutant GFPs can be repaired due to homologous recombination stimulated by CRISPR/Cas9-mediated double strand break at the target site, generating GFP+ cells (Fig. 1C). The GFP expression level should reflect the cutting efficiency of sgRNA-guided CRISPR/Cas9 to the targeted sequence. The data indicated that sgRNA1 had the highest cutting effect for the Fgfr3 mutant sequence, however, it also cut the corresponding WT sequence (one base mismatch to the MT sequence) with high efficiency. sgRNA3 cut much more efficiently than sgRNA1, however, sgRNA3 also cut the WT sequence and generated obvious GFP positive signaling (Fig. S2, A–C). Thus, sgRNA2 had moderate efficiency for cutting the Fgfr3 mutant allele among the three sgRNAs tested, whereas no cutting at its corresponding WT allele (Fig. 1, C and D), reflecting its high specificity to the mutant allele. Because avoiding the off-target effect is most crucial, we selected sgRNA2 for our subsequent in vivo experiments.

      Targeted correction of dwarf phenotype associated with Fgfr3-G374R mutation in mice

      We used the mouse model of dominant achondroplasia caused by a defined Gly to Arg mutation in exon 10 of the Fgfr3 gene (
      • Wang J.M.
      • Du. X.L.
      • Li C.L.
      • Yin L.J.
      • Chen B.
      • Sun J.
      • Su N.
      • Zhao L.
      • Song R.H.
      • Song W.W.
      • Chen L.
      • Deng C.X.
      Gly374Arg mutation in Fgfr3 causes achondroplasia in mice.
      ). As achondroplasia mice grow poorly throughout development and have reduced fertility, a Loxp-neo-Loxp (lnl) cassette was inserted into intron 10 to block the expression of the mutant allele (Fig. 2A). Heterozygous mice for the mutation (Fgfr3G1120A−lnl/+) were normal, whereas homozygous mice (Fgfr3G1120A-lnl/G1120A-lnl) survived to adulthood and were fertile, although exhibited skeleton abnormalities similar to Fgfr3−/− mice generated by gene targeting (
      • Deng C.
      • Wynshaw-Boris A.
      • Zhou F.
      • Kuo A.
      • Leder P.
      Fibroblast growth factor receptor 3 is a negative regulator of bone growth.
      ). Deletion of the lnl cassette from mutant mice by breeding with mice carrying the Ella-Cre transgene will allow the expression of Fgfr3G1120A allele, generating mice with a dwarf phenotype (Fig. 2A). Our data indicated that the cross between homozygous male Fgfr3G1120A-lnl/G1120A-lnl mice with homozygous female Ella-Cre mice could generate dwarf phenotypes in all offspring (Fgfr3G1120A/+; Ella-Cre, or Fgfr3G1120A), thus the fertilized zygotes from this breeding were collected for pronuclear injection.
      Figure thumbnail gr2
      Figure 2Targeted correction of Fgfr3-G374R mutation by CRISPR/Cas9 system. A, flow chat of the breeding procedure. Male mice homozygous for Fgfr3G1120A-ln/lG1120A-lnl, which phenotypically mimic Fgfr3−/− mice, were crossed with homozygous Ella-Cre transgenic female mice to cut off the ploxPneo gene to activate the Fgfr3 point-mutated gene. All offspring were heterozygous for the activated Fgfr3 mutation (Fgfr3G1120A/+) and exhibited achondroplasia phenotypes (an example was shown in the lower right panel). Oocytes were isolated from this cross and Cas9 protein or mRNA, sgRNA, and ssDO were injected into the pronuclear to correction the mutation. B, the sequence of mutant Fgfr3 gene. Fgfr3 gene harbors a G1120A point mutation in exon 10 causing G374R transition. G to A point mutation generated a SfcI restriction site (CTACAG) that can be identified by restriction enzyme SfcI. Primers P1 and P3 were used for SfcI cutting assay, amplicon from primers P1 and P2 were used for the Sanger sequence. C, PCR amplicon (374 bp) was cut with SfcI to identify the indel generation. Without any change, the fragment can be cut into 2 fragments of 198 and 176 bp each. 1–6, pups from injected embryos; 7, parental EIIa-Cre mouse; 8, parental Fgfr3G1120A-lnl/G1120A-lnl. D, sequence around the target from 6 pups revealed by Sanger sequencing. E, further Sanger sequence was conducted to verify the gene correction.
      Cas9 mRNA or protein, sgRNA2 RNA, and single strand donor WT DNA were co-injected into the pronuclear of Fgfr3G1120A zygotes, which were then implanted into pseudopregnant female mice for further development. Of a total of 63 live pups born, 42 mice were phenotypically normal (Table 1), suggesting the mutant allele was either disrupted or corrected.
      Table 1CRISPR/Cas9-mediated phenotype correction
      Experiment groupCas9 mRNA/proteinsgRNAssDO (200 nt)No. of embryos injectedNo. of pups bornNo. of normal pupsTargeted knockin
      ngl
      ImRNA (200 ng/μl)12.5184138 (61.6%)0/8
      IImRNA (5 ng/μl)2.55113150NA
      NA, not applicable.
      IIIProtein (30 ng/μl)2.5103953534 (97.1%)9/34 (26.5%)
      a NA, not applicable.
      We injected Cas9 mRNA at two concentrations, 200 and 5 ng/μl, which generated 13 and 15 pups, respectively. Of 13 pups generated by injection of 200 ng/μl mRNA, 8 (61.6%) where morphologically normal (Table 1). In contrast, none of the 15 mice generated by injection of 5 ng/μl of mRNA were normal, despite the increased concentrations of sgRNA and single strand DNA oligo (ssDO) by 2.5- and 2-fold, respectively, to increase genomic cutting and homologous recombination at the targeting site. This observation indicates that Cas9 mRNA at 5 ng/μl was too low to achieve effective genomic editing to the Fgfr3G1120A allele. Of note, 34 of 35 (97.1%) pups generated by injection of the Cas9 protein exhibited normal phenotype and 9 of them carry precision correction of the point mutation. Collectively, these data suggest that injection of Cas9 protein has a much higher efficiency in genome editing, although more conclusive data may be provided in future studies.

      High frequency of targeted correction of the Fgfr3G1120A point mutation by the injection of Cas9 protein in mice

      The correction of achondroplasia phenotype suggests that the Fgfr3G1120A allele in these mice was edited (either through targeted disruption or point mutation correction). Because the G to A transition at position 1120 (leading to G374R mutation) generates a SfcI restriction site (CTACAG) (Fig. 2B), thus change at the genomic level at this mutation site can be monitored by SfcI digestion. Therefore, to provide a quick analysis of the status of the mutation site, we used a pair of PCR primers flanking the G1120A mutation site to amplify a 374-bp fragment followed by restriction digestion with SfcI, which should cut this fragment into 176 and 198 bp if the CTACAG site is intact, whereas the loss of this site should render the SfcI resistance.
      Of 13 samples isolated from phenotype-corrected mice generated by the injection of Cas9 mRNA, 8 were SfcI resistant (pup 1 in Fig. 2, B and C), indicating the loss of cutting site; and 5 were sensitive (pups 2 and 3, Fig. 2, B and C), suggesting the G1120A mutation was not corrected. To understand the reason why these animals were normal while keeping the mutation, we conducted Sanger’s sequencing and the data indicated that pup 1 carried a 13-bp deletion, including the SfcI site, pup 2 had a deletion of “G,” and pup 3 had an insertion of “T,” respectively, in front of the SfcI site (Fig. 2D). In all three cases, the Fgfr3G1120A allele was effectively mutated due to the frameshift caused by indels, leading to the loss of the mutant phenotype.
      In 34 DNA samples isolated from injection of the Cas9 protein, 25 (73.5%) carried indels that disrupted the Fgfr3G374R allele and the animal appeared normal. In the remaining 9 mice (26.5%), the “A” was replaced by G (pups 5 and 6, Fig. 2, C–E) directed by ssDO template (Correction), all of which destroyed the SfcI site. These data indicate that injection of the Cas9 protein elicits high genome editing efficiency with a significant fraction generated by precise correction. Similar high efficiency was observed when a loxP site was introduced into the EGFP locus through homologous recombination (Fig. 3).
      Figure thumbnail gr3
      Figure 3Knock-in a loxP site into GFP transgene through pronuclear injection of CRISPR/Cas9 system. A, WT FVB mice were mated with homozygous GFP male mice to get pronuclear stage GFP-positive embryos for microinjection. B, Cas9 protein, sgRNA (CAACTACAAGACCCGCGCCG) targeting GFP gene, and a single strand donor DNA, which contains a loxp sequence and HindIII sequence flanked by two 100-nt homology arms for both sides were co-injected. C, after injection, we have checked 110 blastocysts for the fluoresce and knock-in event by using HindIII restriction enzyme to cut the insertion site after PCR amplification. There are 95 (86.4%) embryos that lost the GFP signal, whereas the remaining 15 (13.6%) still expressed the GFP gene. Among 95 GFP negative blastocysts, there are only 21 blastocyst containing a HindIII cutting site that indicates the ssDNA donor gets knocked-in successfully. D, represent picture of GFP negative blastocysts and fetus at embryonic day 14.

      Targeting single strand DNA oligos of 100 and 200 nt corrected Fgfr3G1120A mutation at equal efficiency in the achondroplasia mice

      The correction of the point mutation is primarily mediated by targeting ssDOs, which provide a template for repairing the double strand break created by the cutting of CRISPR/Cas9. The length of the ssDO used in the above experiment is 200 nt. Next, we investigated whether supplying a shorter ssDO would still maintain comparable efficiency of HDR-mediated precise genome repair. Our data indicated that the injection of a ssDO of 101 nt yielded a phenotype correction efficiency of 97.5% (40/41), which is comparable with 97.1% (34/35) obtained by using the ssDO of 200 nt. DNA sequencing revealed that 31 (77.5%) mice carried indels, whereas the remaining 9 (22.5%) mice had their G to A mutation precisely corrected (Table 2).
      Table 2Comparison of ssDO length for CRISPR/Cas9-mediated genomic editing
      ssDONo. of embryos injectedNo. of pups bornNo. of normal pupsNo. of pups (with indels)No. of pups (correction)
      200 nt3953534259 (26.5%)p = 0.86
      101 nt3714140319 (22.5%)
      Taken together, our data suggest that the Cas9 protein have higher efficiency to promote HDR-based genetic repairing, which was directed by exogenous oligos. 101 nt homologous ssDO keep the same efficiency for HDR compared with 200-nt ssDO.

      CRISPR/Cas9-mediated correction or disruption of Fgfr3G374R mutation completely suppresses phenotypes of achondroplasia mice

      The most significant features of achondroplasia mice include reduced body weight, dome-shaped head, short status, reduced lengths of long bone, most promptly in femurs and humerus (
      • Wang J.M.
      • Du. X.L.
      • Li C.L.
      • Yin L.J.
      • Chen B.
      • Sun J.
      • Su N.
      • Zhao L.
      • Song R.H.
      • Song W.W.
      • Chen L.
      • Deng C.X.
      Gly374Arg mutation in Fgfr3 causes achondroplasia in mice.
      ) (Fig. 4A). We carefully examined all 18 mice that carried the corrected point mutation allele using these features of achondroplasia, as measured by bone length and body weight, and the data indicated that all these were absent in the phenotype-corrected mice and they were indistinguishable with WT mice (Fig. 4, A–D). One major concern of the achondroplasia therapy is whether the dwarf phenotype can be inherited to the pups. We therefore mated them with WT mice and then analyzed the pups. Progeny obtained from corrected Fgfr3 point mutation mice were all normal. DNA sequencing using primer against the mutant allele showed that the pups carried the repaired Fgfr3 allele originating from their parents, indicating that the corrected allele could be successfully transmitted to the next generation through the germline.
      Figure thumbnail gr4
      Figure 4Observation of gene corrected and mutant mice. A, whole mount skeleton staining of Fgfr3G1120A/+ dwarf (left) and corrected mice (right) at 1 month of age. B, X-ray images and its statistical analysis (C) showing the length of femur of dwarf and corrected mice. D, comparison of body weight between dwarf mice, corrected mice, and knockout mice.
      The remaining 56 (25 from using of 200 nt ssDO and 31 from using 101-nt ssDO) mice that carry non-homologous end joining-mediated insertions or deletions were also phenotypically normal. Analysis of the nucleotide sequences of the Fgfr3 locus revealed that the mutant Fgfr3G374R allele was disrupted by indels at or near the G1120A mutation site in all these mice, whereas the remaining WT Fgfr3 allele is intact. This is consistent with our previous observation that mice heterozygous for the Fgfr3 knockout mutation is phenotypically normal (
      • Deng C.
      • Wynshaw-Boris A.
      • Zhou F.
      • Kuo A.
      • Leder P.
      Fibroblast growth factor receptor 3 is a negative regulator of bone growth.
      ).

      Correction of G1120A mutation without off-targeting events

      The CRISPR/Cas9 system sometimes suffers an off-target effect by cutting at sites with shared homology (
      • Martin F.
      • Sánchez-Hernandez S.
      • Gutiérrez-Guerrero A.
      • Pinedo-Gomez J.
      • Benabdellah K.
      Biased and unbiased methods for the detection of off-target cleavage by CRISPR/Cas9: an overview.
      ,
      • Rastogi A.
      • Murik O.
      • Bowler C.
      • Tirichine L.
      PhytoCRISP-Ex: a web-based and stand-alone application to find specific target sequences for CRISPR/CAS editing.
      • Haeussler M.
      • Schönig K.
      • Eckert H.
      • Eschstruth A.
      • Mianné J.
      • Renaud J.B.
      • Schneider-Maunoury S.
      • Shkumatava A.
      • Teboul L.
      • Kent J.
      • Joly J.S.
      • Concordet J.P.
      Evaluation of off-target and on-target scoring algorithms and integration into the guide RNA selection tool CRISPOR.
      ). To investigate this, we conducted whole-genome DNA sequencing for 7 mice delivered by one surrogate mother, including 3 mice with targeted mutation correction, 3 mice with indels (these 6 mice are phenotypically normal), and 1 mouse with un-corrected G1120A mutation (achondroplasia phenotype) revealed by SfcI digestion. As shown in Fig. 5, A and B, #1083 still contains the G to A mutation, whereas showing the dwarf phenotype. #1081, #1084, and #1089 only contained WT allele, suggesting the G1120A mutation was corrected. The remaining 3 mice, #1082, #1085, and #1086, exhibited mixed ratios of WT allele and indels. In the #1082 and #1086 mice, the WT allele was presented at 60 and 42%, respectively, which is considered the normal range of heterozygosity. Meanwhile, the #1082 mouse also contained a deletion of 17 bp and an insertion of a T, presented at 16 and 23%, respectively, suggesting this animal carries these two types of mutations in different portions of cells, presumably due to two distinct indels at 2- or 4-cell staged embryos. A similar situation was observed in #1086, which also carried insertion of a C and a T in different cells besides the WT allele. Mouse #1085 contained WT allele at 81.4% and an allele of 27 bp deletion at 18.6%, suggesting a precise correction of G1120A, and a deletion of 27 bp occurred in different cells at early 2-cell or 4-cell staged embryos. All these events (precision correction and indels) are mediated by CRISPR/Cas9 and effectively disrupted the Fgfr3G1120A mutation, resulting in the correction of achondroplasia phenotype of the mutant mice.
      Figure thumbnail gr5
      Figure 5Off-target analysis of sgRNA2 based on whole genome sequence. A, summary of genomic editing at the Fgfr3-G1120A site of 7 mice. #1083 still contains the Fgfr3-G1120A allele at about 40% suggesting no editing in the target region. #1081, #1084, and #1089 mice only contain WT allele, suggesting their Fgfr3-G1120A mutation was corrected. #1082 and #1086 mice contain WT allele with 60 and 42%, respectively, and also contain two different types of indels at lower percentages, suggesting these indels occurred at 2- or 4-cell staged embryos. #1085 mouse contain a WT allele at 81.4% and an allele of 27 bp deletion at 18.6%, suggesting the G1120A mutation may be corrected in some cells. B, types of mutations occurred in the Fgfr3-G1120A region of these 7 mice, including deletions of 17 or 27 bp, and two insertions of single nucleotides near the Fgfr3-G1120A mutation site, which effectively disrupted expression of the mutant allele.
      Next, we examined the entire Fgfr3 locus spanning 15.4 kb and detected no alterations. We then studied all homologous loci (Fig. 1C) for sgRNA2 with 1–4 mismatched nucleotides and detected no indels generated. These analyses indicate that the G1120A mutation was corrected specifically without obvious off-targeting events.

      Discussion

      CRISPR/Cas9 is a powerful tool for genome editing that has been used in multiple organisms (
      • Cong L.
      • Ran F.A.
      • Cox D.
      • Lin S.
      • Barretto R.
      • Habib N.
      • Hsu P.D.
      • Wu X.
      • Jiang W.
      • Marraffini L.A.
      • Zhang F.
      Multiplex genome engineering using CRISPR/Cas systems.
      ,
      • Mali P.
      • Yang L.
      • Esvelt K.M.
      • Aach J.
      • Guell M.
      • DiCarlo J.E.
      • Norville J.E.
      • Church G.M.
      RNA-guided human genome engineering via Cas9.
      ,
      • Chen S.
      • Sun H.
      • Miao K.
      • Deng C.X.
      CRISPR-Cas9: from genome editing to cancer res.
      ). However, CRISPR/Cas9, at its current stands, also suffers some limitations. For example, the CRISPR/Cas9 system frequently suffers an off-target effect by cutting at sites with shared homology (
      • Martin F.
      • Sánchez-Hernandez S.
      • Gutiérrez-Guerrero A.
      • Pinedo-Gomez J.
      • Benabdellah K.
      Biased and unbiased methods for the detection of off-target cleavage by CRISPR/Cas9: an overview.
      ,
      • Rastogi A.
      • Murik O.
      • Bowler C.
      • Tirichine L.
      PhytoCRISP-Ex: a web-based and stand-alone application to find specific target sequences for CRISPR/CAS editing.
      • Haeussler M.
      • Schönig K.
      • Eckert H.
      • Eschstruth A.
      • Mianné J.
      • Renaud J.B.
      • Schneider-Maunoury S.
      • Shkumatava A.
      • Teboul L.
      • Kent J.
      • Joly J.S.
      • Concordet J.P.
      Evaluation of off-target and on-target scoring algorithms and integration into the guide RNA selection tool CRISPOR.
      ). Furthermore, due to its high efficiency cutting at targeting sites, CRISPR/Cas9 preferably generates more indels than precise modifications, casting difficulty to generate the desired mutations. Using Fgfr3G374R achondroplasia mice as a model system, we have addressed some of the potential problems. Our optimized experiment conditions have achieved high specificity of targeted correction of achondroplasia with several navel features.
      Unintended binding, modification, and cleavage of nucleic acids, so called the off-target effect, is a major challenge to the CRISPR/Cas9 system (
      • Knott G.J.
      • Doudna J.A.
      CRISPR-Cas guides the future of genetic engineering.
      ). To overcome this problem, many sgRNA designing tools have been developed to increase specificity of sgRNAs at the silicon level (
      • Cui Y.
      • Xu J.
      • Cheng M.
      • Liao X.
      • Peng S.
      Review of CRISPR/Cas9 sgRNA design tools.
      ,
      • Listgarten J.
      • Weinstein M.
      • Kleinstiver B.P.
      • Sousa A.A.
      • Joung J.K.
      • Crawford J.
      • Gao K.
      • Hoang L.
      • Elibol M.
      • Doench J.G.
      • Fusi N.
      Prediction of off-target activities for the end-to-end design of CRISPR guide RNAs.
      ). A number of methods have also been employed to avoid off-target cutting or identify undesired mutations, such as using single-based editors correct base-pairing to avoid inducing a dsDNA break, or replace Cas9 with Cas13a to target RNA editing (
      • Knott G.J.
      • Doudna J.A.
      CRISPR-Cas guides the future of genetic engineering.
      ). In our study, the off-target effect still occurred even if CRISPR design software was used to design our sgRNAs. Because of the newly designed GFP tester, we found two of three sgRNAs could cut the WT Fgfr3 sequence with a single mismatched base, whereas the other sgRNA cuts only at the MT Fgfr3 sequence. We used this sgRNA for further genome editing in mice and it indeed did not cut the corresponding WT Fgfr3 sequence with single mismatched base.
      More recently, there are some discussions about whether or not the CRISPR/Cas9 system could generate unexpected mutations in vivo (
      • Lareau C.A.
      • Clement K.
      • Hsu J.Y.
      • Pattanayak V.
      • Joung J.K.
      • Aryee M.J.
      • Pinello L.
      Response to “unexpected mutations after CRISPR-Cas9 editing in vivo.”.
      • Nutter L.M.J.
      • Heaney J.D.
      • Lloyd K.C.K.
      • Murray S.A.
      • Seavitt J.R.
      • Skarnes W.C.
      • Teboul L.
      • Brown S.D.M.
      • Moore M.
      Response to “unexpected mutations after CRISPR-Cas9 editing in vivo.”.
      ,
      • Kim S.T.
      • Park J.
      • Kim D.
      • Kim K.
      • Bae S.
      • Schlesner M.
      • Kim J.S.
      Response to “unexpected mutations after CRISPR-Cas9 editing in vivo.”.
      ,
      • Wilson C.J.
      • Fennell T.
      • Bothmer A.
      • Maeder M.L.
      • Reyon D.
      • Cotta-Ramusino C.
      • Fernandez C.A.
      • Marco E.
      • Barrera L.A.
      • Jayaram H.
      • Albright C.F.
      • Cox G.F.
      • Church G.M.
      • Myer V.E.
      Response to “unexpected mutations after CRISPR-Cas9 editing in vivo.”.
      • Lescarbeau R.M.
      • Murray B.
      • Barnes T.M.
      • Bermingham N.
      Response to “unexpected mutations after CRISPR-Cas9 editing in vivo.”.
      ). To investigate this, we conducted whole genome sequencing, and detected no off-target events in the region near the targeting site, the entire Fgfr3 locus, and in all predicted potential off-target sites, which share varying degrees of homology with the target sequence. Altogether, we demonstrated that by using our optimized sgRNA, the CRISPR/Cas9-mediated targeted correction of the Fgfr3-G374R point mutation completely suppresses phenotypes of achondroplasia mice with high fidelity. This result demonstrates in principle the importance of using such a tester system to avoid potential off-target effects prior to the experiment. For example, if sgRNA 1 or 3 was used in our experiment, it might have generated off-target events. The GFP tester is user friendly, and reliable for providing functional evaluation of sgRNAs for their specificity and efficiency. From the long run, highly suspected sequences for off-target predicted by software could also be tested first to avoid any potential problems.
      Both Cas9 mRNA and proteins have been widely used for achieving genome editing with variable frequencies (
      • Chu V.T.
      • Weber T.
      • Graf R.
      • Sommermann T.
      • Petsch K.
      • Sack U.
      • Volchkov P.
      • Rajewsky K.
      • Kühn R.
      Efficient generation of Rosa26 knock-in mice using CRISPR/Cas9 in C57BL/6 zygotes.
      ,
      • Wang L.
      • Shao Y.
      • Guan Y.
      • Li L.
      • Wu L.
      • Chen F.
      • Liu M.
      • Chen H.
      • Ma Y.
      • Ma X.
      • Liu M.
      • Li D.
      Large genomic fragment deletion and functional gene cassette knock-in via Cas9 protein mediated genome editing in one-cell rodent embryos.
      ). Our data revealed significant higher frequency of CRRSPR/Cas9-mediated genome cutting at the target site and the precise correction of the G1120A point mutation by the Cas9 protein than Cas9 mRNA. For the underlying mechanism, we believe that there will be a gap between the time of injection of mRNA and producing the effective amount of protein by translation. However, this gap no longer exists by direct admission of Cas9 protein, as the protein can be immediately used for genome editing. We have also tested the knock-in efficiency by Cas9 protein to EGFP genes, and a comparable efficiency was achieved by introducing a 47-base loxP site into the EGFP gene. In conclusion, Cas9 protein achieves much higher efficiency than Cas9 mRNA of genome editing and introduction of knock-in mutation.
      Conventional gene targeting by homologous recombination has high demand on the length of targeting homology to achieve ideal targeting efficiency (
      • Deng C.
      • Capecchi M.R.
      Reexamination of gene targeting frequency as a function of the extent of homology between the targeting vector and the target locus.
      ). Because the CRISPR/Cas9 system generates double strand breaks at the target sites that greatly stimulates efficiency of homologous recombination, the long homologous arm of targeting constructs is less critical (
      • Deng C.
      • Capecchi M.R.
      Reexamination of gene targeting frequency as a function of the extent of homology between the targeting vector and the target locus.
      ). Donor oligo or targeting constructs with less than 1000-bp homology arms are frequently used for genome editing in mouse zygotes (
      • Yao X.
      • Zhang M.
      • Wang X.
      • Ying W.
      • Hu X.
      • Dai P.
      • Meng F.
      • Shi L.
      • Sun Y.
      • Yao N.
      • Zhong W.
      • Li Y.
      • Wu K.L.
      • Li W.P.
      • Chen Z.
      • Yang H.
      Tild-CRISPR allows for efficient and precise gene knockin in mouse and human cells.
      ). In our study, we first tested targeting ssDO with 200 nt (with around 100 nt on each side of the G to A mutation) and obtained genome editing frequency in 34/35 (97%) mice. In the 34 mice, 9 carry precise G1120A mutation correction (26.5%). To further test the effect of the length of homology, we injected a 101-nt ssDO and obtained comparable editing frequency and mutation correction rate. These data demonstrate that the targeting oligos of 101 and 200 nt precisely corrected the mutation at similar efficiency. As oligos of 101 nt can be generated much more economically, our finding is of great significance when designing targeting ssODs.
      In summary, we have demonstrated that the CRISPR/Cas9 system can be used to cure achondroplasia in mouse by directly correcting the genetic defect through homology-mediated gene editing. Our GFP reporter can be used to select sgRNAs with high specificity and reasonable cutting efficiency to avoid off-target effects. Microinjection of Cas9 protein greatly improves genome editing frequency. Targeting oligo of 101 nt, which is more easily affordable than longer oligos, can be reliably used for the editing. The combination of these factors greatly facilitates the introduction of the desired genome modification with high fidelity.

      Experimental procedures

      Animal model

      As previously described to bypass the lethal effect of Fgfr3G1120A mutation, a lnl cassette is inserted between exons 10 and 11 to knockout the mutant Fgfr3 gene (
      • Wang J.M.
      • Du. X.L.
      • Li C.L.
      • Yin L.J.
      • Chen B.
      • Sun J.
      • Su N.
      • Zhao L.
      • Song R.H.
      • Song W.W.
      • Chen L.
      • Deng C.X.
      Gly374Arg mutation in Fgfr3 causes achondroplasia in mice.
      ). Mice harboring homozygous Fgfr3G1120A-lnl/G1120A-lnl exhibit phenotypes similar to Fgfr3−/− mice (
      • Deng C.
      • Wynshaw-Boris A.
      • Zhou F.
      • Kuo A.
      • Leder P.
      Fibroblast growth factor receptor 3 is a negative regulator of bone growth.
      ) and are maintained to serve as stud mice. Homozygote of EIIa-Cre (JAX, ME) female mice were superovulated by injection with pregnant mare’s serum gonadotropin followed by human chorionic gonadotropin, and then crossed with Fgfr3G1120A-lnl/G1120A-lnl stud mice to delete the lnl cassette and obtain Fgfr3 point mutant embryos for pronuclear injection. All animal procedures were performed under the ethical guidelines of the University of Macau (animal protocol number: UMAEC-037-2015).

      sgRNA reporter plasmid construction

      sgRNA reporter plasmid (pCAG-EGx-xFP-BbsI) was constructed by modifying pCAG-EGxxFP (Addgene). Briefly, a short dual BbsI cassette was amplified from plasmid PX330 (
      • Cong L.
      • Ran F.A.
      • Cox D.
      • Lin S.
      • Barretto R.
      • Habib N.
      • Hsu P.D.
      • Wu X.
      • Jiang W.
      • Marraffini L.A.
      • Zhang F.
      Multiplex genome engineering using CRISPR/Cas systems.
      ) and inserted into the pCAG-EGxxFP digested with EcoRI to obtain pCAG-EGxxFP-BbsI (Fig. 1A). sgRNA tester sequence was designed based on the sgRNA sequence including the NGG PAM sequence, then the adapter was added for further cloning. To clone the target sequence into the reporter plasmid, paired oligos with adaptor were synthesized (Fig. 1A). Each pair of oligos was phosphorylated and annealed followed by ligation with a BbsI pre-digested reporter plasmid.

      Microinjection

      sgRNA(s) were designed by using Optimized CRISPR Design online tools (
      • Hsu P.D.
      • Scott D.A.
      • Weinstein J.A.
      • Ran F.A.
      • Konermann S.
      • Agarwala V.
      • Li Y.
      • Fine E.J.
      • Wu X.
      • Shalem O.
      • Cradick T.J.
      • Marraffini L.A.
      • Bao G.
      • Zhang F.
      DNA targeting specificity of RNA-guided Cas9 nucleases.
      ), then ligation with PX330 plasmid. sgRNAs with a T7 promoter were amplified by PCR and in vitro transcribed using a MEGAshortscript T7 kit (Thermo Fisher Scientific, MA). After transcription, the sgRNAs were purified with a MEGAclear kit (Thermo Fisher Scientific, MA) according to the manufacturer's instructions. The sequences used for preparation of template for in vitro transcription of sgRNA2 was amplified from PX330-sgRNA2 using primers (TTAATACGACTCACTATAGGCGCAGGCGTCCTCAGCTAC and AAAAGCACCGACTCGGTGCC).
      Two single strand donor DNA oligo (ssDO) of the total length of 200 and 101 nt were designed based on WT Fgfr3 genomic DNA. The ssDO-200 nt contains 99 and 100 nt at 5′ and 3′ to the G to A point mutation, respectively; and the ssODN-101 nt contains 50 nt on each side flanking the G to A point mutation as shown; Wt-200nt, AAGGCTGGATGAGGCCCCAAAATTTGTATCTTTGCAGCTGAGGAGGAGCTGATGGAAACTGATGAGGCTGGCAGCGTGTACGCAGGCGTCCTCAGCTACGGGGTGGTCTTCTTCCTCTTCATCCTGGTGGTGGCAGCTGTGATACTCTGCCGCCTGCGCAGTCCCCCAAAGAAGGGCTTGGGCTCGCCCACCGTGCACAA; and Wt-101nt: TGATGGAAACTGATGAGGCTGGCAGCGTGTACGCAGGCGTCCTCAGCTACGGGGTGGTCTTCTTCCTCTTCATCCTGGTGGTGGCAGCTGTGATACTCTGC.
      EIIa-Cre female mice were superovulated and mated with Fgfr3G1120A-lnl/G1120A-lnl males. Pronuclear stage embryos were collected followed by microinjection punched by FemtoJet device (Eppendorf). Cas9 mRNA or proteins were co-injected with sgRNA and ssDO. For Cas9 mRNA injection, different concentrations of Cas9 mRNAs (5 or 200 ng/μl) and sgRNA (2.5 or 1 ng/μl) were mixed with ssDO (5 or 2.5 ng/μl) injected into the embryos. The injected embryos were transferred into pseudopregnant female mice.

      Genotyping, DNA sequencing analysis, and Sanger sequencing

      The target site (sequence around Fgfr3G1120A mutation site) was amplified from genomic DNA of the genetic modified mice by PCR with primers (Fgfr3-s and Fgfr3-a) for SfcI digestion diagnosis, and flanking primers (Fgfr3-s and Loxp-a) for Sanger sequence; Fgfr3-s (P1), CTCTTCTCCAAGTATCCCAGGTCC; Fgfr3-a (P2), CCTGCTGGGACTCTAGGAGACAC; and Loxp-a (P3), CGAAGTTATCTAGAGTCGACCATCG.

      Whole genome sequencing

      Total DNA was extracted from mouse tail using DNeasy Blood & Tissue Kits (Qiagen, Hilden, Germany). Whole genome sequencing libraries were prepared using the DNA Library Prep kit (New England Biolabs) in accordance with the manufacturer’s instructions. Briefly, 1 μg of DNA was sheared using a Covarias sonicator, then followed by end-repair, ligation, and amplification. Whole genome sequencing libraries were evaluated using BioAnalyzer (Agilent, CA) and quantitative PCR. Libraries were sequenced on an Illumina Hiseq ×10 sequencer using 2 × 150 bp cycles to meet coverage of 30 times for further analysis.

      Bioinformatics analyses

      Sequenced reads were mapped to the mouse (Mus musculus) genome (mm10) using Burrows-Wheeler Aligner (version 0.7.15). Unique mapped reads were used for the further analysis. GATK best practices (version 4.0.3) workflows were used to identify SNV(s) and Indel(s). Variants were further filtered to exclude variants in dbSNPv150. Variant annotations were done via SnpEff (version 4.3T) The top potential off-target sites were predicted by Optimized CRISPR Design online tools (
      • Hsu P.D.
      • Scott D.A.
      • Weinstein J.A.
      • Ran F.A.
      • Konermann S.
      • Agarwala V.
      • Li Y.
      • Fine E.J.
      • Wu X.
      • Shalem O.
      • Cradick T.J.
      • Marraffini L.A.
      • Bao G.
      • Zhang F.
      DNA targeting specificity of RNA-guided Cas9 nucleases.
      ).

      X-ray and whole mount skeletal preparation

      For X-ray imaging, animals were euthanatized with CO2. After removing the skin, photographs were taken by using an X-ray machine (Bruker, MA). For skeletal preparation, the carcasses were eviscerated, fixed in 95% ethanol, stained with Alizarin red S and Alcian blue, cleared by KOH treatment, and then stored in glycerol as previously described (
      • Deng C.
      • Wynshaw-Boris A.
      • Zhou F.
      • Kuo A.
      • Leder P.
      Fibroblast growth factor receptor 3 is a negative regulator of bone growth.
      ).

      Author contributions

      K. M., X. Z., J. Z., X. X., and C.-X. D. resources; K. M., X. Z., J. Z., X. X., and C.-X. D. data curation; K. M., X. Z., and J. Z. software; K. M., X. X., and C.-X. D. funding acquisition; K. M., X. Z., and S. M. S. validation; K. M., X. Z., S. M. S., Z. H., U. I. C., X. X., and C.-X. D. investigation; K. M., X. Z., and S. M. S. visualization; K. M., X. Z., S. M. S., J. Z., Z. H., U. I.C., X. X., and C.-X. D. methodology; K. M., X. X., and C.-X. D. writing-original draft; K. M. and C.-X. D. writing-review and editing; X. Z., J. Z., X. X., and C.-X. D. formal analysis; X. X. and C.-X. D. conceptualization; X. X. and C.-X. D. supervision; X. X. and C.-X. D. project administration.

      Acknowledgments

      We thank members of the C. Deng and X. Xu laboratories for helpful advice and discussion.

      Supplementary Material

      References

        • Cong L.
        • Ran F.A.
        • Cox D.
        • Lin S.
        • Barretto R.
        • Habib N.
        • Hsu P.D.
        • Wu X.
        • Jiang W.
        • Marraffini L.A.
        • Zhang F.
        Multiplex genome engineering using CRISPR/Cas systems.
        Science. 2013; 339 (23287718): 819-823
        • Mali P.
        • Yang L.
        • Esvelt K.M.
        • Aach J.
        • Guell M.
        • DiCarlo J.E.
        • Norville J.E.
        • Church G.M.
        RNA-guided human genome engineering via Cas9.
        Science. 2013; 339 (23287722): 823-826
        • Knott G.J.
        • Doudna J.A.
        CRISPR-Cas guides the future of genetic engineering.
        Science. 2018; 361 (30166482): 866-869
        • Chen S.
        • Sun H.
        • Miao K.
        • Deng C.X.
        CRISPR-Cas9: from genome editing to cancer res.
        Int. J. Biol. Sci. 2016; 12 (27994508): 1427-1436
        • Luther D.C.
        • Lee Y.W.
        • Nagaraj H.
        • Scaletti F.
        • Rotello V.M.
        Delivery approaches for CRISPR/Cas9 therapeutics in vivo: advances and challenges.
        Expert Opin. Drug Deliv. 2018; 15: 905-913
        • Chen S.
        • Sanjana N.E.
        • Zheng K.
        • Shalem O.
        • Lee K.
        • Shi X.
        • Scott D.A.
        • Song J.
        • Pan J.Q.
        • Weissleder R.
        • Lee H.
        • Zhang F.
        • Sharp P.A.
        Genome-wide CRISPR screen in a mouse model of tumor growth and metastasis.
        Cell. 2015; 160 (25748654): 1246-1260
        • Wu Y.
        • Liang D.
        • Wang Y.
        • Bai M.
        • Tang W.
        • Bao S.
        • Yan Z.
        • Li D.
        • Li J.
        Correction of a genetic disease in mouse via use of CRISPR-Cas9.
        Cell Stem Cell. 2013; 13 (24315440): 659-662
        • Yin H.
        • Xue W.
        • Chen S.
        • Bogorad R.L.
        • Benedetti E.
        • Grompe M.
        • Koteliansky V.
        • Sharp P.A.
        • Jacks T.
        • Anderson D.G.
        Genome editing with Cas9 in adult mice corrects a disease mutation and phenotype.
        Nat. Biotechnol. 2014; 32 (24681508): 551-553
        • Ou Z.
        • Niu X.
        • He W.
        • Chen Y.
        • Song B.
        • Xian Y.
        • Fan D.
        • Tang D.
        • Sun X.
        The combination of CRISPR/Cas9 and iPSC technologies in the gene therapy of human β-thalassemia in mice.
        Sci. Rep. 2016; 6 (27581487): 32463
        • Tabebordbar M.
        • Zhu K.
        • Cheng J.K.W.
        • Chew W.L.
        • Widrick J.J.
        • Yan W.X.
        • Maesner C.
        • Wu E.Y.
        • Xiao R.
        • Ran F.A.
        • Cong L.
        • Zhang F.
        • Vandenberghe L.H.
        • Church G.M.
        • Wagers A.J.
        In vivo gene editing in dystrophic mouse muscle and muscle stem cells.
        Science. 2016; 351 (26721686): 407-411
        • Nelson C.E.
        • Hakim C.H.
        • Ousterout D.G.
        • Thakore P.I.
        • Moreb E.A.
        • Castellanos Rivera R.M.
        • Madhavan S.
        • Pan X.
        • Ran F.A.
        • Yan W.X.
        • Asokan A.
        • Zhang F.
        • Duan D.
        • Gersbach C.A.
        In vivo genome editing improves muscle function in a mouse model of Duchenne muscular dystrophy.
        Science. 2016; 351 (26721684): 403-407
        • Long C.
        • Amoasii L.
        • Mireault A.A.
        • McAnally J.R.
        • Li H.
        • Sanchez-Ortiz E.
        • Bhattacharyya S.
        • Shelton J.M.
        • Bassel-Duby R.
        • Olson E.N.
        Postnatal genome editing partially restores dystrophin expression in a mouse model of muscular dystrophy.
        Science. 2016; 351 (26721683): 400-403
        • Yang Y.
        • Wang L.
        • Bell P.
        • McMenamin D.
        • He Z.
        • White J.
        • Yu H.
        • Xu C.
        • Morizono H.
        • Musunuru K.
        • Batshaw M.L.
        • Wilson J.M.
        A dual AAV system enables the Cas9-mediated correction of a metabolic liver disease in newborn mice.
        Nat. Biotechnol. 2016; 34 (26829317): 334-338
        • Yin H.
        • Song C.Q.
        • Dorkin J.R.
        • Zhu L.J.
        • Li Y.
        • Wu Q.
        • Park A.
        • Yang J.
        • Suresh S.
        • Bizhanova A.
        • Gupta A.
        • Bolukbasi M.F.
        • Walsh S.
        • Bogorad R.L.
        • Gao G.
        • Weng Z.
        • Dong Y.
        • Koteliansky V.
        • Wolfe S.A.
        • Langer R.
        • Xue W.
        • Anderson D.G.
        Therapeutic genome editing by combined viral and non-viral delivery of CRISPR system components in vivo.
        Nat. Biotechnol. 2016; 34 (26829318): 328-333
        • Martin F.
        • Sánchez-Hernandez S.
        • Gutiérrez-Guerrero A.
        • Pinedo-Gomez J.
        • Benabdellah K.
        Biased and unbiased methods for the detection of off-target cleavage by CRISPR/Cas9: an overview.
        Int. J. Mol. Sci. 2016; 17 (27618019): e1507
        • Rastogi A.
        • Murik O.
        • Bowler C.
        • Tirichine L.
        PhytoCRISP-Ex: a web-based and stand-alone application to find specific target sequences for CRISPR/CAS editing.
        BMC Bioinformatics. 2016; 17 (27363443): 261
        • Haeussler M.
        • Schönig K.
        • Eckert H.
        • Eschstruth A.
        • Mianné J.
        • Renaud J.B.
        • Schneider-Maunoury S.
        • Shkumatava A.
        • Teboul L.
        • Kent J.
        • Joly J.S.
        • Concordet J.P.
        Evaluation of off-target and on-target scoring algorithms and integration into the guide RNA selection tool CRISPOR.
        Genome Biol. 2016; 17 (27380939): 148
        • Ornitz D.M.
        • Legeai-Mallet L.
        Achondroplasia: development, pathogenesis, and therapy.
        Dev. Dyn. 2017; 246 (27987249): 291-309
        • Bellus G.A.
        • Hefferon T.W.
        • Ortiz de Luna R.I.
        • Hecht J.T.
        • Horton W.A.
        • Machado M.
        • Kaitila I.
        • McIntosh I.
        • Francomano C.A.
        Achondroplasia is defined by recurrent G380R mutations of FGFR3.
        Am. J. Hum. Genet. 1995; 56 (7847369): 368-373
        • Stoilov I.
        • Kilpatrick M.W.
        • Tsipouras P.
        A common FGFR3 gene mutation is present in achondroplasia but not in hypochondroplasia.
        Am. J. Med. Genet. 1995; 55 (7702086): 127-133
        • Wang J.M.
        • Du. X.L.
        • Li C.L.
        • Yin L.J.
        • Chen B.
        • Sun J.
        • Su N.
        • Zhao L.
        • Song R.H.
        • Song W.W.
        • Chen L.
        • Deng C.X.
        Gly374Arg mutation in Fgfr3 causes achondroplasia in mice.
        Chinese J. Med. Genet. 2004; 21 (15583977): 537
        • Deng C.
        • Wynshaw-Boris A.
        • Zhou F.
        • Kuo A.
        • Leder P.
        Fibroblast growth factor receptor 3 is a negative regulator of bone growth.
        Cell. 1996; 84 (8601314): 911-921
        • Cui Y.
        • Xu J.
        • Cheng M.
        • Liao X.
        • Peng S.
        Review of CRISPR/Cas9 sgRNA design tools.
        Interdiscip. Sci. 2018; 10 (29644494): 455-465
        • Listgarten J.
        • Weinstein M.
        • Kleinstiver B.P.
        • Sousa A.A.
        • Joung J.K.
        • Crawford J.
        • Gao K.
        • Hoang L.
        • Elibol M.
        • Doench J.G.
        • Fusi N.
        Prediction of off-target activities for the end-to-end design of CRISPR guide RNAs.
        Nat. Biomed. Eng. 2018; 2 (29998038): 38-47
        • Lareau C.A.
        • Clement K.
        • Hsu J.Y.
        • Pattanayak V.
        • Joung J.K.
        • Aryee M.J.
        • Pinello L.
        Response to “unexpected mutations after CRISPR-Cas9 editing in vivo.”.
        Nat. Methods. 2018; 15 (29600992): 238-239
        • Nutter L.M.J.
        • Heaney J.D.
        • Lloyd K.C.K.
        • Murray S.A.
        • Seavitt J.R.
        • Skarnes W.C.
        • Teboul L.
        • Brown S.D.M.
        • Moore M.
        Response to “unexpected mutations after CRISPR-Cas9 editing in vivo.”.
        Nat. Methods. 2018; 15 (29600991): 235-236
        • Kim S.T.
        • Park J.
        • Kim D.
        • Kim K.
        • Bae S.
        • Schlesner M.
        • Kim J.S.
        Response to “unexpected mutations after CRISPR-Cas9 editing in vivo.”.
        Nat. Methods. 2018; 15 (29600990): 239-240
        • Wilson C.J.
        • Fennell T.
        • Bothmer A.
        • Maeder M.L.
        • Reyon D.
        • Cotta-Ramusino C.
        • Fernandez C.A.
        • Marco E.
        • Barrera L.A.
        • Jayaram H.
        • Albright C.F.
        • Cox G.F.
        • Church G.M.
        • Myer V.E.
        Response to “unexpected mutations after CRISPR-Cas9 editing in vivo.”.
        Nat. Methods. 2018; 15 (29600989): 236-237
        • Lescarbeau R.M.
        • Murray B.
        • Barnes T.M.
        • Bermingham N.
        Response to “unexpected mutations after CRISPR-Cas9 editing in vivo.”.
        Nat. Methods. 2018; 15 (29600988): 237
        • Chu V.T.
        • Weber T.
        • Graf R.
        • Sommermann T.
        • Petsch K.
        • Sack U.
        • Volchkov P.
        • Rajewsky K.
        • Kühn R.
        Efficient generation of Rosa26 knock-in mice using CRISPR/Cas9 in C57BL/6 zygotes.
        BMC Biotechnol. 2016; 16 (26772810): 4
        • Wang L.
        • Shao Y.
        • Guan Y.
        • Li L.
        • Wu L.
        • Chen F.
        • Liu M.
        • Chen H.
        • Ma Y.
        • Ma X.
        • Liu M.
        • Li D.
        Large genomic fragment deletion and functional gene cassette knock-in via Cas9 protein mediated genome editing in one-cell rodent embryos.
        Sci. Rep. 2015; 5 (26620761): 17517
        • Deng C.
        • Capecchi M.R.
        Reexamination of gene targeting frequency as a function of the extent of homology between the targeting vector and the target locus.
        Mol. Cell. Biol. 1992; 12 (1321331): 3365-3371
        • Yao X.
        • Zhang M.
        • Wang X.
        • Ying W.
        • Hu X.
        • Dai P.
        • Meng F.
        • Shi L.
        • Sun Y.
        • Yao N.
        • Zhong W.
        • Li Y.
        • Wu K.L.
        • Li W.P.
        • Chen Z.
        • Yang H.
        Tild-CRISPR allows for efficient and precise gene knockin in mouse and human cells.
        Dev. Cell. 2018; 45 (29787711): 526-536
        • Hsu P.D.
        • Scott D.A.
        • Weinstein J.A.
        • Ran F.A.
        • Konermann S.
        • Agarwala V.
        • Li Y.
        • Fine E.J.
        • Wu X.
        • Shalem O.
        • Cradick T.J.
        • Marraffini L.A.
        • Bao G.
        • Zhang F.
        DNA targeting specificity of RNA-guided Cas9 nucleases.
        Nat. Biotechnol. 2013; 31 (23873081): 827-832