Comparative genomics confirms a rare melioidosis human-to-human transmission event and reveals incorrect phylogenomic reconstruction due to polyclonality

Genetic diversity, determinants, and dissemination of Burkholderia pseudomallei lineages implicated in melioidosis in Northeast Thailand

Melioidosis is an often-fatal neglected tropical disease caused by an environmental bacterium Burkholderia pseudomallei. However, our understanding of the disease-causing bacterial lineages, their dissemination, and adaptive mechanisms remains limited. To address this, we conduct a comprehensive genomic analysis of 1,391 B. pseudomallei isolates collected from nine hospitals in northeast Thailand between 2015 and 2018, and contemporaneous isolates from neighbouring countries, representing the most densely sampled collection to date. Our study identifies three dominant lineages, each with unique gene sets potentially enhancing bacterial fitness in the environment. We find that recombination drives lineage-specific gene flow. Transcriptome analyses of representative clinical isolates from each dominant lineage reveal increased expression of lineage-specific genes under environmental conditions in two out of three lineages. This underscores the potential importance of environmental persistence for these dominant lineages. The study also highlights the influence of environmental factors such as terrain slope, altitude, and river direction on the geographical dispersal of B. pseudomallei. Collectively, our findings suggest that environmental persistence may play a role in facilitating the spread of B. pseudomallei, and as a prerequisite for exposure and infection, thereby providing useful insights for informing melioidosis prevention and control strategies.

Latent melioidosis activation presenting with urinary tract infection and bacteremia

We report a rare case of latent melioidosis activation in a patient with a distant travel history to an endemic region. Melioidosis is an infection caused by Burkholderia pseudomallei which is highly endemic in Southeast Asia and northern Australia. The patient exhibited common clinical risk factors, presenting with urinary tract infection and bacteremia. The treatment course was complicated by the adverse effect of trimethoprim/sulfamethoxazole. This case underscores the importance of early detection and appropriate treatment of melioidosis, particularly given its expanding global distribution.

Optimal control and cost-effectiveness analysis for the human melioidosis model

In this work, we formulated and investigated an optimal control problem of the melioidosis epidemic to explain the effectiveness of time-dependent control functions in controlling the spread of the epidemic. The basic reproduction number ( R 0 c ) with control measures is obtained, using the next-generation matrix approach and the impact of the controls on R 0 c is illustrated numerically. The optimal control problem is analyzed using Pontryagin's maximum principle to derive the optimality system. The optimality system is simulated using the forward-backward sweep method based on the fourth-order Runge-Kutta method in the MATLAB program to illustrate the impact of all the possible combinations of the control interventions on the transmission dynamics of the disease. The numerical results indicate that among strategies considered, strategy C is shown to be the most effective in reducing the number of infectious classes compared to both strategy A and strategy B. Furthermore, we carried out a cost-effectiveness analysis to determine the most cost-effective strategy and the result indicated that the strategy B (treatment control strategy) should be recommended to mitigate the spread and impact of the disease regarding the costs of the strategies.

Genomic epidemiology reveals geographical clustering of multidrug-resistant Escherichia coli ST131 associated with bacteraemia in Wales

Antibiotic resistance is a significant global public health concern. Uropathogenic Escherichia coli sequence type (ST)131, a widely prevalent multidrug-resistant clone, is frequently associated with bacteraemia. This study investigates third-generation cephalosporin resistance in bloodstream infections caused by E. coli ST131. From 2013-2014 blood culture surveillance in Wales, 142 E. coli ST131 genomes were studied alongside global data. All three major ST131 clades were represented across Wales, with clade C/H30 predominant (n = 102/142, 71.8%). Consistent with global findings, Welsh strains of clade C/H30 contain β-lactamase genes from the blaCTX-M-1 group (n = 65/102, 63.7%), which confer resistance to third-generation cephalosporins. Most Welsh clade C/H30 genomes belonged to sub-clade C2/H30Rx (58.3%). A Wales-specific sub-lineage, named GB-WLS.C2, diverged around 1996-2000. An introduction to North Wales around 2002 led to a localised cluster by 2009, depicting limited genomic diversity within North Wales. This investigation emphasises the value of genomic epidemiology, allowing the detection of genetically similar strains in local areas, enabling targeted and timely public health interventions.

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.

Chronic Splenic Melioidosis in a Patient with Fever of Unknown Origin Diagnosed by Metagenomics Next-Generation Sequencing: An Emerging Cause and Literature Review

Introduction

Human melioidosis is an emerging infectious disease in tropical areas of China, and chronic melioidosis can be a rare cause of fever of unknown origin (FUO). Timely diagnosis may improve the prognosis of melioidosis.

Case Presentation

We report a case of melioidosis with splenic abscesses caused by Burkholderia pseudomallei in a 57-year-old man, who presented with FUO. Positron emission tomography/computed tomography (PET/CT) revealed multiple hypermetabolic lesions in the spleen. The spleen biopsy was conducted and metagenomics next-generation sequencing (mNGS) of the spleen specimen identified the presence of B. pseudomallei, confirming the diagnosis of melioidosis. Antimicrobial treatment was initiated with intravenous meropenem, followed by oral faropenem. During the follow-up, the patient was in good condition except having a low-grade fever occasionally. A splenectomy was performed, and subsequent culture and mNGS of the spleen pus were both positive for B. pseudomallei. Histopathological characteristics of chronic splenic melioidosis were noted.

Conclusion

Melioidosis is a serious endemic disease, and it is critical to raise awareness about this disease.

Hcp1-loaded staphylococcal membrane vesicle vaccine protects against acute melioidosis

Burkholderia pseudomallei is the causal agent of melioidosis, a deadly tropical infectious disease that lacks a vaccine. On the basis of the attenuated Staphylococcus aureus RN4220-Δagr (RN), we engineered the RN4220-Δagr/pdhB-hcp1 strain (RN-Hcp1) to generate B. pseudomallei hemolysin-coregulated protein 1 (Hcp1)-loaded membrane vesicles (^hcp1MVs). The immunization of BALB/c mice with ^hcp1MVs mixed with adjuvant by a three-dose regimen increased the serum specific IgG production. The serum levels of inflammatory factors, including TNF-α and IL-6, in ^hcp1MV-vaccinated mice were comparable with those in PBS-challenged mice. The partial adjuvant effect of staphylococcal MVs was observed with the elevation of specific antibody titer in ^hcp1MV-vaccinated mice relative to those that received the recombinant Hcp1 protein (rHcp1) or MVs derived from RN strain (^ΔagrMVs). The ^hcp1MVs/adjuvant vaccine protected 70% of mice from lethal B. pseudomallei challenge. Immunization with ^hcp1MVs only protected 60% of mice, whereas vaccination with rHcp1 or ^ΔagrMVs conferred no protection. Moreover, mice that received ^hcp1MVs/adjuvant and ^hcp1MVs immunization had low serum TNF-α and IL-6 levels and no inflammatory infiltration in comparison with other groups. In addition, all surviving mice in ^hcp1MVs/adjuvant and ^hcp1MVs groups exhibited no culturable bacteria in their lungs, livers, and spleens five days postinfection. Overall, our data highlighted a new strategy for developing B. pseudomallei vaccine and showed that Hcp1-incorporated staphylococcal MV is a promising candidate for the prevention of acute melioidosis.

Low risk of acquiring melioidosis from the environment in the continental United States

Melioidosis is an underreported human disease of tropical and sub-tropical regions caused by the saprophyte Burkholderia pseudomallei. Although most global melioidosis cases are reported from tropical regions in Southeast Asia and northern Australia, there are multiple occurrences from sub-tropical regions, including the United States (U.S.). Most melioidosis cases reported from the continental U.S. are the result of acquiring the disease during travel to endemic regions or from contaminated imported materials. Only two human melioidosis cases from the continental U.S. have likely acquired B. pseudomallei directly from local environments and these cases lived only ~7 km from each other in rural Texas. In this study, we assessed the risk of acquiring melioidosis from the environment within the continental U.S. by surveying for B. pseudomallei in the environment in Texas where these two human melioidosis cases likely acquired their infections. We sampled the environment near the homes of the two cases and at additional sampling locations in surrounding counties in Texas that were selected based on ecological niche modeling. B. pseudomallei was not detected at the residences of these two cases or in the surrounding region. These negative data are important to demonstrate that B. pseudomallei is rare in the environment in the U.S. even at locations where locally acquired human cases likely have occurred, documenting the low risk of acquiring B. pseudomallei infection from the environment in the continental U.S.

Genotyping of Burkholderia pseudomallei Isolated From Patients in South-Western Coastal Region of India

Burkholderia pseudomallei is the causative agent of melioidosis. Various tools have been used to determine the genetic diversity in B. pseudomallei isolates. In this study, Random Amplified Polymorphic DNA (RAPD)-PCR and flagellin gene (fliC) based PCR-Restriction Fragment Length Polymorphism (RFLP) were used to genotype Indian clinical B. pseudomallei isolates. A total of 89 clinical isolates could be grouped in 6 groups (A through F) by RAPD-PCR analysis. Some of the isolates in various groups had identical banding pattern suggesting them to be epidemiologically related. The RAPD groups also correlated with MLST sequence types suggesting the utility of this easy to do typing method. The PCR- RFLP analysis suggested Type III to be the predominant type which is different from other RFLP types reported from Southeast Asia. In conclusion, the results of this study show that RAPD-PCR, a simple genotyping method, may be used for analyzing the B. pseudomallei isolates and also establish epidemiological relevant relatedness among them. The results of fliC PCR-RFLP further suggest the Indian isolates are different from other Southeast Asian isolates.

Genomic diversity and antimicrobial resistance of Prevotella species isolated from chronic lung disease airways

Cystic fibrosis (CF) and chronic obstructive pulmonary disease (COPD) are characterized by increasingly frequent acute pulmonary exacerbations that reduce life quality and length. Human airways are home to a rich polymicrobial environment, which includes members of the obligately anaerobic genus Prevotella. Despite their commonness, surprisingly little is known about the prevalence, role, genomic diversity and antimicrobial resistance (AMR) potential of Prevotella species and strains in healthy and diseased airways. Here, we used comparative genomics to develop a real-time PCR assay to permit rapid Prevotella species identification and quantification from cultures and clinical specimens. Assay specificity was validated across a panel of Prevotella and non-Prevotella species, followed by PCR screening of CF and COPD respiratory-derived cultures. Next, 35 PCR-positive isolates were subjected to whole-genome sequencing. Of eight identified Prevotella species, P. histicola, P. melaninogenica, P. nanceiensis, P. salivae and P. denticola overlapped between participant cohorts. Phylogenomic analysis revealed considerable interhost but limited intrahost diversity, suggesting patient-specific lineages in the lower airways, probably from oral cavity aspirations. Correlation of phenotypic AMR profiles with AMR genes identified excellent correlation between tetQ presence and decreased doxycycline susceptibility, and ermF presence and decreased azithromycin susceptibility and clindamycin resistance. AMR rates were higher in the CF isolates, reflecting greater antibiotic use in this cohort. All tested Prevotella isolates were tobramycin-resistant, providing a potential selection method to improve Prevotella culture retrieval rates. Our addition of 35 airway-derived Prevotella genomes to public databases will enhance ongoing efforts to unravel the role of this diverse and enigmatic genus in both diseased and healthy lungs.

Genomic epidemiology links Burkholderia pseudomallei from individual human cases to B. pseudomallei from targeted environmental sampling in Northern Australia

Each case of melioidosis results from a single event when a human is infected by the environmental bacterium Burkholderia pseudomallei. Darwin in tropical northern Australia has the highest incidences of melioidosis globally and the Darwin Prospective Melioidosis Study (DPMS) commenced in 1989, documenting all culture confirmed melioidosis cases. From 2000-2019 we sampled DPMS patient's environments for B. pseudomallei when a specific location was considered to have been where infection occurred. With the aim to use genomic epidemiology to understand B. pseudomallei transmission and infecting scenarios. Environmental sampling was performed at 98 DPMS patient sites, where we collected 975 environmental samples (742 soil; 233 water). Genotyping matched the clinical and epidemiologically linked environmental B. pseudomallei for 19 patients (19%), with the environmental isolates cultured from soil (n=11) or water (n=8) sources. B. pseudomallei isolates from patients and their local environments that matched on genotyping were whole genome sequenced (WGS). Of the 19 patients with a clinical-environmental genotype match, 17 pairs clustered on a Darwin core genome single-nucleotide polymorphism (SNP) phylogeny, later confirmed by single ST phylogenies and pairwise comparative genomics. When related back to patient clinical scenarios, the matched clinical and environmental B. pseudomallei pairs informed likely modes of infection: percutaneous inoculation, inhalation, and ingestion. Targeted environmental sampling for B. pseudomallei can inform infecting scenarios for melioidosis and dangerous occupational and recreational activities and identify hot spots of B. pseudomallei presence. However, WGS and careful genomics are required to avoid overcalling the relatedness between clinical and environmental isolates of B. pseudomallei.

Genetic diversity and transmission patterns of Burkholderia pseudomallei on Hainan island, China, revealed by a population genomics analysis

Burkholderia pseudomallei is a Gram-negative soil-dwelling bacillus that causes melioidosis, a frequently fatal infectious disease, in tropical and subtropical regions. Previous studies have identified the overall genetic and evolutionary characteristics of B. pseudomallei on a global scale, including its origin and transmission routes. However, beyond its known hyperendemicity foci in northern Australia and Southeast Asia, the distribution and genetic characteristics of B. pseudomallei in most tropical regions remain poorly understood, including in southern China. Here, we sequenced the genomes of 122 B. pseudomallei strains collected from Hainan, an island in southern China, in 2002–2018, to investigate the population structure, relationships with global strains, local epidemiology, and virulence and antimicrobial-resistance factors. A phylogenetic analysis and hierarchical clustering divided the Hainan strains into nine phylogenic groups (PGs), 80 % of which were concentrated within five major groups (group 1: corresponding to minor sequence types [STs], 12.3 %; group 3: ST46 and ST50, 31.1 %; group 9: ST58, 13.1 %; group 11: ST55, 8.2 %; group 15: mainly ST658, 15.6%). A phylogenetic analysis that included global strains suggested that B. pseudomallei in Hainan originated from Southeast Asian countries, transmitted in multiple historical importation events. We also identified several mutual transmission events between Hainan and Southeast Asian countries in recent years, including three importation events from Thailand and Singapore to Hainan and three exportation events from Hainan to Singapore, Malaysia, and Taiwan island. A statistical analysis of the temporal distribution showed that the Hainan strains of groups 3, 9, and 15 have dominated the disease epidemic locally in the last 5 years. The spatial distribution of the Hainan strains demonstrated that some PGs are distributed in different cities on Hainan island, and by combining phylogenic and geographic distribution information, we detected 21 between-city transmission events, indicating its frequent local transmission. The detection of virulence factor genes showed that 56 % of the Hainan strains in group 1 encode a B. pseudomallei-specific adherence factor, boaB, confirming the specific pathogenic characteristics of the Hainan strains in group 1. An analysis of the antimicrobial-resistance potential of B. pseudomallei showed that various kinds of alterations were identified in clinically relevant antibiotic resistance factors, such as AmrR, PenA and PBP3, etc. Our results clarify the population structure, local epidemiology, and pathogenic characteristics of B. pseudomallei in Hainan, providing further insight into its regional and global transmission networks and improving our knowledge of its global phylogeography.

Development and Validation of a Burkholderia pseudomallei Core Genome Multilocus Sequence Typing Scheme To Facilitate Molecular Surveillance

Burkholderia pseudomallei causes the severe disease melioidosis. Whole-genome sequencing (WGS)-based typing methods currently offer the highest resolution for molecular investigations of this genetically diverse pathogen. Still, its routine application in diagnostic laboratories is limited by the need for high computing power, bioinformatic skills, and variable bioinformatic approaches, with the latter affecting the results. We therefore aimed to establish and validate a WGS-based core genome multilocus sequence typing (cgMLST) scheme, applicable in routine diagnostic settings. A soft defined core genome was obtained by challenging the B. pseudomallei reference genome K96243 with 469 environmental and clinical genomes, resulting in 4,221 core and 1,351 accessory targets. The scheme was validated with 320 WGS data sets. We compared our novel typing scheme with single nucleotide polymorphism-based approaches investigating closely and distantly related strains. Finally, we applied our scheme for tracking the environmental source of a recent infection. The validation of the scheme detected >95% good cgMLST target genes in 98.4% of the genomes. Comparison with existing typing methods revealed very good concordance. Our scheme proved to be applicable to investigating not only closely related strains but also the global B. pseudomallei population structure. We successfully utilized our scheme to identify a sugarcane field as the presumable source of a recent melioidosis case. In summary, we developed a robust cgMLST scheme that integrates high resolution, maximized standardization, and fast analysis for the nonbioinformatician. Our typing scheme has the potential to serve as a routinely applicable classification system in B. pseudomallei molecular epidemiology.

Molecular Epidemiology of Third-Generation-Cephalosporin-Resistant Enterobacteriaceae in Southeast Queensland, Australia

Third-generation cephalosporin-resistant (3GC-R) Enterobacteriaceae represent a major threat to human health. Here, we captured 288 3GC-R Enterobacteriaceae clinical isolates from 264 patients presenting at a regional Australian hospital over a 14-month period. In addition to routine mass spectrometry and antibiotic sensitivity testing, isolates were examined using rapid (∼40-min) real-time PCR assays targeting the most common extended-spectrum β-lactamases (ESBLs; bla _CTX-M-1 and bla _CTX-M-9 groups, plus bla _TEM, bla _SHV, and an internal 16S rRNA gene control). AmpC CMY β-lactamase (bla _CMY) prevalence was also examined. Escherichia coli (80.2%) and Klebsiella pneumoniae (17.0%) were dominant, with Klebsiella oxytoca, Klebsiella aerogenes, and Enterobacter cloacae infrequently identified. Ceftriaxone and cefoxitin resistance were identified in 97.0% and 24.5% of E. coli and K. pneumoniae isolates, respectively. Consistent with global findings in Enterobacteriaceae, most (98.3%) isolates harbored at least one β-lactamase gene, with 144 (50%) harboring bla _CTX-M-1 group, 92 (31.9%) harboring bla _CTX-M-9 group, 48 (16.7%) harboring bla _SHV, 133 (46.2%) harboring bla _TEM, and 34 (11.8%) harboring bla _CMY genes. A subset of isolates (n = 98) were subjected to whole-genome sequencing (WGS) to identify the presence of cryptic resistance determinants and to verify genotyping accuracy. WGS of β-lactamase-negative or carbapenem-resistant isolates identified uncommon ESBL and carbapenemase genes, including bla _NDM and bla _IMP, and confirmed all PCR-positive genotypes. We demonstrate that our PCR assays enable the rapid and cost-effective identification of ESBLs in the hospital setting, which has important infection control and therapeutic implications.

An Evolutionary Arms Race Between Burkholderia pseudomallei and Host Immune System: What Do We Know?

A better understanding of co-evolution between pathogens and hosts holds promise for better prevention and control strategies. This review will explore the interactions between Burkholderia pseudomallei, an environmental and opportunistic pathogen, and the human host immune system. B. pseudomallei causes "Melioidosis," a rapidly fatal tropical infectious disease predicted to affect 165,000 cases annually worldwide, of which 89,000 are fatal. Genetic heterogeneities were reported in both B. pseudomallei and human host population, some of which may, at least in part, contribute to inter-individual differences in disease susceptibility. Here, we review (i) a multi-host-pathogen characteristic of the interaction; (ii) selection pressures acting on B. pseudomallei and human genomes with the former being driven by bacterial adaptation across ranges of ecological niches while the latter are driven by human encounter of broad ranges of pathogens; (iii) the mechanisms that generate genetic diversity in bacterial and host population particularly in sequences encoding proteins functioning in host-pathogen interaction; (iv) reported genetic and structural variations of proteins or molecules observed in B. pseudomallei-human host interactions and their implications in infection outcomes. Together, these predict bacterial and host evolutionary trajectory which continues to generate genetic diversity in bacterium and operates host immune selection at the molecular level.

Pharmacodynamic Evaluation of Plasma and Epithelial Lining Fluid Exposures of Amikacin against Pseudomonas aeruginosa in a Dynamic In Vitro Hollow-Fiber Infection Model

Given that aminoglycosides, such as amikacin, may be used for multidrug-resistant Pseudomonas aeruginosa infections, optimization of therapy is paramount for improved treatment outcomes. This study aims to investigate the pharmacodynamics of different simulated intravenous amikacin doses on susceptible P. aeruginosa to inform ventilator-associated pneumonia (VAP) and sepsis treatment choices. A hollow-fiber infection model with two P. aeruginosa isolates (MICs of 2 and 8 mg/liter) with an initial inoculum of ∼108 CFU/ml was used to test different amikacin dosing regimens. Three regimens (15, 25, and 50 mg/kg) were tested to simulate a blood exposure, while a 30 mg/kg regimen simulated the epithelial lining fluid (ELF) for potential respiratory tract infection. Data were described using a semimechanistic pharmacokinetic/pharmacodynamic (PK/PD) model. Whole-genome sequencing was used to identify mutations associated with resistance emergence. While bacterial density was reduced by >6 logs within the first 12 h in simulated blood exposures following this initial bacterial kill, there was amplification of a resistant subpopulation with ribosomal mutations that were likely mediating amikacin resistance. No appreciable bacterial killing occurred with subsequent doses. There was less (<5 log) bacterial killing in the simulated ELF exposure for either isolate tested. Simulation studies suggested that a dose of 30 and 50 mg/kg may provide maximal bacterial killing for bloodstream and VAP infections, respectively. Our results suggest that amikacin efficacy may be improved with the use of high-dose therapy to rapidly eliminate susceptible bacteria. Subsequent doses may have reduced efficacy given the rapid amplification of less-susceptible bacterial subpopulations with amikacin monotherapy.