Chapter 5: Network biology approach to complex diseases

Uncovering the Anti-Lung-Cancer Mechanisms of the Herbal Drug FDY2004 by Network Pharmacology

With growing evidence on the therapeutic efficacy and safety of herbal drugs, there has been a substantial increase in their application in the lung cancer treatment. Meanwhile, their action mechanisms at the system level have not been comprehensively uncovered. To this end, we employed a network pharmacology methodology to elucidate the systematic action mechanisms of FDY2004, an anticancer herbal drug composed of Moutan Radicis Cortex, Persicae Semen, and Rhei Radix et Rhizoma, in lung cancer treatment. By evaluating the pharmacokinetic properties of the chemical compounds present in FDY2004 using herbal medicine-associated databases, we identified its 29 active chemical components interacting with 141 lung cancer-associated therapeutic targets in humans. The functional enrichment analysis of the lung cancer-related targets of FDY2004 revealed the enriched Gene Ontology terms, involving the regulation of cell proliferation and growth, cell survival and death, and oxidative stress responses. Moreover, we identified key FDY2004-targeted oncogenic and tumor-suppressive pathways associated with lung cancer, including the phosphatidylinositol 3-kinase-Akt, mitogen-activated protein kinase, tumor necrosis factor, Ras, focal adhesion, and hypoxia-inducible factor-1 signaling pathways. Overall, our study provides novel evidence and basis for research on the comprehensive anticancer mechanisms of herbal medicines in lung cancer treatment.

Automated inference of Boolean models from molecular interaction maps using CaSQ

Motivation

Molecular interaction maps have emerged as a meaningful way of representing biological mechanisms in a comprehensive and systematic manner. However, their static nature provides limited insights to the emerging behaviour of the described biological system under different conditions. Computational modelling provides the means to study dynamic properties through in silico simulations and perturbations. We aim to bridge the gap between static and dynamic representations of biological systems with CaSQ, a software tool that infers Boolean rules based on the topology and semantics of molecular interaction maps built with CellDesigner.

Results

We developed CaSQ by defining conversion rules and logical formulas for inferred Boolean models according to the topology and the annotations of the starting molecular interaction maps. We used CaSQ to produce executable files of existing molecular maps that differ in size, complexity and the use of Systems Biology Graphical Notation (SBGN) standards. We also compared, where possible, the manually built logical models corresponding to a molecular map to the ones inferred by CaSQ. The tool is able to process large and complex maps built with CellDesigner (either following SBGN standards or not) and produce Boolean models in a standard output format, Systems Biology Marked Up Language-qualitative (SBML-qual), that can be further analyzed using popular modelling tools. References, annotations and layout of the CellDesigner molecular map are retained in the obtained model, facilitating interoperability and model reusability.

Availability and implementation

The present tool is available online: https://lifeware.inria.fr/∼soliman/post/casq/ and distributed as a Python package under the GNU GPLv3 license. The code can be accessed here: https://gitlab.inria.fr/soliman/casq.

Supplementary information

Supplementary data are available at Bioinformatics online.

A Network Pharmacology Study on the Molecular Mechanisms of FDY003 for Breast Cancer Treatment

Herbal medicines have drawn considerable attention with regard to their potential applications in breast cancer (BC) treatment, a frequently diagnosed malignant disease, considering their anticancer efficacy with relatively less adverse effects. However, their mechanisms of systemic action have not been understood comprehensively. Based on network pharmacology approaches, we attempted to unveil the mechanisms of FDY003, an herbal drug comprised of Lonicera japonica Thunberg, Artemisia capillaris Thunberg, and Cordyceps militaris, against BC at a systemic level. We found that FDY003 exhibited pharmacological effects on human BC cells. Subsequently, detailed data regarding the biochemical components contained in FDY003 were obtained from comprehensive herbal medicine-related databases, including TCMSP and CancerHSP. By evaluating their pharmacokinetic properties, 18 chemical compounds in FDY003 were shown to be potentially active constituents interacting with 140 BC-associated therapeutic targets to produce the pharmacological activity. Gene ontology enrichment analysis using g:Profiler indicated that the FDY003 targets were involved in the modulation of cellular processes, involving the cell proliferation, cell cycle process, and cell apoptosis. Based on a KEGG pathway enrichment analysis, we further revealed that a variety of oncogenic pathways that play key roles in the pathology of BC were significantly enriched with the therapeutic targets of FDY003; these included PI3K-Akt, MAPK, focal adhesion, FoxO, TNF, and estrogen signaling pathways. Here, we present a network-perspective of the molecular mechanisms via which herbal drugs treat BC.

Cancers in Agreement? Exploring the Cross-Talk of Cancer Metabolomic and Transcriptomic Landscapes Using Publicly Available Data

One of the major hallmarks of cancer is the derailment of a cell's metabolism. The multifaceted nature of cancer and different cancer types is transduced by both its transcriptomic and metabolomic landscapes. In this study, we re-purposed the publicly available transcriptomic and metabolomics data of eight cancer types (breast, lung, gastric, renal, liver, colorectal, prostate, and multiple myeloma) to find and investigate differences and commonalities on a pathway level among different cancer types. Topological analysis of inferred graphical Gaussian association networks showed that cancer was strongly defined in genetic networks, but not in metabolic networks. Using different statistical approaches to find significant differences between cancer and control cases, we highlighted the difficulties of high-level data-merging and in using statistical association networks. Cancer transcriptomics and metabolomics and landscapes were characterized by changed macro-molecule production, however, only major metabolic deregulations with highly impacted pathways were found in liver cancer. Cell cycle was enriched in breast, liver, and colorectal cancer, while breast and lung cancer were distinguished by highly enriched oncogene signaling pathways. A strong inflammatory response was observed in lung cancer and, to some extent, renal cancer. This study highlights the necessity of combining different omics levels to obtain a better description of cancer characteristics.

Protein-protein interaction network analysis applied to DNA copy number profiling suggests new perspectives on the aetiology of Mayer-Rokitansky-Küster-Hauser syndrome

Mayer-Rokitansky-Küster-Hauser (MRKH) syndrome is a rare disease, characterised by the aplasia of vagina and uterus in women with a 46,XX karyotype. Most cases are sporadic, but familial recurrence has also been described. Herein, we investigated an Italian cohort of 36 unrelated MRKH patients to explore the presence of pathogenic copy number variations (CNVs) by array-CGH and MLPA assays. On the whole, aberrations were found in 9/36 (25%) patients. Interestingly, one patient showed a novel heterozygous microduplication at Xp22.33, not yet described in MRKH patients, containing the PRKX gene. Moreover, a novel duplication of a specific SHOX enhancer was highlighted by MLPA. To predict the potential significance of CNVs in MRKH pathogenesis, we provided a network analysis for protein-coding genes found in the altered genomic regions. Although not all of these genes taken individually showed a clear clinical significance, their combination in a computational network highlighted that the most relevant biological connections are related to the anatomical structure development. In conclusion, the results described in the present study identified novel genetic alterations and interactions that may be likely involved in MRKH phenotype determination, so adding new insights into the complex puzzle of MRKH disease.

NetMix: A Network-Structured Mixture Model for Reduced-Bias Estimation of Altered Subnetworks

A classic problem in computational biology is the identification of altered subnetworks: subnetworks of an interaction network that contain genes/proteins that are differentially expressed, highly mutated, or otherwise aberrant compared with other genes/proteins. Numerous methods have been developed to solve this problem under various assumptions, but the statistical properties of these methods are often unknown. For example, some widely used methods are reported to output very large subnetworks that are difficult to interpret biologically. In this work, we formulate the identification of altered subnetworks as the problem of estimating the parameters of a class of probability distributions that we call the Altered Subset Distribution (ASD). We derive a connection between a popular method, jActiveModules, and the maximum likelihood estimator (MLE) of the ASD. We show that the MLE is statistically biased, explaining the large subnetworks output by jActiveModules. Based on these insights, we introduce NetMix, an algorithm that uses Gaussian mixture models to obtain less biased estimates of the parameters of the ASD. We demonstrate that NetMix outperforms existing methods in identifying altered subnetworks on both simulated and real data, including the identification of differentially expressed genes from both microarray and RNA-seq experiments and the identification of cancer driver genes in somatic mutation data.

Dynamic Regulatory Event Mining by iDREM in Large-Scale Multi-omics Datasets During Biotic and Abiotic Stress in Plants

The system-wide complexity of genome regulation encoding the organism phenotypic diversity is well understood. However, a major challenge persists about the appropriate method to describe the systematic dynamic genome regulation event utilizing enormous multi-omics datasets. Here, we describe Interactive Dynamic Regulatory Events Miner (iDREM) which reconstructs gene-regulatory networks from temporal transcriptome, proteome, and epigenome datasets during stress to envisage "master" regulators by simulating cascades of temporal transcription-regulatory and interactome events. The iDREM is a Java-based software that integrates static and time-series transcriptomics and proteomics datasets, transcription factor (TF)-target interactions, microRNA (miRNA)-target interaction, and protein-protein interactions to reconstruct temporal regulatory network and identify significant regulators in an unsupervised manner. The hidden Markov model detects specialized manipulated pathways as well as genes to recognize statistically significant regulators (TFs/miRNAs) that diverge in temporal activity. This method can be translated to any biotic or abiotic stress in plants and animals to predict the master regulators from condition-specific multi-omics datasets including host-pathogen interactions for comprehensive understanding of manipulated biological pathways.

An Investigation of the Molecular Mechanisms Underlying the Analgesic Effect of Jakyak-Gamcho Decoction: A Network Pharmacology Study

Herbal drugs have drawn substantial interest as effective analgesic agents; however, their therapeutic mechanisms remain to be fully understood. To address this question, we performed a network pharmacology study to explore the system-level mechanisms that underlie the analgesic activity of Jakyak-Gamcho decoction (JGd; Shaoyao-Gancao-Tang in Chinese and Shakuyaku-Kanzo-To in Japanese), an herbal prescription consisting of Paeonia lactiflora Pallas and Glycyrrhiza uralensis Fischer. Based on comprehensive information regarding the pharmacological and chemical properties of the herbal constituents of JGd, we identified 57 active chemical compounds and their 70 pain-associated targets. The JGd targets were determined to be involved in the regulation of diverse biological activities as follows: calcium- and cytokine-mediated signalings, calcium ion concentration and homeostasis, cellular behaviors of muscle and neuronal cells, inflammatory response, and response to chemical, cytokine, drug, and oxidative stress. The targets were further enriched in various pain-associated signalings, including the PI3K-Akt, estrogen, ErbB, neurotrophin, neuroactive ligand-receptor interaction, HIF-1, serotonergic synapse, JAK-STAT, and cAMP pathways. Thus, these data provide a systematic basis to understand the molecular mechanisms underlying the analgesic activity of herbal drugs.

Systems Pharmacology Study of the Anticervical Cancer Mechanisms of FDY003

Increasing data support that herbal medicines are beneficial in the treatment of cervical cancer; however, their mechanisms of action remain to be elucidated. In the current study, we used a systems pharmacology approach to explore the pharmacological mechanisms of FDY003, an anticancer herbal formula comprising Lonicera japonica Thunberg, Artemisia capillaris Thunberg, and Cordyceps militaris (Linn.) Link, in the treatment of cervical cancer. Through the pharmacokinetic assessment of absorption-distribution-metabolism-excretion characteristics, we found 18 active compounds that might interact with 106 cervical cancer-related targets responsible for the pharmacological effects. FDY003 targets were significantly associated with gene ontology terms related to the regulation of cellular behaviors, including cell proliferation, cell cycle processes, cell migration, cell apoptosis, cell death, and angiogenesis. The therapeutic targets of the herbal drug were further enriched in various oncogenic pathways that are implicated in the tumorigenesis and progression of cervical cancer, including the phosphatidylinositol 3-kinase, mitogen-activated protein kinase, focal adhesion, human papillomavirus infection, and tumor necrosis factor signaling pathways. Our study provides a systematic approach to explore the anticancer properties of herbal medicines against cervical cancer.

Integrated Co-functional Network Analysis on the Resistance and Virulence Features in Acinetobacter baumannii

Acinetobacter baumannii is one of the most troublesome bacterial pathogens that pose major public health threats due to its rapidly increasing drug resistance property. It is not only derived from clinic setting but also emerges from aquaculture as a fish pathogen, which could pass the resistant genes in the food chain. Understanding the mechanism of antibiotic resistance development and pathogenesis will aid our battle with the infections caused by A. baumannii. In this study, we constructed a co-functional network by integrating multiple sources of data from A. baumannii and then used the k-shell decomposition to analyze the co-functional network. We found that genes involving in basic cellular physiological function, including genes for antibiotic resistance, tended to have high k-shell values and locate in the internal layer of our network. In contrast, the non-essential genes, such as genes associated with virulence, tended to have lower k-shell values and locate in the external layer. This finding allows us to fish out the potential antibiotic resistance factors and virulence factors. In addition, we constructed an online platform ABviresDB (https://acba.shinyapps.io/ABviresDB/) for visualization of the network and features of each gene in A. baumannii. The network analysis in this study will not only aid the study on A. baumannii but also could be referenced for the research of antibiotic resistance and pathogenesis in other bacteria.

Uncertainty Relations and Fluctuation Theorems for Bayes Nets

Recent research has considered the stochastic thermodynamics of multiple interacting systems, representing the overall system as a Bayes net. I derive fluctuation theorems governing the entropy production (EP) of arbitrary sets of the systems in such a Bayes net. I also derive "conditional" fluctuation theorems, governing the distribution of EP in one set of systems conditioned on the EP of a different set of systems. I then derive thermodynamic uncertainty relations relating the EP of the overall system to the precisions of probability currents within the individual systems.

Genome-wide Network-assisted Association and Enrichment Study of Amyloid Imaging Phenotype in Alzheimer's Disease

BACKGROUND

The etiology of Alzheimer's disease remains poorly understood at the mechanistic level, and genome-wide network-based genetics have the potential to provide new insights into the disease mechanisms.

OBJECTIVE

The study aimed to explore the collective effects of multiple genetic association signals on an AV-45 PET measure, which is a well-known Alzheimer's disease biomarker, by employing a network assisted strategy.

METHODS

First, we took advantage of a dense module search algorithm to identify modules enriched by genetic association signals in a protein-protein interaction network. Next, we performed statistical evaluation to the modules identified by dense module search, including a normalization process to adjust the topological bias in the network, a replication test to ensure the modules were not found randomly , and a permutation test to evaluate unbiased associations between the modules and amyloid imaging phenotype. Finally, topological analysis, module similarity tests and functional enrichment analysis were performed for the identified modules.

RESULTS

We identified 24 consensus modules enriched by robust genetic signals in a genome-wide association analysis. The results not only validated several previously reported AD genes (APOE, APP, TOMM40, DDAH1, PARK2, ATP5C1, PVRL2, ELAVL1, ACTN1 and NRF1), but also nominated a few novel genes (ABL1, ABLIM2) that have not been studied in Alzheimer's disease but have shown associations with other neurodegenerative diseases.

CONCLUSION

The identified genes, consensus modules and enriched pathways may provide important clues to future research on the neurobiology of Alzheimer's disease and suggest potential therapeutic targets.

Genome-wide pathogenesis interpretation using a heat diffusion-based systems genetics method and implications for gene function annotation

Background

Genetics is best dedicated to interpreting pathogenesis and revealing gene functions. The past decade has witnessed unprecedented progress in genetics, particularly in genome‐wide identification of disorder variants through Genome‐Wide Association Studies (GWAS) and Phenome‐Wide Association Studies (PheWAS). However, it is still a great challenge to use GWAS/PheWAS‐derived data to elucidate pathogenesis.

Methods

In this study, we used HotNet2, a heat diffusion‐based systems genetics algorithm, to calculate the networks for disease genes obtained from GWAS and PheWAS, with an attempt to get deeper insights into disease pathogenesis at a molecular level.

Results

Through HotNet2 calculation, significant networks for 202 (for GWAS) and 167 (for PheWAS) types of diseases were identified and evaluated, respectively. The GWAS‐derived disease networks exhibit a stronger biomedical relevance than PheWAS counterparts. Therefore, the GWAS‐derived networks were used for pathogenesis interpretation by integrating the accumulated biomedical information. As a result, the pathogenesis for 64 diseases was elucidated in terms of mutation‐caused abnormal transcriptional regulation, and 47 diseases were preliminarily interpreted in terms of mutation‐caused varied protein‐protein interactions. In addition, 3,802 genes (including 46 function‐unknown genes) were assigned with new functions by disease network information, some of which were validated through mice gene knockout experiments.

Conclusions

Systems genetics algorithm HotNet2 can efficiently establish genotype‐phenotype links at the level of biological networks. Compared with original GWAS/PheWAS results, HotNet2‐calculated disease‐gene associations have stronger biomedical significance, hence provide better interpretations for the pathogenesis of genome‐wide variants, and offer new insights into gene functions as well. These results are also helpful in drug development.

Essential Oil Phytocomplex Activity, a Review with a Focus on Multivariate Analysis for a Network Pharmacology-Informed Phytogenomic Approach

Thanks to omic disciplines and a systems biology approach, the study of essential oils and phytocomplexes has been lately rolling on a faster track. While metabolomic fingerprinting can provide an effective strategy to characterize essential oil contents, network pharmacology is revealing itself as an adequate, holistic platform to study the collective effects of herbal products and their multi-component and multi-target mediated mechanisms. Multivariate analysis can be applied to analyze the effects of essential oils, possibly overcoming the reductionist limits of bioactivity-guided fractionation and purification of single components. Thanks to the fast evolution of bioinformatics and database availability, disease-target networks relevant to a growing number of phytocomplexes are being developed. With the same potential actionability of pharmacogenomic data, phytogenomics could be performed based on relevant disease-target networks to inform and personalize phytocomplex therapeutic application.

The self-organization model reveals systematic characteristics of aging

BACKGROUND

Aging is a fundamental biological process, where key bio-markers interact with each other and synergistically regulate the aging process. Thus aging dysfunction will induce many disorders. Finding aging markers and re-constructing networks based on multi-omics data (i.e. methylation, transcriptional and so on) are informative to study the aging process. However, optimizing the model to predict aging have not been performed systemically, although it is critical to identify potential molecular mechanism of aging related diseases.

METHODS

This paper aims to model the aging self-organization system using a series of supervised learning methods, and study complex molecular mechanisms of aging at system level: i.e. optimizing the aging network; summarizing interactions between aging markers; accumulating patterns of aging markers within module; finding order-parameters in the aging self-organization system.

RESULTS

In this work, the normal aging process is modeled based on multi-omics profiles across different tissues. In addition, the computational pipeline aims to model aging self-organizing systems and study the relationship between aging and related diseases (i.e. cancers), thus provide useful indicators of aging related diseases and could help to improve prediction abilities of diagnostics.

CONCLUSIONS

The aging process could be studied thoroughly by modelling the self-organization system, where key functions and the crosstalk between aging and cancers were identified.

Network analysis infers the wilt pathogen invasion associated with non-detrimental bacteria

The microbiota colonizing the root endophytic compartment and surrounding rhizosphere soils contribute to plant growth and health. However, the key members of plant soil and endophytic microbial communities involved in inhibiting or assisting pathogen invasion remain elusive. By utilizing 16S high-throughput sequencing and a molecular ecological network (MEN) approach, we systematically studied the interactions within bacterial communities in plant endophytic compartments (stem and root) and the surrounding soil (bulk and rhizosphere) during bacterial wilt invasion. The endophytic communities were found to be strongly influenced by pathogen invasion according to analysis of microbial diversity and community structure and composition. Endophytic communities of the infected plants were primarily derived from soil communities, as assessed by the SourceTracker program, but with rare migration from soil communities to endophytic communities observed in healthy plants. Soil and endophytic microbiomes from infected plants showed modular topology and greater complexity in network analysis, and a higher number of interactions than those in healthy plants. Furthermore, interactions among microbial members revealed that pathogenic Ralstonia members were positively correlated with several bacterial genera, including Delftia, Stenotrophomonas, Bacillus, Clostridium XlVa, Fontibacillus, Acidovorax, Herminiimonas, and three unclassified bacterial genera, in infected plant roots. Our findings indicated that the pathogen invasion in the rhizosphere and endophytic compartments may be highly associated with bacteria that are normally not detrimental, and sometimes even beneficial, to plants.

Network Pharmacology-Based Investigation of the System-Level Molecular Mechanisms of the Hematopoietic Activity of Samul-Tang, a Traditional Korean Herbal Formula

Hematopoiesis is a dynamic process of the continuous production of diverse blood cell types to meet the body's physiological demands and involves complex regulation of multiple cellular mechanisms in hematopoietic stem cells, including proliferation, self-renewal, differentiation, and apoptosis. Disruption of the hematopoietic system is known to cause various hematological disorders such as myelosuppression. There is growing evidence on the beneficial effects of herbal medicines on hematopoiesis; however, their mechanism of action remains unclear. In this study, we conducted a network pharmacological-based investigation of the system-level mechanisms underlying the hematopoietic activity of Samul-tang, which is an herbal formula consisting of four herbal medicines, including Angelicae Gigantis Radix, Rehmanniae Radix Preparata, Paeoniae Radix Alba, and Cnidii Rhizoma. In silico analysis of the absorption-distribution-metabolism-excretion model identified 16 active phytochemical compounds contained in Samul-tang that may target 158 genes/proteins associated with myelosuppression to exert pharmacological effects. Functional enrichment analysis suggested that the targets of Samul-tang were significantly enriched in multiple pathways closely related to the hematopoiesis and myelosuppression development, including the PI3K-Akt, MAPK, IL-17, TNF, FoxO, HIF-1, NF-kappa B, and p53 signaling pathways. Our study provides novel evidence regarding the system-level mechanisms underlying the hematopoiesis-promoting effect of herbal medicines for hematological disorder treatment.

PANEV: an R package for a pathway-based network visualization

Background

During the last decade, with the aim to solve the challenge of post-genomic and transcriptomic data mining, a plethora of tools have been developed to create, edit and analyze metabolic pathways. In particular, when a complex phenomenon is considered, the creation of a network of multiple interconnected pathways of interest could be useful to investigate the underlying biology and ultimately identify functional candidate genes affecting the trait under investigation.

Results

PANEV (PAthway NEtwork Visualizer) is an R package set for gene/pathway-based network visualization. Based on information available on KEGG, it visualizes genes within a network of multiple levels (from 1 to n) of interconnected upstream and downstream pathways. The network graph visualization helps to interpret functional profiles of a cluster of genes.

Conclusions

The suite has no species constraints and it is ready to analyze genomic or transcriptomic outcomes. Users need to supply the list of candidate genes, specify the target pathway(s) and the number of interconnected downstream and upstream pathways (levels) required for the investigation. The package is available at https://github.com/vpalombo/PANEV.

HerGePred: Heterogeneous Network Embedding Representation for Disease Gene Prediction

The discovery of disease-causing genes is a critical step towards understanding the nature of a disease and determining a possible cure for it. In recent years, many computational methods to identify disease genes have been proposed. However, making full use of disease-related (e.g., symptoms) and gene-related (e.g., gene ontology and protein-protein interactions) information to improve the performance of disease gene prediction is still an issue. Here, we develop a heterogeneous disease-gene-related network (HDGN) embedding representation framework for disease gene prediction (called HerGePred). Based on this framework, a low-dimensional vector representation (LVR) of the nodes in the HDGN can be obtained. Then, we propose two specific algorithms, namely, an LVR-based similarity prediction and a random walk with restart on a reconstructed heterogeneous disease-gene network (RW-RDGN), to predict disease genes with high performance. First, to validate the rationality of the framework, we analyze the similarity-based overlap distribution of disease pairs and design an experiment for disease-gene association recovery, the results of which revealed that the LVR of nodes performs well at preserving the local and global network structure of the HDGN. Then, we apply tenfold cross validation and external validation to compare our methods with other well-known disease gene prediction algorithms. The experimental results show that the RW-RDGN performs better than the state-of-the-art algorithm. The prediction results of disease candidate genes are essential for molecular mechanism investigation and experimental validation. The source codes of HerGePred and experimental data are available at https://github.com/yangkuoone/HerGePred.

Binding site matching in rational drug design: algorithms and applications

Interactions between proteins and small molecules are critical for biological functions. These interactions often occur in small cavities within protein structures, known as ligand-binding pockets. Understanding the physicochemical qualities of binding pockets is essential to improve not only our basic knowledge of biological systems, but also drug development procedures. In order to quantify similarities among pockets in terms of their geometries and chemical properties, either bound ligands can be compared to one another or binding sites can be matched directly. Both perspectives routinely take advantage of computational methods including various techniques to represent and compare small molecules as well as local protein structures. In this review, we survey 12 tools widely used to match pockets. These methods are divided into five categories based on the algorithm implemented to construct binding-site alignments. In addition to the comprehensive analysis of their algorithms, test sets and the performance of each method are described. We also discuss general pharmacological applications of computational pocket matching in drug repurposing, polypharmacology and side effects. Reflecting on the importance of these techniques in drug discovery, in the end, we elaborate on the development of more accurate meta-predictors, the incorporation of protein flexibility and the integration of powerful artificial intelligence technologies such as deep learning.

Biological Network Approaches and Applications in Rare Disease Studies

Network biology has the capability to integrate, represent, interpret, and model complex biological systems by collectively accommodating biological omics data, biological interactions and associations, graph theory, statistical measures, and visualizations. Biological networks have recently been shown to be very useful for studies that decipher biological mechanisms and disease etiologies and for studies that predict therapeutic responses, at both the molecular and system levels. In this review, we briefly summarize the general framework of biological network studies, including data resources, network construction methods, statistical measures, network topological properties, and visualization tools. We also introduce several recent biological network applications and methods for the studies of rare diseases.

Applications of molecular networks in biomedicine

Due to the large interdependence between the molecular components of living systems, many phenomena, including those related to pathologies, cannot be explained in terms of a single gene or a small number of genes. Molecular networks, representing different types of relationships between molecular entities, embody these large sets of interdependences in a framework that allow their mining from a systemic point of view to obtain information. These networks, often generated from high-throughput omics datasets, are used to study the complex phenomena of human pathologies from a systemic point of view. Complementing the reductionist approach of molecular biology, based on the detailed study of a small number of genes, systemic approaches to human diseases consider that these are better reflected in large and intricate networks of relationships between genes. These networks, and not the single genes, provide both better markers for diagnosing diseases and targets for treating them. Network approaches are being used to gain insight into the molecular basis of complex diseases and interpret the large datasets associated with them, such as genomic variants. Network formalism is also suitable for integrating large, heterogeneous and multilevel datasets associated with diseases from the molecular level to organismal and epidemiological scales. Many of these approaches are available to nonexpert users through standard software packages.

Network analysis of aging acceleration reveals systematic properties of 11 types of cancers

Cancers are known to be associated with accelerated aging, but to date, there has been a paucity of systematic and in‐depth studies of the correlation between aging and cancer. DNA methylation (DNAm) profiles can be used as aging markers and utilized to construct aging predictors. In this study, we downloaded 333 paired samples of DNAm, expression and mutation profiles encompassing 11 types of tissues from The Cancer Genome Atlas public access portal. The DNAm aging scores were calculated using the Support Vector Machine regression model. The DNAm aging scores of cancers revealed significant aging acceleration compared to adjacent normal tissues. Aging acceleration‐associated mutation modules and expression modules were identified in 11 types of cancers. In addition, we constructed bipartite networks of mutations and expression, and the differential expression modules related to aging‐associated mutations were selected in 11 types of cancers using the expression quantitative trait locus method. The results of enrichment analyses also identified common functions across cancers and cancer‐specific characteristics of aging acceleration. The aging acceleration interaction network across cancers suggested a core status of thyroid carcinoma and neck squamous cell carcinoma in the aging process. In summary, we have identified correlations between aging and cancers and revealed insights into the biological functions of the modules in aging and cancers.

A dual controllability analysis of influenza virus-host protein-protein interaction networks for antiviral drug target discovery

BACKGROUND

Host factors of influenza virus replication are often found in key topological positions within protein-protein interaction networks. This work explores how protein states can be manipulated through controllability analysis: the determination of the minimum manipulation needed to drive the cell system to any desired state. Here, we complete a two-part controllability analysis of two protein networks: a host network representing the healthy cell state and an influenza A virus-host network representing the infected cell state. In this context, controllability analyses aim to identify key regulating host factors of the infected cell's progression. This knowledge can be utilized in further biological analysis to understand disease dynamics and isolate proteins for study as drug target candidates.

RESULTS

Both topological and controllability analyses provide evidence of wide-reaching network effects stemming from the addition of viral-host protein interactions. Virus interacting and driver host proteins are significant both topologically and in controllability, therefore playing important roles in cell behavior during infection. Functional analysis finds overlap of results with previous siRNA studies of host factors involved in influenza replication, NF-kB pathway and infection relevance, and roles as interferon regulating genes. 24 proteins are identified as holding regulatory roles specific to the infected cell by measures of topology, controllability, and functional role. These proteins are recommended for further study as potential antiviral drug targets.

CONCLUSIONS

Seasonal outbreaks of influenza A virus are a major cause of illness and death around the world each year with a constant threat of pandemic infection. This research aims to increase the efficiency of antiviral drug target discovery using existing protein-protein interaction data and network analysis methods. These results are beneficial to future studies of influenza virus, both experimental and computational, and provide evidence that the combination of topology and controllability analyses may be valuable for future efforts in drug target discovery.

Identification of key regulatory genes connected to NF-κB family of proteins in visceral adipose tissues using gene expression and weighted protein interaction network

Obesity is connected to the activation of chronic inflammatory pathways in both adipocytes and macrophages located in adipose tissues. The nuclear factor (NF)-κB is a central molecule involved in inflammatory pathways linked to the pathology of different complex metabolic disorders. Investigating the gene expression data in the adipose tissue would potentially unravel disease relevant gene interactions. The present study is aimed at creating a signature molecular network and at prioritizing the potential biomarkers interacting with NF-κB family of proteins in obesity using system biology approaches. The dataset GSE88837 associated with obesity was downloaded from Gene Expression Omnibus (GEO) database. Statistical analysis represented the differential expression of a total of 2650 genes in adipose tissues (p = <0.05). Using concepts like correlation, semantic similarity, and theoretical graph parameters we narrowed down genes to a network of 23 genes strongly connected with NF-κB family with higher significance. Functional enrichment analysis revealed 21 of 23 target genes of NF-κB were found to have a critical role in the pathophysiology of obesity. Interestingly, GEM and PPP1R13L were predicted as novel genes which may act as potential target or biomarkers of obesity as they occur with other 21 target genes with known obesity relationship. Our study concludes that NF-κB and prioritized target genes regulate the inflammation in adipose tissues through several molecular signaling pathways like NF-κB, PI3K-Akt, glucocorticoid receptor regulatory network, angiogenesis and cytokine pathways. This integrated system biology approaches can be applied for elucidating functional protein interaction networks of NF-κB protein family in different complex diseases. Our integrative and network-based approach for finding therapeutic targets in genomic data could accelerate the identification of novel drug targets for obesity.

Gene expression changes associated with the evolutionary loss of a metabolic trait: lack of lipogenesis in parasitoids

Background

Trait loss is a pervasive phenomenon in evolution, yet the underlying molecular causes have been identified in only a handful of cases. Most of these cases involve loss-of-function mutations in one or more trait-specific genes. However, in parasitoid insects the evolutionary loss of a metabolic trait is not associated with gene decay. Parasitoids have lost the ability to convert dietary sugars into fatty acids. Earlier research suggests that lack of lipogenesis in the parasitoid wasp Nasonia vitripennis is caused by changes in gene regulation.

Results

We compared transcriptomic responses to sugar-feeding in the non-lipogenic parasitoid species Nasonia vitripennis and the lipogenic Drosophila melanogaster. Both species adjusted their metabolism within 4 hours after sugar-feeding, but there were sharp differences between the expression profiles of the two species, especially in the carbohydrate and lipid metabolic pathways. Several genes coding for key enzymes in acetyl-CoA metabolism, such as malonyl-CoA decarboxylase (mcd) and HMG-CoA synthase (hmgs) differed in expression between the two species. Their combined action likely blocks lipogenesis in the parasitoid species. Network-based analysis showed connectivity of genes to be negatively correlated to the fold change of gene expression. Furthermore, genes involved in the fatty acid metabolic pathway were more connected than the set of genes of all metabolic pathways combined.

Conclusion

High connectivity of lipogenesis genes is indicative of pleiotropic effects and could explain the absence of gene degradation. We conclude that modification of expression levels of only a few little-connected genes, such as mcd, is sufficient to enable complete loss of lipogenesis in N. vitripennis.

Electronic supplementary material

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

The POLR2E rs3787016 polymorphism is strongly associated with the risk of female breast and cervical cancer

The rs3787016 polymorphism, in polymerase II polypeptide E (POLR2E), was previously identified as being associated with the risk for prostate cancer, esophageal cancer, breast cancer, papillary thyroid carcinoma and liver cancer, suggesting that rs3787016 may server as a common genetic factor to affect individual susceptibility to cancer. To prove the hypothesis, we here performed a case-control study to explore the association between rs3787016 and cervical cancer risk, and to confirm the association between rs3787016 and breast cancer in a central Chinese population, which was followed by a meta-analysis to precisely estimate the association between rs3787016 and risk of female breast and cervical cancer. The genotyping was performed by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) and confirmed by sequencing. Our results indicated that rs3787016 was associated with the risk of both breast cancer and cervical cancer, and stratified analysis indicated that the association remained particularly for ≤60 years old females who smoke and drink. Moreover, after grouping breast cancer and cervical cancer together, our meta-analysis demonstrated that rs3787016 was associated with overall cancer risk and breast cancer risk. Collectively, the POLR2E rs3787016 polymorphism may be a valuable biomarker for female breast and cervical cancer predisposition.

Integrative Analysis of Dysfunctional Modules Driven by Genomic Alterations at System Level Across 11 Cancer Types

Aim and Objective:

Integrating multi-omics data to identify driver genes and key biological functions for tumorigenesis remains a major challenge.

Method:

A new computational pipeline was developed to identify the Driver Mutation-Differential Co-Expression (DM-DCE) modules based on dysfunctional networks across 11 TCGA cancers.

Results:

Functional analyses provided insight into the properties of various cancers, and found common cellular signals / pathways of cancers. Furthermore, the corresponding network analysis identified conservations or interactions across different types of cancers, thus the crosstalk between the key signaling pathways, immunity and cancers was found. Clinical analysis also identified key prognostic / survival patterns.

Conclusion:

Taken together, our study sheds light on both cancer-specific and cross-cancer characteristics systematically.

Construction of a Suite of Computable Biological Network Models Focused on Mucociliary Clearance in the Respiratory Tract

Mucociliary clearance (MCC), considered as a collaboration of mucus secreted from goblet cells, the airway surface liquid layer, and the beating of cilia of ciliated cells, is the airways’ defense system against airborne contaminants. Because the process is well described at the molecular level, we gathered the available information into a suite of comprehensive causal biological network (CBN) models. The suite consists of three independent models that represent (1) cilium assembly, (2) ciliary beating, and (3) goblet cell hyperplasia/metaplasia and that were built in the Biological Expression Language, which is both human-readable and computable. The network analysis of highly connected nodes and pathways demonstrated that the relevant biology was captured in the MCC models. We also show the scoring of transcriptomic data onto these network models and demonstrate that the models capture the perturbation in each dataset accurately. This work is a continuation of our approach to use computational biological network models and mathematical algorithms that allow for the interpretation of high-throughput molecular datasets in the context of known biology. The MCC network model suite can be a valuable tool in personalized medicine to further understand heterogeneity and individual drug responses in complex respiratory diseases.

Network spectra for drug-target identification in complex diseases: new guns against old foes

The fundamental understanding of altered complex molecular interactions in a diseased condition is the key to its cure. The overall functioning of these molecules is kind of jugglers play in the cell orchestra and to anticipate these relationships among the molecules is one of the greatest challenges in modern biology and medicine. Network science turned out to be providing a successful and simple platform to understand complex interactions among healthy and diseased tissues. Furthermore, much information about the structure and dynamics of a network is concealed in the eigenvalues of its adjacency matrix. In this review, we illustrate rapid advancements in the field of network science in combination with spectral graph theory that enables us to uncover the complexities of various diseases. Interpretations laid by network science approach have solicited insights into molecular relationships and have reported novel drug targets and biomarkers in various complex diseases.

Personalized Medicine and Molecular Interaction Networks in Amyotrophic Lateral Sclerosis (ALS): Current Knowledge

Multiple genes and mechanisms of pathophysiology have been implicated in amyotrophic lateral sclerosis (ALS), suggesting it is a complex systemic disease. With this in mind, applying personalized medicine (PM) approaches to tailor treatment pipelines for ALS patients may be necessary. The modelling and analysis of molecular interaction networks could represent valuable resources in defining ALS-associated pathways and discovering novel therapeutic targets. Here we review existing omics datasets and analytical approaches, in order to consider how molecular interaction networks could improve our understanding of the molecular pathophysiology of this fatal neuromuscular disorder.

Integrating molecular networks with genetic variant interpretation for precision medicine

More reliable and cheaper sequencing technologies have revealed the vast mutational landscapes characteristic of many phenotypes. The analysis of such genetic variants has led to successful identification of altered proteins underlying many Mendelian disorders. Nevertheless the simple one‐variant one‐phenotype model valid for many monogenic diseases does not capture the complexity of polygenic traits and disorders. Although experimental and computational approaches have improved detection of functionally deleterious variants and important interactions between gene products, the development of comprehensive models relating genotype and phenotypes remains a challenge in the field of genomic medicine. In this context, a new view of the pathologic state as significant perturbation of the network of interactions between biomolecules is crucial for the identification of biochemical pathways associated with complex phenotypes. Seminal studies in systems biology combined the analysis of genetic variation with protein–protein interaction networks to demonstrate that even as biological systems evolve to be robust to genetic variation, their topologies create disease vulnerabilities. More recent analyses model the impact of genetic variants as changes to the “wiring” of the interactome to better capture heterogeneity in genotype–phenotype relationships. These studies lay the foundation for using networks to predict variant effects at scale using machine‐learning or algorithmic approaches. A wealth of databases and resources for the annotation of genotype–phenotype relationships have been developed to support developments in this area. This overview describes how study of the molecular interactome has generated insights linking the organization of biological systems to disease mechanism, and how this information can enable precision medicine.

This article is categorized under:

Translational, Genomic, and Systems Medicine > Translational Medicine

Biological Mechanisms > Cell Signaling

Models of Systems Properties and Processes > Mechanistic Models

Analytical and Computational Methods > Computational Methods

SyNDI: synchronous network data integration framework

BACKGROUND

Systems biology takes a holistic approach by handling biomolecules and their interactions as big systems. Network based approach has emerged as a natural way to model these systems with the idea of representing biomolecules as nodes and their interactions as edges. Very often the input data come from various sorts of omics analyses. Those resulting networks sometimes describe a wide range of aspects, for example different experiment conditions, species, tissue types, stimulating factors, mutants, or simply distinct interaction features of the same network produced by different algorithms. For these scenarios, synchronous visualization of more than one distinct network is an excellent mean to explore all the relevant networks efficiently. In addition, complementary analysis methods are needed and they should work in a workflow manner in order to gain maximal biological insights.

RESULTS

In order to address the aforementioned needs, we have developed a Synchronous Network Data Integration (SyNDI) framework. This framework contains SyncVis, a Cytoscape application for user-friendly synchronous and simultaneous visualization of multiple biological networks, and it is seamlessly integrated with other bioinformatics tools via the Galaxy platform. We demonstrated the functionality and usability of the framework with three biological examples - we analyzed the distinct connectivity of plasma metabolites in networks associated with high or low latent cardiovascular disease risk; deeper insights were obtained from a few similar inflammatory response pathways in Staphylococcus aureus infection common to human and mouse; and regulatory motifs which have not been reported associated with transcriptional adaptations of Mycobacterium tuberculosis were identified.

CONCLUSIONS

Our SyNDI framework couples synchronous network visualization seamlessly with additional bioinformatics tools. The user can easily tailor the framework for his/her needs by adding new tools and datasets to the Galaxy platform.

Integrated network analysis and machine learning approach for the identification of key genes of triple-negative breast cancer

Triple-negative breast cancer (TNBC) has attracted more attention compared with other breast cancer subtypes due to its aggressive nature, poor prognosis, and chemotherapy remains the mainstay of treatment with no other approved targeted therapy. Therefore, the study aimed to discover more promising therapeutic targets and investigating new insights of biological mechanism of TNBC. Six microarray data sets consisting of 463 non-TNBC and 405 TNBC samples were mined from Gene Expression Omnibus. The data sets were integrated by meta-analysis and identified 1075 differentially expressed genes. Protein-protein interaction network was constructed which consists of 486 nodes and 1932 edges, where 29 hub genes were obtained with high topological measures. Further, 16 features (hub genes), 12 upregulated (AURKB, CCNB2, CDC20, DDX18, EGFR, ENO1, MYC, NUP88, PLK1, PML, POLR2F, and SKP2) and four downregulated ( CCND1, GLI3, SKP1, and TGFB3) were selected through machine learning correlation based feature selection method on training data set. A naïve Bayes based classifier built using the expression profiles of 16 features (hub genes) accurately and reliably classify TNBC from non-TNBC samples in the validation test data set with a receiver operating curve of 0.93 to 0.98. Subsequently, Gene Ontology analysis revealed that the hub genes were enriched in mitotic cell cycle processes and Kyoto Encyclopedia of Genes and Genomes pathway analysis showed that they were enriched in cell cycle pathways. Thus, the identified key hub genes and pathways highlighted in the study would enhance the understanding of molecular mechanism of TNBC which may serve as potential therapeutic target.

Integration of Molecular Interactome and Targeted Interaction Analysis to Identify a COPD Disease Network Module

The polygenic nature of complex diseases offers potential opportunities to utilize network-based approaches that leverage the comprehensive set of protein-protein interactions (the human interactome) to identify new genes of interest and relevant biological pathways. However, the incompleteness of the current human interactome prevents it from reaching its full potential to extract network-based knowledge from gene discovery efforts, such as genome-wide association studies, for complex diseases like chronic obstructive pulmonary disease (COPD). Here, we provide a framework that integrates the existing human interactome information with experimental protein-protein interaction data for FAM13A, one of the most highly associated genetic loci to COPD, to find a more comprehensive disease network module. We identified an initial disease network neighborhood by applying a random-walk method. Next, we developed a network-based closeness approach (C_AB) that revealed 9 out of 96 FAM13A interacting partners identified by affinity purification assays were significantly close to the initial network neighborhood. Moreover, compared to a similar method (local radiality), the C_AB approach predicts low-degree genes as potential candidates. The candidates identified by the network-based closeness approach were combined with the initial network neighborhood to build a comprehensive disease network module (163 genes) that was enriched with genes differentially expressed between controls and COPD subjects in alveolar macrophages, lung tissue, sputum, blood, and bronchial brushing datasets. Overall, we demonstrate an approach to find disease-related network components using new laboratory data to overcome incompleteness of the current interactome.

Co-Expression Network Analysis Identifies miRNA⁻mRNA Networks Potentially Regulating Milk Traits and Blood Metabolites

MicroRNAs (miRNA) regulate mRNA networks to coordinate cellular functions. In this study, we constructed gene co-expression networks to detect miRNA modules (clusters of miRNAs with similar expression patterns) and miRNA–mRNA pairs associated with blood (triacylglyceride and nonesterified fatty acids) and milk (milk yield, fat, protein, and lactose) components and milk fatty acid traits following dietary supplementation of cows’ diets with 5% linseed oil (LSO) (n = 6 cows) or 5% safflower oil (SFO) (n = 6 cows) for 28 days. Using miRNA transcriptome data from mammary tissues of cows for co-expression network analysis, we identified three consensus modules: blue, brown, and turquoise, composed of 70, 34, and 86 miRNA members, respectively. The hub miRNAs (miRNAs with the most connections with other miRNAs) were miR-30d, miR-484 and miR-16b for blue, brown, and turquoise modules, respectively. Cell cycle arrest, and p53 signaling and transforming growth factor–beta (TGF-β) signaling pathways were the common gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways enriched for target genes of the three modules. Protein percent (p = 0.03) correlated with the turquoise module in LSO treatment while protein yield (p = 0.003) and milk yield (p = 7 × 10⁻⁰⁴) correlated with the turquoise model, protein and milk yields and lactose percent (p < 0.05) correlated with the blue module and fat percent (p = 0.04) correlated with the brown module in SFO treatment. Several fatty acids correlated (p < 0.05) with the blue (CLA:9,11) and brown (C4:0, C12:0, C22:0, C18:1n9c and CLA:10,12) modules in LSO treatment and with the turquoise (C14:0, C18:3n3 and CLA:9,11), blue (C14:0 and C23:0) and brown (C6:0, C16:0, C22:0, C22:6n3 and CLA:10,12) modules in SFO treatment. Correlation of miRNA and mRNA data from the same animals identified the following miRNA–mRNA pairs: miR-183/RHBDD2 (p = 0.003), miR-484/EIF1AD (p = 0.011) and miR-130a/SBSPON (p = 0.004) with lowest p-values for the blue, brown, and turquoise modules, respectively. Milk yield, protein yield, and protein percentage correlated (p < 0.05) with 28, 31 and 5 miRNA–mRNA pairs, respectively. Our results suggest that, the blue, brown, and turquoise modules miRNAs, hub miRNAs, miRNA–mRNA networks, cell cycle arrest GO term, p53 signaling and TGF-β signaling pathways have considerable influence on milk and blood phenotypes following dietary supplementation of dairy cows’ diets with 5% LSO or 5% SFO.

NetControl4BioMed: a pipeline for biomedical data acquisition and analysis of network controllability

BACKGROUND

Network controllability focuses on discovering combinations of external interventions that can drive a biological system to a desired configuration. In practice, this approach translates into finding a combined multi-drug therapy in order to induce a desired response from a cell; this can lead to developments of novel therapeutic approaches for systemic diseases like cancer.

RESULT

We develop a novel bioinformatics data analysis pipeline called NetControl4BioMed based on the concept of target structural control of linear networks. Our pipeline generates novel molecular interaction networks by combining pathway data from various public databases starting from the user's query. The pipeline then identifies a set of nodes that is enough to control a given, user-defined set of disease-specific essential proteins in the network, i.e., it is able to induce a change in their configuration from any initial state to any final state. We provide both the source code of the pipeline as well as an online web-service based on this pipeline http://combio.abo.fi/nc/net_control/remote_call.php .

CONCLUSION

The pipeline can be used by researchers for controlling and better understanding of molecular interaction networks through combinatorial multi-drug therapies, for more efficient therapeutic approaches and personalised medicine.

Supplementary data for the biological age linked to oxidative stress modifies breast cancer aggressiveness

The data presented in this article are related to the research paper entitled “The biological age linked to oxidative stress modifies breast cancer aggressiveness” (M.M. Sáez-Freire, A. Blanco-Gómez, S. Castillo-Lluva, A. Gómez-Vecino, J.M. Galvis-Jiménez, C. Martín-Seisdedos, M. Isidoro-García, L. Hontecillas-Prieto, M.B. García-Cenador, F.J. García-Criado, M.C. Patino-Alonso, P. Galindo-Villardón, J.H. Mao, C. Prieto, A. Castellanos-Martín, L. Kaderali, J. Pérez-Losada). The data shown were obtained from a population of transgenic mice, MMTV-Erbb2/Neu, with different susceptibility to breast cancer and a mixed genetic background generated by backcrossing. It was observed that the aggressiveness of breast cancer negatively correlates with age, being lower in chronologically old mice, similar to what occurs in humans. Given that oxidative stress is associated with tumour susceptibility and the degree of aging, the association between the aggressiveness of breast cancer and multiple intermediate phenotypes directly or indirectly related to oxidative stress was studied. Using a mathematical model, we defined biological age and the degree of aging as the difference between biological and chronological ages. As a result, we observed that biologically old mice predominated among those that developed the disease early on, that is, those that were chronologically young. We then identified the specific and common genetic components of Quantitative Trait loci or QTL associated with different evolution of breast cancer, the intermediate phenotypes related to oxidative stress studied, the biological age and the degree of aging. Lastly, we showed that the expression pattern in the livers of biologically old mice were enriched in signalling pathways related to inflammation and response to infections; whereas the biologically young mice exhibited enriched pathways related to mitochondrial activity. For the explanation and discussion of these data refer to the research article cited above.

Differential network as an indicator of osteoporosis with network entropy

Osteoporosis is a common skeletal disorder characterized by a decrease in bone mass and density. The peak bone mass (PBM) is a significant determinant of osteoporosis. To gain insights into the indicating effect of PBM to osteoporosis, this study focused on characterizing the PBM networks and identifying key genes. One biological data set with 12 monocyte low PBM samples and 11 high PBM samples was derived to construct protein-protein interaction networks (PPINs). Based on clique-merging, module-identification algorithm was used to identify modules from PPINs. The systematic calculation and comparison were performed to test whether the network entropy can discriminate the low PBM network from high PBM network. We constructed 32 destination networks with 66 modules divided from monocyte low and high PBM networks. Among them, network 11 was the only significantly differential one (P<0.05) with 8 nodes and 28 edges. All genes belonged to precursors of osteoclasts, which were related to calcium transport as well as blood monocytes. In conclusion, based on the entropy in PBM PPINs, the differential network appears to be a novel therapeutic indicator for osteoporosis during the bone monocyte progression; these findings are helpful in disclosing the pathogenetic mechanisms of osteoporosis.

Network biology discovers pathogen contact points in host protein-protein interactomes

In all organisms, major biological processes are controlled by complex protein-protein interactions networks (interactomes), yet their structural complexity presents major analytical challenges. Here, we integrate a compendium of over 4300 phenotypes with Arabidopsis interactome (AI-1MAIN). We show that nodes with high connectivity and betweenness are enriched and depleted in conditional and essential phenotypes, respectively. Such nodes are located in the innermost layers of AI-1MAIN and are preferential targets of pathogen effectors. We extend these network-centric analyses to Cell Surface Interactome (CSILRR) and predict its 35 most influential nodes. To determine their biological relevance, we show that these proteins physically interact with pathogen effectors and modulate plant immunity. Overall, our findings contrast with centrality-lethality rule, discover fast information spreading nodes, and highlight the structural properties of pathogen targets in two different interactomes. Finally, this theoretical framework could possibly be applicable to other inter-species interactomes to reveal pathogen contact points.

The biological age linked to oxidative stress modifies breast cancer aggressiveness

The incidence of breast cancer increases with age until menopause, and breast cancer is more aggressive in younger women. The existence of epidemiological links between breast cancer and aging indicates that both processes share some common mechanisms of development. Oxidative stress is associated with both cancer susceptibility and aging. Here we observed that ERBB2-positive breast cancer, which developed in genetically heterogeneous ERBB2-positive transgenic mice generated by a backcross, is more aggressive in chronologically younger than in older mice (differentiated by the median survival of the cohort that was 79 weeks), similar to what occurs in humans. In this cohort, we estimated the oxidative biological age using a mathematical model that integrated several subphenotypes directly or indirectly related to oxidative stress. The model selected the serum levels of HDL-cholesterol and magnesium and total AKT1 and glutathione concentrations in the liver. The grade of aging was calculated as the difference between the predicted biological age and the chronological age. This comparison permitted the identification of biologically younger and older mice compared with their chronological age. Interestingly, biologically older mice developed more aggressive breast cancer than the biologically younger mice. Genomic regions on chromosomes 2 and 15 linked to the grade of oxidative aging were identified. The levels of expression of Zbp1 located on chromosome 2, a gene related to necroptosis and inflammation, positively correlated with the grade of aging and tumour aggressiveness. Moreover, the pattern of gene expression of genes linked to the inflammation and the response to infection pathways was enriched in the livers of biologically old mice. This study shows part of the complex interactions between breast cancer and aging.

Applications of Bayesian network models in predicting types of hematological malignancies

Network analysis is the preferred approach for the detection of subtle but coordinated changes in expression of an interacting and related set of genes. We introduce a novel method based on the analyses of coexpression networks and Bayesian networks, and we use this new method to classify two types of hematological malignancies; namely, acute myeloid leukemia (AML) and myelodysplastic syndrome (MDS). Our classifier has an accuracy of 93%, a precision of 98%, and a recall of 90% on the training dataset (n = 366); which outperforms the results reported by other scholars on the same dataset. Although our training dataset consists of microarray data, our model has a remarkable performance on the RNA-Seq test dataset (n = 74, accuracy = 89%, precision = 88%, recall = 98%), which confirms that eigengenes are robust with respect to expression profiling technology. These signatures are useful in classification and correctly predicting the diagnosis. They might also provide valuable information about the underlying biology of diseases. Our network analysis approach is generalizable and can be useful for classifying other diseases based on gene expression profiles. Our previously published Pigengene package is publicly available through Bioconductor, which can be used to conveniently fit a Bayesian network to gene expression data.

Advances in metabolome information retrieval: turning chemistry into biology. Part II: biological information recovery

This work reports the second part of a review intending to give the state of the art of major metabolic phenotyping strategies. It particularly deals with inherent advantages and limits regarding data analysis issues and biological information retrieval tools along with translational challenges. This Part starts with introducing the main data preprocessing strategies of the different metabolomics data. Then, it describes the main data analysis techniques including univariate and multivariate aspects. It also addresses the challenges related to metabolite annotation and characterization. Finally, functional analysis including pathway and network strategies are discussed. The last section of this review is devoted to practical considerations and current challenges and pathways to bring metabolomics into clinical environments.

Decipher reliable biomarkers of brain aging by integrating literature-based evidence with interactome data

Aging is an inevitable progressive decline in every physiological function and serves as a primary risk factor for cognitive decline and Alzheimer's disease. Thus, age-dependent impairments in cognitive function must be understood in association with general aging processes with an integrative approach in a systemic manner. An integrative aging gene network was constructed based on mutual molecular interactions using literature-curated interactome data and separated into functionally distinct modules. To investigate key surrogate biomarkers of the aging brain in the context of the general aging process, co-expression networks were built on post-mortem and Alzheimer's brain transcriptome data. In both the normal aging brain and the brain affected by Alzheimer's disease, the immune-related co-expression module was positively correlated with advancing age, whereas the synaptic transmission-related co-expression module was decreased with age. Importantly, the network topology-based analysis indicated that complement system genes were prioritized as a surrogate biomarker in evaluating the process of brain aging. Our public data-centered analysis coupled with experimental validation revealed that the complement system is likely to be a master regulator in initiating and regulating the immune system in the aging brain and could serve as reliable and surrogate biomarkers for the diagnosis of cognitive dysfunction.

Systemic tracking of diagnostic function modules for post-menopausal osteoporosis in a differential co-expression network view

Post-menopausal osteoporosis is one of the most common bone diseases in women. The aim of the present study was to predict the diagnostic function modules from a differential co-expression gene network in order to enhance the current understanding of the biological processes and to promote the early prevention and intervention of post-menopausal osteoporosis. The diagnostic function modules were extracted from a differential co-expression network by the established protein-protein interaction (PPI) network analysis. First, significant genes were identified from the differential co-expression network, which were regarded as seed genes. Starting from the seed genes, the sub-networks in this disease, referred to as diagnostic function modules, were exhaustively searched and prioritized through a snowball sampling strategy to identify genes to accurately predict clinical outcomes. In addition, crucial function inference was performed for each diagnostic function module. Based on the microarray and PPI data, the differential co-expression network was constructed, which contained 1,607 genes and 4,197 interactions. A total of 110 seed genes were identified, and nine diagnostic modules that accurately distinguished post-menopausal osteoporosis from healthy controls were screened out from these seed genes. The diagnostic modules may be associated with five functional pathways with emphasis on metabolism. A total of nine diagnostic functional modules screened in the present study may be considered as potential targets for predicting the clinical outcomes of post-menopausal osteoporosis, and may contribute to the early diagnosis and therapy of osteoporosis.