Development and validation of an autophagy-related prognostic signature in esophageal cancer

Triphenylphosphonium-Conjugated Palmitic Acid for Mitochondrial Targeting of Pancreatic Cancer Cells: Proteomic and Molecular Evidence

Pancreatic ductal adenocarcinoma (PDAC)'s resistance to therapies is mainly attributed to pancreatic cancer stem cells (PCSCs). Mitochondria-impairing agents can be used to hamper PCSC propagation and reduce PDAC progression. Therefore, to develop an efficient vector for delivering drugs to the mitochondria, we synthesized tris(3,5-dimethylphenyl)phosphonium-conjugated palmitic acid. Triphenylphosphonium (TPP) is a lipophilic cationic moiety that promotes the accumulation of conjugated agents in the mitochondrion. Palmitic acid (PA), the most common saturated fatty acid, has pro-apoptotic activity in different types of cancer cells. TPP-PA was prepared by the reaction of 16-bromopalmitic acid with TPP, and its structure was characterized by 1H and 13C NMR and HRMS. We compared the proteomes of TPP-PA-treated and untreated PDAC cells and PCSCs, identifying dysregulated proteins and pathways. Furthermore, assessments of mitochondrial membrane potential, intracellular ROS, cardiolipin content and lipid peroxidation, ER stress, and autophagy markers provided information on the mechanism of action of TPP-PA. The findings showed that TPP-PA reduces PDAC cell proliferation through mitochondrial disruption that leads to increased ROS, activation of ER stress, and autophagy. Hence, TPP-PA might offer a new approach for eliminating both the primary population of cancer cells and PCSCs, which highlights the promise of TPP-derived compounds as anticancer agents for PDAC.

Comprehensive analysis of ZNF family genes in prognosis, immunity, and treatment of esophageal cancer

BACKGROUND

As a common malignant tumor, esophageal carcinoma (ESCA) has a low early diagnosis rate and poor prognosis. This study aimed to construct the prognostic features composed of ZNF family genes to effectively predict the prognosis of ESCA patients.

METHODS

The mRNA expression matrix and clinical data were downloaded from TCGA and GEO database. Using univariate Cox analysis, lasso regression and multivariate Cox analysis, we screened six prognosis-related ZNF family genes to construct the prognostic model. We then used Kaplan-Meier plot, time-dependent receiver operating characteristic (ROC), multivariable Cox regression analysis of clinical information, and nomogram to evaluate the prognostic value within and across sets, separately and combined. We also validated the prognostic value of the six-gene signature using GSE53624 dataset. The different immune status was observed in the single sample Gene Set Enrichment Analysis (ssGSEA). Finally, real-time quantitative PCR was used to detect the expression of six prognostic ZNF genes in twelve pairs of ESCA and adjacent normal tissues.

RESULTS

A six prognosis-related ZNF family genes model consisted of ZNF91, ZNF586, ZNF502, ZNF865, ZNF106 and ZNF225 was identified. Multivariable Cox regression analysis revealed that six prognosis-related ZNF family genes were independent prognostic factors for overall survival of ESCA patients in TCGA and GSE53624. Further, a prognostic nomogram including the riskScore, age, gender, T, stage was constructed, and TCGA/GSE53624-based calibration plots indicated its excellent predictive performance. Drug Sensitivity and ssGSEA analysis showed that the six genes model was closely related to immune cells infiltration and could be used as a potential predictor of chemotherapy sensitivity.

CONCLUSION

We identified six prognosis-related ZNF family genes model of ESCA, which provide evidence for individualized prevention and treatment.

Modeling tumour heterogeneity of PD-L1 expression in tumour progression and adaptive therapy

Although PD-1/PD-L1 inhibitors show potent and durable anti-tumour effects in some refractory tumours, the response rate in overall patients is unsatisfactory, which in part due to the inherent heterogeneity of PD-L1. In order to establish an approach for predicting and estimating the dynamic alternation of PD-L1 heterogeneity during cancer progression and treatment, this study establishes a comprehensive modelling and computational framework based on a mathematical model of cancer cell evolution in the tumour-immune microenvironment, and in combination with epigenetic data and overall survival data of clinical patients from The Cancer Genome Atlas. Through PD-L1 heterogeneous virtual patients obtained by the computational framework, we explore the adaptive therapy of administering anti-PD-L1 according to the dynamic of PD-L1 state among cancer cells. Our results show that in contrast to the continuous maximum tolerated dose treatment, adaptive therapy is more effective for PD-L1 positive patients, in that it prolongs the survival of patients by administration of drugs at lower dosage.

Predicting the prognosis of esophageal cancer based on extensive analysis of new inflammatory response‐related signature

The prognosis of esophageal cancer (ESCA) is very poor, with a 5-year survival rate of less than 20%. On the other hand, inflammation is the characteristic hallmark of ESCA; however, the prognostic relationship between inflammatory response-related genes and ESCA has not been clarified yet. Therefore, in the present manuscript, we intend to investigate the correlation and specific signature of inflammation for the prediction of the prognosis of ESCA. A total of 173 samples were obtained from The Cancer Genome Atlas (TCGA) database, including 162 tumors and 11 normal specimens. The prognostic signature was established by least absolute shrinkage and selection operator Cox regression analysis. The transcription factor regulatory network with genes of the prognostic signature was analyzed from the transcriptional regulatory relationships unravelled by sentence-based text-mining database. Chemotherapy sensitivity and immunotherapy analysis were also performed. Multivariate Cox analysis showed that the signature was an independent prognostic risk factor. The low-risk group had poorer outcomes than the high-risk group. In the high-risk group, the infiltration of most immune cells was high and strongly correlated with the riskScore. In chemotherapeutic drug sensitivity analysis, OSM, AHR, and BTG2 were significantly correlated with the current chemotherapeutic drugs of ESCA. We have demonstrated a valid prognostic signature of inflammatory response-related genes and found strong associations with immune cells, targeted genes, and chemotherapeutic agents.

Identification and Validation of Cuproptosis-Related LncRNA Signatures in the Prognosis and Immunotherapy of Clear Cell Renal Cell Carcinoma Using Machine Learning

(1) Objective: We aimed to mine cuproptosis-related LncRNAs with prognostic value and construct a corresponding prognostic model using machine learning. External validation of the model was performed in the ICGC database and in multiple renal cancer cell lines via qPCR. (2) Methods: TCGA and ICGC cohorts related to renal clear cell carcinoma were included. GO and KEGG analyses were conducted to determine the biological significance of differentially expressed cuproptosis-related LncRNAs (CRLRs). Machine learning (LASSO), Kaplan-Meier, and Cox analyses were conducted to determine the prognostic genes. The tumor microenvironment and tumor mutation load were further studied. TIDE and IC50 were used to evaluate the response to immunotherapy, a risk model of LncRNAs related to the cuproptosis genes was established, and the ability of this model was verified in an external independent ICGC cohort. LncRNAs were identified in normal HK-2 cells and verified in four renal cell lines via qPCR. (3) Results: We obtained 280 CRLRs and identified 66 LncRNAs included in the TCGA-KIRC cohort. Then, three hub LncRNAs (AC026401.3, FOXD2-AS1, and LASTR), which were over-expressed in the four ccRCC cell lines compared with the human renal cortex proximal tubule epithelial cell line HK-2, were identified. In the ICGC database, the expression of FOXD2-AS1 and LASTR was consistent with the qPCR and TCGA-KIRC. The results also indicated that patients with low-risk ccRCC-stratified by tumor-node metastasis stage, sex, and tumor grade-had significantly better overall survival than those with high-risk ccRCC. The predictive algorithm showed that, according to the three CRLR models, the low-risk group was more sensitive to nine target drugs (A.443654, A.770041, ABT.888, AG.014699, AMG.706, ATRA, AP.24534, axitinib, and AZ628), based on the estimated half-maximal inhibitory concentrations. In contrast, the high-risk group was more sensitive to ABT.263 and AKT inhibitors VIII and AS601245. Using the CRLR models, the correlation between the tumor immune microenvironment and cancer immunotherapy response revealed that high-risk patients are more likely to respond to immunotherapy than low-risk patients. In terms of immune marker levels, there were significant differences between the high- and low-risk groups. A high TMB score in the high-risk CRLR group was associated with worse survival, which could be a prognostic factor for KIRC. (4) Conclusions: This study elucidates the core cuproptosis-related LncRNAs, FOXD2-AS1, AC026401.3, and LASTR, in terms of potential predictive value, immunotherapeutic strategy, and outcome of ccRCC.

LASSO Model Better Predicted the Prognosis of DLBCL than Random Forest Model: A Retrospective Multicenter Analysis of HHLWG

Background. Diffuse large B-cell lymphoma (DLBCL) is a heterogeneous non-Hodgkin’s lymphoma with great clinical challenge. Machine learning (ML) has attracted substantial attention in diagnosis, prognosis, and treatment of diseases. This study is aimed at exploring the prognostic factors of DLBCL by ML. Methods. In total, 1211 DLBCL patients were retrieved from Huaihai Lymphoma Working Group (HHLWG). The least absolute shrinkage and selection operator (LASSO) and random forest algorithm were used to identify prognostic factors for the overall survival (OS) rate of DLBCL among twenty-five variables. Receiver operating characteristic (ROC) curve and decision curve analysis (DCA) were utilized to compare the predictive performance and clinical effectiveness of the two models, respectively. Results. The median follow-up time was 43.4 months, and the 5-year OS was 58.5%. The LASSO model achieved an Area under the curve (AUC) of 75.8% for the prognosis of DLBCL, which was higher than that of the random forest model (AUC: 71.6%). DCA analysis also revealed that the LASSO model could augment net benefits and exhibited a wider range of threshold probabilities by risk stratification than the random forest model. In addition, multivariable analysis demonstrated that age, white blood cell count, hemoglobin, central nervous system involvement, gender, and Ann Arbor stage were independent prognostic factors for DLBCL. The LASSO model showed better discrimination of outcomes compared with the IPI and NCCN-IPI models and identified three groups of patients: low risk, high-intermediate risk, and high risk. Conclusions. The prognostic model of DLBCL based on the LASSO regression was more accurate than the random forest, IPI, and NCCN-IPI models.

The association between several autophagy-related genes and their prognostic values in hepatocellular carcinoma: a study on the foundation of TCGA, GEPIA and HPA databases

BACKGROUND

The purpose of this study was to investigate the relationship between the expression of autophagy-related genes and prognosis in hepatocellular carcinoma (HCC).

METHODS AND RESULTS

We selected three autophagy-related genes (ATG3, ATG7, and ATG9A) from gene expression data of liver cancer patients in The Cancer Genome Atlas (TCGA) database by Kaplan-Meier survival analysis, univariate and multivariate Cox regression analysis, and Gene Set Enrichment Analysis (GSEA). Human Protein Atlas (HPA) and Gene Expression Profiling Interactive Analysis (GEPIA) databases were applied to testify the credibility of our results. The expression levels of ATG3, ATG7, and ATG9A were verified by real-time quantitative PCR (RT-qPCR) in normal liver cells (L02) and three HCC cell lines (HepG2, Hep3b, and Li-7). Data analysis results from TCGA showed high ATG3, ATG7, ATG9A expression in HCC tumor tissues. Kaplan-Meier survival analysis showed that the survival rate of the high expression group of ATG3, ATG7, and ATG9A was all significantly lower than the low expression group. GSEA analysis showed that many signaling pathways (such as the regulation of autophagy, glycine serine and threonine metabolism, pathways in cancer, mitogen-activated protein kinase (MAPK) signaling pathway, mammalian target of rapamycin (mTOR) signaling pathway, as well as P53 signaling pathway) were differentially enriched in HCCs with ATG3, ATG7, and ATG9A expression. GEPIA and RT-qPCR also identified that the mRNA expression level of ATG3, ATG7, and ATG9A in normal liver cells were significantly lower than in HCC cells. High protein expression of ATG3, ATG7, and ATG9A was displayed in HCCs from the HPA database.

CONCLUSIONS

The ATG3, ATG7, ATG9A might be utilized as prognostic biomarkers for liver cancer.

Autophagy in gastrointestinal cancers

Gastrointestinal cancers are a group of cancers occurred in gastrointestinal tissues with high morbidity and mortality rate. Although numerous studies were conducted on the investigation of gastrointestinal cancers, the real mechanisms haven’t been discovered, and no effective methods of prevention and treatment of gastrointestinal cancers have been developed. Autophagy, a vital catabolic process in organisms, have been proven to participate in various mechanisms and signaling pathways, thus producing a regulatory effect on various diseases. The role of autophagy in gastrointestinal cancers remains unclear due to its high complexity. In this review, firstly, the biological features of autophagy will be introduced. Secondly, the role of autophagy in three popular gastrointestinal cancers, namely esophageal cancer, gastric cancer, and colorectal cancer will be described and discussed by reviewing the related literature. We aimed to bring novel insights in exploring the real mechanisms for gastrointestinal cancers and developing effective and efficient therapeutic methods to treat gastrointestinal cancers.

Construction of an immune-related gene signature for the prognosis and diagnosis of glioblastoma multiforme

Background

Increasing evidence has suggested that inflammation is related to tumorigenesis and tumor progression. However, the roles of immune-related genes in the occurrence, development, and prognosis of glioblastoma multiforme (GBM) remain to be studied.

Methods

The GBM-related RNA sequencing (RNA-seq), survival, and clinical data were acquired from The Cancer Genome Atlas (TCGA), Genotype-Tissue Expression (GTEx), Chinese Glioma Genome Atlas (CGGA), and Gene Expression Omnibus (GEO) databases. Immune-related genes were obtained from the Molecular Signatures Database (MSigDB). Differently expressed immune-related genes (DE-IRGs) between GBM and normal samples were identified. Prognostic genes associated with GBM were selected by Kaplan–Meier survival analysis, Least Absolute Shrinkage and Selection Operator (LASSO)-penalized Cox regression analysis, and multivariate Cox analysis. An immune-related gene signature was developed and validated in TCGA and CGGA databases separately. The Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analyses were performed to explore biological functions of the signature. The correlation between immune cell infiltration and the signature was analyzed by single-sample gene set enrichment analysis (ssGSEA), and the diagnostic value was investigated. The gene set enrichment analysis (GSEA) was performed to explore the potential function of the signature genes in GBM, and the protein–protein interaction (PPI) network was constructed.

Results

Three DE-IRGs [Pentraxin 3 (PTX3), TNFSF9, and bone morphogenetic protein 2 (BMP2)] were used to construct an immune-related gene signature. Receiver operating characteristic (ROC) curves and Cox analyses confirmed that the 3-gene-based prognostic signature was a good independent prognostic factor for GBM patients. We found that the signature was mainly involved in immune-related biological processes and pathways, and multiple immune cells were disordered between the high- and low-risk groups. GSEA suggested that PTX3 and TNFSF9 were mainly correlated with interleukin (IL)-17 signaling pathway, nuclear factor kappa B (NF-κB) signaling pathway, tumor necrosis factor (TNF) signaling pathway, and Toll-like receptor signaling pathway, and the PPI network indicated that they could interact directly or indirectly with inflammatory pathway proteins. Quantitative real-time PCR (qRT-PCR) indicated that the three genes were significantly different between target tissues.

Conclusion

The signature with three immune-related genes might be an independent prognostic factor for GBM patients and could be associated with the immune cell infiltration of GBM patients.

Prognostic score model based on six m6A-related autophagy genes for predicting survival in esophageal squamous cell carcinoma

BACKGROUND

Prognostic signatures based on autophagy genes have been proposed for esophageal squamous cell carcinoma (ESCC). Autophagy genes are closely associated with m6A genes. Our purpose is to identify m6A-related autophagy genes in ESCC and develop a survival prediction model.

METHODS

Differential expression analyses for m6A genes and autophagy genes were performed based on TCGA and HADd databases followed by constructing a co-expression network. Uni-variable Cox regression analysis was performed for m6A-related autophagy genes. Using the optimal combination of feature genes by LASSO Cox regression model, a prognostic score (PS) model was developed and subsequently validated in an independent dataset.

RESULTS

The differential expression of 13 m6A genes and 107 autophagy genes was observed between ESCC and normal samples. The co-expression network contained 13 m6A genes and 96 autophagy genes. Of the 12 m6A-related autophagy genes that were significantly related to survival, DAPK2, DIRAS3, EIF2AK3, ITPR1, MAP1LC3C, and TP53 were used to construct a PS model, which split the training set into two risk groups with significant different survival ratios (p = 0.015, 1-year, 3-year, and 5-year AUC = 0.873, 0.840, and 0.829). Consistent results of GSE53625 dataset confirmed predictive ability of the model (p = 0.024, 1-year, 3-year, and 5-year AUC = 0.793, 0.751, and 0.744). The six-gene PS score was an independent prognostic factor from clinical factors (HR, 2.362; 95% CI, 1.390-7.064; p-value = 0.012).

CONCLUSION

Our study recommends 6 m6A-related autophagy genes as promising prognostic biomarkers and develops a PS model to predict survival in ESCC.

The role of autophagy in initiation, progression, TME modification, diagnosis, and treatment of esophageal cancers

Autophagy is a highly conserved metabolic process with a cytoprotective function. Autophagy is involved in cancer, infection, immunity, and inflammation and may be a potential therapeutic target. Increasing evidence has revealed that autophagy has primary implications for esophageal cancer, including its initiation, progression, tumor microenvironment (TME) modification, diagnosis, and treatment. Notably, autophagy displayed excellent application potential in radiotherapy combined with immunotherapy. Radiotherapy combined with immunotherapy is a new potential therapeutic strategy for cancers, including esophageal cancer. Autophagy modulators can work as adjuvant enhancers in radiotherapy or immunotherapy of cancers. This review highlights the most recent data related to the role of autophagy regulation in esophageal cancer.

Endoplasmic Reticulum Stress-Related Four-Biomarker Risk Classifier for Survival Evaluation in Esophageal Cancer

Purpose

Esophageal cancer (EC) is a lethal digestive tumor worldwide with a dismal clinical outcome. Endoplasmic reticulum (ER) stress poses essential implications for a variety of tumor malignant behaviors. Here, we set up an ER stress-based risk classifier for assessing patient outcome and exploiting robust targets for medical decision-making of EC cases.

Methods

340 EC cases with transcriptome and survival data from two independent public datasets (TCGA and GEO) were recruited for this project. Cox regression analyses were employed to create a risk classifier based on ER stress-related genes (ERGs) which were strongly linked to EC cases' outcomes. Then, we detected and confirmed the predictive ability of our proposed classifier via a host of statistical methods, including survival analysis and ROC method. In addition, immune-associated algorithm was implemented to analyze the immune activity of EC samples.

Results

Four EGRs (BCAP31, HSPD1, PDHA1, and UBE2D1) were selected to build an EGR-related classifier (ERC). This classifier could distinguish the patients into different risky subgroups. The remarkable differences in patient outcome between the two groups were observed, and similar results were also confirmed in GEO cohort. In terms of the immune analysis, the ERC could forecast the infiltration level of immunocytes, such as Tregs and NK cells.

Conclusion

We created a four-ERG risk classifier which displays the powerful capability of survival evaluation for EC cases.

Feature screening for survival trait with application to TCGA high-dimensional genomic data

Background

In high-dimensional survival genomic data, identifying cancer-related genes is a challenging and important subject in the field of bioinformatics. In recent years, many feature screening approaches for survival outcomes with high-dimensional survival genomic data have been developed; however, few studies have systematically compared these methods. The primary purpose of this article is to conduct a series of simulation studies for systematic comparison; the second purpose of this article is to use these feature screening methods to further establish a more accurate prediction model for patient survival based on the survival genomic datasets of The Cancer Genome Atlas (TCGA).

Results

Simulation studies prove that network-adjusted feature screening measurement performs well and outperforms existing popular univariate independent feature screening methods. In the application of real data, we show that the proposed network-adjusted feature screening approach leads to more accurate survival prediction than alternative methods that do not account for gene-gene dependency information. We also use TCGA clinical survival genetic data to identify biomarkers associated with clinical survival outcomes in patients with various cancers including esophageal, pancreatic, head and neck squamous cell, lung, and breast invasive carcinomas.

Conclusions

These applications reveal advantages of the new proposed network-adjusted feature selection method over alternative methods that do not consider gene-gene dependency information. We also identify cancer-related genes that are almost detected in the literature. As a result, the network-based screening method is reliable and credible.

Development and validation of a survival model for esophageal adenocarcinoma based on autophagy-associated genes

Autophagy is a highly conserved catabolic process which has been implicated in esophageal adenocarcinoma (EAC). We sought to investigate the biological functions and prognostic value of autophagy-related genes (ARGs) in EAC. A total of 21 differentially expressed ARGs were identified between EAC and normal samples. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis were then applied for the differentially expressed ARGs in EAC, and the protein-protein interaction (PPI) network was established. Cox survival analysis and Lasso regression analysis were performed to establish a prognostic prediction model based on nine overall survival (OS)-related ARGs (CAPN1, GOPC, TBK1, SIRT1, ARSA, BNIP1, ERBB2, NRG2, PINK1). The 9-gene prognostic signature significantly stratified patient outcomes in The Cancer Genome of Atlas (TCGA)-EAC cohort and was considered as an independently prognostic predictor for EAC patients. Moreover, Gene set enrichment analysis (GSEA) analyses revealed several important cellular processes and signaling pathways correlated with the high-risk group in EAC. This prognostic prediction model was confirmed in an independent validation cohort (GSE13898) from The Gene Expression Omnibus (GEO) database. We also developed a nomogram with a concordance index of 0.78 to predict the survival possibility of EAC patients by integrating the risk signature and clinicopathological features. The calibration curves substantiated favorable concordance between actual observation and nomogram prediction. Last but not least, Erb-B2 Receptor Tyrosine Kinase 2 (ERBB2), a member of the prognostic gene signature, was identified as a potential therapeutic target for EAC patients. To sum up, we established and verified a novel prognostic prediction model based on ARGs which could optimize the individualized survival prediction in EAC.

Identification and Validation of Pyroptosis-Related Gene Signature to Predict Prognosis and Reveal Immune Infiltration in Hepatocellular Carcinoma

Background: Hepatocellular carcinoma (HCC) is characterized by a poor prognosis and accounts for the fourth common cause of cancer-related deaths. Recently, pyroptosis has been revealed to be involved in the progression of multiple cancers. However, the role of pyroptosis in the HCC prognosis remains elusive.

Methods: The clinical information and RNA-seq data of the HCC patients were collected from the TCGA-LIHC datasets, and the differential pyroptosis-related genes (PRG) were firstly explored. The univariate Cox regression and consensus clustering were applied to recognize the HCC subtypes. The prognostic PRGs were then submitted to the LASSO regression analysis to build a prognostic model in the TCGA training cohort. We further evaluated the predictive model in the TCGA test cohort and ICGC validation cohort (LIRI-JP). The accuracy of prediction was validated using the Kaplan—Meier (K-M) and receiver operating characteristic (ROC) analyses. The single-sample gene set enrichment analysis (ssGSEA) was used to determine the differential immune cell infiltrations and related pathways. Finally, the expression of the prognostic genes was validated by qRT-PCR in vivo and in vitro.

Results: We identified a total of 26 differential PRGs, among which three PRGs comprising GSDME, GPX4, and SCAF11 were subsequently chosen for constructing a prognostic model. This model significantly distinguished the HCC patients with different survival years in the TCGA training, test, and ICGC validation cohorts. The risk score of this model was confirmed as an independent prognostic factor. A nomogram was generated indicating the survival years for each HCC patient. The ssGSEA demonstrated several tumor-infiltrating immune cells to be remarkably associated with the risk scores. The qRT-PCR results also showed the apparent dysregulation of PRGs in HCC. Finally, the drug sensitivity was analyzed, indicating that Lenvatinib might impact the progression of HCC via targeting GSDME, which was also validated in human Huh7 cells.

Conclusion: The PRG signature comprised of GSDME, GPX4, and SCAF11 can serve as an independent prognostic factor for HCC patients, which would provide further evidence for more clinical and functional studies.

Autophagy and Digestive Disorders: Advances in Understanding and Therapeutic Approaches

The gastrointestinal (GI) tract is a series of hollow organs that is responsible for the digestion and absorption of ingested foods and the excretion of waste. Any changes in the GI tract can lead to GI disorders. GI disorders are highly prevalent in the population and account for substantial morbidity, mortality, and healthcare utilization. GI disorders can be functional, or organic with structural changes. Functional GI disorders include functional dyspepsia and irritable bowel syndrome. Organic GI disorders include inflammation of the GI tract due to chronic infection, drugs, trauma, and other causes. Recent studies have highlighted a new explanatory mechanism for GI disorders. It has been suggested that autophagy, an intracellular homeostatic mechanism, also plays an important role in the pathogenesis of GI disorders. Autophagy has three primary forms: macroautophagy, microautophagy, and chaperone-mediated autophagy. It may affect intestinal homeostasis, host defense against intestinal pathogens, regulation of the gut microbiota, and innate and adaptive immunity. Drugs targeting autophagy could, therefore, have therapeutic potential for treating GI disorders. In this review, we provide an overview of current understanding regarding the evidence for autophagy in GI diseases and updates on potential treatments, including drugs and complementary and alternative medicines.

Identification and Validation of Autophagy-Related Gene Nomograms to Predict the Prognostic Value of Patients with Cervical Cancer

Autophagy is a process of engulfing one's own cytoplasmic proteins or organelles and coating them into vesicles, fusing with lysosomes to form autophagic lysosomes, and degrading the contents it encapsulates. Increasing studies have shown that autophagy disorders are closely related to the occurrence of tumors. However, the prognostic role of autophagy genes in cervical cancer is still unclear. In this study, we constructed risk signatures of autophagy-related genes (ARGs) to predict the prognosis of cervical cancer. The expression profiles and clinical information of autophagy gene sets were downloaded from TCGA and GSE52903 queues as training and validation sets. The normal cervical tissue expression profile data from the UCSC XENA website (obtained from GTEx) were used as a supplement to the TCGA normal cervical tissue. Univariate COX regression analysis of 17 different autophagy genes was performed with the consensus approach. Tumor samples from TCGA were divided into six subtypes, and the clinical traits of the six subtypes had different distributions. Further absolute shrinkage and selection operator (LASSO) and multivariable COX regression yielded an autophagy genetic risk model consisting of eight genes. In the training set, the survival rate of the high-risk group was lower than that of the low-risk group (p < 0.0001). In the validation set, the AUC area of the receiver operating characteristic (ROC) curve was 0.772 for the training set and 0.889 for the verification set. We found that high and low risk scores were closely related to TNM stage (p < 0.05). The nomogram shows that the risk score combined with other indicators, such as G, T, M, and N, better predicts 1-, 3-, and 5-year survival rates. Decline curve analysis (DCA) shows that the risk model combined with other indicators produces better clinical efficacy. Immune cells with an enrichment score of 28 showed statistically significant differences related to high and low risk. GSEA enrichment analysis showed the main enrichment being in KRAS activation, genes defining epithelial and mesenchymal transition (EMT), raised in response to the low oxygen level (hypoxia) gene and NF-kB in response to TNF. These pathways are closely related to the occurrence of tumors. Our constructed autophagy risk signature may be a prognostic tool for cervical cancer.