四川动物  2018, Vol. 37 Issue (1): 67-73

扩展功能

文章信息

周传江, 顾钱洪, 孟晓林, 李斌顺, 赵文武, 庞显炳, 张建新, 宋东蓥, 聂国兴
ZHOU Chuanjiang, GU Qianhong, MENG Xiaolin, LI Binshun, ZHAO Wenwu, PANG Xianbing, ZHANG Jianxin, SONG Dongying, NIE Guoxing
基于多来源数据分析河南省刺鱼目新纪录种——中华多刺鱼
Analysis of New Fish Record of Pungitius sinensis in Henan Province Based on Multiple Data Resource
四川动物, 2018, 37(1): 67-73
Sichuan Journal of Zoology, 2018, 37(1): 67-73
10.11984/j.issn.1000-7083.20170037

文章历史

收稿日期: 2017-02-09
接受日期: 2017-11-29
基于多来源数据分析河南省刺鱼目新纪录种——中华多刺鱼
周传江1 , 顾钱洪1 , 孟晓林1 , 李斌顺2 , 赵文武3 , 庞显炳3 , 张建新1 , 宋东蓥1 , 聂国兴1*     
1. 河南省水产动物养殖工程技术研究中心, 河南师范大学水产学院, 河南新乡 453007
2. 林州市水产站, 河南林州 456550
3. 安阳市水产科学研究所, 河南安阳 455000
摘要:2015年10月,河南省鱼类资源调查队在河南省林州市进行鱼类资源普查时,采集到30尾刺鱼目Gasterosteiformes鱼类,综合采用形态学和分子系统学的方法确定其为中华多刺鱼Pungitius sinensis,为河南省新纪录种,标本保存于河南师范大学水产学院鱼类标本室。本文对该鱼的主要鉴定特征、分布区域、可能的来源、生存环境等做了初步分析。该鱼在我国北方(辽宁、黑龙江、吉林、天津、河北等省市)及韩国、日本北方均有分布,记录最南分布区域为北京市,为北京市二级保护鱼类。此次报道进一步扩大了该鱼最南端分布界限,为进一步开展该鱼的基础研究提供了重要参考资料。
关键词刺鱼目     中华多刺鱼     林州     河南省     新纪录    
Analysis of New Fish Record of Pungitius sinensis in Henan Province Based on Multiple Data Resource
ZHOU Chuanjiang 1, GU Qianhong 1, MENG Xiaolin 1, LI Binshun 2, ZHAO Wenwu 3, PANG Xianbing 3, ZHANG Jianxin 1, SONG Dongying 1, NIE Guoxing 1*     
1. College of Fisheries, Engineering Technology Research Center of Henan Province for Aquatic Animal Cultivation, Henan Normal University, Xinxiang, Henan Province 453007, China;
2. Linzhou City Aquatic Station, Linzhou, Henan Province 456550, China;
3. Anyang City Fisheries Science Institute, Anyang, Henan Province 455000, China
Abstract: In this study, a total of 30 Gasterosteiformes fishes were collected when the fish investigation team of Henan province worked in Linzhou city in October, 2015. By using integrated application of morphological and molecular systematic methods, these fish were identified as Pungitius sinensis which is a new record in Henan province. The specimen was preserved in the college of Fisheries of Henan Normal University. Subsequently, preliminarily analyses including mainly identification characteristics, distribution region, possible sources and living environment were performed. This fish is widely distributed in northern China (e.g. Liaoning province, Heilongjiang province, Jilin province, Tianjin city, and Hebei province) and The Republic of Korea, northern Japan. Beijing city is the most southern distribution region of this fish species in previous record, and this fish was characterized as Class Ⅱ municipally protected fish in Beijing. Therefore, the southernmost distribution region of this fish was further expanded by this report, and provided important reference material for further fundamental research of this fish.
Keywords: Gasterosteiformes     Pungitius sinensis     Linzhou     Henan province     new record    

中华多刺鱼Pungitius sinensis隶属于刺鱼目Gasterosteiformes刺鱼科Gasterosteidae多刺鱼属Pungitius。关于该属及近缘种的研究涉及到种群、生态学、系统发育、进化生物学等方面(Chae,1988Chae & Yang, 1988, 1993Takahashi & Goto,2001Watanabe et al., 2003Koizumi et al., 2007卞绍雷等,2008Hubert et al., 2008Ki et al., 2008Aldenhoven et al., 2010Shikano et al., 2010Teacher et al., 2011Hwang et al., 2012McCusker et al., 2013Geiger et al., 2014Denys et al., 2015Rezansoff et al., 2015Thalinger et al., 2016)。据资料记载,我国有4种刺鱼分布,即刺鱼属Gastersteus的三刺鱼Gastersteus aculeatus(郑葆珊等,1980成庆泰,周才武,1997)和多刺鱼属的中华多刺鱼P. sinensis、刺鱼P. pungitius(即八棘多刺鱼;任慕莲等,1981)、多刺鱼P. pungitius sinensis(尼科尔斯基,1960)。中华多刺鱼在我国北方分布广泛,如黑龙江(尼科尔斯基,1960任慕莲,1981张觉民,1995)、图们江(郑葆珊等,1980)、北京(王鸿媛,1984)、辽宁(刘婵馨,秦克静,1987)、天津(天津市水产学会,1990李明德,2011)、河北(李明德,杨竹舫,1992王所安等,2000)、山东(成庆泰,周才武,1997)。《中国鱼类系统检索》(成庆泰,郑葆珊,1987)中记录中华多刺鱼在黑龙江、松花江、图们江、辽河、内蒙古、河北及长江流域有分布,但《长江鱼类》(湖北省水生生物研究所鱼类研究室,1976)中并没有记录该种,所以卞绍雷等(2008)认为北京及其周边地区可能为中华多刺鱼在亚洲大陆的最南分布点。2015年10月,河南省鱼类资源调查队在河南省林州市进行鱼类资源普查时,采集到30尾刺鱼目鱼类,综合采用形态学和分子系统学的方法确定其为中华多刺鱼,本种系河南省新纪录种,标本保存于河南师范大学水产学院鱼类标本室。本文初步分析了该鱼的主要鉴定特征、在河南省的分布区域、可能来源、生存环境等。

1 材料与方法 1.1 材料

30尾实验鱼均采自河南省林州市五龙镇河头村河头大桥(113°56'27.89"E,35°50'14.59"N),地笼诱捕获得,13尾100%乙醇和17尾10%甲醛溶液整体固定,带回实验室进行测量和后续实验。

1.2 形态测量

对17尾甲醛溶液固定标本用数显游标卡尺(单位:mm)测量全长、体长、头长、眼径、吻长、眼后头长、吻宽、体高、尾柄长、尾柄高、吻至背鳍长、第一背鳍长、第二背鳍长、背鳍高、臀鳍基长、臀鳍高、胸鳍基长、腹鳍基长,并计数背鳍棘数目、臀鳍和胸鳍鳍条数目及棘数。

1.3 分子生物学及分子系统学方法鉴定 1.3.1 基因组DNA提取

取适量经酒精固定的标本背部肌肉组织置于1.5 mL离心管中,55 ℃烘干4 h后采用酚-氯仿提取法提取基因组DNA(萨姆布鲁克,拉塞尔,2008)。

1.3.2 PCR扩增获得目的基因

PCR扩增线粒体COⅠ基因的引物对为HCO:5'-AAACYTCWGGGTGVCCAAAGAAYC-3',LCO:5'-TACTAAYCACAAAGATATTGGCAC-3'(Folmer et al., 1994)。30 μL PCR反应体系包含:15 μL mix,上、下游引物各1 μL,基因组DNA 1 μL,灭菌水12 μL。PCR反应条件为:94 ℃预变性5 min;94 ℃变性30 s,54 ℃退火30 s,72 ℃延伸30 s,30个循环;72 ℃延伸7 min。PCR产物分别用上、下游引物测序。

1.4 数据分析

从测序公司返回的测序结果为峰图文件,采用DNASTAR的Seqman组装为一个contig,通过NCBI网站进行BLAST分析(https://blast.ncbi.nlm.nih.gov/Blast.cgi),确认为目的基因COⅠ,同时从NCBI下载刺鱼目及其近缘种的COⅠ基因序列共71条(表 1),用MEGA 7.0(Kumar et al., 2016)构建系统发育树和计算遗传距离。系统发育树构建采用NJ原理Kimura双参法,1 000个重复抽样,其他参数为默认设置。遗传距离估算依据构建的系统发育树将同种的分为同一个群,之后分别计算群内和群间的遗传距离,进而对这些序列及物种之间的关系进行评价。

表 1 样品信息及来源 Table 1 Information of specimen and resource
物种名称 NCBI登录号或样品编号 数据来源
内群 中华多刺鱼 Pungitius sinensis HTU402427~HTU402431、HTU402436~HTU402440 本研究
沼泽多刺鱼 P. hellenicus KJ554130、KJ554255、KJ554514 Geiger et al., 2014
平腹多刺鱼 P. platygaster KJ554195、KJ554199、KJ554378 Geiger et al., 2014
光身多刺鱼 P. laevis KJ554327、KJ554413 Geiger et al., 2014
图们江多刺鱼 P. tymensis KM019403 Geiger et al., 2014
八棘多刺鱼 P. pungitius EU524319~EU524321、EU525105~EU525112 Hubert et al., 2008
八棘多刺鱼 P. pungitius KJ128593、KJ128594 Bergsten et al., unpublished
八棘多刺鱼 P. pungitius KT630717、KT630718 Rodgers & Mock,unpublished
八棘多刺鱼 P. pungitius KR477087~KR477089、KR477243~KR477245 Thalinger et al., 2016
八棘多刺鱼 P. pungitius HQ712742 Mecklenburg et al., 2011
八棘多刺鱼 P. pungitius KC015860、KC015861 McCusker et al., 2013
八棘多刺鱼 P. pungitius KM287043~KM287048 Knebelsberger et al., 2015
八棘多刺鱼 P. pungitius KR862862 Denys et al., 2015
外群 Gasterosteus sp. KJ553578、KJ553600、KJ553622、KJ553623 Geiger et al., 2014
二棘刺鱼 G. wheatlandi EU524640~EU524642 Hubert et al., 2008
三刺鱼 G. aculeatus KP823178、KP823179 Rezansoff et al., 2015
三刺鱼 G. aculeatus KR862765~KR862768 Denys et al., 2015
裸尾刺鱼G. gymnurus KR862829~KR862832 Denys et al., 2015
冰岛刺鱼 G. islandicus KR862833~KR862836 Denys et al., 2015
2 结果 2.1 形态学分析结果

形态鉴定30尾,均符合中华多刺鱼的主要鉴别特征:眼较大,体长,呈侧扁形,尾柄极细,有明显侧棱。侧线上有31~33个骨板,尾截形或微凹入(图版Ⅰ:1),具2个背鳍,第一背鳍为游离的棘,具有交错排列的8~9个棘(图版Ⅰ:2),第二背鳍和臀鳍相对,8~10个鳍条,臀鳍具1棘8~10个鳍条,胸鳍为8~9个鳍条,腹鳍为1棘1~2个鳍条(图版Ⅰ:3)。形态测量显示(n=17):全长49.22±3.90,体长41.79±5.57,头长11.85±1.35,眼径3.57±0.47,吻长3.02±0.70,眼后头长5.83±0.88,尾柄长7.18±0.80,尾柄高1.34±0.16,吻至背鳍长13.28±0.99,第一背鳍长13.67±1.64,第二背鳍长10.50±1.31,背鳍高4.33±0.87,臀鳍基长10.47±1.16,臀鳍高4.75±0.57,胸鳍基长2.48±0.30,腹鳍基长6.93±0.72。体长为体高的4.4~6.9倍,为头长的3.0~4.5倍。头长为吻长的2.4~6.2倍,为眼径的2.5~4.0倍,为尾柄长的1.2~2.0倍,为尾柄高的7.0~12.0倍,尾柄长为尾柄高的4.3~7.3倍。

图版Ⅰ 中华多刺鱼的主要鉴别特征及生境 PlateⅠ The main diagnostic characteristics and habitat of Pungitius sinensis 1.中华多刺鱼整体;2.第一背鳍的棘(白色箭头所指)和第二背鳍(黄色箭头所指);3.腹鳍及臀鳍的棘(白色箭头所指)和骨板(黄色箭头所指);4.中华多刺鱼生活水体环境;比例尺=1 cm 1. Whole body of P. sinensis; 2. Spines of first dorsal fin (white arrow) and the second dorsal fin (yellow arrow); 3. Spines of ventral fin and anal fin (white arrow) and bone plate (yellow arrow); 4. Living water environment of P. sinensis; scale bar=1 cm
2.2 分子生物学及分子系统学结果

PCR扩增共获得中华多刺鱼线粒体COⅠ基因序列10条,与NCBI下载的刺鱼目及近缘种的61条COⅠ基因序列一起比对后发现:589 bp序列中保守位点448个,简约信息位点133个,T、C、A、G含量分别为28.4%、30.2%、23.6%、17.8%,A+T含量高于G+C含量。71条序列构建的NJ树显示(图 1),外群聚为单独的一支,多刺鱼属构成支持率为100%的单系类群,但各支情况有所不同,图们江多刺鱼P. tymensis位于多刺鱼属的基部。光身多刺鱼P. laevis、本研究的中华多刺鱼、沼泽多刺鱼P. hellenicus、平腹多刺鱼P. platygaster的所有个体均聚在一起构成高支持率的单系。而广布种八棘多刺鱼则分为多个分支,可能因其分布范围较广,不同区域的种群已经发生了分化。

图 1 基于线粒体COⅠ基因的多刺鱼属系统发育关系 Fig. 1 Phylogenetic relationship of the genus Pungitius inferred from mitochondrial COⅠ gene

同种鱼的群内遗传距离为0~0.006(如八棘多刺鱼为0.006、仅有一个个体的图们江多刺鱼无法计算、光身多刺鱼为0.002、沼泽多刺鱼为0、平腹多刺鱼为0.001、中华多刺鱼为0),可以判断依据系统发育结果进行群的划分是可靠的,同种个体均聚在一支且遗传距离小于0.02。八棘多刺鱼和光身多刺鱼的遗传距离比较近(0.008<0.02),且在系统发育结果中,光身多刺鱼也是嵌套在八棘多刺鱼中,初步推测八棘多刺鱼和光身多刺鱼可能为同一物种。其他不同物种之间遗传距离均大于0.02,初步推断为不同物种(表 2)。

表 2 不同群(物种)之间遗传距离计算结果 Table 2 The genetic distance between different groups (species) of genus Pungitius
P. pungitius P. tymensis P. laevis P. hellenicus P. platygaster
八棘多刺鱼 Pungitius pungitius
图们江多刺鱼 P. tymensis 0.069
光身多刺鱼 P. laevis 0.008 0.074
沼泽多刺鱼 P. hellenicus 0.059 0.091 0.060
平腹多刺鱼 P. platygaster 0.061 0.089 0.062 0.043
中华多刺鱼 P. sinensis 0.022 0.063 0.026 0.063 0.068
3 讨论 3.1 中华多刺鱼物种的确定

在中国分布的刺鱼类有以下4种,一种为分布于图们江(郑葆珊等,1980)和山东(成庆泰,周才武,1997)的三刺鱼,另一种为刺鱼(任慕莲等, 1981, 即八棘多刺鱼),还有一种是多刺鱼(尼科尔斯基,1960)。尼科尔斯基在《黑龙江流域鱼类》提到除了三刺鱼外,同时也提到了其他多个地方志均对中华多刺鱼有记载(王鸿媛,1984刘婵馨,秦克静,1987天津水产学会,1990李明德,杨竹舫,1992成庆泰,周才武,1997王所安等,2000李明德,2011)。这些文献中记载的中华多刺鱼重要形态鉴定特征基本一致,但对比记录数据,发现形态长度及长度指标比例存在较大偏差,同时本文获得的数据也和已有报道有一定出入,可能是采集地不同及测量误差造成的,如本文中头长为尾柄长的1.2~2.0倍,为尾柄高的7.0~12.0倍,而《图们江鱼类志》(郑葆珊等,1980)记载的中华多刺鱼相应指标为4.6~4.8倍和1.7~1.9倍。结合系统发育研究,本研究样本单独聚为一支,群内遗传距离为0,与其他已知的多刺鱼属鱼类的遗传距离均超过了同种的遗传距离,可以判定为中华多刺鱼。

3.2 中华多刺鱼的分布范围及生境情况

卞绍雷等(2008)认为北京是中华多刺鱼分布的最南边界,此次在河南林州的发现进一步扩大了该鱼的分布范围。《中国鱼类系统检索》(成庆泰,郑葆珊,1987)提到中华多刺鱼在长江流域也有分布,但《长江鱼类》(湖北省水生生物研究所鱼类研究室,1976)并没有该鱼的记载,因此该鱼在长江流域的分布情况需进一步确认。本文采集点属淇河(卫河水系),该样点水质状况良好,水体清澈,流水处有细沙质地质,而静水处有20 cm厚淤泥,水体中有大量水草分布,该鱼在水草较丰富靠近岸边水流较缓处采集较多,这与其有筑巢行为有关。

3.3 林州市中华多刺鱼的来源及保护

在对河南省鱼类资源调查中,仅在林州一个采样点发现了中华多刺鱼,且种群数量有限,这和该处生境水质良好、有大量水草,适宜有特殊产卵习性的中华多刺鱼繁殖有关。同时作者也访问了附近居民及林州市渔政管理人员,他们均反映之前没有发现此种鱼类,也没有该种鱼的引入记录,这也许与该种鱼个体较小,没有太大经济价值有关,因此推测可能为自然分布。此后调查队数次在此地开展复查,均获得了数量不等的样本,说明中华多刺鱼在该地仍保持一定数量。关于该种鱼在林州的分布范围,在调查队扩大采集范围及安阳市水产科学研究所开展全市其他区县调查中,均没有采集到此种鱼类,因此不排除人为带入的可能性。但该种鱼类在该样点已经形成相对稳定的种群,且多次采样均获得该鱼,说明其可以适应该处的水文条件。同时调查队在该样点采集获得了包括中华多刺鱼在内的17种鱼类,该样点鱼类种类多样性较高。中华多刺鱼为北京市二级保护动物(http://www.bjny.gov.cn/nyj/231595/232732/605601/490633/index.html),同时,特殊的体型及繁殖习性让该鱼成为研究鱼类进化生物学、繁殖生物学、行为生态学的优良模型。

参考文献
卞绍雷, 赵亚辉, 张洁, 等. 2008. 温度和盐度对中华多刺鱼胚胎发育过程的影响[J]. 动物学报, 54(2): 282–289.
成庆泰, 郑葆珊. 1987. 中国鱼类系统检索[M]. 北京: 科学出版社.
成庆泰, 周才武. 1997. 山东鱼类志[M]. 济南: 山东科学技术出版社.
湖北省水生生物研究所鱼类研究室. 1976. 长江鱼类[M]. 北京: 科学出版社.
李明德, 杨竹舫. 1992. 河北省鱼类[M]. 北京: 海洋出版社.
李明德. 2011. 天津鱼类志[M]. 天津: 天津科学技术出版社.
刘婵馨, 秦克静. 1987. 辽宁动物志鱼类[M]. 沈阳: 辽宁科学技术出版社.
任慕莲. 1981. 黑龙江鱼类[M]. 哈尔滨: 黑龙江人民出版社.
天津市水产学会. 1990. 天津鱼类[M]. 北京: 海洋出版社.
王鸿媛. 1984. 北京鱼类志[M]. 北京: 北京出版社.
王所安, 王志敏, 李国良, 等. 2000. 河北动物志鱼类[M]. 石家庄: 河北科学技术出版社.
张觉民. 1995. 黑龙江省鱼类志[M]. 哈尔滨: 黑龙江科学技术出版社.
郑葆珊, 黄皓明, 张玉玲, 等. 1980. 图们江鱼类[M]. 长春: 吉林人民出版社.
J. 萨姆布鲁克, D. W. 拉塞尔, 黄培堂译. 2008. 分子克隆实验指南[M]. 北京: 科学出版社.
Г. B. 尼科尔斯基, 高岫译. 1960. 黑龙江流域鱼类[M]. 北京: 科学出版社.
Aldenhoven JT, Miller MA, Corneli PS, et al. 2010. Phylogeography of ninespine sticklebacks (Pungitius pungitius) in north America:glacial refugia and the origins of adaptive traits[J]. Molecular Ecology, 19: 4061–4076. DOI:10.1111/j.1365-294X.2010.04801.x
Chae BS, Yang HJ. 1993. Spawning and growth of eightspine stickleback, Pungitius sinensis kaibarae (Tanaka) (Gasterosteidae, Pisces), in the Chayang stream, a tributary of Kumho River, Korea[J]. Korean Journal of Zoology, 36: 181–192.
Chae BS, Yang HJ. 1988. Comparison of the morphology between two populations of eighspine stickleback, Pungitius sinensis kaibarae (Tanaka)[J]. Korean Journal of Limnology, 21: 79–91.
Chae BS. 1988. Classification on the eightspine stickleback (Pungitius sinensis: Gasterosteidae) to the morphological characteristics[D]. Daegu, South Korea: Kyungpook National University.
Denys GPJ, Geiger MF, Persat H, et al. 2015. Invalidity of Gasterosteus gymnurus (Cuvier, 1829) (Actinopterygii, Gasterosteidae) according to integrative taxonomy[J]. Cybium, 39(1): 37–45.
Folmer O, Black M, Hoeh W, et al. 1994. DNA primers for amplification of mitochondrial cytochrome c oxidase subunit Ⅰ from diverse metazoan invertebrates[J]. Molecular Marine Biology & Biotechnology, 3: 294–299.
Geiger MF, Herder F, Monaghan MT, et al. 2014. Spatial heterogeneity in the Mediterranean biodiversity hotspot affects barcoding accuracy of its freshwater fishes[J]. Molecular Ecology Resources, 14(6): 1210–1221. DOI:10.1111/men.2014.14.issue-6
Hubert N, Hanner R, Holm E, et al. 2008. Identifying Canadian freshwater fishes through DNA barcodes[J]. PLoS ONE, 3(6): e2490. DOI:10.1371/journal.pone.0002490
Hwang DS, Song HB, Lee JS. 2012. Complete mitochondrial genome of the Amur stickleback Pungitius sinensis (Gasterosteiformes, Gasterosteidae)[J]. Mitochondrial DNA, 23(4): 293–294. DOI:10.3109/19401736.2012.674125
Ki JS, Jung SO, Lee YM, et al. 2008. Unusual mitochondrial genome structures of the freshwater goby Odontobutis platycephala:rearrangement of tRNAs and an additional non-coding region[J]. Journal of Fish Biology, 73: 414–428. DOI:10.1111/j.1095-8649.2008.01911.x
Knebelsberger T, Dunz AR, Neumann D, et al. 2015. Molecular diversity of Germany's freshwater fishes and lampreys assessed by DNA barcoding[J]. Molecular Ecology Resources, 15(3): 562–572. DOI:10.1111/1755-0998.12322
Koizumi N, Jinguji H, Takahashi H, et al. 2007. Isolation and characterization of polymorphic microsatellite DNA markers in the Omono type of ninespine stickleback, genus Pungitius[J]. Molecular Ecology Notes, 7: 1315–1318. DOI:10.1111/men.2007.7.issue-6
Kumar S, Stecher G, Tamura K. 2016. MEGA7:molecular evolutionary genetics analysis version 7.0 for bigger datasets[J]. Molecular Biology & Evolution, 33: 1870–1874.
McCusker MR, Denti D, Van Guelpen L, et al. 2013. Barcoding Atlantic Canada's commonly encountered marine fishes[J]. Molecular Ecology Resources, 13(2): 177–188. DOI:10.1111/1755-0998.12043
Mecklenburg CW, Moller PR, Steinke D. 2011. Biodiversity of arctic marine fishes:taxonomy and zoogeography[J]. Marine Biodiversity, 41(1): 109–140. DOI:10.1007/s12526-010-0070-z
Rezansoff AM, Crispo E, Blair C, et al. 2015. Toward the genetic origins of a potentially non-native population of threespine stickleback (Gasterosteus aculeatus) in Alberta[J]. Conservation Genetics, 16(4): 1–15.
Shikano T, Shimada Y, Herczeg G, et al. 2010. History vs. habitat type:explaining the genetic structure of European nine-spined stickleback (Pungitius pungitius) populations[J]. Molecular Ecology, 19: 1147–1161. DOI:10.1111/mec.2010.19.issue-6
Takahashi H, Goto A. 2001. Evolution ofeast Asian ninespine sticklebacks as shown by mitochondrial DNA control region sequences[J]. Molecular Phylogenetics & Evolution, 21: 135–155.
Teacher AG, Shikano T, Karjalainen ME, et al. 2011. Phylogeography and genetic structuring of European nine-spined sticklebacks (Pungitius pungitius)-mitochondrial DNA evidence[J]. PLoS ONE: e19476. DOI:10.1371/journal.pone.0019476
Thalinger B, Oehm J, Mayr H, et al. 2016. Molecular prey identification in central European piscivores[J]. Molecular Ecology Resources, 16(1): 123–137. DOI:10.1111/1755-0998.12436
Watanabe K, Mori S, Nishida M. 2003. Genetic relationships and origin of two geographic groups of the freshwater threespine stickleback, "hariyo"[J]. Zoological Science, 20: 265–274. DOI:10.2108/zsj.20.265