2. 农业农村部淡水渔业健康养殖重点实验室, 浙江 湖州 313001;
3. 中国水产科学研究院东海水产研究所浙江研究中心, 浙江 湖州 313001;
4. 浙江省淡水水产遗传育种重点实验室, 浙江 湖州 313001
梓桐核心区是千岛湖国家级水产种质资源保护区之一,是新安江水系中的一个支流水域。新安江发源于安徽省休宁县六股尖,由歙县街口镇流入浙江省淳安县境内,至建德市梅城镇,旧名徽港、青溪、歙港,主流长365 km,流域面积为11 047 km2,属钱塘江水系[1];而千岛湖(新安江水库)是1960年建成新安江水电站后拦蓄新安江下游而形成的人工湖。目前,千岛湖水域的鱼类研究的报道颇多,主要集中以下几个方面:(1)单个种类的生长[2]、食性[3-4]和遗传[5-6]等方面研究;(2)刺网渔获物[7-8]的相关研究;(3)渔业资源状况及增殖管理[9-12]、保水渔业对渔业生态系统的影响[13-14]等方面研究。但千岛湖国家级水产种质资源保护区——梓桐核心区水域的鱼类资源及群落结构相关调查研究未见报道。由于大坝的建立、水体污染、保水渔业和增殖放流等人为因素造成水体环境的改变,水域生境发生巨大变化,在很大程度上对该水域的鱼类资源产生重要的影响。另外,国内关于水产种质资源保护区鱼类相关的研究[15-18]报道颇多,而针对千岛湖国家级水产种质资源保护区鱼类的相关研究存在空白。因此,本调查于2016年3月至2017年11月期间,针对千岛湖梓桐核心区水域鱼类资源开展了8次采样调查,旨在了解该水域的鱼类资源现状并为水库的渔业资源保护和科学利用提供参考,以填补此处水域鱼类相关研究的空白,也可为水产种质资源保护区的渔业管理提供一定的参考资料。
1 材料与方法 1.1 采样区域和调查时间2016年3月、6月、9月、11月和2017年3月、6月、9月、11月针对千岛湖国家级水产种质资源保护区梓桐核心保护区水域(图 1)内的鱼类资源进行季度调查。
通过下坑坞的渔民在梓桐核心保护区的库区内水域随机下网,每次采样所用网具类型(单层和三层流刺网)、数量基本相同,网目为3、4、5、6、7、8、10和12 cm:3~5 cm网目的网高为1 m,网长为50~100 m,下网3~5片;6~8 cm网目的网高为3 m,网长为100 m,下网1~2片;10和12 cm网目的网高为10~12 m,网长为100~200 m,下网1~2片。捕捞时间和渔获物种类鉴定等参照文献[19]。
1.3 数据处理与分析采用相对重要性指数对保护区鱼类优势种组成进行分析;利用马格利夫多样性指数、香农-维纳多样性指数、皮诺均匀度指数和辛普森多样性指数对调查水域鱼类群落多样性进行分析。
相对重要性指数(relative important index, IRI)[20]:
式中:IRI为某种类IRI值;N为某种类的尾数百分比;W为某种类质量百分比;F为某种类季度出现频率。
优势种划分参考各种类的IRI值,采用IRI >2 000为优势种,2 000>IRI>1 000为重要种,1 000>IRI>500为常见种,500>IRI>300为一般种,IRI < 300为少有种的划分标准[21]。
香农-维纳多样性指数(Shannon-Wener index, H′)[22-23]:
皮诺均匀度指数(Pielou index, J′)[24]:
马格利夫多样性指数(Margalef index, D)[25]:
辛普森多样性指数(Simpson index, λ)[26]:
式中:S为渔获物的种类总数;N为渔获物的总个体数;Pi为第i种的个体体质量占总体质量的百分比;Hmax表示H′的最大值。
据种类质量百分比(W)建立矩阵,以Bray-Curtis相似性系数(Bray-Curtis similarity index)为基础,构建相似性矩阵,利用非加权组平均法(unweighted pair group method arithmetic averaging, UPGMA)进行聚类与非线性多维标序(non-metric multidimensional scaling, NMDS)进行排序,分别构建聚类分析图和NMDS 2D图,两种分析方法的结果可相互验证其正确性[27],以分析调查水域的鱼类群落特征。NMDS分析结果由胁强系数(stress)来衡量优劣,衡量标准参考文献[28]。并运用相似性分析(analysis of similarities, ANOSIM)进行差异显著性检验。
数据统计和作图运用Excel 2007和Primer 5软件。
2 结果与分析 2.1 群落种类组成及生态类型调查共采集到鱼类样本1 528尾,总质量达341.70 kg。种类分析,全部鉴定到种,共鉴定鱼类43种,隶属4目8科30属,如表 1所示,其中:鲤科鱼类28种,占总种数的65.12%,
在8次的水域调查鱼类样本采集中,黄尾鲴(Xenocypris davidi)、鲢(Hypophthalmichthys molitrix)、翘嘴鲌(Culter alburnus)、大眼华鳊(Sinibrama macrops)和华鳈(Sarcocheilichthys sinensis)的出现频率达到了100%(表 2)。按划分方法,梓桐支流水域的鱼类群落优势种有3种,分别为黄尾鲴、银飘鱼(Pseudolaubuca sinensis)和鳙(Aristichthys nobilis);重要的鱼类种类也有3种,分别为鲢、翘嘴鲌和大眼华鳊;常见鱼类也为3种,分别为银
2016年期间:Margalef丰富度指数(D)全年变化范围为3.17~4.49,整体上呈现逐渐减少的趋势;Simpson多样性指数(λ)全年变化范围为0.68~0.89,除6月份值最小外,其他3个月的指数变化无明显差异,变化趋势与Margalef丰富度指数(D)相一致;Pielou均匀度指数(J′)全年变化范围为0.61~0.84,变化趋势呈先减后增的趋势;Shannon-Wiener多样性指数(H′)全年变化范围为1.79~2.55,全年情况与Simpson多样性指数(λ)相一致(表 3)。2017年期间,Margalef丰富度指数(D)全年变化范围为3.66~5.01,呈现先增后减的变化趋势;Simpson多样性指数(λ)全年变化范围为0.84~0.88,4个季度的指数变化差异不明显;Pielou均匀度指数(J′)全年变化范围为0.72~0.79,全年变化与2016年Simpson多样性指数(λ)变化相似;Shannon-Wiener多样性指数(H′)全年变化范围为2.26~2.48,整体上呈现先增后减的趋势。
据Bray-Curtis相似性指数分析结果,在83.19%的相似性水平上,可将8个采样时间划分为2个群组。群组Ⅰ包括2016年3月、9月、11月和2017年9月、11月,群组Ⅱ包括2016年6月和2017年3月、6月(图 2)。NMDS排序的胁强系数(stress)为0,说明NMDS排序结果很好,可以用二维点阵图表示。ANOSIM分析结果R=0.99,P<0.05,表明2个群组间存在显著差异,即鱼类群落结构季节变化显著。
群落结构中的重要种或优势种占据重要的地位,其对群落结构起着重中之重的贡献作用。因此,对该支流的6种重要种或优势种进行Bray-Curtis相似性指数分析,结果显示在69.15%的相似性指数水平上,可将支流水域内的6种重要或优势种鱼类划分为2个类群,类群Ⅰ包括大眼华鳊、银飘鱼、翘嘴鲌和黄尾鲴4种,而类群Ⅱ包含鲢和鳙2种(图 3)。NMDS排序的胁强系数(stress)为0,表明NMDS排序图与相似性分析结果之间可以很好地相互解释。ANOSIM分析结果R=0.714,P>0.05,表明2个类群间无显著差异。
丰度-生物量比较曲线(ABC曲线)可以用来反映群落中种类组成变化情况[29]。调查水域鱼类群落的生物量优势度曲线始终位于丰度优势度曲线的上方的月份有2016年3月和2017年3月(图 4),且W值为正,据李圣法[30]的研究分析,可判定该时期的鱼类群落结构处于稳定状态,未受到较大干扰。2016年和2017年两年中的6月、9月和11月的ABC曲线都出现了交叉变换,说明该时期的群落结构处于中度干扰状态。
本次调查,千岛湖梓桐核心保护区的鱼类群落呈现出以湖泊定居性鱼类为主和鲤科鱼类在种数上占绝对优势的特征,与以往对千岛湖鱼类资源的研究[7-8]结论相一致。本次梓桐核心保护区水域调查采集了43种鱼类,与何喜光等[7]和侯文华等[8]利用流刺网调查的结果相比,种类组成上存在明显的差异性。可见,流刺网对于保护区鱼类种类组成的调查存在局限性,种类难以采集齐全。另外,新安江大坝的形成阻断了长距离溯河洄游性鱼类的产卵通道,此类鱼在调查水域捕获的数量较少,如短颌鲚仅捕获1尾,造成这样的结果可能受到以下两个方面的因素影响:一方面,采样方法单一,捕获的种类和数量受到一定的限制;另一方面,这些鱼类可能演变成定居性鱼类后较难繁殖,从而造成种群资源衰退。大坝的建立改变了水域环境,影响了鱼类群落种类组成[31]。
一方面,外来种对鱼类群落结构存在一定影响。调查发现了外来物种美洲鲥(16尾)、斑点叉尾
水库蓄水前,以马口鱼(Opsariichthys bidens)、宽鳍
鱼类群落多样性可侧面说明群落结构的稳定性,群落物种越丰富,各种类个体数分布越均匀,则丰富度和均匀性指数越高,群落结果越稳定[34]。本次调查结果显示:2016年期间,春季(3月)鱼类群落Shannon-Wiener多样性指数(H′)、Simpson多样性指数(λ)和Pielou均匀度指数(J′)均最大,Margalef丰富度指数(D)相对较大,本研究认为调查水域2016年春季的种类丰富度最大,且种群分布均匀。相反,夏季(6月)4个生物多样性指数均最小,说明夏季鱼类种类丰富度不高、各种类个体分布不均,这可能与水位高低变化[35]、天气等以及流刺网捕获形式等的影响有关。2017年期间,夏季(6月)鱼类群落Shannon-Wiener多样性指数(H′)和Margalef丰富度指数(D)最大,Simpson多样性指数(λ)和Pielou均匀度指数(J′)相对较大,本研究认为调查水域2017年夏季的种类丰富度最大,且种群分布均匀。上述结果与调查方法存在很大的关系[36],可以推测,单一性的设置流刺网会导致结果形成偏差。另外,水位变化、偷钓和旅游等因素也对鱼类群落多样性存在影响。因此,不同季节鱼类群落多样性的差异与水环境改变和调查方法存在相应的关系。
聚类和排序分析结果显示,在83.19%的相似性指数水平上,将8个月份划分为2个类群,ANOSIM分析显示(P<0.05),调查水域鱼类群落结构季节变化差异显著,可能与优势种或者重要种的各月所占生物量百分比变化有关,如鲴亚科鱼类在2016年6月、2017年3月和6月的捕获量占当次总捕获量的百分比较大。类群Ⅰ中的黄尾鲴、翘嘴鲌和大眼华鳊等是当地的主要捕捞对象,而类群Ⅱ中的鲢和鳙是保护对象,这是引起两大类鱼类群落差异显著的主要原因。在69.15%的相似性指数水平上,将优势种和重要种可划分为2个群组,ANOSIM分析显示(P>0.05),但两者之间不存在差异性,这可能与当地主要捕捞对象在数量和生物量上占据优势有关。
丰度-生物量比较曲线(ABC曲线)结果显示:2016年3月(春季)和2017年3月(春季)调查水域鱼类群落结构未受到干扰;其他调查季节鱼类群落结构处于中度干扰状态,这可能与水体环境[34]、水体沿岸周边环境和人为干扰因素等存在关系。
3.4 启示与建议在调查过程中,发现存在外来物种,外来物种一旦形成自然种群,具有形成生物入侵种的潜能,如蓝鳃太阳鱼具有捕食其他鱼类卵粒的习性,这将对鱼类群落产生影响。因此,有必要加强对外来物种的监测,如加强对水域及周边的养殖区引进外来种养殖工作的监管,实时监测水域内的外来种状况,开展预防外来种入侵的相关研究工作,以建立相应的预防应对措施。
[1] |
林衍峰. 新安江上游渔业资源变迁与保护[J]. 中国水产, 2006(1): 81-82. LIN Y F. Changes and protection of fishery resources in the upper reaches of Xin'an River[J]. China Fisheries, 2006(1): 81-82. DOI:10.3969/j.issn.1002-6681.2006.01.046 |
[2] |
刘国栋, 何光喜, 刘其根, 等. 千岛湖大眼华鳊年龄、生长和繁殖的初步研究[J]. 上海海洋大学学报, 2011, 20(3): 382-391. LIU G D, HE G X, LIU Q G, et al. Primary study on the age, growth and reproduction of Sinibrama macrops in Qiandao Lake[J]. Journal of Shanghai University, 2011, 20(3): 382-391. |
[3] |
宋固.千岛湖鲢、鳙鱼食性研究[D].上海: 上海海洋大学, 2013. SONG G. Feeding habits of silver carp and bighead crap in Qiandaohu Lake[D]. Shanghai: Shanghai Ocean University, 2013. http://cdmd.cnki.com.cn/Article/CDMD-10264-1014013545.htm |
[4] |
宋固, 胡梦红, 刘其根. 运用稳定同位素技术研究千岛湖秋季刺网渔获物的食性和营养级[J]. 上海海洋大学学报, 2014, 23(1): 117-122. SONG G, HU M H, LIU Q G. Study on the feeding habits and trophic levels of the gillnet catches of Qiandao Lake in autumn by stable isotope technology[J]. Journal of Shanghai Ocean University, 2014, 23(1): 117-122. |
[5] |
张峻德, 赵良杰, 刘其根. 千岛湖细鳞鲴和黄尾鲴COI种群遗传结构比较的初步研究[J]. 淡水渔业, 2013, 43(5): 3-7. ZHANG J D, ZHAO L J, LIU Q G. Primary study on comparison of population genetic structure between Xenocypris microlepis and X.davidi in Qiandao lake base on mitochondrial COIgene[J]. Freshwater Fisheries, 2013, 43(5): 3-7. DOI:10.3969/j.issn.1000-6907.2013.05.001 |
[6] |
曲宪成, 张勇, 刘其根, 等. 千岛湖细鳞鲴种群线粒体COI基因变异及遗传分化研究[J]. 湖南农业科学, 2011(15): 150-153. QU X C, ZHANG Y, LIU Q G, et al. Gene mutation and genetic differentiation of mitochondrial COI in Xenocypris microlepis in Qiandaohu Lake[J]. Hunan Agricultural Sciences, 2011(15): 150-153. DOI:10.3969/j.issn.1006-060X.2011.15.045 |
[7] |
何光喜, 张峻德, 刘其根, 等. 基于单层刺网渔获物的千岛湖鱼类群落结构初步研究[J]. 水产科技情报, 2011, 38(4): 197-203. HE G X, ZHANG J D, LIU Q G, et al. Preliminary study on fish community structure of Qiandao Lake based on single layer gillnet catch[J]. Fisheries Science & Technology Information, 2011, 38(4): 197-203. DOI:10.3969/j.issn.1001-1994.2011.04.010 |
[8] |
侯文华, 胡梦红, 刘其根. 千岛湖刺网渔获物中鱼类群落结构及多样性特征分析[J]. 上海海洋大学学报, 2014, 23(6): 826-833. HOU W H, HU M H, LIU Q G. Community structure and diversity of fishes in Qiandao Lake[J]. Journal of Shanghai Ocean University, 2014, 23(6): 826-833. |
[9] |
罗仙池, 吴振兴. 千岛湖的渔业资源现状及资源增殖途径探讨[J]. 水产科技情报, 1992, 19(2): 33-37. LUO X C, WU Z X. Present situation of fishery resources in the Thousand Island Lake and discussion on the way of their reproduction[J]. Fisheries Science & Technology Information, 1992, 19(2): 33-37. |
[10] |
徐田祥, 蔡茂胜, 项为民. 千岛湖渔业资源下降原因及对策[J]. 水利渔业, 2000, 20(6): 35-36. XU T X, CAI M S, XIANG W M. Reasons for the decline of fisheries resources in Qiandaohu Lake and strategies[J]. Reservoir Fisheries, 2000, 20(6): 35-36. DOI:10.3969/j.issn.1003-1278.2000.06.021 |
[11] |
方中平, 向国春. 千岛湖渔业生态资源增殖保护的对策与措施[J]. 水利渔业, 2007, 27(4): 50-52. FANG Z P, XIANG G C. Countermeasures and measures for the protection of fishery ecological resources in Qiandao Lake[J]. Reservoir Fisheries, 2007, 27(4): 50-52. DOI:10.3969/j.issn.1003-1278.2007.04.019 |
[12] |
李亚楠, 臧斯颖, 侯文华, 等. 主要外来鱼类对千岛湖生态系统的潜在影响及对策[J]. 生物学通报, 2014, 49(9): 3-5. LI Y N, ZANG S Y, HOU W H, et al. Potential impacts and countermeasures of major exotic fishes on the Qiandao Lake ecosystem[J]. Bulletin of Biology, 2014, 49(9): 3-5. |
[13] |
LIU Q G. Aquatic environmental protection oriented fishery in Lake Qiandaohu and its influences on lake ecosystem[J]. Shanghai:East China Normal University, 2005. |
[14] |
刘其根, 王钰博, 陈立侨, 等. 保水渔业对千岛湖食物网结构及其相互作用的影响[J]. 生态学报, 2010, 30(10): 2774-2783. LIU Q G, WANG Y B, CHEN L Q, et al. Impacts of aquatic environment protection oriented fishery on the structure of food web in Lake Qiandaohu[J]. Acta Ecologica Sinica, 2010, 30(10): 2774-2783. |
[15] |
王开锋, 靳铁治, 孙长铭, 等. 陕西库峪岷县高原鳅水产种质资源保护区的鱼类及其多样性[J]. 淡水渔业, 2011, 41(4): 16-20. WANG K F, JIN T Z, SUN C M, et al. Fish and their diversities in Shaanxi Kuyu aquatic resources reserve[J]. Freshwater Fisheries, 2011, 41(4): 16-20. DOI:10.3969/j.issn.1000-6907.2011.04.003 |
[16] |
何美峰, 袁定清, 崔利峰, 等. 汀江大刺鳅国家级水产种质资源保护区鱼类群落结构及其影响因子[J]. 上海海洋大学学报, 2015, 24(1): 121-129. HE M F, YUAN D Q, CUI L F, et al. Study on fish community and impact factors in Mastacembelusaculedtus national aquatic reserve in Tingjiang River[J]. Journal of Shanghai Fisheries University, 2015, 24(1): 121-129. |
[17] |
李万宝. 福建建溪细鳞鲴国家级水产种质资源保护区的鱼类多样性[J]. 淡水渔业, 2016, 46(4): 31-37. LI W B. Species diversity of fish community of Xenocypris microlepis national aquatic germplasm resource reservein Jianxi River[J]. Freshwater Fisheries, 2016, 46(4): 31-37. DOI:10.3969/j.issn.1000-6907.2016.04.005 |
[18] |
刘凯, 段金荣, 许萌霆, 等. 淮河荆涂峡鲤、长吻鮠国家级水产种质资源保护区冬季渔业群落多样性研究[J]. 安徽农业大学学报, 2015, 42(3): 424-428. LIU K, DUAN J R, XU M T, et al. Winter biodiversity of the fishery community in Huaihe Cyprinus carpio and Leiocassis longirostris national aquatic germplasm resources conservation area[J]. Journal of Anhui Agricultural University, 2015, 42(3): 424-428. |
[19] |
郝雅宾, 刘金殿, 郭爱环, 等. 钱塘江富阳段鱼类资源现状[J]. 淡水渔业, 2017, 47(3): 45-51. HAO Y B, LIU J D, GUO A H, et al. Current status of fishery resources in Fuyang section of Qiantang River[J]. Freshwater Fisheries, 2017, 47(3): 45-51. DOI:10.3969/j.issn.1000-6907.2017.03.007 |
[20] |
PINKAS L, OLIPHANT M S, IVEROSN K. Food habits of albacore, bluefin tuna, and bonito in California waters[J]. Fishery Bulletin, 1971, 152: 9-10. |
[21] |
张建禄, 边坤, 靳铁治, 等. 秦岭黑河流域鱼类资源现状调查[J]. 淡水渔业, 2016, 46(1): 103-108. ZHANG J L, BIAN K, JIN T Z, et al. Investigation on status of fish resources in Qinling Heihe River basin[J]. Freshwater Fisheries, 2016, 46(1): 103-108. DOI:10.3969/j.issn.1000-6907.2016.01.017 |
[22] |
WILHMJ L. Use of biomass units in Shannon's formula[J]. Ecology, 1968, 49(1): 153-156. DOI:10.2307/1933573 |
[23] |
SHANNON C E, WEAVER W W. The mathematical theory of communication[M]. Urbana: University of Illinois Press, 1963: 1-65.
|
[24] |
PIELOU E C. The measurement of diversity in different types of biological collections[J]. Journal of Theoretical Biology, 1966, 13: 131-144. DOI:10.1016/0022-5193(66)90013-0 |
[25] |
MARGALEF R. Information theory in ecology[J]. General Systems, 1958, 3: 36-71. |
[26] |
LUDWIG J A, REYNOLDS J F. Statistical Ecology:a primer on methods and computing[M]. New York: John Wiley & Sons, 1988: 197.
|
[27] |
BRAZNER J C, BEALS E W. Patterns in fish assemblages from coastal wetland and beach habitats in Green Bay, Lake Michigan:a multivariate analysis of abiotic and biotic forcing factors[J]. Canadian Journal of Fisheries and Aquatic Sciences, 1997, 54(8): 1743-1761. DOI:10.1139/f97-079 |
[28] |
CLARKE K R, WARWICK R M. Change in marine communities:an approach to statistical analysis and interpretation[M]. 2nd ed. Plymouth: Primer-E, 2001.
|
[29] |
BLANCHARD F, LELOC'H F, HILY C, et al. Fishing effects on diversity, size and community structure of the benthic invertebrate and fish megafauna on the Bay of Biscay coast of France[J]. MarineEcology Progress Series, 2004, 280: 249-260. DOI:10.3354/meps280249 |
[30] |
李圣法. 以数量生物量比较曲线评价东海鱼类群落的状况[J]. 中国水产科学, 2008, 15(1): 136-144. LI S F. Status of fish community in East China Sea using the method of abundance-biomass comparison (ABC)curve[J]. Journal of Fishery Sciences of China, 2008, 15(1): 136-144. DOI:10.3321/j.issn:1005-8737.2008.01.018 |
[31] |
王崇, 谢山, 王进国, 等. 红水河龙滩水库鱼类资源调查[J]. 水生态学杂志, 2014, 35(2): 39-48. WANG C, XIE S, WANG J G, et al. Investigationof fish resources in Longtan reservoir of Hongshuihe River[J]. Journal of Hydroecology, 2014, 35(2): 39-48. DOI:10.3969/j.issn.1674-3075.2014.02.007 |
[32] |
王敏, 朱峰跃, 刘绍平, 等. 三峡库区汉丰湖鱼类群落结构的季节变化[J]. 湖泊科学, 2017, 29(2): 439-447. WANG M, ZHU F Y, LIU S P, et al. Seasonal variations of fish community structure of Lake Hanfeng in Three Gorges reservoir region[J]. Journal of Lake Science, 2017, 29(2): 439-447. |
[33] |
冯三杰, 王从锋, 谭宏, 等. 北盘江董箐水电站库区鱼类资源调查[J]. 水生态学杂志, 2018, 39(2): 70-76. FENG S J, WANG C F, TAN H, et al. Fish resource investigation of Dongqing hydropower station reservoir in Beipan River basin[J]. Journal of Hydroecology, 2018, 39(2): 70-76. |
[34] |
MATTHEWS W J. Patterns in freshwater fish ecology[M]. New York: Chapman and Hall Publishers, 1998.
|
[35] |
丁庆秋, 彭建华, 杨志, 等. 三峡水库高、低水位下汉丰湖鱼类资源变化特征[J]. 水生态学杂志, 2015, 36(3): 1-9. DING Q Q, PENG J H, YANG Z, et al. Variation in the Hanfeng Lake fishery resource between low and high water levels in Three Gorges reservoir[J]. Journal of Hydroecology, 2015, 36(3): 1-9. |
[36] |
WANG Y G, YE Y M, MILTON D A. Efficient designs for sampling and subsampling in fisheries research based on ranked sets[J]. ICES Journal of Marine Science, 2009, 66(5): 928-934. DOI:10.1093/icesjms/fsp112 |
2. Key Laboratory of Healthy Freshwater Aquaculture, Ministry of Agriculture and Rural Affairs, Huzhou 313001, Zhejiang, China;
3. Zhejiang Research Center of East China Sea Fishery Research Institute, Huzhou 313001, Zhejiang, China;
4. Key Laboratory of Freshwater Aquaculture Genetic and Breeding of Zhejiang Province, Huzhou 313001, Zhejiang, China