Fleas (Siphonaptera) are Cretaceous, and evolved with Theria

The first mitogenome of the subfamily Stenoponiinae (Siphonaptera: Ctenophthalmidae) and implications for its phylogenetic position

Fleas are the most important insect vectors that parasitize warm-blooded animals and are known vectors of zoonotic pathogens. A recent study showed that Stenoponia polyspina parasitizing Eospalax baileyi in Zoige County have carried Bartonella spp. and Spotted fever group Rickettsia (SFGR). Accurate identification and differentiation of fleas are essential for prevention and control of zoonotic pathogens. To understand phylogenetic relationship of the subfamily Stenoponiinae, we described morphological characteristics of S. polyspina and sequenced its mitogenome with 14,933 bp in size and high A + T content (~ 79%). The S. polyspina mitogenome retained the ancestral pattern of mitochondrial gene arrangement of arthropods without rearrangement. The start codons of 13 protein-coding genes (PCGs) are traditional ATN and the stop codons are TAA or TAG. Anticodon loops of all tRNA genes were 7 bp except for trnL2 and trnD had anticodon loops with 9 bp and the abnormal anticodon loops may be associated with frameshifting mutation. Genetic distance and Ka/Ks ratios indicated that all 13 PCGs of S. polyspina were subjected to purifying selection, with cox1 at the slowest rate and atp8 at the fastest rate. The mitogenomes of 24 species representing 7 families in the order Siphonaptera were selected to reconstruct phylogenetic tree based on concatenated nucleotide sequences of two datasets (PCGRNA matrix and PCG12RNA matrix) using Bayesian inference (BI) and Maximum likelihood (ML) methods. Phylogenetic tree supported that the superfamilies Ceratophylloidea, Vermipsylloidea, Pulicoidea were monophyletic and the superfamily Hystrichopsylloidea was paraphyletic. The family Ctenophthalmidae was monophyletic in PCGRNA-ML (codon partition) and paraphyletic in the remain trees. S. polyspina belongs to the subfamily Stenoponiinae was closely more related to the subfamily Rhadinopsyllinae. This paper explored phylogenetic position of diverse clades within the order Siphonaptera based on morphological and mitogenome data of S. polyspina. Our research enriched NCBI database of the order Siphonaptera.

Ectoparasites of wild rodents in forest sites invaded and uninvaded by Maesopsis eminii in Amani nature forest reserve, Tanzania

Parasites are important component of communities in a forest ecosystem with profound effects on trophic interactions such as food web. Modification of the forest structure (e.g. changes in species composition and abundance of key species) can have a strong impact on the occurrence, diversity, and abundance of parasites, with subsequent repercussions for ecosystem functioning. In this study, we compared the occurrence and abundance of wild rodents' ectoparasites from forest sites invaded and uninvaded by an invasive tree, Maesopsis eminii in Amani Nature Forest Reserve, Tanzania. Three large plots (40 m × 100 m) were randomly established in each forest sites invaded and uninvaded by M. eminii. In each plot, 50 Sherman traps were systematically placed at 10 m interval for capturing wild rodents through a capture-mark-recapture technique. Wilcox rank sum test was used to compare for differences in the abundance of infested rodents and ectoparasites between the invaded and uninvaded forest sites. A total of 297 individual rodents were captured and screened for ectoparasites, including 174 rodents from uninvaded forest site and 123 rodents from invaded forest site. The number of infested rodents were significantly (W = 8592, P < 0.001) greater in uninvaded forest site (66.27%) than in the invaded forest site (36.2%). Furthermore, a significant greater number of Echinolaelaps echidninus (W = 1849, P < 0.01) and Dinopsyllus ellobius (W = 2800.5, P < 0.05) ectoparasites were found in uninvaded as compared to the invaded forest sites. The results of this study suggest that the invasion and dominance by, M. eminii in Amani Nature Reserve has created unfavorable conditions for rodents and ectoparasites and therefore impacting the diversity and function of the forest ecosystem. We recommend prevention of further introduction of the M. eminii outside their natural range and mitigating the impact of the established M. eminii in Amani Forest Nature Reserve.

Nudiviruses in free-living and parasitic arthropods: evolutionary taxonomy

The nudiviruses (family: Nudiviridae) are large double-stranded DNA (dsDNA) viruses that infect insects and crustaceans, and have most recently been identified from ectoparasitic members (fleas and lice). This virus family was created in 2014 and has since been expanded via the discovery of multiple novel viral candidates or accepted members, sparking the need for a new taxonomic and evolutionary overview. Using current information (including data from public databases), we construct a new comprehensive phylogeny, encompassing 49 different nudiviruses. We use this novel phylogeny to propose a new taxonomic structure of the Nudiviridae by suggesting two new viral genera (Zetanudivirus and Etanudivirus), from ectoparasitic lice. We detail novel emerging relationships between nudiviruses and their hosts, considering their evolutionary history and ecological role.

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.

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.

Detection of Rickettsia raoultii in Vermipsylla alakurt-Like Fleas of Sheep in Northwestern China

INTRODUCTION

To date, a total of 2574 validated flea species have been discovered. Vermipsyllidae is a family of fleas that comprises at least eight species. Vermipsylla is a genus of the family Vermipsyllidae within the order Siphonaptera of fleas. Here a novel Vermipsylla species was described, and rickettsial agent was also detected in it.

METHODS

A total of 128 fleas were collected directly from 260 pastured sheep in China. Of these, eight representative fleas (four males and four females) were identified by key morphological features. Meanwhile, 120 flea DNAs, including six flea samples for molecular taxonomy, were subjected to Rickettsia spp. DNA detection. The molecular identity of fleas was determined by amplification and sequenmce analysis of four genetic markers (the 28S rDNA genes, the 18S rDNA genes, the mitochondrial cytochrome c oxidase subunit I and subunit II). In addition, five Rickettsia-specific gene fragments were used to identify the species of the rickettsial agents. The amplified products were sequenced and phylogenetically analyzed.

RESULTS

The morphological characteristics of the flea species identified in this study were similar to Vermipsylla alakurt, but presented difference in hair number of the metepimeron, the third tergum, the genitals and the tibiae of hind leg. The 18S rDNA, 28S rDNA and COII genetic markers from fleas showed the highest identity to those of V. alakurt, shared 98.45% (954/969), 95.81% (892/931) and 85.86% (571/665) similarities, respectively. However, the COI sequence showed the highest identity to that of Dorcadia ioffi with 88.48% (576/651) similarity. Rickettsia raoutii tested positive in 14.17% (17/120) flea DNA samples.

CONCLUSION

Our study reports the detection of R. raoultii in V. alakurt-like fleas infesting sheep in China.

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.

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.

Host specificity of flea parasites of mammals from the Andean Biogeographic Region

Host specificity of fleas affects their biodiversity that plays a major role in determining the potential transmission routes by pathogens through vertebrate hosts, including humans. In the Biogeographic Andean region, numerous systematic and ecological studies have been conducted, revealing a high diversity of flea taxa of mammals and the presence of pathogenic organisms transmitted by fleas; however, the degree of preference with which each flea species associates with a mammal host remains poorly understood in this region. Herein, host specificity in mammal fleas from the Andean region was analysed. We employed the number of host species for each flea species and the index of host specificity STD *. Following the literature, 144 species and 13 subspecies of fleas (31 genera and 10 families) have been described in the Andean biogeographic region; 76 taxa are endemic to this region. To carry out the analyses of host specificity, we considered 1759 records of fleas collected from 124 species and 59 genera of wild and domestic mammals, mostly rodent species (85.9%). Our results indicate that typical Andean fleas are genus or family host specific (mostly STD * less than 3.0). More diverse mammal hosts are parasitized by more diverse flea genera and families and these hosts are phylogenetically related. Otherwise, these hosts are associated with different flea lineages, suggesting the interaction of ecological and evolutionary mechanisms (host-switching, ecological adaptations and co-evolutionary alternation). The fields of disease ecology and One Health are considering the host specificity of arthropod vectors as an important point to understand the mechanisms of emergence and re-emergence of diseases. Our results allow us to estimate the risk of diseases involving fleas in the Andean region.

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.

Arthropod Ectoparasites of Two Rodent Species Occurring in Varied Elevations on Tanzania’s Second Highest Mountain

Climate change causes organisms, including species that act as parasite reservoirs and vectors, to shift their distribution to higher altitudes, affecting wildlife infestation patterns. We studied how ectoparasite distributions varied with altitude using two rodent species, Montemys delectorum and Rhabdomys dilectus, at different elevations (1500–3500 m). The ectoparasites infesting the two rodent species were influenced by the host sex, species, and temperature. We expected host density to predict parasite infestation patterns, because hosts in higher densities should have more parasites due to increased contact between individuals. However, temperature, not host density, affected ectoparasite distribution. Since temperatures decrease with elevation, parasite prevalences and abundances were lower at higher elevations, highlighting that the cold conditions at higher elevations limit reproduction and development—this shows that higher elevation zones are ideal for conservation. The rodents and ectoparasite species described in this study have been reported as vectors of diseases of medical and veterinary importance, necessitating precautions. Moreover, Mount Meru is a refuge for a number of endemic and threatened species on the IUCN Red List. Thus, the parasitic infection can also be an additional risk to these critical species as well as biodiversity in general. Therefore, our study lays the groundwork for future wildlife disease surveillance and biodiversity conservation management actions. The study found a previously uncharacterized mite species in the Mesostigmata group that was previously known to be a parasite of honeybees. Further investigations may shed light into the role of this mite species on Mount Meru.

New Insights into the Taxonomy of Malacopsylloidea Superfamily (Siphonaptera) Based on Morphological, Molecular and Phylogenetic Characterization of Phthiropsylla agenoris (Malacopsyllidae) and Polygenis (Polygenis) rimatus (Rhopalopsyllidae)

From a phylogenetic point of view, the Malacopsyllidae family and the Rhopalopsillidae family (comprising Parapsyllinae and Rhopalopsyllinae subfamilies) have been traditionally classified within the Malacopsylloidea superfamily, mostly restricted to South America. The phylogenetic relationships and taxonomic status of Malacopsyllidae and Rhopalopsillidae have never been assessed since no molecular loci of Malacopsyllidae have been sequenced by any authors, and the phylogeny provided so far was not based on any sort of formal quantitative analysis of flea morphology. Based on these precedents, the objective of this study was to carry out a comparative phylogenetic, molecular and morphological study of two different species belonging to each family, Phthiropsylla agenoris (Malacopsylla) and Polygenis (Polygenis) rimatus (Rhopalopsyllidae, Rhopalopsyllinae). In this study, we demonstrated the usefulness of several morphological features as diagnostic characters to differentiate between P. (P.) rimatus and P. agenoris. Using molecular and phylogenetic data, we easily discriminated between the two taxa (P. agenoris and P. (P.) rimatus) by comparing both nuclear and mitochondrial markers. This fact proves the usefulness of ITS2, EF1−α, cox1, cytb and cox2 as molecular diagnostic markers to characterize and identify different Siphonaptera taxa. Additionally, the phylogenetic results confirm, for the first time, the monophyly of the Malacopsyllidae family and suggest a clear paraphyletic position of the Paraspsyllinae subfamily and, consequently, the Rhopalopsyllidae family.

Thermal niche partitioning and phenology of Nearctic and Palearctic flea (Siphonaptera) communities on rodents (Mammalia: Rodentia) from five ecoregions

Seasonality of fleas (Siphonaptera) may be due to species competition, prompting the idea that flea species partition temperature, along with correlated variables such as moisture (thermal-niche partitioning hypothesis). I compared the fleas of five rodent-flea communities described from the literature for thermal-niche optima by fitting non-linear LRF (Lobry-Rosso-Flandrois) curves to examine whether flea species in a community show distinct, partitioned thermal niches. LRF curves estimate physiological parameters of temperature minimum, optimum, maximum, and maximum abundance, and facilitate comparison between species by summarizing seasonal data. Flea-communities were on Nearctic Southern flying squirrel (Glaucomys volans volans), Richardson's ground-squirrel (Urocitellus richardsonii), North American deer-mouse (Peromyscus maniculatus), and Palearctic Midday jird (Meriones meridianus), and Wagner's gerbil (Dipodillus dasyurus). Flea communities appeared to show seasonality consistent with thermal-niche partitioning. Several flea families and genera had characteristic thermal niches: Ceratophyllidae had broad tolerance to extreme temperature, Leptopsyllidae (one species in this study) to cold, and Pulicidae to hot. In contrast, at the local, species level, climatic speciation could be significant in flea diversification. Non-competition hypotheses (environmental filtering, neutrality) require testing, too. Thermal-niche partitioning may increase flea species richness on hosts and could occur in other insect and plant communities. Implications for biodiversity conservation and disease ecology under global warming are wide-ranging.

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.

The Angiosperm Terrestrial Revolution and the origins of modern biodiversity

Biodiversity today has the unusual property that 85% of plant and animal species live on land rather than in the sea, and half of these live in tropical rainforests. An explosive boost to terrestrial diversity occurred from c. 100-50 million years ago, the Late Cretaceous and early Palaeogene. During this interval, the Earth-life system on land was reset, and the biosphere expanded to a new level of productivity, enhancing the capacity and species diversity of terrestrial environments. This boost in terrestrial biodiversity coincided with innovations in flowering plant biology and evolutionary ecology, including their flowers and efficiencies in reproduction; coevolution with animals, especially pollinators and herbivores; photosynthetic capacities; adaptability; and ability to modify habitats. The rise of angiosperms triggered a macroecological revolution on land and drove modern biodiversity in a secular, prolonged shift to new, high levels, a series of processes we name here the Angiosperm Terrestrial Revolution.

Phylogenetic structure of specialization: A new approach that integrates partner availability and phylogenetic diversity to quantify biotic specialization in ecological networks

Biotic specialization holds information about the assembly, evolution, and stability of biological communities. Partner availabilities can play an important role in enabling species interactions, where uneven partner availabilities can bias estimates of biotic specialization when using phylogenetic diversity indices. It is therefore important to account for partner availability when characterizing biotic specialization using phylogenies. We developed an index, phylogenetic structure of specialization (PSS), that avoids bias from uneven partner availabilities by uncoupling the null models for interaction frequency and phylogenetic distance. We incorporate the deviation between observed and random interaction frequencies as weights into the calculation of partner phylogenetic α‐diversity. To calculate the PSS index, we then compare observed partner phylogenetic α‐diversity to a null distribution generated by randomizing phylogenetic distances among the same number of partners. PSS quantifies the phylogenetic structure (i.e., clustered, overdispersed, or random) of the partners of a focal species. We show with simulations that the PSS index is not correlated with network properties, which allows comparisons across multiple systems. We also implemented PSS on empirical networks of host–parasite, avian seed‐dispersal, lichenized fungi–cyanobacteria, and hummingbird pollination interactions. Across these systems, a large proportion of taxa interact with phylogenetically random partners according to PSS, sometimes to a larger extent than detected with an existing method that does not account for partner availability. We also found that many taxa interact with phylogenetically clustered partners, while taxa with overdispersed partners were rare. We argue that species with phylogenetically overdispersed partners have often been misinterpreted as generalists when they should be considered specialists. Our results highlight the important role of randomness in shaping interaction networks, even in highly intimate symbioses, and provide a much‐needed quantitative framework to assess the role that evolutionary history and symbiotic specialization play in shaping patterns of biodiversity. PSS is available as an R package at https://github.com/cjpardodelahoz/pss.

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.

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.

The Acari Hypothesis, II: Interspecies Operability of Pattern Recognition Receptors

Hypersensitivity to galactose-α-1,3-galactose (α-gal) is an informative example of a pathologic IgE-mediated process. By way of their saliva, ticks are able to sensitize humans to tick dietary elements that express α-gal. Mites, which along with ticks constitute the phyletic subclass Acari, feed on proteinaceous foodstuffs that represent most, if not all, human allergens. Given: (1) the gross nature of the pathophysiological reactions of allergy, especially anaphylaxis, (2) the allergenicity of acarian foodstuffs, and (3) the relatedness of ticks and mites, it has been hypothesized that human-acarian interactions are cardinal to the pathogenesis of allergy. In this report, a means by which such interactions contribute to that pathogenesis is proposed.

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.

Combination of nuclear and mitochondrial markers as a useful tool to identify Ctenophthalmus species and subspecies (Siphonaptera: Ctenophthalmidae)

Ctenophthalmus is considered the largest genus within the Order Siphonaptera. From a morphological point of view, only males of this genus can be identified at species and subspecies levels using morphological keys, whereas there are no morphological criteria in order to classify females at these taxonomical levels. Furthermore, the amount of available molecular and phylogenetic data for this genus is quite scarce so far. The main objective of this work was to assess the utility of the combination of nuclear and mitochondrial markers with respect to their ability to differentiate among different subspecies within the Ctenophthalmus genus. With this purpose, we carried out a comparative morphological and molecular study of three different subspecies (Ctenophthalmus baeticus arvernus, Ctenophthalmus nobilis dobyi, and Ctenophthalmus andorrensis catalaniensis) in order to clarify and discuss its taxonomic status. In addition, our study complemented the molecular data previously provided for Ctenophthalmus baeticus boisseauorum and Ctenophthalmus apertus allani subspecies. We sequenced five different molecular markers: EF1-α, ITS1, ITS2, cox1, and cytb. Our results confirmed that morphological data by themselves are not able to discriminate among Ctenophthalmus female taxa; however, the combination of the nuclear marker EF1-α together with mtDNA markers cytb and cox1 constituted a useful taxonomical and phylogenetic tool to solve this issue. Based on these results, we consider that the use of this molecular approach should be gradually used within Ctenophthalmus genus in order to complement its classical taxonomy and clarifying the complex taxonomy of other congeneric species of fleas.

Molecular Identification, Transcriptome Sequencing and Functional Annotation of Pulex irritans

PURPOSE

Pulex irritans are vectors of various zoonotic pathogens. However, molecular studies on P. irritans and flea-borne diseases are limited due to the lack of molecular data. This study aimed to conduct transcriptome sequencing, functional annotation, and pathogen analysis of P. irritans.

METHODS

Fleas collected from a dog were identified morphologically and molecularly. RNA was extracted for transcriptome sequencing and functional annotation. Open reading frames (ORFs) of unigenes were confirmed by employing bioinformatics strategies, and maximum likelihood (ML) trees were reconstructed based on the highly expressed genes of ejaculation globulin-specific 3-like protein, salivary protein, and actin for phylogenetic relationship analysis.

RESULTS

The obtained mitochondrial 16S rRNA gene sequences showed 99.71% of similarity with P. irritans obtained from GenBank database. Transcriptome sequencing generated 74,412 unigenes, of which 53,211 were functionally annotated. A total of 195 unigenes were assigned to fleas, of which 69 contained complete ORFs. Phylogenetic trees of both ejaculatory globulin and salivary protein genes demonstrated that P. irritans first clustered with Pulicidae sp., indicating the reliability of transcriptome data. It is noteworthy that 1070 unigenes were assigned to Hymenolepis microstoma and Dipylidium caninum, of which 62 contained complete ORFs. The phylogenetic tree of the actin gene showed that the unigenes had closer relationships with Echinococcus sp., suggesting the role of P. irritans as intermediate hosts of tapeworms.

CONCLUSION

The results of this study provide the possibility for functional exploration of important genes and lay foundations for the prevention and control of P. irritans and flea-borne diseases.

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.

Functional and phylogenetic uniqueness of helminth and flea assemblages of two South African rodents

The loss of a particular species from a community may have different effects on its functioning, depending on the presence or absence of functionally similar or phylogenetically close species in that community (redundancy). Redundancy is thus defined as the fraction of species diversity not expressed by functional or phylogenetic diversity. We assessed functional and phylogenetic alpha- and beta-redundancy in helminth and flea assemblages of two species of South African rodents, Rhabdomys dilectus and Rhabdomys pumilio, using community uniqueness as the inverse indicator of redundancy. We asked whether patterns of functional and phylogenetic alpha- and beta-uniqueness differed between (i) parasite groups (endo- versus ectoparasites), (ii) host species within parasite groups, and (iii) biomes within host species. We found differences between the two hosts in the functional and phylogenetic alpha-uniqueness (but not beta-uniqueness) of flea, but not helminth, assemblages. Significant correlations between the alpha-uniqueness of parasite assemblages and the total parasite prevalence were found only for phylogenetic uniqueness and only in helminths. Pairwise site-by-site dissimilarities in uniqueness (beta-uniqueness) and pairwise dissimilarity in prevalence were significantly associated (positively) in helminths but not in fleas. A between-biome difference in functional (but not phylogenetic) alpha-uniqueness was found in both helminth and flea assemblages harboured by R. pumilio. We conclude that the resilience of parasite assemblages in terms of the effect on hosts depends not only on their transmission strategy but also on traits of host species and environmental factors.

Bats are key hosts in the radiation of mammal-associated Bartonella bacteria

Bats are notorious reservoirs of several zoonotic diseases and may be uniquely tolerant of infection among mammals. Broad sampling has revealed the importance of bats in the diversification and spread of viruses and eukaryotes to other animal hosts. Vector-borne bacteria of the genus Bartonella are prevalent and diverse in mammals globally and recent surveys have revealed numerous Bartonella lineages in bats. We assembled a sequence database of Bartonella strains, consisting of nine genetic loci from 209 previously characterized Bartonella lineages and 121 new cultured isolates from bats, and used these data to perform a comprehensive phylogenetic analysis of the Bartonella genus. This analysis included estimation of divergence dates using a molecular clock and ancestral reconstruction of host associations and geography. We estimate that Bartonella began infecting mammals 62 million years ago near the Cretaceous-Paleogene boundary. Additionally, the radiation of particular Bartonella clades correlate strongly to the timing of diversification and biogeography of mammalian hosts. Bats were inferred to be the ancestral hosts of all mammal-associated Bartonella and appear to be responsible for the early geographic expansion of the genus. We conclude that bats have had a deep influence on the evolutionary radiation of Bartonella bacteria and their spread to other mammalian orders. These results support a 'bat seeding' hypothesis that could explain similar evolutionary patterns in other mammalian parasite taxa. Application of such phylogenetic tools as we have used to other taxa may reveal the general importance of bats in the ancient diversification of mammalian parasites.

Molecular detection of Rickettsia spp. in ticks and fleas collected from rescued hedgehogs (Erinaceus europaeus) in Portugal

Hedgehogs (e.g., Erinaceus europaeus, E. roumanicus) are wild mammals that frequently are observed near residential areas. The aim of this study was to investigate ticks and fleas collected from European hedgehogs in Portugal and to evaluate the prevalence of Rickettsia in those ectoparasites. Ticks and fleas were identified by morphological and molecular methods, and molecular detection by PCR and genotypic characterization of Rickettsia spp. was performed targeting ompB, ompA and gltA gene fragments. In total, 1892 ticks and 213 fleas were collected from 33 rescued European hedgehogs captured in seven districts of the north and centre of Portugal. Two tick species were identified - Rhipicephalus sanguineus accounted for 91 % (n = 1719) of the total ticks collected and 9 % (n = 173) were Ixodes hexagonus. All fleas were identified as Archaeopsylla erinacei. Regarding pathogen detection, Rickettsia massiliae DNA was found in 22 of the 212 tested Rh. sanguineus. None of the 48 I. hexagonus tested showed to be positive for rickettsiae. Rickettsia asembonensis DNA was identified in 55 A. erinacei fleas tested (n = 117). These results show that European hedgehogs are exposed to R. massiliae transmitted by ticks and to R. asembonensis via fleas suggesting that these mammals might be involved in the natural transmission cycle of these Rickettsia species. This study is the first report of R. asembonensis in fleas in Portugal.

Molecular Differentiation of Four Species of Oropsylla (Siphonaptera: Ceratophyllidae) Using PCR-Based Single Strand Conformation Polymorphism Analyses and DNA Sequencing

It is often difficult to distinguish morphologically between closely related species of fleas (Siphonaptera). Morphological identification of fleas often requires microscopic examination of internal structures in specimens cleared using caustic solutions. This process degrades DNA and/or inhibits DNA extraction from specimens, which limits molecular-based studies on individual fleas and their microbiomes. Our objective was to distinguish between Oropsylla rupestris (Jordan), Oropsylla tuberculata (Baker), Oropsylla bruneri (Baker), and Oropsylla labis (Jordan & Rothschild) (Ceratophyllidae) using PCR-based single strand conformation polymorphism (SSCP) analyses and DNA sequencing. A 446 bp region of the nuclear 28S ribosomal RNA (rRNA) gene was used as the genetic marker. The results obtained for 36 reference specimens (i.e., fleas that were morphologically identified to species) revealed no intraspecific variation in DNA sequence, whereas the DNA sequences of the four species of Oropsylla differed from one another at two to six nucleotide positions. Each flea species also had a unique SSCP banding pattern. SSCP analyses were then used to identify another 84 fleas that had not been identified morphologically. DNA sequencing data confirmed the species identity of fleas subjected to SSCP. This demonstrates that PCR-SSCP combined with DNA sequencing of the 28S rRNA gene is a very effective approach for the delineation of four closely related species of flea.

Ctenophthalmus baeticus boisseauorum (Beaucournu, 1968) and Ctenophthalmus apertus allani (Smit, 1955) (Siphonaptera: Ctenophthalmidae) as synonymous taxa: morphometric, phylogenetic, and molecular characterization

The family Ctenophthalmidae (Order Siphonaptera) has been considered as a 'catchall' for a wide range of divergent taxa showing a paraphyletic origin. In turn, Ctenophthalmus sp. (Ctenophthalmidae) includes 300 valid described taxa. Within this genus, males are easily distinguishable basing on the size, shape, and chaetotaxy of their genitalia; however, females show slight morphological differences with each other. The main objective of this work was to carry out a comparative morphometric, phylogenetic, and molecular study of two different subspecies: Ctenophthalmus baeticus boisseauorum and Ctenophthalmus apertus allani in order to clarify and discuss its taxonomic status. From a morphological and biometrical point of view, we found clear differences between modified abdominal segments of males of both subspecies and slight differences in the margin of sternum VII of all female specimens which did not correspond with molecular and phylogenetic results based on four different molecular markers (Internal Transcribed Spacer 1 and 2 of ribosomal DNA, and the partial cytochrome c oxidase subunit 1 and cytochrome b of mitochondrial DNA). Thus, we observed a phenotypic plasticity between both subspecies, which did not correspond with a real genotypic variability nor different environmental or ecological conditions. Basing on these results, we could consider that there are no solid arguments to consider these two 'morphosubspecies' as two different taxa. We propose that C. b. boisseauorum should be considered as a junior synonym of C. a. allani.

New taxonomic and evolutionary insights relevant to the cat flea, Ctenocephalides felis: A geographic perspective

The cat flea, Ctenocephalides felis, is an obligate haematophagous ectoparasite of wildlife and domestic cats and dogs worldwide. Since cat fleas can affect the health of humans and their pets, an uncertain taxonomy of this taxon can greatly inhibit pest and disease management. To address the evolution and taxonomy of the cat flea, we set out to determine 1) how many genetically distinct taxa exist, 2) whether there is morphological support for the genetically distinct taxa, and 3) the role of host range and paleoclimatic events in speciation. We collected a total of 3352 fleas sampled from 576 domestic cats and dogs as well as 10 wildlife species across 30 localities in South Africa. A total of three flea genera, five species, and three of the currently recognized cat flea subspecies, C. f. damarensis, C. f. strongylus and C. f. felis were obtained. Geometric morphometric analyses on head shape were performed on 68 female and 107 male cat flea individuals. Principal component analysis demonstrated large overlap in head shape variation between C. f. strongylus and C. f. felis, rendering this character not useful for phylogenetic inferences. DNA was extracted from 188 Ctenocephalides spp. and mitochondrial COII and nuclear EF1-α sequences were generated. Bayesian and Maximum Likelihood analyses as well as a TCS parsimony haplotype network of the mitochondrial DNA confirmed the presence of three well supported monophyletic clades. These assemblages did not fully corroborate the existence of the three C. felis subspecies. A single well-supported molecular clade included only C. f. damarensis morphotypes that were mostly collected from wildlife. The recognition of this subspecies as a distinct taxon was further corroborated by sequence distances and also the number of plantar spiniform bristles on fore-tarsi V in males. Despite the overall lack of support for the recognition of C. f. damarensis and C. f. strongylus, a geographic trend was visible whereby one genetic lineage corresponded to the western dryer hot subregion, whereas the other was found throughout the region. Bayesian dating suggested that these two clades diverged during the early Pliocene (4.18 mya), a date that corresponds well with the establishment of a dry hot climate in the west of southern Africa. If so, the off-host environment, particularly temperature and humidity, are important factors to consider in the evolution of the cat flea. The present study rejects recent assertions that the three cat flea subspecies are valid entities and rather point to a situation where more sampling is required before the taxonomic status of C. f. damarensis can be resolved.

Feeding performance on a novel host: no adaptation over generations and differential patterns in two flea species

To model the colonization of a novel host by fleas, Synosternus cleopatrae and Xenopsylla ramesis, we established experimental lines maintained for 15 generations on a principal or a novel host (either co-occurring with a flea or not). We compared the blood meal size and the energy expended for digestion by fleas from the 15th generation of each line on these hosts between hosts within a line and between lines within a host asking (a) whether fleas adapt to a novel host (increased blood consumption/decreased energy expended for digestion); (b) if yes, whether this adaptation leads to the loss of ability to exploit an original host, and (c) whether the success of adaptation to a novel host depends on its ecological co-occurrence with a flea. The blood consumption and digestion energetics of fleas fed on the principal host differed from those on other hosts. The effect of the principal host on feeding performance differed between fleas, with S. cleopatrae consuming less blood and expending more energy for digestion on the principal than on any other host, whereas the opposite was true for X. ramesis. No changes in feeding performance on a novel host over generations were found. We propose several explanations for the lack of adaptation to a novel host over time. We explain the poor performance of S. cleopatrae on its principal host via its immune response mounting pattern. We argue that the principal host of a parasite is not necessarily the host on which the parasite demonstrates the best performance.

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.

Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry as a useful tool for the rapid identification of wild flea vectors preserved in alcohol

Flea identification is a significant issue because some species are considered as important vectors of several human pathogens that have emerged or re-emerged recently, such as Bartonella henselae (Rhizobiales: Bartonellaceae) and Rickettsia felis (Rickettsiales: Rickettsiaceae). Matrix-assisted laser desorption/ionization time-of-flight mass spectrometry (MALDI-TOF MS) has been evaluated in recent years for the identification of multicellular organisms, including arthropods. A preliminary study corroborated the usefulness of this technique for the rapid identification of fleas, creating a preliminary database containing the spectra of five species of flea. However, longterm flea preservation in ethanol did not appear to be an adequate method of storage in the context of specimen identification by MALDI-TOF MS profiling. The goal of the present work was to assess the performance of MALDI-TOF MS in the identification of seven flea species [Ctenocephalides felis (Siphonaptera: Pulicidae), Ctenocephalides canis, Pulex irritans (Siphonaptera: Pulicidae), Archaeopsylla erinacei (Siphonaptera: Pulicidae), Leptopsylla taschenbergi (Siphonaptera: Ceratophyllidae), Stenoponia tripectinata (Siphonaptera: Stenoponiidae) and Nosopsyllus fasciatus (Siphonaptera: Ceratophyllidae)] collected in the field and stored in ethanol for different periods of time. The results confirmed that MALDI-TOF MS can be used for the identification of wild fleas stored in ethanol. Furthermore, this technique was able to discriminate not only different flea genera, but also the two congeneric species C. felis and C. canis.

Origin, evolution, phylogeny and taxonomy of Pulex irritans

The human flea Pulex irritans Linnaeus, 1758 (Siphonaptera: Pulicidae) is one of the most studied species together with the cat flea Ctenocephalides felis Bouché, 1835, because they have a cosmopolitan distribution and are closely related to humans. The present study aimed to carry out a comparative morphometric and molecular study of two different populations of P. irritans (Spain and Argentina). Accordingly, internal transcribed spacer (ITS)1 and ITS2 of rDNA and the partial cytochrome c oxidase subunit 1 (cox1) and cytochrome b (cytb) mtDNA genes of these taxa were sequenced. Furthermore, the taxonomy, origin, evolution and phylogeny of P. irritans was assessed. The morphometric data obtained did not show significant differences between P. irritans specimens from Spain and Argentina, even when these two populations were collected from different hosts; however, there was a considerable degree of molecular divergence between both populations based on nuclear and mitochondrial markers. Thus, it is proposed that P. irritans, in contrast with other generalist fleas, maintains a certain degree of morphological similarity, at least between Western Palearctic and Neotropical areas. Furthermore, two well defined geographical genetic lineages within the P. irritans species are indicated, suggesting the existence of two cryptic species that could be discriminated by a polymerase chain reaction-linked restriction fragment length polymorphism.

Out-of-Africa, human-mediated dispersal of the common cat flea, Ctenocephalides felis: The hitchhiker's guide to world domination

The cat flea (Ctenocephalides felis) is the most common parasite of domestic cats and dogs worldwide. Due to the morphological ambiguity of C. felis and a lack of - particularly largescale - phylogenetic data, we do not know whether global C. felis populations are morphologically and genetically conserved, or whether human-mediated migration of domestic cats and dogs has resulted in homogenous global populations. To determine the ancestral origin of the species and to understand the level of global pervasion of the cat flea and related taxa, our study aimed to document the distribution and phylogenetic relationships of Ctenocephalides fleas found on cats and dogs worldwide. We investigated the potential drivers behind the establishment of regional cat flea populations using a global collection of fleas from cats and dogs across six continents. We morphologically and molecularly evaluated six out of the 14 known taxa comprising genus Ctenocephalides, including the four original C. felis subspecies (Ctenocephalides felis felis, Ctenocephalides felis strongylus, Ctenocephalides felis orientis and Ctenocephalides felis damarensis), the cosmopolitan species Ctenocephalides canis and the African species Ctenocephalides connatus. We confirm the ubiquity of the cat flea, representing 85% of all fleas collected (4357/5123). Using a multigene approach combining two mitochondrial (cox1 and cox2) and two nuclear (Histone H3 and EF-1α) gene markers, as well as a cox1 survey of 516 fleas across 56 countries, we demonstrate out-of-Africa origins for the genus Ctenocephalides and high levels of genetic diversity within C. felis. We define four bioclimatically limited C. felis clusters (Temperate, Tropical I, Tropical II and African) using maximum entropy modelling. This study defines the global distribution, African origin and phylogenetic relationships of global Ctenocephalides fleas, whilst resolving the taxonomy of the C. felis subspecies and related taxa. We show that humans have inadvertently precipitated the expansion of C. felis throughout the world, promoting diverse population structure and bioclimatic plasticity. By demonstrating the link between the global cat flea communities and their affinity for specific bioclimatic niches, we reveal the drivers behind the establishment and success of the cat flea as a global parasite.

Fleas from domestic dogs and rodents in Rwanda carry Rickettsia asembonensis and Bartonella tribocorum

Fleas (Siphonaptera) are ubiquitous blood-sucking parasites that transmit a range of vector-borne pathogens. The present study examined rodents (n = 29) and domestic dogs (n = 7) living in the vicinity of the Volcanoes National Park, Rwanda, for fleas, identified flea species from these hosts, and detected Bartonella (Rhizobiales: Bartonellaceae) and Rickettsia (Rickettsiales: Rickettsiaceae) DNA. The most frequently encountered flea on rodents was Xenopsylla brasiliensis (Siphonaptera: Pulicidae). In addition, Ctenophthalmus (Ethioctenophthalmus) calceatus cabirus (Siphonaptera: Hystrichopsyllidae) and Ctenocephalides felis strongylus (Siphonaptera: Pulicidae) were determined using morphology and sequencing of the cytochrome c oxidase subunit I and cytochrome c oxidase subunit II genes (cox1 and cox2, respectively). Bartonella tribocorum DNA was detected in X. brasiliensis and Rickettsia asembonensis DNA (a Rickettsia felis-like organism) was detected in C. felis strongylus. The present work complements studies that clarify the distributions of flea-borne pathogens and potential role of fleas in disease transmission in sub-Saharan Africa. In the context of high-density housing in central sub-Saharan Africa, the detection of B. tribocorum and R. asembonensis highlights the need for surveillance in both rural and urban areas to identify likely reservoirs.

Morphological, biometrical and molecular characterization of Archaeopsylla erinacei (Bouché, 1835)

In the present work, we carried out a morphological, biometrical and molecular study of the species Archaeopsylla erinacei (Bouché, 1835) and their subspecies: Archaeopsylla erinacei erinacei (Bouché, 1835) and Archaeopsylla erinacei maura (Jordan & Rothschild, 1912) isolated from hedgehogs (Erinaceus europaeus) from different geographical regions (Seville and Corse). We have found morphological differences in females of A. erinacei from the same geographical origin that did not correspond with molecular differences. We suggest that some morphological characters traditionally used to discriminate females of both subspecies should be revised as well as we set the total length of the spermatheca as a valid criterion in order to discriminate between both subspecies. The Internal Transcribed Spacers 1 and 2 (ITS1, ITS2) and partial 18S rRNA gene, and partial cytochrome c-oxidase 1 (cox1) and cytochrome b (cytb) mtDNA gene sequences were determined to clarify the taxonomic status of these taxa and to assess intra-specific and intra-population similarity. In addition, a phylogenetic analysis with other species of fleas using Bayesian and Maximum Likelihood analysis was performed. All molecular markers used, except 18S, showed molecular differences between populations corresponding with geographical origins. Thus, based on the phylogenetic and molecular study of two nuclear markers (ITS1, ITS2) and two mitochondrial markers (cox1 and cytb), as well as concatenated sequences of both subspecies, we reported the existence of two geographical genetic lineages in A. erinacei corresponding with two different subspecies: A. e. erinacei (Corse, France) and A. e. maura (Seville, Spain), that could be discriminated by polymerase chain reaction-linked random-fragment-length polymorphism.

Ectoparasitism and infections in the exoskeletons of large fossil cingulates

Studies on paleopathological alterations in fossil vertebrates, including damages caused by infections and ectoparasites, are important because they are potential sources of paleoecological information. Analyzing exoskeleton material (isolated osteoderms, carapace and caudal tube fragments) from fossil cingulates of the Brazilian Quaternary Megafauna, we identified damages that were attributed to attacks by fleas and dermic infections. The former were compatible with alterations produced by one species of flea of the genus Tunga, which generates well-delimited circular perforations with a patterned distribution along the carapace; the latter were attributable to pathogenic microorganisms, likely bacteria or fungi that removed the ornamentation of osteoderms and, in certain cases, generated craters or pittings. Certain bone alterations observed in this study represent the first record of flea attack and pitting in two species of large glyptodonts (Panochthus and Glyptotherium) and in a non-glyptodontid large cingulate (Pachyarmatherium) from the Quaternary of the Brazilian Intertropical Region. These new occurrences widen the geographic distribution of those diseases during the Cenozoic and provide more evidence for the co-evolutionary interaction between cingulates and parasites registered to date only for a small number of other extinct and extant species.

Body size and ecological traits in fleas parasitic on small mammals in the Palearctic: larger species attain higher abundance

We studied the relationships between body size and (a) abundance and (b) host specificity in fleas parasitic on small mammals (rodents and shrews) in the Palearctic taking into account the confounding effect of phylogeny. We tested these relationships both across 127 flea species and within separate phylogenetic clades, predicting higher abundance and lower host specificity (in terms of the number or diversity of hosts used by a flea) in smaller species. We also tested for the relationships between body size and abundance separately for species that spend most of their lives on a host's body (the "body" fleas) and species that spend most of their lives in a host's burrow or nest (the "nest" fleas). A significant phylogenetic signal in body size was detected across all fleas, as well as in five of six separate clades. Across all fleas and in majority of phylogenetic clades, mean abundance significantly increased with an increase in body size. The same pattern was found for both the "body" and the "nest" fleas, although the slope of the relationship appeared to be steeper in the former than in the latter. Neither measure of host specificity demonstrated a significant correlation with body size regardless of the subset of flea species analysed. We explain higher abundance attained by larger flea species by higher fecundity and/or competitive advantage upon smaller species at larval stage. We conclude that the macroecological patterns reported to date in parasites are far from being universal.

The latitudinal, but not the longitudinal, geographic range positions of haematophagous ectoparasites demonstrate historical signatures

We tested whether geographic range position of fleas parasitic on small mammals in the Palearctic is affected by environmental niche conservatism or geographic range conservatism by measuring phylogenetic signal in range centroids and boundaries. We predicted that stronger phylogenetic signal in latitudinal than longitudinal range positions would indicate the important role of niche conservatism as a driver of the evolution of fleas' geographic ranges. Phylogenetic signals in geographic range positions were measured across 120 species, as well as within five flea lineages (subfamily/family rank) of different evolutionary ages. To investigate the temporal pattern of the geographic range position's evolution, we fitted the phylogenetic patterns in the geographic coordinates of range centroids and border extremes to four models of trait evolution. We consistently detected significant phylogenetic signals in the latitudes of the range centroids and the northern range borders. The latitudes of the southern range borders and the longitudes of the eastern/western borders demonstrated phylogenetic signals less often, whereas no signal was found for the longitudes of the range centroids. The phylogenetic signal in range position was more pronounced in younger lineages. The phylogenetic signal indices mainly suggested the evolution of range positions according to the Brownian motion model, whereas the best fit was often provided by the Orstein-Uhlenbeck model. This contradiction forced us to invoke a parsimonious explanation that the phylogenetic signal in range positions results from the interplay between the footprint of the speciation pattern and limited dispersal from the ancestral ranges.

Phylogenetic heritability of geographic range size in haematophagous ectoparasites: time of divergence and variation among continents

To understand existence, patterns and mechanisms behind phylogenetic heritability in the geographic range size (GRS) of parasites, we measured phylogenetic signal (PS) in the sizes of both regional (within a region) and continental (within a continent) geographic ranges of fleas in five regions. We asked whether (a) GRS is phylogenetically heritable and (b) the manifestation of PS varies between regions. We also asked whether geographic variation in PS reflects the effects of the environment's spatiotemporal stability (e.g. glaciation disrupting geographic ranges) or is associated with time since divergence (accumulation differences among species over time). Support for the former hypothesis would be indicated by stronger PS in southern than in northern regions, whereas support for the latter hypothesis would be shown by stronger PS in regions with a large proportion of species belonging to the derived lineages than in regions with a large proportion of species belonging to the basal lineages. We detected significant PS in both regional and continental GRSs of fleas from Canada and in continental GRS of fleas from Mongolia. No PS was found in the GRS of fleas from Australia and Southern Africa. Venezuelan fleas demonstrated significant PS in regional GRS only. Local Indicators of Phylogenetic Association detected significant local positive autocorrelations of GRS in some clades even in regions in which PS has not been detected across the entire phylogeny. This was mainly characteristic of younger taxa.

High mitochondrial sequence divergence in synanthropic flea species (Insecta: Siphonaptera) from Europe and the Mediterranean

Background

Adult fleas are haematophagous ectoparasites of warm-blooded vertebrates, particularly mammals. Among them, the cat flea (Ctenocephalides felis) and the human flea (Pulex irritans) have high veterinary-medical significance, owing to their cosmopolitan distribution and role in the transmission of important vector-borne pathogens. While the taxonomy of Ct. felis has been investigated on a morphological basis during the past decades, its molecular-phylogenetic analyses have been only recently conducted. This study expands the knowledge on Ct. felis from hitherto less studied geographical regions, and includes representatives from additional flea families, less investigated with molecular approaches.

Methods

Fleas were collected in four countries of the Mediterranean Basin (Croatia, Italy, Malta and Israel), as well as in Hungary, from domestic and wild carnivores, rodents and humans. The DNA extracts of representative fleas (n = 148), belonging to ten species of eight genera, were used for PCR amplification of part of their cytochrome c oxidase subunits 1, 2 (cox1, cox2) and 18S rRNA genes, followed by sequencing and phylogenetic analyses.

Results

The majority (65.6%) of Ct. felis felis cox2 sequences showed 99.4–100% similarity to each other (haplogroup A), whereas those from Malta and Israel had 98.1–98.7% sequence similarity (haplogroup B), and a third sequence from Israel (haplotype C) had as low as 96.3% sequence similarity in comparison with a reference sequence from group “A”. Except for the shape of the head, no consistent morphological differences (e.g. in chaetotaxy) were found between haplogroups “A” and “C”. Haplotypes of Ct. canis were genetically more homogenous, with 99.6–100% sequence similarity to each other. However, when P. irritans collected from humans was compared to those from three species of wild carnivores, these only had 96.6% cox2 similarity. The mouse flea, Leptopsylla segnis and the northern rat flea, Nosopsyllus fasciatus were both shown to have haplotypes with low intraspecific cox2 similarities (96.2 and 94.4%, respectively). Taken together, differences between mitochondrial lineages within four flea species exceeded that observed between two Chaetopsylla spp. (which had 97.3% cox2 similarity). The topologies of cox1 and cox2 phylogenetic trees were in line with relevant sequence comparisons. Conversely, 18S rRNA gene analyses only resolved differences above the species level.

Conclusions

Ctenocephalides felis felis, P. irritans, L. segnis and N. fasciatus were shown to have such a high level of mitochondrial gene heterogeneity, that the uniformity of these flea taxa should be reconsidered. Although the present results are limited (especially in the case of L. segnis and N. fasciatus), there appears to be no geographical or host restriction, which could explain the divergence of these genetic lineages.

Electronic supplementary material

The online version of this article (10.1186/s13071-018-2798-4) contains supplementary material, which is available to authorized users.

The complete mitochondrial genome of the flea Ceratophyllus wui (Siphonaptera: Ceratophyllidae)

The complete mitogenome sequence of the flea, Ceratophyllus wui (Siphonaptera: Ceratophyllidae), was sequenced. The 18,081 bp long genome has the standard metazoan complement of 37 genes. These genes contain 13 protein-coding genes, 22 transfer RNA genes, two ribosomal RNA genes, and one control region. The nucleotide composition of the C. wui mitogenome was A: 37.72%, T: 38.99%, G: 9.51% and C: 13.78%. The A + T content is 76.71%, showing strong AT skew. Phylogenetic analysis indicated that Siphonapteran has sister relationship with Dipteran.

MOLECULAR SURVEILLANCE FOR BARTONELLA, BORRELIA, AND RICKETTSIA SPECIES IN TICKS FROM DESERT BIGHORN SHEEP ( OVIS CANADENSIS) AND MULE DEER ( ODOCOILEUS HEMIONUS) IN SOUTHERN CALIFORNIA, USA

:  Ticks (Acari: Ixodidae) were collected from 44 desert bighorn sheep ( Ovis canadensis) and 10 mule deer ( Odocoileus hemionus) in southern California, US during health inspections in 2015-16. Specimens were identified and screened by PCR analysis to determine the presence and prevalence of Bartonella, Borrelia, and Rickettsia species in ticks associated with these wild ruminants. None of the 60 Dermacentor hunteri and 15 Dermacentor albipictus ticks tested yielded positive PCR results. Additional tick specimens should be collected and tested to determine the prevalence of these confirmed or suspected tickborne pathogens within ruminant populations.

Relationships among different facets of host specificity in three taxa of haematophagous ectoparasites

Host specificity is a fundamental trait of a parasite species. Recently, multiple aspects of host specificity have been recognized, but the relationships between these facets are still poorly understood. Here, we studied pairwise relationships between basic, structural, phylogenetic and geographic host specificity in three taxa of haematophagous ectoparasitic arthropods that differ in tightness of their association with the host. We asked which metrics of host specificity are correlated within each parasite taxon and whether the patterns of the association between different facets of host specificity are similar among parasite taxa. Data on bat flies were taken from published surveys across the Neotropics while data on fleas and mites parasitic on small mammals were compiled from multiple published surveys across the Palaearctic. Basic, structural, phylogenetic and geographic specificity indices were calculated for 18 bat fly species recorded on 40 host species from 15 regions, 109 flea species recorded on 120 host species from 51 regions and 34 mite species recorded on 67 host species from 28 regions. Then, we tested for the correlation between any two measures of host specificity using model II regressions. We found that structural and basic specificity, as well as structural and geographic specificity, exhibited a positive association in all three taxa. However, basic and geographic specificity, as well as basic and phylogenetic specificity, were significantly positively associated in fleas but did not correlate in bat flies or mites. In addition, we found a significant negative association between structural and phylogenetic specificity in bat flies but no association in the remaining taxa. Moreover, geographic and phylogenetic specificity were not associated in any parasite taxon. Our results suggest that different facets of host specificity were shaped differently by natural selection in different taxa.

Revisiting the role of dissimilarity of host communities in driving dissimilarity of ectoparasite assemblages: non-linear vs linear approach

We revisited the role of dissimilarity of host assemblages in shaping dissimilarity of flea assemblages using a non-linear approach. Generalized dissimilarity models (GDMs) were applied using data from regional surveys of fleas parasitic on small mammals in four biogeographical realms. We compared (1) model fit, (2) the relative effects of host compositional and phylogenetic turnover and geographic distance on flea compositional and phylogenetic turnover, and (3) the rate of flea turnover along gradients of host turnover and geographic distance with those from earlier application of a linear approach. GDMs outperformed linear models in explaining variation in flea species turnover and host dissimilarity was the best predictor of flea dissimilarity, irrespective of scale. The shape of the relationships between flea compositional turnovers along host compositional turnover was similar in all realms, whereas turnover along geographic distance differed among realms. In contrast, the rate of flea phylogenetic turnover along gradients of host phylogenetic turnover differed among realms, whereas flea phylogenetic turnover did not depend on geographic distance in any realm. We demonstrated that a non-linear approach (a) explained spatial variation in parasite community composition better than and (b) revealed patterns that were obscured by earlier linear analyses.