《Nature Communications》：Differential methanogenic archaea-induced TLR8-dependent signaling engages gene classes regulated by NF-κB p65 and STAT1/2
编辑推荐：
古菌在人类肠道中普遍存在，但缺乏经典的细菌细胞壁相关免疫配体。对人外周血单个核细胞（PBMCs）进行时间分辨RNA测序（RNA-seq）显示，Methanosphaera stadtmanae和Methanobrevibacter smithii诱导了共享的免
古菌在人类肠道中普遍存在，但缺乏经典的细菌细胞壁相关免疫配体。对人外周血单个核细胞（PBMCs）进行时间分辨RNA测序（RNA-seq）显示，Methanosphaera stadtmanae和Methanobrevibacter smithii诱导了共享的免疫程序，但具有不同的动力学特征：M. stadtmanae触发了更早且更强的反应，而M. smithii的作用则更为渐进。保守的早期上调基因分为Class I基因（优先由M. stadtmanae诱导）和Class II基因（由两种古菌相似诱导）。M. stadtmanae表现出更早且更高效的摄取，而对两种古菌的反应均广泛依赖于TLR8。染色质免疫沉淀-定量PCR（ChIP-qPCR）将M. stadtmanae诱导的更强的Class I基因激活与Class I启动子上更强的NF-κB p65结合相关联，而Class II启动子上的STAT1/2结合则在两种刺激下均发生。以古菌RNA输入进行的剂量-反应分析支持不同的激活阈值，Class II基因在较低输入量下被激活，而Class I基因则需要更高的输入量。综上所述，这些发现支持一种输入量调谐的TLR8依赖性机制，该机制塑造了针对古菌的差异性免疫程序。
论文解读
一、研究背景与科学问题
人类肠道中栖息着多种微生物，其中古菌（Archaea）作为一类独特的微生物群体，长期未受到充分关注。传统观点认为古菌主要存在于极端环境，如热液喷口和高盐湖泊，但近年来的研究逐渐认识到，古菌也广泛分布于温带生态系统，并且是人类皮肤、口腔、呼吸道和消化道微生物群落的稳定成员。在人类肠道中，产甲烷古菌（methanogenic archaea）通过消耗细菌发酵产生的CO?和氢气，在厌氧条件下生长繁殖，参与交叉互养（cross-feeding）相互作用并维持微生态稳定性。尽管在多样性上占比较低，产甲烷古菌可占肠道微生物组的4%，具有重要的生态学意义。
人类肠道中的古菌以甲烷杆菌目（Methanobacteriales）为主，主要包括Methanobrevibacter smithii、Methanobrevibacter intestini和Methanosphaera stadtmanae。这些产甲烷古菌与多种健康和疾病状态相关联：它们与整体健康和长寿相关，但Methanobrevibacter的局部过度丰度与牙周炎、便秘、帕金森病和结直肠癌相关，而在炎症性肠病（IBD）中则观察到其耗竭。M. stadtmanae在IBD中表现出促炎特性，但其早期定植与儿童哮喘风险降低相关。这些看似矛盾的关联表明产甲烷古菌的作用具有背景依赖性。
然而，古菌如何与宿主免疫系统相互作用，其机制仍知之甚少。与细菌不同，古菌缺乏经典的细胞壁相关免疫刺激分子，如脂多糖（LPS）或典型肽聚糖。此前的研究表明，M. smithii和M. stadtmanae可激活人外周血单个核细胞（PBMCs）和单核细胞来源树突状细胞（moDCs），诱导肿瘤坏死因子（TNF）和白细胞介素-1β（IL-1β）分泌。进一步研究揭示，古菌的免疫识别主要由RNA驱动，M. stadtmanae的细胞因子诱导严格依赖吞噬作用和内体酸化，TLR8的敲除或药理学抑制几乎完全消除了反应，确立了古菌RNA作为先天感知的主要触发因子。
尽管已有这些进展，大多数研究仅检测了少数细胞因子读数，未涉及更广泛的转录程序和调节回路。特别是，对于两种重要的肠道产甲烷古菌M. stadtmanae和M. smithii的比较分析尤为缺乏。这两种古菌栖息于相同的宿主环境，但在基础生物学上存在差异：M. stadtmanae具有高度受限的产甲烷代谢和相对紧凑的基因组，而M. smithii是优势肠道产甲烷菌，具有适应肠道栖息地的基因组和代谢特征。此外，两种物种在免疫刺激能力上也存在差异，M. stadtmanae通常比M. smithii引发更强的促炎反应。这一研究空白阻碍了对保守与物种特异性免疫通路如何在人类免疫细胞中被编程的机制性理解。
基于上述背景，研究人员开展了一项全基因组转录组分析，利用时间分辨转录组学解析人类PBMCs及PBMC来源单核细胞对两种重要肠道产甲烷古菌反应的动态变化，提出了细胞进入途径和TLR8介导的古菌RNA感知驱动不同转录程序的模型，为古菌-宿主相互作用提供了机制性见解。该研究发表在《Nature Communications》。
二、主要关键技术方法
研究人员从九名独立人类供体（n=9）中分离PBMCs，采用配对设计进行RNA测序（RNA-seq），细胞分别接受M. stadtmanae刺激、M. smithii刺激或不刺激对照，持续刺激4小时或24小时。利用Sleuth统计模型进行差异表达分析，采用KEGG通路富集分析、基因本体论（GO）分析、STRING蛋白-蛋白相互作用网络分析、基因集富集分析（GSEA）等方法。通过19重U-PLEX多重免疫测定法检测细胞因子分泌，使用RT-qPCR验证基因表达，采用共聚焦显微镜和成像流式细胞术定量古菌摄取，利用辅酶F???自发荧光和物种特异性抗体进行检测，使用TLR8选择性拮抗剂CU-CPT9a进行药理学抑制实验，并在TLR8敲除（KO）的BLaER1单核细胞中进行遗传学功能验证，通过染色质免疫沉淀-定量PCR（ChIP-qPCR）检测转录因子结合。
三、研究结果
共享和刺激特异性免疫转录程序在PBMCs中具有不同的时间动态
通过联合统计模型鉴定出10,261个在整体实验设计中具有显著效应的基因。KEGG富集分析显示，前20条通路均以正向富集为主，包括沙门氏菌、耶尔森菌、志贺菌、致病性大肠杆菌等细菌感染通路，以及EB病毒、巨细胞病毒、卡波西肉瘤相关疱疹病毒等病毒感染通路，这些术语主要表示宿主侧信号和效应模块。同时，Toll样受体、NOD样受体和C型凝集素受体等关键先天免疫信号级联显著富集。此外，非酒精性脂肪肝病、脂质代谢和AGE-RAGE信号等与代谢炎症相关的通路也被富集。
在4小时时，M. stadtmanae诱导的差异表达基因（DEGs）数量多于M. smithii。到24小时时，DEG数量趋于收敛。主成分分析显示刺激样本与对照明显分离，4小时时两种古菌的转录组谱分离更明显，24小时时更相似但仍未完全重叠。重叠分析鉴定出443个在所有刺激条件下均上调的保守核心基因，这些基因富集于先天免疫和炎症程序，病毒防御和先天免疫调节主导生物学过程（BP）条目。STRING蛋白相互作用网络分析鉴定出两个主要簇：一个富集干扰素/抗病毒效应因子（包括STAT1、STAT2、ISG15、RSAD2及IFIT和OAS家族），另一个包含促炎细胞因子和趋化因子（如TNF、IL6、CXCL8、CXCL10、CCL2、CCL3），STAT1作为连接两个簇的桥接枢纽。下调模式与上调模式平行，M. stadtmanae在4小时时对单核细胞亚群标志物的抑制更强且更一致，而M. smithii的抑制较为温和，到24小时时两种刺激趋同。
保守443基因核心中Class I/II模块的早期分化
19重免疫测定显示，4小时时M. stadtmanae对大多数靶标诱导的细胞因子分泌明显强于M. smithii。这一趋势在24小时时基本持续，TNF、IL-6和IL-12B在M. stadtmanae组中仍显著较高。然而，IL-8和CXCL10在两种刺激下达相似水平，CCL2在M. smithii刺激后显著更高。基于RNA-seq的log?FC比较，在4小时时将443个核心基因分为两个模块：Class I（247个基因，由M. stadtmanae更强诱导，Δlog?FC≥1）和Class II（195个基因，由两种古菌相似诱导，?1<Δlog?FC<1）。GSEA显示Class I富集于蛋白质转运/细胞外定位和细胞因子活性/受体结合，而Class II缺乏显著富集。RT-qPCR验证了Class I基因在4小时时由M. stadtmanae显著更强上调，而Class II基因在两种刺激间无显著差异。
M. stadtmanae比M. smithii更高效地被吞噬，在单核细胞和moDCs中产生更强的Class I反应
通过辅酶F???自发荧光筛选和物种特异性抗体免疫荧光显微镜检测，研究人员发现CD14?单核细胞是主要的摄取群体，M. stadtmanae在1小时时摄取频率（47.66% vs 7.32%）、2小时（78.85% vs 13.64%）和4小时（83.08% vs 31.04%）均显著高于M. smithii，摄取负荷也一致更高。在纯化的CD14?单核细胞中，Class I基因由M. stadtmanae更强诱导，Class II基因由两种物种相似诱导。细胞松弛素D抑制吞噬作用后，两种古菌诱导的IL-6和IL-1β释放均被强烈抑制，物种间差异消失。将M. smithii输入量提高十倍可达到与标准剂量M. stadtmanae相当的Class I基因诱导水平。在moDCs中，M. stadtmanae诱导的TNF和IL-6显著更高，而CCL2分泌相当，且M. stadtmanae的摄取负荷显著更高。
时间剖析揭示Class I和Class II基因诱导的动力学
在1、2、4、6小时时间点进行的RT-qPCR分析显示，Class I基因在M. stadtmanae刺激下于1至2小时之间急剧上升并持续高水平，而M. smithii诱导的上升较为渐进且水平始终较低。Class II基因的差异主要是动力学性的：1和2小时时M. stadtmanae诱导的转录水平较高，但M. smithii诱导延迟但更陡峭的上升，在4小时时接近M. stadtmanae水平，CCL2在6小时时甚至超过。所有Class II基因在4至6小时之间持续上升，表明4小时并非Class II反应的最大值。
单核细胞对古菌的反应主要为TLR8依赖性，伴有残余TLR8非依赖性特征
两种古菌均诱导LysoTracker阳性单核细胞，M. stadtmanae刺激后比例更高。巴佛洛霉素A1消除LysoTracker信号并显著降低古菌诱导的TNF和IL-6释放，但不影响古菌摄取，支持内体酸化需求。TLR8选择性拮抗剂CU-CPT9a强烈抑制Class I和Class II基因诱导至接近基线水平。RNA-seq分析显示，CU-CPT9a处理后M. smithii特征几乎完全消除（仅剩4个DEGs），而M. stadtmanae保留了235个显著上调基因。其中171个基因在有无抑制剂条件下均显著上调，根据CU-CPT9a敏感性分为64个抑制剂敏感基因（Δlog?FC≥1）和107个基本不受影响的基因（Δlog?FC<1）。在TLR8敲除的BLaER1单核细胞中，IL6（来自CU-CPT9a敏感类别）在TLR8-KO细胞中反应丧失，而SLC7A11（来自不敏感类别）在WT和TLR8-KO细胞中均被诱导，确认了残余的TLR8非依赖性调控。107个CU-CPT9a不敏感基因的STRING网络分析显示三个相连的簇。
转录因子结合和输入阈值塑造Class I/II基因激活
ChIP-qPCR分析显示，M. stadtmanae刺激后在Class I启动子上NF-κB p65结合更强，而Class II启动子上的STAT1/2结合在两种刺激下均发生。剂量-反应分析表明Class II基因在较低RNA输入量下被激活，而Class I基因需要更高输入量，支持输入量调谐的激活阈值机制。
四、讨论与结论
该研究通过时间分辨转录组学系统解析了人类PBMCs和单核细胞对两种重要肠道产甲烷古菌的免疫反应，揭示了共享核心免疫程序与物种特异性动力学特征的结合。M. stadtmanae更高效的吞噬摄取导致更强的早期Class I基因激活，这与NF-κB p65在Class I启动子上更强的结合相关；而Class II基因的STAT1/2结合在两种刺激下均发生，其激活阈值较低。研究提出了一种输入量调谐的TLR8依赖性机制：细胞进入效率和TLR8介导的古菌RNA感知共同塑造了针对不同古菌的差异性免疫程序。这一发现为理解古菌-宿主免疫互作提供了机制性框架，也为TLR8作为治疗调节靶点提供了进一步的理论支持，同时为未来研究输入量调谐的先天信号传导提供了概念参考。
