CellBRF: a feature selection method for single-cell clustering using cell balance and random forest

Motivation: Single-cell RNA sequencing (scRNA-seq) offers a powerful tool to dissect the complexity of biological tissues through cell sub-population identification in combination with clustering approaches. Feature selection is a critical step for improving the accuracy and interpretability of single-cell clustering. Existing feature selection methods underutilize the discriminatory potential of genes across distinct cell types. We hypothesize that incorporating such information could further boost the performance of single cell clustering.Results: We develop CellBRF, a feature selection method that considers genes' relevance to cell types for single-cell clustering. The key idea is to identify genes that are most important for discriminating cell types through random forests guided by predicted cell labels. Moreover, it proposes a class balancing strategy to mitigate the impact of unbalanced cell type distributions on feature importance evaluation. We benchmark CellBRF on 33 scRNA-seq datasets representing diverse biological scenarios and demonstrate that it substantially outperforms state-of-the-art feature selection methods in terms of clustering accuracy and cell neighborhood consistency. Furthermore, we demonstrate the outstanding performance of our selected features through three case studies on cell differentiation stage identification, non-malignant cell subtype identification, and rare cell identification. CellBRF provides a new and effective tool to boost single-cell clustering accuracy.Availability and implementationAll source codes of CellBRF are freely available at .

Automated summary

Single-cell RNA sequencing (scRNA-seq) is a powerful tool for dissecting biological tissues. This study presents CellBRF, a feature selection method that considers the relevance of genes to cell types, improving single-cell clustering accuracy. CellBRF outperforms existing methods in terms of clustering accuracy and consistency, and has been successfully applied to cell differentiation stage identification, non-malignant cell subtype identification, and rare cell identification.