支气管肺发育不良(bronchopulmonary dysplasia, BPD)是早产儿常见的慢性呼吸系统疾病,发病机制复杂,且缺乏有效的治疗手段,仍以预防为主,寻找一种新型的、非侵入性的生物标志物用于BPD的早期预测仍然是当务之急。肠道菌群被认为是人体的第二个基因组,研究表明,正常的肠道和肺部共生菌有助于维持免疫动态平衡,即存在微生物-肠-肺轴,肠道菌群失调可能影响肺的炎症及免疫反应
[1]。与足月儿相比,早产儿肠道菌群多样性低、定植速度慢、成熟晚,易受外界因素影响,导致菌群的定植与发育模式发生改变
[2]。越来越多的研究表明,哮喘(肠道微生物多样性降低,双歧杆菌、嗜黏蛋白阿克曼菌和普拉梭菌的丰度减少)、呼吸道感染(如结核病患儿肠道放线菌门和变形菌门增加,拟杆菌门减少)等慢性肺部疾病的发生发展与早期肠道内菌群的失调相关联
[3]。BPD亦属于慢性肺部疾病,由于早产儿肠道菌群定植所受影响因素多、标本难以留取等,目前国内外关于早产儿BPD与肠道菌群的研究较少,肠道菌群是否会影响对BPD的易感性目前尚不清楚。本研究应用高通量测序技术分析不同严重程度BPD极早产儿(very preterm infant, VPI)肠道菌群的定植模式以及菌群结构差异,为临床预测、预防及治疗BPD提供了新思路。
1 资料与方法
1.1 研究对象
本研究采用前瞻性单中心队列研究,选取2023年1月—2024年6月厦门大学附属妇女儿童医院新生儿重症监护室收治的胎龄<32周且出生体重<1 500 g的早产儿作为研究对象。随访VPI的临床结局,根据2018年美国NICHD修正后的诊断标准
[4],BPD组分为3个亚组:Ⅰ级、Ⅱ级和Ⅲ级。比较无/Ⅰ级BPD组和Ⅱ/Ⅲ级BPD组肠道菌群组成的差异。
样本量计算:首先用简单随机抽样样本量计算公式计算研究所需样本量,再考虑整群抽样设计效应deff。简单随机抽样的样本量计算公式如下:
上式中,检验水准
范围为0.001~0.05,允许误差
范围为0.001~0.01,本次调查中的BPD发生率
,Scientific Reports的研究结果显示,VPI BPD的发生率
为22.1%
[1]。因此,本研究中,取
=0.05、
=0.1来计算研究最小样本,得出总样本量为66例。考虑到整群抽样设计效应
deff为1.3,样本含量至少为86例。本研究最终收集样本93例,符合最小样本量要求。
入选标准:(1)住院时间≥纠正胎龄36周;(2)生后24 h内入院;(3)单胎;(4)纯母乳喂养或母乳+捐赠母乳。排除标准:(1)存在先天畸形或遗传代谢病;(2)出生时有重症感染、重度窒息、胃肠道功能障碍;(3)围手术期或住院期间行手术者;(4)住院期间死亡、中断治疗或自动出院者;(5)因肠道问题在1周龄内无法开始肠内喂养者;(6)出生后口服过肠道益生菌类药物;(7)母亲分娩前1个月内使用微生态制剂;(8)资料不完整。本研究已通过厦门大学附属妇女儿童医院伦理委员会批准(KY⁃2025⁃100⁃K02),并获得患儿家属的知情同意。
1.2 资料收集
收集研究对象临床资料:(1)围生期及母孕期并发症指标:胎龄、出生体重、产前使用糖皮质激素、延迟脐带结扎、性别、分娩方式、小于胎龄儿(small for gestational age infant, SGA)
[5]、母亲产前抗生素使用等;(2)生长及营养状况指标:开始肠内营养时间、达全肠内营养日龄
[6]、恢复出生体重日龄、平均体重增长速率
[7];(3)住院期间主要治疗指标:抗生素累计使用时间、有创/无创机械通气时间、总用氧时间;(4)住院期间合并症及并发症:新生儿呼吸窘迫综合征(neonatal respiratory distress syndrome, NRDS)
[8]、有血流动力学改变的动脉导管未闭
[9]、≥Ⅱ期新生儿坏死性小肠结肠炎(neonatal necrotizing enterocolitis, NEC)
[8]、早发型败血症(early⁃onset sepsis, EOS)
[8]、Ⅲ~Ⅳ级脑室内出血(intraventricular hemorrhage, IVH)
[8]、宫外生长迟缓(extrauterine growth restriction, EUGR)
[10]。
1.3 粪便标本采集
采用一次性无菌粪便冻存管收集住院VPI在不同时间点[生后第7天(D7)、第21天(D21)、纠正胎龄36周]的粪便样本,编号后在1.5 h内立即放入-80℃冰箱保存待用。
1.4 DNA抽提和PCR扩增
采用 MagPure Soil DNA LQ Kit(Magan)试剂盒提取粪便样本的基因组DNA。利用NanoDrop 2000(Thermo Fisher Scientific, USA)和琼脂糖凝胶电泳检测DNA的浓度和纯度。以提取的基因组DNA为模板,使用带Barcode的特异引物和Takara Ex Taq高保真酶进行细菌16S rRNA基因的聚合酶链反应(polymerase chain reaction, PCR)扩增。采用通用引物343F(5'⁃TACGGRAGGCAGCAG⁃3')和798R(5'⁃AGGGTATCTAATCCT⁃3')扩增16S rRNA基因的V3⁃V4可变区,用于细菌多样性分析。
1.5 文库构建和测序步骤
首先对PCR扩增产物进行琼脂糖凝胶电泳检测,随后使用AMPure XP磁珠进行纯化。将纯化后的产物作为模板,进行第二轮PCR扩增,并对扩增产物再次进行磁珠纯化。取纯化后的二轮产物进行Qubit定量,根据定量结果调整浓度后上机测序。使用Illumina NovaSeq 6000测序平台进行测序,并生成250 bp双端reads。测序由上海欧易生物技术有限公司完成。
1.6 生物信息学分析
数据下机后,首先使用Cutadapt软件,将raw data序列剪切掉引物序列。然后使用DADA2,将上一步合格的双端raw data按照QIIME 2 (2020.11)默认参数进行质量过滤、降噪、拼接及去嵌合体等质控分析,得到代表序列及ASV丰度表格。使用QIIME 2软件包将所有代表序列与Silva(version138)数据库进行比对注释。物种比对注释使用q2⁃feature⁃classifier软件默认参数进行分析。采用QIIME 2软件进行α和β多样性分析。利用R计算的weighted Unifrac距离矩阵,进行weighted Unifrac主坐标分析(PCoA)来评估样本的β多样性。
1.7 统计学分析
基础数据分析采用SPSS 26.0软件进行统计学分析。计数资料采用例数和百分率(%)表示,组间比较采用检验或Fisher确切概率法。正态分布的计量资料采用均数±标准差()表示,组间比较采用两样本t检验;非正态分布的计量资料采用中位数(四分位数间距)[M(P25,P75 )]表示,组间比较采用Mann⁃Whitney U检验。肠道菌群生物信息学分析采用R(v3.5.1)包,P<0.05为差异有统计学意义。
2 结果
2.1 两组基础资料比较
本研究共纳入118例VPI,最终入组93例,其中无BPD组20例,BPD组73例(Ⅰ级32例,Ⅱ级26例和Ⅲ级15例)。男婴41例(44.1%),女婴52例(55.9%),胎龄29.9(28.2,31.1)周,出生体重1 265(969,1 437)g。Ⅱ/Ⅲ级BPD组的胎龄和出生体重显著降低,NRDS的发生率更高,需要更长时间的有创和无创机械通气,更长的用氧时间(
P<0.05)。此外,Ⅱ/Ⅲ级BPD组达全肠内营养时间延长,EUGR的发生率更高(
P<0.05),见
表1。
2.2 多因素logistic回归分析
以分组为因变量(0=无/Ⅰ级BPD,1=Ⅱ/Ⅲ级BPD),单因素分析结果具有统计学意义的变量作为多因素分析的自变量,拟合二分类多因素logistic回归模型,分析Ⅱ/Ⅲ级BPD的影响因素,结果显示,出生体重、NRDS共2个变量进入了模型。具体为:出生体重越高,越不容易出现Ⅱ/Ⅲ级BPD(
OR=0.993,95%
CI:0.991~0.996);NRDS是发生Ⅱ/Ⅲ级BPD的危险因素(
OR=3.535,95%
CI:1.099~11.369)。见
表2。
2.3 粪便样本收集及ASV分析
共收集粪便样本279份,其中无/Ⅰ级BPD组在3个时间点共收集156份,Ⅱ/Ⅲ级BPD组123份。两组VPI在ASV水平上的Venn图结果显示,无/Ⅰ级BPD组3个时间点分别有12 540、2 215、1 786个物种,Ⅱ/Ⅲ级BPD组则有8 849、3 101、1 426个物种,6个亚组所共有的物种数量为229个。
2.4 物种组成及差异分析
门水平上,两组VPI均以变形菌门、厚壁菌门、拟杆菌门、放线菌门为主,占全部菌门的90%以上。随时间推移的自身对照中,Ⅱ/Ⅲ级BPD组变形菌门丰度呈上升趋势,厚壁菌门和拟杆菌门丰度呈下降趋势,无/Ⅰ级BPD组趋势不明显,见
图1A、图
1C和
1D。
属水平上,两组VPI患儿各时点肠道内菌群主要由克雷伯菌属、大肠埃希菌属-志贺菌属、肠球菌属、葡萄球菌属、狭义梭菌属1、未分类
Muribaculaceae、链球菌属、双歧杆菌属组成。随着时间的推移,Ⅱ/Ⅲ级BPD组葡萄球菌属相对丰度逐渐下降,克雷伯菌属相对丰度呈上升趋势,无/Ⅰ级BPD组上述趋势不明显(图
1B、
1E~F)。
以属水平的角度进行差异菌群分析,在丰度前15位的菌属中,两组VPI患儿肠道菌群结构在出生后D7和纠正胎龄36周差异无统计学意义。在出生后D21 Ⅱ/Ⅲ级BPD组链球菌属较无/Ⅰ级BPD组显著减少(
P=0.017),见
图2。
2.5 α多样性分析
与无/Ⅰ级BPD组比较,Ⅱ/Ⅲ级BPD组物种丰富度指数(ACE指数和Chao1指数)和多样性指数(Shannon指数和Simpson指数)在生后D7和D21有下降趋势,但差异无统计学意义(
P>0.05)。然而,在纠正胎龄36周时,Ⅱ/Ⅲ级BPD组物种丰富度指数和多样性指数显著低于无/Ⅰ级BPD组,差异有统计学意义(ACE指数,
P=0.012;Chao1指数,
P=0.016;Shannon指数,
P=0.014;Simpson指数,
P=0.048),见
图3。
2.6 β多样性分析
主坐标分析(principal coordinates analysis, PCoA)显示,两组VPI肠道菌群的组成在生后D7和D21相似(
P=0.762、
P=0.442)。在纠正胎龄36周时,两组早产儿肠道菌群的组成出现明显的差异(
P=0.019),在PCoA图中可区分大部分样本,其中PC1可解释35.7%的差异。见
图4。
3 讨论
近年来,肠道和肺部的微生物群及其代谢物已成为BPD发病机制中潜在的关键因素。肠道菌群失调和随之而来的营养不良可能通过改变表观遗传学、促进炎症和氧化应激、改变肠道通透性、改变血管和淋巴管发育,以及影响微量营养素等机制参与BPD的发生和发展
[11-12]。早产儿肠道菌群的发展过程受到多种围产期因素的影响,其中最重要的是胎龄、分娩方式、喂养方式和抗生素。本研究中,Ⅱ/Ⅲ级BPD组的胎龄和出生体重更低,NRDS和EUGR的发生率更高,无创呼吸机和用氧时间更长,其余围产期因素比较差异均无统计学意义。这意味着虽然两组VPI优势肠道菌群相似,但其丰度存在差异。
早产儿正常厌氧菌定植延迟、潜在病原菌定植增加以及肠道短链脂肪酸组成的改变,可能是由于其免疫系统发育不成熟、肠道氧化环境改变以及新生儿重症监护室环境中的多种医疗干预所致
[13]。本研究结果显示,VPI双歧杆菌属的相对丰度始终较低,均<5%,随着时间推移,拟杆菌门的相对丰度下降,变形菌门、克雷伯氏菌属、大肠埃希菌属-志贺菌属、肠球菌属成为优势菌群较长时间定植在VPI肠道,这与许景林等
[14]和Khan等
[15]的研究结果相似。Chen等
[16]病例对照研究显示,生后第7天至第28天,BPD患儿肠道菌群多样性和相对丰度下降,变形菌门丰度明显增加,拟杆菌门丰度明显下降。Gallacher等
[17]研究显示,接受机械通气的VPI(其中84%发展为中-重度BPD)随着时间推移,肠道微生物区中变形菌门丰富,而厚壁菌门数量下降。本研究中Ⅱ/Ⅲ级BPD组患儿机械通气时间和用氧时间更长,可能会干扰出生后第一个月内肠道内的厌氧环境,进一步造成患儿生长障碍
[18]。其变形菌门(相对丰度分别为30.3%、67.8%、75.7%)和克雷伯氏菌属(相对丰度分别为10.9%、33.1%、40.1%)呈上升趋势,厚壁菌门(相对丰度分别为40.6%、24.6%、20.8%)、拟杆菌门(相对丰度分别为21.2%、5.1%、1.4%)和葡萄球菌属(相对丰度分别为16.7%、2.2%、0.1%)呈下降趋势,无/Ⅰ级BPD组趋势不明显,与上述文献
[16-17]相符。变形菌门属于兼性厌氧菌,厚壁菌门和拟杆菌门属于专性厌氧菌。一般认为变形菌门在新生儿肠道发挥定植作用后,相对丰度会降低,这是微生物初始定植的正常组成部分;而这种定植模式的紊乱与新生儿疾病风险增加有关,比如NEC、新生儿晚发型败血症
[19-20]。变形菌门高丰度代表了菌群失调的微生物特征,反映了肠道微生物群落结构的不稳定性
[21]。克雷伯氏菌属于变形菌门,作为条件致病菌,通过产生脂多糖被肠道中的Toll样受体4受体识别,从而导致促炎性Th17细胞的定向募集和促炎性细胞因子(白细胞介素⁃17、白细胞介素⁃22)的释放,诱导肺泡内皮细胞免疫激活、肺细胞特异性分子信号失调和肺发育畸形,导致BPD
[22-23]。葡萄球菌属于厚壁菌门,葡萄球菌促进了克雷伯氏菌在肠道中的定居,而克雷伯氏菌抑制了葡萄球菌的生长
[23]。此外,拟杆菌门和厚壁菌门是乙酸、丙酸、丁酸等短链脂肪酸(short chain fatty acids, SCFAs)的主要生产者,SCFAs通过免疫调节、屏障保护、表观遗传调控及代谢支持等多途径缓解BPD的病理进程
[12]。Stricker等
[24]发现,肠道中产SCFAs的细菌含量较高的婴儿患BPD的风险较低。这些结果提示肠道菌群的异常定植增加了BPD的严重程度,其机制可能与SCFAs生成减少及促炎信号通路激活有关。
本研究结果发现,两组差异菌群为链球菌属。链球菌属属于厚壁菌门,是兼性厌氧菌。某些链球菌种类,如肺炎链球菌是肺部疾病的重要致病菌,而其他种类则可能在肺部微生物群落中发挥保护性作用。链球菌主要产生乳酸和短链脂肪酸(乙酸、丙酸和丁酸),这些物质可抑制大肠杆菌等病原体在上皮细胞中的生长,还能与宿主的免疫系统协同作用,减少肠道炎症
[25-26]。研究表明,肠道菌群中的链球菌属有助于维持正常免疫反应,如嗜热链球菌可能减轻如哮喘等过敏性肺部疾病的症状
[27]。一项系统性综述指出,早产儿肠道链球菌与早期生长呈正相关
[28]。本研究发现,NRDS是发生Ⅱ/Ⅲ级BPD的危险因素。根据Fu等
[29]的研究,NRDS患儿表现出肠道菌群多样性的增加,通过增加链球菌、双歧杆菌和梭菌属的丰度,显著修复了NRDS婴儿的肠道菌群结构。以上研究均显示了链球菌属对机体有利的一面,故推测链球菌属的减少可能加重BPD的严重程度。
本研究结果显示,Ⅱ/Ⅲ级BPD组肠道微生物群落α多样性在纠正胎龄36周时显著低于无/Ⅰ级BPD组,与Chen等
[16]和Zhang等
[30]的研究结果相似。肠道菌群多样性高有利于建立稳定群落,防止单一物种过度增殖。VPI肠道微生物多样性普遍降低,BPD患儿更甚。多样性降低易出现菌群失调,肠道通透性增加导致条件致病菌的过量定植,同时菌群失调产生多糖A、鞘脂和其他被免疫细胞识别的代谢物,导致局部或全身炎症反应和代谢紊乱,影响肺泡及微血管发育并导致BPD进展
[11,31]。PCoA分析显示,两组菌群组成差异在纠正胎龄36周时明显区分开来。结果表明不同程度BPD早产儿肠道菌群可能存在着不同的定植和发育模式。
本研究局限之处:(1)本研究为单中心前瞻性队列研究,样本量小,且未分析菌群代谢组或宿主免疫指标;(2)肠道菌群可能受到一些混杂因素的影响,受样本量的限制,我们未进一步按胎龄、出生体重、喂养方式和抗生素使用进行分组分析,我们根据现有数据谨慎得出结论,将在今后的研究中扩大样本量进行验证;(3)肠道微生物组的组成可能因种族或地理区域而异,本研究中所有婴儿均来自同一家医院的新生儿重症监护室,这可能是一些常见微生物的来源,未来需要更大规模的多中心研究结合多组学技术(如宏基因组、代谢组)明确关键菌群的功能作用。
综上所述,本研究揭示了不同严重程度VPI的肠道菌群特征,研究结果为BPD的微生态干预提供了理论依据,例如通过靶向调节变形菌门或补充益生菌改善肠道免疫稳态。未来需开展多中心大样本研究,结合代谢组学进一步验证微生物在BPD发展中的作用机制。
厦门市自然科学基金项目(3502Z20227408)
厦门市自然科学基金项目(3502Z20227297)
福建省2023年省级临床重点专科建设项目(新生儿科)(闽卫医政函[2023]1163号)