Weighted Fused Pathway Graphical Lasso for Joint Estimation of Multiple Gene Networks

An exploratory study of CT radiomics using differential network feature selection for WHO/ISUP grading and progression-free survival prediction of clear cell renal cell carcinoma

Objectives

To explore the feasibility of predicting the World Health Organization/International Society of Urological Pathology (WHO/ISUP) grade and progression-free survival (PFS) of clear cell renal cell cancer (ccRCC) using the radiomics features (RFs) based on the differential network feature selection (FS) method using the maximum-entropy probability model (MEPM).

Methods

175 ccRCC patients were divided into a training set (125) and a test set (50). The non-contrast phase (NCP), cortico-medullary phase, nephrographic phase, excretory phase phases, and all-phase WHO/ISUP grade prediction models were constructed based on a new differential network FS method using the MEPM. The diagnostic performance of the best phase model was compared with the other state-of-the-art machine learning models and the clinical models. The RFs of the best phase model were used for survival analysis and visualized using risk scores and nomograms. The performance of the above models was tested in both cross-validated and independent validation and checked by the Hosmer-Lemeshow test.

Results

The NCP RFs model was the best phase model, with an AUC of 0.89 in the test set, and performed superior to other machine learning models and the clinical models (all p <0.05). Kaplan-Meier survival analysis, univariate and multivariate cox regression results, and risk score analyses showed the NCP RFs could predict PFS well (almost all p < 0.05). The nomogram model incorporated the best two RFs and showed good discrimination, a C-index of 0.71 and 0.69 in the training and test set, and good calibration.

Conclusion

The NCP CT-based RFs selected by differential network FS could predict the WHO/ISUP grade and PFS of RCC.

Detecting abnormal connectivity in schizophrenia via a joint directed acyclic graph estimation model

Functional connectivity (FC) between brain region has been widely studied and linked with cognition and behavior of an individual. FC is usually defined as the correlation or partial correlation of fMRI blood oxygen level-dependent (BOLD) signals between two brain regions. Although FC has been effective to understand brain organization, it cannot reveal the direction of interactions. Many directed acyclic graph (DAG) based methods have been applied to study the directed interactions but their performance was limited by the small sample size while high dimensionality of the available data. By enforcing group regularization and utilizing samples from both case and control groups, we propose a joint DAG model to estimate the directed FC. We first demonstrate that the proposed model is efficient and accurate through a series of simulation studies. We then apply it to the case-control study of schizophrenia (SZ) with data collected from the MIND Clinical Imaging Consortium (MCIC). We have successfully identified decreased functional integration, disrupted hub structures and characteristic edges (CtEs) in SZ patients. Those findings have been confirmed by previous studies with some identified to be potential markers for SZ patients. A comparison of the results between the directed FC and undirected FC showed substantial differences in the selected features. In addition, we used the identified features based on directed FC for the classification of SZ patients and achieved better accuracy than using undirected FC or raw features, demonstrating the advantage of using directed FC for brain network analysis.

Inferring Differential Networks by Integrating Gene Expression Data With Additional Knowledge

Evidences increasingly indicate the involvement of gene network rewiring in disease development and cell differentiation. With the accumulation of high-throughput gene expression data, it is now possible to infer the changes of gene networks between two different states or cell types via computational approaches. However, the distribution diversity of multi-platform gene expression data and the sparseness and high noise rate of single-cell RNA sequencing (scRNA-seq) data raise new challenges for existing differential network estimation methods. Furthermore, most existing methods are purely rely on gene expression data, and ignore the additional information provided by various existing biological knowledge. In this study, to address these challenges, we propose a general framework, named weighted joint sparse penalized D-trace model (WJSDM), to infer differential gene networks by integrating multi-platform gene expression data and multiple prior biological knowledge. Firstly, a non-paranormal graphical model is employed to tackle gene expression data with missing values. Then we propose a weighted group bridge penalty to integrate multi-platform gene expression data and various existing biological knowledge. Experiment results on synthetic data demonstrate the effectiveness of our method in inferring differential networks. We apply our method to the gene expression data of ovarian cancer and the scRNA-seq data of circulating tumor cells of prostate cancer, and infer the differential network associated with platinum resistance of ovarian cancer and anti-androgen resistance of prostate cancer. By analyzing the estimated differential networks, we find some important biological insights about the mechanisms underlying platinum resistance of ovarian cancer and anti-androgen resistance of prostate cancer.

A novel probabilistic generator for large-scale gene association networks

MOTIVATION

Gene expression data provide an opportunity for reverse-engineering gene-gene associations using network inference methods. However, it is difficult to assess the performance of these methods because the true underlying network is unknown in real data. Current benchmarks address this problem by subsampling a known regulatory network to conduct simulations. But the topology of regulatory networks can vary greatly across organisms or tissues, and reference-based generators-such as GeneNetWeaver-are not designed to capture this heterogeneity. This means, for example, benchmark results from the E. coli regulatory network will not carry over to other organisms or tissues. In contrast, probabilistic generators do not require a reference network, and they have the potential to capture a rich distribution of topologies. This makes probabilistic generators an ideal approach for obtaining a robust benchmarking of network inference methods.

RESULTS

We propose a novel probabilistic network generator that (1) provides an alternative to address the inherent limitation of reference-based generators and (2) is able to create realistic gene association networks, and (3) captures the heterogeneity found across gold-standard networks better than existing generators used in practice. Eight organism-specific and 12 human tissue-specific gold-standard association networks are considered. Several measures of global topology are used to determine the similarity of generated networks to the gold-standards. Along with demonstrating the variability of network structure across organisms and tissues, we show that the commonly used "scale-free" model is insufficient for replicating these structures.

AVAILABILITY

This generator is implemented in the R package "SeqNet" and is available on CRAN (https://cran.r-project.org/web/packages/SeqNet/index.html).

Estimation of Gene Regulatory Networks from Cancer Transcriptomics Data

Cancer is a genetic disease in which multiple genes are perturbed. Thus, information about the regulatory relationships between genes is necessary for the identification of biomarkers and therapeutic targets. In this review, methods for inference of gene regulatory networks (GRNs) from transcriptomics data that are used in cancer research are introduced. The methods are classified into three categories according to the analysis model. The first category includes methods that use pair-wise measures between genes, including correlation coefficient and mutual information. The second category includes methods that determine the genetic regulatory relationship using multivariate measures, which consider the expression profiles of all genes concurrently. The third category includes methods using supervised and integrative approaches. The supervised approach estimates the regulatory relationship using a supervised learning method that constructs a regression or classification model for predicting whether there is a regulatory relationship between genes with input data of gene expression profiles and class labels of prior biological knowledge. The integrative method is an expansion of the supervised method and uses more data and biological knowledge for predicting the regulatory relationship. Furthermore, simulation and experimental validation of the estimated GRNs are also discussed in this review. This review identified that most GRN inference methods are not specific for cancer transcriptome data, and such methods are required for better understanding of cancer pathophysiology. In addition, more systematic methods for validation of the estimated GRNs need to be developed in the context of cancer biology.

Stabilized Reconstruction of Signaling Networks from Single-Cell Cue-Response Data

Inferring cell-signaling networks from high-throughput data is a challenging problem in systems biology. Recent advances in cytometric technology enable us to measure the abundance of a large number of proteins at the single-cell level across time. Traditional network reconstruction approaches usually consider each time point separately, resulting thus in inferred networks that strongly vary across time. To account for the possibly time-invariant physical couplings within the signaling network, we extend the traditional graphical lasso with an additional regularizer that penalizes network variations over time. ROC evaluation of the method on in silico data showed higher reconstruction accuracy than standard graphical lasso. We also tested our approach on single-cell mass cytometry data of IFNγ-stimulated THP1 cells with 26 phospho-proteins simultaneously measured. Our approach recapitulated known signaling relationships, such as connection within the JAK/STAT pathway, and was further validated in characterizing perturbed signaling network with PI3K, MEK1/2 and AMPK inhibitors.