亞健康狀態(tài)人群腸道菌群結(jié)構(gòu)明顯失調(diào)：基于某醫(yī)院150例護(hù)理人員分析

張星星，蔡文文，廖生武，何璇昱，楊秋玉，白 楊，阮偉清

南方醫(yī)科大學(xué)1南方醫(yī)院惠僑醫(yī)療中心，2護(hù)理學(xué)院，3南方醫(yī)院門急診片區(qū)，4南方醫(yī)院消化內(nèi)科，廣東 廣州510515

腸道菌群（GM）及其代謝產(chǎn)物在維持人體健康穩(wěn)態(tài)中起至關(guān)重要的作用［1,2］，近年來，亞健康狀態(tài)與腸道菌群失衡的關(guān)系越來越受到關(guān)注。亞健康狀態(tài)（SHS）是人體處于健康和疾病之間的一種狀態(tài)［3,4］，常伴有食欲不振、腹部不適、便溏便秘等胃腸道癥狀［5,6］，慢性疲勞綜合征（CFS）為其中一種亞型［7］。Sheedy等［8］基于細(xì)胞培養(yǎng)的技術(shù)，觀察到CFS人群腸道中的革蘭氏陽性需氧菌明顯占優(yōu)勢，尤其是腸球菌和鏈球菌。Fremont等［9］研究來自比利時(shí)和挪威的CFS組和健康對照組腸道菌群的差異，發(fā)現(xiàn)CFS組的類桿菌屬（Bacteroides）減少。Shukla等［10］收集了CFS組和對照組的運(yùn)動(dòng)前、運(yùn)動(dòng)后15 min、48 h、72 h的癥狀報(bào)告、血液和糞便標(biāo)本，發(fā)現(xiàn)CFS人群主要菌群的相對豐度增加更顯著，血液中細(xì)菌清除時(shí)間延長。

上述研究均提供了亞健康狀態(tài)人群腸道菌群失調(diào)的證據(jù)，但由于既往研究的納入標(biāo)準(zhǔn)、測序技術(shù)的不同，導(dǎo)致亞健康狀態(tài)人群具體何種腸道菌群改變及變化的特點(diǎn)尚無定論［11,12］。近年來，16S rRNA測序技術(shù)的發(fā)展，彌補(bǔ)了傳統(tǒng)鑒定方法無法準(zhǔn)確鑒定細(xì)菌種屬的缺點(diǎn)，已廣泛用于腸道菌群的研究［13］，但未進(jìn)一步明確何種腸道菌群失調(diào)在亞健康狀態(tài)發(fā)生發(fā)展中起著重要作用。在本研究中，我們運(yùn)用16S rRNA測序技術(shù)，比較亞健康人群和健康人群的菌群構(gòu)成差異，明確其差異物種。

1 資料和方法

1.1 臨床資料

1.1.1 一般資料 選取2020年11月～2021年5月南方醫(yī)科大學(xué)南方醫(yī)院工作的護(hù)理人員150例作為研究對象，其中亞健康人群75例為亞健康組（SHS組），男性6例，女性69例，年齡20～40歲；選擇75例同期健康人群作為健康對照組（H 組），男性2 例，女性73 例；年齡21～40歲。組間年齡、性別、BMI、吸煙、飲酒差異無統(tǒng)計(jì)學(xué)意義（P均＞0.05），SHS 組亞健康總分明顯高于H 組（P＜0.05）。由于抗生素和激素對腸道菌群影響較大［14］，要求所有志愿者在采樣前2月均未接受過任何抗生素和激素治療。所有參與志愿者均已簽署實(shí)驗(yàn)知情同意書。此項(xiàng)研究通過醫(yī)院倫理委員會(huì)批準(zhǔn)，所有納入的研究對象均簽署知情同意書（倫理編號(hào)：NFEC-2020-278）。

1.1.2 納入標(biāo)準(zhǔn) 亞健康組納入標(biāo)準(zhǔn)：（1）年齡在20～40歲，性別不限；（2）亞健康的判定標(biāo)準(zhǔn)參照《亞健康的中醫(yī)臨床研究指導(dǎo)原則（試行）》［3］，滿足以下條件：入組持續(xù)3個(gè)月以上反復(fù)出現(xiàn)的不適狀態(tài)或適應(yīng)能力顯著減退，但能維持正常工作；無重大器官器質(zhì)性疾病及精神心理疾??；盡管具有明確的非重大器官器質(zhì)性疾病或精神心理疾病診斷但無需用藥維持，且與目前不適狀態(tài)或適應(yīng)能力的減退無因果聯(lián)系；（3）亞健康評定量表診斷為亞健康狀態(tài)；（4）研究對象同意參加本次研究。健康組納入標(biāo)準(zhǔn)：（1）年齡20～40歲，性別不限；（2）無重大器官器質(zhì)性疾病及精神心理疾??；（3）盡管具有明確的非重大器官器質(zhì)性疾病或精神心理疾病診斷但無需用藥維持；（4）亞健康評定量表診斷為健康狀態(tài)；（5）研究對象同意參加本次研究。

1.1.3 排除標(biāo)準(zhǔn)（1）入組前2月內(nèi)曾使用抗生素、微生態(tài)制劑者；（2）過敏體質(zhì)或?qū)λ幬镞^敏者；（3）妊娠或準(zhǔn)備妊娠、哺乳期婦女；（4）懷疑或確有酒精、藥物濫用史者；（5）有消化系統(tǒng)疾病病史及手術(shù)史者；（6）有高血壓、糖尿病、冠心病、乙肝、肺結(jié)核病史者；（7）認(rèn)知功能障礙不能給予充分知情同意者。

1.2 調(diào)查工具

1.2.1 一般情況調(diào)查問卷 調(diào)查內(nèi)容包括年齡、性別、BMI、是否吸煙、飲酒。

1.2.2 亞健康評定量表［15］該量表包括3個(gè)子量表（生理亞健康、心理亞健康和社會(huì)亞健康）、9個(gè)維度（精力、身體運(yùn)動(dòng)功能、社會(huì)適應(yīng)、正向情緒、社會(huì)資源與社會(huì)支持、認(rèn)知功能、心理癥狀、器官功能、身體癥狀）和39個(gè)條目，其中正向計(jì)分20條，反向計(jì)分19條。采用Likert 5級評分，總分為0～100分，各維度、子量表及總量表得分越高，表示健康狀況越好。計(jì)算公式：轉(zhuǎn)化分=（原始分-理論最低分）/（理論最高分-理論最低分）×100。排除疾病人群后，當(dāng)生理子量表低于閾值68，或心理子量表低于閾值67，或社會(huì)子量表低于閾值67時(shí)，被診斷為亞健康狀態(tài)。該量表Cronbach's α系數(shù)為0.917，重測信度為0.644。

1.3 標(biāo)本收集

本研究采取病例對照研究的方法，根據(jù)納入、排除標(biāo)準(zhǔn)招募調(diào)查對象，首先由2名課題組成員（經(jīng)過量表評定一致性培訓(xùn)）對調(diào)查對象進(jìn)行亞健康評定量表評定，并收集調(diào)查對象的一般資料，當(dāng)場回收問卷，后期由雙人錄入電腦中。所有研究對象在清晨（6:00-9:00）用糞便采集管（IGE）收集中間段沒有接觸空氣和地面部分新鮮的糞便標(biāo)本2～3 g，按操作說明混勻后在15分鐘內(nèi)將裝有樣品的采集管置于-80 ℃冰箱保存。

1.4 DNA提取和PCR擴(kuò)增

基因組DNA 通過E.Z.N.A.?DNA Kit（Omega Bio-tek,Norcross,GA,U.S.）試劑盒按照說明書進(jìn)行提取。使用紫外分光光度計(jì)（NanoDrop 2000,Thermo Scientific)檢測DNA的濃度和純度，用2%瓊脂糖凝膠電泳鑒定DNA完整性。使用引物515F（5'-GTGCCA GCMGCCGCGGTAA-3'）和806R6（5'-GGACTACⅤS GGGTATCT AAT-3'）擴(kuò)增細(xì)菌16S rRNA基因的Ⅴ3-Ⅴ4高可變區(qū)。使用qiagen公司提供的膠回收試劑盒回收產(chǎn)物，Tris-HCl 洗脫，2%瓊脂糖電泳檢測。

1.5 文庫構(gòu)建和上機(jī)測序

使用TruSeq?DNA PCR-Free Sample Preparation Kit建庫試劑盒進(jìn)行文庫構(gòu)建，構(gòu)建好的文庫經(jīng)過Qubit和Q-PCR定量，文庫合格后，使用NovaSeq6000進(jìn)行上機(jī)測序。

1.6 生物信息分析

樣本委托北京諾禾致源股技術(shù)有限公司進(jìn)行高通量測序，根據(jù)Barcode序列和PCR擴(kuò)增引物序列從下機(jī)數(shù)據(jù)中拆分出各樣本數(shù)據(jù)，截去Barcode和引物序列后使用FLASH軟件（Ⅴ1.2.7）對每個(gè)樣本的reads進(jìn)行拼接，得到原始Tags數(shù)據(jù)（Raw Tags）。拼接得到的Raw Tags，需要經(jīng)過嚴(yán)格的過濾處理得到高質(zhì)量的Tags數(shù)據(jù)（Clean Tags）。參照Qiime（Ⅴ1.9.1）的Tags質(zhì)量控制流程，進(jìn)行如下操作：a）Tags截?。簩aw Tags從連續(xù)低質(zhì)量值（默認(rèn)質(zhì)量閾值為≤19）堿基數(shù)達(dá)到設(shè)定長度（默認(rèn)長度值為3）的第一個(gè)低質(zhì)量堿基位點(diǎn)截?cái)?；b）Tags長度過濾：Tags經(jīng)過截取后得到的Tags數(shù)據(jù)集，進(jìn)一步過濾掉其中連續(xù)高質(zhì)量堿基長度小于Tags長度75%的Tags。經(jīng)過以上處理后得到的Tags需要進(jìn)行去除嵌合體序列的處理，Tags序列通過與物種注釋數(shù)據(jù)庫進(jìn)行比對檢測嵌合體序列，并最終去除其中的嵌合體序列，得到最終的有效數(shù)據(jù)。通過USEARCH軟件對序列進(jìn)行聚類，將序列相似度≥97%的序列歸為同一OTU，然后通過OTU與數(shù)據(jù)庫比對，對OTU進(jìn)行物種注釋。α多樣性分析用Chao1指數(shù)和Shannon指數(shù)來評估。β多樣性用基于Weighted Unifrac 距離的PCoA 分析?；贚EfSe方法分析，LDA最小值設(shè)為4，進(jìn)行2組之間的比較，尋找在2組之間具有顯著差異的菌屬。

1.7 統(tǒng)計(jì)學(xué)方法

采用SPSS 23.0統(tǒng)計(jì)學(xué)軟件進(jìn)行數(shù)據(jù)分析，計(jì)數(shù)資料采用頻數(shù)和構(gòu)成比描述，組間比較采用卡方檢驗(yàn)；正態(tài)分布的計(jì)量資料以均數(shù)±標(biāo)準(zhǔn)差表示，組間比較采用獨(dú)立樣本t檢驗(yàn)；偏態(tài)分布的計(jì)量資料，采用中位數(shù)和四分位數(shù)M（P25，P75）描述，組間比較采用Mann-WhitneyU檢驗(yàn)。檢驗(yàn)水準(zhǔn)α=0.05。

2 結(jié)果

2.1 亞健康組與健康組一般情況比較

2組年齡、性別、BMI、吸煙、飲酒差異無統(tǒng)計(jì)學(xué)意義（P均＞0.05），SHS組亞健康總分明顯高于H組（P＜0.05）（表1）。

表1 亞健康組與健康組一般情況比較Tab.1 Demographic and clinical characteristics of sub-health status group and control group

2.2 亞健康組與健康組腸道菌群微生物多樣性分析

2.2.1 測序數(shù)據(jù) 2組150個(gè)樣本共獲得平均每樣品測得96 391 條tags。以97%的一致性將序列聚類成為OTUs（Operational Taxonomic Units），共得到3616 個(gè)OTUs。SHS組獨(dú)有494 個(gè)OTU，H組獨(dú)有1784 個(gè)OTU，兩組共有1333個(gè)OTU（圖1）。

圖1 SHS組與H組Ⅴenn圖Fig.1 Venn diagram of OTU in the two groups.H:Healthy group;SHS:Sub-health status group.

2.2.2 α多樣性比較 等級聚類曲線（Rank Abundance曲線）結(jié)果表明2組腸道菌群的群落豐度較高，均勻度較好（圖2）。物種累積箱形結(jié)果顯示箱形圖位置趨于平緩（圖3），當(dāng)前樣本量足夠，現(xiàn)有的測序數(shù)據(jù)量可以反映樣本中菌群的多樣性信息。本研究中菌群豐度用Chao1指數(shù)來衡量；菌群多樣性用Shannon 指數(shù)來衡量，Shannon指數(shù)越大，樣品中的物種多樣性越高。結(jié)果顯示，SHS組的Chao1指數(shù)與H組之間差異無統(tǒng)計(jì)學(xué)意義（P=0.619，圖4），SHS組的Shannon指數(shù)升高且2組之間差異有統(tǒng)計(jì)學(xué)意義（P＜0.001，圖5）。

圖2 等級聚類曲線Fig.2 Rank abundance curves of the two groups.

圖3 物種累積箱型圖Fig.3 Species accumulation boxplots.

圖4 Chao1指數(shù)組間差異箱形圖Fig.4 Boxplot of Chao1 index in the two groups.

圖5 Shannon指數(shù)組間差異箱形圖Fig.5 Boxplot of Shannon index in the two groups.

2.2.3 β多樣性比較 主坐標(biāo)分析PCoA 圖顯示，除個(gè)別樣本存在偏差，亞健康組和健康組基本呈現(xiàn)在不同區(qū)域（圖6）。經(jīng)過Anosim相似性分析發(fā)現(xiàn)，2組間群落結(jié)構(gòu)呈現(xiàn)顯著區(qū)別（R=0.227，P=0.001）。

圖6 基于Unweighted Unifrac距離的PCoA分析Fig.6 PCoA analysis based on Unweighted Unifrac distance.

2.2.4 腸道菌群結(jié)構(gòu)分析 在門水平，SHS組中豐度最高的4 個(gè)門依次為：Firmicutes（59.5%）、Bacteroidota（31.9%）、Actinobacteriota（4.4%），變形菌門占2.5%。而HS組中，Bacteroidota（48.4%）、Firmicutes（41.1%）、變形菌門占5.3%、Actinobacteriota（2.5%）（圖7）。

圖7 門水平上2組的相對豐度Fig.7 Relative abundance of bacterial phyla in the two groups.

在屬水平上，H組中豐度最高的依次是：Prevotella（25.6%）、Bacteroides（19.5%）、Faecalibacterium（12.2%）、Blautia（3.3%）、Megamonas（2.1%）、Lachnospira（1.8%）、Bifidobacterium（1.8% ）、Alloprevotella（1.3% ）、Fusobacterium（1.3%）、Romboutsia（0.9%）。而在SHS組中豐度最高的依次是：Bacteroides（16.4%）、Prevotella（13.5%）、Blautia（11.8%）、Faecalibacterium（10.8%）、Bifidobacterium（3.6%）、Romboutsia（1.4%）、Lachnospira（1.3%）、Megamonas（0.6%）、Fusobacterium（0.5%）、Alloprevotella（0.3%）（圖8）。

圖8 屬水平上2組的相對豐度Fig.8 Relative abundance of bacterial genera in the two groups.

2.2.5 腸道菌群差異分析 在LDA值分布柱狀圖中可以看出，SHS組的優(yōu)勢物種為厚壁菌門，H組的優(yōu)勢物種為擬桿菌門。通過LEfSe分析得出兩組在4個(gè)門3個(gè)綱3個(gè)目3個(gè)科3個(gè)屬的相對豐度差異有統(tǒng)計(jì)學(xué)意義（P＜0.05，圖9）。

圖9 LDA分布柱狀圖Fig.9 Histogram of LDA score distribution.

3 討論

本研究運(yùn)用16S rRNA測序技術(shù)，對比亞健康組與健康組之間的腸道菌群，發(fā)現(xiàn)了亞健康組人群腸道菌群的變化特點(diǎn)，與健康組相比，亞健康組人群的腸道菌群組成結(jié)構(gòu)和多樣性上均存在顯著差異。結(jié)果顯示，在α多樣性分析的Shannon指數(shù)上，SHS組較H組明顯升高，且β多樣性分析也進(jìn)一步顯示2組間腸道菌群組成結(jié)構(gòu)上存在差異，表明亞健康人群腸道菌群多樣性增加，與健康人群相比結(jié)構(gòu)失衡。

本研究發(fā)現(xiàn)SHS組的Shannon指數(shù)較H組明顯升高，這一結(jié)果與Armstrong［16］的研究結(jié)果相反，后者通過CFS組和健康組的腸道菌群對比分析，發(fā)現(xiàn)CFS組的物種多樣性減少。研究結(jié)果差異的原因可能為：（1）與樣本來源的人種、地域等因素造成的差異有關(guān)［17］；（2）使用了不同的微生物鑒定技術(shù)［18］。

在門水平，2組的優(yōu)勢菌門均為擬桿菌門、厚壁菌門、放線菌門，在屬水平，2 組的優(yōu)勢菌屬均為Prevotella、Bacteroides。這與人類腸道微生物群研究的先前報(bào)告一致［19］，人類腸道主要由厚壁菌門或擬桿菌門組成，厚壁菌門的優(yōu)勢屬為Ruminococc、Blautia、Eubacterium、Faecalibacterium，擬桿菌門的優(yōu)勢屬為Bacteroides、Prevotella。

為了進(jìn)一步確定關(guān)鍵菌群，本研究采用LEfSe方法，結(jié)果發(fā)現(xiàn)，亞健康組的厚壁菌門、放線菌門顯著高于健康組，變形菌門、擬桿菌門顯著低于健康組。這與Giloteaux等［20］的部分研究結(jié)果相似。Giloteaux等對一人患有CFS和另一人健康的單卵雙胞胎，進(jìn)行為期兩天的心肺運(yùn)動(dòng)測試，并對運(yùn)動(dòng)前和運(yùn)動(dòng)后的糞便標(biāo)本進(jìn)行了16S rRNA測序分析。在人類腸道微生物中，厚壁菌門與擬桿菌門（F/B）比率升高被廣泛認(rèn)為是腸道失調(diào)的標(biāo)志，并與肥胖［21］、糖尿病［22］和心血管［23］等疾病有關(guān)。我們在本次研究中也得到亞健康組的F/B比率的比值高于健康組，這可能與亞健康狀態(tài)的發(fā)生發(fā)展有密切聯(lián)系。

本研究還發(fā)現(xiàn)，在屬水平，亞健康組的Blautia、Agathobacter顯著高于健康組，Prevotella顯著低于健康組，Giloteaux等［20］也發(fā)現(xiàn)CFS組的Prevotella顯著低于健康組。Prevotella分解蛋白質(zhì)和碳水化合物，發(fā)酵膳食纖維，產(chǎn)生短鏈脂肪酸，與葡萄糖穩(wěn)態(tài)和宿主代謝有關(guān)［24］，在自閉癥兒童中也發(fā)現(xiàn)Prevotella的減少［25］。值得注意的是在既往研究中并未發(fā)現(xiàn)Blautia的變化。Blautia屬于毛螺菌科，可導(dǎo)致人體血漿IL-1β水平升高，睡眠障礙已被證明與IL-1β升高有關(guān)［26］，Smith［27］發(fā)現(xiàn)Blautia的豐度與睡眠質(zhì)量呈負(fù)相關(guān)。李媛媛等［28］也提出，Blautia是區(qū)分慢性失眠患者與健康人的特征細(xì)菌。睡眠障礙是亞健康狀態(tài)發(fā)生發(fā)展的危險(xiǎn)因素，也是亞健康狀態(tài)的重要表現(xiàn)之一［29］，腸道菌群的改變與亞健康狀態(tài)人群的睡眠障礙及疲憊癥狀相關(guān)［30］，提示我們之后的研究可以從這一方面進(jìn)行深入研究。

綜上所述，我們發(fā)現(xiàn)亞健康狀態(tài)人群的腸道菌群具有整體結(jié)構(gòu)失衡的特征，厚壁菌門與擬桿菌門比率升高，并明確了差異物種，亞健康人群的Blautia、Agathobacter顯著增加，Prevotella顯著減少。目前，亞健康狀態(tài)的病因和發(fā)病機(jī)制至今尚未清楚，但菌群穩(wěn)態(tài)失調(diào)可能在亞健康狀態(tài)的發(fā)生發(fā)展中發(fā)揮重要的作用。本研究通過運(yùn)用高通量測序技術(shù)，將亞健康人群與健康人群進(jìn)行比較研究，不僅發(fā)現(xiàn)二者呈現(xiàn)明顯的不同，驗(yàn)證了亞健康狀態(tài)影響腸道菌群多樣性以及結(jié)構(gòu)的組成，而且還將差異物種進(jìn)一步鑒定至菌屬，有望為治療亞健康狀態(tài)提供思路。本研究的不足之處在于，僅初步分析亞健康人群腸道菌群的分布情況，隨著多組學(xué)技術(shù)的發(fā)展，下一步研究將結(jié)合代謝組學(xué)、蛋白質(zhì)組學(xué)等方法確定具有核心功能的亞健康人群腸道菌群種類。