《Genes》：Genetic Evidence for Metabolites Mediated Causal Inference Between Gut Microbiota and Venous Thromboembolism
编辑推荐：
摘要
背景： 细胞异质性限制了复杂疾病中的因果推断。研究人员应用了一套已建立的细胞类型分层孟德尔随机化（cell-type-stratified Mendelian randomization, csMR）框架，该框架最初由Hao等人（2024）开发，用以研究
摘要
背景： 细胞异质性限制了复杂疾病中的因果推断。研究人员应用了一套已建立的细胞类型分层孟德尔随机化（cell-type-stratified Mendelian randomization, csMR）框架，该框架最初由Hao等人（2024）开发，用以研究脂质表型是否通过不同的大脑细胞群影响肾脏疾病，并扩展了全蛋白质组中介分析组分。方法： 利用贝叶斯共定位（后验概率假设4 [posterior probability of hypothesis 4, PPH4]，≥0.8）和多维工具变量筛选，研究人员评估了五种脂质性状在十种大脑细胞或组织分层中与肾脏疾病风险的细胞类型分层孟德尔随机化支持的关联。工具变量来源于已发表的全基因组关联研究（genome-wide association studies, GWAS）和人类大脑单细胞表达数量性状基因座（expression quantitative trait locus, eQTL）图谱。一项探索性蛋白质组中介分析使用差异法（Difference Method）进行，以UKB-PPP作为发现资源，deCODE作为外部验证资源，中介信号被解释为产生假设性的。对于主要的csMR分析，多重检验在每个脂质性状内分别使用Bonferroni阈值3.33 × 10?4进行校正，对应于每个脂质性状的150项细胞/组织-结果检验。结果： 在csMR分析中，胆固醇相关脂质显示出与肾脏疾病的细胞类型分层关联，最显著的是在少突胶质细胞相关分析中（对数优势比尺度上的最大β = 1.025）。低密度脂蛋白胆固醇（low-density lipoprotein cholesterol, LDL-C）相关的广泛慢性肾小球疾病关联在兴奋性神经元相关分析中最为明显。双队列蛋白质组分析优先考虑了血浆蛋白，包括SNAP29和ICAM4，作为探索性中介分析的候选蛋白相关信号。这些发现代表了遗传共定位和孟德尔随机化支持的关联，而非实验确立的机制，应被解释为产生假设性的。结论： 本研究提供了细胞类型分辨的遗传和蛋白质组关联图谱，将脂质性状、大脑细胞分层和肾脏疾病结局联系起来，为未来的实验验证提供了假设。
论文解读
一、研究背景与问题提出
中枢神经系统作为全身稳态调节的关键终端，通过神经内分泌反馈和自主神经调节协调多维生理网络。尽管下丘脑-垂体-肾上腺轴在代谢调节中的作用已得到公认，但新出现的证据表明存在超越经典内分泌通路的更复杂的“脑-外周”串扰。遗传学研究日益将脑脂质稳态与全身性病理联系起来，然而脑相关遗传调控信号与颅外代谢性疾病（尤其是肾脏疾病）之间的关系仍知之甚少。
解读这类跨器官通讯的一个主要挑战是细胞异质性。传统的组织水平分析掩盖了不同脑细胞类型（如神经元、星形胶质细胞和小胶质细胞）的独特遗传结构。单细胞RNA测序和细胞类型特异性表达数量性状基因座（eQTL）图谱的近期进展表明，近一半的表达调控位点（eGene）表现出细胞类型依赖性效应。这提示全基因组关联研究（GWAS）中鉴定的疾病相关变异可能通过特定的细胞亚群而非整个组织发挥致病作用。
多项证据表明，不同的脑细胞类型参与全身代谢和肾脏病理生理过程。中枢神经炎症和交感神经过度活跃是肾脏疾病的既定危险因素，星形胶质细胞和小胶质细胞被认为是这些过程的关键调节因子。少突胶质细胞对髓鞘形成具有独特的高胆固醇需求，通过GPR17–cAMP–PKA轴将中枢脂质代谢与全身能量平衡耦合。兴奋性和抑制性神经元群，特别是室旁核和迷走神经背运动核中的神经元，直接控制肾交感神经和迷走神经输出。内皮细胞和周细胞作为血脑屏障转运的守门人，与慢性肾脏病伴随的脑微血管功能障碍有关。尽管有这些积累的生物学证据，但系统性的群体水平因果推断将特定脑细胞类型与肾脏疾病联系起来的研究尚未开展。
孟德尔随机化（MR）利用遗传变异作为工具变量，最大限度地减少混杂和反向因果，已成为因果推断的稳健框架。然而，传统MR依赖组织水平数据，无法解析遗传风险表现的细胞层机制。以往研究虽然考察了循环脂质与肾脏结局的关系，但将大脑视为黑箱，忽略了脂质相关遗传信号可能注释到脑细胞群体的细胞类型特异性调控景观。因此，缺乏脂质介导的脑-肾轴的系统性、细胞分辨因果框架。脑细胞类型相关遗传调控如何与外周肾脏结局相关仍不清楚。循环蛋白提供了可测量的全身分子性状，可能有助于优先排序候选跨器官通路。整合细胞特异性eQTL与大规模蛋白质组数据，为绘制脑-血浆信号轴提供了机会。
基于此，研究人员应用已建立的csMR框架，考察脂质相关遗传信号注释到人类脑细胞类型后是否与肾脏结局相关，通过整合脑单细胞eQTL图谱、GWAS汇总统计和循环蛋白质组数据，构建了连接脂质性状、脑细胞分层、血浆蛋白和肾脏疾病结局的细胞类型分层遗传和蛋白质组关联框架。
二、主要技术方法
研究人员应用了Hao等人（2024）开发的细胞类型分层孟德尔随机化（csMR）框架，该框架实施三级分层推断系统：一是变异精细定位和遗传共定位，采用SuSiE算法和Coloc贝叶斯框架，识别GWAS与单细胞eQTL数据之间的共享因果变异，PPH4 ≥ 0.8为显著阈值；二是工具变量精炼，涉及多维过滤以选择高置信度遗传工具；三是多算法因果估计，量化暴露对肾脏相关结局的因果效应。脑单细胞eQTL汇总数据来自Bryois等人（2020），涵盖八种主要脑细胞类型，来源于192名欧洲裔个体的前额叶和颞叶皮层。组织水平脑eQTL数据来自GTEx v8项目。脂质性状GWAS数据来源于英国生物银行（Neale实验室，N ≈ 340,000–360,000）和全球脂质遗传学联盟（GLGC）。肾脏结局涵盖十五种疾病，数据来源于FinnGen R10、UK Biobank和公开GWAS荟萃分析。蛋白质组分析以UKB-PPP为发现队列，deCODE为外部验证队列。工具变量筛选遵循MR三大假设，采用LD聚类（r2 < 0.001）、PhenoScanner混杂表型过滤、F统计量 > 10等标准。因果估计以逆方差加权（IVW）为主要方法，MR-Egger、加权中位数、MR-RAPS和加权众数为补充敏感性估计器。蛋白质中介分析采用差异法估计统计中介效应。
三、研究结果
3.1 总胆固醇与脑细胞类型eQTL的共定位
基因组研究初步鉴定出254个与总胆固醇（TC）相关的近独立先导单核苷酸多态性（SNP）。共定位分析观察到259个跨脑细胞分层的SNP-细胞类型共定位事件。17.4–89.2%的共定位SNP和14.6–80.5%的共定位基因表现出细胞类型特异性。跨细胞类型eQTL效应大小分析显示，共定位SNP在其主要作用细胞类型之外的环境中效应一致减弱，支持这些脂质相关变异对基因表达具有细胞类型特异性调控。77.8%的共定位基因表现出细胞类型限制性表达模式。
3.2 总胆固醇的细胞类型分层因果效应
使用全基因组显著变异作为工具变量的主要MR分析确立了TC与高血压肾病、慢性肾脏病（CKD）和IgA肾病的因果关联，LD评分回归进一步证实了这些关联。多重检验校正后，TC在多个脑细胞背景下与肾乳头状腺癌、慢性肾小球肾炎和慢性肾脏病显示出细胞类型特异性关联，涉及内皮细胞、星形胶质细胞、皮质组织、额叶皮质、兴奋性神经元、抑制性神经元、小胶质细胞、少突胶质前体/定向少突胶质细胞谱系细胞、少突胶质细胞和周细胞。少突胶质细胞在TC相关广泛慢性肾小球疾病风险中的显著参与可能反映了髓鞘形成的高胆固醇需求。
3.3 极低密度脂蛋白胆固醇的共定位与因果效应
GWAS分析鉴定出59个与VLDL-C相关的近独立先导SNP，共定位鉴定出46个SNP影响约40个基因，27.2–60.0%的SNP表现出细胞类型特异性。csMR分析显示VLDL-C与慢性肾小球肾炎（CGN）呈细胞类型分层模式关联。周细胞分层显示与CGN呈负相关（β = ?0.688，对数优势比尺度；p = 3.23 × 10?26），而内皮细胞、小胶质细胞和少突胶质细胞显示正对数优势比估计。这表明VLDL-C关联方向在不同脑细胞分层间存在差异。
3.4 高密度脂蛋白胆固醇的共定位与因果效应
GWAS分析鉴定出455个与HDL-C相关的近独立先导SNP，共定位鉴定出160个SNP。csMR分析显示HDL-C在星形胶质细胞、皮质组织和小胶质细胞中与CGN呈正对数优势比关联，代表性估计包括星形胶质细胞（β = 0.479；p = 2.84 × 10?16）、皮质组织（β = 0.329；p = 6.33 × 10?27）和小胶质细胞（β = 0.237；p = 1.83 × 10?6）。与胆固醇相关脂质不同，HDL-C主要通过胶质细胞发挥作用，提示炎症而非代谢机制。
3.5 低密度脂蛋白胆固醇的共定位与因果效应
GWAS分析鉴定出194个与LDL-C相关的近独立先导SNP，共定位鉴定出58个SNP。csMR分析显示LDL-C在多个脑细胞或组织分层中与CGN呈正对数优势比关联，包括额叶皮质（β = 0.985；p = 9.52 × 10?96）、内皮细胞（β = 1.035；p = 1.08 × 10?113）和兴奋性神经元（β = 0.961；p = 2.69 × 10?132）。兴奋性神经元特异性LDL信号提示LDL-C可能激活神经元通路，增强肾交感神经输出。
3.6 甘油三酯的共定位与因果效应
GWAS分析鉴定出392个与甘油三酯（TG）相关的近独立先导SNP，共定位鉴定出145个SNP。csMR分析显示TG升高通过星形胶质细胞（β = 0.264；p = 1.75 × 10?4）和脑皮质（β = 0.167；p = 6.81 × 10?6）增加广泛慢性肾小球疾病风险，效应量小于胆固醇相关脂质。
3.7 甘油三酯对肾脏结局的细胞类型特异性效应
多重检验校正后，TG在星形胶质细胞、皮质组织、额叶皮质和少突胶质前体/定向少突胶质细胞谱系细胞中与CGN显示出细胞类型特异性关联。在这些细胞背景下，遗传预测的TG升高一致地与疾病风险增加相关，支持TG代谢在CGN中收敛但细胞类型依赖的作用。
3.8 肾脏疾病的可重复血浆候选蛋白中介因子
对于慢性肾脏病，UKB-PPP发现队列和deCODE外部验证鉴定出三个可重复血浆蛋白中介因子：TIMD4、IL1RN和B4GALT1。TIMD4显示保护性关联，而IL1RN和B4GALT1显示风险增加关联，涉及免疫调节、炎症信号和细胞外基质重塑通路。对于肾和输尿管结石，鉴定出22个可重复血浆蛋白，其中SNAP29显示显著的保护性关联，可能反映囊泡运输和膜融合相关保护性蛋白特征。对于慢性肾脏病，鉴定出18个可重复血浆蛋白，ICAM4显示显著的风险增加关联，提示细胞黏附和炎症通路参与。对于肾乳头状腺癌，23个疾病相关血浆蛋白中有4个在deCODE中重复，包括CD14、OXT、CTSF和TP53I3，涉及免疫信号、神经内分泌调节、溶酶体蛋白酶活性和p53相关氧化应激反应。
3.9 探索性蛋白质组中介信号
csMR分析鉴定出54个脂质-蛋白关联，涉及TIMD4、ICAM4、OXT、CTSF、TP53I3和SNAP29等关键蛋白。优先化的探索性中介分析鉴定出8个候选统计中介信号，涉及6种蛋白。TIMD4出现在3个候选通路中，连接抑制性神经元、皮质组织或额叶皮质中的胆固醇相关信号与CGN相关结局。ICAM4被纳入连接少突胶质细胞谱系分层总胆固醇信号与慢性肾脏病的候选统计中介信号，估计中介比例约为18.5%。SNAP29与小胶质细胞相关的肾和输尿管结石通路相关。OXT、CTSF和TP53I3在内皮细胞相关的总胆固醇通路中被优先化。多个间接效应为负或与总效应方向相反，被解释为反向或抑制子样统计中介信号。
四、讨论总结与结论翻译
讨论部分指出，本研究将已建立的csMR框架应用于脂质相关脑细胞分层与肾脏疾病的关联研究。传统组织水平分析可能通过平均化不同细胞群的遗传调控效应而掩盖细胞类型特异性信号。研究发现应被解释为细胞类型分辨的遗传证据，而非脂质诱导肾脏疾病物理上起源于大脑的证明。胶质细胞介导的机制方面，少突胶质细胞谱系细胞和周细胞可能参与脂质相关遗传效应对肾脏结局的影响，与既往MR研究一致。神经元分层关联方面，兴奋性神经元和抑制性神经元显示出方向不同的LDL相关关联，支持神经元“双向切换”模型。内皮细胞和周细胞分层分析突出了候选血管或血管周围细胞分层。候选蛋白中介因子方面，TIMD4、IL1RN和B4GALT1提示免疫-炎症和组织重塑特征，SNAP29提示保护性运输相关或自噬相关通路，ICAM4提示细胞黏附相关炎症过程，CD14、OXT、CTSF和TP53I3指向免疫、神经内分泌、溶酶体和肿瘤应激通路。局限性包括潜在样本重叠和赢家诅咒偏倚、祖先偏倚、主要利用顺式eQTL、静态横断面数据以及需要功能验证。
研究结论为：研究人员应用细胞类型分层孟德尔随机化（csMR）框架，整合大规模遗传学与脑单细胞表达数量性状基因座（eQTL）图谱，系统绘制了“脂质→脑细胞→血浆蛋白→肾脏疾病”关联框架。主要发现包括：脂质对肾脏疾病的遗传效应在不同脑细胞和组织分层间存在差异，少突胶质细胞和星形胶质细胞相关分析对胆固醇相关性状显示显著关联，而神经元分层显示方向不同的LDL相关关联；CGN作为反复出现的肾脏结局，但被分析为广泛慢性肾小球疾病表型而非亚型特异性机制；通过整合血浆蛋白质组GWAS，优先化了SNAP29、ICAM4、OXT和TP53I3等候选血浆蛋白，ICAM4与慢性肾脏病的探索性中介信号估计中介比例约为18.5%；csMR框架克服了传统组织均质化分析的分辨率限制，有助于生成关于多器官调控网络的更精细假设。这些发现应被解释为产生假设的遗传和蛋白质组证据，而非确认的生物学机制。
