Identifying Drug Sensitivity Subnetworks with NETPHIX

Phenotypic heterogeneity in cancer is often caused by different patterns of genetic alterations. Understanding such phenotype-genotype relationships is fundamental for the advance of personalized medicine. One of the important challenges in the area is to predict drug response on a personalized level. The pathway-centric view of cancer significantly advanced the understanding of genotype-phenotype relationships. However, most of network identification methods in cancer focus on identifying subnetworks that include general cancer drivers or are associated with discrete features such as cancer subtypes, hence cannot be applied directly for the analysis of continuous features like drug response. On the other hand, existing genome wide association approaches do not fully utilize the complex proprieties of cancer mutational landscape. To address these challenges, we propose a computational method, named NETPHLIX (NETwork-to-PHenotpe mapping LeveragIng eXlusivity), which aims to identify mutated subnetworks that are associated with drug response (or any continuous cancer phenotype). Utilizing properties such as mutual exclusivity and interactions among genes, we formulate the problem as an integer linear program and solve it optimally to obtain a set of genes satisfying the constraints. NETPHLIX identified gene modules significantly associated with many drugs, including interesting response modules to MEK1/2 inhibitors in both directions (increased and decreased sensitivity to the drug) that the previous method, which does not utilize network information, failed to identify. The genes in the modules belong to MAPK/ERK signaling pathway, which is the targeted pathway of the drug.