Abstract
Oncogenic gene fusions drive many human cancers, but tools to more quickly unravel their functional contributions are needed. Here we describe methodology permitting fusion gene construction for functional evaluation. Using this strategy, we engineered the known fusion oncogenes, BCR-ABL1, EML4-ALK, and ETV6-NTRK3, as well as 20 previously uncharacterized fusion genes identified in The Cancer Genome Atlas datasets. In addition to confirming oncogenic activity of the known fusion oncogenes engineered by our construction strategy, we validated five novel fusion genes involving MET, NTRK2, and BRAF kinases that exhibited potent transforming activity and conferred sensitivity to FDA-approved kinase inhibitors. Our fusion construction strategy also enabled domain-function studies of BRAF fusion genes. Our results confirmed other reports that the transforming activity of BRAF fusions results from truncation-mediated loss of inhibitory domains within the N-terminus of the BRAF protein. BRAF mutations residing within this inhibitory region may provide a means for BRAF activation in cancer, therefore we leveraged the modular design of our fusion gene construction methodology to screen N-terminal domain mutations discovered in tumors that are wild-type at the BRAF mutation hotspot, V600. We identified an oncogenic mutation, F247L, whose expression robustly activated the MAPK pathway and sensitized cells to BRAF and MEK inhibitors. When applied broadly, these tools will facilitate rapid fusion gene construction for subsequent functional characterization and translation into personalized treatment strategies. Cancer Res; 77(13); 3502–12. ©2017 AACR.
Introduction
Cancer genome profiling efforts by large consortia such as The Cancer Genome Atlas (TCGA) are cataloging the complex genomic landscape of diverse tumors, which are comprised of numerous somatically acquired genetic alterations ranging from point mutations to ploidy change. Similar to activating mutations in oncogenes, chromosomal rearrangements also lead to oncogene activation through formation of gene fusion transcripts by, for example, increasing oncogene expression via a fused hyperactive transcription promoter or creation of chimeric hypermorphic or neomorphic alleles (1). Indeed, oncogenic fusion genes represent an important class of cancer driver aberrations, some of which have been exploited clinically for cancer therapy. For example, the functional role of BCR-ABL1 in promoting chronic myeloid leukemia led to successful therapies incorporating ABL inhibitors such as imatinib and dasatinib (2, 3). Similarly, use of ALK inhibitors crizotinib and ceritinib has significantly improved clinical outcome in non–small cell lung cancers driven by ALK fusions (4, 5).
The discovery of gene fusions has been accelerated by advances in next-generation sequencing (NGS) technologies (6). While the overall frequency of recurrent fusion transcripts is lower than activating mutations in oncogenes, the oncogenic role of individual fusion genes is suggested by their presence in multiple tumor types as well as the anti-correlation between their presence and that of cancer driver mutations in known oncogenes (7). More importantly, several recent reports describing the oncogenic behavior and therapeutic response of tumors driven by extremely rare fusions highlight their clinical impact. For example, individual cases of myeloid neoplasms driven by fusions involving JAK2 and FLT3 are sensitive to JAK inhibitor (ruxolitinib; ref. 8) and tyrosine kinase inhibitor (sorafenib; ref. 9), respectively. Likewise, we recently reported an oncogenic fusion involving the RET kinase in a single medullary thyroid carcinoma patient whose activity is highly sensitive to multiple tyrosine kinase inhibitors (10). Together, these examples highlight the importance of identifying the subset of rare, oncogenic gene fusions and assessing their sensitivity to therapeutics.
The functional interrogation of fusion genes is complicated given their large number, inability to accurately predict those with driver activity and technical roadblocks preventing efficient fusion gene construction for biological assays. To address these challenges, we report here a method enabling rapid and accurate fusion gene construction using a multi-fragment, recombineering-based strategy. We used this approach to construct known oncogenic fusion genes BCR-ABL1, EML4-ALK, and ETV6-NTRK3, all of which exhibited strong driver activity consistent with their role in cancer. We next scaled our fusion gene construction strategy to successfully build a pilot set of 20 fusion genes identified from TCGA datasets. Functional validation assays revealed that five of these fusion genes, which contained portions of MET, NTRK2, and BRAF kinases, exhibited robust transforming activity and marked responsiveness to inhibitors targeting their activated pathways. To illustrate another use of our fusion gene cloning strategy that leverages its versatility and modular design, we performed domain-function studies of BRAF fusion genes by differentially recombining N-terminal segments/domains of BRAF onto BRAF's C-terminal kinase domain. Data resulting from this work support previous reports indicating that the transforming activity by BRAF fusion genes results from truncation-mediated loss of inhibitory domains located within the N-terminus of BRAF (11–14). Because gene mutations residing within this inhibitory domain might serve as a means to activate BRAF in cancer, we leveraged the modular design of our construction methodology to fuse onto BRAF's kinase domain a set of inhibitory domains, each containing individual patient mutations, to screen for those capable of attenuating kinase inhibition. Using this approach, we identified an oncogenic mutation, F247L, whose expression robustly activates the MAPK pathway and sensitizes cells to inhibitors of BRAF and MEK.
Materials and Methods
Fusion gene construction
The DNA sequences of positive control fusion genes (BCR-ABL1, EML4-ALK, and ETV6-NTRK3) and the 20 uncharacterized fusion genes (Supplementary Tables S1 and S2) were obtained from The Cancer Genome Atlas (TCGA) and other published sources (7, 15–17). PCR templates were from sequence-verified open reading frame (ORF) collections by the ORFeome collaboration (18–20), Mammalian Gene Collection (21), and commercial ORF sources (Life Technologies). On the basis of the sequence of each fusion arm, primers were designed to contain recombination sequences (B1/B2/B2r/B4) followed by 18-nucleotide ORF-specific sequence (Supplementary Table S3). As illustrated in Fig. 1A, the forward primer for ORF 1 carried a B1 site (5′-GGGGACAACTTTGTACAAAAAAGTTGGC) while the reverse primer carried a B2 site (5′-GGGGACCACTTTGTACAAGAAAGCTGGGT). As to ORF 2, the forward primer carried a B2r site (5′-GGGGACCCAGCTTTCTTGTACAAAGTGGTTA) while the reverse primer carried a B4 site (5′-GGGGACAACTTTGTATAGAAAAGTTGGGTG). Depending on the given ORF sequence combined with each recombination site sequence, it is possible that the fusion primer would introduce an unwanted, in frame stop codon that would ultimately reside between the left and right fusion gene arms. One prevents this by manually checking the fusion primer sequence for stop codons, which would be removed by altering the wobble position, or simply substituting the following B2 site and B2r site sequences that eliminates the possibility of stop codons: B2 site (5′-GGGGACCACTTTGTACAAGAAAGCTGGGA); B2r site (5′- GGGGACCCAGCTTTCTTGTACAAAGTGGTTC). The resulting PCR products were incorporated into compatible pDONR vectors (P1/P2 for Fusion fragment 1, P2r/P4 for Fusion fragment 2) through BP recombination (Life Technologies) following the manufacturer's recommendations. The resulting Fusion fragment 1 will be flanked with L1/L2 recombination sites, while Fusion fragment 2 will be flanked with R2/L4 recombination sites. These products were subsequently transferred into pLenti-EF1α-DEST and pHAGE-EF1α-DEST vectors containing R1/R4 sites through multi-site LR recombination reaction (Life Technologies). The reaction mixtures were transformed into STBL3 (Life Technologies) competent bacteria. The 5′ and 3′ pDONR vectors as well as pLenti-EF1α-DEST and pHAGE-EF1α-DEST vectors required for fusion gene construction/expression are available from Addgene (www.addgene.org).
Multifragment recombineering of fusion genes. A, Schematic illustration of fusion gene construction. ATG, translation start site; B1/P1, B2/P2, B2r/P2r, B4/P4, recombination sites for BP recombination; L1/R1, L2/R2, L4/R4, recombination sites for LR recombination; pFor and pRev, PCR detection primers. B, Ba/F3 cell survival assay for BCR-ABL1, EML4-ALK, and ETV6-NTRK3 seven days following IL3 depletion (mean luminescence; error bars, SD; N = 3). C, Immunoblots of BCR-ABL1 and EML4-ALK expression in Ba/F3. Arrow, the correct size of BCR-ABL1. D, PCR detection of the indicated fusion transcripts from Ba/F3 RNA/cDNA extracts. B, fusion DNA backbone (positive control); −, cDNA from GFP-expressing cells as negative control. E, Dose-dependent survival assays of Ba/F3 cells expressing BCR-ABL1 and EML4-ALK treated with imatinib and crizotinib, respectively, for 72 hours (mean percentage of cell survival; error bars, SD; N = 4). F, Endpoint volumes (day 59 postinjection) of xenograft tumors by HMLER cells expressing ETV6-NTRK3 (N = 8) and GFP control (N = 15). Horizontal bars denote mean volumes; error bars, SD. All P values calculated by t test; *, P < 0.05; **, P < 0.01; ****, P < 0.0001.
Multifragment recombineering of fusion genes. A, Schematic illustration of fusion gene construction. ATG, translation start site; B1/P1, B2/P2, B2r/P2r, B4/P4, recombination sites for BP recombination; L1/R1, L2/R2, L4/R4, recombination sites for LR recombination; pFor and pRev, PCR detection primers. B, Ba/F3 cell survival assay for BCR-ABL1, EML4-ALK, and ETV6-NTRK3 seven days following IL3 depletion (mean luminescence; error bars, SD; N = 3). C, Immunoblots of BCR-ABL1 and EML4-ALK expression in Ba/F3. Arrow, the correct size of BCR-ABL1. D, PCR detection of the indicated fusion transcripts from Ba/F3 RNA/cDNA extracts. B, fusion DNA backbone (positive control); −, cDNA from GFP-expressing cells as negative control. E, Dose-dependent survival assays of Ba/F3 cells expressing BCR-ABL1 and EML4-ALK treated with imatinib and crizotinib, respectively, for 72 hours (mean percentage of cell survival; error bars, SD; N = 4). F, Endpoint volumes (day 59 postinjection) of xenograft tumors by HMLER cells expressing ETV6-NTRK3 (N = 8) and GFP control (N = 15). Horizontal bars denote mean volumes; error bars, SD. All P values calculated by t test; *, P < 0.05; **, P < 0.01; ****, P < 0.0001.
BRAF mutagenesis
Site-directed mutagenesis was performed on the basis of the HiTMMoB platform (22, 23). Mutagenesis primers were designed as listed in Supplementary Table S4. Mutant fragments and C-terminal BRAF kinase domain (BRAF-ex9) in pDONR vectors were transferred into pHAGE-EF1α-DEST vectors containing R1/R4 sites through multi-site LR recombination reaction (Life Technologies).
Cell culture and transduction
All cell lines were propagated at 37°C and 5% CO2 in humidified atmosphere. MCF10A cells were obtained from the ATCC and cultured as described previously (24). HMLER cells were provided by S. Mani (MD Anderson Cancer Center, Houston, TX) and maintained in MEGM medium (Lonza) without serum and antibiotics. Parental Ba/F3 cells were cultured in RPMI1640 medium supplemented with final concentration of 5% FBS and 2.5 ng/mL recombinant mouse IL3 (R&D Systems). Lentivirus production and cell transduction were described earlier (25, 26). Cells lines were fingerprinted prior to use on April 21, 2015 by the MD Anderson Cancer Center Characterized Cell Line Core using STR testing platform. Ba/F3 is a mouse-originated cell line, thus STR testing could not be performed.
Ba/F3 cell viability and inhibitor assays
As described previously (25, 27), Ba/F3 cell viability was determined at 7 days after IL3 depletion. Cells were treated with DMSO or respective inhibitors at the indicated concentrations for 72 hours, and cell viability was determined using CellTiter-Glo (Promega). All inhibitor compounds were purchased from Selleck Chemicals.
MCF-10A anchorage-independent growth assay
Soft agar assays were performed in 6-well plates in triplicate. First, bottom layers were prepared at 0.8% noble agar (Affymetrix, Inc.) with complete MCF-10A growth medium. After solidification, 1 × 104 cells were mixed with 0.45% agar in complete growth medium and laid on top of the bottom layer. Two milliliters medium was added in each well after 3 days and medium was refreshed every 3 days. Colonies were counted 2 weeks after seeding.
Immunoblotting
The following antibodies were used to detect protein expression: c-Abl (Cell Signaling Technology), ALK (Cell Signaling Technology), Raf-B (Santa Cruz Biotechnology), MET (Cell Signaling Technology), TrkB (Abcam), Phospho-Stat1 (Y701; Cell Signaling Technology), Stat1 (Cell Signaling Technology), phospho-Stat3 (Y705; Cell Signaling Technology), phospho-AKT (S473 & T308; Cell Signaling Technology), phospho-ERK1/2 (T202/Y204; Cell Signaling Technology), ERK1/2 (Cell Signaling Technology) Vinculin (Cell Signaling Technology), and GAPDH (Santa Cruz Biotechnology).
Immunofluorescence
MCF-10A cells and those expressing MET fusions and wild-type MET were seeded at 2 × 104 cells/well on Millicell EZ Slide (EMD Millipore). Cells are fixed in 2% formalin and permeabilized by 0.5% Triton X-100/PBS. The following antibodies were used to detect protein expression and localization: MET (Cell Signaling Technology) and GM130 (Golgi marker; Abcam). Slides were mounted by SlowFade Gold Antifade Mountant with DAPI (Thermo Fisher Scientific) to label DNA.
Quantitative PCR and fusion detection PCR
Total RNA was isolated from transduced cells for cDNA synthesis using SuperScript IV First-Strand Synthesis System (Life Technologies) as described previously (26, 28). For qPCR, coding regions within wild-type gene transcript were amplified using gene-specific primers (Supplementary Table S5) and SYBR Green PCR Master Mix (Life Technologies). Expression levels were normalized to mouse β-actin and comparative cycle threshold method was used to quantify mRNA copy number. For fusion detection PCR, coding regions in fusion transcripts were amplified using primers (Supplementary Table S5) annealing to either fusion fragment by PCR and the products were run on 1.5% agarose gel.
Animal studies
All studies using mice were performed in accordance with our IACUC-approved animal protocol (AN-5428) at Baylor College of Medicine. For xenograft tumor assays, primary, nontransformed cells of the defined lineages (106 cells per injection site) were suspended in a 1:1 Hank balanced salt solution (Life Technologies) and Matrigel (BD Biosciences) and injected into female athymic mice (CrTac:NCr-Foxn1nu;4–6 weeks in age) subcutaneously at bilateral flanks (28). Mice were monitored twice a week and tumors were measured and calculated by length × width2/2. The assays were terminated (assay endpoint) once the first animal presents a tumor at maximal burden as allowed in our animal protocol.
Results
Fusion gene construction
To address technical challenges related to fusion gene cloning, we developed an approach incorporating a Gateway-based (Life Technologies) phage recombination strategy permitting highly efficient and accurate fusion gene construction (Fig. 1A). This strategy incorporates use of two PCR primer sets, each containing requisite recombination sites (attB1/B2/B2r/B4) adjacent to nucleotides complementary to ORF sequences comprising the target fusion fragment. These primers are entered into a high-fidelity PCR reaction using total cellular cDNA or ORF gene clones as template, and the resulting two PCR amplicons are subsequently recombined into separate entry clones by recombination. The preparation of separate entry clones permits simultaneous assembly of control fusion clones containing optional fluorescent proteins, regulatory elements, etc. (see GFP fusion studies below). Each set of two entry clones are compatible with secondary multi-fragment recombination into a variety of compatible destination vectors, permitting fusion gene expression in desired cell models. Recombination provides fusion fragment assembly in the appropriate order to model gene breakpoints while maintaining fragment reading frames. The strategy results in a fixed 27-nucleotide linker between the two fusion fragments, which when translated leads to an in-frame, 9-amino acid segment (TQLSCTKWL). Inclusion of the short linker did not impact functional studies using positive control fusion constructs (see positive control studies below). Moreover, the modularity of fusion genes whose activity are commonly based on gain/loss of protein functional domains, as well as the observation that breakpoint locations vary widely among recurrent fusion oncoproteins, further supports the limited impact of the linker on the ability to assess the activity of fusion genes.
Functional validation of known fusion oncogenes
Use of sequence-validated ORF collections made available by the ORFeome collaboration (20) and others (21) as fragment PCR template permits the scaling of fusion gene construction. As proof-of-concept, we leveraged the ORFeome collection as well as commercial ORF sources to build three known fusion oncogenes based on published sequences: BCR-ABL1, EML4-ALK, and ETV6-NTRK3. To validate activity of these fusion genes, we first employed murine pro-B Ba/F3 cells, which die in the absence of exogenous IL3 (29), to quantify the ability of driver fusion genes to rescue cell survival and proliferation in the absence of IL3. The ability to assay transfer of Ba/F3 addiction from IL3 to oncogenes has been used by us (23, 25, 27) and others (30) to investigate activity and therapeutic sensitivity by kinase proteins. Viral delivery of BCR-ABL1 promoted robust 245-fold (P < 0.0001) Ba/F3 cell proliferation in the absence of IL3 compared with undetectable growth by parental cells and those expressing GFP (negative control), and this activity was at a level comparable with parental BCR-ABL1 cloned directly from a patient sample (Fig. 1B). Similarly, expression of EML4-ALK and ETV6-NTRK3 led to a 100- (P = 0.0217) and 159- (P = 0.0022) fold increase in Ba/F3 cell growth, respectively, compared with GFP-expressing cells (Fig. 1B). We next confirmed BCR-ABL1 and EML4-ALK fusion gene expression in Ba/F3 cells by immunoblot analysis (Fig. 1C). Given our difficulty identifying a suitable antibody for detection of ETV6-NTRK3, we confirmed mRNA transcription of the chimeric gene by PCR amplifying over the 27-nucleotide linker from cDNA preparations (Fig. 1D) using DNA primers annealing to either arm of the fusion transcript (pFor and pRev, Fig. 1A). This simple transcript detection strategy can be used to confirm fusion gene construct expression across transduced cell lines.
Ba/F3 has been widely used to evaluate kinase inhibitors (30). Therefore, we next used BCR-ABL1- and EML4-ALK–expressing Ba/F3 cells to measure response to clinically approved inhibitors of ABL1 and ALK. Treatment of Ba/F3 cells addicted to BCR-ABL1 and EML4-ALK exhibited robust sensitivity to ABL inhibitor imatinib (IC50 = 0.3012 μmol/L) and ALK inhibitor crizotinib (IC50 = 0.0097 μmol/L), respectively, compared with parental Ba/F3 cells (Fig. 1E).
ETV6-NTRK3 is a recurrent event in secretory breast carcinomas (31). We chose to evaluate ETV6-NTRK3 activity in primary, nontumorigenic human mammary epithelial cells (HMEC) modified to express the telomerase catalytic subunit, SV40 large T and small t antigens and H-RasV12 (HMLER; ref. 32). Expression of ETV6-NTRK3 in HMLER cells led to robust tumor formation (N = 7/8; P = 0.0013) when implanted into athymic mice, whereas control cells stably expressing GFP failed to form significant numbers of tumors (N = 1/15; Fig. 1F).
Construction of previously uncharacterized fusion genes
We resourced recently published literature (7, 15–17) on TCGA datasets to select 20 functionally uncharacterized fusion genes based on mRNA expression profiles and frame status (Supplementary Tables S1 and S2). We scaled fusion PCR reactions to a 96-well format to amplify the unique individual 5′ and 3′ DNA fragments corresponding to each of the 20 fusion genes using the required fusion primers (Fig. 1A; Supplementary Table S3). Each PCR reaction was used for BP recombination into corresponding fusion donor clones (5′ and 3′ fragments recombined into pDONR-P1/P2 and pDONR-P2r/P4, respectively; Fig. 1A) by direct transfer in a 96-well format. We selected, on average, 3 bacterial colonies following each donor recombineering reaction for DNA sequencing, which revealed a 94.7% success rate (both fusion arms) for proper PCR synthesis and recombination (Supplementary Table S6). A sequence-verified plasmid isolate representing each 5′ and 3′ fusion arm was combined with destination vector for multi-fragment LR recombination (Fig. 1A). Sequencing an average of 3.9 colonies per LR reaction indicated an 86.7% success rate for proper recombination across all fusion clones (Supplementary Table S6). In total, we achieved 100% of fusion clones attempted in a single fusion construction run with an overall efficiency of 91.5% considering correctly sequenced clones analyzed throughout each step of the construction strategy.
We next sought to examine each of the 20 fusion genes for driver activity; however, systematic testing for driver activity across groups of genes is confounded by differences in their encoded protein function and differing roles in the recognized hallmarks of cancer (33). Therefore, there are no true “generalizable” driver testing systems applicable to detecting oncogenic or oncogene effector activity across all gene and cancer lineage types. Considering this caveat, we again employed the Ba/F3 cell model acknowledging a potential high false negative rate for detecting activity for bona fide fusion gene drivers of cancer. Nevertheless, viral delivery of the 20 fusion genes revealed five (BAIAP2L1-MET, TFG-MET, AFAP1-NTRK2, SQSTM1-NTRK2, and FAM114A2-BRAF) that promoted robust cell proliferation in the absence of IL3 compared with GFP-expressing control cells (described below). Immunoblotting and fusion transcript PCR using Ba/F3 extracts confirmed expression of all 20 fusion genes (Supplementary Figs. S1–S3 and below).
MET fusion genes
Among the 20 uncharacterized fusion genes selected for construction were four that involved the MET receptor tyrosine kinase (Fig. 2A; ref. 34). BAIAP2L1-MET and TFG-MET, which were identified to be expressed in kidney papillary cell carcinoma (TCGA-BQ-7049) and thyroid carcinoma (TCGA-FK-A3S3), respectively, strongly enhanced Ba/F3 cell survival and proliferation (79- and 213-fold compared with GFP-expressing and wild-type MET-expressing cells (P < 0.0001; Fig. 2B). In contrast, both variants of CAPZA2-MET failed to stimulate Ba/F3 growth (Fig. 2B) despite encoding the complete tyrosine kinase domain of MET and being expressed in Ba/F3 (Supplementary Fig. S2A). Control assays revealed that expression of the individual gene fragments encoding N-terminal gene partners (BAIAP2L1 and TFG) of the active MET fusion genes failed to promote Ba/F3 growth (Supplementary Fig. S2B) despite their verified expression (Supplementary Fig. S2C). We next assayed MET fusion driver activity using normal MCF-10A breast epithelial cells (35), which are widely used for anchorage-independent growth assays to assess oncogene activity. Expression of BAIAP2L1-MET and TFG-MET significantly increased MCF-10A colony formation (39- and 43-fold compared with GFP-expressing cells, respectively; P < 0.0001) similar to PIK3CAH1047R oncogenic control (Fig. 2C; ref. 23). Both variants of CAPZA2-MET failed to stimulate transformation and colony growth in MCF-10A (Fig. 2C) similar to our findings with Ba/F3. Expression of wild-type MET moderately increased MCF-10A cell colony formation by 9-fold (P = 0.005; Fig. 2C) while it failed to relieve IL3 dependency of Ba/F3 cells (Fig. 2B). Immunoblotting of MCF-10A extracts confirmed expression of each fusion gene and revealed heightened phosphorylation of ERK1/2 (T202/Y204), STAT1 (Y701), and STAT3 (Y705) in cells expressing BAIAP2L1-MET and TFG-MET compared with control and both CAPZA2-MET variants (Supplementary Fig. S2D and S2E). Finally, we performed Ba/F3 dose–response assays to examine the response of BAIAP2L1-MET and TFG-MET to crizotinib, which is known to inhibit the MET kinase (36, 37). Ba/F3 cells expressing BAIAP2L1-MET or TFG-MET exhibited marked sensitivity to crizotinib (IC50, BAIAP2L1-MET = 8.42 nmol/L; TFG-MET = 10.53 nmol/L) compared with control cells (IC50 = 675 nmol/L; Fig. 2D).
Oncogenic validation of MET fusions. A, Schematic illustration of MET fusion genes. B, Ba/F3 cell survival assay for MET fusions (mean luminescence; error bars, SD, N = 3) compared with positive control, BCR-ABL1. GFP, negative control. C, MCF-10A anchorage-independent colony formation assays for all MET fusions (mean colony count from 10 random areas; error bars, SD; N = 3). PIK3CAH1047R, positive control; GFP, negative control. D, Dose-dependent survival assays of Ba/F3 cells expressing BAIAP2L1-MET and TFG-MET treated with crizotinib for 72 hours (mean percentage of cell survival; error bars, SD; N = 4). E, MCF-10A cells expressing BAIAP2L1-MET, CAPZA2-MET-2, wild-type MET, and parental were immunostained for MET (red) and Golgi body marker GM130 (green). DNA was labeled with DAPI. Scale bar, 50 μm. All P values calculated by t test; ns, not significant; **, P < 0.01; ****, P < 0.0001.
Oncogenic validation of MET fusions. A, Schematic illustration of MET fusion genes. B, Ba/F3 cell survival assay for MET fusions (mean luminescence; error bars, SD, N = 3) compared with positive control, BCR-ABL1. GFP, negative control. C, MCF-10A anchorage-independent colony formation assays for all MET fusions (mean colony count from 10 random areas; error bars, SD; N = 3). PIK3CAH1047R, positive control; GFP, negative control. D, Dose-dependent survival assays of Ba/F3 cells expressing BAIAP2L1-MET and TFG-MET treated with crizotinib for 72 hours (mean percentage of cell survival; error bars, SD; N = 4). E, MCF-10A cells expressing BAIAP2L1-MET, CAPZA2-MET-2, wild-type MET, and parental were immunostained for MET (red) and Golgi body marker GM130 (green). DNA was labeled with DAPI. Scale bar, 50 μm. All P values calculated by t test; ns, not significant; **, P < 0.01; ****, P < 0.0001.
We next sought to investigate why CAPZA2-MET fusion genes failed to promote Ba/F3 growth and MCF-10A colony formation. The MET protein is normally processed by the Golgi apparatus and related protein processing machinery for deposition at the plasma membrane where it interacts with its ligand, hepatocyte growth factor (HGF; ref. 34). One explanation for CAPZA2-MET inactivity might relate to its inability to reach the cell surface due to altered localization by CAPZA2. To test this possibility, we performed immunofluorescence assays using an antibody specific to the C-terminus of MET. Staining of MCF-10A cells expressing active BAIAP2L1-MET and wild-type MET revealed MET localization at the Golgi apparatus and plasma membrane, whereas signal was undetectable in vector control (parental) cells, suggesting low endogenous MET expression (Fig. 2E). In contrast, cells expressing CAPZA2-MET exhibited MET localization at the nucleus (Fig. 2E), an observation consistent with our hypothesis. However, we cannot eliminate a possible oncogenic role of CAPZA2-MET protein in the nucleus that might be discernable with use of other models to assess oncogenic activities.
NTRK2 fusion genes
Our series of 20 fusion genes also included AFAP1-NTRK2 and SQSTM1-NTRK2, which were identified in low-grade glioma specimens (TCGA-HT-7680 and TCGA-DU-A76L, respectively; Fig. 3A). AFAP1-NTRK2 and SQSTM1-NTRK2 strongly enhanced Ba/F3 cell survival and proliferation (316- and 103-fold, respectively, compared with GFP-expressing cells; P < 0.0001; Fig. 3B) in the absence of IL3 at a level comparable with positive control, ETV6-NTRK3. In contrast, wild-type NTRK2 exhibited much weaker activity in the Ba/F3 model (7-fold, P = 0.306; Fig. 3B). We also observed no activity for N-terminal partners of the two NTRK2 fusions (AFAP1 and SQSTM1; Supplementary Fig. S3A) compared with GFP and BRAFV600E-positive control cells despite verified expression by qPCR (Supplementary Fig. S3B). We observed a similar trend using the MCF-10A colony formation assay, as expression of AFAP1-NTRK2 exhibited greater transforming activity (34-fold, P < 0.0001) than SQSTM1-NTRK2 (7-fold, P = 0.0216; Fig. 3C). Wild-type NTRK2 did not significantly promote colony formation (Fig. 3C) despite verified expression by qPCR (Supplementary Fig. S3C). Immunoblot analysis confirmed fusion gene expression and activation of MAPK signaling through elevated phosphorylation of ERK1/2 (T202/Y204; Supplementary Fig. S3D and S3E). Interestingly, the same anti-NTRK2 (Trk-B) antibody used to detect AFAP1-NTRK2 and SQSTM1-NTRK2 protein expression did not detect wild-type NTRK2 (Supplementary Fig. S3D and S3E) despite the fact all were expressed from the same vector type and qPCR analysis of RNA extracts from the same cells indicated similar transcript levels for all three constructs (Supplementary Fig. S3C). Finally, dose–response assays using Ba/F3 revealed both NTRK2 fusions sensitized cells to a pan-NTRK inhibitor, entrectinib (38) (IC50, AFAP1-NTRK2 = 4.57 nmol/L; SQSTM1-NTRK2 = 2.26 nmol/L; parental > 1 μmol/L; Fig. 3D).
Oncogenic validation of NTRK2 fusions. A, Schematic illustration of NTRK2 fusion genes. B, Ba/F3 cell survival assay for NTRK2 fusions (mean luminescence; error bars, SD; N = 3) compared with positive control, ETV6-NTRK3. GFP, negative control. C, MCF-10A anchorage-independent colony formation assays for NTRK2 fusions (mean colony count from 10 random areas; error bars, SD; N = 3). GFP, negative control. D, Dose-dependent survival assays of Ba/F3 cells expressing AFAP1-NTRK2 and SQSTM1-NTRK2 treated with entrectinib for 72 hours (mean percentage of cell survival; error bars, SD; N = 4). All P values calculated by t test; ns, not significant; *, P < 0.05; ****, P < 0.0001.
Oncogenic validation of NTRK2 fusions. A, Schematic illustration of NTRK2 fusion genes. B, Ba/F3 cell survival assay for NTRK2 fusions (mean luminescence; error bars, SD; N = 3) compared with positive control, ETV6-NTRK3. GFP, negative control. C, MCF-10A anchorage-independent colony formation assays for NTRK2 fusions (mean colony count from 10 random areas; error bars, SD; N = 3). GFP, negative control. D, Dose-dependent survival assays of Ba/F3 cells expressing AFAP1-NTRK2 and SQSTM1-NTRK2 treated with entrectinib for 72 hours (mean percentage of cell survival; error bars, SD; N = 4). All P values calculated by t test; ns, not significant; *, P < 0.05; ****, P < 0.0001.
BRAF fusion genes
We also constructed two BRAF fusion genes discovered in cutaneous melanoma and thyroid carcinoma specimens, AHCYL2-BRAF (TCGA-D3-A3C3) and FAM114A2-BRAF (TCGA-ET-A3BN), both of which contain an intact serine/threonine kinase domain encoded by BRAF (Fig. 4A). Examination in Ba/F3 revealed potent growth-promoting activity for FAM114A2-BRAF (336-fold increase compared with GFP control cells in the absence of IL3; P < 0.0001) similar to ATG7-BRAF previously reported by TCGA and other studies (Fig. 4B; refs. 7, 15). Unlike FAM114A2-BRAF, AHCYL2-BRAF provided no growth advantage in Ba/F3 at the same time point following IL3 removal (Fig. 4B). Consistent with this observation, immunoblot analysis of cell lysates from Ba/F3 cells expressing FAM114A2-BRAF and ATG7-BRAF revealed activation of MAPK signaling evidenced by phosphorylation of ERK1/2 (T202/Y204), and this response was absent in cells expressing AHCYL2-BRAF and GFP (Fig. 4C). Previous studies reported sensitivity of melanocytic tumors harboring BRAF fusion genes to targeted therapies (39, 40). Indeed, dose–response assays using Ba/F3 revealed marked sensitivity to clinically-approved BRAF and MEK inhibitors dabrafenib and trametinib, respectively, by Ba/F3 cells expressing FAM114A2-BRAF (dabrafenib, IC50 = 0.714 μmol/L; trametinib, 1.458 nmol/L) and ATG7-BRAF (dabrafenib, IC50 = 0.794 μmol/L; trametinib, 1.170 nmol/L) as predicted from the fusion's strong activation of MAPK signaling (Fig. 4D).
Oncogenic validation and domain-function studies of BRAF fusion genes. A, Schematic illustration of BRAF fusion genes. B, Ba/F3 cell survival assay for BRAF fusions (mean luminescence; error bars, SD; N = 3, respectively). ATG7-BRAF, positive control; GFP, negative control. C, Immunoblots of BRAF fusions expression and MAPK signaling activation in Ba/F3. D, Dose-dependent survival assays of Ba/F3 cells expressing FAM114A2-BRAF fusions treated with dabrafenib and trametinib for 72 hours (mean percentage of cell survival; error bars, SD; N = 4, respectively). ATG7-BRAF, positive control. E, Ba/F3 cell survival assay for the indicated full-length, wild-type genes (mean luminescence; error bars, SD; N = 3, respectively) compared with BRAFV600E (positive control). GFP, negative control. F, Ba/F3 cell survival assay for BRAF kinase domain (BRAF-ex9, exons 9–18; BRAF-ex11, exons 11–18) and corresponding GFP-BRAF-ex9/11 fusions with and without STOP codon following GFP (mean luminescence; error bars, SD; N = 3) compared with full-length, wild-type BRAF. BRAFV600E, positive control; GFP, negative control. G and H, Schematic illustration of construct structures (G) and activities in Ba/F3 cell survival assay: BRAF kinase domain (exons 9–18) fused to (i) full-length BRAF N-terminus, N-BRAF-ex9; (ii) fragment corresponding to BRAF AA100-345, N100-345-BRAF-ex9; (iii) kinase domain only, BRAF-ex9; (iv) GFP, GFP-BRAF-ex9 (H). Shown mean luminescence; error bars, SD; N = 3. All P values calculated by t test; **, P < 0.01; ***, P < 0.001; ****, P < 0.0001.
Oncogenic validation and domain-function studies of BRAF fusion genes. A, Schematic illustration of BRAF fusion genes. B, Ba/F3 cell survival assay for BRAF fusions (mean luminescence; error bars, SD; N = 3, respectively). ATG7-BRAF, positive control; GFP, negative control. C, Immunoblots of BRAF fusions expression and MAPK signaling activation in Ba/F3. D, Dose-dependent survival assays of Ba/F3 cells expressing FAM114A2-BRAF fusions treated with dabrafenib and trametinib for 72 hours (mean percentage of cell survival; error bars, SD; N = 4, respectively). ATG7-BRAF, positive control. E, Ba/F3 cell survival assay for the indicated full-length, wild-type genes (mean luminescence; error bars, SD; N = 3, respectively) compared with BRAFV600E (positive control). GFP, negative control. F, Ba/F3 cell survival assay for BRAF kinase domain (BRAF-ex9, exons 9–18; BRAF-ex11, exons 11–18) and corresponding GFP-BRAF-ex9/11 fusions with and without STOP codon following GFP (mean luminescence; error bars, SD; N = 3) compared with full-length, wild-type BRAF. BRAFV600E, positive control; GFP, negative control. G and H, Schematic illustration of construct structures (G) and activities in Ba/F3 cell survival assay: BRAF kinase domain (exons 9–18) fused to (i) full-length BRAF N-terminus, N-BRAF-ex9; (ii) fragment corresponding to BRAF AA100-345, N100-345-BRAF-ex9; (iii) kinase domain only, BRAF-ex9; (iv) GFP, GFP-BRAF-ex9 (H). Shown mean luminescence; error bars, SD; N = 3. All P values calculated by t test; **, P < 0.01; ***, P < 0.001; ****, P < 0.0001.
We next investigated why the AHCYL2-BRAF fusion gene failed to promote Ba/F3 growth despite its inclusion of an intact kinase domain. Previous studies suggest that the oncogenic potential of BRAF kinase fusion proteins can be attributed to the kinase domains but not to its N-terminal partner (40–42). Indeed, stable expression of individual gene fragments encoding the N-terminal gene partners of both BRAF fusion genes (AHCYL2 and FAM114A2) failed to promote Ba/F3 growth (Fig. 4E) despite verified expression (Supplementary Fig. S4A). Expression of full-length BRAF, which contains the wild-type kinase domain, similarly failed to influence Ba/F3 growth in contrast to expression of oncogenic BRAFV600E (Fig. 4E). We complemented these studies by expressing BRAF exons included in active BRAF fusion genes ATG7-BRAF and FAM114A2-BRAF: exons 9–18 (BRAF-ex9; as encoded by ATG7-BRAF) and exons 11–18 (BRAF-ex11; as encoded by FAM114A2-BRAF). Expression of these BRAF fragments led to significant Ba/F3 growth (average 145-fold for BRAF-ex9/11) compared with control cells (BRAF-ex9, P = 0.0002; BRAF-ex11, P = 0.0032; Fig. 4F).
We next leveraged the versatility afforded by the modular design of our fusion cloning strategy (Fig. 1A) to examine the consequence of fusing active BRAF-ex9/11 to N-terminal sequence other than FAM114A2, AHCYL2, or ATG7. Expression of each kinase fragment as a fusion to monomeric GFP (GFP-BRAF-ex9/ex11) robustly stimulated Ba/F3 growth (average 350-fold, P < 0.0001), whereas expressing similar constructs engineered with a stop codon following GFP coding sequence (GFP-STOP-BRAF-ex9/11) did not enhance Ba/F3 growth as expected (Fig. 4F). Interestingly, immunoblot analysis and qPCR of outgrowth Ba/F3 indicated lower expression for GFP-BRAF proteins but higher activation of MAPK signaling compared with truncated BRAF-ex9/11 proteins (Supplementary Fig. S4B and S4C) consistent with their greater activity in Ba/F3 (Fig. 4F), suggesting that addition of protein sequence N-terminal to BRAF kinase domain may perturb and, in fact, increase kinase activity.
Our data along with the structure of the inactive AHCYL2-BRAF fusion gene, which preserves not only the kinase domain but also a significant portion of the N-terminus of BRAF, also suggests that BRAF fusion genes like ATG7-BRAF and FAM114A2-BRAF are activated through removal of an inhibitory component within the N-terminus of BRAF. This notion is supported by previous findings (11–13) indicating that RAF kinases are regulated by the N-terminal autoinhibitory component and truncation of the N terminus can lead to activation of RAF kinases. To examine this further, we again resourced the modular design of our fusion construction strategy to fuse segments of BRAF N-terminus onto the active C-terminal (exons 9–18) fragment of BRAF (Fig. 4G). Fusion of full-length N-terminus back to the kinase domain (N-BRAF-ex9) completely abolished BRAF-ex9 activity in Ba/F3 (Fig. 4H) compared with the truncated kinase (BRAF-ex9) and GFP fusion (GFP-BRAF-ex9). Fusion of an N-terminal BRAF fragment (AA100-345), which is present in the inactive fusion gene AHCYL2-BRAF, back to the kinase domain (N100-345-BRAF-ex9) similarly attenuated BRAF-ex9 activity in Ba/F3 (Fig. 4H; versus GFP-BRAF-ex9, P < 0.0001) and showed no evidence of MAPK activation in early passage Ba/F3 cell lysates (Supplementary Fig. S4D) thus supporting earlier findings (11–13) that a portion of BRAF's N-terminus serves to inactivate kinase activity.
BRAF mutation studies
The presence of an inhibitory domain located within BRAF's N-terminus raises the possibility that mutation of this region might serve as another mechanism for BRAF activation during tumorigenesis. To investigate this possibility, we mined cBioPortal for Cancer Genomics (43, 44) for tumors and cell lines that contained a missense mutation within the coding sequence corresponding to AA100-345 but devoid of V600 hotspot mutations (Supplementary Table S4). Among all 1,438 missense mutations identified across 146 tumors, our analysis revealed 11 cases corresponding to 12 individual mutations, which were engineered into the N-terminal AA100-345 inhibitory fragment used for the studies in Fig. 4G and H.
We once again leveraged the modular design of our fusion gene construction strategy to fuse each mutated N-terminal fragment to the C-terminal BRAF kinase domain (BRAF-ex9), followed by testing all mutated fragments for ability to attenuate kinase suppression mediated by A100-345 in Ba/F3. Of the 12 mutations examined, only one (F247L; identified in TCGA case TCGA-AG-A002) promoted Ba/F3 growth in the absence of IL3 (184-fold compared with N100-345-BRAF-ex9, P < 0.0001) to a level comparable with GFP (GFP-BRAF-ex9) and V600E fusion genes [N-BRAF-ex9 (V600E); Fig. 5A and B]. Immunoblot analysis indicated that the F247L mutation in the fusion construct [N-BRAF-ex9 (F247L)] restored activation of MAPK signaling through phosphorylation of ERK1/2 (T202/Y204; Fig. 5C). To confirm these findings, we engineered the F247L mutation into full-length BRAF. Expression of BRAFF247L in Ba/F3 led to a significant IL3-independent growth (116-fold compared with GFP-expressing cells, P = 0.0011), albeit to a lesser extent compared with BRAFV600E (Fig. 5D), and immunoblot analysis confirmed MAPK activation by both BRAFF247L and BRAFV600E (Fig. 5E). Importantly, BRAFF247L sensitized Ba/F3 cells to both dabrafenib (IC50 = 0.6646 μmol/L) and trametinib (IC50 = 0.5412 nmol/L; Fig. 5F), suggesting potential use of these therapeutics in patients whose tumor harbor the BRAFF247L mutation.
BRAF fusion modularity and mutation studies. Schematic illustration (A) and activities (B) in Ba/F3 of BRAF kinase domain (Exons 9–18) fused to BRAF AA100-345 with or without F247L mutation (mean luminescence; error bars, SD; N = 3, respectively). GFP-BRAF-ex9 and N-BRAF-ex9 (V600E), positive control; N-BRAF-ex9 and GFP, negative control. C, Immunoblots of BRAF structural constructs expression and MAPK signaling activation in Ba/F3. D, Ba/F3 cell survival assay of full-length BRAFF247L; shown mean luminescence; error bars, SD; N = 3; BRAFV600E, positive control; GFP, negative control. E, Immunoblots of expression of full-length BRAFF247L, wild-type BRAF, and BRAFV600E in Ba/F3. F, Dose-dependent survival assays of Ba/F3 cells expressing full-length BRAFF247L treated with dabrafenib and trametinib for 72 hours (mean percentage of cell survival; error bars, SD; N = 4, respectively). All P values calculated by t test; **, P < 0.01; ****, P < 0.0001.
BRAF fusion modularity and mutation studies. Schematic illustration (A) and activities (B) in Ba/F3 of BRAF kinase domain (Exons 9–18) fused to BRAF AA100-345 with or without F247L mutation (mean luminescence; error bars, SD; N = 3, respectively). GFP-BRAF-ex9 and N-BRAF-ex9 (V600E), positive control; N-BRAF-ex9 and GFP, negative control. C, Immunoblots of BRAF structural constructs expression and MAPK signaling activation in Ba/F3. D, Ba/F3 cell survival assay of full-length BRAFF247L; shown mean luminescence; error bars, SD; N = 3; BRAFV600E, positive control; GFP, negative control. E, Immunoblots of expression of full-length BRAFF247L, wild-type BRAF, and BRAFV600E in Ba/F3. F, Dose-dependent survival assays of Ba/F3 cells expressing full-length BRAFF247L treated with dabrafenib and trametinib for 72 hours (mean percentage of cell survival; error bars, SD; N = 4, respectively). All P values calculated by t test; **, P < 0.01; ****, P < 0.0001.
Discussion
Over the past few years, cooperative NGS studies have identified more than 9,000 previously uncharacterized gene fusions (http://cgap.nci.nih.gov/Chromosomes/Mitelman); however, the functional consequence of the majority of these events remains to be determined. The functional interrogation of fusion genes is hampered by difficulties cloning such events using traditional molecular biology techniques, which involves obtaining tissues harboring the desired fusion gene or complex multi-step PCR as well as restriction digestion and ligation reactions that are inefficient.
Recent development of CRISPR/Cas9 genome editing technologies have enabled engineering of chromosomal rearrangements within endogenous loci (45–49). While these strategies offer great promise for characterizing fusion genes, inability to cover all types of chromosomal rearrangements (e.g., tandem duplication) and potential off-target complications render it less efficient for functional screening of the thousands of fusion genes identified from NGS studies. Publications (45–49) on modeling chromosomal rearrangements using CRISPR/Cas9 system reported a very low efficiency (1%–8%) for engineering events including EML4-ALK (8% and 4%), KIF5B-RET (1.6%), CD74-ROS1 (1.04%), and BCR-ABL (0.78%) into standard cell lines such as HEK293T and murine embryonic fibroblasts. The potential insertions and deletions resulting from nonhomologous end joining could easily cause mutations and/or frameshift in cases involving exon–exon fusions. Moreover, the incompatibility of mouse chromosomal orientations with human orthologs renders some fusion genes extremely difficult to engineer in animal models. On the other hand, we successfully made all 23 fusion genes in our study and achieved a >91% average success rate across all clones resulting from each recombineering step of our strategy. Resulting fusion genes are present in expression-ready vectors that can be immediately delivered to mouse and human cells for use in diverse downstream applications. DNA synthesis offers another route to fusion gene design, the cost of commercial gene synthesis ranges from $0.20 to $0.50/base, depending on gene sequence length plus additional service fees due to structural complexity thus eliminating this approach for most investigators wanting to build multiple constructs. For example, we were quoted $18,400 by a leading gene synthesis service company to provide the 23 fusion genes highlighted in our study, and this cost would not have covered additional work required for subcloning the synthesized DNA into expression vectors.
We applied our fusion gene construction strategy to a pilot study for construction of 3 bona fide fusion gene drivers and 20 uncharacterized fusion genes identified by TCGA. Our functional investigation revealed 5 fusion genes containing kinase domains encoded by MET, NTRK2, and BRAF. While the majority of the previously uncharacterized fusion genes appeared inactive in our test assays, it is likely that some of the remaining 15 fusion genes represent false negatives given that they were simply inactive in the Ba/F3 platform that is most sensitive to kinases (30). Future use of the described construction technology would perhaps be best applied for scaled production of cancer lineage-specific fusion genes, followed by functional screening in cancer context-specific model systems.
Our mechanistic studies on BRAF fusion genes corroborates other studies suggesting that loss of an N-terminal regulatory domain of wild-type BRAF contributes to the oncogenic activity of BRAF fusion events. Our fusion construction strategy can therefore be used to map inhibitory or activating protein sequences whose removal or ectopic insertion can activate oncogenes through gene fusion. Importantly, we also leveraged the versatility and modular design of our fusion gene construction strategy to discover a single missense mutation (F247L) within the N-terminus of BRAF that could potentially abolish its inhibitory function, leading to activation of the C-terminal kinase domain.
Another interesting observation made during this study relates to our finding that, despite being expressed from the same vector as AFAP1-NTRK2 and SQSTM1-NTRK2, we did not detect wild-type NTRK2 protein expression in our cell models even though its RNA transcript was expressed at levels similar to the NTRK2 fusion constructs. On the basis of this observation, we hypothesize that truncating NTRK2 through the process of fusion genesis results in increasing the overall stability of the fusion gene-encoded protein (and associated NTRK2 kinase domain) compared to wild-type protein, whose normal levels may be tightly regulated by proteasome or related machinery.
Our hypothesized means of NTRK2 activation, if proven, along with our confirmation that BRAF fusion activation occurs through loss of an N-terminal kinase inhibitory domain illustrates the notion that fusion genes can function through diverse mechanisms. Moreover, our observations related to differential subcellular localization and activity of MET fusion gene products further highlight the importance of understanding the biological functions of fusion gene proteins beyond their activity as a driver. While there have been numerous attempts to build computational algorithms to predict the functional impact of mutations on protein function, to our knowledge similar algorithms currently do not exist for predicting fusion proteins that are active in cancer. Such functional prediction algorithms for fusion proteins would greatly facilitate filtering of candidate fusion transcripts from the many identified from NGS datasets. Given the diverse mechanisms by which fusion genes activate, the systematic functional testing of numerous fusion gene events using construction pipelines such as the one presented here would ultimately inform development and enable refinement of functional prediction algorithms.
Disclosure of Potential Conflicts of Interest
G.B. Mills reports receiving commercial research grants from the Adelson Medical Research Foundation, AstraZeneca, Clinical Outcome Technologies, Komen Research Foundation, NanoString, Breast Cancer Research Foundation, Karus, Illumina, Takeda/Millenium Pharmaceuticals, has received speakers bureau honoraria from Symphogen, MedImmune, AstraZeneca, ISIS Pharmaceuticals, Lilly, Novartis, ImmunoMet, Allostery, Tarveda, and Pfizer; and has ownership interest (including patents) in Catena Pharmaceuticals, PTV Ventures, Spindletop Ventures, Myriad Genetics, ImmunoMet, and is a consultant/advisory board member for Adventist Health, AstraZeneca, Provista Diagnostics, Signalchem Lifesciences, Symphogen, Lilly, Novartis, Tarveda, Tau Therapeutics, Allostery, Catena Pharmaceuticals, Critical Outcome Technologies, ISIS Pharmaceuticals, ImmunoMet, Takeda/Millenium Pharmaceuticals, MedImmune, and Precision Medicine. No potential conflicts of interest were disclosed by the other authors.
Authors' Contributions
Conception and design: H. Lu, T. Dogruluk, K. Kong, N.J. Neill, R.G. Verhaak, G.B. Mills, R.S. Kucherlapati, K.L. Scott
Development of methodology: H. Lu, T. Dogruluk, K. Kong, Y.H. Tsang, N.J. Neill, K.L. Scott
Acquisition of data (provided animals, acquired and managed patients, provided facilities, etc.): H. Lu, N. Villafane, T. Dogruluk, C.L. Grzeskowiak, K. Kong, O. Zagorodna, Y.W. Kim, G.B. Mills, K.L. Scott
Analysis and interpretation of data (e.g., statistical analysis, biostatistics, computational analysis): H. Lu, N. Villafane, T. Dogruluk, O. Zagorodna, A. Pantazi, L. Yang, C.J. Creighton, G.B. Mills, P.J. Park, R.S. Kucherlapati, K.L. Scott
Writing, review, and/or revision of the manuscript: H. Lu, T. Dogruluk, C.L. Grzeskowiak, A. Pantazi, N.J. Neill, G.B. Mills, R.S. Kucherlapati, K.L. Scott
Administrative, technical, or material support (i.e., reporting or organizing data, constructing databases): H. Lu, R.G. Verhaak, K.L. Scott
Study supervision: K.L. Scott
Acknowledgments
We thank Zayed Institute for Personalized Cancer Therapy at The University of Texas MD Anderson Cancer Center for reagents and technical support.
Grant Support
This work was supported by the Cancer Prevention and Research Institute of Texas (CPRIT) RP120046 (to K.L. Scott), R21CA198320 from NIH (to K.L. Scott), U01CA168394 from NIH (to G.B. Mills and K.L. Scott), P30CA016672 from NIH (University of Texas MD Anderson Cancer Center,Houston, TX), and P30CA125123 from NIH (Baylor College of Medicine, Houston, TX). H. Lu and N. Villafane were supported by the CPRIT Pre-Doctoral Fellowship (RP140102) and Medical Resident Training Grant (RP14010), respectively.
The costs of publication of this article were defrayed in part by the payment of page charges. This article must therefore be hereby marked advertisement in accordance with 18 U.S.C. Section 1734 solely to indicate this fact.