Identification of liver cancer progenitors whose malignant progression depends on autocrine IL-6 signaling.

Hepatocellular carcinoma (HCC) is a slowly developing malignancy postulated to evolve from premalignant lesions in chronically damaged livers. However, it was never established that premalignant lesions actually contain tumor progenitors that give rise to cancer. Here, we describe isolation and characterization of HCC progenitor cells (HcPCs) from different mouse HCC models. Unlike fully malignant HCC, HcPCs give rise to cancer only when introduced into a liver undergoing chronic damage and compensatory proliferation. Although HcPCs exhibit a similar transcriptomic profile to bipotential hepatobiliary progenitors, the latter do not give rise to tumors. Cells resembling HcPCs reside within dysplastic lesions that appear several months before HCC nodules. Unlike early hepatocarcinogenesis, which depends on paracrine IL-6 production by inflammatory cells, due to upregulation of LIN28 expression, HcPCs had acquired autocrine IL-6 signaling that stimulates their in vivo growth and malignant progression. This may be a general mechanism that drives other IL-6-producing malignancies.

Deciphering cell-cell interactions and communication from gene expression.

Cell-cell interactions orchestrate organismal development, homeostasis and single-cell functions. When cells do not properly interact or improperly decode molecular messages, disease ensues. Thus, the identification and quantification of intercellular signalling pathways has become a common analysis performed across diverse disciplines. The expansion of protein-protein interaction databases and recent advances in RNA sequencing technologies have enabled routine analyses of intercellular signalling from gene expression measurements of bulk and single-cell data sets. In particular, ligand-receptor pairs can be used to infer intercellular communication from the coordinated expression of their cognate genes. In this Review, we highlight discoveries enabled by analyses of cell-cell interactions from transcriptomic data and review the methods and tools used in this context.

Whole-genome CRISPR screening identifies PI3K/AKT as a downstream component of the oncogenic GNAQ-focal adhesion kinase signaling circuitry.

G proteins and G protein-coupled receptors activate a diverse array of signal transduction pathways that promote cell growth and survival. Indeed, hot spot-activating mutations in GNAQ/GNA11, encoding Gαq proteins, are known to be driver oncogenes in uveal melanoma (UM), for which there are limited effective therapies currently available. Focal adhesion kinase (FAK) has been recently shown to be a central mediator of Gαq-driven signaling in UM, and as a result, is being explored clinically as a therapeutic target for UM, both alone and in combination therapies. Despite this, the repertoire of Gαq/FAK-regulated signaling mechanisms have not been fully elucidated. Here, we used a whole-genome CRISPR screen in GNAQ-mutant UM cells to identify mechanisms that, when overactivated, lead to reduced sensitivity to FAK inhibition. In this way, we found that the PI3K/AKT signaling pathway represented a major resistance driver. Our dissection of the underlying mechanisms revealed that Gαq promotes PI3K/AKT activation via a conserved signaling circuitry mediated by FAK. Further analysis demonstrated that FAK activates PI3K through the association and tyrosine phosphorylation of the p85 regulatory subunit of PI3K and that UM cells require PI3K/AKT signaling for survival. These findings establish a novel link between Gαq-driven signaling and the stimulation of PI3K as well as demonstrate aberrant activation of signaling networks underlying the growth and survival of UM and other Gαq-driven malignancies.

Autophagy gene haploinsufficiency drives chromosome instability, increases migration, and promotes early ovarian tumors.

Autophagy, particularly with BECN1, has paradoxically been highlighted as tumor promoting in Ras-driven cancers, but potentially tumor suppressing in breast and ovarian cancers. However, studying the specific role of BECN1 at the genetic level is complicated due to its genomic proximity to BRCA1 on both human (chromosome 17) and murine (chromosome 11) genomes. In human breast and ovarian cancers, the monoallelic deletion of these genes is often co-occurring. To investigate the potential tumor suppressor roles of two of the most commonly deleted autophagy genes in ovarian cancer, BECN1 and MAP1LC3B were knocked-down in atypical (BECN1+/+ and MAP1LC3B+/+) ovarian cancer cells. Ultra-performance liquid chromatography mass-spectrometry metabolomics revealed reduced levels of acetyl-CoA which corresponded with elevated levels of glycerophospholipids and sphingolipids. Migration rates of ovarian cancer cells were increased upon autophagy gene knockdown. Genomic instability was increased, resulting in copy-number alteration patterns which mimicked high grade serous ovarian cancer. We further investigated the causal role of Becn1 haploinsufficiency for oncogenesis in a MISIIR SV40 large T antigen driven spontaneous ovarian cancer mouse model. Tumors were evident earlier among the Becn1+/- mice, and this correlated with an increase in copy-number alterations per chromosome in the Becn1+/- tumors. The results support monoallelic loss of BECN1 as permissive for tumor initiation and potentiating for genomic instability in ovarian cancer.

COVID-19 Outcomes Among Patients With Cancer: Observations From the University of California Cancer Consortium COVID-19 Project Outcomes Registry.

BACKGROUND: The risks associated with severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) and its associated illness, coronavirus disease 2019 (COVID-19), among patients with a cancer diagnosis have not been fully characterized. This study leverages data from a multi-institutional cohort study, the University of California Cancer COVID Consortium, to evaluate outcomes associated with SARS-CoV-2 infection among patients with cancer. METHODS: Clinical data were collected from March to November 2020 and included patient demographics, cancer history and treatment, SARS-CoV-2 exposure and testing, and COVID-19 clinical management and outcomes. Multivariate ordinal logistic regression permitting unequal slopes was used to evaluate the impact of demographic, disease, and treatment factors on SARS-CoV-2 related hospitalization, intensive care unit (ICU) admission, and mortality. FINDINGS: Among all evaluated patients (n = 303), 147 (48%) were male, 118 (29%) were older adults (≥65 years old), and 104 (34%) were non-Hispanic white. A subset (n = 63, 21%) had hematologic malignancies and the remaining had solid tumors. Patients were hospitalized for acute care (n = 79, 26%), ICU-level care (n = 28, 9%), or died (n = 21, 7%) due to COVID-19. Patients with ≥2 comorbidities were more likely to require acute care (odds ratio [OR] 2.09 [95% confidence interval (CI), 1.23-3.55]). Cough was identified as a significant predictor of ICU hospitalization (OR 2.16 [95% CI, 1.03-4.57]). Importantly, mortality was associated with an active cancer diagnosis (OR 3.64 [95% CI, 1.40-9.5]) or advanced age (OR 3.86 [95% CI, 1.2-12.44]). INTERPRETATION: This study observed that patients with active cancer or advanced age are at an increased risk of death from COVID-19. These study observations can inform risk counseling related to COVID-19 for patients with a cancer diagnosis.

Accurate genome-wide genotyping from archival tissue to explore the contribution of common genetic variants to pre-cancer outcomes.

PURPOSE: The contribution of common genetic variants to pre-cancer progression is understudied due to long follow-up time, rarity of poor outcomes and lack of available germline DNA collection. Alternatively, DNA from diagnostic archival tissue is available, but its somatic nature, limited quantity and suboptimal quality would require an accurate cost-effective genome-wide germline genotyping methodology. EXPERIMENTAL DESIGN: Blood and tissue DNA from 10 individuals were used to benchmark the accuracy of Single Nucleotide Polymorphisms (SNP) genotypes, Polygenic Risk Scores (PRS) or HLA haplotypes using low-coverage whole-genome sequencing (lc-WGS) and genotype imputation. Tissue-derived PRS were further evaluated for 36 breast cancer patients (11.7 years median follow-up time) diagnosed with DCIS and used to model the risk of Breast Cancer Subsequent Events (BCSE). RESULTS: Tissue-derived germline DNA profiling resulted in accurate genotypes at common SNPs (blood correlation r2 > 0.94) and across 22 disease-related polygenic risk scores (PRS, mean correlation r = 0.93). Imputed Class I and II HLA haplotypes were 96.7% and 82.5% concordant with clinical-grade blood HLA haplotypes, respectively. In DCIS patients, tissue-derived PRS was significantly associated with BCSE (HR = 2, 95% CI 1.2-3.8). The top and bottom decile patients had an estimated 28% and 5% chance of BCSE at 10 years, respectively. CONCLUSIONS: Archival tissue DNA germline profiling using lc-WGS and imputation, represents a cost and resource-effective alternative in the retrospective design of long-term disease genetic studies. Initial results in breast cancer suggest that common risk variants contribute to pre-cancer progression.

Co-Localization of Gastrointestinal Stromal Tumors (GIST) and Peritoneal Mesothelioma: A Case Series.

PURPOSE: Gastrointestinal stromal tumor (GIST) is associated with increased risk of additional cancers. In this study, synchronous GIST, and peritoneal mesothelioma (PM) were characterized to evaluate the relationship between these two cancers. METHODS: A retrospective chart review was conducted for patients diagnosed with both GIST and PM between July 2010 and June 2021. Patient demographics, past tumor history, intraoperative reports, cross-sectional imaging, peritoneal cancer index (PCI) scoring, somatic next-generation sequencing (NGS) analysis, and histology were reviewed. RESULTS: Of 137 patients who underwent primary GIST resection from July 2010 to June 2021, 8 (5.8%) were found to have synchronous PM, and 4 patients (50%) had additional cancers and/or benign tumors. Five (62.5%) were male, and the median age at GIST diagnosis was 57 years (range: 45-76). Seventy-five percent of GISTs originated from the stomach. Of the eight patients, one patient had synchronous malignant mesothelioma (MM), and the remaining had well-differentiated papillary mesothelioma (WDPM), which were primarily located in the region of the primary GIST (89%). The median PCI score was 2 in the WDPM patients. NGS of GIST revealed oncogenic KIT exon 11 (62.5%), PDGFRA D842V (25%), or SDH (12.5%) mutations, while NGS of the MM revealed BAP1 and PBRM1 alterations. CONCLUSIONS: One in 17 GIST patients undergoing resection in this series have PM, which is significantly higher than expected if these two diseases were considered as independent events. Our results indicate that synchronous co-occurrence of GIST and PM is an underrecognized finding, suggesting a possible relationship that deserves further investigation.

Large and diffuse ductal carcinoma in situ: potentially lethal subtypes of "preinvasive" disease.

PURPOSE: Trials for DCIS have not explored whether outcomes for patients with large disease burden requiring mastectomy are comparable to those of patients with lumpectomy-amenable disease. We aim to identify whether patients with DCIS larger than 5 cm and diffuse-type DCIS differ in breast cancer mortality (BCM) from patients with disease less than 5 cm. METHODS: Patients diagnosed with DCIS in the SEER program were assessed to identify factors prognostic of breast-cancer-specific survival using competing risks regression. RESULTS: 44,849 patients met criteria for the cumulative incidence estimate. On competing risks cumulative incidence approximation, the 10-year estimate for BCM for each group was 1.3%, 1.3%, 2.3%, and 5.1%, respectively, and the difference among groups was significant (p = 0.017). On competing risks regression of patients with known covariates, both diffuse-type disease and disease larger than 5 cm (hazard ratio [HR] = 6.2 and 1.7, p = 0.013 and p = 0.042, respectively) were associated with increased risk of BCM. After matching, DCIS > 5 cm and diffuse disease were associated with increased BCM relative to disease < 5 cm (HR = 1.69, p = 0.04). Among patients undergoing mastectomy for disease larger than 5 cm or diffuse disease, the 10-year cumulative incidence for BCM was 0.5% among patients undergoing bilateral mastectomy and 2.4% for patients undergoing unilateral mastectomy. CONCLUSION: Patients with large and diffuse DCIS represent uncommon but poorly studied DCIS subgroups with worse prognoses than patients with disease smaller than 5 cm. Further studies are needed to elucidate the appropriate treatment for these patients.

Genomics Links Inflammation With Neurocognitive Impairment in Children Living With Human Immunodeficiency Virus Type-1.

BACKGROUND: We identified host single-nucleotide variants (SNVs) associated with neurocognitive impairment (NCI) in perinatally HIV-infected (PHIV) children. METHODS: Whole-exome sequencing (WES) was performed on 217 PHIV with cognitive score for age (CSA) < 70 and 247 CSA ≥ 70 (discovery cohort [DC]). SNVs identified in DC were evaluated in 2 validation cohorts (VC). Logistic regression was used to estimate adjusted odds ratios (ORs) for NCI. A human microglia NLRP3 inflammasome assay characterized the role of identified genes. RESULTS: Twenty-nine SNVs in 24 genes reaching P ≤ .002 and OR ≥ 1.5 comparing CSA < 70 to CSA ≥ 70 were identified in the DC, of which 3 SNVs were identified in VCs for further study. Combining the 3 cohorts, SNV in CCRL2 (rs3204849) was associated with decreased odds of NCI (P < .0001); RETREG1/FAM134B (rs61733811) and YWHAH (rs73884247) were associated with increased risk of NCI (P < .0001 and P < .001, respectively). Knockdown of CCRL2 led to decreased microglial release of IL-1ß following exposure to ssRNA40 while knockdown of RETREG1 and YWHAH resulted in increased IL-1ß release. CONCLUSIONS: Using WES and 2 VCs, and gene silencing of microglia we identified 3 genetic variants associated with NCI and inflammation in HIV-infected children.

Activation of hedgehog signaling associates with early disease progression in chronic lymphocytic leukemia.

Targeted sequencing of 103 leukemia-associated genes in leukemia cells from 841 treatment-naive patients with chronic lymphocytic leukemia (CLL) identified 89 (11%) patients as having CLL cells with mutations in genes encoding proteins that putatively are involved in hedgehog (Hh) signaling. Consistent with this finding, there was a significant association between the presence of these mutations and the expression of GLI1 (χ2 test, P < .0001), reflecting activation of the Hh pathway. However, we discovered that 38% of cases without identified mutations also were GLI1+ Patients with GLI1+ CLL cells had a shorter median treatment-free survival than patients with CLL cells lacking expression of GLI1 independent of IGHV mutation status. We found that GANT61, a small molecule that can inhibit GLI1, was highly cytotoxic for GLI1+ CLL cells relative to that of CLL cells without GLI1. Collectively, this study shows that a large proportion of patients have CLL cells with activated Hh signaling, which is associated with early disease progression and enhanced sensitivity to inhibition of GLI1.

Integrative genomic analysis of mouse and human hepatocellular carcinoma.

Cancer genomics has enabled the exhaustive molecular characterization of tumors and exposed hepatocellular carcinoma (HCC) as among the most complex cancers. This complexity is paralleled by dozens of mouse models that generate histologically similar tumors but have not been systematically validated at the molecular level. Accurate models of the molecular pathogenesis of HCC are essential for biomedical progress; therefore we compared genomic and transcriptomic profiles of four separate mouse models [MUP transgenic, TAK1-knockout, carcinogen-driven diethylnitrosamine (DEN), and Stelic Animal Model (STAM)] with those of 987 HCC patients with distinct etiologies. These four models differed substantially in their mutational load, mutational signatures, affected genes and pathways, and transcriptomes. STAM tumors were most molecularly similar to human HCC, with frequent mutations in Ctnnb1, similar pathway alterations, and high transcriptomic similarity to high-grade, proliferative human tumors with poor prognosis. In contrast, TAK1 tumors better reflected the mutational signature of human HCC and were transcriptionally similar to low-grade human tumors. DEN tumors were least similar to human disease and almost universally carried the Braf V637E mutation, which is rarely found in human HCC. Immune analysis revealed that strain-specific MHC-I genotype can influence the molecular makeup of murine tumors. Thus, different mouse models of HCC recapitulate distinct aspects of HCC biology, and their use should be adapted to specific questions based on the molecular features provided here.

Illuminating the Onco-GPCRome: Novel G protein-coupled receptor-driven oncocrine networks and targets for cancer immunotherapy.

G protein-coupled receptors (GPCRs) are the largest gene family of cell membrane-associated molecules mediating signal transmission, and their involvement in key physiological functions is well-established. The ability of GPCRs to regulate a vast array of fundamental biological processes, such as cardiovascular functions, immune responses, hormone and enzyme release from endocrine and exocrine glands, neurotransmission, and sensory perception (e.g. vision, odor, and taste), is largely due to the diversity of these receptors and the layers of their downstream signaling circuits. Dysregulated expression and aberrant functions of GPCRs have been linked to some of the most prevalent human diseases, which renders GPCRs one of the top targets for pharmaceutical drug development. However, the study of the role of GPCRs in tumor biology has only just begun to make headway. Recent studies have shown that GPCRs can contribute to the many facets of tumorigenesis, including proliferation, survival, angiogenesis, invasion, metastasis, therapy resistance, and immune evasion. Indeed, GPCRs are widely dysregulated in cancer and yet are underexploited in oncology. We present here a comprehensive analysis of GPCR gene expression, copy number variation, and mutational signatures in 33 cancer types. We also highlight the emerging role of GPCRs as part of oncocrine networks promoting tumor growth, dissemination, and immune evasion, and we stress the potential benefits of targeting GPCRs and their signaling circuits in the new era of precision medicine and cancer immunotherapies.

Deep mining heterogeneous networks of biomedical linked data to predict novel drug-target associations.

MOTIVATION: A heterogeneous network topology possessing abundant interactions between biomedical entities has yet to be utilized in similarity-based methods for predicting drug-target associations based on the array of varying features of drugs and their targets. Deep learning reveals features of vertices of a large network that can be adapted in accommodating the similarity-based solutions to provide a flexible method of drug-target prediction. RESULTS: We propose a similarity-based drug-target prediction method that enhances existing association discovery methods by using a topology-based similarity measure. DeepWalk, a deep learning method, is adopted in this study to calculate the similarities within Linked Tripartite Network (LTN), a heterogeneous network generated from biomedical linked datasets. This proposed method shows promising results for drug-target association prediction: 98.96% AUC ROC score with a 10-fold cross-validation and 99.25% AUC ROC score with a Monte Carlo cross-validation with LTN. By utilizing DeepWalk, we demonstrate that: (i) this method outperforms other existing topology-based similarity computation methods, (ii) the performance is better for tripartite than with bipartite networks and (iii) the measure of similarity using network topology outperforms the ones derived from chemical structure (drugs) or genomic sequence (targets). Our proposed methodology proves to be capable of providing a promising solution for drug-target prediction based on topological similarity with a heterogeneous network, and may be readily re-purposed and adapted in the existing of similarity-based methodologies. AVAILABILITY AND IMPLEMENTATION: The proposed method has been developed in JAVA and it is available, along with the data at the following URL: https://github.com/zongnansu1982/drug-target-prediction . CONTACT: nazong@ucsd.edu. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

ADAR1-mediated RNA-editing of 3'UTRs in breast cancer.

BACKGROUND: Whole transcriptome RNA variant analyses have shown that adenosine deaminases acting on RNA (ADAR) enzymes modify a large proportion of cellular RNAs, contributing to transcriptome diversity and cancer evolution. Despite the advances in the understanding of ADAR function in breast cancer, ADAR RNA editing functional consequences are not fully addressed. RESULTS: We characterized A to G(I) mRNA editing in 81 breast cell lines, showing increased editing at 3'UTR and exonic regions in breast cancer cells compared to immortalized non-malignant cell lines. In addition, tumors from the BRCA TCGA cohort show a 24% increase in editing over normal breast samples when looking at 571 well-characterized UTRs targeted by ADAR1. Basal-like subtype breast cancer patients with high level of ADAR1 mRNA expression shows a worse clinical outcome and increased editing in their 3'UTRs. Interestingly, editing was particularly increased in the 3'UTRs of ATM, GINS4 and POLH transcripts in tumors, which correlated with their mRNA expression. We confirmed the role of ADAR1 in this regulation using a shRNA in a breast cancer cell line (ZR-75-1). CONCLUSIONS: Altogether, these results revealed a significant association between the mRNA editing in genes related to cancer-relevant pathways and clinical outcomes, suggesting an important role of ADAR1 expression and function in breast cancer.

Dynamic epigenetic regulation of glioblastoma tumorigenicity through LSD1 modulation of MYC expression.

The available evidence suggests that the lethality of glioblastoma is driven by small subpopulations of cells that self-renew and exhibit tumorigenicity. It remains unclear whether tumorigenicity exists as a static property of a few cells or as a dynamically acquired property. We used tumor-sphere and xenograft formation as assays for tumorigenicity and examined subclones isolated from established and primary glioblastoma lines. Our results indicate that glioblastoma tumorigenicity is largely deterministic, yet the property can be acquired spontaneously at low frequencies. Further, these dynamic transitions are governed by epigenetic reprogramming through the lysine-specific demethylase 1 (LSD1). LSD depletion increases trimethylation of histone 3 lysine 4 at the avian myelocytomatosis viral oncogene homolog (MYC) locus, which elevates MYC expression. MYC, in turn, regulates oligodendrocyte lineage transcription factor 2 (OLIG2), SRY (sex determining region Y)-box 2 (SOX2), and POU class 3 homeobox 2 (POU3F2), a core set of transcription factors required for reprogramming glioblastoma cells into stem-like states. Our model suggests epigenetic regulation of key transcription factors governs transitions between tumorigenic states and provides a framework for glioblastoma therapeutic development.

Pan-cancer analysis reveals technical artifacts in TCGA germline variant calls.

BACKGROUND: Cancer research to date has largely focused on somatically acquired genetic aberrations. In contrast, the degree to which germline, or inherited, variation contributes to tumorigenesis remains unclear, possibly due to a lack of accessible germline variant data. Here we called germline variants on 9618 cases from The Cancer Genome Atlas (TCGA) database representing 31 cancer types. RESULTS: We identified batch effects affecting loss of function (LOF) variant calls that can be traced back to differences in the way the sequence data were generated both within and across cancer types. Overall, LOF indel calls were more sensitive to technical artifacts than LOF Single Nucleotide Variant (SNV) calls. In particular, whole genome amplification of DNA prior to sequencing led to an artificially increased burden of LOF indel calls, which confounded association analyses relating germline variants to tumor type despite stringent indel filtering strategies. The samples affected by these technical artifacts include all acute myeloid leukemia and practically all ovarian cancer samples. CONCLUSIONS: We demonstrate how technical artifacts induced by whole genome amplification of DNA can lead to false positive germline-tumor type associations and suggest TCGA whole genome amplified samples be used with caution. This study draws attention to the need to be sensitive to problems associated with a lack of uniformity in data generation in TCGA data.

Immunodeficiency and autoimmune enterocolopathy linked to NFAT5 haploinsufficiency.

The link between autoimmune diseases and primary immunodeficiency syndromes has been increasingly appreciated. Immunologic evaluation of a young man with autoimmune enterocolopathy and unexplained infections revealed evidence of immunodeficiency, including IgG subclass deficiency, impaired Ag-induced lymphocyte proliferation, reduced cytokine production by CD8(+) T lymphocytes, and decreased numbers of NK cells. Genetic evaluation identified haploinsufficiency of NFAT5, a transcription factor regulating immune cell function and cellular adaptation to hyperosmotic stress, as a possible cause of this syndrome. Inhibition or deletion of NFAT5 in normal human and murine cells recapitulated several of the immune deficits identified in the patient. These results provide evidence of a primary immunodeficiency disorder associated with organ-specific autoimmunity linked to NFAT5 deficiency.

9p21 DNA variants associated with coronary artery disease impair interferon-γ signalling response.

Genome-wide association studies have identified single nucleotide polymorphisms (SNPs) in the 9p21 gene desert associated with coronary artery disease (CAD) and type 2 diabetes. Despite evidence for a role of the associated interval in neighbouring gene regulation, the biological underpinnings of these genetic associations with CAD or type 2 diabetes have not yet been explained. Here we identify 33 enhancers in 9p21; the interval is the second densest gene desert for predicted enhancers and six times denser than the whole genome (P < 6.55 × 10(-33)). The CAD risk alleles of SNPs rs10811656 and rs10757278 are located in one of these enhancers and disrupt a binding site for STAT1. Lymphoblastoid cell lines homozygous for the CAD risk haplotype show no binding of STAT1, and in lymphoblastoid cell lines homozygous for the CAD non-risk haplotype, binding of STAT1 inhibits CDKN2BAS (also known as CDKN2B-AS1) expression, which is reversed by short interfering RNA knockdown of STAT1. Using a new, open-ended approach to detect long-distance interactions, we find that in human vascular endothelial cells the enhancer interval containing the CAD locus physically interacts with the CDKN2A/B locus, the MTAP gene and an interval downstream of IFNA21. In human vascular endothelial cells, interferon-γ activation strongly affects the structure of the chromatin and the transcriptional regulation in the 9p21 locus, including STAT1-binding, long-range enhancer interactions and altered expression of neighbouring genes. Our findings establish a link between CAD genetic susceptibility and the response to inflammatory signalling in a vascular cell type and thus demonstrate the utility of genome-wide association study findings in directing studies to novel genomic loci and biological processes important for disease aetiology.

FGFR1 and NTRK3 actionable alterations in "Wild-Type" gastrointestinal stromal tumors.

BACKGROUND: About 10-15% of adult, and most pediatric, gastrointestinal stromal tumors (GIST) lack mutations in KIT, PDGFRA, SDHx, or RAS pathway components (KRAS, BRAF, NF1). The identification of additional mutated genes in this rare subset of tumors can have important clinical benefit to identify altered biological pathways and select targeted therapies. METHODS: We performed comprehensive genomic profiling (CGP) for coding regions in more than 300 cancer-related genes of 186 GISTs to assess for their somatic alterations. RESULTS: We identified 24 GIST lacking alterations in the canonical KIT/PDGFRA/RAS pathways, including 12 without SDHx alterations. These 24 patients were mostly adults (96%). The tumors had a 46% rate of nodal metastases. These 24 GIST were more commonly mutated at 7 genes: ARID1B, ATR, FGFR1, LTK, SUFU, PARK2 and ZNF217. Two tumors harbored FGFR1 gene fusions (FGFR1-HOOK3, FGFR1-TACC1) and one harbored an ETV6-NTRK3 fusion that responded to TRK inhibition. In an independent sample set, we identified 5 GIST cases lacking alterations in the KIT/PDGFRA/SDHx/RAS pathways, including two additional cases with FGFR1-TACC1 and ETV6-NTRK3 fusions. CONCLUSIONS: Using patient demographics, tumor characteristics, and CGP, we show that GIST lacking alterations in canonical genes occur in younger patients, frequently metastasize to lymph nodes, and most contain deleterious genomic alterations, including gene fusions involving FGFR1 and NTRK3. If confirmed in larger series, routine testing for these translocations may be indicated for this subset of GIST. Moreover, these findings can be used to guide personalized treatments for patients with GIST. Trial registration NCT 02576431. Registered October 12, 2015.