Abstract
The CD19 antigen, expressed on most B-cell acute lymphoblastic leukemias (B-ALL), can be targeted with chimeric antigen receptor–armed T cells (CART-19), but relapses with epitope loss occur in 10% to 20% of pediatric responders. We detected hemizygous deletions spanning the CD19 locus and de novo frameshift and missense mutations in exon 2 of CD19 in some relapse samples. However, we also discovered alternatively spliced CD19 mRNA species, including one lacking exon 2. Pull-down/siRNA experiments identified SRSF3 as a splicing factor involved in exon 2 retention, and its levels were lower in relapsed B-ALL. Using genome editing, we demonstrated that exon 2 skipping bypasses exon 2 mutations in B-ALL cells and allows expression of the N-terminally truncated CD19 variant, which fails to trigger killing by CART-19 but partly rescues defects associated with CD19 loss. Thus, this mechanism of resistance is based on a combination of deleterious mutations and ensuing selection for alternatively spliced RNA isoforms.
Significance: CART-19 yield 70% response rates in patients with B-ALL, but also produce escape variants. We discovered that the underlying mechanism is the selection for preexisting alternatively spliced CD19 isoforms with the compromised CART-19 epitope. This mechanism suggests a possibility of targeting alternative CD19 ectodomains, which could improve survival of patients with B-cell neoplasms. Cancer Discov; 5(12); 1282–95. ©2015 AACR.
See related commentary by Jackson and Brentjens, p. 1238.
This article is highlighted in the In This Issue feature, p. 1225
Introduction
Despite significant advances in the treatment of pediatric B-cell acute lymphoblastic leukemias (B-ALL), children with relapsed or refractory disease still account for a substantial number of all childhood cancer deaths. Adults with B-ALL experience even higher relapse rates and long-term event-free survival of less than 50% (1). Relapsed leukemia is generally not curable with chemotherapy alone, so the prospect of long-term disease control via an immunologic mechanism holds tremendous promise. One of the most innovative approaches involves the use of adoptive T cells expressing chimeric antigen receptors (CAR-T) against CD19 (2, 3). Despite obvious successes, there have been documented relapses in which CART-19 cells were still present, but the leukemia cells lost surface expression of CD19 epitopes, as detected by clinical flow cytometry. According to the recent estimates, epitope loss occurs in 10% to 20% of pediatric B-ALL treated with CD19-directed immunotherapy (4, 5), raising questions about its significance for neoplastic growth.
The cell surface signaling protein CD19 is required for several diverse processes in B-cell development and function. In the bone marrow, CD19 augments pre–B-cell receptor (pre-BCR) signaling (6, 7), thereby promoting the proliferation and differentiation of late pro-B cells bearing functional immunoglobulin heavy chains into pre-B cells. Engaging the CD19 pathway in normal and neoplastic B-lineage cells induces the activation of the growth-promoting kinases PI3K and LYN, which are activated via intracellular interactions with conserved tyrosine residues in the CD19 cytoplasmic tail (8). Significantly, whereas CD19 possesses conserved extracellular domains needed for mature B-cell function (9), the role of CD19 ectodomains in the proliferation and differentiation of normal B-lineage precursors is unknown. Likewise, CD19 is thought to play an essential role in B-cell neoplasm, but it is usually attributed to its ability to recruit intracellular kinases (10–12).
Results
Post–CART-19 Pediatric B-ALL Relapses Retain and Transcribe the CD19 Gene
To study mechanisms and consequences of CD19 loss in vivo, we reanalyzed CD19-positive pre–CART-19 leukemia and relapsed CD19-negative leukemia obtained from the same patient (ref. 13; CHOP101/101R in Fig. 1A, top). We also studied two sequential relapses after CART-19 therapy from patient CHOP105. The first CD19-positive relapse (R1) was due to the loss of CAR-T cells, and the patient achieved complete remission following CART-19 reinfusion. However, the second relapse (R2) was accompanied by loss of the CD19 epitope (Fig. 1A, bottom) and rapid disease progression. Upon successful engraftment in NOD scid gamma (NSG) mice, these four paired leukemia samples were used for molecular analyses. Samples CHOP101/101R were subjected to whole-genome sequencing, and we observed no copy-number variations or focal deletions in the CD19 locus (Fig. 1B). Clinical karyotyping and LOH analysis of samples CHOP105R1/R2 revealed a very large hemizygous deletion within chromosome 16 extending from p13.11 to p11.1 (Fig. 1C) and spanning the entire CD19 locus.
To further characterize the B-ALL samples, we performed whole-exome sequencing (WES) and RNA sequencing as well as copy-number alteration (CNA) analysis. These approaches revealed the existence in relapsed leukemias of de novo genomic alterations primarily, but not exclusively, affecting exon 2. In sample CHOP101R, we observed two independent frameshift mutations (one in exon 2 and one in exon 4); however, they were each subclonal and accounted for less than 50% of tumor cells. In the CHOP105 samples, we identified the insertion of 3 codons in exon 2, which was detectable with very low frequency by RNA sequencing (RNA-seq) in the R1 leukemia but became clonal in the R2 leukemia (Table 1). To better understand the relevance of such mutations, we analyzed three other post–CART-19 relapses: CHOP107Ra/107Rb and CHOP133R, for which matched baseline samples were not available. Neither of the CHOP107R samples (which had been xenografted from the same patient at different times during disease progression) contained mutations. However, leukemia CHOP133R suffered hemizygous loss of the entire chromosome 16, and the remaining allele contained a frameshift mutation also in exon 2 (Table 1), which could have led to nonsense-mediated decay (14). In summary, genetic alterations could have accounted for CD19 protein loss in some (e.g., CHOP105R2 and CHOP133R) but not in other (e.g., CHOP101R and CHOP107a/b) samples, prompting us to investigate transcriptional deregulation.
Using bisulfite-based sequencing, we first showed that there was no increase in methylation of CD19 promoter or enhancer elements, which could have accounted for gene silencing in the two matched relapse samples (Fig. 1D). We also performed qRT-PCR for PAX5, the key regulator of CD19 transcription (15), but observed no consistent downregulation of PAX5 mRNA (Fig. 1E). More surprisingly, CD19 mRNA levels were found to be downregulated only 2- to 3-fold, depending on the choice of primers (Fig. 1F). The discrepancy between mRNA and protein levels suggested that post–CART-19 samples may have altered regulation of transcript processing.
Alternatively Spliced CD19 mRNA Variants Accumulate in Post–CART-19 Relapses
The SIB Genes Track (16) implemented in the UCSC Genome Browser postulates the existence of CD19 mRNA isoforms skipping exons 2 and 5 to 6 (Fig. 1G). To study these isoforms, we confirmed sustained CD19 mRNA expresssion in relapsed tumors using RNA-seq (Fig. 2A) and then aligned CHOP101/101R RNA-seq reads to CD19 exons using the MAJIQ algorithm (ref. 17; Fig. 2B, top). These alignments were used to estimate the relative inclusion [percent spliced in (PSI)] of splicing variants and visualize them in violin plots generated by VOILA (ref. 17; Fig. 2B, bottom). This analysis revealed that in CHOP101 exon 4 is always spliced to exon 5, whereas in CHOP101R 25% to 30% of the observed transcripts skip exon 5 to 6, leading to juxtaposition of exons 4 and 7. We also observed a trend toward fewer reads spanning the ex1/2 and ex2/3 junctions in CHOP101R (Fig. 2B, bottom).
To further validate these changes, we performed additional analyses on CHOP101/101R and CHOP105R1/105R2. The appearance of the Δex5–6 splicing variant in the relapsed samples was detected using very stringent radioactive low-cycle semiquantitative RT-PCR (Fig. 2C) and confirmed by Sanger sequencing of RT-PCR products (Supplementary Fig. S1A). When exon 1–4-specific primers were used, in both samples with CD19 epitope loss, there was 2.5- to 4.5-fold increased abundance of Δex2 and decreased levels of the full-length isoform (Fig. 2D, left and Supplementary Fig. S1B). The Δex2 isoform was also detectable in two other post–CART-19 leukemias for which no matching pretreatment samples were available: CHOP107Ra and CHOP133R (Fig. 2D, right and Supplementary Fig. S1C).
To perform even more stringent quantification, we designed a forward primer spanning the exon1–exon3 junction and thus specific for the Δex2 isoform. By qRT-PCR, we confirmed a sharp increase in exon 2 skipping in CHOP101R leukemia relative to CHOP101 (Fig. 2E). To determine if some Δex2 mRNA species retain exons 5 to 6, we designed another pair of primers to amplify the exon1–exon5 fragment. Using CHOP101R cDNA as template, we observed fragments corresponding to both full-length and Δex2 isoforms (red/green in Fig. 2F). Sanger sequencing of these bands confirmed their makeup and revealed a frameshift mutation present in the full-length (but not Δex2) isoform (Fig. 2G).
Consistent skipping of exon 2 prompted us to reevaluate the seemingly deleterious frameshift mutations in exon 2 found in CHOP101R and CHOP133R (Supplementary Fig. S1D). We used the CRISPR/Cas9 system with a guide RNA homologous to exon 2 to introduce double-stranded breaks in this exon in various B-cell lines and allowed them to repair by nonhomologous end joining. Frameshift events were selected for using sorting for CD19-negative cells and confirmed by sequencing. In all three cell lines tested [697, NALM-6 (both B-ALL), and Raji (Burkitt lymphoma)], frameshifts resulted in expression of a large CD19 protein isoform consistent in size with exon 2 skipping (Supplementary Fig. S1E)—despite alternative stop codons downstream of the mutations site. Exon skipping was confirmed at the mRNA levels by qRT-PCR with exon1–exon 3 junction-specific primers (Fig. 2H and Supplementary S1F). Thus, alternative splicing of exon 2 can override normally deleterious mutations.
The SRSF3 Splicing Factor Binds to and Regulates Inclusion of CD19 Exon2
To understand the mechanism behind CD19 mRNA splicing, we used the AVISPA algorithm, which predicts RNA-binding proteins specific to particular intron–exon cassettes (18). The predictions for exons 2 and 5 to 6 had a considerable overlap, consisting of NOVA, HuD, hnRNP-C, hnRNP-F/H, hnRNP-G, PTBP1/2, SRSF2, SRSF3, and CELF1/2 (Supplementary Figs. S2A and S3A). We then generated T3-transcribed 32P-labeled CD19 RNA containing intron 1/exon 2/intron 2, cross-linked it to protein lysates from B-ALL cells, and separated protein thus labeled by PAGE. The sizes of the observed bands were consistent with molecular weights of AVISPA-predicted as well as additional splice factors (SF), such as hnRNP-M, hnRNP-A1, hnRNP-U, SRSF7, and PSF (Supplementary Fig. S2B; ref. 19). Nevertheless, most of these factors were negative by immunoprecipitation with SF-specific antibodies (Supplementary Fig. S2C). In contrast, SRSF3 and hnRNP-A and –C were all positively confirmed in the same assay (Fig. 3B). To this list we added SRSF2, which is thought to act in concert with SRSF3 (20).
To determine if any of these SFs were involved in exon 2 alternative splicing, we tested four siRNA pools in B-lymphoid P493–6 cells amenable to efficient transfection (21) and in all cases observed efficient knockdown at the mRNA level (Supplementary Fig. S2D). However, only SRSF3 knockdown affected skipping of exon 2, as evidenced by the qRT-PCR assay (Fig. 3C). The knockdown experiment was repeated in NALM-6 B-ALL cells where a 65% to 75% decrease in SRSF2 and SRSF3 mRNA levels was achieved by siRNA transfection (Supplementary Fig. S2E, left). Once again, only SRSF3 but not SRSF2 knockdown affected exon 2 processing (Supplementary Fig. S2E, right). Most importantly, knockdown of SRSF3 resulted in increased abundance of the Δex2 protein isoform in both P493-6 and NALM-6 B-ALL cells, as measured by immunoblotting for CD19 (Fig. 3D). To further confirm the role of SRSF3 in CD19 exon 2 retention, we mined the publicly available GSE52834 dataset where 22 RNA-binding proteins were knocked down in the GM19238 lymphoblastoid cell line. Of note, only knockdown of SRSF3 resulted in a significant increase in CD19 exon 2 skipping (Fig. 3E). We then asked whether any SRSF3 sites are present in exon 2 of CD19. Unfortunately, the commonly used ESE-Finder tool (22) does not include binding motifs for human SRSF3, because the consensus is not well defined. However, the Drosophila homolog of SRSF3, RBP-1, is known to bind to the [A/T]CAAC[A/G] hexamer (23). Of note, this motif is found twice in CD19 exon 2, where it is not directly affected by de novo CD19 mutations (Supplementary Fig. S2D).
To determine how SRSF3 function could be impaired in relapse leukemias, we assessed SRSF3 expression levels in CHOP101/101R and CHOP105R1/105R2 matched sets. In both cases, relapsed leukemias expressed lower amounts of SRSF3. Also, two other post–CART-19 relapses, CHOP107R and CHOP133R (for which matched baseline samples were not available), expressed even lower levels of this protein (Fig. 3F, left). In parallel, we measured the protein levels of hnRNPC1/C2 and hnRNPA1, but there was no consistent pattern of change for either of these splicing factors in paired post– versus pre–CART-19 samples (Fig. 3F, right). Taken together, these results suggest that SRSF3 insufficiency in relapsed leukemias could be at least partly responsible for the abundance of the CD19 Δex2 isoform.
The CD19 Δex2 Isoform Partially Rescues Defects Associated with CD19 Loss
The detected alterations in exon inclusion should result in truncated CD19 variants, with profound implications for both CD19 functionality and CART-19 recognition. Skipping of exon 2 could compromise the FMC63 epitope targeted by the CAR (24, 25), making it invisible to this immunotherapy. Skipping of exons 5 and 6 would result in premature termination and elimination of the transmembrane and the cytosolic domains (Fig. 4A). The expected truncated variants were readily detectable in leukemia cell lysates by immunoblotting using antibodies recognizing either extracellular or cytoplasmic epitopes (Fig. 4B), the hallmark of relapsed leukemias being the lack of the full-length isoform. CD19 Δex2 was also detectable in all human B-cell lines tested (Fig. 4C), and the corresponding mRNA was as stable as the full-length transcript (Fig. 4D), attesting to their possible significance.
To test this hypothesis, we generated a series of CD19-encoding retroviruses (Supplementary Fig. S1G, left) and transduced them into the previously generated murine B-cell line Myc5, which had lost endogenous CD19 expression following silencing of the CD19 transcriptional regulator PAX5 (26, 27). In this system, retrovirally encoded Δex2 and Δex5–6 isoforms were robustly expressed (Fig. 4E). Interestingly, when we measured half-lives of CD19 protein isoforms using treatment with cycloheximide, we observed an increase in Δex2 protein stability compared with the full-length isoform (Fig. 4F). As predicted, the Δex2 isoform was not recognized by the CD19 flow cytometry antibody (Fig. 4G). Importantly, in Myc5 cells restoration of full-length CD19 resulted in enhanced proliferation [consistent with our prior data (10)], and the Δex2 isoform (but not Δex5–6) partly recapitulated this growth phenotype (Fig. 4H).
To establish relevance of this finding to human disease, we generated NALM-6 and 697 B-ALL subclones, in which the endogenous CD19 gene was knocked out using the CRISPR/Cas9 system, resulting in the loss of CD19 expression (Fig. 5A). We then reconstituted it with either full-length or Δex2 CD19 and confirmed robust protein expression by immunoblotting (Fig. 5B and C). We also transduced them with CD19–GFP fusion–encoding retroviruses (Supplementary Fig. S1G, right). Unexpectedly, confocal microscopy revealed that unlike full-length CD19-GFP, which localizes exclusively to the plasma membrane, the CD19 Δex2–GFP isoform is largely cytosolic; however, up to 10% of this isoform can still be found on the membrane (Figs. 5D). Further experiments were performed to validate the relevance of this fraction.
Glycosylation of CD19 is a prerequisite for plasma membrane localization (28), and, unlike its intracellular precursor, plasma membrane–bound CD19 is susceptible to extracellular cleavage by trypsin (29). To determine whether the Δex2 isoform is glycosylated, whole-cell protein lysates obtained from CD19 retrovirus-transduced cultures were treated with a mix of glycosylases followed by Western blotting. Just like its full-length counterpart, the Δex2 isoform was reduced in size upon treatment (Fig. 5E), indicating that it is glycosylated and that some of it could be transported to the plasma membrane. To quantitate the membrane-bound fraction, we incubated reconstituted live cells with trypsin. As expected, almost all of the full-length CD19 were cleaved by trypsin, whereas most of the Δex2 isoform and all of the Δex5–6 isoform retained their original size. However, over 10% of the CD19 Δex2 protein was sensitive to trypsinization, fully consistent with the results of confocal microscopy (Fig. 5F).
To test whether this plasma membrane–associated fraction is functional, we asked whether it contributes to tonic or antigen-driven pre-BCR signaling by directly recruiting PI3K and Src family tyrosine kinases (SFTK), such as LYN (21, 30, 31). In both full-length and Δex2-reconstituted cells, PI3K and LYN coimmunoprecipitated with CD19, albeit less abundantly in the latter case, reflecting a much smaller pool of plasma membrane–associated Δex2. When pre-BCR was ligated with the anti-IgM antibody, there was an increase in CD19 Δex2-LYN binding, although the amount of CD19 Δex2–bound PI3K was reduced (Fig. 5G). Moreover, similar to reconstituted murine Myc5 cells, human Δex2 cells grew in culture almost as rapidly as their full-length CD19 counterparts and significantly faster than control CD19 Δex5–6 cells (Fig. 5H). In principle, the presence of functional Δex2 on the plasma membrane could be sufficient to trigger killing by CART-19 cells. However, when exposed to CART-19, only the full-length CD19 cultures were killed, whereas CD19 Δex2–transduced cells remained fully viable (Fig. 5I), confirming the loss of the cognate CART-19 epitope.
Discussion
Our study addresses the important clinical issue of resistance to CART-19 and establishes a novel combinatorial mechanism by which its cognate epitope could be removed from the cell surface without discarding the target protein entirely. This mechanism involves the clustering of nonsense and missense mutations in exon 2 of CD19. Distributed frameshift mutations would have prevented CD19 protein expression but also left the leukemic cells without the important activator of PI3K and SFTK signaling. In contrast, frameshift mutations clustered in the nonconstitutive exon 2 eliminate full-length CD19, but allow expression of the Δex2 isoform. Not only is this isoform mostly cytosolic and thus hidden from T cells, but expression of its membrane fraction does not trigger killing by CART-19, at least not at physiologic levels. At the same time, it was found to be even more stable than full-length CD19, which could be due to either the presence of a degron within exon 2–encoded amino acid sequence or mislocalization of CD19 Δex2 protein away from its normal degradation pathways. Moreover, this isoform at least partly rescues defects in cell proliferation and pre-BCR signaling associated with CD19 loss. Thus, its retention in relapsed B-ALL is highly advantageous for leukemic cells, whether or not they carry de novo mutations in exon 2. At present, we cannot distinguish between such mutations merely selecting for Δex2 (the “permissive” model) versus actively redirecting the splicing machinery (the “instructive” model). Although the de novo CD19 mutations we have discovered do not interrupt the putative SRSF3 site, they could still compromise its binding, a scenario we plan to address in future studies.
A related question is whether epigenetics plays any role in the recruitment of SRSF3 and other splicing factors that likely cooperate with it. It is well known that splicing occurs co-transcriptionally when pre-mRNA is still in the vicinity of chromatin, which can influence intron removal (32–34). Certain histone modifications are enriched on chromatin associated with exonic sequences (35), and spliceosome machinery is recruited via cofactors recognizing histone modifications and/or associates directly with modified histones. For example, H3K36me3 interacts with PSIP1, which then recruits various SFs, including SRSF3 (36). The underlying mechanism notwithstanding, it is becoming clear that alterations in splicing factors are important drivers of hematologic malignancies, as evidenced by the discovery of acquired mutations in the splicing factor gene SF3B1 in chronic lymphocytic leukemia (37) and myelodysplasia (38).
Whether hematologic malignancies are driven by global deregulation of splicing or by alterations in select target genes is not known. Although our data do not resolve this issue, they underscore the importance of splicing alterations at the level of individual genes, at least in the context of disease progression. Similarly, in the realm of solid tumors, BRAFV600E splicing variants lacking the RAS-binding domain were found in one third of melanomas with acquired resistance to vemurafenib (39). The existence of such splicing-based adaptive mechanisms suggests that future CARs and other antibody-based therapeutics should be designed to target essential exons, as a way to prevent escape (40).
On the other hand, it is conceivable that splicing is globally deregulated in B-ALL, either owing to downregulation of SRSF3 and related splicing factors or due to pervasive epigenetic changes. In that case, it should be possible to define a set of genes that are alternatively spliced in B-ALL versus normal B cells and encode extracellular epitopes. Such epitopes could be targets for completely new chimeric antigen receptors capable of killing B-ALL blasts while sparing normal B cells, with the selectivity CART-19 does not possess.
Methods
Cell Culture, Transfections, Treatments, and Infections
All B-lymphoid cell lines (NALM-6, MYC5, 697, and P493-6) were cultured and maintained in RPMI-1640 medium supplemented with 10% FBS, 2 mmol/L l-glutamine, and penicillin/streptomycin at 37°C and 5% CO2. SMARTpool siRNAs for splicing factors SRSF3, SRSF7, hnRNPC, and hnRNPA (Dharmacon) were transfected at indicated concentrations into B-cell lines by electroporation using the AMAXA system program 0-006 and Reagent V (Lonza). siRNA knockdown efficiency was measured 24 and 48 hours after transfection by RT-qPCR. BCR ligation was performed by incubation of 20 × 106 cells with 10 μg/mL of pre–BCR-specific α-IgM Jackson Immuno antibody (IgM-5μ) or with isotype control goat anti-IgG (Southern Biotech; #0109-01) for indicated time points at room temperature. Cells were lysed in RIPA buffer and loaded onto PAGE gels for immunoblotting analysis. Cleavage of plasma membrane proteins by trypsin was performed by incubation of 1 × 106 cells in 200 μL of 1× trypsin-EDTA solution (Gibco; #15400-054) in PBS for 4 minutes at 37°C. Control cells were incubated under the same conditions in PBS. Trypsinization was stopped by adding 1 ml of ice-cold PBS/10% FBS followed by quick centrifugation and immediate cellular lysis for whole-cell protein extraction. Protein half-life was measured by treating cells with cycloheximide (Sigma) at 50 μg/mL. mRNA half-life was measured by treating the cells with actinomycin D (Sigma) 5 μg/mL.
Cell Line Authentication
NALM-6 and 697 cells were obtained from the Center for Childhood Cancer Research biodepository in 2009. NALM-6 cells were authenticated using fingerprinting based on multiplex PCR of minisatellite markers. P493-6 cells were obtained from Dr. Dirk Eick in 2003 and authenticated in 2010 through targeted resequencing of the transgenic MYC allele. Other cell lines have not been authenticated.
Retroviral and Lentiviral Constructs
Lentiviral vector expressing luciferase and GFP pELNS-CBR-T2A-GFP has been previously described (41). Retroviral constructs expressing full-length CD19 cDNA were generated by digestion of pMX-IRES-CD19-GFP vector (10) with EcoRI/XhoI restriction enzymes, followed by ligation into MSCV-IRES-DsRedFP (Addgene) and pMXs-IRES-Blasticidin (RTV-016; Cell Biolabs) retroviral backbones. To generate CD19 Δex2– and CD19 Δex5–6-expressing vectors, the cDNA fragments (Supplementary Table S1) were synthesized (Genewiz) and cloned into MSCV-CD19-IRES-DsRedFP via EcoRI/BglII or BglII/XhoI, and later moved into pMX-IRES-Blast via EcoRI/XhoI cloning. Retroviral and lentiviral particles were generated by transfection of GP293 cells with Lipofectamine-2000 (Invitrogen). Viral supernatants were harvested 24, 36, and 48 hours after transfection and used to infect B-ALL cell lines in the presence of polybrene (4 μg/mL). Where indicated, selection of infected cells was done with 10 μg/mL Blasticidine over the course of 1 week, or by cell sorting.
CRISPR/Cas9 Genome Editing System
CD19-CRISPR/Cas9-KO plasmid was obtained from Santa Cruz Biotechnologies (sc-400719) and transfected into NALM-6 and 697 cell lines via electroporation using the AMAXA system program 0-006 and Reagent V (Lonza). Cells were stained with α-CD19-phycoerythrin (PE)–conjugated antibody (Beckman Coulter) 4 days after transfection, and CD19-deficient (ΔCD19) cells were sorted and plated individually in 96-well clusters for single-cell clone selection and expansion, or maintained as a pool. CD19 knockdown was confirmed by flow cytometry and by Western blot using antibodies that recognize epitopes in the cytosolic and the extracellular domains of CD19. DNA and RNA were extracted, and the CD19 gene was sequenced to analyze the mutations induced by the CRISPR/Cas9 system. To generate frameshift mutations into CD19 exon 2, a single CRISPR/Cas9 exon 2–gRNA plasmid was transfected by electroporation into 697, NALM-6, and Raji cell lines as described above. Effective insertion of frameshift mutations at expected targeted region was assessed by Sanger sequencing.
Immunofluorescence and Colocalization Studies
The 697 ΔCD19 cells expressing CD19-GFP and CD19 Δex2–GFP fusion proteins were incubated and stained with 5 μg/mL Wheat Germ Agglutinin Alexa Fluor-680 (Molecular Probes; #W32465) according to the manufacturer's instructions. Once fixed, cells were mounted on precharged glass microscope slides with DAPI-containing medium (Vectashield; cat# H1200) and visualized under a Leica STED 3× Super-Resolution Confocal System HC PL APO CS2 63x/1.40 Oil 63× objective. Images were acquired using 4,184 × 4,184 resolution with limited signal saturation. Colocalization was quantified by Pearson correlation coefficient. Six images for each CD19 construct containing 100 cells on average were analyzed with BioImageXD and FIJI Coloc2 plugin software. Statistical Costes P value is 1 for this analysis (42).
Cell Proliferation Assays
Myc5 cells expressing CD19-FL, CD19 Δex2, or empty Blasticidine vector were seeded in 10 mL of medium at 100,000 cells/mL. Daily samples were taken and counted by flow cytometry, and cell density was calculated based on absolute counts. Each cell line was assayed in triplicate, and each assay was repeated two times. The 697 ΔCD19 cells expressing CD19-FL, CD19-Δex2, or CD19-Δex5–6 vector together with pELNS-CBR-T2A-GFP were seeded in triplicate in a standard 96-well plate, 10,000 cells per well in 100 μL of media. GFP fluorescent signal was measured at 485 nm excitation and 528 nm emission in a Synergy2 (Biotek) plate reader daily for 4 days. Each cell line was assayed in triplicate, and each assay was repeated two times. Proliferation rates were statistically compared by the Student t test at indicated times points (*, P ≤ 0.05; **, P < 0.01).
Xenografted Tumor Samples
Xenograft models of tumor samples have been previously described (41). All animal experiments were approved by the Institutional Animal Care and Use Committee of the Children's Hospital of Philadelphia.
Cytotoxicity Assays
NALM-6 ΔCD19 cells expressing CD19-FL, CD19 Δex2, CD19 Δex5–6, or empty vector together with pELNS-CBR-T2A-GFP were used as targets for T-cell cytotoxicity assay as previously described (43). Briefly, target cells (T) were incubated with effector (E) T cells (CART-19) at the indicated E:T ratios for 24 hours. D-luciferin (Goldbio; #LUCK-1G) was then added to the cell culture, and bioluminescence imaging was performed on a Xenogen IVIS-200 Spectrum camera. Target killing was analyzed using the software Living Image 4.3.1 (Caliper LifeSciences).
Flow Cytometry
Live cells were stained with PE-conjugated CD19 antibody (IM1285U; Beckman Coulter) and analyzed in an AccuriC6 cytometer as previously described (10, 13).
DNA Extraction and Sequencing of CD19 Gene
Genomic DNA was obtained from 2 × 106 cells from xenografted pre-CART and post-CART tumor samples using a DNeasy Blood & Tissue Kit (Qiagen). The CD19 gene, expanding 1.2 kb upstream to include the enhancer and promoter regions, was amplified by PCR and sequenced. Primer sets used are described in Supplementary Table S2.
Methylation of Promoter and Enhancer Region
Genomic DNA from xenografted tumor samples was subjected to bisulfate conversion using the Epitect Fast DNA Bisulfite Kit (Qiagen). The CD19 enhancer and promoter regions, as well as the coding region comprising exon1–intron1–exon2–intron2, were PCR amplified using bisulfite specific primer (primer sequences are described in Supplementary Table S3). PCR products were purified (PCR-Purification Kit; Qiagen) and Sanger sequenced. The HOXA3 locus was used as positive control.
Reverse Transcription and Radioactive Semiquantitative PCR
Reverse transcription-PCR (RT-PCR) was performed as previously described in detail (44), using sequence-specific primers. PCR step was performed with radiolabeled primers (Supplementary Table S4) and cycle numbers chosen to provide a signal that is linear with respect to input RNA. Quantification was done by densitometry using a Typhoon PhosphorImager (Amersham Biosciences).
CD19 Mini-Gene Cross-Linking and Pull-Down Assays
The region expanding exon2 and 20nt of its flanking introns was synthesized (Genewiz) and cloned into pBSKii+ (Supplementary Table S5). This mini-gene was transcribed in vitro, radioactively labeled with CTP-32P, and incubated with nuclear lysates from NALM-6 B-ALL cells for 30 minutes at 30°C. Exposure to UV light (254 nm) induced covalent cross-linking of nuclear proteins to RNA, as previously described (45). Immunoprecipitation of cross-linked RNA/protein complexes using antibodies specific for splicing factors (Supplementary Table S6) was performed as previously described (46, 47).
Western Blotting and Coimmunoprecipitation
Whole-cell protein lysates were obtained in RIPA buffer. Protein concentrations were estimated by Bio-Rad colorimetric assay. Immunoblotting was performed as previously described by loading 10 μg of protein onto 7.5% PAGE gels (48). Signals were detected by enhanced chemiluminescence (Pierce) or by Odyssey Infrared Imager (LI-COR Biosciences). Representative blots are shown. Antibodies used are listed in Supplementary Table S7. Coimmunoprecipitations were performed in whole-cell protein lysates from 15 million cells in 500 μL of nondenaturing buffer (150 mmol/L NaCl, 50 mmol/L Tris-pH8, 1% NP-10, 0.25% sodium deoxycholate) and 10 μL of kinase-specific antibodies. After overnight incubation at 4°C, 50 μL of Protein A agarose beads (Invitrogen) were added and incubated for 1 hour at 4°C. Beads were washed 3 times with nondenaturing buffer, and proteins were eluted in Laemmli sample buffer, boiled, and loaded onto PAGE gels.
Deglycosylation Assay
Whole-cell protein lysates were obtained using a nondenaturing buffer (150 mmol/L NaCl, 50 mmol/L Tris-pH8, 1% NP-10, 0.25% sodium deoxycholate) and treated with deglycosylation mix (New England BioLabs; #P6039S) following the manufacturer's instructions in the presence of protease and phosphatase inhibitors (Thermo Scientific; #78446). Control and deglycosylated lysates were loaded onto 8% PAGE gels for Western blot analysis.
Reverse Transcription, Real-Time Quantitative PCR (qRT-PCR), and PCR
Total RNA was isolated using TRIzol reagent (Invitrogen). cDNAs were prepared with random hexamers using the High Capacity cDNA RT Kit (Life Technologies). CD19 mRNA isoforms were visualized in 1% agarose gels after semiquantitative PCR amplification of cDNA using Platinum Taq-polymerase (Invitrogen) following the manufacturer's instructions. Primers used for each CD19 mRNA isoform and expected amplicon sizes are listed in Supplementary Table S8. When required, individual bands were gel-purified (QIAquick Gel Extraction Kit; Qiagen) and Sanger sequenced. qRT-PCR was performed using PowerSYBR Green PCR Master Mix (Life Technologies) and gene-specific oligo pairs (Supplementary Table S9). Reactions were performed on an Applied Biosystems Viia7 machine and analyzed with Viia7 RUO software (Life Technologies).
RNA-seq
RNA-seq reads were aligned using STAR version 2.4.0b (49) with a custom index based on the hg19 reference genome and a splice junction database consisting of all RefSeq isoforms supplemented with the exon 1–3 (Δex2) and exon 4–7 (Δex5–6) junctions for CD19. Aligned read counts per gene were computed using the htseq-count software with “–mode = intersection-strict” and normalized to gene RPKM by the formula: 109 × (read count aligning to gene)/[(mRNA length in bp) × (total aligned read count over all genes)].
WGS/WES Bioinformatic Processing and Point Mutation and LOH Analysis
Read alignment to the hg19 human reference genome for WES and whole-genome sequencing (WGS) samples was performed using the BWA v0.7.7 algorithm with default parameters. Unbiased point mutation calling was performed using SAMtools and BEDtools v0.1.18, and the aligned sequence at CD19 was further manually reviewed in the Integrative Genomics Viewer (IGV) in order to detect subclonal mutations in CD19. For genome-wide LOH analysis based on WES samples, B allele fractions (BAF) were computed for all common germline SNPs in the Single Nucleotide Polymorphism Database (dbSNP) build 142 (obtained from the UCSC Genome Browser) using SAMtools Mpileup v0.1.18. Genomic BAF profiles were then visualized in the R statistical programming language.
AVISPA Splicing Predictions
To find putative regulators of CD19 exon 2 and exon 5–6 skipping, we used AVISPA, which not only predicts if a cassette exon is alternatively spliced, but also gives a list of putative regulatory motifs that contribute to this splicing outcome (18). Hg19 coordinates were extracted for exons 1 through 3 to define the exon 2 triplet. Because AVISPA currently handles only single cassette exon events as inputs, coordinates for exons 4 through 7 were extracted to define two overlapping cassette exon triplets for the tandem skipping of exons 5 and 6 (i.e., an exon 4, 5, 6 triplet and an exon 4, 5, 7 triplet). These three triplets were run, and the top motifs and predicted associated splicing factors for the alternative versus constitutive splicing prediction step were compared. These top motifs were defined by their normalized feature effect, described in ref. 18. Briefly, this value represents the effect on splicing prediction outcome if a motif is removed in silico, normalized by the total effects observed from removing each of the top features in this way. This method has been used previously to detect and experimentally verify novel regulators of cassette exon splicing (17).
MAJIQ and VOILA Splicing Analysis
In order to identify and visualize splicing variations in CD19 from RNA-seq, we applied the MAJIQ and VOILA software (17). Briefly, STAR (49) was run to map the RNA-seq reads. Next, MAJIQ used the junction-spanning reads detected by STAR to construct a splice graph of CD19 and quantitate the PSI of the alternative exons. Finally, the VOILA visualization package was used to plot the resulting splice graph, the alternative splicing variants, and the violin plots representing the PSI estimates.
Disclosure of Potential Conflicts of Interest
S.L. Maude is a consultant/advisory board member for Novartis Pharmaceuticals. S. Gill and S.F. Lacey report receiving commercial research grants from Novartis. J.J. Melenhorst reports receiving commercial research support from Novartis. S.A. Grupp reports receiving a commercial research grant and other commercial research support from Novartis; has ownership interest (including patents) in University of Pennsylvania; and is a consultant/advisory board member for Novartis. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: E. Sotillo, D.M. Barrett, K.L. Black, A. Bagashev, G. Wu, R. Sussman, D. Allman, C. Mackall, K.W. Lynch, J.M. Maris, S.A. Grupp, A. Thomas-Tikhonenko
Development of methodology: E. Sotillo, D.M. Barrett, K.L. Black, A. Bagashev, G. Wu, S. Gill, D. Allman, J.M. Maris
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): E. Sotillo, D.M. Barrett, K.L. Black, A. Bagashev, G. Wu, C. Lanauze, M. Ruella, N.M. Martinez, C.T. Harrington, E.Y. Chung, J. Perazzelli, T.J. Hofmann, S.L. Maude, S.F. Lacey, J.J. Melenhorst, E. Jacoby, J.M. Maris, S.A. Grupp
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): E. Sotillo, D.M. Barrett, K.L. Black, A. Bagashev, D. Oldridge, G. Wu, M. Ruella, M.R. Gazzara, P. Raman, A. Barrera, S. Gill, S.F. Lacey, J.J. Melenhorst, C. Mackall, Y. Barash, K.W. Lynch, J.M. Maris, S.A. Grupp, A. Thomas-Tikhonenko
Writing, review, and/or revision of the manuscript: E. Sotillo, D.M. Barrett, K.L. Black, A. Bagashev, M. Ruella, N.M. Martinez, S.L. Maude, P. Raman, S.F. Lacey, J.J. Melenhorst, T. Fry, C. Mackall, K.W. Lynch, J.M. Maris, S.A. Grupp, A. Thomas-Tikhonenko
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): E. Sotillo, T.J. Hofmann, S.A. Grupp
Study supervision: E. Sotillo, R. Sussman, D. Allman, S.A. Grupp, A. Thomas-Tikhonenko
Other (provided samples): T. Fry
Acknowledgments
The authors thank the members of their laboratories and Drs. Carl June and David Porter (University of Pennsylvania) for stimulating discussions.
Grant Support
This work was supported by the V Foundation for Cancer Research, the William Lawrence and Blanche Hughes Foundation, the Leukemia & Lymphoma Society, the Alex's Lemonade Stand Foundation, the NIH (grants R01 CA102709, T32 GM007229 “Training Program in Cell and Molecular Biology,” and T32 HL007439 “Hematology Clinical Research Training Program”), the Pennsylvania Department of Health (to S.A. Grupp), and Stand Up To Cancer St. Baldrick's Pediatric Dream Team Translational Research Grant SU2C-AACR-DT1113. Stand Up To Cancer is a program of the Entertainment Industry Foundation administered by the American Association for Cancer Research.