BackgroundPhylogeographic reconstruction of some bacterial populations is hindered by low diversity coupled with high levels of lateral gene transfer. A comparison of recombination levels and diversity at seven housekeeping genes for eleven bacterial species, most of which are commonly cited as having high levels of lateral gene transfer shows that the relative contributions of homologous recombination versus mutation for Burkholderia pseudomallei is over two times higher than for Streptococcus pneumoniae and is thus the highest value yet reported in bacteria. Despite the potential for homologous recombination to increase diversity, B. pseudomallei exhibits a relative lack of diversity at these loci. In these situations, whole genome genotyping of orthologous shared single nucleotide polymorphism loci, discovered using next generation sequencing technologies, can provide very large data sets capable of estimating core phylogenetic relationships. We compared and searched 43 whole genome sequences of B. pseudomallei and its closest relatives for single nucleotide polymorphisms in orthologous shared regions to use in phylogenetic reconstruction.ResultsBayesian phylogenetic analyses of >14,000 single nucleotide polymorphisms yielded completely resolved trees for these 43 strains with high levels of statistical support. These results enable a better understanding of a separate analysis of population differentiation among >1,700 B. pseudomallei isolates as defined by sequence data from seven housekeeping genes. We analyzed this larger data set for population structure and allele sharing that can be attributed to lateral gene transfer. Our results suggest that despite an almost panmictic population, we can detect two distinct populations of B. pseudomallei that conform to biogeographic patterns found in many plant and animal species. That is, separation along Wallace's Line, a biogeographic boundary between Southeast Asia and Australia.ConclusionWe describe an Australian origin for B. pseudomallei, characterized by a single introduction event into Southeast Asia during a recent glacial period, and variable levels of lateral gene transfer within populations. These patterns provide insights into mechanisms of genetic diversification in B. pseudomallei and its closest relatives, and provide a framework for integrating the traditionally separate fields of population genetics and phylogenetics for other bacterial species with high levels of lateral gene transfer.

Background: The Gram-negative bacterium Burkholderia pseudomallei (Bp) is the causative agent of the human disease melioidosis. To understand the evolutionary mechanisms contributing to Bp virulence, we performed a comparative genomic analysis of Bp K96243 and B. thailandensis (Bt) E264, a closely related but avirulent relative.

Natural isolates of Burkholderia pseudomallei (Bp), the causative agent of melioidosis, can exhibit significant ecological flexibility that is likely reflective of a dynamic genome. Using whole-genome Bp microarrays, we examined patterns of gene presence and absence across 94 South East Asian strains isolated from a variety of clinical, environmental, or animal sources. 86% of the Bp K96243 reference genome was common to all the strains representing the Bp “core genome”, comprising genes largely involved in essential functions (eg amino acid metabolism, protein translation). In contrast, 14% of the K96243 genome was variably present across the isolates. This Bp accessory genome encompassed multiple genomic islands (GIs), paralogous genes, and insertions/deletions, including three distinct lipopolysaccharide (LPS)-related gene clusters. Strikingly, strains recovered from cases of human melioidosis clustered on a tree based on accessory gene content, and were significantly more likely to harbor certain GIs compared to animal and environmental isolates. Consistent with the inference that the GIs may contribute to pathogenesis, experimental mutation of BPSS2053, a GI gene, reduced microbial adherence to human epithelial cells. Our results suggest that the Bp accessory genome is likely to play an important role in microbial adaptation and virulence.

Burkholderia pseudomallei, the etiological agent of melioidosis, is a facultative intracellular pathogen. As B. pseudomallei is a gram-negative bacterium, its outer membrane contains lipopolysaccharide (LPS) molecules, which have been shown to have low-level immunological activities in vitro. In this study, the biological activities of B. pseudomallei LPS were compared to those of Burkholderia thailandensis LPS, and it was found that both murine and human macrophages produced levels of tumor necrosis factor alpha, interleukin-6 (IL-6), and IL-10 in response to B. pseudomallei LPS that were lower than those in response to B. thailandensis LPS in vitro. In order to elucidate the molecular mechanisms underlying the low-level immunological activities of B. pseudomallei LPS, its lipid A moiety was characterized using mass spectrometry. The major lipid A species identified in B. pseudomallei consists of a biphosphorylated disaccharide backbone, which is modified with 4-amino-4-deoxy-arabinose (Ara4N) at both phosphates and penta-acylated with fatty acids (FA) C 14:0 (3-OH), C 16:0 (3-OH), and either C 14:0 or C 14:0 (2-OH). In contrast, the major lipid A species identified in B. thailandensis was a mixture of tetra-and penta-acylated structures with differing amounts of Ara4N and FA C 14:0 (3-OH). Lipid A species acylated with FA C 14:0 (2-OH) were unique to B. pseudomallei and not found in B. thailandensis. Our data thus indicate that B. pseudomallei synthesizes lipid A species with long-chain FA C 14:0 (2-OH) and Ara4N-modified phosphate groups, allowing it to evade innate immune recognition.

Many microbial pathogens express specific virulence traits at distinct growth phases. To understand the molecular pathways linking bacterial growth to pathogenicity, we have characterized the growth transcriptome of Burkholderia pseudomallei, the causative agent of melioidosis. Using a fine-scale sampling approach, we found approximately 17% of all B. pseudomallei genes displaying regulated expression during growth in rich medium, occurring as broad waves of functionally coherent gene expression tightly associated with distinct growth phases and transition points. We observed regulation of virulence genes across all growth phases and identified serC as a potentially new virulence factor by virtue of its coexpression with other early-phase virulence genes. serC-disrupted B. pseudomallei strains were serine auxotrophs and in mouse infection assays exhibited a dramatic attenuation of virulence compared to wild-type B. pseudomallei. Immunization of mice with serC-disrupted B. pseudomallei also conferred protection against subsequent challenges with different wild-type B. pseudomallei strains. At a genomic level, early-phase genes were preferentially localized on chromosome 1, while stationary-phase genes were significantly biased towards chromosome 2. We detected a significant level of chromosomally clustered gene expression, allowing us to predict ϳ100 potential operons in the B. pseudomallei genome. We computationally and experimentally validated these operons by showing that genes in these regions are preferentially transcribed in the same 533 direction, possess significantly shorter intergenic lengths than the overall genome, and are expressed as a common mRNA transcript. The availability of this transcriptome map provides an important resource for understanding the transcriptional architecture of B. pseudomallei.