METHODS article

Front. Bioinform., 22 June 2023

Sec. Network Bioinformatics

Volume 3 - 2023 | https://doi.org/10.3389/fbinf.2023.1197310

A guide for developing comprehensive systems biology maps of disease mechanisms: planning, construction and maintenance

  • 1. Luxembourg Centre for Systems Biomedicine (LCSB), University of Luxembourg, Esch-sur-Alzette, Luxembourg

  • 2. Independent Researcher, Massy, France

  • 3. Université Paris-Saclay, Laboratoire Européen de Recherche Pour la Polyarthrite Rhumatoïde–Genhotel, University Evry, Evry, France

  • 4. Lifeware Group, Inria Saclay-Ile de France, Palaiseau, France

  • 5. ITTM Information Technology for Translational Medicine, Esch-sur-Alzette, Luxemburg

  • 6. Computer Engineering Department, Bilkent University, Ankara, Türkiye

  • 7. ELIXIR Luxembourg, Belvaux, Luxembourg

Abstract

As a conceptual model of disease mechanisms, a disease map integrates available knowledge and is applied for data interpretation, predictions and hypothesis generation. It is possible to model disease mechanisms on different levels of granularity and adjust the approach to the goals of a particular project. This rich environment together with requirements for high-quality network reconstruction makes it challenging for new curators and groups to be quickly introduced to the development methods. In this review, we offer a step-by-step guide for developing a disease map within its mainstream pipeline that involves using the CellDesigner tool for creating and editing diagrams and the MINERVA Platform for online visualisation and exploration. We also describe how the Neo4j graph database environment can be used for managing and querying efficiently such a resource. For assessing the interoperability and reproducibility we apply FAIR principles.

1 Introduction

Disease maps have been introduced in systems biomedicine as conceptual models of disease mechanisms (; ). In a disease map, interconnected pathway diagrams depict events from the level of molecular processes to higher levels like communication between cells and organs (; ; ). Approaches to the construction of disease maps advanced during past years, supported by the activity of the Disease Maps community (https://disease-maps.io) and formalised into best practices (). However, there is still a lack of clear, formal and well-structured guidelines for main steps of developing a disease map in collaborative projects.

The need for such guidelines became evident during the development of the COVID-19 Disease Map (), a large-scale project supported by an extensive community with more than 200 volunteer participants. The efforts of COVID-19 Disease Map Community would be greatly advanced with step-by-step guidelines to accelerate learning about the approach and creating high-quality resources.

There are a number of publications that can be used as a reference when constructing disease maps (; ). Moreover, brief and informative guideline articles about biocuration () and pathway construction () offer valuable introductory material for constructing disease-oriented resources. In particular, an article on construction of a Reactome-based specific knowledge repository () proposes instructions on reuse of pathway diagrams in a specific context. However, this subject still lacks a comprehensive summary, encompassing the entire lifecycle of a disease map resource.

In this article, we offer guidelines for building diagrammatic systems biology knowledge repositories of disease mechanisms that are explorable online and transformable to computational models. In particular, we discuss planning and design of a disease map including best practices for defining the scope of the resource and communication with domain experts. We discuss the process of biocuration, the required granularity of the map and the required standard languages. Finally, we outline aspects of visual exploration and possible downstream bioinformatic analyses, and the strategy for maintaining the resource following FAIR principles (). We review the literature and examples for a particular area of disease map construction, structuring them into steps of a reference workflow. The methods described in this review are built on the experience of the previously developed disease maps including recently published resources (; ; ; ). Definitions of domain-specific terms are available at https://disease-maps.io/glossary.

2 Disease map development workflow

The major phases of map development are illustrated in Figure 1. The planning phase includes i) determining map purpose and type, ii) defining the coverage of biological functions, iii) choosing a graphical standard (data model) and an editor, and iv) designing map architecture (Figure 1, Steps 1–4). While these steps are shown in a sequence, they are often realised in parallel. Next, in the curation phase, available knowledge from literature and databases is encoded into a diagram with accompanying annotation of map entities (Figure 1, Steps 5 and 6). When the map is ready for publication, it is made available online in a navigable tool and updated based on feedback from domain experts—the publishing phase (Figure 1, Steps 7 and 8). The application phase depends on the goals of the project, expertise and resources available, and can include data visualisation, drug target exploration, computational modelling and hypothesis generation (Figure 1, Step 9).

FIGURE 1

()]. The workflow includes nine steps in four phases (please see in the text). Curators and domain experts are involved in planning and designing the map, then curators develop the map using a diagram editor of choice, advised by domain experts. After publishing online, data analysts, modellers and other users have access to the map for exploration and applications in translational projects, as well as domain experts suggest further evolving the map.

2.1 Determining map purpose and type

At the initial step, it is important to define goals of the project and to determine the map type and its applications. Later, it will influence map’s format, design, size and coverage. The key questions at this stage are:

  • – What is the main focus? Is it one disease or multiple related disorders?

  • – What is the intended use of the map?

  • – What resources are available or needed?

  • – What is the timeline?

  • – Which domain experts (clinicians and biologists) can advise and revise the curation?

Examples of motivation, objectives and primary users for several disease map projects are shown in Table 1. Often they are aligned with the focus of a research centre so the developed mam supports related research. It also can be part of a large-scale translational project for integrating knowledge and future applications in interpreting data generated by the project.

TABLE 1

Primary purposeDetailsExamplesPrimary users
Integrating knowledgeVisualisation, identifying gaps, educational purposesAtlas of Cancer Signalling Networks , Parkinson’s Disease Map , RA-Map Research centres, students, modellers
Combining knowledge with external databases and toolsExploring drug targets, text miningCOVID-19 Disease Map COVID-19 research community, data analysts
Data analysisNetwork-based analysis, omics and multi-omics data visualisation, interpretation and hypothesis generationAsthmaMap IMI U-BIOPRED (No 115010, https://www.imi.europa.eu) community, IMI eTRIKS (No 115446, https://www.etriks.org) community, data analysts
Computational modellingUsing maps as templates for building computational models, refining hypotheses and making predictionsAtherosclerosis Map Research centres, modellers
Computational modelling combined with data analysisUsing maps as templates for building Boolean and hybrid modelling and data analysisRA-Atlas , ; Research centres, modellers, data analysts

Examples of disease map purpose, intended use and expected primary users.

Intended map purpose might include one or several of the following objectives:

  • – integrating available knowledge:

    • – reviewing and visualising disease mechanisms,

    • – identifying gaps in knowledge,

    • – educational purposes;

  • – combining knowledge with external databases and tools:

    • – exploring drug targets,

    • – text mining;

  • – analysis:

    • – network-based analysis,

    • – omics and multi-omics data visualisation, interpretation and hypothesis generation;

  • – computational modelling:

    • – using maps as the basis/template for different types of models,

    • – refining hypothesis,

    • – making predictions.

Clarifying the motivation, objectives and intended applications is needed for the following steps about the coverage of biological functions. This can be documented in a brief project description document (https://disease-maps.io/template), or in a web page (e.g., https://disease-maps.io/projects), or in a project management platform such as FAIRDOMHub (https://fairdomhub.org).

2.2 Defining coverage of biological functions

Exchange with the disease domain experts allows to define the scope and boundaries of the map. Also, key reviews and figures that can help to create an initial skeleton of the map. The main tasks at this stage include:

  • – Identify priority pathways;

  • – Identify cell types, tissues and organs relevant to the pathology mechanisms;

  • – List the molecules known to be involved and drugs used in a given condition

  • – Define disease outcomes on the molecular, cellular and tissue levels, including subdivision into disease stages and subtypes, if applicable.

Such lists can be stored in a document shared and discussed with the domain experts, and maintained throughout the map development, e.g., https://disease-maps.io/template.

Often this information is available in relevant databases, like disease-focused pathways of Reactome (Reactome:R-HSA-1643685) including cancers (Reactome:R-HSA-5663202) and infectious diseases (Reactome:R-HSA-5663205, Reactome:R-HSA-9694516). Disease pathways of KEGG (https://www.genome.jp/kegg/disease), include Parkinson’s disease (KEGG:hsa05012) and COVID-19 (KEGG:hsa05171). The Open Targets Platform (https://platform.opentargets.org) [18] lists molecules relevant for various diseases, e.g., targets associated with allergic asthma (https://platform.opentargets.org/disease/MONDO_0004784/associations).

Outcomes of this step guide the curation efforts in steps 5 and 6 (Figure 1), and the coverage of the biological functions may change depending on available knowledge. In effect, this step is often revisited, so it is important to maintain versioning of the map planning document (https://disease-maps.io/template).

2.3 Choosing a graphical standard and an editor

Majority of disease maps are created in Systems Biology Graphical Notation (SBGN) () using the CellDesigner editor (https://www.celldesigner.org), using the SBGN Process Description (PD) and Activity Flow (AF) languages () (Figure 2). The Entity Relationship (ER) language () is not currently employed in developing disease maps.

FIGURE 2

Process Description (PD) () and Activity Flow (AF) () are the two most used languages of the Systems Biology Graphical Notation (SBGN) standard (www.sbgn.org) () (Figure 2). More information about the SBGN standard are available in the “SBGN Learning” materials at https://sbgn.github.io/learning and the “Figure to SBGN” page at https://sbgn.github.io/figuretosbgn.

CellDesigner (https://www.celldesigner.org) is a diagram editor that is frequently used for developing disease maps. It supports extensive diagrams and follows SBGN logic and makes it possible to draw SBGN-compatible diagrams in Process Description and Activity Flow type of languages using the default CellDesigner’s palette or its Reduced Notation palette. Alternative solutions include the open-source web-based Newt Editor (), the open-source VANTED editor with its SBGN-ED add-on () and the freely available yEd Graph Editor from yWorks (https://www.yworks.com/products/yed). Each tool has its own advantages. Newt offers support for the latest 0.3 SBGN-ML (), SBGN bricks (), advanced automatic layout (), semantic validation of maps, conversion to and from CellDesigner () and SBML formats, and an easy switch to alternative colour schemes. SBGN-ED allows map validation (), conversion from KEGG and SBML formats and PD-to-AF conversion (). yEd has an SBGN palette and SBGN layout tool (), offers intuitive drawing and can support very extensive graphs. The AsthmaMap and the Cystic Fibrosis Map () were developed in the yEd Graph Editor, then converted to the SBGN-ML format via the ySBGN tool () and then to CellDesigner and MINERVA.

The choice of an editor is also connected to the visualisation platform used (Figure 1, step 7). The MINERVA and NaviCell (; ) platforms by design are compatible to CellDesigner’s XML. The MINERVA Platform also accepts SBGN-ML files and SBML files with layout (), and can therefore be used with Newt, SBGN-ED, as well as with yEd in combination with the ySBGN converter. Importantly, the platform also supports export to these formats, allowing to review, improve and reupload disease map diagrams.

2.4 Designing map architecture

Depending on the goals of the project, map coverage and the resources available, there are several possibilities for designing map architecture and for step-by-step progress including starting from PD or AF layer and working with a single diagram or a set of sub-pathways (Table 2).

TABLE 2

#Disease mapArchitectureFormatTop-level view
aAlzPathway , , RA-Map Single diagramProcess DescriptionNone
bFluMap Single diagramProcess DescriptionNone
cCystic Fibrosis Map Single diagramActivity FlowNone
dParkinson’s Disease Map ; ACSN HierarchicalProcess Description, Activity FlowIllustrative
eACSN ; RA-Atlas , COVID-19 Disease Map HierarchicalProcess DescriptionIntegrating diagrama
fAsthmaMap HierarchicalProcess Description, Activity FlowIntegrating diagrama

Examples of map architecture from a single diagram to a hierarchically-organised multi-layered structure.

a

The integrating diagram has elements that represent lower-layer sub-pathways and this way unites them into a single virtual map.

Suggested starting design is a PD diagram of mechanisms with accompanying top-level view of disease mechanisms. This design focuses on the prioritised pathways, allows reuse of the PD layer from other projects, and helps communicating with domain experts via the top-level view (Figure 3). Please note that each map version (A, B, C and D) is a functional, complete resource. Potentially the AF layer can be created automatically from the PD layer using CaSQ () or a PD-to-AF converter (). More options for building a multi-layered structure and example architectures of the published maps are available at https://disease-maps.io/architecture.

FIGURE 3

2.4.1 Modularisation

Modular design of a disease map means developing diagrams that are manageable and reusable () for faster curation, reuse in other projects and downstream use in computational modelling (; ; ). A module is a diagram part that is self-contained and can be developed and maintained independently while being included in a larger network. Modularisation can be applied at different scales including signalling subnetworks, e.g., the MAPK cascade, and on a larger scale, e.g., cell-specific pathways in mast cells, eosinophils and macrophages in asthma (AsthmaMap, https://asthma-map.org/mr).

2.4.2 Examples

AlzPathway (http://alzpathway.org) (; ) catalogues signalling pathways in Alzheimer’s disease as the first comprehensive map of a particular disease. It is a single PD diagram in CellDesigner with an overview image highlighting the included pathways.

The Parkinson’s Disease Map () was originally designed as a single diagram integrating disease-relevant pathways, e.g., mitochondrial dysfunction, calcium homeostasis, synaptic pathobiology, α-synuclein misfolding, failure of protein degradation systems, apoptosis and neuroinflammation. Later a top-level view was added in MINERVA for easier overview and navigation (“SHOW OVERVIEW” button), and a number of pathways were added as underlying PD diagrams (https://pdmap.uni.lu/minerva).

The FluMap () is a comprehensive map of the influenza A virus replication cycle. It is available as a single PD diagram with a simplified version in an AF-compatible format developed in the Reduced Notation in CellDesigner.

The Atlas of Cancer Signalling Networks (ACSN) () is a collection of large-scale PD diagrams that cover mechanisms of cancer and include cell cycle and DNA repair, cell survival, regulated cell death, telomere maintenance, epithelial-mesenchymal transition (EMT) and senescence, invasion and motility, and many more (https://acsn.curie.fr/ACSN2/maps.html). These modules are also available as part of a single extensive map (https://acsn.curie.fr).

The AsthmaMap includes three layers: the detailed Biochemical Mechanisms layer in PD (https://asthma-map.org/bm), the intermediate Molecular Relations layer in AF (https://asthma-map.org/mr) and the top-layer Cellular Interactions view in AF-compatible format (https://asthma-map.org/ci) (). The three layers are interconnected and presented in a hierarchically-organised map in MINERVA (https://asthma.uni.lu).

The Rheumatoid Arthritis Map (RA-Map) was initially designed as a single Process Description map () and then evolved into the Rheumatoid Arthritis Atlas (RA-Atlas), a collection of interconnected diagrams with an access point via a top-level view (). The RA-Atlas includes an updated version of the global RA-Map () covering relevant metabolic pathways and cell-specific molecular interaction maps for CD4+ Th1 cells, fibroblasts, and M1 and M2 macrophages.

The COVID-19 Disease Map is a collection of interconnected Process Description diagrams developed by a community-driven effort in the answer to the SARS-CoV-2 pandemic (https://covid.pages.uni.lu). The map gathers diagrams developed by curators using CellDesigner and VANTED software, describing mechanisms of virus cycle in the host cells, its interactions with host molecular pathways, and the mechanisms of host response to the infection. These diagrams are organised into one hierarchical map via a top-level view diagram (https://covid19map.elixir-luxembourg.org/minerva). Notably, the COVID-19 Disease Map Community involves diagram curators of the Reactome and WikiPathways platforms, hosting their SARS-CoV-2-related diagrams.

2.4.3 Choice and implementation of a map architecture

The choice of a map design should match the biological coverage of the map (

Figure 1

, step 2) and the overall progress of the project. The design should be adapted in case of a change in the coverage, timeline or available resources, as the stages introduced in

Figure 3

are flexible. Considering the discussed examples, a multi-layered and modularised hierarchically-organised map should be considered as a default design (

Figure 3D

). Depending on a project, the focus can be only on one layer, PD or AF, and the top-level view is optional (

Figures 3A, B

). Though, with adequate resources available, it is recommended to work on a multi-layered hierarchical map as the state-of-the-art practice. The main advantages of this hierarchically-organised structure (

Table 3

):

  • – rich content and detailed representation supports different types of visualisation and different types of computational modelling;

  • – a possibility of reusing available resources;

  • – compatibility with similar disease map projects, map parts are easy to reuse in new projects;

  • – a top-level view of disease hallmarks helps communicating with domain experts.

TABLE 3

FeaturesTop layer: top-level viewIntermediate layer: AFDetailed layer: PD
SBGN languageActivity Flow or free formatbActivity FlowProcess Description
ConnectivityInteractionsInteractionsProcesses
GranularityLeast detailed: focus on key molecules and cell typesIntermediate level of detailsMost detailed, the level of processes
Compatible databasesNot availableaKEGG signalling pathways, BioCarta, OmniPath, SIGNORReactome, PANTHER, Recon X, KEGG metabolic pathways
Application/UseCommunication with domain experts for defining disease hallmarksomics visualisation, network analysis, Boolean modelling, integration with machine learning inferred networksomics visualisation, network analysis, ODE modelling, Boolean modelling (via CaSQ)

Comparison of the top-level view vs. Activity Flow vs Process Description layers in a hierarchically-organised map. Resources: KEGG: Kyoto Encyclopaedia of Genes and Genomes (KEGG) (), OmniPath (), SIGnaling Network Open Resource (SIGNOR) (), Reactome (), Recon and ReconMap (; ), PANTHER ().

a

The closest similar visualisations are figures with an overview of disease mechanisms in review papers, representations with cell types, tissues and organs shown.

b

By free-format here we mean flexible representation with one necessary function: all lower-layer sub-pathways in a hierarchically-organised structure are represented by an element on the top-level view and this way integrate them into a single virtual map.

Finally, other design alternatives are possible, e.g., with more high-level layers or using other standard languages. For instance, the disease map approach supports the concept of adverse outcome pathways (AOPs) in the field of immunotoxicology (paper in preparation).

2.5 Formalising knowledge from literature and databases

In the knowledge formalisation step, relevant causal interactions and processes are taken from selected sources including disease-relevant scientific papers and existing representations in pathway databases. In case of scientific papers, causal interactions and processes are manually identified, integrated and encoded into a diagrammatic representation (). In case of pathway databases, relevant causal interactions are manually reviewed and contextualised for a specific disease. Knowledge about molecular processes relevant in the context of a particular disease need to be represented in a standard way, for example, receptor activation, phosphorylation, complex formation, positive and negative regulation of transcription, transport, glycosylation, and more. Knowledge is transformed from textual descriptions and figures in scientific articles into a standardised diagrammatic representation. This representation is also machine-readable and can be further converted into executable computational models.

Knowledge about disease mechanisms is formalised in the SBGN standard within two conceptually different languages, AF and PD. AF is focusing on molecular activities and their sequential connectivity, and PD is focusing on the detailed description of the involved molecular processes (see Table 3; Figure 2). AF is compact and is suitable for overview diagrams and for highlighting key mechanisms. PD is detailed and is better for comprehensive representation of molecular events.

In the curation cycle (Figure 4) information is manually integrated from different sources such as review papers and discussions with domain experts, producing lists of relevant molecules, pathways and cell types. Prioritised pathways are curated in CellDesigner from literature or pathway databases, and visualised in MINERVA for use and review, leading to updates in CellDesigner.

FIGURE 4

), Recon (; ) and PANTHER (), and annotation databases such as UniProt (https://www.uniprot.org) and ChEBI (https://www.ebi.ac.uk/chebi).

The knowledge formalisation and integration works in a similar way to how we build a review paper but, in this case, it is structured and systematic to follow the rules of the graphical standard and to ensure the connectivity in the network. Disease map curation can be seen as a graphical review of relevant disease mechanisms and can be compared to a combination of scoping, rapid and systematic reviews (). A scoping review, corresponding to step 2 of the map development process (Figure 1), determines the coverage of biological functions, defines disease hallmarks, and identifies knowledge gaps. Rapid and systematic reviews, corresponding to steps 5 and 6 (Figure 1), narrow down the focus and collect evidence from literature and databases for drawing and annotating diagrams in an editor. These steps are usually performed by a dedicated curator with advice and feedback from a panel of domain experts.

This “graphical review” requires following diagram-specific knowledge representation formats. It demands consistency in using a standardised graphical language and describing connectivity between diagram elements that may span beyond the currently curated article. For example, if a paper states that IL6 activates epithelial cells, it is important to identify the relevant protein (IL6, UniProt:P05231), identify receptors involved (IL6RA, IL6ST, IL6RB), and describe the following signalling pathway, which may require looking for additional articles.

More details on knowledge formalisation for disease maps are available at https://disease-maps.io/formalisation: the webpage provides examples of text-to-SBGN and figure-to-SBGN transformation, examples of search queries in PubMed and PubMed Central, and describes text-mining capabilities of Europe PMC ().

A map project can include multiple collaboration and exploration pipelines. An extended example of map development and application ecosystem is available for the COVID-19 Disease Map community in Figure 1 in the 2021 publication ().

2.6 Annotating map entities and interactions

Two objectives during the annotation step are identifying map components and providing evidence for interactions. Identification assigns stable identifiers or ontology terms to map components to uniquely define them. This ensures interoperability with systems biology tools and external databases. Providing evidence for interactions defines the context of a given interaction, e.g., a specific cell type, tissue, organism, or relevant pathology. Contrary to components, these annotations are not unique but interaction-specific, and require evidence from a relevant source, e.g., scientific publications, curated pathways, or interaction databases.

2.6.1 Annotating entities

Different types of map entities are shown in Table 4 together with recommended annotation.

TABLE 4

Biochemical entityNamingaIdentifier
ProteinHGNC official symbolbUniProt/HGNC
RNAHGNC official symbolHGNC
GeneHGNC official symbolHGNC
MetaboliteChEBI/PubChem recommended nameChEBI/PubChem
DrugChEBI/PubChem recommended namecChEBI/PubChem/ChEMBL/DrugBank
ComplexSpecific name from literature or listing complex components: Element A:Element BNot required. If available: Complex Portal
CompartmentAppropriate term from Cell Ontology, BRENDA, Cellosaurus, or a specific name from literatureCell Ontology/BRENDA/Cellosaurus
Phenotype (biological processes)Appropriate Gene Ontology (GO) Biological Process (BP) term if availabledGO Biological Process
Phenotype (symptom, disease state)Appropriate MeSH term if availabledMeSH

Recommended annotation of map entities. Resources: HUGO Gene Nomenclature Committee (HGNC, https://www.genenames.org), UniProt (https://www.uniprot.org), Complex Portal (https://www.ebi.ac.uk/complexportal) (), ChEBI (https://www.ebi.ac.uk/chebi), PubChem (https://pubchem.ncbi.nlm.nih.gov), ChEMBL (https://www.ebi.ac.uk/chembl), DrugBank (https://go.drugbank.com), Cell Ontology (https://www.ebi.ac.uk/ols/ontologies/cl), BRENDA (https://www.brenda-enzymes.org), Cellosaurus (https://www.cellosaurus.org), GO (http://geneontology.org), MeSH (https://www.ncbi.nlm.nih.gov/mesh).

a

Flexible if an identifier is provided.

b

This way the protein will be annotated on upload to MINERVA but UniProt will be assigned only if there are no multiple proteins that correspond to that name.

c

Non-standard names are allowed as long as the annotation is provided: for example, “ERK1” instead of “MAPK3” (UniProt:P27361) or “PGE2” instead of “prostaglandin E2” (ChEBI:15551).

d

As in the literature or as suggested by the domain experts.

Proteins, RNAs and genes should be named according to the official symbols of the HUGO Gene Nomenclature Committee (HGNC, https://www.genenames.org). For example, while COX1 is a commonly used name, PTGS1 symbol is recommended. This provides consistency in naming and allows skipping manual addition of other stable identifiers in CellDesigner, because HGNC symbols are recognised by the MINERVA annotation tool on upload.

Annotation of complexes is optional because they can be identified by their components. However, it is recommended to name them by listing the molecules involved, for example, “FCER1A: FCER1G:MS4A2”, or using specific names if available. When annotating, complex-specific annotations are recommended, e.g., GO Cellular Component (CC) or Complex Portal identifiers () (https://www.ebi.ac.uk/complexportal). For instance, GO CC term “Fc-epsilon receptor I complex” (GO:0032998) for the complex above.

Phenotypes or other elements representing higher order processes should be annotated with Gene Ontology (GO) Biological Process (BP) terms () () or MeSH terms (https://www.ncbi.nlm.nih.gov/mesh). GO BP is recommended for annotating components describing biological processes and MeSH terms are recommended for diseases-related terms.

Metabolites (“simple chemical” in SBGN, or “simple molecule” in CellDesigner) are identified via the Chemical Entities of Biological Interest database (ChEBI, https://www.ebi.ac.uk/chebi) (). Alternatively, PubChem (https://pubchem.ncbi.nlm.nih.gov) (), Human Metabolome Database (https://hmdb.ca) () or LIPID MAPS (https://www.lipidmaps.org) () can be used. Automatic MINERVA annotation works for standard ChEBI names, for example, “arachidonic acid” (CHEBI:15843), however in other cases annotations need to be provided manually.

An extended description of annotation for disease maps at offers suggested annotation of diagram information and specific details of annotation fields in CellDesigner.

Other useful sources on the topic of annotation are the recommendations of the logical modelling community (), the Reactome data model and pathway development practices (), the rules suggested for creating reusable pathway models (), the documentation of the MI2CAST format (Minimum Information about a Molecular Interaction CAusal STatement) () and the protocol for high-quality metabolic reconstruction ().

2.6.2 Annotating interactions

Scientific articles contain evidence that confirm the existence of molecular interactions in the context of a specific disease. Each interaction should be annotated with stable identifiers of relevant articles, a PubMed ID or a DOI. If pathway interactions are used, their identifiers can be used instead. Although some guidelines suggest a minimal number of references confirming an interaction (), this should be adapted to the literature coverage of a given disorder.

Importantly, often the relevance of molecular processes in the context of a disease can only be confirmed with a combination of references. For example, one paper can confirm that a specific process is involved in a disease, and another paper will describe details of this process. Only together they offer strong evidence for the interaction.

Some molecular processes related to the disease of interest can be detected in non-human organisms, typically in well-established rodent models for that disease. For parts built based on evidence from cell or animal models, it is important to annotate the corresponding interactions. As proteins, RNA and genes are named using HGNC symbols (human gene symbols), it is recommended to annotate interactions with the NCBI taxon ID of the non-human organism in which the molecular process was determined. For example, if the phosphorylation of STAT3 by JAK2 was determined in mice, the NCBI Taxonomy ID “10090” (NCBI:txid10090) annotation should be added to the interaction.

As most of disease-specific evidence comes from the literature, annotating map interactions requires maintaining a large repository of articles. Such a repository should be stored in a reference manager like Zotero (https://www.zotero.org) or Mendeley (https://www.mendeley.com). Articles should be archived while curating and annotating interactions, preferably in a structure mirroring the organisation of the map. Overall, it is important to check that all interactions in the diagram are annotated with evidence. Automated validation is strongly advised, for instance by using the validators available in the MINERVA Platform.

2.7 Visualising and analysing the maps

Systems biology diagrams of disease maps are developed for online interactive exploration and visual analytics to bridge the gap between the domain experts and bioinformatics experts. The maps can also be used as a part of larger bioinformatic workflows. For example, they can be transformed into a graph database for scalable exploration and analysis, for instance via Neo4j (Section 4), or integrated into a shared network repository like the NDEx platform (). Their content can be transformed with the use of the CaSQ tool for Boolean modelling (Section 2.4; Table 3), or for analysis of network topology with the help of tools such as Cytoscape (). However, the reminder of the section focuses on the visual exploration aspect of disease maps, using a dedicated online platform.

The MINERVA Platform (

;

;

) enables easy online access to and visual interactive exploration of disease maps. Examples are:

Video tutorials are available for introducing key MINERVA functionalities with the AsthmaMap as an example: navigating the maps, adding comments, exploring drug targets and visualising data (https://asthma-map.org/tutorials).

Maps can be uploaded to MINERVA in several formats including CellDesigner SBML, SBML + layout + render, SBGN-ML and GPML (). The conversion module within the MINERVA Platform allows export to these formats as well. On upload users can choose from CellDesigner’s view and SBGN view.

Automatic annotation can be selected on upload. As discussed above, with HGNC names used for proteins, mRNAs and genes, the MINERVA Platform automatically adds links to relevant external databases.

Map browsing and zooming works in MINERVA in the Google Maps manner. For multi-layered maps, interconnected diagrams are organised into a single entity via the submaps function and are visible from the top-level diagram in the SUBMAPS tab (https://asthma-map.org/tutorials/#navigating-the-asthmamap). Search functionality queries all connected maps.

Data can be uploaded as a tab-separated file (https://minerva.uni.lu/doc/examples). This way, omics data (for example, differential expressed genes) can be converted into a gradient of colours on the map (https://asthma-map.org/tutorials/#visualising-data). Multiple datasets can be visualised simultaneously and compared, including a multi-omics combination (e.g., proteomics and metabolomics), or time series data with different time points visualised as separate datasets.

Drug target search connects via online queries to DrugBank (https://www.drugbank.ca) () and ChEMBL (https://www.ebi.ac.uk/chembl) (; ) databases and allows highlighting all entities affected by a drug (https://asthma-map.org/tutorials/#exploring-drug-targets).

Moreover, MINERVA API and plugins () allow the implementation of customised visual exploration and analytics pipelines.

An alternative web-based tool is NaviCell (; ) developed primarily to support the ACSN project. It offers functionalities for visualising omics data for performing functional analysis of the maps, and for computing aggregated values for sample groups and their visualisation on the maps (https://navicell.vincent-noel.fr). A video tutorial on ACSN in NaviCell is available at https://acsn.curie.fr/documentation.html#video.

2.8 Collecting feedback and updating the map

When the first complete version of a map is available online, it is advisable to request feedback from domain experts regarding the coverage of disease hallmarks, the depth and quality of pathways included, accuracy and adequacy of the disease mechanisms representation. This input can be collected in the form of 1) free-format notes, 2) suggested publications, or 3) directs comments on the map in MINERVA (https://asthma-map.org/tutorials/#commenting).

It is also suggested to have an automatic check of any inconsistencies in curation and annotation. This can be done using Neo4j-based tools (see the “Disease map management in the Neo4j graph database” section) or the automated verification in the MINERVA Platform (https://minerva.pages.uni.lu/doc/admin_manual/v16.0/index/#configure-automatic-verification).

2.9 Applying the resource

Map applications can play a role in the development process, as analysis and computational modelling may identify missing feedback loops or analytical endpoints, key pathways, or larger gaps in the network. Therefore, considering the use-cases of a map can help with its curation, revision and refinement ().

One of map applications is omics data visualisation for its interpretation, hypothesis generation, and planning validation experiments. Another is exploration of drug targets in MINERVA using the corresponding functionality (https://asthma-map.org/tutorials). Finally, computational modelling allows to transform a static diagrammatic representation into a dynamic view of disease progression to be used for refining hypotheses and for predictions.

Demonstrated applications of disease maps include:

  • Data visualisation:

    • – omics data visualisation (; ; );

    • – analysis of cell-specific mechanisms using single-cell expression data ();

    • – RNA-Seq-based analysis of the activity of transcription factors ();

  • Network-based analysis:

    • – network analysis and representing map-based molecular signatures for sub-groups of patients in cancer (; );

    • – structural network analysis together with omics data to rationalise the synergistic effects of drugs towards the design of complex disease stage-specific druggable interventions in preclinical studies ();

    • – network-based analysis and prediction of epithelial-to-mesenchymal-like transition mechanisms followed by experimental validation with the use of an animal model, a transgenic mice ();

  • Computational modelling:

    • – creating computational models based on disease maps, for example, for rheumatoid arthritis (; ) and atherosclerosis ();

    • – creating causal-interaction networks based on a disease map ();

  • In preclinical studies with follow-up validation experiments:

    • – in preclinical studies in cancer (; ; ) with validation of the proposed hypothesis by follow-up experiments;

  • Fusing disease maps with other solutions:

    • – integrating disease maps with machine learning inferred networks ();

    • – identifying new crosstalks between pathways from combined analysis of interactions and text mining datasets ().

3 Long-term sustainability of disease maps

Developing a disease map requires a significant amount of work thus it is important to ensure its long-term sustainability and use. First, the map should be deployed on a stable IT infrastructure for uninterrupted access to the content. This can happen via institutional resources or via a third party, such as the open-access map hosting service offered by ELIXIR Luxembourg (https://elixir-luxembourg.org/services/catalog/minerva). Second, a universal programmatic interface should be set up allowing reproducible queries to the disease map contents. The performance of the already available MINERVA API can be further improved by deployment of graph databases, supporting large and complex queries. Finally, a map should be set up as a FAIR resource (https://www.go-fair.org/fair-principles) (), ensuring its long-term reuse. Following sections discuss the topics of graph database deployment and map FAIRification.

4 Disease map management in the Neo4j graph database

Maps offer visual representations primarily intended for human comprehension. For their exchange and edition, they are stored in machine-readable formats, e.g., CellDesigner SBML and SBGN-ML. While these formats are well suited for storing and exchanging graphical layout information, they do not offer a convenient solution for the network-based management (integration, querying, exploration) and analysis of collections of maps. The Neo4j graph database was proposed for an efficient management of large-scale highly-interconnected heterogeneous biological data (), including biological pathway representations in systems biomedicine like Reactome (), IntAct () and Recon2Neo4j (; ). Collections of SBGN maps may also be stored and explored in a Neo4j database using the StonPy software (), which builds a Neo4j graph from an SBGN-ML map or a CellDesigner map converted to SBGN-ML using cd2sbgnml (). Data model of StonPy is close to the visual representation of maps and that is easy to learn. In this model, SBGN nodes (e.g., proteins, states variables) are modelled using Neo4j nodes, while SBGN arcs (e.g., catalyses, production arcs) and relationships between concepts are modelled using Neo4j relationships (edges). Annotations are also stored and can be easily queried. Figure 5 shows an example of a map and an excerpt of its corresponding Neo4j graph as built by StonPy.

FIGURE 5

The StonPy workflow was used to perform a comparative analysis of integrated biological pathways in several major resources including the Atlas of Cancer Signalling Networks (ACSN), PANTHER, and Reactome (; ). The cd2sbgnml () and StonPy tools were also used to develop the dedicated C19DM-Neo4j graph database component (https://c19dm-neo4j.covid.pages.uni.lu), which integrates the biological content of the COVID-19 Disease Map diagrams available in MINERVA.

5 FAIR principles for disease maps available in MINERVA

Disease Maps approaches evolve together with other systems biology communities such as COMBINE (), CoLoMoTo (), SysMod (), SBGN (), SBML () and FAIR () to enhance the accessibility and reusability of the resources ().

The FAIR (Findable, Accessible, Interoperable and Reusable) principles (https://www.go-fair.org/fair-principles) establish criteria to assess and improve the reuse of scientific resources including data, metadata and infrastructure.

The underlying MINERVA Platform provides a set of functionalities that addresses several elements of FAIR:

  • – unified API calls to access the standardised content of all diagrams on the platform;

  • – functionality for translating between systems biology standards;

  • – authentication and authorisation procedures for user management;

  • – explicitly mentioning the reuse licence (for example, the COVID-19 Disease Map diagrams are released in MINERVA under the CC-BY 4.0 licence).

More specifically, MINERVA conforms to FAIR expectations in the following ways:

Findability: The content of the disease maps diagrams in MINERVA is annotated with a range of information elements from initial publications, connected to well-established biomedical resources (e.g., HGNC, UniProt, ChEBI, ChEML, PubMed, etc.) using standard identifiers. The maps are openly available online (e.g., at https://covid19map.elixir-luxembourg.org) and indexed in all major search engines.

Accessibility: The content of biological maps in MINERVA is accessible easily online via standard web protocols. Moreover, MINERVA provides uniform API calls to access the content of the diagrams and can provide required authentication (credentials: user ID and password) and authorisation procedures where necessary.

Interoperability: The biological content in MINERVA follows XML-based systems biology standards such as CellDesigner SBML; translation to other systems biology standard formats is also possible via MINERVA. Moreover, the content of maps is rich in metadata through annotations with information from initial publications connected to biomedical resources as recommended in Table 3, and there are ongoing efforts to connect the content to the Bricks Ontology ().

Reusability: The biological content of disease map diagrams in MINERVA is released under a specific data usage licence, e.g., under the CC-BY 4.0 usage licence for COVID-19 Disease Map (). The biological diagrams are primarily edited using CellDesigner, thus complying with the incorporated rules for molecular mechanism representation, and are annotated with various types of information including references to publications (PMIDs) and standard identifiers for molecular entities (UniProt IDs for proteins, ChEBI IDs for metabolites, etc.).

By conforming to most aspects of the FAIR principles, the utility of disease maps available in MINERVA to the systems biomedicine community is greatly enhanced. While there is always room for improvement, and particularly in the domains of versioning and obsoleting policies, the application of FAIR principles throughout the development process of disease maps ensures the long-term relevance of the resource.

6 Conclusion

Knowledge of disease mechanisms is dispersed across many publications and domains. Systematic reconstructions of disease mechanisms with the help from disease domain experts allow building reference resources where fragmented information is integrated into disease maps as conceptual disease models. These maps are applied for interpreting newly generated omics data and for computational modelling in order to study diseases as complex systems. For this reason, disease maps need to be constructed in a responsible way, based on a formal, structured and detailed protocol enriched with illustrative examples. This will ensure high-quality and consistency, as well as enable efficient coordination in large-scale collaborations and help young researchers to be faster introduced to best practices, tools and methods.

Statements

Data availability statement

The original contributions presented in the study are included in the article/supplementary material, further inquiries can be directed to the corresponding authors.

Author contributions

First draft: AM. Complete draft: AM, MO, MA, AR, IB and DW. Conceptualisation: AM, MO and MA. Project coordination: MO. Online materials: AM, MO and MA. Reviewing and editing: all authors. All authors contributed to the article and approved the submitted version.

Conflict of interest

Authors DR and AK were employed by company ITTM Information Technology for Translational Medicine.

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

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

References

  • 1

    AghakhaniS.SolimanS.NiarakisA. (2022). Metabolic reprogramming in rheumatoid arthritis synovial fibroblasts: A Hybrid modeling approach. PLoS Comput. Biol.18 (12), e1010408. 10.1371/journal.pcbi.1010408

  • 2

    AghamiriS. S.SinghV.NaldiA.HelikarT.SolimanS.NiarakisA. (2020). Automated inference of Boolean models from molecular interaction maps using CaSQ. Bioinforma. Oxf. Engl.36, 4473–4482. 10.1093/bioinformatics/btaa484

  • 3

    BalaurI.MazeinA.SaqiM.LysenkoA.RawlingsC. J.AuffrayC. (2016). Recon2Neo4j: Applying graph database technologies for managing comprehensive genome-scale networks. Bioinforma. Oxf. Engl.33, 1096–1098. 10.1093/bioinformatics/btw731

  • 4

    BalaurI.RoyL.MazeinA.KaracaS. G.DogrusozU.BarillotE.et al (2020). cd2sbgnml: bidirectional conversion between CellDesigner and SBGN formats. Bioinforma. Oxf. Engl.36, 4975. 10.1093/bioinformatics/btaa528

  • 5

    BalaurI.RoyL.TouréV.MazeinA.AuffrayC. (2022). GraphML-SBGN bidirectional converter for metabolic networks. J. Integr. Bioinforma.19, 20220030. 10.1515/jib-2022-0030

  • 6

    BalciH.DogrusozU. (2022). fCoSE: A Fast Compound graph layout algorithm with Constraint support. IEEE Trans. Vis. Comput. Graph.28, 4582–4593. 10.1109/TVCG.2021.3095303

  • 7

    BalciH.SiperM. C.SalehN.SafarliI.RoyL.KilicarslanM.et al (2020). Newt: A comprehensive web-based tool for viewing, constructing and analyzing biological maps. Bioinformatics37, 1475–1477. 10.1093/bioinformatics/btaa850

  • 8

    BentoA. P.GaultonA.HerseyA.BellisL. J.ChambersJ.DaviesM.et al (2014). The ChEMBL bioactivity database: An update. Nucleic Acids Res.42, D1083–D1090. 10.1093/nar/gkt1031

  • 9

    BergmannF. T.CzaudernaT.DogrusozU.RougnyA.DrägerA.TouréV.et al (2020). Systems biology graphical notation markup language (SBGNML) version 0.3. J. Integr. Bioinforma.17, 20200016. 10.1515/jib-2020-0016

  • 10

    BonnetE.ViaraE.KupersteinI.CalzoneL.CohenD. P. A.BarillotE.et al (2015). NaviCell Web Service for network-based data visualization. Nucleic Acids Res.43, W560–W565. 10.1093/nar/gkv450

  • 11

    ChanrionM.KupersteinI.BarrièreC.El MarjouF.CohenD.VignjevicD.et al (2014). Concomitant Notch activation and p53 deletion trigger epithelial-to-mesenchymal transition and metastasis in mouse gut. Nat. Commun.5, 5005. 10.1038/ncomms6005

  • 12

    CoolingM. T.NickersonD. P.NielsenP. M. F.HunterP. J. (2016). Modular modelling with Physiome standards. J. Physiol.594, 6817–6831. 10.1113/JP272633

  • 13

    CzaudernaT.KlukasC.SchreiberF. (2010). Editing, validating and translating of SBGN maps. Bioinforma. Oxf. Engl.26, 2340–2341. 10.1093/bioinformatics/btq407

  • 14

    Del ToroN.ShrivastavaA.RagueneauE.MeldalB.CombeC.BarreraE.et al (2022). The IntAct database: Efficient access to fine-grained molecular interaction data. Nucleic Acids Res.50, D648–D653. 10.1093/nar/gkab1006

  • 15

    DrägerA.HelikarT.BarberisM.BirtwistleM.CalzoneL.ChaouiyaC.et al (2021). SysMod: The ISCB community for data-driven computational modelling and multi-scale analysis of biological systems. Bioinforma. Oxf. Engl.37, 3702–3706. 10.1093/bioinformatics/btab229

  • 16

    FabregatA.KorningerF.ViteriG.SidiropoulosK.Marin-GarciaP.PingP.et al (2018). Reactome graph database: Efficient access to complex pathway data. PLoS Comput. Biol.14, e1005968. 10.1371/journal.pcbi.1005968

  • 17

    FergusonC.AraújoD.FaulkL.GouY.HamelersA.HuangZ.et al (2021). Europe PMC in 2020. Nucleic Acids Res.49, D1507–D1514. 10.1093/nar/gkaa994

  • 18

    FujitaK. A.OstaszewskiM.MatsuokaY.GhoshS.GlaabE.TrefoisC.et al (2013). Integrating pathways of Parkinson’s disease in a molecular interaction map. Mol. Neurobiol.49, 88–102. 10.1007/s12035-013-8489-4

  • 19

    GaultonA.BellisL. J.BentoA. P.ChambersJ.DaviesM.HerseyA.et al (2012). ChEMBL: A large-scale bioactivity database for drug discovery. Nucleic Acids Res.40, D1100–D1107. 10.1093/nar/gkr777

  • 20

    GawronP.OstaszewskiM.SatagopamV.GebelS.MazeinA.KuzmaM.et al (2016). MINERVA-a platform for visualization and curation of molecular interaction networks. NPJ Syst. Biol. Appl.2, 16020. 10.1038/npjsba.2016.20

  • 21

    Gene Ontology ConsortiumDouglassE.GoodB. M.UnniD. R.HarrisN. L.MungallC. J.et al (2021). The gene ontology resource: Enriching a GOld mine. Nucleic Acids Res.49, D325–D334. 10.1093/nar/gkaa1113

  • 22

    GillespieM.JassalB.StephanR.MilacicM.RothfelsK.Senff-RibeiroA.et al (2022). The reactome pathway knowledgebase 2022. Nucleic Acids Res.50, D687–D692. 10.1093/nar/gkab1028

  • 23

    GrantM. J.BoothA. (2009). A typology of reviews: an analysis of 14 review types and associated methodologies. Health Inf. Libr. J.26, 91–108. 10.1111/j.1471-1842.2009.00848.x

  • 24

    HanspersK.KutmonM.CoortS. L.DiglesD.DupuisL. J.EhrhartF.et al (2021). Ten simple rules for creating reusable pathway models for computational analysis and visualization. PLoS Comput. Biol.17, e1009226. 10.1371/journal.pcbi.1009226

  • 25

    HastingsJ.OwenG.DekkerA.EnnisM.KaleN.MuthukrishnanV.et al (2016). ChEBI in 2016: Improved services and an expanding collection of metabolites. Nucleic Acids Res.44, D1214–D1219. 10.1093/nar/gkv1031

  • 26

    HokszaD.GawronP.OstaszewskiM.HausenauerJ.SchneiderR. (2019a). Closing the gap between formats for storing layout information in systems biology. Brief. Bioinform21, 1249–1260. 10.1093/bib/bbz067

  • 27

    HokszaD.GawronP.OstaszewskiM.SmulaE.SchneiderR. (2019b). MINERVA API and plugins: Opening molecular network analysis and visualization to the community. Bioinforma. Oxf. Engl.35, 4496–4498. 10.1093/bioinformatics/btz286

  • 28

    HuckaM.NickersonD. P.BaderG. D.BergmannF. T.CooperJ.DemirE.et al (2015). Promoting coordinated development of community-based information standards for modeling in biology: The COMBINE initiative. Front. Bioeng. Biotechnol.3, 19. 10.3389/fbioe.2015.00019

  • 29

    JdeyW.ThierryS.RussoC.DevunF.Al AboM.Noguiez-HellinP.et al (2016). Drug-driven synthetic lethality: Bypassing tumor cell Genetics with a combination of AsiDNA and PARP Inhibitors. Clin. Cancer Res. Off. J. Am. Assoc. Cancer Res.23, 1001–1011. 10.1158/1078-0432.CCR-16-1193

  • 30

    KanehisaM.GotoS. (2000). Kegg: Kyoto encyclopedia of genes and genomes. Nucleic Acids Res.28, 27–30. 10.1093/nar/28.1.27

  • 31

    KeatingS. M.WaltemathD.KönigM.ZhangF.DrägerA.ChaouiyaC.et al (2020). SBML level 3: An extensible format for the exchange and reuse of biological models. Mol. Syst. Biol.16, e9110. 10.15252/msb.20199110

  • 32

    KimS.ChenJ.ChengT.GindulyteA.HeJ.HeS.et al (2021). PubChem in 2021: New data content and improved web interfaces. Nucleic Acids Res.49, D1388–D1395. 10.1093/nar/gkaa971

  • 33

    KondratovaM.SompairacN.BarillotE.ZinovyevA.KupersteinI. (2018). Signalling maps in cancer research: Construction and data analysis. Database J. Biol. Databases Curation2018, bay036. 10.1093/database/bay036

  • 34

    KupersteinI.CohenD. P. A.PookS.ViaraE.CalzoneL.BarillotE.et al (2013). NaviCell: A web-based environment for navigation, curation and maintenance of large molecular interaction maps. BMC Syst. Biol.7, 100. 10.1186/1752-0509-7-100

  • 35

    KupersteinI.BonnetE.NguyenH.-A.CohenD.ViaraE.GriecoL.et al (2015a). Atlas of cancer signalling network: A systems biology resource for integrative analysis of cancer data with Google maps. Oncogenesis4, e160. 10.1038/oncsis.2015.19

  • 36

    KupersteinI.GriecoL.CohenD. P. A.ThieffryD.ZinovyevA.BarillotE. (2015b). The shortest path is not the one you know: Application of biological network resources in precision oncology research. Mutagenesis30, 191–204. 10.1093/mutage/geu078

  • 37

    Le NovèreN.HuckaM.MiH.MoodieS.SchreiberF.SorokinA.et al (2009). The systems biology graphical notation. Nat. Biotechnol.27, 735–741. 10.1038/nbt.1558

  • 38

    Le NovèreN. (2015). Quantitative and logic modelling of molecular and gene networks. Nat. Rev. Genet.16, 146–158. 10.1038/nrg3885

  • 39

    LicataL.Lo SurdoP.IannuccelliM.PalmaA.MicarelliE.PerfettoL.et al (2020). SIGNOR 2.0, the SIGnaling network open resource 2.0: 2019 update. Nucleic Acids Res.48, D504–D510. 10.1093/nar/gkz949

  • 40

    LysenkoA.RoznovăţI. A.SaqiM.MazeinA.RawlingsC. J.AuffrayC. (2016). Representing and querying disease networks using graph databases. BioData Min.9, 23. 10.1186/s13040-016-0102-8

  • 41

    MatsuokaY.MatsumaeH.KatohM.EisfeldA. J.NeumannG.HaseT.et al (2013). A comprehensive map of the influenza A virus replication cycle. BMC Syst. Biol.7, 97. 10.1186/1752-0509-7-97

  • 42

    MazeinA.OstaszewskiM.KupersteinI.WattersonS.Le NovèreN.LefaudeuxD.et al (2018). Systems medicine disease maps: Community-driven comprehensive representation of disease mechanisms. Npj Syst. Biol. Appl.4, 21. 10.1038/s41540-018-0059-y

  • 43

    MazeinA.IvanovaO.BalaurI.OstaszewskiM.BerzhitskayaV.SerebriyskayaT.et al (2021a). AsthmaMap: An interactive knowledge repository for mechanisms of asthma. J. Allergy Clin. Immunol.147, 853–856. 10.1016/j.jaci.2020.11.032

  • 44

    MazeinA.RougnyA.KarrJ. R.Saez-RodriguezJ.OstaszewskiM.SchneiderR. (2021b). Reusability and composability in process description maps: RAS-RAF-MEK-ERK signalling. Brief. Bioinform.22, bbab103. 10.1093/bib/bbab103

  • 45

    MeldalB. H. M.PerfettoL.CombeC.LubianaT.Ferreira CavalcanteJ. V.Bye-A-JeeH.et al (2022). Complex portal 2022: New curation frontiers. Nucleic Acids Res.50, D578–D586. 10.1093/nar/gkab991

  • 46

    MiH.SchreiberF.MoodieS.CzaudernaT.DemirE.HawR.et al (2015). Systems biology graphical notation: Activity Flow language Level 1 version 1.2. J. Integr. Bioinforma.12, 340–381. 10.1515/jib-2015-265

  • 47

    MiH.MuruganujanA.EbertD.HuangX.ThomasP. D. (2019). PANTHER version 14: More genomes, a new PANTHER GO-slim and improvements in enrichment analysis tools. Nucleic Acids Res.47, D419–D426. 10.1093/nar/gky1038

  • 48

    MiagouxQ.SinghV.de MézquitaD.ChaudruV.ElatiM.Petit-TeixeiraE.et al (2021). Inference of an integrative, executable network for rheumatoid arthritis combining data-driven machine learning approaches and a state-of-the-art mechanistic disease map. J. Pers. Med.11, 785. 10.3390/jpm11080785

  • 49

    MizunoS.IijimaR.OgishimaS.KikuchiM.MatsuokaY.GhoshS.et al (2012). AlzPathway: A comprehensive map of signaling pathways of Alzheimer’s disease. BMC Syst. Biol.6, 52. 10.1186/1752-0509-6-52

  • 50

    Monraz GomezL. C.KondratovaM.RavelJ.-M.BarillotE.ZinovyevA.KupersteinI. (2019). Application of Atlas of cancer signalling network in preclinical studies. Brief. Bioinform.20, 701–716. 10.1093/bib/bby031

  • 51

    NiarakisA.KuiperM.OstaszewskiM.Malik SheriffR. S.Casals-CasasC.ThieffryD.et al (2021). Setting the basis of best practices and standards for curation and annotation of logical models in biology-highlights of the [BC]2 2019 CoLoMoTo/SysMod Workshop. Brief. Bioinform22, 1848–1859. 10.1093/bib/bbaa046

  • 52

    NiarakisA.WaltemathD.GlazierJ.SchreiberF.KeatingS. M.NickersonD.et al (2022). Addressing barriers in comprehensiveness, accessibility, reusability, interoperability and reproducibility of computational models in systems biology. Brief. Bioinform.23, bbac212. 10.1093/bib/bbac212

  • 53

    NoronhaA.DaníelsdóttirA. D.GawronP.JóhannssonF.JónsdóttirS.JarlssonS.et al (2017). ReconMap: An interactive visualization of human metabolism. Bioinforma. Oxf. Engl.33, 605–607. 10.1093/bioinformatics/btw667

  • 54

    OgishimaS.MizunoS.KikuchiM.MiyashitaA.KuwanoR.TanakaH.et al (2016). AlzPathway, an updated map of curated signaling pathways: Towards deciphering Alzheimer’s disease pathogenesis. Methods Mol. Biol. Clifton N. J.1303, 423–432. 10.1007/978-1-4939-2627-5_25

  • 55

    OstaszewskiM.GebelS.KupersteinI.MazeinA.ZinovyevA.DogrusozU.et al (2018). Community-driven roadmap for integrated disease maps. Brief. Bioinform.20, 659–670. 10.1093/bib/bby024

  • 56

    OstaszewskiM.NiarakisA.MazeinA.KupersteinI.PhairR.Orta-ResendizA.et al (2021). COVID19 Disease Map, a computational knowledge repository of virus-host interaction mechanisms. Mol. Syst. Biol.17, e10387. 10.15252/msb.202110387

  • 57

    PartonA.McGilliganV.ChemalyM.O’KaneM.WattersonS. (2019). New models of atherosclerosis and multi-drug therapeutic interventions. Bioinforma. Oxf. Engl.35, 2449–2457. 10.1093/bioinformatics/bty980

  • 58

    PereiraC.MazeinA.FarinhaC. M.GrayM. A.KunzelmannK.OstaszewskiM.et al (2021). CyFi-MAP: An interactive pathway-based resource for cystic fibrosis. Sci. Rep.11, 22223. 10.1038/s41598-021-01618-3

  • 59

    PrattD.ChenJ.WelkerD.RivasR.PillichR.RynkovV.et al (2015). NDEx, the network data exchange. Cell Syst.1, 302–305. 10.1016/j.cels.2015.10.001

  • 60

    RavelJ.-M.Monraz GomezL. C.SompairacN.CalzoneL.ZhivotovskyB.KroemerG.et al (2020). Comprehensive map of the regulated cell death signaling network: A powerful analytical tool for studying diseases. Cancers12, 990. 10.3390/cancers12040990

  • 61

    RougnyA.TouréV.MoodieS.BalaurI.CzaudernaT.BorlinghausH.et al (2019). Systems biology graphical notation: Process Description language Level 1 version 2.0. J. Integr. Bioinforma.16, 20190022. 10.1515/jib-2019-0022

  • 62

    RougnyA.TouréV.AlbaneseJ.WaltemathD.ShirshovD.SorokinA.et al (2021). SBGN Bricks Ontology as a tool to describe recurring concepts in molecular networks. Brief. Bioinform22, bbab049. 10.1093/bib/bbab049

  • 63

    RougnyA.BalaurI.LunaA.MazeinA. (2023). StonPy: A tool to parse and query collections of SBGN maps in a graph database. Bioinforma. Oxf. Engl.39, btad100. 10.1093/bioinformatics/btad100

  • 64

    SatagopamV.GuW.EifesS.GawronP.OstaszewskiM.GebelS.et al (2016). Integration and visualization of translational medicine data for better understanding of human diseases. Big Data4, 97–108. 10.1089/big.2015.0057

  • 65

    ShannonP.MarkielA.OzierO.BaligaN. S.WangJ. T.RamageD.et al (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res.13, 2498–2504. 10.1101/gr.1239303

  • 66

    SiebenhallerM.NielsenS. S.McGeeF.BalaurI.AuffrayC.MazeinA. (2018). Human-like layout algorithms for signalling hypergraphs: Outlining requirements. Brief. Bioinform.21, 62–72. 10.1093/bib/bby099

  • 67

    SinghV.KallioliasG. D.OstaszewskiM.VeyssiereM.PilalisE.GawronP.et al (2020). RA-Map: Building a state-of-the-art interactive knowledge base for rheumatoid arthritis. Database J. Biol. Databases Curation2020, baaa017. 10.1093/database/baaa017

  • 68

    SudM.FahyE.CotterD.BrownA.DennisE. A.GlassC. K.et al (2007). Lmsd: LIPID MAPS structure database. Nucleic Acids Res.35, D527–D532. 10.1093/nar/gkl838

  • 69

    TangY. A.PichlerK.FüllgrabeA.LomaxJ.MaloneJ.Munoz-TorresM. C.et al (2019). Ten quick tips for biocuration. PLOS Comput. Biol.15, e1006906. 10.1371/journal.pcbi.1006906

  • 70

    ThieleI.PalssonB. Ø. (2010). A protocol for generating a high-quality genome-scale metabolic reconstruction. Nat. Protoc.5, 93–121. 10.1038/nprot.2009.203

  • 71

    ThieleI.SwainstonN.FlemingR. M. T.HoppeA.SahooS.AurichM. K.et al (2013). A community-driven global reconstruction of human metabolism. Nat. Biotechnol.31, 419–425. 10.1038/nbt.2488

  • 72

    TouréV.FlobakÅ.NiarakisA.VercruysseS.KuiperM. (2021a). The status of causality in biological databases: Data resources and data retrieval possibilities to support logical modeling. Brief. Bioinform.22, bbaa390. 10.1093/bib/bbaa390

  • 73

    TouréV.VercruysseS.AcencioM. L.LoveringR. C.OrchardS.BradleyG.et al (2021b). The minimum information about a molecular interaction CAusal STatement (MI2CAST). Bioinforma. Oxf. Engl.36, 5712–5718. 10.1093/bioinformatics/btaa622

  • 74

    TüreiD.KorcsmárosT.Saez-RodriguezJ. (2016). OmniPath: Guidelines and gateway for literature-curated signaling pathway resources. Nat. Methods13, 966–967. 10.1038/nmeth.4077

  • 75

    VarusaiT. M.JupeS.SevillaC.MatthewsL.GillespieM.SteinL.et al (2020). Using Reactome to build an autophagy mechanism knowledgebase. Autophagy0, 1543–1554. 10.1080/15548627.2020.1761659

  • 76

    VogtT.CzaudernaT.SchreiberF. (2013). Translation of SBGN maps: Process description to activity Flow. BMC Syst. Biol.7, 115. 10.1186/1752-0509-7-115

  • 77

    WilkinsonM. D.DumontierM.AalbersbergI. J. J.AppletonG.AxtonM.BaakA.et al (2016). The FAIR Guiding Principles for scientific data management and stewardship. Sci. Data3, 160018. 10.1038/sdata.2016.18

  • 78

    WishartD. S.TzurD.KnoxC.EisnerR.GuoA. C.YoungN.et al (2007). Hmdb: The human Metabolome database. Nucleic Acids Res.35, D521–D526. 10.1093/nar/gkl923

  • 79

    WishartD. S.FeunangY. D.GuoA. C.LoE. J.MarcuA.GrantJ. R.et al (2018). DrugBank 5.0: A major update to the DrugBank database for 2018. Nucleic Acids Res.46, D1074–D1082. 10.1093/nar/gkx1037

  • 80

    ZerroukN.AghakhaniS.SinghV.AugéF.NiarakisA. (2022). A mechanistic cellular Atlas of the rheumatic Joint. Front. Syst. Biol.2. 10.3389/fsysb.2022.925791

Summary

Keywords

disease mechanisms, curation, pathway biology, systems biology, translational research

Citation

Mazein A, Acencio ML, Balaur I, Rougny A, Welter D, Niarakis A, Ramirez Ardila D, Dogrusoz U, Gawron P, Satagopam V, Gu W, Kremer A, Schneider R and Ostaszewski M (2023) A guide for developing comprehensive systems biology maps of disease mechanisms: planning, construction and maintenance. Front. Bioinform. 3:1197310. doi: 10.3389/fbinf.2023.1197310

Received

30 March 2023

Accepted

09 June 2023

Published

22 June 2023

Volume

3 - 2023

Edited by

George M. Spyrou, The Cyprus Institute of Neurology and Genetics, Cyprus

Reviewed by

Somnath Tagore, Columbia University, United States

Vincenzo Bonnici, University of Parma, Italy

Updates

Copyright

*Correspondence: Alexander Mazein, ; Marek Ostaszewski,

Disclaimer

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.

Outline

Figures

Cite article

Copy to clipboard


Export citation file


Share article

Article metrics