Integrative network analysis for survival-associated gene-gene interactions across multiple genomic profiles in ovarian cancer

A Drug Repurposing Pipeline Based on Bladder Cancer Integrated Proteotranscriptomics Signatures

Delivering better care for patients with bladder cancer (BC) necessitates the development of novel therapeutic strategies that address both the high disease heterogeneity and the limitations of the current therapeutic modalities, such as drug low efficacy and patient resistance acquisition. Drug repurposing is a cost-effective strategy that targets the reuse of existing drugs for new therapeutic purposes. Such a strategy could open new avenues toward more effective BC treatment. BC patients' multi-omics signatures can be used to guide the investigation of existing drugs that show an effective therapeutic potential through drug repurposing. In this book chapter, we present an integrated multilayer approach that includes cross-omics analyses from publicly available transcriptomics and proteomics data derived from BC tissues and cell lines that were investigated for the development of disease-specific signatures. These signatures are subsequently used as input for a signature-based repurposing approach using the Connectivity Map (CMap) tool. We further explain the steps that may be followed to identify and select existing drugs of increased potential for repurposing in BC patients.

Obstructive Sleep Apnea and Dementia-Common Gene Associations through Network-Based Identification of Common Driver Genes

Background: Obstructive Sleep Apnea (OSA) occurs in 7% of the adult population. The relationship between neurodegenerative diseases such as dementia and sleep disorders have long attracted clinical attention; however, no comprehensive data exists elucidating common gene expression between the two diseases. The objective of this study was to (1) demonstrate the practicability and feasibility of utilizing a systems biology approach called network-based identification of common driver genes (NICD) to identify common genomic features between two associated diseases and (2) utilize this approach to identify genes associated with both OSA and dementia. Methods: This study utilized 2 public databases (PCNet, DisGeNET) and a permutation assay in order to identify common genes between two co-morbid but mutually exclusive diseases. These genes were then linked to their mechanistic pathways through Enrichr, producing a list of genes that were common between the two different diseases. Results: 42 common genes were identified between OSA and dementia which were primarily linked to the G-coupled protein receptor (GPCR) and olfactory pathways. No single nucleotide polymorphisms (SNPs) were identified. Conclusions: This study demonstrates the viability of using publicly available databases and permutation assays along with canonical pathway linkage to identify common gene drivers as potential mechanistic targets for comorbid diseases.

Integrative Survival Analysis of Breast Cancer with Gene Expression and DNA Methylation Data

Motivation

Integrative multi-feature fusion analysis on biomedical data has gained much attention recently. In breast cancer, existing studies have demonstrated that combining genomic mRNA data and DNA methylation data can better stratify cancer patients with distinct prognosis than using single signature. However, those existing methods are simply combining these gene features in series and have ignored the correlations between separate omics dimensions over time.

Results

In the present study, we propose an adaptive multi-task learning method, which combines the Cox loss task with the ordinal loss task, for survival prediction of breast cancer patients using multi-modal learning instead of performing survival analysis on each feature dataset. First, we use local maximum quasi-clique merging (lmQCM) algorithm to reduce the mRNA and methylation feature dimensions and extract cluster eigengenes respectively. Then, we add an auxiliary ordinal loss to the original Cox model to improve the ability to optimize the learning process in training and regularization. The auxiliary loss helps to reduce the vanishing gradient problem for earlier layers and helps to decrease the loss of the primary task. Meanwhile, we use an adaptive weights approach to multi-task learning which weighs multiple loss functions by considering the homoscedastic uncertainty of each task. Finally, we build an ordinal cox hazards model for survival analysis and use long short-term memory (LSTM) method to predict patients’ survival risk. We use the cross-validation method and the concordance index (C-index) for assessing the prediction effect. Stringent cross-verification testing processes for the benchmark dataset and two additional datasets demonstrate that the developed approach is effective, achieving very competitive performance with existing approaches.

Availability and implementation

https://github.com/bhioswego/ML_ordCOX.

Multi-layered network-based pathway activity inference using directed random walks: application to predicting clinical outcomes in urologic cancer

MOTIVATION

To better understand the molecular features of cancers, a comprehensive analysis using multiomics data has been conducted. Additionally, a pathway activity inference method has been developed to facilitate the integrative effects of multiple genes. In this respect, we have recently proposed a novel integrative pathway activity inference approach, iDRW, and demonstrated the effectiveness of the method with respect to dichotomizing two survival groups. However, there were several limitations, such as a lack of generality. In this study, we designed a directed gene-gene graph using pathway information by assigning interactions between genes in multiple layers of networks.

RESULTS

: As a proof-of-concept study, it was evaluated using three genomic profiles of urologic cancer patients. The proposed integrative approach achieved improved outcome prediction performances compared with a single genomic profile alone and other existing pathway activity inference methods. The integrative approach also identified common/cancer-specific candidate driver pathways as predictive prognostic features in urologic cancers. Furthermore, it provides better biological insights into the prioritized pathways and genes in an integrated view using a multi-layered gene-gene network. Our framework is not specifically designed for urologic cancers and can be generally applicable for various datasets.

AVAILABILITY

iDRW is implemented as the R software package. The source codes are available at https://github.com/sykim122/iDRW.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Identification of a Metastasis-Associated Gene Signature of Clear Cell Renal Cell Carcinoma

Clear cell renal cell carcinoma (ccRCC) is one of the most frequent pathological subtypes of kidney cancer, accounting for ~70-75%, and the major cause of mortality is metastatic disease. The difference in gene expression profiles between primary ccRCC tumors and metastatic tumors has not been determined. Thus, we report integrated genomic and transcriptomic analysis for identifying differentially expressed genes (DEGs) between primary and metastatic ccRCC tumors to understand the molecular mechanisms underlying the development of metastases. The microarray datasets GSE105261 and GSE85258 were obtained from the Gene Expression Omnibus (GEO) database, and the R package limma was used for DEG analyses. In summary, the results described herein provide important molecular evidence that metastatic ccRCC tumors are different from primary tumors. Enrichment analysis indicated that the DEGs were mainly enriched in ECM-receptor interaction, platelet activation, protein digestion, absorption, focal adhesion, and the PI3K-Akt signaling pathway. Moreover, we found that DEGs associated with a higher level of tumor immune infiltrates and tumor mutation burden were more susceptible to poor prognosis of ccRCC. Specifically, our study indicates that seven core genes, namely the collagen family (COL1A2, COL1A1, COL6A3, and COL5A1), DCN, FBLN1, and POSTN, were significantly upregulated in metastatic tumors compared with those in primary tumors and, thus, potentially offer insight into novel therapeutic and early diagnostic biomarkers of ccRCC.

An Efficient Algorithm for Identifying Mutated Subnetworks Associated with Survival in Cancer

Protein-protein interaction (PPI) network models interconnections between protein-encoding genes. A group of proteins that perform similar functions are often connected to each other in the PPI network. The corresponding genes form pathways or functional modules. Mutation in protein-encoding genes affect behavior of pathways. This results in initiation, progression, and severity of diseases that propagates through pathways. In this work, we integrate mutation, survival information of patients, and PPI network to identify connected subnetworks associated with survival. We define the computational problem using a fitness function called log-rank statistic to score subnetworks. Log-rank statistic compares the survival between two populations. We propose a novel method, Survival Associated Mutated Subnetwork (SAMS) that adopts genetic algorithm strategy to find the connected subnetwork within the PPI network whose mutation yields highest log-rank statistic. We test on real cancer and synthetic datasets. SAMS generate solutions in negligible time while the state-of-art method in literature takes exponential time. Log-rank statistic of SAMS selected mutated subnetworks are comparable to the method. Our result genesets show significant overlap with well-known cancer driver genes derived from curated datasets and studies in literature, display high text-mining score in terms of number of citations combined with disease-specific keywords in PubMed, and identify pathways having high biological relevance.

Enhanced Permutation Tests via Multiple Pruning

Big multi-omics data in bioinformatics often consists of a huge number of features and relatively small numbers of samples. In addition, features from multi-omics data have their own specific characteristics depending on whether they are from genomics, proteomics, metabolomics, etc. Due to these distinct characteristics, standard statistical analyses using parametric-based assumptions may sometimes fail to provide exact asymptotic results. To resolve this issue, permutation tests can be a way to exactly analyze multi-omics data because they are distribution-free and flexible to use. In permutation tests, p-values are evaluated by estimating the locations of test statistics in an empirical null distribution generated by random shuffling. However, the permutation approach can be infeasible when the number of features increases, because more stringent control of type I error is needed for multiple hypothesis testing, and consequently, much larger numbers of permutations are required to reach significance. To address this problem, we propose a well-organized strategy, “ENhanced Permutation tests via multiple Pruning (ENPP).” ENPP prunes the features in every permutation round if they are determined to be non-significant. In other words, if the feature statistics from the permuted datasets exceed the feature statistics from the original dataset, beyond a predetermined threshold, the feature is determined to be non-significant. If so, ENPP removes the feature and iterates the process without the feature in the next permutation round. Our simulation study showed that the ENPP method could remove about 50% of the features at the first permutation round, and, by the 100th permutation round, 98% of the features had been removed and only 7.4% of the computation time with the original unpruned permutation approach had elapsed. In addition, we applied this approach to a real data set (Korea Association REsource: KARE) of 327,872 SNPs to find association with a non-normally distributed phenotype (fasting plasma glucose), interpreted the results, and discussed the feasibility and advantages of the approach.

New Analysis Framework Incorporating Mixed Mutual Information and Scalable Bayesian Networks for Multimodal High Dimensional Genomic and Epigenomic Cancer Data

We propose a novel two-stage analysis strategy to discover candidate genes associated with the particular cancer outcomes in large multimodal genomic cancers databases, such as The Cancer Genome Atlas (TCGA). During the first stage, we use mixed mutual information to perform variable selection; during the second stage, we use scalable Bayesian network (BN) modeling to identify candidate genes and their interactions. Two crucial features of the proposed approach are (i) the ability to handle mixed data types (continuous and discrete, genomic, epigenomic, etc.) and (ii) a flexible boundary between the variable selection and network modeling stages - the boundary that can be adjusted in accordance with the investigators' BN software scalability and hardware implementation. These two aspects result in high generalizability of the proposed analytical framework. We apply the above strategy to three different TCGA datasets (LGG, Brain Lower Grade Glioma; HNSC, Head and Neck Squamous Cell Carcinoma; STES, Stomach and Esophageal Carcinoma), linking multimodal molecular information (SNPs, mRNA expression, DNA methylation) to two clinical outcome variables (tumor status and patient survival). We identify 11 candidate genes, of which 6 have already been directly implicated in the cancer literature. One novel LGG prognostic factor suggested by our analysis, methylation of TMPRSS11F type II transmembrane serine protease, presents intriguing direction for the follow-up studies.

Prediction of survival risks with adjusted gene expression through risk-gene networks

MOTIVATION

Network-based analysis of biomedical data has been extensively studied over the last decades. As a successful application, gene networks have been used to illustrate interactions among genes and explain the associated phenotypes. However, the gene network approaches have not been actively applied for survival analysis, which is one of the main interests of biomedical research. In addition, a few previous studies using gene networks for survival analysis construct networks mainly from prior knowledge, such as pathways, regulations and gene sets, while the performance considerably depends on the selection of prior knowledge.

RESULTS

In this paper, we propose a data-driven construction method for survival risk-gene networks as well as a survival risk prediction method using the network structure. The proposed method constructs risk-gene networks with survival-associated genes using penalized regression. Then, gene expression indices are hierarchically adjusted through the networks to reduce the variance intrinsic in datasets. By illustrating risk-gene structure, the proposed method is expected to provide an intuition for the relationship between genes and survival risks. The risk-gene network is applied to a low grade glioma dataset, and produces a hypothesis of the relationship between genetic biomarkers of low and high grade glioma. Moreover, with multiple datasets, we demonstrate that the proposed method shows superior prediction performance compared to other conventional methods.

AVAILABILITY AND IMPLEMENTATION

The R package of risk-gene networks is freely available in the web at http://cdal.korea.ac.kr/NetDA/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Topological integration of RPPA proteomic data with multi-omics data for survival prediction in breast cancer via pathway activity inference

BACKGROUND

The analysis of integrated multi-omics data enables the identification of disease-related biomarkers that cannot be identified from a single omics profile. Although protein-level data reflects the cellular status of cancer tissue more directly than gene-level data, past studies have mainly focused on multi-omics integration using gene-level data as opposed to protein-level data. However, the use of protein-level data (such as mass spectrometry) in multi-omics integration has some limitations. For example, the correlation between the characteristics of gene-level data (such as mRNA) and protein-level data is weak, and it is difficult to detect low-abundance signaling proteins that are used to target cancer. The reverse phase protein array (RPPA) is a highly sensitive antibody-based quantification method for signaling proteins. However, the number of protein features in RPPA data is extremely low compared to the number of gene features in gene-level data. In this study, we present a new method for integrating RPPA profiles with RNA-Seq and DNA methylation profiles for survival prediction based on the integrative directed random walk (iDRW) framework proposed in our previous study. In the iDRW framework, each omics profile is merged into a single pathway profile that reflects the topological information of the pathway. In order to address the sparsity of RPPA profiles, we employ the random walk with restart (RWR) approach on the pathway network.

RESULTS

Our model was validated using survival prediction analysis for a breast cancer dataset from The Cancer Genome Atlas. Our proposed model exhibited improved performance compared with other methods that utilize pathway information and also out-performed models that did not include the RPPA data utilized in our study. The risk pathways identified for breast cancer in this study were closely related to well-known breast cancer risk pathways.

CONCLUSIONS

Our results indicated that RPPA data is useful for survival prediction for breast cancer patients under our framework. We also observed that iDRW effectively integrates RNA-Seq, DNA methylation, and RPPA profiles, while variation in the composition of the omics data can affect both prediction performance and risk pathway identification. These results suggest that omics data composition is a critical parameter for iDRW.

ClearF: a supervised feature scoring method to find biomarkers using class-wise embedding and reconstruction

BACKGROUND

Feature selection or scoring methods for the detection of biomarkers are essential in bioinformatics. Various feature selection methods have been developed for the detection of biomarkers, and several studies have employed information-theoretic approaches. However, most of these methods generally require a long processing time. In addition, information-theoretic methods discretize continuous features, which is a drawback that can lead to the loss of information.

RESULTS

In this paper, a novel supervised feature scoring method named ClearF is proposed. The proposed method is suitable for continuous-valued data, which is similar to the principle of feature selection using mutual information, with the added advantage of a reduced computation time. The proposed score calculation is motivated by the association between the reconstruction error and the information-theoretic measurement. Our method is based on class-wise low-dimensional embedding and the resulting reconstruction error. Given multi-class datasets such as a case-control study dataset, low-dimensional embedding is first applied to each class to obtain a compressed representation of the class, and also for the entire dataset. Reconstruction is then performed to calculate the error of each feature and the final score for each feature is defined in terms of the reconstruction errors. The correlation between the information theoretic measurement and the proposed method is demonstrated using a simulation. For performance validation, we compared the classification performance of the proposed method with those of various algorithms on benchmark datasets.

CONCLUSIONS

The proposed method showed higher accuracy and lower execution time than the other established methods. Moreover, an experiment was conducted on the TCGA breast cancer dataset, and it was confirmed that the genes with the highest scores were highly associated with subtypes of breast cancer.

Robust pathway-based multi-omics data integration using directed random walks for survival prediction in multiple cancer studies

Background

Integrating the rich information from multi-omics data has been a popular approach to survival prediction and bio-marker identification for several cancer studies. To facilitate the integrative analysis of multiple genomic profiles, several studies have suggested utilizing pathway information rather than using individual genomic profiles.

Methods

We have recently proposed an integrative directed random walk-based method utilizing pathway information (iDRW) for more robust and effective genomic feature extraction. In this study, we applied iDRW to multiple genomic profiles for two different cancers, and designed a directed gene-gene graph which reflects the interaction between gene expression and copy number data. In the experiments, the performances of the iDRW method and four state-of-the-art pathway-based methods were compared using a survival prediction model which classifies samples into two survival groups.

Results

The results show that the integrative analysis guided by pathway information not only improves prediction performance, but also provides better biological insights into the top pathways and genes prioritized by the model in both the neuroblastoma and the breast cancer datasets. The pathways and genes selected by the iDRW method were shown to be related to the corresponding cancers.

Conclusions

In this study, we demonstrated the effectiveness of a directed random walk-based multi-omics data integration method applied to gene expression and copy number data for both breast cancer and neuroblastoma datasets. We revamped a directed gene-gene graph considering the impact of copy number variation on gene expression and redefined the weight initialization and gene-scoring method. The benchmark result for iDRW with four pathway-based methods demonstrated that the iDRW method improved survival prediction performance and jointly identified cancer-related pathways and genes for two different cancer datasets.

Reviewers

This article was reviewed by Helena Molina-Abril and Marta Hidalgo.

Metabolomics of mammalian brain reveals regional differences

Background

The mammalian brain is organized into regions with specific biological functions and properties. These regions have distinct transcriptomes, but little is known whether they may also differ in their metabolome. The metabolome, a collection of small molecules or metabolites, is at the intersection of the genetic background of a given cell or tissue and the environmental influences that affect it. Thus, the metabolome directly reflects information about the physiologic state of a biological system under a particular condition. The objective of this study was to investigate whether various brain regions have diverse metabolome profiles, similarly to their genetic diversity. The answer to this question would suggest that not only the genome but also the metabolome may contribute to the functional diversity of brain regions.

Methods

We investigated the metabolome of four regions of the mouse brain that have very distinct functions: frontal cortex, hippocampus, cerebellum, and olfactory bulb. We utilized gas- and liquid- chromatography mass spectrometry platforms and identified 215 metabolites.

Results

Principal component analysis, an unsupervised multivariate analysis, clustered each brain region based on its metabolome content, thus providing the unique metabolic profile of each region. A pathway-centric analysis indicated that olfactory bulb and cerebellum had most distinct metabolic profiles, while the cortical parenchyma and hippocampus were more similar in their metabolome content. Among the notable differences were distinct oxidative-anti-oxidative status and region-specific lipid profiles. Finally, a global metabolic connectivity analysis using the weighted correlation network analysis identified five hub metabolites that organized a unique metabolic network architecture within each examined brain region. These data indicate the diversity of global metabolome corresponding to specialized regional brain function and provide a new perspective on the underlying properties of brain regions.

Conclusion

In summary, we observed many differences in the metabolome among the various brain regions investigated. All four brain regions in our study had a unique metabolic signature, but the metabolites came from all categories and were not pathway-centric.

Electronic supplementary material

The online version of this article (10.1186/s12918-018-0644-0) contains supplementary material, which is available to authorized users.

Integrative pathway-based survival prediction utilizing the interaction between gene expression and DNA methylation in breast cancer

Background

Integrative analysis on multi-omics data has gained much attention recently. To investigate the interactive effect of gene expression and DNA methylation on cancer, we propose a directed random walk-based approach on an integrated gene-gene graph that is guided by pathway information.

Methods

Our approach first extracts a single pathway profile matrix out of the gene expression and DNA methylation data by performing the random walk over the integrated graph. We then apply a denoising autoencoder to the pathway profile to further identify important pathway features and genes. The extracted features are validated in the survival prediction task for breast cancer patients.

Results

The results show that the proposed method substantially improves the survival prediction performance compared to that of other pathway-based prediction methods, revealing that the combined effect of gene expression and methylation data is well reflected in the integrated gene-gene graph combined with pathway information. Furthermore, we show that our joint analysis on the methylation features and gene expression profile identifies cancer-specific pathways with genes related to breast cancer.

Conclusions

In this study, we proposed a DRW-based method on an integrated gene-gene graph with expression and methylation profiles in order to utilize the interactions between them. The results showed that the constructed integrated gene-gene graph can successfully reflect the combined effect of methylation features on gene expression profiles. We also found that the selected features by DA can effectively extract topologically important pathways and genes specifically related to breast cancer.

Role of CXCL12-CXCR4 axis in ovarian cancer metastasis and CXCL12-CXCR4 blockade with AMD3100 suppresses tumor cell migration and invasion in vitro

Ovarian cancer (OC) is a lethal gynecologic tumor, which brings its mortality to the head. CXCL12 and its receptor chemokine receptor 4 ( CXCR4) have been found to be highly expressed in OC and contribute to the disease progression by affecting tumor cell proliferation and invasion. Here, in this study, we aim to explore whether the blockade of CXCL12-CXCR4 axis with AMD3100 (a selective CXCR4 antagonist) has effects on the progression of OC. On the basis of the gene expression omnibus database of OC gene expression chips, the OC differentially expressed genes were screened by microarray analysis. OC (nonmetastatic and metastatic) and normal ovarian tissues were collected to determine the expressions of CXCL12 and CXCR4. A series of AMD3100, shRNA against CXCR4, and pCNS-CXCR4 were introduced to treat CAOV3 cells with the highest CXCR4 was assessed. Cell viability, apoptosis, migration, and invasion were all evaluated. The microarray analysis screened out the differential expression of CXCL12-CXCR4 in OC. CXCL12 and CXCR4 expressions were increased in OC tissues, particularly in the metastatic OC tissues. Downregulation of CXCR4 by AMD3100 or shRNA was observed to have a critical role in inhibiting cell proliferation, migration, and invasion of the CAOV3 OC cell line while promoting cell apoptosis. Overexpressed CXCR4 brought significantly promoting effects on the proliferation and invasiveness of OC cells. These results reinforce that the blockade of CXCL12-CXCR4 axis with AMD3100 inhibits the growth of OC cells. The antitumor role of the inhibition of CXCL12-CXCR4 axis offers a preclinical validation of CXCL12-CXCR4 axis as a therapeutic target in OC.

Integrative information theoretic network analysis for genome-wide association study of aspirin exacerbated respiratory disease in Korean population

BACKGROUND

Aspirin Exacerbated Respiratory Disease (AERD) is a chronic medical condition that encompasses asthma, nasal polyposis, and hypersensitivity to aspirin and other non-steroidal anti-inflammatory drugs. Several previous studies have shown that part of the genetic effects of the disease may be induced by the interaction of multiple genetic variants. However, heavy computational cost as well as the complexity of the underlying biological mechanism has prevented a thorough investigation of epistatic interactions and thus most previous studies have typically considered only a small number of genetic variants at a time.

METHODS

In this study, we propose a gene network based analysis framework to identify genetic risk factors from a genome-wide association study dataset. We first derive multiple single nucleotide polymorphisms (SNP)-based epistasis networks that consider marginal and epistatic effects by using different information theoretic measures. Each SNP epistasis network is converted into a gene-gene interaction network, and the resulting gene networks are combined as one for downstream analysis. The integrated network is validated on existing knowledgebase of DisGeNET for known gene-disease associations and GeneMANIA for biological function prediction.

RESULTS

We demonstrated our proposed method on a Korean GWAS dataset, which has genotype information of 440,094 SNPs for 188 cases and 247 controls. The topological properties of the generated networks are examined for scale-freeness, and we further performed various statistical analyses in the Allergy and Asthma Portal (AAP) using the selected genes from our integrated network.

CONCLUSIONS

Our result reveals that there are several gene modules in the network that are of biological significance and have evidence for controlling susceptibility and being related to the treatment of AERD.

Dissecting the relationships of IgG subclasses and complements in membranous lupus nephritis and idiopathic membranous nephropathy

Membranous lupus nephritis (MLN) and idiopathic membranous nephropathy (IMN) are kidney diseases with similar morphology, but distinct etiologies, both producing glomeruli with immune deposits. Immunoglobulins and complements, the main components of the deposits, can be detected by immunofluorescence (IF) microscopy. Previous researches characterized the immune deposits only individually, but not the interactions between them. To study these relationships we analyzed an IF profile of IgG subclasses and complements (IgG1, IgG2, IgG3, IgG4, C3, C1q, and C4) in 53 and 95 cases of biopsy-confirmed MLNs and IMNs, respectively, mainly using information theory and Bayesian networks. We identified significant entropy differences between MLN and IMN for all markers except C3 and IgG1, but mutual information (a measure of mutual dependence) were not significantly different for all the pairs of markers. The entropy differences between MLN and IMN, therefore, were not attributable to the mutual information. These findings suggest that disease type directly and/or indirectly influences the glomerular deposits of most of IgG subclasses and complements, and that the interactions between any pair of the markers were similar between the two diseases. A Markov chain of IgG subclasses was derived from the mutual information about each pair of IgG subclass. Finally we developed an integrated disease model, consistent with the previous findings, describing the glomerular immune deposits of the IgG subclasses and complements based on a Bayesian network using the Markov chain of IgG subclasses as seed. The relationships between the markers were effectively explored by information theory and Bayesian network. Although deposits of IgG subclasses and complements depended on both disease type and the other markers, the interaction between the markers appears conserved, independent from the disease type. The disease model provided an integrated and intuitive representation of the relationships of the IgG subclasses and complements in MLN and IMN.

Use of "big data" in drug discovery and clinical trials

Oncology is undergoing a data-driven metamorphosis. Armed with new and ever more efficient molecular and information technologies, we have entered an era where data is helping us spearhead the fight against cancer. This technology driven data explosion, often referred to as "big data", is not only expediting biomedical discovery, but it is also rapidly transforming the practice of oncology into an information science. This evolution is critical, as results to-date have revealed the immense complexity and genetic heterogeneity of patients and their tumors, a sobering reminder of the challenge facing every patient and their oncologist. This can only be addressed through development of clinico-molecular data analytics that provide a deeper understanding of the mechanisms controlling the biological and clinical response to available therapeutic options. Beyond the exciting implications for improved patient care, such advancements in predictive and evidence-based analytics stand to profoundly affect the processes of cancer drug discovery and associated clinical trials.

Network pharmacology of cancer: From understanding of complex interactomes to the design of multi-target specific therapeutics from nature

Despite massive investments in drug research and development, the significant decline in the number of new drugs approved or translated to clinical use raises the question, whether single targeted drug discovery is the right approach. To combat complex systemic diseases that harbour robust biological networks such as cancer, single target intervention is proved to be ineffective. In such cases, network pharmacology approaches are highly useful, because they differ from conventional drug discovery by addressing the ability of drugs to target numerous proteins or networks involved in a disease. Pleiotropic natural products are one of the promising strategies due to their multi-targeting and due to lower side effects. In this review, we discuss the application of network pharmacology for cancer drug discovery. We provide an overview of the current state of knowledge on network pharmacology, focus on different technical approaches and implications for cancer therapy (e.g. polypharmacology and synthetic lethality), and illustrate the therapeutic potential with selected examples green tea polyphenolics, Eleutherococcus senticosus, Rhodiola rosea, and Schisandra chinensis). Finally, we present future perspectives on their plausible applications for diagnosis and therapy of cancer.