Abstract
BACKGROUND: The mechanism underlying cervical carcinogenesis that is mediated by persistent human papillomavirus (HPV) infection remains elusive. AIMS: Here, for the first time, we deciphered both the temporal transition and spatial distribution of cellular subsets during disease progression from normal cervix tissues to precursor lesions to cervical cancer. MATERIALS & METHODS: We generated scRNA-seq profiles and spatial transcriptomics data from nine patient samples, including two HPV-negative normal, two HPV-positive normal, two HPV-positive HSIL and three HPV-positive cancer samples. RESULTS: We not only identified three 'HPV-related epithelial clusters' that are unique to normal, high-grade squamous intraepithelial lesions (HSIL) and cervical cancer tissues but also discovered node genes that potentially regulate disease progression. Moreover, we observed the gradual transition of multiple immune cells that exhibited positive immune responses, followed by dysregulation and exhaustion, and ultimately established an immune-suppressive microenvironment during the malignant program. In addition, analysis of cellular interactions further verified that a 'homeostasis-balance-malignancy' change occurred within the cervical microenvironment during disease progression. DISCUSSION: We for the first time presented a spatiotemporal atlas that systematically described the cellular heterogeneity and spatial map along the four developmental steps of HPV-related cervical oncogenesis, including normal, HPV-positive normal, HSIL and cancer. We identified three unique HPV-related clusters, discovered critical node genes that determined the cell fate and uncovered the immune remodeling during disease escalation. CONCLUSION: Together, these findings provided novel possibilities for accurate diagnosis, precise treatment and prognosis evaluation of patients with precancer and cervical cancer.