An atlas of epithelial cell states and plasticity in lung adenocarcinoma.

Understanding the cellular processes that underlie early lung adenocarcinoma (LUAD) development is needed to devise intervention strategies1. Here we studied 246,102 single epithelial cells from 16 early-stage LUADs and 47 matched normal lung samples. Epithelial cells comprised diverse normal and cancer cell states, and diversity among cancer cells was strongly linked to LUAD-specific oncogenic drivers. KRAS mutant cancer cells showed distinct transcriptional features, reduced differentiation and low levels of aneuploidy. Non-malignant areas surrounding human LUAD samples were enriched with alveolar intermediate cells that displayed elevated KRT8 expression (termed KRT8+ alveolar intermediate cells (KACs) here), reduced differentiation, increased plasticity and driver KRAS mutations. Expression profiles of KACs were enriched in lung precancer cells and in LUAD cells and signified poor survival. In mice exposed to tobacco carcinogen, KACs emerged before lung tumours and persisted for months after cessation of carcinogen exposure. Moreover, they acquired Kras mutations and conveyed sensitivity to targeted KRAS inhibition in KAC-enriched organoids derived from alveolar type 2 (AT2) cells. Last, lineage-labelling of AT2 cells or KRT8+ cells following carcinogen exposure showed that KACs are possible intermediates in AT2-to-tumour cell transformation. This study provides new insights into epithelial cell states at the root of LUAD development, and such states could harbour potential targets for prevention or intervention.

The microRNA-183/96/182 cluster inhibits lung cancer progression and metastasis by inducing an interleukin-2-mediated antitumor CD8+ cytotoxic T-cell response.

One of the mechanisms by which cancer cells acquire hyperinvasive and migratory properties with progressive loss of epithelial markers is the epithelial-to-mesenchymal transition (EMT). We have previously reported that in different cancer types, including nonsmall cell lung cancer (NSCLC), the microRNA-183/96/182 cluster (m96cl) is highly repressed in cells that have undergone EMT. In the present study, we used a novel conditional m96cl mouse to establish that loss of m96cl accelerated the growth of Kras mutant autochthonous lung adenocarcinomas. In contrast, ectopic expression of the m96cl in NSCLC cells results in a robust suppression of migration and invasion in vitro, and tumor growth and metastasis in vivo. Detailed immune profiling of the tumors revealed a significant enrichment of activated CD8+ cytotoxic T lymphocytes (CD8+ CTLs) in m96cl-expressing tumors, and m96cl-mediated suppression of tumor growth and metastasis was CD8+ CTL-dependent. Using coculture assays with naïve immune cells, we show that m96cl expression drives paracrine stimulation of CD8+ CTL proliferation and function. Using tumor microenvironment-associated gene expression profiling, we identified that m96cl elevates the interleukin-2 (IL2) signaling pathway and results in increased IL2-mediated paracrine stimulation of CD8+ CTLs. Furthermore, we identified that the m96cl modulates the expression of IL2 in cancer cells by regulating the expression of transcriptional repressors Foxf2 and Zeb1, and thereby alters the levels of secreted IL2 in the tumor microenvironment. Last, we show that in vivo depletion of IL2 abrogates m96cl-mediated activation of CD8+ CTLs and results in loss of metastatic suppression. Therefore, we have identified a novel mechanistic role of the m96cl in the suppression of lung cancer growth and metastasis by inducing an IL2-mediated systemic CD8+ CTL immune response.

Mechanical regulation of glycolysis via cytoskeleton architecture.

The mechanics of the cellular microenvironment continuously modulates cell functions such as growth, survival, apoptosis, differentiation and morphogenesis via cytoskeletal remodelling and actomyosin contractility1-3. Although all of these processes consume energy4,5, it is unknown whether and how cells adapt their metabolic activity to variable mechanical cues. Here we report that the transfer of human bronchial epithelial cells from stiff to soft substrates causes a downregulation of glycolysis via proteasomal degradation of the rate-limiting metabolic enzyme phosphofructokinase (PFK). PFK degradation is triggered by the disassembly of stress fibres, which releases the PFK-targeting E3 ubiquitin ligase tripartite motif (TRIM)-containing protein 21 (TRIM21). Transformed non-small-cell lung cancer cells, which maintain high glycolytic rates regardless of changing environmental mechanics, retain PFK expression by downregulating TRIM21, and by sequestering residual TRIM21 on a stress-fibre subset that is insensitive to substrate stiffness. Our data reveal a mechanism by which glycolysis responds to architectural features of the actomyosin cytoskeleton, thus coupling cell metabolism to the mechanical properties of the surrounding tissue. These processes enable normal cells to tune energy production in variable microenvironments, whereas the resistance of the cytoskeleton in response to mechanical cues enables the persistence of high glycolytic rates in cancer cells despite constant alterations of the tumour tissue.

EMT activates exocytotic Rabs to coordinate invasion and immunosuppression in lung cancer.

Epithelial-to-mesenchymal transition (EMT) underlies immunosuppression, drug resistance, and metastasis in epithelial malignancies. However, the way in which EMT orchestrates disparate biological processes remains unclear. Here, we identify an EMT-activated vesicular trafficking network that coordinates promigratory focal adhesion dynamics with an immunosuppressive secretory program in lung adenocarcinoma (LUAD). The EMT-activating transcription factor ZEB1 drives exocytotic vesicular trafficking by relieving Rab6A, Rab8A, and guanine nucleotide exchange factors from miR-148a-dependent silencing, thereby facilitating MMP14-dependent focal adhesion turnover in LUAD cells and autotaxin-mediated CD8+ T cell exhaustion, indicating that cell-intrinsic and extrinsic processes are linked through a microRNA that coordinates vesicular trafficking networks. Blockade of ZEB1-dependent secretion reactivates antitumor immunity and negates resistance to PD-L1 immune checkpoint blockade, an important clinical problem in LUAD. Thus, EMT activates exocytotic Rabs to drive a secretory program that promotes invasion and immunosuppression in LUAD.

Analysis of Immune Intratumor Heterogeneity Highlights Immunoregulatory and Coinhibitory Lymphocytes as Hallmarks of Recurrence in Stage I Non-Small Cell Lung Cancer.

Our understanding of the molecular mechanisms underlying postsurgical recurrence of non-small cell lung cancer (NSCLC) is rudimentary. Molecular and T cell repertoire intratumor heterogeneity (ITH) have been reported to be associated with postsurgical relapse; however, how ITH at the cellular level impacts survival is largely unknown. Here we report the analysis of 2880 multispectral images representing 14.2% to 27% of tumor areas from 33 patients with stage I NSCLC, including 17 cases (relapsed within 3 years after surgery) and 16 controls (without recurrence ≥5 years after surgery) using multiplex immunofluorescence. Spatial analysis was conducted to quantify the minimum distance between different cell types and immune cell infiltration around malignant cells. Immune ITH was defined as the variance of immune cells from 3 intratumor regions. We found that tumors from patients having relapsed display different immune biology compared with nonrecurrent tumors, with a higher percentage of tumor cells and macrophages expressing PD-L1 (P =.031 and P =.024, respectively), along with an increase in regulatory T cells (Treg) (P =.018), antigen-experienced T cells (P =.025), and effector-memory T cells (P =.041). Spatial analysis revealed that a higher level of infiltration of PD-L1+ macrophages (CD68+PD-L1+) or antigen-experienced cytotoxic T cells (CD3+CD8+PD-1+) in the tumor was associated with poor overall survival (P =.021 and P =.006, respectively). A higher degree of Treg ITH was associated with inferior recurrence-free survival regardless of tumor mutational burden (P =.022), neoantigen burden (P =.021), genomic ITH (P =.012) and T cell repertoire ITH (P =.001). Using multiregion multiplex immunofluorescence, we characterized ITH at the immune cell level along with whole exome and T cell repertoire sequencing from the same tumor regions. This approach highlights the role of immunoregulatory and coinhibitory signals as well as their spatial distribution and ITH that define the hallmarks of tumor relapse of stage I NSCLC.

CD30 expression is frequently decreased in relapsed classic Hodgkin lymphoma after anti-CD30 CAR T-cell therapy.

Chimeric antigen receptor (CAR) T-cells anti-CD30 is an innovative therapeutic option that has been used to treat cases of refractory/relapsed (R/R) classic Hodgkin lymphoma (CHL). Limited data are available regarding the CD30 expression status of patients who relapsed after this therapy. This is the first study to show decreased CD30 expression in R/R CHL in patients (n = 5) who underwent CAR T-cell therapy in our institution between 2018 and 2022. Although conventional immunohistochemical assays showed decreased CD30 expression in neoplastic cells in all cases (8/8) the tyramide amplification assay and RNAScope in situ hybridisation detected CD30 expression at different levels in 100% (n = 8/8) and 75% (n = 3/4), respectively. Hence, our findings document that certain levels of CD30 expression are retained by the neoplastic cells. This is not only of biological interest but also diagnostically important, as detection of CD30 is an essential factor in establishing a diagnosis of CHL.

Gene expression profiling and multiplex immunofluorescence analysis of bullous pemphigoid immune-related adverse event reveal upregulation of toll-like receptor 4/complement-induced innate immune response and increased density of TH 1 T-cells.

BACKGROUND: Immune checkpoint inhibitor (ICI)-based cancer therapies cause a variety of cutaneous immune-related adverse events (irAEs) including immunobullous skin eruptions like bullous pemphigoid (BP). However, little is known about the underlying immunopathogenic drivers of these reactions, and understanding the unique gene expression profile and immune composition of BP-irAE remains a critical knowledge gap in the field of oncodermatology/oncodermatopathology. METHODS: BP-irAE (n = 8) and de novo BP control (n = 8) biopsy samples were subjected to gene expression profiling using the NanoString® Technologies nCounter PanCancer Immune Profiling Panel. Multiplex immunofluorescence (mIF) studies using markers for T-cells (CD3 and CD8), T helper 1 (TH 1) cells (Tbet), TH 2 cells (Gata3), TH 17 cells (RORγT), and regulatory T-cells (Tregs; FoxP3) were further evaluated using InForm® image analysis. RESULTS: Compared with de novo BP controls, BP-irAE samples exhibited upregulation of 30 mRNA transcripts (p < 0.025), including toll-like receptor 4 (TLR4) and genes associated with complement activation, and downregulation of 89 mRNA transcripts (p < 0.025), including genes associated with TH 2, TH 17, and B-cell immune response. BP-irAE demonstrated a greater density of Tbet+ (TH 1) cells in the dermis (p = 0.004) and fewer Tregs in the blister floor (p = 0.028) when compared with that of de novo control BP samples. CONCLUSIONS: BP-irAE exhibited activation of the TLR4/complement-driven classical innate immune response pathway, with dermal TH 1 immune cell polarization and decreased Tregs in the blister floor. TLR/complement signaling may underlie the immunopathogenesis of BP-irAE.

Effector memory cytotoxic CD3+/CD8+/CD45RO+ T cells are predictive of good survival and a lower risk of recurrence in triple-negative breast cancer.

Triple-negative breast cancer (TNBC) with high tumour-infiltrating lymphocytes (TILs) has been associated with a promising prognosis. To better understand the prognostic value of immune cell subtypes in TNBC, we characterised TILs and the interaction between tumour cells and immune cell subtypes. A total of 145 breast cancer tissues were stained by multiplex immunofluorescence (mIF), including panel 1 (PD-L1, PD-1, CD3, CD8, CD68 and CK) and panel 2 (Foxp3, Granzyme B, CD45RO, CD3, CD8 and CK). Phenotypes were analysed and quantified by pathologists using InForm software. We found that in the ER-negative (ER <1% and HER2-negative) group and the ER/PR-low positive (ER 1-9% and HER2-negative) group, 11.2% and 7.1% of patients were PD-L1+ by the tumour cell score, 29.0% and 28.6% were PD-L1+ by the modified immune cell score and 30.8% and 32.1% were PD-L1+ by the combined positive score. We combined ER-negative and ER/PR-low positive cases for the survival analysis since a 10% cut-off is often used in clinical practice for therapeutic purposes. The densities of PD-L1+ tumour cells (HR: 0.366, 95% CI: 0.138-0.970; p = 0.043) within the tumour compartment and CD3+ immune cells in the total area (tumour and stromal compartments combined) (HR: 0.213, 95% CI: 0.070-0.642; p = 0.006) were favourable prognostic biomarkers for overall survival (OS) in TNBC. The density of effector/memory cytotoxic T cells (CD3+CD8+CD45RO+) in the tumour compartment was an independent prognostic biomarker for OS (HR: 0.232, 95% CI: 0.086-0.628; p = 0.004) and DFS (HR: 0.183, 95% CI: 0.1301-0.744; p = 0.009) in TNBC. Interestingly, spatial data suggested that patients with a higher density of PD-L1+ tumour cells had shorter cell-cell distances from tumour cells to cytotoxic T cells (p < 0.01). In conclusion, we found that phenotyping tumour immune cells by mIF is highly informative in understanding the immune microenvironment in TNBC. PD-L1+ tumour cells, total T cells and effector/memory cytotoxic T cells are promising prognostic biomarkers in TNBC.

Augmented Lipocalin-2 Is Associated with Chronic Obstructive Pulmonary Disease and Counteracts Lung Adenocarcinoma Development.

Rationale: Early pathogenesis of lung adenocarcinoma (LUAD) remains largely unknown. We found that, relative to wild-type littermates, the innate immunomodulator Lcn2 (lipocalin-2) was increased in normal airways from mice with knockout of the airway lineage gene Gprc5a (Gprc5a-/-) and that are prone to developing inflammation and LUAD. Yet, the role of LCN2 in lung inflammation and LUAD is poorly understood.Objectives: Delineate the role of Lcn2 induction in LUAD pathogenesis.Methods: Normal airway brushings, uninvolved lung tissues, and tumors from Gprc5a-/- mice before and after tobacco carcinogen exposure were analyzed by RNA sequencing. LCN2 mRNA was analyzed in public and in-house data sets of LUAD, lung squamous cancer (LUSC), chronic obstructive pulmonary disease (COPD), and LUAD/LUSC with COPD. LCN2 protein was immunohistochemically analyzed in a tissue microarray of 510 tumors. Temporal lung tumor development, gene expression programs, and host immune responses were compared between Gprc5a-/- and Gprc5a-/-/Lcn2-/- littermates.Measurements and Main Results:Lcn2 was progressively elevated during LUAD development and positively correlated with proinflammatory cytokines and inflammation gene sets. LCN2 was distinctively elevated in human LUADs, but not in LUSCs, relative to normal lungs and was associated with COPD among smokers and patients with LUAD. Relative to Gprc5a-/- mice, Gprc5a-/-/Lcn2-/- littermates exhibited significantly increased lung tumor development concomitant with reduced T-cell abundance (CD4+) and richness, attenuated antitumor immune gene programs, and increased immune cell expression of protumor inflammatory cytokines.Conclusions: Augmented LCN2 expression is a molecular feature of COPD-associated LUAD and counteracts LUAD development in vivo by maintaining antitumor immunity.

Naproxen chemoprevention promotes immune activation in Lynch syndrome colorectal mucosa.

OBJECTIVE: Patients with Lynch syndrome (LS) are at markedly increased risk for colorectal cancer. It is being increasingly recognised that the immune system plays an essential role in LS tumour development, thus making an ideal target for cancer prevention. Our objective was to evaluate the safety, assess the activity and discover novel molecular pathways involved in the activity of naproxen as primary and secondary chemoprevention in patients with LS. DESIGN: We conducted a Phase Ib, placebo-controlled, randomised clinical trial of two dose levels of naproxen sodium (440 and 220 mg) administered daily for 6 months to 80 participants with LS, and a co-clinical trial using a genetically engineered mouse model of LS and patient-derived organoids (PDOs). RESULTS: Overall, the total number of adverse events was not different across treatment arms with excellent tolerance of the intervention. The level of prostaglandin E2 in the colorectal mucosa was significantly decreased after treatment with naproxen when compared with placebo. Naproxen activated different resident immune cell types without any increase in lymphoid cellularity, and changed the expression patterns of the intestinal crypt towards epithelial differentiation and stem cell regulation. Naproxen demonstrated robust chemopreventive activity in a mouse co-clinical trial and gene expression profiles induced by naproxen in humans showed perfect discrimination of mice specimens with LS and PDOs treated with naproxen and control. CONCLUSIONS: Naproxen is a promising strategy for immune interception in LS. We have discovered naproxen-induced gene expression profiles for their potential use as predictive biomarkers of drug activity. TRIAL REGISTRATION NUMBER: gov Identifier: NCT02052908.

CD8+ T cells inhibit metastasis and CXCL4 regulates its function.

BACKGROUND: The mechanism by which immune cells regulate metastasis is unclear. Understanding the role of immune cells in metastasis will guide the development of treatments improving patient survival. METHODS: We used syngeneic orthotopic mouse tumour models (wild-type, NOD/scid and Nude), employed knockout (CD8 and CD4) models and administered CXCL4. Tumours and lungs were analysed for cancer cells by bioluminescence, and circulating tumour cells were isolated from blood. Immunohistochemistry on the mouse tumours was performed to confirm cell type, and on a tissue microarray with 180 TNBCs for human relevance. TCGA data from over 10,000 patients were analysed as well. RESULTS: We reveal that intratumoral immune infiltration differs between metastatic and non-metastatic tumours. The non-metastatic tumours harbour high levels of CD8+ T cells and low levels of platelets, which is reverse in metastatic tumours. During tumour progression, platelets and CXCL4 induce differentiation of monocytes into myeloid-derived suppressor cells (MDSCs), which inhibit CD8+ T-cell function. TCGA pan-cancer data confirmed that CD8lowPlatelethigh patients have a significantly lower survival probability compared to CD8highPlateletlow. CONCLUSIONS: CD8+ T cells inhibit metastasis. When the balance between CD8+ T cells and platelets is disrupted, platelets produce CXCL4, which induces MDSCs thereby inhibiting the CD8+ T-cell function.

CD73 expression defines immune, molecular, and clinicopathological subgroups of lung adenocarcinoma.

INTRODUCTION: CD73 is a membrane-bound enzyme crucial in adenosine generation. The adenosinergic pathway plays a critical role in immunosuppression and in anti-tumor effects of immune checkpoint inhibitors (ICI). Here, we interrogated CD73 expression in a richly annotated cohort of human lung adenocarcinoma (LUAD) and its association with clinicopathological, immune, and molecular features to better understand the role of this immune marker in LUAD pathobiology. MATERIALS AND METHODS: Protein expression of CD73 was evaluated by immunohistochemistry in 106 archived LUADs from patients that underwent surgical treatment without neoadjuvant therapy. Total CD73 (T +) was calculated as the average of luminal (L +) and basolateral (BL +) percentage membrane expression scores for each LUAD and was used to classify tumors into three groups based on the extent of T CD73 expression (high, low, and negative). RESULTS: CD73 expression was significantly and progressively increased across normal-appearing lung tissue, adenomatous atypical hyperplasia, adenocarcinoma in situ, minimally invasive adenocarcinoma, and LUAD. In LUAD, BL CD73 expression was associated with an increase in PD-L1 expression in tumor cells and increase of tumor-associated immune cells. Stratification of LUADs based on T CD73 extent also revealed that tumors with high expression of this enzyme overall exhibited significantly elevated immune infiltration and PD-L1 protein expression. Immune profiling demonstrated that T-cell inflammation and adenosine signatures were significantly higher in CD73-expressing lung adenocarcinomas relative to those lacking CD73. CONCLUSION: Our study suggests that higher CD73 expression is associated with an overall augmented host immune response, suggesting potential implications in the immune pathobiology of early stage lung adenocarcinoma. Our findings warrant further studies to explore the role of CD73 in immunotherapeutic response of LUAD.

Long-circulating siRNA nanoparticles for validating Prohibitin1-targeted non-small cell lung cancer treatment.

RNA interference (RNAi) represents a promising strategy for identification and validation of putative therapeutic targets and for treatment of a myriad of important human diseases including cancer. However, the effective systemic in vivo delivery of small interfering RNA (siRNA) to tumors remains a formidable challenge. Using a robust self-assembly strategy, we develop a unique nanoparticle (NP) platform composed of a solid polymer/cationic lipid hybrid core and a lipid-poly(ethylene glycol) (lipid-PEG) shell for systemic siRNA delivery. The new generation lipid-polymer hybrid NPs are small and uniform, and can efficiently encapsulate siRNA and control its sustained release. They exhibit long blood circulation (t1/2 ∼ 8 h), high tumor accumulation, effective gene silencing, and negligible in vivo side effects. With this RNAi NP, we delineate and validate the therapeutic role of Prohibitin1 (PHB1), a target protein that has not been systemically evaluated in vivo due to the lack of specific and effective inhibitors, in treating non-small cell lung cancer (NSCLC) as evidenced by the drastic inhibition of tumor growth upon PHB1 silencing. Human tissue microarray analysis also reveals that high PHB1 tumor expression is associated with poorer overall survival in patients with NSCLC, further suggesting PHB1 as a therapeutic target. We expect this long-circulating RNAi NP platform to be of high interest for validating potential cancer targets in vivo and for the development of new cancer therapies.

Different Array CGH profiles within hereditary breast cancer tumors associated to BRCA1 expression and overall survival.

BACKGROUND: Array CGH analysis of breast tumors has contributed to the identification of different genomic profiles in these tumors. Loss of DNA repair by BRCA1 functional deficiency in breast cancer has been proposed as a relevant contribution to breast cancer progression for tumors with no germline mutation. Identifying the genomic alterations taking place in BRCA1 not expressing tumors will lead us to a better understanding of the cellular functions affected in this heterogeneous disease. Moreover, specific genomic alterations may contribute to the identification of potential therapeutic targets and offer a more personalized treatment to breast cancer patients. METHODS: Forty seven tumors from hereditary breast cancer cases, previously analyzed for BRCA1 expression, and screened for germline BRCA1 and 2 mutations, were analyzed by Array based Comparative Genomic Hybridization (aCGH) using Agilent 4x44K arrays. Overall survival was established for tumors in different clusters using Log-rank (Mantel-Cox) Test. Gene lists obtained from aCGH analysis were analyzed for Gene Ontology enrichment using GOrilla and DAVID tools. RESULTS: Genomic profiling of the tumors showed specific alterations associated to BRCA1 or 2 mutation status, and BRCA1 expression in the tumors, affecting relevant cellular processes. Similar cellular functions were found affected in BRCA1 not expressing and BRCA1 or 2 mutated tumors. Hierarchical clustering classified hereditary breast tumors in four major, groups according to the type and amount of genomic alterations, showing one group with a significantly poor overall survival (p = 0.0221). Within this cluster, deletion of PLEKHO1, GDF11, DARC, DAG1 and CD63 may be associated to the worse outcome of the patients. CONCLUSIONS: These results support the fact that BRCA1 lack of expression in tumors should be used as a marker for BRCAness and to select these patients for synthetic lethality approaches such as treatment with PARP inhibitors. In addition, the identification of specific alterations in breast tumors associated with poor survival, immune response or with a BRCAness phenotype will allow the use of a more personalized treatment in these patients.

Loss of the novel tumour suppressor and polarity gene Trim62 (Dear1) synergizes with oncogenic Ras in invasive lung cancer.

Deregulation of cell polarity proteins has been linked to the processes of invasion and metastasis. TRIM62 is a regulator of cell polarity and a tumour suppressor in breast cancer. Here, we demonstrate that human non-small cell lung cancer lesions show a step-wise loss of TRIM62 levels during disease progression, which was associated with poor clinical outcomes. To directly examine the role of Trim62 in development of lung cancer, we deleted Trim62 in a mutant K-Ras mouse model of lung cancer. In this context, haploinsufficiency of Trim62 synergized with a K-RasG12D mutation to promote invasiveness and disrupt three-dimensional morphogenesis, both of which are associated with epithelial-mesenchymal transitions. Re-expression of Trim62 reverted these phenotypes in tumour cell lines. Thus, Trim62 loss cooperates with K-Ras mutation in tumourigenesis and metastasis in vivo, indicating that decreased levels of TRIM62 may play an important role in the evolution of lung cancer.

Elucidating immune-related gene transcriptional programs via factorization of large-scale RNA-profiles.

Recent developments in immunotherapy, including immune checkpoint blockade (ICB) and adoptive cell therapy (ACT), have encountered challenges such as immune-related adverse events and resistance, especially in solid tumors. To advance the field, a deeper understanding of the molecular mechanisms behind treatment responses and resistance is essential. However, the lack of functionally characterized immune-related gene sets has limited data-driven immunological research. To address this gap, we adopted non-negative matrix factorization on 83 human bulk RNA sequencing (RNA-seq) datasets and constructed 28 immune-specific gene sets. After rigorous immunologist-led manual annotations and orthogonal validations across immunological contexts and functional omics data, we demonstrated that these gene sets can be applied to refine pan-cancer immune subtypes, improve ICB response prediction and functionally annotate spatial transcriptomic data. These functional gene sets, informing diverse immune states, will advance our understanding of immunology and cancer research.

Elucidating immune-related gene transcriptional programs via factorization of large-scale RNA-profiles.

Recent developments in immunotherapy, including immune checkpoint blockade (ICB) and adoptive cell therapy, have encountered challenges such as immune-related adverse events and resistance, especially in solid tumors. To advance the field, a deeper understanding of the molecular mechanisms behind treatment responses and resistance is essential. However, the lack of functionally characterized immune-related gene sets has limited data-driven immunological research. To address this gap, we adopted non-negative matrix factorization on 83 human bulk RNA-seq datasets and constructed 28 immune-specific gene sets. After rigorous immunologist-led manual annotations and orthogonal validations across immunological contexts and functional omics data, we demonstrated that these gene sets can be applied to refine pan-cancer immune subtypes, improve ICB response prediction and functionally annotate spatial transcriptomic data. These functional gene sets, informing diverse immune states, will advance our understanding of immunology and cancer research.