The mRNA–miRNA–lncRNA Regulatory Network and Factors Associated with Prognosis Prediction of Hepatocellular Carcinoma

The aim of this study was to identify novel prognostic mRNA and microRNA (miRNA) biomarkers for hepatocellular carcinoma (HCC) using methods in systems biology. Differentially expressed mRNAs, miRNAs, and long non-coding RNAs (lncRNAs) were compared between HCC tumor tissues and normal liver tissues in The Cancer Genome Atlas (TCGA) database. Subsequently, a prognosis-associated mRNA co-expression network, an mRNA–miRNA regulatory network, and an mRNA–miRNA–lncRNA regulatory network were constructed to identify prognostic biomarkers for HCC through Cox survival analysis. Seven prognosis-associated mRNA co-expression modules were obtained by analyzing these differentially expressed mRNAs. An expression module including 120 mRNAs was significantly correlated with HCC patient survival. Combined with patient survival data, several mRNAs and miRNAs, including CHST4, SLC22A8, STC2, hsa-miR-326, and hsa-miR-21 were identified from the network to predict HCC patient prognosis. Clinical significance was investigated using tissue microarray analysis of samples from 258 patients with HCC. Functional annotation of hsa-miR-326 and hsa-miR-21-5p indicated specific associations with several cancer-related pathways. The present study provides a bioinformatics method for biomarker screening, leading to the identification of an integrated mRNA–miRNA–lncRNA regulatory network and their co-expression patterns in relation to predicting HCC patient survival.     

Weighted gene correlation network analysis identifies RSAD2, HERC5, and CCL8 as prognostic candidates for breast cancer

As the most commonly diagnosed malignant tumor in female population, the prognosis of breast cancer is affected by complex gene interaction networks. In this research weighted gene co-expression network analysis (WGCNA) would be utilized to build a gene co-expression network to identify potential biomarkers for prediction the prognosis of patients with breast cancer. We downloaded GSE25065 from Gene Expression Omnibus database as the test set. GSE25055 and GSE42568 were utilized to validate findings in the research. Seven modules were established in the GSE25065 by utilizing average link hierarchical clustering. Three hub genes, RSAD2, HERC5, and CCL8 were screened out from the significant module (R ² = 0.44), which were considerably interrelated to worse prognosis. Within test dataset GSE25065, RSAD2, and CCL8 were correlated with tumor stage, grade, and lymph node metastases, whereas HERC5 was correlated with lymph node metastases and tumor grade. In the validation dataset GSE25055 and RSAD2 expression was correlated with tumor grade, stage, and size, whereas HERC5 was related to tumor stage and tumor grade, and CCL8 was associated with tumor size and tumor grade. Multivariable survival analysis demonstrated that RSAD2, HERC5, and CCL8 were independent risk factors. In conclusion, the WGCNA analysis conducted in this study screened out novel prognostic biomarkers of breast cancer. Meanwhile, further in vivo and in vitro studies are required to make the clear molecular mechanisms.     

PNME – A gene-gene parallel network module extraction method

In the domain of gene-gene network analysis, construction of co-expression networks and extraction of network modules have opened up enormous possibilities for exploring the role of genes in biological processes. Through such analysis, one can extract interesting behaviour of genes and would help in the discovery of genes participating in a common biological process. However, such network analysis methods in sequential processing mode often have been found time-consuming even for a moderately sized dataset.It is observed that most existing network construction techniques are capable of handling only positive correlations in gene-expression data whereas biologically-significant genes exhibit both positive and negative correlations. To address these problems, we propose a faster method for construction and analysis of gene-gene network and extraction of modules using a similarity measure which can identify both negatively and positively correlated co-expressed patterns. Our method utilizes General-purpose computing on graphics processing units (GPGPU) to provide fast, efficient and parallel extraction of biologically relevant network modules to support biomarker identification for breast cancer. The modules extracted are validated using p-value and q-value for both metastasis and non-metastasis stages of breast cancer. PNME has been found capable of identifying interesting biomarkers for this critical disease. We identified six genes with the interesting behaviours which have been found to cause breast cancer in homo-sapiens.     

Identification of an extracellular vesicle-related gene signature in the prediction of pancreatic cancer clinical prognosis

Although extracellular vesicles (EVs) in body fluid have been considered to be ideal biomarkers for cancer diagnosis and prognosis, it is still difficult to distinguish EVs derived from tumor tissue and normal tissue. Therefore, the prognostic value of tumor-specific EVs was evaluated through related molecules in pancreatic tumor tissue. NA sequencing data of pancreatic adenocarcinoma (PAAD) were acquired from The Cancer Genome Atlas (TCGA) and International Cancer Genome Consortium (ICGC). EV-related genes in pancreatic cancer were obtained from exoRBase. Protein–protein interaction (PPI) network analysis was used to identify modules related to clinical stage. CIBERSORT was used to assess the abundance of immune and non-immune cells in the tumor microenvironment. A total of 12 PPI modules were identified, and the 3-PPI-MOD was identified based on the randomForest package. The genes of this model are involved in DNA damage and repair and cell membrane-related pathways. The independent external verification cohorts showed that the 3-PPI-MOD can significantly classify patient prognosis. Moreover, compared with the model constructed by pure gene expression, the 3-PPI-MOD showed better prognostic value. The expression of genes in the 3-PPI-MOD had a significant positive correlation with immune cells. Genes related to the hypoxia pathway were significantly enriched in the high-risk tumors predicted by the 3-PPI-MOD. External databases were used to verify the gene expression in the 3-PPI-MOD. The 3-PPI-MOD had satisfactory predictive performance and could be used as a prognostic predictive biomarker for pancreatic cancer.     

Integrated analysis of co‐expression and ceRNA network identifies five lncRNAs as prognostic markers for breast cancer

Long non‐coding RNAs (lncRNAs), which competitively bind miRNAs to regulate target mRNA expression in the competing endogenous RNAs (ceRNAs) network, have attracted increasing attention in breast cancer research. We aim to find more effective therapeutic targets and prognostic markers for breast cancer. LncRNA, mRNA and miRNA expression profiles of breast cancer were downloaded from TCGA database. We screened the top 5000 lncRNAs, top 5000 mRNAs and all miRNAs to perform weighted gene co‐expression network analysis. The correlation between modules and clinical information of breast cancer was identified by Pearson's correlation coefficient. Based on the most relevant modules, we constructed a ceRNA network of breast cancer. Additionally, the standard Kaplan‐Meier univariate curve analysis was adopted to identify the prognosis of lncRNAs. Ultimately, a total of 23 and 5 modules were generated in the lncRNAs/mRNAs and miRNAs co‐expression network, respectively. According to the Green module of lncRNAs/mRNAs and Blue module of miRNAs, our constructed ceRNA network consisted of 52 lncRNAs, 17miRNAs and 79 mRNAs. Through survival analysis, 5 lncRNAs (AL117190.1, COL4A2‐AS1, LINC00184, MEG3 and MIR22HG) were identified as crucial prognostic factors for patients with breast cancer. Taken together, we have identified five novel lncRNAs related to prognosis of breast cancer. Our study has contributed to the deeper understanding of the molecular mechanism of breast cancer and provided novel insights into the use of breast cancer drugs and prognosis.     

Detecting microRNAs of high influence on protein functional interaction networks: a prostate cancer case study

ABSTRACT: BACKGROUND: The use of biological molecular network information for diagnostic and prognostic purposes and elucidation of molecular disease mechanism is a key objective in systems biomedicine. The network of regulatory miRNA-target and functional protein interactions is a rich source of information to elucidate the function and the prognostic value of miRNAs in cancer. The objective of this study is to identify miRNAs that have high influence on target protein complexes in prostate cancer as a case study. This could provide biomarkers or therapeutic targets relevant for prostate cancer treatment. RESULTS: Our findings demonstrate that a miRNA's functional role can be explained by its target protein connectivity within a physical and functional interaction network. To detect miRNAs with high influence on target protein modules, we integrated miRNA and mRNA expression profiles with a sequence based miRNA-target network and human functional and physical protein interactions (FPI). miRNAs with high influence on target protein complexes play a role in prostate cancer progression and are promising diagnostic or prognostic biomarkers. We uncovered several miRNA-regulated protein modules which were enriched in focal adhesion and prostate cancer genes. Several miRNAs such as miR-96, miR-182, and miR-143 demonstrated high influence on their target protein complexes and could explain most of the gene expression changes in our analyzed prostate cancer data set. CONCLUSIONS: We describe a novel method to identify active miRNA-target modules relevant to prostate cancer progression and outcome. miRNAs with high influence on protein networks are valuable biomarkers that can be used in clinical investigations for prostate cancer treatment.     

Coexpression network analysis linked H2AFJ to chemoradiation resistance in colorectal cancer

Neoadjuvant chemoradiotherapy (CRT) resistance is a complex phenomenon and it remains a major problem for patients with a priori resistant tumor. Therefore, there is a strong need to investigate molecular biomarkers which may guide for treatment decision-making. In our study, weighted gene coexpression network analysis was applied to identify CRT-resistance hub modules in 12 colorectal cancer (CRC) cell lines with different CRT sensitivities from GSE20298 data set. The green module and purple module had the highest correlations with CRT resistance. Gene ontology enrichment analysis indicated that the function of these two modules focused on interferon-mediated signaling pathway, immune response, chromatin modulation, Rho GTPases activities, and regulation of apoptotic process. Then, 15 hub genes in both the coexpression and protein-protein interaction networks were selected. Among these hub genes, higher H2A histone family member J (H2AFJ) expression was independently validated in patient cohorts from two testing data sets of GSE46862 and GSE68204 to be related to CRT resistance. The receiver operating characteristic curve showed that H2AFJ could efficiently distinguish CRT-resistance cases from CRT-sensitive cases in another two testing data sets. Furthermore, meta-analysis of 12 Gene Expression Omnibus–sourced data sets showed that H2AFJ messenger RNA levels were significantly higher in CRC tissues than in normal colon tissues. High H2AFJ expression was correlated with a significant worse event- and relapse-free survival by analyzing the data from the R2: Genomics Analysis and Visualization Platform. Gene set enrichment analysis determined that the mechanism of H2AFJ-mediated CRT resistance might involve the ERK5 (MAPK7), human immunodeficiency virus Nef (HIV Nef), and inflammatory pathways. This study is the first, to the best of our knowledge, to implicate and verify H2AFJ as an effective new marker for CRT response prediction.     

Construction of miRNA‐lncRNA‐mRNA co‐expression network affecting EMT‐mediated cisplatin resistance in ovarian cancer

Platinum resistance is one of the major concerns in ovarian cancer treatment. Recent evidence shows the critical role of epithelial–mesenchymal transition (EMT) in this resistance. Epithelial‐like ovarian cancer cells show decreased sensitivity to cisplatin after cisplatin treatment. Our study prospected the association between epithelial phenotype and response to cisplatin in ovarian cancer. Microarray dataset {"type":"entrez-geo","attrs":{"text":"GSE47856","term_id":"47856"}}GSE47856 was acquired from the GEO database. After identifying differentially expressed genes (DEGs) between epithelial‐like and mesenchymal‐like cells, the module identification analysis was performed using weighted gene co‐expression network analysis (WGCNA). The gene ontology (GO) and pathway analyses of the most considerable modules were performed. The protein–protein interaction network was also constructed. The hub genes were specified using Cytoscape plugins MCODE and cytoHubba, followed by the survival analysis and data validation. Finally, the co‐expression of miRNA‐lncRNA‐TF with the hub genes was reconstructed. The co‐expression network analysis suggests 20 modules relating to the Epithelial phenotype. The antiquewhite4, brown and darkmagenta modules are the most significant non‐preserved modules in the Epithelial phenotype and contain the most differentially expressed genes. GO, and KEGG pathway enrichment analyses on these modules divulge that these genes were primarily enriched in the focal adhesion, DNA replication pathways and stress response processes. ROC curve and overall survival rate analysis show that the co‐expression pattern of the brown module''s hub genes could be a potential prognostic biomarker for ovarian cancer cisplatin resistance.     

Identification of prognostic signature in cancer based on DNA methylation interaction network

Background

The identification of prognostic biomarkers for cancer patients is essential for cancer research. These days, DNA methylation has been proved to be associated with cancer prognosis. However, there are few methods which identify the prognostic markers based on DNA methylation data systematically, especially considering the interaction among DNA methylation sites.

Methods

In this paper, we first evaluated the stabilities of microRNA, mRNA, and DNA methylation data in prognosis of cancer. After that, a rank-based method was applied to construct a DNA methylation interaction network. In this network, nodes with the largest degrees (10% of all the nodes) were selected as hubs. Cox regression was applied to select the hubs as prognostic signature. In this prognostic signature, DNA methylation levels of each DNA methylation site are correlated with the outcomes of cancer patients. After obtaining these prognostic genes, we performed the survival analysis in the training group and the test group to verify the reliability of these genes.

Results

We applied our method in three cancers (ovarian cancer, breast cancer and Glioblastoma Multiforme).In all the three cancers, there are more common ones of prognostic genes selected from different samples in DNA methylation data, compared with gene expression data and miRNA expression data, which indicates the DNA methylation data may be more stable in cancer prognosis. Power-law distribution fitting suggests that the DNA methylation interaction networks are scale-free. And the hubs selected from the three networks are all enriched by cancer related pathways. The gene signatures were obtained for the three cancers respectively, and survival analysis shows they can distinguish the outcomes of tumor patients in both the training data sets and test data sets, which outperformed the control signatures.

Conclusions

A computational method was proposed to construct DNA methylation interaction network and this network could be used to select prognostic signatures in cancer.

     

Genomic Heterogeneity of Osteosarcoma - Shift from Single Candidates to Functional Modules

Osteosarcoma (OS), a bone tumor, exhibit a complex karyotype. On the genomic level a highly variable degree of alterations in nearly all chromosomal regions and between individual tumors is observable. This hampers the identification of common drivers in OS biology. To identify the common molecular mechanisms involved in the maintenance of OS, we follow the hypothesis that all the copy number-associated differences between the patients are intercepted on the level of the functional modules. The implementation is based on a network approach utilizing copy number associated genes in OS, paired expression data and protein interaction data. The resulting functional modules of tightly connected genes were interpreted regarding their biological functions in OS and their potential prognostic significance. We identified an osteosarcoma network assembling well-known and lesser-known candidates. The derived network shows a significant connectivity and modularity suggesting that the genes affected by the heterogeneous genetic alterations share the same biological context. The network modules participate in several critical aspects of cancer biology like DNA damage response, cell growth, and cell motility which is in line with the hypothesis of specifically deregulated but functional modules in cancer. Further, we could deduce genes with possible prognostic significance in OS for further investigation (e.g. EZR, CDKN2A, MAP3K5). Several of those module genes were located on chromosome 6q. The given systems biological approach provides evidence that heterogeneity on the genomic and expression level is ordered by the biological system on the level of the functional modules. Different genomic aberrations are pointing to the same cellular network vicinity to form vital, but already neoplastically altered, functional modules maintaining OS. This observation, exemplarily now shown for OS, has been under discussion already for a longer time, but often in a hypothetical manner, and can here be exemplified for OS.     

Identification of key gene modules and pathways of human glioma through coexpression network

Glioma causes great harm to people worldwide. Systemic coexpression analysis of this disease could be beneficial for the identification and development of new prognostic and predictive markers in the clinical management of glioma. In this study, we extracted data sets from the Gene Expression Omnibus data set by using “glioma” as the keyword. Then, a coexpression module was constructed with the help of Weighted Gene Coexpression Network Analysis software. Besides, Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses were performed on the genes in these modules. As a result, the critical modules and target genes were identified. Eight coexpression modules were constructed using the 4,000 genes with a high expression value of the total 141 glioma samples. The result of the analysis of the interaction among these modules showed that there was a high scale independence degree among them. The GO and KEGG enrichment analyses showed that there was a significant difference in the enriched terms and degree among these eight modules, and module 5 was identified as the most important module. Besides, the pathways it was enriched in, hsa04510: Focal adhesion and hsa04610: Complement and coagulation cascades, were determined as the most important pathways. In summary, module 5 and the pathways it was enriched in, hsa04510: Focal adhesion and has 04610: Complement and coagulation cascades, have the potential to serve as biomarkers for patients with glioma.     

Novel key genes in triple-negative breast cancer identified by weighted gene co-expression network analysis

Triple-negative breast cancer (TNBC) is a special subtype of breast cancer (BC) with poor prognosis. Although some molecular mechanisms of TNBC have been elucidated, the efficacy of current treatments is limited. Therefore, it is urgently demanded to screen for novel biomarkers and drug targets for TNBC. In this study, we obtained four independent data sets (GSE76250, GSE31448, GSE43358, and METABRIC) from the Gene Expression Omnibus (GEO) database and the cBioPortal website. In the GSE76250 data set, 890 differentially expressed genes were identified and weighted gene co-expression network analysis was performed based on them. Then, two preserved modules associated with the KI67 score were detected. Gene ontology and pathway enrichment analyses showed genes in the modules participated in some cancer-related biological processes or pathways. Non-SMC condensin I complex subunit G (NCAPG) and ATP-binding cassette subfamily A member 9 (ABCA9) were identified as hub genes of the modules, and the significance of hub genes was validated in the GSE43358 data set. Finally, their prognostic value was assessed by survival analysis. These findings suggested that NCAPG and ABCA9 may be the key genes of TNBC. Moreover, ABCA9 was first reported in TNBC. They deserved further studies.     

High accordance in prognosis prediction of colorectal cancer across independent datasets by multi-gene module expression profiles

A considerable portion of patients with colorectal cancer have a high risk of disease recurrence after surgery. These patients can be identified by analyzing the expression profiles of signature genes in tumors. But there is no consensus on which genes should be used and the performance of specific set of signature genes varies greatly with different datasets, impeding their implementation in the routine clinical application. Instead of using individual genes, here we identified functional multi-gene modules with significant expression changes between recurrent and recurrence-free tumors, used them as the signatures for predicting colorectal cancer recurrence in multiple datasets that were collected independently and profiled on different microarray platforms. The multi-gene modules we identified have a significant enrichment of known genes and biological processes relevant to cancer development, including genes from the chemokine pathway. Most strikingly, they recruited a significant enrichment of somatic mutations found in colorectal cancer. These results confirmed the functional relevance of these modules for colorectal cancer development. Further, these functional modules from different datasets overlapped significantly. Finally, we demonstrated that, leveraging above information of these modules, our module based classifier avoided arbitrary fitting the classifier function and screening the signatures using the training data, and achieved more consistency in prognosis prediction across three independent datasets, which holds even using very small training sets of tumors.     

Identifying breast cancer subtypes associated modules and biomarkers by integrated bioinformatics analysis

Breast cancer is the most common form of cancer afflicting women worldwide. Patients with breast cancer of different molecular classifications need varied treatments. Since it is known that the development of breast cancer involves multiple genes and functions, identification of functional gene modules (clusters of the functionally related genes) is indispensable as opposed to isolated genes, in order to investigate their relationship derived from the gene co-expression analysis. In total, 6315 differentially expressed genes (DEGs) were recognized and subjected to the co-expression analysis. Seven modules were screened out. The blue and turquoise modules have been selected from the module trait association analysis since the genes in these two modules are significantly correlated with the breast cancer subtypes. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment show that the blue module genes engaged in cell cycle, DNA replication, p53 signaling pathway, and pathway in cancer. According to the connectivity analysis and survival analysis, 8 out of 96 hub genes were filtered and have shown the highest expression in basal-like breast cancer. Furthermore, the hub genes were validated by the external datasets and quantitative real-time PCR (qRT-PCR). In summary, hub genes of Cyclin E1 (CCNE1), Centromere Protein N (CENPN), Checkpoint kinase 1 (CHEK1), Polo-like kinase 1 (PLK1), DNA replication and sister chromatid cohesion 1 (DSCC1), Family with sequence similarity 64, member A (FAM64A), Ubiquitin Conjugating Enzyme E2 C (UBE2C) and Ubiquitin Conjugating Enzyme E2 T (UBE2T) may serve as the prognostic markers for different subtypes of breast cancer.     

A Network of Cancer Genes with Co-Occurring and Anti-Co-Occurring Mutations

Certain cancer genes contribute to tumorigenesis in a manner of either co-occurring or mutually exclusive (anti-co-occurring) mutations; however, the global picture of when, where and how these functional interactions occur remains unclear. This study presents a systems biology approach for this purpose. After applying this method to cancer gene mutation data generated from large-scale and whole genome sequencing of cancer samples, a network of cancer genes with co-occurring and anti-co-occurring mutations was constructed. Analysis of this network revealed that genes with co-occurring mutations prefer direct signaling transductions and that the interaction relations among cancer genes in the network are related with their functional similarity. It was also revealed that genes with co-occurring mutations tend to have similar mutation frequencies, whereas genes with anti-co-occurring mutations tend to have different mutation frequencies. Moreover, genes with more exons tend to have more co-occurring mutations with other genes, and genes having lower local coherent network structures tend to have higher mutation frequency. The network showed two complementary modules that have distinct functions and have different roles in tumorigenesis. This study presented a framework for the analysis of cancer genome sequencing outputs. The presented data and uncovered patterns are helpful for understanding the contribution of gene mutations to tumorigenesis and valuable in the identification of key biomarkers and drug targets for cancer.     

Exploring metabolic pathway disruption in the subchronic phencyclidine model of schizophrenia with the Generalized Singular Value Decomposition

Background

One of the major goals in gene and protein expression profiling of cancer is to identify biomarkers and build classification models for prediction of disease prognosis or treatment response. Many traditional statistical methods, based on microarray gene expression data alone and individual genes' discriminatory power, often fail to identify biologically meaningful biomarkers thus resulting in poor prediction performance across data sets. Nonetheless, the variables in multivariable classifiers should synergistically interact to produce more effective classifiers than individual biomarkers.

Results

We developed an integrated approach, namely network-constrained support vector machine (netSVM), for cancer biomarker identification with an improved prediction performance. The netSVM approach is specifically designed for network biomarker identification by integrating gene expression data and protein-protein interaction data. We first evaluated the effectiveness of netSVM using simulation studies, demonstrating its improved performance over state-of-the-art network-based methods and gene-based methods for network biomarker identification. We then applied the netSVM approach to two breast cancer data sets to identify prognostic signatures for prediction of breast cancer metastasis. The experimental results show that: (1) network biomarkers identified by netSVM are highly enriched in biological pathways associated with cancer progression; (2) prediction performance is much improved when tested across different data sets. Specifically, many genes related to apoptosis, cell cycle, and cell proliferation, which are hallmark signatures of breast cancer metastasis, were identified by the netSVM approach. More importantly, several novel hub genes, biologically important with many interactions in PPI network but often showing little change in expression as compared with their downstream genes, were also identified as network biomarkers; the genes were enriched in signaling pathways such as TGF-beta signaling pathway, MAPK signaling pathway, and JAK-STAT signaling pathway. These signaling pathways may provide new insight to the underlying mechanism of breast cancer metastasis.

Conclusions

We have developed a network-based approach for cancer biomarker identification, netSVM, resulting in an improved prediction performance with network biomarkers. We have applied the netSVM approach to breast cancer gene expression data to predict metastasis in patients. Network biomarkers identified by netSVM reveal potential signaling pathways associated with breast cancer metastasis, and help improve the prediction performance across independent data sets.