肝细胞癌患者自噬相关基因的预后作用

构建由自噬相关基因组成的预后模型,预测肝细胞癌(HCC)患者的生存预后情况,为其个性化诊疗和临床研究提供依据.利用TCGA数据库中HCC的测序信息与人类自噬数据库联合,筛选差异表达的自噬相关基因,对其进行GO富集与KEGG通路分析;通过单因素与多因素Cox分析筛选与患者生存预后明显相关的风险基因,构建预后风险评分模型;...     

基于转录组测序的Cd胁迫下芹菜细胞自噬相关基因的筛选

细胞自噬是植物逆境应答过程中最常见的保护机制之一。动物中,自噬相关基因抵御镉(Cd)毒害的功能研究较清楚,但植物却知之甚少。文中以芹菜品种‘皇后’为试材,采用外源Cd(终浓度为0、2、4、8mg/L)添加营养液水培处理,利用转录组测序(RNA-seq)技术筛选细胞自噬相关差异基因并进行q RT-PCR验证。结果表明Cd胁迫对芹菜植株产生了明显的毒害作用,并与浓度间产生了量效关系。在筛选的8个差异表达的自噬相关基因中,ATG8a、ATG8f、ATG13、AMPK-1、AMPK-2基因随Cd浓度升高表达上调,ATG12、VPS30和VPS34则先上调后下降,说明自噬相关基因可能通过表达上调增加了自噬小体结构以抵御Cd毒性作用;而高浓度Cd(8mg/L)可能超出芹菜的耐受范围,导致多个自噬基因又出现表达下调趋势。以上结果有助于后期自噬相关基因的功能研究,为进一步探讨芹菜对Cd胁迫的耐性机制提供参考依据。     

Identification of novel prognostic indicators for oral squamous cell carcinoma based on proteomics and metabolomics

BackgroundThe low 5-year survival rate of oral squamous cell carcinoma (OSCC) suggests that new prognostic indicators need to be identified to aid the clinical management of patients.MethodsSaliva samples from OSCC patients and healthy controls were collected for proteomic and metabolomic sequencing. Gene expressed profiling was downloaded from TCGA and GEO databases. After the differential analysis, proteins with a significant impact on the prognosis of OSCC patients were screened. Correlation analysis was performed with metabolites and core proteins were identified. Cox regression analysis was utilized to stratify OSCC samples based on core proteins. The prognostic predictive ability of the core protein was then evaluated. Differences in infiltration of immune cells between the different strata were identified.ResultsThere were 678 differentially expressed proteins (DEPs), 94 intersected DEPs among them by intersecting with differentially expressed genes in TCGA and GSE30784 dataset. Seven core proteins were identified that significantly affected OSCC patient survival and strongly correlated with differential metabolites (R² > 0.8). The samples were divided into high- and low-risk groups according to median risk score. The risk score and core proteins were well prognostic factor in OSCC patients. Genes in high-risk group were enriched in Notch signaling pathway, epithelial mesenchymal transition (EMT), and angiogenesis. Core proteins were strongly associated with the immune status of OSCC patients.ConclusionsThe results established a 7-protein signatures with the hope of early detection and the capacity for risk assessment of OSCC patient prognosis. Further providing more potential targets for the treatment of OSCC.     

A risk signature with four autophagy-related genes for predicting survival of glioblastoma multiforme

Glioblastoma multiforme (GBM) is a devastating brain tumour without effective treatment. Recent studies have shown that autophagy is a promising therapeutic strategy for GBM. Therefore, it is necessary to identify novel biomarkers associated with autophagy in GBM. In this study, we downloaded autophagy-related genes from Human Autophagy Database (HADb) and Gene Set Enrichment Analysis (GSEA) website. Least absolute shrinkage and selection operator (LASSO) regression and multivariate Cox regression analysis were performed to identify genes for constructing a risk signature. A nomogram was developed by integrating the risk signature with clinicopathological factors. Time-dependent receiver operating characteristic (ROC) curve and calibration plot were used to evaluate the efficiency of the prognostic model. Finally, four autophagy-related genes (DIRAS3, LGALS8, MAPK8 and STAM) were identified and were used for constructing a risk signature, which proved to be an independent risk factor for GBM patients. Furthermore, a nomogram was developed based on the risk signature and clinicopathological factors (IDH1 status, age and history of radiotherapy or chemotherapy). ROC curve and calibration plot suggested the nomogram could accurately predict 1-, 3- and 5-year survival rate of GBM patients. For function analysis, the risk signature was associated with apoptosis, necrosis, immunity, inflammation response and MAPK signalling pathway. In conclusion, the risk signature with 4 autophagy-related genes could serve as an independent prognostic factor for GBM patients. Moreover, we developed a nomogram based on the risk signature and clinical traits which was validated to perform better for predicting 1-, 3- and 5-year survival rate of GBM.     

Prognostic significance of autophagy related genes in estrogen receptor positive tamoxifen treated breast cancer

Resistance to Tamoxifen constitutes a major therapeutic challenge in treating hormone sensitive breast cancer. The induction of autophagy has been shown to be involved as one of the mechanism responsible for Tamoxifen resistance. Autophagy related gene (ATG) members are the regulators and effectors of Macroautophagy process in the cellular systems. In this study, we evaluated the prognostic significance of ATGs in Tamoxifen treated breast cancer. The "Kaplan- Meier plotter" database was utilized to analyze the relevance and significance of ATGs mRNA expression to Relapse Free Survival in breast cancer patients. We used the data of patients who are Estrogen receptor positive and are treated with Tamoxifen. Hazard ratio and log-rank p-value were calculated using KM survival plots for various ATGs. Overexpressed ATG3, ATG 5, ATG 8B and PIK3R4 resulted in a poor prognosis. A gene signature of these ATGs predicts deteriorated RFS (p-value=8.3e-05 and HR=1.84 (1.35-2.51) and Distant Metastasis Free Survival (p value = 0.0027 and HR=2.03 (1.27-3.26). We report the distinct prognostic values of ATGs in patients of breast cancer treated with Tamoxifen. Thus, better understandings of the induction of autophagy pathway may potentially form the basis for use of autophagy inhibitors in the Tamoxifen treated breast cancer.     

Identification and validation of autophagy-related prognostic signature for head and neck squamous cell carcinoma

BackgroundMany studies have demonstrated that autophagy plays a significant role in regulating tumor growth and progression. However, the effect of autophagy-related genes (ARGs) on the prognosis have rarely been analyzed in head and neck squamous cell carcinoma (HNSCC).MethodsWe obtained differentially expressed ARGs from HNSCC mRNA data in The Cancer Genome Atlas (TCGA) database. And then we performed gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analyses to explore the autophagy-related biological functions. The overall survival (OS)-related and disease specific survival (DSS)-related ARGs were identified by univariate Cox regression analyses. With these genes, we established OS-related and DSS-related risk signature by LASSO regression method, respectively. We validated the reliability of the risk signature with receiver operating characteristic (ROC) analysis, Kaplan-Meier survival curves, clinical correlation analysis, and nomogram. Then we analyzed relationships between risk signature and immune cell infiltration.ResultsWe established the prognostic signatures based on 14 ARGs for OS and 12 ARGs for DSS. The ROC curves, survival analysis, and nomogram validated the predictive accuracy of the models. Clinic correlation analysis showed that the risk group was closely related to Stage, pathological T stage, pathological N stage and human papilloma virus (HPV) subtype. Cox regression demonstrated that the risk score was an independent predictor for the prognosis of HNSCC patients. Furthermore, patients in low-risk score group exhibited higher immunescore and distinct immune cell infiltration than high-risk score group. And we further analysis revealed that the copy number alterations (CNAs) of ARGs-based signature affected the abundance of tumor-infiltrating immune cells.ConclusionIn this study, we identified novel autophagy-related signature for the prediction of OS and DSS in patients with HNSCC. Meanwhile, our study provides a novel sight to understand the role of autophagy and elucidate the important role of autophagy in tumor immune microenvironment (TIME) of HNSCC.     

ATIC inhibits autophagy in hepatocellular cancer through the AKT/FOXO3 pathway and serves as a prognostic signature for modeling patient survival

Background: Autophagy regulates many cell functions related to cancer, ranging from cell proliferation and angiogenesis to metabolism. Due to the close relationship between autophagy and tumors, we investigated the predictive value of autophagy-related genes.Methods: Data from patients with hepatocellular carcinoma were obtained from The Cancer Genome Atlas (TCGA) and the International Cancer Genome Consortium (ICGC) databases. A regression analysis of differentially expressed genes was performed. Based on a prognostic model, patients were divided into a high-risk or low-risk group. Kaplan-Meier survival analyses of patients were conducted. The immune landscapes, as determined using single-sample gene set enrichment analysis (ssGSEA), exhibited different patterns in the two groups. The prognostic model was verified using the ICGC database and clinical data from patients collected at Zhongnan Hospital. Based on the results of multivariate Cox regression analysis, 5-aminoimidazole-4-carboxamide ribonucleotide formyltransferase/inosine monophosphate (IMP) cyclohydrolase (ATIC) had the largest hazard ratio, and thus we studied the effect of ATIC on autophagy and tumor progression by performing in vitro and in vivo experiments.Results: Fifty-eight autophagy-related genes were differentially expressed (false discovery rate (FDR)<0.05, log₂ fold change (logFC)>1); 23 genes were related to the prognosis of patients. A prognostic model based on 12 genes (ATG10, ATIC, BIRC5, CAPN10, FKBP1A, GAPDH, HDAC1, PRKCD, RHEB, SPNS1, SQSTM1 and TMEM74) was constructed. A significant difference in survival rate was observed between the high-risk group and low-risk group distinguished by the model (P<0.001). The model had good predictive power (area under the curve (AUC)>0.7). Risk-related genes were related to the terms type II IFN response, MHC class I (P<0.001) and HLA (P<0.05). ATIC was confirmed to inhibit autophagy and promote the proliferation, invasion and metastasis of liver cancer cells through the AKT/Forkhead box subgroup O3 (FOXO3) signaling pathway in vitro and in vivo.Conclusions: The prediction model effectively predicts the survival time of patients with liver cancer. The risk score reflects the immune cell features and immune status of patients. ATIC inhibits autophagy and promotes the progression of liver cancer through the AKT/FOXO3 signaling pathway.     

A novel autophagy-related lncRNA prognostic risk model for breast cancer

Long non-coding RNAs (lncRNAs) are well known as crucial regulators to breast cancer development and are implicated in controlling autophagy. LncRNAs are also emerging as valuable prognostic factors for breast cancer patients. It is critical to identify autophagy-related lncRNAs with prognostic value in breast cancer. In this study, we identified autophagy-related lncRNAs in breast cancer by constructing a co-expression network of autophagy-related mRNAs-lncRNAs from The Cancer Genome Atlas (TCGA). We evaluated the prognostic value of these autophagy-related lncRNAs by univariate and multivariate Cox proportional hazards analyses and eventually obtained a prognostic risk model consisting of 11 autophagy-related lncRNAs (U62317.4, LINC01016, LINC02166, C6orf99, LINC00992, BAIAP2-DT, AC245297.3, AC090912.1, Z68871.1, LINC00578 and LINC01871). The risk model was further validated as a novel independent prognostic factor for breast cancer patients based on the calculated risk score by Kaplan-Meier analysis, univariate and multivariate Cox regression analyses and time-dependent receiver operating characteristic (ROC) curve analysis. Moreover, based on the risk model, the low-risk and high-risk groups displayed different autophagy and oncogenic statues by principal component analysis (PCA) and Gene Set Enrichment Analysis (GSEA) functional annotation. Taken together, these findings suggested that the risk model of the 11 autophagy-related lncRNAs has significant prognostic value for breast cancer and might be autophagy-related therapeutic targets in clinical practice.     

Autophagy is a regulator of TRAIL-induced apoptosis in NSCLC A549 cells

Autophagy, a catabolic process by which cytoplasmic components are degraded in lysosomes, plays an important role in the maintenance of cellular homeostasis. Dysregulation of autophagy is associated with several diseases. However, few studies have addressed the role of autophagy in the lung, and its role in lung diseases remains unclear. In the present study, we examined the effect of tumor necrosis factor-related apoptosis-inducing ligand (TRAIL) on autophagy in A549 cells and explored the underlying mechanisms. We showed that TRAIL promoted autophagosome formation, as detected by the levels of LC3-II, and its effect on promoting autophagy was dependent on the expression of the autophagy related genes (ATGs) Atg5, Atg7, and beclin-1. TRAIL-induced ATG expression was attenuated by JNK silencing or treatment with the JNK inhibitor SP600125, indicating the involvement of the JNK pathway. Crosstalk between autophagy and apoptosis was demonstrated by silencing the autophagy related genes Atg5, Atg7, and beclin-1, and the dependence of TRAIL-induced apoptosis on autophagy-related gene expression. Taken together, our results indicate that TRAIL promotes autophagy in A549 cells via a mechanism involving the modulation of ATG expression through the JNK pathway. Inhibition of autophagy enhanced TRAIL-induced cell proliferative inhibition and apoptosis in A549 cells.     

Identification of a six-gene prognostic signature for oral squamous cell carcinoma

Oral squamous cell carcinoma (OSCC) is one of the most common types of malignancies worldwide, and its morbidity and mortality have increased in the near term. Consequently, the purpose of the present study was to identify the notable differentially expressed genes (DEGs) involved in their pathogenesis to obtain new biomarkers or potential therapeutic targets for OSCC. The gene expression profiles of the microarray datasets GSE85195, GSE23558, and GSE10121 were obtained from the Gene Expression Omnibus (GEO) database. After screening the DEGs in each GEO dataset, 249 DEGs in OSCC tissues were obtained. Kyoto Encyclopedia of Genes and Genomes and Gene Ontology pathway enrichment analysis was employed to explore the biological functions and pathways of the above DEGs. A protein–protein interaction network was constructed to obtain a central gene. The corresponding total survival information was analyzed in patients with oral cancer from The Cancer Genome Atlas (TCGA). A total of six candidate genes (CXCL10, OAS2, IFIT1, CCL5, LRRK2, and PLAUR) closely related to the survival rate of patients with oral cancer were identified, and expression verification and overall survival analysis of six genes were performed based on TCGA database. Time-dependent receiver operating characteristic curve analysis yields predictive accuracy of the patient's overall survival. At the same time, the six genes were further verified by quantitative real-time polymerase chain reaction using samples obtained from the patients recruited to the present study. In conclusion, the present study identified the prognostic signature of six genes in OSCC for the first time via comprehensive bioinformatics analysis, which could become potential prognostic markers for OCSS and may provide potential therapeutic targets for tumors.     

IDUA,NDST1, SAP30L,CRYBA4, and SI as novel prognostic signatures clear cell renal cell carcinoma

Carcinoma of the kidney is one of the most prevalent carcinoma worldwide. The majority types of carcinoma are clear cell renal cell carcinoma (CCRCC), which consist more than 80% of the cases. As a genetically diverse disease, identification of prognosis-related genes has utmost importance in the early diagnosis and prognosis of the CCRCC. In this study, we performed gene expression profiling to identify prognosis-related genes for CCRCC. In addition, we developed and validated a gene signature-based risk score to comprehensively assess the prognostic function of differentially expressed genes. Furthermore, we performed a ROC analysis to identify the optimal cut-off point for classification risk level of the patients. Univariate Cox regression models were used to assess the association between differentially expressed genes in relation to the prognosis of patients with different stages of CCRCC. Five genes were identified significantly differentially expressed in CCRCC and associated with their survival time, namely: IDUA, NDST1, SAP30L, CRYBA4, and SI. A 5-gene signature-based risk score was developed based on the Cox coefficient of the individual genes. The prognostic value of this risk score was validated in an internal testing data set. In summary, a gene-based risk score was identified and validated, which can predict CCRCC patient survival. The potential functions of this gene expression signature and individual differentially expressed gene as prognostic targets of CCRCC were revealed by this study. Furthermore, these findings may have important implications in the understanding of the potential therapeutic method for the CCRCC patients.     

基于转录组测序探究ATP7B基因缺陷小鼠铜累积诱导肝细胞自噬的相关机制*

目的：分析ATP7B基因缺陷(Wilson's disease,WD)小鼠肝脏组织中自噬相关基因的表达和自噬相关蛋白的相互作用方式,探讨铜累积诱导肝内自噬活化的可能机制。方法：对4周龄和12周龄WD小鼠肝组织进行铜含量检测和转录组测序,对差异基因进行GO和KEGG富集分析,筛选自噬相关差异基因做qRT-PCR和Western blot验证,采用GeneMANIA数据库构建自噬相关差异蛋白的互作网络(PPI)并进行功能注释分析,抑制自噬相关蛋白的表达分析其对自噬的影响。结果：与野生型小鼠相比,WD小鼠肝铜含量显著升高,铜累积导致基因表达模式改变;基于GO数据库统计自噬相关差异基因数目,4周龄和12周龄分别有8个、51个,基于KEGG数据库统计,4周龄和12周龄分别有5个、19个;筛选Ulk1、Ddit4、Plk3等9个基因进行qRT-PCR,定量结果与测序结果表达趋势基本一致;其编码的蛋白质通过共表达、共定位等方式互相作用;Western blot结果显示铜累积导致Ulk1、Plk3、Park2蛋白表达显著增加和细胞自噬发生,抑制Ulk1、Plk3、Park2的蛋白质表达可显著下调细胞自噬水平。结论：WD不同阶段的铜累积可调节肝脏多个自噬相关基因的表达,通过其编码的自噬相关蛋白的互相作用共同诱导肝脏自噬活化以缓解肝损伤。     

Autophagy-related circRNA evaluation reveals hsa_circ_0001747 as a potential favorable prognostic factor for biochemical recurrence in patients with prostate cancer

Prostate cancer (PCa) is a common high-incidence malignancy in men, some of whom develop biochemical recurrence (BCR) in the advanced stage. However, there are currently no accurate prognostic indicators of BCR in PCa. The aim of our study was to identify an autophagy-related circular RNA prognostic factor of BCR for patients with PCa. In this study, immunochemistry revealed that the classic autophagy marker MAP1LC3B was positively correlated with Gleason score. Least absolute shrinkage and selector operator regression were conducted to develop a novel prognostic model with tenfold cross-validation and an L1 penalty. Five autophagy-related circRNA signatures were included in the prognostic model. Patients with PCa were ultimately divided into high- and low-risk groups, based on the median risk score. Patients with PCa, who had a high risk score, were more likely to develop BCR in a shorter period of time. Univariate and multivariate Cox regression analyses demonstrated that the risk score was an independent variable for predicting BCR in PCa. In addition, a prognostic nomogram integrated with the risk score and numerous clinicopathological parameters was developed to accurately predict 3- and 5-year BCR of patients with PCa. Finally, the hsa_circ_0001747 signature was selected for further experimental verification in vitro and in vivo, which showed that downregulated hsa_circ_0001747 might facilitate PCa via augmenting autophagy. Our findings indicate that the autophagy-related circRNA signature hsa_circ_0001747 may serve as a promising indicator for BCR prediction in patients with PCa.Subject terms: Tumour biomarkers, Macroautophagy     

Editor's choice: A comprehensive,genome-wide analysis of autophagy-related genes identified in tobacco suggests a central role of autophagy in plant response to various environmental cues

Autophagy is an evolutionarily conserved mechanism in both animals and plants, which has been shown to be involved in various essential developmental processes in plants. Nicotiana tabacum is considered to be an ideal model plant and has been widely used for the study of the roles of autophagy in the processes of plant development and in the response to various stresses. However, only a few autophagy-related genes (ATGs) have been identified in tobacco up to now. Here, we identified 30 ATGs belonging to 16 different groups in tobacco through a genome-wide survey. Comprehensive expression profile analysis reveals an abroad expression pattern of these ATGs, which could be detected in all tissues tested under normal growth conditions. Our series tests further reveal that majority of ATGs are sensitive and responsive to different stresses including nutrient starvation, plant hormones, heavy metal and other abiotic stresses, suggesting a central role of autophagy, likely as an effector, in plant response to various environmental cues. This work offers a detailed survey of all ATGs in tobacco and also suggests manifold functions of autophagy in both normal plant growth and plant response to environmental stresses.     

Processing of ATG8s, ubiquitin-like proteins, and their deconjugation by ATG4s are essential for plant autophagy

Autophagy is an intracellular process for vacuolar degradation of cytoplasmic components. Thus far, plant autophagy has been studied primarily using morphological analyses. A recent genome-wide search revealed significant conservation among autophagy genes (ATGs) in yeast and plants. It has not been proved, however, that Arabidopsis thaliana ATG genes are required for plant autophagy. To evaluate this requirement, we examined the ubiquitination-like Atg8 lipidation system, whose component genes are all found in the Arabidopsis genome. In Arabidopsis, all nine ATG8 genes and two ATG4 genes were expressed ubiquitously and were induced further by nitrogen starvation. To establish a system monitoring autophagy in whole plants, we generated transgenic Arabidopsis expressing each green fluorescent protein-ATG8 fusion (GFP-ATG8). In wild-type plants, GFP-ATG8s were observed as ring shapes in the cytoplasm and were delivered to vacuolar lumens under nitrogen-starved conditions. By contrast, in a T-DNA insertion double mutant of the ATG4s (atg4a4b-1), autophagosomes were not observed, and the GFP-ATG8s were not delivered to the vacuole under nitrogen-starved conditions. In addition, we detected autophagic bodies in the vacuoles of wild-type roots but not in those of atg4a4b-1 in the presence of concanamycin A, a V-ATPase inhibitor. Biochemical analyses also provided evidence that autophagy in higher plants requires ATG proteins. The phenotypic analysis of atg4a4b-1 indicated that plant autophagy contributes to the development of a root system under conditions of nutrient limitation.     

BAG6/BAT3 modulates autophagy by affecting EP300/p300 intracellular localization

We recently reported that BAG6/BAT3 (BCL2-associated athanogene 6) is essential for basal and starvation-induced autophagy in E18.5 bag6^−/− mouse embryos and in mouse embryonic fibroblasts (MEFs) through the modulation of the EP300/p300-dependent acetylation of TRP53 and autophagy-related (ATG) proteins. We observed that BAG6 increases TRP53 acetylation during starvation and pro-autophagic TRP53-target gene expression. BAG6 also decreases the EP300 dependent-acetylation of ATG5, ATG7, and LC3-I, posttranslational modifications that inhibit autophagy. In addition, in the absence of BAG6 or when using a mutant of BAG6 exclusively located in the cytoplasm, autophagy is inhibited, ATG7 is hyperacetylated, TRP53 acetylation is abrogated, and EP300 accumulates in the cytoplasm indicating that BAG6 is involved in the regulation of the nuclear localization of EP300. We also reported that the interaction between BAG6 and EP300 occurs in the cytoplasm rather than the nucleus. Moreover, during starvation, EP300 is transported to the nucleus in a BAG6-dependent manner. We concluded that BAG6 regulates autophagy by controlling the localization of EP300 and its accessibility to nuclear (TRP53) and cytoplasmic (ATGs) substrates.