This article has Open Peer Review reports available.
Representing and querying disease networks using graph databases
- Artem Lysenko†1,
- Irina A. Roznovăţ†2Email author,
- Mansoor Saqi†2,
- Alexander Mazein2,
- Christopher J Rawlings1 and
- Charles Auffray2
© The Author(s). 2016
Received: 16 December 2015
Accepted: 30 June 2016
Published: 25 July 2016
Systems biology experiments generate large volumes of data of multiple modalities and this information presents a challenge for integration due to a mix of complexity together with rich semantics. Here, we describe how graph databases provide a powerful framework for storage, querying and envisioning of biological data.
We show how graph databases are well suited for the representation of biological information, which is typically highly connected, semi-structured and unpredictable. We outline an application case that uses the Neo4j graph database for building and querying a prototype network to provide biological context to asthma related genes.
Our study suggests that graph databases provide a flexible solution for the integration of multiple types of biological data and facilitate exploratory data mining to support hypothesis generation.
KeywordsDisease management platform Graph database Neo4j graph Protein-centric framework Systems medicine Computational approach
A major effort in translational medicine is to understand the molecular basis of disease [1, 2]. Analysis of high-throughput experimental data together with patient phenotypic information has led to the identification of sets of candidate genes, proteins and pathways that may be implicated in many disease conditions.
However, in order to build a higher level picture of the underlying processes involved in the disease pathology, it is necessary to integrate various classes of heterogeneous information, and to explore the complex relationships between entities such as diseases, candidate genes, proteins, interactions and pathways. These relationships could have a variety of different types, like participation (protein A participates in Pathway X), sequence similarity (protein A is sequence similar to protein B) or protein interaction (protein A interacts with protein B), and typically, will be associated with some provenance information. Generally, the relationships will describe complex networks. Some genes or proteins, for example, may be associated with multiple diseases, while others may be implicated in only one or two conditions. Additionally, some candidate genes may encode large families of sequence-similar proteins, and may thus be involved in multiple protein-protein interactions. Identification of these wider relationships between entities can help in providing biological context and in hypothesis generation on disease mechanisms. Traversing paths through the generated network can suggest hitherto unexpected relationships between entities, (e.g. relationships between disease conditions).
Biological information is typically highly connected, semi-structured and unpredictable, and these features are important in deciding an appropriate way to represent and query disease-relevant biomedical networks. Representing multiple types of relationships will generally lead to highly connected networks. Additionally, the amount of information associated with a given entity will not be uniform. Some proteins have been extensively investigated experimentally, are well annotated and associated with multiple complementary pieces of information. For example, a gene may be involved in a protein-protein interaction with another protein, may be a drug target for a particular drug, may be known to be enhanced in a particular tissue type or may participate in a particular pathway. Other proteins can have little or no information associated with them and are annotated only as ‘hypothetical proteins’. Therefore, such datasets do not naturally fit into a classic relational database model which would have been originally designed to handle application cases with a large number of uniform entries associated with a limited number of data types and representing relatively few relationships in each data model. Although highly connected and sparse networks can be stored in a relational database, generally traversal-type queries (joins), which connect data linked by different relationships, become too computationally expensive and cumbersome to design. Traversal type queries are important for hypothesis generation, as they can reveal paths connecting entities that might not have been expected, and might not be apparent from visual inspection of the network. Another issue is that biological data are often semi-structured, i.e. lacking a structured data model, while still possessing some semantic annotation. As relational databases require all data to be transformed to conform to a pre-defined schema, developing parsers for such semi-structured data can be particularly problematic. Biological research is also unpredictable in that new types of information can emerge with little advance notice – for example when new analytical instruments or new data resources are released. In a relational database management system, capturing these new data types often means re-design of a database schema, which is generally a complex and costly process. However, in a research setting this information is often of great importance and can impact on our understanding of the roles of genes and pathways in disease networks. A framework that facilitates the rapid and seamless inclusion of new data types would therefore be particularly well-suited to support the dynamic and ever-changing requirements of biomedical research community. Graph databases are a natural way to represent highly connected, non-uniformly distributed, semi-structured and unpredictable data as found in many biological systems studies [3, 4]. They also offer agile and flexible solutions and easily allow the inclusion of new data types.
Recently, there has been much interest in network representations of biological systems, particularly for network visualisation and network analytics. Several successful frameworks have been established and used to support bioinformatics analysis including, but not limited to, Cytoscape , Gephi  and NetworkX . However, the tools currently available are primarily for visualisation and analysis of network properties and do not offer an in-built support for semantics, graph querying or the ability to effectively work with datasets that are too large to visualise. Graph databases can naturally complement such tools and methods, by offering this much-needed extra functionality.
Several graph-based approaches have been developed for biological data integration such as Biozon , BN++ , Ondex [10, 11], and BioMine ). Some approaches, although providing network views, use an underlying relational database (e.g. BN++). The Ondex platform is a data warehouse using a graph data model and has been applied to plant systems biology, a domain where the data is typically fragmented across a large number of data sources and often poorly annotated. More recently, Bio4j  uses a graph database (Neo4j) to integrate data from several major repositories such as Uniprot KB , Gene Ontology (GO) , RefSeq , NCBI Taxonomy , and Expasy Enzyme DB .
The need to move towards considering graphs as biological knowledge repositories is particularly evident in the biomedical domain. Many important diseases have now been extensively researched for decades and the amounts of knowledge accumulated are both vast and heterogeneous. The large volume of interconnected data means that its exploration using a network visualizer cannot be carried out easily and the heterogeneity of the data can make cost-effective management difficult in a relational database system.
In this work, we have explored the potential of using a graph database to facilitate data management and analysis in order to provide biological context to disease-related genes and proteins. This application case used the well-established Neo4j graph database for building and querying a prototype disease map for a complex disease, such as asthma. Neo4j has a free community edition that can be run as a lightweight desktop database as well as a database server, and an SQL-like versatile query language for graphs called Cypher. The database comes with an immediately available easy-to-use web interface that can display the results as both graphs and tables and also has interfaces to R, Python and Java - core languages of modern bioinformatics research.
We have integrated information from human protein-protein interactions, pathways, sequence similarity, disease-gene, gene-tissue and protein-drug associations. The selected examples are intended to illustrate how the Neo4j Cypher query language can be used to perform complex queries on this integrated dataset. We also show how information from gene expression studies can be easily added to the database and how the database can provide biological context for these differentially expressed genes. The main advantage of presented approach is in its ability to contextualise disease-associated genes by facilitating the identification and visual exploration of their network neighbourhood. In addition, as we have already explained, graph databases excel in traversal-type queries, facilitating the exploration of chains of connected concepts. This ‘link discovery’ can lead to new hitherto unexpected relationships being identified. Therefore, in our application case we have aimed to explore both of these data mining strategies, namely discovering connections between different entities of biomedical interest (queries in Listings 1, 2 and 5) and characterising the relevant context of diseases and biomarkers (queries in Listings 3 and 4).
The knowledgebase we constructed is organised according to the following principles. First, we integrated several typical multi-omics data (sequence similarity, transcriptomics metabolic pathways and protein-protein interactions), which provide a low-level perspective on the relationships between biological entities (genes, proteins and drugs). Then we have added the more general information that relates these lower-level entities to the higher-level processes – e.g. different diseases. The queries presented here are intended to provide illustrative examples about how to leverage connections across these two levels – for example, by discovering links in the multi-omics layer relevant to understanding the differences between different related diseases or using that information for suggesting drug targets.
A graph database offers several options for representation of data. Connections between entities (i.e. nodes) can be encoded directly as edges linking them together (e.g. protein A interacts with protein B, where the edge between the nodes represents the ‘interacts with’ relationship). Additionally, properties can be associated with nodes (e.g. protein A has annotation X, where ‘has annotation X’ is a property of the node representing protein A), and nodes that share a given property could be understood to be linked. And lastly, nodes can be linked via edges to an intermediate node indicating a common property (e.g. membership in the same pathway, for example node i representing protein A is connected to node j, representing pathway P by an edge that describes the relationship ‘participates in’). For this application case, the choice was guided primarily by the concerns of efficiency in capturing the structure of the data. For example, a property is an appropriate choice when annotation is simple, like a category label. An edge can have properties of its own, so can capture situations where the relationship is qualified, like a protein-protein interaction with given methods (e.g. 2-hybrid, anti tag co-immunoprecipitation) and associated with a confidence score.
Data sources imported
All data used in this application case were imported into the Neo4j v2.3.1 graph database. The integration was done from a protein-centric perspective, where the nodes corresponding to the individual entries in the human subset of UniProt/Swissprot served as key anchoring points for all other types of annotation. As many of the databases we used are continuously updated, the point of reference when all data were acquired for this study is 27/03/2015. Sequence similarity was determined by running all versus all Tera-BlastPTM for the UniProt Human reference proteome on the Timelogic DeCypher system. The significance threshold was 1e − 05, with effective database size adjusted to 16.9 bln bases. The results were further processed by retaining bi-directional hits only and modelled in the graph database as sequence similarity edges. The set of sequence similarity data was refined to exclude duplicates based on the UniProt id pair: the similarity between protein1-protein2 and protein2-protein1 was represented by a single, undirected edge. A score was assigned to each edge, calculated as a negative log10 of the average of the two e-values; only hits for which the ratio between the alignment length and the shortest sequence was above 0.60 were included.
Protein-protein interaction data were imported from the EBI IntAct database  (PSI-MI tab format). Each interaction edge was annotated with a list of experimental methods supporting it, as well as a confidence score (as given by the IntAct database). Data about human metabolic pathways were acquired from the Reactome database (see for example ). Individual pathways were modelled as nodes, with edges linking them to proteins to indicate pathway membership. Data about protein association with particular diseases were modelled in a similar way. The diseases, as defined by the curated subset of DisGeNET  (2014 release), were likewise represented as nodes with edges connecting them to corresponding proteins; (mapping gene names to UniProt identifiers was done via the Restful web services API of the UniProt database). Data on approved drug compounds were downloaded from DrugBank  on 27/03/2015 as DrugBank csv files. Where such information was available, those compounds were linked to proteins they target via UniProt identifier references. Data on gene expression-tissue enhanced association have been downloaded from the Human Protein Atlas (HPA), (version 13), , and integrated in the database. The tissue types were represented by nodes in the graph database and the enhanced expression in particular tissue types were modelled as edges, with RNA-TS-FPKM values as property. Information on ENSEMBL ids, (used for gene identification within HPA), was mapped to UniProt ids using the ‘Retrieve/ID mapping’ tool from the UniProt database, and noted as a property of the Protein node.
We have included a set of high-confidence manually curated genes believed to be implicated in one of the respiratory diseases as identified by Kaneko et al. . The review categorized genes into four partially overlapping subsets of Bronchial Asthma (BA), Chronic Obstructive Pulmonary Disease (COPD), Essential hypertension (E-HTN) and Tuberculosis (TB). TB and E-HTN are included in the list as controls . The gene names given in the paper were manually resolved to the UniProt identifiers and imported as a list attribute of corresponding annotations. This resulted in 219 genes (104 BA, 58 COPD, 35 TB, 54 E-HTN); detailed information on these gene sets is given in Additional file 1: Table S4. This type of information, gathered typically from domain experts or from careful study of the literature, is important in disease network construction as it serves as a starting point to explore pathways and processes that may be associated with the disease phenotype.
Overall occurrences and types of nodes and relationships included in the graph database
(a) Occurrences of node types
(b) Types and occurrences of relationships
Data modelling strategy
Unlike the relational database, where tables and relationships are heavily influenced by technical constraints formalized as “normalization rules”, graph databases offer much more flexibility and are more closely guided by the needs of the intended application. For this application case we have chosen the entities that correspond to either logical or real entities of biomedical interest. From the multi-omics perspective, the two physical entities central to our application case are proteins and drugs. These can be members of multiple thematic sets – e.g. the proteins can be grouped by disease involvement, tissue and metabolic pathway and drugs can be linked to particular diseases. To allow us to explore the potentially complex inter-relationships between these sets, they are also represented as nodes and therefore can be referenced directly as part of the query pattern. The edges are used to indicate both i) the membership of an entity in a particular set (e.g. protein participates in a pathway) and ii) existence of relationships between individual entities (e.g. protein interacts with other protein).
The level of the detail was chosen in accordance with the need of our intended application case. For example, in theory it may be valid to model an interaction as a concept rather than an edge and then link it to specific tissues to indicate whether it can occur there. However, by doing so an extra node and two edges will be introduced into the graph, which will increase the computational complexity of querying the database and make query construction more challenging for the user. In this case, due to the applied nature of the work presented here, it was therefore more appropriate to opt for a more generalized representation. However, we recognize that often it may be necessary to develop detailed representations in order to capture provenance and in the case of broad-content rather than thematically focused resources. To demonstrate how such a scenario can be accommodated in Neo4j, we have modelled transcriptomics component of our dataset in a way that captures all information about its provenance (GEOStudy node), experimental design (GEOStudy to GEOComparison relationship) and the analysis used (GEOComparison node) to process the data.
Here we illustrate the use of a graph database for facilitating hypothesis generation on disease mechanisms. We describe four queries related to asthma, a complex respiratory disease whose symptoms include airway inflammation and remodelling. Asthma is associated with multiple phenotypes and co-morbidities and has similarities to other respiratory diseases such as COPD. A significant number of asthma patients do not respond to conventional corticosteroid therapy, suggesting a complex underlying aetiology of this disease at both physiological and molecular levels.
Identification of common proteins between Asthma and COPD, TB and E-HTN
Uniprot id set (Kaneko)
Uniprot id set (Disgenet)
[P29474, P04040, P12821, P07550, P30711]
Pulmonary Dis- Ease, Chronic Obstructive
[P01375, O00206, P01137, P35225, P01584, P29474, Q9BZ11, P07550, P12821, P09211, Q96QV1, P09488, P16410, P05305]
[P35228, P09601, P01137, P14780, P01375]
[P01579, P01920, P01911, P01584, P22301, P29460, P05112, P13501, P11473, P01375, O60603, Q9NR96]
Providing biological context: comparison of the network neighbourhoods of differentially expressed genes from two asthma studies
Retrieval of drug associations with asthma-related proteins and other sequence similar proteins
Results obtained can be used to infer new associations between drugs and asthma biomarkers. For example, although no direct association has been reported for the P08727 protein (asthma biomarker) and Tenecteplase (drug), an indirect relationship can be identified between these two entities given that the P05787 protein (KRT8 gene) has sequence similarity with the P08727 protein, which was reported as a target for Tenecteplase, (Fig. 7). Exploring such indirect relationships can be of interest in drug development research. The ENTREZ gene name and the protein name for Fig. 7 are provided in Additional file 2: Table S5c.
Graph traversal queries for exploring relationships between concepts
Gene symbols and corresponding UniProt identifiers, (resolved via the “Retrieve/ID mapping” tool), for the core circadian components
Durrington et al.  suggest the possibility of the disruption of clock genes affecting airway inflammation through a possible link of clock proteins to the immune response involving REV-ErbA-alpha. The graph network can be queried for P20393 (REV-ErbA-alpha) which does not appear on the shortest paths network in Fig. 8, but lies 3 steps from O00327(ARNTL/BMAL1) and Q99743(NPAS2), (Cypher query given in Listings 6b, in Additional file 5 and output results in Fig. 9b).
Discussion and conclusions
High volume datasets are being routinely generated in systems biomedicine. An on-going concern for the interpretation of these data is that biological information has a challenging mix of complexity and rich semantics. Modelling of data in graph database ideally are guided by the same principles as those used for ontology design. In a way, the ontology may play the same role as a schema plays in a relational database by defining a valid set of relationships and/or properties that can exist for particular entities or their combinations. However, for this application case we have elected not to explore this topic in great detail in order to keep focus firmly on the applied aspects of managing biomedical data using the Neo4j database and to present it at a level accessible to a more general audience. For completeness, we would like to point out that from the ontology-centric data modelling perspective a Neo4j graph can be thought of as a collection of instances of data, where node “Labels” are equivalent to classes, and types of edges – to relationship types. Another important consideration in choosing how to represent the data is the performance constraints. [Note: We would encourage developers to look at the Neo4j features relating to Cartesian product warning and to query flow diagnosis, using the EXPLAIN and PROFILE features provided.] In this respect, as we were specifically interested in exploring the applications of the link discovery type of queries in a biomedical setting, it was necessary to ensure both the compactness of paths and adequate annotation necessary to ensure that the traversal space can be computed, e.g. it does not produce an all-versus-all set of entities as part of computing the solution (a “Cartesian product”).
In traditional relational databases approaches, the technical structure often gets in the way of exploratory data analysis either by visualisation or through data mining techniques. Graph databases (explored here with Neo4j) offer a powerful but lightweight, intuitive and flexible solution that is more readily compatible with a biological network view, a representation that is now near-ubiquitously employed for modelling complex biological data. This study shows how a combination of a powerful query language (to narrow down the search space) with network visualisation can help data exploration and hypothesis generation. The use of graph database approaches has some similarities with previous systems for logic-based modelling of integrated datasets [29, 30].
As we have illustrated in this work, graph databases offer many features that make them a particularly attractive option for a research-based setting, in particular in cases where it might be necessary to dynamically develop and interactively mine heterogeneous and not uniformly inter-related data. However, we would like to stress that in our opinion graph databases are best thought of as a complement to relational database technologies. There will clearly be cases where a relational database solution will be preferable, especially when data are dense and naturally fit the tabular representation where relational databases would offer much better processing performance.
Another important consideration is the place of Neo4j in a wider family of NoSQL and especially graph databases. The concept of graph databases is still relatively novel and therefore currently there are multiple alternative systems under active development. However, one standard – Resource Description Framework (RDF) and its enabling technologies (triple stores and SPARQL query language) have gained particular prominence due to their importance for Semantic Web and Internet-powered federated data solutions. For this reason we would like to further elaborate on the differences between these two standards and specifically the roles they could potentially play in management of large biomedical datasets.
Neo4j database belongs to the family of “Property graph” databases, which in contrast to RDF have conceptually distinct property elements that can be attached to both nodes and edges. In some respects, this syntax may be more familiar to users of object-oriented languages, as in contrast to RDF, there is neither a global uniqueness constraint nor specifically required format for node identifiers. In combination, these features make graph query syntax more compact and easier to read. Given that the Neo4j database is inherently a “closed world” solution, it is possible to collect the entirety of meta-data about the types and number of different entities and relationships between them, which can then inform further exploration. As SPARQL and RDF are primarily designed to operate in an “open world” setting of the Semantic Web, creating such summaries may not always be possible or can be very computationally intensive.
However, we would like to acknowledge that the RDF formalism does have its undisputed advantages in facilitating integration and exchange of distributed data and issues like lack of properties on edges can be addressed by leveraging more advanced features like reification or named graphs. In our view, the matter of preferences between the properties databases like Neo4j or RDF triple stores comes down to whether a particular application case can actually benefit from the extra features available via Semantic Web technologies. As adopting such technologies will inevitably come with considerable costs in terms of technical overhead and increased complexity, property databases may be preferable when speedy and flexible development of a single non-federated resource is the ultimate requirement.
The lack of schema in Neo4j may also be somewhat of a double-edge sword: although it does offer much flexibility, it also removes the interoperability standard from the data, which may make sharing and management of large distributed projects datasets more challenging. One possible solution to this would be to re-introduce some consistent semantics by binding critical attribute names and/or values of nodes and edges to a consensus controlled vocabulary. In a biomedical domain, this controlled vocabulary will likely come in the form of a suitable ontology. The ontologies could then provide a de-facto data modelling standard and many key concepts in the biomedical domain have already been modelled by projects like SNOMED , MeSH , GO  and OBI , with several hundred of other more specialized ontologies available from the OBO foundry  and NCBO Bioportal  repositories. Although at present the availability of advanced ontology-based features, like reasoning and data validation, varies in different graph databases, some notable solutions do offer examples of very close integration - e.g. the already-mentioned graph databases based on the resource description framework (RDF) .
We have explored here the utility of a graph database in providing a powerful yet flexible solution for disease network representations and for exploratory data mining and analysis to support hypothesis generation on disease mechanisms. Our study suggests that the Neo4j system offers a level of performance and an appropriate query and visualisation interface to effectively mine and manipulate these data. In particular, we have found the Cypher graph query language to be of great utility, because, as illustrated by our sample queries, it enabled us to generate a representative selection of common biologically-motivated queries with minimal efforts. As systems medicine projects continue to generate large amounts of heterogeneous datasets, graph database approaches may offer useful solutions for their knowledge management.
BA, bronchial asthma; COPD, chronic obstructive pulmonary disease; DEG, differentially expressed genes; E-HTN, essential hypertension; GEO, Gene Expression Omnibus (database); GSM, samples from the GEO database; HPA, Human Protein Atlas; MiA, mild asthma; MMA, mild-moderate asthma; NC, normal control; RDF, the resource description framework; SA, severe asthma; TB, tuberculosis
The ENTREZ gene names and the protein names are given in Additional file 2.
Authors would like to thank the team of the CNRS/IN2P3 Computing Centre, Mr. B. Guillon, Mr. J. Bussery, Dr. G. Marchetti and Dr. G. Rahal for their support on deploying the Neo4j graph database presented here to a server and making it accessible. Authors would like also to thank Mr. R. Henkel (University of Rostock, Germany) for his valuable comments on the Neo4j aspects and Dr. C. Marshall (BioSciConsulting, Belgium) for his valuable advice on the licenses of the interrogated databases. Authors are also grateful to reviewers for their valuable comments on the manuscript. Authors would like to acknowledge for access to the Neo4j framework and for data resources within the databases interrogated: DisGeNET 2.1, DrugBank, Human Protein Atlas, IntAct, Reactome, UniProt.
The research leading to these results has received support from the Innovative Medicines Initiative Joint Undertaking under grant agreement no. IMI 115010, (U-BIOPRED) and IMI 115446 (eTRIKS) resources of which are composed of financial contribution from the European Union’s Seventh Framework Programme (FP7/2007–2013) and EFPIA companies. AL and CJR acknowledge support from the BBSRC through their strategic funding of Rothamsted Research.
Availability of data and materials
The framework, described above, was developed predominantly in Java, (eclipse), using the Neo4j 2.3.1 functionality. The queries described here can be explored at the following url: https://diseaseknowledgebase.etriks.org/protein/browser/. The framework is freely available for non-commercial purposes and the java code used to parse the data and populate the graph database is available at https://github.com/ibalaur/ProteinFramework.git. For further details please contact firstname.lastname@example.org.
MS and AL conceived of the project, IR and AL implemented the approach. AM provided mappings for gene identifiers, MS, AL, CJR, IR wrote the first draft of the manuscript, CA and CJR provided coordination and guidance; all authors revised and approved the final version of the manuscript.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Auffray C, Charron D, Hood L. Predictive, preventive, personalized and participatory medicine: back to the future. Genome Med. 2010;2:57. doi:10.1186/gm178.View ArticlePubMedPubMed CentralGoogle Scholar
- Hood L, Tian Q. Systems approaches to biology and disease enable translational systems medicine. Genomics Proteomics Bioinformatics. 2012;10:181–5. doi:10.1016/j.gpb.2012.08.004.View ArticlePubMedGoogle Scholar
- Callahan A, Cruz-Toledo J, Ansell P, Dumontier M. Bio2RDF release 2: improved coverage, interoperability and provenance of life science linked data. In: Cimiano P, Corcho O, Presutti V, et al., editors. Semantic Web Semant. Berlin Heidelberg: Big Data. Springer; 2013. p. 200–12.Google Scholar
- Pareja-Tobes P, Tobes R, Manrique M, et al. Bio4j: a high-performance cloud-enabled graph-based data platform. bioRxiv 016758. 2015. doi: http://dx.doi.org/10.1101/016758.
- Smoot ME, Ono K, Ruscheinski J, et al. Cytoscape 2.8: new features for data integration and network visualization. Bioinformatics. 2011;27(3):431–32. doi:110.1093/bioinformatics/btq675.
- Bastian M, Heymann S, Jacomy M. Gephi: an open source software for exploring and manipulating networks. ICWSM. 2009;8:361–62.Google Scholar
- Schult DA, Swart P. Exploring network structure, dynamics, and function using NetworkX. In Proceedings of the 7th Python in Science Conferences (SciPy 2008). 2008;2008:11–6.Google Scholar
- Birkland A, Yona G. BIOZON: a system for unification, management and analysis of heterogeneous biological data. BMC Bioinformatics. 2006;7(1)1. doi:10.1186/1471-2105-7-70.
- Küntzer J, Blum T, Gerasch A, Backes C, Hildebrandt A, Kaufmann M, Lenhof HP. BN++-a biological information system. Journal of Integrative Bioinformatics. 2006;3(2)34.Google Scholar
- Köhler J, Baumbach J, Taubert J, et al. Graph-based analysis and visualization of experimental results with ONDEX. Bioinforma Oxf Engl. 2006;22:1383–90. doi:10.1093/bioinformatics/btl081.View ArticleGoogle Scholar
- Lysenko A, Lysenko A, Hindle MM, et al. Data integration for plant genomics–exemplars from the integration of Arabidopsis thaliana databases. Brief Bioinform. 2009;10:676–93. doi:10.1093/bib/bbp047.View ArticlePubMedGoogle Scholar
- Eronen L, Toivonen H. Biomine: predicting links between biological entities using network models of heterogeneous databases. BMC Bioinformatics. 2012;13:119. doi:10.1186/1471-2105-13-119.View ArticlePubMedPubMed CentralGoogle Scholar
- Consortium TU. UniProt: a hub for protein information. Nucleic Acids Res. 2015;43:D204–12. doi:10.1093/nar/gku989.View ArticleGoogle Scholar
- Consortium TGO. Gene ontology consortium: going forward. Nucleic Acids Res. 2015;43:D1049–56. doi:10.1093/nar/gku1179.View ArticleGoogle Scholar
- Tatusova T, Ciufo S, Fedorov B, et al. RefSeq microbial genomes database: new representation and annotation strategy. Nucleic Acids Res. 2014;42:D553–559. doi:10.1093/nar/gkt1274.View ArticlePubMedGoogle Scholar
- Sayers EW, Barrett T, Benson DA, et al. Database resources of the National Center for Biotechnology Information. Nucleic Acids Res. 2009;37:D5–D15. doi:10.1093/nar/gkn741.View ArticlePubMedGoogle Scholar
- Bairoch A. The ENZYME database in 2000. Nucleic Acids Res. 2000;28:304–5.View ArticlePubMedPubMed CentralGoogle Scholar
- Kerrien S, Aranda B, Breuza L, et al. The IntAct molecular interaction database in 2012. Nucleic Acids Res. 2012;40:D841–846. doi:10.1093/nar/gkr1088.View ArticlePubMedGoogle Scholar
- Croft D, Mundo AF, Haw R, et al. The Reactome pathway knowledgebase. Nucleic Acids Res. 2014;42:D472–477. doi:10.1093/nar/gkt1102.View ArticlePubMedGoogle Scholar
- Bauer-Mehren A, Bundschus M, Rautschka M, et al. Gene-disease network analysis reveals functional modules in mendelian, complex and environmental diseases. PLoS One. 2011;6:e20284. doi:10.1371/journal.pone.0020284.View ArticlePubMedPubMed CentralGoogle Scholar
- Knox C, Law V, Jewison T, et al. DrugBank 3.0: a comprehensive resource for “Omics” research on drugs. Nucleic Acids Res. 2011;39:D1035–41. doi:10.1093/nar/gkq1126.View ArticlePubMedGoogle Scholar
- Uhlén M, Fagerberg L, Hallström BM, et al. Proteomics. Tissue-based map of the human proteome. Science. 2015;347:1260419. doi:10.1126/science.1260419.View ArticlePubMedGoogle Scholar
- Kaneko Y, Yatagai Y, Yamada H, et al. The search for common pathways underlying asthma and COPD. Int J Chron Obstruct Pulmon Dis. 2013;8:65–78. doi:10.2147/COPD.S39617.PubMedPubMed CentralGoogle Scholar
- Voraphani N, Gladwin MT, Contreras AU, et al. An airway epithelial iNOS-DUOX2-thyroid peroxidase metabolome drives Th1/Th2 nitrative stress in human severe asthma. Mucosal Immunol. 2014;7:1175–85. doi:10.1038/mi.2014.6.View ArticlePubMedPubMed CentralGoogle Scholar
- Modena BD, Tedrow JR, Milosevic J, et al. Gene expression in relation to exhaled nitric oxide identifies novel asthma phenotypes with unique biomolecular pathways. Am J Respir Crit Care Med. 2014;190:1363–72. doi:10.1164/rccm.201406-1099OC.View ArticlePubMedPubMed CentralGoogle Scholar
- Ritchie ME, Phipson B, Wu D, et al. (2015) limma powers differential expression analyses for RNA-sequencing and microarray studies. Nucleic Acids Res gkv007. doi: 10.1093/nar/gkv007.
- Durrington HJ, Farrow SN, Loudon AS, Ray DW. The circadian clock and asthma. Thorax. 2014;69:90–2. doi:10.1136/thoraxjnl-2013-203482.View ArticlePubMedGoogle Scholar
- Ko CH, Takahashi JS. Molecular components of the mammalian circadian clock. Hum Mol Genet. 2006;15:R271–7. doi:10.1093/hmg/ddl207.View ArticlePubMedGoogle Scholar
- Lesk V, Taubert J, Rawlings C, et al. WIBL: Workbench for Integrative Biological Learning. J Integr Bioinforma. 2011;8:156. doi:10.2390/biecoll-jib-2011-156.Google Scholar
- Sternberg MJE, Tamaddoni-Nezhad A, Lesk VI, et al. Gene function hypotheses for the Campylobacter jejuni glycome generated by a logic-based approach. J Mol Biol. 2013;425:186–97. doi:10.1016/j.jmb.2012.10.014.View ArticlePubMedPubMed CentralGoogle Scholar
- Côté RA. Architecture of SNOMED. In: Orthner HF, Blum BI, editors. Implement. Health Care Inf. Syst. New York: Springer; 1989. p. 167–79.View ArticleGoogle Scholar
- Rogers FB. Communications to the Editor. Bull Med Libr Assoc. 1963;51:114–6.PubMedPubMed CentralGoogle Scholar
- Ashburner M, Ball CA, Blake JA, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9. doi:10.1038/75556.View ArticlePubMedPubMed CentralGoogle Scholar
- Brinkman R, Courtot M, Derom D, et al. Modeling biomedical experimental processes with OBI. J Biomed Semant. 2010;1:S7+. doi:10.1186/2041-1480-1-s1-s7.View ArticleGoogle Scholar
- Smith B, Ashburner M, Rosse C, et al. The OBO Foundry: coordinated evolution of ontologies to support biomedical data integration. Nat Biotechnol. 2007;25:1251–5. doi:10.1038/nbt1346.View ArticlePubMedPubMed CentralGoogle Scholar
- Noy NF, Shah NH, Whetzel PL, et al. (2009) BioPortal: ontologies and integrated data resources at the click of a mouse. Nucleic Acids Res gkp440. doi: 10.1093/nar/gkp440.
- Lassila O, Swick RR, Wide W, Consortium W. Resource Description Framework (RDF) model and syntax specification. 1998.Google Scholar