条文本
PDF
摘要
客观的先天免疫在胰腺导管腺癌(PDAC)这种非t细胞富集的肿瘤中起着重要作用。中性粒细胞是先天免疫系统的主要参与者。在这里，我们旨在探索PDAC中中性粒细胞的异质性和促肿瘤机制。
设计我们分析了5例PDAC患者的外周血多形核白细胞(PMNs)和肿瘤浸润性免疫细胞的单细胞转录组，并进行了免疫荧光/免疫组化染色、多组学分析在体外验证生物信息学分析发现的实验。
结果对肿瘤相关中性粒细胞(TANs)异质性的探索揭示了一种与预后不良相关的最终分化的前肿瘤亚群(TAN-1)、一种炎症亚群(TAN-2)、一种刚刚迁移到肿瘤微环境的过渡阶段的亚群(TAN-3)和一种优先表达干扰素刺激基因的亚群(TAN-4)。糖酵解信号沿中性粒细胞转变轨迹上调，TAN-1具有糖酵解活性高度激活的特征。通过转录组学、蛋白质组学和代谢组学的综合多组学方法验证了TANs的糖酵解开关。LDHA过表达激活糖酵解活性诱导中性粒细胞样分化HL-60 (dHL-60)细胞的免疫抑制和促肿瘤功能。机制研究显示，BHLHE40在缺氧和内质网应激下游是中性粒细胞向TAN-1表型极化的关键调节因子，通过染色质免疫沉淀试验证实BHLHE40对TAN-1标记基因的直接转录调控。在过表达BHLHE40的dHL-60细胞中观察到促肿瘤和免疫抑制功能。重要的是，PDAC组织的免疫组化分析揭示了BHLHE40的不良预后价值+中性粒细胞。
结论本研究揭示的TANs的动态特性将有助于推进针对先天免疫的PDAC治疗。
- 胰腺癌
- 免疫疗法
- 葡萄糖代谢
数据可用性声明
数据可在一个公共的、开放访问的存储库中获得。本研究产生的单细胞和大块RNA测序数据保存在国家组学数据百科(NODE)，登录代码为OEP003254。
这是一篇开放获取的文章，按照创作共用署名非商业性(CC BY-NC 4.0)许可发布，该许可允许其他人以非商业性的方式发布、混编、改编、构建本作品，并以不同的条款授权他们的衍生作品，前提是原创作品被正确引用，给予适当的荣誉，任何更改都被注明，且使用是非商业性的。看到的:http://creativecommons.org/licenses/by-nc/4.0/．
来自Altmetric.com的统计
关于这个话题我们已经知道了什么
作为非T细胞富集的肿瘤，胰导管腺癌(PDAC)对当前免疫检查点封锁或嵌合抗原受体T的免疫治疗策略反应有限。
在PDAC肿瘤微环境中发现大量肿瘤相关中性粒细胞(TANs)浸润，这被认为是大多数实体肿瘤的一个重要的不利预后因素。
中性粒细胞是先天免疫的重要分子，通过促进细胞增殖、血管生成、组织重塑、免疫抑制和转移等途径促进肿瘤进展。
中性粒细胞是一种脆弱的细胞，在以往的单细胞rna测序研究中未被发现，该研究显示了人类PDAC肿瘤微环境中主要细胞类型的综合基因表达图谱。
这项研究补充了什么
与来自外周血的中性粒细胞相比，TANs是由异质群体组成的。
我们在PDAC肿瘤微环境中观察到了四个具有不同特征的中性粒细胞亚群，首次揭示了中性粒细胞在单细胞水平上的异质性。
在PDAC环境中，具有高活性糖酵解活性和促肿瘤功能的终末分化中性粒细胞亚群与患者较差的预后相关。
PDAC中糖酵解活性沿中性粒细胞转变轨迹上调。
糖酵解开关促进中性粒细胞的免疫抑制和促肿瘤功能。
BHLHE40在缺氧和内质网应激下被激活，是一种关键的调节因子，驱动中性粒细胞向促肿瘤和免疫抑制亚型发展。
这项研究将如何影响研究、实践和/或政策
我们的工作提供了PDAC肿瘤微环境中性粒细胞的全面图谱。
在我们的研究中证明的驱动中性粒细胞向肿瘤前亚型转化的分子机制将促进针对PDAC患者的肿瘤前亚型或TANs免疫代谢的新免疫疗法的开发，为目前主要基于细胞毒性T细胞作用的免疫治疗提供替代选择。
简介
胰腺癌，主要由胰导管腺癌(PDAC)组成，是最致命的恶性肿瘤之一，也是癌症相关死亡的主要原因之一，5年生存率仅为9%。1免疫检查点封锁或嵌合抗原受体T的免疫治疗策略在PDAC中疗效有限，因为高基质密度对T细胞浸润产生了物理屏障。2 3在PDAC的免疫环境中观察到最小的抗肿瘤t细胞浸润。4在非t细胞浸润的肿瘤中，它需要触发先天免疫激活，从而促进效应t细胞运输的信号和通向适应性免疫的桥梁。5个6因此，针对先天免疫系统的治疗干预可能为目前PDAC的免疫治疗提供另一种选择。
与其他实体肿瘤相似，PDAC的特点是在肿瘤微环境中大量浸润免疫抑制细胞。3 7 8纤维炎症滤液被癌细胞培养，提供一个有利的微环境，支持肿瘤细胞的免疫逃逸、恶性转化和进展。3 7 8在各种类型的肿瘤浸润免疫细胞中，中性粒细胞在肿瘤进展中的重要作用在过去十年中引起了广泛的研究兴趣。9中性粒细胞是一种高度可塑性的种群，在各种环境线索的反应中呈现出异质表型。9日10n2极化的中性粒细胞是肿瘤微环境中绝大多数中性粒细胞，通过促进肿瘤细胞增殖、血管生成、组织重塑、免疫抑制和转移支持肿瘤进展。9日10肿瘤相关中性粒细胞浸润(TANs)已被认为是大多数实体性恶性肿瘤的不利预后因素。11然而，仍然有证据表明中性粒细胞可以阻止肿瘤的发展，因为它们对肿瘤细胞表现出直接的细胞毒性12日13或者抗体依赖性细胞毒性。14因此，肿瘤中TANs的特征尚不明确。
在之前的一项研究中，对TCGA队列中大量rna -测序(RNA-seq)数据的生物信息学分析显示，与其他癌症类型相比，PDAC中TAN的浸润相对丰富，这暗示了TANs在PDAC微环境中的重要性。15然而，以往大多数使用单细胞测序方法的研究未能从PDAC组织中检测到人中性粒细胞。月16日我们假设，这是由于中性粒细胞作为纤弱细胞具有较短的寿命所导致的系统性偏差。22中性粒细胞衍生颗粒的高敏感性和多种酶是其易碎性的原因在体外．22因此，中性粒细胞可能无法在组织消化或单细胞捕获过程中存活。使用BD Rhapsody的基于微井的单细胞捕获方法为描述人类中性粒细胞提供了另一种选择。23日24
在这项研究中，我们使用BD Rhapsody成功地从5名PDAC患者中获得了21 972个中性粒细胞的单细胞转录组。我们探讨了PDAC微环境中TANs的异质性，并通过PDAC组织的免疫荧光(IF)和免疫组化(IHC)染色进一步验证了新发现的中性粒细胞亚型的存在和临床相关性。此外，我们发现了驱动TANs向肿瘤前表型转化的潜在机制，并确定了与中性粒细胞相关的潜在治疗靶点。
材料和方法
临床样本的处理
对于外周血，红细胞通过重力沉降消耗，中性粒细胞用cd66b - phycollythrin (PE)单克隆抗体(BioLegend, 305106)标记，根据制造商的协议，用抗PE珠和MACS柱(Miltenyi)分离。用肿瘤分离试剂盒将PDAC肿瘤组织消化成单细胞悬液(Miltenyi, 130-095-929)。浸润性免疫细胞采用CD45微珠和MACS柱(Miltenyi)分离，浸润性中性粒细胞采用与外周血相同的方法分离。患者登记、样本处理和纯度评估的详细信息载于在线补充方法．
单细胞rna测序和数据分析
单细胞捕获由BD Rhapsody系统实现。根据BD Rhapsody单细胞全转录组扩增工作流程制备整个转录组文库，并使用HiSeq Xten (Illumina, San Diego, California, USA)在150 bp对端运行进行测序。原始数据处理使用fastp过滤适配器序列和删除低质量的读取。25利用ui -tools识别小区条码白名单。26利用STAR算法将基于ui的清洁数据映射到人类参考基因组(Ensemble V.91)。27为每个样本生成UMI计数矩阵，并导入Seurat R工具包(V.3.2.3)。28聚类分析由Seurat进行，28伪时间轨迹分析使用Monocle2包(V.2.18.0)进行，29用CellPhoneDB分析细胞与细胞之间的通讯30.利用单细胞调控网络推理聚类(SCENIC)分析构建基因调控网络。31单细胞测序和数据分析的详细描述在在线补充方法．
结果
来自外周血和PDAC肿瘤的中性粒细胞的单细胞图
基于我们中心的PDAC队列，我们观察到PDAC组织中恶性细胞内部或邻近区域有大量TANs浸润(在线补充图S1A-B中性粒细胞浸润较高的PDAC患者预后不良(在线补充图S1C)．为了全面编目肿瘤微环境中中性粒细胞及其与其他免疫滤液的相互作用，我们生成了CD66b的单细胞RNA-seq数据+外周血多形核白细胞(pmn)和CD45+来自5名treatment-naïve PDAC患者的肿瘤浸润免疫细胞。来自健康对照组和慢性胰腺炎患者的pmn的转录特征也与对照组相同(图1一个，左面板和在线补充表S1)．为了验证基于单细胞RNA-seq的发现，我们收集了成对的CD66b+从24例PDAC患者的外周血和肿瘤组织(分别为pmn和TANs)中提取细胞用于多组学分析和定量PCR (qPCR)分析(验证队列1)，并招募另外114例患者进行IHC、IF和空间转录组学研究(验证队列2)(图1一个右面板,在线补充图S1D而且在线补充表S2)．
来自外周血和PDAC肿瘤的单细胞图谱，以及PDAC微环境中中性粒细胞和其他细胞之间的细胞-细胞通讯。(A)实验流程图解，包括初级队列的单细胞测序研究，验证队列1的多组学分析和定量PCR分析，验证队列2的IHC、IF和ST研究。每个亚组的患者/健康供体数量见图。(B) UMAP图描述了通过单细胞测序确定的主要细胞类型，包括CD66b+来自健康对照组、慢性胰腺炎患者和PDAC患者外周血的pmn和CD45+PDAC肿瘤组织的免疫滤液。(C)点图描述PDAC肿瘤微环境中中性粒细胞和其他细胞之间的配体-受体相互作用。平均值和p值由CellPhoneDB (http:// www.cellphonedb.org/)，并分别用圆圈的颜色和大小表示。(D)柱状图描述了PDAC肿瘤微环境中中性粒细胞和其他细胞之间显著的配体-受体相互作用的数量。细胞类型标签写为(表达配体的细胞类型)→(表达相应受体的细胞类型)。(E) PDAC组织中CD66b和CD68的IF染色。(F)显示TANs中TNFα和IL-1通路富集分数的基因集富集分析(GSEA)图。基于PDAC患者pmn和TANs的单细胞RNA测序(RNA-seq)数据进行路径富集分析。CP,慢性胰腺炎;如果,免疫荧光;包含IHC,免疫组织化学;IL,白介素; HC, healthy controls; NES, normalised enrichment score; NF-κB, nuclear factor kappa B; PDAC, pancreatic ductal adenocarcinoma; PMNs, polymorphonuclear leucocytes; ST, spatial transcriptomics; TNF, tumour necrosis factor; UMAP, uniform manifold approximation and projection.
经过初步质量控制，我们从共51 823个细胞中获得96 692 883个独特的转录本，其中33 891个细胞来自外周血，17 932个细胞来自PDAC肿瘤组织(在线补充表S3)．通过主成分分析和基于图的聚类分析，我们根据已知的细胞类型标记基因(图1 b，在线补充图S2A-D，在线补充表S4)．我们发现了两种不同类型的中性粒细胞，其中中性粒细胞1在外周血中强烈富集(主要由pmn组成)，而中性粒细胞2只存在于肿瘤组织中(由TANs组成)(图1 b，在线补充图S2E)．
为了比较来自不同来源的中性粒细胞的表达谱，我们通过t-SNE分析可视化了来自每个样本的中性粒细胞的转录组，该方法将高维数据映射到二维空间，并保留了输入对象的成对相似性。32根据t-SNE图，TANs具有独特的RNA谱，而来自PDAC患者的pmn与来自健康对照或慢性胰腺炎患者的pmn非常相似(在线补充图S3A-B)，表明中性粒细胞被招募到肿瘤组织后发生了显著的转录重编程。途径富集分析显示，与健康对照/慢性胰腺炎患者相比，PDAC患者的pmn中中性粒细胞分化受到抑制，而在TANs中观察到抑制的先天免疫、激活的多种癌症相关信号通路和深刻的代谢变化(在线补充图S3C)．
细胞-细胞通讯分析显示肿瘤微环境中中性粒细胞和巨噬细胞之间存在密切的串扰
为了探索肿瘤微环境中中性粒细胞和其他免疫细胞之间的相互作用，我们应用了一个交互式web应用程序CellPhoneDB，该应用程序根据配体-受体信号数据库推断细胞间的相互作用。30.值得注意的是，与其他类型的免疫细胞相比，巨噬细胞表达的与中性粒细胞的配体对应的受体数量显著增加，相反，巨噬细胞也表达的与中性粒细胞表达的受体对应的配体数量显著增加(图1 c, D)，表明肿瘤微环境中中性粒细胞与巨噬细胞之间相互作用密切。我们的数据表明，肿瘤相关巨噬细胞通过CCL13-CCR1, CCL3-CCR1, CCL3L3-CCR1, CXCL2-CXCR1, CXCL2-CXCR2和CXCL8-CXCR2轴吸引中性粒细胞，而中性粒细胞通过CCL3L3-CCR1轴(图1 c)．与趋化因子和趋化因子受体的表达所暗示的相互招募一致，来自癌症基因组图谱-胰腺癌(TCGA-PAAD)队列的treatment-naïve PDAC中巨噬细胞和中性粒细胞信号正相关(在线补充图S4A)．此外，IF染色证实了PDAC组织中中性粒细胞和巨噬细胞之间的物理相似性(图1 e)．值得注意的是，巨噬细胞通过促炎细胞因子白细胞介素1 (IL-1)和肿瘤坏死因子(TNF)刺激中性粒细胞(图1 c与pmn相比，TANs中TNFα/核因子kappa B (NF-κB)和IL-1信号通路显著激活(图1 f)．有趣的是，PDAC组织的IF染色验证了TANs中NF-κB和MAPK信号通路的激活(由核中p65和c-JUN的存在所揭示)下游的TNFα和IL-1 (在线补充图S4B)， qPCR分析也显示多个TNFα和IL-1靶基因表达上调(在线补充图S4C)，到三十五证实了TNFα和il -1诱导的TANs激活。
来自PDAC肿瘤的TANs是由一个异质群体组成的
为了进一步研究中性粒细胞的异质性，我们对PDAC患者的pmn和TANs进行了降维和聚类。未经批效应校正的数据见在线补充图S5A，其中pmn和tan都根据它们的起源样本聚类。校正批效应后，识别出6个PMN子簇和6个TAN子簇(图2 a, B，在线补充图S5B、C，在线补充表S5、S6)，并且在所有PDAC患者样本中都存在中性粒细胞亚簇(在线补充图S5D)．由于嗜酸性粒细胞标记物Charcot-Leyden晶体凝集素优先表达，PMN-5被排除在外(在线补充图S5C)．TAN-5由低质量细胞组成，也被忽略用于进一步分析(在线补充图S5E)．相关分析显示，所有PMN亚群具有相似的表达谱，而TAN亚群之间的差异更大(图2 c)，表明肿瘤浸润性中性粒细胞是一个异质性种群，这可能归因于肿瘤微环境中对不同环境刺激的响应的不同表型。
来自PDAC肿瘤的TANs是由一个异质群体组成的。(A) UMAP图显示PDAC患者的pmn亚簇。(B)显示PDAC患者TANs亚簇的UMAP图。所有来自PDAC肿瘤组织的注释为中性粒细胞的细胞被分析为TANs，包括中性粒细胞簇1和2图1 b．(C)各亚群中性粒细胞表达谱之间的Pearson相关系数。计算各子聚类中前2000个高变异基因的归一化平均表达量，并进行Pearson相关分析。(D)显示(i) TAN-1亚簇、(ii) TAN-2亚簇、(iii) TAN-3亚簇和(iv) TAN-4亚簇标记基因的小提琴图。(E)不同中性粒细胞亚簇间通路活性的比较。这些通路与中性粒细胞功能相关。通过基因组变异分析对每个细胞的通路活性进行评分。T值和p值的计算基于线性模型，分析来自一个团簇的中性粒细胞与来自所有其他团簇的中性粒细胞之间的差异，并分别用圆圈的颜色和大小表示。(F)中性粒细胞沿二维空间伪时间运动轨迹。每个点对应一个单元格。 (G) Heatmap showing the dynamic changes of gene expression along pseudotime. The differentially expressed genes were clustered hierarchically into three groups, and the representative enriched pathways of each group were shown. (H) Expression of TAN-1 marker genes in neutrophil-like dHL-60 cells stimulated with THG, IL-1β, TNFα and hypoxia for 24 hours, analysed by qPCR. qPCR data were normalised to fold over β-actin (housekeeping gene), and represented as mean with SD. *P<0.05; **p<0.01; ***p<0.001. dHL-60, differentiated HL-60; ER, endoplasmic reticulum; HYP, hypoxia; IL, interleukin; NF-κB, nuclear factor kappa B; PDAC, pancreatic ductal adenocarcinoma; PMN, polymorphonuclear leucocytes; qPCR, quantitative PCR; TAN, tumour-associated neutrophils; THG, thapsigargin; TNF, tumour necrosis factor; UMAP, uniform manifold approximation and projection.
考虑到TAN的异质性，我们进一步探讨了每个TAN子簇的特征。值得注意的是，在TAN-0中没有发现集群特异性的显著特征，因为TAN-0的标记基因在其他TAN亚集群中也有高表达(在线补充图S5C)．TAN-1可以被识别为“促肿瘤亚群”，它高度表达促血管生成因子VEGFA，36PLAU pro-metastatic因素37 38LGALS3是一种具有多种促进肿瘤进展功能的分子，包括增强恶性细胞的增殖和干性、促进血管生成、抑制免疫监视、促进M2巨噬细胞分化和促进耐药39 40（图2 d(我))。TAN-2由一个炎症亚群组成，该亚群强烈表达炎症相关基因NLRP3和PDE4B，41 42中性粒细胞激活标记CD6943 44（图2 d(2))。TAN-2也强烈表达促肿瘤分子IL1RN和肾上腺髓质素(ADM) (图2 d(ii))，已被发现支持其他癌症类型的肿瘤进展。45 46在所有的TAN亚簇中，TAN-3的基因表达谱与pmn最相似(图2 c中性粒细胞跨内皮迁移相关的TAN-3高表达基因(VNN2和SELL，图2 d(3)),47 48这表明TAN-3是一种刚刚迁移到肿瘤微环境的过渡阶段中性粒细胞，从pmn转化为tanan。TAN-4具有独特的转录特征，表达干扰素(IFN)刺激的基因，包括IFIT1、IFIT2、IFIT3、ISG15和RSAD2 (图2 d(四))。49
为了描述每个中性粒细胞亚簇的功能，我们用基因集变异分析分析了通路活性50（图2 e)．pmn高表达与中性粒细胞正常免疫功能相关的基因，包括先天免疫反应、吞噬作用、呼吸爆发和中性粒细胞颗粒合成。同样，TAN-3与吞噬、呼吸爆发和中性粒细胞颗粒相关的基因也高度表达，TAN-4与先天免疫应答相关。然而，与pmn相比，TAN-3和TAN-4优先表现出较高的抗原加工和呈递活性。TAN-1和TAN-2同样表达高水平的趋化因子、细胞因子和血管生成因子。TAN-2也表现出强烈的炎症特征表达，与上述炎症基因表达一致。
轨迹分析显示中性粒细胞最终分化为肿瘤前TAN-1状态
为了进一步研究中性粒细胞从外周血进入肿瘤微环境的动态过渡过程，我们使用了monocle229建立中性粒细胞状态轨迹的伪时间图。轨迹确定为以pmn为起点，经过TAN-3作为pmn和TANs之间的过渡状态，然后是以TAN-0和TAN-4为特征的中间肿瘤浸润状态，最后达到TAN-2和TAN-1的终分化状态(图2 f)．接下来，我们对沿着转录轨迹的单细胞转录组进行了分析，发现了1757个表达发生显著变化的基因，可分为三种表达模式:第1组是沿着转录轨迹表达水平下降的基因。途径富集分析显示，这些基因与IFN信号通路和先天免疫功能相关。第2组基因在肿瘤浸润早期上调，参与了吞噬和抗原呈递。第3组基因在肿瘤微环境晚期被激活，在缺氧、内质网(ER)应激、IL-1和TNF信号通路及糖酵解等条件下富集(图2 g)．为了进一步探讨这些信号通路的激活如何影响中性粒细胞表型，我们将内质网应激诱导剂THG (THG)、促炎因子IL-1β和tnf - α和缺氧处理中性粒细胞样分化HL-60细胞(dHL-60)在体外．有趣的是，我们观察到THG和缺氧都导致TAN-1标记表达显著上调(图2 h)，表明内质网应激和缺氧是TAN-1极化的有效刺激因子。
TAN-1与PDAC患者预后较差相关
为了验证PDAC组织中这些新发现的TAN亚簇，我们进行了IF染色，并证实了每个TAN亚簇(VEGFA, NLRP3, MME和IFIT2)中存在表达标记基因的中性粒细胞(图3一)．为了进一步研究这些TAN亚型的空间分布，我们使用免疫组化技术(IHC)对一系列肿瘤切片进行CD66b和TAN亚型标记基因的染色(图3 b，在线补充图S6A)．VEGFA的中位数百分比+晒得NLRP3,+黝黑色,居里夫人+鞣料和IFIT2+PDAC组织中TANs占细胞总数的比例分别为2.6%、1.4%、4.7%和0.9%。尤其是大部分的VEGFA+TANs在空间上接近恶性细胞，使其能够执行促肿瘤功能，而IFIT2+TANs常出现在纤维化间质组织中(图3 b)．进一步的临床资料分析表明，VEGFA+TANs与美国癌症联合委员会(AJCC)晚期显著相关(在线补充图S6B)和较差的预后(图3 c，在线补充图S6C)。
TAN亚簇的空间分布和临床相关性。(A) PDAC组织中中性粒细胞(CD66b)和TAN亚簇(VEGFA为TAN-1, NLRP3为TAN-2, MME为TAN-3, IFIT2为TAN-4)标记基因的IF染色。(B)连续载片上代表性PDAC组织的免疫组化图像，染色为中性粒细胞(CD66b)标记基因和TAN亚簇(VEGFA为TAN-1, NLRP3为TAN-2, MME为TAN-3, IFIT2为TAN-4)。根据连续载玻片CD66b染色和多核形态鉴定为中性粒细胞。粉色箭头显示表达TAN亚簇标记的中性粒细胞。(C) Kaplan-Meier生存曲线，显示免疫组化分析中PDAC患者的总体生存期。根据VEGFA的百分比将患者平均分为两组+PDAC组织中总细胞的TANs。(D) TCGA-PAAD队列中具有标准PDAC组织学的treatment-naïve例患者的总体生存期Kaplan-Meier生存曲线。根据TAN-1标记基因表达情况将患者平均分为两组，采用gsa对各样本TAN-1标记基因表达进行评估。(E) Kaplan-Meier生存曲线，显示了来自Cao发表的数据集的treatment-naïve例PDAC标准组织学患者的总生存期等．根据TAN-1标记基因表达情况将患者平均分为两组，采用gsa对各样本TAN-1标记基因表达进行评估。GSVA，基因集变异分析;HPF，高功率场;包含IHC,免疫组织化学;如果,免疫荧光;LPF，低功率场;PDAC，胰导管腺癌;TAN，肿瘤相关中性粒细胞。
为了进一步确认基于公共数据库的TAN-1亚群的临床相关性，我们检索了TCGA-PAAD队列的RNA-seq和临床数据，51以及曹发表的数据集等．52与我们队列的结果一致，TAN-1信号是这两个队列中treatment-naïve例PDAC组织学标准患者的不利预后因素(图3 d e)．
代谢分析显示，糖酵解开关在TAN-1中主要上调
考虑到如上所述肿瘤微环境中嗜中性粒细胞的深刻代谢重编程(在线补充图S3C)，以及代谢计划对免疫细胞功能的重大影响，53 54我们分析了每个TAN亚簇的标志性代谢特征，发现TAN-1的糖酵解和缺氧活性明显高于其他中性粒细胞亚簇(图4一，在线补充图S7A,B)，并沿中性粒细胞转变轨迹上调(图4 b)．流式细胞术分析显示，从PDAC组织分离的表达TAN-1标记物LGALS3的中性粒细胞中，葡萄糖转运体GLUT1和糖酵解酶HK2、PFKFB3和LDHA的表达上调(图4 c，在线补充图S7C,D)．为了进一步探索TANs的代谢特征是否与其空间分布相关，我们从PDAC组织切片上的2498个点(在线补充图S7E、F)，平均每个斑点约检测到13 095个uis和3969个独特基因。有趣的是，与间质区中性粒细胞富集点相比，肿瘤区域内或邻近中性粒细胞富集点的糖酵解活性显著上调(图4 d, E)，暗示PDAC肿瘤微环境中恶性细胞诱导的中性粒细胞糖酵解开关。
代谢分析显示，糖酵解开关在TAN-1中主要上调。(A)标记糖酵解和缺氧标记在每个中性粒细胞亚簇中的表达箱图(左图)和标记基因的平均归一化表达热图(右图)。(B)糖酵解特征和缺氧特征沿伪时间轨迹动态表达的二维图。(C)流式细胞术分析GLUT1、HK2、PFKFB3和LDHA在LGALS3中的表达+和LGALS3−来自PDAC组织的晒黑。LGALS3中GLUT1、HK2、PFKFB3和LDHA表达的MFI+和LGALS3−在4例PDAC患者的肿瘤组织中测量TANs，并以SD患者的平均值表示。(D)肿瘤、邻近肿瘤和间质区中性粒细胞富集点的空间特征图。根据单细胞测序鉴定的中性粒细胞标记基因，通过单样本基因集富集分析(ssGSEA)计算各空间点的中性粒细胞特征，将中性粒细胞特征值前10%的点(ssGSEA评分>0.53)定义为中性粒细胞富集点。(E)盒图显示肿瘤、邻近肿瘤和间质区中性粒细胞富集点的糖酵解特征。(F-G) GSEA图显示TANs中HALLMARK糖酵解途径的NES。根据配对CD66b的转录组学(F)和蛋白质组学(G)数据进行路径富集分析+PDAC患者的pmn和TANs。(H) PMN和TAN裂解物中糖酵解中间体的代谢组学分析。数据用均值和标准差表示。(一)糖酵解途径示意图。红色框表示与PDAC患者的pmn相比，TANs中的代谢物显著增加(p<0.05)。MFI，平均荧光强度;NES，归一化浓缩分数;PDAC，胰导管腺癌;中性粒细胞,多形核白细胞;棕褐色,肿瘤相关中性粒细胞。 *P<0.05.
为了进一步验证TANs的糖酵解开关，我们从19例PDAC患者中分离出成对的pmn和TANs (图1一个)，并应用综合多组学方法比较其在转录组、蛋白质组和代谢组水平上的糖酵解活性。对大量RNA-seq数据的分析显示，pmn和TANs之间有3920个基因的差异表达(在线补充表S7)，途径富集分析验证了TANs中糖酵解的上调(图4 f而且在线补充表S8)．同样，根据基于数据独立采集质谱(DIA-MS)的定量蛋白质组学分析所识别的蛋白表达谱(在线补充表S9)，糖酵解途径在TANs中显著富集(图4 g)．对中性粒细胞裂解物的代谢组学分析也揭示了TANs中几种糖酵解中间体的显著上调(图4 h,我)．
糖酵解开关增强TANs的促肿瘤功能
接下来，我们研究了中性粒细胞的糖酵解开关和促肿瘤功能之间的关系。我们关注的是LDHA，它是糖酵解途径(图4一)，也是TAN-1的标记基因之一(图2 d(我))。代谢分析显示LDHA在dHL-60细胞中过表达(在线补充图S8A、B)导致葡萄糖消耗增加(在线补充图S8C)和乳酸生产(在线补充图S8D)，证实了ldha -过表达dHL-60中糖酵解活性的上调。与过表达ldha的dHL-60共培养的PDAC肿瘤细胞(PATU-8988或Aspc-1)增殖能力增强(在线补充图S8E-H)，揭示了中性粒细胞的代谢重编程导致的促肿瘤效应。先前的研究表明，肿瘤微环境中ldha相关乳酸的积累抑制了抗肿瘤免疫监视，55这促使我们研究LDHA过表达对中性粒细胞免疫抑制功能的影响。有趣的是，我们观察到与过表达ldha的dHL-60细胞共培养会降低CD8中IFNγ和TNFα的表达+用phorbol 12-肉豆蔻酸13-醋酸酯(PMA)和离子霉素刺激T细胞，以及对CD3/CD28激活的T细胞增殖轻微抑制(在线补充图S8I,J)．
上游调控因子分析显示BHLHE40驱动中性粒细胞向肿瘤前表型发展
风景优美的31用于评估驱动TANs异质性的上游转录因子(规则)(在线补充图S9A)．我们确定了4个具有独特调控活性的TAN簇，与基于RNA谱的原始TAN子簇划分相关(在线补充图S9B、C)．值得注意的是，TAN-1的标记基因之一BHLHE40是TAN-1中最显著上调的转录因子之一(图5一个，在线补充图S9A)．根据SCENIC构成的基因调控网络，BHLHE40是转录因子HIF1A和XBP1的下游靶基因，BHLHE40调控TAN-1标记VEGFA、PLAU、LGALS3、LDHA和BHLHE40的表达(自调控)，以及TAN-2标记PDE4B和IL1RN (图5 b)，提示在肿瘤微环境中，缺氧和内质网应激诱导的BHLHE40促进中性粒细胞向肿瘤前表型分化。
上游调控因子分析发现BHLHE40是TAN-1的关键转录因子。(A)显示BHLHE40在TANs中的表达的UMAP图(左图)和估计BHLHE40在TANs中的调节活性的曲线下面积(AUC)评分(右图)。(B)单细胞调控网络推理和聚类(SCENIC)揭示的BHLHE40代表性调控网络。(C)通过qPCR分析BHLHE40靶蛋白在对照细胞和BHLHE40过表达dHL-60细胞中的表达情况。每个基因的mRNA表达水平归一化到β-actin(管家基因)的水平。数据用均值和标准差表示。(D) dHL-60细胞中BHLHE40在靶基因启动子区结合，用染色质免疫沉淀- qpcr分析。(E) THG和/或HYP刺激dHL-60细胞24小时后BHLHE40靶点的表达和BHLHE40敲除的dHL-60细胞的表达，用qPCR分析。每个基因的mRNA表达水平归一化到β-actin(管家基因)的水平。数据用均值和标准差表示。 *P<0.05; **p<0.01; ***p<0.001; ****p<0.0001. dHL-60, neutrophil-like differentiated HL-60; HYP, hypoxia; KD, knockdown; mRNA, messenger RNA; NC, negative control; OE, overexpression; THG, thapsigargin; TSS, transcription start site; qPCR, quantitative PCR; TAN, tumour-associated neutrophils; UMAP, uniform manifold approximation and projection.
为了进一步验证BHLHE40的调节作用，我们分析了BHLHE40过表达对中性粒细胞样dHL-60细胞表达谱的影响。qPCR证实，BHLHE40过表达导致VEGFA、PLAU、LGALS3、LDHA和PDE4B的表达显著上调(图5 c)．染色质免疫沉淀- qpcr测定也证实BHLHE40与这些基因的启动子区结合(图5 d)，证明了中性粒细胞中BHLHE40对促肿瘤基因的直接转录调控。有趣的是，与SCENIC推断的调节网络一致，我们观察到缺氧和内质网应激显著诱导TAN-1标记的表达，且这两个因素有很强的协同作用，而在BHLHE40敲低的dHL-60细胞中，TAN-1标记的诱导显著降低(图5 e)，证明BHLHE40在中性粒细胞向TAN-1表型的极化过程中起着关键作用。
接下来，我们试图探索BHLHE40对中性粒细胞促肿瘤功能的影响。有趣的是，我们观察到与过表达bhlhe40的dHL-60共培养的PDAC细胞增殖和迁移能力显著增强(图6 a, B，在线补充图S9D-F)，表明上游调控因子BHLHE40驱动中性粒细胞向肿瘤前亚型发展。此外，免疫学评估显示，bhlhe40过表达的dHL-60细胞对CD8的促炎细胞因子产生有抑制作用+T细胞及淋巴细胞增殖能力(图6 c, D)，提示BHLHE40可能是骨髓源性抑制细胞(MDSCs)的潜在调节因子。核bhlhe40阳性tan经IF染色鉴定(在线补充图S9G)， IHC分析显示BHLHE40共定位+PDAC微环境中表达TAN-1标记基因VEGFA的中性粒细胞和中性粒细胞(图6 e)．值得注意的是，BHLHE40浸润水平较高+中性粒细胞与不良预后相关(图6 f)，进一步证明了TANs中上游调控因子BHLHE40的促肿瘤作用。
BHLHE40驱动中性粒细胞向肿瘤前表型发展。(A)菌落形成实验，对照组和bhlhe40过表达的dHL-60细胞培养于上室，PATU-8988细胞培养于下室。菌落数用SD表示为均值。(B)将PATU-8988细胞与对照细胞和过表达bhlhe40的dHL-60细胞共培养3天，并用transwell法检测其迁移能力。迁移的细胞数用SD表示为平均值。(C) PBMCs与对照组和bhlhe40过表达的dHL-60细胞共培养3天，PMA/离子霉素刺激，IFNγ的百分比+和肿瘤坏死因子α+受刺激的CD8细胞+T细胞用流式细胞术分析。来自四个不同捐赠者的数据汇总在右边的面板中。(D) PBMCs用CellTrace Violet染色，在无dHL-60细胞(对照或bhlhe40过表达)和CD3增殖的情况下，用抗人CD3/CD28刺激4天+流式细胞术分析淋巴细胞。三倍培养物的增殖指数用均值和SD表示。(E)连续载片上代表性PDAC组织的免疫组化图像，分别为CD66b(左)、BHLHE40(中)和VEGFA(右)染色。根据CD66b染色和多核形态鉴定中性粒细胞。粉色箭头显示表达BHLHE40或VEGFA的中性粒细胞。(F)免疫组化分析中显示PDAC患者总体生存期的Kaplan-Meier生存曲线。根据BHLHE40的百分比将患者平均分为两组+PDAC组织中总细胞的TANs。* P < 0.05;* * p < 0.01;* * * p < 0.001。dHL-60，中性粒细胞样分化HL-60;HPF，高功率场;干扰素干扰素;包含IHC,免疫组织化学;数控、消极控制;OE,超表达; PBMCs, peripheral blood mononuclear cells; PDAC, pancreatic ductal adenocarcinoma; TAN, tumour-associated neutrophils; TNF, tumour necrosis factor.
讨论
癌症在复杂的微环境中发展，在这种微环境中，肿瘤相关基质被劫持，创造一个有利的生态位，支持肿瘤的发展。3.另一方面，考虑到基质细胞的功能可塑性，它们也有能力在再教育中抑制肿瘤生长，56这为开发针对恶性肿瘤的新治疗策略提供了见解。在PDAC基因工程小鼠模型中，抑制肿瘤微环境中的中性粒细胞浸润导致肿瘤消退和延长生存时间，15 57证明了TANs在支持PDAC进展中的作用，提示TANs可能是一个潜在的治疗靶点。
先前的单细胞RNA-seq研究显示了PDAC微环境中主要细胞类型的综合基因表达图谱。月16日然而，尽管PDAC中存在大量TAN浸润并与不良预后相关，但中性粒细胞的特征研究较少。在本研究中，我们的目的是在单细胞水平上探索PDAC肿瘤微环境中中性粒细胞的异质性，以确定促进肿瘤进展的特定亚群，并研究驱动中性粒细胞向肿瘤前表型发展的潜在机制。为了描述肿瘤微环境中中性粒细胞的重编程，基于肿瘤微环境中的中性粒细胞在趋化因子的作用下不断从循环中招募这一事实，我们也将外周血中的中性粒细胞作为对照。9中性粒细胞浸润在非癌性组织中少见。58
因为人体中性粒细胞的半衰期很短在活的有机体内(8小时)，非常脆弱在体外，22在肿瘤微环境中很难捕获和破译它们的重编程。值得注意的是，在之前使用10x Genomics Chromium平台进行的大多数人类PDAC单细胞研究中都没有检测到中性粒细胞，月16日尽管只有一个例外，那就是斯蒂尔等成功分析了PDAC组织中表达典型标记基因的粒细胞。59另外，我们用BD Rhapsody成功获得了21 972个中性粒细胞的表达谱。同样，在一项前列腺癌的研究中，BD Rhapsody测序的样本中检测到粒细胞，而10x Genomics没有获得粒细胞。24因此，直接比较这两种单细胞测序平台捕获中性粒细胞的能力是有意义的。
中性粒细胞寿命短，易受伤害，这也给研究它们的功能带来了挑战在体外．在这里，我们利用HL-60细胞系，急性早幼粒细胞白血病细胞，当被二甲亚砜或其他几种刺激处理时，可以分化成中性粒细胞样状态。60 61这使我们能够研究中性粒细胞的促肿瘤和免疫抑制功能，并调查特定因素如何影响它们的表型和功能状态。尽管dHL-60细胞不能完全替代中性粒细胞，但我们相信该模型的实验为中性粒细胞生物学提供了新的见解。
我们的数据为中性粒细胞在PDAC微环境中过渡状态的连续性提供了证据:中性粒细胞经内皮迁移从循环进入PDAC微环境，表现为过渡状态(TAN-3)，然后转向中间状态(TAN-0)，这代表了肿瘤微环境中绝大多数中性粒细胞，没有明显的特征。一小部分中性粒细胞被IFN信号激活并获得n1样表型(TAN-4)。62年9通过肿瘤微环境的教育，中性粒细胞逐渐转化为炎症状态(TAN-2)，导致癌症相关炎症，并最终分化为TAN-1，优先表达肿瘤前分子的“致病亚群”，与PDAC患者的预后不良相关(图7)．Zilionis等破译了肺癌中性粒细胞的多样性，并报道了三个表达模块:(i)表达典型中性粒细胞标记物(S100A8和S100A9)的中性粒细胞，与我们研究中的TAN-3相似;(ii)表达炎症细胞因子和肿瘤前基因的中性粒细胞，对应于TAN-1和TAN-2;(iii)与TAN-4类似，I型ifn反应基因表达强烈的中性粒细胞。63这两个单细胞数据集的比较表明，不同癌症类型的TAN群体结构相似。
本研究主要发现的图表摘要。呃,内质网;PDAC，胰导管腺癌;中性粒细胞,多形核白细胞;棕褐色,肿瘤相关中性粒细胞。
代谢状态影响免疫细胞的表型和极化已被广泛接受。53 54在本研究中，肿瘤前TAN-1亚群表现出高度激活的糖酵解活性(图4一)，糖酵解开关与中性粒细胞中的促肿瘤功能相关(在线补充图S8)．我们的研究结果与“反向Warburg效应”模型一致，即非恶性基质细胞中的糖酵解会产生富含能量的代谢物，如乳酸和丙酮酸，这些代谢物可以被癌细胞作为替代的TCA循环底物，促进能量的产生和肿瘤的生长。64 65糖酵解酶ldha相关乳酸除了能促进恶性细胞增殖外，还能取消肿瘤免疫监视。55 66 67一致地，LDHA在dHL-60细胞中的过表达导致t细胞增殖和激活的抑制(在线补充图S8)．总的来说，糖酵解可能是逆转PDAC中中性粒细胞的促肿瘤作用和免疫抑制活性的潜在治疗靶点。
BHLHE40已成为免疫的关键调节因子。据报道，BHLHE40调节T细胞中细胞因子的产生，以及巨噬细胞的增殖在稳态和2型免疫中。68 69然而，关于它在中性粒细胞中的作用还知之甚少。在我们的研究中，BHLHE40被确定为缺氧和内质网应激的下游靶点，这两种肿瘤微环境中的强大刺激物驱动中性粒细胞向肿瘤前TAN-1表型转变。重要的是，我们的数据表明BHLHE40是促肿瘤中性粒细胞的关键调节因子(图6)．PMN-MDSCs被定义为具有免疫抑制功能的病理激活中性粒细胞。70 71在肿瘤微环境中，缺氧诱导的HIF1α激活和内质网应激被认为是将中性粒细胞转化为免疫抑制的PMN-MDSCs的重要调节因子。71 - 73有趣的是，我们的数据也证明了BHLHE40在调节中性粒细胞免疫抑制功能(图6)，提示BHLHE40可能是PMN-MDSCs的潜在调控因子。综上所述，我们的研究揭示了BHLHE40作为PDAC中TANs治疗靶点的潜力。
总之，我们在PDAC肿瘤微环境中发现了一个嗜中性粒细胞的肿瘤前亚簇，揭示了高糖酵解活性和TANs的肿瘤前功能之间的联系，并证明了缺氧诱导和内质网应激诱导的BHLHE40激活是驱动TANs向肿瘤前表型发展的潜在机制。
数据可用性声明
数据可在一个公共的、开放访问的存储库中获得。本研究产生的单细胞和大块RNA测序数据保存在国家组学数据百科(NODE)，登录代码为OEP003254。
伦理语句
病人同意发表
伦理批准
本研究由上海交通大学医学院附属瑞金医院伦理委员会批准(文献编号2021-161)。参与者在参与研究前获得了知情同意。
致谢
我们感谢NovelBio对单细胞测序的支持。感谢上海免疫学研究所的苏冰博士和王静博士对我们的稿件进行了批判性的阅读，并提出了有见地的建议。
参考文献
补充材料
脚注
LW、YL、YD和XT贡献相当。
贡献者LW, LJ和BS设计了研究并撰写了手稿。LW, YL, XT, MS, JQ, MX, YC, JL和PL进行了实验。LW、YD和JH进行生物信息学分析。TY参与蛋白质组学和代谢组学分析。CWa、TW和LD参与病理工作。CWe、JZ、ZX、HC参与患者招募和样本采集。FB, XD, CP和LJ帮助优化研究和校对论文。SC和BS监督了研究并修改了手稿。所有作者阅读并批准了最终稿件。BS, SC和LJ作为担保人对整体内容负责。
资金国家自然科学基金项目(No. 81871906、82073326 - BS, 81802316 - LJ, 82002460 - LW);中国博士后科学基金项目(No. 2019M661552 - LW);上海交通大学上海市基础研究试点项目(No. 21TQ1400205 - LJ)资助。
相互竞争的利益没有宣布。
患者和公众的参与患者和/或公众未参与本研究的设计、实施、报告或传播计划。
来源和同行评审不是委托;外部同行评议。
补充材料本内容由作者提供。它没有经过BMJ出版集团有限公司(BMJ)的审查，也可能没有经过同行评审。讨论的任何意见或建议仅仅是那些作者(s)和不被BMJ认可。BMJ放弃从放在内容上的任何依赖产生的所有责任和责任。如果内容包含任何翻译材料，BMJ不保证翻译的准确性和可靠性(包括但不限于当地法规、临床指南、术语、药品名称和药物剂量)，并且不对翻译和改编或其他原因引起的任何错误和/或遗漏负责。