Cell
第 187 卷,第 6 期,2024 年 3 月 14 日,第 1422-1439.e24 页
Article 文章Neutrophil profiling illuminates anti-tumor antigen-presenting potency
中性粒细胞分析阐明了抗肿瘤抗原呈递效力
https://doi.org/10.1016/j.cell.2024.02.005IF:64.5 第一季度Get rights and content 获取权利和内容
Graphical abstract 图形概要
Keywords 关键词
肿瘤相关中性粒细胞单细胞RNA测序泛癌分析抗原呈递新抗原免疫疗法
Introduction 介绍
Neutrophils are believed to be the cells that achieve the most rapid defense against pathogens.1 They can sense diverse cancer signals, such as inflammation and wounding, to initiate the chemotaxis module toward the tumor microenvironment (TME).2 However, human neutrophils are usually too short-lived to be profiled (half-life: 6–8 h),3 rendering most single-cell RNA sequencing (scRNA-seq) approaches unable to achieve high-throughput profiling of these cells. Due to their lower mRNA content (neutrophils: 0.33 μg; macrophage and monocytes: 2.55 μg per million cells),4 our understanding of the transcriptional diversity and spatiotemporal heterogeneity of human neutrophils remains rudimentary, despite the universal distribution of these cells across organs and tumors throughout the body.
中性粒细胞被认为是对病原体实现最快速防御的细胞。 1 它们可以感知多种癌症信号,例如炎症和受伤,以启动针对肿瘤微环境(TME)的趋化模块。 2 然而,人类中性粒细胞的寿命通常太短,无法进行分析(半衰期:6-8 小时), 3 导致大多数单细胞 RNA 测序 (scRNA-seq)方法无法实现这些细胞的高通量分析。由于其 mRNA 含量较低(中性粒细胞:0.33 μg;巨噬细胞和单核细胞:每百万细胞 2.55 μg), 4 我们对人类中性粒细胞转录多样性和时空异质性的了解仍然很初级,尽管这些细胞遍布全身的器官和肿瘤。
In the field of cancer immunology, the paradox of whether neutrophils are generally suppressive or protective remains unresolved. Tumor-associated neutrophils are long believed to be immunosuppressive5,6 and to exacerbate patient outcomes.7,8,9 However, these cells were recently shown to kill cancer cells by releasing active elastase,10 nitric oxide synthase,11 or reactive oxygen species (ROS).12 Alternatively, neutrophils also harbored anti-tumor immune phenotype that promotes autologous T cell responses13 or interferon-related immunostimulatory effects.14 These seemingly contradictory data raise critical but poorly understood questions around population composition and which subsets drive pro- or anti-tumor effects. Systematically decoding the cellular diversity of tumor-infiltrating neutrophils will identify their diverse gene expression patterns as well as their niche framework. In this context, single-cell profiling is uniquely suited for the characterization of neutrophil states and provides the opportunity to create a data-driven map of neutrophil ontology.
在癌症免疫学领域,中性粒细胞通常是抑制性的还是保护性的悖论仍未解决。长期以来,肿瘤相关中性粒细胞被认为具有免疫抑制作用 5 6 并会加剧患者的预后。 7 8 9 然而,最近显示这些细胞通过释放活性弹性蛋白酶、 10 一氧化氮合酶、 或活性氧 (ROS)。 12 另外,中性粒细胞还具有抗肿瘤免疫表型,可促进自体 T 细胞反应 13 或干扰素相关的免疫刺激作用。 14 这些看似矛盾的数据提出了关于人口构成以及哪些子集驱动促肿瘤或抗肿瘤作用的关键但人们知之甚少的问题。系统地解码肿瘤浸润中性粒细胞的细胞多样性将识别其不同的基因表达模式及其生态位框架。在这种情况下,单细胞分析特别适合中性粒细胞状态的表征,并提供了创建数据驱动的中性粒细胞本体图的机会。
To address this challenge, we here designed a neutrophil profiling strategy and generated the single neutrophil transcriptomes from 225 samples collected from 143 patients across 17 cancer types, including paired metastases from selected cancers. We found that neutrophils exhibit a complex and diverse transcriptional profile with 10 distinct cell states, among which three were potentially dominant across various cancer types, including inflammation, angiogenesis, and antigen presentation. In particular, antigen-presenting neutrophils showed unique immunophenotypes and metabolic features, which can induce T cell neoantigen reactivity. Our study not only generates cancer neutrophil transcriptomes but also unravels potential therapeutic opportunities such as antigen-presenting neutrophil delivery.
为了应对这一挑战,我们设计了中性粒细胞分析策略,并从 17 种癌症类型(包括来自选定癌症的配对转移)的 143 名患者收集的 225 个样本中生成了单个中性粒细胞转录组。我们发现中性粒细胞表现出复杂多样的转录谱,具有 10 种不同的细胞状态,其中三种在各种癌症类型中可能占主导地位,包括炎症、血管生成和抗原呈递。特别是,抗原呈递中性粒细胞表现出独特的免疫表型和代谢特征,可以诱导 T 细胞新抗原反应。我们的研究不仅产生癌症中性粒细胞转录组,而且揭示了潜在的治疗机会,例如抗原呈递中性粒细胞递送。
Results 结果
Neutrophils preferably infiltrate into certain cancer types
中性粒细胞优先浸润某些癌症类型
The extent of neutrophil infiltration into solid tumors varies widely,15 but there still exists no consensus on the exact cancer types and infiltration level. To explore the infiltration patterns and choose the appropriate cancer types, we tested 8 common immune quantification algorithms and analyzed The Cancer Genome Atlas (TCGA) covering 8,766 samples across 31 solid cancers (STAR Methods). We developed a consensus neutrophil infiltration score based on three algorithms,16,17,18 which revealed strong tissue-selective patterns of neutrophil infiltration that could be clustered into 3 subtypes (29.6% high, 30.0% heterogeneous, and 40.4% low; Figures 1A and S1A). For example, neutrophils show high infiltration in lung and kidney cancers and intermediate infiltration in gastrointestinal cancers, which was consistent in another pan-cancer dataset, Clinical Proteomic Tumor Analysis Consortium (CPTAC)19 (Figures S1B and S1C). By comparing neutrophil levels among immune subtypes,20,21 we observed preferential neutrophil infiltration in inflammatory or fibrotic suppressive TMEs (Figure 1B), in agreement with previously reported neutrophil infiltration patterns and functions.15 Together, these data highlighted the diversity of neutrophil infiltration depending upon the tissue and cancer types, providing a basis for our subsequent neutrophil sampling strategy.
中性粒细胞浸润实体瘤的程度差异很大, 15 但对于确切的癌症类型和浸润水平仍没有达成共识。为了探索浸润模式并选择合适的癌症类型,我们测试了 8 种常见的免疫定量算法,并分析了涵盖 31 种实体癌的 8,766 个样本的癌症基因组图谱 (TCGA)(STAR 方法)。我们基于三种算法制定了一致的中性粒细胞浸润评分, 16 17 18 ,它揭示了中性粒细胞浸润的强组织选择性模式,可分为 3 种亚型(29.6% 高,30.0% 异质性,40.4% 低;图 1A 和 S1A)。例如,中性粒细胞在肺癌和肾癌中表现出高浸润性,在胃肠道癌中表现出中等浸润性,这与另一个泛癌症数据集临床蛋白质组肿瘤分析联盟 (CPTAC) 19 一致(图 S1B 和 S1C) 。通过比较免疫亚型之间的中性粒细胞水平, 20 21 ,我们观察到炎症或纤维化抑制性 TME 中中性粒细胞优先浸润(图 1B),这与之前报道的中性粒细胞浸润模式和功能一致。 15 总之,这些数据突出了中性粒细胞浸润的多样性,具体取决于组织和癌症类型,为我们后续的中性粒细胞采样策略提供了基础。
Following our initial findings, we devised a standardized sampling strategy focused on cancer types with high or medium neutrophil infiltration (sample statistics, Figure S1D), ensuring the inclusion of matched blood, adjacent normal tissues, and metastasis samples when available. Considering the short half-life and data quality, we further designed a neutrophil sorting protocol (Figure S1E) and an in silico strategy (STAR Methods). We successfully sequenced 103 samples from 64 patients that passed quality control (Table S1), including primary and metastatic samples as well as matched normal tissues and blood (STAR Methods). We further applied the standardized pipeline on published datasets and finally generated a neutrophil map of 225 samples from 143 patients across 17 cancer types, among which 12 cancer type data (70.59%) were newly generated or in-house (Figure 1C). After harmonizing the data batches, excluding the low-quality cells, and balancing the RNA dropouts, 1,79,908 single neutrophils finally passed quality control, of which 79.29% data were in-house or freshly released. Together, our single-cell profiling constitutes a potential resource for neutrophil investigation (available at http://www.pancancer.cn/neu/).
根据我们的初步发现,我们设计了一种标准化采样策略,重点关注中性粒细胞浸润程度较高或中等的癌症类型(样本统计数据,图 S1D),确保包含匹配的血液、邻近正常组织和可用的转移样本。考虑到半衰期短和数据质量,我们进一步设计了中性粒细胞分选方案(图 S1E)和计算机策略(STAR 方法)。我们成功对来自 64 名患者的 103 个样本进行了测序,并通过了质量控制(表 S1),包括原发性和转移性样本以及匹配的正常组织和血液(STAR 方法)。我们进一步将标准化流程应用于已发布的数据集,最终生成了来自 17 种癌症类型的 143 名患者的 225 个样本的中性粒细胞图谱,其中 12 种癌症类型数据 (70.59%) 是新生成的或内部的(图 1C)。经过协调数据批次、剔除低质量细胞、平衡RNA dropout等后,最终有1,79,908个中性粒细胞通过质控,其中79.29%的数据为内部数据或新发布的数据。总之,我们的单细胞分析构成了中性粒细胞研究的潜在资源(可在http://www.pancancer.cn/neu/获取)。
Transcriptional signatures across cancers
跨癌症的转录特征
To decode the transcriptional signature, we clustered neutrophils and noted high heterogeneity across cancer and tissue types (Figures 1D and S1F). We observed 10 distinct states composed of S100A12+, HLA-DR+CD74+, VEGFA+SPP1+, TXNIP+, CXCL8+IL1B+, CXCR2+, IFIT1+ISG15+, MMP9+, NFKBIZ+HIF1A+, and ARG1+ neutrophils (Figures 1E and 1F; Table S2). For example, HLA-DR+CD74+ subset showed high expression of major histocompatibility complex (MHC)-II molecules and universal infiltration across cancers. Conforming to the neutrophil biology,15 we also identified clusters potentially representing inflammatory response (CXCL8+IL1B+) and specific chemotaxis (CXCR2+) features (Figure 1F). By decoding neutrophil transcriptome into transcriptional programs (STAR Methods), we consistently confirmed their featured activation modes such as chemotaxis or inflammation (Figure S1G). Given the neutrophil single-cell profiles were only reported in certain cancer types, we computed the correlation with previously defined neutrophil states and subpopulations (STAR Methods) of 5 independent studies.7,22,23,24,25 Some subsets of our data showed strong consistency with published states such as IFIT1+ISG15+ and hNeutro2 (Figure S1H). Together, our neutrophil map not only captured known neutrophil subpopulation features but also revealed potentially uncharacterized neutrophil subsets.
为了解码转录特征,我们对中性粒细胞进行了聚类,并注意到癌症和组织类型之间的高度异质性(图 1D 和 S1F)。我们观察到由 S100A12 + 、 HLA-DR + CD74 + 、 VEGFA + SPP1 + 、TXNIP + 、CXCL8 + IL1B + 、CXCR2 + 、IFIT1 + ISG15 + 、NFKBIZ + HIF1A + 和 ARG1 + 中性粒细胞(图 1E 和 1F;表 S2)。例如,HLA-DR + CD74 + 子集显示主要组织相容性复合体 (MHC)-II 分子的高表达和跨癌症的普遍浸润。符合中性粒细胞生物学, 15 我们还鉴定了可能代表炎症反应(CXCL8 + IL1B + )和特定趋化性(CXCR2 + 22 23 24 25 我们的数据的一些子集与已发布的状态(例如 IFIT1)显示出很强的一致性< b26> ISG15 + 和 hNeutro2(图 S1H)。总之,我们的中性粒细胞图不仅捕获了已知的中性粒细胞亚群特征,而且还揭示了潜在的未表征的中性粒细胞亚群。
Molecular divergence and survival correlation
分子分歧和生存相关性
To explore the principles governing neutrophil clades, we computed the tree structure26 and observed the scattered tree leaves of distinct subsets (Figure 2A). Based on Ro/e analysis (ratio of observed cell number to expected cell number) (Figure 2B), HLA-DR+CD74+ and VEGFA+SPP1+ neutrophils were the most cancer-enriched subsets overall (Figures 2B and S2A) but showed cancer type preferences. HLA-DR+CD74+ neutrophils were enriched in non-small cell lung cancer (NSCLC), bladder cancer (BLCA), and ovarian cancer (OV) while showing decreased infiltration in renal cell carcinoma (RCC) and oral squamous cell carcinoma (OSCC). In contrast, VEGFA+SPP1+ neutrophils showed sparse infiltration in NSCLC, BLCA, and OV but were enriched in RCC and stomach adenocarcinoma (STAD). These data were partly validated by flow cytometry (n = 24, see Table S1) and multiplex immunohistochemistry (mIHC) using an independent multi-cancer-TMA cohort (n = 68, see Table S1) comprising 8 cancer types, including breast invasive carcinoma (BRCA), colon adenocarcinoma (COAD), intrahepatic cholangiocarcinoma (ICC), hepatocellular carcinoma (HCC), STAD, NSCLC, RCC, and pancreatic adenocarcinoma (PAAD) (Figures 2C–2E and S2B).
为了探索控制中性粒细胞分支的原理,我们计算了树结构 26 并观察了不同子集的分散树叶(图 2A)。基于 Ro/e 分析(观察到的细胞数与预期细胞数的比率)(图 2B),HLA-DR + CD74 + 和 VEGFA + SPP1 + 中性粒细胞总体上是癌症最丰富的亚群(图 2B 和 S2A),但表现出癌症类型偏好。 HLA-DR + CD74 + 中性粒细胞在非小细胞肺癌 (NSCLC)、膀胱癌 (BLCA) 和卵巢癌 (OV) 中富集,同时在肾癌中表现出浸润减少细胞癌(RCC)和口腔鳞状细胞癌(OSCC)。相反,VEGFA + SPP1 + 中性粒细胞在 NSCLC、BLCA 和 OV 中表现出稀疏浸润,但在 RCC 和胃腺癌 (STAD) 中富集。这些数据通过流式细胞术(n = 24,参见表 S1)和多重免疫组织化学 (mIHC) 使用独立的多癌症 TMA 队列(n = 68,参见表 S1)进行了部分验证,该队列包括 8 种癌症类型,包括乳腺癌浸润性癌(BRCA)、结肠腺癌 (COAD)、肝内胆管癌 (ICC)、肝细胞癌 (HCC)、STAD、NSCLC、RCC 和胰腺腺癌 (PAAD)(图 2C-2E 和 S2B)。
To explore the association of neutrophil subsets with patient survival, we first analyzed the neutrophil subset signatures based on the TCGA pan-cancer dataset. Among the subsets, VEGFA+SPP1+ subset was linked with the worst patient outcome (Figures 2F and S2C; 8-cancer-TMA cohort), whereas HLA-DR+ was linked with the best. We further corroborated the prognostic relevance of HLA-DR+ neutrophils through mIHC in an independent 8-cancer-TMA cohort (Figure 2G, n = 1,116; HCC, n = 357; COAD, n = 93; NSCLC, n = 90; STAD, n = 85; RCC, n = 150; OV, n = 160; BRCA, n = 129; BLCA, n = 52). These data indicated that HLA-DR+ neutrophils may represent a key anti-tumor neutrophil subset across a majority of cancer types.
为了探索中性粒细胞亚群与患者生存的关联,我们首先根据 TCGA 泛癌数据集分析了中性粒细胞亚群特征。在这些子集中,VEGFA + SPP1 + 子集与最差的患者结果相关(图 2F 和 S2C;8-cancer-TMA 队列),而 HLA-DR + 中性粒细胞的预后相关性(图 2G,n = 1,116;HCC,n = 357;COAD,n = 93;NSCLC ,n = 90;RCC,n = 150;BRCA,n = 129;这些数据表明 HLA-DR + 中性粒细胞可能代表大多数癌症类型中关键的抗肿瘤中性粒细胞亚群。
We next compared cytokine profiles of distinct neutrophil subsets (Figure 2H). HLA-DR+CD74+ neutrophils showed high CCL5, which can recruit T cells.27 IFIT1+ISG15+ neutrophils were associated with high expression of PD-L1 (CD274), indicating its immunosuppressive role. Only HLA-DR+CD74+ neutrophils showed specific enrichment of MHC class II molecules such as HLA-DRA and HLA-DRB1 (Figure 2I). However, almost all neutrophil subsets expressed high levels of MHC class I molecules, aligning with their consensual expression on all nucleated cells.28 In parallel, the immunophenotypic signature29 of neutrophil subsets showed remarkable diversity (Figure 2J), whereas almost all subsets showed high aging signature, supporting the notion that tumor-associated neutrophils are mostly in a maturation state.2 Reciprocally, we also identified the ubiquitous cell cycle gene expression among almost all subsets and the enrichment of interferon-related genes in IFIT1+ISG15+ neutrophils (Figure S2D). Given the known circadian features of neutrophil biology,30,31 we further compared the phenotypic signatures sampled at different times of day and observed higher maturation and chemotaxis levels of neutrophils at daytime (Figure S2E). These data collectively highlighted the subset-specific molecular hallmark of neutrophils.
接下来,我们比较了不同中性粒细胞亚群的细胞因子谱(图 2H)。 HLA-DR + CD74 + 中性粒细胞显示高 CCL5,可以招募 T 细胞。 27 IFIT1 + ISG15 + 中性粒细胞与 PD-L1 (CD274) 高表达相关,表明其具有免疫抑制作用。仅 HLA-DR + CD74 + 中性粒细胞表现出 MHC II 类分子(例如 HLA-DRA 和 HLA-DRB1)的特异性富集(图 2I)。然而,几乎所有中性粒细胞亚群都表达高水平的 MHC I 类分子,这与它们在所有有核细胞上的一致表达一致。 28 与此同时,中性粒细胞亚群的免疫表型特征 29 显示出显着的多样性(图 2J),而几乎所有亚群都表现出高衰老特征,这支持了肿瘤相关中性粒细胞大多是处于成熟状态。 2 相反,我们还鉴定了几乎所有亚群中普遍存在的细胞周期基因表达以及 IFIT1 + ISG15 + 中性粒细胞中干扰素相关基因的富集(图S2D)。鉴于中性粒细胞生物学的已知昼夜节律特征, 30 31 我们进一步比较了一天中不同时间采样的表型特征,并观察到白天中性粒细胞的成熟度和趋化性水平较高(图 S2E) 。这些数据共同突出了中性粒细胞的子集特异性分子标志。
Given neutrophils’ capacity against pathogens, we compared cell states between cancer and inflammatory conditions. As a result, IFIT1+ISG15+ neutrophils showed an expanded proportion in chronic pancreatitis and cholecystitis, which also exhibited high infiltration in pancreatic and gallbladder cancers (Figure S2F). Likewise, neutrophils in COVID-19 lung samples showed a specific spectrum of states, such as enrichment of inflammatory NFKBIZ+HIF1A+ subsets, which was also observed in lung cancer. This result may imply the shared tissue-restricted reprogramming across diverse disease conditions or a causal resemblance between inflammation and cancer in indicated organs. In contrast, certain subsets are highly tumor-specific (i.e., HLA-DR+CD74+ and VEGFA+SPP1+ subsets), indicating the TME-specific stimuli on reprogramming neutrophil states. Taken together, our data revealed the potential preference of neutrophil infiltration (i.e., HLA-DR+CD74+ and VEGFA+SPP1+ subsets) into certain cancer types and indicated the existence of cancer-imprinted transcriptional programs.
鉴于中性粒细胞对抗病原体的能力,我们比较了癌症和炎症条件下的细胞状态。结果,IFIT1 + ISG15 + 中性粒细胞在慢性胰腺炎和胆囊炎中表现出扩大的比例,在胰腺癌和胆囊癌中也表现出高浸润(图 S2F)。同样,COVID-19 肺部样本中的中性粒细胞显示出特定的状态谱,例如炎症 NFKBIZ + HIF1A + 亚群的富集,这在肺癌中也观察到。这一结果可能意味着不同疾病条件下共享的组织限制性重编程或指定器官中炎症和癌症之间的因果相似性。相反,某些亚群具有高度肿瘤特异性(即 HLA-DR + CD74 + 和 VEGFA + SPP1 + 亚群) ,表示对中性粒细胞状态重编程的 TME 特异性刺激。综上所述,我们的数据揭示了中性粒细胞浸润的潜在偏好(即 HLA-DR + CD74 + 和 VEGFA + SPP1 + 子集)进入某些癌症类型,并表明癌症印记转录程序的存在。
Maturation and metabolism states
成熟和代谢状态
Neutrophils have long been regarded as mature and terminally differentiated cells. However, it remains unknown how the diversity of their maturation states is achieved. To address this question, we applied a vector-field-based deep learning algorithm to infer the progressive steps of lineage specification and divergence.32 We observed continuous differentiation along neutrophil states, with the most terminal pseudotime value observed for HLA-DR+CD74+ neutrophils (Figures 3A and S3A); this pattern was replicated by other algorithms, including monocle3,33 CytoTRACE,34 and Slingshot35 (Figures 3A and S3A). We ranked neutrophils according to both their tumor specificity and pseudotime, finding that HLA-DR+CD74+ neutrophils potentially remained at the terminus (Figure 3B). To validate this finding, we first evaluated the maturation markers CD11b and CD1636 using flow cytometry (Figure 3C) in intratumor neutrophils from 24 patients with 8 cancer types (Table S1). As a result, the HLA-DR+ neutrophil subset showed significantly higher CD11b and CD16 expression. We subsequently confirmed the enhanced maturation markers’ (CD11b and MPO)36 expression in HLA-DR+ neutrophils using mIHC across cancers (Figures 3D and 3E). Also, transcription factor RFX5 showed specific activation among the HLA-DR+ subset, consistent with ChIP-seq and knockdown/overexpression assays (Figures S3B–S3E; STAR Methods). Taken together, these data indicated that HLA-DR+ neutrophils were potentially one of the terminally mature neutrophil subsets.
中性粒细胞长期以来被认为是成熟且终末分化的细胞。然而,目前尚不清楚它们的成熟状态的多样性是如何实现的。为了解决这个问题,我们应用了基于向量场的深度学习算法来推断谱系规范和分歧的渐进步骤。 32 我们观察到沿中性粒细胞状态的连续分化,观察到 HLA-DR + CD74 + 中性粒细胞的最终伪时间值(图 3A 和 S3A);其他算法也复制了这种模式,包括 monocle3、 33 CytoTRACE、 34 和 Slingshot 35 (图 3A 和 S3A)。我们根据中性粒细胞的肿瘤特异性和假时间对中性粒细胞进行排名,发现 HLA-DR + CD74 + 中性粒细胞可能保留在终点(图 3B)。为了验证这一发现,我们首先使用流式细胞术(图 3C)评估了来自 24 名患有 8 种癌症类型的患者的瘤内中性粒细胞的成熟标记物 CD11b 和 CD16 36 (表 S1)。结果,HLA-DR + 中性粒细胞亚群表现出显着较高的 CD11b 和 CD16 表达。随后,我们使用 mIHC 跨癌症证实了 HLA-DR + 中性粒细胞中成熟标记物(CD11b 和 MPO) 36 表达的增强(图 3D 和 3E)。此外,转录因子 RFX5 在 HLA-DR + 亚群中表现出特异性激活,与 ChIP-seq 和敲低/过表达测定一致(图 S3B–S3E;STAR 方法)。总而言之,这些数据表明 HLA-DR + 中性粒细胞可能是最终成熟的中性粒细胞亚群之一。
A critical question is how pathway activity varies across tumor neutrophil subsets. To address this issue, we measured the pathway variance and observed the strong diversity of metabolic pathways (Figure S3F), which supports the potential metabolic regulation of neutrophil identity maintenance. Then, we quantified the metabolic pathway activity of the neutrophil subsets (Figures 3F and S3G; Table S3). Notably, HLA-DR+ neutrophils showed remarkable enrichment of amino acid metabolism (i.e., valine, leucine, and isoleucine, Figure 3F). In parallel, the activation of vitamin metabolism and glycan metabolism were dominant among the immunosuppressive VEGFA+SPP1+ neutrophils. Together, these results supported the possibility that tumor neutrophils were metabolically coordinated, raising the idea that amino acid metabolism primes the HLA-DR programs.
一个关键问题是肿瘤中性粒细胞亚群的通路活性如何变化。为了解决这个问题,我们测量了通路方差并观察到代谢通路的强烈多样性(图S3F),这支持中性粒细胞身份维持的潜在代谢调节。然后,我们量化了中性粒细胞亚群的代谢途径活性(图 3F 和 S3G;表 S3)。值得注意的是,HLA-DR + 中性粒细胞显示出氨基酸代谢显着富集(即缬氨酸、亮氨酸和异亮氨酸,图 3F)。同时,维生素代谢和聚糖代谢的激活在免疫抑制性 VEGFA + SPP1 + 中性粒细胞中占主导地位。总之,这些结果支持了肿瘤中性粒细胞代谢协调的可能性,提出了氨基酸代谢启动 HLA-DR 程序的观点。
Leucine metabolism governs the epigenetics of the antigen-presenting machinery
亮氨酸代谢控制抗原呈递机制的表观遗传学
To systematically examine the effect of amino acids on neutrophils, we designed an in vitro screening strategy comprising all 20 amino acids and investigated their impact on antigen presentation (Figure 4A) in circulating neutrophils from healthy donors (STAR Methods). Particularly, leucine upregulated HLA-DR (Figure 4B) and costimulatory molecules such as CD80 (Figure 4C). Although arginine slightly upregulated HLA-DR, it cannot impact on costimulatory molecules (Figure S4A). We further expanded our analysis of leucine on the spectrum of antigen presentation processes by using PCR array analysis. Notably, leucine significantly promoted the gene expression of MHC-II complex assembly (Figure 4D), enhanced antigen processing protease (Figure 4E), and facilitated antigen-loading processes (Figure 4F). We confirmed that intracellular leucine levels did increase upon leucine treatment and that these effects were not due to contamination with other antigen-presenting cells (Figure S4B). RNA-seq analyses also confirmed the impact of leucine on MHC-II but not MHC-I (Figures 4G, S4C, and S4D). We also observed that leucine increased the in vitro survival of neutrophils (Figure S4E). Among the matched clinical samples, the leucine concentration also showed a strong positive correlation with the HLA-DR+ neutrophil signature (Figure S4F; STAR Methods). Collectively, these data supported the notion that leucine potently primes the neutrophil antigen-presenting program.
为了系统地检查氨基酸对中性粒细胞的影响,我们设计了一种包含所有 20 种氨基酸的体外筛选策略,并研究了它们对健康供体循环中性粒细胞中抗原呈递的影响(图 4A)(STAR 方法)。特别是,亮氨酸上调 HLA-DR(图 4B)和共刺激分子,例如 CD80(图 4C)。虽然精氨酸略微上调 HLA-DR,但它不会影响共刺激分子(图 S4A)。通过使用 PCR 阵列分析,我们进一步扩大了对抗原呈递过程范围的亮氨酸分析。值得注意的是,亮氨酸显着促进 MHC-II 复合体组装的基因表达(图 4D),增强抗原加工蛋白酶(图 4E),并促进抗原加载过程(图 4F)。我们证实,亮氨酸处理后细胞内亮氨酸水平确实增加,并且这些影响不是由于其他抗原呈递细胞的污染造成的(图 S4B)。 RNA-seq 分析还证实了亮氨酸对 MHC-II 的影响,但对 MHC-I 没有影响(图 4G、S4C 和 S4D)。我们还观察到亮氨酸增加了中性粒细胞的体外存活率(图 S4E)。在匹配的临床样本中,亮氨酸浓度也显示出与 HLA-DR + 中性粒细胞特征的强正相关性(图 S4F;STAR 方法)。总的来说,这些数据支持亮氨酸有效启动中性粒细胞抗原呈递程序的观点。
To probe how leucine feeds the metabolome and drives the antigen-presenting program, we performed the untargeted metabolomics and observed strong differences in ATP and fatty acid production (Figures 4H, S4G, and S4H; Table S4), partly congruent with the literature.38 Considering that ATP generation mainly occurs in the mitochondria,39,40 we hypothesized that leucine may trigger functional or phenotypic remodeling via mitochondria. Indeed, leucine administration caused mitochondria aggregation and altered phenotypes (i.e., mitochondrial quality, Ca+, and ROS production) (Figures 4I and S4I–S4L). To test the dynamic impact of leucine on mitochondria, we examined the real-time changes in oxygen consumption rate (OCR) and observed that leucine stimulation significantly augmented the mitochondrial OCR (Figure 4J). Transmission electron microscopy (TEM) analysis showed that leucine treatment could induce specific morphological features such as longer mitochondrial length and more pseudopods on the membrane (Figures 4K and S4M), indicating the increased intercellular contact potential.41 These data were partly in line with the reported role of mitochondrial metabolism in regulating professional antigen-presenting cells.42
为了探究亮氨酸如何供给代谢组并驱动抗原呈递程序,我们进行了非靶向代谢组学,并观察到 ATP 和脂肪酸产生的巨大差异(图 4H、S4G 和 S4H;表 S4),与文献部分一致。 38 考虑到ATP的生成主要发生在线粒体中, 39 40 我们假设亮氨酸可能通过线粒体触发功能或表型重塑。事实上,亮氨酸的施用会导致线粒体聚集并改变表型(即线粒体质量、Ca + 和 ROS 产生)(图 4I 和 S4I–S4L)。为了测试亮氨酸对线粒体的动态影响,我们检查了耗氧率 (OCR) 的实时变化,并观察到亮氨酸刺激显着增强了线粒体 OCR(图 4J)。透射电子显微镜 (TEM) 分析表明,亮氨酸处理可以诱导特定的形态特征,例如更长的线粒体长度和膜上更多的伪足(图 4K 和 S4M),表明细胞间接触电位增加。 41 这些数据部分符合报道的线粒体代谢在调节专业抗原呈递细胞中的作用。 42
Given the complexity of the mitochondrial respiration and electron transport chain machinery, we explored the causal effect of leucine on specific mitochondrial subcomponents. By quantifying the mitochondrial respiration signature of single HLA-DR+ neutrophils (Figure 4L), we observed that complex I, which is capable of transferring electrons from reduced nicotinamide adenine dinucleotide (NADH) and producing nicotinamide adenine dinucleotide (NAD),43 showed higher activity upon leucine administration. However, complex III and IV signature activity was weak or moderate. This observation was supported by the stronger output of NAD by leucine-treated neutrophils (Figure S4N). We therefore inhibited the complex I activity and observed decreased mitochondrial membrane state and HLA-DR+ neutrophil proportion (Figures 4M and S4O). Other mitochondria respiration inhibitors also showed coherent results (Figure S4P). Conversely, NAD supplementation caused higher HLA-DR intensity (Figure S4Q), further confirming the mitochondrial respiration-dependent function of leucine in the HLA-DR program.
鉴于线粒体呼吸和电子传递链机制的复杂性,我们探讨了亮氨酸对特定线粒体亚成分的因果影响。通过量化单个 HLA-DR + 中性粒细胞的线粒体呼吸特征(图 4L),我们观察到复合物 I 能够从还原型烟酰胺腺嘌呤二核苷酸 (NADH) 转移电子并产生烟酰胺腺嘌呤二核苷酸( NAD), 43 在亮氨酸给药后表现出更高的活性。然而,复合物 III 和 IV 特征活性较弱或中等。亮氨酸处理的中性粒细胞更强的 NAD 输出支持了这一观察结果(图 S4N)。因此,我们抑制了复合物 I 的活性,并观察到线粒体膜状态和 HLA-DR + 中性粒细胞比例降低(图 4M 和 S4O)。其他线粒体呼吸抑制剂也显示出一致的结果(图 S4P)。相反,补充 NAD 会导致更高的 HLA-DR 强度(图 S4Q),进一步证实了 HLA-DR 程序中亮氨酸的线粒体呼吸依赖性功能。
We further asked how leucine was catabolized and enhanced HLA-DR expression. We initially fed neutrophils from healthy donor blood with 13C-labeled leucine and found that leucine was catalyzed into acetyl-CoA, entered the tricarboxylic acid (TCA) cycle, and generated an increased amount of glutamate and glutamine (Figure 4N), using the reported catabolism route.44 Consistently, CoA biosynthesis signature was primed during leucine treatment (Figure S4R), and the acetyl-CoA activator significantly upregulated HLA-DR (Figure 4O). Acetyl-CoA inhibition significantly reduced HLA-DR level while its restoration rescued HLA-DR, indicating the acetyl-CoA-dependent regulation on HLA-DR. Given the known link between acetyl-CoA and histone H3 lysine 27 acetylation (H3K27ac),45,46,47 we examined the total H3K27ac level and observed the enhanced H3K27ac upon leucine treatment (Figure 4P), without changes in histones H3 themselves (Figure S4S). Further, CUT&Tag showed significantly upregulated H3K27ac on MHC-II genes (Figure 4Q), but not H3K27me3 and H3K4me3 (Figures 4R and 4S). These data supported the notion that leucine can impact H3K27ac and thereby activate MHC-II genes (i.e., HLA-DRA and HLA-DQB1), its transcription factor, and its regulatory element (i.e., MHC-II super-enhancer48) (Figures 4T and S4T). Consistently, the chromatin accessibility of MHC-II genes was also enhanced by leucine treatment (Figure S4U). Overall, our data pointed to the dependency of leucine catabolism for antigen-presenting machinery initiation through mitochondria modeling and metabolism-epigenetic regulation such as the acetyl-CoA/H3K27ac/MHC-II axis.
Antigen-presenting neutrophils spatially link with and fuel T cell responses
Because HLA-DR+ neutrophils favor prognosis, we asked about the mechanisms underlying their potential anti-tumor effects. We first performed bulk RNA-seq of matched tumor samples, deconvoluted18 the immune cell profile, clustered the immune cell proportions, and observed patterned neutrophil-T cell infiltration profiles (Figures 5A and S5A). In particular, HLA-DR+ neutrophils co-localized with a broad spectrum of anti-tumor T cell subsets (i.e., CD4+ effector memory T cells, CD8+ effector memory T cells, and CD8+ central memory T cells). We fetched 50 spatial transcriptomics datasets covering 1,78,330 spots derived from 9 cancer types (BRCA, SKCM, CESC, COAD, OV, PRAD, LGG, HCC, and RCC, Table S5) and computed the correlation between HLA-DR+ neutrophils and major immune lineages (STAR Methods). CD8+ T cells and CD4+ T cells ranked highly among the major lineages (Figure S5B). Of note, in RCC samples (GSM5924040), CD8+ T cells showed strong co-localization with HLA-DR+ neutrophils, with similar results observed in OV (10x) and colorectal cancer (CRC) (OEP001756) samples (Figure 5B). These observations implied that antigen-presenting neutrophils are spatially linked with T cells.
To decode the effect of HLA-DR+ neutrophils on T cells, we cocultured neutrophil subsets (sorted from the tumor) with autologous CD3+ T cells (sorted from PBMCs) for 3 days, indicating that HLA-DR+ neutrophils can promote T cells to express TNFα (Figure 5C). Given their capacity for antigen presentation, we then quantified the allele-specific HLA gene expression at the single-cell level (Figures 5D and S5C). HLA-DRB1, HLA-DPA1, and HLA-DQB1 showed high frequency, where these loci were reported with T cell epitope restricted reactiveness (HLA-DRB1∗09:01)49 or cytokine production (HLA-DPB1∗05:01).50 We further confirmed the antigen uptake ability of HLA-DR+ neutrophils using OVA fluorescent proteins (Figures S5D and S5E). We hence fed the neutrophils with human MHC-II antigens (gp100, 44–59; CMV, 65–71) and cocultured them with autologous T cells (Figure 5E). Notably, sorted HLA-DR+ neutrophils gave rise to the antigen-specific response of autologous T cells, although more weakly than positive controls (DCs). These data naturally raise the possibility that HLA-DR+ neutrophils are likely to present tumor neoantigens and elicit reactive T cell responses. To test this hypothesis, HLA-DR+ neutrophils were incubated with mutation-derived neoantigens (TP53, KRAS, IDH2, and BAP1)51,52,53 for 24 h and were further cocultured with autologous CD3+ T cells (Figure 5F). We observed the stimulated T cell responses of most of the neoantigens, although the responses were variable across donors. For example, KRASG12V neoantigen (MTEYKLVVVGAVGVGKSALTIQLI)53 was associated with strong neutrophil-triggered T cell activation in donors 2 and 3. Our reported KRASG12D neoantigen (LVVVGADGV)52 also led to reactive T cell responses. Focusing on these two peptides, we performed coculture in which the ratio of neutrophil to T cells and the concentration of peptides were controlled in a donor with HLA−DRB1∗07:01 and HLA−A∗02:01 types (Figure 5G). When the ratio of HLA-DR+ neutrophils to T cells was 10:1, a potent neoantigen response of T cells was activated with KRASG12V or KRASG12D neoantigens. By sequencing the TCR repertoire of KRASG12V neoantigen-activated T cells, we observed strong reactive TCR gene rearrangement with HLA-DR+ neutrophil stimulation (negative control, T cell alone; positive control, autologous DCs) (Figures 5H and S5F), which was comparable with DCs (Figure 5I). We finally authenticated the co-localization of reactive T cells (CD39+CXCL13+CD4+ T cells and CD39+CXCL13+CD8+ T cells) and antigen-presenting neutrophils (HLA-DR+CD15+ neutrophils) in the multi-cancer-TMA cohort (Figures 5J and S5G). Together, these data indicated that antigen-presenting neutrophils can effectively generate reactive T cell responses.
We next asked whether neutrophil-elicited T cell activation is dependent on the antigen. By coculturing leucine-stimulated neutrophils with autologous T cells from healthy donors (removing leucine after 24 h of neutrophil stimulation), we observed consistent but slightly weaker T cell activation (Figure S5H) and killing ability for cell lines derived from multiple cancers (Figure S5I), raising the hypothesis that HLA-DR+ neutrophils activate T cells in a nonspecific manner. Subsequently, we tested distinct coculture methods (i.e., in a medium, in a transwell chamber, or directly) and examined the T cell activation levels (Figure S5J). Interestingly, only direct coculture was associated with T cell activation, whereas coculture via medium or transwell system was not. We also neutralized neutrophil-related cytokines but did not observe T cell immunophenotype changes (Figure S5K). These results indicated that HLA-DR+ neutrophils may directly activate T cells mainly via ligand-receptor interactions.
To identify the ligand that drives T cell activation, we ranked the in silico results according to interaction frequency and gene expression percentage and focused on the candidate ligands such as ICAM1 and CXCL10 (Figures S5L and S5M). Upon inhibiting the ICAM1 and its interaction with ITGAL, we observed significantly reduced T cell activation (Figure S5N). Meanwhile, ICAM1 and HLA-DR were co-expressed and co-upregulated by leucine (Figure S5O). In contrast, CXCL10 inhibition showed a negligible effect on T cell activation (Figure S5N). In summary, HLA-DR+ neutrophils are associated with an active TME and can broadly trigger T cell activation, (neo)antigen reactivity, and cytotoxicity, raising the possibility that these cells could be delivered to fuel T cell responses.
Exploring neutrophil-based immunotherapy to fire up TME
Another critical question is whether HLA-DR+ neutrophils have the potential to enhance immunotherapy in vivo. To this end, we first evaluated the conservation of neutrophil subsets within mouse TME by collecting mouse scRNA-seq data from 7 murine cancer types (Table S5). We observed similar subsets such as Cd74+, Spp1+Vegfa+, and Isg15+ neutrophils in mice (Figure 6A) as recently reported (Figure S6A).7,11,14 We further observed H2-Aa (HLA-DQA homolog), H2-Ab1 (HLA-DQB1 homolog), Cd74 (CD74 homolog) expression, and antigen-presenting signature in Cd74+ subsets (Figure 6B) and confirmed this subset using mIHC (Figure 6C). These results not only highlighted the conserved role of antigen-presenting neutrophils but also hinted at the opportunities for investigating these cell subsets in vivo.
Given the role of leucine demonstrated above, we stimulated mouse circulating neutrophils with leucine in vitro and observed the upregulation of Cd74, Cd80, and Cd86 (Figures 6D and S6B). However, circulating neutrophils from knockout mice lacking the leucine transporter or catabolism enzymes (Lat1KO, Bcat2KO, and DbtKO) did not respond to leucine treatment (Figure 6D), again supporting the role of leucine in MHC-II processes.
We subsequently investigated the in vivo association between antigen-presenting neutrophils and tumor phenotypes. Initially, we generated the MHC-IIflox/flox; Ly6GCre-tdTomato mice with specific deletion of antigen-presenting neutrophils (Figures 6E and S6C), subcutaneously injecting murine cancer cells (LLC, MC38, and Hepa 1–6), and observed increased tumor growth. Interestingly, intratumor Cd4 and Cd8 T cells both showed significantly decreased infiltration in MHC-IIflox/flox; Ly6GCre-tdTomato group (Figure 6F). We next attempted to increase the antigen-presenting neutrophils in vivo by giving mice the leucine-rich diet (1.5% leucine in water) and observed increased Cd74+ neutrophils (Figure 6G), but this short-term diet did not impact cancer volume or body weight (Figures S6D and S6E). We also tested other amino acids (arginine, cysteine, glutamine, and tryptophan) but did not find consistently increased Cd74+ neutrophils (Figure S6F).
We further analyzed the TME profile of the leucine diet using scRNA-seq (Figure S6G), which confirmed the expanded Cd74+ proportion upon leucine treatment (Figure S6H). In addition to antigen-presenting, leucine-stimulated neutrophils also exhibited upregulated chemokines and Toll-like receptor pathways (Figure S6I). Within the cancer cells themselves (Figure S6J), the leucine metabolism signature was also upregulated, along with the robust activity of the epithelial-mesenchymal-transition and TNFα. Given that CD74+ neutrophils can activate T cells, we examined T cell infiltration and observed their increased infiltration (Figures 6H and 6I). Together, these data indicated that a short-term leucine diet was beneficial for anti-TME and induced mild changes in cancer cell phenotypes, raising its potential therapeutic usage.
Then, we explored the therapeutic effect of Cd74+ neutrophils in murine cancers. Notably, leucine diet plus anti-PD-1 therapy significantly reduced tumor volumes and enabled the achievement of stable disease (Figure 7A). Conversely, neutrophil deletion diminished the efficacy of this combination therapy (Figure S7A). We also investigated the delivery of antigen-presenting neutrophils into the TME as another therapeutic option. We isolated mouse circulating neutrophils, stimulated them using leucine, and confirmed their Cd74 upregulation. Delivering those Cd74+ neutrophils into the tumors significantly reduced tumor size but still did not generate stable disease outcomes (Figure 7B), even when the number of neutrophils was increased to 1 × 107 (Figure S7B). When transferring Cd74-deficient neutrophils, we found significantly diminished anti-tumor effects (Figure 7B). We also assessed the lifespan of the antigen-presenting neutrophils upon delivering Cd45.1 neutrophils into Cd45.2 mice and estimated their half-life to be potentially longer (Figure S7C). Upon combining anti-PD-1 antibody and Cd74+ neutrophil delivery (Figure 7C), we observed robust anti-tumor efficacy in all tumor models. Notably, in the MC38 and Hepa 1–6 models, a significant proportion of tumors showed complete responses (MC38, 4 of 10; Hepa 1–6, 6 of 10). In contrast, Cd74-knockout neutrophils showed weak efficacy in combination with anti-PD-1.
We finally examined antigen-presenting neutrophils in clinical immunotherapy data. In 8 cohorts receiving immunotherapy spanning SKCM, STAD, HCC, BLCA, and NSCLC patients,54,55,56,57,58,59,60,61 we observed significant positive correlations of the HLA-DR+ neutrophil signature with better survival or responses (Figures 7D and 7E). Next, we tested the adoptive delivery ex vivo in anti-PD-1 immunotherapy-resistant HCC samples following the patient-derived tumor fragment (PDTF) strategy.62 After 3 days of coculture, both CD4 and CD8 T cells displayed upregulation of cytotoxic molecules (IFNγ and TNFα) as well as reactive molecules (4-1BB and CD39) (Figure 7F). Combing antigen-presenting neutrophils and PD-1 antibodies also generated greater T cell reactiveness and cytotoxicity (Figure 7G). Altogether, these in vivo and human-sample-derived data highlighted the synergy of antigen-presenting neutrophils in immunotherapy and suggested therapeutic opportunities such as adoptive delivery.
Discussion
High-dimensional single-cell profiling has revolutionized cancer immunology with its scale and ability to decode complex microenvironment63 and is leading to organ-wide profiling of certain cell types such as T cells or macrophages.64,65 For neutrophils, the opportunity is now here to orchestrate a unified data-driven framework for defining neutrophil ontology. Here, we integrated the neutrophil transcriptomes of 225 samples from 143 patients across 17 cancer types and observed high transcriptional heterogeneity composed of 10 cell states. Notably, the HLA-DR program could be stimulated by leucine treatment, which induced higher mitochondrial respiration, acetyl-coenzyme A output, and epigenetic activation of antigen-presenting genes. As such, a leucine diet or adoptive delivery strategy could boost cancer immunotherapeutic efficacy and may serve as a potential neutrophil-based therapeutic strategy.
Profiling single human neutrophils in an unbiased manner remains a major challenge. Both mouse and human neutrophils are short-lived3,66 and hence have been largely neglected in previous cancer single-cell profiling studies. Our strategy of neutrophil sorting plus single-cell sequencing has made it possible to explore the spectra of neutrophil states, the heterogeneity of cell subsets, and their temporal changes throughout cancer stages. Our cross-species integration of neutrophil data could help establish a hierarchy of neutrophils that is conserved in humans and mice, providing a potential foundation for the application of murine models in exploring neutrophil-based therapies.
The cancer-immune cycle is initiated by antigen presentation processes67; however, impaired antigen-presenting machineries within the TME, such as DC dysfunction and HLA downregulation, can drive immune evasion.68 Our high-resolution mapping approach identified HLA-DR+CD74+ neutrophils as alternative antigen-presenting cells across various cancer types. A similar subset of neutrophils was discovered in early-stage lung cancer in 2016,69 and an independent group recently reported that an anti-FcγRIIIB-antigen conjugate may convert neutrophils into antigen-presenting cells.70 In line with this functional plasticity, our work provides further evidence of metabolic regulation of antigen presentation, demonstrating the role of amino acids in promoting neutrophil antigen presentation, whereas fatty acid or glucose nutrient stress was essential in regulating macrophage antigen-presenting function.71,72 Compared with other professional antigen-presenting cells like dendritic cells, B cells, and certain macrophage subsets, neutrophils expressing HLA-DR may have some advantages in certain contexts. For example, neutrophils are often one of the first responders to sites of inflammation.1 Uniquely, neutrophils possess active phagocytic capabilities and chemotaxis features.1 This attribute potentially empowers them to efficiently migrate and degrade antigens, thus bolstering their role as antigen-presenting cells. Another crucial feature of neutrophils is their short half-life, which potentially minimizes the chance of their being reprogrammed into immunosuppressive cells by the TME. This is a risk that other longer-lived immune cells might face. These distinctive characteristics underscore the potential value of neutrophils in anti-tumor immune response. Further studies may explore direct comparisons of antigen presentation capabilities between HLA-DR+ neutrophils and other professional APCs.
Importantly, our results may offer opportunities for superior immunotherapy. First, different from neutrophil depletion therapies, our neutrophil delivery strategy requires only the sorting of autologous circulating neutrophils and short-term ex vivo stimulation, making it a safer alternative that does not leave patients susceptible to infection. Second, HLA-DR+ neutrophils can be loaded with a broad spectrum of (neo)antigens without the need for complex genetic engineering, which is a potential advantage over CAR-T cells with limited target antigens, complicated purification processes, and costly T cell engineering. However, the adverse events of such therapy are still worth studying. Third, given the short half-life of neutrophils, the side effects of this therapy might be transient and manageable. This characteristic might minimize the duration of any potential side effects, allowing for timely intervention and treatment alteration. Comprehensive preclinical and clinical studies will be needed to fully understand the safety profile of this promising therapy.
In summary, our dataset adds to the growing evidence that the cellular neighborhood within tumors is critical for shaping immune responses. Neutrophil transcriptome profiling allows for simultaneously looking at what cell states are, via the gene program they express, and how they reside across diverse cancer types. Our data will help unravel the cellular circuit of neutrophil subsets as well as potentially bridge the gap separating metabolism, epigenetic modification, and innate immune cell phenotypes. Our study may provide opportunities for modulating (neo)antigen-specific immune responses by developing neutrophil therapy that could potentially complement existing cancer immunotherapies.
Limitations of the study
The exploration of state-switching in tumor-associated neutrophils (e.g., generation of HLA-DR+ neutrophils) continues to be valuable in a variety of contexts, including tumors with diverse genetic backgrounds, immunotherapy-treatment refractory tumors, and metastatic tumors. Our studies covered a relatively limited sample size (i.e., multi-cancer TMA cohort containing 8 cancer types), and further validation of neutrophil states across larger cohorts is needed to support the generalizability of these findings. Although our study demonstrated that leucine can enhance neutrophil-dependent antigen presentation and anti-tumor immunity, the potential adverse effects of a leucine-rich diet remain poorly defined. It also remains elusive why leucine intensity varies among different tumors and whether the antigen-presenting function of neutrophils can be lost in leucine-poor tumors. In addition, the potential metabolic variation among neutrophil subsets remains unclear. For example, although HLA-DR+ neutrophils showed increases in leucine metabolism, the TXNIP+ subset displayed unique enrichment in histidine, arginine, and proline metabolism. Testing whether neutrophils can preferentially take up or utilize these amino acids could provide insights into their functional differences.
We observed increased T cell proliferation with HLA-DR+ neutrophils; further work is needed to detail how they directly activate T cells. Testing antigen presentation to T cells in vitro and analyzing CD4+/CD8+ T cell responses in vivo could elucidate this. Tracking adoptively transferred HLA-DR+ neutrophils would also shed light on their migration patterns and ability to reach draining lymph nodes. Although our data demonstrate HLA-DR+ neutrophils can activate T cells and express high CCL5, which can recruit T cells, the relationship between localized neutrophil-T cell interactions versus broader effects on T cell infiltration warrants further exploration. It is also interesting to explore the junctions and synapses between T cells and antigen-presenting neutrophils to provide further insight into their intercellular communication beyond static snapshot analyses. Tracing the half-life of neutrophils (in the tumor, blood, and lymph nodes) will be informative for their clinical applicability. Future research should focus on developing more effective strategies for reprogramming neutrophils into anti-tumor states such as antigen-presenting states.
STAR★Methods
Key resources table
REAGENT or RESOURCE | SOURCE | IDENTIFIER |
---|---|---|
Antibodies | ||
Biotin anti-human CD66b antibody, Clone G10F5 | BioLegend | Cat#305120, RRID: AB_2566608 |
PE anti-human CD66b antibody, Clone G10F5 | BioLegend | Cat#305105, RRID: AB_10550093 |
HLA-DR anti-human antibody, Clone LN3 | Thermo Fisher | Cat#14-9956, RRID: AB_468638 |
CD15 anti-human antibody | Abnova | Cat#MAB-0015 |
Rabbit Cd74 antibody, reacts with: mouse, Clone EPR25399-94 | Abcam | Cat#ab289885 |
Rabbit Ly6G antibody, reacts with: mouse, Clone EPR22909-135 | Abcam | Cat#ab238132 |
Biotin anti-human CD3 antibody, Clone UCHT1 | Biolegend | Cat#300404, RRID: AB_314058 |
PerCP/Cyanine5.5 anti-human CD66b antibody, Clone G10F5 | Biolegend | Cat#305108, RRID: AB_2077855 |
PE anti-human CD66b antibody, Clone G10F5 | Biolegend | Cat#305105, RRID: AB_10550093 |
FITC anti-human CD66b antibody, Clone G10F5 | Biolegend | Cat#305104, RRID: AB_314496 |
Brilliant Violet 421™ anti-human HLA-DR antibody, Clone L243 | Biolegend | Cat#307636, RRID: AB_2561831 |
eFluor660 anti-human Osteopontin (SPP1) antibody, Clone 2F10 | Invitrogen | Cat#50-9096-41 |
Alexa Fluor(R) 488 anti-human CD182 (CXCR2) antibody, Clone 5E8/CXCR2 | Biolegend | Cat#320712, RRID: AB_492938 |
PEcy5 anti-human CD62L antibody, Clone DREG-56 | Invitrogen | Cat#1946541 |
PE anti-human CD54 (ICAM1) antibody, Clone HCD54 | Biolegend | Cat#322708, RRID: AB_535980 |
APC anti-human CD3 antibody, Clone UCHT1 | Biolegend | Cat#300458, RRID: AB_2564151 |
Brilliant Violet 785™ anti-human CD8a antibody, Clone RPA-T8 | Biolegend | Cat#301046, RRID: AB_2563264 |
APC/Fire™ 750 anti-human CD4 antibody, Clone SK3 | Biolegend | Cat#344638, RRID: AB_2572097 |
PE/Dazzle™ 594 anti-human/mouse Granzyme B Recombinant antibody, Clone QA16A02 | Biolegend | Cat#372216, RRID: AB_2728383 |
Brilliant Violet 711™ anti-human IFN-gamma antibody, Clone 4S.B3 | Biolegend | Cat#502539, RRID: AB_11218602 |
BV650 TNF-α antibody, MAb11 | BD | Cat#563418, RRID: AB_2738194 |
FITC anti-human 4-1BB (CD137) antibody, Clone 4B4 | eBioscience | Cat#11-1379-42 |
BV605 anti-human CD69 antibody, Clone FN50 | BD | Cat#562989, RRID: AB_2737935 |
Histone H3 (D1H2) XP®Rabbit mAb(Alexa Fluor®647 Conjugate) antibody, Clone D1H2 | Cell Signaling | Cat#12230S, RRID: AB_2797852 |
Acetyl-Histone H3 (Lys27) (D5E4) XP® Rabbit mAb (Alexa Fluor® 647 Conjugate) antibody, Clone D5E4 | Cell Signaling | Cat#39030S, RRID: AB_2799145 |
FITC-conjugated OVA | Sangon | Cat#D110528 |
PE anti-mouse Ly-6G antibody, Clone 1A8 | Biolegend | Cat#127608, RRID: AB_1186099 |
Alexa Fluor(R) 488 anti-mouse CD74 (CLIP) antibody, Clone In1/CD74 | Biolegend | Cat#151006, RRID: AB_2750326 |
Alexa Fluor(R) 647 anti-mouse CD74 (CLIP) antibody, Clone In1/CD74 | Biolegend | Cat#151004, RRID: AB_2632609 |
FITC Anti-Mouse Cd45 antibody | Tonbo | Cat#35-0451-U025 |
Alexa Fluor® 700 anti-mouse CD3 antibody | BD | Cat#561388, RRID: AB_10642588 |
BV711 anti-mouse Cd8a antibody, Clone 53-6.7 | BD | Cat#563046, RRID: AB_2737972 |
PEcy5 anti-mouse Cd4 antibody, Clone RM4-5 | BD | Cat#553050, RRID: AB_394586 |
APC anti-mouse CD62L (L-Selectin) antibody, Clone MEL-14 | eBioscience | Cat#17-0621-83, RRID: AB_469411 |
PE/Cyanine7 anti-mouse CD19 antibody, Clone 6D5 | Biolegend | Cat#115519, RRID: AB_313654 |
BV395 anti-mouse Cd11b antibody, Clone M1/70 | BD | Cat#563553, RRID: AB_2738276 |
PerCP-Cyanine5.5 anti-mouse Ly6C antibody, Clone HK1.4 | eBioscience | Cat#45-5932-82, RRID: AB_2723343 |
Brilliant Violet 785™ anti-mouse F4/80 antibody, Clone BM8 | Biolegend | Cat#123141, RRID: AB_2563667 |
Brilliant Violet 421™ anti-mouse CD279 (PD-1) antibody, Clone 29F.1A12 | Biolegend | Cat#135221, RRID: AB_2562568 |
PE-Cyanine5.5 anti-mouse CD11c antibody, Clone N418 | eBioscience | Cat#35-0114-82, RRID: AB_469709 |
Emapalumab (anti-IFNγ) antibody | Selleck | Cat#A2041 |
TNFα neutralizing antibody | Sino Biological | Cat#10602-MM0N1 |
IL-6 neutralizing antibody | Sino Biological | Cat#10395-R508 |
IL-17 neutralizing antibody | Sino Biological | Cat#12047-M237 |
IL-23 neutralizing antibody | Sino Biological | Cat#CT035-mh066 |
Ultra-LEAF™ Purified anti-mouse CD279 (PD-1), Clone 29F.1A12 | BioLegend | Cat#135248 |
Ultra-LEAF™ Purified Rat IgG2a, κ Isotype Ctrl | BioLegend | Cat#400565 |
InVivoPlus anti-mouse Ly6G/Ly6C (Gr-1) antibody, clone RB6-8C5 | Bio X Cell | Cat#BE0075, RRID: AB_10312146 |
Chemicals, peptides, and recombinant proteins | ||
Alanine | Sangon | Cat#A600022-0100 |
Arginine | Sangon | Cat#A600205-0100 |
Asparagine | Sangon | Cat#A694341-0100 |
Aspartate | Sangon | Cat#A600091-0250 |
Cysteine | Sangon | Cat#A600132-0100 |
Glutamine | Sangon | Cat#A100374-0050 |
Glutamate | Sangon | Cat#A600221-0500 |
Glycine | Sangon | Cat#A610235-0500 |
Histidine | Sangon | Cat#A604351-0050 |
Isoleucine | Sangon | Cat#A100803-0050 |
Leucine | Sangon | Cat#A600922-0100 |
Lysine | Sangon | Cat#A602759-0025 |
Methionine | Sangon | Cat#A610346-0100 |
Phenylalanine | Sangon | Cat#A600991-0025 |
Proline | Sangon | Cat#A600923-0100 |
Serine | Sangon | Cat#A601479-0100 |
Threonine | Sangon | Cat#A610919-0100 |
Tryptophan | Sangon | Cat#A601911-0050 |
Tyrosine | Sangon | Cat#A601932-0100 |
Valine | Sangon | Cat#A600172-0025 |
L-Leucine-13C6 | Sigma | Cat#605239 |
Dichloroacetate | Sigma | Cat#2156-56-1 |
ACSS2-IN-2 | MCE | Cat#2332820-04-7 |
TMRE Fluorescent Mitochondrial Probe | Sigma | Cat#87917-25MG |
NAO nonyl bromide | Sigma | Cat#A7847-100MG |
Fluo 3 | Sigma | Cat#73881-1MG |
JC1 | AAT Bioquest | Cat#22200 |
Brite™ HPF ∗Optimized for Detecting Reactive Oxygen Species (ROS) | AAT Bioquest | Cat#16051 |
Trizol | Thermo Fisher Scientific | Cat#15596018 |
NEBNext UltraTM RNA Library Prep Kit | NEB | Cat#E7490 |
Collagenase IV | STEMCELL technologies | Cat#07909_C |
RPMI1640 | Gibco | Cat#11875500BT |
Peptide VVRHCPHHERCSDSD | China Peptides Inc. | N/A |
Peptide QHMTEVVRHCPHHER | China Peptides Inc. | N/A |
Peptide RNTFRHSVVVPCE | China Peptides Inc. | N/A |
Peptide NTFRHSVVVPCEPPE | China Peptides Inc. | N/A |
Peptide HYNYMCNSSCMGSMN | China Peptides Inc. | N/A |
Peptide MTEYKLVVVGAVGVGKSALTIQLI | China Peptides Inc. | N/A |
Peptide LVVVGADGV | China Peptides Inc. | N/A |
Peptide SQEQPRCHY | China Peptides Inc. | N/A |
Peptide RLFERDGLKV | China Peptides Inc. | N/A |
Peptide LVVVGADGV | China Peptides Inc. | N/A |
Peptide gp100 (44-59) | China Peptides Inc. | N/A |
Peptide CMV (6571) | China Peptides Inc. | N/A |
MHC class II antigen presentation Gene Expression PCR Array | Wcgene biotech | Cat#WC-MRNA0283-H |
MHC class I antigen presentation Gene Expression PCR Array | Wcgene biotech | Cat#WC-MRNA0282-H |
DAPI | BioLegend | Cat#422801 |
Critical commercial assays | ||
Chromium™ Single Cell 5′ Library Construction Kit | 10x Genomics | Cat#1000020 |
Chromium™ Next GEM Single Cell 5′ Library and Gel Bead Kit v1.1 | 10x Genomics | Cat#1000165 |
MojoSort™ Whole Blood Human Neutrophil Isolation Kit | BioLegend | Cat#480152 |
Anti-Biotin MicroBeads | Miltenyi Biotec | Cat#130-090-485 |
EasySep™ Direct Human PBMC Isolation Kit | StemCell Technologies | Cata#19654 |
Chromium Next GEM Single Cell 3′ Kit v3.1 | 10x Genomics | Cat#1000268 |
KC-digital™ stranded TCR-seq library prep kit | Seqhealth Technology Co., Ltd | Cat#DT0813-02 |
Experimental models: Cell line | ||
Human: HepG2 cells | Cell Bank of Type Culture Collection Chinese Academy of Sciences (CBTCCCAS) | SCSP-510 |
Human: A549 cells | CBTCCCAS | SCSP-503 |
Human: HCT116 cells | CBTCCCAS | SCSP-5076 |
Human: PANC1 cells | CBTCCCAS | SCSP-535 |
Human: MCF7 cells | CBTCCCAS | SCSP-531 |
Human: dHL-60 cells | Genomeditech Co. Ltd. | N/A |
Mouse: MC38 cells | Shanghai Model Organisms Center | N/A |
Mouse: Hepa 1–6 cells | Shanghai Model Organisms Center | N/A |
Mouse: LLC cells | Shanghai Model Organisms Center | N/A |
Experimental models: Organisms/strains | ||
Mouse: C57BL/6J wildtype | Shanghai Model Organisms Center | SM-001 |
Mouse: C57BL/6J Cd45.1 | Shanghai Model Organisms Center | NM-KI-210226 |
Mouse: C57BL/6J Cd74-KO | Shanghai Model Organisms Center | NM-KO-200715 |
Mouse: C57BL/6J MHC-IIflox/flox | Nanjing GemPharmatech Co. Ltd. | T019085 |
Mouse: C57BL/6J Ly6GCre-tdTomato | Shanghai Model Organisms Center | NM-KI-200219 |
Mouse: C57BL/6J Lat1(Slc7a5)-KO | Nanjing GemPharmatech Co. Ltd. | T031657 |
Mouse: C57BL/6J Bcat2-KO | Nanjing GemPharmatech Co. Ltd. | T049844 |
Mouse: C57BL/6J Dbt-KO | Nanjing GemPharmatech Co. Ltd. | T031427 |
Biological samples | ||
Neutrophil in-house single-cell RNA-seq data (n = 155; 103 samples were newly generated data) | Zhongshan Hospital, Fudan University Summarized in Table S1 | N/A |
Neutrophil single-cell RNA-seq data derived from public data (n = 70, cancer patients; n = 5, healthy donor) | Summarized in Table S1 | N/A |
Spatial transcriptomics data derived from public data (n = 50) | Summarized in Table S5 | N/A |
Tissue microarray (8-Cancer-TMA) cohort with survival information (n = 1,116) | Zhongshan Hospital, Fudan University; Shanghai Outdo Biotech Co. Ltd. Summarized in Table S1 | N/A |
Tissue microarray (Multi-Cancer-TMA) cohort (n = 68) | Shanghai Outdo Biotech Co. Ltd. Summarized in Table S1 | N/A |
Blood from healthy donor (n = 44) | Zhongshan Hospital, Fudan University | N/A |
Hepatocellular carcinoma samples treated with neoadjuvant immunotherapy (n = 5) | Zhongshan Hospital, Fudan University | N/A |
Deposited data | ||
scRNA-Seq of neutrophils | This paper | PRJCA020880; http://pancancer.cn/neu |
scRNA-Seq of neutrophils | Qian, J. et al.73 | E-MTAB-8107, E-MTAB-6149 and E-MTAB-6653 |
scRNA-Seq of neutrophils | Chan, J. et al.74 | Human Tumor Atlas Network (HTAN) |
scRNA-Seq of neutrophils | Yang, L. et al.75 | GSE171145 |
scRNA-Seq of neutrophils | Zilionis et al.23 | GSE127465 |
scRNA-Seq of neutrophils | Wang et al.25 | OEP003254 |
scRNA-Seq of neutrophils | Hu, S. et al.76 | HRA001006 |
scRNA-Seq of neutrophils | Xue et al.7 | PRJCA007744 |
scRNA-Seq of neutrophils | Tabula Sapiens Consortium et al.77 | GSE201333 |
Spatial transcriptomics | Summarized in Table S5 | 10X Genomics website; http://lifeome.net/supp/livercancer-st/data.htm; https://zenodo.org/record/4739739; GSE144239; GSE175540 |
RNA-seq for inferring the neutrophil consensus infiltration | The Cancer Genome Atlas | N/A |
RNA-seq for inferring the neutrophil consensus infiltration | Clinical Proteomic Tumor Analysis Consortium | N/A |
RNA-seq of neutrophils (leucine treatment and control) | This paper | PRJCA020880 |
ATAC-seq of neutrophils (leucine treatment and control) | This paper | PRJCA020880 |
CUT&Tag of neutrophils (leucine treatment and control) | This paper | PRJCA020880 |
TCR-seq (antigen-presenting neutrophil stimulating T-cell response) | This paper | PRJCA020880 |
scRNA-Seq of mouse tumors | Summarized in Table S5 | N/A |
RNA-seq of immunotherapy-treated samples (SKCM) | Gide et al.54 | PRJEB23709 |
RNA-seq of immunotherapy-treated samples (SKCM) | Hugo et al.55 | GSE78220 |
RNA-seq of immunotherapy-treated samples (BLCA) | Mariathasan et al.56 | EGAS00001002556 |
RNA-seq of immunotherapy-treated samples (NSCLC) | Prat et al.57 | GSE93157 |
RNA-seq of immunotherapy-treated samples (SKCM) | Nathanson et al.58 | N/A |
RNA-seq of immunotherapy-treated samples (SKCM) | Lauss et al.59 | GSE100797 |
RNA-seq of immunotherapy-treated samples (STAD) | Kim et al.60 | PRJEB25780 |
RNA-seq of immunotherapy-treated samples (HCC) | Zhu et al.61 | EGAS00001005503 |
Software and algorithms | ||
CellRanger V7 | 10x Genomics | https://10xgenomics.com |
Seurat V4.0.4 | CRAN | https://cran.r-project.org/web/packages/Seurat/index.html |
harmony V0.1.0 | CRAN | https://cran.r-project.org/web/packages/harmony/index.html |
ggplot2 V3.3.5 | CRAN | https://cran.r-project.org/web/packages/ggplot2/index.html |
dittoSeq V1.5.2 | Bioconductor | https://bioconductor.org/packages/dittoSeq/ |
GSVA V1.40.1 | Bioconductor | https://www.bioconductor.org/packages/GSVA/ |
Monocle3 V1.0.0 | Github | https://github.com/cole-trapnell-lab/monocle3 |
shiny V1.6.0 | CRAN | https://cran.r-project.org/package=shiny |
TooManyCells V2.0.0.0 | github | https://github.com/GregorySchwartz/too-many-cells |
SingleR V1.7.1 | Bioconductor | https://bioconductor.org/packages/SingleR |
ggpubr V0.4.0 | CRAN | https://cran.r-project.org/package=ggpubr |
ggsignif V0.6.3 | CRAN | https://cran.r-project.org/web/packages/ggsignif/index.html |
pheatmap V1.0.12 | CRAN | https://cran.r-project.org/web/packages/pheatmap/index.html |
ComplexHeatmap V2.15.1 | Bioconductor | https://bioconductor.org/packages/release/bioc/html/ComplexHeatmap.html |
cowplot V1.1.1 | CRAN | https://cran.r-project.org/web/packages/cowplot/index.html |
sctour V0.1.3 | Pypi | https://pypi.org/project/sctour/ |
xCell V1.1.0 | Github | https://github.com/dviraran/xCell |
dorothea V1.10.0 | Bioconductor | https://bioconductor.org/packages/release/data/experiment/html/dorothea.html |
UCell V1.3.1 | Bioconductor | https://bioconductor.org/packages/release/bioc/html/UCell.html |
ScMetabolism | Github | https://github.com/wu-yc/scMetabolism |
doubletFinder V2.0.3 | Github | https://github.com/chris-mcginnis-ucsf/DoubletFinder |
Other | ||
Code | This paper | https://github.com/wu-yc/neutrophil (https://doi.org/10.5281/zenodo.10531210) |
ScProgram | This paper | https://github.com/wu-yc/scProgram (https://doi.org/10.5281/zenodo.10531218) |
Resource availability
Lead contact
Further information and requests for resources and reagents should be directed to and will be fulfilled by the lead contact, Qiang Gao (gaoqiang@fudan.edu.cn).
Materials availability
This study did not generate new unique reagents.
Data and code availability
- •
Processed gene expression data can be queried and downloaded at http://www.pancancer.cn/neu and raw sequencing data are available at the China National Center for Bioinformation (accession: PRJCA020880) with the permission at Human Genetic Resources Service System of Ministry of Science and Technology. To request access to raw sequencing data, please apply at Human Genetic Resources Service System of Ministry of Science and Technology (https://apply.hgrg.net/) according to the law of Regulations on management of human genetic resources of China. This paper analyzed existing, publicly available data, where the accession numbers are listed in the key resources table.
- •
All original code has been deposited at GitHub and Zenodo and is available at https://github.com/wu-yc/neutrophil (https://doi.org/10.5281/zenodo.10531210) and https://github.com/wu-yc/scProgram (https://doi.org/10.5281/zenodo.10531218).
- •
Any additional information required to reanalyze the data reported in this paper is available from the lead contact upon request.
Experimental model and study participant details
Patient inclusion and sample collection
We collected fresh samples from Zhongshan Hospital Fudan University with written consent and approval from the Institutional Review Board-approved protocols (B2021-381, B2021-084, B2020-348R, B2023-350). After quality control, a total of 103 samples from 64 patients were included in the analyses. The median age of the patients was 59.6 years, and the cohort consisted of 26 females and 38 males, and the detailed information was included in Table S1. Patient inclusion criteria were as follows: Patients with treatment-naïve primary tumors who underwent surgery; patients without major underlying diseases that may seriously affect neutrophils (such as autoimmune diseases and acute infection). The 8-Cancer-TMA (n = 1,116; diameter: 1.5 mm; sample information see Table S1) and Multi-Cancer-TMA (n = 68; diameter: 1.5 mm; tumor and peritumor samples) were from the patients from Zhongshan Hospital Fudan University and Shanghai Outdo Biotech (Table S1) and approval from Shanghai Outdo Biotech Ethics Committee (SHYJS-CP-2210040, SHYJS-CP-1910002, SHYJS-CP-1804011, SHYJS-CP-1510001, SHYJS-CP-1701016, SHYJS-CP-1404018). As for PDTF analysis, samples from 5 HCC patients with pathologic nonresponse (residual viable tumor rates of 95%, 95%, 90%, 90%, and 90% respectively; treatment: anti-PD-1 plus anti-VEGF agent) were included.
Mouse models
5-week-old male C57BL/6 mice, CD74 KO, CD45.1 mice, and Ly6GCre-tdTomato were obtained from the Shanghai Model Organisms Center, Inc. MHC-IIflox/flox mice, LAT1KO, Bcat2KO, DbtKO mice were obtained from Nanjing GemPharmatech Co. Ltd. We housed them under pathogen-free conditions with a maximum of five mice per cage. We strictly adhered to animal care principles and ethics and received approval from the Institutional Animal Care and Use Committee of the Shanghai Model Organisms Center (approval number 2019-0011).
Method details
Standardized cell sorting protocol
Tumor and matched normal samples were dissected into small pieces with a diameter < 1mm, and dissociated using Collagenase IV (STEMCELL technologies; 07909_C) plus 0.4 mg/mL hyaluronidase in RPMI 1640 (Gibco; 11875093) with a GentleMACS Dissociator for 60 minutes. The resulting cells were filtered through a 400 μm filter and washed with DPBS (500g and 10 minutes). To isolate CD66b+ neutrophils, we employed a two-step sorting strategy. First, cells were stained with CD66b Biotin antibody (Biolegend, 305120), and then sorted with MS columns (Miltenyi Biotec, 130-042-201). Next, cells were stained with CD66b PE antibody (Biolegend, 305105) and sorted again by flow cytometry (BD FACS Aria II). Sorted cells were immediately sent for single-cell RNA sequencing. For samples with a sufficient number of cells, we also sorted CD45+ cells and sequenced them. As for the neutrophils derived from blood, we incubated the peripheral blood with red blood cell lysis buffer (Sangon, B541001) for 10 minutes and washed it with DPBS (500g and 10 minutes). The resulting cell suspension was then stained with CD66b PE antibody (Biolegend, 305105) and sorted using flow cytometry (BD FACS Aria II).
Single-cell sequencing
The sorted cells were sequenced using the 10x Chromium single-cell platform with 5′ Reagent Kits following the manufacturer's protocol. Single-cell libraries were then sequenced on the NovaSeq platform from Illumina. To trace the sample source, we used TotalSeq C from BioLegend (399905), which allowed us to distinguish between tumor, adjacent normal, and blood cells. CellRanger V7 was used for processing the barcodes, aligning the data, and generating initial clusters of the raw scRNA-seq profiles.
ScRNA-seq data quality control, processing, annotation, and visualization
Raw fastq files were firstly aligned to human genome (GRCh38, ENSEMBL) by CellRanger V7 following the 10X Genomics neutrophil tutorial (https://support.10xgenomics.com/single-cell-gene-expression/software/pipelines/latest/tutorials/neutrophils), we retained the intronic regions by using the parameter of “--include-introns” and set the “--force-cells=20000” in CellRanger. We then used Seurat (V4.0.4)78 to process the UMI count matrix. We performed doublet removal by using doubletFinder (V2.0.3). The mitochondrial gene percentage was assessed by PercentageFeatureSet(object, pattern = "ˆMT-") function and cells with mitochondrial gene percentage over 10% were removed. Excluding genes in a blacklist as described before,65 the top 5000 genes were identified as highly variable genes (HVG) using FindVariableFeatures function of Seurat.78 We integrated all cells according to sample ID by using harmony (V0.1.0).79 We performed clustering analysis and only reserved the neutrophils defined by markers CSF3R, FPR1, FCGR3B, NAMPT, and MNDA following the 10X Genomics neutrophil tutorial. We further used SingleR (V1.7.1) to confirm the input cells are real neutrophils. In the second round of doublet removal, we observed clusters with high expression of CD3D (T cell), CD79A (B cell), and CD68 (macrophage). We next computed the marker genes using FindAllMarker function of Seurat.78
Single-cell transcriptional program feature gene analysis and signature quantification
We followed a well-established computational strategy designed for decoding single-cell heterogeneity.80 In detail, we first normalized the data matrix using scTransform81 and one-by-one performed NMF analysis (parameter: k = 100). We clustered the Jaccard index of H matrix and visualized it using pheatmap (V1.0.12). We performed the differential gene analysis to find the feature genes of each program. Further gene set enrichment analysis was performed using the highly expressed genes. As for the signature quantification, we used UCell82 to quantify the signature/pathway activity of neutrophils. KEGG,83 hallmark,84 Gene Ontology (GO),85 and neutrophil immunophenotype signatures29 were included for analysis. As for the circadian analysis of neutrophil signature, the sampling time refers to the surgery time of tumor samples.
Single-cell HLA-type quantification
We employed scHLAcount (V0.2.0, available at https://github.com/10XGenomics/scHLAcount) with default parameters to count the molecules of class I genes HLA-A, B, and C; and class II genes DPA1, DPB1, DRA1, DRB1, DQA1, and DQB1.
Single-cell metabolism quantification
To explore the featured pathway among HLA-DR+ neutrophils, we first computed the pathway variance (Gini index, metric entropy, Shannon entropy, and Simpson index). The formula for the Gini index is: G = (2A)/(nB), where G is the Gini index; A is the area between the Lorenz curve and the line of perfect equality of pathway score (the diagonal line from the bottom left to the top right corner of the graph); B is the total area under the line of perfect equality of pathway score; n is the total number of cells. The formula for the Simpson index is: D = 1 - Σ(ni(ni-1))/N(N-1), where: D is the Simpson index; ni is the number of cells in a given pathway. The formula for Shannon entropy is: H = -ΣP(x) log2 P(x), where H is the entropy in bits, P(x) is the probability distribution of pathway scores. The major metabolism subtype quantification was performed using scMetabolism as we developed before37 (parameter: imputation = T, metabolism.type = "KEGG").
Cell-cell communication, cell trajectory, and cell clade analysis
To understand how HLA-DR+ neutrophils interact with T cells, we first performed the down-sampling analysis (to 20,000 cells), impute expression matrix by using ALRA,86 and used CellPhoneDB87 to infer the interactions (parameter: cellphonedb method statistical_analysis --iterations=100 --threads=48). We next ranked the ligand expressed on neutrophils by the gene expression proportion and frequency. We used scTour to infer the differentiation state of neutrophils.32 To validate the trajectory, we next used monocle3,33 CytoTRACE,34 and Slingshot35 to separately infer the pseudotime of neutrophil subsets. To compare the hierarchy difference between neutrophils derived from blood, cancer, and adjacent tissues, we split the expression matrix containing 3000 variable genes and performed clustering analysis using TooManyCells26 (parameter: make-tree PieRing).
Generating the pan-cancer neutrophil infiltration consensus
We download the gene expression matrix of pan-cancer solid tumors by using UCSC Xena88 (data type: HTSeq - FPKM-UQ). We first tested 8 common immune quantification algorithms covering CIBERSORT, ESTIMATE, Quantiseq, MCPCounter, IPS, TIMER, EPIC, and xCell.16,17,18,89,90,91,92,93 We found that 3 of them support the quantification of neutrophil level (MCPCounter, Quantiseq, and xCell). We clustered the pan-cancer samples according to neutrophil level generated by three algorithms and designed the consensus neutrophil score based on the consensus rank of these scores. In detail, as for the samples ranked at the upper quantile or low quantile in three algorithms, we then label the samples as high or low neutrophil consensus samples. While the samples without reaching the consensus among three algorithms are labeled as heterogeneous samples. We ordered the samples according to the consensus score and rank them according to cancer types. We also performed the dimensional reduction analysis using t-SNE embedded in Seurat78 and labeled samples according to their neutrophil consensus status.
Transcription factor activity analysis
For scRNA-seq data, we used dorothea (V1.10.0) to infer the transcription factor activity. As for the ChIP-seq of RFX5 transcription factor, we fetched published RFX5 ChIP-seq data (accession number: SRX150635, SRX150644, SRX150384, SRX186620, SRX186634, SRX150462, and SRX150639) and analyzed it using UCSC Genome Browser.94 The neutrophil-like differentiated HL-60 (dHL-60) cells were obtained by adding 1% DMSO to the HL-60 culture medium for six days.25 We obtained the knockdown and overexpression plasmids of RFX5 from Genomeditech (Shanghai, China). Empty vector was used as the negative control. Each condition was performed with 3 replicates.
Multiplex immunohistochemistry analysis
We performed the immunohistochemistry using Osteopontin/SPP1 (Abcam; ab214050; species reactivity: Human), HLA-DR (thermofisher; 14-9956; species reactivity: Human), CD15 (MAB-0015; species reactivity: Human), Cd74 (Abcam; ab289885; species reactivity: Mouse), Ly6G (Abcam; ab238132; species reactivity: Mouse), CXCL13 (Abcam; ab246518; species reactivity: Human), CD39 (Abcam; ab300065; species reactivity: Human), CD4 (Biolynx; BX50023; species reactivity: Human), CD8 (Dako; M7103; species reactivity: Human), MPO (Abcam; ab300650; species reactivity: Human), CD11b (Abcam; ab133357; species reactivity: Human), DAPI (BioLegend; 422801) antibodies. We scanned the slides using the PerkinElmer Vectra3 platform and quantified the results by using PerkinElmer Vectra3 platform as previously described.52,95,96,97 The 8-Cancer-TMA cohort with matched prognosis metadata (HCC, COAD, NSCLC, STAD, RCC, OV, BRCA, and BLCA) and Multi-Cancer-TMA cohort (PAAD, RCC, HCC, ICC, STAD, NSCLC, BRCA, and COAD) were included. The detailed clinicopathological features were described in Table S1.
Cell culture
We performed Ficoll Paque experiment on blood derived from healthy donors using Ficoll Paque Plus agent (GE, 17-1440-03). Neutrophils were separated by staining with CD66b Biotin antibody (Biolegend, 305120), adding biotin magnetic beads, and sorting with MS columns (Miltenyi Biotec, 130-042-201). T cells were separated by staining with CD3 Biotin antibody (Biolegend, 300404), adding biotin magnetic beads, and sorting with MS columns (Miltenyi Biotec, 130-042-201). Cells (5 ×104) were added to 96-well cell culture plates in a total volume of 200 μL of culture medium. To maintain neutrophil activity, we added Lipopolysaccharides (LPS, MCE, HY-D1056) at a concentration of 100 ng/mL. Additionally, we supplemented the medium with 20 types of amino acids and performed the neutrophil culture for 24 hours (Sangon, A600022-0100, A600205-0100, A694341-0100, A600091-0250, A600132-0100, A100374-0050, A600221-0500, A610235-0500, A604351-0050, A100803-0050, A600922-0100, A602759-0025, A610346-0100, A600991-0025, A600923-0100, A601479-0100, A610919-0100, A601911-0050, A601932-0100, A600172-0025). The concentration of each amino acid was set to match the physiological plasma concentration, as reported by healthmatters.io: Alanine (681 μmol/L), Arginine (137 μmol/L), Asparagine (90 μmol/L), Aspartate (12.6 μmol/L), Cysteine (360 μmol/L), Glutamine (876 μmol/L), Glutamate (214 μmol/L), Glycine (518 μmol/L), Histidine (114 μmol/L), Isoleucine (104 μmol/L), Leucine (196 μmol/L), Lysine (318 μmol/L), Methionine (48 μmol/L), Phenylalanine (95 μmol/L), Proline (363 μmol/L), Serine (172 μmol/L), Threonine (216 μmol/L), Tryptophan (83 μmol/L), Tyrosine (110 μmol/L), Valine (370 μmol/L). Dichloroacetate (mitochondrial acetyl-CoA activator, 0.1 mM; 2156-56-1), ACSS2-IN-2 (acetyl-CoA inhibitor, 5 nM; 2332820-04-7) were used respectively. Each condition was performed with 3-4 replicates. As for the HLA-DR+ neutrophil cocultured with T cells, neutrophils were sorted from samples of hepatocellular carcinoma, colorectal cancer, and non-small cell lung cancer; while T cells were sorted from the PBMC from the same patient.
Flow cytometry
To perform surface staining, we mixed the appropriate antibodies with the cells at room temperature for 15 minutes and washed them with DPBS (500g, 10 minutes). For intracellular staining, we used the Fixation Permeabilization Kit (BD, 554714) to fix and permeabilize the cells, followed by staining with the appropriate antibodies in the Permeabilization buffer for 30 minutes at 4°C, as previously described.98 The following antibodies were used: CD66b PerCP-Cy5.5 (Biolegend 305108), CD66b PE (Biolegend 305105), CD66b FITC (Biolegend 305104), HLA-DR BV421 (Biolegend 307636), Osteopontin (SPP1) eFluor660 (Invitrogen 50-9096-41), CXCR2 A488 (Biolegend 320712), CD62L PEcy5 (invitrogen 1946541), CD54 (ICAM1) PE (Biolegend 322708), CD3 APC (Biolegend 300458), CD8a BV785 (Biolegend 301046), CD4 APC/Fire750 (Biolegend 344638), Granzyme B PE/Dazzle594 (Biolegend 372216), IFN-γ BV711 (Biolegend 502539), TNF-α BV650 (BD 563418), 4-1BB (CD137) FITC (eBioscience 11-1379-42), CD69 BV605 (BD 562989), Histone H3 A647 (Cell Signaling 12230S), Acetyl-Histone H3 Lys27 A647 (Cell Signaling 39030S), OVA FITC (sangon D110528), Ly6G PE (Biolegend 127608), Cd74 A488 (Biolegend 151006), Cd74 A647 (Biolegend 151004), Cd45 FITC (TONBO 35-0451-U025), Cd3 A700 (BD 561388), Cd8a BV711 (BD 563046), Cd4 PEcy5 (BD 553050), Cd62L APC (eBioscience 17-0621-83), Cd19 PEcy7 (Biolegend 115519), Cd11b BV395 (BD 563553), Ly6C PerCPCy5.5 (eBioscience 45-5932-82), F4/80 BV785 (Biolegend 123141), Cd279 BV421 (Biolegend 135221), CD11c PEcy5.5 (eBioscience 35-0114-82), CD16 BV711 (Biolegend 302044). Cells were further analyzed using flow cytometry (BD LSRFortessa) and Flowjo software (BD).
PCR array analysis of MHC class I and II antigen presentation genes
We stimulated neutrophils from healthy donors’ blood with leucine for 24 hours. Gene expression profiling was carried out using the human MHC class I antigen presentation Gene Expression PCR Array (Wcgene Biotech, Shanghai, China) and human MHC class II antigen presentation Gene Expression PCR Array (Wcgene Biotech, Shanghai, China) following the manufacturer's protocol.
Metabolomics LC-MS analysis and leucine tracing
The Thermo Vanquish ultra-high performance liquid phase system (Thermo Fisher Scientific, USA) equipped with an ACQUITY UPLC® HSS T3 column (2.1×150 mm, 1.8 μm) (Waters, Milford, MA, USA) was utilized. The system operated with a flow rate of 0.25 mL/min, a column temperature of 40 °C, and an injection volume of 2 μL. Mass spectral data was collected using the Thermo Orbitrap Exploris 120 mass spectrometer detector (Thermo Fisher Scientific) with electrospray ionization source (ESI). Both positive and negative ion modes were used for data collection, with a positive ion spray voltage of 3.50 kV and a negative ion spray voltage of -2.50 kV. The sheath gas and auxiliary gas were set to 30 arb and 10 arb, respectively. The primary full scan was performed at a resolution of 60,000 over the m/z range of 100-1000, and HCD was utilized for secondary fragmentation with a collision voltage of 30%. The secondary resolution was set to 15,000. The MS data analysis were conducted following a previously established protocol.99 Raw MS data were converted to the mzXML format using ProteoWizard software (http://proteowizard.sourceforge.net). Peaks were extracted using R package XCMS (V3.20.0). The peak table and MS2 files in mgf format (converted using ProteoWizard) were uploaded to the MetDNA web server (http://metdna.zhulab.cn/) for metabolite identification. The identifications were assigned levels 1-3 and unknown, following the MSI (Metabolomics Standard Initiative) guidelines. For the tissue samples, 10 scRNA-seq-matched samples (HCC, n=4; NSCLC, n=2; OV, n=3; STAD, n=1) passed QC and was prepared for LC-MS. As for the neutrophils, autologous neutrophils from blood were separated by staining with CD66b Biotin antibody (Biolegend, 305120), added with biotin magnetic beads, sorted with MS columns (Miltenyi Biotec, 130-042-201), and underwent stimulation with leucine for 24 hours.
For 13C leucine tracing, we treated neutrophils with L-Leucine-13C6 (Sigma 605239) for 24 hours following the published protocol.44 The untreated neutrophils were used as controls. Each condition was performed with 3-4 replicates.
Mitochondrial functional and phenotypic characterization
TMRE (Sigma 87917-25MG), NAO nonyl bromide (Sigma A7847-100MG), Fluo 3 (Sigma 73881-1MG), JC1 (AAT Bioquest 22200), and ROS (AAT Bioquest 16051) were suspended with neutrophils at room temperature for 30min and washed using DPBS (500g, 10 minutes). Cells were further analyzed using flow cytometry (BD LSRFortessa) and Flowjo software (BD). Leica TCS SP5 laser confocal microscope was also used to image neutrophil mitochondria. Each condition was performed with 3 replicates.
Seahorse assays
Autologous neutrophils from blood were separated by staining with CD66b Biotin antibody (Biolegend, 305120), added with biotin magnetic beads, and sorted with MS columns (Miltenyi Biotec, 130-042-201). Post-separation, the neutrophils underwent stimulation with leucine for 24 hours. Cultured plates were used to plate 5 × 105 neutrophils sorted from healthy donors' blood. The neutrophils were then stimulated with leucine for 24 hours, and OCR was measured using an XF24 Seahorse Extracellular Flux Analyzer following the manufacturer's instructions. In the seahorse assays, the neutrophils were treated with oligomycin (0.25 μM), FCCP (0.25 μM), rotenone (0.25 μM), and antimycin A (0.25 μM). Each condition was performed with 3 replicates.
Transmission electron microscopy (TEM)
The neutrophils were fixed in 1% osmium tetroxide in PBS in the dark at room temperature for 2 hours, washed with PBS (pH 7.4) three times for 15 minutes each, and then dehydrated in a series of alcohol concentrations (30%-50%-70%-80%-90%-95%-100%-100%) for 15 minutes each. The cells were then embedded in epoxy resin. The resin blocks were sectioned into ultrathin sections (60-80 nm) using an ultramicrotome, and collected on 150-mesh copper grids. The grids were stained with 2% uranyl acetate in saturated alcohol solution in the dark for 8 minutes, washed three times with 70% alcohol and three times with ultra-pure water, and then stained with 2.6% lead citrate for 8 minutes. For each treatment group, we performed the TEM for 5 cells. Images were obtained using a SUHT7700 electron microscope (hitachi) and subsequently analyzed using ImageJ Fiji (V2.11.0). In detail, we select the "Straight line" tool and measure the mitochondrial length. The exported length was then analyzed using R.
Confocal imaging
To stain the neutrophils with JC1, we incubated them in a staining solution for 30 minutes at 37°C. After removing the staining solution, we washed the cells with assay buffer to remove any unbound dye. We utilized an Olympus SpinSR10 Ixplore confocal microscope to capture images of the cells. Subsequently, we employed ImageJ Fiji (V2.11.0) software to analyze the images.
Bulk RNA-seq
As for the scRNA-seq-matched samples, sufficient remaining fresh specimens were quickly frozen in liquid nitrogen. As for the neutrophil samples, autologous neutrophils from blood were separated by staining with CD66b Biotin antibody (Biolegend, 305120), added with biotin magnetic beads, and sorted with MS columns (Miltenyi Biotec, 130-042-201). Post-separation, the neutrophils underwent stimulation with leucine for 24 hours. Samples were subjected to RNA extraction using Trizol (Thermo Fisher Scientific, 15596018). Library preparation was performed using NEBNext UltraTM RNA Library Prep Kit (NEB #E7490), followed by library purification using beads (AMPure XP system). Finally, sequencing was conducted using NovaSeq 6000 with PE150. For alignment, the STAR software was employed (https://github.com/alexdobin/STAR). The raw expression levels of each gene (based on fragment counts) were calculated using htseq-count (https://htseq.readthedocs.io/en/release_0.11.1). The selection criteria for significantly differentially expressed genes were: |log2FC| > 1 and P-value < 0.05. Neutrophil RNA-seq was performed with 4 replicates. For the immune deconvolution analysis, we use the xCell method.18 First, xCell includes the most types of T cells among the methods we tested, providing us with a more comprehensive view of the T cell landscape in our samples. Secondly, xCell is one of the most widely used immune cell quantification methods in the field.18
H3K27ac, H3K27me3, and H3K4me3 CUT&Tag
Autologous neutrophils from blood were separated by staining with CD66b Biotin antibody (Biolegend, 305120), added with biotin magnetic beads, and sorted with MS columns (Miltenyi Biotec, 130-042-201). Post-separation, the neutrophils underwent stimulation with leucine for 24 hours. After this, procedures of cell permeabilization, antibody incubation, tagmentation, DNA extraction, and sequencing were carried out (Shanghai Jiayin Biotech). For data filtering, the raw reads were processed using Trimmomatic (V0.35, http://www.usadellab.org/cms/?page=trimmomatic). BWA software https://bio-bwa.sourceforge.net/ was used for alignment. Fragment sizes for read pairs were calculated using the BAM file from aligned paired-end sequencing data. The summary statistics on fragment lengths were estimated by sampling several regions, depending on the size of the genome and number of processors. MACS2 (V2.2.7.1, https://pypi.org/project/MACS2/) was used for peak calling in this analysis, Bedtools (V2.30.0, https://bedtools.readthedocs.io/en/latest/) were mainly used for peak annotation analysis.
ATAC-seq
Autologous neutrophils from blood were separated by staining with CD66b Biotin antibody (Biolegend, 305120), added with biotin magnetic beads, and sorted with MS columnsMiltenyi Biotec130-042-201Post-separation, the neutrophils underwent stimulation with leucine for 24 hours. Cell nuclei were subsequently extracted and underwent a transposition reaction via the Tn5 enzyme (Illumina) on 40,000 cell nuclei. Sequencing was performed on an Illumina Novaseq6000using a PE150 sequencing strategy (Shanghai Jiayin Biotech)The raw reads were filtered using Trimmomatic(V0.35, http://www.usadellab.org/cms/?page=trimmomatic)Data alignment was conducted using the BWA softwarehttps://bio-bwa.sourceforge.net/Then the aligned reads in BAM file were used to calculate the fragment sizes for read pairsThe estimation of summary statistics on fragment lengths was done by sampling various regions, with the selection depending on the genome size and number of processorsMACS2(V2.2.7.1, https://pypi.org/project/MACS2/) was used for peak calling in this analysis. Lastly, peak annotation analysis primarily employed Bedtools (V2.30.0, https://bedtools.readthedocs.io/en/latest/).
Spatial transcriptomics
We processed the space ranger output files using Seurat (V4.0.4). The sample information was summarized in Table S5. Afterward, we utilized Seurat's SCTransform function for data normalization, RunPCA function for dimension reduction, and FindNeighbors and FindClusters function for ST spot clustering. To score the cell types, we employed xCell (V1.1.0) and estimated the neutrophil signature using GSVA (V1.40.1).
Neoantigen, T cell, and neutrophil coculture system
Neutrophils were autonomously harvested from the blood of healthy donor through a process involving CD66b Biotin antibody (Biolegend, 305120) staining, the addition of biotin magnetic beads, and sorting using MS columns (Miltenyi Biotec, 130-042-201). Post-harvesting, these neutrophils were stimulated with leucine for 24 hours. Afterward, synthesized peptides were added at a concentration of 10 μg/ml and the antigen-presenting neutrophils were allowed to uptake and process them for 12 hours. Following this, the neutrophils were cocultured with autologous T cells sorted from PBMC for 24 hours. The following neutralizing antibodies were used: TNFα (sinobiological, 10602-MM0N1), IL-6 (sinobiological, 10395-R508), IL-17 (sinobiological, 12047-M237), IL-23 (sinobiological, CT035-mh066), IFNγ (Selleck A2041). Each condition was performed with 3 replicates.
Reactive T cell response evaluated by TCR-seq
Autologous neutrophils from healthy donors’ blood were separated by staining with CD66b Biotin antibody (Biolegend, 305120), added with biotin magnetic beads, and sorted with MS columns (Miltenyi Biotec, 130-042-201). Post-separation, the neutrophils underwent stimulation with leucine for 24 hours. We cocultured the leucine-induced HLA-DR+ neutrophils with autologous T cells (negative control: T cell alone, positive control: DCs) and KRASG12V neoantigen (MTEYKLVVVGAVGVGKSALTIQLI) for 7 days from 4 donors. CD3/CD28 dynabeads were simultaneously added (4:1 to T cells). The T cells were subsequently subjected to TCR-seq. 2 μg of total RNA from each sample was utilized for the preparation of the TCR sequencing library. This process was conducted using the KC-Digital™ Stranded TCR-seq Library Prep Kit (provided by Seqhealth Technology Co., Ltd., Wuhan, China, DT0813-02), in accordance with the manufacturer's guidelines. The library products, with lengths ranging from 250-500 base pairs, were subsequently enriched and quantified. Finally, these samples were sequenced using the NovaSeq platform (Illumina). Sequences during PCR amplification were removed. Subsequently, these sequences were aligned to IMGT by using TRUST4 (https://github.com/liulab-dfci/TRUST4). The data was then analyzed and visualized utilizing the immunarch (V1.0.0, https://immunarch.com/index.html).
Cancer cell, T cell, and neutrophil coculture system
Initially, neutrophils were incubated with leucine for 24 hours. Then, the neutrophils were cocultured with autologous T cells sorted from PBMC for 24 hours. We then add the cancer cell lines (HepG2, A549, HCT116, PANC1, and MCF7) and cultured them together for 48 hours. We finally performed the flow cytometry (PI, AnexinV) to assess the apoptosis level of cancer cells. Cancer cells were gated based on SSC and FSC channels by using Flowjo software (BD). Each condition was performed with 3 replicates.
Mouse model
We obtained 5-week-old male C57BL/6 mice, CD74 KO, CD45.1 mice, and Ly6GCre-tdTomato from the Shanghai Model Organisms Center, Inc. MHC-IIflox/flox mice, LAT1KO, Bcat2KO, DbtKO mice were obtained from Nanjing GemPharmatech Co. Ltd. We housed them under pathogen-free conditions with a maximum of five mice per cage. We strictly adhered to animal care principles and ethics and received approval from the Institutional Animal Care and Use Committee of the Shanghai Model Organisms Center (approval number 2019-0011). MC38, Hepa 1-6, and LLC cells (5 × 106) were injected at day 0 subcutaneously. For the amino acid diet group, 1.5% amino acids were added to the drinking water.100 As for the PD-1 treatment group, Ultra-LEAF™ Purified anti-mouse CD279 (PD-1) (BioLegend, 135248) were injected intraperitoneally (100μg per mouse), and Ultra-LEAF™ Purified Rat IgG2a, κ Isotype Ctrl (BioLegend, 400565) was used in control groups (100μg per mouse). As for the Ly6G antibody treatment group, mice were injected with 50 μg of anti-Ly6G Ab (BE0075, BioXCell).101 As for the neutrophil adoptive delivering group, we split neutrophils from the blood of 5-week-old male C57 mice, stimulated the neutrophils with LPS and leucine for 24 h, and injected the neutrophils inside the tumor (5 × 106). Tumor volume was calculated using the formula length (mm) × widthˆ2 (mm) × 0.5. These leucine-treated and control samples were sent for single-cell RNA-seq. In each group of each cancer type, tumors were merged for sequencing. As for the lifetime of neutrophils staying in the tumor microenvironment, we first bear the Cd45.2 C57BL/6 mouse models with MC38, Hepa 1-6, and LLC cells (5 × 106) subcutaneously. We then delivered neutrophils from Cd45.1 C57BL/6 mice into Cd45.2 C57BL/6 mouse tumors, and assessed the Cd45.1+Ly6G+ cells proportion at day 0 to 5. Each condition was performed with 5 replicates.
Patient-derived tumor fragment
We followed a well-established protocol for the processing and preservation of HCC samples treated with neoadjuvant immunotherapy.62 Tumors were cut into 1–2 mm3 pieces and then gradually frozen using a gradient in cryovials containing 1 mL of freezing media (FBS with 10% DMSO). Fragments from different regions of tumor were mixed to reduce heterogeneity. After pathological assessment, samples from 5 patients with pathologic nonresponse (residual viable tumor rates of 95%, 95%, 90%, 90%, and 90% respectively; treatment: anti-PD-1 plus anti-VEGF agent) were selected for further experimentation. The tumor fragments were mixed with ice-cold matrigel (BD Biosciences; Matrix High Concentration, Phenol Red-Free, 4 mg/mL final concentration) and transferred to a 96-well flat-bottom plate. Autologous neutrophils of each patient from blood were separated by staining with CD66b Biotin antibody (Biolegend, 305120), added with biotin magnetic beads, sorted with MS columns (Miltenyi Biotec, 130-042-201), and stimulated with leucine. Autologous antigen-presenting neutrophils were added to the media and cultured for 3 days. Flow cytometry analysis of T cells was then performed. Each condition was performed with 5 replicates.
Web server
We used Apache and Shiny to construct the web server as previously described.37,102,103 We tested different functions covering common browsers including Chrome, Safari, and IE. Users are not required to register or log in to access features in the web server.
Data visualization
We utilized the following software packages for visualization: Seurat (V4.0.4), ggplot2 (V3.3.5), dittoSeq (V1.5.2), ggrepel (V0.9.1), ggpubr (V0.4.0), ggsignif (V0.6.3), pheatmap (V1.0.12), ComplexHeatmap (V2.15.1), and cowplot (V1.1.1) in R, as well as sctour (V0.1.3) and TooManyCells (V3.0.0) in Python.
Quantification and statistical analysis
We defined statistical significance as P < 0.05 and performed all statistical analyses using R (V4.1.0) and RStudio ("Elsbeth Geranium" Release). Group comparisons were conducted using Student’s t-tests, Wilcoxon rank-sum tests, and ANOVA, while paired t-tests were utilized for paired comparisons. Unless otherwise specified, bar plots were presented as mean± standard deviation. Each experiment was repeated three or more times using biologically independent samples. For correlation analyses, we used Spearman rho or Pearson r. Survival analyses were performed using log-rank tests and the proportion of HLA-DR+CD15+ neutrophils of CD15+ neutrophils and SPP1+CD15+ neutrophils of CD15+ neutrophils were used for analysis. We utilized the ggsurvplot function in the R package survminer (V0.4.9) to determine the cutoff value of proportion and generate Kaplan-Meier survival curves.
Acknowledgments
This study was supported by the National Natural Science Foundation of China (nos. 82130077, 81961128025, 82121002, and 82341008) to Q.G.; the National Natural Science Foundation of China (nos. 82394450, U23A6010) to J.F.; the National Natural Science Foundation of China (31925011) to L.Y.; The Research Projects from the Science and Technology Commission of Shanghai Municipality (nos. 21JC1410100, 21JC1401200, 20JC1418900, and 20YF1407400) to J.F., Q.G., and S.Z.; the Shanghai Municipal Science and Technology Major Project to Q.G.; the Strategic Priority Research Program (XDPB0303) of the Chinese Academy of Sciences, Program of Shanghai Academic Research Leader (23XD1404300), and the Shanghai Municipal Science and Technology Major Project (no. 2019SHZDZX02) to X.Z.; and the Science and Technology Commission of Shanghai Municipality (22ZR1479100) to S.J. We thank Claude Leclerc (Pasteur Institute), Andrew Zhu (Massachusetts General Hospital), Feng Wang, Jing Wang (Shanghai Institute of Immunology), Qunying Lei, Di Zhu, Liangqing Dong, Youpei Lin, Haichao Zhao, Xia Shen, Fanfan Fan, Shuaixi Yang, and Zijian Yang (Fudan University) for their support. We thank Shanghai Biochip Co., Ltd. for technical assistance. We thank the computing platform of the Medical Research Data Center of Shanghai Medical College Fudan University.
Author contributions
Conceptualization, J.F., X.Z., and Q.G.; methodology, Y.W., J.M., X.Y., F.N., T.Z., and S.J.; software, Y.W., J.M., X.Y., F.N., T.Z., and S.J.; validation, Y.W., J.M., X.Y., F.N., T.Z., and S.J.; investigation, Y.W., J.M., X.Y., F.N., T.Z., and S.J.; formal analysis, Y.W., J.M., X.Y., and F.N.; writing – original draft, Y.W., F.N., J.M., L.Y., X.Z., and Q.G.; writing – review & editing, Y.W., J.M., X.Y., F.N., T.Z., S.J., D.R., H.F., K.G., X.G., S.J., G.S., J.P., M.Z., Y.X., S.Z., Y.F., X.W., J.Z., L.Y., J.F., X.Z., and Q.G.; visualization, Y.W., J.M., X.Y., F.N., T.Z., and S.J.; funding acquisition, Y.W., F.N., S.J., S.Z., L.Y., J.F., X. Z., and Q.G.; resources, X.W., J.Z., L.Y., J.F., X.Z., and Q.G.; and supervision, L.Y., J.F., X.Z., and Q.G.
Declaration of interests
The authors declare no potential conflicts of interest.
Supplemental information
References
- 1The NeutrophilImmunity, 54 (2021), pp. 1377-1391
- 2Neutrophils in cancer: heterogeneous and multifacetedNat. Rev. Immunol., 22 (2022), pp. 173-187
- 3Neutrophil kinetics in health and diseaseTrends Immunol., 31 (2010), pp. 318-324
- 4Single-Cell Analysis Reveals the Range of Transcriptional States of Circulating Human NeutrophilsJ. Immunol., 209 (2022), pp. 772-782
- 5Myeloid-Derived Suppressor Cells: A Propitious Road to ClinicCancer Discov., 11 (2021), pp. 2693-2706
- 6Deterministic reprogramming of neutrophils within tumorsScience, 383 (2024), p. eadf6493
- 7Liver tumour immune microenvironment subtypes and neutrophil heterogeneityNature, 612 (2022), pp. 141-147
- 8Multidimensional immune profiling in Gastric Cancer Multiplex Immunohistochemistry Atlas from Peking University Cancer Hospital project informs PD-1/PD-L1 blockade efficacyEur. J. Cancer, 189 (2023), p. 112931
- 9Predicting response to immunotherapy in gastric cancer via multi-dimensional analyses of the tumour immune microenvironmentNat. Commun., 13 (2022), p. 4851
- 10Neutrophil elastase selectively kills cancer cells and attenuates tumorigenesisCell, 184 (2021), pp. 3163-3177.e21
- 11T cell immunotherapies engage neutrophils to eliminate tumor antigen escape variantsCell, 186 (2023), pp. 1432-1447.e17
- 12Innate Immune Training of Granulopoiesis Promotes Anti-tumor ActivityCell, 183 (2020), pp. 771-785.e12
- 13Origin and Role of a Subset of Tumor-Associated Neutrophils with Antigen-Presenting Cell Features in Early-Stage Human Lung CancerCancer Cell, 30 (2016), pp. 120-135
- 14A neutrophil response linked to tumor control in immunotherapyCell, 186 (2023), pp. 1448-1464.e20
- 15Neutrophil diversity and plasticity in tumour progression and therapyNat. Rev. Cancer, 20 (2020), pp. 485-503
- 16Estimating the population abundance of tissue-infiltrating immune and stromal cell populations using gene expressionGenome Biol., 17 (2016), p. 218
- 17Molecular and pharmacological modulators of the tumor immune contexture revealed by deconvolution of RNA-seq dataGenome Med., 11 (2019), p. 34
- 18xCell: digitally portraying the tissue cellular heterogeneity landscapeGenome Biol., 18 (2017), p. 220
- 19Pan-cancer proteogenomics connects oncogenic drivers to functional statesCell, 186 (2023), pp. 3921-3944.e25
- 20The Immune Landscape of CancerImmunity, 48 (2018), pp. 812-830.e14
- 21Conserved pan-cancer microenvironment subtypes predict response to immunotherapyCancer Cell, 39 (2021), pp. 845-865.e7
- 22High-resolution single-cell atlas reveals diversity and plasticity of tissue-resident neutrophils in non-small cell lung cancerCancer Cell, 40 (2022), pp. 1503-1520.e8
- 23Single-Cell Transcriptomics of Human and Mouse Lung Cancers Reveals Conserved Myeloid Populations across Individuals and SpeciesImmunity, 50 (2019), pp. 1317-1334.e10
- 24Atlas of clinically distinct cell states and ecosystems across human solid tumorsCell, 184 (2021), pp. 5482-5496.e28
- 25Single-cell RNA-seq analysis reveals BHLHE40-driven pro-tumour neutrophils with hyperactivated glycolysis in pancreatic tumour microenvironmentGut, 72 (2023), pp. 958-971
- 26TooManyCells identifies and visualizes relationships of single-cell cladesNat. Methods, 17 (2020), pp. 405-413
- 27CCL5 mediates CD40-driven CD4+ T cell tumor infiltration and immunityJCI Insight, 5 (2020), Article e137263
- 28The MHC class I antigen presentation pathway: strategies for viral immune evasionImmunology, 110 (2003), pp. 163-169
- 29Single-cell transcriptome profiling reveals neutrophil heterogeneity in homeostasis and infectionNat. Immunol., 21 (2020), pp. 1119-1133
- 30Programmed 'disarming' of the neutrophil proteome reduces the magnitude of inflammationNat. Immunol., 21 (2020), pp. 135-144
- 31A Neutrophil Timer Coordinates Immune Defense and Vascular ProtectionImmunity, 50 (2019), pp. 390-402.e10
- 32scTour: a deep learning architecture for robust inference and accurate prediction of cellular dynamicsPreprint atbioRxiv (2022)
- 33Reversed graph embedding resolves complex single-cell trajectoriesNat. Methods, 14 (2017), pp. 979-982
- 34Single-cell transcriptional diversity is a hallmark of developmental potentialScience, 367 (2020), pp. 405-411
- 35Slingshot: cell lineage and pseudotime inference for single-cell transcriptomicsBMC Genomics, 19 (2018), p. 477
- 36Human neutrophils in the saga of cellular heterogeneity: insights and open questionsImmunol. Rev., 273 (2016), pp. 48-60
- 37Spatiotemporal Immune Landscape of Colorectal Cancer Liver Metastasis at Single-Cell LevelCancer Discov., 12 (2022), pp. 134-153
- 38Interplay between lipids and branched-chain amino acids in development of insulin resistanceCell Metab., 15 (2012), pp. 606-614
- 39Mitochondrial function - gatekeeper of intestinal epithelial cell homeostasisNat. Rev. Gastroenterol. Hepatol., 15 (2018), pp. 497-516
- 40Mitochondrial TCA cycle metabolites control physiology and diseaseNat. Commun., 11 (2020), p. 102
- 41Extension of chemotactic pseudopods by nonadherent human neutrophils does not require or cause calcium burstsSci. Signal., 11 (2018), p. eaal4289
- 42The circadian clock influences T cell responses to vaccination by regulating dendritic cell antigen processingNat. Commun., 13 (2022), p. 7217
- 43Mammalian Mitochondrial Complex I Structure and Disease-Causing MutationsTrends Cell Biol., 28 (2018), pp. 835-867
- 44SIRT4 is an early regulator of branched-chain amino acid catabolism that promotes adipogenesisCell Rep., 36 (2021), p. 109345
- 45Toll-like Receptor Signaling Rewires Macrophage Metabolism and Promotes Histone Acetylation via ATP-Citrate LyaseImmunity, 51 (2019), pp. 997-1011.e7
- 46Glis1 facilitates induction of pluripotency via an epigenome-metabolome-epigenome signalling cascadeNat. Metab., 2 (2020), pp. 882-892
- 47NAD+ is critical for maintaining acetyl-CoA and H3K27ac in embryonic stem cells by Sirt1-dependent deacetylation of AceCS1Life Med., 1 (2022), pp. 401-405
- 48MHC class II super-enhancer increases surface expression of HLA-DR and HLA-DQ and affects cytokine production in autoimmune vitiligoProc. Natl. Acad. Sci. USA, 113 (2016), pp. 1363-1368
- 49Identification of HLA-DRB1∗09:01-restricted Mycobacterium tuberculosis CD4+ T-cell epitopesFEBS Lett., 590 (2016), pp. 4541-4549
- 50HLA-DPB1∗05: 01-restricted WT1332-specific TCR-transduced CD4+ T lymphocytes display a helper activity for WT1-specific CTL induction and a cytotoxicity against leukemia cellsJ. Immunother., 36 (2013), pp. 159-170
- 51Neoantigen screening identifies broad TP53 mutant immunogenicity in patients with epithelial cancersJ. Clin. Invest., 129 (2019), pp. 1109-1114
- 52Geospatial Immune Heterogeneity Reflects the Diverse Tumor-Immune Interactions in Intrahepatic CholangiocarcinomaCancer Discov., 12 (2022), pp. 2350-2371
- 53Enhanced detection of neoantigen-reactive T cells targeting unique and shared oncogenes for personalized cancer immunotherapyJCI Insight, 3 (2018), Article e122467
- 54Distinct Immune Cell Populations Define Response to Anti-PD-1 Monotherapy and Anti-PD-1/Anti-CTLA-4 Combined TherapyCancer Cell, 35 (2019), pp. 238-255.e6
- 55Genomic and Transcriptomic Features of Response to Anti-PD-1 Therapy in Metastatic MelanomaCell, 165 (2016), pp. 35-44
- 56TGFbeta attenuates tumour response to PD-L1 blockade by contributing to exclusion of T cellsNature, 554 (2018), pp. 544-548
- 57Immune-Related Gene Expression Profiling After PD-1 Blockade in Non-Small Cell Lung Carcinoma, Head and Neck Squamous Cell Carcinoma, and MelanomaCancer Res., 77 (2017), pp. 3540-3550
- 58Somatic Mutations and Neoepitope Homology in Melanomas Treated with CTLA-4 BlockadeCancer Immunol. Res., 5 (2017), pp. 84-91
- 59Mutational and putative neoantigen load predict clinical benefit of adoptive T cell therapy in melanomaNat. Commun., 8 (2017), p. 1738
- 60Comprehensive molecular characterization of clinical responses to PD-1 inhibition in metastatic gastric cancerNat. Med., 24 (2018), pp. 1449-1458
- 61Molecular correlates of clinical response and resistance to atezolizumab in combination with bevacizumab in advanced hepatocellular carcinomaNat. Med., 28 (2022), pp. 1599-1611
- 62An ex vivo tumor fragment platform to dissect response to PD-1 blockade in cancerNat. Med., 27 (2021), pp. 1250-1261
- 63Mapping cell types across human tissuesScience, 376 (2022), pp. 695-696
- 64A pan-cancer single-cell transcriptional atlas of tumor infiltrating myeloid cellsCell, 184 (2021), pp. 792-809.e23
- 65Pan-cancer single-cell landscape of tumor-infiltrating T cellsScience, 374 (2021), p. abe6474
- 66Co-option of Neutrophil Fates by Tissue EnvironmentsCell, 183 (2020), pp. 1282-1297.e18
- 67Oncology meets immunology: the cancer-immunity cycleImmunity, 39 (2013), pp. 1-10
- 68Dysfunction of antigen processing and presentation by dendritic cells in cancerMol. Immunol., 113 (2019), pp. 31-37
- 69Origin and Role of a Subset of Tumor-Associated Neutrophils with Antigen-Presenting Cell Features in Early-Stage Human Lung CancerCancer Cell, 30 (2016), pp. 120-135, 10.1016/j.ccell.2016.06.001IF: 50.3 Q1
- 70FcγR engagement reprograms neutrophils into antigen cross-presenting cells that elicit acquired anti-tumor immunityNat. Commun., 12 (2021), p. 4791
- 71The transcription factor RFX5 coordinates antigen-presenting function and resistance to nutrient stress in synovial macrophagesNat. Metab., 4 (2022), pp. 759-774
- 72CD40 signal rewires fatty acid and glutamine metabolism for stimulating macrophage anti-tumorigenic functionsNat. Immunol., 24 (2023), pp. 452-462
- 73A pan-cancer blueprint of the heterogeneous tumor microenvironment revealed by single-cell profilingCell Res., 30 (2020), pp. 745-762, 10.1038/s41422-020-0355-0IF: 44.1 Q1
- 74Signatures of plasticity, metastasis, and immunosuppression in an atlas of human small cell lung cancerCancer Cell, 39 (2021), pp. 1479-1496.e18, 10.1016/j.ccell.2021.09.008IF: 50.3 Q1
- 75Single-cell transcriptome analysis revealed a suppressive tumor immune microenvironment in EGFR mutant lung adenocarcinomaJ. Immunother. Cancer, 10 (2022), p. e003534, 10.1136/jitc-2021-003534IF: 10.9 Q1
- 76TDO2+ myofibroblasts mediate immune suppression in malignant transformation of squamous cell carcinomaJ. Clin. Invest., 132 (2022), p. e157649, 10.1172/JCI157649IF: 15.9 Q1
- 77The Tabula Sapiens: A multiple-organ, single-cell transcriptomic atlas of humansScience, 376 (2022), p. eabl4896, 10.1126/science.abl4896IF: 56.9 Q1
- 78Integrating single-cell transcriptomic data across different conditions, technologies, and speciesNat. Biotechnol., 36 (2018), pp. 411-420
- 79Fast, sensitive and accurate integration of single-cell data with HarmonyNat. Methods, 16 (2019), pp. 1289-1296
- 80Pan-cancer single-cell RNA-seq identifies recurring programs of cellular heterogeneityNat. Genet., 52 (2020), pp. 1208-1218
- 81Normalization and variance stabilization of single-cell RNA-seq data using regularized negative binomial regressionGenome Biol., 20 (2019), p. 296
- 82UCell: Robust and scalable single-cell gene signature scoringComput. Struct. Biotechnol. J., 19 (2021), pp. 3796-3798
- 83KEGG: new perspectives on genomes, pathways, diseases and drugsNucleic Acids Res., 45 (2017), pp. D353-D361
- 84The Molecular Signatures Database (MSigDB) hallmark gene set collectionCell Syst., 1 (2015), pp. 417-425
- 85Gene ontology: tool for the unification of biology. The Gene Ontology ConsortiumNat. Genet., 25 (2000), pp. 25-29
- 86Zero-preserving imputation of single-cell RNA-seq dataNat. Commun., 13 (2022), p. 192
- 87CellPhoneDB: inferring cell-cell communication from combined expression of multi-subunit ligand-receptor complexesNat. Protoc., 15 (2020), pp. 1484-1506
- 88Visualizing and interpreting cancer genomics data via the Xena platformNat. Biotechnol., 38 (2020), pp. 675-678
- 89Profiling Tumor Infiltrating Immune Cells with CIBERSORTMethods Mol. Biol., 1711 (2018), pp. 243-259
- 90Inferring tumour purity and stromal and immune cell admixture from expression dataNat. Commun., 4 (2013), p. 2612
- 91Pan-cancer Immunogenomic Analyses Reveal Genotype-Immunophenotype Relationships and Predictors of Response to Checkpoint BlockadeCell Rep., 18 (2017), pp. 248-262
- 92Comprehensive analyses of tumor immunity: implications for cancer immunotherapyGenome Biol., 17 (2016), p. 174
- 93Simultaneous enumeration of cancer and immune cell types from bulk tumor gene expression dataeLife, 6 (2017)
- 94The UCSC Genome Browser database: 2022 updateNucleic Acids Res., 50 (2022), pp. D1115-D1122
- 95Multimodule characterization of immune subgroups in intrahepatic cholangiocarcinoma reveals distinct therapeutic vulnerabilitiesJ. Immunother. Cancer, 10 (2022)
- 96Single-cell transcriptomic analysis suggests two molecularly subtypes of intrahepatic cholangiocarcinomaNat. Commun., 13 (2022), p. 1642
- 97Distribution and density of tertiary lymphoid structures predict clinical outcome in intrahepatic cholangiocarcinomaJ. Hepatol., 76 (2022), pp. 608-618
- 98Single-cell profiling reveals distinct adaptive immune hallmarks in MDA5+ dermatomyositis with therapeutic implicationsNat. Commun., 13 (2022), p. 6458
- 99Serum Metabolomics Identifies Dysregulated Pathways and Potential Metabolic Biomarkers for Hyperuricemia and GoutArthritis Rheumatol., 73 (2021), pp. 1738-1748
- 100Leucine-tRNA-synthase-2-expressing B cells contribute to colorectal cancer immunoevasionImmunity, 55 (2022), p. 1748
- 101An Inflammatory Checkpoint Generated by IL1RN Splicing Offers Therapeutic Opportunity for KRAS-Mutant Intrahepatic CholangiocarcinomaCancer Discov., 13 (2023), pp. 2248-2269
- 102SPACE: a web server for linking chromatin accessibility with clinical phenotypes and the immune microenvironment in pan-cancer analysisCell. Mol. Immunol., 17 (2020), pp. 1294-1296
- 103Multi-omics analysis reveals the functional transcription and potential translation of enhancersInt. J. Cancer, 147 (2020), pp. 2210-2224