- 1Department of Thoracic Surgery, National Cancer Center, National Clinical Research Center for Cancer, Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
- 2State Key Laboratory of Molecular Oncology, National Cancer Center, National Clinical Research Center for Cancer, Cancer Hospital, Chinese Academy of Medical Sciences and Peking Union Medical College, Beijing, China
The iron-dependent cell death named ferroptosis has been implicated in the progression and therapeutic response of several tumors. However, potential role of ferroptosis in lung adenocarcinomas (LUAD) remained less well understood. In TCGA-LUAD cohort, unsupervised clustering was first conducted based on ferroptosis regulators extracted from FerrDb database. Comprehensive correlation analysis and comparisons were performed among ferroptosis subtypes. The ferroptosis-related prognostic (FRP) signature was identified based on filtered features and repeated LASSO and was validated in five independent cohorts. The clinical relevance between the risk score and therapeutic response was further explored by multiple algorithms. qPCR was implemented to verify gene expression. A total of 1,168 LUAD patients and 161 ferroptosis regulators were included in this study. Three ferroptosis subtypes were identified and patients in subtype B had the best prognosis among the three subtypes. Significant differences in immune microenvironment and biological function enrichment were illustrated in distinct subtypes. The Boruta algorithm was conducted on 308 common differentially expressed genes for dimensionality reduction. A total of 56 genes served as input for model construction and a six-gene signature with the highest frequencies of 881 was chosen as FRP. The prognostic significance of FRP was validated in five independent cohorts. High FRP risk score was also linked to increased tumor mutation burden, PD-L1 protein expression and number of neoantigens. Of the FRP genes, 83.3% was abnormally expressed in LUAD cell lines. In conclusion, ferroptosis plays a non-negligible role in LUAD. Exploration of the ferroptosis pattern will enhance the prognostic stratification of individual patients and move toward the purpose of personalized treatment.
Introduction
Ferroptosis is a non-apoptotic modality of regulated cell death (RCD) that is induced by excessive phospholipid peroxidation that progresses in an iron-dependent manner (Dixon et al., 2012). The biology of ferroptotic cells is morphologically and mechanistically distinct from that of apoptotic cells, as indicated by ferroptosis induction in BAX/BAK-deficient cells and the lack of detected caspase-3 cleavage during ferroptosis (Stockwell et al., 2017). Extrinsic and intrinsic pathways that result in unrestricted phospholipid peroxidation can trigger ferroptosis (Tang and Kroemer, 2020). The extrinsic pathway includes blockade of the system xc– cystine/glutamate antiporter, which is the main transporter of cystine, the rate-limiting precursor of glutathione, and activator of the iron transporters serotransferrin and lactotransferrin, while the intrinsic pathway is initiated by the inhibition of intracellular antioxidant enzymes, such as glutathione peroxidase GPX4. Current studies have identified ferroptosis as a critical tumor suppression mechanism, and inhibition of ferroptosis, similar to inhibition of apoptosis, promotes tumor development (Gao et al., 2019; Liu et al., 2019). In addition, multiple tumor suppressors can sensitize cancer cells to ferroptosis, and among these suppressors, p53 has been thoroughly investigated (Zhang et al., 2018, 2019; Jiang et al., 2021). p53 can inhibit the transcription of the systemic xc– subunit SLC7A11 (Jiang et al., 2015; Wang et al., 2016), and P47S amino acid mutation of p53 was reported to confer resistance to ferroptosis, thus promoting tumor formation (Jennis et al., 2016).
A number of immune checkpoint inhibitors (ICIs) have been developed and utilized in cancer patients in recent decades (He and Xu, 2020), among which anti-PD-1/PD-L1 therapy has led to great success and has been approved to treat a wide variety of cancer types and has led to intermediate response rates (Ribas and Wolchok, 2018). Anti-PD-1/PD-L1 therapy reactivates the immune system to exert antitumor effects by relieving the suppressive functions of PD-1 signaling (Ribas, 2015). Mechanistically, the antitumor effect of ICIs is strongly associated with the host’s immune system and tumor immune microenvironment (TIME). The interaction between cancer cells and the TIME is complicated due to the coexistence of both protumorigenic and antitumorigenic factors. Induction of the T cell-inflamed TIME induced by JHU083 (Leone et al., 2019), CDK4/6 (Schaer et al., 2018), DCR-BCAT (Ganesh et al., 2018), or other treatments has been demonstrated to enhance the efficacy of ICIs. Thus, deeper insight into the TIME is needed to facilitate the identification of novel predictive biomarkers and the development of novel therapeutic strategies.
Recently, several studies have explored the specific association between ferroptosis and TIME infiltrating immune cells. It has been reported that anti-PD-L1 and ferroptosis inducers [such as erastin, RSL3, and cyst(e)inase] synergistically inhibit tumor development both in vivo and in vitro (Wang et al., 2019). Mechanistically, interferon-γ (IFNγ) released from activated cytotoxic T cells downregulates the expression of SCL7A11 and SLC3A2 through activation of the JAK-STAT1 pathway and thus induces ferroptosis of cancer cells. However, the cellular components released by ferroptotic cancer cells play dual roles in antitumor immunity. KRAS-G12D within exosomes secreted by ferroptotic pancreatic cancer cells can be absorbed by macrophages through AGER, which results in the polarization of macrophages to the M2 phenotype and an immunosuppressive microenvironment for tumor growth (Dai et al., 2020a). In addition, ferroptosis induced by conditional depletion of Gpx4 in the pancreas promotes mutant Kras-driven tumorigenesis in mice through secretion of DNA and activation of the TMEM173/STING-dependent DNA sensor pathway in macrophages (Dai et al., 2020b).
Although success in revealing the molecular mechanisms between ferroptosis and TIME has been achieved, the abovementioned studies focused on limited ferroptosis regulators and were mainly conducted in vitro. Therefore, comprehensive analysis of ferroptosis regulators and TIME in real-world cohorts of patients is urgently needed. In this study, we integrated transcriptome datasets of 500 lung adenocarcinoma (LUAD) samples from The Cancer Genome Atlas (TCGA) and found three distinct ferroptosis subtypes, among which significant survival differences, diverse biological pathways, and different infiltrated immune cells were observed. Furthermore, we constructed a robust ferroptosis-related prognostic (FRP) signature composed of 6 genes (EIF5A, CACYBP, CYCS, ANLN, ARNTL2, andPPM1M) based on defined ferroptosis subtypes and confirmed the FRP signature as an independent prognostic factor for LUAD patients. Notably, patients in the high-risk FRP signature group may potentially benefit from immunotherapy due to the elevation of tumor mutation burden (TMB) and PD-L1 and neoantigen expression. In conclusion, our results revealed that ferroptosis plays critical roles in LUAD and that the FRP signature may be a promising predictive biomarker for prognosis and immunotherapy response.
Materials and Methods
Lung Adenocarcinoma Dataset Acquisition and Preprocessing
Public lung adenocarcinoma transcriptome datasets were obtained from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO). Patients without prognostic information were excluded from the study. A total of six datasets (TCGA-LUAD, GSE13213, GSE30219, GSE31210, GSE3141, and GSE41271) were gathered for further analysis (Bild et al., 2006; Tomida et al., 2009; Okayama et al., 2012; Rousseaux et al., 2013; Sato et al., 2013). For RNA-Seq data from the TCGA, FPKM values were downloaded from the TCGA Data Portal1 and transformed into TPM values before analysis. For datasets in the GEO, normalized probe expression values generated from the microarray were acquired using the R package GEOquery (Davis and Meltzer, 2007). The probes in each dataset were annotated according to the appropriate platform annotation file. When multiple probes mapped to the same gene, the average expression values of these genes were used as the final expression value. To calculate tumor mutation burden (TMB) in TCGA-LUAD, gene-level mutation (VarScan) data were extracted from the TCGA Data Portal. In the present study, TMB was defined as the total number of non-synonymous mutations divided by the size of the exome.
Consensus Clustering of Ferroptosis Regulators
A comprehensive list of ferroptosis drivers and suppressors was acquired from the FerrDb database2 (Zhou and Bao, 2020). After removing duplicated entries and performing univariate Cox regression for each gene, 37 unique prognosis-associated genes were obtained for unsupervised clustering. The identification of ferroptosis subtypes was conducted with the TCGA cohort by applying the partitioning around medoids cluster algorithm in the R package ConsensusClusterPlus (Wilkerson and Hayes, 2010). The number of repetitions was set to 100 to ensure the stability of the result. The other parameters of the cluster algorithm were set to the default values.
Gene Set Variation Analysis and Functional Comparison Between Subtypes
To determine the biological function associated with different ferroptosis subtypes, GSVA was performed to detect subtle pathway activity in each sample using the R package GSVA (Hänzelmann et al., 2013). GSVA is a non-parametric, unsupervised method for calculating gene set enrichment through expression profiles. The gene set file named “c2.cp.kegg.v7.4.symbols.gmt,” including 186 canonical pathways, was downloaded from MSigDB and was chosen for pathway enrichment (Liberzon et al., 2011). A pathway comparison was conducted using the R package limma (Smyth, 2005), and an adjusted p value of less than 0.05 was considered statistically significant.
Exploration of Tumor-Infiltrating Immune Cells in the TCGA-Lung Adenocarcinomas Cohort
Two widely validated algorithms, CIBERSORT and ImmuCellAI (Newman et al., 2015; Miao et al., 2020), were implemented to estimate the abundance of specific immune components. CIBERSORT is a computational approach reinforced by support vector regression to quantify 22 immune cell fractions. ImmuCellAI was used for abundance estimation of 24 immune cell types, including 18 T cell subsets, from gene expression matrices. The absolute fraction and abundance of immune cells were estimated by CIBERSORT and ImmuCellAI, respectively. Only intergroup comparisons were conducted and the differences in immune cell abundance were depicted between risk groups and among subtypes.
Construction of a Ferroptosis-Related Prognostic Model
To identify common DEGs among the three ferroptosis subtypes, gene expression was quantified by TPM value and empirical Bayes statistics was applied for differential expression analysis using limma package (Smyth, 2005). The obtained p-values were further adjusted by the Benjamini–Hochberg method. The DEGs were defined as genes with an adjusted p-value less than 0.001. By taking the intersection of the outputs, common DEGs of the three subtypes were obtained and illustrated using a Venn diagram. Next, the Boruta feature selection algorithm was applied to narrow the genes for model construction using the R package Boruta (Kursa and Rudnicki, 2010). Principal component analysis (PCA) was p to verify the ability to distinguish subtypes using selected genes. Furthermore, genes related to prognosis were screened out through univariate Cox regression and served as input for model training. To generate a robust ferroptosis-related prognostic (FRP) signature, LASSO Cox regression with 10-fold cross-validation was performed using a random seed, and the whole procedure was repeated 1,000 times (Tibshirani, 1997). The model with the highest frequency was selected as the final signature. Patients were subsequently stratified into high-risk and low-risk groups according to the median of the risk score. Both Kaplan–Meier curves and multivariate Cox regression were performed in the training cohort.
Evaluation of Ferroptosis Activity in the Different Risk Groups
To quantify the activity of ferroptosis, an established indicator named the ferroptosis potential index (FPI) was obtained from a previous study and calculated for each sample based on gene expression data (Liu et al., 2020). During computation, key ferroptosis regulator genes were first divided into positive and negative groups according to their potential function. Then, the enrichment score of these two groups was calculated separately using the R package GSVA (Hänzelmann et al., 2013). The FPI was defined as the enrichment score of the positive group minus the negative group and was considered a reliable indicator of ferroptosis level.
Validation of the Ferroptosis-Related Prognostic Signature in Multiple Validation Cohorts
To further validate the prognostic value of the FRP signature, the risk score was first calculated for each sample based on the formula derived from the training cohort. Then, samples were divided into low-risk and high-risk groups according to the median risk score in each cohort. Kaplan–Meier curves were generated, and univariate Cox regression was performed, and the concordance index (C-index) was applied to evaluate the performance of the FRP signature.
Exploring the Implications of the Ferroptosis-Related Prognostic Signature on Immunotherapeutic and Chemotherapeutic Responses
To infer the potential benefit of immunotherapy for each patient, the number of neoantigens and the PD-L1 protein expression of patients in the TCGA-LUAD cohort were extracted from The Cancer Immunome Atlas3 and cBioPortal for Cancer Genomics4, respectively (Gao et al., 2013; Charoentong et al., 2017). To determine the chemotherapeutic implications, experimental data for cancer cell drug sensitivity were obtained from the Genomics of Drug Sensitivity in Cancer (GDSC), and the whole prediction process was conducted with the R package pRRophetic (Geeleher et al., 2014). “allSoldTumours” was selected for the tissue type parameter, and the remaining parameters were set to the default values.
Cell Lines and Cell Culture Conditions
BEAS-2B and HCC827 cell lines were purchased from ATCC. BEAS-2B cells were cultured in BEBM basal medium supplemented with 10% fetal bovine serum (FBS), 1,000 U ml–1 penicillin and 100 μg ml–1 streptomycin. HCC827 cells were maintained in RPMI 1640 supplemented with 10% FBS, 1,000 U ml–1 penicillin and 100 μg ml–1 streptomycin. The cells were cultured with 5% CO2 at 37°C in a humidified incubator.
RNA Extraction and Quantitative RT-PCR Analysis
Total RNA was extracted from cells using TRIzol reagent (Invitrogen, United States) according to the manufacturer’s instructions. Equal amounts of RNA samples were used for cDNA synthesis with a RevertAid First Strand cDNA synthesis kit (Thermo). SYBR Green-based qRT-PCR on a 7900HT fast real-time PCR system (Applied Biosystems/Life Technologies, Waltham, United States) was performed. The relative mRNA expression levels were calculated by the 2–ΔΔCt method with normalization to ACTB; the PCR primers are listed in Supplementary Table 1.
Statistical Analysis
All bioinformatic analyses were conducted in R software (Version 3.6.3). Wilcoxon test and Kruskal–Wallis test were performed to statistically test the difference in continuous data for defined subtypes. The Spearman coefficient was used to evaluate the association between variables. Survival analyses, including Kaplan–Meier curves and Cox regression, were performed using the R package survival. The C-index was calculated by R package survcomp. A P-value less than 0.05 was considered statistically significant.
Results
Identification of Distinct Ferroptosis Subtypes on the Basis of Ferroptosis Regulators
The flow chart of this study is shown in Figure 1. A comprehensive list of genes including 108 that encode ferroptosis drivers and 69 that encode suppressors was extracted from the FerrDb. After removing multirole regulators that would be repeatedly counted, a total of 161 genes remained for screening. Univariate Cox regression was performed, and 37 genes with prognostic significance (p < 0.05) were chosen from the list. Based on the expression level of these genes, three subtypes were identified in the TCGA-LUAD cohort. The clinical information and gene expression patterns are displayed in Figure 1. In subtype A, most of the genes, including both drivers and suppressors, were highly expressed, and in contrast, a substantial portion of ferroptosis regulators was expressed at low levels in subtype B (Figure 2A). The expression pattern in subtype C was somewhere in between the other two subtypes. An association analysis suggested that there was a close connection between the expression levels of the 37 ferroptosis regulators (Figure 2B). Of 24 ferroptosis drivers, nearly one-third (7/24) of the genes was revealed to be favorable prognostic factors, and the remaining genes (17/24) were risk factors for overall survival (OS; Figure 2B). All 13 ferroptosis suppressors were unfavorable prognostic factors (Figure 2B). These results suggested that inhibition of ferroptosis suppressors may improve the survival of LUAD patients. In addition, a significant survival difference (p < 0.0001) was observed among the three ferroptosis subtypes (Figure 2C). Patients with subtype B had the best prognosis compared with patients with the other subtypes, while those with subtype C had the worst prognosis (Figure 2C). The clinicopathological characteristics of the three subtypes were summarized in Table 1.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/743724/fcell-09-743724-HTML/image_m/fcell-09-743724-g002.jpg)
Figure 2. Unsupervised clustering of ferroptosis regulators and characteristics of ferroptosis subtypes. (A) Heat map showing the expression patterns of 37 ferroptosis regulator genes in the TCGA-LUAD cohort. Sex, tumor location, stage, smoking history, age, and survival status were the annotations. (B) Interactions of 37 ferroptosis regulator genes. The size of each circle indicates the survival significance of each gene. Ferroptosis drivers are shown in blue, and ferroptosis suppressors are shown in red. Favorable prognostic factors are indicated in black, and unfavorable prognostic factors are indicated in green. The strength of the relationship between two genes is symbolized by the thickness of the line. (C) Kaplan–Meier curves of the overall survival in the TCGA-LUAD cohort on the basis of the ferroptosis subtypes.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/743724/fcell-09-743724-HTML/image_m/fcell-09-743724-t001.jpg)
Table 1. Correlations between ferroptosis subtypes and clinicopathological parameters in TCGA-LUAD cohort.
Functional Enrichment and Immune Microenvironment Assessment of the Ferroptosis Subtypes
To further explore the biological pathway underlying each distinct ferroptosis subtype, GSVA was performed to determine pathway enrichment (Subtype A vs. Subtype B; Subtype B vs. Subtype C). As shown in Figure 3A, ferroptosis subtype A was significantly enriched in genes involved with RNA polymerase, pyrimidine metabolism, purine metabolism, cell cycle, and mismatch repair. Ferroptosis subtype B was markedly enriched in genes related to fatty acid and drug metabolism, including bile acid biosynthesis, alpha-linolenic acid metabolism, arachidonic acid metabolism, linoleic acid metabolism, and cytochrome P450 drug metabolism (Figure 3B). In addition, subtype B genes are also involved in cell adhesion and molecules (CAMS), intestinal immune network, and JAK-STAT signaling pathway (Figure 3A). Subtype C genes were predominantly related to pathogenic Escherichia coli infection, the p53 signaling pathway, ubiquitin-mediated proteolysis, and basal transcription factors (Figure 3B). Subsequent evaluation of infiltrated immune cells revealed conspicuous differences in the immune microenvironment among subtypes (Figure 4A). Particularly, as the subtype with the best prognosis, the ferroptosis subtype B genes were prominently enriched in Tfh, NKT, NK, CD4+ T, and CD8+ T cells (Figure 4A).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/743724/fcell-09-743724-HTML/image_m/fcell-09-743724-g003.jpg)
Figure 3. Biological characteristics of different ferroptosis subtypes. GSVA analysis revealed highly enriched KEGG pathways in distinct ferroptosis subtypes. The 20 most enriched pathways as identified in each comparison are visualized by heat map. (A) Ferroptosis subtype A vs. subtype B. (B) Ferroptosis subtype B vs. subtype C.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/743724/fcell-09-743724-HTML/image_m/fcell-09-743724-g004.jpg)
Figure 4. Feature selection and signature construction process. (A) Assessing the abundance of infiltrated immune cells in distinct ferroptosis subtypes. **p < 0.01, ***p < 0.001; ns, non-significant. (B) The number of unique and shared DEGs from different comparisons as revealed in a Venn diagram. (C) Principal component analysis of 174 genes selected with the Boruta feature selection algorithm. These genes can effectively discriminate the different ferroptosis subtypes. (D) The distribution of model occurrences after 1,000 iterations of LASSO regression. The six-gene prognostic signature was selected because it had a higher frequency than the other four models. (E) Kaplan–Meier curves and univariate Cox regression revealed that the FRP risk score was markedly related to overall survival in the TCGA-LUAD cohort. (F) Differences in FPI scores between the high-risk and low-risk groups. (G) Alluvial diagram showing the relations of ferroptosis subtype, FRP risk group, FPI, and status.
Construction of an Ferroptosis-Related Prognostic Signature for the TCGA-Lung Adenocarcinomas Cohort
To construct a prognostic model related to the defined ferroptosis subtypes, the intersecting DEGs identified in the pairwise comparisons were identified, and a total of 308 genes were selected as initial features (Figure 4B). After the dimension reduction by the Boruta algorithm, PCA was performed on the residual 174 genes, and the results indicated that the expression of these genes had great discriminatory power for the different subtypes (Figure 4C). Thus, these 174 genes constituted a representative collection of genes in the defined ferroptosis subtypes. These demonstrated that the dimensionality reduction process successfully removed unuseful features (n = 134) and retained relevant genes (n = 174). By applying univariate Cox regression, 56 genes with prognostic significance (p < 0.05) were retained for repeated LASSO Cox regression to generate the best model (Supplementary Table 2). A total of 1,000 iterations were conducted, and the distribution of model occurrences is shown in Figure 4D. As presented in Figure 4D, a prognostic signature with 6 genes (EIF5A, CACYBP, CYCS, ANLN, ARNTL2, and PPM1M) had the highest frequencies, with 881, which greatly exceeded those of the other four models and was, therefore, chosen as the FRP signature. The risk score of the FRP signature was calculated as follows: risk score = (0.0413 × EXPEIF5A) + (0.0068 × EXPCACYBP) + (0.0254 × EXPCYCS) + (0.1300 × EXPANLN) + (0.1503 × EXPARNTL2) – (0.0848 × EXPPPM1M). After calculating the risk score for each sample in the training cohort, samples were stratified into low-risk and high-risk groups based on the median risk score, and the survival analysis demonstrated that patients in the high-risk group had a significantly shorter overall survival time (log-rank p < 0.001, HR = 2.33, 95% CI = 1.72–3.17) than patients in the low-risk group (Figure 4E). The FPI was applied to quantify the activity of ferroptosis in each group, and the results suggested that the FPI in the high-risk group was significantly (p = 0.044) higher than that in the low-risk group (Figure 4F), which indicated a higher ferroptosis activity in samples that with higher risk. The association of attributes including ferroptosis subtype, risk group, FPI, and survival status of individual patients is presented in an alluvial diagram (Figure 4G).
Identification of the Ferroptosis-Related Prognostic Signature as an Independent Prognostic Factor
To decrease the effects of other confounding factors on the relationship between the FRP risk score and OS time, multivariate Cox regression analysis including age, sex, stage, and smoking history was conducted with the TCGA-LUAD cohort. The results proved that age, stage, and risk score (p < 0.001, HR = 2.17, 95% CI = 1.54–3.04) were independent prognostic factors of the OS time for LUAD patients (Figure 5).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/743724/fcell-09-743724-HTML/image_m/fcell-09-743724-g005.jpg)
Figure 5. Multivariate Cox regression analysis including age, sex, stage, smoking history, and FRP risk score in the TCGA-LUAD cohort.
Validation and Evaluation of the Predictive Performance of the Ferroptosis-Related Prognostic Signature in Five Validation Cohorts
The basic information related to the five validation cohorts is described in Table 2. Briefly, a total of 668 cases made up our validation datasets. Consistent and significant differences were observed in validation cohort 1 (log-rank p = 0.001, HR = 2.62, 95% CI = 1.44–4.76), validation cohort 2 (log-rank p = 0.002, HR = 2.62, 95% CI = 1.40–4.89), validation cohort 3 (log-rank p = 0.005, HR = 2.75, 95% CI = 1.32–5.74), validation cohort 4 (log-rank p = 0.014, HR = 2.39, 95% CI = 1.17–4.92), and validation cohort 5 (log-rank p < 0.001, HR = 2.59, 95% CI = 1.57–4.26; Figures 6A–E). In addition, the C-index was calculated to assess predictive performance and the C-index was 0.663 for the training cohort, 0.641 for validation cohort 1, 0.694 for validation cohort 2, 0.666 for validation cohort 3, 0.647 for validation cohort 4, and 0.685 for validation cohort 5 (Figure 6F). Statistical significance (p < 0.05) was observed between the actual value of the C-index and the null hypothesis (C-index = 0.5) in all datasets (Figure 6F).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/743724/fcell-09-743724-HTML/image_m/fcell-09-743724-g006.jpg)
Figure 6. Validation and evaluation of the performance of the FRP risk score. Kaplan–Meier curves and univariate Cox regression were performed with respect to the different risk groups in (A) validation cohort 1 (GSE13213), (B) validation cohort 2 (GSE30219), (C) validation cohort 3 (GSE31210), (D) validation cohort 4 (GSE3141), and (E) validation cohort 5 (GSE41271). (F) Assessing the predictive performance of the FRS by the C-index. The C-index was calculated for each cohort and compared to the null hypothesis. ***p < 0.001, ****p < 0.0001.
Depiction of the Immune Microenvironment and Implications of the Immunotherapeutic Response Between Risk Groups
The abundance of infiltrated immune cells was estimated using CIBERSORT and ImmuCellAI for each sample, and a heat map was generated to illustrate the association between immune cell infiltration and the FRP risk score (Figure 7A). The top five immune cell types with significant changes in abundance between risk groups were CD4+ T cells (p < 0.001), nTregs (p < 0.001), Tfhs (p < 0.001), Tems (p < 0.001), and monocytes (p < 0.001). The detailed results of these differential analyses are presented in Supplementary Table 3. Furthermore, an association analysis of the FRP risk score and potential immunotherapeutic biomarkers, including TMB, PD-L1 protein, and neoantigen, was performed. The results suggested that TMB in the high-risk group was significantly (p < 0.001) higher than that in the low-risk group (Figure 7B). A correlation analysis revealed that the FRP risk score was significantly associated (R = 0.365, p < 0.0001) with TMB (Figure 7C). In addition, the expression level of the PD-L1 protein (p < 0.001) and the numbers of subclonal neoantigens (p = 0.027), clonal neoantigens (p < 0.001), and total neoantigens (p < 0.001) were all significantly higher in the high-risk group than in the low-risk group (Figures 7D–G). These results implied that patients in the high-risk group may potentially benefit from immunotherapy.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/743724/fcell-09-743724-HTML/image_m/fcell-09-743724-g007.jpg)
Figure 7. Depiction of the immune microenvironment and implications of the immunotherapeutic response. (A) Heat map showing the immune microenvironment of the TCGA-LUAD cohort as estimated by CIBERSORT and ImmuCellAI. (B) Comparison of TMB between high-risk and low-risk groups. (C) Scatter plot presenting the correlation between TMB and risk score. Comparisons of (D) the expression of PD-L1 protein, (E) number of subclonal neoantigens, (F) number of clonal neoantigens, and (G) number of total neoantigens between different risk groups.
Extrapolating the Relationship Between the Ferroptosis-Related Prognostic Signature and Chemotherapy Response
Considering that chemotherapy constitutes an important component of LAUD treatment, drug susceptibility data including 138 compounds from the GDSC were applied to the training group to predict drug sensitivities for each sample in the TCGA-LUAD cohort. Based on the significance level, the top 10 compounds that might benefit patients in the high-risk group were identified: A.443654, obatoclax mesylate, epothilone B, GW843682X, BI.2536, TW.37, thapsigargin, ZM.447439, PF.562271, and S-Trityl-L-cysteine (Figures 8A–J). These results provide a novel perspective for exploring new drug candidates. A detailed list showing the significance value for all 138 compounds is available in Supplementary Table 4.
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/743724/fcell-09-743724-HTML/image_m/fcell-09-743724-g008.jpg)
Figure 8. Box plots based on the estimated IC50 of the 10 most significant compounds: (A) A.443654, (B) obatoclax mesylate, (C) epothilone B, (D) GW843682X, (E) BI.2536, (F) TW.37, (G) thapsigargin, (H) ZM.447439, (I) PF.562271, and (J) S-Trityl-L-cysteine. ****p < 0.0001.
Quantification of Gene Expression in Normal and Lung Adenocarcinoma Cell Lines
To further investigate the expression pattern in vitro, six genes that constituted the FRP signature were quantified by qPCR and those in normal bronchial epithelial cells (BEAS-2B cells) and LUAD cells (HCC827 cells) were compared. Of these genes, EIF5A and ARNTL2 were significantly overexpressed in tumor cells (p < 0.05; Figures 9A,E). In contrast, the expression levels of CYCS, ANLN, and PPM1M were significantly elevated in normal lung cells (p < 0.05; Figures 9C,D,F). No significant difference in CACYBP expression was observed between the two cell lines (Figure 9B).
![www.frontiersin.org](https://www.frontiersin.org/files/Articles/743724/fcell-09-743724-HTML/image_m/fcell-09-743724-g009.jpg)
Figure 9. Exploration of the expression of the six genes in normal and LUAD cell lines. Quantification and comparison of mRNA expression in the FRP signature as determined by qPCR: (A) EIF5A, (B) CACYBP, (C) CYCS, (D) ANLN, (E) ARNTL2, and (F) PPM1M. ***p < 0.001, ****p < 0.0001; ns, non-significant.
Discussion
Increasing studies about the association between ferroptosis and the TIME have been performed since the term ferroptosis was coined in 2012, and most of these studies demonstrated that ferroptosis might contribute to the immune surveillance and cytotoxic effects of CD8+ T cells (Wang et al., 2019), while others reported that ferroptosis facilitated the formation of a tumor-promoting immune microenvironment (Wen et al., 2019). As the roles of ferroptosis in shaping the tumor suppressive or promoting microenvironment remain elusive and most studies focus on a single ferroptotic regulator and a single cell type populating the TIME, the overall characteristics of the TIME linked to a specific ferroptotic signature in LUAD need to be elaborated.
Here, we identified three distinct ferroptosis subtypes based on the expression level of 37 ferroptosis driver and suppressor genes with prognostic significance (p < 0.05) in the TCGA-LUAD cohort. Patients with ferroptosis subtype B, in which most of the ferroptosis regulators were expressed at low levels, had the best prognosis. In addition, we evaluated the infiltrated immune cells in these three ferroptosis subtypes and found that subtype B was significantly enriched in Tfh, NKT, NK, CD4+ T, and CD8+ T cells, compared with the other two subtypes, suggesting that the TIME in association with ferroptosis subtype B was enriched in tumor-suppressive immune cells, which can partially explain the reason that patients with subtype B have the best prognosis.
Then, we constructed an FRP signature consisting of 6 genes (EIF5A, CACYBP, CYCS, ANLN, ARNTL2, and PPM1M) based on the defined ferroptosis subtypes and identified the FRP signature as an independent prognostic factor (p < 0.001, HR = 2.17, 95% CI = 1.54–3.04) in the training TCGA-LUAD cohort and five validation cohorts. Patients in the high-risk group had a significantly shorter overall survival time than patients in the low-risk group. In addition, we applied the FPI to quantify the overall activity of ferroptosis in each group and unsurprisingly found that the FPI in the low-risk group patients, who had a better prognosis, was significantly lower than that in the high-risk group, which was in accordance with the abovementioned results indicating that most of the ferroptosis regulators are expressed at low levels in patients with ferroptosis subtype B, who have the best prognosis, and a previous study showing that a high FPI predicts a poor prognosis (Liu et al., 2020).
Among the six genes in the FRP signature, ANLN and ARNTL2 were mechanistically explored in LUAD in previous studies. ANLN, a homolog of anillin whose stability and nuclear localization are regulated by the PI3K/AKT pathway, can bind to and activate RHOA, thus promoting cancer growth and development (Suzuki et al., 2005). In addition, Xu et al. (2019) reported that ANLN is involved in the metastasis of LUAD by promoting epithelial mesenchymal transformation of tumor cells, and Long et al. (2018) pointed out that higher expression of ANLN predicted a relatively poor prognosis. ARNTL2 is a transcription factor that is involved in promoting LUAD metastasis by orchestrating the expression of a complex prometastatic secretome (Brady et al., 2016). With respect to the other genes, EIF5A, a eukaryotic translation initiation factor, enhances the expression of a series of proliferation-associated proteins essential for cancer cells and is upregulated in several cancers, including pancreatic ductal adenocarcinoma, glioblastoma, leukemia, liver, colon, lung, cervical, and ovarian cancer (Mathews and Hershey, 2015; Strnadel et al., 2017). CACYBP, a calcyclin-binding protein, interacts with different partner proteins and, through these interactions, is involved in various cellular processes, such as proliferation, differentiation, cytoskeletal reorganization, and protein ubiquitination (Lian et al., 2019). For example, the Siah1-interacting protein (SIP)/CACYBP complex participates in the ubiquitination and proteasomal degradation of β-catenin, an oncogenic protein that regulates gene transcription (Jiang et al., 2019). Previous studies have reported that CACYBP contributes to the development of a wide range of human malignancies, including gastric, pancreatic, colon, breast, brain, and renal cancer (Topolska-Wos et al., 2016). CYCS, also known as cytochrome c, is a small soluble heme protein that shuttles electrons during the oxidative phosphorylation process and engages in the intrinsic apoptotic pathway under proapoptotic conditions (Ow et al., 2008). PPM1M, protein phosphatase 1M, is a manganese/magnesium ion-dependent serine/threonine phosphatase (Kamada et al., 2020). It was reported that PPM1M inhibits IL-1-induced activation of NF-κB by dephosphorylating IKK β (Henmi et al., 2009).
We also determined the abundance of infiltrated immune cell types in different risk groups using CIBERSORT and ImmuCellAI. The low-risk group with a better prognosis was enriched with CD4+ T cells and Tfhs, which was consistent with the infiltrating cell types associated with ferroptosis subtype B, and the abundance of nTregs, Tems, and monocytes was relatively low. Follicular helper T (Tfh) cells, emerging as a fifth helper T subset, interact with B cells, driving their differentiation into long-lived antibody-secreting plasma cells or memory B cells (King et al., 2008). In breast cancers, Tfh cells recruit immune cells to the TIME and promote tertiary lymphoid structure formation, in which effective antitumor immune responses can be generated and maintained (Gu-Trantien et al., 2013). Certain TIMEs were affected by several factors including exercise, age, diet, adiposity, the microbiome, sex, and tumor-dependent effects. Comprehensive consideration of the above factors during research will be more conducive to dissect the formation mechanism of TIME and lay the foundation for developing intervention measures (Binnewies et al., 2018). Our results suggested that patients in the low-risk group may have a better prognosis partially because of the infiltration of Tfh cells, which facilitate the formation of an antitumor TIME. Furthermore, we found that potential immunotherapeutic biomarkers, including TMB, the PD-L1 protein, and neoantigens, were elevated in the high-risk group, suggesting that patients in the high-risk group may benefit the most from immunotherapy. Our results revealed distinct immune and mutational profiles between different risk groups. Thus, we recommend to include ferroptosis-related genes or ferroptosis signature in future immunotherapy cohort for biomarker study as compared with previous ferroptosis associated prognostic model (Wang et al., 2021), which was directly constructed utilizing LASSO Cox regression model in TCGA-LUAD dataset and validated in one GEO dataset. By contrast, we firstly used ferroptosis regulators to identify three ferroptosis subtypes with prognostic significance and constructed a FRP signature by repeated LASSO Cox using 308 subtype-distinguishing genes. We also validated the signature in five GEO datasets and evaluated its implications of the immunotherapeutic response and chemotherapy, which was of great importance for further clinical application.
Inevitably, several limitations should be acknowledged in our study: (1) A total of 6 cohorts (a TCGA-LUAD cohort and 5 GEO cohorts) used for the construction and further validation of the FRP signature were downloaded from public databases retrospectively. A study with a prospective clinical trial cohort is needed to verify the clinical application value of our FRP signature. (2) Ferroptosis regulators obtained from the FerrDb website were manually extracted from previous literature reports, and therefore, some unselected crucial ferroptosis-associated genes may be missing in our study. In addition, ferroptosis regulators have been reported to regulate ferroptosis in human cells but not specifically in LUAD. (3) The FRP signature in our study was constructed based on DEGs identified in three ferroptosis subtypes, and the six genes in the FRP signature have not yet been identified to be involved in ferroptosis. Further studies are needed to determine the association and mechanisms between the six FRP genes and ferroptosis in LUAD. (4) Algorithmic calculation of immune cell abundance may not be as accurate as immunofluorescence, and further verification is required.
Conclusion
A robust ferroptosis-related prognostic signature was developed and validated with five independent datasets. Different tumor immune microenvironments were observed between the low-risk and high-risk groups. This signature may serve as a potential prognostic biomarker for LUAD in the future.
Data Availability Statement
Publicly available datasets were analyzed in this study. This data can be found here: the datasets collected in the current study are available in the TCGA (https://portal.gdc.cancer.gov/) and GEO repository (https://www.ncbi.nlm.nih.gov/geo/).
Author Contributions
YG and JH directed and designed the study. ZY, JW, RL, and HT extracted the data. SS and YY conducted the data analysis and experiments. SS wrote the manuscript. FT and QX reviewed and edited the manuscript. All authors read and approved the manuscript.
Funding
This study was supported by the National Key R&D Program of China (2018YFC1312100 and 2017YFC1311000), CAMS Initiative for Innovative Medicine (2017-I2M-1-005, 2017-I2M-2-003, and 2019-I2M-2-002), Beijing Municipal Science and Technology Commission (Z191100006619116), Innovation Team Development Project of Ministry of Education (IRT_17R10), and Peking Union Medical College Youth Fund (3332018144).
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
The authors thank all staff in the Department of Thoracic Surgery for their support during the study. The results shown here are in whole or part based upon data generated by the TCGA Research Network: https://www.cancer.gov/tcga.
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fcell.2021.743724/full#supplementary-material
Supplementary Table 1 | Quantitative RT-PCR primers used in this study.
Supplementary Table 2 | The list of 56 genes for model construction.
Supplementary Table 3 | The significance of differences in infiltrated immune cell between risk groups.
Supplementary Table 4 | List of significance for all 138 compounds.
Footnotes
- ^ https://portal.gdc.cancer.gov
- ^ http://www.zhounan.org/ferrdb
- ^ https://tcia.at/home
- ^ http://www.cbioportal.org/
References
Bild, A. H., Yao, G., Chang, J. T., Wang, Q., Potti, A., Chasse, D., et al. (2006). Oncogenic pathway signatures in human cancers as a guide to targeted therapies. Nature 439, 353–357. doi: 10.1038/nature04296
Binnewies, M., Roberts, E. W., Kersten, K., Chan, V., Fearon, D. F., Merad, M., et al. (2018). Understanding the tumor immune microenvironment (TIME) for effective therapy. Nat. Med. 24, 541–550.
Brady, J. J., Chuang, C. H., Greenside, P. G., Rogers, Z. N., Murray, C. W., Caswell, D. R., et al. (2016). An Arntl2-driven secretome enables lung adenocarcinoma metastatic self-sufficiency. Cancer Cell 29, 697–710. doi: 10.1016/j.ccell.2016.03.003
Charoentong, P., Finotello, F., Angelova, M., Mayer, C., Efremova, M., Rieder, D., et al. (2017). Pan-cancer immunogenomic analyses reveal genotype-immunophenotype relationships and predictors of response to checkpoint blockade. Cell Rep. 18, 248–262. doi: 10.1016/j.celrep.2016.12.019
Dai, E., Han, L., Liu, J., Xie, Y., Kroemer, G., Klionsky, D. J., et al. (2020a). Autophagy-dependent ferroptosis drives tumor-associated macrophage polarization via release and uptake of oncogenic KRAS protein. Autophagy 16, 2069–2083. doi: 10.1080/15548627.2020.1714209
Dai, E., Han, L., Liu, J., Xie, Y., Zeh, H. J., Kang, R., et al. (2020b). Ferroptotic damage promotes pancreatic tumorigenesis through a TMEM173/STING-dependent DNA sensor pathway. Nat. Commun. 11:6339.
Davis, S., and Meltzer, P. S. (2007). GEOquery: a bridge between the Gene Expression Omnibus (GEO) and BioConductor. Bioinformatics 23, 1846–1847. doi: 10.1093/bioinformatics/btm254
Dixon, S. J., Lemberg, K. M., Lamprecht, M. R., Skouta, R., Zaitsev, E. M., Gleason, C. E., et al. (2012). Ferroptosis: an iron-dependent form of nonapoptotic cell death. Cell 149, 1060–1072. doi: 10.1016/j.cell.2012.03.042
Ganesh, S., Shui, X., Craig, K. P., Park, J., Wang, W., Brown, B. D., et al. (2018). RNAi-mediated beta-catenin inhibition promotes T cell infiltration and antitumor activity in combination with immune checkpoint blockade. Mol. Ther. 26, 2567–2579. doi: 10.1016/j.ymthe.2018.09.005
Gao, J., Aksoy, B. A., Dogrusoz, U., Dresdner, G., Gross, B., Sumer, S. O., et al. (2013). Integrative analysis of complex cancer genomics and clinical profiles using the cBioPortal. Sci. Signal. 6:l1.
Gao, M., Yi, J., Zhu, J., Minikes, A. M., Monian, P., Thompson, C. B., et al. (2019). Role of mitochondria in ferroptosis. Mol. Cell 73, 354–363.e353.
Geeleher, P., Cox, N., and Huang, R. S. (2014). pRRophetic: an R package for prediction of clinical chemotherapeutic response from tumor gene expression levels. PLoS One 9:e107468. doi: 10.1371/journal.pone.0107468
Gu-Trantien, C., Loi, S., Garaud, S., Equeter, C., Libin, M., De Wind, A., et al. (2013). CD4(+) follicular helper T cell infiltration predicts breast cancer survival. J. Clin. Invest. 123, 2873–2892. doi: 10.1172/jci67428
Hänzelmann, S., Castelo, R., and Guinney, J. (2013). GSVA: gene set variation analysis for microarray and RNA-seq data. BMC Bioinformatics 14:7. doi: 10.1007/1-4020-7788-2_2
He, X., and Xu, C. (2020). Immune checkpoint signaling and cancer immunotherapy. Cell Res. 30, 660–669. doi: 10.1038/s41422-020-0343-4
Henmi, T., Amano, K., Nagaura, Y., Matsumoto, K., Echigo, S., Tamura, S., et al. (2009). A mechanism for the suppression of interleukin-1-induced nuclear factor kappaB activation by protein phosphatase 2Ceta-2. Biochem. J. 423, 71–78. doi: 10.1042/bj20090208
Jennis, M., Kung, C. P., Basu, S., Budina-Kolomets, A., Leu, J. I., Khaku, S., et al. (2016). An African-specific polymorphism in the TP53 gene impairs p53 tumor suppressor function in a mouse model. Genes. Dev. 30, 918–930. doi: 10.1101/gad.275891.115
Jiang, L., Kon, N., Li, T., Wang, S. J., Su, T., Hibshoosh, H., et al. (2015). Ferroptosis as a p53-mediated activity during tumour suppression. Nature 520, 57–62. doi: 10.1038/nature14344
Jiang, T. X., Zou, J. B., Zhu, Q. Q., Liu, C. H., Wang, G. F., Du, T. T., et al. (2019). SIP/CacyBP promotes autophagy by regulating levels of BRUCE/Apollon, which stimulates LC3-I degradation. Proc. Natl. Acad. Sci. U.S.A. 116, 13404–13413. doi: 10.1073/pnas.1901039116
Jiang, X., Stockwell, B. R., and Conrad, M. (2021). Ferroptosis: mechanisms, biology and role in disease. Nat. Rev. Mol. Cell Biol. 22, 266–282. doi: 10.1038/s41580-020-00324-8
Kamada, R., Kudoh, F., Ito, S., Tani, I., Janairo, J. I. B., Omichinski, J. G., et al. (2020). Metal-dependent Ser/Thr protein phosphatase PPM family: Evolution, structures, diseases and inhibitors. Pharmacol. Ther. 215:107622. doi: 10.1016/j.pharmthera.2020.107622
King, C., Tangye, S. G., and Mackay, C. R. (2008). T follicular helper (TFH) cells in normal and dysregulated immune responses. Annu Rev Immunol 26, 741–766. doi: 10.1146/annurev.immunol.26.021607.090344
Kursa, M. B., and Rudnicki, W. R. (2010). Feature selection with the Boruta package. J. Stat. Softw. 36, 1–13.
Leone, R. D., Zhao, L., Englert, J. M., Sun, I. M., Oh, M. H., Sun, I. H., et al. (2019). Glutamine blockade induces divergent metabolic programs to overcome tumor immune evasion. Science 366, 1013–1021. doi: 10.1126/science.aav2588
Lian, Y. F., Huang, Y. L., Zhang, Y. J., Chen, D. M., Wang, J. L., Wei, H., et al. (2019). CACYBP enhances cytoplasmic retention of P27(Kip1) to promote hepatocellular carcinoma progression in the absence of RNF41 mediated degradation. Theranostics 9, 8392–8408. doi: 10.7150/thno.36838
Liberzon, A., Subramanian, A., Pinchback, R., Thorvaldsdóttir, H., Tamayo, P., and Mesirov, J. P. (2011). Molecular signatures database (MSigDB) 3.0. Bioinformatics 27, 1739–1740. doi: 10.1093/bioinformatics/btr260
Liu, T., Jiang, L., Tavana, O., and Gu, W. (2019). The deubiquitylase OTUB1 mediates ferroptosis via stabilization of SLC7A11. Cancer Res. 79, 1913–1924. doi: 10.1158/0008-5472.can-18-3037
Liu, Z., Zhao, Q., Zuo, Z.-X., Yuan, S.-Q., Yu, K., Zhang, Q., et al. (2020). Systematic analysis of the aberrances and functional implications of ferroptosis in cancer. Iscience 23:101302. doi: 10.1016/j.isci.2020.101302
Long, X., Zhou, W., Wang, Y., and Liu, S. (2018). Prognostic significance of ANLN in lung adenocarcinoma. Oncol. Lett. 16, 1835–1840.
Mathews, M. B., and Hershey, J. W. (2015). The translation factor eIF5A and human cancer. Biochim. Biophys. Acta 1849, 836–844. doi: 10.1016/j.bbagrm.2015.05.002
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. 7:1902880. doi: 10.1002/advs.201902880
Newman, A. M., Liu, C. L., Green, M. R., Gentles, A. J., Feng, W., Xu, Y., et al. (2015). Robust enumeration of cell subsets from tissue expression profiles. Nat. Methods 12, 453–457. doi: 10.1038/nmeth.3337
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
Ow, Y. P., Green, D. R., Hao, Z., and Mak, T. W. (2008). Cytochrome c: functions beyond respiration. Nat. Rev. Mol. Cell Biol. 9, 532–542. doi: 10.1038/nrm2434
Ribas, A. (2015). Releasing the brakes on cancer immunotherapy. N. Engl. J. Med. 373, 1490–1492. doi: 10.1056/nejmp1510079
Ribas, A., and Wolchok, J. D. (2018). Cancer immunotherapy using checkpoint blockade. Science 359, 1350–1355.
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.
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
Schaer, D. A., Beckmann, R. P., Dempsey, J. A., Huber, L., Forest, A., Amaladas, N., et al. (2018). The CDK4/6 inhibitor abemaciclib induces a T Cell inflamed tumor microenvironment and enhances the efficacy of PD-L1 checkpoint blockade. Cell Rep. 22, 2978–2994. doi: 10.1016/j.celrep.2018.02.053
Smyth, G. K. (2005). “Limma: linear models for microarray data,” in Bioinformatics and Computational Biology Solutions Using R and Bioconductor, eds R. Gentleman, V. J. Carey, W. Huber, R. A. Irizarry, and S. Dudoit (New York, NY: Springer), 397–420. doi: 10.1007/0-387-29362-0_23
Stockwell, B. R., Friedmann Angeli, J. P., Bayir, H., Bush, A. I., Conrad, M., Dixon, S. J., et al. (2017). Ferroptosis: a regulated cell death nexus linking metabolism. Redox Biol. Dis. Cell 171, 273–285. doi: 10.1016/j.cell.2017.09.021
Strnadel, J., Choi, S., Fujimura, K., Wang, H., Zhang, W., Wyse, M., et al. (2017). eIF5A-PEAK1 signaling regulates YAP1/TAZ protein expression and pancreatic cancer cell growth. Cancer Res. 77, 1997–2007. doi: 10.1158/0008-5472.can-16-2594
Suzuki, C., Daigo, Y., Ishikawa, N., Kato, T., Hayama, S., Ito, T., et al. (2005). ANLN plays a critical role in human lung carcinogenesis through the activation of RHOA and by involvement in the phosphoinositide 3-kinase/AKT pathway. Cancer Res. 65, 11314–11325. doi: 10.1158/0008-5472.can-05-1507
Tibshirani, R. (1997). The lasso method for variable selection in the cox model. Stat. Med. 16, 385–395. doi: 10.1002/(sici)1097-0258(19970228)16:4<385::aid-sim380>3.0.co;2-3
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
Topolska-Wos, A. M., Chazin, W. J., and Filipek, A. (2016). CacyBP/SIP–Structure and variety of functions. Biochim. Biophys. Acta 1860, 79–85. doi: 10.1016/j.bbagen.2015.10.012
Wang, S. J., Li, D., Ou, Y., Jiang, L., Chen, Y., Zhao, Y., et al. (2016). Acetylation Is Crucial for p53-mediated ferroptosis and tumor suppression. Cell Rep. 17, 366–373.
Wang, W., Green, M., Choi, J. E., Gijon, M., Kennedy, P. D., Johnson, J. K., et al. (2019). CD8(+) T cells regulate tumour ferroptosis during cancer immunotherapy. Nature 569, 270–274. doi: 10.1038/s41586-019-1170-y
Wang, Y., Chen, W., Zhu, M., and Xian, L. (2021). Ferroptosis-related gene signature and patterns of immune infiltration predict the overall survival in patients with lung adenocarcinoma. Front. Mol. Biosci. 8:692530.
Wen, Q., Liu, J., Kang, R., Zhou, B., and Tang, D. (2019). The release and activity of HMGB1 in ferroptosis. Biochem. Biophys. Res. Commun. 510, 278–283. doi: 10.1016/j.bbrc.2019.01.090
Wilkerson, M. D., and Hayes, D. N. (2010). ConsensusClusterPlus: a class discovery tool with confidence assessments and item tracking. Bioinformatics 26, 1572–1573. doi: 10.1093/bioinformatics/btq170
Xu, J., Zheng, H., Yuan, S., Zhou, B., Zhao, W., Pan, Y., et al. (2019). Overexpression of ANLN in lung adenocarcinoma is associated with metastasis. Thorac. Cancer 10, 1702–1709. doi: 10.1111/1759-7714.13135
Zhang, Y., Shi, J., Liu, X., Feng, L., Gong, Z., Koppula, P., et al. (2018). BAP1 links metabolic regulation of ferroptosis to tumour suppression. Nat. Cell Biol. 20, 1181–1192. doi: 10.1038/s41556-018-0178-0
Zhang, Y., Zhuang, L., and Gan, B. (2019). BAP1 suppresses tumor development by inducing ferroptosis upon SLC7A11 repression. Mol. Cell Oncol. 6:1536845. doi: 10.1080/23723556.2018.1536845
Keywords: lung adenocarcinoma, ferroptosis, classification, immunotherapy, machine learning
Citation: Sun S, Yang Y, Yang Z, Wang J, Li R, Tian H, Tan F, Xue Q, Gao Y and He J (2021) Ferroptosis Characterization in Lung Adenocarcinomas Reveals Prognostic Signature With Immunotherapeutic Implication. Front. Cell Dev. Biol. 9:743724. doi: 10.3389/fcell.2021.743724
Received: 19 July 2021; Accepted: 22 September 2021;
Published: 20 October 2021.
Edited by:
Geng Chen, GeneCast Biotechnology Co., Ltd., ChinaReviewed by:
Zitong Gao, University of Hawai‘i at Mānoa, United StatesAhreum Kim, CHA University, South Korea
Copyright © 2021 Sun, Yang, Yang, Wang, Li, Tian, Tan, Xue, Gao and He. 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: Yibo Gao, gaoyibo@cicams.ac.cn; Jie He, prof.jiehe@gmail.com
†These authors have contributed equally to this work