Background: Cancer stem cells (CSCs) are biologically characterized by self-renewal, multi-directional differentiation and infinite proliferation, inducing anti-tumor drug resistance and metastasis. In the present study, we attempted to depict the baseline landscape of CSC-mediated biological properties, knowing that it is vital for tumor evolution, anti-tumor drug selection and drug resistance against fatal malignancy. Methods: We performed single-cell RNA sequencing (scRNA-seq) analysis in 15208 cells from a pair of primary and metastatic sites of collecting duct renal cell carcinoma (CDRCC). Cell subpopulations were identified and characterized by t-SNE, RNA velocity, monocle and other computational methods. Statistical analysis of all single-cell sequencing data was performed in R and Python. Results: A CSC population of 1068 cells was identified and characterized, showing excellent differentiation and self-renewal properties. These CSCs positioned as a center of the differentiation process and transformed into CDRCC primary and metastatic cells in spatial and temporal order, and played a pivotal role in promoting the bone destruction process with a positive feedback loop in the bone metastasis microenvironment. In addition, CSC-specific marker genes BIRC5, PTTG1, CENPF and CDKN3 were observed to be correlated with poor prognosis of CDRCC. Finally, we pinpointed that PARP, PIGF, HDAC2, and FGFR inhibitors for effectively targeting CSCs may be the potential therapeutic strategies for CDRCC. Conclusion: The results of the present study may shed new light on the identification of CSCs, and help further understand the mechanism underlying drug resistance, differentiation and metastasis in human CDRCC.