Relationships between functional alpha and beta diversities of flea parasites and their small mammalian hosts.

We studied the relationships between functional alpha and beta diversities of fleas and their small mammalian hosts in 4 biogeographic realms (the Afrotropics, the Nearctic, the Neotropics and the Palearctic), considering 3 components of alpha diversity (functional richness, divergence and regularity). We asked whether (a) flea alpha and beta diversities are driven by host alpha and beta diversities; (b) the variation in the off-host environment affects variation in flea alpha and beta diversities; and (c) the pattern of the relationship between flea and host alpha or beta diversities differs between geographic realms. We analysed alpha diversity using modified phylogenetic generalized least squares and beta diversity using modified phylogenetic generalized dissimilarity modelling. In all realms, flea functional richness and regularity increased with an increase in host functional richness and regularity, respectively, whereas flea functional divergence correlated positively with host functional divergence in the Nearctic only. Environmental effects on the components of flea alpha diversity were found only in the Holarctic realms. Host functional beta diversity was invariantly the best predictor of flea functional beta diversity in all realms, whereas the effects of environmental variables on flea functional beta diversity were much weaker and differed between realms. We conclude that flea functional diversity is mostly driven by host functional diversity, whereas the environmental effects on flea functional diversity vary (a) geographically and (b) between components of functional alpha diversity.

Nestedness of flea assemblages harboured by small mammalian hosts revisited: phylogenetic and functional nestedness do not follow compositional nestedness.

We studied compositional, phylogenetic, and functional nestedness in the flea assemblages of 14 host species across regions. Our main questions were (a) are a host's flea assemblages compositionally, phylogenetically, or functionally nested? (b) Do similar processes drive these nestedness facets? (d) Are a host's biological traits associated with nestedness of its flea assemblages? Rows of host matrices were ordered by decreasing species richness/the sum of the branch lengths of a phylogenetic tree/functional dendrogram or by decreasing region area or by increasing distance from the centre of a host's geographic range. None of the matrices sorted by species richness/sum of branch lengths were nested from a compositional perspective, but they were significantly nested from phylogenetic and functional perspectives. Compositional, phylogenetic, and functional nestedness of matrices sorted by region area or by distance from the host's geographic range centre varied between hosts. In some hosts, flea assemblages were nested from all three perspectives independently of how matrix rows were sorted, whereas in other hosts, the occurrence of significant nestedness depended on the order of the matrix rows. The degree of phylogenetic and functional nestedness for matrices sorted by the sum of branch lengths was associated with a host species' morphoecological traits and the latitude of its geographic range. We conclude that consideration of nestedness based solely on species composition does not allow a comprehensive understanding of the patterns of parasite community structure. Nestedness should also be considered from phylogenetic and functional perspectives.

Environment and traits affect parasite and host species positions but not roles in flea-mammal networks.

We studied spatial variation in the effects of environment and network size on species positions and roles in multiple flea-mammal networks from four biogeographic realms. We asked whether species positions (measured as species strength [SS], the degree of interaction specialization [d'], and the eigenvector centrality [C]) or the roles of fleas and their hosts in the interaction networks: (a) are repeatable/conserved within a flea or a host species; (b) vary in dependence on environmental variables and/or network size; and (c) the effects of environment and network size on species positions or roles in the networks depend on species traits. The repeatability analysis of species position indices for 441 flea and 429 host species, occurring in at least two networks, demonstrated that the repeatability of SS, d', and C within a species was significant, although not especially high, suggesting that the indices' values were affected by local factors. The majority of flea and host species in the majority of networks demonstrated a peripheral role. A value of at least one index of species position was significantly affected by environmental variables or network size in 41 and 36, respectively, of the 52 flea and 52 host species that occurred in multiple networks. In both fleas and hosts, the occurrence of the significant effect of environment or network size on at least one index of species position, but not on a species' role in a network, was associated with some species traits.

Evoregions of fleas and their small mammalian hosts: Do they coincide?

Combining the biogeography and phylogenetic patterns of parasite-host associations allows a better understanding of the history of parasite­host interactions, which can be achieved via biogeographic regionalization incorporating phylogenetic information. Recently, the concepts of evoregions (regions where a majority of species evolved from one or several ancestors inhabiting these regions) and evolutionary transition zones (regions of high phylogenetic turnover) have been proposed, coupled with a classification approach for these concepts. We applied this approach to 206 flea species and 265 host species of the Palearctic and aimed to identify evoregions and evolutionary transition zones for both fleas and hosts and to understand whether these evoregions and transition zones match each other. We identified 5 evoregions with 3 transition zones for either fleas or hosts, but neither the positions and boundaries of the flea and host evoregions nor the transition zones coincided. Indications of multiple geographic centres of diversification of the same flea lineages suggested that (a) the common evolutionary history of fleas and hosts was characterized by multiple events other than codiversification and that (b) dispersal played an important role in flea and host assemblies. Barriers to dispersal could be represented by landscape features (deserts and mountain ranges) and/or climate differences.

Metabolic rate and ecological traits of ectoparasites: a case study with seven flea species from the Negev Desert.

We studied the relationship between fleas' metabolic rate and their ecological traits, using data on standard metabolic rate (SMR), mean abundance, host specificity, and geographic range size in males and females of seven desert flea species. SMR was measured via mass-specific CO2 emission, whereas host specificity was measured as (a) the mean number of host species used by a flea per region in regions where this flea was recorded; (b) the total number of host species a flea exploited across its geographic range; and (c) the phylogenetic diversity of the flea's hosts. To control for confounding effects of phylogeny when analysing data on multiple species, we applied the Phylogenetic Generalised Least Squares (PGLS) model. We found that the only ecological trait significantly correlating with flea SMR was the phylogenetic diversity of hosts utilized by a flea across its geographic range. The strength of the association between SMR and host phylogenetic diversity was higher in male than in female fleas. We explain the relationship between flea SMR and their host specificity by the necessity of host-opportunistic species to compensate for the high energetic cost of neutralizing multiple defences from multiple hosts by increased SMR.

Compositional and phylogenetic nestedness of host assemblages exploited by generalist ectoparasites across their geographic ranges: drivers and associations with ectoparasite traits.

We investigated compositional and phylogenetic nestedness in the host assemblages of 26 host-generalist fleas across regions within the Palearctic. We asked the following questions: (i) are host assemblages exploited by a flea species compositionally or phylogenetically nested (=C-nested and P-nested, respectively) across regions?; (ii) if yes, what are the processes that generate nestedness, and does phylogenetic nestedness follow the same processes as compositional nestedness?; and (iii) are the biological traits of a flea species associated with its host assemblages' degree of nestedness? Nestedness was calculated for matrices with rows ordered either by decreasing region area (=a-matrices) or increasing distance from the centre of a flea's geographic range (d-matrices). Significant C-nestedness was found in either a- (three fleas) or d-matrices (three fleas) or both (10 fleas). Significant P-nestedness was detected in either a- (three fleas) or d-matrices (four fleas) or both (two fleas). In some but not other species, P-nestedness followed C-nestedness. The probability of C-nestedness to be significant, as well as its degree for d-matrices, was associated with a flea's morphoecological traits, whereas this was not the case for either a-matrices or the P-nestedness for either type of ordered matrices. We conclude that compositional, but not phylogenetic, nestedness is (i) generated by similar mechanisms in many flea species and (ii) may be simultaneously driven by different mechanisms in the same flea. In contrast, mechanisms promoting phylogenetic nestedness differ between flea species and seem to act separately.

Phylogenetic patterns in regional flea assemblages from 6 biogeographic realms: strong links between flea and host phylogenetic turnovers and weak effects of phylogenetic originality on host specificity.

We investigated phylogenetic patterns in flea assemblages from 80 regions in 6 biogeographic realms and asked whether (a) flea phylogenetic turnover is driven by host phylogenetic turnover, environmental dissimilarity or geographic distance; (b) the relative importance of these drivers differs between realms; and (c) the environmental drivers of flea phylogenetic turnover are similar to those of host phylogenetic turnover. We also asked whether the phylogenetic originality of a flea species correlates with the degree of its host specificity and whether the phylogenetic originality of a host species correlates with the diversity of its flea assemblages. We found that host phylogenetic turnover was the best predictor of flea phylogenetic turnover in all realms, whereas the effect of the environment was weaker. Environmental predictors of flea phylogenetic turnover differed between realms. The importance of spatial distances as a predictor of the phylogenetic dissimilarity between regional assemblages varied between realms. The responses of host turnover differed from those of fleas. In 4 of the 6 realms, geographic distances were substantially better predictors of host phylogenetic turnover than environmental gradients. We also found no general relationship between flea phylogenetic originality and its host specificity in terms of either host species richness or host phylogenetic diversity. We conclude that flea phylogenetic turnover is determined mainly by the phylogenetic turnover of their hosts rather than by environmental gradients. Phylogenetic patterns in fleas are manifested at the level of regional assemblages rather than at the level of individual species.

Similar adaptative mechanism but divergent demographic history of four sympatric desert rodents in Eurasian inland.

Phenotypes associated with metabolism and water retention are thought to be key to the adaptation of desert species. However, knowledge on the genetic changes and selective regimes on the similar and divergent ways to desert adaptation in sympatric and phylogenetically close desert organisms remains limited. Here, we generate a chromosome level genome assembly for Northern three-toed jerboa (Dipus sagitta) and three other high-quality genome assemblies for Siberian jerboa (Orientallactaga sibirica), Midday jird (Meriones meridianus), and Desert hamster (Phodopus roborovskii). Genomic analyses unveil that desert adaptation of the four species mainly result from similar metabolic pathways, such as arachidonic acid metabolism, thermogenesis, oxidative phosphorylation, insulin related pathway, DNA repair and protein synthesis and degradation. However, the specific evolved genes in the same adaptative molecular pathway often differ in the four species. We also reveal similar niche selection but different demographic histories and sensitivity to climate changes, which may be related to the diversified genomic adaptative features. In addition, our study suggests that nocturnal rodents have evolved some specific adaptative mechanism to desert environments compared to large desert animals. Our genomic resources will provide an important foundation for further research on desert genetic adaptations.

Latitudinal gradients in body size and sexual size dimorphism in fleas: males drive Bergmann's pattern.

We tested for the effects of latitude and geographic range size (GRS) on body size, leg length, and sexual size dimorphism (SSD) across 103 species of fleas, taking into account phylogenetic between-species relationships. When the data on body size were combined for males and females, the positive correlation between body size and latitude, but not GRS, was revealed. When the analysis was restricted to one sex only, the effect of latitude appeared to be non-significant for females, whereas male body size increased with an increase in latitude. Intraspecific body size variation was not associated with either the latitude or the latitudinal span of the geographic range, independently of which data subset was analyzed. No evidence of association between size-independent tibia length and latitude was found for either females, males, or both sexes combined. The degree of SSD decreased with a decrease in latitude but was not affected by GRS. We conclude that macroecological patterns might be manifested differently in males and females. This should be kept in mind while searching for these patterns.

Phylogeography of the common hamster (Cricetus cricetus): paleoclimatic reconstructions of Late Pleistocene colonization.

The phylogeographic structure of Cricetus cricetus is described across the whole range of the species for the first time. The modern phylogenetic structure was formed 70-45 kyr BP. The most favorable conditions for the species are those of interglacial periods rather than of the periods of glaciations. Earlier an opposite hypothesis was suggested. Thus, the present period is generally favorable for the species' existence. The observed decrease in its abundance and fragmentation of its range in natural habitats are clear evidence of the leading role of anthropogenic factors in this process.

Regional flea and host assemblages form biogeographic, but not ecological, clusters: evidence for a dispersal-based mechanism as a driver of species composition.

We used data on the species composition of regional assemblages of fleas and their small mammalian hosts from 6 biogeographic realms and applied a novel method of step-down factor analyses (SDFA) and cluster analyses to identify biogeographic (across the entire globe) and ecological (within a realm across the main terrestrial biomes) clusters of these assemblages. We found that, at the global scale, the clusters of regional assemblage loadings on SDFA axes reflected well the assemblage distribution, according to the biogeographic realms to which they belong. At the global scale, the cluster topology, corresponding to the biogeographic realms, was similar between flea and host assemblages, but the topology of subtrees within realm-specific clusters substantially differed between fleas and hosts. At the scale of biogeographic realms, the distribution of regional flea and host assemblages did not correspond to the predominant biome types. Assemblages with similar loadings on SDFA axes were often situated in different biomes and vice versa. The across-biome, within-realm distributions of flea vs host assemblages suggested weak congruence between these distributions. Our results indicate that dispersal is a predominant mechanism of flea and host community assembly across large regions.

Phylogenetic signals in flea-host interaction networks from four biogeographic realms: differences between interactors and the effects of environmental factors.

The structure of ecological interaction networks is associated with evolutionary histories of the interacting species. This is reflected by the phylogenetic signals (PS) in these networks when closely related species interact with similar partners because some traits inherited from the ancestors may determine ecological interactions. We investigated PS for small mammalian hosts and fleas in 80 regional interaction networks from four biogeographic realms (the Palearctic, the Nearctic, the Afrotropics, and the Neotropics). We asked (i) whether the relative strength of PS in host-flea networks is similar between hosts and fleas and/or between realms; (ii) how environmental variation affects the PS of hosts and fleas in their interaction networks; and (iii) whether the PS for hosts or fleas is affected by the phylogenetic diversity of either hosts or fleas, respectively. We found that the PS for hosts was stronger than that for fleas in all realms. An environmental effect on the PS for hosts, but not for fleas, was found in three of the four realms (except the Neotropics). In the Palearctic and the Nearctic, a stronger PS was characteristic for cooler and/or drier regions, whereas the opposite was the case for the Afrotropics in regard to precipitation. The phylogenetic diversity of regional host and flea assemblages was not associated with the values of the respective PS in any realm. We conclude that the pattern of the relative strength of the PS for hosts and fleas in their interaction networks is similar in different biogeographic realms with vastly different host and flea faunas. However, the environmental effects on the PS are geographically variable and might be associated with the history of host-flea associations, as well as the spatial pattern of environmental variation, within a realm.

Dark host specificity in two ectoparasite taxa: repeatability, parasite traits, and environmental effects.

We applied the concept of dark diversity (species that may potentially inhabit a locality but are absent) to the host spectrum of a parasite and defined it as dark host specificity (DHS). We studied the trait-associated and geographic patterns of dark host specificity in fleas and gamasid mites parasitic on small mammals, asking the following questions: (a) Is dark host specificity repeatable across populations of the same species? (b) Is it associated with morphological and/or ecological species traits? (c) What are the factors associated with geographical variation in the DHS among populations of the same species? The DHS was repeatable within species with a large proportion of variance among samples, accounted for by differences between species. The average DHS of fleas, but not mites, was affected by parasite traits, with the DHS being higher in fleas with larger geographic ranges, higher characteristic abundance levels, and summer reproduction peaks. In the majority of ectoparasites, the regional DHS decreased with an increase in either structural or phylogenetic host specificity. The associations between the DHS and the environmental or host-associated characteristics of a region were revealed in a few species (eight of 22 fleas and three of 12 mites). The DHS decreased with (a) an increase in air temperature in two fleas, (b) a decrease in precipitation in two fleas, and (c) an increase in regional host species richness (in three fleas and three mites). Overall, our results suggest that dark host specificity in arthropod ectoparasites is a species-specific character associated, to a large extent, with the breadth of their host-related niches, while the influences of parasite traits and local environmental conditions are minor.

Phylogenetic relations and range history of jerboas of the Allactaginae subfamily (Dipodidae, Rodentia).

Five-toed jerboas of the subfamily Allactaginae comprise several complex taxa occurring over a wide distribution range covering a large part of the Eurasian arid belt. In this study, we employed current methods of molecular phylogenetics based on 15 nuclear genes and the mitochondrial gene cytb to revise relations and systematics within Allactaginae. We also applied species distribution modelling projected on paleo-environmental data to reconstruct the geographic patterns of speciation in Allactaginae. We elucidated the intergeneric relationships within this subfamily and clarified interspecies relations within the genus Scarturus. Moreover, our results demonstrate the species status of S. caprimulga; outline the currently understudied diversity within Orientallactaga, Allactaga, and Pygeretmus; and improve the divergence estimates of these taxa. Based on our results from modelling of geographic range fragmentation in allactagines, we suggest the dating and location of speciation events and present hypotheses regarding general habitat niche conservatism in small mammals.

Dark diversity of flea assemblages of small mammalian hosts: effects of environment, host traits and host phylogeny.

An assemblage of species in a locality comprises two components, namely (i) species that are present (realised diversity) and (ii) species from the regional pool that may potentially inhabit this locality due to suitable ecological conditions, but that are absent (dark diversity). We investigated factors affecting the dark diversity of component communities of fleas parasitic on small mammals in the northern Palearctic at two scales. First, we considered the dark diversity of flea assemblages of the same host (for 13 host species) across regions and tested for the effects of environmental factors and the number of available host species on the dark diversity of within-region flea assemblages. Second, we considered the dark diversity of fleas across host species within a region (for 20 regions) and asked whether within-host dark diversity is associated with host phylogeny and/or traits. We found that the dark diversity of flea assemblages harboured by small mammals varied substantially (i) within the same host species across space (in 12 of 13 host species) and (ii) between host species within a region (in eight of 20 regions). The size of the dark diversity of flea assemblages of the same host across regions was generally affected by environmental factors (mainly by the amount of green vegetation), whereas the size of the dark diversity of flea assemblages of a host species within a region was affected by host traits (mainly by the degree of host sociality and the structure of its shelter and, to a lesser degree, by its geographic range size) but was not associated with host phylogenetic affinities. We conclude that application of the dark diversity concept to parasite communities across space or hosts allows a better understanding of the factors affecting the species richness and composition of these communities.

Dispersal-based versus niche-based processes as drivers of flea species composition on small mammalian hosts: inferences from species occurrences at large and small scales.

Biological communities may be assembled by both niche-based and dispersal-based (= historic) processes with the relative importance of these processes in community assembly being scale- and context-dependent. To infer whether (a) niche-based or dispersal-based processes play the main role in the assembly of flea communities parasitic on small mammals and whether (b) the main processes of flea community assembly are scale-dependent, we applied a novel permutation-based algorithm (PER-SIMPER) and the dispersal-niche continuum index (DNCI), to data on the species incidence of fleas and their hosts at two spatial scales. At the larger (continental) scale, we analysed flea communities in four biogeographic realms across adjacent continental sections. At the smaller (local) scale, we considered flea communities across two main regions (lowlands and mountains) and seven habitat types within Slovakia. Our analyses demonstrated that species composition of fleas and their small mammalian hosts depended predominantly on historical processes (dispersal) at both scale. This was true for the majority of biogeographic realms at continental scale (except the Nearctic) and both regions at local scale. Nevertheless, strong niche-based assembly mechanism was found in the Nearctic assemblages. At local scale, the intensity of dispersal processes was weaker and niche-driven processes were stronger between habitats within a region than between mountain and lowland regions. We provide historical and ecological explanations for these patterns. We conclude that the assembly of compound flea communities is governed, to a great extent, by the dispersal processes acting on their hosts and, to a lesser extent, by the niche-based processes.

Spatial and temporal variation of compositional, functional, and phylogenetic diversity in ectoparasite infracommunities harboured by small mammals.

We studied patterns of compositional, functional, and phylogenetic α- and ß-diversity in flea and gamasid mite infracommunities of small Siberian mammals, taking into account host-associated (species) and environmental (biome or sampling period) factors. We asked: (a) How do these factors and their interactions affect infracommunity diversity? (b) Does infracommunity composition, in terms of species, traits, and phylogenetic lineages, deviate from random? (c) Are species, traits, and phylogenetic lineages in infracommunities clustered or overdispersed?, and (d) Do patterns of diversity differ between the three diversity facets and/or the two ectoparasite taxa? We found that the α-diversity of infracommunities was strongly affected by host species, biome, and sampling period. The highest proportion of infracommunity diversity in both taxa was associated with the interaction between either host species and biome or host species and sampling period. Infracommunities of both taxa within, as well as between, host species, biomes, and sampling periods were characterized by the clustering of species, traits and lineages. The patterns of the effects of host species, biome, and sampling period on infracommunity diversity were congruent among the three diversity facets in both fleas and mites. We conclude that the assembly patterns in ectoparasite infracommunities mirror those characteristics of component and compound communities.

Harrison's rule scales up to entire parasite assemblages but is determined by environmental factors.

Harrison's rule states that parasite body size and the body size of their hosts tend to be positively correlated. After it was proposed a century ago, a number of studies have investigated this trend, but the support level has varied greatly between parasite/host associations. Moreover, while the rule has been tested at the individual species level, we still lack knowledge on whether Harrison's rule holds at the scale of parasite and host communities. Here, we mapped flea (parasites) and rodent (hosts) body sizes across Mongolia and asked whether Harrison's rule holds for parasite/host assemblages (i.e. whether a parasite's average body size in a locality is positively correlated with its host's average body size). In addition, we attempted to disentangle complex relationships between flea size, host size and environmental factors by testing alternative hypotheses for the determinants of fleas' body size variation. We gathered occurrence data for fleas and rodents from 2,370 sites across Mongolia, constructed incidence matrices for both taxa and calculated the average body sizes of fleas and their hosts over half-degree cells. Then, we applied a path analysis, accounting for spatial autocorrelation, trying to disentangle the drivers of the correlation between parasite and host body sizes. We found a strong positive correlation between average flea and host size across assemblages. Surprisingly though, we found that environmental factors simultaneously affected the body sizes of both fleas and hosts in the same direction, leading to a most likely deceptive correlation between parasite and host size across assemblages. We suggest that environmental factors may, to a great extent, reflect the environmental conditions inside the hosts' burrows where fleas develop and attain their adult body size, thus influencing their larval growth. Similarly, rodent body size is strongly influenced by air temperature, in the direction predicted by Bergmann's rule. If our findings are valid in other host-parasite associations, this may explain the dissenting results of both support and lack thereof for Harrison's rule.

Drivers of compositional turnover are related to species' commonness in flea assemblages from four biogeographic realms: zeta diversity and multi-site generalised dissimilarity modelling.

We investigated drivers of species turnover in fleas parasitic on small mammals in four biogeographic realms using novel methodology (zeta diversity, and Multi-Site Generalised Dissimilarity Modelling). We asked whether (i) flea turnover was better explained by host turnover or environmental variables; (ii) different factors drive the turnover of rare and widespread fleas; (iii) the factors affecting the turnover of rare or widespread fleas differ between realms; and (iv) environmental variables drive flea turnover directly or via their effects on hosts. Dissimilarity in host species composition was the most important factor affecting flea turnover in all realms. In the Afrotropics, the Nearctic, and the Neotropics, this was true mainly for rare species, whereas the zeta diversity of the Palearctic hosts exerted a strong effect on the turnover of both rare and widespread fleas. Dissimilarity in temperature contributed the most to the turnover of rare fleas in the Neotropics and the Palearctic, whereas the turnover of widespread species in these realms was strongly affected by dissimilarity in precipitation. In the Nearctic, dissimilarity in precipitation or temperature mostly affected the turnover of rare fleas or common species, respectively. In the Afrotropics, dissimilarity in the Normalised Difference Vegetation Index and temperature affected the turnover of all species, independently of their level of commonness, while dissimilarity in rainfall was important for the turnover of rare fleas. The responses of flea assemblages to environmental factors represented a combination of direct responses and responses mediated via effects on host turnover. We conclude that host turnover is a more important factor than environmental dissimilarity in its effect on flea species turnover. However, the relative effects of host composition and environment, as well as those of temperature, precipitation and the amount of vegetation, on flea turnover differ (i) between realms and (ii) between rare and common fleas.

The effects of environment, hosts and space on compositional, phylogenetic and functional beta-diversity in two taxa of arthropod ectoparasites.

We studied the effects of variation in environmental, host-associated and spatial factors on variation in compositional, phylogenetic/taxonomic and functional facets of beta-diversity in fleas and gamasid mites parasitic on small mammals and asked whether (a) the importance of these factors as drivers of beta-diversity differs among its multiple facets and (b) the effects of variation in environment, hosts and space on beta-diversity variation differ between the two ectoparasite taxa. To understand the relative effects of each group of predictors, we used a distance-based redundancy analysis and variation partitioning. The greatest portions of variation in the compositional beta-diversity of fleas were equally explained by host-associated and spatial predictors, whereas variation in host species composition contributed the most to variation in the compositional beta-diversity of mites. Variation in the phylogenetic (i.e. based on phylogenetic tree) beta-diversity of fleas was mainly due to variation in the phylogenetic composition of host communities, while the taxonomic (i.e. based on Linnean taxonomy) beta-diversity of mites was influenced by environmental variation. Unique contributions of spatial and environmental variation explained most of the variation in functional beta-diversity and its species replacement (= turnover) component (i.e. beta-diversity explained by replacement of species alone) in fleas and mites, respectively. Variation in the richness difference component (i.e. beta-diversity explained by species loss/gain alone) of functional beta-diversity was mainly affected by either variation in the functional composition of host assemblages (fleas) or its joint action with environmental variables (mites). We conclude that the pattern of the relative effects of environmental, host-associated and spatial factors on beta-diversity is context-dependent and may differ among different facets of beta-diversity, among different beta-diversity components and also among taxa dependent on biological affinities.