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.
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.