 Software
 Open Access
 Published:
SPINNAKER: an Rbased tool to highlight key RNA interactions in complex biological networks
BMC Bioinformatics volume 23, Article number: 166 (2022)
Abstract
Background
Recently, we developed a mathematical model for identifying putative competing endogenous RNA (ceRNA) interactions. This methodology has aroused a broad acknowledgment within the scientific community thanks to the encouraging results achieved when applied to breast invasive carcinoma, leading to the identification of PVT1, a long noncoding RNA functioning as ceRNA for the miR200 family. The main shortcoming of the model is that it is no freely available and implemented in MATLAB®, a proprietary programming platform requiring a paid license for installing, operating, manipulating, and running the software.
Results
Breaking through these model limitations demands to distribute it in an opensource, freely accessible environment, such as R, designed for an ordinary audience of users that are not able to afford a proprietary solution. Here, we present SPINNAKER (SPongeINteractionNetworkmAKER), the opensource version of our widely established mathematical model for predicting ceRNAs crosstalk, that is released as an exhaustive collection of R functions. SPINNAKER has been even designed for providing many additional features that facilitate its usability, make it more efficient in terms of further implementation and extension, and less intense in terms of computational execution time.
Conclusions
SPINNAKER source code is freely available at https://github.com/sportingCode/SPINNAKER.git together with a thoroughgoing PPTbased guideline. In order to help users get the key points more conveniently, also a practical Rstyled plaintext guideline is provided. Finally, a short movie is available to help the user to set the own directory, properly.
Background
microRNAs (miRNAs) are singlestranded short RNAs (about 20–22 nucleotide) that posttranscriptionally regulate gene expression by degradation or translation inhibition of their target messenger RNAs (mRNAs). Most of the biological processes have been shown to entail the regulation orchestared by miRNAs, such as cell proliferation, differentiation, metabolism, development, and apoptosis [1]. A novel mechanism of miRNA regulation regarding the ability of RNAs to compete for the miRNA binding has recently been discovered [2, 3]. Key triggers of this new layer of posttranscriptional regulation are the socalled competing endogenous RNAs (ceRNAs)—or miRNA ‘sponges’, involving both coding and noncoding RNAs, such as pseudogenes [2], circular RNAs [4, 5], and long non coding RNAs (lncRNAs) [6, 7]. ceRNAs exert their decoy activity by recruiting miRNA molecules through basepairing with miRNArecognition elements (MREs) that they share with a target, consequently determining the target release from the miRNA control. The ceRNA crosstalks (i.e., their interactions mediated by miRNAs) have been identified as drivers of most of the pathological conditions, including human cancers [2, 8,9,10].
The research field related to ceRNA mechanism has rapidly grown during the last decade, as observed by looking at the increasing number of published studies over the last years (Additional file 1: Figure S1). Meanwhile, the recent years have seen the development of several computational methodology to build ceRNA networks. Yet, there remain great opportunities as well as challenges to propose computational models helping the generation of hypotheses able to drive wetlab experiments towards the elucidation of the roles of ceRNAs in a particular disease [11].
Since considering all types of miRNA sponges can lead to a high computation complexity, in our recent study we focused on lncRNAs given their acknowledged importance in diverse biological and physiopathological contexts [12] and we developed a new methodology suitable to exploring their potential role as ceRNA regulators [13]. According to a recent review, comparing the most widespread computational models for ceRNAceRNA interactions’ identification [11], our method proved to be the best in terms of the percentage of predicted RNAs acting as ceRNAs related to breast invasive carcinoma. However, the main limitation towards its broader usability is that it was no freely distributed so far, and originally written in MATLAB, a proprietary programming language requiring a paid license to install, exploit, operate, and run the software. The interesting and acknowledged results obtained by applying our model to breast cancer and the necessity to spread it to a broader scientific audience created the call to handle this issue and to design an opensource version of the model for a universal community of nonexpert users.
Here, we present SPINNAKER, the implementation of the originally developed ceRNA model [13] in a exhaustive collection of R functions. Yet, SPINNAKER comes as a simplified and an improved version with respect to the original MATLABbased implementation (Additional file 2), ensuring more efficiency by speeding up the entire pipeline of many orders of magnitude (i.e., from many hours to minutes). In particular, an additional feature of SPINNAKER is the possibility of choosing among different pools of RNAs acting as ceRNAs, as long as the total number of triplets to be tested is within the order of magnitude of \({10}^{6}\), otherwise it collides with a huge computation complexity.
Implementation
SPINNAKER (SPongeINteractionNetworkmAKER) is an Rbased implementation of a methodology for identifying putative ceRNA interactions that we recently published along with its application in breast invasive carcinoma [13]. SPINNAKER takes as input normalized expression levels of RNAs and miRNAs (e.g., FPKM) and predicts the ceRNA interaction network by implementing two modules: 1) data collection and processing, 2) ceRNA network building (Fig. 1). Each module consists of several steps detailed in the following.
Module 1: data collection and processing
The goal of this module is to collect and process data for running SPINNAKER.

(i)
Selection of RNA molecular biotype
To define the molecular entities competing for the miRNA binding, SPINNAKER automatically queries NCBI's Gene database, including information about chromosomal localization, nomenclature, gene products, and their attributes (e.g., molecular biotype). Then, SPINNAKER separates the two selected classes of candidate ceRNAs to be tested (e.g., protein coding versus long noncoding RNAs).

(ii)
Collection of miRNAtarget interactions
SPINNAKER collects miRNA sequences from miRBase (currently release 22.1, October 2018, http://www.miRBase.org). Then, it retrieves the predictions of miRNAmRNA target interactions from TargetScan [14] and the predictions of miRNAlncRNA target interactions from miRWalk [15]. TargetScan appears as the most uptodate database for sequencebased predictions of miRNAtarget interactions, it predicts miRNAs targets by considering the exact matching between the seed region of a miRNA (i.e., positions 2–7 from the miRNA 5′end) and the 3′ UTR of its targets [16]. miRWalk provides the predicted and validated miRNAbinding sites of known genes of human and other species, most importantly including lncRNAsmiRNA interactions, and it is entirely updated more than once a year. miRWalk core is based on TarPmiR [17], a prediction tool that exploits an approach based on randomforest algorithm to look for putative miRNA binding sites within the whole transcript sequence including the 3′UTR, 5′UTR, and CDS [15].

(iii)
Removing missing values and log2 transformation
SPINNAKER applies a logarithmic (log2) transformation to the RNAs and miRNAs expression levels and conducts a processing analysis to remove those genes having too many missing values among the samples (i.e., by default, SPINNAKER filters out entries showing missing values for more than ten percentage of samples).
Module 2: ceRNA network building
The goal of this module is to build the ceRNA interaction network.

(i)
Computation of Pearson correlation
SPINNAKER computes the Pearson correlation coefficients between the expression profiles of the RNA pairs \(({\rho }_{XY})\).

(ii)
Selection of highly correlated pairs
SPINNAKER selects the RNA pairs with \({\rho }_{XY}\) greater than a defined threshold (by default equal to 99th percentile) on the overall correlation distribution. This threshold allows to reduce both the computational effort in evaluating the RNA interactions and the number of false positives.

(iii)
Computation of Pearson correlation with miRNAs
SPINNAKER computes the Pearson correlation coefficient between the expression profiles of RNA X and miRNA Z (\({\rho }_{XZ}\)) and the expression profiles of RNA Y and miRNA Z (\({\rho }_{YZ}\)).

(iv)
Computation of sensitivity correlation
To determine if the Pearson correlation between the RNA pairs is direct or mediated by the miRNA, SPINNAKER implements the following metric, called sensitivity correlation S (Fig. 2):
$$S={\rho }_{XY}{\rho }_{XYZ}$$with \({\uprho }_{XY}\) referring to the Pearson correlation coefficient between RNA X and RNA Y, and \({\uprho }_{XYZ}\) referring to the partial correlation between RNA X and RNA Y controlling for the miRNA Z defined as:
$${\rho }_{XYZ}= \frac{{\rho }_{XY}{\rho }_{XZ}{\rho }_{YZ}}{\sqrt{1{\rho }_{XZ}^{2}}\sqrt{1{\rho }_{YZ}^{2}}}$$with \({\rho }_{XZ}\) (\({\rho }_{YZ}\)) referring to the Pearson correlation between RNA X (RNA Y) and miRNA Z [18]. The partial correlation \({\rho }_{XYZ}\) measures how much the correlation between two variables X (RNA 1 expression profile) and Y (RNA 2 expression profile) remains after removing a third variable Z (miRNA expression profile). Thus, a low value of sensitivity correlation (i.e., partial correlation approaching the Pearson correlation) refers to a direct interaction between the two RNAs competing for the same miRNA, i.e., whose interaction is not mediated by the miRNA (Fig. 2a); while a high value of sensitivity correlation (i.e., partial correlation approaching to zero) refers to an indirect interaction between the two RNAs competing for the same miRNA, i.e., whose interaction is mediated by the miRNA (Fig. 2b).
An example of output of this step, which describes the sensitivity correlation obtained for normal breast tissues [13], is presented in Fig. 2c. The computed sensitivity matrix is rendered as an heatmap, where rows represent the highly correlated RNA pairs, columns refer to all the analysed miRNA, and sensitivity values are colorcoded increasing from red to blue.
The sensitivity correlation calculated in normal breast samples [13] unveiled a general trend of the interactions between RNA pairs marked by highly positive correlations that appear miRNAindependent (S∼0, red background in Fig. 2c), except for a small pool of miRNAs for which the sensitivity appears to be pretty different from zero (bright vertical stripes in Fig. 2c). This observation highlights the presence of particular miRNAs (such as miR200 family) able of generating a crosstalk throughout the whole transcriptome. These miRNAs, corresponding to the bright vertical stripes in the sensitivity heatmap, represent the links of ceRNA network generated by SPINNAKER, whereas the highly correlated RNA pairs with high sensitivity correlation mediated by these miRNAs represent the nodes of the ceRNA network.
It is worth noting that the outcome of SPINNAKER might not be rendered as a single heatmap, since the size of the picture could exceed the memory limits of computer. To overcome this limitation, SPINNAKER segregates the sensitivity matrix into a variable number of heatmaps (depending on the total number of RNA pairs), each one composed of a maximum number of 5000 rows (RNA pairs).

(v)
Computation of ceRNA interaction network.
SPINNAKER selects the XYZ triplets with S greater than a defined threshold (by default equal to the 99^{th} percentile) on the overall distribution of the Svalues. This threshold allows to reduce both the computational effort in evaluating the ceRNA interactions and the number of false positives. The X and Y variables correspond to the topcorrelated RNA pairs. Then, SPINNAKER builds the ceRNA interaction network.
Nodes in the ceRNA network represent ceRNAs marked by a high correlation between their expression profiles; whereas links represent miRNAs that are mediating their interactions. A link between two nodes (ceRNAs) occurs if they fulfilled the following conditions (Fig. 3): (1) showing a high Pearson correlation value; (2) showing a high sensitivity correlation value.

(vi)
[Optional] Search seedmatch for all triplets
SPINNAKER searches for the seedmatch of all the highly correlated pairs with the miRNA mediating their interactions, in order to narrow the above selected triplets (step v.) to those including only ceRNAs that are targets of the shared miRNA.

(vii)
[Optional] Computation of statistical analysis
For each miRNA, SPINNAKER performs a seedmatch enrichment analysis by computing the following statistics (i.e., pvalue resulting from the hypergeometric test) (Fig. 4):
$$p=1\sum_{i=0}^{X1}\frac{\left(\genfrac{}{}{0pt}{}{K}{i}\right)\left(\genfrac{}{}{0pt}{}{UK}{Si}\right)}{\left(\genfrac{}{}{0pt}{}{U}{S}\right)}=\sum_{i=X}^{S}\frac{\left(\genfrac{}{}{0pt}{}{K}{i}\right)\left(\genfrac{}{}{0pt}{}{UK}{Si}\right)}{\left(\genfrac{}{}{0pt}{}{U}{S}\right)}$$where U is the universe dimension, that is the number of the topcorrelated RNA pairs; K is the property, that is the number of RNA pairs sharing the binding site for the miRNA under test; S is the selection, that is the number of RNA pairs with high sensitivity for the miRNA under test; X is the number of RNA pairs with sensitivity correlation exceeding the defined threshold on the Svalues distribution and sharing the binding site for the miRNA under test.
The output of this step is the ceRNA interaction network with the additional information of pvalues, and where a link between two nodes (ceRNAs) occurs if they fulfilled the following conditions (Fig. 3): (1) showing a high Pearson correlation value; (2) showing a high sensitivity correlation value; and (3) sharing the binding sites for the miRNAs mediating their interaction.
An example of output, which describes the ceRNA network built for normal breast tissues [13], is presented in Fig. 5. The ceRNA network released by SPINNAKER can be easily uploaded on Cytoscape software [19], an open source tool for complex networks visualization and integration with other types of attribute data. By using Cytoscape, the ceRNA network can be visualized by adding a different color to each miRNA and analyzed by using several builtin apps, like the “network analyzer” tool. In this way, the user can obtain the total number of connected components of the ceRNA network and other network properties associated to each node, such as betweenness, closeness, clustering coefficient, degree, and thus consequently identified the network hubs (i.e., nodes with more than 5 links [20]).
In particular, the ceRNA network of normal breast (32,375 links and 1738 nodes) is markedly distinguished into two wellconnected components: a smaller one (954 links and 378 nodes) mainly dominated by miR452 whose mRNAs functioning as ceRNAs were functionally enriched in cellular metabolic processes, and a larger one (31,417 links and 1354 nodes) mainly controlled by the miR200 family and functionally enriched in cell–cell adhesion functions [13]. Of note, the first hub of ceRNA normal breast network was the lncRNA PVT1, revealing a net binding preference with the miR200 family and competing with mRNAs mostly associated to cancer development and progression (e.g., CDH1, GATA3, RUNX1, RUNX3, TP53, TP63, TP73).
Results
Comparisons with other methods
The widespread computational approaches developed for modelling ceRNA regulatory mechanism in human cancers can be grouped into two main classes: (i) statisticsbased methods, which exploited statistical analyses (e.g., multivariate analysis) to infer putative evidences of ceRNA crosstalk and construct ceRNA interaction networks by considering miRNA expression levels (Additional file 1: Table S1); (ii) mathematical methods, which make use of stochastic or deterministic models to predict and analyze the behavior of ceRNA crosstalk (Additional file 1: Table S2). The model implemented by SPINNAKER, based on sensitivity metric, is included in the statisticsbased methods.
In order to evaluate the effectiveness of SPINNAKER predictions, we exploited the comparison study conducted in a recent review article [11], where the authors compared the methodology implemented by SPINNAKER [13] with respect to other computational methods, when applied to breast cancer dataset. In particular, to conduct the comparison, they chose to focus on statisticsbased methods, since approaches based on mathematical modelling only quantitatively describe a minimum number of ceRNA interaction network, and hence they can be only used to investigate a small number of ceRNA interactions. Among the statisticsbased methods, they selected as representative those which are based on Conditional Mutual Information (CMI), Positive Correlation (PC), and Partial Pearson Correlation (PPC) [11]. These methods, as SPINNAKER, encompass two principal parts: computing the candidate ceRNA interactions; assessing the ceRNA interactions. The effectiveness of these methods was evaluated in terms of their ability to recover breast cancer related and/or confirmed sponge interactions by exploiting publicly available databases of computationally predicted and/or experimentally validated ceRNA interactions (Additional file 1: Table S3).
Specifically, the authors considered breast cancer related ceRNA interactions, the ones where the two ceRNAs are genes associated to breast cancer, according to the broad range of experimentally validated databases. CMI predicted a total of 509,504 sponge interactions, among which 20,082 associated to breast cancer and 5 of them experimentally validated. PC method predicted a total of 1,274,708 sponge interactions, among which 46,716 associated to breast cancer and 7 of them experimentally validated. PPC method predicted a total of 702,029 sponge interactions, with 26,803 interactions associated to breast cancer, but no one was experimentally confirmed. The model implemented by SPINNAKER predicted a total of 99,662 sponge interactions, with 4223 associated to breast cancer and 2 of them experimentally confirmed in breast cancer (CNOT6LPTEN and ZEB2PTEN). Thus, from this comparative analysis [11], SPINNAKER revealed the highest percentage of identified ceRNA interactions related to breast invasive carcinoma, then resulting as the best method (Fig. 1). It is worth noting that the total number of interactions predicted by SPINNAKER is much lower than other methods shown in the table of Fig. 1. This due to the fact that the rationale behind SPINNAKER methodology is to decrease the number of false positive values as much as possible, by using very high thresholds both in the data collection and processing and ceRNA network building module. This strategy resulted in the higher percentage of diseasespecific interactions, as witnessed by the comparison between SPINNAKER and other statistical methods when applied to breast carcinoma, and by the comparison with other more recently developed methods for studying various diseases [21,22,23,24,25,26,27]. In Table 1, the comparison between SPINNAKER and two of these methods [21, 22] is shown for thyroid carcinoma (thca). Once again, SPINNAKER resulted the best one in predicting the higher percentage of thyroid carcinoma related genes acting as ceRNAs (Additional file 3).
Experimentally validated interactions
The most significant prediction of the computational model implemented by SPINNAKER was the discovery of lncRNA PVT1 acting as ceRNA in breast invasive carcinoma (brca) dataset, where it antagonized the miR200 family to regulate the expression of several messenger RNAs [13]. This finding was confirmed by applying SPINNAKER on the uptodate TCGA brca dataset with an increased number of patients as in the original paper (Additional file 4) [13]. As further proof of the reliability of our methodology, this result was experimentally validated in a recent research study conducted in nonsmall cell lung cancer cells, where the authors showed that PVT1 facilitates the tumor invasion functioning as ceRNA to regulate the MMP9 expression via the competitively binding of miR200 family [7].
Performance evaluation
In order to quantify the speedup reached by SPINNAKER implementation, we compared the elapsed times of running it on MALTAB and R platforms, when applied to brca dataset [13]. We observed that SPINNAKER is from 5x (Module 1) up to 100x (Module 2) faster than its ancestor running on MATLAB (Table 2). Furthermore, we checked that all the results obtained by running the original code (MATLABbased) were confirmed by SPINNAKER (Rbased), when applied to the original brca dataset (Additional file 5). We can conclude that the Rbased solution can lead to gain a significant amount of time, thus greatly affecting the whole analysis process.
Conclusions
In this study, we presented SPINNAKER (SPongeINteractionNetworkmAKER), the R implementation and opensource version of a widely established mathematical model that we published for identifying putative competing endogenous RNA (ceRNA) interactions. According to a recent review [11], the methodology implemented by SPINNAKER resulted as the best one in terms of the percentage of discovered ceRNA interactions associated with breast invasive carcinoma. However, the main limitation towards a broader usability of this methodology is that it was developed in MATLAB®, a proprietary programming environment requiring a paid license for installing, operating, and running the software. To solve this issue, SPINNAKER came as an Rbased, opensource, simplified, and improved version with respect to the original MATLABbased implementation, ensuring a greater efficiency by speeding up the whole process of several orders of magnitude. By comparing SPINNAKER with other statisticalbased methods, once again it resulted as the best one in terms of higher percentage of diseaseassociated genes acting as ceRNAs when applied to thyroid carcinoma dataset.
Availability and requirements
Project name: SPINNAKER
Operating system(s): Windows 10 Pro, Ubuntu 20.04.3 LTS, macOS High Sierra 10.13.6
Programming language: R
Project page: https://github.com/sportingCode/SPINNAKER.git
Other requirements: R version 3.5.1, R 4.1.2 or higher
License: GNU AFFERO GENERAL PUBLIC LICENSE
Any restrictions to use by nonacademics: license needed
Availability of data and materials
SPINNAKER source code is freely available at https://github.com/sportingCode/SPINNAKER.git along with a comprehensive PPTbased guideline. In order to help users get the key points more conveniently, also a practical Rstyled plaintext guideline is provided. Finally, a short movie is available to help the user to set the own directory, properly.
Abbreviations
 ceRNA:

Competing endogenous RNA
 SPINNAKER:

SPongeINteractionNetworkmAKER
 lncRNAs:

Long noncoding RNAs
 brca:

Breast invasive carcinoma
 thca:

Thyroid carcinoma
 TCGA:

The Cancer Genome Atlas
 miRNAs:

MicroRNAs
References
Gebert LFR, MacRae IJ. Regulation of microRNA function in animals. Nat Rev Mol Cell Biol. 2019;20:21–37.
Poliseno L, Salmena L, Zhang J, et al. A codingindependent function of gene and pseudogene mRNAs regulates tumour biology. Nature. 2010;465:1033–8.
Salmena L, Poliseno L, Tay Y, et al. A ceRNA hypothesis: the Rosetta Stone of a hidden RNA language? Cell. 2011;146:353–8.
Memczak S, Jens M, Elefsinioti A, et al. Circular RNAs are a large class of animal RNAs with regulatory potency. Nature. 2013;495:333–8.
Kristensen LS, Jakobsen T, Hager H, et al. The emerging roles of circRNAs in cancer and oncology. Nat Rev Clin Oncol. 2022;19:188–206.
Cesana M, Cacchiarelli D, Legnini I, et al. A long noncoding RNA controls muscle differentiation by functioning as a competing endogenous RNA. Cell. 2011;147:358–69.
Chen W, Zhu H, Yin L, et al. lncRNAPVT1 facilitates invasion through upregulation of MMP9 in nonsmall cell lung cancer cell. DNA Cell Biol. 2017;36:787–93.
Sumazin P, Yang X, Chiu HS, et al. An extensive microRNAmediated network of RNA–RNA interactions regulates established oncogenic pathways in glioblastoma. Cell. 2011;147:370–81.
Yang L, Peng X, Jin H, et al. Long noncoding RNA PVT1 promotes autophagy as ceRNA to target ATG3 by sponging microRNA365 in hepatocellular carcinoma. Gene. 2019;697:94–102.
He Y, Huang H, Jin L, et al. CircZNF609 enhances hepatocellular carcinoma cell proliferation, metastasis, and stemness by activating the Hedgehog pathway through the regulation of miR15a5p/15b5p and GLI2 expressions. Cell Death Dis. 2020;11:1–12.
Le TD, Zhang J, Liu L, et al. Computational methods for identifying miRNA sponge interactions. Brief Bioinform. 2017;18:577–90.
Statello L, Guo CJ, Chen LL, et al. Gene regulation by long noncoding RNAs and its biological functions. Nat Rev Mol Cell Biol. 2021;22:96–118.
Paci P, Colombo T, Farina L. Computational analysis identifies a sponge interaction network between long noncoding RNAs and messenger RNAs in human breast cancer. BMC Syst Biol. 2014;8:83.
Agarwal V, Bell GW, Nam JW, et al. Predicting effective microRNA target sites in mammalian mRNAs. Elife. 2015;4: e05005.
Sticht C, Torre CDL, Parveen A, et al. miRWalk: an online resource for prediction of microRNA binding sites. PLoS ONE. 2018;13: e0206239.
RiffoCampos ÁL, Riquelme I, BrebiMieville P. Tools for sequencebased miRNA target prediction: what to choose? Int J Mol Sci. 2016;17:1987.
Ding J, Li X, Hu H. TarPmiR: a new approach for microRNA target site prediction. Bioinforma Oxf Engl. 2016;32:2768–75.
Conte F, Fiscon G, Sibilio P, et al. An overview of the computational models dealing with the regulatory ceRNACeRNAsMechanism and ceRNACeRNAsDeregulation in cancer. In: Poliseno L, editor., et al., Pseudogenes: functions and protocols. New York: Springer; 2014. p. 149–64.
Shannon P, Markiel A, Ozier O, et al. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13:2498–504.
Han JDJ, Bertin N, Hao T, et al. Evidence for dynamically organized modularity in the yeast protein–protein interaction network. Nature. 2004;430:88–93.
Zhao Y, Wang H, Wu C, et al. Construction and investigation of lncRNAassociated ceRNA regulatory network in papillary thyroid cancer. Oncol Rep. 2018;39:1197–206.
Jiang Y, Wang J, Chen J, et al. Construction and analysis of an aberrant lncRNAmiRNAmRNA network associated with papillary thyroid cancer. Medicine (Baltimore). 2020;99: e22705.
Du J, Li G, Li Y, et al. Identification of prognostic biomarkers related to the tumor microenvironment in thyroid carcinoma. Sci Rep. 2021;11:16239.
Ning P, Wu Z, Hu A, et al. Integrated genomic analyses of lung squamous cell carcinoma for identification of a possible competitive endogenous RNA network by means of TCGA datasets. PeerJ. 2018;6: e4254.
Qi L, Zhang T, Yao Y, et al. Identification of lncRNAs associated with lung squamous cell carcinoma prognosis in the competitive endogenous RNA network. PeerJ. 2019;7: e7727.
Li R, Yang YE, Jin J, et al. Identification of lncRNA biomarkers in lung squamous cell carcinoma using comprehensive analysis of lncRNA mediated ceRNA network. Artif Cells Nanomedicine Biotechnol. 2019;47:3246–58.
Hunt M, Banerjee S, Surana P, et al. Small RNA discovery in the interaction between barley and the powdery mildew pathogen. BMC Genomics. 2019;20:610.
Acknowledgements
Not applicable.
Funding
This work has been partially funded by the BiBiNet project (grant number: H35F21000430002) within the PORLazio FESR 2014–2020, by PRIN 2017—Settore ERC LS2—Codice Progetto 20178L3P38, and by Sapienza University of Rome grant entitled “Sviluppo di un algoritmo basato sulla teoria delle reti per il riposizionamento dei farmaci e la sua applicazione al COVID19” (grant number: RM12117A34663A2C). The funding body played no roles in the design of the study and collection, analysis, and interpretation of data and in writing the manuscript.
Author information
Authors and Affiliations
Contributions
PP conceived and designed the research, and developed the methodology. PP and GF developed the software in R language and wrote the manuscript. All the authors read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Not applicable.
Consent for publication
Not applicable.
Competing interest
The authors declare that they have neither financial or nonfinancial 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:
It includes the Figure S1 and Tables S1–S3.
Additional file 2:
Performance comparison between the original code (MATLAB) versus SPINNAKER (R).
Additional file 3:
ceRNA interaction network obtained by running SPINNAKER on thyroid carcinoma dataset.
Additional file 4:
ceRNA interaction network obtained by running SPINNAKER on breast cancer dataset.
Additional file 5:
Results comparison between the original code (MATLAB) versus SPINNAKER (R).
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.
About this article
Cite this article
Paci, P., Fiscon, G. SPINNAKER: an Rbased tool to highlight key RNA interactions in complex biological networks. BMC Bioinformatics 23, 166 (2022). https://doi.org/10.1186/s1285902204695x
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1285902204695x
Keywords
 Algorithms
 ceRNA network
 Network theory
 miRNA sponge