Access the full text.
Sign up today, get DeepDyve free for 14 days.
Abstract Motivation Moonlighting proteins are a class of proteins having multiple distinct functions, which play essential roles in a variety of cellular and enzymatic functioning systems. Although there have long been calls for computational algorithms for the identification of moonlighting proteins, research on approaches to identify moonlighting long non-coding RNAs (lncRNAs) has never been undertaken. Here, we introduce a novel methodology, MoonFinder, for the identification of moonlighting lncRNAs. MoonFinder is a statistical algorithm identifying moonlighting lncRNAs without a priori knowledge through the integration of protein interactome, RNA–protein interactions and functional annotation of proteins. Results We identify 155 moonlighting lncRNA candidates and uncover that they are a distinct class of lncRNAs characterized by specific sequence and cellular localization features. The non-coding genes that transcript moonlighting lncRNAs tend to have shorter but more exons and the moonlighting lncRNAs have a variable localization pattern with a high chance of residing in the cytoplasmic compartment in comparison to the other lncRNAs. Moreover, moonlighting lncRNAs and moonlighting proteins are rather mutually exclusive in terms of both their direct interactions and interacting partners. Our results also shed light on how the moonlighting candidates and their interacting proteins implicated in the formation and development of cancers and other diseases. Availability and implementation The code implementing MoonFinder is supplied as an R package in the supplementary material. Supplementary information Supplementary data are available at Bioinformatics online. 1 Introduction Protein moonlighting is a common phenomenon in nature involving a protein with a single polypeptide chain that can perform more than one independent cellular function (Boukouris et al., 2016; Monaghan and Whitmarsh, 2015). Enzymes, receptors, ion channels or chaperones are the typical form of moonlighting proteins (MPs). Enzyme is the most common form of moonlighting proteins whose primary function is enzymatic catalysis, but they are also in possession of additional roles such as signal transduction, transcriptional regulation, apoptosis, motility and structural proteins (Jeffery, 2015). For example, crystallins, a class of well-studied moonlighting proteins, function as enzymes when expressed at low levels in many tissues, but are densely packed to form lenses when expressed at high levels in eye tissue (Piatigorsky et al., 1988; Piatigorsky and Wistow, 1989). The genes encoding crystallins need to sustain functions of both catalytic and transparency maintenance. Another example is glycolysis, an ancient universal metabolic pathway, in which a high percentage of proteins are moonlighting proteins (Boukouris et al., 2016; Sriram et al., 2005). Moreover, some proteins work on their moonlighting by being assembled to supramolecules, such as the ribosome, which usually composed of more than a hundred of proteins and RNAs. However, the studies of moonlighting merely concentrated on proteins and the genes coding them, yet the moonlighting of non-coding RNAs has not been investigated, despite the fact that ncRNAs have gained widespread attention due to their functional importance over the last decade (Chen, 2016a,b; Liao et al., 2011; Quinn and Chang, 2016; Zhou et al., 2017). Currently, the information of MPs, such as protein functions, cell localization and primary structures, is scattered across a number of publications, since the MPs perform a variety of functions in different tissues and cell types. Some researchers have summarized the literature about MPs from different aspects of the functional diversity, such as regulation circuits or signaling pathways. The Jeffery lab constructed a manually curated database MoonProt, which consists of over 200 MPs that have been experimentally verified (Mani et al., 2015). The structures and function information about the MPs can aid researchers to understand how proteins function in a moonlighting manner and help in designing proteins with novel functions. Min et al. summarized the MPs from the perspective of a coupled intracellular signaling pathway (Min et al., 2016). Numerous proteins are localized in more than one compartment in cells and the aberrant translocation of proteins may cause cancer or other disorders. Hence, it is necessary to study the localization dynamic and trans localization activity of MPs. Monaghan et al. reviewed several MPs with dual mitochondrial and nuclear functions (Monaghan and Whitmarsh, 2015). It is pointed out that the nuclear moonlighting of mitochondrial proteins is part of a mitochondria-to-nucleus signaling pathway in cells. They also discussed various mechanisms commanding the dual localization of proteins and indicated that the nuclear moonlighters perform as a regulating loop to maintain homeostasis in mitochondria. Boukouris et al. summarized the moonlighting functions of metabolic enzymes in the nucleus (Boukouris et al., 2016). They proposed a new mechanistic connection between metabolic flux and differential expression of genes, which is implemented via nuclear translocation or moonlighting of nuclear metabolic enzymes. This mechanistic link aids cells in adapting a changing environment in normal and disease states, such as cancer, and thus has the potential to be explored for the novel therapeutic target. In parallel to the serendipitous findings of MPs through experiments, some computational approaches have been developed to predict MPs in recent years (Pritykin et al., 2015). Specifically, three algorithms were proposed for moonlighting protein identification, MoonGO (Chapple et al., 2015), MPFit (Khan and Kihara, 2016) and DextMP (Khan et al., 2017), executing statistic, machine learning and text mining techniques, respectively. These studies investigated different aspects of MPs such as conserved sequence domains, structural disorders, protein interaction patterns and network topology. MoonGO first identifies overlapping protein clusters in the human interactome (Chapple et al., 2015). Then, the clusters are annotated to the Gene Ontology (GO) terms of biological process. GO terms annotating more than half of a cluster’s proteins are assigned to the cluster. Each individual protein then inherits the annotations of its clusters in addition to its own. Finally, the proteins shared by dissimilar functional clusters are identified as MPs. MPFit uses a variety of protein features to address the diverse nature of MPs, including functional annotation, protein interactions, gene expression, phylogenetic profiles, genetic interactions, network-based graph properties and protein structural properties (Khan and Kihara, 2016). In general, MPs are assigned in more clusters because they interact with proteins of diverse functions, so the number of clusters that a protein involved is used as an omics feature. For proteins that do not have an available record of certain features, an imputation step using random forest is used to predict the missing features. Eventually, these features are combined with machine learning classifiers to make moonlighting protein prediction. DextMP is a text mining algorithm consisting of four logical steps (Khan et al., 2017). First, it extracts textual information of proteins from literature and functional description in the UniProt database. Next, it constructs a k-dimensional feature vector from each text using three language models, i.e. paragraph vector, Term Frequency-Inverse Document Frequency (TFIDF) in the bag-of-words category and Latent Dirichlet Allocation (LDA) in the topic modeling category. Third, using four machine learning classifiers, a text is classified to either MP or non-MP based on the text features. Finally, the text predictions for each protein are separately summarized to predict which ones are MPs. Long non-coding RNAs (lncRNAs) is a subclass of non-coding RNAs with little coding potential whose transcript consists of no less than 200 nucleotides. lncRNAs are implicated in a variety of biological processes through diverse functional mechanisms such as chromatin remodeling, chromatin interactions and functioning as competing endogenous RNAs (Ferre et al., 2016). Specific expression patterns of lncRNA in cells correspond to certain cell development and disorder. Nuclear and cytoplasmic lncRNAs can regulate gene expression and function in multiple ways, e.g. i) affecting mRNA translation directly, ii) interfering with protein post-translational modifications to disturb signal transduction, iii) functioning as decoys for miRNAs and proteins, iv) acting as miRNA sponges, v) interacting proteins to enhancer regions and vi) encoding small proteins and functional micro peptides, etc. (Cabili et al., 2015; Ferre et al., 2016; Quinn and Chang, 2016; Zhou et al., 2017; Zhu et al., 2016). Many lncRNAs diversely reside in the nucleus and play an essential role as modulators for nuclear functions. Some others are translocated to the cytoplasm to enforce their regulatory roles. Moreover, lncRNAs have a variety of subcellular localization patterns, which are not limited to specific nuclear and cytoplasm localization but also nonspecific localization in both the nucleus and cytoplasm (Barabasi and Oltvai, 2004; Buxbaum et al., 2015). For the lncRNAs localized in multiple compartments, the intercommunication can modulate the interaction pattern or expression abundance, e.g. regulating the lncRNA abundance in one compartment may influence the function of the other cell compartment. Also, inappropriate moonlighting may trigger genetic diseases (Abumrad and Lange, 2006; Espinosa-Cantu et al., 2015; Min et al., 2016). Hence, it is necessary to study the localization dynamic and expression activity of moonlighting lncRNAs (mlncRNAs) and to investigate the mechanism of how the mlncRNAs modulate and switch the functions in the metabolic processes, which is of vital importance for cancer therapeutics and will provide tremendous opportunities for anti-cancer purposes (Du et al., 2013; Liu et al., 2014; Wang et al., 2015; Zhu et al., 2016). We have demonstrated that using clustering algorithms is able to group proteins into functional modules allowing the identification of MPs (Chapple et al., 2015; Khan et al., 2014; Pritykin et al., 2015). A module corresponds to a functional unit, which is composed of several closely interacted proteins involved in specific tasks in the cell. Therefore, it is promising to use the functional module approach to identify mlncRNAs that exhibit multiple but distinct functionalities. Our study is focused on the moonlighting of human lncRNAs, since lncRNAs are pervasively transcribed in the mammalian genome and several of them play the roles as oncogenic or tumor-suppressor genes in multiple cancers (Ning et al., 2016; Quinn and Chang, 2016; Wahlestedt, 2013; Zhu et al., 2016). We first propose a novel algorithm MoonFinder to identify mlncRNAs that have multiple but unrelated functions. Then, we characterize the sequence features and the localization tendency of these mlncRNA candidates. After that, we construct an mlncRNA–module network and topologically analyzed the mlncRNAs with regard to the association with cancer, diseases, drug targets and moonlighting proteins. We also predict two cancer lncRNAs exclusively interacted with functional modules from the network. 2 Material and methods 2.1 The workflow of MoonFinder Moonlighting lncRNAs (mlncRNAs) may perform their distinct functions through interacting with different functional modules and module is a cluster of proteins interconnected densely and associated with specific tasks in cells. Hence, it is promising to identify mlncRNAs from the aspect of whether they are targeting functionally unrelated protein modules. MoonFinder integrates several statistical models based on the RNA and protein interactome as well as the protein functional annotations for the identification of moonlighting macromolecules. The workflow of MoonFinder contains the following six steps, Protein interaction network refinement. Refine the protein interaction network by filtering out protein pairs sharing no identical cell compartments and only the co-localized interactions are considered for further analysis. Protein module identification. Detect protein modules from the co-localized protein interaction network (using ClusterONE by default). Each of the detected modules is highly interconnected and expected to be implicated in a specific cellular process. Functional annotation of modules. Establish the annotation of the modules with GO terms by performing the functional category enrichment using hypergeometric test (HGT). The pairwise association between the modules and the GO terms are constructed. Establish RNA–module interactions. Assess whether the target proteins of an RNA are significantly overrepresented in a module. If yes, we define the RNA functionally regulates the module. Construct similarity matrix of modules. All pairwise functional similarities are pre-calculated by a given semantic similarity measurement using a BMA combination approach and eventually a similarity matrix is produced. Moonlighting determination. Calculate the principal components (PCs) and their contribution weights of the similarity matrix of the modules targeted by an RNA using principal component analysis (PCA). An RNA is determined as moonlighting if more than one PC play an indispensable role (weight > 0.1) and none of these PCs play a dominant role (weight > 0.7). The weight can be customized. The higher the value, the more moonlighting candidates can be identified. For instance, we say an RNA is moonlighting if the weights of the top three PCs of the module similarity matrix are 0.4, 0.3 and 0.2, respectively, because three PCs (more than one) are substantial. On the other hand, an RNA is multifunctional but not moonlighting if the weight of the first PC is 0.85, the weight of the second PC is 0.1, and the rest share the negligible weight of 0.05. The workflow of MoonFinder is explained in more detail in Figure 1. Fig. 1. View largeDownload slide Schematic diagram of MoonFinder. The orange, green and blue boxes represent the input, generated, and output data, respectively. The clustering algorithm and statistical models are shown in the dotted boxes. The inputs are the four datasets, i.e. interactions of proteins, localization annotation of proteins, function annotation of proteins, and interactions between ncRNAs and proteins, respectively. The output is a list of predicted moonlighting lncRNA candidates. The associations between modules and GO terms (or ncRNAs) are established using Hypergeometric Test. The similarity among modules is quantified using semantic similarity. PCA is used to obtain the latent features of the input ncRNAs. Please refer to Section 2.1 for more details Fig. 1. View largeDownload slide Schematic diagram of MoonFinder. The orange, green and blue boxes represent the input, generated, and output data, respectively. The clustering algorithm and statistical models are shown in the dotted boxes. The inputs are the four datasets, i.e. interactions of proteins, localization annotation of proteins, function annotation of proteins, and interactions between ncRNAs and proteins, respectively. The output is a list of predicted moonlighting lncRNA candidates. The associations between modules and GO terms (or ncRNAs) are established using Hypergeometric Test. The similarity among modules is quantified using semantic similarity. PCA is used to obtain the latent features of the input ncRNAs. Please refer to Section 2.1 for more details 2.2 Statistical methods 2.2.1 Module-function and lncRNA–module association The hypergeometric test (HGT) was used to evaluate the statistical significance of the association between lncRNAs and functional modules. A lncRNA is considered to be interacting with a module if the proteins interacting with the lncRNA are significantly enriched in the module, i.e. P-value <0.05, which is shown as follows, p=1-∑i=0k-1ixm-in-xmn (2.1) where n is the total number of all proteins in the protein interaction network, m is the number of proteins analyzed in a module, x is the number of proteins interacting with a specific lncRNA, and i is the number of proteins in the module interacting with the lncRNA. The P-value describes the probability of randomly select no less than k interacted proteins in the module with size m. Similarly, the model was also utilized to carry out the functional annotation for the identified modules by functional categories (or GO terms), where in this case n is still the total number of gene products in the protein interaction network, but m represents the module size and x represents the term size. i is the number of gene products annotated in the term and involved in the module. The P-value shows the probability of randomly choosing no less than k proteins annotated by a GO term for a given module. Monte Carlo simulation was adopted to model the probability of interacting mlncRNAs and MPs. N pairs of mlncRNAs and MPs were randomly extracted from the lncRNA–protein interaction network and then we calculated the interacting pair number, Ti. The P-value is the ratio of the number of simulated interactions ( Ti) that is larger than the number of practical interactions ( T), which is mathematically defined as: p= ∑i=1Nsgn(Ti-T)N (2.2) where sgn is defined as sgnx=1, x≥00, x<0. To gain the statistical significance, all comparisons in this study between two lists of genes or gene products were analyzed using Wilcoxon Ranksum Test (WRT, R-3.4.1). 2.2.2 Decomposition of the functional similarity matrix Principal component analysis (PCA) is a statistical procedure used to reduce the number of features used to represent data. We accomplish the reduction by projecting data from a higher dimension to a lower dimensional manifold such that the error incurred by reconstructing the data in the higher dimension is minimized (Ma and Dai, 2011; Ma and Kosorok, 2009). Mathematically, we want to map the features x∈Rp to x∼∈Rq where q<p. Here, eigenvalue decomposition of the similarity matrix was used to calculate the principal components (PCs) and their weights for the modules interacted by each lncRNA. According to the Spectral Theorem, Avi⃗=λivi⃗ (2.3) we can calculate the eigenvalues λ1+⋯+λm of the Similarity matrix (arranged in decreasing order) and accordingly the trace of the matrix is T=λ1+⋯+λm (2.4) Only the lncRNAs whose similarity matrix can be decomposed into PCs without a dominant role were selected as the candidates, such as the sum of the weights of the first k PCs is less than a threshold τ ( τ=0.7 by default) as follows, argmaxk=0,1,2…∑i=1kλiT<τ (2.5) The threshold τ is a user-defined parameter. The higher the value, the more mlncRNA can be identified. Here, we set τ=0.7 because a PC will play a dominant role with a weight distribution coefficient less than 0.5. The weight distribution coefficient is defined as wdc=1-∑(λiT)2, which captures the distribution of the weights of PCs for a lncRNA. If the weights are extremely similar then the distribution coefficient is close to 1, while if some weights play a dominant role the distribution coefficient approaches 0. Then the maximum k is the number of latent features. The eigenvalue or the weight of a latent feature has minor influence when it is less than 0.1, so we also define λi≥0.1 in Eqt 2.4 and accordingly the number of latent features with key contribution is n=∑i=1msgn(λi-λcutoff) (2.6) where sgn is defined as sgnx=1, x≥00, x<0 and λcutoff=0.1. Consequently, the number of latent functions for an RNA is l= min(k,n) (2.7) An RNA is determined as moonlighting RNA if l is larger than one. Namely, the RNA has more than one latent feature in terms of the interacting functional modules. 3 Results 3.1 Experimental and parametric setups of MoonFinder Moonlighting lncRNAs (mlncRNAs) are assumed to execute multiple distinct functions through interactions with proteins that are localized in different cell compartments, so we identify mlncRNAs from the aspect of whether they are targeting multiple but functionally unrelated protein modules in a co-localized protein interactome. The human protein interactome was first refined and only the co-localized interactions were maintained for module identification, since proteins closely interacted with each other in a module are more likely to reside in the same cell compartment. Eventually, a compartment-specific protein interaction network with 210 410 interactions among 10 111 proteins was constructed. After that, a total of 765 functional modules were identified using ClusterONE, an algorithm that can detect overlapping clusters of proteins highly connected inside but sparse outside (Nepusz et al., 2012). We choose ClusterONE because not only it can detect biological relevant clusters that can be appropriately mapped to modules, but also its ability to softly identify the overlapping modules considering the network topological structure. Functional enrichment analysis was employed to annotate each identified module with specific and significant function categories. We applied the hypergeometric test (HGT) to obtain the enrichment P-values and the FDR adjusted P-values of 0.01 were eventually used as the threshold. To establish the RNA–module interactions, similarly, we used HGT to assess whether there are significant connections between the lncRNAs and the functional modules. The target proteins of a lncRNA are significantly overrepresented in a module if the FDR adjusted enrichment P-values were less than a given threshold of 0.01. Accordingly, we obtained a bipartite network with 2726 interactions among 538 lncRNA and 106 protein modules. The function similarity matrices were calculated using the semantic similarities among the modules of each lncRNA using five semantic similarity measurements, i.e. Resnik, Lin, Jiang and Conrath, Schlicker and Wang (see Supplementary Section S4). Only the intersection of the five sets of mlncRNAs identified using the five measures was determined as the candidates (in total 155), because one measure may outperform the others in different expression or interaction scenarios. Importantly, we utilized eigenvalue decomposition to calculate the number of latent features of each lncRNA. The lncRNAs whose module similarity matrix can be decomposed into the principal components without a dominant role (more than one latent feature) were selected as the candidates. The workflow of MoonFinder is described in more detail in the Materials and methods Section 2.2 as well as in Figure 1. 3.2 Overview of the mlncRNA candidates As shown in the Venn diagram of Figure 2A, among the 1284 lncRNAs with interacted proteins (background lncRNAs), 538 lncRNAs (flncRNAs) are annotated to at least one functional module and eventually 155 out of them are determined as moonlighting lncRNAs (mlncRNAs), whose target modules are functionally unrelated. The mlncRNAs were identified using five semantic similarity (SS) measures to quantify the functional similarity between modules. To obtain more reliable results, only lncRNAs in the intersection of the sets of identified mlncRNAs using the five different measures were defined as mlncRNAs (Fig. 2B). These identified mlncRNA candidates are displayed in Supplementary Table S1, including the number of interacting modules and the number of latent features for different SS measurements. The modules targeted by an identical lncRNA have a higher probability of sharing the same functions than the randomly picked modules. Not surprisingly, significantly more mlncRNAs are detected when simulating the randomly selected modules as the target modules. Specifically, Figure 2C shows that around 40% of the lncRNAs with target module can be identified as mlncRNAs using either SS measures, the ratio is as low as 28% when using Lin’s measure, while the ratios increase to about 60% when the target modules are randomly picked. Fig. 2. View largeDownload slide Overview of the identified mlncRNAs. (A) Venn diagram of the lncRNAs. (B) Venn diagram of the mlncRNAs identified using five semantic similarity measures. (C) The ratio of real and randomly identified mlncRNAs. (D) An example of mlncRNA ANCR. (E–G) Gene Ontology functional enrichment of the three modules regulated by ANCR Fig. 2. View largeDownload slide Overview of the identified mlncRNAs. (A) Venn diagram of the lncRNAs. (B) Venn diagram of the mlncRNAs identified using five semantic similarity measures. (C) The ratio of real and randomly identified mlncRNAs. (D) An example of mlncRNA ANCR. (E–G) Gene Ontology functional enrichment of the three modules regulated by ANCR Here, we take ANCR as an example to illustrate its moonlighting functions. ANCR is an anti-differentiation lncRNA that is required to enforce the undifferentiated state in somatic progenitor populations of epidermis. Using MoonFinder, we observed ANCR mainly interacts with three functional modules with closely linked proteins inside and no proteins are shared between any two modules, as shown in Figure 2D. Specifically, one module was enriched in a variety of metabolic processes, such as estrogen, retinal and steroid, etc., whereas another module was enriched in functions like tissue morphogenesis and development. Signaling pathways for enforcing intracellular receptors like toll-like receptor 5 and 10 were highly represented for the other module (Fig. 2E–G). Consequently, ANCR was determined as a mlncRNA who shows its functional diversity via regulating protein modules taking part in distinct biological processes and it would serve as a highly reliable candidate of moonlighting lncRNA. 3.3 Sequence features of mlncRNAs To investigate whether the mlncRNAs form a distinct group of lncRNAs, we analyzed the sequences of the corresponding non-coding genes to detect common features such as biotype, gene length, transcript length, transcript number, exon length and exon number, as well as the evolutionary conservation and expression correlation with orthologous genomes. The candidates were compared with other two groups of lncRNAs, i.e. the entire set of lncRNAs with protein targets (background lncRNAs) and the functional module related lncRNAs (flncNRAs), which interacts with functionally related or unrelated modules (Fig. 2A). Hence, it is important to identify the unique characteristics of mlncRNAs that the other types of lncRNAs do not exhibit. Although no significant differences of gene length were detected for distinct categories of lncRNAs (Fig. 3A), the candidate mlncRNAs have a significantly longer transcript than the other types of lncRNAs. On average, the transcript length is about 19 500 compared with about 11 000 for the flncRNAs (WRT P-value = 1.27e-3; Fig. 3B). But the number of transcripts in these lncRNA categories are similar to each other, only a marginal significance was tested between the mlncRNAs and the background lncRNAs (WRT P-value = 4.7e-2; Fig. 3C). Importantly, we observed that the number of exon of mlncRNAs is much larger than the regular lncRNAs (10 versus 1, WRT P-value = 2.2-e16; Fig. 3D), but these exons are significantly shorter than the other categories of lncRNAs (WRT P-value = 3.4e-6 and P-value = 8.3e-9; Fig. 3E), probably owing to the transcripts of mlncRNAs are on average much longer than that of the regular lncRNAs. In short, mlncRNAs have short but more exons, which is a potential sequence feature for lncRNAs to moonlight in between multiple biological functions. Additionally, Tilgner et al. observed that the splicing efficiency of an exon is negatively correlated with the length of the downstream intron, supposedly owing to the RNA-seq reads spanning the junction between exon and intron can only be detected when transcription occurs, which is a necessary condition for splicing to be carried out (Tilgner et al., 2012). So, we also explored if there is a difference in the intron length or splicing efficiency compare to background lncRNAs and flncRNAs. As shown in Figure 3F, unfortunately, no significant difference in the average intron length was observed among the three lncRNA categories. Fig. 3. View largeDownload slide The distributions of lncRNAs in different functional groups with regard to distinct sequence features. (A) Gene length. (B) Transcript length. (C) Transcript number. (D) Exon number. (E) Exon length. (F) Intron length. (G) Evolutionary conservation and expression correlation. Outliers are not shown. The bars from left to right correspond to lncRNA, flncRNA, and mlncRNA, respectively Fig. 3. View largeDownload slide The distributions of lncRNAs in different functional groups with regard to distinct sequence features. (A) Gene length. (B) Transcript length. (C) Transcript number. (D) Exon number. (E) Exon length. (F) Intron length. (G) Evolutionary conservation and expression correlation. Outliers are not shown. The bars from left to right correspond to lncRNA, flncRNA, and mlncRNA, respectively Phylogenetic conservation and expression correlation are strong evidence for inferring functions of coding or non-coding genes (Cheng et al., 2016a,b; Park et al., 2014). As lncRNAs are crucial in biological processes if they are evolutionarily conserved or expression-correlated across species, we checked whether the mlncRNAs tend to be conserved among orthologous genomes and whose expression patterns are highly correlated in orthologs. As shown in Figure 3G, the mlncRNAs are prone to be evolutionarily conserved (>12%) than flncRNAs and background lncRNAs (<10%). The ratios of flncRNAs are lower than mlncRNAs but still higher than that of the background lncRNAs. Similarly, mlncRNAs have the largest proportion of expression-conserved ones and the ratio of flncRNAs is second to it. Consequently, the evolution and expression patterns of mlncRNAs are more conserved than those of the other lncRNAs, suggesting that the lncRNAs moonlighting in the cells may play more important biological roles. Besides, RNA species were officially grouped into several biotypes by their transcriptional direction and exosome sensitivity (Hon et al., 2017). Here we also examined the relationship between the functional categories and biotypes, but no significant correlation was detected (Supplementary Fig. S1A). 3.4 Subcellular localization features Next, we aimed to understand how the mlncRNAs behave relative to the other lncRNAs in terms of subcellular localization. To investigate the spatial distribution of lncRNAs at a subcellular level, we applied Relative Concentration Index (RCI) (Mas-Ponte et al., 2017), a ratio of a transcript’s concentration between two cellular compartments, to measure the localization tendency of non-coding RNAs. Essentially, RCI is the log2 transformed ratio of FPKM (fragments per kilobase per million mapped) in two compartments like cytoplasm and nucleus. First, we calculated the cytoplasmic-nuclear RCI to measure the relative concentration of a lncRNA between the cytoplasm and the nucleus in 15 cell lines. Figure 4A illustrates the RCI distributions of mlncRNAs, flncRNAs, and background lncRNAs. It is apparent that mlncRNAs tend to have higher RCI values compared to the other two categories of lncRNAs in almost all these cell lines except SK.N.DZ, SK.MEL.5 and K562, suggesting a more variable localization pattern of mlncRNAs with much higher chance of residing in the cytoplasmic compartment in comparison to the other lncRNAs. Then, we further investigated the localization of mlncRNAs at the sub-compartment level, since LncATLAS also provides information about enrichment in the cytoplasmic and nuclear sub-compartments of the K562 cells. As shown in Figure 4B, the sub-nuclear RCI values of mlncRNAs are higher than that of the other two groups of lncRNAs while the sub-cytoplasmic RCI values are relatively small. Namely, in the nucleus, the mlncRNAs are prone to appear in the sub-compartments of nucleoplasm, nucleolus and chromatin, whereas in the cytoplasm, the mlncRNAs are not likely to reside in insoluble and membrane relative to the other lncRNAs. That is why the cytoplasmic-nuclear RCI values of mlncRNAs are almost the same as the background lncRNAs in the K562 cell line (Fig. 4A). Meanwhile, we also calculated the expression value distribution of lncRNAs in each sub-compartment in the K562 cell line. More importantly, the expression values of the mlncRNAs are significantly higher in all the sub-compartments, revealing that the expression abundance of lncRNAs is crucial in executing the part-time functions (Fig. 4C). Fig. 4. View largeDownload slide The mlncRNA localization and expression features. (A) The boxes from top to bottom correspond to lncRNA, flncRNA, and mlncRNA. (B) Sub-compartment expression value distribution of lncRNAs in the K562 cell line. (C) Sub-compartment RCI distribution of lncRNAs in the K562 cell line. (D) The ratios of different groups of lncRNAs residing in multiple locations. (E–G) The ratios of different groups of lncRNAs residing in the nucleus, the cytosol, and the cytoplasm, respectively. RCI, Relative Concentration Index Fig. 4. View largeDownload slide The mlncRNA localization and expression features. (A) The boxes from top to bottom correspond to lncRNA, flncRNA, and mlncRNA. (B) Sub-compartment expression value distribution of lncRNAs in the K562 cell line. (C) Sub-compartment RCI distribution of lncRNAs in the K562 cell line. (D) The ratios of different groups of lncRNAs residing in multiple locations. (E–G) The ratios of different groups of lncRNAs residing in the nucleus, the cytosol, and the cytoplasm, respectively. RCI, Relative Concentration Index In addition, we also used another RNA localization resource RNALocate, which contains manually-curated localizations classifications, to investigate the localization tendency of mlncRNAs. In this database, the lncRNAs were collected and annotated to different cell compartments, e.g. nucleus, cytosol and cytoplasm. We calculated the ratio of lncRNAs in these compartments separately for each category of lncRNAs. We found that mlncRNAs tend to appear in more than one compartment and localize in the cytoplasmic compartments such as cytosol and cytoplasm (Fig. 4D–G). The ratio of multilocation mlncRNAs is as high as 0.35, which is much higher than that of flncRNAs and the background lncRNAs (about 0.3 and 0.26, respectively; Fig. 4D). More importantly, mlncRNAs were found to be enriched in cytosol and cytoplasm with the ratios of 0.55 and 0.3 (Fig. 4F, G), respectively, whereas the ratio is comparable to the other categories of lncRNAs in the nucleus (Fig. 4E). Consequently, we can draw the same conclusion that mlncRNAs have a higher chance of residing in the cytoplasmic compartment and therefore have a more variable spatial distribution. 3.5 Topological features of mlncRNAs and its roles in cancers lncRNA functions through its interacting partners. Accumulating studies show that the multi-functionality of lncRNAs as interacting hubs with other molecules such as proteins, DNAs and RNAs. Apparently, the candidate mlncRNAs connect a significantly larger number of proteins and modules than the other lncRNAs according to the identification methodology. On average, the number of partner proteins is 36.1 for mlncRNA while less than 20 for the others (WRT P-value = 7.4e-8; Supplementary Fig. S1B). The number of the interacted module is around 6.8 for mlncRNA compared with 5 for the other lncRNAs (WRT P-value = 8.8e-14; Supplementary Fig. S1C). To illustrate the combinatorial regulation and give a systematic description, we constructed an mlncRNA–module regulatory network, in which the edges link the candidate mlncRNAs to their corresponding functional modules. This network contains in total 1055 predicted regulatory interactions between 155 mlncRNAs and 83 modules (Fig. 5A). Some modules connect more candidate mlncRNAs than others, indicating that they might be engaged in a larger number of moonlighting regulations. In particular, the largest module in the center of the network shows the highest degree of 70, suggesting that it is subjected to the regulations of 70 mlncRNAs. The proteins in this module were found to be mainly implicated in biological processes such as nucleic acid metabolic process, gene expression and amniotic stem cell differentiation (Fig. 5B). Fig. 5. View largeDownload slide Association between mlncRNAs and diseases. (A) The mlncRNA–module regulation network. RNAs are represented in squares while modules are in circles. The bars from left to right correspond to lncRNA, flncRNA, and mlncRNA. The circle size corresponds to the module size. (B) Gene Ontology function enrichment of the module with the largest number of regulated mlncRNAs. (C, D) The ratio of cancer and disease lncRNAs among the three lncRNA categories. (E, F) The ratio of lncRNAs associated with cancer and drug target module. (G) mlncRNAs tend not to interact with MPs. T, the number of practical mlncRNA-MP interactions. (H) mlncRNAs and MPs tend to share less interacting partners. T, the number of common partners practically interacted by mlncRNAs and MPs. (I) Summary of the top ten mlncRNAs with the highest OR scores Fig. 5. View largeDownload slide Association between mlncRNAs and diseases. (A) The mlncRNA–module regulation network. RNAs are represented in squares while modules are in circles. The bars from left to right correspond to lncRNA, flncRNA, and mlncRNA. The circle size corresponds to the module size. (B) Gene Ontology function enrichment of the module with the largest number of regulated mlncRNAs. (C, D) The ratio of cancer and disease lncRNAs among the three lncRNA categories. (E, F) The ratio of lncRNAs associated with cancer and drug target module. (G) mlncRNAs tend not to interact with MPs. T, the number of practical mlncRNA-MP interactions. (H) mlncRNAs and MPs tend to share less interacting partners. T, the number of common partners practically interacted by mlncRNAs and MPs. (I) Summary of the top ten mlncRNAs with the highest OR scores To determine whether the moonlighting of lncRNAs is implicated in the formation and development of cancers and other diseases, we associated the mlncRNAs with public available cancer and disease lncRNAs as well as cancer proteins and drug target proteins (see Materials and methods). Around 13% of the mlncRNA candidates have been studied involved in cancer processes, while the ratios of flncRNAs and background lncRNAs are less than 10% (Fig. 5C). When considering other complex diseases, not surprisingly, the ratio is as high as 23% for mlncRNAs, which is still much higher than that the other categories (16 and 17%, respectively; Fig. 5D). From the perspective of regulated functional modules, about 39% of the candidates are included in modules significantly enriching cancer proteins, whereas the ratios decrease to about 29 and 33% for the other two groups (Fig. 5E). Similarly, for the modules enriched of drug targets, they are more likely to interact with mlncRNAs than the flncRNAs and the background lncRNAs (12% versus 7% versus 8%; Fig. 5F). Besides, we can draw the same conclusion when concentrating on the proportion of cancer modules or drug target modules regulated by the mlncRNAs (Supplementary Fig. S1D and E). Therefore, these results demonstrate that the mlncRNAs exercise a great influence on cancer metastasis and progression through pairwise interactions and clustered modules of proteins in the regulatory network. More importantly, we observed that the mlncRNAs and moonlighting proteins (MPs) tend to be mutually exclusive in terms of interactions and interacted partners. Only five MPs directly interacts with the mlncRNAs, whereas on average the value is as high as 18.9 when randomly selected the same number of MPs (Monte Carlo P-value = 1e-04, HGT P-value = 5.4e-18, Fig. 5G). Also, we simulated the interacted partners of both mlncRNAs and MPs 10 000 times and found the number of common partners between them is 691 on average, which is significantly higher than the practical value of 529 (Monte Carlo P-value = 0, HGT P-value = 2.1e-86, Fig. 5H). In other words, the number of common partner proteins that the moonlighting lncRNAs and proteins shared is significantly less than that of randomly selected ones. These results indicate the possible mechanism that the cells make full use of the macromolecules to efficiently and systematically perform cellular tasks avoiding the redundant implementations. Additionally, from the mlncRNA–module network in Figure 5A, we found the mlncRNAs that exclusively interact functional modules tend to be cancer-related. Accordingly, we introduced a score, Interactor Share Rate (ISR), to measure how likely the interactors of a given lncRNA are shared by the other lncRNAs (see Supplementary Section S11). We found that the cancer mlncRNAs have significantly higher ISRs than that of the others (WRT P-value = 3.2e-03, Supplementary Fig. S1F). For the mlncRNAs with the top ten highest ISR scores, six out of them are cancer lncRNAs (Fig. 5I). When strengthening the threshold to 0.5, six out of eight (75%) of the mlncRNAs are cancer genes and the other two, ANCR and LRRC75A-AS1, could be considered as the candidates of cancer mlncRNA, where the dysfunction or inappropriate switching of these RNAs in different cell compartments may result in the biological activity of cancer, although further experimental works are needed to warrant this claim. 4 Discussion In this study, we introduced a computational framework MoonFinder to systematically identify moonlighting lncRNAs (mlncRNAs) based on the integrated lncRNA and protein interaction network as well as the protein functional annotations. In total 155 lncRNAs were determined as candidates with multiple but distinct functions. Also, we characterized them from various aspects of sequence features, evolutionary conservation, expression correlation, expression abundance, localization tendency and interaction patterns, which will facilitate our further understanding of their functions and the mechanism of moonlighting. Remarkably, we observed that the non-coding genes that transcript mlncRNAs tend to have shorter but more exons, which is a potential sequence feature for lncRNAs to moonlight in between multiple biological functions. Also, we found the evolution and expression patterns of mlncRNAs are more conserved than the other lncRNAs, which in contrast with the conventional knowledge that lncRNAs are generally less conserved than mRNAs and proteins (Hon et al., 2017; Park et al., 2014), suggesting that mlncRNAs are central for the homeostasis maintenance of human. More importantly, we found that mlncRNAs have a more variable localization pattern with a high chance of residing in the cytoplasmic compartment, although they display high expression across all the cell compartments. mlncRNAs are expressed significantly higher in all the sub-compartments of the K562 cell lines in comparison with the other lncRNAs, suggesting that the high expression abundance is necessary for executing the part-time functions. We studied the localization tendency and translocation activity of these mlncRNAs because lncRNAs are diversely resided in the cells and play a crucial role as modulators to regulate gene expression in multiple ways (Cabili et al., 2015; Ferre et al., 2016; Quinn and Chang, 2016; Zhou et al., 2017; Zhu et al., 2016). It is well known that the competing endogenous RNAs (ceRNAs) can protect the coding RNAs from repression by binding them to the co-targeted miRNAs. mlncRNAs may function diversely in the same way by competing or shared miRNAs, which mainly reside in the cytoplasm to execute their functions, leading to a more variable spatial distribution (Rashid et al., 2016). In the future, we will investigate whether the mlncRNAs act as ceRNAs to modulate their expression abundance and interaction pattern, e.g. regulating the abundance of lncRNAs in one compartment may influence the function of the other cell compartment. Our result also shows that mlncRNAs and MPs are rather mutually exclusive in terms of their direct interactions and interacting partners. In other words, lncRNAs and proteins with moonlighting functions are not likely to interact with each other and they even tend to share fewer neighbors in the regulatory network. The reason might be that the macromolecules in cells are usually organized to be efficient to perform different cellular tasks without redundancy. According to the mlncRNA–module bipartite network, we also predicted eight cancer lncRNAs and six out of them were previously identified as cancer lncRNAs by different experimental assays. A potential limitation is that the low value of semantic similarity between a pair of modules might due to the incompleteness of GO. For instance, one lncRNA is detected as a moonlighting candidate and it regulates some modules with low semantic similarity, which may partially due to the proteins involved in these modules show no evidence of association with GO terms, but the association might have not been established yet. It is noted that the absence of annotation does not imply the absence of function, because the information of GO is still far from complete, although the creation of new annotations increased dramatically. Another potential bias is that the GO is structured as a graph with unbalance of high and low levels of details among branches (Pritykin et al., 2015). This has a strong implication on measuring the semantic similarity of GO terms. For instance, sibling GO terms can be semantically very similar or extremely distinct in different branches of the GO structure. In this study, we leveraged semantic similarity among the GO terms to identify moonlighting lncRNAs. Since similar processes may be categorized in distinct branches in GO, low values of semantic similarity do not always imply multifunctionality and hence our findings may include some false positive discoveries. Pritykin et al. proposed an approach to identify multifunctional genes based on GO, which defined a multifunctional gene as it is annotated with at least ‘two sufficiently distinct functional terms of comparable specificity’ (Pritykin et al., 2015). The definition of multifunctionality might help mitigate the problem of false positive discovery in detecting moonlighting lncRNAs. In our future work, we will integrate this co-occurrence based definition of multifunctionality into MoonFinder. We believe our observations can aid our and other research groups to understand how they function in a moonlighting manner and help in designing RNAs with novel functions and applications. Moreover, investigating the mechanisms that determine the functional diversity of mlncRNAs has the potential to provide new insights into their biogenesis, physical interaction, subcellular localization and therapeutic application in diseases. In the future, we will investigate the mechanism of how the mlncRNAs modulate and switch the functions in metabolic processes, which is of vital importance for cancer therapeutics and will provide tremendous opportunities for anti-cancer strategies. The moonlighting feature of the other types of RNAs, such as miRNA and circRNA (Chen, 2016a,b), will also be studied and compared and eventually a moonlighting atlas of both RNAs and proteins will be provided. Similar to animal cells, plant cells are heavily compartmentalized. Besides the typical organelles in animal cells, plant cells usually contain a large number of plastids and a central vacuole, which is critical in storing organelles and involved in the regulation of cellular osmoregulation (Gao et al., 2015; Zeng et al., 2015; Zhuang et al., 2017). Hence, we will explore and establish methodologies to identify moonlighting ncRNAs and proteins for in-depth functional studies in plants. Funding This work was supported by The Chinese University of Hong Kong Direct Grant; and the Research Grants Council of Hong Kong GRF Grant . Conflict of Interest: none declared. References Abumrad N.A. , Lange A.J. ( 2006 ) The metabolism of cancer cells: moonlighting proteins and growth control . Curr. Opin. Clin. Nutrition Metab. Care , 9 , 337 – 338 . Google Scholar Crossref Search ADS Barabasi A.L. , Oltvai Z.N. ( 2004 ) Network biology: understanding the cell's functional organization . Nat. Rev. Genet ., 5 , 101 – 113 . Google Scholar Crossref Search ADS PubMed Boukouris A.E. et al. ( 2016 ) Metabolic enzymes moonlighting in the nucleus: metabolic regulation of gene transcription . Trends Biochem. Sci ., 41 , 712 – 730 . Google Scholar Crossref Search ADS PubMed Buxbaum A.R. et al. ( 2015 ) In the right place at the right time: visualizing and understanding mRNA localization . Nat. Rev. Mol. Cell Biol ., 16 , 95 – 109 . Google Scholar Crossref Search ADS PubMed Cabili M.N. et al. ( 2015 ) Localization and abundance analysis of human lncRNAs at single-cell and single-molecule resolution . Genome Biol ., 16 , 20. Google Scholar Crossref Search ADS PubMed Chapple C.E. et al. ( 2015 ) Extreme multifunctional proteins identified from a human protein interaction network . Nat. Commun ., 6 , 7412 . Google Scholar Crossref Search ADS PubMed Chen L.L. ( 2016a ) The biogenesis and emerging roles of circular RNAs . Nat. Rev. Mol. Cell Biol ., 17 , 205 – 211 . Google Scholar Crossref Search ADS Chen L.L. ( 2016b ) Linking long noncoding RNA localization and function . Trends Biochem. Sci ., 41 , 761 – 772 . Google Scholar Crossref Search ADS Cheng L. et al. ( 2016a ) CrossNorm: a novel normalization strategy for microarray data in cancers . Sci. Rep ., 6 , 18898 . Google Scholar Crossref Search ADS Cheng L. et al. ( 2016b ) ICN: a normalization method for gene expression data considering the over-expression of informative genes . Mol. Biosyst ., 12 , 3057 – 3066 . Google Scholar Crossref Search ADS Du Z. et al. ( 2013 ) Integrative genomic analyses reveal clinically relevant long noncoding RNAs in human cancer . Nat. Struct. Mol. Biol ., 20 , 908 – 913 . Google Scholar Crossref Search ADS PubMed Espinosa-Cantu A. et al. ( 2015 ) Gene duplication and the evolution of moonlighting proteins . Front. Genet ., 6 , 227. Google Scholar Crossref Search ADS PubMed Ferre F. et al. ( 2016 ) Revealing protein–lncRNA interaction . Brief. Bioinform ., 17 , 106 – 116 . Google Scholar Crossref Search ADS PubMed Gao C. et al. ( 2015 ) Dual roles of an Arabidopsis ESCRT component FREE1 in regulating vacuolar protein transport and autophagic degradation . Proc. Natl. Acad. Sci. USA , 112 , 1886 – 1891 . Google Scholar Crossref Search ADS Hon C.C. et al. ( 2017 ) An atlas of human long non-coding RNAs with accurate 5' ends . Nature , 543 , 199 – 204 . Google Scholar Crossref Search ADS PubMed Jeffery C.J. ( 2015 ) Why study moonlighting proteins? Front. Genet ., 6 , 211. Google Scholar Crossref Search ADS PubMed Khan I. et al. ( 2014 ) Genome-scale identification and characterization of moonlighting proteins . Biol. Direct , 9 , 30 . Google Scholar Crossref Search ADS PubMed Khan I.K. et al. ( 2017 ) DextMP: deep dive into text for predicting moonlighting proteins . Bioinformatics , 33 , i83 – i91 . Google Scholar Crossref Search ADS PubMed Khan I.K. , Kihara D. ( 2016 ) Genome-scale prediction of moonlighting proteins using diverse protein association information . Bioinformatics , 32 , 2281 – 2288 . Google Scholar Crossref Search ADS PubMed Liao Q. et al. ( 2011 ) Large-scale prediction of long non-coding RNA functions in a coding-non-coding gene co-expression network . Nucleic Acids Res ., 39 , 3864 – 3878 . Google Scholar Crossref Search ADS PubMed Liu W. et al. ( 2014 ) Gene co-expression analysis identifies common modules related to prognosis and drug resistance in cancer cell lines . Int. J. Cancer , 135 , 2795 – 2803 . Google Scholar Crossref Search ADS PubMed Ma S. , Dai Y. ( 2011 ) Principal component analysis based methods in bioinformatics studies . Brief. Bioinform ., 12 , 714 – 722 . Google Scholar Crossref Search ADS PubMed Ma S. , Kosorok M.R. ( 2009 ) Identification of differential gene pathways with principal component analysis . Bioinformatics , 25 , 882 – 889 . Google Scholar Crossref Search ADS PubMed Mani M. et al. ( 2015 ) MoonProt: a database for proteins that are known to moonlight . Nucleic Acids Res ., 43 , D277 – D282 . Google Scholar Crossref Search ADS PubMed Mas-Ponte D. et al. ( 2017 ) LncATLAS database for subcellular localization of long noncoding RNAs . RNA , 23 , 1080 – 1087 . Google Scholar Crossref Search ADS PubMed Min K.W. et al. ( 2016 ) Moonlighting proteins in cancer . Cancer Lett ., 370 , 108 – 116 . Google Scholar Crossref Search ADS PubMed Monaghan R.M. , Whitmarsh A.J. ( 2015 ) Mitochondrial proteins moonlighting in the nucleus . Trends Biochem. Sci ., 40 , 728 – 735 . Google Scholar Crossref Search ADS PubMed Nepusz T. et al. ( 2012 ) Detecting overlapping protein complexes in protein–protein interaction networks . Nat. Methods , 9 , 471 – 472 . Google Scholar Crossref Search ADS PubMed Ning S. et al. ( 2016 ) Lnc2Cancer: a manually curated database of experimentally supported lncRNAs associated with various human cancers . Nucleic Acids Res ., 44 , D980 – D985 . Google Scholar Crossref Search ADS PubMed Park C. et al. ( 2014 ) lncRNAtor: a comprehensive resource for functional investigation of long non-coding RNAs . Bioinformatics , 30 , 2480 – 2485 . Google Scholar Crossref Search ADS PubMed Piatigorsky J. et al. ( 1988 ) Gene sharing by delta-crystallin and argininosuccinate lyase . Proc. Natl. Acad. Sci. USA , 85 , 3479 – 3483 . Google Scholar Crossref Search ADS Piatigorsky J. , Wistow G.J. ( 1989 ) Enzyme/crystallins: gene sharing as an evolutionary strategy . Cell , 57 , 197 – 199 . Google Scholar Crossref Search ADS PubMed Pritykin Y. et al. ( 2015 ) Genome-wide detection and analysis of multifunctional genes . PLoS Comput. Biol ., 11 , e1004467. Google Scholar Crossref Search ADS PubMed Quinn J.J. , Chang H.Y. ( 2016 ) Unique features of long non-coding RNA biogenesis and function . Nat. Rev. Genet ., 17 , 47 – 62 . Google Scholar Crossref Search ADS PubMed Rashid F. et al. ( 2016 ) Long non-coding RNAs in the cytoplasm . Genomics Proteomics Bioinform ., 14 , 73 – 80 . Google Scholar Crossref Search ADS Sriram G. et al. ( 2005 ) Single-gene disorders: what role could moonlighting enzymes play? Am. J. Hum. Genet ., 76 , 911 – 924 . Google Scholar Crossref Search ADS PubMed Tilgner H. et al. ( 2012 ) Deep sequencing of subcellular RNA fractions shows splicing to be predominantly co-transcriptional in the human genome but inefficient for lncRNAs . Genome Res ., 22 , 1616 – 1625 . Google Scholar Crossref Search ADS PubMed Wahlestedt C. ( 2013 ) Targeting long non-coding RNA to therapeutically upregulate gene expression . Nat. Rev. Drug Discov ., 12 , 433 – 446 . Google Scholar Crossref Search ADS PubMed Wang P. et al. ( 2015 ) Identification of lncRNA-associated competing triplets reveals global patterns and prognostic markers for cancer . Nucleic Acids Res ., 43 , 3478 – 3489 . Google Scholar Crossref Search ADS PubMed Zeng Y. et al. ( 2015 ) Unique COPII component AtSar1a/AtSec23a pair is required for the distinct function of protein ER export in Arabidopsis thaliana . Proc. Natl. Acad. Sci. USA , 112 , 14360 – 14365 . Google Scholar Crossref Search ADS Zhou J. et al. ( 2017 ) LncFunNet: an integrated computational framework for identification of functional long noncoding RNAs in mouse skeletal muscle cells . Nucleic Acids Res ., 45 , e108 . Google Scholar Crossref Search ADS PubMed Zhu X. et al. ( 2016 ) A long non-coding RNA signature to improve prognosis prediction of gastric cancer . Mol. Cancer , 15 , 60 . Google Scholar Crossref Search ADS PubMed Zhuang X. et al. ( 2017 ) ATG9 regulates autophagosome progression from the endoplasmic reticulum in Arabidopsis . Proc. Natl. Acad. Sci. USA , 114 , E426 – E435 . Google Scholar Crossref Search ADS © The Author(s) 2018. Published by Oxford University Press. All rights reserved. For permissions, please e-mail: email@example.com This article is published and distributed under the terms of the Oxford University Press, Standard Journals Publication Model (https://academic.oup.com/journals/pages/open_access/funder_policies/chorus/standard_publication_model)
Bioinformatics – Oxford University Press
Published: Oct 15, 2018
Access the full text.
Sign up today, get DeepDyve free for 14 days.