《Gut Microbes》：Longitudinal modelling of microbiome subcommunities reveals parity-dependent dynamics during pregnancy and postpartum
编辑推荐：
背景：由于高孕前体重指数（pre-pregnancy BMI, pBMI）导致的母体对妊娠适应性失调与母儿健康结局恶化相关，而肠道微生物组在这些适应性中的作用尚不清楚。方法：研究人员从纳入“Be Healthy in Pregnancy”研究（NCT01689
背景：由于高孕前体重指数（pre-pregnancy BMI, pBMI）导致的母体对妊娠适应性失调与母儿健康结局恶化相关，而肠道微生物组在这些适应性中的作用尚不清楚。方法：研究人员从纳入“Be Healthy in Pregnancy”研究（NCT01689961）的孕妇（n=52）处收集了第一、二、三 trimester 及产后6个月的粪便样本，以及其婴儿6月龄时的样本。经16S rRNA基因测序后，研究人员实施了时间对齐隐含狄利克雷分布（time-aligned LDA, TALDA），采用指数衰减权重的时加权采样策略构建时间近端队列，对各队列独立应用LDA，并使用alto R包对齐所得主题（topic）；婴儿样本使用标准LDA分析。结果：研究人员识别出7个不同的亚群落（subcommunity），其中之一代表会阴污染并被排除。剩余亚群落的分类学定义截然不同且在妊娠期保持稳定（平均跨队列稳定性0.899），而个体内亚群落比例存在差异。经产妇（Multiparous）在早期妊娠表现出亚群落比例的较大转变，初产妇（Primiparous）则显示渐进性增加的变化，在妊娠向产后过渡时变化最大。高pBMI与微生物组重塑减少相关，尤其在经产妇中。探索性分析显示，晚孕期母体短链脂肪酸（short-chain fatty acid, SCFA）产生型亚群落丰度与婴儿6月龄时双歧杆菌（Bifidobacterium）主导的亚群落比例呈正相关，而pBMI>25独立与较低婴儿双歧杆菌水平相关。结论：TALDA有效区分了稳定的亚群落定义与动态的亚群落比例，揭示微生物组在维持功能组织的同时，亚群落比例贡献在妊娠及产后进程中发生偏移。这些母体动态可能通过影响婴儿肠道定殖产生代际后果。该工作表明母体因素修饰微生物组轨迹，并支持母体肠道微生物组内存在妊娠史的生态记忆（ecological memory）。
研究背景方面，母体代谢适应对于确保妊娠期间适宜的母体-胎儿-胎盘能量分配至关重要，其失调会增加妊娠并发症风险及母儿长期代谢功能障碍。现有证据表明妊娠诱导的代谢适应受肠道微生物组重塑影响，但关于妊娠史即产次（parity）如何塑造妊娠过程中微生物组动态的研究甚少。肠道 microbiota 具有高个体间变异性，不同微生物类群可占据相同功能生态位（niche），传统基于多样性指标或单一类群差异丰度的分析难以捕捉其功能意义转变，且既往聚类方法如肠型（enterotypes）定义在连续分类变异上强加离散边界备受批评。近期统计建模如隐含狄利克雷分布（Latent Dirichlet Allocation, LDA）主题建模虽能识别微生物亚群落，但多数现有方法将纵向样本视为独立或不可互换的重复测量，未模拟时间动态，限制了表征妊娠等渐进生理变化的能力。鉴于此前研究发现孕前体重指数（pre-pregnancy BMI, pBMI）影响妊娠微生物组重塑且受产次调节，并假设高pBMI关联重塑减少及产次调制该效应，经产妇因既往妊娠生态记忆（ecological memory）呈现更早偏移，且母体微生物偏移关联婴儿定殖，故开展此项研究。该研究发表于《Gut Microbes》。
主要关键技术方法中，研究人员基于“Be Healthy in Pregnancy”（BHIP, NCT01689961）队列，纳入52名孕妇于妊娠第一、二、三 trimester 及产后6个月采集粪便，并收集匹配婴儿6月龄粪便。经16S rRNA基因测序得扩增子序列变体（Amplicon Sequence Variants, ASVs）并以Silva 138注释，剔除非细菌及会阴污染菌属后，开发时间对齐LDA（time-aligned LDA, TALDA）：以指数衰减权重（λ=0.75）构建各时间点为中心的四时加权队列，各拟合K=7的LDA模型，借alto R包最优传输（optimal-transport）法跨时间对齐主题；婴儿样本施以标准LDA。借Jensen–Shannon散度算稳定性，欧氏距离量化个体内转变，线性混合效应模型与层次聚类析轨迹，对数-对数线性回归析母婴关联，敏感性与验证比对标准LDA及不同λ值。
研究结果部分如下。
Participant characteristics：研究人员最终分析52名母体样本（T1 n=45, T2 n=48, T3 n=43, PP n=43），多数阴道分娩（37/52），产次与pBMI组间分娩方式无差，经产妇年龄大于初产妇（p=0.001），产次与pBMI对 gestational weight gain（GWG）及孕周无显著影响。
TALDA identifies seven microbial subcommunities with distinct temporal patterns：TALDA识7亚群落具独特分类特征，SC1强占Bacteroides伴Faecalibacterium等；SC5 Bacteroides富集共占Escherichia/Shigella与Akkermansia；SC6为典型丁酸（butyrate）产生类群Faecalibacterium、Roseburia及Lachnospiraceae主导；SC7共享此类但均布Lachnospiraceae多属；SC4 Pseudobutyrivibrio与Bacteroides共占；SC2分类弥散含Lachnospiraceae、Blautia等；SC3为Ezakiella、Porphyromonas等阴道栖居类群判为会阴污染剔除。计算Jensen-Shannon散度转稳定性得分（1—JS divergence），SC1、SC6、SC7均值>0.9最稳；SC5略低>0.85受Escherichia/Shigella增与Akkermansia减驱导；SC2与SC4早期T1-T2最不稳而后增稳，次生贡献者互反波动而核心不变，显功能冗余（functional redundancy）与TALDA分辨率。
Maternal factors shape subcommunity proportions：Bootstrap置信区间示亚群落比例关联pBMI与产次。SC4在初产妇自T1 0.07 [0.05–0.10]升至T2 0.16 [0.10–0.22]，pBMI≥25者T1低而T2追平，经产妇无时态变。SC1 Bacteroides主导之比例关联产次与pBMI呈反向：初产妇pBMI<25趋高（0.22 [0.18–0.27] vs 0.15 [0.09–0.21]），经产妇pBMI≥25趋高（0.24 [0.18–0.31] vs 0.18 [0.12–0.25]），T3反差最大。SC6丁酸产族在初产妇全时点高于经产妇（0.19 [0.14–0.23] vs 0.10 [0.07–0.13]），pBMI<25组间最显（0.19 vs 0.08）。
Within-participant microbiome dynamics reveal BMI and parity effects：算个体内亚群落比例向量欧氏距离，线性混合效应模型示pBMI影响稳定性：pBMI<25者距高于≥25（F=7.08, p=0.011），此效驱于经产妇（p=0.011）而非初产妇（p=0.34）。T1-T2经产妇（multiparous）偏移大于初产妇（p=0.048），后时无异。PCA与层次聚为三轨迹簇：簇1微偏移多pBMI≥25（15/16）；簇2渐增偏移峰于T3-PP多初产妇；簇3早T1-T2峰偏移多经产妇。轨迹隶属关联产次（χ2=9.15, p=0.010）与pBMI（χ2=7.01, p=0.030）。
Maternal parity and pBMI associate with infant microbiome composition at 6 months：46例婴儿6月龄样以标准LDA识6亚群落：SC1 Bifidobacterium、SC2 Veillonella、SC3 Lachnospiraceae、SC4 Escherichia/Shigella、SC5 Bacteroides、SC6 Clostridiaceae。排益生菌与配方奶后，母乳加辅食者SC6高SC1低，益生菌反相；剖宫产（Cesarean section, CS）无论益生菌近无SC5。婴儿微生组关联母体：pBMI≥25婴SC4高（0.47 [0.27–0.66] vs 0.18 [0.11–0.25]）、SC2低（0.03 [0.004–0.05] vs 0.12 [0.06–0.20]）、SC5趋低；pBMI<25下经产妇婴SC6低（0.003 [0.001–0.007] vs 0.11 [0.01–0.23]），SC1经产妇与pBMI<25趋高。限阴道娩、母乳、无益生菌23对母婴，对数-对数回归校pBMI与产次示母体T3时SC7（SCFA产族Lachnospiraceae等）比例正关婴SC1 Bifidobacterium（β=0.73, p=0.0006, FDR=0.02），无共有优势类群示间接生态链；母体pBMI≥25独立预测婴SC1低（β=-1.31, p=0.018）。
讨论部分总结，研究人员指出妊娠母体代谢适应部分受微生物组偏移介导，TALDA示亚群落分类定义在妊娠产后稳而比例变，妊娠系既有母微生群落再平衡非菌群失调（dysbiotic disruption）。高pBMI关联少重塑尤经产妇，最大重塑时机关产次异：经产妇早T1-T2，初产妇峰于产后T3-PP，示微生物组存妊娠史生态记忆促后续妊娠速适。亚群落稳性非均：SC2与SC4早孕不稳由功能相似次生贡献者重组驱，显功能冗余与分类弥散却功能连贯之真生态单元，待宏基因组或代谢组证功能稳。产次调动态两签名：SC4仅初产妇时态变显差异响应，SC6初产妇全时高显持久结构改，生态记忆至少双机制——抗扰缓冲与结构持改，与前德国队列观一脉相承并指特定亚群落。其机制或宿主（免疫记忆T细胞预适应）、微生物（群落生态构型重适）或共适应，需多妊娠纵贯免谱与孕前微生组辨。产次时异重塑或影响婴定殖：经产妇婴6月龄Bifidobacterium主导SC1高，与前经产妇母Bifidobacterium高致垂直传增一致；高pBMI≥25降此SC1尤初产妇婴显，经产妇婴避部高pBMI效，与近期纵贯队列阴道娩婴6月龄内 parity 差 beta 多样性 一致。局限含厕纸拭样引阴道菌已滤除、SC3单隔证效、n=52力限缺早10周前及循环代谢物。结论为妊娠再稳亚群落间平衡，产次经生态记忆调制此程或宿主黏膜免致，母动态借婴肠道定殖具代际意，未来需整合黏膜免谱、代调与婴纵贯采以明母微生轨迹联母代谢健康与婴定殖结局。
结论部分原文翻译为：妊娠重塑稳定微生物亚群落间的平衡，产次通过看似生态记忆的机制调节此过程，可能由宿主黏膜免疫启动介导。这些母体动态可能通过其对婴儿肠道定殖的影响产生代际后果。未来研究整合黏膜免疫分析、母体代谢组学及纵向婴儿采样，可阐明母体微生物组轨迹与母体代谢健康及婴儿肠道定殖结局关联的机理。
