Abstract
Rickettsia felis (Alphaproteobacteria: Rickettsiales) is the causative agent of an emerging flea-borne rickettsiosis with worldwide occurrence. Originally described from the cat flea, Ctenocephalides felis, recent reports have identified R. felis from other flea species, as well as other insects and ticks. This diverse host range for R. felis may indicate an underlying genetic variability associated with host-specific strains. Accordingly, to determine a potential genetic basis for host specialization, we sequenced the genome of R. felis str. LSU-Lb, which is an obligate mutualist of the parthenogenic booklouse Liposcelis bostrychophila (Insecta: Psocoptera). We also sequenced the genome of R. felis str. LSU, the second genome sequence for cat flea-associated strains (cf. R. felis str. URRWXCal2), which are presumably facultative parasites of fleas. Phylogenomics analysis revealed R. felis str. LSU-Lb diverged from the flea-associated strains. Unexpectedly, R. felis str. LSU was found to be divergent from R. felis str. URRWXCal2, despite sharing similar hosts. Although all three R. felis genomes contain the pRF plasmid, R. felis str. LSU-Lb carries an additional unique plasmid, pLbaR (plasmid of L. bostrychophila associated Rickettsia), nearly half of which encodes a unique 23-gene integrative conjugative element. Remarkably, pLbaR also encodes a repeats-in-toxin-like type I secretion system and associated toxin, heretofore unknown from other Rickettsiales genomes, which likely originated from lateral gene transfer with another obligate intracellular parasite of arthropods, Cardinium (Bacteroidetes). Collectively, our study reveals unexpected genomic diversity across three R. felis strains and identifies several diversifying factors that differentiate facultative parasites of fleas from obligate mutualists of booklice.