scFPC-DE: Robust Differential Expression Analysis Along Single Cell Trajectories via Functional Principal Component Analysis
Lopez Candelaria, R. J.; Qian, Y. J.; Law, M.; Chen, F.; Zhang, Y.; Qiu, X.
Show abstract
MotivationIdentifying temporally differentially expressed genes (TDEGs) along pseudotime trajectories from single cell RNA sequencing (scRNA-seq) data helps characterize the cellular states that underlie the dynamic process of cellular development. However, existing tests based on generalized additive models (GAMs) suffer from increased false positive rates under zero inflation caused by high dropouts, a ubiquitous technical artifact of scRNA-seq data. Furthermore, by testing each gene independently, existing tests ignore the variance-covariance structure shared across genes along the trajectory, leading to suboptimal power and reduced interpretability. ResultsWe present scFPC-DE, a trajectory-based differential expression analysis (TDEA) method based on functional data analysis (FDA). It models the gene expression as a function of pseudotime in the L2 space and represents the covariance structure of these functions by eigenfunctions derived from functional principal component (FPC) analysis. This approach effectively captures informative gene expression patterns along the trajectory, while mitigating the influence of zero inflation in both simulation and real data analysis. In simulations, scFPC-DE exhibited superior control of type I error and achieved the highest ROC-AUC among competing methods. When applied to an scRNA-seq dataset of B cell subtypes, scFPC-DE uniquely identified TDEGs enriched for B cell differentiation pathways, outperforming existing methods in biological relevance. These results show that scFPC-DE effectively captures the shared gene expression variation and pseudo-temporal structure along the single cell trajectory for TDEG identification. AvailabilityR package and code vignettes are publicly available at https://github.com/LopezRicardo1/scFPCDE. Contactxing_qiu@urmc.rochester.edu; yun.zhang@nih.gov. Supplementary informationSupplementary data are available at Bioinformatics online
Matching journals
The top 3 journals account for 50% of the predicted probability mass.
Similar papers in this journal
- PseudotimeDE: inference of differential gene expression along cell pseudotime with well-calibrated p-values from single-cell RNA sequencing data 97%
- Heterogeneous pseudobulk simulation enables realistic benchmarking of cell-type deconvolution methods 97%
- scPDA: Denoising Protein Expression in Droplet-Based Single-Cell Data 96%
Similar papers in this journal
- SMNN: Batch Effect Correction for Single-cell RNA-seq data via Supervised Mutual Nearest Neighbor Detection 97%
- SSMD: A semi-supervised approach for a robust cell type identification and deconvolution of mouse transcriptomics data 95%
- A comprehensive comparison on cell type composition inference for spatial transcriptomics data 95%
Similar papers in this journal
- A Markov Random Field Model for Network-based Differential Expression Analysis of Single-cell RNA-seq Data 96%
- scTensor detects many-to-many cell-cell interactions from single cell RNA-sequencing data 96%
- eSVD-DE: Cohort-wide differential expression in single-cell RNA-seq data using exponential-family embeddings 96%
Similar papers in this journal
"Similar papers" are the closest papers from that journal in the model's embedding space. They show what the match is built on, but the ranking comes mostly from a classifier over the whole training set, not from these examples alone.