Bhavjinder K. Dhillon
Maren Smith
Arjun Baghela
Amy H. Y. Lee
Robert E. W. Hancock- 1Centre for Microbial Diseases and Immunity Research, University of British Columbia, Vancouver, BC, Canada
- 2Molecular Biology & Biochemistry Department, Simon Fraser University, Burnaby, BC, Canada
Systems biology is an approach to interrogate complex biological systems through large-scale quantification of numerous biomolecules. The immune system involves >1,500 genes/proteins in many interconnected pathways and processes, and a systems-level approach is critical in broadening our understanding of the immune response to vaccination. Changes in molecular pathways can be detected using high-throughput omics datasets (e.g., transcriptomics, proteomics, and metabolomics) by using methods such as pathway enrichment, network analysis, machine learning, etc. Importantly, integration of multiple omic datasets is becoming key to revealing novel biological insights. In this perspective article, we highlight the use of protein-protein interaction (PPI) networks as a multi-omics integration approach to unravel information flow and mechanisms during complex biological events, with a focus on the immune system. This involves a combination of tools, including: InnateDB, a database of curated interactions between genes and protein products involved in the innate immunity; NetworkAnalyst, a visualization and analysis platform for InnateDB interactions; and MetaBridge, a tool to integrate metabolite data into PPI networks. The application of these systems techniques is demonstrated for a variety of biological questions, including: the developmental trajectory of neonates during the first week of life, mechanisms in host-pathogen interaction, disease prognosis, biomarker discovery, and drug discovery and repurposing. Overall, systems biology analyses of omics data have been applied to a variety of immunology-related questions, and here we demonstrate the numerous ways in which PPI network analysis can be a powerful tool in contributing to our understanding of the immune system and the study of vaccines.
Introduction
In the field of immunology, a systems biology approach is necessary to understanding the immune response to vaccination, infection and diseases, since these involve complex interactions between a large number of genetic, epigenetic, physiological and environmental factors. Systems-level strategies can ultimately be applied to better understand the molecular changes in humans upon exposure to a vaccine or an immunotherapeutic, to understand the mechanisms underlying disease or pathogenesis, and to characterize the effect(s) of specific challenges to the immune system (1–5). Omics technologies offer the ability to measure such aspects in an unbiased way that is high-throughput and cost-effective. Several omics methods have been employed in the context of systems vaccinology (3), including but not limited to, whole genome sequencing (genomics), RNA-Seq for measuring mRNA levels (transcriptomics), high-throughput mass spectrometry for measuring protein levels (proteomics) and metabolite levels (metabolomics), CHiP-Seq for determining transcription factor binding sites, ATAC-Seq to identify DNA modification sites (epigenomics), 16S rRNA sequencing for microbiota profiling (microbiomics), and equivalent omics analyses performed at the single-cell level. Recently, there has also been a growing effort to obtain multiple omics profiles in the same individuals, since shared insights across omics datasets strengthens links between underlying biological mechanisms and responses of interest, and can provide more reliable interpretation of gene function, higher-level changes and novel insights not observed in single-omics studies (6–8). Overall, biological samples can be manipulated to generate numerous omics datasets, and can be applied to study how our immune systems elicit effective, therapeutic and/or pathological responses.
A key challenge in systems biology is building the appropriate bioinformatics tools to integrate omics datasets, ultimately enabling the correlation of global changes with the underlying biological events that drove those changes. Statistical and machine learning approaches have been applied to omics datasets [reviewed previously (9–11)] to identify sets of molecular features that (i) are dysregulated/correlated with observed phenotypes, (ii) can be used as biomarkers to predict observed phenotypes, or (iii) can be targeted by drugs for improved therapies. A wide array of tools are available to run single- or multi-omics analysis pipelines (12), including commercial platforms and more recently published “self-serve” platforms [e.g., OmicsNet (13), OmicsPlayground (14)]. Typically, such methods interrogate information in either a supervised or unsupervised manner; supervised methods identify differences between labeled omics data from different conditions (e.g., responders vs. non-responders or treated vs. untreated) while, unsupervised methods reveal global patterns of gene dysregulation without any labels.
Downstream characterization of dysregulated molecules can further our understanding of underlying biological mechanisms at play. This can be achieved by interrogating curated functional genomics information from databases of gene ontologies (functional descriptions), pathways, known interactors, transcription factor binding sites (TFBS) upstream of dysregulated genes, etc. through various enrichment analyses. However, a large proportion of genes have not been assigned to canonical pathways in pathways databases (such as KEGG or Reactome), so pathway enrichment limits the ability of such approaches to reveal novel insights (15).
The use of biological networks is a powerful approach to integrate multi-omics data to identify novel biological insights (15–18). To characterize the role of individual molecular features in larger cellular processes and global changes using networks involves either overlaying omics data on experimentally-derived known networks [e.g., protein-protein interaction (PPI) networks], or by inferring networks directly from the data [e.g., co-expressed genes (19)], the strengths/limitations of which have been reviewed previously (15). A few commonly used biological networks along with related resources and tools are summarized in Table 1. The application of PPI networks to interrelate dysregulated genes is a very powerful method for revealing the systems-level flow of information through key hubs (highly connected protein nodes) and subnetworks. Because PPIs include direct, metabolic, and regulatory interactions between proteins, they essentially chart potentially biologically relevant, i.e., functional, interconnections. This can enable the determination of emergent properties, which are essentially new biological insights into the processes driving the observed transcriptional differences. The results from a PPI network analysis are always framed as hypotheses rather than knowledge per se, and must be eventually tested using downstream wet lab experiments (15).
Table 1. Examples of functional biological information that can be represented using networks, along with corresponding databases/repositories and supplementary data analysis tools that can be used to assess the functional data in high-throughput omics datasets.
In this article, we provide an overview of the philosophies and methodologies that can be employed in the analysis of omics data, especially with regards to integration of omics datasets using an unsupervised network analysis approach. Examples are provided of how such analyses enable novel hypothesis generation for: (a) immune system development, (b) mechanisms of host-pathogen interactions, (c) discovery of mechanism-based biomarkers, and (d) strategies to define prospective new interventions based on drug repurposing. While the methods are somewhat biased toward the study of innate immune and inflammatory responses, it is worth mentioning that “innate immunity instructs adaptive immunity” (63) in that (i) the effectors of adaptive immunity are often innate immune mechanisms, (ii) many of the pathways involved are the same, and (iii) vaccine adjuvants that improve adaptive immune responses boost innate immunity. Therefore, the tools we describe have value in investigating adaptive immunity as well as human genetic diseases/conditions with an underlying inflammatory pathology.
Systems Tools for Network-Based Analyses Using PPIs
InnateDB (43, 44) and other International Molecular Exchange (IMEx) consortium databases (42) provide the basis for understanding biological connections in cells according to known interactions between molecular elements, such as proteins. InnateDB is a publicly available database (www.innatedb.com) focused on elucidating the genes, proteins, and molecular “interactome” of the innate immune response, with an emphasis on curation of experimentally-validated PPIs and signaling pathways in human, mouse and bovine. The interactome can be used to understand the interplay between multi-omics datasets that measure different parts of a larger system of physical, metabolic, and regulatory networks. For example, human TRAF6 and MyD88 are usually defined as having a role in the major TLR4 to NFκB signaling pathway of innate immunity. However, in InnateDB, they are experimentally documented to interact with 398 and 129 proteins, respectively, in humans. This means that there is a massive potential for these proteins to bridge and/or participate in multiple biological pathways when activated by innate immune stimuli.
InnateDB is an important tool in immunology as evidenced by the >6,000,000 hits from more than 55,000 visitors annually. While all known pathways (>3,500) and molecular interactions (318,000 in human) are present, the emphasis on innate immunity is achieved through the contextual review, curation and annotation of molecular interactions and pathways involved in innate immunity. To date, the InnateDB curation team has reviewed more than 5,200 publications annotating >27,000 molecular interactions of >9,400 separate genes in rich detail including annotation of the cell, cell-line and tissue type; the molecules involved; the interaction detection method; etc. By including interaction and pathway data relevant to all biological processes, a much broader perspective of innate immunity can be achieved, especially since an effective innate immune response requires the coordinated efforts of many important processes including the endocrine, circulatory, and nervous systems (64). Additionally, it becomes possible to investigate any biological signaling process of interest beyond the immune system, as well as inflammation and adaptive immunity.
InnateDB facilitates systems-level analyses by enabling the integration, analysis and visualization of user-supplied quantitative data, such as gene expression data, in the context of molecular interaction networks and pathways. This includes the statistically robust analysis of overrepresented pathways, interactomes, ontologies, TFBS, and networks. One can, for example, refine the network to show only molecular interactions between a list of differentially expressed (DE) genes (and their encoded products) or view all potential interactors regardless of whether they are DE. This can aid in the identification of important nodes that may not be regulated transcriptionally or which are expressed at an earlier or later time. Networks derived from InnateDB can be interactively visualized using the Cerebral plug-in for Cytoscape (65) to generate biologically intuitive, pathway-like layouts of networks, or in a more recently developed tool, NetworkAnalyst (50–52). NetworkAnalyst is an extremely fast network analysis and visualization tool for the analysis of gene expression data in the context of PPI networks. In addition, MetaBridge (21) is a tool that can be used for the integration of metabolite-protein interactions into these existing networks. In combination, these tools can be used to perform multi-omics integration of transcriptomics, proteomics, and metabolomics data in an unsupervised manner.
In addition to these outlined methods, there are bioinformatics tools available for performing other types of network analyses specifically for studying the immune system. Examples include immuneExpresso (66), a data mining tool built as part of Immport to capture inter-cell interactions, and Ontogenet (67), a component of the ImmGen database enabling construction of gene regulatory networks based on sets of co-expressed genes. Such tools can be useful in revealing novel inter-cell interactions or regulatory factors, respectively, but ultimately may be too limited in scope for a systems-level analysis. Thus, we focus here on how PPI-based network analysis tools can be applied to better understand human health and disease.
Mechanistic Insights Into Human Immune Development
Most recently, as a part of the EPIC-HIPC consortium, we published a study that revealed a robust developmental trajectory of immune ontogeny during the first week of life in newborns using a multi-omics integration approach (9). Transcriptomic, proteomic, and metabolomic data were derived from <1 ml of blood collected from West African (The Gambia) neonates at two time points: day of life (DOL) 0 and a second DOL, either 1, 3, or 7.
Importantly, through this study, we were able to show that multi-omics integration using PPI networks (through NetworkAnalyst, InnateDB, and MetaBridge) provided similar biological insights, but greater depth, when compared to data-driven supervised integration approaches [namely, DIABLO (68) and Multifactorial Response Network (MMRN) (69, 70)]. Major observations from this study revealed that the first week of life is highly dynamic; DOL0 and DOL1 were quite similar with few DE genes, but by DOL3, 1,125 DE genes were detected, and 1,864 DE genes by DOL7. These represented several key pathways in immune development, mainly centered around interferon signaling, the complement cascade, and neutrophil activity. These have previously been shown to play a role in the newborn immune response to infection, but until this study were not identified as central to ontogeny in the first week of life. Importantly, these pathways and nearly 60% of transcriptomic changes were confirmed in a second independent cohort of neonates from Papua New Guinea/Australasia, revealing that neonatal immune development is not random, but follows a precise and possibly purposeful age-specific path.
An unsupervised PPI network was used to integrate the transcriptomic, metabolomics, and proteomic data to reveal a single functional network, highlighting that individual omics datasets are complementary, reporting different facets of the same biological processes. For example, both the transcriptomic and proteomic data confirmed the increase in type I interferon-related functions and the regulation of complement cascades. Importantly, this integration also revealed novel nodes in the PPI network that were not identified by any single-omics dataset on its own, representing novel biological insights, including changes in cellular replication machinery, creatinine metabolism, fibrin clotting cascade, adaptive immunity markers and phagosome activity.
Thus, these systems biology approaches allowed novel insights into the immune developmental trajectory during the first week of life in newborns. Further studies are being conducted to provide insights into the mechanistic differences in the susceptibility of neonates to infection-related disease or death during this critical phase of life. Also, in the context of vaccinology, an integrative systems biology approach is being used to reveal mechanistic insights into the molecular determinants of vaccination efficacy, while taking into account this developmental trajectory.
Mechanistic Insights Into Host-Pathogen Interactions
Systems biology methods have also been leveraged to study host-pathogen interactions (71). One example is of infection by the obligate human intracellular pathogen Chlamydia trachomatis, the major cause of bacterial sexually-transmitted diseases (STDs) and preventable blindness worldwide. This involved a study that coupled transcriptomics and proteomics to assess the macrophage responses to infection with C. trachomatis (72). Macrophages were derived from human induced pluripotent stem cells (iPSdMs), which share >95% similarity in terms of gene expression with primary human blood monocyte-derived macrophages, and were able to support the growth of C. trachomatis intracellularly to mimic infection in-vitro.
Pathway analysis of 2,029 DE genes (from transcriptomics) and 307 DE proteins (from proteomics) at 24 h post-infection, revealed strong interferon α, β, and γ responses, and dysregulation of various Toll-like receptor pathways, the endosomal/vacuolar pathway, energy metabolism, and metabolism of amino acids and nucleotides and inhibition of translation. Most significantly upregulated were genes associated with type I interferon signaling, including key transcription factors such as interferon regulatory factors (IRF)-1, 3, and 7, which are known to contribute to the regulation of type I interferons during Chlamydia infection.
Importantly, IRF5 and IL-10RA, not previously characterized for their role in Chlamydia infection, were identified as key players in limiting infection in macrophages. Indeed, IRF5−/− and IL-10RA−/− mutant iPSDM cells were both shown to have increased susceptibility to C. trachomatis infection. These results, along with numerous other published studies [e.g., (73–77)], demonstrate that multi-omics integration using PPI networks can reveal novel insight into the factors that play a significant role in the host immune response to infections.
Mechanism-Based Biomarkers for Disease Diagnosis and Prognosis Prediction
Systems biology analyses have also led to insights into mechanisms underlying disease prognosis and prediction of diagnostic biomarkers. One such study of the enteric pathogen Salmonella enterica sv. Typhimurium (78) involved the use of transcriptomics to compare gene expression in HIV patients with and without severe invasive non-typhoidal Salmonella (iNTS) infections, as well as HIV patients with other acute bacterial infections (including E. coli and Streptococcus pneumoniae). Initially, 1,200 genes were upregulated in HIV patients with iNTS and with other acute bacterial infections, compared to HIV patients without a bacterial infection. However, genes upregulated in patients with non-Salmonella acute infections showed enrichment for pathways typically associated with innate immune/inflammatory responses, while conversely the gene expression response in patients with iNTS could be explained by upregulation of genes that are associated with suppression of inflammation (NFKBIB, PI3K, REL, SIGIRR, SOCS4, SOCS7). This lack of innate immune response and viral signature, which was subsequently shown to be consistent with increased viral load (79), leading to insights into the poor prognosis of HIV patients with iNTS.
These types of analyses were also used to explore immune manipulation using host defense (antimicrobial) peptides. Such peptides selectively modulate the innate immune response and protect against infection, and are produced by many organisms to defend against infections (80). Furthermore, novel small innate defense regulator (IDR) peptides have been shown to be effective in animal models against antibiotic resistant bacteria, tuberculosis, cerebral malaria, pre-term birth and inflammation (81, 82). To better understand the cellular cascade that occurs after these IDR peptides enter the cell, transcriptional changes were assessed in human monocytes and peripheral blood mononuclear cells (83). The biological relevance of these gene expression changes was assessed using pathway over-representation, TFBS analysis, and network analysis with NetworkAnalyst, implicating 11 pathways including the p38, Erk1/2, and JNK mitogen-activated (MAP)-kinases, NFκB, two Src family kinases, and more than 15 transcription factors [including NFκB (most subunits), Creb, IRF4, AP-1, AP-2, Are, E2F1, SP1, Gre, and STAT3]. NetworkAnalyst showed that some of the top connected hub proteins within networks constructed from dysregulated genes were involved in the functioning of MAP kinases and induction of chemokines, anti-inflammatory pathways particularly TGFβ, and type I interferon responses. These highly connected hubs reveal mechanistic insights and could potentially represent diagnostic or treatment biomarkers. Ultimately, a similar approach can be utilized to evaluate any agent perturbing cellular function, including immunomodulators and vaccines, and can define biomarkers differentiating between responders and non-responders.
Drug Discovery and Repurposing
Systems biology techniques have been applied to aid in drug discovery and repurposing of existing agents for the treatment of cancers, bacterial and viral infections, and genetic disorders (84). One such study aimed at finding better therapeutics for cystic fibrosis (CF) utilized transcriptomics to study immortalized CFTR−/− (cystic fibrosis transmembrane regulator) epithelial cells stimulated for hyperinflammation, a state known to lead to deterioration of lung function in CF patients (85). Genes differentially expressed between CFTR−/− cells and corrected variants were submitted to InnateDB for analysis and integration with PPI networks. This revealed the interconnectivity of the CFTR and innate immune networks through the PRKAA1 (AMP kinase)/AKT1 and HSPB1 pathways. Genes within this network were then submitted to DrugBank (86), allowing for the identification of the diabetes drug Metformin as an AMP kinase activator, which was then tested in-vitro and shown to reduce inflammation by ~50%. DE genes between CFTR−/− cells and corrected variants also included 54 genes involved in autophagy. In disease states, autophagy is an adaptive response to stress that favors infection survival and resolution (87). Follow up studies confirmed that CFTR mutant cells demonstrated arrested autophagy. It was then demonstrated that the antimicrobial peptide IDR-1018 resolved this arrested autophagy state and reduced inflammation. These genes also revealed a strong upregulation of ER stress and unfolded protein response pathways, through activation of the IRE-1 pathway (88). Follow up studies showed that salubrinal, an inhibitor of negative regulator GADD34, upregulated this pathway and suppressed inflammation. Thus, through these systems biology-based studies, novel pharmaceuticals (IDR-1018) and 2 existing drugs (Metformin and salubrinal) were identified as potential treatments for CF-related hyperinflammation. As such, along with numerous other studies [e.g., (89–92)], it has been shown that integrating omics datasets using resources such as InnateDB and DrugBank can reveal potential drug targets for improved therapies.
Discussion and the Future
The analyses outlined in this article merely scratch the surface of what is possible using systems biology and high-throughput omics techniques to study the immune system, e.g., the major tools described here (43, 44, 50–52) have been used and cited more than 1,500 times. The above-described examples highlight that using unbiased multi-omics experiments in conjunction with incisive bioinformatics tools, such as PPI network integration, one can go beyond the hypothesis-testing scientific method to use unbiased omics data to generate fundamentally new hypotheses and develop new biological insights. Ultimately such studies should lead to the development of novel diagnostics, individualized therapies for diseases and vaccines. Furthermore, systems biology approaches can provide invaluable insights to inform the stratification of individuals with the same syndrome but different underlying mechanisms, the diagnosis of disease and/or flare-ups, ongoing development of new vaccines and/or adjuvants as well as immune-based therapeutics providing insights into the optimal strategies for delivery of interventions.
Author Contributions
BD, MS, and AB all contributed to writing of the first draft of this article. AL performed data analysis of the ontogeny study and provided valuable feedback through the writing process. RH supervised all authors, edited the manuscript and provided critical insights and feedback. All authors contributed to the article and approved the submitted version.
Funding
Our bioinformatics research is currently supported by a grant from the Canadian Institutes for Health Research FDN-154287, and previously received funding from Genome Canada, Genome BC, and the Foundation for the National Institutes of Health through their Grand Challenges in Global Health Research program. RH holds a Canada Research Chair and a UBC Killam Professorship.
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The authors wish to acknowledge collaborators Fiona Brinkman, David Lynn, Jeff Xia, Tamara Munzner, and previous lab members Erin Gill, Chris Fjell, and Jennifer Gardy as well as the InnateDB curation team for their fantastic and critical contributions.
Abbreviations
CF, cystic fibrosis; CFTR, cystic fibrosis transmembrane regulator; DE, differentially expressed; DOL, day of life; IDR, innate defense regulator; iNTS, invasive non-typhoidal Salmonella; MAP, mitogen-activated; PPI, protein-protein interaction; TFBS, transcription factor binding site.
References
1. Trautmann L, Sekaly R. Solving vaccine mysteries: a systems biology perspective. Nat Immunol. (2011) 12:729. doi: 10.1038/ni.2078
2. Mooney M, McWeeney S, Canderan G, Sékaly R. A systems framework for vaccine design. Curr Opin Immunol. (2013) 25:551–5. doi: 10.1016/j.coi.2013.09.014
3. Pulendran B, Li S, Nakaya HI. Systems vaccinology. Immunity. (2010) 33:516–29. doi: 10.1016/j.immuni.2010.10.006
4. Oberg AL, Kennedy RB, Li P, Ovsyannikova IG, Poland GA. Systems biology approaches to new vaccine development. Curr Opin Immunol. (2011) 23:436–43. doi: 10.1016/j.coi.2011.04.005
5. Kotliarov Y, Sparks R, Martins AJ, Mulè MP, Lu Y, Goswami M, et al. Broad immune activation underlies shared set point signatures for vaccine responsiveness in healthy individuals and disease activity in patients with lupus. Nat Med. (2020) 26:618–29. doi: 10.1038/s41591-020-0769-8
6. Ge H, Walhout AJ, Vidal M. Integrating ‘omic'information: a bridge between genomics and systems biology. Trends Genet. (2003) 19:551–60. doi: 10.1016/j.tig.2003.08.009
7. Lee AH, Shannon CP, Amenyogbe N, Bennike TB, Diray-Arce J, Idoko OT, et al. Dynamic molecular changes during the first week of human life follow a robust developmental trajectory. Nat Commun. (2019) 10:1092. doi: 10.1038/s41467-019-08794-x
8. Ebrahim A, Brunk E, Tan J, O'brien EJ, Kim D, Szubin R, et al. Multi-omic data integration enables discovery of hidden biological regularities. Nat Commun. (2016) 7:1–9. doi: 10.1038/ncomms13091
9. Tolios A, De Las Rivas J, Hovig E, Trouillas P, Scorilas A, Mohr T. Computational approaches in cancer multidrug resistance research: Identification of potential biomarkers, drug targets and drug-target interactions. Drug Resist Updates. (2020) 48:100662. doi: 10.1016/j.drup.2019.100662
10. Kidd BA, Peters LA, Schadt EE, Dudley JT. Unifying immunology with informatics and multiscale biology. Nat Immunol. (2014) 15:118–27. doi: 10.1038/ni.2787
11. Tavassoly I, Goldfarb J, Iyengar R. Systems biology primer: the basic methods and approaches. Essays Biochem. (2018) 62:487–500. doi: 10.1042/EBC20180003
12. Beale DJ, Karpe AV, Ahmed W. Beyond metabolomics: a review of multi-omics-based approaches. Microb Metab. (2016) 289−312. doi: 10.1007/978-3-319-46326-1_10
13. Zhou G, Xia J. OmicsNet: a web-based tool for creation and visual analysis of biological networks in 3D space. Nucleic Acids Res. (2018) 46:W514–22. doi: 10.1093/nar/gky510
14. Akhmedov M, Martinelli A, Geiger R, Kwee I. Omics playground: a comprehensive self-service platform for visualization, analytics and exploration of big omics data. NAR Genom Bioinform. (2020) 2:lqz019. doi: 10.1093/nargab/lqz019
15. Charitou T, Bryan K, Lynn DJ. Using biological networks to integrate, visualize and analyze genomics data. Genet Select Evol. (2016) 48:27. doi: 10.1186/s12711-016-0205-1
16. Barabasi A, Oltvai ZN. Network biology: understanding the cell's functional organization. Nat Rev Genet. (2004) 5:101. doi: 10.1038/nrg1272
17. Saint-Antoine MM, Singh A. Network inference in systems biology: recent developments, challenges, and applications. Curr Opin Biotechnol. (2020) 63:89–98. doi: 10.1016/j.copbio.2019.12.002
18. Mardinoglu A, Boren J, Smith U, Uhlen M, Nielsen J. Systems biology in hepatology: approaches and applications. Nat Rev Gastroenterol Hepatol. (2018) 15:365–77. doi: 10.1038/s41575-018-0007-8
19. Costa RL, Boroni M, Soares MA. Distinct co-expression networks using multi-omic data reveal novel interventional targets in HPV-positive and negative head-and-neck squamous cell cancer. Sci Rep. (2018) 8:1–13. doi: 10.1038/s41598-018-33498-5
20. Kanehisa M, Goto S. KEGG: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res. (2000) 28:27–30. doi: 10.1093/nar/28.1.27
21. Hinshaw SJ, Lee AHY, Gill EE, Hancock REW. MetaBridge: enabling network-based integrative analysis via direct protein interactors of metabolites. Bioinformatics. (2018) 34:3225–7. doi: 10.1093/bioinformatics/bty331
22. Croft D, Mundo AF, Haw R, Milacic M, Weiser J, Wu G. The reactome pathway knowledgebase. Nucleic Acids Res. (2014) 42:D472–7. doi: 10.1093/nar/gkt1102
23. Mi H, Thomas P. PANTHER pathway: an ontology-based pathway database coupled with data analysis tools. Protein Netw Pathw Anal. (2009) 563:123–40. doi: 10.1007/978-1-60761-175-2_7
24. Yu G, He Q. ReactomePA: an R/bioconductor package for Reactome pathway analysis and visualization. Mol BioSyst. (2016) 12:477–9. doi: 10.1039/C5MB00663E
25. Gene Ontology Consortium. The Gene Ontology (GO) database and informatics resource. Nucleic Acids Res. (2004) 32 (Suppl_1):D258–61. doi: 10.1093/nar/gkh036
26. Foroushani AB, Brinkman FS, Lynn DJ. Pathway-GPS and SIGORA: identifying relevant pathways based on the over-representation of their gene-pair signatures. PeerJ. (2013) 1:e229. doi: 10.7717/peerj.229
27. Ma H, Sorokin A, Mazein A, Selkov A, Selkov E, Demin O, et al. The Edinburgh human metabolic network reconstruction and its functional analysis. Mol Syst Biol. (2007) 3:135. doi: 10.1038/msb4100177
28. Zur H, Ruppin E, Shlomi T. iMAT: an integrative metabolic analysis tool. Bioinformatics. (2010) 26:3140–2. doi: 10.1093/bioinformatics/btq602
29. Brunk E, Sahoo S, Zielinski DC, Altunkaya A, Dräger A, Mih N, et al. Recon3D enables a three-dimensional view of gene variation in human metabolism. Nat Biotechnol. (2018) 36:272. doi: 10.1038/nbt.4072
30. Agren R, Bordel S, Mardinoglu A, Pornputtapong N, Nookaew I, Nielsen J. Reconstruction of genome-scale active metabolic networks for 69 human cell types and 16 cancer types using INIT. PLoS Comp Biol. (2012) 8:e1002518. doi: 10.1371/journal.pcbi.1002518
31. Blais EM, Rawls KD, Dougherty BV, Li ZI, Kolling GL, Ye P, et al. Reconciled rat and human metabolic networks for comparative toxicogenomics and biomarker predictions. Nat Commun. (2017) 8:1–15. doi: 10.1038/ncomms14250
32. Wang Y, Eddy JA, Price ND. Reconstruction of genome-scale metabolic models for 126 human tissues using mCADRE. BMC Syst Biol. (2012) 6:153. doi: 10.1186/1752-0509-6-153
33. Malik-Sheriff RS, Glont M, Nguyen TV, Tiwari K, Roberts MG, Xavier A, et al. BioModels - 15 years of sharing computational models in life science. Nucleic Acids Res. (2020) 48:D407–15. doi: 10.1093/nar/gkz1055
34. ENCODE Project Consortium. The ENCODE (ENCyclopedia of DNA elements) project. Science. (2004) 306:636–40. doi: 10.1126/science.1105136
35. Ecker JR, Bickmore WA, Barroso I, Pritchard JK, Gilad Y, Segal E. ENCODE explained. Nature. (2012) 489:52–4. doi: 10.1038/489052a
36. Mathelier A, Fornes O, Arenillas DJ, Chen CY, Denay G, Lee J, et al. JASPAR 2016: a major expansion and update of the open-access database of transcription factor binding profiles. Nucleic Acids Res. (2016) 44:D110–5. doi: 10.1093/nar/gkv1176
37. Matys V, Fricke E, Geffers R, Gößling E, Haubrock M, Hehl R, et al. TRANSFAC®: transcriptional regulation, from patterns to profiles. Nucleic Acids Res. (2003) 31:374–8. doi: 10.1093/nar/gkg108
38. Meyer PE, Kontos K, Lafitte F, Bontempi G. Information-theoretic inference of large transcriptional regulatory networks. EURASIP J Bioinform Syst Biol. (2007) 2007:1–9. doi: 10.1155/2007/79879
39. Margolin AA, Nemenman I, Basso K, Wiggins C, Stolovitzky G, Dalla Favera R, et al. ARACNE: an algorithm for the reconstruction of gene regulatory networks in a mammalian cellular context. BMC Bioinform. (2006) 7:S7. doi: 10.1186/1471-2105-7-S1-S7
40. Verfaillie A, Imrichová H, Van de Sande B, Standaert L, Christiaens V, Hulselmans G, et al. iRegulon: from a gene list to a gene regulatory network using large motif and track collections. PLoS Comp Biol. (2014) 10:e1003731. doi: 10.1371/journal.pcbi.1003731
41. Geurts P. dynGENIE3: dynamical GENIE3 for the inference of gene networks from time series expression data. Sci Rep. (2018) 8:1–12. doi: 10.1038/s41598-018-21715-0
42. Orchard S, Kerrien S, Abbani S, Aranda B, Bhate J, Bidwell S, et al. Protein interaction data curation: the international molecular exchange (IMEx) consortium. Nat Methods. (2012) 9:345–50. doi: 10.1038/nmeth.1931
43. Lynn DJ, Winsor GL, Chan C, Richard N, Laird MR, Barsky A, et al. InnateDB: facilitating systems-level analyses of the mammalian innate immune response. Mol Syst Biol. (2008) 4:218. doi: 10.1038/msb.2008.55
44. Breuer K, Foroushani AK, Laird MR, Chen C, Sribnaia A, Lo R, et al. InnateDB: systems biology of innate immunity and beyond—recent updates and continuing curation. Nucleic Acids Res. (2012) 41:D1228–33. doi: 10.1093/nar/gks1147
45. Bader GD, Betel D, Hogue CW. BIND: the biomolecular interaction network database. Nucleic Acids Res. (2003) 31:248–50. doi: 10.1093/nar/gkg056
46. Salwinski L, Miller CS, Smith AJ, Pettit FK, Bowie JU, Eisenberg D. The database of interacting proteins: 2004 update. Nucleic Acids Res. (2004) 32 (Suppl_1):D449–51. doi: 10.1093/nar/gkh086
47. Licata L, Briganti L, Peluso D, Perfetto L, Iannuccelli M, Galeota E, et al. MINT, the molecular interaction database: 2012 update. Nucleic Acids Res. (2012) 40:D857–61. doi: 10.1093/nar/gkr930
48. Orchard S, Ammari M, Aranda B, Breuza L, Briganti L, Broackes-Carter F, et al. The MIntAct project - IntAct as a common curation platform for 11 molecular interaction databases. Nucleic Acids Res. (2014) 42:D358–63. doi: 10.1093/nar/gkt1115
49. Oughtred R, Stark C, Breitkreutz B, Rust J, Boucher L, Chang C, et al. The BioGRID interaction database: 2019 update. Nucleic Acids Res. (2019) 47:D529–41. doi: 10.1093/nar/gky1079
50. Xia J, Benner MJ, Hancock REW. NetworkAnalyst-integrative approaches for protein–protein interaction network analysis and visual exploration. Nucleic Acids Res. (2014) 42:W167–74. doi: 10.1093/nar/gku443
51. Xia J, Gill EE, Hancock REW. NetworkAnalyst for statistical, visual and network-based meta-analysis of gene expression data. Nat Protoc. (2015) 10:823. doi: 10.1038/nprot.2015.052
52. Zhou G, Soufan O, Ewald J, Hancock REW, Basu N, Xia J. NetworkAnalyst 3.0: a visual analytics platform for comprehensive gene expression profiling and meta-analysis. Nucleic Acids Res. (2019) 47:W234–41. doi: 10.1093/nar/gkz240
53. Liu X, Chang C, Han M, Yin R, Zhan Y, Li C, et al. PPIExp: a web-based platform for integration and visualization of Protein-Protein interaction data and spatiotemporal proteomics data. J Proteome Res. (2018) 18:633–41. doi: 10.1021/acs.jproteome.8b00713
54. Gough NR. Science's signal transduction knowledge environment: the connections maps database. Ann N Y Acad Sci. (2002) 971:585–7. doi: 10.1111/j.1749-6632.2002.tb04532.x
55. Krull M, Pistor S, Voss N, Kel A, Reuter I, Kronenberg D, et al. TRANSPATH®: an information resource for storing and visualizing signaling pathways and their pathological aberrations. Nucleic Acids Res. (2006) 34 (Suppl_1):D546–51. doi: 10.1093/nar/gkj107
56. Law V, Knox C, Djoumbou Y, Jewison T, Guo AC, Liu Y, et al. DrugBank 4.0: shedding new light on drug metabolism. Nucleic Acids Res. (2014) 42:D1091–7. doi: 10.1093/nar/gkt1068
57. Qin C, Zhang C, Zhu F, Xu F, Chen SY, Zhang P, et al. Therapeutic target database update 2014: a resource for targeted therapeutics. Nucleic Acids Res. (2014) 42:D1118–23. doi: 10.1093/nar/gkt1129
58. Hecker N, Ahmed J, von Eichborn J, Dunkel M, Macha K, Eckert A, et al. SuperTarget goes quantitative: update on drug-target interactions. Nucleic Acids Res. (2012) 40:D1113–7. doi: 10.1093/nar/gkr912
59. Kuhn M, Szklarczyk D, Pletscher-Frankild S, Blicher TH, Von Mering C, Jensen LJ, et al. STITCH 4: integration of protein-chemical interactions with user data. Nucleic Acids Res. (2014) 42:D401–7. doi: 10.1093/nar/gkt1207
60. Gaulton A, Bellis LJ, Bento AP, Chambers J, Davies M, Hersey A, et al. ChEMBL: a large-scale bioactivity database for drug discovery. Nucleic Acids Res. (2012) 40:D1100–7. doi: 10.1093/nar/gkr777
61. Liu T, Lin Y, Wen X, Jorissen RN, Gilson MK. BindingDB: a web-accessible database of experimentally determined protein-ligand binding affinities. Nucleic Acids Res. (2007) 35 (Suppl_1):D198–201. doi: 10.1093/nar/gkl999
62. Yamanishi Y, Kotera M, Moriya Y, Sawada R, Kanehisa M, Goto S. DINIES: drug-target interaction network inference engine based on supervised analysis. Nucleic Acids Res. (2014) 42:W39–45. doi: 10.1093/nar/gku337
63. Jain A, Pasare C. Innate control of adaptive immunity: beyond the three-signal paradigm. J Immunol. (2017) 198:3791–800. doi: 10.4049/jimmunol.1602000
64. Gardy JL. Enabling a systems biology approach to immunology: focus on innate immunity. Trends Immunol. (2009) 30:249–62. doi: 10.1016/j.it.2009.03.009
65. Barsky A, Gardy JL, Hancock REW, Munzner T. Cerebral: a cytoscape plugin for layout of and interaction with biological networks using subcellular localization annotation. Bioinformatics. (2007) 23:1040–2. doi: 10.1093/bioinformatics/btm057
66. Bhattacharya S, Dunn P, Thomas CG, Smith B, Schaefer H, Chen J, et al. ImmPort, toward repurposing of open access immunological assay data for translational and clinical research. Sci Data. (2018) 5:180015. doi: 10.1038/sdata.2018.15
67. Shay T, Kang J. Immunological genome project and systems immunology. Trends Immunol. (2013) 34:602–9. doi: 10.1016/j.it.2013.03.004
68. Singh A, Shannon CP, Gautier B, Rohart F, Vacher M, Tebbutt SJ, et al. DIABLO: an integrative approach for identifying key molecular drivers from multi-omics assays. Bioinformatics. (2019) 35:3055–62. doi: 10.1093/bioinformatics/bty1054
69. Chaussabel D, Baldwin N. Democratizing systems immunology with modular transcriptional repertoire analyses. Nat Rev Immunol. (2014) 14:271–80. doi: 10.1038/nri3642
70. Li S, Sullivan NL, Rouphael N, Yu T, Banton S, Maddur MS, et al. Metabolic phenotypes of response to vaccination in humans. Cell. (2017) 169:862–77. e17. doi: 10.1016/j.cell.2017.04.026
71. Yeung A, Hale C, Clare S, Palmer S, Scott JB, Baker S, et al. Using a systems biology approach study host-pathogen interactions. Bacteria Intracell. (2019) 18:337–47. doi: 10.1128/microbiolspec.BAI-0021-2019
72. Yeung AT, Hale C, Lee AH, Gill EE, Bushell W, Parry-Smith D, et al. Exploiting induced pluripotent stem cell-derived macrophages to unravel host factors influencing chlamydia trachomatis pathogenesis. Nat Commun. (2017) 8:1–12. doi: 10.1038/ncomms15013
73. Elmassry MM, Mudaliar NS, Colmer-Hamood JA, et al. New markers for sepsis caused by Pseudomonas aeruginosa during burn infection. Metabolomics. (2020) 16:1–16. doi: 10.1007/s11306-020-01658-2
74. Baschal EE, Larson ED, Bootpetch Roberts TC, et al. Identification of novel genes and biological pathways that overlap in infectious and nonallergic diseases of the upper and lower airways using network analyses. Front Genet. (2020) 10:1352. doi: 10.3389/fgene.2019.01352
75. Sun X, Hua S, Gao C, Blackmer JE, Ouyang Z, Ard K, et al. Immune-profiling of ZIKV-infected patients identifies a distinct function of plasmacytoid dendritic cells for immune cross-regulation. Nat Commun. (2020) 11:1–13. doi: 10.1038/s41467-020-16217-5
76. Smith J. Immunological molecular responses of human retinal pigment epithelial cells to infection with Toxoplasma gondii. Front Immunol. (2019) 10:708. doi: 10.3389/fimmu.2019.00708
77. Mulindwa J, Matovu E, Enyaru J, Clayton C. Blood signatures for second stage human African trypanosomiasis: a transcriptomic approach. BMC Med Genom. (2020) 13:1–12. doi: 10.1186/s12920-020-0666-5
78. Schreiber F, Lynn DJ, Houston A, Peters J, Mwafulirwa G, Finlay BB, et al. The human transcriptome during nontyphoid Salmonella and HIV coinfection reveals attenuated NFκB-mediated inflammation and persistent cell cycle disruption. J Infect Dis. (2011) 204:1237–45. doi: 10.1093/infdis/jir512
79. Preziosi MJ, Kandel SM, Guiney DG, Browne SH. Microbiological analysis of nontyphoidal Salmonella strains causing distinct syndromes of bacteremia or enteritis in HIV/AIDS patients in San Diego, California. J Clin Microbiol. (2012) 50:3598–603. doi: 10.1128/JCM.00795-12
80. Scott MG, Dullaghan E, Mookherjee N, Glavas N, Waldbrook M, Thompson A, et al. An anti-infective peptide that selectively modulates the innate immune response. Nat Biotechnol. (2007) 25:465–72. doi: 10.1038/nbt1288
81. Mansour SC, de la Fuente-Núñez C, Hancock REW. Peptide IDR-1018: modulating the immune system and targeting bacterial biofilms to treat antibiotic-resistant bacterial infections. J Peptide Sci. (2015) 21:323–9. doi: 10.1002/psc.2708
82. Wu BC, Lee AHY, Hancock REW. Mechanisms of the innate defense regulator peptide-1002 anti-inflammatory activity in a sterile inflammation mouse model. J Immunol. (2017) 199:3592–603. doi: 10.4049/jimmunol.1700985
83. Mookherjee N, Hamill P, Gardy J, Blimkie D, Falsafi R, Chikatamarla A, et al. Systems biology evaluation of immune responses induced by human host defence peptide LL-37 in mononuclear cells. Mol BioSyst. (2009) 5:483–96. doi: 10.1039/b813787k
84. Hopkins AL. Network pharmacology: the next paradigm in drug discovery. Nat Chem Biol. (2008) 4:682. doi: 10.1038/nchembio.118
85. Mayer ML, Blohmke CJ, Falsafi R, Fjell CD, Madera L, Turvey SE, et al. Rescue of dysfunctional autophagy attenuates hyperinflammatory responses from Cystic Fibrosis cells. J Immunol. (2013) 190:1227–38. doi: 10.4049/jimmunol.1201404
86. Knox C, Law V, Jewison T, Liu P, Ly S, Frolkis A, et al. DrugBank 3.0: a comprehensive resource for ‘omics' research on drugs. Nucleic Acids Res. (2010) 39 (Suppl_1):D1035–41. doi: 10.1093/nar/gkq1126
87. Saha S, Panigrahi DP, Patil S, Bhutia SK. Autophagy in health and disease: a comprehensive review. Biomed Pharmacother. (2018) 104:485–95. doi: 10.1016/j.biopha.2018.05.007
88. Blohmke CJ, Mayer ML, Tang AC, Hirschfeld AF, Fjell CD, Sze MA, et al. Atypical activation of the unfolded protein response in cystic fibrosis airway cells contributes to p38 MAPK-mediated innate immune responses. J Immunol. (2012) 189:5467–75. doi: 10.4049/jimmunol.1103661
89. Cheng F, Desai RJ, Handy DE, Wang R, Schneeweiss S, Barabási AL, et al. Network-based approach to prediction and population-based validation of in silico drug repurposing. Nat Commun. (2018) 9:1–12. doi: 10.1038/s41467-018-05116-5
90. Cheng F, Murray JL, Zhao J, Sheng J, Zhao Z, Rubin DH. Systems biology-based investigation of cellular antiviral drug targets identified by gene-trap insertional mutagenesis. PLoS Comp Biol. (2016) 12:e1005074. doi: 10.1371/journal.pcbi.1005074
91. Han Z, Xue W, Tao L, Zhu F. Identification of novel immune-relevant drug target genes for Alzheimer's disease by combining ontology inference with network analysis. CNS Neurosci Ther. (2018) 24:1253–63. doi: 10.1111/cns.13051
Keywords: systems biology, multi-omic integration, transcriptomics, innate immunity, immune ontogeny, host-pathogen interaction, drug discovery and repurposing, systems vaccinology
Citation: Dhillon BK, Smith M, Baghela A, Lee AHY and Hancock REW (2020) Systems Biology Approaches to Understanding the Human Immune System. Front. Immunol. 11:1683. doi: 10.3389/fimmu.2020.01683
Received: 11 March 2020; Accepted: 24 June 2020;
Published: 30 July 2020.
Edited by:
Jay Evans, University of Montana, United StatesReviewed by:
Susana Magadan, University of Vigo, SpainGeert Leroux-Roels, Ghent University, Belgium
Copyright © 2020 Dhillon, Smith, Baghela, Lee and Hancock. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Robert E. W. Hancock, bob@hancocklab.com