Trematode diversity reflecting the community structure of Danish freshwater systems: molecular clues
Parasites & Vectors volume 14, Article number: 43 (2021)
Digenean trematodes are parasitic platyhelminths that use several hosts in their life cycles and are thereby embedded in various ecosystems affected by local environmental conditions. Their presence in a habitat will reflect the presence of different host species and, as such, they can serve as ecological indicators. Only limited information on the occurrence of trematodes and their link to other trophic levels in the Danish freshwater ecosystems is currently available.Therefore, the main aim of the present study was to increase our knowledge in this field.
Snails were sampled from 21 freshwater lakes in Denmark, following which shedding procedures were performed, cercariae were recoved and the released parasites were identified using molecular tools (PCR and sequencing).
A total of 5657 snail hosts belonging to ten species were identified, revealing a highly diverse parasite fauna comprising 22 trematode species. The overall trematode prevalence was 12.6%, but large variations occurred between host species. The snail host Lymnaea stagnalis showed the highest prevalence and also exhibited the highest diversity, accounting for 47.6% of the species richness.
This survey contributes updated information on parasite–host relations and compatibility and may assist in describing the ecological structure of the investigated Danish freshwater ecosystems.
Trematodes, dominated by the digeneans, comprise a group of abundant and ubiquitous parasitic metazoans with high diversity and more than 18,000 nominal species . These metazoans utilize multi-host systems to complete their complex life cycles, often with a mollusc as the first intermediate host [2,3,4]. Details on the associative life-cycle patterns between 1000 digenean species and their mollusc hosts have been recorded . Other invertebrates and vertebrates may play a role as the second intermediate host, whereas the definitive hosts are often vertebrates, such as various fishes, amphibians, reptiles, birds and mammals . The trematode life cycles involve intricate transmission patterns adapted to the habitat, the abiotic and biotic conditions and especially the different host species present .
Due to the obligate relation to free-living organisms in the ecosystem, these parasites may serve as bioindicators and used to indicate presence of obligate hosts, environmental conditions and disease risks under various conditions (climate changes, anthropogenic impacts) [2, 8,9,10,11,12]. This potential frames the importance of conducting continuous surveys of trematode occurrence in the freshwater environment in any region. However, due to the characteristics of adult digeneans, often hidden in their hosts and with uneven spatial distributions, it is relatively difficult to detect and monitor that stage. In contrast, sampling of the mollusc host in a habitat creates a method for efficient biomonitoring of their larvae. Incubation of an infected snail in the laboratory allows shedding of cercariae and thereby detection of a parasite at several trophic levels and hosts in the environment. The adults in the final hosts possess characteristic structures, but the larval stages related to the intermediate mollusc host (e.g. cercariae, rediae, sporocysts) and even metacercariae in the second intermediate host are generally more difficult to identify . With the advent of molecular tools, such as PCR and subsequent sequencing of selected genes, it is now possible to elucidate the occurrence and diversity of various trematodes species in freshwater systems simply by sampling the mollusc host, recovering the released cercariae, isolating DNA, performing PCR and obtaining sequence information. Selected parts of the parasite genome, such as the nuclear ribosomal DNA regions comprising 18S, ITS1, 5.8S, ITS2, 28S and the mitochondrial cytochrome c oxidase subunit genes (cox1 and cox2) genes may be used to identify trematode species across a broad taxonomic range . Diagnosis can be done by performing BLAST analysis against well-described species in GenBank.
Almost one century ago, the limnologist Wesenberg-Lund investigated the occurrence of cercariae in snails in a few Danish freshwater systems . However, with a few exceptions, the link between the cercariae and other organisms in the environment was at that time poorly known. In addition, the identity of species of the recovered cercariae and their connection to other hosts in the ecosystem was in most cases not known. By using a molecular approach it is now possible to supplement our knowledge on the fauna of trematodes in Danish freshwater lakes [4, 16]. The present study describes the diversity of trematodes in 21 Danish freshwater lakes (geographically representing a major part of the national territory) by investigating ten species of molluscs and their released cercariae. By using current knowledge on transmission pathways, we obtained information on second intermediate and definitive host occurrence and trophic interactions in the ecosystems.
Samples of snails were recovered from a total of 21 selected freshwater lakes (sø) in Zealand (eastern part of Denmark) and Jutland (Western part of Denmark). The precise locations of the lakes were determined using the global positioning system (Table 1), and samplings were conducted from March 2019 to September 2019. Snails were collected by hand in shallow water areas, from stones and aquatic plants, placed in 1-l containers in cooled transport boxes and transferred (maximum 4 h transport time) to the university laboratory (Frederiksberg, Copenhagen area) for subsequent shedding procedures. In order to evaluate the effect of season (sampling date) on infection, three lakes (Bagsværd sø, Furesø and Esrum sø) were sampled several times between March and September 2019.
Cercarial shedding from snails
Once in the laboratory, the sampled snails were rinsed in de-chlorinated municipal water for removal of surface-adhering organisms (preventing interference of false positives), whereafter they were identified to species by shell morphology . Individual snails were placed at room temperature in 100-ml beakers containing 50 ml of 0.45-µm filtered lake water (Filtropur; Sarstedt, Nümbrecht, Germany)  and incubated for 8 h in daylight. The beakers were subsequently examined under a stereomicroscope (magnification ×7–40; Leica Microsystems GmbH, Wetzlar, Germany) to record presence/absence of cercariae. Infected snails were characterized as “positive” and non-infected snails as “negative”, following which the cercarial type was determined [19,20,21,22]. Recovered cercariae were preserved in 96% ethanol (Kemetyl A/S, Køge, Denmark) for subsequent molecular analysis .
Molecular identification of cercariae
Individual ethanol-fixed cercariae were separated and placed in Eppendorf tubes (Axygen®, Sorenson BioScience, Salt Lake City, UT, USA), following which the ethanol was evaporated (Eppendorf ThermoMixer® Comfort; Eppendorf, Hamburg, Germany). DNA from each specimen was then isolated following lysis (using a guanidine thiocyanate lysis buffer) (Tkach and Pawlowski ). DNA was purified from the PCR mixture using Illustra™ GFX™PCR and Gel band purification kit (VWR International, Søborg, Denmark). DNA concentration and purity were measured using a Nanodrop 2000 spectrophotometer (Saveen & Werner ApS, Jyllinge, Denmark). The complete ITS1–5.8S–ITS2 regions of the rDNA were amplified using forward primer BD1 (5′-GTC GTA ACA AGG TTT CCG TA-3′) and reverse primer BD2 (5′-TAT GCT TAA ATT CAG CGG GT-3′)  under PCR cycling conditions as previously described . Each sequence obtained from the PCR conducted was subjected to BLAST analysis at NCBI (National Center for Biotechnology Information). All sequences were submitted to GenBank and obtained accession numbers.
Evolutionary analysis by maximum likelihood method
Evolutionary analyses were conducted in MEGA X . Clustal W was used to align the sequences, and the model General Time Reversible G+I using the Akaike criterion was the best fitting model of the 24 models tested by MEGA X. Therefore, the evolutionary history was inferred by using the maximum likelihood method and general time-reversible model. The tree with the highest log likelihood (− 31524, 07) is shown in the Electronic Supplementary Material. The percentage of trees in which the associated taxa clustered together is shown next to the branches. Initial tree(s) for the heuristic search were obtained automatically by applying the neighbor-joining and BioNJ algorithms to a matrix of pairwise distances estimated using the maximum composite likelihood (MCL) approach, and then selecting the topology with the superior log likelihood value. A discrete gamma distribution was used to model evolutionary rate differences among sites (5 categories (+ G, parameter = 15143)). The rate variation model allowed for some sites to be evolutionarily invariable ([+ I], 9.91% sites). This analysis involved 89 nucleotide sequences. There were a total of 1960 positions in the final dataset.
Raw data were entered into a Microsoft Excel (Microsoft Corp., Redmond, WA, USA) spreadsheet, and descriptive statistics were used to summarize the data. A Chi-square test was used to assess any significant differences with regard to prevalence in different seasons. Pearson's correlation analysis was used to assess correlation between the number of trematode species and snail species in each lake. The data were analyzed statistically using the RStudio software program (v.3.44; R Foundation for Statistical Computing, Vienna, Austria). The software Tableau (v.2019.3) was applied for community structure (biodiversity) analysis, and Gephi v.0.9.2 was used to illustrate the interconnection between trematode, infected snail and parasite–host transmission pattern. Statistical significance was set at P < 0.05 for all statistical tests.
Sample sizes at the 21 different locations are indicated in Fig. 1. Large samples were recovered in Esrum sø (1030), Furesø (852), Bagsværd (695) and Skanderborg sø (405), whereas relatively small samples were collected in Utterslev mose (81), Almind sø (18) and Mossø (11). Sixteen lakes were positive with regard to trematode infection (cercarial shedding), whereas the remaining five locations were negative for parasites.
Prevalence of infection at different locations
We examined a total of 5657 individual freshwater snails belonging to three families and ten different species: Lymnaeidae (Lymnaea stagnalis, Radix balthica, Radix auricularia, Stagnicola palustris and Omphiscola glabra); Planorbidae (Planorbarius corneus, Planorbis planorbis and Anisus vortex); and Bithyniidae (Bithynia leachii and Bithynia tentaculata). The samples were dominated by Anisus vortex (1402), Radix balthica (1200) and Lymnaea stagnalis (1173). Of the 5657 snails examined, 715 (eight species) released cercariae, resulting in an overall prevalence of 12.6%. The most prevalent snail overall was Lymnaea stagnalis (18.2%), followed by Radix balthica (14.4%). Seven snail species were collected in Furesø and six species in Bagsværd sø, Esrum sø and Skanderborg sø (Table 2), with only one snail species collected from each of the other lakes, such as the artificial lake Ringen sø. Sixteen lakes were positive and five lakes were negative for trematode infection, but some lakes showed a relatively high prevalence of infection. Thus, Ringen sø showed a 46.7% prevalence of trematode infection, Bromme Lille sø 29.0% and Tjele Langsø 27.9%. A number of other lakes showed an intermediate prevalence level of between 5 and 23% (Almind sø, Bagsværd sø, Furesø, Lyngby sø, Tissø, Skanderborg sø), and a few lakes showed a low prevalence of < 5% (Sjælsø, Haraldsted sø, Buresø). Since the 21 lakes were monitored at different sampling dates, it was considered important to test for a possible occurrence of a seasonal variation which could affect the overall data set (Table 3). To this end, we sampled three different lakes as an internal reference (Bagsværd sø, Furesø and Esrum sø) three times from March to September. During these samplings, we collected a total of 119 snails in March, 900 in July and 386 in September. A Chi-square test showed that prevalence of infection on these sampling dates was significantly different (P < 0.05), with the highest prevalence appearing in September (48.96%) and lowest prevalence showed in March (1.68%).
Sequences of the rDNA locus ITS of the different parasite isolates were submitted to GenBank. In order to identify the trematodes, both percentage identity obtained from a BLAST search of GenBank and a phylogenetic analysis were used. Sequenced specimens all belonged to class Digenea (GenBank acc. nos. MW000945 to MW001148), and of these 52 sequences were unique. The ITS1–5.8S–ITS2 part of these 52 sequences were subjected to a phylogenetic analysis together with 34 other digenean species. Three species of the sister class Aspidogastrea were used as an outgroup (marked in bold on Additional file 1: Fig. S1). A total of 22 trematode species belonging to nine families were recorded (Additional file 2: Table S1). The diversity of trematode species differed markedly between snail hosts. Some species of snails exhibited a rich trematode fauna. Lymnaea stagnalis was found to be infected by ten species and showed the highest overall prevalence of infection (18.2%). In some snail species, such as Planorbis planorbis and Bithynia tentaculata, only one trematode species was recovered.
The host snail Lymnaea stagnalis was found to be infected by several types of digeneans belonging to the family Plagiorchiidae. Species within the genus Plagiorchis were the most frequently recovered, showing an overall prevalence of 5.0%. This species was also prevalent in spring samples despite a generally lower infection rate in this season. Three species with a zoonotic potential were recovered, including Isthmiophora melis (found in Skanderborg sø), Trichobilharzia sp. (found in Bagsværd sø, Tjele langsø and Furesø) and Metorchis orientalis (found in Lyngby sø and Skanderborg sø).
The parasites identified utilize a wide variety of second intermediate and definitive hosts (based on published data from the literature; Additional file 2: Table S1) falling within six transmission units. One of the entities are parasites utilizing fishes as the second intermediate host and fish-eating birds as the definitive host (seven species). This group of trematodes is mainly associated with Radix balthica and dominated the digenean fauna in BrommeLille sø. Another group uses molluscs as the second intermediate host and waterfowl as the definitive host (two species). A third group utilizes larval insects as second intermediate host and bats as the definitive host (one species), whereas a fourth group uses frogs (two species) or waterfowl (one genus) as the final host. The fifth group utilizes leeches as the second intermediate host and waterfowl as the definitive host (2 species). Finally, a sixth group comprises species with multiple intermediate and definitive hosts (four species), including larval insects, amphipods, molluscs, leeches and mammals.
Parasites reflecting host occurrence in a lake
Digeneans are relatively host specific, especially at the mollusc host level, and each species utilizes generally one or few related hosts at each stage of its life cycle. Based on published information on hosts and life cycles, we constructed a list of possible animal species present in each lake (Additional file 2: Table S1). Lakes in Zealand and Jutland did not differ in terms of host occurrence inferred from the trematode species present (Fig. 2). Freshwater fish, fish-eating birds, larval insects, amphipods, molluscs, various birds and mammals were host types deducted from the parasite lists in both the eastern and western locations. Primarily BrommeLille sø in Zealand and secondly Tjele langsø in Jutland were populated with parasites, indicating a high abundance of freshwater fish and fish-eating birds. Snails collected from Ringen sø shed parasites, indicating larval insects, amphipods and molluscs, birds and mammals. The presence of waterfowl, tadpoles, leeches and bats were indicated both in Zealand and Jutland. In Zealand (the eastern part of Denmark), a number of host types were indicated, with a higher host diversity in Lyngby sø (9 types), Esrum sø (7 types), Bagsværd sø (7 types) and Furesø (5 types). These were followed, in terms of host diversity, by Haraldsted sø, Tissø and Sjælsø. In Jutland (western part of Denmark) the highest number of host types were seen in Skanderborg sø (7 types) followed by Tjele langsø (6 types).
Investigations on the trematode fauna of freshwater snails in a certain habitat can contribute to our knowledge and understanding of the ecology of trematodes. A basic knowledge of life cycles with the inclusion of obligate intermediate and/or final host(s) will reflect the existence of other faunal elements in the habitat. The most recent study on Danish freshwater systems involving cercariae released by snails was conducted almost one century ago by Wesenberg-Lund  in selected Danish water bodies. Since then considerable research efforts have been invested in the elucidation of life cycles, and this knowledge in this field can be used to elevate our overall understanding of the freshwater biotope. In our study, we conducted a comprehensive and expanded investigation on cercariae released from various snail hosts in Danish lakes. This will serve as an updated record of the trematode fauna in Denmark, link the snail parasites to other trophic levels and facilitate evaluations of the freshwater systems involved. Recent studies in other European regions and in Lake Victoria in Afrtica  provide a valuable basis for comparison. Selbach et al.  examined 5347 snails sampled in German lakes belonging to six species which were infected with a total of 36 trematode species. Similar studies were conducted at five aquatic locations of Europe during the latest decades (Table 4). Our samplings comprised 5657 specimens of snails belonging to ten species within three families, and we found 715 snails infected with parasites, resulting in an overall prevalence of 12.6%. This level of prevalence and the number of species (22) in Danish lakes correspond well with several records from other freshwater bodies in Europe and Africa.
We selected a number of freshwater lakes for our investigation with the aim to cover both the western (represented by Jutland) and eastern (represented by Zealand) parts of Denmark. The intentions was to collect a high number of snails at all locations, but the sample sizes were low at some localities and on certain sampling dates (Almind sø and Mossø). Low sample sizes may be explained by the practical and logistic challenges associated with limited access to private property, physical inaccessibility, elevated water levels and inferior weather conditions, which in turn may bias observations when different lakes are compared. However, the numbers were sufficient to discuss the biological implications when finding a certain type (genus, species) of cercaria released from freshwater snails in a specific lake. Thus, the presence of obligate intermediate and final hosts is indicated if a parasite is found in the ecosystem. We do not at present have sufficient information on each lake’s history to draw conclusions on causative factors (abiotic, biotic, anthropogenic). However, it is known that both abiotic and biotic factors influence the occurrence of snails and parasites, that the pond area of eutrophic water bodies influences the occurrence of freshwater snails  and that Bithynia tentaculata is affected by human recreational transport (proximity to a boat launch) and eutrophication (anthropogenic land use, elevated nutrient concentrations) . Future analyses should therefore include advanced and comprehensive physico-chemical lake data collection to allow advance the conclusions drawn from the present study. This may allow us to implement corresponding studies on how these factors theoretically could affect parasites in the infected snails in Danish lakes as well.
It is noteworthy that the minimal snail richness in Ringen sø was associated with a relatively high prevalence of trematode infection (46.7%). Corresponding high infection levels have been observed previously, such as a 47.1% prevalence of infection in Lymnaea stagnalis in the Czech Republic . It may be that the local environmental conditions Ringen sø account for this high prevalence (46.7%). This leads to the speculation that this artificial and young lake (established 1977), which generally has a low biodiversity, provides optimal conditions for the propagation of a single snail species and a few associated parasite species following their first introduction with local mallards and swans.
We showed that the prevalence of infection showed a seasonal variation—at least based on the three representative Danish lakes sampled from March to September. Therefore, it is not reasonable to directly compare different lakes which were sampled at different time points. Our analyses will therefore not include overall correlations between abiotic conditions and infection levels, although Soldánová  did show that the infection probability of L. stagnalis was associated with levels of dissolved oxygen and pH of lakes.
The most abundant snail species in our study was Anisus vortex (sample size 1402), followed by Radix balthica (1200) and Lymnaea stagnalis (1173). The trematode diversity was highest in Lymnaea stagnalis (10 species), Radix balthica (10 species) and Stagnicola palustris (7 species); these results correspond partly to those of a previous study conducted by Selbach et al.  showing that snails within the genera Radix (Lymnaeidae) and Gyraulus (Planorbidae) accounted for almost 90% of the trematode diversity. Similar studies from Canada and Scandinavia also point to high prevalences and trematode diversity, especially in lymnaeids . When evaluating the importance of snail hosts in Danish lakes, we show here that the three species L. stagnalis, R. balthica and S. palustris play a central role for trematode infections in the ecosystem.
A total of 22 trematode species belonging to nine families were found in ten species of snails. It may be argued that the trematode diversity may be influenced by the host sample size, a possibility referred in previous studies showing that a higher number of parasite species was collected from the more abundant snail species . In the present study we collected almost equal numbers of the snails R. balthica and L. stagnalis, and both of these hosts carried ten digenean species, corresponding to data reported in a recent Norwegian study . The central role of lymnaeids for trematode diversity in freshwater lakes was further framed by a high diversity (23 species) in R. auricularia in Germany , although a lower diversity may be found in the Czech Republic (three species)  and southeast Germany (four trematode species) .
In a German study R. auricularia harbored six trematode species of family Diplostomidae and contributed to a high diversity of trematode species , which is consistent with our study. In addition to identifying Diplostomum baeri, D. mergi and Tylodelphys clavata, the authors of that study also found D. parviventosum and D. spathaceum, compared with Ornithodiplostomum scardinii, Posthodiplostomum cuticola and D. pseudospathaceum released from Radix balthica in our study. Our study is the first to report the presence of O. scardinii cercariae, but the adult specimens have been isolated from infected birds in the Czech Republic and Poland .
Another investigation of biodiversity of trematodes in freshwater ecosystems of Europe  reported that species within the families Diplostomidae, Echinostomatidae and Schistosomatidae were the most frequently recorded. The most frequently recorded trematode species in these studies is Diplostomum spathaceum, which was recorded in the Czech Republic , UK , Finland , Lithuania , Poland , Russia  and Switzerland . Another species, Diplostomum pseudospathaceum, was found in Switzerland , UK , Czech Republic , Finland , Poland , Lithuania  and Germany . A third common species is Echinoparyphium recurvatum found in the Czech Republic [38, 48], Poland , Germany , Finland , Iceland , Ukraine  and Russia .
The findings of our study in Denmark correspond to a certain extent to these trends. Our study shows that Diplostomidae (six species), Echinostomatidae (five species) and Plagiorchiidae (four species) were the most frequently recorded trematode families. Various species within the genus Plagiorchis were frequently isolated in the present study, corresponding to occurrence reported in Poland , Czech Republic , Finland  and Germany . Lecithodendrium linstowi was the most frequently recorded species in our study (Fig. 3; Table 4).
Comparison of the latest Danish records with the results of studies on freshwater snail–trematode associations carried out in Europe since the beginning of the 20th century revealed some trends. Data from the Czech Republic indicate high a similarity in terms of trematode diversity, with 14 species (66.7%) of trematodes overlappingwith our study. This was followed by Germany (10 species, 47.6%), Ukraine (8 species, 38.1%) and Poland (7 species, 33.3%)  (Fig. 4). Three species, Echinostoma revolutum, Australapatemon sp. and Tylodelphys clavata were the most universally recorded parasites (Fig. 4).
In our study, three types of zoonotic trematodes were discovered, Isthmiophora melis, Trichobilharzia sp. and Metorchis orientalis, which parasitized L. stagnalis, R. balthica and S. palustris, respectively, all snail hosts belonging to the family Lymnaeidae (Table 5). The first intermediate host of Isthmiophora melis is L. stagnalis , and amphibians and fish serve as second intermediate hosts [53, 54]. In China, raw loach meat was reported to be a source of human infection, with adult specimens parasitizing the human intestine and causing diarrhea [55, 56]. Isthmiophora melis has not been reported in many countries, with the exception of Germany  and Czech Republic. Various bird schistosome species within the genus Trichobilharzia use only one intermediate host (pulmonate snails) [58,59,60,61]. Species within the genus Radix may serve as intermediate hosts for Trichobilharzia spp., but they are considered less specific hosts compared to L. stagnalis . Cercariae shed from the snail may accidentially penetrate human skin, instead of their natural target (bird skin), and may within a few minutes after elicit cercarial dermatitis [63, 64]. Species within the genus Trichobilharzia have been reported in numerous locations including, among others Czech Republic , Poland , Finland , Russia  and Germany .
Digenean species within the family Opisthorchiidae have been less well studied in Europe, and data from our study adds to our knowledge on the presence and life cycle of Metorchis orientalis. Species of genus Metorchis have been previously recorded in Finland , and numerous human cases are known in Asia (e.g. China and Japan), where Parafossarulus striatulus, Bulimus striatulus [66,67,68] and Bythinia sp.  can act as the first intermediate host. Adult trematodes commonly inhabit the bile duct and gallbladder of mammals, including humans, causing hepatitis and, in severe cases, even death [70, 71]. Infection is acquired by eating raw freshwater fishes , and susceptible fish host species do occur in Denmark. It is therefore relevant to note that the local lymnaeid snails are potential hosts for this type of zoonotic parasites in Danish freshwater lakes.
Trematodes of the family Plagiorchiidae are ubiquitous parasites and widespread in different host species , as reflected by their common occurrence in Danish snails (Fig. 3). We recorded 12 lakes positive for genera within the family Plagiorchiidae. Based on the findings of the present study, it is clear that the genus Plagiorchis is prevalent in Lymnaea stagnalis, with a prevalence of up to 6.4% (for types within the genus Plagorchis which were not identified to the species level). For the two species P. elegans and P. maculosus, we found a prevalence of 1.5 and 2.8%, respectively (Additional file 2: Table S1). The four different species of snail hosts could indicate a relatively low host specificity (Fig. 3), but it must be emphasized that species identification was not achieved for these parasites in all cases. If an improved differential diagnostic technique for the separation of species could be applied, we may achieve a higher species resolution which could reveal a higher specificity. This group of parasites was also found in Ringen sø with only one snail species present. This is an artificial lake constructed 40 years ago and harbors relatively few faunistic elements. It is therefore noteworthy that this group of parasites, even with low specificity, will survive in a young ecosystem as this group of trematodes can utilize molluscs as both the first and second intermediate host, transforming into metacercariae in the snails in which they were formed .
The snail R. balthica harbored the most species-rich trematode faunas, with the trematode family Diplostomidae contributing the most. Members of this family group utilize freshwater fish as the second intermediate host and fish-eating birds as the definitive host (Additional file 2: Table S1). BrommeLille sø is a well-known lake hosting a rich teleost fauna. Correspondingly, the parasite fauna showed the strongest relationship with freshwater fish, emphasizing that species richness of trematodes in a snail population relies on the ecology of the freshwater habitat . BrommeLille sø is a shallow lake, with a water depth of 1.5 m , thereby providing optimal conditions for an efficient transmission between birds, freshwater fish and snails. This was also reflected by the positive relationship between bird species richness and trematode richness in an estuarine snail in a coastal wetland . In our investigation of another lake, Tissø, a well-known bird habitat in Denmark, trematodes utilizing larval insects as the second intermediate host and piscivorous birds as the definitive host were dominant (Fig. 5). Also Furesø, a large waterfowl habitat, is characterized by the snail parasites using birds as final hosts. Tjele langsø harbors parasites that reflect a dominant representation of freshwater fish, and in Lyngby sø, the occurrence of tadpoles and frogs was reflected by the dominant representation of specific trematodes. It has been argued that parasite diversity is positively related to host diversity , and based on theoretical predictions, host and trematode species abundance has been found to be correlated . In our study, Pearson’s correlation analysis between these factors revealed a non-significant positive correlation (r = 0.4245, P = 0.1012). It has also been suggested that a high diversity of snail hosts in itself is not the driving factor of trematode species richness; rather, it is the particular composition of snail species present, the abundance of those species and the degree of host specialization among digeneans which determine the trematode species richness . In this context, it is worthwhile to note that we confirmed that L. stagnalis and R. balthica can be regarded as ‘key host species’ hosting the main proportion of the parasite species present . In our study, the prevalence of infection in August and September was relatively high, compared to the spring time, which indicates that temperature could promote the transmission of parasites and increase their local abundance . A relatively high snail and hosts species abundance was found in Furesø, Bagsværd sø and Esrum sø, whereas Ringen sø, Tissø, BrommeLille sø and Tjele langsø contain few species. Species richness may in that case be influenced by short distances between the ponds/lakes, resulting in additional dispersal modes, which in turn results in high immigration rates and/or low extinction rates .
In this work, we have proposed a method for supplementing biological surveys of Danish freshwater systems. By investigating the occurrence of snail-borne trematodes released as cercariae in 21 freshwater lakes in Denmark, we show the presence of additional hosts in the habitats. Using this approach, we found that trematodes released from the snail host as cercariae may act as bioindicators that elucidate the ecology of the investigated freshwater habitat. We confirmed that L. stagnalis and R. balthica can be regarded as ‘key host species’ hosting the main proportion of the parasite species present. The presence of the parasites in a certain freshwater system indicates the existence of other host organisms that act as secondary and final hosts for the parasites. The presence of these hosts may have implications for fish health in the ecosystem due to the pathogenicity of the trematodes, such as eye flukes, heart flukes and blood flukes, that parasitize fish. We also recovered three zoonotic trematode species from the snail species L. stagnalis, R. balthica and S. palustris, in four different lakes. This finding emphasizes that monitoring of trematode infections in snails may be applied as a tool for surveying occurrence of zoonotic parasites with medical importance. The seasonal occurrence of infection suggests that the climate may have a profound influence on the ecosystems. In this paper, we propose a framework for studying trematode–host interactions and association patterns. We therefore advocate for linking molecular data (for precise and easy diagnostics) to an updated ecological situation of any freshwater habitat. Such a tool could be used as a novel and effective prevention system by farms rearing aquatic animals.
Availability of data and materials
The data from the study are available upon publication. All data supporting the conclusion within this article are included within the article. Representative sequences were submitted to the GenBank database under the accession numbers MW000945-MW001148. All raw sampling data are available in Table 2 and Additional file 2: Table S1.
Bray R, Gibson D, Jones A. Keys to the Trematoda, volume 3. Wallingford: CABI International; 2008.
Morley NJ. The role of the invasive snail Potamopyrgus antipodarum in the transmission of trematode parasites in Europe and its implications for ecotoxicological studies. Aquat Sci. 2008;70:107–14.
Roberts LS, Janovy J. Foundations of parasitology. 8th ed. Boston: McGraw-Hill Higher Education; 2009.
de León GP-P, Hernández-Mena D. Testing the higher-level phylogenetic classification of Digenea (Platyhelminthes, Trematoda) based on nuclear rDNA sequences before entering the age of the ‘next-generation’ Tree of Life. J Helminthol. 2019;93:260–76.
Cribb T, Bray R, Littlewood DTJ. The nature and evolution of the association among digeneans, molluscs and fishes. Int J Parasitol. 2001;31:997–1011.
Zhytova E, Romanchuk L, Guralska S, Andreieva OY, Shvets M. Circulation pathways of trematodes of freshwater gastropod mollusks in forest biocenoses of the Ukrainian Polissia. Vestn Zool. 2019;53:13–22.
Shoop WL. Trematode transmission patterns. J Parasitol. 1988;74:46–59.
Jorge F, Poulin R. Poor geographical match between the distributions of host diversity and parasite discovery effort. Proc Biol Sci. 2018;285:20180072.
Paull SH, Johnson PT. High temperature enhances host pathology in a snail–trematode system: possible consequences of climate change for the emergence of disease. Freshw Biol. 2011;56:767–78.
Hechinger RF, Lafferty KD, Kuris AM. Trematodes indicate animal biodiversity in the chilean intertidal and Lake Tanganyika. J Parasitol. 2008;94:966–8.
Hechinger RF, Lafferty KD, Huspeni TC, Brooks AJ, Kuris AM. Can parasites be indicators of free-living diversity? Relationships between species richness and the abundance of larval trematodes and of local benthos and fishes. Oecologia. 2007;151:82–92.
De Leon GP-P. The diversity of digeneans (Platyhelminthes: Cercomeria: Trematoda) in vertebrates in Mexico. Comp Parasitol. 2001;68:1–8.
Faltýnková A, Sures B, Kostadinova A. Biodiversity of trematodes in their intermediate mollusc and fish hosts in the freshwater ecosystems of Europe. Syst Parasitol. 2016;93:283–93.
McNamara M, Miller T, Cribb T. Evidence for extensive cryptic speciation in trematodes of butterflyfishes (Chaetodontidae) of the tropical Indo-West Pacific. Int J Parasitol. 2014;44:37–48.
Wesenberg-Lund C. Contributions to the development of the Trematoda Digenea. Part II. The biology of the freshwater Cercariae in Danish fresh-waters. Memoires de 1’Academic Royale des Sciences et des Lettres de Danemark. Section des Sciences. 1934;5:1–223.
Selbach C, Soldánová M, Georgieva S, Kostadinova A, Sures BJP. vectors, Integrative taxonomic approach to the cryptic diversity of Diplostomum spp in lymnaeid snails from Europe with a focus on the ‘Diplostomum mergi’ species complex. Parasites Vectors. 2015;8:300.
Glöer P. The freshwater gastropods of the West-Palaearctis. Hetlingen: Biodiversity Research Lab; 2019.
Buchmann K. An introduction to fish parasitological methods: classical and molecular techniques. Biofolia; 2007.
Schell S. Key to Cercariae. Handbook of trematodes of North America north of Mexico. Moscow: University of Idaho Press; 1985.
Frandsen F, Christensen N. Introductory guide to the identification of cercariae from African freshwater snails with special reference to cercariae of trematode species of medical and veterinary importance Taxonomic key. Acta Trop. 1984;41:181–202.
Combes C, Imbert-Establet D. Infectivity in rodents of Schistosoma mansoni cercariae of human and murine origin. J Helminthol. 1980;54:167–71.
Sewell RS. Cercariae indicae. Indian J Med Res. 1922;10:13–70.
Curran JE, Johnson MP, Dyer TD, Göring HH, Kent JW, Charlesworth JC, et al. Genetic determinants of mitochondrial content. Hum Mol Genet. 2007;16:1504–14.
Tkach V, Pawlowski JW. A new method of DNA extraction from the ethanol-fixed parasitic worms. Acta Parasitol. 1999;44:147–8.
Galazzo DE, Dayanandan S, Marcogliese DJ, McLaughlin JD. Molecular systematics of some North American species of Diplostomum (Digenea) based on rDNA-sequence data and comparisons with European congeners. Can J Zool. 2002;80:2207–17.
Christiansen AØ, Olsen A, Buchmann K, Kania PW, Nejsum P, Vennervald BJ. Molecular diversity of avian schistosomes in Danish freshwater snails. Parasitol Res. 2016;115:1027–37.
Kumar S, Stecher G, Li M, Knyaz C, Tamura K. MEGA X: molecular evolutionary genetics analysis across computing platforms. Mol Biol Evol. 2018;35:1547–9.
Outa JO, Sattmann H, Köhsler M, Walochnik J, Jirsa F. Diversity of digenean trematode larvae in snails from Lake Victoria, Kenya: first reports and bioindicative aspects. Acta Trop. 2020;206:105437.
Selbach C, Soldánová M, Feld CK, Kostadinova A, Sures B. Hidden parasite diversity in a European freshwater system. Sci Rep. 2020;10:2694.
Brönmark C. Freshwater snail diversity: effects of pond area, habitat heterogeneity and isolation. Oecologia. 1985;67:127–31.
Schock NT, Reisinger AJ, Reisinger LS, Cooper MJ, Cibrowski JJ, Gehring TM, et al. Relationships between the distribution and abundance of the invasive faucet snail (Bithynia tentaculata) and environmental factors in Laurentian Great Lakes coastal wetlands. Biol Invasions. 2019;21:2613–28.
Faltýnková A. Larval trematodes (Digenea) in molluscs from small water bodies near Šeské Budšjovice Czech Republic. Acta Parasitol. 2005;50:49–55.
Soldánová M, Selbach C, Sures B, Kostadinova A, Pérez-del-Olmo A. Larval trematode communities in Radix auricularia and Lymnaea stagnalis in a reservoir system of the Ruhr River. Parasites Vectors. 2010;3:56.
Gordy MA, Kish L, Tarrabain M, Hanington PC. A comprehensive survey of larval digenean trematodes and their snail hosts in central Alberta Canada. Parasitol Res. 2016;115:3867–80.
Soldánová M, Georgieva S, Roháčová J, Knudsen R, Kuhn JA, Henriksen EH, et al. Molecular analyses reveal high species diversity of trematodes in a sub-Arctic lake. Int J Parasit. 2017;47:327–45.
Faltýnková A, Haas W. Larval trematodes in freshwater molluscs from the Elbe to Danube rivers (Southeast Germany): before and today. Parasitol Res. 2006;99:572–82.
Sitko J, Rząd I. Diplostomum and Ornithodiplostomum scardinii (Diplostomidae, Digenea) species from naturally infected birds (Anatinae) in the Czech Republic and in Poland: morphological, morphometric and ecological features. Helminthologia. 2014;51:215–24.
Zajíček D. Cercariae and further developmental stages of flukes in molluscs of some south Bohemian fishponds. Ceskoslovenska Parasitol. 1963;10:187–205.
Morley NJ, Lewis JW. Anthropogenic pressure on a molluscan–trematode community over a long-term period in the Basingstoke Canal, UK, and its implications for ecosystem health. EcoHealth. 2006;3:269–80.
Karvonen A, Savolainen M, Seppälä O, Valtonen ET. Dynamics of Diplostomum spathaceum infection in snail hosts at a fish farm. Parasitol Res. 2006;99:341–5.
Niewiadomska K, Kiseliene V. Diplostomum cercariae [Digenea] in snails from Lithuania. II. Survey of species. Acta Parasitol. 1994;39:179–86.
Wisniewski W. Characterization of the parasitofauna of an eutrophic lake. Acta Parasitol Pol. 1958;6:1–61.
Ginetsinskaya T. Cercarial fauna of molluscs from the Rybinsk water reservoir. I. Systematic survey of cercariae. In: Polyanski Y I, editor. Ekologicheskaya parazitologiya; 1959. p. 96–149.
Cichy A, Faltynkova A, Zbikowska E. Cercariae (Trematoda, Digenea) in European freshwater snails-a checklist of records from over one hundred years. Folia Malacol. 2011;19:165–89.
Combes C. World atlas of Cercariae. Wallingford: CABI International. 1980.
Faltýnková A, Našincová V, Kablásková L. Larval trematodes (Digenea) of the great pond snail, Lymnaea stagnalis (L.),(Gastropoda, Pulmonata) in Central Europe: a survey of species and key to their identification. Parasite. 2007;14:39–51.
Żbikowska E. Digenea species in chosen populations of freshwater snails in northern and central part of Poland. Wiad Parazytol. 2007;53:301–8.
Zdarska Z. Larval trematodes of freshwater snails in Czechoslovakia. Cesk Parazitol. 1963;10:207–62.
Väyrynen T, Siddall R, Valtonen ET, Taskinen J. Patterns of trematode parasitism in lymnaeid snails from northern and central Finland. Ann Zool Fennici. 2000;37(3):189–99.
Blair D. Observations and experiments on some larval trematodes of freshwater snails and fish from southern Iceland. J Helminthol. 1973;47:409–14.
Chernogorenko M. Trematode fauna of mollusks in the Kremenchug Reservoir. Hydrobiol J. 1977;13:87–94.
Dönges J, Götzelmann M. Isthmiophora melis: Experimental reinfection of Lymnaea stagnalis by implantation of miracidia after implantation of rediae. Exp Parasitol. 1977;42:318–21.
Motzel W, Haas W. Studies on the attachment response of Isthmiophora melis cercariae (Trematoda: Echinostomatidae). Z Parasitenkd. 1985;71:519–26.
Radev V, Kanev I, Khrusanov D, Fried B. Reexamination of the life cycle of Isthmiophora melis (Trematoda: Echinostomatidae) on material from southeast Europe. Parazitologiia. 2009;43:445.
Yu S-H, Mott KE, Organization WH: Epidemiology and morbidity of food-borne intestinal trematode infections. Geneva: World Health Organization. 1994. https://apps.who.int/iris/bitstream/handle/10665/61103/WHO_SCHISTO_94.108.pdf. Accessed 9 November 2020.
Chai J-Y: Human intestinal flukes: From discovery to treatment and control. Berlin Heidelberg: Springer Nature; 2019.
Loy C, Haas W. Prevalence of cercariae from Lymnaea stagnalis snails in a pond system in Southern Germany. Parasitol Res. 2001;87:878–82.
Ferte H, Depaquit JM, Carre S, Villena I, Leger N. Presence of Trichobilharzia szidati in Lymnaea stagnalis and T. franki in Radix auricularia in northeastern France: molecular evidence. Parasitol Res. 2005;95:150–4.
Żbikowska E. Do larvae of Trichobilharzia szidati and Echinostoma revolutum generate behavioral fever in Lymnaea stagnalis individuals? Parasitol Res. 2005;97:68–72.
Soldanova M, Selbach C, Sures B. The early worm catches the bird? Productivity and patterns of Trichobilharzia szidati cercarial emission from Lymnaea stagnalis. PLoS ONE. 2016;11:e0149678.
Rudolfová J, Hampl V, Bayssade-Dufour C, Lockyer A, Littlewood D, Horák P. Validity reassessment of Trichobilharzia species using Lymnaea stagnalis as the intermediate host. Parasitol Res. 2005;95:79–89.
Jouet D, Ferté H, Depaquit J, Rudolfová J, Latour P, Zanella D, et al. Trichobilharzia spp. in natural conditions in Annecy Lake France. Parasitol Res. 2008;103:51.
Haas W, Haeberlein S. Penetration of cercariae into the living human skin: Schistosoma mansoni vs Trichobilharzia szidati. Parasitol Res. 2009;105:1061–6.
Haas W, van de Roemer A. Invasion of the vertebrate skin by cercariae of Trichobilharzia ocellata: penetration processes and stimulating host signals. Parasitol Res. 1998;84:787–95.
Wikgren B-J. Studies on Finnish larval flukes with a list of known Finnish adult flukes (Trematoda: Malacocotylea). Acta Zool Fennica. 1956;91:1–106.
Gao J-F, Gao Y, Qiu J-H, Chang Q-C, Zhang Y, Fang M, et al. De novo assembly and functional annotations of the transcriptome of Metorchis orientalis (trematoda: Opisthorchiidae). Exp Parasitol. 2018;184:90–6.
Li C, Wu H, Sun E, Zhu Y. Investigation on the endemic characteristics of Metorchisorientalis in Huainan area China. Nutr Hosp. 2017;34:675–9.
Bowman DD, Hendrix CM, Lindsay DS, Barr SC. Feline clinical parasitology. 1st ed. Hoboken: Wiley; 2008.
Murell K, Pozio E. The liver flukes: Clonorchis sinensis, Opisthorchis spp, and Metorchis spp. Global Water Pathogen Project. 2017. http://www.waterpathogens.org. Accessed 9 November 2020.
Na L, Gao J-F, Liu G-H, Fu X, Su X, Yue D-M, et al. The complete mitochondrial genome of Metorchis orientalis (Trematoda: Opisthorchiidae): Comparison with other closely related species and phylogenetic implications. Infect Genet Evol. 2016;39:45–50.
Jinxuang L, Youzhu C, Yousong L. The discovery of natural infection of human with Metorchis orientalis and the investigation of its focus. Chin J Zoonoses. 2001;17:19–20.
Tao C. Notes on the study of life-cycle of Metorchis orientalis and M.taiwanensis. Chin Rev Trop Med. 1948;1:9–14.
Zikmundová J, Georgieva S, Faltýnková A, Soldánová M, Kostadinova A. Species diversity of Plagiorchis Lühe, 1899 (Digenea: Plagiorchiidae) in lymnaeid snails from freshwater ecosystems in central Europe revealed by molecules and morphology. Syst Parasitol. 2014;88:37–54.
Styczynska-Jurewicz E. The life cycle of Plagiorchis elegans (Rud, 1802) and the revision of the genus Plagiorchis Lühe, 1889. Acta Parasitol. 1962;10:419–45.
Żbikowska E, Nowak A. One hundred years of research on the natural infection of freshwater snails by trematode larvae in Europe. Parasitol Res. 2009;105:301.
Høy T. Morphometry of 150 Danish Lakes with some considerations on cartometry. Geogr Tidsskr. 1988;88:44–54.
Hechinger RF, Lafferty KD. Host diversity begets parasite diversity: bird final hosts and trematodes in snail intermediate hosts. Proc R Soc B. 2005;272:1059–66.
Thieltges DW, Hof C, Dehling DM, Brändle M, Brandl R, Poulin R. Host diversity and latitude drive trematode diversity patterns in the European freshwater fauna. Glob Ecol Biogeogr. 2011;20:675–82.
Johnson PT, Preston DL, Hoverman JT, LaFonte BE. Host and parasite diversity jointly control disease risk in complex communities. Proc Natl Acad Sci USA. 2013;110:16916–21.
Kamiya T, O’Dwyer K, Nakagawa S, Poulin R. Host diversity drives parasite diversity: Meta-analytical insights into patterns and causal mechanisms. Ecography. 2014;37:689–97.
Poulin R. Global warming and temperature-mediated increases in cercarial emergence in trematode parasites. Parasitology. 2006;132:143–51.
Sannia A, James B. The digenea in marine molluscs from Eyjafjördur North Iceland. Ophelia. 1977;16:97–109.
Manga-Gonzalez Y, González-Lanza C, Kanev I. Lymnaea truncatula, intermediate host of some Plagiorchiidae and Notocotylidae species in Leon NW Spain. J Helminthol. 1994;68:135–41.
Jezewski W. Occurrence of Digenea (Trematoda) in two viviparus species from lakes. Helminthologia. 2004;41:147–50.
Faltýnková A, Niewiadomska K, Santos M, Valtonen E. Furcocercous cercariae (Trematoda) from freshwater snails in Central Finland. Acta Parasitol. 2007;52:310–7.
Schwelm J, Kudlai O, Smit N, Selbach C, Sures B. High parasite diversity in a neglected host: larval trematodes of Bithynia tentaculata in Central Europe. J Helminthol. 2020;94:1–23.
Masahaii N. New records of trematode parasites (Digenea) in the banded frog (Rana camerani) and marsh frog (Rana ridibunda)(Anura: Ranidae), from Southwest of Iran. Iran Fish Sci J. 1999;1:41–8 (in Iranian).
Chikhlyaev IV, Ruchin AB, Fayzulin AI. An overview of the trematodes fauna of pool frog Pelophylax lessonae (Camerano, 1882) in the Volga Basin, Russia: 1 Adult stages. Nusantara Biosci. 2018;10:256–62.
Esteban J, Toledo R, Sánchez L, Muñoz-Antoli C. Life-cycle of Euparyphium albuferensis n. sp. (Trematoda: Echinostomatidae) from rats in Spain. Syst Parasitol. 1997;38:211–9.
Enabulele EE, Lawton SP, Walker AJ, Kirk RS. Molecular and morphological characterization of the cercariae of Lecithodendrium linstowi (Dollfus, 1931), a trematode of bats, and incrimination of the first intermediate snail host Radix balthica. Parasitology. 2018;145:307–12.
Murai E, Meszaros F, Sey O. On parasitic helminths of mammals living in the environs of Lake Balaton. Parasit Hung. 1992;25:23–36.
Moema EB, King PH, Rakgole JN. Phylogenetic studies of larval digenean trematodes from freshwater snails and fish species in the proximity of Tshwane metropolitan South Africa. Onderstepoort J Vet Res. 2019;86:1–7.
Besprozvannykh V, Rozhkovan K, Ermolenko A, Izrailskaya A. Diplodiscus mehrai Pande, 1937 and D. japonicus (Yamaguti, 1936): morphology of developmental stages and molecular data. Helminthologia. 2018;55:60–9.
Nugaraitė D, Mažeika V, Paulauskas A. Molecular and morphological characterization of Isthmiophora melis (Schrank, 1788) Luhe, 1909 (Digenea: Echinostomatidae) from American mink (Neovison vison) and European polecat (Mustela putorius) in Lithuania. Helminthologia. 2017;54:97–104.
Aksenova OV, Bespalaya YV, Bolotov IN, Kondakov AV, Sokolova SE. First molecular identification of Australapatemon burti (Miller, 1923) (Trematoda: Digenea: Strigeidae) from an intermediate host Radix labiata (Rossmaessler) (Gastropoda: Lymnaeidae) in Europe. Zootaxa. 2016;4132:588–90.
Qiu J-H, Zhang Y, Zhang X-X, Gao Y, Li Q, Chang Q-C, et al. Metacercaria infection status of fishborne zoonotic trematodes, except for Clonorchis sinensis in fish from the Heilongjiang Province China. Foodborne Pathog Dis. 2017;14:440–6.
Borkovcová M, Kopřiva JJPR. Parasitic helminths of reptiles (Reptilia) in south Moravia (Czech Republic). Parasitol Res. 2005;95:77–8.
We thank the Chinese Scholarship Council (No. CSC201706910086) for the financial support by a stipend to YJD.
We acknowledge financial support from the Danish Veterinary and Food Administration (FVST) as part of the agreement of commissioned work between the Danish Ministry of Food, Agriculture and Fisheries and the University of Copenhagen.
Ethics approval and consent to participate
All animal procedures were performed observing the general ethical guidelines of the University of Copenhagen. All snail samples were euthanized after the shedding procedures.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Phylogenetic analysis of ITS sequences. Own sequences are shown in red.
Overall prevalence of the trematode species (% of infected snails in the pooled samples from March to September in 2019) infecting the eight snail species in the 21 lakes. Lakes where prevalence was zero are indicated by a minus (−). Intermediate and definitive hosts are provided based on published data from literature [6, 87,88,89,90,91,92,93,94,95,96,97]. *small sample size.
About this article
Cite this article
Duan, Y., Al-Jubury, A., Kania, P.W. et al. Trematode diversity reflecting the community structure of Danish freshwater systems: molecular clues. Parasites Vectors 14, 43 (2021). https://doi.org/10.1186/s13071-020-04536-x
- Molecular identification
- Freshwater lakes