Issue
Knowl. Manag. Aquat. Ecosyst.
Number 421, 2020
Topical Issue on Fish Ecology
Article Number 14
Number of page(s) 5
DOI https://doi.org/10.1051/kmae/2020006
Published online 03 April 2020

© K.A. Tkachuk and D.A. Dunn, Published by EDP Sciences 2020

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License CC-BY-ND (https://creativecommons.org/licenses/by-nd/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. If you remix, transform, or build upon the material, you may not distribute the modified material.

Petromyzon marinus, sea lampreys of the family Petromyzontidae are an invasive semelparous species within the Great Lakes basin of North America. Sea lampreys were first observed in Lake Ontario in the mid-1800s and fully invaded the Great Lakes basin through the early to mid 20th century (Christie and Goddard, 2003; Larson et al., 2003).

The Great Lakes Fishery Commission (GLFC) was established in the 1950's to monitor sea lamprey populations and develop control methods because mortality caused by parasitic lampreys had devastating effects on populations of lake trout, Salvelinus namaycush, and salmonids among other host fish (Fetterolf, 1980; Christie and Goddard, 2003; Larson et al., 2003). Since then, sea lamprey larvae have been identified in 66 tributaries of Lake Ontario (Adair and Sullivan, 2013; Hansen et al., 2016).

Because hatched larvae (ammocoetes) burrow into sandy/silty locations for two to ten-plus years, before metamorphosing into the parasitic phase of their life cycle (Dawson et al., 2015; Hansen et al., 2016), current detection methods are expensive, labor intensive, and can be disruptive to vulnerable populations of several species including native lamprey (Gingera et al., 2016). Methods include direct lamprey counting via electrofishing, trapping and application of granular bayluscide to tributaries (Slade et al., 2003; Steeves et al., 2003). Methods for population control can also be expensive and disruptive. For example, the lampricide 3-trifluoromethyl-4-nitrophenol (TFM) is a mitochondrial uncoupler (Birceanu et al., 2014) and is a primary control method for sea lampreys (Brege et al., 2003; Christie and Goddard, 2003). Acute TFM toxicity is known to occur in several non-target species found in TFM-treated waters including the threatened lake sturgeon, Acipenser fulvescens (Boogaard et al., 2003). In efforts to minimize off-target toxicity, the distribution, and concentration of TFM treatments has significantly decreased in tributaries (Brege et al., 2003). The need to limit TFM treatments in waterways not inhabited by sea lamprey populations, and accurately determine lamprey abundance highlights the importance of developing non-invasive, inexpensive and accurate methods for detection of lamprey populations.

Analysis of environmental DNA (eDNA) is a molecular technique that can be effective in detecting relatively low DNA copy number from environmental samples such as soil, air and water from a wide range of species (Dougherty et al., 2016; Gingera et al., 2016). The ease, comparative low cost and non-invasive nature of conducting eDNA analysis can provide advantages over traditional population survey methods. In recent studies, eDNA techniques were used to detect the invasive rusty crayfish, Orconectes rusticus, in lakes within Gogebic County, Michigan and Vilas County, Wisconsin, USA. These efforts proved to be more successful than trapping detection methods (Dougherty et al., 2016). eDNA was also shown to be more effective and sensitive in determining the threat of invasion of Asian carp species such as the silver carp, Hypophthalmichthys molitrix, and the bighead carp, Hypophthalmichthys nobilis, from the Mississippi River basin into the Great Lakes basin (Jerde et al., 2011). Specific amplification of sea lamprey DNA has been performed from water samples collected from multiple sites across the Great Lakes basin (Gingera et al., 2016). Since then, quantitative PCR (qPCR) analyses have been developed for sea lamprey that offer greater sensitivity (Bracken et al., 2019; Schloesser et al., 2018). With methodological improvements, these techniques are under consideration for implementation as a standard practice for determining lamprey abundance and distribution (Sullivan and Mullett, 2018). Nevertheless, eDNA analysis using end-point PCR is an inexpensive surveying technique that has the advantage of being relatively simple and accessible in many laboratories world-wide. The objective of this study is to evaluate the usefulness of a simple endpoint PCR assay to detect changes in lamprey environmental DNA before and after TFM treatment.

Ninemile Creek, a tributary to Lake Ontario which flows through parts of Cayuga and Oswego Counties, New York, USA, was treated with TFM in late May and early June of 2017. TFM treatment was conducted by employees of the Sea Lamprey Control Program of the GLFC administered by Fisheries and Oceans Canada and the United States Fish and Wildlife Service. Water samples were collected at four locations along Ninemile Creek (Fig. 1) on four separate dates: May 29, 2017 (early morning, pre-TFM treatment), June 2, 2017 (post-TFM treatment), September 22, 2017, and on June 9, 2018 one year following TFM treatment. Three water samples and a negative control sample consisting of commercially bottled drinking water were taken each time a sample was collected. Negative controls were collected and processed using the same protocol as the experimental samples to rule out on-site contamination.

Water sampling and DNA isolation were performed using the methods of Dougherty et al. (2016). Each water sample was collected in 250 mL sterile, individually wrapped Nalgene Analytical Test Filter Funnels (Thermo Fisher Scientific, Waltham, MA, USA). Each of the samples were filtered in the location where they were collected using 47 mm diameter cellulose nitrate (CN) membranes with 0.45 µm pores, a 1 L Nalgene Polypropylene Vacuum Flask (Thermo Fisher) and a CP7839 Hand Vacuum Pump (Actron, Warren, MI, USA). Each filtered CN membrane was removed using clean forceps and placed in 2 mL Eppendorf DNA LoBind microcentrifuge tubes (Thermo Fisher) containing 700 mL Longmire's buffer and placed in a dark, cold icebox (Longmire et al., 1997). Following collection, microcentrifuge tubes were stored in the laboratory at 4 °C until DNA isolation.

DNA isolation from CN membranes was performed via chloroform/isoamyl alcohol extraction followed by isopropyl alcohol precipitation. DNA pellets were resuspended in 100 µl TE with low EDTA (10mM Tris, 0.1mM EDTA, pH 8.0).

Mitochondrial cytochrome oxidase C subunit I gene (COI) PCR primers were used in all analyses. Primer sequences and PCR conditions used were according to Gingera et al. (2016).

Forward: 5'-GGCAACTGACTTGTACCCCTAATACTTGGT-3'

Reverse: 5'-GGCTAAGTGTAAGGAAAAGATTGTTAGGTCGAC-3'

Positive control DNA was obtained by performing PCR amplification of P. marinus genomic DNA using the above primers yielding a 225 bp product. Amplified DNA was gel isolated and quantified using a NanoDrop 2000 spectrophotometer (Thermo Fisher). Molecular mass of the expected 225bp DNA band was estimated (650 Da/bp) and serially diluted to 101 copies/µl. 5µl diluted stock (50 copies/reaction) were used as positive control template and was included in each experimental PCR.

The PCR performed for each eDNA extract consisted of 13.75 µl ultrapure H2O, 2.5 µl BioReady rTaq 10X Reaction Buffer (Bulldog Bio, Inc. Portsmouth, NH, USA), 0.25 µl dNTPs (20mM), 2.5 µl primers (2.0 µM each), 1 µl (5U) BioReady rTaq (Bulldog), 5 µl template DNA in a total reaction volume of 25 µl. Cycling conditions were 95° C for 5 min followed by 55 cycles of 95 °C for 30 s, 55 °C for 30 s, 72 °C for 30 s; 72 °C for 10 min with a terminal 4 °C hold. PCR products (expected band = 225 bp) were viewed on a 1.5% agarose gel stained with ethidium bromide. Each eDNA sample was subjected to a single PCR reaction.

All negative control samples collected during all four collection time points at each site were negative for sea lamprey DNA (n = 16). Seven and six total positive samples were found from the pre-TFM treatment time point and the immediate post-TFM treatment time point, respectively. All samples collected 16 weeks and one year post-TFM treatment from each site were negative for lamprey DNA (Tab. 1). Specificity of lamprey DNA amplification was analyzed via DNA sequencing of two positive PCR products (29-May-2017, site 4, sample 1; 2-June-2017, site 1, sample 1). 225bp PCR products were Sanger sequenced with both the forward and reverse PCR primers listed above. All sequence results aligned 100% with GenBank accession no. KX145312 (Petromyzon marinus COI).

Prior to the May 2017 TFM treatment, Ninemile Creek was treated in May 2014. An assessment of larval populations in the tributary conducted in August, 2016 estimated that residual larval lamprey − survivors from the 2014 treatment − and new larval populations were around 75,000 total larvae (Mullett and Sullivan, 2017). Ninemile Creek was subsequently treated with 182.4 kg liquid TFM, 2.0 kg Solid TFM and 0.2 kg granular bayluscide over a distance of 25.4 km in May, 2017. Our ability to detect sea lamprey presence pre-TFM treatment is consistent with this report on the state of sea lamprey infestation in Ninemile Creek at this timepoint. Our sampling sites spanned an approximate distance of 11.67 km from the opening of the tributary into Lake Ontario to Sample Site 1, the most upstream collection site. All sample sites fell within the range of the tributary that was treated. Exposure to TFM usually occurs over 12 hours for maximum lethality (Clifford et al., 2012) and TFM concentrations for 2017 treatments were designed to eliminate 93% of the larvae in the population (Sullivan and Mullett, 2018).

Our detection of lamprey DNA four days post treatment in three of the four sample sites might be interpreted as detection of residual live larval lamprey. More probably, this result might reflect detection of DNA from decomposing larval carcasses. This finding is in line with the results of Gingera et al., (2016) where sea lamprey eDNA signals were detectable for several weeks following spawning and subsequent death of adults in the absence of TFM treatment. Future work is warranted to correlate TFM treatment efficacy, flow rate, weather, and other environmental parameters with temporal resolution of eDNA detection. Our finding of not detecting lamprey DNA 16 weeks after TFM treatment is evidence of decrease in sea lamprey abundance to below a detectable level. This result alone, however does not yield precision data on the proportion of larval death. Caution should be taken as limitations in sensitivity of this assay led Gingera et al. (2016) to find false negatives in some lamprey infested waters. Nevertheless, this result agrees with the GLFC 2017 report in which larval populations assessed by electrofishing before and after treatment showed significant TFM-related mortality (Sullivan and Mullett, 2018). The August 2017 GLFC assessment found no evidence of residual populations and no recruitment was evident.

Sea lampreys within the Great Lakes usually enter the spawning phase of their life cycles after 12–18 months of parasitizing lake fish species (Hansen et al., 2016). Starting in March, spawning lamprey travel into tributaries, often following a pheromone secreted by larval lampreys. Between March and July, lampreys spawn, and shortly thereafter die (Potter et al., 2015). Our inability to detect lamprey DNA shortly following a spawning season (September, 2017) and a year following treatment (June, 2018) is further indication of population decline following TFM treatment. Since measurable larval populations were eliminated by August 2017 and no recruitment was evident (Sullivan and Mullett, 2018), pheromone-induced attraction of spawning adults to the creek may have been negligible (Bjerselius et al., 2000; Hansen et al., 2016). This may have contributed to our finding of no detectable lamprey DNA one year following TFM treatment. Overall, the results of this study suggest that there was a decline in larval populations post-TFM treatment that was detectable using eDNA surveying techniques.

Here, we took a simple endpoint PCR approach to eDNA analysis for presence/absence of a signal, with species specificity confirmed via DNA sequencing. Other techniques show promise for analytical improvements: Next-generation DNA sequencing/metabarcoding of environmental samples might provide insight into changes in population dynamics of multiple taxa in a single assay or in response to changing weather patterns or other environmental conditions (Valentini et al., 2015; Cannon et al., 2017). Likewise, qPCR assays could offer greater sensitivity (Bracken et al., 2019; Schloesser et al., 2018). Nevertheless, a major value of endpoint PCR is its simplicity, low cost, and availability in a wide range of academic and governmental laboratories. These factors make this technique a robust option for preliminary research and management-based molecular ecological monitoring, particularly in areas of low infrastructure. This study provides evidence for the utility of simple environmental DNA analysis in assessing population decline in aquatic environments.

Acknowledgments

We thank Ayla Busch for assistance in PCR analyses, as well as Seth Dunn and Joel Dunn for assistance in collecting samples. We are grateful to Bruce Morrison at Fisheries and Oceans Canada for communication and coordination regarding dates of TFM treatment at Ninemile Creek. We thank Brian Lantry at the United States Geological Survey Lake Ontario Biological Station for providing a sea lamprey specimen from which DNA was isolated to serve as positive control. Research was supported by a grant from Rice Creek Associates.

References

  • Adair R, Sullivan P. 2013. Sea lamprey control in the Great Lakes 2012, Great Lake Fishery Commission, Ann Arbor, MI USA, Available from http://www.glfc.org/pubs/slcp/annual_reports/ANNUAL_REPORT_2012.pdf (Accessed January 2020) [Google Scholar]
  • Birceanu O, Sorensen LA, Henry M, McClelland GB, Wang YS, Wilkie MP. 2014. The effects of the lampricide 3-trifluoromethyl-4-nitrophenol (TFM) on fuel stores and ion balance in a non-target fish, the rainbow trout (Oncorhynchus mykiss). Comp Biochem Physiol Part C Toxicol Pharmcol 160: 30–41. [CrossRef] [Google Scholar]
  • Bjerselius R, Weiming L, Teeter JH, et al. 2000. Direct behavioral evidence that unique bile acids released by larval sea lamprey (Petromyzon marinus) function as a migratory pheromone. Can J Fish Aquat Sci 57: 557–569. [Google Scholar]
  • Boogaard MA, Bills TD, Johnson DA. 2003. Acute toxicity of TFM and a TFM/Niclosamide mixture to selected species of fish, including lake sturgeon (Acipenser fulvescens) and mudpuppies (Necturus maculosus), in laboratory and field exposures. J Great Lakes Res 29: 529–541. [Google Scholar]
  • Bracken FSA, Rooney SM, Kelly-Quinn M, King JJ, Carlsson J. 2019. Identifying spawning sites and other critical habitat in lotic systems using eDNA “snapshots”: a case study using the sea lamprey Petromyzon marinus L. Ecol Evol 9: 553–567. [PubMed] [Google Scholar]
  • Brege DC, Davis DM, Genovese JH, McAuley TC, Stephens BE, Westman WR. 2003. Factors responsible for the reduction in quantity of the lampricide, TFM, applied annually in streams tributary to the Great Lakes from 1979 to 1999. J Great Lakes Res 29: 500–509. [Google Scholar]
  • Cannon MV, Craine J, Hester J, et al. 2017. Dynamic microbial populations along the Cuyahoga River. PLoS ONE 12: e0186290. [CrossRef] [PubMed] [Google Scholar]
  • Christie GC, Goddard CI. 2003. Sea Lamprey International Symposium (SLIS II): advances in the integrated management of sea lamprey in the Great Lakes. J Great Lakes Res 29: 1–14. [Google Scholar]
  • Clifford AM, Henry M, Bergstedt R, McDonald DG, Smits AS, Wilkie MP. 2012. Recovery of larval sea lampreys from short-term exposure to the pesticide 3-trifluoromethyl-4-nitrophenol: implications for sea lamprey control in the Great Lakes. Trans Am Fish Soc 141: 1697–1710. [Google Scholar]
  • Dawson HA, Quintella BR, Almeida PR, Treble AJ, Jolley JC. 2015. Ecology of larval and metamorphosing lampreys. In Docker MF, ed. Lampreys: Biology, Conservation and Control, Springer, Dordrecht, 75–137. [Google Scholar]
  • Dougherty MM, Larson ER, Renshaw MA, et al. 2016. Environmental DNA (eDNA) detects the invasive rusty crayfish Orconectes rusticus at low abundances. J Appl Ecol 53: 722–732. [PubMed] [Google Scholar]
  • Fetterolf CM. 1980. Why a Great Lakes Fishery Commission and why a Sea Lamprey International Symposium. Can J Fish Aquat Sci 37: 1588–1593. [Google Scholar]
  • Gingera TD, Steeves TB, Boguski DA, Whyard S, Li W, Docker MF. 2016. Detection and identification of lampreys in Great Lakes streams using environmental DNA. J Great Lakes Res 42: 649–659. [Google Scholar]
  • Hansen MJ, Madenjian CP, Slade WJ, Steeves TB, Almeida PR, Quintella BR. 2016. Population ecology of the sea lamprey (Petromyzon marinus) as an invasive species in the Laurentian Great Lakes and an imperiled species in Europe. Rev Fish Biol Fish 26: 509–535. [Google Scholar]
  • Jerde CL, Mahon AR, Chadderton WL, Lodge DM. 2011. “Sight-unseen” detection of rare aquatic species using environmental DNA: eDNA surveillance of rare aquatic species. Conserv Lett 4: 150–157. [Google Scholar]
  • Larson GL, Christie GC, Johnson DA, Koonce JF, Mullett KM, Sullivan WP. 2003. The history of sea lamprey control in Lake Ontario and updated estimates of suppression targets. J Great Lakes Res 29: 637–654. [Google Scholar]
  • Longmire JL, Maltbie M, Baker, RJ. 1997. Use of “lysis buffer” in DNA isolation and its implication for museum collections. Spec Publ Mus Tex Tech Univ 163: 1–3. [Google Scholar]
  • Mullett K, Sullivan P. 2017. Sea lamprey control in the Great Lakes 2016, Great Lakes Fishery Commission, Ann Arbor, MI USA, Available from http://www.glfc.org/pubs/slcp/annual_reports/ANNUAL_REPORT_2016.pdf (accessed January 2020). [Google Scholar]
  • Potter IC, Gill HS, Renaud CB, Haoucher D. 2015. The taxonomy, phylogeny, and distribution of lampreys. In: Docker M (Ed.), Lampreys: Biology, Conservation and Control. Dordrecht: Springer, 35–73. [CrossRef] [Google Scholar]
  • Schloesser NA, Merkes CM, Rees CB, Amberg JJ, Steeves TB, Docker MF. 2018. Correlating sea lamprey density with environmental DNA detections in the lab. Manag Biol Invasions 9: 483–95. [CrossRef] [Google Scholar]
  • Slade JW, Adams JV, Christie GC, et al. 2003. Techniques and methods for estimating abundance of larval and metamorphosed sea lampreys in Great Lakes tributaries, 1995 to 2001. J Great Lakes Res 29: 137–151. [Google Scholar]
  • Steeves TB, Slade JW, Fodale MF, Cuddy DW, Jones ML. 2003. Effectiveness of using backpack electrofishing gear for collecting sea lamprey (Petromyzon marinus) larvae in Great Lakes tributaries. J Great Lakes Res 29: 161–173. [Google Scholar]
  • Sullivan P, Mullett K. 2018. Sea lamprey control in the Great Lakes 2017, Great Lakes Fishery Commission, Ann Arbor, MI USA, Available from http://www.glfc.org/pubs/slcp/annual_reports/ANNUAL_REPORT_2017.pdf (accessed January 2020) [Google Scholar]
  • Valentini A, Taberlet P, Miaud C, et al. 2015. Next-generation monitoring of aquatic biodiversity using environmental DNA metabarcoding. Mol Ecol 25: 929–942. [Google Scholar]

Cite this article as: Tkachuk KA, Dunn DA. 2020. Analysis of sea Lamprey environmental DNA during Lampricide treatment in a tributary of Lake Ontario. Knowl. Manag. Aquat. Ecosyst., 421, 14.

All Tables

Table 1

Detection of sea lamprey DNA at corresponding sampling locations and dates.

All Figures

thumbnail Fig. 1

Sampling sites along Ninemile Creek.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.