Published May 18, 2020 | Version v1

Data from: Detecting aquatic invasive species in bait and pond stores with targeted environmental (e) DNA high-throughput sequencing metabarcode assays: angler, retailer, and manager implications

  • 1. National Oceanic and Atmospheric Administration
  • 2. University of Michigan–Ann Arbor
  • 3. University of Miami
  • 4. University of Toledo

Description

Bait and pond stores comprise potential, yet poorly understood, vectors for aquatic invasive species (AIS). We tested for AIS and illegal native species in 51 bait and 21 pond stores from the central Great Lakes (Lake Erie, Ohio and Lake St. Clair, Michigan) and the adjacent Wabash River (Indiana) using environmental DNA (eDNA) metabarcode assays of water samples and morphological identifications. Retailers were questioned about supply chains, and anglers surveyed about baitfish use and disposal. Assays revealed unadvertised species eDNA in 100% of bait stores, with 61% containing illegal native non-bait (totaling 13 species) and 88% having AIS (11 species). Illegal native non-bait species included juvenile walleye, yellow perch, and white sucker eDNA. AIS eDNA included Eurasian ruffe in seven stores (all states), silver carp in five (including a Lake Erie store in two separate years), and bighead carp in two Lake Erie stores that also had silver carp. Among pond stores, two in Lake St. Clair had bighead carp eDNA, one also contained silver carp, and a Wabash River location showed European ide. Unadvertised invasive snails were discerned in 55% of pond stores. Four contained zebra mussel eDNA and two had invasive bryozoans. Illegal native species and AIS were widespread, but showed little relationship to the retailers' variable and extensive supply chains. Live baitfish releases were reported by 50% of Lake Erie anglers and 35% in Lake St. Clair. Consumer behavior and AIS prevalence in the bait and pond trades thus pose serious risks for introductions and spread.

Notes

MetaTrim.py: See README.md

SeqTabToFasta.pl: script converts SeqTab format files with results from Dada2 (R package) to FASTA files for each sample. Expects an input file named SeqTab.txt in the CWD.

BlastCylce500.pl: script performs BLAST searches for every sequence in files ending with ".fasta" in the CWD. BLAST searches are performed against custom databases for all Actinopterygii on GenBank and only fishes in the Great Lakes. If file name begins with "m" (MiFish), sequences are blasted against 12S databases, otherwise, against Cytochrome b databases.

BlastCycle500Remote.pl: script performs BLAST searches for every sequence in files ending with ".fasta" in the CWD. BLAST is performed remotely against the entirety of GenBank.

SummarizeBlasteDNAEnvironmentSepASVs.pl: script summarizes BLAST results (files ending in "Res.txt" from BlastCycle500.pl"). Species names for hits with >90% query cover and identity are concatenated if the e value is the lowest for all hits. Returns a txt file with summarized hits per sample. Expects an input file named InReads.txt in the CWD. This file is tab delimited. The first value should be the sample name identical to that used as input in the above script, followed by tab, then the number of reads in that sample. This allows calculation of proportion of read represented by the hit.

BOPPSummarize*.py: script removes species hits that do not occur at a portion above the error for a specific marker in a specific sample, unless the occur in multiple markers in that same sample. Hits are then compared against morphological results from bait and pond stores. Several inputs are required and are included in this repository. They include a file showing the error rate calculated from positive controls for each marker on each sequencing run (ErrorCons.txt; if none is supplied, the default value will be 0.001), a shop key file (ShopKeybaitShops.txt), a file indicating the legality of species in each state sampled (SppStatusBait.txt), morphological results from each store at the time of sampling (InMorph.txt), and BLAST results files (summarized by the script above, separated by marker).

Funding provided by: Environmental Protection Agency
Crossref Funder Registry ID: http://dx.doi.org/10.13039/501100001589

Funding provided by: National Science Foundation
Crossref Funder Registry ID: http://dx.doi.org/10.13039/100000001

Files

ErrorCons.txt

Files (106.2 kB)

Name Size Download all
md5:7c02690b4ff0574f3251305980ff0e37
2.1 kB Download
md5:f6c53d77d25d2e42f9eed3946eca0658
835 Bytes Download
md5:71ba4d5dff102af89f6b7dda841e0dcb
16.8 kB Download
md5:5d5e24f88cffe62a3a340fa3e3d64099
6.1 kB Download
md5:be97cd0d880d81dfc2c4e1a76c6c83d8
8.5 kB Download
md5:c2f8607f6b9e0840bd4f0ec58c4ed60c
5.3 kB Preview Download
md5:1f6523066f82d39bf35e8ac161dc4de6
7.2 kB Preview Download
md5:9d01670bced4a78cdc0f80e3d890f628
29.5 kB Download
md5:05b0ec8bb742d87a7503802cf82e3373
7.0 kB Preview Download
md5:01a27e0caf97f2a3472b07f71e874817
1.5 kB Download
md5:9d53f74ca2364d3938f7621901c395cb
1.5 kB Download
md5:80af96cf2c67b96cbd920a38c3d9a4a4
3.5 kB Preview Download
md5:47254f645b18f3dde6866a5626e1788d
410 Bytes Preview Download
md5:3923c2b55cc059c3c3938beca3ae6b9e
3.3 kB Preview Download
md5:b0ca31b613f9737874f7972e3a7df766
3.8 kB Preview Download
md5:a7f71f7cab7e48b435a0fb13eaf2f5cf
8.9 kB Download

Additional details

Related works