Skip to main content

ScLSTM: single-cell type detection by siamese recurrent network and hierarchical clustering



Categorizing cells into distinct types can shed light on biological tissue functions and interactions, and uncover specific mechanisms under pathological conditions. Since gene expression throughout a population of cells is averaged out by conventional sequencing techniques, it is challenging to distinguish between different cell types. The accumulation of single-cell RNA sequencing (scRNA-seq) data provides the foundation for a more precise classification of cell types. It is crucial building a high-accuracy clustering approach to categorize cell types since the imbalance of cell types and differences in the distribution of scRNA-seq data affect single-cell clustering and visualization outcomes.


To achieve single-cell type detection, we propose a meta-learning-based single-cell clustering model called ScLSTM. Specifically, ScLSTM transforms the single-cell type detection problem into a hierarchical classification problem based on feature extraction by the siamese long-short term memory (LSTM) network. The similarity matrix derived from the improved sigmoid kernel is mapped to the siamese LSTM feature space to analyze the differences between cells. ScLSTM demonstrated superior classification performance on 8 scRNA-seq data sets of different platforms, species, and tissues. Further quantitative analysis and visualization of the human breast cancer data set validated the superiority and capability of ScLSTM in recognizing cell types.

Peer Review reports


The human body comprises roughly 40 trillion cells, each exhibiting a remarkable diversity of shapes and functions [1]. Identifying and visualizing cell types offers valuable insight into cellular heterogeneity and reveals specific mechanisms underlying pathological conditions. Remarkable advancements in scRNA-seq technology have made it possible to cost-effectively and efficiently study genome-wide expression at the single-cell level. This allows us to determine cell types by examining the transcriptome status of thousands of individual cells.

Traditional sequencing methods get overall heterogeneity at the transcriptome and phenotypic levels by detecting the total signal of a cell population and calculating the average, which fails to illustrate the individual states and differences among cells [2]. The objective of scRNA-seq technology is to sequence the genome at the single-cell level, offering a distinct advantage in uncovering the subtle variations unique to individual cells and comprehending the diversity in gene expression within biological tissues. Identifying single-cell types serves as the foundation for exploring cellular heterogeneity and developmental processes.

With the ongoing innovation and improvement of scRNA-seq technology, the throughput of single-cell technology has substantially improved, enabling the detection of both common and rare cell types. scRNA-seq technologies consist not only of low-throughput sequencing technologies like Smart-Seq and Celseq but also of high-throughput sequencing technologies like inDrop, 10X-v2, and 10X-v3. The gene expression data generated by various types of scRNA-seq technologies constitute a vast and extensive database of scRNA-seq. Gene Expression Omnibus (GEO), ArrayExpress, PanglaoDB, and Human Cell Atlas (HCA), among others, are frequently used single-cell data repositories. The data used in this study were collected from the GEO database and ArrayExpress.

Among the cell type detection models based on scRNA-seq data, the models based on cell-cell similarity are the most classic. For example, Corr evaluates cell-cell associations from a global analysis of variance based on differential gene expression of cells [3]. POCR selects different kernel embedding techniques for single-cell data sets of varying sizes; Applying a Gaussian kernel for small-scale data sets and a linear kernel for large-scale data sets; On the resulting kernel similarity, spectral clustering analysis is then performed [4]. Based on the assumption that the same cell type is in the same subspace, SinNLRR describes the expression of the same cell type through a cell expression; SinNLRR finds the low-rank and non-negative representation of the expression matrix from all candidate subspaces and optimizes it using the alternating direction multiplier method (ADMM) [5]. SIMLR produces similar symmetric matrices by learning kernel functions with varying weights, then decomposes the similar matrices into approximation diagonal matrices based on the number of cell types [6]. SCENA sorts the gene expression levels in descending order of variance, selects genes to construct cell-cell similarity matrices, and combines the spectral clustering results of all matrices to obtain the final clustering results [7]. However, the above methods cannot reflect the complex nonlinear relationship between genes, affecting the cell type detection results.

Due to its robust feature learning capabilities, deep learning has proven to be advantageous when dealing with large-scale scRNA-seq data with highly complicated features. Previous research has demonstrated that neural networks can extract insights from single-cell gene expression data. The scDCC model proposed by Tian et al. encodes prior single-cell knowledge as constraint information and integrates it into the clustering procedure via a new loss function [8]. scAdapt combines domain adaptation with semi-supervised learning based on virtual adversarial training for accurate cell classification using labeled sources and unlabeled target data [9]. Song et al. proposed to use a graph convolutional network for cell type detection by constructing a mixture of inter and intra-data set cell mapping graphs [10]. The success of deep learning depends on massive datasets and robust computational resources. Despite the widespread application of deep learning frameworks such as graph neural networks and capsule networks to various problems, they continue to face the challenge of limited sample sizes, particularly in the classification of single-cell types when the amount of data varies significantly between cell types [11]. Meta-learning has been proven effective in numerous fields, such as natural language processing and robotics, due to its crucial framework for addressing the sample imbalance problem in deep learning.

Due to the varied experimental sequencing platforms in different laboratories, the distribution features of the single-cell data collected from different laboratories will vary. As indicated in Table 1, the number of samples (cells) in single-cell data sets ranges from tens to hundreds. In addition, the number of cells of various types in a single-cell data set can vary by a factor of ten. The difference in feature distribution of single-cell data sets and the imbalance of the number of cells between cell types make it challenging to identify cell types with fewer cells.

Siamese network can calculate similarity based on the distance between two samples and simultaneously enhance the learning model with the knowledge gathered from several learning events. LSTM is a unique type of recurrent neural network (RNN) that can tackle the problems of vanishing and exploding gradients during lengthy sequence training. We are interested in integrating siamese and LSTM networks to discover relevant feature representations with the potential to enhance cell type detection.

In this work, we propose a single-cell type detection model called ScLSTM, which integrates the siamese and LSTM networks. The network architecture of ScLSTM is shown in Fig. 1. ScLSTM employs an improved sigmoid kernel to calculate a measure of cell similarity. The “siamese” of a siamese LSTM is achieved by sharing weights between two identical LSTMs. Siamese LSTM takes two inputs and maps them to a new space using two LSTMs. ScLSTM learns how to minimize the distance between single-cell data of the same category and maximize the distance between different categories, enabling us to obtain more discriminative features for each cell. Subsequently, the agglomerative clustering algorithm was used to cluster single cells. The superiority of ScLSTM, compared to several state-of-the-art approaches, is validated on 8 different human and mouse scRNA-seq datasets. When comparing the three evaluation indicators, ARI, NMI, ACC, and BAS, the performance of ScLSTM in clustering is superior to that of other methods. Additionally, we apply ScLSTM to the human breast cancer dataset and two 10x Genomics datasets, demonstrating that ScLSTM can effectively categorize cancer cell subtypes. We demonstrate that the semantic representation of cells (semantic correlation between cell pairs) generated by siamese LSTM networks is more suitable for representing cells in the cell type detection problem. ScLSTM, an enhanced version of the meta-learning framework, is more effective at balancing datasets and identifying a small proportion of cell subtypes, while also performing exceptionally well on large-scale data sets.

Fig. 1
figure 1

Network architecture of ScLSTM. The improved sigmoid kernel is used to build the scRNA-seq feature matrix, which is then fed into the siamese LSTM. Using the agglomerative clustering algorithm, determine the cell type of the output of the siamese LSTM


In single-cell-related research, the relationship between different cells can be assessed based on cell similarity calculated from gene expression. Since the scale and distribution of different single-cell transcriptome data sets are different, designing a similarity calculation method for different types of single-cell data sets is one of the significant difficulties. The Euclidean distance, Pearson, and Spearman similarity computation methods are all commonly utilized. The Euclidean distance represents the straight-line distance between two points, the Pearson correlation coefficient is used to test the linear correlation between the data, and the Spearman tests the rank correlation coefficient. We employ an improved sigmoid kernel embedding approach to express cell similarity to fully extract the global information of the single-cell data set structure and the intra-class information of the cell population.

The input of ScLSTM is the expression matrix \({\textbf{A}}\), where rows correspond to cells and columns correspond to genes. Suppose there are n cells in \({\textbf{A}}\), the gene expression of cell i is denoted by \(A_i=(a_{i1},a_{i2},\dots ,a_{im})\), \(i\in (1,2,\dots , n)\), where m represents the number of genes. Next, a logarithmic transformation is executed on each element in matrix \({\textbf{A}}\) to obtain matrix \({\textbf{L}}\), where \({\textbf{L}}=log({\textbf{A}}+1)\). For the matrix \({\textbf{L}}\), the inner product matrix \({\textbf{I}}\) is obtained by multiplying the row vectors, \(I_{ij}=L_i*L_j\), where \(L_i\) and \(L_j\) are both row vectors of the matrix \({\textbf{L}}\). Due to the different sizes and distributions of different single-cell data sets, the inner product matrix is updated, \(I_{ij} = (I_{ij}-I_{min})/(I_{max}-I_{min})\), where \(I_{max}\) is the largest element in matrix \({\textbf{I}}\), and \(I_{min}\) is the the smallest element in matrix \({\textbf{I}}\). Finally, the sigmoid kernel similarity matrix \({\textbf{S}}\) is calculated through the updated inner product matrix, \(S_{ij}=(e^{I_{ij}}-e^{-I_{ij}})/(e^{I_{ij}}+e^{-I_{ij}})\).

Siamese recurrent network architecture of the ScLSTM model

The siamese network two identical neural networks, each encoding different features of individual cells and mapping these features to a new embedding space for comparison. A siamese LSTM network is constructed by combining two identical LSTM networks with the same structure and characteristics. The LSTM structure is illustrated in Fig. 2. LSTM have two transmission states, \(c_t\) (cell state) and \(h_t\) (hidden state). \(h_t\) can be seen as short- term memory for the current knowledge, which is transformed by the tanh function. The purpose of \(c_t\) is to transform and process the memory of the last and current times using a linear transformation. LSTM is composed of three gates:

  • Forget stage. The retention of past information is determined by assessing the importance of current input information. The first step in LSTM is to decide what information to discard from the cell state. This decision is made through the forget gate layer. This gate reads \(h_t-1\) and \(A_t\), and outputs a number between 0 and 1 for each element in the cell state \(c_t-1\). The value of 1 means completely retained, while 0 means completely discarded.

  • Input gate. The retention degree of input information is determined by assessing the importance of the current input information. The decision of how much new information to add to the cell state is made. Firstly, the sigmoid layer of the input gate determines which information needs to be updated. Then, a tanh layer generates a vector, denoted as \(c_t\), which is an alternative for updating. These two parts are then combined to produce an update to the cell state.

  • Output gate. How much the current output depends on the current memory cell. This output is based on the cell state, which is also a filtered version. First, a sigmoid layer is used to determine which part of the cell state will be output. Next, the cell state is processed by tanh (producing a value between −1 and 1), and this result is multiplied by the output of the sigmoid gate, which selects the part that determines to the output.

After performing kernel embedding, we obtain the matrix S. Through a base model consisting of two symmetric LSTMs, the matrix S is mapped to the head model. The structure of the head model is shown in Fig. 1, which consists of a concatenate layer, two dense layers, and two sigmoid layers. The output confusion matrix of the siamese LSTM serves as the input for the agglomerative clustering.

Agglomerative clustering

Agglomerative clustering is a commonly used hierarchical clustering algorithm that can partition data sets into multiple levels to create a tree-like clustering structure. Initially, agglomerative clustering classifies each object as a cluster and subsequently merges these clusters based on a distance measure function until the desired number of clusters is achieved. In the case of ScLSTM, we apply the similarity matrix learned by the siamese LSTM to the agglomerative clustering process. Figure 1 depicts the agglomerative clustering process utilized by our method.

Fig. 2
figure 2

The chain structure of LSTM; Each box represents a neural network layer, consisting of weights, biases, and activation functions; Each circle represents an element-level operation; The arrow indicates the vector flow direction. A forked arrow indicates a copy of the vector


Data sets and performance assessment

Human and mouse cells yielded 8 publically available scRNA-seq data sets (see Table 1 for a complete list) that were used to compare the performance of ScLSTM with other comparative techniques. Five scRNA-seq data sets were obtained from the NCBI GEO database (access IDs: GSE36552, GSE81252, GSE57249, GSE59739, and GSE81861), and one scRNA-seq data set was downloaded from the NCBI SRA database (access ID: SRP041736). Two scRNA-seq data sets were downloaded from the ArrayExpress database (access IDs: E-MTAB-3321 and E-MTAB-2600). The scRNA-seq data sets vary in size from tens to tens of thousands of cells.

Notably, the data sets we collected include those with severely unbalanced cell type,s such as GSE36552, and data sets with extremely limited data, such as GSE57249 and GSE36552. In  the GSE36552 data set, the sample size of the cell type with the largest proportion is 10 times that of the smallest proportion. There are 90 cells of 7 cell types in the GSE36552 data set and 56 cells of 4 cell types in the GSE57249 data set. Additionally, we collected a gold standard labeled data set, the Chung data set [12] (GSE75688), to further validate the robustness and generalization capability of ScLSTM. We also collected two 10x Genomics data sets to assess the performance of ScLSTM on large-scale data. The sc_10x_5cl data set (access ID: GSE118767) five cell lines (HCC827, H1975, H2228, H838, and A549) using the 10x Chromium Genomics protocol, totaling 3919 cells. The second data set is the Jurkat data set, which consist of two cell lines (293T cells and jurkat cells), using UMIs and the droplet-based protocol from 10x Genomics, with a total of 6143 cells (, We applied the same preprocessing step for both 10x Genomics datasets, filtering out genes expressed in fewer  than three cells.

Table 1 The description of data sets used in experiments

We use three metric indicators to measure the performance of ScLSTM: normalized mutual information (NMI), accuracy (ACC), adjusted rand index (ARI), and balaced accuracy score (BAS). In particular, the ranges of NMI and ACC are between 0 and 1, whereas ARI can be negative. A larger value indicates a more remarkable agreement between predicted and actual labels. NMI, ACC, ARI, and BAS are calculated as follows:

$$\begin{aligned} \begin{aligned} NMI(T,P)&= \frac{2I(T,P)}{H(T)+H(P)}, \end{aligned} \end{aligned}$$
$$\begin{aligned} \begin{aligned} ACC&=\frac{\sum _{i=1}^{n} \delta \left( s_{i}, {\text {Map}(\text {r}_\text {i})}\right) }{n}, \end{aligned} \end{aligned}$$
$$\begin{aligned} \begin{aligned} \delta (x, y)&=\left\{ \begin{array}{ll}1, &{} x=y, \\ 0, &{} \text { otherwise, }\end{array}\right. \end{aligned} \end{aligned}$$
$$\begin{aligned} \begin{aligned} ARI(T, P)&=\frac{\sum _{ij}{n_{ij}\atopwithdelims ()2}-[\sum _{i}{a_i\atopwithdelims ()2}\sum _{j}{b_j\atopwithdelims ()2}]/{n \atopwithdelims ()2}}{\frac{1}{2}[\sum _i{a_i\atopwithdelims ()2}+\sum _j{b_i\atopwithdelims ()2}]-[\sum _i{a_i\atopwithdelims ()2}\sum _j{b_j\atopwithdelims ()2}]/{n \atopwithdelims ()2}}, \end{aligned} \end{aligned}$$
$$\begin{aligned} \begin{aligned} BAS&= \frac{(TPR+TNR)}{2} \end{aligned} \end{aligned}$$

where T represents the actual cell cluster and P represents the anticipated cell clusters. In equation (1), H denotes entropy and I(TP) represents the mutual information between T and P. In equation (2), \(r_{i}\) represents the predicted label, \(s_{i}\) represents the true label, n represents the total number of samples, \(Map(\cdot )\) represents the mapping function that maps the predicted labels to the equivalent true label and can be obtained by applying the Hungarian Algorithm [21]; equation (3) is the indicator function. In equation (5), TPR represents the probability that recall correctly predicts the positive class, while TNR represents the probability that the prediction for the negative class is correct.

Influence of network structure on ScLSTM

Comparative analysis of clustering

To evaluate ScLSTM, we applied it to scRNA-seq data sets generated by different sequencing platforms (see Table 1), all of which contained real labels. We choose six state-of-the-art models (Corr [3], POCR [4], SIMLR [6], SinNNLR [5], SNN-cliq [22], and ZIFA [23]) for comparison. Corr evaluates the relationship between cells based on an analysis of variance. POCR and SIMLR use kernel embedding to assess the similarity between cells. Based on similarity learning, SinNLRR imposes a non-negative low-rank structure on the similarity matrix. SNN-Cliq is a clustering algorithm based on a new shared nearest neighbor graph and quasi-clique finding techniques. Each experimental result was obtained through a grid search strategy. Figure 3 displays heatmaps of 7 models across 8 data sets using 4 evaluation metrics. Darker colors indicate closer performance to 1. The results of ScLSTM were obtained using the true number of clusters. It is evident that ScLSTM outperforms all other models in these four metrics: ACC, ARI, BAS, and NMI, followed by SinNNLR, across the 8 scRNA-seq data sets. It is worth noting that there is still a significant gap between SinNNLR and ScLSTM. The results indicate that our proposed ScLSTM performs exceptionally well on clustering tasks.

Fig. 3
figure 3

ACC, ARI, BAS, and NMI heatmaps of Corr, POCR, SIMLR, SinNLRR, SNN-cliq, ZIFA, and ScLSTM on 8 single-cell data sets

Furthermore, we compared ScLSTM with three different clustering algorithms: K-means (KM), spectral clustering (SC), and agglomerative clustering (AG). It should be noted that the overall structure of each variant algorithm still adopts the architecture of ScLSTM. Table 2 shows the results of three indicators (NMI, ACC, ARI, and BAS) on 8 data sets using three clustering methods. We can observe that agglomerative clustering almost achieves the best results, followed by K-means.

Table 2 The influence of KM, SC, and AG on the ScLSTM model is evaluated on eight data sets
Fig. 4
figure 4

A The visualization results of ScLSTM on eight data sets with TSNE and annotated each cell type; B The visualization results of eight data sets directly after agglomerative clustering

Calculate the number of clusters

The number of cell types in biological research is often unknown. Therefore, determining the number of cell types within a data set is crucial for single-cell clustering method. Agglomerative clustering is implemented usingscikit learn’s function (sklearn.cluster.AgglomerativeClustering). Since POCR and SIMLR require the specification of the number of cell types, we supplemented the evaluation of the number of clusters with CORR, SNN-Cliq, and SinNLRR on eight data sets. Table 3 presents the results of four models for predicting the number of cell types across these eight data sets: ScLSTM, SNN-Cliq, SinNLR, and CORR. In addition, we used ScLSTM to predict the number of clusters to evaluate the influence of KM, SC, and AG on the ScLSTM model across eight data sets using the predicted cluster number. The experimental results are shown in Table 4.

Table 3 ScLTM, SNN-Cliq, SinNLRR, and CORR estimate the number of clusters on eight datasets
Table 4 The influence of KM, SC, and AG on the ScLSTM model is evaluated on eight data sets using the predicted cluster number

Compared different siamese networks

To verify the effect of siamese LSTM on the effectiveness of the ScLSTM model, we propose a variant model, sigLSTM, for comparison. The sigLSTM removes one of the two identical LSTMs from the Siamese network portion of ScLSTM. The sigLSTM use the extracted feature matrix to obtain the final results through agglomerative clustering. We compare sigLSTM and ScLSTM on the 8 real scRNA-seq data sets listed in Table 1. These data sets are derived from different species and sequencing platforms, and they vary in data scale. The experimental results are presented in Table 5, which displays the optimal outcomes of ScLSTM and sigLSTM after 50 computation repeats for fourevaluation metrics (ARI, NMI, ACC, and BAS) on 8 distinct data sets. Across the Biase, Yan, Goolam, Pollen, Li, Kolodziejczyk, and Usoskin data sets, the ScLSTM model consistently outperforms the sigLSTM model. In conclusion, the ScLSTM model demonstrates superior performance across all 8 real scRNA-seq data sets generated by various single-cell sequencing platforms.

Table 5 The ARI, NMI, ACC, and BAS of sigLSTM and ScLSTM models are evaluated on eight data sets

The performance of ScLSTM on scRNA-seq data

Visual comparison of ScLSTM on eight datasets

Fig. 5
figure 5

Visualization of ScLSTM results on sc_10x_5cl and Jurkat datasets

Fig. 6
figure 6

ARI, NMI, ACC and BAS results of sc_10x_5cl and Jurkat datasets based on ScLSTM model on three clustering algorithms (KM, SC, and AG)

Biologists utilize scRNA-seq data visualization techniques to identify cell subpopulations,with t-distribution Stochastic Neighbor Embedding (t-SNE) being one of the most commonly employed methods. To demonstrate ScLSTM's ability to distinguish cell types more intuitively, we visualized the results of ScLSTM on 8 data sets with t-SNE and annotated each cell type (Fig. 4A). As a comparison group, we visualized the results of 8 data sets directly after agglomerative clustering (Fig. 4B). It is evident from the visualizations that ScLSTM can effectively distinguish different cell types in each dataset, and the boundaries are clearly defined, allowing for intuitive categorization. In contrast, Fig. 4B shows that agglomerative clustering alone struggles to distinguish different categories. In the 8 data sets, various cell types are mixed to varying degrees, and not all cell categories can be distinguished in the Biase, Usoskin, and Goolam data sets.

Performance of ScLSTM on 10x Genomics dataset

On small-scale data sets (ranging from 50 to 800 cells), ScLSTM consistently demonstrates superior performance. To highlight ScLSTM's capability to identify cell types in large-scale single-cell RNA-seq data containing thousands of cells, we conducted experiments on two publicly available data sets obtained from the 10x Genomics website. The Jurkat data set includes both Jurkat and HEK293T cell lines. The Jurkat cell line comprises 3258 293T cells, and the HEK293T cell line includes 2885 Jurkat cells, both sequenced using UMIs and Drople-based protocols. The sc_10x_5cl dataset was sequenced using the 10x Chromium Genomics protocol. Figure 5 displays the 3D visual classification results of the two data sets after ScLSTM processing,showcasing clear boundaries between different cell categories with no confounding elements. Figure 6 presents the ARI, NMI, ACC, and BAS results of three different clustering methods on sc_10x_5cl and Jurkat data sets. It is evident that AG yields the best results. These findings demonstrate that ScLSTM performs exceptionally well on large-scale data.

Human breast cancer cell type detection with ScLSTM

To demonstrate the robustness and generalizability of ScLSTM, we apply it to the breast cancer single-cell data set supplied by Chung et al. [12]. This dataset consists of 515 cells from 11 breast cancer patients, representing 4 types of breast cancer subtypes: luminal A, luminal B, HER2, and triple negative breast cancer (TNBC). Pathological testing confirmed the following four breast cancer markers: ER-positive (BC01 and BC02; luminal A), ER/HER2-positive (BC03; luminal B), HER2-positive (BC04, BC05, and BC06; HER2), and triple-negative (BC07-BC11; TNBC) invasive ductal carcinoma. In addition, regional lymph nodes were collected from the luminal B (BC03LN) sample and a triple-negative breast cancer (BC07LN) sample.

The results obtained by ScLSTM by Agglomerative Clustering are visualized using t-SNE. As depicted in Fig. 7, ScLSTM generates 6 clusters, each corresponding to a different cell type. ScLSTM separates different cell types into distinct clusters. Importantly, ScLSTM separates BC07LN and BC03LN from other cell types and maintains a clear boundary between them. It can be observed from Fig. 7 that although there are distinct  boundaries between BC03LN and BC03, the two block structures are very close, which is related to the siamese network structure of ScLSTM. Overall, all visualizations demonstrate  that the siamese LSTM structures excel at learning cell-type features and improving clustering performance.

Fig. 7
figure 7

Visualization of the cells in breast cancer dataset (GSE75688) based on the t-SNE learned similarity matrix from ScLSTM


In single-cell transcriptome analysis, grouping individual cells based on gene expression data is valuable for characterizing cellular composition in tissues, distinguishing developmental stages, and understanding the pathological process of tissues. Differences in the distribution of scRNA-seq data and an unequal number of cell types are significant obstacles to detecting single cell types. To address this, we propose a meta-learning model called ScLSTM, composed of an LSTM network with an improved sigmoid kernel framework. ScLSTM describes the features of each cell using various hidden units, and the cell semantic representation obtained by mining gene expression sequences is better appropriate for describing cell types.

We conducted validation experiments and in-depth analyses on 8 scRNA-seq data sets and two 10x Genomics from various single-cell sequencing platforms. Using default hyperparameters, ScLSTM consistently yields the highest ARIs, NMIs, ACCs, and BASs compared to other clustering methods. The experiment results illustrate the viability and efficacy of our proposed ScLSTM for detecting single cell types in datasets with varying scales and distributions. Unlike the Corr, POCR, and SIMLR methods, ScLSTM can directly compare the slight differences between two cells and magnify the difference between cell types.

ScLSTM’s success is contingent on its novel feature extraction and integration style. Unlike other approaches, ScLSTM locates the perspective of feature extraction at the difference between two cells, ensuring that the properties of all individual cells help distinguish between cell types. Compared to conventional single-cell type detection methods, the feature extraction from ScLSTM is more accurate and closer to the target classification. As scRNA-seq demonstrates its distinct benefits in biomedicine, we expect that ScLSTM will provide researchers with accurate single-cell clustering services.

Availability of data and materials

For ScLSTM, we conduct experiments on 8 different real human and mouse cells scRNA-seq data sets. The datasets generated and analysed dyring the current study are available in the Gene Expression Omnibus (GEO) repository, which were taken with URL Yan, Camp, Biase, Usoskin, and Li data sets were downloaded from the NCBI GEO database (access IDs are GSE36552, GSE81252, GSE57249, GSE59739, and GSE81861), and Pollen data set was downloaded from the NCBI SRA database (access IDs is SRP041736). Goolam and Kolodziejczyk data sets were downloaded from the ArrayExpress database (access IDs are E-MTAB-3321 and E-MTAB-2600). The sc_10x_5cl data set GEO accession number is GSE118767. The Jurkat data set downloaded from and The source codes of ScLSTM and the processed data sets are available at


  1. Elowitz MB, Levine AJ, Siggia ED, et al. Stochastic gene expression in a single cell. Science 2002;297(5584):1183–6.

  2. Butler A, Hoffman P, Smibert P, et al. Integrating single-cell transcriptomic data across different conditions, technologies, and species. Genome biology 2018;36(5):411–20.

  3. Jiang H, Sohn LL, Huang H, Chen L. Single cell clustering based on cell-pair differentiability correlation and variance analysis. Bioinformatics. 2018;34(21):3684–94.

    Article  CAS  PubMed  Google Scholar 

  4. Liu Z, Liu F, Hong C, Gao M, Chen Y-PP, Liu S, Zhang W. Detection of cell types from single-cell RNA-seq data using similarity via kernel preserving learning embedding. In: 2019 IEEE international conference on bioinformatics and biomedicine (BIBM). IEEE; 2019. p. 451–457.

  5. Zheng R, Li M, Liang Z, Wu F-X, Pan Y, Wang J. SinNLRR: a robust subspace clustering method for cell type detection by non-negative and low-rank representation. Bioinformatics. 2019;35(19):3642–50.

    Article  CAS  PubMed  Google Scholar 

  6. Wang B, Zhu J, Pierson E, Ramazzotti D, Batzoglou S. Visualization and analysis of single-cell RNA-seq data by kernel-based similarity learning. Nat Methods. 2017;14(4):414–6.

    Article  CAS  PubMed  Google Scholar 

  7. Cui Y, Zhang S, Liang Y, Wang X, Ferraro TN, Chen Y. Consensus clustering of single-cell RNA-seq data by enhancing network affinity. Brief Bioinf. 2021;22(6):bbab236.

    Article  CAS  Google Scholar 

  8. Tian T, Zhang J, Lin X, Wei Z, Hakonarson H. Model-based deep embedding for constrained clustering analysis of single cell RNA-seq data. Nat Commun. 2021;12(1):1–12.

    Article  Google Scholar 

  9. Zhou X, Chai H, Zeng Y, Zhao H, Yang Y. scAdapt: virtual adversarial domain adaptation network for single cell RNA-seq data classification across platforms and species. Brief Bioinf. 2021;22(6):bbab281.

    Article  CAS  Google Scholar 

  10. Song Q, Su J, Zhang W. scGCN is a graph convolutional networks algorithm for knowledge transfer in single cell omics. Nat Commun. 2021;12(1):1–11.

    CAS  Google Scholar 

  11. Jiang H, Huang Y. An effective drug-disease associations prediction model based on graphic representation learning over multi-biomolecular network. BMC Bioinf 2022.

  12. Chung W, Eum HH, Lee H-O, Lee K-M, Lee H-B, Kim K-T, Ryu HS, Kim S, Lee JE, Park YH, et al. Single-cell RNA-seq enables comprehensive tumour and immune cell profiling in primary breast cancer. Nat Commun. 2017;8(1):1–12.

    Article  Google Scholar 

  13. Yan L, Yang M, Guo H, Yang L, Wu J, Li R, Liu P, Lian Y, Zheng X, Yan J, et al. Single-cell RNA-seq profiling of human preimplantation embryos and embryonic stem cells. Nat Struct Mol Biol. 2013;20(9):1131.

    Article  CAS  PubMed  Google Scholar 

  14. Goolam M, Scialdone A, Graham SJ, Macaulay IC, Jedrusik A, Hupalowska A, Voet T, Marioni JC, Zernicka-Goetz M. Heterogeneity in oct4 and sox2 targets biases cell fate in 4-cell mouse embryos. Cell. 2016;165(1):61–74.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  15. Camp JG, Sekine K, Gerber T, Loeffler-Wirth H, Binder H, Gac M, Kanton S, Kageyama J, Damm G, Seehofer D, et al. Multilineage communication regulates human liver bud development from pluripotency. Nature. 2017;546(7659):533–8.

    Article  CAS  PubMed  Google Scholar 

  16. Pollen AA, Nowakowski TJ, Shuga J, Wang X, Leyrat AA, Lui JH, Li N, Szpankowski L, Fowler B, Chen P, et al. Low-coverage single-cell mRNA sequencing reveals cellular heterogeneity and activated signaling pathways in developing cerebral cortex. Nat Biotechnol. 2014;32(10):1053.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  17. Biase FH, Cao X, Zhong S. Cell fate inclination within 2-cell and 4-cell mouse embryos revealed by single-cell RNA sequencing. Genome Res. 2014;24(11):1787–96.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  18. Usoskin D, Furlan A, Islam S, Abdo H, Lönnerberg P, Lou D, Hjerling-Leffler J, Haeggström J, Kharchenko O, Kharchenko PV, et al. Unbiased classification of sensory neuron types by large-scale single-cell RNA sequencing. Nat Neurosci. 2015;18(1):145–53.

    Article  CAS  PubMed  Google Scholar 

  19. Kolodziejczyk AA, Kim JK, Tsang JC, Ilicic T, Henriksson J, Natarajan KN, Tuck AC, Gao X, Bühler M, Liu P, et al. Single cell RNA-sequencing of pluripotent states unlocks modular transcriptional variation. Cell Stem Cell. 2015;17(4):471–85.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  20. Li H, Courtois ET, Sengupta D, Tan Y, Chen KH, Goh JJL, Kong SL, Chua C, Hon LK, Tan WS, et al. Reference component analysis of single-cell transcriptomes elucidates cellular heterogeneity in human colorectal tumors. Nat Genet. 2017;49(5):708–18.

    Article  CAS  PubMed  Google Scholar 

  21. Kuhn HW. The Hungarian method for the assignment problem. Naval Res Logist Q. 1955;2(1–2):83–97.

    Article  Google Scholar 

  22. Xu C, Su Z. Identification of cell types from single-cell transcriptomes using a novel clustering method. Bioinformatics. 2015;31(12):1974–80.

    Article  CAS  PubMed  PubMed Central  Google Scholar 

  23. Pierson E, Yau C. Zifa: dimensionality reduction for zero-inflated single-cell gene expression analysis. Genome Biol. 2015;16(1):1–10.

    Article  Google Scholar 

Download references


The authors would like to thank anonymous reviewers for providing valuable comments for our article.


This work was supported by the National Natural Science Foundation of China (62276191) and the Fundamental Research Funds for the Central Universities (HUST:2022JYCXJJ007).

Author information

Authors and Affiliations



HJ and QL conceived the study. HJ and YH developed the methodology. QL implemented the methodology. YH analysed the results. BF conceived the supplementary experiments. HJ conducted the supplementary experiments and analyzed the experimental results. YH wrote and reviewed the manuscript.

Corresponding author

Correspondence to Yabing Huang.

Ethics declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

Not applicable.

Competing interests

The authors declare that they have no competing interests.

Additional information

Publisher’s Note

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

Rights and permissions

Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Jiang, H., Huang, Y., Li, Q. et al. ScLSTM: single-cell type detection by siamese recurrent network and hierarchical clustering. BMC Bioinformatics 24, 417 (2023).

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: