PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of microgenLink to Publisher's site
 
Microb Genom. 2017 November; 3(11): e000139.
Published online 2017 November 14. doi:  10.1099/mgen.0.000139
PMCID: PMC5729916

Suspected cases of intracontinental Burkholderia pseudomallei sequence type homoplasy resolved using whole-genome sequencing

Abstract

Burkholderia pseudomallei is a Gram-negative environmental bacterium that causes melioidosis, a disease of high mortality in humans and animals. Multilocus sequence typing (MLST) is a popular and portable genotyping method that has been used extensively to characterise the genetic diversity of B. pseudomallei populations. MLST has been central to our understanding of the underlying phylogeographical signal present in the B. pseudomallei genome, revealing distinct populations on both the intra- and the inter-continental level. However, due to its high recombination rate, it is possible for B. pseudomallei isolates to share the same multilocus sequence type (ST) despite being genetically and geographically distinct, with two cases of ‘ST homoplasy’ recently reported between Cambodian and Australian B. pseudomallei isolates. This phenomenon can dramatically confound conclusions about melioidosis transmission patterns and source attribution, a critical issue for bacteria such as B. pseudomallei that are of concern due to their potential for use as bioweapons. In this study, we used whole-genome sequencing to identify the first reported instances of intracontinental ST homoplasy, which involved ST-722 and ST-804 B. pseudomallei isolates separated by large geographical distances. In contrast, a third suspected homoplasy case was shown to be a true long-range (460 km) dispersal event between a remote Australian island and the Australian mainland. Our results show that, whilst a highly useful and portable method, MLST can occasionally lead to erroneous conclusions about isolate origin and disease attribution. In cases where a shared ST is identified between geographically distant locales, whole-genome sequencing should be used to resolve strain origin.

Keywords: genomics, melioidosis, phylogenetics, homoplasy, source tracing

Abbreviations

BTFC, Burkholderia thailandensis-like flagella and chemotaxis; DPMS, Darwin Prospective Melioidosis Study; ML, maximum-likelihood; MLST, multilocus sequence typing; MP, maximum-parsimony; SNP, single-nucleotide polymorphism; ST, sequence type; WGS, whole-genome sequencing; YLF, Yersinia-like fimbriae.

Data Summary

1. Whole-genome sequencing data have been deposited in the NCBI Sequence Read Archive (SRA) under BioProject accession number PRJNA393909 (url – https://www.ncbi.nlm.nih.gov/bioproject/PRJNA393909).

2. The SRA accession number for the raw sequence data of MSHR0052 is SRR5818275.

3. The SRA accession number for the raw sequence data of MSHR0116 is SRR5818274.

4. The SRA accession number for the raw sequence data of MSHR3528 is SRR5818273.

5. The SRA accession number for the raw sequence data of MSHR9076 is SRR5818272.

6. The SRA accession number for the raw sequence data of MSHR4608 is SRR5818271.

Impact Statement

Melioidosis is a high-mortality tropical disease that is acquired following contact with soil or surface waters containing Burkholderia pseudomallei. Characterising the spread and genetic diversity of B. pseudomallei is critical for understanding melioidosis epidemiology. Multilocus sequence typing (MLST) is a popular genetic tool for fingerprinting bacterial populations. However, the relatively low resolution of this method means that isolates can occasionally be assigned the same sequence type (ST) despite being otherwise unrelated, a phenomenon known as ST homoplasy. Using whole-genome sequencing (WGS), two cases of ST homoplasy occurring between Asian and Australian B. pseudomallei isolates have been previously identified. Here, we used WGS to uncover two intracontinental cases of ST homoplasy in isolates from northern Australia, with a third suspected case shown to be a true long-range (460 km) dispersal event. This latter finding points to the relatively recent transmission of a B. pseudomallei clone into a new region.

Introduction

Burkholderia pseudomallei is a Gram-negative soil-dwelling bacterium that causes melioidosis, a potentially fatal tropical infectious disease that is hyper-endemic in northern Australia and parts of Southeast Asia [1]. B. pseudomallei is an opportunistic organism that infects both humans and animals, although person-to-person and zoonotic transmissions are exceedingly rare [2]. B. pseudomallei infections are almost exclusively acquired via contact with a contaminated environmental source [3]. The potential use of B. pseudomallei as a biological weapon [4], combined with limited treatment options, difficulty in diagnosis and current lack of a vaccine [1] make this organism of substantial importance to public health in tropical regions worldwide.

Multilocus sequence typing (MLST) [5] is a widely adopted genotyping method for characterising the population diversity of several pathogenic bacterial species, including B. pseudomallei. The B. pseudomallei MLST scheme [6] uses seven housekeeping loci (ace, gltB, gmhD, lepA, lipA, narK and ndh) encoded solely on chromosome 1, and yields sequence data of a combined total length of 3401 bp. In B. pseudomallei, MLST has been used extensively to identify strain relatedness [6, 7], for point-source attribution [8, 9] and to document the geographical distribution of isolates, particularly those from the ‘Top End’ region of the Northern Territory, Australia [10, 11]. A major finding to arise from this body of work is that, despite severe weather events (e.g. cyclones, floods), animal migration and anthropogenic factors, B. pseudomallei dispersal in the environment is surprisingly limited [10]. Indeed, identical sequence types (STs) are typically found across a maximum linear distance of only 45 km [11].

High rates of genetic recombination [12], coupled with genetic drift over large timescales, can lead to occasional instances where B. pseudomallei strains have converged on the same ST by chance rather than by sharing a recent common ancestor. Such cases are termed ST homoplasy [13]. To date, there have been two documented intercontinental cases of ST homoplasy in B. pseudomallei, which occurred in isolates retrieved from both Cambodia and Australia [13].

Unlike MLST, whole-genome sequencing (WGS) provides high resolution of closely related isolates and can more accurately characterise the origin of an isolate. Phylogenetic reconstruction of B. pseudomallei populations using WGS data is a robust way to identify the geographical origin of strains on a continental level [12–16]. In the present study we use WGS to uncover, for the first time, two intracontinental cases of ST homoplasy in B. pseudomallei, with a third suspected case instead being identified as a true long-distance ST dispersal.

Methods

The ongoing Darwin Prospective Melioidosis Study (DPMS) [17] has documented and collected B. pseudomallei isolates from >1000 human melioidosis cases occurring in the ‘Top End’ region of the Northern Territory since October 1989. This collection also includes isolates from animal infections and from environmental sampling efforts undertaken across the Northern Territory since 1991. Most DPMS isolates have previously been genotyped using a combination of conventional MLST [6] or MLST derived from WGS [18]. GPS data were collected at the time of sampling for the environmental isolates, whereas the provenance of clinical isolates was based on patient travel history taken at the time of diagnosis. Using this dataset, isolates belonging to the same ST that were separated by a linear distance of >100 km were investigated as potential cases of intracontinental ST homoplasy.

Comparative analysis of 145 publicly available B. pseudomallei genomes that represent a global snapshot of isolates [15], plus five new genomes sequenced as part of this study, was carried out to determine the genetic relatedness of the suspected ST homoplasy cases. ST-149, ST-722 or ST-804 isolates (MSHR0052, MSHR0116, MSHR3528, MSHR9076 and MSHR4608) that lacked existing WGS data were sequenced using the Illumina HiSeq2500 platform (Australian Genome Research Facility). Publicly available data belonging to previously sequenced ST-149 isolates (MSHR0503, MSHR0356 and MSHR4300) were retrieved from the NCBI SRA database. Orthologous, biallelic single-nucleotide polymorphisms (SNPs) were identified from the Illumina data using the default settings of SPANDx v3.2 [19]. The closed Australian B. pseudomallei genome MSHR1153 [20] was used as the reference for read mapping for the global phylogeny, whereas the environmental ST-149 isolate MSHR4300 [20] was chosen as the reference for examining ST-149 diversity as this genome is a high-quality assembly, it has the fewest number of contigs (n=3) and it has excellent provenance due to its environmental origin. A draft assembly of an additional ST-149 isolate, MSHR0503, was assembled using the Microbial Genome Assembler Pipeline (MGAP) v1.0 [21]. Both maximum-parsimony (MP) and maximum-likelihood (ML) phylogenetic trees were reconstructed based on 207 694 SNPs identified among the 150 genomes using PAUP* 4.0a151 [22] and RAxML v8.0.X (GTRCAT model) [23]. Recombinogenic SNPs were identified using Gubbins v2.2.0 (default parameters) [24]. In silico MLST was performed using the BIGSdb tool [25], which is available on the B. pseudomallei MLST website (http://pubmlst.org/bpseudomallei/).

Results

We identified eight DPMS isolates belonging to three STs (ST-149, ST-722 and ST-804) as possible cases of intracontinental homoplasy, with isolates of these STs being separated by 460, 300 and 1000 km, respectively (Table 1). Based on a comparative analysis of 150 global genomes (Figs 1 and S1, available in the online Supplementary Material), geographically distant isolates belonging to ST-722 (MSHR0052 and MSHR9076) or ST-804 (MSHR3528 and MSHR4608) did not group closely together on the whole-genome level, differing by 21 211 and 20 567 SNPs, respectively. In contrast, the ST-149 environmental isolates MSHR0503 and MSHR4300 were separated by only 404 SNPs (Fig. 2a), which is comparable to previously observed intra-ST diversity [14]. Two clinical ST-149 isolates, MSHR0356 and MSHR0116, both of which were retrieved from melioidosis patients residing in a remote Australian island community where environmental isolate MSHR0503 was sampled [26], also grouped with MSHR4300, differing from this strain by 408 SNPs each. The three ST-149 isolates from the remote island community (MSHR0503, MSHR0356, MSHR0116) differed by a total of just 19 SNPs, suggesting very recent shared ancestry. When MSHR4300 was excluded and a draft assembly of MSHR0503 was used as the reference for read mapping, 18 SNPs were identified among these three strains, and only three SNPs separated the clinical isolates (MSHR0116 and MSHR0356) (Fig. S2).

Table 1.
Burkholderia pseudomallei isolates used in this study
Fig. 1.
MP phylogeny of 150 Burkholderia pseudomallei genomes constructed using 207 694 orthologous, biallelic SNPs. MSHR1153 was used as a reference, and the tree was rooted with MSHR0668. Suspected cases of ST homoplasy are coloured as follows: red=ST-149 ...
Fig. 2.
MP phylogenetic analysis of Australian ST-149 isolates, reconstructed using SNPs identified from comparative genomic analysis (reference genome: MSHR4300). (a) Based on 416 SNPs identified among the four ST-149 isolates, 404 SNPs separate the two environmental ...

Given that B. pseudomallei is highly recombinogenic, SNP differences separating the ST-149 strains could have arisen due to one or more recombination events, which can confound molecular clock analyses. We therefore assessed the contribution of recombination in the ST-149 isolates. In total, 416 biallelic SNPs were found among the four ST-149 isolates (Fig. 2a). Removal of the 241 recombinogenic SNPs did not affect the topology of the MP or ML phylogenies (Figs 2b and S3). Using an average substitution rate of 5.9×10−7 substitutions per site per year (95 % CI=2.0×10−7 to 9.8×10−7) [16, 27–29] and following the removal of recombinogenic SNPs, divergence of the two geographically disparate ST-149 isolates was estimated to have occurred in 1984 (95 % CI=1947 to 1992).

Finally, we investigated three variable genetic markers with known geographical [30] or clinical [31] associations in B. pseudomallei: the mutually exclusive bimABp and bimABm loci [32], the variably present virulence gene fhaB3 [33], and the mutually exclusive Burkholderia thailandensis-like flagella and chemotaxis (BTFC) cluster and Yersinia-like fimbriae (YLF) cluster loci [30]. These genetic markers were developed in the pre-genomics era to provide a rapid determination of strain relatedness, virulence potential and origin [30, 31]. All isolates were positive for the virulence genes fhaB3 [33] and bimABp [32], except for MSHR9076 (ST-722), which harboured the B. mallei-like neurotropic variant bimABm [34] instead of bimABp. In addition, MSHR0052 (ST-722) was YLF-positive, whereas all other isolates possessed the BTFC cluster, including MSHR9076 (Table 1). These three markers were unable to differentiate the ST-804 strains; however, this result was not surprising as bimABp, fhaB3 and BTFC comprise the most common genotype in Australian isolates [30, 31]. The three variable genetic markers were also identical across the four ST-149 isolates, consistent with their close phylogenomic relatedness.

Discussion

MLST remains a well-adopted genotyping method for bacterial epidemiological studies due to its robustness, portability and data accessibility [5]. Although MLST can provide a snapshot of the evolutionary history and relatedness of an isolate, there are certain instances where MLST does not perform as intended and may in fact provide misleading results, particularly for highly recombinogenic species such as the melioidosis bacterium, B. pseudomallei [35]. WGS has previously been used to resolve two such cases of ST homoplasy that occurred between Australian and Cambodian isolates of B. pseudomallei [13]. Here, we performed WGS analysis on eight Australian B. pseudomallei isolates from three STs (ST-149, ST-722 and ST-804) suspected of also being homoplastic due to shared STs despite wide geographical separation.

Phylogenomic analysis confirmed that the ST-722 and ST-804 isolates represent two new ST homoplasy cases, with isolates of each ST separated by >20 000 SNPs, although all still resided within the Australian clade (Figs 1 and S1). This genetic distance is characteristic of Australian isolates belonging to different STs. In addition, the two ST-722 isolates differed at the YLF/BTFC cluster and bimABp/Bm loci, providing further evidence of their dissimilarity. Our results add to the two previously documented intercontinental ST homoplasy cases reported by De Smet et al. [13]. Although ST homoplasy cases are rare, they represent important instances of genotyping failure that can significantly confound analyses that rely on MLST data, particularly where accurate source attribution is paramount, such as in forensic or outbreak investigations. It is therefore crucial that high-resolution techniques such as WGS are used to confirm unusual or unexpected genotyping results, especially in highly recombinogenic pathogens such as B. pseudomallei.

Unlike ST-722 and ST-804, ST-149 isolates showed evidence of relatively recent diversification. Despite being retrieved from soil samples separated by 460 km, the environmental isolates differed by just 404 SNPs (Fig. 2a), or 165 non-recombinogenic SNPs (Fig. 2b). Using molecular clocks to calculate divergence has gained popularity in recent years. However, such estimates in B. pseudomallei have proven to be non-trivial due to uncertain generation times and difficulties in estimating the contribution of recombination to genetic diversity. Nevertheless, our results estimate that divergence of ST-149 occurred in 1984 (95 % CI=1947 to 1992). Based on this estimate, we hypothesise that an anthropogenic, zoogenic, ecological or severe weather event drove the long-range dispersal of ST-149 between these two locations sometime between 18 and 63 years before MSHR4300 was isolated in 2010.

The ST-149 case is not the first instance of long-range B. pseudomallei dispersal. A possible intracontinental long-range dispersal event is thought to have occurred in Western Australia across a linear distance of 500 km [36], although WGS would be needed to rule out ST homoplasy in this case. The largest known long-range dispersal event has recently been documented between Asia and Australia, with the unprecedented transmission of an Asian B. pseudomallei clone, ST-562, into Darwin, the capital city of the Northern Territory [14]. Although the precise origin and transmission of this clone into Australia has not yet been determined, its restricted geographical locale and very low levels of genomic diversity point to a very recent introduction event [14, 36]. Our current study adds to the growing body of work supporting the potential long-range transmission of B. pseudomallei across large geographical distances. The precise factors driving these dispersal events requires ongoing monitoring and investigation.

Conclusion

We report, for the first time, two instances of intracontinental ST homoplasy in B. pseudomallei. These examples illustrate the occasional limits of MLST for accurate source attribution of highly recombinogenic species. We also report an unusual instance of long-range B. pseudomallei dispersal, which may have been aided by recent anthropogenic, ecological or weather-related factors.

Data bibliography

  1. MSHR0503
  2. Sahl JW, Vazquez AJ, Hall CM et al. BioProject PRJNA285704, Sequence Read Archive Number SRR2975627 (2016) https://www.ncbi.nlm.nih.gov/sra/SRR2975627
  3. MSHR0356
  4. Nandi T, Holden MT, Didelot X, et al. BioProject PRJEB2866 European Nucleotide Archive ERR311049 (2015), https://www.ncbi.nlm.nih.gov/sra/ERX284330
  5. MSHR4300
  6. Johnson SL, Baker AL, Chain PS, et al. NCBI BioProject PRJNA230723 (2015) https://www.ncbi.nlm.nih.gov/Traces/wgs/?val=JPQI01
  7. MSHR1153
  8. Bunnell A, Chain PS, Chertkov O et al. NCBI BioProject PRJNA231860 (2015) https://www.ncbi.nlm.nih.gov/bioproject/231860

Funding information

This study was funded by the National Health and Medical Research Council (NHMRC) through Project Grants 1046812, 1098337 and 1131932 (the HOT NORTH initiative). AA is supported by a Research Training Program Scholarship from the Australian Government and an NHMRC Centres for Research Excellence top-up scholarship (1078557). DSS is supported by an Advance Queensland Fellowship (AQRF13016-17RD2). TMH is supported by an NHMRC Centres for Research Excellence post-doctoral fellowship (1078557). EPP is supported by a University of the Sunshine Coast Fellowship.

Conflicts of interest

The authors declare that there are no conflicts of interest.

Ethical statement

Approval for this study was obtained as described elsewhere [17].

Supplementary Data

506Supplementary File 1:

References

1. Wiersinga WJ, Currie BJ, Peacock SJ. Melioidosis. N Engl J Med. 2012;367:1035–1044. doi: 10.1056/NEJMra1204699. [PubMed] [Cross Ref]
2. Limmathurotsakul D, Peacock SJ. Melioidosis: a clinical overview. Br Med Bull. 2011;99:125–139. doi: 10.1093/bmb/ldr007. [PubMed] [Cross Ref]
3. Barnes JL, Ketheesan N. Route of infection in melioidosis. Emerg Infect Dis. 2005;11:638–639. doi: 10.3201/eid1104.041051. [PMC free article] [PubMed] [Cross Ref]
4. Gilad J, Harary I, Dushnitsky T, Schwartz D, Amsalem Y. Burkholderia mallei and Burkholderia pseudomallei as bioterrorism agents: national aspects of emergency preparedness. Isr Med Assoc J. 2007;9:499–503. [PubMed]
5. Maiden MC, Bygraves JA, Feil E, Morelli G, Russell JE, et al. Multilocus sequence typing: a portable approach to the identification of clones within populations of pathogenic microorganisms. Proc Natl Acad Sci USA. 1998;95:3140–3145. doi: 10.1073/pnas.95.6.3140. [PubMed] [Cross Ref]
6. Godoy D, Randle G, Simpson AJ, Aanensen DM, Pitt TL, et al. Multilocus sequence typing and evolutionary relationships among the causative agents of melioidosis and glanders, Burkholderia pseudomallei and Burkholderia mallei. J Clin Microbiol. 2003;41:2068–2079. doi: 10.1128/JCM.41.5.2068-2079.2003. [PMC free article] [PubMed] [Cross Ref]
7. McCombie RL, Finkelstein RA, Woods DE. Multilocus sequence typing of historical Burkholderia pseudomallei isolates collected in Southeast Asia from 1964 to 1967 provides insight into the epidemiology of melioidosis. J Clin Microbiol. 2006;44:2951–2962. doi: 10.1128/JCM.00725-06. [PMC free article] [PubMed] [Cross Ref]
8. Mayo M, Kaesti M, Harrington G, Cheng AC, Ward L, et al. Burkholderia pseudomallei in unchlorinated domestic bore water, Tropical Northern Australia. Emerg Infect Dis. 2011;17:1283–1285. doi: 10.3201/eid1707.100614. [PMC free article] [PubMed] [Cross Ref]
9. McRobb E, Sarovich DS, Price EP, Kaestli M, Mayo M, et al. Tracing melioidosis back to the source: using whole-genome sequencing to investigate an outbreak originating from a contaminated domestic water supply. J Clin Microbiol. 2015;53:1144–1148. doi: 10.1128/JCM.03453-14. [PMC free article] [PubMed] [Cross Ref]
10. McRobb E, Kaestli M, Price EP, Sarovich DS, Mayo M, et al. Distribution of Burkholderia pseudomallei in northern Australia, a land of diversity. Appl Environ Microbiol. 2014;80:3463–3468. doi: 10.1128/AEM.00128-14. [PMC free article] [PubMed] [Cross Ref]
11. Chapple SN, Price EP, Sarovich DS, McRobb E, Mayo M, et al. Burkholderia pseudomallei genotype distribution in the Northern Territory, Australia. Am J Trop Med Hyg. 2016;94:68–72. doi: 10.4269/ajtmh.15-0627. [PMC free article] [PubMed] [Cross Ref]
12. Pearson T, Giffard P, Beckstrom-Sternberg S, Auerbach R, Hornstra H, et al. Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer. BMC Biol. 2009;7:78. doi: 10.1186/1741-7007-7-78. [PMC free article] [PubMed] [Cross Ref]
13. De Smet B, Sarovich DS, Price EP, Mayo M, Theobald V, et al. Whole-genome sequencing confirms that Burkholderia pseudomallei multilocus sequence types common to both Cambodia and Australia are due to homoplasy. J Clin Microbiol. 2015;53:323–326. doi: 10.1128/JCM.02574-14. [PMC free article] [PubMed] [Cross Ref]
14. Price EP, Sarovich DS, Smith EJ, Machunter B, Harrington G, et al. Unprecedented melioidosis cases in Northern Australia caused by an Asian Burkholderia pseudomallei strain identified by using large-scale comparative genomics. Appl Environ Microbiol. 2016;82:954–963. doi: 10.1128/AEM.03013-15. [PMC free article] [PubMed] [Cross Ref]
15. Sarovich DS, Garin B, de Smet B, Kaestli M, Mayo M, et al. Phylogenomic analysis reveals an Asian origin for African Burkholderia pseudomallei and further supports melioidosis endemicity in Africa. mSphere. 2016;1:e00089-15 doi: 10.1128/mSphere.00089-15. [PMC free article] [PubMed] [Cross Ref]
16. Chewapreecha C, Holden MT, Vehkala M, Välimäki N, Yang Z, et al. Global and regional dissemination and evolution of Burkholderia pseudomallei. Nat Microbiol. 2017;2:16263. doi: 10.1038/nmicrobiol.2016.263. [PMC free article] [PubMed] [Cross Ref]
17. Currie BJ, Ward L, Cheng AC. The epidemiology and clinical spectrum of melioidosis: 540 cases from the 20 year Darwin Prospective Study. PLoS Negl Trop Dis. 2010;4:e900 doi: 10.1371/journal.pntd.0000900. [PMC free article] [PubMed] [Cross Ref]
18. Larsen MV, Cosentino S, Rasmussen S, Friis C, Hasman H, et al. Multilocus sequence typing of total-genome-sequenced bacteria. J Clin Microbiol. 2012;50:1355–1361. doi: 10.1128/JCM.06094-11. [PMC free article] [PubMed] [Cross Ref]
19. Sarovich DS, Price EP. SPANDx: a genomics pipeline for comparative analysis of large haploid whole genome re-sequencing datasets. BMC Res Notes. 2014;7:618. doi: 10.1186/1756-0500-7-618. [PMC free article] [PubMed] [Cross Ref]
20. Johnson SL, Baker AL, Chain PS, Currie BJ, Daligault HE, et al. Whole-genome sequences of 80 environmental and clinical isolates of Burkholderia pseudomallei. Genome Announc. 2015;3:e01282-14 doi: 10.1128/genomeA.01282-14. [PMC free article] [PubMed] [Cross Ref]
21. Sarovich D. Microbial Genome Assembler Pipeline (MGAP): Initial Release. 2017. https://github.com/dsarov/MGAP---Microbial-Genome-Assembler-Pipeline
22. Swofford DL. PAUP*: Phylogenetic Analysis Using Parsimony (*and other methods) 4.0.b5. Sunderland, MA: Sinauer Associates; 2001.
23. Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30:1312–1313. doi: 10.1093/bioinformatics/btu033. [PMC free article] [PubMed] [Cross Ref]
24. Croucher NJ, Page AJ, Connor TR, Delaney AJ, Keane JA, et al. Rapid phylogenetic analysis of large samples of recombinant bacterial whole genome sequences using Gubbins. Nucleic Acids Res. 2015;43:e15 doi: 10.1093/nar/gku1196. [PMC free article] [PubMed] [Cross Ref]
25. Jolley KA, Maiden MC. BIGSdb: scalable analysis of bacterial genome variation at the population level. BMC Bioinformatics. 2010;11:595. doi: 10.1186/1471-2105-11-595. [PMC free article] [PubMed] [Cross Ref]
26. Sarovich DS, Chapple SNJ, Price EP, Mayo M, Holden MTG, et al. Whole-genome sequencing to investigate a non-clonal melioidosis cluster on a remote Australian island. Microb Genom. 2017;3:e0.000117 doi: 10.1099/mgen.0.000117. [PMC free article] [PubMed] [Cross Ref]
27. Viberg LT, Sarovich DS, Kidd TJ, Geake JB, Bell SC, et al. Within-host evolution of Burkholderia pseudomallei during chronic infection of seven Australasian cystic fibrosis patients. MBio. 2017;8:e00356-17 doi: 10.1128/mBio.00356-17. [PMC free article] [PubMed] [Cross Ref]
28. Lieberman TD, Michel JB, Aingaran M, Potter-Bynoe G, Roux D, et al. Parallel bacterial evolution within multiple patients identifies candidate pathogenicity genes. Nat Genet. 2011;43:1275–1280. doi: 10.1038/ng.997. [PMC free article] [PubMed] [Cross Ref]
29. Silva IN, Santos PM, Santos MR, Zlosnik JE, Speert DP, et al. Long-term evolution of Burkholderia multivorans during a chronic cystic fibrosis infection reveals shifting forces of selection. mSystems. 2016;1:e00029-16 doi: 10.1128/mSystems.00029-16. [PMC free article] [PubMed] [Cross Ref]
30. Tuanyok A, Auerbach RK, Brettin TS, Bruce DC, Munk AC, et al. A horizontal gene transfer event defines two distinct groups within Burkholderia pseudomallei that have dissimilar geographic distributions. J Bacteriol. 2007;189:9044–9049. doi: 10.1128/JB.01264-07. [PMC free article] [PubMed] [Cross Ref]
31. Sarovich DS, Price EP, Webb JR, Ward LM, Voutsinos MY, et al. Variable virulence factors in Burkholderia pseudomallei (melioidosis) associated with human disease. PLoS One. 2014;9:e91682 doi: 10.1371/journal.pone.0091682. [PMC free article] [PubMed] [Cross Ref]
32. Stevens MP, Stevens JM, Jeng RL, Taylor LA, Wood MW, et al. Identification of a bacterial factor required for actin-based motility of Burkholderia pseudomallei. Mol Microbiol. 2005;56:40–53. doi: 10.1111/j.1365-2958.2004.04528.x. [PubMed] [Cross Ref]
33. Sim SH, Yu Y, Lin CH, Karuturi RK, Wuthiekanun V, et al. The core and accessory genomes of Burkholderia pseudomallei: implications for human melioidosis. PLoS Pathog. 2008;4:e1000178 doi: 10.1371/journal.ppat.1000178. [PMC free article] [PubMed] [Cross Ref]
34. Morris JL, Fane A, Sarovich DS, Price EP, Rush CM, et al. Increased neurotropic threat from Burkholderia pseudomallei strains with a B. mallei-like variation in the bimA motility gene, Australia. Emerg Infect Dis. 2017;23:eid2305.151417 doi: 10.3201/eid2305.151417. [PMC free article] [PubMed] [Cross Ref]
35. Turner KM, Hanage WP, Fraser C, Connor TR, Spratt BG. Assessing the reliability of eBURST using simulated populations with known ancestry. BMC Microbiol. 2007;7:30. doi: 10.1186/1471-2180-7-30. [PMC free article] [PubMed] [Cross Ref]
36. Inglis TJ, Levy A, Merritt AJ, Hodge M, McDonald R, et al. Melioidosis risk in a tropical industrial environment. Am J Trop Med Hyg. 2009;80:78–84. [PubMed]

Articles from Microbial Genomics are provided here courtesy of Microbiology Society