Skip to main content

Predicted functional interactome of Caenorhabditis elegans and a web tool for the functional interpretation of differentially expressed genes

Abstract

Background

The nematode worm, Caenorhabditis elegans, is a saprophytic species that has been emerging as a standard model organism since the early 1960s. This species is useful in numerous fields, including developmental biology, neurobiology, and ageing. A high-quality comprehensive molecular interaction network is needed to facilitate molecular mechanism studies in C. elegans.

Results

We present the predicted functional interactome of Caenorhabditis elegans (FIC), which integrates functional association data from 10 public databases to infer functional gene interactions on diverse functional perspectives. In this work, FIC includes 108,550 putative functional associations with balanced sensitivity and specificity, which are expected to cover 21.42% of all C. elegans protein interactions, and 29.25% of these associations may represent protein interactions. Based on FIC, we developed a gene set linkage analysis (GSLA) web tool to interpret potential functional impacts from a set of differentially expressed genes observed in transcriptome analyses.

Conclusion

We present the predicted C. elegans interactome database FIC, which is a high-quality database of predicted functional interactions among genes. The functional interactions in FIC serve as a good reference interactome for GSLA to annotate differentially expressed genes for their potential functional impacts. In a case study, the FIC/GSLA system shows more comprehensive and concise annotations compared to other widely used gene set annotation tools, including PANTHER and DAVID. FIC and its associated GSLA are available at the website http://worm.biomedtzc.cn.

Background

More than 50 years ago, Sydney Brenner selected Caenorhabditis elegans as a genetic model to study developmental biology and neurobiology because of its fully mapped genome and nervous system, rapid life cycle (3 days), and ease of laboratory cultivation [1,2,3,4]. The experimental strength and physiological similarities present in C. elegans and higher organisms (e.g., humans) have made it an important organism to explore a variety of subjects, including neuron fate specificity [5], axon guidance [6], and differentiation [7].

Recently, the development of omics technologies for C. elegans has become indispensable to systems biology, which can help to elucidate the mechanisms governing cellular physiology at the molecular level [8, 9]. Although the variety and complexity of omics data provide a global overview of potential mechanisms of physiological change, it also presents unparalleled challenges to describe the underlying design logic of physiological processes from molecular-level descriptions.

In fact, the existing approaches mostly rely on enrichment analysis to obtain high-level biological insights from the observed set of differentially expressed genes (SDEG). These enrichment-based methods evaluate whether an SDEG is enriched or clustered in a defined biological process. To date, many gene set annotation tools based on enrichment analysis have been developed to analyse SDEG, and they show significantly different expression in two or more physiological statuses. Some tools are widely used, including PANTHER [10], KEGG [11], and DAVID [12].

In reality, the observed SDEG is summarized into established biological concepts by the above strategies, which are successful in many cases. However, in practical use, enrichment-based methods frequently report that no annotation term is enriched or only report conceptually general terms (such as GO: 0005634, nucleus), where no established biological concepts can be used to accurately describe the differentially expressed gene set. Clearly, these results provide little assistance for researchers attempting to elucidate molecular mechanisms.

Alternatively, while there are no established concepts to accurately describe the observed SDEG, we still may utilize established biological concepts to interpret the functional impacts of SDEG. For example, observed SDEG may lead collectively to GO: 0051704 (defence response to fungus), even when the SDEGs themselves are not enriched in these terms (for details, please see the Discussion section). In this study, the gene set linkage analysis (GSLA) tool was developed to describe the potential functional impacts of the observed SDEG, especially in cases in which no established concepts or no suitable concepts can describe these changes. If the SDEG has strong functional associations with genes in an established biological process, the observed SDEG will be expected to interfere with this biological function. This strategy is called GSLA and has been successfully used in human and Arabidopsis transcriptome interpretation [13, 14]. The main contribution to the successful interpretations of functional impacts is the high-quality functional association networks in these two species [13, 15]. Worm researchers have made efforts to construct molecular interaction databases, which include WormBase (2390 proteins, 6343 interactions) [16], WormNet (16,122 proteins, 760,116 interactions) [17], mentha (4774 proteins, 12,136 interactions) [18], MIST (4929 proteins, 626,262 interactions) [19], STRING (12,050 proteins, 3,300,700 interactions) [20], and ComPPI (4958 proteins, 13,659 interactions) [21]. These databases have greatly facilitated C. elegans research, but they do not support the gene set linkage analysis (GSLA) algorithm well in our evaluation.

Therefore, in this study, we constructed a high-quality reference functional gene association network for C. elegans, which provides an alternative resource for high-reliability functional gene associations to facilitate mechanism exploration. Additionally, this set of functional gene associations supports the GSLA interpretation of the collective functional impacts of SDEGs in worm. The predicted functional interactome of Caenorhabditis elegans (FIC) integrates six types of functional association data from 10 public databases with timestamps before 2018. To evaluate the accuracy of inferred functional associations in FIC, we used experimentally confirmed protein-protein interactions recently reported after 2018. The current version of FIC includes 108,550 functional gene associations, which are expected to cover 21.42% of protein-protein interactions in C. elegans, and 29.25% of functional gene interactions may represent protein interactions. For users to query functional associations of their genes of interest in C. elegans, we provide an FIC web interface that is easy to operate. We also provide a GSLA web tool so that users can interpret the potential functional impacts of the observed set of differently expressed genes. A case study is also provided to illustrate the way to use the FIC/GSLA system.

Methods

Evidence data

In this work, we selected six types of evidence that suggest functional associations between genes from 7 public databases to build the interaction prediction model. These data were collected before 2018 and include 18,947 expression profiles (Coxpresdb), 44,300 gene annotations (GOC), 35,195 domain interactions (IDDI and Pfam), 14,457 subcellular gene localizations (Compartments), 14,866 phylogenetic profiles (DIOPT), and inparalogue/orthologue relationships between 5137 worm proteins and proteins from A. thaliana, D. melanogaster, H. sapiens, M. musculus, R. norvegicus, S. cerevisiae and S. pombe to compute interologues. Thirty-six features belonging to six categories were computed based on these evidence data, each suggesting a certain type of functional association (Supplementary Table 2) [22, 23]. Detailed methods and equations can be found in website help (the indirect evidence section).

Protein-protein interactions were considered to be evidence of strong functional associations. In this study, we attempt to predict functional associations that are as strong as protein interactions [13]. We collected a total of 26,367 experimentally reported protein-protein interactions of C. elegans from three public databases, including WormBase [16], BioGRID [24], and IntAct [25] (Additional file 1: Table S1). Among 26,367 interactions, we retained 5606 high-quality protein-protein interactions to ensure that these interactions are experimentally confirmed, rather than predicted. The interactions that were reported in less than two independent studies or reported only in high-throughput experiments were removed (Additional file 1: Table S1). To obtain the uniform gene ID for the functional association prediction, UniProt [26] and BioMart [27] software were used to convert different gene IDs to WormBase ID (Fig. 1).

Fig. 1
figure1

Workflow for inferencing functional interactions between C. elegans genes. High-quality experimentally reported protein interactions were integrated from three databases and were used as positive examples. Six types of functional association evidence from 10 databases were collected to infer putative functional interactions. A total of 10 high-quality feature values were selected from 36 feature values that characterize this evidence with different mathematical representations. Random gene pairs after removing positive examples were used as negative examples. The number of negative examples was 100 times that of the positive examples

Computation and assessment of feature values

From six types of functional association evidence, we used 36 mathematical characterizations to compute feature values, which include 3 shared annotation features, 2 co-expression features, 4 subcellular co-localization features, 23 domain interaction features, 3 phylogenetic profile features, and 1 homologous interaction (Additional file 3: Table S2).

To evaluate whether the above feature values can suggest functional associations, the area under the curve (AUC) of the receiver operating characteristic (ROC) test was employed. Each feature value produced a series of sensitivities and specificities by applying different cut-offs when predicting the protein interactions. According to different cut-offs, the ROC curve was plotted with sensitivities and specificities (X-axis, 1-specificity; Y-axis, sensitivity). Assessment of feature values was performed using the training dataset (data before 2018). If the AUC of a feature value was higher than 0.6, a strong functional association was present. A total of 10 features were eventually selected to predict the functional associations (Additional file 4: Table S3 and Additional file 2: Fig. S1).

Website construction

Our developed FIC database runs on LNMP that integrated the running environment, including Linux, Nginx, MySQL, and PHP. Data storage, maintenance, and operation are supported by the MySQL database. The interaction interface is developed with the Laravel framework based on PHP. The front-end of FIC is a Vue.js-based single page application (SPA). Vue.js is an open-source JavaScript framework for interface creation and a web application framework compatible with SPA. The functional association network is visualized with Cytoscape [28].

Microarray data analysis

From the GEO database [29, 30], we retrieved microarray data GSE97678 [31]. In the original article, 35 genes were upregulated between 2- and 8-fold on the E. coli HT115 diet, while 22 genes were upregulated between 2- and 20-fold in the E. coli OP50 diet. Three biological replicates were performed by the authors. These expression profiles were re-analysed based on the online GEO2R tool [32,33,34] using default parameters. The top transcriptionally changed genes were selected by adj. P. Val (P-value after adjustment for multiple testing). The GEO2R tool adjusts the P-values to correct for false positive results. We chose the default Benjamini & Hochberg false discovery rate method to adjust multiple tests.

Results

Data integration for functional association prediction in C. elegans

Six types of evidence for functional association inference between C. elegans genes were collected from seven databases, including Coxpresdb [35], Gene Ontology Consortium (GOC) [36], Compartments [37], IDDI [38], Pfam [39], DIOPT [40] and Inparanoid [41] (Fig. 1). Thirty-six feature values from six types of evidence were used to measure the strength of functional associations (Additional file 3: Table S2).

Not all 36 of these features were suitable to separate protein interactions from random gene pairs. To decrease the noise-to-signal ratio in the following functional associations prediction step, we only kept those features that showed a strong correlation to functional associations. The AUC of the receiver operating characteristic (ROC) curve was used to measure the capability of a feature to indicate protein interactions. In this study, we selected ten features with AUC higher than 0.6 for the subsequent inference of functional gene associations (Additional file 4: Table S3 and Additional file 2: Fig. S1).

In addition, protein-protein interactions reported in experimental studies of C. elegans were collected from three databases, including WormBase [16], BioGRID [24], and IntAct [25] (Fig. 1 and Additional file 1: Table S1). We removed protein-protein interactions based on the supporting evidence provided in each database (Fig. 1), only retaining the experimentally confirmed high-quality protein-protein interactions, which were used as positive examples in prediction model training (Additional file 1: Table S1).

Functional gene association prediction

We used the libSVM package to train and predict functional associations [42, 43] (Fig. 1). Specifically, 5606 high-quality experimentally confirmed protein-protein interactions that were published before 2018 were used as positive examples, which represent examples of strong functional associations between C. elegans genes. Randomly generated gene pairs that do not overlap with positive examples served as negative examples. We considered that two random genes may have functional associations, although the probability is low. To reduce the false positive rate in the negative examples, we set the positive-to-negative ratio in the training dataset to be 1:100 such that only a notably small fraction of gene pairs has functional associations.

We utilized the soft-margin Gaussian kernel SVM algorithm to train the prediction model. The parameters σ (kernel width) and C (soft margin) targeted an optimal harmonic mean of sensitivity and specificity that were optimized with a 5-fold cross-validation. Based on the optimized σ and C, all training data were used to train the prediction model, which was validated with an external validation dataset consisting of protein-protein interactions that were published after December 31, 2017, and randomly generated negative examples. After validation, this model showed a sensitivity of 21.42% and a specificity of 99.95%. In this study, we also evaluated how well the predicted interaction in WormNet, MIST, and STRING covered these new interactions. The results are shown in Supplementary Table S4.

Applying this model to all gene pairs of C. elegans produced 101,727 inferred functional associations. These inferred functional interactions together with the 6823 known protein interactions make the FIC dataset, which consists of 108,550 interactions. A total of 108,550 putative functional interactions were inferred when applying this model to all C. elegans gene pairs. Among the protein-protein interactions, we further estimated the proportion covered by the predicted functional interactome using the following equation:

$$ {N}_{interactome}\times Sensitivity+\left({N}_{all- pairs}-{N}_{interactome}\right)\times \left(1- specificity\right)={N}_{\mathrm{predict}} $$

where Ninteractome is the expected number of all protein-protein interactions in C. elegans, Nall − pairs is the number of all gene pairs in C. elegans, and Npredict is the number of predicted gene associations. The sensitivity and specificity are the accuracy measures produced when the prediction model was validated with newly published protein interactions. Solving this equation yields an estimated size of 1.42 × 105 for the C. elegans protein interactome. Based on the estimated interactome size (1.42 × 105) and the estimated sensitivity (21.42%, the conservative estimation from the training stage sensitivity 21.69% and the evaluation stage sensitivity 21.42%), the predicted interactions in FIC are expected to include 29,755 protein interactions. Therefore, 29.25% of the FIC functional interactions (29,755 out of 101,727) are expected to represent protein interactions.

Assessment of functional gene association network

We evaluated the quality of the functional gene association interactome FIC by its ability to group functionally associated genes together. A gene’s function prediction accuracy based on its network neighbours may be used to measure the quality of our predicted interactome FIC. To this end, the quality of the inferred functional associations was compared with six other public C. elegans interactomes, including WormBase [16], WormNet [17], mentha [18], MIST [19], STRING [20], and ComPPI [21]. In this study, the PATHER term enrichment tool [10] was used to measure the accuracies to predict the new GO biological process annotations on each interactome.

The collection date of data for the inference of FIC gene associations was before 2018 (December 31, 2017). We collected 13,043 genes from GO [44, 45] with new annotations dated up to August 1, 2018. These genes contained a total of 117,848 annotations, 9108 annotations of which were newly added after 2018. We evaluated the prediction performance of gene function based on these genes and their annotations.

The overall prediction accuracy of new annotations across seven interactomes was compared by precision-recall curve. In this instance, precision means whether the annotations reported by PANTHER are consistent with the known annotations (all 117,848 annotations), while recall means the proportion of annotations reported by PANTHER that covered the newly added 9108 annotations. The number of PANTHER reported annotations is correlated with the cut-off on the significant value. Setting a higher cut-off value results in more reported annotations and a higher recall but with a higher false positive rate. In contrast, setting a lower cut-off value results in fewer reported annotations and lower recall but with higher precision. Therefore, the advantage of the precision-recall curve is that it shows the rates of precision and recall on different cut-offs. Independent of the selection cut-offs, the precision-recall curve can provide a more comprehensive view of the quality of the interactome.

Figure 2 shows that the curve of our predicted FIC interactome is above six other interactomes, indicating its superior quality in grouping functionally associated genes together. Compared to others, only the curve of FIC reaches the high-recall region, and it still maintains the highest precision. Although WormBase, mentha, MIST, and ComPPI have similar high-precision regions, none of these curves reached high-recall regions. Alternatively, STRING and WormNet reached the high-recall region, but their precision did not increase considerably and always stayed in the low-recall region. This observation suggested that a high proportion of STRING and WormNet interactions were weak functional gene associations, which may raise the false positive rates during function prediction. In general, FIC showed balanced coverage and accuracy, and its quality exceeds those of other compared interactomes.

Fig. 2
figure2

Assessment of the capabilities of seven interactomes to group functionally associated genes together. The precision-recall curves of gene function prediction using different interactomes are illustrated. Precision estimates the proportion of correct annotations identified by an interactome. Recall estimates the proportion of new annotations that is identified by an interactome

Web interface of FIC/GSLA

We developed a user-friendly interface for FIC, which provides two search modes, a single gene search mode and a multiple gene search mode (Fig. 3a). Both gene name and WormBase ID are acceptable means of searching for a gene. The single gene search option reports putative functional gene-associated interactions involving the query gene, while the multiple search option reports the functional associations between query genes. The functionally associated interactions reported by FIC are provided in a tabular form (Fig. 3b). A graphical view of these reported interactions is also provided at the right side of the query interface. Users can check for feature values for the interaction prediction in our model if they click the edges in the graphical view of their functional associations. Moreover, users can click the nodes in the graphical view to obtain gene information that provides more detailed information about a gene. A full dump of the FIC database is available for download. We also provide a help section on the FIC/GSLA website with more details for users.

Fig. 3
figure3

Interface of FIC and GSLA. a Two search options in FIC. b Search result page. A right click on the edge and node will show interaction details. c Interface of GSLA. d Results of a GSLA analysis job

The GSLA web tool was first developed for the Predicted Arabidopsis Interactome Resource (PAIR) to interpret potential functional impacts of observed SDEG in Arabidopsis [14]. Two testing hypotheses (Q1 and Q2) are used in GSLA to ensure that the reported associations between two C. elegans gene sets are significant (Fig. 4). Q1 examines whether the gene association density of inter-gene-sets between two functionally associated gene sets is higher than the background random gene sets. Q2 examines whether the observed high-density gene associations between two gene sets can only be observed in the biologically correct functional gene association network. In other words, the density observed in FIC is higher than the densities observed in random gene association networks consisting of the same genes, with each gene having the same number of neighbours. From a biological perspective, Q1 tests the strength of a functional association between two gene sets. Q2 confirms that the observed strong inter-gene-set functional interaction is the result of a biologically correct functional interactome (i.e., our knowledge of the molecular mechanisms), rather than the result of the gene set compositions. As is well-known, in an interactome, some hub genes have substantially more neighbours than others. Gene sets that include a number of hubs may easily have many inter-gene-set functional interactions with other gene sets. Therefore, Q2 is used to control this confounding factor of gene set composition. Q1 and Q2 are different but complementary tests, both of which can increase the sensitivity and specificity of GSLA and can ensure the biological significance of the functional associations detected between gene sets. The default significance cut-offs for GSLA to report a gene set interaction are density > 0.01 (Q1) and p < 0.001 (Q2).

Fig. 4
figure4

Two hypothesis tests that GSLA used to identify significant functional associations between two gene sets that are biologically meaningful. Q1 tests whether the density of functional associations between two biologically meaningful gene sets is higher than random gene pairs, while Q2 tests whether the strong functional associations observed between two gene sets can only be observed from the biologically correct network, rather than any random interactomes

The link for GSLA online service is provided in the web interface of FIC, the functional gene association network of which is used to interpret the functional impacts of observed SDEG in C. elegans. The main interface is shown in Fig. 3c. In this instance, we provide five types of C. elegans gene IDs for users to query SDEG, including WormBase ID, gene name, UniProt ID, Ensembl gene ID, Ensembl protein ID, and NCBI Entrez ID. We suggest providing SDEG directly in WormBase IDs because WormBase ID is only recognized by the internal server. Therefore, all of the submitted IDs will be automatically mapped to WormBase ID before further computation (Fig. 3d). To avoid loss of information during ID mapping processes, submitting SDEG lists to WormBase ID is better. The criteria of GSLA to report significant functional gene associations (Q1 and Q2 tests, described above) can be optimized by users (Fig. 3c). Before submission, an email address is needed for the results. Utilizing the top 50–200 changed genes is recommended when querying the observed SDEG to obtain optimal functional impact interactions. The analysis parameters are provided at the top ten lines in the result file (Fig. 3e). Below these parameters, a table is presented to show the functionally associated biological processes, functional associations between genes in reported biological processes, and the genes in the query SDEG.

Using the FIC/GSLA system to re-analyse the E. coli HT115 diet dataset

E. coli OP50 and E. coli HT115 are two standard laboratory diets for C. elegans [31]. E. coli OP50 is the most commonly used food, while E. coli HT115 is typically reared in RNAi-mediated gene knockdown experiments. However, some studies have discovered that E. coli OP50, as the most common lab diet for C. elegans, will cause a mild, chronic vitamin B12 deficiency [31, 46, 47]. Therefore, Revtovich et al. performed a variety of assays and confirmed a diet of E. coli OP50 results in vitamin B12 deficiency, which disrupts mitochondrial homeostasis and decreases the host resistance [31]. Feeding C. elegans E. coli HT115 or overexpression of the B12 transporter improved mitochondrial homeostasis and increased resistance. To further explore the molecular mechanisms underlying the diet-induced difference in stress resistance, these researchers performed a genetic analysis to map this phenotype to the methylmalonyl/succinyl-CoA breakdown pathway, where vitamin B12 serves as a cofactor for MMCM-1/MUT. The authors compared the transcriptomic profiles of C. elegans with different diets, E. coli OP50 and E. coli HT115 (GEO database, GSE97678). The results of the microarray analysis showed that the number of differentially regulated genes was relatively small. Only 22 genes were upregulated in E. coli OP50 (between 2- and 20-fold). Interestingly, more than half (12 genes) of these 22 upregulated gene-encoded proteins were localized on mitochondria.

In this study, we used the FIC/GSLA system to reanalyse the functional impacts of different diet-induced gene changes. To evaluate whether the prediction tools can obtain biological insights from these changed genes of the microarray data (GEO database, GSE97678) [31], we performed FIC/GSLA, DAVID [12], and GO enrichment analysis [44, 45]. As shown in Fig. 5, GO enrichment analysis reported immune response- and defence response-related pathways (Additional file 6: Table S5). Defence response-related GO terms were consistent with the original publication; however, the authors have clearly documented that diet-mediated sensitivity is independent of innate immunity [31]. The term clustering technology-based tool DAVID reported a total of 43 terms in 7 clusters. Among them, there are 22 GO terms. The top 5 GO terms (ranked by PValue) included the biological processes of mitochondrial and acyl-CoA dehydrogenase-related metabolism that are consistent with those reported in the original article (Additional file 7: Table S6). For comparison, we also used FIC/GSLA to annotate the SDEGs (Additional file 8: Table S7). Eighteen terms covered defence response, mitochondria, and propionyl-CoA pathways (Fig. 5), which are known to be the “functional impact” of the E. coli OP50 diet. Moreover, both the DAVID and FIC/GSLA systems reported fatty acid metabolic processes that are consistent with the discoveries by Brooks et al. [48]. In addition, GSLA also found SUMOylation-related biological processes (Fig. 5). In an independent study, Benedetti et al. indeed observed that ubiquitin-like protein 5 (UBL-5) positively regulates chaperone gene expression in response to mitochondrial unfolded proteins [49]. In this case study, the interpretation provided by FIC/GSLA is broader and more accurate, providing new insights for experimental researchers to explore molecular mechanisms, while other widely used enrichment-based tools did not provide similar insights.

Fig. 5
figure5

Functional interpretations produced by FIC/GSLA. Compared to GO enrichment analysis and DAVID, the annotations produced by GSLA are more comprehensive and more accurate

Discussion

Before this work, numerous efforts have been made to build reference interactomes for C. elegans. An accurate and comprehensive reference interactome may facilitate the interpretation of gene transcriptional changes to higher-level biological process changes. To date, many C. elegans interactome databases have emerged. Some of these interactome databases, e.g., WormBase [16], BioGRID [24], and IntAct [25], collect experimentally reported molecular interactions. Others, such as WormNet [17] and STRING [20], provide predicted molecular interactions. In general, although experimentally reported molecular interactions are considered to be more accurate than predicted molecular interactions, the number of molecular interactions reported by experiments is too low. Based on our estimated size of the C. elegans interactome, which is 1.42 × 105, the WormBase interactome database that provides experimental molecular interactions shows the highest coverage with sizes up to 10.35%, which represents 17.15% of the C. elegans protein interactome. In fact, due to the false positive experimental molecular interactions, the coverage of the WormBase database will be lower than 10.35%. Such a low coverage may provide limited support for researchers to explain molecular mechanisms during their research. In contrast to the databases with experimental interactions, the databases with predicted molecular interactions may present a high coverage of the true C. elegans protein interactomes. STRING is a widely used interactome that contains a large number of predicted interactions. In C. elegans, STRING presents 31.89% of predicted interactions, which cover 31.89% of the protein interactome. However, these predicted interactions often have a very high false positive rate, and only 3.46% of them are expected to represent protein interactions. Compared to the currently widely used interactome databases, our newly developed FIC interactome shows better performance with balanced coverage and reliability (21.42% coverage and 29.25% reliability if evaluated as a protein interaction network). In conclusion, FIC is a high-quality reference protein interaction network used to analyse functional gene interactions.

Our high-quality FIC enables GSLA for the interpretation of the observed SDEG in C. elegans. FIC has both a high-precision and high-coverage functional interactome, which helps GSLA to report significant functional associations between gene sets. The assessment of GSLA evaluates the density of functional gene associations between individual genes in two gene sets. The interactome of FIC with balanced accuracy will facilitate the successful application of this strategy, while previous interactomes did not satisfy this requirement. We assessed the GSLA interpretation of FIC compared to other interactomes, and FIC showed the best performance (data not shown). The same observation was made when we previously developed high-quality functional interactomes for humans and Arabidopsis.

As mentioned in the introduction, the FIC/GSLA system can interpret the potential functional impacts of the observed SDEG of C. elegans. Therefore, the FIC/GSLA system extends the availability of current enrichment-based tools to summarize SDEG into known biological processes. In some cases, even when no established biological concept can accurately explain the observed SDEG, the FIC/GSLA system may still be able to interpret the observed gene transcriptional changes and to connect the related physiologies. Moreover, the functional association resource provided in FIC is a useful reference for researchers to elucidate the molecular mechanisms of their genes of interest.

Conclusions

Thus, predicted FIC is a reliable and high-quality resource for querying functional associations between genes of C. elegans, as it can help researchers understand the molecular mechanisms of certain genes. Based on the functional gene association network of FIC, GSLA was developed to facilitate the interpretation of potential functional impacts from the observed differentially expressed genes, especially when no established biological concept is available to describe the observed SDEG. A case study of the FIC/GSLA system shows that the reported annotations are more comprehensive and concise annotations compared to the other widely used annotation tools, including PANTHER and DAVID.

Availability of data and materials

The predicted interactome FIC and its associated GSLA web tool are available at http://worm.biomedtzc.cn.

References

  1. 1.

    Corsi AK, Wightman B, Chalfie M. A transparent window into biology: a primer on Caenorhabditis elegans. Genetics. 2015;200:387–407.

  2. 2.

    Leung MCK, Williams PL, Benedetto A, Au C, Helmcke KJ, Aschner M, et al. Caenorhabditis elegans: an emerging model in biomedical and environmental toxicology. Toxicol Sci. 2008;106:5–28.

  3. 3.

    Hillier LW, Coulson A, Murray JI, Bao Z, Sulston JE, Waterston RH. Genomics in C. elegans: so many genes, such a little worm. Genome Res. 2005;15:1651–60.

  4. 4.

    Brenner S. The genetics of behaviour. Br Med Bull. 1973;29:269–71.

    CAS  Article  Google Scholar 

  5. 5.

    Hobert O. Specification of the nervous system. WormBook. 2005:1–19.

  6. 6.

    Quinn CC, Wadsworth WG. Axon guidance: ephrins at WRK on the midline. Curr Biol. 2006;16:R954–5.

    CAS  Article  Google Scholar 

  7. 7.

    Chisholm AD, Jin Y. Neuronal differentiation in C. elegans. Curr Opin Cell Biol. 2005;17:682–9.

  8. 8.

    Van Assche R, Broeckx V, Boonen K, Maes E, De Haes W, Schoofs L, et al. Integrating -Omics: systems biology as explored through C. elegans research. J Mol Biol. 2015;427:3441–51.

  9. 9.

    Hastings J, Mains A, Virk B, Rodriguez N, Murdoch S, Pearce J, et al. Multi-Omics and Genome-Scale Modeling Reveal a Metabolic Shift During C elegans Aging. Front Mol Biosci. 2019;6:2.

    CAS  Article  Google Scholar 

  10. 10.

    Mi H, Huang X, Muruganujan A, Tang H, Mills C, Kang D, et al. PANTHER version 11: expanded annotation data from gene ontology and Reactome pathways, and data analysis tool enhancements. Nucleic Acids Res. 2017;45:D183–9.

    CAS  Article  Google Scholar 

  11. 11.

    Kanehisa M, Sato Y, Furumichi M, Morishima K, Tanabe M. New approach for understanding genome variations in KEGG. Nucleic Acids Res. 2019;47:D590–5.

    CAS  Article  Google Scholar 

  12. 12.

    Huang DW, Sherman BT, Lempicki RA. Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources. Nat Protoc. 2009;4:44–57.

    CAS  Article  Google Scholar 

  13. 13.

    Zhou X, Chen P, Wei Q, Shen X, Chen X. Human interactome resource and gene set linkage analysis for the functional interpretation of biologically meaningful gene sets. Bioinformatics. 2013;29:2024–31.

    CAS  Article  Google Scholar 

  14. 14.

    Yao H, Wang X, Chen P, Hai L, Jin K, Yao L, et al. Predicted Arabidopsis Interactome Resource and Gene Set Linkage Analysis: A Transcriptomic Analysis Resource. Plant Physiol. 2018;2018:00144.

    Google Scholar 

  15. 15.

    Lin M, Shen X, Chen X. PAIR: the predicted Arabidopsis interactome resource. Nucleic Acids Res. 2011;39:D1134–40.

    CAS  Article  Google Scholar 

  16. 16.

    Lee RYN, Howe KL, Harris TW, Arnaboldi V, Cain S, Chan J, et al. WormBase 2017: molting into a new stage. Nucleic Acids Res. 2018;46:D869–74.

    CAS  Article  Google Scholar 

  17. 17.

    Cho A, Shin J, Hwang S, Kim C, Shim H, Kim H, et al. WormNet v3: a network-assisted hypothesis-generating server for Caenorhabditis elegans. Nucleic Acids Res. 2014;42:W76–82.

  18. 18.

    Calderone A, Castagnoli L, Cesareni G. mentha: a resource for browsing integrated protein-interaction networks. Nat Methods. 2013;10:690–1.

    CAS  Article  Google Scholar 

  19. 19.

    Hu Y, Vinayagam A, Nand A, Comjean A, Chung V, Hao T, et al. Molecular interaction search tool (MIST): an integrated resource for mining gene and protein interaction data. Nucleic Acids Res. 2018;46:D567–74.

    CAS  Article  Google Scholar 

  20. 20.

    Szklarczyk D, Franceschini A, Wyder S, Forslund K, Heller D, Huerta-Cepas J, et al. STRING v10: protein-protein interaction networks, integrated over the tree of life. Nucleic Acids Res. 2015;43:D447–52.

    CAS  Article  Google Scholar 

  21. 21.

    Veres DV, Gyurkó DM, Thaler B, Szalay KZ, Fazekas D, Korcsmáros T, et al. ComPPI: a cellular compartment-specific database for protein-protein interaction network analysis. Nucleic Acids Res. 2015;43:D485–93.

    CAS  Article  Google Scholar 

  22. 22.

    Rhodes DR, Tomlins SA, Varambally S, Mahavisno V, Barrette T, Kalyana-Sundaram S, et al. Probabilistic model of the human protein-protein interaction network. Nat Biotechnol. 2005;23:951–9.

    CAS  Article  Google Scholar 

  23. 23.

    Shoemaker BA, Panchenko AR. Deciphering protein-protein interactions. Part II Computational methods to predict protein and domain interaction partners. PLoS Comput Biol. 2007;3:–e43.

  24. 24.

    Oughtred R, Stark C, Breitkreutz B-J, Rust J, Boucher L, Chang C, et al. The BioGRID interaction database: 2019 update. Nucleic Acids Res. 2019;47:D529–41.

    CAS  Article  Google Scholar 

  25. 25.

    Orchard S, Ammari M, Aranda B, Breuza L, Briganti L, Broackes-Carter F, et al. The MIntAct project--IntAct as a common curation platform for 11 molecular interaction databases. Nucleic Acids Res. 2014;42:D358–63.

    CAS  Article  Google Scholar 

  26. 26.

    UniProt Consortium. UniProt: a worldwide hub of protein knowledge. Nucleic Acids Res. 2019;47:D506–15.

    Article  Google Scholar 

  27. 27.

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

    Article  Google Scholar 

  28. 28.

    Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.

    CAS  Article  Google Scholar 

  29. 29.

    Edgar R, Domrachev M, Lash AE. Gene expression omnibus: NCBI gene expression and hybridization array data repository. Nucleic Acids Res. 2002;30:207–10.

    CAS  Article  Google Scholar 

  30. 30.

    Barrett T, Wilhite SE, Ledoux P, Evangelista C, Kim IF, Tomashevsky M, et al. NCBI GEO: archive for functional genomics data sets--update. Nucleic Acids Res. 2013;41:D991–5.

    CAS  Article  Google Scholar 

  31. 31.

    Revtovich AV, Lee R, Kirienko NV. Interplay between mitochondria and diet mediates pathogen and stress resistance in Caenorhabditis elegans. Garsin DA, editor. PLoS Genet. 2019;15:e1008011.

  32. 32.

    Smyth GK. Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004;3:Article3.

    Article  Google Scholar 

  33. 33.

    Smyth GK. limma: Linear Models for Microarray Data. In: Gentleman R, Carey VJ, Huber W, Irizarry RA, Dudoit S, editors. Bioinformatics and Computational Biology Solutions Using R and Bioconductor [Internet]. New York: Springer-Verlag; 2005. p. 397–420. Available from: http://link.springer.com/10.1007/0-387-29362-0_23, [cited 2020 Mar 23].

    Google Scholar 

  34. 34.

    Davis S, Meltzer PS. GEOquery: a bridge between the gene expression omnibus (GEO) and BioConductor. Bioinformatics. 2007;23:1846–7.

    Article  Google Scholar 

  35. 35.

    Obayashi T, Kagaya Y, Aoki Y, Tadaka S, Kinoshita K. COXPRESdb v7: a gene coexpression database for 11 animal species supported by 23 coexpression platforms for technical evaluation and evolutionary inference. Nucleic Acids Res. 2019;47:D55–62.

    CAS  Article  Google Scholar 

  36. 36.

    Gene Ontology Consortium. Gene ontology consortium: going forward. Nucleic Acids Res. 2015;43:D1049–56.

    Article  Google Scholar 

  37. 37.

    Binder JX, Pletscher-Frankild S, Tsafou K, Stolte C, O’Donoghue SI, Schneider R, et al. COMPARTMENTS: unification and visualization of protein subcellular localization evidence. Database (Oxford). 2014;2014:bau012.

    Article  Google Scholar 

  38. 38.

    Kim Y, Min B, Yi G-S. IDDI: integrated domain-domain interaction and protein interaction analysis system. Proteome Sci. 2012;10(Suppl 1):S9.

    Article  Google Scholar 

  39. 39.

    El-Gebali S, Mistry J, Bateman A, Eddy SR, Luciani A, Potter SC, et al. The Pfam protein families database in 2019. Nucleic Acids Res. 2019;47:D427–32.

    CAS  Article  Google Scholar 

  40. 40.

    Hu Y, Flockhart I, Vinayagam A, Bergwitz C, Berger B, Perrimon N, et al. An integrative approach to ortholog prediction for disease-focused and other functional studies. BMC Bioinformatics. 2011;12:357.

    Article  Google Scholar 

  41. 41.

    O’Brien KP, Remm M, Sonnhammer ELL. Inparanoid: a comprehensive database of eukaryotic orthologs. Nucleic Acids Res. 2005;33:D476–80.

    Article  Google Scholar 

  42. 42.

    Chang C-C, Lin C-J. LIBSVM: a library for support vector machines. ACM Trans Intell Syst Technol. 2011;2:1–27.

    Article  Google Scholar 

  43. 43.

    Winters-Hilt S, Yelundur A, McChesney C, Landry M. Support vector machine implementations for classification & clustering. BMC Bioinformatics. 2006;7(Suppl 2):S4.

    Article  Google Scholar 

  44. 44.

    Ashburner M, Ball CA, Blake JA, Botstein D, Butler H, Cherry JM, et al. Gene ontology: tool for the unification of biology. The Gene Ontology Consortium. Nat Genet. 2000;25:25–9.

    CAS  Article  Google Scholar 

  45. 45.

    The Gene Ontology Consortium. The gene ontology resource: 20 years and still GOing strong. Nucleic Acids Res. 2019;47:D330–8.

    Article  Google Scholar 

  46. 46.

    MacNeil LT, Watson E, Arda HE, Zhu LJ, AJM W. Diet-induced developmental acceleration independent of TOR and insulin in C elegans. Cell. 2013;153:240–52.

    CAS  Article  Google Scholar 

  47. 47.

    Watson E, Yilmaz LS, Walhout AJM. Understanding metabolic regulation at a systems level: metabolite sensing, mathematical predictions, and model organisms. Annu Rev Genet. 2015;49:553–75.

    CAS  Article  Google Scholar 

  48. 48.

    Brooks KK, Liang B. Watts JL. The influence of bacterial diet on fat storage in C elegans PLoS ONE. 2009;4:e7545.

    Google Scholar 

  49. 49.

    Benedetti C, Haynes CM, Yang Y, Harding HP, Ron D. Ubiquitin-like protein 5 positively regulates chaperone gene expression in the mitochondrial unfolded protein response. Genetics. 2006;174:229–39.

    CAS  Article  Google Scholar 

Download references

Acknowledgements

We thank Prof. Xiao-Ming Zhao for offering useful discussion and Taizhou Bigdata AI Research Center for providing computing resources.

Funding

This work was financially supported by the National Key Research and Development Program (2019YFA0905401), the National Natural Science Foundation of China (81830073 and 31571356), and the Breeding program of Taizhou University (grant no. 2018PY010).

Author information

Affiliations

Authors

Contributions

Peng-Cheng Chen and Ruan Li developed the computational workflow and analysed the data with the help of Heng Yao. Jie Jin developed the web database and tools. Yu-Tian Tao wrote the manuscript. Hai-Bo Zhang built the interaction model. Xiao-Bao Ding evaluated the functional association network of FIC. Wen-Ping Guo performed the interaction prediction. Qiao-Lei Yang reanalysed the microarray dataset. Heng Yao provided initial conceptualizations of the FIC/GSLA system. Xin Chen designed and coordinated the project and, together with Yu-Tian Tao, wrote the first draft of the manuscript. All authors reviewed and edited the final manuscript.

Corresponding author

Correspondence to Xin Chen.

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 interests.

Additional information

Publisher’s Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary information

Additional file 1: Table S1.

Number of protein-protein interactions and participating unique proteins from three databases.

Additional file 2: Table S2.

Functional association evidence and computing methods.

Additional file 3: Table S3.

Assessment of feature quality.

Additional file 4: Figure S1.

Receiver operating characteristic (ROC) curves of 36 feature values. Features with areas under the curve (AUC) above 0.6 were selected for use in the SVM model to predict functional interactions between genes.

Additional file 5: Table S4.

Evaluation of the predicted interactions in different datasets.

Additional file 6: Table S5.

Annotations produced by the GO enrichment analysis tool for the transcriptionally changed genes.

Additional file 7: Table S6.

Annotations produced by DAVID for the transcriptionally changed genes.

Additional file 8: Table S7.

Annotations produced by FIC/GSLA for the transcriptionally changed genes.

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 http://creativecommons.org/licenses/by/4.0/. 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 in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Chen, PC., Ruan, L., Jin, J. et al. Predicted functional interactome of Caenorhabditis elegans and a web tool for the functional interpretation of differentially expressed genes. Biol Direct 15, 20 (2020). https://doi.org/10.1186/s13062-020-00271-6

Download citation

Keywords

  • Caenorhabditis elegans
  • Functional interaction
  • Database
  • Gene set linkage analysis
  • Transcriptomic analysis tool