Adaptive evolution and recombination of Rickettsia antigens

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.

Molecular Characterization by Multilocus Sequence Typing and Diversity Analysis of Rickettsia asembonensis in Peru

Despite several reports worldwide documenting the presence of Rickettsia asembonensis in samples derived from ectoparasites, animals and more recently humans, genomic information of these specimens remains scarce, and when available, is usually limited to small genomic fragments of limited value. We generated complete sequences for two conserved (17-kDa antigen gene and gltA) and three variable (sca4, ompB and ompA) genes in five R. asembonensis DNA samples detected in cat and dog fleas in Peru. Complete gene sequences were used to conduct multi-locus sequence typing and phylogenetic analyses to assess diversity and infer relationships among strains and other reference sequences. The 17-kDa antigen gene was highly conserved across Rickettsia species. Of the variable genes ompB was the most variable, but this diversity was not captured through phylogenetics alone even when efforts were made to maximize potential diversity in terms of flea species, animal host and location. Through a combination of de novo and reference-based genome assembly we identified a 75 bp insertion in ompA that encodes a 25 aa repetitive motif found in other Rickettsia species, but not present in the original prototype strain from Kenya. R. asembonensis has only recently been shown to be a bona-fide human pathogen. As such, and compounded by a lack of available genomic information, it remains understudied. Our work directly addresses the lack of genomic information available worldwide for the study of these novel Rickettsia species and specifically contributes to our understanding of the diversity and molecular epidemiology of R. asembonensis in Peru.

Rickettsia association with two Macrolophus (Heteroptera: Miridae) species: A comparative study of phylogenies and within-host localization patterns

Many arthropods host bacterial symbionts, some of which are known to influence host nutrition and diet breadth. Omnivorous bugs of the genus Macrolophus (Heteroptera: Miridae) are mainly predatory, but may also feed on plants. The species M. pygmaeus and M. melanotoma (=M. caliginosus) are key natural enemies of various economically important agricultural pests, and are known to harbor two Rickettsia species, R. bellii and R. limoniae. To test for possible involvement of symbiotic bacteria in the nutritional ecology of these biocontrol agents, the abundance, phylogeny, and distribution patterns of the two Rickettsia species in M. pygmaeus and M. melanotoma were studied. Both of the Rickettsia species were found in 100 and 84% of all tested individuals of M. pygmaeus and M. melanotoma, respectively. Phylogenetic analysis showed that a co-evolutionary process between Macrolophus species and their Rickettsia is infrequent. Localization of R. bellii and R. limoniae has been detected in both female and male of M. pygmaeus and M. melanotoma. FISH analysis of female gonads revealed the presence of both Rickettsia species in the germarium of both bug species. Each of the two Rickettsia species displayed a unique distribution pattern along the digestive system of the bugs, mostly occupying separate epithelial cells, unknown caeca-like organs, the Malpighian tubules and the salivary glands. This pattern differed between the two Macrolophus species: in M. pygmaeus, R. limoniae was distributed more broadly along the host digestive system and R. bellii was located primarily in the foregut and midgut. In contrast, in M. melanotoma, R. bellii was more broadly distributed along the digestive system than the clustered R. limoniae. Taken together, these results suggest that Rickettsia may have a role in the nutritional ecology of their plant-and prey-consuming hosts.

Occurrence and Genetic Diversity of Bartonella spp. (Rhizobiales: Bartonellaceae) and Rickettsia spp. (Rickettsiales: Rickettsiaceae) in Cat Fleas (Siphonaptera: Pulicidae) From Chile

The aim of this study was to perform a molecular survey and determine the genetic diversity of Bartonella spp. (Rhizobiales: Bartonellaceae) and Rickettsia spp. (Rickettsiales: Rickettsiaceae) in cat fleas (Siphonaptera: Pulicidae) from Southern Chile. Fleas (n = 251) were collected from 150 cats in Valdivia city and identified using morphological keys. Fleas belonging to the same cat were pooled (two to seven fleas per pool). DNA was purified from individual (n = 92) and pooled (n = 58) fleas and submitted to conventional polymerase chain reaction assays targeting Bartonella spp. (gltA) and Rickettsia spp. (ompA). Selected positive samples were sequenced for Bartonella gltA (n = 19), Rickettsia ompA (n = 14), and Rickettsia gltA (n = 11) for speciation, phylogenetic, and diversity analyses. All fleas (n = 251) were identified as Ctenocephalides felis felis (Bouché) (Siphonaptera: Pulicidae). Bartonella and Rickettsia occurrences in C. felis felis were 39.3% (59/150) and 76.6% (115/150), respectively. From sequenced Bartonella spp., 47.3% (9/19) were identified as Bartonella clarridgeiae, 42.1% (8/19) as Bartonella henselae, 5.3% (1/19) as Bartonella koehlerae, and 5.3% (1/19) as Bartonella sp. Rickettsia felis was the only Rickettsiaceae species identified in both ompA (14/14) and gltA (11/11) products. B. henselae and B. clarridgeiae presented five genotypes. R. felis ompA sequences presented seven genotypes. On the other hand, R. felis gltA sequences showed only one genotype. Bartonella spp. and R. felis are described for the first time in C. felis felis fleas from Southern Chile, highlighting the importance of these vectors as a source of zoonotic agents.

Molecular Detection of Vector-Borne Pathogens in Rural Dogs and Associated Ctenocephalides felis Fleas (Siphonaptera: Pulicidae) in Easter Island (Chile)

The presence of vector-borne pathogens of veterinary and public health interest have received little attention in Chile. In Easter Island, in particular, a Chilean territory in the southeastern Pacific Ocean, no information is available. To fill this gap, 153 rural dogs were inspected for ectoparasites during a sterilization campaign carried out in 2016. Fleas were observed in 46% of the dogs, and Ctenocephalides felis (Bouché, 1835) was the only species present. Morphological identification of fleas was genetically confirmed using conventional polymerase chain reaction targeting the cox2 gene. No tick was observed in any dog. The presence of DNA of Rickettsia sp. (gltA and ompA fragment genes), Anaplasmataceae (16S rRNA), and Bartonella sp. (16S-23S ribosomal RNA intergenic spacer) was investigated in blood samples of 70 of the dogs and in 126 fleas analyzed in 68 pools that included 1-5 fleas. Rickettsial DNA was detected in 97% (n = 66) of the flea pools. Of these, 57 showed between 99 and 100% identity for both genes with published sequences of Candidatus Rickettsia asemboensis (CRa), six with Rickettsia felis, and one with Candidatus Rickettsia senegalensis. For two pools, gltA amplicons were identical to CRa but ompB amplicions showed 99-100% identity with R. felis. Anaplasmataceae DNA was detected in 16% (n = 11) pools. Sequenced amplicons showed highest identity with the endosymbiont Wolbachia pipientis. Bartonella DNA, showing 99% identity to Bartonella clarridgeiae, was detected in one pool (1.4%). No positive reaction was observed for any dog. This is the first detection of members of the 'R. felis-like' group other than R. felis in Chile.

Uniparental Inheritance Promotes Adaptive Evolution in Cytoplasmic Genomes

Eukaryotes carry numerous asexual cytoplasmic genomes (mitochondria and plastids). Lacking recombination, asexual genomes should theoretically suffer from impaired adaptive evolution. Yet, empirical evidence indicates that cytoplasmic genomes experience higher levels of adaptive evolution than predicted by theory. In this study, we use a computational model to show that the unique biology of cytoplasmic genomes-specifically their organization into host cells and their uniparental (maternal) inheritance-enable them to undergo effective adaptive evolution. Uniparental inheritance of cytoplasmic genomes decreases competition between different beneficial substitutions (clonal interference), promoting the accumulation of beneficial substitutions. Uniparental inheritance also facilitates selection against deleterious cytoplasmic substitutions, slowing Muller's ratchet. In addition, uniparental inheritance generally reduces genetic hitchhiking of deleterious substitutions during selective sweeps. Overall, uniparental inheritance promotes adaptive evolution by increasing the level of beneficial substitutions relative to deleterious substitutions. When we assume that cytoplasmic genome inheritance is biparental, decreasing the number of genomes transmitted during gametogenesis (bottleneck) aids adaptive evolution. Nevertheless, adaptive evolution is always more efficient when inheritance is uniparental. Our findings explain empirical observations that cytoplasmic genomes-despite their asexual mode of reproduction-can readily undergo adaptive evolution.

Comparative genome sequencing of Rickettsia rickettsii strains that differ in virulence

Rickettsia rickettsii is an obligate intracellular pathogen that is the causative agent of Rocky Mountain spotted fever. Strains of R. rickettsii differ dramatically in virulence. In a guinea pig model of infection, the severity of disease as assessed by fever response varies from the most virulent, Sheila Smith, to Iowa, which causes no fever. To identify potential determinants of virulence in R. rickettsii, the genomes of two additional strains were sequenced for comparison to known sequences (comparative genome sequencing [CGS]). R. rickettsii Morgan and R strains were compared to the avirulent R. rickettsii Iowa and virulent R. rickettsii Sheila Smith strains. The Montana strains Sheila Smith and R were found to be highly similar while the eastern strains Iowa and Morgan were most similar to each other. A major surface antigen, rickettsial outer membrane protein A (rOmpA), is severely truncated in the Iowa strain. The region of ompA containing 13 tandem repeats was sequenced, revealing only seven shared SNPs (four nonsynonymous) for R and Morgan strains compared to Sheila Smith, with an additional 17 SNPs identified in Morgan. Another major surface antigen and autotransporter, rOmpB, exhibits a defect in processing in the Iowa strain such that the beta fragment is not cleaved. Sequence analysis of ompB reveals identical sequences between Iowa and Morgan strains and between the R and Sheila Smith strains. The number of SNPs and insertions/deletions between sequences of the two Montana strains and the two eastern strains is low, thus narrowing the field of possible virulence factors.

Secretome of obligate intracellular Rickettsia

The genus Rickettsia (Alphaproteobacteria, Rickettsiales, Rickettsiaceae) is comprised of obligate intracellular parasites, with virulent species of interest both as causes of emerging infectious diseases and for their potential deployment as bioterrorism agents. Currently, there are no effective commercially available vaccines, with treatment limited primarily to tetracycline antibiotics, although others (e.g. josamycin, ciprofloxacin, chloramphenicol, and azithromycin) are also effective. Much of the recent research geared toward understanding mechanisms underlying rickettsial pathogenicity has centered on characterization of secreted proteins that directly engage eukaryotic cells. Herein, we review all aspects of the Rickettsia secretome, including six secretion systems, 19 characterized secretory proteins, and potential moonlighting proteins identified on surfaces of multiple Rickettsia species. Employing bioinformatics and phylogenomics, we present novel structural and functional insight on each secretion system. Unexpectedly, our investigation revealed that the majority of characterized secretory proteins have not been assigned to their cognate secretion pathways. Furthermore, for most secretion pathways, the requisite signal sequences mediating translocation are poorly understood. As a blueprint for all known routes of protein translocation into host cells, this resource will assist research aimed at uniting characterized secreted proteins with their apposite secretion pathways. Furthermore, our work will help in the identification of novel secreted proteins involved in rickettsial 'life on the inside'.

The role of selection in shaping diversity of natural M. tuberculosis populations

Mycobacterium tuberculosis (M.tb), the cause of tuberculosis (TB), is estimated to infect a new host every second. While analyses of genetic data from natural populations of M.tb have emphasized the role of genetic drift in shaping patterns of diversity, the influence of natural selection on this successful pathogen is less well understood. We investigated the effects of natural selection on patterns of diversity in 63 globally extant genomes of M.tb and related pathogenic mycobacteria. We found evidence of strong purifying selection, with an estimated genome-wide selection coefficient equal to -9.5 × 10(-4) (95% CI -1.1 × 10(-3) to -6.8 × 10(-4)); this is several orders of magnitude higher than recent estimates for eukaryotic and prokaryotic organisms. We also identified different patterns of variation across categories of gene function. Genes involved in transport and metabolism of inorganic ions exhibited very low levels of non-synonymous polymorphism, equivalent to categories under strong purifying selection (essential and translation-associated genes). The highest levels of non-synonymous variation were seen in a group of transporter genes, likely due to either diversifying selection or local selective sweeps. In addition to selection, we identified other important influences on M.tb genetic diversity, such as a 25-fold expansion of global M.tb populations coincident with explosive growth in human populations (estimated timing 1684 C.E., 95% CI 1620-1713 C.E.). These results emphasize the parallel demographic histories of this obligate pathogen and its human host, and suggest that the dominant effect of selection on M.tb is removal of novel variants, with exceptions in an interesting group of genes involved in transportation and defense. We speculate that the hostile environment within a host imposes strict demands on M.tb physiology, and thus a substantial fitness cost for most new mutations. In this respect, obligate bacterial pathogens may differ from other host-associated microbes such as symbionts.

Surface proteome analysis and characterization of surface cell antigen (Sca) or autotransporter family of Rickettsia typhi

Surface proteins of the obligate intracellular bacterium Rickettsia typhi, the agent of murine or endemic typhus fever, comprise an important interface for host-pathogen interactions including adherence, invasion and survival in the host cytoplasm. In this report, we present analyses of the surface exposed proteins of R. typhi based on a suite of predictive algorithms complemented by experimental surface-labeling with thiol-cleavable sulfo-NHS-SS-biotin and identification of labeled peptides by LC MS/MS. Further, we focus on proteins belonging to the surface cell antigen (Sca) autotransporter (AT) family which are known to be involved in rickettsial infection of mammalian cells. Each species of Rickettsia has a different complement of sca genes in various states; R. typhi, has genes sca1 thru sca5. In silico analyses indicate divergence of the Sca paralogs across the four Rickettsia groups and concur with previous evidence of positive selection. Transcripts for each sca were detected during infection of L929 cells and four of the five Sca proteins were detected in the surface proteome analysis. We observed that each R. typhi Sca protein is expressed during in vitro infections and selected Sca proteins were expressed during in vivo infections. Using biotin-affinity pull down assays, negative staining electron microscopy, and flow cytometry, we demonstrate that the Sca proteins in R. typhi are localized to the surface of the bacteria. All Scas were detected during infection of L929 cells by immunogold electron microscopy. Immunofluorescence assays demonstrate that Scas 1–3 and 5 are expressed in the spleens of infected Sprague-Dawley rats and Scas 3, 4 and 5 are expressed in cat fleas (Ctenocephalides felis). Sca proteins may be crucial in the recognition and invasion of different host cell types. In short, continuous expression of all Scas may ensure that rickettsiae are primed i) to infect mammalian cells should the flea bite a host, ii) to remain infectious when extracellular and iii) to infect the flea midgut when ingested with a blood meal. Each Sca protein may be important for survival of R. typhi and the lack of host restricted expression may indicate a strategy of preparedness for infection of a new host.

Rickettsia typhi, a member of the typhus group (TG) rickettsia, is the agent of murine or endemic typhus fever – a disease exhibiting mild to severe flu-like symptoms resulting in significant morbidity. It is maintained in a flearodent transmission cycle in urban and suburban environments. The obligate intracellular lifestyle of rickettsiae makes genetic manipulation difficult and impedes progress towards identification of virulence factors. All five Scas were detected on the surface of R.. typhi using a combination of a biotin-labeled affinity assay, negative stain electron microscopy and flow cytometry. Sca proteins are members of the autotransporter (AT) family or type V secretion system (TVSS). We employed detailed bioinformatic analyses and evaluated their transcript abundance in an in vitro infection model where sca transcripts are detected at varying levels over the course of a 5 day in vitro infection. We also observe expression of selected Sca proteins during infection of fleas and rats. Our study provides a proteomic analysis of the bacterial surface and an initial characterization of the Sca family as it exists in R. typhi.

Comparative analysis of Wolbachia surface protein in D. melanoagster, A. tabida and B. malayi

Wolbachia surface protein (WSP) is an eight beta-barrel transmembrane structure which participates in host immune response, cell proliferation, pathogenicity and controlled cell death program. The protein has four extracellular loops containing hyper variable regions separated by conserved regions. The WSP structure is homologous to Neisseria surface protein (Nsp A) which has about 34% similarity including antigenic variation and hydrophilicity. Recombination has a large impact on diversity of this protein including positive selection which is major constraint on protein evolution. The molecular mechanism through which Wolbachia induces various reproductive anomalies is unclear; a key feature observed for such anomalies might be because of Wolbachia undergoing extensive recombination. In Wolbachia, increased recombination is observed in ankyrin proteins, surface proteins and in some hypothetical proteins. Genetic divergence is extensive in the WSP gene, WSP is known to be a chimeric protein involved in host-symbiont interactions. Here we predicted the structural and functional variations in WSP sequences of Wolbachia present in D. melanogaster, A. tabida and in B. malayi.

Accelerated microevolution in an outer membrane protein (OMP) of the intracellular bacteria Wolbachia

Background

Outer membrane proteins (OMPs) of Gram-negative bacteria are key players in the biology of bacterial-host interactions. However, while considerable attention has been given to OMPs of vertebrate pathogens, relatively little is known about the role of these proteins in bacteria that primarily infect invertebrates. One such OMP is found in the intracellular bacteria Wolbachia, which are widespread symbionts of arthropods and filarial nematodes. Recent experimental studies have shown that the Wolbachia surface protein (WSP) can trigger host immune responses and control cell death programming in humans, suggesting a key role of WSP for establishment and persistence of the symbiosis in arthropods.

Results

Here we performed an analysis of 515 unique alleles found in 831 Wolbachia isolates, to investigate WSP structure, microevolution and population genetics. WSP shows an eight-strand transmembrane β-barrel structure with four extracellular loops containing hypervariable regions (HVRs). A clustering approach based upon patterns of HVR haplotype diversity was used to group similar WSP sequences and to estimate the relative contribution of mutation and recombination during early stages of protein divergence. Results indicate that although point mutations generate most of the new protein haplotypes, recombination is a predominant force triggering diversity since the very first steps of protein evolution, causing at least 50% of the total amino acid variation observed in recently diverged proteins. Analysis of synonymous variants indicates that individual WSP protein types are subject to a very rapid turnover and that HVRs can accommodate a virtually unlimited repertoire of peptides. Overall distribution of WSP across hosts supports a non-random association of WSP with the host genus, although extensive horizontal transfer has occurred also in recent times.

Conclusions

In OMPs of vertebrate pathogens, large recombination impact, positive selection, reduced structural and compositional constraints, and extensive lateral gene transfer are considered hallmarks of evolution in response to the adaptive immune system. However, Wolbachia do not infect vertebrates. Here we predict that the rapid turnover of WSP loop motifs could aid in evading or inhibiting the invertebrate innate immune response. Overall, these features identify WSP as a strong candidate for future studies of host-Wolbachia interactions that affect establishment and persistence of this widespread endosymbiosis.

Genome wide evolutionary analyses reveal serotype specific patterns of positive selection in selected Salmonella serotypes

Background

The bacterium Salmonella enterica includes a diversity of serotypes that cause disease in humans and different animal species. Some Salmonella serotypes show a broad host range, some are host restricted and exclusively associated with one particular host, and some are associated with one particular host species, but able to cause disease in other host species and are thus considered "host adapted". Five Salmonella genome sequences, representing a broad host range serotype (Typhimurium), two host restricted serotypes (Typhi [two genomes] and Paratyphi) and one host adapted serotype (Choleraesuis) were used to identify core genome genes that show evidence for recombination and positive selection.

Results

Overall, 3323 orthologous genes were identified in all 5 Salmonella genomes analyzed. Use of four different methods to assess homologous recombination identified 270 genes that showed evidence for recombination with at least one of these methods (false discovery rate [FDR] <10%). After exclusion of genes with evidence for recombination, site and branch specific models identified 41 genes as showing evidence for positive selection (FDR <20%), including a number of genes with confirmed or likely roles in virulence and ompC, a gene encoding an outer membrane protein, which has also been found to be under positive selection in other bacteria. A total of 8, 16, 7, and 5 genes showed evidence for positive selection in Choleraesuis, Typhi, Typhimurium, and Paratyphi branch analyses, respectively. Sequencing and evolutionary analyses of four genes in an additional 42 isolates representing 23 serotypes confirmed branch specific positive selection and recombination patterns.

Conclusion

Our data show that, among the four serotypes analyzed, (i) less than 10% of Salmonella genes in the core genome show evidence for homologous recombination, (ii) a number of Salmonella genes are under positive selection, including genes that appear to contribute to virulence, and (iii) branch specific positive selection contributes to the evolution of host restricted Salmonella serotypes.

Evidence of extensive homologous recombination in the core genome of rickettsia

The important role of homologous recombination has been extensively demonstrated to be fundamental for genetic variation in bacterial genomes. In contrast to extracellular or facultative intracellular bacteria, obligate intracellular bacteria are considered to be less prone to recombination, especially for their core genomes. In Rickettsia, only antigen-related genes were identified to have experienced homologous recombination. In this study, we employed evolutionary genomic approaches to investigate the impact of recombination on the core genome of Rickettsia. Phylogenetic network and phylogenetic compatibility matrix analyses are clearly consistent with the hypothesis that recombination has occurred frequently during Rickettsia evolution. 28% of Rickettsia core genes (194 out of 690) are found to present the evidence of recombination under four independent statistical methods. Further functional classification shows that these recombination events occur across all functional categories, with a significant overrepresentation in the cell wall/membrane/envelope biogenesis, which may provide a molecular basis for the parasite adaptation to host immunity. This evolutionary genomic analysis provides insight into the substantial role of recombination in the evolution of the intracellular pathogenic bacteria Rickettsia.

Evolution and diversity of Rickettsia bacteria

Background

Rickettsia are intracellular symbionts of eukaryotes that are best known for infecting and causing serious diseases in humans and other mammals. All known vertebrate-associated Rickettsia are vectored by arthropods as part of their life-cycle, and many other Rickettsia are found exclusively in arthropods with no known secondary host. However, little is known about the biology of these latter strains. Here, we have identified 20 new strains of Rickettsia from arthropods, and constructed a multi-gene phylogeny of the entire genus which includes these new strains.

Results

We show that Rickettsia are primarily arthropod-associated bacteria, and identify several novel groups within the genus. Rickettsia do not co-speciate with their hosts but host shifts most often occur between related arthropods. Rickettsia have evolved adaptations including transmission through vertebrates and killing males in some arthropod hosts. We uncovered one case of horizontal gene transfer among Rickettsia, where a strain is a chimera from two distantly related groups, but multi-gene analysis indicates that different parts of the genome tend to share the same phylogeny.

Conclusion

Approximately 150 million years ago, Rickettsia split into two main clades, one of which primarily infects arthropods, and the other infects a diverse range of protists, other eukaryotes and arthropods. There was then a rapid radiation about 50 million years ago, which coincided with the evolution of life history adaptations in a few branches of the phylogeny. Even though Rickettsia are thought to be primarily transmitted vertically, host associations are short lived with frequent switching to new host lineages. Recombination throughout the genus is generally uncommon, although there is evidence of horizontal gene transfer. A better understanding of the evolution of Rickettsia will help in the future to elucidate the mechanisms of pathogenicity, transmission and virulence.

Rickettsia phylogenomics: unwinding the intricacies of obligate intracellular life

Background

Completed genome sequences are rapidly increasing for Rickettsia, obligate intracellular α-proteobacteria responsible for various human diseases, including epidemic typhus and Rocky Mountain spotted fever. In light of phylogeny, the establishment of orthologous groups (OGs) of open reading frames (ORFs) will distinguish the core rickettsial genes and other group specific genes (class 1 OGs or C1OGs) from those distributed indiscriminately throughout the rickettsial tree (class 2 OG or C2OGs).

Methodology/Principal Findings

We present 1823 representative (no gene duplications) and 259 non-representative (at least one gene duplication) rickettsial OGs. While the highly reductive (∼1.2 MB) Rickettsia genomes range in predicted ORFs from 872 to 1512, a core of 752 OGs was identified, depicting the essential Rickettsia genes. Unsurprisingly, this core lacks many metabolic genes, reflecting the dependence on host resources for growth and survival. Additionally, we bolster our recent reclassification of Rickettsia by identifying OGs that define the AG (ancestral group), TG (typhus group), TRG (transitional group), and SFG (spotted fever group) rickettsiae. OGs for insect-associated species, tick-associated species and species that harbor plasmids were also predicted. Through superimposition of all OGs over robust phylogeny estimation, we discern between C1OGs and C2OGs, the latter depicting genes either decaying from the conserved C1OGs or acquired laterally. Finally, scrutiny of non-representative OGs revealed high levels of split genes versus gene duplications, with both phenomena confounding gene orthology assignment. Interestingly, non-representative OGs, as well as OGs comprised of several gene families typically involved in microbial pathogenicity and/or the acquisition of virulence factors, fall predominantly within C2OG distributions.

Conclusion/Significance

Collectively, we determined the relative conservation and distribution of 14354 predicted ORFs from 10 rickettsial genomes across robust phylogeny estimation. The data, available at PATRIC (PathoSystems Resource Integration Center), provide novel information for unwinding the intricacies associated with Rickettsia pathogenesis, expanding the range of potential diagnostic, vaccine and therapeutic targets.

Symbiont-mediated protection

Despite the fact that all vertically transmitted symbionts sequester resources from their hosts and are therefore costly to maintain, there is an extraordinary diversity of them in invertebrates. Some spread through host populations by providing their hosts with fitness benefits or by manipulating host sex ratio, but some do not: their maintenance in host lineages remains an enigma. In this review, I explore the evolutionary ecology of vertically transmitted symbionts and their impact on host resistance, and provide an overview of the evidence for the three-way interactions between these symbionts, natural enemies and invertebrate hosts. A number of recent empirical and theoretical studies suggest that vertically transmitted symbionts may protect their hosts from pathogens. If this 'symbiont-mediated protection' is widespread, it is likely that vertically transmitted symbionts contribute significantly to variation in measures of invertebrate resistance to natural enemies.

Investigation of the demographic and selective forces shaping the nucleotide diversity of genes involved in nod factor signaling in Medicago truncatula

Symbiotic nitrogen-fixing rhizobia are able to trigger root deformation in their Fabaceae host plants, allowing their intracellular accommodation. They do so by delivering molecules called Nod factors. We analyzed the patterns of nucleotide polymorphism of five genes controlling early Nod factor perception and signaling in the Fabaceae Medicago truncatula to understand the selective forces shaping the evolution of these genes. We used 30 M. truncatula genotypes sampled in a genetically homogeneous region of the species distribution range. We first sequenced 24 independent loci and detected a genomewide departure from the hypothesis of neutrality and demographic equilibrium that suggests a population expansion. These data were used to estimate parameters of a simple demographic model incorporating population expansion. The selective neutrality of genes controlling Nod factor perception was then examined using a combination of two complementary neutrality tests, Tajima's D and Fay and Wu's standardized H. The joint distribution of D and H expected under neutrality was obtained under the fitted population expansion model. Only the gene DMI1, which is expected to regulate the downstream signal, shows a pattern consistent with a putative selective event. In contrast, the receptor-encoding genes NFP and NORK show no significant signatures of selection. Among the genes that we analyzed, only DMI1 should be viewed as a candidate for adaptation in the recent history of M. truncatula.

Lateral gene transfer between obligate intracellular bacteria: evidence from the Rickettsia massiliae genome

Rickettsia massiliae is a tick-borne obligate intracellular alpha-proteobacteria causing spotted fever in humans. Here, we present the sequence of its genome, comprising a 1.3-Mb circular chromosome and a 15.3-kb plasmid. The chromosome exhibits long-range colinearity with the other Spotted Fever Group Rickettsia genomes, except for a large fragment specific to R. massiliae that contains 14 tra genes presumably involved in pilus formation and conjugal DNA transfer. We demonstrate that the tra region was acquired recently by lateral gene transfer (LGT) from a species related to Rickettsia bellii. Further analysis of the genomic sequences identifies additional candidates of LGT between Rickettsia. Our study indicates that recent LGT between obligate intracellular Rickettsia is more common than previously thought.

Evolution of Chlamydia trachomatis diversity occurs by widespread interstrain recombination involving hotspots

Chlamydia trachomatis is an obligate intracellular bacterium of major public health significance, infecting over one-tenth of the world's population and causing blindness and infertility in millions. Mounting evidence supports recombination as a key source of genetic diversity among free-living bacteria. Previous research shows that intracellular bacteria such as Chlamydiaceae may also undergo recombination but whether this plays a significant evolutionary role has not been determined. Here, we examine multiple loci dispersed throughout the chromosome to determine the extent and significance of recombination among 19 laboratory reference strains and 10 present-day ocular and urogenital clinical isolates using phylogenetic reconstructions, compatibility matrices, and statistically based recombination programs. Recombination is widespread; all clinical isolates are recombinant at multiple loci with no two belonging to the same clonal lineage. Several reference strains show nonconcordant phylogenies across loci; one strain is unambiguously identified as recombinantly derived from other reference strain lineages. Frequent recombination contrasts with a low level of point substitution; novel substitutions relative to reference strains occur less than one per kilobase. Hotspots for recombination are identified downstream from ompA, which encodes the major outer membrane protein. This widespread recombination, unexpected for an intracellular bacterium, explains why strain-typing using one or two genes, such as ompA, does not correlate with clinical phenotypes. Our results do not point to specific events that are responsible for different pathogenicities but, instead, suggest a new approach to dissect the genetic basis for clinical strain pathology with implications for evolution, host cell adaptation, and emergence of new chlamydial diseases.