我國(guó)沿海5個(gè)地區(qū)冬季養(yǎng)殖仿刺參體表微生物多樣性比較分析

潘美利 黨慧鳳 黃宇希 王曉航 黎睿君 丁君

摘 要：為探究仿刺參（Apostichopus japonicus）體表微生物在仿刺參疾病發(fā)生與環(huán)境交互中的作用，采用Illumina HiSeq高通量測(cè)序技術(shù)對(duì)來自大連、營(yíng)口、唐山、威海、寧德等5個(gè)沿海地區(qū)的養(yǎng)殖仿刺參體表微生物進(jìn)行了多樣性分析。結(jié)果顯示，地域差異對(duì)仿刺參體表微生物組成影響明顯，其中福建省寧德市霞浦縣（ND）養(yǎng)殖仿刺參的體表微生物組成與其他4個(gè)地區(qū)差異最明顯。在門水平上，各地區(qū)的最優(yōu)勢(shì)菌門均為變形菌門（Proteobacteria），其中ND組相對(duì)豐富度（76.16%）最低；在目和科水平上，ND組的最優(yōu)勢(shì)微生物分別屬于柄桿菌目（Caulobacterales）和生絲單胞菌科（Hyphomonadaceae），所占比例分別為14.48%、14.34%；在屬水平上，ND組仿刺參體表微生物中不動(dòng)桿菌屬（Acinetobacter）、副球菌屬（Paracoccus）、毛螺菌屬（Lachnospira）等所占的比例與其他各組存在顯著差異（P＜0.05）；山東省威海市乳山市（WH）養(yǎng)殖仿刺參體表微生物多樣性在屬水平上的相對(duì)豐富度最高，該組副球菌屬所占比例與其他地區(qū)均存在顯著差異（P＜0.05）。研究結(jié)果可為進(jìn)一步開展仿刺參體表微生物多樣性與環(huán)境適應(yīng)進(jìn)化研究以及養(yǎng)殖仿刺參病害防治提供參考。

關(guān)鍵詞：仿刺參；冬季；體表微生物；多樣性差異

仿刺參（Apostichopus japonicus）又稱刺參，隸屬于棘皮動(dòng)物門（Echinodermata）、海參綱（Holothuroidea）、刺參科（Stichopodidae）、仿刺參屬（Apostichopus），是我國(guó)海參的代表物種之一，也是我國(guó)重要的棘皮動(dòng)物養(yǎng)殖品種，其產(chǎn)地從北部的渤海灣和遼東半島到南部的南沙群島均有分布［1-2］。仿刺參富含7種必需氨基酸和大量抗衰老的酸性黏多糖，被認(rèn)為是最好的海參品種之一［3-4］。為滿足人們的消費(fèi)需求，我國(guó)自20世紀(jì)50年代開始開展刺參增養(yǎng)殖技術(shù)研究［5］，其養(yǎng)殖在20世紀(jì)80年代迅速發(fā)展，仿刺參也逐漸成為我國(guó)最重要的水產(chǎn)養(yǎng)殖品種之一［6］。我國(guó)仿刺參的大規(guī)模養(yǎng)殖地區(qū)以遼東半島和山東半島為主，其次是河北和福建的沿海地區(qū)。根據(jù)農(nóng)業(yè)農(nóng)村部的數(shù)據(jù)統(tǒng)計(jì)，2022年中國(guó)仿刺參養(yǎng)殖產(chǎn)量達(dá)到248 508 t［7］。

隨著仿刺參養(yǎng)殖規(guī)模和范圍的不斷擴(kuò)大，仿刺參養(yǎng)殖過程中病害頻頻發(fā)生，對(duì)其養(yǎng)殖業(yè)造成了嚴(yán)重經(jīng)濟(jì)損失，其中以腐皮綜合征、幼參爛胃病、體泡綜合征（BVS）等疾病尤為嚴(yán)重［8-9］。Becker等［10］發(fā)現(xiàn)，引起仿刺參表皮潰爛病的病原有弧菌屬（Vibrio sp.）、擬桿菌屬（Bacteroides sp.）和α-變形菌門（α-Proteobacterium）等。王印庚等［11］認(rèn)為，引起養(yǎng)殖刺參腐皮綜合征的主要病原菌為弧菌和假單胞菌。張春云等［12］利用形態(tài)學(xué)、生理生化實(shí)驗(yàn)和16s rRNA分子生物學(xué)鑒定方法，發(fā)現(xiàn)養(yǎng)殖刺參腐皮綜合征的致病菌是燦爛弧菌（Vibrio splendidus）。這些疾病的發(fā)病部位大多在仿刺參體表，因此，研究仿刺參體表微生物的多樣性對(duì)于預(yù)防和控制養(yǎng)殖仿刺參疾病至關(guān)重要。

本研究對(duì)來自中國(guó)沿海5個(gè)養(yǎng)殖地區(qū)的冬季養(yǎng)殖仿刺參體表的微生物多樣性進(jìn)行了分析，以期為進(jìn)一步開展仿刺參體表微生物多樣性與環(huán)境適應(yīng)進(jìn)化研究以及仿刺參疾病防治提供參考。

1 材料和方法

1.1 取樣

2021年11月—2022年1月，分別在我國(guó)沿海5個(gè)大型半開放仿刺參養(yǎng)殖基地采集仿刺參樣品，采樣點(diǎn)具體信息見表1。各基地單個(gè)養(yǎng)殖池塘面積約為6.67 hm2，水深2～3 m，養(yǎng)殖水來自近海海域利用納潮納入的海水。全年不投喂飼料，不投放益生菌，保證仿刺參自然生長(zhǎng)。

從5個(gè)地區(qū)的養(yǎng)殖基地各隨機(jī)采集10頭仿刺參樣品，樣品平均濕體質(zhì)量在55～60 g/頭（排水后體質(zhì)量為35～40 g/頭）。將仿刺參從養(yǎng)殖基地的養(yǎng)殖池中撈起后立即放入加有養(yǎng)殖海水的無菌密封袋中，采用冰袋低溫運(yùn)輸，24 h內(nèi)空運(yùn)至大連市海珍品疾病防控重點(diǎn)實(shí)驗(yàn)室。

1.2 仿刺參樣品準(zhǔn)備

從各采樣點(diǎn)的仿刺參樣品中分別隨機(jī)選取3頭作為3個(gè)平行樣。在超凈工作臺(tái)上，使用滅菌的手術(shù)刀取仿刺參背部中間的表皮0.2 g放入滅菌的小燒杯中，用無菌海水濯洗2～3次，裝入滅菌的2 mL離心管中；在無菌條件下，使用勻漿器將離心管中的仿刺參表皮充分研磨成勻漿，將裝有勻漿的2 mL離心管放入液氮中速凍；使用干冰送交北京諾禾致源科技股份有限公司進(jìn)行擴(kuò)增子高通量測(cè)序分析。

1.3 仿刺參體表微生物基因組提取及檢測(cè)

采用CTAB法提取仿刺參體表細(xì)菌基因組DNA［13］。取適量樣本DNA于離心管中，用無菌水將樣本稀釋至1 ng/μL，利用瓊脂糖凝膠電泳檢測(cè)DNA的純度和濃度。

1.4 PCR擴(kuò)增與純化

分別以DL、YK、WH、TS、ND 5個(gè)地區(qū)的仿刺參體表DNA為模版，用16s rDNA-V4通用引物（上游序列為515F：5-GTGCCAGCMGCCGCGGTAA-3，下游序列為806R：5-GGACTACHVGGGTWTCTAAT-3）［14］對(duì)16s rDNA-V4的基因進(jìn)行PCR擴(kuò)増。取5 μL的PCR產(chǎn)物于2％的瓊脂糖凝膠中電泳檢測(cè)，檢測(cè)條件為5 V/cm，25 min，并用Qiagen Gel Extraction Kit（QIAGEN，Germany）凝膠回收試劑盒回收16S rDNA-V4區(qū)的片段。

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

測(cè)序所得的原始數(shù)據(jù)經(jīng)過FLASH 1.2.7軟件拼接［15］，拼接得到的原始序列數(shù)據(jù)（raw tags）利用Qiime V1.9.1進(jìn)行質(zhì)量過濾［16］，得到高質(zhì)量的標(biāo)簽數(shù)據(jù)（clean tags），通過UCHIME Algorithm去除嵌合體［17］，得到最終的有效數(shù)據(jù)（effective tags）。利用Uparse 7.0.1001軟件［18］對(duì)所有樣本的全部有效數(shù)據(jù)進(jìn)行聚類，以97%的一致性（identity）將序列聚類成為操作分類單元（operational taxonomic units，OTUs），同時(shí)選取OTUs中出現(xiàn)頻次最高的序列作為代表序列。用Mothur軟件和SILVA［19］的SSUrRNA數(shù)據(jù)庫［18］對(duì)OTUs代表序列進(jìn)行物種注釋分析（設(shè)定閾值為0.8～1.0）。此過程委托北京諾禾致源科技股份有限公司完成。

2 結(jié)果和分析

2.1 序列拼接組裝與OTU聚類分析

基于Illumina HiSeq測(cè)序平臺(tái)測(cè)序構(gòu)建PCR-free文庫，然后進(jìn)行雙末端（paired-end）測(cè)序。通過對(duì)讀長(zhǎng)（reads）拼接，平均每個(gè)樣品測(cè)得88 391條tags，經(jīng)過質(zhì)控，平均得到81 547條有效數(shù)據(jù)，質(zhì)控有效數(shù)據(jù)75 292條，質(zhì)控有效率達(dá)85.12%。以97%的一致性將序列聚類成為OTUs，共得到6 319個(gè)OTUs，然后對(duì)OTUs序列與Silva 132數(shù)據(jù)庫進(jìn)行物種注釋。其中能夠注釋到數(shù)據(jù)庫的OTUs數(shù)目為6 300，分類為41個(gè)門、63個(gè)綱、158個(gè)目、315個(gè)科、865個(gè)屬、606個(gè)種。每個(gè)組別分析后的詳細(xì)數(shù)據(jù)見表2。

2.2 仿刺參體表微生物Alpha多樣性分析

Alpha多樣性分析可反映樣本內(nèi)的微生物群落多樣性。冬季沿海5個(gè)地區(qū)仿刺參體表微生物多樣性指數(shù)見表3。由表3可見，5個(gè)地區(qū)養(yǎng)殖仿刺參體表的微生物多樣性存在差異。5組樣品的覆蓋率（goods_coverage）均達(dá)到99%以上，說明樣品中序列檢測(cè)率很高，測(cè)序結(jié)果符合樣品中體表微生物的真實(shí)情況。對(duì)比樣品中香農(nóng)（Shannon）指數(shù)，辛普森（Simpson）指數(shù)，DL組的物種多樣性最高；同時(shí)，DL組的豐富度（Chao1）指數(shù)最大，表明DL組物種豐富度最大。

2.3 仿刺參體表微生物Beta多樣性分析

Beta多樣性可以反映不同地區(qū)冬季仿刺參體表微生物物種的組成差異。采用非度量多維尺度分析（NMDS）方法對(duì)不同組間樣本物種差異進(jìn)行分析，結(jié)果見圖1。樣品間距離越近，表示二者物種差異程度越??；樣品間距離越遠(yuǎn)，表示二者物種差異越大。結(jié)果顯示，DL組和YK組的仿刺參體表微生物物種組成與其他3個(gè)地區(qū)仿刺參體表微生物物種組成差異較大；TS、WH以及ND組這3個(gè)地區(qū)的仿刺參體表微生物物種組成差異較小。

2.4 仿刺參體表微生物組成分析

在門水平上（見圖2-a），5個(gè)地區(qū)養(yǎng)殖仿刺參體表微生物的第一優(yōu)勢(shì)菌門均為變形菌門（Proteobacteria），其中TS組相對(duì)豐度占比最高，達(dá)87.44%；ND組相對(duì)豐度占比最低，為76.16%。YK、WH、TS組這3個(gè)地區(qū)養(yǎng)殖仿刺參體表微生物的次優(yōu)勢(shì)菌門均為厚壁菌門（Firmicutes），相對(duì)豐度占比分別為5.67%、8.94%、5.90%。

在綱水平上（見圖2-b），α-變形菌綱（Alphaproteobacteria）是5個(gè)地區(qū)養(yǎng)殖仿刺參體表微生物的第一優(yōu)勢(shì)菌綱，其中YK組相對(duì)豐度占比最高，達(dá)83.05%；DL組相對(duì)豐度占比最低，為52.16%。DL、WH、TS和ND組這4個(gè)地區(qū)養(yǎng)殖仿刺參體表微生物的次優(yōu)勢(shì)菌綱均為γ-變形菌綱（Gammaproteobacteria），相對(duì)豐度占比分別為22.82%、17.12%、9.42%和20.25%；而YK組養(yǎng)殖仿刺參體表微生物的次優(yōu)勢(shì)菌綱為擬桿菌綱（Bacteroidia）。

在目水平上（見圖2-c），DL、YK、WH和TS組養(yǎng)殖仿刺參體表微生物的第一優(yōu)勢(shì)菌目均為紅螺菌目（Rhodobacterales），相對(duì)豐度占比分別為15.58%、38.03%、22.72%和22.09%；而ND組的第一優(yōu)勢(shì)菌目為柄桿菌目（Caulobacterales），相對(duì)豐度占比為14.48%。在目水平上的第一優(yōu)勢(shì)菌目呈現(xiàn)出明顯的南北地區(qū)差異。DL、ND組養(yǎng)殖仿刺參體表微生物的次優(yōu)勢(shì)菌目為纖維弧菌目（Cellvibrionales），其相對(duì)豐度占比分別為14.58%和12.93%。

在科水平上（見圖2-d），DL、YK、WH和TS組養(yǎng)殖仿刺參體表微生物的第一優(yōu)勢(shì)菌科均為紅螺菌科（Rhodobacteraceae），相對(duì)豐度占比分別為15.58%、38.03%、22.72%和22.09%，而ND組的第一優(yōu)勢(shì)菌科為生絲單胞菌科（Hyphomonadaceae），相對(duì)豐度占比為14.34%。在目、科分類水平上的第一優(yōu)勢(shì)菌均呈現(xiàn)出明顯的南北地區(qū)差異。DL、TS組兩地養(yǎng)殖仿刺參體表微生物的次優(yōu)勢(shì)菌科為生絲單胞菌科，其相對(duì)豐度占比分別為11.17%和12.48%。

在屬水平上（見圖2-e），鑒定到的細(xì)菌總體水平偏低，常見病原菌弧菌屬并不在排名前30的優(yōu)勢(shì)菌內(nèi)，輔證了近年來由弧菌引起的仿刺參發(fā)病率逐年降低，同年度采集的仿刺參病例同樣較少。

由圖3可見，在屬水平上，沿海5個(gè)地區(qū)仿刺參體表微生物物種的相對(duì)豐度依次為WH>TS>DL>ND>YK，其中，WH和TS組乳桿菌屬（Lactobacillus）相對(duì)豐度占比較高，分別為2.70%和2.62%；WH組雙歧桿菌屬（Bifidobacterium）相對(duì)豐度占比較高，為2.26%。

由圖4可見，位于南部的寧德地區(qū)與其他4個(gè)地區(qū)在屬水平上具有顯著差異的微生物有嗜冷菌屬（Colwellia）、不動(dòng)桿菌屬（Acinetobacter）、溶桿菌屬（Lysobacter）、副球菌屬（Paracoccus）、鞘氨醇單胞菌（Sphingomonas）、亞硝化螺菌屬（Nitrosospira）、獨(dú)島桿菌屬（Dokdonella）、熱單胞菌屬（Thermomonas）等。威海地區(qū)在屬水平上的微生物多樣性最為豐富，該地區(qū)與其他4個(gè)地區(qū)存在差異的微生物有不動(dòng)桿菌屬（Acinetobacter）、副球菌屬（Paracoccus）等，其中副球菌屬（Paracoccus）與其他4地均有顯著差異。

3 討論

3.1 仿刺參體表微生物的結(jié)構(gòu)特征

本研究發(fā)現(xiàn)，沿海5個(gè)地區(qū)的仿刺參體表微生物主要由變形菌門（76%～88%）和厚壁菌門（2%～10%）組成，其中變形菌門為第一優(yōu)勢(shì)菌，這與主要魚類體表微生物結(jié)構(gòu)分析的結(jié)果一致［20-21］。有研究發(fā)現(xiàn)，在門水平上，虹鱒（Oncorhynchus mykiss）［20］皮膚微生物以變形菌門和厚壁菌門為主，小頭裸裂尻魚（Herzensteinia microcephalus）［21］皮膚微生物的主要優(yōu)勢(shì)菌為放線菌門（Actinobacteria）和變形菌門。仿刺參屬于底棲生物，其生活的沉積泥中的最優(yōu)勢(shì)菌也是變形菌門［22］。變形菌門是細(xì)菌中最大的一個(gè)門，有著最龐大的微生物類群，在各種復(fù)雜的環(huán)境中都能存活，這類微生物可以參與環(huán)境中的碳循環(huán)，為宿主積累能量，與宿主的生理特性關(guān)系緊密相連［23-24］。仿刺參體表微生物在綱、目水平上的優(yōu)勢(shì)菌分別為α-變形菌綱以及紅螺菌目，均歸屬于變形菌門，其中紅螺菌目是α-變形菌綱中的主要目，這類微生物可以通過光合作用進(jìn)行生長(zhǎng)，還可參與海洋中的碳、氮循環(huán)［25］。厚壁菌門多數(shù)為革蘭氏陽性細(xì)菌，其分布于近海和前海的沉積物中［26］，與刺參活動(dòng)范圍相符［27］。

3.2 仿刺參體表微生物的潛在致病菌

本研究對(duì)5個(gè)不同沿海地區(qū)的仿刺參體表微生物進(jìn)行分析，發(fā)現(xiàn)仿刺參體表微生物存在不動(dòng)桿菌屬，該屬微生物可能會(huì)導(dǎo)致水質(zhì)污染和水生動(dòng)物疾病的發(fā)生。不動(dòng)桿菌屬是一類非發(fā)酵革蘭氏陰性桿菌，在自然環(huán)境的水體和土壤中廣泛存在，是一類毒力較低的條件致病菌［28］。本研究發(fā)現(xiàn)，ND與YK、DL、TS等組的顯著差異菌在屬水平上均有不動(dòng)桿菌。王軼南等［29］發(fā)現(xiàn)，在化皮參苗（Apostichopus japonicus）病灶組織微生物群落結(jié)構(gòu)中存在不動(dòng)桿菌屬，且其屬于第二優(yōu)勢(shì)菌。不動(dòng)桿菌屬作為重要的病原菌可致使魚類發(fā)生疾病，如鮑曼不動(dòng)桿菌對(duì)美洲鰻鱺（Anguilla rostrata）［30］、鱖（Sinipercachuatsi）［31］、異育銀鯽（Carassius auratus gibelio）［32］、斑點(diǎn)叉尾魚回（Ictalurus punctatus）［33］、線鱧（Channastriata）［34］、露斯塔野鯪（Labeo rohita）［35］均具有致病性；瓊氏不動(dòng)桿菌可引起石鰈（Kareius bicoloratus）［36］、錦鯉（Cyprinus carpio）［37］發(fā)生疾病。有研究表明，不動(dòng)桿菌在水生動(dòng)物處于健康狀態(tài)時(shí)不具有致病性［38］，但該菌的一個(gè)重要的感染途徑是在水生動(dòng)物體表損傷時(shí)進(jìn)行感染［［39］。不動(dòng)桿菌屬的一些菌種可以分解有機(jī)質(zhì)，產(chǎn)生氨、亞硝酸鹽和硝酸鹽等營(yíng)養(yǎng)物質(zhì)，為水生動(dòng)物生長(zhǎng)提供養(yǎng)分，但在水生動(dòng)物養(yǎng)殖密度過高時(shí)則會(huì)導(dǎo)致水質(zhì)污染，引起動(dòng)物疾病的發(fā)生［40-41］。

弧菌作為條件致病菌，只有當(dāng)宿主抵抗力降低或環(huán)境極度惡化時(shí)在養(yǎng)殖動(dòng)物體表大量分離。李曉宇等［42］某次在大連地區(qū)發(fā)現(xiàn)仿刺參腐皮綜合征的致病菌為溶藻弧菌（Vibrio alginolyticus），人工回歸感染試驗(yàn)結(jié)果顯示，浸浴感染沒有出現(xiàn)病癥，而腹部注射和體壁創(chuàng)傷浸浴感染則出現(xiàn)病癥并伴有死亡，證實(shí)仿刺參體表黏膜層可抵御外來病原弧菌的侵襲，且溶藻弧菌是一種機(jī)會(huì)致病弧菌。張春云等［12］從患有腐皮綜合征仿刺參的病灶處分離得到優(yōu)勢(shì)菌株——燦爛弧菌，人工回歸感染試驗(yàn)發(fā)現(xiàn)，該菌株引發(fā)的癥狀與自然發(fā)病的癥狀相同，并且對(duì)健康仿刺參有很強(qiáng)的致病及致死作用。在本研究中，在屬水平上被注釋為前30位的優(yōu)勢(shì)屬中并未見到弧菌屬，推測(cè)這與采集的樣本均為健康仿刺參有關(guān)。

3.3 仿刺參體表微生物的有益生物

在本研究中，WH組的優(yōu)勢(shì)菌為乳桿菌屬，此外該組的副球菌屬與其他地區(qū)仿刺參體表微生物均存在顯著差異，以上菌屬都是益生菌，推測(cè)威海地區(qū)的仿刺參疾病發(fā)生較少。乳桿菌屬中的乳酸菌（Lactic acid bacteria）、嗜酸乳桿菌（Lactobacillus acidophilus）等作為水產(chǎn)動(dòng)物的腸道益生菌可改善腸道內(nèi)環(huán)境，增強(qiáng)機(jī)體的免疫力［43］，其中，乳酸菌可作為一種生物屏障抑制病原弧菌的感染；嗜酸乳酸桿菌可分泌抑制物質(zhì)，與病原弧菌競(jìng)爭(zhēng)營(yíng)養(yǎng)從而抑制病原菌的生長(zhǎng)，也可降低病原弧菌的黏附率［44］；乳酸桿菌作為益生菌可以促進(jìn)刺參腸道中益生菌的生長(zhǎng)，增強(qiáng)其免疫能力，降低刺參感染疾病的風(fēng)險(xiǎn)［45］。仿刺參腐皮綜合征等皮膚病是由弧菌引起的，周海平等［46］研究表明，乳酸桿菌對(duì)病原弧菌有很好的抑制作用；雙歧桿菌屬次之，其屬中的兩歧雙歧桿菌（Bifdobacterium bidfium）是水產(chǎn)養(yǎng)殖中常用的益生菌［47］。

副球菌屬隸屬于變形菌門中的紅螺菌目、紅螺菌科，與本研究中在目水平上和科水平上注釋到的微生物相呼應(yīng)。副球菌屬的存在可提高刺參的免疫能力。有研究表明，副球菌屬中的馬氏副球菌（Paracoccus marcusii）作為刺參飼料的添加劑可以有效提高刺參體腔細(xì)胞的呼吸爆發(fā)活力，上調(diào)刺參體腔細(xì)胞非特異性免疫基因的表達(dá)量，增強(qiáng)刺參的非特異性免疫力［48］，還可促進(jìn)其生長(zhǎng)，提高刺參腸道和體腔細(xì)胞的免疫能力。在體外試驗(yàn)中發(fā)現(xiàn)，馬氏副球菌對(duì)引起刺參腐皮綜合征的病原菌燦爛弧菌具有抑制作用［49］。

3.4 不同海域仿刺參體表微生物的異同及其影響因素

本研究發(fā)現(xiàn)，不同海域仿刺參體表微生物的群落結(jié)構(gòu)存在差異，推測(cè)其與仿刺參生存海域環(huán)境的微生物組成有關(guān)。在物種組成上，采自黃海海域（DL和WH）仿刺參體表微生物的多樣性和豐富度總和均高于采自渤海海域（YK和TS）以及東海海域（ND）仿刺參體表微生物的多樣性和豐富度總和。據(jù)報(bào)道，仿刺參腸道微生物與其生活的環(huán)境密切相關(guān)，黃海海域的仿刺參腸道微生物多樣性顯著高于渤海的［50］。這與本研究仿刺參體表微生物多樣性的研究結(jié)果一致。

不同海域仿刺參體表微生物的物種多樣性可能與海域間的水溫有關(guān)。本研究樣品采集時(shí)，黃海海域水溫高于渤海海域水溫，其微生物多樣性也高于后者，然而，東海海域水溫高于黃海海域水溫，其微生物多樣性卻低于黃海海域。有研究指出，溫度是影響環(huán)境中微生物結(jié)構(gòu)重要的物理因子，溫度能夠加快微生物的生長(zhǎng)和繁殖速度，同樣也可以抑制菌群的生長(zhǎng)和繁殖速度，從而通過這種促進(jìn)或抑制作用改變環(huán)境中的菌群結(jié)構(gòu)［51-52］。

4 結(jié)論

仿刺參體表微生物多樣性分析發(fā)現(xiàn)，仿刺參的體表存在益生菌屬，也存在條件致病菌屬。健康仿刺參的體表微生物處于相對(duì)平衡狀態(tài)，養(yǎng)殖環(huán)境惡化會(huì)導(dǎo)致一些有害微生物大量增加進(jìn)而誘發(fā)疾病，甚至導(dǎo)致養(yǎng)殖生物體死亡。仿刺參與其共生微生物的關(guān)系還需進(jìn)一步探索。本研究分析了仿刺參體表微生物與自身環(huán)境的關(guān)聯(lián)性，可為進(jìn)一步開展仿刺參體表微生物多樣性與環(huán)境適應(yīng)進(jìn)化以及仿刺參疾病的防治提供參考。

參考文獻(xiàn)

［1］夏海峰，杜宗軍，陳冠軍.刺參腸道及養(yǎng)殖池塘底泥微生物多樣性的比較研究［J］.海洋湖沼通報(bào)，2015（4）：105-110.

［2］張啟宇，劉峰，陳英義，等.海參病害防治診斷專家系統(tǒng)的研究［J］.江蘇農(nóng)業(yè)科學(xué)，2017，45（18）：226-229.

［3］樊繪曾.海參：海中人參——關(guān)于海參及其成分保健醫(yī)療功能的研究與開發(fā)［J］.中國(guó)海洋藥物，2001，20（4）：37-44.

［4］常亞青，隋錫林，李俊.刺參增養(yǎng)殖業(yè)現(xiàn)狀、存在問題與展望［J］.水產(chǎn)科學(xué)，2006，25（4）：198-201.

［5］黃華偉，王印庚.海參養(yǎng)殖的現(xiàn)狀、存在問題與前景展望［J］.中國(guó)水產(chǎn)，2007（10）：50-53.

［6］ZHAO Z L，JIANG J W，PAN Y J，et al.Proteomic analysis reveals the important roles of alpha-5-collagen and ATP5β during skin ulceration syndrome progression of sea cucumber Apostichopus japonicus［J］.Journal of Proteomics，2018，175：136-143.

［7］農(nóng)業(yè)農(nóng)村部漁業(yè)漁政管理局，全國(guó)水產(chǎn)技術(shù)推廣站，中國(guó)水產(chǎn)學(xué)會(huì).中國(guó)漁業(yè)統(tǒng)計(jì)年鑒［M］.北京：中國(guó)農(nóng)業(yè)出版社，2023.

［8］姜森顥，任貽超，唐伯平，等.我國(guó)刺參養(yǎng)殖產(chǎn)業(yè)發(fā)展現(xiàn)狀與對(duì)策研究［J］.中國(guó)農(nóng)業(yè)科技導(dǎo)報(bào)，2017，19（9）：15-23.

［9］ZHAO Z L，ZHOU Z C，DONG Y，et al.Association of intestinal fungal communities with the body vesicular syndrome：an emerging disease of sea cucumber（Apostichopus japonicus）［J］.Aquaculture，2021，530：735758.

［10］BECKER P，GILLAN D，LANTERBECQ D，et al.The skin ulceration disease in cultivated juveniles of Holothuria scabra （Holothuroidea，Echinodermata）［J］.Aquaculture，2004，242（1/2/3/4）：13-30.

［11］王印庚，榮小軍，張春云，等.養(yǎng)殖刺參暴發(fā)性疾病——“腐皮綜合癥”的初步研究與防治［J］.齊魯漁業(yè)，2004，21（5）：44-47.

［12］張春云，王印庚，榮小軍.養(yǎng)殖刺參腐皮綜合征病原菌的分離與鑒定［J］.水產(chǎn)學(xué)報(bào)，2006，30（1）：118-123.

［13］LUTZ K A，WANG W Q，ZDEPSKI A，et al.Isolation and analysis of high quality nuclear DNA with reduced organellar DNA for plant genome sequencing and resequencing［J］.BMC Biotechnology，2011，11：54.

［14］WU L Y，WEN C Q，QIN Y J，et al.Phasing amplicon sequencing on Illumina Miseq for robust environmental microbial community analysis［J］.BMC Microbiology，2015，15（1）：125.

［15］MAGOC T，SALZBERG S L.FLASH：fast length adjustment of short reads to improve genome assemblies［J］.Bioinformatics，2011，27（21）：2957-2963.

［16］CAPORASO J G，KUCZYNSKI J，STOMBAUGH J，et al.QIIME allows analysis of high-throughput community sequencing data［J］.Nature Methods，2010，7（5）：335-336.

［17］EDGAR R C，HAAS B J，CLEMENTE J C，et al.UCHIME improves sensitivity and speed of chimera detection［J］.Bioinformatics，2011，27（16）：2194-2200.

［18］HAAS B J，GEVERS D，EARL A M，et al.Chimeric 16S rRNA sequence formation and detection in Sanger and 454-pyrosequenced PCR amplicons［J］.Genome Research，2011，21（3）：494-504.

［19］EDGAR R C.UPARSE：highly accurate OTU sequences from microbial amplicon reads［J］.Nature Methods，2013，10（10）：996-998.

［20］CHO K，USHIKI T，ISHIGURO H，et al.Altered microbiota by a high-fat diet accelerates lethal myeloid hematopoiesis associated with systemic SOCS3 deficiency［J］.iScience，2021，24（10）：103117.

［21］郝佳慧，蔡偉杰，李柯懋，等.高通量測(cè)序分析小頭裸裂尻魚皮膚和腸道的微生物多樣性［J］.微生物學(xué)報(bào)，2023，63（1）：233-250.

［22］李曉龍.正常和發(fā)病刺參養(yǎng)殖池塘沉積環(huán)境中細(xì)菌群落變化的研究［D］.保定：河北農(nóng)業(yè)大學(xué)，2014.

［23］程慶，江超，陳蘇芳，等.高通量測(cè)序分析湖南汝城縣熱水圩熱泉微生物多樣性［J］.湘南學(xué)院學(xué)報(bào)，2021，42（5）：22-27.

［24］LU H P，WANG Y B，HUANG S W，et al.Metagenomic analysis reveals a functional signature for biomass degradation by cecal microbiota in the leaf-eating flying squirrel（Petaurista alborufus Lena）［J］.BMC Genomics，2012，13：466.

［25］白潔，李海艷，趙陽國(guó).黃海北部不同站位海洋細(xì)菌群落分布特征［J］.微生物學(xué)報(bào)，2009，49（3）：343-350.

［26］RAVENSCHLAG K，SAHM K，PERNTHALER J，et al.High bacterial diversity in permanently cold marine sediments［J］.Applied and Environmental Microbiology，1999，65（9）：3982-3989.

［27］LIU S L，SUN J C，RU X S，et al.Differences in feeding，intestinal mass and metabolites between a thermotolerant strain and common Apostichopus japonicus under high summer temperature［J］.Aquaculture Research，2018，49（5）：1957-1966.

［28］ELIOPOULOS G M，MARAGAKIS L L，PERL T M.Acinetobacter baumannii：epidemiology，antimicrobial resistance，and treatment options［J］.Clinical Infectious Diseases，2008，46（8）：1254-1263.

［29］王軼南，王俊杰，王姣姣，等.用16S rDNA克隆文庫法分析患病刺參幼苗的菌群結(jié)構(gòu)［J］.大連海洋大學(xué)學(xué)報(bào)，2014，29（3）：205-211.

［30］孫云，何明旺，張盼盼，等.美洲鰻鱺致病性鮑曼不動(dòng)桿菌的分離、鑒定及致病性分析［J］.水產(chǎn)學(xué)報(bào)，2020，44（9）：1561-1571.

［31］顧天釗，陸承平，陳懷青.鮑氏不動(dòng)桿菌──鱖魚暴發(fā)性死亡的新病原［J］.微生物學(xué)通報(bào)，1997，24（2）：104-106.

［32］陸文浩，陳輝，王習(xí)達(dá)，等.異育銀鯽致病性鮑曼不動(dòng)桿菌的鑒定和系統(tǒng)發(fā)育分析［J］.中國(guó)獸醫(yī)科學(xué)，2009，39（4）：303-309.

［33］MARTOJA R，MARTOJA-PIERSON M.Initiation aux techniques de lhistologie animale［R］.Masson，1967.

［34］RAUTA P R，KUMAR K，SAHOO P K.Emerging new multi-drug resistant bacterial pathogen，Acinetobacter baumannii associated with snakehead Channa striatus eye infection［J］.Current Science，2011，101（4）：548-553.

［35］BEHERA B K，PARIA P，DAS A，et al.Molecular characterization and pathogenicity of a virulent Acinetobacter baumannii associated with mortality of farmed Indian Major Carp Labeo rohita（Hamilton 1822）［J］.Aquaculture，2017，471：157-162.

［36］陳翠珍，房海，張曉君.石鰈病原瓊氏不動(dòng)桿菌形態(tài)型Ⅰ的鑒定［J］.微生物學(xué)通報(bào)，2005，32（3）：34-39.

［37］陳翠珍，張曉君，房海，等.錦鯉中瓊氏不動(dòng)桿菌的分離與鑒定［J］.水生態(tài)學(xué)雜志，2009，30（1）：86-90.

［38］張涵，周濤，王巖.綜合養(yǎng)殖池塘中三角帆蚌和魚類腸道細(xì)菌的組成［J］.水生生物學(xué)報(bào)，2013，37（5）：824-835.

［39］吳春艷，任思宇，汪開毓，等.一株虹鱒源烏爾新不動(dòng)桿菌的分離鑒定與其感染的病理損傷［J］.水產(chǎn)科學(xué)，2019，38（6）：797-803.

［40］羅國(guó)芝，陳家捷，于文杰，等.一株新型異養(yǎng)硝化細(xì)菌處理養(yǎng)殖水的效果［J］.環(huán)境工程學(xué)報(bào)，2016，10（8）：4206-4212.

［41］張喜昌，薛明洋，戚曉舟，等.從海泥中分離的約氏不動(dòng)桿菌聚磷特性研究［J］.家畜生態(tài)學(xué)報(bào)，2022，43（5）：60-64.

［42］李曉宇，馮麗娟，徐樂，等.刺參腐皮綜合征病原菌的分離鑒定及特性分析［J］.水產(chǎn)科學(xué)，2022，41（5）：759-768.

［43］孔祎頔，田佳鑫，陳秀梅，等.水產(chǎn)動(dòng)物益生菌及其對(duì)腸道菌群影響的研究進(jìn)展［J］.中國(guó)畜牧雜志，2019，55（1）：29-33.

［44］李建光.刺參菌群結(jié)構(gòu)的分析及益生菌對(duì)刺參的影響［D］.大連：大連理工大學(xué)，2014.

［45］張濤，白嵐，李蕾，等.不同添加量的益生菌組合對(duì)仿刺參消化和免疫指標(biāo)的影響［J］.大連水產(chǎn)學(xué)院學(xué)報(bào)，2009，24（Sup 1）：64-68.

［46］周海平，李卓佳，楊鶯鶯，等.乳酸桿菌LH對(duì)幾種水產(chǎn)養(yǎng)殖病原弧菌的抑制作用［J］.臺(tái)灣海峽，2006，25（3）：388-395.

［47］劉廣，趙陽，李元莉，等.益生菌在水產(chǎn)無脊椎動(dòng)物養(yǎng)殖中的研究與應(yīng)用［J］.飼料與畜牧，2013（4）：46-49.

［48］陽鋼，田相利，董雙林，等.四種不同添加物對(duì)仿刺參非特異性免疫力和腸道免疫基因Aj-p105、Aj-p50、Aj-rel和Aj-lys mRNA表達(dá)的影響［J］.水產(chǎn)學(xué)報(bào)，2015，39（5）：638-647.

［49］YAN F J，TIAN X L，DONG S L，et al.Growth performance，immune response，and disease resistance against Vibrio splendidus infection in juvenile sea cucumber Apostichopus japonicus fed a supplementary diet of the potential probiotic Paracoccus marcusii DB11［J］.Aquaculture，2014，420/421：105-111.

［50］柴英輝，高菲，王金鋒，等.仿刺參（Apostichopus japonicus）腸道菌群的地域性差異與共性研究［J］.海洋與湖沼，2019，50（5）：1127-1137.

［51］李彬，榮小軍，廖梅杰，等.冬季刺參養(yǎng)殖環(huán)境與腸道內(nèi)細(xì)菌菌群的研究［J］.海洋科學(xué)，2010，34（4）：64-69.

［52］李佳霖，汪光義，秦松.秦皇島近海養(yǎng)殖對(duì)潮間帶微生物群落多樣性的影響［J］.生態(tài)環(huán)境學(xué)報(bào)，2011，20（5）：920-926.

Comparative analysis of microbial diversity on the skin of Apostichopus japonicus located at five coastal areas of China in winter

PAN Meili， DANG Huifeng， HUANG Yuxi， WANG Xiaohang， LI Ruijun， DING Jun

（College of Fisheries and Life Sciences，Dalian Ocean University，Dalian 116023，China）

Abstract： Apostichopus japonicus was an important aquatic echinoderm in China，and its surface microorganisms played an important role in the occurrence of diseases of A. japonicus and the interaction with the environment.In this study，Illumina HiSeq high-throughput sequencing technology was used to analyze the diversity of microorganisms on the skin of sea cucumbers in five different coastal areas of China（Dalian，Yingkou，Tangshan，Weihai and Ningde）.The results showed that the regional differences had a significant impact on the surface microbial composition of A. japonicus.The skin microorganisms of A. japonicus collected from Ningde（ND） had the most significant difference compared with the other 4 groups.At the phylum level，Proteobacteria was the most dominant phylum in all groups，with the lowest relative richness in ND group（76.16%）.At the order and family level，the most dominant microorganisms of ND group belonged to Caulobacterales and Hyphomonadaceae respectively，accounting for 14.48% and 14.34% respectively.At the genus level，the percentage of some microorganisms in ND group were significantly different from other 4 groups（P＜0.05），in terms of Acinetobacter，Paracoccus，Lachnospira，etc.In addition，Weihai（WH） group had the highest relative richness at the genus level，with significant difference in Paracoccus percentage（P＜0.05）.The results of four research would provide a reference for the further study on microbial diversity and environmental adaptation of sea cucumbers，as well as the prevention and treatment of diseases of sea cucumbers.

Key words： Apostichopus japonicus; winter; skin microorganisms; diversity difference