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.

Comparative phylogeography of parasitic Laelaps mites contribute new insights into the specialist-generalist variation hypothesis (SGVH)

Background

The specialist-generalist variation hypothesis (SGVH) in parasites suggests that, due to patchiness in habitat (host availability), specialist species will show more subdivided population structure when compared to generalist species. In addition, since specialist species are more prone to local stochastic extinction events with their hosts, they will show lower levels of intraspecific genetic diversity when compared to more generalist.

Results

To test the wider applicability of the SGVH we compared 337 cytochrome oxidase I mitochondrial DNA and 268 nuclear tropomyosin DNA sequenced fragments derived from two co-distributed Laelaps mite species and compared the data to 294 COI mtDNA sequences derived from the respective hosts Rhabdomys dilectus, R. bechuanae, Mastomys coucha and M. natalensis. In support of the SGVH, the generalist L. muricola was characterized by a high mtDNA haplotypic diversity of 0.97 (±0.00) and a low level of population differentiation (mtDNA F_st = 0.56, p < 0.05; nuDNA F_st = 0.33, P < 0.05) while the specialist L. giganteus was overall characterized by a lower haplotypic diversity of 0.77 (±0.03) and comparatively higher levels of population differentiation (mtDNA F_st = 0.87, P < 0.05; nuDNA F_st = 0.48, P < 0.05). When the two specialist L. giganteus lineages, which occur on two different Rhabdomys species, are respectively compared to the generalist parasite, L. muricola, the SGVH is not fully supported. One of the specialist L. giganteus species occurring on R. dilectus shows similar low levels of population differentiation (mtDNA F_st = 0.53, P < 0.05; nuDNA F_st = 0.12, P < 0.05) than that found for the generalist L. muricola. This finding can be correlated to differences in host dispersal: R. bechuanae populations are characterized by a differentiated mtDNA F_st of 0.79 (P < 0.05) while R. dilectus populations are less structured with a mtDNA F_st = 0.18 (P < 0.05).

Conclusions

These findings suggest that in ectoparasites, host specificity and the vagility of the host are both important drivers for parasite dispersal. It is proposed that the SGHV hypothesis should also incorporate reference to host dispersal since in our case only the specialist species who occur on less mobile hosts showed more subdivided population structure when compared to generalist species.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1245-7) contains supplementary material, which is available to authorized users.

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.

Host range and distribution of small mammal fleas in South Africa, with a focus on species of medical and veterinary importance

The host range and distribution of flea species on rodents and insectivores across multiple vegetation types in South Africa were investigated. Habitat suitability for flea species considered as important vectors of disease in humans and domestic animals was modelled. Data originated from fleas that were recovered from small mammals captured at 29 localities during 2009-2013 and published literature searched for flea records. Climate-based predictor variables, widely used in arthropod vector distribution, were selected and habitat suitability modelled for 10 flea vector species. A total of 2469 flea individuals representing 33 species and subspecies were collected from 1185 small mammals. Ten of each of the flea and rodent species are plague vectors and reservoirs, respectively. Multiple novel flea-host associations and locality records were noted. Three vector species were recorded from insectivores. Geographic distributions of flea species ranged from broad, across-biome distributions to narrower distributions within one or two biomes. Habitat suitability models performed excellently for the majority of flea vectors and identified regions of summer and all-year rainfall as representing suitable habitats for most vector species. Current knowledge of vector and disease ecology can benefit from similar sampling approaches that will be important not only for South Africa, but also for the sub-region.

The influence of life history characteristics on flea (Siphonaptera) species distribution models

BACKGROUND

Ectoparasites exhibit pronounced variation in life history characteristics such as time spent on the host and host range. Since contemporary species distribution (SD) modelling does not account for differences in life history, the accuracy of predictions of current and future species' ranges could differ significantly between life history groups.

RESULTS

SD model performance was compared between 21 flea species that differ in microhabitat preferences and level of host specificity. Distribution models generally performed well, with no significant differences in model performance based on either microhabitat preferences or host specificity. However, the relative importance of predictor variables was significantly related to host specificity, with the distribution of host-opportunistic fleas strongly limited by thermal conditions and host-specific fleas more associated with conditions that restrict their hosts' distribution. The importance of temperature was even more pronounced when considering microhabitat preference, with the distribution of fur fleas being strongly limited by thermal conditions and nest fleas more associated with variables that affect microclimatic conditions in the host nest.

CONCLUSIONS

Contemporary SD modelling, that includes climate and landscape variables, is a valuable tool to study the biogeography and future distributions of fleas and other parasites taxa. However, consideration of life history characteristics is cautioned as species may be differentially sensitive to environmental conditions.