Search

BackgroundAntimicrobial resistance (AMR) is an intensifying threat that requires urgent mitigation to avoid a post-antibiotic era. The bacterial pathogen, Pseudomonas aeruginosa, represents one of the greatest AMR concerns due to ever-increasing multi- and pan-drug resistance rates. Shotgun sequencing is gaining traction for AMR profiling due to its unambiguity and data transferability; however, existing tools for in silico AMR prediction perform poorly for this pathogen.MethodsWe curated the most comprehensive database to date of known P. aeruginosa AMR determinants. We then performed comparative genomics and microbial genome-wide association (mGWAS) across a Global isolate Dataset (n=1391) with paired antimicrobial phenotype and genome data to identify novel AMR determinants. Finally, we assessed the predictive performance of our ARDaP-implemented database, along with AMRFinderPlus and ResFinder, for both the Global Dataset and an additional Validation Dataset of 89 isolates.FindingsOur AMR database of 754 chromosomal variants includes 79 newly discovered variants associated with AMR phenotype/s, and 246 variants in AMR-encoding loci that do not confer AMR. Interrogation of the Global Dataset with our AMR database achieved genotype-phenotype concordance (as measured by balanced accuracy [bACC]) of 86%, compared with bACCs of 57% with AMRFinderPlus and 63% with ResFinder. For the Validation Dataset, bACCs of 80%, 54%, and 53%, respectively, were attained.InterpretationOur pipeline provides the most accurate tool yet for AMR prediction in P. aeruginosa. Implementing our ARDaP-compatible AMR determinant database for routine AMR prediction from P. aeruginosa (meta)genomes will vastly improve AMR identification and treatment modalities, addressing a critical facet in combatting this treatment-refractory pathogen.FundingThis work was supported by Advance Queensland (AQIRF0362018 and AQRF13016-17RD2), the Wishlist Sunshine Coast Health Foundation (2019-14), and the National Health and Medical Research Council (GNT455919). The funders had no role in study design, data acquisition, analysis, interpretation, writing, or manuscript submission.Research in contextEvidence before the studySeveral bioinformatic tools have been developed with the purpose of facilitating rapid and accurate AMR prediction from microbial (meta)genomic data, particularly for superbugs. To this end, we searched PubMed on July 29, 2022 using ((((“pseudomonas aeruginosa”) AND (“antimicrobial resistance”)) AND (“genome”)) AND (“predict*”)) NOT (review) for studies published in all languages. Twenty-four results were obtained, of which 20 described the use of one or more bioinformatic methods to predict AMR determinants (8 x ResFinder, 6 x CARD, 4 x comparative genomics, 4 x machine learning, 2 x transcriptomics, 1 x OpGen internal bioinformatic pipeline). Only 1 study created an AMR determinant database, which demonstrated limited ability to predict AMR across a dataset of 100 P. aeruginosa isolates. No study has yet created a publicly available, open-source tool and associated AMR determinant database for widespread use.Added value of this studyOur ARDaP-compatible AMR determinant database, which was built from the largest-ever paired P. aeruginosa AMR phenotype-genome dataset to date (n=1480 isolates), provides a crucial advancement in the way we use WGS data to rapidly and accurately predict AMR in P. aeruginosa. Using ARDaP, systematic analysis of large volumes of WGS data can be used to develop a more detailed understanding of the specific genetic variants that underpin AMR.Implications of all the available evidenceAMR prediction tools will play an essential role in the battle against superbugs as they provide precise information on the genetic basis of a pathogen’s resistance profile – information that can be used to inform timely antimicrobial-treatment shifts, reduce the reliance on broad-spectrum antibiotics, and ultimately improve patient outcomes and slow the spread of AMR. Our study demonstrates that predicting AMR phenotype from genotype in P. aeruginosa can rival the performance of phenotypic AMR sensitivity testing such as the VITEK® 2, which is routinely used in clinical microbiology laboratories worldwide. The rapid detection of emerging AMR and multidrug resistance, and subsequent implementation of personalised treatment regimens, serves as our best chance to extend antibiotic efficacy amidst the AMR crisis.

Between 2010 and 2015, nocardiosis outbreaks caused by

Between 2010 and 2015, nocardiosis outbreaks caused by Nocardia seriolae affected many permit farms throughout Vietnam, causing mass fish mortalities. To understand the biology, origin and epidemiology of these outbreaks, 20 N . seriolae strains collected from farms in four provinces in the South Central Coast region of Vietnam, along with two Taiwanese strains, were analysed using genetics and genomics. PFGE identified a single cluster amongst all Vietnamese strains that was distinct from the Taiwanese strains. Like the PFGE findings, phylogenomic and SNP genotyping analyses revealed that all Vietnamese N. seriolae strains belonged to a single, unique clade. Strains fell into two subclades that differed by 103 SNPs, with almost no diversity within clades (0–5 SNPs). There was no association between geographical origin and subclade placement, suggesting frequent N. seriolae transmission between Vietnamese mariculture facilities during the outbreaks. The Vietnamese strains shared a common ancestor with strains from Japan and China, with the closest strain, UTF1 from Japan, differing by just 220 SNPs from the Vietnamese ancestral node. Draft Vietnamese genomes range from 7.55 to 7.96 Mbp in size, have an average G+C content of 68.2 % and encode 7 602–7958 predicted genes. Several putative virulence factors were identified, including genes associated with host cell adhesion, invasion, intracellular survival, antibiotic and toxic compound resistance, and haemolysin biosynthesis. Our findings provide important new insights into the epidemiology and pathogenicity of N. seriolae and will aid future vaccine development and disease management strategies, with the ultimate goal of nocardiosis-free aquaculture.

Since 2005, the range of Burkholderia pseudomallei sequence type 562 (ST562) has expanded in northern Australia. During 2005–2019, ST562 caused melioidosis in 61 humans and 3 animals. Cases initially occurred in suburbs surrounding a creek before spreading across urban Darwin, Australia and a nearby island community. In urban Darwin, ST562 caused 12% (53/440) of melioidosis cases, a proportion that increased during the study period. We analyzed 2 clusters of cases with epidemiologic links and used genomic analysis to identify previously unassociated cases. We found that ST562 isolates from Hainan Province, China, and Pingtung County, Taiwan, were distantly related to ST562 strains from Australia. Temporal genomic analysis suggested a single ST562 introduction into the Darwin region in ≈1988. The origin and transmission mode of ST562 into Australia remain uncertain.

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

Nocardia seriolae has caused significant fish losses in Asia and the Americas in recent decades, including in Vietnam, which has witnessed devastating economic and social impacts due to this bacterial pathogen. Surveillance strategies are urgently needed to mitigate N. seriolae dissemination in Vietnamese aquaculture and mariculture industries. Whole-genome sequencing (WGS) offers the highest level of resolution to discriminate closely related strains and to determine their putative origin and transmission routes. However, WGS is impractical for epidemiological investigations and pathogen surveillance due to its time-consuming and costly nature, putting this technology out-of-reach for many industry end-users. To overcome this issue, we targeted two previously characterised, phylogenetically informative single-nucleotide polymorphisms (SNPs) in N. seriolae that accurately distinguish: i) Vietnamese from non-Vietnamese strains, and ii) the two Vietnamese subclades. Using the mismatch amplification mutation assay (MAMA) format, we developed assays that genotype strains based on differences in amplicon melting temperature (melt-MAMA) and size (agarose-MAMA). Our MAMA assays accurately genotyped strains both from culture and fish tissues at low cost, using either real-time (~AUD$1/per sample) or conventional (~AUD$0.50/per sample) PCR instrumentation. Our novel assays provide a rapid, reproducible, and cost-effective tool for routine genotyping of this pathogen, allowing faster identification and treatment of nocardiosis-effected permit fish within Vietnamese aquaculture/mariculture facilities, an essential step in mitigating N. seriolae-associated losses.

BackgroundAntimicrobial resistance (AMR) is an ever-increasing global health concern. One crucial facet in tackling the AMR epidemic is earlier and more accurate AMR diagnosis, particularly in the dangerous and highly multi-drug resistant ESKAPE pathogen,Pseudomonas aeruginosa.ObjectivesWe aimed to develop two SYBR Green-based mismatch amplification mutation assays (SYBR-MAMAs) targeting GyrA T83I (gyrA248), and GyrA D87N, D87Y, and D87H (gyrA259). Together, these variants cause the majority of fluoroquinolone (FQ) AMR inP. aeruginosa.MethodsFollowing assay validation, thegyrA248 andgyrA259 SYBR-MAMAs were tested on 84 clinicalP. aeruginosaisolates from Queensland, Australia, 45 of which demonstrated intermediate/full ciprofloxacin resistance according to antimicrobial susceptibility testing.ResultsOur two SYBR-MAMAs correctly predicted an AMR phenotype in the majority (84%) of isolates with intermediate/full FQ resistance. Importantly, all FQ-sensitive strains were predicted to have a sensitive phenotype. Whole-genome sequencing confirmed 100% concordance with SYBR-MAMA genotypes.ConclusionsOur GyrA SYBR-MAMAs provide a rapid and cost-effective method for same-day identification of FQ AMR inP. aeruginosa. An additional SYBR-MAMA targeting the GyrB S466Y/S466F variants would increase FQ AMR prediction to 91%. Clinical implementation of our assays will permit more timely treatment alterations in cases where decreased FQ susceptibility is identified, leading to improved patient outcomes and antimicrobial stewardship.

Between 2010 and 2015, nocardiosis outbreaks caused by Nocardia seriolae affected many permit farms throughout Vietnam, causing mass fish mortalities. To understand the biology, origin, and epidemiology of these outbreaks, 20 N. seriolae strains collected from farms in four provinces in the South-Central Coast of Vietnam, along with two Taiwanese strains, were analysed using genetics and genomics. Pulsed-field gel electrophoresis identified a single cluster amongst all Vietnamese strains that was distinct from the Taiwanese strains. Like the PFGE findings, phylogenomic and single-nucleotide polymorphism (SNP) genotyping analyses revealed that all Vietnamese N. seriolae strains belonged to a single, unique clade. Strains fell into two subclades that differed by 103 SNPs, with almost no diversity within clades (0-2 SNPs). There was no association between geographic origin and subclade placement, suggesting frequent N. seriolae transmission between Vietnamese mariculture facilities during the outbreaks. Vietnamese strains shared a common ancestor with strains from Japan and China, with the closest strain, UTF1 from Japan, differing by just 217 SNPs from the Vietnamese ancestral node. Draft Vietnamese genomes range from 7.55-7.96 Mbp in size, have an average G+C content of 68.2%, and encode 7,602-7,958 predicted genes. Several putative virulence factors were identified, including genes associated with host cell adhesion, invasion, intracellular survival, antibiotic and toxic compound resistance, and haemolysin biosynthesis. Our findings provide important new insights into N. seriolae epidemiology and pathogenicity and will aid future vaccine development and disease management strategies, with the ultimate goal of nocardiosis-free aquaculture.

Cystic fibrosis (CF) and chronic obstructive pulmonary disease (COPD) are characterised by increasingly frequent acute pulmonary exacerbations that reduce life quality and length. Human airways are home to a rich polymicrobial environment, including 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/strains in healthy and diseased airways. Here, we used comparative genomics to develop a real-time PCR assay to permit rapid Prevotella spp. 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 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 gene presence 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 chronic respiratory diseases and healthy lungs.Data summaryThirty-five Prevotella spp. genomes generated in this study are available in the Sequence Read Archive (SRA) and GenBank databases under BioProject accession PRJNA742126.

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.

Background Antimicrobial resistance (AMR) poses a major threat to human health. Whole-genome sequencing holds great potential for AMR identification; however, there remain major gaps in accurately and comprehensively detecting AMR across the spectrum of AMR-conferring determinants and pathogens. Methods Using 16 wild-type Burkholderia pseudomallei and 25 with acquired AMR, we first assessed the performance of existing AMR software (ARIBA, CARD, ResFinder, and AMRFinderPlus) for detecting clinically relevant AMR in this pathogen. B. pseudomallei was chosen due to limited treatment options, high fatality rate, and AMR caused exclusively by chromosomal mutation (i.e. single-nucleotide polymorphisms [SNPs], insertions-deletions [indels], copy-number variations [CNVs], inversions, and functional gene loss). Due to poor performance with existing tools, we developed ARDaP ( A ntimicrobial R esistance D etection a nd P rediction) to identify the spectrum of AMR-conferring determinants in B. pseudomallei. Findings CARD, ResFinder, and AMRFinderPlus failed to identify any clinically-relevant AMR in B. pseudomallei; ARIBA identified AMR encoded by SNPs and indels that were manually added to its database. However, none of these tools identified CNV, inversion, or gene loss determinants, and ARIBA could not differentiate AMR determinants from natural genetic variation. In contrast, ARDaP accurately detected all SNP, indel, CNV, inversion, and gene loss AMR determinants described in B. pseudomallei (n≈50). Additionally, ARDaP accurately predicted three previously undescribed determinants. In mixed strain data, ARDaP identified AMR to as low as ~5% allelic frequency. Interpretation Existing AMR software packages are inadequate for chromosomal AMR detection due to an inability to detect resistance conferred by CNVs, inversions, and functional gene loss. ARDaP overcomes these major shortcomings. Further, ARDaP enables AMR prediction from mixed sequence data down to 5% allelic frequency, and can differentiate natural genetic variation from AMR determinants. ARDaP databases can be constructed for any microbial species of interest for comprehensive AMR detection. Funding National Health and Medical Research Council (BJC, EPP, DSS); Australian Government (DEM, ES); Advance Queensland (EPP, DSS).

Burkholderia pseudomallei is predominantly a tropical pathogen uncommonly found in the environment of temperate climatic regions. It is unclear if introduction into temperate regions is sporadic and temporary or if B. pseudomallei can persist in such environments. B. pseudomallei was identified in the environment of southwest Western Australia with melioidosis cases between 1966 and 1991. We report a new cluster with 23 animal fatalities in the same region from 2017, with B. pseudomallei again being recovered from the environment. Comparison of the isolates from the first and second clusters using genomics revealed a single sequence type, high clonality, and limited recombination, even though the time of recovery of the isolates spanned 51 years. This is a major contrast to the extensive genomic diversity seen in the tropics. Our data support the suggestion that B. pseudomallei has the ability to persist in nontropical environments, potentially in a latent state, and has the ability to activate following favorable conditions (rainfall) and then infect animals and humans., Burkholderia pseudomallei is the causative agent of the high-mortality disease melioidosis. Although melioidosis is classified as a tropical disease, rare autochthonous cases have been reported from temperate climatic regions, with uncertainty as to whether B. pseudomallei is persistent in the local environment and whether specific genetic mechanisms facilitate the survival of B. pseudomallei outside the tropics. Sporadic cases of melioidosis occurred in a valley region (latitude 31.6°S) in southwest Western Australia, Australia, between 1966 and 1992. We report a new melioidosis cluster in the same region following high rainfall in January 2017. More than 20 animals died, and B. pseudomallei was isolated from four alpacas, a parrot, and three environmental samples taken from the farm where the alpacas resided. Epidemiological data and genomics revealed that two locations on the farm were the probable sources of the alpaca infections. We determined that B. pseudomallei isolates from the 2017 cluster belonged to sequence type 284 (ST-284), as did all isolates recovered from 1966 to 1992. Genomic analysis confirmed that the ST-284 isolates were clonal and contained conserved genomic islands and limited evidence of recombination. We identified protein-coding regions unique to these isolates that might influence the persistence of B. pseudomallei in this temperate region. We demonstrate the environmental persistence of B. pseudomallei in a temperate region for over 50 years, with limited genetic changes suggesting a latent state and with activation, potential aerosolization, and local dispersal following unusually high rainfall. IMPORTANCE Burkholderia pseudomallei is predominantly a tropical pathogen uncommonly found in the environment of temperate climatic regions. It is unclear if introduction into temperate regions is sporadic and temporary or if B. pseudomallei can persist in such environments. B. pseudomallei was identified in the environment of southwest Western Australia with melioidosis cases between 1966 and 1991. We report a new cluster with 23 animal fatalities in the same region from 2017, with B. pseudomallei again being recovered from the environment. Comparison of the isolates from the first and second clusters using genomics revealed a single sequence type, high clonality, and limited recombination, even though the time of recovery of the isolates spanned 51 years. This is a major contrast to the extensive genomic diversity seen in the tropics. Our data support the suggestion that B. pseudomallei has the ability to persist in nontropical environments, potentially in a latent state, and has the ability to activate following favorable conditions (rainfall) and then infect animals and humans.

The Elizabethkingia genus has gained global attention in recent years as a sporadic, worldwide, nosocomial pathogen. Elizabethkingia spp. are intrinsically multidrug resistant, primarily infect immunocompromised individuals, and are associated with high mortality (∼20-40%). Yet, gaps remain in our understanding of transmission, global strain relatedness, antimicrobial resistance and effective therapy. Over a 16-year period 22 clinical and six hospital environmental isolates were collected from Queensland, Australia. Identification using the MALDI-TOF MS (VITEK® MS) and whole-genome sequencing was compared with a global strain dataset. Phylogenomic reconstruction robustly identified 22 E. anophelis, three E. miricola, two E. meningoseptica and one E. bruuniana, most of which branched as unique lineages. Global analysisrevealed some Australian E. anophelis isolates are genetically closely related to strains identified from the USA, England and Asia. Comparative genomics of clinical and environmental strains identified evidence of nosocomial transmission in patients, indicating probable infection from a hospital reservoir. Furthermore, broth microdilution against 39 antimicrobials revealed almost ubiquitous resistance to aminoglycosides, carbapenems, cephalosporins and penicillins. Like other international strains, our isolates expressed susceptibility to minocycline and levofloxacin and the less common trimethoprim/sulfamethoxazole. Our study demonstrates important new insights into the genetic diversity, environmental persistence, transmission of and potential effective therapy for Australian Elizabethkingia species.

The Elizabethkingia genus has gained global attention in recent years as containing sporadic, worldwide, nosocomial pathogens. Elizabethkingia spp. are intrinsically multidrug resistant, primarily infect immunocompromised individuals, and are associated with high mortality (∼20 to 40%). As yet, gaps remain in our understanding of transmission, global strain relatedness, antimicrobial resistance, and effective therapy. Over a 16-year period, 22 clinical and 6 hospital environmental isolates were collected from Queensland, Australia. Identification using matrix-assisted laser desorption ionization–time of flight mass spectrometry (MALDI-TOF MS) (Vitek MS) and whole-genome sequencing was compared with a global strain data set. Phylogenomic reconstruction robustly identified 22 Elizabethkingia anophelis, 3 Elizabethkingia miricola, 2 Elizabethkingia meningoseptica, and 1 Elizabethkingia bruuniana isolates, most of which branched as unique lineages. Global analysis revealed that some Australian E. anophelis isolates are genetically closely related to strains from the United States, England, and Asia. Comparative genomics of clinical and environmental strains identified evidence of nosocomial transmission in patients, indicating probable infection from a hospital reservoir. Furthermore, broth microdilution against 39 antimicrobials revealed almost ubiquitous resistance to aminoglycosides, carbapenems, cephalosporins, and penicillins. Like other international strains, our isolates expressed susceptibility to minocycline and levofloxacin and the less common trimethoprim-sulfamethoxazole. Our study demonstrates important new insights into the genetic diversity, environmental persistence, and transmission of and potential effective therapy for Australian Elizabethkingia species.

The Elizabethkingia genus has gained global attention in recent years as a nosocomial pathogen. Elizabethkingia spp. are intrinsically multidrug resistant, primarily infect immunocompromised individuals, and are associated with high mortality (∼20-40%). Although Elizabethkingia infections appear sporadically worldwide, gaps remain in our understanding of transmission, global strain relatedness and patterns of antimicrobial resistance. To address these knowledge gaps, 22 clinical isolates collected in Queensland, Australia, over a 16-year period along with six hospital environmental isolates were examined using MALDI-TOF MS (VITEK® MS) and whole-genome sequencing to compare with a global strain dataset. Phylogenomic reconstruction against all publicly available genomes (n=100) robustly identified 22 E. anophelis, three E. miricola, two E. meningoseptica and one E. bruuniana from our isolates, most with previously undescribed diversity. Global relationships show Australian E. anophelis isolates are genetically related to those from the USA, England and Asia, suggesting shared ancestry. Genomic examination of clinical and environmental strains identified evidence of nosocomial transmission in patients admitted several months apart, indicating probable infection from a hospital reservoir. Furthermore, broth microdilution of the 22 clinical Elizabethkingia spp. isolates against 39 antimicrobials revealed almost ubiquitous resistance to aminoglycosides, carbapenems, cephalosporins and penicillins, but susceptibility to minocycline, levofloxacin and trimethoprim/sulfamethoxazole. Our study demonstrates important new insights into the genetic diversity, environmental persistence and transmission of Australian Elizabethkingia species. Furthermore, we show that Australian isolates are highly likely to be susceptible to minocycline, levofloxacin and trimethoprim/sulfamethoxazole, suggesting that these antimicrobials may provide effective therapy for Elizabethkingia infections.ImportanceElizabethkingia are a genus of environmental Gram-negative, multidrug resistant, opportunistic pathogens. Although an uncommon cause of nosocomial and community-acquired infections, Elizabethkingia spp. are known to infect those with underlying co-morbidities and/or immunosuppression, with high mortality rates of ∼20-40%. Elizabethkingia have a presence in Australian hospitals and patients; however, their origin, epidemiology, and antibiotic resistance profile of these strains is poorly understood. Here, we performed phylogenomic analyses of clinical and hospital environmental Australian Elizabethkingia spp., to understand transmission and global relationships. Next, we performed extensive minimum inhibitory concentration testing to determine antimicrobial susceptibility profiles. Our findings identified a highly diverse Elizabethkingia population in Australia, with many being genetically related to international strains. A potential transmission source was identified within the hospital environment where two transplant patients were infected and three E. anophelis strains formed a clonal cluster within the phylogeny. Furthermore, near ubiquitous susceptibility to tetracyclines, fluoroquinolones and trimethoprim/sulfamethoxazole was observed in clinical isolates. We provide new insights into the origins, transmission and epidemiology of Elizabethkingia spp., in addition to understanding their intrinsic resistance profiles and potential effective treatment options, which has implications to managing infections and detecting outbreaks globally.

SynopsisThird-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. Alongside routine mass spectrometry speciation and antibiotic sensitivity testing, isolates were examined using rapid (∼40 min) real-time PCR assays targeting the most common extended spectrum β-lactamases (ESBLs; CTX-M-1 and CTX-M-9 groups, plus TEM, SHV, and an internal 16S ribosomal DNA control). AmpC CMY β-lactamase 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%) encoding blaCTX-M-1 group, 92 (31.9%) blaCTX-M-9 group, 48 (16.7%) blaSHV, 133 (46.2%) blaTEM, and 34 (11.8%) blaCMY. 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 ESBLs and carbapenemases, including blaNDM and blaIMP, 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.

Several members of the Gram-negative environmental bacterial genus Achromobacter are associated with serious infections, with Achromobacter xylosoxidans being the most common. Despite their pathogenic potential, little is understood about these intrinsically drug-resistant bacteria and their role in disease, leading to suboptimal diagnosis and management. Here, we performed comparative genomics for 158 Achromobacter spp. genomes to robustly identify species boundaries, reassign several incorrectly speciated taxa and identify genetic sequences specific for the genus Achromobacter and for A. xylosoxidans . Next, we developed a Black Hole Quencher probe-based duplex real-time PCR assay, Ac-Ax, for the rapid and simultaneous detection of Achromobacter spp. and A. xylosoxidans from both purified colonies and polymicrobial clinical specimens. Ac-Ax was tested on 119 isolates identified as Achromobacter spp. using phenotypic or genotypic methods. In comparison to these routine diagnostic methods, the duplex assay showed superior identification of Achromobacter spp. and A. xylosoxidans , with five Achromobacter isolates failing to amplify with Ac-Ax confirmed to be different genera according to 16S rRNA gene sequencing. Ac-Ax quantified both Achromobacter spp. and A. xylosoxidans down to ~110 genome equivalents and detected down to ~12 and ~1 genome equivalent(s), respectively. Extensive in silico analysis, and laboratory testing of 34 non- Achromobacter isolates and 38 adult cystic fibrosis sputa, confirmed duplex assay specificity and sensitivity. We demonstrate that the Ac-Ax duplex assay provides a robust, sensitive and cost-effective method for the simultaneous detection of all Achromobacter spp. and A. xylosoxidans and will facilitate the rapid and accurate diagnosis of this important group of pathogens.

The heterogeneous and highly recombinogenic genus Haemophilus comprises several species, some of which are pathogenic to humans. All share an absolute requirement for blood-derived factors during growth. Certain species, such as the pathogen Haemophilus influenzae and the commensal Haemophilus haemolyticus, are thought to require both haemin (X-factor) and nicotinamide adenine dinucleotide (NAD, V-factor), whereas others, such as the informally classified 'Haemophilus intermedius subsp. intermedius', and Haemophilus parainfluenzae, only require V-factor. These differing growth requirements are commonly used for species differentiation, although a number of studies are now revealing issues with this approach. Here, we perform large-scale phylogenomics of 240 Haemophilus spp. genomes, including five 'H. intermedius' genomes generated in the current study, to reveal that strains of the 'H. intermedius' group are in fact haemin-independent H. haemolyticus (hiHh). Closer examination of these hiHh strains revealed that they encode an intact haemin biosynthesis pathway, unlike haemin-dependent H. haemolyticus and H. influenzae, which lack most haemin biosynthesis genes. Our results suggest that the common ancestor of modern-day H. haemolyticus and H. influenzae lost key haemin biosynthesis loci, likely as a consequence of specialized adaptation to otorhinolaryngeal and respiratory niches during their divergence from H. parainfluenzae. Genetic similarity analysis demonstrated that the haemin biosynthesis loci acquired in the hiHh lineage were likely laterally transferred from a H. parainfluenzae ancestor, and that this event probably occurred only once in hiHh. This study further challenges the validity of phenotypic methods for differentiating among Haemophilus species, and highlights the need for whole-genome sequencing for accurate characterization of species within this taxonomically challenging genus.

Background Burkholderia pseudomallei, the causative agent of the high-mortality disease melioidosis, is a gram-negative bacterium that is naturally resistant to many antibiotics. There is no vaccine for melioidosis, and effective eradication is reliant on biphasic and prolonged antibiotic administration. The carbapenem drug meropenem is the current gold standard option for treating severe melioidosis. Intrinsic B. pseudomallei resistance toward meropenem has not yet been documented; however, resistance could conceivably develop over the course of infection, leading to prolonged sepsis and treatment failure. Methods We examined our 30-year clinical collection of melioidosis cases to identify B. pseudomallei isolates with reduced meropenem susceptibility. Isolates were subjected to minimum inhibitory concentration (MIC) testing toward meropenem. Paired isolates from patients who had evolved decreased susceptibility were subjected to whole-genome sequencing. Select agent-compliant genetic manipulation was carried out to confirm the molecular mechanisms conferring resistance. Results We identified 11 melioidosis cases where B. pseudomallei isolates developed decreased susceptibility toward meropenem during treatment, including 2 cases not treated with this antibiotic. Meropenem MICs increased from 0.5-0.75 µg/mL to 3-8 µg/mL. Comparative genomics identified multiple mutations affecting multidrug resistance-nodulation-division (RND) efflux pump regulators, with concomitant overexpression of their corresponding pumps. All cases were refractory to treatment despite aggressive, targeted therapy, and 2 were associated with a fatal outcome. Conclusions This study confirms the role of RND efflux pumps in decreased meropenem susceptibility in B. pseudomallei. These findings have important ramifications for the diagnosis, treatment, and management of life-threatening melioidosis cases.

Bacterial infections are the primary cause of respiratory decline and mortality in cystic fibrosis (CF) patients. In a recent study, Diaz Caballero and colleagues [1] (PLoS Pathog. 2018;14:e1007453) catalogued the molecular adaptation of a decade-long Burkholderia multivorans infection in a Canadian CF patient, which evolved to become resistant towards multiple classes of antibiotics.

Genomics has fundamentally transformed our understanding of infectious diseases. Here, we provide an overview of the insights gained from genomic investigations of the melioidosis-causing pathogen, Burkholderia pseudomallei. B. pseudomallei has a large and complex genome that encodes an impressive array of virulence factors, some of which are variably present. Despite frequent lateral gene transfer, phylogenomics resolves B. pseudomallei populations to the continental level, enabling strain source tracing in non-endemic regions and shedding light on the origin and anthropogenic spread of B. pseudomallei populations across tropical and subtropical regions worldwide. Within-host evolution, outbreak point-sources and antibiotic resistance can now be rapidly elucidated using comparative genomic approaches. Genomics is an essential tool for understanding the phylogeography, transmission, evolution, virulence, epidemiology, and antibiotic resistance of B. pseudomallei. A shift towards more detailed characterization of genetic variation using RNA-Seq, Methyl-Seq, and genome-wide association studies will provide additional insights into this fascinating and deadly bacterium.

These strains have heightened pathogenic potential for rapid dissemination to multiple tissues, including the central nervous system., Neurologic melioidosis is a serious, potentially fatal form of Burkholderia pseudomallei infection. Recently, we reported that a subset of clinical isolates of B. pseudomallei from Australia have heightened virulence and potential for dissemination to the central nervous system. In this study, we demonstrate that this subset has a B. mallei–like sequence variation of the actin-based motility gene, bimA. Compared with B. pseudomallei isolates having typical bimA alleles, isolates that contain the B. mallei–like variation demonstrate increased persistence in phagocytic cells and increased virulence with rapid systemic dissemination and replication within multiple tissues, including the brain and spinal cord, in an experimental model. These findings highlight the implications of bimA variation on disease progression of B. pseudomallei infection and have considerable clinical and public health implications with respect to the degree of neurotropic threat posed to human health.

Background: Extended spectrum beta‐lactamase (ESBL) and AmpC‐producing Enterobacterales are prevalent worldwide. Poor clinical outcomes, including death, are more common in patients both colonised and infected with these organisms1. Australia is thought to have a comparatively low prevalence of such organisms, however no clear continuous surveillance system exists2. Specific ESBL type has shown significant geographic variation, having different therapeutic and infection control implications for different regions3. Little is known about the prevalence, molecular type and spread of ESBL and AmpC‐producing Enterobacterales among hospitalised patients in the Sunshine Coast region. Aims: To (i) Determine the prevalence and molecular epidemiology of ESBL and AmpC‐producing Enterobacterales, (ii) Develop an in‐house multiplex polymerase chain reaction (PCR) assay to detect common ESBL genes, (iii) Compare our PCR assay to the double disc diffusion method for detection of ESBLs, and to (iv) Establish, with whole genome sequencing (WGS), whether transmission was occurring in our intensive care unit (ICU). Methods: Clinical specimens from adult and paediatric patients culturing Enterobacterales resistant to a third‐generation cephalosporin were collected prospectively from July 2017 to August 2018 at the Sunshine Coast University Hospital laboratory. Antibiotic susceptibility testing was determined using VITEK‐2®. A validated in‐house pentaplex PCR assay was developed and used to detect common beta‐lactamase genes (CTX‐1, CTX‐9, SHV, TEM, CMY). All isolates collected from patients admitted to the ICU were identified and underwent WGS analysis using the Illumina® platform. Results: 295 clinical specimens were identified. Escherichia coli (78.8%), Klebsiella pneumoniae (16.8%) with 1‐2% consisting of Enterobacter spp., Serratia marcescens and Citrobacter freundii. Specimens included urine (53.8%), rectal swab (39.4%), pus (3.0%) and blood (2.4%). Carbapenem resistance was observed in 3.3%. ICU patient specimens made up 18.2%. The majority (87.8%) of isolates had at least one gene identified on PCR testing of which 48.5% produced CTX‐1, 31.2% produced CTX‐9, 15.9% produced SHV, 12.5% produced CMY and 40.6% produced TEM. Approximately forty percent of isolates produced more than one ESBL with the most common being CTX‐1 and TEM coproduction. A small amount (

Human-to-human transmission of the melioidosis bacterium,Burkholderia pseudomallei, is exceedingly rare, with only a handful of suspected cases documented to date. Here, we used whole-genome sequencing (WGS) to characterise one such unusualB. pseudomalleitransmission event, which occurred between a breastfeeding mother with mastitis and her child. Two strains corresponding to multilocus sequence types (STs) 259 and 261 were identified in the mother’s sputum from both the primary culture sweep and in purified colonies, confirming an unusual polyclonal infection in this patient. In contrast, primary culture sweeps of the mother’s breast milk and the child’s cerebrospinal fluid and blood samples contained only ST-259, indicating monoclonal transmission to the child. Analysis of purified ST-259 isolates showed no genetic variation between mother and baby isolates, providing the strongest possible evidence ofB. pseudomalleitransmission, probably via breastfeeding. Next, phylogenomic analysis of all isolates, including the mother’s mixed ST-259/261 sputum sample was performed to investigate the effects of mixtures on phylogenetic inference. Inclusion of this mixture caused a dramatic reduction in the number of informative SNPs, resulting in branch collapse of ST-259 and ST-261 isolates, and several instances of incorrect topology in a globalB. pseudomalleiphylogeny, resulting in phylogenetic incongruence. Although phylogenomics can provide clues about the presence of mixtures within WGS datasets, our results demonstrate that this methodology can lead to phylogenetic misinterpretation if mixed genomes are not correctly identified and omitted. Using current bioinformatic tools, we demonstrate a robust method for bacterial mixture identification and strain parsing that avoids these pitfalls.Impact StatementBurkholderia pseudomalleiis the causative agent of melioidosis, a tropical disease of high mortality.B. pseudomalleiinfections occur almost exclusively through contact with contaminated soil and water. Using whole-genome sequencing (WGS), we investigated a rare case of suspectedB. pseudomalleitransmission from mother to child. The mother’s sputum, breast milk and the baby’s blood and cerebrospinal fluid (CSF) specimens were collected, and DNA was extracted from both pure colonies and primary culture sweeps to capture potential strain mixtures. In-depth analysis of genetic variants identified two strains in the mother’s sputum belonging to multilocus sequence types ST-259 and ST-261, whereas the child was infected with only ST-259. Comparative genomics revealed no genetic differences between mother and child ST-259 isolates, providing the strongest possible evidence of transmission to the child via breast milk. The sputum strain mixture was subsequently used to develop a bioinformatic method for identification and quantification of mixtures from WGS data. Using this method, we found ST-259 and ST-261 at an 87%:13% ratio, respectively. Finally, we demonstrate the negative impact that even a single strain mixture event can have on both within-ST and global phylogenomic inferences. Our findings highlight the need for bioinformatic quality control to avoid unintended consequences of phylogenomic incongruence and branch collapse.Data SummaryWhole-genome sequencing data have been deposited in the NCBI Sequence Read Archive (SRA) and GenBank under BioProject accession numberPRJNA559002.The GenBank accession number for MSHR0643 assembly is VXLH00000000.1.The SRA accession numbers for all raw sequence data are listed in Table 1.RepositoriesAll sequencing data generated as part of this study can be found under the NCBI BioProject PRJNA559002 with accession numbers listed in Table 1.

Stenotrophomonas maltophilia is emerging as an important cause of disease in nosocomial and community-acquired settings, including bloodstream, wound and catheter-associated infections. Cystic fibrosis (CF) airways also provide optimal growth conditions for various opportunistic pathogens with high antibiotic tolerance, including S. maltophilia . Currently, there is no rapid, cost-effective and accurate molecular method for detecting this potentially life-threatening pathogen, particularly in polymicrobial specimens, suggesting that its true prevalence is underestimated. Here, we used large-scale comparative genomics to identify a specific genetic target for S. maltophilia , with subsequent development and validation of a real-time PCR assay for its detection. Analysis of 167 Stenotrophomonas spp. genomes identified a conserved 4 kb region in S. maltophilia , which was targeted for Black Hole Quencher assay design. Our assay yielded the positive detection of 89 of 89 (100%) clinical S. maltophilia strains, and no amplification of 23 non- S. maltophilia clinical isolates. S. maltophilia was detected in 10 of 16 CF sputa, demonstrating the assay's utility for direct detection in respiratory specimens. The assay demonstrated good sensitivity, with limits of detection and quantitation on pure culture of ~10 and ~100 genome equivalents, respectively. Our assay provides a highly specific, sensitive and cost-effective method for the accurate identification of S. maltophilia , and will improve the diagnosis and treatment of this under-recognized pathogen by enabling its accurate and rapid detection from polymicrobial clinical and environmental samples.

BackgroundAntimicrobial resistance (AMR) poses a major threat to human health. Whole-genome sequencing holds great potential for AMR identification; however, there remain major gaps in comprehensively detecting AMR across the spectrum of AMR-conferring determinants and pathogens.MethodsUsing 16 wild-type Burkholderia pseudomallei and 25 with acquired AMR, we first assessed the performance of existing AMR software (ARIBA and CARD) for detecting clinically relevant AMR in this pathogen. B. pseudomallei was chosen due to limited treatment options, high fatality rate, and AMR caused exclusively by chromosomal mutation (i.e. single-nucleotide polymorphisms [SNPs], insertions-deletions [indels], copy-number variations [CNVs], and functional gene loss). Due to poor performance with existing tools, we developed ARDaP (Antimicrobial Resistance Detection and Prediction) to identify the spectrum of AMR-conferring determinants in B. pseudomallei.ResultsCARD failed to identify any clinically-relevant AMR in B. pseudomallei, ARIBA cannot differentiate AMR determinants from natural genetic variation, and neither CARD or ARIBA can identify CNV or gene loss determinants. In contrast, ARDaP accurately detected all SNP, indel, CNV, and gene loss AMR determinants described in B. pseudomallei (n≈50). Additionally, ARDaP accurately predicted three previously undescribed determinants. In mixed strain data, ARDaP identified AMR to as low as ~5% allelic frequency.ConclusionsWe demonstrate that existing AMR software are inadequate for comprehensive AMR detection; ARDaP overcomes the shortcomings of existing tools. Further, ARDaP enables AMR prediction from mixed sequence data down to 5% allelic frequency. ARDaP databases can be constructed for any microbial species of interest for comprehensive AMR detection.

Stenotrophomonas maltophilia is emerging as an important cause of disease in nosocomial and community-acquired settings, including bloodstream, wound and catheter-associated infections. Cystic fibrosis airways also provide optimal growth conditions for various opportunistic pathogens with high antibiotic tolerance, including S. maltophilia. Currently, there is no rapid, cost-effective, and accurate molecular method for detecting this potentially life-threatening pathogen, particularly in polymicrobial specimens, suggesting that its true prevalence may be underestimated. Here, we used large-scale comparative genomics to identify a specific genetic target for S. maltophilia, with subsequent development and validation of a real-time PCR assay for its detection. Analysis of 165 Stenotrophomonas spp. genomes identified a 4kb region specific to S. maltophilia, which was targeted for Black Hole Quencher assay design. Our assay yielded the positive detection of 89 of 89 (100%) clinical S. maltophilia strains, and no amplification of 23 non-S. maltophilia clinical isolates. S. maltophilia was detected in 10/16 CF sputa, demonstrating the utility for direct detection in respiratory specimens. The assay demonstrated good sensitivity, with limits of detection and quantitation on pure culture of ~10 and ~100 genome equivalents, respectively. Our assay provides a highly specific, sensitive, and cost-effective method for the accurate identification of S. maltophilia, and will improve the diagnosis and treatment of this under-recognized pathogen by enabling its accurate and rapid detection from polymicrobial clinical and environmental samples.

Although acute melioidosis is the most common outcome of Burkholderia pseudomallei infection, we have documented a case, P314, where disease severity lessened with time, and the pathogen evolved towards a commensal relationship with the host. In the current study, we used whole-genome sequencing to monitor this long-term symbiotic relationship to better understand B. pseudomallei persistence in P314’s sputum despite intensive initial therapeutic regimens. We collected and sequenced 118 B. pseudomallei isolates from P314’s airways over a >16-year period, and also sampled the patient’s home environment, recovering six closely related B. pseudomallei isolates from the household water system. Using comparative genomics, we identified 126 SNPs in the core genome of the 124 isolates or 162 SNPs/indels when the accessory genome was included. The core SNPs were used to construct a phylogenetic tree, which demonstrated a close relationship between environmental and clinical isolates and detailed within-host evolutionary patterns. The phylogeny had little homoplasy, consistent with a strictly clonal mode of genetic inheritance. Repeated sampling revealed evidence of genetic diversification, but frequent extinctions left only one successful lineage through the first four years and two lineages after that. Overall, the evolution of this population is nonadaptive and best explained by genetic drift. However, some genetic and phenotypic changes are consistent with in situ adaptation. Using a mouse model, P314 isolates caused greatly reduced morbidity and mortality compared to the environmental isolates. Additionally, potentially adaptive phenotypes emerged and included differences in the O-antigen, capsular polysaccharide, motility, and colony morphology. The >13-year co-existence of two long-lived lineages presents interesting hypotheses that can be tested in future studies to provide additional insights into selective pressures, niche differentiation, and microbial adaptation. This unusual melioidosis case presents a rare example of the evolutionary progression towards commensalism by a highly virulent pathogen within a single human host., Author summary Pathogens frequently jump between different hosts, and associated adaptation may lead to the emergence of new infectious agents. Such host-jumping evolution is witnessed through endpoint analyses but these cannot capture genetic changes in lineages that have gone extinct. In this study, we have identified and monitored an example of the evolution of a bacterium often deadly to its mammalian host, in an unprecedented case whereby disease lessened through time and the pathogen became a part of the commensal human flora. We used genomic analyses to characterize more than 16 years of this evolutionary process and the stepwise mutations that control pathogen interactions with the patient. Soon after infection, mutational changes occurred that allowed the bacterium to remain in the airways without causing disease. This shift towards avirulence was determined based on clinical data and virulence testing in an animal model. In addition, mutations occurred that contributed to the persistence of the bacteria in the patient's lungs. Finally, we found evidence for the evolutionary emergence and persistence of two distinct lineages of the bacterium over the last 13 years, presenting interesting questions about niche utilization. Bacteria are ubiquitous in the human body and almost all are beneficial or benign. In this study, we document the evolutionary conversion of a normally deadly bacterium towards a commensal.

Non-typeable Haemophilus influenzae (NTHi), an opportunistic pathogen of the upper airways of healthy children, can infect the lower airways, driving chronic lung disease. However, the molecular basis underpinning NTHi transition from a commensal to a pathogen is not clearly understood. Here, we performed comparative genomic and transcriptomic analyses of 12 paired, isogenic NTHi strains, isolated from the nasopharynx (NP) and bronchoalveolar lavage (BAL) of 11 children with chronic lung disease, to identify convergent molecular signatures associated with lung adaptation. Comparative genomic analyses of the 12 NP-BAL pairs demonstrated that five were genetically identical, with the remaining seven differing by only 1 to 3 mutations. Within-patient transcriptomic analyses identified between 2 and 58 differentially expressed genes in 8 of the 12 NP-BAL pairs, including pairs with no observable genomic changes. Whilst no convergence was observed at the gene level, functional enrichment analysis revealed significant under-representation of differentially expressed genes belonging to Coenzyme metabolism, Function unknown, Translation, ribosomal structure and biogenesis Cluster of Orthologous Groups categories. In contrast, Carbohydrate transport and metabolism, Cell motility and secretion, Intracellular trafficking and secretion, and Energy production categories were over-represented. This observed trend amongst genetically-unrelated NTHi strains provides evidence of convergent transcriptional adaptation of NTHi to paediatric airways that deserves further exploration. Understanding the pathoadaptative mechanisms that NTHi employs to infect and persist in the lower paediatric airways is essential for devising targeted diagnostics and treatments aimed at minimising disease severity, and ultimately, preventing NTHi lung infections and subsequent chronic lung disease in children.

The Tier 1 select agentBurkholderia pseudomalleiis an environmental bacterium that causes melioidosis, a high mortality disease. Variably present genetic markers used to elucidate strain origin, relatedness and virulence inB. pseudomalleiinclude theBurkholderiaintracellular motility factor A (bimA) and filamentous hemagglutinin 3 (fhaB3) gene variants. Three lipopolysaccharide (LPS) O-antigen types inB. pseudomalleihave been described, which vary in proportion between Australian and Asian isolates. However, it remains unknown if these LPS types can be used as genetic markers for geospatial analysis within a contiguous melioidosis-endemic region. Using a combination of whole-genome sequencing (WGS), statistical analysis and geographical mapping, we examined if the LPS types can be used as geographical markers in the Northern Territory, Australia. The clinical isolates revealed that LPS A prevalence was highest in the Darwin and surrounds (n = 660; 96% being LPS A and 4% LPS B) and LPS B in the Katherine and Katherine remote and East Arnhem regions (n = 79; 60% being LPS A and 40% LPS B). Bivariate logistics regression of 999 clinicalB. pseudomalleiisolates revealed that the odds of getting a clinical isolate with LPS B was highest in East Arnhem in comparison to Darwin and surrounds (OR 19.5, 95% CI 9.1 – 42.0;pbimAgenotypebimABmhas a similar remote region geographical footprint to that of LPS B. In addition, correlation of LPS type with multi-locus sequence typing (MLST) was strong, and where multiple LPS types were identified within a single sequence type, WGS confirmed homoplasy of the MLST loci. The clinical, sero-diagnostic and vaccine implications of geographically-basedB. pseudomalleiLPS types, and their relationships to regional and global dispersal of melioidosis, require global collaborations with further analysis of larger clinically and geospatially-linked datasets.Author SummaryBurkholderia pseudomalleiis a pathogenic soil bacterium that causes the disease melioidosis, which occurs in many tropical regions globally and in recent years has emerged in non-tropical regions. Melioidosis has been predicted to affect 165,000 people every year resulting in an estimated 89,000 deaths. Person to person transmission is rare with most cases linked to exposure to the bacterium from the environment. The genetic background ofB. pseudomalleihas been well studied and variably present genes have been linked to distinct melioidosis disease states and geographic regions, however we still need a stronger understanding of the association of genes with geography. Three lipopolysaccharide types exist inB. pseudomalleiand the prevalence of the lipopolysaccharide genes vary between melioidosis endemic regions, but it is unknown if the lipopolysaccharide genes can be used as geographical markers in a single melioidosis-endemic region. In this study, we used a combination of whole-genome sequencing, statistics and geographical mapping to elucidate if the three lipopolysaccharide genes can be used as geographical markers within the Northern Territory, Australia. We show that the three LPS types have distinct but overlapping geographical footprints within a single melioidosis region and can be used as geographic markers alongside a number of other important variably presentB. pseudomalleigenes.

Background The causative agent of melioidosis is the Gram-negative bacterium Burkholderia pseudomallei. Clinical presentations of melioidosis are notably diverse, with host risk factors considered central to progression from infection to disease and clinical outcome. Ubiquitous and variably present virulence determinants have been described for B pseudomallei, with several variably present minority genotypes associated with specific disease presentations. The lipopolysaccharide (LPS) O-antigen of B pseudomallei is highly diverse with 3 types described. In vitro data suggest differential virulence between LPS types, but it remains unclear whether this LPS O-antigen diversity influences clinical presentation, severity, and outcomes in patients with melioidosis. Methods Whole-genome sequencing was performed to assign an LPS type to 1005 consecutive B pseudomallei strains, each corresponding to a melioidosis patient enrolled in the 28-year Darwin Prospective Melioidosis study. Correlations of LPS genotype with clinical parameters was then undertaken. Results Bivariate analysis demonstrated that mortality and the rates of bacteremia and septic shock were the same for patients with the 2 predominant B pseudomallei LPS genotypes A (87% of cases) and B (12% of all cases). Mortality was 12% and 12%, bacteremia was 57% and 53%, and septic shock was 22% and 18% for LPS A and LPS B, respectively. Conclusions Lipopolysaccharide genotype was not associated with melioidosis severity or outcome. These findings suggest that in vitro differential virulence between B pseudomallei LPS genotypes does not translate to clinical significance, and this supports the primary role of host risk factors in determining disease severity and outcomes in melioidosis.