scholarly journals lionessR: single-sample network reconstruction in R

2019 ◽  
Author(s):  
Marieke L. Kuijjer ◽  
John Quackenbush ◽  
Kimberly Glass

SummaryWe recently developed LIONESS (Linear Interpolation to Obtain Network Estimates for Single Samples), a method that can be used together with network reconstruction algorithms to extract networks for individual samples in a population. LIONESS was originally made available as a function within the PANDA (Passing Attributes between Networks for Data Assimilation) regulatory network reconstruction framework. In this application note, we describe lionessR, an R implementation of LIONESS that can be applied to any network reconstruction method in R that outputs a complete, weighted adjacency matrix. As an example, we use lionessR to model single-sample co-expression networks on a bone cancer dataset, and show how lionessR can be used to identify differential co-expression between two groups of patients.Availability and implementationThe lionessR open source R package, which includes a vignette of the application, is freely available athttps://github.com/mararie/[email protected]

BMC Cancer ◽  
2019 ◽  
Vol 19 (1) ◽  
Author(s):  
Marieke L Kuijjer ◽  
Ping-Han Hsieh ◽  
John Quackenbush ◽  
Kimberly Glass

Abstract Background In biomedical research, network inference algorithms are typically used to infer complex association patterns between biological entities, such as between genes or proteins, using data from a population. This resulting aggregate network, in essence, averages over the networks of those individuals in the population. LIONESS (Linear Interpolation to Obtain Network Estimates for Single Samples) is a method that can be used together with a network inference algorithm to extract networks for individual samples in a population. The method’s key characteristic is that, by modeling networks for individual samples in a data set, it can capture network heterogeneity in a population. LIONESS was originally made available as a function within the PANDA (Passing Attributes between Networks for Data Assimilation) regulatory network reconstruction framework. However, the LIONESS algorithm is generalizable and can be used to model single sample networks based on a wide range of network inference algorithms. Results In this software article, we describe lionessR, an R implementation of LIONESS that can be applied to any network inference method in R that outputs a complete, weighted adjacency matrix. As an example, we provide a vignette of an application of lionessR to model single sample networks based on correlated gene expression in a bone cancer dataset. We show how the tool can be used to identify differential patterns of correlation between two groups of patients. Conclusions We developed lionessR, an open source R package to model single sample networks. We show how lionessR can be used to inform us on potential precision medicine applications in cancer. The lionessR package is a user-friendly tool to perform such analyses. The package, which includes a vignette describing the application, is freely available at: https://github.com/kuijjerlab/lionessR and at: http://bioconductor.org/packages/lionessR.


2021 ◽  
Author(s):  
Marieke Lydia Kuijjer ◽  
Kimberly Glass

We recently developed LIONESS, a method to estimate sample-specific networks based on the output of an aggregate network reconstruction approach. In this manuscript, we describe how to apply LIONESS to different network reconstruction algorithms and data types. We highlight how decisions related to data preprocessing may affect the output networks, discuss expected outcomes, and give examples of how to analyze and compare single sample networks.


2021 ◽  
Vol 5 (3) ◽  
pp. 83
Author(s):  
Bilgi Görkem Yazgaç ◽  
Mürvet Kırcı

In this paper, we propose a fractional differential equation (FDE)-based approach for the estimation of instantaneous frequencies for windowed signals as a part of signal reconstruction. This approach is based on modeling bandpass filter results around the peaks of a windowed signal as fractional differential equations and linking differ-integrator parameters, thereby determining the long-range dependence on estimated instantaneous frequencies. We investigated the performance of the proposed approach with two evaluation measures and compared it to a benchmark noniterative signal reconstruction method (SPSI). The comparison was provided with different overlap parameters to investigate the performance of the proposed model concerning resolution. An additional comparison was provided by applying the proposed method and benchmark method outputs to iterative signal reconstruction algorithms. The proposed FDE method received better evaluation results in high resolution for the noniterative case and comparable results with SPSI with an increasing iteration number of iterative methods, regardless of the overlap parameter.


2015 ◽  
Vol 2015 ◽  
pp. 1-5 ◽  
Author(s):  
Yang Cao ◽  
Xiaofei Zheng ◽  
Fei Li ◽  
Xiaochen Bo

The human microbiome plays important roles in human health and disease. Previous microbiome studies focused mainly on single pure species function and overlooked the interactions in the complex communities on system-level. A metagenomic approach introduced recently integrates metagenomic data with community-level metabolic network modeling, but no comprehensive tool was available for such kind of approaches. To facilitate these kinds of studies, we developed an R package,mmnet, to implement community-level metabolic network reconstruction. The package also implements a set of functions for automatic analysis pipeline construction including functional annotation of metagenomic reads, abundance estimation of enzymatic genes, community-level metabolic network reconstruction, and integrated network analysis. The result can be represented in an intuitive way and sent to Cytoscape for further exploration. The package has substantial potentials in metagenomic studies that focus on identifying system-level variations of human microbiome associated with disease.


2018 ◽  
Vol 11 (02) ◽  
pp. 1750014 ◽  
Author(s):  
Jingjing Yu ◽  
Qiyue Li ◽  
Haiyu Wang

Bioluminescence tomography (BLT) is an important noninvasive optical molecular imaging modality in preclinical research. To improve the image quality, reconstruction algorithms have to deal with the inherent ill-posedness of BLT inverse problem. The sparse characteristic of bioluminescent sources in spatial distribution has been widely explored in BLT and many L1-regularized methods have been investigated due to the sparsity-inducing properties of L1 norm. In this paper, we present a reconstruction method based on L[Formula: see text] regularization to enhance sparsity of BLT solution and solve the nonconvex L[Formula: see text] norm problem by converting it to a series of weighted L1 homotopy minimization problems with iteratively updated weights. To assess the performance of the proposed reconstruction algorithm, simulations on a heterogeneous mouse model are designed to compare it with three representative sparse reconstruction algorithms, including the weighted interior-point, L1 homotopy, and the Stagewise Orthogonal Matching Pursuit algorithm. Simulation results show that the proposed method yield stable reconstruction results under different noise levels. Quantitative comparison results demonstrate that the proposed algorithm outperforms the competitor algorithms in location accuracy, multiple-source resolving and image quality.


2020 ◽  
Author(s):  
Tobias Rubel ◽  
Anna Ritz

AbstractSignaling pathways drive cellular response, and understanding such pathways is fundamental to molecular systems biology. A mounting volume of experimental protein interaction data has motivated the development of algorithms to computationally reconstruct signaling pathways. However, existing methods suffer from low recall in recovering protein interactions in ground truth pathways, limiting our confidence in any new predictions for experimental validation. We present the Pathway Reconstruction AUGmenter (PRAUG), a higher-order function for producing high-quality pathway reconstruction algorithms. PRAUG modifies any existing pathway reconstruction method, resulting in augmented algorithms that outperform their un-augmented counterparts for six different algorithms across twenty-nine diverse signaling pathways. The algorithms produced by PRAUG collectively reveal potential new proteins and interactions involved in the Wnt and Notch signaling pathways. PRAUG offers a valuable framework for signaling pathway prediction and discovery.


2021 ◽  
Author(s):  
Roozbeh Manshaei

Understanding and interpreting the inherently uncertain nature of complex biological systems, as well as the time to an event in these systems, are notable challenges in the field of bioinformatics. Overcoming these challenges could potentially lead to scientific discoveries, for example paving the path for the design of new drugs to target specific diseases such as cancer, or helping to apply more effective treatment for these diseases. In general, reverse engineering of these types of biological systems using online datasets is difficult. In particular, finding a unique solution to these systems is hard due to their complexity and the small sample size of datasets. This remains an unsolved problem due to such uncertainty, and the often intractable solution space of these systems. The term"uncertainty" describes the application-based margin of significance, validity, and efficiency of inferred or predictive models in their ability to extract characteristic properties and features describing the observed state of a given biological system. In this work, uncertainties within two specific bioinformatics domains are considered, namely "gene regulatory network reconstruction" (in which gene interactions/relationships within a biological entity are inferred from gene expression data); and "cancer survivorship prediction" (in which patient survival rates are predicted based on clinical factors and treatment outcomes). One approach to reduce uncertainty is to apply different constraints that have particular relevance to each application domain. In gene network reconstruction for instance, the consideration of constraints such as sparsity, stability and modularity, can inform and reduce uncertainty in the inferred reconstructions. While in cancer survival prediction, there is uncertainty in determining which clinical features (or feature aggregates) can improve associated prediction models. The inherent lack of understanding of how, why and when such constraints should be applied, however, prompts the need for a radically new approach. In this dissertation, a new approach is thus considered to aid human expert users in understanding and exploring inherent uncertainties associated with these two bioinformatics domains. Specifically, a novel set of tools is introduced and developed to assist in evidence gathering, constraint definition, and refinement of models toward the discovery of better solutions. This dissertation employs computational approaches, including convex optimization and feature selection/aggregation, in order to increase the chances of finding a unique solution. These approaches are realized through three novel interactive tools that employ tangible interaction in combination with graphical visualization to enable experts to query and manipulate the data. Tangible interaction provides physical embodiments of data and computational functions in support of learning and collaboration. Using these approaches, the dissertation demonstrates: (1) a modified stability constraint for reconstructing gene regulatory network that shows improvement in accuracy of predicted networks, (2) a novel modularity constraint (neighbor norm) for extracting available structures in the data which is validated with Laplacian eigenvalue spectrum, and (3) a hybrid method for estimating overall survival and inferring effective prognosis factors for patients with advanced prostate cancer that improves the accuracy of survival analysis.


Sign in / Sign up

Export Citation Format

Share Document