Health of wild fish exposed to contaminants of emerging concern in freshwater ecosystems utilized by a Minnesota Tribal community

Fish serve as indicators of exposure to contaminants of emerging concern (CECs)-chemicals such as pharmaceuticals, hormones, and personal care products-which are often designed to impact vertebrates. To investigate fish health and CECs in situ, we evaluated the health of wild fish exposed to CECs in waterbodies across northeastern Minnesota with varying anthropogenic pressures and CEC exposures: waterbodies with no human development along their shorelines, those with development, and those directly receiving treated wastewater effluent. Then, we compared three approaches to evaluate the health of fish exposed to CECs in their natural environment: a refined fish health assessment index, a histopathological index, and high-throughput (ToxCast) in vitro assays. Lastly, we mapped adverse outcome pathways (AOPs) associated with identified ToxCast assays to determine potential impacts across levels of biological organization within the aquatic system. These approaches were applied to subsistence fish collected from the Grand Portage Indian Reservation and 1854 Ceded Territory in 2017 and 2019. Overall, 24 CECs were detected in fish tissues, with all but one of the sites having at least one detection. The combined implementation of these tools revealed that subsistence fish exposed to CECs had histological and macroscopic tissue and organ abnormalities, although a direct causal link could not be established. The health of fish in undeveloped sites was as poor, or sometimes poorer, than fish in developed and wastewater effluent-impacted sites based on gross and histologic tissue lesions. Adverse outcome pathways revealed potential hazardous pathways of individual CECs to fish. A better understanding of how the health of wild fish harvested for consumption is affected by CECs may help prioritize risk management research efforts and can ultimately be used to guide fishery management and public health decisions. Integr Environ Assess Manag 2024;20:846-863. © 2023 The Authors. Integrated Environmental Assessment and Management published by Wiley Periodicals LLC on behalf of Society of Environmental Toxicology & Chemistry (SETAC).

Modeling the risk of aquatic species invasion spread through boater movements and river connections

Aquatic invasive species (AIS) are one of the greatest threats to the functioning of aquatic ecosystems worldwide. Once an invasive species has been introduced to a new region, many governments develop management strategies to reduce further spread. Nevertheless, managing AIS in a new region is challenging because of the vast areas that need protection and limited resources. Spatial heterogeneity in invasion risk is driven by environmental suitability and propagule pressure, which can be used to prioritize locations for surveillance and intervention activities. To better understand invasion risk across aquatic landscapes, we developed a simulation model to estimate the likelihood of a waterbody becoming invaded with an AIS. The model included waterbodies connected via a multilayer network that included boater movements and hydrological connections. In a case study of Minnesota, we used zebra mussels (Dreissena polymorpha) and starry stonewort (Nitellopsis obtusa) as model species. We simulated the impacts of management scenarios developed by stakeholders and created a decision-support tool available through an online application provided as part of the AIS Explorer dashboard. Our baseline model revealed that 89% of new zebra mussel invasions and 84% of new starry stonewort invasions occurred through boater movements, establishing it as a primary pathway of spread and offering insights beyond risk estimates generated by traditional environmental suitability models alone. Our results highlight the critical role of interventions applied to boater movements to reduce AIS dispersal.

Open water dreissenid mussel control projects: lessons learned from a retrospective analysis

Dreissenid mussels are one of the most problematic aquatic invasive species (AIS) in North America, causing substantial ecological and economic effects. To date, dreissenid mussel control efforts in open water have included physical, biological, and chemical methods. The feasibility of successful dreissenid mussel management or eradication in lakes is relatively undocumented in the freshwater management literature. This review presents information on 33 open water dreissenid mussel control projects in 23 North America lakes. We reviewed data from past dreissenid mussel control projects and identified patterns and knowledge gaps to help inform adaptive management strategies. The three key lessons learned include (1) pre- and post-treatment survey methods that are designed to meet management objectives are beneficial, e.g., by sampling for all life stages and taking into account that no survey method is completely comprehensive; (2) defining the treatment area-particularly ensuring it is sufficiently large to capture all life stages present-is critical to meeting management objectives; and (3) control projects provide an opportunity to collect water chemistry, effects on non-target organisms, and other efficacy-related data that can inform safe and effective adaptive management.

Release of live baitfish by recreational anglers drives fish pathogen introduction risk

Emerging diseases of wildlife are an existential threat to biodiversity, and human-mediated movements of live animals are a primary vector of their spread. Wildlife disease risk analyses offer an appealing alternative to precautionary approaches because they allow for explicit quantification of uncertainties and consideration of tradeoffs. Such considerations become particularly important in high-frequency invasion pathways with hundreds of thousands of individual vectors, where even low pathogen prevalence can lead to substantial risk. The purpose of this study was to examine the landscape-level dynamics of human behavior-mediated pathogen introduction risk in the context of a high-frequency invasion pathway. One such pathway is the use and release of live fish used as bait by recreational anglers. We used a stochastic risk assessment model parameterized by angler survey data from Minnesota, USA, to simulate one year of fishing in Minnesota and estimate the total number of risky trips for each of three pathogens: viral hemorrhagic septicemia virus, the microsporidian parasite Ovipleistophora ovariae, and the Asian fish tapeworm Schizocotyle acheilognathi. We assessed the number of introductions under four scenarios: current/baseline conditions, outbreak conditions (increased pathogen prevalence), source-focused control measures (decreased pathogen prevalence), and angler-focused control measures (decreased rates of release). We found that hundreds of thousands of introduction events can occur per year, even for regulated pathogens at low pathogen prevalence. Reducing the rate of illegal baitfish release had significant impact on risky trips in scenarios where a high number of anglers were involved, but was less impactful in circumstances with limited outbreaks and fewer affected anglers. In contrast, reducing pathogen prevalence in the source populations of baitfish had relatively little impact. In order to make meaningful changes in pathogen introduction risk, managers should focus efforts on containing local outbreaks and reducing illegal baitfish release to reduce pathogen introduction risk. Our study also demonstrates the risk associated with high-frequency invasion pathways and the importance of incorporating human behaviors into wildlife disease models and risk assessments.

Disease ecology and host range of Cyprinid herpesvirus 3 (CyHV-3) in CyHV-3 endemic lakes of North America

Cyprinid herpesvirus-3 (CyHV-3) is an important pathogen of common carp (Cyprinus carpio, carp) causing significant economic and ecological impacts worldwide. The recent emergence of CyHV-3 in the Upper Midwest region of the United States has raised questions related to the disease ecology and host specificity of CyHV-3 in wild carp populations. To determine the prevalence of CyHV-3 in wild populations of fishes in Minnesota, we surveyed five lakes in 2019 in which the virus was known to have caused mass mortality events in carp from 2017 to 2018. A total of 28 species (n = 756 total fish) of native fishes and 730 carp were screened for the presence of CyHV-3 DNA using specific qPCR. None of the native fish tissues tested positive for CyHV-3 although the prevalence of CyHV-3 in carp was 10%-50% in the five lakes. A single lake (Lake Elysian) with a 50% DNA detection rate and evidence of ongoing transmission and CyHV-3-associated mortality was surveyed again in 2020 from April to September. During this period, none of the tissues from 24 species (n = 607 total fish) tested positive for CyHV-3 though CyHV-3 DNA and mRNA (indicating viral replication) was detected in carp tissues during the sampling period. CyHV-3 DNA was detected most often in brain samples without evidence of replication, potentially indicating that brain tissue is a site for CyHV-3 latency. Paired qPCR and ELISA testing for Lake Elysian in 2019-2020 identified young carp (especially males) to be the primary group impacted by CyHV-3-associated mortality and acute infections, but with no positive detections in juvenile carp. Seroprevalence of carp from Lake Elysian was 57% in 2019, 92% in April of 2020 and 97% in September 2020. These results further corroborate the host specificity of CyHV-3 to carp in mixed wild populations of fishes in Minnesota and provide additional insights into the ecological niche of CyHV-3 in shallow lake populations of carp in North America.

A bi-level model for state and county aquatic invasive species prevention decisions

Recreational boats are important vectors of spread of aquatic invasive species (AIS) among waterbodies of the United States. To limit AIS spread, state and county agencies fund watercraft inspection and decontamination stations at lake access points. We present a bi-level model for determining how a state planner can efficiently allocate inspection resources to county managers, who independently decide where to locate inspection stations. In our formulation, each county manager determines a set of optimal plans for the locations of inspection stations under various resource constraints. Each plan maximizes inspections of risky boats that may carry AIS from infested to uninfested lakes within the county. Then, the state planner selects the set of county plans (i.e., one plan for each county) that maximizes the number of risky boats inspected throughout the state subject to a statewide resource constraint. We apply the model using information from Minnesota, USA, including the infestation status of 9182 lakes and estimates of annual numbers of boat movements from infested to uninfested lakes. Comparison of solutions of the bi-level model with solutions of a state-level model where a state planner selects lakes for inspection stations statewide shows that when state and county objectives are not aligned, the loss in efficiency at the state-level can be substantial.

AIS explorer: Prioritization for watercraft inspections-A decision-support tool for aquatic invasive species management

Invasions of aquatic invasive species have caused significant economic and ecological damage to global aquatic ecosystems. Once an invasive population has established in a new habitat, eradication can be financially and logistically impossible, motivating management strategies to rely heavily upon prevention measures to reduce the introduction and spread. To be productive, on-the-ground management of aquatic invasive species requires effective decision-making surrounding the allocation of limited resources. Watercraft inspections play an important role in managing aquatic invasive species by preventing the overland transport of invasive species between waterbodies and providing education to boaters. In this study, we developed and tested an interactive web-based decision-support tool, AIS Explorer: Prioritization for Watercraft Inspections, to guide AIS managers in developing efficient watercraft inspection plans. The decision-support tool is informed by a network-based algorithm that maximized the number of inspected watercraft that move from AIS infested to uninfested lakes within and between counties in Minnesota, USA. It was iteratively built with stakeholder feedback, including consultations with county managers, beta-testing of the web-based application, and workshops to educate and train end-users. The co-development and implementation of data-driven decision support tools demonstrate how interdisciplinary methods can be used to connect science and management to support decision-making. The AIS Explorer: Prioritization for Watercraft Inspections application makes optimized research outputs accessible in multiple dynamic forms that maintain pace with discovery of new infestations and local needs. In addition, the decision support tool has supported improved and closer communication between AIS managers and researchers on this topic.

Investigation of Cyprinid Herpesvirus 3 (CyHV-3) Disease Periods and Factors Influencing CyHV-3 Transmission in A Low Stocking Density Infection Trial

Simple Summary

Pathogens are the primary limitation to aquaculture production of fish and a major issue in consideration of the interface between cultured and wild populations of fishes worldwide. While rapid spread of fish pathogens between populations (wild or farmed) is generally anthropogenic and the result of trade, the mechanisms of transmission once a pathogen has been introduced to a fish population are not well understood. The most widespread pathogen impacting both aquaculture and wild populations of common carp (Cyprinus carpio, carp) is Cyprinid herpesvirus 3 (CyHV-3). To understand how CyHV-3 is transmitted in a population we conducted a series of infection trials, designed to determine the kinetics CyHV-3 infections, identify the contributions of direct and indirect forms of CyHV-3 transmission, and to determine the contributions of contact rate, viral load, pathogenicity, and contact type. We found that direct contact between fish was the primary mechanism of CyHV-3 transmission rather than transmission through contaminated water. Additionally, CyHV-3 transmission occurred primarily during the incubation period of CyHV-3, prior to the appearance of disease signs and disease-associated reduction in contact rate.

Abstract

Cyprinid herpesvirus 3 (CyHV-3) is the etiological agent of koi herpesvirus disease (KHVD) and important pathogen of aquaculture and wild populations of common carp worldwide. Understanding the relative contributions of direct and indirect transmission of CyHV-3 as well as the factors that drive CyHV-3 transmission can clarify the importance of environmental disease vectors and is valuable for informing disease modeling efforts. To study the mechanisms and factors driving CyHV-3 transmission we conducted infection trials that determined the kinetics of KHVD and the contributions of direct and indirect forms of CyHV-3 transmission, as well as the contributions of contact rate, viral load, pathogenicity and contact type. The incubation period of KHVD was 5.88 + 1.75 days and the symptomatic period was 5.31 + 0.87 days. Direct transmission was determined to be the primary mechanism of CyHV-3 transmission (OR = 25.08, 95%CI = 10.73–99.99, p = 4.29 × 10⁻¹⁸) and transmission primarily occurred during the incubation period of KHVD. Direct transmission decreased in the symptomatic period of disease. Transmissibility of CyHV-3 and indirect transmission increased during the symptomatic period of disease, correlating with increased viral loads. Additionally, potential virulence-transmission tradeoffs and disease avoidance behaviors relevant to CyHV-3 transmission were identified.

Optimizing the location of watercraft inspection stations to slow the spread of aquatic invasive species

The accidental spread of aquatic invasive species (AIS) by recreational boaters is a major concern of state and county environmental planners in the USA. While programs for watercraft inspection to educate boaters and slow AIS spread are common practice, large numbers of boats and waterbodies, together with limited budgets, make program design difficult. To facilitate program design, we developed an integer programming model for allocation of scarce inspection resources among lakes. Our model uses species-specific infestation status of lakes and estimates of boat movement between lakes. The objective is to select lakes for inspection stations to maximize the number of risky boats inspected, where risky boats are ones that move from infested to uninfested lakes. We apply our model in Stearns County in central Minnesota, USA, to prioritize lakes for inspection stations and evaluate alternative management objectives. With an objective of protecting uninfested lakes within and outside Stearns County, the optimal policy is to locate stations at infested lakes having the most boats departing for uninfested lakes inside and outside the county. With an objective of protecting only Stearns County lakes, the optimal policy is to locate stations at both infested and uninfested lakes having the riskiest boats arriving from within and outside the county and departing to in-county lakes. The tradeoff between these objectives is significant.

Anthropogenic factors associated with contaminants of emerging concern detected in inland Minnesota lakes (Phase II)

Contaminants of emerging concern (CECs) include a variety of pharmaceuticals, personal care products, and hormones commonly detected in surface waters. Human activities, such as wastewater treatment and discharge, contribute to the distribution of CECs in water, but other sources and pathways are less frequently examined. This study aimed to identify anthropogenic activities and environmental characteristics associated with the presence of CECs, previously determined to be of high priority for further research and mitigation, in rural inland lakes in northeastern Minnesota, United States. The setting for this study consisted of 21 lakes located within both the Grand Portage Indian Reservation and the 1854 Ceded Territory, where subsistence hunting and fishing are important to the cultural heritage of the indigenous community. We used data pertaining to numbers of buildings, healthcare facilities, wastewater treatment plants, impervious surfaces, and wetlands within defined areas surrounding the lakes as potential predictors of the detection of high priority CECs in water, sediment, and fish. Separate models were run for each contaminant detected in each sample media. We used least absolute shrinkage and selection operator (LASSO) models to account for both predictor selection and parameter estimation for CEC detection. Across contaminants and sample media, the percentage of impervious surface was consistently positively associated with CEC detection. Number of buildings in the surrounding area was often negatively associated with CEC detection, though nonsignificant. Surrounding population, presence of wastewater treatment facilities, and percentage of wetlands in surrounding areas were positively, but inconsistently, associated with CECs, while catchment area and healthcare centers were generally not associated. The results of this study highlight human activities and environmental characteristics associated with CEC presence in a rural area, informing future work regarding specific sources and transport pathways. We also demonstrate the utility of LASSO modeling in the identification of these important relationships.

A chemical prioritization process: Applications to contaminants of emerging concern in freshwater ecosystems (Phase I)

Contaminants of emerging concern (CECs), such as pharmaceuticals, personal care products, and hormones, are frequently found in aquatic ecosystems around the world. Information on sublethal effects from exposure to commonly detected concentrations of CECs is lacking and the limited availability of toxicity data makes it difficult to interpret the biological significance of occurrence data. However, the ability to evaluate the effects of CECs on aquatic ecosystems is growing in importance, as detection frequency increases. The goal of this study was to prioritize the chemical hazards of 117 CECs detected in subsistence species and freshwater ecosystems on the Grand Portage Indian Reservation and adjacent 1854 Ceded Territory in Minnesota, USA. To prioritize CECs for management actions, we adapted Minnesota Pollution Control Agency's Aquatic Toxicity Profiles framework, a tool for the rapid assessment of contaminants to cause adverse effects on aquatic life by incorporating chemical-specific information. This study aimed to 1) perform a rapid-screening assessment and prioritization of detected CECs based on their potential environmental hazard; 2) identify waterbodies in the study region that contain high priority CECs; and 3) inform future monitoring, assessment, and potential remediation in the study region. In water samples alone, 50 CECs were deemed high priority. Twenty-one CECs were high priority among sediment samples and seven CECs were high priority in fish samples. Azithromycin, DEET, diphenhydramine, fluoxetine, miconazole, and verapamil were high priority in all three media. Due to the presence of high priority CECs throughout the study region, we recommend future monitoring of particular CECs based on the prioritization method used here. We present an application of a chemical hazard prioritization process and identify areas where the framework may be adapted to meet the objectives of other management-related assessments.

Field and Laboratory Evaluation of the Microsporidian Parasite Heterosporis sutherlandae: Prevalence, Severity, and Transmission

Heterosporis sutherlandae is an invasive microsporidian parasite in the Great Lakes region of North America that infects the skeletal muscle of numerous fish species, rendering the fillet unfit for human consumption. Although H. sutherlandae has been identified as a pathogen of concern by state management agencies, there is little information to inform regulation and intervention. We sampled fishes over 1 year from three lakes in northern Minnesota with known infected populations to determine the importance of host demographic and environmental variables for influencing H. sutherlandae infection prevalence. Heterosporis sutherlandae was present during all sampling periods, ranging in prevalence from 1% to 11%. The prevalence of H. sutherlandae among Yellow Perch Perca flavescens varied significantly according to season, with winter having the lowest prevalence (1%) and summer having the highest prevalence (11%). For other fish species, the prevalence of H. sutherlandae also varied significantly with season: the lowest prevalence occurred during spring (1%) and the highest prevalence occurred in fall (9%). Rates of pathogen transmission were estimated by exposing Fathead Minnows Pimephales promelas in the laboratory. Transmission rates were 23% when naïve fish were fed infected tissues and only 2% when naïve fish were held in cohabitation with tissue-fed fish. Exposure method and exposure duration (d) increased the probability that a fish was infected with H. sutherlandae. These findings suggest that H. sutherlandae transmission is greater when a susceptible host consumes infected tissue than when the fish is exposed to spores present in the water column. The current rates of infection in wild fishes are in stark contrast to the prevalence documented in 2004 (28%), suggesting a reduction in H. sutherlandae prevalence within at least one Yellow Perch population in the Laurentian Great Lakes region since the early 2000s.

Influence of Environmental Stressors on the Microbiota of Zebra Mussels (Dreissena polymorpha)

Host-associated microbiota play a critical role in host fitness by providing nutrition, enhancing digestion capabilities, and by providing protection from pathogens. Here, we investigated the effects of two environmental stressors, temperature, and salinity, on the microbiota associated with zebra mussels (ZMs), a highly invasive bivalve in North America. To examine this in detail, lake-collected ZMs were acclimated to laboratory conditions, and subjected to temperature and salinity stress conditions. The impact of these stressors on the diversity, composition, and dynamics of ZM-associated microbiota were assessed by using amplicon- and shotgun-based sequencing, and qPCR-based approaches. Elevated temperature was found to be the primary driver of ZM mortality, although salinity alone also increased its likelihood. Stressor-induced ZM mortality, which ranged between 53 and 100%, was concomitant with significant increases in the relative abundance of several genera of putative opportunistic pathogens including Aeromonas. These genera were only present in low relative abundance in ZMs obtained from the control tank with 0% mortality. Shotgun sequencing and qPCR analyses indicated that the relative and absolute abundances of pathogenic Aeromonas species (particularly A. veronii) were significantly greater in temperature-induced dead ZMs. Taken together, our results show that environmental stress, especially elevated temperature (> 25 °C), is associated with the rapid mortality of ZMs as well as the proliferation of putative opportunistic bacterial pathogens.

An expert-based risk ranking framework for assessing potential pathogens in the live baitfish trade

As global trade of live animals expands, there is increasing need to assess the risks of invasive organisms, including pathogens, that can accompany these translocations. The movement and release of live baitfish by recreational anglers has been identified as a particularly high‐risk pathway for the spread of aquatic diseases in the United States. To provide risk‐based decision support for preventing and managing disease invasions from baitfish release, we developed a hazard identification and ranking tool to identify the pathogens that pose the highest risk to wild fish via this pathway. We created a screening protocol and semi‐quantitative stochastic risk ranking framework, combining published data with expert elicitation (n = 25) and applied the framework to identify high‐priority pathogens for the bait supply in Minnesota, USA. Normalized scores were developed for seven risk criteria (likelihood of transfer, prevalence in bait supply, likelihood of colonization, current distribution, economic impact if established, ecological impact if established and host species) to characterize a pathogen's ability to persist in the bait supply and cause impacts to wild fish species of concern. The generalist macroparasite Schizocotyle acheilognathi was identified as presenting highest overall threat, followed by the microsporidian Ovipleistophora ovariae, and viral haemorrhagic septicaemia virus. Our findings provide risk‐based decision support for managers charged with maintaining both the recreational fishing industry and sustainable, healthy natural resources. Particularly, the identification of several high‐risk but currently unregulated pathogens suggests that focusing risk management on pathogens of concern in all potential host species could reduce disease introduction risk. The ranking process, implemented here for a single state case study, provides a conceptual framework for integrating expert opinion and sparse available data that could be scaled up and applied across jurisdictions to inform risk‐based management of the live baitfish pathway.

Occurrence of contaminants of emerging concern in aquatic ecosystems utilized by Minnesota tribal communities

Pharmaceuticals, personal care products, hormones, and other chemicals lacking water quality standards are frequently found in surface water. While evidence is growing that these contaminants of emerging concern (CECs) - those previously unknown, unrecognized, or unregulated - can affect the behavior and reproduction of fish and wildlife, little is known about the distribution of these chemicals in rural, tribal areas. Therefore, we surveyed the presence of CECs in water, sediment, and subsistence fish species across various waterbodies, categorized as undeveloped (i.e., no human development along shorelines), developed (i.e., human development along shorelines), and wastewater effluent-impacted (i.e., contain effluence from wastewater treatment plants), within the Grand Portage Indian Reservation and 1854 Ceded Territory in northeastern Minnesota, U.S.A. Overall, in 28 sites across three years (2016-2018), 117 of the 158 compounds tested were detected in at least one form of medium (i.e., water, sediment, or fish). CECs were detected most frequently at wastewater effluent-impacted sites, with up to 83 chemicals detected in one such lake, while as many as 17 were detected in an undeveloped lake. Although there was no statistically significant difference between the number of CECs present in developed versus undeveloped lakes, a range of 3-17 CECs were detected across these locations. Twenty-two CECs were detected in developed and undeveloped sites that were not detected in wastewater effluent-impacted sites. The detection of CECs in remote, undeveloped locations, where subsistence fish are harvested, raises scientific questions about the safety and security of subsistence foods for indigenous communities. Further investigation is warranted so that science-based solutions to reduce chemical risks to aquatic life and people can be developed locally and be informative for indigenous communities elsewhere.

Lessons Learned From the Stakeholder Engagement in Research: Application of Spatial Analytical Tools in One Health Problems

Stakeholder engagement in research is widely advocated as a tool to integrate diverse knowledge and perspectives in the management of health threats while addressing potential conflicts of interest. Although guidelines for stakeholder engagement exist in public health and environmental sciences, the feasibility of actionable decisions based on scientific analyses and the lessons learned from the stakeholder engagement in the process co-creation of knowledge have been rarely discussed in One Health literature and veterinary sciences. Risk maps and risk regionalization using spatiotemporal epidemiological/analytical tools are known to improve risk perception and communication. Risk maps are useful when informing policy and management decisions on quarantine, vaccination, and surveillance intended to prevent or control threats to human, animal, or environmental health interface (i.e., One Health). We hypothesized that researcher-stakeholder engagement throughout the research process could enhance the utility of risk maps; while identifying opportunities to improve data collection, analysis, interpretation, and, ultimately, implementation of scientific/evidence-based management and policy measures. Three case studies were conducted to test this process of co-creation of scientific knowledge, using spatiotemporal epidemiological approaches, all related to One Health problems affecting Minnesota. Our interpretation of the opportunities, challenges, and lessons learned from the process are summarized from both researcher and stakeholder perspectives. By sharing our experience we intend to provide an understanding of the expectations, realizations, and “good practices” we learned through this slow-moving iterative process of co-creation of knowledge. We hope this contribution benefits the planning of future transdisciplinary research related to risk map-based management of One Health problems.

Koi herpesvirus and carp oedema virus: Infections and coinfections during mortality events of wild common carp in the United States

Koi herpesvirus (KHV; cyprinid herpesvirus-3) and carp oedema virus (CEV) are important viruses of common and koi carp (Cyprinus carpio); however, the distribution of these viruses in wild common carp in North America is largely unknown. During the summers of 2017 and 2018, 27 mass mortalities of common carp were reported from four states in the USA (Minnesota, Iowa, Pennsylvania and Wisconsin), the majority of which were distributed across eight major watersheds in southern Minnesota. Samples from 22 of these mortality events and from five clinically healthy nearby carp populations were screened for KHV, CEV and SVCV using real-time polymerase chain reaction (qPCR). KHV was confirmed in 13 mortality events, CEV in two mortality events and coinfections of KHV/CEV in four mortality events. Nucleotide sequence analysis revealed that the KHV and CEV detected here are closely related to European lineages of these viruses. While molecular detection alone cannot conclusively link either virus with disease, the cases described here expand the known range of two important viruses. This is also the first reported detection of KHV and CEV coinfections in wild carp populations.

Antibiotic Resistance Genes in Freshwater Trout Farms in a Watershed in Chile

Point sources such as wastewater treatment plants, terrestrial agriculture, and aquaculture may release antibiotic residues, antibiotic resistant bacteria, and antibiotic resistance genes (ARGs) into aquatic ecosystems. However, there is a lack of quantitative studies attributing environmental ARG abundance to specific sources. The goal of this study was to evaluate the role of freshwater trout farms in the release and dissemination of ARGs into the environment. Sediment samples upstream and downstream from five rainbow trout farms were collected over time in southern Chile. A microfluidic quantitative polymerase chain reaction approach was used to quantify an ARG array covering different mechanisms of resistance, and data were analyzed using principal component analysis (PCA) and linear mixed regression models. Surveys were also conducted to obtain information about management practices, including antibiotic use, at the farms. Florfenicol and oxytetracycline were used at these farms, although at different rates. A total of 93 samples were analyzed. In the PCA, , , , , (A), (B), (C), (W), and grouped together. A statistically significant increase in abundance of , , , and several genes was found downstream from the farms compared with upstream sites, and retention ponds had the highest ARG abundance at each site. Antibiotic resistance gene levels returned to baseline at an average distance of 132.7 m downstream from the farms. Although results from this study indicate an influence of trout farms on the presence of ARGs in the immediate environment, the extent of their contribution to ARG dissemination is unknown and deserves further investigation.

Retrospective and Predictive Investigation of Fish Kill Events

Fish kill investigations are critical to understanding threats to aquatic ecosystems and can serve as a measure of environmental disruption as well as an early indicator of emerging disease. The goal of this study was to analyze historical data related to such events among wild fish populations in Minnesota in order to assess the quality and completeness of the data and potential trends in fish kills. After excluding events with incomplete data (e.g., in which the location was not reported), we analyzed 225 unique fish kills from 2003 to 2013 that were recorded in two Minnesota Department of Natural Resources databases. The most reported fish kills occurred during 2007 (n = 41) and during the month of June (n = 81) across all years. Centrarchid species were present in the most fish kills (138), followed by cyprinid and ictalurid species, which were present in 53 and 40 events, respectively. Environmental factors were the most common cause of death reported. Models of environmental factors revealed that the maximum nighttime land surface temperature was the most critical factor in fish mortality, followed by changes in primary productivity and human disturbances. During the course of this study, data gaps were identified, including underreporting, inconsistent investigation, and the lack of definitive diagnoses, making interpretation of our results challenging. Even so, understanding these historical trends and data gaps can be useful in generating hypotheses and advancing data collection systems for investigating future fish kills. Our study is a primer investigation of fish kills providing information on the plausible areas, seasons, and fish groups at risk that can guide active environmental monitoring and epidemiological surveillance of fishes.

A Probability Co-Kriging Model to Account for Reporting Bias and Recognize Areas at High Risk for Zebra Mussels and Eurasian Watermilfoil Invasions in Minnesota

Zebra mussels (ZMs) (Dreissena polymorpha) and Eurasian watermilfoil (EWM) (Myriophyllum spicatum) are aggressive aquatic invasive species posing a conservation burden on Minnesota. Recognizing areas at high risk for invasion is a prerequisite for the implementation of risk-based prevention and mitigation management strategies. The early detection of invasion has been challenging, due in part to the imperfect observation process of invasions including the absence of a surveillance program, reliance on public reporting, and limited resource availability, which results in reporting bias. To predict the areas at high risk for invasions, while accounting for underreporting, we combined network analysis and probability co-kriging to estimate the risk of ZM and EWM invasions. We used network analysis to generate a waterbody-specific variable representing boater traffic, a known high risk activity for human-mediated transportation of invasive species. In addition, co-kriging was used to estimate the probability of species introduction, using waterbody-specific variables. A co-kriging model containing distance to the nearest ZM infested location, boater traffic, and road access was used to recognize the areas at high risk for ZM invasions (AUC = 0.78). The EWM co-kriging model included distance to the nearest EWM infested location, boater traffic, and connectivity to infested waterbodies (AUC = 0.76). Results suggested that, by 2015, nearly 20% of the waterbodies in Minnesota were at high risk of ZM (12.45%) or EWM (12.43%) invasions, whereas only 125/18,411 (0.67%) and 304/18,411 (1.65%) are currently infested, respectively. Prediction methods presented here can support decisions related to solving the problems of imperfect detection, which subsequently improve the early detection of biological invasions.

Genomic characterization of a novel calicivirus, FHMCV-2012, from baitfish in the USA

During regulatory sampling of fathead minnows (Pimephales promelas), a novel calicivirus was isolated from homogenates of kidney and spleen inoculated into bluegill fry (BF-2) cells. Infected cell cultures exhibiting cytopathic effects were screened by PCR-based methods for selected fish viral pathogens. Illumina HiSeq next generation sequencing of the total RNA revealed a novel calicivirus genome that showed limited protein sequence similarity to known homologs in a BLASTp search. The complete genome of this fathead minnow calicivirus (FHMCV) is 6564 nt long, encoding a polyprotein of 2114 aa in length. The complete polyprotein shared only 21% identity with Atlantic salmon calicivirus,followed by 11% to 14% identity with mammalian caliciviruses. A molecular detection assay (RT-PCR) was designed from this sequence for screening of field samples for FHMCV in the future. This virus likely represents a prototype species of a novel genus in the family Caliciviridae, tentatively named "Minovirus".

Novel Methods in Disease Biogeography: A Case Study with Heterosporosis

Disease biogeography is currently a promising field to complement epidemiology, and ecological niche modeling theory and methods are a key component. Therefore, applying the concepts and tools from ecological niche modeling to disease biogeography and epidemiology will provide biologically sound and analytically robust descriptive and predictive analyses of disease distributions. As a case study, we explored the ecologically important fish disease Heterosporosis, a relatively poorly understood disease caused by the intracellular microsporidian parasite Heterosporis sutherlandae. We explored two novel ecological niche modeling methods, the minimum-volume ellipsoid (MVE) and the Marble algorithm, which were used to reconstruct the fundamental and the realized ecological niche of H. sutherlandae, respectively. Additionally, we assessed how the management of occurrence reports can impact the output of the models. Ecological niche models were able to reconstruct a proxy of the fundamental and realized niche for this aquatic parasite, identifying specific areas suitable for Heterosporosis. We found that the conceptual and methodological advances in ecological niche modeling provide accessible tools to update the current practices of spatial epidemiology. However, careful data curation and a detailed understanding of the algorithm employed are critical for a clear definition of the assumptions implicit in the modeling process and to ensure biologically sound forecasts. In this paper, we show how sensitive MVE is to the input data, while Marble algorithm may provide detailed forecasts with a minimum of parameters. We showed that exploring algorithms of different natures such as environmental clusters, climatic envelopes, and logistic regressions (e.g., Marble, MVE, and Maxent) provide different scenarios of potential distribution. Thus, no single algorithm should be used for disease mapping. Instead, different algorithms should be employed for a more informed and complete understanding of the pathogen or parasite in question.

Forecasting distributions of an aquatic invasive species (Nitellopsis obtusa) under future climate scenarios

Starry stonewort (Nitellopsis obtusa) is an alga that has emerged as an aquatic invasive species of concern in the United States. Where established, starry stonewort can interfere with recreational uses of water bodies and potentially have ecological impacts. Incipient invasion of starry stonewort in Minnesota provides an opportunity to predict future expansion in order to target early detection and strategic management. We used ecological niche models to identify suitable areas for starry stonewort in Minnesota based on global occurrence records and present-day and future climate conditions. We assessed sensitivity of forecasts to different parameters, using four emission scenarios (i.e., RCP 2.6, RCP 4.5, RCP 6, and RCP 8.5) from five future climate models (i.e., CCSM, GISS, IPSL, MIROC, and MRI). From our niche model analyses, we found that (i) occurrences from the entire range, instead of occurrences restricted to the invaded range, provide more informed models; (ii) default settings in Maxent did not provide the best model; (iii) the model calibration area and its background samples impact model performance; (iv) model projections to future climate conditions should be restricted to analogous environments; and (v) forecasts in future climate conditions should include different future climate models and model calibration areas to better capture uncertainty in forecasts. Under present climate, the most suitable areas for starry stonewort are predicted to be found in central and southeastern Minnesota. In the future, suitable areas for starry stonewort are predicted to shift in geographic range under some future climate models and to shrink under others, with most permutations indicating a net decrease of the species’ suitable range. Our suitability maps can serve to design short-term plans for surveillance and education, while future climate models suggest a plausible reduction of starry stonewort spread in the long-term if the trends in climate warming remain.

Potential distribution of the viral haemorrhagic septicaemia virus in the Great Lakes region

Viral haemorrhagic septicaemia virus (VHSV) genotype IVb has been responsible for large-scale fish mortality events in the Great Lakes of North America. Anticipating the areas of potential VHSV occurrence is key to designing epidemiological surveillance and disease prevention strategies in the Great Lakes basin. We explored the environmental features that could shape the distribution of VHSV, based on remote sensing and climate data via ecological niche modelling. Variables included temperature measured during the day and night, precipitation, vegetation, bathymetry, solar radiation and topographic wetness. VHSV occurrences were obtained from available reports of virus confirmation in laboratory facilities. We fit a Maxent model using VHSV-IVb reports and environmental variables under different parameterizations to identify the best model to determine potential VHSV occurrence based on environmental suitability. VHSV reports were generated from both passive and active surveillance. VHSV occurrences were most abundant near shore sites. We were, however, able to capture the environmental signature of VHSV based on the environmental variables employed in our model, allowing us to identify patterns of VHSV potential occurrence. Our findings suggest that VHSV is not at an ecological equilibrium and more areas could be affected, including areas not in close geographic proximity to past VHSV reports.

A multiplex RT-PCR assay for the detection of fish picornaviruses

With the emergence of high profile fish diseases in the Great Lakes region, surveillance and regulatory inspections of fish populations have increased. This has resulted in a better understanding of known pathogens and isolation of many new pathogens of fish. In this study, a multiplex RT-PCR assay was developed for the detection of three newly discovered fish picornaviruses: bluegill picornavirus-1 (BGPV-1), fathead minnow picornavirus (FHMPV), and eel picornavirus-1 (EPV-1). This assay was found to be very sensitive with a detection limit of 81.9pg/μl of extracted RNA from a pool of FHMPV and BGPV-1 and was able to detect 501 and 224 gene copies/μl of BGPV-1 and FHMPV, respectively. The assay was highly reproducible and did not cross react with other closely related pathogens. We believe that this new assay provides a rapid and cost effective tool for confirming cell culture isolates and conducting prevalence studies of these newly detected fish picornaviruses.

Description of the Microsporidian Parasite, Heterosporis sutherlandae n. sp., Infecting Fish in the Great Lakes Region, USA

Heterosporosis is an increasingly important microsporidian disease worldwide, impacting wild and farmed raised fishes in both marine and freshwater environments. A previously undescribed species (Heterosporis sp.), with widespread distribution in the Great Lakes region, was the subject of this study. Three angler-caught fish were submitted to the Minnesota Veterinary Diagnostic Laboratory from 2009–2010 with lesions caused by intracellular proliferation of parasitic spores, resulting in destruction and eventual widespread necrosis of the host skeletal muscles. Mature ovoid (5.8 x 3.5μm) spores of a microsporidian parasite, consistent with the genus Heterosporis, were observed by light and electron microscopy. Molecular identification was performed using primer walking to obtain a near-complete rRNA gene sequence (~3,600 bp). A unique species of Heterosporis was identified, demonstrating less than 96% sequence identity to other published Heterosporis sp. on the basis of partial rRNA gene sequence analysis. Heterosporis sutherlandae n. sp. (formerly Heterosporis sp.) was identified in yellow perch (Perca flavescens), northern pike (Esox lucius) and walleye (Sander vitreus) from inland lakes in Minnesota and Wisconsin. Previous research suggests this species may be even more widespread in the Great Lakes region and should be reexamined using molecular techniques to better understand the distribution of this novel species.

Detection and surveillance of viral hemorrhagic septicemia virus using real-time RT-PCR. I. Initial comparison of four protocols

Eight laboratories worked collectively to evaluate 4 real-time RT-PCR (rRT-PCR) protocols targeting viral hemorrhagic septicemia virus (VHSV) being considered for deployment to a USA laboratory testing network. The protocols utilized previously published primers and probe sets developed for detection and surveillance of VHSV. All participating laboratories received and followed a standard operating protocol for extraction and for each of the rRT-PCR assays. Performance measures specifically evaluated included limit of detection (defined as the smallest amount of analyte in which 95% of the samples are classified as positive), analytical specificity, assay efficiency across genotype representatives, within- and between-plate variation within a laboratory, and variation between laboratories using the same platform, between platforms, and between software versions. This evaluation clearly demonstrated that the TaqMan®-based assay developed by Jonstrup et al. (2013; J Fish Dis 36:9-23) produced the most consistent analytical performance characteristics for detecting all genotypes of VHSV across the 8 participating laboratories.

Detection and surveillance of viral hemorrhagic septicemia virus using real-time RT-PCR. II. Diagnostic evaluation of two protocols

Two real-time reverse transcription polymerase chain reaction (rRT-PCR) assays under consideration for deployment to multiple testing laboratories across the USA were evaluated for diagnostic sensitivity and specificity on tissue homogenates obtained from natural and experimental viral hemorrhagic septicemia (VHS)-infected fish. Estimates for diagnostic specificity using virus isolation as the reference method were similar between laboratories regardless of the assay. Diagnostic sensitivity estimates of 0.96 (95% CI: 0.95, 0.97) for Jonstrup et al. (2013)'s assay (J Fish Dis 36:9-23) exceeded the diagnostic sensitivity of 0.85 (95% CI: 0.83, 0.87) for Phelps et al. (2012)'s assay (J Aquat Anim Health 24:238-243). The Jonstrup rRT-PCR assay is robust as demonstrated by high sensitivity and specificity estimates across laboratories and can be used as a valuable tool for targeted surveillance and for testing of suspect VHSV samples.

Isolation and molecular characterization of a novel picornavirus from baitfish in the USA

During both regulatory and routine surveillance sampling of baitfish from the states of Illinois, Minnesota, Montana, and Wisconsin, USA, isolates (n = 20) of a previously unknown picornavirus were obtained from kidney/spleen or entire viscera of fathead minnows (Pimephales promelas) and brassy minnows (Hybognathus hankinsoni). Following the appearance of a diffuse cytopathic effect, examination of cell culture supernatant by negative contrast electron microscopy revealed the presence of small, round virus particles (∼ 30-32 nm), with picornavirus-like morphology. Amplification and sequence analysis of viral RNA identified the agent as a novel member of the Picornaviridae family, tentatively named fathead minnow picornavirus (FHMPV). The full FHMPV genome consisted of 7834 nucleotides. Phylogenetic analysis based on 491 amino acid residues of the 3D gene showed 98.6% to 100% identity among the 20 isolates of FHMPV compared in this study while only 49.5% identity with its nearest neighbor, the bluegill picornavirus (BGPV) isolated from bluegill (Lepomis macrochirus). Based on complete polyprotein analysis, the FHMPV shared 58% (P1), 33% (P2) and 43% (P3) amino acid identities with BGPV and shared less than 40% amino acid identity with all other picornaviruses. Hence, we propose the creation of a new genus (Piscevirus) within the Picornaviridae family. The impact of FHMPV on the health of fish populations is unknown at present.

Comparison of treatments to inactivate viral hemorrhagic septicemia virus (VHSV-IVb) in frozen baitfish

Current US state and federal fish health regulations target the spread of viral hemorrhagic septicemia virus-IVb (VHSV-IVb) through movement restrictions of live fish; however, they largely ignore the potential for the virus to be spread through commercial distribution and use of frozen baitfish from VHSV-IVb-positive regions. Some state laws do require treatment of frozen baitfish to inactivate VHSV, and additional methods have been proposed, but few scientific studies have examined the efficacy of these treatments. In this study, bluegills Lepomis macrochirus were challenged with VHSV-IVb and frozen to represent standard industry methods, disinfected by various treatments, and tested for infectious VHSV-IVb using virus isolation. The virus was isolated from 70% of fish subjected to 3 freeze/thaw cycles. All other treatment methods were effective in inactivating the virus, including treatment with isopropyl alcohol, mineral oil, salt and borax, and dehydration. Dehydration followed by rehydration is rapid and effective, and therefore, seems to be the best option for inactivating VHSV-IVb present in frozen baitfish while maintaining their usefulness as bait.

Spring viremia of carp virus in Minnehaha Creek, Minnesota

Spring viremia of carp virus (SVCV) causes a highly contagious and serious disease of freshwater cyprinid fishes, generating significant economic and ecological impacts throughout the world. The SVCV is therefore listed as a notifiable pathogen by the International Organization for Animal Health. In June 2011, a significant mortality event of wild common carp Cyprinus carpio occurred in Minnehaha Creek near its confluence with Mississippi River Pool 2 in Minneapolis, Minnesota. Clinical signs of moribund fish included hemorrhagic lesions in the skin, eyes, and internal adipose tissue. The SVCV was isolated from pooled kidney and spleen of the fish. Rhabdovirus particles were seen upon examination of infected cell culture fluid by electron microscopy. The virus was confirmed to be SVCV subtype Ia by reverse transcription PCR and sequencing. This is the first report of SVCV within the state of Minnesota and the ninth documented case in North America.

The use of a one-step real-time reverse transcription polymerase chain reaction (rRT-PCR) for the surveillance of viral hemorrhagic septicemia virus (VHSV) in Minnesota

Viral hemorrhagic septicemia virus (VHSV) is a highly contagious and pathogenic virus of fish. The virus infects more than 70 fish species worldwide, in both fresh and salt water. A new viral strain (VHSV-IVb) has proven both virulent and persistent, spreading throughout the Great Lakes of North America and to inland water bodies in the region. To better understand the geographic distribution of the virus, we used a modified real-time reverse transcription polymerase chain reaction (rRT-PCR) assay for high-throughput testing of fish for VHSV. The assay was shown to be twice as sensitive as the gold standard, virus isolation, and did not cross react with other viruses found in fish. In addition, the diagnostic turnaround time was reduced from 28 to 30 d for virus isolation to 2-4 d for rRT-PCR. To demonstrate the usefulness of the rRT-PCR assay, 115 high-priority water bodies in Minnesota were tested by both methods from April 2010 to June 2011. All survey sites tested negative for VHSV by both methods. The survey results have informed fisheries managers on the absence of VHSV in Minnesota and have better prepared them for the eventual arrival of the disease. In addition, the results demonstrate the value of this rRT-PCR as a surveillance tool to rapidly identify an outbreak so that it can be controlled in a timely manner.

Vertical transmission of Ovipleistophora ovariae (microspora) within the eggs of the golden shiner

Fertilized eggs collected from broodfish infected by Ovipleistophora ovariae were tested by quantitative polymerase chain reaction (PCR) and found to be positive for the O. ovariae genome at 7.77 X 10(2) to 3.26 x 10(7) copies per microgram of host DNA. Fry hatched from these eggs contained from 1.37 X 10(2) to 9.89 X 10(6) copies of the O. ovariae genome per microgram of host DNA. Surface treatments of fertilized eggs with 150 mg formalin/L (used by farms as a fungicide) or a 1.5% solution of sodium sulfite (which removes the adhesive egg matrix) did not reduce vertical transmission to fry. Treatment of eggs with a 10% solution of bleach or a proprietary commercial DNA denaturant did not reduce the number of egg-associated copies of the O. ovariae genome. Histology of ovaries of infected fish demonstrated spores within the oocytes. However, no spores were observed by histology in positive fry hatched from infected eggs. The PCR and histological demonstration of the presence of O. ovariae spores in oocytes and fry, and the failure of strong DNA denaturants to reduce egg-associated copies, give evidence that O. ovariae is vertically transmitted within eggs.

Validation of a quantitative PCR diagnostic method for detection of the microsporidian Ovipleistophora ovariae in the cyprinid fish Notemigonus crysoleucas

Microsporidian parasites are easily detected by light microscopy when infections are heavy and spores are present. However, early infections without spores, or light infections with low numbers of spores, are easily missed. This limitation has made it difficult to conduct investigations into microsporidian prevalence and transmission. In this study, we developed a quantitative TaqMan polymerase chain reaction assay to assess the presence of Ovipleistophora ovariae in the tissues of the cyprinid fish Notemigonus crysoleucas (golden shiner). The efficiency of the primer set was 100.8%, with a correlation coefficient of threshold position to copy number of 0.997 over 9 logs using a plasmid containing the cloned reaction product. No product was produced from other closely related microsporidian species (Nucleospora salmonis, Pseudoloma neurophila, Glugea stephani, Heterosporis sp., and O. mirandella). The coefficient of variation for replicate assays done on different days was 12.4%. The assay detects O. ovariae reliably at less than 10 genomic copies and 0.14 spores per reaction, but maximum sensitivity is only achieved when sonication is included as part of the DNA purification step. Using the assay, we found 4.44 x 10(1) to 7.91 x 10(6) copies microg(-1) host DNA in female golden shiners, with the spore density increasing during the spawning season. The parasite was also detected for the first time in the testes of male golden shiners at 2.60 x 10(1) to 8.62 x 102 copies microg(-1) host DNA.