Rtapas: An R Package to Assess Cophylogenetic Signal between Two Evolutionary Histories.

Cophylogeny represents a framework to understand how ecological and evolutionary process influence lineage diversification. The recently developed algorithm Random Tanglegram Partitions provides a directly interpretable statistic to quantify the strength of cophylogenetic signal and incorporates phylogenetic uncertainty into its estimation, and maps onto a tanglegram the contribution to cophylogenetic signal of individual host-symbiont associations. We introduce Rtapas, an R package to perform Random Tanglegram Partitions. Rtapas applies a given global-fit method to random partial tanglegrams of a fixed size to identify the associations, terminals, and internal nodes that maximize phylogenetic congruence. This new package extends the original implementation with a new algorithm that examines the contribution to phylogenetic incongruence of each host-symbiont association and adds ParaFit, a method designed to test for topological congruence between two phylogenies, to the list of global-fit methods than can be applied. Rtapas facilitates and speeds up cophylogenetic analysis, as it can handle large phylogenies (100+ terminals) in affordable computational time as illustrated with two real-world examples. Rtapas can particularly cater for the need for causal inference in cophylogeny in two domains: (i) Analysis of complex and intricate host-symbiont evolutionary histories and (ii) assessment of topological (in)congruence between phylogenies produced with different DNA markers and specifically identify subsets of loci for phylogenetic analysis that are most likely to reflect gene-tree evolutionary histories. [Cophylogeny; cophylogenetic signal; gene tree incongruence; phylogenetic congruence; phylogenomics.].

Dissimilarity in flea and host assemblages and their interaction networks along a spatial distance gradient: different patterns revealed by different network dissimilarity metrics.

We investigated the distance-decay pattern (an increase in dissimilarity with increasing geographic distance) in regional assemblages of fleas and their small mammalian hosts, as well as their interaction networks, in four biogeographic realms. Dissimilarity of assemblages (ßtotal) was partitioned into species richness differences (ßrich) and species replacement (ßrepl) components. Dissimilarity of networks was assessed using two metrics: (a) whole network dissimilarity (ßWN) partitioned into species replacement (ßST) and interaction rewiring (ßOS) components and (b) D statistics, measuring dissimilarity in the pure structure of the networks, without using information on species identities and calculated for hosts-shared-by-fleas networks (Dh) and fleas-shared-by-hosts networks (Df). We asked whether the distance-decay pattern (a) occurs among interactor assemblages or their interaction networks; (b) depends on the network dissimilarity metric used; and (c) differs between realms. The ßtotal and ßrepl of flea and host assemblages increased with distance in all realms except for host assemblages in the Afrotropics. ßrich for flea and host assemblages increased with distance in the Nearctic only. In networks, ßWN and ßST demonstrated a distance-decay pattern, whereas ßOS was mainly spatially invariant except in the Neotropics. Correlations of Dh or Df and geographic distance were mostly non-significant. We conclude that investigations of dissimilarity in interaction networks should include both types of dissimilarity metrics (those that consider partner identities and those that consider the pure structure of networks). This will allow elucidating the predictability of some facets of network dissimilarity and the unpredictability of other facets.

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.

Functional similarity affects similarity in partner composition in flea-mammal networks.

Functional signal in an interaction network is a phenomenon in which species resembling each other in their traits interact with similar partners. We tested the functional signal concept in realm-specific and regional flea-host networks from four biogeographic realms and asked whether the species composition of (a) host spectra and (b) flea assemblages is similar between functionally similar flea and host species, respectively. Analogously to testing for phylogenetic signal, we applied Mantel tests to investigate the correlation between flea or host functional distances calculated from functional dendrograms and dissimilarities in sets of interacting partners. In all realm-specific networks, functionally similar fleas tended to exploit similar hosts often belonging to the same genus, whereas functionally similar hosts tended to harbour similar fleas, again often belonging to the same genus. The strength of realm-specific functional signals and the frequency of detecting a significant functional signal in the regional networks differed between realms. The frequency of detecting a significant functional signal in the regional networks correlated positively with the network size for fleas and with the number of hosts in a network for hosts. A functional signal in the regional networks was more frequently found for hosts than for fleas. We discuss the mechanisms behind the functional signal in both fleas and their hosts, relate geographic functional signal patterns to the historic biogeography of fleas and conclude that functional signals in the species composition of host spectra for fleas and of flea assemblages for hosts result from the interplay of evolutionary and ecological processes.

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.

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.

Nestedness and beta diversity of gastrointestinal helminth communities in common warthogs, Phacochoerus africanus (Suidae), at 2 localities in South Africa.

Few studies have investigated the ecological interactions between wild species of Suidae and their parasites, leaving our knowledge concerning this host­parasite system fragmented. In the present study, we applied network studies to analyse community nestedness in helminth assemblages of common warthogs, Phacochoerus africanus (Gmelin) (Suidae). Helminth data were compiled from 95 warthogs, including young and adult males and females, from 2 different conservation areas in Mpumalanga and Limpopo Provinces, South Africa, collected monthly over a period of 1 year each. The aim was to study the effect of host sex, age and season of sampling on the structure of helminth infracommunities harboured by the warthogs and to search for non-random structural patterns in the warthog­helminth interaction networks. Furthermore, we investigated the influence of a warthog's age, sex and season of sampling on beta diversity and dark diversity of their helminth infracommunities. Lastly, we asked whether the effects of host sex, age and sampling season on helminth communities differed between the 2 localities. We found that helminth communities of warthogs were nested and host­parasite interactions were influenced by all 3 factors as well as combinations thereof. However, the resulting patterns differed at the 2 localities, indicating that local environmental processes are important drivers of community structure.

Ectoparasites associated with the Bushveld gerbil (Gerbilliscus leucogaster) and the role of the host and habitat in shaping ectoparasite diversity and infestations.

Rodents are known hosts for various ectoparasite taxa such as fleas, lice, ticks and mites. South Africa is recognized for its animal diversity, yet little is published about the parasite diversity associated with wild rodent species. By focusing on a wildlife-human/domestic animal interface, the study aims to record ectoparasite diversity and levels of infestations of the Bushveld gerbil, Gerbilliscus leucogaster, and to establish the relationship between ectoparasite infestation parameters and host- and habitat factors. Rodents (n = 127) were trapped in 2 habitat types (natural and agricultural) during 2014­2020. More than 6500 individuals of 32 epifaunistic species represented by 21 genera and belonging to 5 taxonomic groups (fleas, sucking lice, ticks, mesostigmatan mites and trombiculid mites) were collected. Mesostigmatan mites and lice were the most abundant and fleas and mesostigmatan mites the most prevalent groups. Flea and mesostigmatan mite numbers and mesostigmatan mite species richness was significantly higher on reproductively active male than female rodents. Only ticks were significantly associated with habitat type, with significantly higher tick numbers and more tick species on rodents in the natural compared to the agricultural habitat. We conclude that the level of infestation by ectoparasites closely associated with the host (fleas and mites) was affected by host-associated factors, while infestation by ectoparasite that spend most of their life in the external environment (ticks) was affected by habitat type.

Documenting the microbiome diversity and distribution in selected fleas from South Africa with an emphasis on the cat flea, Ctenocephalides f. felis.

The factors that influence parasite associated bacterial microbial diversity and the geographic distributions of bacteria are not fully understood. In an effort to gain a deeper understanding of the relationship between the bacterial diversity of Ctenocephalides fleas and host species and the external environment, we conducted a metagenetic analysis of 107 flea samples collected from 8 distinct sampling sites in South Africa. Pooled DNA samples mostly comprising of 2 or 3 individuals sampled from the same host, and belonging to the same genetic cluster, were sequenced using the Ion PGM™ Hi-Q™ Kit and the Ion 316™ Chip v2. Differences were detected in the microbiome compositions between Ctenocephalides felis, Ctenocephalides canis and Ctenocephalides connatus. Although based on a small sample, C. connatus occurring on wildlife harboured a higher bacterial richness when compared to C. felis on domestic animals. Intraspecific differences in the microbial OTU diversity were detected within C. f. felis that occurred on domestic cats and dogs. Different genetic lineages of C. f. felis were similar in microbial compositions but some differences exist in the presence or absence of rare bacteria. Rickettsia and Bartonella OTU's identified in South African cat fleas differ from those identified in the USA and Australia. Intraspecific microbial compositions also differ across geographic sampling sites. Generalized dissimilarity modelling showed that temperature and humidity are potentially important environmental factors explaining the pattern obtained.

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.

The species composition of local flea assemblages at a small scale in two South American regions is predominantly driven by niche-based mechanisms.

We applied a step-down factor analysis (SDFA) and multi-site generalised dissimilarity modelling (MS-GDM) to local flea communities harboured by small mammals (i.e., collected at small sampling sites over a short time period) in two South American regions (Patagonia and the Northwestern Argentina) with the aim of understanding whether these communities were assembled via niche-based or dispersal-based processes. The SDFA allows us to determine whether clusters of flea assemblages across different types of climates, vegetation and soils can be distinguished (suggesting niche-based assembly). MS-GDM allows us to determine whether a substantial proportion of the variation in flea species turnover is explained by specific climate-associated, vegetation-associated and soil-associated variables (indicating niche-based assembly) or host turnover (indicating dispersal-based assembly). Mapping of assemblages on climate, vegetation and soil maps, according to their loadings on axis 1 or axis 2 of the SDFA, did not provide clear-cut results. Clusters of similar loadings could be recognized within some, but not other, climate, vegetation and soil types. However, MS-GDM demonstrated that the effect of environmental variables (especially air temperature) on flea compositional turnover was much stronger than that of host turnover, indicating the predominance of niche-based processes in local community assembly. A comparison of our results with those on the mechanisms that drive species assembly in regional communities allows us to conclude that local and regional communities result from the joint action of niche-based and dispersal-based processes, with the former more important at a smaller spatial scale and the latter at a larger spatial scale.

Similarity in ixodid tick communities harboured by wildlife and livestock in the Albany Thicket Biome of South Africa.

South Africa has a diverse fauna of ixodid tick species, several of which are of medical or veterinary importance. Elucidation of which host or environmental characteristics determine the composition and abundance of tick assemblages is, therefore, highly relevant to disease management and wildlife conservation efforts. Here, we analysed the similarity in ixodid tick assemblages of three wildlife and three livestock species in a natural (Great Fish River Reserve) and anthropogenic (an adjacent farm) habitat in South Africa. We compared tick infracommunities of three wild host species between the reserve and the farm; of three wild host species within the reserve and of wild and livestock species on the farm (considering body size). Hosts examined in this study harboured the adults and immature stages of 11 tick species belonging to five ixodid genera. Notably, several tick species of South African wildlife have successfully made the switch to livestock and thus both wild hosts and livestock now contribute to the pool of maintenance hosts for these ticks as well as their associated pathogens. This is an important consideration when translocating wildlife or livestock as part of farming or conservation activities.

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.

Host phylogeny and ecology, but not host physiology, are the main drivers of (dis)similarity between the host spectra of fleas: application of a novel ordination approach to regional assemblages from four continents.

We investigated the patterns of phylogenetic and functional (dis)similarity in the species composition of host spectra between co-habitating generalist flea species in regional assemblages from four continents (Europe, Asia, North America and Africa) using a recently developed ordination approach (Double Similarity Principal Component Analysis). From the functional perspective, we considered physiological [body mass and basal metabolic rate (BMR)] and ecological (shelter depth and complexity) host traits. We asked (a) whether host phylogeny, physiology or ecology is the main driver of (dis)similarities between flea host spectra and (b) whether the patterns of phylogenetic and functional (dis)similarity in host spectra vary between flea assemblages from different continents. Phylogenetic similarity between the host spectra was highest in Africa, lowest in North America and moderate in Europe and Asia. In each assemblage, phylogenetic clusters of hosts dominating in the host spectra could be distinguished. The functional similarity between the host spectra of co-occurring fleas was low for shelter structure in all assemblages and much higher for body mass and BMR in three of the four assemblages (except North America). We conclude that host phylogeny and shelter structure are the main drivers of (dis)similarity between the host spectra of co-habitating fleas. However, the effects of these factors on the patterns of (dis)similarity varied across continents.

Fitness consequences of host colonization in two generalist fleas: Context-dependency and the effect of spatial co-occurrence.

We studied the fitness consequences of colonizing a novel host by experimental lines of fleas (Synosternus cleopatrae and Xenopsylla ramesis) maintained for 18-22 generations on the principal or novel (sympatric or allopatric) hosts via number, developmental success and size of the offspring of the fleas exploiting these hosts. We asked whether (a) fitness on non-principal hosts increases after prolonged maintenance; (b) the colonization success depends on the spatial co-occurrence of a flea and a host and (c) colonization of a novel host is accompanied by a decreased ability to exploit an original host. The ability of fleas to colonize novel hosts differed between species, with S. cleopatrae, but not X. ramesis, increasing its offspring production on novel hosts. Spatial co-occurrence did not affect colonization success. Maintenance on an alternative host was not accompanied by decreased adaptation to the original host. When fleas returned to the original host, their reproductive output was higher than that of their ancestors. We conclude that the success of colonizing a novel host is (a) context-dependent and varies between flea and host species and (b) not accompanied by the loss of ability to exploit an ancestral host but may lead to an increase in this ability.

Temporal variation of metacommunity structure in arthropod ectoparasites harboured by small mammals: the effects of scale and climatic fluctuations.

We applied the elements of metacommunity structure (EMS) approach and studied the temporal dynamics of metacommunity structure in arthropod ectoparasites (fleas, gamasid mites and ixodid ticks) harboured by six small mammalian hosts sampled for three decades in the same locality in Western Siberia at three hierarchical scales (inframetacommunities, component metacommunities and a compound metacommunity). All metacommunities were positively coherent. Inframetacommunity structures varied across sampling periods in all host species. The main structural pattern in an inframetacommunity of the same host varied across sampling times but was mostly characterized by clumped species distributions (Clementsian, Gleasonian and their quasi-versions). Component metacommunities in five of the six host species were characterized by either a Clementsian or a quasi-Clementsian distribution. In four of the six host species, this pattern was driven by mite distribution. The temporal structure of compound metacommunity was characterized by a Clementsian pattern. In contrast to the majority of component metacommunities, this pattern was driven by fleas, whereas the temporal structure of gamasid mite compound metacommunities demonstrated a Gleasonian distribution. The temporal gradient in infracommunity composition was not associated with temporal changes in either air temperature or precipitation, whereas the precipitation gradient was positively correlated with the structure of component (in five host species) and compound metacommunities. In conclusion, the best-fit metacommunity structure of ectoparasites varies temporally due to temporal changes in distribution patterns that can be associated with year-to-year climatic variation, affecting both hosts and parasites.

Impact of host sex and age on the diversity of endoparasites and structure of individual-based host-parasite networks in nyalas (Tragelaphus angasii Angas) from three game reserves in KwaZulu-Natal province, South Africa.

In recent years, numerous studies have examined the effect of host sex and age on the structure of parasite communities in several host taxa under various environmental conditions and in different geographic regions. However, the influence of such factors on the structure of host-parasite networks has received less attention, and remarkably few studies have been carried out on large terrestrial mammals. In this study, we investigated the effects of host age and sex on the parasite infra- and component communities of nyalas (Tragelaphus angasii) and on the structure of individual-based nyala-endoparasite networks. We also aimed to evaluate to what extent these effects vary spatially and if they are mediated by conservation management. Based on a data set of internal macroparasites of 74 nyalas from three game reserves in KwaZulu-Natal province, we found that host age strongly influenced parasite community structure as well as the structure of parasite-nyala networks, whereas host sex played a minor role. However, the effects of both host sex and age were mediated by environmental conditions and thus led to different patterns at the three localities. Our findings highlight that host-parasite communities from different localities should not be pooled when conducting host-parasite network and community studies as this may bias results and mask patterns that are typical for a given locality.

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.

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.