Skip to main content

Advertisement

CRISPR/Cas9-mediated genome editing induces exon skipping by alternative splicing or exon deletion

Genome Biology volume 18, Article number: 108 (2017) Cite this article

Abstract

CRISPR is widely used to disrupt gene function by inducing small insertions and deletions. Here, we show that some single-guide RNAs (sgRNAs) can induce exon skipping or large genomic deletions that delete exons. For example, CRISPR-mediated editing of β-catenin exon 3, which encodes an autoinhibitory domain, induces partial skipping of the in-frame exon and nuclear accumulation of β-catenin. A single sgRNA can induce small insertions or deletions that partially alter splicing or unexpected larger deletions that remove exons. Exon skipping adds to the unexpected outcomes that must be accounted for, and perhaps taken advantage of, in CRISPR experiments.

Background

CRISPR/Cas9 genome editing has transformed the study of gene function in many organisms [1,2,3,4,5]. Guide RNAs direct the Cas9 nuclease to create double-strand DNA breaks at complementary target sites in the genome. Repair of these double-strand DNA breaks by non-homologous end-joining (NHEJ) often introduces small insertions or deletions (indels) that shift the open reading frame, thereby inactivating the target gene. CRISPR therefore provides a simple way to generate loss-of-function (LOF) mutations in virtually any gene in the mammalian genome [1]. Nonetheless, CRISPR can also induce off-target editing at genomic positions that imperfectly match the single-guide RNA (sgRNA) sequence, which calls for the implementation of strategies to reduce off-target effects [6, 7]. Besides off-target editing, it remains unknown whether CRISPR-mediated editing has unintended consequence at the post-transcriptional level of the target gene.

We have previously used in vivo delivery of CRISPR to inactivate tumor suppressor genes in mice [8,9,10]. We also showed that CRISPR can edit oncogenes or disease genes through homolog-directed repair [8, 11, 12]. Here we show that CRISPR-mediated editing of mammalian exons can induce exon skipping. Exon skipping can result from alternative splicing or from genomic deletions that remove exons. Moreover, exon skipping can produce messenger RNAs (mRNAs) with intact reading frames that encode functional proteins.

Results

We recently used CRISPR to disrupt the Kras oncogene in two independent lung adenocarcinoma cell lines [13], which were derived from Kras G12D; p53 fl/fl (KP) mice [14, 15]. We isolated two single-cell clones each carrying frameshifting deletions in exon 2 (Fig. 1a and Additional file 1: Figure S1a): KP1 carries a 2-nt “-CG” deletion in the G12D allele and a 1-nt “-C” deletion in the otherwise wild-type (WT) Kras allele; and KP2 carries a 2-nt “-GG” deletion. Neither clone produces full length Kras protein [13], indicating that all three deletions disrupt the Kras reading frame.

Fig. 1
figure1

sgRNA targeting Kras induces exon skipping in single cell clones. a Schematic of an sgRNA targeting exon 2 of the mouse Kras gene (sgKras). The red arrowhead denotes the Cas9 cleavage site. KP1 and KP2 cell lines were transduced with lentivirus that encodes Cas9 and sgKras. Two single-cell clones (KP1 clone and KP2 clone) harbor frameshift deletions. Black arrows indicate the positions of reverse transcription polymerase chain reaction (RT-PCR) primers. The G12D codon is underlined. b Normalized Kras read counts from RNA-sequencing (RNA-seq) analysis of KP parental cells (blue) and KP clones (red). RNA-seq was done twice for KP2 clone and three times for the other groups. “+” denotes WT allele. c RNA-seq showing partial exon 2 skipping in KP1 clones. RNA-seq numbers indicate reads spanning the indicated exon junctions. Two representative biological replicates are shown. d RT-PCR analysis of Kras mRNA detects an exon 2 skipped band. The expected band sizes are 331 bp and 209 bp. M, molecular marker. “*” denotes indels in PCR products from clones. e Scatter plot showing 22 exon events that change in both KP1 and KP2 clones. Exclusion of Kras exon 2 is the most frequent event. Ψ, Percentage Splicing Index

Full size image

Frameshift mutations in early exons are known to trigger nonsense-mediated decay (NMD) [16], which eliminates mRNAs with premature termination codons. When we analyzed mRNA-sequencing (RNA-seq) data, however, we found that apparent Kras mRNA levels (i.e. total normalized mRNA reads) were only reduced by 19% in KP1 cells and 47% in KP2 cells, compared with parental KP cells (Fig. 1b). Both clones produced fewer exon 2 reads, but normal levels of exon 1 and 3 reads (Fig. 1c), suggesting that exon 2 might be skipped in the KP1 and KP2 clones. Indeed, we detected exon 1-3 junction reads, indicating that exon 2 was skipped (Fig. 1c and Additional file 1: Figure S1b). Calculating the ratio between exon 2 reads and total reads, we found that exon 2 is included in only 64.0 ± 9.1% of Kras reads from KP1-clone (Fig. 1c and Table 1). Similar exon 2 skipping was observed in KP2-clones (Additional file 1: Figure S1c). Concordantly, reverse transcription of Kras mRNA followed by polymerase chain reaction (RT-PCR) yielded two products: one corresponded to intact Kras complementary DNA (cDNA) and the other corresponded to the exon 1-3 isoform (Fig. 1d). The exon 1-3 isoform retains a partial Kras open reading frame that could initiate translation from an ATG codon in exon 3 (Additional file 1: Figure S2) and produce a severely truncated Kras protein.

Table 1 Genomic lesion and mRNA splicing results of single cell clones
Full size table

Editing of Kras did not induce alternative splicing genome-wide. We identified 97 alternatively spliced exons in KP1 cells and 177 events in KP2 cells. KP1 and KP2 clones shared 22 cassette inclusion or exclusion events, with the exclusion of Kras exon 2 being the greatest change in both clones (Fig. 1e and Additional file 1: Table S3). Thus, editing of Kras exon 2 specifically induced skipping of Kras exon 2. Notably, whereas mouse Kras G12D (GGU to GAU) transcripts do not skip exon 2 in parental KP cells, we found that ~15% of human KRAS G12S (codon 12 GGU to AGU) transcripts skip exon 2 in the A549 human lung cancer cell line (Additional file 1: Figure S1d). We were unable to predict the gain or loss of exon splice enhancers or silencers [17], but our data suggest that sequences near Kras codon 12 promote exon 2 inclusion in mouse and human Kras. Exon skipping induced by CRISPR editing was not limited to Kras or to mouse KP cells. A recent study showed that CRISPR editing of FLOT1 exon 3 in HeLa cells can cause skipping of exon 3, exon 4, or exons 3, 4, and 5 [18]. We also detected infrequent exon skipping when we targeted exon 11 of LMNA in human HCT116 cells (Additional file 1: Figure S3). Skipping LMNA exon 11 produces an in-frame transcript that could be translated into a neomorphic protein.

To further explore the idea that exon skipping could produce a functional in-frame transcript, we asked whether CRISPR-mediated editing of Ctnnb1 exon 3 might induce exon skipping and cause a gain-of-function phenotype. Exon 3 of Ctnnb1 encodes phosphoacceptor residues that promote degradation of the β-Catenin transcription factor [19]; genetic excision of Ctnnb1 exon 3—which is in frame with exon 4—stabilizes a constitutively active β-Catenin that accumulates in the nucleus [20, 21]. We designed 11 sgRNAs that target regions along Ctnnb1 exon 3 (Ctnnb1-sg1 to -sg11), transduced individual sgRNAs into KP cells, and used high-throughput sequencing to analyze the extent of editing at the sgRNA target site in each line (Fig. 2b x-axis, Additional file 1: Figure S4 and Additional file 2: Table S4). Three sgRNAs (sg6, sg9, and sg10) inefficiently targeted Ctnnb1. Eight of the Ctnnb1 sgRNAs (sg1 to sg5, sg7, sg8, and sg11), however, induced indels at their target sites with frequencies that exceeded 20%. For example, Ctnnb1-sg1 generated + T insertions in about 65% of reads (Fig. 2c). In each population targeted by a strong Ctnnb1 sgRNA, we detected three RT-PCR products that span exons 2 to 5 (Fig. 2d). The major product corresponds to the normally spliced transcript that includes exon 3. The other two products correspond to alternatively spliced transcripts: one that skips exon 3 (i.e. exon 2-4 splicing, Fig. 2e) and one that skips both exons 3 and 4 (i.e., exon 2-5 splicing, Fig. 2f). Ctnnb1 sgRNAs targeting either DNA strand induced exon skipping and Cas9 nuclease activity was essential for exon skipping (Fig. 3a).

Fig. 2
figure2

Ctnnb1 sgRNAs targeting exon 3 induces exon skipping. a Schematic of the Ctnnb1 gene. The in-frame exon 3 encodes an inhibitory domain: phosphorylation amino acids 33, 37, 41, and 45 promotes degradation of the β-Catenin protein. Loss of exon 3 stabilizes β-Catenin. Eleven sgRNAs were designed to target exon 3: strong sgRNAs in red and weak sgRNAs in black, respectively. sgRNAs that use “NGG” PAM are shown above exon 3 and those that use “CCN” PAM are shown below exon 3. b Correlation between exon 3 skipping and sgRNA efficiency. Genomic indels were measured by deep sequencing. KP cells were infected with lentivirus. Exon 3 skipping efficiencies are from (d). Indels of sg11 were not determined. sgRNAs that induce > 20% indels are marked in red. c Distribution of sg1 indels shows that a T insertion (+T) at the Cas9 cleavage site nucleotide 97 of exon 3 (red arrowhead) was the most frequent. PAM sequence is in blue. d RT-PCR using primers spanning exons 2 and 5 shows partial exon skipping. M molecular marker. sgGFP is a control sgRNA. Exon 3 skipping bands were quantified using ImageQuant TL software and normalized to full length cDNA bands. sg4 showed visible weak bands that could not be quantified. e, f TOPO cloning and Sanger sequencing confirmed that the two major lower RT-PCR bands in (c) are alternative splicing of exon 2-4 and exon 2-5, respectively. g Western blot analysis of β-Catenin. Full length β-Catenin is ~86 kD. β-Catenin without exon 3 (delta exon 3) is ~77 kDa. Actin served as a loading control

Full size image
Fig. 3
figure3

Cas9 nuclease activity required for skipping of one or more exons. a RT-PCR analysis of Ctnnb1 mRNA in KP cells transduced with lentiviruses that encode sgCtnnb1.2 and nuclease-defective Cas9 (dCas9), dCas9-KRAB fusion, or WT Cas9. RT-PCR was performed using primers in exons 2 and 7 on transduced KP cell populations after puromycin selection and FACS sorting. The exon length and reading frame phase are shown. Only the exon 2-4 splice product retains an in-frame β-Catenin coding sequence. b RT-PCR analysis of Ctnnb1 mRNA in KP cells transduced with lentiviruses that encode Cas9 and sgGFP, sg3, or sg5. “–”, untreated

Full size image

Western blot analysis revealed that cell populations transduced with the strong sgRNAs produce a smaller ~74 kD β-Catenin protein that corresponds in size to that expected from the exon 2-4 splice product (Fig. 2g). The full length β-Catenin protein was not significantly depleted four days after transduction. To test whether the alternative splicing is dependent on the continuous expression of Cas9 or sgRNA in the lentiviral vectors, we co-transfected Cas9 and Ctnnb1-sg1 or a non-targeting sgRNA control. Seven days after transfection, when transfected Cas9 and guide RNAs should be depleted, we examined β-Catenin localization by immunofluorescence. In mouse fibroblast cells transfected with a non-targeting control sgRNA, β-Catenin localized to cell junctions (Additional file 1: Figure S5a). By contrast, in many cells transfected with Ctnnb1-sg1, we detected β-Catenin in the nucleus (Additional file 1: Figure S5a). These results suggest that continuous editing is not required for exon skipping and that exon 3 skipping induced by CRISPR-mediated editing of Ctnnb1 exon 3 produces a gain-of-function β-Catenin isoform.

We further analyzed transcripts spanning exons 2 to 7 in cell populations treated with Ctnnb1-sg2, -sg3, and -sg5. In addition to the full-length isoform, we detected four transcripts with exon 2 apparently spliced to each downstream exon (i.e. exon 2-4, exon 2-5, exon 2-6, and exon 2-7; Fig. 3a, b). We do not understand the mechanism of this apparently promiscuous exon skipping induced by Ctnnb1 exon 3 editing, nor have we been able to correlate promiscuous exon skipping with specific target sites or indel mutations in exon 3. Nevertheless, we isolated a Ctnnb1-sg3 edited clone that suggests a potential mechanism (Additional file 1: Figure S6a). This biallelic clone contains a 3-bp in-frame deletion on one allele and a large 832-bp deletion on the other; the 832-bp deletion fuses the 5’ end of intron 2 to the 3’ end of exon 4 (Additional file 1: Figure S6). We detected two transcripts in these cells: the properly spliced transcript that includes the 3-bp deletion and a transcript that includes intron 2 fused to exon 4 (Additional file 1: Figure S6c and Table 1). These results suggest that apparent exon skipping detected in populations of edited cells could reflect genome rearrangements that remove exons.

Two experiments support the idea that a single sgRNA can induce large genomic deletions that remove exons. For example, we isolated 15 clones from mouse 3T3 cells transiently transfected with Cas9 and Ctnnb1-sg1, and found that four clones (i.e. clones 4, 5, 13, and 15) showed apparent exon skipping by RT-PCR. Genomic PCR revealed genome rearrangements in three of these clones: large deletions (>500 bp) and smaller deletions (~100 bp) in clones 4 and 15, and large insertions in clones 13 and 15 (Additional file 1: Figure S7). Moreover, after targeting exon 6 of p65/RelA, we isolated a biallelic p65 clone (#15): one allele harbors a 1-nt “+A” insertion and the other harbors a 2268-bp deletion that removes exons 5, 6, and 7 (Additional file 1: Figure S8a, c–e). In p65 clone #15, we detected the fully spliced transcript and an exon 4-8 splice product (Additional file 1: Figure S8c). Both alleles encode frameshifted transcripts and both p65 transcripts are present at lower levels than WT (Additional file 1: Figure S8b). We also isolated an edited p65 clone (#31) homozygous for the same + A insertion as in clone #15, but clone #31 does not produce alternatively spliced transcripts. Thus, the exon 4-8 spliced transcript in clone #15 results from the deletion of exons 5, 6, and 7. These large exon deletion events were unexpected and would be missed using typical PCR-based screening assays.

The ability to cause a gain-of-function activity by inducing exon skipping or exon excision suggested that CRISPR-meditated editing using a single sgRNA might be a useful way to partially rescue function to a disease gene that requires low-level rescue. CRISPR-mediated homologous DNA repair has been used to correct premature stop codon mutations in the Dmd gene in a mouse model of DMD [22] and several groups have used CRISPR to delete Dmd exons and partially restore Dmd expression [23,24,25,26]. We designed four sgRNA/Cas9 lentiviruses that target different sites in exon 23 of the Dmd gene (Fig. 4a, b) and transduced mouse C2C12 myoblasts, a cell line widely used as a model for Duchenne muscular dystrophy (DMD) [27]. In C2C12 cells transduced with Dmd sgRNAs, we detected an RT-PCR product that corresponds to the normal splice product containing exon 23. Sequencing these RT-PCR products revealed that only Dmd-sg2 efficiently edited Dmd exon 23, as evidenced by mixed sequence peaks beyond the sgRNA target site (Additional file 1: Figure S9). In cells transduced with Dmd-sg2, we also detected an RT-PCR product corresponding to exon 22 spliced to exon 24 (Fig. 4c, d). Thus targeting exon 23 with one sgRNA might be sufficient to induce partial exon skipping and produce an intact dystrophin open reading frame. DMD is a classic example of a disease in which a small amount of functional restoration can provide substantial clinical benefit [28].

Fig. 4
figure4

An sgRNA targeting exon 23 of Dmd can partially restore in-frame dystrophin mRNA. a Schematic of sgRNA targeting and skipping of mouse Dmd exon 23 and location of primers for RT-PCR analysis. Skipping of exon 23 will generate in-frame mRNA. b sgRNA target sites in Dmd exon 23. c RT-PCR analysis of C2C12 mouse myoblast cells transduced with lentiviruses that encode Cas9 and sgDmd1, 2, 3, or 4. The expected band sizes are 353 bp and 140 bp. M molecular marker. d Sequence analysis of the 140-bp cDNA band from sgDmd2-treated cells confirmed splicing of exon 22 to exon 24

Full size image

Discussion

Whereas gene inactivation is most often the goal of CRISPR-mediated editing, our findings identify exon skipping as an unintended consequence of genome editing. We also show that exon skipping can result from indels that cause alternative splicing or from larger deletions that remove exons. Novel splice isoforms could encode proteins that retain partial function and should be carefully considered when interpreting phenotypes that result from CRISPR-induced mutations.

The frequency with which CRISPR-induced indels cause exon skipping is difficult to predict. Nevertheless, exon skipping caused by point mutations—including nonsense, missense, and translationally silent mutations—is well documented [29,30,31,32] and our results complement a recent study, which showed that CRISPR-mediated editing of the human FLOT1 gene can cause exon skipping by alternative splicing [18]. Roles for nonsense-mediated decay or cis-acting regulatory elements have been proposed, but mechanisms remain elusive. DNA damage has also been shown to regulate exon skipping [33]. Our data do not resolve whether the DNA damage, the indel, or the premature stop codon induces exon skipping, but they are consistent with the model that some indel mutations disrupt cis-acting sequences that promote splicing [29]. Future studies are needed to determine how CRISPR-induced indels cause alternative splicing and identify rules for predicting when exon skipping will occur.

We detected an unexpectedly high frequency of large deletions induced by CRISPR using a single sgRNA. We and others previously showed that two sgRNAs can generate large genomic deletion or inversion [34, 35]. However, large deletions induced by a single sgRNA have not been systematically analyzed in the literature. We initially missed these large deletions with the short-range PCR assays typically used to genotype CRISPR clones. We therefore recommend that, whenever possible, long-range PCR be used to genotype CRISPR clones. In many cases, large deletions will disrupt gene function and accomplish the goal of a CRISPR-mediated genome editing experiment. But our findings warrant careful analysis of editing events, because the aberrant juxtaposition and splicing of exons could result in neomorphic alleles.

Although exon skipping is an unintended consequence of CRISPR-mediated editing, we have shown that exon skipping can produce mRNAs that encode gain-of-function or partially functional proteins. Thus, exon skipping induced by CRISPR-mediated editing might be harnessed as a way to restore partial function to disease genes, in much the way that exon skipping induced by antisense oligonucleotides is being explored as a therapeutic to treat genetic diseases that result from splicing mutations [36].

Methods

CRISPR vectors

sgRNAs (Additional file 1: Table S1) were cloned into the lentiV2 (Addgene 52961) or pX330 (Addgene 42230) vectors using standard protocols [37].

Cell culture and infection

Cell culture conditions were as described [34]. A total of 293 fs cells were used to package lentiviruses encoding individual sgRNA and Cas9. KP cells or C2C12 cells were infected with lentiV2 lentiviruses and selected with puromycin. For Fig. 3a, cells were transduced with sgCtnnb1.2 cloned into lentiGuide-Puro vector (Addgene 52963), lenti Cas9-Blast (Addgene 52962), dCas9-BFP (Addgene 46910), or dCas9-KRAB-BFP (Addgene 46911). Cells were selected with puromycin, blasticidin, or FACS sorted for BFP.

Isolation of single-cell clones

KP or NIH-3T3 cells were transduced with lentiviruses Cas9 and sgRNAs targeting Kras, Ctnnb1, or p65 and were selected with puromycin for four days. For each transduction, 500 puromycin-resistant cells were seeded into a 100-mm dish and cultured until cell colonies were observed under a microscope. Individual colonies were transferred to 12-well plates—one colony per well—and grown to confluence. Genomic DNA and total RNA was isolated and PCR or RT-PCR was used to identify clones with indels, deletions, or insertions and exon skipping. Genomic PCR products were cloned into a TOPO vector to sequence alleles with indels or deletions.

CRISPR-induced insertion/deletion detection

Genomic DNA from cells was harvested by quick extraction buffer (Epibio), sgRNA target sites were amplified by PCR, and the products were sequenced on an Illumina NextSeq 500 [34]. We mapped the reads to the reference sequence using BWA (version 0.7.5) and SAMtools (version 0.1.19). VarScan2 (version 2.3) was used to identify insertions and deletions with the “pileup2indel” mode and parameters “--min-var-freq,” “--min-avg-qual,” and “--p-value.”

RNA-seq and bioinformatics analysis

RNA-seq libraries were generated using Illumina TruSeq kit, as described [38]. Paired-end 75-nt sequencing was performed using NextSeq. Reads were trimmed and primer sequences were removed using Trimmomatic (v 0.30). Reads were aligned to the mm10 genome using STAR (version 2.3.0e) with default parameters and uniquely mapping reads were selected. Redundant read pairs were removed using Samtools (version 0.0.19). For each gene annotated in GENCODE M7, the number of reads per gene was calculated using HTSeq. Percent exon 2 inclusion (Percentage Splicing Index, PSI or Ψ) for Kras was calculated as: (exon 1-2 + exon 2-3)/(exon 1-2 + exon 2-3 + exon 1-3). For global alternative splicing analysis, alternatively spliced exons were called using MISO 0.5.3 with default settings [39] and filtered with stringent cutoffs (ΔΨ ≥ 0.2, total reads ≥ 10, and Bayes factor ≥ 10).

Western blot analysis

Protein lysates from cultured cells were prepared in RIPA buffer with proteinase and phosphatase inhibitors. Proteins were separated on 4–12% NuPage Bis-Tris gels (Life Technologies, NP0321), transferred to nitrocellulose membrane, and probed with 1:1000 anti-β-Catenin antibody (BD 610154) or 1:5000 anti-Actin antibody (CST 8457).

RT-PCR and TOPO cloning

RNA was purified using RNeasy Mini Kit (Qiagen). First strand cDNA was synthesized using Superscript (ABI) and target sequences were amplified using LA-Taq (Clontech) or Herculase II (Agilent). Primers were listed in Additional file 1: Table S2. Gel bands were quantified using the ImageQuant TL software. Exon skipping products were gel purified, re-amplified using the same PCR protocol to increase yield, and cloned into a TOPO vector. TOPO clones were submitted to Genewiz for sequencing. Representative results from two lentiviral infections are shown.

References

  1. 1.

    Doudna JA, Charpentier E. Genome editing. The new frontier of genome engineering with CRISPR-Cas9. Science. 2014;346:1258096.

  2. 2.

    Cong L, Ran FA, Cox D, Lin S, Barretto R, Habib N, et al. Multiplex genome engineering using CRISPR/Cas systems. Science. 2013;339:819–23.

  3. 3.

    Mali P, Yang L, Esvelt KM, Aach J, Guell M, DiCarlo JE, et al. RNA-guided human genome engineering via Cas9. Science. 2013;339:823–6.

  4. 4.

    Mali P, Esvelt KM, Church GM. Cas9 as a versatile tool for engineering biology. Nat Methods. 2013;10:957–63.

  5. 5.

    Sander JD, Joung JK. CRISPR-Cas systems for editing, regulating and targeting genomes. Nat Biotechnol. 2014;32:347–55.

  6. 6.

    Fu Y, Foden JA, Khayter C, Maeder ML, Reyon D, Joung JK, et al. High-frequency off-target mutagenesis induced by CRISPR-Cas nucleases in human cells. Nat Biotechnol. 2013;31:822–6.

  7. 7.

    Hsu PD, Scott DA, Weinstein JA, Ran FA, Konermann S, Agarwala V, et al. DNA targeting specificity of RNA-guided Cas9 nucleases. Nat Biotechnol. 2013;31:827–32.

  8. 8.

    Xue W, Chen S, Yin H, Tammela T, Papagiannakopoulos T, Joshi NS, et al. CRISPR-mediated direct mutation of cancer genes in the mouse liver. Nature. 2014;514:380–5.

  9. 9.

    Wang D, Mou H, Li S, Li Y, Hough S, Tran K, et al. Adenovirus-mediated somatic genome editing of Pten by CRISPR/Cas9 in mouse liver in spite of Cas9-specific immune responses. Hum Gene Ther. 2015;26:432–42.

  10. 10.

    Song CQ, Li Y, Mou H, Moore J, Park A, Pomyen Y, et al. Genome-wide CRISPR screen identifies regulators of mitogen-activated protein kinase as suppressors of liver tumors in mice. Gastroenterology. 2016;152:1161–1173.e1161.

  11. 11.

    Yin H, Xue W, Chen S, Bogorad RL, Benedetti E, Grompe M, et al. Genome editing with Cas9 in adult mice corrects a disease mutation and phenotype. Nat Biotechnol. 2014;32:551–3.

  12. 12.

    Yin H, Song CQ, Dorkin JR, Zhu LJ, Li Y, Wu Q, et al. Therapeutic genome editing by combined viral and non-viral delivery of CRISPR system components in vivo. Nat Biotechnol. 2016;34:328–33.

  13. 13.

    Mou H, Moore J, Malonia SK, Li Y, Ozata DM, Hough S, et al. Genetic disruption of oncogenic Kras sensitizes lung cancer cells to Fas receptor-mediated apoptosis. Proc Natl Acad Sci U S A. 2017;114:3648–53.

  14. 14.

    Tuveson DA, Shaw AT, Willis NA, Silver DP, Jackson EL, Chang S, et al. Endogenous oncogenic K-ras(G12D) stimulates proliferation and widespread neoplastic and developmental defects. Cancer Cell. 2004;5:375–87.

  15. 15.

    Xue W, Meylan E, Oliver TG, Feldser DM, Winslow MM, Bronson R, et al. Response and resistance to NF-kappaB inhibitors in mouse models of lung adenocarcinoma. Cancer Discov. 2011;1:236–47.

  16. 16.

    Popp MW, Maquat LE. Leveraging rules of nonsense-mediated mRNA decay for genome engineering and personalized medicine. Cell. 2016;165:1319–22.

  17. 17.

    Cartegni L, Wang J, Zhu Z, Zhang MQ, Krainer AR. ESEfinder: A web resource to identify exonic splicing enhancers. Nucleic Acids Res. 2003;31:3568–71.

  18. 18.

    Kapahnke M, Banning A, Tikkanen R. Random splicing of several exons caused by a single base change in the target exon of CRISPR/Cas9 mediated gene knockout. Cells. 2016;5:45.

  19. 19.

    Moon RT, Kohn AD, De Ferrari GV, Kaykas A. WNT and beta-catenin signalling: diseases and therapies. Nat Rev Genet. 2004;5:691–701.

  20. 20.

    Iwao K, Nakamori S, Kameyama M, Imaoka S, Kinoshita M, Fukui T, et al. Activation of the beta-catenin gene by interstitial deletions involving exon 3 in primary colorectal carcinomas without adenomatous polyposis coli mutations. Cancer Res. 1998;58:1021–6.

  21. 21.

    Harada N, Tamai Y, Ishikawa T, Sauer B, Takaku K, Oshima M, et al. Intestinal polyposis in mice with a dominant stable mutation of the beta-catenin gene. EMBO J. 1999;18:5931–42.

  22. 22.

    Long C, McAnally JR, Shelton JM, Mireault A, Bassel-Duby R, Olson EN. Prevention of muscular dystrophy in mice by CRISPR/Cas9-mediated editing of germline DNA. Science. 2014;345:1184–8.

  23. 23.

    Long C, Amoasii L, Mireault AA, McAnally JR, Li H, Sanchez-Ortiz E, et al. Postnatal genome editing partially restores dystrophin expression in a mouse model of muscular dystrophy. Science. 2016;351:400–3.

  24. 24.

    Nelson CE, Hakim CH, Ousterout DG, Thakore PI, Moreb EA, Castellanos Rivera RM, et al. In vivo genome editing improves muscle function in a mouse model of Duchenne muscular dystrophy. Science. 2016;351:403–7.

  25. 25.

    Tabebordbar M, Zhu K, Cheng JK, Chew WL, Widrick JJ, Yan WX, et al. In vivo gene editing in dystrophic mouse muscle and muscle stem cells. Science. 2016;351:407–11.

  26. 26.

    Bengtsson NE, Hall JK, Odom GL, Phelps MP, Andrus CR, Hawkins RD, et al. Muscle-specific CRISPR/Cas9 dystrophin gene editing ameliorates pathophysiology in a mouse model for Duchenne muscular dystrophy. Nat Commun. 2017;8:14454.

  27. 27.

    Yaffe D, Saxel O. Serial passaging and differentiation of myogenic cells isolated from dystrophic mouse muscle. Nature. 1977;270:725–7.

  28. 28.

    Robinson-Hamm JN, Gersbach CA. Gene therapies that restore dystrophin expression for the treatment of Duchenne muscular dystrophy. Hum Genet. 2016;135:1029–40.

  29. 29.

    Liu HX, Cartegni L, Zhang MQ, Krainer AR. A mechanism for exon skipping caused by nonsense or missense mutations in BRCA1 and other genes. Nat Genet. 2001;27:55–8.

  30. 30.

    Cartegni L, Chew SL, Krainer AR. Listening to silence and understanding nonsense: exonic mutations that affect splicing. Nat Rev Genet. 2002;3:285–98.

  31. 31.

    Wang Z, Rolish ME, Yeo G, Tung V, Mawson M, Burge CB. Systematic identification and analysis of exonic splicing silencers. Cell. 2004;119:831–45.

  32. 32.

    Singh G, Pratt G, Yeo GW, Moore MJ. The clothes make the mRNA: past and present trends in mRNP fashion. Annu Rev Biochem. 2015;84:325–54.

  33. 33.

    Kalantari R, Chiang CM, Corey DR. Regulation of mammalian transcription and splicing by nuclear RNAi. Nucleic Acids Res. 2016;44:524–37.

  34. 34.

    Li Y, Park A, Mou H, Anderson DG, Jacks T, Weng ZP, et al. A versatile reporter system for CRISPR-mediated chromosomal rearrangements. Genome Biol. 2015;16:111.

  35. 35.

    Canver MC, Bauer DE, Dass A, Yien YY, Chung J, Masuda T, et al. Characterization of genomic deletion efficiency mediated by clustered regularly interspaced palindromic repeats (CRISPR)/Cas9 nuclease system in mammalian cells. J Biol Chem. 2014;289:21312–24.

  36. 36.

    Burghes AH, McGovern VL. Antisense oligonucleotides and spinal muscular atrophy: skipping along. Genes Dev. 2010;24:1574–9.

  37. 37.

    Sanjana NE, Shalem O, Zhang F. Improved vectors and genome-wide libraries for CRISPR screening. Nat Methods. 2014;11:783–4.

  38. 38.

    Shao DD, Xue W, Krall EB, Bhutkar A, Piccioni F, Wang X, et al. KRAS and YAP1 converge to regulate EMT and tumor survival. Cell. 2014;158:171–84.

  39. 39.

    Katz Y, Wang ET, Airoldi EM, Burge CB. Analysis and design of RNA sequencing experiments for identifying isoform regulation. Nat Methods. 2010;7:1009–15.

Download references

Acknowledgements

We thank T. Jacks, C. Mello, P. Zamore, T. Flotte, V. Ambros, S. Wolfe, Z. Kennedy, T. Partridge, J. Chen, and T. Pederson for their insightful comments and reagents. We thank D. Conte for critically editing this manuscript. We thank E. Kittler in the UMass Medical School Deep Sequencing Core for support.

Funding

HY was supported by the Skoltech Center and the NCI in the MIT-Harvard Center of Cancer Nanotechnology Excellence (5-U54-CA151884-04). YL was supported by the China Scholarship Council (201506260151) and the Thousand Talent Plan funding to ZW from the Chinese government. This work was supported by grants from the NIH (DP2HL137167, P01HL131471, and R00CA169512 to WX; GM115911 to EJS), the American Cancer Society (129056-RSG-16-093), and the Lung Cancer Research Foundation to WX.

Availability of data and materials

Deep sequencing data are available under BioProject ID: PRJNA375870.

Authors’ contributions

HM, JS, MJM, ZW, and WX designed the study. HM, JS, LP, HY, JM, XZ, CQS, AS, QW, DMO, and YL performed experiments or analyzed data. HM, JS, DGA, CPE, EJS, MJM, ZW, and WX wrote the manuscript with comments from all authors. All authors read and approved the final manuscript.

Competing interests

The authors declare that they have no competing interests.

Ethics approval and consent to participate

Ethics approval was not needed for the study.

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Author information

Author notes

    Affiliations

    1. RNA Therapeutics Institute, University of Massachusetts Medical School, Worcester, MA, 01605, USA
      • Haiwei Mou
      • , Jordan L. Smith
      • , Lingtao Peng
      • , Chun-Qing Song
      • , Ankur Sheel
      • , Deniz M. Ozata
      • , Erik J. Sontheimer
      • , Melissa J. Moore
      •  & Wen Xue
    2. David H. Koch Institute for Integrative Cancer Research, Massachusetts Institute of Technology, Cambridge, MA, 02142, USA
      • Hao Yin
      • , Qiongqiong Wu
      • , Daniel G. Anderson
      •  & Zhiping Weng
    3. Program in Bioinformatics and Integrative Biology, University of Massachusetts Medical School, Worcester, MA, 01605, USA
      • Jill Moore
      • , Xiao-Ou Zhang
      • , Yingxiang Li
      •  & Zhiping Weng
    4. Department of Bioinformatics, School of Life Science and Technology, Tongji University, Shanghai, People’s Republic of China
      • Yingxiang Li
    5. Department of Chemical Engineering, Massachusetts Institute of Technology, Cambridge, MA, 02142, USA
      • Daniel G. Anderson
    6. Harvard-MIT Division of Health Sciences & Technology, Cambridge, MA, 02139, USA
      • Daniel G. Anderson
    7. Institute for Medical Engineering and Science, Massachusetts Institute of Technology, Cambridge, MA, USA
      • Daniel G. Anderson
    8. Wellstone Muscular Dystrophy Program, Department of Neurology, University of Massachusetts Medical School, Worcester, MA, 01605, USA
      • Charles P. Emerson
    9. Program in Molecular Medicine, University of Massachusetts Medical School, Worcester, MA, 01605, USA
      • Erik J. Sontheimer
      •  & Wen Xue
    10. Department of Biochemistry and Molecular Pharmacology, Howard Hughes Medical Institute, University of Massachusetts Medical School, Worcester, MA, 01605, USA
      • Melissa J. Moore
    11. Department of Molecular, Cell and Cancer Biology, University of Massachusetts Medical School, 368 Plantation Street, Worcester, MA, 01605, USA
      • Wen Xue
    Authors
    1. Search for Haiwei Mou in:
    2. Search for Jordan L. Smith in:
    3. Search for Lingtao Peng in:
    4. Search for Hao Yin in:
    5. Search for Jill Moore in:
    6. Search for Xiao-Ou Zhang in:
    7. Search for Chun-Qing Song in:
    8. Search for Ankur Sheel in:
    9. Search for Qiongqiong Wu in:
    10. Search for Deniz M. Ozata in:
    11. Search for Yingxiang Li in:
    12. Search for Daniel G. Anderson in:
    13. Search for Charles P. Emerson in:
    14. Search for Erik J. Sontheimer in:
    15. Search for Melissa J. Moore in:
    16. Search for Zhiping Weng in:
    17. Search for Wen Xue in:

    Corresponding authors

    Correspondence to Melissa J. Moore or Zhiping Weng or Wen Xue.

    Additional files

    Additional file 1:

    Supplemental Figures and Tables 1–3. (PDF 4407 kb)

    Additional file 2:

    Supplemental Table 4. (XLSX 32 kb)

    Rights and permissions

    Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

    Reprints and Permissions

    About this article

    Verify currency and authenticity via CrossMark

    Cite this article

    Mou, H., Smith, J.L., Peng, L. et al. CRISPR/Cas9-mediated genome editing induces exon skipping by alternative splicing or exon deletion. Genome Biol 18, 108 (2017) doi:10.1186/s13059-017-1237-8

    Download citation

    Genome Biology

    ISSN: 1474-760X

    Contact us