Abstract
Metastasis is the major cause of death in patients with pancreatic ductal adenocarcinoma (PDAC), and circulating tumor cells (CTCs) play an important role in the development of metastasis. However, few studies have uncovered the metastasis mechanism of PDAC based on CTCs. In this study, the existing bulk RNA-sequencing (bulk RNA-seq) and single-cell sequencing (scRNA-seq) data for CTCs in pancreatic cancer were obtained from the Gene Expression Omnibus (GEO) database. Analysis of tumor-infiltrating immune cells (TIICs) by CIBERSORT showed that the CTCs enriched from the peripheral blood of metastatic PDAC were found to contain a high proportion of T cell regulators (Tregs) and macrophages, while the proportion of dendritic cells (DCs) was lower than that enriched from localized PDAC. Through weighted gene co-expression network analysis (WGCNA) and the result of scRNA-seq, we identified the hub module (265 genes) and 87 marker genes, respectively, which were highly associated with metastasis. The results of functional enrichment analysis indicated that the two gene sets mentioned above are mainly involved in cell adhesion and cytoskeleton and epithelial-mesenchymal transition (EMT). Finally, we found that HMGB3 was the hub gene according to the Venn diagram. The expression of HMGB3 in PDAC was significantly higher than that in normal tissues (protein and mRNA levels). HMGB3 expression was significantly positively correlated with both EMT-related molecules and CTC cluster-related markers. Furthermore, it was also found that HMGB3 mutations were favorably related to tumor-associated immune cells through the TIMER2.0 online tool. We further demonstrated that PDAC patients with higher HMGB3 expression had significantly worse overall survival (OS) in multiple datasets. In summary, our study suggests that HMGB3 is a hub gene associated with EMT in CTCs, the formation of CTC clusters, and infiltration patterns of immune cells favorable for tumor progression and metastasis to distant organs.