Antimicrobial spectrum, resistance gene detection and ERIC-PCR genotyping of Vibrio scophthalmi
-
摘要: 2020年6—10月,从辽宁地区患病大菱鲆 (Scophthalmus maximus) 体内共分离到394株菌株。基于16S rRNA序列鉴定分离株,采用微量稀释法分析随机挑选的18株大菱鲆弧菌 (Vibrio scophthalmi) 对8种抗生素的敏感性,检测其耐药基因携带情况,并基于ERIC-PCR进行分型研究。结果显示,分离株以弧菌属为主,共232株 (58.88%),其中大菱鲆弧菌117株 (29.70%)。经rpoD序列比对确定的18株大菱鲆弧菌对硫酸新霉素、氟甲喹和盐酸多西环素不耐药,对其他5种抗菌药物有不同的耐药率,多重耐药率为66.7%。18株大菱鲆弧菌共检测到酰胺醇类耐药基因floR (61.11%) 和cmlA (66.67%)、磺胺类耐药基因sul2 (55.56%)、喹诺酮类耐药基因qnrA (50%) 和qnrS (5.56%) 5种耐药基因,未检测到氨基糖苷类和四环素类耐药基因。酰胺醇类耐药基因与耐药表型符合率为88.89%,有一定的相关性。磺胺类耐药基因与耐药表型符合率次之 (50%),喹诺酮类符合率最低 (33.33%)。ERIC-PCR将该18株大菱鲆弧菌分为4个亚型,I型 (44.44%) 和III型 (44.44%) 为主要谱型,与耐药表型或耐药基因无明显相关性。根据药敏结果,防治辽宁地区大菱鲆弧菌建议首选盐酸多西环素。Abstract: From June to October of 2020, we had isolated 394 strains from diseased turbot (Scophthalmus maximus) in Liaoning Province, and identified them based on the 16S rRNA sequence. Then we randomly selected 18 Vibrio scophthalmi strains, which were further typed by enterobacterial repetitive intergenic consensus PCR (ERIC-PCR), so as to analyze the susceptibility to 8 antibiotics and the presence of drug-resistance genes. The results show that most of the isolates were Vibrio with 232 strains (58.88%), among which 117 strains (29.70%) were V. scophthalmi. The 18 strains identified by rpoD sequence alignment were not resistant to neomycin sulfate, flumequine and doxycycline hydrochloride, but resistant to other five antibiotics with different ratios. The multiple drug resistance rate was 66.7%. Five drug resistance genes in the 18 strains of V. scophthalmi were detected (Amide alcohol resistance genes floR of 61.11% and cmlA of 66.67%, sulfonamide resistance gene sul2 of 55.56%, quinolone resistance genes qnrA of 50% and qnrS of 5.56%). No aminoglycosides or tetracycline resistant genes were detected. The coincidence rate of amide alcohol resistance genes and drug resistance phenotype was 88.89%, indicating a certain correlation between them. Sulfonamide resistance genes had a coincidence rate of 50% with the resistance phenotype, and quinolones resistance genes had the lowest coincidence rate of 33.33%. Finally, we clustered the 18 strains into four types by using ERIC-PCR fingerprinting. Type I (44.44%) and Type III (44.44%) were the main clusters. We observed apparent absence of correlation between ERIC-PCR fingerprinting and the resistance phenotype or resistance genes. Thus, it is conluded that doxycycline hydrochloride can be the first choice for the control of V. scophthalmi in Liaoning aquaculture farms.
-
Keywords:
- Vibrio scophthalmi /
- Drug resistance /
- Microdilution method /
- Drug resistance genes /
- ERIC-PCR
-
少鳞鱚 (Sillago japonica)是鲈形目、鱚科、鱚属的一种小型海洋经济鱼类,主要分布在红海、印度洋、太平洋西部和南部,盛产于我国南海和东海的浅水区或河口区[1-2]。少鳞鱚肉质鲜美、营养丰富,深受人们喜爱,也是近岸渔业捕捞对象和游钓鱼种。20世纪末开始,日本学者围绕少鳞鱚渔业管理[3]、繁殖发育学[4-7]、组织胚胎学[8-9]和资源生物学[10-11]领域开展了大量基础性研究工作,近年来国内有关其形态学和遗传学方面的研究主要集中在外部形态特征描述[12]、耳石形态多元统计分析[13]、分子系统学研究[14-15]、线粒体控制区[16]和微卫星分子标记开发[17]。遗传多样性反映了生物对环境变化的适应能力,是物种多样性的基础和决定性因素,体现在表型、染色体、蛋白质和DNA多个层次水平的变异上,开展遗传多样性研究对于探讨物种的系统进化、遗传分化和资源保护都具有十分重要的意义[18]。
线粒体DNA (Mitochondrial DNA,mtDNA)具有分子量小、结构简单、排列紧凑,几乎不发生重组且遵循母系遗传的特性,目前已广泛应用在鱼类分子系统学、生物地理学和种群遗传学等诸多领域[19-20]。在mtDNA包含的13个蛋白质编码基因中,NADH脱氢酶亚基2 (NADH dehydrogenase subunit 2,ND2)基因由于进化速率相对较快,能较好地反映种内与群体间的遗传变异及系统分类地位,正逐渐应用于鱼类系统发育和种群划分研究中[21-25]。本研究采集了中国沿海6个地理种群的少鳞鱚样本,基于线粒体ND2基因序列开展其群体遗传多样性的研究,以期为少鳞鱚种质资源的保护提供理论依据。
1. 材料与方法
1.1 实验材料
本研究所用少鳞鱚样本分别于2017—2018年间采自莱州(LZ,19尾)、胶南(JN,20尾)、舟山(ZS,20尾)、厦门(XM,20尾)、汕头(ST,20尾)、北海(BH,20尾),共计119尾样品(图1)。所有样品取背部新鲜肌肉组织存放于5 mL离心管,加入无水乙醇固定,4 ℃保存备用。
1.2 基因组DNA提取及PCR扩增
剪取适量已浸泡在无水乙醇中的背部肌肉组织,采用传统的酚-氯仿法[26]提取少鳞鱚基因组DNA,用无水乙醇沉淀和75%乙醇润洗后经自然风干的DNA溶解于100 μL灭菌水中,置于4 ℃保存备用。
参考少鳞鱚线粒体基因组全序列(GenBank登录号KR363149),采用Primer primer 6.0软件设计扩增ND2基因的引物,分别为上游引物ND2-1F (5'-CACG AACGCCCCTATACTCA-3')和下游引物ND2-1R (5'-CCTGGTAGGTTGTTAGGGGA-3')。PCR总反应体系为25 μL,包括0.25 μL浓度为5 U·μL−1的Taq酶,各1 μL正反向引物(10 μmol·L−1),2 μL的dNTPs (2.5 mmol·L−1),2.5 μL的10×buffer缓冲液(含Mg2+),1 μL的模板DNA (50~100 ng),以及17.25 μL的灭菌双蒸水。PCR反应程序为95 ℃预变性5 min;95 ℃变性0.5 min,52 ℃退火0.5 min,72 ℃延伸1 min,总计循环35次;72 ℃延伸5 min。用1.5%琼脂糖凝胶电泳检测PCR扩增产物,凝胶成像系统拍照后,挑选条带亮度较高的PCR产物送至上海美吉生物有限公司进行测序。
1.3 数据分析
所获序列均由DNAstar[27]软件包进行编辑、校对和排序;采用DnaSP 6.0[28]和Arlequin[29]软件进行多态位点数目(number of polymorphic sites)、单倍型数目(number of haplotype)、核苷酸多样性(nucleotide diversity)和单倍型多样性(haplotype diversity)等计算;使用MEGA 7.0[30]软件基于Kimura 2-parameter model参数计算少鳞鱚群体间平均遗传距离、碱基组成和碱基转换/颠换,并采用邻接法(neighbor-joining,NJ)基于1 000次Bootstrap重抽样构建单倍型系统发育树;运用Network软件基于中间连接法(median-joining)绘制少鳞鱚单倍型网络图。使用Arlequin[29]软件统计少鳞鱚群体间的遗传分化水平,分子变异水平分析(analysis of molecular variance analysis,AMOVA)种群间遗传变异情况;使用Tajima's D和Fu's Fs中性检验以及核苷酸不配对分布(mismatch distribution)来估算种群的历史动态变化。使用公式τ=2ut[31]估算种群扩张时间,其中τ为扩张时间参数(Tau),u为进化速率,u=μk,其中μ为每个碱基的变异速率,k为所分析序列的长度,t表示自扩张以来所经历的代数,扩张时间T=t×代时,μ参考Bermingham和Mccafferty[32]的研究,采用每百万年2%的突变速率。
2. 结果
2.1 序列变异及遗传多样性组成
使用DNAstar软件包中SeqMan程序进行序列比对,得到长度为450 bp的ND2基因片段。119尾少鳞鱚样本平均碱基组成为A=20.96%,T=30.38%,G=14.05%,C=34.6%,A+T含量(51.34%)略高于G+C (48.66%),其中C含量最高,G含量最低,碱基组成偏向于嘧啶碱基。所获得序列中多态位点数目77个,占所获得序列长度的19.84%,没有发现碱基的插入和缺失。119条序列共检测到61个单倍型,从各群体所含单倍型分布来看,厦门群体检测到的单倍型数量最多,而莱州群体检测到的单倍型数量最少。遗传多样性参数显示,厦门群体的单倍型多样性最高(0.984 2±0.020 5),北海群体单倍型多样性最低(0.910 5±0.053 8);汕头群体核苷酸多样性最高(0.014 225±0.007 953),莱州群体核苷酸多样性最低(0.005 788±0.003 723,表1)。
表 1 基于线粒体ND2基因的少鳞鱚6个种群遗传多样性Table 1. Genetic diversity in seven S. japonica populations based on mtDNA ND2 gene种群
population样本数
number of individuals单倍型数目
number of haplotypes多态位点数目
number of polymorphic sites单倍型多样性
Hd核苷酸多样性
π莱州 LZ 19 12 16 0.924 0±0.045 8 0.005 788±0.003 723 胶南 JN 20 13 17 0.957 9±0.025 5 0.007 800±0.004 739 舟山 ZS 20 13 34 0.947 4±0.032 3 0.013 019±0.007 349 厦门 XM 20 17 21 0.984 2±0.020 5 0.010 282±0.005 993 汕头 ST 20 14 48 0.915 8±0.054 6 0.014 225±0.007 953 北海 BH 20 13 19 0.910 5±0.053 8 0.007 678±0.004 677 总计 total 119 51 77 0.945 3±0.015 5 0.009 718±0.005 445 2.2 群体遗传结构和聚类分析
在MEGA软件中采用Kimura 2-parameter模型计算少鳞鱚6个群体间的遗传距离(表2),种群间的遗传距离在0.006 85~0.010 02,其中厦门和汕头的遗传距离最大 (0.010 02);莱州和北海的遗传距离最小 (0.006 85),遗传分化尚未达到种级水平[33]。
表 2 少鳞鱚不同群体间平均遗传距离Table 2. Average genetic distances among different S. japonica populations种群
population北海
BH胶南
JN莱州
LZ汕头
ST厦门
XM胶南 JN 0.007 74 莱州 LZ 0.006 85 0.006 88 汕头 ST 0.008 64 0.008 51 0.007 66 厦门 XM 0.009 11 0.009 02 0.008 50 0.010 02 舟山 ZS 0.008 04 0.007 95 0.007 42 0.008 98 0.009 17 少鳞鱚不同群体间遗传多样性指数(FST)分析发现两两群体间的FST介于−0.014 16 ~0.044 38,除莱州和厦门、舟山群体之间差异显著外,其余群体间P均大于0.05 (表3)。采用分子方差分析(AMOVA)对6个少鳞鱚群体遗传变异来源和结构进行统计(表4)。结果表明,少鳞鱚遗传变异主要来自群体内个体间(99.96%),其P检验也不显著(P=0.426 20±0.015 26),表明少鳞鱚群体间存在广泛的基因交流,使得遗传分化程度较低、遗传差异不显著。
表 3 基于ND2序列单倍型频率的少鳞鱚群体FST (对角线下) 和相应P (对角线上)Table 3. F-Statistics (below diagonal) and FST P value (above diagonal) from haplotype frequencies of S. japonica种群
population北海
BH胶南
JN莱州
LZ汕头
ST厦门
XM舟山
ZSBH 0.625 00 0.201 17 0.898 44 0.535 16 0.564 45 JN −0.008 44 0.434 57 0.693 36 0.457 03 0.517 58 LZ 0.000 93 0.012 84 0.638 67 0.018 55 0.021 48 ST −0.005 06 −0.011 29 −0.002 92 0.355 47 0.654 30 XM −0.001 54 −0.004 82 0.044 38 0.003 94 0.826 17 ZS −0.002 41 −0.005 85 0.035 13 −0.013 66 −0.014 16 表 4 少鳞鱚群体线粒体ND2基因的AMOVA分析Table 4. AMOVA analysis of S. japonica populations based on mtDNA ND2 gene变异来源
source of variation自由度
df平方和
SS方差分量
variance components变异百分比
percentage of variation种群间 among populations 5 9.720 0.000 78Va 0.04 种群内 within populations 113 217.911 1.928 41Vb 99.96 总计 total 118 227.630 1.929 19 固定值 fixation index (FST) 0.000 41 注:Va. 组群间方差、Vb. 组群内群体间方差
Note: Va. variance between groups; Vb. variance between groups within the group从GenBank数据库下载近缘种多鳞鱚 (S. sihama,GenBank登录号KR363150)和中国鱚 (S. sinica,GenBank登录号KR363151) 作为外群,采用邻接法(neighbor-joining,NJ) 构建系统发育树 (图2),聚类树中不同群体的个体相互交织,没有形成明显的分支,种群之间存在广泛的基因交流。
2.3 群体历史动态及分化时间
在1 000次模拟抽样的情况下,用Arlequin软件对6个群体的少鳞鱚进行中性检验(表5和表6)和核苷酸不配对分析(图3)。当Tajima's D为负,且在统计学上达到显著水平,表明种群分化偏离中性突变理论模型,预示着群体在进化过程中不仅受到随机漂变因素影响,还有可能经历过大规模的群体扩张或受到过选择压力、瓶颈效应的影响[34]。Fu's Fs检验中,如果Fs大于零,表明种群趋于稳定,反之则表明种群趋于扩张[35]。结果表明,无论是从群体还是个体角度来看,Tajima's D和Fu's Fs均为负值且显著偏离中性(P<0.05),表明少鳞鱚在历史进化过程中经历过群体扩张事件。此外,偏离方差(SSD)和Raggedness index (Rg)均较小,统计检验不显著,进一步表明群体扩增假说的成立[29]。核苷酸不配对分布图呈单峰且中性检验差异显著,表明少鳞鱚群体近期可能受到了瓶颈效应的打击。根据核苷酸不配对分布得到的Tau,少鳞鱚性成熟时间为1年,代时取1,计算得到少鳞鱚群体发生扩张的时间大约在 (12~29) 万年前的第四纪更新世(Pleistocene)晚期。
表 5 少鳞鱚群体中性检验及分化时间Table 5. Neutral test and differentiation time of fit for S. japonica中性检验
neutral test莱州
LZ胶南
JN舟山
ZS厦门
XM汕头
ST北海
BHTajima's D −1.918 06 −1.678 23 −0.721 48 −1.244 16 −2.098 64 −1.378 17 P 0.012 00 0.033 00 0.272 00 0.090 00 0.008 00 0.072 00 Fu's Fs −7.368 12 −6.792 29 −3.610 54 −12.466 49 −4.380 38 −6.687 65 P 0.000 00 0.001 00 0.050 00 0.000 00 0.030 00 0.003 00 Tau 2.2 3.7 3.5 4.3 5.3 3.4 T (Mya) 0.12 0.21 0.19 0.24 0.29 0.18 表 6 少鳞鱚种群中性检验Table 6. Neutral test of population for S. japonica数量
numberTajima's D Fu's Fs 偏离方差
SSDRg D P Fs P 119 −2.223 54 0 −25.916 43 0 0.001 89 0.020 24 3. 讨论
单倍型多样性(h)和核苷酸多样性(π)是衡量物种遗传多样性丰富度的重要指标,也是物种生存适应和发展进化的前提。物种的遗传多样性越高,对环境的适应能力越强,对变化也有更强的调整能力[36]。本研究表明中国沿海少鳞鱚群体遗传多样性较高(h=0.945 3,π=0.009 718),与薛泰强等[14]和王林燕[17]分别采用控制区序列和微卫星分子标记对不同地理群体少鳞鱚遗传结构的分析结果相似,都呈现出相对丰富的遗传多样性水平。这可能与ND2基因本身具有较高的突变速率有关,另外种群的高单倍型多样性与种群的数量、所处环境的变化以及生活习性有密切联系[37],少鳞鱚种群繁殖时间较长,单次后代产生数量大,能够弥补因捕食、捕捞等带来的种群缺失,使其存在较高的遗传多样性。
莱州群体核苷酸多样性指数最低,且与厦门、舟山群体间的遗传差异显著。由于莱州位于采样点最北端,受渤海沿岸流和黄海冷水团的共同作用,所处海域水温相对较低,而浙闽沿岸受亚热带海洋性季风气候影响,水温普遍偏高,推测温度变化可能对少鳞鱚遗传变异产生一定影响。鱼类作为变温性动物,温度变化对其分布和生长发育都具有重要影响。温度变化限制了鱼类成体的扩散并促进其适应局限的生存环境,使得彼此间交流的概率大幅降低,这种差异经长期积累并最终导致遗传分化。有文献资料表明,一些海洋鱼类的遗传多样性分布受到温度和纬度变化的影响,如王秀亮[38]对西北太平洋玉筋鱼(Ammodytes personatus) Cyt b基因研究发现,玉筋鱼单倍型频率在地理分布上与水温变化具有相关性,因而推测水域环境温度的改变以及洋流系统对玉筋鱼种群扩散及种群地理分布格局可能有着重要影响;Xu等[39]对中日沿海5个群体褐菖鲉(Sebastiscus marmoratus) ATP6和Cyt b基因所编码的氨基酸替换率进行比较时发现,温度可能是影响群体遗传结构和遗传分化的潜在环境因素。
遗传多样性指数FST是评价物种群体间遗传分化尺度的标准。Wright[40]的研究发现,若FST>0.25,表明群体间存在高度分化,0.15<FST<0.25,表示群体间存在中度分化,0.05<FST<0.15表示群体间存在低度分化,而FST<0.05则表明群体间没有分化。本研究发现少鳞鱚各群体间FST均小于0.05,结合平均遗传距离的大小来看,暗示群体间未发生明显的遗传分化,这与Gao等[16]的研究结果一致。推测少鳞鱚这类小型鱼类群体数量大、分布范围广,加之早期营浮游生活时间长、本身具备较强迁移能力,基因交流频繁造成遗传的同质化,使得种群间遗传分化水平普遍偏低,阻碍了种群的分化。此外,基于邻接法构建的系统关系树显示来自不同群体的个体彼此相互交错分布,没有呈现显著的谱系结构,也反映了少鳞鱚遗传分化水平较低这一结论。
Hewitt[41]通过研究证实,更新世的剧烈气候波动对地球上动植物遗传多样性和分布格局产生了重要影响。基于线粒体控制区序列估算出中日韩不同地理种群少鳞鱚发生群体扩张的时间大约在 (0.423~1.06)百万年前的晚更新世时期[16],本研究对ND2基因的核苷酸不配对分析得出少鳞鱚群体分化时间大约在 (0.12~0.29)百万年前的第四纪冰川期-间冰期旋回[42]。尽管线粒体不同序列变异速率存在一定差异,但推算出扩张时间均处于第四纪更新世晚期,这一时期气候的急剧变化,全球性大幅度气温变冷,中高纬形成大面积冰盖,大气环流和洋流的变化直接影响动植物生长、演化和分布,导致大量生物种群的迁徙或者灭绝[43],使得少鳞鱚等海洋鱼类的空间分布格局以及遗传结构发生了较大变化。
鱼类种群结构的时空分布和动态变化与其所处水域的环境因子有着密切关系,过度捕捞、水体污染、工程建设以及外来物种入侵等都会使鱼类资源和遗传多样性遭到破坏。基于线粒体ND2基因分析的结果显示,我国沿海少鳞鱚种群遗传多样性丰富程度较高,种质资源尚处于相对安全状态,具有一定的资源开发与利用潜力,但由于多遗传标记结合对于揭示种群进化历史和谱系地理格局的形成往往更全面、更准确。因此,在今后的渔业资源管理过程中,还需要综合采用多种分子生物学标记技术来监测少鳞鱚的遗传多样性水平,以实现其资源多样性的合理保护和可持续利用。
-
表 1 PCR引物
Table 1 Primers used for PCR
目的基因
Target gene引物序列 (5'—3')
Primer sequence (5'—3')产物长度
Product size/bp参考文献
Reference16S rRNA 27F: AGAGTTTGATCCTGGCTCAG
1492R: GGTTACCTTGTTACGACTT约1 500 [8] rpoD F:ATAGAAATAACCAGACGTAAGTTNGCYTCNACCATYTCYTTYT
R:ACGACTGACCCGGTACGCATGTAYATGMGNGARATGGGNACNGT780 [16] 氨基糖苷类耐药基因
Aminoglycoside resistance geneaph (3')-Ⅰa F: TGACTGGGCACAACAGACAA
R: CGGCGATACCGTAAAGCAC677 [17] aac (6') -Ⅰb F: ATGACCTTGCCATGCTCTATGA
R: CGAATGCCTGGCGTGTTT486 [17] aac (3')-Ⅰb F: ACCCTACGAGGAGACTCTGAATG
R: CCAAGCATCGGCATCTCATA384 [17] 酰胺醇类耐药基因
Amide alcohol resistance genefloR F: GCTTCACTGGCGATGGATATTTA
R: CAAAGTAATGAATATCGCCTGCC450 [18] catA F: AAAAATTATATCCGACTCTCTTA
R: CTTGAATCGATAAGGGAATATAG420 [18] cmlA F: TGCCAGCAGTGCCGTTTAT
R: CACCGCCCAAGCAGAAGTA500 [18] 四环素类耐药基因
Tetracycline resistance genetetA F: TTTCGGGTTCGGGATGGT
R: CAGGCAGAGCAAGTAGAGGG915 [19] tetC F: CTGGGCTGCTTCCTAATGC
R: AGCTGTCCCTGATGGTCGT480 [19] tetM F: GAGGTCCGTCTGAACTTTGCG
R: AGAAAGGATTTGGCGGCACT580 [19] 磺胺类耐药基因
Sulfonamide resistance genesul1 F: CATTGCCTGGTTGCTTCAT
R: ATCCGACTCGCAGCATTT238 [19] sul2 F: CATCATTTTCGGCATCGTC
R: TCTTGCGGTTTCTTTCAGC793 [19] sul3 F: AGATGTGATTGATTTGGGAGC
R: TAGTTGTTTCTGGATTAGAGCCT443 [19] 喹诺酮类耐药基因
Quinolone resistance geneqnrA F: ATTTCTCACGCCAGGATTTG
R: GATCGGCAAAGGTCAGGTCA516 [20] qnrB F: GATCGTGAAAGCCAGAAAGG
R: ACGATGCCTGGTAGTTGTCC469 [20] qnrS F: ACGACATTCGTCAACTGCAA
R: TAAATTGGCACCCTGTAGGC417 [20] aac (6')-Ⅰb-cr F: TTGCGATGCTCTATGAGTGGCTA
R: CTCGAATGCCTGGCGTGTTT523 [20] oqxAB F: CCCTGGACCGCACATAAAG
R: AAAGAACAAGATTCACCGCAAC482 [20] 表 2 基于16S rRNA测序的394株分离株鉴定结果
Table 2 Identification result of 394 strains based on 16S rRNA sequencing
属
Genus种
Species6月
Jun.7月
Jul.8月
Aug.9月
Sep.10月
Oct.合计
Total弧菌属 Vibrio 大菱鲆弧菌 V. scophthalmi 13 10 35 26 33 117 (29.70%) 溶藻弧菌 V. alginolyticus 0 1 23 11 0 35 (8.89%) V. atlanticus 2 3 1 6 3 15 (3.81%) 嗜环弧菌 V. cyclitrophicus 0 0 0 4 6 10 (2.54%) V. toranzoniae 1 2 1 7 4 15 (3.81%) 其他 Other 7 3 2 17 11 40 (10.15%) 假单胞菌属 Pseudomonas 8 3 8 11 16 46 (11.68%) 假交替单胞菌属 Pseudoalteromonas 0 5 12 18 0 35 (8.88%) 其他 Other 14 17 10 33 7 81 (20.56%) 合计 Total 45 44 92 133 80 394 表 3 18株大菱鲆弧菌的药物敏感性
Table 3 Drug sensitivity of 18 V. scophthalmi strains
分离株
Isolate恩诺沙星
ENR硫酸新霉素
NEO甲砜霉素
TAP氟苯尼考
FFC盐酸多西环素
DO氟甲喹
FLU磺胺间甲氧嘧啶钠
SMM磺胺甲恶唑+甲氧苄啶
SMZ+TMPVs1 S I S S S S S S Vs2 S S R R S S R S Vs3 S S S S S S R R Vs4 S S S S S S S S Vs5 S S S S S S R S Vs6 S I S S S S S S Vs7 S S S S S S R R Vs8 S S S S S S S S Vs9 R I R R S S S S Vs10 R S R R S S S S Vs11 S I R R S S S S Vs12 R S R R S S S S Vs13 S I S S S S R R Vs14 S I S S S S S S Vs15 I S R R S S R R Vs16 R S R R S S S S Vs17 R I R R S S S S Vs18 R I R R S S S S 注:S. 敏感;I. 中敏;R. 耐药。 Note: S. Sensitive; I. Intermediary; R. Resistant. 表 4 大菱鲆弧菌耐药基因携带情况
Table 4 Drug resistance genes status of V. scophthalmi
抗生素类型
Type of antibiotics耐药基因
Drug resistance gene携带耐药基因的分离株 (N=18)
Isolates carrying resistance gene氨基糖苷类 Aminoglycosides aac (6')-Ⅰb, aac (6')-Ⅰb, aac (3')-Ⅰb 0 (0%) 四环素类 Tetracyclines tetA, tetC, tetM 0 (0%) 酰胺醇类 Amido alcohols floR Vs2, Vs3, Vs4, Vs7, Vs9, Vs11, Vs12, Vs15, Vs16, Vs17, Vs18 (61.11%) catA 0 (0%) cmlA Vs1, Vs3, Vs4, Vs5, Vs6, Vs7, Vs8, Vs9, Vs12, Vs15, Vs16, Vs18 (66.67%) 磺胺类 Sulfonamides sul1 0 (0%) sul2 Vs3, Vs4, Vs7, Vs9, Vs11, Vs12, Vs15, Vs16, Vs17, Vs18 (55.56%) sul3 0 (0%) 喹诺酮类 Quinolones qnrA Vs1, Vs2, Vs3, Vs5, Vs7, Vs8, Vs9, Vs14, Vs18 (50%) qnrB 0 (0%) qnrS Vs11 (5.56%) oqxAB 0 (0%) aac (6')-Ⅰb-cr 0 (0%) 表 5 大菱鲆弧菌分离株耐药基因与耐药表型符合率
Table 5 Coincidence rate of drug resistance genes and phenotype of V. scophthalmi strains
抗生素类型
Type of antibiotics耐药表型株数
Number of isolates with
drug resistance phenotype有耐药表型且携带耐药基因株数
Number of strains with drug
resistance phenotype and
drug resistance genes符合率
Coincidence rate/%氨基糖苷类 Aminoglycosides 0 0 100 四环素类 Tetracyclines 0 0 100 酰胺醇类 Amido alcohols 9 8 88.89 磺胺类 Sulfonamides 6 3 50 喹诺酮类 Quinolones 6 2 33.33 表 6 大菱鲆弧菌分离株耐药谱、耐药基因谱和ERIC基因型
Table 6 Drug resistance spectrum, resistance gene spectrum and ERIC genotype of V. scophthalmi strains
耐药谱
Antibiotic resistance spectrum耐药基因谱
Drug resistance genes spectrumERIC基因型
ERIC genotypeVs1 无 cmlA+qnrA Ⅰ Vs2 甲砜霉素、氟苯尼考、磺胺间甲氧嘧啶钠 floR+qnrA Ⅰ Vs6 无 cmlA Ⅰ Vs7 磺胺间甲氧嘧啶钠、磺胺甲恶唑甲氧苄啶复合物 floR+cmlA+sul2+qnrA Ⅰ Vs10 甲砜霉素、氟苯尼考、恩诺沙星 无 Ⅰ Vs13 磺胺间甲氧嘧啶钠、磺胺甲恶唑甲氧苄啶复合物 无 Ⅰ Vs15 甲砜霉素、氟苯尼考、磺胺间甲氧嘧啶钠、磺胺甲恶唑甲氧苄啶复合物 floR+cmlA+sul2 Ⅰ Vs18 甲砜霉素、氟苯尼考、恩诺沙星 floR+cmlA+sul2+qnrA Ⅰ Vs11 甲砜霉素、氟苯尼考 floR+cmlA+sul2+qnrS Ⅱ Vs3 磺胺间甲氧嘧啶钠、磺胺甲恶唑甲氧苄啶复合物 floR+cmlA+sul2+qnrA Ⅲ Vs4 无 floR+cmlA+sul2 Ⅲ Vs5 磺胺间甲氧嘧啶钠 cmlA+qnrA Ⅲ Vs8 无 cmlA+qnrA Ⅲ Vs9 甲砜霉素、氟苯尼考、恩诺沙星 floR+cmlA+sul2+qnrA Ⅲ Vs14 无 qnrA Ⅲ Vs16 甲砜霉素、氟苯尼考、恩诺沙星 floR+cmlA+sul2 Ⅲ Vs17 甲砜霉素、氟苯尼考、恩诺沙星 floR+sul2 Ⅲ Vs12 甲砜霉素、氟苯尼考、恩诺沙星 floR+cmlA+sul2 Ⅳ -
[1] 崔惠敬, 孟玉霞, 徐永平, 等. 植物精油对大菱鲆弧菌的体外和体内抗菌活性[J]. 微生物学通报, 2017, 44(2): 274-284. [2] 关长涛, 杨正勇, 王启要, 等. 大菱鲆产业发展报告[J]. 中国水产, 2021(1): 22-35. [3] 于新然, 姚洪, 叶仕根, 等. 养殖大菱鲆感染迟缓爱德华氏菌的分离、毒力基因及ERIC-PCR分析[J]. 大连海洋大学学报, 2018, 33(2): 169-174. [4] 崔惠敬, 孟玉霞, 冯文倩, 等. 养殖大菱鲆肠道中大菱鲆弧菌的分离鉴定及药敏试验[J]. 水产科学, 2017, 36(2): 125-131. [5] ZHANG Z, YU Y X, WANG Y G, et al. Complete genome analysis of a virulent Vibrio scophthalmi strain VSc190401 isolated from diseased marine fish half-smooth tongue sole, Cynoglossus semilaevis[J]. BMC Microbiol, 2020, 20(1): 1-2. doi: 10.1186/s12866-019-1672-7
[6] WANG X A, MA A J, SUN Z B. Genetic parameters of seven immune factors in turbot (Scophthalmus maximus) infected with Vibrio anguillarum[J]. J Fish Dis, 2021, 44(3): 263-271. doi: 10.1111/jfd.13320
[7] 刘朝阳, 王印庚, 孙晓庆. 颗粒饲料携带细菌与大菱鲆疾病发生的相关性[J]. 南方水产, 2009, 5(4): 13-21. [8] 王岚, 王印庚, 张正, 等. 养殖大菱鲆 (Scophthalmus maximus) 腹水病的病原多样性及其耐药性分析[J]. 渔业科学进展, 2017, 38(4): 17-24. doi: 10.11758/yykxjz.20160417001 [9] QIAO G, LEE D C, WOO S H, et al. Microbiological characteristics of Vibrio scophthalmi isolates from diseased olive flounder Paralichthys olivaceus[J]. Fish Sci, 2012, 78(4): 853-863. doi: 10.1007/s12562-012-0502-8
[10] 郑怀东, 李小进, 黄萍萍, 等. 大菱鲆弧菌耐药性初探[J]. 中国水产, 2017(9): 65-68. [11] 王鹤, 王力勇, 刘朋, 等. 养殖大菱鲆病原菌耐药现状与用药建议[J]. 中国水产, 2019(8): 32-36. [12] FU H Y, YU P, LIANG W L, et al. Virulence, resistance, and genomic fingerprint traits of Vibrio cholerae isolated from 12 species of aquatic products in Shanghai, China[J]. Microb Drug Resist, 2020, 26(12): 1526-1539. doi: 10.1089/mdr.2020.0269
[13] ZHANG X J, YANG W M, ZHANG D F, et al. Does the gastrointestinal tract serve as the infectious route of Aeromonas hydrophila in crucian carp (Carassius carassius)?[J]. Aquacult Res, 2015, 46(1): 141-154. doi: 10.1111/are.12168
[14] 张海强, 安伟, 肖雨. 虾源副溶血弧菌致病基因检测与ERIC-PCR分型[J]. 上海海洋大学学报, 2019, 28(6): 848-856. doi: 10.12024/jsou.20181102433 [15] 张永政, 李聪, 王世锋, 等. 海南罗非鱼源维氏气单胞菌的分离、毒力基因检测及ERIC-PCR分型[J]. 海南大学学报(自然科学版), 2017, 35(3): 273-281. [16] TARAZONA E, PÉREZ-CATALUÑA A, LUCENA T, et al. Multilocus sequence analysis of the redefined clade Scophthalmi in the genus Vibrio[J]. Syst Appl Microbiol, 2015, 38(3): 169-175. doi: 10.1016/j.syapm.2015.03.005
[17] 杨明伟, 韦慕兰, 罗福广, 等. 罗非鱼源无乳链球菌对氨基糖苷类耐药性及耐药基因检测[J]. 西南农业学报, 2018, 31(11): 2438-2444. [18] 赵敏. 四川地区斑点叉尾鮰源维氏气单胞菌的耐药性以及耐药基因检测分析[D]. 成都: 四川农业大学, 2013: 23-24. [19] 乔毅. 江苏省沿海地区水产养殖主要致病菌耐药性研究 [D].上海: 上海海洋大学, 2015: 57-58. [20] 魏文娟, 赵姝, 王元, 等. 副溶血弧菌养殖对虾分离株耐药性及耐药基因分析[J]. 南方水产科学, 2020, 16(1): 9-16. [21] CLSI. Methods for antimicrobial dilution and disk susceptibility testing of infrequently isolated or fastidious bacteria[M]. 2nd ed. Wayne, PA: Clinical and Laboratory Standards Institute, 2010: 70-96. [22] 杨晓霞, 阙国鹰. 16S核糖体RNA基因测序技术在口腔微生物多样性中的研究进展[J]. 中南大学学报(医学版), 2020, 45(7): 849-855. doi: 10.11817/j.issn.1672-7347.2020.190236 [23] 黄华, 张建柏, 徐晓莹, 等. 养殖大菱鲆细菌性腹水病及其防治方法研究进展[J]. 水产科技情报, 2021, 48(2): 112-116. [24] PASCUAL J, MACIAN M C, ARAHAL D R, et al. Multilocus sequence analysis of the central clade of the genus Vibrio by using the 16S rRNA, recA, pyrH, rpoD, gyrB, rctB and toxR genes[J]. Int J Syst Evol Micr, 2010, 60(Pt 1): 154-165.
[25] MO W Y, CHEN Z, LEUNG H M, et al. Application of veterinary antibiotics in China's aquaculture industry and their potential human health risks[J]. Environ Sci Pollut R, 2017, 24(10): 8978-8989. doi: 10.1007/s11356-015-5607-z
[26] XU N, LI M, AI X, et al. Determination of pharmacokinetic and pharmacokinetic-pharmacodynamic parameters of doxycycline against Edwardsiella ictaluri in yellow catfish (Pelteobagrus fulvidraco)[J]. Antibiot, 2021, 10(3): 1-2.
[27] 孙朋辉, 吕利群. 硫酸新霉素控制肠型嗜水气单胞菌性草鱼肠炎病的给药方案[J]. 水产科学, 2016, 35(1): 54-59. [28] 张刚, 冯婕. 细菌固有耐药的研究进展[J]. 遗传, 2016, 38(10): 872-880. [29] 刘蔚雯, 汪洋. 细菌对氟苯尼考的耐药机制研究进展[J]. 中国动物传染病学报, 2018, 26(1): 1-6. [30] 孙仲琳, 崔锦, 王锐, 等. 细菌耐药的外排机制及其抑制剂[J]. 复旦学报(自然科学版), 2019, 58(2): 135-143. [31] 高晓华, 安伟, 张明辉. 上海凡纳滨对虾源副溶血弧菌耐药性及其耐药基因检测分析[J]. 南方农业学报, 2021, 52(3): 827-836. doi: 10.3969/j.issn.2095-1191.2021.03.031 [32] 吴甘林, 邓玉婷, 姜兰, 等. 水产动物源细菌质粒介导的喹诺酮类耐药研究概况[J]. 水产科学, 2020, 39(4): 631-638. [33] MARSHALL B M, LEVY S B. Food animals and antimicrobials: impacts on human health[J]. Clin Microbiol Rev, 2011, 24(4): 718-733. doi: 10.1128/CMR.00002-11
[34] NIU G, WONGSATHEIN D, BOONYAYATRA S, et al. Occurrence of multiple antibiotic resistance and genotypic characterization in Edwardsiella tarda isolated from cage-cultured hybrid red tilapia (Oreochromis sp.) in the Ping River, Northern Thailand[J]. Aquacult Res, 2019, 50(12): 3643-3652. doi: 10.1111/are.14322
[35] BARTOSZEWICZ M, CZYZEWSKA U. Spores and vegetative cells of phenotypically and genetically diverse Bacillus cereus sensu lato are common bacteria in fresh water of northeastern Poland[J]. Can J Microbiol, 2017, 63(12): 939-950. doi: 10.1139/cjm-2017-0337
-
期刊类型引用(3)
1. 张文博,黄洪辉,巩秀玉,刘华雪. 南海西部不同粒级浮游动物碳氮稳定同位素研究. 海洋环境科学. 2025(01): 89-96 . 百度学术
2. 贡艺,梁茜,李云凯,刘必林,牛成功. 西北太平洋桡足类稳定同位素特征及其影响因素. 生态学杂志. 2024(02): 505-513 . 百度学术
3. 李纯厚,齐占会. 中国渔业生态环境学科研究进展与展望. 水产学报. 2023(11): 132-147 . 百度学术
其他类型引用(2)