基于高通量测序的营养不良儿童肠道菌群多样性研究.pdf

175中国现代药物应用2024年1月第18卷第1期Chin J Mod Drug Appl,Jan 2024,Vol.18,No.1基于高通量测序的营养不良儿童肠道菌群多样性研究吴慧丹黄峥阳丽君江美花李扬帆梅桂斌【摘要】目的使用高通量测序模式和方法,对比分析营养不良儿童和正常儿童中肠道菌群物种存在的多样性和实际组成特点存在的差异。方法选取 20 例营养不良患儿作为 A 组,选择 20 例身体正常的儿童作为 B 组,采集分析粪便样本,提取细菌中的核糖体 DNA(RDNA),通过两组聚合酶链式反应(PCR)模式和方法,增加 V4V5 区域实现高通量的测序设计,并进行生物信息学的分析研究；构建相应的稀释曲线,使用 MOTHUR 软件实现 APLHA 的多样化的计算分析,通过 R 语言实现热图(HEAT MAP)分析研究,实现典型关联分析(CCA),之后分别在样本的门和属两个列表中实现数据的整合和归纳,进行统计分析。结果两组年龄、性别、分娩方式等一般资料比较,差异无统计学意义(P0.05)。两组生物多样性比较差异无统计学意义(P0.05)。两组生物多样性比较差异无统计学意义(P0.05)。两组样本中共发现微生物 154 种。A 组肠道菌群的操作性分类单元(OTU)为 320,B 组大肠产甲烷菌群的 OTU 为 348,以及其共同 OTU 为 290,两组肠细菌群门的细菌中共检出 14 个。为更好地展现检测结果,本研究依托柱状图优势,将各个样品进行分类学层面的比较和研究。A 组硬壁菌门及未分级菌(乳酸菌)的 OTU 少于 B 组,拟杆菌(多形状杆菌、罗斯菌)的 OTU 多于 B 组,差异有统计学意义(P0.05),其他菌群丰度比较差异有统计学意义(P0.05).There was no statistically significant difference in biodiversity between the two groups(P0.05).A total of 154 microbial species were found in the two groups of samples.The operational taxonomic unit(OTU)of intestinal flora of group A was 320,the OTU of group B colorectal methanogens was 348,and their common OTU was 290.A total of 14 species of intestinal flora were detected in the two groups of children.In order to better display the test results,this study relied on the advantages of bar charts to compare and study each sample at the taxonomic level.OTU of Firmicutes and ungraded bacteria(lactic acid bacteria)in group A was less than that in group B,and OTU of Bacteroides 基金项目：福建省卫生健康科技计划项目(项目编号：2021QNA034)作者单位：350012福建省疾病预防控制中心食品营养研究室实验研究176中国现代药物应用2024年1月第18卷第1期Chin J Mod Drug Appl,Jan 2024,Vol.18,No.1营养不足是当各种营养素由于摄食过度或摄入量不够或消化吸收障碍等各种因素引起产生的某种营养素失调症状1。营养不足是一种全球性的普遍性问题,在欧盟和美洲等发达国家以超重为主,非洲及东南亚等发达国家以身材矮小及生长发育迟缓为主,而南亚和东非等区域主要为出现 5 岁内的低体重儿,其中有多达 33%的儿童发生营养素不足和肠道的通透性显著提高,出现腹泻及吸收障碍等一系列问题,导致儿童的生长发育减缓2。就目前来说,营养素干预及相应的饮食调整方法无法切实可行的治疗营养不良对身体正常生长、发育产生阻碍作用,提升疾病发生几率。近期研究成果表明,微生物群落作为人类第二基因组,在很多层面影响了人类健康水平3。以肠道菌群为例,其共生在肠道中,能够提取食物中的营养物质,增强食物消化和吸收能力,尤其对肠道功能尚不健全的儿童来说,肠道菌群的完整性和丰富性,对其生长发育有着最为直接的影响。结合以往经验,人体内的肠道菌群结构组成和数量分布总体上处于动态平衡的状态,如果倡导微生态的平衡出现一定程度的紊乱,出现菌群失调等问题,对人体肠道运动能力、营养吸收能力和能量换化能力产生消极影响4。本文就主要使用高通量测序方法和模式,对比研究营养不良患儿和正常儿童中肠道菌群物种的多样性及差异情况,分析儿童营养不良和肠道微生物群落的具体联系性。1资料与方法1.1一般资料本研究选择 2022 年 16 月在本院门诊就诊中出现营养不良的患儿共 20 例作为 A 组,其中男 8 例,女 12 例。纳入标准：由于喂养不当、饮食搭配不合理等原因出现营养不良的情况,如挑食、偏食及没有科学的喂食方法等；排除标准：近半年出现传染性腹泻及抗生素应用。选取20例健康儿童视为B组,其中男 9 例,女 11 例。1.2诊断标准按照国际通用标准判定儿童营养不足有关情况,并在此基础上开展有效估测和预算。根据患儿体重判定。根据同年龄段孩子的体重标准作为参考,如果较正常体重低 15%25%之间,属于轻度营养不良；较正常体重低 26%40%则为中度(Bacteroides thetaiotaomicron and Rothia)was more than that in group B.The difference was statistically significant(P0.05).The abundance of other bacterial communities showed statistical significance(P0.05).Conclusion The number of Firmicutes in the intestinal flora of malnourished children is low,and the content of Bacteroides thetaiotaomicron and Rothia is high.Therefore,in order to ensure the stability and balance of the flora,it is necessary to adjust the colony by taking advantage of the actual characteristics of different bacteria genera to ensure that children have good absorption and physical development.【Key words】High-throughput sequencing;Malnutrition;Childrens intestines;Diversity of flora营养不良；较正常体重低 40%以上为重度营养不良5。根据皮下脂肪厚度情况判定：腹壁皮下脂肪厚度为0.40.8 cm为轻度营养不良,腹壁皮下脂肪厚度0.4 cm 为中度营养不良,全身皮下脂肪几乎消失为重度营养不良6。1.3方法1.3.1样本采集研究过程中严格按照规定的采样流程,使用采集盒开展采集工作。为防范样本污染情况,采集环节应当科学运用采集勺,采集重量为 5 g 左右的粪便,并放入到具有保存液的保存管中,并马上转入到-80的环境保存7。1.3.2DNA 的提取和测序采集和使用 QLAAMP 试剂盒进行粪便样本的抽取提炼,使用浓度为 12%的琼脂凝胶进行抽取,得到基因组的 DNA,通过 PCR 的模式实现文档构建,利用 525 和 926 引物进行扩增,之后使用相关检测平台,将接头、测序引物和其他标签列加入到两端位置,使用美国 AXYGEN 公司生产的凝胶回收试剂盒进行回收,使用实时荧光定量 PCR 仪 FTC-3000 技术进行荧光的定量分析,在均一化混合完成后通过文库进行构建研究,使用 ILLUMIAN 进行文库的测序。1.3.3稀释性曲线对研究过程中获得的系列数据,借助 TRIMMOMATIC 软件对低质量序列进行去除,对去除后的序列使用 FLASH 软件开展序列的二次拼接,在此基础上,借助 MOTHUR 对序列开展质量控制,从而获得最优化的序列8。从实际效果来看,MOTHUR软件对数据的处理可以有效排除冗余数据、干扰数据的影响,形成最具代表性的 OTU 稀释性曲线,以更加真实地反映出质量系列变化情况,最大程度地保证测序数据的丰富程度,为后续研究工作高质量开展提供便利9。如果曲线最终达到平稳后表明测序科学合理,对优化的序列通过 97%的相似性进行对比分析,并进行相应的信息标注处理。1.3.4多样性分析本研究借助 MOTHUR 软件,对ALOHA 开展多样性分析和计算,从最终研究获得的ALPHA 指标中,评估分析 CHAO 和 ACE 的指数变化情况10。结合过往经验,CHAO 和 ACE 指数越大,表177中国现代药物应用2024年1月第18卷第1期Chin J Mod Drug Appl,Jan 2024,Vol.18,No.1明 SIMPSON 实际数值越小,检测样本中的微生物种类越多、类型越丰富。使用 R 语言进行分析,箱线图实现展示样品组中存在的多样化分布,实现数据的分布情况,具体差异使用 P 值进行表示,箱体中也就证明该组数据中50%的分布区间和中间线都确定为中位数,如果范围 1.5 箱体则为最大值或最小值11。1.4统计学方法采用 SPSS20.0 统计学软件对研究数据进行统计分析。计量资料以均数 标准差(x-s)表示,采用 t 检验；计数资料以率(%)表示,采用2、q、w 检验或 Fisher 确切概率法。P0.05),具有可比性。见表 1。2.2两组物种丰度及生物多样性分析两组粪便中共得到1394074个有效序列,优化序列数量为1268032个,平均长度 410。通过对上述序列分析,实现序列数当代表的 OTU 进行稀释性的结构曲线分析,可见大部分样本呈平坦趋势达到饱和,但部分样本曲线没有达到平坦,证明两组微生物种类都较为丰度。通过箱线图像进行显示,两组生物多样性比较差异无统计学意义(P0.05)。证明两组肠道菌群的丰度和均匀度无显著差异。见图 1。表 1两组一般资料比较n,x-s,n(%)组别例数性别平均月龄(个月)剖宫产早产儿母乳喂养男女A 组2081218154(20)3(15)4(20)B 组2091122122(10)2(10)2(10)2/t0.1020.9310.7840.2290.784P0.7490.3580.3760.6330.376注：两组比较,P0.05图 1OTU 水平系数曲线和 Shannon 指数图序列数OTU 数值序列数Shannon 多样性Shannon 指数图曲线多样本稀释曲线2.3两组肠道菌群物种分析两组样本中共发现微生物 154 种。A 组肠道菌群的 OTU 为 320,B 组大肠产甲烷菌群的 OTU 为 348,以及其共同 OTU 为 290,两组肠细菌群门的细菌中共检出 14 个。为更好地展现检测结果,本研究依托柱状图优势,将各个样品进行分类学层面的比较和研究。A 组硬壁菌门(Firmicutes)及未分级菌(乳酸菌)的OTU少于B组,拟杆菌(Bacteroidetes)(多形状杆菌、罗斯菌)的 OTU 多于 B 组,差异有统计学意义(P0.05),其他菌群丰度比较差异有统计学意义(P0.05)。见表 3。178中国现代药物应用2024年1月第18卷第1期Chin J Mod Drug Appl,Jan 2024,Vol.18,No.1表 2两组属水平肠道菌群 OTU 比较组别例数多形状杆菌乳酸菌罗斯菌x-方差标准误x-方差标准误x-方差标准误A 组200.4734070.0197520.0314260.00007400.0000030.0126630.0004010.004476B 组200.2555370.0441700.0469950.0142620.0035750.0133700.0025610.0000410.001424q0.0170.0610.095P0.0030.0160.033注：表 3两组属水平肠道菌群比较(n,%)组别例数多形杆状菌未分类菌隐秘杆菌属克雷伯杆菌韦荣球菌属OTU分布(%)OTU分布(%)OTU分布(%)OTU分布(%)OTU分布(%)A 组2027315448.786534814.27332585.57241162.85219862.49B 组2048653637.7425054219.43663035.14439343.40396983.08w363121228168175P0.0010.0340.4570.3940.508注：2.5聚类分析通过研究不同样本 OTU(97%相似性)组成反映样本之间的差别与距离,主成分分析(PCA)通过方差分析,将各组样本差别表示在二维坐标图上,坐标轴为可以最大地表示方数差的两个特征值,如样品成分越是接近,反映到 PCA 曲线上的时间越 近12。不同条件中的样品可以显示出离散与集中的分布关系13。对两个样品的主要组分进行 PCA 研究,根据样品聚集现象,两组尽管高度接近,但仍可发现两组聚集。根据国民营养计划(2017-2030 年)等相关政策文件,为最大程度地保障人体生命健康,在生命开始最早的 100 d 就应当开展保健活动,通过营养充分供给,稳步增强孕妇、婴幼儿的营养供给水平。“生命早期1000 d”是指280 d胎期,在婴儿出生后2周岁(720 d)。这1000 d是新生儿神经系统等人体组织的关键发育期,也是肠道菌群形成的关键时期。这一时期的肠道菌群极易受到营养状况、病情变化等因素影响14,15。为更好地应对上述问题,需要加强身体健康状况监测,一旦发现问题要及时进行应急处理,使得肠道菌群活性、丰度达到预期,更好地服务于婴幼儿生长发育16。从本次研究过程中获得数据来看,借助高通量的16S rDNA 测序技术,对肠道菌群进行系统分析,开展多样性研究,能够准确发现营养不足儿童和营养正常儿童体内肠道菌群状况,为后续研究活动有序开展提供数据支持17。本研究中,在门水平上,A 组含有较多的拟球菌,B 组含有较多硬壁菌门；在属水平上,两组 24 个菌株比较差异有统计学意义(P0.05),考虑可能与样本数量低有关,但仍有深入探讨的价值。高通量 16S rDNA 的主要优点在于为研究生物组织多样性和研究生命系统发育分类信息学资源,主成分 PCA 研究结果显示 A 组与B 组尽管距离较接近,但仍可发现两组聚集,表明两组间在物质结构方面仍存在着相当的差异,具有进一步深入研究探索的重要意义。但对 16S rDNA 的基因组研究,对推断生物的结构仍具有一定局限性。因为对肠道菌群的功能串扰及对数据的冗余处理,仅通过测定产甲烷菌的胆固醇含量就很难合理、精确地推测其生理机制,因此目前尚无法确定与营养不足患儿的相关致病微生物。本研究方法不足就是对缺乏儿童疾病机理的研究,为进一步认识营养不足儿童肠道产甲烷菌与疾病之间的相互关联和发生机理,本研究进一步对儿童肠道细菌中的基因组、蛋白质组和代谢组成等方面进行相应的研究和探索。系统化分析儿童肠道内甲烷菌与食物消化、人体代谢、免疫作用发挥等方面的作用,更加明确关键性菌群在营养不良患儿治疗中的作用和意义。结合过往经验,对营养不良患儿整个治疗过程中可以采取菌群治疗方案,通过采取有效举措,平衡菌群关系,在体内形成良好的微生物环境。本次研究选择 2022 年 16 月在本院儿童营养门诊的 20 例营养不良患儿作为 A 组,选择 20 例身体正常的儿童作为 B 组,采集分析粪便样本,提取细菌中的 RDNA,通过两组 PCR 的模式和方法,增加 V4V5区域实现高通量的测序设计,且进行生物信息学的分析研究,构建相应的稀释曲线,使用 MOTHUR 软件实现 APLHA 的多样化的计算分析,通过 R 语言实现HEATMAP 分析研究,实现典型关联分析 CCA,之后分别在样本的门和属两个列表中实现数据的整合和归纳,进行统计分析。两组年龄、性别、分娩方式等一般资料比较,差异无统计学意义(P0.05)。两组肠道菌群的丰度和均匀度比较差异无统计学意义(P0.05)。两组肠道菌群共检测到 14 种,但 A 组硬壁菌门和未分级菌的 OTU 明显低于 B 组,且拟杆菌中的 OTU 数量超过 B 组,差异有统计学意义(P0.05)。两组肠道菌群共检测出 154 种微生物,具体分析后发现,乳酸菌、罗斯菌等 24 个菌群的 OTU 差异明显,其中 A 组杆状菌和罗斯菌的 OTU 超过 B 组,乳酸菌 OTU 低于 B 组。两组样本成分进行 PCA 分析,根据具体样本聚合情况,两组虽然较为接近,但也具有差异。出现营养不良的儿童肠道菌群中硬壁菌的数量较低,其中多形杆状菌和罗斯菌的含量较高,为保证菌群的稳定和平衡,需要利用不同菌属的实际特点进行菌落的调整,保证儿童具有良好的吸收和体格发育。肠道菌群深刻地关系着宿主的营养加工、能量平衡、免疫活动、胃肠道发育和其他许多关键的生命行为。人体肠道菌群的形成在婴幼儿阶段就开始进行训练,和健康青少年一样,婴幼儿的消化道内微生物菌群形成方式截然不同,婴幼儿消化道菌群定植和成长都是系统的动态变化,因此消化道菌群对婴幼儿生长发育和免疫形成很重要,婴幼儿时期消化道菌群失衡除患正常的代谢类疾病外,也可罹患一些细菌类疾病,如流感、肺炎等。研究证实,在婴幼儿阶段,营养素的正常摄入量对保证肠道产甲烷菌的生长平衡有重要影响,而矿物质(如钙、铁、镁、锌)的正常营养利用率则对宿主肠道产甲烷菌的形成起重要作用,肠道产甲烷菌与维生素代谢相互作用,提高了机体的身心健康度。而营养包是促进中国贫困地区 624 月龄婴幼儿健康成长重要的营养补充品,含有蛋白质、脂肪、矿物质、微量元素等营养素,因此本文将运用高通量测序方法研究营养包对婴幼儿阶段肠道产甲烷菌的作用。肠道菌群多样性指肠道中的一些微生物,如肠道中有乳酸杆菌、双歧杆菌等,这些微生物都能够促进肠道的蠕动,能够防止出现便秘的情况。这些益生菌还能增强抵抗能力,出现肠道菌群失调后,患者可能会出现腹泻及便秘情况,对患者身体造成的影响是较大的。因此,平时一定要做好相关防护措施与调整,防止出现肠道疾病。患上肠道疾病后应当及时采取对症治疗,才有助于病情控制。对人体肠道中具有数以亿计大量的共生菌群,其自身的种类数量较多,可以达到上千种,且自身数量甚至达到人体细胞总量的10倍以上,至今为止,科学发展至今还有超过 80%的微生物没有180中国现代药物应用2024年1月第18卷第1期Chin J Mod Drug Appl,Jan 2024,Vol.18,No.1检测到,对此类肠道微生物其自身和人体存在互利共生的情况和关系,可以维持和保证人体健康的发展具有积极作用,肠道菌群在肠道中一直保证这动态化的平衡,可以实现微生物的合成,帮助人体更好地吸收营养情况,维持肠道免疫功能得到更好的维持,可以抵抗有害微生物的具体侵害问题,对此类平衡也会因为一部分因素导致动态平衡紊乱,进而可能出现肥胖与和肠炎等一系列疾病,为可以更好地强化和加深人类疾病的发生原因和药理学认知,自从人类基因组计划后,科学化开始了大量的原基因组计划体系,让肠道菌群研究取得更好的发展。在传统肠道微生物群落多样性研究的过程中,科学家主要使用 DGGE 之类的分子生物法实现对肠道菌群的研究分析,但对此类的方法其自身都存在一定缺陷,如 DGGE 只能简单地检测环境样品中存在的 10 多种优势菌群,但对弱势的微生物检测不够精确,且电泳条带中存在着多种 16DNA 的排列序列。为更好地了解到具体菌种的详细信息还需要进行克隆检测、实验室实验等多种方法,且随着科学技术的不断提高,对微生物检测的精准性要求不断提高,无法实现对微生物丰度的检测。但生物芯片其自身可以实现通过固定的探针得到微生物的多样性信息,通过信号的情况了解到微生物的丰度,但准确性依然不足,高通量测序技术最早发布于 2005 年,具有稳定性强、检测率高等一系列优势,当下被大量使用在肠道菌群的检测中,学术界中大量学者开展了详细的研究,其中大部分研究结果都是在国际顶尖杂志刊登,是当下进行微生物多样性研究最好的科学方法。参考文献1 曹丹丹,马灌楠,董梦醒,等.肠道菌群与儿童营养不良关系的研究进展.发育医学电子杂志,2017,5(1):45-49,59.2 杨光军,惠宇,陈欣亦,等.海口市学龄前儿童营养不良现状及影响因素分析.中国公共卫生管理,2022,38(3):369-371.3 胡黎舒,魏兵,廖世峨,等.先天性心脏病患儿术后合并营养不良影响因素及其生存质量分析.中国儿童保健杂志,2021,29(8):824-828.4 蒲吉慧.营养不良小儿的饮食应怎样要求.人人健康,2022(14):55.5 杨萍.及时干预儿童营养不良.家庭医药(快乐养生),2020(4):14.6 晏建东.谨防“营养不良”症.江西教育,1996(12):7.7 郑惠槟,陈峥嵘,何翠娟,等.儿童营养不良的病因、病原体和诱发因素年代变迁的研究.现代临床医学生物工程学杂志,2001,7(4):263-266.8 林华锦.75 例 3 岁以下营养不良儿血铅、锌、铜、铁、钙含量的变化.广东微量元素科学,1996(4):25-27.9 叶美琴.205 例 06 岁营养不良儿童调查分析.浙江预防医学,1996(4):29-30.10 陈绮霞.婴幼儿重症营养不良 138 例治疗体会.广东医学,1995(12):847-848.11 吴丹,任彩霞,徐文,等.基于肠道菌群探讨益生菌对儿童营养不良的影响.药物生物技术,2023,30(2):216-220.12 丁姮月,朱惠萍,梁国强,等.基于高通量测序技术的脾虚腹泻型肠易激综合征小鼠与健康小鼠肠道菌群的差异研究.中国中医基础医学杂志,2021,27(2):260-324.13 杨国平,叶敏霞,张荣,等.基于高通量测序的鼠类肠道细菌群落多样性分析.中华卫生杀虫药械,2021,27(2):148-153.14 宫幼喆,钟雪梅,朱丹,等.基于 16S rDNA 高通量测序技术对克罗恩病患儿肠道菌群变化的研究.中国微生态杂志,2021,33(9):993-999.15 曾佩佩,曾婷,邓梁琼,等.基于高通量测序分析孤独症谱系障碍儿童肠道菌群的构成变化.检验医学,2021,36(2):153-161.16 黄焕迪,秦华珍,何瑞坤,等.基于高通量测序对中药治疗胃肠道疾病中肠道菌群研究进展.辽宁中医药大学学报,2021,23(8):65-69.17 刘智程,吴德,屈爱娜.孤独症谱系障碍儿童肠道菌群多样性研究及功能预测分析.中国当代儿科杂志,2022(12):1356-1364.18 凌宝殿,谢志军,张文娟,等.基于高通量测序技术研究原发性胃癌患者胃肠道菌群特征.赣南医学学院学报.2023(3):228-234.19 陈文辉,于龙,刘文韬,等.基于高通量测序技术的胆囊癌患者肠道菌群变化特点.河南医学研究,2022(24):4431-4437.20 钱香,王莉,齐艳,等.高通量测序技术对腹泻患儿使用抗生素前后肠道菌群变化的研究.中文科技期刊数据库(全文版)医药卫生,2021(7):229-231.收稿日期：2023-08-03