Spatial patterns and MRI-based radiomic prediction of high peritumoral tertiary lymphoid structure density in hepatocellular carcinoma: a multicenter study

Tertiary lymphoid structures (TLS) within the tumor microenvironment have been associated with cancer prognosis and therapeutic response. However, the immunological pattern of a high peritumoral TLS (pTLS) density and its clinical potential in hepatocellular carcinoma (HCC) remain poor. This study aimed to elucidate biological differences related to pTLS density and develop a radiomic classifier for predicting pTLS density in HCC, offering new insights for clinical diagnosis and treatment. Spatial transcriptomics (n=4) and RNA sequencing data (n=952) were used to identify critical regulators of pTLS density and evaluate their prognostic significance in HCC. Baseline MRI images from 660 patients with HCC who had undergone surgery treatment between October 2015 and January 2023 were retrospectively recruited for model development and validation. This included training (n=307) and temporal validation (n=76) cohorts from Xiangya Hospital, and external validation cohorts from three independent hospitals (n=277). Radiomic features were extracted from intratumoral and peritumoral regions of interest and analyzed using machine learning algorithms to develop a predictive classifier. The classifier's performance was evaluated using the area under the curve (AUC), with prognostic and predictive value assessed across four independent cohorts and in a dual-center outcome cohort of 41 patients who received immunotherapy. Patients with HCC and a high pTLS density experienced prolonged median overall survival (p<0.05) and favorable immunotherapy response (p=0.03). Moreover, immune infiltration by mature B cells was observed in the high pTLS density region. Spatial pseudotime analysis and immunohistochemistry staining revealed that expansion of pTLS in HCC was associated with elevated CXCL9 and CXCL10 co-expression. We developed an optimal radiomic-based classifier with excellent discrimination for predicting pTLS density, achieving an AUC of 0.91 (95% CI 0.87, 0.94) in the external validation cohort. This classifier also exhibited promising stratification ability in terms of overall survival (p<0.01), relapse-free survival (p<0.05), and immunotherapy response (p<0.05). We identified key regulators of pTLS density in patients with HCC and proposed a non-invasive radiomic classifier capable of assisting in stratification for prognosis and treatment.