Skip to main content
  • Research article
  • Open access
  • Published:

Predicted protein-protein interactions in the moss Physcomitrella patens: a new bioinformatic resource



Physcomitrella patens, a haploid dominant plant, is fast becoming a useful molecular genetics and bioinformatics tool due to its key phylogenetic position as a bryophyte in the post-genomic era. Genome sequences from select reference species were compared bioinformatically to Physcomitrella patens using reciprocal blasts with the InParanoid software package. A reference protein interaction database assembled using MySQL by compiling BioGrid, BIND, DIP, and Intact databases was queried for moss orthologs existing for both interacting partners. This method has been used to successfully predict interactions for a number of angiosperm plants.


The first predicted protein-protein interactome for a bryophyte based on the interolog method contains 67,740 unique interactions from 5,695 different Physcomitrella patens proteins. Most conserved interactions among proteins were those associated with metabolic processes. Over-represented Gene Ontology categories are reported here.


Addition of moss, a plant representative 200 million years diverged from angiosperms to interactomic research greatly expands the possibility of conducting comparative analyses giving tremendous insight into network evolution of land plants. This work helps demonstrate the utility of “guilt-by-association” models for predicting protein interactions, providing provisional roadmaps that can be explored using experimental approaches. Included with this dataset is a method for characterizing subnetworks and investigating specific processes, such as the Calvin-Benson-Bassham cycle.


A compilation of the physical interactions between the proteins of an organism is known as the protein-protein interactome. It is through these interactions that most biological processes within a living cell are performed [1]. While interactomes can be a composed from many individual published experiments, several high-throughput methods have been developed to rapidly identify interactions between proteins even those with no known function. In doing so, high-throughput methods reveal much of the signalling and communication within the proteomes of yeast, human, fruit fly, and nematode worm [2-5]. With 84% of the proteins interactions experimentally determined, the yeast interactome map is nearly complete [6,7]. However the same is not yet true for the other model organisms, especially plants where high-throughput methods have only recently begun to be employed [8].

To fill in these gaps, predicted interactome maps have been constructed for Arabidopsis thaliana (herein referred to as Arabidopsis), rice, coffee, yeast, human, mouse, and fruit fly at the organelle and whole cell levels using computational interolog methodology [6,7,9-13]. The interolog method of predicting protein interactions is based on the hypothesis that conserved orthologous genes encode functionally similar proteins forming similar complexes and signalling pathways [1,14]. This method assumes that functionally similar proteins retain their interacting partners. Using interologs, a predicted protein interactome can be constructed without the extreme expense or time requirement of high-throughput methods, which currently have only been done for selected model organisms. This provisional roadmap of the predicted interactions is by nature of the methodology incomplete, and only contains interactions among proteins that are highly conserved.

The model moss Physcomitrella patens is fast becoming a tool for bioinformatics and molecular work due to its key phylogenetic position as sister to the other land plant lineages. Physcomitrella patens has a protein-coding genome similar in size to Arabidopsis thaliana, but is similar to yeast in efficiency of gene targeting experiments and haploid-dominance of the life cycle, making an interesting and useful molecular genetic tool for plants [15-17]. This moss genome has been sequenced and annotated based on sequence homology to known genes and domains [18,19]. The annotated moss genome is a key tool for research into the evolution of all plant functions. As a bryophyte, this model moss is the placeholder for early land plants on the tree of life, making evolutionary comparisons of biological pathways at the protein interaction level a useful new avenue of investigation.

Presented here is the first predicted PPI for a bryophyte based on the interolog method. Bias and enrichment of gene functions in the predicted interactome were analysed in order to help the user evaluate the utility and interpret the results of this tool. Although derived from many reference organisms, only a few plant and cyanobacterial reference interactions were available for orthology mapping. A confidence value (CV) was assigned to each interaction in order to help determine credibility of predicted interactions. Resolution of the plant specific protein networks, the Calvin-Benson-Bassham (CBB) cycle network, is shown to illustrate the bait and prey methods of capturing functional subnetworks. This tool provides only a conserved eukaryotic skeleton of biological pathways and interactions, but can aid current research and provide the framework for future avenues of interactomic research in plants.


Construction of MySQL database

Creating the merged interaction database was a four-step process beginning with the development of standardized identifiers for each interactor. Secondly, conversions from non-chosen identifiers that exist in four interaction databases (BIND, DIP, BioGrid, and IntAct) to the selected identifiers were found and placed inside SQL databases [20-23].

Third, a Universal Translator program for reference interactomes (see Additional file 1 for source code) was developed to create mappings between the three chosen identifiers, Swiss/Uniprot ID, Entrez Gene ID, and Ensembl Peptide ID. These were chosen because there are a large number of already existing mappings of these identifiers. BioMart [24] was relied upon to create the mappings because of the amount and quality of information found in the databases available at the BioMart website ( The data were output into a highly configurable format that was input into our local database. The Universal Translator program used data from the Uniprot website as well as the four interaction databases to create a mappings database between Uniprot, Entrez Gene, and Ensembl IDs. This database allows for the Interactome Merger application to operate much quicker.

Finally, an Interactome Merger application (source code also in Additional file 1) was developed to pull data from the four interaction databases and, using the database created by the Universal Translator and the mappings, make a merged interaction database. Interactome Merger was responsible for using the data in the Universal Translator database and mappings databases to convert interactions contained in DIP, BIND, BioGrid, and IntAct into a standardized format where each interaction had at least one of the three chosen identifiers associated to each interactor.

The reference database was queried for interactions for which moss orthologs existed for both interacting partners. Output included source interactions, referenced metadata including placeholder ID and moss predicted interactions. This data spreadsheet was sorted for unique interactions and contains predicted orthologous protein IDs, type of experiment, organism the interaction was found, source database the interaction data were retrieved, and the PubMed ID for each interaction (Additional file 2: Table S1). The general process for assembling this database can be outlined visually using a flow-chart (Figure 1).

Figure 1
figure 1

A flow-chart can be used to visualize the process of generating the predicted interactome. The predicted interactome of Physcomitrella patens was derived from orthologs of Arabidopsis, nematode, fruitfly, bacteria, mouse, rat, human and yeast using the InParanoid algorithm (See Methods). One to one orthology was used to query a MySQL database containing interaction data from BioGrid, BIND, DIP and Intact databases. A spreadsheet of orthologous interactions and their supporting information was generated and input into Cytoscape for the predicted interactome visualization.

Construction of predicted moss interactome

Interologs were determined using InParanoid ortholog predicting algorithm that compares, in pairwise manner, protein sequences derived from the gene models of Physcomitrella patens genome assembly V1.2 and annotation V1.6 from COSMOSS [19] against the gene models of E. coli, yeast, nematode worm, fruit fly, rat, mouse, and human [25,26]. The program separated outparalogs and orthologs into different clusters. These orthologs were loaded into a MySQL database that housed interactome datasets from BIND, BIOGRID, DIP, and IntAct databases [20-23]. One pair of proteins (moss, reference species) from each outparalog and ortholog cluster was selected to form the ortholog/outparalog one to one match (herein referred to as ortholog for brevity). A complete table of orthologs for all reference species is included as Additional file 3: Table S2.

The raw interolog data, unique interactions and calculation matrix for confidence value was maintained as a spreadsheet (Additional file 2: Table S1). This file contains 104,392 raw predicted interactions for Physcomitrella, reference species, type of methodology, and PubMed ID for the reference interaction. Moss protein IDs were provided for COSMOSS V1.2 and V1.6 annotation [19,27]. Current Uniprot, Genbank and other IDs (where available) can be obtained through Genonaut ( Duplicate interactions were removed to produce a dataset of 67,740 predicted interactions with the number of interactions for each protein, this constituted the unique moss interactome.

The confidence value (CV) for each unique interaction was calculated in the same way as Geisler-Lee et al. [11]. Essentially the CV is the total number of times the interaction is predicted from different experimental references, with a multiplier added for different experimental methods, and different reference species. Thus the CV of an interaction with evidence from 4 different references using 2 different experimental methods and 3 different species would be 4 × 2 × 3 = 24, while one with 4 references all using the same experimental method in the same species would be 4 × 1 × 1 = 4. This metric thus favours different experimental methodology and evolutionary conservation of the interaction.

Added to this spreadsheet were the KOG descriptions and definitions for each protein, which comprise the node attributes for the moss interactome. These data were loaded into Cytoscape 3.0.2 [28,29] and the predicted protein-protein interactome was constructed and displayed in the organic layout, thus rendering the ball and stick interactome model.

Interactome analysis

Analysis of the overrepresented GO categories was conducted using the BiNGO v2.41 [30] plugin available through the Cytoscape web interface. The following settings were used to generate the network visualization; get the cluster from network and visualize overrepresented clusters after correction, perform hypergeometric test against the whole annotation with the Benjamini & Hochberg False Discovery Rate correction at the 0.05 significance level. Since the main concern was with the categories associated with plants, the GO Slim Plant ontology file was selected for cross-referencing against our custom Moss GO Annotation file (Additional file 4: Table S3). The same settings were employed to determine the underrepresented clusters. Hub distribution analysis and comparison used the Arabidopsis interactome from Geisler-Lee et al. [11] and the human and yeast interactomes downloaded from BioGrid [23].


Identification of moss gene orthologs

Protein sequences from whole genomes of reference species were compared to those in the genome of Physcomitrella patens in a pairwise fashion using reciprocal BLASTs to separate inparalogs from orthologs and outparalogs with the InParanoid software package [25,26,31]. Of the nearly 28,000 genes in moss, only a few had true orthologs in yeast and animal reference species (Additional file 3: Table S2). The majority of orthologs present in the predicted interactome were from animal or fungus. Where multiple possible orthologs occurred, only one was selected for strict one-to-one orthology. The inparalogs were deliberately left out of the interaction prediction. The reasoning is that the last common species ancestor for most reference genomes is very early in the evolution of eukaryotes, thus each lineage has had enormous amount of time to diverge. Where such gene duplication has occurred, the resulting lack of selection pressure was thought likely lead to strong protein/gene neofunctionalization or subfunctionalization (or gene loss). Given the long evolutionary time, it was considered unlikely that a completely redundant function will remain, and thus cannot expect that such divergent inparalogs would have all the same interacting partners. Each moss gene was thus partnered with a single ortholog from each of the reference species where they occurred; this was assembled in MySQL and constituted the orthology database.

Assembly of the reference database and interolog prediction

More than 104,000 total interactions were predicted from different references including 67,740 unique interactions from 5,694 different Physcomitrella patens proteins (Additional file 2: Table S1). When visualized in its entirety in Cytoscape v3.2.0, the interactome looks like a ball of densely tangled circles and lines (Figure 2a,b). Like a map of roads and towns, this interactome is most useful when plotting connections between genes or pathways of interest and their neighbours, which is made possible in Cytoscape by the select node and zoom tools. Individual interactions were systematically evaluated using an evidence-based confidence value to assess the quality of the predicted interactions. 3936 high confidence (CV > 10), 10,318 medium confidence (CV between 2 and 10), and 53,400 low confidence (CV = 1) interactions were found in moss. The prediction efficiency of different confidence levels was evaluated by comparing the current experimentally determined interactions accumulated at the BioArray Resource [32] to the predicted interactions using the same interlog method and CV calculation for Arabidopsis by Geisler-Lee et al. [11]. The BAR resource combined the work of several high throughput experimental studies [8,33] with individual experiments culled from over 1190 publications (full list available at The 37645 experimentally determined interactions in Arabidopsis were compared to the 72266 predicted interactions in Arabidopsis from non-plant reference genomes, and 1450 matched. The expected overlap between these datasets by chance alone (equal number of random protein pairs) was 91. This did not account for non-overlapping protein sets (proteins present in one set but absent in the other), nor did it account for the non-overlapping bias, as experimental interactomes tend to be biased towards plant specific genes, of which none would have occurred in the Arabidopsis predicted interactome. When broken down by confidence value, there was only slight improvement comparing interactions with low CV to those with medium CV, and there was 2.1 fold enrichment when comparing interactions with low and high confidence values.

Figure 2
figure 2

The Physcomitrella patens network was viewed in Cytoscape and the hub distribution of proteins in the network were analysed. (A) Large hairy ball of the 67,740 non-redundant interactions in the organic view of Cytoscape. (B) Detailed view showing layers of interacting proteins. (C) Pie chart of hub distribution by node class; “Free ends” are defined as have a single interaction, “pipes” have two interactions, and grouped hubs based on numbers of interactions. Numbers in parentheses indicate number of interacting proteins used to determine hub size.

The topology of the whole interactome was evaluated by dividing it into groups of proteins based on the number interactions (connectivity in Figure 2c). Greater than half the proteins are major hubs with between 51 and 100 interactions. In comparison, hub distributions in yeast and Arabidopsis are primarily medium in size between 11 and 50 interactions while the human protein interactome has more minor hubs (Figure 3a,c,e,g). Regardless of hub distributions in the interactomes of moss, Arabidopsis, yeast, and human, protein connectivity follows a scale free power law distribution (Figure 3b,d,f,h).

Figure 3
figure 3

Frequency distributions of hub sizes among Physcomitrella patens , Arabidopsis thaliana , Saccharomyces cerevisiae , and Homo sapiens . Graphs for the frequency distribution of hub (protein) connectivity in shown both as a bar chart by node classes (left) and binned every 10 connections (right) in A, B: Physcomitrella patens (predicted, this work), C, D: Arabidopsis thaliana (predicted by [11]), E, F: Saccharomyces cerevisiae (experimentally determined; from BioGrid), and G, H: Homo sapiens (experimentally determined; from BioGrid).

Large hub sizes for proteins in both predicted and experimentally determined interactomes may be attributed to several factors. The most common high throughput experimental techniques used in reference organisms included Affinity Capture and Yeast Two Hybrid (Figure 4a,b), which can generate experimental false positive results [34]. This prompted the use of an experiment multiplier (E) when calculating confidence values. Additionally, conserved proteins tend to have more interacting partners and tend to be more essential [35], and this might skew the structure of the network composed of mostly conserved proteins towards mid to large hub sizes. Finally, it might indeed be true that this many potential interactions are possible, but that only a subset of these interactions actually occur in any given living cell or tissue due to differential expression of genes encoding the proteins. The distribution of species investigated was determined by counting the number experiments for each species, most coming from fungi (yeast) or animal references (Figure 4c). Approximately 5,130 moss predicted interactions that come from plant and cyanobacterial reference organisms (i.e. Arabidopsis and Synechocystis) were detected. Many of these interactions are highly conserved across eukaryotes (Additional file 5: Table S4), making evolutionary comparisons of plant networks feasible.

Figure 4
figure 4

A frequency distribution of interactions by experiment type and organism was analyzed. Histogram showing numbers of interactions by type of experimental system. Major contributions shown in (A) and a smaller scale of contributions shown in (B). Code abbreviations are as follows: NR, Not reported; AC, Affinity Capture; NG, Negative Genetic; Y2H, Yeast-two-hybrid; PE, Phenotypic Enhancement; SL, Synthetic Lethality; SGD, Synthetic Growth Defect; DR, Dosage Rescue; PS, Phenotypic Suppression; PG, Positive Genetic; RC, Reconstituted Complex; Co-P, Co-purification; SR, Synthetic Rescue; BA, Biochemical Activity. PCA, Principal Components Analysis; Co-F, Co-fractionation; Co-CS, Co-crystal Structure; Co-L, Co-localization; DL, Dosage Lethality; FW, Far Western; SH, Synthetic Haploinsufficiency; DGD, Dosage Growth Defect; Prot, Protein/Peptide; FRET. (C) Histogram of interactions by organism type.

Highly-connected nodes and conserved eukaryotic orthologs

A search for the most connected proteins within the moss interactome revealed the largest hubs (greatest number of interacting partners) are linked to each other with high degrees of confidence (Figure 5). Many of these hubs are involved with the ribosome, nuclear DNA repair, ubiquitin, proteosome, and the cytoskeleton. These likely represent strongly conserved pathways that have not altered significantly during eukaryotic evolution. Heat shock proteins (HSP 90 and HSP 60) were also among the most highly connected proteins, possibly due to their ubiquitous involvement in forming dimers with numerous other proteins. An interesting node is Pp1s323_19V6.1, a GDP-mannose pyrophosphorylase enzyme that converts mannose to GDP-mannose used in cell wall biosynthesis, protein glycosylation and ascorbic acid biosynthesis. This protein appears to have many more predicted interacting partners (237; Additional file 6: Table S5) in moss than the ortholog for Arabidopsis CYT1/AT2G39770 (88 partners) [11,32]. Of the additional interacting partners in moss most came from references not used in the construction of the Arabidopsis predicted interactome, but 37 interacting partners had no Arabidopsis ortholog. These represent moss genes with no matching ortholog in Arabidopsis, but with orthologs other reference species (i.e. yeast, human, E. coli).

Figure 5
figure 5

The top twenty hubs in the Physcomitrella interactome were elucidated from the network. The hubs represent the most highly connected proteins found within the dataset.

The moss interactome, like other partially completed interactomes for Arabidopsis, human, drosophila and rice is surprisingly well connected, with most proteins assembled into a single interacting ball. One early theory for the possible global organisation of biological pathways is that of densely connected subnetworks that are only connected to each other due to the presence of highly connected hubs such as HSPs; this is referred to as the “party and date” hub model. A second model, observed in yeast, is that there is dense interconnection throughout the entire proteome and few or no separate subnetworks connected by date hubs; this is known as the stratus (cloud) model [36-39]. The robustness of the interactome was tested by deleting the top 20 connected proteins to determine if the moss interactome core fits the party-date model or the stratus model. The deletion of these proteins had no overall effect on the interactome topology. All remaining proteins stayed connected and the distribution of hub size was not affected significantly except for the top twenty.

The top one hundred most conserved interactions were examined across reference eukaryotes. Like the most conserved hubs, the conserved interactions were found to represent core biological processes such as DNA transcription and polymerization, RNA polymerization, cell cycle, and vesicle associated interactions (Additional file 5: Table S4). However, within this list of fifty interactions there are ten proteins with completely unknown functions and among those, two of the unknown proteins interact with each other. This occurs due to the high throughput experimental methodologies tend to investigate interactions without knowing what the proteins do (i.e. somewhat double blind). Thus there are still core-conserved processes that are not understood thus proteins involved in these are good candidates for future research of significance to all of biological science.

Analysis of gene ontology

The conservation of proteins using the interlog method will potentially bias the resulting network towards conserved processes. To learn what these conserved processes were, genes in the interactome were annotated by Gene Ontology, using a reference list of approximately 17,000 GO annotations for Physcomitrella patens [19]. An analysis for enriched and depleted processes revealed that proteins involved in intracellular and cytoplasmic metabolic and catalytic processes are overrepresented likely due to their conserved nature. Additionally the protein binding category was enriched due to physical interaction requirement for inclusion in the interactome (Figure 6, Additional file 7: Table S6). Within the cellular and metabolic processes, a majority of overrepresented proteins involved DNA and protein metabolism, likely due to these processes being conserved in all life and consequently being the most studied in the reference organisms. Also cellular component organization and biogenesis are overrepresented in our interactome. There were relatively few, but significant, representation of proteins involved in cell cycle, growth, embryonic development, cell differentiation, and behavior. The catalytic and binding activities most represented in the interactome relate to hydrolase and transferase enzymes, as well as nucleic acid, protein, and nucleotide binding. A complete protein list from the analysis is provided in Additional file 5: Table S4.

Figure 6
figure 6

A BiNGO analysis of the Physcomitrella patens predicted interactome was performed in order to determine over-represented cellular processes. Node size is related to the number of times a protein in that category is found in the predicted interactome. Node color indicates level of significance (i.e. p-value) at the 0.05 level. (See Additional file 5: Table S4 for complete protein list and their p-values).

Orthology and comparison of moss and Arabidopsis genomes

Physcomitrella appears to have more orthologs in general, although this value does not represent a significant increase. For example: moss has 2,999 orthologous proteins to human and 1,793 orthologous proteins to yeast, while Arabidopsis has 2,708 and 1,708 proteins respectively. Of the 2,999 human orthologs moss and Arabidopsis share 1,485 proteins, but moss has 896 proteins with unique orthology to human proteins. Most interesting though, is the orthology of arachidonate 5-lipoxygenase (5-LOX) in moss and human. This enzyme in humans adds oxygen to fatty acids for form arachidonic acids (animal fat) that are implicated in prostate tumor growth and asthma [40]. Arachidonic acids are significantly less abundant in angiosperms with respect to linoleic and α-linoleic acids, while the inverse is the condition in algae and bryophytes [41]. However arachidonate 12-lipoxygenases (12-LOX), not 5-LOX, significantly prefer arachidonic acid as a substrate. Therefore we speculate that the 5-LOX enzyme activity is reduced due to purifying selection for 12-LOX in Physcomitrella patens. This discovery highlights that mosses and other bryophytes may possess enzymes and pathways that are shared with non-plant eukaryotes, but are reduced or absent in angiosperms or vascular plants. As each sister lineage of land plants diverged from their last common ancestor, they may have selectively lost different conserved pathways. These differences are detected by ortholog and interolog analysis at the whole genome level.


Topology and features of the predicted interactome

A predicted protein interactome for Physcomitrella patens was constructed using the interolog method based on orthologous genes across 14 reference organisms. Nearly 68,000 unique interactions were recovered from 5,695 proteins in Physcomitrella patens. The increase in the number of reference species more than doubled the number of interactions, and added just over 400 more proteins than were found in Arabidopsis using the same methods [11]. Increased taxon sampling resulted in increased CV for the most highly conserved interactions. Interactions with the highest CV are typically those involving proteins from highly conserved processes such as DNA replication and polymerization, RNA polymerization, endomembrane trafficking, metabolism and mitosis. Consequently accurate prediction of protein interaction pathways is possible for these processes in moss. Additionally proteins involved in formation of the proteosome complexes are overrepresented in the moss interactome. This both points out the conserved nature of this complex and illustrates the limitation of this methodology. As more interactomic data become available from reference species, there is greater confidence in the predicted proteins involved in deeply conserved processes [1,2,6,7,9,11,13].

Functional annotation of most (>99%) proteins in genome databases is sequence homology based without any experimental data. This work adds to this annotation by inference based on the presence of known interacting partners. When two unknown proteins interact with each other such inferences are not possible, however the co-occurrence of these proteins across several taxa makes it interesting because whole processes, completely unknown as to their function, can be conserved across eukaryotes. This implies that there are still fundamental biological pathways as yet unexplored. The top hubs (i.e. proteins most interacting partners) represent well-known fundamental biological processes (Figure 5). Given their importance to the interactome one would expect removal of these proteins would lead to a breakdown into several disconnected networks. Yet removal of these highly connected nodes did not result in such a breakdown. These results suggest the moss predicted protein interactome fits a stratus structure of connectivity [36]. This type of network is composed of densely connected hubs that have high degrees of overlap of interacting proteins. This differs from the view that protein interaction networks exhibit a party/date hub or altocumulus structure where subnetworks interact through a very few (date) hubs [39]. Protein interactomes that fit a stratus structure rather than an altocumulus structure support data demonstrating that many proteins have multiple functions or participate in multiple complexes [37,38]. However high density connectedness found in stratus-like interactomes are less amenable to finding specific modules within the structure [36].

When the general features were looked at for the moss, Arabidopsis, yeast and human interactomes a key difference was discovered in the distribution of hubs. Physcomitrella patens has more major hubs (51–100 interactions) while yeast and human have more minor, small, and medium hubs. The Arabidopsis interactome connectivity was more similar to the moss genome than it was to either yeast or human. Moss has fewer free-ends and pipes when compared to Arabidopsis, but this could due to the fewer number of reference species used to generate the Arabidopsis dataset [11]. Since hub size is directly associated with numbers of interacting proteins, the increased number of reference species increases the number of predicted interacting proteins. The increased number of interactions would result in a shift in the types of hubs represented in the moss interactome. Regardless of this shift, hub connectivity follows a scale-free power law distribution (Figure 3), a feature found in other interactomes, metabolic networks, and neural networks [42].

Another feature of Physcomitrella patens is that it is a paleopolyploid as result from an ancient genome duplication that preferentially maintained a large complement of metabolic genes [27,43] encoding for proteins that are overrepresented in the predicted interactome. While not significantly different, in general moss appears to have more orthologs with the reference species than does Arabidopsis. This is expected given that the genomes of Physcomitrella patens and Arabidopsis have very similar numbers of protein coding sequences. However the common ancestor to the Brassicales also underwent a genome duplication event followed by subsequent genome reductions in Arabidopsis resulting in biased retention toward transcriptional regulation and signal transduction [44].

“Guilt-by-association” and analysis of subnetworks

The utility of predicted interactomes lies in the “guilt-by-association” model of predicting proteins in a pathway under the assumption that orthologous proteins have similar functions [2,4,11,45]. Roadmaps are provided for all the conserved pathways in moss and assign pathway position and function through interactions. This was demonstrated, with high confidence, in predicting disease genes in humans [46,47].

This Physcomitrella interactome can be used to analyze specific subnetworks and further characterize key interacting proteins. Using the interactome, one can view interactions between proteins that are known to be involved in a certain processes (bait) and others that directly interact with the bait (prey). Of all known CBB-associated proteins, 23 were found to exist in the dataset and were chosen as bait resulting in a prey capture 242 proteins (Figure 7, Additional file 8: Table S7). Finding key prey proteins can be accomplished by look at the ratio of that protein’s connectivity within the subnetwork to its overall connectivity. Prey proteins that exclusively occur in this subnetwork include a carbohydrate transmembrane transporter (Pp1s88_88V6.1), isocitrate dehydrogenase (Pp1s78_143V6.1), glycosyltransferase (Pp1s90_132V6.1), UDP-acetylglucosamine acyltransferase (Pp1s377_21V6.1), and a protein with a domain of unknown function (Pp1s30_41V6.1). The average connectivity ratio in this subnetwork was 0.33 suggesting that many of these proteins are specific to this process. Reconstruction of subnetworks is aided through this “bait and prey” model and can help elucidate mechanisms that emerge from the network’s activity.

Figure 7
figure 7

A bait and prey method was used to reconstruct a sub-network using the Physcomitrella predicted interactome. Two proteins known to be involved in the Calvin-Benson-Bassham cycle were used as bait (blue color nodes) and their first neighbors are the captured prey (green color nodes). These two known CBB-cycle proteins do not directly interact, but interact with several other proteins that serve as intermediates. Additional file 8: Table S7 contains the complete protein list and annotations used for this analysis.


We developed the moss interactome to increase the community resource base by improving the annotation of moss proteins by adding predicted interactions and allowing guilt-by-association functional annotation. The ability to view an interactome in its entirety facilitates the verification of known interacting proteins and discovery of new proteins involved in a particular pathway or process. The availability of an annotated interactome for the moss provides an invaluable resource and focal point for system-scale analysis of this bryophyte model organism. The addition of moss, a plant representative of about 450 million years divergent evolution from seed plants like Arabidopsis, to interactomic research greatly expands the possibility of conducting comparative analyses thus giving additional insight into network evolution of land plants.



Protein-protein interaction




Heat-shock protein


Calvin-Benson-Bassham cycle


  1. Wiles AM, Doderer M, Ruan J, Gu TT, Ravi D, Blackman B, et al. Building and analyzing protein interactome networks by cross-species comparisons. BMC Sys Bio. 2010;4:36.

    Article  CAS  Google Scholar 

  2. Gandhi TKB, Zhong J, Mathivanan S, Karthick L, Chandrika KN, Mohan SS, et al. Analysis of the human protein interactome and comparison with yeast, worm and fly interaction datasets. Nature Gen. 2006;38:285–93.

    Article  CAS  Google Scholar 

  3. Giot L, Bader JS, Brouwer C, Chaudhuri A, Kuang B, Li Y, et al. A protein interaction map of Drosophila melanogaster. Science. 2003;302:1727–36.

    Article  CAS  PubMed  Google Scholar 

  4. Li S, Armstrong CM, Bertin N, Ge H, Milstein S, Boxem M, et al. A map of the interactome network of the metazoan C. elegans. Science. 2004;303:540–3.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  5. Uetz P, Giot L, Cagney G, Mansfield TA, Judson RS, Knight JR, et al. A comprehensive analysis of protein-protein interactions in Saccharomyces cerevisiae. Nature. 2000;403:623–7.

    Article  CAS  PubMed  Google Scholar 

  6. Yu H, Yildirim MA, Lemmens I, Venkateson K, Sahalie J, Hirozane-Kishikawa T, et al. High-quality binary protein interaction map of the yeast interactome network. Science. 2008;322:104–10.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  7. Yu QB, Li G, Wang G, Sun JC, Wang PC, Wang C, et al. Construction of a chloroplast protein interaction network and functional mining of photosynthetic proteins in Arabidopsis thaliana. Cell Res. 2008;18:1007–19.

    Article  CAS  PubMed  Google Scholar 

  8. Arabidopsis Interactome Mapping Consortium. Evidence for Network Evolution in an Arabidopsis Interactome Map. Science. 2011;333(6042):601–7.

    Article  CAS  PubMed Central  Google Scholar 

  9. DeBodt S, Proost S, Vandepoele K, Rouzé P, Van de Peer Y. Predicting protein-protein interactions in Arabidopsis thaliana through integrations of orthology, gene ontology and co-expression. BMC Genom. 2009;10:288.

    Article  CAS  Google Scholar 

  10. Geisler M, Fitzek E. A predicted interactome for coffee (Coffea canephora var robusta). Plant Mol Biol Biotechol. 2011;2(2):34–46.

    Google Scholar 

  11. Geisler-Lee J, O’Toole N, Ammar R, Provart NJ, Millar AH, Geisler M. A predicted interactome for Arabidopsis. Plant Physiol. 2007;145:317–29.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  12. Ho CL, Wu Y, Shen HB, Provart NJ, Geisler M. A predicted protein interactome for rice. Rice (N Y). 2012;5(1):15. doi: 10.1186/1939-8433-5-15.

    Google Scholar 

  13. von Mering C, Krause R, Cornell M, Oliver SG, Fields S, Bork P. Comparative assessment of large-scale data sets of protein-protein interactions. Nature. 2002;417:399–403.

    Article  CAS  Google Scholar 

  14. Walhout AJ, Sordella R, Lu X, Hartley JL, Temple GF, Brasch MA, et al. Protein interaction mapping in C. Elegans using proteins involved in vulval development. Science. 2000;287(5450):116–22.

    Article  CAS  PubMed  Google Scholar 

  15. Lang D, Zimmer AD, Rensing SA, Reski R. Exploring plant biodiversity: the Physcomitrella genome and beyond. Tr Plant Sci. 2008;13:542–9.

    Article  CAS  Google Scholar 

  16. Schaefer D. A new moss genetics: targeted mutagenesis in Physcomitrella patens. Annu Rev Plant Bio. 2002;53:477–501.

    Article  CAS  Google Scholar 

  17. Schuette S, Wood AJ, Geisler M, Geisler-Lee J, Ligrone R, Renzaglia KS. Novel localization of callose in the spores of Physcomitrella patens and phylogenomics of the callose synthase gene family. Anna Bot. 2009;103:749–56.

    Article  CAS  Google Scholar 

  18. Rensing SA, Lang D, Zimmer AD, Terry A, Salamov A, Shapiro H, et al. The Physcomitrella genome reveals evolutionary insights into the conquest of land by plants. Science. 2008;319:64–9.

    Article  CAS  PubMed  Google Scholar 

  19. Zimmer AD, Lang D, Buchta K, Rombauts S, Nishiyama T, Hasebe M, et al. Reannotation and extended community resources of the non-seed plant Physcomitrella patens provide insights into the evolution of plant gene structures and functions. BMC Genomics. 2013;14:498.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Aranda B, Achuthan P, Alam-Faruque Y, Armean I, Bridge A, Derow C, et al. The IntAct molecular interaction database. Nucl Acids Res. 2010;38(1):D525–31.

    Article  CAS  PubMed  Google Scholar 

  21. Bader GD, Hogue CW. An automated method for finding molecular complexes in large protein interaction networks. BMC Bioinform. 2003;4:2.

    Article  Google Scholar 

  22. Salwinski L, Miller CS, Smith AJ, Pettit FK, Bowie JU, Eisenberg D. The database of interacting proteins: 2004 update. Nucl Acids Res. 2004;32(1):D449–51.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Stark C, Breitkreutz BJ, Chatr-arvamontri A, Boucher L, Oughtred R, Livstone MS, et al. The BioGRID interaction database: 2001 update. Nucl Acids Res. 2011;39(1):D674–98.

    Google Scholar 

  24. Smedley D, Haider S, Ballester B, Holland R, London D, Thorisson G, et al. BioMart - biological queries made easy. BMC Genomics. 2009;10:22. doi: 10.1186/1471-2164-10-22.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  25. O’Brien KP, Remm M, Sonnhammer E. Inparanoid: a comprehensive database of eukaryotic orthologs. Nucl Acids Res. 2005;33(Database issue):D476–80.

    Article  PubMed  Google Scholar 

  26. Remm M, Storm C, Sonnhammer E. Automatic clustering of orthologs and in-paralogs from pairwise species comparisons. J Mol Bio. 2001;314:1041–52.

    Article  CAS  Google Scholar 

  27. Lang D, Eisinger J, Reski R, Rensing SA. Representation and high-quality annotation of the Physcomitrella patens transcriptome demonstrates a high proportion of proteins involved in metabolism among mosses. Plant Bio. 2005;7:228–37.

    Article  CAS  Google Scholar 

  28. Cline MS, Smoot M, Cerami E, Kuchinsky A, Landys N, Workman C, et al. Integration of biological networks and gene expression data using Cytoscape. Nature Protoc. 2007;2:2366–82.

    Article  CAS  Google Scholar 

  29. Saito R, Smoot ME, Ono K, Ruscheinski J, Wang PL, Lotia S, et al. A travel guide to Cytoscape plugins. Nat Methods. 2012;9(11):1069–76.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  30. Maere Heymans S, Kuiper M. BiNGO: a Cytoscape plugin to assess overrepresentation of Gene Ontology categories in biological networks. Bioinformatics. 2005;21(16):3449–88.

    Google Scholar 

  31. Sonnhammer ELL, Östlund G. InParanoid 8: orthology analysis between 273 proteomes, mostly eukaryotic. Nucl Acids Res. 2015;43(Database issue):D234–39.

    Article  PubMed  Google Scholar 

  32. Toufighi K, Brady SM, Austin R, Ly E, Provart NJ. The botany array resource: e-Northerns: expression angling, and promoter analyses. Plant J. 2005;43:153–63.

    Article  CAS  PubMed  Google Scholar 

  33. Jones AM, Xuan Y, Xu M, Wang R-S, Ho C-H, Lalonde S, et al. Border control—a membrane-linked interactome of Arabidopsis. Science. 2014;344:711–6.

    Article  CAS  PubMed  Google Scholar 

  34. Huang H, Jedynak BM, Bader JS. Where have all the interactions gone? Estimating the coverage of two-hybrid protein interaction maps. PLoS Comput Biol. 2007;3(11):e214. doi: 10.1371/journal.pcbi.0030214.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  35. He X, Zhang J. Why do hubs tend to be essential in protein networks? PLoS Genet. 2006;2(6):e88. doi:10.1371/journal.pgen.0020088.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  36. Batada NN, Reguly T, Breitkreutz A, Boucher L, Breitkreutz BJ, Hurst LD, et al. Stratus not altocumulus: a new view of the yeast protein interaction network. PLoS Bio. 2006;4:1720–31.

    Article  CAS  Google Scholar 

  37. Gavin AC, Aloy P, Grandi P, Krause R, Boesche M. Proteome survey reveals modularity of the yeast cell machinery. Nature. 2006;440:631–7.

    Article  CAS  PubMed  Google Scholar 

  38. Krause R, von Mering C, Jensen LJ, Bork P. Shared components of protein complexes – versatile building blocks of biochemical artefacts? Bioess. 2004;26:1333–43.

    Article  CAS  Google Scholar 

  39. Maslov S, Sneppen K. Specificity and stability in topology of protein networks. Science. 2002;296:910–3.

    Article  CAS  PubMed  Google Scholar 

  40. Moreno JJ. New aspects of the role of hydroxyeicosatertraenoic acids in cell growth and cancer development. Biochem Pharm. 2009;77:1–10.

    Article  CAS  PubMed  Google Scholar 

  41. Anterola A, Göbel C, Hornung E, Sellhorn G, Feussner I, Grimes H. Physcomitrella patens has lipoxygenases for both eicosanoid and octadecanoid pathways. Phytochem. 2009;70:45–57.

    Article  CAS  Google Scholar 

  42. Wolf YI, Karev G, Koonin EV. Scale-free networks in biology: new insights into the fundamentals of evolution? Bioess. 2002;24:105–9.

    Article  Google Scholar 

  43. Rensing SA, Ick J, Fawcett JA, Lang D, Zimmer A, Van de Peer Y, et al. An ancient genome duplication contributed to the abundance of metabolic genes in the moss Physcomitrella patens. BMC Evo Bio. 2007;7:130.

    Article  CAS  Google Scholar 

  44. Blanc G, Wolfe KH. Functional divergence of duplicated genes formed by polyploidy during Arabidopsis evolution. Plant Cell. 2004;16:1679–91.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  45. Brown KR, Jurisica I. Online predicted human interaction database. Bioinform. 2005;21:2076–82.

    Article  CAS  Google Scholar 

  46. Goh KI, Cusick ME, Valle D, Childs B, Vidal M, Barabasi AL. The human disease network. Proc Nat Acad Sci USA. 2007;104:8685–90.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  47. Oti M, Snel B, Huynen MA, Brunner HG. Predicting disease genes using protein-protein interactions. J Med Gen. 2006;43:691–8.

    Article  CAS  Google Scholar 

Download references


The authors extend thanks to Elisabeth Fitzek for valuable comments on and assistance with assembling earlier versions of this manuscript. This research supported by the NSF GK-12 Fellowship Program (#0638722) at Southern Illinois University.

Author information

Authors and Affiliations


Corresponding author

Correspondence to Matt Geisler.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

BP performed analysis of the data, created the table of conserved eukaryotic orthologs, created/edited figures, and contributed to preparing the manuscript. SS constructed the interactome, performed data analysis and contributed to preparing the manuscript. AC0 eveloped software to aid in interactome construction and database integration. DL aided with Physcomitrella patens genome sequences and gene annotation. MG supervised the project and aided in interactome construction and preparing/editing the manuscript. All authors read and approved the final manuscript.

Scott Schuette and Brian Piatkowski contributed equally to this work.

Additional files

Additional file 1:

Software package used in generating the interactome from databases.

Additional file 2: Table S1.

Moss Interolog Database.

Additional file 3: Table S2.

Moss Ortholog Database.

Additional file 4: Table S3.

Moss GO Annotation.

Additional file 5: Table S4.

List of interactions conserved in at least 5 eukaryotes.

Additional file 6: Table S5.

GDP-Mannose PPase interaction comparison.

Additional file 7: Table S6.

BinGO Analysis.

Additional file 8: Table S7.

CBB Subnetwork.

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit

The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Schuette, S., Piatkowski, B., Corley, A. et al. Predicted protein-protein interactions in the moss Physcomitrella patens: a new bioinformatic resource. BMC Bioinformatics 16, 89 (2015).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: