Acknowledgments:
We thank Hokkaido Aquaculture Promotion Corporation, Michio Fukushima and Hiroki Mizumoto for providing a part of the water samples, Fumito Tashiro for providing a tissue sample, Ryukoku University Center for Biodiversity Science for analyzing a part of the eDNA metabarcoding samples, and Kazusa DNA Research Institute for DNA extraction of a part of the eDNA samples. We also thank Akira Goto and the members in Animal Ecology lab in Faculty of Agriculture, Hokkaido University for many suggestions in the development of this study. We used the DNA Sequencing Facility of Research Faculty of Agriculture, Hokkaido University.
Funding:
Environment Research and Technology Development Fund of the Environmental Restoration and Conservation Agency Provided by the Ministry of Environment of Japan grant no. JPMEERF20204004 (HA)
JSPS KAKENHI grant no. 20H03005 (HA)
JSPS KAKENHI grant no. 23H00329 (HA)
Grant-in-Aid for JSPS Fellow grant no. JP21J12688 (TY)
Author contributions:Conceptualization: TY, TK, HA Methodology: TY, TK, SI Investigation: TY, KF, HA Visualization: TY Funding acquisition: TY, HA Project administration: TY, HA Supervision: HA Writing – original draft: TY, TK, KF, SI, HA Writing – review & editing: TY, TK, KF, SI, HAData availability:
Raw sequence reads are deposited in DDBJ Sequence Read Archive (DRA). For the MiFish sequence data, accession numbers are DRR440535–DRR440728 and DRR488036–DRR488183. For Barbatula toni ’s cyt-bsequence data, accession numbers are DRR490086–DRR490140. Full details of the methods, results, supporting tables and figures are available in Supplementary Information.