Density variation of Vibrio in ablone aquaculture water and its relationship with environmental factors
-
摘要:
采用培养法对皱纹盘鲍(Haliotis discus hannai)养殖水体中的弧菌(Vibrio)密度进行了测定, 监测了主要环境因子(水温、溶解氧、pH、盐度、氧化还原电位、亚硝酸盐、硝酸盐、氨盐、磷酸盐), 并探讨了弧菌密度与环境因子的关系。结果表明, 养殖水中弧菌密度全年变化范围为0.2×102~6.0×103 cfu·mL-1, 平均为1.7×103 cfu·mL-1。水温≤25 ℃时弧菌密度与各环境因子均无相关性。水温＞25 ℃时弧菌密度与水温显著相关(P＜0.05), 与氧化还原电位、硝酸盐、氨盐极显著相关(P＜0.01), 与其他环境因子相关性不明显。
Abstract:We measured the Vibrio density in aquaculture water of Haliotis discus hannai by cultivation method, and investigated the main environmental factors including temperature, dissolved oxygen, potential of hydrogen, salinity, redox potential, nitrite, nitrate, ammonium and phosphate, so as to study the relationship between Vibrio quantity and main environmental factors. The results show that Vibrio density throughout the whole year varied from 0.2×102 cfu·mL-1 to 6.03×103 cfu·mL-1(average 1.7×103 cfu·mL-1). When the temperature was less than or equal to 25 ℃, the Vibrio density had no correlation with any environmental factor; when the temperature was greater than 25 ℃, the Vibrio density had significant correlation with temperature (P < 0.05), redox potential (P < 0.01), nitrate (P < 0.01) and ammonium (P < 0.01), but had no correlation with the others.
-
Keywords:
- Haliotis discus hannai /
- Vibrio /
- environmental factors
-
细胞色素P450(cytochrome P450，CYPs)是一类广泛分布于生物体内的超基因家族，其酶系能够催化内源物质的合成、降解以及外源性物质的代谢[1]，在解毒代谢阶段I起到重要作用[2]。CYPs基因具有结构多样性，这决定了细胞色素P450酶系可能是自然界最具催化多样性的酶系之一[3]。NELSON等[2]基于481个细胞色素P450基因和22个假基因的趋异进化关系以及氨基酸序列相似度将CYPs归为不同的“家族(family)”、“亚家族(subfamily)”，同一家族序列同源性大于40%，同一亚家族序列同源性大于55%[2, 4]。目前无脊椎动物中已发现的CYPs家族数已超过了70多个[5]，而软体动物中CYPs基因的相关研究基础十分薄弱，仅在栉孔扇贝(Chlamys farreri)[6]、菲律宾蛤仔(Ruditapes philippinarum)[7]、翡翠贻贝(Perna viridis)[8]以及长牡蛎(Crassostrea gigas)[9]等少数物种有相关基因的研究报道，且缺乏系统的命名。
福寿螺(Pomacea canaliculata)具有适应性强、繁殖力高等特点，严重危害水域生态、农作物及人类健康，是中国恶性外来入侵种之一[10]。目前普遍采用四聚乙醛(methaldehyde)为有效成分的药物对其进行杀灭防治，然而四聚乙醛的持续重复施用诱导了福寿螺对其产生微弱抗药性[11-12]。细胞色素P450作为生物体内三大解毒酶系之一，已经被证明是害虫对杀虫剂产生抗性和交互抗性的主要原因之一[13-14]，研究认为海洋贝类CYPs酶活力与水体有毒有机物代谢密切相关[7]，而福寿螺耐药机理还不明了，相关研究十分薄弱。基于此，文章克隆分析了福寿螺CYP3192A1，采用荧光定量PCR(real time quantitative PCR，RT-qPCR)检测了CYP3192A1基因表达水平的组织及性别差异性，进而研究浓度及时间维度上四聚乙醛处理与CYP3192A1表达水平的相关性，初步探索CYP3192A1在四聚乙醛代谢中的作用地位，以为研究福寿螺CYPs基因结构与解毒功能，细胞色素P450功能及底物特异性奠定理论基础。
1. 材料与方法
1.1 材料处理
福寿螺采集于中国水产科学研究院珠江水产研究所外来水生生物入侵风险评估中心，螺体质量为9.32~15.28 g，壳高3.5~4.5 cm，取自相同世代群体。分别取雌、雄螺的肝、心、中肠、鳃、胃及腹足等6个组织，每个组织设3个样重复，采用RT-qPCR测定CYP3192A1基因表达水平，每样3个重复。
分别以0.6 mg·L-1、2.0 mg·L-1的亚致死浓度四聚乙醛处理福寿螺雌性群体，在用药后第0、第2、第6、第12、第24小时采集样品，选取鳃、肝、肠、胃4个组织，每个组织设3个样重复，RT-qPCR测定CYP3192A1表达水平，每样3个重复。
1.2 总RNA的提取
应用Mollusc RNA Kit R6875-01(OMEGA Bio-TEK)试剂盒提取福寿螺组织中总RNA，以电泳仪(BG-submidi)和核酸蛋白分析仪(东南D30)检测RNA质量与浓度后，置于-80 ℃保存备用。
1.3 细胞色素P450核心区cDNA克隆
应用PrimeScriptTMⅡ1st Strand cDNA synthesis Kit(TaKaRa)试剂盒进行反转录获得cDNA，-20 ℃保存。Primer Premier 5.0软件设计目的片段扩增引物P3-F和P3-R(表 1)，以TaKaRa LA Taq®(RR02MA)进行PCR，测序；纯化目的片段、连接至pMD19-T载体，转染感受态细胞E.coli DH5α，菌液涂布于SOC培养基上12 h后挑菌落，经克隆、测序、鉴定后与已知序列进行比对。
1.4 cDNA末端快速扩增
以核心区序列为模板设计5′RACE-PCR套式引物P3-5′-1、P3-5′-2(表 1)，应用Clontech SMAR-Ter® RACE cDNA Amplification Kit试剂盒进行反转录和Outer-PCR，通用引物为UPM；利用LA Taq®和Ex Taq®(TaKaRa)进行Inner-PCR，通用引物为NUP，依次进行PCR测序和克隆测序。设计引物P3-3′-1、P3-3′-2、P3-3′-3、P3-3′-4进行3′ RACE-PCR，利用试剂盒3′-Full RACE core Set with Prime ScriptTM RTase、LA Taq®(TaKaRa)、HiFi酶(Tran-sgene bio-tech)获得3′末端序列，通用引物为3′ Outer Primer、3′ Inner Primer，克隆测序。PCR主程序为94 ℃，30 s；退火30 s；72 ℃，30~60 s；30个循环；退火温度值见表 1。
表 1 细胞色素CYP3192A1扩增、RACE-PCR引物及荧光定量PCR引物Table 1. Primers used for PCR, RACE-PCR and quatitative real-time PCR引物名称
primer引物序列(5′→3′)
primer sequence扩增片段/bp
amplified fragment退火温度/℃
annealing temperatureP3-F
P3-R
P3-5′-1
P3-5′-2
P3-3′-1
P3-3′-2
P3-3′-3
P3-3′-4
UPM
NUP
3′ Outer Primer
3′ Inner Primer
18S-f
18S-r
P3-QT-F
P3-QT-R5′-ATCATCCACACTCTCACTGC-3′
5′-GAATCGTTACGCCGTTGA-3′
5′-ATCCTCTGCCTTGGCGTTTTGTTCC-3′
5′-TAACCGAATGACCTTCCACACAGCC-3′
5′-ATCCACCCGACAGACTCTCACAGGC-3′
5′-CTCAGCAAGAGGTTGTCTCCCAGTG-3′
5′-GCTACCCACGACCCCACCAACTACC-3′
5′-CTTCCCTTCGGTCTGGGTCCTCGTC-3′
5′-CTAATACGACTCACTATAGGGCAAGCAGTGGTATCAACGCAGAGT-3′
5′-CTAATACGACTCACTATAGGGC-3′
5′-AAGCAGTGGTATCAACGCAGAGT-3′Control Reagents
5′-TACCGTCGTTCCACTAGTGATTT-3′
5′-CGCGGATCCTCCACTAGTGATTTCACTATAGG-3′
5′-AATACATGCAAACCAGCTCC-3′
5′-ATTTTTCGTCACTACCTCCC-3′
5′-TTCGGTCTGGGTCCTCGTCA-3′
5′-TACAGCCCTGGTGGAAGTGG-3′1 247
1 247
758
672
1 225
1 094
1 010
922
5′ universal Primer mix
nested universal primer
universal
universal
329
329
322
32256
56
60
62
60
60
62
62
57
57
57
571.5 实时荧光定量PCR
预实验基础上以18S rDNA为荧光定量表达内参基因。参考福寿螺18S(EU520453.1)设计引物18S-f、18S-r；设计福寿螺CYP3A序列克隆片段引物P3-QT-F、P3-QT-R(表 1)，分别克隆、连接载体、转染感受态细胞制备18S rDNA和CYP3192A1质粒。分别将质粒以10-1数量级稀释，共稀释6个浓度梯度；分别与SYBR® Select Master Mix体系混合，置于LightCycler480 Software Setup(Roche罗氏)进行荧光定量PCR，获得荧光信号强度与浓度的标准曲线，结果显示18S内参基因表达的标准曲线为y=-3.404x+38，R2=0.993 3，E=98.4%；CYP3192A1基因的为y=-3.404 262x+34.09，R2=0.990 7，E=101.3%。说明实验数据可靠，可继续进行实验数据的检测。同实验组每个组织3个样品重复，然后每个cDNA样品分别进行荧光定量PCR，测定每个样CYP3192A1绝对定量表达，每个样品3次检测重复。
1.6 序列分析与数据统计
序列拼接与组装，获得CYP3192A1基因cDNA全序列；ORF finder在线工具预测开放阅读框(open reading frame，ORF)序列；NCBI数据库对蛋白序列进行Blastp和SmartBLAST搜索相似序列，MEGA 5.0进行ClustalW比对并构建Neighbor-Joining树；用ProtParam在线工具预测蛋白质的理论等电点pI、分子质量和亲水性等氨基酸特性；Motif Scan工具预测序列的模体结构；TMHMM server.V.2.0预测跨膜结构区；Interpro预测蛋白结构域并归类；SWISS-MODEL Workspace预测蛋白三维结构并分析底物识别位点(substrate recognition sites，SRS)；相似序列比对分析CYP3192A1序列保守位点。
用2-ΔΔCT法处理分析荧光定量PCR数据[15]，采用软件Excel 2010作图，SPSS 19.0统计分析性别、组织表达差异性及药物处理与福寿螺CYP3192A1表达水平的相关性。
2. 结果
2.1 CYP3192A1的克隆与鉴定分析
通过PCR和RACE-PCR从福寿螺鳃组织中克隆获得一条长度为2 523 bp的CYPs基因cDNA全序列(GenBank号KU923379)；包含了260 bp的5′非编码区(UTR)和712 bp的3′UTR序列，以及1 551 bp的开放阅读框序列，预测显示其编码517个氨基酸(图 1)，将氨基酸序列提交至细胞色素P450命名委员会，定名为CYP3192A1。预测蛋白质分子量约为58.0 kD，理论等电点为8.61，亲水性较弱，亮氨酸(Leu)含量高达12.6%。
![]()
图 1 福寿螺CYP3192A1克隆及氨基酸序列分析CYP3192A1的核酸及编码氨基酸序列，不同背景颜色氨基酸序列为P450超家族E class的标签序列；蓝色字体标注的序列为跨膜螺旋；红色字体的序列依次为C螺旋(W***R)、I螺旋(AG*ET)、K螺旋(E**R)和血红素结合区F**G*** C*G。Figure 1. Cloning and amino acid sequence analysis of P.canaliculata CYP3192A1The CYP3A-like sequences of nucleic acids and amino acids encoding; the sequence amino acids were marked by different background colors for P450 superfamily E class groupⅠsignature. Blue font tags sequence represent transmembrane helix. Red fonts represent helix C (W***R), I helix (AG*ET), helix K (E**R) and heme combination zone F**G***C*G.跨膜结构域预测显示1~11号氨基酸位于生物膜以内，35~517氨基酸位于细胞质内，而12~34氨基酸为一个跨膜螺旋(图 1)，其中12~31氨基酸镶嵌于生物膜内部；蛋白质结构分析显示其归类于细胞色素P450酶，具有P450超家族E类(E class)的标签(signature)序列(PR00463)：69~88氨基酸(amino acid，AA)、93~114 AA、380~398 AA、420~444 AA、452~462 AA以及462~485 AA(图 1)；还具有CYP450高度保守序列血红素结合区F** G***C*G、K螺旋(E**R)、I螺旋(AG*ET)、C螺旋(W***R)(图 1)等，其中亚铁血红素结合区视为CYPs的特征基序，与I螺旋、K螺旋等参与稳定蛋白核心结构[16]，空间结构如图 2。
![]()
图 2 参照CYP3A4构建的福寿螺CYP3192A1蛋白质空间结构图Figure 2. Three-dimensional structure model of P.canaliculata CYP3192A1 protein basing on CYP3A42.2 CYP3192A1的序列比对与系统发育分析
以福寿螺CYP3192A1氨基酸进行Blastp同源比对搜索，获得不同物种同源性序列，均为CYP3A亚家族序列；基于不同物种的CYP3A亚家族氨基酸序列同源性构建了系统发育树(图 3)，聚类分析显示福寿螺CYP3192A1与长牡蛎的CYP3A24归为一支，鱼类的CYP3A亚家族序列归为一支，人(Homo sapiens)和印度野牛(Bos taurus)等哺乳动物CYP3A亚家族归为一支，各种鸟的CYP3A亚家族序列归为一支，物种亲缘关系的亲疏符合传统的动物分类学中分类法则。与福寿螺CYP3192A1基因氨基酸序列同源性最高的为长牡蛎的CYP3A24序列，同源性仅为34%；推测福寿螺该CYPs属于细胞色素P450新家族的一个成员；将其序列提交至细胞色素P450命名委员会，正式定名为CYP3192A1。
![]()
表 2 图 3中基因、物种名及对应的登录号Table 2. Gene, Latin name, Chinese name of species and their Accession No. in Fig. 3基因gene 物种species 登录号Accession no. CY3A9 白鹭Egretta garzetta KFP21783.1 CYP3A21 红蜂虎Merops nubicus KFQ25408.1 CYP3A24 红冠蕉鹃Tauraco erythrolophus KFV20621.1 CYP3A24 鹃鴗Leptosomus discolor KFQ05822.1 CYP3A9 马来犀鸟Buceros rhinoceros silvestris KFO86946.1 CYP3A9 金领娇鹟Manacus vitellinus KFW84447.1 CYP3A9 短嘴鸦Corvus brachyrhynchos KFO59023.1 CYP3A9 刺鹩Acanthisitta chloris KFP82045.1 CYP3A24 双领鸻Charadrius vociferous KGL97267.1 CYP3A37 鸵鸟Struthio camelus BAN33784.1 CYP3A24 南非鸵鸟Struthio camelus australis KFV75188.1 CYP3A13 小家鼠Mus musculus NP 031845.1 CYP3A43 智人Homo sapiens NP 476436.1 CYP3A4 普通牛Bos taurus NP 001092837.1 CYP3A56 底鳉Fundulus heteroclitus NP 001296866.1 CYP3A27 安大略鲑Salmo salar ACI33861.1 CYP3A45 虹鳟鱼Oncorhynchus mykiss AAK58569.1 CYP3A24 长牡蛎Crassostrea gigas EKC24008.1 2.3 CYP3192A1基因表达的性别及组织差异性
福寿螺CYP3192A1表达具有性别差异性，性别对福寿螺肝、肠、鳃和心组织中CYP3192A1表达水平具有极显著影响(P < 0.01)，而对胃及腹足中CYP3192A1表达水平无显著影响(P＞0.05)。肠组织CYP3192A1的表达量表现为雄性高于雌性，而在其他组织中CYP3192A1的表达水平均为雄性低于雌性(图 4)。
![]()
图 4 福寿螺CYP3192A1基因在雌雄各组织中的表达情况A、B表示福寿螺组织CYP3192A1基因表达水平受性别的影响；a、b、c表示福寿螺不同组织CYP3192A1基因表达差异。Figure 4. Sex-dependent and tissue-dependent differences of P.canaliculata CYP3192A1 gene expressionA and B represent CYP3192A1 expression level affected by gender in different tissues; a, b and c represent tissue-dependent differences of CYP3192A1 expression in male and female individuals.福寿螺CYP3192A1表达具有组织差异性，且在雌、雄个体中的表达有差异。雄螺CYP3192A1在肠组织表达量最高，显著高于其他组织，而其他组织间无显著性差异；雌螺CYP3192A1在肠和鳃组织表达量最丰富，其次是肝和胃组织，心和腹足组织内表达量最低。
2.4 四聚乙醛处理对CYP3192A1基因表达的影响
分别以0.6 mg·L-1、2.0 mg·L-1的四聚乙醛处理雌性福寿螺24 h，RT-qPCR检测不同时刻肝、鳃、肠和胃组织的表达情况。结果如图 5所示，四聚乙醛的处理质量浓度为0.6 mg·L-1时，肝和胃组织CYP3192A1表达量在0~6 h呈逐渐上升趋势；鳃和肠组织CYP3192A1表达量在2~6 h呈逐渐上升趋势。四聚乙醛处理质量浓度为2.0 mg·L-1时，肝和胃组织CYP3192A1表达量与0.6 mg·L-1时变化趋势类似，但是表达量上升时间缩短至2 h，表达量峰值降低；鳃和肠组织CYP3192A1表达量实验期间呈持续下降趋势，未检测到CYP3192A1的诱导表达。四聚乙醛诱导福寿螺CYP3192A1表达的持续时间随浓度变化在时间维度上表现出不同的趋势：低浓度诱导期更长，高浓度诱导期缩短甚至在实验时间内未检测到诱导作用。
![]()
图 5 CYP3192A1基因在不同浓度四聚乙醛诱导下各组织中表达变化a.肝；b.鳃；c.肠；d.胃Figure 5. CYP3192A1 expression in different tissues of P.canaliculata exposed to various concentrations of methaldehydea. liver; b. gill; c. intestine; d. stomach3. 讨论
该研究利用在线工具对克隆基因的编码蛋白序列进行分析鉴定，在氨基酸序列中未发现明显的跨膜信号肽，表明该蛋白质不属于分泌蛋白；跨膜结构区分析结果显示CYP3192A1具有由12~34号氨基酸组成的跨膜螺旋，镶嵌于生物膜内，35~517号氨基酸位于细胞质基质内，据此推测CYP3192A1可能依赖于自身的跨膜螺旋锚定于内质网上[1]。
细胞色素P450是一类含有亚铁血红素的超家族单加氧酶，参与多类底物的不同催化反应类型，但是不同形态的细胞色素P450具有位点特异性、立体结构特异性和底物特异性[17-18]，引起这种特异性的结构基础归结于氨基酸和核酸序列的差异性。研究认为这种差异性是细胞色素P450超家族通过基因扩增和适应多样化而形成的，另外细胞色素P450也具有特定的保守序列[19]。细胞色素P450结构的一般特征为C端富集螺旋状物且高度保守，N端富集β折叠物；CYPs标志性的区域为F**G***C*G基序，最保守的部分是围绕血红素的蛋白质中心，该保守中心空间上接近螺旋体I及螺旋K[20]；CYPs的可变区通常与N端锚定、底物结合和识别有关，具有底物识别位点(SRSs)[21]。该研究中的蛋白质序列具有Cyt_P450_ E_ grp-I(Cytochorome P450 E class group I)(IPR002401)的特征结构域，由于大多数微粒体细胞色素P450均属于此分类[2]研究基础，推测CYP3192A1蛋白序列属于微粒体细胞色素P450酶。另外CYP3192A1蛋白序列具有细胞色素P450酶的高度保守区域，空间结构呈球形，具有血红素结合区及蛋白中心区域；N端结构变异性较大，空间上没有结构域包围，易于三维结构的变化，符合细胞色素的一般结构特征。系统发育分析显示，与福寿螺CYP3192A1同源性最高的序列是长牡蛎的CYP3A24，同源性仅为34%，并未达到同家族序列同源性所规定的40%；但GenBank数据库中与此同源性较高的其他物种CYPs序列均为CYP3A亚家族的成员，且空间结构高度保守，而研究认为细胞色素P450超家族的多样性是基因扩增进化而来的[4]，推测福寿螺CYP3192A1基因可能是由CYP3A进化而来，其在蛋白结构上呈现高度相似性，因此CYP3192A1在结构及功能研究上可以参考CYP3A亚族基因的研究成果。研究认为CYP3A亚家族酶系是药物代谢中最重要且唯一具多个底物识别位点的CYP亚家族同工酶，其代谢底物种类繁多[13]；目前，在人、老鼠、鱼类等多种生物体内已发现与药物代谢相关的CYP3A亚家族成员[22-24]。因此福寿螺CYP3192A1序列结构特征的克隆与鉴定，与CYP3A亚家族结构与功能研究相互促进，为福寿螺细胞色素P450酶学特性及药物代谢研究提供重要理论基础。
细胞色素P450广泛存在于不同的生物个体中，且同种基因在不同生物、不同组织中分布具有差异性[25]；水生动物研究显示鳃是异源物质吸收、生物转化及代谢物排泄的重要器官[26]，消化器官因其细胞色素P450的高水平表达而作为主要的解毒器官研究[27]，因此该实验以鳃和消化系统作为实验研究对象组织。研究认为细胞色素P450基因的表达是在体内的信号或外界异生物质的复杂而精确地调控下进行的，其表达水平除了与物种和组织有关，还应考虑其性别及其所处环境条件等情况[19]。该研究中，性别和组织特异性对福寿螺CYP3192A1的表达水平有显著影响。整体表现为雌性表达量高于雄性，这可能与体内激素分泌相关；不同性别个体的CYPs的表达量不同，并且特定CYPs在体内的表达及其活性也存在性别差异性[28-29]；而在组织间，表现为表达量最丰富的组织是肠和鳃，肝和胃次之，在心和腹足中表达最少。这与常见水生动物组织特异性研究有相似之处，但也有诸多差异性；如异育银鲫(Carassius auratus gibelio)CYP3A亚家族基因在肝和肠组织中转录水平最高，在肾和鳃中次之，在其余组织较低[30]。
检测代谢酶系的cDNA表达水平是研究化合物对CYPs的抑制或诱导效应的常用方法之一[31]。为了解福寿螺CYPs表达水平与四聚乙醛处理的相关性，文章在四聚乙醛对福寿螺急性攻毒实验的基础上，参考韩文素等[32]的实验方法选取不同倍数LC50的亚致死浓度的四聚乙醛处理福寿螺，检测了四聚乙醛处理后CYP3192A1基因表达的变化。研究认为胁迫环境下福寿螺几乎不活动，新陈代谢维持在很低水平[33]，而该实验结果显示在福寿螺受到四聚乙醛胁迫后活力下降的同时，不同组织均出现了CYP3192A1表达量一定时期内升高的现象，说明四聚乙醛能增强福寿螺组织CYP3192A1的表达，总体表现为“低浓度诱导，高浓度缩短诱导期”；另外，随着四聚乙醛浓度的升高，福寿螺CYP3192A1表达水平也发生显著的变化，证明了四聚乙醛处理与福寿螺CYP3192A1表达量变化具有密切相关性。可能由于四聚乙醛在福寿螺体内存在药物的“首过效应”[34]，最新接触到药物的组织对药物存在吸收和稀释作用，鳃和肠诱导期更早出现；最后所有组织CYP3192A1表达量均表现为持续下降的趋势，可能是由于福寿螺自身解毒系统并不能完全消除四聚乙醛的毒害作用，机体生命活动状态趋向衰竭，高浓度四聚乙醛能够增强四聚乙醛作用效果。
总体来说四聚乙醛能够诱导福寿螺组织内CYP3192A1表达水平的增加，研究认为由表达增加或活力增强等引起的代谢酶酶活性增加一般会促使药物代谢增加[35]，因此CYP3192A1在药物胁迫下的诱导表达可能在增强四聚乙醛代谢中起到作用，然而这需要对四聚乙醛在福寿螺体内代谢途径及CYP3192A1底物特异性进行深入研究验证。
-
![]()
图 2 养殖水温度(a)、溶解氧(a)、盐度(b)、pH(b)、氧化还原电位(c)和主要营养盐(d)随时间变化趋势
Figure 2. Variation trends of temperature (a), dissolved oxygen (a), salinity (b), potential of hydrogen (b), redox potential (c) and environmental factors (d) in aquaculture water along with time
表 1 弧菌密度与环境因子的Pearson相关系数
Table 1 Pearson correlation coefficients of Vibrio density and environmental parameters
弧菌密度Vibrio density 高温期(T>25 ℃)
high-temperature period低温期(T≤25 ℃)
low-temperature period温度/℃ temperature 0.667* 0.211 ρ(溶氧)/mg·L-1 dissolved oxygen -0.182 0.007 酸碱度potential of hydrogen 0.137 -0.002 盐度salinity -0.484 -0.180 氧化还原电位/mV redox potential 0.832** 0.011 ρ(亚硝酸盐)/μg·L-1 nitrite 0.216 -0.316 ρ(硝酸盐)/μg·L-1 nitrate 0.771** -0.214 ρ(氨盐)/μg·L-1 ammonium -0.781** -0.153 ρ(磷酸盐)/μg·L-1 phosphate -0.072 -0.399 注：*. 显著性相关(P＜0.05)；* *. 极显著性相关(P＜0.01)
Note: *. significant at 0.05 level (P＜0.05)；* *. very significant at 0.01 level (P＜0.01)
下载: 导出CSV
表 2 弧菌密度和水温的模型概述和参数估计值
Table 2 Model summary and parameter estimates of Vibrio density and water temperature
方程
equation模型概述model summary 参数估计值estimated value of parameter R2 F df1 df2 Sig. α α1 α2 α3 线性方程linear 0.445 8.024 1 10 0.018 -15 289.182 590.942 对数方程logarithmic 0.450 8.173 1 10 0.017 -54 643.412 16 776.839 倒数方程inverse 0.454 8.325 1 10 0.016 18 268.378 -476 129.359 二次方程quadratic 0.445 8.024 1 10 0.018 -15 289.182 590.942 0 三次方程cubic 0.705 10.74 2 9 0.004 -645 009.140 34 045.683 0 -13.968 复合方程compound 0.455 8.352 1 10 0.016 0.002 1.615 幂方程power 0.460 8.507 1 10 0.015 2.385E-17 13.602 增长方程growth 0.455 8.352 1 10 0.016 -6.369 0.479 指数方程exponential 0.455 8.352 1 10 0.016 0.002 0.479 逻辑方程logistic 0.455 8.352 1 10 0.016 583.579 0.619 注：R2表示决定系数；F表示方差齐性检验；df1为F统计量分子的自由度，df2为F统计量分母的自由度；Sig.表示统计学意义；后表同此；α表示方程常数；α1、α2、α3表示方程系数
Note：R2 is coefficient of determination；F is homogeneity of variance；df1 is freedom of F statistics molecules；df2 is freedom of F statistics denominator；Sig. means statistical significance；the same case in the following tables；α is equation constant；α1，α2 and α3 are equation coefficients.
下载: 导出CSV
表 3 弧菌密度和氧化还原电位的模型概述和参数估计值
Table 3 Model summary and parameter estimates of Vibrio density and redox potential
方程
equation模型概述model summary 参数估计值estimated value of parameter R2 F df1 df2 Sig. β β1 β2 β3 线性方程linear 0.692 22.497 1 10 0.001 474.345 33.940 对数方程logarithmic 0.614 15.938 1 10 0.003 -596.567 631.913 倒数方程inverse 0.491 9.645 1 10 0.011 1 838.289 -8 618.092 二次方程quadratic 0.706 10.831 2 9 0.004 721.987 4.869 0.603 三次方程cubic 0.717 6.758 3 8 0.014 1 126.761 -68.851 4.153 -0.049 复合方程compound 0.700 23.283 1 10 0.001 611.936 1.028 幂方程power 0.630 17.029 1 10 0.002 255.211 0.513 增长方程growth 0.700 23.283 1 10 0.001 6.417 0.027 指数方程exponential 0.700 23.283 1 10 0.001 611.936 0.027 逻辑方程logistic 0.700 23.283 1 10 0.001 0.002 0.973 注：β表示方程常数；β1、β2、β3表示方程系数
Note: β is equation constant；β1，β2 and β3 are equation coefficients.
下载: 导出CSV
表 4 弧菌密度和硝酸盐的模型概述和参数估计值
Table 4 Model summary and parameter estimates of Vibrio density and nitrate
方程
equation模型概述model summary 参数估计值estimated value of parameter R2 F df1 df2 Sig. γ γ1 γ2 γ3 线性方程linear 0.595 14.692 1 10 0.003 -242.475 21.662 对数方程logarithmic 0.545 11.990 1 10 0.006 -4 256.247 1 315.419 倒数方程inverse 0.491 9.660 1 10 0.011 2 452.585 -75 065.398 二次方程quadratic 0.695 10.276 2 9 0.005 2 726.633 -76.043 0.734 三次方程cubic 0.697 10.370 2 9 0.005 1 210.742 0 -0.477 0.006 复合方程compound 0.611 15.712 1 10 0.003 339.887 1.018 幂方程power 0.565 12.990 1 10 0.005 12.758 1.074 增长方程growth 0.611 15.712 1 10 0.003 5.829 0.018 指数方程exponential 0.611 15.712 1 10 0.003 339.887 0.018 逻辑方程logistic 0.611 15.712 1 10 0.003 0.003 0.983 注：γ表示方程常数；γ1、γ2、γ3表示方程系数
Note: γ is equation constant；γ1，γ2 and γ3 are equation coefficients.
下载: 导出CSV
表 5 弧菌密度和氨盐的模型概述和参数估计值
Table 5 Model summary and parameter estimates of Vibrio density and ammonium
方程
equation模型概述model summary 参数估计值estimated value of parameter R2 F df1 df2 Sig. k k1 k2 k3 线性方程linear 0.610 15.660 1 10 0.003 2 783.214 -15.626 对数方程logarithmic 0.661 19.542 1 10 0.001 8 511.565 -1 597.177 倒数方程inverse 0.702 23.609 1 10 0.001 -450.977 154 462.864 二次方程quadratic 0.757 14.028 2 9 0.002 6 901.759 -102.574 0.425 三次方程cubic 0.757 14.028 2 9 0.002 6 901.759 -102.574 0.425 0 复合方程compound 0.627 16.827 1 10 0.002 3 976.426 0.987 幂方程power 0.674 20.629 1 10 0.001 407 596.730 -1.292 增长方程growth 0.627 16.827 1 10 0.002 8.288 -0.013 指数方程exponential 0.627 16.827 1 10 0.002 3 976.426 -0.013 逻辑方程logistic 0.627 16.827 1 10 0.002 0 1.013 注：k表示方程常数；k1、k2、k3表示方程系数
Note: k is equation constant；k1，k2 and k3 are equation coefficients.
下载: 导出CSV
-
[1] JENNIFER L, HSIEH J, FRIES J S, et al. Vibrio and phytoplankton dynamics during the summer of 2004 in a eutrophying estuary[J]. Ecol Appl, 2007, 17(5): 102-109. doi: 10.1890/05-1274.1
[2] NUMATA K, DOIY. Biosynthesis of polyhydroxyalkanaotes by a novel facultatively anaerobic Vibrio sp. under marine conditions[J]. Mar Biotechnol, 2012, 14(13): 323-331. doi: 10.1007/s10126-011-9416-1
[3] RAMAIAH N, HILL R T, CHUN J, et al. Use of a chiA probe for detection of chitinase genes in bacteria from the Chesapeake Bay[J]. Fems Microbiol Ecol, 2000, 34(1): 63-71. https://academic.oup.com/femsec/article/34/1/63/617866
[4] SHERR E B, SHERR B F. Significance of predation by protists in aquatic microbial food webs[J]. Anton Leeuw Int J G, 2002, 81(1/2/3/4): 293-308. doi: 10.1023/A:1020591307260
[5] KARL D, MICHAELS B, BERGMAN D, et al. Dinitrogen fixation in the world's oceans[J]. Biogeochemistry, 2002, 57(1): 47-98. doi: 10.1023/A:1015798105851
[6] 陈志胜, 吕军仪, 吴金英, 等. 杂色鲍(Haliotis diversicolor)溃疡症病原菌的研究[J]. 热带海洋学报, 2000, 19(3): 72-77. doi: 10.3969/j.issn.1009-5470.2000.03.011 [7] NICOLAS J L, BASUYAUX O, MAZURIE J, et al. Vibrio carchariae, a pathogen of the abalone Haliotis tuberculata[J]. Dis Aquat Org, 2002, 50(1): 35-43. https://pubmed.ncbi.nlm.nih.gov/12152903/
[8] 刘广锋, 周世宁, 徐力文, 等. 杂色鲍幼苗"急性死亡脱落症"病原菌分析[J]. 中国水产科学, 2006, 13(4): 655-661. doi: 10.3321/j.issn:1005-8737.2006.04.023 [9] CAI J P, HAN Y, WANG Z. Isolation of Vibrio parahaemolyticus from abalone (Haliotis diversicolor supertexta L.) postlarvae associated with mass mortalities[J]. Aquaculture, 2006, 257(3): 161-166. https://www.sciencedirect.com/science/article/pii/S0044848606002110
[10] BALCÁZAR J L, GALLO-BUENO A, PLANAS M, et al. Isolation of Vibrio alginolyticus and Vibrio splendidus from captive-bred seahorses with disease symptoms[J]. Anton Leeuw Int J G, 2010, 97(2): 207-210. doi: 10.1007/s10482-009-9398-4
[11] MAHMUD Z H, NEOGI S B, KASSU A, et al. Occurrence, seasonality and genetic diversity of Vibrio vulnificus in coastal seaweeds and water along the Kii Channel, Japan[J]. Fems MicroBiol Ecol, 2008, 64(2): 209-218. doi: 10.1111/j.1574-6941.2008.00460.x
[12] SCHETS F M, MARCHESE A, GARBOM S. Potentially human pathogenic vibrios in marine and fresh bathing waters related to environmental conditions and disease outcome[J]. Int J Hyg Environ Health, 2011, 214(5): 399-406. doi: 10.1016/j.ijheh.2011.05.003
[13] 王振华. 循环水养鱼系统水质参数关系分析[J]. 中国农学通报, 2014, 30(8): 57-62. doi: 10.11924/j.issn.1000-6850.2013-2633 [14] MAUGERI T L, CACCAMO D, GUGLIANDOLO C. Potentially pathogenic vibrios in brackish waters and mussels[J]. J Gen Appl Microbiol, 2000, 89(2): 261-266. doi: 10.1046/j.1365-2672.2000.01096.x
[15] LAFISCA A, PEREIRA C S, GIACCONE V, et al. Enzymatic characterization of Vibrio alginolyticus strains isolated from bivalves harvested at Venice Lagoon (Italy) and Guanabara Bay (Brazil)[J]. Rev Inst Med Trop S Paulo, 2008, 50(4): 199-202. doi: 10.1590/S0036-46652008000400002
[16] 吴玲, 张景隆, 付利军. 中国沿海十省、市、自治区水域中致病性弧菌的分布调查[J]. 中华航海医学与高气压医学杂志, 2005, 12(3): 171-172. doi: 10.3760/cma.j.issn.1009-6906.2005.03.020 [17] 李和阳, 陈明霞, 郑天凌, 等. 深圳海域水体和九龙江口沉积物中TCBS菌群与弧菌相关性的研究[J]. 台湾海峡, 2011, 30(3): 394-399. doi: 10.3969/J.ISSN.1000-8160.2011.03.015 [18] 陈金玲, 赖秋明, 苏树叶, 等. 精养虾池主要生态因子变化特点与相关性分析[J]. 南方水产科学, 2012, 8(4): 49-56. doi: 10.3969/j.issn.2095-0780.2012.04.008 [19] 杜虹, 黄显兵, 郑兵, 等. 粤东深澳湾养殖区域异养细菌和弧菌的动态分布[J]. 热带海洋学报, 2010, 29(6): 110-117. doi: 10.3969/j.issn.1009-5470.2010.06.017 [20] 李卓佳, 李烁寒, 杨莺莺, 等. 凡纳滨对虾高位池养殖水体细菌变动及其与理化因子的关系[J]. 南方水产, 2010, 6(4): 6-12. doi: 10.3969/j.issn.1673-2227.2010.04.002 [21] 张新中, 谢珍玉, 王世锋, 等. 海南万宁小海近岸表层海水弧菌组成及其变化[J]. 海洋学报, 2011, 33(6): 142-150. https://www.cnki.com.cn/Article/CJFDTOTAL-SEAC201106018.htm [22] 倪纯志, 林燕顺, 叶德赞, 等. 海水养虾池的几种致病弧菌生态[J]. 台湾海峡, 1992, 14(1): 73-79. https://www.cnki.com.cn/Article/CJFDTotal-TWHX199501012.htm [23] DEPAOLA A, HOPKINS L H, PEELER J T, et al. Incidence of Vibrio parahaemolyticus in U.S. coastal waters and Ostrea[J]. Appl Environ Microb, 1990, 56(8): 2299-2302. doi: 10.1128/aem.56.8.2299-2302.1990
[24] DEPAOLA A, NORDSTROM J L, BOWERS J C, et al. Seasonal abundance of total and pathogenic Vibrio parahaemolyticus in Alabama oysters[J]. Appl Environ Microb, 2003, 69(3): 1521-1526. doi: 10.1128/AEM.69.3.1521-1526.2003
[25] 李涛, 宁喜斌. 副溶血性弧菌温度-盐度双因素预测模型的建立[J]. 微生物学通报, 2009, 36(8): 1200-1205. https://www.cnki.com.cn/Article/CJFDTotal-WSWT200908023.htm [26] 王璐华, 宁喜斌. 副溶血性弧菌的温度预测模型[J]. 食品与生物技术学报, 2009, 28(2): 263-266. doi: 10.3321/j.issn:1673-1689.2009.02.025 [27] 黄和, 田金玲, 励建荣. 波动温度模拟蟹肉副溶血性弧菌生长模型[J]. 食品研究与开发, 2011, 32(11): 1-3. doi: 10.3969/j.issn.1005-6521.2011.11.001 [28] 林强, 李宁求, 付小哲, 等. 牡蛎养殖过程中副溶血弧菌与水质因子间的关系[J]. 水产学报, 2012, 36(3): 415-421. doi: 10.3724/SP.J.1231.27656 [29] 吴建平, 蔡创华, 周毅频, 等. 大亚湾网箱养殖区异养细菌和弧菌的数量动态[J]. 湛江海洋大学学报, 2006, 26(3): 21-25. doi: 10.3969/j.issn.1673-9159.2006.03.006 [30] 章洁香, 张瑜斌, 张才学, 等. 高位虾池水体细菌和弧菌的密度变化及影响因素[J]. 集美大学学报: 自然科学版, 2010, 15(5): 12-17. https://lib.cqvip.com/Qikan/Article/Detail?id=35312122 -
期刊类型引用(3)
1. 凌慧，赵旺，王瑞旋，黄晟，刘晓靖，吴开畅，王江勇. 方斑东风螺养殖环境中细菌数量变动与环境因子的相关性和主成分分析. 生态科学. 2019(02): 130-137 . 
百度学术
2. 孙永婵，王瑞旋，赵曼曼，曹超，鲁义善，姚托，王江勇. 鲍消化道及其养殖水体异养菌的耐药性研究. 南方水产科学. 2017(03): 58-65 . 本站查看
3. 唐姝，熊金波，王凯，胡常巨，郑嘉来，张德民. 象山港海水网箱养殖区弧菌群落动态变化. 生物学杂志. 2017(04): 76-81 . 百度学术
其他类型引用(7)
计量
- 文章访问数: 3056
- HTML全文浏览量: 122
- PDF下载量: 1543
- 被引次数: 10
粤公网安备 44010502001741号
