Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Biosci., 07 March 2022
Sec. Molecular Diagnostics and Therapeutics
This article is part of the Research Topic Insights in Molecular Diagnostics and Therapeutics: 2021 View all 14 articles

Systematic Analysis and Validation of the Prognosis, Immunological Role and Biology Function of the Ferroptosis-Related lncRNA GSEC/miRNA-101-3p/CISD1 Axis in Lung Adenocarcinoma

Xiulin Jiang,&#x;Xiulin Jiang1,2Yixiao Yuan&#x;Yixiao Yuan1Lin TangLin Tang1Juan WangJuan Wang1Dahang ZhangDahang Zhang1Lincan Duan
Lincan Duan1*
  • 1The Department of Thoracic Surgery, The Third Affiliated Hospital of Kunming Medical University (Yunnan Tumor Hospital), Kunming, China
  • 2Key Laboratory of Animal Models and Human Disease Mechanisms of Chinese Academy of Sciences, Kunming Institute of Zoology, Kunming, China

Lung adenocarcinoma (LUAD) is the most common type of lung cancer, accounting for approximately 85% of pulmonary malignancies. Emerging evidence has demonstrated that ferroptosis plays a central role in both immunities as well as tumor proliferation. However, the clinical significance, immunological function, and upstream modulatory mechanism of ferroptosis-related genes in LUAD remain unclear. Here, we utilized various bioinformatics data to identify differentially expressed (DEGs) and prognosis-related ferroptosis (FRGs) genes in LUAD. Based upon identified DEGs, FRG, and ceRNA modulatory networks were constructed. Pearson’s correlation analysis was used to evaluate the correlation between FRGs and the tumor mutational burden, microsatellite instability, tumor-infiltrating immunity, cellular checkpoint control, and drug sensitivity in LUAD. A loss-of-function analysis was performed to verify the function of CISD1 in LUAD progression. Our findings revealed that certain FRGs (CISD1, ATP5MC3, PGD, SLC7A11, ACSL3, and FANCD2) are significantly upregulated in LUAD and that their elevated expression is associated with both advanced tumor stage and unfavorable prognosis. Furthermore, Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment results revealed these FRGs to be primarily involved in ferroptosis and glutathione metabolism in LUAD. We constructed a prognostic FRG-based model capable of accurately predicting LUAD patient overall survival with high specificity. The upstream lncRNA GSEC/miRNA-101-3p regulatory axis involving CISD1, ATP5MC3, and PGD was identified to be relevant in tumor progression. We also found GSEC, CISD1, ATP5MC3, and PGD to be upregulated, with miRNA-101-3p downregulated, in the setting of LUAD. Immunohistochemical analysis revealed CISD1, ATP5MC3, and PGD overexpression in LUAD tissue samples; CISD1 knockdown was noted to significantly inhibit LUAD proliferation and migration. In summary, this study characterizes relevant functional roles of the lncRNA GSEC/miR-101-3p axis in the setting of LUAD and suggests diagnostic and therapeutic biomarkers potentially useful in the clinical management of this illness.

Introduction

Lung cancer remains the commonest fatal condition globally, with an estimated 2.09 million new cases and 1.76 million deaths annually (Rivera and Wakelee, 2016). Non-small cell lung cancer, the most frequently encountered disease subtype, typically results in adverse clinical outcomes and carries a 5-year survival rate of only 18% (Schwartz and Cote, 2016). Despite advances in early detection and treatment standardization, management strategies remain varied and range from chemo-radiotherapy to immunotherapy. Urgent identification of novel, specific biomarkers with clinicopathological and prognostic significance is thus necessary for developing successful lung adenocarcinoma (LUAD) management strategies.

Ferroptosis describes a recently discovered type of programmed cell death distinct from apoptosis that is iron-dependent and characterized by both lipid peroxidation and reactive oxygen species production (Wang et al., 2019; Luo et al., 2021; Xu et al., 2021). Ferroptosis has been reported to play an essential role in lung cancer progression; the interaction of LINC00336 with ELAVL1 was found to result in the inhibition of ferroptosis in LUAD (Wang et al., 2019). Similarly, the interaction of G3BP1 with lncRNA-P53RRA was reported to facilitate both ferroptosis and apoptosis in lung cancer via regulation of p53 expression (Mao et al., 2018). Ferroptosis was additionally reported to be regulated via cellular redox and cell cycle signaling pathways (Yang and Stockwell, 2016). Prior studies have suggested several oncogenic signaling pathways related to the ferroptosis process to inhibit tumor progression via modulation of ferroptosis (Kang et al., 2019; Mou et al., 2019). For instance, P53, a well-studied tumor-suppressor gene, was reported to suppress cystine/glutamate antiporter expression and thus inhibit the ferroptosis process (Mou et al., 2019). However, a comprehensive analysis has yet to establish the prognostic value of ferroptosis-related genes (FRGs) and their relevant upstream regulatory axis in LUAD.

Analysis of gene expression, mutation, DNA methylation, and prognostic value in LUAD revealed CISD1, ATP5MC3, PGD, SLC7A11, ACSL3, and FANCD2 to be significantly upregulated in LUAD and elevated expression of these FRGs to be associated with advanced tumor stage and poor prognosis. Furthermore, KEGG enrichment findings revealed these FRGs to primarily be involved in both ferroptosis and glutathione metabolism in the setting of LUAD. Six prognosis-related genes were also used to construct an FRG-based prognostic model capable of accurately predicting LUAD patient overall survival with high specificity. Finally, CISD1, ATP5MC3, PGD, SLC7A11, ACSL3, and FANCD2 expression was confirmed to significantly correlate with tumor mutational burden (TMB), microsatellite instability (MSI), immune cell infiltration, cellular checkpoint dysfunction, and cancer sensitivity to drugs. Additionally, we identified the upstream regulatory (namely, the lncRNA GSEC/miRNA-101-3p) axis of relevant FRGs (CISD1, ATP5MC3, PGD) and found GSEC, CISD1, ATP5MC3, and PGD to be upregulated, with miRNA-101-3p downregulated, in the setting of LUAD. Knockdown of CISD1 was noted to significantly inhibit the proliferative and migratory capabilities of LUAD cells. In summary, by characterizing the functional roles of the lncRNA GSEC/miR-101-3p axis in LUAD, we suggest diagnostic and therapeutic biomarkers potentially useful in the future clinical management of this condition.

Materials and Methods

Data Collection

The Cancer Genome Atlas (TCGA)-LUAD cohort data and corresponding clinical information of 535 LUAD patients were downloaded from the TCGA website up to November 14, 2020 (https://portal.gdc.cancer.gov/repository). The gene expression profiles were normalized using the scale method provided in the “limma” R package. Data analysis was performed with the R (version 3.6.3) and ggplot2 (3.3.3) packages. The expression data were normalized to transcripts per kilobase million (TPM) values before further analysis.

Identification of Differentially Expressed 34 FRGs

A total of 34 FRGs were obtained from previous reviews (Su et al., 2020; Feng et al., 2020), which are shown in Supplementary Table S1. The difference in FRG expression in LUAD and normal tissues was identified using the “limma” packages. We then constructed a gene–gene interaction network for 34 FRGs using the GeneMANIA (http://www.genemania.org) database (Warde-Farley et al., 2010).

Gene Mutation Analysis of FRGs

The mutation frequency and oncoplot waterfall plot of 34 FRGs in LUAD patients were analyzed by the Gene Set Cancer Analysis (http://bioinfo.life.hust.edu.cn/web/GSCALite/) database (Liu et al., 2018).

Functional Enrichment Analysis

Gene Ontology (GO), consisting of the biological process (BP), cellular component (CC), and molecular function (MF) categories, was conducted with the “ggplot2” package in R software. Similarly, this package was also utilized to perform Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis (Ito and Murphy, 2013).

Development of the Ferroptosis-Related Gene Prognostic Model

Development of the FRG prognostic model was performed as documented (Lin et al., 2021). Briefly, we performed the Cox regression analysis to examine the prognostic significance of the FRGs. For Kaplan–Meier curves, p-values and hazard ratios (HRs) with 95% confidence intervals (CIs) were generated by log-rank tests and univariate Cox proportional hazard regression. FRGs with a significant prognostic value was selected for further analysis. Based on these prognostic FRGs, LASSO Cox regression analysis was then used to construct the prognostic model. The TCGA LUAD patients were divided into low- and high-risk subgroups according to the median risk score, and the overall survival time was compared between the two subgroups via Kaplan–Meier analysis. The predictive accuracy of each gene and the risk score was evaluated by performing time receiver-operating characteristic (ROC) analysis. Considering the clinical characteristics, a predicted nomogram was developed to predict the 1-, 3-, and 5-year overall survival. A forest was used to show the p-value, HR, and 95% CI of each variable through the “forest plot” R package.

Prediction of lncRNA and ceRNA Network Construction

We used starbase (http://starbase.sysu.edu.cn/) to predict the potential upstream miRNAs of FRGs and examine the expression, prognosis, and correlation between miRNA-101-3p and lncRNA, as well as to predict the binding sites among miRNA, mRNA, and lncRNA (Li et al., 2014). We utilized lncLocator (www.csbio.sjtu.edu.cn/bioinf/lncLocator.) and CPC2 (http://cpc2.cbi.pku.edu.cn) to explore the subcellular localization and the protein-coding ability of lncRNAs (Kang et al., 2017; Cao et al., 2018).

Analysis of the Immunological Roles of FRGs in Lung Cancer

We utilize TIMER (https://cistrome.shinyapps.io/timer/) and XCELL tools (https://xcell.ucsf.edu/) to examine the immunological roles of FRGs (Aran et al., 2017; Li et al., 2017), including the correlation between diverse immune cells and immune regulators. The TMB and MSI scores were obtained from TCGA. A correlation analysis between the FRG expression and TMB or MSI was performed using Spearman’s method.

Analysis of the Correlation Between the FRGs Expression and Drug Sensitivity

We utilize the Genomics of Drug Sensitivity in Cancer (GDSC) (www.cancerRxgene.org) and the Cancer Therapeutics Response Portal (CTRP) databases to analyze the correlation between FRG expression and drug sensitivity (Basu et al., 2013; Yang et al., 2013).

Cells and Cell Culture Conditions

The BEAS-2B cell line was purchased from the cell bank of Kunming Institute of Zoology and cultured in BEGM media (Lonza, Morrisville, NC, USA, CC-3170). Lung cancer cell lines, including A549 and H1975, were purchased from Cobioer (Nanjing, China), with the STR document; A549 and H1975 cells were all cultured in RPMI 1640 medium (Corning, Tewksbury, MA, USA) supplemented with 10% fetal bovine serum (Cat# 10099141C, Gibco, Grand Island, NY, USA) and 1% penicillin/streptomycin.

Quantitative Real-Time PCR

The qRT-PCR assay was performed as documented (Jiang et al., 2018). Total RNA was extracted according to the manufacturer’s protocol, and then reverse-transcribed using RT Reagent Kit (Takara Bio, Beijing, China, Cat# RR047A; Tiangen Biotech, Beijing, China, Cat# KR211-02). Real-time PCR was performed by FastStart Universal SYBR Green Master Mix (Roche, Cat# 04194194001; Tiangen Biotech, Beijing, China, Cat# FP411-02) using an Applied Biosystems 7500 machine. The primer sequences are listed as follows: GSEC-F: TTC​CAA​TTA​ACC​TGG​CCG​GAG, GSEC-R: GTC​AGC​CAA​CCC​ATT​GCA​AC, PGD-F: ATG​GCC​CAA​GCT​GAC​ATC​G, PGD-R: AAA​GCC​GTG​GTC​ATT​CAT​GTT, CISD1-F: CTG​ACT​TCC​AGT​TCC​AGC​GT, CISD1-R:TGATCAGAGGGCCCACATTG, ATP5G3-F: ATG​TTC​GCC​TGC​GCC​AAG, ATP5G3-R: GGC​AAA​CAA​GAT​CAA​GAA​CGC​A, miRNA-101-3p: TAC​AGT​ACT​GTG​ATA​ACT​GAA. β-Actin was used to normalize expression levels: β-actin-F: CTTCGCGGGCGACGAT, β-actin-R: CCA​TAG​GAA​TCC​TTC​TGA​CC. The expression quantification was obtained with the 2−ΔΔCt method.

Cell Proliferation and Cell Migration Assays

Cell proliferation assay was performed as previously described (Xu et al., 2020). Briefly, the indicated cells were plated onto 12-well plates, and the cell numbers were subsequently counted each day using an automated cell analyzer Countstar (Shanghai Ruiyu Biotech Co., Shanghai, China, IC1000). Cell migration assay was performed as previously described (Xiong et al., 2021). For transwell assay, 1-2×104 cells in 100 l serum-free medium were plated in an 8.0-µm, 24-well plate chamber insert (Corning Life Sciences, catalog no. 3422), with a medium containing 10% FBS at the bottom of the insert. Cells were incubated for 24 h and then fixed with 4% paraformaldehyde for 20 min. After washing, cells were stained with 0.5% crystal violet-blue. The positively stained cells were examined under the microscope.

Immunohistochemistry Staining

The immunohistochemistry staining assay was performed as documented (Dixon et al., 2012). Briefly, cancer tissues were collected; the primary antibody was incubated overnight and the second antibody incubated. Finally, the instrument was developed. Detailed information of the primers used in this study is as follows: PGD, Catalog number: #13389, dilution, 1:100, CST, Shanghai, China; CISD1, Catalog number: #83775, dilution, 1:200, CST, Shanghai, China; and ATP5MC3, Catalog number: ab129742, dilution, 1:100, Abcam, Shanghai, China.

Statistical Analysis

Analysis of the FRG expression lung cancer was estimated using t-tests. For survival analysis, the HR and p-value were calculated employing univariate Cox regression analysis. Kaplan–Meier analysis was used to examine the survival time of patients stratified according to high or low levels of the FRG expression. p-values less than 0.05 were considered statistically significant. For all figures, ∗, ∗∗, and ∗∗∗ indicate p < 0.05, p < 0.01, and p < 0.001, respectively.

Results

Analysis of FRG Expression and Mutation in LUAD

The expression of 34 FRGs and their function in the setting of LUAD was analyzed using data obtained from the TCGA database (Supplementary Table S1). Findings confirmed 26 genes to be upregulated and seven genes to be downregulated in lung cancer; no other significant differences were noted (Figures 1A–D). Furthermore, analysis of FRG mutations in LUAD revealed 97 of 127 (76.38%) samples to exhibit genetic mutations; ZEB1 exhibited the highest mutation frequency and was followed by ACACA and ABCC1 (Figure 1E). Finally, analysis of correlations among copy number variation (CNV), DNA methylation, and FRG expression in LUAD revealed FDFT1, IREB2, NFS1, TFRC, HSBP1, ACO1, LPCAT3, CARS, GSS, ACACA, ACSL3, CISD1, SLC1A5, and NCOA4 expression to positively correlate with their respective CNV (Supplementary Figure S1A). Meanwhile, DNA methylation was found to negatively associate with CBS, FDFT1, STEAP3, FADS2, SCL1A5, GCLM, GOT1, SLC7A11, CD44, ACSL3, and G6PD expression in LUAD (Supplementary Figure S1B). Collectively, these data indicate that CNV and DNA methylation exert a critical influence on FRG expression in LUAD.

FIGURE 1
www.frontiersin.org

FIGURE 1. Analysis of the expression and gene mutation of FRGs in LUAD. (A) The gene–gene interaction networks of FRG analysis by the GeneMANIA database. (B–D) The expression of FRGs in LUAD is examined by the TCGA database. (E) The mutation frequency and classification of FRGs in LUAD were examined by GSCA tools. For all figures, ∗, ∗∗, and ∗∗∗ indicate p < 0.05, p < 0.01, and p < 0.001, respectively.

Analysis of FRG Function in LUAD

To explore the function of FRGs in LUAD, GO and KEGG enrichment analyses were performed. Results revealed 34 FRGs to be primarily involved in cellular responses to oxidative and other chemical stress, reactive oxygen species metabolism, extracellular stimuli, nutrient levels, biosynthetic processes, and iron ion binding (Figure 2A). Furthermore, KEGG pathway enrichment analysis revealed these FRGs to be primarily involved in vascular fluid dynamics and pathogenesis of atherosclerosis, HIF-1 signaling, lipid metabolism, ferroptosis, amino acid metabolism, necroptosis, thyroid hormone synthesis, serotonergic communication, and inflammatory bowel disease pathogenesis, as well as the metabolism of arachidonic acid, aspartate, alanine, glutamate, and 2-monocarboxylic acid (Figure 2B). Collectively, these data confirm that FRGs influence cancer cell metabolism, with ferroptosis playing an essential role in LUAD progression.

FIGURE 2
www.frontiersin.org

FIGURE 2. Analysis of the functions of FRGs in LUAD. (A) The biological process involved by FRGs in LUAD. (B) The KEGG pathway is involved by FRGs in LUAD.

The Prognostic Value of FRGs in LUAD

An investigation of the correlation between FRG expression and tumor stage revealed only ATP5G3, CISD1, and PGD expression to positively correlate with the LUAD stage (Figure 3A). Evaluation of the prognostic value of FRGs revealed elevated SLC7A11, PGD, FANCD2, CISD1, ATP5G3, and ASCL3 expression to closely correlate with adverse clinical outcomes in LUAD (Figure 3B). The diagnostic value of FRGs was evaluated using ROC analysis. Findings confirmed SLC7A11, FANCD2, CISD1, and ATP5G3 expression to accurately predict the LUAD prognosis area under the curve (AUC > 0.8) (Figure 3C). Higher levels of SLC7A11, FANCD2, SIDS1, ATP5G3, and ASCL3 expression were similarly found to correlate with disease-specific survival in LUAD (Supplementary Figure S2A). Elevated SLC7A11, ATP5G3, and ASCL3 expressions were found to correlate with progression-free survival in LUAD (Supplementary Figure S2B). Our findings thus confirm that FRG plays an essential role in LUAD progression.

FIGURE 3
www.frontiersin.org

FIGURE 3. Analysis of the prognosis of FRGs in LUAD. (A) The pathologic stage of FRGs in LUAD is examined by the GEPIA database. (B) The overall survival of FRGs in LUAD was examined by the GEPIA database. (C) The ROC curve of FRGs in LUAD is examined by the TCGA database.

Construction of an FRG-Based Prognostic Model

A Lasso Cox regression prognostic model based on the six aforementioned FRGs was constructed (Figures 4A,B) with the risk score=(0.0387)*PGD+(0.082)*ATP5MC3+(0.2118)*CISD1+(0.0029)*SLC7A11+(0.2104)*ACSL3+(0.1528)*FANCD2. Patients suffering LUAD were divided into two groups based on the risk score. Risk score distribution, survival status, and gene expression data are presented in Figure 4C. Patient risk of death increased and survival time decreased as the risk score increased (Figure 4C). Kaplan–Meier survival analysis suggested that LUAD patients with high-risk scores possess a worse overall survival probability as compared to those with low-risk scores (median time = 3.3 vs. 4.2 years, p = 0.00198; Figure 4D), with AUCs of 0.653, 0.617, and 0.588 noted on 1-, 3-, and 5-year receiver operating characteristic (ROC) curves, respectively (Figure 4E). The above results were verified by the GEO dataset; Kaplan–Meier survival analysis suggested that LUAD patients with high-risk scores possess a worse overall survival probability as compared to those with low-risk scores (median time = 2.9 vs. 4.9 years, p = 6.52e-08), with AUCs of 0.66, 0.693, and 0.664 noted on 1-, 3-, and 5-year ROC curves, respectively (Supplementary Figures S6A–D). A nomogram predictive of patient survival state was constructed. Univariate and multivariate analyses revealed ATP5MC3 expression and TNM stage to be independent factors affecting LUAD prognosis (Figures 5A,B). The nomogram could accurately predict 1-, 3-, and 5‐year overall survival rates as compared with an ideal model of the entire cohort (Figures 5C,D).

FIGURE 4
www.frontiersin.org

FIGURE 4. Construction of a prognostic FRG model in LUAD. (A) LASSO coefficient profiles of six FRGs. (B) Plots of the ten-fold cross-validation error rates. (C) Distribution of risk score, survival status, and the expression of six prognostics FRGs in LUAD. (D,E) Overall survival curves for LUAD patients in the high-/low-risk group and the ROC curve of measuring the predictive value. *p < 0.05, **p < 0.01, ***p < 0.001.

FIGURE 5
www.frontiersin.org

FIGURE 5. Construction of a predictive nomogram in LUAD. (A,B) Hazard ratio and p‐value of the constituents involved in univariate and multivariate Cox regression considering clinical parameters and six prognostics FRGs in LUAD. (C,D) Nomogram predicts the 1-, 3-, and 5-year overall survival rate of LUAD patients and calibration curve for the overall survival nomogram model in the discovery group. A dashed diagonal line represents the ideal nomogram. *p < 0.05, **p < 0.01, ***p < 0.001.

Analysis of the Correlation Between FRG Expression and TMB, MSI, and Drug Sensitivity

Emerging evidence has suggested that TMB and MSI may serve as potential predictive biomarkers for immunotherapy efficacy in the setting of lung cancer (Goodman et al., 2019). As such, several studies have reported that the expression of FRGs significantly correlates with tumor immune infiltration. Correlation analysis revealed FANCD2, ATP5MC3, and PGD expression to positively correlate with TMB (Supplementary Figure S3A); FANCD2 expression was also noted to positively correlate with MSI. Meanwhile, CISD1 expression was found to negatively correlate with MSI in LUAD (Supplementary Figure S3B). To explore potential therapeutic targets, gene-set co-expression analysis was utilized to explore the relationship between FRG expression and drug sensitivity. Results revealed SLC7A11 expression to positively correlate with sensitivity to PRIMA-1, PX-12, necrosulfonamide, methylstat, piperlongumine, SMER-3, NSC95397, manumycin A, ML162, PL-DI1S, 3R-RSL-3, pifithrin-mu, and cerulenin (r > 0.36). Expression of ACSL3 was found to positively correlate with sensitivity to tozasertib, PRIMA-1, PX-12, manumycin A, BRD-K30748066, ML210, BRD-A94377914, 1S, 3R-RSL-3, and docetaxel (r > 0.30). On the contrary, PGD expression was found to negatively correlate with sensitivity to tivantinib, SCH-79797, BI-2536, GW-405833, GSK461364, nakiterpiosin, docetaxel, SB-743921, bafilomycin A1, linifanib, ceranib-2, BRD-K70511574, BRD-K01737880, FQI-2, and BRD-K30748066 (r < -0.22). Expression of ATP5MC3 was found to negatively correlate with sensitivity to BRD-K30748066, marinopyrrole A, COL-3, dinaciclib, BI-2536, alvocidib, methotrexate, and bafilomycin A1 (r < -0.22). The expression of FANCD2 was found to negatively correlate with sensitivity to GSK-J4, BRD-K30748066, COL-3, docetaxel, GSK461364, BRD-K66453893, tivantinib, BI-2536, narciclasine, BRD-K70511574, and triazolothiadiazine (r < -0.34; Supplementary Figure S3D). These findings thus underscore that FRG expression either positively or negatively correlates with LUAD sensitivity to drugs.

FRG Expression Associated With LUAD Immune Infiltration

Analysis of CISD1, FANCD2, PGD, ASCL1, ATP5MC3, and SLC7A11 expression in the setting of C1 and C2 LUAD subtypes revealed high CISD1, FANCD2, and ATP5MC3 expression in the C2 subtype and high PGD and SLC7A11 expression in the C1 subtype (Supplementary Figures S4A,B). Analysis of data obtained from the TIMER database revealed FRG somatic copy number alterations to significantly correlate with levels of immune cell infiltration in LUAD (Supplementary Figures S5A–F). As ferroptosis plays crucial roles in both immunity and pulmonary carcinogenesis (Huang et al., 2021), the correlation between CISD1, FANCD2, PGD, ASCL3, ATP5MC3, and SLC7A11 expression and tumor immune infiltration in LUAD was evaluated utilizing data obtained from the TIMER database. Expression of ASCL3 was found to positively correlate with the level of infiltration by central memory T cells, T helper cells, T follicular helper cells, B cells, natural killer cells, CD56 bright natural killer cells, activated dendritic cells, cytotoxic cells, and mast cells. On the contrary, ATP5MC3 expression was found to negatively correlate with the level of infiltration by cytotoxic cells, dendritic cells, interstitial dendritic cells, T cells, eosinophils, macrophages, T helper one cells, B cells, mast cells, plasmacytoid dendritic cells, T follicular helper cells, natural killer cells, and central and effector memory T cells. The expression of CISD1 was found to negatively correlate with the level of infiltration by mast cells, T cells, cytotoxic cells, B cells, eosinophils, natural killer cells, T helper 17 cells, CD56 bright natural killer cells, CD8 T cells, plasmacytoid dendritic cells, T follicular helper cells, and central and effector memory T cells (Figure 6A). The expression of FANCD2 was found to negatively correlate with the level of infiltration by cytotoxic cells, macrophages, CD56 bright natural killer cells, neutrophils, B cells, T helper 17 cells, natural killer cells, plasmacytoid dendritic cells, T follicular helper cells, dendritic cells, CD8 T cells, eosinophils, interstitial dendritic cells, and mast cells. Similarly, the PGD expression was found to negatively correlate with the level of infiltration by mast cells, eosinophils, interstitial dendritic cells, natural killer cells, cytotoxic cells, macrophages, dendritic cells, plasmacytoid dendritic cells, CD8 T cells, B cells, T helper one cell, T cells, T follicular helper cells, and effector and memory T cells. The expression of SLC7A11 was found to negatively correlate with the level of infiltration by eosinophils, mast cells, cytotoxic cells, activated dendritic cells, B cells, CD8 T cells, natural killer cells, regulatory T cells, T cells, plasmacytoid dendritic cells, T helper one cell, macrophages, interstitial dendritic cells, dendritic cells, and T follicular helper cells (Figure 6B).

FIGURE 6
www.frontiersin.org

FIGURE 6. Analysis of the association between FRG expression and immune infiltration level in LUAD. (A,B) The association between FRG expression and immune infiltration level in LUAD.

As immune checkpoints play a crucial role in tumor immunosuppression, we analyzed the distribution of a relevant gene expression in stage I–IV LUAD tissues (Figure 7A). We additionally evaluated the relationship between FRG and immune checkpoint-related gene (including CD274, CTLA4, HAVCR2, LAG3, PDCD1, PDCD1LG12, TIGIT, SIGLEC15) expression in LUAD via Pearson correlation analysis; results revealed the expression of SLC7A11, PGD, CISD1, ATP5MC3, and ACSL3 to negatively correlate with that of CD274, CTLA4, HAVCR2, LAG3, PDCD1, PDCD1LG12, TIGIT, and SIGLEC15. Meanwhile, FANCD2 expression was found to positively correlate with checkpoint-related gene expression (Figure 7B). These findings confirm that FRG expression significantly correlates with that of immune checkpoint-related genes in LUAD.

FIGURE 7
www.frontiersin.org

FIGURE 7. Analysis of the correlation between FRG expression and diverse immune modulators. (A) The expression of immune checkpoint-related genes in diverse pathologic stages and normal samples. (B) The correlation between FRGs expression and diverse immune checkpoint-related genes in LUAD. For all figures, ∗, ∗∗, and ∗∗∗ indicate p < 0.05, p < 0.01, and p < 0.001, respectively.

Analysis of the Upstream FRG Molecular Regulatory Axis

The aforementioned evidence suggests that ATP5MC3, PGD, and CISD1 expression correlates with both stage and progression of LUAD. To further explore the relevant upstream FRG regulatory axis, data obtained from various public databases were used to construct a network of mRNA–miRNA–lncRNA interactions. We found miRNA-101-3p to most significantly bind the 3′UTR of ATP5MC3, PGD, and CISD1 (Figure 8A). Further analysis revealed miRNA-101-3p to be decreased in lung cancer based on TCGA-LUAD and GSE74190 datasets (Figures 8B,C). Moreover, downregulated miRNA-101-3p expression correlated with the manifestation of poorer clinical features and prognosis among lung cancer patients (Figures 8D,E). Analysis of miRNA-101-3p ROC curve data revealed an AUC value of 0.864 among lung cancer patients (Figure 8F). We also found miRNA-101-3p expression to significantly negatively correlate with ATP5MC3, PGD, and CISD1 expression in LUAD (Figure 8G). To predict and obtain the PVT1 and GSEC lncRNAs, starbase and lncbase were utilized (Figures 8H–J). According to the ceRNA theory, lncRNA negatively and positively correlates with miRNA and mRNA expression, respectively. We found only GSEC to negatively correlate with miRNA-101-3p expression (Figure 8K) and positively correlate with ATP5MC3, PGD, and CISD1 expression in LUAD (Figure 9A). Importantly, GSEC was found to be highly expressed in the setting of this malignancy based on the GSE81089 dataset (Figure 9B) and primarily localized to the cytoplasm as determined using the lncLocator and AnnoLnc RNA tools (Figure 9C). However, GSEC was not found to possess coding potential (Figure 9D). High levels of GSEC expression were found to correlate not only with tumor stage but also with poor prognosis in LUAD. Analysis of the GSEC ROC curve data revealed an AUC value of 0.813 in LUAD patients (Figures 9E–G). In summary, our findings confirm that the lncRNA GSEC/miRNA-101-3p axis modulates ATP5MC3, PGD, and CISD1 expression in LUAD.

FIGURE 8
www.frontiersin.org

FIGURE 8. Prediction and analysis of the upstream miRNAs of FRGs in LUAD. (A) Prediction of the potential miRNAs of FRGs in LUAD examined by starbase. (B,C) The expression of miRNA-101-3p in LUAD was examined by TCGA/GEO databases. (D) The correlation between FRGs and clinical features in LUAD. (E) The prognosis value of miRNA-101-3p in LUAD. (F) ROC curve analyses and AUC values for miRNA-101-3p in LUAD. (G) Pearson’s correlation analysis determined the correlation between miRNA-101-3p and ATP5MC3, CISD1, and PGD expression in LUAD examined by starbase. (H) Explored the potential lncRNAs of miRNA-101-3p in LUAD by starbase. (I,J) The expression of GESC and PVT1 in LUAD was examined by starbase. (K) Pearson’s correlation analysis determined the correlation between miRNA-101-3p and GSEC in LUAD. *p < 0.05, **p < 0.01, ***p < 0.001.

FIGURE 9
www.frontiersin.org

FIGURE 9. Construction of lncRNA/miRNA/FRG interaction network. (A) Pearson’s correlation analysis determined the correlation between GSEC and ATP5MC3, CISD1, and PGD in LUAD. (B) The expression of GESC in LUAD is examined by the GEO dataset. (C) The subcellular location of GESC is examined by lncLocator and Annolnc2 databases. (D) The coding ability of GESC is examined by coding potential calculator databases. (E) The pathologic stage of GESC in LUAD is examined by GEPIA databases. (F) ROC curve analyses and AUC values for GSEC in LUAD. (G) The prognosis of GESC in LUAD is examined by GEPIA databases. *p < 0.05, **p < 0.01, ***p < 0.001.

Knockdown of CISD1 Inhibits LUAD Proliferation and Migration

To verify the aforementioned findings, a qRT-PCR assay was used to detect GSEC, ATP5MC3, PGD, CISD1, and miRNA-101-3p expression in LUAD cell lines. While we found GSEC, ATP5MC3, PGD, and CISD1 expression to be upregulated in LUAD cells (Figure 10A), miRNA-101-3p expression was noted to be decreased (Figure 10B). The immunohistochemical evaluation further revealed PGD, ATP5MC3, and CISD1 overexpression in LUAD tissues (Figures 10C–E). As prior literature has detailed the roles of PGD and ATP5MC3 in LUAD, we focused on further studying the influence of CISD1 on LUAD cell proliferation and migration via the utilization of target siRNA transfected into H1975 cells. We found that CISD1 expression was decreased after CISD1 knockdown in this cell line (Figure 10F). Growth curve and colony formation findings revealed that CISD1 knockdown significantly inhibits cell proliferation (Figures 10G,H). Transwell assays indicated that the migratory capability of H1975 cells was significantly suppressed via CISD1 downregulation (Figure 10I). These findings confirm that CISD1 functions as an oncogene, promoting cell growth and migration among LUAD cells.

FIGURE 10
www.frontiersin.org

FIGURE 10. Depletion of CISD1 inhibits LUAD cell proliferation and migration. (A,B) The expression of GSEC, FRGs, and miRNA-101-3p in LUAD cell lines was examined by qRT-PCR assay. (C–E) The expression of PGD, ATP5MC3, and CISD1 in LUAD tissue was examined by IHC assay. (F) The expression of CISD1 in H1975 cells after knockdown of CISD1 was examined by qRT-PCR assay. (G,H) The growth curve and colony formation assays were utilized to detect the silencing of CISD1 on the growth of H1975 cells. (I) The transwell assay was utilized to detect the silencing of CISD1 on the migration of H1975 cells. *p < 0.05, **p < 0.01, ***p < 0.001.

Discussion

An iron-dependent type of programmed cell death, ferroptosis differs from other forms of cell death such as apoptosis, necrosis, and autophagy. The ferroptosis process generally results in excessive lipid peroxidation and is induced by abnormalities in cellular redox processes. Importantly, ferroptosis is understood to play an indispensable role in the elimination of carcinogenic cells (Hirschhorn and Stockwell, 2019). Activation of Ras/mitogen-activated protein kinase signaling increases the sensitivity of cancer cells to ferroptosis and thus leads to the maintenance of relative iron abundance in cancer via control of transferrin receptor and ferritin expression (Dixon et al., 2012). Knockdown of FMS-like tyrosine kinase three was reported to reduce lipid peroxidation and inhibit p22phox activity, resulting in the inhibition of ferroptosis among malignant cells (Song et al., 2018). Recently, 5-monophosphate was reported to increase Beclin1 phosphorylation and inhibit system Xc activity, thus facilitating ferroptosis (Kang et al., 2018). However, the modulatory mechanism of the ferroptosis process and associated FRG expression correlated with LUAD progression remains to be elucidated.

Here, we first investigated the expression and prognostic value of FRGs in LUAD. We found a total of 26 genes to be upregulated; seven genes were found to be downregulated in LUAD as compared with normal tissue. Meanwhile, the expression of FDFT1, IREB2, NFS1, TFRC, HSBP1, ACO1, LPCAT3, CARS, GSS, ACACA, ACSL3, CISD1, SLC1A5, and NCOA4 was found to positively correlate with respective CNV in LUAD. Levels of DNA methylation were found to negatively associate with CBS, FDFT1, STEAP3, FADS2, SCL1A5, GCLM, GOT1, SLC7A11, CD44, ACSL3, and G6PD expression in LUAD. Collectively, these data indicate that CNV and DNA methylation exert critical effects on FRG expression in LUAD.

Functional enrichment analysis revealed 33 FRGs to be primarily involved in vascular fluid dynamics and atherosclerosis, HIF-1a signaling, lipid metabolism, ferroptosis, cysteine and methionine metabolism, glutathione metabolism, amino acid biosynthesis, necroptosis, arachidonic acid metabolism, thyroid hormone synthesis, serotonergic communication, alanine metabolism, aspartate metabolism, glutamate metabolism, inflammatory bowel disease, 2-monocarboxylic acid metabolism, and arginine biosynthesis. Our findings demonstrate that FRGs play an important role in both LUAD pathogenesis and illness progression.

Prognosis analysis revealed increased SLC7A11, PGD, FANCD2, CISD1, ATP5G3, and ASCL3 expression to closely correlate with poor LUAD clinical prognosis. Importantly, levels of ATP5G3, CISD1, and PGD expression were found to positively correlate with tumor stage. Analysis ROC curve data revealed SLC7A11, FANCD2, CISD1, and ATP5G3 to be highly accurate markers in LUAD (AUC > 0.8). These findings suggest that FRG expression correlates with LUAD prognosis. A highly accurate five-gene prognostic model capable of predicting overall LUAD patient survival was constructed based on SLC7A11, PGD, FANCD2, SIDS1, ATP5G3, and ASCL3 using Lasso Cox regression analysis. A predictive nomogram useful in predicting 3- and 5‐year overall survival rates as compared to an ideal model in the studied cohort was similarly developed. Here, we identified a pyroptosis-related prognostic gene signature for LUAD, further advancing potential prognostic methods relevant to LUAD.

Here, SLC7A11, PGD, FANCD2, CIDS1, ATP5G3, and ASCL3 were found to serve as effective gene signatures for predicting LUAD prognosis. An important cystine/glutamate antiporter regulated by p53, SLC7A11 plays suppressive roles in the ferroptosis process. The SLC7A11-mediated importation of cystine enhances glutathione biosynthesis and increases GPX4-mediated detoxification of lipid peroxides, thereby inhibiting ferroptosis (Stockwell and Jiang, 2020). A prior study reported PGD capable of predicting overall survival in the setting of papillary thyroid carcinoma (Yang et al., 2021). In addition, FANCD2 was reported to inhibit ferroptosis and protect from bone marrow injury (Song et al., 2016).

Importantly, we found FRG expression to significantly correlate with tumor immune infiltration. Furthermore, we found that SLC7A11, PGD, CISD1, ATP5MC3, and ACSL3 expression significantly negatively correlates with CD274, CTLA4, HAVCR2, LAG3, PDCD1, PDCD1LG12, TIGIT, and SIGLEC15 expression. However, FANCD2 expression was found to positively correlate with checkpoint-related gene expression in LUAD. As such, these data confirm that FRG expression is associated with immune checkpoint-related gene expression in LUAD.

The lncRNA GSEC was previously reported to upregulate EIF5A2 expression via miR-588 sponging and subsequently facilitate cell proliferation, migration, and invasion capabilities in osteosarcoma (Liu et al., 2020). High levels of GSEC expression were reported in triple-negative breast cancer tissue and cell lines, with GSEC knockdown significantly decreasing these proliferative, migratory, and invasive capabilities. The sponging of miR-202-5p by GSEC, along with resultant increased AXL expression, was reported to enhance LUAD progression (Zhang et al., 2021). Similarly, miRNA-101-3p was reported to inhibit both proliferative and migratory capabilities of gastric carcinoma cells and enhance apoptosis via PIM-1 expression modulation (Wu et al., 2019). Indeed, in an LUAD setting of miRNA-101-3p downregulation, overexpression of miR-101-3p mimics was reported to reduce rates of tumor growth and progression in vitro (Liu et al., 2021). In this study, data obtained from various public databases were used to illustrate novel mRNA–miRNA–lncRNA interactions (namely, the lncRNA GSEC/miR-101-3p/ATP5MC3/PGD/CISD1 axis) detailing modulation of FRG expression following the ceRNA hypothesis (Figure 11). Although we found GSEC, ATP5MC3, PGD, and CISD1 to be upregulated in LUAD cell lines, miRNA-101-3p expression was found to be decreased. The immunohistochemical evaluation also confirmed increased expression of ATP5MC3, PGD, and CISD1 in LUAD tissue. Knockdown of CISD1 was found to significantly inhibit the proliferative and migratory capabilities of LUAD cells.

FIGURE 11
www.frontiersin.org

FIGURE 11. A working model for lncRNA GSEC/miR-101-3p/CISD1 axis in LUAD.

Conclusion

Here, we found that GSEC, CISD1, ATP5MC3, and PGD are upregulated, and miRNA-101-3p is downregulated, in LUAD. Knockdown of CISD1 was found to significantly inhibit the proliferative and migratory capabilities of LUAD cells. In summary, this study characterizes the functional roles of the lncRNA GSEC/miR-101-3p/CISD1 axis in LUAD and suggests potential diagnostic and therapeutic biomarkers for future clinical application in the management of this illness.

Data Availability Statement

The original contributions presented in the study are included in the article/Supplementary Material; further inquiries can be directed to the corresponding author.

Author Contributions

XJ and YY designed this work, performed the related assay, and analyzed the data. LT, JW, and DZ contributed to the study materials. LD supervised and wrote the manuscript. All authors have read and approved the final version of the manuscript.

Funding

This work was supported by the National Nature Science Foundation of China (82160508) and Yunnan Applied Basic Research Projects (YNWRMY-2019-067, 2019FE001) and Yunnan Province Specialized Training Grant for High-Level Healthcare Professionals (D-201614).

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 would like to thank support from the Department of Thoracic Surgery, The Third Affiliated Hospital of Kunming Medical University (Yunnan Tumor Hospital), Kunming, China.

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2021.793732/full#supplementary-material

References

Aran, D., Hu, Z., and Butte, A. J. (2017). xCell: Digitally Portraying the Tissue Cellular Heterogeneity Landscape. Genome Biol. 18 (1), 220. doi:10.1186/s13059-017-1349-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Basu, A., Bodycombe, N. E., Cheah, J. H., Price, E. V., Liu, K., Schaefer, G. I., et al. (2013). An Interactive Resource to Identify Cancer Genetic and Lineage Dependencies Targeted by Small Molecules. Cell 154 (5), 1151–1161. doi:10.1016/j.cell.2013.08.003

PubMed Abstract | CrossRef Full Text | Google Scholar

Cao, Z., Pan, X., Yang, Y., Huang, Y., and Shen, H.-B. (2018). The lncLocator: a Subcellular Localization Predictor for Long Non-coding RNAs Based on a Stacked Ensemble Classifier. Bioinformatics 34 (13), 2185–2194. doi:10.1093/bioinformatics/bty085

PubMed Abstract | CrossRef Full Text | Google Scholar

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 (5), 1060–1072. doi:10.1016/j.cell.2012.03.042

PubMed Abstract | CrossRef Full Text | Google Scholar

Feng, H., Schorpp, K., Jin, J., Yozwiak, C. E., Hoffstrom, B. G., Decker, A. M., et al. (2020). Transferrin Receptor Is a Specific Ferroptosis Marker. Cel Rep. 30 (10), 3411–3423. doi:10.1016/j.celrep.2020.02.049

PubMed Abstract | CrossRef Full Text | Google Scholar

Goodman, A. M., Sokol, E. S., Frampton, G. M., Lippman, S. M., and Kurzrock, R. (2019). Microsatellite-Stable Tumors with High Mutational Burden Benefit from Immunotherapy. Cancer Immunol. Res. 7 (10), 1570–1573. doi:10.1158/2326-6066.cir-19-0149

PubMed Abstract | CrossRef Full Text | Google Scholar

Hirschhorn, T., and Stockwell, B. R. (2019). The Development of the Concept of Ferroptosis. Free Radic. Biol. Med. 133, 130–143. doi:10.1016/j.freeradbiomed.2018.09.043

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, J., Liu, D., Wang, Y., Liu, L., Li, J., Yuan, J., et al. (2021). Ginseng Polysaccharides Alter the Gut Microbiota and Kynurenine/tryptophan Ratio, Potentiating the Antitumour Effect of Antiprogrammed Cell Death 1/programmed Cell Death Ligand 1 (Anti-PD-1/pd-L1) Immunotherapy. Gut 2021 [Epub ahead of print], 321031. doi:10.1136/gutjnl-2020-321031

CrossRef Full Text | Google Scholar

Ito, K., and Murphy, D. (2013). Application of Ggplot2 to Pharmacometric Graphics. CPT Pharmacometrics Syst. Pharmacol. 2 (10), e79. doi:10.1038/psp.2013.56

PubMed Abstract | CrossRef Full Text | Google Scholar

Jiang, L.-P., Fan, S.-Q., Xiong, Q.-X., Zhou, Y.-C., Yang, Z.-Z., Li, G.-F., et al. (2018). GRK5 Functions as an Oncogenic Factor in Non-small-cell Lung Cancer. Cell Death Dis 9 (3), 295. doi:10.1038/s41419-018-0299-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, R., Kroemer, G., and Tang, D. (2019). The Tumor Suppressor Protein P53 and the Ferroptosis Network. Free Radic. Biol. Med. 133, 162–168. doi:10.1016/j.freeradbiomed.2018.05.074

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, R., Zhu, S., Zeh, H. J., Klionsky, D. J., and Tang, D. (2018). BECN1 Is a New Driver of Ferroptosis. Autophagy 14 (12), 2173–2175. doi:10.1080/15548627.2018.1513758

PubMed Abstract | CrossRef Full Text | Google Scholar

Kang, Y.-J., Yang, D.-C., Kong, L., Hou, M., Meng, Y.-Q., Wei, L., et al. (2017). CPC2: a Fast and Accurate Coding Potential Calculator Based on Sequence Intrinsic Features. Nucleic Acids Res. 45 (W1), W12–w16. doi:10.1093/nar/gkx428

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, J. H., Liu, S., Zhou, H., Qu, L. H., and Yang, J. H. (2014). starBase v2.0: Decoding miRNA-ceRNA, miRNA-ncRNA and Protein-RNA Interaction Networks from Large-Scale CLIP-Seq Data. Nucleic Acids Res. 42, D92–D97. doi:10.1093/nar/gkt1248

PubMed Abstract | CrossRef Full Text | Google Scholar

Li, T., Fan, J., Wang, B., Traugh, N., Chen, Q., Liu, J. S., et al. (2017). TIMER: A Web Server for Comprehensive Analysis of Tumor-Infiltrating Immune Cells. Cancer Res. 77 (21), e108–e110. doi:10.1158/0008-5472.can-17-0307

PubMed Abstract | CrossRef Full Text | Google Scholar

Lin, W., Chen, Y., Wu, B., chen, Y., and Li, Z. (2021). Identification of the Pyroptosis-related P-rognostic G-ene S-ignature and the A-ssociated R-egulation axis in L-ung A-denocarcinoma. Cell Death Discov. 7 (1), 161. doi:10.1038/s41420-021-00557-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, C.-J., Hu, F.-F., Xia, M.-X., Han, L., Zhang, Q., and Guo, A.-Y. (2018). GSCALite: a Web Server for Gene Set Cancer Analysis. Bioinformatics 34 (21), 3771–3772. doi:10.1093/bioinformatics/bty411

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, R., Ju, C., Zhang, F., Tang, X., Yan, J., Sun, J., et al. (2020). LncRNA GSEC Promotes the Proliferation, Migration and Invasion by Sponging miR-588/EIF5A2 axis in Osteosarcoma. Biochem. Biophysical Res. Commun. 532 (2), 300–307. doi:10.1016/j.bbrc.2020.08.056

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, S.-H., Hsu, K.-W., Lai, Y.-L., Lin, Y.-F., Chen, F.-H., Peng, P.-H., et al. (2021). Systematic Identification of Clinically Relevant miRNAs for Potential miRNA-Based Therapy in Lung Adenocarcinoma. Mol. Ther. - Nucleic Acids 25, 1–10. doi:10.1016/j.omtn.2021.04.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Luo, L., Zheng, Y., Lin, Z., Li, X., Li, X., Li, M., et al. (2021). Identification of SHMT2 as a Potential Prognostic Biomarker and Correlating with Immune Infiltrates in Lung Adenocarcinoma. J. Immunol. Res. 2021, 6647122. doi:10.1155/2021/6647122

CrossRef Full Text | Google Scholar

Mao, C., Wang, X., Liu, Y., Wang, M., Yan, B., Jiang, Y., et al. (2018). A G3BP1-Interacting lncRNA Promotes Ferroptosis and Apoptosis in Cancer via Nuclear Sequestration of P53. Cancer Res. 78 (13), 3484–3496. doi:10.1158/0008-5472.CAN-17-3454

PubMed Abstract | CrossRef Full Text | Google Scholar

Mou, Y., Wang, J., Wu, J., He, D., Zhang, C., Duan, C., et al. (2019). Ferroptosis, a New Form of Cell Death: Opportunities and Challenges in Cancer. J. Hematol. Oncol. 12 (1), 34. doi:10.1186/s13045-019-0720-y

CrossRef Full Text | Google Scholar

Rivera, G. A., and Wakelee, H. (2016). Lung Cancer in Never Smokers. Adv. Exp. Med. Biol. 893, 43–57. doi:10.1007/978-3-319-24223-1_3

PubMed Abstract | CrossRef Full Text | Google Scholar

Schwartz, A. G., and Cote, M. L. (2016). Epidemiology of Lung Cancer. Adv. Exp. Med. Biol. 893, 21–41. doi:10.1007/978-3-319-24223-1_2

PubMed Abstract | CrossRef Full Text | Google Scholar

Song, X., Xie, Y., Kang, R., Hou, W., Sun, X., Epperly, M. W., et al. (2016). FANCD2 Protects against Bone Marrow Injury from Ferroptosis. Biochem. Biophysical Res. Commun. 480 (3), 443–449. doi:10.1016/j.bbrc.2016.10.068

CrossRef Full Text | Google Scholar

Song, X., Zhu, S., Chen, P., Hou, W., Wen, Q., Liu, J., et al. (2018). AMPK-mediated BECN1 Phosphorylation Promotes Ferroptosis by Directly Blocking System Xc- Activity. Curr. Biol. 28 (15), 2388–2399. doi:10.1016/j.cub.2018.05.094

PubMed Abstract | CrossRef Full Text | Google Scholar

Stockwell, B. R., and Jiang, X. (2020). The Chemistry and Biology of Ferroptosis. Cel Chem. Biol. 27 (4), 365–375. doi:10.1016/j.chembiol.2020.03.013

PubMed Abstract | CrossRef Full Text | Google Scholar

Su, Y., Zhao, B., Zhou, L., Zhang, Z., Shen, Y., Lv, H., et al. (2020). Ferroptosis, a Novel Pharmacological Mechanism of Anti-cancer Drugs. Cancer Lett. 483, 127–136. doi:10.1016/j.canlet.2020.02.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Mao, C., Ouyang, L., Liu, Y., Lai, W., Liu, N., et al. (2019). Long Noncoding RNA LINC00336 Inhibits Ferroptosis in Lung Cancer by Functioning as a Competing Endogenous RNA. Cell Death Differ 26 (11), 2329–2343. doi:10.1038/s41418-019-0304-y

PubMed Abstract | CrossRef Full Text | Google Scholar

Warde-Farley, D., Donaldson, S. L., Comes, O., Zuberi, K., Badrawi, R., Chao, P., et al. (2010). The GeneMANIA Prediction Server: Biological Network Integration for Gene Prioritization and Predicting Gene Function. Nucleic Acids Res. 38, W214–W220. doi:10.1093/nar/gkq537

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, F., Huang, W., Yang, L., and Xu, F. (2019). MicroRNA-101-3p Regulates Gastric Cancer Cell Proliferation, Invasion and Apoptosis by Targeting PIM 1 Expression. Cel Mol Biol (Noisy-le-grand) 65 (7), 118–122. doi:10.14715/cmb/2019.65.7.20

CrossRef Full Text | Google Scholar

Xiong, Q., Jiang, L., Liu, K., Jiang, X., Liu, B., Shi, Y., et al. (2021). miR-133b Targets NCAPH to Promote β-catenin Degradation and Reduce Cancer Stem Cell Maintenance in Non-small Cell Lung Cancer. Sig Transduct Target. Ther. 6 (1), 252. doi:10.1038/s41392-021-00555-x

CrossRef Full Text | Google Scholar

Xu, G., Wang, H., Li, X., Huang, R., and Luo, L. (2021). Recent Progress on Targeting Ferroptosis for Cancer Therapy. Biochem. Pharmacol. 190, 114584. doi:10.1016/j.bcp.2021.114584

PubMed Abstract | CrossRef Full Text | Google Scholar

Xu, P., Jiang, L., Yang, Y., Wu, M., Liu, B., Shi, Y., et al. (2020). PAQR4 Promotes Chemoresistance in Non-small Cell Lung Cancer through Inhibiting Nrf2 Protein Degradation. Theranostics 10 (8), 3767–3778. doi:10.7150/thno.43142

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, D., Wang, J., Li, C., Shi, L., and Zhang, M. (2021). Ferroptosis-related Gene Model to Predict Overall Survival of Papillary Thyroid Carcinoma. Am. J. Otolaryngol. 42 (6), 103163. doi:10.1016/j.amjoto.2021.103163

CrossRef Full Text | Google Scholar

Yang, W., Soares, J., Greninger, P., Edelman, E. J., Lightfoot, H., Forbes, S., et al. (2013). Genomics of Drug Sensitivity in Cancer (GDSC): a Resource for Therapeutic Biomarker Discovery in Cancer Cells. Nucleic Acids Res. 41, D955–D961. doi:10.1093/nar/gks1111

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, W. S., and Stockwell, B. R. (2016). Ferroptosis: Death by Lipid Peroxidation. Trends Cel Biol. 26 (3), 165–176. doi:10.1016/j.tcb.2015.10.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, J., Du, C., Zhang, L., Wang, Y., Zhang, Y., and Li, J. (2021). lncRNA GSEC Promotes the Progression of Triple Negative Breast Cancer (TNBC) by Targeting the miR-202-5p/AXL Axis. Ott Vol. 14, 2747–2759. doi:10.2147/ott.s293832

CrossRef Full Text | Google Scholar

Keywords: NSCLC, ferroptosis, ceRNA, immune cell infiltration, cell proliferation, cell migration

Citation: Jiang X, Yuan Y, Tang L, Wang J, Zhang D and Duan L (2022) Systematic Analysis and Validation of the Prognosis, Immunological Role and Biology Function of the Ferroptosis-Related lncRNA GSEC/miRNA-101-3p/CISD1 Axis in Lung Adenocarcinoma. Front. Mol. Biosci. 8:793732. doi: 10.3389/fmolb.2021.793732

Received: 12 October 2021; Accepted: 22 December 2021;
Published: 07 March 2022.

Edited by:

Matteo Becatti, University of Firenze, Italy

Reviewed by:

Anirban Chakraborty, University of Texas Medical Branch at Galveston, United States
Qiang Guo, Hubei University of Medicine, China
Rui-Qing Zhai, Harbin Medical University, China

Copyright © 2022 Jiang, Yuan, Tang, Wang, Zhang and Duan. 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: Lincan Duan, duanmumuhuosan@163.com

These authors have contributed equally to this work

Disclaimer: 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.