Establishment and validation of an aging-related risk signature associated with prognosis and tumor immune microenvironment in breast cancer.

Breast cancer (BC) is a highly malignant and heterogeneous tumor which is currently the cancer with the highest incidence and seriously endangers the survival and prognosis of patients. Aging, as a research hotspot in recent years, is widely considered to be involved in the occurrence and development of a variety of tumors. However, the relationship between aging-related genes (ARGs) and BC has not yet been fully elucidated. The expression profiles and clinicopathological data were acquired in the Cancer Genome Atlas (TCGA) and the gene expression omnibus (GEO) database. Firstly, the differentially expressed ARGs in BC and normal breast tissues were investigated. Based on these differential genes, a risk model was constructed composed of 11 ARGs via univariate and multivariate Cox analysis. Subsequently, survival analysis, independent prognostic analysis, time-dependent receiver operating characteristic (ROC) analysis and nomogram were performed to assess its ability to sensitively and specifically predict the survival and prognosis of patients, which was also verified in the validation set. In addition, functional enrichment analysis and immune infiltration analysis were applied to reveal the relationship between the risk scores and tumor immune microenvironment, immune status and immunotherapy. Finally, multiple datasets and real-time polymerase chain reaction (RT-PCR) were utilized to verify the expression level of the key genes. An 11-gene signature (including FABP7, IGHD, SPIB, CTSW, IGKC, SEZ6, S100B, CXCL1, IGLV6-57, CPLX2 and CCL19) was established to predict the survival of BC patients, which was validated by the GEO cohort. Based on the risk model, the BC patients were divided into high- and low-risk groups, and the high-risk patients showed worse survival. Stepwise ROC analysis and Cox analyses demonstrated the good performance and independence of the model. Moreover, a nomogram combined with the risk score and clinical parameters was built for prognostic prediction. Functional enrichment analysis revealed the robust relationship between the risk model with immune-related functions and pathways. Subsequent immune microenvironment analysis, immunotherapy, etc., indicated that the immune status of patients in the high-risk group decreased, and the anti-tumor immune function was impaired, which was significantly different with those in the low-risk group. Eventually, the expression level of FABP7, IGHD, SPIB, CTSW, IGKC, SEZ6, S100B, CXCL1, IGLV6-57 and CCL19 was identified as down-regulated in tumor cell line, while CPLX2 up-regulated, which was mostly similar with the results in TCGA and Human Protein Atlas (HPA) via RT-PCR. In summary, our study constructed a risk model composed of ARGs, which could be used as a solid model for predicting the survival and prognosis of BC patients. Moreover, this model also played an important role in tumor immunity, providing a new direction for patient immune status assessment and immunotherapy selection.

Wang Z ，Liu H ，Gong Y ，Cheng Y ... -  《-》

Derivation and Comprehensive Analysis of Aging Patterns in Patients with Bladder Cancer.

Aging is an essential risk factor for cancer. However, aging-related genes (ARGs) have not been comprehensively analyzed in bladder cancer (BC). Therefore, the study is aimed at derivating a risk stratification system for BC patients based on ARGs. Public databases were used to acquire ARGs sets, transcriptome files, and clinical data. The "limma" package was then used to screen for differential ARGs while also using univariate Cox regression analysis to explore for prognostic ARGs. The "ConsensusClusterPlus" package was used to perform aging patterns in BC patients based on the above prognostic ARGs. Subsequently, aging patterns were investigated in survival prediction, mutation landscape, immunotherapy, immunological checkpoints, and immune microenvironment. We likewise utilized gene enrichment analysis to explore the biological functions that were behind the findings. To construct a risk signature and nonogram for prognostic prediction, we used LASSO and Cox regression analysis based on differential genes in aging patterns. In addition, we plotted a nomogram and validate the accuracy of the risk signature in GEO and TCGA cohorts. We explored the possible biological mechanism using GSEA analysis and preliminarily identified a hub gene using PPI network. Finally, we validated the expression of hub gene in BC cell lines. We screened 84 downregulated ARGs, 74 upregulated ARGs, and 32 prognostic ARGs in the human aging genome resource. The aging patterns based on prognostic genes had excellent survival prediction (p < 0.001) and discriminatory ability in 405 BC patients. In addition, we found no significant differences in aging patterns in mutation analysis, which were all characterized by TP53, TTN, and KMT2D mutations. It is worth noting that cluster B in the aging patterns has a better response to immunotherapy and a more active immune microenvironment (p < 0.05). In addition, gene enrichment analysis showed that aging patterns may be related to biological processes such as Staphylococcus aureus infection, phagosome, and cytokine-cytokine receptor interaction. Subsequently, we constructed a risk signature based on 16 differential genes from different aging patterns and had good survival prediction ability in both GEO and TCGA cohort. Specifically, survival analysis revealed a significantly shorter survival time in the high-risk group than in the low-risk group (TCGA and GEO, p < 0.001). In addition, AUC values in the ROC analysis predicted 1, 3, and 5 years in TCGA cohort that are 0.713, 0.714, and 0.738, respectively. AUC values predicted 1, 3, and 5 years in GEO cohort that are 0.606, 0.663, and 0.718, respectively. There is no doubt that risk score was an independent prognostic factor from results of multivariate Cox regression analysis in BC patients (p < 0.001). There were also significant differences in immune cell infiltration, immune checkpoint, and immune score between the two groups (p < 0.05), but it should not be ignored that the correlation with the HLA expression was weak. Finally, we identified and validated CLIC3 as a hub gene that may be involved in the Wnt signaling pathway, etc. We provided robust evidences that aging patterns based on ARGs can guide targeted therapy and survival prediction in BC patients.

Wang B ，Tong F ，Zhai C ，Wang L ，Liu Y ，Wang J ... -  《-》

Anoikis-related genes in breast cancer patients: reliable biomarker of prognosis.

Tang M ，Rong Y ，Li X ，Pan H ，Tao P ，Wu Z ，Liu S ，Tang R ，Liu Z ，Cai H ... -  《BMC CANCER》

Construction and validation of a bladder cancer risk model based on autophagy-related genes.

Autophagy has an important association with tumorigenesis, progression, and prognosis. However, the mechanism of autophagy-regulated genes on the risk prognosis of bladder cancer (BC) patients has not been fully elucidated yet. In this study, we created a prognostic model of BC risk based on autophagy-related genes, which further illustrates the value of genes associated with autophagy in the treatment of BC. We first downloaded human autophagy-associated genes and BC datasets from Human Autophagy Database and The Cancer Genome Atlas (TCGA) database, and finally obtained differential prognosis-associated genes for autophagy by univariate regression analysis and differential analysis of cancer versus normal tissues. Subsequently, we downloaded two datasets from Gene Expression Omnibus (GEO), GSE31684 and GSE15307, to expand the total number of samples. Based on these genes, we distinguished the molecular subtypes (C1, C2) and gene classes (A, B) of BC by consistent clustering analysis. Using the genes merged from TCGA and the two GEO datasets, we conducted least absolute shrinkage and selection operator (LASSO) and multivariate Cox regression analysis to obtain risk genes and construct autophagy-related risk prediction models. The accuracy of this risk prediction model was assessed by receiver operating characteristic (ROC) and calibration curves, and then nomograms were constructed to predict the survival of bladder cancer patients at 1, 3, and 5 years, respectively. According to the median value of the risk score, we divided BC samples into the high- and low-risk groups. Kaplan-Meier (K-M) survival analysis was performed to compare survival differences between subgroups. Then, we used single sample gene set enrichment analysis (ssGSEA) for immune cell infiltration abundance, immune checkpoint genes, immunotherapy response, gene ontology (GO) and Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis, and tumor mutation burden (TMB) analysis for different subgroups. We also applied quantitative real-time polymerase chain reaction (PCR) and immunohistochemistry (IHC) techniques to verify the expression of these six genes in the model. Finally, we chose the IMvigor210 dataset for external validation. Six risk genes associated with autophagy (SPOCD1, FKBP10, NAT8B, LDLR, STMN3, and ANXA2) were finally screened by LASSO regression algorithm and multivariate Cox regression analysis. ROC and calibration curves showed that the model established was accurate and reliable. Univariate and multivariate regression analyses were used to verify that the risk model was an independent predictor. K-M survival analysis indicated that patients in the high-risk group had significantly worse overall survival than those in the low-risk group. Analysis by algorithms such as correlation analysis, gene set variation analysis (GSVA), and ssGSEA showed that differences in immune microenvironment, enrichment of multiple biologically active pathways, TMB, immune checkpoint genes, and human leukocyte antigens (HLAs) were observed in the different risk groups. Then, we constructed nomograms that predicted the 1-, 3-, and 5-year survival rates of different BC patients. In addition, we screened nine sensitive chemotherapeutic drugs using the correlation between the obtained expression status of risk genes and drug sensitivity results. Finally, the external dataset IMvigor210 verified that the model is reliable and efficient. We established an autophagy-related risk prognostic model that is accurate and reliable, which lays the foundation for future personalized treatment of bladder cancer.

Shen C ，Yan Y ，Yang S ，Wang Z ，Wu Z ，Li Z ，Zhang Z ，Lin Y ，Li P ，Hu H ... -  《-》

Identification of Novel Anoikis-Related Gene Signatures to Predict the Prognosis, Immune Microenvironment, and Drug Sensitivity of Breast Cancer Patients.

Liu J ，Wu H ，Wang Q ，Jin S ，Hou S ，Shen Z ，Zhao L ，Xu S ，Pang D ... -  《-》