深入微观战场的全景探查
为了厘清遗传与环境的复杂纠葛，研究团队设计了一项规模宏大且设计巧妙的实验。他们并没有满足于传统的批量测序，而是通过流式细胞术（FACS）分选出了七种主要的免疫细胞类型，并结合了单核甲基化测序（snmC-seq2）和单核染色质开放性测序（snATAC-seq）。
这种技术组合意味着什么？
意味着研究人员可以精确到每一个细胞核，查看DNA上的甲基化修饰（Methylation）——那些像“开关”一样决定基因沉默或活跃的化学标记，同时还能看到染色质（Chromatin）是紧密缠绕还是处于开放状态。
这并不是一个小样本的探索。研究人员从110名个体中收集了171份外周血单个核细胞（PBMC）样本。更关键的是，这些样本涵盖了七种截然不同的环境暴露情况：
病毒感染 艾滋病病毒（HIV-1）、甲型流感病毒（IAV）、新冠病毒（SARS-CoV-2）
细菌与疫苗 耐甲氧西林金黄色葡萄球菌（MRSA）、甲氧西林敏感金黄色葡萄球菌（MSSA）、炭疽疫苗接种
化学暴露 有机磷农药（OPs）
最终，经过严苛的质量控制，研究人员筛选出了104,000个高质量的单细胞数据，并以5kb的基因组窗口为单位进行了聚类分析。这十万个细胞，每一个都成为了一个独立的数据点，共同构建了一个高分辨率的免疫细胞表观基因组图谱。在这个图谱中，我们不仅看到了不同免疫细胞的身份标识，更看到了它们在面对不同环境压力时，那令人惊叹的动态变化。
环境的烙印：增强子上的风暴
在这个研究中，一个最令人着迷的发现是：环境对表观基因组的重塑并非漫无目的，而是具有高度的特异性和靶向性。
为了量化环境的影响，研究人员引入了一个概念：暴露相关差异甲基化区域（exposure-associated Differentially Methylated Regions，简称eDMRs）。你可以把它理解为环境在基因组上留下的“伤痕”或“记号”。
数据显示，在所有暴露类型和细胞类型中，研究人员共识别出了高达756,575个eDMRs。这其中，有517,698个表现为去甲基化（Hypomethylated），即原本关闭的开关被打开了；还有238,877个表现为高甲基化（Hypermethylated），即原本开启的区域被关闭了。平均而言，每一种特定的暴露和细胞类型的组合，都能识别出大约10,000个这样的特异性标记。
这其中，SARS-CoV-2（新冠病毒）、有机磷农药（OP）以及MRSA/MSSA感染对细胞的影响尤为剧烈，它们在大多数细胞类型中都引发了最丰富的eDMRs。这意味着，这些强烈的环境刺激在细胞的DNA层面引发了巨大的震荡。
那么，这些环境烙印究竟打在了哪里？
当研究人员将这些eDMRs与基因组特征进行比对时，一个清晰的规律浮现出来：环境相关的甲基化改变，绝大多数（尤其是去甲基化区域）富集在增强子（Enhancers）区域。
增强子是基因组中一段能够调控基因表达的序列，它们往往距离基因本身很远，但像远程指挥官一样控制着基因的活性。研究发现，eDMRs显著富集在组蛋白修饰H3K4me1和H3K27ac标记的区域，这正是活性增强子的典型特征。
这告诉我们一个深刻的生物学逻辑：环境因素并不直接改变基因的“躯干”（基因体），而是通过改变基因的“开关”（增强子）来快速调节细胞的功能。这是一种极其高效且灵活的策略，允许免疫细胞在保留身份（我是谁）的同时，迅速改变状态（我该做什么）以应对外界危机。
HIV-1：重塑免疫系统的时光胶囊
在所有研究的暴露类型中，HIV-1的数据尤为珍贵，因为它包含了一组极其难得的纵向样本：来自同一供体的感染前（Pre）、诊断当天（Acute）和治疗后200天（Chronic）的样本。这相当于为我们提供了一个观察免疫系统如何随时间演变的“时光胶囊”。
数据揭示了HIV-1感染对免疫细胞的深刻重塑。在流式细胞术和甲基化聚类的分析中，研究人员发现HIV-1感染显著改变了细胞在不同聚类中的比例。特别是自然杀伤细胞（NK cells）、CD8记忆T细胞和CD8初始T细胞，它们的全球甲基化组（Global Methylome）和功能状态被彻底改写。
一个具体的例子发生在CD8记忆T细胞中。研究人员观察到一个位于DGKH基因内含子区域的变化。在HIV-1感染的“急性期”和“慢性期”，这个区域经历了显著的甲基化缺失（Loss of Methylation），与此同时，其染色质开放性（Accessibility）显著增加。
DGKH基因编码二酰甘油激酶，参与细胞内的信号转导。更有趣的是，此前的研究曾报道该基因在HIV精英控制者（能够自发控制病毒的人）和接受抗逆转录病毒治疗的个体之间存在差异甲基化。本研究的数据进一步证实，这种表观遗传的改变是病毒感染直接诱导的，并且在染色质结构上也发生了协同变化。
此外，研究人员还发现，不同阶段的细胞分布极不均匀。通过计算，研究人员识别出了特异性富集在“急性期”的CD8记忆T细胞亚群，这些细胞富集了与“免疫反应正向调节”相关的功能基因。这意味着，在感染初期，免疫系统动员了一支特定的突击队，而它们的动员令就写在表观基因组的甲基化密码中。
新冠风暴眼中的单核细胞
如果说HIV是潜伏的刺客，那么SARS-CoV-2就是正面的强攻。研究人员在分析新冠患者的样本时，发现了一个独特的现象：单核细胞（Monocytes）发生了剧烈的表观遗传分化。
在聚类图上，几乎一半来自重症和非重症新冠患者的单核细胞被分离到了两个独特的亚群中，分别被标记为“Monocyte1”和“Monocyte2”。这两个亚群在对照组样本中几乎不存在，统计学差异极度显著（P值达到了惊人的2.05 × 10⁻²³⁷，这是一个在生物学统计中几乎代表“绝对确定”的数字）。
不仅如此，这两个新冠特有的单核细胞亚群在功能上也大相径庭。通过分析它们之间的差异甲基化基因（DMGs），研究人员发现，在这两个亚群中低甲基化（即被激活）的基因，主要富集在促炎功能相关的通路中，例如“IL-18信号通路”和“白细胞介素-1（IL-1）产生的调节”。
IL-1和IL-18是著名的促炎细胞因子，它们在诱导“细胞因子风暴”中扮演着核心角色。细胞因子风暴正是导致新冠重症患者多器官衰竭的主要原因之一。此外，在“Monocyte2”亚群中，差异甲基化基因还富集了“吞噬作用”和“内吞作用”的功能，暗示这个亚群可能在抗原提呈方面具有特殊功能。
这一发现从表观遗传的角度解释了为什么新冠感染会引发如此强烈的炎症反应：病毒暴露直接打开了单核细胞中促炎基因的甲基化开关，让它们处于一种“亢奋”的战备状态。
基因的锚点：坚守在基因体
当我们把目光从环境转向遗传时，风景截然不同。
为了探究遗传变异如何影响甲基化，研究人员利用一种名为Biscuit的工具，直接从甲基化测序数据中识别出了每个个体的单核苷酸多态性（SNPs）。这是一种非常巧妙的策略，使得他们无需额外的全基因组测序就能获得遗传信息。
通过将遗传变异（Genotype）与甲基化水平（Methylation）进行关联分析（meQTL），研究人员在九种细胞类型中识别出了275,283个基因型相关的差异甲基化区域（genotype-associated DMRs，简称gDMRs）。
与环境留下的印记（eDMRs）不同，遗传留下的印记（gDMRs）在基因组上的分布有着完全不同的偏好。gDMRs主要位于基因体（Gene Body）区域，特别是由H3K36me3这种组蛋白修饰标记的区域。
这是一个极具深意的对比：
eDMRs（环境）： 偏好增强子（Enhancers），这是基因调控的灵敏开关，负责应对瞬息万变的外界信号。
gDMRs（遗传）： 偏好基因体（Gene Bodies），这是基因转录的结构基础，负责维持细胞功能的稳定性。
这种“分工”揭示了生命演化的智慧。遗传变异更倾向于通过影响基因转录的延伸或剪接等基础过程来发挥作用，这种影响往往是持久且稳定的；而环境因素则通过拨动增强子这个“快慢机”，对基因表达进行微调，以适应短期的环境压力。
当基因遇上疾病：解码GWAS的盲区
全基因组关联分析（GWAS）已经为我们找到了成千上万个与疾病相关的遗传位点（SNPs），但其中绝大多数位点并不位于基因的编码区，而是落在所谓的“非编码荒漠”中。这让研究人员经常感到困惑：这个突变到底是怎么导致疾病的？
本研究构建的单细胞表观基因组图谱，为解开这个谜题提供了绝佳的钥匙。研究人员将识别出的meQTL（甲基化数量性状位点）与已知的疾病GWAS位点进行了共定位分析。结果令人振奋：许多与免疫疾病相关的遗传变异，正好与特定细胞类型中的甲基化调控位点重合。
让我们看一个具体的例子：湿疹（Eczema）。湿疹是一种常见的免疫介导皮肤病。GWAS研究发现了一个与之高度相关的SNP位点：rs10791824。然而，光看这个SNP，我们不知道它在什么细胞里起作用，也不知道它影响了哪个基因。
利用本研究的数据，研究人员发现，这个SNP恰好与一个meQTL共定位，而且这个共定位具有极高的细胞特异性——它主要发生在CD8初始T细胞和CD4记忆T细胞中。进一步的SMR（基于摘要数据的孟德尔随机化）分析显示，这个SNP相关的甲基化区域与EFEMP2基因的表达水平显著相关。
这就串起了一条完整的致病逻辑链条：
01. 患者携带了rs10791824这个遗传变异；
02. 这个变异特异性地改变了T细胞中特定区域的DNA甲基化水平；
03. 甲基化的改变进而影响了下游EFEMP2基因的表达；
04. EFEMP2基因表达的异常最终导致了湿疹的易感性。
研究还发现，胆结石病（Gallstone disease）相关的遗传位点也主要富集在T细胞的meQTL中。这提示T细胞的表观遗传调控在这些免疫相关疾病的遗传风险中扮演了核心角色。这种细胞特异性的解析能力是传统的大块组织（Bulk tissue）研究无法企及的。
染色质开放性：DNA甲基化的镜像
为了验证DNA甲基化变化的真实功能后果，研究人员还将目光投向了染色质的物理状态——开放性。利用单核ATAC-seq数据，他们发现DNA甲基化与染色质开放性之间存在显著的负相关关系，尤其是在单核细胞中，相关性最高。
然而，更有趣的是细节中的差异。在HIV-1感染的样本中，研究人员将两种数据进行了整合分析。结果显示，大约25.6%的CD8初始T细胞中的去甲基化区域（hypo-eDMRs）同时也伴随着染色质开放性的增加（ATAC-seq peaks gained）。
这说明，虽然大趋势是一致的，但在很多情况下，DNA甲基化的改变可能早于染色质结构的改变，或者独立于染色质结构发生。这种“不同步”或许暗示了表观遗传调控具有多重保险机制，或者是时间尺度上的滞后效应。
种族差异的表观遗传回响
在探索环境影响时，研究人员还观察到了一个不可忽视的现象：遗传背景（Ancestry）对环境反应的修饰作用。在剔除了年龄和性别的影响后，研究人员发现，在MRSA/MSSA感染和SARS-CoV-2感染中，非洲裔个体（African genetic ancestry）表现出了比欧洲裔个体更强烈的甲基化反应。
具体来说，在MRSA/MSSA暴露组中，研究人员识别出了13,409个非洲裔特有的eDMRs，而欧洲裔特有的仅有5,471个。同样，在新冠暴露组中，非洲裔特有的eDMRs高达16,393个，远超欧洲裔的9,471个。
这些数据与流行病学观察到的现象——即非洲裔美国人在新冠等传染病中往往表现出更严重的疾病后果——不仅在统计上相符，更为其提供了潜在的分子机制解释。这表明，遗传背景不仅决定了基线状态，还可能决定了免疫系统对环境刺激的反应阈值或反应强度。
精密的双重奏
当我们回溯这项长达数年、涵盖十万个细胞、七种环境暴露的宏大研究时，一个清晰的生命图景展现在我们面前。
我们的免疫系统并非刻板的机器，而是一个能够感知环境并做出适应性改变的智能网络。在这个网络中，基因（Genetics）和环境（Environment）并非简单的加减法关系，而是分工明确的合作伙伴：
基因构建框架：
遗传变异主要通过影响基因体（Gene Bodies）的甲基化，特别是富集H3K36me3的区域，来确立细胞功能的基本架构和稳定性。这是生命的“底色”。
环境微调细节：
病原体和化学物质暴露主要通过重塑增强子（Enhancers）的甲基化，特别是富集H3K4me1和H3K27ac的区域，来快速调节基因表达的开关。这是生命的“应变”。
这种机制的发现，不仅回答了“先天与后天”的哲学之问，更为医学提供了极具操作性的启示。
对于临床医生和药物研发人员来说，这意味着如果我们想干预某种环境因素导致的疾病（如病毒感染后的炎症风暴），我们应该瞄准那些特异性的增强子区域；而如果我们想攻克某种遗传性免疫疾病，可能需要深入理解基因体区域的调控机制。
更重要的是，这项研究发现的成千上万个暴露特异性的eDMRs，有望成为未来的生物标志物（Biomarkers）。也许在不久的将来，只需一滴血，通过检测这些特定的甲基化位点，医生就能知道你的免疫系统曾经经历过怎样的战斗——是曾与新冠病毒殊死搏斗，还是正受到某种化学物质的隐秘侵袭。
在这场基因与环境的精密双重奏中，每一个甲基化修饰都是一个音符。通过这项研究，我们终于开始读懂这本写在细胞核里的无字天书，听懂了那首关于生存与适应的生命乐章。
参考文献