《International Journal of Molecular Sciences》：Longitudinal Gut Microbiota Changes in Women with Breast Cancer-Related Lymphedema Undergoing a Lifestyle Intervention
编辑推荐：
：乳腺癌相关淋巴水肿（BCRL）与肥胖密切相关，但该人群中肠道菌群的特征仍未得到充分阐明。研究人员对48名超重或肥胖且患有BCRL的女性进行了肠道菌群谱的纵向评估，这些受试者被随机分配至多模式生活方式干预组（结合限能量地中海式饮食与运动，n = 28）或对照组
：乳腺癌相关淋巴水肿（BCRL）与肥胖密切相关，但该人群中肠道菌群的特征仍未得到充分阐明。研究人员对48名超重或肥胖且患有BCRL的女性进行了肠道菌群谱的纵向评估，这些受试者被随机分配至多模式生活方式干预组（结合限能量地中海式饮食与运动，n = 28）或对照组（n = 20）。在基线、3个月和6个月时采集粪便样本（共144份），并进行16S rRNA基因测序。干预组在体重、身体质量指数（BMI）和脂肪量方面随时间出现更大程度的下降（组别×时间交互作用p < 0.001、p = 0.001和p = 0.002）。基于Bray–Curtis相异度的α多样性和总体微生物群落组成未观察到显著的组别×时间交互作用（PERMANOVA R2 = 0.006，F = 0.459，p = 0.324）。探索性DESeq2差异丰度分析识别出经FDR校正后显著的类群特异性差异，包括干预组中Barnesiellaceae和Desulfovibrionaceae丰度较高，对照组中Prevotella_7和Lachnospiraceae_UCG_010丰度较高（FDR校正后p < 0.05）。这些DESeq2结果未在FDR校正后的混合效应敏感性分析中得到证实。Prevotella_7在基线和6个月时在对照组中丰度更高。具有推断优先权的受试者分块组别×时间模型在基线至3个月的交互作用检验中得出Izemoplasmatales和Frisingicoccus的FDR显著结果，尽管相应的收缩效应估计值接近零。基于补充剂使用状态、补充剂依从性和临床相关减重反应的次要假设生成分析识别出与补充剂使用和依从性相关的额外类群特异性模式，包括补充剂使用者中涉及两歧双歧杆菌（Bifidobacterium bifidum）的初步物种水平关联，而按减重反应分组未观察到显著的菌群差异。总体而言，这些发现未提供肠道微生物群落广泛重构的证据，但识别出值得在更大规模纵向研究中采用菌株分辨和功能方法加以验证的探索性类群特异性信号。
论文解读
一、研究背景与目的
乳腺癌相关淋巴水肿（BCRL）是乳腺癌治疗后因淋巴引流受损和组织进行性重塑导致的慢性并发症，以持续炎症、脂肪组织沉积和纤维化为特征，造成肢体肿胀、功能障碍和生活质量下降。超重和肥胖与BCRL发生风险增加密切相关，脂肪组织扩张、慢性低度炎症和淋巴转运受损可能进一步损害淋巴功能；反之，慢性淋巴功能障碍也可能促进局部脂肪堆积和炎症性组织重塑，提示淋巴功能障碍与肥胖之间存在双向关系。肠道菌群作为饮食、脂肪组织、代谢和免疫调节之间的重要界面，其组成和代谢产物通过短链脂肪酸、胆汁酸、脂多糖等途径影响能量稳态、肠道屏障完整性及局部和全身免疫应答。肠道菌群改变与肥胖及其相关代谢紊乱存在关联，而地中海饮食模式和体力活动均可影响肠道菌群组成。然而，此前尚无研究直接利用16S rRNA基因测序表征BCRL女性的肠道菌群特征，肠道菌群在该人群生活方式干预生物学应答中的潜在作用亦未被阐明。基于此，本研究作为一项已发表的随机生活方式干预试验的探索性肠道菌群亚研究，旨在纵向表征BCRL合并超重或肥胖女性的肠道菌群多样性和分类组成，并探讨这些模式是否随干预组别、补充剂使用状态、补充剂依从性及临床相关减重反应而存在差异。
二、研究设计与方法概述
本研究为一项在西班牙马拉加地区医院开展的随机先导试验（ClinicalTrials.gov注册号：NCT04974268）的肠道菌群亚研究，采用开放标签、双臂平行设计，纳入18–79岁、BMI为25–40 kg/m2的BCRL女性。最终分析队列包括48名受试者，其中干预组28人、对照组20人，均需在基线、3个月和6个月提供粪便样本且各样本经质量过滤后获得至少3000条序列读数。干预组接受结构化生活方式方案，包括监督下抗阻和有氧运动联合限能量地中海式饮食；干预组受试者可自愿每日服用一份含初榨橄榄油、亮氨酸、omega-3脂肪酸和热灭活形式的动物双歧杆菌乳亚种CECT 8145（BPL1）的口服营养补充剂。主要技术方法包括：粪便样本DNA提取后采用Ion 16S宏基因组试剂盒扩增16S rRNA基因多个可变区，在Ion S5平台测序；QIIME 2流程经DADA2去噪生成扩增子序列变异（ASV），采用Greengenes2数据库进行分类学注释；α多样性以Chao1、Shannon、Simpson和逆Simpson指数评估，β多样性基于Bray–Curtis相异度经主坐标分析（PCoA）和PERMANOVA评估；差异丰度分析采用DESeq2，辅以混合效应模型敏感性分析。
三、主要研究结果
3.1 研究人群基线特征
干预组与对照组在基线时总体可比，唯一显著差异为对照组血脂异常的患病率高于干预组（p = 0.007）。其余人口学、生活方式、癌症相关、合并症、人体测量、身体成分、淋巴水肿相关、饮食和生化变量均无显著组间差异。
3.2 临床结局变化
与对照组相比，干预组在3个月和6个月时体重下降幅度更大（p = 0.022和p = 0.002），组别×时间交互作用显著（p < 0.001）；BMI和脂肪量亦呈现类似模式，均存在显著组别×时间交互作用（p = 0.001和p = 0.002）。腰围、去脂体重和相位角无显著组间差异或交互作用。所有淋巴水肿相关结局均无显著组间差异。能量摄入在3个月时干预组降幅更大（p = 0.045），存在显著组别×时间交互作用（p = 0.040）；地中海饮食依从性变化两组相当。
3.3 肠道菌群分类组成与多样性
经质量过滤和抽平至每样本3000条序列后，共鉴定出255个属水平或最近解析分类水平的细菌类群，平均Good's覆盖率为99%。所有α多样性指数均未显示显著的组别×时间交互作用（所有p ≥ 0.430）。β多样性方面，PCoA显示各组-时间组合间存在大量重叠，PERMANOVA未检测到显著的组别×时间效应（R2 = 0.006，F = 0.459，p = 0.324），PERMDISP亦无显著差异（p = 0.223）。
3.4 按研究组别的差异丰度类群
总体DESeq2分析（校正采样时间）识别出4个经FDR校正后显著差异的类群：Barnesiellaceae（收缩log2FoldChange = 2.248，FDR校正p < 0.001）和Desulfovibrionaceae（收缩log2FoldChange = 1.443，FDR校正p = 0.019）在干预组中富集；Prevotella_7（收缩log2FoldChange = ?1.946，FDR校正p < 0.001）和Lachnospiraceae_UCG_010（收缩log2FoldChange = ?1.419，FDR校正p = 0.014）在对照组中富集。混合效应敏感性分析虽保留了方向但均未在FDR校正后保持显著。受试者分块纵向分析中，基线至3个月的组别×时间交互作用对Izemoplasmatales（收缩交互log2FoldChange = 0.042，FDR校正p < 0.001）和Frisingicoccus（收缩交互log2FoldChange = ?0.127，FDR校正p = 0.042）达到FDR显著，但收缩后效应估计值接近零。时间分层分析显示Prevotella_7在基线和6个月时均在对照组中富集，6个月时Escherichia-Shigella亦在对照组中富集（收缩log2FoldChange = ?0.169，FDR校正p = 0.011）。采用≥20%患病率阈值的敏感性分析中，Barnesiellaceae、Desulfovibrionaceae和Lachnospiraceae_UCG_010仍保持显著，Prevotella_7因患病率低于阈值被排除。校正基线血脂异常的事后敏感性分析中，4个类群均保持FDR显著且方向一致。干预组内纵向分析未发现任何FDR显著的类群变化。
3.5 按补充剂使用状态的探索性菌群分析
补充剂使用为干预组内的自愿行为，未随机化，故此类比较属次要假设生成分析。PERMANOVA未检测到显著的补充剂状态×时间效应（R2 = 0.009，F = 0.368，p = 0.548）。总体分析中Prevotella_9（收缩log2FoldChange = 3.626，FDR校正p < 0.001）和Klebsiella（收缩log2FoldChange = 0.932，FDR校正p = 0.003）在补充剂使用者中丰度更高，但时间分层分析显示这些差异在基线时已存在，且所有未使用补充剂的7名受试者在基线时均未检出这两个类群。属水平Bifidobacterium未因补充剂状态而显著差异（FDR校正p = 0.740）；物种水平上，两歧双歧杆菌（Bifidobacterium bifidum）在补充剂使用者中丰度更高（收缩log2FoldChange = 1.707，FDR校正p < 0.001），而所施用BPL1菌株所属的动物双歧杆菌（Bifidobacterium animalis）未显示显著差异（FDR校正p = 1.000）。
3.6 按补充剂依从性的探索性菌群分析
依从性比较亦属探索性分析。PERMANOVA未检测到显著的依从性组×时间效应（R2 = 0.012，F = 0.377，p = 0.537）。总体分析中，Rikenellaceae_RC9_gut_group（收缩log2FoldChange = 0.149，FDR校正p < 0.001）和Escherichia-Shigella（收缩log2FoldChange = 2.521，FDR校正p < 0.001）在高依从组中丰度更高；Prevotella_9（收缩log2FoldChange = ?2.080，FDR校正p < 0.001）、Parasutterella（收缩log2FoldChange = ?2.572，FDR校正p = 0.003）和Lachnospiraceae_AC2044_group（收缩log2FoldChange = ?3.235，FDR校正p = 0.015）在低依从组中丰度更高。6个月时Clostridium_sensu_stricto_1在高依从组中丰度更高（收缩log2FoldChange = 0.445，FDR校正p = 0.012）。混合效应敏感性分析均未在FDR校正后保持显著。
3.7 按减重反应的探索性菌群分析
按是否达到≥5%减重分类的PERMANOVA未检测到显著的减重反应×时间效应（R2 = 0.007，F = 0.496，p = 0.099）。总体和分时间点的DESeq2分析均未发现任何FDR显著的类群差异，混合效应敏感性分析亦未识别出显著关联。
四、讨论与结论总结
讨论部分指出，尽管干预组在体重、BMI和脂肪量方面显著下降，但未检测到整体菌群多样性和群落结构的显著变化，这与既往研究显示减重可在α和β多样性变化甚微或不一致的情况下发生的发现一致。探索性差异丰度分析识别出若干类群特异性信号，包括干预组中Barnesiellaceae和Desulfovibrionaceae富集、对照组中Prevotella_7和Lachnospiraceae_UCG_010富集，其中Prevotella_7的差异在基线时已存在，提示可能反映预先存在的菌群配置而非干预诱导效应。补充剂相关分析中Prevotella_9和Klebsiella的差异亦在基线时已存在，限制了因果归因。Bifidobacterium bifidum在补充剂使用者中的富集属探索性发现，因BPL1以热灭活形式施用，预期不会直接定植或增加所施用菌株的丰度。依从性相关分析识别出多个类群差异，但效应量在收缩后大幅衰减，且混合效应敏感性分析未确认这些发现。研究局限性包括样本量有限、补充剂使用非随机化、基线血脂异常组间不平衡、多次比较问题及16S rRNA测序在菌株水平分辨率上的局限。研究结论为：这是首项在BCRL合并超重或肥胖女性中纵向表征肠道菌群的研究，尽管干预带来显著的体重和脂肪量下降，但未伴随可检测的整体菌群多样性或群落结构变化；探索性分析识别出若干类群特异性和时间性模式，包括Izemoplasmatales和Frisingicoccus的组别×时间交互作用检验达到FDR显著但效应量衰减，以及补充剂和依从性相关的多个类群差异；按临床相关减重反应未观察到显著菌群差异。这些发现不提供群落水平广泛重构的证据，但识别出需要独立验证的探索性类群特异性信号，未来需采用菌株分辨测序和功能方法在更大规模纵向研究中加以阐明。
