An order independent algorithm for inferring gene regulatory network using quantile value for conditional independence tests

CNNGRN: A Convolutional Neural Network-Based Method for Gene Regulatory Network Inference From Bulk Time-Series Expression Data

Gene regulatory networks (GRNs) participate in many biological processes, and reconstructing them plays an important role in systems biology. Although many advanced methods have been proposed for GRN reconstruction, their predictive performance is far from the ideal standard, so it is urgent to design a more effective method to reconstruct GRN. Moreover, most methods only consider the gene expression data, ignoring the network structure information contained in GRN. In this study, we propose a supervised model named CNNGRN, which infers GRN from bulk time-series expression data via convolutional neural network (CNN) model, with a more informative feature. Bulk time series gene expression data imply the intricate regulatory associations between genes, and the network structure feature of ground-truth GRN contains rich neighbor information. Hence, CNNGRN integrates the above two features as model inputs. In addition, CNN is adopted to extract intricate features of genes and infer the potential associations between regulators and target genes. Moreover, feature importance visualization experiments are implemented to seek the key features. Experimental results show that CNNGRN achieved competitive performance on benchmark datasets compared to the state-of-the-art computational methods. Finally, hub genes identified based on CNNGRN have been confirmed to be involved in biological processes through literature.

Gene Regulatory Network-Classifier: Gene Regulatory Network-Based Classifier and Its Applications to Gastric Cancer Drug (5-Fluorouracil) Marker Identification

The complex mechanisms of diseases involve the disturbance of the molecular network, rather than disorder in a single gene, implying that single gene-based analysis is insufficient to understand these mechanisms. Gene regulatory networks (GRNs) have attracted a lot of interest and various approaches have been developed for their statistical inference and gene network-based analysis. Although various computational methods have been developed, relatively little attention has been paid to incorporation of biological knowledge into the computational approaches. Furthermore, existing studies on network-based analysis perform prediction/classification of status of cell lines based on preconstructed GRNs, implying that we cannot extract prediction/classification-specific gene networks, leading to difficulty in interpretation of biological mechanisms and marker identification related to the status of cancer cell lines. We developed a novel strategy to build a GRN-based classifier, called a GRN-classifier. The proposed GRN-classifier estimates GRNs and classifies cell lines simultaneously, where the gene network is estimated to minimize error in gene network estimation and the negative log-likelihood for classifying cell lines. Thus, we can identify biological status-specific gene regulatory systems, enabling us to achieve biologically reliable interpretation of the classification. We also propose an algorithm to implement the GRN-classifier based on coordinate descent update. Monte Carlo simulations were conducted to examine performance of the GRN-classifier. Results: Our strategy provides effective results in feature selection in the classification model and edge selection in gene network estimation. The GRN-classifier also shows outstanding classification accuracy. We apply the GRN-classifier to classify cancer cell lines into anticancer drug-related status, that is, 5-fluorouracil (5-FU)-sensitive/resistant and 5-FU target/nontarget cancer cell lines. We then identified 5-FU markers based on 5-FU-related status classification-specific gene networks. The mechanisms of the identified markers were verified through literature survey. Our results suggest that the molecular interplay between MYOF and AHNAK2 may play a crucial role in drug resistance and can provide information on the chemotherapy efficiency of 5-FU. It is also suggested that suppression of the identified 5-FU markers, including MYOF/AHNAK2 and AKR1C1/AKR1C3 may improve 5-FU resistance of cancer cell lines.

An approach of gene regulatory network construction using mixed entropy optimizing context-related likelihood mutual information

MOTIVATION

The question of how to construct gene regulatory networks has long been a focus of biological research. Mutual information can be used to measure nonlinear relationships, and it has been widely used in the construction of gene regulatory networks. However, this method cannot measure indirect regulatory relationships under the influence of multiple genes, which reduces the accuracy of inferring gene regulatory networks.

APPROACH

This work proposes a method for constructing gene regulatory networks based on mixed entropy optimizing context-related likelihood mutual information (MEOMI). First, two entropy estimators were combined to calculate the mutual information between genes. Then, distribution optimization was performed using a context-related likelihood algorithm to eliminate some indirect regulatory relationships and obtain the initial gene regulatory network. To obtain the complex interaction between genes and eliminate redundant edges in the network, the initial gene regulatory network was further optimized by calculating the conditional mutual inclusive information (CMI2) between gene pairs under the influence of multiple genes. The network was iteratively updated to reduce the impact of mutual information on the overestimation of the direct regulatory intensity.

RESULTS

The experimental results show that the MEOMI method performed better than several other kinds of gene network construction methods on DREAM challenge simulated datasets (DREAM3 and DREAM5), three real Escherichia coli datasets (E.coli SOS pathway network, E.coli SOS DNA repair network and E.coli community network) and two human datasets.

AVAILABILITY AND IMPLEMENTATION

Source code and dataset are available at https://github.com/Dalei-Dalei/MEOMI/ and http://122.205.95.139/MEOMI/.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

A decomposition structure learning algorithm in Bayesian network based on a two-stage combination method

Decomposition hybrid algorithms with the recursive framework which recursively decompose the structural task into structural subtasks to reduce computational complexity are employed to learn Bayesian network (BN) structure. Merging rules are commonly adopted as the combination method in the combination step. The direction determination rule of merging rules has problems in using the idea of keeping v-structures unchanged before and after combination to determine directions of edges in the whole structure. It breaks down in one case due to appearances of wrong v-structures, and is hard to operate in practice. Therefore, we adopt a novel approach for direction determination and propose a two-stage combination method. In the first-stage combination method, we determine nodes, links of edges by merging rules and adopt the idea of permutation and combination to determine directions of contradictory edges. In the second-stage combination method, we restrict edges between nodes that do not satisfy the decomposition property and their parent nodes by determining the target domain according to the decomposition property. Simulation experiments on four networks show that the proposed algorithm can obtain BN structure with higher accuracy compared with other algorithms. Finally, the proposed algorithm is applied to the thickening process of gold hydrometallurgy to solve the practical problem.