RNA-sequencing to discover genes and signaling pathways associated with venous thromboembolism in glioblastoma patients: A case-control study

Background: Glioblastoma patients are at high risk of developing venous thromboembolism (VTE). Tumor-intrinsic features are considered to play a role, but the underlying pathophysiological mechanisms remain incompletely understood. Objectives: To identify tumor-expressed genes and signaling pathways that associate with glioblastoma-related VTE by using next generation RNA-sequencing (RNA-Seq). Methods: The tumor gene expression profile of 23 glioblastoma patients with VTE and 23 glioblastoma patients without VTE was compared using an unpaired analysis. Ingenuity Pathway Analysis (IPA) core analysis was performed on the top 50 differentially expressed genes to explore associated functions and pathways. Based on full RNA-Seq data, molecular glioblastoma subtypes were determined by performing cluster analysis. Results: Of the 19,327 genes, 1246 (6.4 %) were differentially expressed between glioblastoma patients with and without VTE (unadjusted P < 0.05). The most highly overexpressed gene was GLI1, a classical target gene in the Sonic Hedgehog (Shh) signaling pathway (log2 fold change: 3.7; unadjusted P < 0.0001, adjusted P = 0.219). In line, Shh signaling was among the top canonical pathways and processes associated with VTE. The proportion of patients with the proneural/neural glioblastoma subtype was higher among those with VTE than controls. Conclusion: Shh signaling may be involved in the development of glioblastoma-related VTE.