- Department of Respiratory and Critical Care Medicine, Union Hospital, Tongji Medical College, Huazhong University of Science and Technology, Wuhan, China
The role of autophagy in lung cancer is context-dependent and complex. Recent studies have reported the important role of autophagy in tumor immune escape. However, the association between autophagy and tumor-infiltrating lymphocytes (TILs) in early-stage lung adenocarcinoma (LUAD) remains unclear. In this study, we aimed to develop and validate the autophagy-related gene pair index (ATGPI) and autophagy clinical prognostic index (ACPI) in multiple LUAD cohorts, including The Cancer Genome Atlas (TCGA) cohort, Gene Expression Omnibus cohorts, and one cohort from Union Hospital, Wuhan (UH cohort), using a Cox proportional hazards regression model with the least absolute shrinkage and selection operator. Multivariate Cox regression analysis demonstrated that there was a significant difference in overall survival (OS) between patients with high and low ATGPI in the testing [hazard ratio (HR) = 1.97; P < 0.001] and TCGA validation (HR = 2.25; P < 0.001) cohorts. Time-dependent receiver operating characteristic curve analysis was also performed. We found that high ATGPI could accurately identify patients with early-stage LUAD with shorter OS, with the areas under the curve of 0.703 and 0.676 in the testing and TCGA validation cohorts, respectively. Concordance index (C-index) was used to evaluate the efficiency of ATGPI and ACPI. The C-index of ACPI was higher than that of ATGPI in the testing (0.71 vs. 0.66; P < 0.001), TCGA validation (0.69 vs. 0.65; P = 0.028), and UH (0.80 vs. 0.70; P = 0.015) cohorts. TIL analysis demonstrated that the proportions of tumor-infiltrating CD4+ T cells were lower in the high-ATGPI group than in the low-ATGPI group in both the TCGA validation and UH cohorts. These results indicate the potential clinical use of ATG signatures which are associated with TILs, in identifying patients with early-stage LUAD with different OS.
Introduction
Macroautophagy (referred to here as autophagy) is a highly conserved catabolic process that involves the formation of autophagosomes, growth of double membranes, merging with lysosomes, and disintegration of engulfed cellular proteins and organelles (Levine and Kroemer, 2008; Levy et al., 2017). Autophagy plays a context-dependent role in cancer development. The enhancement and inhibition of autophagy have both been suggested as therapeutic strategies in cancers (Levy and Thorburn, 2011; Amaravadi et al., 2016). The dynamic process of autophagy in cancers is coordinated by a series of proteins encoded by autophagy-related (ATG) genes and their various activation factors such as the tumor type and status of oncogenes and tumor suppressors (Kimmelman and White, 2017). The role of autophagy in lung cancer has been previously demonstrated in multiple genetically engineered mouse models. Inactivation of the essential autophagy gene ATG5 accelerates the early phases of non-small cell lung cancer (NSCLC) oncogenesis (Rao et al., 2014). Additionally, autophagy maintains the functioning of mitochondria to support the metabolism and growth of KRAS-driven NSCLC (Guo and White, 2013). Moreover, multiple factors, including circular RNA, REV-ERB agonists, and copper, can function as autophagy regulators to modulate the proliferation, migration, and invasion of lung cancer (Chen et al., 2020; Shen et al., 2020; Tsang et al., 2020). These findings demonstrate the biological significance of autophagy in lung cancer. However, few studies have focused on the potential of ATG genes in predicting survival of patients with early-stage lung adenocarcinoma (LUAD).
Lung cancer is the most frequently diagnosed cancer and the most common cause of cancer-related mortality worldwide, and LUAD is the most common histologic subtype of NSCLC (Torre et al., 2015). LUAD is currently considered a cluster of discrete molecular subtypes that are defined by high-throughput genomic profiling. Therefore, it is crucial to assess the complex genomic alterations of LUAD using next-generation sequencing technologies The Cancer Genome Atlas Research Network (2014). Combining sequencing profiles with clinical characteristics will improve the prognostic evaluation of NSCLC. To date, multiple prognostic biomarkers, including long non-coding RNA signatures, microRNA signatures, circulating tumor DNA signatures, DNA methylation signatures, and immune signatures, have been constructed for NSCLC (Ajona et al., 2013; Hao et al., 2017; Li B. et al., 2017; Li Y. et al., 2017; Ooki et al., 2017; Martínez-Terroba et al., 2018; Dong et al., 2019). However, the clinical application of these signatures is limited because of technical biases across different platforms.
Thus, in this study, we aimed to develop and validate the autophagy-related gene pair index (ATGPI) and autophagy clinical prognostic index (ACPI) in multiple LUAD cohorts. ATGPI and ACPI could accurately distinguish early-stage LUAD with different overall survival (OS) in independent cohorts. Compared with ATGPI, the combination with clinical characteristics improved the accuracy of ACPI. The proportions of tumor-infiltrating CD4+ T cells were lower in the high-ATGPI group than in the low-ATGPI group.
Materials and Methods
Data Collection and LUAD Cohorts
Gene expression data from 13 public LUAD datasets and one LUAD cohort from Union Hospital, Wuhan, China were used in this retrospective study (Bhattacharjee et al., 2001; Beer et al., 2002; Takeuchi et al., 2006; Kuner et al., 2009; Tomida et al., 2009; Okayama et al., 2012; Wilkerson et al., 2012; Rousseaux et al., 2013; Sato et al., 2013; Tang et al., 2013; Der et al., 2014). In total, 68 samples of LUAD tissues were obtained from surgical resections of patients without preoperative treatment from Union Hospital, Wuhan, China. Specimens were collected from July 2017 to August 2020 and were stored at −80°C. Human specimens were obtained with the consent of patients; the study was approved by Ethics Committees of Tongji Medical College, Huazhong University of Science and Technology (2019-S671). RNA-Seq data and clinical information of LUAD patients were downloaded from Genomic Data Commons Data Portal1. Data from 12 datasets (GSE11969, GSE13213, GSE14814, GSE26939, GSE30219, GSE31210, GSE31547, GSE37745, GSE41271, GSE42127, GSE50081, and GSE68465) were downloaded from the Gene Expression Omnibus (GEO)2. Patients who received neoadjuvant therapy, adjuvant chemotherapy, or other immune-modulating therapy were excluded. All microarray datasets were merged into one dataset and then randomly divided into training and testing cohorts. The Cancer Genome Atlas (TCGA) and Union Hospital (UH) cohorts were included as independent validation cohorts. The information about all platforms included in this study is shown in Table 1. Raw median intensity values of tumor samples in Agilent datasets were extracted using the R package agilp (v3.4.0) (Stafford and Brun, 2007). For Affymetrix microarrays, gene expression profiles were normalized using the R package affy (version 1.50.0) (Gautier et al., 2004). For Illumina datasets, model-based background correction processed data were downloaded.
Construction of ATGPI as a Robust Prognostic Biomarker for Early-Stage LUAD
The method used to construct the gene pair index was based on that used in a previous study (Li B. et al., 2017). The ATG genes were downloaded from The Human Autophagy Database3. Among them, only 142 ATG genes detected in all the datasets were included. To construct the autophagy-related gene pairs (ATGPs), the ATG genes were compared in pairs. To construct one gene pair, the ATGP score was 0 if the gene expression level of the ATG gene1 was more than that of the ATG gene2; otherwise, the ATGP score was 1. The ATGPs with constant values (0 or 1) in each dataset were excluded considering the probable biases. The log-rank test was used to select ATGPs that were associated with the OS of patients. ATGPs with p values < 0.05 were included in ATGPI in the training cohort. A Cox proportional hazards regression model with the least absolute shrinkage and selection operator (LASSO) was used to select appropriate ATGPs. The optimal model parameter λ was evaluated by 10-fold cross-validation at 1 standard error as previously recommended. The final ATGPI was constructed using 10 ATGPs that contained 15 unique ATG genes. Complete information of the final 10 ATGPs is provided in Table 2. Time-dependent receiver operating characteristic (ROC) curves were used to determine the optimal cutoff of ATGPI. The nearest neighbor estimation method was used to estimate the ROC curve. Patients were divided into low- and high-risk groups based on the optimal cutoff value of ATGPI.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719011/fcell-09-719011-HTML/image_m/fcell-09-719011-t002.jpg)
Table 2. The autophagy-related (ATG) genes used to construct the autophagy-related gene pair index (ATGPI).
Validation of the Prognostic Value of ATGPI
To assess the prognostic value of ATGPI, univariate Cox regression analysis was performed in the training, testing, and TCGA validation cohorts. Age and tumor stage were considered continuous variables. The tumor stage was considered as I = 1 and II = 2. Male and high ATGPI were risk factors in the univariate analysis. The prognostic accuracy of ATGPI was assessed using time-dependent ROC curves.
Association Between ATGPI and Immune Infiltration
To explore the association between ATGPI and immune infiltration, we used the ImmuCellAI method4 (Miao et al., 2020) to analyze tumor-infiltrating lymphocytes (TILs) abundance in the TCGA validation cohort. The proportions of TILs were compared between the high- and low-ATGPI groups.
Combination of ATGPI and Clinical Characteristics
Multivariate Cox regression analysis was conducted by combining ATGPI with clinical factors that were significantly associated with OS in the univariate analysis. ATGPI, age, and tumor stage were combined to fit a Cox proportional hazards regression in the training cohort. The cutoff value of ACPI was estimated using time-dependent ROC curves. The prognostic accuracy of ACPI was compared with that of ATGPI using restricted mean survival (RMS) curves.
Total RNA Isolation and Quantitative Real-Time Polymerase Chain Reaction (qRT-PCR)
Total RNAs were isolated from fresh LUAD tissues using Trizol reagent (Takara Bio Inc., China). HiScript II One Step qRT-PCR SYBR Green Kit (Vazyme, China) was used in qRT-PCR analysis. The sequences of all primers used are shown in Table 3. All primers were designed and purchased from Tsingke (Wuhan, China).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719011/fcell-09-719011-HTML/image_m/fcell-09-719011-t003.jpg)
Table 3. The primer sequences of 15 autophagy-related (ATG) genes used in quantitative real-time polymerase chain reaction (qRT-PCR) analysis.
Immunohistochemical (IHC) Analyses
The tumor tissues collected from LUAD patients were fixed, dehydrated, and embedded in paraffin. The 4-μm-thick sections were then stained for IHC analyses using recombinant anti-CD4 antibodies (ab133616, Abcam, United Kingdom).
Statistical Analysis
The log-rank test was used to select ATGPs associated with OS. The LASSO Cox proportional hazards regression model was used to select ideal prognostic ATGPs to construct ATGPI. The time-dependent ROC curve was used to determine the optimal cutoff of ATGPI and ACPI. The discrimination of the model was evaluated using the areas under the curve (AUCs). The Kaplan–Meier method was used to estimate survival curves between the low- and high-risk groups. A multivariate Cox proportional hazards regression model that included variates significantly associated with OS in the univariate analysis was used. The concordance index (C-index) was calculated using survcomp (version 1.22.0) and was compared using compareC (version 1.3.1) packages (Kang et al., 2015). The RMS curve and RMS time ratio were estimated using survival (version 2.44-1.1) and survRM2 (version 1.0-2) packages. All statistical analyses were performed using R (version 3.6.1), P values < 0.05 were considered statistically significant.
Results
Clinical Characteristics of Patients
Detailed patient characteristics are shown in Table 4. A total of 1888 LUAD patients from multiple datasets were included. No statistically significant difference was observed between training and testing cohorts with respect to clinical characteristics. The analysis process of this study was demonstrated in a flow chart (Figure 1).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719011/fcell-09-719011-HTML/image_m/fcell-09-719011-g001.jpg)
Figure 1. Study design for the construction and validation of the autophagy-related gene pair index (ATGPI) and autophagy clinical prognostic index (ACPI) in early-stage lung adenocarcinoma.
Construction of ATGPI
Among the ATG genes from The Human Autophagy Database, 142 ATG genes were detected in all public datasets enrolled in this study. A total of 10011 ATGPs were built. Among 5402 ATGPs with differential expression in all datasets, 96 ATGPs that were associated with OS in the training dataset were selected to fit the LASSO regression model. Finally, 10 ATGPs that contained 15 unique ATG genes were selected to construct the ATGPI (Figures 2A,B). We then built the ATGPI using Cox proportional hazards regression. Based on the time-dependent ROC curve analysis, the optimal cutoff value for ATGPI was 1.35 (Figure 2C). Risk curve and scatterplot were used to determine ATGPI and the vital status of each patient. Patients in the high-ATGPI group had higher mortality than those in the low-ATGPI group (Figures 2D,E). Survival curves of different groups were estimated using the Kaplan–Meier method and were compared using the log-rank test (P < 0.001) (Figure 2F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719011/fcell-09-719011-HTML/image_m/fcell-09-719011-g002.jpg)
Figure 2. Construction of the autophagy-related gene pair index (ATGPI), a classifier comprising 15 unique autophagy-related genes. (A) The penalty parameter was estimated by 10-fold cross-validation in the training data set at 1 SE beyond the minimum partial likelihood deviance. (B) The LASSO coefficient profiles of the 10 prominent autophagy-related gene pairs. (C) The optimal cutoff of the ATGPI for stratifying patients based on time-dependent receiver operator characteristic curve analysis in the training cohort. (D) The risk curve of the ATGPI in the training cohort. (E) The scatterplot of vital status of each patient in the training cohort. (F) The Kaplan-Meier survival curves of low and high ATGPI groups.
Validation of ATGPI in the Testing and TCGA Cohorts
Autophagy-related gene pair index in the testing and TCGA validation cohorts was calculated using the same model in the training cohort. ATGPI and the vital status of each patient with LUAD were demonstrated by risk curves and scatterplots in the testing (Figures 3A,B) and validation (Figures 3E,F) cohorts. The Kaplan–Meier method and log-rank test demonstrated that ATGPI stratified early-stage LUAD into different prognostic groups in the testing (P < 0.001) (Figure 3C) and TCGA validation (P < 0.001) (Figure 3G) cohorts. In the testing cohort, the AUCs at 1, 3, and 5 years were 0.745, 0.703, and 0.687, respectively (Figure 3D). In the TCGA validation cohort, the AUCs at 1, 3, and 5 years were 0.668, 0.676, and 0.623, respectively (Figure 3H).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719011/fcell-09-719011-HTML/image_m/fcell-09-719011-g003.jpg)
Figure 3. The risk curves, Kaplan-Meier (KM) curves, and time-dependent receiver operator characteristic (ROC) curves of patients in the testing and The Cancer Genome Atlas (TCGA) validation cohorts based on the autophagy-related gene pair index (ATGPI). (A,B) Risk curves and scatterplots of vital status for patients in the testing cohort. (C) The KM curves of low and high ATGPI groups in the testing cohort. (D) The area under the curves (AUC) of time-dependent ROC curves in the testing cohort. (E,F) Risk curves and scatterplots of vital status for patients in the TCGA validation cohort. (G) The KM curves of low and high ATGPI groups in the TCGA validation cohort. (H) The AUC of time-dependent ROC curves in the TCGA validation cohort.
ATGPI as an Independent Prognostic Biomarker Related to Immune Infiltration for Early-Stage LUAD
In early-stage (stage I and II) LUAD, ATGPI stratified patients into different prognostic groups in the testing [hazard ratio (HR) = 2.07; 95% confidence interval (CI), 1.64–2.62; P < 0.001] and TCGA validation (HR = 2.15; 95% CI, 1.46–3.14; P < 0.001) cohorts. In stage II LUAD, ATGPI retained its high prognostic accuracy in the testing (HR = 2.32; 95% CI, 1.56–3.45; P < 0.001) and the TCGA validation (HR = 2.63; 95% CI, 1.35–5.10; P = 0.003) cohorts (Figure 4). Multivariate analysis was performed with variables that were significantly associated with OS in the univariate analysis, including age, ATGPI, and stage. Multivariate Cox proportional hazards regression analyses revealed that high ATGPI was an independent risk factor for poor prognosis of early-stage LUAD in the training (HR = 1.95; 95%CI, 1.52–2.50; P < 0.001), testing (HR = 1.97; 95%CI, 1.56–2.50; P < 0.001) and TCGA validation (HR = 2.25; 95%CI, 1.51–3.36; P < 0.001) cohorts (Table 5). ImmuCellAI analysis demonstrated that the percentages of CD4+ T cells infiltration were significantly different between the high- and low- ATGPI groups (Figure 5). The mean level of CD4+ TILs in the low-ATGPI group was 1.3-fold that in the high-ATGPI group (P < 0.0001).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719011/fcell-09-719011-HTML/image_m/fcell-09-719011-g004.jpg)
Figure 4. Hazard ratio (HR) between high- and low-risk groups based on the autophagy-related gene pair index (ATGPI) was evaluated using univariate cox proportional hazards regression. The length of the horizontal line corresponds to the confidence interval, and the size of the blue square is proportional to the number of patients. The vertical line indicates that HR is equal to 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719011/fcell-09-719011-HTML/image_m/fcell-09-719011-t005.jpg)
Table 5. Univariate and multivariate analyses of autophagy-related gene pair index (ATGPI) and clinical characteristics in the training, testing, and The Cancer Genome Atlas (TCGA) validation data sets.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719011/fcell-09-719011-HTML/image_m/fcell-09-719011-g005.jpg)
Figure 5. The immune infiltration status of high- and low-autophagy-related gene pair index (ATGPI) groups in the TCGA validation cohort. The percentages of different tumor-infiltrating lymphocytes (TILs) were compared by Wilcoxon rank sum test. ns, P ≥ 0.05; *P < 0.05; **P < 0.01; ***P < 0.001; ****P < 0.0001.
Integrated Prognostic Model by Combining ATGPI With Clinical Characteristics
ATGPI, age, and tumor stage were independent prognostic factors for OS of patients with early-stage LUAD; thus, these variables were combined to fit a Cox proportional hazards regression model in the training cohort, and ACPI was built (ACPI = 0.6657 × ATGPI + 0.4445 × stage + 0.0199 × age). The optimal cutoff value of ACPI was determined to be 1.31 based on the time-dependent ROC curve analysis in the training cohort. The RMS curves of continuous ATGPI and ACPI demonstrated that ACPI was significantly better in predicting OS than ATGPI in the training (C-index: 0.66, 95% CI, 0.63–0.70 vs. C-index: 0.68, 95% CI, 0.65–0.71; P = 0.032), testing (C-index: 0.66, 95% CI, 0.62–0.69 vs. C-index: 0.71, 95% CI, 0.68–0.74; P < 0.001), and TCGA validation (C-index: 0.65, 95% CI, 0.60–0.71 vs. C-index: 0.69, 95% CI, 0.63–0.74; P = 0.028) cohorts (Figure 6).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719011/fcell-09-719011-HTML/image_m/fcell-09-719011-g006.jpg)
Figure 6. The restricted mean survival (RMS) curves for autophagy-related gene pair index (ATGPI) and autophagy clinical prognostic index (ACPI) in the training, testing, and The Cancer Genome Atlas (TCGA) validation cohorts. (A–C) The RMS curves for the ATGPI and ACPI in the training cohort, the testing cohort, and the TCGA validation cohort. Each point demonstrates the RMS time of corresponding ATGPI and ACPI value. The slope of RMS curves was larger, the estimation for overall survival was better. Concordance index (C-index) was used to represent the difference between the ATGPI and ACPI.
Validation of ATGPI and ACPI in the UH Cohort
Detailed qRT-PCR data of patients are shown in Supplementary Table 1. ATGPI and ACPI of the UH cohort were calculated using the same model as that in the training cohort. The Kaplan–Meier method and log-rank test demonstrated that high ATGPI was significantly associated with short OS of patients with early-stage LUAD in the UH cohort (P = 0.038) (Figure 7A). The RMS curves of continuous ATGPI and ACPI in the UH cohort demonstrated that significant improvement in predicting OS was achieved by ACPI (C-index: 0.70, 95% CI, 0.54–0.86 vs. C-index: 0.80, 95% CI, 0.68–0.92; P = 0.015) (Figure 7B). IHC analysis demonstrated that the proportions of tumor-infiltrating CD4+ T cells were higher in the low-ATGPI group that in the high-ATGPI group (Figure 7C).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/719011/fcell-09-719011-HTML/image_m/fcell-09-719011-g007.jpg)
Figure 7. The Kaplan–Meier curves for the autophagy-related gene pair index (ATGPI) and restricted mean survival (RMS) curves for the ATGPI and autophagy clinical prognostic index (ACPI) in the Union hospital (UH) cohort. (A) The Kaplan–Meier curves showed that patients were stratified into different prognostic groups by the ATGPI. (B) The RMS curves revealed that a superior estimation for overall survival was achieved by the ACPI in the UH cohort. (C) Immunohistochemistry (IHC) analysis showed the proportions of tumor-infiltrating CD4+ T cells were higher in the low-ATGPI group.
Discussion
Complete surgical resection of early-stage LUAD significantly improves prognosis, but some patients are still at a high risk of recurrence and death. Previous studies have not demonstrated a survival benefit of adjuvant therapy in patients with early-stage LUAD. To identify high-risk early-stage LUAD, clinicians urgently need reliable prognostic biomarkers. In this study, we first constructed a TIL-related prognostic biomarker based on ATGPs for patients with early-stage LUAD by combining multiple gene expression datasets. To construct a robust signature, data from 12 GEO LUAD datasets were randomly divided into training and testing datasets. TCGA and UH cohorts were included as independent validation datasets. The optimal prediction model was established using the Cox proportional hazards regression model with the LASSO in the training cohort and was validated in independent cohorts. Moreover, we integrated ATG signatures with clinical characteristics to build the ACPI, which improved the estimation of prognosis for patients with early-stage LUAD. The expression levels of genes in each sample were compared in pairs to generate the ATGPI. This gene pair-based approach was based on the relative gene expression of each sample and did not require normalization.
Autophagy is a dynamic process that captures and degrades intracellular components such as proteins and organelles to sustain metabolism and homeostasis (Mizushima and Komatsu, 2011). This survival-promoting pathway mediated by evolutionarily conserved ATG genes, recycles intracellular materials in lysosomes, prevents the accumulation of toxic cellular wastes, and involves multiple signaling pathways, such as the RAS, p53/DRAM, PI3K/AKT/mTOR, and JAK-STAT signaling pathways (Levy et al., 2017). Autophagy plays an important role in lung cancer. For example, PAQR3 enhances autophagy induced by the EGFR inhibitor erlotinib to inhibit the growth of NSCLC cells, and a combination therapy, including lipophilic bisphosphonates, facilitates autophagy, and prevents KRAS-driven LUAD cell proliferation (Xia et al., 2014). Because of the important role of autophagy in lung cancer, including ATG gene signatures in survival prediction for lung cancer is crucial.
To construct a robust ATG signature, we identified ATG genes from The Human Autophagy Database and profiled the mRNA expression of 142 ATG genes that were detected in all datasets used in this study. The final ATGPI was constructed using 10 ATGPs consisting of 15 ATG genes (ARSB, BAK1, BAX, BCL2, BIRC5, CCL2, CCR2, CDKN2A, FADD, CX3CL1, SPHK1, MAP2K7, TSC1, MAPK9, and HSPB8). Dysregulation of the expression of anti-apoptotic BCL-2 family proteins is observed in multiple cancers. Combining BCL-2/BCL-XL inhibitors with kinase inhibitors enhanced the efficacy of EGFR inhibitors against EGFR mutant LUAD (Gong et al., 2007). Lung cancer cells can induce macrophage infiltration by increasing the production of CCL2 and CXCL3 (Schmall et al., 2015). Inhibition of the CCL2-CCR2 signaling pathway blocks the recruitment of inflammatory monocytes and prolongs survival in a mouse model (Qian et al., 2011). Increased expression of FADD is associated with adverse clinical outcomes in patients with LUAD (Chen et al., 2005). SPHK1 regulates tumor cell apoptosis and promotes NSCLC development (Ma et al., 2021). Loss of TSC1 synergizes with KRAS mutation to enhance the development of lung cancer in mice (Liang et al., 2010). Inactivation of MAP2K7 in KRAS-driven lung cancer accelerates tumorigenesis and reduces survival (Schramek et al., 2011). The autophagy process in the NSCLC microenvironment can be modulated by variable immune signals, such as extracellular damage-associated molecular patterns and cytokines. For example, the activation of Toll-like receptors (TLR3 and TLR4) with polyinosinic:polycytidylic acid and lipopolysaccharide can induce autophagy and upregulate metastasis-related cytokines (Zhan et al., 2014). Few studies focused on the influence of the relative expression of ATGPs on cancer progression. For example, overexpression of SPHK1 induced the expression of BCL2 to promote the proliferation and migration of NSCLC cells, suggesting that the expression of SPHK1 is positively associated with the expression of BCL2 (Ma et al., 2021). Additionally, the coexpression of CCL2 and CX3CL1 was strongly associated with prognosis of patients with soft tissue sarcoma, particularly female patients (Kehlen et al., 2014). Further studies are needed to explore the potential value of the relative expression of ATGP signatures. Bioinformatic and IHC analyses indicated that the proportions of CD4+ TILs were lower in the high-ATGPI group than in the low-ATGPI group. The increased infiltration of CD4+ TILs was associated with better prognosis of patients with esophageal squamous cell carcinoma (Chen et al., 2017) and triple-negative breast cancer (Gao et al., 2020). Furthermore, CD4+ TILs demonstrated a functional heterogeneity in NSCLC (Oja et al., 2018). The vital role of these ATG genes in the development and progression of NSCLC and the association between ATGPI and TILs might explain the robust prediction ability of ATGPI and ACPI.
Limitations
Although batch effects were reduced by the pairwise comparison of ATG genes, some were unavoidable because of the combination of different platforms (Leek et al., 2010). Additionally, sampling bias caused by intratumor genetic heterogeneity was inevitable for gene expression-based signatures.
In summary, ATGPI is a prospective prognostic biomarker for LUAD, especially early-stage LUAD. Combining ATGPI with clinical characteristics improved the accuracy of the model. Further studies are needed to validate the prediction accuracy of ATGPI and ACPI and explore the clinical utility of ATGPI and ACPI in the individualized management of LUAD.
Data Availability Statement
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.
Ethics Statement
The studies involving human participants were reviewed and approved by the Ethics Committees of Tongji Medical College, Huazhong University of Science and Technology (2019-S671). The patients/participants provided their written informed consent to participate in this study. Written informed consent was obtained from the individual(s) for the publication of any potentially identifiable images or data included in this article.
Author Contributions
QZ conceived the idea, designed, and supervised the study, had full access to all data and took responsibility for the integrity of the data. X-SW, YL, L-LY, W-BP, S-YZ, and Y-RN recorded and sorted the data. XX and Q-QX analyzed the data. Z-HW and PZ interpreted the data and wrote the manuscript. All authors contributed significantly to this work, agreed to be accountable for the work, and read and approved the final manuscript.
Funding
This work was supported in part by grants from National Natural Science Foundation of China (Nos. 81973990, 81900096, and 81770090).
Conflict of Interest
The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s Note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
Acknowledgments
We would like to thank Shi-Peng Guo and Zhen-Hua Zhu for the bioinformatic skills. We would also like to thank Editage for English language editing.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.719011/full#supplementary-material
Abbreviations
ATG, autophagy-related; ATGPs, autophagy-related gene pairs; ATGPI, autophagy-related gene pair index; ACPI, autophagy clinical prognostic index; AUC, area under the curve; C-index, concordance index; FP, false positive rate; GEO, gene expression omnibus; LASSO, least absolute shrinkage and selection operator; LUAD, lung adenocarcinoma; NSCLC, non-small cell lung cancer; OS, overall survival; ROC, receiver operating characteristics; TCGA, The Cancer Genome Atlas; TP, true positive rate; TILs, tumor-infiltrating lymphocytes.
Footnotes
- ^ https://portal.gdc.cancer.gov/
- ^ https://www.ncbi.nlm.nih.gov/geo/
- ^ http://www.autophagy.lu/
- ^ http://bioinfo.life.hust.edu.cn/ImmuCellAI#!/
References
Ajona, D., Pajares, M. J., Corrales, L., Perez-Gracia, J. L., Agorreta, J., Lozano, M. D., et al. (2013). Investigation of complement activation product c4d as a diagnostic and prognostic biomarker for lung cancer. J. Natl. Cancer Inst. 105, 1385–1393. doi: 10.1093/jnci/djt205
Amaravadi, R., Kimmelman, A. C., and White, E. (2016). Recent insights into the function of autophagy in cancer. Genes Dev. 30, 1913–1930. doi: 10.1101/gad.287524.116
Beer, D. G., Kardia, S. L., Huang, C. C., Giordano, T. J., Levin, A. M., Misek, D. E., et al. (2002). Gene-expression profiles predict survival of patients with lung adenocarcinoma. Nat. Med. 8, 816–824. doi: 10.1038/nm733
Bhattacharjee, A., Richards, W. G., Staunton, J., Li, C., Monti, S., Vasa, P., et al. (2001). Classification of human lung carcinomas by mRNA expression profiling reveals distinct adenocarcinoma subclasses. Proc. Natl. Acad. Sci. U.S.A. 98, 13790–13795. doi: 10.1073/pnas.191502998
Chen, G., Bhojani, M. S., Heaford, A. C., Chang, D. C., Laxman, B., Thomas, D. G., et al. (2005). Phosphorylated FADD induces NF-kappaB, perturbs cell cycle, and is associated with poor outcome in lung adenocarcinomas. Proc. Natl. Acad. Sci. U.S.A. 102, 12507–12512. doi: 10.1073/pnas.0500397102
Chen, K., Zhu, Z., Zhang, N., Cheng, G., Zhang, F., Jin, J., et al. (2017). Tumor-infiltrating CD4+ lymphocytes predict a favorable survival in patients with operable esophageal squamous cell carcinoma. Med. Sci. Monit. 23, 4619–4632. doi: 10.12659/msm.904154
Chen, X., Mao, R., Su, W., Yang, X., Geng, Q., Guo, C., et al. (2020). Circular RNA circHIPK3 modulates autophagy via MIR124-3p-STAT3-PRKAA/AMPKα signaling in STK11 mutant lung cancer. Autophagy 16, 659–671. doi: 10.1080/15548627.2019.1634945
Der, S. D., Sykes, J., Pintilie, M., Zhu, C. Q., Strumpf, D., Liu, N., et al. (2014). Validation of a histology-independent prognostic gene signature for early-stage, non-small-cell lung cancer including stage IA patients. J. Thorac. Oncol. 9, 59–64. doi: 10.1097/jto.0000000000000042
Dong, S., Li, W., Wang, L., Hu, J., Song, Y., Zhang, B., et al. (2019). Histone-related genes are hypermethylated in lung cancer and hypermethylated HIST1H4F could serve as a pan-cancer biomarker. Cancer Res. 79, 6101–6112. doi: 10.1158/0008-5472.Can-19-1019
Gao, G., Wang, Z., Qu, X., and Zhang, Z. (2020). Prognostic value of tumor-infiltrating lymphocytes in patients with triple-negative breast cancer: a systematic review and meta-analysis. BMC Cancer 20:179. doi: 10.1186/s12885-020-6668-z
Gautier, L., Cope, L., Bolstad, B. M., and Irizarry, R. A. (2004). affy–analysis of Affymetrix GeneChip data at the probe level. Bioinformatics 20, 307–315. doi: 10.1093/bioinformatics/btg405
Gong, Y., Somwar, R., Politi, K., Balak, M., Chmielecki, J., Jiang, X., et al. (2007). Induction of BIM is essential for apoptosis triggered by EGFR kinase inhibitors in mutant EGFR-dependent lung adenocarcinomas. PLoS Med. 4:e294. doi: 10.1371/journal.pmed.0040294
Guo, J. Y., and White, E. (2013). Autophagy is required for mitochondrial function, lipid metabolism, growth, and fate of KRAS(G12D)-driven lung tumors. Autophagy 9, 1636–1638. doi: 10.4161/auto.26123
Hao, X., Luo, H., Krawczyk, M., Wei, W., Wang, W., Wang, J., et al. (2017). DNA methylation markers for diagnosis and prognosis of common cancers. Proc. Natl. Acad. Sci. U.S.A. 114, 7414–7419. doi: 10.1073/pnas.1703577114
Kang, L., Chen, W., Petrick, N. A., and Gallas, B. D. (2015). Comparing two correlated C indices with right-censored survival outcome: a one-shot nonparametric approach. Stat. Med. 34, 685–703. doi: 10.1002/sim.6370
Kehlen, A., Greither, T., Wach, S., Nolte, E., Kappler, M., Bache, M., et al. (2014). High coexpression of CCL2 and CX3CL1 is gender-specifically associated with good prognosis in soft tissue sarcoma patients. Int. J. Cancer 135, 2096–2106. doi: 10.1002/ijc.28867
Kimmelman, A. C., and White, E. (2017). Autophagy and tumor metabolism. Cell Metab. 25, 1037–1043. doi: 10.1016/j.cmet.2017.04.004
Kuner, R., Muley, T., Meister, M., Ruschhaupt, M., Buness, A., Xu, E. C., et al. (2009). Global gene expression analysis reveals specific patterns of cell junctions in non-small cell lung cancer subtypes. Lung Cancer 63, 32–38. doi: 10.1016/j.lungcan.2008.03.033
Leek, J. T., Scharpf, R. B., Bravo, H. C., Simcha, D., Langmead, B., Johnson, W. E., et al. (2010). Tackling the widespread and critical impact of batch effects in high-throughput data. Nat. Rev. Genet. 11, 733–739. doi: 10.1038/nrg2825
Levine, B., and Kroemer, G. (2008). Autophagy in the pathogenesis of disease. Cell 132, 27–42. doi: 10.1016/j.cell.2007.12.018
Levy, J. M. M., Towers, C. G., and Thorburn, A. (2017). Targeting autophagy in cancer. Nat. Rev. Cancer 17, 528–542. doi: 10.1038/nrc.2017.53
Levy, J. M., and Thorburn, A. (2011). Targeting autophagy during cancer therapy to improve clinical outcomes. Pharmacol. Ther. 131, 130–141. doi: 10.1016/j.pharmthera.2011.03.009
Li, B., Cui, Y., Diehn, M., and Li, R. (2017). Development and validation of an individualized immune prognostic signature in early-stage nonsquamous non-small cell lung cancer. JAMA Oncol. 3, 1529–1537. doi: 10.1001/jamaoncol.2017.1609
Li, Y., Sun, N., Lu, Z., Sun, S., Huang, J., Chen, Z., et al. (2017). Prognostic alternative mRNA splicing signature in non-small cell lung cancer. Cancer Lett. 393, 40–51. doi: 10.1016/j.canlet.2017.02.016
Liang, M. C., Ma, J., Chen, L., Kozlowski, P., Qin, W., Li, D., et al. (2010). TSC1 loss synergizes with KRAS activation in lung cancer development in the mouse and confers rapamycin sensitivity. Oncogene 29, 1588–1597. doi: 10.1038/onc.2009.452
Ma, Y., Xing, X., Kong, R., Cheng, C., Li, S., Yang, X., et al. (2021). SphK1 promotes development of non-small cell lung cancer through activation of STAT3. Int. J. Mol. Med. 47, 374–386. doi: 10.3892/ijmm.2020.4796
Martínez-Terroba, E., Behrens, C., de Miguel, F. J., Agorreta, J., Monsó, E., Millares, L., et al. (2018). A novel protein-based prognostic signature improves risk stratification to guide clinical management in early-stage lung adenocarcinoma patients. J. Pathol. 245, 421–432. doi: 10.1002/path.5096
Miao, Y. R., Zhang, Q., Lei, Q., Luo, M., Xie, G. Y., Wang, H., et al. (2020). ImmuCellAI: a unique method for comprehensive T-cell subsets abundance prediction and its application in cancer immunotherapy. Adv. Sci. (Weinh) 7:1902880. doi: 10.1002/advs.201902880
Mizushima, N., and Komatsu, M. (2011). Autophagy: renovation of cells and tissues. Cell 147, 728–741. doi: 10.1016/j.cell.2011.10.026
Oja, A. E., Piet, B., van der Zwan, D., Blaauwgeers, H., Mensink, M., de Kivit, S., et al. (2018). Functional heterogeneity of CD4(+) tumor-infiltrating lymphocytes with a resident memory phenotype in NSCLC. Front. Immunol. 9:2654. doi: 10.3389/fimmu.2018.02654
Okayama, H., Kohno, T., Ishii, Y., Shimada, Y., Shiraishi, K., Iwakawa, R., et al. (2012). Identification of genes upregulated in ALK-positive and EGFR/KRAS/ALK-negative lung adenocarcinomas. Cancer Res. 72, 100–111. doi: 10.1158/0008-5472.Can-11-1403
Ooki, A., Maleki, Z., Tsay, J. J., Goparaju, C., Brait, M., Turaga, N., et al. (2017). A panel of novel detection and prognostic methylated DNA markers in primary non-small cell lung cancer and serum DNA. Clin. Cancer Res. 23, 7141–7152. doi: 10.1158/1078-0432.Ccr-17-1222
Qian, B. Z., Li, J., Zhang, H., Kitamura, T., Zhang, J., Campion, L. R., et al. (2011). CCL2 recruits inflammatory monocytes to facilitate breast-tumour metastasis. Nature 475, 222–225. doi: 10.1038/nature10138
Rao, S., Yang, H., Penninger, J. M., and Kroemer, G. (2014). Autophagy in non-small cell lung carcinogenesis: a positive regulator of antitumor immunosurveillance. Autophagy 10, 529–531. doi: 10.4161/auto.27643
Rousseaux, S., Debernardi, A., Jacquiau, B., Vitte, A. L., Vesin, A., Nagy-Mignotte, H., et al. (2013). Ectopic activation of germline and placental genes identifies aggressive metastasis-prone lung cancers. Sci. Transl. Med. 5:186ra166. doi: 10.1126/scitranslmed.3005723
Sato, M., Larsen, J. E., Lee, W., Sun, H., Shames, D. S., Dalvi, M. P., et al. (2013). Human lung epithelial cells progressed to malignancy through specific oncogenic manipulations. Mol. Cancer Res. 11, 638–650. doi: 10.1158/1541-7786.Mcr-12-0634-t
Schmall, A., Al-Tamari, H. M., Herold, S., Kampschulte, M., Weigert, A., Wietelmann, A., et al. (2015). Macrophage and cancer cell cross-talk via CCR2 and CX3CR1 is a fundamental mechanism driving lung cancer. Am. J. Respir. Crit. Care Med. 191, 437–447. doi: 10.1164/rccm.201406-1137OC
Schramek, D., Kotsinas, A., Meixner, A., Wada, T., Elling, U., Pospisilik, J. A., et al. (2011). The stress kinase MKK7 couples oncogenic stress to p53 stability and tumor suppression. Nat. Genet. 43, 212–219. doi: 10.1038/ng.767
Shen, W., Zhang, W., Ye, W., Wang, H., Zhang, Q., Shen, J., et al. (2020). SR9009 induces a REV-ERB dependent anti-small-cell lung cancer effect through inhibition of autophagy. Theranostics 10, 4466–4480. doi: 10.7150/thno.42478
Stafford, P., and Brun, M. (2007). Three methods for optimization of cross-laboratory and cross-platform microarray expression data. Nucleic Acids Res. 35:e72. doi: 10.1093/nar/gkl1133
Takeuchi, T., Tomida, S., Yatabe, Y., Kosaka, T., Osada, H., Yanagisawa, K., et al. (2006). Expression profile-defined classification of lung adenocarcinoma shows close relationship with underlying major genetic changes and clinicopathologic behaviors. J. Clin. Oncol. 24, 1679–1688. doi: 10.1200/jco.2005.03.8224
Tang, H., Xiao, G., Behrens, C., Schiller, J., Allen, J., Chow, C. W., et al. (2013). A 12-gene set predicts survival benefits from adjuvant chemotherapy in non-small cell lung cancer patients. Clin. Cancer Res. 19, 1577–1586. doi: 10.1158/1078-0432.Ccr-12-2321
The Cancer Genome Atlas Research Network (2014). Comprehensive molecular profiling of lung adenocarcinoma. Nature 511, 543–550. doi: 10.1038/nature13385
Tomida, S., Takeuchi, T., Shimada, Y., Arima, C., Matsuo, K., Mitsudomi, T., et al. (2009). Relapse-related molecular signature in lung adenocarcinomas identifies patients with dismal prognosis. J. Clin. Oncol. 27, 2793–2799. doi: 10.1200/jco.2008.19.7053
Torre, L. A., Bray, F., Siegel, R. L., Ferlay, J., Lortet-Tieulent, J., and Jemal, A. (2015). Global cancer statistics, 2012. CA Cancer J. Clin. 65, 87–108. doi: 10.3322/caac.21262
Tsang, T., Posimo, J. M., Gudiel, A. A., Cicchini, M., Feldser, D. M., and Brady, D. C. (2020). Copper is an essential regulator of the autophagic kinases ULK1/2 to drive lung adenocarcinoma. Nat. Cell Biol. 22, 412–424. doi: 10.1038/s41556-020-0481-4
Wilkerson, M. D., Yin, X., Walter, V., Zhao, N., Cabanski, C. R., Hayward, M. C., et al. (2012). Differential pathogenesis of lung adenocarcinoma subtypes involving sequence mutations, copy number, chromosomal instability, and methylation. PLoS One 7:e36530. doi: 10.1371/journal.pone.0036530
Xia, Y., Liu, Y. L., Xie, Y., Zhu, W., Guerra, F., Shen, S., et al. (2014). A combination therapy for KRAS-driven lung adenocarcinomas using lipophilic bisphosphonates and rapamycin. Sci. Transl. Med. 6:263ra161. doi: 10.1126/scitranslmed.3010382
Keywords: autophagy, autophagy-related genes, gene signatures, lung adenocarcinoma, prognosis, tumor-infiltrating lymphocytes
Citation: Wang Z-H, Li Y, Zhang P, Xiang X, Wei X-S, Niu Y-R, Ye L-L, Peng W-B, Zhang S-Y, Xue Q-Q and Zhou Q (2021) Development and Validation of a Prognostic Autophagy-Related Gene Pair Index Related to Tumor-Infiltrating Lymphocytes in Early-Stage Lung Adenocarcinoma. Front. Cell Dev. Biol. 9:719011. doi: 10.3389/fcell.2021.719011
Received: 01 June 2021; Accepted: 30 July 2021;
Published: 20 September 2021.
Edited by:
Guillermo Velasco, Complutense University of Madrid, SpainReviewed by:
Emenike K. Onyido, Swansea University Medical School, United KingdomMarina Mendiburu-Eliçabe, Centro de Investigación del Cancer, Spain
Copyright © 2021 Wang, Li, Zhang, Xiang, Wei, Niu, Ye, Peng, Zhang, Xue and Zhou. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Qiong Zhou, zhouqiongtj@126.com
†These authors have contributed equally to this work