关于天然杀伤细胞和CD8+ T细胞免疫记忆表观遗传调控的学术研究报告
第一， 研究作者、机构、发表期刊与时间 本研究由来自纪念斯隆凯特琳癌症中心（Memorial Sloan Kettering Cancer Center）免疫学项目和计算与系统生物学项目的 Colleen M. Lau 和 Nicholas M. Adams 等人，以及威尔康奈尔医学院（Weill Cornell Medical College）免疫学与微生物发病机制系的 Joseph C. Sun 等人合作完成。通讯作者为 Joseph C. Sun。该研究以题为 “Epigenetic control of innate and adaptive immune memory” 的论文形式，于2018年9月发表在 Nature Immunology 期刊上。
第二， 学术背景与研究目的 本研究属于免疫学，特别是天然免疫与适应性免疫记忆的交叉领域。长期以来，克隆性扩增和产生免疫记忆被认为是适应性免疫系统（如T细胞和B细胞）的独有特征。然而，近年的研究颠覆了这一传统观念，证明天然杀伤（Natural Killer， NK）细胞在病毒（如小鼠巨细胞病毒，MCMV）感染过程中，也能表现出类似适应性的特征，即由特异性受体（如Ly49H，识别病毒蛋白m157）介导的克隆扩增、收缩，并形成能够产生保护性回忆应答的长寿命记忆NK细胞。尽管先前的研究已经描述了MCMV感染过程中NK细胞独特的转录组特征，但人们对NK细胞从初始状态到效应状态，再转变为记忆状态的过程中，转录是如何在表观遗传水平被调控的，尚不清楚。
因此，本研究的主要目标是：第一，阐明病毒特异性NK细胞在感染过程中染色质可及性的动态变化，以及这种变化如何与其转录谱相关联；第二，探究关键的炎症信号通路（特别是通过STAT4和STAT1的细胞因子信号）如何影响早期的表观遗传重塑；第三，通过平行分析经典适应性免疫细胞——病毒特异性CD8+ T细胞在相同感染模型中的表观遗传景观，探究NK细胞与T细胞在生成免疫记忆过程中是否共享共同的表观遗传程序。
第三， 详细研究流程与方法 本研究设计了一个系统性的多组学分析流程，核心在于对Ly49H+ NK细胞在MCMV感染时间进程中进行平行的高通量染色质可及性测序（Assay for Transposase-Accessible Chromatin using high-throughput sequencing， ATAC-seq）和转录组测序（RNA-seq）。具体流程如下：
实验模型与样本制备：
- 研究对象： 主要研究对象是C57BL/6背景小鼠中，对MCMV特异性的Ly49H+ NK细胞。作为对比，研究还分离了针对同一病毒MCMV抗原（M45）特异性的CD8+ T细胞。
- 感染与时间点： 小鼠感染MCMV，在感染后第0天（初始，D0）、第2天（D2）、第4天（D4）、第7天（D7）、第14天（D14）和第35天（记忆，D35）获取脾脏细胞。对于CD8+ T细胞，则分析了初始（D0）、效应期（D7）和记忆期（D35）的细胞。
- 细胞分选： 通过流式细胞分选技术，从脾脏中精确分选纯化的Ly49H+ NK细胞（TCRβ- CD19- CD3ε- F4/80- NK1.1+ Ly49H+）以及M45四聚体阳性的CD8+ T细胞。每个时间点通常包含2-4个生物学重复样本。
多组学数据生成：
- ATAC-seq： 对分选的细胞进行ATAC-seq实验，以全基因组尺度检测染色质的开放区域。该方法利用转座酶切割开放的染色质区域并插入测序接头，从而揭示DNA的可及性状态。研究人员对每个样本生成平均约5100万条双端测序读数。
- RNA-seq： 平行地，从相同分选的细胞群中提取RNA，构建cDNA文库进行RNA-seq测序，以获取全基因组的转录表达谱。每个样本平均生成约4700万条双端读数。
- ChIP-seq（染色质免疫沉淀测序）： 为了研究特定转录因子（STAT4和STAT1）的结合情况，研究人员用IL-12 + IL-18（激活STAT4）或IFN-α（激活STAT1）体外刺激纯化的NK细胞，然后分别使用抗STAT4和抗STAT1抗体进行ChIP-seq实验，以确定这些STAT蛋白在全基因组的结合位点。
遗传缺陷模型验证：
- 为了直接验证STAT4和STAT1在早期感染中对表观遗传和转录的影响，研究人员构建了混合骨髓嵌合体小鼠，其中一半骨髓细胞来自野生型小鼠，另一半来自Stat4或Stat1基因敲除（KO）小鼠。这使得在同一个宿主内，可以同时比较野生型和KO型Ly49H+ NK细胞在感染MCMV后（特别是D2）的表观遗传和转录状态。对这些细胞也进行了ATAC-seq和RNA-seq分析。
生物信息学与数据分析工作流：
- 数据比对与标准化： 所有测序读数经过质量修剪后，比对到小鼠参考基因组（mm10）。使用软件如MACS2进行ATAC-seq和ChIP-seq的峰值（peak） calling，并结合Irreproducible Discovery Rate (IDR)分析确保重复样本间的可重复性。最终合并所有样本的峰值区域，构建一个综合的“染色质可及性图谱”用于后续比较。
- 差异分析： 使用DESeq2软件包进行差异分析。对于ATAC-seq，鉴定不同时间点之间差异可及的染色质区域（differentially accessible， DA peaks）；对于RNA-seq，鉴定差异表达基因。设定严格的统计阈值（错误发现率FDR < 0.05）。
- 注释与功能富集： 将DA peaks根据其相对于基因的位置（启动子区、内含子区、基因间区等）进行注释。使用GREAT工具对DA peaks相关的基因集合进行通路富集分析，以揭示生物学过程。
- 整合分析与关联： 对ATAC-seq和RNA-seq数据进行主成分分析和分层聚类，观察其动态模式的异同。计算特定通路中基因的染色质可及性变化与表达变化之间的相关性。
- 转录因子结合模体分析： 使用HOMER软件在差异可及区域中进行从头（de novo）转录因子结合模体搜索，以识别可能调控这些区域开放或关闭的关键转录因子家族。
- NK细胞与CD8+ T细胞比较： 创建一个包含NK细胞和CD8+ T细胞共有染色质可及区域的合并图谱。通过对细胞类型内部归一化后的可及性值进行分析，比较两种细胞在从初始到记忆转变过程中的共同和独特表观遗传特征。
第四， 主要研究结果 本研究通过上述严谨的实验和分析流程，获得了一系列层次分明的关键结果：
NK细胞在病毒感染过程中经历剧烈的染色质动态重塑：
- ATAC-seq时间进程分析显示，病毒特异性NK细胞在整个感染过程中，染色质可及性发生了显著变化。这些变化主要发生在感染第一周（D0到D7），而在后期（D14到D35）变化很小，表明表观遗传状态在记忆细胞形成早期就已趋于稳定。
- 变化最剧烈的区域多位于内含子和基因间区（推测为增强子区域），而启动子区域的变化相对温和。聚类分析揭示了不同“波次”的可及性变化：有些区域在记忆细胞中持续开放或关闭（稳定变化，如与Socs3， Pdcd1等基因相关的区域），而另一些则只在感染早期短暂改变后恢复到基线水平（瞬时变化，如与Tbx21， Ifng等基因相关的区域）。
- 通路富集分析表明，这些动态可及区域关联的基因涉及NK细胞介导的细胞毒性、细胞因子信号（如IL-2、IL-12、JAK-STAT通路）以及T细胞抗原受体信号通路，这共同反映了NK细胞兼具的天然和适应性免疫功能。
早期感染中转录与表观遗传的协调与分化：
- 主成分分析显示，ATAC-seq和RNA-seq的轨迹在感染早期（D0-D7）高度相似，但在后期（D7-D35）分道扬镳。染色质重塑在D14基本完成，而基因表达在D14至D35间仍在持续调整，提示表观遗传层面的“记忆命运决定”可能先于转录组层面的最终稳定。
- 早期时间点（特别是D0-D2和D2-D4）染色质可及性与基因表达呈显著正相关的通路包括I型和II型干扰素（IFN）信号通路、NFAT信号通路等，这些都是抗病毒NK细胞应答的关键通路。例如，Ifng基因座的可及性和表达同步上调。
STAT4和STAT1以不同模式介导染色质重塑：
- ChIP-seq分析发现，IL-12/IL-18诱导的STAT4主要结合在增强子区域（内含子/基因间区），而IFN-α诱导的STAT1主要结合在启动子区域。
- 在感染早期（D2），STAT4结合的区域显示出染色质可及性的强烈瞬时增加，而STAT1结合的区域则表现出温和的、持续的下降趋势。
- Stat4敲除导致NK细胞在早期感染时，非启动子区域可及性显著下降，且与STAT4结合、可及性变化和基因表达下调三者重叠的基因（如Fyn， Ifng）被识别出来，表明STAT4可能直接促进增强子区域开放以驱动转录。相反，Stat1敲除的影响模式不同，且三者重叠的基因很少，提示STAT1的作用可能更间接或依赖于其他信号。
- 交叉调控分析发现，Stat4敲除上调了Stat1和Stat2的表达，增强了I型IFN应答基因；而Stat1敲除则增强了Ifng的表达，表明两条通路相互拮抗，共同微调NK细胞应答。
记忆NK细胞具有独特的“预备就绪”表观遗传状态：
- 与初始NK细胞相比，记忆NK细胞约有30%的染色质可及区域发生了显著改变（log2FC>0.5）。
- 在记忆细胞中可及性更高的区域，富集了NK细胞介导的细胞毒性通路（如Prf1基因座），同时这些基因的表达也更高，表明记忆NK细胞在表观遗传上预备好快速执行杀伤功能。
- 在记忆细胞中可及性降低的区域，富集了MAP激酶信号等通路。
- 模体分析揭示，记忆细胞中开放的区域富含干扰素刺激反应元件样模体；而关闭的区域则富含TCF-LEF（与Wnt信号相关）和NF-κB转录因子家族的结合模体。TCF7和NFKB1等基因在记忆阶段表达下调，与其基因座可及性降低一致，这可能限制了某些信号通路在记忆细胞中的过度激活。
NK细胞与CD8+ T细胞共享免疫记忆的表观遗传特征：
- 通过对比两种细胞在感染过程中的染色质可及性，研究者发现一个有趣的现象：经过细胞类型内归一化后，初始NK细胞的表观遗传状态与记忆CD8+ T细胞最为接近，这为NK细胞能快速启动效应功能提供了表观遗传学证据。
- 进一步分析记忆相对于初始状态的共同变化，发现尽管两种细胞整体的表观遗传景观各异，但仍有21%的差异可及区域表现出相同的趋势（在记忆中都升高或都降低），而只有9%的差异表达基因共享相同趋势。这表明存在大量不立即导致转录改变的共同表观遗传变化，可能代表了“预备”状态。
- 这些共同调控的基因区域关联到一些已知对T细胞记忆至关重要的基因，如Bach2， Tcf7， Zeb2， 以及效应分子（Gzmb， Prf1）和表面标志物（Sell， Klrg1）。
- 模体分析发现，在两类细胞记忆状态中共同变得更开放的区域，富集了AP-1转录因子家族（如Jun蛋白）的结合模体；而共同关闭的区域则再次富集了TCF模体。这表明AP-1和TCF-LEF家族可能共同调控着天然与适应性淋巴细胞记忆的生成。
第五， 研究结论与价值 本研究系统揭示了NK细胞在获得免疫记忆过程中，其染色质景观经历了与CD8+ T细胞类似的、动态且有序的表观遗传重编程。关键结论如下： 1. NK细胞的记忆形成伴随着广泛且稳定的染色质可及性变化，这些变化在感染早期（约两周内）基本完成，使记忆细胞处于一种功能特化的“表观遗传预备”状态。 2. 不同的细胞因子信号（IL-12/STAT4 vs. IFN/STAT1）通过结合于不同的基因组区域（增强子 vs. 启动子），以截然不同的方式塑造着早期的表观遗传和转录程序，并存在交叉调控。 3. 尽管NK细胞与T细胞在受体和发育上不同，但它们在形成长期保护性记忆的过程中，共享了核心的表观遗传调控逻辑和关键转录因子网络（如AP-1和TCF）。
这项研究的科学价值在于：首次在表观基因组层面深入剖析了天然淋巴细胞免疫记忆的分子基础，填补了该领域的知识空白；通过平行比较NK细胞与CD8+ T细胞，揭示了天然与适应性免疫记忆在深层调控机制上的共性与特性，为理解免疫记忆的进化提供了新视角。其应用潜力包括：为通过表观遗传手段干预或增强免疫记忆（例如在疫苗设计或免疫治疗中）提供了新的理论依据和潜在的靶点；发现的共享调控基因和通路（如Bach2， Tcf7， AP-1等）是未来研究免疫记忆形成机制的宝贵线索。
第六， 研究亮点 1. 研究视角新颖： 将表观遗传学分析首次系统性地应用于NK细胞适应性免疫特征的研究，特别是其记忆形成过程，选题具有前沿性和创新性。 2. 实验设计系统全面： 采用了时间进程的多组学（ATAC-seq + RNA-seq）平行分析，并结合了ChIP-seq和遗传缺陷模型验证，构成了一个从观察到机制验证的完整证据链。 3. 比较生物学意义深刻： 创新性地将NK细胞与经典的CD8+ T细胞记忆模型进行头对头的表观遗传学比较，发现了超越细胞类型的保守记忆特征，结论更具普遍性和深度。 4. 发现重要机制： 明确了STAT4在直接重塑增强子可及性中的关键作用，以及记忆NK细胞中特定转录因子模体（ISRE， TCF， AP-1）的富集状态，提出了“表观遗传预备”这一重要概念。 5. 提供宝贵资源： 本研究产生的全面基因组学数据集，为免疫学和表观遗传学领域研究者进一步探索淋巴细胞分化与记忆提供了宝贵的数据资源。
第七， 其他有价值内容 本研究的补充材料（Supplementary Information）包含了详尽的实验方法细节、数据分析的扩展图表、以及STAT缺陷型小鼠的表型验证数据（如细胞扩增缺陷），这些内容为论文结论的可靠性提供了坚实的支撑。此外，论文最后讨论部分不仅总结了发现，还就增强子与启动子调控的差异、TCF因子的潜在作用、STAT信号的复杂性以及未来研究方向（如组蛋白修饰分析、候选基因功能验证）进行了深入的思考和展望，启发了后续研究。