Identifying immune infiltration by deep learning to assess the prognosis of patients with hepatocellular carcinoma

The treatment situation for hepatocellular carcinoma remains critical. The use of deep learning algorithms to assess immune infiltration is a promising new diagnostic tool.

Patient data and whole slide images (WSIs) were obtained for the Xijing Hospital (XJH) cohort and TCGA cohort. We wrote programs using Visual studio 2022 with C# language to segment the WSI into tiles. Pathologists classified the tiles and later trained deep learning models using the ResNet 101V2 network via ML.NET with the TensorFlow framework. Model performance was evaluated using AccuracyMicro versus AccuracyMacro. Model performance was examined using ROC curves versus PR curves. The percentage of immune infiltration was calculated using the R package survminer to calculate the intergroup cutoff, and the Kaplan‒Meier method was used to plot the overall survival curve of patients. Cox regression was used to determine whether the percentage of immune infiltration was an independent risk factor for prognosis. A nomogram was constructed, and its accuracy was verified using time-dependent ROC curves with calibration curves. The CIBERSORT algorithm was used to assess immune infiltration between groups. Gene Ontology was used to explore the pathways of differentially expressed genes.

There were 100 WSIs and 165,293 tiles in the training set. The final deep learning models had an AccuracyMicro of 97.46% and an AccuracyMacro of 82.28%. The AUCs of the ROC curves on both the training and validation sets exceeded 0.95. The areas under the classification PR curves exceeded 0.85, except that of the TLS on the validation set, which might have had poor results (0.713) due to too few samples. There was a significant difference in OS between the TIL classification groups (p < 0.001), while there was no significant difference in OS between the TLS groups (p = 0.294). Cox regression showed that TIL percentage was an independent risk factor for prognosis in HCC patients (p = 0.015). The AUCs according to the nomogram were 0.714, 0.690, and 0.676 for the 1-year, 2-year, and 5-year AUCs in the TCGA cohort and 0.756, 0.797, and 0.883 in the XJH cohort, respectively. There were significant differences in the levels of infiltration of seven immune cell types between the two groups of samples, and gene ontology showed that the differentially expressed genes between the groups were immune related. Their expression levels of PD-1 and CTLA4 were also significantly different.

We constructed and tested a deep learning model that evaluates the immune infiltration of liver cancer tissue in HCC patients. Our findings demonstrate the value of the model in assessing patient prognosis, immune infiltration and immune checkpoint expression levels.