基于高通量测序与培养基鉴定法分析不同省份不同养殖模式下小龙虾菌群的多样性

范鹏飞，冯 武，肖 瑶，张月洁，杜红英，陈加平

（华中农业大学食品科学与技术学院，湖北武汉 430070）

克氏原螯虾俗称小龙虾（Procambarus clarkii），原产北美地区，20世纪80年代起国内将之作为养殖对象引至各地，现分布在国内十多个省市[1]。近年来，小龙虾产业发展迅速，2020年全国小龙虾养殖总产量达239.37万吨，主要养殖省份沿长江分布，包括湖北、江苏、安徽等[2]。小龙虾的养殖及食品安全卫生问题在这几年随之增多，常见情况如养殖细菌病害[3]和横纹肌溶解综合征，2020年广州市白云区共发生25例小龙虾相关横纹肌溶解综合征[4]。由于小龙虾的市场逐渐增长，因此需要对小龙虾的养殖质量和在保藏运输中的存活与污染投入更多的关注，而这一过程中重要的影响因素就是小龙虾所携带的菌群情况[5-6]。目前，小龙虾养殖模式以稻虾互作为主，2020年稻虾互作占86.61%，小龙虾营养丰富且主要生活在稻塘水体环境中[7-8]。鲜活小龙虾自身携带多种微生物，极不耐运输储藏[9-10]。因此除应季小龙虾外，多以龙虾制品的形式进行运输，包括冻虾尾[11]、虾仁[12]和熟制小龙虾等形式[13-14]。小龙虾生长偏好湿热条件的环境，不同产地的气候、水土情况与饲养原料都存在一定差别，小龙虾出塘后携带的微生物也由此存在差别，需要对多种影响因素进行研究[15]。

培养基鉴定法是使用根据营养成分和抑制物质表现出选择性的培养基对样品中的微生物进行筛选，从而分离出多种细菌菌株，通过提取和扩增细菌DNA进行测序，使用NCBI Blast查询并确定细菌种类[16]。该方法在实施过程中主观影响较大，仅凭实验人员自主挑选菌株，容易出现重复与漏选的情况，影响菌群丰富度和菌种占比的测定。此外，存在多种细菌无法在培养基上进行培养的情况，自然条件下生长的细菌可通过培养基筛选的仅占总数的1%[17]。

高通量测序法不存在受限于人为培养的缺点，作为一种高效且全面的测定方式，高通量测序法（二代测序法）所需时长远低于传统培养的时长，在较短的时间内能同时对上百万条DNA分子进行检测，由于不需要进行培养增菌，无法培养或丰度较低的微生物同样会被检出，相对培养基鉴定法更为全面[18]。现在已在分析肉制品和水生生物等的微生物组成有广泛应用，Li等[19]应用高通量技术分析得冷藏后期猪肉中的优势菌属为假单胞菌属、不动杆菌属和发光杆菌属。该方法同样被逐渐应用于鉴定水产品携带或腐败微生物方面，如Xing等[20]利用高通量测序法分析得毛尾鱼微生物群落的丰富度和多样性在冷藏条件下先下降后短暂上升后下降的动态变化。汤纯等[15]通过高通量测序法得浦口小龙虾携带的优势菌为柠檬酸杆菌和气单胞菌，太湖东山小龙虾携带优势菌为蜂房哈夫尼菌。

本研究使用培养基鉴定法和高通量测序法共同实验，对湖北荆州、江苏盱眙和安徽霍邱三地的小龙虾的微生物情况进行研究，其中荆州两组样本分别来自稻虾互作塘与养殖塘，分析鲜活小龙虾整体和虾尾肉部位的菌群组成，同时研究分析两种养殖模式下的养殖环境中的微生物，确认小龙虾携带的菌群情况，分析小龙虾菌群与养殖环境的关系，为小龙虾在流通过程中的微生物防治措施研究提供参考。

1 材料与方法

1.1 材料与仪器

本实验选取六月份可上市生鲜小龙虾 在荆州、盱眙和霍邱三个地区的当地养殖塘通过五点取样法取小龙虾，每个取样点选取大小相近（25±2 g）的5尾小龙虾冷藏打包，通过快递送至实验室备用，其中荆州小龙虾捕捞自稻虾塘与养殖池塘两种养殖模式，盱眙和霍邱小龙虾均取自养殖池塘。采用同样方法，使用灭菌封口袋对荆州稻虾塘和养殖池塘的土壤和水体进行取样，置于冰盒送至实验室备用。PCA菌落计数培养基、结晶紫中性红胆盐葡萄糖琼脂（VRBGA）、甘露醇氯化钠琼脂（MSA）、乳酸菌琼脂基础培养基（MRS）、假单胞菌琼脂基础培养基（CN琼脂）、LB肉汤（LB）、细菌DNA提取试剂盒山东青岛海博生物技术有限公司。

SW-CJ-1D型无菌操作台 浙江苏净净化有限公司；ME204E电子分析天平 梅特勒-托利多仪器（上海）有限公司；SHP-250型生化培养箱 上海精宏实验设备有限公司；Microfuge 20R高速冷冻离心机德国Beckman coulter公司；超纯水系统 上海和泰仪器有限公司。

1.2 实验方法

1.2.1 样品处理 参考Özogul等[21]的方法并略有改动。对三产地的四组小龙虾分别做出如下2种处理：a. 称取75 g鲜活小龙虾，均质破碎混匀后取任意25 g加入225 mL生理盐水得到小龙虾整体样品；b. 无菌条件下取出不带虾线的虾尾肉碾碎混匀，取25 g加入225 mL生理盐水均质充分得到小龙虾尾肉（不带肠）样品处理液。每组样品处理3个平行。对荆州两组土样均质后放入-80 ℃冷冻保藏，水样则用0.22 μm滤膜过滤，在滤膜上得到菌群样本，将滤膜放入-80 ℃冰箱保存。

1.2.2 菌落计数与优势菌种保藏 根据GB 4789.2-2016[22]测定每组样品整体与虾尾肉的菌落总数。将1.2.1节中的样品稀释液通过平板涂布的方式涂布于PCA琼脂表面培养30 ℃培养24 h后计数，同时涂布在VRBGA、CN琼脂、MSA和MRS琼脂平板上，30 ℃培养24 h后挑取菌落。在4种培养基平板上选取不同的菌落，在对应的培养基上进行划线培养直至纯化得单菌落，挑取单菌落至10 mL LB肉汤中，30 ℃，120 r/min振荡培养至菌液混浊（此时OD600nm=1）时，取1 mL菌液与0.2 mL无菌甘油至2 mL无菌离心管中混匀密封，放入-80 ℃冰箱保存。

1.2.3 样品总DNA提取 取1.2.1节中的小龙虾整体样品和虾尾肉均质液各50 mL于100 mL无菌离心管中，4 ℃，2000 r/min条件下离心10 min，取上清液在4 ℃，12000 r/min条件下再次离心10 min，倒掉上清液，取沉淀置于无菌袋中。挑取少量沉淀使用细菌基因组DNA提取试剂盒提取DNA进行检测，将检测合格的沉淀样品放于-80 ℃下保存。

1.2.4 高通量测序和数据分析 将保存在-80 ℃冰箱中的样品处理产物置于干冰箱中送至诺禾致源有限公司的Illumina PE100高通量测序平台进行检测。使用引物27F（5’-AGAGTTTGATCCTGGCTCAG-3’）和1492R（5’-TACGACTTAACCCCAATCGC-3’）对样品菌群进行扩增，对扩增产物在16S rDNA V3～V4区的DNA信息进行鉴定[23]。完成后用High Performance Computing平台处理数据，分析获取可用序列[24]。以各样品的序列作为基础单元，以相似性在98%以上为标准，将各个单元分类得操作分类单元（Operational Taxonomic Units，OTU）。将各个OTU比对High Performance Computing平台数据库，获得相应物种信息。Alpha多样性分析选用Mothur软件（version 1.31.2）。汇总比较三个产地四组生鲜小龙虾携带的菌群，统计微生物群落结构。

1.3 数据处理

方差分析使用SPSS 16.0软件，进行显著性分析。所有实验重复3次。

2 结果与分析

2.1 三地区小龙虾初始菌落数分析

图1为三个产地四组小龙虾初始菌落数。由图1a可知，荆州稻虾塘和养殖塘两种小龙虾整体样品相比江苏盱眙小龙虾整体样品、安徽霍邱小龙虾整体样品的菌落总数较少，分别为5.85 lg（CFU/g）和5.75 lg（CFU/g）；而后两者菌落总数相近，分别为6.32 lg（CFU/g）和6.35 lg（CFU/g）。由于在捕捞后盱眙和霍邱小龙虾运输至实验室所用时间相近且均高于荆州小龙虾运输时间，而同为池塘养殖的荆州小龙虾样本菌落总数显著低于前两者，因此推测可能菌落总数差别与运输时间与运输温度相关，这也说明活体小龙虾在运输中仍存在微生物菌群生长情况[25]。

从图1b可知，荆州稻虾塘小龙虾尾肉的菌落数为5.27 lg（CFU/g），荆州养殖塘小龙虾尾肉的菌落数为5.21 lg（CFU/g），盱眙地区小龙虾尾肉菌落数为5.51 lg（CFU/g），霍邱地区小龙虾尾肉菌落数为5.49 lg（CFU/g）。三地小龙虾尾肉初始菌落数均高于5 lg（CFU/g），这可能与小龙虾整体的菌落数有关，同时受到小龙虾出塘是在环境温度较高的6月的影响[26]。

结合图1a与图1b两图可以看出，小龙虾虾尾肉的菌落数相比整体样品略低。这可能是由于小龙虾表面的甲壳中含有致密的几丁质，可以减少环境微生物入侵体内，也可能是由于虾头与甲壳被去除而导致微生物含量的减少。

图 1 三地四组小龙虾整体和虾尾肉的菌落数差异Fig.1 Differences in bacterial counts of whole crayfish and tail meat of four groups of crayfish in three places

2.2 优势携带菌的分离与鉴定

挑取特征菌落在同类培养基上进行划线分离得单菌株，共分离209株菌。提取各个菌株的DNA后通过PCR扩增后检测产物的质量。扩增的样品的特异性条带在1500 bp附近，对扩增产物使用NCBI Blast进行比对，通过同源性高于98%的已知序列来确定菌株。209株菌株中包括荆州稻虾塘29种，荆州养殖塘30种，盱眙74种，霍邱76种。其中大多菌株属于肠杆菌科，包括柠檬酸杆菌、哈夫尼亚菌和粘质沙雷氏菌，有22株属于柠檬酸杆菌科。Pan等[27]从患病鲫鱼中提取出致病性弗氏柠檬酸杆菌。Liu等[28]发现有弗氏柠檬酸杆菌株对小龙虾存在高毒性会造成死亡。因此对柠檬酸杆菌这类条件致病菌，需要注意防范。

2.3 Alpha多样性指数分析

Alpha多样性指数（包括Chao1、PD whole tree、observed species及Shannon指数）分析通过不同指数的指标表现各个样品菌群的丰度和均匀度等数据。其中Chao1及observed species指数表现各个菌群的丰度水平（Species Richness）。而Shannon指数与菌群多样性（Species Diversity）呈正相关，PD whole tree指数表现各个样品对进化历史的留存差别，如图2所示。

由于Shannon指数反映出各样品微生物菌落的均匀度和丰度，受各样品微生物菌落的多样性的直接影响，因此本研究通过分析Shannon指数确认各样品的菌群多样性[29]。根据图2可以看出，荆州两处小龙虾整体和虾肉Shannon指数的色块位置相对低于另外两组。这表明荆州小龙虾整体和虾尾肉的菌群多样性均较低。而盱眙和霍邱两处小龙虾整体Shannon指数的色块位置都较高，霍邱样品的尾肉部分Shannon指数色块高于盱眙样品虾尾肉部分。这表明盱眙和霍邱的样品菌群多样性均高于荆州的两组样品，四组样本虾尾肉的微生物多样性又以霍邱的最多，盱眙样本虾尾肉相较略低。

而就养殖环境微生物多样性而言，荆州稻虾塘水样与土样的Shannon指数色块均高于荆州养殖塘对应的Shannon指数色块，即稻虾塘养殖环境的微生物多样性高于养殖塘。

2.4 物种组成和聚类分析

分类并归纳各个OTU的物种所属信息，对样品中的微生物在科和属两个水平进行菌株种属分类。

2.4.1 基于科水平的分析 由图3可得，三个产地的四组小龙虾样本的携带微生物菌落的主要菌株组成属于气单胞菌科（Aeromonadaceae）、肠杆菌科（Enterobacteriaceae）、莫拉菌科（Moraxellaceae）和丛毛单胞菌科（Comamonadaceae）4种菌科，不同产地小龙虾携带菌相有一定差别。气单胞菌科（Aeromonadaceae）存在于三个产地小龙虾菌群中并占据优势，但在荆州两组整虾样本中占比接近或低于肠杆菌科（Enterobacteriaceae）。而莫拉菌科（Moraxellaceae）和丛毛单胞菌科（Comamonadaceae）也是盱眙和霍邱两组小龙虾整体的优势菌群。

图 2 Alpha多样性指数的箱图Fig.2 The box plots of Alpha diversity indexes

就龙虾与虾尾肉的菌相关系来看，相对于在小龙虾整体菌群的占比，虾尾肉中肠杆菌科（Enterobacteriaceae）占比均较低，丛毛单胞菌科（Comamonadaceae）的占比则相近或相对较大。气单胞菌科（Aeromonadaceae）在小龙虾和虾尾肉中均为优势菌群。莫拉菌科（Moraxellaceae）也在虾尾肉中占一定比例。

就养殖环境的菌群来看，荆州稻虾塘的土壤菌群丰富度明显高于荆州养殖塘，而水体菌群组成则比较相似，优势菌占比较为均匀，与小龙虾和虾肉的菌群组成差别较大，其中丛毛单胞菌科（Comamonadaceae）是养殖环境和小龙虾都存在的优势菌种。

2.4.2 基于属水平的分析 从图4可以发现，各样品微生物多样性丰富，丰度较高的有气单胞菌属（Aeromonas）、乳球菌属（Lactococcus）、不动杆菌属（Acinetobacter）、金黄杆菌属（Chryseobacterium）、哈夫尼亚-肥胖杆菌（Hafnia-Obesumbacterium）、柠檬酸杆菌属（Citrobacter）、卟啉单胞菌属（Porphyromonas）和志贺氏菌属（Escherichia-Shigella）等，三产地的四组样品菌相大体类似，但菌群丰度差异大，其中荆州稻虾塘小龙虾整体携带菌为气单胞菌属、乳球菌属、柠檬酸杆菌属和志贺氏菌属；荆州养殖塘小龙虾整体携带菌主要有气单胞菌属、不动杆菌属、哈夫尼亚-肥胖杆菌属和乳球菌属；盱眙小龙虾整体携带菌为不动杆菌属、气单胞菌属、微小杆菌属（Exiguobacterium）和柠檬酸杆菌属；霍邱小龙虾整体携带菌主要有气单胞菌属、不动杆菌属、微小杆菌属和柠檬酸杆菌属。这些差异主要与养殖所在地区环境有关。可以看出，不同地区不同养殖模式下小龙虾携带的优势菌株较为一致，在菌群占比中存在差异，而菌群中其他菌株的差别可能受养殖环境和模式的影响存在一定差异，推测可能与所在水系有关，荆州两组小龙虾均处于长江干流地区水域，不同养殖模式导致菌种占比存在差异；而盱眙和霍邱则处于下游的支流淮河水系，因此两组小龙虾携带菌群相似度更高[30]。

图 3 基于科水平的菌种组成分析柱状图Fig.3 Column chart of bacterial composition analysis based on family level

图 4 基于属水平的菌种组成分析柱状图Fig.4 Column chart of bacterial composition analysis based on genus level

而小龙虾尾肉部分，三产地四组样本的携带菌均以气单胞菌属、不动杆菌属为主，菌群占比略有差别。其中荆州两组样本菌群组成比较接近，含有的微生物有气单胞菌属、不动杆菌属、金黄杆菌属及希瓦氏菌属（Shewanella），其中虾尾肉与小龙虾整体的菌群差别主要在黄金杆菌属和希瓦氏菌属的菌群占比上，而这两种菌属于小龙虾的常见腐败菌和条件致病菌，因此对虾尾肉产品，应更加注意这两种菌的防治。而盱眙和霍邱的两组样本小龙虾整体菌相与虾尾肉菌相相近，虾尾肉中的气单胞杆菌属的占比相对小龙虾整体则较高。本实验结果验证了汤纯等[15]的研究的结论：淡水小龙虾的携带菌主要包括希瓦氏菌和柠檬酸杆菌等。

2.5 样品菌群多样性变化的热图分析

将表格中数据的差异与信息在一张图中用颜色的不同与深浅进行表达即可得到热图，图5中对12组样品数据表达的信息以不同深浅的颜色进行呈现，较为直接清晰。图中共12组样品以携带菌的丰度差异进行了划分与对比的标准，可以看出，各组小龙虾样品携带的菌群：荆州稻虾塘小龙虾的柠檬酸杆菌属、乳球菌属和志贺氏菌属等丰度较高，而荆州养殖塘小龙虾样本的哈夫尼亚菌属丰度较高。盱眙和霍邱样品中微小杆菌属的丰度较荆州地区样品高。虾尾肉的情况，荆州两组样品中希瓦氏菌属的相对丰度比盱眙和霍邱地区样品高。其中希瓦氏菌被验证为是低温贮藏的大黄鱼中的优势腐败菌[31]。而腐败希瓦氏菌会引起南美白对虾黑变及品质劣变[32]。李婷婷等[33]发现哈夫尼亚菌是一种耐寒的腐败菌。样品携带的菌株的差异可能是水产品的养殖方式、地理环境和以及运输方式等的因素导致优势携带菌的差异[34]。

图 5 OTU 及属水平热图Fig.5 OUT and heat map of genus level

3 结论

对不同产地不同养殖模式下的小龙虾进行研究，运用Illumina二代高通量测序法研究得出生鲜小龙虾整体和虾尾肉样品及其生长环境的菌相构成情况。

培养基鉴定法的结果表明不同来源小龙虾的主要携带菌种均包括柠檬酸杆菌、哈夫尼亚菌和粘质沙雷氏菌等。主要菌种多数为肠杆菌科，其中柠檬酸杆菌作为一类条件致病菌需要注意防范。

宏基因组测序显示，四组小龙虾所携带的优势菌种较为一致，但菌群比例存在一定差异，在科水平上主要属于气单胞菌科（Aeromonadaceae）、肠杆菌科（Enterobacteriaceae）、莫拉菌科（Moraxellaceae）和丛毛单胞菌科（Comamonadaceae）4种菌科。在属水平上，四组小龙虾携带菌中均存在气单胞菌属和不动杆菌属等菌属，具有较高一致性，但也存在一定差异。荆州不同养殖模式的两组小龙虾携带优势菌以气单胞菌属、乳球菌属和哈夫尼亚菌属为主，其中稻虾塘小龙虾的优势菌株占比较高；盱眙和霍邱两地小龙虾主要携带菌株较为相似，以气单胞杆菌属和不动杆菌属为主，出现这些差异可能是养殖模式、地理位置和所在水系等因素造成的。

此外，小龙虾整体的菌群多样性高于虾尾肉，这可能是因为整体菌群包括虾尾肉中的菌群，且在生长运输过程中受到环境的影响，因此相对虾尾肉菌群更为丰富，而虾尾肉中优势菌的占比则更大，四组小龙虾的虾尾肉携带菌均以气单胞菌属和不动杆菌属为主，在虾尾肉的加工和保鲜中应重点关注。

而不同养殖模式下的小龙虾的菌群，稻虾互作模式下，土壤微生物种类更丰富更均衡，养殖塘土壤微生物种类略低，优势菌占比较为明显。而两种模式下水体微生物菌群存在差异较小。

本研究对不同地区不同养殖环境下生鲜小龙虾及加工品携带微生物的进行研究，对携带菌的情况进行了探索，为携带优势菌种的防治等研究提供理论依据。