Molecular phylogeny of Burkholderia pseudomallei from a remote region of Papua New Guinea

A systematic review and meta-analysis of the global prevalence and relationships among Burkholderia pseudomallei sequence types isolated from humans, animals, and the environment

Background and Aim

Burkholderia pseudomallei, a highly pathogenic bacterium responsible for melioidosis, exhibits ecological ubiquity and thrives within soil and water reservoirs, posing significant infection risks to humans and animals through direct contact. The aim of this study was to elucidate the genetic diversity and prevalence patterns of B. pseudomallei sequence types (STs) across a global spectrum and to understand the relationships between strains isolated from different sources.

Materials and Methods

We performed a systematic review and meta-analysis in this study. Extensive research was carried out across three comprehensive databases, including PubMed, Scopus, and ScienceDirect with data collected from 1924 to 2023.

Results

A total of 40 carefully selected articles contributed 2737 B. pseudomallei isolates attributed to 729 distinct STs and were incorporated into the systematic review. Among these, ST46 emerged as the most prominent, featuring in 35% of the articles and demonstrating a dominant prevalence, particularly within Southeast Asia. Moreover, ST51 consistently appeared across human, animal, and environmental studies. Subsequently, we performed a meta-analysis, focusing on nine specific STs: ST46, ST51, ST54, ST70, ST84, ST109, ST289, ST325, and ST376. Surprisingly, no statistically significant differences in their pooled prevalence proportions were observed across these compartments for ST46, ST70, ST289, ST325, and ST376 (all p > 0.69). Conversely, the remaining STs, including ST51, ST54, ST84, and ST109, displayed notable variations in their prevalence among the three domains (all p < 0.04). Notably, the pooled prevalence of ST51 in animals and environmental samples surpassed that found in human isolates (p < 0.01).

Conclusion

To the best of our knowledge, this study is the first systematic review and meta-analysis to investigate the intricate relationships between STs and their sources and contributes significantly to our understanding of B. pseudomallei diversity within the One Health framework.

Phylogeographic characterization of Burkholderia pseudomallei isolated from Bangladesh

BACKGROUND

Burkholderia pseudomallei possesses a diverse set of genes which encode a vast array of biological functions reflecting its clinical, ecological and phenotypic diversity. Strain variation is linked to geographic location as well as pattern of land uses. This soil-dwelling Gram-negative pathogen causes melioidosis, a tropical disease endemic in northern Australia and Southeast Asian regions including Bangladesh. Phylogeographic analyses of B. pseudomallei isolates by molecular typing techniques could be used to examine the diversity of this organism as well as to track melioidosis epidemics.

METHODS

In this study, 22 B. pseudomallei isolates, of which 20 clinical and two soil isolates were analyzed, utilizing Real-time PCR assay and multilocus sequence typing (MLST). The sequences were then submitted to PubMLST database for analysis and construction of phylogenetic tree.

FINDINGS

A total of 12 different sequence types (STs) that includes four novel STs were identified for the first time. Strains having STs 1005, 1007 and 56 were the most widespread STs frequently isolated in Bangladesh. ST 1005, ST 56, ST 1007 and ST 211 have been detected not only in Bangladesh but are also present in many Southeast Asian countries.

SIGNIFICANCE

ST 1005 was detected in both soil and clinical samples of Gazipur. Most prevalent, ST 56 has been previously reported from Myanmar, Thailand, Cambodia and Vietnam, confirming the persistence of the genotype over the entire continent. Further large-scale study is necessary to find out the magnitude of the infection and its different reservoirs in the environment along with phylogeographic association.

Antimicrobial Resistance in Papua New Guinea: A Narrative Scoping Review

Antimicrobial-resistant bacterial infections are a known threat to the public health of low-income countries and are undercharacterized in Papua New Guinea. A scoping literature review of scientific peer-reviewed publications on antimicrobial resistance in Papua New Guinea was conducted, and their results were summarized. Many of the available data on resistant bacteria in Papua New Guinea have come from Port Moresby and Goroka and have been focused on Staphylococcus aureus, as well as important pediatric pathogens such as Streptococcus pneumoniae and Haemophilus influenzae. Progressive resistance to the commonly used antibiotics penicillin and chloramphenicol among most clinically important bacterial pathogens has prompted healthcare workers to adopt expensive broad-spectrum antibiotics. There is already evidence of resistance to newly adopted antibiotics among several Gram-negative organisms. Drivers of antimicrobial resistance in Papua New Guinea include a high burden of infectious diseases, inappropriate antibiotic prescription practices, poor regulation of antibiotics, incomplete adherence, substandard drug quality, and overcrowding of healthcare facilities. There is a lack of information on antimicrobial resistance among priority pathogens and from several important regions of Papua New Guinea.

In silico analyses of penicillin binding proteins in Burkholderia pseudomallei uncovers SNPs with utility for phylogeography, species differentiation, and sequence typing

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. 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. 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.

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.

Epidemiological and Clinical Characteristics of Melioidosis Caused by Gentamicin-Susceptible Burkholderia pseudomallei in Sarawak, Malaysia

Background

Burkholderia pseudomallei, the causative agent of melioidosis, is intrinsically resistant to a broad range of antibiotics, including aminoglycosides. In Sarawak, Malaysia, a high proportion of melioidosis cases are caused by gentamicin-susceptible isolates. There are limited epidemiological and clinical data on these infections.

Methods

We conducted a retrospective study of culture-confirmed melioidosis among adults admitted to Bintulu Hospital in Sarawak, Malaysia, from January 2011 until December 2016.

Results

One hundred forty-eight adults with culture-confirmed melioidosis were identified. Of 129 (87%) tested, 84 (65%) had gentamicin-susceptible B pseudomallei. The average annual incidence of melioidosis was 12.3 per 100 000 population, with marked variation between districts ranging from 5.8 to 29.3 per 100 000 population. Rural districts had higher incidences of melioidosis and overwhelmingly larger proportions of gentamicin-susceptible B pseudomallei infection. Significantly more patients with gentamicin-susceptible infection had no identified risk factors, with diabetes less frequently present in this group. Ninety-eight percent had acute presentations. Pneumonia, reported in 71%, was the most common presentation. Splenic abscesses were found in 54% of those imaged. Bacteremia was present in 88%; septic shock occurred in 47%. Forty-five (35%) patients died. No differences in clinical, laboratory, or outcome characteristics were noted between gentamicin-susceptible and gentamicin-resistant infections.

Conclusions

Gentamicin-susceptible B pseudomallei infections are common in Sarawak and dominate in the high-incidence rural interior regions. Clinical manifestations and outcomes are the same as for gentamicin-resistant B pseudomallei infections. Further studies are required to determine if all gentamicin-susceptible B pseudomallei infections in Sarawak are clonal and to ascertain their environmental drivers and niches.

Selective whole genome amplification and sequencing of Coxiella burnetii directly from environmental samples

Whole genome sequencing (WGS) is a widely available, inexpensive means of providing a wealth of information about an organism's diversity and evolution. However, WGS for many pathogenic bacteria remain limited because they are difficult, slow and/or dangerous to culture. To avoid culturing, metagenomic sequencing can be performed directly on samples, but the sequencing effort required to characterize low frequency organisms can be expensive. Recently developed methods for selective whole genome amplification (SWGA) can enrich target DNA to provide efficient sequencing. We amplified Coxiella burnetii (a bacterial select agent and human/livestock pathogen) from 3 three environmental samples that were overwhelmed with host DNA. The 68- to 147-fold enrichment of the bacterial sequences provided enough genome coverage for SNP analyses and phylogenetic placement. SWGA is a valuable tool for the study of difficult-to-culture organisms and has the potential to facilitate high-throughput population characterizations as well as targeted epidemiological or forensic investigations.

Burkholderia pseudomallei distribution in Australasia is linked to paleogeographic and anthropogenic history

Burkholderia pseudomallei is the environmental bacillus that causes melioidosis; a disease clinically significant in Australia and Southeast Asia but emerging in tropical and sub-tropical regions around the globe. Previous studies have placed the ancestral population of the organism in Australia with a single lineage disseminated to Southeast Asia. We have previously characterized B. pseudomallei isolates from New Guinea and the Torres Strait archipelago; remote regions that share paleogeographic ties with Australia. These studies identified regional biogeographical boundaries. In this study, we utilize whole-genome sequencing to reconstruct ancient evolutionary relationships and ascertain correlations between paleogeography and present-day distribution of this bacterium in Australasia. Our results indicate that B. pseudomallei from New Guinea fall into a single clade within the Australian population. Furthermore, clades from New Guinea are region-specific; an observation possibly linked to limited recent anthropogenic influence in comparison to mainland Australia and Southeast Asia. Isolates from the Torres Strait archipelago were distinct yet scattered among those from mainland Australia. These results provide evidence that the New Guinean and Torres Strait lineages may be remnants of an ancient portion of the Australian population. Rising sea levels isolated New Guinea and the Torres Strait Islands from each other and the Australian mainland, and may have allowed long-term isolated evolution of these lineages, providing support for a theory of microbial biogeography congruent with that of macro flora and fauna. Moreover, these findings indicate that contemporary microbial biogeography theories should consider recent and ongoing impacts of globalisation and human activity.

Phylogenetic Diversity of Burkholderia pseudomallei isolated from veterinary cases and the environments in Peninsular Malaysia

This study was designed to determine the genotype and the phylogeny of Burkholderia pseudomallei isolated from veterinary cases and from the animal environments in Peninsular Malaysia. The Malaysian B. pseudomallei population were then compared to those found elsewhere. A total of 113 isolates from veterinary cases (35) and the environment (56 from soil and 22 from water) were characterized using multilocus sequence typing (MLST). Two novel alleles, allele 97 and 69 of the gene locus ace and lepA respectively were recovered. Isolates were resolved into 12 distinct sequence types (STs) out of which five were novel, namely ST1130, ST1131, ST1338, ST1339 and ST1367. The isolates from veterinary cases co-clustered with those from the environment. B. pseudomallei isolates in this study were highly clonal and have descended from a common ancestor clonal complex (CC) 48 found in Southeast Asia. This study shows that veterinary case isolates are often caused by similar STs, with similar populations found in the direct animal environment and those previously reported to cause human infections in Malaysia and elsewhere. Isolates of B. pseudomallei from human infections have been given more attention, with a comparatively lower focus on isolates from animals and the farm environment. This study highlighted the genotype and phylogeny of B. pseudomallei isolated from animals and the environment and their relations to the isolates from human cases reported in Malaysia and elsewhere. Most STs reported in this study, from veterinary cases and animal environment are similar to those previously reported as causing human infections in Malaysia and elsewhere. Therefore, even though direct zoonosis is uncommon, monitoring melioidosis occurrences in animals can provide insights on the bacterial strains infecting humans.

Melioidosis in Papua New Guinea and Oceania

Melioidosis has only been sporadically reported throughout Melanesia and the Pacific region since the first report from Guam in 1946; therefore, its contribution to the disease burden in this region is largely unknown. However, the outcome of a small number of active surveillance programs, serological surveys, and presumptive imported cases identified elsewhere provide an insight into its epidemiology and potential significance throughout the region. Both clinical cases and environmental reservoirs have been described from the rural district of Balimo in the Western Province of Papua New Guinea and from the Northern Province of New Caledonia. In both these locations the incidence of disease is similar to that described in tropical Australia and Burkholderia pseudomallei isolates are also phylogenetically linked to Australian isolates. Serological evidence and presumptive imported cases identified elsewhere suggest that melioidosis exists in other countries throughout the Pacific. However, the lack of laboratory facilities and clinical awareness, and the burden of other infections of public health importance such as tuberculosis, contribute to the under-recognition of melioidosis in this region.

Phylogeographic, genomic, and meropenem susceptibility analysis of Burkholderia ubonensis

The bacterium Burkholderia ubonensis is commonly co-isolated from environmental specimens harbouring the melioidosis pathogen, Burkholderia pseudomallei. B. ubonensis has been reported in northern Australia and Thailand but not North America, suggesting similar geographic distribution to B. pseudomallei. Unlike most other Burkholderia cepacia complex (Bcc) species, B. ubonensis is considered non-pathogenic, although its virulence potential has not been tested. Antibiotic resistance in B. ubonensis, particularly towards drugs used to treat the most severe B. pseudomallei infections, has also been poorly characterised. This study examined the population biology of B. ubonensis, and includes the first reported isolates from the Caribbean. Phylogenomic analysis of 264 B. ubonensis genomes identified distinct clades that corresponded with geographic origin, similar to B. pseudomallei. A small proportion (4%) of strains lacked the 920kb chromosome III replicon, with discordance of presence/absence amongst genetically highly related strains, demonstrating that the third chromosome of B. ubonensis, like other Bcc species, probably encodes for a nonessential pC3 megaplasmid. Multilocus sequence typing using the B. pseudomallei scheme revealed that one-third of strains lack the "housekeeping" narK locus. In comparison, all strains could be genotyped using the Bcc scheme. Several strains possessed high-level meropenem resistance (≥32 μg/mL), a concern due to potential transmission of this phenotype to B. pseudomallei. In silico analysis uncovered a high degree of heterogeneity among the lipopolysaccharide O-antigen cluster loci, with at least 35 different variants identified. Finally, we show that Asian B. ubonensis isolate RF23-BP41 is avirulent in the BALB/c mouse model via a subcutaneous route of infection. Our results provide several new insights into the biology of this understudied species.

Molecular evidence of Burkholderia pseudomallei genotypes based on geographical distribution

Background. Central intermediary metabolism (CIM) in bacteria is defined as a set of metabolic biochemical reactions within a cell, which is essential for the cell to survive in response to environmental perturbations. The genes associated with CIM are commonly found in both pathogenic and non-pathogenic strains. As these genes are involved in vital metabolic processes of bacteria, we explored the efficiency of the genes in genotypic characterization of Burkholderia pseudomallei isolates, compared with the established pulsed-field gel electrophoresis (PFGE) and multilocus sequence typing (MLST) schemes.

Methods. Nine previously sequenced B. pseudomallei isolates from Malaysia were characterized by PFGE, MLST and CIM genes. The isolates were later compared to the other 39 B. pseudomallei strains, retrieved from GenBank using both MLST and sequence analysis of CIM genes. UniFrac and hierachical clustering analyses were performed using the results generated by both MLST and sequence analysis of CIM genes.

Results. Genetic relatedness of nine Malaysian B. pseudomallei isolates and the other 39 strains was investigated. The nine Malaysian isolates were subtyped into six PFGE profiles, four MLST profiles and five sequence types based on CIM genes alignment. All methods demonstrated the clonality of OB and CB as well as CMS and THE. However, PFGE showed less than 70% similarity between a pair of morphology variants, OS and OB. In contrast, OS was identical to the soil isolate, MARAN. To have a better understanding of the genetic diversity of B. pseudomallei worldwide, we further aligned the sequences of genes used in MLST and genes associated with CIM for the nine Malaysian isolates and 39 B. pseudomallei strains from NCBI database. Overall, based on the CIM genes, the strains were subtyped into 33 profiles where majority of the strains from Asian countries were clustered together. On the other hand, MLST resolved the isolates into 31 profiles which formed three clusters. Hierarchical clustering using UniFrac distance suggested that the isolates from Australia were genetically distinct from the Asian isolates. Nevertheless, statistical significant differences were detected between isolates from Malaysia, Thailand and Australia.

Discussion. Overall, PFGE showed higher discriminative power in clustering the nine Malaysian B. pseudomallei isolates and indicated its suitability for localized epidemiological study. Compared to MLST, CIM genes showed higher resolution in distinguishing those non-related strains and better clustering of strains from different geographical regions. A closer genetic relatedness of Malaysian isolates with all Asian strains in comparison to Australian strains was observed. This finding was supported by UniFrac analysis which resulted in geographical segregation between Australia and the Asian countries.

Multilocus sequence typing of 102 Burkholderia pseudomallei strains isolated from China

The phylogenetic and epidemiological relationships of 102 Burkholderia pseudomallei clinical isolates from different geographical and population sources in China were investigated by multilocus sequence typing (MLST). The MLST data were analysed using the e-BURST algorithm, and an unweighted pair-group method with arithmetic mean dendrogram was constructed based on the pair-wise differences in the allelic profiles of the strains. Forty-one sequence types (STs) were identified, of which eight were novel (ST1341, ST1345, ST1346, ST1347, ST1348, ST1349, ST1350, ST1351). No geographical-specific or host population-specific phylogenetic lineages were identified. ST46, ST50, ST55, ST58, ST70 and ST1095 predominated, but ~44% of isolates were assigned to 45 STs illustrating high genetic diversity in the strain collection. Additionally, the phylogenetic relationships of the dominant STs in China showed significant linkeage with B. pseudomallei isolates from Thailand. Analysis of the gmhD allele suggests high genetic variation in B. pseudomallei in China.

Melioidosis in New Caledonia: a dominant strain in a transmission hotspot

Melioidosis is an infectious disease caused by Burkholderia pseudomallei, a bacterium endemic in Southeast Asia and northern Australia. In New Caledonia, sporadic cases were first described in 2005; since then, more cases have been identified. To improve our understanding of melioidosis epidemiology in New Caledonia, we compared the local cases and B. pseudomallei isolates with those from endemic areas. Nineteen melioidosis cases have been diagnosed in New Caledonia since 1999, mostly severe and with frequent bacteraemia, leading to three (16%) fatalities. All but one occurred in the North Province. Besides sporadic cases caused by non-clonal strains, we also identified a hotspot of transmission related to a clonal group of B. pseudomallei that is phylogenetically related to Australian strains.

Environmental Attributes Influencing the Distribution of Burkholderia pseudomallei in Northern Australia

Factors responsible for the spatial and temporal clustering of Burkholderia pseudomallei in the environment remain to be elucidated. Whilst laboratory based experiments have been performed to analyse survival of the organism in various soil types, such approaches are strongly influenced by alterations to the soil micro ecology during soil sanitisation and translocation. During the monsoonal season in Townsville, Australia, B. pseudomallei is discharged from Castle Hill (an area with a very high soil prevalence of the organism) by groundwater seeps and is washed through a nearby area where intensive sampling in the dry season has been unable to detect the organism. We undertook environmental sampling and soil and plant characterisation in both areas to ascertain physiochemical and macro-floral differences between the two sites that may affect the prevalence of B. pseudomallei. In contrast to previous studies, the presence of B. pseudomallei was correlated with a low gravimetric water content and low nutrient availability (nitrogen and sulphur) and higher exchangeable potassium in soils favouring recovery. Relatively low levels of copper, iron and zinc favoured survival. The prevalence of the organism was found to be highest under the grasses Aristida sp. and Heteropogon contortus and to a lesser extent under Melinis repens. The findings of this study indicate that a greater variety of factors influence the endemicity of melioidosis than has previously been reported, and suggest that biogeographical boundaries to the organisms' distribution involve complex interactions.

Burkholderia pseudomallei isolates from Sarawak, Malaysian Borneo, are predominantly susceptible to aminoglycosides and macrolides

Melioidosis is a potentially fatal disease caused by the saprophytic bacterium Burkholderia pseudomallei. Resistance to gentamicin is generally a hallmark of B. pseudomallei, and gentamicin is a selective agent in media used for diagnosis of melioidosis. In this study, we determined the prevalence and mechanism of gentamicin susceptibility found in B. pseudomallei isolates from Sarawak, Malaysian Borneo. We performed multilocus sequence typing and antibiotic susceptibility testing on 44 B. pseudomallei clinical isolates from melioidosis patients in Sarawak district hospitals. Whole-genome sequencing was used to identify the mechanism of gentamicin susceptibility. A novel allelic-specific PCR was designed to differentiate gentamicin-sensitive isolates from wild-type B. pseudomallei. A reversion assay was performed to confirm the involvement of this mechanism in gentamicin susceptibility. A substantial proportion (86%) of B. pseudomallei clinical isolates in Sarawak, Malaysian Borneo, were found to be susceptible to the aminoglycoside gentamicin, a rare occurrence in other regions where B. pseudomallei is endemic. Gentamicin sensitivity was restricted to genetically related strains belonging to sequence type 881 or its single-locus variant, sequence type 997. Whole-genome sequencing identified a novel nonsynonymous mutation within amrB, encoding an essential component of the AmrAB-OprA multidrug efflux pump. We confirmed the role of this mutation in conferring aminoglycoside and macrolide sensitivity by reversion of this mutation to the wild-type sequence. Our study demonstrates that alternative B. pseudomallei selective media without gentamicin are needed for accurate melioidosis laboratory diagnosis in Sarawak. This finding may also have implications for environmental sampling of other locations to test for B. pseudomallei endemicity.

Biogeography of Burkholderia pseudomallei in the Torres Strait Islands of Northern Australia

It has been hypothesized that biogeographical boundaries are a feature of Burkholderia pseudomallei ecology, and they impact the epidemiology of melioidosis on a global scale. This study examined the relatedness of B. pseudomallei sourced from islands in the Torres Strait of Northern Australia to determine if the geography of isolated island communities is a determinant of the organisms' dispersal. Environmental sampling on Badu Island in the Near Western Island cluster recovered a single clone. An additional 32 clinical isolates from the region were sourced. Isolates were characterized using multilocus sequence typing and a multiplex PCR targeting the flagellum gene cluster. Gene cluster analysis determined that 69% of the isolates from the region encoded the ancestral Burkholderia thailandensis-like flagellum and chemotaxis gene cluster, a proportion significantly lower than that reported from mainland Australia and consistent with observations of isolates from southern Papua New Guinea. A goodness-of-fit test indicated that there was geographic localization of sequence types throughout the archipelago, with the exception of Thursday Island, the economic and cultural hub of the region. Sequence types common to mainland Australia and Papua New Guinea were identified. These findings demonstrate for the first time an environmental reservoir for B. pseudomallei in the Torres Strait, and multilocus sequence typing suggests that the organism is not randomly distributed throughout this region and that seawater may provide a barrier to dispersal of the organism. Moreover, these findings support an anthropogenic dispersal hypothesis for the spread of B. pseudomallei throughout this region.