Introduction
Latitudinal gradients in species diversity are among the best-studied biogeographical patterns (Willig et al., Reference Willig, Kaufman and Stevens2003). Several studies to date have investigated the relationship between latitude and parasite species diversity (Rohde & Heap, Reference Rohde and Heap1998; Choudhury & Dick, Reference Choudhury and Dick2000; Poulin, Reference Poulin2001; Poulin & Morand, Reference Poulin and Morand2004). However, there have been no serious attempts to document how the composition, not the richness, of parasite communities changes as a function of latitude. All else being equal, the types of parasites making up the helminth community of a vertebrate host in temperate areas may differ from those in a similar vertebrate in the tropics, because many environmental variables correlate with latitude, including the diversity and availability of intermediate hosts.
As an example, we can consider how latitude might correlate with the contribution of trematodes to the parasite communities of vertebrates. Trematodes require a mollusc, usually a gastropod, as first intermediate host, and often another invertebrate as second intermediate host before they reach a vertebrate host. Generally speaking, the diversity of molluscs increases toward the equator (Macpherson, Reference Macpherson2002; Witman & Roy, Reference Witman and Roy2009). Given the strict host specificity of trematodes for their mollusc host, we might expect that the greater the number of gastropods in an area, the higher the total number of trematode species present. Thus, all else being equal, the proportion of trematodes in the parasite communities of tropical vertebrates could be higher than in similar hosts living at higher latitudes. The validity of such predictions, however, relies on several assumptions, and empirical data are necessary to determine whether the proportion of any taxonomic group of parasites changes as a function of latitude.
Other factors can also influence the composition of parasite communities. For instance, species-rich communities may tend to be dominated by small-bodied parasite taxa, since space saturation might constrain the body size of new species that could colonize an existing community (Poulin et al., Reference Poulin, Beveridge and Spratt2008). For instance, in a relatively species-rich helminth community, further diversification may only be possible via the addition of small parasites like trematodes and not large ones like cestodes. In addition, the nature of the host can affect the composition of its parasite community. Hosts and parasites often share a long co-evolutionary history, such that the parasite taxa found in an extant host are, to a large degree, inherited from the host's ancestor (Brooks & McLennan, Reference Brooks and McLennan1993). This legacy, combined with the host's ecological traits that determine what kinds of parasites it is exposed to, shape its parasite fauna. Furthermore, marked diversification of particular parasites is often only achieved in specific types of hosts. For example, some nematodes have radiated extensively into ‘species flocks’ in tortoises and kangaroos (Schad, Reference Schad1963; Beveridge & Spratt, Reference Beveridge and Spratt1996; Beveridge et al., Reference Beveridge, Chilton and Spratt2002), but this sort of diversification has not been reported in fish or bird hosts. Therefore, in addition to latitude, the type of host used and the species richness of the community may also influence the relative contribution of major taxa to species composition of parasite communities.
Here, we investigate the determinants of taxonomic composition of parasite communities in vertebrate hosts. We limit our analyses to adult and larval stages of helminth endoparasites (trematodes, cestodes, nematodes and acanthocephalans) of fish, birds and mammals, since these have been the subjects of numerous investigations. We address the following questions: (1) Does the taxonomic composition of helminth communities show any latitudinal gradients? (2) Is taxonomic composition different in species-rich helminth communities than in those with fewer species? (3) What are the quantitative differences in taxonomic composition between the helminth communities of fish, birds and mammals? Our analyses take advantage of a very large dataset to extract general patterns that may shed light on the factors governing the assembly of parasite communities over evolutionary time scales.
Material and methods
We used a dataset compiled by Poulin & Leung (Reference Poulin and Leung2010) from published surveys of endohelminth communities in vertebrates. This compilation provides a large, representative and unbiased sample derived from 545 studies published between 1936 and 2009. The dataset only includes studies that provided a clear list of all species of trematodes, cestodes, nematodes and acanthocephalans found in one sample of individual hosts from the same species. Therefore, each entry in the dataset corresponds to one parasite survey for one host species from one study, i.e. one local parasite community; some host species have been the subject of a few separate studies, and thus for those host species there are two or more entries in the dataset. For each entry, the following information was recorded: (1) the name of the host species; (2) the host type, i.e. fish, bird or mammal, and the taxonomic order to which the host species belonged (based on information at www.fishbase.org); (3) the total number of endohelminth species, as well as the numbers of trematode, cestode, nematode and acanthocephalan species; and (4) the latitude at which the host sample was collected, or the latitude at the centre of the area where the hosts were obtained if they came from a region rather than a precise locality. Most available surveys came from the northern hemisphere. Thus, in preliminary analyses latitude was treated in three different ways: by using latitude regardless of whether it was north or south, by assigning negative values to southern latitudes to differentiate them from northern ones, and by using only northern latitudes. Since the effect of latitude emerging from the analyses was the same whichever of these three approaches was used, below we only present results from analyses that used latitude regardless of whether it was north or south.
Studies of parasite species richness generally correct for sampling effort, since the more hosts are examined in a particular survey, the more parasite species are found (Walther et al., Reference Walther, Cotgreave, Price, Gregory and Clayton1995). However, here sampling effort was not treated as a confounding variable since the analyses did not focus on how many species were found, but on what taxon the species that were found belonged to; preliminary analyses showed that sampling effort did not affect the proportions of species belonging to different groups (correlations with proportions of species belonging to each of the four groups: all P>0.23).
The proportion of parasite species in a community that belong to higher groups like trematodes or nematodes is, in part, a biological property of the host species, because of host–parasite co-evolutionary history, and, in part, a property of the environment, which determines the presence of essential intermediate hosts and the survivorship of infective stages. Host phylogeny within fish, birds and mammals may thus influence the composition of parasite communities. Therefore, in the absence of a comprehensive species-level phylogeny of vertebrates, we included host taxonomic order as a variable in the analyses to account for possible higher-order phylogenetic influences. Furthermore, some host species contribute two or more entries in the dataset, and to avoid potential pseudoreplication, the analysis was conducted both at the survey level (across all entries in the dataset) and at the host species level. For the latter analysis, for each host species with two or more entries in the dataset, only the survey with the highest number of parasite species was retained as the most representative of the parasite fauna of that host species; other entries for that host were excluded from the host-species-level analyses.
There were therefore two levels of analysis (across all surveys, and across host species) for each of four measures of taxonomic composition (proportion of trematodes, of cestodes, of nematodes and of acanthocephalans per community), for a total of eight versions of the analysis of the determinants of the taxonomic composition of parasite communities. For each of these eight analyses, we used generalized linear models (GLM) implemented with the statistical software JMP 7.0 (SAS Institute Inc., Cary, North Carolina, USA). Since the response variable was always a proportion, i.e. the proportion of species belonging to a particular parasite group, all models had a binomial error structure and logit link function. The explanatory variables included in the starting models were: host type (fish, bird or mammal), host order (nested within host type), parasite species richness (or the log10-transformed total number of helminth taxa in the community), latitude, and the interaction between host type and latitude. After starting with a full model, significance levels were based on the deviance explained by each factor, based on χ2 statistics, following backward stepwise elimination of unimportant (P>0.10) terms. Only the final models are shown in the results.
Results
The dataset included a total of 950 helminth endoparasite communities from 650 different host species (see Poulin & Leung, Reference Poulin and Leung2010, for the full dataset). The surveys consisted of 419 on fish hosts, 161 on bird hosts and 370 on mammal hosts. The parasite species richness in these surveys ranged from 1 to 52 helminth species per community, with a median of eight species (0–28 trematodes, 0–45 cestodes, 0–40 nematodes and 0–10 acanthocephalans). The surveys were conducted between 78°N and 80°S, though the vast majority (828 out of 950) were from the northern hemisphere, and most were from temperate latitudes (fig. 1). We found no overall latitudinal gradient in parasite species richness (latitude versus log no. helminth species per survey: r = 0.046, P = 0.153); when the four helminth groups were treated separately, only very weak correlations were obtained (all r < 0.13). In other words, the number of parasite species per community showed no consistent variation as a function of latitude, and therefore any change in taxonomic composition (see below) is independent of latitudinal gradients in richness.
Frequency distribution of latitudes at which helminth communities were surveyed, across 950 surveys in fish, bird and mammal hosts. Negative values represent southern latitudes, and positive values represent northern ones.
In the GLMs testing for the determinants of species composition in helminth community, the effect of host taxonomic order (nested within host type) only appeared influential in analyses of the proportion of trematode species per community. The influence of host order only emerged in preliminary tests using specific subsets of factors, and it was not retained in our final model which started with all factors prior to backward elimination. Therefore, the results below are unlikely to be masked by any influence of higher-order host taxonomic affiliation on the inclusion of trematodes in helminth communities. Host taxonomic order had no detectable effect on the proportions of other types of helminths (cestodes, nematodes, acanthocephalans) among the surveys included in our analyses.
Almost identical results were obtained in the GLMs whether they were carried out across all surveys, or across host species, i.e. including a single survey per host species (table 1). Parasite species richness was not retained in any of the final models. However, both host type and latitude significantly correlated with the proportion of helminth species belonging to various groups (table 1). In general, helminth communities in fish included a higher proportion of trematodes and acanthocephalans than those in other host types, whereas communities in birds and especially mammals comprised a higher proportion of nematodes than those in fish (fig. 2). The proportion of trematodes in helminth communities tended to be higher at higher latitudes, though this effect was not quite significant (table 1, fig. 3). This pattern was subject to a weak interaction between host type and latitude (table 1), resulting mainly from the fact that the proportion of trematodes per community shows practically no variation in fish hosts, but tends toward its lowest values in the tropics for bird and mammal hosts. However, there was a clear and significant increase in the proportion of nematodes per community toward lower latitudes (fig. 3). Indeed, nematodes make up nearly half of the species in helminth communities in the tropics, independently of host type.
Results of generalized linear models of the factors influencing taxonomic composition in helminth communities of vertebrates. Eight analyses are shown, at two levels (across all surveys, across host species) and for four measures of composition (proportion of trematodes, cestodes, nematodes and acanthocephalans). The χ2 tests assess the significance of the deviance explained by each factor. df, Degrees of freedom.

Mean ( ± SE) composition of parasite communities with respect to the proportions of trematodes (shaded bars), cestodes (hatched bars), nematodes (white bars) and acanthocephalans (black bars), in three types of vertebrate hosts, fish (419 surveys), birds (161 surveys) and mammals (370 surveys).
Mean ( ± SE) proportions of trematodes (shaded bars) and nematodes (white bars) in parasite communities of vertebrates, as a function of the latitude at which the community occurs. Numbers above the bars indicate the number of communities surveyed in each latitudinal band.
Discussion
As reported in other studies (Rohde & Heap, Reference Rohde and Heap1998; Choudhury & Dick, Reference Choudhury and Dick2000; Poulin, Reference Poulin2001), the number of endohelminth parasite species per community showed no consistent variation as a function of latitude; in other words, there was no general gradient in parasite diversity among the 950 surveys included in our study. However, we found some relatively clear patterns in the taxonomic composition of helminth communities: (1) trematodes and acanthocephalans account for a higher proportion of species in helminth communities of fish, whereas nematodes represent a higher proportion of the species in communities of bird and especially mammal hosts; (2) the proportion of trematodes in helminth communities of birds and mammals increased toward higher latitudes; and (3) the proportion of nematodes per community increased toward lower latitudes regardless of the type of host.
Trematodes use snails and often a second intermediate host before reaching their vertebrate definitive host, whereas acanthocephalans use an arthropod, commonly a crustacean, as intermediate host. Most of these intermediate host taxa are aquatic, which may account for the fact that these two parasite groups account for higher proportions of helminth communities in fish hosts than in birds and especially mammals. In contrast, nematodes use a variety of life cycles, many of which require no intermediate hosts. Anderson (Reference Anderson1996) noted that fish have generally few nematode parasites compared to terrestrial vertebrates. He attributed this to the terrestrial origins of parasitic nematodes, and their relatively slow colonization of aquatic habitats. Our large-scale quantitative evidence supports Anderson's (Reference Anderson1996) qualitative observation. Nematodes are particularly dominant in communities of mammalian hosts. It is worth noting that radiations and species flocks of parasitic nematodes have only been reported from large grazing mammals, such as elephants, horses and kangaroos (Inglis, Reference Inglis1971; Bucknell et al., Reference Bucknell, Hoste, Gasser and Beveridge1996; Beveridge & Spratt, Reference Beveridge and Spratt1996; Beveridge et al., Reference Beveridge, Chilton and Spratt2002), and never from birds or fish. Such radiations can result in nematodes accounting for a substantial proportion of the parasites in a community.
Because the diversity of molluscs generally increases toward the equator (Macpherson, Reference Macpherson2002; Witman & Roy, Reference Witman and Roy2009), one obvious prediction was that they may harbour a richer fauna of trematodes, which could translate into a high proportion of trematodes in the helminth communities of definitive hosts. In fact, we found the opposite: trematodes account for a slightly higher proportion of species in helminth communities at higher latitudes, and not in the tropics. Interestingly, a study of trematode parasitism in their snail intermediate hosts has found that snails from high latitude areas tend to harbour slightly richer trematode communities than those of low latitude areas (Poulin & Mouritsen, Reference Poulin and Mouritsen2003). In contrast, small crustaceans serving as second intermediate hosts for trematodes harbour richer trematode communities at low latitudes than in high latitude areas (Thieltges et al., Reference Thieltges, Fredensborg, Studer and Poulin2009). These contrasting latitudinal patterns in different hosts, i.e. for different stages in the life cycles of trematodes, are difficult to reconcile. Although trematode species richness in crustacean intermediate hosts is higher at low than at high latitudes, perhaps a disproportionate number of trematode species use other types of second intermediate hosts at high latitudes. It has been hypothesized that truncated life cycles, in which the second intermediate host is bypassed, are more common at high latitudes (Galaktionov & Skirnisson, Reference Galaktionov and Skirnisson2007); this may explain the congruence between the present results and those of Poulin & Mouritsen (Reference Poulin and Mouritsen2003) and the mismatch between the present results and those of Thieltges et al. (Reference Thieltges, Fredensborg, Studer and Poulin2009). Latitudinal variation in host specificity may provide an alternative explanation. For certain types of parasites, host specificity is relaxed at higher latitudes (see Krasnov et al., Reference Krasnov, Shenbrot, Khokhlova, Mouillot and Poulin2008), and lower host specificity for the definitive hosts at high latitudes could certainly account for our results, independently of latitudinal gradients in species richness found in intermediate hosts.
We can only speculate on possible explanations for the third pattern identified by our analyses, i.e. the finding that the proportion of nematodes per community increased toward lower latitudes regardless of the type of host. One possibility is that there may be a more diverse fauna of intermediate hosts for nematodes in the tropics, and also more complex food webs, leading to a greater potential for the diversification of nematodes and their trophic acquisition by many hosts. Alternatively, warmer conditions, or less variable temperatures, at lower latitudes may favour the survival of the free-living stages of nematodes (many nematodes do not use intermediate hosts), which could promote the maintenance of a rich nematode fauna in tropical and subtropical areas, and their overrepresentation in helminth communities of vertebrates in those areas. These and other possible explanations will require further investigation. Although several questions remain unanswered, our study has highlighted the kind of novel insights that can be gained by searching for latitudinal gradients not only in parasite species richness, which appear inconsistent (Rohde & Heap, Reference Rohde and Heap1998; Choudhury & Dick, Reference Choudhury and Dick2000; Poulin, Reference Poulin2001; Poulin & Morand, Reference Poulin and Morand2004), but also in the taxonomic compositions of parasite assemblages.
Acknowledgements
We are grateful to Isabel Blasco-Costa and Haseeb Randhawa for useful comments on an earlier draft of this manuscript. Financial support came in part via a grant from the Marsden Fund.

