Length polymorphisms within the 16S-23S ribosomal DNA internal transcribed spacer (ITS) have been described as stable genetic markers for studying bacterial phylogenetics. In this study, we used these genetic markers to investigate phylogenetic relationships in Burkholderia pseudomallei and its near-relative species. B. pseudomallei is known as one of the most genetically recombined bacterial species. In silico analysis of multiple B. pseudomallei genomes revealed approximately four homologous rRNA operons and ITS length polymorphisms therein. We characterized ITS distribution using PCR and analyzed via a high-throughput capillary electrophoresis in 1,191 B. pseudomallei strains. Three major ITS types were identified, two of which were commonly found in most B. pseudomallei strains from the endemic areas, whereas the third one was significantly correlated with worldwide sporadic strains. Interestingly, mixtures of the two common ITS types were observed within the same strains, and at a greater incidence in Thailand than Australia suggesting that genetic recombination causes the ITS variation within species, with greater recombination frequency in Thailand. In addition, the B. mallei ITS type was common to B. pseudomallei, providing further support that B. mallei is a clone of B. pseudomallei. Other B. pseudomallei near-neighbors possessed unique and monomorphic ITS types. Our data shed light on evolutionary patterns of B. pseudomallei and its near relative species.
Citation: Liguori AP, Warrington SD, Ginther JL, Pearson T, Bowers J, Glass MB, et al. (2011) Diversity of 16S-23S rDNA Internal Transcribed Spacer (ITS) Reveals Phylogenetic Relationships in Burkholderia pseudomallei and Its Near-Neighbors. PLoS ONE 6(12): e29323. doi:10.1371/journal.pone.0029323
Editor: Paul J. Planet, Columbia University, United States of America
Received: July 4, 2011; Accepted: November 24, 2011; Published: December 14, 2011
Copyright: © 2011 Liguori et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
Funding: This work was supported by National Institutes of Health (NIH) - National Institute of Allergy and Infectious Diseases (NIAID) grant nos. U01 AI075568 (PI: Keim) and U54 AI065359 (PI: Tuanyok), and US Department of Homeland Security contract no. HSHQDC-10-C-00149 (PI: Tuanyok). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Competing interests: The authors have declared that no competing interests exist.
Utilizing length polymorphisms within the 16S-23S ribosomal DNA internal transcribed spacer (ITS) as genetic markers, aka PCR ribotyping, was developed by Kostman and colleagues . This technique was first used as an alternative method to traditional 16S ribosomal DNA (rDNA) ribotyping for detecting epidemiological differences between Burkholderia cepacia strains . The technique involves PCR of the 16S-23S internal transcribed spacer (ITS), aka the 16S-23S intergenic spacer region, and assigns phylogenetic groupings based on length polymorphisms, sequence differences or restriction patterns in the ITS , , . Detailed sequence analyses of the ITS in many bacteria reveal its functional role in ribosomal RNA (rRNA) maturation and transcriptional anti-termination , , . In many bacterial species the ITS contains coding sequences for tRNA genes. Consequently, the region is under selective pressure and is fairly conserved, even across bacterial species. Multiple rRNA operons (rrn) have been confirmed within single bacterium, sometimes as many as 15, which often display intragenomic heterogeneity in ITS type , , . Copy number of the rrn operon has been positively associated with bacterial resource utilization and growth rate . ITS variation has been attributed to intragenic evolution as well as the insertion or deletion of sequence blocks; because of this, the ITS has been referred to as a mosaic sequence , . High intra- and interspecies homology of conserved blocks has led to the conclusion that the primary source of ITS variation is homologous recombination , .
ITS typing has since been applied successfully to a number of bacterial systematic inquiries, e.g., subtyping of Clostridium difficile , Prochlorococcus and Synechococcus , and Acinetobacter calcoaceticus-A. baumannii complex . Widely applied elsewhere, this technique has only been minimally investigated in Burkholderia pseudomallei , , and only for the purposes of clinical diagnostics.
B. pseudomallei is the causative agent of melioidosis. Melioidosis is primarily acquired environmentally via inhalation, percutaneous inoculation or ingestion, and is associated with a variety of clinical outcomes, ranging from latent or subclinical infection to pneumonia, which often progresses to systemic sepsis and death . The disease is endemic in tropical regions of northern Australia and northeastern Thailand and is associated with mortality rates of 20% and 50%, respectively , . Advances in B. pseudomallei diagnostics indicate greater geographic distribution than previously thought, and the bacterium has been categorized as sporadic in other tropical regions, such as in Africa and South America . Additionally, the bacterium is categorized by the Centers for Disease Control and Prevention (CDC) as a Category B Select Agent, because of its amenability to deliberate large-scale dissemination and high morbidity and mortality .
A high degree of genomic variability exists between strains, complicating genetic and phylogenetic studies. The B. pseudomallei genome is composed of ∼7.2 Mbp on two chromosomes and exhibits a relatively high GC-content of 68% . The four completed B. pseudomallei whole genome sequences (strains K96243, 1710b, 1106a, and MSHR668) had four rrn operons. This genome is described as an “open genome,” with variability attributed mostly to genomic islands and islets, paralogous genes and indels , . Such variation is mediated by the interaction of insertion sequence elements and DNA mobilizing genes such as integrases and recombinases, or by horizontal gene transfer and subsequent homologous recombination , , . B. pseudomallei has been identified as one of the most genetically recombined bacterial species . Because of this high level of variation, finding reliable molecular targets to identify strains and their relationships remains difficult. High rates of recombination have complicated multi-locus sequence typing (MLST) and multi-locus variable number tandem repeat analysis (MLVA) schemes in the species ,  and may also be responsible for ribotype variability observed in restriction fragment length polymorphism (RFLP)  and pulsed-field gel electrophoresis (PFGE)  studies. In contrast, sequence analysis of the 16S rRNA gene showed nearly 100% conservation, often differing by only one single nucleotide polymorphism, making it difficult to resolve isolates using this locus .
The objective of this study was to elucidate phylogenetic relationship within B. pseudomallei and among closely related species using a high-throughput ITS typing scheme across a large, geographically and epidemiologically diverse strain panel. This adds to the scientific understanding of bacterial ITS diversity, as well as the ever-increasing knowledge of B. pseudomallei diversity, and may be used in conjunction with other typing studies to inform accurate phylogenetic conclusions within the species. Here, we also report that ITS typing can be complementary to MLST, and used in phylogenetic analysis of B. pseudomallei and its near neighbors.
NCBI's BLASTN algorithm was used in order to discover variability in ITS size and in rrn operon copy number in many Burkholderia species' whole genome sequences (Fig. 1; Table S1). Variability in rrn genomic location was assessed in ten B. pseudomallei genomes, in both completed genomes and whole genome shotgun sequences (Table 1) that had been assembled using the Artemis Comparison Tool . To discover the causative elements of length polymorphisms between ITS types, ITS sequences were aligned using the ClustalW multiple alignment method in BioEdit . To deduce ITS sequence identity among types, a sequence distance matrix was generated using MegAlign software in the Lasergene suite version 8.0.2 (DNASTAR, Inc), from which average pair-wise sequence similarity values (d-values) were produced. To deduce sequence identity between types, non-variable regions from each ITS type were extracted, concatenated, and average pair-wise similarities calculated as above. For generation of the ITS alignment, consensus ITS sequences for each type were generated using BioEdit, and subsequently aligned using the slow/accurate multiple alignment tool in the CLC Sequence Viewer version 6.0 (CLC bio A/S).
Panel A shows diagram of bacterial rRNA operon in Burkholderia species, demonstrating the specific location of the 16S-23S ITS. Panel B shows low resolution genomic alignment generated from consensus sequences of known ITS sequences, to highlight 1) sequence similarity among ITS types and 2) regions of difference between ITS types. Sequence labels include ITS code followed by species designation as follows: Bp (B. pseudomallei), Bm (B. mallei), Bt (B. thailandensis), Bt-like (B. thailandensis-like, referring to the reference strain MSMB43) and Bok (B. oklahomensis). ITS sizes are described at the right. Black bars at the bottom highlight major regions of variability I, II, III and IV that contribute to ITS fragment length polymorphisms, which are further demonstrated in the sequence conservation plot. Note that peaks outside of regions I, II, III and IV correspond to rare frame-shift and single nucleotide mutations, and are primarily the result of alignment with B. oklahomensis, whose ITS sequence is most divergent. Dotted lines represent missing sequences.
DNA samples and bacterial strains
DNA samples were obtained from different sources (see Table S2) and were prepared by several different methods. Essentially any genomic DNA preparation sufficient in yield and quality was interrogated. All putative B. pseudomallei isolates were cultured on selective Ashdown agar prior to DNA extraction, and were further confirmed as B. pseudomallei using species-specific SYBR-Green (Applied Biosystems, Inc.) real-time PCR assays, including the TTS1  and the multiplex assay targeting either the B. thailandensis-like flagellar gene cluster (BTFC) or Yersinia-like fimbrial (YLF) gene cluster . The identity of B. thailandensis isolates used were confirmed by targeting the cheB gene in the B. thailandensis homolog of BTFC , and B. mallei isolates were confirmed by targeting the bimA gene . For the B. thailandensis-like MSMB43 strains and B. oklahomensis strains used, characterization and species confirmation was performed at the United States CDC laboratory and included biochemical testing, 16S rDNA sequencing, MLST typing and DNA/DNA hybridization , .
To interrogate ITS fragments sizes in many Burkholderia species, PCR primers were designed to conserved regions of the 16S and 23S rRNA genes which flank the ITS. The ITS was amplified using approximately 500 picograms of template DNA in 10 µL reactions of the HotStarTaq Master Mix (QIAGEN) and the PCR primers GTW_For (5′- GTGAAGTCGTAACAAGGTAGCCGT-3′) and, to facilitate capillary electrophoresis, a 6-carboxyfluorescein labeled reverse primer, GTW_Rev_FAM (5′-/56-FAM/ AAGGCATCCACCACATGCACTT-3′) at a final concentration of 0.2 µM each. Thermal cycling was performed using an MJ Research PTC-0200 DNA Engine (BioRad) by the following parameters: 95°C for 15 min, followed by 35 cycles of 95°C for 30 sec, 62°C for 30 sec, 72°C for 1 min and a final extension of 72°C for 10 min.
To prepare PCR products for capillary electrophoresis, amplicons were diluted 1∶150 in Type I water and 1 µL each was denatured in reactions of 13.875 µL Hi-Di Formamide (Applied Biosystems, Inc.) and 0.125 µL GeneScan – 1200 LIZ Size Standard (Applied Biosystems, Inc.) at 95°C for 4 min. Products were then electrophoresed on an Applied Biosystems 3730xl DNA Analyzer, and ITS fragment sizes were auto-analyzed and binned using GeneMapper Software Version 4.0 (Applied Biosystems, Inc.).
To facilitate rapid screening of B. pseudomallei ITS types, TaqMan real-time PCR assays were developed and validated against a subset of those isolates screened (Text S1).
We used chi-square tests of association to examine relationships of ITS incidence in two endemic populations, Australia and Thailand. Three degrees of freedom were used in the analysis.
Twenty-two representative isolates (Table S1) were selected for ITS sequencing in order to confirm 1) sequence identity among identical ITS fragment sizes and 2) sequence conservation regardless of the geographic origin of the allele or its presence in a mixed array. For sequencing, PCR of the ITS was performed as before with the exception that an unlabelled reverse primer GTW_Rev (5′-AAGGCATCCACCACATGCACTT-3′) was used. PCR products were purified using ExoSAP-IT according to the manufacturer's protocol, with the difference that half of the recommended solution was used (USB Corporation; Affymetrix, Inc.). Cycle sequencing was carried out using the Big Dye Terminator (BDT) version 3.1 kit (Applied Biosystems, Inc.), according to the manufacturer's protocol. For isolates that possessed a single ITS type, direct cycle sequencing was performed with purified amplicons as template DNA.
GenBank ITS sequence submission
In order to assess rrn genomic distribution, we compared genomes from nine B. pseudomallei strains to the reference strain K96243, and found rrn genomic locations to be conserved (Table 1). We found four homologous rrn operons in all completed Burkholderia genomes analyzed although in some unfinished B. pseudomallei genomes, only three loci were observed. This may be due to incomplete genomic coverage or imperfect assembly. The only differences in location were explained by rare large-scale genomic rearrangement events (Table 1).
In silico analysis of various Burkholderia genomes revealed variation in the size of ITS regions within and among species. To facilitate comparisons, ITS fragment sizes were assigned alpha codes (Table 2). When searched, four completed B. pseudomallei genomes were found to possess one of two ITS types (type C or E) or had a mixture of the two ITS types distributed on the four rrn operons (type CE).
ITS types were interspecific in B. pseudomallei near-neighbors. All four B. mallei genomes possessed ITS type C that was virtually identical to the B. pseudomallei homolog. Genomes of B. thailandensis E264, B. thailandensis-like strain MSMB43, and B. oklahomensis strains C6786 and E0147 possessed an ITS fragment size distinct from any other Burkholderia species analyzed, designated as ITS types A, B and D, respectively (Table 2).
PCR screening of Burkholderia ITS types confirmed trends observed through in silico analysis of whole genome sequences and facilitated the discovery of a third major B. pseudomallei ITS type (ITS type G). Comprehensive typing data is presented in Table S2. All tested B. pseudomallei isolates possessed either a single homogeneous ITS type (either C, E, or G) or a heterogeneous mixture of two ITS types arrayed in unknown orientation in their four rrn operons (Table 2). The predominant ITS types found in 1,191 tested B. pseudomallei strains were C and E, accounting for 30.31% and 43.91% of ITS types, respectively. ITS type G represented a third major type of B. pseudomallei and was rare, accounting for 1.26% of isolates tested. Pair-wise mixtures were observed of all single ITS types, including CE (24.27%), GC (0.08%) and GE (0.17%).
Additionally, ITS typing confirmed species classifications of B. pseudomallei phylogenetic near-neighbors (Table S2). ITS types were not species specific based upon testing B. thailandensis (n = 88), B. mallei (n = 47), B. oklahomensis (n = 4) and Australian B. thailandensis-like (n = 2) strains. All B. mallei strains possessed only ITS type C, a type common to B. pseudomallei. Two strains of B. cepacia (Table S2) were found to possess ITS type F, but these are likely clonal members of the Burkholderia cepacia complex . Other members of the complex were screened and a high degree of diversity was found (data not shown). Therefore, allele F should not be considered as a marker for the diverse “Pseudomallei” phylogenetic group.
Worldwide distribution of B. pseudomallei ITS types and other correlates
B. pseudomallei ITS types correlated with their geographic origins (Fig. 2). In endemic populations such as northern Australia and Southeast Asia, ITS types C, E and hybrid CE dominated, whereas the less common type G was associated with isolates found in sporadic melioidosis regions such as Africa and South America. In fact, eight out of nine samples isolated from non-endemic locations possessed a G allele. The G allele was rarely observed in endemic populations, occurring 0.88% and 0.68% in Australian and Thai isolates tested, respectively. Interestingly, isolates possessing type CE were much more common in Thailand compared to Australia, whereas ITS type E was greatest in the Australian population (Fig. 3).
ITS types are color coded. Pie chart sizes reflect the number of isolates characterized. Samples taken from sporadic melioidosis regions, as in Africa, Central and South America, frequently possessed ITS type G, shown in purple. Endemic melioidosis areas show a high degree of ITS diversity, but possessed very few isolates with the G allele (0.88% of Australian isolates and 0.68% of Thai isolates).
ITS types C, E and CE were observed most frequently in these regions, while other ITS types were rarely observed. ITS type CE was observed 24% more often in Thailand as compared to Australia.
Other ITS correlations were investigated. For example, many isolates that were ITS typed had also been MLST typed, and single sequence types (STs) were rarely composed of single ITS types (Table S2). In some cases, however, MLST sequence types had conserved ITS types, as with the dominant Australian sequence types ST109 and ST132, which exclusively possessed ITS types E and C, respectively. Such a result was rarely seen in Thailand.
Epidemiological correlations with ITS type were investigated, but no consistent pattern was observed when stratifying isolates by their sources (clinical versus environmental). This may be due to the recombining nature of this locus and a lack of cause-and-effect relationships between virulence and these rrn alleles (data not shown). It is likely that ITS alleles are frequently mobilized in the environment, therefore confounding epidemiological and phylogenetic investigations using this particular marker.
ITS sequence analysis
DNA sequencing of representatives from each ITS type (Table S1) confirmed the near identical DNA sequence similarity among ITS types regardless of the strain's geographic origin or mixed ITS status. Pair-wise sequence similarities (d) among each ITS type ranged from d = 0.994 (ITS type G) to d = 0.999 (ITS type C) (Table 2). Alignment of consensus sequences from each ITS type by the multiple alignment method revealed high sequence similarity between ITS types (d = 0.97) with the exception of four variable regions (Var. I –IV), which are the sole contributors to ITS fragment length polymorphism (Fig. 1B).
Burkholderia species differed in specific ways within these variable regions (Fig. 1B). Variable regions I, II and III represent variation between different Burkholderia species, whereas variable region IV represents variation within B. pseudomallei. For example, all B. pseudomallei possessed conserved sequences for Var. I – II, lacked Var. III, and Var. IV variation was due to presence, absence, or length variation. However, the B. thailandensis ITS type A lacked sequences for all variable regions except Var. IV. ITS type B in Australian B. thailandensis-like MSMB43 appeared to be very similar to B. thailandensis, except that the lengths of missing sequence in Var. I and II were less than in ITS type A. The B. oklahomensis ITS type D was most divergent, lacking Var. I sequence, but possessing Var. II and III sequences, with an additional small deletion in Var. IV. Note that Var. IV sequence in B. oklahomensis was divergent from Var. IV sequence in B. thailandensis, B. thailandensis-like MSMB43 or B. pseudomallei (data not observable in Fig. 1).
Insights from Burkholderia ITS typing
In this study, we sought to characterize the distribution of ITS differences in B. pseudomallei and its near-neighbors. Our ITS data indicate the existence of three major ITS types of B. pseudomallei. A high rate of recombination , ,  together with mobilizing genetic elements ,  has made the B. pseudomallei genome extremely diverse, and thus phylogenetic assignments are complicated by homoplasy. Standard 16S rDNA typing has been useful in identifying the genus and species, but the gene is too conserved to resolve B. pseudomallei populations . ITS polymorphisms observed within B. pseudomallei appear to be a result of genetic recombination. Therefore, we hypothesize that single ITS types represent ancestral states on which genomic diversity has been built in the species.
Strains with mixed ITS types are most likely the result of contemporary recombination events between strains rather than being generated within a single bacterium by current in situ mutation events. ITS mobilization, a lateral gene transfer, has been hypothesized in other bacteria, and has been associated with homologous recombination between identical sequence blocks , , including potentially conserved tRNA genes. In fact, tRNA genes serve as “hot spots” for site-specific recombination events in other loci . The species itself is prone to recombination, as found in MLST studies in which new sequence types were more commonly the result of mixing of existing alleles rather than intragenic mutations , . The homogenous nature of lesser-recombining B. pseudomallei near-neighbor ITS types also supports the recombination hypothesis (Table 2).
ITS typing and molecular epidemiology
ITS mixed types were found significantly more often in Thailand than in Australia, as evidenced by the 24% greater incidence of ITS type CE (Fig. 3). This observed high incidence of mixed types is consistent with a study suggesting higher rates of recombination in relation to mutation in Southeast Asian strains . Increased soil bacterial loads likely increase the rate of contact and therefore genomic sharing between strains, and Thai soil has been found to harbor significantly more viable B. pseudomallei than Australian soil , , . Additionally, B. pseudomallei incidence has been linked with soil disturbance and agriculture , which are commonplace activities in the Thai rice farming industry. Such ITS mixing may help to explain the low levels of association between ITS types and MLST types, especially in the Thai population.
ITS types C, E and CE represent most endemic isolates, whereas type G encompasses less than 1% of Thai or Australian isolates but is the primary type in sporadic isolates (Fig. 2). It remains unclear whether the former are simply more fit than the latter to survive in endemic regions, or whether a sampling bias exists in our isolate collection. Further ecological and functional studies, in conjunction with more thorough sampling, will be needed to elucidate this. Less clear is why sporadic B. pseudomallei isolates almost always possessed ITS type G. Although our sample size of sporadic isolates was small (n = 9), sampling bias is unlikely because of the high degree of association, even from multiple locations. This dispersal pattern appears to be a genetic bottleneck of unknown cause.
Dispersal by human trade is one explanation for the disparate distribution of B. pseudomallei ITS type G. Because genomic evidence indicates Australia as the founding population , , , , worldwide seeding from this or another endemic area is plausible. Indeed, the G allele is present in the Australian population, albeit minimally. The genotype may have accessed sporadic melioidosis locations via multiple exports from the same location, or via a single exporter that travelled to several ports worldwide. Trading of potted tropical plants may have served as the perfect temporary habitat for B. pseudomallei. If true, it is possible that many ITS types were exported, and that ITS type G was most fit to survive in sporadic locations. Such isolates, then, would be less fit to survive in endemic areas, which is consistent with the low incidence of allele G in Australia and Southeast Asia. Surely, comparative genomic and functional studies between endemic and non-endemic isolates must be performed in order to answer questions about bacterial fitness within and outside endemic areas.
Conservation of the ITS sequences
ITS sequence conservation is not random (Fig. 1B). Among individual ITS types, sequence identity neared 100% in all cases (Table 2). This finding makes ITS sequencing unnecessary in future studies, unless novel ITS types are discovered. Such high levels of conservation indicate selection and a functional role. In support of this hypothesis, other studies have characterized clear roles of the ITS in rRNA processing, maturation and transcriptional anti-termination , , . Additionally, tRNA gene sequences within the ITS were conserved (Fig. 1; tRNA-Ile at bp 83..142; tRNA-Ala at bp 194..253). Because sequence between ITS types was conserved when variable regions I-IV were neglected, it is likely that these conserved regions contain the functional motifs and that variable regions are extraneous. It is also possible that variable regions do serve functional roles and influence the efficiency of transcription, affecting environmental fitness and thus explaining observed bacterial niches. Indeed, rrn copy number has been associated with bacterial fitness , , and so it is plausible that the ITS may also play a role. Further structural and functional studies will be necessary to confirm functional roles of the Burkholderia ITS.
In addition, ITS mutational patterns can be traced when analyzed in light of a Burkholderia MLST phylogeny (Fig. 4). Interspecific ITS length differences between the ITS types are attributed to specific insertion and deletion events, possibly as the result of homologous recombination between strains, which may be evolutionarily traced when compared to this phylogeny. Variable region I, for example, did not appear in Burkholderia species until B. pseudomallei speciated from a common ancestor. It is therefore likely that all B. pseudomallei ancestors lacked Var. I. An alternative hypothesis that a basal Burkholderia ancestor possessed the region, which was lost in three speciation events but kept in B. pseudomallei, is considerably less probable.
The presence (+), absence (-) or partial absence (indicated by the Δ symbol) of variable ITS sequence blocks are described for major species-level Burkholderia clades, based on MLST sequence analysis. Because B. mallei is a clone of B. pseudomallei and possesses an identical ITS type, its clade was collapsed into that of B. pseudomallei.
The most parsimonious evolutionary hypothesis for Var. II is one in which a clonal ancestor possessed sequence at that locus, which was differentially partially deleted in three speciation events: MSMB43 (Δ17 bp), B. thailandensis (Δ21 bp) and B. pseudomallei (Δ4 bp). This hypothesis is in contrast to one that suggests that the basal ancestor lacked sequence, and that identical sequence was inserted in four distinct speciation events.
ITS and the Pseudomallei group's evolutionary patterns
ITS typing supported confirmed species classifications among species in the “Pseudomallei group” and provides further support of evolutionary hypotheses. The Pseudomallei group is one defined by NCBI's taxonomy database, and specifically includes the species studied here (http://www.ncbi.nlm.nih.gov/Taxonomy/). B. mallei, for example, is an obligate pathogen and has been called a clone of B. pseudomallei, said to have become specialized through genome reduction to occupy only the intracellular niche, rather than soil niches , . All B. mallei possessed ITS type C, which was identical to the B. pseudomallei ortholog, supporting the species as a clone of B. pseudomallei, and more specifically, of the basal group of B. pseudomallei that possessed only ITS type C. This species showed no variation in the ITS type and had no incidence of ITS mixing (d = 1.00 for B. mallei's ITS type C), a result consistent with our previous hypothesis that inter-strain contact and subsequent horizontal gene transfer were the sources of ITS mixing. In other words, the bacterium rarely, if ever, encounters suitable gene donors with which to recombine. Because of this, B. mallei may be considered a clonal bacterium, while B. pseudomallei displays high levels of homoplastic evolution and thus cannot be considered clonal.
B. thailandensis, on the other hand, is more puzzling. A close genetic relative and soil commensal of B. pseudomallei, all B. thailandensis possessed ITS type A, supporting its confirmed species status , . Like B. mallei, the species was completely monomorphic with regard to ITS allele, but in contrast, the bacterium certainly encounters suitable gene donors in the environment. Such aplasticity in the ITS may indicate greater overall genomic stability compared to B. pseudomallei, which had three possible ITS states and mixing therein. Further comparative studies in a larger B. thailandensis population will be necessary to confirm this hypothesis.
MSMB43, an Australian soil bacterial strain that had been described as similar to B. thailandensis, by 16S rDNA sequencing, MLST typing and DNA/DNA hybridization , and was found to possess ITS type B in silico. ITS sequence comparisons confirmed this relationship, as the two species possess very similar ITS types (Fig. 1B). When screened, one other isolate possessed ITS type B, MSHR1554, a strain whose MLST sequence type was found to group phylogenetically with MSMB43 . The same result was confirmed in B. oklahomensis , a pathogenic soil bacterium found in the southern USA, in which all putative strains screened were found to have ITS type D. Thus, at least in non-B. pseudomallei isolates, ITS typing is consistent with MLST typing (Fig. 4). Because ITS type E possessed sequence at Var. IV that is most similar to the ITS types of B. thailandensis and B. thailandensis-like MSMB43 (Fig 1), we hypothesize that ITS type E represents the most ancient B. pseudomallei group. This inference is consistent with the result that ITS type E was most prevalent in the Australian population, which is suspected to be the founding population for all B. pseudomallei . From ITS type E, two derivatives evolved: C and G. It is impossible to determine whether each evolved independently from ITS type E, or whether one evolved from the other as each scenario is equally parsimonious. If the latter hypothesis is correct, then type G would have evolved from type C.
In summary, this study fills an important knowledge gap in Burkholderia genomics and confirms trends observed in other research. By applying the widely-used technique of 16S-23S ITS typing to Burkholderia species we were able to demonstrate several important points. We were able to show that ancient B. pseudomallei was composed of three genomic groups, represented by ITS types C, E and G, and hypothesize that ITS type E is the most ancestral type. This research demonstrates that B. pseudomallei is highly variable in comparison to other species in the Pseudomallei group. Moreover, we demonstrate a greater level of recombination in Thai B. pseudomallei isolates compared to isolates from other locations, which agrees with previous findings. Regarding geographic distribution, we identified one B. pseudomallei group (represented by ITS type G) that correlated strongly with worldwide sporadic distribution of melioidosis, and hypothesize human travel and trade as the mechanism of this dispersal pattern. With regard to phylogenetic near-neighbors, we were able to support B. mallei as a clone of B. pseudomallei. We show that other species are evolutionarily distinct and display less genomic plasticity than B. pseudomallei. Increased sampling, in conjunction with comparative functional and ecological studies, will be necessary to confirm hypotheses that we present here.
List of ITS sequences or whole-genome sources of ITS sequences used for bioinformatics analysis in this study.
All ITS typing data generated in this study, including isolate contributor, geographic origin, epidemiological data and MLST data, when known.
Validation of PCR-ribotyping using TaqMan real-time PCR assays.
Thank you to Drs. David DeShazer, Steven P. Harvey, Abeera Naureen and Muhammad Saqib for providing DNA samples. Thank you to Linda Ward for providing epidemiological data on strains used in this study.
Conceived and designed the experiments: AT. Performed the experiments: APL SDW JLG JB MM DE. Analyzed the data: APL TP DMW AT. Contributed reagents/materials/analysis tools: MBG VW SJP. Wrote the paper: APL AT BJC DMW SJP PK.
- 1. Kostman JR, Edlind TD, LiPuma JJ, Stull TL (1992) Molecular epidemiology of Pseudomonas cepacia determined by polymerase chain reaction ribotyping. J Clin Microbiol 30: 2084–2087.
- 2. Fisher MM, Triplett EW (1999) Automated approach for ribosomal intergenic spacer analysis of microbial diversity and its application to freshwater bacterial communities. Appl Environ Microbiol 65: 4630–4636.
- 3. Nielsen XC, Justesen US, Dargis R, Kemp M, Christensen JJ (2009) Identification of Clinically Relevant Nonhemolytic Streptococci on the Basis of Sequence Analysis of 16S-23S Intergenic Spacer Region and Partial gdh Gene. J Clin Microbiol 47: 932–939.
- 4. Rachman C, Kabadjova P, Valcheva R, Prevost H, Dousset X (2004) Identification of Carnobacterium species by restriction fragment length polymorphism of the 16S-23S rRNA gene intergenic spacer region and species-specific PCR. Appl Environ Microbiol 70: 4468–4477.
- 5. Apirion D, Miczak A (1993) RNA processing in prokaryotic cells. BioEssays 15: 113–120.
- 6. Berg KL, Squires C, Squires CL (1989) Ribosomal RNA operon anti-termination. Function of leader and spacer region box B-box A sequences and their conservation in diverse micro-organisms. J Mol Biol 209: 345–358.
- 7. Rocap G, Distel DL, Waterbury JB, Chisholm SW (2002) Resolution of Prochlorococcus and Synechococcus ecotypes by using 16S-23S ribosomal DNA internal transcribed spacer sequences. Appl Environ Microbiol 68: 1180–1191.
- 8. Acinas SG, Marcelino LA, Klepac-Ceraj V, Polz MF (2004) Divergence and Redundancy of 16S rRNA Sequences in Genomes with Multiple rrn Operons. J Bacteriol 186: 2629–2635.
- 9. Gürtler V (1999) The role of recombination and mutation in 16S-23S rDNA spacer rearrangements. Gene 238: 241–252.
- 10. Tambong JT, Xu R, Bromfield ES (2009) Intercistronic heterogeneity of the 16S-23S rRNA spacer region among Pseudomonas strains isolated from subterranean seeds of Hog Peanut (Amphicarpa bracteata). Microbiology.
- 11. Klappenbach JA, Dunbar JM, Schmidt TM (2000) rRNA operon copy number reflects ecological strategies of bacteria. Appl Environ Microbiol 66: 1328–1333.
- 12. Sadeghifard N, Gurtler V, Beer M, Seviour RJ (2006) The mosaic nature of intergenic 16S-23S rRNA spacer regions suggests rRNA operon copy number variation in Clostridium difficile strains. Appl Environ Microbiol 72: 7311–7323.
- 13. Lan R, Reeves PR (1998) Recombination between rRNA operons created most of the ribotype variation observed in the seventh pandemic clone of Vibrio cholerae. Microbiology 144(Pt 5): 1213–1221.
- 14. Chang HC, Wei YF, Dijkshoorn L, Vaneechoutte M, Tang CT, et al. (2005) Species-level identification of isolates of the Acinetobacter calcoaceticus-Acinetobacter baumannii complex by sequence analysis of the 16S-23S rRNA gene spacer region. J Clin Microbiol 43: 1632–1639.
- 15. Tyler SD, Strathdee CA, Rozee KR, Johnson WM (1995) Oligonucleotide primers designed to differentiate pathogenic pseudomonads on the basis of the sequencing of genes coding for 16S-23S rRNA internal transcribed spacers. Clin Diagn Lab Immunol 2: 448–453.
- 16. Kunakorn M, Markham RB (1995) Clinically practical seminested PCR for Burkholderia pseudomallei quantitated by enzyme immunoassay with and without solution hybridization. J Clin Microbiol 33: 2131–2135.
- 17. Cheng AC, Currie BJ (2005) Melioidosis: epidemiology, pathophysiology, and management. Clin Microbiol Rev 18: 383–416.
- 18. Cheng AC, Hanna JN, Norton R, Hills SL, Davis J, et al. (2003) Melioidosis in northern Australia, 2001-02. Commun Dis Intell 27: 272–277.
- 19. White NJ (2003) Melioidosis. Lancet 361: 1715–1722.
- 20. Currie BJ, Dance DA, Cheng AC (2008) The global distribution of Burkholderia pseudomallei and melioidosis: an update. Trans R Soc Trop Med Hyg 102: Suppl 1S1–4.
- 21. Rotz LD, Khan AS, Lillibridge SR, Ostroff SM, Hughes JM (2002) Public health assessment of potential biological terrorism agents. Emerg Infect Dis 8: 225–230.
- 22. Holden MT, Titball RW, Peacock SJ, Cerdeno-Tarraga AM, Atkins T, et al. (2004) Genomic plasticity of the causative agent of melioidosis, Burkholderia pseudomallei. Proc Natl Acad Sci U S A 101: 14240–14245.
- 23. Sim SH, Yu Y, Lin CH, Karuturi RK, Wuthiekanun V, et al. (2008) The core and accessory genomes of Burkholderia pseudomallei: implications for human melioidosis. PLoS Pathog 4: e1000178.
- 24. Tuanyok A, Leadem BR, Auerbach RK, Beckstrom-Sternberg SM, Beckstrom-Sternberg JS, et al. (2008) Genomic islands from five strains of Burkholderia pseudomallei. BMC Genomics 9: 566.
- 25. Tumapa S, Holden MT, Vesaratchavest M, Wuthiekanun V, Limmathurotsakul D, et al. (2008) Burkholderia pseudomallei genome plasticity associated with genomic island variation. BMC Genomics 9: 190.
- 26. Pearson T, Giffard P, Beckstrom-Sternberg S, Auerbach R, Hornstra H, et al. (2009) Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer. BMC Biology 7: 14.
- 27. Pearson T, U'Ren JM, Schupp JM, Allan GJ, Foster PG, et al. (2007) VNTR analysis of selected outbreaks of Burkholderia pseudomallei in Australia. Infect Genet Evol 7: 416–423.
- 28. Vesaratchavest M, Tumapa S, Day NP, Wuthiekanun V, Chierakul W, et al. (2006) Nonrandom distribution of Burkholderia pseudomallei clones in relation to geographical location and virulence. J Clin Microbiol 44: 2553–2557.
- 29. Sermswan RW, Wongratanacheewin S, Trakulsomboon S, Thamlikitkul V (2001) Ribotyping of Burkholderia pseudomallei from clinical and soil isolates in Thailand. Acta Trop 80: 237–244.
- 30. Pitt TL, Trakulsomboon S, Dance DA (2000) Molecular phylogeny of Burkholderia pseudomallei. Acta Trop 74: 181–185.
- 31. Gee JE, Sacchi CT, Glass MB, De BK, Weyant RS, et al. (2003) Use of 16S rRNA gene sequencing for rapid identification and differentiation of Burkholderia pseudomallei and B. mallei. J Clin Microbiol 41: 4647–4654.
- 32. Carver TJ, Rutherford KM, Berriman M, Rajandream MA, Barrell BG, et al. (2005) ACT: the Artemis Comparison Tool. Bioinformatics 21: 3422–3423.
- 33. Hall TA (1999) BioEdit: a user-friendly biological sequence alignment editor and analysis program for Windows 95/98/NT. Nucl Acids Symp Ser 41: 95–98.
- 34. Novak RT, Glass MB, Gee JE, Gal D, Mayo MJ, et al. (2006) Development and evaluation of a real-time PCR assay targeting the type III secretion system of Burkholderia pseudomallei. J Clin Microbiol 44: 85–90.
- 35. Tuanyok A, Auerbach RK, Brettin TS, Bruce DC, Munk AC, et al. (2007) A Horizontal Gene Transfer Event Defines Two Distinct Groups within Burkholderia pseudomallei That Have Dissimilar Geographic Distributions. J Bacteriol 189: 9044–9049.
- 36. Ulrich RL, Ulrich MP, Schell MA, Kim HS, DeShazer D (2006) Development of a polymerase chain reaction assay for the specific identification of Burkholderia mallei and differentiation from Burkholderia pseudomallei and other closely related Burkholderiaceae. Diagn Microbiol Infect Dis 55: 37–45.
- 37. Gee JE, Glass MB, Novak RT, Gal D, Mayo MJ, et al. (2008) Recovery of a Burkholderia thailandensis-like isolate from an Australian water source. BMC Microbiol 8: 54.
- 38. Glass MB, Steigerwalt AG, Jordan JG, Wilkins PP, Gee JE (2006) Burkholderia oklahomensis sp. nov., a Burkholderia pseudomallei-like species formerly known as the Oklahoma strain of Pseudomonas pseudomallei. Int J Syst Evol Microbiol 56: 2171–2176.
- 39. Mahenthiralingam E, Baldwin A, Dowson CG (2008) Burkholderia cepacia complex bacteria: opportunistic pathogens with important natural biology. J Appl Microbiol 104: 1539–1551.
- 40. Cheng AC, Ward L, Godoy D, Norton R, Mayo M, et al. (2008) Genetic diversity of Burkholderia pseudomallei isolates in Australia. J Clin Microbiol 46: 249–254.
- 41. Kaestli M, Mayo M, Harrington G, Ward L, Watt F, et al. (2009) Landscape Changes Influence the Occurrence of the Melioidosis Bacterium Burkholderia pseudomallei in Soil in Northern Australia. PLoS Negl Trop Dis 3: e364.
- 42. Smith MD, Wuthiekanun V, Walsh AL, White NJ (1995) Quantitative recovery of Burkholderia pseudomallei from soil in Thailand. Trans R Soc Trop Med Hyg 89: 488–490.
- 43. Wuthiekanun V, Smith MD, Dance DA, White NJ (1995) Isolation of Pseudomonas pseudomallei from soil in north-eastern Thailand. Trans R Soc Trop Med Hyg 89: 41–43.
- 44. Currie BJ (2008) Advances and remaining uncertainties in the epidemiology of Burkholderia pseudomallei and melioidosis. Transactions of the Royal Society of Tropical Medicine and Hygiene 102: 225–227.
- 45. Shrestha PM, Noll M, Liesack W (2007) Phylogenetic identity, growth-response time and rRNA operon copy number of soil bacteria indicate different stages of community succession. Environmental Microbiology 9: 2464–2474.
- 46. Godoy D, Randle G, Simpson AJ, Aanensen DM, Pitt TL, et al. (2003) Multilocus sequence typing and evolutionary relationships among the causative agents of melioidosis and glanders, Burkholderia pseudomallei and Burkholderia mallei. J Clin Microbiol 41: 2068–2079.
- 47. Nierman WC, DeShazer D, Kim HS, Tettelin H, Nelson KE, et al. (2004) Structural flexibility in the Burkholderia mallei genome. Proc Natl Acad Sci U S A 101: 14246–14251.
- 48. Brett PJ, DeShazer D, Woods DE (1998) Burkholderia thailandensis sp. nov., a Burkholderia pseudomallei-like species. Int J Syst Bacteriol 48(Pt 1): 317–320.