《Biomedicines》：Identification of Genes Associated with Gut Microbiota Metabolites in Type 2 Diabetes Mellitus Complicated by Acute Cerebral Infarction Based on Bioinformatics and Mendelian Randomization
编辑推荐：
摘要
背景： 2型糖尿病（Type 2 diabetes mellitus, T2DM）与急性脑梗死（Acute cerebral infarction, ACI）是全球流行的疾病，造成了沉重的公共卫生负担。尽管肠道菌群（gut microbiota, GM）
摘要
背景： 2型糖尿病（Type 2 diabetes mellitus, T2DM）与急性脑梗死（Acute cerebral infarction, ACI）是全球流行的疾病，造成了沉重的公共卫生负担。尽管肠道菌群（gut microbiota, GM）日益被认为可调节与这两种疾病相关的代谢和免疫稳态，但其共享的遗传预测微生物类群以及T2DM与ACI中潜在的代谢物-基因调控机制仍不清楚，凸显了系统性研究的必要性。因此，本研究旨在阐明其调控机制。方法： 从公共数据库获取转录组数据集、全基因组关联研究（genome-wide association study, GWAS）数据和GM数据。采用两样本孟德尔随机化（Mendelian randomization, MR）和16S rRNA测序鉴定关键微生物类群。通过公共数据库预测代谢物-靶点基因（metabolite-target genes, MTGs）。基于Lipinski规则和低毒性筛选关键代谢物，并通过分子对接优先确定关键基因。利用公共数据库预测调控网络和潜在药物。结果： 以FinnGen糖尿病终点作为T2DM相关遗传易感性的替代指标，MR分析鉴定出两种疾病共有的12种候选微生物类群。进一步与16S rRNA测序的前30个微生物类群取交集，确认了两个关键微生物类群（Alistipes 和 Bifidobacterium）。随后，针对源自关键微生物类群的代谢物，鉴定出602个MTGs。经过成药性评估和毒性分析，确定了9种具有良好成药性和低毒性的关键代谢物。分子对接预测了这些关键代谢物与编码蛋白的结合亲和力，表明ANPEP和CYP1B1是候选关键基因。此外，ANPEP与12个转录因子（transcription factors, TFs）相互作用，而CYP1B1受6个TFs和13个微小RNA调控，后者进一步与26个长链非编码RNA相互作用。分别鉴定出68种和32种靶向ANPEP和CYP1B1的潜在药物，为转化研究提供了候选药物。结论： 本研究系统性地将ANPEP和CYP1B1鉴定为T2DM相关遗传易感性与ACI共病状态下GM相关调控网络中的候选关键基因，为开发基于GM的靶向治疗提供了候选分子靶点和候选化合物。
论文解读：基于多组学整合策略探究2型糖尿病并发急性脑梗死的肠道菌群代谢物调控机制
一、研究背景与科学问题
2型糖尿病（T2DM）与急性脑梗死（ACI）是全球范围内患病率持续攀升的重大慢性疾病，二者常合并存在，显著增加患者的致残率和死亡率。流行病学证据表明，糖尿病患者发生ACI的风险约为非糖尿病患者的2倍，且血糖控制虽能减少微血管并发症，却未能有效降低卒中风险，提示传统降糖策略在脑血管保护方面存在局限性。近年来，肠道菌群（GM）及其代谢产物在T2DM与ACI发病中的作用受到广泛关注。研究表明，T2DM相关的肠道菌群失调可通过破坏肠黏膜屏障、促进脂多糖入血、损害血脑屏障等途径加剧脑损伤。然而，目前尚缺乏系统性研究从遗传预测层面揭示介导T2DM与ACI共病的共享肠道微生物类群及其代谢物-基因调控网络。为此，研究人员整合了转录组学、孟德尔随机化（MR）、16S rRNA测序及分子对接等多维度分析方法，旨在阐明GM代谢物在该共病状态下的分子机制，并筛选潜在的治疗靶点与候选药物。
二、主要技术方法概述
研究人员首先利用两样本孟德尔随机化（MR）分析，以GM为暴露因素、T2DM与ACI为结局变量，从MiBioGen数据库获取GM的全基因组关联研究（GWAS）数据，并以FinnGen数据库的糖尿病终点和UK Biobank的脑梗死数据作为疾病GWAS来源。同时，在中国昆明地区第九二〇医院招募了20例ACI患者和20例T2DM合并ACI患者，采集粪便样本进行16S rRNA测序（Illumina NovaSeq 6000平台），经QIIME2流程处理获得扩增子序列变异（ASV）丰度表，并基于SILVA数据库进行物种注释。随后，利用随机森林（RF）分析筛选区分两组的关键菌属。通过gutMGene数据库获取关键菌属的代谢物信息，并利用SEA和SwissTargetPrediction平台预测代谢物靶点基因（MTGs）。差异表达分析基于GEO数据库的GSE20966（T2DM胰腺β细胞标本）和GSE16561（ACI全血样本）两个转录组数据集。最后，采用SwissADME和ADMETlab平台评估代谢物的类药性与毒性，并通过CB-Dock2进行分子对接验证。
三、主要研究结果
3.1 候选微生物群落的遗传预测关联
IVW分析显示，34种GM与ACI存在遗传预测关联，24种GM与T2DM存在遗传预测关联。其中，Bifidobacterium（OR=0.999）被识别为ACI的保护因素，而Alistipes（OR=1.006）为ACI的危险因素；对于T2DM，Alistipes（OR=0.997）为保护因素，Bifidobacterium（OR=1.001）为危险因素。取交集后获得12种共享候选微生物类群。敏感性分析（Cochran's Q检验、MR-Egger截距检验、留一法分析）及Steiger方向性检验均支持结果的稳健性与因果方向的正确性。
3.2 T2DM-ACI患者的肠道菌群特征
α多样性分析显示ACI组与T2DM-ACI组间无显著差异，但β多样性分析（基于未加权UniFrac距离的主坐标分析PCoA）及置换多元方差分析（PERMANOVA）证实两组间菌群结构存在显著分离（R2=0.039，p=0.016）。线性判别分析效应大小（LEfSe）鉴定出10个在ACI组富集的菌属（如Muribaculaceae_g_Unclassified、Haemophilus）和4个在T2DM-ACI组富集的菌属（如Eisenbergiella、Slackia）。随机森林分析显示Bifidobacterium具有最高的变量重要性评分，是区分两组的最关键菌属。
3.3 关键微生物群落的确定
将MR分析获得的12种候选GM与随机森林分析排名前30的菌属取交集，最终确定Alistipes和Bifidobacterium为关键肠道微生物类群。
3.4 候选基因的筛选与分析
基于关键菌属预测得到43种代谢物和602个MTGs。差异表达分析在GSE16561中获得676个差异表达基因（DEGs），在GSE20966中获得1061个DEGs。取两数据集共有的27个一致性DEGs（23个上调、4个下调）与602个MTGs的交集，获得4个候选共病基因：ANPEP、CTSS、CYP1B1和DYRK2。基因本体论（GO）富集分析显示候选基因主要参与“小梁网发育”、“内溶酶体腔”和“类固醇羟化酶活性”等条目；京都基因与基因组百科全书（KEGG）通路富集分析提示涉及“吞噬作用”、“色氨酸代谢”、“类固醇激素生物合成”等通路。
3.5 关键基因的筛选与验证
对与4个候选基因相关的63种代谢物进行类药性评估，54种符合Lipinski五规则；进一步毒性分析筛选出9种无或低毒的核心代谢物。将核心代谢物对应的359个靶基因与4个候选基因取交集，获得3个核心基因：ANPEP、CTSS和CYP1B1。分子对接结果显示，ANPEP编码蛋白与5种代谢物（如3-（3-羟基苯基）丙酸、瓜氨酸、蛋氨酸等）具有预测的良好结合亲和力（结合能<-4.25 kcal/mol），CYP1B1编码蛋白与3-羟基-4-甲氧基苯丙酸（HMCA）表现出最佳结合活性，而CTSS蛋白与代谢物的结合力较弱。因此，ANPEP和CYP1B1被优先确定为候选关键基因。
3.6 关键基因的调控网络与药物预测
调控网络分析显示，CYP1B1可与13个微小RNA（miRNA）和26个长链非编码RNA（lncRNA）相互作用，构建了lncRNA-miRNA-mRNA调控网络。转录因子（TF）分析表明，ANPEP与12个TFs互作，CYP1B1与6个TFs互作，其中SPI1、STAT3、CEBPB、ESR1和FOXA1为两者共同的TFs。转录因子结合位点（TFBS）预测发现ANPEP启动子区有170个潜在结合TFs，CYP1B1有201个。GeneMANIA分析显示功能相似基因主要富集于“氨肽酶活性”、“外肽酶活性”等功能。药物预测共鉴定出68种靶向ANPEP的药物（如白杨素Chrysin）和32种靶向CYP1B1的药物（如秋水仙碱Colchicine）。
四、讨论与结论
讨论总结： 研究人员指出，CYP1B1编码的单加氧酶参与代谢、炎症、血管生成等过程，其异常上调可能通过产生活性氧（ROS）、激活NF-κB通路、促进色氨酸犬尿氨酸代谢等途径参与T2DM和ACI的发生发展。分子对接提示HMCA（由Eubacterium ramulus和Flavonifractor plautii产生）可能与CYP1B1存在潜在相互作用，但具体效应需实验验证。ANPEP参与蛋白质加工、氨基酸代谢及肾素-血管紧张素系统（RAS）和谷胱甘肽（GSH）代谢通路，可能在T2DM中通过维持GSH平衡发挥抗氧化保护作用，但其在ACI中的角色尚不明确。研究还讨论了Bifidobacterium-蛋氨酸-ANPEP轴的潜在关联，但强调其为待验证的假设模型而非确证的保护性级联反应。此外，共同上游TFs（如STAT3、ESR1）的预测为转录调控机制提供了线索，而白杨素和秋水仙碱等候选药物的预测为后续转化研究奠定了基础。研究人员也指出了多项局限性，包括临床16S队列样本量小且缺乏健康对照和单纯T2DM对照、人群和组织的异质性、非T2DM特异性GWAS终点的使用、以及分析结果的假设生成性质，强调未来需结合多中心多族群队列、靶向代谢组学、酶活测定和动物模型进行验证。
结论翻译： 总之，本研究结合生物信息学分析和16S rRNA测序，将核心基因和关键肠道微生物类群鉴定为T2DM相关遗传易感性与ACI的候选分子特征，并进一步表征了它们与疾病进展相关的多种生物学功能和信号通路。本工作为解析T2DM并发ACI的分子机制建立了系统的生物信息学基础。尽管如此，仍需额外的实验研究来验证所提出的调控关系在疾病起始和进展中的作用。预测的ANPEP/CYP1B1相关网络和候选药物需要酶活测定、靶向代谢组学和体内验证。这些发现应被解释为支持T2DM相关遗传易感性，而非严格的T2DM特异性因果关系。
