Continuing evolution of Burkholderia mallei through genome reduction and large-scale rearrangements

SENCA: A Multilayered Codon Model to Study the Origins and Dynamics of Codon Usage

Gene sequences are the target of evolution operating at different levels, including the nucleotide, codon, and amino acid levels. Disentangling the impact of those different levels on gene sequences requires developing a probabilistic model with three layers. Here we present SENCA (site evolution of nucleotides, codons, and amino acids), a codon substitution model that separately describes 1) nucleotide processes which apply on all sites of a sequence such as the mutational bias, 2) preferences between synonymous codons, and 3) preferences among amino acids. We argue that most synonymous substitutions are not neutral and that SENCA provides more accurate estimates of selection compared with more classical codon sequence models. We study the forces that drive the genomic content evolution, intraspecifically in the core genome of 21 prokaryotes and interspecifically for five Enterobacteria. We retrieve the existence of a universal mutational bias toward AT, and that taking into account selection on synonymous codon usage has consequences on the measurement of selection on nonsynonymous substitutions. We also confirm that codon usage bias is mostly driven by selection on preferred codons. We propose new summary statistics to measure the relative importance of the different evolutionary processes acting on sequences.

Antibiotic resistance in Burkholderia species

The genus Burkholderia comprises metabolically diverse and adaptable Gram-negative bacteria, which thrive in often adversarial environments. A few members of the genus are prominent opportunistic pathogens. These include Burkholderia mallei and Burkholderia pseudomallei of the B. pseudomallei complex, which cause glanders and melioidosis, respectively. Burkholderia cenocepacia, Burkholderia multivorans, and Burkholderia vietnamiensis belong to the Burkholderia cepacia complex and affect mostly cystic fibrosis patients. Infections caused by these bacteria are difficult to treat because of significant antibiotic resistance. The first line of defense against antimicrobials in Burkholderia species is the outer membrane penetration barrier. Most Burkholderia contain a modified lipopolysaccharide that causes intrinsic polymyxin resistance. Contributing to reduced drug penetration are restrictive porin proteins. Efflux pumps of the resistance nodulation cell division family are major players in Burkholderia multidrug resistance. Third and fourth generation β-lactam antibiotics are seminal for treatment of Burkholderia infections, but therapeutic efficacy is compromised by expression of several β-lactamases and ceftazidime target mutations. Altered DNA gyrase and dihydrofolate reductase targets cause fluoroquinolone and trimethoprim resistance, respectively. Although antibiotic resistance hampers therapy of Burkholderia infections, the characterization of resistance mechanisms lags behind other non-enteric Gram-negative pathogens, especially ESKAPE bacteria such as Acinetobacter baumannii, Klebsiella pneumoniae and Pseudomonas aeruginosa.

Phenotypic Characterization of a Novel Virulence-Factor Deletion Strain of Burkholderia mallei That Provides Partial Protection against Inhalational Glanders in Mice

Burkholderia mallei (Bm) is a highly infectious intracellular pathogen classified as a category B biological agent by the Centers for Disease Control and Prevention. After respiratory exposure, Bm establishes itself within host macrophages before spreading into major organ systems, which can lead to chronic infection, sepsis, and death. Previously, we combined computational prediction of host-pathogen interactions with yeast two-hybrid experiments and identified novel virulence factor genes in Bm, including BMAA0553, BMAA0728 (tssN), and BMAA1865. In the present study, we used recombinant allelic exchange to construct deletion mutants of BMAA0553 and tssN (ΔBMAA0553 and ΔTssN, respectively) and showed that both deletions completely abrogated virulence at doses of >100 times the LD₅₀ of the wild-type Bm strain. Analysis of ΔBMAA0553- and ΔTssN-infected mice showed starkly reduced bacterial dissemination relative to wild-type Bm, and subsequent in vitro experiments characterized pathogenic phenotypes with respect to intracellular growth, macrophage uptake and phagosomal escape, actin-based motility, and multinucleated giant cell formation. Based on observed in vitro and in vivo phenotypes, we explored the use of ΔTssN as a candidate live-attenuated vaccine. Mice immunized with aerosolized ΔTssN showed a 21-day survival rate of 67% after a high-dose aerosol challenge with the wild-type Bm ATCC 23344 strain, compared to a 0% survival rate for unvaccinated mice. However, analysis of histopathology and bacterial burden showed that while the surviving vaccinated mice were protected from acute infection, Bm was still able to establish a chronic infection. Vaccinated mice showed a modest IgG response, suggesting a limited potential of ΔTssN as a vaccine candidate, but also showed prolonged elevation of pro-inflammatory cytokines, underscoring the role of cellular and innate immunity in mitigating acute infection in inhalational glanders.

A Unique Set of the Burkholderia Collagen-Like Proteins Provides Insight into Pathogenesis, Genome Evolution and Niche Adaptation, and Infection Detection

Burkholderia pseudomallei and Burkholderia mallei, classified as category B priority pathogens, are significant human and animal pathogens that are highly infectious and broad-spectrum antibiotic resistant. Currently, the pathogenicity mechanisms utilized by Burkholderia are not fully understood, and correct diagnosis of B. pseudomallei and B. mallei infection remains a challenge due to limited detection methods. Here, we provide a comprehensive analysis of a set of 13 novel Burkholderia collagen-like proteins (Bucl) that were identified among B. pseudomallei and B. mallei select agents. We infer that several Bucl proteins participate in pathogenesis based on their noncollagenous domains that are associated with the components of a type III secretion apparatus and membrane transport systems. Homology modeling of the outer membrane efflux domain of Bucl8 points to a role in multi-drug resistance. We determined that bucl genes are widespread in B. pseudomallei and B. mallei; Fischer’s exact test and Cramer’s V² values indicate that the majority of bucl genes are highly associated with these pathogenic species versus nonpathogenic B. thailandensis. We designed a bucl-based quantitative PCR assay which was able to detect B. pseudomallei infection in a mouse with a detection limit of 50 CFU. Finally, chromosomal mapping and phylogenetic analysis of bucl loci revealed considerable genomic plasticity and adaptation of Burkholderia spp. to host and environmental niches. In this study, we identified a large set of phylogenetically unrelated bucl genes commonly found in Burkholderia select agents, encoding predicted pathogenicity factors, detection targets, and vaccine candidates.

Microbial species delineation using whole genome sequences

Increased sequencing of microbial genomes has revealed that prevailing prokaryotic species assignments can be inconsistent with whole genome information for a significant number of species. The long-standing need for a systematic and scalable species assignment technique can be met by the genome-wide Average Nucleotide Identity (gANI) metric, which is widely acknowledged as a robust measure of genomic relatedness. In this work, we demonstrate that the combination of gANI and the alignment fraction (AF) between two genomes accurately reflects their genomic relatedness. We introduce an efficient implementation of AF,gANI and discuss its successful application to 86.5M genome pairs between 13,151 prokaryotic genomes assigned to 3032 species. Subsequently, by comparing the genome clusters obtained from complete linkage clustering of these pairs to existing taxonomy, we observed that nearly 18% of all prokaryotic species suffer from anomalies in species definition. Our results can be used to explore central questions such as whether microorganisms form a continuum of genetic diversity or distinct species represented by distinct genetic signatures. We propose that this precise and objective AF,gANI-based species definition: the MiSI (Microbial Species Identifier) method, be used to address previous inconsistencies in species classification and as the primary guide for new taxonomic species assignment, supplemented by the traditional polyphasic approach, as required.

Efflux pump-mediated drug resistance in Burkholderia

Several members of the genus Burkholderia are prominent pathogens. Infections caused by these bacteria are difficult to treat because of significant antibiotic resistance. Virtually all Burkholderia species are also resistant to polymyxin, prohibiting use of drugs like colistin that are available for treatment of infections caused by most other drug resistant Gram-negative bacteria. Despite clinical significance and antibiotic resistance of Burkholderia species, characterization of efflux pumps lags behind other non-enteric Gram-negative pathogens such as Acinetobacter baumannii and Pseudomonas aeruginosa. Although efflux pumps have been described in several Burkholderia species, they have been best studied in Burkholderia cenocepacia and B. pseudomallei. As in other non-enteric Gram-negatives, efflux pumps of the resistance nodulation cell division (RND) family are the clinically most significant efflux systems in these two species. Several efflux pumps were described in B. cenocepacia, which when expressed confer resistance to clinically significant antibiotics, including aminoglycosides, chloramphenicol, fluoroquinolones, and tetracyclines. Three RND pumps have been characterized in B. pseudomallei, two of which confer either intrinsic or acquired resistance to aminoglycosides, macrolides, chloramphenicol, fluoroquinolones, tetracyclines, trimethoprim, and in some instances trimethoprim+sulfamethoxazole. Several strains of the host-adapted B. mallei, a clone of B. pseudomallei, lack AmrAB-OprA, and are therefore aminoglycoside and macrolide susceptible. B. thailandensis is closely related to B. pseudomallei, but non-pathogenic to humans. Its pump repertoire and ensuing drug resistance profile parallels that of B. pseudomallei. An efflux pump in B. vietnamiensis plays a significant role in acquired aminoglycoside resistance. Summarily, efflux pumps are significant players in Burkholderia drug resistance.

Mining host-pathogen protein interactions to characterize Burkholderia mallei infectivity mechanisms

Burkholderia pathogenicity relies on protein virulence factors to control and promote bacterial internalization, survival, and replication within eukaryotic host cells. We recently used yeast two-hybrid (Y2H) screening to identify a small set of novel Burkholderia proteins that were shown to attenuate disease progression in an aerosol infection animal model using the virulent Burkholderia mallei ATCC 23344 strain. Here, we performed an extended analysis of primarily nine B. mallei virulence factors and their interactions with human proteins to map out how the bacteria can influence and alter host processes and pathways. Specifically, we employed topological analyses to assess the connectivity patterns of targeted host proteins, identify modules of pathogen-interacting host proteins linked to processes promoting infectivity, and evaluate the effect of crosstalk among the identified host protein modules. Overall, our analysis showed that the targeted host proteins generally had a large number of interacting partners and interacted with other host proteins that were also targeted by B. mallei proteins. We also introduced a novel Host-Pathogen Interaction Alignment (HPIA) algorithm and used it to explore similarities between host-pathogen interactions of B. mallei, Yersinia pestis, and Salmonella enterica. We inferred putative roles of B. mallei proteins based on the roles of their aligned Y. pestis and S. enterica partners and showed that up to 73% of the predicted roles matched existing annotations. A key insight into Burkholderia pathogenicity derived from these analyses of Y2H host-pathogen interactions is the identification of eukaryotic-specific targeted cellular mechanisms, including the ubiquitination degradation system and the use of the focal adhesion pathway as a fulcrum for transmitting mechanical forces and regulatory signals. This provides the mechanisms to modulate and adapt the host-cell environment for the successful establishment of host infections and intracellular spread.

Burkholderia species need to manipulate many host processes and pathways in order to establish a successful intracellular infection in eukaryotic host organisms. Burkholderia mallei uses secreted virulence factor proteins as a means to execute host-pathogen interactions and promote pathogenesis. While validated virulence factor proteins have been shown to attenuate infection in animal models, their actual roles in modifying and influencing host processes are not well understood. Here, we used host-pathogen protein-protein interactions derived from yeast two-hybrid screens to study nine known B. mallei virulence factors and map out potential virulence mechanisms. From the data, we derived both general and specific insights into Burkholderia host-pathogen infectivity pathways. We showed that B. mallei virulence factors tended to target multifunctional host proteins, proteins that interacted with each other, and host proteins with a large number of interacting partners. We also identified similarities between host-pathogen interactions of B. mallei, Yersinia pestis, and Salmonella enterica using a novel host-pathogen interactions alignment algorithm. Importantly, our data are compatible with a framework in which multiple B. mallei virulence factors broadly influence key host processes related to ubiquitin-mediated proteolysis and focal adhesion. This provides B. mallei the means to modulate and adapt the host-cell environment to advance infection.

The Identification and Differentiation between Burkholderia mallei and Burkholderia pseudomallei Using One Gene Pyrosequencing

The etiologic agents for melioidosis and glanders, Burkholderia mallei and Burkholderia pseudomallei respectively, are genetically similar making identification and differentiation from other Burkholderia species and each other challenging. We used pyrosequencing to determine the presence or absence of an insertion sequence IS407A within the flagellin P (fliP) gene and to exploit the difference in orientation of this gene in the two species. Oligonucleotide primers were designed to selectively target the IS407A-fliP interface in B. mallei and the fliP gene specifically at the insertion point in B. pseudomallei. We then examined DNA from ten B. mallei, ten B. pseudomallei, 14 B. cepacia, eight other Burkholderia spp., and 17 other bacteria. Resultant pyrograms encompassed the target sequence that contained either the fliP gene with the IS407A interruption or the fully intact fliP gene with 100% sensitivity and 100% specificity. These pyrosequencing assays based upon a single gene enable investigators to reliably identify the two species. The information obtained by these assays provides more knowledge of the genomic reduction that created the new species B. mallei from B. pseudomallei and may point to new targets that can be exploited in the future.

Cross-species comparison of the Burkholderia pseudomallei, Burkholderia thailandensis, and Burkholderia mallei quorum-sensing regulons

Burkholderia pseudomallei, Burkholderia thailandensis, and Burkholderia mallei (the Bptm group) are close relatives with very different lifestyles: B. pseudomallei is an opportunistic pathogen, B. thailandensis is a nonpathogenic saprophyte, and B. mallei is a host-restricted pathogen. The acyl-homoserine lactone quorum-sensing (QS) systems of these three species show a high level of conservation. We used transcriptome sequencing (RNA-seq) to define the quorum-sensing regulon in each species, and we performed a cross-species analysis of the QS-controlled orthologs. Our analysis revealed a core set of QS-regulated genes in all three species, as well as QS-controlled factors shared by only two species or unique to a given species. This global survey of the QS regulons of B. pseudomallei, B. thailandensis, and B. mallei serves as a platform for predicting which QS-controlled processes might be important in different bacterial niches and contribute to the pathogenesis of B. pseudomallei and B. mallei.

Genome sequencing of Giardia lamblia genotypes A2 and B isolates (DH and GS) and comparative analysis with the genomes of genotypes A1 and E (WB and Pig)

Giardia lamblia (syn G. intestinalis, G. duodenalis) is the most common pathogenic intestinal parasite of humans worldwide and is a frequent cause of endemic and epidemic diarrhea. G. lamblia is divided into eight genotypes (A-H) which infect a wide range of mammals and humans, but human infections are caused by Genotypes A and B. To unambiguously determine the relationship among genotypes, we sequenced GS and DH (Genotypes B and A2) to high depth coverage and compared the assemblies with the nearly completed WB genome and draft sequencing surveys of Genotypes E (P15; pig isolate) and B (GS; human isolate). Our results identified DH as the smallest Giardia genome sequenced to date, while GS is the largest. Our open reading frame analyses and phylogenetic analyses showed that GS was more distant from the other three genomes than any of the other three were from each other. Whole-genome comparisons of DH_A2 and GS_B with the optically mapped WB_A1 demonstrated substantial synteny across all five chromosomes but also included a number of rearrangements, inversions, and chromosomal translocations that were more common toward the chromosome ends. However, the WB_A1/GS_B alignment demonstrated only about 70% sequence identity across the syntenic regions. Our findings add to information presented in previous reports suggesting that GS is a different species of Giardia as supported by the degree of genomic diversity, coding capacity, heterozygosity, phylogenetic distance, and known biological differences from WB_A1 and other G. lamblia genotypes.

Leptospiral pathogenomics

Leptospirosis, caused by pathogenic spirochetes belonging to the genus Leptospira, is a zoonosis with important impacts on human and animal health worldwide. Research on the mechanisms of Leptospira pathogenesis has been hindered due to slow growth of infectious strains, poor transformability, and a paucity of genetic tools. As a result of second generation sequencing technologies, there has been an acceleration of leptospiral genome sequencing efforts in the past decade, which has enabled a concomitant increase in functional genomics analyses of Leptospira pathogenesis. A pathogenomics approach, by coupling of pan-genomic analysis of multiple isolates with sequencing of experimentally attenuated highly pathogenic Leptospira, has resulted in the functional inference of virulence factors. The global Leptospira Genome Project supported by the U.S. National Institute of Allergy and Infectious Diseases to which key scientific contributions have been made from the international leptospirosis research community has provided a new roadmap for comprehensive studies of Leptospira and leptospirosis well into the future. This review describes functional genomics approaches to apply the data generated by the Leptospira Genome Project towards deepening our knowledge of virulence factors of Leptospira using the emerging discipline of pathogenomics.

Comparative phylogenomics and evolution of the Brucellae reveal a path to virulence

Brucella species include important zoonotic pathogens that have a substantial impact on both agriculture and human health throughout the world. Brucellae are thought of as "stealth pathogens" that escape recognition by the host innate immune response, modulate the acquired immune response, and evade intracellular destruction. We analyzed the genome sequences of members of the family Brucellaceae to assess its evolutionary history from likely free-living soil-based progenitors into highly successful intracellular pathogens. Phylogenetic analysis split the genus into two groups: recently identified and early-dividing "atypical" strains and a highly conserved "classical" core clade containing the major pathogenic species. Lateral gene transfer events brought unique genomic regions into Brucella that differentiated them from Ochrobactrum and allowed the stepwise acquisition of virulence factors that include a type IV secretion system, a perosamine-based O antigen, and systems for sequestering metal ions that are absent in progenitors. Subsequent radiation within the core Brucella resulted in lineages that appear to have evolved within their preferred mammalian hosts, restricting their virulence to become stealth pathogens capable of causing long-term chronic infections.

The Type VI secretion system spike protein VgrG5 mediates membrane fusion during intercellular spread by pseudomallei group Burkholderia species

Pseudomallei group Burkholderia species are facultative intracellular parasites that spread efficiently from cell to cell by a mechanism involving the fusion of adjacent cell membranes. Intercellular fusion requires the function of the cluster 5 type VI secretion system (T6SS-5) and its associated valine-glycine repeat protein, VgrG5. Here we show that VgrG5 alleles are conserved and functionally interchangeable between Burkholderia pseudomallei and its relatives B. mallei, B. oklahomensis, and B. thailandensis. We also demonstrate that the integrity of the VgrG5 C-terminal domain is required for fusogenic activity, and we identify sequence motifs, including two hydrophobic segments, that are important for fusion. Mutagenesis and secretion experiments using B. pseudomallei strains engineered to express T6SS-5 in vitro show that the VgrG5 C-terminal domain is dispensable for T6SS-mediated secretion of Hcp5, demonstrating that the ability of VgrG5 to mediate membrane fusion can be uncoupled from its essential role in type VI secretion. We propose a model in which a unique fusogenic activity at the C terminus of VgrG5 facilitates intercellular spread by B. pseudomallei and related species following injection across the plasma membranes of infected cells.

Novel Burkholderia mallei virulence factors linked to specific host-pathogen protein interactions

Burkholderia mallei is an infectious intracellular pathogen whose virulence and resistance to antibiotics makes it a potential bioterrorism agent. Given its genetic origin as a commensal soil organism, it is equipped with an extensive and varied set of adapted mechanisms to cope with and modulate host-cell environments. One essential virulence mechanism constitutes the specialized secretion systems that are designed to penetrate host-cell membranes and insert pathogen proteins directly into the host cell's cytosol. However, the secretion systems' proteins and, in particular, their host targets are largely uncharacterized. Here, we used a combined in silico, in vitro, and in vivo approach to identify B. mallei proteins required for pathogenicity. We used bioinformatics tools, including orthology detection and ab initio predictions of secretion system proteins, as well as published experimental Burkholderia data to initially select a small number of proteins as putative virulence factors. We then used yeast two-hybrid assays against normalized whole human and whole murine proteome libraries to detect and identify interactions among each of these bacterial proteins and host proteins. Analysis of such interactions provided both verification of known virulence factors and identification of three new putative virulence proteins. We successfully created insertion mutants for each of these three proteins using the virulent B. mallei ATCC 23344 strain. We exposed BALB/c mice to mutant strains and the wild-type strain in an aerosol challenge model using lethal B. mallei doses. In each set of experiments, mice exposed to mutant strains survived for the 21-day duration of the experiment, whereas mice exposed to the wild-type strain rapidly died. Given their in vivo role in pathogenicity, and based on the yeast two-hybrid interaction data, these results point to the importance of these pathogen proteins in modulating host ubiquitination pathways, phagosomal escape, and actin-cytoskeleton rearrangement processes.

Bifidobacterium animalis subsp. lactis ATCC 27673 is a genomically unique strain within its conserved subspecies

Many strains of Bifidobacterium animalis subsp. lactis are considered health-promoting probiotic microorganisms and are commonly formulated into fermented dairy foods. Analyses of previously sequenced genomes of B. animalis subsp. lactis have revealed little genetic diversity, suggesting that it is a monomorphic subspecies. However, during a multilocus sequence typing survey of Bifidobacterium, it was revealed that B. animalis subsp. lactis ATCC 27673 gave a profile distinct from that of the other strains of the subspecies. As part of an ongoing study designed to understand the genetic diversity of this subspecies, the genome of this strain was sequenced and compared to other sequenced genomes of B. animalis subsp. lactis and B. animalis subsp. animalis. The complete genome of ATCC 27673 was 1,963,012 bp, contained 1,616 genes and 4 rRNA operons, and had a G+C content of 61.55%. Comparative analyses revealed that the genome of ATCC 27673 contained six distinct genomic islands encoding 83 open reading frames not found in other strains of the same subspecies. In four islands, either phage or mobile genetic elements were identified. In island 6, a novel clustered regularly interspaced short palindromic repeat (CRISPR) locus which contained 81 unique spacers was identified. This type I-E CRISPR-cas system differs from the type I-C systems previously identified in this subspecies, representing the first identification of a different system in B. animalis subsp. lactis. This study revealed that ATCC 27673 is a strain of B. animalis subsp. lactis with novel genetic content and suggests that the lack of genetic variability observed is likely due to the repeated sequencing of a limited number of widely distributed commercial strains.

Less is more: Burkholderia pseudomallei and chronic melioidosis

The Gram-negative bacterium Burkholderia pseudomallei is the causative agent of melioidosis, a serious infectious disease of humans and animals. Once considered an esoteric tropical disease confined to Southeast Asia and northern Australia, research on B. pseudomallei has recently gained global prominence due to its classification as a potential bioterrorism agent by countries such as the United States and also by increasing numbers of case reports from regions where it is not endemic. An environmental bacterium typically found in soil and water, assessing the true global prevalence of melioidosis is challenged by the fact that clinical symptoms associated with B. pseudomallei infection are extremely varied and may be confused with diverse conditions such as lung cancer, tuberculosis, or Staphyloccocus aureus infection. These diagnostic challenges, coupled with lack of awareness among clinicians, have likely contributed to underdiagnosis and the high mortality rate of melioidosis, as initial treatment is often either inappropriate or delayed. Even after antibiotic treatment, relapses are frequent, and after resolution of acute symptoms, chronic melioidosis can also occur, and the symptoms can persist for months to years. In a recent article, Price et al. [mBio 4(4):e00388-13, 2013, doi:10.1128/mBio.00388-13] demonstrate how comparative genomic sequencing can reveal the repertoire of genetic changes incurred by B. pseudomallei during chronic human infection. Their results have significant clinical ramifications and highlight B. pseudomallei’s ability to survive in a wide range of potential niches within hosts, through the acquisition of genetic adaptations that optimize fitness and resource utilization.

When outgroups fail; phylogenomics of rooting the emerging pathogen, Coxiella burnetii

Rooting phylogenies is critical for understanding evolution, yet the importance, intricacies and difficulties of rooting are often overlooked. For rooting, polymorphic characters among the group of interest (ingroup) must be compared to those of a relative (outgroup) that diverged before the last common ancestor (LCA) of the ingroup. Problems arise if an outgroup does not exist, is unknown, or is so distant that few characters are shared, in which case duplicated genes originating before the LCA can be used as proxy outgroups to root diverse phylogenies. Here, we describe a genome-wide expansion of this technique that can be used to solve problems at the other end of the evolutionary scale: where ingroup individuals are all very closely related to each other, but the next closest relative is very distant. We used shared orthologous single nucleotide polymorphisms (SNPs) from 10 whole genome sequences of Coxiella burnetii, the causative agent of Q fever in humans, to create a robust, but unrooted phylogeny. To maximize the number of characters informative about the rooting, we searched entire genomes for polymorphic duplicated regions where orthologs of each paralog could be identified so that the paralogs could be used to root the tree. Recent radiations, such as those of emerging pathogens, often pose rooting challenges due to a lack of ingroup variation and large genomic differences with known outgroups. Using a phylogenomic approach, we created a robust, rooted phylogeny for C. burnetii. [Coxiella burnetii; paralog SNPs; pathogen evolution; phylogeny; recent radiation; root; rooting using duplicated genes.].

Within-host evolution of Burkholderia pseudomallei over a twelve-year chronic carriage infection

Burkholderia pseudomallei causes the potentially fatal disease melioidosis. It is generally accepted that B. pseudomallei is a noncommensal bacterium and that any culture-positive clinical specimen denotes disease requiring treatment. Over a 23-year study of melioidosis cases in Darwin, Australia, just one patient from 707 survivors has developed persistent asymptomatic B. pseudomallei carriage. To better understand the mechanisms behind this unique scenario, we performed whole-genome analysis of two strains isolated 139 months apart. During this period, B. pseudomallei underwent several adaptive changes. Of 23 point mutations, 78% were nonsynonymous and 43% were predicted to be deleterious to gene function, demonstrating a strong propensity for positive selection. Notably, a nonsense mutation inactivated the universal stress response sigma factor RpoS, with pleiotropic implications. The genome underwent substantial reduction, with four deletions in chromosome 2 resulting in the loss of 221 genes. The deleted loci included genes involved in secondary metabolism, environmental survival, and pathogenesis. Of 14 indels, 11 occurred in coding regions and 9 resulted in frameshift mutations that dramatically affected predicted gene products. Disproportionately, four indels affected lipopolysaccharide biosynthesis and modification. Finally, we identified a frameshift mutation in both P314 isolates within wcbR, an important component of the capsular polysaccharide I locus, suggesting virulence attenuation early in infection. Our study illustrates a unique clinical case that contrasts a high-consequence infectious agent with a long-term commensal infection and provides further insights into bacterial evolution within the human host.

Some bacterial pathogens establish long-term infections that are difficult or impossible to eradicate with current treatments. Rapid advances in genome sequencing technologies provide a powerful tool for understanding bacterial persistence within the human host. Burkholderia pseudomallei is considered a highly pathogenic bacterium because infection is commonly fatal. Here, we document within-host evolution of B. pseudomallei in a unique case of human infection with ongoing chronic carriage. Genomic comparison of isolates obtained 139 months (11.5 years) apart showed a strong signal of adaptation within the human host, including inactivation of virulence and immunogenic factors, and deletion of pathways involved in environmental survival. Two global regulatory genes were mutated in the 139-month isolate, indicating extensive regulatory changes favoring bacterial persistence. Our study provides insights into B. pseudomallei pathogenesis and, more broadly, identifies parallel evolutionary mechanisms that underlie chronic persistence of all bacterial pathogens.

Reductive evolution in Streptococcus agalactiae and the emergence of a host adapted lineage

BACKGROUND

During host specialization, inactivation of genes whose function is no more required is favored by changes in selective constraints and evolutionary bottlenecks. The Gram positive bacteria Streptococcus agalactiae (also called GBS), responsible for septicemia and meningitis in neonates also emerged during the seventies as a cause of severe epidemics in fish farms. To decipher the genetic basis for the emergence of these highly virulent GBS strains and of their adaptation to fish, we have analyzed the genomic sequence of seven strains isolated from fish and other poikilotherms.

RESULTS

Comparative analysis shows that the two groups of GBS strains responsible for fish epidemic diseases are only distantly related. While strains belonging to the clonal complex 7 cannot be distinguished from their human CC7 counterparts according to their gene content, strains belonging to the ST260-261 types probably diverged a long time ago. In this lineage, specialization to the fish host was correlated with a massive gene inactivation and broad changes in gene expression. We took advantage of the low level of sequence divergence between GBS strains and of the emergence of sublineages to reconstruct the different steps involved in this process. Non-homologous recombination was found to have played a major role in the genome erosion.

CONCLUSIONS

Our results show that the early phase of genome reduction during host specialization mostly involves accumulation of small and likely reversible indels, followed by a second evolutionary step marked by a higher frequency of large deletions.

φX216, a P2-like bacteriophage with broad Burkholderia pseudomallei and B. mallei strain infectivity

BACKGROUND

Burkholderia pseudomallei and B. mallei are closely related Category B Select Agents of bioterrorism and the causative agents of the diseases melioidosis and glanders, respectively. Rapid phage-based diagnostic tools would greatly benefit early recognition and treatment of these diseases. There is extensive strain-to-strain variation in B. pseudomallei genome content due in part to the presence or absence of integrated prophages. Several phages have previously been isolated from B. pseudomallei lysogens, for example φK96243, φ1026b and φ52237.

RESULTS

We have isolated a P2-like bacteriophage, φX216, which infects 78% of all B. pseudomallei strains tested. φX216 also infects B. mallei, but not other Burkholderia species, including the closely related B. thailandensis and B. oklahomensis. The nature of the φX216 host receptor remains unclear but evidence indicates that in B. mallei φX216 uses lipopolysaccharide O-antigen but a different receptor in B. pseudomallei. The 37,637 bp genome of φX216 encodes 47 predicted open reading frames and shares 99.8% pairwise identity and an identical strain host range with bacteriophage φ52237. Closely related P2-like prophages appear to be widely distributed among B. pseudomallei strains but both φX216 and φ52237 readily infect prophage carrying strains.

CONCLUSIONS

The broad strain infectivity and high specificity for B. pseudomallei and B. mallei indicate that φX216 will provide a good platform for the development of phage-based diagnostics for these bacteria.

Developing insights into the mechanisms of evolution of bacterial pathogens from whole-genome sequences

Evolution of bacterial pathogen populations has been detected in a variety of ways including phenotypic tests, such as metabolic activity, reaction to antisera and drug resistance and genotypic tests that measure variation in chromosome structure, repetitive loci and individual gene sequences. While informative, these methods only capture a small subset of the total variation and, therefore, have limited resolution. Advances in sequencing technologies have made it feasible to capture whole-genome sequence variation for each sample under study, providing the potential to detect all changes at all positions in the genome from single nucleotide changes to large-scale insertions and deletions. In this review, we focus on recent work that has applied this powerful new approach and summarize some of the advances that this has brought in our understanding of the details of how bacterial pathogens evolve.

The genetic integrity of bacterial species: the core genome and the accessory genome, two different stories

Strains within a bacterial species typically have a set of conserved core genes and a variable set of accessory genes. The accessory genes often appear to move laterally between strains, thereby forming new trait combinations. Sometimes, genetic material also moves laterally between species, thereby resulting in diffuse borders between them. The growing number of genome sequences offers new possibilities to study these processes. Ten species for which abundant genomic data exists were here selected for analysis of the species border integrity. The average core genome similarities and relative core genome sizes (RCGSs) were determined for strain pairs within the species and for strain pairs crossing the species border. The variability within the species as well as the border integrity varies for different bacterial species. Some have very distinct borders while others are more or less indefinable. From the growing amount of genomic data, it becomes even clearer that the concept of bacterial species is, in many cases, far from absolute.

Investigation of the bacterial communities associated with females of Lutzomyia sand fly species from South America

Phlebotomine sand flies are vectors of Leishmania that are acquired by the female sand fly during blood feeding on an infected mammal. Leishmania parasites develop exclusively in the gut lumen during their residence in the insect before transmission to a suitable host during the next blood feed. Female phlebotomine sand flies are blood feeding insects but their life style of visiting plants as well as animals, and the propensity for larvae to feed on detritus including animal faeces means that the insect host and parasite are exposed to a range of microorganisms. Thus, the sand fly microbiota may interact with the developing Leishmania population in the gut. The aim of the study was to investigate and identify the bacterial diversity associated with wild adult female Lutzomyia sand flies from different geographical locations in the New World. The bacterial phylotypes recovered from 16S rRNA gene clone libraries obtained from wild caught adult female Lutzomyia sand flies were estimated from direct band sequencing after denaturing gradient gel electrophoresis of bacterial 16 rRNA gene fragments. These results confirm that the Lutzomyia sand flies contain a limited array of bacterial phylotypes across several divisions. Several potential plant-related bacterial sequences were detected including Erwinia sp. and putative Ralstonia sp. from two sand fly species sampled from 3 geographically separated regions in Brazil. Identification of putative human pathogens also demonstrated the potential for sand flies to act as vectors of bacterial pathogens of medical importance in addition to their role in Leishmania transmission.

Glanders in animals: a review on epidemiology, clinical presentation, diagnosis and countermeasures

Glanders or farcy, caused by Burkholderia mallei, is an infectious and zoonotic disease of solipeds. Horses, donkeys and mules are the only known natural reservoir of B. mallei. Although glanders has been eradicated from most countries, it has regained the status of a re-emerging disease because of the numerous recent outbreaks. Pre-symptomatic or carrier animals are the potential source of infection for the healthy equine population and play a crucial role in the spreading of the infectious agent. Glanders is characterized by ulcerating nodular lesions of the skin and mucous membrane. Generalized symptoms include fever, malaise, depression, cough, anorexia and weight loss. Burkholderia mallei can invade its host through mucous membranes, gastrointestinal tract and the integument. Its virulence mechanisms and pathogenesis are not yet completely understood. A major problem when using serological tests for diagnosing glanders is the occurrence of false-positive and false-negative results leading to difficulties in international trade with equids and to the spread of glanders to disease-free regions. Moreover, poor tests critically result in poor control of disease. These tests are not only incapable of discriminating between B. mallei and B. pseudomallei antibodies, they are also unable to differentiate between malleinized and naturally infected animals. Combined use of both serological and molecular detection methods increases the detection rate of glanders. Countermeasures against glanders include early detection of disease in susceptible animals, stringent quarantine measures, testing and safe destruction of infected carcasses, adequate compensation to the animal owners, disinfection of infected premises and awareness about glanders and the zoonotic implications through veterinary extension services. An account of the clinical picture and successful experimental therapy of spontaneous equine glanders is also given.

Insights from genomic comparisons of genetically monomorphic bacterial pathogens

Some of the most deadly bacterial diseases, including leprosy, anthrax and plague, are caused by bacterial lineages with extremely low levels of genetic diversity, the so-called 'genetically monomorphic bacteria'. It has only become possible to analyse the population genetics of such bacteria since the recent advent of high-throughput comparative genomics. The genomes of genetically monomorphic lineages contain very few polymorphic sites, which often reflect unambiguous clonal genealogies. Some genetically monomorphic lineages have evolved in the last decades, e.g. antibiotic-resistant Staphylococcus aureus, whereas others have evolved over several millennia, e.g. the cause of plague, Yersinia pestis. Based on recent results, it is now possible to reconstruct the sources and the history of pandemic waves of plague by a combined analysis of phylogeographic signals in Y. pestis plus polymorphisms found in ancient DNA. Different from historical accounts based exclusively on human disease, Y. pestis evolved in China, or the vicinity, and has spread globally on multiple occasions. These routes of transmission can be reconstructed from the genealogy, most precisely for the most recent pandemic that was spread from Hong Kong in multiple independent waves in 1894.

Environmental mycobacteria: a threat to human health?

In many cases, bacterial pathogens are close relatives to nonpathogens. Pathogens seem to be limited lineages within nonpathogenic bacteria. Nonpathogenic isolates are generally more diverse and widespread in the environment and it is generally considered that environmental bacteria do not pose a risk to human health as clinical isolates do; this may not be the case with mycobacteria, but environmental mycobacteria have not been well studied. It is documented that several environmental mycobacteria constitute a source for human infections. Diverse mycobacterial environmental isolates are rarely involved in human disease. Environmental mycobacteria may have a role in degradation of different compounds. Environmental mycobacteria have had a long interaction with humans, maybe as long as the human species, and may have contributed to human evolution.

Outer membrane proteome of Burkholderia pseudomallei and Burkholderia mallei from diverse growth conditions

Burkholderia mallei and Burkholderia pseudomallei are closely related, aerosol-infective human pathogens that cause life-threatening diseases. Biochemical analyses requiring large-scale growth and manipulation at biosafety level 3 under select agent regulations are cumbersome and hazardous. We developed a simple, safe, and rapid method to prepare highly purified outer membrane (OM) fragments from these pathogens. Shotgun proteomic analyses of OMs by trypsin shaving and mass spectrometry identified >155 proteins, the majority of which are clearly outer membrane proteins (OMPs). These included: 13 porins, 4 secretins for virulence factor export, 11 efflux pumps, multiple components of a Type VI secreton, metal transport receptors, polysaccharide exporters, and hypothetical OMPs of unknown function. We also identified 20 OMPs in each pathogen that are abundant under a wide variety of conditions, including in serum and with macrophages, suggesting these are fundamental for growth and survival and may represent prime drug or vaccine targets. Comparison of the OM proteomes of B. mallei and B. pseudomallei showed many similarities but also revealed a few differences, perhaps reflecting evolution of B. mallei away from environmental survival toward host-adaptation.

Genome sequencing and analysis of Yersina pestis KIM D27, an avirulent strain exempt from select agent regulation

Yersinia pestis is the causative agent of the plague. Y. pestis KIM 10+ strain was passaged and selected for loss of the 102 kb pgm locus, resulting in an attenuated strain, KIM D27. In this study, whole genome sequencing was performed on KIM D27 in order to identify any additional differences. Initial assemblies of 454 data were highly fragmented, and various bioinformatic tools detected between 15 and 465 SNPs and INDELs when comparing both strains, the vast majority associated with A or T homopolymer sequences. Consequently, Illumina sequencing was performed to improve the quality of the assembly. Hybrid sequence assemblies were performed and a total of 56 validated SNP/INDELs and 5 repeat differences were identified in the D27 strain relative to published KIM 10+ sequence. However, further analysis showed that 55 of these SNP/INDELs and 3 repeats were errors in the KIM 10+ reference sequence. We conclude that both 454 and Illumina sequencing were required to obtain the most accurate and rapid sequence results for Y. pestis KIMD27. SNP and INDELS calls were most accurate when both Newbler and CLC Genomics Workbench were employed. For purposes of obtaining high quality genome sequence differences between strains, any identified differences should be verified in both the new and reference genomes.

Macrophage and Galleria mellonella infection models reflect the virulence of naturally occurring isolates of B. pseudomallei, B. thailandensis and B. oklahomensis

BACKGROUND

Burkholderia pseudomallei is the causative agent of melioidosis, a tropical disease of humans with a variable and often fatal outcome. In murine models of infection, different strains exhibit varying degrees of virulence. In contrast, two related species, B. thailandensis and B. oklahomensis, are highly attenuated in mice. Our aim was to determine whether virulence in mice is reflected in macrophage or wax moth larvae (Galleria mellonella) infection models.

RESULTS

B. pseudomallei strains 576 and K96243, which have low median lethal dose (MLD) values in mice, were able to replicate and induce cellular damage in macrophages and caused rapid death of G. mellonella. In contrast, B. pseudomallei strain 708a, which is attenuated in mice, showed reduced replication in macrophages, negligible cellular damage and was avirulent in G. mellonella larvae. B. thailandensis isolates were less virulent than B. pseudomallei in all of the models tested. However, we did record strain dependent differences. B. oklahomensis isolates were the least virulent isolates. They showed minimal ability to replicate in macrophages, were unable to evoke actin-based motility or to form multinucleated giant cells and were markedly attenuated in G. mellonella compared to B. thailandensis.

CONCLUSIONS

We have shown that the alternative infection models tested here, namely macrophages and Galleria mellonella, are able to distinguish between strains of B. pseudomallei, B. thailandensis and B. oklahomensis and that these differences reflect the observed virulence in murine infection models. Our results indicate that B. oklahomensis is the least pathogenic of the species investigated. They also show a correlation between isolates of B. thailandensis associated with human infection and virulence in macrophage and Galleria infection models.

Genetic and phenotypic diversity in Burkholderia: contributions by prophage and phage-like elements

BACKGROUND

Burkholderia species exhibit enormous phenotypic diversity, ranging from the nonpathogenic, soil- and water-inhabiting Burkholderia thailandensis to the virulent, host-adapted mammalian pathogen B. mallei. Genomic diversity is evident within Burkholderia species as well. Individual isolates of Burkholderia pseudomallei and B. thailandensis, for example, carry a variety of strain-specific genomic islands (GIs), including putative pathogenicity and metabolic islands, prophage-like islands, and prophages. These GIs may provide some strains with a competitive advantage in the environment and/or in the host relative to other strains.

RESULTS

Here we present the results of analysis of 37 prophages, putative prophages, and prophage-like elements from six different Burkholderia species. Five of these were spontaneously induced to form bacteriophage particles from B. pseudomallei and B. thailandensis strains and were isolated and fully sequenced; 24 were computationally predicted in sequenced Burkholderia genomes; and eight are previously characterized prophages or prophage-like elements. The results reveal numerous differences in both genome structure and gene content among elements derived from different species as well as from strains within species, due in part to the incorporation of additional DNA, or 'morons' into the prophage genomes. Implications for pathogenicity are also discussed. Lastly, RNAseq analysis of gene expression showed that many of the genes in varphi1026b that appear to contribute to phage and lysogen fitness were expressed independently of the phage structural and replication genes.

CONCLUSIONS

This study provides the first estimate of the relative contribution of prophages to the vast phenotypic diversity found among the Burkholderiae.