Genome-scale analysis of DNA methylation in lung adenocarcinoma and integration with mRNA expression

Integration of osimertinib-targeted EGFR gene-associated differential gene expression in constructing a prognostic model for lung adenocarcinoma

Lung adenocarcinoma (LUAD) is one of the deadliest cancers. Epidermal growth factor receptor-tyrosine kinase inhibitor (EGFR-TKI)-targeted therapy is an important approach for treating LUAD. However, the development of acquired resistance poses a serious clinical challenge. Our objective was to explore the differentially expressed genes (DEGs) associated with EGFR and detect biomarkers for diagnosing and treating osimertinib resistance in LUAD patients. LUAD datasets were downloaded from public databases. Differential expression analysis was performed to screen DEGs, and prognostic modules were constructed by Cox regression. Enrichment analysis, gene regulatory network analysis and immune microenvironment analysis were employed to explore the underlying mechanisms in LUAD. Finally, the expression of prognosis module genes (PMGs) was validated in 8 LUAD tissue specimens and 5 cell lines by qRT-PCR. In total, 13 differential module genes (BIRC3, CCT6A, CPLX2, GLCCI1, GSTA1, HLA-DQB2, ID1, KCTD12, MUC15, NOTUM, NT5E, TCIM, and TM4SF4) were screened for the construction of a prognostic module. Notably, CCT6A and KCTD12 demonstrated excellent accuracy in the diagnosis of LUAD. Immune dysregulation and BIRC3, HLA-DQB2, KCTD12, and NT5E expression were significantly associated with invasive immune cells in LUAD patients. The expression level of CCT6A was highest in PC9-OR and H1975-OR cells, while the expression level of KCTD12 was highest in paracancerous tissue and HBE cells. The constructed prognostic model showed promise in predicting the survival of LUAD patients. Notably, KCTD12 and CCT6A might be candidate biomarkers for improving diagnostic performance and guiding individualized therapy for EGFR-TKI-resistant LUAD patients.

A differentially-methylated-region signature predicts the recurrence risk for patients with early stage lung adenocarcinoma

Predicting prognosis in lung cancer patients is important in establishing future treatment and monitoring plans. Lung adenocarcinoma (LUAD) is the most common and aggressive type of lung cancer with dismal prognosis and prognostic stratification would help to guide treatment. Aberrant DNA methylation in tumors occurs earlier than clinical variations, and keeps accumulating as cancer progresses. Preliminary studies have given us some clues that DNA methylation might serve as a promising biomarker for prognosis prediction. Herein, we aimed to study the potential utility of DNA methylation pattern in predicting the recurrence risk of early stage resectable LUAD and to develop a risk-modeling signature based on differentially methylated regions (DMRs). This study consisted of three cohorts of 244 patients with stage I-IIIA LUAD, including marker discovery cohort (n = 39), prognostic model training cohort (n = 117) and validation cohort (n = 80). 468 DMRs between LUAD tumor and adjacent tissues were screened out in the marker discovery cohort (adjusted P < 0.05), and a prognostic signature was developed based on 15 DMRs significantly related to disease-free survival in early stage LUAD patients. The DMR signature showed commendable performance in predicting the recurrence risk of LUAD patients both in model training cohort (P < 0.001; HR = 4.32, 95% CI = 2.39-7.80) and model validation cohort (P = 0.009; HR = 9.08, 95% CI = 1.20-68.80), which might be of great utility both for understanding the molecular basis of LUAD relapse, providing risk stratification of patients, and establishing future monitoring plans.

Development of a genome atlas for discriminating benign, preinvasive, and invasive lung nodules

To tackle misdiagnosis in lung cancer screening with low-dose computed tomography (LDCT), we aimed to compile a genome atlas for differentiating benign, preinvasive, and invasive lung nodules and characterize their molecular pathogenesis. We collected 432 lung nodule tissue samples from Chinese patients, spanning benign, atypical adenomatous hyperplasia (AAH), adenocarcinoma in situ (AIS), minimally invasive adenocarcinoma (MIA), and invasive adenocarcinoma (IA). We performed comprehensive sequencing, examining somatic variants, gene expressions, and methylation levels. Our findings uncovered EGFR and TP53 mutations as key drivers in - early lung cancer development, with EGFR mutation frequency increasing with disease progression. Both EGFR mutations and EGF/EGFR hypo-methylation activated the EGFR pathway, fueling cancer growth. Transcriptome analysis identified four lung nodule subtypes (G1-4) with distinct molecular features and immune cell infiltrations: EGFR-driven G1, EGFR/TP53 co-mutation G2, inflamed G3, stem-like G4. Estrogen/androgen response was associated with the EGFR pathway, proposing a new therapy combining tyrosine kinase inhibitors with antiestrogens. Preinvasive nodules exhibited stem cell pathway enrichment, potentially hindering invasion. Epigenetic regulation of various genes was essential for lung cancer initiation and development. This study provides insights into the molecular mechanism of neoplastic progression and identifies potential diagnostic biomarkers and therapeutic targets for lung cancer.

Pan-cancer prognosis, immune infiltration, and drug resistance characterization of lung squamous cell carcinoma tumor microenvironment-related genes

Background

The tumor microenvironment (TME) plays an important role in cancer development; however, its implications in lung squamous cell carcinoma (LUSC) and pan-cancer have been poorly understood.

Methods

In this study, The Cancer Genome Atlas (TCGA) and Estimation of Stromal and Immune cells in Malignant Tumor tissue using Expression Data (ESTIMATE) datasets were applied to identify differentially expressed genes. Additionally, online public databases were utilized for in-depth bioinformatics analysis of pan-cancer datasets to investigate the prognostic implications of TME-related genes further.

Results

Our study demonstrated a significant association between stromal scores, immune scores, and specific clinical characteristics in LUSC patients. C3AR1, CSF1R, CCL2, CCR1, and CD14 were identified as prognostic genes related to the TME. All TME-related prognostic genes demonstrated varying degrees of correlation with immune infiltration subtypes and tumor cell stemness. Moreover, our study revealed that TME-related prognostic genes, particularly C3AR1 and CCR1, might contribute to drug resistance in cancer cells.

Conclusions

The identified TME-related prognostic genes, particularly C3AR1 and CCR1, have potential implications for understanding and targeting drug resistance mechanisms in cancer cells.

TERC promotes non-small cell lung cancer progression by facilitating the nuclear localization of TERT

The core of telomerase consists of the protein subunit telomerase reverse transcriptase (TERT) and the telomerase RNA component (TERC). So far, the role of TERC in cancer development has remained elusive. Here, we found TERC expression elevated in non-small cell lung cancer (NSCLC) tissues, which was associated with disease progression and poor prognosis in patients. Using NSCLC cell lines and xenograft models, we showed that knockdown of TERC caused cell cycle arrest, and inhibition of cell proliferation and migration. Mechanistically, TERC was exported to the cytoplasm by nuclear RNA export factor 1 (NXF1), where it mediated the interaction of TERT with other telomerase subunits. Depletion of TERC hindered the assembly and subsequent nuclear localization of the telomerase complex, preventing TERT from functioning in telomere maintenance and transcription regulation. Our findings suggest that TERC is a potential biomarker for NSCLC diagnosis and prognosis and can be a target for NSCLC treatment.

Knockdown of NCAPD3 inhibits the tumorigenesis of non-small cell lung cancer by regulation of the PI3K/Akt pathway

BACKGROUND

Accumulating evidence indicates that aberrant non-SMC condensin II complex subunit D3 (NCAPD3) is associated with carcinogenesis of various cancers. Nevertheless, the biological role of NCAPD3 in the pathogenesis of non-small cell lung cancer (NSCLC) remains unclear.

METHODS

Immunohistochemistry and Western blot were performed to assess NCAPD3 expression in NSCLC tissues and cell lines. The ability of cell proliferation, invasion, and migration was evaluated by CCK-8 assays, EdU assays, Transwell assays, and scratch wound healing assays. Flow cytometry was performed to verify the cell cycle and apoptosis. RNA-sequence and rescue experiment were performed to reveal the underlying mechanisms.

RESULTS

The results showed that the expression of NCAPD3 was significantly elevated in NSCLC tissues. High NCAPD3 expression in NSCLC patients was substantially associated with a worse prognosis. Functionally, knockdown of NCAPD3 resulted in cell apoptosis and cell cycle arrest in NSCLC cells as well as a significant inhibition of proliferation, invasion, and migration. Furthermore, RNA-sequencing analysis suggested that NCAPD3 contributes to NSCLC carcinogenesis by regulating PI3K/Akt/FOXO4 pathway. Insulin-like growth factors-1 (IGF-1), an activator of PI3K/Akt signaling pathway, could reverse NCAPD3 silence-mediated proliferation inhibition and apoptosis in NSCLC cells.

CONCLUSION

NCAPD3 suppresses apoptosis and promotes cell proliferation via the PI3K/Akt/FOXO4 signaling pathway, suggesting a potential use for NCAPD3 inhibitors as NSCLC therapeutics.

Prognostic Model for High-Grade Neuroendocrine Carcinoma of the Lung Incorporating Genomic Profiling and Poly (ADP-ribose) Polymerase-1 Expression

PURPOSE

High-grade neuroendocrine carcinoma (HGNEC) of the lung is an aggressive cancer with a complex biology. We aimed to explore the prognostic value of genetic aberrations and poly(ADP-ribose) polymerase-1 (PARP1) expression in HGNEC and to establish a novel prognostic model.

MATERIALS AND METHODS

We retrospectively enrolled 191 patients with histologically confirmed HGNEC of the lung. Tumor tissues were analyzed using PARP1 immunohistochemistry (IHC; N = 191) and comprehensive cancer panel sequencing (n = 102). Clinical and genetic data were used to develop an integrated Cox hazards model.

RESULTS

Strong PARP1 IHC expression (intensity 3) was observed in 153 of 191 (80.1%) patients, and the mean PARP1 H-score was 285 (range, 5-300). To develop an integrated Cox hazard model, our data set included information from 357 gene mutations and 19 clinical profiles. When the targeted mutation profiles were combined with clinical profiles, 12 genes (ATRX, CCND2, EXT2, FGFR2, FOXO1, IL21R, MAF, TGM7, TNFAIP3, TP53, TSHR, and DDR2) were identified as prognostic factors for survival. The integrated Cox hazard model, which combines mutation profiles with a baseline model, outperformed the baseline model (incremental area under the curve 0.84 v 0.78; P = 8.79e-12). The integrated model stratified patients into high- and low-risk groups with significantly different disease-free and overall survival (integrated model: hazard ratio, 7.14 [95% CI, 4.07 to 12.54]; P < .01; baseline model: 4.38 [2.56 to 7.51]; P < .01).

CONCLUSION

We introduced a new prognostic model for HGNEC that combines genetic and clinical data. The integrated Cox hazard model outperformed the baseline model in predicting the survival of patients with HGNEC.

Mapping cancer biology in space: applications and perspectives on spatial omics for oncology

Technologies to decipher cellular biology, such as bulk sequencing technologies and single-cell sequencing technologies, have greatly assisted novel findings in tumor biology. Recent findings in tumor biology suggest that tumors construct architectures that influence the underlying cancerous mechanisms. Increasing research has reported novel techniques to map the tissue in a spatial context or targeted sampling-based characterization and has introduced such technologies to solve oncology regarding tumor heterogeneity, tumor microenvironment, and spatially located biomarkers. In this study, we address spatial technologies that can delineate the omics profile in a spatial context, novel findings discovered via spatial technologies in oncology, and suggest perspectives regarding therapeutic approaches and further technological developments.

UBE2C is a diagnosis and therapeutic biomarker involved in immune infiltration of cancers including lung adenocarcinoma

The mechanism of action of UBE2C in lung adenocarcinoma (LUAD) and its significance in cancer diagnosis, targeted therapy and immunotherapy, even in pan-cancer, are still unclear. Several large public databases and online analysis tools were used for big data mining analysis. RNA interference technology, CCK8 assay, flow cytometry and apoptosis detection, and western blot were used for in vitro experiments. UBE2C was found to be overexpressed in various of tumors, including LUAD, and its expression level was found to be significantly related to gender, weight, tumor stage, grade and prognosis in LUAD. Downregulation of UBE2C expression induced proliferation suppression and G2/M phase arrest and cell apoptosis in LUAD cells and suppressed LUAD cell growth through inhibiting the Akt-mTOR signaling pathway. Expression level of UBE2C was negatively correlated with B cells and CD4+ T cell, and also with immune checkpoint genes in LUAD. Pan-cancer assay shown that UBE2C was significantly overexpressed in 28 cancers and was correlated with Ki-67 index in many cancers. Overexpression of UBE2C in BRCA, LUAD and MESO indicated worse Overall Survival (OS). UBE2C expression levels were positively associated with immunocyte infiltration, immune regulatory genes, immune checkpoints, TMB, MSI and MMRs in some cancers. Additionally, Single-cell functional analysis showed that UBE2C was positively correlated with cell cycle, proliferation, DNA damage, EMT, DNA repair, invasion and differentiation in some cancers. These findings suggested that UBE2C could be regarded as a latent diagnosis and prognostic biomarker and a new target for immunological therapy of cancers including LUAD.

Carcinogenic Role and Clinical Significance of Histone H3-H4 Chaperone Anti-silencing Function 1 B (ASF1B) in Lung Adenocarcinoma

Histone H3-H4 chaperone anti-silencing function 1 (ASF1) plays an important role in the polymerization, transport, and modification of histones. However, the significance of ASF1B in lung adenocarcinoma (LUAD) is largely overlooked. We investigated the aberrant expression of ASF1B in LUAD and its potential link to patient survival using multiple databases. ASF1B-overexpressing and knockdown cell lines were constructed to explore its effects on the biological behavior of lung cancer cells. ssGSEA, TMB, TIDE and IMvigor210 cohort were used to explore and validate the association of ASF1B to tumor immunity. Our data suggested that ASF1B was overexpressed in LUAD, and was associated with poor prognosis. ASF1B promoted the proliferation, migration, and invasion of lung cancer cells by regulating the phosphorylation of AKT in vitro. ASF1B was associated with tumor immunity. In summary, ASF1B may promote malignant behavior of LUAD cells, and its overexpression correlates with worse prognosis and better immunotherapy effect.

Comprehensive analysis of mitophagy-related genes in NSCLC diagnosis and immune scenery: based on bulk and single-cell RNA sequencing data

Lung cancer is the main cause of cancer-related deaths, and non-small cell lung cancer (NSCLC) is the most common type. Understanding the potential mechanisms, prognosis, and treatment aspects of NSCLC is essential. This study systematically analyzed the correlation between mitophagy and NSCLC. Six mitophagy-related feature genes (SRC, UBB, PINK1, FUNDC1, MAP1LC3B, and CSNK2A1) were selected through machine learning and used to construct a diagnostic model for NSCLC. These feature genes are closely associated with the occurrence and development of NSCLC. Additionally, NSCLC was divided into two subtypes using unsupervised consensus clustering, and their differences in clinical characteristics, immune infiltration, and immunotherapy were systematically analyzed. Furthermore, the interaction between mitophagy-related genes (MRGs) and immune cells was analyzed using single-cell sequencing data. The findings of this study will contribute to the development of potential diagnostic biomarkers for NSCLC and the advancement of personalized treatment strategies.

FOXM1 and CENPF are associated with a poor prognosis through promoting proliferation and migration in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is a clinically challenging disease due to its poor prognosis and limited therapeutic methods. The aim of the present study was to identify prognosis-related genes and therapeutic targets for LUAD. Raw data from the GSE32863, GSE41271 and GSE42127 datasets were downloaded from the Gene Expression Omnibus database. Following normalization, the data were merged into a matrix, which was first used to identify differentially expressed genes (DEGs). Weighted gene co-expression network analysis (WGCNA) and survival analysis were performed to screen potential prognosis-related genes. Gene overlaps among DEGs, survival-related genes and WGCNA genes were finally constructed to obtain candidate genes. An analysis with the STRING database was performed to construct a protein-protein interaction network and hub genes were selected using Cytoscape. The candidate genes were finally identified by univariate and multivariate Cox regression analysis. Furthermore, in vivo and in vitro experiments, including immunohistochemistry, immunofluorescence, Cell Counting Kit-8, colony-formation and migration assays, were performed to validate the potential mechanism of these genes in LUAD. Two genes, namely forkhead box M1 (FOXM1) and centromere protein F (CENPF), were identified as unfavorable indicators of prognosis in patients with LUAD. High expression of FOXM1 and CENPF were associated with poor survival. Furthermore, LUAD cells with FOXM1 and CENPF knockdown showed a significant reduction in proliferation and migration (P<0.05). FOXM1 and CENPF may have an essential role in the prognosis of patients with LUAD by influencing cell proliferation and migration, and they provide potential molecular targets for LUAD therapy.

Deep neural network for discovering metabolism-related biomarkers for lung adenocarcinoma

Introduction

Lung cancer is a major cause of illness and death worldwide. Lung adenocarcinoma (LUAD) is its most common subtype. Metabolite-mRNA interactions play a crucial role in cancer metabolism. Thus, metabolism-related mRNAs are potential targets for cancer therapy.

Methods

This study constructed a network of metabolite-mRNA interactions (MMIs) using four databases. We retrieved mRNAs from the Tumor Genome Atlas (TCGA)-LUAD cohort showing significant expressional changes between tumor and non-tumor tissues and identified metabolism-related differential expression (DE) mRNAs among the MMIs. Candidate mRNAs showing significant contributions to the deep neural network (DNN) model were mined. Using MMIs and the results of function analysis, we created a subnetwork comprising candidate mRNAs and metabolites.

Results

Finally, 10 biomarkers were obtained after survival analysis and validation. Their good prognostic value in LUAD was validated in independent datasets. Their effectiveness was confirmed in the TCGA and an independent Clinical Proteomic Tumor Analysis Consortium (CPTAC) dataset by comparison with traditional machine-learning models.

Conclusion

To summarize, 10 metabolism-related biomarkers were identified, and their prognostic value was confirmed successfully through the MMI network and the DNN model. Our strategy bears implications to pave the way for investigating metabolic biomarkers in other cancers.

Anchoring Filament Protein Ladinin-1 is an Immunosuppressive Microenvironment and Cold Tumor Correlated Prognosticator in Lung Adenocarcinoma

Anchoring filament protein ladinin-1 (LAD1) codes for an anchor filament protein in the basement membrane. Here, we have aimed to determine its potential role in LUAD. According to the comprehensive analyses conducted in this study, we studied the expression, prognostic significance, function, methylation, copy number variations, and the immune cell infiltration of LAD1 in LUAD. A higher level of LAD1 gene expression was observed in the LUAD tumor tissues compared to the normal lung tissues (p < 0.001). Furthermore, the multivariate analysis indicated that a higher LAD1 gene expression level was the independent prognostic factor. Additionally, the DNA methylation level of the LAD1 was inversely linked to its expression (p < 0.001). We noted that the patients affected due to LAD1 hypomethylation showed a very low overall survival rate compared to the patients with a higher LAD1 methylation score (p < 0.05). Moreover, the results of the immunity analysis indicated that the LAD1 expression might be inversely linked to the immune cell infiltration degree, expression of the infiltrated immune cells, and the PD-L1 levels. Lastly, we supplemented some verification to increase the rigor of the study. The results suggested that high expression of LAD1 may be related to cold tumors. Hence, this indirectly reflects that the immunotherapy effect of LUAD patients with high LAD1 expression might be worse. Based on the role played by the LAD1 in the tumor immune microenvironment, it can be considered a potential biomarker for predicting the immunotherapy response to LUAD.

Identification of TLRs as potential prognostic biomarkers for lung adenocarcinoma

Lung adenocarcinoma (LUAD) is one of the most common tumors with the highest cancer-related death rate worldwide. Early diagnosis of LUAD can improve survival. Abnormal expression of the Toll-like receptors (TLRs) is related to tumorigenesis and development, inflammation and immune infiltration. However, the role of TLRs as an immunotherapy target and prognostic marker in lung adenocarcinoma is not well understood and needs to be analyzed. Relevant data obtained from databases such as ONCOMINE, UALCAN, GEPIA, and the Kaplan-Meier plotter, GSCALite, GeneMANIA, DAVID 6.8, Metascape, LinkedOmics and TIMER, to compare transcriptional TLRs and survival data of patients with LUAD. The expression levels of TLR1/2/3/4/5/7/8 in LUAD tissues were significantly reduced while the expression levels of TLR6/9/10 were significantly elevated. LUAD patients having low expression of TLR1/2/3/5/8 and high expression of TLR9 had a poor overall survival while patients with low expression of TLR2/3/7 presented with worse first progress. TLR4, TLR7 and TLR8 are the 3 most frequently mutated genes in the TLR family. Correlation suggested a low to moderate correlation among TLR family. TLR family was also involved in the activation or inhibition of the famous cancer related pathways. Analysis of immune infiltrates analysis suggested that TLR1/2/7/8 levels significantly correlated with immune infiltration level. Enrichment analysis revealed that TLRs were involved in TLR signaling pathway, immune response, inflammatory response, primary immunodeficiency, regulation of IL-8 production and PI3K-Akt signaling pathway. Our results provided information on TLRs expression and potential regulatory networks in LUAD. Moreover, our results suggested TLR2/7/8 as a potential prognostic biomarker for lung adenocarcinoma.

Immune Infiltration in Tumor and Adjacent Non-Neoplastic Regions Codetermines Patient Clinical Outcomes in Early-Stage Lung Cancer

INTRODUCTION

In recent years, the proportion of patients with NSCLC diagnosed at an early stage has increased continuously.

METHODS

In this study, we analyzed samples and data collected from 119 samples from 67 early stage patients with NSCLC, including 52 pairs of tumor and adjacent non-neoplastic samples, and performed RNA-sequencing analysis with high sequencing depth.

RESULTS

We found that immune-related genes were highly enriched among the differentially expressed genes and observed significantly higher inferred immune infiltration levels in adjacent non-neoplastic samples than in tumor samples. In survival analysis, the infiltration of certain immune cell types in tumor, but not adjacent non-neoplastic, samples were associated with overall patient survival, and excitingly, the differential infiltration between paired samples (tumor minus non-neoplastic) was more prognostic than expression in either non-neoplastic or tumor tissues. We also performed B cell receptor (BCR) and T cell receptor (TCR) repertoire analysis and observed more BCR/TCR clonotypes and increased BCR clonality in tumor than in non-neoplastic samples. Finally, we carefully quantified the fraction of the five histologic subtypes in our adenocarcinoma samples and found that higher histologic pattern complexity was associated with higher immune infiltration and low TCR clonality in the tumor-proximal regions.

CONCLUSIONS

Our results indicated significantly differential immune characteristics between tumor and adjacent non-neoplastic samples and suggested that the two regions provided complementary prognostic values in early-stage NSCLCs.

Exploring the methylation status of CFTR and PKIA genes as potential biomarkers for lung adenocarcinoma

BACKGROUND

One of the most prevalent cancers in the world is lung cancer, with adenocarcinoma (LUAD) making up a significant portion of cases. According to the National Cancer Institute (NCI), there are new cases and fatality rates per 100,000 individuals as follows: New instances of lung and bronchial cancer occur annually at a rate of 50.0 per 100,000 persons. The yearly death rate for men and women is 35.0 per 100,000. DNA methylation is one of the earliest discovered and widely studied epigenetic regulatory mechanisms, and its abnormality is closely related to the occurrence and development of cancer. However, the prognostic value of DNA methylation and LUAD needs to be further explored to improve the survival prediction of LUAD patients.

METHODS

The transcriptome data and clinical data of LUAD were downloaded from TCGA and GEO databases, and the Illumina Human Methylation450 array (450k array) data were downloaded from the TCGA database. Firstly, the intersection of the expressed genes of the two databases is corrected, the differential analysis is performed, and the methylation data is evaluated by the MethylMix package to obtain differentially methylated genes. Independent prognostic genes were screened out using univariate and multivariate Cox regression analysis, and a methylation prognostic model was developed using univariate Cox analysis and validated with the GSE30219 dataset in the GEO database. Survival analysis between methylation high-risk and low-risk groups was performed and a methylation-based gene prognostic model was constructed. Finally, the prediction of potential drugs associated with the LUAD gene signature using Drug Sensitivity Genomics in Cancer (GDSC).

RESULTS

In this study, a total of 555 samples from the TCGA database and 307 samples from GSE30219 were included, and a total of 24 differential methylation driver genes were identified. Univariate and multivariate Cox regression analyzes were used to screen out independent prognostic genes, involving 2 genes: CFTR, PKIA. Survival analysis was different between the methylation high-risk group and the low-risk group, the CFTR high methylation group and the low methylation group were poor, and the opposite was true for PKIA.

CONCLUSIONS

Our study revealed that the methylation status of CFTR and PKIA can serve as potential prognostic biomarkers and therapeutic targets in lung cancer.

Systematic analysis of histone acetylation regulators across human cancers

BACKGROUND

Histone acetylation (HA) is an important and common epigenetic pathway, which could be hijacked by tumor cells during carcinogenesis and cancer progression. However, the important role of HA across human cancers remains elusive.

METHODS

In this study, we performed a comprehensive analysis at multiple levels, aiming to systematically describe the molecular characteristics and clinical relevance of HA regulators in more than 10000 tumor samples representing 33 cancer types.

RESULTS

We found a highly heterogeneous genetic alteration landscape of HA regulators across different human cancer types. CNV alteration may be one of the major mechanisms leading to the expression perturbations in HA regulators. Furthermore, expression perturbations of HA regulators correlated with the activity of multiple hallmark oncogenic pathways. HA regulators were found to be potentially useful for the prognostic stratification of kidney renal clear cell carcinoma (KIRC). Additionally, we identified HDAC3 as a potential oncogene in lung adenocarcinoma (LUAD).

CONCLUSION

Overall, our results highlights the importance of HA regulators in cancer development, which may contribute to the development of clinical strategies for cancer treatment.

Comprehensive analysis of suppressor of cytokine signaling 2 protein in the malignant transformation of NSCLC

Suppressor of cytokine signaling 2 (SOCS2) plays an essential role in a number of physiological phenomena and functions as a tumor suppressor. Understanding the predictive effects of SOCS2 on non-small cell lung cancer (NSCLC) is urgently needed. The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) databases were used to assess SOCS2 gene expression levels in NSCLC. The clinical significance of SOCS2 was evaluated through Kaplan-Meier curve analysis and the analysis of related clinical factors. Gene Set Enrichment Analysis (GSEA) was used to identify the biological functions of SOCS2. Subsequently proliferation, wound-healing, colony formation and Transwell assays, and carboplatin drug experiments were used for verification. The results revealed that SOCS2 expression was low in the NSCLC tissues of patients in TCGA and GEO database analyses. Downregulated SOCS2 was associated with poor prognosis, as determined by Kaplan-Meier survival analysis (HR 0.61, 95% CI 0.52-0.73; P<0.001). GSEA showed that SOCS2 was involved in intracellular reactions, including epithelial-mesenchymal transition (EMT). Cell experiments indicated that knockdown of SOCS2 caused the malignant progression of NSCLC cell lines. Furthermore, the drug experiment showed that silencing of SOCS2 promoted the resistance of NSCLC cells to carboplatin. In conclusion, low expression of SOCS2 was associated with poor clinical prognosis by effecting EMT and causing drug resistance in NSCLC cell lines. Furthermore, SOCS2 could act as a predictive indicator for NSCLC.

The Prognostic Significance of KIAA1522 Expression in Non-Small-Cell Lung Cancer Patients

The majority of lung cancers belong to the non-small-cell lung cancer (NSCLC) category, which is linked to a high mortality rate despite significant progress in diagnosis and treatment. Therefore, there is a need for novel prognostic NSCLC biomarkers to improve prognosis which currently remains poor. Recent studies and analyses of gene expression data of NSCLC revealed that high expression of KIAA1522 was significantly associated with poor prognosis and decreased overall survival. We identified 98 patients who underwent radical curative surgical resections or metastasectomy for pulmonary adenocarcinoma and squamous cell carcinoma at our institution or the pathological diagnosis confirmed by our pathologists. Following the latest data, we utilized immunohistochemistry to assess the expression of KIAA1522 and investigated its association with various clinic-demographic parameters, pathological stages, recurrence rates, overall survival, and disease-free survival in patients who achieved complete remission. Notably, there were no significant differences in the expression profiles of KIAA1522 between adenocarcinoma and squamous cell carcinoma samples (p=0.6). Survival analysis was conducted using log-rank tests and a multivariate Cox proportional hazard model. Of the 98 samples, 54 (55.1%) exhibited high expression of KIAA1522, and patients with high KIAA1522 expression had a significantly shorter overall survival than the low-expression group (p=0.01). Multivariate Cox proportional hazard models in which metastatic patients were included revealed that along with older age, higher TNM stage (tumor, node, metastasis system), and Eastern Cooperative Oncology Group (ECOG) performance status, high expression of KIAA1522 served as an independent prognostic factor. A high expression profile was not significantly associated with relapses in those whose complete remission had been achieved. Still, those patients with high expression of KIAA1522 tended to exhibit a shorter disease-free survival rate. In conclusion, our findings suggest that KIAA1522 expression is an independent factor for predicting overall survival and may serve as a valuable prognostic indicator for relapse and disease-free survival in NSCLC patients.

Bioinformatics analysis and single-cell RNA sequencing: elucidating the ubiquitination pathways and key enzymes in lung adenocarcinoma

BACKGROUND

Lung adenocarcinoma (LUAD) is a prevalent subtype of lung cancer associated with high mortality rates. We aimed to utilize single-cell multiomics analysis to identify the key molecules involved in ubiquitination modification, which plays a role in LUAD development and progression.

METHODS

We use a systematic approach to analyze LUAD-related single-cell and bulk transcriptome datasets from Gene Expression Omnibus (GEO) and The Cancer Genome Atlas (TCGA) databases. Single-cell RNA sequencing (scRNA-seq) data were normalized, clustered, and annotated with the Seurat package in R. InferCNV was used to distinguish malignant from epithelial cells, and AUCell evaluated the area under the curve (AUC) score of ubiquitination-related enzymes. Survival and differential analyses identified significant molecular markers associated with ubiquitination. PSMD14 expression was confirmed using reverse-transcription quantitative polymerase chain reaction (RT-qPCR) and Western blot assays, and its knockdown cell lines were assessed for effects on cellular processes and tumor formation in mice. PSMD14's interacting proteins were predicted, and its impact on AGR2 protein half-life and ubiquitination was evaluated. Rescue experiments involving PSMD14 overexpression and AGR2 silencing assessed their impact on malignant behaviors.

RESULTS

By means of single-cell sequencing analysis, we probed the ubiquitination modification landscape in the LUAD microenvironment. Malignant cells had elevated scores for enzymes and ubiquitin-binding domains compared to normal epithelial cells, with 53 ubiquitination-related molecules showing prognostic disparities. FGR, PSMD14, and ZBTB16 were identified as genes with prognostic significance, with PSMD14 showing higher expression in epithelial and malignant cells. Two missense mutation sites were identified in PSMD14, which had a high copy number amplification ratio and positive correlation with messenger RNA (mRNA) expression. PSMD14 expression and tumor stage were found to be independent prognostic factors, and interfering with PSMD14 expression reduced the malignant behavior of LUAD cells. PSMD14 was found to bind to AGR2 protein and reduce its ubiquitination, leading to increased AGR2 stability. Knockdown of AGR2 inhibited the enhancement of cell viability, invasion, and migration resulting from PSMD14 overexpression.

CONCLUSIONS

This study examined ubiquitination modifications in LUAD using sequencing data, identifying PSMD14's critical role in malignancy regulation and its potential as a prognostic and therapeutic biomarker. These insights enhance understanding of LUAD mechanisms and treatment.

CAVIN2/SDPR Functioned as a Tumor Suppressor in Lung Adenocarcinoma from Systematic Analysis of Caveolae-Related Genes and Experimental Validation

Background: Caveolae-Related Genes include caveolins and cavins, which are the main component of the fossa and, play important roles in a variety of physiological and pathological processes. Although increasing evidence indicated that caveolins (CAVs) and cavins (CAVINs) are involved in carcinogenesis and progression, their clinical significance and biological function in lung cancer are still limited. Methods: We investigated the expression of CAVs and CAVINs at transcriptional levels using Oncomine and Gene Expression Profiling Interactive Analysis. The protein and mRNA expression levels of CAVs and CAVINs were determined by the human protein atlas website and our surgically resected samples, respectively. The clinical value of prognostic prediction based on the expression of CAVs and CAVINs was also assessed. cBioPortal, GeneMANIA and STRING were used to analyze the molecular characteristics of CAVs and CAVINs in lung adenocarcinoma (LUAD) comprehensively. Finally, we investigated the effect of CAVIN2/SDPR (serum deprivation protein response) on LUAD cells with biological experiments in vitro. Results: The expression of CAV1/2 and CAVIN1/2/3 were significantly downregulated in LUAD and lung squamous cell carcinoma (LUSC). The patients with high expression of CAV1, CAV2, CAV3, CAVIN1 and CAVIN2/SDPR were tightly correlated with a better prognosis in LUAD, while no statistical significances in LUSC. Further, our results found that CAVIN2/SDPR can be identified as a prognostic biomarker independent of other CAVINs in patients with LUAD. Mechanically, the overexpression of CAVIN2/SDPR inhibited cell proliferation and migration owing to the cell apoptosis induction and cell cycle arrest at S phase in LUAD cells. Conclusions: CAVIN2/SDPR functioned as a tumor suppressor, and was able to serve as prognostic biomarkers in precision medicine of LUAD. Mechanically, overexpression of CAVIN2/SDPR inhibited cell proliferation by inducing cell apoptosis and S phase arrest in LUAD cells.

MBD2 facilitates tumor metastasis by mitigating DDB2 expression

Despite past extensive studies, the pathoetiologies underlying tumor metastasis remain poorly understood, which renders its treatment largely unsuccessful. The methyl-CpG-binding domain 2 (MBD2), a "reader" to interpret DNA methylome-encoded information, has been noted to be involved in the development of certain types of tumors, while its exact impact on tumor metastasis remains elusive. Herein we demonstrated that patients with LUAD metastasis were highly correlated with enhanced MBD2 expression. Therefore, knockdown of MBD2 significantly attenuated the migration and invasion of LUAD cells (A549 and H1975 cell lines) coupled with attenuated epithelial-mesenchymal transition (EMT). Moreover, similar results were observed in other types of tumor cells (B16F10). Mechanistically, MBD2 selectively bound to the methylated CpG DNA within the DDB2 promoter, by which MBD2 repressed DDB2 expression to promote tumor metastasis. As a result, administration of MBD2 siRNA-loaded liposomes remarkably suppressed EMT along with attenuated tumor metastasis in the B16F10 tumor-bearing mice. Collectively, our study indicates that MBD2 could be a promising prognostic marker for tumor metastasis, while administration of MBD2 siRNA-loaded liposomes could be a viable therapeutic approach against tumor metastasis in clinical settings.

Association of frequent hypermethylation with high grade histological subtype in lung adenocarcinoma

Lung adenocarcinoma is classified morphologically into five histological subtypes according to the WHO classification. While each histological subtype correlates with a distinct prognosis, the molecular basis has not been fully elucidated. Here we conducted DNA methylation analysis of 30 lung adenocarcinoma cases annotated with the predominant histological subtypes and three normal lung cases using the Infinium BeadChip. Unsupervised hierarchical clustering analysis revealed three subgroups with different methylation levels: high-, intermediate-, and low-methylation epigenotypes (HME, IME, and LME). Micropapillary pattern (MPP)-predominant cases and those with MPP components were significantly enriched in HME (p = 0.02 and p = 0.03, respectively). HME cases showed a significantly poor prognosis for recurrence-free survival (p < 0.001) and overall survival (p = 0.006). We identified 365 HME marker genes specifically hypermethylated in HME cases with enrichment of "cell morphogenesis" related genes; 305 IME marker genes hypermethylated in HME and IME, but not in LME, with enrichment "embryonic organ morphogenesis"-related genes; 257 Common marker genes hypermethylated commonly in all cancer cases, with enrichment of "regionalization"-related genes. We extracted surrogate markers for each epigenotype and designed pyrosequencing primers for five HME markers (TCERG1L, CXCL12, FAM181B, HOXA11, GAD2), three IME markers (TBX18, ZNF154, NWD2) and three Common markers (SCT, GJD2, BARHL2). DNA methylation profiling using Infinium data was validated by pyrosequencing, and HME cases defined by pyrosequencing results also showed the worse recurrence-free survival. In conclusion, lung adenocarcinomas are stratified into subtypes with distinct DNA methylation levels, and the high-methylation subtype correlated with MPP-predominant cases and those with MPP components and showed a poor prognosis.

CircSCN8A suppresses malignant progression and induces ferroptosis in non-small cell lung cancer by regulating miR-1290/ACSL4 axis

Circular RNAs (CircRNAs) are reported to exert vital regulatory roles in the occurrence and development of various human malignancies, including non-small cell lung cancer (NSCLC). Bioinformatics methods identified the down-regulation of circSCN8A (circBase ID: hsa_circ_0026337) in NSCLC tissues. However, its biological functions and molecular mechanisms in NSCLC remain unknown. In this study, we found that circSCN8A expression was down-regulated in NSCLC tissues and cells. Low circSCN8A expression was positively associated with aggressive clinicopathological characteristics and poor prognosis in NSCLC patients. CircSCN8A suppressed cell proliferation, migration, invasion, and epithelial-mesenchymal transition (EMT) in vitro and blocked tumor growth in vivo. Moreover, circSCN8A promoted cell ferroptosis in NSCLC. Mechanistically, circSCN8A acted as a competing endogenous RNA (ceRNA) by sponging miR-1290 to enhance the expression of long-chain acyl-CoA synthetase-4 (ACSL4). Furthermore, the knockdown of ACSL4 or overexpression of miR-1290 reversed the effect of circSCN8A on facilitating ferroptosis and inhibiting cell proliferation and metastasis. In summary, circSCN8A represses cell proliferation and metastasis in NSCLC by regulating the miR-1290/ACSL4 axis to induce ferroptosis. Thus, circSCN8A may represent a promising therapeutic target against NSCLC.

High Expression of COA6 Is Related to Unfavorable Prognosis and Enhanced Oxidative Phosphorylation in Lung Adenocarcinoma

Mitochondrial metabolism plays an important role in the occurrence and development of cancers. Cytochrome C oxidase assembly factor six (COA6) is essential in mitochondrial metabolism. However, the role of COA6 in lung adenocarcinoma (LUAD) remains unknown. Here we report that the expression of COA6 mRNA and protein were upregulated in LUAD tissues compared with lung normal tissues. We found that COA6 had high sensitivity and specificity to distinguish LUAD tissues from normal lung tissues shown by a receiver operating characteristic (ROC) curve. In addition, our univariate and multivariate Cox regression analysis indicated that COA6 was an independent unfavorable prognostic factor for LUAD patients. Furthermore, our survival analysis and nomogram showed that a high expression of COA6 mRNA was related to the short overall survival (OS) of LUAD patients. Notably, our weighted correlation network analysis (WGCNA) and functional enrichment analysis revealed that COA6 may participate in the development of LUAD by affecting mitochondrial oxidative phosphorylation (OXPHOS). Importantly, we demonstrated that depletion of COA6 could decrease the mitochondrial membrane potential (MMP), nicotinamide adenine dinucleotide (NAD) + hydrogen (H) (NADH), and adenosine triphosphate (ATP) production in LUAD cells (A549 and H1975), hence inhibiting the proliferation of these cells in vitro. Together, our study strongly suggests that COA6 is significantly associated with the prognosis and OXPHOS in LUAD. Hence, COA6 is highly likely a novel prognostic biomarker and therapeutic target of LUAD.

Core genes in lung adenocarcinoma identified by integrated bioinformatic analysis

This study aims to identify potential core genes of lung adenocarcinoma (LUAD). Three datasets (GSE32863, GSE43458, and GSE116959) were retrieved from the Gene Expression Omnibus (GEO) database. Differentially expressed genes (DEGs) between LUAD and normal tissues were filtrated by GEO2R tool. Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway analysis were performed via Metascape database. The protein-protein interaction (PPI) network was constructed and core genes were identified using STRING and Cytoscape. Core genes expressions and their relevant clinical characteristics were performed via Oncomine and UALCAN databases respectively. The correlation between core genes and immune infiltrates was investigated by TIMER database. Kaplan-Meier plotter was performed for survival analysis. The signal pathway network of core genes was mapped by KEGG Mapper analysis tool. In this study, ten core genes were significantly related to overall survival (OS) of LUAD patients, which can provide clues for prognosis of LUAD.

MTSS1 curtails lung adenocarcinoma immune evasion by promoting AIP4-mediated PD-L1 monoubiquitination and lysosomal degradation

Immune checkpoint blockade (ICB) therapy targeting PD-1/PD-L1 has shown durable clinical benefits in lung cancer. However, many patients respond poorly to ICB treatment, underscoring an incomplete understanding of PD-L1 regulation and therapy resistance. Here, we find that MTSS1 is downregulated in lung adenocarcinoma, leading to PD-L1 upregulation, impairment of CD8⁺ lymphocyte function, and enhanced tumor progression. MTSS1 downregulation correlates with improved ICB efficacy in patients. Mechanistically, MTSS1 interacts with the E3 ligase AIP4 for PD-L1 monoubiquitination at Lysine 263, leading to PD-L1 endocytic sorting and lysosomal degradation. In addition, EGFR-KRAS signaling in lung adenocarcinoma suppresses MTSS1 and upregulates PD-L1. More importantly, combining AIP4-targeting via the clinical antidepressant drug clomipramine and ICB treatment improves therapy response and effectively suppresses the growth of ICB-resistant tumors in immunocompetent mice and humanized mice. Overall, our study discovers an MTSS1-AIP4 axis for PD-L1 monoubiquitination and reveals a potential combinatory therapy with antidepressants and ICB.

When "No-Smoking" is not enough: Hypoxia and nicotine acetylcholine receptor signaling may drive lung adenocarcinoma progression in never-smokers

The question of how lung cancer progresses in never-smokers remains largely unanswered. In our analysis of data from 1727 lung cancer patients, we observed a difference of only 47 days in the overall survival between lung adenocarcinoma patients who were smokers vis-a-vis never-smokers - the disease has a poor prognosis irrespective of the smoking status, or gender. We have investigated the possible collaboration between the nAChR and hypoxia signaling pathway to explicate a mechanism of disease progression in never-smokers using patient-derived tumor cells. We found a previously unidentified increase in both acetylcholine and nAChR-α7 levels in non-small cell lung cancer cells in hypoxia. A similar increase in ubiquitously expressed nAChR-α7 transcripts was also observed in other cancer lines and primary tumor tissues. A direct binding of HIF-1α with the hypoxia-response element (HRE) present at -48 position preceding the transcriptional start site in nAChR-α7 promoter region was established. Crucially, the increased acetylcholine levels in hypoxia drove a feedback loop via modulation of PI3K/AKT pathway to stabilize HIF-1α in hypoxia. Further, hypoxia-mediated metastasis and induction of HIF-1α in these cells was significantly reversed by bungarotoxin, an antagonist of nAChR-α7. The nAChR-AKT-HIF network needs to be further investigated to conclusively prove its mechanism and to explore its therapeutic potential. Our study gives a plausible explanation for the equally worse prognosis of lung adenocarcinoma in never-smokers wherein the nAChR signaling is enhanced in hypoxia by acetylcholine in the absence of nicotine.

STRIP2 motivates non-small cell lung cancer progression by modulating the TMBIM6 stability through IGF2BP3 dependent

BACKGROUND

Striatin interacting protein 2 (STRIP2) is a core component of the striatin-interacting phosphatase and kinase (STRIPAK) complexes, which is involved in tumor initiation and progression via the regulation of cell contractile and metastasis. However, the underlying molecular mechanisms of STRIP2 in non-small cell lung cancer (NSCLC) progression remain largely unknown.

METHODS

The expressions of STRIP2 and IGF2BP3 in human NSCLC specimens and NSCLC cell lines were detected using quantitative RT-PCR, western blotting, and immunohistochemistry (IHC) analyses. The roles and molecular mechanisms of STRIP2 in promoting NSCLC progression were investigated in vitro and in vivo.

RESULTS

Here, we found that STRIP2 expression was significantly elevated in NSCLC tissues and high STRIP2 expression was associated with a poor prognosis. Knockdown of STRIP2 suppressed tumor growth and metastasis in vitro and in vivo, while STRIP2 overexpression obtained the opposite effect. Mechanistically, P300/CBP-mediated H3K27 acetylation activation in the promoter of STRIP2 induced STRIP2 transcription, which interacted with insulin-like growth factor 2 mRNA-binding protein 3 (IGF2BP3) and upregulated IGF2BP3 transcription. In addition, STRIP2-IGF2BP3 axis stimulated m6A modification of TMBIM6 mRNA and enhanced TMBIM6 stability. Consequently, TMBIM6 involved NSCLC cell proliferation, migration and invasion dependent on STRIP2 and IGF2BP3. In NSCLC patients, high co-expression of STRIP2, IGF2BP3 and TMBIM6 was associated with poor outcomes.

CONCLUSIONS

Our findings indicate that STRIP2 interacts with IGF2BP3 to regulate TMBIM6 mRNA stability in an m6A-dependent manner and may represent a potential prognostic biomarker and therapeutic target for NSCLC.

Xanthine dehydrogenase rewires metabolism and the survival of nutrient deprived lung adenocarcinoma cells by facilitating UPR and autophagic degradation

Xanthine dehydrogenase (XDH) is the rate-limiting enzyme in purine catabolism by converting hypoxanthine to xanthine and xanthine to uric acid. The altered expression and activity of XDH are associated with the development and prognosis of multiple types of cancer, while its role in lung adenocarcinoma (LUAD) remains unknown. Herein, we demonstrated that XDH was highly expressed in LUAD and was significantly correlated with poor prognosis. Though inhibition of XDH displayed moderate effect on the viability of LUAD cells cultured in the complete medium, it significantly attenuated the survival of starved cells. Similar results were obtained in XDH-knockout cells. Nucleosides supplementation rescued the survival of starved LUAD cells upon XDH inhibition, while inhibition of purine nucleoside phosphorylase abrogated the process, indicating that nucleoside degradation is required for the XDH-mediated survival of LUAD cells. Accordingly, metabolic flux revealed that ribose derived from nucleoside fueled key carbon metabolic pathways to sustain the survival of starved LUAD cells. Mechanistically, down-regulation of XDH suppressed unfolded protein response (UPR) and autophagic flux in starved LUAD cells. Inhibition of XDH decreased the level of amino acids produced by autophagic degradation, which was accompanied with down-regulation of mTORC1 signaling. Supplementation of amino acids including glutamine or glutamate rescued the survival of starved LUAD cells upon knockout or inhibition of XDH. Finally, XDH inhibitors potentiated the anti-cancer activity of 2-deoxy-D-glucose that induced UPR and/or autophagy in vitro and in vivo. In summary, XDH plays a crucial role in the survival of starved LUAD cells and targeting XDH may improve the efficacy of drugs that induce UPR and autophagy in the therapy of LUAD.

Penetrating Exploration of Prognostic Correlations of the FKBP Gene Family with Lung Adenocarcinoma

The complexity of lung adenocarcinoma (LUAD), the development of which involves many interacting biological processes, makes it difficult to find therapeutic biomarkers for treatment. FK506-binding proteins (FKBPs) are composed of 12 members classified as conservative intracellular immunophilin family proteins, which are often connected to cyclophilin structures by tetratricopeptide repeat domains and have peptidyl prolyl isomerase activity that catalyzes proline from residues and turns the trans form into the cis form. Since FKBPs belong to chaperone molecules and promote protein folding, previous studies demonstrated that FKBP family members significantly contribute to the degradation of damaged, misfolded, abnormal, and foreign proteins. However, transcript expressions of this gene family in LUAD still need to be more fully investigated. In this research, we adopted high-throughput bioinformatics technology to analyze FKBP family genes in LUAD to provide credible information to clinicians and promote the development of novel cancer target drugs in the future. The current data revealed that the messenger (m)RNA levels of FKBP2, FKBP3, FKBP4, FKBP10, FKBP11, and FKBP14 were overexpressed in LUAD, and FKBP10 had connections to poor prognoses among LUAD patients in an overall survival (OS) analysis. Based on the above results, we selected FKBP10 to further conduct a comprehensive analysis of the downstream pathway and network. Through a DAVID analysis, we found that FKBP10 was involved in mitochondrial electron transport, NADH to ubiquinone transport, mitochondrial respiratory chain complex I assembly, etc. The MetaCore pathway analysis also indicated that FKBP10 was involved in "Ubiquinone metabolism", "Translation_(L)-selenoaminoacid incorporation in proteins during translation", and "Transcription_Negative regulation of HIF1A function". Collectively, this study revealed that FKBP family members are both significant prognostic biomarkers for lung cancer progression and promising clinical therapeutic targets, thus providing new targets for treating LUAD patients.

An evaluation of KIF20A as a prognostic factor and therapeutic target for lung adenocarcinoma using integrated bioinformatics analysis

The identification of prognostic and therapeutic biomarkers is essential to reduce morbidity and mortality from lung adenocarcinoma (LUAD). This study aimed to identify a reliable prognostic and therapeutic biomarker for LUAD using integrated bioinformatics. Based on the cancer genome atlas (TCGA) and genome-tissue expression (GTEx) analyses, KIF20A has been identified as the hub gene. Following validation using a series of cohorts, survival analysis, meta-analysis, and univariate Cox analysis was conducted. ESTIMATE and CIBERSORT algorithms were then used to study the association of KIF20A with the tumor microenvironment (TME) and the percentage of tumor-infiltrating immune cells (TICs). In vitro experiments were conducted to determine the function of KIF20A. Finally, there was a negative association between the expression of the KIF20A and overall survival, progression-free survival, and disease-free survival, which was confirmed by meta-analysis and COX analysis. Furthermore, KIF20A also had a potential role of altering the TME and TICs proportions in LUAD. Validations in vitro were performed on A549 and PC-9 cell lines, and we found that the knockdown of KIF20A exhibited inhibitory effects on cell proliferation, resulted in cell cycle arrest during the G2/M phase, and induced cellular apoptosis. Our study demonstrated that KIF20A could be utilized as a reliable prognostic marker and treatment target for LUAD. However, further studies are required to validate these findings.

Comprehensive analysis of the significance of METTL7A gene in the prognosis of lung adenocarcinoma

Background

The most common subtype of lung cancer, called lung adenocarcinoma (LUAD), is also the largest cause of cancer death in the world. The aim of this study was to determine the importance of the METTL7A gene in the prognosis of patients with LUAD.

Methods

This particular study used a total of four different LUAD datasets, namely TCGA-LUAD, GSE32863, GSE31210 and GSE13213. Using RT-qPCR, we were able to determine METTL7A expression levels in clinical samples. Univariate and multivariate Cox regression analyses were used to identify factors with independent effects on prognosis in patients with LUAD, and nomograms were designed to predict survival in these patients. Using gene set variation analysis (GSVA), we investigated differences in enriched pathways between METTL7A high and low expression groups. Microenvironmental cell population counter (MCP-counter) and single-sample gene set enrichment analysis (ssGSEA) methods were used to study immune infiltration in LUAD samples. Using the ESTIMATE technique, we were able to determine the immune score, stromal score, and estimated score for each LUAD patient. A competing endogenous RNA network, also known as ceRNA, was established with the help of the Cytoscape program.

Results

We detected that METTL7A was down-regulated in pan-cancer, including LUAD. The survival study indicates that METTL7A was a protective factor in the prognosis of LUAD. The univariate and multivariate Cox regression analyses revealed that METTL7A was a robust independent prognostic indicator in survival prediction. Through the use of GSVA, several immune-related pathways were shown to be enriched in both the high-expression and low-expression groups of METTL7A. Analysis of the tumor microenvironment revealed that the immune microenvironment of the group with low expression was suppressed, which may be connected to the poor prognosis. To explore the ceRNA regulatory mechanism of METTL7A, we finally constructed a regulatory network containing 1 mRNA, 2 miRNAs, and 5 long non-coding RNAs (lncRNAs).

Conclusion

In conclusion, we presented METTL7A as a potential and promising prognostic indicator of LUAD. This biomarker has the potential to offer us with a comprehensive perspective of the prediction of prognosis and treatment for LUAD patients.

Glutathione peroxidase 2 overexpression promotes malignant progression and cisplatin resistance of KRAS‑mutated lung cancer cells

Kirsten rat sarcoma viral oncogene homolog (KRAS) aberrations frequently occur in patients with lung cancer. Oncogenic KRAS is characterized by excessive reactive oxygen species (ROS) accumulation, thus, ROS detoxification may contribute to KRAS‑driven lung tumorigenesis. In the present study, the influence of glutathione peroxidase 2 (GPX2) on malignant progression and cisplatin resistance of KRAS‑driven lung cancer was explored. The RNA sequencing data from TCGA lung cancer samples and GEO database were downloaded and analyzed. The effects of GPX2 on KRAS‑driven lung tumorigenesis were evaluated by western blotting, cell viability assay, soft agar assay, Transwell assay, tumor xenograft model, flow cytometry, BrdU incorporation assay, transcriptome RNA sequencing, luciferase reporter assay and RNA immunoprecipitation. In the present study, GPX2 was upregulated in patients with non‑small cell lung carcinoma (NSCLC), and positively correlated with poor overall survival. Ectopic GPX2 expression facilitated malignant progression of KRASG12C‑transformed BEAS‑2B cells. Moreover, GPX2 overexpression promoted growth, migration, invasion, tumor xenograft growth and cisplatin resistance of KRAS‑mutated NSCLC cells, while GPX2 knockdown exhibited the opposite effects. GPX2 overexpression reduced ROS accumulation and increased matrix metalloproteinase‑1 (MMP1) expression in KRAS‑mutated NSCLC cells. In addition, GPX2 was directly targeted by miR‑325‑3p, while MMP1 knockdown or miR‑325‑3p overexpression partially abrogated the effects of GPX2 in NSCLC cells. In conclusion, the results indicated that GPX2 facilitated malignant progression and cisplatin resistance of KRAS‑driven lung cancer, and inhibition of GPX2 may be a feasible strategy for lung cancer treatment, particularly in patients with active KRAS mutations.

The roles of CPSF6 in proliferation, apoptosis and tumorigenicity of lung adenocarcinoma

Cleavage and polyadenylation specific factor 6 (CPSF6), a member of serine/arginine-rich protein family, is implicated in HIV-1-infection and replication. Overexpression of CPSF6 predicts poor prognostic outcomes of breast cancer. However, the expression and possible function of CPSF6 in lung adenocarcinoma (LUAD) still needs to be explored. Here, we found that CPSF6 is significantly higher expressed in tumor tissues than normal tissues in multiple cancer types. Besides, CPSF6 plays a significant risky role in LUAD that is associated with overall survival (HR=1.337, P=0.051) and disease-specific survival (HR=1.4739, P=0.042). CPSF6 mRNA was up-regulated in LUAD tissues by analyzing publicly available datasets from Gene Expression Omnibus (GEO). Further survival analysis on The Cancer Genome Atlas (TCGA) dataset suggested a close correlation between CPSF6 expression and overall survival, and disease-free survival of LUAD patients. Inhibition of CPSF6 expression by lentivirus-mediated RNA interference (RNAi) in two LUAD cell lines (A549 and NCH-H1299) caused a significant reduction in cell proliferation, colony formation and a notable induction in apoptotic rate. CPSF6 knockdown in xenograft tumors inhibited LUAD cell growth in vivo. Moreover, we identified differentially expressed genes with CPSF6 inhibition by Microarray analysis, and pathway analyses revealed that CPSF6 knockdown resulted in the dysregulation of the phosphatidylinositol 3-kinase/protein kinase B (PI3K/AKT) pathway. Collectively, our results are the first to demonstrate that CPSF6 functions as an oncoprotein by regulating cancer-related pathways in LUAD.

Association of the Timeless Gene with Prognosis and Clinical Characteristics of Human Lung Cancer

(1) Background: As the most common malignant tumor type worldwide, it is necessary to identify novel potential prognostic biomarkers to improve the poor prognosis of lung cancer. The Timeless gene, a circadian rhythm-related gene, is associated with several types of cancer. However, studies analyzing the clinical significance of the Timeless gene in patients with lung cancer are currently limited. (2) Methods: In the present study, the expression levels and prognostic potential of the Timeless gene and its co-expressed genes in different subtypes of lung cancer were explored using multiple bioinformatics approaches. The correlations between the Timeless gene and its co-expressed genes were validated using A549 and NCI-H226 cells by transfecting them with expression vectors and analyses using Western blot and reverse transcription-quantitative PCR. (3) Results: The Oncomine and GEPIA database analyses indicated that the expression of the Timeless gene was significantly higher in lung cancer as compared to that in the normal tissue. Using the UALCAN database, significant differences in Timeless gene expression were determined among different stages of lung cancer and between genders. A Kaplan-Meier plotter analysis indicated that high expression of the Timeless gene was associated with poor overall survival (OS) and progression-free survival (PFS) of patients with lung cancer. In the cBioPortal and GEPIA database analyses, extra spindle pole bodies like 1 (ESPL1) was the top correlated gene of Timeless in patients with lung cancer. Similar to the Timeless gene, high expression of the ESPL1 gene was also associated with poor OS and PFS. Of note, overexpression of the Timeless gene increased the expression level of ESPL1 at both the mRNA and protein levels. (4) Conclusion: The present study explored the clinical significance of the Timeless gene and its correlated gene ESPL1 in patients with lung cancer, thereby providing a potential therapeutic target for lung cancer.

Analysis on methylation and expression of PSMB8 and its correlation with immunity and immunotherapy in lung adenocarcinoma

Aim: To find biomarkers for immunity and immunotherapy in lung adenocarcinoma (LUAD) through multiomics analysis. Materials & methods: The multiomics data of patients with LUAD were downloaded from the TCGA and GEO databases. CIBERSORT, quanTIseq, ESTIMATEScore, k-means clustering, gene set enrichment analysis, gene set variation analysis, immunophenoscore and logistic regression were used in this study. Results: PSMB8 HypoMet-HighExp group patients have more active immune-related pathways, more antitumor immune cells, less protumor immune cells, higher immunophenoscore and longer progression-free survival of immune checkpoint inhibitor therapy than HyperMet-LowExp group. In multivariate analysis, PSMB8 showed an independent value. Conclusion: The combination of DNA methylation and mRNA expression of PSMB8 could independently distinguish types of tumor immune microenvironment and predict programmed cell death protein 1/programmed cell death-ligand 1 inhibitors' effects in patients with LUAD.

Prognostic Value of GIMAP4 and Its Role in Promoting Immune Cell Infiltration into Tumor Microenvironment of Lung Adenocarcinoma

GIMAPs are recognized as an important regulator in the carcinogenesis and development of lung cancer, but the function of GIMAP4 in the tumor microenvironment (TME) of lung cancers is unclear. In this study, we investigated the expression and variation of GIMAP4 in lung adenocarcinoma (LUAD), to explore its association with infiltration of immune cells. The Cancer Genome Atlas (TCGA) data and Gene Expression Omnibus (GEO) data were analyzed. Infiltration of immune cells was identified with TIMER (Tumor Immune Estimation Resource) and TISIDB (an integrated repository portal for tumor-immune system interactions). GIMAP4 expression declined in non-small-cell lung cancer (NSCLC), correlated with a poor overall survival (OS) in LUAD, indicating that GIMAP4 was a promising prognostic biomarker in LUAD. GIMAP4 mutation frequency was 1.76% in TCGA cohort and was relevant to the expression of immune components. TIMER and CIBERSORT analysis further confirmed that high GIMAP4 expression possibly promoted immune cell infiltration into the TME, with low GIMAP4 impairing the efficacy of immunotherapies targeting common immune check point inhibitors (ICI). GO (Gene Ontology) and KEGG (Kyoto Encyclopedia of Genes and Genomes) analyses were performed to provide insights into biological processes involved in LUAD. GIMAP4 was expected to be a prognostic biomarker in LUAD and provides potential adjuvant or neoadjuvant therapeutic strategies for targeting ICIs.

Immune infiltration landscape on prognosis and therapeutic response and relevant epigenetic and transcriptomic mechanisms in lung adenocarcinoma

Objective

Lung adenocarcinoma (LUAD) is the most prevalent lung cancer subtype, but its immune infiltration features are not comprehensively understood. To address the issue, the present study was initiated to describe the immune infiltrations across LUAD from cellular compositional, functional, and mechanism perspectives.

Methods

We adopted five LUAD datasets (GSE32863, GSE43458, GSE75037, TCGA-LUAD, and GSE72094). Differentially expressed genes between LUAD and controls were selected for co-expression network analysis. Risky immune cell types were determined for classifying LUAD patients as diverse subtypes, followed by a comparison of antitumor immunity and therapeutic response between subtypes. Then, LUAD- and subtype-related key module genes affected by DNA methylation were determined for quantifying a scoring scheme. EXO1 was chosen for functional analysis via in vitro assays.

Results

Two immune cell infiltration-based subtypes (C1 and C2) were established across LUAD, with poorer prognostic outcomes and lower infiltration of immune cell types in C1. Additionally, C1 presented higher responses to immune checkpoint blockade and targeted agents (JNK inhibitor VIII, BI-D1870, RO-3306, etc.). The scoring system (comprising GAPDH, EXO1, FYN, CFTR, and KLF4) possessed higher accuracy in estimating patients’ prognostic outcomes. EXO1 upregulation contributed to the growth, migration, and invasion of LUAD cells. In addition, EXO1 facilitated PD-L1 and sPD-L1 expression in LUAD cells.

Conclusion

Altogether, our findings offer a comprehensive understanding of the immune infiltration landscape on prognosis and therapeutic response of LUAD as well as unveil potential epigenetic and transcriptomic mechanisms, which might assist personalized treatment.

Smoking-associated Downregulation of FILIP1L Enhances Lung Adenocarcinoma Progression Through Mucin Production, Inflammation, and Fibrosis

Lung adenocarcinoma (LUAD) is the major subtype in lung cancer, and cigarette smoking is essentially linked to its pathogenesis. We show that downregulation of Filamin A interacting protein 1-like (FILIP1L) is a driver of LUAD progression. Cigarette smoking causes its downregulation by promoter methylation in LUAD. Loss of FILIP1L increases xenograft growth, and, in lung-specific knockout mice, induces lung adenoma formation and mucin secretion. In syngeneic allograft tumors, reduction of FILIP1L and subsequent increase in its binding partner, prefoldin 1 (PFDN1) increases mucin secretion, proliferation, inflammation, and fibrosis. Importantly, from the RNA-sequencing analysis of these tumors, reduction of FILIP1L is associated with upregulated Wnt/β-catenin signaling, which has been implicated in proliferation of cancer cells as well as inflammation and fibrosis within the tumor microenvironment. Overall, these findings suggest that down-regulation of FILIP1L is clinically relevant in LUAD, and warrant further efforts to evaluate pharmacologic regimens that either directly or indirectly restore FILIP1L-mediated gene regulation for the treatment of these neoplasms.

Significance

This study identifies FILIP1L as a tumor suppressor in LUADs and demonstrates that downregulation of FILIP1L is a clinically relevant event in the pathogenesis and clinical course of these neoplasms.

Identification of HMGB2 associated with proliferation, invasion and prognosis in lung adenocarcinoma via weighted gene co-expression network analysis

Background

High mobility group protein B2 (HMGB2) is a multifunctional protein that plays various roles in different cellular compartments. Moreover, HMGB2 serves as a potential prognostic biomarker and therapeutic target for lung adenocarcinoma (LUAD).

Methods

In this study, the expression pattern, prognostic implication, and potential role of HMGB2 in LUAD were evaluated using the integrated bioinformatics analyses based on public available mRNA expression profiles from The Cancer Genome Atlas and Gene Expression Omnibus databases, both at the single-cell level and the tissue level. Further study in the patient-derived samples was conducted to explore the correlation between HMGB2 protein expression levels with tissue specificity, (tumor size-lymph node-metastasis) TNM stage, pathological grade, Ki-67 status, and overall survival. In vitro experiments, such as CCK-8, colony-formation and Transwell assay, were performed with human LUAD cell line A549 to investigate the role of HMGB2 in LUAD progression. Furthermore, xenograft tumor model was generated with A549 in nude mice.

Results

The results showed that the HMGB2 expression was higher in the LUAD samples than in the adjacent normal tissues and was correlated with high degree of malignancy in different public data in this study. Besides, over-expression of HMGB2 promoted A549 cells proliferation and migration while knocking down of HMGB2 suppressed the tumor promoting effect.

Conclusions

Our study indicated that HMGB2 was remarkably highly expressed in LUAD tissues, suggesting that it is a promising diagnostic and therapeutic marker for LUAD in the future.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12890-022-02110-y.

Epigenetic underpinnings of inflammation: Connecting the dots between pulmonary diseases, lung cancer and COVID-19

Inflammation is an essential component of several respiratory diseases, such as chronic obstructive pulmonary disease (COPD), asthma and acute respiratory distress syndrome (ARDS). It is central to lung cancer, the leading cancer in terms of associated mortality that has affected millions of individuals worldwide. Inflammation and pulmonary manifestations are also the major causes of COVID-19 related deaths. Acute hyperinflammation plays an important role in the COVID-19 disease progression and severity, and development of protective immunity against the virus is greatly sought. Further, the severity of COVID-19 is greatly enhanced in lung cancer patients, probably due to the genes such as ACE2, TMPRSS2, PAI-1 and furin that are commonly involved in cancer progression as well as SAR-CoV-2 infection. The importance of inflammation in pulmonary manifestations, cancer and COVID-19 calls for a closer look at the underlying processes, particularly the associated increase in IL-6 and other cytokines, the dysregulation of immune cells and the coagulation pathway. Towards this end, several reports have identified epigenetic regulation of inflammation at different levels. Expression of several key inflammation-related cytokines, chemokines and other genes is affected by methylation and acetylation while non-coding RNAs, including microRNAs as well as long non-coding RNAs, also affect the overall inflammatory responses. Select miRNAs can regulate inflammation in COVID-19 infection, lung cancer as well as other inflammatory lung diseases, and can serve as epigenetic links that can be therapeutically targeted. Furthermore, epigenetic changes also mediate the environmental factors-induced inflammation. Therefore, a better understanding of epigenetic regulation of inflammation can potentially help develop novel strategies to prevent, diagnose and treat chronic pulmonary diseases, lung cancer and COVID-19.

CCT6A and CHCHD2 Are Coamplified with EGFR and Associated with the Unfavorable Clinical Outcomes of Lung Adenocarcinoma

Chaperonin containing TCP1 subunit 6A (CCT6A) and coiled-coil-helix-coiled-coil-helix domain containing 2 (CHCHD2) are located at the chromosome 7p11 region proximal to epidermal growth factor receptor (EGFR). However, the amplifications, expressions, and the prognostic effects of CCT6A and CHCDH2 in lung adenocarcinoma (LUAD) are unclear. Here, using The Cancer Genome Atlas (TCGA) and Gene Expression Omnibus (GEO) datasets, we found that CCT6A was coamplified and coexpressed with EGFR in LUAD patients. CCT6A amplification was correlated with the unfavorable outcomes of LUAD. Moreover, CCT6A was upregulated in LUAD tissues, and CCT6A overexpression was correlated with the unfavorable relapse free survival or overall survival of LUAD. On the contrary, CCT6A was hypomethylated in LUAD, and CCT6A hypermethylation was correlated with the favorable overall survival of LUAD. Similar expression and methylation profiling of CCT6A were obtained in 479 lung normal tissues and 544 LUAD tissues collected from 11 independent datasets. In 1,462 LUAD patients from eight independent cohorts, CCT6A was also correlated with LUAD relapse-free survival or overall survival. Furthermore, CCT6A overexpression promoted the cell growth and invasion of LUAD. Identification of genes differentially expressed in CCT6A highly expressed LUAD patients revealed that CHCHD2 was the most correlated with CCT6A expression. CHCHD2 was coamplified with CCT6A. CHCHD2 was upregulated in LUAD tissues, and overexpression of CHCHD2 was correlated with the shorted relapse-free survival or overall survival of LUAD. Overall, our results revealed that CCT6A and CHCHD2 were coamplifying and coexpressing with EGFR and were correlated with the unfavorable clinical outcomes of LUAD.

Validation and analysis of expression, prognosis and immune infiltration of WNT gene family in non-small cell lung cancer

Early diagnosis and prognosis prediction of non-small cell lung cancer (NSCLC) have been challenging. Signaling cascades involving the Wingless-type (WNT) gene family play important biological roles and show prognostic value in various cancers, including NSCLC. On this basis, this study aimed to investigate the significance of WNTs in the prognosis and tumor immunity in NSCLC by comprehensive analysis. Expression and methylation levels of WNTs were obtained from the ONCOMINE, TIMER, and UALCAN. The dataset obtained from The Cancer Genome Atlas (TCGA) was utilized for prognostic analysis. cBioPortal was used to perform genetic alterations and correlation analysis of WNTs. R software was employed for functional enrichment and pathway analysis, partial statistics, and graph drawing. TRRUST was used to find key transcription factors. GEPIA was utilized for the analysis of expression, pathological staging, etc. Correlative analysis of immune infiltrates from TIMER. TISIDB was used for further immune infiltration validation analysis. Compared with that of normal tissues, WNT2/2B/3A/4/7A/9A/9B/11 expressions decreased, while WNT3/5B/6/7B/8B/10A/10B/16 expressions increased in lung adenocarcinoma (LUAD); WNT2/3A/7A/11 expressions were lessened, while WNT2B/3/5A/5B/6/7B/10A/10B/16 expressions were enhanced in squamous cell lung cancer (LUSC). Survival analysis revealed that highly expressed WNT2B and lowly expressed WNT7A predicted better prognostic outcomes in LUAD and LUSC. In the study of immune infiltration levels, WNT2, WNT9B, and WNT10A were positively correlated with six immune cells in LUAD; WNT1, WNT2, and WNT9B were positively correlated with six immune cells in LUSC, while WNT7B was negatively correlated. Our study indicated that WNT2B and WNT7A might have prognostic value in LUAD, and both of them might be important prognostic factors in LUSC and correlated to immune cell infiltration in LUAD and LUSC to a certain extent. Considering the prognostic value of WNT2B and WNT7A in NSCLC, we validated their mRNA and protein expression levels in NSCLC by performing qRT-PCR, western blot, and immunohistochemical staining on NSCLC pathological tissues and cell lines. This study may provide some direction for the subsequent exploration of the prognostic value of the WNTs and their role as biomarkers in NSCLC.

Identifying General Tumor and Specific Lung Cancer Biomarkers by Transcriptomic Analysis

The bioinformatic pipeline previously developed in our research laboratory is used to identify potential general and specific deregulated tumor genes and transcription factors related to the establishment and progression of tumoral diseases, now comparing lung cancer with other two types of cancer. Twenty microarray datasets were selected and analyzed separately to identify hub differentiated expressed genes and compared to identify all the deregulated genes and transcription factors in common between the three types of cancer and those unique to lung cancer. The winning DEGs analysis allowed to identify an important number of TFs deregulated in the majority of microarray datasets, which can become key biomarkers of general tumors and specific to lung cancer. A coexpression network was constructed for every dataset with all deregulated genes associated with lung cancer, according to DAVID’s tool enrichment analysis, and transcription factors capable of regulating them, according to oPOSSUM´s tool. Several genes and transcription factors are coexpressed in the networks, suggesting that they could be related to the establishment or progression of the tumoral pathology in any tissue and specifically in the lung. The comparison of the coexpression networks of lung cancer and other types of cancer allowed the identification of common connectivity patterns with deregulated genes and transcription factors correlated to important tumoral processes and signaling pathways that have not been studied yet to experimentally validate their role in lung cancer. The Kaplan–Meier estimator determined the association of thirteen deregulated top winning transcription factors with the survival of lung cancer patients. The coregulatory analysis identified two top winning transcription factors networks related to the regulatory control of gene expression in lung and breast cancer. Our transcriptomic analysis suggests that cancer has an important coregulatory network of transcription factors related to the acquisition of the hallmarks of cancer. Moreover, lung cancer has a group of genes and transcription factors unique to pulmonary tissue that are coexpressed during tumorigenesis and must be studied experimentally to fully understand their role in the pathogenesis within its very complex transcriptomic scenario. Therefore, the downstream bioinformatic analysis developed was able to identify a coregulatory metafirm of cancer in general and specific to lung cancer taking into account the great heterogeneity of the tumoral process at cellular and population levels.

Systematic analysis of expression profiles of HMGB family members for prognostic application in non-small cell lung cancer

Background: Lung cancer is a significant challenge to human health. Members of the high mobility group (HMG) superfamily (HMGB proteins) are implicated in a wide variety of physiological and pathophysiological processes, but the expression and prognostic value of HMGB family members in non-small cell lung cancer (NSCLC) have not been elucidated.

Methods: In this study, ONCOMINE, UALCAN, GEPIA, Kaplan–Meier Plotter, starBase, OncomiR databases, and GeneMANIA were utilized to evaluate the prognostic significance of HMGB family members in NSCLC.

Results: HMGB2/3 expression levels were higher in NSCLC patients. HMGB1 expression was higher in lung squamous cell carcinoma (LUSC) and was lower in lung adenocarcinoma (LUAD) tissue than in normal lung tissue. HMGB2 expression was related to cancer stage. Increased HMGB1 mRNA expression levels were associated with improved lung cancer prognosis, including overall survival (OS), first-progression survival (FP), and post-progression survival (PPS). There was no significant association between HMGB2 levels and prognostic indicators. HMGB3 expression was associated with poorer OS. GeneMANIA and GO/KEGG pathway analysis showed that HMGB family members mainly associated with chromosome condensation, regulation of chromatin organization, and nucleosome binding in NSCLC. HMGBs expression were closely correlated with infiltrating levels of specific types of immune cells in NSCLC, especially Th2 cells, Th17 cells, and mast cells. hsa-miR-25-3p, hsa-miR-374a-3p, and hsa-miR-93-5p were significantly positively correlated with HMGB1, HMGB2, and HMGB3, respectively. However, hsa-miR-30a-5p was predicted to significantly negatively regulate HMGB3 expression.

Conclusion: Our study revealed that HMGB1 is positively related to the improved prognosis in NSCLC, and demonstrate that HMGB3 might be a risk factor for poorer survival of NSCLC patients.

A novel prognostic signature of metastasis-associated genes and personalized therapeutic strategy for lung adenocarcinoma patients

Lung adenocarcinoma (LUAD) is a highly invasive and metastatic malignant tumor with high morbidity and mortality. This study aimed to construct a prognostic signature for LUAD patients based on metastasis-associated genes (MAGs). RNA expression profiles were downloaded from the Cancer Genome Atlas (TCGA) database. RRA method was applied to identify differentially expressed MAGs. A total of 192 significantly robust MAGs were determined among seven GEO datasets. MAGs were initially selected through the Lasso Cox regression analysis and 6 MAGs were included to construct a prognostic signature model. Transcriptome profile, patient prognosis, correlation between the risk score and clinicopathological features, immune cell infiltration characteristics, immunotherapy sensitivity and chemotherapy sensitivity differed between low- and high-risk groups after grouping according to median risk score. The reliability and applicability of the signature were further validated in the GSE31210, GSE50081 and GSE68465 cohort. CMap predicted 62 small molecule drugs on the base of the prognostic MAGs. Targeted drug staurosporine had hydrogen bonding with Gln-172 of SLC2A1, which is one of MAGs. Staurosporine could inhibit cell migration in A549 and H1299. We further verified mRNA and protein expression of 6 MAGs in A549 and H1299. The signature can serve as a promising prognostic tool and may provide a novel personalized therapeutic strategy for LUAD patients.

Expressional regulation of NKG2DLs is associated with the tumor development and shortened overall survival in lung adenocarcinoma

Natural killer group 2D ligands (NKG2DLs) are expressed on tumor cells as a ligand for Natural killer group 2D (NKG2D) receptors. NKG2DLs interact with NKG2D to induce immune cell-mediated cytotoxicity for eliminating tumors. Studies demonstrated that tumor cells can reduce NKG2DLs' expression to escape from anti-tumor immunity, leading to an aggressive cancer phenotype and poor prognosis in some cancers. However, these studies are limited and there is no comprehensive work on the regulation of NKG2DLs in lung adenocarcinoma (LUAD) which is one of the deadliest cancers worldwide. Here, we conducted an in silico analysis to evaluate the changes in NKG2DLs in LUAD by analyzing The Cancer Genome Atlas and the Gene Expression Omnibus datasets including tumor vs. normal comparisons, TNM stages, survival and infiltrating immune estimation profile. Results indicated that some members of NKG2DL were downregulated in LUAD as compared to normal samples. We determined that MICA (MHC class I polypeptide-related sequence A) was the most and significantly downregulated ligand among others and the results were nearly consistent with the different datasets which we used. Furthermore, survival analysis revealed that down-regulated MICA transcript expression might be one of the prognostic indicators of LUAD. Interestingly, according to the immune cell infiltrating analysis, there wasn't a direct correlation between the MICA transcript expression and immune cell infiltration, while for MICB there was. In addition, in genetic alteration, DNA methylation and miRNA analyses, we did not observe critical outcomes that would clarify the down-regulated MICA expression in detail. Regardless, this study is highly comprehensive and contributes valuable suggestions to further functional studies about the regulation of NKG2DLs and promising immunotherapeutic approaches in LUAD.