Skip to main content
Download PDF

The metazoan parasite communities of the shoal flounder (Syacium gunteri) as bioindicators of chemical contamination in the southern Gulf of Mexico

Parasites & Vectors volume 7, Article number: 541 (2014) Cite this article

Abstract

Background

Because agriculture and offshore oil extraction are significant economic activities in the southern Gulf of Mexico, high concentrations of nutrients and hydrocarbons are expected. As parasite communities are sensitive to environmental impacts, these contaminants should have an effect on metrics such as species richness, relative abundance and similarity. Consequently, these community metrics can be used as indicators of aquatic environmental health. Our objectives were to describe the parasite communities of the shoal flounder Syacium gunteri and to determine potential thresholds above which environmental contaminants become major controlling factors of parasite community metrics.

Methods

The study area included 33 sampling sites in the southern Gulf of Mexico, where benthic sediments, water and shoal flounder individuals were collected. Data on ecto- and endo-parasites from flounder and nutrients, contaminants and physicochemical variables from the water and sediments were obtained. The statistical associations of the parasite community metrics at the component and infracommunity levels and the environmental data were analysed using redundancy analysis (RDA).

Results

Overall, 203 shoal flounder were examined for parasites, recovering 13 metazoan parasite species, and 48 physicochemical (e.g. temperature, nutrients) and contaminant (e.g. hydrocarbons, heavy metals) variables were obtained. The larval stages of the cestode Oncomegas wageneri and the nematodes Pseudoterranova decipiens and Hysterothylacium sp. were numerically dominant at the component and infracommunity levels. The parasite community metrics had significant negative statistical associations with both nitrate and total PAHs. With the exception of these two chemicals, which exceeded the threshold effect levels (TELs), no other environmental variable exceeded the range considered safe for marine organisms.

Conclusions

The community metrics chosen generally had robust statistically significant associations with both physicochemical and contaminant variables, which supports the ecological relevance of these parameters as indicators of aquatic environmental health. Within the study area, the shoal flounder and their parasites live in a polluted environment with relatively high levels of hydrocarbons and nitrate. Regarding nitrate, we emphasise that if uncontrolled sewage discharge continues in the southern Gulf of Mexico, hypoxic conditions similar to those caused by the Mississippi river can be expected in the near future.

Background

Marine biologists currently recognise that both marine free-living benthic organisms and their parasites are useful as bioindicators of aquatic environmental health, i.e., as species or communities used to assess the quality of the environment and how it changes over time [1]-[6]. Parasites of aquatic organisms are good bioindicators of environmental impact because their populations and communities are sensitive to environmental insults [5],[7]-[9]. For example, parasite abundance in Finnish lakes that were highly polluted with pulp and paper effluents (PPE) was lower than in lakes with lower PPE levels [10]. After environmental conditions in the lakes improved, parasites returned to levels similar to those in less polluted lakes. At the community level, Huspeni and Lafferty [11] used larval digenean infections of snails as bioindicators of environmental recovery in salt marshes. After six years, the number of infected snails in the impacted zones recovered to at least as high as in the control zones. Recently, two meta-analyses have noted the usefulness of parasite communities as a bioindicator of the recovery of impacted habitats [5],[12]. Vidal-Martínez [5] showed that parasite communities have a positive interaction term with eutrophication, meaning that this environmental impact increases several community metrics such as species richness or relative abundance.

To monitor the effect of the contamination in the southern Gulf of Mexico produced by the oil industry and other environmental insults such as unregulated sewage discharge or pesticide transport by river run-off, studies have been performed on species of parasitic helminths and protozoans of the Mayan catfish Hexanematichthys assimilis[13], the pink shrimp Farfantepenaeus duorarum[14], and the pufferfish Spheroides testudineus[15]. These studies have concluded that chemical contamination influences parasite infection parameters.

During studies in the southern Gulf of Mexico to determine the environmental quality of sediments and water for the Mexican Petroleum company (PEMEX), we obtained data on the helminth parasites infecting the shoal flounder Syacium gunteri. Given that agriculture and offshore petroleum extraction are significant economic activities in the southern Gulf of Mexico [16], high concentrations of nutrients from river run-off, hydrocarbons, and other contaminants such as heavy metals or pesticides would be expected. Because parasite communities are sensitive to environmental impacts, we hypothesise that these contaminants should have a direct effect on parasite community metrics such as species richness, relative abundance or similarity within and among sampling sites. As the parasite communities of S. gunteri have yet to be described for the region, the objectives of this study were twofold: 1) to describe the parasite communities of the shoal flounder at the component and infracommunity levels and 2) to determine whether there are thresholds above which environmental contaminant levels become major controlling factors of parasite community metrics.

Methods

Study area and sampling procedures for sediments

The study area included 33 sampling stations in the southern Gulf of Mexico (Figure 1). Benthic sediments were collected between September 6 and October 8, 2005, at depths between 1.30 and 253.5 m from the oceanographic vessel (OV) Justo Sierra using a 0.25 m2 Hessler Sandia MK-III box corer. Water samples were taken at 5 m depth intervals using 1-gallon amber glass bottles, which were closed under water to avoid contamination with surface mixtures. We obtained 48 physicochemical parameters from water and sediments including oxygen (mg/L), salinity (UPS), pH, and nitrogen concentrations (see Additional file 1 of the online supporting material for a complete list per sampling site). Sediment samples were placed in high-density polythene (HDPE) bags and maintained at 4 °C for transport to CINVESTAV-IPN, Unidad Mérida. Hydrocarbon sampling procedures have been described elsewhere [14],[17],[18]. The physicochemical characteristics, hydrocarbon and metal concentrations of the sediment were determined in the Marine Chemistry and Geochemistry laboratories, respectively, by applying standardised methods [19]-[21]. Principal coordinates of neighbour matrices (PCNM) analyses were used to generate a set of spatial variables from the geographical position of each station [22],[23]. These spatial variables (called PCNM vectors) are spectral decompositions of the spatial relationships among stations corresponding to all the spatial scales that can be perceived from the data [23]. To calculate these variables, we followed the procedure of Santana-Piñeros [24]. Both environmental and biological data were obtained within a very narrow window of time (one month). Therefore, the temporal variability in both types of data is expected to be minimal, and thus, most of the variability in the data should be spatial.

Figure 1
figure 1

Sampling sites in Campeche sound, Gulf of Mexico, where sediment and water samples were obtained. The shoal flounders Syacium gunteri were caught in trawls around each station. Symbols were as follows: sites of group 1; sites of group 2; sites of group 3.

Full size image

Sampling procedures for flatfishes and helminth parasites

Fish were collected using 1 h trawls with 20 m shrimp nets. The trawls lasted 50–60 minutes at 0.6-0.7 knots around each station. The shoal flounder Syacium gunteri were kept in isolated plastic bags in coolers in the cold storage of the vessel and transported to CINVESTAV-IPN Unidad Mérida for parasitological examination. The total length, standard length, maximum height (cm) and weight (g) were recorded for each individual fish. Subsequently, the body surface, cavities and all internal organs were individually examined using a dissection microscope. The metazoan parasites were counted in situ and preserved in 70% alcohol (digeneans) or 4% formalin (nematodes) in labelled vials for subsequent taxonomic identification. Digeneans, cestodes and acanthocephalans were stained using the Mayer paracarmine technique, and the nematodes were cleared using increasing concentrations of glycerine [25]. Voucher specimens were deposited in the National Collections, Universidad Nacional Autónoma de México (CNHE).

Prevalence and mean abundance were calculated for each metazoan parasite species following Bush et al. [26]. Host-specialist helminth species were those previously reported from fish species of the same genus or fish family, whereas generalist species were those previously reported in fish species from different fish families. Species richness and similarity were calculated separately at the component and infracommunity level [27] for total communities (i.e., including every possible infected organ). Autogenic species were those completing their life cycles within aquatic environments. In contrast, allogenic species were those reaching sexual maturity in piscivorous birds or terrestrial mammals [28]. An infracommunity was defined as all of the metazoan parasites infecting an individual fish. At the infracommunity level, the mean ± SD (standard deviation) of the number of species and individual helminths per host examined were calculated. Brillouin’s diversity index [29] was calculated for all infracommunities and expressed as the mean ± SD for each host sample. The Berger-Parker dominance index [29] was calculated at the infracommunity level for the parasites of each individual shoal flounder and presented as the mean ± SD of the index per sampling site. The component community included all of the metazoan parasites infecting a sample of flatfishes of the same species at a specific sampling site. At the component community level, the total number of species and the total number of individual metazoan parasites per sampling site were obtained. The Berger-Parker dominance index [29] was calculated at the component community level for the parasite fauna of the sample of shoal flounder at each sampling site. The qualitative and quantitative faunal similarities were determined using Jaccard’s similarity index and the percentage of similarity [29] for all possible combinations within and among sampling sites at both the component and infracommunity levels. Data normality was determined by using rankit plots [30]. Variables with large deviations from normality were natural log-transformed. One-way ANOVA was used to determine potential differences among the mean numbers of species and individuals per sampling site. Where transformation did not improve normality, the non-parametric Kruskal-Wallis (KW) test was used. The significance of all statistical analyses was established at α <0.05, unless otherwise stated. We focused on community metrics to obtain a more holistic view of the response of parasites rather than concentrating on individual parasite species. Due to the large number of sampling sites, we examined whether these sites would group based on the parasite species composition. Thus, we used cluster analysis with UPGMA (Unweighted Pair Group Method with Arithmetic Mean) as a hierarchical clustering method and Jaccard’s similarity index at the component community level to determine the possible presence of groups of sampling sites that were similar due to their parasite species composition. The cluster analysis was undertaken using MVSP V. 3.22. Then, the tree obtained was then arbitrarily cut at 0.5 of Jaccard’s similarity. Thus, the multivariate analysis (see below) was applied to the original 33 sampling sites and to the groups of sampling sites produced by the cluster analysis.

By using detrended correspondence analysis (DCA), the length of the ordination axes scale (gradient) in standard deviation units (SDU) for the species data was found to be 0.06 SDU. This SDU value indicates that our data only include a limited portion of the environmental range of the shoal flounder, and this analysis assumes linear associations between the environmental and biological variables. Because the SDU value was less than 3, the recommended option for analysis is redundancy analysis (RDA) [31]. RDA is a constrained form of both PCA and multivariate multiple regression and was applied using CANOCO [31],[32] to determine possible associations between environmental variables and the infracommunity and component community metrics. Covariable analysis was used to control for confounding variables such as standard length and weight and to determine the actual influence of independent environmental variables on the infracommunity and component community metrics. Monte Carlo tests were used to determine the significance of the canonical axes for both parasite community metrics and environmental variables. The use of constrained models decreased the percentage of explained variance, although the models were still statistically significant [31].

Results

A total of 203 shoal flounder from 33 sampling sites (Figure 1) were examined for parasites, and 13 metazoan parasite species were recovered. The cluster analysis showed the presence of three groups of sampling sites, called groups 1, 2 and 3 (Figure 2). There were no significant differences in fish standard length (ANOVA ONE-WAY, F[2, 200] = 1.61, p = 0.20) or weight (ANOVA ONE-WAY, F[2, 200] = 3.09, p = 0.05) among the three groups. In contrast, Fulton’s condition factor for the fish from group 1 was significantly higher than for groups 2 and 3 (ANOVA ONE-WAY, F[2, 200] = 4.36, p =0.01, n = 203). The metazoan parasites infecting the shoal flounder were two digeneans (one larval and one adult), five nematodes (four larvae and one adult), four cestodes (all in the larval stage), one adult acanthocephalan and one adult parasitic copepod (Table 1). The prevalence and mean abundance of the metazoan parasites infecting the shoal flounder are summarised in Table 1 for all 33 sampling stations in the overall category and for groups 1 through 3. In all of the groups of sampling sites, the most frequent and abundant species was the larval cestode Oncomegas wageneri, followed by the larval nematodes Pseudoterranova decipiens and Hysterothylacium sp. Overall, there was no significant difference among groups (ANOVA ONE-WAY, F[2, 200] = 3.06, p = 0.05) in the mean abundance of O. wageneri. However, there were significant differences in the mean abundance of P. decipiens among the groups (ANOVA ONE-WAY, F[2, 200] = 4.85, p = 0.008), with Tukey’s test indicating that group 2 had significantly more P. decipiens than did groups 1 and 3 (p <0.05). The remainder of the metazoan parasites in Table 1 had a low prevalence and low mean abundance, which in turn precluded meaningful statistical analysis.

Figure 2
figure 2

Groupings based on the parasite species of Syacium gunteri shared among sampling sites. The dendogram was produced using the unweighted pair-group method (UPGMA) and Jaccard’s similarity index. Note that the symbols for groups 1 (), 2 (), and 3 () are the same as in Figure 1.

Full size image
Table 1 The metazoan parasites of the shoal flounder Syacium gunteri in the southern Gulf of Mexico
Full size table

Component communities

The metazoan parasite fauna of the shoal flounder consisted of 13 species, with groups 1 and 2 containing 12 and 10 species, respectively, and group 3 containing the lowest number of species (4) (Table 2). There were significant differences in the total number of species per sampling site among groups (ANOVA ONE-WAY, F[2, 30] = 3.76, p = 0.03), with Tukey’s test indicating that group 3 had significantly fewer species than groups 1 and 2 (p <0.05). The number of individuals was variable among groups, with group 1 having the lowest number of individuals (Table 2). However, there was no significant difference among groups in the mean number of individuals per sampling site (ANOVA ONE-WAY, F[2, 30] = 1.43, p = 0.25). The overall diversity of the metazoan parasite component communities of S. gunteri is presented in Table 2. There were significant differences among groups in Simpson’s diversity index (ANOVA ONE-WAY, F[2, 30] = 5.50, p = 0.009), with Tukey’s test indicating that group 3 had a significantly lower Simpson’s diversity index than groups 1 and 2 (p <0.05). There were also significant differences among groups in the numerical dominance expressed by the Berger-Parker dominance index (ANOVA ONE-WAY, F[2, 30] = 5.97, p = 0.006), with Tukey’s test indicating that group 3 had a significantly higher numerical dominance than groups 1 and 2 (p <0.05).

Table 2 The component communities of the metazoan parasites of the shoal flounder Syacium gunteri
Full size table

The qualitative similarity of the component community within the groups was variable, with group 3 having the highest similarity (0.81 ± 0.17) within its sampling sites (Table 2). There were significant differences among groups in Jaccard’s index of qualitative similarity (ANOVA ONE-WAY, F[2, 169] = 97.36, p <0.001), with Tukey’s test indicating that group 3 had a significantly higher qualitative similarity than groups 1 and 2 (p <0.05). The quantitative similarity of the component communities within groups was also variable, with group 3 having the highest similarity (53.72 ± 24.32) within its sampling sites (Table 2). There were also significant differences among groups in the percentage of similarity index of the component communities (ANOVA ONE-WAY, F[2, 169] = 20.93, p <0.001), with Tukey’s test indicating that group 3 had significantly higher quantitative similarity than groups 1 and 2 (p <0.05).

Infracommunities

The metazoan parasite infracommunities of the shoal flounder for the overall group and groups 1 through 3 are presented in Table 3, with group 3 having the lowest mean number of species. There were significant differences in the mean number of species per fish among groups (ANOVA ONE-WAY, F[2, 200] = 18.40, p <0.001), with Tukey’s test indicating that group 3 had significantly fewer species than groups 1 and 2 (p <0.05). The mean number of individual parasites per fish among groups was variable. However, no significant difference among groups in the mean number of individuals per fish was detected (ANOVA ONE-WAY, F[2, 200] = 0.17, p = 0.85). The overall values of the Brillouin’s diversity index for the parasite infracommunities of S. gunteri can be seen in Table 3. There were significant differences in the mean Brillouin’s diversity index among groups (ANOVA ONE-WAY, F[2, 200] = 305.69, p <0.001), with Tukey’s test indicating that group 3 had significantly lower values of the diversity index than groups 1 and 2 (p <0.05). There were also significant differences among groups in the numerical dominance expressed by the mean Berger-Parker dominance index per individual fish (Kruskal-Wallis test, H = 17.28, p <0.001, n = 203), with the rank test indicating that group 2 had significantly lower numerical dominance than groups 1 and 3 (p <0.05).

Table 3 The infracommunities of the metazoan parasites of the shoal flounder Syacium gunteri
Full size table

The mean qualitative similarity of the infracommunity within groups was variable, with group 3 having the highest similarity (0.86 ± 0.10) within its sampling sites (Table 3). There were significant differences among groups in Jaccard’s index of qualitative similarity (ANOVA ONE-WAY, F[2, 7194] = 546.56, p <0.001, n = 7197 pair comparisons), with Tukey’s test indicating that group 3 had significantly higher qualitative similarity than groups 1 and 2 (p <0.05). The quantitative similarity of the infracommunity within groups was also variable, with group 3 having the highest similarity (66.89 ± 17.17) within its sampling sites (Table 3). There were significant differences among groups in the quantitative similarity of the infracommunity (ANOVA ONE-WAY, F[2, 7536] = 128.61, p <0.001), with Tukey’s test indicating that group 3 had a significantly higher numerical dominance than groups 1 and 2 (p <0.05).

Qualitative and quantitative similarity among the groups of sampling sites

The qualitative and quantitative similarities between pairs of groups of sampling sites are presented in Table 4. The highest qualitative and quantitative similarities at the component community level were between groups 1 and 3 (qualitative: 0.64 ± 0.21; quantitative: 43.11 ± 25.55). However, there were significant differences among the groups in terms of qualitative similarity (ANOVA ONE-WAY, F[2, 169] = 97.36, p <0.001, n = 172), with Tukey’s test indicating significant differences among the three groups of sampling sites (p <0.05). For the quantitative similarity at the component community level, there were significant differences among the groups (ANOVA ONE-WAY, F[2, 169] = 20.93, p <0.001, n = 172), with Tukey’s test indicating significant differences among the three groups of sampling sites (p <0.05).

Table 4 Similarity of the metazoan parasite communities of the shoal flounder Syacium gunteri
Full size table

At the infracommunity level, there were significant differences in the qualitative similarity among the groups (ANOVA ONE-WAY, F[2, 7194] = 546.56, p <0.001), with Tukey’s test indicating that group 3 had significantly higher similarity values than groups 1 and 2 (p <0.05). For the quantitative similarity, there were also significant differences among the groups (ANOVA ONE-WAY, F[2, 7536] = 128.61, p <0.001), with Tukey’s test indicating that group 3 had significantly higher similarity values than groups 1 and 2 (p <0.05).

Statistical associations between environmental and biological variables

Table 5 shows the environmental variables and contaminants in the sediment and water selected by the forward stepwise procedure of CANOCO for all 33 sampling sites and for each of the three groups of sampling sites in Figure 2. These variables explained the greatest percentage of the variance in the redundancy analysis (Figure 3A-D). Figure 3A shows the statistical associations between the five environmental variables explaining the greatest amount of variance and the parasite infracommunity metrics. The most relevant pattern in Figure 3A was the negative association between the infracommunity metrics (especially Brillouin’s diversity index and the number of parasites per fish) and the concentrations of aliphatic hydrocarbons, aluminium and phosphorus. In contrast, the number of parasites per fish had a positive association with the concentration of clay in the sediments and an unknown environmental variable acting at a spatial scale of 58 km (PCNM58KM). Figure 3B-D show the statistical associations between the environmental variables explaining the greatest amount of variance and the infracommunity metrics for groups 1, 2 and 3, respectively. Figure 3B shows the negative statistical associations of nitrate, vanadium and the unresolved complex mixture with Brillouin’s diversity index and the number of parasite species per fish. For these two community metrics, the concentration of cobalt was apparently irrelevant. The number of parasites per fish and the Berger Parker dominance index were both associated positively with clay and negatively with nitrate and vanadium. Figure 3C shows a negative statistical association between total PAHs, porosity and phosphorus and most of the infracommunity metrics. The exception was the Berger-Parker dominance index, which apparently increased with the concentration of phosphorus. Figure 3D shows that for group 3, the number of individuals per fish had a negative association with the alkalinity and pH of the sediment. Additionally, with an increasing silt concentration in the sediments, the number of parasite species per fish and the Brillouin’s diversity index increased, and Berger-Parker index values decreased.

Table 5 Selected environmental variables from the sampling sites in the Gulf of Mexico
Full size table
Figure 3
figure 3

A-D Redundancy analysis of the environmental variables and infracommunity metrics of the parasites of Syacium gunteri . (A) The RDA for all sampling sites accounted for 36% of the total variance, and was highly significant for the first and all four canonical axes (F = 51.46; P-value = 0.0002; 4999 permutations) for the 33 sampling stations. (B) For group 1, the RDA accounted for 35% of the total variance, and was highly significant for the first and all four canonical axes (F = 51.46; P-value = 0.0002; 4999 permutations). (C) For group 2, the RDA accounted for 23% of the total variance, and was highly significant for the first and all four canonical axes (F = 14.10; P-value = 0.03; 4999 permutations). (D) For group 3, the RDA accounted for 68% of the total variance, and was highly significant for the first and all four canonical axes (F = 38.24; P-value = 0.03; 4999 permutations).

Full size image

Discussion

The results show that there were significant statistical associations between environmental variables and contaminants and the parasite community metrics of the shoal flounder. Therefore, the environmental impacts primarily produced by river run-off and polyaromatic hydrocarbons in the study area apparently have a direct effect on the parasite community metrics of S. gunteri. The sampling sites in groups 1 and 2 had a relatively high number of species, high diversity values and low qualitative and quantitative similarity at the component and infracommunity levels (Tables 2 and 3). Thus, despite the presence of pollutants (Table 5), the environmental conditions in both groups of sites were apparently not bad enough to prevent the life cycles from being completed. Consequently, we considered these groups of sampling sites as mildly disturbed. In contrast, as detailed below, several of the component and infracommunity metrics of the sites of group 3 presented values suggesting alterations related to environmental impact, including a low number of species, high numerical dominance and high similarity within sites (Tables 2 and 3). These characteristics suggest that parasite community metrics, and hence the probability of completion of the life cycles, may have been affected. Consequently, we considered the sampling sites in group 3 as strongly disturbed.

Prevalence and abundance of the most frequent and abundant species

All of the metazoan parasites recovered were generalist autogenic parasites that have been reported in marine and brackish fishes from other families in the region, namely, the red grouper Epinephelus morio[33], the grey snapper Lutjanus griseus[34], the tonguefish Symphurus plagiusa[35], the Florida pompano Trachinotus carolinus[36], the Mayan sea catfish Ariopsis assimilis[13] and the Mexican flounder Cyclopsetta chittendeni[37].

The most frequent and abundant helminth species infecting the shoal flounder across all 33 sampling sites were the larval cestode O. wageneri and the larval nematode P. decipiens (Table 1). The most likely explanation for the presence of these larval species throughout the study area is that they use copepods, amphipods, shrimp, and isopods as intermediate hosts, which in turn form part of the diet of S. gunteri[38]-[41]. In the case of P. decipiens, the adults of this nematode have been reported in pygmy whale (Kogia breviceps) along the coast of the Yucatan Peninsula [42]. In the case of O. wageneri, the putative definitive hosts are the stingrays, Dasyatis centroura[40].

Table 1 shows the low number of metazoan parasite species found at the sites of group 3, where only 4 species were recovered, versus sites in groups 1 and 2, where 10 and 12 parasite species were present, respectively. The most likely explanation for this difference is the absence of suitable conditions for the completion of the life cycle of all other parasite species in the sites of group 3. This is evidently the case for adult helminths such as Lecithochirium sp., the nematode Capillariidae gen. sp., the acanthocephalan Serrasentis sagittifer, the larval digenean Stephanostomum sp. and the copepod Caligus pomacentrus. These results suggest that a number of environmental characteristics prevent the presence of intermediate hosts in the sites of group 3. Similar effects of anthropogenic environmental impacts on the parasite populations of flatfishes have been reported for European waters [43],[44] and for the Mexican flounder Cyclopsetta chittendeni from the Campeche Sound, Gulf of Mexico [37].

Component communities

The parasite species composition of the component communities of the shoal flounder was formed by a large number of species in the larval stage (8/13 = 61%), which is most likely because this fish species is a benthic predator feeding on shrimps and copepods [41]. The shoal flounder had a lower number of parasite species compared with C. chittendeni (17 species; [37]) but a higher number compared with S. plagiusa in the Gulf of Mexico (8; [35]). In comparison with flatfishes from other regions, the number of parasite species of the shoal flounder was in the range reported for Spain and Germany for Platichthys flesus (13 and 17 species, respectively) [45],[46]. However, there were fewer parasite species in S. gunteri than in Hippoglossina macrops and Paralichthys adspersus from the coasts of Brazil and Chile, which harbour 20 and 22 helminth species, respectively [47],[48].

The overall diversity of the parasite component communities of S. gunteri (Table 2) is similar to that of S. plagiusa (0.24 ± 0.11) reported in the Gulf of Mexico [35]. However, the diversity of the component community for the sites of group 3 is the lowest obtained in the region to date for any flatfish species (see Table six in [37]). The statistically significant low values of diversity, the mean number of metazoan parasite species per site and the significantly high values of numerical dominance and similarity for the sites of group 3 (Table 2) were most likely due to the large number of larval parasites and the unsuitable conditions for the completion of the life cycle of the parasites using intermediate hosts. Further support for this argument comes from the study of benthic bioindicators, in which benthic fauna affected by oil spills or organic matter deposition undergo changes in their community composition from contaminant-sensitive species to contaminant-resistant species [49],[50]. Because the intermediate hosts (e.g., mollusks, shrimps) of the metazoan parasites are also residents of the benthic zone, their removal could negatively affect the parasite life cycles at the sites of group 3. In contrast, the completion of the life cycle of all other parasites in Table 1 was apparently not affected at the sites of the less environmentally disturbed groups 1 and 2. With increasing levels of environmental contamination, similar patterns of decreasing species numbers and parasite community metrics in flounder (Platichthys flesus) and of the infection parameters of digeneans in Mayan sea catfish (Ariopsis assimilis) have been reported by Schmidt et al. [46] and Vidal-Martínez et al. [13] , respectively.

Infracommunities

The mean number of metazoan parasite species per fish of the shoal flounder (Table 3) was within the range previously described for the Mexican flounder C. chittendeni (1.19 ± 1.17 to 3.27 ± 1.64 species per fish; [37]) and for the tonguefish S. plagiusa (2.00 ± 0.55 to 3.00 ± 0.66; [35]). These results suggest that the shoal flounder individuals are acquiring similar numbers of metazoan parasite species compared to C. chittendeni and S. plagiusa in the Gulf of Mexico, which is not surprising since these flatfish species share the same regional pool of parasites. However, once again, sites of group 3 had the fewest species per fish in the present study (Table 3). The mean number of individual parasites per fish (Table 3) in the shoal flounder was also within the range found for S. plagiusa in the Campeche Sound (47.00 ± 31.81) [35]. However, the identity of the numerically dominant species was different, with O. wageneri for S. gunteri and the larval digenean Stephanostomum sp. for S. plagiusa. This difference in the identity of the numerically dominant species reflects the oceanic versus coastal character of the parasite fauna of the shoal flounder and the tonguefish, respectively.

The diversity values for the parasite infracommunities of S. gunteri fall within the range commonly observed in the Campeche Sound, but the diversity values in sites of group 3 were the lowest obtained in the present study (see [35], and Table six in [37]). The statistically significant low diversity and mean numbers of metazoan parasite species per fish and the significantly high values of numerical dominance and similarity of sites of group 3 (Table 3) were most likely due to the large numbers of larval parasites and the unsuitable conditions for the parasites’ life cycles. Similar numerical dominance and high similarity values at the infracommunity level were observed in the parasite communities of the Mexican flounder C. chittendeni[37].

The coincidence in the numerical dominance of larval parasites of the infracommunities of the three flatfish species studied in the Campeche Sound, namely, the larval cestode O. wageneri in S. gunteri, the larval nematode Ascarididae gen. sp. in C. chittendeni and the larval digenean Stephanostomum sp. in S. plagiusa, is remarkable. Whether this pattern reflects a transition in the parasite community composition of the three flatfishes from contaminant-sensitive parasite species to contaminant-resistant parasite species is difficult to assess. However, because the three flatfish species were captured in the Campeche Sound, a zone influenced by oil extraction activities and the discharge of the Grijalva-Usumacinta River, the community metrics of their metazoan parasite fauna likely reflect the effect of living in a highly disturbed environment. To undertake such a comparison, the parasite communities of flatfishes in other, less-impacted environments in the Gulf of Mexico must be assessed.

Qualitative and quantitative similarity among groups at the component and infracommunity levels

At the component community level, the highest qualitative and quantitative similarities were observed between sites of groups 1 and 3 (Table 4), although the similarity values were rather low (Table 4). The most likely explanation for these similarities is that there were only a few parasite species that were shared between the sites. Clearly, this is the case for the sites of groups 1 and 3 since there were more parasite species in the former than in the latter (Table 1). These low similarities at the component community level are typical of diverse parasite communities. For example, flatfish species such as Pleuronectes flesus, Hippoglossoides platessoides, Gyptocephalus cynoglossoides and Hippoglossus hippoglossus had quantitative similarities between 24 ± 25 to 53 ± 19 for intestinal parasite communities in Norway [43],[51].

At the infracommunity level, the qualitative and quantitative similarities were very high among all three groups of sampling sites, especially for Jaccard’s similarity index (Table 4). The most likely explanation for this pattern is that at this hierarchical level, the individual fish were sharing the most frequent and abundant species of larval parasites (Table 1). This pattern can be explained based on the concept of “seed rain” [52], in which the larvae of the parasites become dispersed by ocean or tidal coastal currents, infecting intermediate hosts, which are eventually eaten by the shoal flounder at all of the study sites. This larval dispersion process has been well documented for free-living stages of digeneans (cercariae) infecting benthic intermediate hosts (crabs) [52]-[54].

Statistical associations between environmental variables and biological variables

The positive statistical associations in Figure 3A between both the number of parasites per fish and the Berger-Parker dominance index of the 33 sampling sites and the concentrations of aliphatic hydrocarbons, aluminium and phosphorus suggest a positive effect of these chemicals on these infracommunity metrics. The concentrations of aliphatic hydrocarbons (range 0–1.28 μg.g−1) and aluminium (range 0.24-18.52 μg.g−1; Table 5) found in the present study did not exceed the normal average concentrations in deep-sea sediments (10 μg.g−1 for aliphatics [55]; 95 μg.g−1 for aluminium [56]). A potential explanation for these associations is that the presence of aliphatic hydrocarbons and phosphorus, transported from the continent through rivers and then discharged into the marine sediments, enhances the growth of hydrocarbonoclastic bacteria [57]. The increase in these bacteria would in turn enhance primary and secondary productivity in the area, with a consequent increase in the number of intermediate hosts, as has been suggested for regions affected by oil spills such as the Prestige [8],[9]. The Brillouin’s diversity index had a negative statistical association with all of the environmental variables in Figure 3A. Whether some of these parasite species are sensitive to contamination remains an open question, however the combination of the environmental variables in Figure 3A is apparently having negative effects on this parasite community metric. The number of parasites per fish also had a positive association with the concentration of clay in the sediments and other unknown environmental variables acting at a spatial scale of 58 km (PCNM58KM), which could presumably be ocean currents transporting infective stages. Thus, an environment with a low concentration of aliphatic hydrocarbons, aluminium and clay seems to allow the intermediate hosts to survive, which in turn enhances transmission within these parasite communities at a spatial scale of 58 km.

The negative statistical associations in group 1 (Figure 3B) between nitrate, vanadium and the unresolved complex mixture, and the values of the Brillouin’s diversity index and the numbers of parasites and parasite species per fish suggest a deleterious effect of these chemicals on the infracommunity metrics. The nitrate concentrations in sites of group 1 (Figure 3B) were between 0.06 and 15.76 μM (Table 5), and concentrations of this nutrient greater than 0.59 μM are generally considered harmful for aquatic organisms [58]. This value was exceeded at 8 out of 14 sampling sites in group 1 (data not shown). The most likely negative effect would be on intermediate hosts such as mollusks or benthic crustaceans. The problem is apparently not extreme at present because the parasite community still has a relatively large number of species and values of diversity similar to those of parasite communities in the region (see Table six in [37]). However, if this uncontrolled sewage discharge continues, hypoxia similar to that caused by the nutrients deposited in the Mississippi river [59] can be expected. For the sites in group 1, the concentrations of vanadium, the unresolved complex mixture (UCM) and cobalt were between 0.02 and 4.0 μg.g−1, 0.72 and 30.98 μg.g−1 and 0.05 and 0.49 μg.g−1, respectively (Table 5). None of these concentrations exceeded the sediment benchmarks for aquatic life proposed for vanadium (57 μg.g−1[60]), UCM (see Table one in [61]), or cobalt (50 mg.g−1 for aquatic organisms) by the EPA [62]. The percentage of clay in the sediments (0-55%) is apparently adequate for the completion of parasite life cycles.

In the sites of group 2 (Figure 3C), there were negative statistical associations between porosity and total concentrations of hydrocarbons (total PAHs), phosphorus and zinc and all of the parasite infracommunity metrics (with the exception of the Berger Parker dominance index). The increase in porosity could theoretically have a positive effect on the probability of the presence of invertebrates [63]. If this is the case, predators of the potential intermediate hosts could reduce the probability of parasite transmission. In addition, the total PAH concentrations for the sites of group 2 ranged from 0.68 to 2.15 μg.g−1 (Table 5), whereas the threshold effect level (TEL) of total PAHs reported by NOAA (National Oceanic and Atmospheric Administration) is 1.68 μg.g−1 for marine organisms. In group 2, 4 out of 10 sampling sites (40%) slightly exceeded this threshold. The phosphorus concentrations ranged from 1.09 to 6.50 μM, and the zinc concentrations ranged from 0 to 1.86 μg−1 (Table 5). These values were low; the Canadian Water Quality Guidance states that water with 129.15 μM of phosphorus is oligotrophic [64]. In the same way, the zinc concentrations were below the threshold effect level for zinc in sediments (123 μg.g−1; [65]). However, a possible synergistic effect, especially with porosity and total PAHs, cannot be ruled out because higher porosity indicates more interstitial space for PAHs, phosphorus and zinc. Our results are in accordance with those of Vidal-Martínez et al. [14], who found negative correlations between the number of symbionts in the gills of pink shrimp (Penaeus duorarum) and the UCM concentrations in sediments from the Campeche Sound. UCM contains both aliphatic and polyaromatic hydrocarbons, which in turn are related to the total PAHs present in the sediments [55].

The results for the sites of group 3 (Figure 3D) suggest that there was a decrease in the number of individual parasites per fish with decreasing alkalinity and silt proportion. The range of alkalinity was between 0.84 and 2.92 meq.L-1 (Table 5), within the normal range for the Gulf of Mexico (2.46-2.54 meq.L-1; [66]) and other saline-alkaline environments (2.5 meq.L-1; [67]). The proportion of silt was also within the range found in the southern Gulf of Mexico [68]. Therefore, the alkalinity and the silt proportions for the study area do not seem to pose any threat to the parasites. Other community metrics in Figure 3D, for example the parasite species per fish, increased with increasing silt proportion. Brillouin’s diversity index also increased with the alkalinity values. The pH range in the sediments was small, between 7.45 and 8.33 (Table 5). However, our results suggest that there were very few parasites with complex life cycles in the sites of group 3 (Table 1). Thus, apparently this narrow range of pH could be enough to create an inhospitable environment for parasites with complex life cycles. In further support of this hypothesis, pH values between 7.7 and 7.9 increased the mortality of larval brittlestar Ophiothrix fragilis in comparison with the control pH value of 8.1 [69]. For parasites, Koprivnikar et al. [70] reported an interaction between pH (range = 7.8 to 8.2), salinity and time affecting the survival of the larval stages of the digenean Acanthoparyphium spinulosum. Thus, the potential interaction between pH and the decrease in the number of parasites cannot be disregarded without experimental evidence.

An equally possible explanation is that there were other contaminants that were not included in the present study and that could be acting alone or synergistically to negatively affect the parasite community metrics of group 3. This unknown contaminant could be co-varying strongly with both alkalinity and pH in the sediment, which could be the reason that these two variables appear significant in the analysis. Potential candidates for this unknown contaminant are PCBs and pesticides, as negative statistical associations have been found between these contaminants and the number of gill symbionts of pink shrimp in the Campeche Sound [14].

Conclusions

Our results suggest that different combinations of environmental and contaminant variables affect the three groups of sites detected and that this influence is spatially highly patchy, acting at different spatial scales. This finding poses a very complex scenario in which parasites could be responding specifically to different combinations of contaminants. Most of the environmental variables did not exceed the limits established by different agencies (EPA, NOAA), with the exception of nitrate and total PAHs, which were above the values established in the literature, although the number of sites above the threshold was low in most cases. Thus, both the shoal flatfish and its parasites apparently live in tolerable levels of these hydrocarbons and nutrients in most cases. In the specific case of nitrate, if this uncontrolled sewage discharge continues, the occurrence of hypoxia similar to that caused by the Mississippi river [59] can be expected. In examining sites of group 3 specifically, other contaminants such as PCBs and pesticides should be included to obtain a more complete picture of the potential contaminants affecting the parasite community metrics. The unexpected combined effect of alkalinity and pH upon the parasites is very suggestive, and we speculate that the potential interaction between pH and the decrease in the number of parasites cannot be disregarded. The community metrics chosen generally had robust statistically significant associations with both physicochemical and contaminant variables, which supports the ecological relevance of these parameters as indicators of aquatic environmental health [5],[71].

Additional file

References

  1. Lafferty KD: Environmental parasitology: what can parasites tell us about human impacts on the environment?. Parasitol Today. 1997, 13: 251-255. 10.1016/S0169-4758(97)01072-7.

    Article  CAS  PubMed  Google Scholar 

  2. Pietrock M, Marcogliese DJ: Free-living endohelminth stages: at the mercy of environmental conditions. Trends Parasitol. 2003, 19: 293-299. 10.1016/S1471-4922(03)00117-X.

    Article  PubMed  Google Scholar 

  3. Khan RA: Parasites of fish as biomarkers of environmental degradation: a field study. Bull Environ Contam Toxicol. 2004, 72: 394-400. 10.1007/s00128-003-9092-6.

    Article  CAS  PubMed  Google Scholar 

  4. Hudson PJ, Dobson AP, Lafferty KD: Is a healthy ecosystem one that is rich in parasites?. Trends Ecol Evol. 2006, 21: 381-385. 10.1016/j.tree.2006.04.007.

    Article  PubMed  Google Scholar 

  5. Vidal-Martinez VM, Pech D, Sures B, Purucker ST, Poulin R: Can parasites really reveal environmental impacts?. Trends Parasitol. 2010, 26: 44-51. 10.1016/j.pt.2009.11.001.

    Article  PubMed  Google Scholar 

  6. Holt EA, Miller SW: Bioindicators: using organisms to measure environmental impacts. Nat Educ Knowl. 2011, 3: 8-

    Google Scholar 

  7. Vidal-Martinez VM: Helminths and protozoans of aquatic organisms as bioindicators of chemical pollution. Parassitologia. 2007, 49: 77-84.

    Google Scholar 

  8. Pérez-del Olmo A, Raga JA, Kostadinova A, Fernández M: Parasite communities in boops boops (L.) (Sparidae) after the prestige oil-spill: detectable alterations. Mar Pollut Bull. 2007, 54: 266-276. 10.1016/j.marpolbul.2006.10.003.

    Article  PubMed  Google Scholar 

  9. Pérez-del-Olmo A, Fernández M, Raga JA, Kostadinova A, Morand S: Follow-up trends of parasite community alteration in a marine fish after the prestige oil-spill: shifting baselines?. Environ Pollut. 2009, 157: 221-228. 10.1016/j.envpol.2008.07.010.

    Article  PubMed  Google Scholar 

  10. Valtonen ET, Holmes JC, Koskivaara M: Eutrophication, pollution and fragmentation: effects on parasite communities in roach (rutilus rutilus) and perch (Perca fluviatilis) in four lakes in central Finland. Can J Fish Aquat Sci. 1997, 54: 572-585. 10.1139/f96-306.

    Article  Google Scholar 

  11. Huspeni TC, Lafferty KD: Using larval trematodes that parasitize snails to evaluate a saltmarsh restoration project. Ecol Appl. 2014, 14: 795-804. 10.1890/01-5346.

    Article  Google Scholar 

  12. Blanar CA, Munkittrick KR, Houlahan J, Maclatchy DL, Marcogliese DJ: Pollution and parasitism in aquatic animals: a meta-analysis of effect size. Aquat Toxicol. 2009, 93: 18-28. 10.1016/j.aquatox.2009.03.002.

    Article  CAS  PubMed  Google Scholar 

  13. Vidal-Martínez VM, Aguirre-Macedo ML, Noreña- Barroso E, Gold-Bouchot G, Caballero-Pinzón PI: Potential interactions between the metazoan parasites of the Mayan catfish Ariopsis assimilis and chemical pollution in Chetumal Bay, Mexico. J Helminthol. 2003, 77: 173-184. 10.1079/JOH2002158.

    Article  PubMed  Google Scholar 

  14. Vidal-Martinez VM, Aguirre-Macedo ML, Del Rio-Rodríguez R, Gold-Bouchot G, Rendón-von Osten J, Miranda-Rosas GA: The pink shrimp Farfantepenaeus duorarum, its symbionts and helminths as bioindicators of chemical pollution in Campeche sound, Mexico. J Helminthol. 2006, 80: 159-174. 10.1079/JOH2006358.

    Article  CAS  PubMed  Google Scholar 

  15. Pech D, Vidal-Martínez VM, Aguirre-Macedo ML, Gold-Bouchot G, Herrera-Silveira J, Zapata-Pérez O, Marcogliese D: The checkered puffer (Spheroides testudineus) and its helminths as bioindicators of chemical pollution in Yucatan coastal lagoons. Sci Total Environ. 2009, 407: 2315-2324. 10.1016/j.scitotenv.2008.11.054.

    Article  CAS  PubMed  Google Scholar 

  16. García-Cuellar A, Arreguín-Sánchez F, Hernández-Vázquez S, Lluch-Cota DB y: Impacto ecológico de la industria petrolera en la Sonda de Campeche, México, tras tres décadas de actividad: una revisión. Interciencia. 2004, 29: 311-319.

    Google Scholar 

  17. Gold-Bouchot G, Barroso-Noreña E, Zapata-Perez O: Hydrocarbon concentrations in the american oyster (Crassostrea virginica) in Laguna de Terminos, Campeche, Mexico. B Environ Contam Tox. 1995, 53: 222-227.

    Google Scholar 

  18. Gold-Bouchot G, Zavala-Coral M, Zapata-Pérez O, Ceja-Moreno V: Hydrocarbon concentrations in oysters (Crassostrea virginica) and recent sediments from three coastal lagoons in Tabasco, Mexico. B Environ Contam Tox. 1997, 59: 430-437. 10.1007/s001289900496.

    Article  CAS  Google Scholar 

  19. Strickland JDH, Parsons TR: A practical Handbook of Seawater Analysis. 1972, Fisheries Research Board of Canada Bulletin, Otawa

    Google Scholar 

  20. Comisión Oceanográfica Intergubernamental: Determinacion de los Hidrocarburos del Petróleo en los Sedimentos. Manuales y Guías 11. 1982, Comisión Oceanográfica Intergubernamental de la UNESCO (IOC-UNESCO), Paris

    Google Scholar 

  21. UNEP/FAO/IAEA/IOC: Sampling of Selected Marine Organisms and Sample Preparation for the Analysis of Chlorinated Hydrocarbons. Reference Methods for Marine Pollution Studies No. 12 Rev. 2. 1991, Intergovernmental Oceanographic Commission of UNESCO (IOC-UNESCO), Paris

    Google Scholar 

  22. Borcard D, Legendre P, Drapeau P: Partialling out the spatial component of ecological variation. Ecology. 1992, 73: 1045-1055. 10.2307/1940179.

    Article  Google Scholar 

  23. Borcard D, Legendre P, Avois-Jacquet C, Tuomisto H: Dissecting the spatial structure of ecological data at multiple scales. Ecology. 2004, 85: 1826-1832. 10.1890/03-3111.

    Article  Google Scholar 

  24. Santana-Piñeiros AM, Pech D, Vidal-Martínez VM: Spatial structure of the helminth parasite communities of the tonguefish, Symphurus plagiusa, from the Campeche coast, southern Mexico. Int J Parasitol. 2012, 42: 911-920. 10.1016/j.ijpara.2012.07.008.

    Article  Google Scholar 

  25. Vidal-Martínez VM, Aguirre-Macedo ML, Scholz T, González-Solís D, Mendoza-Franco E: Atlas of the helminth parasites of cichlid fishes of Mexico. 2001, Academia, Praha

    Google Scholar 

  26. Bush AO, Lafferty KD, Lotz JM, Shostak AW: Parasitology meets ecology on its own terms: Margolis et al. Revisited. J Parasitol. 1997, 83: 575-583. 10.2307/3284227.

    Article  CAS  PubMed  Google Scholar 

  27. Holmes JC, Price PW: Communities of Parasites. Community Ecology: Patterns and Processes. Edited by: Anderson DJ, Kikkawa J. 1986, Blackwell, Scientific Publications, Oxford, 187-213.

    Google Scholar 

  28. Esch GW, Fernandez J: A Functional Biology of Parasitism. 1993, Chapman & Hall, London

    Book  Google Scholar 

  29. Magurran AE: Ecological Diversity and its Measurements. 1988, Princeton, Princeton University Press

    Book  Google Scholar 

  30. Sokal RR, Rohlf FJ: Biometry. 1995, W.H. Freeman and Company, San Francisco, 2

    Google Scholar 

  31. Jongman RHG, Ter Braak CJF, van Tongeren OFR: Data Analysis in Community and Landscape Ecology. 1995, Cambridge University Press, Cambridge

    Book  Google Scholar 

  32. ter Braak CJF, Šmilauer P: CANOCO Reference manual and users guide to Canoco for Windows: software for Canonical Community Ordination (version 4). 1998, Microcomputer Power, New York

    Google Scholar 

  33. Moravec F, Vidal-Martínez VM, Vargas-Vazquez J: Helminth parasites of Epinephelus morio (Pisces: Serranidae) of the Yucatan peninsula, Mexico. Folia Parasit. 1997, 44: 255-266.

    CAS  Google Scholar 

  34. Argáez-García N, Guillén-Hernández S, Aguirre-Macedo ML: Intestinal helminths of Lutjanus griseus (Perciformes: Lutjanidae) from three environments in Yucatán (Mexico), with a checklist of its parasites in the Gulf of Mexico and Caribbean region. Rev Mex Biodivers. 2010, 81: 903-912.

    Google Scholar 

  35. Rodríguez-Gonzales A, Vidal-Martínez VM: Las comunidades de helmintos del lenguado (Symphurus plagiusa) en la costa de Campeche, México. Rev Mex Biodivers. 2008, 79: 159-173.

    Google Scholar 

  36. Sánchez-Ramírez C, Vidal-Martínez VM: Metazoan parasite infracommunities of Florida pompano (Trachinotus carolinus) from the coast of the Yucatan Peninsula, Mexico. J Parasitol. 2002, 88: 1087-1094. 10.1645/0022-3395(2002)088[1087:MPIOFP]2.0.CO;2.

    Article  PubMed  Google Scholar 

  37. Centeno-Chalé OA: Efectos potenciales del impacto ambiental en la Sonda de Campeche sobre las comunidades de helmintos parásitos del “lenguado Mexicano” Cyclopsetta chittendeni Bean, 1985. MSc Thesis. 2012, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Unidad Mérida, Departamento de Recursos del Mar, Yucatán, México, ᅟ

    Google Scholar 

  38. Marcogliese DJ: Transmission of the sealworm, Pseudoterranova decipiens (Krabbe) from invertebrates to fish in an enclosed brackish pond. J Exp Mar Biol Ecol. 1996, 205: 205-219. 10.1016/S0022-0981(96)02613-5.

    Article  Google Scholar 

  39. Palm HW: Ecology of Pseudoterranova decipiens Krabbe, 1878 (Nematoda: Anisakidae) from Antarctic waters. Parasitol Res. 1999, 85: 638-646. 10.1007/s004360050608.

    Article  CAS  PubMed  Google Scholar 

  40. Toth LM, Campbell RA, Schmidt GD: A revision of Oncomegas Dollfus, 1929 (Cestoda: Trypanorhyncha; Eutetrarhynchidae), the description of two new species and comments on its classification. Syst Parasitol. 1992, 22: 167-187. 10.1007/BF00009664.

    Article  Google Scholar 

  41. Sanchez-Gil P, Arreguin-Sánchez F, García-Abad MC: Ecological strategies and recruitment of Syacium gunteri (Pisces: Bothidae) in the southern Gulf of Mexico shelf. Neth J Sea Res. 1994, 32: 433-439. 10.1016/0077-7579(94)90019-1.

    Article  Google Scholar 

  42. González-Solís D, Vidal-Martínez VM, Antochiw-Alonso DM, Ortega-Argueta A: Anisakid Nematodes From Stranded Pygmy Sperm Whales, Kogia breviceps (Kogiidae), in Three Localities of the Yucatan Peninsula, Mexico. J Parasitol. 2006, 92: 1120-1122. 10.1645/GE-3553RN.1.

    Article  Google Scholar 

  43. Lile NK: Alimentary tract helminths of four pleuronectid flatfish in relation to host phylogeny and ecology. J Fish Biol. 1998, 53: 945-953. 10.1111/j.1095-8649.1998.tb00455.x.

    Article  Google Scholar 

  44. Marques JF, Santos MJ, Cabral HNI: Soleidae macroparasites along the Portuguese coast: latitudinal variation and host-parasite associations. Mar Biol. 2006, 150: 285-298. 10.1007/s00227-006-0339-8.

    Article  Google Scholar 

  45. Álvarez F, Iglesias R, Parama AI, Leiro J, Sanmatin M: Abdominal macroparasites of commercially important flatfish (Teleostei: Scophthalmidae, Pleuronectidae, Soleidae) in northwest Spain (iCES IXa). Aquaculture. 2002, 213: 31-53. 10.1016/S0044-8486(02)00025-X.

    Article  Google Scholar 

  46. Schmidt V, Zander S, Körting W, Broeg K, Westernhagen V, Dizer H, Hansen H, Skouras A, Steinhagen D: Parasites of flounder (Platichthys flesus) from the German Bight, North Sea, and their potential use in ecosystem monitoring. C. Pollution effects on the parasite community and a comparison to biomarker responses. Helgoland Mar Res. 2003, 57: 262-271. 10.1007/s10152-003-0159-x.

    Article  Google Scholar 

  47. Oliva ME, Castro R, Burgos R: Parasites of the flatfish Paralichthys adspersus (Steindachner, 1867) (Pleuronectiformes) from northern Chile. Mem Inst Oswaldo Cruz. 1996, 91: 301-306. 10.1590/S0074-02761996000300009.

    Article  Google Scholar 

  48. González MT, Acuña E, Oliva ME: Metazoan parasite fauna of the bigeye flounder, Hippoglossina macrops, from northern Chile. Influence of host age and sex. Mem I Oswaldo Cruz. 2001, 96: 1049-1054.

    Google Scholar 

  49. Dauvin JC, Alizier S, Rolet C, Bakalemc A, Belland G, Gomez Gesteira JL, Grimes S, de-la-Ossa-Carreterog JA, Del-Pilar-Ruso Y: Response of different benthic indices to diverse human pressures. Ecol Indic. 2012, 12: 143-153. 10.1016/j.ecolind.2011.03.019.

    Article  Google Scholar 

  50. Rombouts I, Beaugrand G, Dauvin JC: Potential changes in benthic macrofaunal distributions from the English Channel simulated under climate change scenarios. Estuar Coast Shelf S. 2013, 99: 153-161. 10.1016/j.ecss.2011.12.026.

    Article  Google Scholar 

  51. National Aquaculture Legislation Overview. Norway: National Aquaculture Legislation Overview (NALO) Fact Sheets; http://www.fao.org/fishery/legalframework/nalo_norway/en.

  52. Zimmer RK, Fingerut JT, Zimmer CA: Dispersal pathways, seed rains, and the dynamics of larval behavior. Ecology. 2009, 90: 1933-1947. 10.1890/08-0786.1.

    Article  PubMed  Google Scholar 

  53. Fingerut JT, Zimmer CA, Zimmer RK: Patterns and processes of larval emergence in an estuarine parasite system. Biol Bull. 2003, 205: 110-120. 10.2307/1543232.

    Article  PubMed  Google Scholar 

  54. Fingerut JT, Zimmer CA, Zimmer RK: Larval swimming overpowers turbulent mixing and facilitates transmission of a marine parasite. Ecology. 2003, 84: 2502-2515. 10.1890/02-4035.

    Article  Google Scholar 

  55. Venturini N, Tommasi LR, Mícego MC, Martins CC: Characterization of the benthic environment of a coastal area adjacent to an oil refinery, Todos os Santos Bay (ne-brazil). Braz J Oceanogr. 2004, 52: 123-134. 10.1590/S1679-87592004000200004.

    Article  Google Scholar 

  56. Martin JM, Whitfield M: The Significance of the River Inputs to the Ocean. Trace Metals in Seawater. Edited by: Wong CS, Boyle E, KBruland W, Burton JD, Goldberg ED. 1983, Plenum Press, New York, 265-269. 10.1007/978-1-4757-6864-0_16.

    Chapter  Google Scholar 

  57. Juhasz AL, Naidu R: Bioremediation of high molecular weight polycyclic aromatic hydrocarbons: a review of the microbial degradation of benzo[a ]pyrene. Int Biodeter Biodegr. 2000, 45: 57-88. 10.1016/S0964-8305(00)00052-4.

    Article  CAS  Google Scholar 

  58. Camargo JA, Alonso A: Ecological and toxicological effects of inorganic nitrogen pollution in aquatic ecosystems: A global assessment. Environ Int. 2006, 32: 831-849. 10.1016/j.envint.2006.05.002.

    Article  CAS  PubMed  Google Scholar 

  59. Rabalais N: Nitrogen in aquatic ecosystems. Ambio. 2000, 31: 102-112.

    Article  Google Scholar 

  60. Sediment Benchmarks for Aquatic Life: http://www.epa.gov/bpspill/sediment-benchmarks.html.

  61. Scarlett A, Galloway TS, Rowland SJ: Chronic toxicity of unresolved complex mixtures (UCM) of hydrocarbons in marine sediments. J Soils Sediments. 2007, 7: 200-206. 10.1065/jss2007.06.232.

    Article  CAS  Google Scholar 

  62. Mid-Atlantic Risk Assessment: http://www.epa.gov/reg3hwmd/risk/eco/btag/sbv/fwsed/screenbench.htm.

  63. Bo T, Fenoglio S, Malacarne G, Pessino M, Sgariboldi F: Effects of clogging on stream macroinvertebrates: an experimental approach. Limnologica. 2007, 37: 186-192. 10.1016/j.limno.2007.01.002.

    Article  Google Scholar 

  64. Canadian Council of Ministers of the Environment: Canadian water quality guidelines for the protection of aquatic life: Phosphorus: Canadian Guidance Framework for the Management of Freshwater Systems. 2004, Canadian Council of Ministers of the Environment, Winnipeg

    Google Scholar 

  65. Burton GAJ: Sediment quality criteria in use around the world. Limnology. 2002, 3: 65-75. 10.1007/s102010200008.

    Article  CAS  Google Scholar 

  66. Environmental Protection Agency: Designation of a site in the Gulf of Mexico for incineration of chemical wastes. Final Environmental Impact Statement. EPA-EIS-WA 76X-054. 1976, Environmental Protection Agency, Texas

    Google Scholar 

  67. Lin T, Lai Q, Yao S, Lu J, Zhou K, Wang H: Combined effects of carbonate alkalinity and pH on survival, growth and haemocyte parameters of the Venus clam Cyclina sinensis. Fish Shellfish Immun. 2013, 35: 525-531. 10.1016/j.fsi.2013.05.006.

    Article  CAS  Google Scholar 

  68. Yáñez-Arancibia A, Lara-Domínguez AL, Rojas JL, Villalobos GJ, Zárate D, Sánchez-Gil P: Integrated coastal zone management plan for Terminos Lagoon, Campeche Mexico. The Gulf of Mexico Large Marine Ecosystem: Assessment, Sustainability and Management. Edited by: Kumpf H, Steidinger K, Sherman K. 1999, Blackwell Science, Massachusetts, 565-592.

    Google Scholar 

  69. Dupont S, Havenhand J, Thorndyke W, Peck L, Thorndyke M: Near-future level of CO2-driven ocean acidification radically affects larval survival and development in the brittlestar Ophiothrix fragilis. Mar Ecol-Prog Ser. 2008, 373: 285-294. 10.3354/meps07800.

    Article  CAS  Google Scholar 

  70. Koprivnikar J, Lim D, Fu C, Brack SHM: Effects of temperature, salinity, and pH on the survival and activity of marine cercariae. Parasitol Res. 2010, 106: 1167-1177. 10.1007/s00436-010-1779-0.

    Article  PubMed  Google Scholar 

  71. Adams SM: Biological Indicators of Aquatic Ecosystem Stress. 2002, American Fisheries Society, Bethesda, Maryland

    Google Scholar 

Download references

Acknowledgements

The authors thank the Gerencia de Seguridad Industrial y Protección Ambiental RMNE, PEMEX Exploración y Producción for financial support through contract no. 428 818 804 “Determination of potential biotic effects for contingency Kab 121 well”. They also thank Dr. David Valdes-Lozano for allowing us to use the physicochemical data of Xcambo 2 oceanographic campaign, and to Clara Vivas-Rodríguez, Gregory Arjona-Torres, and Francisco Puc-Itza of CINVESTAV-Mérida for support with the field and laboratory work.

Author information

Authors and Affiliations

  1. Laboratorios de Parasitología y Geoquímica, Centro de Investigación y de Estudios Avanzados del Instituto Politécnico Nacional, Unidad Mérida, Km 6 Carretera Antigua a Progreso, Cordemex, Mérida, 97310, Yucatán, México

    Víctor Manuel Vidal-Martínez, Oscar A Centeno-Chalé, Edgar Torres-Irineo, Juan Sánchez-Ávila, Gerardo Gold-Bouchot & M Leopoldina Aguirre-Macedo

Authors
  1. Víctor Manuel Vidal-Martínez
    View author publications

    You can also search for this author in PubMed Google Scholar

  2. Oscar A Centeno-Chalé
    View author publications

    You can also search for this author in PubMed Google Scholar

  3. Edgar Torres-Irineo
    View author publications

    You can also search for this author in PubMed Google Scholar

  4. Juan Sánchez-Ávila
    View author publications

    You can also search for this author in PubMed Google Scholar

  5. Gerardo Gold-Bouchot
    View author publications

    You can also search for this author in PubMed Google Scholar

  6. M Leopoldina Aguirre-Macedo
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Víctor Manuel Vidal-Martínez.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

VMVM conceived of the study, performed part of the statistical analysis and drafted the manuscript. ACC carried out the parasitological examinations, produced the figures and helped to draft the manuscript. ETI carried out part of the statistical analysis and helped to draft the manuscript. JSA helped to draft the manuscript. GGB provided the data on hydrocarbons and heavy metals, and helped to draft the manuscript. MLAM participated in the design of the study and coordination and helped to draft the manuscript. All authors read and approved the final manuscript.

Electronic supplementary material

13071_2014_541_MOESM1_ESM.xlsx

Additional file 1: Environmental variables of sediments and water from the sampling sites in the Gulf of Mexico. Codes are as follows: BPM/APM = ratio of PAHBPM and PAHAPM, PAHAPM = polycyclic aromatic hydrocarbons of high molecular weight (4–5 benzene rings), PAHBPM = polycyclic aromatic hydrocarbons of low molecular weight (2–3 benzene rings), PAHTOT = total polycyclic aromatic hydrocarbons, PCNM2, PCNM15, PCNM 21 and PCNM58K correspond to spatial variables acting at the spatial scales of 2, 15, 21 and 58 km respectively, UCM = unresolved complex mixture. (XLSX 24 KB)

Authors’ original submitted files for images

Rights and permissions

Open Access  This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made.

The images or other third party material in this article are included in the article’s Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder.

To view a copy of this licence, visit https://creativecommons.org/licenses/by/4.0/.

The Creative Commons Public Domain Dedication waiver (https://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Vidal-Martínez, V.M., Centeno-Chalé, O.A., Torres-Irineo, E. et al. The metazoan parasite communities of the shoal flounder (Syacium gunteri) as bioindicators of chemical contamination in the southern Gulf of Mexico. Parasites Vectors 7, 541 (2014). https://doi.org/10.1186/s13071-014-0541-3

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/s13071-014-0541-3

Keywords