Single-cell RNA-seq (scRNA-seq) allows the analysis of gene expression in each cell, which enables the detection of highly variable genes (HVG) that contribute to cell-to-cell variation within a homogeneous cell population. HVG detection is necessary for clustering analysis to improve the clustering result. scRNA-seq includes some genes that are expressed with a certain probability in all cells which make the cells indistinguishable. These genes are referred to as background noise. To remove the background noise and select the informative genes for clustering analysis, in this paper, we propose an effective HVG detection method based on principal component analysis (PCA). The proposed method utilizes PCA to evaluate the genes (features) on the sample space. The distortion-free principal components are selected to calculate the distance from the origin to gene as the weight of each gene. The genes that have the greatest distances to the origin are selected for clustering analysis. Experimental results on both synthetic and gene expression datasets show that the proposed method not only removes the background noise to select the informative genes for clustering analysis, but also outperforms the existing HVG detection methods.
A Oshlack , M D Robinson , M D Young . From RNA-seq reads to differential expression results. Genome Biology, 2010, 11( 12): 220–
2
X Ye , W Zhang , T Sakurai . Adaptive unsupervised feature learning for gene signature identification in non-small-cell lung cancer. IEEE Access, 2020, 8 : 154354– 154362
3
F Ozsolak , P M Milos . RNA sequencing: advances, challenges and opportunities. Nature Reviews Genetics, 2011, 12( 2): 87– 98
4
A Wagner , A Regev , N Yosef . Revealing the vectors of cellular identity with single-cell genomics. Nature Biotechnology, 2016, 34( 11): 1145– 1160
5
V Y Kiselev , T S Andrews , M Hemberg . Challenges in unsupervised clustering of single-cell RNA-seq data. Nature Reviews Genetics, 2019, 20( 5): 273– 282
6
X Ye , W Zhang , Y Futamura , T Sakurai . Detecting interactive gene groups for single-cell RNA-Seq data based on co-expression network analysis and subgraph learning. Cells, 2020, 9( 9): 1938–
7
X Ye , T Sakurai . Robust similarity measure for spectral clustering based on shared neighbors. ETRI Journal, 2016, 38( 3): 540– 550
8
F Emmert-Streib , M Dehmer , B Haibe-Kains . Gene regulatory networks and their applications: understanding biological and medical problems in terms of networks. Frontiers in Cell and Developmental Biology, 2014, 2 : 38–
9
D Thompson , A Regev , S Roy . Comparative analysis of gene regulatory networks: from network reconstruction to evolution. Annual Review of Cell and Developmental Biology, 2015, 31 : 399– 428
10
X Ye , T Sakurai . Spectral clustering with adaptive similarity measure in Kernel space. Intelligent Data Analysis, 2018, 22( 4): 751– 765
11
S H Yip , P C Sham , J Wang . Evaluation of tools for highly variable gene discovery from single-cell RNA-seq data. Briefings in Bioinformatics, 2019, 20( 4): 1583– 1589
12
G Finak , A McDavid , M Yajima , J Deng , V Gersuk , A K Shalek , C K Slichter , H W Miller , M J McElrath , M Prlic , P S Linsley , R Gottardo . MAST: a flexible statistical framework for assessing transcriptional changes and characterizing heterogeneity in single-cell RNA sequencing data. Genome Biology, 2015, 16 : 278–
13
S H Yip , P Wang , J P A Kocher , P C Sham , J Wang . Linnorm: improved statistical analysis for single cell RNA-seq expression data. Nucleic Acids Research, 2017, 45( 22): e179–
14
C W Law , Y Chen , W Shi , G K Smyth . Voom: precision weights unlock linear model analysis tools for RNA-seq read counts. Genome Biology, 2014, 15( 2): R29–
15
C A Vallejos , J C Marioni , S Richardson . BASiCS: bayesian analysis of single-cell sequencing data. PLoS Computational Biology, 2015, 11( 6): e1004333–
16
A T L Lun , K Bach , J C Marioni . Pooling across cells to normalize single-cell RNA sequencing data with many zero counts. Genome Biology, 2016, 17 : 75–
17
A T L Lun , D J McCarthy , J C Marioni . A step-by-step workflow for low-level analysis of single-cell RNA-seq data with Bioconductor. F1000Research, 2016, 5 : 2122–
18
P Brennecke , S Anders , J K Kim , A A Kolodziejczyk , X W Zhang , V Proserpio , B Baying , V Benes , S A Teichmann , J C Marioni , M G Heisler . Accounting for technical noise in single-cell RNA-seq experiments. Nature Methods, 2013, 10( 11): 1093– 1095
19
H I H Chen , Y Jin , Y Huang , Y Chen . Detection of high variability in gene expression from single-cell RNA-seq profiling. BMC Genomics, 2016, 17( S7): 508–
20
I Guyon, J Weston, S Barnhill, V Vapnik. Gene selection for cancer classification using support vector machines. Machine Learning, 2002, 46(1−3): 389−422
21
R Díaz-Uriarte , Andrés S A de . Gene selection and classification of microarray data using random forest. BMC Bioinformatics, 2006, 7 : 3–
22
R Satija , J A Farrell , D Gennert , A F Schier , A Regev . Spatial reconstruction of single-cell gene expression data. Nature Biotechnology, 2015, 33( 5): 495– 502
23
T Stuart , A Butler , P Hoffman , C Hafemeister , E Papalexi , III W M Mauck , Y Hao , M Stoeckius , P Smibert , R Satija . Comprehensive integration of single-cell data. Cell, 2019, 177( 7): 1888– 1902.e21
24
C Mayer , C Hafemeister , R C Bandler , R Machold , R B Brito , X Jaglin , K Allaway , A Butler , G Fishell , R Satija . Developmental diversification of cortical inhibitory interneurons. Nature, 2018, 555( 7697): 457– 462
25
H Hotelling . Analysis of a complex of statistical variables into principal components. Journal of Educational Psychology, 1993, 24( 6): 417– 441
26
I T Jolliffe. Principal Component Analysis. Springer, 1986
27
K Pearson. On the criterion that a given system of deviations from the probable in the case of a correlated system of variables is such that it can be reasonably supposed to have arisen from random sampling. In: Kotz S, Johnson N L, eds. Breakthroughs in Statistics. New York: Springer, 1992
28
Heckert N A, Filliben J J. NIST/SEMATECH e-Handbook of statistical methods; Chapter 1: Exploratory Data Analysis. 2003
29
T M Gierahn , II M H Wadsworth , T K Hughes , B D Bryson , A Butler , R Satija , S Fortune , J C Love , A K Shalek . Seq-Well: portable, low-cost RNA sequencing of single cells at high throughput. Nature Methods, 2017, 14( 4): 395– 398
30
A H Liu , T J Nowakowski , A A Pollen , J H Lui , M A Horlbeck , F J Attenello , D He , J S Weissman , A R Kriegstein , A A Diaz , D A Lim . Single-cell analysis of long non-coding RNAs in the developing human neocortex. Genome Biology, 2016, 17 : 67–
31
A A Pollen , T J Nowakowski , J Shuga , X Wang , A A Leyrat . Low-coverage single-cell mRNA sequencing reveals cellular heterogeneity and activated signaling pathways in developing cerebral cortex. Nature Biotechnology, 2014, 32( 10): 1053– 1058
32
C Hafemeister , R Satija . Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regression. Genome Biology, 2019, 20( 1): 296–
33
W M Rand . Objective criteria for the evaluation of clustering methods. Journal of the American Statistical Association, 1971, 66( 336): 846– 850
34
L McInnes , J Healy , N Saul , L Großberger . UMAP: uniform manifold approximation and projection. The Journal of Open Source Software, 2018, 3( 29): 861–
