Molecular Pathway Analysis of Mutation Data for Biomarkers Discovery and Scoring of Target Cancer Drugs

Quantified pathway mutations associate epithelial-mesenchymal transition and immune escape with poor prognosis and immunotherapy resistance of head and neck squamous cell carcinoma

BACKGROUND

Pathway mutations have been calculated to predict the poor prognosis and immunotherapy resistance in head and neck squamous cell carcinoma (HNSCC). To uncover the unique markers predicting prognosis and immune therapy response, the accurate quantification of pathway mutations are required to evaluate epithelial-mesenchymal transition (EMT) and immune escape. Yet, there is a lack of score to accurately quantify pathway mutations.

MATERIAL AND METHODS

Firstly, we proposed Individualized Weighted Hallmark Gene Set Mutation Burden (IWHMB, https://github.com/YuHongHuang-lab/IWHMB ) which integrated pathway structure information and eliminated the interference of global Tumor Mutation Burden to accurately quantify pathway mutations. Subsequently, to further elucidate the association of IWHMB with EMT and immune escape, support vector machine regression model was used to identify IWHMB-related transcriptomic features (IRG), while Adversarially Regularized Graph Autoencoder (ARVGA) was used to further resolve IRG network features. Finally, Random walk with restart algorithm was used to identify biomarkers for predicting ICI response.

RESULTS

We quantified the HNSCC pathway mutation signatures and identified pathway mutation subtypes using IWHMB. The IWHMB-related transcriptomic features (IRG) identified by support vector machine regression were divided into 5 communities by ARVGA, among which the Community 1 enriching malignant mesenchymal components promoted EMT dynamically and regulated immune patterns associated with ICI responses. Bridge Hub Gene (BHG) identified by random walk with restart was key to IWHMB in EMT and immune escape, thus, more predictive for ICI response than other 70 public signatures.

CONCLUSION

In summary, the novel pathway mutation scoring-IWHMB suggested that the elevated malignancy mediated by pathway mutations is a major cause of poor prognosis and immunotherapy failure in HNSCC, and is capable of identifying novel biomarkers to predict immunotherapy response.

OncoboxPD: human 51 672 molecular pathways database with tools for activity calculating and visualization

OncoboxPD (Oncobox pathway databank) available at https://open.oncobox.com is the collection of 51 672 uniformly processed human molecular pathways. Superposition of all pathways formed interactome graph of protein–protein interactions and metabolic reactions containing 361 654 interactions and 64 095 molecular participants. Pathways are uniformly classified by biological processes, and each pathway node is algorithmically functionally annotated by specific activator/repressor role. This enables online calculation of statistically supported pathway activation levels (PALs) with the built-in bioinformatic tool using custom RNA/protein expression profiles. Each pathway can be visualized as static or dynamic graph, where vertices are molecules participating in a pathway and edges are interactions or reactions between them. Differentially expressed nodes in a pathway can be visualized in two-color mode with user-defined color scale. For every comparison, OncoboxPD also generates a graph summarizing top up- and downregulated pathways.

RNA Sequencing for Personalized Treatment of Metastatic Leiomyosarcoma: Case Report

Uterine leiomyosarcoma (UL) is a rare malignant tumor that develops from the uterine smooth muscle tissue. Due to the low frequency and lack of sufficient data from clinical trials there is currently no effective treatment that is routinely accepted for UL. Here we report a case of a 65-years-old female patient with metastatic UL, who progressed on ifosfamide and doxorubicin therapy and developed severe hypertensive crisis after administration of second line pazopanib, which lead to treatment termination. Rapid progression of the tumor stressed the need for the alternative treatment options. We performed RNA sequencing and whole exome sequencing profiling of the patient's biopsy and applied Oncobox bioinformatic algorithm to prioritize targeted therapeutics. No clinically relevant mutations associated with drug efficiencies were found, but the Oncobox transcriptome analysis predicted regorafenib as the most effective targeted treatment option. Regorafenib administration resulted in a complete metabolic response which lasted for 10 months. In addition, RNA sequencing analysis revealed a novel cancer fusion transcript of YWHAE gene with fusion partner JAZF1. Several chimeric transcripts for YWHAE and JAZF1 genes were previously found in uterine neoplasms and some of them were associated with tumor prognosis. However, their combination was detected in this study for the first time. Taken together, these findings evidence that RNA sequencing may complement analysis of clinically relevant mutations and enhance management of oncological patients by suggesting putative treatment options.

Using proteomic and transcriptomic data to assess activation of intracellular molecular pathways

Analysis of molecular pathway activation is the recent instrument that helps to quantize activities of various intracellular signaling, structural, DNA synthesis and repair, and biochemical processes. This may have a deep impact in fundamental research, bioindustry, and medicine. Unlike gene ontology analyses and numerous qualitative methods that can establish whether a pathway is affected in principle, the quantitative approach has the advantage of exactly measuring the extent of a pathway up/downregulation. This results in emergence of a new generation of molecular biomarkers-pathway activation levels, which reflect concentration changes of all measurable pathway components. The input data can be the high-throughput proteomic or transcriptomic profiles, and the output numbers take both positive and negative values and positively reflect overall pathway activation. Due to their nature, the pathway activation levels are more robust biomarkers compared to the individual gene products/protein levels. Here, we review the current knowledge of the quantitative gene expression interrogation methods and their applications for the molecular pathway quantization. We consider enclosed bioinformatic algorithms and their applications for solving real-world problems. Besides a plethora of applications in basic life sciences, the quantitative pathway analysis can improve molecular design and clinical investigations in pharmaceutical industry, can help finding new active biotechnological components and can significantly contribute to the progressive evolution of personalized medicine. In addition to the theoretical principles and concepts, we also propose publicly available software for the use of large-scale protein/RNA expression data to assess the human pathway activation levels.

System, Method and Software for Calculation of a Cannabis Drug Efficiency Index for the Reduction of Inflammation

There are many varieties of Cannabis sativa that differ from each other by composition of cannabinoids, terpenes and other molecules. The medicinal properties of these cultivars are often very different, with some being more efficient than others. This report describes the development of a method and software for the analysis of the efficiency of various cannabis extracts to detect the anti-inflammatory properties of the various cannabis extracts. The method uses high-throughput gene expression profiling data but can potentially use other omics data as well. According to the signaling pathway topology, the gene expression profiles are convoluted into the signaling pathway activities using a signaling pathway impact analysis (SPIA) method. The method was tested by inducing inflammation in human 3D epithelial tissues, including intestine, oral and skin, and then exposing these tissues to various extracts and then performing transcriptome analysis. The analysis showed a different efficiency of the various extracts in restoring the transcriptome changes to the pre-inflammation state, thus allowing to calculate a different cannabis drug efficiency index (CDEI).

Intratumoral Heterogeneity and Longitudinal Changes in Gene Expression Predict Differential Drug Sensitivity in Newly Diagnosed and Recurrent Glioblastoma

Background: Inevitable recurrence after radiochemotherapy is the major problem in the treatment of glioblastoma, the most prevalent type of adult brain malignancy. Glioblastomas are notorious for a high degree of intratumor heterogeneity manifest through a diversity of cell types and molecular patterns. The current paradigm of understanding glioblastoma recurrence is that cytotoxic therapy fails to target effectively glioma stem cells. Recent advances indicate that therapy-driven molecular evolution is a fundamental trait associated with glioblastoma recurrence. There is a growing body of evidence indicating that intratumor heterogeneity, longitudinal changes in molecular biomarkers and specific impacts of glioma stem cells need to be taken into consideration in order to increase the accuracy of molecular diagnostics still relying on readouts obtained from a single tumor specimen. Methods: This study integrates a multisampling strategy, longitudinal approach and complementary transcriptomic investigations in order to identify transcriptomic traits of recurrent glioblastoma in whole-tissue specimens of glioblastoma or glioblastoma stem cells. In this study, 128 tissue samples of 44 tumors including 23 first diagnosed, 19 recurrent and 2 secondary recurrent glioblastomas were analyzed along with 27 primary cultures of glioblastoma stem cells by RNA sequencing. A novel algorithm was used to quantify longitudinal changes in pathway activities and model efficacy of anti-cancer drugs based on gene expression data. Results: Our study reveals that intratumor heterogeneity of gene expression patterns is a fundamental characteristic of not only newly diagnosed but also recurrent glioblastomas. Evidence is provided that glioblastoma stem cells recapitulate intratumor heterogeneity, longitudinal transcriptomic changes and drug sensitivity patterns associated with the state of recurrence. Conclusions: Our results provide a transcriptional rationale for the lack of significant therapeutic benefit from temozolomide in patients with recurrent glioblastoma. Our findings imply that the spectrum of potentially effective drugs is likely to differ between newly diagnosed and recurrent glioblastomas and underscore the merits of glioblastoma stem cells as prognostic models for identifying alternative drugs and predicting drug response in recurrent glioblastoma. With the majority of recurrent glioblastomas being inoperable, glioblastoma stem cell models provide the means of compensating for the limited availability of recurrent glioblastoma specimens.

Mutation Enrichment and Transcriptomic Activation Signatures of 419 Molecular Pathways in Cancer

Carcinogenesis is linked with massive changes in regulation of gene networks. We used high throughput mutation and gene expression data to interrogate involvement of 278 signaling, 72 metabolic, 48 DNA repair and 47 cytoskeleton molecular pathways in cancer. Totally, we analyzed 4910 primary tumor samples with individual cancer RNA sequencing and whole exome sequencing profiles including ~1.3 million DNA mutations and representing thirteen cancer types. Gene expression in cancers was compared with the corresponding 655 normal tissue profiles. For the first time, we calculated mutation enrichment values and activation levels for these pathways. We found that pathway activation profiles were largely congruent among the different cancer types. However, we observed no correlation between mutation enrichment and expression changes both at the gene and at the pathway levels. Overall, positive median cancer-specific activation levels were seen in the DNA repair, versus similar slightly negative values in the other types of pathways. The DNA repair pathways also demonstrated the highest values of mutation enrichment. However, the signaling and cytoskeleton pathways had the biggest proportions of representatives among the outstandingly frequently mutated genes thus suggesting their initiator roles in carcinogenesis and the auxiliary/supporting roles for the other groups of molecular pathways.