Genome scale modeling in systems biology: algorithms and resources

PathIN: an integrated tool for the visualization of pathway interaction networks

PathIN is a web-service that provides an easy and flexible way for rapidly creating pathway-based networks at several functional biological levels: genes, compounds and reactions. The tool is supported by a database repository of reference pathway networks across a large set of species, developed through the freely available information included in the KEGG, Reactome and Wiki Pathways database repositories. PathIN provides networks by means of five diverse methodologies: (a) direct connections between pathways of interest, (b) direct connections as well as the first neighbours of the given pathways, (c) direct connections, the first neighbours and the connections in between them, and (d) two additional methodologies for creating complementary pathway-to-pathway networks that involve additional (missing) pathways that interfere in-between pathways of interest. PathIN is expected to be used as a simple yet informative reference tool for understanding networks of molecular mechanisms related to specific diseases.

PathExNET: A tool for extracting pathway expression networks from gene expression statistics

A fundamental issue related to the understanding of the molecular mechanisms, is the way in which common pathways act across different biological experiments related to complex diseases. Using network-based approaches, this work aims to provide a numeric characterization of pathways across different biological experiments, in the prospect to create unique footprints that may characterise a specific disease under study at a pathway network level. In this line we propose PathExNET, a web service that allows the creation of pathway-to-pathway expression networks that hold the over- and under expression information obtained from differential gene expression analyses. The unique numeric characterization of pathway expression status related to a specific biological experiment (or disease), as well as the creation of diverse combination of pathway networks generated by PathExNET, is expected to provide a concrete contribution towards the individualization of disease, and further lead to a more precise personalised medicine and management of treatment. PathExNET is available at: https://bioinformatics.cing.ac.cy/PathExNET and at https://pathexnet.cing-big.hpcf.cyi.ac.cy/.

A Systems Bioinformatics Approach to Interconnect Biological Pathways

Signal transduction tasks as well as other complex biological processes involve many different changes in groups of genes, proteins, and metabolites linked together in chains or networks called pathways or networks of pathways. In a classical functional analysis, the biomolecules found to play a role in the biological status under investigation are members of a group of pathways that are not necessarily interconnected. However, interconnectivity is a critical factor for functionality. Thus, it is necessary to be able to construct "connected functional stories" to understand better the complex biological processes. PathwayConnector is a recently introduced web-tool that facilitates the construction of complementary pathway-to-pathway networks, bringing to our attention missing pathways that are crucial links towards the understanding of the molecular mechanisms related to complex diseases. Current version of the web-tool draws from an expanded pathway reference network and provides information deriving from 19 different organisms and 2 different pathway repositories: the KEGG and the REACTOME. Novel genes, proteins, and pathways derived from any experimental/computational method either in large-scale (omics) or even in smaller scale (specific laboratory experiments) can potentially be projected and analyzed through PathwayConnector. This chapter describes in details the pipeline and methodologies used for the latest updated version of PathwayConnector, providing an easy way for rapidly relating human or other organism's pathways together. Recent studies have shown that pathway networks and subnetworks, generated by PathwayConnector, are an integral part towards the individualization of disease, leading to a more precise and personalized management of the treatment.

Integrating Single-Step GWAS and Bipartite Networks Reconstruction Provides Novel Insights into Yearling Weight and Carcass Traits in Hanwoo Beef Cattle

Simple Summary

Hanwoo is an indigenous cattle breed in Korea and popular for meat production owing to its rapid growth and high-quality meat. Its yearling weight and carcass traits (backfat thickness, carcass weight, eye muscle area, and marbling score) are economically important for the selection of young and proven bulls. In recent decades, the advent of high throughput genotyping technologies has made it possible to perform genome-wide association studies (GWAS) for the detection of genomic regions associated with traits of economic interest in different species. In this study, we conducted a weighted single-step genome-wide association study which combines all genotypes, phenotypes and pedigree data in one step (ssGBLUP). It allows for the use of all SNPs simultaneously along with all phenotypes from genotyped and ungenotyped animals. Our results revealed 33 relevant genomic regions related to the traits of interest. Gene set enrichment analysis indicated that the identified candidate genes were related to biological processes and functional terms that were involved in growth and lipid metabolism. In conclusion, these results suggest that the incorporation of GWAS results and network analysis can help us to better understand the genetic bases underlying growth and carcass traits.

Abstract

In recent years, studies on the biological mechanisms underlying complex traits have been facilitated by innovations in high-throughput genotyping technology. We conducted a weighted single-step genome-wide association study (WssGWAS) to evaluate backfat thickness, carcass weight, eye muscle area, marbling score, and yearling weight in a cohort of 1540 Hanwoo beef cattle using BovineSNP50 BeadChip. The WssGWAS uncovered thirty-three genomic regions that explained more than 1% of the additive genetic variance, mostly located on chromosomes 6 and 14. Among the identified window regions, seven quantitative trait loci (QTL) had pleiotropic effects and twenty-six QTL were trait-specific. Significant pathways implicated in the measured traits through Gene Ontology (GO) term enrichment analysis included the following: lipid biosynthetic process, regulation of lipid metabolic process, transport or localization of lipid, regulation of growth, developmental growth, and multicellular organism growth. Integration of GWAS results of the studied traits with pathway and network analyses facilitated the exploration of the respective candidate genes involved in several biological functions, particularly lipid and growth metabolism. This study provides novel insight into the genetic bases underlying complex traits and could be useful in developing breeding schemes aimed at improving growth and carcass traits in Hanwoo beef cattle.

Advances in Astrocyte Computational Models: From Metabolic Reconstructions to Multi-omic Approaches

The growing importance of astrocytes in the field of neuroscience has led to a greater number of computational models devoted to the study of astrocytic functions and their metabolic interactions with neurons. The modeling of these interactions demands a combined understanding of brain physiology and the development of computational frameworks based on genomic-scale reconstructions, system biology, and dynamic models. These computational approaches have helped to highlight the neuroprotective mechanisms triggered by astrocytes and other glial cells, both under normal conditions and during neurodegenerative processes. In the present review, we evaluate some of the most relevant models of astrocyte metabolism, including genome-scale reconstructions and astrocyte-neuron interactions developed in the last few years. Additionally, we discuss novel strategies from the multi-omics perspective and computational models of other glial cell types that will increase our knowledge in brain metabolism and its association with neurodegenerative diseases.

ProtExA: A tool for post-processing proteomics data providing differential expression metrics, co-expression networks and functional analytics

ProTExA is a web-tool that provides a post-processing workflow for the analysis of protein and gene expression datasets. Using network-based bioinformatics approaches, ProTExA facilitates differential expression analysis and co-expression network analysis as well as pathway and post-pathway analysis. Specifically, for a given set of protein-gene expression data across samples, ProTExA: (1) performs statistical analysis and filtering to highlight the differentially expressed proteins-genes, (2) performs enrichment analysis to identify top-scored pathways, (3) generates pathway-to-pathway and pathway-to-gene networks (4) generates protein and gene co-expression networks using a variety of methodologies, and (5) applies clustering methodologies to identify sub-networks of co-expressed proteins-genes. The proposed web-tool is a simple yet informative tool, towards understanding and exploitation of protein and gene expression datasets, especially for those that do not have the expertise and local resources to replicate specific analyses in the context of collaborative and scientific data exchanging.

Multiple Pathways Involved in Palmitic Acid-Induced Toxicity: A System Biology Approach

Inflammation is a complex biological response to injuries, metabolic disorders or infections. In the brain, astrocytes play an important role in the inflammatory processes during neurodegenerative diseases. Recent studies have shown that the increase of free saturated fatty acids such as palmitic acid produces a metabolic inflammatory response in astrocytes generally associated with damaging mechanisms such as oxidative stress, endoplasmic reticulum stress, and autophagic defects. In this aspect, the synthetic neurosteroid tibolone has shown to exert protective functions against inflammation in neuronal experimental models without the tumorigenic effects exerted by sexual hormones such as estradiol and progesterone. However, there is little information regarding the specific mechanisms of tibolone in astrocytes during inflammatory insults. In the present study, we performed a genome-scale metabolic reconstruction of astrocytes that was used to study astrocytic response during an inflammatory insult by palmitate through Flux Balance Analysis methods and data mining. In this aspect, we assessed the metabolic fluxes of human astrocytes under three different scenarios: healthy (normal conditions), induced inflammation by palmitate, and tibolone treatment under palmitate inflammation. Our results suggest that tibolone reduces the L-glutamate-mediated neurotoxicity in astrocytes through the modulation of several metabolic pathways involved in glutamate uptake. We also identified a set of reactions associated with the protective effects of tibolone, including the upregulation of taurine metabolism, gluconeogenesis, cPPAR and the modulation of calcium signaling pathways. In conclusion, the different scenarios studied in our model allowed us to identify several metabolic fluxes perturbed under an inflammatory response and the protective mechanisms exerted by tibolone.

Identification of biomarkers in common chronic lung diseases by co-expression networks and drug-target interactions analysis

Background

asthma, chronic obstructive pulmonary disease (COPD), and idiopathic pulmonary fibrosis (IPF) are three serious pulmonary diseases that contain common and unique characteristics. Therefore, the identification of biomarkers that differentiate these diseases is of importance for preventing misdiagnosis. In this regard, the present study aimed to identify the disorders at the early stages, based on lung transcriptomics data and drug-target interactions.

Methods

To this end, the differentially expressed genes were found in each disease. Then, WGCNA was utilized to find specific and consensus gene modules among the three diseases. Finally, the disease-disease similarity was analyzed, followed by determining candidate drug-target interactions.

Results

The results confirmed that the asthma lung transcriptome was more similar to COPD than IPF. In addition, the biomarkers were found in each disease and thus were proposed for further clinical validations. These genes included RBM42, STX5, and TRIM41 in asthma, CYP27A1, GM2A, LGALS9, SPI1, and NLRC4 in COPD, ATF3, PPP1R15A, ZFP36, SOCS3, NAMPT, and GADD45B in IPF, LRRC48 and CETN2 in asthma-COPD, COL15A1, GIMAP6, and JAM2 in asthma-IPF and LMO7, TSPAN13, LAMA3, and ANXA3 in COPD-IPF. Finally, analyzing drug-target networks suggested anti-inflammatory candidate drugs for treating the above mentioned diseases.

Conclusion

In general, the results revealed the unique and common biomarkers among three chronic lung diseases. Eventually, some drugs were suggested for treatment purposes.

Systems Bioinformatics: increasing precision of computational diagnostics and therapeutics through network-based approaches

Systems Bioinformatics is a relatively new approach, which lies in the intersection of systems biology and classical bioinformatics. It focuses on integrating information across different levels using a bottom-up approach as in systems biology with a data-driven top-down approach as in bioinformatics. The advent of omics technologies has provided the stepping-stone for the emergence of Systems Bioinformatics. These technologies provide a spectrum of information ranging from genomics, transcriptomics and proteomics to epigenomics, pharmacogenomics, metagenomics and metabolomics. Systems Bioinformatics is the framework in which systems approaches are applied to such data, setting the level of resolution as well as the boundary of the system of interest and studying the emerging properties of the system as a whole rather than the sum of the properties derived from the system's individual components. A key approach in Systems Bioinformatics is the construction of multiple networks representing each level of the omics spectrum and their integration in a layered network that exchanges information within and between layers. Here, we provide evidence on how Systems Bioinformatics enhances computational therapeutics and diagnostics, hence paving the way to precision medicine. The aim of this review is to familiarize the reader with the emerging field of Systems Bioinformatics and to provide a comprehensive overview of its current state-of-the-art methods and technologies. Moreover, we provide examples of success stories and case studies that utilize such methods and tools to significantly advance research in the fields of systems biology and systems medicine.

Metabolic network-based stratification of hepatocellular carcinoma reveals three distinct tumor subtypes

Hepatocellular carcinoma (HCC) is one of the most frequent forms of liver cancer, and effective treatment methods are limited due to tumor heterogeneity. There is a great need for comprehensive approaches to stratify HCC patients, gain biological insights into subtypes, and ultimately identify effective therapeutic targets. We stratified HCC patients and characterized each subtype using transcriptomics data, genome-scale metabolic networks and network topology/controllability analysis. This comprehensive systems-level analysis identified three distinct subtypes with substantial differences in metabolic and signaling pathways reflecting at genomic, transcriptomic, and proteomic levels. These subtypes showed large differences in clinical survival associated with altered kynurenine metabolism, WNT/β-catenin-associated lipid metabolism, and PI3K/AKT/mTOR signaling. Integrative analyses indicated that the three subtypes rely on alternative enzymes (e.g., ACSS1/ACSS2/ACSS3, PKM/PKLR, ALDOB/ALDOA, MTHFD1L/MTHFD2/MTHFD1) to catalyze the same reactions. Based on systems-level analysis, we identified 8 to 28 subtype-specific genes with pivotal roles in controlling the metabolic network and predicted that these genes may be targeted for development of treatment strategies for HCC subtypes by performing in silico analysis. To validate our predictions, we performed experiments using HepG2 cells under normoxic and hypoxic conditions and observed opposite expression patterns between genes expressed in high/moderate/low-survival tumor groups in response to hypoxia, reflecting activated hypoxic behavior in patients with poor survival. In conclusion, our analyses showed that the heterogeneous HCC tumors can be stratified using a metabolic network-driven approach, which may also be applied to other cancer types, and this stratification may have clinical implications to drive the development of precision medicine.

Metabolic Network-Based Identification and Prioritization of Anticancer Targets Based on Expression Data in Hepatocellular Carcinoma

Hepatocellular carcinoma (HCC) is a deadly form of liver cancer with high mortality worldwide. Unfortunately, the large heterogeneity of this disease makes it difficult to develop effective treatment strategies. Cellular network analyses have been employed to study heterogeneity in cancer, and to identify potential therapeutic targets. However, the existing approaches do not consider metabolic growth requirements, i.e., biological network functionality, to rank candidate targets while preventing toxicity to non-cancerous tissues. Here, we developed an algorithm to overcome these issues based on integration of gene expression data, genome-scale metabolic models, network controllability, and dispensability, as well as toxicity analysis. This method thus predicts and ranks potential anticancer non-toxic controlling metabolite and gene targets. Our algorithm encompasses both objective-driven and-independent tasks, and uses network topology to finally rank the predicted therapeutic targets. We employed this algorithm to the analysis of transcriptomic data for 50 HCC patients with both cancerous and non-cancerous samples. We identified several potential targets that would prevent cell growth, including 74 anticancer metabolites, and 3 gene targets (PRKACA, PGS1, and CRLS1). The predicted anticancer metabolites showed good agreement with existing FDA-approved cancer drugs, and the 3 genes were experimentally validated by performing experiments in HepG2 and Hep3B liver cancer cell lines. Our observations indicate that our novel approach successfully identifies therapeutic targets for effective treatment of cancer. This approach may also be applied to any cancer type that has tumor and non-tumor gene or protein expression data.

Integrating multi-source information on a single network to detect disease-related clusters of molecular mechanisms

The abundance of available information for each disease from multiple sources (e.g. as genetic, regulatory, metabolic, and protein-protein interaction) constitutes both an advantage and a challenge in identifying disease-specific underlying mechanisms. Integration of multi-source data is a rising topic and a great challenge in precision medicine and is crucial in enhancing disease understanding, identifying meaningful clusters of molecular mechanisms and increasing precision and personalisation towards the goal of Predictive, Preventive and Personalised Medicine (PPPM). The overall aim of this work was to develop a novel network-based integration methodology with the following characteristics: (i) maximise the number of data sources, (ii) utilise holistic approaches to integrate these sources (iii) be simple, flexible and extendable, (iv) be conclusive. Here, we present the case of Alzheimer's disease as a paradigm for illustrating our novel approach.

SIGNIFICANCE

In this work we present an integration methodology, which aggregates a large number of the available data sources and types by exploiting the holistic nature of network approaches. It is simple, flexible and extendable generating solid conclusions regarding the molecular mechanisms that underlie the input data. We have illustrated the strength of our proposed methodology using Alzheimer's disease as a paradigm. This method is expected to serve as a stepping-stone for further development of integration methods of multi-source omic-data and to contribute to progress towards the goal of Predictive, Preventive and Personalised Medicine (PPPM). The output of this methodology may act as a reference map of implicated pathways in the disease under investigation, where pathways related to additional omics data from any kind of experiment may be projected. This will increase the precision in the understanding of the disease and may contribute to personalised approaches for patients with different disease-related pathway profile, leading to a more precise, personalised and ideally preventive management of the disease.

Dynamic modeling of folliculogenesis signaling pathways in the presence of miRNAs expression

Background

TEK signaling plays a very important role in folliculogenesis. It activates Ras/ERK/MYC, PI3K/AKT/mTORC1 and ovarian steroidogenesis activation pathways. These are the main pathways for cell growth, differentiation, migration, adhesion, proliferation, survival and protein synthesis.

Results

TEK signaling on each of the two important pathways where levels of pERK, pMYC, pAkt, pMCL1 and pEIF4EBP1 are increased in dominant follicles and pMYC is decreased in dominant follicles. Over activation of ERK and MYC which are the main cell growth and proliferation and over activation of Akt, MCl1, mTORC1 and EIF4EBP1 which are the main cell survival and protein synthesis factors act as promoting factors for folliculogenesis. In case of over expression of hsa-miR-30d-3p and hsa-miR-451a, MYC activity level is considerably increased in subordinate follicles. Our simulation results show that in the presence of has-miR-548v and bta-miR-22-3p, downstream factors of pathways are inhibited.

Conclusions

Our work offers insight into the design of natural biological procedures and makes predictions that can guide further experimental studies on folliculogenesis pathways. Moreover, it defines a simple signal processing unit that may be useful for engineering synthetic biology and genes circuits to carry out cell-based computation.

Electronic supplementary material

The online version of this article (doi:10.1186/s13048-017-0371-y) contains supplementary material, which is available to authorized users.

MAMMOTh: A new database for curated mathematical models of biomolecular systems

MOTIVATION

Living systems have a complex hierarchical organization that can be viewed as a set of dynamically interacting subsystems. Thus, to simulate the internal nature and dynamics of the entire biological system, we should use the iterative way for a model reconstruction, which is a consistent composition and combination of its elementary subsystems. In accordance with this bottom-up approach, we have developed the MAthematical Models of bioMOlecular sysTems (MAMMOTh) tool that consists of the database containing manually curated MAMMOTh fitted to the experimental data and a software tool that provides their further integration.

RESULTS

The MAMMOTh database entries are organized as building blocks in a way that the model parts can be used in different combinations to describe systems with higher organizational level (metabolic pathways and/or transcription regulatory networks). The tool supports export of a single model or their combinations in SBML or Mathematica standards. The database currently contains 110 mathematical sub-models for Escherichia coli elementary subsystems (enzymatic reactions and gene expression regulatory processes) that can be combined in at least 5100 complex/sophisticated models concerning more complex biological processes as de novo nucleotide biosynthesis, aerobic/anaerobic respiration and nitrate/nitrite utilization in E. coli. All models are functionally interconnected and sufficiently complement public model resources.

AVAILABILITY

http://mammoth.biomodelsgroup.ru.

miRNA-mRNA network involved in folliculogenesis interactome: systems biology approach

At later phases of folliculogenesis, the mammalian ovarian follicle contains layers of granulosa cells surrounding an antral cavity. To better understand the molecular basis of follicular growth and granulosa cell maturation, we study transcriptome profiling of granulosa cells from small (<5 mm) and large (>10 mm) bovine follicles using simultaneous method of Affymetrix microarrays (24,128 probe sets) and RNA-Seq data sets. This study proposes a computational method to discover the functional miRNA-mRNA regulatory modules, that is, groups of miRNAs and their target mRNAs that are believed to take part cooperatively in post-transcriptional gene regulation under specific conditions. The reconstructed network was named Integrated miRNA-mRNA Bipartite Network. 277 genes and 6 key modules were disclosed through clustering for mRNA master list. The 66 genes are among the genes that belong to at least two modules. All these genes, being involved in at least one of the phenomena, namely cell survival, proliferation, metastasis and apoptosis, have an overexpression pattern (P < 0.01). For miRNA master list, a total of 172 sequences were differentially expressed (P < 0.01) between dominant (large) and each of subordinate (small) follicles. Within the follicle, these miRNAs were predominantly expressed in mural granulosa cells. Finally, predicted and validated targets of these miRNAs enriched in dominant (large) follicles were identified, which are mapped to signaling pathways involved in follicular cell proliferation, steroidogenesis, PI3K/AKT/mTOR and Ras/Raf/MEK/ERK. The identification of miRNAs and their target mRNAs and the construction of their regulatory networks may give new insights into biological procedures.

Genome-Scale Reconstruction of the Human Astrocyte Metabolic Network

Astrocytes are the most abundant cells of the central nervous system; they have a predominant role in maintaining brain metabolism. In this sense, abnormal metabolic states have been found in different neuropathological diseases. Determination of metabolic states of astrocytes is difficult to model using current experimental approaches given the high number of reactions and metabolites present. Thus, genome-scale metabolic networks derived from transcriptomic data can be used as a framework to elucidate how astrocytes modulate human brain metabolic states during normal conditions and in neurodegenerative diseases. We performed a Genome-Scale Reconstruction of the Human Astrocyte Metabolic Network with the purpose of elucidating a significant portion of the metabolic map of the astrocyte. This is the first global high-quality, manually curated metabolic reconstruction network of a human astrocyte. It includes 5,007 metabolites and 5,659 reactions distributed among 8 cell compartments, (extracellular, cytoplasm, mitochondria, endoplasmic reticle, Golgi apparatus, lysosome, peroxisome and nucleus). Using the reconstructed network, the metabolic capabilities of human astrocytes were calculated and compared both in normal and ischemic conditions. We identified reactions activated in these two states, which can be useful for understanding the astrocytic pathways that are affected during brain disease. Additionally, we also showed that the obtained flux distributions in the model, are in accordance with literature-based findings. Up to date, this is the most complete representation of the human astrocyte in terms of inclusion of genes, proteins, reactions and metabolic pathways, being a useful guide for in-silico analysis of several metabolic behaviors of the astrocyte during normal and pathologic states.

Neuroblastoma, a Paradigm for Big Data Science in Pediatric Oncology

Pediatric cancers rarely exhibit recurrent mutational events when compared to most adult cancers. This poses a challenge in understanding how cancers initiate, progress, and metastasize in early childhood. Also, due to limited detected driver mutations, it is difficult to benchmark key genes for drug development. In this review, we use neuroblastoma, a pediatric solid tumor of neural crest origin, as a paradigm for exploring "big data" applications in pediatric oncology. Computational strategies derived from big data science-network- and machine learning-based modeling and drug repositioning-hold the promise of shedding new light on the molecular mechanisms driving neuroblastoma pathogenesis and identifying potential therapeutics to combat this devastating disease. These strategies integrate robust data input, from genomic and transcriptomic studies, clinical data, and in vivo and in vitro experimental models specific to neuroblastoma and other types of cancers that closely mimic its biological characteristics. We discuss contexts in which "big data" and computational approaches, especially network-based modeling, may advance neuroblastoma research, describe currently available data and resources, and propose future models of strategic data collection and analyses for neuroblastoma and other related diseases.

Airway remodeling: Systems biology approach, from bench to bedside

Airway Remodeling, a patho-physiologic process, is considered as a key feature of chronic airway diseases. In recent years, our understanding of the complex diseases has increased significantly by the use of combined approaches, including systems biology, which may contribute to the development of personalized and predictive medicine approaches. Integrative analysis, along with the cooperation of clinicians, computer scientists, research scientists, and bench scientists, has become an important part of the experimental design and therapeutic strategies in the era of omics. The airway remodeling process is the result of the dysregulation of several signaling pathways that modulate the airway regeneration; therefore, high-throughput experiments and systems biology approach can help to understand this process better. The study reviews related literature and is consistent with the existing clinical evidence.