Mutation-Driven Divergence and Convergence Indicate Adaptive Evolution of the Intracellular Human-Restricted Pathogen, Bartonella bacilliformis

Single-nucleotide polymorphisms in ialB, gltA and rpoB genes of Bartonella bacilliformis isolated from patients in endemic Peruvian regions

Bartonella bacilliformis is a Gram-negative, aerobic bacterium and the known causal agent of Carrion's disease, still considered a neglected disease. There is limited information about the nucleotide sequences of this bacterium in international databases, and few studies have addressed the genetic diversity of B. bacilliformis. We analyzed a total of 20 isolates of B. bacilliformis from the Peruvian regions of Ancash and Cajamarca. Three genes (ialB, gltA, and rpoB) were sequenced in each isolate and nucleotide sequences retrieved from GenBank (16 B. bacilliformis genomes) were also included in the study. All this information was merged in order to obtain clearer evidence of the phylogenetic relationships of B. bacilliformis. In the phylogenetic analysis conducted with the concatenated markers, four isolates (B.b-1, B. b-3, B. b- 7, B.b-8) from the Ancash region were observed to form a subgroup different from B. bacilliformis type strain KC583, showing dissimilarity levels of 5.96% (ialB), 3.69% (gltA) and 3.04% (rpoB). Our results suggest that B. bacilliformis consists of two different subgroups. Future investigations are needed to establish the taxonomic status of these subgroups.

Immunogenic Peptides from Pap31 and SCS-α of Bartonella bacilliformis: One Step Closer to a Rapid Diagnostic Tool for Carrion's Disease

Bartonella bacilliformis is the causal agent of Carrion’s disease, an overlooked illness endemic in the Andean Mountains with Peru being the most affected country. The diagnostic of this illness is a challenge due to the limited resources and the common symptomatology with other infectious diseases. The goal of this study was to identify immunogenic peptides from Pap31 and succinyl-CoA synthetase α (SCS-α) of B. bacilliformis that might be suitable for developing a serologic tool. The immunodominant character of Pap31 and SCS-α was determined by Western blotting and in-silico analysis. Subsequently, 35 peptides were selected for epitope mapping and their immunoreactivity was tested by enzyme-linked immunosorbent assay (ELISA). A total of 30 sera were tested including pre-exposed people with high IgM levels for Pap31/SCS-α (23 sera), patients (2 sera) as well as 5 sera with no reactivity to Pap31/SCS-α. The results indicate that Pap31-8 (₁₈₇QAIGSAILKGTKDTGT₂₀₂) and SCS-α-12 (₅₉IFASVAEGKEKTGANA₇₄) are the most immunogenic peptides, with Pap31-8 showing potential to discriminate between B. bacilliformis and the remaining Bartonella spp., and SCS-α-12 differentiating Bartonella spp. from other microorganisms.

Molecular Characterization of Fluoroquinolone-Resistant Bartonella bacilliformis

The presence of amino acid changes in GyrA, GyrB, ParC, ParE, and in a proposed chromosomal chloramphenicol acetyl transferase (CAT), as well as mutations at 23S rRNA, were established by PCR and sequencing in 38 B. bacilliformis clinical isolates from four different areas in Peru. Eighteen out of 24 (75%) isolates showing ciprofloxacin resistance for both disk-diffusion and e-test presented amino acid substitutions in GyrA (G₈₉C, six isolates, A₉₁V, 1 isolate) GyrB (S₄₇₄F, 10 isolates) or both (GyrA D₉₅N and GyrB S₄₇₄F, one isolate). Two out of 14 susceptible isolates presented amino acid substitutions at GyrB (S₄₇₄F) or a double substitution GyrA D₉₅N and GyrB S₄₇₄F. Of note, ciprofloxacin-resistant isolates were recovered in the four areas studied. No amino acid change was observed at ParC or ParE. Only one isolate showed chloramphenicol resistance, but no alteration was present in either 23S rRNA or CAT. B. bacilliformis resistant to quinolones are extended throughout Peru, with amino acid substitutions at GyrA or GyrB as the main, albeit not exclusive, cause. B. bacilliformis seems to have an apparent facility to develop mutations on GyrB outside the classical positions 91, 95 of GyrA and 85, 88 of ParC.

Gain and loss of antibiotic resistant genes in multidrug resistant bacteria: One Health perspective

The emergence of multidrug resistance (MDR) has become a global health threat due to the increasing unnecessary use of antibiotics. Multidrug resistant bacteria occur mainly by accumulating resistance genes on mobile genetic elements (MGEs), made possible by horizontal gene transfer (HGT). Humans and animal guts along with natural and engineered environments such as wastewater treatment plants and manured soils have proven to be the major reservoirs and hotspots of spreading antibiotic resistance genes (ARGs). As those environments support the dissemination of MGEs through the complex interactions that take place at the human-animal-environment interfaces, a growing One Health challenge is for multiple sectors to communicate and work together to prevent the emergence and spread of MDR bacteria. However, maintenance of ARGs in a bacterial chromosome and/or plasmids in the environments might place energy burdens on bacterial fitness in the absence of antibiotics, and those unnecessary ARGs could eventually be lost. This review highlights and summarizes the current investigations into the gain and loss of ARG genes in MDR bacteria among human-animal-environment interfaces. We also suggest alternative treatments such as combinatory therapies or sequential use of different classes of antibiotics/adjuvants, treatment with enzyme-inhibitors, and phage therapy with antibiotics to solve the MDR problem from the perspective of One Health issues.

In silico analysis of Pap31 from Bartonella bacilliformis and other Bartonella spp

Pap31 is an outer membrane protein of Bartonella bacilliformis which is considered to be a potential antigenic candidate for the development of diagnostic tools. The present study aimed to compare Pap31 from B. bacilliformis with that of other Bartonella spp. The results showed the presence of at least 5 different B. bacilliformis Pap31 alleles, with the strain Ver097 being the most divergent (89.7% of identity with the reference strain KC583). The most significant finding was the presence of a variable number (1 to 3) of 6 amino acid tandem repeats (GTEGGG) in the different B. bacilliformis Pap31 alleles, with no similar structure in other established Bartonella spp., except for Bartonella ancashensis, another Bartonella spp. isolated from chronic cases of Carrion's disease. In both B. bacilliformis and B. ancashensis this repetitive region was coincident with the most predicted immunogenic region of the protein. In other microorganisms, the presence of amino acid tandem repeats has been related to the development of poorly functional antibodies. The findings of this study also suggest a utility of Pap31 amino acid tandem repeats as potential contributors to the immune evasion of Carrion's disease-related Bartonella spp. and the establishment of asymptomatic B. bacilliformis / B. ancashensis infections.

Gene duplication and deletion, not horizontal transfer, drove intra-species mosaicism of Bartonella henselae

Bartonella henselae is a facultative intracellular pathogen that occurs worldwide and is responsible primarily for cat-scratch disease in young people and bacillary angiomatosis in immunocompromised patients. The principal source of genome-level diversity that contributes to B. henselae's host-adaptive features is thought to be horizontal gene transfer events. However, our analyses did not reveal the acquisition of horizontally-transferred islands in B. henselae after its divergence from other Bartonella. Rather, diversity in gene content and genome size was apparently acquired through two alternative mechanisms, including deletion and, more predominantly, duplication of genes. Interestingly, a majority of these events occurred in regions that were horizontally transferred long before B. henselae's divergence from other Bartonella species. Our study indicates the possibility that gene duplication, in response to positive selection pressures in specific clones of B. henselae, might be linked to the pathogen's adaptation to arthropod vectors, the cat reservoir, or humans as incidental host-species.

Mutational convergence acts as a major player in adaptive parallel evolution of Shigella spp

Shigella spp., emerging from multiple origins of Escherichia coli, poses a significant health threat as a causative agent of bacillary dysentery. While multiple serotypes of four different species have evolved via independent lineages, Shigella spp. are designated as a single pathotype, primarily because of their common mode of pathogenesis. Convergent horizontal transfer events have so far been attributed to the commonalities in the evolution of virulence across diverse lineages. However, the role of mutational convergence in such parallel evolution is not yet well understood. Here we have carried out a genome-wide analysis of Shigella strains from all four species to detect the core genes (i.e. the ones present in all analyzed strains) acquiring convergent mutations of evolutionarily recent origin. Simulation studies show non-neutral accumulation of these convergent mutations across species, suggesting their adaptive role in the evolution of Shigella virulence. S. dysenteriae strain 197, representing highly virulent type 1 (Sd1) clone, carries excessively high number of core genes with recent convergent mutations compared to other analyzed strains. We propose that this high frequency of adaptive convergence in S. dysenteriae strain 197 could be linked to recent re-emergence of the Sd1 clone and its increased resistance to antimicrobials.

Recombination-independent rapid convergent evolution of the gastric pathogen Helicobacter pylori

BACKGROUND

Helicobacter pylori is a human stomach pathogen, naturally-competent for DNA uptake, and prone to homologous recombination. Extensive homoplasy (i.e., phylogenetically-unlinked identical variations) observed in H. pylori genes is considered a hallmark of such recombination. However, H. pylori also exhibits a high mutation rate. The relative adaptive role of homologous recombination and mutation in species diversity is a highly-debated issue in biology. Recombination results in homoplasy. While convergent mutation can also account for homoplasy, its contribution is thought to be minor. We demonstrate here that, contrary to dogma, convergent mutation is a key contributor to Helicobacter pylori homoplasy, potentially driven by adaptive evolution of proteins.

RESULTS

Our present genome-wide analysis shows that homoplastic nonsynonymous (amino acid replacement) changes are not typically accompanied by homoplastic synonymous (silent) variations. Moreover, the majority of the codon positions with homoplastic nonsynonymous changes also contain different (i.e. non-homoplastic) nonsynonymous changes arising from mutation only. This indicates that, to a considerable extent, nonsynonymous homoplasy is due to convergent mutations. High mutation rate or limited availability of evolvable sites cannot explain this excessive convergence, as suggested by our simulation studies. Rather, the genes with convergent mutations are overrepresented in distinct functional categories, suggesting possible selective responses to conditions such as distinct micro-niches in single hosts, and to differences in host genotype, physiology, habitat and diet.

CONCLUSIONS

We propose that mutational convergence is a key player in H. pylori's adaptation and extraordinary persistence in human hosts. High frequency of mutational convergence could be due to saturation of evolvable sites capable of responding to selection pressures, while the number of mutable residues is far from saturation. We anticipate a similar scenario of mutational vs. recombinational genome dynamics or plasticity for other naturally competent microbes where strong positive selection could favor frequent convergent mutations in adaptive protein evolution.

Prophage-Driven Genomic Structural Changes Promote Bartonella Vertical Evolution

Bartonella is a genetically diverse group of vector-borne bacteria. Over 40 species have been characterized to date, mainly from mammalian reservoirs and arthropod vectors. Rodent reservoirs harbor one of the largest Bartonella diversity described to date, and novel species and genetic variants are continuously identified from these hosts. Yet, it is still unknown if this significant genetic diversity stems from adaptation to different niches or from intrinsic high mutation rates. Here, we explored the vertical occurrence of spontaneous genomic alterations in 18 lines derived from two rodent-associated Bartonella elizabethae-like strains, evolved in nonselective agar plates under conditions mimicking their vector- and mammalian-associated temperatures, and the transmission cycles between them (i.e., 26 °C, 37 °C, and alterations between the two), using mutation accumulation experiments. After ∼1,000 generations, evolved genomes revealed few point mutations (average of one-point mutation per line), evidencing conserved single-nucleotide mutation rates. Interestingly, three large structural genomic changes (two large deletions and an inversion) were identified over all lines, associated with prophages and surface adhesin genes. Particularly, a prophage, deleted during constant propagation at 37 °C, was associated with an increased autonomous replication at 26 °C (the flea-associated temperature). Complementary molecular analyses of wild strains, isolated from desert rodents and their fleas, further supported the occurrence of structural genomic variations and prophage-associated deletions in nature. Our findings suggest that structural genomic changes represent an effective intrinsic mechanism to generate diversity in slow-growing bacteria and emphasize the role of prophages as promoters of diversity in nature.

PanGFR-HM: A Dynamic Web Resource for Pan-Genomic and Functional Profiling of Human Microbiome With Comparative Features

The conglomerate of microorganisms inhabiting various body-sites of human, known as the human microbiome, is one of the key determinants of human health and disease. Comprehensive pan-genomic and functional analysis approach for human microbiome components can enrich our understanding about impact of microbiome on human health. By utilizing this approach we developed PanGFR-HM (http://www.bioinfo.iicb.res.in/pangfr-hm/) – a novel dynamic web-resource that integrates genomic and functional characteristics of 1293 complete microbial genomes available from Human Microbiome Project. The resource allows users to explore genomic/functional diversity and genome-based phylogenetic relationships between human associated microbial genomes, not provided by any other resource. The key features implemented here include pan-genome and functional analysis of organisms based on taxonomy or body-site, and comparative analysis between groups of organisms. The first feature can also identify probable gene-loss events and significantly over/under represented KEGG/COG categories within pan-genome. The unique second feature can perform comparative genomic, functional and pathways analysis between 4 groups of microbes. The dynamic nature of this resource enables users to define parameters for orthologous clustering and to select any set of organisms for analysis. As an application for comparative feature of PanGFR-HM, we performed a comparative analysis with 67 Lactobacillus genomes isolated from human gut, oral cavity and urogenital tract, and therefore characterized the body-site specific genes, enzymes and pathways. Altogether, PanGFR-HM, being unique in its content and functionality, is expected to provide a platform for microbiome-based comparative functional and evolutionary genomics.

Transmission in the Origins of Bacterial Diversity, From Ecotypes to Phyla

Any two lineages, no matter how distant they are now, began their divergence as one population splitting into two lineages that could coexist indefinitely. The rate of origin of higher-level taxa is therefore the product of the rate of speciation times the probability that two new species coexist long enough to reach a particular level of divergence. Here I have explored these two parameters of disparification in bacteria. Owing to low recombination rates, sexual isolation is not a necessary milestone of bacterial speciation. Rather, irreversible and indefinite divergence begins with ecological diversification, that is, transmission of a bacterial lineage to a new ecological niche, possibly to a new microhabitat but at least to new resources. Several algorithms use sequence data from a taxon of focus to identify phylogenetic groups likely to bear the dynamic properties of species. Identifying these newly divergent lineages allows us to characterize the genetic bases of speciation, as well as the ecological dimensions upon which new species diverge. Speciation appears to be least frequent when a given lineage has few new resources it can adopt, as exemplified by photoautotrophs, C1 heterotrophs, and obligately intracellular pathogens; speciation is likely most rapid for generalist heterotrophs. The genetic basis of ecological divergence may determine whether ecological divergence is irreversible and whether lineages will diverge indefinitely into the future. Long-term coexistence is most likely when newly divergent lineages utilize at least some resources not shared with the other and when the resources themselves will coexist into the remote future.

Revisiting Bartonella bacilliformis MLST

All the studies published including Bartonella bacilliformis MLST data, as well as all B. bacilliformis genomes present in GenBank were analyzed. Overall 64 isolates and their geographical distribution were analyzed, and 14 different MLST patterns were observed. The results highlight the need for expanding the MLST studies and adding a higher number of isolates from all endemic areas.

Carrion's Disease: the Sound of Silence

Carrion's disease (CD) is a neglected biphasic vector-borne illness related to Bartonella bacilliformis. It is found in the Andean valleys and is transmitted mainly by members of the Lutzomyia genus but also by blood transfusions and from mother to child. The acute phase, Oroya fever, presents severe anemia and fever. The lethality is high in the absence of adequate treatment, despite the organism being susceptible to most antibiotics. Partial immunity is developed after infection by B. bacilliformis, resulting in high numbers of asymptomatic carriers. Following infection there is the chronic phase, Peruvian warts, involving abnormal proliferation of the endothelial cells. Despite potentially being eradicable, CD has been expanded due to human migration and geographical expansion of the vector. Moreover, in vitro studies have demonstrated the risk of the development of antimicrobial resistance. These findings, together with the description of new Bartonella species producing CD-like infections, the presence of undescribed potential vectors in new areas, the lack of adequate diagnostic tools and knowledge of the immunology and bacterial pathogenesis of CD, and poor international visibility, have led to the risk of increasing the potential expansion of resistant strains which will challenge current treatment schemes as well as the possible appearance of CD in areas where it is not endemic.