A Gold Standard-Derived Modular Barcoding Approach to Cancer Transcriptomics.

A challenge with studying cancer transcriptomes is in distilling the wealth of information down into manageable portions of information. In this resource, we develop an approach that creates and assembles cancer type-specific gene expression modules into flexible barcodes, allowing for adaptation to a wide variety of uses. Specifically, we propose that modules derived organically from high-quality gold standards such as The Cancer Genome Atlas (TCGA) can accurately capture and describe functionally related genes that are relevant to specific cancer types. We show that such modules can: (1) uncover novel gene relationships and nominate new functional memberships, (2) improve and speed up analysis of smaller or lower-resolution datasets, (3) re-create and expand known cancer subtyping schemes, (4) act as a "decoder" to bridge seemingly disparate established gene signatures, and (5) efficiently apply single-cell RNA sequencing information to other datasets. Moreover, such modules can be used in conjunction with native spreadsheet program commands to create a powerful and rapid approach to hypothesis generation and testing that is readily accessible to non-bioinformaticians. Finally, we provide tools for users to create and interpret their own modules. Overall, the flexible modular nature of the proposed barcoding provides a user-friendly approach to rapidly decoding transcriptome-wide data for research or, potentially, clinical uses.

Comprehensive Immunogenomic Profiling of IDH1-/2-Altered Cholangiocarcinoma.

PURPOSE: Isocitrate dehydrogenase (IDH)1/2 genomic alterations (GA) occur in 20% of intrahepatic cholangiocarcinoma (iCCA); however, the immunogenomic landscape of IDH1-/2-mutated iCCA is largely unknown. METHODS: Comprehensive genomic profiling (CGP) was performed on 3,067 cases of advanced iCCA. Tumor mutational burden (TMB), PD-L1 expression (Dako 22C3), microsatellite instability (MSI), and genomic loss of heterozygosity (gLOH) as a surrogate marker for homologous recombination deficiency were examined. RNA sequencing of 73 patient samples was analyzed for differences in stromal/immune cell infiltration, immune marker expression, and T-cell inflammation. Tissue microarray arrays were subjected to multiplex immunohistochemistry and colocalization analysis in 100 surgical samples. Retrospective clinical data were collected for 501 patients with cholangiocarcinoma to examine median overall survival (mOS) in IDH1/2+ versus IDHwt. RESULTS: Of 3,067 iCCA cases subjected to CGP, 426 (14%) were IDH1+ and 125 (4%) were IDH2+. IDH1 GA included R132C (69%) and R132L/G/S/H/F (16%/7%/4%/3%/<1%). IDH2 GA occurred at R172 (94.4%) and R140 (6.6%). No significant difference was seen in median gLOH between IDH1+ versus IDHwt iCCA (P = .37), although patterns of comutations differed. MSI-High (P = .009), TMB ≥10 mut/Mb (P < .0001), and PD-L1 positivity were lower in IDH1/2+ versus IDHwt iCCA. Resting natural killer cell population, CD70, and programmed cell death 1 expression were significantly higher in non-IDH1-mutated cases, whereas V-set domain containing T-cell activation inhibitor 1 (B7-H4) expression was significantly higher in IDH1+. No significant difference in mOS was observed between IDH1/2+ versus IDHwt patients. CONCLUSION: Significant differences in GA and immune biomarkers are noted between IDH1/2+ and IDHwt iCCA. IDH1-/2-mutated tumors appear immunologically cold without gLOH. These immunogenomic data provide insight for precision targeting of iCCA with IDH alterations.

CRISPR screening identifies BET and mTOR inhibitor synergy in cholangiocarcinoma through serine glycine one carbon.

Patients with cholangiocarcinoma have poor clinical outcomes due to late diagnoses, poor prognoses, and limited treatment strategies. To identify drug combinations for this disease, we have conducted a genome-wide CRISPR screen anchored on the bromodomain and extraterminal domain (BET) PROTAC degrader ARV825, from which we identified anticancer synergy when combined with genetic ablation of members of the mTOR pathway. This combination effect was validated using multiple pharmacological BET and mTOR inhibitors, accompanied by increased levels of apoptosis and cell cycle arrest. In a xenograft model, combined BET degradation and mTOR inhibition induced tumor regression. Mechanistically, the 2 inhibitor classes converged on H3K27ac-marked epigenetic suppression of the serine glycine one carbon (SGOC) metabolism pathway, including the key enzymes PHGDH and PSAT1. Knockdown of PSAT1 was sufficient to replicate synergy with single-agent inhibition of either BET or mTOR. Our results tie together epigenetic regulation, metabolism, and apoptosis induction as key therapeutic targets for further exploration in this underserved disease.

Convergent MAPK pathway alterations mediate acquired resistance to FGFR inhibitors in FGFR2 fusion-positive cholangiocarcinoma.

BACKGROUND & AIMS: There is a knowledge gap in understanding mechanisms of resistance to fibroblast growth factor receptor (FGFR) inhibitors (FGFRi) and a need for novel therapeutic strategies to overcome it. We investigated mechanisms of acquired resistance to FGFRi in patients with FGFR2-fusion-positive cholangiocarcinoma (CCA). METHODS: A retrospective analysis of patients who received FGFRi therapy and underwent tumor and/or cell-free DNA analysis, before and after treatment, was performed. Longitudinal circulating tumor DNA samples from a cohort of patients in the phase I trial of futibatinib (NCT02052778) were assessed. FGFR2-BICC1 fusion cell lines were developed and secondary acquired resistance mutations in the mitogen-activated protein kinase (MAPK) pathway were introduced to assess their effect on sensitivity to FGFRi in vitro. RESULTS: On retrospective analysis of 17 patients with repeat sequencing following FGFRi treatment, new FGFR2 mutations were detected in 11 (64.7%) and new alterations in MAPK pathway genes in nine (52.9%) patients, with seven (41.2%) patients developing new alterations in both the FGFR2 and MAPK pathways. In serially collected plasma samples, a patient treated with an irreversible FGFRi tested positive for previously undetected BRAF V600E, NRAS Q61K, NRAS G12C, NRAS G13D and KRAS G12K mutations upon progression. Introduction of a FGFR2-BICC1 fusion into biliary tract cells in vitro sensitized the cells to FGFRi, while concomitant KRAS G12D or BRAF V600E conferred resistance. MEK inhibition was synergistic with FGFRi in vitro. In an in vivo animal model, the combination had antitumor activity in FGFR2 fusions but was not able to overcome KRAS-mediated FGFRi resistance. CONCLUSIONS: These findings suggest convergent genomic evolution in the MAPK pathway may be a potential mechanism of acquired resistance to FGFRi. CLINICAL TRIAL NUMBER: NCT02052778. IMPACT AND IMPLICATIONS: We evaluated tumors and plasma from patients who previously received inhibitors of fibroblast growth factor receptor (FGFR), an important receptor that plays a role in cancer cell growth, especially in tumors with abnormalities in this gene, such as FGFR fusions, where the FGFR gene is fused to another gene, leading to activation of cancer cell growth. We found that patients treated with FGFR inhibitors may develop mutations in other genes such as KRAS, and this can confer resistance to FGFR inhibitors. These findings have several implications for patients with FGFR2 fusion-positive tumors and provide mechanistic insight into emerging MAPK pathway alterations which may serve as a therapeutic vulnerability in the setting of acquired resistance to FGFRi.

USP1 promotes cholangiocarcinoma progression by deubiquitinating PARP1 to prevent its proteasomal degradation.

Despite its involvement in various cancers, the function of the deubiquitinase USP1 (ubiquitin-specific protease 1) is unexplored in cholangiocarcinoma (CCA). In this study, we provide evidence that USP1 promotes CCA progression through the stabilization of Poly (ADP-ribose) polymerase 1 (PARP1), consistent with the observation that both USP1 and PARP1 are upregulated in human CCA. Proteomics and ubiquitylome analysis of USP1-overexpressing CCA cells nominated PARP1 as a top USP1 substrate. Indeed, their direct interaction was validated by a series of immunofluorescence, co-immunoprecipitation (CO-IP), and GST pull-down assays, and their interaction regions were identified using deletion mutants. Mechanistically, USP1 removes the ubiquitin chain at K197 of PARP1 to prevent its proteasomal degradation, with the consequent PARP1 stabilization being necessary and sufficient to promote the growth and metastasis of CCA in vitro and in vivo. Additionally, we identified the acetyltransferase GCN5 as acetylating USP1 at K130, enhancing the affinity between USP1 and PARP1 and further increasing PARP1 protein stabilization. Finally, both USP1 and PARP1 are significantly associated with poor survival in CCA patients. These findings describe PARP1 as a novel deubiquitination target of USP1 and a potential therapeutic target for CCA.

Multiplatform Analysis of Intratumoral PTEN Heterogeneity in Melanoma.

Loss of protein expression of the tumor suppressor PTEN is associated with increased cancer aggressiveness, decreased tumor immune infiltration, and resistance to immune and targeted therapies in melanoma. We assessed a unique cohort of eight melanoma samples with focal loss of PTEN protein expression to understand the features and mechanisms of PTEN loss in this disease. We compared the PTEN-negative (PTEN[-]) areas to their adjacent PTEN-positive (PTEN[+]) areas using DNA sequencing, DNA methylation, RNA expression, digital spatial profiling, and immunohistochemical platforms. Variations or homozygous deletions of PTEN were identified in PTEN(-) areas that were not detected in the adjacent PTEN(+) areas in three cases (37.5%), but no clear genomic or DNA methylation basis for loss was identified in the remaining PTEN(-) samples. RNA expression data from two independent platforms identified a consistent increase in chromosome segregation gene expression in PTEN(-) versus adjacent PTEN(+) areas. Proteomic analysis showed a relative paucity of tumor-infiltrating lymphocytes in PTEN(-) versus adjacent PTEN(+) areas. The findings add to our understanding of potential molecular intratumoral heterogeneity in melanoma and the features associated with the loss of PTEN protein in this disease.

Microparticle-Delivered Cxcl9 Prolongs Braf Inhibitor Efficacy in Melanoma.

Patients with BRAF-mutant melanoma show substantial responses to combined BRAF and MEK inhibition, but most relapse within 2 years. A major reservoir for drug resistance is minimal residual disease (MRD), comprised of drug-tolerant tumor cells laying in a dormant state. Towards exploiting potential therapeutic vulnerabilities of MRD, we established a genetically engineered mouse model of BrafV600E-driven melanoma MRD wherein genetic BrafV600E extinction leads to strong but incomplete tumor regression. Transcriptional time-course analysis after BrafV600E extinction revealed that after an initial surge of immune activation, tumors later became immunologically "cold" after MRD establishment. Computational analysis identified candidate T-cell recruiting chemokines as strongly upregulated initially and steeply decreasing as the immune response faded. Therefore, we hypothesized that sustaining chemokine signaling could impair MRD maintenance through increased recruitment of effector T cells. We found that intratumoral administration of recombinant Cxcl9 (rCxcl9), either naked or loaded in microparticles, significantly impaired MRD relapse in BRAF-inhibited tumors, including several complete pathologic responses after microparticle-delivered rCxcl9 combined with BRAF and MEK inhibition. Our experiments constitute proof of concept that chemokine-based microparticle delivery systems are a potential strategy to forestall tumor relapse and thus improve the clinical success of first-line treatment methods.

Targeting RAS Mutant Colorectal Cancer with Dual Inhibition of MEK and CDK4/6.

KRAS and NRAS mutations occur in 45% of colorectal cancers, with combined MAPK pathway and CDK4/6 inhibition identified as a potential therapeutic strategy. In the current study, this combinatorial treatment approach was evaluated in a co-clinical trial in patient-derived xenografts (PDX), and safety was established in a clinical trial of binimetinib and palbociclib in patients with metastatic colorectal cancer with RAS mutations. Across 18 PDX models undergoing dual inhibition of MEK and CDK4/6, 60% of tumors regressed, meeting the co-clinical trial primary endpoint. Prolonged duration of response occurred predominantly in TP53 wild-type models. Clinical evaluation of binimetinib and palbociclib in a safety lead-in confirmed safety and provided preliminary evidence of activity. Prolonged treatment in PDX models resulted in feedback activation of receptor tyrosine kinases and acquired resistance, which was reversed with a SHP2 inhibitor. These results highlight the clinical potential of this combination in colorectal cancer, along with the utility of PDX-based co-clinical trial platforms for drug development. SIGNIFICANCE: This co-clinical trial of combined MEK-CDK4/6 inhibition in RAS mutant colorectal cancer demonstrates therapeutic efficacy in patient-derived xenografts and safety in patients, identifies biomarkers of response, and uncovers targetable mechanisms of resistance.

Multi-modal molecular programs regulate melanoma cell state.

Melanoma cells display distinct intrinsic phenotypic states. Here, we seek to characterize the molecular regulation of these states using multi-omic analyses of whole exome, transcriptome, microRNA, long non-coding RNA and DNA methylation data together with reverse-phase protein array data on a panel of 68 highly annotated early passage melanoma cell lines. We demonstrate that clearly defined cancer cell intrinsic transcriptomic programs are maintained in melanoma cells ex vivo and remain highly conserved within melanoma tumors, are associated with distinct immune features within tumors, and differentially correlate with checkpoint inhibitor and adoptive T cell therapy efficacy. Through integrative analyses we demonstrate highly complex multi-omic regulation of melanoma cell intrinsic programs that provide key insights into the molecular maintenance of phenotypic states. These findings have implications for cancer biology and the identification of new therapeutic strategies. Further, these deeply characterized cell lines will serve as an invaluable resource for future research in the field.

IDH1 Inhibition Reawakens the Immune Response against Cholangiocarcinoma.

Isocitrate dehydrogenase 1 mutations (mIDH1) are common in cholangiocarcinoma, but their exact mechanisms in cholangiocarcinoma initiation and maintenance are unclear. In this issue of Cancer Discovery, Wu and colleagues identify immune suppression via TET2 inactivation as the primary means by which mIDH1 maintains cholangiocarcinoma survival, leading to an efficacious new combination of mIDH1 inhibitors and immune checkpoint blockade targeting regulatory T cells. See related article by Wu et al., p. 812 (9).

Limitations and opportunities of technologies for the analysis of cell-free DNA in cancer diagnostics.

Cell-free DNA (cfDNA) in the circulating blood plasma of patients with cancer contains tumour-derived DNA sequences that can serve as biomarkers for guiding therapy, for the monitoring of drug resistance, and for the early detection of cancers. However, the analysis of cfDNA for clinical diagnostic applications remains challenging because of the low concentrations of cfDNA, and because cfDNA is fragmented into short lengths and is susceptible to chemical damage. Barcodes of unique molecular identifiers have been implemented to overcome the intrinsic errors of next-generation sequencing, which is the prevailing method for highly multiplexed cfDNA analysis. However, a number of methodological and pre-analytical factors limit the clinical sensitivity of the cfDNA-based detection of cancers from liquid biopsies. In this Review, we describe the state-of-the-art technologies for cfDNA analysis, with emphasis on multiplexing strategies, and discuss outstanding biological and technical challenges that, if addressed, would substantially improve cancer diagnostics and patient care.

Monitoring of Dynamic Changes and Clonal Evolution in Circulating Tumor DNA From Patients With IDH-Mutated Cholangiocarcinoma Treated With Isocitrate Dehydrogenase Inhibitors.

PURPOSE: IDH mutations occur in about 30% of patients with cholangiocarcinoma. Analysis of mutations in circulating tumor DNA (ctDNA) can be performed by droplet digital polymerase chain reaction (ddPCR). The analysis of ctDNA is a feasible approach to detect IDH mutations. METHODS: We isolated ctDNA from the blood of patients with IDH-mutated advanced cholangiocarcinoma collected at baseline, on therapy, and at progression to isocitrate dehydrogenase (IDH) inhibitors. RESULTS: Of 31 patients with IDH1R132 (n = 26) or IDH2R172 mutations (n = 5) in the tumor, IDH mutations were detected in 84% of ctDNA samples analyzed by ddPCR and in 83% of ctDNA samples analyzed by next-generation sequencing (NGS). Patients with a low variant allele frequency of ctDNA detected by NGS at baseline had a longer median time to treatment failure compared to patients with high variant allele frequency of ctDNA (3.6 v 1.5 months; P = .008). Patients with a decrease in IDH-mutated ctDNA on therapy by ddPCR compared with no change/increase had a trend to a longer median survival (P = .07). Most frequent emergent alterations in ctDNA by NGS at progression were ARID1A (n = 3) and TP53 mutations (n = 3). CONCLUSION: Detection of IDH mutations in ctDNA in patients with advanced cholangiocarcinoma is feasible, and dynamic changes in ctDNA can correspond with the clinical course and clonal evolution.

The immunogenomic landscape of resected intrahepatic cholangiocarcinoma.

BACKGROUND AND AIMS: Cholangiocarcinoma (CCA) is a deadly and highly therapy-refractory cancer of the bile ducts, with early results from immune checkpoint blockade trials showing limited responses. Whereas recent molecular assessments have made bulk characterizations of immune profiles and their genomic correlates, spatial assessments may reveal actionable insights. APPROACH AND RESULTS: Here, we have integrated immune checkpoint-directed immunohistochemistry with next-generation sequencing of resected intrahepatic CCA samples from 96 patients. We found that both T-cell and immune checkpoint markers are enriched at the tumor margins compared to the tumor center. Using two approaches, we identify high programmed cell death protein 1 or lymphocyte-activation gene 3 and low CD3/CD4/inducible T-cell costimulator specifically in the tumor center as associated with poor survival. Moreover, loss-of-function BRCA1-associated protein-1 mutations are associated with and cause elevated expression of the immunosuppressive checkpoint marker, B7 homolog 4. CONCLUSIONS: This study provides a foundation on which to rationally improve and tailor immunotherapy approaches for this difficult-to-treat disease.

Neural Crest-Like Stem Cell Transcriptome Analysis Identifies LPAR1 in Melanoma Progression and Therapy Resistance.

Metastatic melanoma is challenging to clinically address. Although standard-of-care targeted therapy has high response rates in patients with BRAF-mutant melanoma, therapy relapse occurs in most cases. Intrinsically resistant melanoma cells drive therapy resistance and display molecular and biologic properties akin to neural crest-like stem cells (NCLSC) including high invasiveness, plasticity, and self-renewal capacity. The shared transcriptional programs and vulnerabilities between NCLSCs and cancer cells remains poorly understood. Here, we identify a developmental LPAR1-axis critical for NCLSC viability and melanoma cell survival. LPAR1 activity increased during progression and following acquisition of therapeutic resistance. Notably, genetic inhibition of LPAR1 potentiated BRAFi ± MEKi efficacy and ablated melanoma migration and invasion. Our data define LPAR1 as a new therapeutic target in melanoma and highlights the promise of dissecting stem cell-like pathways hijacked by tumor cells. SIGNIFICANCE: This study identifies an LPAR1-axis critical for melanoma invasion and intrinsic/acquired therapy resistance.

Cholangiocarcinoma Risk Factors Open the Floodgates for Gut Microbes and Immunosuppressive Myeloid Cells.

Primary sclerosing cholangitis and colitis are known predisposing factors for bile duct cancer, but their exact pro-oncogenic mechanisms are unclear. In this issue of Cancer Discovery, Zhang and colleagues identify intestinal barrier impairment as a key mechanism, resulting in gut microbes spilling into the portal vein, in turn recruiting immunosuppressive myeloid-derived suppressor cells and promoting cholangiocarcinoma.See related article by Zhang et al. p. 1248.

High sensitivity sanger sequencing detection of BRAF mutations in metastatic melanoma FFPE tissue specimens.

Mutations in the BRAF gene at or near the p. V600 locus are informative for therapy selection, but current methods for analyzing FFPE tissue DNA generally have a limit of detection of 5% variant allele frequency (VAF), or are limited to the single variant (V600E). These can result in false negatives for samples with low VAFs due to low tumor content or subclonal heterogeneity, or harbor non-V600 mutations. Here, we show that Sanger sequencing using the NuProbe VarTrace BRAF assay, based on the Blocker Displacement Amplification (BDA) technology, is capable of detecting BRAF V600 mutations down to 0.20% VAF from FFPE lymph node tissue samples. Comparison experiments on adjacent tissue sections using BDA Sanger, immunohistochemistry (IHC), digital droplet PCR (ddPCR), and NGS showed 100% concordance among all 4 methods for samples with BRAF mutations at ≥ 1% VAF, though ddPCR did not distinguish the V600K mutation from the V600E mutation. BDA Sanger, ddPCR, and NGS (with orthogonal confirmation) were also pairwise concordant for lower VAF mutations down to 0.26% VAF, but IHC produced a false negative. Thus, we have shown that Sanger sequencing can be effective for rapid detection and quantitation of multiple low VAF BRAF mutations from FFPE samples. BDA Sanger method also enabled detection and quantitation of less frequent, potentially actionable non-V600 mutations as demonstrated by synthetic samples.

Genomic profiling reveals high frequency of DNA repair genetic aberrations in gallbladder cancer.

DNA repair gene aberrations (GAs) occur in several cancers, may be prognostic and are actionable. We investigated the frequency of DNA repair GAs in gallbladder cancer (GBC), association with tumor mutational burden (TMB), microsatellite instability (MSI), programmed cell death protein 1 (PD-1), and its ligand (PD-L1) expression. Comprehensive genomic profiling (CGP) of 760 GBC was performed. We investigated GAs in 19 DNA repair genes including direct DNA repair genes (ATM, ATR, BRCA1, BRCA2, FANCA, FANCD2, MLH1, MSH2, MSH6, PALB2, POLD1, POLE, PRKDC, and RAD50) and caretaker genes (BAP1, CDK12, MLL3, TP53, and BLM) and classified patients into 3 groups based on TMB level: low (< 5.5 mutations/Mb), intermediate (5.5-19.5 mutations/Mb), and high (≥ 19.5 mutations/Mb). We assessed MSI status and PD-1 & PD-L1 expression. 658 (86.6%) had at least 1 actionable GA. Direct DNA repair gene GAs were identified in 109 patients (14.2%), while 476 (62.6%) had GAs in caretaker genes. Both direct and caretaker DNA repair GAs were significantly associated with high TMB (P = 0.0005 and 0.0001, respectively). Tumor PD-L1 expression was positive in 119 (15.6%), with 17 (2.2%) being moderate or high. DNA repair GAs are relatively frequent in GBC and associated with coexisting actionable mutations and a high TMB.

Intrahepatic Cholangiocarcinoma: Genomic Heterogeneity Between Eastern and Western Patients.

PURPOSE: Intrahepatic cholangiocarcinoma (IHCCA), a global health problem, is increasing in incidence and has differing etiologies worldwide. Next-generation sequencing (NGS) is rapidly being incorporated into the clinical management of biliary cancers. IHCCA is enriched with actionable mutations, and there are several promising targeted therapies under development. NGS data from Asia, where IHCCA is most prevalent, are limited. METHODS: Comprehensive genomic profiling of formalin-fixed paraffin-embedded tumor tissue from 164 Asian and 283 Western patients with IHCCA was performed using NGS. We measured the distribution of DNA repair genetic aberrations (GAs) in IHCCA, along with actionable mutations. Also, we evaluated the association between DNA repair GAs and tumor mutation burden (TMB). Based on the TMB status, patients were distinguished into 3 levels: low (< 6 mut/Mb), intermediate (6-10 mut/Mb), and high (TMB-H; ≥ 10 mut/Mb). RESULTS: Seventy-two percent of Asian patients had ≥ 1 actionable GA, with a significantly higher frequency in KMT2C , BRCA1/2, and DDR2 compared with Western patients (P = .02, .003, and .003, respectively); 60.9% of Western patients had ≥ 1 actionable GA and higher frequency of CDKN2A/B and IDH1/2 GAs (P = .0004 and < .001, respectively). GAs in nuclear factor kappa B pathway regulators and DNA repair genes occurred more frequently in Asian patients (P = .006 and .001, respectively). There was a higher frequency of TMB-H in Asian compared with the Western cohort (12.2% v 5.9%; P = .07). CONCLUSION: A higher burden of DNA repair mutations and frequency of patients with TMB-H in the Asian IHCCA cohort compared with the Western patients suggests a potential role for DNA repair and immune checkpoint inhibitors in the Asian population. Future clinical trials should account for this genetic heterogeneity.

Same Name, Different Game: EGFR Drives Intrinsic KRASG12C Inhibitor Resistance in Colorectal Cancer.

The discovery of covalent inhibitors of KRASG12C has led to promising clinical results in lung cancer, but disappointing response rates in colon cancer. In this issue of Cancer Discovery, Misale and colleagues identify high endogenous EGFR activity as the underlying mechanism of intrinsic resistance, which can be overcome by anti-EGFR antibody coadministration.See related article by Amodio et al., p. 1129.

Generation of An Endogenous FGFR2-BICC1 Gene Fusion/58 Megabase Inversion Using Single-Plasmid CRISPR/Cas9 Editing in Biliary Cells.

Fibroblast growth factor receptor 2 (FGFR2) gene fusions are bona fide oncogenic drivers in 10-15% of intrahepatic cholangiocarcinoma (CCA), yet currently there are no cell lines publically available to study endogenous FGFR2 gene fusions. The ability of clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 to generate large yet precise chromosomal rearrangements has presented the possibility of engineering endogenous gene fusions for downstream studies. In this technical report, we describe the generation of an endogenous FGFR2-Bicaudal family RNA binding protein 1 (BICC1) fusion in multiple independent cholangiocarcinoma and immortalized liver cell lines using CRISPR. BICC1 is the most common FGFR2 fusion partner in CCA, and the fusion arises as a consequence of a 58-megabase-sized inversion on chromosome 10. We replicated this inversion to generate a fusion product that is identical to that seen in many human CCA. Our results demonstrate the feasibility of generating large megabase-scale inversions that faithfully reproduce human cancer aberrations.