Identification of m6A modification patterns and development of m6A-hypoxia prognostic signature to characterize tumor microenvironment in triple-negative breast cancer.

N6-methylation (m6A) modification of RNA has been found to have essential effects on aspects of the tumor microenvironment (TME) including hypoxia status and mobilization of immune cells. However, there are no studies to explore the combined effect of m6A modification and hypoxia on molecular heterogeneity and TME of triple-negative breast cancer (TNBC). We collected The Cancer Genome Atlas (TCGA-TNBC, N=139), the Molecular Taxonomy of Breast Cancer International Consortium (METABRIC-TNBC, N=297), the GSE103091, GSE21653, and GSE135565 series from the Gene Expression Omnibus (GEO-TNBC, N=247), and FUSCCTNBC (N=245) for our study. The non-negative matrix factorization algorithm was used to cluster TNBC samples. Immune cell infiltration was analyzed by the CIBERSORT algorithm. The enrichment scores were calculated by single-sample gene set enrichment analysis(ssGSEA) to characterize TME in TNBC samples. Immunohistochemistry (IHC) and qRT-PCR were performed to detect the gene expression. Based on the expression of m6A-related genes, we identified three distinct m6A clusters (denoted A, B, and C) in TNBC samples. Comparing the TME characteristics among the three clusters, we observed that cluster C was strongly related to hypoxia status and immune suppression, whereas clusters A and B displayed more immune cell infiltration. Therefore, we combine m6A and hypoxia related genes to classify two m6A-hypoxia clusters of TNBC and screened six prognostic genes by LASSO-Cox regression to construct a m6A-hypoxia signature(MHPS), which divided TNBC samples into high- and low-risk groups. We identified different TME features, immune cell infiltration between the two groups, and a better immunotherapy response was observed in the low-risk group. A nomogram was constructed with tumor size, lymph node, and risk score to improve clinical application of MHPS. We identified distinct TME characteristics of TNBC based on three different m6A modification patterns. Then, we constructed a specific m6A-hypoxia signature for TNBC to evaluate risk and predict immunotherapy response of patients, to enable more accurate treatment in the future.

Shen X ，Zhong J ，He J ，Han J ，Chen N ... -  《Frontiers in Immunology》

A combined hypoxia and immune gene signature for predicting survival and risk stratification in triple-negative breast cancer.

Yang X ，Weng X ，Yang Y ，Zhang M ，Xiu Y ，Peng W ，Liao X ，Xu M ，Sun Y ，Liu X ... -  《Aging-US》

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

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

The significance of m6A RNA methylation modification in prognosis and tumor microenvironment immune infiltration of cervical cancer.

Recent studies have highlighted that N6-methyladenosine (m6A) plays a significant role in tumorigenicity and progression. However, the mechanism of m6A modifications in the tumor microenvironment (TME) immune cell infiltration in cervical cancer (CC) remains unclear. Clinical and RNA sequencing data of 25 m6A RNA methylation regulators were acquired from the Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) database. LASSO Cox regression analysis was used to generate a prognostic risk signature. m6A modification patterns were identified based on the expression of 25 m6A regulators, and their correlation with TME immune cell-infiltrating characterization was analyzed. Principal component analysis was used to construct an m6A-scoring signature (m6A score) to evaluate the m6A modification patterns of individual CC samples and guide the selection of more effective immunotherapeutic strategies. Genetic and expression alterations of 25 m6A regulators were highly heterogeneous between CC and normal tissues. METTL14 and IGF2BP1 were selected to conduct the prognostic risk signature. Three m6A modification patterns were identified in 659 CC samples, which were associated with distinct clinical outcomes and biological pathways. The TME immune cell-infiltrating characterization of the three m6A modification patterns was highly consistent with 3 tumor immune phenotypes, including immune-excluded, immune-inflamed, and immune-desert phenotypes. Due to the heterogeneity of m6A modification patterns, an m6A scoring signature was established to evaluate the m6A modification patterns of individual CC samples. Univariate and multivariate Cox regression analyses revealed that the m6A score is a robust and independent prognostic biomarker for assessing the prognosis of CC patients. A low m6A score, characterized by higher somatic mutation and higher expression of proliferation-related and DNA repair-related genes, indicated poor overall survival. Activation of immune infiltration was exhibited by the high m6A score, which was likely to have a good response and clinical benefits to antiPD-1/L1 immunotherapy. This study highlights the prognostic value of 25 m6A regulators in CC. The m6A modification is related to immune regulation and the formation of TME heterogeneity and complexity. An m6A scoring signature to clarify the individual m6A modification pattern could enhance our understanding of TME immune cell-infiltrating characterization and guide immunotherapy strategies.

Guo Y ，Bai Y ，Wang L ，Xu Z ，Wang X ，Wang W ... -  《-》

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》