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

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

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

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

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

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

1 材料和方法

1.1 取樣

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

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

1.2 仿刺參樣品準備

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

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

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

1.4 PCR擴增與純化

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

1.5 文庫構建和上機測序

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

2 結果和分析

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

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

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

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

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

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

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

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

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

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

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

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

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

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

3 討論

3.1 仿刺參體表微生物的結構特征

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

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

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

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

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

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

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

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

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

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

4 結論

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

參考文獻

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

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

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

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

［5］黃華偉，王印庚.海參養(yǎng)殖的現(xiàn)狀、存在問題與前景展望［J］.中國水產，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］農業(yè)農村部漁業(yè)漁政管理局，全國水產技術推廣站，中國水產學會.中國漁業(yè)統(tǒng)計年鑒［M］.北京：中國農業(yè)出版社，2023.

［8］姜森顥，任貽超，唐伯平，等.我國刺參養(yǎng)殖產業(yè)發(fā)展現(xiàn)狀與對策研究［J］.中國農業(yè)科技導報，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］.水產學報，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］郝佳慧，蔡偉杰，李柯懋，等.高通量測序分析小頭裸裂尻魚皮膚和腸道的微生物多樣性［J］.微生物學報，2023，63（1）：233-250.

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

［23］程慶，江超，陳蘇芳，等.高通量測序分析湖南汝城縣熱水圩熱泉微生物多樣性［J］.湘南學院學報，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］白潔，李海艷，趙陽國.黃海北部不同站位海洋細菌群落分布特征［J］.微生物學報，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克隆文庫法分析患病刺參幼苗的菌群結構［J］.大連海洋大學學報，2014，29（3）：205-211.

［30］孫云，何明旺，張盼盼，等.美洲鰻鱺致病性鮑曼不動桿菌的分離、鑒定及致病性分析［J］.水產學報，2020，44（9）：1561-1571.

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

［32］陸文浩，陳輝，王習達，等.異育銀鯽致病性鮑曼不動桿菌的鑒定和系統(tǒng)發(fā)育分析［J］.中國獸醫(yī)科學，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］陳翠珍，房海，張曉君.石鰈病原瓊氏不動桿菌形態(tài)型Ⅰ的鑒定［J］.微生物學通報，2005，32（3）：34-39.

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

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

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

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

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

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

［43］孔祎頔，田佳鑫，陳秀梅，等.水產動物益生菌及其對腸道菌群影響的研究進展［J］.中國畜牧雜志，2019，55（1）：29-33.

［44］李建光.刺參菌群結構的分析及益生菌對刺參的影響［D］.大連：大連理工大學，2014.

［45］張濤，白嵐，李蕾，等.不同添加量的益生菌組合對仿刺參消化和免疫指標的影響［J］.大連水產學院學報，2009，24（Sup 1）：64-68.

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

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

［48］陽鋼，田相利，董雙林，等.四種不同添加物對仿刺參非特異性免疫力和腸道免疫基因Aj-p105、Aj-p50、Aj-rel和Aj-lys mRNA表達的影響［J］.水產學報，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)境與腸道內細菌菌群的研究［J］.海洋科學，2010，34（4）：64-69.

［52］李佳霖，汪光義，秦松.秦皇島近海養(yǎng)殖對潮間帶微生物群落多樣性的影響［J］.生態(tài)環(huán)境學報，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