Tractable mouse TNBC models capture the heterogeneous tumor immune microenvironment and adaptation to PD-L1 blockade

Triple-negative breast cancer (TNBC) patients exhibit variable responses to programmed death (PD)-ligand (L)1 blockade, largely determined by the 'hot' versus 'cold' state of the tumor immune microenvironment (TIME). We here characterized 9 mouse TNBC models, relying on intraductal mammary gland inoculation of established mouse TNBC cell lines, with a heterogeneous TIME to study anti-PD-L1 resistance mechanisms. Complementary in vitro and in vivo screening classified the 4T1-hot-based model, a highly inflamed control through its immunogenic luciferase tag expression compared to the untagged 4T1-cold-based model, as displaying the 'hottest' TIME. However, both 4T1-based counterparts did not respond to anti-PD-L1, which was attributed to their immunosuppressive myeloid cell content as well as upregulation of cancer-associated fibroblasts in the 4T1-hot and high PD-L1-expressing CXCL10+ tumor-associated macrophages in 4T1-cold primary tumors. These anti-PD-L1 adaptation mechanisms across TIME states as captured by mouse TNBC models highlight specific cellular targets for future studies. Overall design: 9 Different TNBC models were characterized on different levels, with as a first step genomic characterization of the 9 selected TNBC celllines (4T1-hot, 4T1-cold, EMT6, D2A1, 66CL4, PY230, AT3, MVT1 and MET1) in vitro using single-cell RNA-sequencing. Creating a single-cell ATLAS of these 9 cell lines, analyzing transcriptomic differences between these cell lines identifying their heterogeneity based on their transcriptome on the in vitro level. The scRNAseq was performed in collaboration with the VIB Single-cell core.