Abstract

Tumor evolution is a key feature of tumorigenesis and plays a pivotal role in driving intratumor heterogeneity, treatment failure and patients prognosis. Here we performed single-cell transcriptome profiling of 46 primary liver cancers from 37 patients enrolled for interventional studies. We surveyed the landscape of ~57,000 malignant and non-malignant cells and determined tumor cell clonality by developing a machine learning-based consensus clustering method. We found evidence of tumor cell branching evolution using hierarchical clustering, RNA velocity as well as reverse graph embedding methods. Interestingly, an increasing tumor cell clonality was tightly linked to patients prognosis, accompanied by a polarized immune cell landscape. We identified osteopontin as a key player for tumor cell evolution and microenvironmental reprogramming. Our study offers insight into the collective behavior of tumor cell communities in liver cancer as well as potential drivers for tumor evolution in response to therapy.