Development of a machine learning-based radiomics signature for estimating breast cancer TME phenotypes and predicting anti-PD-1/PD-L1 immunotherapy response

Abstract Backgrounds Since breast cancer patients respond diversely to immunotherapy, there is an urgent need to explore novel biomarkers to precisely predict clinical responses and enhance therapeutic efficacy. The purpose of our present research was to construct and independently validate a biomarker of tumor microenvironment (TME) phenotypes via a machine learning-based radiomics way. The interrelationship between the biomarker, TME phenotypes and recipients’ clinical response was also revealed. Methods In this retrospective multi-cohort investigation, five separate cohorts of breast cancer patients were recruited to measure breast cancer TME phenotypes via a radiomics signature, which was constructed and validated by integrating RNA-seq data with DCE-MRI images for predicting immunotherapy response. Initially, we constructed TME phenotypes using RNA-seq of 1089 breast cancer patients in the TCGA database. Then, parallel DCE-MRI images and RNA-seq of 94 breast cancer patients obtained from TCIA were applied to develop a radiomics-based TME phenotypes signature using random forest in machine learning. The repeatability of the radiomics signature was then validated in an internal validation set. Two additional independent external validation sets were analyzed to reassess this signature. The Immune phenotype cohort ( n = 158) was divided based on CD8 cell infiltration into immune-inflamed and immune-desert phenotypes; these data were utilized to examine the relationship between the immune phenotypes and this signature. Finally, we utilized an Immunotherapy-treated cohort with 77 cases who received anti-PD-1/PD-L1 treatment to evaluate the predictive efficiency of this signature in terms of clinical outcomes. Results The TME phenotypes of breast cancer were separated into two heterogeneous clusters: Cluster A, an "immune-inflamed" cluster, containing substantial innate and adaptive immune cell infiltration, and Cluster B, an "immune-desert" cluster, with modest TME cell infiltration. We constructed a radiomics signature for the TME phenotypes ([AUC] = 0.855; 95% CI 0.777–0.932; p < 0.05) and verified it in an internal validation set (0.844; 0.606–1; p < 0.05). In the known immune phenotypes cohort, the signature can identify either immune-inflamed or immune-desert tumor (0.814; 0.717–0.911; p < 0.05). In the Immunotherapy-treated cohort, patients with objective response had higher baseline radiomics scores than those with stable or progressing disease ( p < 0.05); moreover, the radiomics signature achieved an AUC of 0.784 (0.643–0.926; p < 0.05) for predicting immunotherapy response. Conclusions Our imaging biomarker, a practicable radiomics signature, is beneficial for predicting the TME phenotypes and clinical response in anti-PD-1/PD-L1-treated breast cancer patients. It is particularly effective in identifying the "immune-desert" phenotype and may aid in its transformation into an "immune-inflamed" phenotype.