Review of biological network data and its applications

Special Issue on “Biological Network Approaches and Applications”

Biological phenomena comprise various interactions between genes and molecules [...]

Intersection of network medicine and machine learning towards investigating the key biomarkers and pathways underlying amyotrophic lateral sclerosis: a systematic review

BACKGROUND

Network medicine is an emerging area of research that focuses on delving into the molecular complexity of the disease, leading to the discovery of network biomarkers and therapeutic target discovery. Amyotrophic lateral sclerosis (ALS) is a complicated rare disease with unknown pathogenesis and no available treatment. In ALS, network properties appear to be potential biomarkers that can be beneficial in disease-related applications when explored independently or in tandem with machine learning (ML) techniques.

OBJECTIVE

This systematic literature review explores recent trends in network medicine and implementations of network-based ML algorithms in ALS. We aim to provide an overview of the identified primary studies and gather details on identifying the potential biomarkers and delineated pathways.

METHODS

The current study consists of searching for and investigating primary studies from PubMed and Dimensions.ai, published between 2018 and 2022 that reported network medicine perspectives and the coupling of ML techniques. Each abstract and full-text study was individually evaluated, and the relevant studies were finally included in the review for discussion once they met the inclusion and exclusion criteria.

RESULTS

We identified 109 eligible publications from primary studies representing this systematic review. The data coalesced into two themes: application of network science to identify disease modules and promising biomarkers in ALS, along with network-based ML approaches. Conclusion This systematic review gives an overview of the network medicine approaches and implementations of network-based ML algorithms in ALS to determine new disease genes, and identify critical pathways and therapeutic target discovery for personalized treatment.

Basics on network theory to analyze biological systems: a hands-on outlook

Biological processes result from interactions among molecules and cell-to-cell communications. In the last 50 years, network theory has empowered advances in understanding molecular networks' structure and dynamics that regulate biological systems. Adopting a network data analysis point of view at more laboratories might enrich their research capacity to generate forward working hypotheses. This work briefly describes network theory origins and provides basic graph analysis principles in biological systems, specific centrality measurements, and the main models for network structures. Also, we describe a workflow employing user-friendly free platforms to process, construct, and analyze transcriptome data from a network perspective. With this assay, we expect to encourage the implementation of network theory analysis on biological data in everyday laboratory research.

Modulation of RNA splicing associated with Wnt signaling pathway using FD-895 and pladienolide B

Alterations in RNA splicing are associated with different malignancies, including leukemia, lymphoma, and solid tumors. The RNA splicing modulators such as FD-895 and pladienolide B have been investigated in different malignancies to target/modulate spliceosome for therapeutic purpose. Different cell lines were screened using an RNA splicing modulator to test in vitro cytotoxicity and the ability to modulate RNA splicing capability via induction of intron retention (using RT-PCR and qPCR). The Cignal Finder Reporter Array evaluated [pathways affected by the splice modulators in HeLa cells. Further, the candidates associated with the pathways were validated at protein level using western blot assay, and gene-gene interaction studies were carried out using GeneMANIA. We show that FD-895 and pladienolide B induces higher apoptosis levels than conventional chemotherapy in different solid tumors. In addition, both agents modulate Wnt signaling pathways and mRNA splicing. Specifically, FD-895 and pladienolide B significantly downregulates Wnt signaling pathway-associated transcripts (GSK3β and LRP5) and both transcript and proteins including LEF1, CCND1, LRP6, and pLRP6 at the transcript, total protein, and protein phosphorylation's levels. These results indicate FD-895 and pladienolide B inhibit Wnt signaling by decreasing LRP6 phosphorylation and modulating mRNA splicing through induction of intron retention in solid tumors.

Computational analyses of mechanism of action (MoA): data, methods and integration

The elucidation of a compound's Mechanism of Action (MoA) is a challenging task in the drug discovery process, but it is important in order to rationalise phenotypic findings and to anticipate potential side-effects. Bioinformatic approaches, advances in machine learning techniques and the increasing deposition of high-throughput data in public databases have significantly contributed to recent advances in the field, but it is not straightforward to decide which data and methods are most suitable to use in a given case. In this review, we focus on these methods and data and their applications in generating MoA hypotheses for subsequent experimental validation. We discuss compound-specific data such as -omics, cell morphology and bioactivity data, as well as commonly used supplementary prior knowledge such as network and pathway data, and provide information on databases where this data can be accessed. In terms of methodologies, we discuss both well-established methods (connectivity mapping, pathway enrichment) as well as more developing methods (neural networks and multi-omics integration). Finally, we review case studies where the MoA of a compound was successfully suggested from computational analysis by incorporating multiple data modalities and/or methodologies. Our aim for this review is to provide researchers with insights into the benefits and drawbacks of both the data and methods in terms of level of understanding, biases and interpretation - and to highlight future avenues of investigation which we foresee will improve the field of MoA elucidation, including greater public access to -omics data and methodologies which are capable of data integration.

Network theory reveals principles of spliceosome structure and dynamics

Cryoelectron microscopy has revolutionized spliceosome structural biology, and structures representing much of the splicing process have been determined. Comparison of these structures is challenging due to extreme dynamics of the splicing machinery and the thousands of changing interactions during splicing. We have used network theory to analyze splicing factor interactions by constructing structure-based networks from protein-protein, protein-RNA, and RNA-RNA interactions found in eight different spliceosome structures. Our analyses reveal that connectivity dynamics result in step-specific impacts of factors on network topology. The spliceosome's connectivity is focused on the active site, in part due to contributions from nonglobular proteins. Many essential factors exhibit large shifts in centralities during splicing. Others show transiently high betweenness centralities at certain stages, thereby suggesting mechanisms for regulating splicing by briefly bridging otherwise poorly connected network nodes. These observations provide insights into organizing principles of the spliceosome and provide frameworks for comparative analysis of other macromolecular machines.

IMEx Databases: Displaying Molecular Interactions into a Single, Standards-Compliant Dataset

Molecular interaction databases aim to systematically capture and organize the experimental interaction information described in the scientific literature. These data can then be used to perform network analysis, to assign putative roles to uncharacterized proteins and to investigate their involvement in cellular pathways.This chapter gives a brief overview of publicly available molecular interaction databases and focuses on the members of the IMEx Consortium, on their curation policies and standard data formats. All of the goals achieved by IMEx databases over the last 15 years, the data types provided and the many different ways in which such data can be utilized by the research community, are described in detail. The IMEx databases curate molecular interaction data to the highest caliber, following a detailed curation model and supplying rich metadata by employing common curation rules and harmonized standards. The IMEx Consortium provides comprehensively annotated molecular interaction data integrated into a single, non-redundant, open access dataset.

Novel cancer subtyping method based on patient-specific gene regulatory network

The identification of cancer subtypes is important for the understanding of tumor heterogeneity. In recent years, numerous computational methods have been proposed for this problem based on the multi-omics data of patients. It is widely accepted that different cancer subtypes are induced by different molecular regulatory networks. However, only a few incorporate the differences between their molecular systems into the identification processes. In this study, we present a novel method to identify cancer subtypes based on patient-specific molecular systems. Our method realizes this by quantifying patient-specific gene networks, which are estimated from their transcriptome data, and by clustering their quantified networks. Comprehensive analyses of The Cancer Genome Atlas (TCGA) datasets applied to our method confirmed that they were able to identify more clinically meaningful cancer subtypes than the existing subtypes and found that the identified subtypes comprised different molecular features. Our findings also show that the proposed method can identify the novel cancer subtypes even with single omics data, which cannot otherwise be captured by existing methods using multi-omics data.

Phenotype correlations reveal the relationships of physiological systems underlying human ageing

Ageing is characterized by degeneration and loss of function across multiple physiological systems. To study the mechanisms and consequences of ageing, several metrics have been proposed in a hierarchical model, including biological, phenotypic and functional ageing. In particular, phenotypic ageing and interconnected changes in multiple physiological systems occur in all ageing individuals over time. Recently, phenotypic age, a new ageing measure, was proposed to capture morbidity and mortality risk across diverse subpopulations in US cohort studies. Although phenotypic age has been widely used, it may overlook the complex relationships among phenotypic biomarkers. Considering the correlation structure of these phenotypic biomarkers, we proposed a composite phenotype analysis (CPA) strategy to analyse 71 biomarkers from 2074 individuals in the Rugao Longitudinal Ageing Study. CPA grouped these biomarkers into 18 composite phenotypes according to their internal correlation, and these composite phenotypes were mostly consistent with prior findings. In addition, compared with prior findings, this strategy exhibited some different yet important implications. For example, the indicators of kidney and cardiovascular functions were tightly connected, implying internal interactions. The composite phenotypes were further verified through associations with functional metrics of ageing, including disability, depression, cognitive function and frailty. Compared to age alone, these composite phenotypes had better predictive performances for functional metrics of ageing. In summary, CPA could reveal the hidden relationships of physiological systems and identify the links between physiological systems and functional ageing metrics, thereby providing novel insights into potential mechanisms underlying human ageing.

Learning representation for multiple biological networks via a robust graph regularized integration approach

Learning node representation is a fundamental problem in biological network analysis, as compact representation features reveal complicated network structures and carry useful information for downstream tasks such as link prediction and node classification. Recently, multiple networks that profile objects from different aspects are increasingly accumulated, providing the opportunity to learn objects from multiple perspectives. However, the complex common and specific information across different networks pose challenges to node representation methods. Moreover, ubiquitous noise in networks calls for more robust representation. To deal with these problems, we present a representation learning method for multiple biological networks. First, we accommodate the noise and spurious edges in networks using denoised diffusion, providing robust connectivity structures for the subsequent representation learning. Then, we introduce a graph regularized integration model to combine refined networks and compute common representation features. By using the regularized decomposition technique, the proposed model can effectively preserve the common structural property of different networks and simultaneously accommodate their specific information, leading to a consistent representation. A simulation study shows the superiority of the proposed method on different levels of noisy networks. Three network-based inference tasks, including drug-target interaction prediction, gene function identification and fine-grained species categorization, are conducted using representation features learned from our method. Biological networks at different scales and levels of sparsity are involved. Experimental results on real-world data show that the proposed method has robust performance compared with alternatives. Overall, by eliminating noise and integrating effectively, the proposed method is able to learn useful representations from multiple biological networks.

CeTF: an R/Bioconductor package for transcription factor co-expression networks using regulatory impact factors (RIF) and partial correlation and information (PCIT) analysis

BACKGROUND

Finding meaningful gene-gene interaction and the main Transcription Factors (TFs) in co-expression networks is one of the most important challenges in gene expression data mining.

RESULTS

Here, we developed the R package "CeTF" that integrates the Partial Correlation with Information Theory (PCIT) and Regulatory Impact Factors (RIF) algorithms applied to gene expression data from microarray, RNA-seq, or single-cell RNA-seq platforms. This approach allows identifying the transcription factors most likely to regulate a given network in different biological systems - for example, regulation of gene pathways in tumor stromal cells and tumor cells of the same tumor. This pipeline can be easily integrated into the high-throughput analysis. To demonstrate the CeTF package application, we analyzed gastric cancer RNA-seq data obtained from TCGA (The Cancer Genome Atlas) and found the HOXB3 gene as the second most relevant TFs with a high regulatory impact (TFs-HRi) regulating gene pathways in the cell cycle.

CONCLUSION

This preliminary finding shows the potential of CeTF to list master regulators of gene networks. CeTF was designed as a user-friendly tool that provides many highly automated functions without requiring the user to perform many complicated processes. It is available on Bioconductor ( http://bioconductor.org/packages/CeTF ) and GitHub ( http://github.com/cbiagii/CeTF ).

Relationship between gene regulation network structure and prediction accuracy in high dimensional regression

The least absolute shrinkage and selection operator (lasso) and principal component regression (PCR) are popular methods of estimating traits from high-dimensional omics data, such as transcriptomes. The prediction accuracy of these estimation methods is highly dependent on the covariance structure, which is characterized by gene regulation networks. However, the manner in which the structure of a gene regulation network together with the sample size affects prediction accuracy has not yet been sufficiently investigated. In this study, Monte Carlo simulations are conducted to investigate the prediction accuracy for several network structures under various sample sizes. When the gene regulation network is a random graph, a sufficiently large number of observations are required to ensure good prediction accuracy with the lasso. The PCR provided poor prediction accuracy regardless of the sample size. However, a real gene regulation network is likely to exhibit a scale-free structure. In such cases, the simulation indicates that a relatively small number of observations, such as [Formula: see text], is sufficient to allow the accurate prediction of traits from a transcriptome with the lasso.

CLEP: A Hybrid Data- and Knowledge- Driven Framework for Generating Patient Representations

Summary

As machine learning and artificial intelligence increasingly attain a larger number of applications in the biomedical domain, at their core, their utility depends on the data used to train them. Due to the complexity and high dimensionality of biomedical data, there is a need for approaches that combine prior knowledge around known biological interactions with patient data. Here, we present CLinical Embedding of Patients (CLEP), a novel approach that generates new patient representations by leveraging both prior knowledge and patient-level data. First, given a patient-level dataset and a knowledge graph containing relations across features that can be mapped to the dataset, CLEP incorporates patients into the knowledge graph as new nodes connected to their most characteristic features. Next, CLEP employs knowledge graph embedding models to generate new patient representations that can ultimately be used for a variety of downstream tasks, ranging from clustering to classification. We demonstrate how using new patient representations generated by CLEP significantly improves performance in classifying between patients and healthy controls for a variety of machine learning models, as compared to the use of the original transcriptomics data. Furthermore, we also show how incorporating patients into a knowledge graph can foster the interpretation and identification of biological features characteristic of a specific disease or patient subgroup. Finally, we released CLEP as an open source Python package together with examples and documentation.

Availability and implementation

CLEP is available to the bioinformatics community as an open source Python package at https://github.com/hybrid-kg/clep under the Apache 2.0 License.

Supplementary information

Supplementary data are available at Bioinformatics online.

Probabilistic model by Bayesian network for the prediction of antibody glycosylation in perfusion and fed-batch cell cultures

Glycosylation is a critical quality attribute of therapeutic monoclonal antibodies (mAbs). The glycan pattern can have a large impact on the immunological functions, serum half-life and stability. The medium components and cultivation parameters are known to potentially influence the glycosylation profile. Mathematical modelling provides a strategy for rational design and control of the upstream bioprocess. However, the kinetic models usually contain a very large number of unknown parameters, which limit their practical applications. In this article, we consider the metabolic network of N-linked glycosylation as a Bayesian network (BN) and calculate the fluxes of the glycosylation process as joint probability using the culture parameters as inputs. The modelling approach is validated with data of different Chinese hamster ovary cell cultures in pseudo perfusion, perfusion, and fed batch cultures, all showing very good predictive capacities. In cases where a large number of cultivation parameters is available, it is shown here that principal components analysis can efficiently be employed for a dimension reduction of the inputs compared to Pearson correlation analysis and feature importance by decision tree. The present study demonstrates that BN model can be a powerful tool in upstream process and medium development for glycoprotein productions.

SIGNAL: A web-based iterative analysis platform integrating pathway and network approaches optimizes hit selection from genome-scale assays

Hit selection from high-throughput assays remains a critical bottleneck in realizing the potential of omic-scale studies in biology. Widely used methods such as setting of cutoffs, prioritizing pathway enrichments, or incorporating predicted network interactions offer divergent solutions yet are associated with critical analytical trade-offs. The specific limitations of these individual approaches and the lack of a systematic way by which to integrate their rankings have contributed to limited overlap in the reported results from comparable genome-wide studies and costly inefficiencies in secondary validation efforts. Using comparative analysis of parallel independent studies as a benchmark, we characterize the specific complementary contributions of each approach and demonstrate an optimal framework to integrate these methods. We describe selection by iterative pathway group and network analysis looping (SIGNAL), an integrated, iterative approach that uses both pathway and network methods to optimize gene prioritization. SIGNAL is accessible as a rapid user-friendly web-based application (https://signal.niaid.nih.gov). A record of this paper's transparent peer review is included in the Supplemental information.

Acute Inflammation in Traumatic Brain Injury and Polytrauma Patients Using Network Analysis

INTRODUCTION

Traumatic brain injury (TBI) is associated with secondary injury to the central nervous system (CNS) via inflammatory mechanisms. The combination of polytrauma and TBI further exacerbates the inflammatory response to injury; however, combined injury phenomena have not been thoroughly studied. In this study, we examined the inflammatory differences between patients with TBI versus patients with polytrauma, but no TBI (polytrauma). We hypothesize that patients with TBI have a heightened early inflammatory response compared with polytrauma.

METHODS

We conducted a single-center retrospective study of a cohort of patients with polytrauma, who were enrolled in the PROPPR study. These patients had blood samples prospectively collected at eight time points in the first 3 days of admission. Using radiological data to determine TBI, our polytrauma cohort was dichotomized into TBI (n = 30) or polytrauma (n = 54). Inflammatory biomarkers were measured using ELISA. Data across time were compared for TBI versus polytrauma groups using Wilcoxon rank-sum test. Network analysis techniques were used to systematically characterize the inflammatory responses at admission.

RESULTS

Patients with TBI (51.6%) had a higher 30-day mortality compared with polytrauma (16.9%) (P <0.001). Expression levels of IL6, IL8, and CCL2 were elevated from the 2-h through 24-h time points, becoming significant at the 6-h time point (IL6, IL8, and CCL2; P <0.05) (). CSF3 showed a similar pattern, but did not attain significance. TBI and polytrauma networks underwent diverging trends from admission to the 6-h time point.

CONCLUSION

Patients with TBI demonstrated upregulations in proinflammatory cytokines IL6, IL8, and CCL2. Utilizing informatics methods, we were able to identify temporal differences in network trends, as well as uncharacterized cytokines and chemokines in TBI. These data suggest TBI induces a distinct inflammatory response and pathologically heightened inflammatory response in the presence of polytrauma and may propagate worsened patient outcomes including mortality.

An Efficient and Easy-to-Use Network-Based Integrative Method of Multi-Omics Data for Cancer Genes Discovery

Identifying personalized driver genes is essential for discovering critical biomarkers and developing effective personalized therapies of cancers. However, few methods consider weights for different types of mutations and efficiently distinguish driver genes over a larger number of passenger genes. We propose MinNetRank (Minimum used for Network-based Ranking), a new method for prioritizing cancer genes that sets weights for different types of mutations, considers the incoming and outgoing degree of interaction network simultaneously, and uses minimum strategy to integrate multi-omics data. MinNetRank prioritizes cancer genes among multi-omics data for each sample. The sample-specific rankings of genes are then integrated into a population-level ranking. When evaluating the accuracy and robustness of prioritizing driver genes, our method almost always significantly outperforms other methods in terms of precision, F1 score, and partial area under the curve (AUC) on six cancer datasets. Importantly, MinNetRank is efficient in discovering novel driver genes. SP1 is selected as a candidate driver gene only by our method (ranked top three), and SP1 RNA and protein differential expression between tumor and normal samples are statistically significant in liver hepatocellular carcinoma. The top seven genes stratify patients into two subtypes exhibiting statistically significant survival differences in five cancer types. These top seven genes are associated with overall survival, as illustrated by previous researchers. MinNetRank can be very useful for identifying cancer driver genes, and these biologically relevant marker genes are associated with clinical outcome. The R package of MinNetRank is available at https://github.com/weitinging/MinNetRank.

Detecting Common Pathways and Key Molecules of Neurodegenerative Diseases from the Topology of Molecular Networks

MotivationNeurodegenerative diseases (NDs), including amyotrophic lateral sclerosis, Parkinson's disease, Alzheimer's disease, and Huntington's disease, occur as a result of neurodegenerative processes. Thus, it has been increasingly appreciated that many neurodegenerative conditions overlap at multiple levels. However, traditional clinicopathological correlation approaches to better classify a disease have met with limited success. Discovering this overlap offers hope for therapeutic advances that could ameliorate many ND simultaneously. In parallel, in the last decade, systems biology approaches have become a reliable choice in complex disease analysis for gaining more delicate biological insights and have enabled the comprehension of the higher order functions of the biological systems.ResultsToward this orientation, we developed a systems biology approach for the identification of common links and pathways of ND, based on well-established and novel topological and functional measures. For this purpose, a molecular pathway network was constructed, using molecular interactions and relations of four main neurodegenerative diseases (Alzheimer's disease, Parkinson's disease, amyotrophic lateral sclerosis, and Huntington's disease). Our analysis captured the overlapped subregions forming molecular subpathways fully enriched in these four NDs. Also, it exported molecules that act as bridges, hubs, and key players for neurodegeneration concerning either their topology or their functional role.ConclusionUnderstanding these common links and central topologies under the perspective of systems biology and network theory and greater insights are provided to uncover the complex neurodegeneration processes.

A biochemical network modeling of a whole-cell

All cellular processes can be ultimately understood in terms of respective fundamental biochemical interactions between molecules, which can be modeled as networks. Very often, these molecules are shared by more than one process, therefore interconnecting them. Despite this effect, cellular processes are usually described by separate networks with heterogeneous levels of detail, such as metabolic, protein-protein interaction, and transcription regulation networks. Aiming at obtaining a unified representation of cellular processes, we describe in this work an integrative framework that draws concepts from rule-based modeling. In order to probe the capabilities of the framework, we used an organism-specific database and genomic information to model the whole-cell biochemical network of the Mycoplasma genitalium organism. This modeling accounted for 15 cellular processes and resulted in a single component network, indicating that all processes are somehow interconnected. The topological analysis of the network showed structural consistency with biological networks in the literature. In order to validate the network, we estimated gene essentiality by simulating gene deletions and compared the results with experimental data available in the literature. We could classify 212 genes as essential, being 95% of them consistent with experimental results. Although we adopted a relatively simple organism as a case study, we suggest that the presented framework has the potential for paving the way to more integrated studies of whole organisms leading to a systemic analysis of cells on a broader scale. The modeling of other organisms using this framework could provide useful large-scale models for different fields of research such as bioengineering, network biology, and synthetic biology, and also provide novel tools for medical and industrial applications.

Turning genome-wide association study findings into opportunities for drug repositioning

Drug development is a very costly and lengthy process, while repositioned or repurposed drugs could be brought into clinical practice within a shorter time-frame and at a much reduced cost. Numerous computational approaches to drug repositioning have been developed, but methods utilizing genome-wide association studies (GWASs) data are less explored.

The past decade has observed a massive growth in the amount of data from GWAS; the rich information contained in GWAS has great potential to guide drug repositioning or discovery. While multiple tools are available for finding the most relevant genes from GWAS hits, searching for top susceptibility genes is only one way to guide repositioning, which has its own limitations.

Here we provide a comprehensive review of different computational approaches that employ GWAS data to guide drug repositioning. These methods include selecting top candidate genes from GWAS as drug targets, deducing drug candidates based on drug-drug and disease-disease similarities, searching for reversed expression profiles between drugs and diseases, pathway-based methods as well as approaches based on analysis of biological networks. Each method is illustrated with examples, and their respective strengths and limitations are discussed. We also discussed several areas for future research.

Molecular networks in Network Medicine: Development and applications

Network Medicine applies network science approaches to investigate disease pathogenesis. Many different analytical methods have been used to infer relevant molecular networks, including protein-protein interaction networks, correlation-based networks, gene regulatory networks, and Bayesian networks. Network Medicine applies these integrated approaches to Omics Big Data (including genetics, epigenetics, transcriptomics, metabolomics, and proteomics) using computational biology tools and, thereby, has the potential to provide improvements in the diagnosis, prognosis, and treatment of complex diseases. We discuss briefly the types of molecular data that are used in molecular network analyses, survey the analytical methods for inferring molecular networks, and review efforts to validate and visualize molecular networks. Successful applications of molecular network analysis have been reported in pulmonary arterial hypertension, coronary heart disease, diabetes mellitus, chronic lung diseases, and drug development. Important knowledge gaps in Network Medicine include incompleteness of the molecular interactome, challenges in identifying key genes within genetic association regions, and limited applications to human diseases. This article is categorized under: Models of Systems Properties and Processes > Mechanistic Models Translational, Genomic, and Systems Medicine > Translational Medicine Analytical and Computational Methods > Analytical Methods Analytical and Computational Methods > Computational Methods.

Target Identification Using Homopharma and Network-Based Methods for Predicting Compounds Against Dengue Virus-Infected Cells

Drug target prediction is an important method for drug discovery and design, can disclose the potential inhibitory effect of active compounds, and is particularly relevant to many diseases that have the potential to kill, such as dengue, but lack any healing agent. An antiviral drug is urgently required for dengue treatment. Some potential antiviral agents are still in the process of drug discovery, but the development of more effective active molecules is in critical demand. Herein, we aimed to provide an efficient technique for target prediction using homopharma and network-based methods, which is reliable and expeditious to hunt for the possible human targets of three phenolic lipids (anarcardic acid, cardol, and cardanol) related to dengue viral (DENV) infection as a case study. Using several databases, the similarity search and network-based analyses were applied on the three phenolic lipids resulting in the identification of seven possible targets as follows. Based on protein annotation, three phenolic lipids may interrupt or disturb the human proteins, namely KAT5, GAPDH, ACTB, and HSP90AA1, whose biological functions have been previously reported to be involved with viruses in the family Flaviviridae. In addition, these phenolic lipids might inhibit the mechanism of the viral proteins: NS3, NS5, and E proteins. The DENV and human proteins obtained from this study could be potential targets for further molecular optimization on compounds with a phenolic lipid core structure in anti-dengue drug discovery. As such, this pipeline could be a valuable tool to identify possible targets of active compounds.

Similarities in biological processes can be used to bridge ecology and molecular biology

Much of the research in biology aims to understand the origin of diversity. Naturally, ecological diversity was the first object of study, but we now have the necessary tools to probe diversity at molecular scales. The inherent differences in how we study diversity at different scales caused the disciplines of biology to be organized around these levels, from molecular biology to ecology. Here, we illustrate that there are key properties of each scale that emerge from the interactions of simpler components and that these properties are often shared across different levels of organization. This means that ideas from one level of organization can be an inspiration for novel hypotheses to study phenomena at another level. We illustrate this concept with examples of events at the molecular level that have analogs at the organismal or ecological level and vice versa. Through these examples, we illustrate that biological processes at different organization levels are governed by general rules. The study of the same phenomena at different scales could enrich our work through a multidisciplinary approach, which should be a staple in the training of future scientists.

scNPF: an integrative framework assisted by network propagation and network fusion for preprocessing of single-cell RNA-seq data

Background

Single-cell RNA-sequencing (scRNA-seq) is fast becoming a powerful tool for profiling genome-scale transcriptomes of individual cells and capturing transcriptome-wide cell-to-cell variability. However, scRNA-seq technologies suffer from high levels of technical noise and variability, hindering reliable quantification of lowly and moderately expressed genes. Since most downstream analyses on scRNA-seq, such as cell type clustering and differential expression analysis, rely on the gene-cell expression matrix, preprocessing of scRNA-seq data is a critical preliminary step in the analysis of scRNA-seq data.

Results

We presented scNPF, an integrative scRNA-seq preprocessing framework assisted by network propagation and network fusion, for recovering gene expression loss, correcting gene expression measurements, and learning similarities between cells. scNPF leverages the context-specific topology inherent in the given data and the priori knowledge derived from publicly available molecular gene-gene interaction networks to augment gene-gene relationships in a data driven manner. We have demonstrated the great potential of scNPF in scRNA-seq preprocessing for accurately recovering gene expression values and learning cell similarity networks. Comprehensive evaluation of scNPF across a wide spectrum of scRNA-seq data sets showed that scNPF achieved comparable or higher performance than the competing approaches according to various metrics of internal validation and clustering accuracy. We have made scNPF an easy-to-use R package, which can be used as a versatile preprocessing plug-in for most existing scRNA-seq analysis pipelines or tools.

Conclusions

scNPF is a universal tool for preprocessing of scRNA-seq data, which jointly incorporates the global topology of priori interaction networks and the context-specific information encapsulated in the scRNA-seq data to capture both shared and complementary knowledge from diverse data sources. scNPF could be used to recover gene signatures and learn cell-to-cell similarities from emerging scRNA-seq data to facilitate downstream analyses such as dimension reduction, cell type clustering, and visualization.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5747-5) contains supplementary material, which is available to authorized users.

Fast Subgraph Matching Strategies Based on Pattern-Only Heuristics

Many scientific applications entail solving the subgraph isomorphism problem, i.e., given an input pattern graph, find all the subgraphs of a (usually much larger) target graph that are structurally equivalent to that input. Because subgraph isomorphism is NP-complete, methods to solve it have to use heuristics. This work evaluates subgraph isomorphism methods to assess their computational behavior on a wide range of synthetic and real graphs. Surprisingly, our experiments show that, among the leading algorithms, certain heuristics based only on pattern graphs are the most efficient.

Integrating Multiple Interaction Networks for Gene Function Inference

In the past few decades, the number and variety of genomic and proteomic data available have increased dramatically. Molecular or functional interaction networks are usually constructed according to high-throughput data and the topological structure of these interaction networks provide a wealth of information for inferring the function of genes or proteins. It is a widely used way to mine functional information of genes or proteins by analyzing the association networks. However, it remains still an urgent but unresolved challenge how to combine multiple heterogeneous networks to achieve more accurate predictions. In this paper, we present a method named ReprsentConcat to improve function inference by integrating multiple interaction networks. The low-dimensional representation of each node in each network is extracted, then these representations from multiple networks are concatenated and fed to gcForest, which augment feature vectors by cascading and automatically determines the number of cascade levels. We experimentally compare ReprsentConcat with a state-of-the-art method, showing that it achieves competitive results on the datasets of yeast and human. Moreover, it is robust to the hyperparameters including the number of dimensions.

A guide to gene regulatory network inference for obtaining predictive solutions: Underlying assumptions and fundamental biological and data constraints

The study of biological systems at a system level has become a reality due to the increasing powerful computational approaches able to handle increasingly larger datasets. Uncovering the dynamic nature of gene regulatory networks in order to attain a system level understanding and improve the predictive power of biological models is an important research field in systems biology. The task itself presents several challenges, since the problem is of combinatorial nature and highly depends on several biological constraints and also the intended application. Given the intrinsic interdisciplinary nature of gene regulatory network inference, we present a review on the currently available approaches, their challenges and limitations. We propose guidelines to select the most appropriate method considering the underlying assumptions and fundamental biological and data constraints.

Construction of a Transcriptome-Driven Network at the Early Stage of Infection with Influenza A H1N1 in Human Lung Alveolar Epithelial Cells

We aimed to understand the molecular changes in host cells that accompany infection by the seasonal influenza A H1N1 virus because the initial response rapidly changes owing to the fact that the virus has a robust initial propagation phase. Human epithelial alveolar A549 cells were infected and total RNA was extracted at 30 min, 1 h, 2 h, 4 h, 8 h, 24 h, and 48 h post infection (h.p.i.). The differentially expressed host genes were clustered into two distinct sets of genes as the infection progressed over time. The patterns of expression were significantly different at the early stages of infection. One of the responses showed roles similar to those associated with the enrichment gene sets to known 'gp120 pathway in HIV.' This gene set contains genes known to play roles in preventing the progress of apoptosis, which infected cells undergo as a response to viral infection. The other gene set showed enrichment of 'Drug Metabolism Enzymes (DMEs).' The identification of two distinct gene sets indicates that the virus regulates the cell's mechanisms to create a favorable environment for its stable replication and protection of gene metabolites within 8 h.

Bipartite graphs in systems biology and medicine: a survey of methods and applications

The latest advances in high-throughput techniques during the past decade allowed the systems biology field to expand significantly. Today, the focus of biologists has shifted from the study of individual biological components to the study of complex biological systems and their dynamics at a larger scale. Through the discovery of novel bioentity relationships, researchers reveal new information about biological functions and processes. Graphs are widely used to represent bioentities such as proteins, genes, small molecules, ligands, and others such as nodes and their connections as edges within a network. In this review, special focus is given to the usability of bipartite graphs and their impact on the field of network biology and medicine. Furthermore, their topological properties and how these can be applied to certain biological case studies are discussed. Finally, available methodologies and software are presented, and useful insights on how bipartite graphs can shape the path toward the solution of challenging biological problems are provided.

Systematic Evaluation of Molecular Networks for Discovery of Disease Genes

Gene networks are rapidly growing in size and number, raising the question of which networks are most appropriate for particular applications. Here, we evaluate 21 human genome-wide interaction networks for their ability to recover 446 disease gene sets identified through literature curation, gene expression profiling, or genome-wide association studies. While all networks have some ability to recover disease genes, we observe a wide range of performance with STRING, ConsensusPathDB, and GIANT networks having the best performance overall. A general tendency is that performance scales with network size, suggesting that new interaction discovery currently outweighs the detrimental effects of false positives. Correcting for size, we find that the DIP network provides the highest efficiency (value per interaction). Based on these results, we create a parsimonious composite network with both high efficiency and performance. This work provides a benchmark for selection of molecular networks in human disease research.

Biological network border detection

Complex networks have been widely used to model biological systems. The concept of accessibility has been proposed recently as a means to organize the nodes of complex networks as belonging to its border or center. Such an approach paves the way to investigating how the functional and structural properties of nodes vary with their respective position in the networks. In this work, we approach such a problem in a biological context applying border detection to Protein-Protein Interaction networks from four organisms of the Mycoplasma genus. We found evidence that the borderness of proteins bears a relation with their spatial organization and molecular function specificity.

NFP: An R Package for Characterizing and Comparing of Annotated Biological Networks

Large amounts of various biological networks exist for representing different types of interaction data, such as genetic, metabolic, gene regulatory, and protein-protein relationships. Recent approaches on biological network study are based on different mathematical concepts. It is necessary to construct a uniform framework to judge the functionality of biological networks. We recently introduced a knowledge-based computational framework that reliably characterized biological networks in system level. The method worked by making systematic comparisons to a set of well-studied "basic networks," measuring both the functional and topological similarities. A biological network could be characterized as a spectrum-like vector consisting of similarities to basic networks. Here, to facilitate the application, development, and adoption of this framework, we present an R package called NFP. This package extends our previous pipeline, offering a powerful set of functions for Network Fingerprint analysis. The software shows great potential in biological network study. The open source NFP R package is freely available under the GNU General Public License v2.0 at CRAN along with the vignette.

Graph-theoretical model of global human interactome reveals enhanced long-range communicability in cancer networks

Malignant transformation is known to involve substantial rearrangement of the molecular genetic landscape of the cell. A common approach to analysis of these alterations is a reductionist one and consists of finding a compact set of differentially expressed genes or associated signaling pathways. However, due to intrinsic tumor heterogeneity and tissue specificity, biomarkers defined by a small number of genes/pathways exhibit substantial variability. As an alternative to compact differential signatures, global features of genetic cell machinery are conceivable. Global network descriptors suggested in previous works are, however, known to potentially be biased by overrepresentation of interactions between frequently studied genes-proteins. Here, we construct a cellular network of 74538 directional and differential gene expression weighted protein-protein and gene regulatory interactions, and perform graph-theoretical analysis of global human interactome using a novel, degree-independent feature—the normalized total communicability (NTC). We apply this framework to assess differences in total information flow between different cancer (BRCA/COAD/GBM) and non-cancer interactomes. Our experimental results reveal that different cancer interactomes are characterized by significant enhancement of long-range NTC, which arises from circulation of information flow within robustly organized gene subnetworks. Although enhancement of NTC emerges in different cancer types from different genomic profiles, we identified a subset of 90 common genes that are related to elevated NTC in all studied tumors. Our ontological analysis shows that these genes are associated with enhanced cell division, DNA replication, stress response, and other cellular functions and processes typically upregulated in cancer. We conclude that enhancement of long-range NTC manifested in the correlated activity of genes whose tight coordination is required for survival and proliferation of all tumor cells, and, thus, can be seen as a graph-theoretical equivalent to some hallmarks of cancer. The computational framework for differential network analysis presented herein is of potential interest for a wide range of network perturbation problems given by single or multiple gene-protein activation-inhibition.

Seasonal Dynamics and Metagenomic Characterization of Marine Viruses in Goseong Bay, Korea

Viruses are the most abundant biological entities in the oceans, and account for a significant amount of the genetic diversity of marine ecosystems. However, there is little detailed information about the biodiversity of viruses in marine environments. Rapid advances in metagenomics have enabled the identification of previously unknown marine viruses. We performed metagenomic profiling of seawater samples collected at 6 sites in Goseong Bay (South Sea, Korea) during the spring, summer, autumn, and winter of 2014. The results indicated the presence of highly diverse virus communities. The DNA libraries from samples collected during four seasons were sequenced using Illumina HiSeq 2000. The number of viral reads was 136,850 during March, 70,651 during June, 66,165 during September, and 111,778 during December. Species identification indicated that Pelagibacter phage HTVC010P, Ostreococcus lucimarinus OIV5 and OIV1, and Roseobacter phage SIO1 were the most common species in all samples. For viruses with at least 10 reads, there were 204 species during March, 189 during June, 170 during September, and 173 during December. Analysis of virus families indicated that the Myoviridae was the most common during all four seasons, and viruses in the Polyomaviridae were only present during March. Viruses in the Iridoviridae were only present during three seasons. Additionally, viruses in the Iridoviridae, Herpesviridae, and Poxviridae, which may affect fish and marine animals, appeared during different seasons. These results suggest that seasonal changes in temperature contribute to the dynamic structure of the viral community in the study area. The information presented here will be useful for comparative analyses with other marine viral communities.

Computational Methods for Annotation Transfers from Sequence

Surveys of public sequence resources show that experimentally supported functional information is still completely missing for a considerable fraction of known proteins and is clearly incomplete for an even larger portion. Bioinformatics methods have long made use of very diverse data sources alone or in combination to predict protein function, with the understanding that different data types help elucidate complementary biological roles. This chapter focuses on methods accepting amino acid sequences as input and producing GO term assignments directly as outputs; the relevant biological and computational concepts are presented along with the advantages and limitations of individual approaches.

On the Variable Ordering in Subgraph Isomorphism Algorithms

Graphs are mathematical structures to model several biological data. Applications to analyze them require to apply solutions for the subgraph isomorphism problem, which is NP-complete. Here, we investigate the existing strategies to reduce the subgraph isomorphism algorithm running time with emphasis on the importance of the order with which the graph vertices are taken into account during the search, called variable ordering, and its incidence on the total running time of the algorithms. We focus on two recent solutions, which are based on an effective variable ordering strategy. We discuss their comparison both with the variable ordering strategies reviewed in the paper and the other algorithms present in the ICPR2014 contest on graph matching algorithms for pattern search in biological databases.

Compact Integration of Multi-Network Topology for Functional Analysis of Genes

The topological landscape of molecular or functional interaction networks provides a rich source of information for inferring functional patterns of genes or proteins. However, a pressing yet-unsolved challenge is how to combine multiple heterogeneous networks, each having different connectivity patterns, to achieve more accurate inference. Here, we describe the Mashup framework for scalable and robust network integration. In Mashup, the diffusion in each network is first analyzed to characterize the topological context of each node. Next, the high-dimensional topological patterns in individual networks are canonically represented using low-dimensional vectors, one per gene or protein. These vectors can then be plugged into off-the-shelf machine learning methods to derive functional insights about genes or proteins. We present tools based on Mashup that achieve state-of-the-art performance in three diverse functional inference tasks: protein function prediction, gene ontology reconstruction, and genetic interaction prediction. Mashup enables deeper insights into the structure of rapidly accumulating and diverse biological network data and can be broadly applied to other network science domains.

Correlation-Based Network Generation, Visualization, and Analysis as a Powerful Tool in Biological Studies: A Case Study in Cancer Cell Metabolism

In the last decade vast data sets are being generated in biological and medical studies. The challenge lies in their summary, complexity reduction, and interpretation. Correlation-based networks and graph-theory based properties of this type of networks can be successfully used during this process. However, the procedure has its pitfalls and requires specific knowledge that often lays beyond classical biology and includes many computational tools and software. Here we introduce one of a series of methods for correlation-based network generation and analysis using freely available software. The pipeline allows the user to control each step of the network generation and provides flexibility in selection of correlation methods and thresholds. The pipeline was implemented on published metabolomics data of a population of human breast carcinoma cell lines MDA-MB-231 under two conditions: normal and hypoxia. The analysis revealed significant differences between the metabolic networks in response to the tested conditions. The network under hypoxia had 1.7 times more significant correlations between metabolites, compared to normal conditions. Unique metabolic interactions were identified which could lead to the identification of improved markers or aid in elucidating the mechanism of regulation between distantly related metabolites induced by the cancer growth.

Molecular signatures of ovarian diseases: Insights from network medicine perspective

Dysfunctions and disorders in the ovary lead to a host of diseases including ovarian cancer, ovarian endometriosis, and polycystic ovarian syndrome (PCOS). Understanding the molecular mechanisms behind ovarian diseases is a great challenge. In the present study, we performed a meta-analysis of transcriptome data for ovarian cancer, ovarian endometriosis, and PCOS, and integrated the information gained from statistical analysis with genome-scale biological networks (protein-protein interaction, transcriptional regulatory, and metabolic). Comparative and integrative analyses yielded reporter biomolecules (genes, proteins, metabolites, transcription factors, and micro-RNAs), and unique or common signatures at protein, metabolism, and transcription regulation levels, which might be beneficial to uncovering the underlying biological mechanisms behind the diseases. These signatures were mostly associated with formation or initiation of cancer development, and pointed out the potential tendency of PCOS and endometriosis to tumorigenesis. Molecules and pathways related to MAPK signaling, cell cycle, and apoptosis were the mutual determinants in the pathogenesis of all three diseases. To our knowledge, this is the first report that screens these diseases from a network medicine perspective. This study provides signatures which could be considered as potential therapeutic targets and/or as medical prognostic biomarkers in further experimental and clinical studies. Abbreviations DAVID: Database for Annotation, Visualization and Integrated Discovery; DEGs: differentially expressed genes; GEO: Gene Expression Omnibus; KEGG: Kyoto Encyclopedia of Genes and Genomes; LIMMA: Linear Models for Microarray Data; MBRole: Metabolite Biological Role; miRNA: micro-RNA; PCOS: polycystic ovarian syndrome; PPI: protein-protein interaction; RMA: Robust Multi-Array Average; TF: transcription factor.

CyTRANSFINDER: a Cytoscape 3.3 plugin for three-component (TF, gene, miRNA) signal transduction pathway construction

Background

Biological research increasingly relies on network models to study complex phenomena. Signal Transduction Pathways are molecular circuits that model how cells receive, process, and respond to information from the environment providing snapshots of the overall cell dynamics. Most of the attempts to reconstruct signal transduction pathways are limited to single regulator networks including only genes/proteins. However, networks involving a single type of regulator and neglecting transcriptional and post-transcriptional regulations mediated by transcription factors and microRNAs, respectively, may not fully reveal the complex regulatory mechanisms of a cell. We observed a lack of computational instruments supporting explorative analysis on this type of three-component signal transduction pathways.

Results

We have developed CyTRANSFINDER, a new Cytoscape plugin able to infer three-component signal transduction pathways based on user defined regulatory patterns and including miRNAs, TFs and genes. Since CyTRANSFINDER has been designed to support exploratory analysis, it does not rely on expression data. To show the potential of the plugin we have applied it in a study of two miRNAs that are particularly relevant in human melanoma progression, miR-146a and miR-214.

Conclusions

CyTRANSFINDER supports the reconstruction of small signal transduction pathways among groups of genes. Results obtained from its use in a real case study have been analyzed and validated through both literature data and preliminary wet-lab experiments, showing the potential of this tool when performing exploratory analysis.

Electronic supplementary material

The online version of this article (doi:10.1186/s12859-016-0964-2) contains supplementary material, which is available to authorized users.

Tumor suppressor genes and their underlying interactions in paclitaxel resistance in cancer therapy

Objectives

Paclitaxel (PTX) is frequently used in the clinical treatment of solid tumors. But the PTX-resistance is a great obstacle in cancer treatment. Exploration of the mechanisms of drug resistance suggests that tumor suppressor genes (TSGs) play a key role in the response of chemotherapeutic drugs. TSGs, a set of genes that are often inactivated in cancers, can regulate various biological processes. In this study, an overview of the contribution of TSGs to PTX resistance and their underlying relationship in cancers are reported by using GeneMANIA, a web-based tool for gene/protein function prediction.

Methods

Using PubMed online database and Google web site, the terms “paclitaxel resistance” or “taxol resistance” or “drug resistance” or “chemotherapy resistance”, and “cancer” or “carcinoma”, and “tumor suppressor genes” or “TSGs” or “negative regulated protein” or “antioncogenes” were searched and analyzed. GeneMANIA data base was used to predict gene/protein interactions and functions.

Results

We identified 22 TSGs involved in PTX resistance, including BRCA1, TP53, PTEN, APC, CDKN1A, CDKN2A, HIN-1, RASSF1, YAP, ING4, PLK2, FBW7, BLU, LZTS1, REST, FADD, PDCD4, TGFBI, ING1, Bax, PinX1 and hEx. The TSGs were found to have direct and indirect relationships with each other, and thus they could contribute to PTX resistance as a group. The varied expression status and regulation function of the TSGs on cell cycle in different cancers might play an important role in PTX resistance.

Conclusion

A further understanding of the roles of tumor suppressor genes in drug resistance is an important step to overcome chemotherapy tolerance. Tumor suppressor gene therapy targets the altered genes and signaling pathways and can be a new strategy to reverse chemotherapy resistance.

Interlog protein network: an evolutionary benchmark of protein interaction networks for the evaluation of clustering algorithms

BACKGROUND

In the field of network science, exploring principal and crucial modules or communities is critical in the deduction of relationships and organization of complex networks. This approach expands an arena, and thus allows further study of biological functions in the field of network biology. As the clustering algorithms that are currently employed in finding modules have innate uncertainties, external and internal validations are necessary.

METHODS

Sequence and network structure alignment, has been used to define the Interlog Protein Network (IPN). This network is an evolutionarily conserved network with communal nodes and less false-positive links. In the current study, the IPN is employed as an evolution-based benchmark in the validation of the module finding methods. The clustering results of five algorithms; Markov Clustering (MCL), Restricted Neighborhood Search Clustering (RNSC), Cartographic Representation (CR), Laplacian Dynamics (LD) and Genetic Algorithm; to find communities in Protein-Protein Interaction networks (GAPPI) are assessed by IPN in four distinct Protein-Protein Interaction Networks (PPINs).

RESULTS

The MCL shows a more accurate algorithm based on this evolutionary benchmarking approach. Also, the biological relevance of proteins in the IPN modules generated by MCL is compatible with biological standard databases such as Gene Ontology, KEGG and Reactome.

CONCLUSION

In this study, the IPN shows its potential for validation of clustering algorithms due to its biological logic and straightforward implementation.

Deciphering complex patterns of class-I HLA-peptide cross-reactivity via hierarchical grouping

T-cell responses in humans are initiated by the binding of a peptide antigen to a human leukocyte antigen (HLA) molecule. The peptide-HLA complex then recruits an appropriate T cell, leading to cell-mediated immunity. More than 2000 HLA class-I alleles are known in humans, and they vary only in their peptide-binding grooves. The polymorphism they exhibit enables them to bind a wide range of peptide antigens from diverse sources. HLA molecules and peptides present a complex molecular recognition pattern, as many peptides bind to a given allele and a given peptide can be recognized by many alleles. A powerful grouping scheme that not only provides an insightful classification, but is also capable of dissecting the physicochemical basis of recognition specificity is necessary to address this complexity. We present a hierarchical classification of 2010 class-I alleles by using a systematic divisive clustering method. All-pair distances of alleles were obtained by comparing binding pockets in the structural models. By varying the similarity thresholds, a multilevel classification was obtained, with 7 supergroups, each further subclassifying to yield 72 groups. An independent clustering performed based only on similarities in their epitope pools correlated highly with pocket-based clustering. Physicochemical feature combinations that best explain the basis of clustering are identified. Mutual information calculated for the set of peptide ligands enables identification of binding site residues contributing to peptide specificity. The grouping of HLA molecules achieved here will be useful for rational vaccine design, understanding disease susceptibilities and predicting risk of organ transplants.

A genetically modified protein-based hydrogel for 3D culture of AD293 cells

Hydrogels have strong application prospects for drug delivery, tissue engineering and cell therapy because of their excellent biocompatibility and abundant availability as scaffolds for drugs and cells. In this study, we created hybrid hydrogels based on a genetically modified tax interactive protein-1 (TIP1) by introducing two or four cysteine residues in the primary structure of TIP1. The introduced cysteine residues were crosslinked with a four-armed poly (ethylene glycol) having their arm ends capped with maleimide residues (4-armed-PEG-Mal) to form hydrogels. In one form of the genetically modification, we incorporated a peptide sequence ‘GRGDSP’ to introduce bioactivity to the protein, and the resultant hydrogel could provide an excellent environment for a three dimensional cell culture of AD293 cells. The AD293 cells continued to divide and displayed a polyhedron or spindle-shape during the 3-day culture period. Besides, AD293 cells could be easily separated from the cell-gel constructs for future large-scale culture after being cultured for 3 days and treating hydrogel with trypsinase. This work significantly expands the toolbox of recombinant proteins for hydrogel formation, and we believe that our hydrogel will be of considerable interest to those working in cell therapy and controlled drug delivery.

Toward a systems-level understanding of gene regulatory, protein interaction, and metabolic networks in cyanobacteria

Cyanobacteria are essential primary producers in marine ecosystems, playing an important role in both carbon and nitrogen cycles. In the last decade, various genome sequencing and metagenomic projects have generated large amounts of genetic data for cyanobacteria. This wealth of data provides researchers with a new basis for the study of molecular adaptation, ecology and evolution of cyanobacteria, as well as for developing biotechnological applications. It also facilitates the use of multiplex techniques, i.e., expression profiling by high-throughput technologies such as microarrays, RNA-seq, and proteomics. However, exploration and analysis of these data is challenging, and often requires advanced computational methods. Also, they need to be integrated into our existing framework of knowledge to use them to draw reliable biological conclusions. Here, systems biology provides important tools. Especially, the construction and analysis of molecular networks has emerged as a powerful systems-level framework, with which to integrate such data, and to better understand biological relevant processes in these organisms. In this review, we provide an overview of the advances and experimental approaches undertaken using multiplex data from genomic, transcriptomic, proteomic, and metabolomic studies in cyanobacteria. Furthermore, we summarize currently available web-based tools dedicated to cyanobacteria, i.e., CyanoBase, CyanoEXpress, ProPortal, Cyanorak, CyanoBIKE, and CINPER. Finally, we present a case study for the freshwater model cyanobacteria, Synechocystis sp. PCC6803, to show the power of meta-analysis, and the potential to extrapolate acquired knowledge to the ecologically important marine cyanobacteria genus, Prochlorococcus.