Ectoparasites found on bats are known to contain important microbes. However, the viruses hosted by these obligate parasites are understudied. This has led to the near oversight of the potential role of these ectoparasites in virus maintenance and transmission from bats to other interacting species and the environment. Here, we sampled bat ectoparasites parasitizing a diverse selection of bat species in the families Rhinolophidae, Vespertilionidae, Megadermatidae, Hipposideridae and Pteropodidae in Yunnan Province, China. We show that the ectoparasite prevalence was generally higher in male compared to female bats. Most ectoparasites were found to fall within the Nycteribiidae, Spinturnicidae and Streblidae bat ectoparasite families. We subsequently applied a non-biased sequencing of libraries prepared from the pooled ectoparasites, followed by an in-silico virus-centric analysis of the resultant reads. We show that ectoparasites hosted by the sampled families of bats are found to carry, in addition to a diverse set of phages, vertebrate and insect viruses in the families Aliusviridae, Ascoviridae, Chuviridae, Circoviridae, Flaviviridae, Hepadnaviridae, Hepeviridae, Herpesviridae, Iridoviridae, Marseilleviridae, Nairoviridae, Orthomyxoviridae, Parvoviridae, Poxviridae, Reoviridae, Retroviridae, and Rhabdoviridae. We further report a partial Parvovirus VP1/VP2 gene and partial Poxvirus ubiquitin-like gene predicted by two independent next generation sequencing data analysis pipelines. This study describes the natural virome of bat ectoparasites, providing a platform for understanding the role these ectoparasites play in the maintenance and spread of viruses to other animals.

In recent years, the global pandemic of bat-associated pathogens has led to increasing attention on bat ectoparasites. Numerous studies have identified human-associated pathogens in Nycteribiidae, indicating their potential as vectors. In this study, the first complete sequencing of the mitochondrial genome of Nycteribia allotopa Speiser, 1901 was sequenced and analyzed. We also compared the mitochondrial sequences of N. allotopa with those available in the database for other Nycteribiidae species. The complete mitochondrial genome of N. allotopa was found to be 15,161 bp in size with an A + T content of 82.49%. Nucleotide polymorphism analysis of 13 protein-coding genes from five species of Nycteribiidae showed that nad6 exhibited the most significant variation, while cox1 was the most conserved. Furthermore, selection pressure analysis revealed cox1 to exhibit the strongest purifying selection, while atp8, nad2, nad4L, and nad5 showed slightly looser purifying selection. Pairwise genetic distances indicated that cox1 and cox2 were evolving comparatively slowly, whereas atp8, nad2, and nad6 were evolving comparatively quickly. Phylogenetic trees constructed using Bayesian inference and maximum likelihood methods demonstrated that all four families within the superfamily Hippoboscoidea clustered into one branch each, indicating their monophyly. N. allotopa was found to be most closely related to the same genus N. parvula. This study significantly enriches the molecular database for Nycteribiidae and provides invaluable reference data for future species identification, phylogenetic analysis, and exploration of their potential as vectors for human-associated pathogens.

In recent years, bat-associated pathogens, such as 2019 novel coronavirus, have been ravaging the world, and ectoparasites of bats have received increasing attention. Penicillidia jenynsii is a member of the family Nycteribiidae which is a group of specialized ectoparasites of bats. In this study, the complete mitochondrial genome of P. jenynsii was sequenced for the first time and a comprehensive phylogenetic analysis of the superfamily Hippoboscoidea was conducted. The complete mitochondrial genome of P. jenynsii is 16 165 base pairs (bp) in size, including 13 protein-coding genes (PCGs), 22 transfer RNA genes, 2 ribosomal RNA genes and 1 control region. The phylogenetic analysis based on 13 PCGs of the superfamily Hippoboscoidea known from the NCBI supported the monophyly of the family Nycteribiidae, and the family Nycteribiidae was a sister group with the family Streblidae. This study not only provided molecular data for the identification of P. jenynsii, but also provided a reference for the phylogenetic analysis of the superfamily Hippoboscoidea.

Species of the family Nycteribiidae are blood-sucking ectoparasites that parasitize bats. To further enrich the molecular data of species in the family Nycteribiidae, the complete mitochondrial genome of Nycteribia parvula was sequenced for the first time in this study. The complete mitochondrial genome of N. parvula is 16,060 base pairs (bp) in size, including 13 protein-coding genes (PCGs), 22 transfer RNA genes, two ribosomal RNA genes, and a control region. The nucleotide contents of A, T, G, and C are respectively 40.86%, 42.19%, 6.51%, and 10.44%. The phylogenetic analysis based on 13 PCGs supports the monophyly of the family Nycteribiidae, and N. parvula is the closest relative to Phthiridium szechuanum.

Phthiridium szechuanum is a bat surface parasite under the family Nycteribiidae that prefers to roost in the hair of bats to feed on their blood. In this study, the complete mitochondrial genome of P. szechuanum was studied for the first time using Illumina sequencing technology. The mitochondrial genome was 14,896 bp in size and was predicted to encode 37 genes including 13 protein-coding genes, 22 transfer RNA genes, and 2 ribosomal RNA genes. Phylogenetic trees were constructed using the IQ-TREE web server and phylogenetic analysis was performed using the maximum likelihood method, and P. szechuanum was found to be phylogenetically closest to Basilia ansifera. These data will provide a molecular biological approach to the species identification of P. szechuanum and provide a new reference for further studies on the population genetics and phylogeny of the family Nycteribiidae.

BACKGROUND: Understanding the interactions between bat flies and host bats offer us fundamental insights into the coevolutionary and ecological processes in host-parasite relationships. Here, we investigated the identities, host specificity, and patterns of host association of bat flies in a subtropical region in East Asia, which is an understudied region for bat fly research.
METHODS: We used both morphological characteristics and DNA barcoding to identify the bat fly species found on 11 cavernicolous bat species from five bat families inhabiting Hong Kong. We first determined the phylogenetic relationships among bat fly species. Then, we elucidated the patterns of bat-bat fly associations and calculated the host specificity of each bat fly species. Furthermore, we assembled the mitogenomes of three bat fly species from two families (Nycteribiidae and Streblidae) to contribute to the limited bat fly genetic resources available.
RESULTS: We examined 641 individuals of bat flies and found 20 species, of which many appeared to be new to science. Species of Nycteribiidae included five Nycteribia spp., three Penicillidia spp., two Phthiridium spp., one Basilia sp., and one species from a hitherto unknown genus, whereas Streblidae included Brachytarsina amboinensis, three Raymondia spp., and four additional Brachytarsina spp. Our bat-bat fly association network shows that certain closely related bat flies within Nycteribiidae and Streblidae only parasitized host bat species that are phylogenetically more closely related. For example, congenerics of Raymondia only parasitized hosts in Rhinolophus and Hipposideros, which are in two closely related families in Rhinolophoidea, but not other distantly related co-roosting species. A wide spectrum of host specificity of these bat fly species was also revealed, with some bat fly species being strictly monoxenous, e.g. nycteribiid Nycteribia sp. A, Phthiridium sp. A, and streblid Raymondia sp. A, while streblid B. amboinensis is polyxenous.
CONCLUSIONS: The bat fly diversity and specificity uncovered in this study have shed light on the complex bat-bat fly ecology in the region, but more bat-parasite association studies are still needed in East Asian regions like China as a huge number of unknown species likely exists. We highly recommend the use of DNA barcoding to support morphological identification to reveal accurate host-ectoparasite relationships for future studies.

Bartonella species are facultative intracellular bacteria and recognized worldwide as emerging zoonotic pathogens. Bartonella were isolated or identified by polymerase chain reaction (PCR) in bats and their ectoparasites worldwide, whereas the association between them was scarce, especially in Asia. In this study, a retrospective analysis with frozen samples was carried out to identify the genetic diversity of Bartonella in bats and their ectoparasites and to investigate the relationships of Bartonella carried by bats and their ectoparasites. Bats and their ectoparasites (bat flies and bat mites) were collected from caves in Hubei Province, Central China, from May 2018 to July 2020. Bartonella were screened by PCR amplification and sequencing of three genes (gltA, rpoB, and ftsZ). Bats, bat flies, and bat mites carried diverse novel Bartonella genotypes with a high prevalence. The sharing of some Bartonella genotypes between bats and bat flies or bat mites indicated a potential role of bat flies and bat mites as vectors of bartonellae, while the higher genetic diversity of Bartonella in bat flies than that in bats might be due to the vertical transmission of this bacterium in bat flies. Therefore, bat flies might also act as reservoirs of Bartonella. In addition, human-pathogenic B. mayotimonesis was identified in both bats and their ectoparasites, which expanded our knowledge on the geographic distribution of this bacterium and suggested a potential bat origin with bat flies and bat mites playing important roles in the maintenance and transmission of Bartonella.