《Genes》：Cross-Organ Transcriptomic Integration Reveals Cross-Disease Transcriptional Concordance and Shared Cell-State Patterns Between Metabolic Dysfunction-Associated Steatotic Liver Disease and Chronic Atrophic Gastritis
编辑推荐：
背景/目的：跨器官整合转录组证据受队列间异质性限制，因此可重复的跨器官信号通常以基因列表而非经过验证的候选基因形式报告。为此，研究人员开发并应用了一个跨器官整合框架，其中发现阶段与裁决阶段相分离。方法：将八个公共批量转录组队列（六个肝脏队列、两个胃队列）通过稳
背景/目的：跨器官整合转录组证据受队列间异质性限制，因此可重复的跨器官信号通常以基因列表而非经过验证的候选基因形式报告。为此，研究人员开发并应用了一个跨器官整合框架，其中发现阶段与裁决阶段相分离。方法：将八个公共批量转录组队列（六个肝脏队列、两个胃队列）通过稳健秩聚合（RRA）进行合并，并将共享候选基因带入预先指定的稳健性审计，该审计结合了留一数据集重算、效应量荟萃分析和细胞组成校正。候选基因集在发现后不再重新选择，保留的候选基因在两个独立的、非发现阶段的外部队列中进行检验。结果：应用于代谢功能障碍相关脂肪性肝病（MASLD）和慢性萎缩性胃炎（CAG）时，该框架返回四个基因：CDHR2和ANXA4定义了上皮模式，IL32定义了上皮界面的免疫相关模式，CADM2定义了肝星状细胞模式；LGALS3作为组成敏感性对照单独保留。每个保留的候选基因在两个保留队列中均保持其发现阶段的方向。在胃复制队列中（该检验跨十个混合文库运行），CDHR2和CADM2在Benjamini–Hochberg校正后仍显著（两者q = 0.020）；对于其他三个候选基因，胃部对比仅为方向性。在单细胞和单核数据中，候选基因解析为三种细胞状态模式，与独立的蛋白水平注释相匹配，探索性孟德尔随机化未返回一致的、FDR支持的遗传证据。结论：五个候选基因中有三个（CDHR2、CADM2和IL32）满足较高证据层级的所有三项预先指定标准，两个（ANXA4和LGALS3）满足较低层级。数据所确立的是沿肠-肝轴的跨疾病转录一致性，通过共享细胞状态定位来解释，而非共享病因、因果关系或共同生物标志物。该审计框架可直接迁移至其他跨器官比较，五个基因作为表达支持的候选基因报告，准备进行实验随访。
研究背景方面，代谢功能障碍相关脂肪性肝病（MASLD）曾被称为非酒精性脂肪性肝病，目前估计影响全球38%的成年人，预计到2040年患病率将超过55%。其组织学从单纯性脂肪变性经脂肪性肝炎进展为进行性纤维化、肝硬化和肝细胞癌。在所有组织学特征中，纤维化分期最能预测该人群的肝脏相关和全因死亡率。治疗选择开始扩大，但用于常规活动度或进展分级的分子特征仍然缺乏。慢性萎缩性胃炎（CAG）在胃癌发生中占据类似位置。在Correa模型中，长期黏膜炎症（最常由幽门螺杆菌驱动）从萎缩性胃炎经肠化生进展为异型增生和腺癌。人群队列显示萎缩性胃炎的年胃癌风险约为0.1%，约为正常黏膜的三倍。根除幽门螺杆菌可降低胃癌发病率，但仅部分降低（合并发病率比0.53，95%置信区间0.44–0.64）。即便如此，2020年胃癌在全球造成超过75万人死亡。一些高发地区使用血清学筛查，但组织学仍是风险评估的参考标准。这两种疾病倾向于被孤立研究，但它们在多个层面交叉。一项纳入34项研究和218,573名参与者的荟萃分析将幽门螺杆菌感染与MASLD联系起来（合并比值比1.26，95%置信区间1.17–1.36）。肠-肝轴是明显的解剖和功能通道。当胃和肠道功能障碍时，营养物质、胆汁酸和微生物产物向肝脏的门静脉输送发生变化。一个器官的慢性炎症可以维持全身低度炎症，反过来影响另一个器官。这两种疾病还汇聚于相似的组织水平反应，即上皮损伤、免疫细胞浸润和成纤维细胞活化，这些反应驱动两个器官的进行性重塑。这种汇聚是否反映肝和胃真正共享的分子模式迄今尚未被检验。用转录组学回答这个问题并不简单。表达研究通常在一个队列中分析一个器官；小样本、平台异质性和组织细胞组成的变化使基因列表难以重现。基于秩的聚合提供了一种跨平台整合差异表达结果而不合并原始数据的方法。随后可以通过三项互补检查来 interrogate 批量组织信号：在保留队列中的复制、在单细胞数据中定位到特定细胞类型，以及使用大规模表达数量性状位点资源的遗传工具进行因果检验。联合应用于MASLD和CAG，该框架应将可重复的跨器官信号与队列特异性噪声分开。研究人员将这一框架应用于八个批量转录组队列（六个肝脏、两个胃）。稳健秩聚合（RRA）随后识别出在MASLD和CAG中同向变化的基因。然后研究人员将14个共享候选基因置于预先制定的稳健性审计中，结合留一数据集分析、异质性荟萃分析和细胞组成校正。五个领先基因（IL32、CDHR2、ANXA4、CADM2和LGALS3）和一个观察列表基因（RPS6KA1）从审计中产生；另外八个候选基因被降级。随后研究人员在外部数据中检验这五个基因，不再重新选择。所有五个基因均保持其发现阶段方向，在216例活检肝脏系列中显示与NAFLD活动评分（NAS）的梯度，其中三个在外部胃队列中达到显著性。在两个器官的单细胞和单核数据集中，领先基因分为三种组织或细胞状态模式：上皮模式（CDHR2、ANXA4）、上皮界面的免疫相关模式（IL32）和星状细胞相关模式（CADM2）。这些定位与人体组织中的细胞类型免疫组化相匹配。探索性孟德尔随机化（MR）分析在可用分析中未返回一致的FDR支持的遗传证据，因此研究人员将这些基因框定为优先候选基因而非确认的因果驱动因素。据研究人员所知，这是MASLD和CAG的首个跨器官转录组综合，产生了用于机制和转化随访的优先候选基因集。在研究方法上，研究人员使用了以下关键技术方法：八个公共批量转录组队列的差异表达分析（使用limma-voom处理计数数据和limma处理归一化表达矩阵）；器官特异性稳健秩聚合（RRA）整合队列水平基因排名；稳健性审计包括留一数据集（LODO）精确RRA重算（闭式p值，N=20,000基因背景）、随机效应荟萃分析（I2统计量）和替代细胞组成校正；两个独立外部队列验证（GSE135251肝脏队列n=216；GSE153224胃队列40例患者10个混合文库）；三个公共单细胞数据集的患者水平伪批量分析（GSE202379肝脏单核RNA测序47例患者99,687细胞；GSE115469肝脏单细胞RNA测序5例移植供体仅用于细胞类型确认；GSE134520胃窦黏膜单细胞RNA测序9例患者13例活检）；Human Protein Atlas（版本25.1）免疫组化细胞类型注释和抗体可靠性评分；MuSiC版本1.0.0标记限制性反卷积敏感性分析（四个批量肝脏RNA测序队列：GSE126848、GSE130970、GSE162694和GSE174478）；探索性两样本孟德尔随机化使用eQTLGen血液顺式表达数量性状位点工具（全基因组显著性p < 5 × 10?8，r2 < 0.001对1000 Genomes EUR参考面板clumping）和七个肝脏及胃相关结局全基因组关联研究（包括肝脏MRI衍生质子密度脂肪分数、NAFLD、肝纤维化/肝硬化、胃炎和胃癌），产生69个估计。所有分析在R版本4.6.0中运行。
研究结果方面，3.1 跨队列筛选和稳健性审计：发现集包括六个肝脏队列（Ahrens GSE48452；Arendt GSE89632；Lefebvre GSE83452；Suppli GSE126848；Hoang GSE130970；Pantano GSE162694）跨越从脂肪变性到晚期MASLD纤维化，以及两个胃队列比较CAG与非萎缩胃黏膜（GSE116312 CAG对比滤泡性胃炎；GSE27411萎缩对比非感染胃）。器官特异性稳健秩聚合返回14个在两个器官中一致排名靠前的基因。由于基于秩的整合可能被单个队列主导，所有14个候选基因经过审计，包括留一数据集精确RRA重算（闭式精确p值对N=20,000基因背景）、随机效应荟萃分析（I2）的每队列效应量估计、标记衍生组成评分的替代细胞组成校正和批次诊断。审计保留五个领先基因（IL32、CDHR2、ANXA4、CADM2和LGALS3）和一个观察列表基因（RPS6KA1）；其他八个基因（GOLM1、ANO10、SLC6A16、KIAA1958、FGA、FGB、LEPR和TSPAN3）被降级，其信号在队列间不一致或由单个队列驱动。保留基因之一LGALS3在单个发现队列（Arendt；log2倍数变化?0.24，p = 0.013）中方向不一致，而其余队列和外部验证支持它；该例外记录在补充表S4中。3.2 保留数据集中的外部复制：五个领先候选基因从发现筛选中推进，研究人员在两个外部数据集中检验它们而不重新选择。原始14个候选基因中，13个在GSE135251（肝脏；n=216含NAS和纤维化分级）中重现其发现阶段变化方向。对于五个领先基因，病例对照平均差异（log2 CPM）为IL32 +1.26（p = 1.5 × 10?3）、CDHR2 +1.43（p = 1.2 × 10?5）、ANXA4 +0.13（p = 0.48）、CADM2 ?0.50（p = 8.2 × 10?6）和LGALS3 +0.21（p = 0.39），NAS梯度相关性为ρ = 0.48（IL32；95% CI 0.37–0.58）、0.47（LGALS3；0.36–0.57）、0.34（ANXA4；0.22–0.46）、0.25（CDHR2；0.12–0.37）和?0.17（CADM2，p = 0.013；?0.30至?0.04）。作为辅助判别指标，CDHR2的接收者操作特征曲线下面积为0.911（bootstrap 95% CI 0.860–0.956），IL32为0.798（0.684–0.900），这些估计基于小对照组（n=10）。在GSE153224中，一个独立的CAG对比CNAG全转录组数据集覆盖40例患者10个混合文库（每组5个，每池4例患者），所有五个候选基因均按预期方向变化（log2(FPKM + 1)单位）：CDHR2 +2.16（p = 0.008）、CADM2 ?1.65（p = 0.008）、LGALS3 +1.53（p = 0.032）、ANXA4 +0.66（p = 0.15）和IL32 +0.49（p = 0.55）。CADM2作为低丰度基因，可在70%的混合文库中检测到。原始GSE153224研究本身通过其自己的共表达和蛋白相互作用流程将CDHR2优先纳入八个枢纽基因；由于这一汇聚来自同一数据集，研究人员不将其计为独立队列复制。3.3 细胞状态定位定义三种组织或细胞状态模式：三个单细胞数据集的患者水平伪批量分析（GSE202379，47例患者；GSE115469仅用于细胞类型分配确认；GSE134520，9例患者）将五个候选基因分为三种组织或细胞状态模式。第一种模式是上皮性。CDHR2定位到肝细胞和胃黏膜上皮，具有近乎排他的谱；ANXA4在两个肝脏数据集中97–100%的胆管细胞中检测到，在74%的胃上皮细胞中检测到。第二种模式位于上皮-免疫界面：IL32在肝脏淋巴细胞和肝细胞中表达，在胃黏膜中在T细胞（71%）和上皮细胞（63%）中表达；LGALS3（肝脏巨噬细胞；胃上皮93%、髓系细胞70%）作为组成敏感性对照，其批量上调与上皮内在下调相反。第三种模式属于星状细胞：CADM2在两个肝脏数据集中独家定位到肝星状细胞。胃CADM2在批量CAG组织中可低丰度检测到，但在胃单细胞图谱中缺失，因此其胃细胞来源仍未解决。3.4 蛋白水平与公共组织注释的一致性：研究人员评估了Human Protein Atlas（HPA）中所有五个领先候选基因的细胞类型免疫组化注释与预先指定的 concordance 标准（一致、部分一致、不一致或不可评估；补充表S37）。ANXA4具有最强的蛋白水平支持：肝胆管细胞（三个抗体中两个）和胃腺细胞中中至高染色，所有抗体评级为Enhanced并与胆管细胞-上皮模式一致。CDHR2染色肝细胞（官方注释为中等；抗体间强度变异见补充表S37）和胃腺细胞，与其肝细胞-上皮定位匹配。LGALS3在胃腺细胞中高染色，HPA肝脏单细胞数据显示巨噬细胞和Kupffer细胞主导表达，与其肝脏巨噬细胞定位一致（HPA肝脏免疫组化不单独注释肝脏巨噬细胞）。对于IL32，无蛋白注释可用，两个抗体评级为Uncertain且低染色；即便如此，HPA单细胞数据与研究人员数据集中观察到的T细胞主导模式一致。CADM2在肝脏或胃免疫组化中未检测到，如其低转录丰度（0.8–1.7 nTPM）和限于稀有细胞类型所预期；HPA单细胞数据中唯一的肝脏信号落在肝星状细胞，这反映了研究人员分析中的独家星状定位。HPA分析正常而非疾病组织，因此该分析涉及细胞类型归属；它不是疾病阶段验证。3.4.1 患者水平纤维化阶段梯度：纤维化阶段来自GSE202379样本元数据中的SAF评分。五例终末期移植肝和两例健康供体无SAF评分；研究人员分别将其分配为F4和F0。47例患者覆盖F0（n=5）、F1（n=9）、F2（n=12）、F3（n=12）和F4（n=9）。在预先指定的每患者和细胞类型至少20个细胞的阈值后，肝细胞IL32梯度利用所有47例患者；星状细胞CADM2梯度利用37例患者，胆管细胞ANXA4梯度利用30例。在GSE202379内（纤维化阶段0–4），星状细胞伪批量中CADM2表达与纤维化阶段强烈负相关（Spearman ρ = ?0.761，p = 4.5 × 10?8，n=37），通过Bonferroni校正阈值（p = 4.8 × 10?4；13基因×八细胞类型），呼应其在晚期疾病中的批量下调。该估计与星状细胞转录的细胞内在变化、星状细胞活化状态转变或星状细胞亚群相对丰度变化相容，此处使用的参考注释单一星状细胞区室，因此无法区分三者。肝细胞伪批量中IL32与纤维化阶段正相关（ρ = 0.639，p = 1.4 × 10?6，n=47）。相比之下，胆管细胞伪批量中ANXA4与阶段负相关（ρ = ?0.585，n=30）。调整年龄和性别后梯度方向改变，三者均保持显著（补充表S46）。调整后，系数为星状细胞CADM2 ρ = ?0.75（p = 2.6 × 10?7）、肝细胞IL32 ρ = 0.62（p = 5.4 × 10?6）和胆管细胞ANXA4 ρ = ?0.50（p = 0.0069）。调整前通过Bonferroni校正阈值的两个梯度（CADM2和IL32）调整后仍低于阈值。与晚期疾病中胆管细胞比例上升一起解读，这看起来像胆管反应动态，更多胆道细胞每个表达更少ANXA4；这也调和了批量审计中ANXA4在细胞组成校正后丧失显著性。由于47例患者中七例携带推断而非记录的阶段，研究人员在排除这七例后重复每个梯度，留下40例SAF分级患者（补充表S49）。三个显著关联保持方向并仍显著：星状细胞CADM2 ρ = ?0.666（p = 5.9 × 10?5，n=30）、肝细胞IL32 ρ = 0.587（p = 6.8 × 10?5，n=40）和胆管细胞ANXA4 ρ = ?0.437（p = 0.033，n=24）。没有关联从显著变为不显著，两个不显著的仍不显著。记录的F4阶段包括四例患者，因此限制性分析在范围顶部效力较低。3.4.2 定位候选基因的组成敏感性：为检查候选信号是否由细胞类型比例变化携带，研究人员比较了四个批量队列（GSE126848、GSE130970、GSE162694和GSE174478）中未调整和MuSiC调整的效应。全基因运行未通过质量控制标准（平台基因组成不匹配；比例崩溃），因此标记限制性估计取而代之（补充表S35）。估计比例在每个队列中看起来生物学合理（平均肝细胞分数0.61–0.87；星状细胞分数0.05–0.28），研究人员重建标记基因表达的中位Spearman相关系数为0.44–0.50（补充表S33和S35）；GSE126848具有最低Pearson重建相关性（0.18），研究人员将其标记为反卷积敏感。对于CADM2，调整在三个具有可测量信号的队列中将下调减弱不超过2–11%，效应在每个队列中保持显著（GSE126848，p = 3.1 × 10?3；GSE130970，p = 1.5 × 10?3；GSE174478每阶段斜率，p = 2.0 × 10?3），这表明CADM2下调不能由估计细胞类型比例变化解释。在GSE162694中，发现阶段CADM2效应已经不显著（p = 0.28），两臂均无信号。ANXA4在病例对照比较中保持稳定（从27%增强到32%衰减；调整p ≤ 0.022），但在GSE174478纤维化阶段梯度上丧失86%效应并降至显著性以下（p = 0.59），与胆管反应相关表达一致。LGALS3作为预先指定的组成敏感性对照，在两个信号均可测量的地方显示比CADM2更大的衰减（GSE126848，16%对11%；GSE130970，36%对4%）并在纤维化阶段梯度上方向反转（两臂均不显著），这支持部分组成依赖性。未调整模型的独立重算返回与发现审计中标记评分校正相同的估计（GSE126848中三个基因的log2倍数变化差异≤0.03）。三个分析臂（发现主分析、标记评分校正和MuSiC调整重算）的完整结果，以及每样本比例和质量控制记录，见补充表S32–S36和图S2。3.4.3 探索性遗传三角验证：在69个两样本孟德尔随机化估计（10基因×七个肝脏和胃相关结局）中，没有一个在错误发现率（FDR）校正后以一致方向存活；唯一FDR显著的估计（ANO10对胃炎）方向不一致，因此不支持。IL32由三个血液顺式eQTL变异工具化（平均F统计量682），给出方向一致的胃轴估计但未达显著性。对于细胞状态限制性候选基因（CADM2、ANXA4和IL32），GTEx v10肝脏、胃、脾或全血中不存在组织匹配的顺式eQTL汇总数据。空间验证同样无法从公共处理数据实现：Human MASLD Spatial Multi-omics Atlas的处理矩阵无法直接下载，GSE312698定制面板未披露。这些失败反映公共数据可用性缺口，而非基因本身的特性。因此研究人员将遗传分析报告为探索性，并基于跨队列表达水平复制和细胞状态定位进行优先级排序（补充表S20–S24）。3.4.4 针对内部比较者的证据分级：审计产生不同地位的基因，因此研究人员按分析前固定的标准而非仅按排名分级。候选基因满足三个条件进入Tier 1：在肝脏筛选的所有六个留一数据集重复中排名前200；在至少一个独立保留队列中以预期方向变化且p < 0.05；第二条正交证据线支持它，形式为蛋白水平细胞类型归属、对细胞组成调整的稳健性或患者水平阶段梯度。Tier 2分配给通过稳定性标准但在复制或组成稳健性上不足的候选基因。三个基因达到Tier 1。CDHR2在GSE135251中变化log2倍数变化1.43（p = 1.2 × 10?5），在GSE153224中变化2.16（p = 0.008），bootstrap曲线下面积0.911（95% CI 0.860–0.956），具有一致的肝细胞和胃腺蛋白归属。CADM2变化?0.50（p = 8.2 × 10?6）和?1.65（p = 0.008），MuSiC调整后丧失不超过11%效应，并携带任何候选基因中最陡的阶段梯度（星状细胞ρ = ?0.76）。IL32在肝脏队列中变化1.26（p = 1.5 × 10?3），在胃队列中变化0.49（p = 0.55），在十个队列中保持方向，曲线下面积0.798（95% CI 0.684–0.900）。两个基因被分配到Tier 2。ANXA4在两个保留队列中均不显著（p = 0.48和p = 0.15），组成调整后丧失86%阶段梯度（调整p = 0.59），尽管其蛋白水平证据是面板中最直接的；LGALS3仅在胃队列中达到显著性（p = 0.032），调整后阶段梯度方向反转，并预先指定为组成敏感性对照。补充表S48给出所有14个候选基因的分级。八个降级基因提供分级所依赖的内部比较，因为它们每个都经过相同的审计。它们中没有因缺乏名义显著性而失败。TSPAN3和GOLM1在所有六个肝脏队列中保留闭式p值低于0.001，但分别在六个留一数据集重复中仅两个和一个保持前200排名；FGA、FGB和SLC6A16在任何重复中均未保持该排名。GOLM1（GP73）是有启发性的案例，因为其作为肝细胞癌和晚期肝病血清标志物有广泛文献，此处其跨队列排名依赖于个别队列而非可重复效应。RPS6KA1被置于观察列表而非降级：它在六个重复中四个排名稳定，但未被单核参考面板捕获，因此无法分配来源细胞。分级陈述从公共批量 and 单细胞数据可获得的证据强度；它不是生物学重要性排名，两个Tier 2基因仍是表达支持的候选基因。
讨论部分总结：研究人员着手询问MASLD和CAG是否共享可重复的转录变化。这两种疾病在其器官的逐步进展中占据类似位置，因此寻找共享分子足迹是合理的。五个基因（IL32、CDHR2、ANXA4、CADM2和LGALS3）通过预先指定的稳健性审计，在两个保留队列中复制，并分为三种组织或细胞状态模式。探索性孟德尔随机化在研究人员能运行的所有分析中未给出一致的FDR支持遗传证据，因此研究人员将这些基因讨论为优先候选基因而非确认驱动因素。第六个基因RPS6KA1保持在预先指定的观察列表上。其排名跨队列保持，但单核参考面板未捕获它，因此无法分配来源细胞。总之，三种模式对应于两个器官中发生的组织反应：上皮损伤伴代偿性分化（CDHR2、ANXA4）、上皮界面的免疫反应（IL32）和星状细胞重塑（CADM2）。LGALS3是浸润携带信号的指定对照。五个基因中四个映射到在另一器官中有明确对应物的细胞区室；例外CADM2没有分配胃细胞类型。汇聚在于模式而非任何单个基因，共享跨器官生物学的论据就在于此。上皮损伤、免疫浸润和成纤维细胞活化也是肠-肝轴被认为运作所通过的反应。肝脏表达研究（大多在旧NAFLD命名下进行）已在单队列中绘制严重程度相关特征。与该文献最清晰的接触点是IL32。在一个活检系列中，它排名为严重NAFLD中最强上调转录本，与脂肪变性分级、小叶炎症、气球样变和纤维化相关。一个独立系列将肝脏IL-32与NAS联系起来，并显示重组IL-32损害原代人肝细胞中的胰岛素信号。研究人员的数据增加两个观察。第一，上调超出肝脏：在CAG中，IL32定位到T细胞和上皮细胞。这符合幽门螺杆菌感染胃上皮中IL-32的NF-κB依赖性诱导。第二，肝细胞IL-32在患者水平随纤维化阶段上升（Spearman ρ = 0.64），这将先前的病例对照关联推进一步。CAG的转录组工作更薄弱，大多数分析依赖单一数据集。GSE153224的原始分析通过共表达和蛋白相互作用网络在八个基因集中提名CDHR2。CDHR2此处从六个肝脏和两个胃队列中出现，未参考该流程，并在同一数据集的混合文库再分析中保持。信号在不同分析路线中存活。然而该汇聚来自单一数据集，因此不能作为独立队列复制。CDHR2（PCDH24）在此背景外已知是分化肝、肾和结肠上皮的钙黏蛋白，在相应癌中沉默。研究人员的数据约束解释。如果MASLD和CAG中的上调是随损伤规模化的代偿反应，人们会期望纤维化阶段梯度，但肝细胞CDHR2没有显示（图5）。增加似乎早期饱和，或反映病因而非严重程度。已建立癌中报告的丧失则会在后来出现，作为转化相关事件。ANXA4先前被描述为幽门螺杆菌诱导的胃癌相关蛋白。在研究人员的肝脏数据中，它几乎完全位于胆管细胞，随着纤维化进展每个细胞表达下降而胆管细胞比例上升。这种组合符合胆管反应，其中胆道和祖细胞区室在脂肪性肝炎纤维化进展中扩张。ANXA4的批量上调将反映更大的胆道区室，其个别细胞表达更少该基因。同一模型解释为什么组成校正消除纤维化阶段梯度而保持病例对照信号完整。一个限定适用：胆管细胞定位、上升比例和反卷积面板均来自同一单核参考（GSE202379）。胆管反应解释是内部一致的而非独立确认的，需要在无关胆道数据集中检验。LGALS3做了指定组成敏感性对照应该做的：在批量组织中上调，主要定位到巨噬细胞，反卷积调整后衰减，并相对纤维化阶段方向反转。Galectin-3在实验模型中驱动肌成纤维细胞活化和肝纤维化，它是belapectin的药理靶点，后者已在NASH肝硬化中完成2b期试验。一个众所周知的浸润携带基因以完全预期的方式衰减支持审计的校准。它也警告不要将任一器官中的批量LGALS3上调解读为上皮事件。胃单细胞数据中LGALS3的上皮内在下调与其批量变化相反，仍未解释。一个不一致也值得评论：LGALS3在单个发现队列（Arendt；log2倍数变化?0.24，p = 0.013）中方向不一致。该队列在病例组合和对照确定方面与其他肝脏队列不同，组成携带信号正是预期对此类差异敏感的信号。研究人员将此例外解读为该基因浸润驱动行为的结果；它不反对整体模式，整体模式主要基于IL32。CADM2是五个中最少被表征的。它最著名的是与体重指数的遗传关联，在小鼠中删除Cadm2减少肥胖并改善胰岛素敏感性。据研究人员所知，没有早期工作将其与肝纤维化联系起来。其限于肝星状细胞、与纤维化阶段的陡峭负相关（ρ = ?0.76）及其对反卷积调整的不敏感性使其成为研究人员识别的最细胞特异性信号。它也是机制上最难放置的。表达数据不能说明星状细胞CADM2丧失是伴随、允许还是对抗活化。注释星状细胞区室内的亚型转变（从静息向活化细胞）是研究人员数据无法排除的替代解释，胃细胞来源仍未解决。如所述，探索性孟德尔随机化未产生一致的FDR支持遗传证据。最接近支持信号来自IL32，由三个强血液顺式eQTL变异工具化（平均F = 682）。其胃轴估计方向一致但不显著。这些数据不能区分有限效力与真实无效应。三个进一步考虑反对对零结果的字面解读。第一，工具是血液顺式eQTL，而五个领先基因中三个是细胞状态限制性的，根本没有组织匹配工具。表达的遗传调控强烈依赖组织和细胞组成，因此血液工具至多是星状细胞或胆管细胞转录本的弱代理。第二，孟德尔随机化估计表达的终生平均效应，而此处报告的信号以疾病状态为条件。第三，CAG没有足够效力的GWAS，因此胃炎和胃癌必须作为胃轴代理。暴露工具和所有七个结局GWAS均来自欧洲祖先样本，因此估计可能不推广到其他祖先群体。这些是数据和工具可用性的缺口，不是反对因果性的证据。 underlying队列的两个特征为如何解读这些共享信号设定了界限。第一，三个胃对比中两个是针对炎症黏膜而非健康黏膜：GSE116312中的滤泡性胃炎和GSE153224中的慢性非萎缩性胃炎。第三个GSE27411较小，比较萎缩与非感染胃。因此共享信号主要位于炎症背景上。该设计使信号保守，但不能声称对萎缩特异。其中一些可能只是标记慢性黏膜炎症。第二，肝脏队列的幽门螺杆菌状态不可用，胃队列仅部分解决。当前感染常随萎缩进展而丢失。在一个经典系列中，86%的萎缩性胃体胃炎受试者血清阳性，而仅33%染色可检测到微生物。因此胃信号既不能归因于也不能干净地与活动性感染分离。鉴于幽门螺杆菌感染与MASLD之间的荟萃分析关联，在幽门螺杆菌分层队列中复制将是必要的。在此类队列存在之前，共享基因的两种解读都保持开放：它们可能反映器官自主模式，或可能记录共享暴露的足迹。未测量的幽门螺杆菌感染值得单独处理，因为它是可能生成整个模式的唯一混杂因素。感染在研究人员使用的对比中分布不均。萎缩队列按构造富集过去或现在感染；非萎缩比较者不是，且三个肝脏队列均未记录血清学。这种大小的患病率差异转化为研究人员模型中任何项都无法吸收的偏倚，研究人员没有调整它：审计以队列为条件，而非感染。三个观察限制它可解释多少一致性，尽管没有一个消除担忧。第一，候选基因没有统一响应。LGALS3丧失其阶段梯度并在组成调整后反转，RPS6KA1和GOLM1未通过排名稳定性，ANXA4在两个保留队列中均未复制，而单一感染驱动炎症轴预期会一起移动一组相关基因。第二，探索性孟德尔随机化对感染状态不敏感，因为它基于种系工具，它未返回一致的FDR支持效应，因此共享信号不由感染会混杂的遗传路径携带。第三，最强候选基因还在来自不同群体和平台的两个进一步队列中复制，这是稳定效应所做的而单队列暴露伪影不会做的。明确区分两种解读需要具有CagA或粪便抗原血清状态的队列，通过设计或多变量调整分层；研究人员认为这是最有信息量的下一个数据集。本研究的优势包括审计的预先指定，降级候选基因被报告而非静默丢弃，以及五个基因的外部验证而不重新选择。患者水平伪批量设计避免将细胞视为独立观察。降级本身是有信息量的。GOLM1（GP73）作为肝细胞癌和晚期肝病血清标志物有广泛文献支持；即便如此，其信号此处由个别队列驱动，未通过审计。如此确立的标志物基于这些理由降级显示单队列肝脏发现泛化有多差，并相应给予通过相同审计存活的五个基因更多权重。局限性同样具体。设计是横断面的和基于批量的，因此不能建立沿Correa级联或MASLD谱的方向性。胃发现队列较小，外部胃队列只能作为每组五个混合文库分析。保留肝脏队列中CDHR2的判别估计仅基于十个对照。胃单细胞分析仅描述性，全基因反卷积未通过自己的质量控制标准，必须限制于标记基因。因为HPA免疫组化分析正常组织，蛋白水平分析涉及细胞类型归属而非疾病，对于IL-32抗体证据评级为不确定。从公共数据不可能进行空间验证。最后，IL-32在小鼠中不表达，这限制标准动物模型中的机制测试。五个中没有一个在疾病组织中在蛋白水平验证。胃保留队列中的复制必须同样谨慎解读。所有五个候选基因的方向一致性不同于统计支持的复制：每组五个混合文库，只有CDHR2、CADM2和LGALS3达到显著性，而ANXA4和IL32仅方向一致，文库而非患者是该比较中的统计单位。研究人员还分离独立证据与支持证据。独立复制仅声称用于两个保留批量队列和探索性孟德尔随机化的种系工具。伪批量纤维化阶段分析和反卷积敏感性分析均来自GSE202379，反卷积参考本身源自该数据集，因此它们之间的一致是内部的而非外部的；Human Protein Atlas结果则是正常组织中的细胞类型归属，不是任一疾病中的差异蛋白表达。两个分析均报告因为它们约束如何解读一致性，但两者都不计为复制。直接随访从这些缺口生长：在来自两个器官的独立、组织学分期活检系列中对五基因面板进行qPCR和免疫组化确认。空间转录组学可解决CADM2的胃来源，IL32将需要基于人类的实验系统。循环IL-32水平已被报告跟踪NAFLD严重程度，因此面板的血清基础扩展是合理的转化路线。共享特征是否在任一器官中分层进展风险是纵向队列的问题，不是此处分析的横断面数据的问题。机制步骤最重要，它现在在人类系统中是可处理的。在气液界面生长的患者衍生胃类器官保留腺体结构，可跨transwell膜与肝或胆管细胞类器官共培养，这将允许CDHR2和ANXA4逐个区室扰动，通过一次一个类器官系的CRISPR干扰，同时两个区室交换可溶性因子。CADM2可在含星状细胞的肝脏类器官或精密切割肝切片中接近，其中它跟踪的纤维化梯度被保留。IL32决定模型选择，因为该基因没有小鼠直系同源物：用人造血细胞重建的人源化小鼠，或用原代或iPSC衍生髓系和T细胞的类器官-免疫共培养，是它可被测试的系统。这些设计中的每一个 interrogate 此处报告模式的一个元素，在兩個器官可单独操纵的系统中。结论部分翻译：据研究人员所知，没有早期研究在转录组水平跨器官比较代谢功能障碍相关脂肪性肝病和慢性萎缩性胃炎。整合八个发现队列使研究人员能够优先考虑两种疾病共享的五个候选基因。这些基因追踪三种组织和细胞状态模式：上皮模式（CDHR2、ANXA4）、上皮界面的免疫相关模式（IL32）和星状细胞相关模式（CADM2）。所有五个在保留队列中保持其发现阶段方向；其中三个（CDHR2、CADM2和IL32）满足较高证据层级的所有三项标准，并且每个在来自肝和胃的单细胞数据中映射到特定细胞类型。探索性孟德尔随机化再次未产生一致的FDR支持遗传证据。因此，五个基因仍是表达支持的候选基因，等待机制确认而非确认驱动因素。研究人员所做的声称相应狭窄：这些数据确立的是沿肠-肝轴的跨疾病转录一致性，通过共享细胞状态定位解读，而非共享病因、因果关系或共同生物标志物。审计策略应可重用于其他共病疾病的跨器官比较。
