Linking signaling pathways to transcriptional programs in breast cancer

STAN, a computational framework for inferring spatially informed transcription factor activity

Transcription factors (TFs) drive significant cellular changes in response to environmental cues and intercellular signaling. Neighboring cells influence TF activity and, consequently, cellular fate and function. Spatial transcriptomics (ST) captures mRNA expression patterns across tissue samples, enabling characterization of the local microenvironment. However, these datasets have not been fully leveraged to systematically estimate TF activity governing cell identity. Here, we present STAN ( S patially informed T ranscription factor A ctivity N etwork), a linear mixed-effects computational method that predicts spot-specific, spatially informed TF activities by integrating curated TF-target gene priors, mRNA expression, spatial coordinates, and morphological features from corresponding imaging data. We tested STAN using lymph node, breast cancer, and glioblastoma ST datasets to demonstrate its applicability by identifying TFs associated with specific cell types, spatial domains, pathological regions, and ligand‒receptor pairs. STAN augments the utility of STs to reveal the intricate interplay between TFs and spatial organization across a spectrum of cellular contexts.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a new optimization procedure to identify the top network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

Metabolic Heterogeneity and Potential Immunotherapeutic Responses Revealed by Single-Cell Transcriptomics of Breast Cancer

BACKGROUND

Breast cancer (BC) exhibits remarkable heterogeneity. However, the transcriptomic heterogeneity of BC at the single-cell level has not been fully elucidated.

METHODS

We acquired BC samples from 14 patients. Single-cell RNA sequencing (scRNA-seq), bioinformatic analyses, along with immunohistochemistry (IHC) and immunofluorescence (IF) assays were carried out.

RESULTS

According to the scRNA-seq results, 10 different cell types were identified. We found that Cancer-Associated Fibroblasts (CAFs) exhibited distinct biological functions and may promote resistance to therapy. Metabolic analysis of tumor cells revealed heterogeneity in glycolysis, gluconeogenesis, and fatty acid synthetase reprogramming, which led to chemotherapy resistance. Furthermore, patients with multiple metastases and progression were predicted to benefit from immunotherapy based on a heterogeneity analysis of T cells and tumor cells.

CONCLUSIONS

Our findings provide a comprehensive understanding of the heterogeneity of BC, provide comprehensive insight into the correlation between cancer metabolism and chemotherapy resistance, and enable the prediction of immunotherapy responses based on T-cell heterogeneity.

Mapping drug biology to disease genetics to discover drug impacts on the human phenome

Motivation

Medications can have unexpected effects on disease, including not only harmful drug side effects, but also beneficial drug repurposing. These effects on disease may result from hidden influences of drugs on disease gene networks. Then, discovering how biological effects of drugs relate to disease biology can both provide insight into the mechanism of latent drug effects, and can help predict new effects.

Results

Here, we develop Draphnet, a model that integrates molecular data on 429 drugs and gene associations of nearly 200 common phenotypes to learn a network that explains drug effects on disease in terms of these molecular signals. We present evidence that our method can both predict drug effects, and can provide insight into the biology of unexpected drug effects on disease. Using Draphnet to map a drug's known molecular effects to downstream effects on the disease genome, we put forward disease genes impacted by drugs, and we suggest a new grouping of drugs based on shared effects on the disease genome. Our approach has multiple applications, including predicting drug uses and learning drug biology, with implications for personalized medicine.

Availability and implementation

Code to reproduce the analysis is available at https://github.com/RDMelamed/drug-phenome.

Data-driven modeling of core gene regulatory network underlying leukemogenesis in IDH mutant AML

Acute myeloid leukemia (AML) is characterized by uncontrolled proliferation of poorly differentiated myeloid cells, with a heterogenous mutational landscape. Mutations in IDH1 and IDH2 are found in 20% of the AML cases. Although much effort has been made to identify genes associated with leukemogenesis, the regulatory mechanism of AML state transition is still not fully understood. To alleviate this issue, here we develop a new computational approach that integrates genomic data from diverse sources, including gene expression and ATAC-seq datasets, curated gene regulatory interaction databases, and mathematical modeling to establish models of context-specific core gene regulatory networks (GRNs) for a mechanistic understanding of tumorigenesis of AML with IDH mutations. The approach adopts a novel optimization procedure to identify the optimal network according to its accuracy in capturing gene expression states and its flexibility to allow sufficient control of state transitions. From GRN modeling, we identify key regulators associated with the function of IDH mutations, such as DNA methyltransferase DNMT1, and network destabilizers, such as E2F1. The constructed core regulatory network and outcomes of in-silico network perturbations are supported by survival data from AML patients. We expect that the combined bioinformatics and systems-biology modeling approach will be generally applicable to elucidate the gene regulation of disease progression.

COVID-19db linkage maps of cell surface proteins and transcription factors in immune cells

The highly contagious SARS-CoV-2 and its associated disease (COVID-19) are a threat to global public health and economies. To develop effective treatments for COVID-19, we must understand the host cell types, cell states and regulators associated with infection and pathogenesis such as dysregulated transcription factors (TFs) and surface proteins, including signaling receptors. To link cell surface proteins with TFs, we recently developed SPaRTAN (Single-cell Proteomic and RNA-based Transcription factor Activity Network) by integrating parallel single-cell proteomic and transcriptomic data based on Cellular Indexing of Transcriptomes and Epitopes by sequencing (CITE-seq) and gene cis-regulatory information. We apply SPaRTAN to CITE-seq data sets from patients with varying degrees of COVID-19 severity and healthy controls to identify the associations between surface proteins and TFs in host immune cells. Here, we present COVID-19db of Immune Cell States (https://covid19db.streamlit.app/), a web server containing cell surface protein expression, SPaRTAN-inferred TF activities, and their associations with major host immune cell types. The data include four high-quality COVID-19 CITE-seq data sets with a toolset for user-friendly data analysis and visualization. We provide interactive surface protein and TF visualizations across major immune cell types for each data set, allowing comparison between various patient severity groups for the discovery of potential therapeutic targets and diagnostic biomarkers.

Linking Expression of Cell-Surface Receptors with Transcription Factors by Computational Analysis of Paired Single-Cell Proteomes and Transcriptomes

Complex signaling and transcriptional programs control the development and physiology of specialized cell types. Genetic perturbations in these programs cause human cancers to arise from a diverse set of specialized cell types and developmental states. Understanding these complex systems and their potential to drive cancer is critical for the development of immunotherapies and druggable targets. Pioneering single-cell multi-omics technologies that analyze transcriptional states have been coupled with the expression of cell-surface receptors. This chapter describes SPaRTAN (Single-cell Proteomic and RNA-based Transcription factor Activity Network), a computational framework, to link transcription factors with cell-surface protein expression. SPaRTAN uses CITE-seq (cellular indexing of transcriptomes and epitopes by sequencing) data and cis-regulatory sites to model the effect of interactions between transcription factors and cell-surface receptors on gene expression. We demonstrate the pipeline for SPaRTAN using CITE-seq data from peripheral blood mononuclear cells.

Deregulated transcription factors and poor clinical outcomes in cancer patients

Transcription factors are a group of proteins, which possess DNA-binding domains, bind to DNA strands of promoters or enhancers, and initiate transcription of genes with cooperation of RNA polymerase and other co-factors. They play crucial roles in regulating transcription during embryogenesis and development. Their physiological status in different cell types is also important to maintain cellular homeostasis. Therefore, any deregulation of transcription factors will lead to the development of cancer cells and tumor progression. Based on their functions in cancer cells, transcription factors could be either oncogenic or tumor suppressive. Furthermore, transcription factors have been shown to modulate cancer stem cells, epithelial-mesenchymal transition (EMT) and drug response; therefore, measuring deregulated transcription factors is hypothesized to predict treatment outcomes of patients with cancers and targeting deregulated transcription factors could be an encouraging strategy for cancer therapy. Here, we summarize the current knowledge of major deregulated transcription factors and their effects on causing poor clinical outcome of patients with cancer. The information presented here will help to predict the prognosis and drug response and to design novel drugs and therapeutic strategies for the treatment of cancers by targeting deregulated transcription factors.

Interpretable deep learning for chromatin-informed inference of transcriptional programs driven by somatic alterations across cancers

Cancer is a disease of gene dysregulation, where cells acquire somatic and epigenetic alterations that drive aberrant cellular signaling. These alterations adversely impact transcriptional programs and cause profound changes in gene expression. Interpreting somatic alterations within context-specific transcriptional programs will facilitate personalized therapeutic decisions but is a monumental task. Toward this goal, we develop a partially interpretable neural network model called Chromatin-informed Inference of Transcriptional Regulators Using Self-attention mechanism (CITRUS). CITRUS models the impact of somatic alterations on transcription factors and downstream transcriptional programs. Our approach employs a self-attention mechanism to model the contextual impact of somatic alterations. Furthermore, CITRUS uses a layer of hidden nodes to explicitly represent the state of transcription factors (TFs) to learn the relationships between TFs and their target genes based on TF binding motifs in the open chromatin regions of tumor samples. We apply CITRUS to genomic, transcriptomic, and epigenomic data from 17 cancer types profiled by The Cancer Genome Atlas. CITRUS predicts patient-specific TF activities and reveals transcriptional program variations between and within tumor types. We show that CITRUS yields biological insights into delineating TFs associated with somatic alterations in individual tumors. Thus, CITRUS is a promising tool for precision oncology.

Expression and prognostic potential of TMEM204: a pan-cancer analysis

TMEM204 (Transmembrane Protein 204) is a member of the TMEM family that regulates cell function and angiogenesis. Previous studies showed that TMEM204 is related to pancreatic cancer, but its roles in other cancers remain unknown. To reveal this relationship, we conducted a pan-cancer analysis by several online databases. The expression of TMEM204 was analyzed by Oncomine and Tumor Immune Estimation Resource2.0 (TIMER2.0). The prognostic potential of TMEM204 was evaluated by the GEPIA2, UALCAN, and Oncolnc. The methylation level of gene expression was analyzed by UALCAN, and the relationship between cancer and immune invasion was displayed by TIMER2.0. The Protein-Protein Interactions Network and functional analysis of TMEM204 and its related genes were conducted by STRING and Webgestalt. We found that TMEM204 expression was up-regulated and correlated with prognosis in multiple cancers. In liver hepatocellular carcinoma (LIHC), high TMEM204 expression was associated with a good prognosis, and with high infiltrating levels of CD8⁺ T and CD4⁺ T cells, macrophages, neutrophils, and myeloid dendritic cells. In addition, the methylation level in LIHC was higher than in normal tissues. p53 signaling pathway and Fanconi anemia pathway were implicated by KEGG pathway analysis. These results indicate that TMEM204 is associated with the prognosis, methylation, and immune invasion of cancers, especially LIHC. TMEM204 may act as a prognostic marker of LIHC and its role in other cancers should be studied.

Network modeling predicts personalized gene expression and drug responses in valve myofibroblasts cultured with patient sera

Aortic valve stenosis (AVS) patients experience pathogenic valve leaflet stiffening due to excessive extracellular matrix (ECM) remodeling. Numerous microenvironmental cues influence pathogenic expression of ECM remodeling genes in tissue-resident valvular myofibroblasts, and the regulation of complex myofibroblast signaling networks depends on patient-specific extracellular factors. Here, we combined a manually curated myofibroblast signaling network with a data-driven transcription factor network to predict patient-specific myofibroblast gene expression signatures and drug responses. Using transcriptomic data from myofibroblasts cultured with AVS patient sera, we produced a large-scale, logic-gated differential equation model in which 11 biochemical and biomechanical signals were transduced via a network of 334 signaling and transcription reactions to accurately predict the expression of 27 fibrosis-related genes. Correlations were found between personalized model-predicted gene expression and AVS patient echocardiography data, suggesting links between fibrosis-related signaling and patient-specific AVS severity. Further, global network perturbation analyses revealed signaling molecules with the most influence over network-wide activity, including endothelin 1 (ET1), interleukin 6 (IL6), and transforming growth factor β (TGFβ), along with downstream mediators c-Jun N-terminal kinase (JNK), signal transducer and activator of transcription (STAT), and reactive oxygen species (ROS). Lastly, we performed virtual drug screening to identify patient-specific drug responses, which were experimentally validated via fibrotic gene expression measurements in valvular interstitial cells cultured with AVS patient sera and treated with or without bosentan-a clinically approved ET1 receptor inhibitor. In sum, our work advances the ability of computational approaches to provide a mechanistic basis for clinical decisions including patient stratification and personalized drug screening.

SPaRTAN, a computational framework for linking cell-surface receptors to transcriptional regulators

The identity and functions of specialized cell types are dependent on the complex interplay between signaling and transcriptional networks. Recently single-cell technologies have been developed that enable simultaneous quantitative analysis of cell-surface receptor expression with transcriptional states. To date, these datasets have not been used to systematically develop cell-context-specific maps of the interface between signaling and transcriptional regulators orchestrating cellular identity and function. We present SPaRTAN (Single-cell Proteomic and RNA based Transcription factor Activity Network), a computational method to link cell-surface receptors to transcription factors (TFs) by exploiting cellular indexing of transcriptomes and epitopes by sequencing (CITE-seq) datasets with cis-regulatory information. SPaRTAN is applied to immune cell types in the blood to predict the coupling of signaling receptors with cell context-specific TFs. Selected predictions are validated by prior knowledge and flow cytometry analyses. SPaRTAN is then used to predict the signaling coupled TF states of tumor infiltrating CD8+ T cells in malignant peritoneal and pleural mesotheliomas. SPaRTAN enhances the utility of CITE-seq datasets to uncover TF and cell-surface receptor relationships in diverse cellular states.

Bioinformatics Approaches for Anti-cancer Drug Discovery

Drug discovery is important in cancer therapy and precision medicines. Traditional approaches of drug discovery are mainly based on in vivo animal experiments and in vitro drug screening, but these methods are usually expensive and laborious. In the last decade, omics data explosion provides an opportunity for computational prediction of anti-cancer drugs, improving the efficiency of drug discovery. High-throughput transcriptome data were widely used in biomarkers' identification and drug prediction by integrating with drug-response data. Moreover, biological network theory and methodology were also successfully applied to the anti-cancer drug discovery, such as studies based on protein-protein interaction network, drug-target network and disease-gene network. In this review, we summarized and discussed the bioinformatics approaches for predicting anti-cancer drugs and drug combinations based on the multi-omic data, including transcriptomics, toxicogenomics, functional genomics and biological network. We believe that the general overview of available databases and current computational methods will be helpful for the development of novel cancer therapy strategies.

Integrative network analysis identifies potential targets and drugs for ovarian cancer

BACKGROUND

Though accounts for 2.5% of all cancers in female, the death rate of ovarian cancer is high, which is the fifth leading cause of cancer death (5% of all cancer death) in female. The 5-year survival rate of ovarian cancer is less than 50%. The oncogenic molecular signaling of ovarian cancer are complicated and remain unclear, and there is a lack of effective targeted therapies for ovarian cancer treatment.

METHODS

In this study, we propose to investigate activated signaling pathways of individual ovarian cancer patients and sub-groups; and identify potential targets and drugs that are able to disrupt the activated signaling pathways. Specifically, we first identify the up-regulated genes of individual cancer patients using Markov chain Monte Carlo (MCMC), and then identify the potential activated transcription factors. After dividing ovarian cancer patients into several sub-groups sharing common transcription factors using K-modes method, we uncover the up-stream signaling pathways of activated transcription factors in each sub-group. Finally, we mapped all FDA approved drugs targeting on the upstream signaling.

RESULTS

The 427 ovarian cancer samples were divided into 3 sub-groups (with 100, 172, 155 samples respectively) based on the activated TFs (with 14, 25, 26 activated TFs respectively). Multiple up-stream signaling pathways, e.g., MYC, WNT, PDGFRA (RTK), PI3K, AKT TP53, and MTOR, are uncovered to activate the discovered TFs. In addition, 66 FDA approved drugs were identified targeting on the uncovered core signaling pathways. Forty-four drugs had been reported in ovarian cancer related reports. The signaling diversity and heterogeneity can be potential therapeutic targets for drug combination discovery.

CONCLUSIONS

The proposed integrative network analysis could uncover potential core signaling pathways, targets and drugs for ovarian cancer treatment.

Knowledge-primed neural networks enable biologically interpretable deep learning on single-cell sequencing data

BACKGROUND

Deep learning has emerged as a versatile approach for predicting complex biological phenomena. However, its utility for biological discovery has so far been limited, given that generic deep neural networks provide little insight into the biological mechanisms that underlie a successful prediction. Here we demonstrate deep learning on biological networks, where every node has a molecular equivalent, such as a protein or gene, and every edge has a mechanistic interpretation, such as a regulatory interaction along a signaling pathway.

RESULTS

With knowledge-primed neural networks (KPNNs), we exploit the ability of deep learning algorithms to assign meaningful weights in multi-layered networks, resulting in a widely applicable approach for interpretable deep learning. We present a learning method that enhances the interpretability of trained KPNNs by stabilizing node weights in the presence of redundancy, enhancing the quantitative interpretability of node weights, and controlling for uneven connectivity in biological networks. We validate KPNNs on simulated data with known ground truth and demonstrate their practical use and utility in five biological applications with single-cell RNA-seq data for cancer and immune cells.

CONCLUSIONS

We introduce KPNNs as a method that combines the predictive power of deep learning with the interpretability of biological networks. While demonstrated here on single-cell sequencing data, this method is broadly relevant to other research areas where prior domain knowledge can be represented as networks.

Inference of Bacterial Small RNA Regulatory Networks and Integration with Transcription Factor-Driven Regulatory Networks

Small noncoding RNAs (sRNAs) are key regulators of bacterial gene expression. Through complementary base pairing, sRNAs affect mRNA stability and translation efficiency. Here, we describe a network inference approach designed to identify sRNA-mediated regulation of transcript levels. We use existing transcriptional data sets and prior knowledge to infer sRNA regulons using our network inference tool, the Inferelator. This approach produces genome-wide gene regulatory networks that include contributions by both transcription factors and sRNAs. We show the benefits of estimating and incorporating sRNA activities into network inference pipelines using available experimental data. We also demonstrate how these estimated sRNA regulatory activities can be mined to identify the experimental conditions where sRNAs are most active. We uncover 45 novel experimentally supported sRNA-mRNA interactions in Escherichia coli, outperforming previous network-based efforts. Additionally, our pipeline complements sequence-based sRNA-mRNA interaction prediction methods by adding a data-driven filtering step. Finally, we show the general applicability of our approach by identifying 24 novel, experimentally supported, sRNA-mRNA interactions in Pseudomonas aeruginosa, Staphylococcus aureus, and Bacillus subtilis. Overall, our strategy generates novel insights into the functional context of sRNA regulation in multiple bacterial species.

IMPORTANCE Individual bacterial genomes can have dozens of small noncoding RNAs with largely unexplored regulatory functions. Although bacterial sRNAs influence a wide range of biological processes, including antibiotic resistance and pathogenicity, our current understanding of sRNA-mediated regulation is far from complete. Most of the available information is restricted to a few well-studied bacterial species; and even in those species, only partial sets of sRNA targets have been characterized in detail. To close this information gap, we developed a computational strategy that takes advantage of available transcriptional data and knowledge about validated and putative sRNA-mRNA interactions for inferring expanded sRNA regulons. Our approach facilitates the identification of experimentally supported novel interactions while filtering out false-positive results. Due to its data-driven nature, our method prioritizes biologically relevant interactions among lists of candidate sRNA-target pairs predicted in silico from sequence analysis or derived from sRNA-mRNA binding experiments.

Logic-based analysis of gene expression data predicts association between TNF, TGFB1 and EGF pathways in basal-like breast cancer

For breast cancer, clinically important subtypes are well characterized at the molecular level in terms of gene expression profiles. In addition, signaling pathways in breast cancer have been extensively studied as therapeutic targets due to their roles in tumor growth and metastasis. However, it is challenging to put signaling pathways and gene expression profiles together to characterize biological mechanisms of breast cancer subtypes since many signaling events result from post-translational modifications, rather than gene expression differences. We designed a logic-based computational framework to explain the differences in gene expression profiles among breast cancer subtypes using Pathway Logic and transcriptional network information. Pathway Logic is a rewriting-logic-based formal system for modeling biological pathways including post-translational modifications. Our method demonstrated its utility by constructing subtype-specific path from key receptors (TNFR, TGFBR1 and EGFR) to key transcription factor (TF) regulators (RELA, ATF2, SMAD3 and ELK1) and identifying potential association between pathways via TFs in basal-specific paths, which could provide a novel insight on aggressive breast cancer subtypes. Codes and results are available at http://epigenomics.snu.ac.kr/PL/.

KnockTF: a comprehensive human gene expression profile database with knockdown/knockout of transcription factors

Transcription factors (TFs) and their target genes have important functions in human diseases and biological processes. Gene expression profile analysis before and after knockdown or knockout is one of the most important strategies for obtaining target genes of TFs and exploring TF functions. Human gene expression profile datasets with TF knockdown and knockout are accumulating rapidly. Based on the urgent need to comprehensively and effectively collect and process these data, we developed KnockTF (http://www.licpathway.net/KnockTF/index.html), a comprehensive human gene expression profile database of TF knockdown and knockout. KnockTF provides a number of resources for human gene expression profile datasets associated with TF knockdown and knockout and annotates TFs and their target genes in a tissue/cell type-specific manner. The current version of KnockTF has 570 manually curated RNA-seq and microarray datasets associated with 308 TFs disrupted by different knockdown and knockout techniques and across multiple tissue/cell types. KnockTF collects upstream pathway information of TFs and functional annotation results of downstream target genes. It provides details about TFs binding to promoters, super-enhancers and typical enhancers of target genes. KnockTF constructs a TF-differentially expressed gene network and performs network analyses for genes of interest. KnockTF will help elucidate TF-related functions and potential biological effects.

Investigating the effect of dependence between conditions with Bayesian Linear Mixed Models for motif activity analysis

MOTIVATION

Cellular identity and behavior is controlled by complex gene regulatory networks. Transcription factors (TFs) bind to specific DNA sequences to regulate the transcription of their target genes. On the basis of these TF motifs in cis-regulatory elements we can model the influence of TFs on gene expression. In such models of TF motif activity the data is usually modeled assuming a linear relationship between the motif activity and the gene expression level. A commonly used method to model motif influence is based on Ridge Regression. One important assumption of linear regression is the independence between samples. However, if samples are generated from the same cell line, tissue, or other biological source, this assumption may be invalid. This same assumption of independence is also applied to different yet similar experimental conditions, which may also be inappropriate. In theory, the independence assumption between samples could lead to loss in signal detection. Here we investigate whether a Bayesian model that allows for correlations results in more accurate inference of motif activities.

RESULTS

We extend the Ridge Regression to a Bayesian Linear Mixed Model, which allows us to model dependence between different samples. In a simulation study, we investigate the differences between the two model assumptions. We show that our Bayesian Linear Mixed Model implementation outperforms Ridge Regression in a simulation scenario where the noise, which is the signal that can not be explained by TF motifs, is uncorrelated. However, we demonstrate that there is no such gain in performance if the noise has a similar covariance structure over samples as the signal that can be explained by motifs. We give a mathematical explanation to why this is the case. Using four representative real datasets we show that at most ∼​40% of the signal is explained by motifs using the linear model. With these data there is no advantage to using the Bayesian Linear Mixed Model, due to the similarity of the covariance structure.

AVAILABILITY & IMPLEMENTATION

The project implementation is available at https://github.com/Sim19/SimGEXPwMotifs.

An integrative method to predict signalling perturbations for cellular transitions

Induction of specific cellular transitions is of clinical importance, as it allows to revert disease cellular phenotype, or induce cellular reprogramming and differentiation for regenerative medicine. Signalling is a convenient way to accomplish such transitions without transfer of genetic material. Here we present the first general computational method that systematically predicts signalling molecules, whose perturbations induce desired cellular transitions. This probabilistic method integrates gene regulatory networks (GRNs) with manually-curated signalling pathways obtained from MetaCore from Clarivate Analytics, to model how signalling cues are received and processed in the GRN. The method was applied to 219 cellular transition examples, including cell type transitions, and overall correctly predicted experimentally validated signalling molecules, consistently outperforming other well-established approaches, such as differential gene expression and pathway enrichment analyses. Further, we validated our method predictions in the case of rat cirrhotic liver, and identified the activation of angiopoietins receptor Tie2 as a potential target for reverting the disease phenotype. Experimental results indicated that this perturbation induced desired changes in the gene expression of key TFs involved in fibrosis and angiogenesis. Importantly, this method only requires gene expression data of the initial and desired cell states, and therefore is suited for the discovery of signalling interventions for disease treatments and cellular therapies.

Inferring the transcriptional regulatory mechanism of signal-dependent gene expression via an integrative computational approach

Eukaryotic transcription factors (TFs) coordinate different upstream signals to regulate the expression of their target genes. To unveil this regulatory network in B-cell receptor signaling, we developed a computational pipeline to systematically analyze the extracellular signal-regulated kinase (ERK)- and IκB kinase (IKK)-dependent transcriptome responses. We combined a bilinear regression method and kinetic modeling to identify the signal-to-TF and TF-to-gene dynamics, respectively. We input a set of time-course experimental data for B cells and concentrated on transcriptional activators. The results show that the combination of TFs differentially controlled by ERK and IKK could contribute divergent expression dynamics in orchestrating the B-cell response. Our findings provide insights into the regulatory mechanisms underlying signal-dependent gene expression in eukaryotic cells.

Chromatin-informed inference of transcriptional programs in gynecologic and basal breast cancers

Chromatin accessibility data can elucidate the developmental origin of cancer cells and reveal the enhancer landscape of key oncogenic transcriptional regulators. We develop a computational strategy called PSIONIC (patient-specific inference of networks informed by chromatin) to combine chromatin accessibility data with large tumor expression data and model the effect of enhancers on transcriptional programs in multiple cancers. We generate a new ATAC-seq data profiling chromatin accessibility in gynecologic and basal breast cancer cell lines and apply PSIONIC to 723 patient and 96 cell line RNA-seq profiles from ovarian, uterine, and basal breast cancers. Our computational framework enables us to share information across tumors to learn patient-specific TF activities, revealing regulatory differences between and within tumor types. PSIONIC-predicted activity for MTF1 in cell line models correlates with sensitivity to MTF1 inhibition, showing the potential of our approach for personalized therapy. Many identified TFs are significantly associated with survival outcome. To validate PSIONIC-derived prognostic TFs, we perform immunohistochemical analyses in 31 uterine serous tumors for ETV6 and 45 basal breast tumors for MITF and confirm that the corresponding protein expression patterns are also significantly associated with prognosis.

m1A Regulated Genes Modulate PI3K/AKT/mTOR and ErbB Pathways in Gastrointestinal Cancer

BACKGROUND: Gene expression can be posttranscriptionally regulated by a complex network of proteins. N1-methyladenosine (m1A) is a newly validated RNA modification. However, little is known about both its influence and biogenesis in tumor development. METHODS: This study analyzed TCGA data of patients with five kinds of gastrointestinal (GI) cancers. Using data from cBioPortal, molecular features of the nine known m1A-related enzymes in GI cancers were investigated. Using a variety of bioinformatics approach, the impact of m1A regulators on its downstream signaling pathway was studied. To further confirm this regulation, the effect of m1A writer ALKBH3 knockdown was studied using RNA-seq data from published database. RESULTS: Dysregulation and multiple types of genetic alteration of putative m1A-related enzymes in tumor samples were observed. The ErbB and mTOR pathways with ErbB2, mTOR, and AKT1S1 hub genes were identified as being regulated by m1A-related enzymes. The expression of both ErbB2 and AKT1S1 was decreased after m1A writer ALKBH3 knockdown. Furthermore, Gene Ontology analysis revealed that m1A downstream genes were associated with cell proliferation, and the results showed that m1A genes are reliably linked to mTOR. CONCLUSION: This study demonstrated for the first time the dysregulation of m1A regulators in GI cancer and its signaling pathways and will contribute to the understanding of RNA modification in cancer.

Oncogenic MSH6-CXCR4-TGFB1 Feedback Loop: A Novel Therapeutic Target of Photothermal Therapy in Glioblastoma Multiforme

Glioblastoma multiforme (GBM) has been considered the most aggressive glioma type. Temozolomide (TMZ) is the main first-line chemotherapeutic agent for GBM. Decreased mutS homolog 6 (MSH6) expression is clinically recognized as one of the principal reasons for GBM resistance to TMZ. However, the specific functions of MSH6 in GBM, in addition to its role in mismatch repair, remain unknown.

Methods: Bioinformatics were employed to analyze MSH6 mRNA and protein levels in GBM clinical samples and to predict the potential cancer-promoting functions and mechanisms of MSH6. MSH6 levels were silenced or overexpressed in GBM cells to assess its functional effects in vitro and in vivo. Western blot, qRT-PCR, and immunofluorescence assays were used to explore the relevant molecular mechanisms. Cu₂(OH)PO₄@PAA nanoparticles were fabricated through a hydrothermal method. Their MRI and photothermal effects as well as their effect on restraining the MSH6-CXCR4-TGFB1 feedback loop were investigated in vitro and in vivo.

Results: We demonstrated that MSH6 is an overexpressed oncogene in human GBM tissues. MSH6, CXCR4 and TGFB1 formed a triangular MSH6-CXCR4-TGFB1 feedback loop that accelerated gliomagenesis, proliferation (G1 phase), migration and invasion (epithelial-to-mesenchymal transition; EMT), stemness, angiogenesis and antiapoptotic effects by regulating the p-STAT3/Slug and p-Smad2/3/ZEB2 signaling pathways in GBM. In addition, the MSH6-CXCR4-TGFB1 feedback loop was a vital marker of GBM, making it a promising therapeutic target. Notably, photothermal therapy (PTT) mediated by Cu₂(OH)PO₄@PAA + near infrared (NIR) irradiation showed outstanding therapeutic effects, which might be associated with a repressed MSH6-CXCR4-TGFB1 feedback loop and its downstream factors in GBM. Simultaneously, the prominent MR imaging (T1WI) ability of Cu₂(OH)PO₄@PAA could provide visual guidance for PTT.

Conclusions: Our findings indicate that the oncogenic MSH6-CXCR4-TGFB1 feedback loop is a novel therapeutic target for GBM and that PTT is associated with the inhibition of the MSH6-CXCR4-TGFB1 loop.

The Proliferative and Apoptotic Landscape of Basal-like Breast Cancer

Basal-like breast cancer (BLBC) is an aggressive molecular subtype that represents up to 15% of breast cancers. It occurs in younger patients, and typically shows rapid development of locoregional and distant metastasis, resulting in a relatively high mortality rate. Its defining features are that it is positive for basal cytokeratins and, epidermal growth factor receptor and/or c-Kit. Problematically, it is typically negative for the estrogen receptor and human epidermal growth factor receptor 2 (HER2), which means that it is unsuitable for either hormone therapy or targeted HER2 therapy. As a result, there are few therapeutic options for BLBC, and a major priority is to define molecular subgroups of BLBC that could be targeted therapeutically. In this review, we focus on the highly proliferative and anti-apoptotic phenotype of BLBC with the goal of defining potential therapeutic avenues, which could take advantage of these aspects of tumor development.

HMGA1a induces alternative splicing of estrogen receptor alpha in MCF-7 human breast cancer cells

The high-mobility group A protein 1a (HMGA1a) protein is known as an oncogene whose expression level in cancer tissue correlates with the malignant potential, and known as a component of senescence-related structures connecting it to tumor suppressor networks in fibroblasts. HMGA1 protein binds to DNA, but recent studies have shown it exerts novel functions through RNA-binding. Our previous studies have shown that sequence-specific RNA-binding of HMGA1a induces exon-skipping of Presenilin-2 exon 5 in sporadic Alzheimer disease. Here we show that HMGA1a induced exon-skipping of the estrogen receptor alpha (ERα) gene and increased ERα46 mRNA expression in MCF-7 breast cancer cells. An RNA-decoy of HMGA1a efficiently blocked this event and reduced ERα46 protein expression. Blockage of HMGA1a RNA-binding property consequently induced cell growth through reduced ERα46 expression in MCF-7 cells and increased sensitivity to tamoxifen in the tamoxifen-resistant cell line, MCF-7/TAMR1. Stable expression of an HMGA1a RNA-decoy in MCF-7 cells exhibited decreased ERα46 protein expression and increased estrogen-dependent tumor growth when these cells were implanted in nude mice. These results show HMGA1a is involved in alternative splicing of the ERα gene and related to estrogen-related growth as well as tamoxifen sensitivity in MCF-7 breast cancer cells.

Core signaling pathways in ovarian cancer stem cell revealed by integrative analysis of multi-marker genomics data

Tumor recurrence occurs in more than 70% of ovarian cancer patients, and the majority eventually becomes refractory to treatments. Ovarian Cancer Stem Cells (OCSCs) are believed to be responsible for the tumor relapse and drug resistance. Therefore, eliminating ovarian CSCs is important to improve the prognosis of ovarian cancer patients. However, there is a lack of effective drugs to eliminate OCSCs because the core signaling pathways regulating OCSCs remain unclear. Also it is often hard for biologists to identify a few testable targets and infer driver signaling pathways regulating CSCs from a large number of differentially expression genes in an unbiased manner. In this study, we propose a straightforward and integrative analysis to identify potential core signaling pathways of OCSCs by integrating transcriptome data of OCSCs isolated based on two distinctive markers, ALDH and side population, with regulatory network (Transcription Factor (TF) and Target Interactome) and signaling pathways. We first identify the common activated TFs in two OCSC populations integrating the gene expression and TF-target Interactome; and then uncover up-stream signaling cascades regulating the activated TFs. In specific, 22 activated TFs are identified. Through literature search validation, 15 of them have been reported in association with cancer stem cells. Additionally, 10 TFs are found in the KEGG signaling pathways, and their up-stream signaling cascades are extracted, which also provide potential treatment targets. Moreover, 40 FDA approved drugs are identified to target on the up-stream signaling cascades, and 15 of them have been reported in literatures in cancer stem cell treatment. In conclusion, the proposed approach can uncover the activated up-stream signaling, activated TFs and up-regulated target genes that constitute the potential core signaling pathways of ovarian CSC. Also drugs and drug combinations targeting on the core signaling pathways might be able to eliminate OCSCs. The proposed approach can also be applied for identifying potential activated signaling pathways of other types of cancers.

Impact of suppression of tumorigenicity 14 (ST14)/serine protease 14 (Prss14) expression analysis on the prognosis and management of estrogen receptor negative breast cancer

To elucidate the role of a type II transmembrane serine protease, ST14/Prss14, during breast cancer progression, we utilized publically accessible databases including TCGA, GEO, NCI-60, and CCLE. Survival of breast cancer patients with high ST14/Prss14 expression is significantly poor in estrogen receptor (ER) negative populations regardless of the ratios of ST14/Prss14 to its inhibitors, SPINT1 or SPINT2. In a clustering of 1085 selected EMT signature genes, ST14/Prss14 is located in the same cluster with CDH3, and closer to post-EMT markers, CDH2, VIM, and FN1 than to the pre-EMT marker, CDH1. Coexpression analyses of known ST14/Prss14 substrates and transcription factors revealed context dependent action. In cell lines, paradoxically, ST14/Prss14 expression is higher in the ER positive group and located closer to CDH1 in clustering. This apparent contradiction is not likely due to ST14/Prss14 expression in a cancer microenvironment, nor due to negative regulation by ER. Genes consistently coexpressed with ST14/Prss14 include transcription factors, ELF5, GRHL1, VGLL1, suggesting currently unknown mechanisms for regulation. Here, we report that ST14/Prss14 is an emerging therapeutic target for breast cancer where HER2 is not applicable. In addition we suggest that careful conclusions should be drawn not exclusively from the cell line studies for target development.

A signal-based method for finding driver modules of breast cancer metastasis to the lung

Tumor metastasis is mainly caused by somatic genomic alterations (SGAs) that perturb pathways regulating metastasis-relevant activities and thus help the primary tumor to adapt to the new microenvironment. Identifying drivers of metastasis, i.e. SGAs, sheds light on the metastasis mechanism and provides guidance for targeted therapy. In this paper, we introduce a novel method to search for SGAs driving breast cancer metastasis to the lung. First, we search for transcriptomic modules with genes that are differentially expressed in breast cell lines with strong metastatic activities to the lung and co-expressed in a large number of breast tumors. Then, for each transcriptomic module, we search for a set of SGA genes (driver modules) such that genes in each driver module carry a common signal regulating the transcriptomic module. Evaluations indicate that many genes in driver modules are indeed related to metastasis, and our methods have identified many new driver candidates. We further choose two novel metastatic driver genes, BCL2L11 and CDH9, for in vitro verification. The wound healing assay reveals that inhibiting either BCL2L11 or CDH9 will enhance the migration of cell lines, which provides evidence that these two genes are suppressors of tumor metastasis.

Overexpression of the E2F target gene CENPI promotes chromosome instability and predicts poor prognosis in estrogen receptor-positive breast cancer

During cell division, chromosome segregation is facilitated by the mitotic checkpoint, or spindle assembly checkpoint (SAC), which ensures correct kinetochore-microtubule attachments and prevents premature sister-chromatid separation. It is well established that misexpression of SAC components on the outer kinetochores promotes chromosome instability (CIN) and tumorigenesis. Here, we study the expression of CENP-I, a key component of the HIKM complex at the inner kinetochores, in breast cancer, including ductal, lobular, medullary and male breast carcinomas. CENPI mRNA and protein levels are significantly elevated in estrogen receptor-positive (ER+) but not in estrogen receptor-negative (ER-) breast carcinoma. Well-established prognostic tests indicate that CENPI overexpression constitutes a powerful independent marker for poor patient prognosis and survival in ER+ breast cancer. We further demonstrate that CENPI is an E2F target gene. Consistently, it is overexpressed in RB1-deficient breast cancers. However, CENP-I overexpression is not purely due to cell cycle-associated expression. In ER+ breast cancer cells, CENP-I overexpression promotes CIN, especially chromosome gains. In addition, in ER+ breast carcinomas the degree of CENPI overexpression is proportional to the level of aneuploidy and CENPI overexpression is one of the strongest markers for CIN identified to date. Our results indicate that overexpression of the inner kinetochore protein CENP-I promotes CIN and forecasts poor prognosis for ER+ breast cancer patients. These observations provide novel mechanistic insights and have important implications for breast cancer diagnostics and potentially therapeutic targeting.

Seed-effect modeling improves the consistency of genome-wide loss-of-function screens and identifies synthetic lethal vulnerabilities in cancer cells

BACKGROUND

Genome-wide loss-of-function profiling is widely used for systematic identification of genetic dependencies in cancer cells; however, the poor reproducibility of RNA interference (RNAi) screens has been a major concern due to frequent off-target effects. Currently, a detailed understanding of the key factors contributing to the sub-optimal consistency is still a lacking, especially on how to improve the reliability of future RNAi screens by controlling for factors that determine their off-target propensity.

METHODS

We performed a systematic, quantitative analysis of the consistency between two genome-wide shRNA screens conducted on a compendium of cancer cell lines, and also compared several gene summarization methods for inferring gene essentiality from shRNA level data. We then devised novel concepts of seed essentiality and shRNA family, based on seed region sequences of shRNAs, to study in-depth the contribution of seed-mediated off-target effects to the consistency of the two screens. We further investigated two seed-sequence properties, seed pairing stability, and target abundance in terms of their capability to minimize the off-target effects in post-screening data analysis. Finally, we applied this novel methodology to identify genetic interactions and synthetic lethal partners of cancer drivers, and confirmed differential essentiality phenotypes by detailed CRISPR/Cas9 experiments.

RESULTS

Using the novel concepts of seed essentiality and shRNA family, we demonstrate how genome-wide loss-of-function profiling of a common set of cancer cell lines can be actually made fairly reproducible when considering seed-mediated off-target effects. Importantly, by excluding shRNAs having higher propensity for off-target effects, based on their seed-sequence properties, one can remove noise from the genome-wide shRNA datasets. As a translational application case, we demonstrate enhanced reproducibility of genetic interaction partners of common cancer drivers, as well as identify novel synthetic lethal partners of a major oncogenic driver, PIK3CA, supported by a complementary CRISPR/Cas9 experiment.

CONCLUSIONS

We provide practical guidelines for improved design and analysis of genome-wide loss-of-function profiling and demonstrate how this novel strategy can be applied towards improved mapping of genetic dependencies of cancer cells to aid development of targeted anticancer treatments.

Network-based approaches that exploit inferred transcription factor activity to analyze the impact of genetic variation on gene expression

Over the past decade, a number of methods have emerged for inferring protein-level transcription factor activities in individual samples based on prior information about the structure of the gene regulatory network. We discuss how this has enabled new methods for dissecting trans-acting mechanisms that underpin genetic variation in gene expression.

Integrated systems biology analysis of KSHV latent infection reveals viral induction and reliance on peroxisome mediated lipid metabolism

Kaposi’s Sarcoma associated Herpesvirus (KSHV), an oncogenic, human gamma-herpesvirus, is the etiological agent of Kaposi’s Sarcoma the most common tumor of AIDS patients world-wide. KSHV is predominantly latent in the main KS tumor cell, the spindle cell, a cell of endothelial origin. KSHV modulates numerous host cell-signaling pathways to activate endothelial cells including major metabolic pathways involved in lipid metabolism. To identify the underlying cellular mechanisms of KSHV alteration of host signaling and endothelial cell activation, we identified changes in the host proteome, phosphoproteome and transcriptome landscape following KSHV infection of endothelial cells. A Steiner forest algorithm was used to integrate the global data sets and, together with transcriptome based predicted transcription factor activity, cellular networks altered by latent KSHV were predicted. Several interesting pathways were identified, including peroxisome biogenesis. To validate the predictions, we showed that KSHV latent infection increases the number of peroxisomes per cell. Additionally, proteins involved in peroxisomal lipid metabolism of very long chain fatty acids, including ABCD3 and ACOX1, are required for the survival of latently infected cells. In summary, novel cellular pathways altered during herpesvirus latency that could not be predicted by a single systems biology platform, were identified by integrated proteomics and transcriptomics data analysis and when correlated with our metabolomics data revealed that peroxisome lipid metabolism is essential for KSHV latent infection of endothelial cells.

Kaposi’s Sarcoma herpesvirus (KSHV) is the etiologic agent of Kaposi’s Sarcoma, the most common tumor of AIDS patients. KSHV modulates host cell signaling and metabolism to maintain a life-long latent infection. To unravel the underlying cellular mechanisms modulated by KSHV, we used multiple global systems biology platforms to identify and integrate changes in both cellular protein expression and transcription following KSHV infection of endothelial cells, the relevant cell type for KS tumors. The analysis identified several interesting pathways including peroxisome biogenesis. Peroxisomes are small cytoplasmic organelles involved in redox reactions and lipid metabolism. KSHV latent infection increases the number of peroxisomes per cell and proteins involved in peroxisomal lipid metabolism are required for the survival of latently infected cells. In summary, through integration of multiple global systems biology analyses we were able to identify novel pathways that could not be predicted by one platform alone and found that lipid metabolism in a small cytoplasmic organelle is necessary for the survival of latent infection with a herpesvirus.

Systematic Identification and Assessment of Therapeutic Targets for Breast Cancer Based on Genome-Wide RNA Interference Transcriptomes

With accumulating public omics data, great efforts have been made to characterize the genetic heterogeneity of breast cancer. However, identifying novel targets and selecting the best from the sizeable lists of candidate targets is still a key challenge for targeted therapy, largely owing to the lack of economical, efficient and systematic discovery and assessment to prioritize potential therapeutic targets. Here, we describe an approach that combines the computational evaluation and objective, multifaceted assessment to systematically identify and prioritize targets for biological validation and therapeutic exploration. We first establish the reference gene expression profiles from breast cancer cell line MCF7 upon genome-wide RNA interference (RNAi) of a total of 3689 genes, and the breast cancer query signatures using RNA-seq data generated from tissue samples of clinical breast cancer patients in the Cancer Genome Atlas (TCGA). Based on gene set enrichment analysis, we identified a set of 510 genes that when knocked down could significantly reverse the transcriptome of breast cancer state. We then perform multifaceted assessment to analyze the gene set to prioritize potential targets for gene therapy. We also propose drug repurposing opportunities and identify potentially druggable proteins that have been poorly explored with regard to the discovery of small-molecule modulators. Finally, we obtained a small list of candidate therapeutic targets for four major breast cancer subtypes, i.e., luminal A, luminal B, HER2+ and triple negative breast cancer. This RNAi transcriptome-based approach can be a helpful paradigm for relevant researches to identify and prioritize candidate targets for experimental validation.

Pancancer modelling predicts the context-specific impact of somatic mutations on transcriptional programs

Pancancer studies have identified many genes that are frequently somatically altered across multiple tumour types, suggesting that pathway-targeted therapies can be deployed across diverse cancers. However, the same ‘actionable mutation' impacts distinct context-specific gene regulatory programs and signalling networks—and interacts with different genetic backgrounds of co-occurring alterations—in different cancers. Here we apply a computational strategy for integrating parallel (phospho)proteomic and mRNA sequencing data across 12 TCGA tumour data sets to interpret the context-specific impact of somatic alterations in terms of functional signatures such as (phospho)protein and transcription factor (TF) activities. Our analysis predicts distinct dysregulated transcriptional regulators downstream of somatic alterations in different cancers, and we validate the context-specific differential activity of TFs associated to mutant PIK3CA in isogenic cancer cell line models. These results have implications for the pancancer use of targeted drugs and potentially for the design of combination therapies.

Cancer genomic data sets contain a wealth of data that can be used to predict prognosis and further understand disease. Here, the authors integrate multiple genomics data types to identify transcriptional dysregulation in response to somatic mutations.

Transcriptome- and proteome-oriented identification of dysregulated eIF4G, STAT3, and Hippo pathways altered by PIK3CA H1047R in HER2/ER-positive breast cancer

PURPOSE

Phosphatidylinositol 3-kinase (PI3K)/AKT pathway aberrations are common in human breast cancer. Furthermore, PIK3CA mutations are commonly associated with resistance to anti-epidermal growth factor receptor 2 (HER2) or anti-estrogen receptor (ER) agents in HER2 or ER positive (HER2⁺/ER⁺) breast cancer. Hence, deciphering the underlying mechanisms of PIK3CA mutations in HER2⁺/ER⁺ breast cancer would provide novel insights into elucidating resistance to anti-HER2/ER therapies.

METHODS

In this study, we systematically investigated the biological consequences of PIK3CA ^H1047R in HER2⁺/ER⁺ breast cancer by uniquely incorporating mRNA transcriptomic data from The Cancer Genome Atlas and proteomic data from reverse-phase protein arrays.

RESULTS

Our integrative bioinformatics analyses revealed that several important pathways such as STAT3 and VEGF/hypoxia were selectively altered by PIK3CA ^H1047R in HER2⁺/ER⁺ breast cancer. Protein differential expression analysis indicated that an elevated eIF4G might promote tumor angiogenesis and growth via regulation of the hypoxia-activated switch in HER2⁺ PIK3CA ^H1047R breast cancer. We observed hypo-phosphorylation of EGFR in HER2⁺ PIK3CA ^H1047R breast cancer versus HER2⁺PIK3CA^wild-type (PIK3CA ^WT). In addition, ER and PIK3CA ^H1047R might cooperate to activate STAT3, MAPK, AKT, and Hippo pathways in ER⁺ PIK3CA ^H1047R breast cancer. A higher YAP_pS127 level was observed in ER⁺ PIK3CA ^H1047R patients than that in an ER⁺ PIK3CA ^WT subgroup. By examining breast cancer cell lines having both microarray gene expression and drug treatment data from the Genomics of Drug Sensitivity in Cancer and the Stand Up to Cancer datasets, we found that the elevated YAP1 mRNA expression was associated with the resistance of BCL-2 family inhibitors, but with the sensitivity to MEK/MAPK inhibitors in breast cancer cells.

CONCLUSIONS

In summary, these findings shed light on the functional consequences of PIK3CA ^H1047R-driven breast tumorigenesis and resistance to the existing therapeutic agents in HER2⁺/ER⁺ breast cancer.

Integrative proteomic and gene expression analysis identify potential biomarkers for adjuvant trastuzumab resistance: analysis from the Fin-her phase III randomized trial

Trastuzumab is a remarkably effective therapy for patients with human epidermal growth factor receptor 2 (HER2) - positive breast cancer (BC). However, not all women with high levels of HER2 benefit from trastuzumab.

By integrating mRNA and protein expression data from Reverse-Phase Protein Array Analysis (RPPA) in HER2-positive BC, we developed gene expression metagenes that reflect pathway activation levels. Next we assessed the ability of these metagenes to predict resistance to adjuvant trastuzumab using gene expression data from two independent datasets.

10 metagenes passed external validation (false discovery rate [fdr] < 0.05) and showed biological relevance with their pathway of origin. These metagenes were further screened for their association with trastuzumab resistance. An association with trastuzumab resistance was observed and validated only for the AnnexinA1 metagene (ANXA1). In the randomised phase III Fin-her study, tumours with low levels of the ANXA1 metagene showed a benefit from trastuzumab (multivariate: hazard ratio [HR] for distant recurrence = 0.16[95%CI 0.05–0.5]; p = 0.002; fdr = 0.03), while high expression levels of the ANXA1 metagene were associated with a lack of benefit to trastuzmab (HR = 1.29[95%CI 0.55–3.02]; p = 0.56). The association of ANXA1 with trastuzumab resistance was successfully validated in an independent series of subjects who had received trastuzumab with chemotherapy (Log Rank; p = 0.01).

In conclusion, in HER2-positive BC, some proteins are associated with distinct gene expression profiles. Our findings identify the ANXA1metagene as a novel biomarker for trastuzumab resistance.

A Bayesian algorithm for detecting differentially expressed proteins and its application in breast cancer research

Presence of considerable noise and missing data points make analysis of mass-spectrometry (MS) based proteomic data a challenging task. The missing values in MS data are caused by the inability of MS machines to reliably detect proteins whose abundances fall below the detection limit. We developed a Bayesian algorithm that exploits this knowledge and uses missing data points as a complementary source of information to the observed protein intensities in order to find differentially expressed proteins by analysing MS based proteomic data. We compared its accuracy with many other methods using several simulated datasets. It consistently outperformed other methods. We then used it to analyse proteomic screens of a breast cancer (BC) patient cohort. It revealed large differences between the proteomic landscapes of triple negative and Luminal A, which are the most and least aggressive types of BC. Unexpectedly, majority of these differences could be attributed to the direct transcriptional activity of only seven transcription factors some of which are known to be inactive in triple negative BC. We also identified two new proteins which significantly correlated with the survival of BC patients, and therefore may have potential diagnostic/prognostic values.

Integrating Transcriptomic and Proteomic Data Using Predictive Regulatory Network Models of Host Response to Pathogens

Mammalian host response to pathogenic infections is controlled by a complex regulatory network connecting regulatory proteins such as transcription factors and signaling proteins to target genes. An important challenge in infectious disease research is to understand molecular similarities and differences in mammalian host response to diverse sets of pathogens. Recently, systems biology studies have produced rich collections of omic profiles measuring host response to infectious agents such as influenza viruses at multiple levels. To gain a comprehensive understanding of the regulatory network driving host response to multiple infectious agents, we integrated host transcriptomes and proteomes using a network-based approach. Our approach combines expression-based regulatory network inference, structured-sparsity based regression, and network information flow to infer putative physical regulatory programs for expression modules. We applied our approach to identify regulatory networks, modules and subnetworks that drive host response to multiple influenza infections. The inferred regulatory network and modules are significantly enriched for known pathways of immune response and implicate apoptosis, splicing, and interferon signaling processes in the differential response of viral infections of different pathogenicities. We used the learned network to prioritize regulators and study virus and time-point specific networks. RNAi-based knockdown of predicted regulators had significant impact on viral replication and include several previously unknown regulators. Taken together, our integrated analysis identified novel module level patterns that capture strain and pathogenicity-specific patterns of expression and helped identify important regulators of host response to influenza infection.

An important challenge in infectious disease research is to understand how the human immune system responds to different types of pathogenic infections. An important component of mounting proper response is the transcriptional regulatory network that specifies the context-specific gene expression program in the host cell. However, our understanding of this regulatory network and how it drives context-specific transcriptional programs is incomplete. To address this gap, we performed a network-based analysis of host response to influenza viruses that integrated high-throughput mRNA- and protein measurements and protein-protein interaction networks to identify virus and pathogenicity-specific modules and their upstream physical regulatory programs. We inferred regulatory networks for human cell line and mouse host systems, which recapitulated several known regulators and pathways of the immune response and viral life cycle. We used the networks to study time point and strain-specific subnetworks and to prioritize important regulators of host response. We predicted several novel regulators, both at the mRNA and protein levels, and experimentally verified their role in the virus life cycle based on their ability to significantly impact viral replication.

Network-based approaches for analysis of complex biological systems

Cells function and respond to changes in their environment by the coordinated activity of their molecular components, including mRNAs, proteins and metabolites. At the heart of proper cellular function are molecular networks connecting these components to process extra-cellular environmental signals and drive dynamic, context-specific cellular responses. Network-based computational approaches aim to systematically integrate measurements from high-throughput experiments to gain a global understanding of cellular function under changing environmental conditions. We provide an overview of recent methodological developments toward solving two major computational problems within this field in the past two years (2013-2015): network reconstruction and network-based interpretation. Looking forward, we envision development of methods that can predict phenotypes with high accuracy as well as provide biologically plausible mechanistic hypotheses.

An integrative somatic mutation analysis to identify pathways linked with survival outcomes across 19 cancer types

MOTIVATION

Identification of altered pathways that are clinically relevant across human cancers is a key challenge in cancer genomics. Precise identification and understanding of these altered pathways may provide novel insights into patient stratification, therapeutic strategies and the development of new drugs. However, a challenge remains in accurately identifying pathways altered by somatic mutations across human cancers, due to the diverse mutation spectrum. We developed an innovative approach to integrate somatic mutation data with gene networks and pathways, in order to identify pathways altered by somatic mutations across cancers.

RESULTS

We applied our approach to The Cancer Genome Atlas (TCGA) dataset of somatic mutations in 4790 cancer patients with 19 different types of tumors. Our analysis identified cancer-type-specific altered pathways enriched with known cancer-relevant genes and targets of currently available drugs. To investigate the clinical significance of these altered pathways, we performed consensus clustering for patient stratification using member genes in the altered pathways coupled with gene expression datasets from 4870 patients from TCGA, and multiple independent cohorts confirmed that the altered pathways could be used to stratify patients into subgroups with significantly different clinical outcomes. Of particular significance, certain patient subpopulations with poor prognosis were identified because they had specific altered pathways for which there are available targeted therapies. These findings could be used to tailor and intensify therapy in these patients, for whom current therapy is suboptimal.

AVAILABILITY AND IMPLEMENTATION

The code is available at: http://www.taehyunlab.org

CONTACT

jhcheong@yuhs.ac or taehyun.hwang@utsouthwestern.edu or taehyun.cs@gmail.com

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Affinity regression predicts the recognition code of nucleic acid-binding proteins

Predicting the affinity profiles of nucleic acid-binding proteins directly from the protein sequence is a major unsolved problem. We present a statistical approach for learning the recognition code of a family of transcription factors (TFs) or RNA-binding proteins (RBPs) from high-throughput binding assays. Our method, called affinity regression, trains on protein binding microarray (PBM) or RNA compete experiments to learn an interaction model between proteins and nucleic acids, using only protein domain and probe sequences as inputs. By training on mouse homeodomain PBM profiles, our model correctly identifies residues that confer DNA-binding specificity and accurately predicts binding motifs for an independent set of divergent homeodomains. Similarly, learning from RNA compete profiles for diverse RBPs, our model can predict the binding affinities of held-out proteins and identify key RNA-binding residues. More broadly, we envision applying our method to model and predict biological interactions in any setting where there is a high-throughput ‘affinity’ readout.

Targeted Therapies for Triple-Negative Breast Cancer: Combating a Stubborn Disease

Triple-negative breast cancers (TNBCs) constitute a heterogeneous subtype of breast cancers that have a poor clinical outcome. Although no approved targeted therapy is available for TNBCs, molecular-profiling efforts have revealed promising molecular targets, with several candidate compounds having now entered clinical trials for TNBC patients. However, initial results remain modest, thereby highlighting challenges potentially involving intra- and intertumoral heterogeneity and acquisition of therapy resistance. We present a comprehensive review on emerging targeted therapies for treating TNBCs, including the promising approach of immunotherapy and the prognostic value of tumor-infiltrating lymphocytes. We discuss the impact of pathway rewiring in the acquisition of drug resistance, and the prospect of employing combination therapy strategies to overcome challenges towards identifying clinically-viable targeted treatment options for TNBC.

Breast cancer intrinsic subtype classification, clinical use and future trends

Breast cancer is composed of multiple subtypes with distinct morphologies and clinical implications. The advent of microarrays has led to a new paradigm in deciphering breast cancer heterogeneity, based on which the intrinsic subtyping system using prognostic multigene classifiers was developed. Subtypes identified using different gene panels, though overlap to a great extent, do not completely converge, and the avail of new information and perspectives has led to the emergence of novel subtypes, which complicate our understanding towards breast tumor heterogeneity. This review explores and summarizes the existing intrinsic subtypes, patient clinical features and management, commercial signature panels, as well as various information used for tumor classification. Two trends are pointed out in the end on breast cancer subtyping, i.e., either diverging to more refined groups or converging to the major subtypes. This review improves our understandings towards breast cancer intrinsic classification, current status on clinical application, and future trends.

Constitutive phosphorylated STAT3-associated gene signature is predictive for trastuzumab resistance in primary HER2-positive breast cancer

BACKGROUND

The likelihood of recurrence in patients with breast cancer who have HER2-positive tumors is relatively high, although trastuzumab is a remarkably effective drug in this setting. Signal transducer and activator of transcription 3 protein (STAT3), a transcription factor that is persistently tyrosine-705 phosphorylated (pSTAT3) in response to numerous oncogenic signaling pathways, activates downstream proliferative and anti-apoptotic pathways. We hypothesized that pSTAT3 expression in HER2-positive breast cancer will confer trastuzumab resistance.

METHODS

We integrated reverse phase protein array (RPPA) and gene expression data from patients with HER2-positive breast cancer treated with trastuzumab in the adjuvant setting.

RESULTS

We show that a pSTAT3-associated gene signature (pSTAT3-GS) is able to predict pSTAT3 status in an independent dataset (TCGA; AUC = 0.77, P = 0.02). This suggests that STAT3 induces a characteristic set of gene expression changes in HER2-positive cancers. Tumors characterized as high pSTAT3-GS were associated with trastuzumab resistance (log rank P = 0.049). These results were confirmed using data from the prospective, randomized controlled FinHer study, where the effect was especially prominent in HER2-positive estrogen receptor (ER)-negative tumors (interaction test P = 0.02). Of interest, constitutively activated pSTAT3 tumors were associated with loss of PTEN, elevated IL6, and stromal reactivation.

CONCLUSIONS

This study provides compelling evidence for a link between pSTAT3 and trastuzumab resistance in HER2-positive primary breast cancers. Our results suggest that it may be valuable to add agents targeting the STAT3 pathway to trastuzumab for treatment of HER2-positive breast cancer.

Integrated network analysis and logistic regression modeling identify stage-specific genes in Oral Squamous Cell Carcinoma

BACKGROUND

Oral squamous cell carcinoma (OSCC) is associated with substantial mortality and morbidity but, OSCC can be difficult to detect at its earliest stage due to its molecular complexity and clinical behavior. Therefore, identification of key gene signatures at an early stage will be highly helpful.

METHODS

The aim of this study was to identify key genes associated with progression of OSCC stages. Gene expression profiles were classified into cancer stage-related modules, i.e., groups of genes that are significantly related to a clinical stage. For prioritizing the candidate genes, analysis was further restricted to genes with high connectivity and a significant association with a stage. To assess predictive power of these genes, a classification model was also developed and tested by 5-fold cross validation and on an independent dataset.

RESULTS

The identified genes were enriched for significant processes and functional pathways, and various genes were found to be directly implicated in OSCC. Forward and stepwise, multivariate logistic regression analyses identified 13 key genes whose expression discriminated early- and late-stage OSCC with predictive accuracy (area under curve; AUC) of ~0.81 in a 5-fold cross-validation strategy.

CONCLUSIONS

The proposed network-driven integrative analytical approach can identify multiple genes significantly related to an OSCC stage; the classification model that is developed with these genes may help to distinguish cancer stages. The proposed genes and model hold promise for monitoring of OSCC stage progression, and our findings may facilitate cancer detection at an earlier stage, resulting in improved treatment outcomes.