Regulatory T cells are associated with the tumor immune microenvironment and immunotherapy response in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is the most aggressive subtype of breast cancer with a high risk of distant metastasis, an extremely poor prognosis, and a high risk of death. Regulatory T cells (Tregs) contribute to the formation of a tumor immunosuppressive microenvironment, which plays an important role in the progression and treatment resistance of TNBC. A public single-cell sequencing dataset demonstrated increased infiltration of Tregs in TNBC tissues relative to normal breast tissue. Weighted gene co-expression network analysis was used to identify Treg infiltration-related modules for METABRIC TNBC samples. Subsequently, we obtained two Treg infiltration-associated clusters of TNBC by applying consensus clustering and further constructed a prognostic model based on this Treg infiltration-associated gene module. The ability of the selected gene in the prognostic model, thymidine kinase-1 (TK1), to promote the progression of TNBC was evaluated in vitro. We concluded that two Treg infiltration-associated clusters had different prognoses and sensitivities to drugs commonly used in breast cancer treatment, and multi-omics analysis revealed that the two clusters had different copy number variations of key tumor progression genes. The 7-gene risk score based on TNBC Treg infiltration was a reliable prognostic indicator both in the training and validation cohorts. Moreover, patients with TNBC with high Treg infiltration-related scores lacked the activation of immune activation pathways and exhibited resistance to anti-PD1 immunotherapy. Knocking down TK1 led to impaired proliferation, migration, and invasion of TNBC cells in vitro. In addition, specimens from patients with TNBC with high TK1 expression showed significantly higher Treg infiltration in tumors. Results of spatial transcriptome analysis showed that TK1 positive cells mainly localize in tumor area, and Treg cell infiltration in TNBC tissues was associated with high expression of TK1. Pan-cancer analysis also demonstrated that TK1 is associated with poor prognosis and activation of proliferation pathways in multiple cancers. We established a prognostic model related to Treg infiltration and this model can be used to establish a clinically relevant classification of TNBC progression. Additionally, our work revealed the underestimable potential of TK1 as a tumor biomarker and immunotherapeutic target.

Huang P ，Zhou X ，Zheng M ，Yu Y ，Jin G ，Zhang S ... -  《Frontiers in Immunology》

3D Collagen Fiber Concentration Regulates Treg Cell Infiltration in Triple Negative Breast Cancer.

Triple negative breast cancer (TNBC) is characterized by poor prognosis and a lack of effective therapeutic agents owing to the absence of biomarkers. A high abundance of tumor-infiltrating regulatory T cells (Tregs) was associated with worse prognosis in malignant disease. Exploring the association between Treg cell infiltration and TNBC will provide new insights for understanding TNBC immunosuppression and may pave the way for developing novel immune-based treatments. Patients from TCGA were divided into Treg-high (Treg-H) and Treg-low (Treg-L) groups based on the abundance of Tregs according to CIBERSORT analysis. The association between expression level of Tregs and the clinical characteristics as well as prognosis of breast cancer were evaluated. Next, a Treg-related prognostic model was established after survival-dependent univariate Cox and LASSO regression analysis, companied with an external GEO cohort validation. Then, GO, KEGG and GSEA analyses were performed between the Treg-H and Treg-L groups. Masson and Sirius red/Fast Green staining were applied for ECM characterization. Accordingly, Jurkat T cells were encapsulated in 3D collagen to mimic the ECM microenvironment, and the expression levels of CD4, FOXP3 and CD25 were quantified according to immunofluorescence staining. The expression level of Tregs is significantly associated with the clinical characteristics of breast cancer patients, and a high level of Treg cell expression indicates a poor prognosis in TNBC. To further evaluate this, a Treg-related prognostic model was established that accurately predicted outcomes in both TCGA training and GEO validation cohorts of TNBC patients. Subsequently, ECM-associated signaling pathways were identified between the Treg-H and Treg-L groups, indicating the role of ECM in Treg infiltration. Since we found increasing collagen concentrations in TNBC patients with distant migration, we encapsulated Jurkat T cells within a 3D matrix with different collagen concentrations and observed that increasing collagen concentrations promoted the expression of Treg biomarkers, supporting the regulatory role of ECM in Treg infiltration. Our results support the association between Treg expression and breast cancer progression as well as prognosis in the TNBC subtype. Moreover, increasing collagen density may promote Treg infiltration, and thus induce an immunosuppressed TME.

Gao H ，Tian Q ，Zhou Y ，Zhu L ，Lu Y ，Ma Y ，Feng J ，Jiang Y ，Wang B ... -  《Frontiers in Immunology》

Research and experimental verification on the mechanisms of cellular senescence in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) is an aggressive breast cancer subtype with high heterogeneity, poor prognosis, and a low 10-year survival rate of less than 50%. Although cellular senescence displays extensive effects on cancer, the comprehensions of cellular senescence-related characteristics in TNBC patients remains obscure. Single-cell RNA sequencing (scRNA-seq) data were analyzed by Seurat package. Scores for cellular senescence-related pathways were computed by single-sample gene set enrichment analysis (ssGSEA). Subsequently, unsupervised consensus clustering was performed for molecular cluster identification. Immune scores of patients in The Cancer Genome Atlas (TCGA) dataset and associated immune cell scores were calculated using Estimation of STromal and Immune cells in MAlignantTumours using Expression data (ESTIMATE) and Microenvironment Cell Populations-counter (MCP-counter), Tumor Immune Estimation Resource (TIMER) and Estimating the Proportion of Immune and Cancer cells (EPIC) methods, respectively. Immunotherapy scores were assessed using TIDE. Furthermore, feature genes were identified by univariate Cox and Least Absolute Shrinkage and Selection Operator (LASSO) regression analyses; these were used to construct a risk model. Additionally, quantitative reverse transcription-polymerase chain reaction (qRT-PCR) and transwell assay were conducted for in vitro validation of hub genes. TNBC was classified into three subtypes based on cellular senescence-related pathways as clusters 1, 2, and 3. Specifically, cluster 1 showed the best prognosis, followed by cluster 2 and cluster 3. The levels of gene expression in cluster 2 were the lowest, whereas these were the highest in cluster 3. Moreover, clusters 1 and 3 showed a high degree of immune infiltration. TIDE scores were higher for cluster 3, suggesting that immune escape was more likely in patients with the cluster 3 subtype who were less likely to benefit from immunotherapy. Next, the TNBC risk model was constructed and validated. RT-qPCR revealed that prognostic risk genes (MMP28, ACP5 and KRT6A) were up-regulated while protective genes (CT83) were down-regulated in TNBC cell lines, validating the results of the bioinformatics analysis. Meanwhile, cellular experiments revealed that ACP5 could promote the migration and invasion abilities in two TNBC cell lines. Finally, we evaluated the validity of prognostic models for assessing TME characteristics and TNBC chemotherapy response. In conclusion, these findings help to assess the efficacy of targeted therapies in patients with different molecular subtypes, have practical applications for subtype-specific treatment of TNBC patients, and provide information on prognostic factors, as well as guidance for the revelation of the molecular mechanisms by which senescence-associated genes influence TNBC progression.

Cao T ，Huang M ，Huang X ，Tang T ... -  《PeerJ》

A novel fatty-acid metabolism-based classification for triple negative breast cancer.

Yang X ，Tang W ，He Y ，An H ，Wang J ... -  《Aging-US》

Development and verification of a novel immunogenic cell death-related signature for predicting the prognosis and immune infiltration in triple-negative breast cancer.

Insufficient understanding of the pathogenesis and tumor immunology of triple-negative breast cancer (TNBC) has limited the development of immunotherapy. The importance of tumor microenvironment (TME) in immunotyping, prognostic assessment and immunotherapy efficacy of cancer has been emphasized, however, potential immunogenic cell death (ICD) related genes function in TME of TNBC has been rarely investigated. To initially explore the role and related mechanisms of ICD in TNBC, especially the role played in the TME of TNBC, and to identify different relevant subtypes based on ICD, and then develop an ICD-related risk score to predict each TNBC patient TME status, prognosis and immunotherapy response. In this study, we identified distinct ICD-related modification patterns based on 158 TNBC cases in the TCGA-TNBC cohort. We then investigated the possible correlation between ICD-related modification patterns and TME cell infiltration characteristics in TNBC. By using univariate Cox and least absolute shrinkage and selection operator (LASSO) regression analysis, we created a risk scoring system (ICD score) to quantifiably evaluate the impact of ICD-related modification patterns in individual TNBC patient. Two different ICD-related modification patterns were found with significant differences in immune infiltration. Lower ICD score was correlated with higher immune infiltration, tumor mutational burden and significantly enriched in immune-related pathways, indicating a strong ability to activate immune response, which might account for relatively favorable prognosis of TNBC patients and could serve as a predictor to select suitable candidates for immunotherapy. We used two independent cohorts, GSE58812 cohort and Metabric cohort to validate prognosis and immunohistochemistry for preliminary in vitro validation. This study evidenced that the ICD-related modification patterns might exert pivotal roles in the immune infiltration landscape of TNBC and ICD score might act as potential predictors of prognostic assessment and immunotherapy response. This research provides unique insights for individualize immune treatment strategies and promising immunotherapy candidates screening.

Li J ，Li Z ，Yang W ，Pan J ，You H ，Yang L ，Zhang X ... -  《Cancer Reports》