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.

The evolution of lung adenocarcinoma precursors is associated with chromosomal instability and transition from innate to adaptive immune response/evasion.

Studying lung adenocarcinoma (LUAD) early carcinogenesis is challenging, primarily due to the lack of LUAD precursors specimens. We amassed multi-omics data from 213 LUAD and LUAD precursors to identify molecular features underlying LUAD precancer evolution. We observed progressively increasing mutations, chromosomal aberrations, whole genome doubling and genomic instability from precancer to invasive LUAD, indicating aggravating chromosomal instability (CIN). Telomere shortening, a crucial genomic alteration linked to CIN, emerged at precancer stage. Moreover, later-stage lesions demonstrated increasing cancer stemness and decreasing alveolar identity, suggesting epithelial de-differentiation during early LUAD carcinogenesis. The innate immune cells progressively diminished from precancer to invasive LUAD, concomitant with a gradual recruitment of adaptive immune cells (except CD8+ and gamma-delta T cells that decreased in later stages) and upregulation of numerous immune checkpoints, suggesting LUAD precancer evolution is associated with a shift from innate to adaptive immune response and immune evasion mediated by various mechanisms.

eEF1A2 promotes PTEN-GSK3ß-SCF complex-dependent degradation of Aurora kinase A and is inactivated in breast cancer.

The translation elongation factor eEF1A promotes protein synthesis. Its methylation by METTL13 increases its activity, supporting tumor growth. However, in some cancers, a high abundance of eEF1A isoforms is associated with a good prognosis. Here, we found that eEF1A2 exhibited oncogenic or tumor-suppressor functions depending on its interaction with METTL13 or the phosphatase PTEN, respectively. METTL13 and PTEN competed for interaction with eEF1A2 in the same structural domain. PTEN-bound eEF1A2 promoted the ubiquitination and degradation of the mitosis-promoting Aurora kinase A in the S and G2 phases of the cell cycle. eEF1A2 bridged the interactions between the SKP1-CUL1-FBXW7 (SCF) ubiquitin ligase complex, the kinase GSK3ß, and Aurora-A, thereby facilitating the phosphorylation of Aurora-A in a degron site that was recognized by FBXW7. Genetic ablation of Eef1a2 or Pten in mice resulted in a greater abundance of Aurora-A and increased cell cycling in mammary tumors, which was corroborated in breast cancer tissues from patients. Reactivating this pathway using fimepinostat, which relieves inhibitory signaling directed at PTEN and increases FBXW7 expression, combined with inhibiting Aurora-A with alisertib, suppressed breast cancer cell proliferation in culture and tumor growth in vivo. The findings demonstrate a therapeutically exploitable, tumor-suppressive role for eEF1A2 in breast cancer.

Increased expression of CD70 in relapsed acute myeloid leukemia after hypomethylating agents.

Acute myeloid leukemia (AML) is the most common acute leukemia in adults. While induction chemotherapy leads to remission in most patients, a significant number will experience relapse. Therefore, there is a need for novel therapies that can improve remission rates in patients with relapsed and refractory AML. CD70 is the natural ligand for CD27 (a member of the TNF superfamily) and appears to be a promising therapeutic target. Consequently, there is considerable interest in developing chimeric antigen receptor (CAR) T-cell therapy products that can specifically target CD70 in various neoplasms, including AML. In this study, we employed routine diagnostic techniques, such as immunohistochemistry and flow cytometry, to investigate the expression of CD70 in bone marrow samples from treatment-naïve and relapsed AML patients after hypomethylating agents (HMA). Also, we evaluated the impact of HMA on CD70 expression and examined CD70 expression in various leukemic cell subsets and normal hematopoietic progenitors.

Brief Report: Comprehensive Clinicogenomic Profiling of Small Cell Transformation From EGFR-Mutant NSCLC Informs Potential Therapeutic Targets.

Introduction: NSCLC transformation to SCLC has been best characterized with EGFR-mutant NSCLC, with emerging case reports seen in ALK, RET, and KRAS-altered NSCLC. Previous reports revealed transformed SCLC from EGFR-mutant NSCLC portends very poor prognosis and lack effective treatment. Genomic analyses revealed TP53 and RB1 loss of function increase the risk of SCLC transformation. Little has been reported on the detailed clinicogenomic characteristics and potential therapeutic targets for this patient population. Methods: In this study, we conducted a single-center retrospective analysis of clinical and genomic characteristics of patients with EGFR-mutant NSCLC transformed to SCLC. Demographic data, treatment course, and clinical molecular testing reports were extracted from electronic medical records. Kaplan-Meier analyses were used to estimate survival outcomes. Next generation sequencing-based assays was used to identify EGFR and co-occurring genetic alterations in tissue or plasma before and after SCLC transformation. Single-cell RNA sequencing (scRNA-seq) was performed on a patient-derived-xenograft model generated from a patient with EGFR-NSCLC transformed SCLC tumor. Results: A total of 34 patients were identified in our study. Median age at initial diagnosis was 58, and median time to SCLC transformation was 24.2 months. 68% were female and 82% were never smokers. 79% of patients were diagnosed as stage IV disease, and over half had brain metastases at baseline. Median overall survival of the entire cohort was 38.3 months from initial diagnoses and 12.4 months from time of SCLC transformation. Most patients harbored EGFR exon19 deletions as opposed to exon21 L858R alteration. Continuing EGFR tyrosine kinase inhibitor post-transformation did not improve overall survival compared with those patients where tyrosine kinase inhibitor was stopped in our cohort. In the 20 paired pretransformed and post-transformed patient samples, statistically significant enrichment was seen with PIK3CA alterations (p = 0.04) post-transformation. Profiling of longitudinal liquid biopsy samples suggest emergence of SCLC genetic alterations before biopsy-proven SCLC, as shown by increasing variant allele frequency of TP53, RB1, PIK3CA alterations. ScRNA-seq revealed potential therapeutic targets including DLL3, CD276 (B7-H3) and PTK7 were widely expressed in transformed SCLC. Conclusions: SCLC transformation is a potential treatment resistance mechanism in driver-mutant NSCLC. In our cohort of 34 EGFR-mutant NSCLC, poor prognosis was observed after SCLC transformation. Clinicogenomic analyses of paired and longitudinal samples identified genomic alterations emerging post-transformation and scRNA-seq reveal potential therapeutic targets in this population. Further studies are needed to rigorously validate biomarkers and therapeutic targets for this patient population.

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.

Targeting eIF4A triggers an interferon response to synergize with chemotherapy and suppress triple-negative breast cancer.

Protein synthesis is frequently dysregulated in cancer and selective inhibition of mRNA translation represents an attractive cancer therapy. Here, we show that therapeutically targeting the RNA helicase eIF4A with zotatifin, the first-in-class eIF4A inhibitor, exerts pleiotropic effects on both tumor cells and the tumor immune microenvironment in a diverse cohort of syngeneic triple-negative breast cancer (TNBC) mouse models. Zotatifin not only suppresses tumor cell proliferation but also directly repolarizes macrophages toward an M1-like phenotype and inhibits neutrophil infiltration, which sensitizes tumors to immune checkpoint blockade. Mechanistic studies revealed that zotatifin reprograms the tumor translational landscape, inhibits the translation of Sox4 and Fgfr1, and induces an interferon (IFN) response uniformly across models. The induction of an IFN response is partially due to the inhibition of Sox4 translation by zotatifin. A similar induction of IFN-stimulated genes was observed in breast cancer patient biopsies following zotatifin treatment. Surprisingly, zotatifin significantly synergizes with carboplatin to trigger DNA damage and an even heightened IFN response, resulting in T cell-dependent tumor suppression. These studies identified a vulnerability of eIF4A in TNBC, potential pharmacodynamic biomarkers for zotatifin, and provide a rationale for new combination regimens consisting of zotatifin and chemotherapy or immunotherapy as treatments for TNBC.

Targeting EIF4A triggers an interferon response to synergize with chemotherapy and suppress triple-negative breast cancer.

Protein synthesis is frequently dysregulated in cancer and selective inhibition of mRNA translation represents an attractive cancer therapy. Here, we show that therapeutically targeting the RNA helicase eIF4A by Zotatifin, the first-in-class eIF4A inhibitor, exerts pleiotropic effects on both tumor cells and the tumor immune microenvironment in a diverse cohort of syngeneic triple-negative breast cancer (TNBC) mouse models. Zotatifin not only suppresses tumor cell proliferation but also directly repolarizes macrophages towards an M1-like phenotype and inhibits neutrophil infiltration, which sensitizes tumors to immune checkpoint blockade. Mechanistic studies revealed that Zotatifin reprograms the tumor translational landscape, inhibits the translation of Sox4 and Fgfr1, and induces an interferon response uniformly across models. The induction of an interferon response is partially due to the inhibition of Sox4 translation by Zotatifin. A similar induction of interferon-stimulated genes was observed in breast cancer patient biopsies following Zotatifin treatment. Surprisingly, Zotatifin significantly synergizes with carboplatin to trigger DNA damage and an even heightened interferon response resulting in T cell-dependent tumor suppression. These studies identified a vulnerability of eIF4A in TNBC, potential pharmacodynamic biomarkers for Zotatifin, and provide a rationale for new combination regimens comprising Zotatifin and chemotherapy or immunotherapy as treatments for TNBC.

Targeting immunosuppressive Ly6C+ classical monocytes reverses anti-PD-1/CTLA-4 immunotherapy resistance.

Introduction: Despite significant clinical advancement with the use of immune checkpoint blockade (ICB) in non-small cell lung cancer (NSCLC) there are still a major subset of patients that develop adaptive/acquired resistance. Understanding resistance mechanisms to ICB is critical to developing new therapeutic strategies and improving patient survival. The dynamic nature of the tumor microenvironment and the mutational load driving tumor immunogenicity limit the efficacy to ICB. Recent studies indicate that myeloid cells are drivers of ICB resistance. In this study we sought to understand which immune cells were contributing to resistance and if we could modify them in a way to improve response to ICB therapy. Results: Our results show that combination anti-PD-1/CTLA-4 produces an initial antitumor effect with evidence of an activated immune response. Upon extended treatment with anti-PD-1/CTLA-4 acquired resistance developed with an increase of the immunosuppressive populations, including T-regulatory cells, neutrophils and monocytes. Addition of anti-Ly6C blocking antibody to anti-PD-1/CTLA-4 was capable of completely reversing treatment resistance and restoring CD8 T cell activity in multiple KP lung cancer models and in the autochthonous lung cancer KrasLSL-G12D/p53fl/fl model. We found that there were higher classical Ly6C+ monocytes in anti-PD-1/CTLA-4 combination resistant tumors. B7 blockade illustrated the importance of dendritic cells for treatment efficacy of anti-Ly6C/PD-1/CTLA-4. We further determined that classical Ly6C+ monocytes in anti-PD-1/CTLA-4 resistant tumors are trafficked into the tumor via IFN-γ and the CCL2-CCR2 axis. Mechanistically we found that classical monocytes from ICB resistant tumors were unable to differentiate into antigen presenting cells and instead differentiated into immunosuppressive M2 macrophages or myeloid-derived suppressor cells (MDSC). Classical Ly6C+ monocytes from ICB resistant tumors had a decrease in both Flt3 and PU.1 expression that prevented differentiation into dendritic cells/macrophages. Conclusions: Therapeutically we found that addition of anti-Ly6C to the combination of anti-PD-1/CTLA-4 was capable of complete tumor eradication. Classical Ly6C+ monocytes differentiate into immunosuppressive cells, while blockade of classical monocytes drives dendritic cell differentiation/maturation to reinvigorate the anti-tumor T cell response. These findings support that immunotherapy resistance is associated with infiltrating monocytes and that controlling the differentiation process of monocytes can enhance the therapeutic potential of ICB.

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.

Naproxen chemoprevention induces proliferation of cytotoxic lymphocytes in Lynch Syndrome colorectal mucosa.

Background: Recent clinical trial data from Lynch Syndrome (LS) carriers demonstrated that naproxen administered for 6-months is a safe primary chemoprevention that promotes activation of different resident immune cell types without increasing lymphoid cellularity. While intriguing, the precise immune cell types enriched by naproxen remained unanswered. Here, we have utilized cutting-edge technology to elucidate the immune cell types activated by naproxen in mucosal tissue of LS patients. Methods: Normal colorectal mucosa samples (pre- and post-treatment) from a subset of patients enrolled in the randomized and placebo-controlled 'Naproxen Study' were obtained and subjected to a tissue microarray for image mass cytometry (IMC) analysis. IMC data was processed using tissue segmentation and functional markers to ascertain cell type abundance. Computational outputs were then used to quantitatively compare immune cell abundance in pre- and post-naproxen specimens. Results: Using data-driven exploration, unsupervised clustering identified four populations of immune cell types with statistically significant changes between treatment and control groups. These four populations collectively describe a unique cell population of proliferating lymphocytes within mucosal samples from LS patients exposed to naproxen. Conclusions: Our findings show that daily exposure of naproxen promotes T-cell proliferation in the colonic mucosa, which paves way for developing combination of immunoprevention strategies including naproxen for LS patients.

Immune cellular patterns of distribution affect outcomes of patients with non-small cell lung cancer.

Studying the cellular geographic distribution in non-small cell lung cancer is essential to understand the roles of cell populations in this type of tumor. In this study, we characterize the spatial cellular distribution of immune cell populations using 23 makers placed in five multiplex immunofluorescence panels and their associations with clinicopathologic variables and outcomes. Our results demonstrate two cellular distribution patterns-an unmixed pattern mostly related to immunoprotective cells and a mixed pattern mostly related to immunosuppressive cells. Distance analysis shows that T-cells expressing immune checkpoints are closer to malignant cells than other cells. Combining the cellular distribution patterns with cellular distances, we can identify four groups related to inflamed and not-inflamed tumors. Cellular distribution patterns and distance are associated with survival in univariate and multivariable analyses. Spatial distribution is a tool to better understand the tumor microenvironment, predict outcomes, and may can help select therapeutic interventions.

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.

PD-L1+ macrophages are associated with favorable features in primary mediastinal (thymic) large B-cell lymphoma.

Primary mediastinal (thymic) large B-cell lymphoma (PMBCL) is a rare, aggressive subtype of non-Hodgkin lymphoma and has a complex inflammatory microenvironment. Although most patients can be cured with standard-of-care immunochemotherapy, patients who have disease relapse have an unfavorable prognosis. Pre-treatment prognostic biomarkers in PMBCL are needed. In this retrospective study, we analyzed the clinical features and outcomes of PMBCL patients and their association with immune cell subpopulations identified by multiplex immunofluorescence at initial diagnosis. Two different antibody panels were used to assess macrophages in tissue biopsy specimens collected before the initiation of induction therapy. Twelve PMBCL patients, including five patients who had disease relapse, were included in the analysis. At a median follow-up time of 32.2 months, the median progression-free and overall survival durations were not reached. Our findings suggest that a high density of PD-L1+ macrophages is associated with favorable features, such as early disease stage and the absence of B-symptoms, and indicate that a high percentage of PD-L1+ macrophages and high densities of CD30+PD-L1+ cells and CD30+ cells might be associated with a lower risk of relapse within 12 months of therapy initiation. Further studies are needed to develop a biomarker signature predictive of treatment response with therapeutic consequences for patients with newly diagnosed PMBCL.

Neoadjuvant chemotherapy plus nivolumab with or without ipilimumab in operable non-small cell lung cancer: the phase 2 platform NEOSTAR trial.

Neoadjuvant ipilimumab + nivolumab (Ipi+Nivo) and nivolumab + chemotherapy (Nivo+CT) induce greater pathologic response rates than CT alone in patients with operable non-small cell lung cancer (NSCLC). The impact of adding ipilimumab to neoadjuvant Nivo+CT is unknown. Here we report the results and correlates of two arms of the phase 2 platform NEOSTAR trial testing neoadjuvant Nivo+CT and Ipi+Nivo+CT with major pathologic response (MPR) as the primary endpoint. MPR rates were 32.1% (7/22, 80% confidence interval (CI) 18.7-43.1%) in the Nivo+CT arm and 50% (11/22, 80% CI 34.6-61.1%) in the Ipi+Nivo+CT arm; the primary endpoint was met in both arms. In patients without known tumor EGFR/ALK alterations, MPR rates were 41.2% (7/17) and 62.5% (10/16) in the Nivo+CT and Ipi+Nivo+CT groups, respectively. No new safety signals were observed in either arm. Single-cell sequencing and multi-platform immune profiling (exploratory endpoints) underscored immune cell populations and phenotypes, including effector memory CD8+ T, B and myeloid cells and markers of tertiary lymphoid structures, that were preferentially increased in the Ipi+Nivo+CT cohort. Baseline fecal microbiota in patients with MPR were enriched with beneficial taxa, such as Akkermansia, and displayed reduced abundance of pro-inflammatory and pathogenic microbes. Neoadjuvant Ipi+Nivo+CT enhances pathologic responses and warrants further study in operable NSCLC. (ClinicalTrials.gov registration: NCT03158129 .).

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.

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.

Spatial molecular and cellular determinants of STAT3 activation in liver fibrosis progression in non-alcoholic fatty liver disease.

Background & Aims: The prevalence of non-alcoholic fatty liver disease (NAFLD) and its severe form, non-alcoholic steatohepatitis (NASH), is increasing. Individuals with NASH often develop liver fibrosis and advanced liver fibrosis is the main determinant of mortality in individuals with NASH. We and others have reported that STAT3 contributes to liver fibrosis and hepatocellular carcinoma in mice. Methods: Here, we explored whether STAT3 activation in hepatocyte and non-hepatocyte areas, measured by phospho-STAT3 (pSTAT3), is associated with liver fibrosis progression in 133 patients with NAFLD. We further characterized the molecular and cellular determinants of STAT3 activation by integrating spatial distribution and transcriptomic changes in fibrotic NAFLD livers.Results: pSTAT3 scores in non-hepatocyte areas progressively increased with fibrosis severity (r = 0.53, p <0.001). Correlation analyses between pSTAT3 scores and expression of 1,540 immune- and cancer-associated genes revealed a large effect of STAT3 activation on gene expression changes in non-hepatocyte areas and confirmed a major role for STAT3 activation in fibrogenesis. Digital spatial transcriptomic profiling was also performed on 13 regions selected in hepatocyte and non-hepatocyte areas from four NAFLD liver biopsies with advanced fibrosis, using a customized panel of markers including pSTAT3, PanCK+CK8/18, and CD45. The regions were further segmented based on positive or negative pSTAT3 staining. Cell deconvolution analysis revealed that activated STAT3 was enriched in hepatic progenitor cells (HPCs) and sinusoidal endothelial cells. Regression of liver fibrosis upon STAT3 inhibition in mice with NASH resulted in a reduction of HPCs, demonstrating a direct role for STAT3 in HPC expansion. Conclusion: Increased understanding of the spatial dependence of STAT3 signaling in NASH and liver fibrosis progression could lead to novel targeted treatment approaches. Impact and implications: Advanced liver fibrosis is the main determinant of mortality in patients with NASH. This study showed using liver biopsies from 133 patients with NAFLD, that STAT3 activation in non-hepatocyte areas is strongly associated with fibrosis severity, inflammation, and progression to NASH. STAT3 activation was enriched in hepatic progenitor cells (HPCs) and sinusoidal endothelial cells (SECs), as determined by innovative technologies interrogating the spatial distribution of pSTAT3. Finally, STAT3 inhibition in mice resulted in reduced liver fibrosis and depletion of HPCs, suggesting that STAT3 activation in HPCs contributes to their expansion and fibrogenesis in NAFLD.