 Review
 Open Access
 Open Peer Review
 Published:
Grasping frequent subgraph mining for bioinformatics applications
BioData Mining volume 11, Article number: 20 (2018)
Abstract
Searching for interesting common subgraphs in graph data is a wellstudied problem in data mining. Subgraph mining techniques focus on the discovery of patterns in graphs that exhibit a specific network structure that is deemed interesting within these data sets. The definition of which subgraphs are interesting and which are not is highly dependent on the application. These techniques have seen numerous applications and are able to tackle a range of biological research questions, spanning from the detection of common substructures in sets of biomolecular compounds, to the discovery of network motifs in largescale molecular interaction networks. Thus far, information about the bioinformatics application of subgraph mining remains scattered over heterogeneous literature. In this review, we provide an introduction to subgraph mining for life scientists. We give an overview of various subgraph mining algorithms from a bioinformatics perspective and present several of their potential biomedical applications.
Introduction
Graphs or networks are ubiquitous data types, pervasive in multiple domains, from social sciences to medicine, biology and chemistry. One common problem in graph theory consists of finding the underlying subgraph patterns in graphs, which are also referred to as network motifs or graphlets. In this review, we present a survey of frequent subgraph mining applications that deal with biomolecular graphs, such as interaction networks, protein graph structures and chemical structures.
By definition, a graph consists of nodes and edges. In the molecular biology context, each node corresponds to a biomolecular entity and each edge denotes a certain association or interaction between such entities. For example, a protein structure can be represented as a single graph, in which each node corresponds to an amino acid and each peptide bond is represented as an edge. We can consider any chemical structure as a graph, for example by representing each of its constituent atoms as a node and each molecular bond as an edge. On a cellular scale, the interactions between all proteins can be represented as a graph. In this case each protein is a node, and each physicochemical interaction is represented by an edge.
Frequent subgraph mining deals with identifying frequently occurring subgraphs in a given graph dataset. A subgraph can be considered frequent if the number of its occurrences within the dataset exceeds a specified threshold s. However, how the occurrences of a subgraph are counted depends on the graph setting. For example, subgraph counting will be different for a single graph or for multiple graphs. As an illustration of counting within a single graph, let us represent a protein built up out of n=7 amino acids as a single graph, in which each atom is a node and each covalent bond is an edge (Fig. 1a). Figure 1a shows the protein structure at three levels of granularity: (i) as a chain of amino acids; (ii) as a chain of amino acid residues; and (iii) at the atomic level. If we look for subgraphs that appear in the atomic representation graph at least n times (so that s≥n), some very common small twonode subgraphs will be found, such as CC or C=O, and the protein backbone, i.e. the chemical structure that all amino acids have in common: NCC =O, which will occur exactly n times (Fig. 1b). Among the subgraphs that occur less than n times (s<n), either individual amino acids will be found or common patterns that occur across amino acids such as CCCCC (Fig. 1c). It should be noted that the nonoverlapping counting scheme was adopted when counting appearances of subgraphs in a large graph, more on which can be found in “In a single graph” section.
Subgraph mining can be applied to various biological data sets and has a wealth of applications, ranging from finding patterns (i.e. network motifs) that explain functional wiring in protein–protein interaction networks to finding shared properties in molecular compounds, relevant for example in the context of drug discovery [1–4]. However, despite its ability to tackle different biological research questions, a straightforward introduction into the bioinformatics applications and context of subgraph mining has thus far been lacking.
With this paper, we aim to cover theoretical concepts and provide a summary of diverse bioinformatics applications of subgraph mining. It must be noted that due to the myriad use cases of subgraph mining this review will not be able to cover the field completely. Instead we provide a general introduction into subgraph mining and its applications without assuming a specific computational background.
Definitions
In this section, terms are defined that will be used throughout the paper. More detailed definitions are given in Additional file 1. For a more exhaustive introduction into the field of subgraph mining and its algorithms, we refer to [5–7].
Graphs and subgraphs
A graph G is defined as a pair G(V,E) consisting of two sets, a set of nodes V and a set of edges E⊆V×V. If the nodes and/or edges of a graph have labels, then such a graph is considered a labeled graph. A graph is deemed directed if every edge in the graph represents an ordered pair of nodes. If there is no edge orientation in a graph, the graph is deemed undirected. Undirected graphs are common for molecular structures as there is typically no specified direction in the chemical bonds. A graph is deemed connected if there is a path along the edges that links each pair of nodes; otherwise it is deemed unconnected. While most complete chemical structures are connected graphs, many protein interaction networks are for example unconnected. If there is a numeric value assigned to each edge in the graph, then the graph is considered to be weighted and the assigned value is called the weight. A graph whose edges have no weight is considered unweighted. Weights can be used to denote the certainty of the edge, as estimated by experimental or computational determination. It can also imply the strength of the interaction (for example affinity in interactions between biomolecules).
We say that graph G_{s} is a subgraph of graph G if the set of all nodes and the set of all edges of graph G_{s} are subsets of the set of all nodes and the set of all edges of graph G, respectively. A subgraph G_{s} is an induced subgraph of a graph G if its set of nodes is a subset of the set of nodes of graph G, V_{s}⊆V, and its set of edges, E_{s}⊆E, consists of all edges that connect nodes in V_{s} in G. In other words, all edges between the selected nodes are preserved. If at least a single edge is different, yet all nodes can be uniquely mapped from the subgraph G_{s} to the graph G, it is a noninduced subgraph.
If there exists a mapping between two graphs such that if two nodes are connected in one graph by an edge, they will be connected in the other graph as well, then such graphs are considered isomorphic. In other words, if two graphs are isomorphic, it means that those graphs are considered to be equal. If there are two graphs of different sizes, it is clear they cannot be graph isomorphic, i.e. they cannot be equal. However, if the smaller graph is completely within the larger graph, then the graphs are subgraph isomorphic.
Frequent subgraph mining
In the most straightforward application of subgraph mining, the goal is to find those frequent subgraphs that occur more often than a given threshold. We can distinguish the case of multiple graphs and single graphs when defining how the occurrences of a subgraph are counted, which is sometimes referred to as the support of the subgraph. The support will precisely be defined in “Interestingness measures” section. While other definitions of interesting subgraphs do exist, the concept of counting the number of subgraph occurrences is often the first and indeed the most common step in almost all subgraph mining procedures.
Subgraph mining algorithms
In this section, we will discuss the algorithmic aspect of subgraph mining. The general procedure of finding subgraphs of interest is shown in Fig. 2.
Before mining for frequent subgraphs the input graph data needs to be properly encoded. Graphs should be represented so that each subgraph has a unique encoding for easier detection of isomorphisms. After the graph(s) have been encoded, the first step towards finding the frequent subgraphs is to generate a set of candidate subgraphs. These candidates are the collection of subgraphs that could be frequent in the graph and need to be checked. After obtaining this set of candidate subgraphs, the next step is to count the number of occurrences of these subgraphs in the graph dataset. Frequency counting is usually the most computationally intensive part of subgraph mining algorithms. For this reason, it is essential to have as few candidate subgraphs as possible. Therefore the counting is commonly preceded by search space pruning and removal of redundant graphs in order to achieve a reduction in the number of candidates. The counting procedure will output how many times each subgraph that has been checked occurs, which can be used to calculate the interestingness of a given subgraph. This can be as straightforward as checking against a fixed threshold. Alternatively, it can include more advanced methods such as calculating statistical enrichment compared to a background distribution.
Subgraph representation
Graphs and subgraphs can be represented in several ways. The two most commonly used graph representations are the adjacency matrix and the adjacency list.
Adjacency matrix The graph is represented as an n×n binary matrix; n being the total number of nodes in the graph. The position (i,j) in the matrix can have the value of 0 or 1 for an unweighted graph, depending on whether or not there is an edge between nodes v_{i} and v_{j},(v_{i},v_{j})∈E (Fig. 3a). For weighted graphs, instead of a binary matrix, entries in the adjacency matrix will contain the weights assigned to the corresponding edges.
Adjacency list The adjacency list of a graph is an array A of length n, where n is the number of nodes in a graph. Each entry A[ i] in the array is linked to a list of all nodes connected to the node v_{i} (Fig. 3b). For directed graphs, each edge (v_{i},v_{j}) is stored exactly once, while for undirected graphs, each edge (v_{i},v_{j}) is stored twice; once in a list connected to the v_{i} node and once in the v_{j} related list.
Canonical labeling Neither the adjacency matrix nor the adjacency list take into account graph isomorphism discovery, which means that it is possible for two graphs which are isomorphic not to share the same adjacency matrix/adjacency list. Canonical labeling ensures that if two graphs are isomorphic they share the same canonical representation. There are several possible ways for canonical labeling. Here, we will present two of the most common approaches.
Canonical adjacency matrix (CAM) The CAM is a unique adjacency matrix [8, 9]. Given an undirected graph G of size n and its adjacency matrix A, the CAM code is obtained by concatenating the entries of the upper triangular part of the adjacency matrix A in a columnwise fashion. The canonical code of a graph G is obtained by comparing all possible codes of the graph and choosing the one with the minimal (maximal) lexicographic value (Fig. 3c).
Depthfirst search (DFS) code The DFS procedure assigns a unique label to each node while traversing the graph in a depthfirst fashion (more about graph traversal strategies in “Subgraph searching and matching” section) [10]. This type of search strategy, specifically the way that subgraphs are extended, results in each subgraph having a unique canonical label.
De novo subgraph candidate generation
There are several methods for generating subgraph candidates de novo, most of which can be possible induced undirected unlabeled candidate generation and extendbased candidate generation.
Joinbased candidate generation With this method, a new subgraph candidate of size (k+1) is created by joining two subgraphs of size k that were discovered to be frequent in a previous iteration. Two subgraphs are joined if and only if they share the same (k−1)–size subgraph. Depending on the definition of the subgraph size, joinbased candidate generation can be either nodebased (Fig. 4a) or edgebased (Fig. 4b). Thus k will either be equal to the number of nodes, or the number of edges. The main issue with joinbased candidate generation is that multiple candidates can be generated in one join process and equivalent candidates can be generated in a number of different ways, which can lead to the occurrence of duplicate candidates.
Extendbased candidate generation A subgraph candidate of size (k+1) is generated by extending a frequent ksize subgraph with an additional node. The issue with this approach is that there are multiple ways to add the node to a subgraph (it can be added to multiple nodes) which leads to the generation of redundant candidates. One way to solve this is to only add extensions that meet specific criteria. The most common approach in this manner is to use the rightmost path extension, i.e. to add the extra node only to a node in the rightmost branch (Fig. 4c) [11, 12].
Subgraph searching and matching
If the goal of subgraph mining is to find all subgraphs that occur more than a set number of times, the subgraph candidates can be pruned based on the Apriori principle. This states that the number of occurrences of a subgraph can never be higher than those of the subgraphs that it contains. Therefore if a subgraph does not pass the frequency threshold, any future candidate graph containing this subgraph need no longer be checked. Two types of search strategies exist that are based on this principle: breadthfirst search (BFS) and depthfirst search (DFS).
The breadthfirst search (BFS) strategy The BFS strategy checks the support of all candidates of a certain size, before moving to the next level; i.e. first all possible candidate subgraphs of size k will be generated and checked for support, subsequently the frequent subgraphs will be retained and used to generate the candidate subgraphs of size (k+1). A BFS is necessary if the subgraph candidates are generated by the joinbased generation method. For example, to generate the (k+1)size subgraph candidates two ksize frequent graphs are needed, which means that all frequent subgraphs of size k need to be determined first (Fig. 5a). This approach accounts for effective candidate pruning, but at the cost of a high memory usage.
The depthfirst search (DFS) strategy The DFS strategy first checks the support of a candidate subgraph of size k; if this subgraph is frequent it will be extended to size (k+1) and checked for support again (Fig. 5b). The subgraph will continue to be extended until it is no longer frequent. Compared to the BFS this approach requires less memory but at the cost of less effective pruning.
There is another generation and search strategy that is common in bioinformatic approaches, where all possible subgraphs are generated in advance and each is subsequently tested in turn. This is only feasible if the candidate space is fairly small so that there are few possibilities, for example only unlabeled subgraphs of a fixed small size, or if there is prior knowledge about what kind of subgraphs need to be tested, such as feedforward loop structures. The most common of these deal with induced subgraphs called graphlets, which are very common for the analysis of interaction networks [13–15]. Graphlets are all possible induced undirected unlabeled connected subgraphs up to a given size, usually 5 nodes so that there are 30 unique graphlets (Fig. 6). In these cases, all subgraphs that need to be checked can be predefined. Therefore, no pruning is needed and the occurrence of every defined subgraph is counted.
Interestingness measures
After obtaining a set of candidate subgraphs, the frequency of these subgraphs needs to be counted. In most cases, frequency counting is the most computationally intensive part of the subgraph mining algorithms.
Frequent subgraph mining. Let D be a graph dataset (a graph database or a single graph) and let σ be a userdefined minimal support threshold. Frequent subgraph mining deals with finding all subgraphs G_{S} in D such that:
In multiple graphs
When searching across multiple graphs, the count of a subgraph is usually defined as the number of graphs that contain it, independent of the number of times that the subgraph actually occurs within each of these graphs (Fig. 7).
Support of a subgraph (in multiple graphs) Let D be a graph database. The support of a subgraph G_{S} is defined as the number of graphs in D that contain G_{S} as a subgraph:
Frequency of a subgraph (in multiple graphs) Given a database of graphs D the frequency of a subgraph G_{S} is defined as the fraction of all graphs in D that contain subgraph G_{S}:
Statistical significance
The most frequent subgraphs are not necessarily the most relevant for a given research problem. The statistical significance of a subgraph is usually determined by checking how often that subgraph appears in a background graph database, for example one where the edges have been randomized within each graph. Statistical significance can then be used as a postprocessing step to filter out insignificant frequent subgraphs or to mine for subgraphs that are significantly associated with a certain subset of the database (i.e. subgroup discovery) [16–18] (Fig. 9).
In a single graph
Support of a subgraph (in a single graph) Given a graph G the support of a subgraph G_{S} is defined as the measure of occurrence of subgraph G_{S} within graph G. For a single graph, the support and frequency are the same. However, counting the support within a single graph is not as straightforward as for multiple graphs [19, 20]. First of all, different occurrences of a subgraph can overlap and this can lead to several difficulties. When dealing with overlapping subgraphs there are two approaches to count the occurrences of a subgraph within a single graph: the first approach takes into account only nonoverlapping subgraphs; the second approach allows an overlap between various subgraphs (Fig. 8a). In the first approach, only those subgraphs which share no edges (or no nodes) are considered nonidentical, while in the second approach two subgraphs are deemed nonidentical if they differ by at least one edge (or one node).
When considering how to deal with overlapping subgraphs within a single graph, one must also take into account the Apriori property. This states that the count of the larger set (graph) cannot be greater than the count of any of its subsets (subgraphs), and allows pruning of a large portion of the search space. However, for overlapping subgraphs simply counting the number of occurrences of a subgraph within a graph might not adhere to this property because parts of overlapping subgraphs are shared. This is shown in the example in Fig. 8b. If g_{1} is a subgraph of g_{2}, it should have an equal or higher count. However, it is clear that g_{2} has a higher count. To this end, other overlap support measures have been proposed, a prominent one being the size of the maximum independent set (MIS) measure of graphs [21–23].
Statistical significance Sometimes the subgraphs of interest are not the ones that occur most frequently, but the ones that are significantly enriched in certain nodes. The significance of a subgraph can be determined based on a statistical background distribution (Fig. 9). Significant subgraphs can then be those that occur more frequently than expected in a random graph [24].
Sampling Sometimes the graph is simply too large to be efficiently analyzed and thus the use of sampling approaches is required. Graph sampling considers only a representative sample of the graph instead of the entire graph. This means that when searching for frequent subgraphs, instead of enumerating the entire graph, only a representative sample is taken into account. More on sampling in large graphs can be found in [25–27].
Biomedical subgraph mining
Subgraph mining methods have been applied to a wide variety of biomedical problems. The reason for this is that many biological features and datasets can be represented as a graph or a collection of graphs, which enables usage of these methods. Often this requires conversion from the typical format into a graph representation, and in some cases this incurs a loss of information. However, in these cases, the advantage of being able to use subgraph mining approaches far outweighs any disadvantages. A broad distinction can be made between the problems that deal with multiple graphs, and those that deal with single graphs. The following is a nonexhaustive description of some of the most common usages of subgraph mining in bioinformatics, together with frequent subgraph mining algorithms and available implementations. A flowchart detailing the differences in various subgraph mining implementations, as featured here, can be found in Fig. 10.
In multiple graphs
Application in bioinformatics
As introduced in “Interestingness measures” section, the goal of subgraph mining in multiple graphs typically concerns finding those subgraphs that occur at least once in a sufficiently large number of these graphs. As such, these applications focus on finding common subgraph patterns that are ubiquitous in the graph collection, which may indicate their biological relevance.
Molecular data
One of the archetypical examples is subgraph mining of a database of (bio)molecules for common subgraph patterns. Indeed such datasets are often used as benchmarks for the development of novel subgraph mining approaches [10, 28–31]. In this instance, each graph in the collection represents a single molecule. The atoms of the molecule are then represented as nodes, and the covalent bonds between the atoms as edges. The edges here are therefore undirected, however they are sometimes labeled to denote the type of bond, or weighted to denote the strength or length of the bond. The goal of the subgraph mining often involves the discovery of common molecular structures, such as benzene rings or disulfide bonds that can be related to a chemical, physical or biological feature. This is for example useful for the discovery of new drug compounds. By mining those molecules that are known to have a specific drug activity for common patterns, one can identify potential new candidates by searching for unstudied molecules that contain the identified subgraphs [32–34].
Protein structures
The threedimensional structure of a protein is often represented as a graph for mining purposes. In these cases, this will usually not be at the level of individual atoms, as described for other biomolecules in the previous section, but at the level of amino acid residues. The graph representation will then typically not concern itself with the chemical bonds that exist to build up the protein structure, but will describe those amino acids that are in close proximity, regardless of whether or not an actual bond exists between them. The nodes of the graph thus represent the different amino acid residues of the protein, and two nodes are connected with an edge if the residues are in close spatial proximity in the threedimensional protein structure. Common criteria that define close proximity include a fixed maximum distance, often 7 or up to 11.5 angstrom, between the αcarbon atom of the amino acids, i.e. the carbon atom that the residue chain is attached to, or another fixed position in the amino acid structure [4, 35, 36]. These distances are defined such that it can be supposed that if two amino acids are sufficiently close they likely interact or can influence each other in some manner. The nodes are then labeled with the type of amino acid, e.g. alanine, leucine and so on. The edges are undirected and usually left unweighted. The mined subgraphs therefore consist of patterns of amino acid residues that occur with high frequency in the collection of proteins. While the goal of this mining can be fairly diverse, it most commonly involves correlating the found subgraphs with protein functions. Indeed, these patterns often represent evolutionary conserved threedimensional structures or domains, which may indicate that they have a biological function in these proteins [35, 37], and can in turn be used as features for functional prediction.
Phylogenetic trees The construction of phylogenetic trees that represent similarity between sequences is a common staple in bioinformatics. As a tree is a specific type of graph, many of the approaches detailed here can be readily applied to collections of trees to find socalled subtrees. In this case, the nodes represent the different genes, proteins, species or other biological entities of interest that were contained in the original dataset for comparison, supplemented with the branch nodes signifying inferred ancestors. In the case of a rooted tree, one branch node will be labeled as the root node. The edges denote the relationship between the studied entities and their ancestors. These edges can be weighted with the branch length, however these are often not used in the mining process. The actual mining proceeds in a manner very similar to that of normal subgraph mining as already outlined in this review. However, trees do have specific properties that can be utilized to speed up the mining process. The mining of phylogenetic subtrees is most relevant in the context of the discovery of robust phylogenetic relationships in a larger collection of phylogenetic trees. For example, many methods used to study the evolutionary relationship between sequences generate a large number of different trees, each representing a potential phylogenetic possibility. A common approach is therefore to search for the consensus tree, i.e. the single tree that contains those phylogenetic relationships that most of the generated trees agree on. By mining the different solutions, those frequent subtrees can be found, and combined into the consensus solution [38–40].
Subgraph mining tools available for bioinformatics use
Here we present a short overview of the algorithms and available implementations to mine subgraphs in multiple graph settings. Table 1 summarizes the information on the algorithms and tools for multiple graphs. A more detailed overview of these frequent subgraph mining algorithms and their implementations can be found in [41].
AGM/AcGM The Apriori graph mining algorithm (AGM) [9] represents graphs by an adjacency matrix and searches for frequent subgraphs in a BFS manner. Candidate subgraphs are created by joinbased candidate generation, which in this case is nodebased (i.e. in each iteration one node is added to the subgraph). AcGM [42] is a version of AGM that only looks for frequent connected subgraphs. The AcGM executable for Linux is available.
FSG FSG [43, 44] represents graphs with adjacency lists and uses a BFS strategy to find frequent subgraphs. As opposed to AGM (and AcGM), FSG generates candidate subgraphs by adding a new edge in each iteration, i.e. it uses edgebased candidate generation. An FSG implementation can be found as a part of the PAFI software, which contains executables for Linux (PAFI is written in ANSI C and C++).
GASTON The Graph/Sequence/Tree extraction (GASTON) [30] algorithm divides the subgraph mining procedure into different subtasks by evaluating subgraphs from a lower to a higher complexity. It starts by looking for paths first, continues with free trees and in the final stage it looks for cyclic graphs. Here, a path is a sequence of nodes connected by edges, which includes every node only once; a free tree is a connected graph which does not contain any cycle; and a cyclic graph is a graph which contains at least one cycle. The source code for the Gaston algoritm is available in C++ and Java (Parallel and Sequential Mining Suite (ParSeMiS), ParMol package).
gSpan gSpan [10] is one of the most well known algorithms for frequent subgraph mining. This algorithm traverses the search space in a DFS manner and creates new subgraphs using the rightmost path extension with minimum DFS code. This way it combines the creation of new subgraph candidates with subgraph isomorphism testing. The source codes of gSpan is available in C++ (gBoost toolbox) and Java (Parallel and Sequential Mining Suite (ParSeMiS), ParMol package).
FFSM The Fast Frequent Subgraph Mining (FFSM) [29] algorithm uses a CAM representation of graphs and traverses the search space in a DFS manner. It stores an embedding set for each discovered frequent subgraph and thus avoids subgraph isomorphism testing. Both binaries and Java source code (part of ParMol package) for FFSM are available.
MoFa The Molecular Substructure Miner (MoSS)/Molecular Fragment Miner (MoFa) [34] is an algorithm inspired by eclat [45], designed for finding frequent fragments (frequent connected subgraphs) in molecular data. It uses adjacency matrices for graph representation and traverses the search space in a DFS manner. It is implemented as part of the ParMol package and in the Molecular Substructure Miner (MoSS). Although it was initially designed for molecular data, MoSS can also work with other types of graphs. It supports input for molecular files in the SMILES, SLN or SDfile formats.
Comparison of tools
To give an overview of the performance of the various subgraph mining algorithms we evaluated each on a protein structure graph dataset derived from all RCSB PDB nonredundant structures, similar to those described in “Application in bioinformatics” section. The graph statistics are summarized in Table 2. We compared the performance for support thresholds of 50, 30, 20 and 10%. Tables 3 and 4 show CPU time in seconds and memory consumption in MB, respectively. There is no information about memory consumption for AcGM and FSG, since they are binary executables. These values are only indicative as they result from a single run on a single computer.
The ParMol package contains implementations of four subgraph mining algorithms and is convenient in case the user wants to test several algorithms on the same dataset. However, MoFa/MoSS is more userfriendly than the ParMol package, and accepts a wide variety of input formats. Testing our example graph database, AcGM had a running time one order of magnitude slower than its competitors. FSG was the second slowest algorithm with a running time one order of magnitude slower than GASTON, gSpan, FFSM and MoFa/MoSS for a support threshold of 20% and lower. For a lower support threshold (10%) GASTON and MoFa/MoSS had a twofold faster running time than gSpan and FFSM.
Interpretation of the results
The output of the tools for mining subgraphs in multiple graph settings is given in Section 10.2. All methods, except for FSG, find the same patterns for all support thresholds. FSG returns a lower number of frequent subgraphs than the remaining algorithms because it does not count one node subgraphs. As the input concerned a large number of protein molecular structures encoded as graphs, the output graph patterns concern configurations of amino acids that are common within the molecular structure. For example, one of the threenode subgraphs that was found with high frequency (50.88%) was LEUGLUARG. The close proximity of glutamate and arginine in this subgraph are indicative of a salt bridge, a common configuration in protein structures. Another example with high frequency (40.23%) was THRVALPHE, which may be indicative of a βsheet as all three amino acids are known to be preferred in the middle of βsheets and thus would be in close proximity.
In a single graph
Application in bioinformatics
Much of today’s biomedical relationship data is represented as a large single graph or network representation. Mining these single graphs involves identifying those subgraphs, substructures, or network motifs that are common, and specifically those that are more common than a given background distribution would suggest.
Interaction and association networks The interactions and associations between the different biomolecules in cells, such as genes, proteins, and metabolites, are most commonly conceptualized as a graph. The biomolecules are featured as the graph nodes. The edges in these graphs can then represent direct interactions, such as protein–protein complexes or kinase–target relationships, but also indirect interactions, such as between genes coding for transcription factor proteins and the gene targets of these transcription factors. For association graphs, the edges are biomolecules that are cooccurring, such as in coexpression networks, or sharing similar functionality. Despite representing different types of relationships, these graphs have very similar properties: they are often very large, featuring several hubs of nodes with a very high degree, and are therefore often described as hairballs, as can be seen in Fig. 11. Depending on the relationship type, the edges are either directed for those with a strict regulatortarget relationship, or undirected for associations and complexes. Most of the studies in this field are done on graphs that have unlabeled nodes and edges. The objective of subgraph mining interaction and association networks often concerns the identification of those subgraphs, termed graphlets or network motifs, that are present in a much higher (or lower) quantity than in a collection of similarlysized random networks that function as the background for the relevant statistical test. Given the size of these networks, typically only induced unlabeled subgraphs are examined. These are then generally described as the building blocks of the network [46, 47], and have been related to specific biological functions [48]. The identified network motifs can then be used in a number of different contexts, such as as a measure for graph similarity [14], network alignment [49, 50], and network prediction [51]. Several reviews are available that focus on the identification and usage of network motifs in these types of graphs [52–55].
Other biomedical graphs
A few other single graph problems are worth mentioning as they demonstrate that subgraph mining in bioinformatics is not limited to molecular data. The first is the mining of brain connectivity maps or neuronal networks [56]. These are graphs that represent the synaptic connections between different neurons. The neurons are thus represented as unlabeled nodes, with directed edges as connections from one neuron to the next. The second involves the mining of food webs that document the relationships between species in an ecosystem [47, 56, 57]. Each node represents a species and directed edges denote predatorprey relationships. In both cases, the goal of the mining is typically the identification of common subgraphs that may represent network building blocks with a specific biological function.
Subgraph mining tools available for bioinformatics use
We present a short overview of the algorithms and tools available for mining subgraphs in a single graph. Summarized information can be found in Table 5. It should be noted that we only considered algorithms for which implementations were available. A more detailed overview of the algorithms and tools for motif detection can be found in [58].
Mfinder
Mfinder is a tool for network motif detection. It contains implementations of algorithms for both exhaustive motif enumeration and motif sampling [46, 56]. Exhaustive motif enumeration looks for all possible subgraphs of a certain size (expressed as a number of edges) that appear more frequently than expected in a random network. Motif sampling for subgraph counting estimates the frequency of subgraphs by sampling within a subset of the whole graph. Random subgraphs are sampled by randomly choosing an edge, and then extending that edge until a subgraph of the required size is obtained. Both algorithms look for subgraphs in a DFS manner and generate candidate subgraphs by adding a new edge in each iteration, i.e. they use edgebased candidate generation. Mfinder is compatible with MDraw, a network visualization tool. Mfinder was the first tool for motif detection and it works with both directed and undirected unlabeled graphs and supports a motif size between 3 and 8. Unfortunately, currently it is no longer under active development.
MAVisto
MAVisto is a Java Web Start application that contains an implementation of the frequent pattern finder (FPF) algorithm [59] for motif detection. This algorithm searches for all motifs of a certain size (which can be defined either as a number of nodes or a number of edges) that appear more frequently than in an ensemble of random networks by using one of three different concepts of motif frequency. These frequency concepts differ in the way they allow overlapping matches. FPF builds a pattern tree with the simple oneedge, twonodegraph at its root and expands this by iteratively adding a new edge. It utilizes the downward closure property and spends time only on the promising branches. The MAVisto tool works with both directed and undirected unlabeled graphs, but it is currently no longer under active development.
FANMOD
FANMOD is a tool that implements the RANDESU algorithm for enumerating and sampling subgraphs [60] as well as the full enumeration algorithm. RANDESU was designed to address the bias of the sampling method for subgraph counting implemented in Mfinder, as Mfinder’s sampling method is prone to sampling certain subgraphs more often than others. RANDESU fixes this bias and is faster. It enumerates all subgraphs of a certain size, although during the execution it will ignore some of these to achieve an unbiased sampling. FANMOD works with both directed and undirected, labeled and unlabeled graphs and it can detect motifs of up to a size of 8 nodes (min size: 3).
Kavosh
The Kavosh algorithm [61] is designed with the goal of finding motifs with less memory and CPU time. It enumerates all subgraphs of a certain size in the graph first by finding all subgraphs that contain a certain node, then removing that node from the graph and repeating the process. Kavosh works with both directed and undirected unlabeled graphs and poses no restriction on the motif size (defined as the number of nodes) to be searched.
NetMODE
NetMODE [62] is a software package for motif detection designed specifically as an improvement over Kavosh (upon which it was built) and FANMOD. As opposed to FANMOD, which contains implementations of both sampling and full enumeration algorithms, NetMODE exclusively performs full enumeration. In contrast to Kavosh it searches only for motifs of size 3–6 nodes. Initially much faster than both Kavosh and FANMOD, NetMODE can also be run in parallel. It works with both directed and undirected unlabeled graphs.
accMotif
The accelerated motif (accmotif) program [63] implements a motif detection algorithm that counts motifs of size k using the set of induced subgraphs of size k−2. Here, the size of the motif is equal to the number of nodes. This algorithm is efficient only for smaller motifs (up to size 6). It allows multithreaded execution and works with both directed and undirected unlabeled graphs.
Significant Subgraph Miner
Significant Subgraph Miner (SSM) [24] is an implementation of the algorithm for finding subgraphs associated with a certain set of nodes in a single graph, as per subgroup discovery. The algorithm generates candidate significant subgraphs in a DFS manner, starting with the node of interest as the root node. If significantly more nodes of interest are root nodes for a subgraph compared to the other nodes of the graph (which are not of interest), the subgraph will be deemed as significantly enriched. SSM works with both directed and undirected, labeled and unlabeled graphs.
Comparison of tools
To compare these tools we used the protein–protein interaction graph of Escherichia coli as present in the Bacteriome database [64] with a confidence cutoff of 0.8. The graph was left unlabeled as only SSM supports labels. The graph statistics are summarized in Table 6 and it is visualized in Fig. 11a. We compared the performance for motif sizes of 3, 4 and 5. The number of randomized networks was 100 for motifs of size 3, 10 for 4size motifs and 5 for 5size motifs. All tools were compared in full enumeration mode and a single thread. Table 7 shows CPU time in seconds, while Table 8 shows memory consumption in MB. These are only indicative values as the measurements were made only once on a single computer.
Mfinder and MAVisto exhibited the longest running times, as was shown in [61]. FANMOD is two orders of magnitude faster than Mfinder and is overall the best tool in terms of userfriendliness. FANMOD supports both a full enumeration and a sampling mode, and works with labeled data as well. As expected Kavosh is slightly faster than FANMOD and consumes less memory. Kavosh also has one big advantage — it does not pose a restriction on the motif size. NetMODE and accMotif have similar performances when it comes to speed, they are both one order of magnitude faster than Kavosh with accMotif being slightly faster. They both have multithreading capability. This allows them to be run on graphs that can be too large for FANMOD to handle. However, as opposed to Kavosh they are restricted to motifs of size up to 6 nodes and do not work with labeled graphs. In terms of memory consumption, accMotif requires more memory. SSM is quite unique since it does not look for motifs of a specific size (it only has a restriction for maximum edges allowed), but for all subgraphs which are associated with certain nodes. However, it can also enumerate all subgraphs within a graph. In terms of running time, it has roughly a running time of the same order of magnitude as Kavosh. Its memory consumption is one order of magnitude larger than accMotif’s.
Interpretation of the results
The output of the tools for mining subgraphs in single graph is given in Additional file 2. As the starting graph was a proteinprotein interaction network, all methods returned subgraphs that represent configurations of proteins that are frequent or significantly enriched compared to a randomized graph of the same size. Two configurations are common throughout all of the results, independent of the subgraph size that was being investigated. The first are densely connected subgraphs where all nodes have edges that connect to (almost) all other nodes. These are indicative of protein complexes, where many proteins interact with each other. These will coelute in pulldown experiments and create dense clusters within the graph. For example, MFINDER subgraph 238, featured in the supplementary materials, is a fully connected graph between three nodes. This threenode graph occurred 4683 times in the proteinprotein interaction network, but only on average 3429 times in random graphs of a similar size. This graph was therefore enriched with a Pvalue less than 1×10^{−32} in the proteinprotein interaction network. The second type of common enriched configuration is where a single protein interacts with all other proteins, but they (hardly) interact with each other. These subgraphs are indicative of hub proteins, a well known feature of proteinprotein interaction networks where a single protein has a large number of interactions. An example is the FANMOD subgraph 4382, which contains a single node that is connected to three other nodes with no other shared edges. However, this subgraph was found to occur less than expected when compared to a random graph, indicating that while hub proteins are present in the graph, they are relatively rare. The differences between the output of the various methods are due to their different search parameters. Mfinder finds fewer motifs than the remaining tools because it calculates motif sizes based on their number of edges instead of their number of nodes. SSM finds more subgraphs because it filters the subgraphs based on a subgroup discovery criterion and to directly compare it to the other approaches the entire graph was added as a single subgroup. The remaining methods found the same number of motifs.
Conclusions
We have given an overview of the basic principles of subgraph mining and how it can be used for biomedical applications. We made a distinction based on those methods and problems that feature finding frequent subgraphs across multiple graphs and those that find them within a single graph. These approaches are very similar algorithmically. The basic steps thus remain: first graph encoding, then candidate generation, and finally subgraph counting. However they use different criteria to define which graphs are actually interesting.
Despite its broad applicability there are still some enduring issues that prevent subgraph mining in a biomedical setting from reaching its full potential. First, labeled single graph solutions are very rare, yet very relevant for biomedical research [65–69]. The introduction of labels into existing solutions is however nontrivial as it cause a massive explosion of the number of subgraphs that need to be evaluated. In addition, many biological graphs allow multiple labels per node. Another issue when it comes to subgraph mining approaches for single graphs is that many of them struggle with large graphs, especially those with many hubs, i.e. highly connected nodes. As the input graph gets larger, so does the number of potential subgraphs of interests and for a lower frequency threshold, the number of patterns potentially explodes. Recently, there have been solutions that mine for frequent subgraphs in large graphs at low frequency thresholds [70]. There is also a question of the right interestingness measure, i.e. how to define an interesting subgraph from a bioinformatics point of view. The most frequent subgraphs are not necessarily the ones of interest and different alternative measures have been devised from a theoretical computer science perspective [71, 72]. Unfortunately, many of these interestingness measures are not as relevant for biomedical applications. There is thus still a need for novel algorithms which can find biologically relevant subgraphs in a flexible manner.
Abbreviations
 AGM:

Apriori graph mining
 BFS:

Breadthfirst search
 CAM:

Canonical adjacency matrix
 DFS:

Depthfirst search
 FPF:

Frequent pattern finder
 MIS:

Maximum independent set
References
 1
Koyutürk M, Grama A, Szpankowski W. An efficient algorithm for detecting frequent subgraphs in biological networks. Bioinformatics. 2004; 20(suppl 1):200–7.
 2
Hu H, Yan X, Huang Y, Han J, Zhou XJ. Mining coherent dense subgraphs across massive biological networks for functional discovery. Bioinformatics. 2005; 21(suppl 1):213–21.
 3
Cakmak A, Ozsoyoglu G. Mining biological networks for unknown pathways. Bioinformatics. 2007; 23(20):2775–83.
 4
Meysman P, Zhou C, Cule B, Goethals B, Laukens K. Mining the entire protein databank for frequent spatially cohesive amino acid patterns. BioData Min. 2015; 8(1):1.
 5
Jiang C, Coenen F, Zito M. A survey of frequent subgraph mining algorithms. Knowl Eng Rev. 2013; 28(01):75–105.
 6
Han J, Cheng H, Xin D, Yan X. Frequent pattern mining: current status and future directions. Data Min Knowl Disc. 2007; 15(1):55–86.
 7
Washio T, Motoda H. State of the art of graphbased data mining. SIGKDD Explor Newsl. 2003; 5(1):59–68.
 8
Fortin S. The graph isomorphism problem. Technical report, Technical Report 9620. Edomonton: University of Alberta; 1996.
 9
Inokuchi A, Washio T, Motoda H. An aprioribased algorithm for mining frequent substructures from graph data. In: Proceedings of the 4th European Conference on Principles of Data Mining and Knowledge Discovery. PKDD ’00: 2000. p. 13–23.
 10
Yan X, Han J. gspan: Graphbased substructure pattern mining. In: Proceedings of the 2002 IEEE International Conference on Data Mining. ICDM ’02: 2002. p. 721.
 11
Zaki MJ. Efficiently mining frequent trees in a forest. In: Proceedings of the Eighth ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. KDD ’02: 2002. p. 71–80.
 12
Asai T, Abe K, Kawasoe S, Arimura H, Sakamoto H, Arikawa S. Efficient substructure discovery from large semistructured data. In: Proceedings of the 2002 SIAM International Conference on Data Mining: 2002. p. 158–74.
 13
Pržulj N, Corneil DG, Jurisica I. Efficient estimation of graphlet frequency distributions in protein–protein interaction networks. Bioinformatics. 2006; 22(8):974–80.
 14
Pržulj N. Biological network comparison using graphlet degree distribution. Bioinformatics. 2007; 23(2):177–83.
 15
Hočevar T, Demšar J. A combinatorial approach to graphlet counting. Bioinformatics. 2013; 30(4):559–65.
 16
He H, Singh AK. Graphrank: Statistical modeling and mining of significant subgraphs in the feature space. In: Proceedings of the Sixth International Conference on Data Mining. ICDM ’06: 2006. p. 885–90.
 17
Ranu S, Singh AK. Mining statistically significant molecular substructures for efficient molecular classification. J Chem Inf Model. 2009; 49(11):2537–50.
 18
Ozdemir B, AbdAlmageed W, Roessler S, Wang XW. iSubgraph: integrative genomics for subgroup discovery in hepatocellular carcinoma using graph mining and mixture models. PloS one. 2013; 8(11):78624.
 19
Fiedler M, Borgelt C. Support computation for mining frequent subgraphs in a single graph. In: Proceedings of the 5th International Workshop on Mining and Learning with Graphs. MLG ’07: 2007.
 20
Bringmann B, Nijssen S. What is frequent in a single graph? In: Proceedings of the 12th PacificAsia Conference on Advances in Knowledge Discovery and Data Mining. PAKDD’08: 2008. p. 858–63.
 21
Vanetik N, Gudes E, Shimony SE. Computing frequent graph patterns from semistructured data. In: Proceedings of the 2002 IEEE International Conference on Data Mining. ICDM ’02: 2002. p. 458–65.
 22
Vanetik N, Shimony SE, Gudes E. Support measures for graph data. J Data Min Knowl Discov. 2006; 13(2):243–60.
 23
Calders T, Ramon J, Van Dyck D. Antimonotonic overlapgraph support measures. In: Proceedings of the Eighth IEEE International Conference on Data Mining. ICDM ’08: 2008. p. 73–82.
 24
Meysman P, Saeys Y, Sabaghian E, Bittremieux W, Van de Peer Y, Goethals B, Laukens K. Discovery of significantly enriched subgraphs associated with selected vertices in a single graph. In: Proceedings of the 14th International Workshop on Data Mining in Bioinformatics. BIOKDD ’15: 2015. p. 1–8.
 25
Leskovec J, Faloutsos C. Sampling from large graphs. In: Proceedings of the 12th ACM SIGKDD International Conference on Knowledge Discovery and Data Mining. KDD ’06: 2006. p. 631–6.
 26
Zou R, Holder LB. Frequent subgraph mining on a single large graph using sampling techniques. In: Proceedings of the Eighth Workshop on Mining and Learning with Graphs. MLG ’10: 2010. p. 171–8.
 27
Hu P, Lau WC. A survey and taxonomy of graph sampling. CoRR abs/1308.5865; 2013.
 28
Wörlein M, Meinl T, Fischer I, Philippsen M. A quantitative comparison of the subgraph miners mofa, gspan, ffsm, and gaston. In: Proceeding of the 9th European Conference on Principles and Practice of Knowledge Discovery in Database. PKDD ’05: 2005. p. 392–403.
 29
Huan J, Wang W, Prins J. Efficient mining of frequent subgraphs in the presence of isomorphism. In: Proceedings of the Third IEEE International Conference on Data Mining. ICDM ’03: 2003. p. 549–52.
 30
Nijssen S, Kok JN. The gaston tool for frequent subgraph mining. Electronic Notes Theor Comput Sci. 2005; 127(1):77–87.
 31
Dehaspe L, Toivonen H, King RD. Finding frequent substructures in chemical compounds. In: Proceedings of the Fourth International Conference on Knowledge Discovery and Data Mining. KDD’98: 1998. p. 30–6.
 32
Hassan M, Brown RD, VarmaO’Brien S, Rogers D. Cheminformatics analysis and learning in a data pipelining environment. Mol Divers. 2006; 10(3):283–99.
 33
Takigawa I, Mamitsuka H. Graph mining: procedure, application to drug discovery and recent advances. Drug Discov Today. 2013; 18(1):50–7.
 34
Borgelt C, Berthold MR. Mining molecular fragments: finding relevant substructures of molecules. In: Proceedings of the 2002 IEEE International Conference on Data Mining. ICDM’02: 2002. p. 51–8.
 35
Huan J, Wang W, Bandyopadhyay D, Snoeyink J, Prins J, Tropsha A. Mining protein family specific residue packing patterns from protein structure graphs. In: Proceedings of the Eighth Annual International Conference on Resaerch in Computational Molecular Biology. RECOMB ’04: 2004. p. 308–15.
 36
Dhifli W, Saidi R, Nguifo EM. Smoothing 3d protein structure motifs through graph mining and amino acid similarities. J Comput Biol. 2014; 21(2):162–72.
 37
Bandyopadhyay D, Huan J, Liu J, Prins J, Snoeyink J, Wang W, Tropsha A. Structurebased function inference using protein familyspecific fingerprints. Protein Sci. 2006; 15(6):1537–43.
 38
Ramu A, Kahveci T, Burleigh JG. A scalable method for identifying frequent subtrees in sets of large phylogenetic trees. BMC Bioinforma. 2012; 13(1):256.
 39
Zhang S, Wang JT. Discovering frequent agreement subtrees from phylogenetic data. IEEE Trans Knowl Data Eng. 2008; 20(1):68–82.
 40
Deepak A, FernándezBaca D, Tirthapura S, Sanderson MJ, McMahon MM. Evominer: frequent subtree mining in phylogenetic databases. Knowl Inf Syst. 2014; 41(3):559–90.
 41
Ayed R, Hacid MS, Haque R, Jemai A. An updated dashboard of Complete Search FSM Implementations in Centralized Graph Transaction Databases: BDA 2016 Gestion de DonnéesPrincipes, Technologies et Applications 32 e anniversaire 1518 novembre 2016, Poitiers, Futuroscope; 2016, p. 57.
 42
Inokuchi A. A fast algorithm for mining frequent connected subgraphs. IBM Research Report RT0448. 2002.
 43
Kuramochi M, Karypis G. Frequent subgraph discovery. In: Data Mining, 2001. ICDM 2001, Proceedings IEEE International Conference On: 2001. p. 313–20.
 44
Kuramochi M, Karypis G. An efficient algorithm for discovering frequent subgraphs. IEEE Trans Knowl Data Eng. 2004; 16(9):1038–51.
 45
Zaki MJ, Parthasarathy S, Ogihara M, Li W, et al. New algorithms for fast discovery of association rules. In: KDD: 1997. p. 283–6.
 46
Milo R, ShenOrr S, Itzkovitz S, Kashtan N, Chklovskii D, Alon U. Network motifs: simple building blocks of complex networks. Science. 2002; 298(5594):824–7.
 47
Huang CY, Cheng CY, Sun CT. Bridge and brick network motifs: Identifying significant building blocks from complex biological systems. Artif Intell Med. 2007; 41(2):117–27.
 48
Davis D., Yaverğlu ÖN, MalodDognin N, Stojmirovic A, Pržulj N. Topologyfunction conservation in protein–protein interaction networks. Bioinformatics. 2015; 31(10):1632–9.
 49
Cootes AP, Muggleton SH, Sternberg MJ. The identification of similarities between biological networks: application to the metabolome and interactome. J Mol Biol. 2007; 369(4):1126–39.
 50
Flannick J, Novak A, Srinivasan BS, McAdams HH, Batzoglou S. Graemlin: general and robust alignment of multiple large interaction networks. Genome Res. 2006; 16(9):1169–81.
 51
Ruyssinck J, Demeester P, Dhaene T, Saeys Y. Netter: reranking gene network inference predictions using structural network properties. BMC Bioinforma. 2016; 17(1):1.
 52
Ciriello G, Guerra C. A review on models and algorithms for motif discovery in protein–protein interaction networks. Brief Funct Genom Proteomics. 2008; 7(2):147–156.
 53
Aittokallio T, Schwikowski B. Graphbased methods for analysing networks in cell biology. Brief Bioinform. 2006; 7(3):243–55.
 54
Alon U. Network motifs: theory and experimental approaches. Nat Rev Genet. 2007; 8(6):450–61.
 55
Meysman P, Titeca K, Eyckerman S, Tavernier J, Goethals B, Martens L, Valkenborg D, Laukens K. Protein complex analysis: From raw protein lists to protein interaction networks. Mass Spectrom Rev. 2015; 36(5):600–14.
 56
Kashtan N, Itzkovitz S, Milo R, Alon U. Efficient sampling algorithm for estimating subgraph concentrations and detecting network motifs. Bioinformatics. 2004; 20(11):1746–58.
 57
Krause AE, Frank KA, Mason DM, Ulanowicz RE, Taylor WW. Compartments revealed in foodweb structure. Nature. 2003; 426(6964):282–5.
 58
Tran NTL, Mohan S, Xu Z, Huang CH. Current innovations and future challenges of network motif detection. Brief Bioinform. 2014; 16(3):497–525.
 59
Schreiber F, Schwöbbermeyer H. Towards motif detection in networks: Frequency concepts and flexible search. In: in Proceedings of the International Workshop on Network Tools and Applications in Biology (NETTAB04): 2004. p. 91–102.
 60
Wernicke S. A faster algorithm for detecting network motifs. In: Proceedings of WABI ’05: 2005. p. 165–77.
 61
Kashani ZR, Ahrabian H, Elahi E, NowzariDalini A, Ansari ES, Asadi S, Mohammadi S, Schreiber F, MasoudiNejad A. Kavosh: a new algorithm for finding network motifs. BMC Bioinforma. 2009; 10(1):318.
 62
Li X, Stones DS, Wang H, Deng H, Liu X, Wang G. NetMODE: Network Motif Detection without Nauty. PloS one. 2012; 7(12):50093.
 63
Meira LA, Máximo VR, Fazenda ÁL, Da Conceição AF. AccMotif: Accelerated Network Motif Detection. IEEE/ACM Trans Comput Biol Bioinforma (TCBB). 2014; 11(5):853–62.
 64
PeregrínAlvarez JM, Xiong X, Su C, Parkinson J. The modular organization of protein interactions in escherichia coli. PLoS Comput Biol. 2009; 5(10):1000523.
 65
Hellal A, Romdhane LB. Nodar: mining globally distributed substructures from a single labeled graph. J Intel Inf Syst. 2013; 40(1):1–15.
 66
Jiang X, Xiong H, Wang C, Tan AH. Mining globally distributed frequent subgraphs in a single labeled graph. Data Knowl Eng. 2009; 68(10):1034–58.
 67
Wörlein M, Dreweke A, Meinl T, Fischer I, Philippsen M. Edgar: the embeddingbased graph miner. In: Proceedings of the International Workshop on Mining and Learning with Graphs : in Conjunction with ECML / PKDD, MLG ’06: 2006. p. 1034–58.
 68
Ketkar NS, Holder LB, Cook DJ. Subdue: Compressionbased frequent pattern discovery in graph data. In: Proceedings of the 1st International Workshop on Open Source Data Mining: Frequent Pattern Mining Implementations. OSDM ’05: 2005. p. 71–6.
 69
Ghazizadeh S, Chawathe SS. Seus: Structure extraction using summaries. In: Proceedings of the 5th International Conference on Discovery Science. DS ’02: 2002. p. 71–85.
 70
Elseidy M, Abdelhamid E, Skiadopoulos S, Kalnis P. Grami: Frequent subgraph and pattern mining in a single large graph. Proc VLDB Endowment. 2014; 7(7):517–28.
 71
Cule B, Goethals B, Hendrickx T. Mining interesting itemsets in graph datasets. In: Advances in Knowledge Discovery and Data Mining. 17th PacificAsia Conference, PAKDD 2013: 2013. p. 237–48.
 72
van Leeuwen M, De Bie T, Spyropoulou E, Mesnage C. Subjective interestingness of subgraph patterns. Mach Learn. 2016; 105(1):41–75.
 73
Su C, PeregrinAlvarez JM, Butland G, Phanse S, Fong V, Emili A, Parkinson J. Bacteriome.org–an integrated protein interaction database for e. coli. Nucleic Acids Res. 2008; 36(suppl 1):632–6.
 74
GamaCastro S, Salgado H, SantosZavaleta A, LedezmaTejeida D, MuñizRascado L, GarcíaSotelo JS, AlquiciraHernández K, MartínezFlores I, Pannier L, CastroMondragón JA, et al. RegulonDB version 9.0: highlevel integration of gene regulation, coexpression, motif clustering and beyond. Nucleic Acids Res. 2016; 44:133–43.
 75
Moretto M, Sonego P, Dierckxsens N, Brilli M, Bianco L, LedezmaTejeida D, GamaCastro S, Galardini M, Romualdi C, Laukens K, et al. Colombos v3. 0: leveraging gene expression compendia for crossspecies analyses. Nucleic Acids Res. 2016; 44(D1):620–3.
 76
Inokuchi A. AcGM. http://ist.ksc.kwansei.ac.jp/~inokuchi/acgm.zip. Accessed 24 Nov 2017.
 77
Karypis G. FSG. http://glaros.dtc.umn.edu/gkhome/pafi/overview. Accessed 24 Nov 2017.
 78
Parmol. http://en.verysource.com/parmol_1346_200608100787.html. Accessed 24 Nov 2017.
 79
Borgelt C. MOSS. http://www.borgelt.net/moss.html. Accessed 24 Nov 2017.
 80
Alon U. Mfinder. http://www.weizmann.ac.il/mcb/UriAlon/download/networkmotifsoftware. Accessed 24 Nov 2017.
 81
Schreiber F, Schwöbbermeyer H. Mavisto: a tool for the exploration of network motifs. Bioinformatics. 2005; 21(17):3572–4.
 82
Schreiber F, Schwöbbermeyer H. MAVisto. http://mavisto.ipkgatersleben.de/. Accessed 24 Nov 2017.
 83
Wernicke S, Rasche F. Fanmod: a tool for fast network motif detection. Bioinformatics. 2006; 22(9):1152–3.
 84
Wernicke S, Rasche F. FANMOD. http://theinf1.informatik.unijena.de/motifs/. Accessed 24 Nov 2017.
 85
Kavosh. https://github.com/shmohammadi86/Kavosh. Accessed 24 Nov 2017.
 86
NetMODE. https://sourceforge.net/projects/netmode/. Accessed 24 Nov 2017.
 87
accMotif. http://www.ft.unicamp.br/docentes/meira/accmotifs/. Accessed 24 Nov 2017.
 88
Meysman P. SSM. http://adrem.ua.ac.be/sigsubgraph. Accessed 24 Nov 2017.
Funding
This work was supported by the Fund for Scientific Research – Flanders (FWOVlaanderen) project “Evolving graph patterns”. W.B. is a postdoctoral researcher of the Research Foundation – Flanders (FWO).
Availability of data and materials
The datasets analysed during the current study are available in the Git repository, https://bitbucket.org/ligataaida/data.
Author information
Affiliations
Contributions
PMe and KL developed the concept, AM and PMo performed the analysis, AM, PMe, WB wrote the manuscript. WB, BC, BG gave conceptual advice. All authors read and approved the final manuscript.
Corresponding author
Correspondence to Kris Laukens.
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interests
The authors declare that they have no competing interest.
Publisher’s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional files
Additional file 1
More detailed definitions of graphs and subgraphs. (PDF 385 kb)
Additional file 2
Output of tools for both single and multiple graph settings. (PDF 286 kb)
Rights and permissions
Open Access This 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.
About this article
Received
Accepted
Published
DOI
Keywords
 Subgraph mining
 Frequent subgraphs
 Graph motifs
 Biological networks
 Pattern discovery
 Pattern mining