广西银滩南部海域海洋牧场鱼类群落结构特征及其与环境因子的关系

于杰, 邹建伟, 陈国宝

于杰, 邹建伟, 陈国宝. 广西银滩南部海域海洋牧场鱼类群落结构特征及其与环境因子的关系[J]. 南方水产科学, 2024, 20(5): 81-90. DOI: 10.12131/20240102
引用本文: 于杰, 邹建伟, 陈国宝. 广西银滩南部海域海洋牧场鱼类群落结构特征及其与环境因子的关系[J]. 南方水产科学, 2024, 20(5): 81-90. DOI: 10.12131/20240102
YU Jie, ZOU Jianwei, CHEN Guobao. Characteristics of fish community structure and its relationship with environmental factors in marine ranching zone in southern area of Yintan in Guangxi[J]. South China Fisheries Science, 2024, 20(5): 81-90. DOI: 10.12131/20240102
Citation: YU Jie, ZOU Jianwei, CHEN Guobao. Characteristics of fish community structure and its relationship with environmental factors in marine ranching zone in southern area of Yintan in Guangxi[J]. South China Fisheries Science, 2024, 20(5): 81-90. DOI: 10.12131/20240102

广西银滩南部海域海洋牧场鱼类群落结构特征及其与环境因子的关系

基金项目: 北海市政府采购项目 (BHZC2022-C3-00009-XYZJ);广东省重点领域研发计划项目 (2020B1111030002);农业农村部渔业发展补助资金项目 (SCSFRI2024A012)
详细信息
    作者简介:

    于 杰 (1980—),女,副研究员,博士,研究方向为海洋渔业资源和海洋遥感。E-mail: yujiescs@aliyun.com

    通讯作者:

    陈国宝 (1975—),男,研究员,硕士,研究方向为渔业资源及渔业声学。E-mail: chengb98@163.com

  • 中图分类号: S 931.9

Characteristics of fish community structure and its relationship with environmental factors in marine ranching zone in southern area of Yintan in Guangxi

  • 摘要:

    适时评价国家级海洋牧场示范区的鱼类资源状况,对推进海洋牧场建设具有重要意义。基于2022年4月广西银滩南部海域海洋牧场渔业资源调查数据,分析了鱼类资源特征、群落结构及其与环境因子的相关性,以为北部湾海洋牧场生态效应研究提供参考。结果显示,示范区鱼类数量密度、生物量密度和种类数均高于对照区。示范区和对照区鱼类物种组成在种级水平上为中等不相似,且2个区域鱼类群落均受到一定程度的干扰。研究区鱼类群落可划分为4个具有显著性差异的类群 (p<0.05),分别位于示范区B区和对照区南部站位 (类群I)、对照区北部站位 (类群II)、示范区A区和B区交界处 (类群III) 和示范区A区 (类群IV)。类群I和III的平均相异性指数为53.88%,贡献率排前3位的分歧种为吉打副叶鲹 (Alepes djedaba)、蓝圆鲹 (Decapterus maruadsi)和竹荚鱼 (Trachurus japonicus);类群III和IV的平均相异性指数为48.85%,贡献率排前3位的分歧种为竹荚鱼、短吻鲾 (Leiognathus brevirostris) 和花斑蛇鲻 (Saurida undosquamis)。水温和化学需氧量是显著影响研究区鱼类群落结构的主要环境因子,溶解氧和浮游植物密度的影响稍弱。

    Abstract:

    Timely assessment of fish resources in state-level marine ranching demonstration zone is important in promoting the construction of marine ranching. Based on the fishery resources survey data in the southern area of Guangxi Yintan in April 2022, we analyzed the characteristics of fish resources, community stability, fish community structure and their correlation with environmental factors to provide references for research on the ecological effects of marine ranching in Beibu Gulf. The results show that the fish individual density, biomass density and number of species in the demonstration zone were higher than those in the control zone. The fish species composition in the demonstration and control zones was moderately dissimilar at species level. In addition, fish communities in both zones were somewhat disturbed. The fish community in the survey area was divided into four significantly different clusters (p<0.05), which were located in the southern station of the demonstration zone B and the control zone (Cluster I), the northern station of the control zone (Cluster II), the junction of the demonstration zone A and B (Cluster III), and the demonstration zone A (Cluster IV). The average dissimilarity index of Cluster I and III was 53.88%, and the divergent species that ranked in the top three in terms of their contribution were Alepes djedaba, Decapterus marauds and Trachurus japonicus, while the average dissimilarity index of Cluster III and IV was 48.85%, with the top three divergent species (T. japonicus, Leiognathus brevirostri and Saurida undosquamis) contributing to the dissimilarity index. Water temperature and chemical oxygen demand were the main environmental factors significantly influencing the structure of the fish community in the survey area, while dissolved oxygen and phytoplankton density had a slightly weaker effect.

  • 鳗败血假单胞菌 (Pseudomonas anguilliseptica, PA) 于1972年由 Wakabayashi和Egusa[1]在日本鳗鲡 (Anguilla japonica) 上首次分离并报道,其后发现该菌主要影响欧洲地区的欧洲鳗鲡 (A. anguilla)[2]、大西洋鲷 (Sparus aurata)[3]、大西洋鳕 (Gadus morhua)[4]、圆鳍鱼 (Cyclopterus lumpus)[5]等海洋鱼类。因其感染的流行季节主要在水温低于12 ℃的冬季,且被感染宿主呈现体表出血性红斑的特征,因此被称为“冬季病”[6]和“红斑病”[7]。近年来,笔者研究发现PA是大口黑鲈 (Micropterus salmoides) 的一种新病原菌 (另文报道),对大口黑鲈的健康养殖造成威胁,但目前其在我国大口黑鲈中的流行情况尚不完全清楚。

    PA的生长条件苛刻、培养时间长[8],同时其感染的临床症状与嗜水气单胞菌 (Aeromonas hydrophila)[9]、维氏气单胞菌 (A. veronii)[10]等相似,通过病原菌分离鉴定及临床症状的传统方法难以进行快速诊断,往往出现漏检、错检、误诊的情况,影响了PA感染的有效防控。建立快速、特异的检测方法不仅对PA的早期临床诊断和防控具有重要意义,还能为开展流行病学调查提供条件。目前针对水生动物细菌性病原的检测方法主要有环介导等温扩增技术 (LAMP)[11]、核酸序列依赖性扩增检测技术 (NASBA)[12]、免疫磁珠分离技术 (IMS)[13]、酶联免疫吸附测定技术 (ELISA)[14]、重组酶介导等温扩增技术 (RAA)[15]、实时荧光定量PCR (qPCR)[16]等,综合考虑经济性、操作简便性、时效性等因素后,本研究采用SYBR Green I 实时荧光定量PCR (SYBR Green I real-time quantitative PCR) 和重组酶介导等温扩增技术结合侧流层析试纸条(Recombinase-mediated isothermal amplification combined with lateral flow dipstick, RAA-LFD),以PA的管家基因recA为靶基因建立2种检测方法,以期为PA的早期、快速检测提供方法参考。

    健康大口黑鲈 [外部无异常、PA特异性检测结果为阴性,体质量为 (18.4±1.5) g)]购自四川某苗种生产企业;各菌株 (表1) 由本实验室保存。主要试剂为:pUC18-T载体 (宝生物工程有限公司,大连);细菌基因组DNA提取试剂盒、胶回收试剂盒及质粒提取试剂盒 (天根生化科技有限公司,北京);Taq Pro Universal SYBR qPCR Master Mix、2×Rapid Taq Master Mix (诺唯赞生物科技有限公司,南京);RAA核酸扩增试剂盒及侧流层析试纸条 (乐尚生物科技有限公司,无锡)。

    表  1  研究所用的细菌菌株
    Table  1  Bacterial strains used in this study
    菌株名称
    Bacterial strain
    菌株编号
    Strain No.
    鳗败血假单胞菌 Pseudomonas anguilliseptica SICAU PA001
    无乳链球菌 Streptococcus agalactiae GY104
    拟态弧菌 Vibrio mimicus SCCF01
    嗜水气单胞菌 Aeromonas hydrophila YYL
    鮰爱德华菌 Edwardsiella ictaluri CNHYC01
    鰤诺卡菌 Nocardia seriolae DY01
    荧光假单胞菌 P. fluorescens ATCC®13525
    柱状黄杆菌 Flavobacterium cloumnare ATCC®23463
    维氏气单胞菌 A. veronii ATCC®35624
    下载: 导出CSV 
    | 显示表格

    根据Mjølnerød等[17]对53株PA管家基因的分析结果,选择其中具有完全同源性的recA基因 (登录号:MW685241.1—MW685293.1) 作为建立2种检测方法的靶基因,并采用验证recA基因的特异性引物 (recA-PCR-F/R) 作为本研究PCR和构建标准品质粒的引物。使用Primer Premier 5.0软件,根据SYBR Green I real-time quantitative PCR和RAA引物设计原则设计多对引物,并设计RAA探针。成功筛选出引物recA-qPCR-F/R、recA-RAA-F/R和探针recA-RAA-Probe,并在recA-RAA-R的5' 端标记生物素。本研究所用引物和探针 (表2) 均由生工生物工程 (上海) 有限公司合成。将recA-RAA-F/R配合recA-RAA-Probe的RAA扩增产物和recA-qPCR-F/R的PCR 扩增产物送至生工生物工程 (上海) 有限公司测序,在GenBank中获取相应序列并与测序结果进行比对,以验证引物的特异性。

    表  2  研究所用引物及探针
    Table  2  Primers used in this study
    检测方法
    Detection method
    引物名称
    Primer name
    序列 (5'—3')
    Sequence (5'−3')
    产物大小
    Length of product/bp
    PCRrecA-PCR-F/RF: CTCGCTTGGTCTGGACATCG
    R: CCTTGCCTTGGCCGATCTT
    740
    qPCRrecA-qPCR-F/R F: GGCAGCTCTTACACCCAAAG
    R: GGCTCCCAAACATAACACCTAT
    179
    RAA-LFDrecA-RAA-F/RF: CATGCACTTGACCCAATCTACGCAGCAAAG
    R: Biotin CCCGTAATCTTACGCATTGCCTGACTCATC
    251
    recA-RAA-ProbeFAM AAGCTGAAATCGAAGGTGAAATTGGGGATA/idSp/
    GCACGTTGGGTTGGC C3 Spacer
    下载: 导出CSV 
    | 显示表格

    以PA的基因组DNA为模板,用recA-PCR-F/R引物扩增recA基因目的片段,与pUC18-T载体连接后转化于DH5α感受态细胞,获得重组标准品质粒pUC18-recA。用微量核酸蛋白检测仪测定质粒浓度和纯度,并根据摩尔定律计算出单位体积质粒所含的DNA拷贝数。

    2种方法均采用2.816×106 拷贝·μL−1的标准品质粒进行优化。qPCR引物终浓度梯度设置为0.2~1.0 μmol·L−1,退火温度梯度为 58~64 ℃,以获得qPCR最小的平均Ct值和标准差;RAA引物终浓度梯度设置为 6.25~200 μmol·L−1,RAA扩增反应温度梯度为29~54 ℃,并通过观察不同条件下侧流层析试纸条第30分钟时的显色情况确定RAA扩增的最佳反应体系和反应条件。

    qPCR方法的建立:将标准品质粒10倍稀释成9组浓度梯度 (2.816×109~2.816×101 拷贝·μL−1),采用优化后的反应条件 (反应体系总体积10 μL:recA-qPCR-F/R各0.2 μL,2×Taq Pro Universal SYBR qPCR Mix 5 μL,ddH2O 3.6 μL,模板DNA 1 μL;反应条件为:95 ℃预变性30 s;95 ℃变性10 s,60 ℃退火20 s,72 ℃延伸10 s,40个循环) 进行qPCR检测,建立以质粒拷贝数对数值为x轴和以Ct值为y轴的qPCR标准曲线;每个浓度梯度的标准品质粒设置3次qPCR检测,以进行组内重复性分析,并以相同的标准品质粒为模板,在不同时间点开展3次独立重复实验,进行组间重复性分析;通过Ct值变异系数评价方法的重复性。

    RAA-LFD方法的建立:采用优化后的反应条件 (反应体系总体积50 μL:recA-RAA-Probe 0.6 μL,recA-RAA-LF/LR各2.1 μL,基础缓冲液25 μL,模板DNA 2 μL,ddH2O 15.2 μL;各组分混匀后,加入含冻干酶粉的反应管中溶解,然后加入280 mmol·L−1的醋酸镁溶液3 μL作为启动剂混匀,39 ℃水浴20 min,立即置于冰上5 min) 对标准品质粒 (2.816×109 拷贝·μL−1) 进行RAA扩增,将侧流层析试纸条样品垫末端垂直插入含有50 μL稀释5倍的RAA核酸扩增产物的EP管中,观察试纸条C线和T线在30 min内的颜色变化,并参考Saxena等[18]的标准判定结果。

    各取101~108 拷贝·μL−1的标准品质粒1 μL作为模板,分别进行qPCR、RAA和普通PCR检测,分析3种方法的灵敏性。

    取质量浓度为10 ng·μL−1的各株细菌DNA 1 μL作为模板,同时进行qPCR、RAA和普通PCR扩增,分析3种方法的特异性。

    通过建立的qPCR和RAA-LFD方法以及普通PCR方法,对10份疑似感染PA的大口黑鲈肾组织样本与30份腹腔注射感染PA的大口黑鲈肾组织样本进行检测并比较检出率的差异。

    选取15尾健康大口黑鲈 (共3组,每组5尾),腹腔注射PA,于感染后的第7天分别采集每组病鱼的肝、脾、肾、心、肌肉、脑、眼和鳃组织并提取DNA,采用qPCR方法和标准曲线方程测定每组病鱼组织中的平均PA载量。

    PA提取的DNA通过qPCR引物recA-qPCR-F/R扩增后得到179 bp的单一目的条带 (图1-a),通过RAA引物recA-RAA-F/R配合recA-RAA-Probe扩增后得到251 bp的单一目的条带 (图1-b)。将目的条带的测序结果与GenBank中获取的相应序列进行比对,结果显示2段产物序列与其他序列的同源性均为100% (图2)。

    图  1  基于鳗败血假单胞菌recA基因的引物特异性与标准品质粒鉴定结果
    注:a. qPCR引物特异性结果;b. RAA引物及探针特异性结果;c. 标准品质粒鉴定结果。
    Fig. 1  Results based on primers specificity and standard plasmid identification of recA gene of P. anguilliseptica
    Note: a. Primer specificity results for qPCR; b. Primer and probe specificity results for RAA; c. Results of standard plasmid identification.

    标准品质粒扩增获得的目的片段大小为740 bp (图1-c),与预期大小相符。重组质粒pUC18-recA的质量浓度为1.057×102 ng·μL−1,经质粒拷贝数计算公式换算后为2.816×1010 拷贝·μL−1

    qPCR扩增条件的优化结果显示,qPCR引物浓度为0.6 μmol·L−1、退火温度为60 ℃时平均Ct值 (13.259) 和标准差 (0.045) 最小 (表3),故选择该引物浓度和退火温度为最佳反应条件。

    表  3  鳗败血假单胞菌qPCR扩增条件的优化
    Table  3  Optimization of qPCR amplification conditions for P. anguilliseptica
    退火温度
    Annealing temperature/℃
    引物浓度
    Primer concentration/(μmol·L−1)
    Ct值
    Ct value
    平均Ct值
    Mean Ct value
    标准差
    Standard deviation
    58 0.2 14.583 15.418 14.845 14.949 0.349
    0.4 13.876 14.356 14.230 14.154 0.203
    0.6 13.336 13.576 13.429 13.447 0.099
    0.8 13.745 13.401 13.313 13.486 0.187
    1.0 13.756 13.444 13.895 13.698 0.189
    60 0.2 14.419 15.074 14.355 14.616 0.325
    0.4 14.079 13.832 13.786 13.899 0.129
    0.6 13.204 13.259 13.314 13.259 0.045
    0.8 13.349 13.608 13.573 13.510 0.115
    1.0 13.749 14.023 13.988 13.920 0.122
    62 0.2 16.812 16.033 16.267 16.370 0.326
    0.4 15.305 15.466 15.077 15.283 0.160
    0.6 14.275 14.119 14.439 14.277 0.131
    0.8 14.206 14.401 15.046 14.551 0.359
    1.0 14.748 14.207 14.805 14.587 0.269
    64 0.2 15.745 16.784 15.695 16.075 0.502
    0.4 14.139 14.916 14.878 14.644 0.358
    0.6 15.061 15.328 15.457 15.282 0.165
    0.8 15.740 15.112 15.557 15.470 0.264
    1.0 15.963 15.369 15.844 15.725 0.257
    下载: 导出CSV 
    | 显示表格

    RAA扩增条件的优化结果如图3所示。试纸条显色第30分钟时的T线灰度值与引物浓度成正比 (图3-a—3-b),由于引物浓度过高可能会产生大量引物二聚体,且在25~100 μmol·L−1时T线灰度值差异不显著 (p>0.05),故选择50 μmol·L−1为最佳引物浓度。T线灰度值在反应温度介于29~54 ℃时呈正态分布,在39~44 ℃时达到最高,差异不显著 (p>0.05,图3-c—3-d),故选择其中获得最高T线灰度值的39 ℃为最适反应温度。

    图  2  鳗败血假单胞菌qPCR及RAA扩增序列同源性对比
    Fig. 2  Comparison of homology between qPCR and RAA amplified sequences of P. anguilliseptica
    图  3  鳗败血假单胞菌RAA扩增条件的优化
    注:a. 不同引物浓度的RAA-LFD显色结果;b. 灰度分析结果 (对应3-a);c. 不同反应温度的RAA-LFD显色结果;d. 灰度分析结果 (对应3-c);不同字母表示差异显著 (p<0.05)。
    Fig. 3  Optimization of RAA amplification conditions for P. anguilliseptica
    Note: a. Results of RAA-LFD color development with different primer concentrations; b. Grayscale analysis results (Corresponding to 3-a); c. Results of RAA-LFD color development at different reaction temperatures; d. Grayscale analysis results (Corresponding to 3-c); different letters represent significant differences (p<0.05).

    采用优化后的反应条件对9组10倍梯度稀释的标准品质粒 (2.816×109~2.816×101 拷贝·μL−1) 进行qPCR检测,建立质粒拷贝数对数值与Ct值对应关系的定量标准曲线 (图4-a),标准曲线介于2.816×107~2.816×102 拷贝·μL−1之间Ct值线性良好。标准曲线方程为y= –3.158x+32.213,r2=0.999 2。熔解曲线显示,6个浓度的标准品质粒在熔解温度为 (84±0.5) ℃时均出现特异性细窄单峰 (图4-b),表明扩增产物具有良好的特异性和均一性。标准曲线重复性检测结果表明,组内变异系数≤0.730%,组间变异系数≤1.443% (表4),该方法重复性和稳定性较好。

    图  4  鳗败血假单胞菌qPCR方法建立结果
    注:a. qPCR标准曲线;b. qPCR熔解曲线,1—6. 质粒标准品浓度为2.816×107~2.816×102 拷贝·μL−1
    Fig. 4  Results established by qPCR for P. anguilliseptica
    Note: a. Standard curve for qPCR; b. Melting curves for qPCR, 1−6. Concentration of standard plasmid was 2.816×107−2.816×102 copies·μL−1.

    采用优化后的反应条件对标准品质粒 (2.816×109 拷贝·μL−1) 进行RAA扩增,并通过侧流层析试纸条对扩增产物进行检测。结果显示,T线灰度值与显色时长成正比,在15 min后试纸条T线显色情况较好 (T线灰度值≥40.0),且各时间点无显著性差异 (p>0.05,图5-a—5-b),表明该方法最快15 min可读取检测结果,且显色较为稳定。

    表  4  鳗败血假单胞菌qPCR方法重复性试验结果
    Table  4  Results of qPCR reproducibility test for P. anguilliseptica
    标准品质粒浓度
    Standard quality plasmid
    concentration/
    (拷贝·μL−1)
    组内重复性试验
    Repeatability test within group
    组间重复性试验
    Repeatability test between group
    平均Ct值±标准差
    Mean Ct value±standard deviation
    变异系数
    Variable coefficient/%
    平均Ct值±标准差
    Mean Ct value±standard deviation
    变异系数
    Variable coefficient/%
    2.816×107 10.300±0.032 0.311 10.328±0.149 1.443
    2.816×106 13.213±0.042 0.318 13.257±0.122 0.920
    2.816×105 16.710±0.122 0.730 16.606±0.212 1.277
    2.816×104 19.511±0.102 0.523 19.523±0.111 0.569
    2.816×103 22.609±0.081 0.358 22.618±0.168 0.743
    2.816×102 26.083±0.130 0.498 26.120±0.210 0.804
    下载: 导出CSV 
    | 显示表格
    图  5  鳗败血假单胞菌RAA-LFD反应结果
    注:a. 不同反应时间的RAA-LFD显色;b. 灰度分析 (对应5-a),不同字母表示差异显著 (p<0.05)。
    Fig. 5  Results of RAA-LFD reaction of P. anguilliseptica
    Note: a. Results of RAA-LFD color development at different reaction time; b. Grayscale analysis results (Corresponding to 5-a), different letters represent significant differences (p<0.05).

    由普通PCR灵敏性试验结果可知,模板浓度低于2.816×105 拷贝·μL−1时,无法产生扩增条带 (图6-a),故普通PCR的最低检测浓度为2.816×105 拷贝·μL−1。由2.816×107 拷贝·μL−1浓度起始的6个10倍梯度稀释的标准品质粒扩增曲线可知,模板浓度在低于2.816×102 拷贝·μL−1时,无法产生扩增曲线 (图6-b),故qPCR的最低检测浓度为2.816×102 拷贝·μL−1。由RAA-LFD灵敏性试验结果可知,模板浓度低于2.816×104 拷贝·μL−1时扩增产物在试纸条上显色极不明显 (T线灰度值≤5.0,图6-c—6-d),故RAA-LFD的最低检测浓度为2.816×104 拷贝·μL−1。综上,构建的2种方法灵敏性均高于普通PCR,且qPCR方法的灵敏性最高。

    图  6  鳗败血假单胞菌灵敏性试验结果
    注:a. 普通PCR灵敏性电泳图,M. 标准分子量DL 2000,NC. 阴性对照,1. 2.816×108 拷贝·μL−1,2. 2.816×107 拷贝·μL−1,3. 2.816×106拷贝·μL−1,4. 2.816×105 拷贝·μL−1,5. 2.816×104 拷贝·μL−1;b. qPCR标准质粒扩增曲线,1—6. 质粒标准品浓度为2.816×107~2.816×102 拷贝·μL−1;c. 不同浓度标准质粒的RAA-LFD显色;d. 灰度分析 (对应6-c),不同字母表示差异显著 (p<0.05)。
    Fig. 6  Results of sensitivity test for P. anguilliseptica
    Note: a. Sensitive electrophoresis pattern of PCR, M. Standard molecular weight is DL 2000, NC. Negative control, 1. 2.816×108 copies·μL−1, 2. 2.816×107 copies·μL−1, 3. 2.816×106copies·μL−1, 4. 2.816×105 copies·μL−1, 5. 2.816×104 copies·μL−1; b. Amplification curves of qPCR standard plasmids, 1−6. Concentration of standard plasmid was 2.816×107−2.816×102 copies·μL−1; c. Results of RAA-LFD color development with different concentrations of standard plasmids; d. Grayscale analysis results (Corresponding to 6-c), different letters represent significant differences (p<0.05).

    qPCR特异性扩增曲线结果显示,PA的DNA扩增结果为阳性,而其他8株菌的结果为阴性 (图7-a)。RAA-LFD特异性试验结果显示,PA的DNA扩增产物显色较为明显 (T线灰度值为27.0),而其他8株菌显色极不明显 (T线灰度值≤3.0,图7-b—7-c)。普通PCR特异性试验结果显示,PA的DNA扩增后均出现单一条带(740 bp),而其他8株菌DNA均未出现扩增条带 (图7-d)。综上,3种检测方法均出现明显的特异性扩增结果。

    图  7  鳗败血假单胞菌特异性试验结果
    注:a. qPCR特异性扩增曲线;b. 不同菌株的RAA-LFD显色结果;c. 灰度分析结果 (对应7-b),不同字母表示差异显著 (p<0.05);d. 普通PCR特异性电泳图。
    Fig. 7  Results of specificity test for P. anguilliseptica
    Note: a. Specific amplification curves of qPCR; b. RAA-LFD chromogenic results of different strains; c. Grayscale analysis results (Corresponding to 7-b), different letters represent significant differences (p<0.05); d. Specific electrophoresis pattern of PCR.

    通过建立的2种方法及普通PCR方法对采集的10份疑似感染PA的大口黑鲈肾组织样本与30份腹腔注射感染PA的大口黑鲈肾组织样本进行检测。普通PCR方法的阳性样本检出率为62.50%,而qPCR和RAA-LFD方法的阳性样本检出率分别为87.50%和85.00% (表5),较普通PCR方法有明显提高。且qPCR与RAA-LFD方法具有高一致性,较普通PCR方法更可靠。

    表  5  鳗败血假单胞菌3种检测方法的应用结果对比
    Table  5  Comparison of application results of three detection methods for P. anguilliseptica
    检测方法
    Detection method
    样本总数
    Total sample number
    阳性样本检出率
    Detection rate of positive samples
    共同阳性/阴性样本数
    Number of common positive/negative samples
    Kappa值
    Kappa value
    qPCR/RAA-LFD 40/40 87.50%/85.00% 34/5 0.895
    qPCR/PCR 40/40 87.50%/62.50% 25/5 0.385
    RAA-LFD/PCR 40/40 85.00%/62.50% 25/6 0.455
    注:Kappa值的取值范围为0~1.00,其中1.00表示完全一致,0表示完全不一致。常见标准为:Kappa值在0~0.20表示一致性极低,0.21~0.40表示一致性较低,0.41~0.60表示一致性中等,0.61~0.80表示一致性较高,0.81~1.00表示一致性极高。 Note: The Kappa values ranged from 0 to 1.00, with 1.00 indicating complete agreement and 0 indicating complete disagreement. Common criteria are as follows: Kappa of 0−0.20 represents very low consistency; 0.21−0.40 represents low consistency; 0.41−0.60 represents medium consistency; 0.61−0.80 represents high consistency; 0.81−1.00 represents extremely high consistency.
    下载: 导出CSV 
    | 显示表格

    使用qPCR方法测得PA载量 (图8),PA在大口黑鲈肾中载量最高,达3.533×107 拷贝·ng−1,在PA总载量中占比达70.84%;在肝、脾、眼和脑等组织器官中的载量均较高,但在鳃中载量极低。

    图  8  受感染大口黑鲈各部位的鳗败血假单胞菌载量结果
    注:组1、组2、组3表示3个平行的PA感染组,感染量为1×107 CFU·尾−1
    Fig. 8  Results of P. anguilliseptica load in each part of infected M. salmoides
    Note: Group 1, 2 and 3 represent three parallel PA infected groups at 1×107 CFU·ind−1.

    PA曾被认为是海水鱼类的重要病原菌,近年来其感染谱有向淡水鱼类扩展的趋势。Ajengaekanurmaningdyah和Kurniasih[19]于2018年首次报道了PA感染淡水鱼类金鱼 (Cyprinus carpio) 引起死亡;本实验室近年发现PA已成为危害大口黑鲈的新细菌性病原,给大口黑鲈养殖产业造成严重威胁 (另文报道)。PA感染后病鱼处于无明显症状的亚临床期较长,发展到明显症状的临床期往往导致高死亡率,由于PA仅在TSA和BHI等营养丰富的平板上生长[8],且出现菌落的时间超过1周,采用传统病原菌分离鉴定及临床症状进行诊断的难度大、效率低,建立快速、特异、灵敏的PA检测方法在该病的防控中就显得尤为重要。目前,国内外对于PA检测方法的研究较少,除了传统的分离鉴定外,有采用普通PCR进行检测的报道[20-21],但在笔者实验室的应用中发现普通PCR存在假阴性率高的缺点。

    qPCR方法具有重复性好、特异性强和灵敏性高等优点[22]。相较于TaqMan探针法,SYBR Green I法只需在反应混合液中加入引物和待测样品的DNA,无需设计合成探针,简单高效且成本较低[23],目前已应用于鼠伤寒沙门氏菌 (Salmonella typhimurium)[24]、金黄色葡萄球菌 (Staphylococcus aureus)[25]等病原菌的检测。RAA方法对比NASBA、LAMP等温扩增方法有扩增温度较低、仅需1对引物和探针、反应时间较短的优势[26-27],已在铜绿假单胞菌 (Pseudomonas aeruginosa)[28]、鰤诺卡菌 (Nocardia seriolae)[29]等病原菌的检测上得到了应用。RAA-LFD 方法则更利于脱离实验室环境,应用到现场检测中,目前该技术已在副结核分枝杆菌 (Mycobacterium avium subsp. paratuberculosis)[30]、霍乱弧菌 (Vibrio cholerae)[31]、停乳链球菌 (Streptococcus dysgalactiae)[32]等病原菌上开展了研究,为这些病原菌的现场快速检测提供了重要的支持。本研究基于PA管家基因recA成功建立了SYBR Green I real-time quantitative PCR和RAA-LFD 2种快速、特异的检测方法,检出率相较PCR方法提高了至少22.50%。qPCR方法的最低检测浓度为2.816×102 拷贝·μL−1,RAA-LFD方法的最低检测浓度为2.816×104 拷贝·μL−1,尽管建立的qPCR方法的灵敏性高于RAA-LFD方法,但2种方法检测结果的一致性较高,且RAA-LFD方法更具有开展现场检测的优势。这2种方法的建立为PA的早期快速检测与流行病学调查提供了重要的方法依据。

    据报道,PA感染金鱼致肾、肝、脾、肌肉、心、脑等组织器官病变,且肾的病变最为严重,表现为严重坏死[19];PA感染比目鱼 (Psetta maxima) 可致肾、肝、脾、肌肉、脑、眼等组织器官出现病变[33]。由此可见PA感染宿主致全身多组织器官损伤,但对于这些组织器官出现的病变与病原菌的分布及载量之间的关系尚不清楚。本研究建立的SYBR Green I real-time quantitative PCR方法可对PA 进行定量检测,使用该方法对PA感染大口黑鲈后病原菌在各组织的分布与载量进行分析,结果发现PA感染后菌体在肾、肝、脾、眼、脑、心等全身多组织器官分布,其中在肾中的载量最高,其感染致宿主组织器官的病理损伤程度与病原菌的载量有直接关系,表明肾可能是PA感染最重要的侵袭靶器官。

  • 图  1   调查站位示意图

    Figure  1.   Location of survey stations

    图  2   示范区 (a) 和对照区 (b) 鱼类群落丰度-生物量比较曲线

    Figure  2.   Abundance biomass comparison curve of fish community in demonstration zone (a) and control zone (b)

    图  3   银滩南示范区鱼类聚类分析 (a) 和非度量多维标度分析 (b)

    Figure  3.   Cluster analysis (a) and NMDS (b) of fish community in southern station of demonstration zone in Yintan

    图  4   银滩南示范区主要鱼种与环境因子的RDA分析

    注:TEMP. 水温;SAL. 盐度;TRANS. 透明度;DO. 溶解氧;COD. 化学需氧量;TN. 总氮;TP. 总磷;PHYTO. 浮游植物密度;Ala. 侧带天竺鲷;Dma. 蓝圆鲹;Lbe. 细纹鲾;Lbr. 短吻鲾;Llu. 月尾兔头鲀;Lru. 鹿斑仰口鲾;Ped. 二长棘鲷;Pin. 鲬;Sca. 黄斑篮子鱼;Sja. 少鳞鱚;Sne. 魔拟鲉;Sun. 花斑蛇鲻;Tja. 竹荚鱼;Tka. 赤鼻棱鳀。

    Figure  4.   Redundancy analysis between fish species and environmental factors in southern station of demonstration zone in Yintan

    Note: TEMP. Temperature; SAL. Salinity; TRANS. Transparency; DO. Dissolved oxygen; COD. Chemical oxygen demand; TN. Total nitrogen; TP. Total phosphate; PHYTO. Phytoplankton density; Al. A. lateralis; Dma. D. maruadsi; Lbe. L. berbis; Lbr. L. brevirostris; Llu. L. lunaris; Lru. L. ruconius; Ped. P. edita; Pin. P. indicus; Sca. S. canaliculatus; Sja. S. japonica; Sne. S. neglecta; Sun. S. undosquamis; Tja: T. japonicas; Tka. T. kammalensis.

    表  1   示范区和对照区环境因子和渔获参数的均值、最小值和最大值

    Table  1   Mean, minimum and maximum values of environmental factor and catch parameters in demonstration area and control area

    环境因子
    Environment factor
    示范区
    Demonstration zone
    对照区
    Control zone
    显著性
    Significance
    均值 (最小值~最大值)
    Mean (Minimum–Maximum)
    均值 (最小值~最大值)
    Mean (Minimum–Maximum)
    水深 Water depth/m 11.69 (10.70~13.00) 11.87 (11.20~12.70) n.s.
    透明度 Transparency/m 3.76 (3.50~4.20) 3.83 (3.80~3.90) n.s.
    水温 Water temperature/℃ 21.33 (20.90~21.80) 20.47 (20.40~20.60) **
    盐度Salinity/‰ 32.34 (32.09~32.56) 32.35 (32.30~32.41) n.s.
    溶解氧Dissolved oxygen/(mg·L−1) 8.52 (8.45~8.62) 8.29 (8.10~8.50) *
    化学需氧量Chemical oxygen demand/(mg·L−1) 1.69 (0.60~3.00) 2.85 (2.68~2.96) *
    总磷Total phosphate/(mmol·L−1) 0.67 (0.27~3.31) 0.42 (0.35~0.49) n.s.
    总氮Total nitrogen/(mmol·L−1) 38.64 (32.00~40.40) 41.40 (39.00~43.60) n.s.
    浮游植物密度Phytoplankton density/(103 个·m−3) 595 (250~977) 241 (215~266) *
    数量密度 Individual density/(尾·km−2) 19 137 (14 839~32 997) 16 852 (14 439~20 038) n.s.
    生物量密度Biomass density/(kg·km−2) 1 897 (1 504~2 630) 1 244 (1 006~1 493) *
    种类数Number of species 16 (11~21) 12 (11~13) *
    注:根据Mann-Whitney U检验,n.s.. 无显著性差异;*. p<0.05;**. p<0.01。 Note: According to Mann-Whitney U test: n.s.. Insignificant difference; *. p<0.05; **. p<0.01.
    下载: 导出CSV

    表  2   鱼类不同类群组间差异和主要分歧种及其贡献率 (>4%)

    Table  2   Differences in community structure and contribution rates of different species (>4%)

    类群
    Cluster
    平均相异性
    Average dissimilarity/%
    种类
    Species
    数量密度变化
    Change in individual density
    贡献率
    Contribution rate/%
    类群I & III
    Cluster I & III
    53.88 吉打副叶鲹 Alepes djedaba 增加 6.8
    蓝圆鲹 Decapterus maruadsi 减少 6.1
    竹荚鱼 Trachurus japonicus 减少 5.6
    魔拟鲉 Scorpaenopsis neglecta 增加 5.5
    黄斑篮子鱼 Siganus canaliculatus 减少 5.2
    短吻鲾 Leiognathus brevirostris 减少 4.9
    鹿斑仰口鲾 Secutor ruconius 增加 4.4
    纤羊舌鲆 Arnoglossus tenuis 增加 4.4
    细纹鲾 Leiognathus berbis 增加 4.3
    毒拟鲉 Scorpaenopsis diabolus 减少 4.3
    侧带天竺鲷 Apogon lateralis 增加 4.1
    斑鰶 Konosirus punctatus 减少 4.0
    类群III & IV
    Cluster III & IV
    48.85 竹荚鱼 Trachurus japonicus 增加 8.7
    短吻鲾 Leiognathus brevirostris 增加 6.6
    花斑蛇鲻 Saurida undosquamis 增加 5.8
    蓝圆鲹 Decapterus maruadsi 增加 5.4
    魔拟鲉 Scorpaenopsis neglecta 减少 5.2
    吉打副叶鲹 Alepes djedaba 减少 5.1
    Platycephalus indicus 增加 4.9
    细纹鲾 Leiognathus berbis 减少 4.7
    少鳞鱚 Sillago japonica 增加 4.4
    纤羊舌鲆 Arnoglossus tenuis 减少 4.4
    鹿斑仰口鲾 Secutor ruconius 增加 4.1
    赤鼻棱鳀 Thryssa kammalensis 增加 4.0
    下载: 导出CSV

    表  3   RDA分析结果显著性检验

    Table  3   Significance test of redundancy analysis results

    环境因子
    Environment factor
    F p
    水温 Water temperature 3.336 3 0.002**
    化学需氧量 Chemical oxygen demand 2.806 3 0.005**
    溶解氧 Dissolved oxygen 2.173 1 0.032*
    浮游植物密度 Phytoplankton density 2.149 1 0.041*
    总磷 Total phosphate 1.925 8 0.060
    盐度 Salinity 1.325 2 0.225
    透明度 Transparency 1.282 9 0.252
    总氮 Total nitrogen 1.186 3 0.307
    注:**. 在0.01水平上显著相关;*. 在0.05水平上显著相关。 Note: **. Significant correlation at 0.01 level; *. Significant correlation at 0.05 level.
    下载: 导出CSV
  • [1] 陈丕茂, 舒黎明, 袁华荣, 等. 国内外海洋牧场发展历程与定义分类概述[J]. 水产学报, 2019, 43(9): 1981-1869.
    [2] 杨红生, 章守宇, 张秀梅, 等. 中国现代化海洋牧场建设的战略思考[J]. 水产学报, 2019, 43(4): 1255-1262.
    [3] 罗惠桂, 汪佳仪, 谢珍玉, 等. 三亚蜈支洲岛毗邻海域鱼类物种多样性及群落结构特征[J]. 海洋科学, 2023, 47(7): 74-86.
    [4] 陈丕茂, 袁华荣, 贾晓平, 等. 大亚湾杨梅坑人工鱼礁区渔业资源变动初步研究[J]. 南方水产科学, 2013, 9(5): 100-108. doi: 10.3969/j.issn.2095-0780.2013.05.016
    [5]

    HU C, SHUI Y, TIAN K, et al. Functional group classification and niche identification of major fish species in the Qixing Islands Marine Reserve, Zhejiang Province[J]. Biodiv Sci, 2016, 24(2): 175-184. doi: 10.17520/biods.2015225

    [6] 何倩, 刘淑德, 唐衍力, 等. 山东琵琶岛海域人工鱼礁区鱼类群落物种及功能多样性[J]. 中国水产科学, 2023, 30(12): 1479-1495. doi: 10.12264/JFSC2023-0195
    [7] 冯雪, 戴小杰, 袁华荣, 等. 外伶仃海洋牧场附近海域游泳生物群落结构及其与环境因子的关系[J]. 上海海洋大学学报, 2024, 33(1): 186-201. doi: 10.12024/jsou.20230304115
    [8] 张亚洲, 蒋日进, 梁君. 嵊泗马鞍列岛海域国家级海洋牧场渔业资源增殖养护效果评析[J]. 浙江海洋大学学报 (自然科学版), 2022, 41(5): 466-472.
    [9] 李寇军, 邱永松, 王跃中, 等. 自然环境变动对北部湾渔业资源的影响[J]. 南方水产, 2007, 3(1): 7-13.
    [10] 陈国宝, 李永振, 陈丕茂, 等. 南海北部陆架区海域蓝圆产卵场的研究[J]. 热带海洋学报, 2003, 22(6): 22-28. doi: 10.3969/j.issn.1009-5470.2003.06.004
    [11] 张公俊, 杨长平, 孙典荣, 等. 北部湾中北部海域鱼类群落的季节变化特征[J]. 南方农业学报, 2021, 52(10): 2861-2871. doi: 10.3969/j.issn.2095-1191.2021.10.027
    [12] 张曼, 王雪辉, 蔡研聪, 等. 北部湾带鱼空间聚散变化特征[J]. 中国水产科学, 2022, 29(11): 1647-1658. doi: 10.12264/JFSC2021-0609
    [13] 侯刚, 冯钰婷 , 陈妍颖, 等. 北部湾二长棘犁齿鲷时空分布及其与环境因子的关系[J]. 广东海洋大学学报, 2021, 41(4): 8-16.
    [14] 黎树式, 戴志军, 葛振鹏, 等. 北部湾北部生态环境灾害变化研究[J]. 灾害学, 2014, 29(4): 43-47. doi: 10.3969/j.issn.1000-811X.2014.04.009
    [15]

    JACCARD P. The distribution of flora in the alpine zone[J]. New Phytol, 1912, 11: 37-50. doi: 10.1111/j.1469-8137.1912.tb05611.x

    [16]

    WARWICK R M. A new method for detecting pollution effects on marine macrobenthic communities[J]. Mar Biol, 1986, 92: 557-562. doi: 10.1007/BF00392515

    [17] 裴精花, 陈清华, 范金金, 等. 珠江口海域游泳动物群落结构及多样性特征[J]. 南方农业学报, 2023, 54(12): 3727-3738. doi: 10.3969/j.issn.2095-1191.2023.12.027
    [18] 王雪辉, 林昭进, 王跃中, 等. 南海西北部陆架区鱼类的种类组成与群落格局[J]. 生态学报, 2013, 33(7): 2225-2235.
    [19] 郭彪, 于莹, 张博伦, 等. 天津大神堂海域人工鱼礁区游泳动物群落特征变化[J]. 海洋渔业, 2015, 37(5): 409-418. doi: 10.3969/j.issn.1004-2490.2015.05.003
    [20] 徐浩, 曾晓起, 顾炎斌, 等. 人工鱼礁对山东莱州朱旺港海区游泳动物的群落结构及季节变化的影响[J]. 中国海洋大学学报, 2012, 42(5): 47-54.
    [21] 罗峥力, 杨长平, 王良明, 等. 北部湾北部沿岸海域鱼类资源时空分布特征及多样性变化[J]. 南方农业学报, 2023, 54(6): 1847-1857. doi: 10.3969/j.issn.2095-1191.2023.06.026
    [22] 冯雪, 范江涛, 孙晓, 等. 珠海外伶仃人工鱼礁对鱼类资源养护效果初步评估[J]. 南方农业学报, 2021, 52(12): 3228-3236. doi: 10.3969/j.issn.2095-1191.2021.12.005
    [23] 印瑞, 周永东, 梁君, 等. 中街山列岛海洋牧场大黄鱼时空分布与环境因子的关系[J]. 浙江海洋大学学报 (自然科学版), 2022, 41(6): 483-489.
    [24] 韩业越, 殷蕊, 孙桂清, 等. 北戴河国家级海洋牧场示范区人工鱼礁建设效果评价[J]. 河北农业, 2020, 316(4): 24-28.
    [25] 王欢欢, 毕福洋, 曹敏, 等. 獐子岛海洋牧场秋季渔业资源声学调查与评估[J]. 大连海洋大学学报, 2018, 33(6): 802-807.
    [26] 崔晨, 张云岭, 张秀文, 等. 唐山祥云湾海洋牧场渔业资源增殖效果评估[J]. 河北农业, 2021(1): 25-31.
    [27] 李建都, 赵祺, 刘晋冀, 等. 黄渤海不同人工鱼礁区渔业生物群落结构特征及生物增量影响要素[J]. 中国水产科学, 2023, 30(3): 371-383. doi: 10.12264/JFSC2022-0406
    [28]

    HYKEMA A, DEBROT A O, OSINGA R, et al. Fish assemblages of three common artificial reef designs during early colonization[J]. Ecol Eng, 2020, 157: 105994.

    [29] 王宇, 侯纯强, 汪笑宇, 等. 天津近岸海域人工鱼礁区渔业资源增殖效果初步评估[J]. 海洋湖沼通报, 2018(2): 137-145.
    [30] 张荣良, 刘辉, 孙东洋, 等. 烟台近岸人工鱼礁与自然岩礁底层渔业生物群落特征对比分析[J]. 海洋与湖沼, 2021, 52(3): 697-707. doi: 10.11693/hyhz20201000296
    [31] 王新萌, 唐衍力, 孙晓梅, 等. 威海小石岛人工鱼礁海域渔获物群落结构特征及其与环境因子相关性[J]. 海洋科学, 2016, 40(11): 34-39. doi: 10.11759/hykx20151123001
    [32]

    ROSAS-ALAYOLA J, HEMANDES-HERRERA A, GALVAN-MAGANA F, et al. Diet composition of sailfish (Istiophorus platypterus) from the southern Gulf of California, Mexico[J]. Fish Res, 2002, 57: 185-195. doi: 10.1016/S0165-7836(01)00344-7

    [33]

    KADYE W T, MAGADZA C H D, MOYO N A G, et al. Stream fish assemblages in relation to environmental factors on a montane plateau (Nyika, Malawi)[J]. Enriron Biol Fish, 2008, 83(4): 417-428. doi: 10.1007/s10641-008-9364-4

    [34]

    SANDER R M J, CHANDLER C R, LANDRY A M J. Hydrological, diel and lunar factors affecting fishes on artificial reefs off Panama City, Florida[J]. Bull Mar Sci, 1985, 37(1): 318-328.

    [35]

    OROPEZE O A, BALART E F. Community structure of reef fish in several habitats of a rocky reef in the Gulf of California[J]. Mar Ecol, 2001, 22(4): 283-305.

    [36] 薛燕妮, 王洁英, 别又才, 等. 清远市南部淡水养殖鱼塘水质抽样调查与污染评价[J]. 黑龙江环境通报, 2024, 37(4): 40-42. doi: 10.3969/j.issn.1674-263X.2024.04.014
    [37] 杨斌, 钟秋平, 鲁栋梁, 等. 钦州湾海域COD时空分布及对富营养化贡献分析[J]. 海洋科学, 2014, 38(3): 20-25.
    [38] 刘子萌, 冯民权, 杨锐婧, 等. 漳泽水库浮游植物群落结构特征及其影响因子[J]. 黑龙江大学工程学报, 2024, 15(1): 76-84.
    [39] 王奕杨, 毛志刚, 谷孝鸿, 等. 高宝邵伯湖鱼类群落结构及其多样性特征[J/OL]. 生态学杂志. https://link.cnki.net/urlid/ 21.1148.Q.20240205.1147.002
    [40]

    MAES J, STEVENS M, BREINE J. Modelling the migration opportunities of diadromous fish species along a gradient of dissolved oxygen concentration in a European tidal watershed[J]. Estuar Coast Shelf Sci, 2007, 75(1/2): 151-162.

    [41] 张雷敏, 梁化亮, 米海峰, 等. 低氧胁迫对鱼类的影响以及营养调控策略的研究进展[J]. 动物营养学报, 2023, 35(10): 6207-6219. doi: 10.12418/CJAN2023.570
    [42] 曾嘉维, 林坤, 王学锋, 等. 雷州湾附近海域鱼类群落结构及其与环境因子的关系[J]. 中国水产科学, 2019, 26(1): 108-117.
    [43] 曾雷, 陈国宝, 李纯厚, 等. 大亚湾湾口游泳生物群落季节异质特征与生态效应分析[J]. 南方水产科学, 2019, 15(3): 22-32. doi: 10.12131/20180246
    [44] 胡翠林, 张洪亮, 张亚洲, 等. 带鱼保护区春秋季鱼类群落特征及与环境因子的关系[J]. 水产学报, 2018, 42(5): 694-703.
    [45] 唐广隆, 刘永, 吴鹏, 等. 珠江口万山群岛海域春季渔业资源群落结构特征及其与环境因子的关系[J]. 中国水产科学, 2022, 29(8): 1198-1209.
    [46] 张雪, 徐晓甫, 戴媛媛, 等. 天津近岸人工鱼礁海域浮游植物群落及其变化特征[J]. 渔业科学进展, 2018, 39(6): 1-10.
    [47] 谢斌, 张硕, 李莉, 等. 海州湾海洋牧场浮游植物群落结构特征及其与水质参数的关系[J]. 环境科学学报, 2017, 37(1): 121-129.
    [48] 刘懂, 陈晨, 王莉, 等. 象山港海洋牧场示范区浮游植物落特征及其与环境因子的关系[J]. 海洋与湖沼, 2016, 47(5): 1024-1032. doi: 10.11693/hyhz20160500096
推荐阅读
基于环境dna技术的西江珍稀鱼类省级自然保护区鱼类多样性研究
钟占友 et al., 南方水产科学, 2025
三峡库区夏季和秋季鱼类群落结构及环境影响因子研究
吴凡 et al., 南方水产科学, 2025
雷州半岛东部近岸海域春秋季浮游植物群落结构特征及环境影响因子分析
侯迪遥 et al., 南方水产科学, 2025
广西银滩南部海域海洋牧场渔业资源评估
牛麓连 et al., 南方水产科学, 2024
黄海北部近岸海域渔业生物群落结构
崔培东 et al., 水产学报, 2024
龙羊峡水库鱼类群落结构变化
刘旭 et al., 中国水产科学, 2024
The seas around china in a warming climate
Wang, Fan et al., NATURE REVIEWS EARTH & ENVIRONMENT, 2023
The influencing factors and mechanisms for urban flood resilience in china: from the perspective of social-economic-natural complex ecosystem
Zhu, Shiyao et al., ECOLOGICAL INDICATORS, 2023
Dependence of daily precipitation and wind speed over coastal areas: evidence from china's coastline
HYDROLOGY RESEARCH, 2023
Exploring water quality as a determinant of the existence of soft shell crab (scylla olivacea) in different habitats of the large marine ecosystem in the bay of bengal.
WATER CONSERVATION AND MANAGEMENT, 2024
Powered by
图(4)  /  表(3)
计量
  • 文章访问数:  1247
  • HTML全文浏览量:  39
  • PDF下载量:  94
  • 被引次数: 0
出版历程
  • 收稿日期:  2024-05-09
  • 修回日期:  2024-05-29
  • 录用日期:  2024-06-26
  • 网络出版日期:  2024-06-28
  • 刊出日期:  2024-10-04

目录

/

返回文章
返回