Cell Trajectory-Related Genes of Lung Adenocarcinoma Predict Tumor Immune Microenvironment and Prognosis of Patients

Human epididymis protein 4, a novel potential biomarker for diagnostic and prognosis monitoring of lung cancer

OBJECTIVE

This study aimed to explore the application value of human epididymis protein 4 (HE4) in diagnosing and monitoring the prognosis of lung cancer.

METHODS

First, TCGA (The Cancer Genome Atlas) databases were used to analyze whey-acidic-protein 4-disulfide bond core domain 2 (WFDC2) gene expression levels in lung cancer tissues. Then, a total of 160 individuals were enrolled, categorized into three groups: the lung cancer group (n = 80), the benign lesions group (n = 40), and the healthy controls group (n = 40). Serum HE4 levels and other biomarkers were quantified using an electro-chemiluminescent immunoassay. Additionally, the expression of HE4 in tissues was analyzed through immunohistochemistry (IHC). In vitro cultures of human airway epithelial (human bronchial epithelial [HBE]) cells and various lung cancer cell lines (SPC/PC9/A594/H520) were utilized to detect HE4 levels via western blot (WB).

RESULTS

Analysis of the TCGA and UALCAN (The University of Alabama at Birmingham Cancer Data Analysis Portal) databases showed that WFDC2 gene expression levels were upregulated in lung cancer tissues (p < 0.01). Compared with the control group and the benign group, HE4 was significantly higher in the serum of patients with lung cancer (p < 0.001). Receiver operating characteristic (ROC) analysis confirmed that HE4 had better diagnostic efficacy than classical markers in the differential diagnosis of lung cancer and benign lesions and had the highest diagnostic value in lung adenocarcinoma (area under the ROC curve [AUC] = 0.826). HE4 increased in early lung cancer and positively correlated with poor prognosis (p < 0.001). Moreover, the results of WB and IHC revealed that the expression of HE4 was increased in lung cancer cells (SPC/A549/H520) and lung cancer tissues but decreased in PC9 cells with a lack of exon EGFR19 (p < 0.05).

CONCLUSION

Serum HE4 emerges as a promising novel biomarker for the diagnosis and prognosis assessment of lung cancer.

Unraveling the tumor immune microenvironment of lung adenocarcinoma using single-cell RNA sequencing

Tumor immune microenvironment (TIME) and its indications for lung cancer patient prognosis and therapeutic response have become new hotspots in cancer research in recent years. Tumor cells, immune cells, various regulatory factors, and their interactions in the TIME have been suggested to commonly influence lung cancer development and therapeutic outcome. The heterogeneity of TIME is composed of dynamic immune-related components, including various cancer cells, immune cells, cytokine/chemokine environments, cytotoxic activity, or immunosuppressive factors. The specific composition of cell subtypes may facilitate or hamper the response to immunotherapy and influence patient prognosis. Various markers have been found to stratify the patient prognosis or predict the therapeutic outcome. In this article, we systematically reviewed the recent advancement of TIME studies in lung adenocarcinoma (LUAD) using single-cell RNA sequencing (scRNA-seq) techniques, with specific focuses on the roles of TIME in LUAD development, TIME heterogeneity, indications of TIME in patient prognosis and therapeutic response during immunotherapy and drug resistance. The main findings in TIME heterogeneity and relevant markers or models for prognosis stratification and response prediction have been summarized. We hope that this review provides an overview of TIME status in LUAD and an inspiration for future development of strategies and biomarkers in LUAD treatment.

A novel iTreg-related signature for prognostic prediction in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the most common subtype of lung cancer. Most patients are diagnosed at an advanced stage, therefore it is crucial to identify novel prognostic biomarkers for LUAD. As important regulatory cells, inducible regulatory T cells (iTregs) play a vital role in immune suppression and are important for the maintenance of immune homeostasis. This study explored the prognostic value and therapeutic effects of iTreg-related genes in LUAD. Data for LUAD patients, including immune infiltration data, RNA sequencing data, and clinical features, were acquired from The Cancer Genome Atlas, Gene Expression Omnibus, and Tumor Immune Single-cell Hub 2 databases. Immune-related subgroups with different infiltration patterns and iTreg-related genes were identified through univariate and multivariate Cox regression analyses and weighted correlation network analysis. Functional enrichment analyses were performed to explore the underlying mechanisms of iTreg-related genes. A prognostic risk signature was constructed using Cox regression analysis with the least absolute shrinkage and selection operator penalty. The ESTIMATE algorithm was applied to determine the immune status of LUAD patients. We applied the constructed signature to predict chemosensitivity and performed single-cell RNA sequencing analysis. The infiltration of iTregs was identified as an independent factor for predicting patient outcomes. We constructed a prognostic signature based on seven iTreg-related genes (GIMAP5, SLA, MS4A7, ZNF366, POU2AF1, MRPL12, and COL5A1), which was applied to subdivide patients into high- and low-risk subgroups. Our results revealed that patients in the iTreg-related low-risk subgroup had a better prognosis and possibly greater sensitivity to traditional chemotherapy. Our study provides a novel iTreg-related signature to elucidate the mechanisms underlying LUAD prognosis and promote individualized chemotherapy treatment.

Unraveling dynamic interactions between tumor-associated macrophages and consensus molecular subtypes in colorectal cancer: An integrative analysis of single-cell and bulk RNA transcriptome

Background

Accumulating research substantiated that tumor-associated macrophages (TAMs) have a significant impact on the tumorigenesis, progression, and distant metastasis, representing a novel target for various cancers. However, the underlying dynamic changes and interactions between TAMs and tumor cells remain largely elusive in colorectal cancer (CRC).

Methods

We depicted the dynamic changes of macrophages using sing-cell RNA-seq data and extracted TAM differentiation-related genes. Next, we utilized the weighted gene co-expression network analysis (WGCNA) to acquire CMS-related modular genes using bulk RNA-seq data. Finally, we utilized univariate Cox and Lasso Cox regression analyses to identify TAM differentiation-related biomarkers and established a novel risk signature model. We employed quantitative real-time polymerase chain reaction (qRT-PCR) on CRC tissue samples and used immunohistochemistry (IHC) data frome the HPA database to validate the mRNA and protein expression of prognostic genes. The interaction of TAMs and each consensus molecular subtype (CMS) subpopulation was analyzed at the cellular level.

Results

A total of 47,285 cells from single-cell dataset and 1197 CRC patients from bulk dataset were obtained. Among those, 6400 myeloid cells were re-clustered and annotated. RNASE1, F13A1, DAPK1, CLEC10A, RPN2, REG4 and RGS19 were identified as prognostic genes and the risk signature model was established based on the above genes. The qRT-PCR analysis indicated that the expression of RNASE1 and DAPK1 were significantly up-regulated in CRC tumor tissues. The cell-cell communication analysis demonstrated complex interactions between TAMs and CMS malignant cell subpopulations.

Conclusion

This study presents an in-depth dissection of the dynamic features of TAMs in the tumor microenvironment and provides promising therapeutic targets for CRC.

Intratumor microbiome derived glycolysis-lactate signatures depicts immune heterogeneity in lung adenocarcinoma by integration of microbiomic, transcriptomic, proteomic and single-cell data

Introduction

Microbiome plays roles in lung adenocarcinoma (LUAD) development and anti-tumor treatment efficacy. Aberrant glycolysis in tumor might promote lactate production that alter tumor microenvironment, affecting microbiome, cancer cells and immune cells. We aimed to construct intratumor microbiome score to predict prognosis of LUAD patients and thoroughly investigate glycolysis and lactate signature's association with LUAD immune cell infiltration.

Methods

The Cancer Genome Atlas-LUAD (TCGA-LUAD) microbiome data was downloaded from cBioPortal and analyzed to examine its association with overall survival to create a prognostic scoring model. Gene Set Enrichment Analysis (GSEA) was used to find each group's major mechanisms involved. Our study then investigated the glycolysis and lactate pattern in LUAD patients based on 19 genes, which were correlated with the tumor microenvironment (TME) phenotypes and immunotherapy outcomes. We developed a glycolysis-lactate risk score and signature to accurately predict TME phenotypes, prognosis, and response to immunotherapy.

Results

Using the univariate Cox regression analysis, the abundance of 38 genera were identified with prognostic values and a lung-resident microbial score (LMS) was then developed from the TCGA-LUAD-microbiome dataset. Glycolysis hallmark pathway was significantly enriched in high-LMS group and three distinct glycolysis-lactate patterns were generated. Patients in Cluster1 exhibited unfavorable outcomes and might be insensitive to immunotherapy. Glycolysis-lactate score was constructed for predicting prognosis with high accuracy and validated in external cohorts. Gene signature was developed and this signature was elevated in epithelial cells especially in tumor mass on single-cell level. Finally, we found that the glycolysis-lactate signature levels were consistent with the malignancy of histological subtypes.

Discussion

Our study demonstrated that an 18-microbe prognostic score and a 19-gene glycolysis-lactate signature for predicting prognosis of LUAD patients. Our LMS, glycolysis-lactate score and glycolysis-lactate signature have potential roles in precision therapy of LUAD patients.

Role of the membrane-spanning 4A gene family in lung adenocarcinoma

Lung adenocarcinoma, which is the second most prevalent cancer in the world, has a poor prognosis and a low 5-year survival rate. The MS4A protein family is crucial to disease development and progression, particularly for cancers, allergies, metabolic disorders, autoimmune diseases, infections, and neurodegenerative disorders. However, its involvement in lung adenocarcinoma remains unclear. In this study, we found that 11 MS4A family genes were upregulated or downregulated in lung adenocarcinoma. Furthermore, we described the genetic variation landscape of the MS4A family in lung adenocarcinoma. Notably, through functional enrichment analysis, we discovered that the MS4A family is involved in the immune response regulatory signaling pathway and the immune response regulatory cell surface receptor signaling pathway. According to the Kaplan-Meier curve, patients with lung adenocarcinoma having poor expression of MS4A2, MS4A7, MS4A14, and MS4A15 had a low overall survival rate. These four prognostic genes are substantially associated with immune-infiltrating cells, and a prognosis model incorporating them may more accurately predict the overall survival rate of patients with lung adenocarcinoma than current models. The findings of this study may offer creative suggestions and recommendations for the identification and management of lung adenocarcinoma.

HLA-DQA1 expression is associated with prognosis and predictable with radiomics in breast cancer

BACKGROUND

High HLA-DQA1 expression is associated with a better prognosis in many cancers. However, the association between HLA-DQA1 expression and prognosis of breast cancer and the noninvasive assessment of HLA-DQA1 expression are still unclear. This study aimed to reveal the association and investigate the potential of radiomics to predict HLA-DQA1 expression in breast cancer.

METHODS

In this retrospective study, transcriptome sequencing data, medical imaging data, clinical and follow-up data were downloaded from the TCIA ( https://www.cancerimagingarchive.net/ ) and TCGA ( https://portal.gdc.cancer.gov/ ) databases. The clinical characteristic differences between the high HLA-DQA1 expression group (HHD group) and the low HLA-DQA1 expression group were explored. Gene set enrichment analysis, Kaplan‒Meier survival analysis and Cox regression were performed. Then, 107 dynamic contrast-enhanced magnetic resonance imaging features were extracted, including size, shape and texture. Using recursive feature elimination and gradient boosting machine, a radiomics model was established to predict HLA-DQA1 expression. Receiver operating characteristic (ROC) curves, precision-recall curves, calibration curves, and decision curves were used for model evaluation.

RESULTS

The HHD group had better survival outcomes. The differentially expressed genes in the HHD group were significantly enriched in oxidative phosphorylation (OXPHOS) and estrogen response early and late signalling pathways. The radiomic score (RS) output from the model was associated with HLA-DQA1 expression. The area under the ROC curves (95% CI), accuracy, sensitivity, specificity, positive predictive value, and negative predictive value of the radiomic model were 0.866 (0.775-0.956), 0.825, 0.939, 0.7, 0.775, and 0.913 in the training set and 0.780 (0.629-0.931), 0.659, 0.81, 0.5, 0.63, and 0.714 in the validation set, respectively, showing a good prediction effect.

CONCLUSIONS

High HLA-DQA1 expression is associated with a better prognosis in breast cancer. Quantitative radiomics as a noninvasive imaging biomarker has potential value for predicting HLA-DQA1 expression.

Expression patterns of platinum resistance-related genes in lung adenocarcinoma and related clinical value models

The purpose of this study was to explore platinum resistance-related biomarkers and mechanisms in lung adenocarcinoma. Through the analysis of gene expression data of lung adenocarcinoma patients and normal patients from The Cancer Genome Atlas, Gene Expression Omnibus database, and A database of genes related to platinum resistance, platinum resistance genes in lung adenocarcinoma and platinum resistance-related differentially expressed genes were obtained. After screening by a statistical significance threshold, a total of 252 genes were defined as platinum resistance genes with significant differential expression, of which 161 were up-regulated and 91 were down-regulated. The enrichment results of up-regulated gene Gene Ontology (GO) showed that TOP3 entries related to biological processes (BP) were double-strand break repair, DNA recombination, DNA replication, the down-regulated gene GO enriches the TOP3 items about biological processes (BP) as a response to lipopolysaccharide, muscle cell proliferation, response to molecule of bacterial origin. Gene Set Enrichment Analysis showed that the top three were e2f targets, g2m checkpoint, and rgf beta signaling. A prognostic model based on non-negative matrix factorization classification showed the characteristics of high- and low-risk groups. The prognostic model established by least absolute shrinkage and selection operator regression and risk factor analysis showed that genes such as HOXB7, NT5E, and KRT18 were positively correlated with risk score. By analyzing the differences in m6A regulatory factors between high- and low-risk groups, it was found that FTO, GPM6A, METTL3, and YTHDC2 were higher in the low-risk group, while HNRNPA2B1, HNRNPC, TGF2BP1, IGF2BP2, IGF2BP3, and RBM15B were higher in the high-risk group. Immune infiltration and drug sensitivity analysis also showed the gene characteristics of the platinum-resistant population in lung adenocarcinoma. ceRNA analysis showed that has-miR-374a-5p and RP6-24A23.7 were lower in the tumor expression group, and that the survival of the low expression group was worse than that of the high expression group. In conclusion, the results of this study show that platinum resistance-related differentially expressed genes in lung adenocarcinoma are mainly concentrated in biological processes such as DNA recombination and response to lipopolysaccharide. The validation set proved that the high-risk group of our prognostic model had poor survival. M6A regulatory factor analysis, immune infiltration, and drug sensitivity analysis all showed differences between high and low-risk groups. ceRNA analysis showed that has-miR-374a-5p and RP6-24A23.7 could be protective factors. Further exploration of the potential impact of these genes on the risk and prognosis of drug-resistant patients with lung adenocarcinoma would provide theoretical support for future research.

The dynamic dysregulated network identifies stage-specific markers during lung adenocarcinoma malignant progression and metastasis

Brain metastasis occurs in approximately 30% of patients with lung adenocarcinoma (LUAD) and is closely associated with poor prognosis, recurrence, and death. However, dynamic gene regulation and molecular mechanism driving LUAD progression remain poorly understood. In this study, we performed a comprehensive single-cell transcriptome analysis using data from normal, early stage, advanced stage, and brain metastasis LUAD. Our single-cell-level analysis reveals the cellular composition heterogeneity at different stages during LUAD progression. We identified stage-specific risk genes that could contribute to LUAD progression and metastasis by reprogramming immune-related and metabolic-related functions. We constructed an early advanced metastatic dysregulated network and revealed the dynamic changes in gene regulations during LUAD progression. We identified 6 early advanced (HLA-DRB1, HLA-DQB1, SFTPB, SFTPC, PLA2G1B, and FOLR1), 8 advanced metastasis (RPS15, RPS11, RPL13A, RPS24, HLA-DRB5, LYPLA1, KCNJ15, and PSMA3), and 2 common risk genes in different stages (SFTPD and HLA-DRA) as prognostic markers in LUAD. Particularly, decreased expression of HLA-DRA, HLA-DRB1, HLA-DQB1, and HLA-DRB5 refer poor prognosis in LUAD by controlling antigen processing and presentation and T cell activation. Increased expression of PSMA3 and LYPLA1 refer poor prognosis by reprogramming fatty acid metabolism and RNA catabolic process. Our findings will help further understanding the pathobiology of brain metastases in LUAD.