Human lung cancer harbors spatially organized stem-immunity hubs associated with response to immunotherapy.

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens and found an association with beneficial response to PD-1 blockade. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcome. This hub is distinct from mature tertiary lymphoid structures and is enriched for stem-like TCF7+PD-1+CD8+ T cells, activated CCR7+LAMP3+ dendritic cells and CCL19+ fibroblasts as well as chemokines that organize these cells. Within the stem-immunity hub, we find preferential interactions between CXCL10+ macrophages and TCF7-CD8+ T cells as well as between mature regulatory dendritic cells and TCF7+CD4+ and regulatory T cells. These results provide a picture of the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

Therapy-associated remodeling of pancreatic cancer revealed by single-cell spatial transcriptomics and optimal transport analysis.

In combination with cell intrinsic properties, interactions in the tumor microenvironment modulate therapeutic response. We leveraged high-plex single-cell spatial transcriptomics to dissect the remodeling of multicellular neighborhoods and cell-cell interactions in human pancreatic cancer associated with specific malignant subtypes and neoadjuvant chemotherapy/radiotherapy. We developed Spatially Constrained Optimal Transport Interaction Analysis (SCOTIA), an optimal transport model with a cost function that includes both spatial distance and ligand-receptor gene expression. Our results uncovered a marked change in ligand-receptor interactions between cancer-associated fibroblasts and malignant cells in response to treatment, which was supported by orthogonal datasets, including an ex vivo tumoroid co-culture system. Overall, this study demonstrates that characterization of the tumor microenvironment using high-plex single-cell spatial transcriptomics allows for identification of molecular interactions that may play a role in the emergence of chemoresistance and establishes a translational spatial biology paradigm that can be broadly applied to other malignancies, diseases, and treatments.

Spatial analysis of human lung cancer reveals organized immune hubs enriched for stem-like CD8 T cells and associated with immunotherapy response.

The organization of immune cells in human tumors is not well understood. Immunogenic tumors harbor spatially-localized multicellular 'immunity hubs' defined by expression of the T cell-attracting chemokines CXCL10/CXCL11 and abundant T cells. Here, we examined immunity hubs in human pre-immunotherapy lung cancer specimens, and found that they were associated with beneficial responses to PD-1-blockade. Immunity hubs were enriched for many interferon-stimulated genes, T cells in multiple differentiation states, and CXCL9/10/11 + macrophages that preferentially interact with CD8 T cells. Critically, we discovered the stem-immunity hub, a subtype of immunity hub strongly associated with favorable PD-1-blockade outcomes, distinct from mature tertiary lymphoid structures, and enriched for stem-like TCF7+PD-1+ CD8 T cells and activated CCR7 + LAMP3 + dendritic cells, as well as chemokines that organize these cells. These results elucidate the spatial organization of the human intratumoral immune response and its relevance to patient immunotherapy outcomes.

A single-nucleus and spatial transcriptomic atlas of the COVID-19 liver reveals topological, functional, and regenerative organ disruption in patients.

The molecular underpinnings of organ dysfunction in acute COVID-19 and its potential long-term sequelae are under intense investigation. To shed light on these in the context of liver function, we performed single-nucleus RNA-seq and spatial transcriptomic profiling of livers from 17 COVID-19 decedents. We identified hepatocytes positive for SARS-CoV-2 RNA with an expression phenotype resembling infected lung epithelial cells. Integrated analysis and comparisons with healthy controls revealed extensive changes in the cellular composition and expression states in COVID-19 liver, reflecting hepatocellular injury, ductular reaction, pathologic vascular expansion, and fibrogenesis. We also observed Kupffer cell proliferation and erythrocyte progenitors for the first time in a human liver single-cell atlas, resembling similar responses in liver injury in mice and in sepsis, respectively. Despite the absence of a clinical acute liver injury phenotype, endothelial cell composition was dramatically impacted in COVID-19, concomitantly with extensive alterations and profibrogenic activation of reactive cholangiocytes and mesenchymal cells. Our atlas provides novel insights into liver physiology and pathology in COVID-19 and forms a foundational resource for its investigation and understanding.

Single-nucleus and spatial transcriptome profiling of pancreatic cancer identifies multicellular dynamics associated with neoadjuvant treatment.

Pancreatic ductal adenocarcinoma (PDAC) is a highly lethal and treatment-refractory cancer. Molecular stratification in pancreatic cancer remains rudimentary and does not yet inform clinical management or therapeutic development. Here, we construct a high-resolution molecular landscape of the cellular subtypes and spatial communities that compose PDAC using single-nucleus RNA sequencing and whole-transcriptome digital spatial profiling (DSP) of 43 primary PDAC tumor specimens that either received neoadjuvant therapy or were treatment naive. We uncovered recurrent expression programs across malignant cells and fibroblasts, including a newly identified neural-like progenitor malignant cell program that was enriched after chemotherapy and radiotherapy and associated with poor prognosis in independent cohorts. Integrating spatial and cellular profiles revealed three multicellular communities with distinct contributions from malignant, fibroblast and immune subtypes: classical, squamoid-basaloid and treatment enriched. Our refined molecular and cellular taxonomy can provide a framework for stratification in clinical trials and serve as a roadmap for therapeutic targeting of specific cellular phenotypes and multicellular interactions.

System-wide transcriptome damage and tissue identity loss in COVID-19 patients.

The molecular mechanisms underlying the clinical manifestations of coronavirus disease 2019 (COVID-19), and what distinguishes them from common seasonal influenza virus and other lung injury states such as acute respiratory distress syndrome, remain poorly understood. To address these challenges, we combine transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues to define body-wide transcriptome changes in response to COVID-19. We then match these data with spatial protein and expression profiling across 357 tissue sections from 16 representative patient lung samples and identify tissue-compartment-specific damage wrought by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection, evident as a function of varying viral loads during the clinical course of infection and tissue-type-specific expression states. Overall, our findings reveal a systemic disruption of canonical cellular and transcriptional pathways across all tissues, which can inform subsequent studies to combat the mortality of COVID-19 and to better understand the molecular dynamics of lethal SARS-CoV-2 and other respiratory infections.

High-plex expression profiling reveals that implants drive spatiotemporal protein production and innate immune activation for tissue repair.

Surprisingly little clarity exists concerning effects of biomaterial properties on spatially localized protein expression, which drives implant success. Wound healing and tissue regeneration must be optimally supported by the implant, adsorbed proteins, immune cells, and fibroblasts; cells determine repair and functional recovery through protein production and regulation. However, not yet fully understood is how implants differentially drive spatial quantities of individual proteins both within the implant interior and the tissue surrounding it. Here we apply GeoMxⓇ digital spatial profiling to site-specifically investigate protein production in porous implants. Data is collected on the location and quantity of 40+ proteins from formalin-fixed, paraffin-embedded tissue slides of anisotropic tissue scaffolds (n = 18) with differing pore sizes (35 µm, 53 µm) and implantation durations (2, 14, 28 days); matching bulk gene expression data (700+ genes) is measured for identical implants. Notably, we discover fundamental spatial relationships in protein localization that in both the implant interior and the exterior are either uniquely independent or dependent of implant microstructure: dendritic cell marker CD11c and fibronectin significantly dominate the scaffold interior, while cell-to-cell adhesion marker CD34 and anti-inflammatory M2 polarization marker CD163 localize in the exterior. Lastly, collating spatial and bulk information, unique spatiotemporal expression patterns are identified for markers such as fibronectin, which are only uncoverable through spatial profiling and are otherwise hidden in bulk expression results. Together, these discoveries illustrate the critical importance of quantifying spatial expression patterns for implants, facilitating a paradigm shift in the iterative design, mechanistic understanding, and rapid assessment of biomaterials. STATEMENT OF SIGNIFICANCE: Spatial localization and expression of proteins, which determine implant success, are not fully understood because quantitative high-plex profiling is challenging. Applying GeoMxⓇ digital spatial profiling to site-specifically investigate protein production in porous implants, data is collected on the location and quantity of 40+ protein targets from tissue scaffolds with differing pore sizes (35 µm, 53 µm) and implantation durations (2, 14, 28 days). Collecting in parallel matched bulk gene expression data (700+ genes) for identical implants, we discover significant spatiotemporal expression patterns that remain otherwise hidden in differential bulk results. This new approach for the rapid assessment of biomaterials offers an enhanced mechanistic understanding and enables the tailoring of implants for superior regenerative outcomes.

Visualizing in deceased COVID-19 patients how SARS-CoV-2 attacks the respiratory and olfactory mucosae but spares the olfactory bulb.

Anosmia, the loss of smell, is a common and often the sole symptom of COVID-19. The onset of the sequence of pathobiological events leading to olfactory dysfunction remains obscure. Here, we have developed a postmortem bedside surgical procedure to harvest endoscopically samples of respiratory and olfactory mucosae and whole olfactory bulbs. Our cohort of 85 cases included COVID-19 patients who died a few days after infection with SARS-CoV-2, enabling us to catch the virus while it was still replicating. We found that sustentacular cells are the major target cell type in the olfactory mucosa. We failed to find evidence for infection of olfactory sensory neurons, and the parenchyma of the olfactory bulb is spared as well. Thus, SARS-CoV-2 does not appear to be a neurotropic virus. We postulate that transient insufficient support from sustentacular cells triggers transient olfactory dysfunction in COVID-19. Olfactory sensory neurons would become affected without getting infected.

Best Practices for Spatial Profiling for Breast Cancer Research with the GeoMx® Digital Spatial Profiler.

Breast cancer is a heterogenous disease with variability in tumor cells and in the surrounding tumor microenvironment (TME). Understanding the molecular diversity in breast cancer is critical for improving prediction of therapeutic response and prognostication. High-plex spatial profiling of tumors enables characterization of heterogeneity in the breast TME, which can holistically illuminate the biology of tumor growth, dissemination and, ultimately, response to therapy. The GeoMx Digital Spatial Profiler (DSP) enables researchers to spatially resolve and quantify proteins and RNA transcripts from tissue sections. The platform is compatible with both formalin-fixed paraffin-embedded and frozen tissues. RNA profiling was developed at the whole transcriptome level for human and mouse samples and protein profiling of 100-plex for human samples. Tissue can be optically segmented for analysis of regions of interest or cell populations to study biology-directed tissue characterization. The GeoMx Breast Cancer Consortium (GBCC) is composed of breast cancer researchers who are developing innovative approaches for spatial profiling to accelerate biomarker discovery. Here, the GBCC presents best practices for GeoMx profiling to promote the collection of high-quality data, optimization of data analysis and integration of datasets to advance collaboration and meta-analyses. Although the capabilities of the platform are presented in the context of breast cancer research, they can be generalized to a variety of other tumor types that are characterized by high heterogeneity.

The CD155/TIGIT axis promotes and maintains immune evasion in neoantigen-expressing pancreatic cancer.

The CD155/TIGIT axis can be co-opted during immune evasion in chronic viral infections and cancer. Pancreatic adenocarcinoma (PDAC) is a highly lethal malignancy, and immune-based strategies to combat this disease have been largely unsuccessful to date. We corroborate prior reports that a substantial portion of PDAC harbors predicted high-affinity MHC class I-restricted neoepitopes and extend these findings to advanced/metastatic disease. Using multiple preclinical models of neoantigen-expressing PDAC, we demonstrate that intratumoral neoantigen-specific CD8+ T cells adopt multiple states of dysfunction, resembling those in tumor-infiltrating lymphocytes of PDAC patients. Mechanistically, genetic and/or pharmacologic modulation of the CD155/TIGIT axis was sufficient to promote immune evasion in autochthonous neoantigen-expressing PDAC. Finally, we demonstrate that the CD155/TIGIT axis is critical in maintaining immune evasion in PDAC and uncover a combination immunotherapy (TIGIT/PD-1 co-blockade plus CD40 agonism) that elicits profound anti-tumor responses in preclinical models, now poised for clinical evaluation.

Spatially organized multicellular immune hubs in human colorectal cancer.

Immune responses to cancer are highly variable, with mismatch repair-deficient (MMRd) tumors exhibiting more anti-tumor immunity than mismatch repair-proficient (MMRp) tumors. To understand the rules governing these varied responses, we transcriptionally profiled 371,223 cells from colorectal tumors and adjacent normal tissues of 28 MMRp and 34 MMRd individuals. Analysis of 88 cell subsets and their 204 associated gene expression programs revealed extensive transcriptional and spatial remodeling across tumors. To discover hubs of interacting malignant and immune cells, we identified expression programs in different cell types that co-varied across tumors from affected individuals and used spatial profiling to localize coordinated programs. We discovered a myeloid cell-attracting hub at the tumor-luminal interface associated with tissue damage and an MMRd-enriched immune hub within the tumor, with activated T cells together with malignant and myeloid cells expressing T cell-attracting chemokines. By identifying interacting cellular programs, we reveal the logic underlying spatially organized immune-malignant cell networks.

Models that combine transcriptomic with spatial protein information exceed the predictive value for either single modality.

Immunotherapy has reshaped the field of cancer therapeutics but the population that benefits are small in many tumor types, warranting a companion diagnostic test. While immunohistochemistry (IHC) for programmed death-ligand 1 (PD-L1) or mismatch repair (MMR) and polymerase chain reaction (PCR) for microsatellite instability (MSI) are the only approved companion diagnostics others are under consideration. An optimal companion diagnostic test might combine the spatial information of IHC with the quantitative information from RNA expression profiling. Here, we show proof of concept for combination of spatially resolved protein information acquired by the NanoString GeoMx® Digital Spatial Profiler (DSP) with transcriptomic information from bulk mRNA gene expression acquired using NanoString nCounter® PanCancer IO 360™ panel on the same cohort of immunotherapy treated melanoma patients to create predictive models associated with clinical outcomes. We show that the combination of mRNA and spatially defined protein information can predict clinical outcomes more accurately (AUC 0.97) than either of these factors alone.

Systemic Tissue and Cellular Disruption from SARS-CoV-2 Infection revealed in COVID-19 Autopsies and Spatial Omics Tissue Maps.

The Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) virus has infected over 115 million people and caused over 2.5 million deaths worldwide. Yet, the molecular mechanisms underlying the clinical manifestations of COVID-19, as well as what distinguishes them from common seasonal influenza virus and other lung injury states such as Acute Respiratory Distress Syndrome (ARDS), remains poorly understood. To address these challenges, we combined transcriptional profiling of 646 clinical nasopharyngeal swabs and 39 patient autopsy tissues, matched with spatial protein and expression profiling (GeoMx) across 357 tissue sections. These results define both body-wide and tissue-specific (heart, liver, lung, kidney, and lymph nodes) damage wrought by the SARS-CoV-2 infection, evident as a function of varying viral load (high vs. low) during the course of infection and specific, transcriptional dysregulation in splicing isoforms, T cell receptor expression, and cellular expression states. In particular, cardiac and lung tissues revealed the largest degree of splicing isoform switching and cell expression state loss. Overall, these findings reveal a systemic disruption of cellular and transcriptional pathways from COVID-19 across all tissues, which can inform subsequent studies to combat the mortality of COVID-19, as well to better understand the molecular dynamics of lethal SARS-CoV-2 infection and other viruses.

The spatial landscape of lung pathology during COVID-19 progression.

Recent studies have provided insights into the pathology of and immune response to COVID-191-8. However, a thorough investigation of the interplay between infected cells and the immune system at sites of infection has been lacking. Here we use high-parameter imaging mass cytometry9 that targets the expression of 36 proteins to investigate the cellular composition and spatial architecture of acute lung injury in humans (including injuries derived from SARS-CoV-2 infection) at single-cell resolution. These spatially resolved single-cell data unravel the disordered structure of the infected and injured lung, alongside the distribution of extensive immune infiltration. Neutrophil and macrophage infiltration are hallmarks of bacterial pneumonia and COVID-19, respectively. We provide evidence that SARS-CoV-2 infects predominantly alveolar epithelial cells and induces a localized hyperinflammatory cell state that is associated with lung damage. We leverage the temporal range of fatal outcomes of COVID-19 in relation to the onset of symptoms, which reveals increased macrophage extravasation and increased numbers of mesenchymal cells and fibroblasts concomitant with increased proximity between these cell types as the disease progresses-possibly as a result of attempts to repair the damaged lung tissue. Our data enable us to develop a biologically interpretable landscape of lung pathology from a structural, immunological and clinical standpoint. We use this landscape to characterize the pathophysiology of the human lung from its macroscopic presentation to the single-cell level, which provides an important basis for understanding COVID-19 and lung pathology in general.

Shotgun transcriptome, spatial omics, and isothermal profiling of SARS-CoV-2 infection reveals unique host responses, viral diversification, and drug interactions.

In less than nine months, the Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) killed over a million people, including >25,000 in New York City (NYC) alone. The COVID-19 pandemic caused by SARS-CoV-2 highlights clinical needs to detect infection, track strain evolution, and identify biomarkers of disease course. To address these challenges, we designed a fast (30-minute) colorimetric test (LAMP) for SARS-CoV-2 infection from naso/oropharyngeal swabs and a large-scale shotgun metatranscriptomics platform (total-RNA-seq) for host, viral, and microbial profiling. We applied these methods to clinical specimens gathered from 669 patients in New York City during the first two months of the outbreak, yielding a broad molecular portrait of the emerging COVID-19 disease. We find significant enrichment of a NYC-distinctive clade of the virus (20C), as well as host responses in interferon, ACE, hematological, and olfaction pathways. In addition, we use 50,821 patient records to find that renin-angiotensin-aldosterone system inhibitors have a protective effect for severe COVID-19 outcomes, unlike similar drugs. Finally, spatial transcriptomic data from COVID-19 patient autopsy tissues reveal distinct ACE2 expression loci, with macrophage and neutrophil infiltration in the lungs. These findings can inform public health and may help develop and drive SARS-CoV-2 diagnostic, prevention, and treatment strategies.

Temporal and spatial heterogeneity of host response to SARS-CoV-2 pulmonary infection.

The relationship of SARS-CoV-2 pulmonary infection and severity of disease is not fully understood. Here we show analysis of autopsy specimens from 24 patients who succumbed to SARS-CoV-2 infection using a combination of different RNA and protein analytical platforms to characterize inter-patient and intra-patient heterogeneity of pulmonary virus infection. There is a spectrum of high and low virus cases associated with duration of disease. High viral cases have high activation of interferon pathway genes and a predominant M1-like macrophage infiltrate. Low viral cases are more heterogeneous likely reflecting inherent patient differences in the evolution of host response, but there is consistent indication of pulmonary epithelial cell recovery based on napsin A immunohistochemistry and RNA expression of surfactant and mucin genes. Using a digital spatial profiling platform, we find the virus corresponds to distinct spatial expression of interferon response genes demonstrating the intra-pulmonary heterogeneity of SARS-CoV-2 infection.

Comprehensive analysis of cutaneous and uveal melanoma liver metastases.

BACKGROUND: The profound disparity in response to immune checkpoint blockade (ICB) by cutaneous melanoma (CM) and uveal melanoma (UM) patients is not well understood. Therefore, we characterized metastases of CM and UM from the same metastatic site (liver), in order to dissect the potential underlying mechanism in differential response on ICB. METHODS: Tumor liver samples from CM (n=38) and UM (n=28) patients were analyzed at the genomic (whole exome sequencing), transcriptional (RNA sequencing) and protein (immunohistochemistry and GeoMx Digital Spatial Profiling) level. RESULTS: Comparison of CM and UM metastases from the same metastatic site revealed that, although originating from the same melanocyte lineage, CM and UM differed in somatic mutation profile, copy number profile, tumor mutational burden (TMB) and consequently predicted neoantigens. A higher melanin content and higher expression of the melanoma differentiation antigen MelanA was observed in liver metastases of UM patients. No difference in B2M and human leukocyte antigen-DR (HLA-DR) expression was observed. A higher expression of programmed cell death ligand 1 (PD-L1) was found in CM compared with UM liver metastases, although the majority of CM and UM liver metastases lacked PD-L1 expression. There was no difference in the extent of immune infiltration observed between CM and UM metastases, with the exception of a higher expression of CD163 (p<0.0001) in CM liver samples. While the extent of immune infiltration was similar for CM and UM metastases, the ratio of exhausted CD8 T cells to cytotoxic T cells, to total CD8 T cells and to Th1 cells, was significantly higher in UM metastases. CONCLUSIONS: While TMB was different between CM and UM metastases, tumor immune infiltration was similar. The greater dependency on PD-L1 as an immune checkpoint in CM and the identification of higher exhaustion ratios in UM may both serve as explanations for the difference in response to ICB. Consequently, in order to improve current treatment for metastatic UM, reversal of T cell exhaustion beyond programmed cell death 1 blockade should be considered.

Temporal and Spatial Heterogeneity of Host Response to SARS-CoV-2 Pulmonary Infection.

The relationship of SARS-CoV-2 lung infection and severity of pulmonary disease is not fully understood. We analyzed autopsy specimens from 24 patients who succumbed to SARS-CoV-2 infection using a combination of different RNA and protein analytical platforms to characterize inter- and intra- patient heterogeneity of pulmonary virus infection. There was a spectrum of high and low virus cases that was associated with duration of disease and activation of interferon pathway genes. Using a digital spatial profiling platform, the virus corresponded to distinct spatial expression of interferon response genes and immune checkpoint genes demonstrating the intra-pulmonary heterogeneity of SARS-CoV-2 infection.

Phase Ib/II trial testing combined radiofrequency ablation and ipilimumab in uveal melanoma (SECIRA-UM).

Approximately, 50% of patients with uveal melanoma develop distant metastasis for which no standard therapy is established. In contrast to cutaneous melanoma, the anti-CTLA-4 antibody ipilimumab showed no clinical activity in uveal melanoma. Liver directed therapies improve local control, but fail to show overall survival (OS) benefit. Preclinical experiments demonstrated that radiofrequency ablation (RFA) induced durable responses in combination with anti-CTLA-4. The aim of this phase Ib/II study was to assess safety and efficacy of RFA plus ipilimumab in uveal melanoma. Patients underwent RFA of one liver lesion and subsequently received four courses ipilimumab 0.3, 3 or 10 mg/kg every 3 weeks in a 3 + 3 design. Primary endpoints were safety in terms of dose limiting toxicities per cohort to define the recommended phase II dose (RP2D) in the phase Ib part and confirmed the objective response rate and disease control rate (DCR) of non-RFA lesions in the phase II part. Secondary endpoints were progression-free survival (PFS) and OS. Ipilimumab 10 mg/kg + RFA was initially defined as the RP2D. However, after 19 patients, the study was amended to adjust the RP2D to ipilimumab 3 mg/kg + RFA, because 47% of patients treated with 10 mg/kg had developed grade 3 colitis. In the 3 mg/kg cohort, also 19 patients have been treated. Immunotherapy-related grade ≥3 adverse events were observed in 53% of patients in the 10 mg/kg cohort versus 32% in the 3 mg/kg cohort. No confirmed objective responses were observed; the confirmed DCR was 5% in the 10 mg/kg cohort and 11% in the 3 mg/kg cohort. Median PFS was 3 months and comparable for both cohorts, median OS was 14.2 months for the 10 mg/kg cohort versus 9.7 months for the 3 mg/kg cohort. Combining RFA with ipilimumab 3 mg/kg was well tolerated, but showed very limited clinical activity in uveal melanoma.

Glutaminase is essential for the growth of triple-negative breast cancer cells with a deregulated glutamine metabolism pathway and its suppression synergizes with mTOR inhibition.

Tumor cells display fundamental changes in metabolism and nutrient uptake in order to utilize additional nutrient sources to meet their enhanced bioenergetic requirements. Glutamine (Gln) is one such nutrient that is rapidly taken up by tumor cells to fulfill this increased metabolic demand. A vital step in the catabolism of glutamine is its conversion to glutamate by the mitochondrial enzyme glutaminase (GLS). This study has identified GLS a potential therapeutic target in breast cancer, specifically in the basal subtype that exhibits a deregulated glutaminolysis pathway. Using inducible shRNA mediated gene knockdown, we discovered that loss of GLS function in triple-negative breast cancer (TNBC) cell lines with a deregulated glutaminolysis pathway led to profound tumor growth inhibition in vitro and in vivo. GLS knockdown had no effect on growth and metabolite levels in non-TNBC cell lines. We rescued the anti-tumor effect of GLS knockdown using shRNA resistant cDNAs encoding both GLS isoforms and by addition of an α-ketoglutarate (αKG) analog thus confirming the critical role of GLS in TNBC. Pharmacological inhibition of GLS with the small molecule inhibitor CB-839 reduced cell growth and led to a decrease in mammalian target of rapamycin (mTOR) activity and an increase in the stress response pathway driven by activating transcription factor 4 (ATF4). Finally, we found that GLS inhibition synergizes with mTOR inhibition, which introduces the possibility of a novel therapeutic strategy for TNBC. Our study revealed that GLS is essential for the survival of TNBC with a deregulated glutaminolysis pathway. The synergistic activity of GLS and mTOR inhibitors in TNBC cell lines suggests therapeutic potential of this combination for the treatment of vulnerable subpopulations of TNBC.