鄭 健，楊鎖柱，周艷萍，向泓銓，趙 津*

(1.玉溪農(nóng)業(yè)職業(yè)技術(shù)學(xué)院，云南玉溪 653100；2.玉溪華大高原農(nóng)業(yè)基因測(cè)序中心，云南玉溪 653100)

螺螄是一種軟體動(dòng)物，隸屬田螺科(Viviparidae)螺螄屬(Margarya)，僅分布在我國(guó)云南省的高原湖泊中，為我國(guó)特有物種[1-2]。螺螄不僅是高原湖泊物種多樣性的重要組成部分，而且被列為《國(guó)家重點(diǎn)保護(hù)動(dòng)物名錄》二級(jí)保護(hù)動(dòng)物。除了在物種多樣性方面具有重要的生態(tài)地位外，螺螄在重金屬污染、水體有機(jī)污染等方面也有著獨(dú)特的生態(tài)指示功能[3-4]。滇池是云南省最大的高原淡水湖泊。王丑明等[5]2011年調(diào)查結(jié)果表明，滇池中的螺螄屬包括螺螄(Margaryamelanioides)、牟氏螺螄(M.monodi)、光肋螺螄(M.mansugi)、長(zhǎng)螺螄(M.elongata)、 乳頂螺螄(M.tropidophora)等。其中，牟氏螺螄(Margaryamonodi)是螺螄屬中的一個(gè)特殊物種，僅分布在滇池。20世紀(jì)80年代以來(lái)水體污染、過(guò)度捕撈越來(lái)越嚴(yán)重，包含牟氏螺螄在內(nèi)的螺螄種群數(shù)量急劇減少，2009年螺螄被世界自然保護(hù)聯(lián)盟列為極危物種。牟氏螺螄的種群變化不僅是螺螄屬物種多樣性變化的指征，而且是滇池水環(huán)境的重要指征。然而，目前對(duì)牟氏螺螄在內(nèi)的螺螄屬生物的相關(guān)研究還相對(duì)較少，特別是對(duì)螺螄生存的微生物環(huán)境等方面缺乏了解。

近年來(lái)，隨著測(cè)序技術(shù)的進(jìn)步，學(xué)者對(duì)人和其他動(dòng)物腸道微生物的研究證實(shí)了腸道微生物對(duì)宿主健康、行為有著巨大的影響[6-7]。對(duì)于水生生物及軟體動(dòng)物，腸道微生物同樣與它們的生存環(huán)境及健康密切相關(guān)[8-10]。腸道菌群作為螺螄體內(nèi)微生物生態(tài)的一部分，其群落構(gòu)建的機(jī)制及其對(duì)螺螄的影響尚不清楚。筆者通過(guò)使用3代全長(zhǎng)擴(kuò)增子測(cè)序技術(shù)對(duì)牟氏螺螄腸道菌群的16S rRNA基因進(jìn)行高通量測(cè)序，比較野生和人工飼養(yǎng)牟氏螺螄腸道微生物組成和多樣性的差異，旨在揭示螺螄腸道菌群的形成機(jī)制，為螺螄的保護(hù)、馴養(yǎng)提供一定的參考依據(jù)。

1 材料與方法

1.1 樣品搜集2021年6月，從云南省昆明市的滇池水域采集并鑒定了8個(gè)活體牟氏螺螄。選取其中4個(gè)進(jìn)行解剖，將其腸道內(nèi)容物洗脫至5 mL無(wú)菌水中，-80 ℃下冷凍保存，作為野生組(W組)用于后續(xù)測(cè)序分析；其余4個(gè)牟氏螺螄放置在水箱中，使用25 ℃滇池水，以蒸煮后的海菜花和南瓜以及市售硅藻粉和螺旋藻粉為食物，在通氣條件下喂養(yǎng)2個(gè)月，然后進(jìn)行解剖，作為飼養(yǎng)組(C)，取樣方法與W組相同。具體樣品信息見(jiàn)表1。

表1 采樣信息Table 1 The sample information

1.2 DNA提取及3代全長(zhǎng)擴(kuò)增子測(cè)序含有腸道內(nèi)容物的液體混勻后，吸取250 μL，使用ZymoBIOMICS DNA Miniprep Kit進(jìn)行總DNA提取。提取后的DNA使用0.8%瓊脂糖凝膠電泳檢測(cè)DNA完整性，隨后使用Tecan F200進(jìn)行核酸濃度檢測(cè)(PicoGreen染料法)。質(zhì)量檢測(cè)合格的DNA使用引物對(duì)8F(5′-AGAGTTTGATCATGGCTCAG-3′)和1492R(5′-CGGTTACCTTGTTACGACTT-3′)進(jìn)行全長(zhǎng)16S rRNA基因擴(kuò)增。PCR程序如下：預(yù)變性98 ℃ 30 s，1個(gè)循環(huán)；98 ℃變性5 s，54 ℃退火15 s，72 ℃延伸45 s，28個(gè)循環(huán)；72 ℃ 2 min，1個(gè)循環(huán)；4 ℃保溫。PCR產(chǎn)物使用1% 瓊脂糖凝膠進(jìn)行目的片段電泳檢測(cè)。檢測(cè)合格的樣品取目的條帶使用Zymoclean Gel Recovery Kit(D4008)進(jìn)行回收，并使用Tecan F200定量；將所有樣品的待測(cè)產(chǎn)物等摩爾量混合后，使用Nanopore Ligation Sequencing Kit(SQK-LSK109)以及Nanopore Flow Cell Priming Kit(EXP-FLP002)構(gòu)建測(cè)序文庫(kù)。最后，使用Nanopore GridION Mk1測(cè)序儀及MinION Flow Cell R9.4.1芯片，送至成都羅寧生物科技有限公司完成3代測(cè)序工作，單個(gè)樣本測(cè)序量保證3萬(wàn)條以上。

1.3 數(shù)據(jù)分析原始Fast5序列使用Guppy軟件進(jìn)行堿基解析(base calling)，獲得Fastq格式的序列文件。Fastq文件使用qcat進(jìn)行樣品拆分。對(duì)拆分后的數(shù)據(jù)進(jìn)行質(zhì)量控制，保留序列平均質(zhì)量 >Q10、長(zhǎng)度1 300～1 600 bp的序列。使用Kraken2軟件結(jié)合SILVA數(shù)據(jù)庫(kù)v138對(duì)質(zhì)量控制后的序列進(jìn)行物種注釋，并合并每個(gè)物種的序列數(shù)目，構(gòu)建物種豐度表(OTU表)。隨機(jī)選取每個(gè)物種內(nèi)的一條序列作為代表性序列，使用MUSCLE進(jìn)行多重比對(duì)，隨后使用FastTree構(gòu)建Maximum-likelihood系統(tǒng)發(fā)育樹(shù)。

使用R語(yǔ)言V4.1.0進(jìn)行數(shù)據(jù)統(tǒng)計(jì)分析。使用vegan包函數(shù)diversity計(jì)算樣品的α多樣性指數(shù)Chao1和Simpson。α多樣性的組間Kruskal-Wallis秩和檢驗(yàn)使用kruskal.test函數(shù)。使用GUniFrac包分別計(jì)算Weighted(豐度加權(quán))UniFrac和Unweighted(非豐度加權(quán))UniFrac β多樣性矩陣。NMDS(非度量多維尺度)分析使用vegan包中的metaMDS函數(shù)。β多樣性矩陣的組間ANOSIM分析使用vegan包中的anosim函數(shù)進(jìn)行。random forest差異物種分析使用random forest包。

2 結(jié)果與分析

2.1 野生和飼養(yǎng)牟氏螺螄的腸道菌群群落組成從群落組成角度，該研究重點(diǎn)關(guān)注了門水平和種水平的群落組成及組間差異。從門水平來(lái)看，野生組和飼養(yǎng)組的主要物種組成差異不大(圖1)。相對(duì)豐度在1%以上的門為Proteobacteria、Firmicutes和Bacteroidetes。從種水平來(lái)看，主要使用韋恩圖分析了組間獨(dú)有和共有物種的組成情況。2組間共有物種670個(gè)，這些物種占總體物種豐度的96.97%。 W組獨(dú)有物種為262個(gè)、C組獨(dú)有物種為159個(gè)，分別占總體物種豐度的1.93%和1.10%。由此可見(jiàn)，野生和飼養(yǎng)牟氏螺螄的微生物組成總體區(qū)別不大，雖然2組存在多種稀有物種，但占比極低。由此可見(jiàn)，野生牟氏螺螄在人工飼養(yǎng)馴化過(guò)程中，核心腸道菌群保持了相對(duì)的魯棒性。

圖1 野生(W組)和人工飼養(yǎng)(C組)牟氏螺螄門水平(A)和種水平(B)組成分析Fig.1 The composition analysis of M. monodi at phylum level(A)and species level(B)in wild group (W) and culture group (C)

2.2 野生和飼養(yǎng)牟氏螺螄的腸道菌群多樣性差異群落多樣性分別從α多樣性和β多樣性2個(gè)角度進(jìn)行了分析。α多樣性表征的是樣品的多樣性指標(biāo)，圖2分別展示了Chao1指數(shù)和Simpson指數(shù)的多樣性結(jié)果?？傮w來(lái)看，2個(gè)α多樣性指數(shù)都表現(xiàn)出W組比C組有著更高的多樣性。其中Chao1指數(shù)W組和C組間存在顯著差異(P<0.05)。Simpson指數(shù)雖然2組間存在一定差異，但2組間差異不顯著。究其原因，Chao1指數(shù)對(duì)稀有物種較Simpson指數(shù)有更多的加權(quán)。由此可見(jiàn)，從α多樣性角度來(lái)看，W組和C組間的差異更多地體現(xiàn)在稀有物種的丟失導(dǎo)致的群落豐富度和均勻度的降低。

組間的多樣性差異使用β多樣性來(lái)表征，這里采用基于UniFrac距離矩陣的β多樣性NMDS排序分析(圖3)?；诩訖?quán)的UniFrac距離的NMDS分析(圖3A)中，W組和C組并沒(méi)有出現(xiàn)明顯的分隔；基于非加權(quán)的UniFrac距離的NMDS分析(圖3B)中，W組和C組存在明顯的各自成簇趨勢(shì)?；?種UniFrac距離矩陣的ANOSIM的定量差異統(tǒng)計(jì)分析也證實(shí)了這一結(jié)果：加權(quán)的UniFrac組間ANOSIM分析結(jié)果中r=0.281 3、P=0.135，表示不存在顯著差異；非加權(quán)的UniFrac組間ANOSIM分析結(jié)果中r=0.708 3、P=0.031，P<0.05表示存在顯著差異。這表明群落構(gòu)建的組間差異主要是由于物種組成的變化而不是物種豐度的變化造成的。

圖2 野生(W組)和人工飼養(yǎng)(C組)牟氏螺螄腸道菌群α多樣性比較Fig.2 Alpha diversity comparison of intestinal flora in M. monodi between wild group(W) and culture group(C)

注：A.Weighted UniFrac；B.Unweighted UniFrac。圖3 基于UniFrac距離的NMDS排序圖Fig.3 NMDS ranking based on UniFrac distance

上述基于多樣性的分析結(jié)果表明，野生和飼養(yǎng)牟氏螺螄的腸道微生物存在差異。這一差異是由于稀有菌群在馴化飼養(yǎng)過(guò)程中的丟失造成的。由于飼養(yǎng)過(guò)程中主要是食物產(chǎn)生了變化，因此推測(cè)造成這一結(jié)果的原因有2個(gè)：一是由于在飼養(yǎng)過(guò)程中食物是經(jīng)過(guò)高溫滅菌的，所以食物所攜帶的微生物變少，飼養(yǎng)螺螄失去了來(lái)自食物的微生物在其腸道中的持續(xù)定殖，因此腸道微生物呈現(xiàn)多樣性降低的過(guò)程。二是食物種類的變化對(duì)宿主腸道微生物有較大的影響。食物的類型會(huì)顯著影響宿主腸道微生物的多樣性，例如Bolnick等[11]研究發(fā)現(xiàn)三刺魚(yú)和河鱸的腸道微生物多樣性與食物的多樣性呈負(fù)相關(guān)，而牟氏螺螄的食物多樣性降低帶來(lái)了腸道菌群多樣性的降低。

2.3 野生和飼養(yǎng)牟氏螺螄的腸道菌群差異物種分析從組成和多樣性角度來(lái)看，W組和C組間的差異主要是由于低豐度物種的丟失造成的。該研究使用基于隨機(jī)森林的機(jī)器學(xué)習(xí)方法來(lái)鑒定存在組間差異的物種(圖4)，平均基尼指數(shù)下降度(Mean Decrease Gini)越大，說(shuō)明該物種對(duì)組間差異的貢獻(xiàn)度越大。該研究共篩選出12個(gè)存在組間差異的物種(P<0.05)，其中10個(gè)為相對(duì)豐度均低于1%的低豐度物種，其余2個(gè)物種的相對(duì)豐度也低于5%。這一結(jié)果進(jìn)一步驗(yàn)證了多樣性研究的結(jié)論。12個(gè)物種中8個(gè)物種在培養(yǎng)過(guò)程中被富集(C組相對(duì)豐度大于W組)，分別為Streptococcusinfantarius、Ottowiasp.oral taxon 894、Pseudomonassilesiensis、Pseudomonasprotegens、Bacteroidesfragilis、Pseudomonassyringae、Actinobacilluspleuropneumoniae以及Bartonellasp.Raccoon 60。4個(gè)在培養(yǎng)過(guò)程中被選擇性過(guò)濾的物種分別為Mycoplasmamobile、Variovoraxsp.PMC12、Stenotrophomonasmaltophilia、Sphingomonassp.LK11。無(wú)論是富集還是被過(guò)濾掉的差異微生物，絕大多數(shù)都屬于變形菌門。這些微生物也是環(huán)境中常見(jiàn)的物種，通常分布在人、動(dòng)植物的生存環(huán)境中，很多是宿主的病原菌[12-14]。

組間差異物種中，M.mobile相對(duì)于其他差異物種而言比較特殊。雖然該菌在野生組中的平均相對(duì)豐度很低(0.27%)，但在人工飼養(yǎng)過(guò)程中直接消失。該菌是1987年首次被報(bào)道從魚(yú)腸道中分離得到的一種病原菌[15]。North等[16]通過(guò)高通量測(cè)序的方法發(fā)現(xiàn)了一種未分類的Mycoplasmasp.，該物種以極低的豐度存在于Biomphalariaglabrata、Cipangopaludinajaponica等田螺科物種的腸道中。Mycoplasmasp.是一類寄生在細(xì)胞間的微生物，其與軟體動(dòng)物的結(jié)合能力被認(rèn)為與動(dòng)物唾液酸的殘留有關(guān)[17-18]。由于該研究采用3代測(cè)序技術(shù)測(cè)定了全長(zhǎng)16S rRNA基因，因此在前人研究的基礎(chǔ)上進(jìn)一步確定了在牟氏螺螄中這一微生物為M.mobile。由于M.mobile在野生和人工飼養(yǎng)牟氏螺螄內(nèi)的分布特點(diǎn)，其可以作為區(qū)分牟氏螺螄來(lái)源的指征菌。綜上所述，差異物種分析驗(yàn)證了多樣性分析結(jié)果，低豐度的稀有物種是造成組間差異的重要原因。同時(shí)，筆者發(fā)現(xiàn)M.mobile在人工馴化飼養(yǎng)的過(guò)程中會(huì)完全消失，可以作為區(qū)分野生和飼養(yǎng)牟氏螺螄的指征微生物。

注：A.Weighted UniFrac；B.Unweighted UniFrac。圖4 基于隨機(jī)森林方法的組間差異物種分析Fig.4 Differential species analysis between different groups based on random forest method

3 結(jié)論

牟氏螺螄作為滇池特有的螺螄屬物種和極危物種，具有很高的研究?jī)r(jià)值。該研究首次通過(guò)使用3代全長(zhǎng)16S rRNA基因高通量測(cè)序，對(duì)野生和飼養(yǎng)牟氏螺螄的腸道微生物進(jìn)行了比較。該研究發(fā)現(xiàn)野生和飼養(yǎng)牟氏螺螄腸道的核心菌群相似性很高，主要由Proteobacteria、Firmicutes、Bacteroidetes等門組成。多樣性和差異物種分析表明，占總體物種豐度較小的稀有物種在人工飼養(yǎng)過(guò)程中的丟失是造成牟氏螺螄腸道菌群組間差異的主要原因。此外，M.mobile可以作為野生和飼養(yǎng)牟氏螺螄的指征菌，在飼養(yǎng)過(guò)程中該菌會(huì)消失。筆者首次報(bào)道了牟氏螺螄在野生和人工飼養(yǎng)過(guò)程中腸道微生物的組成及多樣性的變化，為牟氏螺螄的保護(hù)和繁育提供了一定的理論依據(jù)。