2022年8月19日北京中醫(yī)藥大學與電子科技大學合作，在國際期刊JCR Q1區(qū)《Frontiers in Pharmacology》(IF:5.988)在線發(fā)表了題為“Revealing the Mechanism of Huazhi Rougan Granule in the Treatment of Nonalcoholic Fatty Liver Through Intestinal Flora Based on 16S rRNA, Metagenomic Sequencing and Network Pharmacology”研究成果。該研究通過宏基因組測序，16S rRNA基因測序、網(wǎng)絡藥理學和分子對接等手段，探討了HRG（化滯柔肝顆粒）通過腸道菌群治療非酒精性脂肪肝(NAFL)的作用機制。百邁客為該研究提供了16S rRNA基因測序和宏基因組測序和部分分析工作。

研究背景

非酒精性脂肪肝(NAFL)的發(fā)病率呈逐年上升趨勢，越來越多的證據(jù)表明腸道菌群在NAFL中起著致病作用。化脂柔肝顆粒是臨床上常用的治療NAFL的藥物。據(jù)報道，它可以降脂保肝，但還沒有研究證實該藥物的作用是否與腸道菌群有關。因此，研究者從腸道菌群的角度，研究其作用是否與腸道菌群調節(jié)有關，以進一步探討其治療NAFL的機制。

技術路線

主要結果

1.高脂飼料引起的小鼠體內脂質蓄積

如圖1所示，與BC組相比，模型組大鼠肝臟脂肪堆積增多，肝臟體積增大，顏色變白，邊緣變鈍。給藥后，各組大鼠肝臟紅化程度均有不同程度的改善，尤以TH組為佳。HE和油紅O染色結果顯示，模型組大鼠肝組織內可見較多的球狀變性肝細胞和較多的紅色脂肪，而給藥組較少。MC組病理程度重，TH組輕，TL組與PC組之間無顯著差異。因此，HRG可改善HFD誘導的肝臟脂肪堆積和體重增加，并且可能以劑量依賴的方式改善肝臟損傷。

圖1.HRG對小鼠肝臟形態(tài)、肝組織病理學的影響

2.血清分析

與空白組比較，模型組大鼠肝酶、血脂指標升高，其中TC、HDL-C、LDL-C、AKP、ALT、AST差異有統(tǒng)計學意義,而甘油三酯則無顯著差異。給藥后肝酶和血脂均低于模型組。丙氨酸氨基轉移酶(ALT)、谷草轉氨酶(AST)和高密度脂蛋白膽固醇(HDL-C)在MC組和給藥組之間有顯著差異。各給藥組間TC、ALT、HDLc差異無統(tǒng)計學意義。HRG組與PC組比較，AKP、AST、TG、LDL-C差異有統(tǒng)計學意義。

TM組和TH組降低AKP效果優(yōu)于TL組，TG組降低AKP效果優(yōu)于TM組。PC組降低AKP的效果優(yōu)于TL組，但差于TM組和TH組。PC組降低低密度脂蛋白膽固醇的效果優(yōu)于TM組和TH組，表明高脂飼料可導致大鼠血脂和肝酶代謝異常，高脂高脂顆粒可改善高脂大鼠血脂代謝和肝酶活性。與陽性藥物相比，HRG在改善高密度脂蛋白膽固醇、甘油三酯、堿性磷酸酶方面具有優(yōu)勢，且改善程度可能呈劑量依賴性。各組間器官指數(shù)無顯著差異，HRG對小鼠生化指標和器官指數(shù)的影響如圖2所示。

圖2.HRG對小鼠生化指標及器官指標的影響

3.微生物多樣性分析

對糞便樣本進行16s rRNA基因測序分析，在門、綱、目、科和屬的水平上對微生物群進行了鑒定和分析，最終共鑒定到17門、26綱、56目、86科、187屬。主要物種為厚壁菌門和擬桿菌門，其次是放線菌門、疣微菌門和變形菌門。此外，物種主要分布在uncultured_bacterium_f_Muribaculacea屬和lachnospiraceae NK4A136 group，其次是uncultured_bacterium_f_Lachnospirae、Lactobacillus和Akkermansia。

在門水平（圖 3A、B），MC 組的厚壁菌門、變形菌門和放線菌門的水平高于 BC 組，而擬桿菌門和疣微菌門的水平較低。與BC組相比，給藥組厚壁菌門和變形菌門水平降低，尤其是TH組；放線菌在 TH 和 TM 組中減少，在 TM 組中最顯著，但在 PC 和 TL 組中增加。擬桿菌在 TH 和 TM 組中顯著高于 MC 組，Verrucomicrobia 在 TH 和 TM 組中增加最顯著。與BC組相比，MC組在屬水平上乳酸桿菌和脫硫弧菌的比例更高。相比之下，Akkermansia 和 Ruminococcaceae UGG-014 所占比例較低。給藥后乳桿菌（TH組下降多）和脫硫弧菌（TH組下降多）減少，Akkermansia（TH組增加多）和瘤胃球菌UGG-014（PC組增加多）增加。給藥后乳桿菌和脫硫弧菌減少，阿克曼氏菌增多，尤其是TH組； Ruminococcaceae UGG-014 增加，特別是在 PC 組。因此，初步推測HRG可增加擬桿菌門、疣微菌門、未培養(yǎng)菌f Lachnospiceae和Akkermansia，減少厚壁菌門和變形菌門、乳酸桿菌和脫硫弧菌。HRG可通過改變NAFL小鼠腸道菌群結構改善NAFL。

圖3 I物種組成分析

腸道菌群聚類顯示BC組和MC組樣本微生物群落差異明顯，給藥組、BC組和MC組微生物群落差異顯著。此外，Bray-Curtis PCoA 證實了這些發(fā)現(xiàn)（圖 3J）Curtis PCoA 證實了這些發(fā)現(xiàn)（圖 3J）。 PC1（15.21%）能較好地區(qū)分BC和MC組、MC和TM、TH組微生物群落。這些數(shù)據(jù)進一步表明BC和MC組在微生物群落結構上存在較大差異，并且HRG施用組和對照組在微生物群落結構上存在較大差異。

進一步通過LEfSe分析尋找表明組間差異biomarker。結果表明，在門水平上差異物種為擬桿菌門、放線菌門、厚壁菌門和疣微菌門在是顯著差異物種。在屬水平上，差異物種為雙歧桿菌、擬桿菌、Alloprevotella、Lactobacillus、Lachnospiraceae_NK4A136_group、Dubosiella、Faecalibaculum和Akkermansiaceae。AFL小鼠和HRG組的優(yōu)勢菌群主要為厚壁菌門、變形菌門、放線菌門和擬桿菌門。在屬水平上，差異物種主要為Adlercreutzia和Ruminococcaceae_UCG-013、Sphingomonas等。表明HRG可能通過改變厚壁菌門、變形菌門、放線菌門和擬桿菌門的菌群來改善NAFL。

對腸道菌群豐度和血清指標參數(shù)進行 RDA 分析，結果表明，AKP、HDL-C、TC、TG與物種分布高度相關，且呈同向分布。在屬水平上，它們與Faecalibaculum、Desulfovibrio、Dubosiella和Lactobacillus呈正相關，與Bacteroides、Akkermansia和Alloprevotella呈負相關。在門水平上，AKP、HDL-C和TG與物種分布高度相關，與厚壁菌門和放線菌門呈正相關，與Epsilonbacteraeota和Verrucomicrobia呈負相關。 Acidobacteria 和 Proteobacteria 與 AKP 和 HDLC 呈負相關，與 TG 呈正相關。擬桿菌與HDL-C呈正相關，與AKP、TG呈負相關。結果表明，這些菌落豐度的變化與NAFL小鼠給藥前后腸道菌群的結構變化密切相關。

圖3 beta多樣性分析，LEfSe分析和RDA分析

3.網(wǎng)絡藥理學

通過對藥材中重復成分的篩選和去除，得到249個有效成分。將上述化合物用于靶標預測、校正和重復靶標的刪除，共獲得1186個靶標的標準基因名稱。HRG的“草藥-化合物靶點”網(wǎng)絡如圖4A所示(黃色節(jié)點代表藥物靶點，綠色節(jié)點代表化合物，紅色節(jié)點代表草藥)，在20個節(jié)點中，高度值的14個化合物被選為關鍵化合物(表3)。NAFL的PPI網(wǎng)絡如圖4B所示，IFD的PPI網(wǎng)絡如圖4C所示。
將與NAFL相關的靶點、與IFD相關的靶點、與化合物對應的靶點同時合并，得到了HRG通過腸道菌群治療NAFL的靶點(圖4D)。對合并后的網(wǎng)絡進行模塊分析和細胞Hubba分析。模塊分析顯示模塊1(圖4E，Score=7.2)和模塊2(圖4F，Score=7.0)，而cell Hubba分析顯示前10個目標(圖4G)。最終獲得了10個潛在的治療靶點：CXCL10、CXCL8、ICAM1、IFNG、IL10、IL1B、IL2、IL4、IL6、TNF。圖4A-F中的節(jié)點大小與度值呈正相關。圖4B-F中的節(jié)點顏色與度值呈正相關。

圖4.網(wǎng)絡構建和關聯(lián)分析。

(A)HRG藥材靶點網(wǎng)絡。(B)與NAFL有關的PPI網(wǎng)絡。(C)與IFD有關的PPI網(wǎng)絡。(D)HRG-IFD-NAFL合并目標的PPI網(wǎng)絡。(E)模塊1(F)模塊2(G)Hub基因。(H)草藥-關鍵化合物-潛在靶向途徑網(wǎng)絡。(I)關鍵靶點可能導致的關鍵生物進展圖例。

利用10個關鍵治療靶點、14個關鍵化合物和相關通路，通過Cytoscape構建了“藥草-關鍵化合物-潛在靶點-通路”的網(wǎng)絡圖（圖4H）。排名前三的化合物被選為潛在化合物（槲皮素、木犀草素、山柰酚）。選擇前三個信號通路作為關鍵通路。繪制了 HRG 通過腸道菌群治療 NAFL 的通路機制（圖 4I）。將10個關鍵靶點與3個潛在化合物進行分子對接，共獲得18對對接結果（表4），使用 Pymol 可視化結果（圖 5）。

圖5. 分子對接模式圖

4.宏基因組測序

通過宏基因組測序研究基因功能，基于KEGG(圖6A、B)、GO(圖6C)、EGNOG(圖6D)、CARD(圖6E)和CAZY(圖6F)數(shù)據(jù)庫進行了功能注釋，共注釋到了6055個KO(KEGG Ontology)和2104個EC(Ease)。結果表明，相應的復制、推薦、修復和一般功能預測功能基因只占較高的比例。細胞壁/膜/被膜的生物發(fā)生和碳水化合物的運輸和代謝位居第二。CARD數(shù)據(jù)庫的標注結果表明，四環(huán)素類和多藥耐藥對應的耐藥基因相對含量較高，其次為氟喹諾酮類和氨基糖苷類。CAZy數(shù)據(jù)庫的注釋結果表明，糖苷水解酶(GH)、糖基轉移酶(GT)和非催化糖結合模塊(CBM)是前三位的碳水化合物酶。

圖6G-L顯示了在每個組和樣本的不同水平上的KEGG通路的組成。與空白組和模型組相比，模型組3個層級通路的豐度均有不同程度的增加。推測HRG可能通過增加與代謝調節(jié)相關的代謝途徑來改善NAFL。

進一步用MetagenomeSeq分析檢測第 3層級 KEGG的豐度差異(p<0.05)，豐度差異圖如圖6M-X所示。結果表明，MC組蛋白酶體豐度顯著高于BC組，MC組倍半萜和三萜生物合成顯著低于BC組。與MC組相比，PC組改變了這一趨勢。與BC組相比，HRG組降解苯甲酸氟化物和降解甲苯的豐度更高。與MC組相比，HRG組豐度較高的代謝途徑包括氨基苯甲酸降解、類胡蘿卜素生物合成等。因此，HRG可能通過代謝和細胞過程改善NAFL，代謝相關途徑在其中起重要作用。氨基苯甲酸的降解、類固醇的生物合成、倍半萜和三萜的生物合成是重要的代謝途徑。

圖5.宏基因組測序結果。 （Ⅰ：功能注釋分析。Ⅱ：KEGG通路組成和豐度直方圖。Ⅲ：metagenomeSeq差異通路豐度熱圖）。

總結

綜上所述，HRG可能通過改變微生物多樣性、結構和功能來改善高脂飲食誘導的NAFL。還有可能通過調節(jié)腸道相關的代謝途徑、炎癥反應和免疫反應來改善NAFL相關的脂肪堆積和肝臟損傷。這些結果有力地表明，HRG可能通過預防腸道菌群失調來減輕NAFL。此外，HRG可通過腸道菌群治療多種成分、多代謝途徑、多靶點的NAFL。本研究為NAFL的治療提供了新的思路，證實了腸道菌群與NAFL之間的關系，提示HRG具有良好的治療效果。

如果您對宏基因組測序技術感興趣，歡迎點擊下方按鈕聯(lián)系我們，我們將免費為您設計文章思路方案。

百邁客微生物組學服務

百邁客生物在微生物領域深耕多年， 合作文章見刊于《Nature》、《Cell Host & Microbe》、《Water Research》、《Journal of Hazardous Materials》、《Environmental science & Technology》、《Bioresource Technology》等高水平雜志。百邁客可提供Illumina二代宏基因組、ONT三代宏基因組、PacBio全長微生物多樣性、二代微生物多樣性、微生物絕對定量等全套測序分析服務，項目經(jīng)驗豐富，累計處理樣本數(shù)15萬+，合作文章數(shù)600+，影響因子2500+，合作單位覆蓋31個省及直轄市、400+家單位，百邁客致力于提供高質量的組學測序服務，歡迎垂詢！