A closer look at "social" boundary genes reveals knowledge to gene expression profiles

As social network analysis is gaining popularity in modeling real world problems, the task of applying the social network model concepts and notions to biological data is still one of the most attractive research problems to be addressed. According, our work described in this paper focuses on a particular set of genes that reside on the community boundaries in gene co-expression networks. Stemmed from community mining problem in social networks, peripheries of communities (i.e., boundaries) can be used to aid certain biological analysis. The proposed method consists of three parts: 1) Finding communities of gene co-expression networks through clustering. 2) Analyzing stability of community structures by Monte Carlo method. 3) Designing of dynamic adoption of boundaries using geometric convexity. We validated our findings using breast cancer gene expression data from various studies. Our approach contributes to the new branch of applying social network mechanisms in biological data analysis, leading to new data mining strategies implied by witnessing social behaviors in gene expression analysis.     

Recurrent neural network based hybrid model for reconstructing gene regulatory network

One of the exciting problems in systems biology research is to decipher how genome controls the development of complex biological system. The gene regulatory networks (GRNs) help in the identification of regulatory interactions between genes and offer fruitful information related to functional role of individual gene in a cellular system. Discovering GRNs lead to a wide range of applications, including identification of disease related pathways providing novel tentative drug targets, helps to predict disease response, and also assists in diagnosing various diseases including cancer. Reconstruction of GRNs from available biological data is still an open problem. This paper proposes a recurrent neural network (RNN) based model of GRN, hybridized with generalized extended Kalman filter for weight update in backpropagation through time training algorithm. The RNN is a complex neural network that gives a better settlement between biological closeness and mathematical flexibility to model GRN; and is also able to capture complex, non-linear and dynamic relationships among variables. Gene expression data are inherently noisy and Kalman filter performs well for estimation problem even in noisy data. Hence, we applied non-linear version of Kalman filter, known as generalized extended Kalman filter, for weight update during RNN training. The developed model has been tested on four benchmark networks such as DNA SOS repair network, IRMA network, and two synthetic networks from DREAM Challenge. We performed a comparison of our results with other state-of-the-art techniques which shows superiority of our proposed model. Further, 5% Gaussian noise has been induced in the dataset and result of the proposed model shows negligible effect of noise on results, demonstrating the noise tolerance capability of the model.     

Node-based differential network analysis in genomics

Gene dependency networks often undergo changes in response to different conditions. Understanding how these networks change across two conditions is an important task in genomics research. Most previous differential network analysis approaches assume that the difference between two condition-specific networks is driven by individual edges. Thus, they may fail in detecting key players which might represent important genes whose mutations drive the change of network. In this work, we develop a node-based differential network analysis (N-DNA) model to directly estimate the differential network that is driven by certain hub nodes. We model each condition-specific gene network as a precision matrix and the differential network as the difference between two precision matrices. Then we formulate a convex optimization problem to infer the differential network by combing a D-trace loss function and a row-column overlap norm penalty function. Simulation studies demonstrate that N-DNA provides more accurate estimate of the differential network than previous competing approaches. We apply N-DNA to ovarian cancer and breast cancer gene expression data. The model rediscovers known cancer-related genes and contains interesting predictions.     

Perturbations of pathway co-expression network identify a core network in metastatic breast cancer

Metastases are the main cause of death in advanced breast cancer (BC) patients. Although chemotherapy and hormone therapy are current treatment strategies, drug resistance is frequent and still not completely understood.In this study, a bioinformatics analysis was performed on BC patients to explore the molecular mechanisms associated with BC metastasis. Microarray gene expression profiles of metastatic and non metastatic BC patients were downloaded from Gene Expression Omnibus (GEO) dataset. Raw data were normalized and merged using the Combat tool. Pathways enriched with differently expressed genes were identified and a pathway co-expression network was generated using Pearson’s correlation. We identified from this network, which includes 17 pathways and 128 interactions, the pathways that most influence the network efficiency. Moreover, protein interaction network was investigated to identify hub genes of the pathway network. The prognostic role of the network was evaluated with a survival analysis using an independent dataset.In conclusion, the pathway co-expression network could contribute to understanding the mechanism and development of BC metastases.     

Characterizing mutation–expression network relationships in multiple cancers

BackgroundData made available through large cancer consortia like The Cancer Genome Atlas make for a rich source of information to be studied across and between cancers. In recent years, network approaches have been applied to such data in uncovering the complex interrelationships between mutational and expression profiles, but lack direct testing for expression changes via mutation. In this pan-cancer study we analyze mutation and gene expression information in an integrative manner by considering the networks generated by testing for differences in expression in direct association with specific mutations. We relate our findings among the 19 cancers examined to identify commonalities and differences as well as their characteristics.ResultsUsing somatic mutation and gene expression information across 19 cancers, we generated mutation–expression networks per cancer. On evaluation we found that our generated networks were significantly enriched for known cancer-related genes, such as skin cutaneous melanoma (p < 0.01 using Network of Cancer Genes 4.0). Our framework identified that while different cancers contained commonly mutated genes, there was little concordance between associated gene expression changes among cancers. Comparison between cancers showed a greater overlap of network nodes for cancers with higher overall non-silent mutation load, compared to those with a lower overall non-silent mutation load.ConclusionsThis study offers a framework that explores network information through co-analysis of somatic mutations and gene expression profiles. Our pan-cancer application of this approach suggests that while mutations are frequently common among cancer types, the impact they have on the surrounding networks via gene expression changes varies. Despite this finding, there are some cancers for which mutation-associated network behaviour appears to be similar: suggesting a potential framework for uncovering related cancers for which similar therapeutic strategies may be applicable. Our framework for understanding relationships among cancers has been integrated into an interactive R Shiny application, PAn Cancer Mutation Expression Networks (PACMEN), containing dynamic and static network visualization of the mutation–expression networks. PACMEN also features tools for further examination of network topology characteristics among cancers.     

Identifying cancer specific signaling pathways based on the dysregulation between genes

A large collection of studies has shown that the occurrence of cancer is related to the functional dysfunction of the pathways. Identification of cancer-related pathways could help researchers understand the mechanisms of complex diseases well. Whereas, most current signaling pathway analysis methods take no account of the gene interaction variations within pathways. Furthermore, considering that some pathways have connection with two or more cancer types, while some are likely to be cancer-type specific pathways. Identifying cancer-type specific pathways contributes to interpreting the different mechanisms of different cancer types. In this study, we first proposed a pathway analysis method named Pathway Analysis of Intergenic Regulation (PAIGR) to identify pathways with dysregulation between genes and compared the performance of this method with four existing methods on four colorectal cancer (CRC) datasets. The results showed that PAIGR could find cancer-related pathways more accurately. Moreover, in order to explore the relationship between the identified pathways and the cancer type, we constructed a pathway interaction network, in which nodes and edges represented pathways and interactions between pathways respectively. Highly connected pathways were considered to play a central role in an extensive range of biological processes, while sparsely connected pathways are considered to have certain specificity. Our results showed that pathways identified by PAIGR had a low nodal degree (i.e., a few numbers of interactions), which suggested that most of these pathways were cancer-type specific.     

NFκB pathway analysis: An approach to analyze gene co-expression networks employing feedback cycles

The genes of the NFκB pathway are involved in the control of a plethora of biological processes ranking from inhibition of apoptosis to metastasis in cancer. It has been described that Gliobastoma multiforme (GBM) patients carry aberrant NFκB activation, but the molecular mechanisms are not completely understood. Here, we present a NFκB pathway analysis in tumor specimens of GBM compared to non-neoplasic brain tissues, based on the different kind of cycles found among genes of a gene co-expression network constructed using quantized data obtained from the microarrays. A cycle is a closed walk with all vertices distinct (except the first and last). Thanks to this way of finding relations among genes, a more robust interpretation of gene correlations is possible, because the cycles are associated with feedback mechanisms that are very common in biological networks. In GBM samples, we could conclude that the stoichiometric relationship between genes involved in NFκB pathway regulation is unbalanced. This can be measured and explained by the identification of a cycle. This conclusion helps to understand more about the biology of this type of tumor.     

Gene network and canonical pathway analysis in prostate cancer: a microarray study

The molecular mechanism playing a role in the development of prostate cancer (PCA) is not well defined. We decided to determine the changes in gene expression in PCA tissues and to compare them to those in non-cancerous samples. Prostate tissue samples were collected by needle biopsy from 21 PCA and 10 benign prostate hyperplasic (BPH) patients. Total RNA was isolated, cDNA was synthesized, and gene expression levels were determined by microarray method. In the progression to PCA, 738 up-regulated and 515 down-regulated genes were detected in samples. Analysis using Ingenuity Pathway Analysis (IPA) software revealed that 466 network and 423 functions-pathways eligible genes were up-regulated, and 363 network and 342 functions-pathways eligible genes were down-regulated. Up-regulated networks were identified around IL-1beta and insulin-like growth factor-1 (IGF-1) genes. The NFKB gene was centered around two up- and down-regulated networks. Up-regulated canonical pathways were assigned and four of them were evaluated in detail: acute phase response, hepatic fibrosis, actin cytoskeleton, and coagulation pathways. Axonal guidance signaling was the most significant down-regulated canonical pathway. Our data provide not only networks between the genes for understanding the biologic properties of PCA but also useful pathway maps for future understanding of disease and the construction of new therapeutic targets.     

Gene selection of rat hepatocyte proliferation using adaptive sparse group lasso with weighted gene co-expression network analysis

Grouped gene selection is the most important task for analyzing the microarray data of rat liver regeneration. Many existing gene selection methods cannot outstand the interactions among the selected genes. In the process of rat liver regeneration, one of the most important events involved in many biological processes is the proliferation of rat hepatocytes, so it can be used as a measure of the effectiveness of the method. Here we proposed an adaptive sparse group lasso to select genes in groups for rat hepatocyte proliferation. The weighted gene co-expression networks analysis was used to identify modules corresponding to gene pathways, based on which a strategy of dividing genes into groups was proposed. A strategy of adaptive gene selection was also presented by assessing the gene significance and introducing the adaptive lasso penalty. Moreover, an improved blockwise descent algorithm was proposed. Experimental results demonstrated that the proposed method can improve the classification accuracy, and select less number of significant genes which act jointly in groups and have direct or indirect effects on rat hepatocyte proliferation. The effectiveness of the method was verified by the method of rat hepatocyte proliferation.     

Inferring gene regulatory networks from temporal expression profiles under time-delay and noise

Ordinary differential equations (ODE) have been widely used for modeling and analysis of dynamic gene networks in systems biology. In this paper, we propose an optimization method that can infer a gene regulatory network from time-series gene expression data. Specifically, the following four cases are considered: (1) reconstruction of a gene network from synthetic gene expression data with noise, (2) reconstruction of a gene network from synthetic gene expression data with time-delay, (3) reconstruction of a gene network from synthetic gene expression data with noise and time-delay, and (4) reconstruction of a gene network from experimental time-series data in budding yeast cell cycle.     

Genome-wide identification and characterization of long non-coding RNAs involved in acquired resistance to gefitinib in non-small-cell lung cancer

Acquired resistance is a major obstacle to the therapeutic efficacy of gefitinib in non-small-cell lung cancer (NSCLC). Current knowledge about the role of long non-coding RNAs (lncRNAs) in this phenomenon is insufficient. In this study, we searched RNA sequencing data for lncRNAs associated with acquired resistance to gefitinib in NSCLC, and constructed a functional lncRNA-mRNA co-expression network and protein-protein interaction (PPI) network to analyze their putative target genes and biological functions. The expression levels of 14 outstanding dysregulated lncRNAs and mRNA were verified using real-time PCR. Changes in the expression levels of 39 lncRNAs and 121 mRNAs showed common patterns in our two pairs of gefitinib-sensitive and gefitinib-resistant NSCLC cell lines. The co-expression network included 1235 connections among these common differentially expressed lncRNAs and mRNAs. The significantly enriched signaling pathways based on dysregulated mRNAs were mainly involved in the Hippo signaling pathway; proteoglycans in cancer; and valine, leucine, and isoleucine biosynthesis. The results show that LncRNAs play an important part in acquired gefitinib resistance in NSCLC by regulating mRNA expression and function, and may represent potential new molecular biomarkers and therapeutic targets for gefitinib-resistant NSCLC.     

Network regularised Cox regression and multiplex network models to predict disease comorbidities and survival of cancer

In cancer genomics, gene expression levels provide important molecular signatures for all types of cancer, and this could be very useful for predicting the survival of cancer patients. However, the main challenge of gene expression data analysis is high dimensionality, and microarray is characterised by few number of samples with large number of genes. To overcome this problem, a variety of penalised Cox proportional hazard models have been proposed. We introduce a novel network regularised Cox proportional hazard model and a novel multiplex network model to measure the disease comorbidities and to predict survival of the cancer patient. Our methods are applied to analyse seven microarray cancer gene expression datasets: breast cancer, ovarian cancer, lung cancer, liver cancer, renal cancer and osteosarcoma. Firstly, we applied a principal component analysis to reduce the dimensionality of original gene expression data. Secondly, we applied a network regularised Cox regression model on the reduced gene expression datasets. By using normalised mutual information method and multiplex network model, we predict the comorbidities for the liver cancer based on the integration of diverse set of omics and clinical data, and we find the diseasome associations (disease–gene association) among different cancers based on the identified common significant genes. Finally, we evaluated the precision of the approach with respect to the accuracy of survival prediction using ROC curves. We report that colon cancer, liver cancer and renal cancer share the CXCL5 gene, and breast cancer, ovarian cancer and renal cancer share the CCND2 gene. Our methods are useful to predict survival of the patient and disease comorbidities more accurately and helpful for improvement of the care of patients with comorbidity. Software in Matlab and R is available on our GitHub page: https://github.com/ssnhcom/NetworkRegularisedCox.git.     

A network-based pharmacology study of active compounds and targets of Fritillaria thunbergii against influenza

Seasonal and pandemic influenza infections are serious threats to public health and the global economy. Since antigenic drift reduces the effectiveness of conventional therapies against the virus, herbal medicine has been proposed as an alternative. Fritillaria thunbergii (FT) have been traditionally used to treat airway inflammatory diseases such as coughs, bronchitis, pneumonia, and fever-based illnesses. Herein, we used a network pharmacology-based strategy to predict potential compounds from Fritillaria thunbergii (FT), target genes, and cellular pathways to better combat influenza and influenza-associated diseases. We identified five compounds, and 47 target genes using a compound-target network (C-T). Two compounds (beta-sitosterol and pelargonidin) and nine target genes (BCL2, CASP3, HSP90AA1, ICAM1, JUN, NOS2, PPARG, PTGS1, PTGS2) were identified using a compound-influenza disease target network (C-D). Protein-protein interaction (PPI) network was constructed and we identified eight proteins from nine target genes formed a network. The compound-disease-pathway network (C-D-P) revealed three classes of pathways linked to influenza: cancer, viral diseases, and inflammation. Taken together, our systems biology data from C-T, C-D, PPI and C-D-P networks predicted potent compounds from FT and new therapeutic targets and pathways involved in influenza.