Abstract
Retinoblastoma is a childhood retinal tumour that is the most common primary malignant intraocular tumour. However, it has been challenging to identify the cell types associated with genetic complexity. Here, we performed single-cell RNA sequencing on 14,739 cells from two retinoblastoma samples to delineate the heterogeneity and the underlying mechanism of retinoblastoma progression. Using a multiresolution network-based analysis, we identified two major cell types in human retinoblastoma. Cell trajectory analysis yielded a total of 5 cell states organized into two main branches, and the cell cycle-associated cone precursors were the cells of origin of retinoblastoma that were required for initiating the differentiation and malignancy process of retinoblastoma. Tumour cells differentiation reprogramming trajectory analysis revealed that cell-type components of multiple tumour-related pathways and predominantly expressed UBE2C were associated with an activation state in the malignant progression of the tumour, providing a potential novel "switch gene" marker during early critical stages in human retinoblastoma development. Thus, our findings improve our current understanding of the mechanism of retinoblastoma progression and are potentially valuable in providing novel prognostic markers for retinoblastoma.