量化人類腸道菌種對艱難梭菌長期生長的影響（一）
摘要Summary
艱難梭菌(Clostridioides difficile)能夠短暫或長期定殖于人類腸道，并可能引發(fā)感染風(fēng)險。這一過程受到人類腸道微生物群中復(fù)雜的分子與生態(tài)相互作用的深刻影響。通過研究艱難梭菌在人體腸道群落中經(jīng)歷數(shù)百代的動態(tài)變化，我們發(fā)現(xiàn)其定殖模式包括穩(wěn)定共存、不穩(wěn)定性或競爭性排除。降低碳水化合物濃度使包含艱難梭菌和常見腸道共生菌普通擬桿菌的群落從競爭性排除轉(zhuǎn)變?yōu)楣泊?，而這種轉(zhuǎn)變得益于交互喂養(yǎng)的增強(qiáng)。在這種環(huán)境下，艱難梭菌發(fā)生了兩處關(guān)鍵突變，使其代謝生態(tài)位從依賴脯氨酸轉(zhuǎn)向依賴葡萄糖。這些代謝改變顯著影響了艱難梭菌與腸道微生物群的種間相互作用，同時降低了小鼠的疾病的嚴(yán)重程度。綜上所述，與普通擬桿菌的相互作用在塑造艱難梭菌的長期生長動態(tài)和進(jìn)化適應(yīng)中發(fā)揮了關(guān)鍵作用，為制定針對艱難梭菌的防治策略提供了重要啟示。
結(jié)果Results
人類腸道菌種在群落中對艱難梭菌的長期生長具有差異性影響
艱難梭菌(C.difficile)可在不同生命階段的個體中定殖，持續(xù)時間各不相同，而人們對決定這些動態(tài)的因素仍然知之甚少。由于腸道微生物生態(tài)是艱難梭菌定殖的重要決定因素，我們采用了自下而上的微生物群落實(shí)驗(yàn)，研究種間相互作用如何塑造艱難梭菌的長期生長動態(tài)。為此，我們將艱難梭菌DSM 27147(流行性核型027的參考菌株R20291)與9種多樣的人類腸道菌種共同培養(yǎng)在支持其生長的特定培養(yǎng)基中(圖1A和S1A；表S1和S2)。這些人類腸道菌種在不同個體中高度普遍存在，并涵蓋了人類腸道微生物群的系統(tǒng)發(fā)育多樣性。該群落包括以下物種：閃爍梭菌(Clostridium scindens,CS)：此前研究表明其在無菌小鼠中可抑制艱難梭菌的生長；梭菌(Clostridium hiranonis,CH)：通過代謝生態(tài)位重疊抑制艱難梭菌，并導(dǎo)致其代謝發(fā)生顯著變化；擬桿菌屬物種(包括單形擬桿菌(Bacteroides uniformis,BU)、多形擬桿菌(Bacteroides thetaiotaomicron,BT)和普通擬桿菌(Phocaeicola vulgatus,PV)：在不同環(huán)境下可抑制或促進(jìn)艱難梭菌的生長。此外，脫硫弧菌(Desulfovibrio piger,DP)、遲緩埃格特菌(Eggerthella lenta,EL)、產(chǎn)氣柯林斯菌(Collinsella aerofaciens,CA)和普雷沃菌(Prevotella copri,PC)與艱難梭菌的相互作用也已在短時間尺度(約5–10代)內(nèi)得到廣泛研究。
圖1|艱難梭菌在人類腸道群落中的長期生長動態(tài)
(A)艱難梭菌在人類腸道群落中長期生長實(shí)驗(yàn)示意圖。根據(jù)每個物種的16S rRNA序列生成了系統(tǒng)發(fā)育樹。
(B)利用時間序列絕對豐度數(shù)據(jù)推斷gLV模型參數(shù)。
(C)利用35次傳代的豐度數(shù)據(jù)推斷出艱難梭菌與9種腸道菌種之間的長期種間相互作用網(wǎng)絡(luò)。節(jié)點(diǎn)大小代表物種的環(huán)境容納量，邊的寬度代表長期種間相互作用系數(shù)(aij)的大小。腸道菌種之間的完整aij見圖S3A。
(D)艱難梭菌與腸道細(xì)菌在成對(左)和三成員(右)群落中的代表性群落動態(tài)。由彩色線條連接的數(shù)據(jù)點(diǎn)表示每個物種的相對豐度(左y軸)，而帶有95%置信區(qū)間陰影的灰色線條表示群落OD600(右y軸)。完整的群落動態(tài)見圖S1C(成對)和圖S2A(三成員)。穩(wěn)定共存、不穩(wěn)定共存和競爭排斥的分類依據(jù)圖S2B。(E)使用最后一個時間點(diǎn)相對于初始時間點(diǎn)的香農(nóng)多樣性來測量所有包含x軸所示特定物種的成對和三成員群落中共存程度(ΔH)。
(F)在所有包含x軸所示特定物種的成對和三成員群落中，35次傳代過程中各時間點(diǎn)測量值之間歐幾里得距離的箱線圖。對于(E)和(F)，由于組間重疊群落的存在，某些數(shù)據(jù)點(diǎn)跨越了較大的多樣性空間，從而導(dǎo)致了變異性(示例數(shù)據(jù)點(diǎn)已突出顯示)。圖上方的星號表示該組與另一組之間雙向曼-惠特尼U檢驗(yàn)的統(tǒng)計學(xué)顯著p值，星號顏色與另一組對應(yīng)。*p < 0.05，**p < 0.01，***p < 0.001。確切的 p 值見表 S4。
我們組裝了所有可能的雙物種組合(共9種)和三物種組合(共36種)，其中包含艱難梭菌(C.difficile)和腸道細(xì)菌。這些菌群培養(yǎng)了24小時后，取少量樣品以40倍稀釋接種至新鮮培養(yǎng)基，連續(xù)傳代35次(約186代)(詳見STAR方法)。通過16S rRNA測序確定每種物種的相對豐度，并乘以各時間點(diǎn)的總生物量，估算每種物種的絕對豐度。此外，我們單獨(dú)培養(yǎng)了每種物種，采用相同的實(shí)驗(yàn)設(shè)計(圖S1B)。結(jié)果顯示，艱難梭菌在60%的群落中持續(xù)存在(相對豐度>1%，直到最后一次傳代)，包括7/9的雙物種群落和20/36的三物種群落(圖S1C和S2A)。然而，在84%的群落中，艱難梭菌的絕對豐度從初代到最后一次傳代顯著下降(圖S1D)。某些物種(如PC和BT)在明顯抑制艱難梭菌的群落中被顯著富集。所有擬桿菌屬(Bacteroides)物種在單一培養(yǎng)或與艱難梭菌共培養(yǎng)時均顯著降低了培養(yǎng)基的pH值(培養(yǎng)基pH從7.1降至約4.9–5.3；圖S1E)。雖然其他腸道菌在單一培養(yǎng)或共培養(yǎng)中未顯著降低pH值，但PC僅在與艱難梭菌共培養(yǎng)時顯著酸化了培養(yǎng)基(pH從7.1降至5.1)，這可能與PC對艱難梭菌的抑制效應(yīng)有關(guān)(圖S1D)。
為了量化人類腸道菌種對艱難梭菌長期生長的影響，我們采用廣義Lotka-Volterra(gLV)模型，對單一培養(yǎng)、雙物種和三物種群落在35次傳代中的物種豐度時間序列數(shù)據(jù)進(jìn)行了擬合(圖1B；表S3 DATASET001；詳見STAR方法)。gLV模型是一種動態(tài)生態(tài)模型，可根據(jù)各物種的生長和與群落成員的雙物種相互作用預(yù)測群落動態(tài)。該模型廣泛用于解析物種間相互作用并預(yù)測不同環(huán)境中的群落組裝過程。推導(dǎo)出的gLV物種間相互作用系數(shù)量化了特定物種對另一物種長期生長的影響(圖S3A)。模型與實(shí)驗(yàn)數(shù)據(jù)的擬合度較高(Pearson相關(guān)系數(shù)R=0.89–0.91，p=3.6E?42至2.3E?51；圖S3B)。為了評估模型對驗(yàn)證數(shù)據(jù)的預(yù)測性能，我們進(jìn)行了留一交叉驗(yàn)證(圖S3C；詳見STAR方法)。由于用于訓(xùn)練模型的數(shù)據(jù)有限，特別是腸道菌的相關(guān)數(shù)據(jù)，模型僅對艱難梭菌的豐度預(yù)測表現(xiàn)出中等一致性(艱難梭菌的Pearson相關(guān)系數(shù)R=0.36–0.60，總體相關(guān)系數(shù)R=0.63–0.70，p=1.5E?13至1.9E?18)。此外，gLV模型未能捕捉物種間相互作用的進(jìn)化變化，這可能是預(yù)測性能較低的原因之一。
相關(guān)新聞推薦
1、源自丹麥奶酪鹽水的天然漢森德巴利菌菌株生長曲線及抗真菌活性
2、呼寧散對雞肺源大腸桿菌生長曲線、細(xì)胞壁的影響及抑制效果——材料與方法
3、基于生長速率法評估植物乳桿菌Bx6-2的發(fā)酵液對甜瓜疫霉菌抑制效果（二）