# Equine strongyle communities are constrained by horse sex and species dipersal-fecundity trade-off

Guillaume Salle, Slawomir Kornaś, Marta Basiaga

► To cite this version:

Guillaume Salle, Slawomir Kornaś, Marta Basiaga. Equine strongyle communities are constrained by horse sex and species dipersal-fecundity trade-off. Parasites & Vectors, 2018, 11 (1), pp.1-11. 10.1186/s13071-018-2858-9 . hal-02624249

HAL Id: hal-02624249

https://hal.inrae.fr/hal-02624249

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>
RESEARCH

Open Access

<figure><img src="image_2.png" /></figure>

CrossMark

# Equine strongyle communities are constrained by horse sex and species dipersal-fecundity trade-off

Guillaume Sallé¹*, Slawomir Kornaś² and Marta Basiaga²

## Abstract

**Background:** Equine strongyles are a major health issue. Large strongyles can cause death of horses while cyathostomins (small strongyles) have shown increased resistance to anthelmintics worldwide. Description of strongyle communities have accumulated but little is known about the diversity of these communities and underpinning environmental factors.

**Methods:** Strongyles were recovered after ivermectin treatment from 48 horses located in six premises in Poland. Correlation between previously published species fecundity and the observed relative abundance and prevalence were estimated. Significance of horse sex was determined at the species level (prevalence, relative abundance) and at the community level (species richness and dissimilarity between communities).

**Results:** Strongyle species fell into two groups, contrasted by their prevalence and relative abundance. Six to nine horses were necessary to sample at least 90% of strongyle community diversity, providing a minimal cut-off to implement sampling trial in the field. Strongyle communities entertained a network of mostly positive interactions and species co-occurrence was found more often than expected by chance. In addition, species fecundity and prevalence were negatively correlated (Pearson's *r* = -0.71), suggesting functional trade-offs between species dispersal abilities and fecundity. This functional trade-off may underpin species coexistence. Horse sex was also a significant constraint shaping strongyle communities. Indeed, mares generally displayed more similar strongyle communities than stallions (*P* = 0.003) and *Cylicostephanus calicatus* was more abundant in stallions suggesting sex-specific interactions (*P* = 0.006).

**Conclusions:** While niche partitioning is likely to explain some of the positive interactions between equine strongyle species, coexistence may also result from a functional trade-off between dispersal ability and fecundity. There is significant evidence that horse sex drives strongylid community structure, which may require differential control strategies between mares and stallions.

**Keywords:** Horse, Cyathostomes, Strongyles, Ecology, Sex, Nematode, Diversity

## Background

Grazing horses are infected by a wide variety of intestinal strongyle parasites belonging to the subfamilies Strongylinae and Cyathostominae [1]. Among the Strongylinae, *Strongylus vulgaris* has been identified as a major cause of colic in horses [2]. This results from the accumulation of larval stages in the cranial mesenteric artery that leads to an arteritis responsible for fatal intestinal infarction [2].

However, modern anthelmintic treatments have dramatically decreased its prevalence [3]. Cyathostominae (cyathostomes or small strongyles) form a wider tribe than equine Strongylinae, encompassing 51 species, 40 of which infect horses [1]. Unlike the Strongylinae, these nematodes can reach prevalences up to 100% in horses [4]. Cyathostomes have been associated with milder clinical signs like poor hair coat and weight loss [5], although the massive emergence of larval stages from the large intestine mucosa where they reside can lead to larval cyathostominosis cases, characterised by diarrhoea and protein-losing enteropathy [6]. Failure to control larval

<footnote>* Correspondence: Guillaume.Salle@inra.fr
¹UMR1282 INRA/Université de Tours Infectiologie et Santé Publique, F-37380 Nouzilly, France
Full list of author information is available at the end of the article</footnote>
cyathostominosis can also lead to the death of affected horses in at least a third of cases, as suggested by a report on 15 clinical cases in the UK [7]. More importantly, cyathostomes have been involved in multiple reports of anthelmintic failures throughout the world thus representing a major issue in equine medicine [8–11].

Despite the wide variety of species, cyathostome communities are consistently described as a core assemblage of 10–13 species [4, 12–16]. While this has been frequently reported, only scant knowledge has been gathered about the drivers underpinning strongyle biodiversity or community assemblage. Nonetheless, some efforts have been made to characterise the consequences of deworming treatments on strongyle diversity as an attempt to better understand widespread anthelmintic failures. First, a reduced diversity in strongyle communities was observed in Ukrainian farms applying frequent deworming [17]. Secondly, investigations of strongyle community structure in horse populations after anthelmintic treatment found that the first species to reappear represented a limited subset of total community [18, 19]. *Cylicocyclus nassatus* and *Cylicostephanus longibursatus* were consistently found among the major species after a combination of pyrantel followed by fenbendazole treatment [18] or after ivermectin treatment [19]. Another study also suggested that these two species were driving the drug resistance phenotype [20, 21]. However, these observations might be confounded by the reported variation in the number of eggs found in *utero* across strongyle species [22], dispersal ability (measured by the prevalence of infection) or relative abundance [12, 13, 16]. Indeed, these latter properties may affect the post-treatment community structure, e.g. the most abundant species contributing most to the post-treatment community, without any genuine association with anthelmintic susceptibility.

Some insights into the structure of equine strongyle communities have been provided by a report focusing on pair-wise correlations between strongyle species abundance [23]. This analysis revealed negative interactions between major species, suggesting competition might be involved [23]. In addition, an attempt to bring together environmental factors with equine strongyle diversity in Ukraine concluded that foals and horses older than 16 years displayed higher species richness [24]. However, species richness is only a limited facet of the actual biodiversity as it does not account for the relative proportions of every species found. Furthermore, the impact of host sex was not considered [24].

Host sex effect is a major driver of helminth community assemblages, male mammals being generally more frequently and more heavily parasitized [25, 26], although some discrepancies exist to this general pattern [27]. In horses, no clear consensus has been reached so far. The

analysis of necropsy reports in Australia found differential prevalence and abundance of 13 endoparasites between the two sexes [28]. However, no general trend emerged as females were more heavily infected by *Cyathostomum pateratum* but less frequently infected by *Cylicodonto-phorus bicoronatus*, *Cylicocyclus insigne* and *Cylicocyclus elongatus* than stallions and geldings [28]. Study of faecal egg count in a feral horse population reported a spatially contrasted polarity of the sex effect, stallions being more affected than females on one side of the island while the opposite trend was found on the other side [29]. Other studies reported higher faecal egg count in geldings or could not evidence any difference [30, 31].

In this study, we analysed strongyle communities recovered from 48 horses in Poland to investigate species community structure and diversity, the link between previously published fecundity estimates and species dispersal ability and to assess how they were affected by horse sex.

## Methods

### Study population

Strongyles were collected in spring 2011 and 2012 from 48 horses scattered across six premises from Poland (Additional file 1: Table S1). Premises had been enrolled for an ivermectin efficacy test. Horses were selected based on their faecal egg count before treatment to allow for strongyle recovery. The study population was composed of 37 mares and 11 stallions. The Małopolska breed accounted for most of the sampled horses (n = 25), in addition to some Hucul (n = 6) and pure blood Arabian horses (n = 17) that were found in only two sites (Additional file 1: Table S1). Two thirds of the individuals were born before 2007 (Additional file 1: Table S1).

Within each premise, the number of sampled horses ranged from 3 to 18 individuals (Table 1, Additional file 1: Table S1).

Parasite control strategy in every location relied on ivermectin treatments administered in March (or April) and October, before and after pasture season, respectively.

### Strongyle collection and identification procedure

The species composition of strongylid nematodes was performed by collecting 500 g of faecal material expelled 24 h after ivermectin treatment following the strategy of other trials [14, 32]. Subsequently, samples were rinsed with tap water on 250 µm mesh sieves and examined under stereomicroscope for the presence of strongylids. Collected nematodes were fixed in 70% ethyl alcohol with 5% glycerine additive to clear their body structure. Strongylids were identified according to previously published morphological keys [1].

Faecal egg count before treatment was performed following a modified McMaster technique [33] on 5 g of
<table><caption>**Table 1** Premise features and ecological properties of recovered strongyle communities. Average species richness (number of recovered strongyle species) and Shannon's index of diversity are provided for each premise</caption><thead><tr><th>Premise</th><th>Type</th><th>No. of horses (n)</th><th>Mares</th><th>Stallions</th><th>Species richness ± SD</th><th>Shannon index ± SD</th></tr></thead><tbody><tr><td>1</td><td>Stud farm</td><td>18</td><td>18</td><td>0</td><td>12.17 ± 2.04</td><td>2.07 ± 0.17</td></tr><tr><td>2</td><td>Riding school</td><td>14</td><td>8</td><td>4</td><td>9.86 ± 3.74</td><td>1.59 ± 0.48</td></tr><tr><td>3</td><td>Riding school</td><td>3</td><td>2</td><td>1</td><td>6.00 ± 0.01</td><td>1.46 ± 0.23</td></tr><tr><td>4</td><td>Riding school</td><td>3</td><td>1</td><td>2</td><td>6.67 ± 1.53</td><td>1.43 ± 0.12</td></tr><tr><td>5</td><td>Riding school</td><td>3</td><td>0</td><td>3</td><td>11.67 ± 6.51</td><td>1.75 ± 0.38</td></tr><tr><td>6</td><td>Riding school</td><td>7</td><td>6</td><td>1</td><td>11.00 ± 2.16</td><td>1.90 ± 0.14</td></tr><tr><td>Total</td><td></td><td>48</td><td>35</td><td>11</td><td>10.6 ± 3.4</td><td>1.8 ± 0.38</td></tr></tbody></table>

Abbreviation: SD standard deviation

faecal material, diluted in 70 ml of a flotation solution (NaCl, d = 1.18), with a sensitivity of 50 eggs/g.

## Community structure analyses

Statistical analyses were implemented with R v 3.4.1 [34]. Unless stated otherwise, any P-value below 5% was regarded as significant. Strongyles recovered from each horse are hereafter referred to as a community. Within each community, species intensity refers to the number of worms recovered from a given horse, while species relative abundance in a given community corresponds to the proportion of worms found relative to the total worm burden. Species prevalence was defined at the premise level as the proportion of horses infected by a given strongyle species.

Species fecundity may contribute to observed relative abundance and prevalence. To investigate this point, the correlation between previously reported number of eggs *in utero* recovered from strongyle females [22] and observed species abundance and prevalence were estimated (*rcorr* function of the *Hmisc* package v.4.0-3 [35]). Species fecundity estimates were log-transformed to obtain a linear relationship between variables.

Our sampling procedure involved a limited amount of faecal material which may have biased our prevalence and abundance estimates. Therefore, our observations were compared against a dataset [36] of previously published equine strongyle species prevalence (n = 29 papers) and relative abundances (n = 14 papers) provided as Additional file 2: Table S2. These studies spanned all continents and involved various methodologies (necropsy, deworming), thus limiting confounding factors linked to climatic conditions or sampling methods. For each paper, species relative abundance and prevalence were computed from the data if not already provided. These individual estimates were subsequently averaged across studies for comparison with the observed values from our 48 horses.

Community diversity analyses were run using the *vegan* package v. 2.4-3 [37]. Species richness was defined

as the total number of species recovered from each horse. Species diversity within each horse, corresponding to Whittaker's α-diversity [38], was measured by the Shannon's index [39] defined as:

$$H = \sum_{i}^{n} p_{i} \log(p_{i})$$

where $p_i$ is the proportional abundance of species i. This was implemented with the *diversity* function of the *vegan* package [37]. Normality of species richness was tested by a Shapiro-Wilk test.

Species turnover can occur along environmental gradients, hence resulting in variation among individual communities from a given group. This measure of dissimilarities between communities was originally coined as β-diversity [38], and will contribute only marginally to the global diversity if individual communities show little or no variation among groups. Difference in β-diversity among environmental factors (location, horse age and horse sex) was tested following a multivariate framework [40, 41], implemented in the *adonis* function of the *vegan* package [37]. This framework measures β-diversity as the distance of every community to the average group community, determined as the centroid of a principal components analysis of community structures [40].

Non-random species co-occurrence was tested following a probabilistic framework that compares observed species frequencies to the expected frequencies under species independence (discarding species pairs when the expected number of co-occurrence is less than 1) [42, 43]. Species interactions could also result in positive or negative feedbacks that could impact on species relative abundances. To test for these interactions, species counts were regressed upon other species intensity in a pairwise manner after accounting for environmental variables (horse sex and sampling site) and by means of a Poisson regression. For this analysis, the 13 species with more than 200 individuals across horses were considered, as failure to do
so results in strong positive associations for the minor species present. Because of the high number of statistical tests (n = 105), a Bonferroni correction was applied and a P-value below 0.05/105 = 5 × 10⁻⁴ was considered significant. Significant pairwise regression estimates (P < 5 × 10⁻⁴) were used to generate a species network, visualized with the ggnet package [44]. The network uses a force-directed algorithm from the matrix of interactions. Nodes act as repelling objects that are organized to minimize forces in the network.

## Host sex effect on species relative abundance and prevalence

The impact of host sex on strongyle abundance (worm counts) and prevalence was investigated for the 11 core strongyle species, i.e. most abundant and most prevalent. Cyathostomum catinatum showed the highest prevalence (47 horses infected) and was not included in the prevalence data analysis as not enough variability was present to assess factors of variation. To account for the data structure (Additional file 1: Table S1), these effects were investigated using only the data from premises where horses of each gender were present (premises 2, 3, 4 and 6; n = 27 horses). Worm counts and species prevalence were modelled using a negative binomial and logistic regression, respectively.

Samples were recovered in 2011 and 2012, hence introducing a putative confounding annual variation. This was tested in two ways. First, Pearson's correlation between prevalence and relative abundance measured in 2011 and 2012 were estimated to support their congruence. Secondly, worm counts and species prevalence were regressed on a universal year effect shared across species and a species-specific inter-annual variation.

After ruling out a significant inter-annual variation, worm counts and species prevalence were modelled fitting the premise, strongyle species, horse sex and their interaction as fixed effects, and horse age as a covariable:

$$y \sim \text{Strongyle species} + \text{Horse age} + \text{Horse sex} + \text{Strongyle species} \times \text{Horse sex} + \text{Premise}$$

where y stands for species count or prevalence accordingly.

A linear regression analysis was applied to Shannon's index values to evaluate variation between premises and horse ages and sexes. For every model, the sex effect on species relative abundance and prevalence was estimated across premises (19 mares and 8 stallions). The premise effect accounted for both management and breed type differences.

## Results

### Strongyle dispersal abilities were contrasted and inversely correlated with their known fecundity estimates

Considered parasite populations were all ivermectin-susceptible as post-treatment FEC were back to 0 in every horse. On average, 608.7 worms per horse were collected, ranging from 111 to 2884 (Additional file 1: Table S1). Overall, 23 species were found across premises (Additional file 1: Table S1), with average species richness ranging from 6 to 12 species (Table 1).

Between 4 and 20 species could be recovered at the horse level (Additional file 1: Table S1). The male to female worm ratio was 0.58 on average (Additional file 3: Table S3). Total worm count was not significantly correlated with pre-treatment faecal egg count (Spearman's r(48) = 0.22, P = 0.13).

Cyathostomum catinatum and C. biconatus were the most and the least abundant species, respectively; worms of each species summing up to 5462 and one worms recovered across communities, respectively. Strongylus vulgaris was encountered in only 3 horses. Relative species abundance across communities could distinguish between two groups of species, being either highly abundant and prevalent (Fig. 1) or with lower dispersal ability and minor contribution to communities (Fig. 1). Cyathostomum catinatum, C. nassatus and C. longibursatus were the three most abundant species (accounting for 21.1, 17.8 and 13.7% of the total worm burden, respectively) and were found in three quarters of horses (n = 47, 37 and 41 horses, respectively; Fig. 1). These observations showed strong congruency with previously published estimates (Pearson's r(19) = 0.83, P < 0.0001 and r(21)= 0.90, P < 0.0001 for prevalence and relative abundance estimates, respectively). This independent dataset also supported the outlying contribution of C. catinatum, C. nassatus and C. longibursatus, that infected 67.3, 68.9 and 72.1% of horses and accounted for 17.4, 18.2 and 16.8% of total worm burden across studies.

Pearson's correlations between published estimates of worm fecundity and measured species abundance and prevalence were estimated to determine whether fecundity could underpin better dispersal abilities. These correlations suggested that fecundity was negatively correlated with species prevalence (Pearson's r(14) = -0.71, P = 0.004; Fig. 2) and that the same trend held for species relative abundance (Pearson's r(14) = -0.48, P = 0.08; Additional file 4: Figure S1). A significant correlation was found between fecundity and relative abundance when considering only female worm counts (Pearson's r(14) = -0.57, P = 0.03).

A similar pattern was found when using the 29 previously published prevalence estimates (Pearson's
<figure><img src="image_3.png" /><figcaption>**Fig. 1** Relationship between strongyle species prevalence and relative abundance. Mean relative abundance is expressed as the fraction of a given species relative to the total number of strongyles recovered across the 48 horses involved in this study. Mean prevalence represents the percentage of horses infected by a given species. Each dot stands for one species. The black line represents the regression line between the two variables, with associated 95% confidence interval shown as a grey area</figcaption></figure>

$r_{(13)} = -0.86, P < 0.0001$), but the relationship with 14 published species relative abundance was not significant (Pearson's $r_{(10)} = -0.45, P = 0.19$).

Species richness accumulation curves suggested that the minimal sampling effort had to be comprised between six and nine horses to detect 90% of total species present in one location (Fig. 3).

**Communities were structured by a network of mostly positive interactions**

Species co-occurrence analysis was dominated by random co-occurrences that represented 87% of the 200 tested combinations (Fig. 4, Additional file 5: Table S4). Nonetheless, a few species pairs seemed to positively interact with each other ($n = 26$, Fig 4, Additional file 5:

<figure><img src="image_4.png" /><figcaption>**Fig. 2** Relationship between estimated strongyle prevalence and known fecundity estimates. Previously published estimates of strongyle species fecundity (measured as the number of eggs found in female strongyle *in utero*) were plotted against species prevalence across monitored horses (measured as the fraction of horses infected by a given species). Each dot stands for a given strongyle species. The regression line between the two variables is shown as well the 95% confidence interval (grey area)</figcaption></figure>
<figure><img src="image_5.png" /></figure>

<figure><img src="image_6.png" /></figure>

**Fig. 3** Species accumulation curves for the two most populated premises. For these two premises, the sampling effort (number of horses to sample) required to detect a given number of species (y-axis) is plotted as a blue line with 95% confidence interval (blue area). Box-plots stand for random variation measured by permutations. Vertical lines assess the number of horses needed to detect 90% of total species in a given premise

<figure><img src="image_7.png" /></figure>

**Fig. 4** Strongyle species co-occurrence. A pairwise matrix of co-occurrence between species co-occurring more than once. Grey boxes indicate that species pairs occurred randomly, whereas 26 pairs co-occurred significantly more often than expected (blue boxes).

Table S4), including some core species, e.g. *C. longibursatus*-*C. nassatus* (*P* = 0.04) or *Cyathostomum pateratum*-*Cylicostephanus goldi* (*P* = 0.04).

Significant pairwise interactions (Additional file 6: Table S5) were used to build a network between *Coronocyclus*, *Cyathostomum*, *Cylicostephanus* and *Cylicocyclus* species (Fig. 5) that was also mostly driven by positive interactions between species intensities (Fig. 5). Notably, *C. pateratum* was the only species whose intensity varied negatively with other species, except for *C. longibursatus* and *Cylicocyclus lepostomum* (Fig. 5). Three of the best dispersers (*C. catinatum*, *C. nassatus* and *C. longibursatus*) seated on the edge of the network, with weak links to other species, suggesting that variation in their intensity was not impacted by other species. On the contrary, *Coronocyclus* species, i.e. *C. coronatus*, *C. labratus* and *C. labiatus*, were more tightly connected, supporting similar behaviors.

### Horse gender affects dissimilarity between strongyles communities and species abundance

Little inter-annual variation was found as the correlation for prevalence and relative abundances across the two
<figure><img src="image_8.png" /></figure>

years were 0.89 and 0.91, respectively (Pearson's $r_{(23)}$, $P < 0.0001$ in both cases). In addition, neither a universal inter-annual effect shared across species (Wald test $\chi^2 = 1.35$, $df = 1$, $P = 0.25$ for prevalence; Wald test $\chi^2 = 0.28$, $df = 1$, $P = 0.59$ for relative abundance) nor a species-specific annual variation (Wald test $\chi^2 = 14.48$, $df = 10$, $P = 0.15$ for prevalence; Wald test $\chi^2 = 11.01$, $df = 10$, $P = 0.36$ for relative abundance) were found.

Strongyle community diversity ranged from 0.6 to 2.6 at the horse level but was more consistent at the premise level with average values between 1.43–2.07 (Table 1). No significant variation in species α-diversity was found between premises ($F_{(3, 21)} = 2$, $P = 0.15$).

Neither horse age ($F_{(1, 21)} = 2.33$, $P = 0.14$) nor horse sex affected significantly community diversity ($F_{(1, 21)} = 0.33$, $P = 0.57$), suggesting that α-diversity was equivalent between mares and stallions.

An analysis of community dissimilarity was performed to assess the extent of species turnover across horse gender, horse age and premises. Strongyle communities structure was similar across horse ages ($F_{(1, 21)} = 1.60$, $P = 0.24$) but significantly varied between premises ($F_{(3, 21)} = 2.38$, $P = 0.05$) and were also significantly more homogeneous in mares than in stallions ($F_{(1, 21)} = 8.72$, $P = 0.006$).

No systematic horse sex driven differences in strongyle species abundance was identified across strongyle species by the binomial regression approach (Wald test $\chi^2 = 0.8$, $df = 1$, $P = 0.37$, Fig. 5). However, a significant

interaction was found between horse gender and strongyle species abundance (Wald test $\chi^2 = 23.4$, $df = 10$, $P = 0.009$). Notably, *C. calicatus* was significantly more abundant in stallions than in mares ($t_{(271)} = 2.77$, $P = 0.006$), whereas the opposite trend (non-significant) was generally found for other species (Fig. 6, Additional file 7: Table S6). This pattern was also conserved when comparing *C. calicatus* abundance between premise 1 (14.8 worms/horse on average; mares only) and premise 5 (107 worms/horse on average; stallions only) even though this difference could not be disentangled from a premise-associated effect.

Strongyle dispersal ability measured by the prevalence of each species did not vary significantly between mares and stallions, although *C. nassatus* was found slightly more often in mares ($t_{(271)} = -1.94$, $P = 0.05$, Fig. 6, Additional file 7: Table S6).

## Discussion

This study focused on the structure of equine strongyle communities under eastern Europe climatic conditions. The data suggested that reported species fecundity, approximated by the number of eggs *in utero*, was negatively correlated with species prevalence. This relationship was found both in our dataset and across previously published prevalence reports. In addition, horse sex had a statistically significant effect on strongyle community structure and diversity.

Numerous reports on equine strongyle community structure have already evidenced the bimodal pattern
<figure><img src="image_9.png" /><figcaption>**Fig. 6** Strongyle abundance variation according to horse gender. Boxplots of strongyle counts are plotted against horse gender (red for females and blue for males) across premises. Within each box, the midline shows the median, and the lower and upper box edges delimit the 25 and 75% quantiles, respectively. Dots represent outliers (1.5-fold higher than the 75% quantile). Of note, is the higher abundance for *C. calicatus* in males</figcaption></figure>

of their assemblages, with five to ten highly prevalent and abundant species and additional minor species [4, 12–16]. Our data also highlighted the same split between core and satellite species [45], with a set of three successful dispersers, i.e. *C. cati-natum*, *C. nassatus* and *C. longibursatus*, which were found in 75% of horses and represented between 13–21% of the total worm burden. This pattern was also corroborated by estimates from available literature.

Not only corroborating the known literature, this finding also supported the chosen sampling strategy [14, 32]. First, only a limited amount of faecal material was taken from each horse. This could have biased downwards worm counts, eventually resulting in increasing the missing rate of some species. However, it seems that our prevalence and relative abundance fit the pattern observed across previous reports that relied on necropsy or deworming for worm collection. Secondly, a single sampling 24 h after treatment might have biased the sampling efficiency downwards as it has been suggested that a significant proportion of parasite material could still be recovered after 36 and up to 60 h [14]. From our results, some species reported as minor in the present study (*S. edentatus*, *S. vulgaris*, *T. serratus* and *C. elongatus*) may have suffered from the sampling strategy as most of them were expelled 48 h after treatment onwards in a previous trial [14]. Notably, this prevents any firm conclusion regarding the minor contribution of *S. vulgaris* to the strongyle communities, even though displayed prevalence levels were in agreement with other reports from regularly dewormed horses [46, 47].

Nevertheless, community structure and species richness estimates were similar to previous reports under eastern European conditions [14, 17], supporting the implemented sampling method.

Under the considered conditions, a limited number of horses were required to sample 90% of the strongyle community diversity. This is of particular relevance for the implementation of field sampling where necropsy is not an option and minimal interference with premise activity are expected.

Although the overall structure of strongyle communities has been well characterized in equids [4, 12–15], little is known about the mechanisms and variation factors underpinning their stability. Remarkably, sampled communities were entertaining positive interactions, both in terms of presence/absence patterns or in their intensity covariation. This finding was in contradiction with previously reported results from an abattoir survey, which highlighted negative interactions between major species [23]. However, the negative interactions mostly occurred between species that were less represented in our present report, i.e. *S. vulgaris* and *S. edentatus* or *Triodontophorus* species [23]. In addition, our analyses were restricted to the best dispersers (species with highest prevalence); hence preventing spurious positive associations that can occur when correlating low counts [48]. Our results were in line with known positive associations between core species [49]. Reports of helminths species co-occurring more often than expected by chance have also been made for human and rodent systems [50]. These positive interactions could arise from niche partitioning between species within horses. The latter could explain for instance, positive
feedbacks between *C. goldi* and *C. nassatus*, which mostly occupy the dorsal and ventral colon, respect-
ively [16, 23], but does not support positive interac-
tions between *C. minutus* and *C. calicatus*, which
tend to share the same organs [16, 23]. Observed
positive interactions may hence suggest that some of
the equine strongyle species share mutualistic interac-
tions as already reported for rodent models [50].
However, the negative interactions of *C. pateratum*
with other species remains unexplained. It could be
the result of competition between species for space or
nutrients (exploitation competition); it may also sug-
gest interference between species due to local inflam-
mation mediated by one species and affecting others
or cross-reactivity of the host immune system towards
multiple species [51].

The existence of functional trade-offs can result in positive associations between species, as demonstrated in plants [52] or fig wasps [53]. Because of the sparsity of hosts and environmental constraints put on eggs and larvae on pasture, strongyle species with best competing abilities, like higher fecundity, may not be able to colonize every available horse. As a result, horses can be colonized by species with lower competing abilities, e.g. lower fecundity. This ultimately leads to species co-existence despite highly differential competing abilities and an apparent trade-off between species fecundity and prevalence. However, our results may suffer from the proxy used for species fecundity. These estimates are based on the number of eggs found *in utero* in female worms, and were generated for a worm population from the USA [22]. A precise assessment of each species egg production rate by means of barcode sequencing may help to better resolve this point.

The negative correlation between species female rela-
tive abundance and their inferred fecundity may be re-
lated to some crowding effect. It has been described in
sheep trichostrongylids [54] and a negative feedback be-
tween mature worms and incoming larvae has been sug-
gested in horses [55, 56].

Our results suggested that horse sex was not exert-
ing a universal influence on worm relative abundance
across species, but that mares had more consistent
strongyle communities. In addition, higher *C. calica-
tus* worm counts were found in stallions. This finding
was already reported in a previous Australian survey,
which also evidenced additional host sex interactions
for *C. coronatus* and *C. catinatum* [28]. Variation in
the qualitative composition of strongyle structure
might hence contribute to the lack of clear differ-
ences between stallion and mare excretion patterns.
However, mechanistic explanations to support a
higher species turnover among stallions remain
unclear and deserve further exploration.

Conclusions

This survey provided insights into applied and more
fundamental knowledge about equine strongyle commu-
nities and how host sex drives their structure and diver-
sity. The bimodal structure of equine strongyle
communities split into core and satellite species was
confirmed, underscoring the overwhelming contribution
of *C. nassatus*, *C. catinatum* and *C. longibursatus*. We
also found that less than ten horses were needed to sam-
ple most of the diversity, which is particularly relevant
for logistically constrained surveys. More fundamentally,
it seems that equine strongyle communities are largely
structured around a network of positive, if at all, interac-
tions, with the only exception of *C. pateratum*. While
mutualistic interactions and niche partitioning could
explain this pattern, the observed trade-offs between
fecundity and dispersal abilities may also explain this co-
existence. Horse sex did not appear to influence the
overall worm counts, but a higher species turnover was
found among stallions leading to more dissimilar
strongyle communities for this gender. While our data
provided new insights into the structuring equine
strongyle communities, only a limited set of the poten-
tial factors shaping these communities have been consid-
ered. Further research in this area should make use of
recent advances in sequencing technology to increase
the sensitivity of worm sampling and also integrate puta-
tive interactions with the host gut microbiota.

Additional files

**Additional file 1: Table S1.** Individual worm counts and associated metadata relative to each of the 48 horses considered: sex, breed, year of birth, premise, premise type, pre-treatment faecal egg count and worm counts for each nematode species. (XLS 72 kb)

**Additional file 2: Table S2.** List of previously published papers on prevalence and relative abundance for equine strongylids. The list of papers (with digital object identifier and link to PubMed reference) used for comparing observed strongyle species prevalence and relative abundance with previous results obtained across different countries is provided. (XLS 60 kb)

**Additional file 3: Table S3.** Species worm counts distinguishing adult female and adult male worms for every horse. Female and male worm counts are provided for each horse (in rows) and strongyle species (in column). (XLS 79 kb)

**Additional file 4: Figure S1.** Relationship between observed strongyle abundance and published fecundity estimates. Previously published estimates of strongyle species fecundity (measured as the number of eggs found in female strongyles *in utero*) are plotted against species relative abundance (measured as the fraction of the total strongyle community within each horse) across monitored horses. Each dot stands for a given strongyle species. The regression line between the two variables is shown as well the 95% confidence interval (grey area). (TIF 737 kb)

**Additional file 5: Table S4.** Co-occurrence probabilities between species. For each strongyle species pair, the observed and expected co-occurences are reported as well as the probability of co-occurrence. In each case, the probability of negative interaction or positive interaction are provided and significant interactions are highlighted in red. (XLS 49 kb)
**Additional file 6: Table S5.** Pairwise regression coefficients between the 13 most abundant species. To test for species interactions in terms of abundance, species counts were regressed upon other species intensity in a pair-wise manner after accounting for environmental variables (horse sex and sampling site) and by means of a Poisson regression. The pair-wise regression estimates for the 13 species with more than 200 individuals across horses are provided. (XLS 71 kb)

**Additional file 7: Table S6.** Species intensity and prevalence modelling output. Model estimates of the regression of horse age, premise, strongyle species, horse sex and interaction between horse sex and species upon species abundance and prevalence are given. Regression coefficient estimates are provided with associated standard errors and P-values. Within each variable, regression coefficients are expressed as the departure from a reference level (C. *catinatum*, mare, premise 2, for the species, horse sex and premise variable, respectively). Interaction estimates correspond to the deviation of a given species in stallions in comparison to mares. Age was fitted as a covariate. (XLS 68 kb)

### Acknowledgements

The authors are particularly indebted to the premise managers who took part to this survey.

### Availability of data and materials

The datasets supporting the conclusions of this article are included within the article and its additional files.

### Authors' contributions

SK designed the study. SK and MB collected samples and performed cyathostomin species identification. GS analysed the data. GS and SK drafted the manuscript. All authors read and approved the final manuscript.

### Ethics approval and consent to participate

Not applicable.

### Competing interests

The authors declare that they have no competing interests.

### Publisher's Note

Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

### Author details

¹UMR1282 INRA/Université de Tours Infectiologie et Santé Publique, F-37380 Nouzilly, France. ²Department of Environmental Zoology, Institute of Animal Sciences, University of Agriculture in Krakow, al. Mickiewicza 24/28, 30-059, Krakow, Poland.

Received: 10 January 2018 Accepted: 18 April 2018

Published online: 02 May 2018

### References

1. Lichtenfels JR, Kharchenko VA, Dvojnos GM. Illustrated identification keys to strongylid parasites (Strongylidae: Nematoda) of horses, zebras and asses (Equidae). Vet Parasitol. 2008;156:4–161.

2. Duncan JL, Campbell JR. Further observations on the maintenance of a monospecific infection of *Strongylus vulgaris* in the horse. Vet Rec. 1973; 92:533.

3. Nielsen MK, Vidyashankar AN, Olsen SN, Monrad J, Thamsborg SM. *Strongylus vulgaris* associated with usage of selective therapy on Danish horse farms-is it reemerging? Vet Parasitol. 2012;189:260–6.

4. Lyons ET, Tolliver SC, Drudge JH. Historical perspective of cyathostomes: prevalence, treatment and control programs. Vet Parasitol. 1999;85:97–112.

5. Murphy D, Love S. The pathogenic effects of experimental cyathostome infections in ponies. Vet Parasitol. 1997;70:99–110.

6. Love S, Murphy D, Mellor D. Pathogenicity of cyathostome infection. Vet Parasitol. 1999;85:113–22.

7. Giles CJ, Urquhart KA, Longstaffe JA. Larval cyathostomiasis (immature trichonema-induced enteropathy): a report of 15 clinical cases. Equine Vet J. 1985;17:196–201.

8. Canever RJ, Braga PR, Boeckh A, Grycajuck M, Bier D, Molento MB. Lack of *Cyathostomin* sp. reduction after anthelmintic treatment in horses in Brazil. Vet Parasitol. 2013;194:35–9.

9. Traversa D, Castagna G, von Samson-Himmelstjerna G, Meloni S, Bartolini R, Geurden T, et al. Efficacy of major anthelmintics against horse cyathostomins in France. Vet Parasitol. 2012;188:294–300.

10. Relf VE, Lester HE, Morgan ER, Hodgkinson JE, Matthews JB. Anthelmintic efficacy on UK Thoroughbred stud farms. Int J Parasitol. 2014;44:507–14.

11. Geurden T, van Doorn D, Claerebout E, Kooyman F, De Keersmaecker S, Vercruysse J, et al. Decreased strongyle egg re-appearance period after treatment with ivermectin and moxidectin in horses in Belgium, Italy and The Netherlands. Vet Parasitol. 2014;204:291–6.

12. Mfitilodze MW, Hutchinson GW. Prevalence and abundance of equine strongyles (Nematoda: Strongyloidea) in tropical Australia. J Parasitol. 1990; 76:487–94.

13. Ogbourne CP. The prevalence, relative abundance and site distribution of nematodes of the subfamily Cyathostominae in horses killed in Britain. J Helminthol. 1976;50:203–14.

14. Kuzmina TA, Kharchenko VA, Starov AI, Dvojnos GM. Analysis of the strongylid nematodes (Nematoda: Strongylidae) community after deworming of brood horses in Ukraine. Vet Parasitol. 2005;131:283–90.

15. Traversa D, Milillo P, Barnes H, von Samson-Himmelstjerna G, Schurmann S, Demeler J, et al. Distribution and species-specific occurrence of cyathostomins (Nematoda, Strongylida) in naturally infected horses from Italy, United Kingdom and Germany. Vet Parasitol. 2010;168:84–92.

16. Collobert-Laugier C, Hoste H, Sevin C, Dorchies P. Prevalence, abundance and site distribution of equine small strongyles in Normandy, France. Vet Parasitol. 2002;110:77–83.

17. Kuzmina TA, Kharchenko VO. Anthelmintic resistance in cyathostomins of brood horses in Ukraine and influence of anthelmintic treatments on strongylid community structure. Vet Parasitol. 2008;154:277–88.

18. Hodgkinson JE, Freeman KL, Lichtenfels JR, Palfreman S, Love S, Matthews JB. Identification of strongyle eggs from anthelmintic-treated horses using a PCR-ELISA based on intergenic DNA sequences. Parasitol Res. 2005;95:287–92.

19. Ionita M, Howe DK, Lyons ET, Tolliver SC, Kaplan RM, Mitrea IL, et al. Use of a reverse line blot assay to survey small strongyle (Strongylida: Cyathostominae) populations in horses before and after treatment with ivermectin. Vet Parasitol. 2010;168:332–7.

20. Kooyman FN, van Doorn DC, Geurden T, Mughini-Gras L, Ploeger HW, Wagenaar JA. Species composition of larvae cultured after anthelmintic treatment indicates reduced moxidectin susceptibility of immature *Cylicocylus* species in horses. Vet Parasitol. 2016;227:77–84.

21. van Doorn DC, Ploeger HW, Eysker M, Geurden T, Wagenaar JA, Kooyman FN. *Cylicocylus* species predominate during shortened egg reappearance period in horses after treatment with ivermectin and moxidectin. Vet Parasitol. 2014;206:246–52.

22. Kuzmina TA, Lyons ET, Tolliver SC, Dzeverin II, Kharchenko VA. Fecundity of various species of strongylids (Nematoda: Strongylidae)-parasites of domestic horses. Parasitol Res. 2012;111:2265–71.

23. Stancampiano L, Mughini Gras L, Poglayen G. Spatial niche competition among helminth parasites in horse's large intestine. Vet Parasitol. 2010;170: 88–95.

24. Kuzmina TA, Dzeverin I, Kharchenko VA. Strongylids in domestic horses: Influence of horse age, breed and deworming programs on the strongyle parasite community. Vet Parasitol. 2016;227:56–63.

25. Poulin R. Sexual inequalities in helminth infections: a cost of being a male? Am Nat. 1996;147:287–95.

26. Zuk M, McKean KA. Sex differences in parasite infections: patterns and processes. Int J Parasitol. 1996;26:1009–23.

27. Grzybek M, Bajer A, Behnke-Borowczyk J, Al-Sarraf M, Behnke JM. Female host sex-biased parasitism with the rodent stomach nematode *Mastophorus muris* in wild bank voles (*Myodes glareolus*). Parasitol Res. 2015;114:523–33.

28. Bucknell DG, Gasser RB, Beveridge I. The prevalence and epidemiology of gastrointestinal parasites of horses in Victoria, Australia. Int J Parasitol. 1995;25:711–24.

29. Debeffe L, McLoughlin PD, Medill SA, Stewart K, Andres D, Shury T, et al. Negative covariance between parasite load and body condition in a population of feral horses. Parasitology. 2016;143:983–97.

30. Kornaš S, Sallé G, Skalska M, David I, Richard A, Cabaret J. Estimation of genetic parameters for resistance to gastro-intestinal nematodes in pure blood Arabian horses. Int J Parasitol. 2015;45:237–42.
31. Kornaś S, Cabaret J, Skalska M, Nowosad B. Horse infection with intestinal helminths in relation to age, sex, access to grass and farm system. Vet Parasitol. 2010;174:285-91.

32. Osterman Lind E, Eysker M, Nilsson O, Uggla A, Hoglund J. Expulsion of small strongyle nematodes (*Cyathostomum* spp.) following deworming of horses on a stud farm in Sweden. Vet Parasitol. 2003;115:289-99.

33. Raynaud JP. Study of the efficiency of a quantitative coproscopic technic for the routine diagnosis and control of parasitic infestations of cattle, sheep, horses and swine. Ann Parasitol Hum Comp. 1970;45:321-42.

34. R Core Team. R: A Language and Environment for Statistical Computing. Vienna: R Foundation for Statistical Computing; 2016.

35. Harrell FE, Dupont C. Hmisc: Harrell Miscellaneous. 2017.

36. Sallé G, Cabaret J. Meta-analysis of cyathostomins community structure and diversity. Liverpool: International Conference of the World Association for the Advancement of Veterinary Parasitology; 2015.

37. Oksanen J, Blanchet FG, Friendly M, Kindt R, Legendre P, McGlinn D, et al. vegan: Community Ecology Package. R package version 2.4-3. 2017.

38. Withthaker RH. Vegetation of Siskiyou mountains, Oregon and California. Ecol Monogr. 1960;30:279-338.

39. Shannon CE. A mathematical theory of communication. Bell Syst Tech J. 1949;27:379-423.

40. Anderson MJ, Ellingsen KE, McArdle BH. Multivariate dispersion as a measure of beta diversity. Ecol Lett. 2006;9:683-93.

41. Anderson MJ. Distance-based tests for homogeneity of multivariate dispersions. Biometrics. 2006;62:245-53.

42. Griffith DM, Veech JA, Marsh CJ. cooccur: Probabilistic Species Co-Occurrence Analysis in R. J Stat Softw Code Snips. 2016;69(2):1-17. https://doi.org/10.18637/jss.v069.c02.

43. Veech JA. A probabilistic model for analysing species co-occurrence. Glob Ecol Biogeogr. 2013;22:252-60.

44. Briatte F. ggnet: Functions to plot networks with ggplot2. 2018. https://briatte.github.io/ggnet/.

45. Hanski I. Dynamics of regional distribution - the core and satellite species hypothesis. Oikos. 1982;38:210-21.

46. Boxell AC, Gibson KT, Hobbs RP, Thompson RCA. Occurrence of gastrointestinal parasites in horses in metropolitan Perth, Western Australia. Aust Vet J. 2004;82:91-5.

47. Lind EO, Hoglund J, Ljungstrom BL, Nilsson O, Uggla A. A field survey on the distribution of strongyle infections of horses in Sweden and factors affecting faecal egg counts. Equine Vet J. 1999;31:68-72.

48. Poulin R. Interactions between species and the structure of helminth communities. Parasitology. 2001;122(Suppl):S3-11.

49. Holmes JC. The structure of helminth communities. Int J Parasitol. 1987; 17:203-8.

50. Behnke JM. Structure in parasite component communities in wild rodents: predictability, stability, associations and interactions ... or pure randomness? Parasitology. 2008;135:751-66.

51. Bottomley C, Isham V, Basanez MG. Population biology of multispecies helminth infection: competition and coexistence. J Theor Biol. 2007;244:81-95.

52. Angert AL, Huxman TE, Chesson P, Venable DL. Functional tradeoffs determine species coexistence via the storage effect. Proc Natl Acad Sci USA. 2009;106:11641-5.

53. Duthie AB, Abbott KC, Nason JD. Trade-offs and coexistence in fluctuating environments: evidence for a key dispersal-fecundity trade-off in five nonpollinating fig wasps. Am Nat. 2015;186:151-8.

54. Stear MJ, Bishop SC. The curvilinear relationship between worm length and fecundity of *Teladorsagia circumcincta*. Int J Parasitol. 1999;29:777-80.

55. Smith HJ. Strongyle infections in ponies. I. Response to intermittent thiabendazole treatments. Can J Comp Med. 1976;40:327-33.

56. Gibson TE. The effect of repeated anthelmintic treatment with phenothiazine on the faecal egg counts of housed horses, with some observations on the life cycle of *Trichonema* spp. in the horse. J Helminthol. 1953;27:29-40.

Ready to submit your research? Choose BMC and benefit from:

• fast, convenient online submission

• thorough peer review by experienced researchers in your field

• rapid publication on acceptance

• support for research data, including large and complex data types

• gold Open Access which fosters wider collaboration and increased citations

• maximum visibility for your research: over 100M website views per year

At BMC, research is always in progress.

Learn more biomedcentral.com/submissions

<figure><img src="image_10.png" /></figure>