Parasite species co-occurrence patterns on Peromyscus: Joint species distribution modelling

Hosts are often infested by multiple parasite species, but it is often unclear whether patterns of parasite co-occurrence are driven by parasite habitat requirements or parasite species interactions. Using data on infestation patterns of ectoparasitic arthropods (fleas, trombiculid mites, cuterebrid botflies) from deer mice (Peromyscus maniculatus), we analyzed species associations using joint species distribution modelling. We also experimentally removed a flea (Orchopeas leucopus) from a subset of deer mice to examine the effect on other common ectoparasite species. We found that the mite (Neotrombicula microti) and botfly (Cuterebra sp.) had a negative relationship that is likely a true biotic species interaction. The flea had a negative association with the mite and a positive association with the botfly species, both of which appeared to be influenced by host traits or parasite life-history traits. Furthermore, experimental removal of the flea did not have a significant effect on ectoparasite prevalence of another species. Overall, these findings suggest that complex parasite species associations can be present among multiple parasite taxa, and that aggregation is not always the rule for ectoparasite communities of small mammals.