
==== Front
Front Mol Biosci
Front Mol Biosci
Front. Mol. Biosci.
Frontiers in Molecular Biosciences
2296-889X
Frontiers Media S.A.

1452142
10.3389/fmolb.2024.1452142
Molecular Biosciences
Original Research
LncSTPred: a predictive model of lncRNA subcellular localization and decipherment of the biological determinants influencing localization
Hu et al.
10.3389/fmolb.2024.1452142
Hu Si-Le 1
Chen Ying-Li 1
Zhang Lu-Qiang 1

Bai Hui 1
Yang Jia-Hong 1
Li Qian-Zhong 1 2 *

1 School of Physical Science and Technology, Inner Mongolia University, Hohhot, China
2 The State Key Laboratory of Reproductive Regulation and Breeding of Grassland Livestock, Inner Mongolia University, Hohhot, China
Edited by: Vikram Dalal, Washington University in St. Louis, United States

Reviewed by: Bhavika Mam, Independent Researcher, Palo Alto, CA, United States

Preeti Verma, University of Virginia, United States

Gunjan Saini, Purdue University, United States

*Correspondence: Qian-Zhong Li, qzli@imu.edu.cn
05 9 2024
2024
11 145214220 6 2024
19 8 2024
Copyright © 2024 Hu, Chen, Zhang, Bai, Yang and Li.
2024
Hu, Chen, Zhang, Bai, Yang and Li
https://creativecommons.org/licenses/by/4.0/ This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
Introduction

Long non-coding RNAs (lncRNAs) play crucial roles in genetic markers, genome rearrangement, chromatin modifications, and other biological processes. Increasing evidence suggests that lncRNA functions are closely related to their subcellular localization. However, the distribution of lncRNAs in different subcellular localizations is imbalanced. The number of lncRNAs located in the nucleus is more than ten times that in the exosome.

Methods

In this study, we propose a new oversampling method to construct a predictive dataset and develop a predictive model called LncSTPred. This model improves the Adaboost algorithm for subcellular localization prediction using 3-mer, 3-RF sequence, and minimum free energy structure features.

Results and Discussion

By using our improved Adaboost algorithm, better prediction accuracy for lncRNA subcellular localization was obtained. In addition, we evaluated feature importance by using the F-score and analyzed the influence of highly relevant features on lncRNAs. Our study shows that the ANA features may be a key factor for predicting lncRNA subcellular localization, which correlates with the composition of stems and loops in the secondary structure of lncRNAs.

lncRNA
subcellular localization
oversampling method
algorithm improvements
model biological interpretation
National Natural Science Foundation of China 10.13039/501100001809 Nos. 32160216 Nos. 62361047 Nos. 62161033 The author(s) declare that financial support was received for the research, authorship, and/or publication of this article. This work was supported by the National Natural Science Foundation of China (Nos. 32160216, 62361047, 62161033).section-at-acceptanceBiological Modeling and Simulation
==== Body
pmc1 Introduction

Messenger RNAs (mRNAs) encode proteins that underlie various organismal phenomena, although they only represent about 2% of the total RNA. The remaining 98% consists of non-coding RNAs (ncRNAs), whose functions are still poorly understood (Birney et al., 2007; Wang and Li, 2013). Studies suggest that ncRNAs, especially long non-coding RNAs (lncRNAs) exceeding 200 base pairs in length, play vital roles in regulating biological activities (Atianand et al., 2017; Batista and Chang, 2013). The research on lncRNAs has expanded beyond their traditional biological functions to disease studies, particularly in cancer development and diagnosis (Esguerra and Eliasson, 2014; Flynn and Chang, 2014; Kameswaran and Kaestner, 2014; Li et al., 2014; Wang and Chang, 2011; Yan et al., 2015). Some lncRNAs are considered biomarkers for carcinogenesis and show differential expression between cancer and normal tissues (Harries, 2012; Kitagawa et al., 2012). For instance, overexpression of H19 can accelerate bladder and prostate cancer metastasis (Luo et al., 2013; Zhu et al., 2014a; Zhu et al., 2014b), while AFAP1-AS1 is highly expressed in esophageal adenocarcinoma (Wu et al., 2013), and Gas5 is underexpressed in breast cancer (Mourtada-Maarabouni et al., 2009; Wu et al., 2013). Even the Oncotype Dx genes has been used as a strong evidence for the risk classification in NCCN clinical guidelines (Sparano et al., 2018).

Based on the relative abundance of lncRNAs (Derrien et al., 2012), 17% of lncRNAs are enriched in the nucleus and 4% in the cytoplasm, each with distinct functions. In the nucleus, lncRNAs can act as molecular scaffolds (Clemson et al., 2009), assist in alternative splicing (Gonzalez et al., 2015), regulate chromatin remodeling (Jiang et al., 2015; Kugel and Goodrich, 2012; Martens et al., 2004; Melé and Rinn, 2016; Saxena and Carninci, 2011), and modify DNA/RNA methylation (Yang et al., 2015). In the cytoplasm, they can regulate translation, promote or inhibit mRNA degradation (Gong and Maquat, 2011; Wilusz, 2016), and affect gene expression by binding to miRNAs (Lauressergues et al., 2015; Paraskevopoulou et al., 2013; Wang et al., 2010; K; Wang et al., 2015; Winter et al., 2009). Therefore, studying the subcellular localization and functions of lncRNAs is crucial (Bridges et al., 2021; Miao et al., 2019). Although experimental methods like lncRNA-FISH have been widely used, their utilization is hindered by their time-consuming nature, high cost, and low efficiency (L. Wang et al., 2015; Xiao et al., 2015). Consequently, many researchers have directed their efforts toward developing more reliable and efficient predictive models to address these shortcomings.

In recent decades, several predictive models for lncRNA subcellular localization have been proposed. For instance, lncLocation incorporates sequential, physicochemical, and structural features in its predictive model construction. It leverages the SVM algorithm along with binomial distribution and iterative feature selection techniques to create predictive models (Feng et al., 2020). In iLoc-lncRNA 2.0, the predictive model is built using 8-mer features and the mRMR method, which results in 1407 features and subsequently submitted to the SVM algorithms for model construction (Zhang et al., 2022). DeepLncLoc transforms sequence information into a matrix using word2vec, then uses a CNN to construct the DeepLncLoc model (Zeng et al., 2022).

In this paper, we presented the LncSTPred, an Adaboost-based model for predicting and interpreting lncRNA subcellular localization based on primary RNA sequences in the RNALocate database. Our model supports five localization types including nucleus, cytoplasm, cytosol, ribosome, and exosome, accommodating both sequence and structure features. To solve the imbalance of categories, we employed Borderline-SMOTE, ADASYN, and UNCERTAIN WEIGHT on the training set. We enhanced Adaboost using maximum likelihood estimation and selected key features through the F-score. Subsequently, we conducted bioinformatics analyses of sequence and structure distributions of these features. An overview of the research process is depicted in Figure 1.

FIGURE 1 The flowchart of predictive model construction.

2 Materials and methods

2.1 Collection and preprocessing of dataset

The data used in this study related to subcellular localization of LncRNA in mammals, including Homo sapiens and Mus musculus, and were extracted from the RNALocate database (https://www.rnalocate.org/download) (Cui et al., 2022). Researchers can download all raw lncRNAs in the “Download and API” page. Subsequently, we retained lncRNAs related to H. sapiens and M. musculus, and excluded lncRNAs without sequence annotation information and those with multiple subcellular localizations. Redundant sequences were then removed using the CD-HIT program (Huang et al., 2010; Li and Godzik, 2006) with a threshold set at 80%. Subsequently, categories with sample sizes below 10 were excluded, resulting in 1342 lncRNAs across five distinct categories. These categories included 673 lncRNAs located in the nucleus, 407 in the cytoplasm, 152 associated with ribosomal localization, 94 within the cytosol, and 16 originating from exosomes, as depicted in Table 1.

TABLE 1 Number of lncRNA in each subcellular localization.

Subcellular localization	Mammals	Homo sapiens	Li’s Dataset Zeng et al. (2022)	Lin’s Dataset Su et al. (2018)	
Nucleus	673	342	325	156	
Cytoplasm	407	71	328	426	
Cytosol	94	69	88		
Ribosome	152	132	88	43	
Exosome	16	16	28	30	

2.2 Nucleotide composition features

The nucleotide compositional features from lncRNA are significant features used to characterize the biological function of RNA and their species. Sequences are represented by Equation 1. The traditional method of sequence feature extraction is K-mer. Additionally, we aimed to extract sufficient information from the sequences. Therefore, we used reading frame (RF) features to further characterize the sequences, following the methodology outlined by Rainey et al. (Rainey and Repka, 2013). For convenience of description, we defined the 3-RF by Equation 2. Sequence=N1,N2,N3,N4,N5,…,NM‐4,NM‐3,NM‐2,NM‐1,NM (1)

3−RF1y=N1⁢,N2,N3,N4,N5,N6…,NM−5,NM−4⁢,NM−3,NM−2,NM−1,NM

3−RF2y=N1,N2,N3,N4,N5,N6,N7,…,NM−4,NM−3,NM−2,NM−1,NM (2)

3−RF3y=N1,N2,N3,N4,N5,N6,N7,N8,…,NM‐3⁢,NM‐2,NM‐1,NM

Where the sequence described a lncRNA with a length of M base pairs. N i denoted the type of nucleotide at position i. {N 1 , N 2 , N 3 , … , N K} was called k-mer, and there were 4k combinations in total. 3−RFxy represented the combination of three reading frames in different starting points. x = {1,2,3} represented the first, second and third position respectively, y = {1,2, 3, 64} were the 64 combinations of 3-mer respectively.

2.3 Minimum free energy

The formation of base pairs can reduce the energy of RNA molecules and make the structure more stable. Therefore, based on the core idea of the minimum free energy (MFE) and the Zuker algorithm, we defined the global minimum value of the overall energy as Equation 3 (Zuker and Stiegler, 1981; Zuker and Sankoff, 1984). Eij=Ei+1,j−1+αij,⁡min⁡Ei+k+βk,⁡min⁡Ei+k,j−l+γk+l,⁡min⁡Ei+k,j+Ei,j−l+εk+j+i−l,δj−l (3)

Where, αij represents the stacking energy when i and j are paired. βk,γk,εk,and δk describe the energy of the bulge loop, interior loop, multi-branched loop, and hairpin loop, respectively. In the actual calculation, Zuker’s algorithm uses four free energy functions and five dynamic programming matrices. The minimum free energy of the RNA sequence is similar to the backtracking process of the Nussinov base pair maximization algorithm.

2.4 Feature importance

The F-score is a simple and effective feature selection method which measures the discriminative power of features across categories (Xie et al., 2010). The F-score of the ith feature in a multi-classification problem can be defined as Equation 4: Fi=∑j=1lx¯ij−x¯i2∑jl1nj−1∑k=1njxk,ij−x¯i2 (4)

Where l denotes the total number of categories. nj represents the number of samples in the jth class. The x¯i and x¯ij are the average values of the ith feature across the entire dataset and the jth data subset, respectively. xk,ij is the kth observation of the ith feature in the jth class

2.5 Sample balance process

Machine learning algorithms generally assume that positive and negative datasets are balanced. However, when the ratio of positive to negative sets exceeds 1:3, the results will be affected. To address this imbalance, two main approaches can be employed. One is adjusting sample weights in the algorithm, such as using the class-weight parameter in the XGBoost algorithm to balance different categories. The other is oversampling the minority categories to equalize sample numbers across categories. In this study, the ratio of nuclear to ribosome samples exceeds 1:4, and the ratio of nuclear to exosome samples is more than 1:40, making it challenging to find a suitable class-weight to characterize this complex distribution. Therefore, we constructed a predictive dataset using two oversampling methods, Borderline-SMOTE (Han et al., 2005) and ADASYN (He et al., 2008), for sample balancing. By focusing on borderline instances, Borderline-SMOTE generally produces better classification results compared to SMOTE, particularly when the minority class is at high risk of misclassification. It also reduces the risk of overfitting by concentrating on the most informative samples. ADASYN focuses on the more difficult minority class samples. This targeted approach ensures that the classifier is better trained on challenging examples, improving its robustness and generalizability. Finally, we filtered with the UNCERTAIN WEIGHT (Kendall et al., 2018) method.

2.6 Predictive model algorithm

The Boosting algorithm constructs high-accuracy classifiers by combining several base classifiers, each with moderate accuracy. Adaboost exemplifies this strategy and is known for its high accuracy and ability to model complex split interfaces through nonlinear combination. In the Adaboost algorithm, each base classifier generates a predicted classification result and a self-correction factor to estimate the reliability of the classification (Schapire and Singer, 1999).

In the original binary classification problem solved in the Adaboost algorithm, the coefficients αt corrected for the basic classifier during iteration are shown in Equation 5: αt=12log1−εtεt (5)

εt is the classification error rate of the base classifier on the training dataset. This coefficient ensures that, in each round, the classifier’s accuracy is at least greater than random probability, which is more than 1/2 in a binary classification problem. To extend this to a multi-classification problem, after ensuring the training data is balanced, the accuracy of the base classifier must be at least 1/k, where k is the number of categories. In this paper, k = 5. To ensure that each round prioritizes minimizing the classification error of the base classifier with the highest weight in the final classifier, we refined it as Equation 6: αt=12log1−εtεt+⁡logk−1 (6)

In multi-classification Adaboost, we update the sample weights and decrease the weight of the previously classified base classifier. The weight of correctly classified samples is shown as Equation 7: ωt+1,i′=ωt,i×⁡exp−αt (7)

The weight of wrongly classified samples is shown as Equation 8: ωt+1,i′=ωt,i×⁡expαt (8)

To ensure the weights sum to 1, the weights ωt+1,i of round t+1 in the training set after the round t are shown as Equation 9: ωt+1,i=ωt+1,i′∑iNωt+1,i′ (9)

2.7 Performance evaluation

We use the Specificity (Sp), Sensitivity (Sn), Accuracy (ACC), and Matthews Correlation Coefficient (MCC) to measure the performance of the predictive model. Evaluation indicators can be written as Equation 10: Sn=TPTP+FN ; Sp=TNTN+FPACC=TP+FNTP+TN+FP+FN

MCC=TP×TN−FP×FNTP+FPTP+FNTN+FPTN+FN (10)

In the context of classification issue, TP represents the number of correctly recognized positives, FN represents the number of positives recognized as negatives, FP represents the number of negatives recognized as positives, and TN represents the number of correctly recognized negatives. Additionally, the ROC (Receiver Operating Characteristic) curve is established to evaluate the model’s robustness. The AUC value, ranging from 0 to 1, represents the area under the ROC curve. A larger AUC value indicates better model performance.

3 Results

3.1 Description of predictive dataset

In our original dataset, the ratio of ribosome samples, cytosol samples, and exosome samples to nucleus samples exceeded 1:3, the distribution of original dataset is shown in Figure 2A. To address this imbalance, we oversampled cytoplasm, ribosomes, and exosomes. The balanced dataset was analyzed using the K-means clustering method with all features, including k-mer, k-RF, and MFE. Figure 2B illustrates the K-means clustering results after data balancing. Clusters three and five corresponded to nuclear and cytoplasmic localization, while clusters 1, 2, and four represented ribosome, cytosol, and exosome localizations, respectively. After preprocessing, which involved removing outliers and de-linearizing the dataset, the revised K-means clustering results were shown in Figure 2C. In this updated dataset, clusters 1 to five denoted ribosome, cytoplasm, exosome, cytosol, and nuclear localizations, respectively, the processing of the dataset is shown in Table 2.

FIGURE 2 Demonstration of sample distribution at various stages of the balancing process and predictive model results (A). The distribution of the original dataset through K-means clustering. X-axis and Y-axis were projection of the data onto the first dimension of a reduced-dimensional space after using PCA (Principal Component Analysis) (B) The distribution of the balanced dataset through K-means clustering. (C) The distribution of the predictive dataset through K-means Clustering. (D) ROC curve of all categories and averaged cases.

TABLE 2 LncRNA subcellular localization dataset.

Name	Origin Dataset	Balanced Dataset	Predictive Dataset	
Nucleus	673	673	673	
Cytoplasm	407	407	407	
Cytosol	94	407	394	
Ribosome	152	407	363	
Exosome	16	320	302	

3.2 Predictive modelling process

The predicted results, both before and after sample balancing, were shown in Table 3. From Table 4, we found that exosomes, cytosols, and ribosomes were well-identified in the new predictive dataset. Notably, despite the limited data in the sample, these categories exhibited improved predictive results compared to those that were not oversampled. This result suggested that oversampling may enhance the model’s ability to capture more biological characteristics, thus improving prediction accuracy. This aligns with the concept of biological diversity. Furthermore, the ROC curve was shown in Figure 2D. From Figure 2D, we found that the AUC values for exosome, ribosome, cytoplasm, cytosol, and nucleus are 0.99, 0.94, 0.93, 0.84, and 0.77, respectively. These precise predictive outcomes indicate that the model has a commendable generalization ability across diverse subcellular localization categories.

TABLE 3 The predictive results of each subcellular localization.

Sample Balance	Subcellular localization	S n (%)	S p (%)	MCC	ACC (%)	
Before	Nucleus	98.36	89.28	0.885	86.58	
Cytoplasm	95.33	94.97	0.891	
Cytosol	47.38	98.51	0.604	
Ribosome	32.97	98.87	0.456	
Exosome	6.25	99.92	0.173	
After	Nucleus	58.54	90.22	0.606	94.14	
Cytoplasm	75.81	93.25	0.684	
Cytosol	98.39	99.59	0.965	
Ribosome	99.57	98.34	0.885	
Exosome	99.19	99.79	0.989	

TABLE 4 Predictive results using different feature combinations in 10-fold cross validation.

Feature	S n (%)	S p (%)	ACC (%)	MCC	
3-mer	54.31	70.23	68.35	0.467	
3-RF	62.53	81.23	80.35	0.673	
3-mer+3-RF	84.74	95.05	92.81	0.771	
3-mer+3-RF + MFE	86.38	96.67	94.14	0.829	

3.3 Predictive modelling process

The predicted results, both before and after sample balancing, were shown in Table 3. In Table 3, we found that exosomes, cytosols, and ribosomes were well-identified in the new predictive dataset. Notably, despite the limited data in the sample, these categories exhibited improved predictive results compared to those not oversampled. This result suggested that oversampling may enhance the model’s ability to capture more biological characteristics, thus improving prediction accuracy. This aligns with the concept of biological diversity. Furthermore, ROC curve analysis is conducted, and the results were shown in Figure 2D. In Figure 2D, we found that the AUC values for exosome, ribosome, cytoplasm, cytosol, and nucleus are 0.91, 0.89, 0.99, 0.93, and 0.77, respectively. These precise predictive outcomes indicate that the model has a commendable generalization ability across diverse subcellular localization categories.

To identify which feature categories most influence the predictive model’s results, we conducted separate experiments on 3-mer, 3-RF, MFE, and various combinations of these features. Table 4 shown the predictive results of different feature combinations in the LncSTPred model. We observed that the model’s performance using sequential features closely resembles that of other theoretical models. However, a substantial improvement in efficacy was achieved by adding the MFE feature.

3.4 Comparison with other researchers’ methods

Over the past few decades, numerous lncRNA predictive models have emerged. In this study, we conducted a comparative analysis between LncSTPred and existing theoretical algorithms, as well as other scholars’ predictive models to assess LncSTPred’s performance. Our comparison included Support Vector Machine (SVM), Random Forest, and XGBoost, with results presented in Table 5. MCC is a comprehensive performance metric that effectively reflects the classification ability of the model, especially when dealing with unbalanced datasets. Despite both XGBoost and LncSTPred employing the boosting integration strategy, LncSTPred displayed higher accuracy and generalization ability, particularly when using class-weight parameters for predicting lncRNA subcellular localization. This improvement is attributed to optimized error recognition and adjustments in error recognition point weights within LncSTPred.

TABLE 5 Comparison with existing theoretical algorithm.

Method	Subcellular Localization	S p (%)	S n (%)	MCC	ACC (%)	
SVM	Exosome	66.92	65.53	0.262	73.37	
Cytosol	73.74	64.17	0.306	
Cytoplasm	92.15	48.32	0.445	
Ribosome	88.47	59.23	0.451	
Nucleus	94.11	25.85	0.274	
Random Forest	Exosome	85.91	75.95	0.550	84.77	
Cytosol	85.91	72.85	0.427	
Cytoplasm	95.85	76.70	0.604	
Ribosome	91.69	68.06	0.557	
Nucleus	96.85	48.32	0.553	
Xgboost	Exosome	96.97	94.44	0.906	90.22	
Cytosol	92.31	93.54	0.810	
Cytoplasm	92.78	89.47	0.804	
Ribosome	91.26	78.13	0.680	
Nucleus	93.88	51.35	0.519	
LncSTPred	Exosome	99.79	99.19	0.989	94.14	
Cytosol	99.58	98.39	0.965	
Cytoplasm	93.25	75.81	0.684	
Ribosome	98.34	99.57	0.885	
Nucleus	90.22	58.54	0.606	

We compared LncSTPred with machine learning-based models, specifically iLoc-lncRNA 2.0 and LightGBM-LNCLOC, using datasets from these models for performance validation (Table 6). Both iLoc-lncRNA 2.0 and LncSTPred demonstrated high precision in ribosome and exosome categories, with superior sensitivity in predicted cytoplasmic and nuclear localizations. In contrast, lncLocation, which did not use oversampling, showed comparable accuracy in nuclear and cytoplasmic categories but lower performance in ribosomal and exosomal predictions. This indicated that oversampling enhanced lncRNA subcellular localization prediction, particularly improved sensitivity in the cytoplasmic category compared to lncLocator. The analysis of classification datasets highlighted LncSTPred’s robust adaptability across various datasets, affirming its reliability as a classification model.

TABLE 6 Comparison with previous state-of-the-art methods.

Method	Subcellular Localization	S p (%)	S n (%)	MCC	ACC (%)	
iLoc-lncRNA
2.0	Nucleus	95.59	91.03	86.59	91.60	
Cytoplasm	98.96	94.37	94.59	
Ribosome	99.01	83.72	85.71	
Exosome	99.36	66.67	83.33	
lncLocation	Nucleus	—	74.19	95.83	87.78	
Cytoplasm	—	100	85.00	
Ribosome	—	55.56	100	
Exosome	—	33.33	100	
LncSTPred in Lin’s dataset	Nucleus	78.74	90.06	60.48	90.76	
Cytoplasm	95.20	75.12	67.06	
Ribosome	99.84	97.67	97.53	
Exosome	99.84	96.67	92.90	
LightGBM-LNCLOC	Nucleus	—	90.00	42.90	70.60	
Cytoplasm	—	40.00	66.67	
Cytosol	—	50.00	100	
Ribosome	—	40.00	40.00	
Exosome	—	42.90	100	
LncSTPred in Li’s dataset	Nucleus	81.25	81.29	61.39	87.68	
Cytoplasm	87.82	60.49	50.92	
Cytosol	93.58	53.62	45.31	
Ribosome	94.51	47.69	43.03	
Exosome	99.83	6.25	75.43	

3.5 Analysis of feature importance

We computed the feature importance in LncSTPred using the F-score described in Section 2.4, and the results were depicted in Figure 3. The top 10 combinations were “ACA” and “AUA” in the third RF, “AUA”, “AGA”, and “UCA” in the second RF, and “GGU” in the first RF. Additionally, “AUA”, “UAU”, “CGG”, and “GGC” were prominent in the 3-mer. The triplex nucleotide “ANA” had a significant impact on predictive modeling. To explore the potential effects of these 10 triplexes on the subcellular localization of lncRNAs, we analyzed their distribution in lncRNA sequences and secondary structure substructures.

FIGURE 3 Top 10 features of the LncSTPred model.

Figure 4A displayed the frequency of triplex nucleotides in the original dataset. Interestingly, ACA and AGA frequencies exceeded the random probability, whereas AUA and UAU frequencies were below the random probability in the lncRNA sequences.

FIGURE 4 Biological analysis of the 10 most important features: (A) Frequency distribution of triplex nucleotides in lncRNA sequences. The X-axis represents 64 types of triplexes nucleotides, and the Y-axis represents the frequency of triplexes. (B) Three types of substructures in the secondary structure of lncRNAs.

We used the RNAfold software to predict the secondary structure of lncRNA (Gruber et al., 2008). The output was in dot-bracket format, where '...' denoted loop structures, ')))' or '(((' represented stem structures, and '..(', '..)',').. ', or '(..' indicated junction structures. Figure 4B displayed three distinct substructures observed in the secondary structure. A frequency analysis of these substructure types was conducted, and the results were presented in Table 7. Analysis of Table 7 revealed specific enrichments: ACA in the junction, and AGA and TAT in the stem structures. While the frequency of ATA in the predictive dataset was lower than that of random combinations, ATA, ACA, and AGA predominantly appeared at critical positions within the junction and stem substructures of the lncRNA secondary structure. These findings suggested the significant contributory role of ANA in the construction of RNA secondary structures.

TABLE 7 Frequency distribution of triplex nucleotide in secondary structure components.

	AUA (%)	ACA (%)	AGA (%)	UAU (%)	Random combination (%)	
Junction	1.632	3.081	2.998	1.382	1.563	
Loop	1.414	2.691	1.938	1.415	1.563	
Stem	1.613	2.232	3.347	2.101	1.563	

4 Discussion

In recent years, the recognition of lncRNA subcellular localization has garnered increasing attention, as researchers have realized its potential for discovering the function of LncRNA. In this paper, we propose an improved algorithm for predicting lncRNA subcellular localizations, called LncSTPred.

During the establishment of the predictive model, we recognized the significant impact of the predictive dataset on the results. To address this, we utilized three oversampling methods—Borderline-SMOTE, ADASYN, and UNCERTAIN WEIGHT—to oversample the sample sets. This ensured that the predictive model for each category of samples could be sufficiently trained to achieve the best predictive results. Constructing LncSTPred using the improved Adaboost algorithm, we achieved 94.14% accuracy in the 5-categorical dataset and 90.76% accuracy in Lin’s 4-categorical dataset. This demonstrates that our improvements to the Adaboost algorithm, combined with data balancing, can provide better results than using the class-weight parameter in other algorithms.

Several studies have explored the impact of specific nucleotide combinations on RNA secondary structure. For example, the predictions of Smith et al. accurately identify mascRNA and a conserved hairpin upstream of Evolutionarily Conserved Structures (ECS). They observed that “ANA” triplex nucleotides predominantly appear at the stem-loop junction in ECS (Smith et al., 2013). Novikova et al., through biochemical probing, delineated a complex, two-dimensional structure comprising distinct sub-domains, including helical segments, terminal loops, internal loops, and linker regions. This study underscores that purine-rich sequences are highly conserved and often situated in single-stranded regions such as terminal and internal loops (Novikova et al., 2012). These findings corroborate the involvement of “ANA” triplex nucleotide composition in lncRNA secondary structure. Additionally, we performed a quantitative analysis of feature importance to identify the most significant features. By analyzing the frequency of triplex nucleotides and the stem-loop structures of lncRNA, we aimed to understand the relationship between significant features and lncRNA subcellular localization. Our analysis revealed a bias in the frequency of ANA nucleotide combinations within triplex nucleotides and the substructural frequency of stem-loop structures. These findings suggest that ANA nucleotide combinations play key roles in the composition of lncRNA secondary structures. In previous studies, Constanty et al. found that conserved U-rich and A-rich motifs were associated with specific processing and localization functions of lncRNAs like NEAT1 and MALAT1 (Constanty and Shkumatava, 2021). Furthermore, Cai et al. provided evidence that specific triplexes, including ACA, ATA, and AGA, significantly influence localization patterns by analyzing various sequences (Cai et al., 2023). Lyu et al. emphasized the relevance of trinucleotide propensity and position-specific features in recognizing lncRNA subcellular localization, demonstrating that specific triplexes like UAU could play a role in these predictions (Lyu et al., 2023).

Although LncSTPred has achieved better results in predicting lncRNA subcellular localization, we still face some challenges. On the one hand, the training process of AdaBoost is more complex, leading to significantly higher computing time compared to other prediction models. On the other hand, LncSTPred currently only accurately predicts lncRNAs only localized to a single subcellular localization, whereas many lncRNAs are localized in multiple subcellular localizations. Therefore, we will focus on lncRNA subcellular localization prediction in the future, further enhancing the accuracy of predicting subcellular localization of lncRNAs and the prediction of multi-localized lncRNAs using deep learning algorithms.

Data availability statement

The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

Author contributions

S-LH: Data curation, Methodology, Software, Writing–original draft, Writing–review and editing. Y-LC: Funding acquisition, Writing–review and editing. L-QZ: Writing–review and editing. HB: Writing–review and editing. J-HY: Writing–review and editing. Q-ZL: Conceptualization, Funding acquisition, Writing–review and editing.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2024.1452142/full#supplementary-material
==== Refs
References

Atianand M. K. Caffrey D. R. Fitzgerald K. A. (2017). Immunobiology of long noncoding RNAs. Annu. Rev. Immunol. 35 , 177–198. 10.1146/annurev-immunol-041015-055459 28125358
Batista P. J. Chang H. Y. (2013). Long noncoding RNAs: cellular address codes in development and disease. Cell 152 (6 ), 1298–1307.23498938
Birney E. Stamatoyannopoulos J. A. Dutta A. Guigó R. Gingeras T. R. Margulies E. (2007). Identification and analysis of functional elements in 1% of the human genome by the ENCODE pilot project. Nature 447 (7146 ), 799–816. 10.1038/nature05874 17571346
Bridges M. C. Daulagala A. C. Kourtidis A. (2021). LNCcation: lncRNA localization and function. J. Cell Biol. 220 , 2020090455–e202009117. 10.1083/jcb.202009045
Cai J. Wang T. Deng X. Tang L. Liu L. (2023). GM-lncLoc: LncRNAs subcellular localization prediction based on graph neural network with meta-learning. BMC Genomics 24 (52 ), 52–14. 10.1186/s12864-022-09034-1 36709266
Clemson C. M. Hutchinson J. N. Sara S. A. Ensminger A. W. Fox A. H. Chess A. (2009). An architectural role for a nuclear noncoding RNA: NEAT1 RNA is essential for the structure of paraspeckles. Mol. Cell. 33 , 717–726. 10.1016/j.molcel.2009.01.026 19217333
Constanty F. Shkumatava A. (2021). lncRNAs in development and differentiation: from sequence motifs to functional characterization. Develop 148 (1 ), dev182741–11. 10.1242/dev.182741
Cui T. Dou Y. Tan P. Ni Z. Liu T. Wang D. (2022). RNALocate v2.0: an updated resource for RNA subcellular localization with increased coverage and annotation. Nucleic Acids Res. 50 , D333–D339. 10.1093/nar/gkab825 34551440
Derrien T. Johnson R. Bussotti G. Tanzer A. Djebali S. Tilgner H. (2012). The GENCODE v7 catalog of human long noncoding RNAs: analysis of their gene structure, evolution, and expression. Genome Res. 22 , 1775–1789. 10.1101/gr.132159.111 22955988
Esguerra J. L. S. Eliasson L. (2014). Functional implications of long non-coding RNAs in the pancreatic islets of Langerhans. Front. Genet. 5 , 209–9. 10.3389/fgene.2014.00209 25071836
Feng S. Liang Y. Du W. Lv W. Li Y. (2020). LncLocation: efficient subcellular location prediction of long non-coding RNA-based multi-source heterogeneous feature fusion. Int. J. Mol. Sci. 21 , 7271. 10.3390/ijms21197271 33019721
Flynn R. A. Chang H. Y. (2014). Long noncoding RNAs in cell-fate programming and reprogramming. Cell Stem Cell 14 , 752–761. 10.1016/j.stem.2014.05.014 24905165
Gong C. Maquat L. E. (2011). lncRNAs transactivate STAU1-mediated mRNA decay by duplexing with 3' UTRs via Alu elements. Nature 470 , 284–288. 10.1038/nature09701 21307942
Gonzalez I. Munita R. Agirre E. Dittmer T. A. Gysling K. Misteli T. (2015). A lncRNA regulates alternative splicing via establishment of a splicing-specific chromatin signature. Nat. Struct. Mol. Biol. 22 , 370–376. 10.1038/nsmb.3005 25849144
Gruber A. R. Lorenz R. Bernhart S. H. Neuböck R. Hofacker I. L. (2008). The Vienna RNA websuite. Nucleic Acids Res. 36 , 70–74. 10.1093/nar/gkn188
Han H. Wang W. Y. Mao B. H. (2005). Borderline-SMOTE: a new over-sampling method in imbalanced data sets learning. Int. Conf. Int. Comp., 878–887. 10.1007/11538059_91
Harries L. W. (2012). Long non-coding RNAs and human disease. Biochem. Soc. Trans. 40 , 902–906. 10.1042/BST20120020 22817756
He H. Bai Y. Garcia E. A. Li S. (2008). “ADASYN: adaptive synthetic sampling approach for imbalanced learning,” in Ieee int. Jt. Conf. Neural netw. (WCCI) (Hong Kong), 1322–1328.
Huang Y. Niu B. Gao Y. Fu L. Li W. (2010). CD-HIT Suite: a web server for clustering and comparing biological sequences. Bioinformatics 26 (5 ), 680–682. 10.1093/bioinformatics/btq003 20053844
Jiang W. Liu Y. Liu R. Zhang K. Zhang Y. (2015). The lncRNA DEANR1 facilitates human endoderm differentiation by activating FOXA2 expression. Cell Rep. 11 , 137–148. 10.1016/j.celrep.2015.03.008 25843708
Kameswaran V. Kaestner K. H. (2014). The missing lnc(RNA) between the pancreatic β-cell and diabetes. Front. Genet. 5 , 200–210. 10.3389/fgene.2014.00200 25071830
Kendall A. Gal Y. Cipolla R. (2018). Multi-task learning using uncertainty to weigh losses for scene geometry and semantics. Proc. IEEE Conf. Comp. Vis. Pat. Reco, 7482–7491.
Kitagawa M. Kotake Y. Ohhata T. (2012). Long non-coding RNAs involved in cancer development and cell fate determination. Curr. Drug Targets 13 , 1616–1621. 10.2174/138945012803530026 22974399
Kugel J. F. Goodrich J. A. (2012). Non-coding RNAs: key regulators of mammalian transcription. Trends Biochem. Sci. 37 , 144–151. 10.1016/j.tibs.2011.12.003 22300815
Lauressergues D. Couzigou J. M. Clemente H. S. Martinez Y. Dunand C. Bécard G. (2015). Primary transcripts of microRNAs encode regulatory peptides. Nature 520 , 90–93. 10.1038/nature14346 25807486
Li W. Godzik A. (2006). Cd-hit: a fast program for clustering and comparing large sets of protein or nucleotide sequences. Bioinformatics 22 (13 ), 1658–1659. 10.1093/bioinformatics/btl158 16731699
Li Z. Li X. Wu S. Xue M. Chen W. (2014). Long non-coding RNA UCA1 promotes glycolysis by upregulating hexokinase 2 through the mTOR-STAT3/microRNA143 pathway. Cancer Sci. 105 , 951–955. 10.1111/cas.12461 24890811
Luo M. Li Z. Wang W. Zeng Y. Liu Z. Qiu J. (2013). Long non-coding RNA H19 increases bladder cancer metastasis by associating with EZH2 and inhibiting E-cadherin expression. Cancer Lett. 333 , 213–221. 10.1016/j.canlet.2013.01.033 23354591
Lyu J. Zheng P. Qi Y. Huang G. (2023). LightGBM-LncLoc: a LightGBM-based computational predictor for recognizing long non-coding RNA subcellular localization. Mathematics 11 (3 ), 602–614. 10.3390/math11030602
Martens J. A. Laprade L. Winston F. (2004). Intergenic transcription is required to repress the Saccharomyces cerevisiae SER3 gene. Nature 429 , 571–574. 10.1038/nature02538 15175754
Melé M. Rinn J. L. (2016). “Cat’s Cradling” the 3D genome by the act of LncRNA transcription. Mol. Cell 62 , 657–664. 10.1016/j.molcel.2016.05.011 27259198
Miao H. Wang L. Zhan H. Dai J. Chang Y. Wu F. (2019). A long noncoding RNA distributed in both nucleus and cytoplasm operates in the PYCARD-regulated apoptosis by coordinating the epigenetic and translational regulation. PLoS Genet. 15 , 10081444–e1008224. 10.1371/journal.pgen.1008144
Mourtada-Maarabouni M. Pickard M. R. Hedge V. L. Farzaneh F. Williams G. T. (2009). GAS5, a non-protein-coding RNA, controls apoptosis and is downregulated in breast cancer. Oncogene 28 , 195–208. 10.1038/onc.2008.373 18836484
Novikova I. V. Hennelly S. P. Sanbonmatsu K. Y. (2012). Sizing up long non-coding RNAs: do lncRNAs have secondary and tertiary structure? Bioarchitecture 2 (6 ), 189–199. 10.4161/bioa.22592 23267412
Paraskevopoulou M. D. Georgakilas G. Kostoulas N. Reczko M. Maragkakis M. Dalamagas T. M. (2013). DIANA-LncBase: experimentally verified and computationally predicted microRNA targets on long non-coding RNAs. Nucleic Acids Res. 41 , 239–245. 10.1093/nar/gks1246
Rainey S. Repka J. (2013). Quantitative sequence and open reading frame analysis based on codon bias. J. Syst. Cyberne Inf. 4 (1 ), 65–72.
Saxena A. Carninci P. (2011). Long non-coding RNA modifies chromatin: epigenetic silencing by long non-coding RNAs. BioEssays News Rev. Mol. Cell Dev. Biol. 33 , 830–839. 10.1002/bies.201100084
Schapire R. E. Singer Y. (1999). Improved boosting algorithms using confidence-rated predictions. Mach. Learn 37 , 297–336. 10.1023/a:1007614523901
Smith M. A. Gesell T. Stadler P. F. Mattick J. S. (2013). Widespread purifying selection on RNA structure in mammals. Nucleic Acids Res. 41 (17 ), 8220–8236. 10.1093/nar/gkt596 23847102
Sparano J. A. Gray R. J. Makower D. F. Pritchard K. I. Albain K. S. Hayes D. F. (2018). Adjuvant chemotherapy guided by a 21-gene expression assay in breast cancer. N. Engl. J. Med. 379 , 111–121. 10.1056/NEJMoa1804710 29860917
Su Z. D. Huang Y. Zhang Z. Y. Zhao Y. W. Wang D. Chen W. (2018). iLoc-lncRNA: predict the subcellular location of lncRNAs by incorporating octamer composition into general PseKNC. Bioinforma. Oxf Engl. 34 , 4196–4204. 10.1093/bioinformatics/bty508
Wang J. Liu X. Wu H. Ni P. Gu Z. Qiao Y. (2010). CREB up-regulates long non-coding RNA, HULC expression through interaction with microRNA-372 in liver cancer. Nucleic Acids Res. 38 , 5366–5383. 10.1093/nar/gkq285 20423907
Wang K. Liu C. Y. Zhou L. Y. Wang J. X. Wang M. Zhao B. (2015). APF lncRNA regulates autophagy and myocardial infarction by targeting miR-188-3p. Nat. Commun. 6 (1 ), 6779. 10.1038/ncomms7779 25858075
Wang K. C. Chang H. Y. (2011). Molecular mechanisms of long noncoding RNAs. Mol. Cell. 43 , 904–914. 10.1016/j.molcel.2011.08.018 21925379
Wang Z. Li X. (2013). The role of noncoding RNA in hepatocellular carcinoma. Gland. Surg. 2 , 25–29. 10.3978/j.issn.2227-684X.2013.02.07 25083452
Wilusz J. E. (2016). Long noncoding RNAs: re-writing dogmas of RNA processing and stability. Biochim. Biophys. Acta 1859 , 128–138. 10.1016/j.bbagrm.2015.06.003 26073320
Wu W. Bhagat T. D. Yang X. Song J. H. Cheng Y. Agarwal R. (2013). Hypomethylation of noncoding DNA regions and overexpression of the long noncoding RNA, AFAP1-AS1, in Barrett’s esophagus and esophageal adenocarcinoma. Gastroenterology 144 , 956–966. 10.1053/j.gastro.2013.01.019 23333711
Xiao T. Liu L. Li H. Sun Y. Luo H. Li T. (2015). Long noncoding RNA ADINR Regulates adipogenesis by transcriptionally activating C/EBPα. Stem Cell Rep. 5 , 856–865. 10.1016/j.stemcr.2015.09.007
Xie J. Y. Wang C. Jiang S. Zhang Y. (2010). Feature selection method combing improved F-score and support vector machine. J. Comput. Appl. 30 , 993–996.
Yan L. Zhou J. Gao Y. Ghazal S. Lu L. Bellone S. (2015). Regulation of tumor cell migration and invasion by the H19/let-7 axis is antagonized by metformin-induced DNA methylation. Oncogene 34 , 3076–3084. 10.1038/onc.2014.236 25088204
Yang F. Deng X. Ma W. Berletch J. B. Rabaia N. Wei G. (2015). The lncRNA Firre anchors the inactive X chromosome to the nucleolus by binding CTCF and maintains H3K27me3 methylation. Genome Biol. 16 , 52–17. 10.1186/s13059-015-0618-0 25887447
Zeng M. Wu Y. Lu C. Zhang F. Wu F. X. Li M. (2022). DeepLncLoc: a deep learning framework for long non-coding RNA subcellular localization prediction based on subsequence embedding. Brief. Bioinform 23 , bbab360–10. 10.1093/bib/bbab360 34498677
Zhang Z. Y. Sun Z. J. Yang Y. H. Lin H. (2022). Towards a better prediction of subcellular location of long non-coding RNA. Front. Comput. Sci. 16 , 165903–165917. 10.1007/s11704-021-1015-3
Zhu M. Chen Q. Liu X. Sun Q. Zhao X. Deng R. (2014a). LncRNA H19/miR-675 axis represses prostate cancer metastasis by targeting TGFBI. FEBS J. 281 , 3766–3775. 10.1111/febs.12902 24988946
Zhu Y. P. Bian X. J. Ye D. W. Yao X. D. Zhang S. L. Dai B. (2014b). Long noncoding RNA expression signatures of bladder cancer revealed by microarray. Oncol. Lett. 7 , 1197–1202. 10.3892/ol.2014.1843 24944692
Zuker M. Sankoff D. (1984). RNA secondary structures and their prediction. Bull. Math. Biol. 46 , 591–621. 10.1016/s0092-8240(84)80062-2
Zuker M. Stiegler P. (1981). Optimal computer folding of large RNA sequences using thermodynamics and auxiliary information. Nucleic Acids Res. 9 , 133–148. 10.1093/nar/9.1.133 6163133
