Single-cell transcriptional profiling reveals the heterogeneity in embryonal rhabdomyosarcoma

Rhabdomyosarcoma is the most common soft tissue sarcoma in children, and embryonal rhabdomyosarcoma is the most typical type of rhabdomyosarcoma. The heterogeneity, etiology, and origin of embryonal rhabdomyosarcoma remain unknown. After obtaining the gene expression data of every cell in the tumor tissue by single-cell RNA sequencing, we used the Seurat package in R studio for quality control, analysis, and exploration of the data. All cells are divided into tumor cells and non-tumor cells, and we chose tumor cells by marker genes. Then, we repeated the process to cluster the tumor cells and divided the subgroups by their differentially expressed genes and gene ontology/Kyoto Encyclopedia of Genes and Genomes analysis. Additionally, Monocle 2 was used for pseudo-time analysis to obtain the evolution trajectory of cells in tumor tissues. Tumor cells were divided into 5 subgroups according to their functions, which were characterized by high proliferation, sensing and adaptation to oxygen availability, enhanced epigenetic modification, enhanced nucleoside phosphonic acid metabolism, and ossification. Evolution trajectory of cells in tumor tissues is obtained. We used pseudo-time analysis to distinguish between mesenchymal stem cells and fibroblasts, proved that embryonal rhabdomyosarcoma in the pelvic originated from skeletal muscle progenitor cells, showed the evolutionary trajectory of embryonal rhabdomyosarcoma, and improved the method of evaluating the degree of malignancy of embryonal rhabdomyosarcoma.