Published online by Cambridge University Press: 13 May 2004
Recent studies of the forces behind the diversification of parasite assemblages have shed light on many aspects of parasite biodiversity. By using only parasite species richness as their measure of diversity, however, previous investigations have ignored the relatedness among parasite species and the taxonomic structure of the assemblages, which contain much information about their evolutionary origins. Here, we performed a comparative analysis across 50 species of fish from the coast of Brazil; we evaluated the effects of several host traits (body size, social behaviour, feeding habits, preference for benthic vs. pelagic habitats, depth range, and ability to enter brackish waters) on the diversity of their assemblages of metazoan parasites. As measures of diversity, we used parasite species richness, as well as the average taxonomic distinctness of the assemblage and its variance; the latter measures are based on the average taxonomic distance between any two parasite species in an assemblage. Unlike parasite species richness, taxonomic distinctness was unaffected by the number of host individuals examined per species. Fish body length proved to be the main predictor of parasite species richness, even when controlling for the confounding influences of host phylogeny and sampling effort, although it did not correlate with measures of parasite taxonomic distinctness. Predatory fish also had higher parasite species richness than planktivores, but this trend could not be confirmed using phylogenetically independent contrasts between host taxa. The main host feature associated with the taxonomic diversity of parasites was schooling behaviour, with schooling fish having more taxonomically diverse parasite assemblages than those of their non-schooling relatives. When focusing on endoparasite species only, both predatory feeding habits and a broad depth range were associated with the taxonomic distinctness of parasites. Our results suggest that certain host traits (i.e. body size) determine how many parasite species a host can accumulate over evolutionary time, whereas different host features influence the processes causing the taxonomic diversification of parasite assemblages.
Parasite assemblages are playing an increasingly important role as models for the study of biodiversity and biogeography (Poulin & Morand, 2000). The events leading to new species joining parasite assemblages are relatively well understood, and it is possible to search for the key factors that have caused certain parasite assemblages to diversify more than others over evolutionary time (Poulin, 1998; Page, 2003). Yet, the results currently available are inconsistent (Poulin, 1997; Morand, 2000). In fish, for instance, some comparative studies have found that host body size is a good predictor of parasite species richness (Price & Clancy, 1983; Bell & Burt, 1991; Poulin, 1995; Gregory, Keymer & Harvey, 1996; Sasal, Morand & Guégan, 1997) whereas others have found no effect of host size (e.g., Morand et al. 2000, though this study focused on a single fish family). The same is true for a range of other host features (see Poulin, 1997; Morand, 2000). It is therefore difficult to assess the relative importance of different host traits for the evolution of parasite diversity in general.
One reason for the inconsistent results may be that most data sets used in previous analyses have been compiled from many different sources. The differences in the methods used to detect and identify parasites among different sources can generate too much error and statistical noise for any existing signal to emerge. In addition, the different host species included in the same data set often come from different geographical areas; availability of (and exposure to) different parasite species varies across regions, adding another element of variability to such data sets. Ideally, one would examine interspecific variation in parasite diversity among a group of host species from the same area, all examined in a standard way by the same researcher(s).
Another reason why previous results have been inconsistent may be that the measure of parasite diversity used in earlier studies is not the most appropriate one. To date, all have used species richness as their sole measure of the diversity of parasite assemblages. Richness is a convenient measure, but it does not capture all facets of diversity (Purvis & Hector, 2000). It ignores the evolutionary relationships among species coexisting in an assemblage (see Shimatani, 2001). Applied to parasite assemblages, measures of diversity that incorporate information on the relationships among parasite species can shed light on how the assemblage has formed. New species join a parasite assemblage in one of two ways (Poulin, 1998; Page, 2003). First, they can originate from within the assemblage by intra-host speciation, when one parasite species undergoes speciation without the host also speciating. The large numbers of congeneric parasite species occurring in certain assemblages indicate that parasite lineages can radiate within single host species over time (Schad, 1963; Kennedy & Bush, 1992; Beveridge, Chilton & Spratt, 2002). Second, new parasite species can arrive when a parasite from a sympatric host species colonizes the assemblage by switching hosts (Paterson & Gray, 1997). Studies of parasite diversity that focus on species richness cannot distinguish between the different origins of parasite species in an assemblage. If the parasites in an assemblage form a narrow taxonomic grouping, intra-host speciation may have been common, or the host has only been colonized by parasites from a limited taxonomic spectrum. If, in contrast, the parasites in an assemblage are completely unrelated, then the assemblage has formed by repeated colonizations from a broad taxonomic range of parasites, indicating that the host is a ‘magnet’ for many types of parasites. We need a measure of diversity that goes beyond mere species richness, a measure that takes into account the relationships between parasite species in an assemblage if we are to figure out which host features attract parasites.
Predictions about which host features may promote the diversification of parasite assemblages have come from two different theoretical frameworks. First, following from island biogeographical theory (MacArthur & Wilson, 1967; Kuris, Blaustein & Alió, 1980), host features that promote high rates of parasite speciation or colonization by new parasite species should be associated with high parasite diversity. For instance, larger-bodied hosts should be able to accommodate more parasite species than small ones; they may also incur higher exposure to internal parasites because of the quantities of food they ingest, and to external parasites because their larger surface area facilitates contact with infective stages. Other host features that increase exposure to different kinds of parasites should also lead to higher colonization rates. These include a broad geographical range that overlaps with the distributions of several other host species, from which new parasite species can be acquired, as well as the type and breadth of habitats used by the host. Empirical support exists for most of these variables, but with the inconsistencies alluded to above (see reviews in Poulin, 1997; Morand, 2000).
The second theoretical source of predictions regarding parasite diversity has been epidemiological modelling (Dobson & Roberts, 1994; Roberts et al. 2002). The models indicate that host population density, which regulates the contact rate between parasite infective stages and hosts, is the key factor determining whether a parasite species can invade and persist in a host population. In comparisons among different host species, those occurring at higher population density should harbour more species of parasites, because they exceed the persistence threshold of more parasite species than hosts with low population density. Comparative studies on mammals have found empirical support for this prediction (Morand & Poulin, 1998; Arneberg, 2002). In hosts such as fish, social behaviour determines host density. Many workers have found that schooling fish species are used by more species of parasites than solitary species, for both external parasites (Caro, Combes & Euzet, 1997; Sasal & Morand, 1998; Raibaut, Combes & Benoit, 1998) and all parasites combined (Ranta, 1992).
All the host features mentioned above are expected to influence parasite species richness; what is not known is how they should affect the taxonomic structure of the parasite assemblage. Certain host traits may promote the acquisition of a broad taxonomic range of parasite species, whereas other host traits may lead to the acquisition of the same number of parasite species, but from a narrower taxonomic spectrum. In the present study, we examine the relationship between different features of host species and the diversity of metazoan parasite assemblages across species of fish hosts from Brazilian coastal waters. We use two measures of diversity, parasite species richness and the taxonomic distinctness of species within an assemblage. The latter measure is based on the average taxonomic distances between parasite species (Clarke & Warwick, 1998, 1999; Warwick & Clarke, 2001). Our study is the first to use indices of taxonomic distances as measures of parasite diversity.
All fish were collected by local fishermen from the coastal waters off the state of Rio de Janeiro, Brazil (latitude 21–23 °S), during the period 1991 to 2002. The fishes were identified according to Figueiredo & Menezes (1978, 1980, 2000) and Menezes & Figueiredo (1980, 1985). Each individual fish was measured (total length) and examined for both external and internal metazoan parasites, using standard parasitological methods. All internal organs were searched for endoparasites; washings from the gills and the lumen of the gut were passed through a sieve (154 μm mesh size) to recover even the smallest parasites (further methodological details provided in the references at the end of this paragraph). The metazoan parasites recovered consisted of platyhelminths (trematodes, cestodes, monogeneans), nematodes, acanthocephalans, hirudineans and crustaceans. All fish dissections and collection of the parasites were made using the same methods, and all parasite identifications were carried out or confirmed by the same person (J.L.L.). Thus the data do not suffer from the problems associated with data sets compiled from different sources and based on different methods. Quantitative information on the prevalence and intensity of infection by different parasites and on parasite community descriptors has been published previously for about half of the fish species in our data set (Luque et al. 1996, 2000, 2002, 2004; Takemoto, Luque & Amato, 1996; Knoff, Luque & Amato, 1997; Cezar & Luque, 1999; Chaves & Luque, 1999; Luque & Chaves, 1999; Silva et al. 2000; Luque & Alves, 2001; Tavares et al. 2001, 2004a,b; Alves & Luque, 2001; Alves et al. 2002a,b; Paraguassú, Luque & Alves, 2002; Cordeiro & Luque, 2004; Tavares & Luque, 2004a,b). The present analysis is the first one that combines all these data; the entire data set is shown in Table 1.
For each host species, the total number of fish examined for parasites and the total number of parasite species (i.e. parasite species richness) found were recorded. The number of hosts examined, or sampling effort, is often a key determinant of the number of parasites found in a survey (Walther et al. 1995), and must therefore be included as a potential confounding variable. For the computations of taxonomic distinctness (see below), we used the proposed taxonomies of Brooks & McLennan (1993) and Gibson, Jones & Bray (2002) for trematodes, Boeger & Kritsky (1993) for monogeneans, Khalil, Jones & Bray (1994) for cestodes, Amin (1985) for acanthocephalans, Anderson (2000) for nematodes, Apakupakul, Siddall & Burreson (1999) for hirudineans, Boxshall & Montú (1997) for copepods, Yamaguti (1963) for branchiurans, and Bunkley-Williams & Williams, Jr. (1998) for isopods.
In addition to data on average host body length obtained from the individual fish examined, data on a range of variables were also obtained from Haimovici, Martins & Vieira (1996), Carvalho Filho (1999), Bizerril & Costa (2001), and Froese & Pauly (2003), giving priority to the former three as they deal specifically with fish in Brazilian waters. The following variables were recorded for each fish species: (1) whether the fish species forms schools or not, with species adopting schooling only in some parts of the year (e.g. during the reproductive period) classified as schooling; (2) whether its geographical distribution extends to the Atlantic coast of South America (mainly Brazil), the Atlantic coast of both North and South America, or the whole world; (3) whether the fish's habitat is benthic, benthopelagic or pelagic; (4) whether or not it occasionally enters brackish or estuarine waters; (5) whether it is a predator (the majority of species) or a planktivore; and (6) its depth range, measured as the difference between the deepest and shallowest depths at which it occurs.
For each fish species, we computed the average taxonomic distinctness (Δ+) and the variance in taxonomic distinctness (Λ+) of the parasite species found. This was done in three ways for each fish species: for all parasite species, for endoparasites (trematodes, cestodes, acanthocephalans and nematodes) only, and for ectoparasites (monogeneans, hirudineans and crustaceans) only. Taxonomic distinctness measures for endoparasites and ectoparasites were only computed for fish species harbouring at least 4 parasites of either type, because values obtained for fewer parasite species are not reliable.
When parasite species are placed within a taxonomic hierarchy, based on the Linnean classification into kingdom, phyla, classes, orders, families, genera and species, the average taxonomic distinctness, Δ+, is simply the mean number of steps up the hierarchy that must be taken to reach a taxon common to two parasite species, computed across all possible pairs of parasite species in an assemblage (Clarke & Warwick, 1998, 1999; Warwick & Clarke, 2001). Thus, if two species are congeners, one step (species-to-genus) is necessary to reach a common node in the taxonomic tree; if the two species belong to different genera but the same family, two steps will be necessary (species-to-genus and genus-to-family); and so on, with these numbers of steps averaged across all species pairs. For any given species pair, the number of steps corresponds to half the path length connecting two species in the taxonomic tree, with equal step lengths being postulated between each level in the taxonomic hierarchy. Step lengths are standardized so that the distinctness of two species connected at the highest taxonomic level is set equal to 100 (Clarke & Warwick, 1999); with 6 levels above the species in the taxonomy we used, each step length was thus equal to 16.67. The greater the taxonomic distinctness between parasite species, the higher the number of steps needed, and the higher the value of Δ+. A high value means that on average the parasites in a host population are not closely related. Formally, Δ+ is computed as follows (see Clarke & Warwick, 1998):
where s is the number of parasite species, the double summation is over the set {i=1,…, s; j=1,…, s, such that i<j}, and ωij is the taxonomic distinctness between parasite species i and j, or the number of taxonomic steps required to reach a node common to both.
The index Δ+ measures the average taxonomic distinctness between species, and does not capture all of the taxonomic structure of a set of parasite species. It is possible to have two host species, each harbouring the same number of parasite species and each characterized by an identical value of Δ+, but with one host species clearly supporting a broader taxonomic range of parasites. Asymmetries in the taxonomic distribution of parasite species across higher taxa can sometimes be missed by Δ+, which is only the average taxonomic distinctness; in these situations complementary information can be obtained by examining the variance in taxonomic distinctness (see Clarke & Warwick, 2001; Warwick & Clarke, 2001):
where ϖ is simply the average taxonomic distinctness, or Δ+. The variance Λ+ conveys separate information of how much taxonomic heterogeneity there is among a group of parasite species. Note, however, that Λ+ can only be computed when at least 3 parasite species are found in a host species (it always equals zero with 2 species). To calculate Δ+ and Λ+, a computer program was developed (by D.M.) using borland C++ Builder 5.0.
Analyses were performed separately on all parasite species, on endoparasites only, and on ectoparasites only. We analysed our data in two ways, one that highlighted any existing pattern in the distribution of parasite diversity among host species, and one that emphasized which factors played a role in the evolutionary diversification of parasite assemblages. In the first series of analyses, fish species were treated as independent observations, using standard parametric tests on log-transformed continuous variables. In analyses using categorical variables (e.g. schooling behaviour or feeding habits), it was not possible to use multifactorial ANOVAs because there were too few species in some categories and thus too many empty cells in the factorial matrix. Nevertheless, our analyses allowed us to determine which features of host species are associated with either many parasite species or with taxonomically diverse parasite assemblages. In the second round of analyses, we took host phylogeny into account. Closely related host species are likely to harbour a similar number of parasite species, and possibly taxonomically related parasite species, because these were inherited from a recent common ancestor; this means that they do not represent truly independent statistical observations. We must therefore control for phylogenetic influences when evaluating the effects of host features (body size, schooling or feeding habits, etc.) on the evolution of parasite assemblages. To achieve this, we used the phylogenetically independent contrasts method (Felsenstein, 1985; Harvey & Pagel, 1991), implemented with the CAIC version 2.0 program (Purvis & Rambaut, 1994). Contrasts were derived from a host tree constructed from published studies on the phylogenetic relationships of fish (Nelson, 1994; Carpenter, Collette & Russo, 1995; Reed, Carpenter & deGravelle, 2002; Chen, Bonillo & Lecointre, 2003; Miya et al. 2003). Contrasts were computed on log-transformed data and all regression analyses were forced through the origin (Garland, Harvey & Ives, 1992). We obtained contrasts corrected for the influence of one or more confounding variables (e.g., sampling effort) by taking the residuals of regressions of a selected variable against the potential confounding variables. For dichotomous variables (schooling behaviour, feeding habit, tolerance of brackish waters), contrasts were computed following Burt (1989). The mean value of these contrasts was compared with zero, as expected from the null hypothesis, using one-group two-tailed t-tests.
Because our aim was to expose associations between host features and measures of parasite diversity, we did not apply a Bonferroni correction to our results. Applying the correction could possibly mask interesting trends worthy of further investigation (Moran, 2003).
In total, data from 50 fish host species were obtained (see Table 1). These data were derived from the examination of 3067 individual fish, for an average of 61 hosts per species (range 30 to 162). The number of host species (or sets of independent contrasts) varies among the analyses reported below, because of some missing data on certain variables for a few fish species.
Using data on all metazoan parasites across all fish species, the number of hosts examined per species correlated positively with parasite species richness (r=0·453, N=50, P=0·0009), but not with either Δ+ (r=0·196, N=50, P=0·1723) or Λ+ (r=−0·194, N=48, P=0·1874). Measures of taxonomic distinctness are thus independent of sampling effort, as shown in earlier studies (Rogers, Clarke & Reynolds, 1999). However, parasite species richness correlated positively with Δ+ (r=0·679, N=50, P=0·0001) and negatively with Λ+ (r=−0·374, N=48, P=0·0087), indicating that measures of taxonomic distinctness are influenced by the number of parasite species in a sample. Also, Δ+ and Λ+ covaried negatively (r=−0·774, N=48, P=0·0001), such that an increase in average taxonomic diversity is associated with a decrease in its variance. Identical patterns were also found using phylogenetic contrasts, and the relationships generally hold using either only endoparasites or only ectoparasites, with minor exceptions. Thus, in all comparative analyses below, parasite species richness is corrected for sampling effort, and Δ+ and Λ+ are always corrected for parasite species richness, when appropriate.
Using species values as independent statistical observations, we found that fish length correlated positively with parasite species richness (r=0·445, N=50, P=0·0012): larger fish species tended to harbour more parasite species, independently of sampling effort (Fig. 1). In fact, as seen in Fig. 1, the relationship is better explained by a curvilinear regression with a quadratic term (r2=0·236, P=0·0018) than by a linear relationship (r2=0·198, P=0·0012), suggesting there may be a threshold size beyond which parasite species richness stops increasing. No such correlation with fish length was found when using either Δ+ (r=−0·122, N=50, P=0·3984) or Λ+ (r=0·178, N=48, P=0·2266) as measure of parasite diversity. Host depth range did not correlate with parasite species richness (r=0·076, N=49, P=0·6057), but it covaried negatively, though weakly, with Δ+ (r=−0·272, N=49, P=0·0582) and Λ+ (r=−0·318, N=47, P=0·0293). In other words, fish species occurring over a narrow range of depths tend to possess a taxonomically more diverse parasite fauna than fish with a broad depth range. The effect of other host features on parasite diversity was evaluated with ANOVAs, given that fish species fell into different categories. Only two significant results were obtained: parasite species richness varied significantly among fish species inhabiting different parts of the water column (benthic, benthopelagic and pelagic) (F2,47=3·665, P=0·0332) and between predatory and planktivorous fish species (F1,48=8·893, P=0·0045). Benthic and predatory fish harboured more parasite species than pelagic or planktivorous fish (Fig. 2), even after controlling for the effects of sampling effort and fish length. The extent of the geographical distribution of the fish, and whether or not they can enter brackish habitats, had no effect on any of the measures of parasite diversity. The same was true of schooling behaviour, although non-schooling fish species tended to have lower values of Λ+ than schooling species (F1,46=3·582, P=0·0647).
Results obtained when using phylogenetically independent contrasts show some similarities to those obtained from analyses of raw species values. Again, fish length correlated positively with parasite species richness (Fig. 3; r=0·381, N=40 sets of contrasts, P=0·0153), but not with either Δ+ or Λ+ (both P>0·40). Host depth range correlated negatively with only one of the three measures of parasite diversity, Λ+ (r=−0·452, N=37 sets of contrasts, P=0·0051), but once an outlier was excluded from the analysis, the relationship disappeared (P=0·2664). It was not possible to confirm or refute the effect of feeding habits on parasite diversity, as shown in Fig. 2, because there were only two independent contrasts between fish taxa with different feeding modes. There were no effects, however, of either geographical distribution, use of benthic vs. pelagic habitats, or ability to enter brackish waters, on parasite diversity (all P>0·19). Interestingly, in the contrasts between schooling and non-schooling fish taxa, there was evidence that schooling fish had higher values of Λ+ than non-schooling taxa (t=2·377, N=8 sets of contrasts, P=0·0491), as seen in the analysis using fish species values instead of contrasts. In other words, the variance in taxonomic distinctness of parasite assemblages of schooling fish was higher than that of their non-schooling relatives.
Using host species values as independent observations, we found no correlation between fish length and either endoparasite species richness, Δ+, or Λ+ (all P[ges ]0·19). Host depth range correlated positively with endoparasite species richness (r=0·482, N=45, P=0·0008), but not with either Δ+ or Λ+ (both P[ges ]0·13); this pattern changes when host phylogeny is taken into account (see below). Of the categorical variables considered, only one was associated with endoparasite diversity: endoparasite Δ+ varied significantly between predatory and planktivorous fish species (F1,44=6·957, P=0·0115). Despite very unequal sample sizes, predatory fish harboured endoparasite species that are more taxonomically diverse than those of planktivorous fish (Fig. 4), even after controlling for the effects of species richness.
Results obtained when using phylogenetically independent contrasts confirmed the absence of a relationship between fish length and any of the three measures of endoparasite diversity (all P[ges ]0·12). Host depth range did not relate with endoparasite species richness (P=0·8145), but it correlated positively with Δ+ (r=0·567, N=35 sets of contrasts, P=0·0002); 26 of the 35 sets of contrasts had the same sign, suggesting that the two variables tend to covary (Fig. 5). Thus, as a fish species extends its depth range, it also acquires a broader taxonomic range of endoparasites. There was no significant association between any of the categorical variables and either of the measures of endoparasite diversity, although, in the contrasts between schooling and non-schooling fish taxa, there was a hint that schooling fish had higher values of Λ+ than non-schooling taxa (t=2·062, N=8 sets of contrasts, P=0·0781), as seen in the analysis using all parasite species.
No significant associations were found between any of the host features investigated and any of the three measures of ectoparasite diversity, using either raw species values or phylogenetically independent contrasts. There were, however, too few sets of contrasts to test the relationship between ectoparasite diversity and either feeding habits or schooling behaviour.
The search for the key determinants of parasite biodiversity has blossomed in recent years (see reviews in Poulin, 1997; Morand, 2000; Poulin & Morand, 2000). At the same time, phylogenetic studies have revealed the historical and contemporary determinants of the biogeographical distribution of parasites (Hoberg & Klassen, 2002). Yet, there is no consensus view regarding the role, if any, of various host traits in the evolutionary diversification of parasite faunas. This may in part be due to the nature of previous studies, which have often suffered from important limitations. Here, we found that parasite biodiversity was not distributed randomly among marine fish species with respect to host characteristics. We have addressed some features long thought to be associated with parasite acquisition in ecological time by individual fish, and over evolutionary time by fish species (see Dogiel, Petrushevski & Polyanski, 1961). Certain host features, i.e. host body size and possibly host diet, appear to influence the total number of parasite species exploiting a host species, whereas other host features, i.e. depth range, appear to influence the variety of parasite taxa that a host will accumulate over evolutionary time. Our study includes two major improvements on earlier studies of this kind, making its results particularly interesting.
First, many earlier studies of the diversity of fish parasites have used data from fish species that do not occur in the same geographical areas (e.g., Bell & Burt, 1991; Aho & Bush, 1993; Poulin, 1995; Rohde, Hayward & Heap, 1995; Gregory et al. 1996). While relationships between parasite diversity and host features may still be detected in such studies, they are confounded by the fact that the different fish species come from different areas and are therefore not exposed to the same pool of parasite species. The size of the pool of available parasite species must differ from one geographical area to the next, and thus it can limit how many parasite species a host can acquire over time, regardless of the characteristics of this host species. Our study focused on a set of fish species from the same general area (the coastal waters off the state of Rio de Janeiro, Brazil), thus minimizing any differences in parasite availability. In addition, our study was the first investigation of parasite biodiversity based on a representative data set of marine fishes from the Neotropical Region.
Second, previous studies have almost exclusively focused on parasite species richness as a measure of parasite diversity. While convenient, species richness does not capture all facets of biodiversity (Purvis & Hector, 2000). Our analysis also included measures of the taxonomic diversity of parasite species in the different assemblages. Looking at Table 1, it is easy to find pairs of fish species with similar parasite species richness, but with different measures of parasite taxonomic distinctness, or vice versa. We can therefore look at whether certain host features have allowed some host species to ‘capture’ a narrow or a broad taxonomic range of parasite species, regardless of exactly how many parasite species they have acquired.
Fish size proved to be the main predictor of total parasite species richness in the marine fish we investigated. Based on arguments borrowed from island biogeography theory (Kuris et al. 1980), one might predict that the equilibrium number of parasite species that a host can harbour, when rates of colonization by new parasite species are balanced by extinction rates, should be determined mainly by its body size. Host body size is the best measure we can use for the total amount of nutrient or energy available in a host for parasite exploitation. Our results also suggest that host feeding habits may influence parasite species richness. Predatory fish should be exposed to more infective helminth larvae in their diet than planktivores; over evolutionary time, this should translate in higher parasite colonization rates in predatory fish than in planktivorous fish, pushing up the equilibrium parasite species richness.
Host feeding habits were also associated with the taxonomic diversity of endoparasite species. Clearly, predatory fish may not only acquire more parasites via their diet, they may also be exposed to a broader range of different parasite taxa. In contrast, fish schooling behaviour had no effect on parasite species richness, but it was associated, though weakly, with the taxonomic diversity of the parasite species exploiting a host species. Similarly, host depth range did not correlate with parasite species richness, but it covaried positively with the taxonomic diversity of endoparasites. This relationship was observed even when correcting for endoparasite species richness; thus, for a given number of parasite species, those of fish with broad depth range have wider taxonomic origins than those of fish with a limited depth range. A possibility may be that fish living at many depths feed on a wider array of prey species, and thus face colonization by a more diverse group of parasites, than fish restricted to a narrow depth zone.
On a larger scale, other determinants can influence the diversity of fish parasites in the region. The coastal area of Rio de Janeiro is strongly influenced by upwelling systems and by the subtropical convergence (23 °–30 °S), where the warm Brazilian Current (originating at the southernmost limit of the South Equatorial Current) meets the cold Falkland Current which flows northward (Braga, 2001). The water circulation pattern in this region impedes the exchanges of fish between the Caribbean Sea and the Brazilian coast. These conditions affect the fish species biodiversity in the area (Braga, 2001), and may also influence parasite biodiversity. Moreover, salinity levels, which drop near the Amazon River, may also limit the dispersal of fish and their parasites. Santos & Carbonel (2000) proposed that between 10 °S to 10 °N along the South American coast, low salinity and water current patterns are limiting parasite dispersal. It is thus possible that the biodiversity of fish parasites around Rio de Janeiro may be constrained because physical barriers limit the arrival of further taxa.
All previous studies of parasite diversity in fish and other host taxa have focused exclusively on parasite species richness, the most convenient and widely available measure of biodiversity. Species richness conveys no information on the taxonomic structure of an assemblage. The measures of taxonomic distinctness used here provide a summary of this structure, and they suggest possible evolutionary scenarios. For instance, certain feeding habits or living over a broad depth range may facilitate host-switches leading to completely novel host–parasite associations. For a given parasite species richness, this would lead to a taxonomically diverse assemblage. In contrast, parasites in host species with different features (narrow diet and depth range), may have accumulated in part via intra-host speciation, leading to a taxonomically closer set of species. These are only possibilities, of course. As shown by our results, however, shifting the focus from species richness to taxonomic diversity can cast a different light on the evolution of parasite biodiversity.
J.L.L. was partially supported by a Fellowship from CNPq (Conselho Nacional de Pesquisa e Desenvolvimento Tecnológico) and by two grants from FAPERJ (Fundação de Amparo à Pesquisa do Estado do Rio de Janeiro). R.P. was supported by a James Cook Research Fellowship from the Royal Society of New Zealand during this study.