Next Article in Journal
The Metabolic Regulation of Amino Acid Synthesis Counteracts Reactive Nitrogen Stress via Aspergillus nidulans Cross-Pathway Control
Previous Article in Journal
First Record of Flavocillium subprimulinum (Cordycipitaceae, Hypocreales) in Mexico: Morphological and Molecular Characterisation, Nematocidal Activity of Its Liquid Culture Filtrates against Haemonchus contortus and Protease Activity
Previous Article in Special Issue
Establishment and Application of Real-Time Fluorescence Quantitative PCR Detection Technology for Metschnikowia bicuspidata Disease in Eriocheir sinensis
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
JoF
Volume 10
Issue 1
10.3390/jof10010057
jof-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Specific Phylotypes of Saprolegnia parasitica Associated with Atlantic Salmon Freshwater Aquaculture

by
Kypher Varin Shreves
1,
Marcia Saraiva
1,
Tahmina Ruba
1,2,
Claire Miller
3,
E. Marian Scott
3,
Debbie McLaggan
1 and
Pieter van West
1,*
1
International Centre for Aquaculture Research and Development (ICARD), Aberdeen Oomycete Laboratory, Institute of Medical Sciences, University of Aberdeen, Foresterhill, Aberdeen AB25 2ZD, UK
2
Department of Pathology, Faculty of Veterinary Science, Bangladesh Agricultural University, Mymensingh 2202, Bangladesh
3
School Mathematics and Statistics, University of Glasgow, Glasgow G12 8TA, UK
*
Author to whom correspondence should be addressed.
J. Fungi 2024, 10(1), 57; https://doi.org/10.3390/jof10010057
Submission received: 16 October 2023 / Revised: 13 December 2023 / Accepted: 28 December 2023 / Published: 9 January 2024
(This article belongs to the Special Issue Fungal Infections in Fishes and Aquatic Invertebrates, Second Edition)
Browse Figures
Review Reports Versions Notes

Abstract

:
Saprolegniosis is a major destructive disease in freshwater aquaculture. The destructive economic impact of saprolegniosis on freshwater aquaculture necessitates further study on the range of Saprolegnia species within Atlantic salmon fish farms. This study undertook a thorough analysis of a total of 412 oomycete and fungal isolates that were successfully cultured and sequenced from 14 aquaculture sites in Scotland across a two-year sampling period. An ITS phylogenetic analysis of all isolates was performed according to whether they were isolated from fish or water samples and during enzootic or epizootic periods. Several genera of oomycetes were isolated from sampling sites, including Achlya, Leptolegnia, Phytophthora, and Pythium, but by far the most prevalent was Saprolegnia, accounting for 66% of all oomycetes isolated. An analysis of the ITS region of Saprolegnia parasitica showed five distinct phylotypes (S2–S6); S1 was not isolated from any site. Phylotype S2 was the most common and most widely distributed phylotype, being found at 12 of the 14 sampling sites. S2 was overwhelmingly sampled from fish (93.5%) and made up 91.1% of all S. parasitica phylotypes sampled during epizootics, as well as 67.2% of all Saprolegnia. This study indicates that a single phylotype may be responsible for Saprolegnia outbreaks in Atlantic salmon fish farms, and that water sampling and spore counts alone may be insufficient to predict Saprolegnia outbreaks in freshwater aquaculture.

1. Introduction

The total annual production losses in freshwater aquaculture due to saprolegniosis remain consistently high, with losses of 50% having been reported for over two decades, and with Saprolegnia responsible for at least 10% of all annual salmonid economic loss worldwide [1,2,3,4,5,6]. Past publications have identified species within the genus Saprolegnia as the cause of these infections [3,4,5,7,8,9,10]; however, the taxonomic classification of Saprolegnia spp., and oomycetes in general, has often been ambiguous.
The classical identification of Saprolegnia species was based upon the morphological characteristics of sexual structures, such as the oogonia, oospores, and antheridia [11,12,13]. Observations of such structures can often be difficult, particularly as the in vitro growth of many Saprolegnia isolates fails to produce these sexual structures. Following the introduction of molecular identification techniques, it became apparent that these methods could be ambiguous and misleading, resulting in the taxonomical misclassification of species (see [14] for a full list). The phylogeny of Saprolegnia has been re-evaluated and corrected [8,14], resulting in a robust model taxonomy which groups Saprolegnia species into 23 clusters based on genetic sequencing of the internal transcribed spacer (ITS) region of the nuclear ribosomal DNA (nrDNA) [14]. One finding of particular note was the re-classification of S. hypogyna, S. parasitica, and S. salmonis, (three species previously thought to be associated with saprolegniosis), into a single species, S. parasitica [8,14].
Saprolegniosis is caused by several species within the genus Saprolegnia, in particular, the species S. australis, S. diclina, and S. parasitica [3,5,7,8,9,10,14]; however, only a few published studies have investigated the range of Saprolegnia species within Atlantic salmon fish farms. Across farms in Chile and Canada, differences were reported in the species diversity of Saprolegnia, with S. parasitica making up the overall majority of isolated species [4,15]. Whether this may also be true of Scottish sites has not been determined. For Saprolegnia host–pathogen specificity, species-level identification is typically considered; however, it has been suggested that the interaction may actually be strain-specific [15]. Furthermore, whilst several studies have investigated Saprolegnia pathogenicity through the inference of molecular markers or through infection trials, this is the first study that has determined the range of Saprolegnia species present in aquaculture during epizootic and non-epizootic periods across several Scottish Atlantic salmon aquaculture sites. The distribution preference of species and strains between host and environment was also established; comparisons were also made between species isolated from fish epidermal tissue and environmental tank water. This work provides insight into the distribution profiles of Saprolegnia species that are present in Scottish freshwater aquaculture. The variation in Saprolegnia species abundance between water and fish suggests adaptive strategies to increase the probability of attachment to target hosts.

2. Materials and Methods

2.1. Sample Collection and Culturing

Samples from tank water and the epidermis of fish were taken on a monthly basis across fourteen Atlantic salmon freshwater aquaculture sites (anonymized as A to N). Sampling was performed by on-site personnel trained in proper sampling techniques and provided with all necessary sampling equipment.
Water samples were collected, whereby 1 mL of water from the sampling tank was placed into each well of a sterile 24 well plate. A single sterilized hemp seed was introduced to each well, acting as a growth medium for any oomycetes present in the water sample. The plates were sealed and sent by courier to the Aberdeen laboratory, where they were incubated at 12 °C. The plates were checked regularly over a 2-week period; any hemp seeds displaying signs of growth were transferred onto selective agar plates: Potato Dextrose Agar (60 mm plates) supplemented with antibiotics (0.1% (w/v) vancomycin (Alfa Aesar, Haverhill, MA, USA) and 0.5% (w/v) ampicillin (Sigma-Aldrich, St. Louis, MO, USA) and the anti-fungal compound 0.02% (w/v) pimaricin (Sigma-Alrich) (PDA-VAP).
To sample the epidermis of fish, two biopsies were taken from 9 fish per monthly sampling period. Fish from the same tank were selected randomly. The first biopsy was taken from behind the left pectoral fin; a 1 cm2 skin sample was transferred directly to 6-well PDA-VAP plates, skin down. The second biopsy comprised the left pectoral fin, transferred directly to 6-well PDA-VAP plates with the fin fanned out. The plates were sealed and sent by courier to the Aberdeen laboratory, where they were incubated at 12 °C. Purification was achieved through repeated sub-culturing by transferring agar plugs to fresh PDA-VAP plates. The purification of the cultures was confirmed by culturing in Pea Broth [16] with no antimicrobials at 24 °C. A clear broth after 24 h confirmed no bacterial contamination.

2.2. DNA Extraction

Pure cultures were grown in Pea Broth for ca. 48 h. DNA extraction was performed on collected, washed mycelia in accordance with the phenol/chloroform protocol described in [17]. Essentially, the mycelia were transferred to a 2 mL screw-cap tube with sterile glass beads. A volume of 800 µL of an extraction buffer (10 mM tris-HCI (pH 8), 50 mM EDTA, 0.5% (w/v) SDS, 0.5% (v/v) Triton X-100, and 0.5% (v/v) Tween 20) was added to the tube, with 2 µL of proteinase K (20 mg/mL, Promega, Madison, WI, USA). The samples were snap-frozen in liquid nitrogen, placed into a Fastprep24 5G (MP BiomedicalsTM, Santa Ana, CA, USA), and run for four cycles of 45 s at 6.5 m/s to break open the cells by mechanical force. The samples were incubated at 55 °C for 30 min, with the tubes inverted every 10 min. The samples were cooled to room temperature (20 °C); then, 800 µL of equilibrated phenol/chloroform/isoamyl alcohol (25:24:1) (ACROS OrganicsTM, Geel, Belgium) was added. The tubes were vortexed and then centrifuged at 10,000× g for 10 min to separate the phases. Up to 1 mL of the upper phase was transferred to another 2 mL tube containing 1 mL of molecular-grade isopropanol (Fisher Scientific, Hampton, NH, USA) and incubated at −20 °C for at least 15 min. The samples were then centrifuged at 4 °C at 10,000× g for 10 min to precipitate genomic DNA. The supernatant was removed and the DNA pellet was washed with 70% (v/v) ethanol (molecular grade, Sigma Aldrich), centrifuged (1 min at 10,000× g), then washed with 100% molecular-grade ethanol, and air dried. The concentration and purity of the genomic DNA were assessed through Nanodrop and gel electrophoresis.

2.3. Species Identification

Taxonomic identification was performed by targeting the internal transcribed spacer (ITS) region using the primer set ITS_4_ALT (Forward: 5′-TCCTCCGCTTATTGATATG-3′) and ITS_5_ALT (Reverse: 5′-TGAAAAGTCGTAACAAGGTT-3′). The primers were selected for their universal binding, known to amplify a wide range of oomycetes and fungi [18,19]), as well as their large amplicon coverage, covering the ITS1, 5.8S rDNA, and ITS2 regions. The PCR reaction mix (25 µL per sample) contained a 1 µL DNA template (5 ng), 5 µL of 5X Phusion® High Fidelity Buffer, 5 µL of MgCL2 (25 mM), 2 µL of dNTP’s (2.5 mM), 1 µL of ITS_5_ALT (10 µM), 1 µL of ITS_4_ALT (10 µM), 0.25 µL of Phusion® High-Fidelity DNA Polymerase (0.5 U) (Thermo Scientific, Waltham, MA, USA), and 4.75 µL ddH2O. The PCR thermal cycle consisted of an initial denaturation of 30 s at 95 °C, followed by 35 cycles of 35 s at 95 °C, 15 s at 59 °C, 75 s at 72 °C, and a final step of 7 min at 72 °C. The PCR products were sequenced by Sanger sequencing (Eurofins, Luxembourg). The sequences were assessed and edited using Geneious Prime v2020.1.2 The chromatogram of each sequence was assessed for the quality of the reads, with any sequences that appeared to be contaminated or of low quality discarded from the dataset. In cases where both paired-end reads of a single sample were considered to be of good quality, both reads were aligned into a single consensus sequence using de novo assembly with the Geneious assembler at the highest sensitivity. In the event that assembly failed, the read with the highest quality was kept as the representative read for the sample. An initial trimming of both de novo assembled contigs and single reads was performed automatically, with an error probability limit of 0.01. The species identification of each sample was based on the results of the NCBI Basic Local Alignment Search Tool (BLAST) [20]. In the cases where isolates were assigned species-level assignment based on their blast results, the isolates were grouped based on their NCBI blast results using the most frequent nomenclature of the top 10–20 results. To ensure that the isolate sequences were then correctly grouped, an independent alignment was performed within each genus-level group. If it was found that an isolate assigned to a species did not closely match other sequences within the same assigned species group, then the sequence was re-assessed. Re-assessment typically resulted in the group being grouped into a species which had a much more similar sequence or assigned to an ‘sp’ group if no similar sequences could be found. A reference dataset attained from [14,15] was used for the resolution of the S. parasitica cluster, with a single reference sequence selected for each S. parasitica phylotype within the dataset (Appendix A Table A1). All alignments were performed within Geneious Prime v2020.1.2. using MUSCLE v2.5 [21] with the default settings and a maximum of 8 iterations. Consensus sequences were produced for all isolates with identical ITS sequences and used as a representative sequence for the respective species/phylotype. Alignment trimming was performed in Geneious Prime v2020.1.2.

2.4. Phylogenetic Analysis

The isolates were binned into water and fish based upon their sampling origin. The isolates (water and fish together) were also binned into epizootic and enzootic based on the output of a modified z-score model (Z3), which identified epizootic events from saprolegniosis mortalities (Appendix A, Figure A1). Phylogenetic trees were produced using TOPALi v.2.5 [22]. Model selection for each tree was performed using the model selection tool for Phy maximum likelihood (ML) and MrBayes v3.0. The optimal model was selected based on the resulting values of the Bayesian information criterion (BIC), with the lowest-value model selected. Phylogenetic relationships were assessed using both ML analysis with PhyML [23] and Bayesian inference with MrBayes [24]. Whilst model selection varied between trees, all ML models ran with 1000 bootstraps, and all MrBayes models were run with 5 × 106 generations and a 35% burn-in. Treegraph2 (v2.15.0) [25] was used to visualize the ML trees and integrate support values from the Bayesian inference.

3. Results

3.1. Phylogenetic Analysis of All Isolates Sampled across Scottish Atlantic Salmon Aquaculture Sites

A total of 412 isolates were successfully cultured and sequenced from 14 aquaculture sites across the two-year sampling period (Supplementary Materials Table S1, for complete isolate list); of those, 336 isolates were identified as oomycetes and 76 isolates were identified as true fungi. The 336 isolates were identified as belonging to the oomycete genera Saprolegnia (n = 228, 69.9%), Pythium (n = 95, 23%), Achlya (n = 10, 2.4%), Phytophthora (n = 2, 0.5%), and Leptolegnia (n = 1, 0.2%) (Figure 1). Saprolegnia was the most frequently isolated oomycete genus, followed by Pythium. The 228 Saprolegnia samples isolated consisted exclusively of five Saprolegnia species: S. parasitica (n = 151, 36.7%), S. ferax (n = 33, 8%), S. diclina (n = 20, 4.9%), S. australis (n = 14, 3.4%), and S. delica (n = 10, 2.4%), with S. parasitica being the most frequently sampled species overall. Of the 95 Pythium isolates sampled, P. coloratum was the most abundant (n = 24, 5.8%), followed by Pythium sp.5 (n = 21, 5.1%), and P. flevoense (n = 19, 4.6%). The oomycetes formed two distinct clusters, with Achlya, Leptolegnia, and Saprolegnia grouping in one, and Phytophthora and Pythium grouping in the other (Figure 1).
There was a clear differentiation between the fungal and oomycete isolates, although Mortierella formed a cluster distinct from other fungal samples that was marginally closer to the oomycete cluster (Figure 1). Despite the use of selective media containing pimaricin as an anti-fungal compound, 19% of the total isolates sampled were identified as true fungi (see Section 4).

3.2. Phylogenetic Analysis of Saprolegnia Genus

Isolates of the Saprolegnia genus were clearly grouped with their respective cluster, as designated by Sandoval-Sierra et al. [14] (Figure 2). All Saprolegnia isolates were grouped in a distinct phylogenetic branch. This branch consisted of S. delica, S. australis, Saprolegnia sp.1, S. diclina, S. ferax, S. brachydanionis, and S. parasitica (Figure 2), and almost all members of the branch were isolated from aquaculture sites in this study, with only two members, Saprolegnia sp.1 and S. brachydanionis, not isolated. S. delica and S. ferax displayed no variation in the ITS region. The site dispersal of S. delica and S. ferax was very narrow (Table 1), with 70% and 95% isolated from two different sites, respectively, both recirculating. Isolates of S. australis displayed enough variation in the ITS region to separate them into three different phylotypes, although the majority of S. australis isolates belonged to a single phylotype. The two more abundant phylotypes for S. australis were grouped with the reference sequences (Figure 2); the third phylotype had one isolate but no reference sequence, indicating that this phylotype is either rare, or the reference samples were simply not numerous enough to provide full phylotype diversity. S. australis was found at only two of the fourteen sites (Table 1), both recirculating systems. Isolates of S. diclina were clearly divided into two separate groups based on variation in the ITS, with the number of isolates in each group being twelve and eight (Figure 2). The reference sequences only represented one of these S. diclina groups; however, the NCBI BLAST results confirmed that the unrepresented group was not novel. S. diclina was sampled across just over half of the sites (eight of fourteen), with the total number of isolates sampled being relatively low (20) (Table 1).
The only species of Saprolegnia isolated across all sampling sites was S. parasitica (Table 1). The number of S. parasitica isolates varied greatly between sites, with site C sampling the highest number of all S. parasitica isolates (n = 51, 33.8%).

3.3. Phylogenetic Analysis of Saprolegnia Parasitica

This current analysis showed that the distribution of all 151 S. parasitica isolates were divided into five different phylotypes based on variations in the ITS (Figure 2). For clarity and consistency with the literature, the nomenclature used for the phylotypes in this study is adopted from the study of Sarowar et al. [15], who identified four S. parasitica phylotypes/strains using ITS and Cox1 regions and designated them as S1–S4. S1 was not isolated from any Scottish fish farm; the five isolated phylotypes were designated as S2–S6 (Table 2).
The resolution of the S. parasitica phenotypes isolated was achieved using eight reference sequences (four from [14] and four from [15]). The Scottish isolates were grouped with the five phylotypes as designated [14,15] (Appendix A, Figure A2).
S. parasitica phylotype S2 was by far the most abundant phylotype sampled (70.9%) on Scottish fish farms, and it was also the most widely distributed S. parasitica phylotype, sampled across all but two sites. The second most abundant S. parasitica phylotype was S6 (23.8%) (Table 2). S6 was isolated from only three sites, one recirculating site and two sites geographically isolated from the mainland UK.

3.4. Analysis of Nucleotide Changes of S. parasitica Phylotypes

The ITS consensus sequences from each isolated S. parasitica phylotype (S2–S6), and also S1 from Sarowar et al. [15], were compared with the reference sequence KF717839.1 (SAP0208) [14]. The nucleotides are numbered in accordance with this reference sequence. Differences are seen at 10 different bp positions (Figure 3).

3.5. Phylogenetic Comparison of Saprolegnia between Water and Fish Isolates

When Saprolegnia isolates were grouped into those sampled from the epidermis of Atlantic salmon and those isolated from tank water, more than three times the number of isolates were sampled from fish (n = 172, 75.4%) than from water (n = 56, 25.6%) overall. The isolates in fish and in water were S. australis, S. diclina, and S. delica (Table 3), whereas the number of S. ferax isolated from fish (n = 23) was more than double the number found in water (n = 10). The largest difference in the abundance of water and fish isolates was observed with S. parasitica, with the vast majority of S. parasitica isolates sampled from fish (n = 130, 86.1%). S2, the most abundant S. parasitica phylotype, was overwhelmingly sampled from fish (n = 100), with only seven samples isolated from water (Table 3). The second most abundant phylotype, S6, had 23 isolates sampled from fish and 13 isolates sampled from water.

3.6. Phylogenetic Comparison of Saprolegnia between Enzootic and Epizootic Periods

At the time of sampling, 225 of a total of 228 Saprolegnia isolates could be divided based upon their epizootic status. It was found that the majority (164 isolates (78.5%)) were sampled during enzootic periods, rather than epizootic periods. The only S. parasitica phylotype sampled during both enzootic and epizootic periods was S. parasitica phylotype S2. Phylotypes S4 and S6 were found exclusively during enzootic periods, and S3 and S5 were found exclusively during epizootic periods (Table 4). During epizootic periods, almost all S. parasitica isolates were S2 (91.1%), with the remainder being S5 (6.7%) and S3 (2.2%). Both phylotypes of S. diclina were isolated during enzootic (n = 16) and epizootic (n = 4) periods (Table 4). S. delica and S. australis were also found during both enzootic and epizootic periods; the numbers were low and relatively similar for both periods, whereas two of the S. australis phylotypes were only found during epizootic periods. In contrast, S. ferax was exclusively found during enzootic periods (Table 4).
This was also the case for S. parasitica S6; it was the second most frequently sampled S. parasitica phylotype, and was isolated exclusively during enzootic periods, as was S4, although the numbers for S4 were low (Table 4).

4. Discussion

This study has demonstrated that Saprolegnia is ubiquitous across Scottish freshwater Atlantic salmon aquaculture, with the majority of Saprolegnia isolated during enzootic periods. This study also showed that, at the species level, S. parasitica is the primary Saprolegnia species sampled during saprolegniosis epizootics.
Phylogenetically, all Saprolegnia isolates were grouped within the species-level clusters as described by Sandoval-Sierra et al. [14]. In Chile, Sandoval-Sierra et al. [4] reported the isolation of all Saprolegnia species identified in this study, in addition to Saprolegnia sp.1 and Saprolegnia sp.2. In contrast, in Nova Scotia, no isolates of S. australis, S. delica, or S. diclina were identified, although the addition of S. torulosa was reported [15]. S. parasitica and S. ferax were the only species consistently found within Atlantic salmon aquaculture across all three studies. Although there are some commonalities, it seems that Saprolegnia species richness is not uniform within Atlantic salmon aquaculture at a global or even national level, but instead may be unique to individual aquaculture sites.
In the current study, the majority of isolates of S. australis, S. delica, and S. ferax were found at two sites, and these sites were both recirculating sites. These sites recycle most of their water and have greater biocontrol in comparison to flow-through sites. Recirculating sites are designed to reduce the introduction of pathogenic microorganisms, with, for example, controlled decontamination entrances and extensive water treatment facilities. However, it is possible that such infrastructure also prevents microorganisms from being removed. If treatments inhibit, rather than kill, Saprolegnia, zoospores may encyst and reside in biofilms [26]. Environmental conditions such as temperature are typically more stable within recirculating sites than flow-through sites; this stability may be more favorable for some species of Saprolegnia. The extensive water treatment found within recirculating sites may also reduce the ability of other microorganisms to survive which compete with or inhibit the growth of specific Saprolegnia species, which may partially explain why species such as S. australis, S. delica, and S. ferax are almost exclusively found in recirculating sites and in such abundance.
S. parasitica was by far the most abundant Saprolegnia species isolated, and the only species to be found ubiquitously across 12 out of 14 Scottish aquaculture sites tested. In Atlantic salmon studies, an increased abundance of S. parasitica, relative to other species of Saprolegnia, is a common phenomenon seen across different geographic locations, in both wild and captive populations [4,15,27,28]. This pattern may suggest a strong host–pathogen inter-relationship, with S. parasitica adapting to host environmental conditions and proliferating at a much greater rate than other Saprolegnia species. Studies by Matthews [28] found evidence of high genetic variation among S. parasitica isolates, indicating their adaptive potential.
The phylogenetic analysis of all S. parasitica species isolated in this study showed five distinct phylotypes, more phylotypes than any other species of Saprolegnia isolated. No isolates were obtained of phylotype S1. Phylotype S2 was the most common and widely distributed phylotype, being found at 12 of the 14 sampling sites. Isolates with ITS regions matching the S. parasitica S2 ITS phylotype seem to be the most widely reported S. parasitica phylotype worldwide. S. parasitica isolates that can be designated as S2 are found in Canada, Chile, Croatia, Egypt, Norway, Poland, Scotland, Spain, and Switzerland [4,14,15,29,30,31,32] (Table 5). In a study of oomycete isolates from trout farms in Croatia, 86% of all S. parasitica isolates could be designated as S2, with a single isolate of S6 sampled from water [29]. During an outbreak of saprolegniosis in Nile tilapia in Egypt, Saprolegnia spp. were isolated [30]. Two out of six isolates of S. parasitica recovered from naturally infected O. niloticus could be designated as the S2 phylotype based on ITS sequencing (Table 5) [4,15]. In the study by Ravasi et al. [32], isolates were collected from different locations in Switzerland, including fish hatcheries, fish farms, rivers, and lakes. The isolates were categorized using a multilocus sequence typing (MLST) scheme, using seven housekeeping genes; a diploid sequence type (DST) was assigned to each unique combination of alleles. This allowed a higher resolution of phylotyping, and there were a total of 10 different DSTs, with DST3 being the most common [32]. Based on ITS sequencing, we could assign eight of the ten DSTs to the S2 phylotype. The other two DSTs could be assigned to the S1 phylotype. Thus, from studies reported to date, it seems that a single phylotype dominates other phylotypes in salmonid aquaculture. The S2 ITS sequence displayed a 100% match with the S. parasitica strain CBS 223.65 (C65), a strain originally sampled from Northern pike (Esox lucius) in 1965, and also S. parasitica N12 (VI02736), originally sampled from the parr of Atlantic salmon in Scotland in 2002 [33,34], both of which are widely used as model strains for S. parasitica studies.
The S6 ITS phylotype was the second most frequently isolated phylotype of S. parasitica found at three of the fourteen sampling sites. The highest frequency of S6 was found at a recirculating site; this site had relatively high numbers of both S. ferax and S. parasitica, with S. parasitica phylotypes S6 and S2 being relatively equal in numbers. The other two sites that S6 was isolated from are unique in this study, as they are separated from mainland Scotland, and only S. parasitica S6 was isolated from these sites. The separation of land masses often creates geographic isolation in terrestrial species, as the water between the land masses acts as a barrier for migration. Although reports of the S6 phylotype are much less common than those of S2, isolates have been found in Chile, Ecuador, Poland, Russia, Spain, Scotland [4,8,14,31,33], and Korea [35], and a single isolate was found in Croatia [29] (Table 5). A recent report of the diversity and distribution of culturable Saprolegniaceae species in freshwater ecosystems in Korea obtained a total of nine Saprolegnia strains [35]. There were only four isolates of Saprolegnia parasitica isolates found, all of which could be designated as the S6 phylotype based on the ITS sequence analysis.
S. parasitica phylotypes S3, S4, and S5 were isolated far less frequently than S2 or S6, suggesting that they may be less well adapted to an aquaculture environment. The S3 phylotype was found at two sites and has only been previously reported in Canada [15], though an unpublished submission on GenBank has also identified the phylotype in Scotland (accession number: MW356889.1). The S4 phylotype was found only at one site. From the isolates reported in GenBank, S4 has been found in Argentina, Canada, and the USA [14,15,36]. The S5 phylotype was isolated from only two of the fourteen sampling sites, both flow-through sites. From the isolates reported on GenBank, the S5 phylotype is found in the Czech Republic, Ecuador, Spain, and the USA [14,36,37].
Table 5. Global distribution of S. parasitica phylotypes S1–S6 designated by 600 bp ITS region. The GenBank comparison of ITS sequences illustrates the different countries in which phylotypes S1–S6 are found. The relative numbers of isolates are not depicted but discussed in text.
Table 5. Global distribution of S. parasitica phylotypes S1–S6 designated by 600 bp ITS region. The GenBank comparison of ITS sequences illustrates the different countries in which phylotypes S1–S6 are found. The relative numbers of isolates are not depicted but discussed in text.
S. parasitica Isolate Country of OriginHost/HabitatITS PhylotypeGenBankRef.
SAP1091ArgentinaRiver waterS4KF717864[14]
Isolate #5CanadaSalmo salar 1S1MK849947[15]
Isolate #4CanadaSalmo salar 2S2MK849946[15]
Isolate #49CanadaSalvelinus fontinalisS3MK849963[15]
Isolate #1CanadaWaterS4MK849943[15]
Li16Czech RepublicAstacus astacusS5KF386710[37]
SAP0522ChileSalmo salar 3S2KF717845[4]
SAP0530ChileSalmo salarS6KM095949[4]
Isolate B1L1CroatiaOncorhynchus mykiss 4S2MT555893[29]
Isolate 122CroatiaWater S6MT555889[29]
SAP1230EcuadorRiver waterS5KF717872[14]
SAP1381EcuadorRiver waterS6KF717876[14]
SA221013EgyptOreochromis niloticusS2ON797303[30]
W9KoreaWaterS6ON075413[35]
CBS 223.65NetherlandsEsox luciusS2KF717879[33]
SAP1484NorwaySalmo salarS2KF717880[14]
SAP0254PolandWaterS2KF717840[14]
SAP0257PolandWaterS6KF717842[14]
CBS 113187RussiaN/AS6HQ644005[33]
DD.37.04ScotlandSalmo salarS2OQ678594This study
VI0 5977ScotlandSalmo salar eggsS2HG329736[31]
N12 (VI-02736)ScotlandSalmo salarS2NCBI: txid983306[34]
AA.35.02.04ScotlandAquaculture tank waterS3OQ678418This study
DD.57.10ScotlandSalmo salarS4OQ678598This study
GG.48.02ScotlandSalmo salarS5OQ678653This study
EE.19.06ScotlandSalmo salarS6OQ678612This study
VI0 6009ScotlandSalmo salar eggsS6HG329739[31]
SAP0601SpainWaterS2KF717849[14]
SAP1203SpainWaterS5KF717870[14]
SAP26SpainWaterS6AM228725[8]
S026SwitzerlandSalmo trutta 5S1MH030519[32]
S001SwitzerlandSalmo trutta 6S2MH030499[32]
UNCW314USAN/AS4DQ353545[36]
UNCW373USAN/AS5DQ393557[36]
1 S1 also isolated from brown trout (Salmo trutta), brook trout (Salvelinus fontinalis), Arctic charr (Salvelinus alpinus), and striped bass (Salvelinus alpinus). 2 S2 also isolated from brook trout (S. fontinalis), striped bass (S. alpinus), and eggs of Atlantic salmon (Salmo salar) and brook trout (S. fontinalis). 3 S2 also isolated from rainbow trout (Oncorhynchus mykiss) and from water. 4 S2 also isolated from brown trout (S. trutta), and rainbow trout (O. mykiss) eggs. 5 S2 also isolated from common minnow (Phoxinus phoxinus), grayling (Thymallus thymallus), whitefish (Coregonus clupeaformis), and Atlantic salmon (Salmo salar). 6 S2 also isolated from marble trout (Salmo marmoratus).
Regarding the distribution of isolates between water and fish, 75% of all Saprolegnia isolates were sampled from the epidermis of Atlantic salmon, with the remaining from tank water. The numbers of S. delica, S. diclina, and S. australis, isolates were relatively even between water and fish, whereas a higher number of S. ferax was isolated from fish (70%). The greatest difference between water and fish was seen for the S. parasitica isolates. Most notable was the distribution of the S. parasitica S2 phylotype, with 93.5% of the 107 isolates being sampled from fish; thus, during epizootics, virtually no S. parasitica S2 was detected in water samples, indicating that assessing the risk of saprolegniosis by sampling water for Saprolegnia or Saprolegnia zoospores would not provide a true reflection of saprolegniosis risk.
When an analysis was carried out on the distribution of isolates during enzootic and epizootic periods, all species of Saprolegnia isolated in this study were found during both enzootic and epizootic periods, with the exception of S. ferax, which was found only during enzootic periods. Although S. ferax is a known pathogenic species, it is primarily pathogenic to amphibians [38]. S. ferax has been isolated from infected Atlantic salmon at aquaculture sites in both Chile and Canada [4,15]; however, cases were rare (5.6%, and 7.1% of isolates, respectively), with both studies concluding that S. ferax was not a primary infector. The apparent absence of S. ferax during epizootics was not anticipated, as it was the second most frequently sampled Saprolegnia species and was isolated more often from fish than from water. However, a caveat to this observation is that all S. ferax samples were isolated almost exclusively from a single site (97%), which was a recirculating site. S. australis and S. delica were low in number and relatively evenly distributed between enzootic and epizootic periods. S. diclina was primarily found during enzootic periods (80%). The majority of Saprolegnia isolates sampled during epizootics were identified as S. parasitica (73.8%), which was anticipated, as S. parasitica is known as the primary cause of saprolegniosis in salmon [5]. S6 was never isolated during an epizootic period, supporting suggestions that pathogenicity varies between isolates of S. parasitica [39,40]. It is important to note that, whilst the highest number of isolates of S. ferax and S. parasitica phylotypes S2 and S6 were sampled from site C than any other site, the isolates found across the five epizootics that occurred at site C during the sampling period were exclusively S2. Overall, the S2 phylotype made up 91.1% of all S. parasitica phylotypes sampled during epizootics in Scotland, 97.3% of all S. parasitica strains sampled from infected Atlantic salmon in Chile [4], 63.2% of all S. parasitica strains sampled from infected Atlantic salmon in Canada [15], and 86% of all S. parasitica strains sampled from infected trout in trout farms in Croatia [29].
In addition to Saprolegnia, several other genera of oomycetes were isolated from the sampling sites, including Achlya, Leptolegnia, Phytophthora, and Pythium. Both Leptolegnia and Phytophthora are known to contain highly pathogenic species [6,41,42]; however, neither are associated with Atlantic salmon, and the infrequence with which these isolates were obtained in this study indicated their low abundance. A greater species richness of Achlya was present across farm sites, relative to Leptolegnia and Phytophthora. Whilst infection of Achlya species has been reported in several fish species [6,43,44], there are no reports of Achlya infection in Atlantic salmon. Pythium was the second most frequently isolated genera of oomycete, after Saprolegnia. Of the three most abundantly isolated Pythium species, only P. coloratum and P. flevoense could be identified with confidence via NCBI BLAST. Whilst many species of Pythium are pathogenic, the majority of pathogenic species are plant pathogens, with few species known to infect animals [6,45,46,47,48]. P. coloratum is associated with the infection of carrots [49], and P. flevoense has been associated with infection in Ayu fish larvae (Plecoglossus altivelis) [50]. Pythium species have been isolated from fish; however, very few studies have associated Pythium with infections [45,47,50,51], and no study to date has associated Pythium with infection in Atlantic salmon. Pythium are ubiquitous, found in freshwater worldwide [6,48]. Considering the lack of association with fish infection, the abundance and richness of Pythium species found in aquaculture sites may be influenced by their respective water sources, and it is likely that they have no notable effect on fish. However, an indirect effect of Pythium on the manifestation of saprolegniosis cannot be ruled out.
We also isolated 76 fungal isolates contributing to 19% of the total isolates obtained during this study. The true fungi primarily belonged to the genus Mortierella (n = 64, 15.5%), with Mortierella hyalina being the most frequently isolated fungal species (n = 24, 5.8%). Despite the use of selective media containing pimaricin, an anti-fungal compound, the most common fungal isolates belonged to the genus Mortierella, which is consistent with the finding of past studies that have shown reduced inhibitory effects of pimaricin on the growth of Mortierella [52,53,54]. Furthermore, some Mortierella produce white, cotton-like colonies, which are very similar in appearance to oomycetes such as Pythium and Saprolegnia [55,56,57]. It is likely that the anti-fungal resistance and physical appearance of Mortierella caused a bias during sample selection, resulting in an overabundance of these isolates. For information, we chose to list all isolates found under the sampling regime used; however, no further interpretation regarding true fungi abundance can be made.

5. Conclusions

This study established the diversity of Saprolegnia species within Scottish freshwater Atlantic salmon aquaculture. The variation in Saprolegnia species abundance between water and fish suggest adaptive strategies to increase the probability of attachment to target hosts. This study reports that of the current six different S. parasitica phylotypes in the literature, phylotype S1 was not detected in Scottish freshwater aquaculture; phylotypes S3, S4, and S5 were all very low in abundance. Phylotype S6 was the second most abundant, but by far the most ubiquitous across the farm sites was the S. parasitica S2 phylotype, present in high abundance, particularly on fish, suggesting that this particular phylotype is the primary cause of Atlantic salmon saprolegniosis, with similar observations found in both Canada and Chile.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/jof10010057/s1: Table S1: Details of all 412 isolates from Scottish Atlantic salmon aquaculture sites. Figure S1: Phylogenetic summary of all Saprolegnia species divided between water and fish samples. Figure S2: Phylogenetic summary of all Saprolegnia species divided between enzootic and epizootic sampling periods.

Author Contributions

Conceptualization, P.v.W. and M.S.; methodology, K.V.S., M.S., E.M.S., C.M. and T.R.; software, K.V.S., E.M.S. and C.M.; validation, K.V.S., D.M. and P.v.W.; formal analysis, K.V.S., E.M.S. and C.M.; investigation, K.V.S.; resources, P.v.W.; data curation, K.V.S.; writing—original draft preparation, K.V.S.; writing—review and editing, D.M.; visualization, K.V.S. and D.M.; supervision, P.v.W.; project administration, D.M., M.S. and P.v.W.; funding acquisition, P.v.W. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the BBSRC LINK Research Grant BB/P020224/1 Project: Risk factors for escalating saprolegniosis outbreaks in salmon farms (RIFE-SOS), and The APC was funded by BB/P020224/1. Financial support for Tahmina Ruba was provided by the Prime minister fellowship, Government of the Peoples Republic of Bangladesh.

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The nucleotide sequence of the ITS region for all 412 isolates from the Scottish Atlantic salmon aquaculture sites are deposited in GenBank (accession numbers OQ6784-08 to OQ6787-89).

Acknowledgments

The authors would like to acknowledge the dedicated work of all aquaculture personnel involved in this study who took the time to provide regular sampling.

Conflicts of Interest

The authors declare no conflicts of interest.

Appendix A

Table A1. Reference sequences used in Saprolegnia phylogenetic analysis. Sequences cover the ITS region and are taken from Sandoval-Sierra et al. [14]. Each sample is listed with its respective GenBank accession number, species, phylogenetic cluster (as designated by [14]), and geographic origin.
Table A1. Reference sequences used in Saprolegnia phylogenetic analysis. Sequences cover the ITS region and are taken from Sandoval-Sierra et al. [14]. Each sample is listed with its respective GenBank accession number, species, phylogenetic cluster (as designated by [14]), and geographic origin.
Accession N°Species *Cluster *Origin
AB219386S. diclina1Japan
AM228844S. diclina1United Kingdom
JX087968S. diclina1Canada
AM228782S. sp.12Spain
FN186016S. sp.12Norway
EF460349S. sp.12Germany
AM228725S. parasitica3Spain
FN186030S. parasitica3United Kingdom
FN186045S. parasitica3Chile
DQ393530S. ferax4United Kingdom
HQ644001S. ferax4Netherlands
JX087989S. ferax4Canada
DQ393560S. australis5United Kingdom
AM228829S. australis5Chile
DQ393546S. australis5USA
AM228813S. delica6Spain
AM228816S. delica6Spain
AM228819S. delica6Spain
AM228824S. litoralis7United Kingdom
DQ393558S. litoralis7United Kingdom
HQ643991S. litoralis7United Kingdom
FN186033S. sp.28United Kingdom
EU240098S. sp.28Poland
GU014276S. sp.28USA
AB219378S. subterranea10Japan
DQ393534S. subterranea10USA
HQ644009S. subterranea10Ukraine
AB219379S. torulosa11Japan
DQ393567S. torulosa11Japan
JQ974999S. torulosa11USA
DQ393559S. monilifera12United Kingdom
DQ393562S. monilifera12USA
HQ643996S. monilifera12United Kingdom
AB219396S. terrestris13Japan
AM228851S. terrestris13Norway
DQ393523S. terrestris13USA
DQ393556S. eccentrica14United Kingdom
DQ393563S. eccentrica14Japan
HQ643981S. eccentrica14United Kingdom
AB219388S. furcata16United Kingdom
DQ393561S. furcata16United Kingdom
EU544191Saprolegnia sp. 517USA
EU124761Saprolegnia sp. 618USA
JQ974993Saprolegnia sp. 618USA
JQ974994Saprolegnia sp. 618USA
AB219398S. asterophora19United Kingdom
HQ643975S. asterophora19United Kingdom
HQ643406S. megasperma21Japan
HQ643407S. megasperma21Japan
HQ643993S. megasperma21United Kingdom
AB219397S. turfosa22Japan
HQ644012S. turfosa22United Kingdom
HQ644013S. turfosa22Germany
AB219384S. anisospora23Japan
EU124748S. anisospora23USA
EU240104S. anisospora23Poland
EU292729S. brachydanionis24China
AB219392Saprolegnia sp. 825Japan
AM228847Saprolegnia sp. 926Sweden
DQ393533Saprolegnia sp. 926USA
DQ393538Saprolegnia sp. 926USA
AB219391Saprolegnia sp. 1027Japan
HQ643992Saprolegnia sp. 1027Japan
HQ438029Saprolegnia sp. 1128Antarctica
HQ643998S. monoica29United Kingdom
HQ643999S. monoica29United Kingdom
EU124749Leptolegnia sp.OutgroupUSA
KF718204Aphanomyces astaciOutgroupSpain
KF717839.1S. parasitica United Kingdom
KF717842.1S. parasitica Poland
KF717860.1S. parasitica Ecuador
KF717864.1S. parasitica Argentina
* As designated by Sandoval-Sierra et al. [14].
Table A2. Reference sequences used in Saprolegnia parasitica phylogenetic analysis. Sequences cover the ITS region and are taken from Sandoval-Sierra et al. [14] and Sarowar et al. [15] Each sample is listed with its respective GenBank accession number, species, geographic origin, and * adopted phylotype nomenclature from [15].
Table A2. Reference sequences used in Saprolegnia parasitica phylogenetic analysis. Sequences cover the ITS region and are taken from Sandoval-Sierra et al. [14] and Sarowar et al. [15] Each sample is listed with its respective GenBank accession number, species, geographic origin, and * adopted phylotype nomenclature from [15].
Accession N°SpeciesOriginPhylotype *
KF717839.1S. parasiticaUnited KingdomS2
KF717842.1S. parasiticaPolandS6
KF717860.1S. parasiticaEcuadorS5
KF717864.1S. parasiticaArgentinaS4
MK849947.1S. parasiticaCanadaS1
MK849946.1S. parasiticaCanadaS2
MK849963.1S. parasiticaCanadaS3
MK849943.1S. parasiticaCanadaS4
Jof 10 00057 g0a1
Figure A1. Z3 epizootic identification model, where (xi) is the observed value, ( X ~ ) is the population median, and (MAD) is the median absolute distribution: all 0 count observations were disregarded from X ~ and MAD.
Figure A1. Z3 epizootic identification model, where (xi) is the observed value, ( X ~ ) is the population median, and (MAD) is the median absolute distribution: all 0 count observations were disregarded from X ~ and MAD.
Jof 10 00057 g0a1
Jof 10 00057 g0a2
Figure A2. Phylogenetic summary of all S. parasitica phylotypes isolated across Scottish Atlantic salmon aquaculture sites alongside model reference sequences. A single consensus sequence was used for isolates with identical ITS sequences, with the number of isolates for each consensus in brackets next to the species name. Isolates from this study are boldened and underlined. Reference sequences were taken from [14], and [15] are listed with their NCBI accession number and species names. Phy maximum-likelihood tree and Bayesian inference were both produced using the K80 + G model. Supporting values for each branch are displayed with maximum-likelihood bootstrap values below and Bayesian support values above (bold).
Figure A2. Phylogenetic summary of all S. parasitica phylotypes isolated across Scottish Atlantic salmon aquaculture sites alongside model reference sequences. A single consensus sequence was used for isolates with identical ITS sequences, with the number of isolates for each consensus in brackets next to the species name. Isolates from this study are boldened and underlined. Reference sequences were taken from [14], and [15] are listed with their NCBI accession number and species names. Phy maximum-likelihood tree and Bayesian inference were both produced using the K80 + G model. Supporting values for each branch are displayed with maximum-likelihood bootstrap values below and Bayesian support values above (bold).
Jof 10 00057 g0a2

References

  1. van West, P. Molecular Tools to Unravel the Role of Genes from Phytophthora infestans; Wageningen University and Research: Wageningen, The Netherlands, 2000. [Google Scholar]
  2. Bruno, D.W.; van West, P.; Beakes, G.W. Saprolegnia and Other Oomycetes. In Fish Diseases and Disorders. Volume 3: Viral, Bacterial and Fungal Infections; CABI International: Wallingford, UK, 2011; pp. 669–720. [Google Scholar]
  3. van den Berg, A.H.; McLaggan, D.; Diéguez-Uribeondo, J.; van West, P. The Impact of the Water Moulds Saprolegnia diclina and Saprolegnia parasitica on Natural Ecosystems and the Aquaculture Industry. Fungal Biol. Rev. 2013, 27, 33–42. [Google Scholar] [CrossRef]
  4. Sandoval-Sierra, J.V.; Latif-Eugenin, F.; Martín, M.P.; Zaror, L.; Diéguez-Uribeondo, J. Saprolegnia Species Affecting the Salmonid Aquaculture in Chile and Their Associations with Fish Developmental Stage. Aquaculture 2014, 434, 462–469. [Google Scholar] [CrossRef]
  5. van West, P. Saprolegnia parasitica, an Oomycete Pathogen with a Fishy Appetite: New Challenges for an Old Problem. Mycologist 2006, 20, 99–104. [Google Scholar] [CrossRef]
  6. Phillips, A.J.; Anderson, V.L.; Robertson, E.J.; Secombes, C.J.; van West, P. New Insights into Animal Pathogenic Oomycetes. Trends Microbiol. 2008, 16, 13–19. [Google Scholar] [CrossRef] [PubMed]
  7. Willoughby, L.G. Saprolegnias of Salmonid Fish in Windermere: A Critical Analysis. J. Fish Dis. 1978, 1, 51–67. [Google Scholar] [CrossRef]
  8. Diéguez-Uribeondo, J.; Fregeneda-Grandes, J.M.; Cerenius, L.; Pérez-Iniesta, E.; Aller-Gancedo, J.M.; Tellería, M.T.; Söderhäll, K.; Martín, M.P. Re-Evaluation of the Enigmatic Species Complex Saprolegnia Diclina–Saprolegnia Parasitica Based on Morphological, Physiological and Molecular Data. Fungal Genet. Biol. 2007, 44, 585–601. [Google Scholar] [CrossRef] [PubMed]
  9. Rezinciuc, S.; Sandoval-Sierra, J.-V.; Diéguez-Uribeondo, J. Molecular Identification of a Bronopol Tolerant Strain of Saprolegnia australis Causing Egg and Fry Mortality in Farmed Brown Trout, Salmo trutta. Fungal Biol. 2014, 118, 591–600. [Google Scholar] [CrossRef] [PubMed]
  10. Rezinciuc, S.; Sandoval-Sierra, J.V.; Ruiz-León, Y.; van West, P.; Diéguez-Uribeondo, J. Specialized Attachment Structure of the Fish Pathogenic Oomycete Saprolegnia parasitica. PLoS ONE 2018, 13, e0190361. [Google Scholar] [CrossRef]
  11. Coker, W.C. The Saprolegniaceae, with Notes on Other Water Molds, by William Chambers Coker; University of North Carolina Press: Chapel Hill, NC, USA, 1923. [Google Scholar]
  12. Seymore, R.L. The Genus Saprolegnia. Nova Hedwig. 1970, 19, 1–124. [Google Scholar]
  13. Johnson, T.W.; Seymore, R.L.; Padgett, D.E. Biology and Systematics of the Saprolegniaceae; University of North Carolina at Wilmington, Department of Biological Sciences: Wilmington, NC, USA, 2022. [Google Scholar]
  14. Sandoval-Sierra, J.V.; Martín, M.P.; Diéguez-Uribeondo, J. Species Identification in the Genus Saprolegnia (Oomycetes): Defining DNA-Based Molecular Operational Taxonomic Units. Fungal Biol. 2014, 118, 559–578. [Google Scholar] [CrossRef]
  15. Sarowar, M.N.; Cusack, R.; Duston, J. Saprolegnia Molecular Phylogeny among Farmed Teleosts in Nova Scotia, Canada. J. Fish Dis. 2019, 42, 1745–1760. [Google Scholar] [CrossRef] [PubMed]
  16. Van West, P.; De Bruijn, I.; Minor, K.L.; Phillips, A.J.; Robertson, E.J.; Wawra, S.; Bain, J.; Anderson, V.L.; Secombes, C.J. The Putative RxLR Effector Protein SpHtp1 from the Fish Pathogenic Oomycete Saprolegnia parasitica Is Translocated into Fish Cells. FEMS Microbiol. Lett. 2010, 310, 127–137. [Google Scholar] [CrossRef] [PubMed]
  17. Zelaya-Molina, L.X.; Ortega, M.A.; Dorrance, A.E. Easy and Efficient Protocol for Oomycete DNA Extraction Suitable for Population Genetic Analysis. Biotechnol. Lett. 2011, 33, 715–720. [Google Scholar] [CrossRef]
  18. White, T.J.; Bruns, T.; Lee, S.; Taylor, J. Amplification and Direct Sequencing of Fungal Ribosomal RNA Genes for Phylogenetics. In PCR Protocols; Academic Press: Cambridge, MA, USA, 1990; pp. 315–322. [Google Scholar]
  19. Sarowar, M.N.; Hossain, M.J.; Nasrin, T.; Naznin, T.; Hossain, Z.; Rahman, M.M. Molecular Identification of Oomycete Species Affecting Aquaculture in Bangladesh. Aquac. Fish 2019, 4, 105–113. [Google Scholar] [CrossRef]
  20. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic Local Alignment Search Tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef] [PubMed]
  21. Edgar, R.C. MUSCLE: A Multiple Sequence Alignment Method with Reduced Time and Space Complexity. BMC Bioinform. 2004, 5, 113. [Google Scholar] [CrossRef] [PubMed]
  22. Milne, I.; Lindner, D.; Bayer, M.; Husmeier, D.; McGuire, G.; Marshall, D.F.; Wright, F. TOPALi v2: A Rich Graphical Interface for Evolutionary Analyses of Multiple Alignments on HPC Clusters and Multi-Core Desktops. Bioinformatics 2009, 25, 126–127. [Google Scholar] [CrossRef]
  23. Guindon, S.; Gascuel, O. A Simple, Fast, and Accurate Algorithm to Estimate Large Phylogenies by Maximum Likelihood. Syst. Biol. 2003, 52, 696–704. [Google Scholar] [CrossRef]
  24. Ronquist, F.; Huelsenbeck, J.P. MrBayes 3: Bayesian Phylogenetic Inference under Mixed Models. Bioinformatics 2003, 19, 1572–1574. [Google Scholar] [CrossRef]
  25. Stöver, B.C.; Müller, K.F. TreeGraph 2: Combining and Visualizing Evidence from Different Phylogenetic Analyses. BMC Bioinform. 2010, 11, 7. [Google Scholar] [CrossRef]
  26. Ali, S.E.; Thoen, E.; Vrålstad, T.; Kristensen, R.; Evensen, Ø.; Skaar, I. Development and Reproduction of Saprolegnia Species in Biofilms. Vet. Microbiol. 2013, 163, 133–141. [Google Scholar] [CrossRef] [PubMed]
  27. Beckmann, M.J. Saprolegnia Infections after Vaccination in Scottish Salmon Aquaculture: A Host-Pathogen Interaction Influenced by Stress. Ph.D. Thesis, University of Aberdeen, Aberdeen, UK, 2016. [Google Scholar]
  28. Matthews, E. Environmental Factors Impacting Saprolegnia Infections in Wild Fish Stocks. Ph.D. Thesis, Cardiff University, Cardiff, Wales, 2019. [Google Scholar]
  29. Pavić, D.; Miljanović, A.; Grbin, D.; Šver, L.; Vladušić, T.; Galuppi, R.; Tedesco, P.; Bielen, A. Identification and Molecular Characterization of Oomycete Isolates from Trout Farms in Croatia, and Their Upstream and Downstream Water Environments. Aquaculture 2021, 540, 736652. [Google Scholar] [CrossRef]
  30. Elgendy, M.Y.; Ali, S.E.; Abdelsalam, M.; Abd El-Aziz, T.H.; Abo-Aziza, F.; Osman, H.A.; Authman, M.M.N.; Abbas, W.T. Onion (Allium cepa) Improves Nile Tilapia (Oreochromis niloticus) Resistance to Saprolegniasis (Saprolegnia parasitica) and Reduces Immunosuppressive Effects of Cadmium. Aquac. Int. 2023, 31, 1457–1481. [Google Scholar] [CrossRef]
  31. Songe, M.M.; Willems, A.; Wiik-Nielsen, J.; Thoen, E.; Evensen, Ø.; West, P.; Skaar, I. Saprolegnia Diclina IIIA and S. parasitica Employ Different Infection Strategies When Colonizing Eggs of Atlantic Salmon, Salmo salar L. J. Fish Dis. 2016, 39, 343–352. [Google Scholar] [CrossRef] [PubMed]
  32. Ravasi, D.; De Respinis, S.; Wahli, T. Multilocus Sequence Typing Reveals Clonality in Saprolegnia parasitica Outbreaks. J. Fish Dis. 2018, 41, 1653–1665. [Google Scholar] [CrossRef] [PubMed]
  33. Robideau, G.P.; De Cock, A.W.A.M.; Coffey, M.D.; Voglmayr, H.; Brouwer, H.; Bala, K.; Chitty, D.W.; Désaulniers, N.; Eggertson, Q.A.; Gachon, C.M.M.; et al. DNA Barcoding of Oomycetes with Cytochrome c Oxidase Subunit I and Internal Transcribed Spacer. Mol. Ecol. Resour. 2011, 11, 1002–1011. [Google Scholar] [CrossRef] [PubMed]
  34. Stueland, S.; Hatai, K.; Skaar, I. Morphological and Physiological Characteristics of Saprolegnia Spp. Strains Pathogenic to Atlantic Salmon, Salmo salar L. J. Fish Dis. 2005, 28, 445–453. [Google Scholar] [CrossRef]
  35. Nam, B.; Nguyen, T.T.T.; Lee, H.B.; Park, S.K.; Choi, Y.-J. Uncharted Diversity and Ecology of Saprolegniaceae (Oomycota) in Freshwater Environments. Mycobiology 2022, 50, 326–344. [Google Scholar] [CrossRef]
  36. Hulvey, J.P.; Padgett, D.E.; Bailey, J.C. Species Boundaries within Saprolegnia (Saprolegniales, Oomycota) Based on Morphological and DNA Sequence Data. Mycologia 2007, 99, 421–429. [Google Scholar] [CrossRef]
  37. Kozubíková-Balcarová, E.; Koukol, O.; Martín, M.P.; Svoboda, J.; Petrusek, A.; Diéguez-Uribeondo, J. The Diversity of Oomycetes on Crayfish: Morphological vs. Molecular Identification of Cultures Obtained While Isolating the Crayfish Plague Pathogen. Fungal Biol. 2013, 117, 682–691. [Google Scholar] [CrossRef]
  38. Kiesecker, J.M.; Blaustein, A.R.; Miller, C.L. Transfer of a Pathogen from Fish to Amphibians. Conserv. Biol. 2001, 15, 1064–1070. [Google Scholar] [CrossRef]
  39. Yuasa, K.; Hatai, K. Relationship between Pathogenicity of Saprolegnia Spp. Isolates to Rainbow Trout and Their Biological Characteristics. Fish Pathol. 1995, 30, 101–106. [Google Scholar] [CrossRef]
  40. Matthews, E.; Ellison, A.; Cable, J. Saprolegnia parasitica Zoospore Activity and Host Survival Indicates Isolate Variation in Host Preference. Fungal Biol. 2021, 125, 260–268. [Google Scholar] [CrossRef] [PubMed]
  41. Haas, B.J.; Kamoun, S.; Zody, M.C.; Jiang, R.H.Y.; Handsaker, R.E.; Cano, L.M.; Grabherr, M.; Kodira, C.D.; Raffaele, S.; Torto-Alalibo, T.; et al. Genome Sequence and Analysis of the Irish Potato Famine Pathogen Phytophthora Infestans. Nature 2009, 461, 393–398. [Google Scholar] [CrossRef] [PubMed]
  42. Ruthig, G. Water Molds of the Genera Saprolegnia and Leptolegnia Are Pathogenic to the North American Frogs Rana catesbeiana and Pseudacris crucifer, Respectively. Dis. Aquat. Organ. 2009, 84, 173–178. [Google Scholar] [CrossRef] [PubMed]
  43. Khulbe, R.D.; Bisht, G.S.; Joshi, C. Epizootic Infection Due to Achlya debaryana in a Catfish. Mycoses 1994, 37, 61–63. [Google Scholar] [CrossRef] [PubMed]
  44. Hussein, M.; Hassan, W.H.; Mahmoud, M.A. Pathogenicity of Achlya proliferoides and Saprolegnia diclina (Saprolegniaceae) Associated with Saprolegniosis Outbreaks in Cultured Nile Tilapia (Oreochromis niloticus). World J. Fish Mar. Sci. 2013, 5, 188–193. [Google Scholar]
  45. van der Plaats-Niterink, A.J. Monograph of the Genus Pythium. Stud. Mycol. 1981, 21, 1–239. [Google Scholar]
  46. Martin, F.N.; Loper, J.E. Soilborne Plant Diseases Caused by Pythium Spp.: Ecology, Epidemiology, and Prospects for Biological Control. CRC Crit. Rev. Plant Sci. 1999, 18, 111–181. [Google Scholar] [CrossRef]
  47. Mendoza, L.; Vilela, R. The Mammalian Pathogenic Oomycetes. Curr. Fungal Infect. Rep. 2013, 7, 198–208. [Google Scholar] [CrossRef]
  48. Parveen, T.; Sharma, K. Pythium Diseases, Control and Management Strategies: A Review. Int. J. Plant Anim. Environ. Sci. 2015, 5, 244–2547. [Google Scholar]
  49. El-Tarabily, K.A.; Hardy, G.E.S.J.; Sivasithamparam, K.; Hussein, A.M.; Kurtböke, D.I. The Potential for the Biological Control of Cavity-spot Disease of Carrots, Caused by Pythium coloratum, by Streptomycete and Non-streptomycete Actinomycetes. New Phytol. 1997, 137, 495–507. [Google Scholar] [CrossRef]
  50. Miura, M.; Hatai, K.; Tojo, M.; Wada, S.; Kobayashi, S.; Okazaki, T. Visceral Mycosis in Ayu Plecoglossus Altivelis Larvae Caused by Pythium flevoense. Fish Pathol. 2010, 45, 24–30. [Google Scholar] [CrossRef]
  51. Czeczuga, B. Species of Pythium Isolated from Eggs of Fresh-Water Fish. Acta Mycol. 2014, 31, 151–161. [Google Scholar] [CrossRef]
  52. Mircetich, S.M.; Kraft, J.M. Efficiency of Various Selective Media in Determining Pythium Population in Soil. Mycopathol. Mycol. Appl. 1973, 50, 151–161. [Google Scholar] [CrossRef] [PubMed]
  53. Tsao, P.H. Inhibition of Mortierella and Pythium in a Phytophthora-Isolation Medium Containing Hymexazol. Phytopathology 1977, 77, 796. [Google Scholar] [CrossRef]
  54. Robertson, G.I. The Genus Pythium in New Zealand. N. Z. J. Bot. 1980, 18, 73–102. [Google Scholar] [CrossRef]
  55. Vadivelan, G.; Venkateswaran, G. Production and Enhancement of Omega-3 Fatty Acid from Mortierella alpina CFR-GV15: Its Food and Therapeutic Application. BioMed Res. Int. 2014, 2014, 657414. [Google Scholar] [CrossRef]
  56. Ozimek, E.; Jaroszuk-Ściseł, J.; Bohacz, J.; Korniłłowicz-Kowalska, T.; Tyśkiewicz, R.; Słomka, A.; Nowak, A.; Hanaka, A. Synthesis of Indoleacetic Acid, Gibberellic Acid and ACC-Deaminase by Mortierella Strains Promote Winter Wheat Seedlings Growth under Different Conditions. Int. J. Mol. Sci. 2018, 19, 3218. [Google Scholar] [CrossRef]
  57. Johnson, J.M.; Ludwig, A.; Furch, A.C.U.; Mithöfer, A.; Scholz, S.; Reichelt, M.; Oelmüller, R. The Beneficial Root-Colonizing Fungus Mortierella hyalina Promotes the Aerial Growth of Arabidopsis and Activates Calcium-Dependent Responses That Restrict Alternaria brassicae–Induced Disease Development in Roots. Mol. Plant-Microbe Interact. 2019, 32, 351–363. [Google Scholar] [CrossRef]
Jof 10 00057 g001
Figure 1. Phylogenetic summary of all oomycete and fungal species isolated using oomycete-selective media, across Scottish Atlantic salmon aquaculture sites. A single consensus ITS sequence was used for isolates belonging to the same species, with the number of isolates for each species in brackets next to the species name. Phy maximum-likelihood tree was produced using the TrN + I + G model with Bayesian inference using the GTR + I + G model. Supporting values for each branch are displayed with maximum-likelihood bootstrap values below and Bayesian support values above (bold). Conflicting Bayesian topologies are indicated with posterior probabilities in square brackets.
Figure 1. Phylogenetic summary of all oomycete and fungal species isolated using oomycete-selective media, across Scottish Atlantic salmon aquaculture sites. A single consensus ITS sequence was used for isolates belonging to the same species, with the number of isolates for each species in brackets next to the species name. Phy maximum-likelihood tree was produced using the TrN + I + G model with Bayesian inference using the GTR + I + G model. Supporting values for each branch are displayed with maximum-likelihood bootstrap values below and Bayesian support values above (bold). Conflicting Bayesian topologies are indicated with posterior probabilities in square brackets.
Jof 10 00057 g001
Jof 10 00057 g002
Figure 2. Phylogenetic summary of all Saprolegnia species isolated across Scottish Atlantic salmon aquaculture sites alongside model reference sequences (Appendix A, Table A1). A single consensus sequence was used for isolates with identical ITS sequences, with the number of isolates for each consensus in brackets next to the species name. Isolates from this study are boldened and underlined. Reference sequences were taken from [14] and are listed with their NCBI accession number, species names, and cluster number as designated by the reference study. Phy maximum-likelihood tree and Bayesian inference were both produced using the GTR + G model. Supporting values for each branch are displayed with maximum-likelihood bootstrap values below and Bayesian support values above (bold).
Figure 2. Phylogenetic summary of all Saprolegnia species isolated across Scottish Atlantic salmon aquaculture sites alongside model reference sequences (Appendix A, Table A1). A single consensus sequence was used for isolates with identical ITS sequences, with the number of isolates for each consensus in brackets next to the species name. Isolates from this study are boldened and underlined. Reference sequences were taken from [14] and are listed with their NCBI accession number, species names, and cluster number as designated by the reference study. Phy maximum-likelihood tree and Bayesian inference were both produced using the GTR + G model. Supporting values for each branch are displayed with maximum-likelihood bootstrap values below and Bayesian support values above (bold).
Jof 10 00057 g002
Jof 10 00057 g003
Figure 3. Nucleotide sequences of the ITS region for S. parasitica phylotypes. Nucleotides are numbered in accordance with the reference sequence KF717839.1 [14]. Nucleotides identical to the reference sequence for each position are represented by dots; the nucleotides which deviate from the reference are shown. Deviations are seen at the 53, 56, 152, 158, 362, 388, 423, 465, 589, and 598 bp positions.
Figure 3. Nucleotide sequences of the ITS region for S. parasitica phylotypes. Nucleotides are numbered in accordance with the reference sequence KF717839.1 [14]. Nucleotides identical to the reference sequence for each position are represented by dots; the nucleotides which deviate from the reference are shown. Deviations are seen at the 53, 56, 152, 158, 362, 388, 423, 465, 589, and 598 bp positions.
Jof 10 00057 g003
Table 1. Number of isolates of each Saprolegnia species obtained from Atlantic salmon farms in Scotland. Isolates were identified to species level via NCBI BLAST and confirmed through comparison to model sequences identified via molecular operational taxonomic units [14].
Table 1. Number of isolates of each Saprolegnia species obtained from Atlantic salmon farms in Scotland. Isolates were identified to species level via NCBI BLAST and confirmed through comparison to model sequences identified via molecular operational taxonomic units [14].
Aquaculture Site
Saprolegnia SpeciesABCDEFGHIJKLMNTotal
S. australis00000000012000214
S. delica1000000017001010
S. diclina3200403241000120
S. ferax00320000000001033
S. parasitica3135162719713119118151
Table 2. Saprolegnia species phylotypes obtained from Atlantic salmon farms in Scotland. S. parasitica phylotypes were identified via MUSCLE alignment and confirmed via comparison with model sequences published in [14]; S1 was not isolated from any Scottish fish farm.
Table 2. Saprolegnia species phylotypes obtained from Atlantic salmon farms in Scotland. S. parasitica phylotypes were identified via MUSCLE alignment and confirmed via comparison with model sequences published in [14]; S1 was not isolated from any Scottish fish farm.
Aquaculture Site
Saprolegnia PhylotypesABCDEFGHIJKLMNTotal
S21132440015713119118107
S3100000100000002
S4000200000000002
S5100000300000004
S600270270000000036
Table 3. Number of isolates of each Saprolegnia species obtained from water and from fish in Atlantic salmon farms in Scotland. A phylogenetic analysis of these isolates from water and from fish is shown in Supplementary Figure S1.
Table 3. Number of isolates of each Saprolegnia species obtained from water and from fish in Atlantic salmon farms in Scotland. A phylogenetic analysis of these isolates from water and from fish is shown in Supplementary Figure S1.
Saprolegnia SpeciesWaterFishTotal
S. australis8614
S. delica5510
S. diclina12820
S. ferax102333
S. parasitica21130151
S. parasitica phylotypes
S27100107
S3112
S4022
S5044
S6132336
Table 4. Total number of isolates of each Saprolegnia species obtained from enzootic and epizootic periods in Atlantic salmon farms in Scotland. A phylogenetic analysis of these isolates from water and from fish is shown in Supplementary Figure S2.
Table 4. Total number of isolates of each Saprolegnia species obtained from enzootic and epizootic periods in Atlantic salmon farms in Scotland. A phylogenetic analysis of these isolates from water and from fish is shown in Supplementary Figure S2.
Saprolegnia SpeciesEnzooticEpizooticTotal
S. australis8614
S. delica4610
S. diclina16420
S. ferax33033
S. parasitica10345148
S. parasitica phylotypes
S26541106
S3033
S4 202
S5011
S636036
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Shreves, K.V.; Saraiva, M.; Ruba, T.; Miller, C.; Scott, E.M.; McLaggan, D.; van West, P. Specific Phylotypes of Saprolegnia parasitica Associated with Atlantic Salmon Freshwater Aquaculture. J. Fungi 2024, 10, 57. https://doi.org/10.3390/jof10010057

AMA Style

Shreves KV, Saraiva M, Ruba T, Miller C, Scott EM, McLaggan D, van West P. Specific Phylotypes of Saprolegnia parasitica Associated with Atlantic Salmon Freshwater Aquaculture. Journal of Fungi. 2024; 10(1):57. https://doi.org/10.3390/jof10010057

Chicago/Turabian Style

Shreves, Kypher Varin, Marcia Saraiva, Tahmina Ruba, Claire Miller, E. Marian Scott, Debbie McLaggan, and Pieter van West. 2024. "Specific Phylotypes of Saprolegnia parasitica Associated with Atlantic Salmon Freshwater Aquaculture" Journal of Fungi 10, no. 1: 57. https://doi.org/10.3390/jof10010057

APA Style

Shreves, K. V., Saraiva, M., Ruba, T., Miller, C., Scott, E. M., McLaggan, D., & van West, P. (2024). Specific Phylotypes of Saprolegnia parasitica Associated with Atlantic Salmon Freshwater Aquaculture. Journal of Fungi, 10(1), 57. https://doi.org/10.3390/jof10010057

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop