In silico analyses of penicillin binding proteins in Burkholderia pseudomallei uncovers SNPs with utility for phylogeography, species differentiation, and sequence typing
Authors: McLaughlin Heather P. , Gulvik Christopher A. , Sue David 
Published: October 2021

Background

Burkholderia pseudomallei causes melioidosis. Sequence typing this pathogen can reveal geographical origin and uncover epidemiological associations. Here, we describe B. pseudomallei genes encoding putative penicillin binding proteins (PBPs) and investigate their utility for determining phylogeography and differentiating closely related species.

Methodology & Principal Findings

We performed in silico analysis to characterize 10 PBP homologs in B. pseudomallei 1026b. As PBP active site mutations can confer β-lactam resistance in Gram-negative bacteria, PBP sequences in two resistant B. pseudomallei strains were examined for similar alterations. Sequence alignments revealed single amino acid polymorphisms (SAAPs) unique to the multidrug resistant strain Bp1651 in the transpeptidase domains of two PBPs, but not directly within the active sites. Using BLASTn analyses of complete assembled genomes in the NCBI database, we determined genes encoding PBPs were conserved among B. pseudomallei (n=101) and Burkholderia mallei (n=26) strains. Within these genes, single nucleotide polymorphisms (SNPs) useful for predicting geographic origin of B. pseudomallei were uncovered. SNPs unique to B. mallei were also identified. Based on 11 SNPs identified in two genes encoding predicted PBP-3s, a dual-locus sequence typing (DLST) scheme was developed. The robustness of this typing scheme was assessed using 1,523 RefSeq genomes from B. pseudomallei (n=1,442) and B. mallei (n=81) strains, resulting in 32 sequence types (STs). Compared to multi-locus sequence typing (MLST), the DLST scheme demonstrated less resolution to support the continental separation of Australian B. pseudomallei strains. However, several STs were unique to strains originating from a specific country or region. The phylogeography of Western Hemisphere B. pseudomallei strains was more highly resolved by DLST compared to internal transcribed spacer (ITS) typing, and all B. mallei strains formed a single ST.

Significance

Conserved genes encoding PBPs in B. pseudomallei are useful for strain typing, can enhance predictions of geographic origin, and differentiate strains of closely related Burkholderia species.

Author Summary

Burkholderia pseudomallei causes the life-threatening disease melioidosis and is considered a biological threat and select agent by the United States government. This soil-dwelling bacterium is commonly found in regions of southeast Asia and northern Australia, but it is also detected in other tropical and sub-tropical areas around the world. With a predicted global burden of 165,000 annual cases and mortality rate that can exceed 40% without prompt and appropriate antibiotic treatment, understanding the epidemiology of melioidosis and mechanisms of antibiotic resistance in B. pseudomallei can benefit public health and safety. Recently, we identified ten conserved genes encoding putative penicillin binding proteins (PBPs) in B. pseudomallei. Here, we examined B. pseudomallei PBP sequences for amino acid mutations that may contribute to β-lactam resistance. We also uncovered nucleotide mutations with utility to predict the geographical origin of B. pseudomallei strains and to differentiate closely related Burkholderia species. Based on 11 informative single nucleotide polymorphisms in two genes each encoding a PBP-3, we developed a simple, targeted dual-locus typing approach.