Prognostic value of an autophagy-related gene expression signature for endometrial cancer patients

doi:10.21203/rs.3.rs-24691/v3

Download PDF

Primary research

Prognostic value of an autophagy-related gene expression signature for endometrial cancer patients

https://doi.org/10.21203/rs.3.rs-24691/v3

This work is licensed under a CC BY 4.0 License

Journal Publication

published 13 Jul, 2020

Read the published version in Cancer Cell International →

You are reading this latest preprint version

Background: Autophagy is associated with cancer development. Autophagy-related genes play significant roles in endometrial cancer (EC), a major gynecological malignancy worldwide, but little was known about their value as prognostic markers. Here we evaluated the value of a prognostic signature based on autophagy-related genes for EC.

Methods: First, various autophagy-related genes were obtained via the Human Autophagy Database and their expression profiles were downloaded from The Cancer Genome Atlas. Second, key prognostic autophagy-related genes were identified via univariat, LASSO and multivariate Cox regression analyses. Finally, a risk score to predict the prognosis of EC was calculated and validated by using the test and the entire data sets. Besides, the key genes mRNA expression were validated using quantitative real-time PCR in clinical tissue samples.

Results: A total of 40 differentially expressed autophagy-related genes in EC were screened and five of them were prognosis-related (CDKN1B, DLC1, EIF4EBP1, ERBB2 and GRID1). A prognostic signature was constructed based on these five genes using the train set, which stratified EC patients into high-risk and low-risk groups (P<0.05). In terms of overall survival, the analyses of the test set and the entire set yielded consistent results (test set: p < 0.05; entire set: p < 0.05). Time-dependent ROC analysis suggested that the risk score predicted EC prognosis accurately and independently (0.674 at 1 year, 0.712 at 3 years and 0.659 at 5 years). A nomogram with clinical utility was built. Patients in the high-risk group displayed distinct mutation signatures compared with those in the low-risk group. For clinical sample validation, we found that EIF4EBP1and ERBB2 had higher level in EC than that in normal tissues while CDKN1B, DLC1 and GRID1 had lower level, which was consistent with the results predicted.

Conclusions: Based on five autophagy-related genes (CDKN1B, DLC1, EIF4EBP1, ERBB2 and GRID1), our model can independently predict the OS of EC patients by combining molecular signature and clinical characteristics.

Cancer Biology

Oncology

autophagy-related genes

endometrial cancer

prognosis

As an evolutionarily ancient and highly conserved biological behavior, autophagy plays a cytoprotective role in eukaryotic cells through degrading unnecessary or dysfunctional organelles under the condition of hypoxia, starvation, nutrition deficiency or high PH [1]. However, excessive self-degradation can be poisonous [2]. Accordingly, autophagic dysfunction involved in numerous human pathologies, such as liver and heart disease, neurodegeneration, muscle disease, and cancer [3]. However, the mechanisms of autophagy in cancer remain stagnantly understood. A clearer eyesight into these mechanisms may facilitate the invention of autophagy-focused therapeutic interventions.

Endometrial cancer (EC) is one of the most common malignant diseases in women worldwide, with an estimated 61,880 newly-diagnosed cases and 12,160 deaths in the United States in 2019 [4]. Primarily treated with surgical resection, the prognosis of women with early-stage EC is excellent, with a 5-year survival rate of over 70% [5]. However, 15-20% of cases still show post-operative vaginal or pelvic EC relapse, often metastasizing to distant sites [6]. Despite the advance in surgical techniques, chemotherapy, or radiotherapy, data show that the EC mortality of EC has bypassed EC incidence [7]. To further improve the outcomes of EC treatments, physicians need to identify high-risk patients and tailor precise treatment.

Recent studies revealed associations between autophagy and pathophysiological processes in EC. The so-called "stone-like" structure (SLS) has emerged as the hallmark of autophagic activity. SLSs are detected exclusively in EC and rarely in atypical hyperplasias, suggesting that autophagy is more active in EC cells than in normal or hyperplastic endometrial cells [8]. Targeting autophagy enhances sorafenib cytotoxicity and inhibits tumor growth and pulmonary metastasis, which can help us gain insights into the unopposed resistance of advanced EC to sorafenib and provide a possibility to develop a new therapeutic intervention for recurrent EC [9]. However, little literature is available regarding the prediction of EC prognosis based on its molecular mechanisms or biological behavior.

In the present study, we examined the relevance between autophagy-related genes (ARGs) expression profiles and clinical outcomes in 552 EC patients from The Cancer Genome Atlas (TCGA). Then, using ARGs as an independent risk factor, a prognosis-prediction model was established for EC overall survival (OS). Results from this study could provide an autophagy-targeted strategy for predicting and monitoring the prognosis of patients with EC.

Collection of data

The Human Autophagy Database (http://www.autophagy.lu/index.html) is a public repository containing information about the human genes involved in autophagy. Various ARGs were obtained via this database. RNA sequencing (RNA-seq) data of ARGs from each individual and their clinical information were downloaded from The Cancer Genome Atlas (TCGA) data portal (http://cancergenome.nih.gov/).

Identification of differentially expressed ARGs

The limma package in R software was used to search for differentially expressed ARGs between EC and non-tumor samples. Genes exhibiting at least 1-fold change and an FDR< 0.05 were regarded as obviously differentially expressed. Then, their most involved biological functions and pathways were revealed by a series of gene functional enrichment analyses, including Gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) at Enrichr (http://amp.pharm.mssm.edu/Enrichr/) [10].

Calculation of a risk score based on ARG model

ARGs expression profiles were downloaded from TCGA and normalized by [log2(data+1)]. The entire set was randomly separated into the train set and the test set in a ratio of 6 to 4. Using the gene expression profile of the train set, ARGs were submitted to univariate, LASSO Cox regression and multivariate analyses, to remove the genes which might not be an independent indicator in EC prognosis. Time-dependent receiver operating characteristic (ROC) analysis was used to investigate the prognosis-prediction accuracy and risk score of each ARG. The area under the curve (AUC) with different cutoffs was used to measure prognosis-prediction accuracy. Finally, several prognostic ARGs were screened out, and the coefficient of each prognostic gene in the risk model was assessed with corresponding multivariate analysis. EC patients were then divided into high- and low-risk groups by the cutoff medians, which was calculated according to the expression level of ARGs and the estimated regression coefficient. To investigate whether the risk score could be used as an independent predictor of OS in the TCGA cohort of EC patients, the univariate and multivariate Cox regression analyses were conducted. Clinicopathologic characteristics of EC patients were downloaded from TCGA database (Table 1). The risk score, age, clinical stage, grade, peritoneal wash, radiation therapy, surgical approach, cancer status, and histological type were used as covariates. The test set and the entire set were used to validate this signature and a nomogram was constructed. This nomogram was corrected by calibration curves.

Gene Set Enrichment analysis

Gene Set Enrichment Analysis (GSEA) was set up in TCGA series using a molecular signature database (MSigDB) c2.cp.kegg.v6.2.symbols.gmt gene sets [11]. The GSEA, visualized in Enrichment Map software and Cytoscape [12], was applied to determine whether the members of a given gene set were significantly correlated with our risk score.

Somatic mutation analysis

Package “maftools” was used to summarize, analyze, annotate and visualize mutation annotation format (MAF) files [13]. The T test was applied to identify the differentially mutated genes between high-risk and low-risk groups. The plotmaf Summary function was employed to plot the numbers and types of various MAF files. OncoPlot was used to plot the top ten most obviously mutated genes.

Preparation for human tissue samples

This research was approved by the Institutional Review Board of Nanjing Medical University and informed consent was obtained from all subjects. Tissue samples were collected from patients after surgery at the First Affiliated Hospital of Nanjing Medical University from 2017 to 2019 and then stored at -80°C, including fifteen EC tissues and fifteen normal endometrial tissues.

RNA isolation and quantitative real-time PCR ( qRT-PCR)

Total RNA from human tissues were extracted separately using Trizol reagent (Invitrogen, USA). The integrity of isolated RNA was assessed by NanoDrop 2000 Spectrophotometer (Thermo Scientific, Wilmington, DE, USA). RNA was reversely transcribed into cDNA using the high capacity reverse transcription kits (TaKaRa, Shiga, Japan) and qRT-PCR was performed in the LightCycler480II (Roche) using SYBR Green technology (Takara). The relative expression level was calculated by comparing the 2-ΔΔCt method. The primers for GAPDH, CDKN1B, DLC1, EIF4EBP1, ERBB2 and GRID1 were purchased from Tsingke (NanJing, China). The forward primer of GAPDH is CGCTCTCTGCTCCTCCTGTT. The reverse primer of GAPDH is CATGGGTGGAATCATATTGG. The forward primer of CDKN1B is AACGTGCGAGTGTCTAACGG. The reverse primer of CDKN1B is CCCTCTAGGGGTTTGTGATTCT. The forward primer of DLC1 is CCACGGACCTCCCATCTTC. The reverse primer of DLC1 is GCTGTGCATACTGGGGGAA. The forward primer of EIF4EBP1 is CTATGACCGGAAATTCCTGATGG. The reverse primer of EIF4EBP is CCCGCTTATCTTCTGGGCTA. The forward primer of ERBB2 is TGCAGGGAAACCTGGAACTC. The reverse primer of ERBB2 is ACAGGGGTGGTATTGTTCAGC. The forward primer of GRID1 is TGGCTGTGCATCTGCCAAT. The reverse primer of GRID1 is CGTAGAACATGACGAACTTCTGC. All procedures for qRT-PCR were performed according to the manufacturer’s protocol and all experiments were repeated three times.

Statistical analysis

Univariate and multivariate Cox regression analyses were used to select the survival-related ARGs and construct the risk score model. The survival curves were plotted by Kaplan–Meier (K–M) method and differences in the survival rates between high- and low-risk groups were evaluated by the log-rank test. ROC curve and AUC were plotted by SurvivalROC package in R. P-value less than 0.05 was considered statistically significant.

Identification of differently expressed ARGs

RNA-seq and clinical data of 552 EC tissue samples and 23 non-tumor samples were downloaded from TCGA. With |log2(Fold Change)|>1 and FDR <0.05, we finally obtained 22 up-regulated and 18 down-regulated ARGs between EC and non-tumor tissues (Figure 1A, 1B). Box plots were visualized to display the expression patterns of the 40 differentially expressed ARGs, including 22 up-regulated genes (APOL1, ATG4D, BAK1, BAX, BID, BIRC5, BNIP3, CASP3, CDKN2A, CTSB, DAPK2, EIF4EBP1, ERBB2, GAPDH, IFNG, IKBKE, P4HB, PARP1, PTK6, SERPINA1, TP63, WIPI1 ) and 18 down-regulated genes (BCL2, CALCOCO2, CDKN1B, DLC1, FOS, FOXO1, GABARAPL1, GRID1, HSPB8, ITPR1, MYC, NRG2, PPP1R15A, PRKAR1A, RAB33B, ST13, TUSC1, VAMP3) (Figure 1C).

Functional annotation of the differentially expressed ARGs

According to the results of Enrichr database, in the aspect of biological processes, the genes were mostly enriched in apoptotic process, intrinsic apoptotic signaling pathway and intrinsic apoptotic signaling pathway in response to endoplasmic reticulum stress. In the aspect of cellular components, the genes were mostly enriched in mitochondrial outer membrane, mitochondrion, and autophagosome. In the aspect of molecular function, the genes were mostly enriched in protein heterodimerization activity, protein serine/threonine kinase inhibitor activity, and protein phosphatase binding. Besides, we also found the differentially expressed ARGs were notably associated with apoptosis, pathways in cancer, autophagy, Kaposi sarcoma-associated herpesvirus infection and measles in the KEGG pathway enrichment analysis (Additional file 1).

Construction of a risk score formula in the train set

The relationships between the expression level of the 40 differentially expressed ARGs and overall survival (OS) were evaluated based on the data obtained from TCGA. Firstly, we integrated autophagy-related gene expression profiles and clinical follow-up information to screen out 516 EC samples. Then the entire set (516 samples) were randomly separated into the train set and the test set in a ratio of 6 to 4. Hence, the train set contained 312 samples and the test set contained 204 samples. Using the expression files of the train set, the univariate Cox regression, LASSO Cox regression and multivariate Cox regression analyses were performed (Additional file 2). A prognostic model using the expression of CDKN1B, DLC1, EIF4EBP1, ERBB2 and GRID1 was constructed (Table 2). The risk score of EC prognosis was quantified by the following formula: risk score= CDKN1B * (0.016941) + DLC1 * (-0.33079) + EIF4EBP1 * (0.003577) + ERBB2 * (0.001271) + GRID1 * (0.733501). It was noticed that among these five ARGs, just DLC1 had a negative coefficient. Totally, the expressions of CDKN1B, EIF4EBP1, ERBB2 and GRID1 were negatively related to the OS of EC patients, while that of DLC1 was positively related to OS.

Relationships between the risk score and OS in the train set

Based on these five ARGs, we could calculate the risk score for each patient. Survival analysis was performed and a dichotomous score was adopted. Patients in the train set were then stratified into high-risk (n=156) and low-risk groups (n=156) based on the median risk score. The risk score and survival status predicted by the prognostic model were displayed in Figure 2A-C. The survival rate of the patients in the low-risk group was significantly higher than that in the high-risk group (P= 4.351e-04) (Figure 2D). Time-dependent ROC analysis demonstrated that the prognostic accuracy of the risk score was 0.706 at 1 year, 0.742 at 3 years, and 0.655 at 5 years, indicative of its good performance in predicting EC prognosis (Figure 2E). Principal component analysis displayed a different distribution pattern of high and low risk according to 5 autophagy-related gene expression based on the train set (Figure 2F). Univariate and multivariate Cox regression analyses showed that the model can efficiently predict EC prognosis with ARG risk score as an independent indicator (Figure 5A, B).

Validation of the risk score in the test set

The test set was applied to validate the predictive power of the five ARGs. Based on the cut-offs in the train set, 204 samples in the test set were divided into low-risk (n=107) and high-risk groups (n=97) by using the same risk score formula. Figure 3A-C also showed the risk score and survival status. In consistent with the findings described above, low-risk patients had significantly longer overall survival than high-risk patients (P=4.398e-02) (Figure 3D). ROC curve analysis showed that the specificity and sensitivity were highest when the risk score was 0.627 at 1 year, 0.66 at 3 years, 0.67 at 5 years, according to the value of the area under the receiver operating characteristic curve (AUC) (Figure 3E). Principal component analysis displayed a different distribution pattern of high and low risk according to 5 autophagy-related gene expression based on the test set (Figure 3F). Likewise, the independency of the prognostic model was also confirmed by univariate and multivariate Cox regression analyses (Figure 5C, D).

Validation of the risk score in the entire set

We further validated our risk score in the entire set (516 samples). Based on the cut-offs in the train set, the samples were stratified into low-risk group (n=263) and high-risk group (n=253). Figure 4A-C also showed the risk score and survival status. KM plot indicated that patients in high-risk group and low-risk group showed significantly different outcomes (P= 3.335e-05) (Figure 4D). The 1‐year, 3‐year and 5‐year prognostic accuracy of the prognostic model was 0.674, 0.712 and 0.659, respectively (Figure 4E). Principal component analysis displayed a different distribution pattern of high and low risk according to 5 autophagy-related gene expression based on the entire set (Figure 4F). The risk score and clinical factors were incorporated into univariate and multivariate Cox regression analyses. The results demonstrated that the ARG risk score in the model can serve as an independent prognostic indicator (Figure 5E, F). These conclusions supported our speculation that autophagy affects the prognosis of EC.

Relationships between the risk score and clinicopathologic factors

We finally built a complete prognostic model based on the entire set. This model showed a good performance in stratification in clinical stage I-II & III-IV, grade G3&G4, peritoneal wash, cancer status, radiation therapy, minimally invasive surgical approach and endometrioid histological type (Figure 6). In parallel to the results above, high‐risk group in both subgroups was inclined to have worse OS. Meanwhile, relationships were analyzed between the ARGs-based risk score and clinicopathologic factors, including age, grade, clinical stage, peritoneal wash, radiation therapy, surgical approach, cancer status and histological type. Risk score was significantly higher in seniors, poorly differentiated cases, and mix or serous histological type cases (Additional file 3). We also explored the clinical significance of included genes (Table 3).

Nomogram and its clinical utility

To provide clinicians with a quantitative method to predict the overall survival of EC patients, a nomogram was constructed incorporating the risk score and clinical factors (Figure 7A). The calibration curve of 1-year, 3-year and 5-year survival indicated that the nomogram was almost an ideal model (Figure 7B). ROC curve analysis also showed that in this model, the risk score AUC was 0.794, and the clinical factors AUC was 0.851, both much significantly higher than that of patient age (AUC=0.621), clinical stage (AUC=0.701), grade (AUC=0.691), peritoneal wash (AUC=0.726), cancer status (AUC=0.760) and histological type (AUC=0.746). Interestingly, the ROC curve of risk score combined with clinical factors was much higher than that of each alone (AUC=0.859) (Figure 7C). These further validated the reliability of our previous prognostic model.

External validation of the risk score

The representative protein expression of CDKN1B, DLC1, EIF4EBP1 and ERBB2 was explored in the Human Protein Profiles (Figure 8A). However, we failed to find GRID1 on the website. ERBB2 possessed the most genetic alterations (11%), amplification being the commonest (Figure 8B). Moreover, ROC curve analysis found that the key genes had strong capacity to distinguish tumor from normal tissues, and the combined use of these ARGs was more effective (AUC=0.997) (Figure 8C).

Gene set enrichment analysis

We performed GSEA to identify ARGs related pathways. The most enriched pathways in high-risk groups included “Cell cycle”, “Citrate cycle tca cycle”, “DNA replication”, “Homologous recombination”, “Mismatch repair” and “Oxidative phosphorylation”, which implied that the five ARGs are related to EC progress (Additional file 4).

Somatic mutation analysis

The top ten mutated genes in different risk groups were displayed in the Figure 12. The order of somatic mutations in the high-risk group was as follows: TP53 >PTEN >PIK3CA >TTN >ARID1A >KMT2D >PIK3R1 >MUC16 >ZFHX3 > FAT4 (Figure 9A). In the low-risk group, the order was: PTEN >ARID1A > PIK3CA >TTN >PIK3R1 >CTNNB1 >CTCF >KMT2D >MUC16 >CSMD3 (Figure 9B). Besides, we found that the patients in the high-risk group showed significant mutation signatures, compared with those in the low-risk group (Figure 9C, D). It was noticed that TP53, PTEN, CTNNB1, DOCK3, ARID1A, ARHGAP35, FAT4, CTCF and FBXW7 were mutated differentially between the high- and low-risk groups (P<0.05) (Figure 9E).

Expression levels of ARGs in clinical tissue samples

Finally, the qRT-PCR results from our own specimens verified that the expression of EIF4EBP1 and ERBB2 mRNA were significantly higher in EC tissues than that in normal endometrial tissues, while the expression of CDKN1B, DLC1 and GRID1 mRNA were lower. The p values were 0.002, 0.037, 0.021, 0.012 and 0.003 respectively (Figure 10).

Endometrial cancer is a major gynecological malignancy worldwide. Prognosis-related high risk factors include disease stage, tumor size, grade, histological type, myometrial invasion and lymph node metastasis [14]. Nevertheless, these risk factors could only be determined after extensive surgery. Consequently, preoperative prognostic factors are needed, especially for those who are not suitable for surgery or wish to preserve their fertility function. Similar to other tumors, the occurrence, development and metastasis of EC also involve complex molecular mechanism [15]. Thus, reliable molecular biomarkers predicting the prognosis of EC are significant for selecting patients who might be sensitive to additional targeted therapy [16]. Autophagy is a highly conserved catabolic and energy-generating process that maintains cell homeostasis through degrading damaged organelles or intracellular components [7]. It was reported that autophagy plays complicated and momentous roles in EC genesis and progression. In type I endometrial cancer, as autophagy is repressed by activated PI3K/AKT/mTOR signaling pathway, the proliferation of tumor cells becomes facilitated, indicating that autophagy inhibits the development of tumors. Hua Liu, et al. showed that activation of the PI3K/AKT pathway could result in activation of mTOR and the proliferation of endometrial cells by inhibiting autophagy [18]. Currently, autophagy is considered as a "double-edged sword", since it not only inhibits tumor progression by depriving cells of nutrients and energy, but also promotes tumor progression by conferring cells with nutrients and energy[19, 20]. Therefore, autophagy-related research may open a new window into EC pathogenesis and prognosis. Several autophagy-related prognostic indicators had already been proposed. For example, An Y, et al. developed an eight-gene autophagy-related signature, which could serve as an independent prognostic indicator for serous ovarian cancer [21]. Wang SS, et al. identified three prognostic ARGs (JUN, MYC, and ITGA3) and validated a clinical autophagy-based index in predicting overall survival of bladder cancer [22]. However, the prognostic significance of autophagy genes in endometrial cancer has not been reported.

In the present study, we first screened differentially expressed ARGs in EC. Subsequently, several gene functional enrichment analyses were performed to define their biological process, molecular function and signaling pathways. Via univariate and multivariate survival analysis, we identified that CDKN1B, DLC1, EIF4EBP1, ERBB2 and GRID1 were significantly related to EC patients’ overall survival. A prognostic model was constructed based on them, all of which showed potential of being therapeutic targets. CDKN1B encodes p27Kip1 protein, a cyclin-dependent kinase inhibitor that binds to and prevents the activation of a broad range of cyclin CDK complexes, thereby controlling the cell cycle progression at G1 [23]. Recent studies have demonstrated that CDKN1B can be genetically inactivated, particularly in luminal breast cancer, prostate cancer, and small intestine neuroendocrine neoplasms [24]. DLC1 gene encodes GTPase-activating protein (GAP), a member of the Rho-GAP family of proteins, and plays a role in the regulation of small GTP-binding proteins [25]. GAP proteins are involved in signaling pathways that regulate cellular processes in cytoskeletal changes, functioning as a tumor suppressor in a variety of common cancers, such as lung cancer, liver cancer, colorectal cancer, breast cancer and prostate cancer [26]. This role is consistent with our findings. EIF4EBP1 encodes one translation repressor protein, which directly interacts with eukaryotic translation initiation factor 4E (eIF4E), to inhibit complex assembly and repress translation. This protein is phosphorylated in response to various signals, resulting in its dissociation from eIF4E and activation of mRNA translation [27, 28]. High level of EIF4EBP1 phosphorylation is associated with malignant progression and adverse prognosis in several malignancies, such as breast, ovary and prostate cancers [29]. ERBB2 encodes an epidermal growth factor (EGF) receptor and binds tightly to other ligand-bound EGF receptors to form a heterodimer, stabilizing ligand binding and enhancing kinase-mediated activation of downstream signaling pathways, such as phosphatidylinositol-3 kinase and mitogen-activated protein kinase. Amplification and/or overexpression of this gene have been reported in lots of cancers, including ovarian and breast cancers [30]. ERBB2 overexpression shows association with EC prognosis [31, 32]. GRID1 is a member of glutamate receptors, which were thought to exclusively play a role in the central nervous system (CNS). However, it was recently suggested that glutamate receptor subunits are also involved outside the CNS in malignant processes, such as pancreatic cancer, breast cancer and ovarian cancer [33, 34]. From a clinical perspective elevated expression in high-risk endometrial cancers makes the gene, which is reported in the present study for the first time attractive as a potential tumor marker.

Our gene enrichment analysis also confirmed their involvement in cancer-related biological processes. Besides, higher rates of TP53, ZFHX3 and FAT4 mutation were found in high-risk patients than low-risk patients. TP53 is a critical tumor suppressor that regulates cell cycle progression, apoptosis, cell senescence and many others. Its mutation evokes diverse physiological effects, including numerous cancers [35]. Therefore, TP53 mutation can be used to guide the design of new EC treatments. Our GSEA analysis suggested that high risk score tended to be accompanied by a number of up-regulated networks, including tumorigenesis and tumor progression associated pathways. For example, nstability is a hallmark of cancer, most likely caused by abnormal DNA replication [36]. Homologous recombination enables the cell to obtain and replicate intact DNA sequence information，especially to repair DNA damage. Patients with homologous recombination deficiencies were more responsive to treatments, such as platinum-based chemotherapy and poly (ADP-ribose) polymerase (PARP) inhibitors [37, 38].

Using TCGA, we screened out autophagy-related genes and used them to build up a prognostic risk score model. The model displayed powerful ability to predict the survival in patients stratified according to ages, grades, clinical stages and histological types. Besides, the model was validated using the test set and the entire set. Thus, it might be converted into clinical application.

However, there are some limitations in this study. First, the underlying molecular mechanism of the key ARGs in the EC pathogenesis were undefined. Second, other potential prognostic variables associated with EC, such as tumor size, myometrial invasion and lymph node metastasis, should be investigated. Third, prospective studies are required to further confirm the clinical utility and the biological function of this signature.

In conclusion, based on five autophagy-related genes (CDKN1B, DLC1, EIF4EBP1, ERBB2 and GRID1), our model can independently predict the OS of EC patients by combining molecular signature and clinical characteristics.

Ethics approval and consent to participate

This research was approved by the Institutional Review Board of Nanjing Medical University and informed consent was obtained from all subjects.

Consent for publication

Not applicable

Availability of data and materials

The data sets used and/or analyzed during the current study are publicly available data from The Cancer Genome Atlas (TCGA). The figures and materials supporting the conclusions of this article are included within the article and its additional files).

Competing interests

The authors declare that there is no conflict of interests.

Funding

This work was supported by the National Nature Science Foundation of China (81872119 and 81472442) and the Jiangsu province medical innovation team (CXTDA2017008).

Authors' contributions

Yankai Xia and Wenjun Cheng conceived and designed the model. Hui Wang, Xiaoling Ma and Jinhui Liu analyzed the data and wrote the manuscript. Yicong Wan and Yi Jiang contributed significantly to the data analyses and provided important suggestions. All authors provided critical feedback and helped shape the research, analyses and manuscript.

Acknowledgements

We would like to acknowledge the Human Autophagy Database and TCGA for free use. The authors received no financial support for the research and authorship of this article.

ARGs: Autophagy-related genes

AUC: Area under the curve

CNS: Central nervous system

EC: Endometrial cancer

GO: Gene ontology

GSEA: Gene Set Enrichment Analysis

KEGG: Kyoto Encyclopedia of Genes and Genomes

MAF: Mutation annotation format

MSigDB: Molecular signature database

OS: Overall survival

qRT-PCR: Quantitative real-time PCR

RNA-seq: RNA sequencing

ROC: Receiver operating characteristic

SLS: Stone-like structure

TCGA: The Cancer Genome Atlas

[1] J.M.M. Levy, C.G. Towers, A. Thorburn. Targeting autophagy in cancer. Nat Rev Cancer. 2017; 17(9):528-542.

[2] K.R. Parzych, D.J. Klionsky. An overview of autophagy: morphology, mechanism, and regulation. Antioxid Redox Signal. 2014; 20(3): 460-473.

[3] B. Levine, G. Kroemer. Autophagy in the pathogenesis of disease. Cell. 2008; 132(1):27-42.

[4] Siegel RL, Miller KD, Jemal A. Cancer stastics, 2019. CA Cancer J Clin. 2019; 69(1):7-34.

[5] Stefano Uccella, Matteo Bonzini, Stefano Palomba, Francesco Fanfani, Marcello Ceccaroni, Renato Seracchioli, Enrico Vizza , Annamaria Ferrero, Giovanni Roviglione, Paolo Casadio, Giacomo Corrado, Giovanni Scambia, Fabio Ghezzi. Impact of obesity on surgical treatment for endometrial cancer: a multicenter study comparing laparoscopy vs open surgery, with propensity-matched analysis. J Minim Invasive Gynecol. 2016; 23 (1):53-61.

[6] Elizabeth V Connor, Peter G Rose. Management Strategies for Recurrent Endometrial Cancer. Expert Rev Anticancer Ther. 2018; 18(9):873-885.

[7] Balch C, Matei DE, Huang TH, Nephew KP. Role of epigenomics in ovarian and endometrial cancers. Epigenomics. 2010; 2 (3):419-447.

[8] E. Sivridis, A. Giatromanolaki, V. Liberis, M.I. Koukourakis. Autophagy in endometrial carcinomas and prognostic relevance of 'stone-like' structures (SLS): what is destined for the atypical endometrial hyperplasia?. Autophagy. 2011; 7(1):74-82.

[9] N. Eritja, B.J. Chen, R. Rodriguez-Barrueco, M. Santacana, S. Gatius, A. Vidal, M.D. Marti, J. Ponce, L. Bergada, A. Yeramian, M. Encinas, J. Ribera, J. Reventos, J. Boyd, A. Villanueva, X. Matias-Guiu, X. Dolcet, D. Llobet-Navas. Autophagy orchestrates adaptive responses to targeted therapy in endometrial cancer. Autophagy. 2017; 13(3):608-624.

[10] Maxim V. Kuleshov , Matthew R. Jones, Andrew D. Rouillard, Nicolas F. Fernandez, Qiaonan Duan, Zichen Wang, Simon Koplev, Sherry L. Jenkins, Kathleen M. Jagodnik, Alexander Lachmann, Michael G. McDermott, Caroline D. Monteiro, Gregory W. Gundersen, Avi Ma’ayan. Enrichr: a comprehensive gene set enrichment analysis web server 2016 update. Nucleic Acids Res. 2016; 90-97.

[11] Aravind Subramanian, Pablo Tamayo, Vamsi K. Mootha, Sayan Mukherjee, Benjamin L. Ebert, Michael A. Gillette, Amanda Paulovich, Scott L. Pomeroy, Todd R. Golub, Eric S. Lander, and Jill P. Mesirov. Gene Set Enrichment Analysis: A Knowledge-Based Approach for Interpreting Genome-Wide Expression Profiles. Proc Natl Acad Sci U S A. 2015;102 (43):15545-15550.

[12] Daniele Merico, Ruth Isserlin, Oliver Stueker, Andrew Emili, Gary D. Bader. Enrichment Map: A Network-Based Method for Gene-Set Enrichment Visualization and Interpretation. PLoS One. 2010; 5 (11): 13984.

[13] Anand Mayakonda, De-Chen Lin, Yassen Assenov, Christoph Plass, H. Phillip Koeffler. Maftools: efficient and comprehensive analysis of somatic variants in cancer. Genome Res. 2018; 28 (11) :1747-1756.

[14] T. Van Nyen, C.P. Moiola, E. Colas, D. Annibali, F. Amant. Modeling Endometrial Cancer: Past, Present, and Future. Int J Mol Sci. 2018; 19(8).

[15] Albiona Stampoliou, Petroula Arapantoni-Dadioti, Kitty Pavlakis. Epigenetic mechanisms in endometrial cancer. J BUON. 2016; 21 (2):301-306.

[16] Nadia Taoussi, Ali Alghamdi, Konrad Futyma, Tomasz Rechberger. Biological markers with potential clinical value in endometrial cancer — review of the literature. Ginekol Pol. 2017; 88 (6):331-336.

[17] Manuela Antonioli, Martina Di Rienzo, Mauro Piacentini, Gian Maria Fimia. Emerging Mechanisms in Initiating and Terminating Autophagy. Trends Biochem Sci. 2017;42 (1):28-41.

[18] H. Liu, L. Zhang, X. Zhang, Z. Cui. PI3K/AKT/mTOR pathway promotes progestin resistance in endometrial cancer cells by inhibition of autophagy. Onco Targets Ther. 2017; 10:2865-2871.

[19] S.K. Bhutia, S. Mukhopadhyay, N. Sinha, D.N. Das, P.K. Panda, S.K. Patra, T.K. Maiti, M. Mandal, P. Dent, X.Y. Wang, S.K. Das, D. Sarkar, P.B. Fisher. Autophagy: cancer's friend or foe?. Adv Cancer Res. 2013; 118:61-95.

[20] S. Lorin, A. Hamai, M. Mehrpour, P. Codogno. Autophagy regulation and its role in cancer. Semin Cancer Biol. 2013; 23(5):361-379.

[21] Yuanyuan An, Fangfang Bi, Yue You, Xinhui Liu, Qing Yang. Development of a Novel Autophagy-related Prognostic Signature for Serous Ovarian Cancer. J Cancer. 2018; 9 (21):4058-4071.

[22] S.S. Wang, G. Chen, S.H. Li, J.S. Pang, K.T. Cai, H.B. Yan, Z.G. Huang, R.Q. He. Identification and validation of an individualized autophagy-clinical prognostic index in bladder cancer patients. Onco Targets Ther. 2019; 12:3695-3712.

[23] Kornelia Polyak, Jun-ya Kato, z Mark J. Solomon, Charles J. Sherr, z Joan Massague, James M. Roberts, Andrew Koff. p27 Kip1, a cyclin-Cdk inhibitor, links transforming growth factor-13 and contact inhibition to cell cycle arrest. Genes Dev. 1994; 8 (1):9-22.

[24] M. Cusan, G. Mungo, M. De Marco Zompit, I. Segatto, B. Belletti, G. Baldassarre. Landscape of CDKN1B Mutations in Luminal Breast Cancer and Other Hormone-Driven Human Tumors. Front Endocrinol (Lausanne). 2018; 9:393.

[25] B Z Yuan 1, M J Miller, C L Keck, D B Zimonjic, S S Thorgeirsson, N C Popescu. Cloning, Characterization, and Chromosomal Localization of a Gene Frequently Deleted in Human Liver Cancer (DLC-1) Homologous to Rat RhoGAP. Cancer Res. 1998; 58 (10):2196-2199.

[26] Dunrui Wang, Xiaolan Qian, Megha Rajaram, Marian E. Durkin, Douglas R.

Lowy. DLC1 is the principal biologically-relevant down-regulated DLC family member in several cancers. Oncotarget. 2016; 7 (29):45144-45157.

[27] Ivan Topisirovic, Melisa Ruiz-Gutierrez, Katherine L. B. Borden. Phosphorylation of the Eukaryotic Translation Initiation Factor eIF4E Contributes to Its Transformation and mRNA Transport Activities. Cancer Res. 2004; 64 (23): 8639-8642.

[28] E. Petroulakis, A. Parsyan, R.J. Dowling, O. LeBacquer, Y. Martineau, M. Bidinosti, O. Larsson, T. Alain, L. Rong, Y. Mamane, M. Paquet, L. Furic, I. Topisirovic, D. Shahbazian, M. Livingstone, M. Costa-Mattioli, J.G. Teodoro, N. Sonenberg. p53-dependent translational control of senescence and transformation via 4E-BPs. Cancer Cell. 2009; 16(5):439-446.

[29] X. Qin, B. Jiang, Y. Zhang. 4E-BP1, a multifactor regulated multifunctional protein. Cell Cycle. 2016; 15(6):781-786.

[30] M.M. Moasser. The oncogene HER2: its signaling and transforming functions and its role in human cancer pathogenesis. Oncogene. 2007; 26(45):6469-6487.

[31] W. Xiao, X. Dong, H. Zhao, S. Han, R. Nie, X. Zhang, R. An. Expression of MIF and c-erbB-2 in endometrial cancer. Mol Med Rep. 2016; 13(5)3828-3834.

[32] M.K. Halle, I.L. Tangen, H.F. Berg, E.A. Hoivik, K.K. Mauland, K. Kusonmano, A. Berg, A. Hurtado, K.H. Kalland, A.M. Oyan, I. Stefansson, O.K. Vintermyr, H.M. Werner, I.S. Haldorsen, J. Trovik, H.B. Salvesen, C. Krakstad. HER2 expression patterns in paired primary and metastatic endometrial cancer lesions. Br J Cancer. 2018; 118(3):378-387.

[33] Kristian Unger, Johannes Wienberg, Andrew Riches, Ludwig Hieber, Axel Walch, Andreas Brown, Patricia C M O’Brien, Ca¨cilia Briscoe, Lindsey Gray, Elke Rodriguez, Gerhard Jack, Jeroen Knijnenburg, Giovanni Tallini, Malcolm Ferguson-Smith, Horst Zitzelsberger. Novel Gene Rearrangements in Transformed Breast Cells Identified by High-Resolution Breakpoint Analysis of Chromosomal Aberrations. Endocr Relat Cancer. 2010; 17 (1): 87-98.

[34] Andrzej Stepulak, Radoslaw Rola, Krzysztof Polberg, Chrysanthy Ikonomidou. Glutamate and Its Receptors in Cancer. J Neural Transm (Vienna). 2014; 121 (8):933-944.

[35] James A. McCubrey, Kvin Lertpiriyapong, Timothy L. Fitzgerald, Alberto M. Martelli, Lucio Cocco, Dariusz Rakus, Agnieszka Gizak, Massimo Libra, Melchiorre Cervello, Guiseppe Montalto, Li V. Yang , Stephen L. Abrams, Linda S. Steelman. Roles of TP53 in determining therapeutic sensitivity, growth, cellular senescence, invasion and metastasis. Adv Biol Regul. 2017; 63:32-48.

[36] Hélène Gaillard, Tatiana García-Muse, Andrés Aguilera. Replication Stress and Cancer. Nat Rev Cancer. 2015; 15 (5):276-289.

[37] William Douglass Wright, Shanaya Shital Shah, Wolf-Dietrich Heyer. Homologous recombination and the repair of DNA double-strand breaks. J Biol Chem. 2018; 293 (27):10524-10535.

[38] Renata Rodrigues da Cunha Colombo Bonadio, Rodrigo Nogueira Fogace, Vanessa Costa Miranda, Maria del Pilar Estevez Diz. Homologous recombination deficiency in ovarian cancer:a review of its epidemiology and management. Clinics. 2018;73(suppl 1):e450s.

Table 1 Clinicopathologic characteristics of 342 EC patients from TCGA database.

Characteristics	Entire set (n=342)
N (%)
Age, years
≤60	141 (41.2%)
＞60	201 (58.8%)
Stage
I	221 (64.6%)
II	31 (9.1%)
III	77 (22.5%)
IV	13 (3.8%)
Grade
G1	80 (23.4%)
G2	74 (21.6%)
G3	188 (55.0%)
Peritoneal wash
negative	293 (85.7%)
positive	49 (14.3%)
Person neoplasm cancer status
tumor free	291 (85.1%)
with tumor	51 (14.9%)
Radiation therapy
yes	140 (40.9%)
no	202 (59.1%)
Histological type
endometrioid endometrial adenocarcinoma	268 (78.4%)
serous endometrial adenocarcinoma	61 (17.8%)
mixed serous and endometrioid	13 (3.8%)
Surgical approach
minimally Invasive	150 (43.9%)
open	192 (56.1%)

Table 2. The most prognosis-related five autophagy-related genes.

gene	coef	HR	HR.95L	HR.95H	pvalue
CDKN1B	0.016941	1.017085	0.994197	1.040501	0.144618
DLC1	-0.33079	0.718355	0.542355	0.951468	0.02106
EIF4EBP1	0.003577	1.003584	1.000417	1.00676	0.026523
ERBB2	0.001271	1.001272	1.000166	1.002379	0.024188
GRID1	0.733501	2.082357	1.20684	3.59303	0.008402

Download PDF

Journal Publication

published 13 Jul, 2020

Read the published version in Cancer Cell International →

Editorial decision: Accept
08 Jul, 2020
Review #2 received at journal
06 Jul, 2020
Review #1 received at journal
03 Jul, 2020
Reviewer #3 agreed at journal
01 Jul, 2020
Reviewers invited by journal
30 Jun, 2020
Reviewer #1 agreed at journal
30 Jun, 2020
Reviewer #2 agreed at journal
30 Jun, 2020
Editor assigned by journal
17 Jun, 2020
Submission checks completed at journal
16 Jun, 2020
Editor invited by journal
16 Jun, 2020

You are reading this latest preprint version

Prognostic value of an autophagy-related gene expression signature for endometrial cancer patients

Status:

Journal Publication

Version 3

Abstract

Figures

Background

Methods

Results

Discussion

Conclusions

Declarations

Abbreviations

References

Tables

Supplementary Files

Status:

Journal Publication

Version 3