參附注射液對阿霉素誘導(dǎo)的心力衰竭大鼠腸道細菌多樣性的影響

趙震宇， 胡炎芝， 董振華， 劉嘉豪， 張曦寧， 李 琳,2*

(1.湖南中醫(yī)藥大學(xué)，湖南 長沙 410208；2.湖南中醫(yī)藥大學(xué)中醫(yī)診斷研究所，湖南 長沙 410208)

心力衰竭是一種復(fù)雜的臨床綜合征，是許多心血管疾病發(fā)展的晚期[1]。很多疾病可以引起心衰，包括心肌病、心臟瓣膜病、糖尿病、肥胖等[2-3]。據(jù)統(tǒng)計，全球約6 300萬人患有此病，約50%的患者在最初診斷后的5年內(nèi)死亡，超過了大多數(shù)類型的癌癥，已成為一種全球健康的負擔(dān)[4]。近年來心力衰竭以西醫(yī)治療為主，但副作用明顯，如應(yīng)用利尿劑引起的低鈉血癥常導(dǎo)致神志不清、跌倒和痙攣等虛弱癥狀[5]。鑒于目前西醫(yī)治療心衰的困境，人們開始在中醫(yī)藥領(lǐng)域?qū)で蟆敖馑帯?，包括參附注射液在?nèi)的多種中藥注射劑已經(jīng)廣泛用于臨床，療效明確[6]。近年來許多心力衰竭與腸道菌群的相關(guān)研究證實，腸道菌群可通過介導(dǎo)三甲胺(TMA)/三甲胺氧化氮(TMAO)、短鏈脂肪酸、膽汁酸等代謝途徑參與心力衰竭的發(fā)生與發(fā)展[7]。參附注射液的治療效果是否與調(diào)節(jié)腸道細菌菌群結(jié)構(gòu)有關(guān)，目前尚未明確。因此，本實驗以腸道菌群為切入點，利用16S rRNA高通量測序的方法，明確參附注射液對心衰大鼠腸道菌群的影響，探討參附注射液對心力衰竭大鼠療效的作用機制。

1 材料

1.1 動物 15只5周齡雄性Sprague-Dawley大鼠，由湖南斯萊克景達實驗動物有限公司提供，實驗動物生產(chǎn)許可證號SCXK(湘)2019-0004。實驗方案和動物使用已通過湖南中醫(yī)藥大學(xué)動物實驗管理中心動物使用與管理倫理委員會的批準(zhǔn)，批號LL2020060401，實驗設(shè)計和實施依據(jù)《湖南實驗動物條例》和國家衛(wèi)生研究院動物倫理使用指南進行。實驗開始前，所有大鼠適應(yīng)培養(yǎng)1周。

1.2 試劑、藥物與儀器 參附注射液[華潤三九(雅安)藥業(yè)有限公司，批號20030101005]。阿霉素(美國APExBIO公司，批號A39662337769)；生理鹽水(湖南康源制藥有限公司，批號C20090105d)；DNA抽提試劑盒(美國Bio-Tek公司，型號E.Z.N.A.? soil DNA kit)。MiSeq PE300型高通量測序儀(美國Illumina公司)；PCR儀(美國ABI公司，型號ABI GeneAmp?9700)；建庫試劑盒(美國Bioo Scientific公司，型號NEXTFLEX?Rapid DNA-Seq Kit)。

2 方法

2.1 造模 所有大鼠適應(yīng)性培養(yǎng)1周后，按照隨機數(shù)字表法分為空白組、模型組、參附組，每組5只，其中空白組大鼠每周腹腔注射生理鹽水，其余大鼠每周腹腔注射阿霉素(生理鹽水配制，2 mg/mL)1次，劑量為1.5 mL/kg，連續(xù)7周，總劑量21 mg/kg，依據(jù)文獻[8]報道制備心衰模型。

2.2 給藥 7周后，參附組大鼠按6 mL/kg劑量連續(xù)15 d腹腔注射參附注射液，空白組、模型組大鼠腹腔注射等量生理鹽水[9]。

2.3 樣本收集 末次給藥24 h后收集大鼠糞便，放入凍存管中保存，液氮速凍后放入-80 ℃冰箱中保存?zhèn)溆肹10]。

2.4 DNA提取 依據(jù)E.Z.N.A.?soil DNA kit說明書對收集到的大鼠糞便樣本進行微生物群落總DNA抽提，瓊脂糖凝膠電泳檢測DNA提取質(zhì)量，NanoDrop2000測定DNA濃度和純度。

2.5 腸道細菌基因聚合酶鏈反應(yīng)(PCR)和16S rRNA測序 使用338 bp正向5′-ACTCCTACGGGAGGCAGCAG-3′和806 bp反向 5′-GGACTACHVGGGTWTCTAAT-3′序列對16S rRNA基因V3～V4可變區(qū)進行PCR擴增，參數(shù)為95 ℃ 3 min；95 ℃ 30 s，55 ℃ 30 s，72 ℃ 45 s，27個循環(huán)；72 ℃ 10 min，在4 ℃下保存。PCR產(chǎn)物經(jīng)瓊脂糖凝膠電泳檢測，純化后的產(chǎn)品采用上海美吉生物醫(yī)藥科技有限公司提供的Illumina MiSeq PE300測序系統(tǒng)進行測序。

2.6 測序數(shù)據(jù)處理 ①過濾reads尾部質(zhì)量值20以下的堿基，設(shè)置50 bp的窗口，如果窗口內(nèi)的平均質(zhì)量值低于20，從窗口開始截去后端堿基，過濾質(zhì)控后50 bp以下的reads，去除含N堿基的reads；②根據(jù)PE reads之間的overlap關(guān)系，將成對reads拼接(merge)成1條序列，最小overlap長度為10 bp；③拼接篩選，最大錯配比率為0.2；④修正，barcode允許的錯配數(shù)為0，最大引物錯配數(shù)為2；⑤分析，使用UPARSE 7.1軟件，按照97%的相似度對序列進行OTU聚類并剔除嵌合體。利用RDP classifier對每條序列進行物種分類注釋，比對Silva 16S rRNA數(shù)據(jù)庫(版本號132)，比對閾值是70%。

3 結(jié)果

3.1 參附注射液對心衰大鼠腸道菌群多樣性的影響 15個樣本通過測序、質(zhì)控過濾后，共獲得413 910條有效序列數(shù)目，單個樣本平均產(chǎn)生27 594條，稀疏曲線隨測序量的增加而趨于平緩，一直延伸到X軸右端，說明目前的測序量可用于后續(xù)分析，見圖1。

通過對端測序可知，15個樣本共獲得1 028個OTUs，空白組、模型組、參附組OTU數(shù)目分別是844、871、840，共有OTUs 653個，每組樣本的Chao、ACE、Simpson、Shannon等指數(shù)見表1。由此可知，與空白組比較，模型組Simpson指數(shù)值下降(P<0.01)，Chao、ACE、Shannon指數(shù)值升高，以Shannon值更明顯(P<0.05)，說明心衰后大鼠腸道細菌多樣性增加，細菌豐度也呈現(xiàn)增長趨勢；與模型組比較，參附組Chao、ACE、Simpson、Shannon指數(shù)值回調(diào)，說明參附注射液有助于心衰大鼠腸道細菌豐度和多樣性的恢復(fù)。

表1 各樣本Alpha多樣性分析

3.2 參附注射液對心衰大鼠腸道菌群結(jié)構(gòu)的影響 基于unweighted unifrac距離算法的主坐標(biāo)分析(principal coordinates analysis，PCoA)，考察大鼠腸道菌群β多樣性的差異，見圖2。由此可知，PC1對總體菌群的代表性為27.35%，PC2為17.07%，3組組間差異明顯；基于PC2軸可見，模型組的點與空白組的點偏差較大，空白組與參附組的點距離較近，說明空白組與參附組菌群結(jié)構(gòu)相似度高，參附注射液有助于恢復(fù)心衰大鼠腸道細菌的結(jié)構(gòu)。

3.3 參附注射液的心衰大鼠腸道細菌相對豐度的影響 腸道微生物群落結(jié)構(gòu)的直方圖揭示了微生物種類及其相對豐度，共檢測出7個門，見圖3，可知3個樣本組的優(yōu)勢門主要來源于厚壁菌門(Firmicutes)、擬桿菌門(Bacteroidota)，其次是放線菌門(Actinobacteria)、疣微菌門(Verrucomicrobia)、脫硫桿菌門(Desulfobacterota)、變形菌門(Proteobacteria)、螺旋體門(Spirochaetota)。與空白組比較，模型組擬桿菌門(Bacteroidota)、放線菌門(Actinobacteria)、脫硫桿菌門(Desulfobacterota)、螺旋體門(Spirochaetota)豐度最高，厚壁菌門(Firmicutes)豐度最低，模型組厚壁菌門相對豐度(70.88%)低于空白組(81.06%)(P<0.01)，模型組擬桿菌門相對豐度(15.07%)高于空白組(6.87%)(P<0.05)，參附組變形菌門相對豐度(2.95%)高于模型組(0.16%)(P<0.01)，見圖4～5。與空白組比較，模型組厚壁菌門豐度減小(P<0.01)，擬桿菌門豐度增加(P<0.05)，F(xiàn)/B值減??；與模型組比較，參附組厚壁菌門豐度增加，擬桿菌門豐度減小，F(xiàn)/B值增加，差異無統(tǒng)計學(xué)意義(P>0.05)，可認(rèn)為參附組大鼠腸道中厚壁菌門和擬桿菌門豐度及其比值有回調(diào)趨勢，見表2。

所有標(biāo)本共鑒定出30個屬，見圖6。由此可知，空白組的優(yōu)勢菌依次為乳桿菌屬Lactobacillus39.78%，羅姆布茨菌屬Romboutsia7.28%，未知的S24-7科6.36%，未知的梭狀芽胞桿菌UCG-014家族5.85%，艾克曼菌屬Akkermansia5.83%，未分類的毛螺旋菌科unclassifiedLachnospiraceae4.13%；模型組的優(yōu)勢菌依次為乳桿菌屬21.91%，未知的S24-7科11.63%，羅姆布茨菌屬5.85%，雙歧桿菌屬Bifidobacterium5.12%，未分類的毛螺旋菌科unclassifiedLachnospiraceae4.09%；參附組的優(yōu)勢菌依次為異桿菌屬14.98%，未知的S24-7科12.84%，乳桿菌屬9.83%，Turicibacter7.27%，未知的梭狀芽胞桿菌UCG-014家族6.36%，即空白組和模型組乳桿菌屬比例最大，參附組異桿菌屬為優(yōu)勢菌屬。

本實驗根據(jù)3組樣本在屬水平的物種注釋及豐度信息，選取豐度排名前35的屬，根據(jù)其在每組樣本中的豐度信息按組別進行聚類，見圖7。由此可知，這35個屬的細菌分屬于6個優(yōu)勢的細菌門，其中26個屬于厚壁菌門，3個屬于放線菌門，1個屬于疣微菌門，2個屬于擬桿菌門，2個屬于脫硫桿菌門，1個屬于螺旋體門，表明在屬水平上檢測到的細菌中通過相對豐度進行了明顯的分離。

然后，篩選出具有維持腸道結(jié)構(gòu)完整性和(或)改善心衰作用的7個細菌屬，并進行組間差異檢驗，見圖8～9。由此可知，與空白組比較，模型組羅姆布茨菌屬、布勞特氏菌屬Blautia、異桿菌屬、糞球帶菌屬Faecalibaculum的相對豐度降低，密螺旋體菌屬Treponema、Intestinimonas、普雷沃氏菌屬Prevotella的相對豐度升高，以Intestinimonas、普雷沃氏菌屬更明顯(P<0.05)；與模型組比較，參附組羅姆布茨菌屬、布勞特氏菌屬、異桿菌屬、糞球帶菌屬、密螺旋體菌屬、Intestinimonas、普雷沃氏菌屬的相對豐度均呈現(xiàn)回調(diào)趨勢，以異桿菌屬(P<0.05)、Intestinimonas(P<0.01)更明顯，說明參附注射液可有效恢復(fù)菌群豐度的改變。

3.4 線性判別分析效應(yīng)量(LEfSe)確定豐度差異顯著的菌群 本實驗篩選出3組大鼠腸道中具有差異的微生物(LDA閾值為4)，采用LEfSe分析得出多級物種分類樹圖和柱狀圖，見圖10。由此可知，空白組腸道菌群主要包括4種細菌，即乳桿菌屬Lactobacillus、乳桿菌科Lactobacillaceae、乳桿菌目Lactobacillales、厚壁菌門Firmicutes，以益生菌為主；模型組中富含1種細菌，即普雷沃氏菌科Prevotellaceae，為條件致病菌；參附組特征性菌群有9種，即丹毒絲菌科Erysipelotrichaceae、丹毒絲菌目Erysipelotrichales、異桿菌屬Allobaculum、蘇黎世桿菌屬Turicibacter、嚴(yán)格梭菌屬Clostridium_sensu_stricto_1、梭菌科Clostridiaceae、梭菌目Clostridiales、γ-變形菌綱Gammaproteobacteria、變形菌門Proteobacteria。

4 討論

腸道是人體最大的消化器官，與人類的生命活動密切相關(guān)[11]。生活在腸道內(nèi)的腸道菌群在人類日常飲食和口服藥物中扮演著“代謝物過濾器”的重要角色[12]，其對機體物質(zhì)的吸收、消化、代謝功能與中醫(yī)小腸“受盛化物”的功能相契合，對有害物質(zhì)進行識別、抵抗、吸收的功能與中醫(yī)小腸“泌別清濁”的功能相關(guān)[13]，可認(rèn)為腸道菌群是腸道和其他器官之間的一個重要橋梁。

腸道菌群微生態(tài)環(huán)境的動態(tài)平衡才能維持正常的腸道微生態(tài)功能，如免疫功能、營養(yǎng)作用等，一旦平衡被打破，便會伴發(fā)多種疾病，如肥胖、胰島素抵抗、糖尿病等[14-15]。已有研究證實，心衰發(fā)生發(fā)展過程中可引起腸道黏膜損傷，腸道菌群移位，引發(fā)全身炎癥反應(yīng)和免疫反應(yīng)，失調(diào)的菌群通過介導(dǎo)炎癥反應(yīng)進一步影響心衰的病理生理過程，誘發(fā)心衰晚期惡病質(zhì)等不良反應(yīng)[16-17]。大量臨床研究證實，參附注射液治療后心衰患者心功能明顯改善，生存率明顯提高[18-19]。結(jié)合中醫(yī)對腸道功能的認(rèn)知和腸道菌群功能的重要性，推測腸道菌群可能是參附注射液治療心力衰竭的關(guān)鍵，故本實驗采用16S rRNA高通量測序技術(shù)探究參附注射液對心力衰竭大鼠腸道菌群結(jié)構(gòu)的影響。

本研究中，與空白組比較，模型組的OTUs數(shù)量增加，結(jié)合Alpha多樣性分析結(jié)果，可認(rèn)為心衰大鼠腸道細菌的豐度和多樣性較空白組增加，其腸道菌群豐度和多樣性的增加可能是由于各種有害細菌的過度繁殖導(dǎo)致的，這種觀點也被以往研究認(rèn)同[20]。給予參附注射液后，參附組的OTUs數(shù)量減少，Alpha多樣性分析指數(shù)回調(diào)，與模型組比較無統(tǒng)計學(xué)意義，此結(jié)果與本課題組前期研究結(jié)果類似[21]。結(jié)合PCoA圖結(jié)果，PC2軸空白組與參附組的距離較空白組與模型組的距離明顯縮短。提示，參附注射液有助于恢復(fù)心衰大鼠腸道細菌的豐度和多樣性，使菌群結(jié)構(gòu)接近空白組水平。

厚壁菌門是優(yōu)勢菌門，其次是擬桿菌門、放線菌門等[22]。厚壁菌門和擬桿菌門的比值(F/B值)常用來測量腸道內(nèi)穩(wěn)態(tài)[23]。研究發(fā)現(xiàn)，與空白組比較，模型組厚壁菌門相對豐度下降，擬桿菌門相對豐度上升，F(xiàn)/B值減小，說明心衰后大鼠腸道內(nèi)環(huán)境紊亂。給藥后，參附組厚壁菌門相對豐度上升，擬桿菌門相對豐度下降，F(xiàn)/B值增大，與以往研究相符，說明參附注射液對腸道內(nèi)穩(wěn)態(tài)的維持具有積極作用[24]。He等[25]發(fā)現(xiàn)，厚壁菌門和變形菌門是產(chǎn)生乳糖酶的主要菌門，模型組變形菌門的相對豐度低于空白組，參附注射液治療后變形菌門的相對豐度增加，推測參附注射液會促進某些產(chǎn)乳糖酶細菌的增殖來維持腸道內(nèi)環(huán)境穩(wěn)定，改善心衰狀態(tài)。

從屬水平分析3組大鼠腸道菌群結(jié)構(gòu)差異，結(jié)果顯示，3組大鼠具有特征性的高豐度細菌，說明不同組別大鼠腸道菌群的物種結(jié)構(gòu)具有差異性。乳桿菌屬在空白組大鼠腸道菌群中所占比重最大，心衰大鼠乳桿菌屬相對豐度較空白組下降，給藥后乳桿菌屬豐度進一步下降。乳桿菌屬是益生菌，有增強機體免疫力和抑制病原菌侵襲的作用，但給予參附注射液后的大鼠腸道菌群中乳桿菌屬豐度進一步下降，說明參附注射液不能通過增加乳桿菌屬的數(shù)量達到修復(fù)腸黏膜、改善心衰的作用[23,26]。與空白組比較，模型組羅姆布茨菌屬[27]、布勞特氏菌屬[28]、異桿菌屬[29]、糞球帶菌屬[30]等有益菌豐度減少，以往研究表明這些細菌在維持腸道健康方面有重要作用，密螺旋體屬[31]、Intestinimonas[32]、普雷沃氏菌屬[33]等致病菌屬在模型組的豐度增加。給藥后上述細菌豐度回調(diào)，說明參附注射液可能是通過增加有益菌數(shù)量，減少致病菌數(shù)量來維持腸道的完整性。Mao等[30]研究發(fā)現(xiàn)糞球帶菌屬可通過發(fā)酵作用促進丁酸的生成。丁酸是短鏈脂肪酸的一種，近期研究表明，腸道微生物介導(dǎo)的短鏈脂肪酸代謝途徑與心衰有關(guān)，可通過多種機制抑制炎癥反應(yīng)，在改善心衰發(fā)生發(fā)展方面起到重要作用[26]，提示參附注射液可能是通過促進糞球帶菌屬等產(chǎn)短鏈脂肪酸細菌的增殖來改善大鼠心衰狀態(tài)。

LEfSe結(jié)果表明，空白組的特征性細菌是以乳桿菌屬為主的益生菌；模型組以致病菌普雷沃氏菌科為主；給予參附注射液后，普雷沃氏菌科豐度下降，異桿菌屬等有益菌成為參附組的特征性細菌，可認(rèn)為參附注射液具有促進心衰大鼠腸道有益菌增殖，抑制致病菌增殖的作用。

綜上所述，參附注射液對心力衰竭大鼠腸道微生物環(huán)境有積極作用，可改善因心衰導(dǎo)致的大鼠腸道細菌結(jié)構(gòu)紊亂，具有抑制致病菌生長，促進有益菌生長的作用，其機制可能與增加SCFA的產(chǎn)生菌(如糞球帶菌屬)有關(guān)。