Single-cell transcriptome analysis profiling lymphatic invasion-related TME in colorectal cancer

Lymphatic invasion (LI) is extremely aggressive and induces worse prognosis among patients with colorectal cancer (CRC). Thus, it is critical to characterize the cellular and molecular mechanisms underlying LI in order to establish novel and efficacious therapeutic targets that enhance the prognosis of CRC patients. RNA-seq data, clinical and survival information of colon adenocarcinoma (COAD) patients were obtained from the TCGA database. In addition, three scRNA-seq datasets of CRC patients were acquired from the GEO database. Data analyses were conducted with the R packages. We assessed the tumor microenvironment (TME) differences between LI+ and LI- based scRNA-seq data, LI+ cells exhibited augmented abundance of immunosuppression and invasive subset. Marked extracellular matrix network activation was also observed in LI+ cells within SPP1+ macrophages. We revealed that an immunosuppressive and pro-angiogenic TME strongly enhanced LI, as was evidenced by the CD4+ Tregs, CD8+ GZMK+, SPP1+ macrophages, e-myCAFs, and w-myCAFs subcluster infiltrations. Furthermore, we identified potential LI targets that influenced tumor development, metastasis, and immunotherapeutic response. Finally, a novel LIRS model was established based on the expression of 14 LI-related signatures, and in the two testing cohorts, LIRS was also proved to have accurate prognostic predictive ability. In this report, we provided a valuable resource and extensive insights into the LI of CRC. Our conclusions can potentially benefit the establishment of highly efficacious therapeutic targets as well as diagnostic biomarkers that improve patient outcomes.