Recent advances in preclinical models for lung squamous cell carcinoma

HTR2B as a novel biomarker of chronic obstructive pulmonary disease with lung squamous cell carcinoma

Chronic obstructive pulmonary disease (COPD) is often associated with lung squamous cell carcinoma (LUSC), which has the same etiology (smoking, inflammation, oxidative stress, microenvironmental changes, and genetics). Smoking, inflammation, and airway remodeling are the most important and classical mechanisms of COPD comorbidity in LUSC patients. Cancer can occur during repeated airway damage and repair (airway remodeling). Changes in the inflammatory and immune microenvironments, which can cause malignant transformation of some cells, are currently being revealed in both LUSC and COPD patients. We obtained the GSE76925 dataset from the Gene Expression Omnibus database. Screening for possible COPD biomarkers was performed using the LASSO regression model and a random forest classifier. The compositional patterns of the immune cell fraction in COPD patients were determined using CIBERSORT. HTR2B expression was analyzed using validation datasets (GSE47460, GSE106986, and GSE1650). HTR2B expression in COPD cell models was determined via real-time quantitative PCR. Epithelial-mesenchymal transition (EMT) marker expression levels were determined after knocking down or overexpressing HTR2B. HTR2B function and mechanism in LUSC were analyzed with the Kaplan‒Meier plotter database. HTR2B expression was inhibited to detect changes in LUSC cell proliferation. A total of 1082 differentially expressed genes (DEGs) were identified in the GSE76925 dataset (371 genes were significantly upregulated, and 711 genes were significantly downregulated). Kyoto Encyclopedia of Genes and Genomes pathway enrichment analysis indicated that the DEGs were mainly enriched in the p53 signaling and β-alanine metabolism pathways. Gene Ontology enrichment analysis indicated that the DEGs were largely related to transcription initiation from the RNA polymerase I promoter and to the regulation of mononuclear cell proliferation. The LASSO regression model and random forest classifier results revealed that HTR2B, DPYS, FRY, and CD19 were key COPD genes. Immune cell infiltration analysis indicated that these genes were closely associated with immune cells. Analysis of the validation sets suggested that HTR2B was upregulated in COPD patients. HTR2B was significantly upregulated in COPD cell models, and its upregulation was associated with increased EMT marker expression. Compared with that in bronchial epithelial cells, HTR2B expression was upregulated in LUSC cells, and inhibiting HTR2B expression led to the inhibition of LUSC cell proliferation. In conclusions, HTR2B might be a new biomarker and therapeutic target in COPD patients with LUSC.

Cellular Origins and Lineage Plasticity in Cancer

All cancers arise from normal cells whose progeny acquire the cancer-initiating mutations and epigenetic modifications leading to frank tumorigenesis. The identity of those "cells-of-origin" has historically been a source of controversy across tumor types, as it has not been possible to witness the dynamic events giving rise to human tumors. Genetically engineered mouse models (GEMMs) of cancer provide an invaluable substitute, enabling researchers to interrogate the competence of various naive cellular compartments to initiate tumors in vivo. Researchers using these models have relied on lineage-specific promoters, knowledge of preneoplastic disease states in humans, and technical advances allowing more precise manipulations of the mouse germline. These approaches have given rise to the emerging view that multiple lineages within a given organ may generate tumors with similar histopathology. Here, we review some of the key studies leading to this conclusion in solid tumors and highlight the biological and clinical ramifications.

Significance of NKX2-1 as a biomarker for clinical prognosis, immune infiltration, and drug therapy in lung squamous cell carcinoma

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.

Identification of SLC2A1 as a predictive biomarker for survival and response to immunotherapy in lung squamous cell carcinoma

BACKGROUND

As one of the common subtypes of non-small lung cancer, lung squamous cell carcinoma (LUSC) patients with advanced stage have few choices of treatment strategies. Therefore, it is urgent to discover genes that are associated with the survival and efficacy of immunotherapies.

METHOD

Differential gene expression analyses were conducted using TCGA LUSC bulk-sequencing and single-cell RNA-sequencing data. Prognostic genes were identified from the TCGA LUSC cohort. Protein expression validation and survival analyses were performed. Experiments were conducted to explore the underlying mechanisms. In addition, the correlation between gene expression and pathological response to adjuvant immunochemotherapy was also investigated.

RESULTS

After a series of bioinformatic analyses, solute carrier family 2 member 1(SLC2A1), encoding glucose transporter-1 (GLUT1), was found to be differentially expressed between tumor and normal tissues. GLUT1 was subsequently identified as an independent prognostic factor for LUSC. GSEA analysis revealed the glycolysis metabolism pathway of KEGG enriched in SLC2A1high tumor tissues. LASSO analyses revealed that tumor tissues with high expression of SLC2A1 were associated with high levels of protein lactylation. We found that SLC2A1 was preferentially expressed by SPP1+ macrophages in the tumor microenvironment, and the expression of SLC2A1 was associated with the abundance of SPP1+ macrophages. Immunofluorescence demonstrated GLUT1 and HIF1α colocalization in tumor-infiltrating macrophages. In vitro experiments showed HIF-1α-induced macrophage polarization under hypoxia, and GLUT1 inhibition blocked this polarization. In addition, SLC2A1 was negatively associated with the common immune checkpoint molecules, such as programmed cell death 1(PD-1), T cell immunoreceptor with Ig and ITIM domains (TIGIT), cytotoxic T-lymphocyte associated protein 4 (CTLA4) and lymphocyte activating 3 (LAG3), while showed a positive association with CD44. Finally, we observed that there was a significant correlation between pre-adjuvant-treatment GLUT1 expression and the pathological response.

CONCLUSION

SLC2A1 expression was differentially upregulated in tumor tissues, and elevated GLUT1 expression was associated with worse survival and poor pathological response to adjuvant immunochemotherapy. Upregulation of GLUT1 promoted macrophage polarization into the M2 phenotype. The findings will contribute to guiding the treatment selection for LUSC patients and providing personalized immunotherapy strategies.

USP13 drives lung squamous cell carcinoma by switching lung club cell lineage plasticity

Lung squamous cell carcinoma (LUSC) is associated with high mortality and limited targeted therapies. USP13 is one of the most amplified genes in LUSC, yet its role in lung cancer is largely unknown. Here, we established a novel mouse model of LUSC by overexpressing USP13 on KrasG12D/+; Trp53flox/flox background (KPU). KPU-driven lung squamous tumors faithfully recapitulate key pathohistological, molecular features, and cellular pathways of human LUSC. We found that USP13 altered lineage-determining factors such as NKX2-1 and SOX2 in club cells of the airway and reinforced the fate of club cells to squamous carcinoma development. We showed a strong molecular association between USP13 and c-MYC, leading to the upregulation of squamous programs in murine and human lung cancer cells. Collectively, our data demonstrate that USP13 is a molecular driver of lineage plasticity in club cells and provide mechanistic insight that may have potential implications for the treatment of LUSC.

MYH16 upregulation is associated with lung adenocarcinoma aggressiveness and immune infiltration

Myosin heavy chain 16 (MYH16) may significantly affect cell cycle progression. Nevertheless, there is a lack of evidence about the clinical relevance of MYH16 upregulation in pan cancers, including lung adenocarcinoma (LUAD). MYH16 expression patterns were evaluated in various bioinformatics databases using The Cancer Genome Atlas data set. Clinical and pathological factor data were employed to risk-stratify patients. The Kaplan-Meier plotter approach was used to estimate survival rates. Tumor immune infiltration was explored via the TIMER tool, and gene set enrichment analysis (GSEA) was used to identify the pathways involved in MYH16 upregulation. The results showed that MYH16 was abnormally upregulated in pan cancers, including LUAD. MYH16 expression induction in LUAD was found to be related to the tumor stage. Furthermore, MYH16 upregulation was correlated with LUAD development and worse overall survival, particularly in women. Notably, MYH16 overexpression in LUAD tissues corresponded to the amount of immune infiltration in the tumor. Additionally, univariate Cox hazard regression analysis revealed that MYH16 may be an independent prognostic indicator for LUAD. Furthermore, a nomogram was constructed according to MYH16 expression and clinical characteristics. BMP6 expression deficiency may be a key factor contributing to MYH16 upregulation in LUAD. Finally, GSEA demonstrated that MYH16 might mediate meiosis and gene silencing through RNA signaling pathways. This study, for the first time, showed that MYH16 upregulation in LUAD is associated with various risk factors, increased cancer aggressiveness, enhanced infiltration of tumor immune cells, and reduced survival rates.

Identification of prognosis-related lncRNAs and cell validation in lung squamous cell carcinoma based on TCGA data

Objective

To discern long non-coding RNAs (lncRNAs) with prognostic relevance in the context of lung squamous cell carcinoma (LUSC), we intend to predict target genes by leveraging The Cancer Genome Atlas (TCGA) repository. Subsequently, we aim to investigate the proliferative potential of critical lncRNAs within the LUSC milieu.

Methods

DESeq2 was employed to identify differentially expressed genes within the TCGA database. Following this, we utilized both univariate and multivariate Cox regression analyses to identify lncRNAs with prognostic relevance. Noteworthy lncRNAs were selected for validation in cell lines. The intracellular localization of these lncRNAs was ascertained through nucleocytoplasmic isolation experiments. Additionally, the impact of these lncRNAs on cellular proliferation, invasion, and migration capabilities was investigated using an Antisense oligonucleotides (ASO) knockdown system.

Results

Multivariate Cox regression identified a total of 12 candidate genes, consisting of seven downregulated lncRNAs (BRE-AS1, CCL15-CCL14, DNMBP-AS1, LINC00482, LOC100129034, MIR22HG, PRR26) and five upregulated lncRNAs (FAM83A-AS1, LINC00628, LINC00923, LINC01341, LOC100130691). The target genes associated with these lncRNAs exhibit significant enrichment within diverse biological pathways, including metabolic processes, cancer pathways, MAPK signaling, PI3K-Akt signaling, protein binding, cellular components, cellular transformation, and other functional categories. Furthermore, nucleocytoplasmic fractionation experiments demonstrated that LINC00923 and LINC01341 are predominantly localized within the cellular nucleus. Subsequent investigations utilizing CCK-8 assays and colony formation assays revealed that the knockdown of LINC00923 and LINC01341 effectively suppressed the proliferation of H226 and H1703 cells. Additionally, transwell assays showed that knockdown of LINC00923 and LINC01341 significantly attenuated the invasive and migratory capacities of H226 and H1703 cells.

Conclusion

This study has identified 12 candidate lncRNA associated with prognostic implications, among which LINC00923 and LINC01341 exhibit potential as markers for the prediction of LUSC outcomes.

Machine learning analysis of lung squamous cell carcinoma gene expression datasets reveals novel prognostic signatures

BACKGROUND

Lung squamous cell carcinoma (LUSC) patients are often diagnosed at an advanced stage and have poor prognoses. Thus, identifying novel biomarkers for the LUSC is of utmost importance.

METHODS

Multiple datasets from the NCBI-GEO repository were obtained and merged to construct the complete dataset. We also constructed a subset from this complete dataset with only known cancer driver genes. Further, machine learning classifiers were employed to obtain the best features from both datasets. Simultaneously, we perform differential gene expression analysis. Furthermore, survival and enrichment analyses were performed.

RESULTS

The kNN classifier performed comparatively better on the complete and driver datasets' top 40 and 50 gene features, respectively. Out of these 90 gene features, 35 were found to be differentially regulated. Lasso-penalized Cox regression further reduced the number of genes to eight. The median risk score of these eight genes significantly stratified the patients, and low-risk patients have significantly better overall survival. We validated the robust performance of these eight genes on the TCGA dataset. Pathway enrichment analysis identified that these genes are associated with cell cycle, cell proliferation, and migration.

CONCLUSION

This study demonstrates that an integrated approach involving machine learning and system biology may effectively identify novel biomarkers for LUSC.

Subtype classification based on t cell proliferation-related regulator genes and risk model for predicting outcomes of lung adenocarcinoma

Background

Lung adenocarcinoma (LUAD), the major lung cancer histotype, represents 40% lung cancers. Currently, outcomes are remarkably different in LUAD patients with similar AJCC/UICC-TNM features. T cell proliferation-related regulator genes (TPRGs) relate to the proliferation, activity and function of T cells and tumor progression. The values of TPRGs in classifying LUAD patients and predicting outcomes remain unknown.

Methods

Gene expression profile and corresponding clinical data were downloaded from TCGA and the GEO databases. We systematically analyzed the expression profile characteristics of 35 TPRGs in LUAD patients and investigated the differences in overall survival (OS), biology pathway, immunity and somatic mutation between different TPRGs-related subtypes. Subsequently, we constructed a TPRGs-related risk model in TCGA cohort to quantify risk scores using LASSO cox regression analysis and then validated this risk model in two GEO cohorts. LUAD patients were divided into high- and low-risk subtypes according to the median risk score. We systematically compared the biology pathway, immunity, somatic mutation and drug susceptibility between the two risk subtypes. Finally, we validate biological functions of two TPRGs-encoded proteins (DCLRE1B and HOMER1) in LUAD cells A549.

Results

We identified different TPRGs-related subtypes (including cluster 1/cluster A and its counterpart cluster 2/cluster B). Compared to the cluster 1/cluster A subtype, cluster 2/cluster B subtype tended to have a prominent survival advantage with an immunosuppressive microenvironment and a higher somatic mutation frequency. Then, we constructed a TPRGs-related 6-gene risk model. The high-risk subtype characterized by higher somatic mutation frequency and lower immunotherapy response had a worse prognosis. This risk model was an independent prognostic factor and showed to be reliable and accurate for LUAD classification. Furthermore, subtypes with different risk scores were significantly associated with drug sensitivity. DCLRE1B and HOMER1 suppressed cell proliferation, migration and invasion in LUAD cells A549, which was in line with their prognostic values.

Conclusion

We construed a novel stratification model of LUAD based on TPRGs, which can accurately and reliably predict the prognosis and might be used as a predictive tool for LUAD patients.

Using patient-derived organoids to predict locally advanced or metastatic lung cancer tumor response: A real-world study

Predicting the clinical response to chemotherapeutic or targeted treatment in patients with locally advanced or metastatic lung cancer requires an accurate and affordable tool. Tumor organoids are a potential approach in precision medicine for predicting the clinical response to treatment. However, their clinical application in lung cancer has rarely been reported because of the difficulty in generating pure tumor organoids. In this study, we have generated 214 cancer organoids from 107 patients, of which 212 are lung cancer organoids (LCOs), primarily derived from malignant serous effusions. LCO-based drug sensitivity tests (LCO-DSTs) for chemotherapy and targeted therapy have been performed in a real-world study to predict the clinical response to the respective treatment. LCO-DSTs accurately predict the clinical response to treatment in this cohort of patients with advanced lung cancer. In conclusion, LCO-DST is a promising precision medicine tool in treating of advanced lung cancer.

KMT2D deficiency drives lung squamous cell carcinoma and hypersensitivity to RTK-RAS inhibition

Lung squamous cell carcinoma (LUSC) represents a major subtype of lung cancer with limited treatment options. KMT2D is one of the most frequently mutated genes in LUSC (>20%), and yet its role in LUSC oncogenesis remains unknown. Here, we identify KMT2D as a key regulator of LUSC tumorigenesis wherein Kmt2d deletion transforms lung basal cell organoids to LUSC. Kmt2d loss increases activation of receptor tyrosine kinases (RTKs), EGFR and ERBB2, partly through reprogramming the chromatin landscape to repress the expression of protein tyrosine phosphatases. These events provoke a robust elevation in the oncogenic RTK-RAS signaling. Combining SHP2 inhibitor SHP099 and pan-ERBB inhibitor afatinib inhibits lung tumor growth in Kmt2d-deficient LUSC murine models and in patient-derived xenografts (PDXs) harboring KMT2D mutations. Our study identifies KMT2D as a pivotal epigenetic modulator for LUSC oncogenesis and suggests that KMT2D loss renders LUSC therapeutically vulnerable to RTK-RAS inhibition.

Identification of tumour-infiltrating myeloid subsets associated with overall survival in lung squamous cell carcinoma

Lung squamous cell carcinoma (LUSC) is a primary subtype of lung cancer with limited therapeutic options and poor prognosis, and tumour-infiltrating myeloid cells (TIMs) are key regulators of LUSC. However, the correlation between the abundance of TIM subtypes and clinical outcomes of LUSC remains unexplored. This study aimed to develop and validate a prognostic model for low- and high-risk patients with LUSC based on myeloid cell microenvironments. TIM markers in the tumoural (T) and stromal (S) regions were quantified using immunohistochemistry for 502 LUSC patients. L1-penalized Cox regression was used to develop a myeloid survival score (MSS) model based on the training cohort, followed by validation in distinct cohorts from multiple centres. RNA sequencing and immunostaining were used to examine the mechanisms of myeloid cells in LUSC progression and predict potential drug targets and therapeutic agents. Of the 12 myeloid markers, CD163T, CD163S, and S100A12T were highly associated with overall survival (OS) in LUSC patients. The MSS of the three myeloid signatures accurately categorized LUSC patients into risk categories, with an observable difference in OS between the training and validation cohorts. Tumours with high MSS were associated with enhanced antioxidative ability and hedgehog signalling and a shift to a more pro-tumorigenic microenvironment, accompanied by a reduced tumour cell immunogenicity and increased CD8⁺ T cell exhaustion patterns. Additionally, in high-risk patients, potential drug targets and compounds regulating hedgehog signalling were identified. Our study provides the first prognostic myeloid signature for LUSC, which may help advance precision medicine. © 2022 The Authors. The Journal of Pathology published by John Wiley & Sons Ltd on behalf of The Pathological Society of Great Britain and Ireland.

Involvement of CHRNA6 in the Immune Response in Lung Squamous Cell Carcinoma and its Potential as a Drug Target for the Disease

BACKGROUND

Lung squamous cell carcinoma (LUSC) is a subtype of lung cancer with a poor prognosis and limited treatment options. Previous studies show that some components of the cholinergic pathway may play important roles in the tumorigenesis of lung cancer, including LUSC.

OBJECTIVE

The purpose of this study is to investigate the involvement of cholinergic genes in immune infiltration in LUSC, and identify the key genes in the pathway and analyze their potential as targets for LUSC treatment and novel drugs.

METHODS

We first screened the cholinergic genes associated with immune infiltration in LUSC based on transcriptomic samples and explored the correlation between the key genes and immune infiltrating cells and immune pathways. Then, we assessed the effect of immunotherapeutic response in the high and low-expression groups of key genes in vitro. And finally, we screened potential drugs for the treatment of LUSC.

RESULTS

We found that the expression of CHRNA6, the gene encoding the α6 subunit of nicotinic acetylcholine receptors (nAChR), was significantly correlated with the proportion of immune infiltrating cells in LUSC, and the high expression level of the gene was associated with poor prognosis of the disease. Also, the proportion of Tregs, M1 macrophages, and resting mast cells was correlated with the expression of CHRNA6. In addition, LUSC patients with higher CHRNA6 expression levels had better immunotherapy responses. Furthermore, we found that the drugs, i.e., adavosertib, varbulin and pyrazoloacridine, had a strong affinity with CHRNA6, with adavosertib binding most stably with the protein.

CONCLUSION

CHRNA6 may be associated with immune infiltration in LUSC and affects patient prognosis and immunotherapeutic response by regulating immune cells and immune pathways. In addition, adavosertib may be a potential drug for the treatment of LUSC.

miR-139-5p and miR-451a as a Diagnostic Biomarker in LUSC

Background

Lung squamous cell carcinoma (LUSC) is a type of lung cancer that originates from segmental or subsegmental bronchial mucosa. There is evidence that miRNA plays an important role in the occurrence and progression of tumors.

Methods

In this study, plasma samples of patients with early LUSC and healthy volunteers were subjected to miRNA sequencing, and the levels of differentially expressed miRNAs (DEMs) in LUSC tissues were analyzed using R language. Cox regression and Kaplan-Meier (K-M) survival curve analyses were performed to determine the relationship between DEMs and prognosis in LUSC, and PCR method was verified for the plasma expression level of DEMs in patients with LUSC. The levels of CYFRA21-1 and SCC-Ag in plasma were measured, and area under curve (AUC) was used to evaluate the diagnostic value of the DEMs.

Results

A total of 21 DEMs were screened out by sequencing. The expression levels of DEMs in tissue samples in the TCGA database were analyzed, and four DEMs with consistent expression levels were further screened from plasma and tissue samples. Regression analysis and K-M curve were performed to select two DEMs (miR-139-5p, miR-451a) that were correlated with the prognosis. PCR verification results showed that the levels of miR-451a and miR-139-5p were low in patients, and the level of miR-139-5p in late stages III & IV with the patients of LUSC was higher than that in stages I & II. The AUC values of the four indicators (SCC-Ag, CYFRA21-1, miR-451a and miR-139-5p) in the diagnosis of LUSC, early and late cases were 0.884, 0.935 and 0.778, respectively.

Conclusion

The detection of miR-139-5p and miR-451a levels in plasma has a certain potential in the non-invasive diagnosis, especially in patients with early stages of LUSC.

Dissecting the Immunological Profiles in NSD3-Amplified LUSC through Integrative Multi-Scale Analyses

Simple Summary

Although cancer immunotherapy has become a “game changer” in treating LUSC patients, challenges still prevail due to the heterogeneous response and insufficient predictive biomarkers. Recently, NSD3, a neighbor gene of FGFR1, was identified as a key genetic driver of LUSC tumorigenesis. While previous studies have reported the relevance of NSD3 in innate antiviral immunity, the association of NSD3 with the TIME remains an open question requiring further investigation. We first show that NSD3 gene amplification is associated with an immune-desert phenotype and correlated with a worse immunotherapy outcome. Further molecular characterizations pinpoint that the high activity of UPR signaling might be a pivotal mediator for the non-inflamed TME phenotype of NSD3-amplified LUSC. Concordantly, the pharmaco-transcriptomic correlation analysis indicated that the intervention of UPR signaling could be a promising synthetic lethality target for NSD3-amplified LUSC. Our findings reveal a previously unappreciated immunological role for NSD3 in LUSC pathobiology and provide therapeutic rationales for this daunting disease.

Abstract

The histone H3 lysine 36 (H3K36) methyltransferase NSD3, a neighboring gene of FGFR1, has been identified as a critical genetic driver of lung squamous cell carcinoma (LUSC). However, the molecular characteristics, especially the immunological roles of NSD3 in driving carcinogenesis, are poorly understood. In this study, we systematically integrated multi-omics data (e.g., genome, transcriptome, proteome, and TMA array) to dissect the immunological profiles in NSD3-amplified LUSC. Next, pharmaco-transcriptomic correlation analysis was implemented to identify the molecular underpinnings and therapeutic vulnerabilities in LUSC. We revealed that NSD3-amplified LUSC presents a non-inflamed tumor immune microenvironment (TIME) state in multiple independent LUSC patient cohorts. Predictably, elevated NSD3 expression was correlated with a worse immunotherapy outcome. Further molecular characterizations revealed that the high activity of unfolded protein response (UPR) signaling might be a pivotal mediator for the non-immunogenic phenotype of NSD3-amplified LUSC. Concordantly, we showed that NSD3-amplified LUSCs exhibited a more sensitive phenotype to compounds targeting UPR branches than the wild-type group. In brief, our multi-level analyses point to a previously unappreciated immunological role for NSD3 and provide therapeutic rationales for NSD3-amplified squamous lung cancer.

Identification of adenoid subtype characterized with immune-escaped phenotype in lung squamous carcinoma based on transcriptomics

Non-small cell lung cancer (NSCLC) is a heterogeneous disease, and its demarcation contributes to various therapeutic outcomes. However, a small subset of tumors shows different molecular features that are in contradiction with pathological classification. Unsupervised clustering was performed to subtype NSCLC using the transcriptome data from the TCGA database. Next, immune microenvironment features of lung adenocarcinoma (LUAD), lung squamous carcinoma (LUSC), and lung adenoid squamous carcinoma (LASC) were characterized. In addition, diagnostic biomarkers to demarcate LASC among LUSC were screened using weighted gene co-expression network analysis (WGCNA) and validated by the in-house cohort. LASC was identified as a novel subtype with adenoid transcriptomic features in LUSC, which exhibited the most immuno-escaped phenotype among all NSCLC subtypes. In addition, FOLR1 was identified as a biomarker for LASC discrimination using the WGCNA analysis, and its diagnostic value was validated by the in-house cohort. Moreover, FOLR1 was related to immuno-escaped tumors in LUSC but not in LUAD. Overall, we proposed a novel typing strategy in NSCLC and identified FOLR1 as a biomarker for LASC discrimination.

Serum-derived piR-hsa-164586 of extracellular vesicles as a novel biomarker for early diagnosis of non-small cell lung cancer

Non-small cell lung cancer (NSCLC) is a major cause of death in those with malignant tumors. To achieve the early diagnosis of NSCLC, we investigated serum-derived Piwi-interacting RNA (piRNA) of extracellular vesicles to filter diagnostic biomarkers for NSCLC. High-throughput sequencing from cancerous tissues and adjacent noncancerous tissues in patients with NSCLC was first applied to recognize candidate piRNAs as diagnostic biomarkers. These screened piRNAs were further validated in 115 patients (including 95 cases in stage I) and 47 healthy individuals using quantitative real-time PCR (qRT-PCR). We showed that piR-hsa-164586 was significantly upregulated compared with paracancerous tissues and extracellular vesicles from the serum samples of healthy individuals. Moreover, the area under the curve (AUC) value of piR-hsa-164586 was 0.623 and 0.624 to distinguish patients with all stages or stage I of NSCLC, respectively, from healthy individuals. The diagnostic performance of piR-hsa-164586 was greatly improved compared with the cytokeratin-19-fragment (CYFRA21-1). Additionally, piR-hs-164586 was associated with the clinical characteristics of patients with NSCLC. Its expression was associated with the age and TNM stage of patients with NSCLC, indicating that it can serve as an effective and promising biomarker for the early diagnosis of NSCLC.

Comprehensive analysis of molecular features, prognostic values, and immune landscape association of m6A-regulated immune-related lncRNAs in smoking-associated lung squamous cell carcinoma

Lung squamous cell carcinoma (LUSC) is the second most common histopathological subtype of lung cancer, and smoking is the leading cause of this type of cancer. However, the critical factors that directly affect the survival rate and sensitivity to immunotherapy of smoking LUSC patients are still unknown. Previous studies have highlighted the role of N6-methyladenosine (m6A) RNA modification, the most common epigenetic modification in eukaryotic species, together with immune-related long non-coding RNAs (lncRNAs) in promoting the development and progression of tumors. Thus, elucidating m6A-modified immune lncRNAs in LUSC patients with smoking history is vital. In this study, we described the expression and mutation features of the 24 m6A-related regulators in the smoking-associated LUSC cohort from The Cancer Genome Atlas (TCGA) database. Then, two distinct subtypes based on the expression levels of the prognostic m6A-regulated immune lncRNAs were defined, and differentially expressed genes (DEGs) between the subtypes were identified. The distributions of clinical characteristics and the tumor microenvironment (TME) between clusters were analyzed. Finally, we established a lncRNA-associated risk model and exhaustively clarified the clinical features, prognosis, immune landscape, and drug sensitivity on the basis of this scoring system. Our findings give insight into potential mechanisms of LUSC tumorigenesis and development and provide new ideas in offering LUSC patients with individual and effective immunotherapies.

A Novel Small-Molecule Inhibitor of SREBP-1 Based on Natural Product Monomers Upregulates the Sensitivity of Lung Squamous Cell Carcinoma Cells to Antitumor Drugs

The transcription factor, sterol regulatory element binding protein 1 (SREBP-1), plays important roles in modulating the proliferation, metastasis, or resistance to antitumor agents by promoting cellular lipid metabolism and related cellular glucose-uptake/Warburg Effect. However, the underlying mechanism of SREBP-1 regulating the proliferation or drug-resistance in lung squamous cell carcinoma (LUSC) and the therapeutic strategies targeted to SREBP-1 in LUSC remain unclear. In this study, SREBP-1 was highly expressed in LUSC tissues, compared with the paired non-tumor tissues (the para-tumor tissues). A novel small-molecule inhibitor of SREBP-1, MSI-1 (Ma’s inhibitor of SREBP-1), based on natural product monomers, was identified by screening the database of natural products. Treatment with MSI-1 suppressed the activation of SREBP-1-related pathways and the Warburg effect of LUSC cells, as indicated by decreased glucose uptake or glycolysis. Moreover, treatment of MSI-1 enhanced the sensitivity of LUSC cells to antitumor agents. The specificity of MSI-1 on SREBP-1 was confirmed by molecular docking and point-mutation of SPEBP-1. Therefore, MSI-1 improved our understanding of SREBP-1 and provided additional options for the treatment of LUSC.

Multiple datasets to explore the tumor microenvironment of cutaneous squamous cell carcinoma

BACKGROUND

Cutaneous squamous cell carcinoma (cSCC) is one of the most frequent types of cutaneous cancer. The composition and heterogeneity of the tumor microenvironment significantly impact patient prognosis and the ability to practice precision therapy. However, no research has been conducted to examine the design of the tumor microenvironment and its interactions with cSCC.

MATERIAL AND METHODS

We retrieved the datasets GSE42677 and GSE45164 from the GEO public database, integrated them, and analyzed them using the SVA method. We then screened the core genes using the WGCNA network and LASSO regression and checked the model's stability using the ROC curve. Finally, we performed enrichment and correlation analyses on the core genes.

RESULTS

We identified four genes as core cSCC genes: DTYMK, CDCA8, PTTG1 and MAD2L1, and discovered that RORA, RORB and RORC were the primary regulators in the gene set. The GO semantic similarity analysis results indicated that CDCA8 and PTTG1 were the two most essential genes among the four core genes. The results of correlation analysis demonstrated that PTTG1 and HLA-DMA, CDCA8 and HLA-DQB2 were significantly correlated.

CONCLUSIONS

Examining the expression levels of four primary genes in cSCC aids in our understanding of the disease's pathophysiology. Additionally, the core genes were found to be highly related with immune regulatory genes, suggesting novel avenues for cSCC prevention and treatment.

ATP8B1 Knockdown Activated the Choline Metabolism Pathway and Induced High-Level Intracellular REDOX Homeostasis in Lung Squamous Cell Carcinoma

Simple Summary

We found that low expression of ATP8B1 was associated with poor prognosis, and involved in the dysregulation of glutathione (GSH) synthesis and choline metabolism in lung squamous cell carcinoma (LUSC) samples of The Cancer Genome Atlas (TCGA) and Tianjin Medical University Cancer Institute and Hospital (TJMUCH) cohort. We further constructed ATP8B1 knockdown of LUSC cell lines H520^SH-ATP8B1 and SK-MES-1^SH-ATP8B1 to investigate how ATP8B1 knockdown promoted cell proliferation, migration, and invasion in vitro and in vivo via upregulation of the CHKA-dependent choline metabolism pathway. We identified that ATP8B1 knockdown and CHKA upregulation can lead to mitochondrial dysfunction and high reduction-oxidation (REDOX) homeostasis, which may be involved in the roles of cardiolipin in maintaining mitochondrial dynamics and phospholipid homeostasis. Therefore, we proposed ATP8B1 as a novel predictive biomarker in LUSC and targeting ATP8B1-driven specific metabolic disorder might be a promising therapeutic strategy for LUSC.

Abstract

The flippase ATPase class I type 8b member 1 (ATP8B1) is essential for maintaining the stability and polarity of the epithelial membrane and can translocate specific phospholipids from the outer membrane to the inner membrane of the cell. Although ATP8B1 plays important roles in cell functions, ATP8B1 has been poorly studied in tumors and its prognostic value in patients with lung squamous cell carcinoma (LUSC) remains unclear. By investigating the whole genomic expression profiles of LUSC samples from The Cancer Genome Atlas (TCGA) database and Tianjin Medical University Cancer Institute and Hospital (TJMUCH) cohort, we found that low expression of ATP8B1 was associated with poor prognosis of LUSC patients. The results from cellular experiments and a xenograft-bearing mice model indicated that ATP8B1 knockdown firstly induced mitochondrial dysfunction and promoted ROS production. Secondly, ATP8B1 knockdown promoted glutathione synthesis via upregulation of the CHKA-dependent choline metabolism pathway, therefore producing and maintaining high-level intracellular REDOX homeostasis to aggravate carcinogenesis and progression of LUSC. In summary, we proposed ATP8B1 as a novel predictive biomarker in LUSC and targeting ATP8B1-driven specific metabolic disorder might be a promising therapeutic strategy for LUSC.

KRAS Mutations in Squamous Cell Carcinomas of the Lung

KRAS is one of the most commonly mutated oncogenes in cancer, enabling tumor proliferation and maintenance. After various approaches to target KRAS have failed over the past decades, the first specific inhibitor of the p.G12C mutation of KRAS was recently approved by the FDA after showing promising results in adenocarcinomas of the lung and other solid tumors. Lung cancer, the most common cancer worldwide, is a promising use case for these new therapies, as adenocarcinomas in particular frequently harbor KRAS mutations. However, in squamous cell carcinoma (SCC) of the lung, KRAS mutations are rare and their impact on clinical outcome is poorly understood. In this review, we discuss the current knowledge on the prevalence and prognostic and predictive significance of KRAS mutations in the context of SCC.

Low-Coverage Whole Genome Sequencing Using Laser Capture Microscopy with Combined Digital Droplet PCR: An Effective Tool to Study Copy Number and Kras Mutations in Early Lung Adenocarcinoma Development

Defining detailed genomic characterization of early tumor progression is critical to identifying key regulators and pathways in carcinogenesis as potentially druggable targets. In human lung cancer, work to characterize early cancer development has mainly focused on squamous cancer, as the earliest lesions are more proximal in the airways and often accessible by repeated bronchoscopy. Adenocarcinomas are typically located distally in the lung, limiting accessibility for biopsy of pre-malignant and early stages. Mouse lung cancer models recapitulate many human genomic features and provide a model for tumorigenesis with pre-malignant atypical adenomatous hyperplasia and in situ adenocarcinomas often developing contemporaneously within the same animal. Here, we combined tissue characterization and collection by laser capture microscopy (LCM) with digital droplet PCR (ddPCR) and low-coverage whole genome sequencing (LC-WGS). ddPCR can be used to identify specific missense mutations in Kras (Kirsten rat sarcoma viral oncogene homolog, here focused on Kras Q61) and estimate the percentage of mutation predominance. LC-WGS is a cost-effective method to infer localized copy number alterations (CNAs) across the genome using low-input DNA. Combining these methods, the histological stage of lung cancer can be correlated with appearance of Kras mutations and CNAs. The utility of this approach is adaptable to other mouse models of human cancer.