尚玮璇,刘 璐,雷素珍,郑宝东,张 怡,曾红亮
(福建农林大学食品科学学院,福建 福州 350002)
人体肠道是一个大型、复杂的生态系统,肠道菌群的数量、比例、功能等常受到外在因素的影响,如饮食、生活方式等,菌群失调可能引起机体天然免疫系统受损、肠内稳态失衡,有诱发肠道炎症、动脉粥样硬化或胰岛素抵抗的风险,导致2型糖尿病、高血糖等病症的发生。小肠是人体消化吸收的重要场所,主要通过酶类物质作用糖类、胆固醇、蛋白质等营养物质,小肠微生物数量远低于大肠,然而有研究表明,小肠菌群的失调在各类疾病的发生发展中有着独特的作用机制:1)小肠菌群失调可能导致小肠绒毛排列不规则并发生不同程度的断裂,使黏膜通透性增加,屏障功能下降;2)小肠上皮细胞中大量表达钠-钾-ATP酶,该酶在菌群失调时含量下降,导致细胞渗透压改变,细胞结构受到破坏,这两种机制可能协同破坏小肠对糖类、蛋白质、脂肪等物质的吸收作用,引发糖尿病等疾病的发生。因此通过改善饮食等途径增加小肠内益生菌数量、控制有害菌繁殖可作为维持小肠稳态、调节体内代谢、预防疾病的重要举措。
抗性淀粉(resistant starch,RS)可作为底物被结肠微生物发酵并利用,因此具有促进益生菌生长及增强其活性的益生元作用,且能够与其他益生元如膳食纤维相互作用发挥其益生元效应。近年来,已有许多研究报道了RS的功能特性,如张亚楠等研究发现,RS不仅可以增加大鼠体内的、等益生菌相对丰度,还有利于预防肥胖。莲子补脾益肾,对控制2型糖尿病人血糖和胰岛素水平有显著影响,且莲子中直链淀粉含量高达42%,是提供RS的良好来源。经过压热、微波等处理后的莲子抗性淀粉(lotus seed resistant starch,LRS)表面具有特殊的沟壑状结构,能提高益生菌在不良肠道环境的适应性,使其更好地利用碳源,促进自身增殖,并能通过参与胰岛素分泌、调节抗氧化活性等各种关键因素的表达水平发挥预防炎症等益生作用。Qin Renbing等研究发现,III型RS有利于肠道菌群发酵产生乳酸。乳酸盐是RS经肠道菌群发酵后产生的一种后生元,其益生机制主要是调节宿主菌群、改善肠上皮屏障功能、调节全身代谢反应等。目前对LRS和乳酸钠(sodium lactate,SL)的生理功能已进行了一定研究,但对于LRS和SL的协同作用尚鲜有报道,为研究两者是否存在协同益生作用,本实验以LRS和SL为饮食干预,研究大鼠小肠菌群、代谢谱的影响,旨在探讨LRS及SL影响小肠菌群及代谢的潜在机制,从饮食干预角度为调节人体代谢、预防糖尿病等疾病提供科学依据。
SPF级雄性SD大鼠购自上海斯莱克实验动物有限责任公司,实验动物使用许可证编号:SYXK(闽)2016-0005。动物实验在福建省福建中医药大学实验室进行,实验经福建省中医科学院动物伦理委员会批准,批准编号为FJATCM-IAEC2021002。
大鼠饲料 北京华阜康生物科技股份有限公司;新鲜莲子 福建绿地有限公司;SL(纯度≥98%)上海阿拉丁生化科技有限公司;生理盐水 江西科伦药业有限公司;无水乙醇 上海国药集团化学试剂有限公司;异丙醇 美国Merck公司;苯丙氨酸 上海阿达玛斯试剂有限公司;AxyPrepDNA凝胶回收试剂盒美国AXYGEN公司。
Centrifuge 5430R型高速冷冻离心机 德国Eppendorf公司;QuantiFluor-ST微型荧光计 美国Promega公司;Vanquish Horizon system超高效液相色谱系统串联Q-Exactive HF-X复合四极杆-轨道阱质谱仪美国Thermo Fisher Scientific公司;HSS T3色谱柱(100 mm×2.1 mm,1.8 μm) 美国Waters公司;JXDC-20型氮气吹扫仪 上海净信实业发展有限公司。
1.3.1 LRS制备和动物分组
LRS参考文献[14]制备。24 只体质量(70±5)g SPF级雄性SD大鼠饲养在标准的聚丙烯笼中,动物房温度(22±2)℃、相对湿度(55±5)%、12 h昼夜循环。将大鼠随机分为4 组:正常对照(normal control,NC)组、LRS、SL组以及LRS+SL组。依照文献[15]的方法进行灌胃,NC组灌胃2.0 mL/100 g的生理盐水;LRS组灌胃0.25 g/100 g的III型LRS,SL组灌胃0.018 g/100 gSL;LRS+SL组大鼠每100 g体质量灌胃0.25 g的LRS和0.018 g的SL,各组自由摄食饮水干预4 周,每周记录体质量。对大鼠按剂量40 mg/kg进行腹腔注射戊巴比妥钠,麻醉处死后切开腹腔,取3~5 cm回肠肠段,将内容物挤入15 mL无菌离心管中。将含有大鼠小肠内容物的离心管置于液氮中,全部处理结束后转移置于-80 ℃冰箱中冻存。
1.3.2 高通量测序
委托上海美吉生物医药科技有限公司完成高通量测序,流程如图1所示。采用DNA提取试剂盒并按照操作说明对大鼠小肠内容物中的总DNA进行提取,并用1%(质量分数)琼脂糖凝胶电泳检测,按指定测序区域(16S V3~V4区)合成带有barcode的特异引物。随后取特异引物进行聚合酶链式反应(polymerase chain reaction,PCR)扩增,随机选取具有代表性的样本进行预实验,确保在最低循环数中使绝大多数样本能够扩增出浓度合适的产物。混合同一样本的PCR产物并用2%(质量分数)琼脂糖凝胶电泳检测,使用AxyPrepDNA凝胶回收试剂盒切胶回收PCR产物,用Tris-HCl洗脱,再用2%(质量分数)琼脂糖电泳初步定量。将PCR产物用QuantiFluor™-ST蓝色荧光定量系统定量,之后按照每个样本的测序量要求,进行相应比例的混合。最后进行Miseq文库构建和Miseq测序并计算Shannon指数和Chao指数。
图1 测序实验流程Fig. 1 Flow chart of gene sequencing
1.3.3 超高效液相色谱串联复合四极杆-轨道阱质谱分析
参照文献[15]的方法,取1.3.1节中所得样品50 mg,置于2 mL离心管中,加入400 μL含0.02 mg/mL-2-氯苯丙氨酸(内标)溶液(溶剂为甲醇-水(体积比4∶1)溶液)。然后-10 ℃下用冷冻组织研磨机50 Hz研磨6 min,接着低温超声(5 ℃、40 kHz)30 min。样品于-20 ℃条件下放置30 min后,4 ℃、13 000×离心15 min,收集上清液备用,并经0.22 μm滤膜过滤。采用Vanquish Horizon system超高效液相色谱系统串联Q-Exactive HF-X复合四极杆-轨道阱质谱仪进行分析。色谱条件:HSS T3色谱柱(100 mm×2.1 mm,1.8 μm);流动相A为95%水+5%乙腈(含体积分数0.1%甲酸),流动相B为47.5%乙腈+47.5%异丙醇+5%水(含体积分数0.1%甲酸);进样量2 μL;柱温40 ℃。质谱条件:样品经电喷雾电离,分别采用正、负离子扫描模式采集质谱信号。每个样品取20 μL上清液,混匀后作为质控(quality control,QC)样本。在分析的过程中,每检测5~15 个分析样本后测定1 个QC样本,以考察整个检测过程的稳定性。代谢组学分析得到各化合物的保留时间、质荷比和峰强度的数据矩阵,同时将质谱和串联质谱信息与人类代谢组数据库(human metabolome database,HMDB)METLIN数据库以及美吉自建数据库进行比对分析。
实验设置3 个平行,结果用平均值±标准差表示。采用SPSS 20.0软件进行单因素方差分析,采用Tukey-Kramer检验进行显著性分析,<0.05表示差异显著。采用OriginPro 8.5软件和GraphPad Prism 8.0软件作图。物种组间重叠数据采用fastp、FLASH软件进行主坐标分析(principal co-ordinates analysis,PCoA)。采用R-3.3.1(stat)软件分析组间物种差异并绘制Venn图,通过假设检验评估差异的显著性水平。采用python-2.7软件进行菌群群落分析、群落Circos分析并作图;采用python-2.7软件、R-3.3.1(stat)软件进行菌群聚类分析。代谢谱差异数据采用ropls(R packages)软件进行正交偏最小方差判别分析(orthogonal partial least squares-discriminant analysis,OPLS-DA)。差异代谢物差异倍数(fold change,FC)和变量投影重要性值(variable projection importance value,VIP)采用ropls(R packages)、scipy(Python)软件计算。采用ropls(R packages)软件计算某一代谢物在两组间的FC。筛选VIP>1,<0.05、FC<1或FC>1的差异代谢物与京都基因与基因组百科全书(Kyoto Encyclopedia of Genes and Genomes,KEGG)数据库进行比对分析,功能通路数据采用scipy(Python)软件进行拓扑学分析。采用scipy(Python)软件对菌群多样性与代谢谱进行皮尔森相关性分析。
为描述肠道菌群多样性和丰富度,分别计算Shannon指数和Chao指数,结果表明,相较于NC组,LRS组(<0.05)、SL组(<0.01)、LRS+SL组(<0.05)物种丰富度均有显著提高,且LRS组多样性降低(<0.05),这与林珊的研究结果一致,SL组无显著变化(>0.05),LRS+SL组物种多样性显著增加(<0.05)(图2A、B),且LRS+SL组相比LRS组(<0.001)及SL组(<0.05)物种多样性均显著增加,说明LRS和SL协同有利于增加大鼠小肠菌群的多样性。各组样本的PCoA分析如图2C所示,该差异检验结果为0.021,样本间差异显著。NC组与LRS组和SL组间重叠少,与LRS+SL组重叠较多,差异性小。如图2D所示,4 组中共有的物种数为118 个,NC组有20 个特有物种,LRS组有33 个特有物种,SL组有64 个特有物种,LRS+SL组有145 个特有物种。结果表明,当添加LRS和SL时特有物种数均有增加,而在LRS与SL协同作用时增加效果最为显著。
图2 LRS、SL干预对大鼠肠道菌群多样性的影响(n=6)Fig. 2 Effects of LRS and SL intervention on intestinal flora diversity in rats (n = 6)
在不同分类学水平上对物种组成进行分析,如图3A所示,从门水平来看,厚壁菌门(Firmicutes)、放线菌门(Actinobacteriota)和变形菌门(Proteobacteria)是各组中的优势菌群,相比NC组,LRS组Firmicutes相对丰度增加,Actinobacteriota相对丰度减小;SL组Firmicutes相对丰度减小;LRS+SL组拟杆菌门(Bacteroidota)相对丰度增加。其中,Firmicutes代谢产物丁酸盐具有抗炎、改善胰岛素抵抗的作用,其相对丰度增加有利于预防肥胖和2型糖尿病的发生发展,Firmicutes利用或分解乳酸的能力较弱,在SL组中相对丰度下降,这与Policastro等研究结果一致。科水平以乳酸杆菌科(Lactobacillaceae)为主要物种,其次是棒状杆菌科(Corynebacteriaceae)和消化链球菌科(Peptostreptococcaceae),相较于NC组,LRS组乳酸杆菌科、消化链球菌科相对丰度上调,且棒状杆菌科相对丰度下调,SL组和LRS+SL组均表现为消化链球菌科相对丰度下调(图3B)。属水平上,乳酸杆菌属()相对丰度最高,棒状杆菌属()、次之,LRS组上调而相对丰度下调,LRS组相对丰度增加,在其余两组相对丰度降低(图3C)。对图3C中相对丰度前15的物种进行筛选,其中具有显著差异的物种如图3D所示,结合图3E和图3F可知,在LRS组中的相对丰度低于其余3 组,漫游球菌属()在LRS组、SL组和LRS+SL组中的水平均极显著低于NC组(<0.01)。LRS+SL组相较其余各组表现出显著促进生长(<0.01)和抑制异杆菌属()繁殖(<0.05)的效果。其中,被认为是保护肠道屏障、调控免疫应答、抗病原菌感染、调节代谢、预防炎症的益生菌;可通过产生天然类胡萝卜素(角黄素)发挥抗氧化、预防癌症的作用;与人类肠道感染疾病有关,类似于肠杆菌属()和肉食杆菌属()等致病菌。
图3 群落组成及多组间差异分析结果Fig. 3 Analysis of intestinal bacterial community composition and differential species among groups
在属水平上,与NC组相比,LRS组(<0.01)、(<0.05)、(<0.05)、相对丰度显著或极显著降低,、相对丰度显著增加(<0.05)(图4A)。有研究表明,被认为是人体肠道致病菌之一,其参与了肠道内未消化蛋白的发酵,可能破坏肠道稳态,其产物氨、腐胺等与结直肠癌的发生有密切联系;被证明在病态肥胖个体肠道菌群中相对丰度较高,可能具有潜在的致代谢紊乱作用。SL组相对丰度显著降低(<0.05),、相对丰度显著增加(<0.05)(图4B)。LRS+SL组、相对丰度极显著降低(<0.01),(<0.01)、(<0.05)相对丰度极显著或显著增加(图4C)。此外,Lachnospiraceae和Ruminococcaceae同属Firmicutes,是人类肠道中常见的菌群,能水解淀粉和糖类从而产生丁酸盐及短链脂肪酸。同时,LRS组和SL组相对LRS+SL组均表现出显著较低的相对丰度(<0.05)和极显著较高的相对丰度(<0.01)(图4D、E)。的相对丰度与ANGPTL4(即一种脂质代谢的关键调节因子,也是肠道微生物群和脂肪沉积的循环介质)的表达呈正相关,且有实验表明,高脂饮食处理也使小鼠相对丰度显著增加。
图4 两组间群落组成差异分析结果Fig. 4 Analysis of differences in intestinal bacterial community composition between two groups
由此可见,LRS和SL的饮食干预有利于小肠菌群结构的优化、维持肠道稳态,这与Zeng Hongliang等研究结果具有一致性。另外,LRS和SL协同作用促进了的增殖并抑制了的增殖,具有预防代谢紊乱的潜在功效。
代谢谱数据的OPLS-DA结果如图5所示。在阴离子模式下,LRS组与NC组、LRS+SL组与NC组、LRS组与LRS+SL组、SL组与LRS+SL组代谢模式均可完全分离,表明LRS及LRS与SL协同作用均可有效干预代谢模式的变化。
图5 组间代谢差异OPLS-DA评分图Fig. 5 OPLS-DA score plots for differential metabolites between groups
在VIP>1,<0.05、FC<1或FC>1的条件下,将所筛选出的差异代谢物与KEGG数据库比对分类,结果如图6所示,具体差异代谢物见表1。
图6 组间KEGG差异代谢物分析结果Fig. 6 KEGG differential metabolite analysis between groups
由表1可知,相对于NC组,LRS极显著或显著地促进了()-硫辛酸(<0.01)、-姜黄素(<0.05)、17,18-乙醇酰胺(<0.05)、银杏内酯A(<0.01)、17--雌二醇3-硫酸盐-17-(--葡萄糖醛酸)(<0.01)、油酰乙醇酰胺(<0.05)、类固醇激素(雌三醇、雌激素、皮质甾酮)(<0.05)等代谢物的上调,以及磷脂(LysoPC(17:0))、氨基酸(-丝氨酸)等的显著下调(<0.05),表明LRS对体内脂类、氨基酸代谢有调节作用。其中,-姜黄素是改善脂质代谢、保护肠道健康,降低血液胆固醇水平的重要代谢物。LysoPC(17:0)可能通过氧化应激在机体损伤过程中发挥重要作用。-丝氨酸为非必需氨基酸,丝氨酸的从头合成在提高机体免疫功能方面有积极影响。在SL作用下,大鼠体内磷脂代谢有明显变化(LysoPC(18:1(9))、LysoPC(20:2(11,14))显著下调(<0.05),LysoPC(24:1(15))显著上调(<0.05)),油酰乙醇酰胺、类固醇激素(脱氢表雄酮、皮质醇、糖皮质类固醇激素、皮质甾酮、雄烯二酮)水平显著上调(<0.05)。LRS与SL协同作用时,大鼠体内氨基酸代谢变化(-组氨酸、-色氨酸下调,瓜氨酸上调)显著(<0.05),同时表现出3a,7a-二羟基-5b-胆甾烷、维生素(泛酸)水平显著上调(<0.05)和磷脂(LysoPC(17:0))水平下调。此外,与LRS单独作用相比,该协同作用对维生素代谢的调节(烟酸、抗坏血酸显著下调(<0.05))有更显著的影响;相较于SL单独作用,二者协同对17-羟基亚麻酸、LysoPC(16:1(9)/0:0)水平有上调作用,对脱氧雪腐镰刀菌烯醇、山柰酚等代谢物有下调作用。其中,瓜氨酸在机体氮素循环中发挥重要作用,在体内吸收后可转化为精氨酸,具有提高机体的免疫力、保持正常的胆固醇水平、平衡正常的血糖水平、清除自由基等作用。脱氧雪腐镰刀菌烯醇是一种引起机体免疫抑制和肠道菌群紊乱的毒素,对肠道屏障功能有破坏作用。综上所述,LRS和SL均有利于类固醇激素的合成,且LRS单独作用时可促进大鼠体内()-硫辛酸、银杏内酯A、-姜黄素等物质水平的增加,有利于调节氨基酸和脂质代谢;SL干预可调节类固醇激素及磷脂等代谢物水平;二者共同作用时,对氨基酸代谢和维生素代谢的调节作用更为明显,并可能提高机体对有害代谢物的抵抗作用。
表1 LRS组、SL组、LRS+SL组与NC组及LRS组、SL组与LRS+SL组小肠差异代谢物对比Table 1 Differential metabolites in small intestine of LRS, SL and LRS + SL versus NC groups, and LRS and SL versus LRS + SL groups
续表1
LRS组、SL组、LRS+SL组与NC组间及LRS组、SL组与LRS+SL组之间差异显著的功能通路如图7所示。其中,横轴表示通路中代谢物在通路中的相对重要性,气泡大小代表通路重要性值。NC组与LRS组有2 条差异较大的功能通路:类固醇激素(steroid hormone)、半胱氨酸和蛋氨酸代谢(cysteine and methionine metabolism)通路(图7A)。类固醇激素功能通路与雌激素、皮质甾酮、脱氢表雄酮等类固醇及其衍生物有关,并与一种调节胆汁酸代谢及脂质代谢的细菌丰度呈正相关。NC与SL组有4 条差异较明显的功能通路:类固醇降解(steroid degradation)、类固醇激素生物合成(steroid hormone biosynthesis)、初级胆汁酸合成(primary bile acid synthesis)、谷胱甘肽代谢(glutathione metabolism)通路(图7B)。其中,初级胆汁酸合成对抗胆汁淤积有重要作用。谷胱甘肽代谢是机体内重要的抗氧化和代谢调节的途径,能够清除氮氧自由基、预防肝脏疾病。对照组与LRS+SL组有3 条差异较明显的功能通路:氨酰tRNA生物合成(aminoacyl-tRNA biosynthesis)、-丙氨酸代谢(-alanine metabolism)、星孢菌素生物合成(staurosporine biosynthesis)通路(图7C)。LRS组与LRS+SL组间有2 条差异较大的功能通路:谷胱甘肽代谢(Glutathione metabolism)、抗坏血酸和醛酸代谢(ascorbate and aldarate metabolism)(图7D)。SL组与LRS+SL组间有3 条差异较大的功能通路:黄酮和黄酮醇(flavone and flavonol)、类黄酮生物合成(flavonoid biosynthesis)和甘油磷脂代谢(glycerophospholipid metabolism)通路(图7E)。-丙氨酸代谢等氨基酸代谢通路及甘油磷脂代谢等脂质代谢通路均与预防代谢综合征等疾病有关。综上,LRS及SL单独作用时均表现出氨基酸代谢和激素代谢调节的重要作用,SL还在优化肝功能、增强肝防御系统方面有促进作用;二者协同作用时有利于维持脂质代谢、维生素代谢稳定并促进生物活性物质的合成和表达。
图7 组间功能通路差异拓扑学分析结果Fig. 7 Topological analysis of differential functional pathways between groups
图8 肠道菌群-代谢相关性分析结果Fig. 8 Correlation analysis between intestinal flora and metabolism
对小肠菌群差异菌属与差异代谢物进行相关性分析,所得结果如图8所示,球链菌属()、罗斯氏菌属()、红球菌()、、、不动杆菌()、与LysoPC(24:1(15))(<0.05)、去氢催吐萝鞭木醇(<0.001)、辣椒素、油酰乙醇酰胺、皮质甾酮、雌激素、()-硫辛酸、雌三醇、银杏内酯A(<0.05)呈显著正相关,与-丝氨酸、降胭脂木酯、LysoPC(20:2(11,14))、LysoPC(18:1(9))呈显著负相关(<0.05)。与LysoPC(24:1(15))(<0.05)、去氢催吐萝鞭木醇(<0.001)、辣椒素(<0.05)、油酰乙醇酰胺(<0.01)、皮质甾酮(<0.01)呈显著正相关。与雌激素、()-硫辛酸、银杏内酯A呈显著负相关(<0.01)。与瓜氨酸呈显著负相关(<0.01),与LysoPC(17:0)、组氨酸等维生素呈显著正相关(<0.05)。其中,银杏内酯A参与血清和肝组织TC与TG含量的调控,降低肝细胞脂肪变性和脂质蓄积程度,有效预防非酒精性脂肪性肝病。有研究表明,()-硫辛酸能调节线粒体功能和脂肪细胞胰岛素敏感性,是改善胰岛素抵抗和相关代谢紊乱的潜在有效营养素。肝脏是雌激素代谢的主要部位,也是雌激素作用的主要非生殖靶器官;雌激素调节肝脏蛋白质的合成、肝脏能量稳态,参与脂质和葡萄糖代谢的肝脏基因以及肝细胞的增殖和凋亡,对肝功能有重要调节作用。辣椒素是从辣椒中提取的主要活性成分,研究表明其不仅具有抗氧化、清除活性氧自由基和改善线粒体活性等功能,还在改善肝细胞膨大、气球样变及脂肪空泡等病理损伤方面发挥积极作用。油酰乙醇酰胺作为一种天然存在的脂质信号分子,主要通过改善脂肪酸诱导的氧化应激,预防非酒精性脂肪性肝病。与银杏内酯A呈正相关的可以降解脂肪酶、腈水合酶和胆固醇氧化酶为丁酸、丙酸等短链脂肪酸,对促进脂肪氧化、改善代谢健康有重要作用。也是一种产丁酸的细菌,在调节肠道能量代谢、防止肠道微生态失衡方面有重要作用。结合小肠菌群和代谢谱差异可得,LRS单独作用时不仅有利于促进等有益菌的增殖,增加银杏内酯A、()-硫辛酸、油酰乙醇酰胺、辣椒素等有益代谢物的产生,促进和产短链脂肪酸细菌的繁殖,降低等有害菌相对丰度,起到调节氨基酸、激素、脂类代谢的作用。LRS和SL协同作用时,可通过降低相对丰度,调节维生素代谢,可能对维持体内代谢平衡、预防代谢综合征有潜在效果。
本研究探讨了LRS及SL对大鼠的小肠道菌群结构及代谢的影响。通过小肠菌群结构、代谢物组成及其相关性的分析,结果表明,LRS的饮食干预促进了的增殖并降低了相对丰度,与辣椒素、油酰乙醇酰胺、雌激素、()-硫辛酸、银杏内酯A正相关,影响了类固醇激素和半胱氨酸和蛋氨酸代谢代谢通路,可能在调节氨基酸、激素代谢等方面发挥积极作用;SL的干预提高了、的相对丰度,影响了初级胆汁酸合成及谷胱甘肽代谢通路,具有预防胆汁酸堆积、增强肝功能的潜在效果;当二者协同作用时,能够增加小肠菌群多样性,同时促进的增殖并降低的相对丰度,与瓜氨酸呈负相关并与LysoPC(17:0)呈正相关,能够富集维生素、脂质、氨基酸代谢通路,与LRS和SL单独作用相比,其在维持机体稳态方面有着更为显著的效果。研究结果可为益生元及其后生元协同效应的研究提供一定的理论依据。