The incidence and mortality of colorectal cancer (CRC) are increasing year by year. The accurate classification of CRC can realize the purpose of personalized and precise treatment for patients. The tumor microenvironment (TME) plays an important role in the malignant progression and immunotherapy of CRC. An in-depth understanding of the clusters based on the TME is of great significance for the discovery of new therapeutic targets for CRC. We extracted data on CRC, including gene expression profile, DNA methylation array, somatic mutations, clinicopathological information, and copy number variation (CNV), from The Cancer Genome Atlas (TCGA), Gene Expression Omnibus (GEO) (four datasets—GSE14333, GSE17538, GSE38832, and GSE39582), cBioPortal, and FireBrowse. The MCPcounter was utilized to quantify the abundance of 10 TME cells for CRC samples. Cluster repetitive analysis was based on the Hcluster function of the Pheatmap package in R. The ESTIMATE package was applied to compute immune and stromal scores for CRC patients. PCA analysis was used to remove batch effects among different datasets and transform genome-wide DNA methylation profiling into methylation of tumor-infiltrating lymphocyte (MeTIL). We evaluated the mutation differences of the clusters using MOVICS, DeconstructSigs, and GISTIC packages. As for therapy, TIDE and SubMap analyses were carried out to forecast the immunotherapy response of the clusters, and chemotherapeutic sensibility was estimated based on the pRRophetic package. All results were verified in the TCGA and GEO data. Four immune clusters (ImmClust-CS1, ImmClust-CS2, ImmClust-CS3, and ImmClust-CS4) were identified for CRC. The four ImmClusts exhibited distinct TME compositions, cancer-associated fibroblasts (CAFs), functional orientation, and immune checkpoints. The highest immune, stromal, and MeTIL scores were observed in CS2, in contrast to the lowest scores in CS4. CS1 may respond to immunotherapy, while CS2 may respond to immunotherapy after anti-CAFs. Among the four ImmClusts, the top 15 markers with the highest mutation frequency were acquired, and CS1 had significantly lower CNA on the focal level than other subtypes. In addition, CS1 and CS2 patients had more stable chromosomes than CS3 and CS4. The most sensitive chemotherapeutic agents in these four ImmClusts were also found. IHC results revealed that CD29 stained significantly darker in the cancer samples, indicating that their CD29 was highly expressed in colon cancer. This work revealed the novel clusters based on TME for CRC, which would guide in predicting the prognosis, biological features, and appropriate treatment for patients with CRC.

We downloaded gene expression data, clinical data, and somatic mutation data of cervical squamous cell carcinoma (CSCC) patients from The Cancer Genome Atlas (TCGA) and Genotype-Tissue Expression (GTEx) databases. Predictive lncRNAs were screened using univariate analysis and lasso regression, and risk score of each patient were calculated according to the expression levels of lncRNAs and regression coe cients to establish a risk model that could be a novel signature. We assessed the correlation between immune in ltration status, chemotherapeutics sensitivity, immune checkpoint proteins (ICP) and the signature. Therefore, we selected 11 immune-related lncRNAs (WWC2,

Background. Ferroptosis, a form of cell death driven by iron-dependent lipid peroxidation, may be a potential treatment for many cancers, including cervical cancer (CC). However, the regulation of long noncoding RNAs (lncRNAs) in the process of ferroptosis and whether ferroptosis inducers could increase the cytotoxicity of conventional chemotherapy drugs remain to be further elucidated. Methods. We analyzed the variation of 55 differentially ferroptosis-related genes (FRGs) and the influence of mutations in CC patients. The patients with CC were classified into two ferroptosis clusters by the non-negative matrix factorization (NMF) algorithm. The principal components analysis (PCA) was used to measure the ferroptosis score (FerroScore) in patients with CC. Besides, FerroScore was used to predict the sensitivity of chemotherapy and responses to immunotherapy in patients with CC. Finally, experiments were performed to verify the regulatory effect of AC026790.1 on erastin-induced ferroptosis, as well as the effect of erastin in combination with cisplatin on the toxicity of CC cells (SiHa, HeLa). Result. There were significant differences in the overall survival and immune cell infiltration between the two ferroptosis clusters. Patients with low FerroScore were more sensitive to chemotherapy drugs such as cisplatin and docetaxel. The low-FerroScore group had higher CD8+ T cell infiltration and immune checkpoint expression, demonstrating that patients with lower FerroScores were more sensitive to immunotherapy, which was consistent with the result of the submap method. In vitro, overexpression of AC026790.1 could promote erastin-induced ferroptosis, and the combination of erastin and cisplatin could increase the toxicity of CC cells. Conclusion. FerroScore has a potential prognostic value for CC that may provide a reference for chemotherapy and immunotherapy. LncRNA AC026790.1 can influence ferroptosis, and the combination of ferroptosis inducers and chemotherapy drugs can provide a new perspective on cancer treatment.

Breast cancer (BRCA) remains the most prevalent cancer worldwide and the tumor microenvironment (TME) has been discovered to exert a wide influence on the overall survival and therapeutic response. Numerous lines of evidence reported that the effects of immunotherapy of BRCA were manipulated by TME. Immunogenic cell death (ICD) is a form of regulated cell death (RCD) that is capable of fueling adaptive immune responses and aberrant expression of ICD-related genes (ICDRGs) can govern the TME system by emitting danger signals or damage-associated molecular patterns (DAMPs). In the current study, we obtained 34 key ICDRGs in BRCA. Subsequently, using the transcriptome data of BRCA from the TCGA database, we constructed a risk signature based on 6 vital ICDRGs, which had a good performance in predicting the overall survival of BRCA patients. We also examined the efficacy of our risk signature in the validation dataset (GSE20711) in the GEO database and it performed excellently. According to the risk model, patients with BRCA were divided into high-risk and low-risk groups. Also, the unique immune characteristics and TME between the two subgroups and 10 promising small molecule drugs targeting BRCA patients with different ICDRGs risk have been investigated. The low-risk group had good immunity indicated by T cell infiltration and high immune checkpoint expression. Moreover, the BRCA samples could be divided into three immune subtypes according to immune response severity (ISA, ISB, and ISC). ISA and ISB predominated in the low-risk group and patients in the low-risk group exhibited a more vigorous immune response. In conclusion, we developed an ICDRGs-based risk signature that can predict the prognosis of BRCA patients and offer a novel therapeutic strategy for immunotherapy, which would be of great significance in the BRCA clinical setting.

Background. Breast cancer (BC) is a highly heterogeneous disease with high morbidity and mortality. Its subtypes may have distinctly different biological behaviors, clinical outcomes, and therapeutic responses. The metabolic status of BC tissue is closely related to its progress. Therefore, we comprehensively characterized the function of metabolic genes in BC and identified new biomarkers to predict BC patients’ prognoses. Methods. Metabolic genes were identified by intersecting genes obtained from two published pieces of literature. The function of metabolic genes in BC was determined by extracting differentially expressed genes (DEGs), performing functional enrichment analyses, analyzing the infiltrating proportion of immune cells, and conducting metabolic subgroup analyses. A risk score model was constructed to assess the prognoses of BC patients by performing the univariate Cox regression, LASSO algorithm, multivariate Cox regression, Kaplan-Meier survival analyses, and ROC curve analyses in the training set. The prognostic model was then validated on the testing dataset, external dataset, the whole TCGA-BC database, and our clinical specimens. Finally, a nomogram was constructed for clinical prognostic prediction based on the risk score model and other clinicopathological parameters. Results. 955 metabolic genes were obtained. Among these, 157 metabolic DEGs were identified between BC and normal tissues for subsequent GO and KEGG pathway enrichment analyses. 5 metabolic genes were negatively correlated with CD8+ T cells, while 49 genes were positively correlated with CD8+ T cells. Furthermore, 5 metabolic subgroups with varying proportions of PAM50 subtypes, TNM classification, and immune cell infiltration were obtained. Finally, a risk score model was constructed to predict the prognoses of BC patients, and a nomogram incorporating the risk score model was established for clinical application. Conclusion. In this study, we elucidated tumor heterogeneity from metabolite profiling of BC. The roles of metabolic genes in the occurrence of BC were comprehensively characterized, clarifying the relationship between the tumor microenvironment (TME) and metabolic genes. Meanwhile, a concise prediction model was also constructed based on metabolic genes, providing a convenient and precise method for the individualized diagnosis and treatment of BC patients.