Metazoan parasite communities: support for the biological invasion of Barbus barbus and its hybridization with the endemic Barbus meridionalis
Parasites & Vectors volume 9, Article number: 588 (2016)
Recently, human intervention enabled the introduction of Barbus barbus from the Rhône River basin into the Barbus meridionalis habitats of the Argens River. After an introduction event, parasite loss and lower infection can be expected in non-native hosts in contrast to native species. Still, native species might be endangered by hybridization with the incomer and the introduction of novel parasite species. In our study, we aimed to examine metazoan parasite communities in Barbus spp. populations in France, with a special emphasis on the potential threat posed by the introduction of novel parasite species by invasive B. barbus to local B. meridionalis.
Metazoan parasite communities were examined in B. barbus, B. meridionalis and their hybrids in three river basins in France. Microsatellites were used for the species identification of individual fish. Parasite abundance, prevalence, and species richness were compared. Effects of different factors on parasite infection levels and species richness were tested using GLM.
Metazoan parasites followed the expansion range of B. barbus and confirmed its introduction into the Argens River. Here, the significantly lower parasite number and lower levels of infection found in B. barbus in contrast to B. barbus from the Rhône River supports the enemy release hypothesis. Barbus barbus × B. meridionalis hybridization in the Argens River basin was confirmed using both microsatellites and metazoan parasites, as hybrids were infected by parasites of both parental taxa. Trend towards higher parasite diversity in hybrids when compared to parental taxa, and similarity between parasite communities from the Barbus hybrid zone suggest that hybrids might represent “bridges” for parasite infection between B. barbus and B. meridionalis. Risk of parasite transmission from less parasitized B. barbus to more parasitized B. meridionalis indicated from our study in the Argens River might be enhanced in time as higher infection levels in B. barbus from the Rhône River were revealed. Hybrid susceptibility to metazoan parasites varied among the populations and is probably driven by host-parasite interactions and environmental forces.
Scientific attention should be paid to the threatened status of the endemic B. meridionalis, which is endangered by hybridization with the invasive B. barbus, i.e. by genetic introgression and parasite transmission.
Concerns over the ecological implications of the introduction of an alien species into new environments are increasing. From the conservation point of view, there is elevated apprehension if the native species represents an endemic or endangered species, which is often characterized by small population sizes, fragmented distribution and low genetic variability (e.g. [1, 2]). A new incomer may, therefore, represent a serious problem for the resident if it predates the indigenous species, exploits the same resources, or alters its native habitat [3–6]. Following the introduction event, many native species are endangered through hybridization with the closely related alien species [7–9]. At the same time, the new invader serves as a source of novel parasites to which the local species may display a different degree of susceptibility [10, 11]. Consequently, exposing susceptible local hosts to new parasite species carried by introduced individuals may result in accelerated mortality in native populations. For instance, the parasitic nematode Anguillicoloides crassus was imported to Europe probably as a result of the introduction of the Japanese eel Anguilla japonica and was, subsequently, disseminated in the populations of the European eel Anguilla anguilla (reviewed in ). While A. crassus is not highly pathogenic in the Japanese eel probably due to low-intensity infection rates , high infection and more serious pathology connected with high mortalities may be detected in wild European eels [14–17].
Populations of invaders established in new habitats typically exhibit fewer parasite species, and a smaller number of host individuals are parasitized (i.e. there is a lower prevalence of infection) when compared to the source populations . This could be the result of the introduction of a restricted number of individuals carrying only a proportion of the original parasite fauna, new and unsuitable environmental conditions for parasites, and the absence or low abundance of suitable hosts required for the parasite life-cycles [19, 20]. Such a release from co-evolved parasites may therefore provide an advantage for the performance of a novel host species in new habitats . However, this advantage is often of a temporal nature and parasite species richness and prevalence rates can return to original levels or even be multiplied in a short time by transmission from the local hosts . Co-evolutionary relationships which evolve between hosts and their parasites during their co-existence may, therefore, be steered by ecological forces . Hybridization might even alter the composition of metazoan parasite communities of the two interacting host species, since hybrid individuals are often vulnerable to parasites infecting both parental species [24, 25] and may, therefore, represent “bridges” for parasite infection .
In France, two congeneric Barbus species co-exist in several rivers in the Mediterranean basin. The low level of mitochondrial DNA and allozyme variability indicates recent colonization of the French rivers by the common barbel Barbus barbus after the last glaciation [27, 28], where the Mediterranean barbel Barbus meridionalis was already present, probably from the Miocene . Nowadays, the widely distributed European species B. barbus has been found in almost all French river basins and prefers medium-sized to large rivers. By contrast, the occurrence of endemic B. meridionalis is restricted to the Languedoc-Roussillon, Rhône-Alpes, and Provence-Alpes-Côte d’Azur regions in France , where it inhabits mainly upper and middle streams of mountain rivers, probably as a result of competition with B. barbus . Nevertheless, hybridization between these two species has previously been reported from the Hérault, Garonne, Orb and Rhône river basins . Moreover, an irrigation canal which supplies the Var department with water collected from the River Durance (Rhône River basin) enabled the very recent immigration of B. barbus individuals into the habitats of the native B. meridionalis approximately 30 years ago [31, 32]. While Kiener et al.  in 1981 rejected the presence of B. barbus in the River Argens, hybridization between B. barbus and B. meridionalis is already occurring in this river .
In our study, we aimed to examine the composition of metazoan parasite communities of B. barbus and B. meridionalis in (i) allopatric areas, (ii) sympatric areas of late origin of the Rhône River basin, and (iii) sympatric areas of recent origin of the Argens River basin. In the Argens River basin, we further intended to confront the metazoan parasite communities of parental taxa with those found in hybrids which resulted from the biological invasion of B. barbus into this watershed. The metazoan parasite abundance, prevalence and species richness in Barbus populations collected from the three areas were analyzed. We focused on the possible threat posed by the introduction of the widely distributed B. barbus to native and endemic B. meridionalis with respect to the transmission of non-native parasite species and the role of hybrids in facilitating parasite transmission between parental taxa.
From 2007 to 2014, 349 B. barbus (BB), B. meridionalis (BM) and hybrid (H) individuals were collected in France from the allopatric BB and BM populations (site 1 on the Loire River basin and site 14 on the Argens River basin, respectively), the Rhône late sympatric populations (sites 2–3 on the River Ardèche and sites 4–8 on the River Durance of the Rhône River basin), and the Argens recent sympatric populations (sites 9–13 on the Argens River basin); (Fig. 1, Table 1). The full information on the site names and their coordinates are shown in Additional file 1: Table S1. To remove the effect of temporal variation, fish were sampled only in the summer period (i.e. July-August) when the highest parasite diversity and high abundance of many metazoan parasite species are expected. The water temperature was measured (in °C) in each locality. Fish were measured (standard body length in mm), transported to the laboratory, and subsequently examined for metazoan parasites.
Microsatellite genotyping of fish individuals
Microsatellite markers developed for the Barbus species [34–36] were used to identify the fish species in our study. Genomic DNA was isolated from fin clip samples stored in 96% ethanol using DNeasy Blood & Tissue Kit (Qiagen GmbH, Hilden, Germany). Further diluted DNA (approx. 10 ng/l) served as a template in the following multiplex PCR analysis of microsatellite loci following the protocol described in . Amplicons were analyzed on an ABI PRISM 3130 Genetic Analyzer (Applied Biosystems) using 500LIZ® Size Standard (Applied Biosystems) and Hi-Di™ Formamid (Applied Biosystems), and genotypes were finally scored using GeneMapper Software version 4.0 (Applied Biosystems). Due to the tetraploid genome of the investigated Barbus spp., only selected microsatellite loci with probable disomic inheritance  were applied in our study.
The MICRO-CHECKER program  was used to check for microsatellite null alleles in the Barbus populations and, subsequently, locus Barb21 was excluded from our study. Overall, 19 loci (Barbus2, Barbus22, Barbus26, Barbus28, Barbus31, Barbus32, Barbus36, Barbus37, Barbus40, Barbus41, Barbus47, Barbus49, Barbus50, Barbus56, Barbus62, Barbus63, Barbus65, Barb59 and Barb79) were used afterwards following a Bayesian clustering approach implemented in the STRUCTURE software . The program was run for five independent runs assuming an admixture model and the model of correlated allele frequencies, using 1,000,000 iterations after a burn-in period of 100,000 iterations for K = 2 clusters. The Introgress package  implemented in the R statistical software was used to calculate the hybridization index (h-index) in sympatric populations. Allopatric populations of the Loire River and the Argens River basins (i.e. sites 1 and 14) were set in this software as parental BB and BM populations, respectively. First, the interspecific differentiation index (D) between allelic frequencies in parental populations was computed. Subsequently, 14 microsatellites with D ≥ 0.80 (Barbus22, Barbus26, Barbus31, Barbus32, Barbus37, Barbus47, Barbus49, Barbus50, Barbus56, Barbus62, Barbus65, Barbus63, Barb59 and Barb79) were selected for estimation of the h-index in sympatric populations, as applied in Andrés et al. . In general, the resultant h-indices of zero and one should be used to determine pure individuals. Since the occurrence of BB in the Argens River basin is a result of introduction from the River Durance, we applied two approaches for BB designation (i) using an h-index of zero and (ii) using an h-index of up to 0.11. As similar results were obtained using both approaches (Additional file 2: Table S2; Additional file 3: Table S3; Additional file 4: Table S4), we presented only the results when the individuals with an h-index up to 0.11 were treated as BB and those with an h-index between 0.11 and 1 were treated as H individuals.
Quantitative and qualitative comparisons of metazoan parasite communities
Fish dissection was performed following Ergens & Lom . In our study, fins, gills, eyes, heart, kidney, spleen, hepathopancreas, intestine, gonads, gall-bladder, and swim bladder were examined for all metazoan parasites using a stereo microscope Olympus SZX7. Parasites were fixed as described in Lamková et al.  and, subsequently, identified using a light microscope (Olympus BX50) equipped with phase-contrast, differential interference contrast, and Olympus Stream Motion 1.9.2 digital image analysis software. Parasites were identified using the available identification keys and publications providing keys to the identification of metazoan parasites, e.g. [43–46]. Measures of parasite infection, i.e. prevalence and abundance of metazoan parasites were calculated according to Bush et al. . The individual abundance of the myxozoan parasites was not taken into consideration (because these parasites cannot be quantified as in the case of other metazoan parasites) and only the prevalence of this parasitic group was taken into account in the analyses. The effect of sampling effort on parasite species richness was corrected using the Chao1 estimator , and was calculated using the EstimateS program  on the basis of abundance data excluding data on myxozoan parasites. Similarities between metazoan parasite communities based on presence/absence data (Jaccard index) were computed in PAST .
Spearman’s rank correlation was computed between individual admixture q-values and h-index obtained by STRUCTURE software and Introgress package software, respectively. Metazoan parasite abundance and averaged prevalence were log-transformed prior to statistical analyses. Kolmogorov-Smirnov test was used for normality data assessment. Subsequently, Bonferroni post-hoc tests following General Linear Model (GLM) were used to compare the estimated marginal means of total abundance, prevalence and species richness adjusted for fish body length, water temperature, and sampling years between BB, BM and H from different river basins. Since all fish individuals were infected with at least one parasitic species (i.e. the overall prevalence was 100% in each fish population), the average of prevalences for each parasite species across the fish populations was used in our study (further referred to as averaged prevalence).
For the River Argens, GLM analyses were conducted to investigate the potential effects of different factors, i.e. host (BB, BM or H), locality (site), sampling year, water temperature, and host body length on the abundance, averaged prevalence, and species richness of metazoan parasites found in Barbus individuals.
Genetic composition of Barbus spp. populations
Based on microsatellite markers, STRUCTURE analysis confirmed the existence of one allopatric BB population in the Loire River basin (Site 1) and one allopatric BM population in the Argens River basin (Site 14). A low level of admixture between the populations of two Barbus species was revealed in the Rhône River basin (q < 0.12). In contrast, the extent of population admixture in the Argens River basin was high except for the abovementioned the Argens allopatric BM population (Fig. 2). The values of h-index obtained by the Introgress software package resembled the overall picture of the individual admixture obtained by STRUCTURE (Fig. 2) and correlated significantly with the q-values (Spearman’s rank correlation, r (373) = 0.97, P < 0.001). An h-index up to 0.11 was revealed in populations of the Rhône River basin. Using this value as an upper limit for B. barbus individuals, 25 B. barbus, 42 B. meridionalis and 81 hybrids were detected within the Argens River basin (i.e. a system with a very recent introduction of B. barbus), where we considered Barbus populations with the occurrence of hybrids as sympatric populations. In the Argens River basin, the co-existence of both pure species was documented only in Site 9. The Barbus sample originating from Site 13 was composed entirely of hybrid individuals (Table 1, Fig. 2).
Composition of metazoan parasite communities in Barbus spp. populations from different river basins
Examination of the studied Barbus spp. populations for metazoan parasites revealed the presence of parasites belonging to different parasitic groups (Myxozoa, Trematoda, Monogenea, Cestoda, Nematoda, Acanthocephala, Mollusca, Chelicerata and Crustacea). The prevalence and mean abundance of each metazoan parasite species are listed in Table 2. Detailed information on the composition of parasite communities per site is shown in Additional file 5: Table S5. Monogeneans (especially Dactylogyrus spp.) were the most dominant groups of parasites in the Loire allopatric BB and the Rhône late sympatric BB, while they were rare in the Argens recent sympatric BB. Above all, monogeneans and acanthocephalans (specifically Gyrodactylus spp. and Pomphorhynchus tereticollis, respectively) were dominant in the Argens recent sympatric BM and H. In a single Argens allopatric BM population, nematodes represented the most dominant and abundant parasite group (Table 2, Fig. 3).
Similarity in metazoan parasite communities between Barbus spp. populations
High Jaccard index value was found between the Argens recent sympatric BB and H individuals (0.71), while lower similarity in parasite communities (0.57) was found between the Argens recent sympatric BM and H. Parasite communities in the Argens sympatric BB were more similar to parasite communities in the Rhône late sympatric BB (0.52) than to those in the Loire allopatric BB (0.43). The parasite communities in the Argens recent sympatric BB and the Argens recent sympatric BM were more similar (0.50) than those in the Argens recent sympatric BB and the Argens allopatric BM (0.43; Table 3).
Total abundance, averaged prevalence and species richness of metazoan parasites in Barbus spp. populations
Significant differences in abundance, averaged prevalence, and species richness (GLM, abundance: whole model F (4,340) = 18.05, P < 0.001; averaged prevalence: whole model F (4,1241) = 1.70, P = 0.029; species richness: whole model F (4,340) = 23.24, P < 0.001) of metazoan parasites between Barbus groups were revealed (see Additional file 3: Table S3 and Additional file 4: Table S4 for detailed statistics). After controlling for the covariates, Bonferroni post-hoc tests revealed no significant differences in metazoan parasite abundance and averaged prevalence between the Loire allopatric BB and the Rhône late sympatric BB (P > 0.05). However, significantly lower species richness was revealed in the Rhône late sympatric BB than in the Loire allopatric BB (P = 0.008). Significantly lower values of abundance and species richness (P < 0.001), and lower but not significantly different averaged prevalence (P > 0.05) of metazoan parasites were found in the Argens recent sympatric BB in comparison with the Rhône late sympatric BB. In the Argens River basin, significantly lower metazoan abundance (P < 0.001) and species richness (P = 0.004), and lower but not significantly different averaged prevalence (P > 0.05) were found in BB when compared to BM. Metazoan parasite abundance, averaged prevalence, and species richness in H tended to be intermediate between pure species of the Argens River basin. However, statistically significant difference was only revealed in the case of metazoan parasite abundance between H and BM (P < 0.001); (Fig. 4, Additional file 4: Table S4).
Factors influencing parasitism in the B. barbus × B. meridionalis recent hybrid zone
In the Argens River basin, significant effects of locality (F (5,174) = 11.42, P < 0.001), sampling year (F (3,174) = 7.99, P < 0.001), water temperature (F (1,174) = 6.33, P = 0.013), and fish length (F (1,174) = 20.15, P < 0.001) on metazoan parasite abundance were revealed, while effect of host was not significant in GLM (whole model R 2 = 0.45, F (12,174) = 11.55, P < 0.001). Host (F (2,174) = 3.85, P = 0.012), locality (F (5,124) = 20.87, P < 0.001), sampling year (F (3,174) = 20.22, P < 0.001), water temperature (F (1,174) =5.32, P = 0.030), and fish length (F (1,174) = 4.16, P = 0.020) significantly affected the species richness of metazoan parasites (whole model R 2 = 0.74, F (12,174) = 32.88, P < 0.001). Neither host, locality, sampling year, water temperature, nor host body length effected significantly averaged prevalence of metazoan parasites (whole model R 2 = 0.02, F (12,709) = 1.01, P = 0.442).
This study explored composition of metazoan parasite communities in Barbus spp. populations of three river basins. No significant differences in the intensity of metazoan parasite infection was revealed between the Rhône late sympatric BB and the Loire allopatric BB. At the same time, lower number of parasitic taxa was revealed in the Rhône late sympatric BB than in the Loire allopatric BB. Low introgression of BM revealed by microsatellite markers in the Rhône late sympatric BB, therefore, does not provide a disadvantage in terms of high levels of parasite infection or high species richness in contrast to the Loire allopatric BB. Yet, a degree to which parasites have a negative impact on host vigour and fitness components is dependent on a particular parasite species, parasite genotype, or co-infection with other parasite species [51, 52]. This is, however, beyond the scope of our study. Several parasite species were shared between the Loire allopatric BB, the Rhône late sympatric BB, and the Argens recent sympatric BB, which indicates that these parasite species followed the expansion of BB. On the other hand, we detected considerably lower diversity and intensity of infection of metazoan parasites in the Argens recent sympatric BB in contrast to the Rhône late sympatric BB. Still, all parasite species except for Holostephanus sp. found in our study in the Argens recent sympatric BB were also present in the Rhône late sympatric BB (Table 2). These findings support the view that BB of the Argens River basin originates from the Durance River system. However, parasite species infecting the Argens recent sympatric BB represented only a small proportion of the metazoan parasite fauna found in the Rhône late sympatric BB (Table 2). The enemy release hypothesis suggests that individuals introduced outside their natural ranges may benefit from enemy release, e.g. predators or pathogens [53–55]. Our results concerning the significantly lower levels of parasite infection and lower number of metazoan parasites in the Argens recent sympatric BB when compared to the Rhône late sympatric BB are, therefore, in congruence with the general scenario of parasite loss after host introduction into a new environments . Kennedy & Bush  revealed that the parasite communities of native Onchorhynchus mykiss populations were dominated by specialist helminth parasites, while the number of specialist helminths declined with the increasing distance of translocated host populations from their original heartland. In fish, monogeneans are considered to be the most host-specific parasites . In our study, five of six monogenean species documented in the Rhône late sympatric BB were found also in the Argens recent sympatric BB. However, the absence of Dactylogyrus carpathicus, a highly abundant parasite species of the Rhône late sympatric BB, and the decrease in the abundance and prevalence of Dactylogyrus malleus in the Argens sympatric BB resulted in a shift in parasite communities from the dominance of monogenean parasites to a higher proportion of endoparasitic groups in total parasite numbers. Ondračková et al.  showed that the reduction in parasite numbers is dependent on the time after a colonization event and, therefore, a very late arrival of B. barbus into the Argens River tributaries can be expected. The first events involving the migration of BB from the River Durance to the waters of the Argens River system probably occurred in the period 1980–1990 . Since an upstream migration range of up to several dozens of kilometres has been documented for BB individuals , colonization of the Argens River tributaries might thus have taken place at the turn of this century or even more recently.
Parasites often exhibit a shorter generation time, larger population size, and higher migration and mutation rates than their hosts. Consequently, as a result of co-evolutionary host-parasite interactions, local parasite adaptation, i.e. the better performance of parasites in their local hosts than in foreign ones, is expected . Individuals introduced into novel areas may, therefore, benefit from the higher ability of parasites to adapt to their local hosts. By contrast, lower ability of parasites to infect their local hosts than alien ones or no differences in degrees of host resistance/susceptibility to parasites between local and non-indigenous hosts (i.e. no local adaptation), were previously documented [60, 61]. In our study, with respect to the Argens River basin, the lower parasite diversity and lower parameters of metazoan parasite infection found in introduced BB when compared to the local BM (Fig. 4) indicate that parasites are adapted to their local BM populations. On the other hand, elevated metazoan parasite infection in native BM might be a result of parasite transmission from introduced BB. Out of all the parasite species that infected both the Rhône late sympatric BB and the Argens recent sympatric BB individuals, and might, therefore, have been carried along with host introductions, almost 70% were also found in BM (Table 2). To our knowledge, and with the exception of this study, no investigation of metazoan parasite communities in fish from the Argens River system has so far been undertaken. We, therefore, cannot exclude the possibility that the parasite species that we reported from both BB and BM were already present in the Argens River basin before the introduction of BB. Myxobolus spp. parasites as well as endoparasites such as Pseudocapillaria tomentosa, Caryophyllaeus brachycollis and Pomphorhynchus terreticollis, or parasites belonging to the genus Apharyngostrigea, infect a wide range of fish hosts beside barbels [62–65]. On the other hand, Bathybothrium rectangulum and Rhabdochona hellichi, which are commonly found in barbels and rarely found in other fish , were already reported from both BB and BM [62, 65]. Similarly, Gyrodactylus katharineri, G. hemibarbi and G. markewitchi found in both Barbus species in our study, have already been documented from other cyprinids besides barbels, such as Cyprinus carpio, Alburnus alburnus, Gobio gobio, Leuciscus cephalus and Gymnocephalus cornua [66–68]; these cyprinids were already reported from the Argens River basin . Finally, in our study, the monogenean parasite Paradiplozoon homoion occurred in very low abundance and prevalence on both BB and BM. Le Brun et al.  documented that Paradiplozoon gracile never parasitized BB, while it is present on BM and H individuals of these two species (P. gracile was proposed as a synonym of P. homoion ). Slightly higher values of P. homoion infection in BM than in BB of the Argens River basin and its absence in the Loire allopatric BB population may even indicate the opposite direction of transmission, i.e. from BM to BB. However, this parasite species was already reported in BB and exhibits a wide range of cyprinid hosts [43, 62, 63, 70]. Even if all these parasites that were, in our study, shared between barbel hosts were indeed present in the Argens River basin, contact zones between these two Barbus species may act as a reservoir for the reinforcement of infection of local BM from introduced BB. On the other hand, the exposure of susceptible non-native host species to native parasites might lead to parasite spillback, i.e. an increase in the impact of a disease on a local host via transmission from a competent alien host . This may be supported by higher similarity in parasite communities between the Argens recent sympatric BB and the Argens recent sympatric BM than between the Argens recent sympatric BB and the Argens allopatric BM. In general, similar parasite communities are expected in hosts living in sympatric areas as a result of host contact, which facilitates parasite transmission . Our findings, therefore, highlight the risk of possible parasite transmission from invasive BB to endemic BM. Šimková et al.  indicated that local and endemic Parachondrostoma toxostoma became infected by Dactylogyrus parasites after contact with the invasive Chondrostoma nasus in southern France. In our study, Dactylogyrus parasites were very rare in BM. Dactylogyrus malleus, the only representative of the genus which followed the expansion of BB, was not transmitted to BM, even though H individuals are susceptible to infection with this species (Table 2). The higher infection levels of these parasites, as was found in the Loire allopatric BB and the Rhône late sympatric BB, and subsequently their capacity to be transmitted to local BM might, therefore, be expected in the very near future.
A continuum of hybrid genotypes and phenotypes between parental species may fill the gap between the two hybridizing species and facilitate parasite transmission . As a result, we might expect that mixed genotypes arising by introgressive hybridization should enable more opportunities for infection by parasite species which naturally track parental species, i.e. hybrids will be infected by parasites of both parental taxa [24, 25]. In our study, hybrids of the Argens River basin were, overall, parasitized by a number of parasite species intermediate between pure species (Fig. 4). However, in localities where at least one parental species was present, hybrid individuals tended to be parasitized by more parasite species (Tables 1 and 2). Moreover, the effect of host on the number of parasite species which infected Barbus populations was confirmed by the results of GLM analysis and may be a result of differences in host ethology and ecological preference . Thus, our findings suggest that hybrids represent “bridges” for parasite infection between invasive and endemic species. From this point of view, it seems that BB represents a potential threat to local BM in terms of the transmission or increased impact of metazoan parasites on local BM via hybridization. In fact, all parasites that infected both parental taxa were also detected in hybrids of the Argens River basin. Consequently, we may conclude that metazoan parasites represent important biomarkers of BB and BM hybridization in the Argens River basin, as was also reported for Dactylogyrus parasites in Alburnus alburnus × Rutilus rutilus hybrids (Lake Mikri Prespa, Greece; ). Šimková et al.  reported that Cyprinus carpio × Carassius gibelio hybrids were also parasitized by a greater variety of parasite species than pure hosts; however, they remained less susceptible to metazoan parasites, which could be a result of so-called hybrid vigour . In our study, hybrids displayed higher fitness overall than BM in terms of lower parasite abundance and resembled BB more in this respect (Fig. 4). Slightly greater similarity in parasite communities was also revealed between H and invasive BB than between H and endemic BM, based on presence/absence data. Phillipart & Berrebi  showed that experimental crossing between female BB and male BM resulted in F1 hybrids with a similar size structure to BB and favored female F1 hybrids. As a result, the maternal effect may influence offspring (in our case, hybrid) susceptibility to parasites . Yet, hybrid individuals displayed different levels of susceptibility to metazoan parasites (i.e. higher or lower abundance), in contrast to the parental taxa in localities where parental species were also present (Table 1). These findings are consistent with the results of GLM, which showed no effect of host on metazoan parasite abundance, while significant effects of host body size, locality, sampling year, and water temperature on metazoan parasite abundance were demonstrated. Since parasite infection in hosts is unstable in the time and space driven by host-parasite interactions and environmental forces [76–78], our findings indicate that both individual host characteristics (i.e. host genotype) and environmental factors together significantly influence spatio-temporal distribution of metazoan parasite communities in the BB × BM hybrid system. Since metazoan parasites with the complex life-cycle constitute a substantial fraction of metazoan parasite communities in Barbus spp. of the Argens River basin, spatial differences in diversity of parasite communities found in our study might be also partially shaped by the availability and abundance of intermediate hosts, i.e. the presence of Nematode sp. 3 in the Argens allopatric BM. The effect of intermediate host abundance on the parasite diversity was previously shown in eels .
On the basis of our results, we may conclude that metazoan parasites extend along the expansion range of invasive B. barbus. While similar levels of metazoan parasite infection were revealed in B. barbus of the Loire River basin (i.e. in the absence of B. meridionalis), and B. barbus of the Rhône River basin (i.e. where there is a low level of microsatellite introgression from B. meridionalis), B. barbus recently introduced into the Argens River was shown to profit from enemy release after its arrival from the River Durance by displaying lower susceptibility to metazoan parasites in contrast to the source populations. Concerning the Argens River basin, lower levels of parasite infection in populations of B. barbus in comparison to B. meridionalis and a similarity in metazoan parasite communities with those found in the Rhône River basin support the idea that the introduction of B. barbus into the watersheds of the River Argens from the River Durance is of very recent origin. The infection of hybrids by metazoan parasites found in both parental species in the River Argens supports the existence of hybridization between B. barbus and B. meridionalis and indicates that parasites along with molecular markers may be used as powerful tools for detecting recent hybridization events. The transmission of parasites via introgressive hybridization and higher parasite infection in B. barbus from the Loire River basin and the Rhône River basin indicated in our study may highlight the potential risk of non-native B. barbus having an increased disease impact on endangered B. meridionalis.
General linear model
Qiu YX, Hong DY, Fu CHX, Cameron KM. Genetic variation in the endangered and endemic species Changium smyrnioides (Apiaceae). Biochem Sys Ecol. 2004;32:583–96.
Juan A, Crespo MB, Cowan RS, Lexer C, Fay MF. Patterns of variability and gene flow in Medicago citrina, an endangered endemic of islands in the western Mediterranean, as revealed by amplified fragment length polymorphism (AFLP). Mol Ecol. 2004;13:2679–90.
Ogutu-Ohwayo R. The effects of predation by nile perch, Lates niloticus L., on the fish of Lake Nabugabo, with suggestions for conservation of endangered endemic cichlids. Conserv Biol. 1993;7:701–11.
Ogutu-Ohwayo R. Nile perch in Lake Victoria: the balance between benefits and negative impacts of aliens. In: Sanlund OT, Schei PJ, editors. Invasive species and biodiversity management. Dordrecht: Kluwer Academic Publishers; 1999. p. 47–63.
Fisher RN, Bradley SH. The decline of amphibians in California's Great Central Valley. Conserv Biol. 1996;1:1387–97.
Vilà M, Weiner J. Are invasive plant species better competitors than native plant species? Evidence from pair-wise experiments. Oikos. 2004;105:229–38.
Perry WL, Lodge DM, Feder JL. Importance of hybridization between indigenous and nonindigenous freshwater species: an overlooked threat to North American biodiversity. Syst Biol. 2002;51:255–75.
De la Rúa P, Jaffé R, Dall'Olio R, Muñoz I, Serrano J. Biodiversity, conservation and current threats to European honeybees. Apidologie. 2009;40:263–84.
Biedrzycka A, Solarz W, Okarma H. Hybridization between native and introduced species of deer in Eastern Europe. J Mammal. 2012;93:1331–41.
Knopf K, Mahnke M. Differences in susceptibility of the European eel (Anguilla anguilla) and the Japanese eel (Anguilla japonica) to the swimbladder nematode Anguillicola crassus. Parasitology. 2004;129:491–6.
Woodworth BL, Atkinson CT, LaPointe DA, Hart PJ, Spiegel CS, Tweed EJ, et al. Host population persistence in the face of introduced vector-borne diseases: Hawaii amakihi and avian malaria. Proc Natl Acad Sci U S A. 2005;102:1531–6.
Kirk RS. The impact of Anguillicola crassus on European eels. Fish Manag Ecol. 2003;10:385–94.
Egusa S. Nematode diseases. In: Egusa S, editor. Infectious Diseases of Fish. Rotterdam/Brookfield: A.A. Balkema; 1992. p. 643–57.
van Banning P, Haenen OLM. Effects of the swimbladder nematode Anguillicola crassus in wild and farmed eel, Anguilla anguilla. In: Perkins FO, Cheng TC, editors. Pathology in Marine Science, New York: Academic Press, vol. 1990. 1990. p. 317–30.
Molnár K, Székely C, Baska F. Mass mortality of eel in Lake Balaton due to Anguillicola crassus infection. Bull Eur Assoc Fish Pathol. 1991;11:211–2.
Molnár K, Baska F, Csaba G, Glávits R, Székely C. Pathological and histopathological studies of the swimbladder of eels Anguilla anguilla infected by Anguillicola crassus (Nematoda: Dracunculoidea). Dis Aquat Organ. 1993;15:41–50.
Baruš V, Moravec F, Prokeš M. Anguillicolosis of the European eel (Anguilla anguilla) in the Czech Republic. Czech J Anim Sci. 1999;44:423–31.
Torchin ME, Lafferty KD, Dobson AP, McKenzie VJ, Kuris AM. Introduced species and their missing parasites. Nature. 2003;421:628–30.
Keane RM, Crawley MJ. Exotic plant invasions and the enemy release hypothesis. Trends Ecology Evol. 2002;27:164–70.
Blakeslee AMH, Altman I, Whitman Miller A, Byers JE, Hamer CE, Ruiz GM. Parasites and invasions: a biogeographic examination of parasites and hosts in native and introduced ranges. J Biogeogr. 2012;39:609–22.
Torchin ME, Lafferty KD, Kuris AM. Release from parasites as natural enemies: increased performance of a globally introduced marine crab. Biol Invasions. 2001;3:333–45.
Ondračková M, Valová Z, Hudcová I, Michálková V, Šimková A, Borcherding J, Jurajda P. Temporal effects on host-parasite associations in four naturalized goby species living in sympatry. Hydrobiologia. 2015;746:233–43.
Poulin R, Mouillot D. Host introductions and the geography of parasite taxonomic diversity. J Biogeogr. 2003;30:837–45.
Dupont F, Crivelli AJ. Do parasites confer a disadvantage to hybrids? A case study of Alburnus alburnus × Rutilus rubilio, a natural hybrid of Lake Mikri Prespa, Northern Greece. Oecologia. 1988;75:587–92.
Šimková A, Dávidová M, Papoušek I, Vetešník L. Does interspecies hybridization affect the host specificity of parasites in cyprinid fish? Parasit Vectors. 2013;6:95.
Floate KD, Whitham TG. The "hybrid bridge" hypothesis: host shifting via plant hybrid swarms. Am Nat. 1993;141:651–62.
Persat H, Berrebi P. Relative ages of present populations of Barbus barbus and Barbus meridionalis (Cyprinidae) in southern France: preliminary considerations. Aquat Living Resour. 1990;3:253–63.
Kotlík P, Berrebi P. Phylogeography of the barbel (Barbus barbus) assessed by mitochondrial DNA variation. Mol Ecol. 2001;10:2177–86.
Berrebi P, Lamy G, Cattaneo-Berrebi G, Renno JF. Variabilité génétique de Barbus meridionalis Risso (Cyprinidae): une espèce quasi monomorphe. B Fr Peche Piscic. 1988;310:77–84.
Crespin L, Berrebi P. L’hybridation naturelle entre le barbeau commun et le barbeau méridional en France. Compte rendu de dix années de recherche. B Fr Peche Piscic. 1994;334:177–89.
Agence de l'eau Rhône Méditerranée et Corse. Etude de détermination des volumes prélevables. Bassin versant de l’Argens. Phases 4, 5 et 6. 2013. http://www.rhone-mediterranee.eaufrance.fr/docs/gestion-quantitative/EEVPG/argens/EVP_argens_rapport_ph4%265%266_juin2013.pdf. Accessed 27 Feb 2016.
BIOTOPE – Agence Paca. Document d'objectifs du site Natura FR9301627 « Embouchure de L´Argens » Tome 1, Fiches Espèces. 2011. http://embouchure-argens.n2000.fr/sites/embouchure-argens.n2000.fr/files/documents/page/DOCOB20-20Fiches20especes.pdf. Accessed 01 Mar 2016.
Kiener A, Delize A, Belkior P. Aspects piscicoles du flevue Argens (Var). Bull Fr Piscic. 1981;282:43–74.
Molecular Ecology Resources Database. http://tomato.bio.trinity.edu/manuscripts/13-5/mer-13-0090.pdf. Accessed 27 Feb 2016.
Agostini C, Albaladejo RG, Aparicio A, Arthofer W, Berrebi P, Boag PT, et al. Permanent genetic resources added to molecular ecology resources database 1 April 2013–31 May 2013. Mol Ecol Resour. 2013;13:966–8.
Chenuil A, Desmarais E, Pouyaud L, Berrebi B. Does polyploidy lead to fewer and shorter microsatellites in Barbus (Teleostei: Cyprinidae)? Mol Ecol. 1997;6:169–78.
Van Oosterhout C, Hutchinson WF, Derek P, Wills M, Shipley P. MICRO-CHECKER: software for identifying and correcting genotyping errors in microsatellite data. Mol Ecol Notes. 2004;4:535–9.
Pritchard JK, Stephens M, Donnelly P. Inference of population structure using multilocus genotype data. Genetics. 2000;155:945–59.
Gompert Z, Buerkle CA. Introgress: a software package for mapping components of isolation in hybrids. Mol Ecol Resour. 2010;10:378–84.
Andrés JA, Larson EL, Bogdanowicz SM, Harrison RG. Patterns of transcriptome divergence in the male accessory gland of two closely related species of field crickets. Genetics. 2013;193:501–13.
Ergens R, Lom J. Causative agents of fish diseases. Prague: Academia; 1970 (in Czech).
Lamková K, Šimková A, Palíková M, Jurajda P, Lojek A. Seasonal changes of immunocompetence and parasitism in chub (Leuciscus cephalus), a freshwater cyprinid fish. Parasitol Res. 2007;101:775–89.
Pugachev ON, Gerasev PI, Gussev AV, Ergens R, Khotenowsky I. Guide to Monogenoidea of freshwater fish of Palearctic and Amur regios. Milan: Ledizione-Ledi Publishing; 2009.
Špakulová M, Perrot-Minnot MJ, Neuhaus B. Resurrection of Pomphorhynchus tereticollis (Rudolphi, 1809) (Acanthocephala: Pomphorhynchidae) based on new morphological and molecular data. Helminthologia. 2011;48:268–77.
Moravec F. Parasitic nematodes of freshwater fishes of Europe. Prague: Academia; 2013.
Barčák D, Oros M, Hanzelová V, Scholz T. Phenotypic plasticity in Caryophyllaeus brachycollis Janiszewska, 1953 (Cestoda: Caryophyllidea): does fish host play a role? Syst Parasitol. 2014;88:153–66.
Bush AO, Lafferty KD, Lotz JM, Shostack AW. Parasitology meets ecology on its own terms: Margolis et al. revised. J Parasitol. 1997;83:575–83.
Chao A. Species estimation and applications. In: Kotz S, Balakrishnan N, Read CB, Vidakovic B, editors. Encyclopedia of Statistical Sciences, vol. 12. 2nd ed. New York: Wiley; 2005. p. 7907–16.
Colwell RK, Elsensohn JE. EstimateS turns 20: statistical estimation of species richness and shared species from samples, with non-parametric extrapolation. Ecography. 2014;37:609–13.
Hammer Ø, Harper DAT, Ryan PD. PAST: Paleontological statistics software package for education and data analysis. Palaeontol Electronica. 2001;4:1–9.
Ebert D. The effects of Daphnia parasites on host fitness. In: Ebert D, editor. Ecology, epidemiology, and evolution of parasitism in Daphnia. Bethesda: National Center for Biotechnology Information; 2005. p. 41–8.
Johnson PT, Hoverman JT. Parasite diversity and coinfection determine pathogen infection success and host fitness. Proc Natl Acad Sci U S A. 2012;109:9006–11.
Elton CS. The Ecology of invasions by animals and plants. London: Methuen; 1958.
De Bach P. Biological control by natural enemies. London: Cambridge University Press; 1974.
Dobson AP. Restoring island ecosystems: the potential of parasites to control introduced mammals. Conserv Biol. 1998;2:31–9.
Kennedy CR, Bush AO. The relationship between pattern and scale in parasite communities: a stranger in a strange land. Parasitology. 1994;109:187–96.
Rohde K. A critical evaluation of intrinsic and extrinsic factors responsible for niche restriction in parasites. Am Nat. 1979;114:648–71.
Lucas MC, Batley E. Seasonal movements and behaviour of adult barbel Barbus barbus, a riverine cyprinid fish: implications for river management. J Appl Ecol. 1996;33:1345–58.
Gandon S, Michalakis Y. Local adaptation, evolutionary potential and host-parasite coevolution: interactions between migration, mutation, population size and generation time. J Evol Biol. 2002;15:451–62.
Kaltz O, Shykoff JA. Local adaptation in host-parasite systems. Heredity. 1998;81:361–70.
Francová K, Ondračková M. Host-parasite interactions in sympatric and allopatric populations of European bitterling. Parasitol Res. 2011;109:801–8.
Moravec F, Konečný R, Baska F, Rydlo M, Scholz T, Molnár K, Schiemer F. Endohelminth fauna of barbel, Barbus barbus (L.), under ecological conditions of the Danube basin in Central Europe. Prague: Academia; 1997.
Moravec F. Checklist of metazoan parasites of fishes of the Czech Republic and the Slovak Republic (1873–2000). Prague: Academia; 2001.
Molnár K, Eszterbauer E, Marton S, Székely C, Eiras JC. Comparison of the Myxobolus fauna of common barbel from Hungary and Iberian barbel from Portugal. Dis Aquat Organ. 2012;100:231–48.
Guttiérrez-Galindo JF, Lacasa-Millán MI, Castellá-Espuny J, Muñoz-López E. Helminths of Barbus meridionalis meridionalis Risso, 1826 in northeastern Spain. Acta Parasitol. 1995;40:140–1.
Ergens R. A survey of the results of studies on Gyrodactylus katharineri Malmberg, 1964 (Gyrodactylidae: Monogenea). Folia Parasitol. 1983;30:319–27.
Gelnar M, Sebelova S, Dusek L, Koubkova B, Jurajda P, Zahradkova S. Biodiversity of parasites in freshwater environment in relation to pollution. Parassitologia. 1997;39:189–99.
Angelescu N. Contribution to the knowledge of the parasitofauna of the Danube fishes from the reservoir of Iron Gates. Trav Mus Nat His Nat Gr Antipa. 1974;XV:13–22.
Le Brun N, Renaud F, Berrebi P, Lambert A. Hybrid zones and host-parasite relationships: effect on the evolution of parasitic specificity. Evolution. 1992;46:56–61.
Matejusová I, Koubková B, Gelnar M, Cunningham CO. Paradiplozoon homoion Bychowsky & Nagibina, 1959 versus P. gracile Reichenbach-Klinke, 1961 (Monogenea): two species or phenotypic plasticity? Syst Parasitol. 2002;53:39–47.
Mastisky SE, Veres JK. Field evidence for parasite spillback caused by exotic mollusc Dreissena polymorpha in an invaded lake. Parasitol Res. 2010;106:667–75.
Šimková A, Navrátilová P, Dávidová M, Ondračková M, Sinama M, Chappaz R, Gilles A, Costedoat C. Does invasive Chondrostoma nasus shift the parasite community structure of endemic Parachondrostoma toxostoma in sympatric zones? Parasit Vectors. 2012;5:200.
Moulia C, Le Brun N, Loubes C, Marin R, Renaud F. Hybrid vigor in parasites of interspecific crosses between two mice species. Heredity. 1995;74:48–52.
Philippart JC, Berrebi P. Experimental hybridization of Barbus barbus and Barbus meridionalis: physiological, morphological, and genetic aspects. Aquat Living Resour. 1990;3:325–32.
Boots M, Roberts KE. Maternal effects in disease resistance: poor maternal environment increases offspring resistance to an insect virus. Proc R Soc Lond B Biol Sci. 2012;279:4009–14.
Mitchell SE, Rogers ES, Little TJ, Read AF. Host‐parasite and genotype‐by‐environment interactions: temperature modifies potential for selection by a sterilizing pathogen. Evolution. 2005;59:70–80.
Grech K, Watt K, Read AF. Host-parasite interactions for virulence and resistance in a malaria model system. J Evol Biol. 2006;19:1620–30.
Wolinska J, Keller B, Manca M, Spaak P. Parasite survey of a Daphnia hybrid complex: host-specificity and environment determine infection. J Anim Ecol. 2007;76:191–200.
Sures B, Streit B. Eel parasite diversity and intermediate host abundance in the River Rhine, Germany. Parasitology. 2001;123:185–91.
This study was funded by GACR project No. P505/12/0375. We thank Radim Blažek, Rémi Chappaz, Jean-Pierre Cornec, Caroline Costadoat, Martina Dávidová, Kateřina Francová, Milan Gelnar, Lucka Havlátová, Šárka Mašová, Petra Navrátilová, Markéta Ondračková, Jarek Piňos, Eva Řehulková, Tomáš Scholz, and the representatives of Maison Régionale de l’Eau, France for help in the field and with parasite determination, and Kristína Civáňová and Kristýna Koukalová for their assistance in the laboratory. We kindly thank Matthew Nicholls for English revision of the final draft.
This study was funded by the Czech Science Foundation, Project No. P505/12/0375.
Availability of data and material
All datasets supporting the main conclusions are included within the manuscript and its additional files. Individual level data are available from the authors on reasonable request.
AŠ designed and supervised the study; AŠ, AG and LG conducted the field study; LG performed the data analysis; LG wrote the manuscript; AŠ critically revised the draft. All authors read and approved the final manuscript.
The authors declare that they have no competing interests.
Ethics approval and consent to participate
This study was approved by animal care and use committee of Masaryk University in Brno, Czech Republic (approval number CZ01308).
An erratum to this article is available at http://dx.doi.org/10.1186/s13071-017-2214-5.
Site names and their geographical coordinates. (XLSX 10 kb)
Metazoan parasite abundance, prevalence and species richness in Barbus spp. populations. (XLSX 11 kb)
Results of ANCOVA tests. (XLSX 16 kb)
P-values of Bonferroni post-hoc tests. (XLSX 13 kb)
Metazoan parasite communities in Barbus spp. populations. (XLSX 21 kb)
About this article
Cite this article
Gettová, L., Gilles, A. & Šimková, A. Metazoan parasite communities: support for the biological invasion of Barbus barbus and its hybridization with the endemic Barbus meridionalis . Parasites Vectors 9, 588 (2016). https://doi.org/10.1186/s13071-016-1867-9