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.

This study aims to highlight SARS-COV-2 mutations which are associated with increased or decreased viral virulence. We utilize genetic data from all strains available from GISAID and countries’ regional information, such as deaths and cases per million, as well as COVID-19-related public health austerity measure response times. Initial indications of selective advantage of specific mutations can be obtained from calculating their frequencies across viral strains. By applying modelling approaches, we provide additional information that is not evident from standard statistics or mutation frequencies alone. We therefore, propose a more precise way of selecting informative mutations. We highlight two interesting mutations found in genes N (P13L) and ORF3a (Q57H). The former appears to be significantly associated with decreased deaths and cases per million according to our models, while the latter shows an opposing association with decreased deaths and increased cases per million. Moreover, protein structure prediction tools show that the mutations infer conformational changes to the protein that significantly alter its structure when compared to the reference protein.

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. 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.

Huntington’s disease is a rare neurodegenerative disease caused by a cytosine–adenine–guanine (CAG) trinucleotide expansion in the Huntingtin (HTT) gene. Although Huntington’s disease (HD) is well studied, the pathophysiological mechanisms, genes and metabolites involved in HD remain poorly understood. Systems bioinformatics can reveal synergistic relationships among different omics levels and enables the integration of biological data. It allows for the overall understanding of biological mechanisms, pathways, genes and metabolites involved in HD. The purpose of this study was to identify the differentially expressed genes (DEGs), pathways and metabolites as well as observe how these biological terms differ between the pre-symptomatic and symptomatic HD stages. A publicly available dataset from the Gene Expression Omnibus (GEO) was analyzed to obtain the DEGs for each HD stage, and gene co-expression networks were obtained for each HD stage. Network rewiring, highlights the nodes that change most their connectivity with their neighbors and infers their possible implication in the transition between different states. The CACNA1I gene was the mostly highly rewired node among pre-symptomatic and symptomatic HD network. Furthermore, we identified AF198444 to be common between the rewired genes and DEGs of symptomatic HD. CNTN6, DEK, LTN1, MST4, ZFYVE16, CEP135, DCAKD, MAP4K3, NUPL1 and RBM15 between the DEGs of pre-symptomatic and DEGs of symptomatic HD and CACNA1I, DNAJB14, EPS8L3, HSDL2, SNRPD3, SOX12, ACLY, ATF2, BAG5, ERBB4, FOCAD, GRAMD1C, LIN7C, MIR22, MTHFR, NABP1, NRG2, OTC, PRAMEF12, SLC30A10, STAG2 and Y16709 between the rewired genes and DEGs of pre-symptomatic HD. The proteins encoded by these genes are involved in various biological pathways such as phosphatidylinositol-4,5-bisphosphate 3-kinase activity, cAMP response element-binding protein binding, protein tyrosine kinase activity, voltage-gated calcium channel activity, ubiquitin protein ligase activity, adenosine triphosphate (ATP) binding, and protein serine/threonine kinase. Additionally, prominent molecular pathways for each HD stage were then obtained, and metabolites related to each pathway for both disease stages were identified. The transforming growth factor beta (TGF-β) signaling (pre-symptomatic and symptomatic stages of the disease), calcium (Ca²⁺) signaling (pre-symptomatic), dopaminergic synapse pathway (symptomatic HD patients) and Hippo signaling (pre-symptomatic) pathways were identified. The in silico metabolites we identified include Ca²⁺, inositol 1,4,5-trisphosphate, sphingosine 1-phosphate, dopamine, homovanillate and L-tyrosine. The genes, pathways and metabolites identified for each HD stage can provide a better understanding of the mechanisms that become altered in each disease stage. Our results can guide the development of therapies that may target the altered genes and metabolites of the perturbed pathways, leading to an improvement in clinical symptoms and hopefully a delay in the age of onset.

PathwayConnector is a web-tool that facilitates the construction of complementary pathway-to-pathway networks and subnetworks of them, based on a reference pathway network derived from the rich information available either in KEGG or Reactome database for pathway mapping. Specifically, for a given set of pathways, PathwayConnector (i) finds all the direct connections between them, (ii) adds a minimum set of complementary pathways required to achieve connectivity between the pathways, leading to informative fully connected networks and (ii) provides a series of clustering methods for the further grouping of pathways in to sub-clusters. The proposed web-tool is a simple yet informative tool towards identifying connected groups of pathways that are significantly related to specific diseases. http://bioinformatics.cing.ac.cy/PathwayConnector Supplementary data are available at Bioinformatics online.