Metagenomics of culture isolates and insect tissue illuminate the evolution ofWolbachia,RickettsiaandBartonellasymbionts inCtenocephalidesspp. fleas

While fleas are often perceived simply as a biting nuisance and cause of allergic dermatitis, they represent important disease vectors worldwide, especially for bacterial zoonoses such as plague (transmitted by rodent fleas) and some of the rickettsioses and bartonelloses. The cosmopolitan cat ( Ctenocephalides felis ) and dog ( C. canis ) fleas, as well as C. orientis (restricted to tropical and subtropical Asia), breed in human dwellings and are vectors of cat-scratch fever (caused by Bartonella spp.) and Rickettsia spp. of the so-called "transitional group". The latter includes R. felis (agent of flea-borne spotted fever) and R. asembonensis , an emerging pathogen. The relatively depauperate flea microbiome can also contain arthropod-specific endosymbionts, including a diverse range of Wolbachia strains. Here, we present circularized genome assemblies for two C. orientis -associated pathogens ( Bartonella clarridgeiae and R. asembonensis ) from Malaysia, a novel Wolbachia strain ( w Cori), and the C. orientis mitochondrion; all obtained by direct metagenomic sequencing of flea tissues. Moreover, we isolated two Wolbachia strains from Malaysian C. felis into tick cell culture and recovered circularized genome assemblies for both, one of which ( w CfeF) is newly sequenced. We demonstrate that the three Wolbachia strains are representatives of different major clades ("supergroups"), two of which appear to be flea-specific. These Wolbachia genomes exhibit unique combinations of features associated with reproductive parasitism or mutualism, including prophage WO, cytoplasmic incompatibility factors, and the biotin operon of obligate intracellular microbes. The first circularized assembly for R. asembonensis includes a plasmid with a markedly different structure and gene content compared to the published plasmid; moreover, this novel plasmid was also detected in cat flea metagenomes from the US. Analysis of loci under positive selection in the transitional group revealed genes involved in host-pathogen interactions that may facilitate host switching. Finally, the first B. clarridgeiae genome from Asia exhibited largescale genome stability compared to isolates from other continents, except for SNPs in regions predicted to mediate interactions with the vertebrate host. These findings highlight the paucity of data on the genomic diversity of Ctenocephalides -associated bacteria and raise questions regarding how interactions between members of the flea microbiome might influence vector competence.