Emerg Infect DisEmerging Infect. DisEIDEmerging Infectious Diseases1080-60401080-6059Centers for Disease Control and Prevention24655932396637013-096010.3201/eid2004.130960DispatchDispatchBurkholderia pseudomallei Type G in Western HemisphereB.pseudomallei Type G in Western HemisphereGeeJay E.AllenderChristopher J.TuanyokApichaiElrodMindy G.HoffmasterAlex R.Centers for Disease Control and Prevention, Atlanta, Georgia, USA (J.E. Gee, M.G. Elrod, A.R. Hoffmaster); Northern Arizona University, Flagstaff, Arizona, USA (C.J. Allender, A. Tuanyok)Address for correspondence: Jay E. Gee, Centers for Disease Control and Prevention, 1600 Clifton Rd NE, Mailstop G34, Atlanta, GA 30333, USA; email: xzg4@cdc.gov42014204682684

Burkholderia pseudomallei isolates from the Western Hemisphere are difficult to differentiate from those from regions in which melioidosis is traditionally endemic. We used internal transcribed spacer typing to determine that B. pseudomallei isolates from the Western Hemisphere are consistently type G. Knowledge of this relationship might be useful for epidemiologic investigations.

Keywords: melioidosisWestern Hemisphereinternal transcribed spacermolecular typingBurkholderiapseudomalleitype Gbacteria

Burkholderia pseudomallei is the causative agent of the disease melioidosis. Melioidosis is considered endemic to Southeast Asia and northern Australia. However, sporadic cases do occur elsewhere in the world, especially in tropical areas (1).

The predominant method of molecular subtyping of B. pseudomallei is multilocus sequence typing (MLST), which is based on a comparison of the alleles of 7 housekeeping genes to generate a sequence type (ST) (2). These data can then be analyzed by using tools such as eBURST, which is used to infer phylogenetic patterns (3).

As of May 30, 2013, a total of 3,028 B. pseudomallei isolates were listed in the MLST database (www.mlst.net); these isolates are predominantly from Southeast Asia (1,036 isolates) and Australia (1,776). Some entries are from other Pacific areas (e.g., New Caledonia [9 isolates] and Hong Kong [39]) and other parts of the world (e.g., Africa [8], Europe [15], and the Western Hemisphere [30]) or of unknown origin (32).

For a population study, Pearson et al. analyzed the STs in the B. pseudomallei MLST database along with other data, such as single-nucleotide polymorphisms from whole-genome sequencing data (4). Their study indicated some population structures associated with geographic origin, in particular, clades associated with isolates from Southeast Asia or northern Australia. The data were used to support a hypothesis that B. pseudomallei originated on the Australian continent, spread to Southeast Asia, and from there spread throughout the world (4).

For an isolate of unknown origin, however, MLST alone may not provide information about geographic origin. For example, isolates from the Western Hemisphere have yielded STs unique to that hemisphere, but when these isolates were analyzed by eBURST or other methods, no distinct clade was found. In the clusters in which these STs appear, STs are intermixed with those from regions in which melioidosis is traditionally endemic, such as Southeast Asia (2,4). Therefore, a method for determining whether an isolate of unknown origin is from the Western Hemisphere is desirable.

Recently, Ligouri et al. developed a typing scheme by measuring the length polymorphisms in the 16S–23S internal transcribed spacer (ITS) of Burkholderia spp. (5). The typing scheme consists of 10 types: A, B, C, D, E, F, G, CE, GE, and GC. Ligouri et al. found that some types are unique to a given Burkholderia species (i.e., A = B. thailandensis, B = B. humptydooensis, D = B. oklahomensis, and F = B. cepacia). They also determined that type C could be found in B. mallei and in B. pseudomallei. The remaining 5 types were exclusive to B. pseudomallei (5). Ligouri et al. determined that types C, E, GE, and CE were the predominant types for isolates from northern Australia and Southeast Asia and that type G was rare in Australia (4 isolates) and Southeast Asia (3). They noted, on the basis of a limited number of strains from these regions, that type G was overrepresented in isolates from other parts of the world: Madagascar (1 isolate), Ecuador (2), Puerto Rico (2), Venezuela (1), and Kenya (1). They hypothesized that a genetic bottleneck occurred during the dispersal of type G to regions outside of Southeast Asia and Australia (5).

ITS typing of B. pseudomallei might be a powerful tool for linking cases of melioidosis to regions outside of those in which melioidosis in highly endemic, such as Southeast Asia and northern Australia. To further investigate this trend, we assessed the ITS types of B. pseudomallei from the Western Hemisphere.

The Study

All tested Western Hemisphere isolates from our collection were ITS type G (Table). In addition to those isolates tested, we performed in silico analysis of whole-genome sequencing data from other B. pseudomallei isolates in our collection or publicly available data with origins in the Western Hemisphere and found them to also be ITS type G (Table). As expected, eBURST analysis of the STs from the type G strains that originated in the Western Hemisphere did not yield a discrete clade. These STs are interspersed with other STs that are predominantly from Southeast Asia (Figure). Our results support the findings of Ligouri et al. (5). The predominance of type G in isolates in our panel of isolates from the Western Hemisphere is consistent with a hypothesis that these isolates are derived from a bottleneck that occurred during dispersal to the rest of the world from Southeast Asia.

Analysis results for <italic>Burkholderia pseudomallei</italic> isolates from Western Hemisphere*
Identification no.ITS typeYLF geneMLSTSource
Swiss2010G+92Switzerland ex Martinique 2010
PR1998G†+†92Puerto Rico 1998 (GenBank accession no. pending)
4900CFG†+†92Brazil: cystic fibrosis patient 2007 (GenBank accession no. ARZE00000000)
CA2009G+95USA (California) ex Mexico 2009
FL2009G+297USA (Florida) ex Puerto Rico 2009
FL2012G+297USA (Florida) ex Trinidad 2012
PR2012G+297Puerto Rico 2012
MX2013GND297Mexico 2013
PB08298010G†+426USA (Arizona) locally acquired 2008 (GenBank accession no. ARZO00000000)
CA2012dG+436USA (California) ex Guatemala 2012
2002734728G+518USA (California) Iguana 2007
PB 1007001G+518USA (Arizona) ex Costa Rica 2010
CA2013aGND518USA (California) Iguana 2013
NY2010G+698USA (New York) ex Aruba 2010
724644G+698USA (Massachusetts) ex Aruba 2010
BCC215G†+†ND‡Brazil (Ceara) 2003 (GenBank accession no. ABBR00000000)

*ITS, internal transcribed spacer; YLF, Yersinia-like fimbrial; MLST, multilocus sequence type; ex, diagnosis made in first location but infection acquired in second location; ND, not determined.
†Determined by analysis of whole-genome sequence.
‡Allele for gltB not available.

Diagram of eBURST (3) analysis of multilocus sequence types (STs) of Burkholderia pseudomallei isolates; default settings were used. Dots represent sequence types. Blue dot represents the calculated primary founder, yellow dots represent calculated subgroup founders, and black dots represent the remaining STs. The sizes of dots are proportional to the number of isolates in the database representing a given ST. A single line between sequence types indicates that they are single-locus variants of each other. Red arrows indicate positions of STs associated with isolates from the Western Hemisphere and that are type G according to internal transcribed spacer typing.

Other supporting evidence is provided by testing for the Yersinia-like fimbrial (YLF) gene by YLF PCR. The presence of the YLF gene is associated with B. pseudomallei isolates from Southeast Asia (6). For those strains tested, all were positive for YLF or the YLF gene was present according to in silico analysis of whole-genome sequencing data when available (Table).

Conclusions

In melioidosis cases for which the origin is unclear, such as for patients with unknown exposure histories or histories of travel to multiple regions to which melioidosis is endemic, ITS typing along with other molecular epidemiologic tools might be useful for assessing the B. pseudomallei origins. More insight into the relationship of Western Hemisphere isolates to isolates from regions to which melioidosis is highly endemic might come from whole-genome sequencing. As more isolates are analyzed, such studies might enable higher confidence in the geographic origin of isolates.

Suggested citation for this article: Gee JE, Allender CJ, Tuanyok A, Elrod MG, Hoffmaster AR. Burkholderia pseudomallei type G in Western Hemisphere. Emerg Infect Dis [Internet]. 2014 Apr [date cited]. http://dx.doi.org/10.3201/eid2004.130960

A.T. is supported by the United States Department of Homeland Security grant no. HSHQDC-10-C-00135. This study used the Multi Locus Sequence Typing website (www.mlst.net) at Imperial College London, developed by David Aanensen and funded by the Wellcome Trust.

Dr Gee is a research biologist in the Bacterial Special Pathogens Branch, Division of High-Consequence Pathogens and Pathology, National Center for Emerging and Zoonotic Infectious Diseases, Centers for Disease Control and Prevention. His research interests focus on the molecular epidemiology of B. pseudomallei and B. mallei and on other emerging bacterial pathogens.

ReferencesWiersinga WJ, Currie BJ, Peacock SJ. Melioidosis. N Engl J Med. 2012;367:103544 . 10.1056/NEJMra120469922970946Godoy D, Randle G, Simpson AJ, Aanensen DM, Pitt TL, Kinoshita R, Multilocus sequence typing and evolutionary relationships among the causative agents of melioidosis and glanders, Burkholderia pseudomallei and Burkholderia mallei. J Clin Microbiol. 2003;41:206879 10.1128/JCM.41.5.2068-2079.200312734250Feil EJ, Li BC, Aanensen DM, Hanage WP, Spratt BG. eBURST: inferring patterns of evolutionary descent among clusters of related bacterial genotypes from multilocus sequence typing data. J Bacteriol. 2004;186:151830 10.1128/JB.186.5.1518-1530.200414973027Pearson T, Giffard P, Beckstrom-Sternberg S, Auerbach R, Hornstra H, Tuanyok A, Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer. BMC Biol. 2009;7:78 . 10.1186/1741-7007-7-7819922616Liguori AP, Warrington SD, Ginther JL, Pearson T, Bowers J, Glass MB, Diversity of 16S–23S rDNA internal transcribed spacer (ITS) reveals phylogenetic relationships in Burkholderia pseudomallei and its near-neighbors. PLoS ONE. 2011;6:e29323 10.1371/journal.pone.002932322195045Tuanyok A, Auerbach RK, Brettin TS, Bruce DC, Munk AC, Detter JC, A horizontal gene transfer event defines two distinct groups within Burkholderia pseudomallei that have dissimilar geographic distributions. J Bacteriol. 2007;189:90449 10.1128/JB.01264-0717933898