Discovering novel cancer bio-markers in acquired lapatinib resistance using Bayesian methods

Signalling transduction pathways (STPs) are commonly hijacked by many cancers for their growth and malignancy, but demystifying their underlying mechanisms is difficult. Here, we developed methodologies with a fully Bayesian approach in discovering novel driver bio-markers in aberrant STPs given high-throughput gene expression (GE) data. This project, namely 'PathTurbEr' (Pathway Perturbation Driver) uses the GE dataset derived from the lapatinib (an EGFR/HER dual inhibitor) sensitive and resistant samples from breast cancer cell lines (SKBR3). Differential expression analysis revealed 512 differentially expressed genes (DEGs) and their pathway enrichment revealed 13 highly perturbed singalling pathways in lapatinib resistance, including PI3K-AKT, Chemokine, Hippo and TGF-$\beta $ singalling pathways. Next, the aberration in TGF-$\beta $ STP was modelled as a causal Bayesian network (BN) using three MCMC sampling methods, i.e. Neighbourhood sampler (NS) and Hit-and-Run (HAR) sampler that potentially yield robust inference with lower chances of getting stuck at local optima and faster convergence compared to other state-of-art methods. Next, we examined the structural features of the optimal BN as a statistical process that generates the global structure using $p_1$-model, a special class of Exponential Random Graph Models (ERGMs), and MCMC methods for their hyper-parameter sampling. This step enabled key drivers identification that drive the aberration within the perturbed BN structure of STP, and yielded 34, 34 and 23 perturbation driver genes out of 80 constituent genes of three perturbed STP models of TGF-$\beta $ signalling inferred by NS, HAR and MH sampling methods, respectively. Functional-relevance and disease-relevance analyses suggested their significant associations with breast cancer progression/resistance.