运动如何通过肠道菌群调节2型糖尿病患者的胰岛素敏感性？//田浩冬，等
《上海体育学院学报》2023年第4期
论文导读
肠道菌群被证明是2型糖尿病（T2DM）发生与发展的重要因素，且在运动改善T2DM中发挥重要中介作用，本文假设：T2DM患者的运动响应可能与其肠道微生态存在密切联系，宿主肠道菌群的组成关系及其代谢特征可能不仅影响运动响应的出现，还在不同响应患者的运动过程中发挥重要作用。本文基于微生物组学和粪代谢组学，通过分析两类T2DM患者的肠道菌群及其代谢特征的差异，探讨了T2DM患者运动响应与其肠道菌群代谢间的关系。受试者为16名中老年T2DM患者，对其进行为期6周，每周3次的高强度间歇运动干预。采用HOMA-IR指数作为患者运动响应的判别依据，响应阈值为2倍HOMA-IR典型误差，将患者分为两组（R组与NR组）。收集受试者粪便样本用于组学检测分析，结果发现：①干预后R组α多样性出现升高，显著高于NR组。②R组菌群关系更复杂、稳定，其关系网络特征在干预后进一步增强。③2组优势菌群差异显著，R组主要包含克里斯滕森菌、青春双歧杆菌等，NR组主要包含韦荣氏球菌、 大芬戈尔德菌等；组内的显著差异仅在R组中出现，主要包含瘤胃球菌、芽孢杆菌等。④R组优势菌群的关联代谢物主要包含花生四烯酸、正缬氨酸等；NR组优势菌群的关联代谢物主要包含L-谷氨酸、2-乙基乙酸等。得出结论：T2DM患者的运动响应与其肠道微生态联系密切。运动响应者肠道微生态更为稳定，包含更多与糖脂代谢、胰岛素代谢、肠道炎症改善正相关的优势菌群，这些菌群的丰度变化可能是运动响应的重要原因。
第一作者简介
田浩冬，男，四川江油人，西南大学硕士研究生。主要研究方向为慢性代谢性疾病的运动康复、运动促进健康的代谢机制。
通信作者简介
彭莉，重庆人，博士，西南大学体育学院教授，博士生导师。主要研究方向为运动促进健康。主持国家级、省部级科研项目10余项，以第一作者或通信作者身份发表学术论文60余篇，出版学术专著4部、译著1部。
本文引用格式
田浩冬，黄丽，刘昊为，等.T2DM患者运动响应与肠道菌群代谢：关系探究与机制推断[J].上海体育学院学报,2023, 47(4):12-23
T2DM患者运动响应与肠道菌群代谢：关系探究与机制推断
人类进行阶段性运动后机体所出现的变化被称为运动响应。在传统意义上，运动响应被用于反映运动能力指标和人体测量学指标改变的异质性。目前，伴随诸多疾病防治中运动疗法的普及，运动疗效的异质性也开始显现，运动疗法对核心临床指标的改善被称为疾病人群的运动响应。
2型糖尿病（T2DM）为临床最常见的慢性疾病之一，患者基数大、并发症多，是危害全球公共卫生安全的重大问题。对于T2DM患者，血糖紊乱仅为发病表象，胰岛素代谢紊乱才是致病根源，因此反映胰岛素敏感性的胰岛素抵抗指数（HOMA-IR）常被作为判断病情发展的重要指标。运动疗法则可通过改善胰岛β细胞功能、增加其数量来促进胰岛素分泌和提升胰岛素敏感性。现有证据表明，运动疗法对胰岛素敏感性的改善存在异质性，6%～69%的糖尿病或肥胖症患者接受运动干预后胰岛素敏感性并未产生有利变化，被称为运动不响应患者。因此，T2DM运动疗法的疗效问题仍是当前研究的重要靶点。
目前研究大多聚焦运动干预方式的疗效差异，对患者自身的运动响应却鲜有关注，微生物组学技术的发展为其提供了全新的解释途径。诸多研究证明，肠道菌群不仅在运动干预T2DM中发挥了重要中介作用，还独立影响着T2DM的产生与发展。宿主肠道菌群的传统干预以特定饮食或药物补给为主，但外服的益生菌无法在肠道重新定植，疗效短暂。相比之下，运动则能在增强体质的同时持续调节宿主的肠道微生态，并以此改善胰岛细胞功能和炎症反应，帮助控制T2DM。然而，个体间肠道菌群的组成受饮食、生活习惯、体成分和基因等诸多因素的影响而存在普遍差异，这或许会影响T2DM患者机体对运动的应答，影响运动疗效，导致异质性的出现。本文比较了同种干预下运动响应者与不响应者肠道微生物的种群特征、关系网络、相关代谢物等要素，试图探究并推断T2DM患者的运动响应与其肠道微生态之间的关系。
1
研究对象与方法
1.1 研究对象
1.1.1 受试者纳入与剔除标准
受试者为16名T2DM患者。所有受试者均提前签署知情同意书，完成基本信息、运动风险筛查问卷以及《国际体力活动量表短问卷》，确保其符合纳入标准。同时制定剔除标准，确保实验过程的安全性和实验结果的可靠性。
纳入标准：①符合中华医学会糖尿病分会T2DM诊断标准；②1年内病情稳定且未发生糖尿病急性代谢紊乱；③近1个月内血糖稳定且饮食有所控制；④每周运动不高于2次，单次运动时间不超过15min；⑤每周吸烟不多于10支，近1个月未饮酒；⑥通过运动风险筛查。
剔除标准：①干预期间出现严重不良反应或意外，如晕厥、骨折等，无法继续训练；②干预期间病情出现加重或其他突发事件，需紧急就医治疗；③未按照研究制定的方案接受干预、检查、随访，或由于资料过少不能满足疗效和安全性评价；④研究期间因其他个人原因主动退出。
1.1.2 受试者基本情况
根据受试者的基本问卷与血液检查结果，统计受试者的基本信息，包含性别、年龄、患病年限、空腹血糖（FPG）、空腹胰岛素（FINS）、HOMA-IR以及服药种类。根据运动干预后HOMA-IR变化对受试者进行分组：根据3周前胰岛素抵抗的重复测试结果计算典型误差（Typical Error，TE=SD_diff/√2），根据公式HOMA"-" IR=(FPG×FINS)/22.5对受试者干预前后HOMA-IR进行计算，以HOMA-IR降低为响应方向，若受试者干预前后HOMA-IR降低量大于2倍TE则将其归入R组，反之则将其归入NR组。在本文中，SD_diff（HOMA-IR）=0.28，TE=0.2，参照HOMA-IR降低量获得R组（n=9）和NR组（n=7），2组受试者基本信息如表1所示。
1.2 研究方法
1.2.1 实验法
所有受试者均接受6周高强度间歇训练（HIIT）。实验流程主要包含适应阶段与正式干预阶段，在适应阶段结束后第3天清晨进行前测样本收集，包含血液样本与粪便样本。随后进行时长6周、每周3次、每次约35min的HIIT干预，干预结束后第3天清晨，再次收集受试者血液样本与粪便样本，测定与前测相同的指标。实验期间受试者保持其原有药物服用及饮食起居等生活方式不变，并进行定期线上随访。详细实验流程见图1。
有研究认为，运动强度越大越有利于病情改善。HIIT强度较高、干预效果较好，同时更节省运动时间。中长期的HIIT干预在改善T2DM患者血糖控制、心肺功能和降低血管风险方面是安全有效的，效果优于中强度运动，可作为一种有效的T2DM干预手段。《WHO关于身体活动和久坐行为的指南》指出，T2DM患者应尽可能地参与持续性、节律性运动，且以跑跳类运动为宜，因此本文以跳绳为主要运动方式。考虑到T2DM患者群体的特殊性且受试者年龄普遍较大，为规避因牵绊导致的训练风险，采用无绳跳绳的方式提升训练安全性。目前针对T2DM患者的HIIT方案总时间大多在8～28min，训练组数一般在4～8组。通常以心率作为强度监测手段，训练强度通常为70%～95%HRmax，间歇强度通常为50%～60%HRmax，单组训练时间∶间歇时间通常为1∶1、4∶3、2∶1、5∶1等。
考虑到患者训练过程中体能的增强，为保证训练强度基本一致，本文训练采用渐进式负荷（表2）。每次干预主要包含8～10min的准备练习（活动下肢关节，原地踏步热身）；16～18min的正式练习（8组跳绳练习，单组时长90s，恢复30～45s，恢复形式随实验阶段变化），全程佩戴心率表监控强度，并在每组训练及间歇结束后询问RPE[跳绳心率控制在70%～80% HRmax（HRmax= 220−年龄），RPE为16～19，恢复期心率控制在50%～60%HRmax，RPE为11～15]；3～5min的放松练习（拉伸主要运动肌肉）组成。
1.2.2 指标获取及检测方法
（1）糖代谢指标
糖代谢测试要求受试者提前禁食8h以上，于次日清晨空腹采集肘部静脉血3～5mL置于无抗凝剂的试管中，采用葡萄糖氧化酶法检测FPG，采用放射免疫分析法（RIA）检测FINS，并计算HOMA-IR。
（2）肠道菌群及粪便代谢指标
在体检前一天向受试者发放独立密封的粪便采集器、一次性便盆（直径略小于家用马桶）以及无菌保鲜膜，告知受试者于次日清晨排泄前对外阴及肛门区域进行清理，防止尿液或体液污染粪便样本。嘱咐受试者在排便前将无菌保鲜膜贴至无菌便盆表面，随后将其放置于便池或马桶上，进行排便。排便结束后将粪便采集器内小勺插入粪便中段，从新鲜晨便中段内部挖取3～5g样本放入粪便采集器。随后立即转交实验人员，实验人员将预先准备好的包含受试者及其粪便样本对应信息的标签贴至采集器表面，并用生物冰袋包裹，迅速转移至−80℃的冰箱冻存，随后尽快送检。若受试者未能按时采集粪便样本，或粪便样本受到腹泻、肛裂、痔疮破裂等意外情况的污染，则视为不合格样本，不予纳入。
肠道菌群DNA测序流程如下。①DNA 抽提和PCR扩增。根据E.Z.N.A.® soil DNA Kit（Omega Bio-tek, Norcross, GA, 美国）说明书进行微生物群落总 DNA抽提，使用1%的琼脂糖凝胶电泳检测DNA的提取质量，使用NanoDrop2000测定DNA浓度和纯度；使用338F（5′-ACTCCTACGGGAGGCAGCAG-3′）和806R（5′-GGACTACHVGGGTWTCTAAT-3′）对16S rRNA基因V3-V4可变区进行PCR扩增。扩增程序为95℃预变性3min，27个循环（95℃变性30s，55℃退火30s，72℃延伸30s），然后72℃稳定延伸10min，最后在4℃进行保存（PCR仪为ABI GeneAmp® 9700型）。PCR反应体系为5×TransStart FastPfu 缓冲液4μL，2.5mMdNTPs 2μL，上游引物（5uM） 0.8μL，下游引物（5uM）0.8μL，TransStart FastPfu DNA聚合酶0.4μL，模板DNA 10ng，补足至20μL，每个样本重复3次。②Illumina Miseq测序。将同一样本的PCR产物混合后使用2%琼脂糖凝胶回收PCR产物，利用AxyPrep DNA Gel Extraction Kit（Axygen Biosciences, Union City, CA, 美国）进行回收产物纯化，2%琼脂糖凝胶电泳检测，并用Quantus™ Fluorometer（Promega, 美国）对回收产物进行检测定量。使用NEXTflexTM Rapid DNA-Seq Kit（Bioo Scientific，美国）进行建库：接头链接；使用磁珠筛选去除接头自连片段；利用PCR扩增进行文库模板的富集；磁珠回收PCR产物得到最终的文库。利用Illumina公司的Miseq PE300平台进行测序（上海美吉生物医药科技有限公司），原始数据上传至NCBI SRA数据库。③数据处理。使用fastp软件（version 0.20.0）对原始测序序列进行质控，使用FLASH软件（version 1.2.7）进行拼接：过滤reads尾部质量值20以下的碱基，设置50bp的窗口，如果窗口内的平均质量值低于20，从窗口开始截去后端碱基，过滤质控后50bp以下的reads，去除含N碱基的reads；根据PE reads之间的overlap关系，将成对reads拼接（merge）成1条序列，最小overlap长度为10bp；拼接序列的overlap区允许的最大错配比率为0.2，筛选不符合序列；根据序列首尾两端的barcode和引物区分样品，调整序列方向，barcode允许的错配数为0，最大引物错配数为2。
基于默认参数，使用Qiime2流程中的DADA2插件对质控拼接后的优化序列进行降噪处理。DADA2降噪处理后的序列通常被称为ASVs（扩增子序列变体）。为了降低测序深度对后续α多样性和β多样性等分析结果的影响，将每个样本降噪后的序列数抽平到4000条，抽平后每个样本的平均覆盖度（Good's coverage）依然能达到97.90%。基于Sliva 16S rRNA数据库（v138），使用Qiime2中的Naive Bayes（或Vsearch或Blast）分类器对ASVs进行物种分类学分析。
粪便代谢组学分析如下。① 采集受试者粪便样本3～5g置于−80℃的冰箱，待样本收集完全后进行样本前处理，处理流程：精确称取50mg样本到2mL离心管里，加入0.5mL甲醇-水溶液（体积比CH3OH∶H2O =4∶1，含0.02mg/mL的内标L-2-氯-苯丙氨酸）；加入1颗钢珠，在−20 ℃研磨机中研磨（50Hz，3min）；加入200μL的氯仿，在研磨机中研磨（50Hz，3min），超声提取30min，20 ℃静置30min；低温离15min（13000rcf，4℃），取上清液装入玻璃衍生瓶中，氮气吹干；向玻璃衍生瓶中加入80μL的甲氧胺盐酸吡啶溶液（15mg/mL），涡旋振荡2min后，于振荡培养箱中37 ℃肟化反应90min；取出后再加入80μL的BSTFA（含1%TMCS）衍生试剂，涡旋振荡2min后，于70℃反应60min；取出样本，室温放置30min。②完成基本处理后采用气相色谱质谱联用仪进行GC/MS分析、代谢组学分析，处理色谱条件：衍生化后样本用分流模式注入GC/MS系统进行分析，进样量1μL，分流比10∶1。样品经DB-5MS 毛细管柱（40m×0.25mm×0.25μm，Agilent 122-5532G）分离后进入质谱检测。进样口温度260℃，载气为高纯氦气，载气流速1mL/min，隔垫吹扫流速3mL/min，溶剂延迟5.5min。升温程序：初始温度60℃，平衡0.5min，然后以8℃/min的速度升至310℃并维持6min；处理质谱条件为：电子轰击离子源（EI），传输线温度为310℃，离子源温230℃，四极杆温度150℃，电子能70eV。扫描方式为全扫描模式（SCAN），质量扫描范围：50～500m/z，扫描频率为3.2scan/s。
1.2.3 数理统计
采用IBM SPSS 26.0 计算典型误差，对运动干预前后糖代谢指标变化进行配对样本t检验。对粪便样本进行16sRNA全基因测序，采用MassHunter workstation Quantitative Analysis（v10.0.707.0）软件进行峰提取、对齐等数据预处理操作，最终得到代谢物鉴定结果及数据矩阵。采用Wilcoxon rank-sum检验分析不同响应组的肠道菌群α多样性，采用Bray_curtis距离算法对有无响应组的肠道菌群β多样性进行PCoA分析；限定LDA>2.5，P<0.05对各组肠道菌群组成进行LEfSe多级物种差异判别分析；结合Spearman相关系数绘制菌群间单因素相关性网络；提取聚类前50的代谢物与各组丰度前50的菌属进行关联分析，随后对各组优势菌属关联代谢物进行ROC分析，提取预测能力前12的关键代谢物。
2
研究结果
2.1 R组与NR组的糖代谢变化比较
R组与NR组接受干预后糖代谢变化结果（表3）显示，运动干预后2组受试者糖代谢变化存在较大差异：R组FPG、FINS、HOMA-IR均出现了显著降低；NR组FINS和HOMA-IR均出现了显著升高，FPG无显著变化。
2.2 R组与NR组的肠道菌群特征比较
2.2.1 R组与NR组肠道菌群多样性比较
对提取的粪便样本进行全基因组测序，获得了平均每个样本1051118148的碱基对含量，对其进行菌群多样性分析。α多样性分析纳入反映群落内部丰度的Chao指数和多样性的Shannon指数，结果（图2）显示：2组前测样本α多样性均无显著差异，但R组后测样本菌群Shannon指数（3.30±0.68 vs 2.62±0.67，P=0.034）与Chao指数（270.56±75.04 vs 180.46±54.66，P=0.026）均显著高于NR组的；运动干预后R组受试者α多样性比干预前出现了一定提高，NR组则出现了一定下降。采用PCoA分析对2组受试者肠道菌群表型的总体差异进行分析，结果（图3）显示干预前后的组内及组间比较均不存在显著差异（R组干预前后差异，R=−0.09，P=0.99；NR组干预前后差异，R=−0.09，P=0.79；2组前测差异，R=0.05，P=0.29；2组后测差异，R=0.10，P=0.17）。
2.2.2 R组与NR组肠道菌群共存关系网络比较
在属水平下提取丰度前30的菌属，绘制关系网络对比2组受试者肠道菌群的共存关系，结果（图4、图5）显示：R组受试者运动干预后菌群网络传递性增强（0.48→0.52），网络直径（10→7）和最短路径长度（3.61→2.84）均缩短，且干预后围绕厚壁菌门下布劳特氏菌属（g_Blautia_hallii_group）、克里斯滕森菌属（g_Christensenellaceae_R-7_group）、毛螺菌属（g_Lachnoclostridium)、瘤胃球菌属（g_Ruminococcus_torques_group）等出现了更多强相关线段；干预前后的NR组菌群都存在非连通部分，网络直径与最短路径长度均无有效值，且运动干预后的菌群网络传递性降低（0.51→0.44）。这表明相较于NR组，运动干预更显著地增强了R组菌群间的共存关系。
2.2.3 R组与NR组肠道菌群差异物种比较
上述分析显示2组受试者菌群关系的差异伴随着物种差异，因此采用LEfSe多级物种差异判别分析筛选2组具有显著差异的菌群，通过LDA估算各菌群丰度对差异影响的大小。LEfSe分析中LDA的选取一般不低于2，本文选取LDA=2.5的判别阈值，P<0.05的显著性水平对2组菌群从门水平到种水平提取优势差异菌群。结果（图6）显示，运动干预前的样本共定位到28个优势菌群，其中，23种在R组出现了显著富集， 5种在NR组出现了显著富集。干预后共定位到33个优势菌群，其中，31种在R组出现了显著富集，2种在NR组出现了显著富集。组间主要差异菌群见表4。
分别对干预前后2组肠道菌群进行差异菌的组内比较，结果（图7）显示，NR组干预前后不存在差异菌群，R组干预后36种菌群丰度出现了显著上升。组内主要菌群差异见表5。
2.3 R组与NR组的菌群代谢物分析
2.3.1 肠道菌群与代谢物关联分析
分别对2组干预前后的菌群及粪便代谢物进行层级聚类，取丰度前50的菌属变化与聚类前50的代谢物变化进行关联分析，观测差异菌群诱导的代谢物差异。R组菌群代谢物关联分析结果（图8）显示，与R组优势菌属呈显著正相关的代谢物主要包含与高血糖和胰岛素抵抗密切相关的花生四烯酸（arachidonic acid）、反式13-十八烯酸（trans-13-octadecenoic）、正缬氨酸（norvaline）、L-缬氨酸（L-valine）、2-羟基喹啉（2-hydroxyquinoline）、草酸（oxalic acid）及天然抗氧化剂对伞花烃（p-cymene）等；呈显著负相关的代谢物主要包含与糖尿病视网膜并发症有关的2-哌啶酮（2-piperidone），与脂肪异常代谢有关的5-氨基戊酸（5-Aminovaleric Acid）、肉豆蔻酸（myristic acid），与心血管疾病风险和胰岛素抵抗呈正相关的反油酸（elaidic acid）、糖尿肾病潜在标志物1-十六烷醇（l-hexadecanol）及抗氧化剂柠康酸（citraconic acid）。NR组菌群代谢物关联分析结果（图9）显示，与NR组优势菌属呈显著正相关的代谢物主要包含1-辛烯-3-醇（1-Octen-3-ol）、L-谷氨酸（L-glutamic acid）；呈显著负相关的代谢物为肝脏线粒体增加的重要生物标志物2-乙基己酸（2-ethylhexanoic acid）、2-羟基喹啉及黄毒素（xanthotoxin）等。
2.3.2 相关代谢物ROC分析
采用监督学习对样本进行线性判别和分组建模，对伴随各组出现变化的代谢物进行ROC分析，通过AUC（Area Under the Curve）衡量该代谢物的变化对运动响应的预测能力，评估其对受试者运动响应的关键程度，一般认为AUC大于0.7表示代谢物预测能力较为理想。结果（图10）显示，共12种代谢物体现出了一定的预测能力，其中：5种有较为理想（AUC>0.7）的预测能力，包括5-氨基戊酸（AUC=0.81）、正缬氨酸（AUC=0.75）、柠康酸（AUC=0.73）、肉豆蔻酸（AUC=0.71）、2-羟基喹啉（AUC=0.71）；剩余7种代谢物主要包含黄毒素、L-谷氨酸、反油酸、2-乙基己酸、花生四烯酸、1-辛烯-3-醇、反13-十八烯酸，其AUC区间为0.60～0.70。
3
分析与讨论
3.1 菌群多样性及共存关系与运动响应的关系
目前国内探讨T2DM患者运动响应与肠道菌群关系的研究十分少见，本文可在一定程度上填补此空白。本文首先以菌群多样性为切入点，探究运动响应与整体群落特征间的关系，结果显示受试者的β多样性均未在运动后出现显著改变，R组菌群的α多样性随胰岛素抵抗的降低出现了升高，这与其他研究一致，而NR组α多样性却出现了降低。尽管干预后2组的α多样性出现了显著差异，但干预前并未出现，这提示α多样性并非导致运动响应差异的直接诱因。
有研究报道，代谢疾病与宿主肠道微生态的稳定性密切相关，为此本文从菌群共存关系出发进一步探究。对比干预前后2组的菌群关系网络发现，R组的肠道菌群存在更稳固复杂的共存关系。干预后R组在布劳特氏菌属、克里斯滕森菌属、毛螺菌属等有益菌的外围建立了更强的共存关系，而NR组关系网络的传递性在干预后出现下降。这与2组受试者胰岛素的变化一致，提示肠道内有益菌群间的共存关系可能是影响运动响应的重要因素，运动亦能反过来强化该共存关系。
3.2 肠道菌群组成与运动响应的关系
2组患者肠道菌群组成的差异在前后测中均十分明显。尚未干预前，运动响应者的肠道便出现了诸多益生菌及潜在益生菌的显著富集，这些菌群作用十分广泛，且存在明显共性。以青春双歧杆菌、克里斯滕森菌、粪便罗斯氏菌及毛螺旋杆菌为例的优势菌群均对调节血糖、降低胰岛素抵抗具有良好助力。从作用机制上看：上述菌群能从多方面发挥作用，青春双歧杆菌能有效保护胰岛β细胞超微结构，维系其功能；克里斯滕森菌能帮助建立健康的葡萄糖代谢循环；后两者则主要通过直接或间接产生短链脂肪酸（SCFA），改善机体炎症。在干预后，运动响应者的肠道菌群的变化主要体现在两方面：首先，原有菌群优势的延续，克里斯滕森菌、粪便罗斯氏菌、毛螺旋杆菌在后测中仍旧显著富集；其次，经过运动干预，普拉梭菌和芽孢杆菌等潜在益生菌出现了显著富集。普拉梭菌是人体丰度最高也是最重要的细菌之一，同时是丁酸的重要生产者，具备良好的抗炎能力，且普拉梭菌还可诱导CD4+/CD8α+调节性T细胞在结肠固有层中发挥保护作用，抵御诸多肠道病原体的侵害；芽孢杆菌则能抑制肠道致病菌（如金黄色葡萄球菌）的定植、减少炎症诱导因子诱导的一氧化氮合酶（iNOS）。综上，推测运动具有保护和改善运动响应患者肠道微生态的潜力。
相比之下，无论是干预前还是干预后，运动不响应者优势菌群均较为匮乏，甚至出现了韦荣氏球菌和小杆菌等有害菌群的显著富集。二者与肠道炎症、血糖受损、肥胖形成、肠胃疾病发生等多个病理过程关系密切。NR组的组内分析也并未检测到前后富集的显著差异，提示运动可能无法改善不响应者的肠道微生态。综上，运动响应者与不响应者在原生肠道菌群方面存在显著差异，尤其体现在相关益生菌与潜在益生菌的富集方面。运动响应患者肠道中存在多种利于T2DM改善的菌群，这可能是其对运动产生有利应答的重要原因，而对于运动不响应者，运动无法有效改善其肠道菌群组成。特定食物和药物的补给或许是解决方案之一。例如：可发酵型膳食纤维（如植物种子类食物）的摄入或在食物中添加阿卡波糖等淀粉抑制剂（可将升血糖淀粉转化为抗性淀粉）可增加肠道内毛螺旋杆菌的丰度 ；采用芽孢杆菌发酵的绿茶可减少脂肪消化吸收、增加能耗、降低机体厚壁菌门比例以调节脂代谢；青春双歧杆菌与小檗碱联用可改善患者肠道结构与功能，调节菌群组成以叠加降糖效果。
3.3 优势菌群相关代谢物与有无运动响应的关系
菌群与代谢物相互作用，在各种代谢疾病改善中扮演重要角色，在ROC分析中，5-氨基戊酸对运动响应最为关键（AUC=0.81），且与R组中多个优势菌属的富集存在显著正相关。多项基于高脂小鼠的实验证实了其对糖代谢改善的良好效果，指出其对高血糖的缓解主要依赖于对C2C12细胞摄取和转运葡萄糖的促进作用。同样较为关键的正缬氨酸（AUC=0.75）、肉豆蔻酸（AUC=0.71）则均被证实对脂代谢的改善效果显著，它们仅与运动响应组的优势菌属UCG-002菌属和克里斯滕森菌属显著正相关。正缬氨酸在2项针对高血糖大鼠的实验中显著降低了血液总胆固醇、低密度脂蛋白及氧化应激标志物丙二醛；肉豆蔻酸能促进褐色脂肪活化和米色脂肪形成，帮助耗能、平衡体脂。花生四烯酸（AUC=0.65）也与UCG-002和NK4A214两个优势菌属显著正相关，作为一种常见的不饱和长链脂肪酸，它有助于增强胰岛素敏感性，响应基础葡萄糖的浓度变化。上述4种代谢物的增加，提示这些优势菌属可能促进了运动响应组在胰岛素敏感性、血糖控制、有益脂肪形成等方面的应答能力，印证了菌群对运动响应的重要作用。
与R组相比，NR组并未呈现完全相反的代谢特征，如NR组出现了1-辛烯-3-醇（AUC=0.65）、L-谷氨酸（AUC=0.68）等有益代谢物。1-辛烯-3-醇对金黄色葡萄球菌、大肠杆菌、三星镰刀菌等有害菌群具有较强的抗菌活性，可有效抑制其生长；L-谷氨酸可通过产生谷胱甘肽抑制氧化应激的产生，降低对胰岛β细胞的损伤、增强外周胰岛素敏感性。尽管1-辛烯-3-醇与R组瘤胃球菌属、毛螺菌属等优势菌属显著正相关，但其并未在运动响应组出现富集。由此可见，2组患者的代谢特征并未出现两极分化，差异明显弱于菌群特征。值得注意的是，肠道菌群和粪便代谢物的关系并非简明单一，个体生活习性的差异和检测手段有限都是重要局限，现有证据仅指出有益代谢间的程度差异可能是运动响应的重要影响因素。
4
结 论
T2DM患者的运动响应与肠道菌群及其代谢具有密切联系。青春双歧杆菌、克里斯滕森菌、毛螺旋杆菌等有利于糖脂代谢、胰岛素代谢、肠道抗炎的益生菌或潜在益生菌在运动响应者的肠道中显著富集，并建立了稳固的共存关系。这或许是T2DM患者拥有运动响应的重要条件。
运动也能延续或增强运动响应者的肠道微生态优势，主要体现为有益菌群代谢的增强与共存关系的巩固，但无法改善不响应者的肠道微生态，甚至降低了其菌群多样性。因此，对不具有运动响应的患者采用饮食调控和药物干预或许是更好的选择，而是否应避免运动或事先改善其肠道微生态来诱导运动响应值得进一步研究证实。
本文刊于《上海体育学院学报》2023年第4期。为方便阅读，此处删去原文注释，如果其他媒体或机构转载，请标明文章出处。
上海体育学院学报
长按扫码关注我们