Gene landscape and correlation between B-cell infiltration and programmed death ligand 1 expression in lung adenocarcinoma patients from The Cancer Genome Atlas data set

The genomic landscape of the immune system in lung cancer: present insights and continuing investigations

Lung cancer is one of the most prevalent malignancies worldwide, contributing to over a million cancer-related deaths annually. Despite extensive research investigating the genetic factors associated with lung cancer susceptibility and prognosis, few studies have explored genetic predispositions regarding the immune system. This review discusses the most recent genomic findings related to the susceptibility to or protection against lung cancer, patient survival, and therapeutic responses. The results demonstrated the effect of immunogenetic variations in immune system-related genes associated with innate and adaptive immune responses, cytokine, and chemokine secretions, and signaling pathways. These genetic diversities may affect the crosstalk between tumor and immune cells within the tumor microenvironment, influencing cancer progression, invasion, and prognosis. Given the considerable variability in the individual immunegenomics profiles, future studies should prioritize large-scale analyses to identify potential genetic variations associated with lung cancer using highthroughput technologies across different populations. This approach will provide further information for predicting response to targeted therapy and promotes the development of new measures for individualized cancer treatment.

Single-cell transcriptome analysis reveals intratumoral heterogeneity in lung adenocarcinoma

INTRODUCTION

Lung adenocarcinoma (LUAD) is a major health concern worldwide. Single-cell RNA-sequencing (scRNA-seq) provides a valuable platform for exploring the intratumoral heterogeneity in LUAD and holds great potential for facilitating the development and application of personalized therapeutic approaches.

METHODS

The TCGA-LUAD (n = 503), GSE68465 (n = 442), GSE72094 (n = 398), and GSE26939 (n = 115) datasets were retrieved for prognostic assessment. Subgroup analysis was performed for the epithelial cells, endothelial cells, immune cells, and fibroblasts, and the transcription factors and tumor-related pathways enriched in each subgroup were analyzed using PROGENy and DoRothEA package. The InferCNV software was used to calculate the copy number variations (CNVs) in tumor cell subgroups with normal epithelial cells as the reference. The association between the annotated cell types and survival was analyzed using the Scissor software.

RESULTS

We identified eight major cell types in LUAD, namely epithelial cells, NK cells, T and B cells, endothelial cells, mast cells, myeloid cells, and fibroblasts, of which the epithelial cells and B cells showed a marked increase in the tumor samples. In addition, we also detected an intense signal transduction network from the cancer-associated fibroblasts (CAFs) to malignant cells, mainly involving the DCN/MET, COLA1/DDR1, COL1A1/SDC1, and COL1A2/SDC1 pathways. The tumor differentiation trajectory consisted of state 1 and state 2, which were enriched in HIF1A, and state 4. Furthermore, only a few B cells originated from the normal tissue, suggesting significant recruitment and infiltration of B cells in LUAD. Based on differentially upregulated genes in the cells positively and negatively associated with survival, we established a prognostic model that showed satisfactory predictive performance in three different cohorts. States 3 and 2 of epithelial cells included the majority of cells with KRAS mutation, whereas state 2 showed high frequency of EGFR mutations.

CONCLUSION

We analyzed intra-tumor heterogeneity of LUAD at the single-cell level and developed a prognostic index that was highly effective across multiple cohorts.

Neglected no more: B cell-mediated anti-tumor immunity

Immuno-oncology has traditionally focused on the cellular arm of the adaptive immune response, while attributing tumor-promoting activity to humoral responses in tumor-bearing hosts. This view stems from mouse models that do not necessarily recapitulate the antibody response process consistently observed in most human cancers. In recent years, the field has reconsidered the coordinated action of T and B cell responses in the context of anti-tumor immunity, as in any other immune response. Thus, recent studies in human cancer identify B cell responses with better outcome, typically in association with superior T cell responses. An area of particular interest is tertiary lymphoid structures, where germinal centers produce isotype switched antibodies and B cells and T lymphocytes interact with other immune cell types. The presence of these lymphoid structures is associated with better immunotherapeutic responses and remain poorly understood. Here, we discuss recent discoveries on how coordination between humoral and cellular responses is required for effective immune pressure against malignant progression, providing a perspective on the role of tertiary lymphoid structures and interventions to elicit their formation in unresectable tumors.

LINC02126 is a potential diagnostic, prognostic and immunotherapeutic target for lung adenocarcinoma

Background

Adenocarcinoma has long been an independent histological class of lung cancer, which leads to high morbidity and mortality. We aimed to investigate the contribution of LINC02126 in lung adenocarcinoma.

Methods

RNA sequencing data and clinical information were downloaded. Diagnostic efficiency and survival analysis of LINC02126 were performed, followed by functional analysis of genes co-expressed with LINC02126 and differentially expressed genes (DEGs) in different LINC02126 expression groups. Tumor immune microenvironment (TIME) cell infiltration and correlation analysis of tumor mutation burden were performed in different LINC02126 expression groups.

Results

In lung adenocarcinoma, the expression level of LINC02126 was significantly decreased. Significant expression differences of LINC02126 were found in some clinical variables, including T staging, M staging, sex, stage, and EGFR mutation. LINC02126 had potential diagnostic and prognostic value for patients. In the low LINC02126 expression group, the infiltration degree of most immune cells was significantly lower than that in the high LINC02126 expression group. Tumor mutation burden level and frequency of somatic mutation in patients with low LINC02126 expression group were significantly higher than in patients with high LINC02126 expression group.

Conclusions

LINC02126 could be considered as a diagnostic, prognostic and immunotherapeutic target for lung adenocarcinoma.

B cell deficiency promotes the initiation and progression of lung cancer

Currently commercialized CAR-T cell therapies targeting CD19 and BCMA show great efficacy to cure B cell malignancies. However, intravenous infusion of these CAR-T cells severely destroys both transformed and normal B cells in most tissues and organs, in particular lung, leading to a critical question that what the impact of normal B cell depletion on pulmonary diseases and lung cancer is. Herein, we find that B cell frequency is remarkably reduced in both smoking carcinogen-treated lung tissues and lung tumors, which is associated with advanced cancer progression and worse patient survival. B cell depletion by anti-CD20 antibody significantly accelerates the initiation and progression of lung tumors, which is mediated by repressed tumor infiltration of T cells and macrophage elimination of tumor cells. These findings unveil the overall antitumor activity of B cells in lung cancer, providing novel insights into both mechanisms underlying lung cancer pathogenesis and clinical prevention post CAR-T cell therapy.

Upregulation of EMID1 Accelerates to a Favorable Prognosis and Immune Infiltration in Lung Adenocarcinoma

Lung cancer is a difficult-to-treat cancer. Lung adenocarcinoma (LUAD) is the main subtype of lung cancer. Although there are many ways to treat lung cancer, the survival rate of patients is low. Therefore, novel molecules need to be identified to diagnose and treat LUAD. This study utilized The Cancer Genome Atlas (TCGA) LUAD data to analyze and validate the value of EMID1 as a LUAD diagnostic surface marker and overall survival prognostic marker. Differential expression analysis formally confirmed that decreased EMID1 expression was significantly associated with advanced stage and metastasis of lung cancer. Kaplan-Meier survival analysis showed that the patients with low EMID expression are dismal. The relationship between clinicopathological features and EMID1 was scored using Wilcoxon signed-rank test and R (v.3.5.1) logistic regression and suggested that patients with low EMID1 expression had a worse prognosis than patients with high EMID1 expression. (Gene Ontology) GO, Kyoto Encyclopedia of Genes and Genomes(KEGG), and gene set enrichment analysis (GSEA) were performed to investigate the potential mechanism of EMID1 expression on the prognosis of LUAD and suggested that Notch signaling pathway may be an important biological pathway for EMID1 to play a role in LUAD. Further, combined with univariate and multivariate Cox regression analysis, it was speculated that high and low levels of EMID1 expression and the logistic regression analysis of related clinical variables had significant clinical significance to verify the underlying mechanism of LUAD focus and prognosis. EMID1 plays an important role in the immune milieu of LUAD. Meanwhile, the correlation between tumor-infiltrating immune cells and genes was assessed using CIBERSORT, and it was found that the level of B cell infiltration was positively correlated with the expression of EMID1, all of which were validated in the GEO and GEPIA databases. In all, this study helps to understand the immune microenvironment of LUAD and improve the survival of patients with LUAD. Thus, EMID1 may be a novel immune-related prognostic marker of LUAD.

Integrated analysis of single-cell and bulk RNA-sequencing identifies a signature based on B cell marker genes to predict prognosis and immunotherapy response in lung adenocarcinoma

As an essential component of the tumor microenvironment, B cells exist in all stages of tumor and exert important roles in anti-tumor immunity and shaping tumor development. We aimed to explore the expression profile of B cell marker genes and construct a prognostic signature based on these genes in Lung adenocarcinoma (LUAD). A total of 1268 LUAD patients from different cohorts were enrolled in this study. We performed an analysis of single-cell RNA-sequencing (scRNA-seq) data from Gene expression omnibus (GEO) database to identify B cell marker genes in LUAD. TCGA database was used to construct signature, and six cohorts from GEO database were used for validation. We also investigated the association between this signature and immunotherapy response. Based on 258 B cell marker genes identified by scRNA-seq analysis, a nine-gene signature was constructed for prognostic prediction in TCGA dataset, which classified patients into high-risk and low-risk groups according to overall survival. The multivariate analysis demonstrated that the signature was an independent prognostic factor. The signature's predictive power was verified in other six independent cohorts and different clinical subgroups. Analysis of immune profiles showed that high-risk groups presented discriminative immune-cell infiltrations and immune-suppressive states. More importantly, risk scores of the signature were closely correlated with PD-L1, tumor mutation burden, neoantigens, and tumor immune dysfunction and exclusion score. Our study proposed a novel prognostic signature based on B cell marker genes for LUAD patients. The signature could effectively indicate LUAD patients' survival and serve as a predictor for immunotherapy.

Exploration of immune-related cells and ceRNA in squamous cell lung cancer

The treatment for squamous cell lung cancer (SqCLC) is limited, and the prognosis of SqCLC is poor. In this article, we aimed to analyze and identify immune-related cells and competition endogenous RNA (ceRNA) that influence the prognosis of SqCLC. SqCLC and lung adenocarcinoma data were downloaded from TCGA-GDC. A total of 22 types of immune cell fractions were estimated using CIBERSORT. R software was used to identify any significantly different transcriptome data, including mRNA, LncRNA, and miRNA. The univariate cox regression method was applied to screen for prognosis-related lncRNA, miRNA, mRNA and tumor-infiltrating immune cells. There were 504 patients included in this study. There was a higher proportion of memory activated CD4+ T cells and CD8+ T cells in younger women. Follicular helper T (Tfh) cells were predictive of a good prognosis and reflected immune activation in SqCLC. The SFTA1P/NKX2-1-AS1, hsa-mir-503, GREM2 ceRNA axes and NKX2-1-AS1, hsa-mir-96, PROK2 ceRNA axes were found to be important for the immune function, pathogenesis, and prognosis of SqCLC. Collectively, the immune-related ceRNA and tumor-infiltrating immune cells in SqCLC are likely important determinants of SqCLC pathogenesis, prognosis, and immune status.

B cells in lung cancer-not just a bystander cell: a literature review

Metastatic lung cancer represents a significant global issue where it is responsible for the most cancer diagnoses and deaths worldwide. Treatment for advanced lung cancer has undergone a series of paradigm shifts from chemotherapy to targeted molecular agents to the most recent immunotherapy strategies. The most successful of the latter involves antibodies that block inhibitory receptors on tumor infiltrating T cells, thereby enhancing T cell activity against tumor cells. However, only a subset of patients demonstrate durable responses to these drugs and treatment resistance is common. Emerging evidence suggests that a critical role exists for B cells as more than a bystander immune cell in the tumor microenvironment (TME). However, this role is likely context-specific where B cells comprise distinct subtypes with unique effector functions that may result in anti- or pro-tumor effects. As such, the balance between various B cell subtypes affects the net B cell impact upon tumor immunity. To date, the factors needed to polarize B cell function toward anti-tumor activity are unclear. Understanding B cell biology in the lung cancer setting will help redefine and refine treatment strategies to augment anti-tumor immunity. This article presents a review of the literature describing the current knowledge of the development and function of B cells, and explores their role in lung cancer and potential as an immunotherapeutic strategy and as a predictive marker for response to immune checkpoint blockade.

Tumor infiltrated immune cell types support distinct immune checkpoint inhibitor outcomes in patients with advanced non-small cell lung cancer

The evaluation of PD‐L1 expression alone has limitations in predicting clinical outcome in immune‐checkpoint inhibitors (ICI). This study aimed to evaluate the predictive and prognostic effects of the presence of various immune cells in pretreatment tissue samples and to identify determinants associated with response in patients with advanced non‐small cell lung cancer (NSCLC) treated with PD‐1 blockade. Immune cell distribution was heterogeneous and the most dominant immune cell type was T cells. Patients with durable clinical benefit (DCB) showed significantly higher PD‐L1 expression. The ratio of tumor/stroma region of T cell, B cell, and macrophage was significantly higher in patient with DCB. High intratumoral T‐ and B‐cell density (≥median) was associated with DCB in the low PD‐L1 expression (<50%) group. In univariate analyses, the overall survival (OS) benefit was shown according to intratumoral B‐cell density (p = 0.0337). The incidence of hyperprogressive disease (HPD) was 13.0%. The Chi‐square test revealed that HPD was significantly associated with intratumoral B‐cell density but not T‐cell or macrophage density. Our results demonstrate different predictive and prognostic values for infiltrating immune cells in tumor tissue, which may help in selecting patients for ICI.

Development of an immune-related gene pairs signature for predicting clinical outcome in lung adenocarcinoma

Lung adenocarcinoma (LUAD) is the main pathological subtype of Non-small cell lung cancer. We downloaded the gene expression profile and immune-related gene set from the TCGA and ImmPort database, respectively, to establish immune-related gene pairs (IRGPs). Then, IRGPs were subjected to univariate Cox regression analysis, LASSO regression analysis, and multivariable Cox regression analysis to screen and develop an IRGPs signature. The receiver operating characteristic curve (ROC) was applied for evaluating the predicting accuracy of this signature by calculating the area under ROC (AUC) and data from the GEO set was used to validate this signature. The relationship of 22 tumor-infiltrating immune cells (TIICs) to the immune risk score was also investigated. An IRGPs signature with 8 IRGPs was constructed. The AUC for 1- and 3-year overall survival in the TCGA set was 0.867 and 0.870, respectively. Similar results were observed in the AUCs of GEO set 1, 2 and 3 (GEO set 1 [1-year: 0.819; 3-year: 0.803]; GEO set 2 [1-year: 0.834; 3-year: 0.870]; GEO set 3 [1-year: 0.955; 3-year: 0.827]). Survival analysis demonstrated high-risk LUAD patients exhibited poorer prognosis. The multivariable Cox regression indicated that the risk score was an independent prognostic factor. The immune risk score was highly associated with several TIICs (Plasma cells, memory B cells, resting memory CD4 T cells, and activated NK cells). We developed a novel IRGPs signature for predicting 1- and 3- year overall survival in LUAD, which would be helpful for prognosis assessment of LUAD.

Heterogeneous Tumor-Immune Microenvironments between Primary and Metastatic Tumors in a Patient with ALK Rearrangement-Positive Large Cell Neuroendocrine Carcinoma

Evolution of tumor-immune microenviroments (TIMEs) occurs during tumor growth and dissemination. Understanding inter-site tumor-immune heterogeneity is essential to harness the immune system for cancer therapy. While the development of immunotherapy against lung cancer with driver mutations and neuroendocrine tumors is ongoing, little is known about the TIME of large cell neuroendocrine carcinoma (LCNEC) or anaplastic lymphoma kinase (ALK) rearrangement-positive lung cancer. We present a case study of a 32-year-old female patient with ALK-rearrangement-positive LCNEC, who had multiple distant metastases including mediastinal lymph-node, bilateral breasts, multiple bones, liver and brain. Multiple biopsy samples obtained from primary lung and three metastatic tumors were analyzed by fluorescent multiplex immunohistochemistry. Tissue localizations of tumor-infiltrating lymphocytes in the tumor nest and surrounding stroma were evaluated. T cell and B cell infiltrations were decreased with distance from primary lung lesion. Although each tumor displayed a unique TIME, all tumors exhibited concomitant regression after treatment with an ALK-inhibitor. This study provides the first evidence of the coexistence of distinct TIME within a single individual with ALK-rearrangement-positive LCNEC. The present study contributes to our understanding of heterogeneous TIMEs between primary and metastatic lesions and provides new insights into the complex interplay between host-immunity and cancer cells in primary and metastatic lesions.

Duality of B Cell-CXCL13 Axis in Tumor Immunology

Tumor immunity is a rapidly evolving area of research consisting of many possible permutations of immune cell tumor interactions that are dependent upon cell type, tumor type, and stage in tumor progression. At the same time, the majority of cancer immunotherapies have been focused on modulating the T cell-mediated antitumor immune response and have largely ignored the potential utility that B cells possess with respect to tumor immunity. Therefore, this motivated an exploration into the role that B cells and their accompanying chemokine, CXCL13, play in tumor immunity across multiple tumor types. Both B cells and CXCL13 possess dualistic impacts on tumor progression and tumor immunity which is furthered detail in this review. Specifically, various B cells subtypes are able to suppress or enhance several important immunological functions. Paradoxically, CXCL13 has been shown to drive several pro-growth and invasive signaling pathways across multiple tumor types, while also, correlating with improved survival and immune cell tumor localization in other tumor types. Potential tools for better elucidating the mechanisms by which B cells and CXCL13 impact the antitumor immune response are also discussed. In addition, multiples strategies are proposed for modulating the B cell-CXCL13 axis for cancer immunotherapies.

Comprehensive analyses of the heterogeneity and prognostic significance of tumor-infiltrating immune cells in non-small-cell lung cancer: Development and validation of an individualized prognostic model

Understanding the role of tumor-infiltrating immune cells (TIICs) in non-small cell lung cancer (NSCLC) is critical to finding new prognostic biomarkers and improving prognostic evaluation. Herein, we aimed to comprehensively analyze tumor-infiltrating pattern of TIICs in NSCLC and build a TIICs-associated, risk-stratification prognostic model for clinical practice. We applied CIBERSORT and ESTIMATE computational methods to analyze RNA-seq samples of 852 NSCLC patients from The Cancer Genome Atlas (TCGA). Prognotic factors were identified by univariate and multivariate Cox regression analyses for overall survival (OS). A novel model was developed to predict the 1-, 3- and 5-year OS of NSCLC based on the TCGA cohort, validated by external validation cohorts (GSE31210, GSE37745), and then evaluated by C-indexes and calibration plots. Significant heterogeneity in the infiltrating patterns of TIICs was shown among various pathological subtypes of NSCLC and between different genders. Further analyses showed that abundances of naive B cells (NBCs), T cells and mast cells (MCs) were positively correlated with prognosis. Tumor samples with high T cells abundances tended to have higher expression levels of immune checkpoint genes (PD-1, PD-L1, CTLA-4). A new immune-gene related index (IGRI) was built by five immune-related differentially expressed genes (DEGs) including BTK, CCR2, CLEC10A, NCR3 and PRKCB, which were closely correlated with TIICs abundances and prognosis. Tumor stage, IGRI, abundances of NBCs, T cells, MCs and NK cells were significant independent prognostic factors and were included in the nomogram as predictors. The internal and external calibration plots of the nomogram were in excellent agreement. This study reveals that TIICs are significantly correlated with clinicopathological features and prognosis in NSCLC and thus can be potential prognostic biomarker or therapeutic target. The remarkable heterogeneity of TIICs suggests that specific infiltrating patterns of TIICs should also be taken into consideration when determining individualized immunotherapy strategies for NSCLC patients.

Humoral immune responses: Unsung heroes of the war on cancer

Solid cancers progress from primordial lesions through complex interactions between tumor-promoting and anti-tumor immune cell types, ultimately leading to the orchestration of humoral and T cell adaptive immune responses, albeit in an immunosuppressive environment. B cells infiltrating most established tumors have been associated with a dual role: Some studies have associated antibodies produced by tumor-associated B cells with the promotion of regulatory activities on myeloid cells, and also with direct immunosuppression through the production of IL-10, IL-35 or TGF-β. In contrast, recent studies in multiple human malignancies identify B cell responses with delayed malignant progression and coordinated T cell protective responses. This includes the elusive role of Tertiary Lymphoid Structures identified in many human tumors, where the function of B cells remains unknown. Here, we discuss emerging data on the dual role of B cell responses in the pathophysiology of human cancer, providing a perspective on future directions and possible novel interventions to restore the coordinated action of both branches of the adaptive immune response, with the goal of maximizing immunotherapeutic effectiveness.

The Role of Tumor-Infiltrating B Cells in Tumor Immunity

Earlier studies on elucidating the role of lymphocytes in tumor immunity predominantly focused on T cells. However, the role of B cells in tumor immunity has increasingly received better attention in recent studies. The B cells that infiltrate tumor tissues are called tumor-infiltrating B cells (TIBs). It is found that TIBs play a multifaceted dual role in regulating tumor immunity rather than just tumor inhibition or promotion. In this article, latest research advances focusing on the relationship between TIBs and tumor complexity are reviewed, and light is shed on some novel ideas for exploiting TIBs as a possible tumor biomarker and potential therapeutic target against tumors.

The Risks and Benefits of Immune Checkpoint Blockade in Anti-AChR Antibody-Seropositive Non-Small Cell Lung Cancer Patients

Background: Anti-programmed cell death 1 (PD-1) monoclonal antibodies (Abs) unleash an immune response to cancer. However, a disruption of the immune checkpoint function by blocking PD-1/PD-ligand 1(PD-L1) signaling may trigger myasthenia gravis (MG) as a life-threatening immune-related adverse event. MG is a neuromuscular disease and is closely associated with being positive for anti-acetylcholine receptor (anti-AChR) Abs, which are high specific and diagnostic Abs for MG. Methods: A 72-year-old man was diagnosed with chemotherapy-refractory lung squamous cell carcinoma and nivolumab was selected as the third-line regimen. We describe the first report of an anti-AChR Ab-seropositive lung cancer patient achieving a durable complete response (CR) to an anti-PD-1 antibody therapy. To further explore this case, we performed multiplex immunofluorescence analysis on a pretreatment tumor. Results: The patient achieved a durable CR without developing MG. However, the levels of anti-AChR Abs were elevated during two years of anti-PD-1 antibody therapy. The tumor of the subclinical MG patient had high PD-L1 expression and an infiltrated–inflamed tumor immune microenvironment. Conclusions: This study suggests that immune checkpoint inhibitors can be safely used and provide the benefits for advanced cancer patients with immunologically ‘hot’ tumor even if anti-AChR Abs are positive. Although careful monitoring clinical manifestation in consultation with neurologist is needed, immune checkpoint inhibitors should be considered as a treatment option for asymptomatic anti-AChR Ab-seropositive cancer patients.