Integrative multi-omics characterization reveals sex differences in glioblastoma.

BACKGROUND: Glioblastoma (GBM) is the most common and lethal primary brain tumor in adults, with limited treatment modalities and poor prognosis. Recent studies have highlighted the importance of considering sex differences in cancer incidence, prognosis, molecular disparities, and treatment outcomes across various tumor types, including colorectal adenocarcinoma, lung adenocarcinoma, and GBM. METHODS: We performed comprehensive analyses of large-scale multi-omics data (genomic, transcriptomic, and proteomic data) from TCGA, GLASS, and CPTAC to investigate the genetic and molecular determinants that contribute to the unique clinical properties of male and female GBM patients. RESULTS: Our results revealed several key differences, including enrichments of MGMT promoter methylation, which correlated with increased overall and post-recurrence survival and improved response to chemotherapy in female patients. Moreover, female GBM exhibited a higher degree of genomic instability, including aneuploidy and tumor mutational burden. Integrative proteomic and phosphor-proteomic characterization uncovered sex-specific protein abundance and phosphorylation activities, including EGFR activation in males and SPP1 hyperphosphorylation in female patients. Lastly, the identified sex-specific biomarkers demonstrated prognostic significance, suggesting their potential as therapeutic targets. CONCLUSIONS: Collectively, our study provides unprecedented insights into the fundamental modulators of tumor progression and clinical outcomes between male and female GBM patients and facilitates sex-specific treatment interventions. Highlights Female GBM patients were characterized by increased MGMT promoter methylation and favorable clinical outcomes compared to male patients. Female GBMs exhibited higher levels of genomic instability, including aneuploidy and TMB. Each sex-specific GBM is characterized by unique pathway dysregulations and molecular subtypes. EGFR activation is prevalent in male patients, while female patients are marked by SPP1 hyperphosphorylation.

Glioblastoma (GBM) is the most common and lethal primary brain tumor in adults with limited treatment modalities and dismal prognosis. A thorough understanding of sex differences could facilitate personalized therapeutic strategies in GBM. In this study, we conducted a comprehensive multi-omics analysis from TCGA, CPTAC, and GLASS studies, revealing distinct molecular and clinical disparities between male and female GBM patients. We discovered that female GBM patients exhibited enrichments of MGMT promoter methylation and high genomic instability, including aneuploidy and TMB. While male GBMs were characterized by activation of EGFR protein and phosphorylation activities, female GBM patients demonstrated hyperphosphorylation of SPP1. Notably, these proteins demonstrated prognostic significance, highlighting their potential as therapeutic targets. Our findings provide unprecedented insights into the fundamental modulators of tumor progression and clinical outcomes in male and female GBM patients, offering valuable opportunities for sex-specific treatment interventions.

Integrated proteogenomic characterization of glioblastoma evolution.

The evolutionary trajectory of glioblastoma (GBM) is a multifaceted biological process that extends beyond genetic alterations alone. Here, we perform an integrative proteogenomic analysis of 123 longitudinal glioblastoma pairs and identify a highly proliferative cellular state at diagnosis and replacement by activation of neuronal transition and synaptogenic pathways in recurrent tumors. Proteomic and phosphoproteomic analyses reveal that the molecular transition to neuronal state at recurrence is marked by post-translational activation of the wingless-related integration site (WNT)/ planar cell polarity (PCP) signaling pathway and BRAF protein kinase. Consistently, multi-omic analysis of patient-derived xenograft (PDX) models mirror similar patterns of evolutionary trajectory. Inhibition of B-raf proto-oncogene (BRAF) kinase impairs both neuronal transition and migration capability of recurrent tumor cells, phenotypic hallmarks of post-therapy progression. Combinatorial treatment of temozolomide (TMZ) with BRAF inhibitor, vemurafenib, significantly extends the survival of PDX models. This study provides comprehensive insights into the biological mechanisms of glioblastoma evolution and treatment resistance, highlighting promising therapeutic strategies for clinical intervention.

Pharmacogenomic profiling reveals molecular features of chemotherapy resistance in IDH wild-type primary glioblastoma.

BACKGROUND: Although temozolomide (TMZ) has been used as a standard adjuvant chemotherapeutic agent for primary glioblastoma (GBM), treating isocitrate dehydrogenase wild-type (IDH-wt) cases remains challenging due to intrinsic and acquired drug resistance. Therefore, elucidation of the molecular mechanisms of TMZ resistance is critical for its precision application. METHODS: We stratified 69 primary IDH-wt GBM patients into TMZ-resistant (n = 29) and sensitive (n = 40) groups, using TMZ screening of the corresponding patient-derived glioma stem-like cells (GSCs). Genomic and transcriptomic features were then examined to identify TMZ-associated molecular alterations. Subsequently, we developed a machine learning (ML) model to predict TMZ response from combined signatures. Moreover, TMZ response in multisector samples (52 tumor sectors from 18 cases) was evaluated to validate findings and investigate the impact of intra-tumoral heterogeneity on TMZ efficacy. RESULTS: In vitro TMZ sensitivity of patient-derived GSCs classified patients into groups with different survival outcomes (P = 1.12e-4 for progression-free survival (PFS) and 3.63e-4 for overall survival (OS)). Moreover, we found that elevated gene expression of EGR4, PAPPA, LRRC3, and ANXA3 was associated to intrinsic TMZ resistance. In addition, other features such as 5-aminolevulinic acid negative, mesenchymal/proneural expression subtypes, and hypermutation phenomena were prone to promote TMZ resistance. In contrast, concurrent copy-number-alteration in PTEN, EGFR, and CDKN2A/B was more frequent in TMZ-sensitive samples (Fisher's exact P = 0.0102), subsequently consolidated by multi-sector sequencing analyses. Integrating all features, we trained a ML tool to segregate TMZ-resistant and sensitive groups. Notably, our method segregated IDH-wt GBM patients from The Cancer Genome Atlas (TCGA) into two groups with divergent survival outcomes (P = 4.58e-4 for PFS and 3.66e-4 for OS). Furthermore, we showed a highly heterogeneous TMZ-response pattern within each GBM patient using in vitro TMZ screening and genomic characterization of multisector GSCs. Lastly, the prediction model that evaluates the TMZ efficacy for primary IDH-wt GBMs was developed into a webserver for public usage ( http://www.wang-lab-hkust.com:3838/TMZEP ). CONCLUSIONS: We identified molecular characteristics associated to TMZ sensitivity, and illustrate the potential clinical value of a ML model trained from pharmacogenomic profiling of patient-derived GSC against IDH-wt GBMs.

IGFBP5 is an ROR1 ligand promoting glioblastoma invasion via ROR1/HER2-CREB signaling axis.

Diffuse infiltration is the main reason for therapeutic resistance and recurrence in glioblastoma (GBM). However, potential targeted therapies for GBM stem-like cell (GSC) which is responsible for GBM invasion are limited. Herein, we report Insulin-like Growth Factor-Binding Protein 5 (IGFBP5) is a ligand for Receptor tyrosine kinase like Orphan Receptor 1 (ROR1), as a promising target for GSC invasion. Using a GSC-derived brain tumor model, GSCs were characterized into invasive or non-invasive subtypes, and RNA sequencing analysis revealed that IGFBP5 was differentially expressed between these two subtypes. GSC invasion capacity was inhibited by IGFBP5 knockdown and enhanced by IGFBP5 overexpression both in vitro and in vivo, particularly in a patient-derived xenograft model. IGFBP5 binds to ROR1 and facilitates ROR1/HER2 heterodimer formation, followed by inducing CREB-mediated ETV5 and FBXW9 expression, thereby promoting GSC invasion and tumorigenesis. Importantly, using a tumor-specific targeting and penetrating nanocapsule-mediated delivery of CRISPR/Cas9-based IGFBP5 gene editing significantly suppressed GSC invasion and downstream gene expression, and prolonged the survival of orthotopic tumor-bearing mice. Collectively, our data reveal that IGFBP5-ROR1/HER2-CREB signaling axis as a potential GBM therapeutic target.

Mutation-specific non-canonical pathway of PTEN as a distinct therapeutic target for glioblastoma.

PTEN is one of the most frequently altered tumor suppressor genes in malignant tumors. The dominant-negative effect of PTEN alteration suggests that the aberrant function of PTEN mutation might be more disastrous than deletion, the most frequent genomic event in glioblastoma (GBM). This study aimed to understand the functional properties of various PTEN missense mutations and to investigate their clinical relevance. The genomic landscape of PTEN alteration was analyzed using the Samsung Medical Center GBM cohort and validated via The Cancer Genome Atlas dataset. Several hotspot mutations were identified, and their subcellular distributions and phenotypes were evaluated. We established a library of cancer cell lines that overexpress these mutant proteins using the U87MG and patient-derived cell models lacking functional PTEN. PTEN mutations were categorized into two major subsets: missense mutations in the phosphatase domain and truncal mutations in the C2 domain. We determined the subcellular compartmentalization of four mutant proteins (H93Y, C124S, R130Q, and R173C) from the former group and found that they had distinct localizations; those associated with invasive phenotypes ('edge mutations') localized to the cell periphery, while the R173C mutant localized to the nucleus. Invasive phenotypes derived from edge substitutions were unaffected by an anti-PI3K/Akt agent but were disrupted by microtubule inhibitors. PTEN mutations exhibit distinct functional properties regarding their subcellular localization. Further, some missense mutations ('edge mutations') in the phosphatase domain caused enhanced invasiveness associated with dysfunctional cytoskeletal assembly, thus suggesting it to be a potent therapeutic target.

Modulation of Nogo receptor 1 expression orchestrates myelin-associated infiltration of glioblastoma.

As the clinical failure of glioblastoma treatment is attributed by multiple components, including myelin-associated infiltration, assessment of the molecular mechanisms underlying such process and identification of the infiltrating cells have been the primary objectives in glioblastoma research. Here, we adopted radiogenomic analysis to screen for functionally relevant genes that orchestrate the process of glioma cell infiltration through myelin and promote glioblastoma aggressiveness. The receptor of the Nogo ligand (NgR1) was selected as the top candidate through Differentially Expressed Genes (DEG) and Gene Ontology (GO) enrichment analysis. Gain and loss of function studies on NgR1 elucidated its underlying molecular importance in suppressing myelin-associated infiltration in vitro and in vivo. The migratory ability of glioblastoma cells on myelin is reversibly modulated by NgR1 during differentiation and dedifferentiation process through deubiquitinating activity of USP1, which inhibits the degradation of ID1 to downregulate NgR1 expression. Furthermore, pimozide, a well-known antipsychotic drug, upregulates NgR1 by post-translational targeting of USP1, which sensitizes glioma stem cells to myelin inhibition and suppresses myelin-associated infiltration in vivo. In primary human glioblastoma, downregulation of NgR1 expression is associated with highly infiltrative characteristics and poor survival. Together, our findings reveal that loss of NgR1 drives myelin-associated infiltration of glioblastoma and suggest that novel therapeutic strategies aimed at reactivating expression of NgR1 will improve the clinical outcome of glioblastoma patients.

Therapeutic Efficacy of GC1118, a Novel Anti-EGFR Antibody, against Glioblastoma with High EGFR Amplification in Patient-Derived Xenografts.

We aimed to evaluate the preclinical efficacy of GC1118, a novel anti-epidermal growth factor receptor (EGFR) monoclonal antibody (mAb), against glioblastoma (GBM) tumors using patient-derived xenograft (PDX) models. A total of 15 distinct GBM PDX models were used to evaluate the therapeutic efficacy of GC1118. Genomic data derived from PDX models were analyzed to identify potential biomarkers associated with the anti-tumor efficacy of GC1118. A patient-derived cell-based high-throughput drug screening assay was performed to further validate the efficacy of GC1118. Compared to cetuximab, GC1118 exerted comparable growth inhibitory effects on the GBM tumors in the PDX models. We confirmed that GC1118 accumulated within the tumor by crossing the blood-brain barrier in in vivo specimens and observed the survival benefit in GC1118-treated intracranial models. Genomic analysis revealed high EGFR amplification as a potent biomarker for predicting the therapeutic efficacy of GC1118 in GBM tumors. In summary, GC1118 exerted a potent anti-tumor effect on GBM tumors in PDX models, and its therapeutic efficacy was especially pronounced in the tumors with high EGFR amplification. Our study supports the importance of patient stratification based on EGFR copy number variation in clinical trials for GBM. The superiority of GC1118 over other EGFR mAbs in GBM tumors should be assessed in future studies.

Ethnic delineation of primary glioblastoma genome.

Glioblastoma (GBM) is the most malignant primary brain tumor in adults with substantial genomic alterations. The median survival is approximately 14.6 months, despite aggressive therapeutic intervention, which comprised of surgical resection, radiotherapy, and chemotherapy. Recent studies on cancer genomic have revealed crucial insights into dynamic molecular subgroups within GBM, which govern distinct clinical response and sensitivity of each individual to therapy. In the present study, we analyzed genomic composition of primary GBMs between two ethnic groups [IRCR (Institute of Refractory Cancer Research), and TCGA (The Cancer Genome Atlats)] to explore genomic and molecular features that constitute malignant behavior of glioblastoma based on distinct ethnicity. We identified enrichments of MAPK and p53 pathways in IRCR patients, while aberrant activation of Receptor Tyrosine Kinases (RTKs) were predominant in TCGA cohort. We also discovered differential clinical prognosis between two groups and explored essential features that present such diversity.

Multi-Habitat Radiomics Unravels Distinct Phenotypic Subtypes of Glioblastoma with Clinical and Genomic Significance.

We aimed to evaluate the potential of radiomics as an imaging biomarker for glioblastoma (GBM) patients and explore the molecular rationale behind radiomics using a radio-genomics approach. A total of 144 primary GBM patients were included in this study (training cohort). Using multi-parametric MR images, radiomics features were extracted from multi-habitats of the tumor. We applied Cox-LASSO algorithm to build a survival prediction model, which we validated using an independent validation cohort. GBM patients were consensus clustered to reveal inherent phenotypic subtypes. GBM patients were successfully stratified by the radiomics risk score, a weighted sum of radiomics features, corroborating the potential of radiomics as a prognostic biomarker. Using consensus clustering, we identified three distinct subtypes which significantly differed in the prognosis ("heterogenous enhancing", "rim-enhancing necrotic", and "cystic" subtypes). Transcriptomic traits enriched in individual subtypes were in accordance with imaging phenotypes summarized by radiomics. For example, rim-enhancing necrotic subtype was well described by radiomics profiling (T2 autocorrelation and flat shape) and highlighted by the inflammatory genomic signatures, which well correlated to its phenotypic peculiarity (necrosis). This study showed that imaging subtypes derived from radiomics successfully recapitulated the genomic underpinnings of GBMs and thereby confirmed the feasibility of radiomics as an imaging biomarker for GBM patients with comprehensible biologic annotation.

Sphere-Forming Culture for Expanding Genetically Distinct Patient-Derived Glioma Stem Cells by Cellular Growth Rate Screening.

Diffusely infiltrating gliomas (DIGs) are difficult to completely resect and are associated with a high rate of tumor relapse and progression from low- to high-grade glioma. In particular, optimized short-term culture-enriching patient-derived glioma stem cells (GSCs) are essential for customizing the therapeutic strategy based on clinically feasible in vitro drug screening for a wide range of DIGs, owing to the high inter-tumoral heterogeneity. Herein, we constructed a novel high-throughput culture condition screening platform called 'GFSCAN', which evaluated the cellular growth rates of GSCs for each DIG sample in 132 serum-free combinations, using 13 previously reported growth factors closely associated with glioma aggressiveness. In total, 72 patient-derived GSCs with available genomic profiles were tested in GFSCAN to explore the association between cellular growth rates in specific growth factor combinations and genomic/molecular backgrounds, including isocitrate dehydrogenase 1 (IDH1) mutation, chromosome arm 1p and 19q co-deletion, ATRX chromatin remodeler alteration, and transcriptional subtype. GSCs were clustered according to the dependency on epidermal growth factor and basic fibroblast growth factor (E&F), and isocitrate dehydrogenase 1 (IDH1) wild-type GSCs showed higher E&F dependencies than IDH1 mutant GSCs. More importantly, we elucidated optimal combinations for IDH1 mutant glioblastoma and lower grade glioma GSCs with low dependencies on E&F, which could be an aid in clinical decision-making for these DIGs. Thus, we demonstrated the utility of GFSCAN in personalizing in vitro cultivation to nominate personalized therapeutic options, in a clinically relevant time frame, for individual DIG patients, where standard clinical options have been exhausted.

Comprehensive pharmacogenomic characterization of gastric cancer.

BACKGROUND: Gastric cancer is among the most lethal human malignancies. Previous studies have identified molecular aberrations that constitute dynamic biological networks and genomic complexities of gastric tumors. However, the clinical translation of molecular-guided targeted therapy is hampered by challenges. Notably, solid tumors often harbor multiple genetic alterations, complicating the development of effective treatments. METHODS: To address such challenges, we established a comprehensive dataset of molecularly annotated patient derivatives coupled with pharmacological profiles for 60 targeted agents to explore dynamic pharmacogenomic interactions in gastric cancers. RESULTS: We identified lineage-specific drug sensitivities based on histopathological and molecular subclassification, including substantial sensitivities toward VEGFR and EGFR inhibition therapies in diffuse- and signet ring-type gastric tumors, respectively. We identified potential therapeutic opportunities for WNT pathway inhibitors in ALK-mutant tumors, a significant association between PIK3CA-E542K mutation and AZD5363 response, and transcriptome expression of RNF11 as a potential predictor of response to gefitinib. CONCLUSIONS: Collectively, our results demonstrate the feasibility of drug screening combined with tumor molecular characterization to facilitate personalized therapeutic regimens for gastric tumors.

Clinical Targeted Next-Generation sequencing Panels for Detection of Somatic Variants in Gliomas.

PURPOSE: Targeted next-generation sequencing (NGS) panels for solid tumors have been useful in clinical framework for accurate tumor diagnosis and identifying essential molecular aberrations. However, most cancer panels have been designed to address a wide spectrum of pan-cancer models, lacking integral prognostic markers that are highly specific to gliomas. MATERIALS AND METHODS: To address such challenges, we have developed a glioma-specific NGS panel, termed "GliomaSCAN," that is capable of capturing single nucleotide variations and insertion/deletion, copy number variation, and selected promoter mutations and structural variations that cover a subset of intron regions in 232 essential glioma-associated genes. We confirmed clinical concordance rate using pairwise comparison of the identified variants from whole exome sequencing (WES), immunohistochemical analysis, and fluorescence in situ hybridization. RESULTS: Our panel demonstrated high sensitivity in detecting potential genomic variants that were present in the standard materials. To ensure the accuracy of our targeted sequencing panel, we compared our targeted panel to WES. The comparison results demonstrated a high correlation. Furthermore, we evaluated clinical utility of our panel in 46 glioma patients to assess the detection capacity of potential actionable mutations. Thirty-two patients harbored at least one recurrent somatic mutation in clinically actionable gene. CONCLUSION: We have established a glioma-specific cancer panel. GliomaSCAN highly excelled in capturing somatic variations in terms of both sensitivity and specificity and provided potential clinical implication in facilitating genome-based clinical trials. Our results could provide conceptual advance towards improving the response of genomically guided molecularly targeted therapy in glioma patients.

Pharmacogenomic analysis of patient-derived tumor cells in gynecologic cancers.

BACKGROUND: Gynecologic malignancy is one of the leading causes of mortality in female adults worldwide. Comprehensive genomic analysis has revealed a list of molecular aberrations that are essential to tumorigenesis, progression, and metastasis of gynecologic tumors. However, targeting such alterations has frequently led to treatment failures due to underlying genomic complexity and simultaneous activation of various tumor cell survival pathway molecules. A compilation of molecular characterization of tumors with pharmacological drug response is the next step toward clinical application of patient-tailored treatment regimens. RESULTS: Toward this goal, we establish a library of 139 gynecologic tumors including epithelial ovarian cancers (EOCs), cervical, endometrial tumors, and uterine sarcomas that are genomically and/or pharmacologically annotated and explore dynamic pharmacogenomic associations against 37 molecularly targeted drugs. We discover lineage-specific drug sensitivities based on subcategorization of gynecologic tumors and identify TP53 mutation as a molecular determinant that elicits therapeutic response to poly (ADP-Ribose) polymerase (PARP) inhibitor. We further identify transcriptome expression of inhibitor of DNA biding 2 (ID2) as a potential predictive biomarker for treatment response to olaparib. CONCLUSIONS: Together, our results demonstrate the potential utility of rapid drug screening combined with genomic profiling for precision treatment of gynecologic cancers.

Novel Semi-Replicative Retroviral Vector Mediated Double Suicide Gene Transfer Enhances Antitumor Effects in Patient-Derived Glioblastoma Models.

As glioblastomas are mostly localized infiltrative lesions, gene therapy based on the retroviral replicating vector (RRV) system is considered an attractive strategy. Combinations of multiple suicide genes can circumvent the limitations associated with each gene, achieving direct and synergistic cytotoxic effects, along with bystander cell killing. In this study, we constructed a semi-and pseudotyped-RRV (sp-RRV) system harboring two suicide genes-herpes simplex virus type 1 thymidine kinase (TK) and yeast cytosine deaminase (CD)-to verify the dissemination and antitumor efficacy of our sp-RRV system (spRRVe-sEF1α-TK/sRRVgp-sEF1α-CD) in seven patient-derived glioblastoma stem-like cells (GSCs). Flow cytometry and high-content analysis revealed a wide range of transduction efficiency and good correlation between the delivery of therapeutic genes and susceptibility to the prodrugs ganciclovir and 5-fluorocytosine in patient-derived GSCs in vitro. Intra-tumoral delivery of spRRVe-sEF1α-TK/sRRVgp-sEF1α-CD, combined with prodrug treatment, synergistically inhibited cell proliferation and angiogenesis while increasing apoptosis and the depletion of tumor-associated macrophages in orthotopic glioblastoma xenografts. Genomic profiling of patient-derived GSCs revealed that the key genes preventing sp-RRV infection and transmission were associated with cell adhesion, migration, development, differentiation, and proliferation. This is the first report demonstrating that a novel sp-RRV-mediated TK/CD double suicide gene transfer system has high oncolytic power against extremely heterogeneous and treatment-refractory glioblastomas.

Elucidation of Novel Therapeutic Targets for Acute Myeloid Leukemias with RUNX1-RUNX1T1 Fusion.

The RUNX1-RUNX1T1 fusion is a frequent chromosomal alteration in acute myeloid leukemias (AMLs). Although RUNX1-RUNX1T1 fusion protein has pivotal roles in the development of AMLs with the fusion, RUNX1-RUNX1T1, fusion protein is difficult to target, as it lacks kinase activities. Here, we used bioinformatic tools to elucidate targetable signaling pathways in AMLs with RUNX1-RUNX1T1 fusion. After analysis of 93 AML cases from The Cancer Genome Atlas (TCGA) database, we found expression of 293 genes that correlated to the expression of the RUNX1-RUNX1T1 fusion gene. Based on these 293 genes, the cyclooxygenase (COX), vascular endothelial growth factor receptor (VEGFR), platelet-derived growth factor receptor (PDGFR), and fibroblast growth factor receptor (FGFR) pathways were predicted to be specifically activated in AMLs with RUNX1-RUNX1T1 fusion. Moreover, the in vitro proliferation of AML cells with RUNX1-RUNX1T1 fusion decreased significantly more than that of AML cells without the fusion, when the pathways were inhibited pharmacologically. The results indicate that novel targetable signaling pathways could be identified by the analysis of the gene expression features of AMLs with non-targetable genetic alterations. The elucidation of specific molecular targets for AMLs that have a specific genetic alteration would promote personalized treatment of AMLs and improve clinical outcomes.

Phenotypic Plasticity of Invasive Edge Glioma Stem-like Cells in Response to Ionizing Radiation.

Unresectable glioblastoma (GBM) cells in the invading tumor edge can act as seeds for recurrence. The molecular and phenotypic properties of these cells remain elusive. Here, we report that the invading edge and tumor core have two distinct types of glioma stem-like cells (GSCs) that resemble proneural (PN) and mesenchymal (MES) subtypes, respectively. Upon exposure to ionizing radiation (IR), GSCs, initially enriched for a CD133+ PN signature, transition to a CD109+ MES subtype in a C/EBP-ß-dependent manner. Our gene expression analysis of paired cohorts of patients with primary and recurrent GBMs identified a CD133-to-CD109 shift in tumors with an MES recurrence. Patient-derived CD133-/CD109+ cells are highly enriched with clonogenic, tumor-initiating, and radiation-resistant properties, and silencing CD109 significantly inhibits these phenotypes. We also report a conserved regulation of YAP/TAZ pathways by CD109 that could be a therapeutic target in GBM.

Identification of genomic and molecular traits that present therapeutic vulnerability to HGF-targeted therapy in glioblastoma.

BACKGROUND: Cancer is a complex disease with profound genomic alterations and extensive heterogeneity. Recent studies on large-scale genomics have shed light on the impact of core oncogenic pathways, which are frequently dysregulated in a wide spectrum of cancer types. Aberrant activation of the hepatocyte growth factor (HGF) signaling axis has been associated with promoting various oncogenic programs during tumor initiation, progression, and treatment resistance. As a result, HGF-targeted therapy has emerged as an attractive therapeutic approach. However, recent clinical trials involving HGF-targeted therapies have demonstrated rather disappointing results. Thus, an alternative, in-depth assessment of new patient stratification is necessary to shift the current clinical course. METHODS: To address such challenges, we have evaluated the therapeutic efficacy of YYB-101, an HGF-neutralizing antibody, in a series of primary glioblastoma stem cells (GSCs) both in vitro and in vivo. Furthermore, we performed genome and transcriptome analysis to determine genetic and molecular traits that exhibit therapeutic susceptibility to HGF-mediated therapy. RESULTS: We have identified several differentially expressed genes, including MET, KDR, and SOX3, which are associated with tumor invasiveness, malignancy, and unfavorable prognosis in glioblastoma patients. We also demonstrated the HGF-MET signaling axis as a key molecular determinant in GSC invasion, and we discovered that a significant association in HGF expression existed between mesenchymal phenotype and immune cell recruitment. CONCLUSIONS: Upregulation of MET and mesenchymal cellular state are essential in generating HGF-mediated therapeutic responses. Our results provide an important framework for evaluating HGF-targeted therapy in future clinical settings.

Identification of transcriptome signature for predicting clinical response to bevacizumab in recurrent glioblastoma.

Glioblastomas are among the most fatal brain tumors. Although no effective treatment option is available for recurrent glioblastomas (GBMs), a subset of patients evidently derived clinical benefit from bevacizumab, a monoclonal antibody against vascular endothelial growth factor. We retrospectively reviewed patients with recurrent GBM who received bevacizumab to identify biomarkers for predicting clinical response to bevacizumab. Following defined criteria, the patients were categorized into two clinical response groups, and their genetic and transcriptomic results were compared. Angiogenesis-related gene sets were upregulated in both responders and nonresponders, whereas genes for each corresponding angiogenesis pathway were distinct from one another. Two gene sets were made, namely, the nonresponder angiogenesis gene set (NAG) and responder angiogenesis gene set (RAG), and then implemented in independent GBM cohort to validate our dataset. A similar association between the corresponding gene set and survival was observed. In NAG, COL4A2 was associated with a poor clinical outcome in bevacizumab-treated patients. This study demonstrates that angiogenesis-associated gene sets are composed of distinct subsets with diverse biological roles and they represent different clinical responses to anti-angiogenic therapy. Enrichment of a distinct angiogenesis pathway may serve as a biomarker to predict patients who will derive a clinical benefit from bevacizumab.