Apilactobacillus kunkeei releases RNA-associated membrane vesicles and proteinaceous nanoparticles.

Extracellularly released particles, including membrane vesicles, have increasingly been recognized as important for bacterial community functions and host-interaction processes, but their compositions and functional roles differ between species and also between strains of the same species. In this study, we have determined the composition of membrane vesicles and protein particles identified in the cell-free pellets of two strains of Apilactobacillus kunkeei, a defensive symbiont of honeybees. The membrane vesicles were separated from the extracellular particles using density gradient ultracentrifugation. The peaks of the RNA and protein distributions were separated from each other and the highest concentration of RNA was observed in the fractions that contained the membrane vesicles while the highest protein concentration coincided with the fractions that contained extracellular particles. A comparative proteomics analysis by LC-MS/MS showed that 37 proteins with type-I signal peptides were consistently identified across the fractionated samples obtained from the cell-free pellets, of which 29 were orthologs detected in both strains. Functional predictions of the extracellular proteins revealed the presence of glycoside hydrolases, glycosyltransferases, giant proteins and peptidases. The extracellular transcriptomes mapped to a broad set of genes with a similar functional profile as the whole cell transcriptome. This study provides insights into the composition of membrane vesicles and extracellular proteins of a bee-associated symbiont.

Genome Evolution of a Symbiont Population for Pathogen Defense in Honeybees.

The honeybee gut microbiome is thought to be important for bee health, but the role of the individual members is poorly understood. Here, we present closed genomes and associated mobilomes of 102 Apilactobacillus kunkeei isolates obtained from the honey crop (foregut) of honeybees sampled from beehives in Helsingborg in the south of Sweden and from the islands Gotland and Åland in the Baltic Sea. Each beehive contained a unique composition of isolates and repeated sampling of similar isolates from two beehives in Helsingborg suggests that the bacterial community is stably maintained across bee generations during the summer months. The sampled bacterial population contained an open pan-genome structure with a high genomic density of transposons. A subset of strains affiliated with phylogroup A inhibited growth of the bee pathogen Melissococcus plutonius, all of which contained a 19.5 kb plasmid for the synthesis of the antimicrobial compound kunkecin A, while a subset of phylogroups B and C strains contained a 32.9 kb plasmid for the synthesis of a putative polyketide antibiotic. This study suggests that the mobile gene pool of A. kunkeei plays a key role in pathogen defense in honeybees, providing new insights into the evolutionary dynamics of defensive symbiont populations.

The genome of Rhizobiales bacteria in predatory ants reveals urease gene functions but no genes for nitrogen fixation.

Gut-associated microbiota of ants include Rhizobiales bacteria with affiliation to the genus Bartonella. These bacteria may enable the ants to fix atmospheric nitrogen, but no genomes have been sequenced yet to test the hypothesis. Sequence reads from a member of the Rhizobiales were identified in the data collected in a genome project of the ant Harpegnathos saltator. We present an analysis of the closed 1.86 Mb genome of the ant-associated bacterium, for which we suggest the species name Candidatus Tokpelaia hoelldoblerii. A phylogenetic analysis reveals a relationship to Bartonella and Brucella, which infect mammals. Novel gene acquisitions include a gene for a putative extracellular protein of more than 6,000 amino acids secreted by the type I secretion system, which may be involved in attachment to the gut epithelium. No genes for nitrogen fixation could be identified, but genes for a multi-subunit urease protein complex are present in the genome. The urease genes are also present in Brucella, which has a fecal-oral transmission pathway, but not in Bartonella, which use blood-borne transmission pathways. We hypothesize that the gain and loss of the urease function is related to transmission strategies and lifestyle changes in the host-associated members of the Rhizobiales.

A gene transfer agent and a dynamic repertoire of secretion systems hold the keys to the explosive radiation of the emerging pathogen Bartonella.

Gene transfer agents (GTAs) randomly transfer short fragments of a bacterial genome. A novel putative GTA was recently discovered in the mouse-infecting bacterium Bartonella grahamii. Although GTAs are widespread in phylogenetically diverse bacteria, their role in evolution is largely unknown. Here, we present a comparative analysis of 16 Bartonella genomes ranging from 1.4 to 2.6 Mb in size, including six novel genomes from Bartonella isolated from a cow, two moose, two dogs, and a kangaroo. A phylogenetic tree inferred from 428 orthologous core genes indicates that the deadly human pathogen B. bacilliformis is related to the ruminant-adapted clade, rather than being the earliest diverging species in the genus as previously thought. A gene flux analysis identified 12 genes for a GTA and a phage-derived origin of replication as the most conserved innovations. These are located in a region of a few hundred kb that also contains 8 insertions of gene clusters for type III, IV, and V secretion systems, and genes for putatively secreted molecules such as cholera-like toxins. The phylogenies indicate a recent transfer of seven genes in the virB gene cluster for a type IV secretion system from a cat-adapted B. henselae to a dog-adapted B. vinsonii strain. We show that the B. henselae GTA is functional and can transfer genes in vitro. We suggest that the maintenance of the GTA is driven by selection to increase the likelihood of horizontal gene transfer and argue that this process is beneficial at the population level, by facilitating adaptive evolution of the host-adaptation systems and thereby expansion of the host range size. The process counters gene loss and forces all cells to contribute to the production of the GTA and the secreted molecules. The results advance our understanding of the role that GTAs play for the evolution of bacterial genomes.

Comparative genomics of Wolbachia and the bacterial species concept.

The importance of host-specialization to speciation processes in obligate host-associated bacteria is well known, as is also the ability of recombination to generate cohesion in bacterial populations. However, whether divergent strains of highly recombining intracellular bacteria, such as Wolbachia, can maintain their genetic distinctness when infecting the same host is not known. We first developed a protocol for the genome sequencing of uncultivable endosymbionts. Using this method, we have sequenced the complete genomes of the Wolbachia strains wHa and wNo, which occur as natural double infections in Drosophila simulans populations on the Seychelles and in New Caledonia. Taxonomically, wHa belong to supergroup A and wNo to supergroup B. A comparative genomics study including additional strains supported the supergroup classification scheme and revealed 24 and 33 group-specific genes, putatively involved in host-adaptation processes. Recombination frequencies were high for strains of the same supergroup despite different host-preference patterns, leading to genomic cohesion. The inferred recombination fragments for strains of different supergroups were of short sizes, and the genomes of the co-infecting Wolbachia strains wHa and wNo were not more similar to each other and did not share more genes than other A- and B-group strains that infect different hosts. We conclude that Wolbachia strains of supergroup A and B represent genetically distinct clades, and that strains of different supergroups can co-exist in the same arthropod host without converging into the same species. This suggests that the supergroups are irreversibly separated and that barriers other than host-specialization are able to maintain distinct clades in recombining endosymbiont populations. Acquiring a good knowledge of the barriers to genetic exchange in Wolbachia will advance our understanding of how endosymbiont communities are constructed from vertically and horizontally transmitted genes.

A genome-wide study of recombination rate variation in Bartonella henselae.

BACKGROUND: Rates of recombination vary by three orders of magnitude in bacteria but the reasons for this variation is unclear. We performed a genome-wide study of recombination rate variation among genes in the intracellular bacterium Bartonella henselae, which has among the lowest estimated ratio of recombination relative to mutation in prokaryotes. RESULTS: The 1.9 Mb genomes of B. henselae strains IC11, UGA10 and Houston-1 genomes showed only minor gene content variation. Nucleotide sequence divergence levels were less than 1% and the relative rate of recombination to mutation was estimated to 1.1 for the genome overall. Four to eight segments per genome presented significantly enhanced divergences, the most pronounced of which were the virB and trw gene clusters for type IV secretion systems that play essential roles in the infection process. Consistently, multiple recombination events were identified inside these gene clusters. High recombination frequencies were also observed for a gene putatively involved in iron metabolism. A phylogenetic study of this gene in 80 strains of Bartonella quintana, B. henselae and B. grahamii indicated different population structures for each species and revealed horizontal gene transfers across Bartonella species with different host preferences. CONCLUSIONS: Our analysis has shown little novel gene acquisition in B. henselae, indicative of a closed pan-genome, but higher recombination frequencies within the population than previously estimated. We propose that the dramatically increased fixation rate for recombination events at gene clusters for type IV secretion systems is driven by selection for sequence variability.

Genome dynamics of Bartonella grahamii in micro-populations of woodland rodents.

BACKGROUND: Rodents represent a high-risk reservoir for the emergence of new human pathogens. The recent completion of the 2.3 Mb genome of Bartonella grahamii, one of the most prevalent blood-borne bacteria in wild rodents, revealed a higher abundance of genes for host-cell interaction systems than in the genomes of closely related human pathogens. The sequence variability within the global B. grahamii population was recently investigated by multi locus sequence typing, but no study on the variability of putative host-cell interaction systems has been performed. RESULTS: To study the population dynamics of B. grahamii, we analyzed the genomic diversity on a whole-genome scale of 27 B. grahamii strains isolated from four different species of wild rodents in three geographic locations separated by less than 30 km. Even using highly variable spacer regions, only 3 sequence types were identified. This low sequence diversity contrasted with a high variability in genome content. Microarray comparative genome hybridizations identified genes for outer surface proteins, including a repeated region containing the fha gene for filamentous hemaggluttinin and a plasmid that encodes a type IV secretion system, as the most variable. The estimated generation times in liquid culture medium for a subset of strains ranged from 5 to 22 hours, but did not correlate with sequence type or presence/absence patterns of the fha gene or the plasmid. CONCLUSION: Our study has revealed a geographic microstructure of B. grahamii in wild rodents. Despite near-identity in nucleotide sequence, major differences were observed in gene presence/absence patterns that did not segregate with host species. This suggests that genetically similar strains can infect a range of different hosts.

Run-off replication of host-adaptability genes is associated with gene transfer agents in the genome of mouse-infecting Bartonella grahamii.

The genus Bartonella comprises facultative intracellular bacteria adapted to mammals, including previously recognized and emerging human pathogens. We report the 2,341,328 bp genome sequence of Bartonella grahamii, one of the most prevalent Bartonella species in wild rodents. Comparative genomics revealed that rodent-associated Bartonella species have higher copy numbers of genes for putative host-adaptability factors than the related human-specific pathogens. Many of these gene clusters are located in a highly dynamic region of 461 kb. Using hybridization to a microarray designed for the B. grahamii genome, we observed a massive, putatively phage-derived run-off replication of this region. We also identified a novel gene transfer agent, which packages the bacterial genome, with an over-representation of the amplified DNA, in 14 kb pieces. This is the first observation associating the products of run-off replication with a gene transfer agent. Because of the high concentration of gene clusters for host-adaptation proteins in the amplified region, and since the genes encoding the gene transfer agent and the phage origin are well conserved in Bartonella, we hypothesize that these systems are driven by selection. We propose that the coupling of run-off replication with gene transfer agents promotes diversification and rapid spread of host-adaptability factors, facilitating host shifts in Bartonella.

The mosaic genome structure of the Wolbachia wRi strain infecting Drosophila simulans.

The obligate intracellular bacterium Wolbachia pipientis infects around 20% of all insect species. It is maternally inherited and induces reproductive alterations of insect populations by male killing, feminization, parthenogenesis, or cytoplasmic incompatibility. Here, we present the 1,445,873-bp genome of W. pipientis strain wRi that induces very strong cytoplasmic incompatibility in its natural host Drosophila simulans. A comparison with the previously sequenced genome of W. pipientis strain wMel from Drosophila melanogaster identified 35 breakpoints associated with mobile elements and repeated sequences that are stable in Drosophila lines transinfected with wRi. Additionally, 450 genes with orthologs in wRi and wMel were sequenced from the W. pipientis strain wUni, responsible for the induction of parthenogenesis in the parasitoid wasp Muscidifurax uniraptor. The comparison of these A-group Wolbachia strains uncovered the most highly recombining intracellular bacterial genomes known to date. This was manifested in a 500-fold variation in sequence divergences at synonymous sites, with different genes and gene segments supporting different strain relationships. The substitution-frequency profile resembled that of Neisseria meningitidis, which is characterized by rampant intraspecies recombination, rather than that of Rickettsia, where genes mostly diverge by nucleotide substitutions. The data further revealed diversification of ankyrin repeat genes by short tandem duplications and provided examples of horizontal gene transfer across A- and B-group strains that infect D. simulans. These results suggest that the transmission dynamics of Wolbachia and the opportunity for coinfections have created a freely recombining intracellular bacterial community with mosaic genomes.

Genome rearrangements, deletions, and amplifications in the natural population of Bartonella henselae.

Cats are the natural host for Bartonella henselae, an opportunistic human pathogen and the agent of cat scratch disease. Here, we have analyzed the natural variation in gene content and genome structure of 38 Bartonella henselae strains isolated from cats and humans by comparative genome hybridizations to microarrays and probe hybridizations to pulsed-field gel electrophoresis (PFGE) blots. The variation in gene content was modest and confined to the prophage and the genomic islands, whereas the PFGE analyses indicated extensive rearrangements across the terminus of replication with breakpoints in areas of the genomic islands. We observed no difference in gene content or structure between feline and human strains. Rather, the results suggest multiple sources of human infection from feline B. henselae strains of diverse genotypes. Additionally, the microarray hybridizations revealed DNA amplification in some strains in the so-called chromosome II-like region. The amplified segments were centered at a position corresponding to a putative phage replication initiation site and increased in size with the duration of cultivation. We hypothesize that the variable gene pool in the B. henselae population plays an important role in the establishment of long-term persistent infection in the natural host by promoting antigenic variation and escape from the host immune response.

Characterization of the genome composition of Bartonella koehlerae by microarray comparative genomic hybridization profiling.

Bartonella henselae is present in a wide range of wild and domestic feline hosts and causes cat-scratch disease and bacillary angiomatosis in humans. We have estimated here the gene content of Bartonella koehlerae, a novel species isolated from cats that was recently identified as an agent of human endocarditis. The investigation was accomplished by comparative genomic hybridization (CGH) to a microarray constructed from the sequenced 1.93-Mb genome of B. henselae. Control hybridizations of labeled DNA from the human pathogen Bartonella quintana with a reduced genome of 1.58 Mb were performed to evaluate the accuracy of the array for genes with known levels of sequence divergence. Genome size estimates of B. koehlerae by pulsed-field gel electrophoresis matched that calculated by the CGH, indicating a genome of 1.7 to 1.8 Mb with few unique genes. As in B. quintana, sequences in the prophage and the genomic islands were reported absent in B. koehlerae. In addition, sequence variability was recorded in the chromosome II-like region, where B. koehlerae showed an intermediate retention pattern of both coding and noncoding sequences. Although most of the genes missing in B. koehlerae are also absent from B. quintana, its phylogenetic placement near B. henselae suggests independent deletion events, indicating that host specificity is not solely attributed to genes in the genomic islands. Rather, the results underscore the instability of the genomic islands even within bacterial populations adapted to the same host-vector system, as in the case of B. henselae and B. koehlerae.