Apichai Tuanyok

Global and regional dissemination and evolution of Burkholderia pseudomallei

The environmental bacterium Burkholderia pseudomallei causes an estimated 165,000 cases of human melioidosis per year worldwide and is also classified as a biothreat agent. We used whole genome sequences of 469 B. pseudomallei isolates from 30 countries collected over 79 years to explore its geographic transmission. Our data point to Australia as an early reservoir, with transmission to Southeast Asia followed by onward transmission to South Asia and East Asia. Repeated reintroductions were observed within the Malay Peninsula and between countries bordered by the Mekong River. Our data support an African origin of the Central and South American isolates with introduction of B. pseudomallei into the Americas between 1650 and 1850, providing a temporal link with the slave trade. We also identified geographically distinct genes/variants in Australasian or Southeast Asian isolates alone, with virulence-associated genes being among those over-represented. This provides a potential explanation for clinical manifestations of melioidosis that are geographically restricted.

Identification of Burkholderia pseudomallei Near-Neighbor Species in the Northern Territory of Australia

Identification and characterization of near-neighbor species are critical to the development of robust molecular diagnostic tools for biothreat agents. One such agent, Burkholderia pseudomallei, a soil bacterium and the causative agent of melioidosis, is lacking in this area because of its genomic diversity and widespread geographic distribution. The Burkholderia genus contains over 60 species and occupies a large range of environments including soil, plants, rhizospheres, water, animals and humans. The identification of novel species in new locations necessitates the need to identify the true global distribution of Burkholderia species, especially the members that are closely related to B. pseudomallei. In our current study, we used the Burkholderia-specific recA sequencing assay to analyze environmental samples from the Darwin region in the Northern Territory of Australia where melioidosis is endemic. Burkholderia recA PCR negative samples were further characterized using 16s rRNA sequencing for species identification. Phylogenetic analysis demonstrated that over 70% of the bacterial isolates were identified as B. ubonensis indicating that this species is common in the soil where B. pseudomallei is endemic. Bayesian phylogenetic analysis reveals many novel branches within the B. cepacia complex, one novel B. oklahomensis-like species, and one novel branch containing one isolate that is distinct from all other samples on the phylogenetic tree. During the analysis with recA sequencing, we discovered 2 single nucleotide polymorphisms in the reverse priming region of B. oklahomensis. A degenerate primer was developed and is proposed for future use. We conclude that the recA sequencing technique is an effective tool to classify Burkholderia and identify soil organisms in a melioidosis endemic area.

Direct detection of Burkholderia pseudomallei and biological factors in soil

BACKGROUND Burkholderia pseudomallei, a Gram-negative saprophytic bacillus, is a severe infectious agent that causes melioidosis and soil is the most important reservoir. METHODS One hundred and forty soil samples were tested for pH, moisture content and total C and N measurements and used for DNA extraction and culture for B. pseudomallei. The quantitative real-time PCR (qPCR) targeting wcbG, a putative capsular polysaccharide biosynthesis protein gene of B. pseudomallei, was developed to detect the bacterium, and random amplified polymorphic DNA (RAPD) was used to detect the microbial diversity in soil. RESULTS The acidic pH was correlated with the presence of the bacterium. Forty-four soil sites (44/140, 31.4%) were positive for B. pseudomallei by qPCR, of which 21 were positive by culture. The limit of detection is 32 fg of DNA (about 4 genomes). The RAPD method could classify the soil samples into low diversity (LD) and high diversity (HD) sites. The trend of LD was found with B. pseudomallei positive soil sites. CONCLUSIONS The acidity of the soil or metabolites from organisms in the sites may contribute to the presence of the bacterium. Further investigation of microbes by a more robust method should elucidate biological factors that promote the presence of B. pseudomallei and may be used for controlling the bacterium in the environment.

Genomic Characterization of Burkholderia pseudomallei Isolates Selected for Medical Countermeasures Testing: Comparative Genomics Associated with Differential Virulence

Burkholderia pseudomallei is the causative agent of melioidosis and a potential bioterrorism agent. In the development of medical countermeasures against B. pseudomallei infection, the US Food and Drug Administration (FDA) animal Rule recommends using well-characterized strains in animal challenge studies. In this study, whole genome sequence data were generated for 6 B. pseudomallei isolates previously identified as candidates for animal challenge studies; an additional 5 isolates were sequenced that were associated with human inhalational melioidosis. A core genome single nucleotide polymorphism (SNP) phylogeny inferred from a concatenated SNP alignment from the 11 isolates sequenced in this study and a diverse global collection of isolates demonstrated the diversity of the proposed Animal Rule isolates. To understand the genomic composition of each isolate, a large-scale blast score ratio (LS-BSR) analysis was performed on the entire pan-genome; this demonstrated the variable composition of genes across the panel and also helped to identify genes unique to individual isolates. In addition, a set of ~550 genes associated with pathogenesis in B. pseudomallei were screened against the 11 sequenced genomes with LS-BSR. Differential gene distribution for 54 virulence-associated genes was observed between genomes and three of these genes were correlated with differential virulence observed in animal challenge studies using BALB/c mice. Differentially conserved genes and SNPs associated with disease severity were identified and could be the basis for future studies investigating the pathogenesis of B. pseudomallei. Overall, the genetic characterization of the 11 proposed Animal Rule isolates provides context for future studies involving B. pseudomallei pathogenesis, differential virulence, and efficacy to therapeutics.

Pangenome Analysis of Burkholderia pseudomallei: Genome Evolution Preserves Gene Order despite High Recombination Rates

The pangenomic diversity in Burkholderia pseudomallei is high, with approximately 5.8% of the genome consisting of genomic islands. Genomic islands are known hotspots for recombination driven primarily by site-specific recombination associated with tRNAs. However, recombination rates in other portions of the genome are also high, a feature we expected to disrupt gene order. We analyzed the pangenome of 37 isolates of B. pseudomallei and demonstrate that the pangenome is ‘open’, with approximately 136 new genes identified with each new genome sequenced, and that the global core genome consists of 4568±16 homologs. Genes associated with metabolism were statistically overrepresented in the core genome, and genes associated with mobile elements, disease, and motility were primarily associated with accessory portions of the pangenome. The frequency distribution of genes present in between 1 and 37 of the genomes analyzed matches well with a model of genome evolution in which 96% of the genome has very low recombination rates but 4% of the genome recombines readily. Using homologous genes among pairs of genomes, we found that gene order was highly conserved among strains, despite the high recombination rates previously observed. High rates of gene transfer and recombination are incompatible with retaining gene order unless these processes are either highly localized to specific sites within the genome, or are characterized by symmetrical gene gain and loss. Our results demonstrate that both processes occur: localized recombination introduces many new genes at relatively few sites, and recombination throughout the genome generates the novel multi-locus sequence types previously observed while preserving gene order.

Structural diversity of Burkholderia pseudomallei lipopolysaccharides affects innate immune signaling

Burkholderia pseudomallei (Bp) causes the disease melioidosis. The main cause of mortality in this disease is septic shock triggered by the host responding to lipopolysaccharide (LPS) components of the Gram-negative outer membrane. Bp LPS is thought to be a weak inducer of the host immune system. LPS from several strains of Bp were purified and their ability to induce the inflammatory mediators TNF-α and iNOS in murine macrophages at low concentrations was investigated. Innate and adaptive immunity qPCR arrays were used to profile expression patterns of 84 gene targets in response to the different LPS types. Additional qPCR validation confirmed large differences in macrophage response. LPS from a high-virulence serotype B strain 576a and a virulent rough central nervous system tropic strain MSHR435 greatly induced the innate immune response indicating that the immunopathogenesis of these strains is different than in infections with strains similar to the prototype strain 1026b. The accumulation of autophagic vesicles was also increased in macrophages challenged with highly immunogenic Bp LPS. Gene induction and concomitant cytokine secretion profiles of human PBMCs in response to the various LPS were also investigated. MALDI-TOF/TOF was used to probe the lipid A portions of the LPS, indicating substantial structural differences that likely play a role in host response to LPS. These findings add to the evolving knowledge of host-response to bacterial LPS, which can be used to better understand septic shock in melioidosis patients and in the rational design of vaccines.

Antibiotic Resistance Markers in Burkholderia pseudomallei Strain Bp1651 Identified by Genome Sequence Analysis.

ABSTRACT Burkholderia pseudomallei Bp1651 is resistant to several classes of antibiotics that are usually effective for treatment of melioidosis, including tetracyclines, sulfonamides, and β-lactams such as penicillins (amoxicillin-clavulanic acid), cephalosporins (ceftazidime), and carbapenems (imipenem and meropenem). We sequenced, assembled, and annotated the Bp1651 genome and analyzed the sequence using comparative genomic analyses with susceptible strains, keyword searches of the annotation, publicly available antimicrobial resistance prediction tools, and published reports. More than 100 genes in the Bp1651 sequence were identified as potentially contributing to antimicrobial resistance. Most notably, we identified three previously uncharacterized point mutations in penA, which codes for a class A β-lactamase and was previously implicated in resistance to β-lactam antibiotics. The mutations result in amino acid changes T147A, D240G, and V261I. When individually introduced into select agent-excluded B. pseudomallei strain Bp82, D240G was found to contribute to ceftazidime resistance and T147A contributed to amoxicillin-clavulanic acid and imipenem resistance. This study provides the first evidence that mutations in penA may alter susceptibility to carbapenems in B. pseudomallei. Another mutation of interest was a point mutation affecting the dihydrofolate reductase gene folA, which likely explains the trimethoprim resistance of this strain. Bp1651 was susceptible to aminoglycosides likely because of a frameshift in the amrB gene, the transporter subunit of the AmrAB-OprA efflux pump. These findings expand the role of penA to include resistance to carbapenems and may assist in the development of molecular diagnostics that predict antimicrobial resistance and provide guidance for treatment of melioidosis.

The Effects of Signal Erosion and Core Genome Reduction on the Identification of Diagnostic Markers

ABSTRACT Whole-genome sequence (WGS) data are commonly used to design diagnostic targets for the identification of bacterial pathogens. To do this effectively, genomics databases must be comprehensive to identify the strict core genome that is specific to the target pathogen. As additional genomes are analyzed, the core genome size is reduced and there is erosion of the target-specific regions due to commonality with related species, potentially resulting in the identification of false positives and/or false negatives. IMPORTANCE A comparative analysis of 1,130 Burkholderia genomes identified unique markers for many named species, including the human pathogens B. pseudomallei and B. mallei. Due to core genome reduction and signature erosion, only 38 targets specific to B. pseudomallei/mallei were identified. By using only public genomes, a larger number of markers were identified, due to undersampling, and this larger number represents the potential for false positives. This analysis has implications for the design of diagnostics for other species where the genomic space of the target and/or closely related species is not well defined. A comparative analysis of 1,130 Burkholderia genomes identified unique markers for many named species, including the human pathogens B. pseudomallei and B. mallei. Due to core genome reduction and signature erosion, only 38 targets specific to B. pseudomallei/mallei were identified. By using only public genomes, a larger number of markers were identified, due to undersampling, and this larger number represents the potential for false positives. This analysis has implications for the design of diagnostics for other species where the genomic space of the target and/or closely related species is not well defined.

A Burkholderia pseudomallei Colony Variant Necessary for Gastric Colonization

ABSTRACT  Diverse colony morphologies are a hallmark of Burkholderia pseudomallei recovered from infected patients. We observed that stresses that inhibit aerobic respiration shifted populations of B. pseudomallei from the canonical white colony morphotype toward two distinct, reversible, yet relatively stable yellow colony variants (YA and YB). As accumulating evidence supports the importance of B. pseudomallei enteric infection and gastric colonization, we tested the response of yellow variants to hypoxia, acidity, and stomach colonization. Yellow variants exhibited a competitive advantage under hypoxic and acidic conditions and alkalized culture media. The YB variant, although highly attenuated in acute virulence, was the only form capable of colonization and persistence in the murine stomach. The accumulation of extracellular DNA (eDNA) was a characteristic of YB as observed by 4′,6-diamidino-2-phenylindole (DAPI) staining of gastric tissues, as well as in an in vitro stomach model where large amounts of eDNA were produced without cell lysis. Transposon mutagenesis identified a transcriptional regulator (BPSL1887, designated YelR) that when overexpressed produced the yellow phenotype. Deletion of yelR blocked a shift from white to the yellow forms. These data demonstrate that YB is a unique B. pseudomallei pathovariant controlled by YelR that is specifically adapted to the harsh gastric environment and necessary for persistent stomach colonization. IMPORTANCE  Seemingly uniform populations of bacteria often contain subpopulations that are genetically identical but display unique characteristics which offer advantages when the population is faced with infrequent but predictable stresses. The pathogen Burkholderia pseudomallei is capable of forming several reversible colony types, and it interconverted between one white type and two yellow types under certain environmental stresses. The two yellow forms exhibited distinct advantages in low-oxygen and acidic environments. One yellow colony variant was the only form capable of chronic stomach colonization. Areas of gastric infection were marked by bacteria encased in a DNA matrix, and the yellow forms were able to produce large amounts of extracellular DNA in vitro. We also identified the regulator in control of yellow colony variant formation. These findings demonstrate a role in infection for colony variation and provide a mechanism for chronic stomach colonization—a frequently overlooked niche in melioidosis. Seemingly uniform populations of bacteria often contain subpopulations that are genetically identical but display unique characteristics which offer advantages when the population is faced with infrequent but predictable stresses. The pathogen Burkholderia pseudomallei is capable of forming several reversible colony types, and it interconverted between one white type and two yellow types under certain environmental stresses. The two yellow forms exhibited distinct advantages in low-oxygen and acidic environments. One yellow colony variant was the only form capable of chronic stomach colonization. Areas of gastric infection were marked by bacteria encased in a DNA matrix, and the yellow forms were able to produce large amounts of extracellular DNA in vitro. We also identified the regulator in control of yellow colony variant formation. These findings demonstrate a role in infection for colony variation and provide a mechanism for chronic stomach colonization—a frequently overlooked niche in melioidosis.

Multitarget Quantitative PCR Improves Detection and Predicts Cultivability of the Pathogen Burkholderia pseudomallei

ABSTRACT Burkholderia pseudomallei is present in the environment in many parts of the world and causes the often-fatal disease melioidosis. The sensitive detection and quantification of B. pseudomallei in the environment are a prerequisite for assessing the risk of infection. We recently reported the direct detection of B. pseudomallei in soil samples using a quantitative PCR (qPCR) targeting a single type three secretion system 1 (TTSS1) gene. Here, we extend the qPCR-based analysis of B. pseudomallei in soil by validating novel qPCR gene targets selected from a comparative genomic analysis. Two hundred soil samples from two rice paddies in northeast Thailand were evaluated, of which 47% (94/200) were B. pseudomallei culture positive. The TTSS1 qPCR and two novel qPCR assays that targeted open reading frames (ORFs) BPSS0087 and BPSS0745 exhibited detection rates of 76.5% (153/200), 34.5% (69/200), and 74.5% (150/200), respectively. The combination of TTSS1 and BPSS0745 qPCR increased the detection rate to 90% (180/200). Combining the results of the three qPCR assays and the BPSS1187 nested PCR previously published, all 200 samples were positive by at least one PCR assay. Samples positive by either TTSS1 (n = 153) or BPSS0745 (n = 150) qPCR were more likely to be direct-culture positive, with odds ratios of 4.0 (95% confidence interval [CI], 1.7 to 9.5; P < 0.001) and 9.0 (95% CI, 3.1 to 26.4; P < 0.001), respectively. High B. pseudomallei genome equivalents correlated with high CFU counts by culture. In conclusion, multitarget qPCR improved the B. pseudomallei detection rate in soil samples and predicted culture positivity. This approach has the potential for use as a sensitive environmental screening method for B. pseudomallei. IMPORTANCE The worldwide environmental distribution of the soil bacterium Burkholderia pseudomallei remains to be determined. So far, most environmental studies have relied on culture-based approaches to detect this pathogen. Since current culture methods are laborious, are time consuming, and have limited sensitivity, culture-independent and more sensitive methods are needed. In this study, we show that a B. pseudomallei-specific qPCR approach can detect significantly higher numbers of B. pseudomallei-positive soil samples from areas where it is endemic compared with that from culture. The use of multiple independent B. pseudomallei-specific qPCR targets further increased the detection rate of B. pseudomallei compared with that from single targets. Samples with a high molecular B. pseudomallei load were more likely to be culture positive. We conclude that our quantitative multitarget approach might be useful in defining areas where there is a risk of B. pseudomallei infections in different parts of the world.