腸道菌群與結直腸癌進展關系的臨床研究探討*

張默涵 葛振宇 譚楊 谷金華 史立宏 張小茜

近年來，隨著中國人群飲食結構與生活方式逐漸西方化，相關癌癥（結直腸癌、前列腺癌、膀胱癌）發(fā)病率呈上升趨勢。結直腸癌近年來已成為全球，尤其是中國的常見惡性腫瘤[1]。大多數(shù)散發(fā)的結直腸癌遵循常規(guī)的“正常-腺瘤-癌”序列，與在特定進展階段發(fā)生的特定突變相關[2]。在分子水平上，有研究發(fā)現(xiàn)，微環(huán)境中的促炎細胞因子白細胞介素-6（interleukin 6，IL-6）可通過激活STAT3信號通路促進上皮間質轉化和侵襲促進結直腸癌的進展[3-4]。

大量的證據(jù)表明結直腸癌發(fā)生與腸道菌群、遺傳和其他因素密切相關[5-6]。有研究顯示，多種腸道菌群與結直腸癌進展有關，但特異菌仍未明確。有研究表明，在結直腸癌患者的糞便中Prevotellaceae、Fusobacterium和Peptostreptococcus等菌群顯著升高[7]。

本研究旨在尋找與結直腸癌進展相關的特定腸道菌群，通過16srRNA測序技術對健康對照組、腺瘤性息肉組和結直腸癌組的腸道菌群進行比較，并檢測3組樣本腸道組織的IL-6的表達，以評估腸道菌群對結直腸癌進展的影響，為進一步臨床治療策略的確定提供新思路。

1 材料與方法

1.1 臨床資料

回顧性選取2020年1月至2021年5月于濰坊醫(yī)學院附屬醫(yī)院就診的40例患者，根據(jù)疾病類型分為C組（20例結直腸癌患者）、A組（20例腺瘤性息肉患者）。另外，篩選于濰坊醫(yī)學院附屬醫(yī)院進行常規(guī)胃腸鏡健康查體的就診者20例，設為N組作為對照。納入標準：1）所有結直腸癌與腺瘤性息肉的診斷均由至少2名專業(yè)醫(yī)師依據(jù)結腸鏡表現(xiàn)與病理結果進行診斷。2）年齡在18～70歲，過去近6個月無抗生素與益生菌服用史。3）臨床資料完整。排除標準：1）具有其他胃腸道疾病或具有糖尿病、高血壓、肥胖癥、結石及各種存在肝、臟嚴重原發(fā)性疾病患者、精神疾病障礙不能耐受腸鏡檢查者；2）標本收集前3個月應用抗生素或益生菌的；3）依從性差，難以隨訪的患者。4）具有3年及以上酗酒史的。所有參加研究的受試者都被告知研究性質，并簽署知情同意書。

1.2 方法

1.2.1 糞便樣本采集和DNA提取 所有糞便樣本均采用無菌糞便盒收集，收集后分類編號并于2 h內置于-80℃冰箱儲存，以準備進一步提取DNA。使用Power Soil DNA提取試劑盒（MO Bio Laboratories公司，購自北京）按照說明書的操作規(guī)程從樣品中提取細菌總DNA。以260 nm/280 nm和260 nm/230 nm的比值評價DNA的質量和數(shù)量。然后將DNA保存在-80℃，等待下一步處理。

用通用引物對（正向引物：5’-ACTCCTACGGG AGGCAGCA-3’；反向引物：5‘-GGACTACHVGGG TWTCTAAT-3’）結合適配序列和條碼序列，擴增細菌16S rRNA基因的V3-V4區(qū)。擴增總體積為50 μL，包括10 μL Buffer，0.2 μL Q5高保真DNA聚合酶，10 μL高GC增強子，1 μL dNTP，每個引物10 μM，60 ng基因組DNA。熱循環(huán)條件為：在95℃下初始變性5 min，然后在95℃下變性1 min，在50℃下變性1 min，在72℃下加熱1 min，最后在72℃下延伸7 min。第1輪的PCR產物通過VahtsTM DNA清潔珠純化。第2輪PCR在40 μL反應中進行，該反應包含20 μL 2×Phμsion HF MM、8 μL ddH2O、每種引物10 μM和第1輪的PCR產物10 μL。熱循環(huán)條件如下：首先在98℃下變性30 s，然后在98°C下變性10 s，65℃變性30 s，72℃變性30 s，在72℃下延長5 min。最后用Quant-itTMdsdna-hs試劑對所有的PCR產物進行定量，并匯集在一起。利用中國北京Biomarker Technologies Corporation的Illumina HiSeq 2 500平臺（2×250對末端）對純化的混合樣品進行細菌rRNA基因的高通量測序分析。物種鑒定和分類細菌 16S：使用Silva（Release128，http：//www.arb-silva.de）數(shù)據(jù)庫對物種鑒定和分析。

1.2.2 免疫組織化學法檢測IL-6在3組腸道組織的表達 C組取材于結直腸癌患者術后的腫瘤組織，A組取材于腸鏡下所取的腺瘤性息肉患者的息肉組織，N組取材于健康體檢者的腸壁組織。將收集的C組、A組、N組的組織進行固定、包埋，利用免疫組織化學法分別檢測3組組織樣本IL-6的表達情況，IL-6定位于細胞質，以棕黃色顆粒為染色陽性。染色結果判定采用半定量結果判斷：每例切片隨機選取5個高倍視野，按染色強度及陽性細胞數(shù)所占百分比綜合計分。染色強度：細胞無染色為0分，淺黃色為1分，黃色為2分，棕褐色為3分；陽性細胞百分比：陽性細胞＜5%為0分， 5%～25%為1分，26%～50%為2分，＞50%為3分。將陽性細胞百分比和陽性強度兩項相乘，得到該視野的最終得分。每張切片的得分為該切片選取的5個視野得分的平均值。最后該切片按照得分為下述4個等級：0分為陰性（-），1～4分為弱陽性（+），5～8分為陽性（++），8分以上為強陽性（+++）。統(tǒng)計學上把強陽性（+++）與陽性（++）所占比率稱為為陽性率。

1.2.3 觀察指標 1）α多樣性分析：α多樣性反映的是單個樣品物種豐度及物種多樣性，其中度量菌群豐度的指標：Chao1指數(shù)；度量菌群多樣性的指標：香農指數(shù)，香農指數(shù)越大說明群落多樣性越高。

2）β多樣性分析：β多樣性分析用于比較不同樣品在物種多樣性方面存在的相似程度。基于binary_jaccard等算法呈現(xiàn)物種多樣性的矩陣；基于R語言平臺繪制主坐標分析（PCoA），通過主坐標分析實現(xiàn)多個樣品的分類，進一步展示樣品間物種多樣性差異；UPGMA分析，樣品越靠近，枝長越短，說明兩個樣品的物種組成越相似；Anosim分析，可以對不同分組的樣品之間beta多樣性是否有顯著性差異進行檢驗。

3）LefSe和Metastats分析：LefSe分析，即組間差異顯著物種分析（或biomarkers分析），采用線性判別分析（LDA）來估算每個組分（物種）豐度對差異效果影響的大小，該分析目的為找到組間在豐度上有顯著性差異的物種。Metastats分析，對組間的物種豐度數(shù)據(jù)進行t檢驗得到P值，并對P值進行校正得到Q值；最后根據(jù)P值（或Q值）篩選出導致兩組樣品組成差異的物種，P＜0.05，Q＜0.05為差異具有統(tǒng)計學意義。

1.3 統(tǒng)計學分析

采用R軟件用于所有統(tǒng)計分析和圖表構建。使用t檢驗來評估微生物分類群、臨床參數(shù)和多樣性指數(shù)差異的顯著性，計數(shù)資料采用百分率表示，多樣本率的比較采用KrusKal-Wallis秩和檢驗法（H檢驗法）。以P＜0.05為差異具有統(tǒng)計學意義。

2 結果

2.1 不同樣本間菌群的豐度差異

隨著測序條數(shù)的加大，檢測到的OTU數(shù)目逐漸增多，稀釋性曲線先急劇上升后趨于平緩（圖1A）。熱圖展示了不同樣本在屬的水平所包含優(yōu)勢菌群的豐度差異，呈現(xiàn)組內聚集現(xiàn)象（圖1B）。物種分布柱狀圖展示了不同樣本在屬的水平特征菌群所占的比例（圖1C）。綜上所述，在“健康—腺瘤性息肉—結直腸癌”3組中，每個樣本的菌群的豐度并不相同。

圖1 不同樣本間菌群的豐度差異

2.2 不同樣本間菌群的多樣性差異

反映物種多樣性的香農指數(shù)曲線隨著測序條數(shù)的加大先急劇上升后趨于平緩（圖2A）。基于binary_jaccard算法的PCoA顯示，3組樣本的多樣性同組之間相互聚集，不同組之間相互分離。（圖2B）。基于binary_jaccard算法的Anosim分析箱形圖可見R=0.322，P=0.001（圖2C）。

圖2 不同樣本間菌群的多樣性差異

2.3 部分菌群與結直腸癌的進展相關

C組、A組與N組的Chao1指數(shù)具有顯著性差異，且C組與A組的差異較大（P＜0.01），而A組與N組的差異較小（P＜0.05）（圖3A）。物種分布柱狀圖顯示在N組中Bifidobacterium、Faecalibacterium豐度較高、A組中Escherichia-Shigella、Enterobacteriaceae豐度較高、C組中Lachnospiraceae豐度明顯降低，Escherichia-Shigella、Enterobacteriaceae豐 度 較 高（圖3B）。分支圖和LEfSe分析顯示差異最顯著的菌為結直腸癌組中Prevotellaceae；腺瘤性息肉組中的Enterobacteriaceae、Erysipelotrichaceae和正常組中的Ruminococcaceae、Bifidobacteriaceae（LDA得分＞4，P＜0.05，圖3C，3D）。

2.4 IL-6在3組不同組織樣本中的表達水平比較

利用免疫組織化學法檢測3組不同腸道組織的IL-6的表達，IL-6的陽性部位在細胞質，呈棕黃色染色。結果顯示IL-6蛋白在C組中的陽性率為95%（19/20），其中強陽性為60%（12/20），陽性為35%（7/20），弱陽性為5%（1/20），陰性為0（0/20），在A組中的陽性率為60%（12/20），其中強陽性為5%（1/20），陽性為55%（11/20），弱陽性為40%（8/20），陰性為0（0/20），在N組中的陽性率為0（0/20），弱陽性為30%（6/20），陰性為70%（14/20）。IL-6在C組中的陽性率高于A組高于N組，且3組之間差異具有統(tǒng)計學意義（x2=43.570，P＜0.01）（表1，圖4）。

圖4 免疫組織化學檢測各組織中IL-6的表達水平 （IHC×400）

表1 IL-6蛋白在3組腸道組織中的表達情況

3 討論

近年來研究發(fā)現(xiàn)，腸道微生物可以通過多種方式影響結直腸癌的進展，例如誘導產生炎癥因子[8]、抑制NK細胞對腫瘤細胞的殺傷[9]等，也有部分腸道微生物可通過代謝食物成分產生短鏈脂肪酸抑制結直腸癌的進展[10]。與健康人相比，結直腸癌患者的糞便微生物豐度以及優(yōu)勢菌群有很大差異[11]。據(jù)報道，腸道菌群的差異可能會影響結直腸癌的發(fā)生和發(fā)展[12]。目前，16srRNA測序作為一種分析腸道微生物多樣性的工具已在全球范圍內廣泛應用[13]。盡管腸道菌群與腸道疾病之間的關系已有大量研究，但與結直腸癌進展相關的特異性菌群仍不明確。

本研究中采用16srRNA測序技術，探究在“健康-腺瘤性息肉-結直腸癌”的進展過程中受試者腸道菌群的數(shù)量和豐度差異。有研究表明與健康對照組相比，結直腸癌組的糞便腸道菌群多樣性增加[14-15]，本研究也證實了這一點。但也有其他研究表明，與健康對照組相比，結直腸癌組的腸道微生物豐度和多樣性降低[2,14]，這可能與地域差異、樣本數(shù)量或疾病的進展時期有關。通過α多樣性及β多樣性分析，發(fā)現(xiàn)結直腸癌組腸道微生物群落的組成及多樣性與腺瘤性息肉組及健康對照組之間有顯著性差異：即3組中，每個樣本的菌群豐度與多樣性并不相同，且呈現(xiàn)出組內樣本差異小，組間樣本差異大的現(xiàn)象。有研究表明，Prevotellaceae可能與結直腸癌有密切關系，其在結直腸癌患者的腸腔內明顯富集[16]。Prevotellaceae在既往研究中也被證明在AOM/DSS小鼠中的表達高于對照組，并可以增加DSS誘導的小鼠結腸炎的嚴重程度[17]。Bifidobacteriaceae豐度較高的一組可檢測到結腸中較高的短鏈脂肪酸水平[18]，另外有研究表明，Ruminococcaceae也產生短鏈脂肪酸[19]。在本研究中，初步探討了在“健康-腺瘤性息肉-結直腸癌”的進展過程中受試者腸道菌群的變化，并通過LEfSe分析進一步明確了結直腸癌組中差異最顯著的菌為Prevotellaceae，這與既往研究一致，說明該菌群可能參與了結直腸癌的進展；而腺瘤性息肉組中差異最顯著的菌為Enterobacteriaceae、Erysipelotrichaceae，提示上述菌群可能參與腺瘤性息肉的發(fā)生。而正常組中，差異最顯著的菌為Ruminococcaceae、Bifidobacteriaceae，說明這兩種菌可能抑制結直腸癌的進展。簡而言之，LEfSe分析揭示3組間差異最顯著的特定菌群。最后，本研究還驗證了大腸癌進展過程中健康受試者的腸壁組織、腺瘤性息肉患者的息肉組織和結直腸癌患者的腫瘤組織中IL-6的表達，IL-6在腫瘤組織中的陽性率高于腺瘤性息肉組高于健康組，且3組之間差異具有統(tǒng)計學意義（P＜0.01），說明隨著結直腸癌的進展，IL-6的表達逐漸增強，即慢性炎癥可能通過產生IL-6等炎性因子導致癌前病變，進而進展為結直腸癌。

綜上所述，本研究結果表明，腸道菌群在結直腸癌進展的不同階段可能作為評估患者病情的生物標志物，并可能通過IL-6影響癌癥的進展。鑒別癌前病變或腫瘤病變特有的細菌群落可提供更有效的診斷策略，后續(xù)將分析長期處于不同結直腸癌進展風險中的患者，并繼續(xù)加大樣本數(shù)量、深入探索更多通路，探究潛在的促進結直腸癌進展的生物標志物及其促癌作用的機制。