
==== Front
Genome Biol
Genome Biol
Genome Biology
1474-7596
1474-760X
BioMed Central London

39237934
3367
10.1186/s13059-024-03367-8
Method
DeepKINET: a deep generative model for estimating single-cell RNA splicing and degradation rates
Mizukoshi Chikara m-chikara@nagoya-u.ac.jp

12
Kojima Yasuhiro yakojim@ncc.go.jp

34
Nomura Satoshi 1
Hayashi Shuto 4
Abe Ko 4
Shimamura Teppei shimamura.csb@tmd.ac.jp

14
1 https://ror.org/04chrp450 grid.27476.30 0000 0001 0943 978X Division of Systems Biology, Graduate School of Medicine, Nagoya University, Aichi, Japan
2 https://ror.org/008zz8m46 grid.437848.4 0000 0004 0569 8970 Nagoya University Hospital, Aichi, Japan
3 grid.272242.3 0000 0001 2168 5385 Laboratory of Computational Life Science, National Cancer Center Research Institute, Tokyo, Japan
4 https://ror.org/051k3eh31 grid.265073.5 0000 0001 1014 9130 Department of Computational and Systems Biology, Medical Research Institute, Tokyo Medical and Dental University, Tokyo, Japan
6 9 2024
6 9 2024
2024
25 22914 12 2023
4 8 2024
© The Author(s) 2024
2024
https://creativecommons.org/licenses/by-nc-nd/4.0/ Open Access This article is licensed under a Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License, which permits any non-commercial use, sharing, 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 you modified the licensed material. You do not have permission under this licence to share adapted material derived from this article or parts of it. The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by-nc-nd/4.0/.
Messenger RNA splicing and degradation are critical for gene expression regulation, the abnormality of which leads to diseases. Previous methods for estimating kinetic rates have limitations, assuming uniform rates across cells. DeepKINET is a deep generative model that estimates splicing and degradation rates at single-cell resolution from scRNA-seq data. DeepKINET outperforms existing methods on simulated and metabolic labeling datasets. Applied to forebrain and breast cancer data, it identifies RNA-binding proteins responsible for kinetic rate diversity. DeepKINET also analyzes the effects of splicing factor mutations on target genes in erythroid lineage cells. DeepKINET effectively reveals cellular heterogeneity in post-transcriptional regulation.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-024-03367-8.

Keywords

Single-cell RNA sequencing (scRNA-seq)
RNA splicing
RNA degradation
Splicing kinetics
Transcriptome dynamics
RNA-binding proteins
RNA velocity
Neural network
Variational autoencoder (VAE)
Deep generative model
Dimensionality reduction
Cell differentiation
Metabolic labeling
http://dx.doi.org/10.13039/501100001691 Japan Society for the Promotion of Science 22H04925 23H04938 Shimamura Teppei http://dx.doi.org/10.13039/100009619 Japan Agency for Medical Research and Development JP22ama221215 JP22ama221501 JP22wm0425007 JP23wm0325068 Shimamura Teppei http://dx.doi.org/10.13039/100009619 Japan Agency for Medical Research and Development P23tm0424226 Shimamura Teppei http://dx.doi.org/10.13039/501100020963 Moonshot Research and Development Program JPMJMS2025 Shimamura Teppei http://dx.doi.org/10.13039/501100007951 Tokyo Medical and Dental University issue-copyright-statement© BioMed Central Ltd., part of Springer Nature 2024
==== Body
pmcBackground

Messenger RNA (mRNA) splicing and degradation play essential roles in precise gene expression regulation. These processes are vital for accurate utilization of genetic information within cells. Inappropriate splicing can lead to production of dysfunctional proteins, potentially resulting in severe implications for fundamental cellular functions. Recent studies have established that abnormal mRNA splicing and degradation are closely associated with development and progression of diseases such as cancer [1, 2].

Several methodologies are available to estimate mRNA splicing and degradation rates, each with its own limitations and challenges. Metabolic labeling methods [3, 4] are used to estimate the synthesis and degradation rates in genome-wide RNA metabolism by integrating RNA metabolic labeling with cell-specific splicing kinetics. However, owing to the necessity of specific metabolic labeling, this approach is limited and cannot be applied as readily as conventional scRNA-seq data. Combination of scRNA-seq data with the RNA velocity theory [5] was introduced to model the dynamic processes of mRNA in individual cells. However, this approach has been criticized for assuming uniform kinetic rates across cells, which may cause misrepresentation of true biological variation. While transcription rates have been modeled to account for cell-to-cell variability, methods such as scVelo [6] and VeloVI [7] have assumed uniform splicing and degradation rates for each gene. A novel relay velocity model [8] utilizes neighboring cell information and leverages deep neural networks to estimate cell-specific kinetic rates. However, its primary intention is to refine the RNA velocity, leaving questions regarding the accuracy of the kinetic rates for each cell.

In light of these challenges, we introduced DeepKINET (a deep generative model with single-cell RNA kinetics), an advanced analysis framework based on deep generative modeling. This framework uses deep generative model-driven cell states in scRNA-seq data to accurately estimate single-cell splicing and degradation kinetics. Our method aims to reveal the heterogeneity in splicing and degradation rates across cells, enabling to elucidate post-transcriptional regulatory mechanisms mediated by factors such as RNA-binding proteins.

We demonstrate that DeepKINET can estimate mRNA splicing and degradation rates with greater precision than existing methods, as evidenced by simulated and metabolic labeling experimental data. Moreover, we demonstrate its robustness against dropouts. By applying DeepKINET to a forebrain dataset, we analyzed whether genes governed by the same RNA-binding proteins have equivalent trends in their splicing and degradation rates, and we identified the biological functions of these RNA-binding proteins. Furthermore, when applied to breast cancer data, DeepKINET revealed splicing and degradation anomalies related to cancer metastasis; we provide specific examples. In addition, we analyzed the effects of mutations in a splicing factor on the target genes in erythroid lineage cells. The results enhance our understanding of mRNA splicing and degradation processes and help to elucidate underlying molecular mechanisms and potential therapeutic targets.

Results

Conceptual view of DeepKINET

Figure 1 presents a clear overview of the conceptual framework of DeepKINET. This method processes both spliced and unspliced mRNA counts from scRNA-seq data and subsequently generates comprehensive kinetic rates across genes, including splicing and degradation rates, at single-cell resolution. DeepKINET addresses heterogeneity in kinetic rates spanning genes and cells, which is ignored by existing methods [5, 6].Fig. 1 Overview of DeepKINET. a Overview of our method for estimating single-cell transcriptome dynamics from latent variables. DeepKINET receives scRNA-seq data that have unspliced and spliced counts and outputs kinetic rates at the single-cell level. DeepKINET provides biologically meaningful insights by accounting for cellular heterogeneity in kinetic rates, which is ignored by existing methods. For example, DeepKINET can be used to classify genes by their kinetic rates, find genes that show significant rate variation among cell populations, and identify RNA-binding proteins involved in splicing and degradation. b Estimated RNA velocity by DeepKINET in the mouse pancreas dataset visualized on Uniform Manifold Approximation and Projection (UMAP) embedding. The direction of transition in latent space is plotted in 2D coordinates in the same way as scvelo. Trajectory inference by PAGA [10] was performed using RNA velocity from DeepKINET. c Expression, splicing rate, and degradation rate at the single-cell level projected on UMAP embedding. DeepKINET estimates splicing and degradation rates for each cell based on the RNA velocity equation and cell states. The colors of the points indicate the gene expression, the splicing rate, and the degradation rate per cell

DeepKINET uses a deep generative model of mature and immature transcripts based on an RNA velocity equation. This enables optimization in which the splicing and degradation rates are adjusted according to the cell state. First, we use a variational autoencoder (VAE) to model stochastic transitions within the latent cell state space, similar to that in our previous study [9]. DeepKINET assumes that the kinetic parameters for each cell are obtained from transformation of the latent cell state by the neural network. We optimized both cell state dynamics and kinetic parameters to align with the observed mature and immature transcript levels, following the RNA velocity equation.

Beyond kinetic rate heterogeneity estimation across genes and cells, DeepKINET offers the following: (1) gene clustering based on kinetic rates, which enables identification of genes with analogous rate patterns; (2) identification of genes exhibiting significant rate variations by comparing different cell populations; and (3) detection of RNA-binding proteins that influence splicing and degradation rates of their associated targets.

DeepKINET not only delivers refined insights into RNA kinetics but also serves as a springboard for in-depth molecular studies, promising deeper comprehension and demystification of the complex regulatory mechanisms guiding cellular kinetics. It is accessible as a user-friendly open-source Python package with comprehensive documentation at https://github.com/3254c/DeepKINET.

Simulated data to demonstrate accuracy and superiority of DeepKINET

We used simulated data to evaluate the accuracy of the kinetic rates estimated using the DeepKINET software. Simulated data were generated using SERGIO [11], which uses gene regulatory networks and RNA velocity equations to generate the scRNA-seq data. We generated scRNA-seq count data for each cell cluster with different splicing and degradation rates.

We applied DeepKINET to each simulated dataset and confirmed that it predicted the correct direction of differentiation (Additional file 1: Fig. S1a). We then estimated the kinetic rates for each single cell (Fig. 2a), averaged them over each cell cluster, and calculated the correlation coefficient using the set value (Fig. 2b). We found positive correlations across various dropout scenarios. Therefore, we concluded that our method is robust against data sparsity. In addition, we compared the accuracy of our method with that of cellDancer [8] and DeepVelo [12], existing methods for estimating kinetic rates at the single-cell level. Although Velocyto [5] and scVelo [6] are widely used in RNA velocity analysis, we did not include them in our comparison as these methods assume uniform kinetic rates across cells. cellDancer showed positive correlations in splicing rates, whereas DeepVelo showed negative correlations in splicing rates. Both methods were less accurate than DeepKINET. Furthermore, cellDancer showed negative correlations in degradation rates. DeepVelo showed positive correlations with respect to degradation rates, but DeepKINET had higher correlations. Subsequent simulations were conducted using varying numbers of cells. For these simulations, we used the default dropout rates. We applied DeepKINET to each simulated dataset and computed correlation coefficients for the set values. DeepKINET could accurately estimate the kinetic rates, even for small numbers of cells (Fig. 2c). On the other hand, cellDancer also showed positive correlations in splicing rate estimation accuracy, but it was less accurate than DeepKINET and required more cells until the estimation accuracy stabilized. Furthermore, DeepVelo was unable to make accurate estimates even when the number of cells was increased. In degradation rates, cellDancer consistently failed to make correct estimates, and the accuracy decreased as the number of cells increased. DeepVelo showed positive correlations, but the estimation precision was still lower than DeepKINET.Fig. 2 DeepKINET is robust to dropout rates and cell numbers in simulated data, and its performance exceeds that of cellDancer and DeepVelo. a Visualization of the UMAP embedding of the expression, set kinetic rates, and estimated kinetic rates. The gene with the highest correlation in splicing rate and the gene with the highest correlation in degradation rate are shown. To prevent extreme values from affecting the visualization, the minimum or maximum value of the top 1% was forced to the 1% and 99% quantile values. b Scatter plot of correlation coefficient averages of splicing rates and degradation rates for each dataset. Ten datasets were generated for each of the 20 different generation conditions. We applied DeepKINET, cellDancer and DeepVelo once to each dataset and calculated the correlation coefficient between the set rates and the estimated rates by each method. DeepKINET’s accuracy exceeds that of cellDancer and DeepVelo. c Box plot of correlation coefficient averages when varying the number of cells in a cluster. Ten datasets were generated for each of the 14 different generation conditions. DeepKINET always had a positive correlation coefficient and outperformed cellDancer and DeepVelo

These validations confirmed the accuracy of the splicing and degradation rates estimated by DeepKINET, marking a clear advancement over the kinetic parameter estimation capabilities of cellDancer and DeepVelo. Notably, the accuracy of splicing rate estimation by cellDancer appeared to increase slowly as the number of cells increased, implying a requirement for larger datasets than those required by DeepKINET for accurate predictions. A detailed exposition of genes that were successfully estimated and those that were not is shown in Additional file 1: Fig. S1b, S1c.

Due to the presence of multiple unknowns in the RNA velocity equation for spliced mRNA, the solution of the splicing and degradation rates for each cell may be underdetermined, which could lead to correlation between the estimated kinetic rates. To assess whether each method can estimate these parameters independently, we investigated the correlation between the estimated splicing rates and degradation rates in each simulation dataset. We observed that the correlation between splicing and degradation rates estimated by DeepKINET was lower than that of cellDancer and DeepVelo (Additional file 1: Fig. S1d). This suggests that DeepKINET can estimate splicing and degradation rates more independently compared to other methods. In contrast, cellDancer and DeepVelo exhibited relatively high correlations, indicating that these methods have difficulty in separately considering splicing and degradation processes.

Accuracy of DeepKINET for real data evaluated using metabolic labeling data

We next evaluated the accuracy of DeepKINET for real data using multicellular-level kinetic rates derived from metabolic labeling experimental data. The values obtained from the metabolic labeling experiments depended on the assumptions of the mathematical model used and did not represent the perfect ground truth. Nevertheless, the temporal resolution inherent in the metabolic experimental data lost in scRNA-seq provides a benchmark from which to assess the similarity to extrapolated kinetic rates. Li et al. [8] used single-cell EU-labeled RNA sequencing (scEU-seq) [3] to qualitatively assess the accuracy of cellDancer, which is limited to cell cycle genes.

Using the same scEU-seq cell cycle dataset, we evaluated the accuracy of DeepKINET. scEU-seq methodology can be used to estimate multicellular-level kinetic rates by observing temporal variations in the fraction of 5-ethynyl-uridine(EU)-labeled mRNA. Battich et al. did not differentiate between unspliced and spliced mRNAs when modeling mRNA metabolism. Conversely, Dynamo [13] can estimate kinetic rates, including splicing rates, by accounting for splicing events in the scEU-seq data. We partitioned the cell cycle dataset into PULSE and CHASE experimental categories, each distinctly modeling mRNA metabolism. In the Pulse experiment, the EU incubation time differs for each cell. In the Chase experiment, EU incubation is performed under the same conditions, followed by washing time with uridine, which varies depending on the cells. We divided the cells into six to thirty clusters, each containing an equal number of cells across the cell cycle trajectory. Dynamo was used to determine the splicing and degradation rates for each cluster.

Next, we estimated the RNA velocity using DeepKINET and confirmed that the estimated future states of individual cells followed the order of the cell cycle (Fig. 3a, Additional file 1: S2a). We then estimated the single-cell splicing and degradation rates, averaged them across clusters and calculated the correlation coefficient using the kinetic rates determined using Dynamo. Our method showed positive correlations in both PULSE and CHASE experiments, outperformed cellDancer and DeepVelo in terms of accuracy in degradation rates, and demonstrated comparable performance in accuracy in splicing rates (Fig. 3b, Additional file 1: S2b). Notably, the PULSE experimental data were considered more reliable because the proportion of cells in different cell cycles was constant. Regarding the degradation rates, cellDancer showed negative correlation in both experiments, and DeepVelo showed negative correlation in the CHASE experiment. DeepKINET showed no negative correlation in any setting.Fig. 3 DeepKINET is also accurate for real data and outperforms cellDancer and DeepVelo. a The clusters and velocities of the PULSE data were visualized on the pre-defined embedding based on the Geminin-GFP and Cdt1-RFP signals. Geminin-GFP and Cdt1-RFP signals were used in Battich et al.’s paper [3] to estimate cell cycle score. The cells were divided into cell clusters based on cell cycle position. The PULSE data showed an even distribution of cells with respect to cell cycle positions compared to the CHASE data (S2a). DeepKINET is able to estimate the correct direction along the cell cycle. b Box plot of correlation coefficient averages between estimated rates by Dynamo and estimated rates by DeepKINET, cellDancer, and DeepVelo using the PULSE experimental data. A total of 100 estimations were performed by each of DeepKINET, cellDancer, and DeepVelo. DeepKINET showed positive correlations, outperformed cellDancer and DeepVelo in terms of accuracy in degradation rates, and demonstrated comparable performance in accuracy in splicing rates. cellDancer showed negative correlations in degradation rates. c Heatmaps of splicing rates (left) and degradation rates (right). To prevent extreme values from affecting the visualization, the minimum or maximum value of the top 1% was forced to the 1% and 99% quantile values. The genes were clustered by splicing and degradation rates and sorted by their clusters. The cells were sorted by cell cycle positions. The correlation coefficients between splicing and degradation rates for each gene are indicated by colored bars. The genes related to the cell cycle [14] are also shown in white color. Cluster 9 has a large number of genes related to the cell cycle. d Gene Ontology (GO) terms enriched in the gene list belonging to cluster 9 obtained by g:Profiler. One-thousand genes in this analysis were used as background. e Genes with different correlations between splicing and degradation rates. DeepKINET can extract genes by the value of the correlation between splicing and degradation rates. The minimum or maximum value of the top 1% was forced to the 1% and 99% quantile values

Using the PULSE experimental data, we estimated the splicing and degradation rates for each cell and clustered the genes using these rates (Fig. 3c). We then derived the correlation coefficients between the splicing and degradation rates. Genes related to the cell cycle were concentrated in one cluster, and related terms were detected using Gene Ontology (GO) analysis (Fig. 3d). Finally, we classified the genes using the correlation coefficients between splicing and degradation rates (Fig. 3e).

Additionally, we evaluated the accuracy of DeepKINET using another metabolic labeling dataset. Single-cell metabolically labeled new RNA tagging sequencing (scNT-seq) [4] enables the estimation of transcription and degradation rates by distinguishing between old and new transcriptomes in the same cell. We applied DeepKINET and other methods to the hematopoiesis dataset [13]. DeepKINET estimated the differentiation trajectory of hematopoietic cells, which was consistent with the results in the Dynamo paper (Additional file 1: Fig. S3a). We applied Dynamo to estimate the degradation rates of two cell populations (Additional file 1: Fig. S3b) and compared these estimates with those obtained from other methods. The degradation rates estimated by DeepKINET correlated with the values estimated by Dynamo and demonstrated superior performance compared to cellDancer and DeepVelo (Additional file 1: Fig. S3c). These results further support the accuracy of DeepKINET’s estimations.

DeepKINET to investigate functions of RNA-binding proteins and RNA-binding proteins that regulate gene clusters

We applied DeepKINET to a forebrain dataset [5] to examine the functions of RNA-binding proteins. DeepKINET can classify genes based on their kinetic rates and identify RNA-binding proteins that govern these clusters. Additionally, DeepKINET can determine whether an RNA-binding protein regulates the splicing or degradation of its target genes.

First, we confirmed that the direction of RNA velocity estimated by DeepKINET was consistent with the known trajectories of cell differentiation (Fig. 4a). We then used DeepKINET to estimate the single-cell splicing and degradation rates and used these rates separately to cluster the genes. By clustering genes using either splicing rates or degradation rates independently, it becomes possible to discern which process an RNA-binding protein contributes to, thereby providing insights into its functional role in post-transcriptional regulation. We examined whether the gene clusters by kinetic rates matched the gene list of RNA-binding protein targets using Fisher’s exact test (Fig. 4b). We found clusters that matched the target gene lists, indicating that genes regulated by the same RNA-binding protein have similar splicing and degradation rate changes.Fig. 4 RNA-binding protein analysis of the forebrain dataset by DeepKINET. a PAGA trajectory inference of forebrain dataset using DeepKINET’s velocity estimates. b Dot heatmap showing the association of each RNA-binding protein (RBP) targets with each gene cluster. The genes were clustered using splicing and degradation rates separately, and Fisher’s exact test was used to determine if a list of RNA-binding protein target genes were enriched in a particular cluster. The colors indicate the corrected p-values for Fisher’s exact test. The circle size indicates the ratio of the proportion of RNA-binding protein targets in the cluster to the proportion of RNA-binding protein targets in all genes. c Joint plot of the mean correlation coefficient between RNA-binding protein expression levels and the splicing and degradation rates of each target or non-target. Compared with non-target genes, target genes have higher correlations with the expression of RNA-binding proteins. A significant difference was indicated by the Levene’s test. d Box plots show correlation coefficients between RBFOX1 and RBFOX2 expression and the splicing rates of each target or non-target gene. The green dot represents the average value. A significant difference was indicated by the one-sided unpaired t-test. e Visualization of the UMAP embedding of the expression of RBFOX1 and RBFOX2 and the splicing rates of target genes that are highly correlated with RBFOX1 and RBFOX2 expression

Next, we examined the relationship between the expression levels of each RNA-binding protein and the splicing and degradation rates of the target genes. We calculated the average correlation coefficients for both target and non-target genes for all remaining RNA-binding proteins from expression preprocessing and the observed absolute values of correlation coefficients of target genes were significantly greater than those of non-target genes (Fig. 4c). This suggests that DeepKINET accurately reflects the regulatory roles of RNA-binding proteins with respect to their target genes. Further analysis of the highly variable genes that substantially affected the kinetic rates of their targets revealed that the expression levels of RBFOX1 and RBFOX2 correlated with the target splicing rates (Fig. 4d), which is in agreement with established research identifying RBFOX1 and RBFOX2 as regulators of mRNA splicing [15]. Therefore, DeepKINET demonstrated proficiency in deducing the contributions of RNA-binding proteins to splicing and degradation within the dataset, as well as in identifying genes that are potentially regulated by specific RNA-binding proteins (Fig. 4e).

DeepKINET reveals heterogeneity in cancer cell populations

Next, we applied DeepKINET to breast cancer data to identify genes with significant changes in kinetic rates and RNA-binding proteins that exhibit distinct functions across different cell populations. Previous studies have highlighted the critical roles of splicing and degradation abnormalities in cancer development and progression [1, 2]. Additionally, the significant involvement of RNA-binding proteins in cancer has been well documented [16, 17]. Cell Ranger [18] and Velocyto [5] were used to create matrices of the spliced and unspliced breast cancer data [19].

We applied DeepKINET to malignant epithelial cells from the breast cancer data (Fig. 5a) and confirmed that the estimated velocities were in the direction from primary cells to metastatic cells (Additional file 1: Fig. S4a). We then estimated the single-cell kinetic rates and identified genes that exhibited marked differences in their splicing or degradation rates when primary cells were compared with metastatic cells (Fig. 5b, c). Among these, KDM6A [20], PGR [21], PIK3CA [22], PRKAA1 [23], TPM2 [24], TP63 [25], USP9X [26], and TIMP2 [27] have been implicated in breast cancer metastasis. These variations in the kinetic rates may play a pivotal role in metastasis.Fig. 5 DeepKINET can identify kinetics changes involved in metastasis using breast cancer data. a Visualization of UMAP embedding of the velocities estimated by DeepKINET and pre-defined classifications for malignant epithelial cells of breast cancer. There are 15,269 primary cells and 642 metastatic cells. The velocities indicate the direction from the primary cancer to the metastatic cancer. b Genes with large changes in their splicing rates (left) and degradation rates (right) between primary and metastatic cells as determined using t-test. These genes include those involved in cancer metastasis and breast cancer. c Visualization of UMAP embedding of expression levels and kinetic rates of the genes with the largest changes in their splicing or degradation rates. These genes are involved in breast cancer metastasis. To prevent the effect of extreme values in the visualization, the minimum or maximum value of the top 1% was forced to the 1% and 99% quantile values. d Bar plot of corrected p-values for genes whose slopes changed significantly between primary and metastatic cells when linear regression was performed using RBM47 expression levels and target splicing rates. Several of these genes are involved in breast cancer metastasis and metastasis of other cancers. e Scatter plot of RBM47 expression and splicing rates of CD2AP, the gene with the most slope change

Furthermore, we explored the correlation coefficients between the expression of highly variable RNA-binding proteins and the kinetic rates of their target genes. Within this dataset, the effect of RBM47 on the splicing rate of its target genes was significant (Additional file 1: Fig. S4b, S4c). Because RBM47 is involved in RNA splicing and metastasis, including that of breast cancer [28–30], this result indicates the capacity of DeepKINET to accurately reflect authentic biological processes. We also investigated whether the relationship between RBM47 and its target genes differed significantly between primary and metastatic cells. We performed linear regression on the expression of RBM47 and the splicing rates of its targets, and examined whether the slope of the regression varied significantly between primary and metastatic cells. We corrected the p-values using multiple testing corrections and extracted significantly altered genes (Fig. 5d, e). Among these genes, CTSC [31], PSD3 [32], TGFBR3 [33], and USP53 [34] are involved in breast cancer metastasis. CD2AP [35], GFRA1 [36], and EPB41 [37] are implicated in the metastasis of other cancers, but no findings on breast cancer metastasis have been reported. These findings imply that changes in the effect of RBM47 expression on the splicing rates of its target genes are associated with cellular transitions critical for cancer metastasis.

DeepKINET reveals changes in splicing due to mutations in splicing factors

Finally, we investigated changes in target splicing rates due to mutations in a splicing factor using erythroid lineage cells. The subtype of myelodysplastic syndrome (MDS), MDS-RS (MDS with ringed sideroblasts), has mutations in the splicing factor SF3B1 and is characterized by severe anemia and the accumulation of erythrocyte progenitor cells in the bone marrow. SF3B1 is responsible for connecting immature mRNAs to spliceosomes, and mutations in it lead to aberrant splicing, particularly the use of alternative 3′ splice sites [38], resulting in reduced standard transcripts [39].

We extended our model using the conditional VAE framework [40] to integrate and analyze multiple samples (see more details in the “Methods” section). We used data from Adema et al.’s [41] bone marrow mononuclear cells from seven MDS-RS patients with SF3B1 mutations and two age-matched healthy donors. We analyzed erythroid lineage cells, which are known to be damaged by MDS-RS (Fig. 6a). The results of trajectory inference were consistent with the known erythroid differentiation process (Fig. 6b).Fig. 6 Analysis of a mutated splicing factor by DeepKINET. a Patient label (left) and SF3B1 mutation status (right) on the UMAP embedding from the low-dimensional latent cell state estimated by DeepKINET. b PAGA trajectory inference using DeepKINET’s velocity estimates. DeepKINET can accurately estimate the differentiation pathway of erythroid lineage cells. c UMAP coordinates colored by the expression of SF3B1. d Scatter plot of the average splicing rate by SF3B1 mutation status (wild-type on the x-axis and mutant on the y-axis) for each gene. The p-values shown in the figure are derived from a one-sided paired t-test, which shows that cells with the SF3B1 mutation have significantly lower target splicing rates compared to cells without the mutation. The top three genes with the smallest ratio of splicing rates in SF3B1 mutated cells compared to wild-type cells (WAS, APOE, and CHAC1) were highlighted. e Box plot showing the correlation between SF3B1 expression and target or non-target gene splicing rates. In the target genes, the correlation value were significantly higher than those of the non-target genes. f Box plot showing the correlation between SF3B1 expression and target gene splicing rate. In cells with SF3B1 mutations, the correlation value was significantly lower than in cells without mutations

Our analysis revealed that the splicing rate of SF3B1 target genes in SF3B1 mutant cells was significantly lower than in healthy cells (Fig. 6d). This result suggests that DeepKINET effectively captures the changes in splicing kinetics caused by SF3B1 mutations. It should be noted that in the RNA velocity model, the splicing rate is defined as the amount of change from unspliced mRNA to spliced mRNA per unit time. In the standard quantification method such as Velocyto [5], reads are annotated as spliced mRNA if they map only to exon regions, and even if only a small amount maps to intron regions, they are annotated as unspliced mRNA. It is known that in the case of SF3B1 mutations, the usage of alternative 3′ splice sites results in the shifting of splicing sites tens of base pairs upstream compared to the canonical 3′ splice sites, causing an insertion of intronic sequence at the authentic exon junction [39]. From these considerations, it can be inferred that the use of alternative 3′ splicing sites in SF3B1 mutants reduces the amount of transcripts classified as spliced mRNA, which in turn leads to decreases the splicing rates in the RNA velocity model. The results from DeepKINET suggest that this model captures the changes in kinetics underlying these biological processes.

Additionally, the genes that show a large reduction in splicing rate due to the SF3B1 mutation included WAS, APOE, and CHAC1. It has been suggested that the deficiency of WAS, whose splicing rate was most substantially decreased by SF3B1 mutation, promotes the development of hematopoietic malignancies including MDS [42, 43]. APOE is involved in both tumor promotion and suppression [44], and CHAC1 exerts antitumor effects when its expression is increased [45]. The results from DeepKINET suggest that SF3B1 mutations may decrease the expression of these genes through the reduction of the splicing rates, potentially contributing to the development of MDS.

Furthermore, the average correlation between the expression of SF3B1 and the splicing rates of target genes was significantly higher than the average correlation between the expression of SF3B1 and the splicing rates of non-target genes (Fig. 6e), consistent with the functional characteristics of SF3B1. In addition, the average correlation between the expression of SF3B1 and the splicing rates of target genes in SF3B1 mutant cells was significantly lower than that in SF3B1 non-mutant cells (Fig. 6f). This results suggests that mutations in SF3B1 make it difficult to produce the target standard spliced mRNA and that the mutated SF3B1 does not contribute to normal splicing. These findings demonstrate that DeepKINET can capture changes in splicing of targets due to mutations in splicing factors.

Discussion

In this study, we introduced DeepKINET, a groundbreaking method for accurately estimating splicing and degradation rates at single-cell resolution. By harnessing cell state information and RNA velocity, DeepKINET advances beyond conventional models that assign static splicing and degradation rates to genes, offering dynamic and cell-specific analysis. This innovation marks the first instance in which such kinetic rates have been estimated and validated for accuracy at the single-cell level using both simulated and metabolic labeling data, thereby enabling a more nuanced understanding of gene expression regulation. Our approach facilitates a variety of biological analyses, including clustering by the kinetic rate, identifying genes with highly variable kinetics across cell types, and detecting RNA-binding proteins that influence splicing and degradation processes. Importantly, DeepKINET utilizes readily available scRNA-seq data, avoids the need for complex metabolic labeling, and paves the way for novel investigations of gene expression kinetics. Using this method, one can gain insights into the regulatory mechanisms of gene expression and uncover potential therapeutic targets for diseases in which splicing and degradation are dysregulated, such as cancer. These insights will be critical in elucidating variations in gene expression among cells and populations, bringing to light complex regulatory networks.

Despite its advantages, DeepKINET has several inherent limitations. It employs a unified model to estimate splicing and degradation rates, which can lead to correlation trends among these rates (Additional file 1: Fig. S1d, S2c, S5). Nonetheless, the fidelity of our estimates was supported by simulated and metabolic labeling data. In addition, the correlation between splicing and degradation rates estimated by DeepKINET was the lowest among the three methods, while cellDancer and DeepVelo exhibited high correlations. These high correlations suggest that cellDancer and DeepVelo may have limited ability to disentangle the effects of splicing and degradation. While kinetic rate estimation at the single-cell level improves the details of RNA velocity calculations [8], the simultaneous estimation of RNA velocity and kinetic rates presents a challenge, indicating the need for further methodological enhancements and additional constraints for improved accuracy in estimating kinetic rates.

It is worth noting that by extending DeepKINET, the assumption of fixed transcription rates for each gene can also be eliminated. However, this would increase the number of parameters, and further investigation is required to determine whether the estimation of the transcription, splicing, and degradation rates would all remain stable in such a setting. An existing method MultiVelo [46] uses multi-omics data (gene expression and chromatin accessibility) as input and assumes that transcription rates are determined based on chromatin accessibility near the gene. Considering the fact that transcription factors bind almost exclusively to open chromatin and provide dynamic regulation of transcription [47], MultiVelo’s modeling is more realistic than estimating transcription rates using only scRNA-seq data and may allow for more accurate estimation of transcription rates.

A notable challenge lies in the current limitations of RNA velocity analysis in distinguishing mRNA isoforms [48], with implications particularly relevant to diseases such as cancer, where alternative splicing is prevalent. Addressing this issue in future versions of DeepKINET could provide deeper biological insights and a more authentic portrayal of variations in mRNA splicing.

In summary, DeepKINET is a significant contributor to the field of single-cell biology, offering a novel analytical framework that not only advances the current understanding but also sets the stage for future innovations that will further elucidate the complexities of cellular kinetics.

Methods

In DeepKINET, the cell states and RNA velocity were first estimated, as in VICDYF [9], and the learned parameters of the encoders and decoders were fixed. Subsequently, decoders are created that take the cell states as the input and output the splicing and degradation rates at the single-cell level. These decoders are trained to better reconstruct unspliced mRNA amounts.

Derivation of single-cell splicing and degradation rates

The cell state and RNA velocity were estimated as described in the previous VICDYF method. The standard normal distribution is used as a prior for the low-dimensional latent cell state zn∈RD of cell n and the direction of small change dn∈RD on the low-dimensional latent cell space. D is the dimension of the latent cell space and the default value is 20.p(zn)=N(0,I),

p(dn)=N(0,ρI)

where ρ is a scaling factor, and I is the identity matrix. The direction of the small change dn needs to have a small scale with respect to zn; thus, we set ρ=0.01 to be the same as in VICDYF. Unspliced and spliced transcriptomes of a single cell are indicated by un∈Rg and sn∈Rg, where g is the number of genes. Poisson or negative binomial distributions were assumed for the distributions of un and the distribution of sn given zn. A Poisson distribution was assumed for all analyses in this research.sn~=lsnλθ(zn),

p(sn|zn)=Poisson(sn~)

where lsn∈Rg is the mean of spliced counts across all genes in the single cell, and λθ(zn)∈Rg is the decoding neural network of the latent cell states with 100 hidden units, one hidden layer, and layer normalization. sn~ is the reconstructed spliced mRNA counts. We derived the approximate time change in the mean parameter of the spliced transcriptome by decoding a small change in the latent cell state. In VICDYF, only s is used as input for the VAE to quantify the uncertainty of u given s. However, in DeepKINET, both u and s are used as inputs because we do not focus on the uncertainty of u. Moreover, to determine the small change in s, we differentiate the decoder transformation from z to s~ by z using a functorch instead of using the central difference approximation in VICDYF.1 vn=∂λθ(zn)∂zndn.

Here, we assumed that the mean parameter of the abundance of spliced and unspliced transcriptomes was represented by the differential equation of splicing kinetics as an RNA velocity estimation.2 un~≈lunvn+dtγsn~dtβ.

where β∈Rg is a vector of gene-specific splicing rates of unspliced transcripts and γ∈Rg is a vector of gene-specific degradation rates of spliced transcripts. Here, β and γ are the same value for each cell. lun∈Rg is the mean of unspliced counts across all genes in the single cell. un~ is the reconstructed unspliced mRNA counts. By combining (1) and (2), we can approximate the mean parameter of the abundance of unspliced transcripts as follows:un~≈lun∂λθ(zn)∂zndn+dtγsn~dtβ

We assumed that the abundance of unspliced transcriptomes u has a Poisson distribution, as follows: where dt is the small interval and is set to 1.p(un|zn,dn)=Poisson(un~)

We assume the variational posterior distribution of zn is a Gaussian distribution that depends on the raw counts of spliced and unspliced mRNA and the variational posterior distribution of dn is a Gaussian distribution that depends on zn.q(zn|sn,un)=N(μϕ(sn,un),diag(σϕ(sn,un)))

q(dn|zn)=N(μϕ′(zn),diag(σϕ′(zn)))

where μϕ() and σϕ() are the encoding neural networks with 100 hidden units, two hidden layers, and layer normalization [49]. μϕ′(sn,un) and σϕ′(sn,un) are the encoding neural network with 100 hidden units, one hidden layers, and layer normalization.

The generative model and variational posterior distribution were optimized by minimizing the following loss function: Minimizing this loss function is equivalent to maximizing the variational lower bound (ELBO) of transcriptome distribution. This minimization allowed us to learn about the variational autoencoder of the spliced transcriptome, RNA velocity, and the reconstruction of the unspliced transcriptome.L(θ,ϕ)=-Eq(zn,dn|sn,un)logp(sn,un,zn,dn)q(zn,dn|sn,un)≥-logp(sn|zn′)-logp(un|zn′,dn′)+DKL(q(zn|sn,un)‖p(zn))+DKL(q(dn|zn′)‖p(dn))

where zn′ and dn′ are derived through reparametrized sampling from q(zn|sn,un) and q(dn|zn′), and Ep(x)[f(x)] represents the expectation of f(x) given x∼p(x). To minimize the loss function, the Adam W optimizer was used with a learning rate of 0.001 and a mini-batch size of 100. Learning ended when the average loss of the 10 epochs was not been updated for 10 epochs.

After learning the VAE and RNA velocity, and reconstructing the unspliced transcripts as described above, all encoder and decoder parameters were fixed. Next, we create decoders that take latent variables as inputs and output splicing rate βn and degradation rate γn at the single cell level. When reconstructing unspliced transcripts, they were substituted for the previous splicing and degradation rates. By estimating the splicing and degradation rates as cell-state-dependent values, the rates for cells with similar cell states will be similar, weakening the indeterminacy of the solution.un~≈lun∂λθ(zn)∂zndn+dtγnsn~dtβn

The same loss function described above was used to learn the splicing and degradation rates at the single-cell level.

Conditional model that handles multiple samples

To address data with multiple samples, we extended DeepKINET with a conditional VAE framework [40]. The prior distribution of zn and dn is assumed to be the same distribution as in the previous model. We assume a variational posterior distribution of cell state zn with raw mRNA counts sn, un and batches bn∈{0,1}B as inputs. B is the total number of the experimental batches and bn,k=1 denote cell n belongs to experimental batch k.q(zn|sn,un,bn)=N(μϕ(sn,un,bn),diag(σϕ(sn,un,bn)))

We then reconstruct spliced mRNA as follows.sn~=lsnλθ(zn,bn)

The rest of the model is the same as the model described in the previous section. After training the VAE and RNA velocity, we fix encoders and decoders and use the latent state and batch information of each cell as input to estimate the splicing rate and degradation rate of each cell.

Creating simulated datasets

We used SERGIO (version 1.0.0) to generate the scRNA-seq count data with varying splicing and degradation rates per cluster. We used the DS6 differentiation process and the gene network from SERGIO. The SERGIO source code was rewritten to allow the splicing and degradation rates to change on a cluster-by-cluster basis. The base rate for each cell cluster was set by multiplying the SERGIO default splicing and degradation rate values by values sampled from a uniform distribution of 0.5 to 1.5. The base kinetic rates were then multiplied by values sampled from a uniform distribution of 0.75 to 1.25 for each cluster to establish different rates for each cluster. Each cluster contained 100 cells. In experiments with varying dropout rates, the dropout indicator dynamics function was used. Twenty dropout rate conditions were set with shape=1 and five increments from percentile=0 to percentile=95. For the experiments in which the number of cells was varied, 13 conditions were set for the number of cells using a default dropout rate of shape=1 and percentile=65. Ten datasets were created for each condition using different splicing and degradation rates. DeepKINET and cellDancer were used once for each dataset.

Validation using metabolic labeling experimental dataset

Using the scEU-seq cell cycle dataset, we determined the splicing and degradation rates for each cluster using Dynamo [13] and compared the estimates with those from DeepKINET and cellDancer. We split the cell cycle into PULSE and CHASE data and performed default gene filtering using Dynamo to extract 1000 genes. We divided each dataset into cell clusters based on the cell cycle position, with each cluster containing the same number of cells. We then modified the dynamo.tl.recipe_kin_data and dynamo.tl.recipe_deg_data functions to calculate the kinetic rates for each cluster. Using other parameters and following the default values of Dynamo, we derived the splicing and degradation rates for each cluster. We then applied DeepKINET, cellDancer, and DeepVelo to the PULSE and CHASE data 100 times each and derived the correlation coefficients each time. In cellDancer, the seed value used for training was fixed in the source code, so the estimation was performed without setting this seed value.

For the scNT-seq hematopoietic dataset, we estimated the kinetic rates in two separate batches, each containing cells collected at different time points, as in the Dynamo tutorial. We filtered out genes exhibiting a high correlation (> 0.7, 75 genes) between the moments of unspliced and spliced mRNA. We then compared the ratio of degradation rates between the two batches between the estimated values of Dynamo and the estimated values of DeepKINET, cellDancer, and DeepVelo. In DeepKINET, we used a conditional VAE framework to address batch effects by using the time batches as batch labels.

Clustering by kinetic rates

The splicing and degradation rates of each cell were estimated using DeepKINET and Z-transformation. Principal component analysis was then performed using the rates. Leiden clustering was performed on the principal components to cluster the genes (Fig. 3c).

Functional enrichment analysis

We performed gene clustering using kinetic rates on the cell cycle PULSE data. GO analysis was performed on each gene cluster (Fig. 3d). We used g:Profiler [50] for the analysis. When conducting GO term analysis, we used all genes used to estimate splicing and degradation rates as the background.

Enrichment test of RNA-binding protein targets

Using the forebrain dataset, we performed gene clustering based on the kinetic rates. We examined whether the genes in each cluster were enriched for RNA-binding protein targets (Fig. 4b). We selected RNA-binding proteins that were included in the 1000 highly variable genes selected by preprocessing, for which eCLIP data were available in the CLIPdb [51]. Genes with at least one binding site in the eCLIP data were considered as targets. After performing Fisher’s exact test, we used the Benjamini-Hochberg method for multiple testing correction.

Analysis of the relationship between expression of RNA-binding proteins and kinetic rates of their targets

As a preprocessing step, we used scvelo.pp.filter_and_normalize() with min_shared_counts = 20 for the forebrain dataset and min_shared_counts = 100 for the breast cancer dataset. To ensure accuracy, we estimated the kinetic rates of genes with high variability. When all the remaining RNA-binding proteins from the expression preprocessing were used in the analysis, the expression was averaged over the neighborhoods. For the forebrain dataset, we used n_neighbors=30. When analyzing only the RNA-binding proteins in the highly variable genes, we used the expression reconstructed from the latent variables. The top 1000 genes in the forebrain dataset and the top 2000 genes in the breast cancer dataset were used as highly variable genes. When comparing the expression of a specific RNA-binding protein to its target or non-target kinetic rates, we used a t-test to determine any significant difference in the correlation coefficients between targets and non-targets.

Preparation of breast cancer data

We downloaded the FASTQ files from the public data of Liu et al. We then created BAM files using Cell Ranger [18]. Next, Velocyto [5] was used to create count matrices for unspliced and spliced mRNA. We used EPCAM and KRT19 as markers of epithelial cells, following the method described by Liu et al. We used inferCNVpy to extract the cancer cells. Among the seven patients, cells from patient 5 were selected and used for further analysis because the other patients contained few metastatic cells or, conversely, too many metastatic cells or a low number of breast cancer epithelial cells. Because tumor epithelial cells are highly heterogeneous from patient to patient [52], we did not perform an integrated analysis. Cells with at least 100 expressed genes and at least 500 unique molecular identifier counts were used. As a preprocessing step, we used scvelo.pp.filter_and_normalize() with min_shared_counts = 100 and n_top_genes = 2000 to extract genes with high expression variability.

Preparation of bone marrow mononuclear cell data

We downloaded the FASTQ files from the public data of Adema et al. [41]. We then created BAM files using Cell Ranger [18]. Next, we used Velocyto [5] to create count matrices for unspliced and spliced mRNA. As in the original paper, cells with at least 100 expressed genes and at least 500 unique molecular identifier counts were used. We annotated cells using CellTypist [53]. Then, we extracted erythroid lineage cells. As a preprocessing step, we used scvelo.pp.filter_and_normalize() with min_shared_counts = 20 and n_top_genes = 1000 to extract genes with high expression variability.

Identification of targets differentially regulated by different cell populations

We performed the following linear regression using the expression levels of RNA-binding proteins and the kinetic rates of their targets. We then examined whether the slope of the regression line differed significantly among the cell populations.y=β0+β1x1+β2x0x1+β3x0+ϵ

where x0 is the label of the cell population, 0 for primary cells and 1 for metastatic cells, x1 is the expression of a RNA-binding protein, and β0 to β3 are the regression coefficients. We set β2=0 as the null hypothesis and used statsmodels.regression.linear_model.OLS() to perform regression and testing. We corrected the p-values using the Benjamini-Hochberg method.

Two-dimensional embedding of velocity

We projected the transitions in the latent space onto two-dimensional coordinates following the method described by Bergen et al. [6]. We used zj-zi as the change in the latent space of cell i to cell j and di as the velocity in the latent space of cell i. We computed a neighborhood graph, calculated the transition probabilities, and projected them onto two-dimensional coordinates using Scvelo’s functions.

Supplementary Information

Additional file 1. Supplementary Figures S1-S5

Additional file 2. Review history

Acknowledgements

Supercomputing resources were provided by the SHIROKANE supercomputer at the Human Genome Center of the University of Tokyo, the TSUBAME3.0 supercomputer at the Tokyo Institute of Technology, and the AI Bridging Cloud Infrastructure (ABCI) at the National Institute of Advanced Industrial Science and Technology (AIST).

Review history

The review history is available as Additional file 2.

Peer review information

Kevin Pang and Veronique van den Berghe were the primary editors of this article and managed its editorial process and peer review in collaboration with the rest of the editorial team.

Authors’ contributions

Y.K. conceived the concept of the method. K.A. conceived the idea of validation through simulation. C.M. designed the source code for this method, conducted experiments to verify its validity, designed the analysis using this method, and performed the analysis under the supervision of Y.K. and T.S. S.N. and S.H. made minor modifications to the theory of this method. All the authors have read and approved the final version of the manuscript.

Funding

This research was funded by multiple sources. The Grant-in-Aid for Transformative Research Areas (platforms for Advanced Technologies and Research Resources) (grant no. 22H04925) and Grant-in-Aid for Transformative Research Areas (A) (grant no. 23H04938) were provided by the Japan Society for the Promotion of Science (JSPS). Additional support was received from the Project for P-PROMOTE (grant nos. JP22ama221215 and JP22ama221501), Brain/MINDS Health and Diseases (grant no. JP22wm0425007), the Interdisciplinary Cutting-edge Research (grant no. JP23wm0325068), and the Advanced Genome Research and Bioinformatics Study to Facilitate Medical Innovation (GRIFIN) (grant no. JP23tm0424226) from the Japan Agency for Medical Research and Development (AMED). The Moonshot R&D program (grant no. JPMJMS2025) also contributed, through the Japan Science and Technology Agency (JST). Further support came from the Medical Research Center Initiative for High Depth Omics and Multilayered Stress Diseases at Tokyo Medical and Dental University.

Availability of data and materials

The DeepKINET implementation and the modified SERGIO code used to create the simulation data are available at https://github.com/3254c/DeepKINET [54], which has also been deposited via Zenodo (https://zenodo.org/doi/10.5281/zenodo.13054695) [55]. The DeepKINET package are released with MIT license.

The pancreas dataset is available at scvelo.datasets.pancreas() or the original work [56] under the Gene Expression Omnibus (GEO) accession number GSE132188 [57]. The cell cycle dataset is available at dynamo.dataset.cellcycle() or the original work [3] under the GEO accession number GSE128365 [58]. The forebrain dataset is available at scvelo.datasets.forebrain() or the original work [5] under the Sequence Read Archive accession code SRP129388 [59]. The breast cancer data, including the raw data, can be obtained from the original work [19] under the GEO accession number GSE167036 [60].

Declarations

Ethics approval and consent to participate

Not applicable.

Consent for publication

All authors have reviewed and consented to the publication of this manuscript.

Competing interests

The authors declare no competing interests.

Publisher's Note

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

1. Bradley RK Anczuków O RNA splicing dysregulation and the hallmarks of cancer Nat Rev Cancer. 2023 23 3 135 155 10.1038/s41568-022-00541-7 36627445
Bradley RK, Anczuków O. RNA splicing dysregulation and the hallmarks of cancer. Nat Rev Cancer. 2023;23(3):135–55.36627445 10.1038/s41568-022-00541-7
2. Fang Z Mei W Qu C Lu J Shang L Cao F Li F Role of m6A writers, erasers and readers in cancer Exp Hematol Oncol. 2022 11 1 45 10.1186/s40164-022-00298-7 35945641
Fang Z, Mei W, Qu C, Lu J, Shang L, Cao F, Li F. Role of m6A writers, erasers and readers in cancer. Exp Hematol Oncol. 2022;11(1):45.35945641 10.1186/s40164-022-00298-7
3. Battich N Beumer J Barbanson BD Krenning L Baron CS Tanenbaum ME Clevers H Oudenaarden AV Sequencing metabolically labeled transcripts in single cells reveals mRNA turnover strategies Science. 2020 367 6482 1151 1156 10.1126/science.aax3072 32139547
Battich N, Beumer J, Barbanson BD, Krenning L, Baron CS, Tanenbaum ME, Clevers H, Oudenaarden AV. Sequencing metabolically labeled transcripts in single cells reveals mRNA turnover strategies. Science. 2020;367(6482):1151–6.32139547 10.1126/science.aax3072
4. Qiu Q Hu P Qin X Govek KW Cámara PG Wu H Massively parallel and time-resolved RNA sequencing in single cells with scNT-seq Nat Methods. 2020 17 10 991 1001 10.1038/s41592-020-0935-4 32868927
Qiu Q, Hu P, Qin X, Govek KW, Cámara PG, Wu H. Massively parallel and time-resolved RNA sequencing in single cells with scNT-seq. Nat Methods. 2020;17(10):991–1001.32868927 10.1038/s41592-020-0935-4
5. La Manno G Soldatov R Zeisel A Braun E Hochgerner H Petukhov V Lidschreiber K RNA velocity of single cells Nature. 2018 560 7719 494 498 10.1038/s41586-018-0414-6 30089906
La Manno G, Soldatov R, Zeisel A, Braun E, Hochgerner H, Petukhov V, Lidschreiber K, et al. RNA velocity of single cells. Nature. 2018;560(7719):494–8.30089906 10.1038/s41586-018-0414-6
6. Bergen V Lange M Peidli S Wolf FA Theis FJ Generalizing RNA velocity to transient cell states through dynamical modeling Nat Biotechnol. 2020 38 12 1408 1414 10.1038/s41587-020-0591-3 32747759
Bergen V, Lange M, Peidli S, Wolf FA, Theis FJ. Generalizing RNA velocity to transient cell states through dynamical modeling. Nat Biotechnol. 2020;38(12):1408–14.32747759 10.1038/s41587-020-0591-3
7. Gayoso A Weiler P Lotfollahi M Klein D Hong J Streets A Theis FJ Yosef N Deep generative modeling of transcriptional dynamics for RNA velocity analysis in single cells Nat Methods. 2024 21 1 50 59 10.1038/s41592-023-01994-w 37735568
Gayoso A, Weiler P, Lotfollahi M, Klein D, Hong J, Streets A, Theis FJ, Yosef N. Deep generative modeling of transcriptional dynamics for RNA velocity analysis in single cells. Nat Methods. 2024;21(1):50–9.37735568 10.1038/s41592-023-01994-w
8. Li S, Zhang P, Chen W, Ye L, Brannan KW, Le N, Abe J, et al. A relay velocity model infers cell-dependent RNA velocity. Nat Biotechnol. 2023. 10.1038/s41587-023-01728-5.
9. Nagaharu K Kojima Y Hirose H Minoura K Hinohara K Minami H Kageyama Y A bifurcation concept for B-lymphoid/plasmacytoid dendritic cells with largely fluctuating transcriptome dynamics Cell Rep. 2022 40 9 111260 10.1016/j.celrep.2022.111260 36044861
Nagaharu K, Kojima Y, Hirose H, Minoura K, Hinohara K, Minami H, Kageyama Y, et al. A bifurcation concept for B-lymphoid/plasmacytoid dendritic cells with largely fluctuating transcriptome dynamics. Cell Rep. 2022;40(9):111260.36044861 10.1016/j.celrep.2022.111260
10. Wolf FA Hamey FK Plass M Solana J Dahlin JS Göttgens B Rajewsky N PAGA: graph abstraction reconciles clustering with trajectory inference through a topology preserving map of single cells Genome Biol. 2019 20 1 59 10.1186/s13059-019-1663-x 30890159
Wolf FA, Hamey FK, Plass M, Solana J, Dahlin JS, Göttgens B, Rajewsky N, et al. PAGA: graph abstraction reconciles clustering with trajectory inference through a topology preserving map of single cells. Genome Biol. 2019;20(1):59.30890159 10.1186/s13059-019-1663-x
11. Dibaeinia P Sinha S SERGIO: A Single-Cell Expression Simulator Guided by Gene Regulatory Networks Cell Syst. 2020 11 3 252 271.e11 10.1016/j.cels.2020.08.003 32871105
Dibaeinia P, Sinha S. SERGIO: a single-cell expression simulator guided by gene regulatory networks. Cell Syst. 2020;11(3):252-271.e11.32871105 10.1016/j.cels.2020.08.003
12. Cui H Maan H Vladoiu MC Zhang J Taylor MD Wang B DeepVelo: deep learning extends RNA velocity to multi-lineage systems with cell-specific kinetics Genome Biol. 2024 25 1 27 10.1186/s13059-023-03148-9 38243313
Cui H, Maan H, Vladoiu MC, Zhang J, Taylor MD, Wang B. DeepVelo: deep learning extends RNA velocity to multi-lineage systems with cell-specific kinetics. Genome Biol. 2024;25(1):27.38243313 10.1186/s13059-023-03148-9
13. Qiu X Zhang Y Martin-Rufino JD Weng C Hosseinzadeh S Yang D Pogson AN Mapping transcriptomic vector fields of single cells Cell. 2022 185 4 690 711.e45 10.1016/j.cell.2021.12.045 35108499
Qiu X, Zhang Y, Martin-Rufino JD, Weng C, Hosseinzadeh S, Yang D, Pogson AN, et al. Mapping transcriptomic vector fields of single cells. Cell. 2022;185(4):690-711.e45.35108499 10.1016/j.cell.2021.12.045
14. Tirosh I Izar B Prakadan SM Wadsworth MH 2nd Treacy D Trombetta JJ Rotem A Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq Science. 2016 352 6282 189 96 10.1126/science.aad0501 27124452
Tirosh I, Izar B, Prakadan SM, Wadsworth MH 2nd, Treacy D, Trombetta JJ, Rotem A, et al. Dissecting the multicellular ecosystem of metastatic melanoma by single-cell RNA-seq. Science. 2016;352(6282):189–96.27124452 10.1126/science.aad0501
15. Conboy JG. Developmental regulation of RNA processing by Rbfox proteins. Wiley Interdiscip Rev RNA. 2017;8(2). 10.1002/wrna.1398.
16. Pereira B Billaud M Almeida R RNA-Binding Proteins in Cancer: Old Players and New Actors Trends Cancer. 2017 3 7 506 528 10.1016/j.trecan.2017.05.003 28718405
Pereira B, Billaud M, Almeida R. RNA-binding proteins in cancer: old players and new actors. Trends Cancer. 2017;3(7):506–28.28718405 10.1016/j.trecan.2017.05.003
17. Qin H Ni H Liu Y Yuan Y Xi T Li X Zheng L RNA-binding proteins in tumor progression J Hematol Oncol. 2020 13 1 90 10.1186/s13045-020-00927-w 32653017
Qin H, Ni H, Liu Y, Yuan Y, Xi T, Li X, Zheng L. RNA-binding proteins in tumor progression. J Hematol Oncol. 2020;13(1):90.32653017 10.1186/s13045-020-00927-w
18. Zheng G Terry J Belgrader P Massively parallel digital transcriptional profiling of single cells Nat Commun. 2017 8 14049 10.1038/ncomms14049 28091601
Zheng G, Terry J, Belgrader P, et al. Massively parallel digital transcriptional profiling of single cells. Nat Commun. 2017;8:14049.28091601 10.1038/ncomms14049
19. Liu T Liu C Yan M Zhang J Xiao M Li Z Wei X Zhang H Single cell profiling of primary and paired metastatic lymph node tumors in breast cancer patients Nat Commun. 2022 13 1 6823 10.1038/s41467-022-34581-2 36357424
Liu T, Liu C, Yan M, Zhang J, Xiao M, Li Z, Wei X, Zhang H. Single cell profiling of primary and paired metastatic lymph node tumors in breast cancer patients. Nat Commun. 2022;13(1):6823.36357424 10.1038/s41467-022-34581-2
20. Xiao JF, Kua LF, Ding LW, Sun QY, Myint KN, Chia XR, Venkatachalam N, et al. KDM6A depletion in breast epithelial cells leads to reduced sensitivity to anticancer agents and increased TGF activity. Mol Cancer Res. 2022;20(4):637–49.
21. Fowler AM Salem K DeGrave M Ong IM Rassman S Powers GL Kumar M Progesterone Receptor Gene Variants in Metastatic Estrogen Receptor Positive Breast Cancer Horm Cancer. 2020 11 2 63 75 10.1007/s12672-020-00377-3 31942683
Fowler AM, Salem K, DeGrave M, Ong IM, Rassman S, Powers GL, Kumar M, et al. Progesterone receptor gene variants in metastatic estrogen receptor positive breast cancer. Horm Cancer. 2020;11(2):63–75.31942683 10.1007/s12672-020-00377-3
22. Fusco N Malapelle U Fassan M Marchiò C Buglioni S Zupo S Criscitiello C PIK3CA Mutations as a Molecular Target for Hormone Receptor-Positive, HER2-Negative Metastatic Breast Cancer Front Oncol. 2021 11 644737 10.3389/fonc.2021.644737 33842357
Fusco N, Malapelle U, Fassan M, Marchiò C, Buglioni S, Zupo S, Criscitiello C, et al. PIK3CA mutations as a molecular target for hormone receptor-positive, HER2-negative metastatic breast cancer. Front Oncol. 2021;11:644737.33842357 10.3389/fonc.2021.644737
23. Yi Y, Chen D, Ao J, Zhang W, Yi J, Ren X, Fei J, et al. Transcriptional suppression of AMPK1 promotes breast cancer metastasis upon oncogene activation. Proc Natl Acad Sci U S A. 2020;117(14):8013–21.
24. Zhang J Zhang J Xu S Zhang X Wang P Wu H Xia B Hypoxia-Induced TPM2 Methylation is Associated with Chemoresistance and Poor Prognosis in Breast Cancer Cell Physiol Biochem. 2018 45 2 692 705 10.1159/000487162 29414807
Zhang J, Zhang J, Xu S, Zhang X, Wang P, Wu H, Xia B, et al. Hypoxia-induced TPM2 methylation is associated with chemoresistance and poor prognosis in breast cancer. Cell Physiol Biochem. 2018;45(2):692–705.29414807 10.1159/000487162
25. Gatti V Bongiorno-Borbone L Fierro C Annicchiarico-Petruzzelli M Melino G Peschiaroli A p63 at the Crossroads between Stemness and Metastasis in Breast Cancer Int J Mol Sci. 2019 20 11 2683 10.3390/ijms20112683 31159154
Gatti V, Bongiorno-Borbone L, Fierro C, Annicchiarico-Petruzzelli M, Melino G, Peschiaroli A. p63 at the Crossroads between Stemness and Metastasis in Breast Cancer. Int J Mol Sci. 2019;20(11):2683.31159154 10.3390/ijms20112683
26. Guan T Yang X Liang H Chen J Chen Y Zhu Y Liu T Deubiquitinating enzyme USP9X regulates metastasis and chemoresistance in triple-negative breast cancer by stabilizing Snail1 J Cell Physiol. 2022 237 7 2992 3000 10.1002/jcp.30763 35506169
Guan T, Yang X, Liang H, Chen J, Chen Y, Zhu Y, Liu T. Deubiquitinating enzyme USP9X regulates metastasis and chemoresistance in triple-negative breast cancer by stabilizing Snail1. J Cell Physiol. 2022;237(7):2992–3000.35506169 10.1002/jcp.30763
27. Peeney D Jensen SM Castro NP Kumar S Noonan S Handler C Kuznetsov A TIMP-2 suppresses tumor growth and metastasis in murine model of triple-negative breast cancer Carcinogenesis. 2020 41 3 313 325 10.1093/carcin/bgz172 31621840
Peeney D, Jensen SM, Castro NP, Kumar S, Noonan S, Handler C, Kuznetsov A, et al. TIMP-2 suppresses tumor growth and metastasis in murine model of triple-negative breast cancer. Carcinogenesis. 2020;41(3):313–25.31621840 10.1093/carcin/bgz172
28. Kim Y Wom M Chari T Lee S Park C Son C Kim KK RBM47-regulated alternative splicing of TJP1 promotes actin stress fiber assembly during epithelial-to-mesenchymal transition Oncogene. 2019 38 38 6521 6536 10.1038/s41388-019-0892-5 31358901
Kim Y, Wom M, Chari T, Lee S, Park C, Son C, Kim KK. RBM47-regulated alternative splicing of TJP1 promotes actin stress fiber assembly during epithelial-to-mesenchymal transition. Oncogene. 2019;38(38):6521–36.31358901 10.1038/s41388-019-0892-5
29. Vanharanta S Marney CB Shu W Valiente M Zou Y Mele A Darnell RB Massagué J Loss of the multifunctional RNA-binding protein RBM47 as a source of selectable metastatic traits in breast cancer Elife. 2014 3 e02734 10.7554/eLife.02734 24898756
Vanharanta S, Marney CB, Shu W, Valiente M, Zou Y, Mele A, Darnell RB, Massagué J. Loss of the multifunctional RNA-binding protein RBM47 as a source of selectable metastatic traits in breast cancer. Elife. 2014;3:e02734.24898756 10.7554/eLife.02734
30. Guo T You K Chen X Sun Y Wu Y Wu P Jiang Y RBM47 inhibits hepatocellular carcinoma progression by targeting UPF1 as a DNA/RNA regulator Cell Death Discov. 2022 8 1 320 10.1038/s41420-022-01112-3 35831298
Guo T, You K, Chen X, Sun Y, Wu Y, Wu P, Jiang Y. RBM47 inhibits hepatocellular carcinoma progression by targeting UPF1 as a DNA/RNA regulator. Cell Death Discov. 2022;8(1):320.35831298 10.1038/s41420-022-01112-3
31. Xiao Y Cong M Li J He D Wu Q Tian P Wang Y Cathepsin C promotes breast cancer lung metastasis by modulating neutrophil infiltration and neutrophil extracellular trap formation Cancer Cell. 2021 39 3 423 437.e7 10.1016/j.ccell.2020.12.012 33450198
Xiao Y, Cong M, Li J, He D, Wu Q, Tian P, Wang Y, et al. Cathepsin C promotes breast cancer lung metastasis by modulating neutrophil infiltration and neutrophil extracellular trap formation. Cancer Cell. 2021;39(3):423-437.e7.33450198 10.1016/j.ccell.2020.12.012
32. Jin L Zheng D Bhandari A Chen D Xia E Guan Y Wen J Wang O PSD3 is an oncogene that promotes proliferation, migration, invasion, and G1/S transition while inhibits apoptotic in papillary thyroid cancer J Cancer 2021 12 18 5413 5422 10.7150/jca.60885 34405004
Jin L, Zheng D, Bhandari A, Chen D, Xia E, Guan Y, Wen J, Wang O. PSD3 is an oncogene that promotes proliferation, migration, invasion, and G1/S transition while inhibits apoptotic in papillary thyroid cancer. J Cancer. 2021;12(18):5413–22.34405004 10.7150/jca.60885
33. Jovanović B Beeler JS Pickup MW Chytil A Gorska AE Ashby WJ Lehmann BD Transforming growth factor beta receptor type III is a tumor promoter in mesenchymal-stem like triple negative breast cancer Breast Cancer Res. 2014 16 4 R69 10.1186/bcr3684 24985072
Jovanović B, Beeler JS, Pickup MW, Chytil A, Gorska AE, Ashby WJ, Lehmann BD, et al. Transforming growth factor beta receptor type III is a tumor promoter in mesenchymal-stem like triple negative breast cancer. Breast Cancer Res. 2014;16(4):R69.24985072 10.1186/bcr3684
34. Liu Y Tang W Yao F USP53 Exerts Tumor-Promoting Effects in Triple-Negative Breast Cancer by Deubiquitinating CRKL Cancers (Basel). 2023 15 20 5033 10.3390/cancers15205033 37894400
Liu Y, Tang W, Yao F. USP53 exerts tumor-promoting effects in triple-negative breast cancer by deubiquitinating CRKL. Cancers (Basel). 2023;15(20):5033.37894400 10.3390/cancers15205033
35. Xie W Chen C Han Z Huang J Liu X Chen H Zhang T CD2AP inhibits metastasis in gastric cancer by promoting cellular adhesion and cytoskeleton assembly Mol Carcinog. 2020 59 4 339 352 10.1002/mc.23158 31989722
Xie W, Chen C, Han Z, Huang J, Liu X, Chen H, Zhang T, et al. CD2AP inhibits metastasis in gastric cancer by promoting cellular adhesion and cytoskeleton assembly. Mol Carcinog. 2020;59(4):339–52.31989722 10.1002/mc.23158
36. Ma WR Xu P Liu ZJ Zhou J Gu LK Zhang J Deng DJ Impact of GFRA1 gene reactivation by DNA demethylation on prognosis of patients with metastatic colon cancer World J Gastroenterol. 2020 26 2 184 198 10.3748/wjg.v26.i2.184 31988584
Ma WR, Xu P, Liu ZJ, Zhou J, Gu LK, Zhang J, Deng DJ. Impact of GFRA1 gene reactivation by DNA demethylation on prognosis of patients with metastatic colon cancer. World J Gastroenterol. 2020;26(2):184–98.31988584 10.3748/wjg.v26.i2.184
37. Yuan J, Xing H, Li Y, Song Y, Zhang N, Xie M, Liu J, et al. EPB41 suppresses the Wnt/-catenin signaling in non-small cell lung cancer by sponging ALDOC. Cancer Lett. 2021;499:255–64.
38. Ochi T Fujiwara T Ono K Suzuki C Nikaido M Inoue D Kato H Exploring the mechanistic link between SF3B1 mutation and ring sideroblast formation in myelodysplastic syndrome Sci Rep. 2022 12 1 14562 10.1038/s41598-022-18921-2 36028755
Ochi T, Fujiwara T, Ono K, Suzuki C, Nikaido M, Inoue D, Kato H, et al. Exploring the mechanistic link between SF3B1 mutation and ring sideroblast formation in myelodysplastic syndrome. Sci Rep. 2022;12(1):14562.36028755 10.1038/s41598-022-18921-2
39. Shiozawa Y Malcovati L Gallì A Sato-Otsubo A Kataoka K Sato Y Watatani Y Aberrant splicing and defective mRNA production induced by somatic spliceosome mutations in myelodysplasia Nat Commun. 2018 9 1 3649 10.1038/s41467-018-06063-x 30194306
Shiozawa Y, Malcovati L, Gallì A, Sato-Otsubo A, Kataoka K, Sato Y, Watatani Y, et al. Aberrant splicing and defective mRNA production induced by somatic spliceosome mutations in myelodysplasia. Nat Commun. 2018;9(1):3649.30194306 10.1038/s41467-018-06063-x
40. Kingma DP, Mohamed S, Rezende DJ, Welling M. Semi-supervised learning with deep generative models. Neural Inf Process Syst. 2014;27:3581–89.
41. Adema V Ma F Kanagal-Shamanna R Thongon N Montalban-Bravo G Yang H Peslak SA Targeting the EIF2AK1 Signaling Pathway Rescues Red Blood Cell Production in SF3B1-Mutant Myelodysplastic Syndromes With Ringed Sideroblasts Blood Cancer Discov. 2022 3 6 554 567 10.1158/2643-3230.BCD-21-0220 35926182
Adema V, Ma F, Kanagal-Shamanna R, Thongon N, Montalban-Bravo G, Yang H, Peslak SA, et al. Targeting the EIF2AK1 signaling pathway rescues red blood cell production in SF3B1-mutant myelodysplastic syndromes with ringed sideroblasts. Blood Cancer Discov. 2022;3(6):554–67.35926182 10.1158/2643-3230.BCD-21-0220
42. Keszei M, Kritikou JS, Sandfort D, He M, Oliveira MMS, Wurzer H, Kuiper RV, Westerberg LS. Wiskott-Aldrich syndrome gene mutations modulate cancer susceptibility in the p53 murine model. Oncoimmunology. 2018;7(9):e1468954.
43. Biber G Ben-Shmuel A Noy E Joseph N Puthenveetil A Reiss N Levy O Targeting the actin nucleation promoting factor WASp provides a therapeutic approach for hematopoietic malignancies Nat Commun. 2021 12 1 5581 10.1038/s41467-021-25842-7 34552085
Biber G, Ben-Shmuel A, Noy E, Joseph N, Puthenveetil A, Reiss N, Levy O, et al. Targeting the actin nucleation promoting factor WASp provides a therapeutic approach for hematopoietic malignancies. Nat Commun. 2021;12(1):5581.34552085 10.1038/s41467-021-25842-7
44. Miao G Zhuo D Han X Yao W Liu C Liu H Cao H From degenerative disease to malignant tumors: Insight to the function of ApoE Biomed Pharmacother. 2023 158 114127 10.1016/j.biopha.2022.114127 36516696
Miao G, Zhuo D, Han X, Yao W, Liu C, Liu H, Cao H, et al. From degenerative disease to malignant tumors: insight to the function of ApoE. Biomed Pharmacother. 2023;158:114127.36516696 10.1016/j.biopha.2022.114127
45. Xue Y Lu F Chang Z Li J Gao Y Zhou J Luo Y Intermittent dietary methionine deprivation facilitates tumoral ferroptosis and synergizes with checkpoint blockade Nat Commun. 2023 14 1 4758 10.1038/s41467-023-40518-0 37553341
Xue Y, Lu F, Chang Z, Li J, Gao Y, Zhou J, Luo Y, et al. Intermittent dietary methionine deprivation facilitates tumoral ferroptosis and synergizes with checkpoint blockade. Nat Commun. 2023;14(1):4758.37553341 10.1038/s41467-023-40518-0
46. Li C Virgilio MC Collins KL Welch JD Multi-omic single-cell velocity models epigenome-transcriptome interactions and improves cell fate prediction Nat Biotechnol. 2023 41 3 387 398 10.1038/s41587-022-01476-y 36229609
Li C, Virgilio MC, Collins KL, Welch JD. Multi-omic single-cell velocity models epigenome-transcriptome interactions and improves cell fate prediction. Nat Biotechnol. 2023;41(3):387–98.36229609 10.1038/s41587-022-01476-y
47. Klemm SL Shipony Z Greenleaf WJ Chromatin accessibility and the regulatory epigenome Nat Rev Genet. 2019 20 207 220 10.1038/s41576-018-0089-8 30675018
Klemm SL, Shipony Z, Greenleaf WJ. Chromatin accessibility and the regulatory epigenome. Nat Rev Genet. 2019;20:207–20.30675018 10.1038/s41576-018-0089-8
48. Gorin G Fang M Chari T Pachter L RNA velocity unraveled PLoS Comput Biol. 2022 18 9 e1010492 10.1371/journal.pcbi.1010492 36094956
Gorin G, Fang M, Chari T, Pachter L. RNA velocity unraveled. PLoS Comput Biol. 2022;18(9):e1010492.36094956 10.1371/journal.pcbi.1010492
49. Ba JL, Kiros JR, Hinton GE. Layer normalization. 2016. arXiv preprint arXiv:1607.06450. Available from: http://arxiv.org/abs/1607.06450.
50. Raudvere U, Kolberg L, Kuzmin I, Arak1 T, Adler P, Peterson H, Vilo J. g:Profiler: a web server for functional enrichment analysis and conversions of gene lists (2019 update). Nucleic Acids Res. 2019;47(Web Server issue):W191–8.
51. Yang YT Di C Hu B Zhou M Liu Y Song N Li Y Umetsu J Lu ZJ CLIPdb: a CLIP-seq database for protein-RNA interactions BMC Genomics. 2015 16 1 51 10.1186/s12864-015-1273-2 25652745
Yang YT, Di C, Hu B, Zhou M, Liu Y, Song N, Li Y, Umetsu J, Lu ZJ. CLIPdb: a CLIP-seq database for protein-RNA interactions. BMC Genomics. 2015;16(1):51.25652745 10.1186/s12864-015-1273-2
52. Turashvili G Brogi E Tumor Heterogeneity in Breast Cancer Front Med (Lausanne). 2017 4 227 10.3389/fmed.2017.00227 29276709
Turashvili G, Brogi E. Tumor heterogeneity in breast cancer. Front Med (Lausanne). 2017;4:227.29276709 10.3389/fmed.2017.00227
53. Domínguez Conde C Xu C Jarvis LB Rainbow DB Wells SB Gomes T Howlett SK Cross-tissue immune cell analysis reveals tissue-specific features in humans Science. 2022 376 6594 eabl5197 10.1126/science.abl5197 35549406
Domínguez Conde C, Xu C, Jarvis LB, Rainbow DB, Wells SB, Gomes T, Howlett SK, et al. Cross-tissue immune cell analysis reveals tissue-specific features in humans. Science. 2022;376(6594):eabl5197.35549406 10.1126/science.abl5197
54. Mizukoshi C, Kojima Y, Nomura S, Hayashi S, Abe K, Shimamura T. DeepKINET. Github. 2024. https://github.com/3254c/DeepKINET. Accessed 28 Jul 2024.
55. Mizukoshi C, Kojima Y, Nomura S, Hayashi S, Abe K, Shimamura T. DeepKINET v0.2.0. Zenodo. 2024. https://zenodo.org/doi/10.5281/zenodo.13054695. Accessed 28 Jul 2024.
56. Bastidas-Ponce A Tritschler S Dony L Scheibner K Tarquis-Medina M Salinno C Schirge S Comprehensive single cell mRNA profiling reveals a detailed roadmap for pancreatic endocrinogenesis Development. 2019 146 12 dev173849 10.1242/dev.173849 31160421
Bastidas-Ponce A, Tritschler S, Dony L, Scheibner K, Tarquis-Medina M, Salinno C, Schirge S, et al. Comprehensive single cell mRNA profiling reveals a detailed roadmap for pancreatic endocrinogenesis. Development. 2019;146(12):dev173849.31160421 10.1242/dev.173849
57. Bastidas-Ponce A, Tritschler S, Leander D, Scheibner K, Tarquis-Medina M, Salinno C, Schirge S, et al. Comprehensive single-cell mRNA profiling reveals a detailed roadmap for pancreatic endocrinogenesis. Datasets. Gene Expr Omnibus. 2019. http://identifiers.org/geo:GSE132188. Accessed 3 Oct 2023.
58. Battich N, Beumer J, de Barbanson B, Krenning L, Baron CS, Tanenbaum ME, Clevers H, van Oudenaarden A. Sequencing of metabolically labeled transcripts in single cells from RPE1-FUCCI cells and murine intestinal organoids. Datasets. Gene Expr Omnibus. 2020. http://identifiers.org/geo:GSE128365. Accessed 3 Oct 2023.
59. La Manno G, Soldatov R, Zeisel A, Braun E, Hochgerner H, Petukhov V, Lidschreiber K, et al. RNA velocity of single cells. Datasets. Seq Read Arch. 2018. https://www.ncbi.nlm.nih.gov/sra/SRP129388. Accessed 3 Oct 2023.
60. Liu T, Liu C, Zhang J, Wei X, Zhang H. Single cell profiling of primary and paired metastatic lymph node tumors in breast cancer patients(10x RNA and TCR). Datasets. Gene Expr Omnibus. 2022. http://identifiers.org/geo:GSE167036. Accessed 3 Oct 2023.
