Quantitative modeling of signaling in aggressive B cell lymphoma unveils conserved core network.

B cell receptor (BCR) signaling is required for the survival and maturation of B cells and is deregulated in B cell lymphomas. While proximal BCR signaling is well studied, little is known about the crosstalk of downstream effector pathways, and a comprehensive quantitative network analysis of BCR signaling is missing. Here, we semi-quantitatively modelled BCR signaling in Burkitt lymphoma (BL) cells using systematically perturbed phosphorylation data of BL-2 and BL-41 cells. The models unveiled feedback and crosstalk structures in the BCR signaling network, including a negative crosstalk from p38 to MEK/ERK. The relevance of the crosstalk was verified for BCR and CD40 signaling in different BL cells and confirmed by global phosphoproteomics on ERK itself and known ERK target sites. Compared to the starting network, the trained network for BL-2 cells was better transferable to BL-41 cells. Moreover, the BL-2 network was also suited to model BCR signaling in Diffuse large B cell lymphoma cells lines with aberrant BCR signaling (HBL-1, OCI-LY3), indicating that BCR aberration does not cause a major downstream rewiring.

Modeling unveils sex differences of signaling networks in mouse embryonic stem cells.

For a short period during early development of mammalian embryos, both X chromosomes in females are active, before dosage compensation is ensured through X-chromosome inactivation. In female mouse embryonic stem cells (mESCs), which carry two active X chromosomes, increased X-dosage affects cell signaling and impairs differentiation. The underlying mechanisms, however, remain poorly understood. To dissect X-dosage effects on the signaling network in mESCs, we combine systematic perturbation experiments with mathematical modeling. We quantify the response to a variety of inhibitors and growth factors for cells with one (XO) or two X chromosomes (XX). We then build models of the signaling networks in XX and XO cells through a semi-quantitative modeling approach based on modular response analysis. We identify a novel negative feedback in the PI3K/AKT pathway through GSK3. Moreover, the presence of a single active X makes mESCs more sensitive to the differentiation-promoting Activin A signal and leads to a stronger RAF1-mediated negative feedback in the FGF-triggered MAPK pathway. The differential response to these differentiation-promoting pathways can explain the impaired differentiation propensity of female mESCs.

Neuroblastoma signalling models unveil combination therapies targeting feedback-mediated resistance.

Very high risk neuroblastoma is characterised by increased MAPK signalling, and targeting MAPK signalling is a promising therapeutic strategy. We used a deeply characterised panel of neuroblastoma cell lines and found that the sensitivity to MEK inhibitors varied drastically between these cell lines. By generating quantitative perturbation data and mathematical modelling, we determined potential resistance mechanisms. We found that negative feedbacks within MAPK signalling and via the IGF receptor mediate re-activation of MAPK signalling upon treatment in resistant cell lines. By using cell-line specific models, we predict that combinations of MEK inhibitors with RAF or IGFR inhibitors can overcome resistance, and tested these predictions experimentally. In addition, phospho-proteomic profiling confirmed the cell-specific feedback effects and synergy of MEK and IGFR targeted treatment. Our study shows that a quantitative understanding of signalling and feedback mechanisms facilitated by models can help to develop and optimise therapeutic strategies. Our findings should be considered for the planning of future clinical trials introducing MEKi in the treatment of neuroblastoma.

SPEED2: inferring upstream pathway activity from differential gene expression.

Extracting signalling pathway activities from transcriptome data is important to infer mechanistic origins of transcriptomic dysregulation, for example in disease. A popular method to do so is by enrichment analysis of signature genes in e.g. differentially regulated genes. Previously, we derived signatures for signalling pathways by integrating public perturbation transcriptome data and generated a signature database called SPEED (Signalling Pathway Enrichment using Experimental Datasets), for which we here present a substantial upgrade as SPEED2. This web server hosts consensus signatures for 16 signalling pathways that are derived from a large number of transcriptomic signalling perturbation experiments. When providing a gene list of e.g. differentially expressed genes, the web server allows to infer signalling pathways that likely caused these genes to be deregulated. In addition to signature lists, we derive 'continuous' gene signatures, in a transparent and automated fashion without any fine-tuning, and describe a new algorithm to score these signatures.

Comparative Network Reconstruction using mixed integer programming.

Motivation: Signal-transduction networks are often aberrated in cancer cells, and new anti-cancer drugs that specifically target oncogenes involved in signaling show great clinical promise. However, the effectiveness of such targeted treatments is often hampered by innate or acquired resistance due to feedbacks, crosstalks or network adaptations in response to drug treatment. A quantitative understanding of these signaling networks and how they differ between cells with different oncogenic mutations or between sensitive and resistant cells can help in addressing this problem. Results: Here, we present Comparative Network Reconstruction (CNR), a computational method to reconstruct signaling networks based on possibly incomplete perturbation data, and to identify which edges differ quantitatively between two or more signaling networks. Prior knowledge about network topology is not required but can straightforwardly be incorporated. We extensively tested our approach using simulated data and applied it to perturbation data from a BRAF mutant, PTPN11 KO cell line that developed resistance to BRAF inhibition. Comparing the reconstructed networks of sensitive and resistant cells suggests that the resistance mechanism involves re-establishing wild-type MAPK signaling, possibly through an alternative RAF-isoform. Availability and implementation: CNR is available as a python module at https://github.com/NKI-CCB/cnr. Additionally, code to reproduce all figures is available at https://github.com/NKI-CCB/CNR-analyses. Supplementary information: Supplementary data are available at Bioinformatics online.

Modelling signalling networks from perturbation data.

Motivation: Intracellular signalling is realized by complex signalling networks, which are almost impossible to understand without network models, especially if feedbacks are involved. Modular Response Analysis (MRA) is a convenient modelling method to study signalling networks in various contexts. Results: We developed the software package STASNet (STeady-STate Analysis of Signalling Networks) that provides an augmented and extended version of MRA suited to model signalling networks from incomplete perturbation schemes and multi-perturbation data. Using data from the Dialogue on Reverse Engineering Assessment and Methods challenge, we show that predictions from STASNet models are among the top-performing methods. We applied the method to study the effect of SHP2, a protein that has been implicated in resistance to targeted therapy in colon cancer, using a novel dataset from the colon cancer cell line Widr and a SHP2-depleted derivative. We find that SHP2 is required for mitogen-activated protein kinase signalling, whereas AKT signalling only partially depends on SHP2. Availability and implementation: An R-package is available at https://github.com/molsysbio/STASNet. Supplementary information: Supplementary data are available at Bioinformatics online.

An immediate-late gene expression module decodes ERK signal duration.

The RAF-MEK-ERK signalling pathway controls fundamental, often opposing cellular processes such as proliferation and apoptosis. Signal duration has been identified to play a decisive role in these cell fate decisions. However, it remains unclear how the different early and late responding gene expression modules can discriminate short and long signals. We obtained both protein phosphorylation and gene expression time course data from HEK293 cells carrying an inducible construct of the proto-oncogene RAF By mathematical modelling, we identified a new gene expression module of immediate-late genes (ILGs) distinct in gene expression dynamics and function. We find that mRNA longevity enables these ILGs to respond late and thus translate ERK signal duration into response amplitude. Despite their late response, their GC-rich promoter structure suggested and metabolic labelling with 4SU confirmed that transcription of ILGs is induced immediately. A comparative analysis shows that the principle of duration decoding is conserved in PC12 cells and MCF7 cells, two paradigm cell systems for ERK signal duration. Altogether, our findings suggest that ILGs function as a gene expression module to decode ERK signal duration.

Hypoxia-induced gene expression results from selective mRNA partitioning to the endoplasmic reticulum.

Protein synthesis is a primary energy-consuming process in the cell. Therefore, under hypoxic conditions, rapid inhibition of global mRNA translation represents a major protective strategy to maintain energy metabolism. How some mRNAs, especially those that encode crucial survival factors, continue to be efficiently translated in hypoxia is not completely understood. By comparing specific transcript levels in ribonucleoprotein complexes, cytoplasmic polysomes and endoplasmic reticulum (ER)-bound ribosomes, we show that the synthesis of proteins encoded by hypoxia marker genes is favoured at the ER in hypoxia. Gene expression profiling revealed that transcripts particularly increased by the HIF-1 transcription factor network show hypoxia-induced enrichment at the ER. We found that mRNAs favourably translated at the ER have higher conservation scores for both the 5'- and 3'-untranslated regions (UTRs) and contain less upstream initiation codons (uAUGs), indicating the significance of these sequence elements for sustained mRNA translation under hypoxic conditions. Furthermore, we found enrichment of specific cis-elements in mRNA 5'- as well as 3'-UTRs that mediate transcript localization to the ER in hypoxia. We conclude that transcriptome partitioning between the cytoplasm and the ER permits selective mRNA translation under conditions of energy shortage.

Computer-assisted curation of a human regulatory core network from the biological literature.

MOTIVATION: A highly interlinked network of transcription factors (TFs) orchestrates the context-dependent expression of human genes. ChIP-chip experiments that interrogate the binding of particular TFs to genomic regions are used to reconstruct gene regulatory networks at genome-scale, but are plagued by high false-positive rates. Meanwhile, a large body of knowledge on high-quality regulatory interactions remains largely unexplored, as it is available only in natural language descriptions scattered over millions of scientific publications. Such data are hard to extract and regulatory data currently contain together only 503 regulatory relations between human TFs. RESULTS: We developed a text-mining-assisted workflow to systematically extract knowledge about regulatory interactions between human TFs from the biological literature. We applied this workflow to the entire Medline, which helped us to identify more than 45 000 sentences potentially describing such relationships. We ranked these sentences by a machine-learning approach. The top-2500 sentences contained ∼900 sentences that encompass relations already known in databases. By manually curating the remaining 1625 top-ranking sentences, we obtained more than 300 validated regulatory relationships that were not present in a regulatory database before. Full-text curation allowed us to obtain detailed information on the strength of experimental evidences supporting a relationship. CONCLUSIONS: We were able to increase curated information about the human core transcriptional network by >60% compared with the current content of regulatory databases. We observed improved performance when using the network for disease gene prioritization compared with the state-of-the-art. AVAILABILITY AND IMPLEMENTATION: Web-service is freely accessible at http://fastforward.sys-bio.net/. CONTACT: leser@informatik.hu-berlin.de or nils.bluethgen@charite.de SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Consequences of feedback in signal transduction for targeted therapies.

Over the last two decades, many small-molecule inhibitors that target kinase signalling have been developed. More than 20 of these inhibitors are FDA (U.S. Food and Drug Administration)-approved and are now being used in the clinics to treat tumours; even more have entered clinical trials. However, resistance to these inhibitors, either intrinsic to the tumour or acquired during treatment, remains a major problem in targeted therapeutics. One of the mechanisms by which tumours become resistant is the rewiring of the signalling networks via feedback, by which the tumour cells re-activate signalling or activate alternative signalling pathways. In the present article, we review insights from recent quantitative signalling studies combining mathematical modelling and experiments that revealed how feedback rewires MAPK (mitogen-activated protein kinase)/PI3K (phosphoinositide 3-kinase) signalling upon treatment and how that affects drug sensitivity.

Network quantification of EGFR signaling unveils potential for targeted combination therapy.

The epidermal growth factor receptor (EGFR) signaling network is activated in most solid tumors, and small-molecule drugs targeting this network are increasingly available. However, often only specific combinations of inhibitors are effective. Therefore, the prediction of potent combinatorial treatments is a major challenge in targeted cancer therapy. In this study, we demonstrate how a model-based evaluation of signaling data can assist in finding the most suitable treatment combination. We generated a perturbation data set by monitoring the response of RAS/PI3K signaling to combined stimulations and inhibitions in a panel of colorectal cancer cell lines, which we analyzed using mathematical models. We detected that a negative feedback involving EGFR mediates strong cross talk from ERK to AKT. Consequently, when inhibiting MAPK, AKT activity is increased in an EGFR-dependent manner. Using the model, we predict that in contrast to single inhibition, combined inactivation of MEK and EGFR could inactivate both endpoints of RAS, ERK and AKT. We further could demonstrate that this combination blocked cell growth in BRAF- as well as KRAS-mutated tumor cells, which we confirmed using a xenograft model.

Reverse engineering a hierarchical regulatory network downstream of oncogenic KRAS.

RAS mutations are highly relevant for progression and therapy response of human tumours, but the genetic network that ultimately executes the oncogenic effects is poorly understood. Here, we used a reverse-engineering approach in an ovarian cancer model to reconstruct KRAS oncogene-dependent cytoplasmic and transcriptional networks from perturbation experiments based on gene silencing and pathway inhibitor treatments. We measured mRNA and protein levels in manipulated cells by microarray, RT-PCR and western blot analysis, respectively. The reconstructed model revealed complex interactions among the transcriptional and cytoplasmic components, some of which were confirmed by double pertubation experiments. Interestingly, the transcription factors decomposed into two hierarchically arranged groups. To validate the model predictions, we analysed growth parameters and transcriptional deregulation in the KRAS-transformed epithelial cells. As predicted by the model, we found two functional groups among the selected transcription factors. The experiments thus confirmed the predicted hierarchical transcription factor regulation and showed that the hierarchy manifests itself in downstream gene expression patterns and phenotype.

Combination of tumor asphericity and an extracellular matrix-related prognostic gene signature in non-small cell lung cancer patients.

One important aim of precision oncology is a personalized treatment of patients. This can be achieved by various biomarkers, especially imaging parameters and gene expression signatures are commonly used. So far, combination approaches are sparse. The aim of the study was to independently validate the prognostic value of the novel positron emission tomography (PET) parameter tumor asphericity (ASP) in non small cell lung cancer (NSCLC) patients and to investigate associations between published gene expression profiles and ASP. This was a retrospective evaluation of PET imaging and gene expression data from three public databases and two institutional datasets. The whole cohort comprised 253 NSCLC patients, all treated with curative intent surgery. Clinical parameters, standard PET parameters and ASP were evaluated in all patients. Additional gene expression data were available for 120 patients. Univariate Cox regression and Kaplan-Meier analysis was performed for the primary endpoint progression-free survival (PFS) and additional endpoints. Furthermore, multivariate cox regression testing was performed including clinically significant parameters, ASP, and the extracellular matrix-related prognostic gene signature (EPPI). In the whole cohort, a significant association with PFS was observed for ASP (p < 0.001) and EPPI (p = 0.012). Upon multivariate testing, EPPI remained significantly associated with PFS (p = 0.018) in the subgroup of patients with additional gene expression data, while ASP was significantly associated with PFS in the whole cohort (p = 0.012). In stage II patients, ASP was significantly associated with PFS (p = 0.009), and a previously published cutoff value for ASP (19.5%) was successfully validated (p = 0.008). In patients with additional gene expression data, EPPI showed a significant association with PFS, too (p = 0.033). The exploratory combination of ASP and EPPI showed that the combinatory approach has potential to further improve patient stratification compared to the use of only one parameter. We report the first successful validation of EPPI and ASP in stage II NSCLC patients. The combination of both parameters seems to be a very promising approach for improvement of risk stratification in a group of patients with urgent need for a more personalized treatment approach.

Discovering causal signaling pathways through gene-expression patterns.

High-throughput gene-expression studies result in lists of differentially expressed genes. Most current meta-analyses of these gene lists include searching for significant membership of the translated proteins in various signaling pathways. However, such membership enrichment algorithms do not provide insight into which pathways caused the genes to be differentially expressed in the first place. Here, we present an intuitive approach for discovering upstream signaling pathways responsible for regulating these differentially expressed genes. We identify consistently regulated signature genes specific for signal transduction pathways from a panel of single-pathway perturbation experiments. An algorithm that detects overrepresentation of these signature genes in a gene group of interest is used to infer the signaling pathway responsible for regulation. We expose our novel resource and algorithm through a web server called SPEED: Signaling Pathway Enrichment using Experimental Data sets. SPEED can be freely accessed at http://speed.sys-bio.net/.

Mutation-specific effects of NRAS oncogenes in colorectal cancer cells.

In colorectal cancer (CRC), the prevalence of NRAS mutations (5-9%) is inferior to that of KRAS mutations (40-50%). NRAS mutations feature lately during tumour progression and drive resistance to anti-EGFR therapy in KRAS wild-type tumours. To elucidate specific functions of NRAS mutations in CRC, we expressed doxycycline-inducible G12D and Q61K mutations in the CRC cell line Caco-2. A focused phospho-proteome analysis based on the Bio-Plex platform, which interrogated the activity of MAPK, PI3K, mTOR, STAT, p38, JNK and ATF2, did not reveal significant differences between Caco-2 cells expressing NRASG12D, NRASQ61K and KRASG12V. However, phenotypic read-outs were different. The NRAS Q61K mutation promoted anchorage-independent proliferation and tumorigenicity, similar to features driven by canonical KRAS mutations. In contrast, expression of NRASG12D resulted in reduced proliferation and apoptosis. At the transcriptome level, we saw upregulation of cytokines and chemokines. IL1A, IL11, CXCL8 (IL-8) and CCL20 exhibited enhanced secretion into the culture medium. In addition, RNA sequencing results indicated activation of the IL1-, JAK/STAT-, NFκB- and TNFα signalling pathways. These results form the basis for an NRASG12D-driven inflammatory phenotype in CRC.

Early senescence and production of senescence-associated cytokines are major determinants of radioresistance in head-and-neck squamous cell carcinoma.

Resistance against radio(chemo)therapy-induced cell death is a major determinant of oncological treatment failure and remains a perpetual clinical challenge. The underlying mechanisms are manifold and demand for comprehensive, cancer entity- and subtype-specific examination. In the present study, resistance against radiotherapy was systematically assessed in a panel of human head-and-neck squamous cell carcinoma (HNSCC) cell lines and xenotransplants derived thereof with the overarching aim to extract master regulators and potential candidates for mechanism-based pharmacological targeting. Clonogenic survival data were integrated with molecular and functional data on DNA damage repair and different cell fate decisions. A positive correlation between radioresistance and early induction of HNSCC cell senescence accompanied by NF-κB-dependent production of distinct senescence-associated cytokines, particularly ligands of the CXCR2 chemokine receptor, was identified. Time-lapse microscopy and medium transfer experiments disclosed the non-cell autonomous, paracrine nature of these mechanisms, and pharmacological interference with senescence-associated cytokine production by the NF-κB inhibitor metformin significantly improved radiotherapeutic performance in vitro and in vivo. With regard to clinical relevance, retrospective analyses of TCGA HNSCC data and an in-house HNSCC cohort revealed that elevated expression of CXCR2 and/or its ligands are associated with impaired treatment outcome. Collectively, our study identifies radiation-induced tumor cell senescence and the NF-κB-dependent production of distinct senescence-associated cytokines as critical drivers of radioresistance in HNSCC whose therapeutic targeting in the context of multi-modality treatment approaches should be further examined and may be of particular interest for the subgroup of patients with elevated expression of the CXCR2/ligand axis.

Neuronal activity regulates alternative exon usage.

Neuronal activity-regulated gene transcription underlies plasticity-dependent changes in the molecular composition and structure of neurons. A large number of genes regulated by different neuronal plasticity inducing pathways have been identified, but altered gene expression levels represent only part of the complexity of the activity-regulated transcriptional program. Alternative splicing, the differential inclusion and exclusion of exonic sequence in mRNA, is an additional mechanism that is thought to define the activity-dependent transcriptome. Here, we present a genome wide microarray-based survey to identify exons with increased expression levels at 1, 4 or 8 h following neuronal activity in the murine hippocampus provoked by generalized seizures. We used two different bioinformatics approaches to identify alternative activity-induced exon usage and to predict alternative splicing, ANOSVA (ANalysis Of Splicing VAriation) which we here adjusted to accommodate data from different time points and FIRMA (Finding Isoforms using Robust Multichip Analysis). RNA sequencing, in situ hybridization and reverse transcription PCR validate selected activity-dependent splicing events of previously described and so far undescribed activity-regulated transcripts, including Homer1a, Homer1d, Ania3, Errfi1, Inhba, Dclk1, Rcan1, Cda, Tpm1 and Krt75. Taken together, our survey significantly adds to the comprehensive understanding of the complex activity-dependent neuronal transcriptomic signature. In addition, we provide data sets that will serve as rich resources for future comparative expression analyses.

Reduced replication origin licensing selectively kills KRAS-mutant colorectal cancer cells via mitotic catastrophe.

To unravel vulnerabilities of KRAS-mutant CRC cells, a shRNA-based screen specifically inhibiting MAPK pathway components and targets was performed in CaCo2 cells harboring conditional oncogenic KRASG12V. The custom-designed shRNA library comprised 121 selected genes, which were previously identified to be strongly regulated in response to MEK inhibition. The screen showed that CaCo2 cells expressing KRASG12V were sensitive to the suppression of the DNA replication licensing factor minichromosome maintenance complex component 7 (MCM7), whereas KRASwt CaCo2 cells were largely resistant to MCM7 suppression. Similar results were obtained in an isogenic DLD-1 cell culture model. Knockdown of MCM7 in a KRAS-mutant background led to replication stress as indicated by increased nuclear RPA focalization. Further investigation showed a significant increase in mitotic cells after simultaneous MCM7 knockdown and KRASG12V expression. The increased percentage of mitotic cells coincided with strongly increased DNA damage in mitosis. Taken together, the accumulation of DNA damage in mitotic cells is due to replication stress that remained unresolved, which results in mitotic catastrophe and cell death. In summary, the data show a vulnerability of KRAS-mutant cells towards suppression of MCM7 and suggest that inhibiting DNA replication licensing might be a viable strategy to target KRAS-mutant cancers.

SFPQ Depletion Is Synthetically Lethal with BRAFV600E in Colorectal Cancer Cells.

Oncoproteins such as the BRAFV600E kinase endow cancer cells with malignant properties, but they also create unique vulnerabilities. Targeting of BRAFV600E-driven cytoplasmic signaling networks has proved ineffective, as patients regularly relapse with reactivation of the targeted pathways. We identify the nuclear protein SFPQ to be synthetically lethal with BRAFV600E in a loss-of-function shRNA screen. SFPQ depletion decreases proliferation and specifically induces S-phase arrest and apoptosis in BRAFV600E-driven colorectal and melanoma cells. Mechanistically, SFPQ loss in BRAF-mutant cancer cells triggers the Chk1-dependent replication checkpoint, results in decreased numbers and reduced activities of replication factories, and increases collision between replication and transcription. We find that BRAFV600E-mutant cancer cells and organoids are sensitive to combinations of Chk1 inhibitors and chemically induced replication stress, pointing toward future therapeutic approaches exploiting nuclear vulnerabilities induced by BRAFV600E.