Study on monitoring technology of fish motion trajectories based on Acoustic Tag System
-
摘要:
声学标签监测技术是声波监测技术的一种主动识别方法,其通过接收、处理移植或捆绑于鱼类身上的声学标签发射的声波信号来获得鱼类二维、三维运动轨迹,具有原位观察、定位准确、数据处理简单和数据连续性好等优点。文章介绍了声学标签监测技术的组成及工作原理,并根据异常数据的不同特点总结了一套完整的数据去噪、清洗的处理方法,该法可快速得到鱼类的二维、三维运动轨迹。通过实例分析了鱼类运动轨迹的分布情况以及变化规律。根据鱼类在不同水生态环境中的运动行为响应情况,可以直接或间接判断水中污染物或有毒物质对鱼类的潜在影响以及实际毒性,为渔业养殖水质环境评价、水质监测及预警、水生态健康评价以及水生态修复效果评价的研究提供参考依据。
Abstract:Acoustic tag monitoring technology, an active identification method of acoustic detection technology, can obtain 2D and 3D motion trajectories of fish by receiving and processing the acoustic signals emitting from acoustic tags which are transplanted or tied to fish. It has the advantages of in-situ observation, accurate positioning, simple data processing and good data continuity. The paper introduces the composition and working principle of acoustic tag monitoring technology. A complete data processing method of denoising and cleaning is proposed according to the different characteristics of abnormal data. With this method, the real-time 2D and 3D motion trajectories can be quickly obtained. The distribution of fish motion trajectories and changing rules are analyzed by an example. According to the behavioral response of fish in different water ecological environments, the potential effects and actual toxicity of pollutants or toxic substances in water can be judged directly or indirectly. It provides the references for water quality assessment in aquaculture, monitoring and early warning of water environment, health evaluation of aquatic ecological, as well as evaluation of aquatic ecological restoration.
-
鱼类在水体中的生命活动会受到不同环境的影响,水中的溶解氧浓度变化对鱼类有至关重要的影响。水中溶解氧质量浓度在4.0 mg·L−1以上时,鱼类可以正常生长发育;低于1.0 mg·L−1时,大部分鱼类会出现浮头现象[1]。水中溶解氧浓度降低,不仅使鱼类呼吸和摄氧能力下降[2],还会影响鱼类细胞的存活和信号传递[3],并直接影响其产卵、交配、生长、发育等一系列生命活动[4],严重时还会出现死亡,破坏种群内部的动态平衡。
斑马鱼 (Danio rerio) 具有易繁殖、发育快等优点[5],一直作为模式生物用于生物医学研究[6-7]。作为模式生物,斑马鱼的应激调节功能、心血管系统功能等与人类高度相似[8-9]。鳃是鱼类重要的黏膜免疫器官和呼吸器官,在鱼类低氧适应研究中常作为主要的研究对象,如低氧胁迫下鳃组织中热休克蛋白的研究[10],草鱼 (Ctenopharyngodon idellus) 在低氧胁迫下鳃的差异蛋白质组学及热休克诱导[11],低氧胁迫对卵形鲳鲹 (Trachinotus ovatus) 鱼体鳃器官的影响[12]等。
miRNAs可以通过调控其靶基因的表达水平从而参与细胞的各类进程,并参与许多关键的生理进程与病理过程[13]。已有多项研究表明,miRNAs是植物和动物面临环境胁迫作出响应的关键调节剂,一些miRNAs可以通过调控基因表达,恢复或重建新的表达程序,从而增强细胞对胁迫的耐受性[14]。低氧是生物常常面临的一种胁迫,目前已有一些应对低氧胁迫miRNAs的报道。如缺氧性神经胶质瘤来源的外泌体通过信号转导和转录活化因子3 (STAT3) 和核因子κB (NF-κB) 途径靶向端粒重复结合因子2 (terf2ip),传递microRNA-1246诱导M2巨噬细胞极化[15]。miR-204可以作用于血管内皮生长因子 (vegf) 的3'-UTR区,是一个内源性的vegf表达调控因子。miR-204通过基因网络调控应答低氧胁迫[16]。miR-126-5p可以作为一种新型的miRNA,在缺氧条件下靶向白细胞介素17 (IL-17A) 来调节大鼠心肌细胞 (H9c2) 的活力和凋亡[17]。
热休克蛋白 (Heat shock proteins, HSPs) 是生物在面对环境中的物理、化学、生物等刺激发生应激反应后大量产生的,常被称为应激蛋白[18]。热休克蛋白是生物体内最古老的分子之一。生物体为抵御环境变化所带来的刺激,会减少其他正常蛋白的合成,同时增加HSP的合成以应对环境的变化,帮助生物体恢复正常[10]。因此本文在对常氧和低氧胁迫下的斑马鱼鳃组织进行小RNA组测序的基础上,进一步筛选可能与低氧胁迫相关的micRNAs,并用其对热休克蛋白基因进行了靶基因预测,以期进一步挖掘斑马鱼的低氧适应机制。
1. 材料与方法
1.1 材料
1.1.1 实验试剂
small RNA分离试剂盒,DEPC水 (生物生工有限公司),氯仿 (24∶1),Trizol 试剂,异丙醇,无水乙醇 (吉泰生物公司)。
1.1.2 实验器材
研磨棒、EP管、超净工作台、冰块、移液器、各类枪头、剪刀、镊子、锡纸均购买于上海生物生工有限公司,低氧驯化箱 (长沙华晓电子科技有限公司定制),离心机,旋涡仪 (德国Eppendorf有限公司)。
1.2 方法
1.2.1 样本低氧处理
将本实验室培养的多代斑马鱼分别置于1.0 mg·L−1的低氧驯化箱中,保持该浓度进行低氧胁迫,同时保留6.7 mg·L−1的溶解氧质量浓度作为常氧对照组。在使用1.0 mg·L−1的溶解氧质量浓度进行低氧胁迫时,发现2周后斑马鱼不再出现浮头的现象,推测经2周低氧胁迫后其逐渐适应了低氧环境。故采用低氧驯化2周后的斑马鱼与常氧下的斑马鱼进行比对。低氧胁迫2周后,取出常氧/低氧条件下的斑马鱼,每个溶解氧质量浓度下取样本15尾 (体长3~4 cm),解剖取出两组样本的鳃组织,备用。
1.2.2 高通量测序和生信分析
提取RNA样本,按照NucleoZOL试剂盒进行。将所提取的RNA进行电泳检测。检测后将EP管放入−80 ℃冰箱进行保存。将上述常氧和低氧下斑马鱼鳃组织提取后检测合格的RNA[19],按照small RNA分离试剂盒的方法,分别进行small RNA的分离。将提取的small RNA进行纯化,运用荧光光度计进行定量,上机检测测序得到原始测序结果。为保证数据的质量,对原始数据进行处理。去除低质量的reads (N比例大于10%的reads、有5'接头污染的reads、无3'接头序列和插入片段的reads、3'接头序列以及polyA/T/G/C的reads) 后得到的clean small RNA reads数。其中大多数小RNA的长度为21~23 nt。对于该物种的ncRNA注释,若有该物种小RNA的注释信息,就用该物种ncRNA注释所测的small RNA。若没有该物种的信息,则选择Rfam数据库中rRNA、tRNA、snRNA和 snoRNA来注释测序所得的small RNA。
1.2.3 靶基因预测
本实验室前期通过对低氧、常氧条件下鳃组织的转录组比较分析发现,常氧低氧条件下斑马鱼鳃中一共筛选获得28个显著差异表达的热休克蛋白基因,包括表达量显著下调基因12个,显著上调基因16个[12]。对miRNAs 测序和斑马鱼鳃转录组进行关联分析。针对低氧胁迫和常氧条件下斑马鱼鳃中显著差异表达的miRNAs,对实验室前期筛选获得的28个差异热休克蛋白基因进行靶基因的预测分析。miRNAs的靶基因预测使用TargetScanFish (http://www.targetscan.org/fish_62/)、miRanda (http://www.microrna.org/microrna/home.do) 2个网站同时进行。
2. 结果
2.1 低氧与常氧下斑马鱼鳃的小RNA测序
常氧和低氧斑马鱼鳃样本的小RNA组测序分别产生6 995 009和6 662 504 bp的原始数据。去除低质量的数据后分别得到6 585 748和5 941 304 bp的clean data。
常氧与低氧小RNA序列比对后获得了相应结果,获得饼状图 (图1)。根据饼图比例分析,可以看出,其中常氧的miRNA约占整个small RNA总数量的40%。而低氧的约占20%。提示与常氧相比,低氧胁迫下的miRNA有降低的趋势 (图1)。
2.2 低氧/常氧下鳃组织的差异miRNAs分析
根据测序结果,首先排除tRNA、rRNA等小分子的RNA。利用归一化法对比低氧胁迫斑马鱼与常氧斑马鱼中同一个miRNA的表达差异量。利用火山图呈现差异miRNAs的整体分布情况。结果显示,低氧/常氧条件下的斑马鱼鳃间一共筛选出差异表达的miRNAs共32个 (图2)。使用校正后的显著水平 (P) 和差异倍数 (Fold change) 2个水平进行评估,设置显著差异表达miRNAs的筛选条件为P<0.01和 |log2(fold change)|>1。
针对32个差异miRNAs,同时使用 |log2FC|≥1,P<0.05,且表达量≥50作为临界值,从低氧和常氧的比较中,一共鉴定获得低氧与常氧条件下显著差异表达的15个miRNAs。其中,13个miRNAs在低氧胁迫斑马鱼鳃中的表达量显著上调、2个miRNAs (miR-455-3p、miR-125b-5p) 的表达量显著下调。图3为显著差异表达miRNAs的聚类。
2.3 差异miRNAs靶向热休克蛋白基因的预测结果
针对前期筛选获得的低氧胁迫与常氧条件下显著差异表达的28个热休克蛋白基因 (包括12个表达量显著下调的hsp基因、16个显著上调的hsp基因)[11]进行靶基因预测,结果显示,7个显著差异表达的miRNAs可以靶向9个热休克蛋白基因。图4显示出单个miRNA靶向热休克蛋白基因的数量。
表达量显著下调的miR-455-3p可以靶向2个显著上调的热休克蛋白基因 (表1)。表达量显著上调的miRNAs (dre-miR-194a、dre-miR-155、dre-miR-130c、dre-miR-9、dre-miR-29a、dre-miR-96-5p) 可以靶向7个显著下调的热休克蛋白基因 (表2)。
表 1 斑马鱼低氧与常氧鳃中下调的miRNAs靶基因预测Table 1. Prediction of down-regulated miRNAs target genes in hypoxic and normoxic gills of D. reriomiRNA名称
miRNA name靶基因
Gene binding结合位点
Site目标区域的预测配对
Predicted pairing of target regiondre-miR-455-3p hspa14 935—941 dnajb6b 1157—1163 1829—1835 表 2 斑马鱼低氧与常氧鳃中上调miRNAs靶基因预测Table 2. Prediction of up-regulated miRNAs target genes in hypoxic and normoxic gills of D. reriomiRNA名称
miRNA name靶基因
Gene binding结合位点
Site目标区域的预测配对
Predicted pairing of target regiondre-miR-194a hspa12a 2262—2268 dnajc5aa 4114—4120 hspb7 1644—1650 hsp70.3 245—251 dnajb2 3834—3840 dre-miR-155 hspa12a 2771—2777 hspg2 3611—3617 hspa13 308—314 903—910 dnajb2 1484—1490 3420—3426 dre-miR-130c hspa12a 3219—3225 4709—4715 dnajb2 920—926 dre-miR-9 dnajc5aa 2348—2354 dnajb2 5033—5039 dre-miR-29a hspg2 2832—2839 dre-miR-96-5p dnajb2 2290—2296 5045—5051 5094—5100 2.4 差异miRNAs靶向的热休克蛋白基因富集分析
针对差异miRNAs靶向的热休克蛋白基因进行富集分析发现,富集到的生物过程功能前6条通路主要与发育相关,包括dnajb6b和hspg2两个靶基因 (图5)。低氧胁迫的斑马鱼鳃中表达量显著上调的miR-455-3p靶向dnajb6b,而hspg2同时受到2个上调的miRNAs (miR-155和miR-29a) 调控。差异miRNAs靶向的热休克蛋白基因富集的前20条通路主要涉及的基因除了上面提及的2个miRNAs外,还包括miR-194a和miR-130c同时靶向的hspb7。
3. 讨论
本研究发现了在缺氧反应中差异表达的15个miRNAs,13个miRNAs在低氧胁迫的斑马鱼鳃中表达量上调,2个miRNAs的表达显著下调。也有研究表明,在下调的miRNAs中,miR-125b-5p的靶基因rps3a通过在翻译机制中发挥调节作用以抑制细胞凋亡[20];而在上调的miRNAs中,mir-192可以通过靶向E盒结合锌指蛋白基因 (Zeb2) 来保护肝脏免受氧化应激诱导的损伤,增强肝脏的低氧耐受性[21-22]。在人类细胞中,miR-29b的上调可以靶向TNF受体相关因子5 (TRAF5) 保护心肌细胞免受缺氧诱导的细胞凋亡[23]。miR-29b在调节细胞凋亡中显示出重要作用[24]。另外,小鼠中miR-216b可以作用于自噬相关蛋白13基因 (Atg13),且使缺氧条件下细胞的自噬减少,并减少细胞的凋亡[25]。
高海拔人群中的miR-210-3p水平与红细胞计数以及血红蛋白和血细胞比容显示出强正相关性,被认为是人类适应高海拔地区生活的重要miRNA[26]。miR-194过表达可以保护缺氧诱导的人肾皮质近曲小管上皮细胞 (HK-2) 损伤[27]。miR-155被发现在低氧条件下促进内皮细胞的血管生成[28]。人参皂甙 (GS-Rb1) 可以增加mir-29a的表达量,保护缺氧的心肌细胞[29]。大鼠中Hif-1α诱导的miR-9上调在缺氧期间有助于肺动脉平滑肌细胞的表型调节[30],miR-96-5p已知有抑制细胞凋亡的功能[31]。miR-1可能在转录后水平直接或间接调控HSP90aa1和HSP90b1。过表达miR-1后缺氧复氧的HSP90蛋白及其亚型90aa1和90b1表达水平更低,结合之前的结果提示miR-1可能在心肌缺氧复氧中调控HSP90[32]。在高血压心肌肥厚的早期代偿阶段,心肌miR-378表达的下降使其对内源性HSF1转录后抑制作用减弱,进而对HSF1的代偿性升高发挥了重要的调控作用[33]。可见,miRNAs在低氧条件下在其他器官起着抑制细胞凋亡,增强器官的低氧耐受性,保护细胞免受缺氧带来的损伤等一系列功能。因此,本研究筛选出的低氧和常氧之间差异表达的miRNAs很可能在低氧适应机制中起重要作用。已知生物体在面对环境变化时会通过改变蛋白或者调节mRNAs的翻译以适应环境,本实验室已经发表过低氧胁迫下鳃组织相关热休克差异蛋白基因,验证了一些热休克蛋白基因应对低氧胁迫的作用[11]。因此,本研究利用筛选到的斑马鱼低氧与常氧状态下差异表达的15条miRNAs对28个热休克蛋白基因[11]进行靶基因预测,再结合miRNAs与预测的mRNAs的表达呈负相关这一特性,对预测的靶基因进行分析。
热休克蛋白被认为与正常和异常的胚胎发育密切相关。低氧胁迫下,低表达的miR-455-3p通过同时靶向提高hspa14和dnajb6b的表达,有可能增强了生物体的发育和机体保护,进而增强对低氧的适应。本研究发现,miR-194a同时靶向5个热休克蛋白基因 (hspa12a、dnajc5aa、hspb7、hsp70.3、dnajb2);而miR-155可以同时靶向4个热休克蛋白基因 (hspa12a、hspg2、hspa13、dnajb2)。本研究还发现,热休克蛋白基因dnajb2同时受到5个miRNAs调控。基因dnajb2是HSP70的伴侣调节剂,主要在神经系统中表达[34]。等距遗传性运动神经病 (dHMN) 是一组罕见的遗传性神经肌肉疾病,其特征是在没有感觉症状的情况下会影响腓骨肌肉的萎缩,dnajb2是23个与dHMN有关的基因之一,主要从dnajb2起始[35]。可以推测,5个miRNAs靶向抑制dnajb2基因的表达,有可能减轻了对低氧环境下生物体神经系统的伤害。
Hspa14可能是肢体发育的相关基因[36],在成年斑马鱼中进行无偏见的诱变筛选,确定了dnajb6b是心肌病的新型遗传修饰剂[37]。缺乏热休克蛋白基因hspg2可减轻在低氧中引起的动脉高压[38]。另外,在靶向的热休克蛋白基因GO功能富集前20条通路中,均有hspg2的参与。miR-155和miR-29a同时靶向hspg2,很可能通过抑制hspg2的表达来降低低氧胁迫下引起的动脉高压,增强对低氧环境的适应。
综上,本研究结合常氧与低氧下差异表达的miRNAs和热休克蛋白基因的关联分析,为探究鱼类低氧适应机制提供了新的研究思路。
-
-
[1] KANG I J, MOROISHI J, NAKAMURA A, et al. Biological monitoring for detection of toxic chemicals in water by the swimming behavior of small freshwater fish[J]. J Fac Agr Kyushu U, 2009, 54(1): 209-214.
[2] 王丽娜, 方景龙. 基于视频理解的活体鱼水质监测实验[J]. 杭州电子科技大学学报, 2010, 30(6): 36-38. doi: 10.3969/j.issn.1001-9146.2010.06.010 [3] 徐盼麟, 韩军. 基于视频的水箱内鱼类运动行为观测系统[J]. 湖南农业科学, 2011(19): 140-142. doi: 10.3969/j.issn.1006-060X.2011.19.043 [4] 江丹丹, 桂福坤. 基于视频图像的鱼类运动行为轨迹追踪[J]. 浙江海洋学院学报(自然科学版), 2015, 34(2): 112-118. [5] 范伟康. 基于鱼群行为分析的水质监测系统研究[D]. 杭州: 浙江工业大学, 2015: 13-14. [6] 余祥瑞. 基于双摄像机视频的生物水质预警系统的研究[D]. 宁波: 宁波大学, 2017: 5-6. [7] 颜鹏东, 谭均军, 高柱, 等. 基于视频跟踪的竖缝式鱼道内鱼类运动行为分析[J]. 水生生物学报, 2018, 42(2): 250-254. doi: 10.7541/2018.031 [8] 于铭, 杨士莪. 被动式渔探仪开发设计与试验研究[J]. 应用声学, 2007, 26(4): 231-238. doi: 10.3969/j.issn.1000-310X.2007.04.008 [9] 于铭, 杨士莪, 牟冬英. 被动式渔探仪的信号检测与方位估计[J]. 哈尔滨工业大学学报, 2008, 40(1): 147-151. doi: 10.3321/j.issn:0367-6234.2008.01.037 [10] STANTON, TIMOTHY K. 30 years of advances in active bioacoustics: a personal perspective[J]. Meth Oceanogr, 2012, s1/2: 49-77.
[11] GODLEWSKA M, SWIERZOWSKI A. Hydroacoustical parameters of fish in reservoirs with contrasting levels of eutrophication[J]. Aquat Living Resour, 2003, 16(3): 167-173. doi: 10.1016/S0990-7440(03)00014-7
[12] EVERSON I, TAABU M A, KAYANDA R. Acoustic estimates of commercial fish species in Lake Victoria: moving towards ecosystem-based fisheries management[J]. Fish Res, 2013, 139: 65-75. doi: 10.1016/j.fishres.2012.09.019
[13] 武智, 谭细畅, 李新辉, 等. 珠江首次禁渔西江段鱼类资源声学跟踪监测分析[J]. 南方水产科学, 2014, 10(3): 24-28. doi: 10.3969/j.issn.2095-0780.2014.03.004 [14] 武智, 李新辉, 李捷, 等. 红水河岩滩水库鱼类资源声学评估[J]. 南方水产科学, 2017, 13(3): 20-25. doi: 10.3969/j.issn.2095-0780.2017.03.003 [15] SZCZUCKA J, HOPPE T, SCHMIDT B, et al. Acoustical estimation of fish distribution and abundance in two Spitsbergen fjords[J]. Oceanologia, 2017, 59(4): 585-591.
[16] 张翔, 沈蔚, 童剑锋, 等. 基于DIDSON双频识别声纳技术的青草沙水库鱼类资源量评估[J]. 上海海洋大学学报, 2017, 26(4): 561-569. doi: 10.12024/jsou.20161201919 [17] CONTI S G, MAURER B D, ROUX P, et al. Acoustical monitoring of fish behavior in a tank[J]. J Acoust Soc Am, 2004, 116(4): 2489-2489.
[18] 张赞, 蔺丹清, 汤勇, 等. 基于声学测量方法的大伙房水库鱼类资源季节变动特征[J]. 南方水产科学, 2014, 10(6): 12-19. doi: 10.3969/j.issn.2095-0780.2014.06.002 [19] 孙明波, 谷孝鸿, 曾庆飞, 等. 基于水声学方法的天目湖鱼类季节和昼夜空间分布研究[J]. 生态学报, 2015, 23(9): 5597-5605. [20] 连玉喜, 叶少文, 黄耿, 等. 三峡库区香溪河鱼类资源水声学探测效果的昼夜差异研究[J]. 水生生物学报, 2015, 39(5): 1041-1045. doi: 10.7541/2015.136 [21] 李斌, 陈国宝, 于杰, 等. 海南陵水湾口海域不同季节鱼类资源声学探查[J]. 水产学报, 2018, 42(4): 544-555. [22] 郭杰. 航道整治透水框架群对鱼类集群影响的水声学探测[J]. 水生态学杂志, 2015, 36(5): 29-35. [23] 马燕芹, 司纪锋. 基于水声技术的黄海近海鱼类活动定点监测研究[J]. 渔业现代化, 2016, 43(4): 70-75. doi: 10.3969/j.issn.1007-9580.2016.04.013 [24] EHRENBERG J E, STEIG T W. A study of the relationship between tag-signal characteristics and achievable performances in acoustic fish-tag studies[J]. ICES J Mar Sci, 2009, 66(6): 1278-1283. doi: 10.1093/icesjms/fsp051
[25] 王志超, 陈国宝, 曾雷. 基于声学标志和无线跟踪方法的鱼类行为研究[J]. 南方水产科学, 2018, 14(2): 51-59. doi: 10.3969/j.issn.2095-0780.2018.02.007 [26] STEIG T W, TIMKO M A. Innovative techniques for analyzing the three-dimensional behavioral results from acoustically tagged fish[J]. J Acoust Soc Am, 2005, 117(4): 2367-2367.
[27] SEMMENS B X. Acoustically derived fine-scale behaviors of juvenile Chinook salmon (Oncorhynchus tshawytscha) associated with intertidal benthic habitats in an estuary[J]. Can J Fish Aquat Sci, 2008, 65(9): 2053-2062. doi: 10.1139/F08-107
[28] SKALSKI J R, STEIG T W, HEMSTROM S L. Assessing compliance with fish survival standards: a case study at Rock Island dam, Washington[J]. Environ Sci Policy, 2012, 18: 45-51. doi: 10.1016/j.envsci.2012.01.001
[29] 李荣辉.鱼类栖息迁徒习性及其监测技术研究[D]. 南宁: 广西大学, 2013: 75-94 [30] ABAID N, PORFIRI M. Collective behavior of fish shoals in one-dimensional annular domains[C]. Riverside: Institute of Electrical and Electronics Engineers, 2010 American Control Conference, Riverside: Institute of Electrical and Electronics Engineers,2010: 63-68.
[31] SCHULTZ A A, KUMAGAI K K, BRIDGES B B. Methods to evaluate gut evacuation rates and predation using acoustic telemetry in the tracy fish collection facility primary channel[J]. Anim Biotel, 2015, 3(1): 1-9. doi: 10.1186/s40317-014-0021-8
[32] HOLLO T, WATSON B M, JOHNSTON S V, et al. Behaviour of growth hormone transgenic coho salmon Oncorhynchus kisutch in marine mesocosms assessed by acoustic tag telemetry[J]. J Fish Biol, 2017, 90(4): 1660-1667. doi: 10.1111/jfb.13233
[33] ROMINE J G, PERRY R W, JOHNSTON S V, et al. Identifying when tagged fishes have been consumed by piscivorous predators: application of multivariate mixture models to movement parameters of telemetered fishes[J]. Anim Biotel, 2014, 2(1): 3. doi: 10.1186/2050-3385-2-3
[34] 黄月群, 蔡德所, 李明泉, 等. 大坝泄水对鱼类洄游能力的影响研究[J]. 人民长江, 2019, 50(8): 74-80. [35] CAI D S, LI R H. Research on fish habitat based on acoustic fish tracking system[J]. JCIT, 2012, 7(23): 195-201. doi: 10.4156/jcit.vol7.issue23.23
[36] HUANG Y Q, CAI D S, LI M Q, et al. Influence of changes in dissolved oxygen content on fish behavioral trajectories during water eutrophication[J]. Appl Ecol Env Res, 2018, 17(1): 653-666.
[37] 赵方旭. 不同水力条件下鱼类运动行为轨迹特性研究[D]. 南宁: 广西大学, 2016: 48-55 [38] KRISTIANSEN T S, FEMO A, HOLM J C, et al. Swimming behavior as an indicator of low growth rate and impaired welfare in Atlantic halibut (Hippoglossus hippoglossus L.) reared at three stocking densities[J]. Aquaculture, 2004, 230(1/2/3/4): 137-151.
[39] MCFARLANE W J, CUBITT K F, WILLIAMS H, et al. Can feeding status and stress level be assessed by analyzing patterns of muscle activity in free swimming rainbow trout (Oncorhynchus mykiss Walbaum)[J]. Aquaculture, 2004, 239(4): 467-484.