Burkholderia pseudomallei genome plasticity associated with genomic island variation
© Tumapa et al; licensee BioMed Central Ltd. 2008
Received: 01 October 2007
Accepted: 25 April 2008
Published: 25 April 2008
Burkholderia pseudomallei is a soil-dwelling saprophyte and the cause of melioidosis. Horizontal gene transfer contributes to the genetic diversity of this pathogen and may be an important determinant of virulence potential. The genome contains genomic island (GI) regions that encode a broad array of functions. Although there is some evidence for the variable distribution of genomic islands in B. pseudomallei isolates, little is known about the extent of variation between related strains or their association with disease or environmental survival.
Five islands from B. pseudomallei strain K96243 were chosen as representatives of different types of genomic islands present in this strain, and their presence investigated in other B. pseudomallei. In silico analysis of 10 B. pseudomallei genome sequences provided evidence for the variable presence of these regions, together with micro-evolutionary changes that generate GI diversity. The diversity of GIs in 186 isolates from NE Thailand (83 environmental and 103 clinical isolates) was investigated using multiplex PCR screening. The proportion of all isolates positive by PCR ranged from 12% for a prophage-like island (GI 9), to 76% for a metabolic island (GI 16). The presence of each of the five GIs did not differ between environmental and disease-associated isolates (p > 0.05 for all five islands). The cumulative number of GIs per isolate for the 186 isolates ranged from 0 to 5 (median 2, IQR 1 to 3). The distribution of cumulative GI number did not differ between environmental and disease-associated isolates (p = 0.27). The presence of GIs was defined for the three largest clones in this collection (each defined as a single sequence type, ST, by multilocus sequence typing); these were ST 70 (n = 15 isolates), ST 54 (n = 11), and ST 167 (n = 9). The rapid loss and/or acquisition of gene islands was observed within individual clones. Comparisons were drawn between isolates obtained from the environment and from patients with melioidosis in order to examine the role of genomic islands in virulence and clinical associations. There was no reproducible association between the individual or cumulative presence of five GIs and a range of clinical features in 103 patients with melioidosis.
Horizontal gene transfer of mobile genetic elements can rapidly alter the gene repertoire of B. pseudomallei. This study confirms the utility of a range of approaches in defining the presence and significance of genomic variation in natural populations of B. pseudomallei.
Burkholderia pseudomallei is a soil dwelling Gram-negative bacterium and the cause of melioidosis, a serious human infection commonly reported in northeast Thailand and northern Australia [1, 2]. The genome of this biothreat agent contains genomic island (GI) regions that have properties indicative of horizontal gene transfer, i.e. DNA displaying anomalies in % G+C content or dinucleotide frequency signature, and/or the presence of coding sequences (CDSs) with similarities to genes associated with mobile genetic elements such as insertion sequence (IS) elements, bacteriophages and plasmids. Sixteen GIs were identified in the genome of B. pseudomallei strain K96243 , comprising approximately 6% of the genome. Screening of 40 clinical and environmental isolates of B. pseudomallei for the presence of selected K96243 GIs found that there was variable distribution , suggesting that horizontal gene transfer contributes to genome diversity. Genomic islands in other bacterial species encode many different functions, and selection may favour the maintenance of islands that increase fitness in a specialized environmental niche. Genomic islands also play a pivotal role in virulence of a large number of bacterial pathogens, carrying clusters of virulence genes encoding a wide range of functions including iron uptake systems, adhesins, superantigens, and genes that alter the antibiotic resistance phenotype .
B. pseudomallei is found in the soil and surface water of endemic areas and usually enters the host via cuts and abrasions in the skin or by inhalation following contact with contaminated water or soil. Once inside the body, the time taken for disease symptoms to appear can vary from hours to years. In some cases, B. pseudomallei infection does not cause overt disease around the time of exposure, but rather results in a state of bacterial dormancy in the host for many years before causing disease in later life. The most remarkable example is that of a second world war veteran who developed melioidosis 62 years after returning to the US from SE Asia .
The clinical symptoms of melioidosis are so varied that the disease has been termed "the great mimicker" . Disease manifestations are extremely wide ranging and vary from acute, fulminant sepsis to localized disease. The most frequent picture is a septicemic illness, often associated with bacterial dissemination to distant sites such as the liver and spleen . Pneumonia occurs in around 50% of patients. Other sites of infection include bone, joints, skin, soft tissue, prostate and the central nervous system. Overall mortality is around 50% in northeast Thailand (35% in children) , and 19% in Australia . The basis for marked variability in clinical presentation and disease severity is unknown. Host determinants such as diabetes mellitus and renal failure represent risk factors for infection with B. pseudomallei [1, 7], but no host factors have been associated with specific clinical features. The role of bacterial factors in determining disease variability and severity is poorly understood.
In this study we examined the distribution of a representative sample of five GIs previously described for B. pseudomallei K96243. In silico analysis of 10 genome sequences revealed the micro-evolutionary processes responsible for the short-term diversification of these regions. Multiplex PCR screening was used to investigate GI presence in B. pseudomallei isolates from NE Thailand. Rapid loss and/or acquisition of gene islands was observed by assaying the presence of islands within isolates belonging to single clones (as defined by multilocus sequence typing (MLST)) . Comparisons were drawn between isolates obtained from the environment and from patients with melioidosis in order to examine the role of genomic islands in virulence.
Results and Discussion
Overview of five genomic islands in B. pseudomallei strain K96243
The five genomic islands studied here include examples of a prophage (GI 2; φK96243) which has been shown to be mobile in this strain, two prophage-like islands (GI 6 and GI 9), a putative integrated plasmid (GI 11), and a putative metabolic island (GI 16) that lacks any obvious genes for mobilisation or integration . The five islands are a representative sample of the islands present in a single strain, and provide a snapshot of GI diversity. The mosaic structure and variable composition of GIs means that it is impractical to attempt to capture the full gamut of elements circulating in a large natural bacterial population, and we therefore focused our analysis on a representative sample of GIs identified in K96243.
A summary of the CDS content of each of the five GIs in B. pseudomallei K96243 is provided in Additional file 1, and a comprehensive list of the CDSs contained in these regions is provided in Additional File 2. The overwhelming majority of CDSs that are not phage or plasmid related are of unknown function. GIs 2, 6, 9 and 11 contain CDSs with no homology to known metabolic or putative virulence proteins. GIs 2, 6 and 11 contain low % GC regions relative to the rest of the island, indicating that the islands are mosaic and some genes may have been acquired recently; new gene acquisitions may have atypical sequence characteristics that are maintained over a relatively short time frame. GI 16 contains several CDSs with similarity to known virulence determinants, including a putative haemagglutinin and processing protein, and genes involved in the acquisition and utilization of nutrients. The abundance of CDSs on this island that encode functions that potentially broaden the metabolic repertoire of K96243 (for example possible sugar utilization and amino acid catabolism gene clusters and accompanying regulators), has led to this island being previously designated a putative metabolic island . CDSs in this island region exhibited varying levels of similarity to proteins from a range of taxonomically diverse organisms, and it was not possible to speculate on the likely source of this island region from similarity searches.
Comparative genomic analysis of genomic islands
Distribution of islands in sequenced B. pseudomallei strains
+/- (BPSS2060 to BPSS2076)
+/- (BPSS2057 to BPSS2076)
+/- (BPSS2057 to BPSS2076)
+/- (BPSS2057 to BPSS2076)
+/- (BPSS2057 to BPSS2076)
+/- (BPSS2057 to BPSS2076)
+/- (BPSS2057 to BPSS2076)
We examined the structure, composition and insertion site of the five GIs using the ten available genome sequences (complete and assemblies from unfinished genomes). GIs 2, 6, 9 and 11 all had flanking repeats, integrases, and mobile element-like genes, consistent with the suggestion that these were recently acquired mobile genetic elements. There were two instances where the GI was absent at the predicted site of insertion, but DNA with extended similarity to the island was detected elsewhere in the genome (one example each for GI 2 and GI 9). In B. pseudomallei strain 1655 a prophage with extended similarity to GI 2 (95–99% DNA identity over ~95% of the length of GI 2; data not shown), is integrated at an alternative tRNA gene on chromosome 1 (tRNA-Arg as opposed to tRNA-Phe). In the genome of K96243, GI 12 is found at this alternative locus. In 1106a an island with similarity to GI 2 is inserted at an alternative locus (tRNA-Ser) on chromosome 1. An identical island is present in 1106b. The 1106a island has a mosaic structure; approximately half of the island is highly similar (95–99% DNA identity) to the K96243 GI 9, with the other half exhibiting no similarity (see Additional file 3). The location of these sequences at alternative loci may be due to genomic rearrangement following divergence from a common ancestor carrying the island, or may reflect the presence of multiple insertion sites for related mobile genetic elements in the B. pseudomallei genome.
Presence of five genomic islands in B. pseudomallei
Detection by PCR of five genomic islands in 186 B. pseudomallei isolates from northeast Thailand
All isolates (n = 186)
Soil (n = 83)
Invasive (n = 103)
Cumulative number of genomic islands
Variability in genomic islands between clones of B. pseudomallei
Presence of genomic islands in the three largest clones as defined by multilocus sequence type
Number of isolates
These findings suggest that horizontal gene transfer can alter the gene repertoire ("accessory genome") of this species very rapidly, i.e. before any observable changes in the vertically inherited "core" genome, and is consistent with the hypothesis that isolates that are identical by MLST might exploit different ecological niches . The MLST data for B. pseudomallei has revealed a very high frequency of homologous recombination resulting in limited clustering and low levels of linkage disequilibrium [8, 11, 12]. Our results also suggest that the loss and acquisition of genomic islands by horizontal gene transfer significantly contributes to the dynamism of the genome. However, it should be noted that these two processes are distinct, and that a high rate of island transfer does not necessarily predict a high rate of homologous recombination or vice versa. For example, MLST data for the Gram-positive pathogen Staphylococcus aureus has revealed only modest rates of homologous recombination, yet the "accessory" genome of this species is known to change very rapidly [10, 13–16]. B. pseudomallei represents an alternative scenario, the data indicating both high rates of homologous recombination together with rapid changes in GI repertoire.
Comparison of the presence of genomic islands in soil and invasive isolates
We examined the hypothesis that the GIs play a role in the ability to cause human disease by testing whether the presence of GIs was significantly different in clinical and environmental isolates. We found that there was no significant difference in the presence of any of the GIs between these two groups (p > 0.05 for all five islands) (Table 2). We considered the possibility that combinations of GIs could alter the ability to cause disease through either a cumulative or synergistic effect. We defined the cumulative number of GIs per isolate for the 186 Thai isolates; this ranged from 0 to 5 (median 2, IQR 1 to 3) (Table 2). The distribution of cumulative GI number did not differ between environmental and disease-associated isolates (p = 0.27) (Table 2). These data indicate that none of the islands examined were numerically associated either individually or on a cumulative basis with the ability to cause melioidosis.
Relationship between clinical features and genomic islands in disease-associated isolates
The lack of association between the presence of the five GIs and the ability to cause human disease does not exclude the possibility that these GIs are associated with specific clinical features or outcome. To study this, we undertook a detailed statistical analysis of the relationship between a range of clinical features in 103 patients with melioidosis and the presence of GIs in their associated B. pseudomallei isolates. Univariate analysis was used to compare the presence of each GI and their cumulative number with the presence of blood cultures positive for B. pseudomallei, presence of ultrasound confirmed abscess(es) in liver and/or spleen, pneumonia, multiple organ involvement, the presence on admission or subsequent development of hypotension, impaired renal function, impaired liver function and death during hospital admission. GI 6 was negatively associated with positive blood culture (OR 0.30; 95%CI 0.11–0.82, p = 0.02), and GI 11 was negatively associated with impaired renal function (OR 0.36; 95%CI 0.14–0.94, p = 0.04). Cumulative number of GIs was negatively associated with impaired renal function (OR 0.54; 95%CI 0.33–0.86, p = 0.01). These significant associations could be a chance finding related to multiple comparisons. To test this hypothesis, the presence of GIs were determined for an independent set of 255 invasive isolates obtained from patients with melioidosis presenting to Sappasithiprasong Hospital between 2002–2003. There was no significant association between GI 6 and positive blood cultures (OR 1.12; 95%CI 0.51–2.43, p= 0.78), or GI 11 and impaired renal function (OR 0.67; 95%CI 0.34–1.32, p = 0.25). Cumulative number of GIs was not negatively associated with impaired renal function. We conclude that there were no significant associations between clinical features and any of the five GIs examined.
This survey of divergent putative GIs in a well-defined and clinically significant B. pseudomallei population has demonstrated that genetic exchange is widespread. Our in silico analysis provided evidence that GI 16 may represent a composite region rather than a genomic island, and defined a novel GI at the predicted insertion site for GI 11. A detailed analysis of the relationship between five GIs and both the ability to cause disease and a range of clinical features failed to show an association between the two. However, a previous study that utilized MLST to examine the bacterial genotype of 266 B. pseudomallei isolates (83 soil and 183 invasive) obtained from northeast Thailand reported that isolates from patients with melioidosis were significantly over-represented in the 10 largest clones . A classification index applied to examine differences in the frequency of STs, and the frequency of alleles at specific loci, also demonstrated a significant difference between the two bacterial populations. This indicates that bacterial genotypes were not distributed randomly between soil and invasive isolate groups, an observation that could be accounted for by variations in the accessory genome of some strains that leads to heightened virulence. Further evidence for this suggestion comes from a study in Australia, in which variability in genome content between two Australian B. pseudomallei isolates was utilized to define sequence variability and to develop a variable amplicon typing scheme designed to score the presence of 14 PCR amplicons in Thai and Australian isolates . Isolates clustered into groups, one of which was mostly associated with severe disease.
Our study confirms the utility of a range of approaches in defining the presence and significance of genomic variation in natural populations of B. pseudomallei. Further studies are required to determine the full extent of variability in GIs, and their relationship to biological fitness in the environment and to disease pathogenesis in the human host.
A total of 186 B. pseudomallei isolates were obtained from northeast Thailand. Of these, 83 were isolated from the environment in the province of Ubon Ratchathani, northeast Thailand between 1990 and 2003. All sites were flooded rice paddies, the majority of which were sampled after ploughing but before planting. Environmental samples were processed for the presence of B. pseudomallei, as previously described . A further 103 isolates were cultured from consecutive patients presenting with melioidosis during 2001 to Sappasithiprasong hospital, Ubon Ratchathani. Clinical manifestations of infection were varied; 63 patients had positive blood cultures with or without involvement of one or more organs or tissues, and 40 patients had negative blood cultures but one or more organs or tissues involved. All isolates were maintained at -70°C in TSB with 15% glycerol.
Detection of genomic islands
A single bacterial colony was inoculated into TSB and incubated overnight in air at 37°C, after which genomic DNA was extracted using the Wizard Genomic DNA purification kit (Promega, WI, USA). The presence or absence of five GIs (GI 2, prophage φK96243; GI 6, prophage-like; GI 9, prophage-like; GI 11, putative integrated plasmid; and GI 16, metabolic island) , was defined by PCR. Two target genes were selected for each island, as follows: BPSL0130 and BPSL0135 (GI 2), BPSL 1138 and BPSL1155 (GI 6), BPSL2578 and BPSL2579 (GI 9), BPSL3258 and BPSL3260 (GI 11), and BPSS2053 and BPSS2061 (GI 16). Targets encoded predominantly putative hypothetical proteins with unknown functions. Genes that encoded functions necessary for the transfer and maintenance of the island such as integrases were avoided as similar sequences may be found on other unrelated islands. Primers and cycling conditions are as described by Holden et al.  with the exception of the following: BPSL1030-f (5'-GCGCCGCTCGACTTCCTTCTCT), BPSL1030-r (5'-GAGGGGCCGGACTGCTACTTCAC); BPSL1138-f (5'-GATTTGGTTGGCGTCCGTGTTT), BPSL1138-r (5'-CGACCTTGGCCGAATTATGTGA); BPSS2061-f (5'-AACGCTCGCGCCCTTTAC), BPSS2061-r (5'-AATGCCCTTCCGAATCCTTTATG). A positive control (B. pseudomallei K96243) and a negative control (reaction mixture without DNA) were included in each PCR run. To verify the absence of the genomic islands, a second round of PCR was performed for all isolates to amplify the putative insertion site by using primers situated within genes immediately abutting each putative island, as previously described . The expected result was either no PCR product (interpreted as being consistent with the presence of the island) or a product whose size was consistent with amplification across the insertion site minus island PCR. Amplifications were performed using a PTC-0200 DNA engine (MJ Research, Cambridge, MA) with Taq polymerase (Promega), and aliquots of reaction mixtures were analyzed by 1% agarose gel electrophoresis. An island was classified as present if either or both target genes were positive, and negative if no amplification products were obtained for either target gene but a product was amplified across the putative insertion site.
Multiple organ involvement was defined when there was >1 non-contiguous focus of infection, not including blood. Pneumonia was defined as the presence of clinical features plus radiographic changes and/or sputum culture positive for B. pseudomallei. Impaired renal function was defined as an estimated glomerular filtration rate below 60 mL/min/1.73 m2 during admission. Glomerular filtration rate was estimated using an abbreviated form of the Modification of Diet in Renal Disease study equation . Impaired hepatic function was defined as an elevation of aminotransferase more than 5 times the upper normal limit, or the presence of jaundice during admission. Hypotension was defined as a systolic BP <90 mmHg on or during admission.
MLST has been performed by us for all of the Thai isolates used in this study , MLST profiles for which are maintained on the B. pseudomallei MLST web site . Genomic islands were previously identified in the genome of B. pseudomallei K96243 (accession numbers BX571965 and BX571966) . Comparison of genome sequences was performed with ACT (Artemis Comparison Tool)  using BLASTN and TBLASTX . Comparisons were made using the 3 other complete B. pseudomallei genome sequences 1106a (accession numbers CP000572 and CP000573), 1710b (accession numbers CP000124 and CP000125) and 668 (accession numbers CP000570 and CP000571), the 6 B. pseudomallei strains currently undergoing whole genome sequencing (406e, 1106b, 1710a, Pasteur 6068, S13, and 1655) by JCVI , and B. thailandensis E264 (accession numbers CP000086 and CP000085) .
Statistical tests were performed using the statistical program STATA/SE, version 9.0 (StataCorp LP, College Station, Tx.). Proportions were compared using the Chi-square test, or Fisher's exact test where appropriate. Comparison of continuous data was performed using the Student's t test. Unifactorial analysis was used to examine the association between presence of each GI, cumulative number of GIs and clinical factors.
We are grateful to staff at Sappasithiprasong Hospital, in particular to Professor Wipada Chaowagul, and to the Mahidol-Oxford Tropical Medicine Research Unit for their support. SJP was supported by a Wellcome Trust Career Development Award in Clinical Tropical Medicine. This study was funded by the Wellcome Trust.
- Currie BJ, Fisher DA, Howard DM, Burrow JNC, Lo D, Selva-Nayagam S, Anstey NM, Huffam SE, Snelling PL, Marks PJ, Stephens DP, Lum GD, Jacups SP, Krause VL: Endemic melioidosis in tropical northern Australia: A 10-year prospective study and review of the literature. Clinical Infectious Diseases. 2000, 31 (4): 981-986. 10.1086/318116.PubMedView ArticleGoogle Scholar
- White NJ: Melioidosis. Lancet. 2003, 361 (9370): 1715-1722. 10.1016/S0140-6736(03)13374-0.PubMedView ArticleGoogle Scholar
- Holden MT, Titball RW, Peacock SJ, Cerdeno-Tarraga AM, Atkins T, Crossman LC, Pitt T, Churcher C, Mungall K, Bentley SD, Sebaihia M, Thomson NR, Bason N, Beacham IR, Brooks K, Brown KA, Brown NF, Challis GL, Cherevach I, Chillingworth T, Cronin A, Crossett B, Davis P, DeShazer D, Feltwell T, Fraser A, Hance Z, Hauser H, Holroyd S, Jagels K, Keith KE, Maddison M, Moule S, Price C, Quail MA, Rabbinowitsch E, Rutherford K, Sanders M, Simmonds M, Songsivilai S, Stevens K, Tumapa S, Vesaratchavest M, Whitehead S, Yeats C, Barrell BG, Oyston PC, Parkhill J: Genomic plasticity of the causative agent of melioidosis, Burkholderia pseudomallei. Proc Natl Acad Sci U S A. 2004, 101 (39): 14240-14245. 10.1073/pnas.0403302101.PubMedPubMed CentralView ArticleGoogle Scholar
- Schmidt H, Hensel M: Pathogenicity islands in bacterial pathogenesis. Clinical Microbiology Reviews. 2004, 17 (1): 14-+. 10.1128/CMR.17.1.14-56.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Ngauy V, Lemeshev Y, Sadkowski L, Crawford G: Cutaneous melioidosis in a man who was taken as a prisoner of war by the Japanese during World War II. J Clin Microbiol. 2005, 43 (2): 970-972. 10.1128/JCM.43.2.970-972.2005.PubMedPubMed CentralView ArticleGoogle Scholar
- Yee KC, Lee MK, Chua CT, Puthucheary SD: Melioidosis, the Great Mimicker - a Report of 10 Cases from Malaysia. J Trop Med Hyg. 1988, 91 (5): 249-254.PubMedGoogle Scholar
- Currie BJ, Jacups SP, Cheng AC, Fisher DA, Anstey NM, Huffam SE, Krause VL: Melioidosis epidemiology and risk factors from a prospective whole-population study in northern Australia. Tropical Medicine & International Health. 2004, 9 (11): 1167-1174. 10.1111/j.1365-3156.2004.01328.x.View ArticleGoogle Scholar
- Godoy D, Randle G, Simpson AJ, Aanensen DM, Pitt TL, Kinoshita R, Spratt BG: Multilocus sequence typing and evolutionary relationships among the causative agents of melioidosis and glanders, Burkholderia pseudomallei and Burkholderia mallei. J Clin Microbiol. 2003, 41 (5): 2068-2079. 10.1128/JCM.41.5.2068-2079.2003.PubMedPubMed CentralView ArticleGoogle Scholar
- Turner KME, Feil EJ: The secret life of the multilocus sequence type. International Journal of Antimicrobial Agents. 2007, 29 (2): 129-135. 10.1016/j.ijantimicag.2006.11.002.PubMedView ArticleGoogle Scholar
- Cheng AC, Godoy D, Mayo M, Gal D, Spratt BG, Currie BJ: Isolates of Burkholderia pseudomallei from northern Australia are distinct by multilocus sequence typing, but strain types do not correlate with clinical presentation. Journal of Clinical Microbiology. 2004, 42 (12): 5477-5483. 10.1128/JCM.42.12.5477-5483.2004.PubMedPubMed CentralView ArticleGoogle Scholar
- Hanage WP, Fraser C, Spratt BG: The impact of homologous recombination on the generation of diversity in bacteria. Journal of Theoretical Biology. 2006, 239 (2): 210-219. 10.1016/j.jtbi.2005.08.035.PubMedView ArticleGoogle Scholar
- Feil EJ, Cooper JE, Grundmann H, Robinson DA, Enright MC, Berendt T, Peacock SJ, Smith JM, Murphy M, Spratt BG, Moore CE, Day NPJ: How clonal is Staphylococcus aureus?. J Bacteriol. 2003, 185 (11): 3307-3316. 10.1128/JB.185.11.3307-3316.2003.PubMedPubMed CentralView ArticleGoogle Scholar
- Holden MTG, Feil EJ, Lindsay JA, Peacock SJ, Day NPJ, Enright MC, Foster TJ, Moore CE, Hurst L, Atkin R, Barron A, Bason N, Bentley SD, Chillingworth C, Chillingworth T, Churcher C, Clark L, Corton C, Cronin A, Doggett J, Dowd L, Feltwell T, Hance Z, Harris B, Hauser H, Holroyd S, Jagels K, James KD, Lennard N, Line A, Mayes R, Moule S, Mungall K, Ormond D, Quail MA, Rabbinowitsch E, Rutherford K, Sanders M, Sharp S, Simmonds M, Stevens K, Whitehead S, Barrell BG, Spratt BG, Parkhill J: Complete genomes of two clinical Staphylococcus aureus strains: Evidence for the rapid evolution of virulence and drug resistance. Proc Natl Acad Sci U S A. 2004, 101 (26): 9786-9791. 10.1073/pnas.0402521101.PubMedPubMed CentralView ArticleGoogle Scholar
- Lindsay JA, Moore CE, Day NP, Peacock SJ, Witney AA, Stabler RA, Husain SE, Butcher PD, Hinds J: Microarrays reveal that each of the ten dominant lineages of Staphylococcus aureus has a unique combination of surface-associated and regulatory genes. J Bacteriol. 2006, 188 (2): 669-676. 10.1128/JB.188.2.669-676.2006.PubMedPubMed CentralView ArticleGoogle Scholar
- Peacock SJ, Moore CE, Justice A, Kantzanou M, Story L, Mackie K, O'Neil G, Day NPJ: Virulent combinations of adhesin and toxin genes in natural populations of Staphylococcus aureus. Infect Immun. 2002, 70 (9): 4987-4996. 10.1128/IAI.70.9.4987-4996.2002.PubMedPubMed CentralView ArticleGoogle Scholar
- Vesaratchavest M, Tumapa S, Day NPJ, Wuthiekanun V, Chierakul W, Holden MTG, White NJ, Currie BJ, Spratt BG, Feil EJ, Peacock SJ: Nonrandom distribution of Burkholderia pseudomallei clones in relation to geographical location and virulence. Journal of Clinical Microbiology. 2006, 44 (7): 2553-2557. 10.1128/JCM.00629-06.PubMedPubMed CentralView ArticleGoogle Scholar
- Duangsonk K, Gal D, Mayo M, Hart CA, Currie BJ, Winstanley C: Use of a variable amplicon typing scheme reveals considerable variation in the accessory genomes of isolates of Burkholderia pseudomallei. Journal of Clinical Microbiology. 2006, 44 (4): 1323-1334. 10.1128/JCM.44.4.1323-1334.2006.PubMedPubMed CentralView ArticleGoogle Scholar
- Smith MD, Wuthiekanun V, Walsh AL, White NJ: Quantitative Recovery of Burkholderia pseudomallei from Soil in Thailand. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1995, 89 (5): 488-490. 10.1016/0035-9203(95)90078-0.PubMedView ArticleGoogle Scholar
- K/DOQI clinical practice guidelines for chronic kidney disease: evaluation, classification, and stratification. American Journal of Kidney Disease. 2002, 39: S1-266.Google Scholar
- Burkholderia pseudomallei MLST web site. [http://www.mlst.net]
- Carver TJ, Rutherford K, Berriman M, Rajandream MA, Barrell B, Parkhill J: ACT: the Artemis comparison tool. Bioinformatics. 2005, 21: 3422-3423. 10.1093/bioinformatics/bti553.PubMedView ArticleGoogle Scholar
- Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ: Basic Local Alignment Search Tool. J Mol Biol. 1990, 215 (3): 403-410.PubMedView ArticleGoogle Scholar
- The Institute of Genome Research. [http://www.tigr.org]
- Yu Y, Kim HS, Chua HH, Lin CH, Sim SH, Lin D, Derr A, Engels R, DeShazer D, Birren B, Nierman WC, Tan P: Genomic patterns of pathogen evolution revealed by comparison of Burkholderia pseudomallei, the causative agent of melioidosis, to avirulent Burkholderia thailandensis. BMC Microbiol. 2006, 6: 46-10.1186/1471-2180-6-46.PubMedPubMed CentralView ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.