Optimizing network propagation for multi-omics data integration.

Network propagation refers to a class of algorithms that integrate information from input data across connected nodes in a given network. These algorithms have wide applications in systems biology, protein function prediction, inferring condition-specifically altered sub-networks, and prioritizing disease genes. Despite the popularity of network propagation, there is a lack of comparative analyses of different algorithms on real data and little guidance on how to select and parameterize the various algorithms. Here, we address this problem by analyzing different combinations of network normalization and propagation methods and by demonstrating schemes for the identification of optimal parameter settings on real proteome and transcriptome data. Our work highlights the risk of a 'topology bias' caused by the incorrect use of network normalization approaches. Capitalizing on the fact that network propagation is a regularization approach, we show that minimizing the bias-variance tradeoff can be utilized for selecting optimal parameters. The application to real multi-omics data demonstrated that optimal parameters could also be obtained by either maximizing the agreement between different omics layers (e.g. proteome and transcriptome) or by maximizing the consistency between biological replicates. Furthermore, we exemplified the utility and robustness of network propagation on multi-omics datasets for identifying ageing-associated genes in brain and liver tissues of rats and for elucidating molecular mechanisms underlying prostate cancer progression. Overall, this work compares different network propagation approaches and it presents strategies for how to use network propagation algorithms to optimally address a specific research question at hand.

Deciphering the signaling network of breast cancer improves drug sensitivity prediction.

One goal of precision medicine is to tailor effective treatments to patients' specific molecular markers of disease. Here, we used mass cytometry to characterize the single-cell signaling landscapes of 62 breast cancer cell lines and five lines from healthy tissue. We quantified 34 markers in each cell line upon stimulation by the growth factor EGF in the presence or absence of five kinase inhibitors. These data-on more than 80 million single cells from 4,000 conditions-were used to fit mechanistic signaling network models that provide insight into how cancer cells process information. Our dynamic single-cell-based models accurately predicted drug sensitivity and identified genomic features associated with drug sensitivity, including a missense mutation in DDIT3 predictive of PI3K-inhibition sensitivity. We observed similar trends in genotype-drug sensitivity associations in patient-derived xenograft mouse models. This work provides proof of principle that patient-specific single-cell measurements and modeling could inform effective precision medicine strategies.

Convergent network effects along the axis of gene expression during prostate cancer progression.

BACKGROUND: Tumor-specific genomic aberrations are routinely determined by high-throughput genomic measurements. It remains unclear how complex genome alterations affect molecular networks through changing protein levels and consequently biochemical states of tumor tissues. RESULTS: Here, we investigate the propagation of genomic effects along the axis of gene expression during prostate cancer progression. We quantify genomic, transcriptomic, and proteomic alterations based on 105 prostate samples, consisting of benign prostatic hyperplasia regions and malignant tumors, from 39 prostate cancer patients. Our analysis reveals the convergent effects of distinct copy number alterations impacting on common downstream proteins, which are important for establishing the tumor phenotype. We devise a network-based approach that integrates perturbations across different molecular layers, which identifies a sub-network consisting of nine genes whose joint activity positively correlates with increasingly aggressive tumor phenotypes and is predictive of recurrence-free survival. Further, our data reveal a wide spectrum of intra-patient network effects, ranging from similar to very distinct alterations on different molecular layers. CONCLUSIONS: This study uncovers molecular networks with considerable convergent alterations across tumor sites and patients. It also exposes a diversity of network effects: we could not identify a single sub-network that is perturbed in all high-grade tumor regions.

Proteomics-Based Monitoring of Pathway Activity Reveals that Blocking Diacylglycerol Biosynthesis Rescues from Alpha-Synuclein Toxicity.

Proteinaceous inclusions containing alpha-synuclein (α-Syn) have been implicated in neuronal toxicity in Parkinson's disease, but the pathways that modulate toxicity remain enigmatic. Here, we used a targeted proteomic assay to simultaneously measure 269 pathway activation markers and proteins deregulated by α-Syn expression across a panel of 33 Saccharomyces cerevisiae strains that genetically modulate α-Syn toxicity. Applying multidimensional linear regression analysis to these data predicted Pah1, a phosphatase that catalyzes conversion of phosphatidic acid to diacylglycerol at the endoplasmic reticulum membrane, as an effector of rescue. Follow-up studies demonstrated that inhibition of Pah1 activity ameliorates the toxic effects of α-Syn, indicate that the diacylglycerol branch of lipid metabolism could enhance α-Syn neuronal cytotoxicity, and suggest a link between α-Syn toxicity and the biology of lipid droplets.

Comparative analysis of mRNA and protein degradation in prostate tissues indicates high stability of proteins.

Deterioration of biomolecules in clinical tissues is an inevitable pre-analytical process, which affects molecular measurements and thus potentially confounds conclusions from cohort analyses. Here, we investigate the degradation of mRNA and protein in 68 pairs of adjacent prostate tissue samples using RNA-Seq and SWATH mass spectrometry, respectively. To objectively quantify the extent of protein degradation, we develop a numerical score, the Proteome Integrity Number (PIN), that faithfully measures the degree of protein degradation. Our results indicate that protein degradation only affects 5.9% of the samples tested and shows negligible correlation with mRNA degradation in the adjacent samples. These findings are confirmed by independent analyses on additional clinical sample cohorts and across different mass spectrometric methods. Overall, the data show that the majority of samples tested are not compromised by protein degradation, and establish the PIN score as a generic and accurate indicator of sample quality for proteomic analyses.

Multi-region proteome analysis quantifies spatial heterogeneity of prostate tissue biomarkers.

It remains unclear to what extent tumor heterogeneity impacts on protein biomarker discovery. Here, we quantified proteome intra-tissue heterogeneity (ITH) based on a multi-region analysis of prostate tissues using pressure cycling technology and SWATH mass spectrometry. We quantified 6,873 proteins and analyzed the ITH of 3,700 proteins. The level of ITH varied depending on proteins and tissue types. Benign tissues exhibited more complex ITH patterns than malignant tissues. Spatial variability of ten prostate biomarkers was validated by immunohistochemistry in an independent cohort (n=83) using tissue microarrays. PSA was preferentially variable in benign prostatic hyperplasia, while GDF15 substantially varied in prostate adenocarcinomas. Further, we found that DNA repair pathways exhibited a high degree of variability in tumorous tissues, which may contribute to the genetic heterogeneity of tumors. This study conceptually adds a new perspective to protein biomarker discovery: it suggests that recent technological progress should be exploited to quantify and account for spatial proteome variation to complement biomarker identification and utilization.

Several immune escape patterns in non-Hodgkin's lymphomas.

Follicular Lymphomas (FL) and diffuse large B cell lymphomas (DLBCL) must evolve some immune escape strategy to develop from lymphoid organs, but their immune evasion pathways remain poorly characterized. We investigated this issue by transcriptome data mining and immunohistochemistry (IHC) of FL and DLBCL lymphoma biopsies. A set of genes involved in cancer immune-evasion pathways (Immune Escape Gene Set, IEGS) was defined and the distribution of the expression levels of these genes was compared in FL, DLBCL and normal B cell transcriptomes downloaded from the GEO database. The whole IEGS was significantly upregulated in all the lymphoma samples but not in B cells or other control tissues, as shown by the overexpression of the PD-1, PD-L1, PD-L2 and LAG3 genes. Tissue microarray immunostainings for PD-1, PD-L1, PD-L2 and LAG3 proteins on additional biopsies from 27 FL and 27 DLBCL patients confirmed the expression of these proteins. The immune infiltrates were more abundant in FL than DLBCL samples, and the microenvironment of FL comprised higher rates of PD-1+ lymphocytes. Further, DLBCL tumor cells comprised a higher proportion of PD-1+, PD-L1+, PD-L2+ and LAG3+ lymphoma cells than the FL tumor cells, confirming that DLBCL mount immune escape strategies distinct from FL. In addition, some cases of DLBCL had tumor cells co-expressing both PD-1, PD-L1 and PD-L2. Among the DLBCLs, the activated B cell (ABC) subtype comprised more PD-L1+ and PD-L2+ lymphoma cells than the GC subtype. Thus, we infer that FL and DLBCL evolved several pathways of immune escape.

Weighted Kolmogorov Smirnov testing: an alternative for Gene Set Enrichment Analysis.

Gene Set Enrichment Analysis (GSEA) is a basic tool for genomic data treatment. Its test statistic is based on a cumulated weight function, and its distribution under the null hypothesis is evaluated by Monte-Carlo simulation. Here, it is proposed to subtract to the cumulated weight function its asymptotic expectation, then scale it. Under the null hypothesis, the convergence in distribution of the new test statistic is proved, using the theory of empirical processes. The limiting distribution needs to be computed only once, and can then be used for many different gene sets. This results in large savings in computing time. The test defined in this way has been called Weighted Kolmogorov Smirnov (WKS) test. Using expression data from the GEO repository, tested against the MSig Database C2, a comparison between the classical GSEA test and the new procedure has been conducted. Our conclusion is that, beyond its mathematical and algorithmic advantages, the WKS test could be more informative in many cases, than the classical GSEA test.

A six gene expression signature defines aggressive subtypes and predicts outcome in childhood and adult acute lymphoblastic leukemia.

Abnormal gene expression in cancer represents an under-explored source of cancer markers and therapeutic targets. In order to identify gene expression signatures associated with survival in acute lymphoblastic leukemia (ALL), a strategy was designed to search for aberrant gene activity, which consists of applying several filters to transcriptomic datasets from two pediatric ALL studies. Six genes whose expression in leukemic blasts was associated with prognosis were identified:three genes predicting poor prognosis (AK022211, FASTKD1 and STARD4) and three genes associated with a favorable outcome (CAMSAP1, PCGF6 and SH3RF3). Combining the expression of these 6 genes could successfully predict prognosis not only in the two discovery pediatric ALL studies, but also in two independent validation cohorts of adult patients, one from a publicly available study and one consisting of 62 newly recruited Chinese patients. Moreover, our data demonstrate that our six gene based test is particularly efficient in stratifying MLL or BCR.ABL negative patients. Finally, common biological traits characterizing aggressive forms of ALL in both children and adults were found, including features of dormant hematopoietic stem cells, suggesting new therapeutic strategies.