Beyond the Barrier: Unraveling the Mechanisms of Immunotherapy Resistance.

Immune checkpoint blockade (ICB) induces a remarkable and durable response in a subset of cancer patients. However, most patients exhibit either primary or acquired resistance to ICB. This resistance arises from a complex interplay of diverse dynamic mechanisms within the tumor microenvironment (TME). These mechanisms include genetic, epigenetic, and metabolic alterations that prevent T cell trafficking to the tumor site, induce immune cell dysfunction, interfere with antigen presentation, drive heightened expression of coinhibitory molecules, and promote tumor survival after immune attack. The TME worsens ICB resistance through the formation of immunosuppressive networks via immune inhibition, regulatory metabolites, and abnormal resource consumption. Finally, patient lifestyle factors, including obesity and microbiome composition, influence ICB resistance. Understanding the heterogeneity of cellular, molecular, and environmental factors contributing to ICB resistance is crucial to develop targeted therapeutic interventions that enhance the clinical response. This comprehensive overview highlights key mechanisms of ICB resistance that may be clinically translatable.

RNA Modification in the Immune System.

Characterization of RNA modifications has identified their distribution features and molecular functions. Dynamic changes in RNA modification on various forms of RNA are essential for the development and function of the immune system. In this review, we discuss the value of innovative RNA modification profiling technologies to uncover the function of these diverse, dynamic RNA modifications in various immune cells within healthy and diseased contexts. Further, we explore our current understanding of the mechanisms whereby aberrant RNA modifications modulate the immune milieu of the tumor microenvironment and point out outstanding research questions.

B Cell Function in the Tumor Microenvironment.

The tumor microenvironment (TME) is a heterogeneous, complex organization composed of tumor, stroma, and endothelial cells that is characterized by cross talk between tumor and innate and adaptive immune cells. Over the last decade, it has become increasingly clear that the immune cells in the TME play a critical role in controlling or promoting tumor growth. The function of T lymphocytes in this process has been well characterized. On the other hand, the function of B lymphocytes is less clear, although recent data from our group and others have strongly indicated a critical role for B cells in antitumor immunity. There are, however, a multitude of populations of B cells found within the TME, ranging from naive B cells all the way to terminally differentiated plasma cells and memory B cells. Here, we characterize the role of B cells in the TME in both animal models and patients, with an emphasis on dissecting how B cell heterogeneity contributes to the immune response to cancer.

Insights Gained from Single-Cell Analysis of Immune Cells in the Tumor Microenvironment.

Understanding tumor immune microenvironments is critical for identifying immune modifiers of cancer progression and developing cancer immunotherapies. Recent applications of single-cell RNA sequencing (scRNA-seq) in dissecting tumor microenvironments have brought important insights into the biology of tumor-infiltrating immune cells, including their heterogeneity, dynamics, and potential roles in both disease progression and response to immune checkpoint inhibitors and other immunotherapies. This review focuses on the advances in knowledge of tumor immune microenvironments acquired from scRNA-seq studies across multiple types of human tumors, with a particular emphasis on the study of phenotypic plasticity and lineage dynamics of immune cells in the tumor environment. We also discuss several imminent questions emerging from scRNA-seq observations and their potential solutions on the horizon.

Neoadjuvant PARPi or chemotherapy in ovarian cancer informs targeting effector Treg cells for homologous-recombination-deficient tumors.

Homologous recombination deficiency (HRD) is prevalent in cancer, sensitizing tumor cells to poly (ADP-ribose) polymerase (PARP) inhibition. However, the impact of HRD and related therapies on the tumor microenvironment (TME) remains elusive. Our study generates single-cell gene expression and T cell receptor profiles, along with validatory multimodal datasets from >100 high-grade serous ovarian cancer (HGSOC) samples, primarily from a phase II clinical trial (NCT04507841). Neoadjuvant monotherapy with the PARP inhibitor (PARPi) niraparib achieves impressive 62.5% and 73.6% response rates per RECIST v.1.1 and GCIG CA125, respectively. We identify effector regulatory T cells (eTregs) as key responders to HRD and neoadjuvant therapies, co-occurring with other tumor-reactive T cells, particularly terminally exhausted CD8+ T cells (Tex). TME-wide interferon signaling correlates with cancer cells upregulating MHC class II and co-inhibitory ligands, potentially driving Treg and Tex fates. Depleting eTregs in HRD mouse models, with or without PARP inhibition, significantly suppresses tumor growth without observable toxicities, underscoring the potential of eTreg-focused therapeutics for HGSOC and other HRD-related tumors.

A TCF4-dependent gene regulatory network confers resistance to immunotherapy in melanoma.

To better understand intrinsic resistance to immune checkpoint blockade (ICB), we established a comprehensive view of the cellular architecture of the treatment-naive melanoma ecosystem and studied its evolution under ICB. Using single-cell, spatial multi-omics, we showed that the tumor microenvironment promotes the emergence of a complex melanoma transcriptomic landscape. Melanoma cells harboring a mesenchymal-like (MES) state, a population known to confer resistance to targeted therapy, were significantly enriched in early on-treatment biopsies from non-responders to ICB. TCF4 serves as the hub of this landscape by being a master regulator of the MES signature and a suppressor of the melanocytic and antigen presentation transcriptional programs. Targeting TCF4 genetically or pharmacologically, using a bromodomain inhibitor, increased immunogenicity and sensitivity of MES cells to ICB and targeted therapy. We thereby uncovered a TCF4-dependent regulatory network that orchestrates multiple transcriptional programs and contributes to resistance to both targeted therapy and ICB in melanoma.

Macrophage-mediated myelin recycling fuels brain cancer malignancy.

Tumors growing in metabolically challenged environments, such as glioblastoma in the brain, are particularly reliant on crosstalk with their tumor microenvironment (TME) to satisfy their high energetic needs. To study the intricacies of this metabolic interplay, we interrogated the heterogeneity of the glioblastoma TME using single-cell and multi-omics analyses and identified metabolically rewired tumor-associated macrophage (TAM) subpopulations with pro-tumorigenic properties. These TAM subsets, termed lipid-laden macrophages (LLMs) to reflect their cholesterol accumulation, are epigenetically rewired, display immunosuppressive features, and are enriched in the aggressive mesenchymal glioblastoma subtype. Engulfment of cholesterol-rich myelin debris endows subsets of TAMs to acquire an LLM phenotype. Subsequently, LLMs directly transfer myelin-derived lipids to cancer cells in an LXR/Abca1-dependent manner, thereby fueling the heightened metabolic demands of mesenchymal glioblastoma. Our work provides an in-depth understanding of the immune-metabolic interplay during glioblastoma progression, thereby laying a framework to unveil targetable metabolic vulnerabilities in glioblastoma.

Pan-cancer proteogenomics expands the landscape of therapeutic targets.

Fewer than 200 proteins are targeted by cancer drugs approved by the Food and Drug Administration (FDA). We integrate Clinical Proteomic Tumor Analysis Consortium (CPTAC) proteogenomics data from 1,043 patients across 10 cancer types with additional public datasets to identify potential therapeutic targets. Pan-cancer analysis of 2,863 druggable proteins reveals a wide abundance range and identifies biological factors that affect mRNA-protein correlation. Integration of proteomic data from tumors and genetic screen data from cell lines identifies protein overexpression- or hyperactivation-driven druggable dependencies, enabling accurate predictions of effective drug targets. Proteogenomic identification of synthetic lethality provides a strategy to target tumor suppressor gene loss. Combining proteogenomic analysis and MHC binding prediction prioritizes mutant KRAS peptides as promising public neoantigens. Computational identification of shared tumor-associated antigens followed by experimental confirmation nominates peptides as immunotherapy targets. These analyses, summarized at https://targets.linkedomics.org, form a comprehensive landscape of protein and peptide targets for companion diagnostics, drug repurposing, and therapy development.

TULIPs decorate the three-dimensional genome of PFA ependymoma.

Posterior fossa group A (PFA) ependymoma is a lethal brain cancer diagnosed in infants and young children. The lack of driver events in the PFA linear genome led us to search its 3D genome for characteristic features. Here, we reconstructed 3D genomes from diverse childhood tumor types and uncovered a global topology in PFA that is highly reminiscent of stem and progenitor cells in a variety of human tissues. A remarkable feature exclusively present in PFA are type B ultra long-range interactions in PFAs (TULIPs), regions separated by great distances along the linear genome that interact with each other in the 3D nuclear space with surprising strength. TULIPs occur in all PFA samples and recur at predictable genomic coordinates, and their formation is induced by expression of EZHIP. The universality of TULIPs across PFA samples suggests a conservation of molecular principles that could be exploited therapeutically.

Early rhombic lip Protogenin+ve stem cells in a human-specific neurovascular niche initiate and maintain group 3 medulloblastoma.

We identify a population of Protogenin-positive (PRTG+ve) MYChigh NESTINlow stem cells in the four-week-old human embryonic hindbrain that subsequently localizes to the ventricular zone of the rhombic lip (RLVZ). Oncogenic transformation of early Prtg+ve rhombic lip stem cells initiates group 3 medulloblastoma (Gr3-MB)-like tumors. PRTG+ve stem cells grow adjacent to a human-specific interposed vascular plexus in the RLVZ, a phenotype that is recapitulated in Gr3-MB but not in other types of medulloblastoma. Co-culture of Gr3-MB with endothelial cells promotes tumor stem cell growth, with the endothelial cells adopting an immature phenotype. Targeting the PRTGhigh compartment of Gr3-MB in vivo using either the diphtheria toxin system or chimeric antigen receptor T cells constitutes effective therapy. Human Gr3-MBs likely arise from early embryonic RLVZ PRTG+ve stem cells inhabiting a specific perivascular niche. Targeting the PRTGhigh compartment and/or the perivascular niche represents an approach to treat children with Gr3-MB.

Pan-cancer single-cell dissection reveals phenotypically distinct B cell subtypes.

Characterizing the compositional and phenotypic characteristics of tumor-infiltrating B cells (TIBs) is important for advancing our understanding of their role in cancer development. Here, we establish a comprehensive resource of human B cells by integrating single-cell RNA sequencing data of B cells from 649 patients across 19 major cancer types. We demonstrate substantial heterogeneity in their total abundance and subtype composition and observe immunoglobulin G (IgG)-skewness of antibody-secreting cell isotypes. Moreover, we identify stress-response memory B cells and tumor-associated atypical B cells (TAABs), two tumor-enriched subpopulations with prognostic potential, shared in a pan-cancer manner. In particular, TAABs, characterized by a high clonal expansion level and proliferative capacity as well as by close interactions with activated CD4 T cells in tumors, are predictive of immunotherapy response. Our integrative resource depicts distinct clinically relevant TIB subsets, laying a foundation for further exploration of functional commonality and diversity of B cells in cancer.

A developmental constraint model of cancer cell states and tumor heterogeneity.

Cancer is a disease that stems from a fundamental liability inherent to multicellular life forms in which an individual cell is capable of reneging on the interests of the collective organism. Although cancer is commonly described as an evolutionary process, a less appreciated aspect of tumorigenesis may be the constraints imposed by the organism's developmental programs. Recent work from single-cell transcriptomic analyses across a range of cancer types has revealed the recurrence, plasticity, and co-option of distinct cellular states among cancer cell populations. Here, we note that across diverse cancer types, the observed cell states are proximate within the developmental hierarchy of the cell of origin. We thus posit a model by which cancer cell states are directly constrained by the organism's "developmental map." According to this model, a population of cancer cells traverses the developmental map, thereby generating a heterogeneous set of states whose interactions underpin emergent tumor behavior.

Pan-cancer proteogenomics characterization of tumor immunity.

Despite the successes of immunotherapy in cancer treatment over recent decades, less than <10%-20% cancer cases have demonstrated durable responses from immune checkpoint blockade. To enhance the efficacy of immunotherapies, combination therapies suppressing multiple immune evasion mechanisms are increasingly contemplated. To better understand immune cell surveillance and diverse immune evasion responses in tumor tissues, we comprehensively characterized the immune landscape of more than 1,000 tumors across ten different cancers using CPTAC pan-cancer proteogenomic data. We identified seven distinct immune subtypes based on integrative learning of cell type compositions and pathway activities. We then thoroughly categorized unique genomic, epigenetic, transcriptomic, and proteomic changes associated with each subtype. Further leveraging the deep phosphoproteomic data, we studied kinase activities in different immune subtypes, which revealed potential subtype-specific therapeutic targets. Insights from this work will facilitate the development of future immunotherapy strategies and enhance precision targeting with existing agents.

Neutrophil profiling illuminates anti-tumor antigen-presenting potency.

Neutrophils, the most abundant and efficient defenders against pathogens, exert opposing functions across cancer types. However, given their short half-life, it remains challenging to explore how neutrophils adopt specific fates in cancer. Here, we generated and integrated single-cell neutrophil transcriptomes from 17 cancer types (225 samples from 143 patients). Neutrophils exhibited extraordinary complexity, with 10 distinct states including inflammation, angiogenesis, and antigen presentation. Notably, the antigen-presenting program was associated with favorable survival in most cancers and could be evoked by leucine metabolism and subsequent histone H3K27ac modification. These neutrophils could further invoke both (neo)antigen-specific and antigen-independent T cell responses. Neutrophil delivery or a leucine diet fine-tuned the immune balance to enhance anti-PD-1 therapy in various murine cancer models. In summary, these data not only indicate the neutrophil divergence across cancers but also suggest therapeutic opportunities such as antigen-presenting neutrophil delivery.

Glioblastoma evolution and heterogeneity from a 3D whole-tumor perspective.

Treatment failure for the lethal brain tumor glioblastoma (GBM) is attributed to intratumoral heterogeneity and tumor evolution. We utilized 3D neuronavigation during surgical resection to acquire samples representing the whole tumor mapped by 3D spatial coordinates. Integrative tissue and single-cell analysis revealed sources of genomic, epigenomic, and microenvironmental intratumoral heterogeneity and their spatial patterning. By distinguishing tumor-wide molecular features from those with regional specificity, we inferred GBM evolutionary trajectories from neurodevelopmental lineage origins and initiating events such as chromothripsis to emergence of genetic subclones and spatially restricted activation of differential tumor and microenvironmental programs in the core, periphery, and contrast-enhancing regions. Our work depicts GBM evolution and heterogeneity from a 3D whole-tumor perspective, highlights potential therapeutic targets that might circumvent heterogeneity-related failures, and establishes an interactive platform enabling 360° visualization and analysis of 3D spatial patterns for user-selected genes, programs, and other features across whole GBM tumors.

Small-molecule GSDMD agonism in tumors stimulates antitumor immunity without toxicity.

Gasdermin-mediated inflammatory cell death (pyroptosis) can activate protective immunity in immunologically cold tumors. Here, we performed a high-throughput screen for compounds that could activate gasdermin D (GSDMD), which is expressed widely in tumors. We identified 6,7-dichloro-2-methylsulfonyl-3-N-tert-butylaminoquinoxaline (DMB) as a direct and selective GSDMD agonist that activates GSDMD pore formation and pyroptosis without cleaving GSDMD. In mouse tumor models, pulsed and low-level pyroptosis induced by DMB suppresses tumor growth without harming GSDMD-expressing immune cells. Protection is immune-mediated and abrogated in mice lacking lymphocytes. Vaccination with DMB-treated cancer cells protects mice from secondary tumor challenge, indicating that immunogenic cell death is induced. DMB treatment synergizes with anti-PD-1. DMB treatment does not alter circulating proinflammatory cytokine or leukocyte numbers or cause weight loss. Thus, our studies reveal a strategy that relies on a low level of tumor cell pyroptosis to induce antitumor immunity and raise the possibility of exploiting pyroptosis without causing overt toxicity.

Sex-dependent effects in the aged melanoma tumor microenvironment influence invasion and resistance to targeted therapy.

There is documented sex disparity in cutaneous melanoma incidence and mortality, increasing disproportionately with age and in the male sex. However, the underlying mechanisms remain unclear. While biological sex differences and inherent immune response variability have been assessed in tumor cells, the role of the tumor-surrounding microenvironment, contextually in aging, has been overlooked. Here, we show that skin fibroblasts undergo age-mediated, sex-dependent changes in their proliferation, senescence, ROS levels, and stress response. We find that aged male fibroblasts selectively drive an invasive, therapy-resistant phenotype in melanoma cells and promote metastasis in aged male mice by increasing AXL expression. Intrinsic aging in male fibroblasts mediated by EZH2 decline increases BMP2 secretion, which in turn drives the slower-cycling, highly invasive, and therapy-resistant melanoma cell phenotype, characteristic of the aged male TME. Inhibition of BMP2 activity blocks the emergence of invasive phenotypes and sensitizes melanoma cells to BRAF/MEK inhibition.

ARID1A suppresses R-loop-mediated STING-type I interferon pathway activation of anti-tumor immunity.

Clinical trials have identified ARID1A mutations as enriched among patients who respond favorably to immune checkpoint blockade (ICB) in several solid tumor types independent of microsatellite instability. We show that ARID1A loss in murine models is sufficient to induce anti-tumor immune phenotypes observed in ARID1A mutant human cancers, including increased CD8+ T cell infiltration and cytolytic activity. ARID1A-deficient cancers upregulated an interferon (IFN) gene expression signature, the ARID1A-IFN signature, associated with increased R-loops and cytosolic single-stranded DNA (ssDNA). Overexpression of the R-loop resolving enzyme, RNASEH2B, or cytosolic DNase, TREX1, in ARID1A-deficient cells prevented cytosolic ssDNA accumulation and ARID1A-IFN gene upregulation. Further, the ARID1A-IFN signature and anti-tumor immunity were driven by STING-dependent type I IFN signaling, which was required for improved responsiveness of ARID1A mutant tumors to ICB treatment. These findings define a molecular mechanism underlying anti-tumor immunity in ARID1A mutant cancers.

Circadian tumor infiltration and function of CD8+ T cells dictate immunotherapy efficacy.

The quality and quantity of tumor-infiltrating lymphocytes, particularly CD8+ T cells, are important parameters for the control of tumor growth and response to immunotherapy. Here, we show in murine and human cancers that these parameters exhibit circadian oscillations, driven by both the endogenous circadian clock of leukocytes and rhythmic leukocyte infiltration, which depends on the circadian clock of endothelial cells in the tumor microenvironment. To harness these rhythms therapeutically, we demonstrate that efficacy of chimeric antigen receptor T cell therapy and immune checkpoint blockade can be improved by adjusting the time of treatment during the day. Furthermore, time-of-day-dependent T cell signatures in murine tumor models predict overall survival in patients with melanoma and correlate with response to anti-PD-1 therapy. Our data demonstrate the functional significance of circadian dynamics in the tumor microenvironment and suggest the importance of leveraging these features for improving future clinical trial design and patient care.

Alanyl-tRNA synthetase, AARS1, is a lactate sensor and lactyltransferase that lactylates p53 and contributes to tumorigenesis.

Lysine lactylation is a post-translational modification that links cellular metabolism to protein function. Here, we find that AARS1 functions as a lactate sensor that mediates global lysine lacylation in tumor cells. AARS1 binds to lactate and catalyzes the formation of lactate-AMP, followed by transfer of lactate to the lysince acceptor residue. Proteomics studies reveal a large number of AARS1 targets, including p53 where lysine 120 and lysine 139 in the DNA binding domain are lactylated. Generation and utilization of p53 variants carrying constitutively lactylated lysine residues revealed that AARS1 lactylation of p53 hinders its liquid-liquid phase separation, DNA binding, and transcriptional activation. AARS1 expression and p53 lacylation correlate with poor prognosis among cancer patients carrying wild type p53. ß-alanine disrupts lactate binding to AARS1, reduces p53 lacylation, and mitigates tumorigenesis in animal models. We propose that AARS1 contributes to tumorigenesis by coupling tumor cell metabolism to proteome alteration.