Conjugative Transposons and Their Cargo Genes Vary across Natural Populations of Rickettsia buchneri Infecting the Tick Ixodes scapularis

Genomic and Proteomic Analyses of Bacterial Communities of Ixodes scapularis Ticks from Broome County, New York

The microbial communities of Ixodes scapularis, the primary vector of Lyme disease in North America, exhibit regional variations that may affect pathogen transmission and vector competence. We analyzed bacterial communities in I. scapularis ticks collected from Broome County, New York, using 16S rRNA gene sequencing (18 ticks) as well as mass spectrometry-based proteomics (36 ticks). According to the 16S rRNA analysis, the endosymbiont Rickettsia buchneri was the most abundant species, with significantly higher (p = 0.0011) abundance in females (54.76%) compared to males (31.15%). We detected Borreliella burgdorferi in 44.44% of ticks and Anaplasma phagocytophilum in two nymphs but in high relative abundances (12.73% and 46.46%). Male ticks exhibited higher bacterial diversity, although the community composition showed no significant clustering by sex or life stage. Co-occurrence analysis revealed negative associations between R. buchneri and Pseudomonas (p = 0.0245), but no associations with B. burgdorferi. Proteomic analysis identified 12 R. buchneri-specific proteins, additionally detecting the protozoan pathogen Babesia microti in 18.18% of females. These findings provide the first comprehensive characterization of I. scapularis microbiomes in the Southern Tier region of New York and suggest broader distribution of R. buchneri across tick life stages than previously recognized, with potential implications for pathogen transmission dynamics.

Development of a Recombinase-Mediated Cassette Exchange System for Gene Knockout and Expression of Non-Native Gene Sequences in Rickettsia

BACKGROUND/OBJECTIVES

Incidence of vector-borne diseases, including rickettsioses and anaplasmosis, has been increasing in many parts of the world. The obligate intracellular nature of rickettsial pathogens has hindered the development of robust genetic tools for the study of gene function and the identification of therapeutic targets. Transposon mutagenesis has contributed to recent progress in the identification of virulence factors in this important group of pathogens.

METHODS

Combining the efficiency of the himar1 transposon method with a recombinase-mediated system, we aimed to develop a genetic tool enabling the exchange of the transposon with a cassette encoding non-native sequences.

RESULTS

This approach was used in Rickettsia parkeri to insert a himar1 transposon encoding fluorescent protein and antibiotic resistance genes for visualization and selection, flanked by mismatched loxP sites to enable subsequent recombinase-mediated cassette exchange (RMCE). RMCE mediated by a plasmid-encoded Cre recombinase was then employed to replace the transposon with a different cassette containing alternate fluorescent and selection markers and epitopes of Anaplasma phagocytophilum antigens. The resulting genetically modified R. parkeri was trialed as a live-attenuated vaccine against spotted fever rickettsiosis and anaplasmosis in mice.

CONCLUSIONS

The use of this system provides a well-established and relatively efficient way of inserting non-native sequences into the rickettsial genome, with applications for the study of gene function and vaccine development.

Bacterial reprogramming of tick metabolism impacts vector fitness and susceptibility to infection

Arthropod-borne pathogens are responsible for hundreds of millions of infections in humans each year. The blacklegged tick, Ixodes scapularis, is the predominant arthropod vector in the United States and is responsible for transmitting several human pathogens, including the Lyme disease spirochete Borrelia burgdorferi and the obligate intracellular rickettsial bacterium Anaplasma phagocytophilum, which causes human granulocytic anaplasmosis. However, tick metabolic response to microbes and whether metabolite allocation occurs upon infection remain unknown. Here we investigated metabolic reprogramming in the tick ectoparasite I. scapularis and determined that the rickettsial bacterium A. phagocytophilum and the spirochete B. burgdorferi induced glycolysis in tick cells. Surprisingly, the endosymbiont Rickettsia buchneri had a minimal effect on bioenergetics. An unbiased metabolomics approach following A. phagocytophilum infection of tick cells showed alterations in carbohydrate, lipid, nucleotide and protein metabolism, including elevated levels of the pleiotropic metabolite β-aminoisobutyric acid. We manipulated the expression of genes associated with β-aminoisobutyric acid metabolism in I. scapularis, resulting in feeding impairment, diminished survival and reduced bacterial acquisition post haematophagy. Collectively, we discovered that metabolic reprogramming affects interspecies relationships and fitness in the clinically relevant tick I. scapularis.

Metagenome diversity illuminates the origins of pathogen effectors

Recent metagenome-assembled genome (MAG) analyses have profoundly impacted Rickettsiology systematics. The discovery of basal lineages (novel families Mitibacteraceae and Athabascaceae) with predicted extracellular lifestyles exposed an evolutionary timepoint for the transition to host dependency, which seemingly occurred independent of mitochondrial evolution. Notably, these basal rickettsiae carry the Rickettsiales vir homolog (rvh) type IV secretion system and purportedly use rvh to kill congener microbes rather than parasitize host cells as described for later-evolving rickettsial pathogens. MAG analysis also substantially increased diversity for the genus Rickettsia and delineated a sister lineage (the novel genus Tisiphia) that stands to inform on the emergence of human pathogens from protist and invertebrate endosymbionts. Herein, we probed Rickettsiales MAG and genomic diversity for the distribution of Rickettsia rvh effectors to ascertain their origins. A sparse distribution of most Rickettsia rvh effectors outside of Rickettsiaceae lineages illuminates unique rvh evolution from basal extracellular species and other rickettsial families. Remarkably, nearly every effector was found in multiple divergent forms with variable architectures, indicating profound roles for gene duplication and recombination in shaping effector repertoires in Rickettsia pathogens. Lateral gene transfer plays a prominent role in shaping the rvh effector landscape, as evinced by the discovery of many effectors on plasmids and conjugative transposons, as well as pervasive effector gene exchange between Rickettsia and Legionella species. Our study exemplifies how MAGs can yield insight into pathogen effector origins, particularly how effector architectures might become tailored to the discrete host cell functions of different eukaryotic hosts.IMPORTANCEWhile rickettsioses are deadly vector-borne human diseases, factors distinguishing Rickettsia pathogens from the innumerable bevy of environmental rickettsial endosymbionts remain lacking. Recent metagenome-assembled genome (MAG) studies revealed evolutionary timepoints for rickettsial transitions to host dependency. The rvh type IV secretion system was likely repurposed from congener killing in basal extracellular species to parasitizing host cells in later-evolving pathogens. Our analysis of MAG diversity for over two dozen rvh effectors unearthed their presence in some non-pathogens. However, most effectors were found in multiple divergent forms with variable architectures, indicating gene duplication and recombination-fashioned effector repertoires of Rickettsia pathogens. Lateral gene transfer substantially shaped pathogen effector arsenals, evinced by the discovery of effectors on plasmids and conjugative transposons, as well as pervasive effector gene exchanges between Rickettsia and Legionella species. Our study exemplifies how MAGs yield insight into pathogen effector origins and evolutionary processes tailoring effectors to eukaryotic host cell biology.

Single-cell genomics revealed Candidatus Grellia alia sp. nov. as an endosymbiont of Eutreptiella sp. (Euglenophyceae)

Though endosymbioses between protists and prokaryotes are widespread, certain host lineages have received disproportionate attention what may indicate either a predisposition to such interactions or limited studies on certain protist groups due to lack of cultures. The euglenids represent one such group in spite of microscopic observations showing intracellular bacteria in some strains. Here, we perform a comprehensive molecular analysis of a previously identified endosymbiont in the Eutreptiella sp. CCMP3347 using a single cell approach and bulk culture sequencing. The genome reconstruction of this endosymbiont allowed the description of a new endosymbiont Candidatus Grellia alia sp. nov. from the family Midichloriaceae. Comparative genomics revealed a remarkably complete conjugative type IV secretion system present in three copies on the plasmid sequences of the studied endosymbiont, a feature missing in the closely related Grellia incantans. This study addresses the challenge of limited host cultures with endosymbionts by showing that the genomes of endosymbionts reconstructed from single host cells have the completeness and contiguity that matches or exceeds those coming from bulk cultures. This paves the way for further studies of endosymbionts in euglenids and other protist groups. The research also provides the opportunity to study the diversity of endosymbionts in natural populations.

Orientia tsutsugamushi: comprehensive analysis of the mobilome of a highly fragmented and repetitive genome reveals the capacity for ongoing lateral gene transfer in an obligate intracellular bacterium

IMPORTANCE

Obligate intracellular bacteria, or those only capable of growth inside other living cells, have limited opportunities for horizontal gene transfer with other microbes due to their isolated replicative niche. The human pathogen Ot, an obligate intracellular bacterium causing scrub typhus, encodes an unusually high copy number of a ~40 gene mobile genetic element that typically facilitates genetic transfer across microbes. This proliferated element is heavily degraded in Ot and previously assumed to be inactive. Here, we conducted a detailed analysis of this element in eight Ot strains and discovered two strains with at least one intact copy. This implies that the element is still capable of moving across Ot populations and suggests that the genome of this bacterium may be even more dynamic than previously appreciated. Our work raises questions about intracellular microbial evolution and sounds an alarm for gene-based efforts focused on diagnosing and combatting scrub typhus.

Metagenomics of culture isolates and insect tissue illuminate the evolution of Wolbachia, Rickettsia and Bartonella symbionts in Ctenocephalides spp. fleas

While fleas are often perceived simply as a biting nuisance and a 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 (Ctenocephalides canis ) fleas, as well as Ctenocephalides 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., including Rickettsia felis (agent of flea-borne spotted fever) and Rickettsia asembonensis , a suspected pathogen. These Rickettsia spp. are members of a phylogenetic clade known as the ‘transitional group’, which includes both human pathogens and arthropod-specific endosymbionts. The relatively depauperate flea microbiome can also contain other endosymbionts, including a diverse range of Wolbachia strains. Here, we present circularized genome assemblies for two C. orientis -derived pathogens (Bartonella clarridgeiae and R. asembonensis ) from Malaysia, a novel Wolbachia strain (w Cori), and the C. orientis mitochondrion; all were 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 USA. 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 large-scale 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.

Metabolic disruption impacts tick fitness and microbial relationships

Arthropod-borne microbes rely on the metabolic state of a host to cycle between evolutionarily distant species. For instance, arthropod tolerance to infection may be due to redistribution of metabolic resources, often leading to microbial transmission to mammals. Conversely, metabolic alterations aids in pathogen elimination in humans, who do not ordinarily harbor arthropod-borne microbes. To ascertain the effect of metabolism on interspecies relationships, we engineered a system to evaluate glycolysis and oxidative phosphorylation in the tick Ixodes scapularis. Using a metabolic flux assay, we determined that the rickettsial bacterium Anaplasma phagocytophilum and the Lyme disease spirochete Borrelia burgdorferi, which are transstadially transmitted in nature, induced glycolysis in ticks. On the other hand, the endosymbiont Rickettsia buchneri, which is transovarially maintained, had a minimal effect on I. scapularis bioenergetics. Importantly, the metabolite β-aminoisobutyric acid (BAIBA) was elevated during A. phagocytophilum infection of tick cells following an unbiased metabolomics approach. Thus, we manipulated the expression of genes associated with the catabolism and anabolism of BAIBA in I. scapularis and detected impaired feeding on mammals, reduced bacterial acquisition, and decreased tick survival. Collectively, we reveal the importance of metabolism for tick-microbe relationships and unveil a valuable metabolite for I. scapularis fitness.

Orientia tsutsugamushi: analysis of the mobilome of a highly fragmented and repetitive genome reveals ongoing lateral gene transfer in an obligate intracellular bacterium

The rickettsial human pathogen Orientia tsutsugamushi (Ot) is an obligate intracellular Gram-negative bacterium with one of the most highly fragmented and repetitive genomes of any organism. Around 50% of its ~2.3 Mb genome is comprised of repetitive DNA that is derived from the highly proliferated Rickettsiales amplified genetic element (RAGE). RAGE is an integrative and conjugative element (ICE) that is present in a single Ot genome in up to 92 copies, most of which are partially or heavily degraded. In this report, we analysed RAGEs in eight fully sequenced Ot genomes and manually curated and reannotated all RAGE-associated genes, including those encoding DNA mobilisation proteins, P-type (vir) and F-type (tra) type IV secretion system (T4SS) components, Ankyrin repeat- and tetratricopeptide repeat-containing effectors, and other piggybacking cargo. Originally, the heavily degraded Ot RAGEs led to speculation that they are remnants of historical ICEs that are no longer active. Our analysis, however, identified two Ot genomes harbouring one or more intact RAGEs with complete F-T4SS genes essential for mediating ICE DNA transfer. As similar ICEs have been identified in unrelated rickettsial species, we assert that RAGEs play an ongoing role in lateral gene transfer within the Rickettsiales. Remarkably, we also identified in several Ot genomes remnants of prophages with no similarity to other rickettsial prophages. Together these findings indicate that, despite their obligate intracellular lifestyle and host range restricted to mites, rodents and humans, Ot genomes are highly dynamic and shaped through ongoing invasions by mobile genetic elements and viruses.

Orientia and Rickettsia: different flowers from the same garden

Recent discoveries of basal extracellular Rickettsiales have illuminated divergent evolutionary paths to host dependency in later-evolving lineages. Family Rickettsiaceae, primarily comprised of numerous protist- and invertebrate-associated species, also includes human pathogens from two genera, Orientia and Rickettsia. Once considered sister taxa, these bacteria form distinct lineages with newly appreciated lifestyles and morphological traits. Contrasting other rickettsial human pathogens in Family Anaplasmataceae, Orientia and Rickettsia species do not reside in host-derived vacuoles and lack glycolytic potential. With only a few described mechanisms, strategies for commandeering host glycolysis to support cytosolic growth remain to be discovered. While regulatory systems for this unique mode of intracellular parasitism are unclear, conjugative transposons unique to Orientia and Rickettsia species provide insights that are critical for determining how these obligate intracellular pathogens overtake eukaryotic cytosol.

Metagenome diversity illuminates origins of pathogen effectors

Recent metagenome assembled genome (MAG) analyses have profoundly impacted Rickettsiology systematics. Discovery of basal lineages (Mitibacteraceae and Athabascaceae) with predicted extracellular lifestyles reveals an evolutionary timepoint for the transition to host dependency, which occurred independent of mitochondrial evolution. Notably, these basal rickettsiae carry the Rickettsiales vir homolog (rvh) type IV secretion system (T4SS) and purportedly use rvh to kill congener microbes rather than parasitize host cells as described for derived rickettsial pathogens. MAG analysis also substantially increased diversity for genus Rickettsia and delineated a basal lineage (Tisiphia) that stands to inform on the rise of human pathogens from protist and invertebrate endosymbionts. Herein, we probed Rickettsiales MAG and genomic diversity for the distribution of Rickettsia rvh effectors to ascertain their origins. A sparse distribution of most Rickettsia rvh effectors outside of Rickettsiaceae lineages indicates unique rvh evolution from basal extracellular species and other rickettsial families. Remarkably, nearly every effector was found in multiple divergent forms with variable architectures, illuminating profound roles for gene duplication and recombination in shaping effector repertoires in Rickettsia pathogens. Lateral gene transfer plays a prominent role shaping the rvh effector landscape, as evinced by the discover of many effectors on plasmids and conjugative transposons, as well as pervasive effector gene exchange between Rickettsia and Legionella species. Our study exemplifies how MAGs can provide incredible insight on the origins of pathogen effectors and how their architectural modifications become tailored to eukaryotic host cell biology.

One to host them all: genomics of the diverse bacterial endosymbionts of the spider Oedothorax gibbosus

Bacterial endosymbionts of the groups Wolbachia, Cardinium and Rickettsiaceae are well known for their diverse effects on their arthropod hosts, ranging from mutualistic relationships to reproductive phenotypes. Here, we analysed a unique system in which the dwarf spider Oedothorax gibbosus is co-infected with up to five different endosymbionts affiliated with Wolbachia, 'Candidatus Tisiphia' (formerly Torix group Rickettsia), Cardinium and Rhabdochlamydia. Using short-read genome sequencing data, we show that the endosymbionts are heterogeneously distributed among O. gibbosus populations and are frequently found co-infecting spider individuals. To study this intricate host-endosymbiont system on a genome-resolved level, we used long-read sequencing to reconstruct closed genomes of the Wolbachia, 'Ca. Tisiphia' and Cardinium endosymbionts. We provide insights into the ecology and evolution of the endosymbionts and shed light on the interactions with their spider host. We detected high quantities of transposable elements in all endosymbiont genomes and provide evidence that ancestors of the Cardinium, 'Ca. Tisiphia' and Wolbachia endosymbionts have co-infected the same hosts in the past. Our findings contribute to broadening our knowledge about endosymbionts infecting one of the largest animal phyla on Earth and show the usefulness of transposable elements as an evolutionary 'contact-tracing' tool.

Coinfection of Two Rickettsia Species in a Single Tick Species Provides New Insight into Rickettsia-Rickettsia and Rickettsia-Vector Interactions

Rickettsiae are obligate intracellular bacteria that can cause life-threatening illnesses. There is an ongoing debate as to whether established infections by one Rickettsia species preclude the maintenance of the second species in ticks. Here, we identified two Rickettsia species in inoculum from Haemaphysalis montgomeryi ticks and subsequently obtained pure isolates of each species by plaque selection. The two isolates were classified as a transitional group and spotted fever group rickettsiae and named Rickettsia hoogstraalii str CS and Rickettsia rhipicephalii str EH, respectively. The coinfection of these two Rickettsia species was detected in 25.6% of individual field-collected H. montgomeryi. In cell culture infection models, R. hoogstraalii str CS overwhelmed R. rhipicephalii str EH with more obvious cytopathic effects, faster plaque formation, and increased cellular growth when cocultured, and R. hoogstraalii str CS seemed to polymerize actin tails differently from R. rhipicephalii str EH in vitro. This work provides a model to investigate the mechanisms of both Rickettsia-Rickettsia and Rickettsia-vector interactions. IMPORTANCE The rickettsiae are a group of obligate intracellular Gram-negative bacteria that include human pathogens causing an array of clinical symptoms and even death. There is an important question in the field, that is whether one infection can block the superinfection of other rickettsiae. This work demonstrated the coinfection of two Rickettsia species in individual ticks and further highlighted that testing the rickettsial competitive exclusion hypothesis will undoubtedly be a promising area as methods for bioengineering and pathogen biocontrol become amenable for rickettsiae.

Comparative Whole Genome Analysis of an Anaplasma phagocytophilum Strain Isolated from Norwegian Sheep

Anaplasma phagocytophilum is a Gram-negative obligate intracellular tick-borne alphaproteobacteria (family Anaplasmatacea, order Rickettsiales) with a worldwide distribution. In Norway, tick borne fever (TBF), caused by A. phagocytophilum, presents a major challenge in sheep farming. Despite the abundance of its tick vector, Ixodes ricinus, and A. phagocytophilum infections in wild and domestic animals, reports of infections in humans are low compared with cases in the U.S. Although A. phagocytophilum is genetically diverse and complex infections (co-infection and superinfection) in ruminants and other animals are common, the underlying genetic basis of intra-species interactions and host-specificity remains unexplored. Here, we performed whole genome comparative analysis of a newly cultured Norwegian A. phagocytophilum isolate from sheep (ApSheep_NorV1) with 27 other A. phagocytophilum genome sequences derived from human and animal infections worldwide. Although the compared strains are syntenic, there is remarkable genetic diversity between different genomic loci including the pfam01617 superfamily that encodes the major, neutralization-sensitive, surface antigen Msp2/p44. Blast comparisons between the msp2/p44 pseudogene repertoires from all the strains showed high divergence between U. S. and European strains and even between two Norwegian strains. Based on these comparisons, we concluded that in ruminants, complex infections can be attributed to infection with strains that differ in their msp2/p44 repertoires, which has important implications for pathogen evolution and vaccine development. We also present evidence for integration of rickettsial DNA into the genome of ISE6 tick cells.

Genomic diversity across the Rickettsia and 'Candidatus Megaira' genera and proposal of genus status for the Torix group

Members of the bacterial genus Rickettsia were originally identified as causative agents of vector-borne diseases in mammals. However, many Rickettsia species are arthropod symbionts and close relatives of 'Candidatus Megaira', which are symbiotic associates of microeukaryotes. Here, we clarify the evolutionary relationships between these organisms by assembling 26 genomes of Rickettsia species from understudied groups, including the Torix group, and two genomes of 'Ca. Megaira' from various insects and microeukaryotes. Our analyses of the new genomes, in comparison with previously described ones, indicate that the accessory genome diversity and broad host range of Torix Rickettsia are comparable to those of all other Rickettsia combined. Therefore, the Torix clade may play unrecognized roles in invertebrate biology and physiology. We argue this clade should be given its own genus status, for which we propose the name 'Candidatus Tisiphia'.

Cryptic Genes for Interbacterial Antagonism Distinguish Rickettsia Species Infecting Blacklegged Ticks From Other Rickettsia Pathogens

Background

The genus Rickettsia (Alphaproteobacteria: Rickettsiales) encompasses numerous obligate intracellular species with predominantly ciliate and arthropod hosts. Notable species are pathogens transmitted to mammals by blood-feeding arthropods. Mammalian pathogenicity evolved from basal, non-pathogenic host-associations; however, some non-pathogens are closely related to pathogens. One such species, Rickettsia buchneri, is prevalent in the blacklegged tick, Ixodes scapularis. While I. scapularis transmits several pathogens to humans, it does not transmit Rickettsia pathogens. We hypothesize that R. buchneri established a mutualism with I. scapularis, blocking tick superinfection with Rickettsia pathogens.

Methods

To improve estimates for assessing R. buchneri infection frequency in blacklegged tick populations, we used comparative genomics to identify an R. buchneri gene (REIS_1424) not present in other Rickettsia species present throughout the I. scapularis geographic range. Bioinformatic and phylogenomics approaches were employed to propose a function for the hypothetical protein (263 aa) encoded by REIS_1424.

Results

REIS_1424 has few analogs in other Rickettsiales genomes and greatest similarity to non-Proteobacteria proteins. This cohort of proteins varies greatly in size and domain composition, possessing characteristics of Recombination hotspot (Rhs) and contact dependent growth inhibition (CDI) toxins, with similarity limited to proximal C-termini (~145 aa). This domain was named CDI-like/Rhs-like C-terminal toxin (CRCT). As such proteins are often found as toxin-antidote (TA) modules, we interrogated REIS_1423 (151 aa) as a putative antidote. Indeed, REIS_1423 is similar to proteins encoded upstream of CRCT domain-containing proteins. Accordingly, we named these proteins CDI-like/Rhs-like C-terminal toxin antidotes (CRCA). R. buchneri expressed both REIS_1423 and REIS_1424 in tick cell culture, and PCR assays showed specificity for R. buchneri over other rickettsiae and utility for positive detection in three tick populations. Finally, phylogenomics analyses uncovered divergent CRCT/CRCA modules in varying states of conservation; however, only R. buchneri and related Tamurae/Ixodes Group rickettsiae carry complete TA modules.

Conclusion

We hypothesize that Rickettsia CRCT/CRCA modules circulate in the Rickettsia mobile gene pool, arming rickettsiae for battle over arthropod colonization. While its functional significance remains to be tested, R. buchneri CRCT/CRCA serves as a marker to positively identify infection and begin deciphering the role this endosymbiont plays in the biology of the blacklegged tick.

Genomic evolution and adaptation of arthropod-associated Rickettsia

Rickettsia species are endosymbionts hosted by arthropods and are known to cause mild to fatal diseases in humans. Here, we analyse the evolution and diversity of 34 Rickettsia species using a pangenomic meta-analysis (80 genomes/41 plasmids). Phylogenomic trees showed that Rickettsia spp. diverged into two Spotted Fever groups, a Typhus group, a Canadensis group and a Bellii group, and may have inherited their plasmids from an ancestral plasmid that persisted in some strains or may have been lost by others. The results suggested that the ancestors of Rickettsia spp. might have infected Acari and/or Insecta and probably diverged by persisting inside and/or switching hosts. Pangenomic analysis revealed that the Rickettsia genus evolved through a strong interplay between genome degradation/reduction and/or expansion leading to possible distinct adaptive trajectories. The genus mainly shared evolutionary relationships with α-proteobacteria, and also with γ/β/δ-proteobacteria, cytophagia, actinobacteria, cyanobacteria, chlamydiia and viruses, suggesting lateral exchanges of several critical genes. These evolutionary processes have probably been orchestrated by an abundance of mobile genetic elements, especially in the Spotted Fever and Bellii groups. In this study, we provided a global evolutionary genomic view of the intracellular Rickettsia that may help our understanding of their diversity, adaptation and fitness.

The Ixodes scapularis Symbiont Rickettsia buchneri Inhibits Growth of Pathogenic Rickettsiaceae in Tick Cells: Implications for Vector Competence

Ixodes scapularis is the primary vector of tick-borne pathogens in North America but notably does not transmit pathogenic Rickettsia species. This tick harbors the transovarially transmitted endosymbiont Rickettsia buchneri, which is widespread in I. scapularis populations, suggesting that it confers a selective advantage for tick survival such as providing essential nutrients. The R. buchneri genome includes genes with similarity to those involved in antibiotic synthesis. There are two gene clusters not found in other Rickettsiaceae, raising the possibility that these may be involved in excluding pathogenic bacteria from the tick. This study explored whether the R. buchneri antibiotic genes might exert antibiotic effects on pathogens associated with I. scapularis. Markedly reduced infectivity and replication of the tick-borne pathogens Anaplasma phagocytophilum, R. monacensis, and R. parkeri were observed in IRE11 tick cells hosting R. buchneri. Using a fluorescent plate reader assay to follow infection dynamics revealed that the presence of R. buchneri in tick cells, even at low infection rates, inhibited the growth of R. parkeri by 86-100% relative to R. buchneri-free cells. In contrast, presence of the low-pathogenic species R. amblyommatis or the endosymbiont R. peacockii only partially reduced the infection and replication of R. parkeri. Addition of host-cell free R. buchneri, cell lysate of R. buchneri-infected IRE11, or supernatant from R. buchneri-infected IRE11 cultures had no effect on R. parkeri infection and replication in IRE11, nor did these treatments show any antibiotic effect against non-obligate intracellular bacteria E. coli and S. aureus. However, lysate from R. buchneri-infected IRE11 challenged with R. parkeri showed some inhibitory effect on R. parkeri infection of treated IRE11, suggesting that challenge by pathogenic rickettsiae may induce the antibiotic effect of R. buchneri. This research suggests a potential role of the endosymbiont in preventing other rickettsiae from colonizing I. scapularis and/or being transmitted transovarially. The confirmation that the observed inhibition is linked to R. buchneri's antibiotic clusters requires further investigation but could have important implications for our understanding of rickettsial competition and vector competence of I. scapularis for rickettsiae.

Comparative hologenomics of two Ixodes scapularis tick populations in New Jersey

Tick-borne diseases, such as those transmitted by the blacklegged tick Ixodes scapularis, are a significant and growing public health problem in the US. There is mounting evidence that co-occurring non-pathogenic microbes can also impact tick-borne disease transmission. Shotgun metagenome sequencing enables sampling of the complete tick hologenome—the collective genomes of the tick and all of the microbial species contained therein, whether pathogenic, commensal or symbiotic. This approach simultaneously uncovers taxonomic composition and allows the detection of intraspecific genetic variation, making it a useful tool to compare spatial differences across tick populations. We evaluated this approach by comparing hologenome data from two tick samples (N = 6 ticks per location) collected at a relatively fine spatial scale, approximately 23 km apart, within a single US county. Several intriguing variants in the data between the two sites were detected, including polymorphisms in both in the tick’s own mitochondrial DNA and that of a rickettsial endosymbiont. The two samples were broadly similar in terms of the microbial species present, including multiple known tick-borne pathogens (Borrelia burgdorferi, Babesia microti, and Anaplasma phagocytophilum), filarial nematodes, and Wolbachia and Babesia species. We assembled the complete genome of the rickettsial endosymbiont (most likely Rickettsia buchneri) from both populations. Our results provide further evidence for the use of shotgun metagenome sequencing as a tool to compare tick hologenomes and differentiate tick populations across localized spatial scales.

Tick Cell Culture Analysis of Growth Dynamics and Cellular Tropism of Rickettsia buchneri, an Endosymbiont of the Blacklegged Tick, Ixodes scapularis

The blacklegged tick, Ixodes scapularis, a species of significant importance to human and animal health, harbors an endosymbiont Rickettsia buchneri sensu stricto. The symbiont is largely restricted to the ovaries, but all life stages can harbor various quantities or lack R. buchneri entirely. The endosymbiont is cultivable in cell lines isolated from embryos of Ixodes ticks. Rickettsia buchneri most readily grows and is maintained in the cell line IRE11 from the European tick, Ixodes ricinus. The line was characterized by light and electron microscopy and used to analyze the growth dynamics of wildtype and GFPuv-expressing R. buchneri. qPCR indicated that the genome copy doubling time in IRE11 was >7 days. Measurements of fluorescence using a plate reader indicated that the amount of green fluorescent protein doubled every 11 days. Two 23S rRNA probes were tested via RNA FISH on rickettsiae grown in vitro and adapted to evaluate the tissue tropism of R. buchneri in field-collected female I. scapularis. We observed strong positive signals of R. buchneri in the ovaries and surrounding the nucleus of the developing oocytes. Tissue tropism in I. scapularis and in vitro growth dynamics strengthen the contemporary understanding of R. buchneri as a transovarially transmitted, non-pathogenic endosymbiont.

Unpacking the intricacies of Rickettsia-vector interactions

Although Rickettsia species are molecularly detected among a wide range of arthropods, vector competence becomes an imperative aspect of understanding the ecoepidemiology of these vector-borne diseases. The synergy between vector homeostasis and rickettsial invasion, replication, and release initiated within hours (insects) and days (ticks) permits successful transmission of rickettsiae. Uncovering the molecular interplay between rickettsiae and their vectors necessitates examining the multifaceted nature of rickettsial virulence and vector infection tolerance. Here, we highlight the biological differences between tick- and insect-borne rickettsiae and the factors facilitating the incidence of rickettsioses. Untangling the complex relationship between rickettsial genetics, vector biology, and microbial interactions is crucial in understanding the intricate association between rickettsiae and their vectors.

Lipid A Structural Divergence in Rickettsia Pathogens

Spikes in rickettsioses occur as deforestation, urbanization, and homelessness increase human exposure to blood-feeding arthropods. Still, effective Rickettsia vaccines remain elusive.

Species of Rickettsia (Alphaproteobacteria: Rickettsiales) are obligate intracellular parasites of a wide range of eukaryotes, with recognized arthropod-borne human pathogens belonging to the transitional group (TRG), typhus group (TG), and spotted fever group (SFG) rickettsiae. Growing in the host cytosol, rickettsiae pilfer numerous metabolites to make a typical Gram-negative bacterial cell envelope. The O-antigen of rickettsial lipopolysaccharide (LPS) is immunogenic and has been shown to tether the S-layer to the rickettsial surface; however, little is known about the structure and immunogenicity of the Rickettsia lipid A moiety. The structure of lipid A, the membrane anchor of LPS, affects the ability of this molecule to interact with components of the host innate immune system, specifically the MD-2/TLR4 receptor complex. To dissect the host responses that can occur during Rickettsia in vitro and in vivo infection, structural analysis of Rickettsia lipid A is needed. Lipid A was extracted from four Rickettsia species and structurally analyzed. R. akari (TRG), R. typhi (TG), and R. montanensis (SFG) produced a similar structure, whereas R. rickettsii (SFG) altered the length of a secondary acyl group. While all structures have longer acyl chains than known highly inflammatory hexa-acylated lipid A structures, the R. rickettsii modification should differentially alter interactions with the hydrophobic internal pocket in MD2. The significance of these characteristics toward inflammatory potential as well as membrane dynamics between arthropod and vertebrate cellular environments warrants further investigation. Our work adds lipid A to the secretome and O-antigen as variable factors possibly correlating with phenotypically diverse rickettsioses.

IMPORTANCE Spikes in rickettsioses occur as deforestation, urbanization, and homelessness increase human exposure to blood-feeding arthropods. Still, effective Rickettsia vaccines remain elusive. Recent studies have determined that Rickettsia lipopolysaccharide anchors the protective S-layer to the bacterial surface and elicits bactericidal antibodies. Furthermore, growing immunological evidence suggests vertebrate sensors (MD-2/TLR4 and noncanonical inflammasome) typically triggered by the lipid A portion of lipopolysaccharide are activated during Rickettsia infection. However, the immunopotency of Rickettsia lipid A is unknown due to poor appreciation for its structure. We determined lipid A structures for four distinct rickettsiae, revealing longer acyl chains relative to highly inflammatory bacterial lipid A. Surprisingly, lipid A of the Rocky Mountain spotted fever agent deviates in structure from other rickettsiae. Thus, lipid A divergence may contribute to variable disease phenotypes, sounding an alarm for determining its immunopotency and possible utility (i.e., as an adjuvant or anti-inflammatory) for development of more prudent rickettsiacidal therapies.

Cells within cells: Rickettsiales and the obligate intracellular bacterial lifestyle

The Rickettsiales are a group of obligate intracellular vector-borne Gram-negative bacteria that include many organisms of clinical and agricultural importance, including Anaplasma spp., Ehrlichia chaffeensis, Wolbachia, Rickettsia spp. and Orientia tsutsugamushi. This Review provides an overview of the current state of knowledge of the biology of these bacteria and their interactions with host cells, with a focus on pathogenic species or those that are otherwise important for human health. This includes a description of rickettsial genomics, bacterial cell biology, the intracellular lifestyles of Rickettsiales and the mechanisms by which they induce and evade the innate immune response.

Growth Dynamics and Antibiotic Elimination of Symbiotic Rickettsia buchneri in the Tick Ixodes scapularis (Acari: Ixodidae)

Rickettsia buchneri is the principal symbiotic bacterium of the medically significant tick Ixodes scapularis This species has been detected primarily in the ovaries of adult female ticks and is vertically transmitted, but its tissue tropism in other life stages and function with regard to tick physiology is unknown. In order to determine the function of R. buchneri, it may be necessary to produce ticks free from this symbiont. We quantified the growth dynamics of R. buchneri naturally occurring in I. scapularis ticks throughout their life cycle and compared it with bacterial growth in ticks in which symbiont numbers were experimentally reduced or eliminated. To eliminate the bacteria, we exposed ticks to antibiotics through injection and artificial membrane feeding. Both injection and membrane feeding of the antibiotic ciprofloxacin were effective at eliminating R. buchneri from most offspring of exposed females. Because of its effectiveness and ease of use, we have determined that injection of ciprofloxacin into engorged female ticks is an efficient means of clearing R. buchneri from the majority of progeny.IMPORTANCE This paper describes the growth of symbiotic Rickettsia buchneri within Ixodes scapularis through the life cycle of the tick and provides methods to eliminate R. buchneri from I. scapularis ticks.

Evolution of Wolbachia mutualism and reproductive parasitism: insight from two novel strains that co-infect cat fleas

Wolbachiae are obligate intracellular bacteria that infect arthropods and certain nematodes. Usually maternally inherited, they may provision nutrients to (mutualism) or alter sexual biology of (reproductive parasitism) their invertebrate hosts. We report the assembly of closed genomes for two novel wolbachiae, wCfeT and wCfeJ, found co-infecting cat fleas (Ctenocephalides felis) of the Elward Laboratory colony (Soquel, CA, USA). wCfeT is basal to nearly all described Wolbachia supergroups, while wCfeJ is related to supergroups C, D and F. Both genomes contain laterally transferred genes that inform on the evolution of Wolbachia host associations. wCfeT carries the Biotin synthesis Operon of Obligate intracellular Microbes (BOOM); our analyses reveal five independent acquisitions of BOOM across the Wolbachia tree, indicating parallel evolution towards mutualism. Alternately, wCfeJ harbors a toxin-antidote operon analogous to the wPip cinAB operon recently characterized as an inducer of cytoplasmic incompatibility (CI) in flies. wCfeJ cinB and three adjacent genes are collectively similar to large modular toxins encoded in CI-like operons of certain Wolbachia strains and Rickettsia species, signifying that CI toxins streamline by fission of large modular toxins. Remarkably, the C. felis genome itself contains two CI-like antidote genes, divergent from wCfeJ cinA, revealing episodic reproductive parasitism in cat fleas and evidencing mobility of CI loci independent of WO-phage. Additional screening revealed predominant co-infection (wCfeT/wCfeJ) amongst C. felis colonies, though fleas in wild populations mostly harbor wCfeT alone. Collectively, genomes of wCfeT, wCfeJ, and their cat flea host supply instances of lateral gene transfers that could drive transitions between parasitism and mutualism.

The Mosaic Architecture of NRPS-PKS in the Arbuscular Mycorrhizal Fungus Gigaspora margarita Shows a Domain With Bacterial Signature

As obligate biotrophic symbionts, arbuscular mycorrhizal fungi (AMF) live in association with most land plants. Among them, Gigaspora margarita has been deeply investigated because of its peculiar features, i.e., the presence of an intracellular microbiota with endobacteria and viruses. The genome sequencing of this fungus revealed the presence of some hybrid non-ribosomal peptide synthases-polyketide synthases (NRPS-PKS) that have been rarely identified in AMF. The aim of this study is to describe the architecture of these NRPS-PKS sequences and to understand whether they are present in other fungal taxa related to G. margarita. A phylogenetic analysis shows that the ketoacyl synthase (KS) domain of one G. margarita NRPS-PKS clusters with prokaryotic sequences. Since horizontal gene transfer (HGT) has often been advocated as a relevant evolutionary mechanism for the spread of secondary metabolite genes, we hypothesized that a similar event could have interested the KS domain of the PKS module. The bacterial endosymbiont of G. margarita, Candidatus Glomeribacter gigasporarum (CaGg), was the first candidate as a donor, since it possesses a large biosynthetic cluster involving an NRPS-PKS. However, bioinformatics analyses do not confirm the hypothesis of a direct HGT from the endobacterium to the fungal host: indeed, endobacterial and fungal sequences show a different evolution and potentially different donors. Lastly, by amplifying a NRPS-PKS conserved fragment and mining the sequenced AMF genomes, we demonstrate that, irrespective of the presence of CaGg, G. margarita, and some other related Gigasporaceae possess such a sequence.

Oxford Nanopore MinION Sequencing Enables Rapid Whole Genome Assembly of Rickettsia typhi in a Resource-Limited Setting

The infrastructure challenges and costs of next-generation sequencing have been largely overcome, for many sequencing applications, by Oxford Nanopore Technologies' portable MinION sequencer. However, the question remains open whether MinION-based bacterial whole genome sequencing is by itself sufficient for the accurate assessment of phylogenetic and epidemiological relationships between isolates and whether such tasks can be undertaken in resource-limited settings. To investigate this question, we sequenced the genome of an isolate of Rickettsia typhi, an important and neglected cause of fever across much of the tropics and subtropics, for which only three genomic sequences previously existed. We prepared and sequenced libraries on a MinION in Vientiane, Lao PDR, using v9.5 chemistry, and in parallel, we sequenced the same isolate on the Illumina platform in a genomics laboratory in the United Kingdom. The MinION sequence reads yielded a single contiguous assembly, in which the addition of Illumina data revealed 226 base-substitution and 5,856 indel errors. The combined assembly represents the first complete genome sequence of a human R. typhi isolate collected in the last 50 years and differed from the genomes of existing strains collected over a 90-year time period at very few sites, with no rearrangements. Filtering based on the known error profile of MinION data improved the accuracy of the nanopore-only assembly. However, the frequency of false-positive errors remained greater than true sequence divergence from recorded sequences. Although nanopore-only sequencing cannot yet recover phylogenetic signals in R. typhi, such an approach may be applicable for more diverse organisms.

A Metabolic Dependency for Host Isoprenoids in the Obligate Intracellular Pathogen Rickettsia parkeri Underlies a Sensitivity to the Statin Class of Host-Targeted Therapeutics

Obligate intracellular pathogens, which include viruses as well as certain bacteria and eukaryotes, are a subset of infectious microbes that are metabolically dependent on and unable to grow outside an infected host cell because they have lost or lack essential biosynthetic pathways. In this study, we describe a metabolic dependency of the bacterial pathogen Rickettsia parkeri on host isoprenoid molecules that are used in the biosynthesis of downstream products, including cholesterol, steroid hormones, and heme. Bacteria make products from isoprenoids, such as an essential lipid carrier for making the bacterial cell wall. We show that bacterial metabolic dependency can represent a potential Achilles’ heel and that inhibiting host isoprenoid biosynthesis with the FDA-approved statin class of drugs inhibits bacterial growth by interfering with the integrity of the cell wall. This work supports the potential to treat infections by obligate intracellular pathogens through inhibition of host biosynthetic pathways that are susceptible to parasitism.

Gram-negative bacteria in the order Rickettsiales have an obligate intracellular growth requirement, and some species cause human diseases such as typhus and spotted fever. The bacteria have evolved a dependence on essential nutrients and metabolites from the host cell as a consequence of extensive genome reduction. However, it remains largely unknown which nutrients they acquire and whether their metabolic dependency can be exploited therapeutically. Here, we describe a genetic rewiring of bacterial isoprenoid biosynthetic pathways in the Rickettsiales that has resulted from reductive genome evolution. Furthermore, we investigated whether the spotted fever group Rickettsia species Rickettsia parkeri scavenges isoprenoid precursors directly from the host. Using targeted mass spectrometry, we found that infection caused decreases in host isoprenoid products and concomitant increases in bacterial isoprenoid metabolites. Additionally, we report that treatment of infected cells with statins, which inhibit host isoprenoid synthesis, prohibited bacterial growth. We show that growth inhibition correlates with changes in bacterial size and shape that mimic those caused by antibiotics that inhibit peptidoglycan biosynthesis, suggesting that statins lead to an inhibition of cell wall synthesis. Altogether, our results describe a potential Achilles’ heel of obligate intracellular pathogens that can potentially be exploited with host-targeted therapeutics that interfere with metabolic pathways required for bacterial growth.

IMPORTANCE Obligate intracellular pathogens, which include viruses as well as certain bacteria and eukaryotes, are a subset of infectious microbes that are metabolically dependent on and unable to grow outside an infected host cell because they have lost or lack essential biosynthetic pathways. In this study, we describe a metabolic dependency of the bacterial pathogen Rickettsia parkeri on host isoprenoid molecules that are used in the biosynthesis of downstream products, including cholesterol, steroid hormones, and heme. Bacteria make products from isoprenoids, such as an essential lipid carrier for making the bacterial cell wall. We show that bacterial metabolic dependency can represent a potential Achilles’ heel and that inhibiting host isoprenoid biosynthesis with the FDA-approved statin class of drugs inhibits bacterial growth by interfering with the integrity of the cell wall. This work supports the potential to treat infections by obligate intracellular pathogens through inhibition of host biosynthetic pathways that are susceptible to parasitism.

Rickettsia buchneri, symbiont of the deer tick Ixodes scapularis, can colonise the salivary glands of its host

Vertically-transmitted bacterial symbionts are widespread in ticks and have manifold impacts on the epidemiology of tick-borne diseases. For instance, they may provide essential nutrients to ticks, affect vector competence, induce immune responses in vertebrate hosts, or even evolve to become vertebrate pathogens. The deer or blacklegged tick Ixodes scapularis harbours the symbiont Rickettsia buchneri in its ovarian tissues. Here we show by molecular, proteomic and imaging methods that R. buchneri is also capable of colonising the salivary glands of wild I. scapularis. This finding has important implications for the diagnosis of rickettsial infections and for pathogen-symbiont interactions in this notorious vector of Lyme borreliosis.