Improved multilocus sequence typing of Burkholderia pseudomallei and closely related species

Multi-locus sequence typing and genetic diversity of antibiotic-resistant genes and virulence-associated genes in Burkholderia pseudomallei: Insights from whole genome sequencing of animal and environmental isolates in Thailand

Burkholderia pseudomallei is a Gram-negative bacillus and the etiological agent of melioidosis in humans and animals. The disease is highly endemic in northern Australia and Southeast Asia. Comprehensive genomic data are essential for understanding the bacteria's dissemination and genetic relationships among strains from different geographical regions. In this study, we conducted antimicrobial susceptibility testing and whole-genome sequencing of 54 B. pseudomallei isolates obtained from environmental and animal sources in southern Thailand between 2011 and 2018. Their genomics were determined of antibiotic-resistant genes (ARGs), virulence-associated genes, mobile genetic elements (MGEs), sequence types (STs), and single nucleotide polymorphisms (SNPs) to evaluate their epidemiological relatedness. Remarkably, all 54 isolates displayed sensitivity to antimicrobial agents typically used for melioidosis treatment. We identified nine distinct sequence types: ST392, ST51, ST409, ST508, ST376, ST1721, ST389, ST395, and ST289. Oxacillinase genes and the resistance nodulation family of efflux pumps (RND) were identified as contributors to antimicrobial resistance. Phylogenetic analysis demonstrated close genetic relations with other strains isolated from Southeast Asia. Furthermore, 172 virulence-associated genes were identified among the isolates, suggesting variations in clinical presentations. These findings underscore the importance of ongoing molecular genetic surveillance of B. pseudomallei for effective healthcare management and reducing melioidosis mortality.

Virulome and phylogenomic profiling of a novel Burkholderia pseudomallei strain from an Indian clinical isolate

Highly pathogenic Burkholderia pseudomallei is the causative agent of melioidosis, a neglected tropical disease endemic in Southeast Asian tropical region. This bacterium encompasses diverse virulence factors which further undergo dynamic gene-expression flux as it transits through distinct environmental niches within the host which may lead to manifestation of differential clinical symptoms. B. pseudomallei, is classified as a Tier 1 select agent in the United States and regarded as a risk group 3 organism in India with the potential to be used as bioweapon. Considering these facts, it is vital to uncover both physiological and genetic heterogeneity of B. pseudomallei, particularly to identify any novel virulence factors that may contribute to pathogenicity. B. pseudomallei strain CM000113 was isolated from a clinical case in India, characterized it for its physiological, biochemical, and prominently genetic traits through WGS. It has a type 2 morphotype with faster doubling time and high biofilm producing capacity as compared to Pseudomonas aeruginosa. The genome size is 7.3 Mbp and it is phylogenetically close to B. pseudomallei strain Mahidol 1106a and Burkholderia mallei Turkey 2. We observed genetic heterogeneity, as key virulence factors that were identified shows sequence dissimilarity with reference strains. Additionally, presence of genomic islands, harbouring two virulence factors, GmhA and GmhB2, associated with pathogenesis indicates possibility of horizontal gene transfer. These results emphasize the need for an extensive study focusing the genome of B. pseudomallei and its associated heterogeneity, to identify molecular biomarkers aiding to develop point-of-care diagnostic kits for early diagnosis of melioidosis.

The Darwin Prospective Melioidosis Study: a 30-year prospective, observational investigation

BACKGROUND

The global distribution of melioidosis is under considerable scrutiny, with both unmasking of endemic disease in African and Pacific nations and evidence of more recent dispersal in the Americas. Because of the high incidence of disease in tropical northern Australia, The Darwin Prospective Melioidosis Study commenced in October, 1989. We present epidemiology, clinical features, outcomes, and bacterial genomics from this 30-year study, highlighting changes in the past decade.

METHODS

The present study was a prospective analysis of epidemiological, clinical, and laboratory data for all culture-confirmed melioidosis cases from the tropical Northern Territory of Australia from Oct 1, 1989, until Sept 30, 2019. Cases were identified on the basis of culture-confirmed melioidosis, a laboratory-notifiable disease in the Northern Territory of Australia. Patients who were culture-positive were included in the study. Multivariable analysis determined predictors of clinical presentations and outcome. Incidence, survival, and cluster analyses were facilitated by population and rainfall data and genotyping of Burkholderia pseudomallei, including multilocus sequence typing and whole-genome sequencing.

FINDINGS

There were 1148 individuals with culture-confirmed melioidosis, of whom 133 (12%) died. Median age was 50 years (IQR 38-60), 48 (4%) study participants were children younger than 15 years of age, 721 (63%) were male individuals, and 600 (52%) Indigenous Australians. All but 186 (16%) had clinical risk factors, 513 (45%) had diabetes, and 455 (40%) hazardous alcohol use. Only three (2%) of 133 fatalities had no identified risk. Pneumonia was the most common presentation occurring in 595 (52%) patients. Bacteraemia occurred in 633 (56%) of 1135 patients, septic shock in 240 (21%) patients, and 180 (16%) patients required mechanical ventilation. Cases correlated with rainfall, with 80% of infections occurring during the wet season (November to April). Median annual incidence was 20·5 cases per 100 000 people; the highest annual incidence in Indigenous Australians was 103·6 per 100 000 in 2011-12. Over the 30 years, annual incidences increased, as did the proportion of patients with diabetes, although mortality decreased to 17 (6%) of 278 patients over the past 5 years. Genotyping of B pseudomallei confirmed case clusters linked to environmental sources and defined evolving and new sequence types.

INTERPRETATION

Melioidosis is an opportunistic infection with a diverse spectrum of clinical presentations and severity. With early diagnosis, specific antimicrobial therapy, and state-of-the-art intensive care, mortality can be reduced to less than 10%. However, mortality remains much higher in the many endemic regions where health resources remain scarce. Genotyping of B pseudomallei informs evolving local and global epidemiology.

FUNDING

The Australian National Health and Medical Research Council.

Genome-Wide Metabolic Reconstruction of the Synthesis of Polyhydroxyalkanoates from Sugars and Fatty Acids by Burkholderia Sensu Lato Species

Burkholderia sensu lato (s.l.) species have a versatile metabolism. The aims of this review are the genomic reconstruction of the metabolic pathways involved in the synthesis of polyhydroxyalkanoates (PHAs) by Burkholderia s.l. genera, and the characterization of the PHA synthases and the pha genes organization. The reports of the PHA synthesis from different substrates by Burkholderia s.l. strains were reviewed. Genome-guided metabolic reconstruction involving the conversion of sugars and fatty acids into PHAs by 37 Burkholderia s.l. species was performed. Sugars are metabolized via the Entner-Doudoroff (ED), pentose-phosphate (PP), and lower Embden-Meyerhoff-Parnas (EMP) pathways, which produce reducing power through NAD(P)H synthesis and PHA precursors. Fatty acid substrates are metabolized via β-oxidation and de novo synthesis of fatty acids into PHAs. The analysis of 194 Burkholderia s.l. genomes revealed that all strains have the phaC, phaA, and phaB genes for PHA synthesis, wherein the phaC gene is generally present in ≥2 copies. PHA synthases were classified into four phylogenetic groups belonging to class I II and III PHA synthases and one outlier group. The reconstruction of PHAs synthesis revealed a high level of gene redundancy probably reflecting complex regulatory layers that provide fine tuning according to diverse substrates and physiological conditions.

Discovery of a Novel Lineage Burkholderia cepacia ST 1870 Endophytically Isolated from Medicinal Polygala paniculata Which Shows Potent In Vitro Antileishmanial and Antimicrobial Effects

In this study, we report the isolation and identification of an endophytic strain of Burkholderia cepacia (COPS strain) associated with Polygala paniculata roots. Polygala plants are rich sources of promising microbiomes, of which the literature reports several pharmacological effects, such as trypanocidal, antinociceptive, anesthetic, anxiolytics, and anticonvulsant activities. B. cepacia COPS belongs to a new sequence type (ST 1870) and harbors a genome estimated in 8.3 Mbp which exhibits the aminoglycosides and beta-lactams resistance genes aph(3′)-IIa and bla_TEM-116, respectively. Analysis performed using MLST, average nucleotide identity, and digital DNA-DNA hybridization support its species-level identification and reveals its novel housekeeping genes alleles gyrB, lepA, and phaC. The root endophyte B. cepacia COPS drew our attention from a group of 14 bacterial isolates during the primary screening for being potentially active against Staphylococcus aureus ATCC 29213, Enterococcus faecalis ATCC 29212, Micrococcus luteus ATCC 9341, Escherichia coli ATCC 25922, and Candida albicans ATCC 10231 and exhibited the broad-spectrum activity against phytopathogenic fungi. In addition, COPS strain showed production of protease, lipase, and esterase in solid media, and its natural product extract showed potent inhibition against fungal plant pathogens, such as Moniliophthora perniciosa, whose antagonism index (89.32%) exceeded the positive control (74.17%), whereas Sclerotinia sclerotiorum and Ceratocystis paradoxa showed high percentages of inhibition (85.53% and 82.69%, respectively). COPS crude extract also significantly inhibited S. epidermidis ATCC 35984, E. faecium ATCC 700221 (MIC values of 32 μg/mL for both), E. faecalis ATCC 29212 (64 μg/mL), and S. aureus ATCC 25923 (128 μg/mL). We observed moderate antagonistic activity against A. baumannii ATCC 19606 and E. coli ATCC 25922 (both at 512 μg/mL), as well as potent cytotoxic effects on Leishmania infantum and Leishmania major promastigote forms with 78.25% and 57.30% inhibition. In conclusion, this study presents for the first time the isolation of an endophytic B. cepacia strain associated with P. paniculata and enough evidence that these plants may be considered a rich source of microbes for the fight against neglected diseases.

Paraburkholderia acidiphila sp. nov., Paraburkholderia acidisoli sp. nov. and Burkholderia guangdongensis sp. nov., isolated from forest soil, and reclassification of Burkholderia ultramafica as Paraburkholderia ultramafica comb. nov

Three Gram-stain-negative, aerobic, motile and rod-shaped bacterial strains, 7Q-K02^T, DHF22^T and DHOM02^T, were isolated from forest soil sampled at Dinghushan Biosphere Reserve, Guangdong Province, China. Strains 7Q-K02^T, DHF22^T and DHOM02^T grew at 4-37, 4-42 and 12-37 °C, pH 3.0-8.5, 3.5-8.5 and 5.0-8.0, and in the presence of 0-3.0, 0-3.5 and 0-2.5 % (w/v) NaCl; with optima at 28-33, 28 and 28-33 °C, pH 3.5-6.5, 4.0-5.5 and 6.5-7.0, and 0-1.5, 0-1.5 and 0.5-1.5 % (w/v) NaCl, respectively. Strains 7Q-K02^T and DHF22^T have the highest 16S rRNA gene sequence similarities of 99.0 and 98.0 % to Paraburkholderia sacchari LMG 19450^T and 97.7 % between themselves, while strain DHOM02^T shares the highest similarity of 98.4 % to 'Burkholderia rinojensis' A396^T followed by 98.3 % to Burkholderia plantarii ATCC 43733^T. In the 16S rRNA gene sequence phylogram, strain 7Q-K02^T formed a sister branch with Paraburkholderia sacchari, Paraburkholderia oxyphila and Paraburkholderia paradisi, and strain DHF22^T was separated from all other species within the genus Paraburkholderia, while strain DHOM02^T formed a separated clade with members of the genus Burkholderia. The DNA G+C contents of strains 7Q-K02^T, DHF22^T and DHOM02^T wwe 64.3, 65.4 and 66.6 %, respectively. Digital DNA-DNA hybridization (dDDH) and average nucleotide identity (ANI) values of strains 7Q-K02^T, DHF22^T and closely related Paraburkholderia strains were in the ranges of 25.5-43.7 % and 81.5-91.3 %, respectively. While dDDH and ANI values between strain DHOM02^T and Burkholderia strains with genome sequence data were in the ranges of 22.4-31.0 % and 78.2-86.1 %, respectively. These three strains have the same major respiratory quinone: ubiquinone-8. Strains 7Q-K02^T, DHF22^T and DHOM02^T have C_16 : 0, C_17 : 0 cyclo, C_19 : 0 cyclo ω8c and summed feature 8 (C_18 : 1  ω7c/C_18 : 1  ω6c) as their major fatty acid compositions. The major polar lipids were phosphatidylethanolamine, phosphatidylglycerol and diphosphatidylglycerol. On the basis of phenotypic, phylogenetic, genomic analyses and chemotaxonomic data, strains 7Q-K02^T and DHF22^T represent two novel species of the genus Paraburkholderia, for which the names Paraburkholderia acidiphila sp. nov. (type strain 7Q-K02^T=CGMCC 1.15433^T=KCTC 62472^T=LMG 29209^T) and Paraburkholderia acidisoli sp. nov. (type strain DHF22^T=GDMCC 1.1448^T=LMG 30262^T) are proposed, while strain DHOM02^T represents a novel species in the genus Burkholderia, for which the name Burkholderia guangdongensis sp. nov. (type strain DHOM02^T=KCTC 42625^T=LMG 28843^T) is proposed. We also propose to transfer Burkholderia ultramafica to the genus Paraburkholderia as Paraburkholderia ultramafica comb. nov. based mainly on the results of phylogenomic analysis.

Mapping of the Denitrification Pathway in Burkholderia thailandensis by Genome-Wide Mutant Profiling

Burkholderia thailandensis is a soil saprophyte that is closely related to the pathogen Burkholderia pseudomallei, the etiological agent of melioidosis in humans. The environmental niches and infection sites occupied by these bacteria are thought to contain only limited concentrations of oxygen, where they can generate energy via denitrification. However, knowledge of the underlying molecular basis of the denitrification pathway in these bacteria is scarce. In this study, we employed a transposon sequencing (Tn-Seq) approach to identify genes conferring a fitness benefit for anaerobic growth of B. thailandensis Of the 180 determinants identified, several genes were shown to be required for growth under denitrifying conditions: the nitrate reductase operon narIJHGK2K1, the aniA gene encoding a previously unknown nitrite reductase, and the petABC genes encoding a cytochrome bc ₁, as well as three novel regulators that control denitrification. Our Tn-Seq data allowed us to reconstruct the entire denitrification pathway of B. thailandensis and shed light on its regulation. Analyses of growth behaviors combined with measurements of denitrification metabolites of various mutants revealed that nitrate reduction provides sufficient energy for anaerobic growth, an important finding in light of the fact that some pathogenic Burkholderia species can use nitrate as a terminal electron acceptor but are unable to complete denitrification. Finally, we demonstrated that a nitrous oxide reductase mutant is not affected for anaerobic growth but is defective in biofilm formation and accumulates N₂O, which may play a role in the dispersal of B. thailandensis biofilms.IMPORTANCE Burkholderia thailandensis is a soil-dwelling saprophyte that is often used as surrogate of the closely related pathogen Burkholderia pseudomallei, the causative agent of melioidosis and a classified biowarfare agent. Both organisms are adapted to grow under oxygen-limited conditions in rice fields by generating energy through denitrification. Microoxic growth of B. pseudomallei is also considered essential for human infections. Here, we have used a Tn-Seq approach to identify the genes encoding the enzymes and regulators required for growth under denitrifying conditions. We show that a mutant that is defective in the conversion of N₂O to N₂, the last step in the denitrification process, is unaffected in microoxic growth but is severely impaired in biofilm formation, suggesting that N₂O may play a role in biofilm dispersal. Our study identified novel targets for the development of therapeutic agents to treat meliodiosis.

Molecular Characteristics of Burkholderia pseudomallei Collected From Humans in Hainan, China

Melioidosis is a common infectious disease in Southeast Asia and Northern Australia. In Hainan, several cases have been reported, but no systematic study has yet been done on the molecular epidemiology profiles of the organism. An investigation of the molecular epidemiology links and population structure of Burkholderia pseudomallei would help to better understand the clonally of the isolates and differences among them. In this study, multilocus variable-number tandem repeat analysis (MLVA), and multilocus sequence typing (MLST) were applied to examine the epidemiological relatedness and population structure of 166 B. pseudomallei isolates obtained during 2002–2014 in Hainan, China. Both the MLVA_4 and MLST approaches had high discriminatory power for this population, with diversity indices of 0.9899 and 0.9457, respectively. However, the MLVA_4 assay showed a higher discriminatory power than the MLST approach, and a variable-number tandem repeat (VNTR3 933) found by the MLVA approach was the most useful in discriminating strains from this province. A total of 166 strains yielded 99 MLVA_4 genotypes, of which 34 genotypes were shared by 101 isolates, for a clustering rate of 60.8% (101/166), which suggested that some cases may have a common source. Additionally, 65 isolates showed distinct genotypes, indicating that more than 39.2% (65/166) of melioidosis cases in Hainan had epidemiologically unrelated or sporadic characteristics. The 166 isolates were resolved into 48 STs, of which five STs (ST55, -70, -46, -50, and -58) were here found to be predominant. Phylogenetic analysis of 116 isolates conducted using the eBURST v3 segregated the 48 STs into eight groups with ST50 as predicted founder, and 21 STs were found to be singletons, which suggest that the strains in the Hainan region represent a high diversity of ST clones, indicating that many B. pseudomallei clone groups are endemic to this region. Moreover, ST50 had 5 SLV, 7 DLV, 6 TLV, and 29 satellite STs and formed a radial expansion pattern, suggesting that the melioidosis epidemic in this study was mainly caused by the clonal expansion of ST 50. Phylogenetic analysis on global scale suggests that China’s isolates are closely related to isolates from Southeast Asia, particularly from Thailand and Malaysia.

Human Melioidosis

The causative agent of melioidosis, Burkholderia pseudomallei, a tier 1 select agent, is endemic in Southeast Asia and northern Australia, with increased incidence associated with high levels of rainfall. Increasing reports of this condition have occurred worldwide, with estimates of up to 165,000 cases and 89,000 deaths per year. The ecological niche of the organism has yet to be clearly defined, although the organism is associated with soil and water. The culture of appropriate clinical material remains the mainstay of laboratory diagnosis. Identification is best done by phenotypic methods, although mass spectrometric methods have been described. Serology has a limited diagnostic role. Direct molecular and antigen detection methods have limited availability and sensitivity. Clinical presentations of melioidosis range from acute bacteremic pneumonia to disseminated visceral abscesses and localized infections. Transmission is by direct inoculation, inhalation, or ingestion. Risk factors for melioidosis include male sex, diabetes mellitus, alcohol abuse, and immunosuppression. The organism is well adapted to intracellular survival, with numerous virulence mechanisms. Immunity likely requires innate and adaptive responses. The principles of management of this condition are drainage and debridement of infected material and appropriate antimicrobial therapy. Global mortality rates vary between 9% and 70%. Research into vaccine development is ongoing.

Melioidosis fatalities in captive slender-tailed meerkats (Suricata suricatta): combining epidemiology, pathology and whole-genome sequencing supports variable mechanisms of transmission with one health implications

Background

Melioidosis is a tropical infectious disease which is being increasingly recognised throughout the globe. Infection occurs in humans and animals, typically through direct exposure to soil or water containing the environmental bacterium Burkholderia pseudomallei. Case clusters of melioidosis have been described in humans following severe weather events and in exotic animals imported into melioidosis endemic zones. Direct transmission of B. pseudomallei between animals and/or humans has been documented but is considered extremely rare. Between March 2015 and October 2016 eight fatal cases of melioidosis were reported in slender-tailed meerkats (Suricata suricatta) on display at a Wildlife Park in Northern Australia. To further investigate the melioidosis case cluster we sampled the meerkat enclosure and adjacent park areas and performed whole-genome sequencing (WGS) on all culture-positive B. pseudomallei environmental and clinical isolates.

Results

WGS confirmed that the fatalities were caused by two different B. pseudomallei sequence types (STs) but that seven of the meerkat isolates were highly similar on the whole-genome level. Used concurrently with detailed pathology data, our results demonstrate that the seven cases originated from a single original source, but routes of infection varied amongst meerkats belonging to the clonal outbreak cluster. Moreover, in some instances direct transmission may have transpired through wounds inflicted while fighting.

Conclusions

Collectively, this study supports the use of high-resolution WGS to enhance epidemiological investigations into transmission modalities and pathogenesis of melioidosis, especially in the instance of a possible clonal outbreak scenario in exotic zoological collections. Such findings from an animal outbreak have important One Health implications.

Application of qPCR assays based on haloacids transporter gene dehp2 for discrimination of Burkholderia and Paraburkholderia

Background

A major facilitator superfamily transporter Dehp2 was recently shown to be playing an important role in transport and biodegradation of haloacids in Paraburkholderia caribensis MBA4, and Dehp2 is phylogenetically conserved in Burkholderia sensu lato.

Results

We designed both Burkholderia sensu stricto-specific and Paraburkholderia-specific qPCR assays based on dehp2 and 16S rRNA, and validated the qPCR assays in 12 bacterial strains. The qPCR assays could detect single species of Burkholderia sensu stricto or Paraburkholderia with high sensitivity and discriminate them in mixtures with high specificity over a wide dynamic range of relative concentrations. At relatively lower cost compared with sequencing-based approach, the qPCR assays will facilitate discrimination of Burkholderia sensu stricto and Paraburkholderia in a large number of samples.

Conclusions

For the first time, we report the utilization of a haloacids transporter gene for discriminative purpose in Burkholderia sensu lato. This enables not only quick decision on proper handling of putative pathogenic samples in Burkholderia sensu stricto group but also future exploitation of relevant species in Paraburkholderia group for haloacids biodegradation purposes.

Electronic supplementary material

The online version of this article (10.1186/s12866-019-1411-0) contains supplementary material, which is available to authorized users.

Open-access bacterial population genomics: BIGSdb software, the PubMLST.org website and their applications

The PubMLST.org website hosts a collection of open-access, curated databases that integrate population sequence data with provenance and phenotype information for over 100 different microbial species and genera.  Although the PubMLST website was conceived as part of the development of the first multi-locus sequence typing (MLST) scheme in 1998 the software it uses, the Bacterial Isolate Genome Sequence database (BIGSdb, published in 2010), enables PubMLST to include all levels of sequence data, from single gene sequences up to and including complete, finished genomes.  Here we describe developments in the BIGSdb software made from publication to June 2018 and show how the platform realises microbial population genomics for a wide range of applications.  The system is based on the gene-by-gene analysis of microbial genomes, with each deposited sequence annotated and curated to identify the genes present and systematically catalogue their variation.  Originally intended as a means of characterising isolates with typing schemes, the synthesis of sequences and records of genetic variation with provenance and phenotype data permits highly scalable (whole genome sequence data for tens of thousands of isolates) means of addressing a wide range of functional questions, including: the prediction of antimicrobial resistance; likely cross-reactivity with vaccine antigens; and the functional activities of different variants that lead to key phenotypes.  There are no limitations to the number of sequences, genetic loci, allelic variants or schemes (combinations of loci) that can be included, enabling each database to represent an expanding catalogue of the genetic variation of the population in question.  In addition to providing web-accessible analyses and links to third-party analysis and visualisation tools, the BIGSdb software includes a RESTful application programming interface (API) that enables access to all the underlying data for third-party applications and data analysis pipelines.

Profiling microbial strains in urban environments using metagenomic sequencing data

Background

The microbial communities populating human and natural environments have been extensively characterized with shotgun metagenomics, which provides an in-depth representation of the microbial diversity within a sample. Microbes thriving in urban environments may be crucially important for human health, but have received less attention than those of other environments. Ongoing efforts started to target urban microbiomes at a large scale, but the most recent computational methods to profile these metagenomes have never been applied in this context. It is thus currently unclear whether such methods, that have proven successful at distinguishing even closely related strains in human microbiomes, are also effective in urban settings for tasks such as cultivation-free pathogen detection and microbial surveillance. Here, we aimed at a) testing the currently available metagenomic profiling tools on urban metagenomics; b) characterizing the organisms in urban environment at the resolution of single strain and c) discussing the biological insights that can be inferred from such methods.

Results

We applied three complementary methods on the 1614 metagenomes of the CAMDA 2017 challenge. With MetaMLST we identified 121 known sequence-types from 15 species of clinical relevance. For instance, we identified several Acinetobacter strains that were close to the nosocomial opportunistic pathogen A. nosocomialis. With StrainPhlAn, a generalized version of the MetaMLST approach, we inferred the phylogenetic structure of Pseudomonas stutzeri strains and suggested that the strain-level heterogeneity in environmental samples is higher than in the human microbiome. Finally, we also probed the functional potential of the different strains with PanPhlAn. We further showed that SNV-based and pangenome-based profiling provide complementary information that can be combined to investigate the evolutionary trajectories of microbes and to identify specific genetic determinants of virulence and antibiotic resistances within closely related strains.

Conclusion

We show that strain-level methods developed primarily for the analysis of human microbiomes can be effective for city-associated microbiomes. In fact, (opportunistic) pathogens can be tracked and monitored across many hundreds of urban metagenomes. However, while more effort is needed to profile strains of currently uncharacterized species, this work poses the basis for high-resolution analyses of microbiomes sampled in city and mass transportation environments.

Reviewers

This article was reviewed by Alexandra Bettina Graf, Daniel Huson and Trevor Cickovski.

Electronic supplementary material

The online version of this article (10.1186/s13062-018-0211-z) contains supplementary material, which is available to authorized users.

Phylogeography of Burkholderia pseudomallei Isolates, Western Hemisphere

The bacterium Burkholderia pseudomallei causes melioidosis, which is mainly associated with tropical areas. We analyzed single-nucleotide polymorphisms (SNPs) among genome sequences from isolates of B. pseudomallei that originated in the Western Hemisphere by comparing them with genome sequences of isolates that originated in the Eastern Hemisphere. Analysis indicated that isolates from the Western Hemisphere form a distinct clade, which supports the hypothesis that these isolates were derived from a constricted seeding event from Africa. Subclades have been resolved that are associated with specific regions within the Western Hemisphere and suggest that isolates might be correlated geographically with cases of melioidosis. One isolate associated with a former World War II prisoner of war was believed to represent illness 62 years after exposure in Southeast Asia. However, analysis suggested the isolate originated in Central or South America.

Comparative Genomics of Burkholderia singularis sp. nov., a Low G+C Content, Free-Living Bacterium That Defies Taxonomic Dissection of the Genus Burkholderia

Four Burkholderia pseudomallei-like isolates of human clinical origin were examined by a polyphasic taxonomic approach that included comparative whole genome analyses. The results demonstrated that these isolates represent a rare and unusual, novel Burkholderia species for which we propose the name B. singularis. The type strain is LMG 28154T (=CCUG 65685T). Its genome sequence has an average mol% G+C content of 64.34%, which is considerably lower than that of other Burkholderia species. The reduced G+C content of strain LMG 28154T was characterized by a genome wide AT bias that was not due to reduced GC-biased gene conversion or reductive genome evolution, but might have been caused by an altered DNA base excision repair pathway. B. singularis can be differentiated from other Burkholderia species by multilocus sequence analysis, MALDI-TOF mass spectrometry and a distinctive biochemical profile that includes the absence of nitrate reduction, a mucoid appearance on Columbia sheep blood agar, and a slowly positive oxidase reaction. Comparisons with publicly available whole genome sequences demonstrated that strain TSV85, an Australian water isolate, also represents the same species and therefore, to date, B. singularis has been recovered from human or environmental samples on three continents.

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.

A Quadruplex Real-Time PCR Assay for the Rapid Detection and Differentiation of the Most Relevant Members of the B. pseudomallei Complex: B. mallei, B. pseudomallei, and B. thailandensis

The Burkholderia pseudomallei complex classically consisted of B. mallei, B. pseudomallei, and B. thailandensis, but has now expanded to include B. oklahomensis, B. humptydooensis, and three unassigned Burkholderia clades. Methods for detecting and differentiating the B. pseudomallei complex has been the topic of recent research due to phenotypic and genotypic similarities of these species. B. mallei and B. pseudomallei are recognized as CDC Tier 1 select agents, and are the causative agents of glanders and melioidosis, respectively. Although B. thailandensis and B. oklahomensis are generally avirulent, both display similar phenotypic characteristics to that of B. pseudomallei. B. humptydooensis and the Burkholderia clades are genetically similar to the B. pseudomallei complex, and are not associated with disease. Optimal identification of these species remains problematic, and PCR-based methods can resolve issues with B. pseudomallei complex detection and differentiation. Currently, no PCR assay is available that detects the major species of the B. pseudomallei complex. A real-time PCR assay in a multiplex single-tube format was developed to simultaneously detect and differentiate B. mallei, B. pseudomallei, and B. thailandensis, and a common sequence found in B. pseudomallei, B. mallei, B. thailandensis, and B. oklahomensis. A total of 309 Burkholderia isolates and 5 other bacterial species were evaluated. The assay was 100% sensitive and specific, demonstrated sensitivity beyond culture and GC methods for the isolates tested, and is completed in about an hour with a detection limit between 2.6pg and 48.9pg of gDNA. Bioinformatic analyses also showed the assay is likely 100% specific and sensitive for all 84 fully sequenced B. pseudomallei, B. mallei, B. thailandensis, and B. oklahomensis strains currently available in GenBank. For these reasons, this assay could be a rapid and sensitive tool in the detection and differentiation for those species of the B. pseudomallei complex with recognized clinical and practical significance.