《Food and Agricultural Immunology》：Bioinformatics prioritisation of butyrate-associated targets and a four-gene immune-barrier score using public atopic dermatitis transcriptomic cohorts
编辑推荐：
研究人员将丁酸盐相关基因与GSE32924特应性皮炎（atopic dermatitis，AD）转录组整合，鉴定出一个探索性四基因转录组特征（TP53、JAK2、H2AX、ALDH1A1），在GSE16161、GSE120721和GSE121212队列中的样本
研究人员将丁酸盐相关基因与GSE32924特应性皮炎（atopic dermatitis，AD）转录组整合，鉴定出一个探索性四基因转录组特征（TP53、JAK2、H2AX、ALDH1A1），在GSE16161、GSE120721和GSE121212队列中的样本水平AUC（area under the curve，曲线下面积）分别为0.951、0.888和0.900。这些发现并未确立临床可转化性或疾病特异性诊断性能。免疫分析显示出一致的T细胞和中性粒细胞信号。在GSE147424中，基因分布于皮肤各区室，但样本水平汇总未显示清晰分离；供者水平映射不可用。分子对接和100 ns分子动力学（molecular dynamics，MD）轨迹仅为探索性分析，并未确立HDAC3结合。由于没有公共队列包含膳食、微生物组、粪便、组织丁酸盐或干预数据，这些产生假说的发现需通过供者水平、生化和功能学研究进行独立验证。
论文解读
一、研究背景与科学问题
特应性皮炎（atopic dermatitis，AD）是一种慢性复发性炎症性皮肤病，以瘙痒、湿疹样皮损为特征，在不同年龄组和分子内型间表现出显著异质性。当代模型将上皮屏障损伤、免疫激活和环境暴露置于相互强化的网络中，而非简单的线性级联。高通量皮肤分析已使这种异质性可被测量，代表性胶带剥离、单细胞和空间转录组研究显示，皮损AD包含协调的免疫、角质形成细胞和基质程序。
在食品与营养研究领域，肠-皮肤轴提供了实用切入点：饮食塑造肠道微生物代谢，微生物代谢产物可影响免疫和上皮屏障生理。丁酸盐（butyrate）尤为重要，因为它连接可发酵膳食纤维、后生元概念和HDAC（组蛋白去乙酰化酶）相关的免疫调节，包括促进FoxP3+调节性T细胞分化。然而，尚不明确的问题是：由正式丁酸盐相关靶标来源支持的基因是否与可重复的AD皮损皮肤程序重叠，并能否产生在独立队列中验证的固定转录组特征。由于公共队列缺乏膳食、微生物组、丁酸盐暴露和一致的临床元数据，研究仅限于靶标和通路优先排序，而非因果或疗效推断。
二、研究设计与主要结论
研究人员开展了一项多队列生物信息学研究，使用公共GEO转录组数据集。GSE32924用于微阵列发现，GSE16161和GSE120721用于独立微阵列验证，GSE121212作为补充RNA-seq平台验证队列，GSE147424作为公共单细胞细胞情境资源。主要分析包括正式靶标整合、仅发现队列的评分开发和锁定评分在独立批量队列中的验证。蛋白互作、富集和批量免疫分析为次要解释性分析。
研究发现，将正式丁酸盐相关靶标证据与AD皮损转录组整合，将3,991个发现差异表达基因（differentially expressed genes，DEGs）缩减为74个靶标支持的候选基因，并产生由TP53、JAK2、H2AX和ALDH1A1组成的固定符号四基因评分。该评分在两个独立微阵列队列和补充RNA-seq队列中显示出队列水平区分能力。STRING网络密度分析显示候选网络比度匹配背景更密集。富集分析支持免疫-屏障关联。免疫分析显示一致的T细胞和中性粒细胞相关信号。该研究确立了证据层级：转录组优先排序和验证是主要证据，网络、富集和免疫分析为次要解释，单细胞分析和HDAC3对接及分子动力学模拟为探索性、产生假说的分析。
三、关键技术与方法
研究人员整合了CTD（Comparative Toxicogenomics Database，比较毒理基因组学数据库）、SEA（Similarity Ensemble Approach，相似性集成方法）和SwissTargetPrediction三个正式丁酸盐靶标来源，与GSE32924发现队列的limma差异表达分析结果取交集。STRING v12.0蛋白互作网络分析结合MCODE模块检测，采用1,000次度匹配随机置换检验。LASSO逻辑回归和随机森林仅在发现队列中进行特征选择，锁定四基因评分后在独立队列验证。GSE147424单细胞分析使用作者分发的17个处理矩阵。HDAC3-丁酸盐分子对接使用AutoDock Vina，辅以单条100 ns GROMACS分子动力学轨迹和MM-PBSA（Molecular Mechanics Poisson-Boltzmann Surface Area）端点分析。
四、研究结果
正式靶标整合鉴定出丁酸盐相关、靶标优先排序的AD候选基因。 发现差异表达分析鉴定出3,991个调控基因。CTD贡献229个独特基因，SEA贡献38个，SwissTargetPrediction贡献89个，并集含337个独特丁酸盐相关基因。与发现DEGs交集产生74个靶标优先排序的候选基因。四个最终评分基因均仅由CTD层支持。
候选基因形成比度和映射匹配STRING背景更密集的网络。 STRING候选网络含261条诱导边。1,000次探索性DEG大小匹配随机集平均19.24条边；精确映射状态匹配和半八度全局度匹配将零均值提高至146.55条边，观察值仍更高（经验P=0.000999），表明枢纽相关拓扑解释了大部分未调整对比，残余差异为候选优先排序提供拓扑支持。
富集分析支持免疫-屏障关联。 聚焦富集突出NFE2L2调控的炎症基因、免疫系统细胞因子信号、细胞-细胞连接组织、白细胞激活和RUNX1调控的紧密连接组分。方向性GSEA显示AD中免疫应答程序富集，而表皮发育、表皮分化、角化和细胞连接过程负富集。
仅发现队列的机器学习优先排序出四基因评分。 LASSO和随机森林均支持TP53、JAK2、H2AX和ALDH1A1。锁定公式为Score=[z(H2AX) – z(TP53) – z(JAK2) – z(ALDH1A1)]/4。1,000次bootstrap重复全部完成，LASSO-随机森林交集频率为H2AX 93.2%、ALDH1A1 55.0%、JAK2 53.3%、TP53 44.4%。
固定符号四基因评分在独立微阵列队列中显示队列水平区分。 GSE16161 AUC为0.951，GSE120721为0.888。标签置换检验支持与随机分组分配的分离。GSE120721区室特异性分析中，全皮肤AUC为0.733，表皮为0.980，真皮为0.933，全皮肤置信区间包含0.5。
补充RNA-seq验证降低了微阵列平台局限性。 GSE121212中评分AUC为0.9。H2AX、JAK2和ALDH1A1重现发现方向，TP53未显示预期下降。留一基因分析显示省略TP53保留相似AUC（0.891），省略H2AX产生较大下降（0.794）。探索性疾病特异性分析显示皮损银屑病评分高于皮损AD。
多种免疫评分方法观察到一致的T细胞和中性粒细胞相关信号。 ssGSEA显示可重复的中性粒细胞相关增加，T细胞评分在多个队列中方向一致。xCell对Treg、Th2、巨噬细胞、树突状细胞和角质形成细胞/屏障特征的估计一致性较低。
评分基因与免疫-屏障通路评分相关联。 JAK2与皮肤屏障和紧密连接评分相关，H2AX与Th2炎症和T细胞/中性粒细胞相关评分相关，ALDH1A1与T细胞、中性粒细胞相关和Treg相关信号呈负相关。重叠审计确认四个评分基因未出现在自定义免疫、通路或标记参考集中。
探索性GSE147424分析提供描述性细胞情境证据。 评分基因表达分布于角质形成细胞、基质和免疫区室。样本文件水平评分比较未显示清晰皮损-健康分离，因无法恢复样本到供者映射，未作供者水平推断。
探索性HDAC3对接显示低分口袋容纳。 最高排名Vina评分为?4.2 kcal/mol，不确立结合。受体为孤立无金属HDAC3链，缺乏辅阻遏物、肌醇四磷酸和催化Zn2?，无对照配体或重复运行。
探索性单轨迹分析描述轨迹水平模式。 单条100 ns轨迹中蛋白RMSD均值随25 ns块递增，最终20 ns斜率为负。100帧PB端点总数为?14.29 kcal/mol，因省略熵且主要项相互对立，未解释为绝对亲和力估计。
五、讨论与结论
讨论部分强调，四基因评分是数据驱动的转录组优先排序结果，其符号编码相对表达方向而非基因功能或通路激活。?z(JAK2)表示较低JAK2转录本丰度提高了评分，不定义低JAK2为致病性、保护性或临床可转化的风险标志。TP53的负号反映发现表达方向而非p53蛋白活性降低，其RNA-seq方向与微阵列方向不一致，被处理为跨平台敏感评分组分而非稳定机制标志。免疫分析最清晰地收敛于T细胞和中性粒细胞相关激活。结构分析仅提供有限的、产生假说的HDAC3情境。与Trompette等人的小鼠模型比较是方向性的而非机制性的，因为公共人类队列不含饮食、微生物组或丁酸盐测量。
研究结论为：本研究在公共AD队列中鉴定出一个探索性四基因、丁酸盐相关的转录组特征。其临床可转化性和疾病特异性诊断性能尚未确立，发现需独立队列和实验验证。研究未证明丁酸盐暴露、四个基因的直接调控、HDAC3结合、微生物群介导的因果肠-皮肤通路，或膳食纤维、功能食品、后生元的功效。
