Geographic variation in composition of metazoan parasite infracommunities in Galaxias maculatus Jenyns 1842 (Osmeriformes: Galaxiidae) in southern Chile (38-47° S)

Abstract

Galaxias maculatus is an abundant freshwater fish species in Chilean continental waters where it plays important ecological functions, yet few parasitological records of this species exist in Chile and all of them cover a very limited geographic range. The objective of this study was to assess large scale geographic variation in composition of parasite infracommunities of Galaxias maculatus. Specifically, parasite infracommunities of this species were compared among 11 locations across 9 degrees of latitude and 3 ecosystem types (lake, river and estuary). Most taxa found had been previously reported in Chile and Argentina. However, this is the first report for Tylodelphys sp. in this host in Chile. Furthermore, the cranial parasite Tylodelphys sp. had the highest overall prevalence and abundance compared to other parasite species. Despite the fact that the abundance of Tylodelphys sp. was not significantly correlated with Fulton’s condition factor of fish, infected fish seem to have a better body condition compared to uninfected ones. The most important source of variation in composition of infracommunities was the sampling location. Furthermore, fish from lakes have a different composition of parasite infracommunities mainly due to higher abundances of Tylodelphys sp.

Background

Puye Galaxias maculatus Jenyns 1842 (Osmeriformes: Galaxiidae) is one of the most widely distributed freshwater fish species in the Southern Hemisphere [1]. In Chile, Galaxias maculatus is widespread and often abundant in rivers of central Chile and Patagonia, where it plays an important ecological role [2, 3]. There are several parasitological records for this species in Argentina, Australia and New Zealand [4,5,6,7], but just five in Chile [8,9,10,11,12] and these covered very limited geographic range (one to three sampling locations in each study). The aim of this study was to assess geographic variation in composition of infracommunities of metazoan parasites of G. maculatus in Chile (excluding Myxozoa). This was done by comparison of parasite communities among 11 locations across 9 degrees of latitude and 3 ecosystem types (lake, river and estuary, Fig. 1).

Fig. 1

The study area indicating fish sampling locations and ecosystems (rivers, lakes and estuaries)

Methods

One-hundred and sixty-six individuals of G. maculatus were collected during 2013 using a beach seine (5 m long, 1.5 m high and 10 mm stretched mesh size). The beach seine was hauled over a distance of 10–20 m in the following river basins in southern Chile (Fig. 1: Imperial (river, 38°44′52.42″S, 73°24′27.19”O), Valdivia (lake, river and estuary, 39°51′55.33″S, 73°21′21.89”O), Maullín (lake and river, 41°12′38.98″S, 73°2′3.86”O), Puelo (estuary, 41°39′08.75″S, 72°18′16.78″W), Aysén (lake and river, 45°24′17.87″S, 72°47′41.62″W) and Baker (lake and estuary, 47°47′25.62″S, 73°29′57.42″W). Sample sizes ranged from 6 to 20 host fish per sampling location. Small sample size in some locations was due to low host abundances (catch per unit effort). Each fish specimen was measured (total length, TL; in mm) and weighed (total weight, W; in g), then preserved in 70% ethanol and stored individually. Subsequently, each specimen was thoroughly scraped on the body surface, fins and gills to look for ectoparasites. For endoparasites, we examined eyes, brain, heart, intestine, stomach, liver as well as the cranial and body cavity under a stereomicroscope. Finally, longitudinal cuts were made in the musculature to verify presence of metacercariae.

The variation in TL of G. maculatus among locations and ecosystem types was assessed using general linear models. Subsequently, the size effect was analyzed using partial eta-squared statistic. These analyses were performed in SAS/STAT/PROC GLM software, Version 9.3 (2011). A posteriori tests were run using GT-2 because of unequal sample sizes. The variation in the composition of parasite infracommunities was analysed using permutational multivariate analysis of variance (PERMANOVA) [13] based on fourth-root transformed data and Bray Curtis + 1 similarity index. In this analysis, the ecosystem type, sampling locations (nested within ecosystem type) were considered as factors, with total body length as covariate. A Principal Coordinates Analysis (PCoA) was used to visually display the distribution of parasite infracommunities among sampling locations and ecosystems using the vectors whose Spearman correlation coefficients were higher than 0.5. Further insight on this ordination was carried out with a One-way ANOVA on the scores of the first PCoA axis with the type of ecosystem as a factor. Furthermore, distance-based test for homogeneity of multivariate dispersions (PERMDISP) was used to compare composition of parasite infracommunities among ecosystems. These analyses were performed in PRIMER statistical program Version 7 [14]. In addition, a Spearman correlation coefficient was used to assess the association between the distances of sampling locations (in km) with the similarity in composition of parasite infracommunities. These analyses were performed in XLSTAT® software. Finally, the eventual effect of Tylodelphys sp. on Fulton’s condition factor of fishes (a quantitative indicator of host’s well-being) was assessed in two ways: by examining the statistical significance of the Spearman correlation coefficient between the abundance of parasites and the condition factor, and by performing a Mann-Whitney test to compare condition factor of fishes with and without Tylodelphys sp.

Results

Total length of 166 G. maculatus specimens showed medium size effect and significant variations among sampling locations (overall mean ± S.D. = 3.62 ± 0.56 cm, F _{10, 155} = 15.04, R² = 0.49, P < 0.0001, size effect partial eta-squared 90% confidence limits = 0.37–0.54). Specifically, the mean TL of specimens of the small sample from Puelo Estuary (n = 6) differed from other locations. Furthermore, fish collected from the Baker River basin (Lbak) were the smallest (Table 1). TL size effect showed significant though small variations among ecosystems (F _{2, 163} = 4.13, P = 0.018, R² = 0.049, size effect partial eta-squared 90% confidence limits = 0.00–0.10). More specifically, fish from rivers were larger than those from lakes (GT-2 Test, P < 0.05).

Table 1 Number of fish examined, total length (mm), weight (g) standard deviation, and number of parasites per sampling location and taxon. Total number of taxa per location (last row) and total number of parasites per location as well as the range of locations in which each taxon occur (right column)

Sixty-four percent of the specimens harbored at least one parasite. A total of 12 metazoan taxa of parasites was recorded. Most taxa had low prevalence and abundance, except for a diplostomulum type larva found free (unencysted) in the meningeal space (overall prevalence = 49%). This parasite accounted for 89% of all 1214 parasitic individuals collected, and was the most prevalent and widespread parasite across locations (Table 1). This parasite was determined as Tylodelphys sp., in the absence of molecular confirmation. However, it is morphometrically similar to Tylodelphys barilochensis Quaggiotto and Valverde 1992 (Digenea: Diplostomatidae) (14), based on its body length and width (in microns) (body length: average ± S.D. = 456.3 ± 126.3, maximum body width: 127.3 ± 15.0, n = 10). There was no significant correlation between the Fulton’s condition factor and the abundance of Tylodelphys sp. (r_s = − 0.196, P = 0.081). However, condition factor of non-infected fish was lower than that of fishes infected with Tylodelphys sp. (Mann-Whitney U one-tailed test, U = 2871, P = 0.032).

Multivariate dispersion of parasite infracommunities did not differ among ecosystems (Pseudo-F _{2, 163} = 2.09, P (perm) = 0.188). The largest part of the variation was associated with sampling locations (GT-2 test, P < 0.05, Table 2). Two first axes of the PCoA analysis of the composition of parasite infracommunities accounted for 81.5% of the variation (Fig. 2). The infracommunities of the Maullín River separated from all other locations, mostly due to the absence of Tylodelphys sp. in this basin. Two other parasite taxa were relatively prevalent: the branchiuran Argulus sp. (mostly found in Maullín River), and the larval nematode Contracaecum sp. The one-way ANOVA on the scores of the first PCoA axis with the type of ecosystem as a factor revealed significant differences in composition of infracommunities in fish coming from lakes (F _{7, 158} = 9.438, P < 0.0001) due to higher abundance of Tylodelphys sp. in lakes than in rivers or estuaries (Table 2). The number of taxa per sampling location ranged from 1 to 8, and was not correlated with the number of fish examined (Spearman correlation, r_s = 0.01, n = 11, P = 0.96, Table 1). Finally, the distance between sampling locations was not correlated with the similarity of parasite infracommunities (Spearman correlation, r_s = 0.007, n = 55 pairs of distances between 11 sampling locations, P = 0.55).

Table 2 Summary of PERMANOVA of the composition of metazoan parasite infracommunities in 166 G. maculatus according to their total length (TL), type of ecosystem (Ec) and sampling location (Lo)

Fig. 2

Principal Coordinates Analysis (PCo) of the composition of parasite infracommunities in G. maculatus according to sampling locations and ecosystems. Vectors indicate species that had Spearman’s rank correlations > 0.5

Discussion

This report differs from previous ones in Chile, mainly due the broader geographical range and variety of ecosystems sampled. Previously, the only study on parasites of G. maculatus that considered more than two sampling locations in Chile dealt with Stephanostomum sp. metacercariae only [8]. Another difference of this report is that it is the first carried out at the infracommunity level. Major taxonomic differences found with previous studies was the lack of reports on Tylodelphys sp. for this host in Chile. In neighbouring Argentina, Tylodelphys spp. have been widely reported from several fish species [5, 6, 15,16,17,18,19]. This is also a frequent worm reported for many fish species and places elsewhere (see [20, 21]). Austrodiplostomum mordax Szidat and Nani, 1951 and Tylodelphys destructor Szidat and Nani 1951 were reported in another freshwater fish, Basilichthys australis Eigenmann 1927, from Lake Riñihue in the South of Chile [22, 23]. It is not surprising that fish infected with Tylodelphys sp. show a higher condition factor compared to uninfected fish, as previously reported for other endoparasites in freshwater fish [24]. Furthermore, this better condition is expected to be associated with strategies of transmission to definitive hosts [24].

Differences found in parasite infracommunity composition among sampled locations may be associated with the availability and vagility of intermediate and definitive hosts. For example, the absence of records of Tylodelphys sp. in G. maculatus in Chile may be due to the low abundance of its first intermediate host Chilina dombeyana (Bruguiére, 1789) in some sampling locations. In addition, spatial proximity of sampling locations was not a good predictor of similarity of parasite infracommunities, what has also been shown in other freshwater fish species [25]. With a mixture of larval and adult stages of parasites, there is likely a wide variety of hosts involved in parasite transmission, many of which are invertebrates that have limited dispersal ability, potentially resulting in clumped distributions. Although some results indicated that type of ecosystem was not a relevant factor to account for variations in composition of parasite infracommunities, the ANOVA on the scores of the first PCoA axis revealed that fish from lakes were different to those from rivers and estuaries, mainly due to higher abundances of Tylodelphys sp.. This pattern may be an effect of differential influence of autogenic and allogenic parasite species on patterns of parasite distribution and community composition [26, 27].