PHGDH heterogeneity potentiates cancer cell dissemination and metastasis.

Cancer metastasis requires the transient activation of cellular programs enabling dissemination and seeding in distant organs1. Genetic, transcriptional and translational heterogeneity contributes to this dynamic process2,3. Metabolic heterogeneity has also been observed4, yet its role in cancer progression is less explored. Here we find that the loss of phosphoglycerate dehydrogenase (PHGDH) potentiates metastatic dissemination. Specifically, we find that heterogeneous or low PHGDH expression in primary tumours of patients with breast cancer is associated with decreased metastasis-free survival time. In mice, circulating tumour cells and early metastatic lesions are enriched with Phgdhlow cancer cells, and silencing Phgdh in primary tumours increases metastasis formation. Mechanistically, Phgdh interacts with the glycolytic enzyme phosphofructokinase, and the loss of this interaction activates the hexosamine-sialic acid pathway, which provides precursors for protein glycosylation. As a consequence, aberrant protein glycosylation occurs, including increased sialylation of integrin αvß3, which potentiates cell migration and invasion. Inhibition of sialylation counteracts the metastatic ability of Phgdhlow cancer cells. In conclusion, although the catalytic activity of PHGDH supports cancer cell proliferation, low PHGDH protein expression non-catalytically potentiates cancer dissemination and metastasis formation. Thus, the presence of PHDGH heterogeneity in primary tumours could be considered a sign of tumour aggressiveness.

Single-cell Atlas of Penile Cancer Reveals TP53 Mutations as a Driver of an Aggressive Phenotype, Irrespective of Human Papillomavirus Status, and Provides Clues for Treatment Personalization.

BACKGROUND AND OBJECTIVE: TP53 loss-of-function (TP53LOF) mutations might be a driver of poor prognosis and chemoresistance in both human papillomavirus (HPV)-independent (HPV-) and HPV-associated (HPV+) penile squamous cell carcinoma (PSCC). Here, we aim to describe transcriptomic differences in the PSCC microenvironment stratified by TP53LOF and HPV status. METHODS: We used single-cell RNA sequencing (scRNA-seq) and T-cell receptor sequencing to obtain a comprehensive atlas of the cellular architecture of PSCC. TP53LOF and HPV status were determined by targeted next-generation sequencing and sequencing HPV-DNA reads. Six HPV+ TP53 wild type (WT), six HPV- TP53WT, and four TP53LOF PSCC samples and six controls were included. Immunohistochemistry and hematoxylin-eosin confirmed the morphological context of the observed signatures. Prognostic differences between patient groups were validated in 541 PSCC patients using Kaplan-Meier survival estimates. KEY FINDINGS AND LIMITATIONS: Patients with aberrant p53 staining fare much worse than patients with either HPV- or HPV+ tumors and WT p53 expression. Using scRNA-seq, we revealed 65 cell subtypes within 83 682 cells. TP53LOF tumors exhibit a partial epithelial-to-mesenchymal transition, immune-excluded, angiogenic, and morphologically invasive environment, underlying their aggressive phenotype. HPV- TP53WT tumors show stemness and immune exhaustion. HPV+ TP53WT tumors mirror normal epithelial maturation with upregulation of antibody-drug-conjugate targets and activation of innate immunity. Inherent to the scRNA-seq analysis, low sample size is a limitation and validation of signatures in large PSCC cohorts is needed. CONCLUSIONS AND CLINICAL IMPLICATIONS: This first scRNA-seq atlas offers unprecedented in-depth insights into PSCC biology underlying prognostic differences based on TP53 and HPV status. Our findings provide clues for testing novel biomarker-driven therapies in PSCC. PATIENT SUMMARY: Here, we analyzed tissues of penile cancer at the level of individual cells, which helps us understand why patients who harbor a deactivating mutation in the TP53 gene do much worse than patients lacking such a mutation. Such an analysis may help us tailor future therapies based on TP53 gene mutations and human papillomavirus status of these tumors.

PD-1- CD45RA+ effector-memory CD8 T cells and CXCL10+ macrophages are associated with response to atezolizumab plus bevacizumab in advanced hepatocellular carcinoma.

The combination of atezolizumab plus bevacizumab (atezo/bev) has dramatically changed the treatment landscape of advanced HCC (aHCC), achieving durable responses in some patients. Using single-cell transcriptomics, we characterize the intra-tumoural and peripheral immune context of patients with aHCC treated with atezo/bev. Tumours from patients with durable responses are enriched for PDL1+ CXCL10+ macrophages and, based on cell-cell interaction analysis, express high levels of CXCL9/10/11 and are predicted to attract peripheral CXCR3+ CD8+ effector-memory T cells (CD8 TEM) into the tumour. Based on T cell receptor sharing and pseudotime trajectory analysis, we propose that CD8 TEM preferentially differentiate into clonally-expanded PD1- CD45RA+ effector-memory CD8+ T cells (CD8 TEMRA) with pronounced cytotoxicity. In contrast, in non-responders, CD8 TEM remain frozen in their effector-memory state. Finally, in responders, CD8 TEMRA display a high degree of T cell receptor sharing with blood, consistent with their patrolling activity. These findings may help understand the possible mechanisms underlying response to atezo/bev in aHCC.

Obesity-associated changes in molecular biology of primary breast cancer.

Obesity is associated with an increased risk of developing breast cancer (BC) and worse prognosis in BC patients, yet its impact on BC biology remains understudied in humans. This study investigates how the biology of untreated primary BC differs according to patients' body mass index (BMI) using data from >2,000 patients. We identify several genomic alterations that are differentially prevalent in overweight or obese patients compared to lean patients. We report evidence supporting an ageing accelerating effect of obesity at the genetic level. We show that BMI-associated differences in bulk transcriptomic profile are subtle, while single cell profiling allows detection of more pronounced changes in different cell compartments. These analyses further reveal an elevated and unresolved inflammation of the BC tumor microenvironment associated with obesity, with distinct characteristics contingent on the estrogen receptor status. Collectively, our analyses imply that obesity is associated with an inflammaging-like phenotype. We conclude that patient adiposity may play a significant role in the heterogeneity of BC and should be considered for BC treatment tailoring.

A single-cell map of intratumoral changes during anti-PD1 treatment of patients with breast cancer.

Immune-checkpoint blockade (ICB) combined with neoadjuvant chemotherapy improves pathological complete response in breast cancer. To understand why only a subset of tumors respond to ICB, patients with hormone receptor-positive or triple-negative breast cancer were treated with anti-PD1 before surgery. Paired pre- versus on-treatment biopsies from treatment-naive patients receiving anti-PD1 (n = 29) or patients receiving neoadjuvant chemotherapy before anti-PD1 (n = 11) were subjected to single-cell transcriptome, T cell receptor and proteome profiling. One-third of tumors contained PD1-expressing T cells, which clonally expanded upon anti-PD1 treatment, irrespective of tumor subtype. Expansion mainly involved CD8+ T cells with pronounced expression of cytotoxic-activity (PRF1, GZMB), immune-cell homing (CXCL13) and exhaustion markers (HAVCR2, LAG3), and CD4+ T cells characterized by expression of T-helper-1 (IFNG) and follicular-helper (BCL6, CXCR5) markers. In pre-treatment biopsies, the relative frequency of immunoregulatory dendritic cells (PD-L1+), specific macrophage phenotypes (CCR2+ or MMP9+) and cancer cells exhibiting major histocompatibility complex class I/II expression correlated positively with T cell expansion. Conversely, undifferentiated pre-effector/memory T cells (TCF7+, GZMK+) or inhibitory macrophages (CX3CR1+, C3+) were inversely correlated with T cell expansion. Collectively, our data identify various immunophenotypes and associated gene sets that are positively or negatively correlated with T cell expansion following anti-PD1 treatment. We shed light on the heterogeneity in treatment response to anti-PD1 in breast cancer.

Discriminating mild from critical COVID-19 by innate and adaptive immune single-cell profiling of bronchoalveolar lavages.

How the innate and adaptive host immune system miscommunicate to worsen COVID-19 immunopathology has not been fully elucidated. Here, we perform single-cell deep-immune profiling of bronchoalveolar lavage (BAL) samples from 5 patients with mild and 26 with critical COVID-19 in comparison to BALs from non-COVID-19 pneumonia and normal lung. We use pseudotime inference to build T-cell and monocyte-to-macrophage trajectories and model gene expression changes along them. In mild COVID-19, CD8+ resident-memory (TRM) and CD4+ T-helper-17 (TH17) cells undergo active (presumably antigen-driven) expansion towards the end of the trajectory, and are characterized by good effector functions, while in critical COVID-19 they remain more naïve. Vice versa, CD4+ T-cells with T-helper-1 characteristics (TH1-like) and CD8+ T-cells expressing exhaustion markers (TEX-like) are enriched halfway their trajectories in mild COVID-19, where they also exhibit good effector functions, while in critical COVID-19 they show evidence of inflammation-associated stress at the end of their trajectories. Monocyte-to-macrophage trajectories show that chronic hyperinflammatory monocytes are enriched in critical COVID-19, while alveolar macrophages, otherwise characterized by anti-inflammatory and antigen-presenting characteristics, are depleted. In critical COVID-19, monocytes contribute to an ATP-purinergic signaling-inflammasome footprint that could enable COVID-19 associated fibrosis and worsen disease-severity. Finally, viral RNA-tracking reveals infected lung epithelial cells, and a significant proportion of neutrophils and macrophages that are involved in viral clearance.

A pan-cancer blueprint of the heterogeneous tumor microenvironment revealed by single-cell profiling.

The stromal compartment of the tumor microenvironment consists of a heterogeneous set of tissue-resident and tumor-infiltrating cells, which are profoundly moulded by cancer cells. An outstanding question is to what extent this heterogeneity is similar between cancers affecting different organs. Here, we profile 233,591 single cells from patients with lung, colorectal, ovary and breast cancer (n = 36) and construct a pan-cancer blueprint of stromal cell heterogeneity using different single-cell RNA and protein-based technologies. We identify 68 stromal cell populations, of which 46 are shared between cancer types and 22 are unique. We also characterise each population phenotypically by highlighting its marker genes, transcription factors, metabolic activities and tissue-specific expression differences. Resident cell types are characterised by substantial tissue specificity, while tumor-infiltrating cell types are largely shared across cancer types. Finally, by applying the blueprint to melanoma tumors treated with checkpoint immunotherapy and identifying a naïve CD4+ T-cell phenotype predictive of response to checkpoint immunotherapy, we illustrate how it can serve as a guide to interpret scRNA-seq data. In conclusion, by providing a comprehensive blueprint through an interactive web server, we generate the first panoramic view on the shared complexity of stromal cells in different cancers.

Phenotype molding of stromal cells in the lung tumor microenvironment.

Cancer cells are embedded in the tumor microenvironment (TME), a complex ecosystem of stromal cells. Here, we present a 52,698-cell catalog of the TME transcriptome in human lung tumors at single-cell resolution, validated in independent samples where 40,250 additional cells were sequenced. By comparing with matching non-malignant lung samples, we reveal a highly complex TME that profoundly molds stromal cells. We identify 52 stromal cell subtypes, including novel subpopulations in cell types hitherto considered to be homogeneous, as well as transcription factors underlying their heterogeneity. For instance, we discover fibroblasts expressing different collagen sets, endothelial cells downregulating immune cell homing and genes coregulated with established immune checkpoint transcripts and correlating with T-cell activity. By assessing marker genes for these cell subtypes in bulk RNA-sequencing data from 1,572 patients, we illustrate how these correlate with survival, while immunohistochemistry for selected markers validates them as separate cellular entities in an independent series of lung tumors. Hence, in providing a comprehensive catalog of stromal cells types and by characterizing their phenotype and co-optive behavior, this resource provides deeper insights into lung cancer biology that will be helpful in advancing lung cancer diagnosis and therapy.