Single-Cell Profiling of Sarcomas from Archival Tissue Reveals Programs Associated with Resistance to Immune Checkpoint Blockade.

PURPOSE: Sarcoma encompasses a diverse group of cancers that are typically resistant to current therapies, including immune checkpoint blockade (ICB), and underlying mechanisms are poorly understood. The contexture of sarcomas limits generation of high-quality data using cutting-edge molecular profiling methods, such as single-cell RNA-sequencing, thus hampering progress in understanding these understudied cancers. EXPERIMENTAL DESIGN: Here, we demonstrate feasibility of producing multimodal single-cell genomics and whole-genome sequencing data from frozen tissues, profiling 75,716 cell transcriptomes of five undifferentiated pleomorphic sarcoma and three intimal sarcoma samples, including paired specimens from two patients treated with ICB. RESULTS: We find that genomic diversity decreases in patients with response to ICB, and, in unbiased analyses, identify cancer cell programs associated with therapy resistance. Although interactions of tumor-infiltrating T lymphocytes within the tumor ecosystem increase in ICB responders, clonal expansion of CD8+ T cells alone was insufficient to predict drug responses. CONCLUSIONS: This study provides a framework for studying rare tumors and identifies salient and treatment-associated cancer cell intrinsic and tumor microenvironmental features in sarcomas.

Mapping variant effects on anti-tumor hallmarks of primary human T cells with base-editing screens.

Single-nucleotide variants (SNVs) in key T cell genes can drive clinical pathologies and could be repurposed to improve cellular cancer immunotherapies. Here, we perform massively parallel base-editing screens to generate thousands of variants at gene loci annotated with known or potential clinical relevance. We discover a broad landscape of putative gain-of-function (GOF) and loss-of-function (LOF) mutations, including in PIK3CD and the gene encoding its regulatory subunit, PIK3R1, LCK, SOS1, AKT1 and RHOA. Base editing of PIK3CD and PIK3R1 variants in T cells with an engineered T cell receptor specific to a melanoma epitope or in different generations of CD19 chimeric antigen receptor (CAR) T cells demonstrates that discovered GOF variants, but not LOF or silent mutation controls, enhanced signaling, cytokine production and lysis of cognate melanoma and leukemia cell models, respectively. Additionally, we show that generations of CD19 CAR T cells engineered with PIK3CD GOF mutations demonstrate enhanced antigen-specific signaling, cytokine production and leukemia cell killing, including when benchmarked against other recent strategies.

SMARCAL1 is a dual regulator of innate immune signaling and PD-L1 expression that promotes tumor immune evasion.

Genomic instability can trigger cancer-intrinsic innate immune responses that promote tumor rejection. However, cancer cells often evade these responses by overexpressing immune checkpoint regulators, such as PD-L1. Here, we identify the SNF2-family DNA translocase SMARCAL1 as a factor that favors tumor immune evasion by a dual mechanism involving both the suppression of innate immune signaling and the induction of PD-L1-mediated immune checkpoint responses. Mechanistically, SMARCAL1 limits endogenous DNA damage, thereby suppressing cGAS-STING-dependent signaling during cancer cell growth. Simultaneously, it cooperates with the AP-1 family member JUN to maintain chromatin accessibility at a PD-L1 transcriptional regulatory element, thereby promoting PD-L1 expression in cancer cells. SMARCAL1 loss hinders the ability of tumor cells to induce PD-L1 in response to genomic instability, enhances anti-tumor immune responses and sensitizes tumors to immune checkpoint blockade in a mouse melanoma model. Collectively, these studies uncover SMARCAL1 as a promising target for cancer immunotherapy.

Loss of Pip4k2c confers liver-metastatic organotropism through insulin-dependent PI3K-AKT pathway activation.

Liver metastasis (LM) confers poor survival and therapy resistance across cancer types, but the mechanisms of liver-metastatic organotropism remain unknown. Here, through in vivo CRISPR-Cas9 screens, we found that Pip4k2c loss conferred LM but had no impact on lung metastasis or primary tumor growth. Pip4k2c-deficient cells were hypersensitized to insulin-mediated PI3K/AKT signaling and exploited the insulin-rich liver milieu for organ-specific metastasis. We observed concordant changes in PIP4K2C expression and distinct metabolic changes in 3,511 patient melanomas, including primary tumors, LMs and lung metastases. We found that systemic PI3K inhibition exacerbated LM burden in mice injected with Pip4k2c-deficient cancer cells through host-mediated increase in hepatic insulin levels; however, this circuit could be broken by concurrent administration of an SGLT2 inhibitor or feeding of a ketogenic diet. Thus, this work demonstrates a rare example of metastatic organotropism through co-optation of physiological metabolic cues and proposes therapeutic avenues to counteract these mechanisms.

From patient tissue correlates to molecular mechanisms of cancer immune evasion: the emerging role of CD58 and PD-L1 co-regulation via CMTM6.

Immune evasion is a hallmark of cancer, yet the underlying mechanisms are often unknown in many patients. Using single-cell transcriptomics analysis, we previously identified the co-stimulator CD58 as part of a cancer cell-intrinsic immune checkpoint resistance signature in patient melanoma tissue. We subsequently validated CD58 loss as a driver of immune evasion using a patient-derived co-culture model of cancer and cytotoxic tumor-infiltrating lymphocytes in a pooled single-cell perturbation experiment, where we additionally observed concurrent upregulation of PD-L1 protein expression in melanoma cells with CD58 loss. In our most recent study, we uncovered the mechanisms of immune evasion mediated by CD58 loss, including impaired T cell activation and infiltration within tumors, as well as inhibitory signaling by PD-L1 via a shared regulator, CMTM6. Thus, cancer cell-intrinsic reduction of CD58 represents a multi-faceted determinant of immune evasion. Furthermore, its reciprocal interaction with PD-L1 via CMTM6 provides critical insights into how co-inhibitory and co-stimulatory immune cues are regulated.

Neoadjuvant durvalumab plus radiation versus durvalumab alone in stages I-III non-small cell lung cancer: survival outcomes and molecular correlates of a randomized phase II trial.

We previously reported the results of a randomized phase II trial (NCT02904954) in patients with early-stage non-small cell lung cancer (NSCLC) who were treated with either two preoperative cycles of the anti-PD-L1 antibody durvalumab alone or combined with immunomodulatory doses of stereotactic radiation (DRT). The trial met its primary endpoint of major pathological response, which was significantly higher following DRT with no new safety signals. Here, we report on the prespecified secondary endpoint of disease-free survival (DFS) regardless of treatment assignment and the prespecified exploratory analysis of DFS in each arm of the trial. DFS at 2 and 3 years across patients in both arms of the trial were 73% (95% CI: 62.1-84.5) and 65% (95% CI: 52.5-76.9) respectively. For the exploratory endpoint of DFS in each arm of the trial, three-year DFS was 63% (95% CI: 46.0-80.4) in the durvalumab monotherapy arm compared to 67% (95% CI: 49.6-83.4) in the dual therapy arm. In addition, we report post hoc exploratory analysis of progression-free survival as well as molecular correlates of response and recurrence through high-plex immunophenotyping of sequentially collected peripheral blood and gene expression profiles from resected tumors in both treatment arms. Together, our results contribute to the evolving landscape of neoadjuvant treatment regimens for NSCLC and identify easily measurable potential biomarkers of response and recurrence.

The CD58-CD2 axis is co-regulated with PD-L1 via CMTM6 and shapes anti-tumor immunity.

The cell-autonomous balance of immune-inhibitory and -stimulatory signals is a critical process in cancer immune evasion. Using patient-derived co-cultures, humanized mouse models, and single-cell RNA-sequencing of patient melanomas biopsied before and on immune checkpoint blockade, we find that intact cancer cell-intrinsic expression of CD58 and ligation to CD2 is required for anti-tumor immunity and is predictive of treatment response. Defects in this axis promote immune evasion through diminished T cell activation, impaired intratumoral T cell infiltration and proliferation, and concurrently increased PD-L1 protein stabilization. Through CRISPR-Cas9 and proteomics screens, we identify and validate CMTM6 as critical for CD58 stability and upregulation of PD-L1 upon CD58 loss. Competition between CD58 and PD-L1 for CMTM6 binding determines their rate of endosomal recycling over lysosomal degradation. Overall, we describe an underappreciated yet critical axis of cancer immunity and provide a molecular basis for how cancer cells balance immune inhibitory and stimulatory cues.

Multimodal single-cell and whole-genome sequencing of small, frozen clinical specimens.

Single-cell genomics enables dissection of tumor heterogeneity and molecular underpinnings of drug response at an unprecedented resolution1-11. However, broad clinical application of these methods remains challenging, due to several practical and preanalytical challenges that are incompatible with typical clinical care workflows, namely the need for relatively large, fresh tissue inputs. In the present study, we show that multimodal, single-nucleus (sn)RNA/T cell receptor (TCR) sequencing, spatial transcriptomics and whole-genome sequencing (WGS) are feasible from small, frozen tissues that approximate routinely collected clinical specimens (for example, core needle biopsies). Compared with data from sample-matched fresh tissue, we find a similar quality in the biological outputs of snRNA/TCR-seq data, while reducing artifactual signals and compositional biases introduced by fresh tissue processing. Profiling sequentially collected melanoma samples from a patient treated in the KEYNOTE-001 trial12, we resolved cellular, genomic, spatial and clonotype dynamics that represent molecular patterns of heterogeneous intralesional evolution during anti-programmed cell death protein 1 therapy. To demonstrate applicability to banked biospecimens of rare diseases13, we generated a single-cell atlas of uveal melanoma liver metastasis with matched WGS data. These results show that single-cell genomics from archival, clinical specimens is feasible and provides a framework for translating these methods more broadly to the clinical arena.

Massively parallel base editing screens to map variant effects on anti-tumor hallmarks of primary human T cells.

Base editing enables generation of single nucleotide variants, but large-scale screening in primary human T cells is limited due to low editing efficiency, among other challenges 1 . Here, we developed a high-throughput approach for high-efficiency and massively parallel adenine and cytosine base-editor screening in primary human T cells. We performed multiple large-scale screens editing 102 genes with central functions in T cells and full-length tiling mutagenesis of selected genes, and read out variant effects on hallmarks of T cell anti-tumor immunity, including activation, proliferation, and cytokine production. We discovered a broad landscape of gain- and loss-of-function mutations, including in PIK3CD and its regulatory subunit encoded by PIK3R1, LCK , AKT1, CTLA-4 and JAK1 . We identified variants that affected several (e.g., PIK3CD C416R) or only selected (e.g. LCK Y505C) hallmarks of T cell activity, and functionally validated several hits by probing downstream signaling nodes and testing their impact on T cell polyfunctionality and proliferation. Using primary human T cells in which we engineered a T cell receptor (TCR) specific to a commonly presented tumor testis antigen as a model for cellular immunotherapy, we demonstrate that base edits identified in our screens can tune specific or broad T cell functions and ultimately improve tumor elimination while exerting minimal off-target activity. In summary, we present the first large-scale base editing screen in primary human T cells and provide a framework for scalable and targeted base editing at high efficiency. Coupled with multi-modal phenotypic mapping, we accurately nominate variants that produce a desirable T cell state and leverage these synthetic proteins to improve models of cellular cancer immunotherapies.

Dissecting the treatment-naive ecosystem of human melanoma brain metastasis.

Melanoma brain metastasis (MBM) frequently occurs in patients with advanced melanoma; yet, our understanding of the underlying salient biology is rudimentary. Here, we performed single-cell/nucleus RNA-seq in 22 treatment-naive MBMs and 10 extracranial melanoma metastases (ECMs) and matched spatial single-cell transcriptomics and T cell receptor (TCR)-seq. Cancer cells from MBM were more chromosomally unstable, adopted a neuronal-like cell state, and enriched for spatially variably expressed metabolic pathways. Key observations were validated in independent patient cohorts, patient-derived MBM/ECM xenograft models, RNA/ATAC-seq, proteomics, and multiplexed imaging. Integrated spatial analyses revealed distinct geography of putative cancer immune evasion and evidence for more abundant intra-tumoral B to plasma cell differentiation in lymphoid aggregates in MBM. MBM harbored larger fractions of monocyte-derived macrophages and dysfunctional TOX+CD8+ T cells with distinct expression of immune checkpoints. This work provides comprehensive insights into MBM biology and serves as a foundational resource for further discovery and therapeutic exploration.

Structural variants shape driver combinations and outcomes in pediatric high-grade glioma.

We analyzed the contributions of structural variants (SVs) to gliomagenesis across 179 pediatric high-grade gliomas (pHGGs). The most recurrent SVs targeted MYC isoforms and receptor tyrosine kinases (RTKs), including an SV amplifying a MYC enhancer in 12% of diffuse midline gliomas (DMG), indicating an underappreciated role for MYC in pHGG. SV signature analysis revealed that tumors with simple signatures were TP53 wild type (TP53WT) but showed alterations in TP53 pathway members PPM1D and MDM4. Complex signatures were associated with direct aberrations in TP53, CDKN2A and RB1 early in tumor evolution and with later-occurring extrachromosomal amplicons. All pHGGs exhibited at least one simple-SV signature, but complex-SV signatures were primarily restricted to subsets of H3.3K27M DMGs and hemispheric pHGGs. Importantly, DMGs with complex-SV signatures were associated with shorter overall survival independent of histone mutation and TP53 status. These data provide insight into the impact of SVs on gliomagenesis and the mechanisms that shape them.

Dmp1Cre-directed knockdown of parathyroid hormone-related protein (PTHrP) in murine decidua is associated with a life-long increase in bone mass, width, and strength in male progeny.

Parathyroid hormone-related protein (PTHrP, gene name Pthlh) is a pleiotropic regulator of tissue homeostasis. In bone, Dmp1Cre-targeted PTHrP deletion in osteocytes causes osteopenia and impaired cortical strength. We report here that this outcome depends on parental genotype. In contrast to our previous report using mice bred from heterozygous (flox/wild type) Dmp1Cre.Pthlhf/w parents, adult (16-week-old and 26-week-old) flox/flox (f/f) Dmp1Cre.Pthlhf/f mice from homozygous parents (Dmp1Cre.Pthlhf/f(hom) ) have stronger bones, with 40% more trabecular bone mass and 30% greater femoral width than controls. This greater bone size was observed in Dmp1Cre.Pthlhf/f(hom) mice as early as 12 days of age, when greater bone width was also found in male and female Dmp1Cre.Pthlhf/f(hom) mice compared to controls, but not in gene-matched mice from heterozygous parents. This suggested a maternal influence on skeletal size prior to weaning. Although Dmp1Cre has previously been reported to cause gene recombination in mammary gland, milk PTHrP protein levels were normal. The wide-bone phenotype was also noted in utero: Dmp1Cre.Pthlhf/f(hom) embryonic femurs were more mineralized and wider than controls. Closer examination revealed that Dmp1Cre caused PTHrP recombination in placenta, and in the maternal-derived decidual layer that resides between the placenta and the uterus. Decidua from mothers of Dmp1Cre.Pthlhf/f(hom) mice also exhibited lower PTHrP levels by immunohistochemistry and were smaller than controls. We conclude that Dmp1Cre leads to gene recombination in decidua, and that decidual PTHrP might, through an influence on decidual cells, limit embryonic bone radial growth. This suggests a maternal-derived developmental origin of adult bone strength. © 2021 American Society for Bone and Mineral Research (ASBMR).

A molecular single-cell lung atlas of lethal COVID-19.

Respiratory failure is the leading cause of death in patients with severe SARS-CoV-2 infection1,2, but the host response at the lung tissue level is poorly understood. Here we performed single-nucleus RNA sequencing of about 116,000 nuclei from the lungs of nineteen individuals who died of COVID-19 and underwent rapid autopsy and seven control individuals. Integrated analyses identified substantial alterations in cellular composition, transcriptional cell states, and cell-to-cell interactions, thereby providing insight into the biology of lethal COVID-19. The lungs from individuals with COVID-19 were highly inflamed, with dense infiltration of aberrantly activated monocyte-derived macrophages and alveolar macrophages, but had impaired T cell responses. Monocyte/macrophage-derived interleukin-1ß and epithelial cell-derived interleukin-6 were unique features of SARS-CoV-2 infection compared to other viral and bacterial causes of pneumonia. Alveolar type 2 cells adopted an inflammation-associated transient progenitor cell state and failed to undergo full transition into alveolar type 1 cells, resulting in impaired lung regeneration. Furthermore, we identified expansion of recently described CTHRC1+ pathological fibroblasts3 contributing to rapidly ensuing pulmonary fibrosis in COVID-19. Inference of protein activity and ligand-receptor interactions identified putative drug targets to disrupt deleterious circuits. This atlas enables the dissection of lethal COVID-19, may inform our understanding of long-term complications of COVID-19 survivors, and provides an important resource for therapeutic development.

Multimodal pooled Perturb-CITE-seq screens in patient models define mechanisms of cancer immune evasion.

Resistance to immune checkpoint inhibitors (ICIs) is a key challenge in cancer therapy. To elucidate underlying mechanisms, we developed Perturb-CITE-sequencing (Perturb-CITE-seq), enabling pooled clustered regularly interspaced short palindromic repeat (CRISPR)-Cas9 perturbations with single-cell transcriptome and protein readouts. In patient-derived melanoma cells and autologous tumor-infiltrating lymphocyte (TIL) co-cultures, we profiled transcriptomes and 20 proteins in ~218,000 cells under ~750 perturbations associated with cancer cell-intrinsic ICI resistance (ICR). We recover known mechanisms of resistance, including defects in the interferon-γ (IFN-γ)-JAK/STAT and antigen-presentation pathways in RNA, protein and perturbation space, and new ones, including loss/downregulation of CD58. Loss of CD58 conferred immune evasion in multiple co-culture models and was downregulated in tumors of melanoma patients with ICR. CD58 protein expression was not induced by IFN-γ signaling, and CD58 loss conferred immune evasion without compromising major histocompatibility complex (MHC) expression, suggesting that it acts orthogonally to known mechanisms of ICR. This work provides a framework for the deciphering of complex mechanisms by large-scale perturbation screens with multimodal, single-cell readouts, and discovers potentially clinically relevant mechanisms of immune evasion.

CloneSifter: enrichment of rare clones from heterogeneous cell populations.

BACKGROUND: Many biological processes, such as cancer metastasis, organismal development, and acquisition of resistance to cytotoxic therapy, rely on the emergence of rare sub-clones from a larger population. Understanding how the genetic and epigenetic features of diverse clones affect clonal fitness provides insight into molecular mechanisms underlying selective processes. While large-scale barcoding with NGS readout has facilitated cellular fitness assessment at the population level, this approach does not support characterization of clones prior to selection. Single-cell genomics methods provide high biological resolution, but are challenging to scale across large populations to probe rare clones and are destructive, limiting further functional analysis of important clones. RESULTS: Here, we develop CloneSifter, a methodology for tracking and enriching rare clones throughout their response to selection. CloneSifter utilizes a CRISPR sgRNA-barcode library that facilitates the isolation of viable cells from specific clones within the barcoded population using a sequence-specific retrieval reporter. We demonstrate that CloneSifter can measure clonal fitness of cancer cell models in vitro and retrieve targeted clones at abundance as low as 1 in 1883 in a heterogeneous cell population. CONCLUSIONS: CloneSifter provides a means to track and access specific and rare clones of interest across dynamic changes in population structure to comprehensively explore the basis of these changes.

Mechanisms and therapeutic implications of hypermutation in gliomas.

A high tumour mutational burden (hypermutation) is observed in some gliomas1-5; however, the mechanisms by which hypermutation develops and whether it predicts the response to immunotherapy are poorly understood. Here we comprehensively analyse the molecular determinants of mutational burden and signatures in 10,294 gliomas. We delineate two main pathways to hypermutation: a de novo pathway associated with constitutional defects in DNA polymerase and mismatch repair (MMR) genes, and a more common post-treatment pathway, associated with acquired resistance driven by MMR defects in chemotherapy-sensitive gliomas that recur after treatment with the chemotherapy drug temozolomide. Experimentally, the mutational signature of post-treatment hypermutated gliomas was recapitulated by temozolomide-induced damage in cells with MMR deficiency. MMR-deficient gliomas were characterized by a lack of prominent T cell infiltrates, extensive intratumoral heterogeneity, poor patient survival and a low rate of response to PD-1 blockade. Moreover, although bulk analyses did not detect microsatellite instability in MMR-deficient gliomas, single-cell whole-genome sequencing analysis of post-treatment hypermutated glioma cells identified microsatellite mutations. These results show that chemotherapy can drive the acquisition of hypermutated populations without promoting a response to PD-1 blockade and supports the diagnostic use of mutational burden and signatures in cancer.

Brief exposure to full length parathyroid hormone-related protein (PTHrP) causes persistent generation of cyclic AMP through an endocytosis-dependent mechanism.

Parathyroid hormone (PTH)-related protein (PTHrP) (gene name Pthlh) was discovered as the factor responsible for the humoral hypercalcemia of malignancy. It shares such sequence similarity with PTH in the amino-terminal region that the two are equally able to act through a single G protein-coupled receptor, PTH1R. A number of biological activities are ascribed to domains of PTHrP beyond the amino-terminal domain. PTH functions as a circulating hormone, but PTHrP is generated locally in many tissues including bone, where it acts as a paracrine factor on osteoblasts and osteocytes. The present study compares how PTH and PTHrP influence cyclic AMP (cAMP) formation through adenylyl cyclase, the first event in cell activation through PTH1R. Brief exposure to full length PTHrP(1-141) in several osteoblastic cell culture systems was followed by sustained adenylyl cyclase activity for more than an hour after ligand washout. This effect was dose-dependent and was not found with shorter PTHrP or PTH peptides even though they were fully able to activate adenylyl cyclase with acute treatment. The persistent activation response to PTHrP(1-141) was seen also with later events in the cAMP/PKA pathway, including persistent activation of CRE-luciferase and sustained regulation of several CREB-responsive mRNAs, up to 24 h after the initial exposure. Pharmacologic blockade of endocytosis prevented the persistent activation of cAMP and gene responses. We conclude that full length PTHrP, the likely local physiological effector in bone, differs in intracellular action to PTH by undergoing endosomal translocation to induce a prolonged adenylyl cyclase activation in its target cells.

Neuronal differentiation and cell-cycle programs mediate response to BET-bromodomain inhibition in MYC-driven medulloblastoma.

BET-bromodomain inhibition (BETi) has shown pre-clinical promise for MYC-amplified medulloblastoma. However, the mechanisms for its action, and ultimately for resistance, have not been fully defined. Here, using a combination of expression profiling, genome-scale CRISPR/Cas9-mediated loss of function and ORF/cDNA driven rescue screens, and cell-based models of spontaneous resistance, we identify bHLH/homeobox transcription factors and cell-cycle regulators as key genes mediating BETi's response and resistance. Cells that acquire drug tolerance exhibit a more neuronally differentiated cell-state and expression of lineage-specific bHLH/homeobox transcription factors. However, they do not terminally differentiate, maintain expression of CCND2, and continue to cycle through S-phase. Moreover, CDK4/CDK6 inhibition delays acquisition of resistance. Therefore, our data provide insights about the mechanisms underlying BETi effects and the appearance of resistance and support the therapeutic use of combined cell-cycle inhibitors with BETi in MYC-amplified medulloblastoma.