基于網絡藥理學聯合16S rDNA高通量測序技術分析丹參對感染大腸桿菌小鼠腸道菌群的影響

王 樂，陳泓岑，張永紅，吳 瓊，侯佳佳，王天祎，盧天航，黃傳發，張 華，崔德鳳

(北京農學院動物科學技術學院/獸醫學(中醫藥)北京市重點實驗室，北京 102206)

隨著畜牧業和飼料工業的發展，飼料添加劑的發展進入了一個新的階段，中草藥飼料添加劑以其能預防動物疾病、提高畜禽產品產量和繁殖性能等優點而日益受到關注，應用前景極為廣闊。相對于抗生素而言，中草藥具有促進動物生長，提高生產性能，改善動物機體免疫，抗氧化、抗炎、抗菌，減少有害氣體排放，降低抗性基因蓄積及轉移等作用，可作為抗生素替代產品廣泛應用于畜禽健康養殖。不論從健康角度還是從經濟效益角度來看，替代抗生素產品的生產和研發顯得尤為重要，而中草藥以其獨有的優勢成為眾多學者的優先選擇。因此，開發安全性高的中藥飼料添加劑代替原有的抗生素及化學藥物將是大勢所趨。

丹參又名赤參、紫丹參、紅根等，為唇形科植物丹參的干燥根和根狀莖。味苦，微寒，歸心、肝經，有活血祛瘀、通經止痛、清心除煩、涼血消癰的功效，用于胸痹心痛，脘腹脅痛，癥瘕積聚，熱痹疼痛，心煩不眠，月經不調，痛經經閉，瘡瘍腫痛。據報道，丹參的化學成分主要分為兩大類：脂溶性二萜醌類化合物和水溶性酚酸類成分。其中脂溶性包括丹參酮Ⅰ、丹參酮Ⅱ A、丹參酮Ⅱ B、二氫丹參酮、隱丹參酮等，具有抗氧化、消炎、抑菌等作用；而水溶性成分包括丹參素、丹酚酸等，具有抑制血小板凝聚、擴張冠狀動脈等作用。李曉明等研究表明，在衣霉素(1 mg·kg) 灌胃復制小鼠急性肝損傷模型中，經丹參酮 Ⅱ 處理后，其可通過抑制肝組織Grp78/caspase-12 通路活化，降低內質網應激水平，改善衣霉素所致的小鼠急性肝損傷。張東濤等用丹參水提物治療乳鼠心肌細胞缺氧-復氧損傷，發現丹參水提物對缺氧-復氧損傷的心肌細胞具有保護作用，這可能與 PI3K/Akt 信號通路的激活有關。由此可見，丹參對緩解應激具有良好的作用。

然而，目前關于丹參是否對腸道菌群失衡具有調節作用，尚未見報道。因此，本研究以丹參提取物為研究對象，應用網絡藥理學與高通量測序16S rDNA技術研究丹參對小鼠大腸桿菌感染模型的腸道菌群結構的影響，以期闡明丹參提取物調節菌群結構的作用機制，為丹參提取物的臨床應用提供理論依據。

1 材料與方法

1.1 丹參活性成分與相應靶點的收集

在TCMSP 數據庫(https://tcmspw.com/tcmsp.php)搜索丹參活性成分，根據TCMCP，口服利用度(oral bioavailability，OB)≥30%、類藥性(drug-likeness，DL)≥0.18表示藥物成分口服吸收良好并且適合于藥物開發，藥物對小腸上皮細胞通透性(Caco-2)≥0.4表示藥物對小腸上皮細胞通透性良好。因此，根據OB≥30%、DL≥ 0.18、Caco-2≥0.4這3個參數對丹參活性成分進行初步篩選。同時收集活性成分所對應的蛋白靶點，利用 STRING 數據庫(https://cn.string-db.org)，轉換靶點蛋白為對應的基因名稱。

1.2 E. coli靶點基因的獲取

在 Genecards數 據 庫(https://www.genecards.org/) 中輸入關鍵詞“”，查詢并收集其靶點基因；將已轉換成基因名的丹參活性成分作用靶點與靶點基因取交集(https://bioinfogp.cnb.csic.es/tools/venny/)，從而得到丹參治療感染的潛在作用靶點。

1.3 丹參治療E. coli感染的“活性成分-潛在作用靶點”網絡構建與分析

將活性成分與對應潛在作用靶點基因導入 Cytoscape 3.8.2軟件，構建丹參治療感染的“活性成分-潛在作用靶點”網絡，并應用軟件中“Network Analyzer”功能進行網絡拓撲結構分析，可獲得活性成分及作用靶點的連接度(Degree)、緊密度(Closenesss)及介度(Betweenness)等網絡拓撲參數。選取丹參成分連接度≥平均數為關鍵活性成分，靶點連接度≥平均值為關鍵靶點。

1.4 丹參成分蛋白互作網絡KEGG通路與GO 功能富集分析

利用DAVID 數 據 庫(https://david.ncifcrf.gov/)，物種選擇“”，對潛在作用靶點進行GO功能分析，利用Metascape數據庫(https://metascape.org/)進行 KEGG 通路富集分析，其中GO功能分析包括生物學過程(Biological process，BP)、分子功能(Molecular function，MF)及細胞成分(Cellular component，CC)分析。

1.5 動物分組與處理

本試驗以安徽地區精品丹參經煎熬回旋蒸發濃縮后獲得的丹參提取藥液為研究對象，以110只體重為(20±5)g的昆明小鼠作為試驗對象，隨機分為5組：對照組(CON)、大腸桿菌感染組()、飼喂低劑量丹參提取物大腸桿菌感染組(EL)、飼喂中劑量丹參提取物大腸桿菌感染組(EM)、飼喂高劑量丹參提取物大腸桿菌感染組(EH)，每組22只。參考人推薦用量的1、5、10倍，設丹參提取物的藥物濃度分別為0.0195、0.039、0.078 g·mL，分別為本次試驗的低、中、高劑量丹參提取物大腸桿菌感染組給藥劑量。小鼠體重以20 g計，按390 mg·kg量進行灌胃，即每次每只試驗鼠灌胃0.2 mL。每日定時定量灌胃1次，連續28 d，各丹參組灌胃不同濃度藥物0.2 mL，對照組與大腸桿菌感染組均灌胃0.2 mL同等劑量的生理鹽水。

1.6 菌種的復蘇

將檢查過菌種純度的O菌種進行傳代復蘇。在滅菌的試管當中進行將該菌種的復蘇，將凍存菌液加入到5 mL左右的營養肉湯液體培養基，置于搖床上，轉速設定成為220 r·min，溫度設定成為37 ℃，培養12 h，傳三代復蘇至鏡下觀察菌形態良好即可。將菌液用比濁管比濁，將菌液濃度大概培養至 6×10CFU·mL。

1.7 大腸桿菌的腹腔注射

本試驗的5組小鼠連續灌胃28 d后，將新鮮培養所得的O按半數致死量對灌胃低、中、高劑量丹參組和大腸桿菌感染組進行腹腔注射0.2 mL·只，而對照組腹腔注射0.2 mL·只的無菌生理鹽水，觀察24 h。

1.8 腸道菌群高通量測序

丹參組和大腸桿菌感染組腹腔注射O及對照組注射無菌生理鹽水24 h后取小鼠的糞便，進行標記，保存于-80 ℃。使用MagPure Stool DNA KF試劑盒B(Magen，中國)，按照制造商的說明書提取糞便樣本的微生物群落DNA。使用Qubit dsDNA BR Assay kit (Invitrogen公司，美國)對DNA進行定量，并用1%瓊脂糖凝膠進行均一檢測。使用引物341 F(5′-ACTCCTACGGGAGGCAGCAG-3′)和806R (5′-GGACTACHVGGGT-WTCTAAT-3′)擴增細菌16S rRNA基因的可變區V3-V4。正向和反向引物均用Illumina接頭、pad和接頭序列標記。在含有30 ng模板、融合PCR引物和PCR主混合物的50 μL反應中進行PCR富集。PCR循環條件：94 ℃ 3 min；94 ℃ 30 s，56 ℃ 45 s，72 ℃ 45 s，30個循環；最后在72 ℃延伸10 min。PCR產物用AmpureXP磁珠純化，在洗脫緩沖液中洗脫。庫通過Agilent 2100生物分析儀(美國Agilent公司)鑒定。使用經過驗證的文庫在Illumina MiSeq平臺(中國深圳BGI)上按照Illumina標準管道進行測序，并產生2×300 bp 的雙端讀數。

1.9 生物信息學分析

美吉生物一站式科研服務平臺(http://www.majorbio.com/)的物種組成分析選項功能中物種Venn圖分析，統計多組或多個樣本中所共有和獨有的物種(如OTU)數目，可以比較直觀地展現不同環境樣本中物種(如OTU)組成相似性及重疊情況；物種組成分析選項功能中群落組成分析，在不同分類學水平上統計各樣本的物種豐度。根據線性判別分析(liner discriminant analysis，LDA)進行LEfSe 多級物種差異判別分析，找出組間在豐度上有顯著差異的物種。基因組功能預測與分析應用PICRUSt(PICRUSt流程存儲了Greengene ID對應的COG信息和KO信息)對OTU豐度表進行標準化，即去除16S marker gene在物種基因組中的copy數目的影響；然后通過每個OTU 對應的Greengene ID，獲得OTU對應的COG家族信息和KEGG Ortholog (KO)信息；并計算各COG的豐度和KO豐度。根據COG數據庫的信息，可以從eggNOG數據庫中解析到各個COG的描述信息，以及其功能信息，從而得到功能豐度譜；根據KEGG數據庫的信息，可以獲得KO、Pathway、EC信息，并能根據OTU豐度計算各功能類別的豐度。此外，針對Pathway，運用PICRUSt可獲得代謝通路的3個水平信息，并分別得到各個水平的豐度表。

1.10 數據分析

數據采用Excel 2021進行初步處理，采用SPSS 21.0軟件ANOVA過程進行單因素方差分析，試驗數據用“平均值±標準差”表示，<0.05為差異顯著，<0.01為差異極顯著。

2 結 果

2.1 丹參活性成分及相應靶點

通過TCMSP平臺，輸入“丹參”搜索。依據設定的OB值、DL值、Caco-2值，從202個化合物中篩選出41個符合條件的化合物(表1)。將符合條件但靶點數為0的化合物，如α-amyrin、sclareol、microstegiol、NSC 122421、tanshinone iia活性成分去除后，最終剩余36個活性成分及669個蛋白質靶點；進入STRING 數據庫(https://cn.string-db.org)，選擇Multiple protein，Lists of Name輸入669個蛋白質靶點，物種選擇“Human”，最終轉換成對應的基因名。

表1 丹參活性成分Table 1 Active components of Radix Salviae Miltiorrhizae

(續表1 Continued)

2.2 丹參治療E. coli感染的潛在作用靶點

從Genecards(https://www.genecards.org/)數據庫中共獲取基因8 902個。同時利用Venny 2.1軟件將丹參活性成分的作用靶點與基因靶點取交集，共獲得51個潛在作用靶點。

2.3 丹參活性成分蛋白互作網絡的構建

將獲得的51個交集靶點導入STRING(https://string-db.org/)數據庫，通過Cytoscape 3.8.2軟件進行蛋白網絡結構可視化處理，并將丹參活性成分與潛在作用靶點進行對應，得到活性成分-潛在作用靶點網絡(圖1)。該網絡包含87個節點，457條邊，87個節點包括 36個活性成分與51個潛在作用靶點，457條邊代表丹參活性成分與潛在作用靶點的相互作用關系，長方形代表丹參活性成分，六邊形代表潛在作用靶點。通過 Cytoscape 3.8.2軟件的內置“network analyze”插件分析網絡的拓撲屬性，得到活性成分及潛在作用靶點的介度、緊密度、連接度。經計算，活性成分的連接度平均數為12.69，故選取連接度≥12.69為關鍵活性成分。關鍵活性成分作用的靶點個數較多，可能是丹參治療感染的重要物質基礎。因此，選取連接度>10為關鍵靶點，共有18個靶點符合，分別為PTGS2、SCN5A、ADRB2、OPRM1、CHRM3、NCOA1、PTGS1、ACHE、CA2、RXRA、OPRD1、ADRA1B、CHRM2、NCOA2、SERPIND1、ESR1、AR、DPP4。

圖1 丹參治療E. coli感染的活性成分-潛在作用靶點網絡Fig.1 Active components-potential action target network of Radix Salviae Miltiorrhizae response to E. coli infection

2.4 丹參成分蛋白互作網絡通路與GO功能富集分析

通過 Metascape數據庫(https://metascape.org/)的51個潛在作用靶點進行 KEGG 通路富集分析，得到<0.01的通路82條，選取<0.01，Count≥6通路可視化(圖2)。同時，找出18個關鍵靶點所對應的通路為關鍵信號通路，共47條。主要包含神經活性配體-受體相互作用、5-羥色胺能突觸、雌激素信號通路、鈣信號通路、cAMP信號通路、PI3K-Akt信號通路、甲狀腺激素信號通路、cGMP-PKG信號通路等途徑。

圖2 KEGG信號通路富集分析Fig.2 KEGG signaling pathway enrichment analysis

通過 DAVID數據庫對丹參治療感染的51個潛在作用靶點進行 GO 富集分析。以<0.05為篩選條件，獲得 BP、MF、CC條目各 174、76、45條。值從小到大排名前10可視化(圖3)。結果顯示，丹參成分作用的BP主要富集在細胞內類固醇激素受體信號通路、突觸囊泡胞吐調節、信號轉導、凋亡過程的正調控、對藥物的反應、對雌二醇的反應、轉錄的陽性調節，DNA模板、細胞增殖調節、核糖核酸聚合酶ⅱ啟動子轉錄正調控、血壓調節；MF主要富集在RNA聚合酶II轉錄因子活性，配體激活的序列特異性DNA結合、酶結合、轉錄因子結合、蛋白結合、相同蛋白結合、血清素結合、類固醇激素受體活性、高分子復合物結合、蛋白激酶結合、類固醇結合；CC主要富集在神經元投射、突觸前膜的組成成分、樹突、高分子復合物、細胞質、質膜的整體部件、突觸后膜的組成成分、軸突、谷氨酸能突觸、質膜。

圖3 GO功能富集分析Fig.3 GO function enrichment analysis

2.5 關鍵成分-關鍵靶點-關鍵信號通路網絡圖

綜合網絡藥理學拓撲參數、KEGG、GO 富集分析所得到的關鍵成分、關鍵靶點、關鍵信號通路，構建丹參治療“關鍵成分-關鍵靶點-關鍵信號通路”網絡(圖4)。結果顯示，關鍵成分為丹參“酮、醌”型結構的脂溶性成分。由圖4可知，丹參的1種活性成分可對應1個或多個靶標，多種活性成分可對應相同的靶標，提示丹參治療感染具有多成分、多靶標的特點。

圖4 丹參治療E. coli感染的“關鍵成分-關鍵靶點-信號通路”網絡Fig.4 Network of “key components-key targets-signaling pathway” of E.coli infection treated with Radix Salviae Miltiorrhizae

2.6 測序數據處理和稀釋曲線

2.6.1 測序數據處理 小鼠糞便樣本中，剔除低質量序列，最終檢出了950 855個有效序列進行后續的分析，去除掉Barcode后，平均長度為460 bp左右。

2.6.2 稀釋曲線 按試驗數據在各個測序深度下檢測到的OTU數目構建稀釋曲線。在曲線的起始段，每條曲線的斜率很大，隨著測序的增加曲線趨于平緩，顯示大部分樣品達到穩定狀態，說明測序的充分性足夠，滿足后續分析要求(圖5)。

圖5 OTU稀釋曲線Fig.5 OTU dilution curve

2.7 菌群相似度分析

通過聚類，將序列按照彼此的相似性分歸為許多小組，每個小組為一個分類操作單元(OTU)，通常在97%的相似水平下的OTU進行生物信息統計分析。由圖6可知，各組之間共有和單獨含有的OTUs數量。其中對照組本身含有519個OTUs，大腸桿菌組本身含有174個OTUs，低劑量丹參提取物大腸桿菌感染組本身含有581個OTUs，中劑量丹參提取物大腸桿菌感染組本身含有565個OTUs，高劑量丹參提取物大腸桿菌感染組本身含有698個OTUs，對照組與大腸桿菌組共有的OTUs數量為130個。丹參低劑量大腸桿菌感染組與大腸桿菌感染組共有的OTUs數量為156個，丹參中劑量大腸桿菌感染組與大腸桿菌感染組共有的OTUs數量為135個，丹參高劑量大腸桿菌感染組與大腸桿菌感染組共有的OTUs數量為136個；丹參低劑量大腸桿菌感染組與對照組共有的OTUs數量為409個，丹參中劑量大腸桿菌感染組與對照組共有的OTUs數量為429個，丹參高劑量大腸桿菌感染組與對照組共有的OTUs數量為403個。說明單獨感染大腸桿菌減弱了與正常組小鼠菌群的相似度，飼喂低、中、高劑量丹參的大腸桿菌感染組與對照組共有的OTUs數量分別增加了253、294和267個，說明丹參可以提高各組小鼠腸道內菌群的相似度，緩解大腸桿菌對小鼠腸道菌群的損害。

CON. 對照組；E. coli. 大腸桿菌感染組；EL. 低劑量丹參大腸桿菌感染組；EM. 中劑量丹參大腸桿菌感染組；EH. 高劑量丹參大腸桿菌感染組。下同CON. Blank control group; E. coli. E. coli infection group; EL. Low dose Radix Salvia miltiorrhiza infected with E. coli group; EM. Medium dose Radix Salvia miltiorrhiza infected with E. coli group; EH. High dose Radix Salvia miltiorrhiza infected with E. coli group. The same as below圖6 五組小鼠糞便的韋恩圖分析Fig.6 Venn diagram analysis of feces from five groups of mice

2.8 Alpha多樣性分析

Alpha多樣性指數中Chao指數和ACE指數反映樣品中群落的豐富度，Shannon指數和Simpson指數反映樣品中群落的多樣性，其中的Chao指數、ACE指數和Shannon指數越大，Simpson指數越小，說明樣品中的物種越豐富。本研究發現，與對照組相比，大腸桿菌感染組極顯著降低了腸道菌群豐富度(<0.01)，而丹參低、中、高劑量大腸桿菌感染組極顯著提高了腸道菌群豐富度(<0.01)，其中丹參高劑量組與其他兩個劑量組相比，效果更佳(表2)。與對照組相比，大腸桿菌組腸道菌群多樣性顯著降低(<0.05)，而丹參低、中、高劑量大腸桿菌感染組與對照組和大腸桿菌感染組相比，則顯著增加了腸道菌群多樣性(<0.05)。

表2 Alpha多樣性結果Table 2 The results of alpha diversity %

2.9 菌群相對豐度變化

2.9.1 門分類學水平上的差異 如表3與圖7所示，在門水平上相對排名前三的菌群依次為Bacteroidetes(擬桿菌門)、Firmicutes(厚壁菌門)、Proteobacteria(變形菌門)。與對照組相比，大腸桿菌感染組中的Bacteroidetes極顯著減少(<0.01)，Firmicutes豐度顯著減少(<0.05)，Proteobacteria豐度極顯著增高(<0.01)；與大腸桿菌感染組相比，添加丹參低，中，高劑量大腸桿菌組中的Bacteroidetes豐度極顯著增高(<0.01)、Firmicutes豐度顯著增高(<0.05)，Proteobacteria豐度極顯著降低(<0.01)。說明低、中、高劑量丹參提取物可以緩解大腸桿菌對腸道菌群的破壞，從而維護腸道菌群平衡。

表3 門水平排名前三菌群相對豐度的變化Table 3 Changes in relative abundance of top three microflora at phylum level %

不同顏色的柱子代表不同的物種，柱子的長短代表該物種所占比例的大小。下同The columns of different colors represent different species, and the length of the column represents the proportion of the species. The same as below圖7 門水平的樣本物種組成分析Fig.7 Analysis of sample species composition at phylum level

2.9.2 屬分類學水平上的差異分析 如表4與圖8所示，在屬水平上相對排名前五的菌群分別為(擬桿菌屬)、(螺桿菌屬)、(乳桿菌屬)、(普氏菌屬)和Unclassified(未知菌屬)。與對照組相比，大腸桿菌感染組中的、、、豐度極顯著降低(<0.01)，Unclassified豐度無顯著差異(>0.05)；與大腸桿菌組相比，添加低劑量丹參大腸桿菌組中的、、Unclassified豐度極顯著增加(<0.01)，豐度顯著增加(<0.05)，豐度無顯著差異(>0.05)；添加中劑量丹參大腸桿菌組中的、、豐度極顯著增加(<0.01)，、Unclassified豐度無顯著差異(>0.05)；添加高劑量丹參大腸桿菌組中的、、Unclassified豐度極顯著增加(<0.01)，、豐度顯著增加(<0.05)。

表4 屬水平排名前五菌群相對豐度的變化Table 4 Changes in relative abundance of top five microflora at genus level %

圖8 屬水平的樣本物種組成分析Fig.8 Analysis of sample species composition at genus level

2.9.3 差異菌群篩選結果 LEfSe分析展示了組間在豐度上有顯著差異的物種。從圖9可以看出，共發現了14個差異菌群(LDA>4，<0.05)。其中，空白對照組占1種，屬水平上有___；大腸桿菌感染組占6種，科水平上有Enterobacteriaceae，屬水平上有-、___1，__，，；丹參高劑量大腸桿菌感染組占7種，綱水平上有unclassified_p__Firmicutes，目水平上有unclassified_p__Firmicutes、Lachnospirales，科水平上有unclassified_p__Firmicutes、Lachnospiraceae，屬水平上有___、。

圖中不同顏色節點表示各級分類水平上的微生物群落，淡黃色節點表示在兩個時間點對比均無顯著差異的菌群；紅色表示兩個時間點對比差異顯著，且在術前中顯著富集的菌群；藍色表示兩個時間點對比差異顯著The nodes with different colors indicated the microbial communities at all levels of classification, and the nodes with pale yellow color indicated the microflora with no significant difference compared at two time points. The red color indicated the flora with significant difference in comparison between the two time points and significant enrichment in pre-surgery; Blue indicates significant differences between the two time points圖9 LEfSe 分析比較對照組、大腸桿菌感染組、丹參高劑量大腸桿菌感染組差異菌群Fig.9 Differential gut bacterial in control group, E. coli infection group and high-dose Radix Salvia miltiorrhiza infected with E. coli group determined by LEfSe analysis

2.10 功能預測分析

2.10.1 COG 數據庫比對結果 基于美吉生物一站式科研服務平臺PICRUSt分析平臺的功能預測分析將獲得的 OTU 豐度表與Greengene 數據庫比對，獲得相應的 COG 及 KO功能信息及其豐度。分析發現，EL、EM、EH組小鼠腸道菌群功能主要集中于氨基酸、碳水化合物運輸與代謝，未知功能，翻譯、核糖體結構和生物轉化，細胞壁/膜/生物合成等方面，此外還有與能量產生和轉換，復制、重組和修復，轉錄，無機離子運輸和代謝等相關的功能基因，隨著丹參提取物劑量的增加能夠增強這些功能趨勢(圖10)。

圖10 COG功能分類Fig.10 COG functional classification

2.10.2 KEGG 功能預測 根據KEGG數據庫中得到的代謝通路Level 2豐度分析結果，各組共有 44 種代謝通路，根據Level 1可以分為6類，其中新陳代謝12種、遺傳信息處理4、環境信息處理3種、細胞生理過程4種、人類疾病12種和組織系統9種，去除人類疾病作旭日圖。由圖11可知，一級代謝通路分別為新陳代謝37.50%、環境信息處理9.38%、遺傳信息處理 12.50%、細胞過程 12.50%、組織系統28.13%。二級主要代謝通路有Carbohydrate metabolism-碳水化合物代謝(3.99%)、Amino acid metabolism-氨基酸代謝(3.72%)、Signal transduction-信號轉導(7.40%)、Replication and repair-復制和修復(4.41%)、Cell growth and death-細胞生長和死亡(5.56%)、Endocrine system-內分泌系統(9.14%)、Immune system-免疫系統(5.63%)。

圖11 KEGG 一級和二級代謝通路Fig.11 KEGG of category level 1 and 2

2.10.3 腸道菌群基因功能預測 通過美吉PICRUt 分析發現，在KEGG L2水平，組與EH組在代謝途徑上存在明顯差異(圖12)。在EH組存在顯著差異的功能通路有11個(<0.01)，主要包括肽聚糖生物合成、錯配修復、DNA復制、蛋白質輸出、同源重組、核糖體、氨酰tRNA生物合成、光合作用、萜類骨架生物合成、賴氨酸生物合成、谷氨酰胺與谷氨酸代謝，在組有2個差異極顯著的功能通路(<0.01)，主要包括生物膜形成-銅綠假單胞菌、C5-支鏈二元酸代謝。

圖12 E. coli組與EH組在KEGG L2水平代謝途徑差異分析圖Fig.12 Analysis chart of metabolic pathway difference between E. coli group and EH group at KEGG L2 level

3 討 論

腸道菌群數量多、種類龐大而且功能復雜多樣，正常的菌群微生態在動物的發育、代謝、免疫、感染的預防及治療等方面發揮了舉足輕重的作用，而這種平衡的關系一旦被打破，則會引起疾病的發生。在過去的生產中，主要使用抗生素來解決腸道菌群失調等問題，其作用機理是殺滅或抑制畜禽體內(特別是消化道內)的有害微生物，調節畜禽胃腸道微生態平衡，從而改善動物機體健康，促進動物生長。然而長期使用抗生素會導致耐藥菌產生進而降低藥效，使疾病難控，此外抗生素殘留在動物性食品中也將危害人類健康。在此環境下，中藥以其很少產生耐藥性、毒副作用小、且能增強機體免疫功能等特點，具有廣闊的應用前景。臨床上中藥的有效成分與腸道菌群相互作用，大多數藥物經腸內菌群吸收代謝后發揮藥效。因此，本研究利用網絡藥理學理論，通過OB、DL 和Caco-2參數來篩選丹參的有效成分，為其在調節腸道菌群平衡中的應用提供理論依據。

本研究發現，通過篩選丹參的36種有效成分的相關靶點，經過生物信息學分析，丹參成分作用的生物學過程主要富集在細胞內類固醇激素受體信號通路、信號轉導、凋亡過程的正調控、細胞增殖調節、血壓調節等；這與先前的研究結果相一致，丹參具有通經止痛、活血化瘀及涼血消癰等功效，而且還有抗氧化、抗菌消炎、保護心臟、抗肝損傷等作用。可用于臨床上肝炎、腫瘤、冠心病等治療。羅曉梅等研究發現，在膿毒癥大鼠模型中，丹參注射液能改善心肌功能，減少心肌細胞凋亡，調節凋亡相關蛋白的表達，上調抗氧化通路，下調炎癥通路，進而對膿毒癥大鼠心肌損傷起到治療作用。此外，本研究揭示了丹參有效成分的分子功能富集分析主要在RNA聚合酶II轉錄因子活性，配體激活的序列特異性DNA結合、血清素結合、類固醇激素受體活性等方面；細胞成分主要富集在神經元投射、突觸前膜的組成成分、質膜的整體部件、谷氨酸能突觸等方面。值得注意的是，Luo等研究表明，用丹參提取物治療卵清蛋白(OVA)誘導的過敏性哮喘小鼠，可減少氣道炎癥細胞浸潤、Th1/Th2細胞因子量和杯狀細胞增生，減輕OVA致敏哮喘小鼠的氣道炎癥。Liu等采用丹酚酸B (SAB)治療實驗性肺纖維化，發現SAB通過保護內皮細胞免受氧化應激損傷而在肺纖維化中發揮抗炎作用，減少氧化應激損傷主要通過抑制內皮通透性和通過MAPK和NF-κB信號通路表達促炎細胞因子介導。由此可見，丹參具有抗菌、抗炎和減少機體氧化損傷的作用。此外，本研究通過“成分-靶點-信號通路”網絡互作分析結果發現，丹參治療感染的關鍵成分為丹參“酮、醌”型結構的脂溶性成分。據報道，從總丹參酮中分離得到丹參酮 Ⅰ、丹參酮 ⅡA、丹參酮 ⅡB、隱丹參酮、羥基丹參酮、丹參酸甲酯；另外還有4個新成分：二氫丹參酮 Ⅰ、丹參新鯤甲、乙、丙，所分離到的化學成分經體外抑菌試驗表明，具有二氫呋喃結構的隱丹參酮和二氫丹參酮 Ⅰ有較強的抑菌作用，而其他具有呋喃環結構的成分抑菌作用弱或沒有表現出抑菌作用。陳姍姍等研究發現，丹參新酮能抑制 THP-1 細胞的增殖并誘導其凋亡，其機制可能與下調2、1 基因的表達，Bax/Bcl-2 的值增大，激活促凋亡因子caspase-3有關。Wang等研究發現，丹參新酮以劑量依賴的方式下調-、-1b、-6、-8基因的表達，表明丹參新酮對炎癥性腸病的抗炎作用可能是通過調TLR4/NF-κB/IQGAP2信號通路實現的。由此可推測，丹參可能通過抑制大腸桿菌所引起的細胞凋亡以及機體的炎癥反應而發揮主要作用。本研究的丹參提取物雖然表現出了良好的改善腸道菌群結構和抗大腸桿菌感染的作用，但具體發揮作用的有效成分作用機制還有待進一步的研究。

此外，丹參抗感染的關鍵靶點KEGG 富集分析表明，cAMP信號通路和PI3K-Akt信號通路具有重要的作用。據報道，環磷酸腺苷(cyclic adenosine monophosphate, cAMP)是作為哺乳動物體內重要且保守的第二信使之一，可通過轉導細胞外信號參與調節多種器官和組織的發育及生理功能。劉陽研究發現，和胃理氣方通過激活 cAMP/PKA 信號通路，提高 cAMP和PKA的表達，能夠改善肝郁氣滯型 FD 模型大鼠的胃腸排空，調節胃腸功能，進而促進胃腸動力。值得注意的是，PI3K/Akt 信號通路由磷脂酰肌醇3-激酶(PI3K)和其下游分子絲氨酸/蘇氨酸蛋白激B(PKB/Akt)兩部分組成，參與多種生長因子、細胞因子和細胞外基質等的信號轉導，同時還參與細胞增殖、分化、凋亡和葡萄糖轉運等多種細胞功能的調節。本研究采用網絡藥理學技術得到丹參抗感染的主要通路為cAMP信號通路和PI3K-Akt信號通路，在前人的研究基礎上，得知這兩條通路可以作用胃腸和調節細胞功能，從而為丹參治療大腸桿菌感染所致的菌群失調提供了相關理論依據。

本研究OTU聚類分析和Alpha多樣性結果表明，感染大腸桿菌的小鼠腸道內菌群的豐富度與均勻度都呈下降趨勢，而灌服了低、中、高劑量丹參素提取物之后再感染大腸桿菌的小鼠糞便中，菌群豐富度與均勻度與正常組相似，說明灌服丹參提取物可以提高大腸桿菌感染小鼠腸道中菌群的多樣性，其中高劑量丹參提取物組效果最好。此外，通過菌群相似度分析可知，單獨感染大腸桿菌減弱了與正常組小鼠菌群的相似度，但飼喂低、中、高劑量丹參組與對照組共有的OTUs數量皆有所增加，說明丹參可以提高小鼠腸道內菌群OTUs數量，緩解大腸桿菌對小鼠腸道菌群的損害性。運用分類學組成，通過不同組分的菌群構成比例繪制柱狀圖發現，丹參主要影響門水平中的Bacteroidetes(擬桿菌門)、Firmicutes(厚壁菌門)、Proteobacteria(變形菌門)，前兩者是健康宿主共有的腸道微生物優勢菌群。擬桿菌門和厚壁菌門在大腸桿菌感染組中含量較低，在丹參低、中劑量大腸桿菌感染組依次升高，丹參高劑量大腸桿菌感染組有所降低，但整體為增長趨勢，因此丹參可以緩解大腸桿菌對腸道菌群的破壞，使腸道內菌群接近對照組，維護腸道菌群平衡，丹參中劑量組藥效最佳。在屬水平上，丹參主要影響的是(擬桿菌屬)、(螺桿菌屬)，擬桿菌屬是擬桿菌門中重要的糖化菌，可以產生豐富的糖苷水解酶，將三萜糖苷等轉化為極性較小且親脂性更高的小分子，從而更好地被腸道吸收，說明丹參誘導腸道菌群中擬桿菌屬的相對豐度增加有利于小鼠機體更好地維持生理狀態。此外，通過COG、KEGG代謝通路對比分析發現，丹參對抗大腸桿菌感染的功能主要集中在氨基酸運輸和代謝、碳水化合物轉運及代謝上，說明可能含有豐富的蛋白分解、轉運及代謝酶相關基因，但具體的作用機制還有待進一步的驗證。

4 結 論

本研究基于網絡藥理學發現，丹參治療感染作用機制可能是通過18種關鍵藥物活性成分作用于動物機體關鍵靶點PTGS2、SCN5A、ADRB2、OPRM1、CHRM3、NCOA1、PTGS1、ACHE、CA2、RXRA、OPRD1、ADRA1B、CHRM2、NCOA2、SERPIND1、ESR1、AR、DPP4，影響cAMP信號通路和PI3K-Akt信號通路等。丹參主要通過影響腸道菌群的Bacteroidetes(擬桿菌門)、Firmicutes(厚壁菌門)的豐度來緩解大腸桿菌感染對小鼠腸道菌群的損害作用，調控菌群平衡。