Search

Background: This study was performed to determine the biological processes in which NKX2-1 is involved and thus its role in the development of lung squamous cell carcinoma (LUSC) toward improving the prognosis and treatment of LUSC. Methods: Raw RNA sequencing (RNA-seq) data of LUSC from The Cancer Genome Atlas (TCGA) were used in bioinformatics analysis to characterize NKX2-1 expression levels in tumor and normal tissues. Survival analysis of Kaplan–Meier curve, the time-dependent receiver operating characteristic (ROC) curve, and a nomogram were used to analyze the prognosis value of NKX2-1 for LUSC in terms of overall survival (OS) and progression-free survival (PFS). Then, differentially expressed genes (DEGs) were identified, and Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Gene Set Enrichment Analysis (GSEA) were used to clarify the biological mechanisms potentially involved in the development of LUSC. Moreover, the correlation between the NKX2-1 expression level and tumor mutation burden (TMB), tumor microenvironment (TME), and immune cell infiltration revealed that NKX2-1 participates in the development of LUSC. Finally, we studied the effects of NKX2-1 on drug therapy. To validate the protein and gene expression levels of NKX2-1 in LUSC, we employed immunohistochemistry(IHC) datasets, The Gene Expression Omnibus (GEO) database, and qRT-PCR analysis. Results: NKX2-1 expression levels were significantly lower in LUSC than in normal lung tissue. It significantly differed in gender, stage and N classification. The survival analysis revealed that high expression of NKX2-1 had shorter OS and PFS in LUSC. The multivariate Cox regression hazard model showed the NKX2-1 expression as an independent prognostic factor. Then, the nomogram predicted LUSC prognosis. There are 51 upregulated DEGs and 49 downregulated DEGs in the NKX2-1 high-level groups. GO, KEGG and GSEA analysis revealed that DEGs were enriched in cell cycle and DNA replication.The TME results show that NKX2-1 expression was positively associated with mast cells resting, neutrophils, monocytes, T cells CD4 memory resting, and M2 macrophages but negatively associated with M1 macrophages. The TMB correlated negatively with NKX2-1 expression. The pharmacotherapy had great sensitivity in the NKX2-1 low-level group, the immunotherapy is no significant difference in the NKX2-1 low-level and high-level groups. The analysis of GEO data demonstrated concurrence with TCGA results. IHC revealed NKX2-1 protein expression in tumor tissues of both LUAD and LUSC. Meanwhile qRT-PCR analysis indicated a significantly lower NKX2-1 expression level in LUSC compared to LUAD. These qRT-PCR findings were consistent with co-expression analysis of NKX2-1. Conclusion: We conclude that NKX2-1 is a potential biomarker for prognosis and treatment LUSC. A new insights of NKX2-1 in LUSC is still needed further research. [ABSTRACT FROM AUTHOR]

Purpose: To discover the biological function and potential mechanism of LINC01936 in the development of lung squamous cell carcinoma (LUSC). Methods: Transcriptome data of LUSC from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases were used to analyze the differentially expressed lncRNAs in LUSC and normal tissues by R "DEseq2", "edgeR" and "limma" packages. The subcellular localization of LINC01936 was predicted by lncLocator. Cell proliferation and apoptosis were measured by CCK-8, MTT assay and Hoechst fluorescence staining. The migration and invasion were detected by Transwell assay. The function and pathway enrichment analysis were performed by Gene Ontology (GO) terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and gene set variation analysis (GSVA). The downstream targets of LINC01936 were predicted using RNA-Protein Interaction Prediction (RPISeq) program. The effect of LINC01936 on tumor immune infiltration was analyzed using Pearson Correlation Analysis using R "ggpubr" package. Results: Based on the gene expression data of LUSC from TCGA database, 1,603, 1,702 and 529 upregulated and 536, 436 and 630 downregulated lncRNAs were obtained by DEseq2, edgeR and limma programs, respectively. For GSE88862 dataset, we acquired 341 differentially expressed lncRNAs (206 upregulated and 135 downregulated). Venn plot for the intersection of above differential expressed lncRNAs showed that there were 29 upregulated and 23 downregulated genes. LINC01936 was one of downregulated lncRNAs in LUSC tissues. The biological analysis showed that the overexpression of LINC01936 significantly reduced proliferation, migration and invasion of LUSC cells, and promoted cell apoptosis. The knockdown of LINC01936 promoted cell proliferation and metastasis. Pathway and GSVA analysis indicated that LINC01936 might participated in DNA repair, complement, cell adhesion and EMT, etc. LINC01936 was predicted to interact with TCF21, AOC3, RASL12, MEOX2 or HSPB7, which are involved in EMT and PI3KAKT-MTOR pathway, etc. The expression of LINC01936 was also positively correlated with the infiltrating immune cells in LUSC. [ABSTRACT FROM AUTHOR]

Background This study was performed to determine the biological processes in which NKX2-1 is involved and thus its role in the development of lung squamous cell carcinoma (LUSC) toward improving the prognosis and treatment of LUSC. Methods Raw RNA sequencing (RNA-seq) data of LUSC from The Cancer Genome Atlas (TCGA) were used in bioinformatics analysis to characterize NKX2-1 expression levels in tumor and normal tissues. Survival analysis of Kaplan–Meier curve, the time-dependent receiver operating characteristic (ROC) curve, and a nomogram were used to analyze the prognosis value of NKX2-1 for LUSC in terms of overall survival (OS) and progression-free survival (PFS). Then, differentially expressed genes (DEGs) were identified, and Kyoto Encyclopedia of Genes and Genomes (KEGG), Gene Ontology (GO), and Gene Set Enrichment Analysis (GSEA) were used to clarify the biological mechanisms potentially involved in the development of LUSC. Moreover, the correlation between the NKX2-1 expression level and tumor mutation burden (TMB), tumor microenvironment (TME), and immune cell infiltration revealed that NKX2-1 participates in the development of LUSC. Finally, we studied the effects of NKX2-1 on drug therapy. To validate the protein and gene expression levels of NKX2-1 in LUSC, we employed immunohistochemistry(IHC) datasets, The Gene Expression Omnibus (GEO) database, and qRT-PCR analysis. Results NKX2-1 expression levels were significantly lower in LUSC than in normal lung tissue. It significantly differed in gender, stage and N classification. The survival analysis revealed that high expression of NKX2-1 had shorter OS and PFS in LUSC. The multivariate Cox regression hazard model showed the NKX2-1 expression as an independent prognostic factor. Then, the nomogram predicted LUSC prognosis. There are 51 upregulated DEGs and 49 downregulated DEGs in the NKX2-1 high-level groups. GO, KEGG and GSEA analysis revealed that DEGs were enriched in cell cycle and DNA replication.The TME results show that NKX2-1 expression was positively associated with mast cells resting, neutrophils, monocytes, T cells CD4 memory resting, and M2 macrophages but negatively associated with M1 macrophages. The TMB correlated negatively with NKX2-1 expression. The pharmacotherapy had great sensitivity in the NKX2-1 low-level group, the immunotherapy is no significant difference in the NKX2-1 low-level and high-level groups. The analysis of GEO data demonstrated concurrence with TCGA results. IHC revealed NKX2-1 protein expression in tumor tissues of both LUAD and LUSC. Meanwhile qRT-PCR analysis indicated a significantly lower NKX2-1 expression level in LUSC compared to LUAD. These qRT-PCR findings were consistent with co-expression analysis of NKX2-1. Conclusion We conclude that NKX2-1 is a potential biomarker for prognosis and treatment LUSC. A new insights of NKX2-1 in LUSC is still needed further research.

Background. Worldwide, lung squamous cell carcinoma (LUSC) has wreaked havoc on humanity. Matrix metallopeptidase 12 (MMP12) plays an essential role in a variety of cancers. This study aimed to reveal the expression, clinical significance, and potential molecular mechanisms of MMP12 in LUSC. Methods. There were 2,738 messenger RNA (mRNA) samples from several multicenter databases used to detect MMP12 expression in LUSC, and 125 tissue samples were validated by immunohistochemistry (IHC) experiments. Receiver operator charac- teristic (ROC) curves, Kaplan-Meier curves, and univariate and multivariate Cox regression analyses were used to assess the clinical value of MMP12 in LUSC. The potential molecular mechanisms of MMP12 were explored by gene enrichment analysis and immune correlation analysis. Furthermore, single-cell sequencing was used to determine the distribution of MMP12 in multiple tumor microenvironment cells. Results. MMP12 was significantly overexpressed at the mRNA level (p<0:05, SMD D 3.13, 95% CI [2.51 -3.75]), which was verified at the protein level (p<0:001) by internal IHC experiments. MMP12 expression could be used to differentiate LUSC samples from normal samples, and overexpression of MMP12 itself implied a worse clinical prognosis and higher levels of immune cell infiltration in LUSC patients. MMP12 was involved in cancer development and progression through two immune-related signaling pathways. The high expression of MMP12 in LUSC might act as an antigen- presenting cell-associated tumor neoantigen and activate the body's immune response. Conclusions. MMP12 expression is upregulated in LUSC and high expression of MMP12 serves as a risk factor for LUSC patients. MMP12 may be involved in cancer development by participating in immune-related signaling pathways and elevating the level of immune cell infiltration. [ABSTRACT FROM AUTHOR]

Background: Recurrence is a risk factor for the prognosis of lung squamous carcinoma (LUSC). DNA methylation levels of RNAs are also associated with LUSC prognosis. This study aimed to construct a prognostic model with high performance in predicting LUSC prognosis using the methylation levels of lncRNAs and genes. Methods: The differentially expressed RNAs (DERs) and differentially methylated RNAs (DMRs) between the recurrent and non-recurrent LUSC tissues in The Cancer Genome Atlas (TCGA; training dataset) were identified. Weighted correlation network analysis was performed to identify co-methylation networks. Differentially methylated genes and lncRNAs with opposite expression-methylation levels were used for the screening of prognosis-associated RNAs. The prognostic model was constructed and its performance was validated in the GSE39279 dataset. Results: A total of 664 DERs and 981 DMRs (including 972 genes) in recurrent LUSC tissues were identified. Three co-methylation modules, including 226 differentially methylated genes, were significantly associated with LUSC. Among prognosis-associated RNAs, 18 DERs/DMRs with opposite methylation-expression levels were included in the methylation prognostic risk model. LUSC patients with high risk scores had a poor prognosis compared with patients who had low risk scores (TCGA: HR = 3.856, 95% CI [2.297-6.471]; GSE39279: HR = 3.040, 95% CI [1.435-6.437]). This model had a high accuracy in predicting the prognosis (AUC = 0.903 and 0.800, respectively), equivalent to the nomogram model inclusive of clinical variables. Conclusions: Referring to the methylation levels of the 16-RNAs might help to predict the survival outcomes in LUSC. [ABSTRACT FROM AUTHOR]

As an innate feature of human beings, gender differences have an influence on various biological phenotypes, yet it does not attract enough attention in genomics studies. The prognosis of multiple carcinomas usually exhibits a favorable ending for female patients, but the neglect of gender differences can cause serious bias in survival analysis. Enhancer RNAs (eRNAs) are mostly downstream of androgens or estrogen. The present study was aimed to screen eRNAs in patients with non-smallcell lung cancer. The findings revealed that eRNA TBX5-AS1 was expressed differently between female and male patients. Meanwhile, its prognostic significance appeared only in male patients with squamous cell carcinoma (SCC) type. The Gene Set Enrichment Analysis proved that the expression level of TBX5-AS1 increased following the activation of the androgen signaling pathway. In pan-cancer analysis, the prognostic prediction based on gender grouping obtained more meaningful results, and the synergy between TBX5-AS1 and its homologous target was more consistent. Furthermore, immunity variations between sexes prompted us to explore the role that TBX5-AS1 played in tumor microenvironment and immunotherapy. The robust evidence proved that male patients with high expression of TBX5-AS1 possessed a malignant immune microenvironment and urgently needed immune checkpoint inhibitor treatment. In conclusion, TBX5-AS1 may be one of the strongest candidates to predict prognosis for male patients with SCC and provide a reference for immunotherapy. [ABSTRACT FROM AUTHOR]

Background Squamous cell carcinomas of the lung are an extremely common and deadly form of non-small cell lung cancers. Clinical management of the disease is dependent on staging and metastatic status. Metastasis to the lymph node is especially crucial to diagnose as it occurs at an earlier stage. However, lymphadenectomies are invasive and tumor cells may be overlooked during evaluation.There are limited approved biomarkers for predicting lymph node metastasis with squamous cell carcinomas of the lung (LSCC). Methods Genome data of 60 tumor-adjacent samples were downloaded from Genome Expression Omnibus. We identified over-expressed HUB genes using Cytoscape as key prognostic markers. The selected markers were further evaluated based on gene ontology and overall expression levels compared to normal tissue using The Cancer Genome Atlas. We further validated these results using clinical biopsy tissue taken from squamous cell carcinoma patients. Results Analysis of the genome expression data resulted in 13 relevant hub genes that were differentially expressed in cancerous samples. All of these genes are associated with collagen biosynthesis within the tumor microenvironment. We chose Collagen Type 1 Alpha 1 (COL1A1) as the most relevant prognostic marker due to its high number of pathway connections and over expression in the tumor microenvironment compared to the other 12 genes. Additionally, based on analysis of The Cancer Genome Atlas, tumors with higher levels of COL1A1 expression are associated with poorer overall survival. Finally, evaluation of clinical biopsy samples suggests that overexpression of COL1A1 in the LSCC microenvironment highly correlates with lymph node metastasis. These results suggest COL1A1 is a clinically relevant marker that should be used to justify lymphadenectomies.

Background Lung squamous cell carcinoma (LUSC) is a major subtype of lung cancer with limited therapeutic options and poor clinical prognosis. Methods Three datasets (GSE19188, GSE33532 and GSE33479) were obtained from the gene expression omnibus (GEO) database. Differentially expressed genes (DEGs) between LUSC and normal tissues were identified by GEO2R, and functional analysis was employed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) online tool. Protein–protein interaction (PPI) and hub genes were identified via the Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape software. Hub genes were further validated in The Cancer Genome Atlas (TCGA) database. Subsequently, survival analysis was performed using the Kapla–Meier curve and Cox progression analysis. Based on univariate and multivariate Cox progression analysis, a gene signature was established to predict overall survival. Receiver operating characteristic curve was used to evaluate the prognostic value of the model. Results A total of 116 up-regulated genes and 84 down-regulated genes were identified. These DEGs were mainly enriched in the two pathways: cell cycle and p53 signaling way. According to the degree of protein nodes in the PPI network, 10 hub genes were identified. The mRNA expression levels of the 10 hub genes in LUSC were also significantly up-regulated in the TCGA database. Furthermore, a novel seven-gene signature (FLRT3, PPP2R2C, MMP3, MMP12, CAPN8, FILIP1 and SPP1) from the DEGs was constructed and acted as a significant and independent prognostic signature for LUSC. Conclusions The 10 hub genes might be tightly correlated with LUSC progression. The seven-gene signature might be an independent biomarker with a significant predictive value in LUSC overall survival.

Purpose There is plenty of evidence showing that autophagy plays an important role in the biological process of cancer. The purpose of this study was to establish a novel autophagy-related prognostic marker for lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC). Methods The mRNA microarray and clinical data in The Cancer Genome Atlas (TCGA) were analyzed by using a univariate Cox proportional regression model to select candidate autophagy-related prognostic genes. Bioinformatics analysis of gene function using the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) platforms was performed. A multivariate Cox proportional regression model helped to develop a prognostic signature from the pool of candidate genes. On the basis of this prognostic signature, we could divide LUAD and LUSC patients into high-risk and low-risk groups. Further survival analysis demonstrated that high-risk patients had significantly shorter disease-free survival (DFS) than low-risk patients. The signature which contains six autophagy-related genes (EIF4EBP1, TP63, BNIP3, ATIC, ERO1A and FADD) showed good performance for predicting the survival of LUAD and LUSC patients by having a better Area Under Curves (AUC) than other clinical parameters. Its efficacy was also validated by data from the Gene Expression Omnibus (GEO) database. Conclusion Collectively, the prognostic signature we proposed is a promising biomarker for monitoring the outcomes of LUAD and LUSC.

Background Long noncoding RNAs (lncRNAs) play a role in the formation, development, and prognosis of various cancers. Our study aimed to identify prognostic-related lncRNAs in lung squamous cell carcinoma (LUSC), which may provide new perspectives for individualized treatment of patients. Materials and Methods The RNA sequencing (lncRNA, microRNA (miRNA), mRNA) data and clinical information related to LUSC were obtained from The Cancer Genome Atlas (TCGA) database. Differentially expressed RNA sequences were used to construct the competitive endogenous RNA (ceRNA) network. In present study, we mainly used two prognostic verification methods, Cox analysis and survival analysis, to identify the prognostic relevance of specific lncRNAs and construct prognostic model of lncRNA. Results Datasets on 551 samples of lncRNA and mRNA and 523 miRNA samples were retrieved from the TCGA database. Analysis of the normal and LUSC samples identified 170 DElncRNAs, 331 DEmiRNAs, and 417 DEmRNAs differentially expressed RNAs. The ceRNA network contained 27 lncRNAs, 43 miRNAs, and 11 mRNAs. Furthermore, we identified seven specific lncRNAs (ERVH48-1, HCG9, SEC62-AS1, AC022148.1, LINC00460, C5orf17, LINC00261) as potential prognostic factors after correlation analysis, and five of the seven lncRNAs (AC022148.1, HCG9, LINC00460, C5orf17, LINC00261) constructed a prognostic model of LUSC. Conclusion In present study, we identified seven lncRNAs in the ceRNA network that are associated with potential prognosis in LUSC patients, and constructed a prognostic model of LUSC which can be used to assess the prognosis risk of clinical patients. Further biological experiments are needed to elucidate the specific molecular mechanisms underlying them.

Purpose: To discover the biological function and potential mechanism of LINC01936 in the development of lung squamous cell carcinoma (LUSC)., Methods: Transcriptome data of LUSC from The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases were used to analyze the differentially expressed lncRNAs in LUSC and normal tissues by R "DEseq2", "edgeR" and "limma" packages. The subcellular localization of LINC01936 was predicted by lncLocator. Cell proliferation and apoptosis were measured by CCK-8, MTT assay and Hoechst fluorescence staining. The migration and invasion were detected by Transwell assay. The function and pathway enrichment analysis were performed by Gene Ontology (GO) terms, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis and gene set variation analysis (GSVA). The downstream targets of LINC01936 were predicted using RNA-Protein Interaction Prediction (RPISeq) program. The effect of LINC01936 on tumor immune infiltration was analyzed using Pearson Correlation Analysis using R "ggpubr" package., Results: Based on the gene expression data of LUSC from TCGA database, 1,603, 1,702 and 529 upregulated and 536, 436 and 630 downregulated lncRNAs were obtained by DEseq2, edgeR and limma programs, respectively. For GSE88862 dataset, we acquired 341 differentially expressed lncRNAs (206 upregulated and 135 downregulated). Venn plot for the intersection of above differential expressed lncRNAs showed that there were 29 upregulated and 23 downregulated genes. LINC01936 was one of downregulated lncRNAs in LUSC tissues. The biological analysis showed that the overexpression of LINC01936 significantly reduced proliferation, migration and invasion of LUSC cells, and promoted cell apoptosis. The knockdown of LINC01936 promoted cell proliferation and metastasis. Pathway and GSVA analysis indicated that LINC01936 might participated in DNA repair, complement, cell adhesion and EMT, etc . LINC01936 was predicted to interact with TCF21, AOC3, RASL12, MEOX2 or HSPB7, which are involved in EMT and PI3K-AKT-MTOR pathway, etc . The expression of LINC01936 was also positively correlated with the infiltrating immune cells in LUSC., Conclusions: LINC01936 is downregulated in LUSC. LINC01936 affected proliferation, migration and invasion of LUSC cells probably by EMT and immune infiltration, which might serve as a new target for the treatment of LUSC., Competing Interests: The authors declare that they have no competing interests., (© 2023 Tian et al.)

The etiology of cancer includes aberrant cellular homeostasis where a compromised RNA regulatory network is a prominent contributing factor. In particular, noncoding RNAs including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) were recently shown to play important roles in the initiation, progression, and metastasis of human cancers. Nonetheless, a mechanistic understanding of noncoding RNA functions in lung squamous cell carcinoma (LUSC) is lacking. To fill this critical gap in knowledge, we obtained mRNA, miRNA, and lncRNA expression data on patients with LUSC from the updated Cancer Genome Atlas (TCGA) database (2016). We successfully identified 3,366 mRNAs, 79 miRNAs, and 151 lncRNAs as key contributing factors of a high risk of LUSC. Furthermore, we hypothesized that the lncRNA–miRNA–mRNA regulatory axis positively correlates with LUSC and constructed a competitive endogenous RNA (ceRNA) network of LUSC by targeting interrelations with significantly aberrant expression data between miRNA and mRNA or lncRNA. Six ceRNAs (PLAU, miR-31-5p, miR-455-3p, FAM83A-AS1, MIR31HG, and MIR99AHG) significantly correlated with survival (P < 0.05). Finally, real-time quantitative PCR analysis showed that PLAU is significantly upregulated in SK-MES-1 cells compared with 16-BBE-T cells. Taken together, our findings represent new knowledge for a better understanding the ceRNA network in LUSC biology and pave the way to improved diagnosis and prognosis of LUSC.

Background: Worldwide, lung squamous cell carcinoma (LUSC) has wreaked havoc on humanity. Matrix metallopeptidase 12 ( MMP12 ) plays an essential role in a variety of cancers. This study aimed to reveal the expression, clinical significance, and potential molecular mechanisms of MMP12 in LUSC., Methods: There were 2,738 messenger RNA (mRNA) samples from several multicenter databases used to detect MMP12 expression in LUSC, and 125 tissue samples were validated by immunohistochemistry (IHC) experiments. Receiver operator characteristic (ROC) curves, Kaplan-Meier curves, and univariate and multivariate Cox regression analyses were used to assess the clinical value of MMP12 in LUSC. The potential molecular mechanisms of MMP12 were explored by gene enrichment analysis and immune correlation analysis. Furthermore, single-cell sequencing was used to determine the distribution of MMP12 in multiple tumor microenvironment cells., Results: MMP12 was significantly overexpressed at the mRNA level ( p  < 0.05, SMD = 3.13, 95% CI [2.51-3.75]), which was verified at the protein level ( p  < 0.001) by internal IHC experiments. MMP12 expression could be used to differentiate LUSC samples from normal samples, and overexpression of MMP12 itself implied a worse clinical prognosis and higher levels of immune cell infiltration in LUSC patients. MMP12 was involved in cancer development and progression through two immune-related signaling pathways. The high expression of MMP12 in LUSC might act as an antigen-presenting cell-associated tumor neoantigen and activate the body's immune response., Conclusions: MMP12 expression is upregulated in LUSC and high expression of MMP12 serves as a risk factor for LUSC patients. MMP12 may be involved in cancer development by participating in immune-related signaling pathways and elevating the level of immune cell infiltration., Competing Interests: Gang Chen is an Academic Editor for PeerJ., (©2023 Zhang et al.)

Background: Recurrence is a risk factor for the prognosis of lung squamous carcinoma (LUSC). DNA methylation levels of RNAs are also associated with LUSC prognosis. This study aimed to construct a prognostic model with high performance in predicting LUSC prognosis using the methylation levels of lncRNAs and genes., Methods: The differentially expressed RNAs (DERs) and differentially methylated RNAs (DMRs) between the recurrent and non-recurrent LUSC tissues in The Cancer Genome Atlas (TCGA; training dataset) were identified. Weighted correlation network analysis was performed to identify co-methylation networks. Differentially methylated genes and lncRNAs with opposite expression-methylation levels were used for the screening of prognosis-associated RNAs. The prognostic model was constructed and its performance was validated in the GSE39279 dataset., Results: A total of 664 DERs and 981 DMRs (including 972 genes) in recurrent LUSC tissues were identified. Three co-methylation modules, including 226 differentially methylated genes, were significantly associated with LUSC. Among prognosis-associated RNAs, 18 DERs/DMRs with opposite methylation-expression levels were included in the methylation prognostic risk model. LUSC patients with high risk scores had a poor prognosis compared with patients who had low risk scores (TCGA: HR = 3.856, 95% CI [2.297-6.471]; GSE39279: HR = 3.040, 95% CI [1.435-6.437]). This model had a high accuracy in predicting the prognosis (AUC = 0.903 and 0.800, respectively), equivalent to the nomogram model inclusive of clinical variables., Conclusions: Referring to the methylation levels of the 16-RNAs might help to predict the survival outcomes in LUSC., Competing Interests: The authors declare that they have no competing interests., (© 2022 Wang et al.)

As an innate feature of human beings, gender differences have an influence on various biological phenotypes, yet it does not attract enough attention in genomics studies. The prognosis of multiple carcinomas usually exhibits a favorable ending for female patients, but the neglect of gender differences can cause serious bias in survival analysis. Enhancer RNAs (eRNAs) are mostly downstream of androgens or estrogen. The present study was aimed to screen eRNAs in patients with non-small-cell lung cancer. The findings revealed that eRNA TBX5-AS1 was expressed differently between female and male patients. Meanwhile, its prognostic significance appeared only in male patients with squamous cell carcinoma (SCC) type. The Gene Set Enrichment Analysis proved that the expression level of TBX5-AS1 increased following the activation of the androgen signaling pathway. In pan-cancer analysis, the prognostic prediction based on gender grouping obtained more meaningful results, and the synergy between TBX5-AS1 and its homologous target was more consistent. Furthermore, immunity variations between sexes prompted us to explore the role that TBX5-AS1 played in tumor microenvironment and immunotherapy. The robust evidence proved that male patients with high expression of TBX5-AS1 possessed a malignant immune microenvironment and urgently needed immune checkpoint inhibitor treatment. In conclusion, TBX5-AS1 may be one of the strongest candidates to predict prognosis for male patients with SCC and provide a reference for immunotherapy., Competing Interests: The authors declare that they have no competing interests., (© 2021 Yan et al.)

Background: Squamous cell carcinomas of the lung are an extremely common and deadly form of non-small cell lung cancers. Clinical management of the disease is dependent on staging and metastatic status. Metastasis to the lymph node is especially crucial to diagnose as it occurs at an earlier stage. However, lymphadenectomies are invasive and tumor cells may be overlooked during evaluation.There are limited approved biomarkers for predicting lymph node metastasis with squamous cell carcinomas of the lung (LSCC)., Methods: Genome data of 60 tumor-adjacent samples were downloaded from Genome Expression Omnibus. We identified over-expressed HUB genes using Cytoscape as key prognostic markers. The selected markers were further evaluated based on gene ontology and overall expression levels compared to normal tissue using The Cancer Genome Atlas. We further validated these results using clinical biopsy tissue taken from squamous cell carcinoma patients., Results: Analysis of the genome expression data resulted in 13 relevant hub genes that were differentially expressed in cancerous samples. All of these genes are associated with collagen biosynthesis within the tumor microenvironment. We chose Collagen Type 1 Alpha 1 (COL1A1) as the most relevant prognostic marker due to its high number of pathway connections and over expression in the tumor microenvironment compared to the other 12 genes. Additionally, based on analysis of The Cancer Genome Atlas, tumors with higher levels of COL1A1 expression are associated with poorer overall survival. Finally, evaluation of clinical biopsy samples suggests that overexpression of COL1A1 in the LSCC microenvironment highly correlates with lymph node metastasis. These results suggest COL1A1 is a clinically relevant marker that should be used to justify lymphadenectomies., Competing Interests: The authors declare there are no competing interests., (©2020 Dong et al.)

Background: Lung squamous cell carcinoma (LUSC) is a major subtype of lung cancer with limited therapeutic options and poor clinical prognosis., Methods: Three datasets (GSE19188, GSE33532 and GSE33479) were obtained from the gene expression omnibus (GEO) database. Differentially expressed genes (DEGs) between LUSC and normal tissues were identified by GEO2R, and functional analysis was employed using the Database for Annotation, Visualization and Integrated Discovery (DAVID) online tool. Protein-protein interaction (PPI) and hub genes were identified via the Search Tool for the Retrieval of Interacting Genes (STRING) and Cytoscape software. Hub genes were further validated in The Cancer Genome Atlas (TCGA) database. Subsequently, survival analysis was performed using the Kapla-Meier curve and Cox progression analysis. Based on univariate and multivariate Cox progression analysis, a gene signature was established to predict overall survival. Receiver operating characteristic curve was used to evaluate the prognostic value of the model., Results: A total of 116 up-regulated genes and 84 down-regulated genes were identified. These DEGs were mainly enriched in the two pathways: cell cycle and p53 signaling way. According to the degree of protein nodes in the PPI network, 10 hub genes were identified. The mRNA expression levels of the 10 hub genes in LUSC were also significantly up-regulated in the TCGA database. Furthermore, a novel seven-gene signature (FLRT3, PPP2R2C, MMP3, MMP12, CAPN8, FILIP1 and SPP1) from the DEGs was constructed and acted as a significant and independent prognostic signature for LUSC., Conclusions: The 10 hub genes might be tightly correlated with LUSC progression. The seven-gene signature might be an independent biomarker with a significant predictive value in LUSC overall survival., Competing Interests: The authors declare that they have no competing interests., (© 2020 Ma et al.)

Purpose: There is plenty of evidence showing that autophagy plays an important role in the biological process of cancer. The purpose of this study was to establish a novel autophagy-related prognostic marker for lung adenocarcinoma (LUAD) and lung squamous cell carcinoma (LUSC)., Methods: The mRNA microarray and clinical data in The Cancer Genome Atlas (TCGA) were analyzed by using a univariate Cox proportional regression model to select candidate autophagy-related prognostic genes. Bioinformatics analysis of gene function using the Gene Ontology (GO) and the Kyoto Encyclopedia of Genes and Genomes (KEGG) platforms was performed. A multivariate Cox proportional regression model helped to develop a prognostic signature from the pool of candidate genes. On the basis of this prognostic signature, we could divide LUAD and LUSC patients into high-risk and low-risk groups. Further survival analysis demonstrated that high-risk patients had significantly shorter disease-free survival (DFS) than low-risk patients. The signature which contains six autophagy-related genes (EIF4EBP1, TP63, BNIP3, ATIC, ERO1A and FADD) showed good performance for predicting the survival of LUAD and LUSC patients by having a better Area Under Curves (AUC) than other clinical parameters. Its efficacy was also validated by data from the Gene Expression Omnibus (GEO) database., Conclusion: Collectively, the prognostic signature we proposed is a promising biomarker for monitoring the outcomes of LUAD and LUSC., Competing Interests: The authors declare there are no competing interests., (©2020 Zhu et al.)

Background: Long noncoding RNAs (lncRNAs) play a role in the formation, development, and prognosis of various cancers. Our study aimed to identify prognostic-related lncRNAs in lung squamous cell carcinoma (LUSC), which may provide new perspectives for individualized treatment of patients., Materials and Methods: The RNA sequencing (lncRNA, microRNA (miRNA), mRNA) data and clinical information related to LUSC were obtained from The Cancer Genome Atlas (TCGA) database. Differentially expressed RNA sequences were used to construct the competitive endogenous RNA (ceRNA) network. In present study, we mainly used two prognostic verification methods, Cox analysis and survival analysis, to identify the prognostic relevance of specific lncRNAs and construct prognostic model of lncRNA., Results: Datasets on 551 samples of lncRNA and mRNA and 523 miRNA samples were retrieved from the TCGA database. Analysis of the normal and LUSC samples identified 170 DElncRNAs, 331 DEmiRNAs, and 417 DEmRNAs differentially expressed RNAs. The ceRNA network contained 27 lncRNAs, 43 miRNAs, and 11 mRNAs. Furthermore, we identified seven specific lncRNAs (ERVH48-1, HCG9, SEC62-AS1, AC022148.1, LINC00460, C5orf17, LINC00261) as potential prognostic factors after correlation analysis, and five of the seven lncRNAs (AC022148.1, HCG9, LINC00460, C5orf17, LINC00261) constructed a prognostic model of LUSC., Conclusion: In present study, we identified seven lncRNAs in the ceRNA network that are associated with potential prognosis in LUSC patients, and constructed a prognostic model of LUSC which can be used to assess the prognosis risk of clinical patients. Further biological experiments are needed to elucidate the specific molecular mechanisms underlying them., Competing Interests: The authors declare that they have no competing interests., (© 2019 Qi et al.)

The etiology of cancer includes aberrant cellular homeostasis where a compromised RNA regulatory network is a prominent contributing factor. In particular, noncoding RNAs including microRNAs (miRNAs) and long noncoding RNAs (lncRNAs) were recently shown to play important roles in the initiation, progression, and metastasis of human cancers. Nonetheless, a mechanistic understanding of noncoding RNA functions in lung squamous cell carcinoma (LUSC) is lacking. To fill this critical gap in knowledge, we obtained mRNA, miRNA, and lncRNA expression data on patients with LUSC from the updated Cancer Genome Atlas (TCGA) database (2016). We successfully identified 3,366 mRNAs, 79 miRNAs, and 151 lncRNAs as key contributing factors of a high risk of LUSC. Furthermore, we hypothesized that the lncRNA-miRNA-mRNA regulatory axis positively correlates with LUSC and constructed a competitive endogenous RNA (ceRNA) network of LUSC by targeting interrelations with significantly aberrant expression data between miRNA and mRNA or lncRNA. Six ceRNAs (PLAU, miR-31-5p, miR-455-3p, FAM83A-AS1, MIR31HG, and MIR99AHG) significantly correlated with survival ( P < 0.05). Finally, real-time quantitative PCR analysis showed that PLAU is significantly upregulated in SK-MES-1 cells compared with 16-BBE-T cells. Taken together, our findings represent new knowledge for a better understanding the ceRNA network in LUSC biology and pave the way to improved diagnosis and prognosis of LUSC., Competing Interests: The authors declare that they have no competing interests.