Abstract
Single-cell RNA sequencing presents the sophisticated delineation of cell transcriptomes in many cancer types and highlights the tumor heterogeneity at higher resolution, which provides a new chance to explore the molecular mechanism in a minority of cells. In this study, we utilized publicly available single-cell RNA-seq data to discover and comprehensively dissect rare genes existing in few glioblastoma (GBM) cells. Moreover, we designed a framework to systematically identify 51 rare protein-coding genes (PCGs) and 47 rare long non-coding RNAs (lncRNAs) in GBM. Patients with high expression levels of rare genes like CYB5R2 and TPPP3 had worse overall survival and disease-free survival, implying their potential implication in GBM progression and prognosis. We found that these rare genes tended to be specifically expressed in GBM cancer stem cells, which emphasized their ability to characterize stem-like cancer cells and implied their contribution to GBM growth. Furthermore, rare genes were enriched in a 17-cell subset, which was located in an individual branch of the pseudotime trajectory of cancer progression and exhibited high cell cycle activity and invasive potential. Our study captures the rare genes highly expressed in few cells, deepens our understanding of special states during GBM tumorigenesis and progression such as cancer stemness and invasion, and proposes potential targets for cancer therapy.