Prediction of therapeutic response to transarterial chemoembolization plus systemic therapy regimen in hepatocellular carcinoma using pretreatment contrast-enhanced MRI based habitat analysis and Crossformer model.

Purpose: To develop habitat and deep learning (DL) models from multi-phase contrast-enhanced magnetic resonance imaging (CE-MRI) habitat images categorized using the K-means clustering algorithm. Additionally, we aim to assess the predictive value of identified regions for early evaluation of the responsiveness of hepatocellular carcinoma (HCC) patients to treatment with transarterial chemoembolization (TACE) plus molecular targeted therapies (MTT) and anti-PD-(L)1.

Methods: A total of 102 patients with HCC from two institutions (A, n = 63 and B, n = 39) who received TACE plus systemic therapy were enrolled from September 2020 to January 2024. Multiple CE-MRI sequences were used to outline 3D volumes of interest (VOI) of the lesion. Subsequently, K-means clustering was applied to categorize intratumoral voxels into three distinct subgroups, based on signal intensity values of images. Using data from institution A, the habitat model was built with the ExtraTrees classifier after extracting radiomics features from intratumoral habitats. Similarly, the Crossformer model and ResNet50 model were trained on multi-channel data in institution A, and a DL model with Transformer-based aggregation was constructed to predict the response. Finally, all models underwent validation at institution B.

Results: The Crossformer model and the habitat model both showed high area under the receiver operating characteristic curves (AUCs) of 0.869 and 0.877 (training cohort). In validation, AUC was 0.762 for the Crossformer model and 0.721 for the habitat model.

Conclusion: The habitat model and DL model based on CE-MRI possesses the capability to non-invasively predict the efficacy of TACE plus systemic therapy in HCC patients, which is critical for precision treatment and patient outcomes.