Network-based survival-associated module biomarker and its crosstalk with cell death genes in ovarian cancer

SIG: Graph-Based Cancer Subtype Stratification With Gene Mutation Structural Information

Somatic tumors have a high-dimensional, sparse, and small sample size nature, making cancer subtype stratification based on somatic genomic data a challenge. Current methods for improving cancer clustering performance focus on dimension reduction, integrating multi-omics data, or generating realistic samples, yet ignore the associations between mutated genes within the patient-gene matrix. We refer to these associations as gene mutation structural information, which implicitly includes cancer subtype information and can enhance subtype clustering. We introduce a novel method for cancer subtype clustering called SIG(Structural Information within Graph). As cancer is driven by a combination of genes, we establish associations between mutated genes within the same patient sample, pair by pair, and use a graph to represent them. An association between two mutated genes corresponds to an edge in the graph. We then merge these associations among all mutated genes to obtain a structural information graph, which enriches the gene network and improves its relevance to cancer clustering. We integrate the somatic tumor genome with the enriched gene network and propagate it to cluster patients with mutations in similar network regions. Our method achieves superior clustering performance compared to SOTA methods, as demonstrated by clustering experiments on ovarian and LUAD datasets.

Extraction and Prioritization of a Gene-Cancer-By-Survival Network Involved in Homeostasis of Intracellular Calcium Concentrations Using TCGA PANCAN Data

Regulation of intracellular calcium concentrations, [Ca++]i is important in maintaining the viability of normal as well as cancer cells and can be mediated by tumor microenvironment. Calcium release-activated calcium channel protein (ORAI) calcium channels on the plasma membrane (PM) become physically connected by stromal interaction molecules (STIMs) to the endoplasmic reticulum (ER), on which paralogous receptors of inositol phosphate and of ryanodine are also present along with ATP2A/SERCA (sarco/endoplasmic reticulum calcium ATPases) subunits (also known as PM-ER geneset). Proper expression of this functionally and physically interconnected geneset is essential for the maintenance of [Ca++]i , yet has not been interrogated as a whole for its role in cancer prognosis using multivariable Cox regression. In the present study, we examined whether the expression profile of the PM-ER geneset exhibited prognostic significance across different cancers found in The Cancer Genome Atlas (TCGA) by generating gene-cancer-by-survival networks, in which the nodes represented either genes or cancers and the edges, the logarithmically transformed hazard ratios for overall survival (OS). We then applied network clustering to identify the gene-cancer subnetworks with high connectivity, among which uveal melanoma (UVM) emerged exhibiting the highest degree of genes (k = 10). BAP1, a well-known [Ca++]i regulator and a tumor suppressor, was not found to be significant in predicting OS by PM-ER geneset for UVM, yet it was for several others, including mesothelioma (MESO). Moreover, the best subset of the PM-ER geneset obtained by lasso predicted OS in the TCGA UVM cohort with an area under the receiver operating characteristics (AUC) of 91.4%, comparable to or better than previous prognostic signatures in the literature. Our findings indicate that homeostasis of [Ca++]i is an essential determinant of prognosis in multiple cancers and particularly in UVM. The proposed gene-cancer-by-survival network approach can be extended with other gene sets as well as different survival types.

Survival-Based Biomarker Module Identification Associated with Oral Squamous Cell Carcinoma (OSCC)

Simple Summary

In this study, four OSCC-specific hub genes were identified using high-throughput RNA-Seq data from TCGA cohort. The significant genes within tumor and normal samples were used for weighted PPI network construction based on survival of patients along with their expression profiles. The analysis revealed the most significant module in the training and test datasets. The genes from this module were used for pathway enrichment analysis followed by hub gene selection. These novel biomarkers might have clinical utility for diagnosis and prognosis prediction in OSCC, providing diagnosis at a very early stage. Moreover, a combination of all these biomarkers might distinguish the OSCC patients with low risk and high risk for cancer progression and recurrence, which will provide useful guidance for personalized and precision therapy. However, the results in the present study were obtained by integrative theoretical analysis, and the findings remain to be confirmed by further experimental validations.

Abstract

Head and neck squamous cell carcinoma (HNSC) is one of the most common malignant tumors worldwide with a high rate of morbidity and mortality, with 90% of predilections occurring for oral squamous cell carcinoma (OSCC). Cancers of the mouth account for 40% of head and neck cancers, including squamous cell carcinomas of the tongue, floor of the mouth, buccal mucosa, lips, hard and soft palate, and gingival. OSCC is the most devastating and commonly occurring oral malignancy, with a mortality rate of 500,000 deaths per year. This has imposed a strong necessity to discover driver genes responsible for its progression and malignancy. In the present study we filtered oral squamous cell carcinoma tissue samples from TCGA-HNSC cohort, which we followed by constructing a weighted PPI network based on the survival of patients and the expression profiles of samples collected from them. We found a total of 46 modules, with 18 modules having more than five edges. The KM and ME analyses revealed a single module (with 12 genes) as significant in the training and test datasets. The genes from this significant module were subjected to pathway enrichment analysis for identification of significant pathways and involved genes. Finally, the overlapping genes between gene sets ranked on the basis of weighted PPI module centralities (i.e., degree and eigenvector), significant pathway genes, and DEGs from a microarray OSCC dataset were considered as OSCC-specific hub genes. These hub genes were clinically validated using the IHC images available from the Human Protein Atlas (HPA) database.

Differential Co-Expression Analyses Allow the Identification of Critical Signalling Pathways Altered during Tumour Transformation and Progression

Biological systems respond to perturbations through the rewiring of molecular interactions, organised in gene regulatory networks (GRNs). Among these, the increasingly high availability of transcriptomic data makes gene co-expression networks the most exploited ones. Differential co-expression networks are useful tools to identify changes in response to an external perturbation, such as mutations predisposing to cancer development, and leading to changes in the activity of gene expression regulators or signalling. They can help explain the robustness of cancer cells to perturbations and identify promising candidates for targeted therapy, moreover providing higher specificity with respect to standard co-expression methods. Here, we comprehensively review the literature about the methods developed to assess differential co-expression and their applications to cancer biology. Via the comparison of normal and diseased conditions and of different tumour stages, studies based on these methods led to the definition of pathways involved in gene network reorganisation upon oncogenes' mutations and tumour progression, often converging on immune system signalling. A relevant implementation still lagging behind is the integration of different data types, which would greatly improve network interpretability. Most importantly, performance and predictivity evaluation of the large variety of mathematical models proposed would urgently require experimental validations and systematic comparisons. We believe that future work on differential gene co-expression networks, complemented with additional omics data and experimentally tested, will considerably improve our insights into the biology of tumours.

Kinase-substrate Edge Biomarkers Provide a More Accurate Prognostic Prediction in ER-negative Breast Cancer

The estrogen receptor (ER)-negative breast cancer subtype is aggressive with few treatment options available. To identify specific prognostic factors for ER-negative breast cancer, this study included 705,729 and 1034 breast invasive cancer patients from the Surveillance, Epidemiology, and End Results (SEER) and The Cancer Genome Atlas (TCGA) databases, respectively. To identify key differential kinase-substrate node and edge biomarkers between ER-negative and ER-positive breast cancer patients, we adopted a network-based method using correlation coefficients between molecular pairs in the kinase regulatory network. Integrated analysis of the clinical and molecular data revealed the significant prognostic power of kinase-substrate node and edge features for both subtypes of breast cancer. Two promising kinase-substrate edge features, CSNK1A1-NFATC3 and SRC-OCLN, were identified for more accurate prognostic prediction in ER-negative breast cancer patients.

EMDIP: An Entropy Measure to Discover Important Proteins in PPI networks

Discovering important proteins in Protein-Protein Interaction (PPI) networks has attracted a lot of attention in recent years. Most of the previous work applies different network centrality measures such as Closeness, Betweenness, PageRank and many others to discover the most influential proteins in PPI networks. Although entropy is a well-known graph-based method in computer science, according to our knowledge, it is not used in the biology domain for this purpose. In this paper, first, we annotate the human PPI network with available annotation data. Second, we introduce a new concept called annotation-context that describes each protein according to annotation data of its neighbors. Third, we apply an entropy measure to discover proteins with varied annotation-context. Empirical results indicate that our proposed method succeeded in (1) differentiating essential and non-essential proteins in PPI networks with annotation data; (2) outperforming centrality measures in the task of discovering essential nodes; (3) predicting new annotated proteins based on existing annotation data.

Identification of molecular biomarkers for ovarian cancer using computational approaches

Ovarian cancer is one of the major causes of mortality among women. This is partly because of highly asymptomatic nature, lack of reliable screening techniques and non-availability of effective biomarkers of ovarian cancer. The recent availability of high-throughput data and consequently the development of network medicine approach may play a key role in deciphering the underlying global mechanism involved in a complex disease. This novel approach in medicine will pave the way in translating the new molecular insights into an effective drug therapy applying better diagnostic, prognostic and predictive tests for a complex disease. In this study, we performed reconstruction of gene co-expression networks with a query-based method in healthy and different stages of ovarian cancer to identify new potential biomarkers from the reported biomarker genes. We proposed 17 genes as new potential biomarkers for ovarian cancer that can effectively classify a disease sample from a healthy sample. Most of the predicted genes are found to be differentially expressed between healthy and diseased states. Moreover, the survival analysis showed that these genes have a significantly higher effect on the overall survival rate of the patient than the established biomarkers. The comparative analyses of the co-expression networks across healthy and different stages of ovarian cancer have provided valuable insights into the dynamic nature of ovarian cancer.

Analyzing the regulation of miRNAs on protein-protein interaction network in Hodgkin lymphoma

Background

Hodgkin Lymphoma (HL) is a type of aggressive malignancy in lymphoma that has high incidence in young adults and elderly patients. Identification of reliable diagnostic markers and efficient therapeutic targets are especially important for the diagnosis and treatment of HL. Although many HL-related molecules have been identified, our understanding on the molecular mechanisms underlying the disease is still far from complete due to its complex and heterogeneous characteristics. In such situation, exploring the molecular mechanisms underlying HL via systems biology approaches provides a promising option. In this study, we try to elucidate the molecular mechanisms related to the disease and identify potential pharmaceutical targets from a network-based perspective.

Results

We constructed a series of network models. Based on the analysis of these networks, we attempted to identify the biomarkers and elucidate the molecular mechanisms underlying HL. Initially, we built three different but related protein networks, i.e., background network, HL-basic network and HL-specific network. By analyzing these three networks, we investigated the connection characteristic of the HL-related proteins. Subsequently, we explored the miRNA regulation on HL-specific network and analyzed three kinds of simple regulation patterns, i.e., co-regulation of protein pairs, as well as the direct and indirect regulation of triple proteins. Finally, we constructed a simplified protein network combined with the regulation of miRNAs on proteins to better understand the relation between HL-related proteins and miRNAs.

Conclusions

We find that the HL-related proteins are more likely to connect with each other compared to other proteins. Moreover, the HL-specific network can be further divided into five sub-networks and 49 proteins as the backbone of HL-specific network make up and connect these 5 sub-networks. Thus, they may be closely associated with HL. In addition, we find that the co-regulation of protein pairs is the main regulatory pattern of miRNAs on the protein network in the HL-specific network. According to the regulation of miRNA on protein network, we have identified 5 core miRNAs as the potential biomarkers for diagnostic of HL. Finally, several protein pathways have been identified to closely associated with HL, which provides deep insights into underlying mechanism of HL.

Electronic supplementary material

The online version of this article (10.1186/s12859-019-3041-9) contains supplementary material, which is available to authorized users.

Classifying tumors by supervised network propagation

Motivation

Network propagation has been widely used to aggregate and amplify the effects of tumor mutations using knowledge of molecular interaction networks. However, propagating mutations through interactions irrelevant to cancer leads to erosion of pathway signals and complicates the identification of cancer subtypes.

Results

To address this problem we introduce a propagation algorithm, Network-Based Supervised Stratification (NBS²), which learns the mutated subnetworks underlying tumor subtypes using a supervised approach. Given an annotated molecular network and reference tumor mutation profiles for which subtypes have been predefined, NBS² is trained by adjusting the weights on interaction features such that network propagation best recovers the provided subtypes. After training, weights are fixed such that mutation profiles of new tumors can be accurately classified. We evaluate NBS² on breast and glioblastoma tumors, demonstrating that it outperforms the best network-based approaches in classifying tumors to known subtypes for these diseases. By interpreting the interaction weights, we highlight characteristic molecular pathways driving selected subtypes.

Availability and implementation

The NBS² package is freely available at: https://github.com/wzhang1984/NBSS.

Supplementary information

Supplementary data are available at Bioinformatics online.

A Multi-Cohort and Multi-Omics Meta-Analysis Framework to Identify Network-Based Gene Signatures

Although massive amounts of condition-specific molecular profiles are being accumulated in public repositories every day, meaningful interpretation of these data remains a major challenge. In an effort to identify the biomarkers that describe the key biological phenomena for a given condition, several approaches have been developed over the past few years. However, the majority of these approaches either (i) do not consider the known intermolecular interactions, or (ii) do not integrate molecular data of multiple types (e.g., genomics, transcriptomics, proteomics, epigenomics, etc.), and thus potentially fail to capture the true biological changes responsible for complex diseases (e.g., cancer). In addition, these approaches often ignore the heterogeneity and study bias present in independent molecular cohorts. In this manuscript, we propose a novel multi-cohort and multi-omics meta-analysis framework that overcomes all three limitations mentioned above in order to identify robust molecular subnetworks that capture the key dynamic nature of a given biological condition. Our framework integrates multiple independent gene expression studies, unmatched DNA methylation studies, and protein-protein interactions to identify methylation-driven subnetworks. We demonstrate the proposed framework by constructing subnetworks related to two complex diseases: glioblastoma and low-grade gliomas. We validate the identified subnetworks by showing their ability to predict patients' clinical outcome on multiple independent validation cohorts.

EMT network-based feature selection improves prognosis prediction in lung adenocarcinoma

Various feature selection algorithms have been proposed to identify cancer prognostic biomarkers. In recent years, however, their reproducibility is criticized. The performance of feature selection algorithms is shown to be affected by the datasets, underlying networks and evaluation metrics. One of the causes is the curse of dimensionality, which makes it hard to select the features that generalize well on independent data. Even the integration of biological networks does not mitigate this issue because the networks are large and many of their components are not relevant for the phenotype of interest. With the availability of multi-omics data, integrative approaches are being developed to build more robust predictive models. In this scenario, the higher data dimensions create greater challenges. We proposed a phenotype relevant network-based feature selection (PRNFS) framework and demonstrated its advantages in lung cancer prognosis prediction. We constructed cancer prognosis relevant networks based on epithelial mesenchymal transition (EMT) and integrated them with different types of omics data for feature selection. With less than 2.5% of the total dimensionality, we obtained EMT prognostic signatures that achieved remarkable prediction performance (average AUC values >0.8), very significant sample stratifications, and meaningful biological interpretations. In addition to finding EMT signatures from different omics data levels, we combined these single-omics signatures into multi-omics signatures, which improved sample stratifications significantly. Both single- and multi-omics EMT signatures were tested on independent multi-omics lung cancer datasets and significant sample stratifications were obtained.

Identification of a noncoding RNA‑mediated gene pair‑based regulatory module in Alzheimer's disease

Alzheimer's disease (AD) is the most common type of neurological disorder that results from brain cell death; however, not all brain regions are simultaneously affected to the same extent. Despite single biomarkers for AD having been determined on a genome-wide scale, the differential co-expression in gene pairs between regions and interactions with other types of cellular molecules, particularly non-coding (nc)RNAs, are often overlooked in studies investigating the underlying mechanisms associated with AD. In the present study, based on 1,548 samples obtained from a cohort of 90 patients with AD spanning 19 brain regions, a gene-pair based method was established for the classification of 19 brain regions into seven different groups, including marked disparate groupings of six single regions and a cluster of another 13 regions as revealed by principal component analysis (PCA). To further investigate the different underlying mechanisms associated with each group, five highly interconnected functional modules of the protein-protein interaction network were demonstrated to characterize the seven region groups containing six single groups and 13 clustered regions based on 4,731 gene-pairs. Genes in two of the functional modules exhibited a strong association with pathways associated with the nervous system, including cholinergic synapses, circadian entrainment and dopaminergic synapses. Notably, following integration of these two modules with a ncRNA-mediated network, one module demonstrated a strong association with micro (mi)RNAs, which were revealed to interact with numerous long non-coding (lnc)RNAs associated with AD, such as metastasis associated lung adenocarcinoma transcript 1 and taurine upregulated 1. This suggested that mRNAs and lncRNAs may represent competing endogenous RNAs for binding with miRNAs. Thus, these results indicated that the ncRNA-mediated gene regulatory module detected by the established gene pair-based method may further the understanding of underlying mechanisms associated with AD as well as aid the development of novel therapeutic strategies for the treatment of patients with AD.

Optimizing miRNA-module diagnostic biomarkers of gastric carcinoma via integrated network analysis

Several microRNAs (miRNAs) have been suggested as novel biomarkers for diagnosing gastric cancer (GC) at an early stage, but the single-marker strategy may ignore the co-regulatory relationships and lead to low diagnostic specificity. Thus, multi-target modular diagnostic biomarkers are urgently needed. In this study, a Z_summary and NetSVM-based method was used to identify GC-related hub miRNAs and activated modules from clinical miRNA co-expression networks. The NetSVM-based sub-network consisting of the top 20 hub miRNAs reached a high sensitivity and specificity of 0.94 and 0.82. The Z_summary algorithm identified an activated module (miR-486, miR-451, miR-185, and miR-600) which might serve as diagnostic biomarker of GC. Three members of this module were previously suggested as biomarkers of GC and its 24 target genes were significantly enriched in pathways directly related to cancer. The weighted diagnostic ROC AUC of this module was 0.838, and an optimized module unit (miR-451 and miR-185) obtained a higher value of 0.904, both of which were higher than that of individual miRNAs. These hub miRNAs and module have the potential to become robust biomarkers for early diagnosis of GC with further validations. Moreover, such modular analysis may offer valuable insights into multi-target approaches to cancer diagnosis and treatment.

A computational framework for complex disease stratification from multiple large-scale datasets

BACKGROUND

Multilevel data integration is becoming a major area of research in systems biology. Within this area, multi-'omics datasets on complex diseases are becoming more readily available and there is a need to set standards and good practices for integrated analysis of biological, clinical and environmental data. We present a framework to plan and generate single and multi-'omics signatures of disease states.

METHODS

The framework is divided into four major steps: dataset subsetting, feature filtering, 'omics-based clustering and biomarker identification.

RESULTS

We illustrate the usefulness of this framework by identifying potential patient clusters based on integrated multi-'omics signatures in a publicly available ovarian cystadenocarcinoma dataset. The analysis generated a higher number of stable and clinically relevant clusters than previously reported, and enabled the generation of predictive models of patient outcomes.

CONCLUSIONS

This framework will help health researchers plan and perform multi-'omics big data analyses to generate hypotheses and make sense of their rich, diverse and ever growing datasets, to enable implementation of translational P4 medicine.

Salidroside induces apoptosis in human ovarian cancer SKOV3 and A2780 cells through the p53 signaling pathway

Salidroside is one of the most potent compounds extracted from the plant Rhodiola rosea, and its cardiovascular protective effects have been studied extensively. However, the role of salidroside in human ovarian carcinoma remains unknown. The aim of the current study was to investigate the effects of salidroside on the proliferation and apoptosis of SKOV3 and A2780 cells using MTT assay and acridine orange/ethidium bromide staining. Salidroside activated caspase-3 and upregulated the levels of apoptosis-inducing factor, Bcl-2-associated X and Bcl-2-associated death promoter (Bad) proteins. Furthermore, salidroside downregulated the levels of Bcl-2, p-Bad and X-linked inhibitor of apoptosis proteins. Salidroside activated the caspase-dependent pathway in SKOV3 and A2780 cells, upregulating p53, p21^Cip1/Waf1 and p16^INK4a. These results suggest that the p53/p21^Cip1/Waf1/p16^INK4a pathway may serve a key function in salidroside-mediated effects on SKOV3 and A2780 cells. The current findings indicate that salidroside may be a promising novel drug candidate for ovarian cancer therapy.

Dynamic sporulation gene co-expression networks for Bacillus subtilis 168 and the food-borne isolate Bacillus amyloliquefaciens: a transcriptomic model

Sporulation is a survival strategy, adapted by bacterial cells in response to harsh environmental adversities. The adaptation potential differs between strains and the variations may arise from differences in gene regulation. Gene networks are a valuable way of studying such regulation processes and establishing associations between genes. We reconstructed and compared sporulation gene co-expression networks (GCNs) of the model laboratory strain Bacillus subtilis 168 and the food-borne industrial isolate Bacillus amyloliquefaciens. Transcriptome data obtained from samples of six stages during the sporulation process were used for network inference. Subsequently, a gene set enrichment analysis was performed to compare the reconstructed GCNs of B. subtilis 168 and B. amyloliquefaciens with respect to biological functions, which showed the enriched modules with coherent functional groups associated with sporulation. On basis of the GCNs and time-evolution of differentially expressed genes, we could identify novel candidate genes strongly associated with sporulation in B. subtilis 168 and B. amyloliquefaciens. The GCNs offer a framework for exploring transcription factors, their targets, and co-expressed genes during sporulation. Furthermore, the methodology described here can conveniently be applied to other species or biological processes.

Comprehensive analysis of lncRNA expression profiles reveals a novel lncRNA signature to discriminate nonequivalent outcomes in patients with ovarian cancer

There is growing evidence of dysregulated long non-coding RNAs (lncRNAs) serving as potential biomarkers for cancer prognosis. However, systematic efforts of searching for an expression-based lncRNA signature for prognosis prediction in ovarian cancer (OvCa) have not been made yet. Here, we performed comprehensive analysis for lncRNA expression profiles and clinical data of 544 OvCa patients from The Cancer Genome Atlas (TCGA), and identified an eight-lncRNA signature with ability to classify patients of the training cohort into high-risk group showing poor outcome and low-risk group showing significantly improved outcome, which was further validated in the validation cohort and entire TCGA cohort. Multivariate Cox regression analysis and stratified analysis demonstrated that the prognostic value of this signature was independent of other clinicopathological factors. Associating the outcome prediction with BRCA1 and/or BRCA2 mutation revealed a superior prognosis performance both in BRCA1/2-mutated and BRCA1/2 wild-type tumors. Finally, a significantly correlation was found between the lncRNA signature and the complete response rate of chemotherapy, suggesting that this eight-lncRNA signature may be a measure to predict chemotherapy response and identify platinum-resistant patients who might benefit from other more efficacious therapies. With further prospective validation, this eight-lncRNA signature may have important implications for outcome prediction and therapy decisions.

Identifying pathway modules of tuberculosis in children by analyzing multiple different networks

Tuberculosis (TB), which is caused by the mycobacterium TB, is the major cause of human death worldwide. The aim of this study was to identify the biomarkers involved in child TB. Gene expression data were obtained from the Array Express Archive of Functional Genomics Data. Gene expression data and protein-protein interaction (PPI) data were downloaded to construct differential gene co-expression networks (DCNs). The Benjamini-Hochberg algorithm was used to correct the P-value. In total, 3,820 edges (PPIs) and 1,359 nodes (genes) were obtained from the human-related PPIs data and gene expression data at the criteria of absolute value of Pearson's correlation coefficient >0.8. The DCNs were formed by these edges and nodes. Thirteen seed genes were obtained by ranging z-scores. Eight significant multiple different modules were identified from DCNs using the statistical significant test. In conclusion, the seed genes and significant modules constitute potential biomarkers that reveal the underlying mechanisms in child TB. The new identified biomarkers may contribute to an understanding of TB and provide a new therapeutic method for the treatment of TB.

RNA-based ovarian cancer research from 'a gene to systems biomedicine' perspective

Ovarian cancer remains the leading cause of death from a gynecologic malignancy, and treatment of this disease is harder than any other type of female reproductive cancer. Improvements in the diagnosis and development of novel and effective treatment strategies for complex pathophysiologies, such as ovarian cancer, require a better understanding of disease emergence and mechanisms of progression through systems medicine approaches. RNA-level analyses generate new information that can help in understanding the mechanisms behind disease pathogenesis, to identify new biomarkers and therapeutic targets and in new drug discovery. Whole RNA sequencing and coding and non-coding RNA expression array datasets have shed light on the mechanisms underlying disease progression and have identified mRNAs, miRNAs, and lncRNAs involved in ovarian cancer progression. In addition, the results from these analyses indicate that various signalling pathways and biological processes are associated with ovarian cancer. Here, we present a comprehensive literature review on RNA-based ovarian cancer research and highlight the benefits of integrative approaches within the systems biomedicine concept for future ovarian cancer research. We invite the ovarian cancer and systems biomedicine research fields to join forces to achieve the interdisciplinary caliber and rigor required to find real-life solutions to common, devastating, and complex diseases such as ovarian cancer.

ABBREVIATIONS

CAF: cancer-associated fibroblasts; COG: Cluster of Orthologous Groups; DEA: disease enrichment analysis; EOC: epithelial ovarian carcinoma; ESCC: oesophageal squamous cell carcinoma; GSI: gamma secretase inhibitor; GO: Gene Ontology; GSEA: gene set enrichment analyzes; HAS: Hungarian Academy of Sciences; lncRNAs: long non-coding RNAs; MAPK/ERK: mitogen-activated protein kinase/extracellular signal-regulated kinases; NGS: next-generation sequencing; ncRNAs: non-coding RNAs; OvC: ovarian cancer; PI3K/Akt/mTOR: phosphatidylinositol-3-kinase/protein kinase B/mammalian target of rapamycin; RT-PCR: real-time polymerase chain reaction; SNP: single nucleotide polymorphism; TF: transcription factor; TGF-β: transforming growth factor-β.

Identification of oral cancer related candidate genes by integrating protein-protein interactions, gene ontology, pathway analysis and immunohistochemistry

In the recent years, bioinformatics methods have been reported with a high degree of success for candidate gene identification. In this milieu, we have used an integrated bioinformatics approach assimilating information from gene ontologies (GO), protein–protein interaction (PPI) and network analysis to predict candidate genes related to oral squamous cell carcinoma (OSCC). A total of 40973 PPIs were considered for 4704 cancer-related genes to construct human cancer gene network (HCGN). The importance of each node was measured in HCGN by ten different centrality measures. We have shown that the top ranking genes are related to a significantly higher number of diseases as compared to other genes in HCGN. A total of 39 candidate oral cancer target genes were predicted by combining top ranked genes and the genes corresponding to significantly enriched oral cancer related GO terms. Initial verification using literature and available experimental data indicated that 29 genes were related with OSCC. A detailed pathway analysis led us to propose a role for the selected candidate genes in the invasion and metastasis in OSCC. We further validated our predictions using immunohistochemistry (IHC) and found that the gene FLNA was upregulated while the genes ARRB1 and HTT were downregulated in the OSCC tissue samples.

A powerful weighted statistic for detecting group differences of directed biological networks

Complex disease is largely determined by a number of biomolecules interwoven into networks, rather than a single biomolecule. Different physiological conditions such as cases and controls may manifest as different networks. Statistical comparison between biological networks can provide not only new insight into the disease mechanism but statistical guidance for drug development. However, the methods developed in previous studies are inadequate to capture the changes in both the nodes and edges, and often ignore the network structure. In this study, we present a powerful weighted statistical test for group differences of directed biological networks, which is independent of the network attributes and can capture the changes in both the nodes and edges, as well as simultaneously accounting for the network structure through putting more weights on the difference of nodes locating on relatively more important position. Simulation studies illustrate that this method had better performance than previous ones under various sample sizes and network structures. One application to GWAS of leprosy successfully identifies the specific gene interaction network contributing to leprosy. Another real data analysis significantly identifies a new biological network, which is related to acute myeloid leukemia. One potential network responsible for lung cancer has also been significantly detected. The source R code is available on our website.

An integrative and comparative study of pan-cancer transcriptomes reveals distinct cancer common and specific signatures

To investigate the commonalities and specificities across tumor lineages, we perform a systematic pan-cancer transcriptomic study across 6744 specimens. We find six pan-cancer subnetwork signatures which relate to cell cycle, immune response, Sp1 regulation, collagen, muscle system and angiogenesis. Moreover, four pan-cancer subnetwork signatures demonstrate strong prognostic potential. We also characterize 16 cancer type-specific subnetwork signatures which show diverse implications to somatic mutations, somatic copy number aberrations, DNA methylation alterations and clinical outcomes. Furthermore, some of them are strongly correlated with histological or molecular subtypes, indicating their implications with tumor heterogeneity. In summary, we systematically explore the pan-cancer common and cancer type-specific gene subnetwork signatures across multiple cancers, and reveal distinct commonalities and specificities among cancers at transcriptomic level.

Differential Regulatory Analysis Based on Coexpression Network in Cancer Research

With rapid development of high-throughput techniques and accumulation of big transcriptomic data, plenty of computational methods and algorithms such as differential analysis and network analysis have been proposed to explore genome-wide gene expression characteristics. These efforts are aiming to transform underlying genomic information into valuable knowledges in biological and medical research fields. Recently, tremendous integrative research methods are dedicated to interpret the development and progress of neoplastic diseases, whereas differential regulatory analysis (DRA) based on gene coexpression network (GCN) increasingly plays a robust complement to regular differential expression analysis in revealing regulatory functions of cancer related genes such as evading growth suppressors and resisting cell death. Differential regulatory analysis based on GCN is prospective and shows its essential role in discovering the system properties of carcinogenesis features. Here we briefly review the paradigm of differential regulatory analysis based on GCN. We also focus on the applications of differential regulatory analysis based on GCN in cancer research and point out that DRA is necessary and extraordinary to reveal underlying molecular mechanism in large-scale carcinogenesis studies.

A VEGF-dependent gene signature enriched in mesenchymal ovarian cancer predicts patient prognosis

We have previously reported surrogate biomarkers of VEGF pathway activities with the potential to provide predictive information for anti-VEGF therapies. The aim of this study was to systematically evaluate a new VEGF-dependent gene signature (VDGs) in relation to molecular subtypes of ovarian cancer and patient prognosis. Using microarray profiling and cross-species analysis, we identified 140-gene mouse VDGs and corresponding 139-gene human VDGs, which displayed enrichment of vasculature and basement membrane genes. In patients who received bevacizumab therapy and showed partial response, the expressions of VDGs (summarized to yield VDGs scores) were markedly decreased in post-treatment biopsies compared with pre-treatment baselines. In contrast, VDGs scores were not significantly altered following bevacizumab treatment in patients with stable or progressive disease. Analysis of VDGs in ovarian cancer showed that VDGs as a prognostic signature was able to predict patient outcome. Correlation estimation of VDGs scores and molecular features revealed that VDGs was overrepresented in mesenchymal subtype and BRCA mutation carriers. These findings highlighted the prognostic role of VEGF-mediated angiogenesis in ovarian cancer, and proposed a VEGF-dependent gene signature as a molecular basis for developing novel diagnostic strategies to aid patient selection for VEGF-targeted agents.

A Five-Gene Expression Signature Predicts Clinical Outcome of Ovarian Serous Cystadenocarcinoma

Ovarian serous cystadenocarcinoma is a common malignant tumor of female genital organs. Treatment is generally less effective as patients are usually diagnosed in the late stage. Therefore, a well-designed prognostic marker provides valuable data for optimizing therapy. In this study, we analyzed 303 samples of ovarian serous cystadenocarcinoma and the corresponding RNA-seq data. We observed the correlation between gene expression and patients' survival and eventually established a risk assessment model of five factors using Cox proportional hazards regression analysis. We found that the survival time in high-risk patients was significantly shorter than in low-risk patients in both training and testing sets after Kaplan-Meier analysis. The AUROC value was 0.67 when predicting the survival time in testing set, which indicates a relatively high specificity and sensitivity. The results suggest diagnostic and therapeutic applications of our five-gene model for ovarian serous cystadenocarcinoma.

A network-based method for identifying prognostic gene modules in lung squamous carcinoma

Similarities in gene expression between both developing embryonic and precancerous tissues and cancer tissues may help identify much-needed biomarkers and therapeutic targets in lung squamous carcinoma. In this study, human lung samples representing ten successive time points, from embryonic development to carcinogenesis, were used to construct global gene expression profiles. Differentially expressed genes with similar expression in precancerous and cancer samples were identified. Using a network-based greedy searching algorithm to analyze the training cohort (n = 69) and three independent testing cohorts, we successfully identified a significant 22-gene module in which expression levels were correlated with overall survival in lung squamous carcinoma patients.

Predicting cancerlectins by the optimal g-gap dipeptides

The cancerlectin plays a key role in the process of tumor cell differentiation. Thus, to fully understand the function of cancerlectin is significant because it sheds light on the future direction for the cancer therapy. However, the traditional wet-experimental methods were money- and time-consuming. It is highly desirable to develop an effective and efficient computational tool to identify cancerlectins. In this study, we developed a sequence-based method to discriminate between cancerlectins and non-cancerlectins. The analysis of variance (ANOVA) was used to choose the optimal feature set derived from the g-gap dipeptide composition. The jackknife cross-validated results showed that the proposed method achieved the accuracy of 75.19%, which is superior to other published methods. For the convenience of other researchers, an online web-server CaLecPred was established and can be freely accessed from the website http://lin.uestc.edu.cn/server/CalecPred. We believe that the CaLecPred is a powerful tool to study cancerlectins and to guide the related experimental validations.

Revealing potential molecular targets bridging colitis and colorectal cancer based on multidimensional integration strategy

Chronic inflammation may play a vital role in the pathogenesis of inflammation-associated tumors. However, the underlying mechanisms bridging ulcerative colitis (UC) and colorectal cancer (CRC) remain unclear. Here, we integrated multidimensional interaction resources, including gene expression profiling, protein-protein interactions (PPIs), transcriptional and post-transcriptional regulation data, and virus-host interactions, to tentatively explore potential molecular targets that functionally link UC and CRC at a systematic level. In this work, by deciphering the overlapping genes, crosstalking genes and pivotal regulators of both UC- and CRC-associated functional module pairs, we revealed a variety of genes (including FOS and DUSP1, etc.), transcription factors (including SMAD3 and ETS1, etc.) and miRNAs (including miR-155 and miR-196b, etc.) that may have the potential to complete the connections between UC and CRC. Interestingly, further analyses of the virus-host interaction network demonstrated that several virus proteins (including EBNA-LP of EBV and protein E7 of HPV) frequently inter-connected to UC- and CRC-associated module pairs with their validated targets significantly enriched in both modules of the host. Together, our results suggested that multidimensional integration strategy provides a novel approach to discover potential molecular targets that bridge the connections between UC and CRC, which could also be extensively applied to studies on other inflammation-related cancers.

Stroma-associated master regulators of molecular subtypes predict patient prognosis in ovarian cancer

High-grade serous ovarian carcinoma (HGS-OvCa) has the lowest survival rate among all gynecologic cancers and is hallmarked by a high degree of heterogeneity. The Cancer Genome Atlas network has described a gene expression-based molecular classification of HGS-OvCa into Differentiated, Mesenchymal, Immunoreactive and Proliferative subtypes. However, the biological underpinnings and regulatory mechanisms underlying the distinct molecular subtypes are largely unknown. Here we showed that tumor-infiltrating stromal cells significantly contributed to the assignments of Mesenchymal and Immunoreactive clusters. Using reverse engineering and an unbiased interrogation of subtype regulatory networks, we identified the transcriptional modules containing master regulators that drive gene expression of Mesenchymal and Immunoreactive HGS-OvCa. Mesenchymal master regulators were associated with poor prognosis, while Immunoreactive master regulators positively correlated with overall survival. Meta-analysis of 749 HGS-OvCa expression profiles confirmed that master regulators as a prognostic signature were able to predict patient outcome. Our data unraveled master regulatory programs of HGS-OvCa subtypes with prognostic and potentially therapeutic relevance, and suggested that the unique transcriptional and clinical characteristics of ovarian Mesenchymal and Immunoreactive subtypes could be, at least partially, ascribed to tumor microenvironment.