Search

Abstract Background Antimicrobial resistance (AMR) is an intensifying threat that requires urgent mitigation to avoid a post-antibiotic era. Pseudomonas aeruginosa represents one of the greatest AMR concerns due to increasing multi- and pan-drug resistance rates. Shotgun sequencing is gaining traction for in silico AMR profiling due to its unambiguity and transferability; however, accurate and comprehensive AMR prediction from P. aeruginosa genomes remains an unsolved problem. Methods We first curated the most comprehensive database yet of known P. aeruginosa AMR variants. Next, we performed comparative genomics and microbial genome-wide association study analysis across a Global isolate Dataset (n = 1877) with paired antimicrobial phenotype and genomic data to identify novel AMR variants. Finally, the performance of our P. aeruginosa AMR database, implemented in our AMR detection and prediction tool, ARDaP, was compared with three previously published in silico AMR gene detection or phenotype prediction tools—abritAMR, AMRFinderPlus, ResFinder—across both the Global Dataset and an analysis-naïve Validation Dataset (n = 102). Results Our AMR database comprises 3639 mobile AMR genes and 728 chromosomal variants, including 75 previously unreported chromosomal AMR variants, 10 variants associated with unusual antimicrobial susceptibility, and 281 chromosomal variants that we show are unlikely to confer AMR. Our pipeline achieved a genotype-phenotype balanced accuracy (bACC) of 85% and 81% across 10 clinically relevant antibiotics when tested against the Global and Validation Datasets, respectively, vs. just 56% and 54% with abritAMR, 58% and 54% with AMRFinderPlus, and 60% and 53% with ResFinder. ARDaP’s superior performance was predominantly due to the inclusion of chromosomal AMR variants, which are generally not identified with most AMR identification tools. Conclusions Our ARDaP software and associated AMR variant database provides an accurate tool for predicting AMR phenotypes in P. aeruginosa, far surpassing the performance of current tools. Implementation of ARDaP for routine AMR prediction from P. aeruginosa genomes and metagenomes will improve AMR identification, addressing a critical facet in combatting this treatment-refractory pathogen. However, knowledge gaps remain in our understanding of the P. aeruginosa resistome, particularly the basis of colistin AMR.

Historically, Boswellia resin, commonly known as frankincense, has been used for a variety of applications, including medicinal purposes and in religious ceremonies. These practices continued into the modern era, in which the trade of resin has increased dramatically. The genus Boswellia contains many species that are not currently protected by the Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES), but research has found that some of these populations are predicted to decrease by 50 % in the next two decades due to overexploitation (Bongers et al., 2019). This highlights the need for a method that can identify Boswellia resins. In this research, resin from nine different species of Boswellia were collected and analyzed in positive ion mode by Direct Analysis in Real Time Time-of-Flight mass spectrometry (DART TOFMS) by dissolving the resin in methanol to obtain a chemotype devoid of background noise. Chemometric analysis of the data indicated that positive ion spectra from the resins separated the species effectively, indicating that DART TOFMS can be used to identify resins in trade.

Because of the rich omega-3 fatty acids content, harp seal (Pagophilus groenlandicus) oil is a popular supplement that is packaged as pills in Canada and sold for medicinal purposes, although this practice is banned in the United States. Due to US regulations, it is important to be able to distinguish between fish oil and seal oil, but the taxonomic determination of oils provenance has been a difficult problem to solve. In this study, Direct Analysis in Real Time time-of-flight mass spectrometry (DART TOFMS) was used to analyze the chemotypes of blubber samples collected from seven species of marine mammals, including seals, sea lions, and a porpoise. Results indicated that the chemotype profiles found in negative-ion mode could be used to separate all of the species using Discriminant Analysis of Principal Components (DAPC). Consequently, this study suggests that it may be possible to identify the taxonomic source of marine mammal oils based on chemical chemotypes.

The Convention on International Trade in Endangered Species of Wild Fauna and Flora (CITES) has listed over 100 species of trees whose trade is primarily in timber, and whose existence may be endangered due to either indiscriminate or illegal logging. Species determination of timber is essential for enforcing logging quotas and species protections laws, but it is a challenging endeavor. Here we present a technique that relies on mass spectra chemotypes for species assignment and demonstrate that it provides high reliability.The reliability of timber identification by mass spectrometry relies on multiple factors: 1) access to a robust wood spectra database (e.g., ForeST ©: Forensic Spectra of Trees) derived from a curated xylarium, 2) analysis of the correct wood tissue (i.e., heartwood vs sapwood), 3) reproducibility of the spectra, 4) repeatability of the analysis and 5) skill in interpreting the data. The goal of timber identification is to assign wood products to species, but in some cases the precision of the analysis will be to the genus or even family level. This may be due, for example, to shared chemotypes among closely related tree species. These higher taxa assignments do not represent an error (false positive or false negative), but rather reflect the complexity of wood chemistry in relation to phylogeny. Such higher taxa determinations are often sufficient for enforcement purposes, as many timber trees are protected at the genus or above.

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.

In Pterocarpus, three out of 46 tree species (P. erinaceus, P. santalinus, and P. tinctorius) are listed as endangered and protected under Appendix II of CITES, and the status of several Pterocarpus species suggests that more will find a place on CITES in the future. As unsustainable forest exploitation has increased, regulations for timber traceability have also increased with the creation of such laws as the 2005 Lacey Act in the United States (Lacey Act, 2005) and the European Union Timber Regulation (EUTR) in 2013. These laws were implemented to ensure that members of the timber trade operated in accordance with guidelines and practiced due diligence, with the goal of halting illegal timber imports (FAO, 2016). Unfortunately, the illegal timber trade remains prevalent and has progressed into the third largest transnational crime (Mavrellis, 2017). To combat the prevalence of illegal timber, researchers have utilized Direct Analysis in Real Time Time-of-Flight Mass Spectrometry (DART-TOFMS). With the success found using DART-TOFMS for species identification, we investigated whether the method could be used to identify the original geographic region of timber through analysis of 211 samples taken from 12 Pterocarpus species taken from the Neotropics, Africa, South Asia, and Southeast Asia. We found that analysis using DART-TOFMS could identify the broad geographical provenances of the timber with an accuracy of 98%. This research demonstrated that DART-TOFMS can be used to detect the phytochemical trends that exist between timber species occupying broad geographic regions.

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 (Antimicrobial Resistance Detection and Prediction) 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).

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

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.

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.

The Tier 1 select agent Burkholderia pseudomallei is an environmental bacterium that causes melioidosis, a high mortality disease. Variably present genetic markers used to elucidate strain origin, relatedness and virulence in B. pseudomallei include the Burkholderia intracellular motility factor A (bimA) and filamentous hemagglutinin 3 (fhaB3) gene variants. Three lipopolysaccharide (LPS) O-antigen types in B. pseudomallei have 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 clinical B. pseudomallei isolates 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; p

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 pediatric airways that deserves further exploration. Understanding the pathoadaptative mechanisms that NTHi employs to infect and persist in the lower pediatric airways is essential for devising targeted diagnostics and treatments aimed at minimizing disease severity, and ultimately, preventing NTHi lung infections and subsequent chronic lung disease in children.

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

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.

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.

ABSTRACT Cystic fibrosis (CF) is a genetic disorder characterized by progressive lung function decline. CF patients are at an increased risk of respiratory infections, including those by the environmental bacterium Burkholderia pseudomallei, the causative agent of melioidosis. Here, we compared the genomes of B. pseudomallei isolates collected between ~4 and 55 months apart from seven chronically infected CF patients. Overall, the B. pseudomallei strains showed evolutionary patterns similar to those of other chronic infections, including emergence of antibiotic resistance, genome reduction, and deleterious mutations in genes involved in virulence, metabolism, environmental survival, and cell wall components. We documented the first reported B. pseudomallei hypermutators, which were likely caused by defective MutS. Further, our study identified both known and novel molecular mechanisms conferring resistance to three of the five clinically important antibiotics for melioidosis treatment. Our report highlights the exquisite adaptability of microorganisms to long-term persistence in their environment and the ongoing challenges of antibiotic treatment in eradicating pathogens in the CF lung. Convergent evolution with other CF pathogens hints at a degree of predictability in bacterial evolution in the CF lung and potential targeted eradication of chronic CF infections in the future. IMPORTANCE Burkholderia pseudomallei, the causative agent of melioidosis, is an environmental opportunistic bacterium that typically infects immunocompromised people and those with certain risk factors such as cystic fibrosis (CF). Patients with CF tend to develop chronic melioidosis infections, for reasons that are not well understood. This report is the first to describe B. pseudomallei evolution within the CF lung during chronic infection. We show that the pathways by which B. pseudomallei adapts to the CF lung are similar to those seen in better-studied CF pathogens such as Pseudomonas aeruginosa, Staphylococcus aureus, and Burkholderia cepacia complex species. Adaptations include the accumulation of antibiotic resistance, loss of nonessential genes, metabolic alterations, and virulence factor attenuation. Known and novel mechanisms of resistance to three of the five antibiotics used in melioidosis treatment were identified. Similar pathways of evolution in CF pathogens, including B. pseudomallei, provide exciting avenues for more-targeted treatment of chronic, recalcitrant infections.

The frequency with which melioidosis results from inhalation rather than percutaneous inoculation or ingestion is unknown. We recovered Burkholderia pseudomallei from air samples at the residence of a patient with presumptive inhalational melioidosis and used whole-genome sequencing to link the environmental bacteria to B. pseudomallei recovered from the patient.

After heavy rains and flooding during early 2011 in the normally arid interior of Australia, melioidosis was diagnosed in 6 persons over a 4-month period. Although the precise global distribution of the causal bacterium Burkholderia pseudomallei remains to be determined, this organism can clearly survive in harsh and even desert environments outside the wet tropics.

ABSTRACT Whole-genome sequence (WGS) data are commonly used to design diagnostic targets for the identification of bacterial pathogens. To do this effectively, genomics databases must be comprehensive to identify the strict core genome that is specific to the target pathogen. As additional genomes are analyzed, the core genome size is reduced and there is erosion of the target-specific regions due to commonality with related species, potentially resulting in the identification of false positives and/or false negatives. IMPORTANCE A comparative analysis of 1,130 Burkholderia genomes identified unique markers for many named species, including the human pathogens B. pseudomallei and B. mallei. Due to core genome reduction and signature erosion, only 38 targets specific to B. pseudomallei/mallei were identified. By using only public genomes, a larger number of markers were identified, due to undersampling, and this larger number represents the potential for false positives. This analysis has implications for the design of diagnostics for other species where the genomic space of the target and/or closely related species is not well defined.

ABSTRACT Burkholderia pseudomallei, an environmental bacterium that causes the deadly disease melioidosis, is endemic in northern Australia and Southeast Asia. An increasing number of melioidosis cases are being reported in other tropical regions, including Africa and the Indian Ocean islands. B. pseudomallei first emerged in Australia, with subsequent rare dissemination event(s) to Southeast Asia; however, its dispersal to other regions is not yet well understood. We used large-scale comparative genomics to investigate the origins of three B. pseudomallei isolates from Madagascar and two from Burkina Faso. Phylogenomic reconstruction demonstrates that these African B. pseudomallei isolates group into a single novel clade that resides within the more ancestral Asian clade. Intriguingly, South American strains reside within the African clade, suggesting more recent dissemination from West Africa to the Americas. Anthropogenic factors likely assisted in B. pseudomallei dissemination to Africa, possibly during migration of the Austronesian peoples from Indonesian Borneo to Madagascar ~2,000 years ago, with subsequent genetic diversity driven by mutation and recombination. Our study provides new insights into global patterns of B. pseudomallei dissemination and adds to the growing body of evidence of melioidosis endemicity in Africa. Our findings have important implications for melioidosis diagnosis and management in Africa. IMPORTANCE Sporadic melioidosis cases have been reported in the African mainland and Indian Ocean islands, but until recently, these regions were not considered areas where B. pseudomallei is endemic. Given the high mortality rate of melioidosis, it is crucial that this disease be recognized and suspected in all regions of endemicity. Previous work has shown that B. pseudomallei originated in Australia, with subsequent introduction into Asia; however, the precise origin of B. pseudomallei in other tropical regions remains poorly understood. Using whole-genome sequencing, we characterized B. pseudomallei isolates from Madagascar and Burkina Faso. Next, we compared these strains to a global collection of B. pseudomallei isolates to identify their evolutionary origins. We found that African B. pseudomallei strains likely originated from Asia and were closely related to South American strains, reflecting a relatively recent shared evolutionary history. We also identified substantial genetic diversity among African strains, suggesting long-term B. pseudomallei endemicity in this region.

In December 2009, two unusual cases of anthrax were diagnosed in heroin users in Scotland. A subsequent anthrax outbreak in heroin users emerged throughout Scotland and expanded into England and Germany, sparking concern of nefarious introduction of anthrax spores into the heroin supply. To better understand the outbreak origin, we used established genetic signatures that provided insights about strain origin. Next, we sequenced the whole genome of a representative Bacillus anthracis strain from a heroin user (Ba4599), developed Ba4599-specific single-nucleotide polymorphism assays, and genotyped all available material from other heroin users with anthrax. Of 34 case-patients with B. anthracis–positive PCR results, all shared the Ba4599 single-nucleotide polymorphism genotype. Phylogeographic analysis demonstrated that Ba4599 was closely related to strains from Turkey and not to previously identified isolates from Scotland or Afghanistan, the presumed origin of the heroin. Our results suggest accidental contamination along the drug trafficking route through a cutting agent or animal hides used to smuggle heroin into Europe.

Melioidosis, caused by the bacterium Burkholderia pseudomallei, is an infectious disease of humans or animals, and the specific environmental conditions that are present in western Indian Ocean islands are particularly suitable for the establishment/survival of B. pseudomallei. Indeed, an increasing number of new cases have been reported in this region (Madagascar, Mauritius, Réunion (France), and Seychelles, except Comoros and Mayotte (France)), and are described in this review. Our review clearly points out that further studies are needed in order to investigate the real incidence and burden of melioidosis in the western Indian Ocean and especially Madagascar, since it is likely to be higher than currently reported. Thus, research and surveillance priorities were recommended (i) to improve awareness of melioidosis in the population and among clinicians; (ii) to improve diagnostics, in order to provide rapid and effective treatment; (iii) to implement a surveillance and reporting system in the western Indian Ocean; and (iv) to investigate the presence of B. pseudomallei in environmental samples, since we have demonstrated its presence in soil samples originating from the yard of a Madagascan case.

We collected epidemiologic and molecular data from Burkholderia mallei isolates from equines in Punjab, Pakistan from 1999 through 2007. We show that recent outbreaks are genetically distinct from available whole genome sequences and that these genotypes are persistent and ubiquitous in Punjab, probably due to human-mediated movement of equines.

Burkholderia pseudomallei is a Gram-negative environmental bacterium that causes melioidosis, a potentially life-threatening infectious disease affecting mammals, including humans. Melioidosis symptoms are both protean and diverse, ranging from mild, localized skin infections to more severe and often fatal presentations including pneumonia, septic shock with multiple internal abscesses and occasionally neurological involvement. Several ubiquitous virulence determinants in B. pseudomallei have already been discovered. However, the molecular basis for differential pathogenesis has, until now, remained elusive. Using clinical data from 556 Australian melioidosis cases spanning more than 20 years, we identified a Burkholderia mallei-like actin polymerization bimA(Bm) gene that is strongly associated with neurological disease. We also report that a filamentous hemagglutinin gene, fhaB3, is associated with positive blood cultures but is negatively correlated with localized skin lesions without sepsis. We show, for the first time, that variably present virulence factors play an important role in the pathogenesis of melioidosis. Collectively, our study provides a framework for assessing other non-ubiquitous bacterial virulence factors and their association with disease, such as candidate loci identified from large-scale microbial genome-wide association studies.

ABSTRACT Burkholderia pseudomallei causes the potentially fatal disease melioidosis. It is generally accepted that B. pseudomallei is a noncommensal bacterium and that any culture-positive clinical specimen denotes disease requiring treatment. Over a 23-year study of melioidosis cases in Darwin, Australia, just one patient from 707 survivors has developed persistent asymptomatic B. pseudomallei carriage. To better understand the mechanisms behind this unique scenario, we performed whole-genome analysis of two strains isolated 139 months apart. During this period, B. pseudomallei underwent several adaptive changes. Of 23 point mutations, 78% were nonsynonymous and 43% were predicted to be deleterious to gene function, demonstrating a strong propensity for positive selection. Notably, a nonsense mutation inactivated the universal stress response sigma factor RpoS, with pleiotropic implications. The genome underwent substantial reduction, with four deletions in chromosome 2 resulting in the loss of 221 genes. The deleted loci included genes involved in secondary metabolism, environmental survival, and pathogenesis. Of 14 indels, 11 occurred in coding regions and 9 resulted in frameshift mutations that dramatically affected predicted gene products. Disproportionately, four indels affected lipopolysaccharide biosynthesis and modification. Finally, we identified a frameshift mutation in both P314 isolates within wcbR, an important component of the capsular polysaccharide I locus, suggesting virulence attenuation early in infection. Our study illustrates a unique clinical case that contrasts a high-consequence infectious agent with a long-term commensal infection and provides further insights into bacterial evolution within the human host. IMPORTANCE Some bacterial pathogens establish long-term infections that are difficult or impossible to eradicate with current treatments. Rapid advances in genome sequencing technologies provide a powerful tool for understanding bacterial persistence within the human host. Burkholderia pseudomallei is considered a highly pathogenic bacterium because infection is commonly fatal. Here, we document within-host evolution of B. pseudomallei in a unique case of human infection with ongoing chronic carriage. Genomic comparison of isolates obtained 139 months (11.5 years) apart showed a strong signal of adaptation within the human host, including inactivation of virulence and immunogenic factors, and deletion of pathways involved in environmental survival. Two global regulatory genes were mutated in the 139-month isolate, indicating extensive regulatory changes favoring bacterial persistence. Our study provides insights into B. pseudomallei pathogenesis and, more broadly, identifies parallel evolutionary mechanisms that underlie chronic persistence of all bacterial pathogens.

Burkholderia ubonensis is an environmental bacterium belonging to the Burkholderia cepacia complex (Bcc), a group of genetically related organisms that are associated with opportunistic but generally nonfatal infections in healthy individuals. In contrast, the near-neighbour species Burkholderia pseudomallei causes melioidosis, a disease that can be fatal in up to 95% of cases if left untreated. B. ubonensis is frequently misidentified as B. pseudomallei from soil samples using selective culturing on Ashdown's medium, reflecting both the shared environmental niche and morphological similarities of these species. Additionally, B. ubonensis shows potential as an important biocontrol agent in B. pseudomallei-endemic regions as certain strains possess antagonistic properties towards B. pseudomallei. Current methods for characterising B. ubonensis are laborious, time-consuming and costly, and as such this bacterium remains poorly studied. The aim of our study was to develop a rapid and inexpensive real-time PCR-based assay specific for B. ubonensis. We demonstrate that a novel B. ubonensis-specific assay, Bu550, accurately differentiates B. ubonensis from B. pseudomallei and other species that grow on selective Ashdown's agar. We anticipate that Bu550 will catalyse research on B. ubonensis by enabling rapid identification of this organism from Ashdown's-positive colonies that are not B. pseudomallei.

Burkholderia pseudomallei is a gram-negative bacterium that causes the serious human disease, melioidosis. There is no vaccine against melioidosis and it can be fatal if not treated with a specific antibiotic regimen, which typically includes the third-generation cephalosporin, ceftazidime (CAZ). There have been several resistance mechanisms described for B. pseudomallei, of which the best described are amino acid changes that alter substrate specificity in the highly conserved class A β-lactamase, PenA. In the current study, we sequenced penA from isolates sequentially derived from two melioidosis patients with wild-type (1.5 µg/mL) and, subsequently, resistant (16 or ≥256 µg/mL) CAZ phenotypes. We identified two single-nucleotide polymorphisms (SNPs) that directly increased CAZ hydrolysis. One SNP caused an amino acid substitution (C69Y) near the active site of PenA, whereas a second novel SNP was found within the penA promoter region. In both instances, the CAZ resistance phenotype corresponded directly with the SNP genotype. Interestingly, these SNPs appeared after infection and under selection from CAZ chemotherapy. Through heterologous cloning and expression, and subsequent allelic exchange in the native bacterium, we confirmed the role of penA in generating both low-level and high-level CAZ resistance in these clinical isolates. Similar to previous studies, the amino acid substitution altered substrate specificity to other β-lactams, suggesting a potential fitness cost associated with this mutation, a finding that could be exploited to improve therapeutic outcomes in patients harboring CAZ resistant B. pseudomallei. Our study is the first to functionally characterize CAZ resistance in clinical isolates of B. pseudomallei and to provide proven and clinically relevant signatures for monitoring the development of antibiotic resistance in this important pathogen.

Single nucleotide polymorphisms (SNPs) are abundant in genomes of all species and biologically informative markers extensively used across broad scientific disciplines. Newly identified SNP markers are publicly available at an ever-increasing rate due to advancements in sequencing technologies. Efficient, cost-effective SNP genotyping methods to screen sample populations are in great demand in well-equipped laboratories, but also in developing world situations. Dual Probe TaqMan assays are robust but can be cost-prohibitive and require specialized equipment. The Mismatch Amplification Mutation Assay, coupled with melt analysis (Melt-MAMA), is flexible, efficient and cost-effective. However, Melt-MAMA traditionally suffers from high rates of assay design failures and knowledge gaps on assay robustness and sensitivity. In this study, we identified strategies that improved the success of Melt-MAMA. We examined the performance of 185 Melt-MAMAs across eight different pathogens using various optimization parameters. We evaluated the effects of genome size and %GC content on assay development. When used collectively, specific strategies markedly improved the rate of successful assays at the first design attempt from ~50% to ~80%. We observed that Melt-MAMA accurately genotypes across a broad DNA range (~100 ng to ~0.1 pg). Genomic size and %GC content influence the rate of successful assay design in an independent manner. Finally, we demonstrated the versatility of these assays by the creation of a duplex Melt-MAMA real-time PCR (two SNPs) and conversion to a size-based genotyping system, which uses agarose gel electrophoresis. Melt-MAMA is comparable to Dual Probe TaqMan assays in terms of design success rate and accuracy. Although sensitivity is less robust than Dual Probe TaqMan assays, Melt-MAMA is superior in terms of cost-effectiveness, speed of development and versatility. We detail the parameters most important for the successful application of Melt-MAMA, which should prove useful to the wider scientific community.

The bacterium Burkholderia pseudomallei causes melioidosis, a rare but serious illness that can be fatal if untreated or misdiagnosed. Species-specific PCR assays provide a technically simple method for differentiating B. pseudomallei from near-neighbor species. However, substantial genetic diversity and high levels of recombination within this species reduce the likelihood that molecular signatures will differentiate all B. pseudomallei from other Burkholderiaceae. Currently available molecular assays for B. pseudomallei detection lack rigorous validation across large in silico datasets and isolate collections to test for specificity, and none have been subjected to stringent quality control criteria (accuracy, precision, selectivity, limit of quantitation (LoQ), limit of detection (LoD), linearity, ruggedness and robustness) to determine their suitability for environmental, clinical or forensic investigations. In this study, we developed two novel B. pseudomallei specific assays, 122018 and 266152, using a dual-probe approach to differentiate B. pseudomallei from B. thailandensis, B. oklahomensis and B. thailandensis-like species; other species failed to amplify. Species specificity was validated across a large DNA panel (>2,300 samples) comprising Burkholderia spp. and non-Burkholderia bacterial and fungal species of clinical and environmental relevance. Comparison of assay specificity to two previously published B. pseudomallei-specific assays, BurkDiff and TTS1, demonstrated comparable performance of all assays, providing between 99.7 and 100% specificity against our isolate panel. Last, we subjected 122018 and 266152 to rigorous quality control analyses, thus providing quantitative limits of assay performance. Using B. pseudomallei as a model, our study provides a framework for comprehensive quantitative validation of molecular assays and provides additional, highly validated B. pseudomallei assays for the scientific research community.