barCoder: a tool to generate unique, orthogonal genetic tags for qPCR detection.

BACKGROUND: Tracking dispersal of microbial populations in the environment requires specific detection methods that discriminate between the target strain and all potential natural and artificial interferents, including previously utilized tester strains. Recent work has shown that genomic insertion of short identification tags, called "barcodes" here, allows detection of chromosomally tagged strains by real-time PCR. Manual design of these barcodes is feasible for small sets, but expansion of the technique to larger pools of distinct and well-functioning assays would be significantly aided by software-guided design. RESULTS: Here we introduce barCoder, a bioinformatics tool that facilitates the process of creating sets of uniquely identifiable barcoded strains. barCoder utilizes the genomic sequence of the target strain and a set of user-specified PCR parameters to generate a list of suggested barcode "modules" that consist of binding sites for primers and probes, and appropriate spacer sequences. Each module is designed to yield optimal PCR amplification and unique identification. Optimal amplification includes metrics such as ideal melting temperature and G+C content, appropriate spacing, and minimal stem-loop formation; unique identification includes low BLAST hits against the target organism, previously generated barcode modules, and databases (such as NCBI). We tested the ability of our algorithm to suggest appropriate barcodes by generating 12 modules for Bacillus thuringiensis serovar kurstaki-a simulant for the potential biowarfare agent Bacillus anthracis-and three each for other potential target organisms with variable G+C content. Real-time PCR detection assays directed at barcodes were specific and yielded minimal cross-reactivity with a panel of near-neighbor and potential contaminant materials. CONCLUSIONS: The barCoder algorithm facilitates the generation of synthetically barcoded biological simulants by (a) eliminating the task of creating modules by hand, (b) minimizing optimization of PCR assays, and (c) reducing effort wasted on non-unique barcode modules.

Publisher Correction: Shared strategies for ß-lactam catabolism in the soil microbiome.

In the version of the article originally published, the x axis of the graph in Fig. 4d was incorrectly labeled as "Retention time (min)". It should read "Reaction time (min)". The 'deceased' footnote was also formatted incorrectly when published. The footnote text itself should include the name of co-author Tara A. Gianoulis in addition to the previous link to her name in the author list through footnote number 10. The errors have been corrected in the HTML and PDF versions of the article.

Shared strategies for ß-lactam catabolism in the soil microbiome.

The soil microbiome can produce, resist, or degrade antibiotics and even catabolize them. While resistance genes are widely distributed in the soil, there is a dearth of knowledge concerning antibiotic catabolism. Here we describe a pathway for penicillin catabolism in four isolates. Genomic and transcriptomic sequencing revealed ß-lactamase, amidase, and phenylacetic acid catabolon upregulation. Knocking out part of the phenylacetic acid catabolon or an apparent penicillin utilization operon (put) resulted in loss of penicillin catabolism in one isolate. A hydrolase from the put operon was found to degrade in vitro benzylpenicilloic acid, the ß-lactamase penicillin product. To test the generality of this strategy, an Escherichia coli strain was engineered to co-express a ß-lactamase and a penicillin amidase or the put operon, enabling it to grow using penicillin or benzylpenicilloic acid, respectively. Elucidation of additional pathways may allow bioremediation of antibiotic-contaminated soils and discovery of antibiotic-remodeling enzymes with industrial utility.

The Essential Role of Hypermutation in Rapid Adaptation to Antibiotic Stress.

A common outcome of antibiotic exposure in patients and in vitro is the evolution of a hypermutator phenotype that enables rapid adaptation by pathogens. While hypermutation is a robust mechanism for rapid adaptation, it requires trade-offs between the adaptive mutations and the more common "hitchhiker" mutations that accumulate from the increased mutation rate. Using quantitative experimental evolution, we examined the role of hypermutation in driving the adaptation of Pseudomonas aeruginosa to colistin. Metagenomic deep sequencing revealed 2,657 mutations at ≥5% frequency in 1,197 genes and 761 mutations in 29 endpoint isolates. By combining genomic information, phylogenetic analyses, and statistical tests, we showed that evolutionary trajectories leading to resistance could be reliably discerned. In addition to known alleles such as pmrB, hypermutation allowed identification of additional adaptive alleles with epistatic relationships. Although hypermutation provided a short-term fitness benefit, it was detrimental to overall fitness. Alarmingly, a small fraction of the colistin-adapted population remained colistin susceptible and escaped hypermutation. In a clinical population, such cells could play a role in reestablishing infection upon withdrawal of colistin. We present here a framework for evaluating the complex evolutionary trajectories of hypermutators that applies to both current and emerging pathogen populations.

A Standard Method To Inactivate Bacillus anthracis Spores to Sterility via Gamma Irradiation.

In 2015, a laboratory of the United States Department of Defense (DoD) inadvertently shipped preparations of gamma-irradiated spores of Bacillus anthracis that contained live spores. In response, a systematic evidence-based method for preparing, concentrating, irradiating, and verifying the inactivation of spore materials was developed. We demonstrate the consistency of spore preparations across multiple biological replicates and show that two different DoD institutions independently obtained comparable dose-inactivation curves for a monodisperse suspension of B. anthracis spores containing 3 × 1010 CFU. Spore preparations from three different institutions and three strain backgrounds yielded similar decimal reduction (D10) values and irradiation doses required to ensure sterility (DSAL) to the point at which the probability of detecting a viable spore is 10-6 Furthermore, spores of a genetically tagged strain of B. anthracis strain Sterne were used to show that high densities of dead spores suppress the recovery of viable spores. Together, we present an integrated method for preparing, irradiating, and verifying the inactivation of spores of B. anthracis for use as standard reagents for testing and evaluating detection and diagnostic devices and techniques.IMPORTANCE The inadvertent shipment by a U.S. Department of Defense (DoD) laboratory of live Bacillus anthracis (anthrax) spores to U.S. and international destinations revealed the need to standardize inactivation methods for materials derived from biological select agents and toxins (BSAT) and for the development of evidence-based methods to prevent the recurrence of such an event. Following a retrospective analysis of the procedures previously employed to generate inactivated B. anthracis spores, a study was commissioned by the DoD to provide data required to support the production of inactivated spores for the biodefense community. The results of this work are presented in this publication, which details the method by which spores can be prepared, irradiated, and tested, such that the chance of finding residual living spores in any given preparation is 1/1,000,000. These irradiated spores are used to test equipment and methods for the detection of agents of biological warfare and bioterrorism.

Mutations in global regulators lead to metabolic selection during adaptation to complex environments.

Adaptation to ecologically complex environments can provide insights into the evolutionary dynamics and functional constraints encountered by organisms during natural selection. Adaptation to a new environment with abundant and varied resources can be difficult to achieve by small incremental changes if many mutations are required to achieve even modest gains in fitness. Since changing complex environments are quite common in nature, we investigated how such an epistatic bottleneck can be avoided to allow rapid adaptation. We show that adaptive mutations arise repeatedly in independently evolved populations in the context of greatly increased genetic and phenotypic diversity. We go on to show that weak selection requiring substantial metabolic reprogramming can be readily achieved by mutations in the global response regulator arcA and the stress response regulator rpoS. We identified 46 unique single-nucleotide variants of arcA and 18 mutations in rpoS, nine of which resulted in stop codons or large deletions, suggesting that subtle modulations of ArcA function and knockouts of rpoS are largely responsible for the metabolic shifts leading to adaptation. These mutations allow a higher order metabolic selection that eliminates epistatic bottlenecks, which could occur when many changes would be required. Proteomic and carbohydrate analysis of adapting E. coli populations revealed an up-regulation of enzymes associated with the TCA cycle and amino acid metabolism, and an increase in the secretion of putrescine. The overall effect of adaptation across populations is to redirect and efficiently utilize uptake and catabolism of abundant amino acids. Concomitantly, there is a pronounced spread of more ecologically limited strains that results from specialization through metabolic erosion. Remarkably, the global regulators arcA and rpoS can provide a "one-step" mechanism of adaptation to a novel environment, which highlights the importance of global resource management as a powerful strategy to adaptation.

Whole-Genome Sequencing in Microbial Forensic Analysis of Gamma-Irradiated Microbial Materials.

Effective microbial forensic analysis of materials used in a potential biological attack requires robust methods of morphological and genetic characterization of the attack materials in order to enable the attribution of the materials to potential sources and to exclude other potential sources. The genetic homogeneity and potential intersample variability of many of the category A to C bioterrorism agents offer a particular challenge to the generation of attributive signatures, potentially requiring whole-genome or proteomic approaches to be utilized. Currently, irradiation of mail is standard practice at several government facilities judged to be at particularly high risk. Thus, initial forensic signatures would need to be recovered from inactivated (nonviable) material. In the study described in this report, we determined the effects of high-dose gamma irradiation on forensic markers of bacterial biothreat agent surrogate organisms with a particular emphasis on the suitability of genomic DNA (gDNA) recovered from such sources as a template for whole-genome analysis. While irradiation of spores and vegetative cells affected the retention of Gram and spore stains and sheared gDNA into small fragments, we found that irradiated material could be utilized to generate accurate whole-genome sequence data on the Illumina and Roche 454 sequencing platforms.

Pathosphere.org: pathogen detection and characterization through a web-based, open source informatics platform.

BACKGROUND: The detection of pathogens in complex sample backgrounds has been revolutionized by wide access to next-generation sequencing (NGS) platforms. However, analytical methods to support NGS platforms are not as uniformly available. Pathosphere (found at Pathosphere.org) is a cloud - based open - sourced community tool that allows for communication, collaboration and sharing of NGS analytical tools and data amongst scientists working in academia, industry and government. The architecture allows for users to upload data and run available bioinformatics pipelines without the need for onsite processing hardware or technical support. RESULTS: The pathogen detection capabilities hosted on Pathosphere were tested by analyzing pathogen-containing samples sequenced by NGS with both spiked human samples as well as human and zoonotic host backgrounds. Pathosphere analytical pipelines developed by Edgewood Chemical Biological Center (ECBC) identified spiked pathogens within a common sample analyzed by 454, Ion Torrent, and Illumina sequencing platforms. ECBC pipelines also correctly identified pathogens in human samples containing arenavirus in addition to animal samples containing flavivirus and coronavirus. These analytical methods were limited in the detection of sequences with limited homology to previous annotations within NCBI databases, such as parvovirus. Utilizing the pipeline-hosting adaptability of Pathosphere, the analytical suite was supplemented by analytical pipelines designed by the United States Army Medical Research Insititute of Infectious Diseases and Walter Reed Army Institute of Research (USAMRIID-WRAIR). These pipelines were implemented and detected parvovirus sequence in the sample that the ECBC iterative analysis previously failed to identify. CONCLUSIONS: By accurately detecting pathogens in a variety of samples, this work demonstrates the utility of Pathosphere and provides a platform for utilizing, modifying and creating pipelines for a variety of NGS technologies developed to detect pathogens in complex sample backgrounds. These results serve as an exhibition for the existing pipelines and web-based interface of Pathosphere as well as the plug-in adaptability that allows for integration of newer NGS analytical software as it becomes available.

Protein export by the mycobacterial SecA2 system is determined by the preprotein mature domain.

At the core of the bacterial general secretion (Sec) pathway is the SecA ATPase, which powers translocation of unfolded preproteins containing Sec signal sequences through the SecYEG membrane channel. Mycobacteria have two nonredundant SecA homologs: SecA1 and SecA2. While the essential SecA1 handles "housekeeping" export, the nonessential SecA2 exports a subset of proteins and is required for Mycobacterium tuberculosis virulence. Currently, it is not understood how SecA2 contributes to Sec export in mycobacteria. In this study, we focused on identifying the features of two SecA2 substrates that target them to SecA2 for export, the Ms1704 and Ms1712 lipoproteins of the model organism Mycobacterium smegmatis. We found that the mature domains of Ms1704 and Ms1712, not the N-terminal signal sequences, confer SecA2-dependent export. We also demonstrated that the lipid modification and the extreme N terminus of the mature protein do not impart the requirement for SecA2 in export. We further showed that the Ms1704 mature domain can be efficiently exported by the twin-arginine translocation (Tat) pathway. Because the Tat system exports only folded proteins, this result implies that SecA2 substrates can fold in the cytoplasm and suggests a putative role of SecA2 in enabling export of such proteins. Thus, the mycobacterial SecA2 system may represent another way that bacteria solve the problem of exporting proteins that can fold in the cytoplasm.

Effect of the Bacillus atrophaeus subsp. globigii Spo0F H101R mutation on strain fitness.

Sporulation is a critical developmental process in Bacillus spp. that, once initiated, removes the possibility of further growth until germination. Therefore, the threshold conditions triggering sporulation are likely to be subject to evolutionary constraint. Our previous studies revealed two spontaneous hypersporulating mutants of Bacillus atrophaeus subsp. globigii, both containing point mutations in the spo0F gene. One of these strains (Detrick-2; contains the spo0F101 allele with a C:T [His101Arg] substitution) had been deliberately selected in the early 1940s as an anthrax surrogate. To determine whether the experimental conditions used during the selection of the "military" strains could have supported the emergence of hypersporulating variants, the relative fitness of strain Detrick-2 was measured in several experimental settings modeled on experimental conditions employed during its development in the 1940s as a simulant. The congenic strain Detrick-1 contained a wild-type spo0F gene and sporulated like the wild-type strain. The relative fitness of Detrick-1 and Detrick-2 was evaluated in competition experiments using quantitative single nucleotide polymorphism (SNP)-specific real-time PCR assays directed at the C:T substitution. The ancestral strain Detrick-1 had a fitness advantage under all conditions tested except when competing cultures were subjected to frequent heat shocks. The hypersporulating strain gained the maximum fitness advantage when cultures were grown at low oxygen tension and when heat shock was applied soon after the formation of the first heat-resistant spores. This is interpreted as gain of fitness by the hypersporulating strain in fast-changing fluctuating environments as a result of the increased rate of switching to the sporulating phenotype.

Genetic barcodes for improved environmental tracking of an anthrax simulant.

The development of realistic risk models that predict the dissemination, dispersion and persistence of potential biothreat agents have utilized nonpathogenic surrogate organisms such as Bacillus atrophaeus subsp. globigii or commercial products such as Bacillus thuringiensis subsp. kurstaki. Comparison of results from outdoor tests under different conditions requires the use of genetically identical strains; however, the requirement for isogenic strains limits the ability to compare other desirable properties, such as the behavior in the environment of the same strain prepared using different methods. Finally, current methods do not allow long-term studies of persistence or reaerosolization in test sites where simulants are heavily used or in areas where B. thuringiensis subsp. kurstaki is applied as a biopesticide. To create a set of genetically heterogeneous yet phenotypically indistinguishable strains so that variables intrinsic to simulations (e.g., sample preparation) can be varied and the strains can be tested under otherwise identical conditions, we have developed a strategy of introducing small genetic signatures ("barcodes") into neutral regions of the genome. The barcodes are stable over 300 generations and do not impact in vitro growth or sporulation. Each barcode contains common and specific tags that allow differentiation of marked strains from wild-type strains and from each other. Each tag is paired with specific real-time PCR assays that facilitate discrimination of barcoded strains from wild-type strains and from each other. These uniquely barcoded strains will be valuable tools for research into the environmental fate of released organisms by providing specific artificial detection signatures.

Detection and tracking of a novel genetically tagged biological simulant in the environment.

A variant of Bacillus thuringiensis subsp. kurstaki containing a single, stable copy of a uniquely amplifiable DNA oligomer integrated into the genome for tracking the fate of biological agents in the environment was developed. The use of genetically tagged spores overcomes the ambiguity of discerning the test material from pre-existing environmental microflora or from previously released background material. In this study, we demonstrate the utility of the genetically "barcoded" simulant in a controlled indoor setting and in an outdoor release. In an ambient breeze tunnel test, spores deposited on tiles were reaerosolized and detected by real-time PCR at distances of 30 m from the point of deposition. Real-time PCR signals were inversely correlated with distance from the seeded tiles. An outdoor release of powdered spore simulant at Aberdeen Proving Ground, Edgewood, MD, was monitored from a distance by a light detection and ranging (LIDAR) laser. Over a 2-week period, an array of air sampling units collected samples were analyzed for the presence of viable spores and using barcode-specific real-time PCR assays. Barcoded B. thuringiensis subsp. kurstaki spores were unambiguously identified on the day of the release, and viable material was recovered in a pattern consistent with the cloud track predicted by prevailing winds and by data tracks provided by the LIDAR system. Finally, the real-time PCR assays successfully differentiated barcoded B. thuringiensis subsp. kurstaki spores from wild-type spores under field conditions.

Detection of Burkholderia pseudomallei O-antigen serotypes in near-neighbor species.

BACKGROUND: Burkholderia pseudomallei is the etiological agent of melioidosis and a CDC category B select agent with no available effective vaccine. Previous immunizations in mice have utilized the lipopolysaccharide (LPS) as a potential vaccine target because it is known as one of the most important antigenic epitopes in B. pseudomallei. Complicating this strategy are the four different B. pseudomallei LPS O-antigen types: A, B, B2, and rough. Sero-crossreactivity is common among O-antigens of Burkholderia species. Here, we identified the presence of multiple B. pseudomallei O-antigen types and sero-crossreactivity in its near-neighbor species. RESULTS: PCR screening of O-antigen biosynthesis genes, phenotypic characterization using SDS-PAGE, and immunoblot analysis showed that majority of B. mallei and B. thailandensis strains contained the typical O-antigen type A. In contrast, most of B. ubonensis and B. thailandensis-like strains expressed the atypical O-antigen types B and B2, respectively. Most B. oklahomensis strains expressed a distinct and non-seroreactive O-antigen type, except strain E0147 which expressed O-antigen type A. O-antigen type B2 was also detected in B. thailandensis 82172, B. ubonensis MSMB108, and Burkholderia sp. MSMB175. Interestingly, B. thailandensis-like MSMB43 contained a novel serotype B positive O-antigen. CONCLUSIONS: This study expands the number of species which express B. pseudomallei O-antigen types. Further work is required to elucidate the full structures and how closely these are to the B. pseudomallei O-antigens, which will ultimately determine the efficacy of the near-neighbor B serotypes for vaccine development.

Putative Phenotypically Neutral Genomic Insertion Points in Prokaryotes.

The barriers to effective genome editing in diverse prokaryotic organisms have been falling at an accelerated rate. As editing becomes easier in more organisms, quickly identifying genomic locations to insert new genetic functions without disrupting organism fitness becomes increasingly useful. When the insertion is noncoding DNA for applications such as information storage or barcoding, a neutral insertion point can be especially important. Here we describe an approach to identify putatively neutral insertion sites in prokaryotes. An algorithm (targetFinder) finds convergently transcribed genes with gap sizes within a specified range, and looks for annotations within the gaps. We report putative editing targets for 10 common synthetic biology chassis organisms, including coverage of available RNA-seq data, and provide software to apply to others. We further experimentally evaluate the neutrality of six identified targets in Escherichia coli through insertion of a DNA barcode. We anticipate this information and the accompanying tool will prove useful for synthetic biologists seeking neutral insertion points for genome editing.

The ChpR transcriptional regulator of Sinorhizobium meliloti senses 3,5,6-trichloropyridinol, a degradation product of the organophosphate pesticide chlorpyrifos.

The global use of organophosphate insecticides (OPPs) and the growing concern of off-target side effects due to OPP exposure has prompted the need for sensitive and economical detection methods. Here we set out to engineer a previously identified OPP responsive transcription factor, ChpR, from Sinorhizobium melilotii to respond to alternative OPPs and generate a repertoire of whole-cell biosensors for OPPs. The ChpR transcription factor and cognate promoter P chpA, have been shown to activate transcription in the presence of the OPP chlorpyrifos (CPF). Utilizing a GFP reporter regulated by ChpR in a whole-cell biosensor we found that the system responds significantly better to 3,5,6-trichloro-2-pyridinol (TCP), the main degradation product of CPF, compared to CPF itself. This biosensor was able to respond to TCP at 390 nM within 4 h compared to 50 µM of CPF in 7 h. The ChpR-P chpA , and the activating ligand TCP, were able to regulate expression of a kanamycin resistance/sucrose sensitivity (kan/sacB) selection/counterselection module suitable for high throughput mutagenesis screening studies. The ability to control both GFP and the kan/sacB module demonstrates the utility of this reporter for the detection of CPF affected areas. The ChpR-P chpA system serves as an additional positive regulator switch to add to the growing repertoire of controllers available within synthetic biology.

ATPase activity of Mycobacterium tuberculosis SecA1 and SecA2 proteins and its importance for SecA2 function in macrophages.

The Sec-dependent translocation pathway that involves the essential SecA protein and the membrane-bound SecYEG translocon is used to export many proteins across the cytoplasmic membrane. Recently, several pathogenic bacteria, including Mycobacterium tuberculosis, were shown to possess two SecA homologs, SecA1 and SecA2. SecA1 is essential for general protein export. SecA2 is specific for a subset of exported proteins and is important for M. tuberculosis virulence. The enzymatic activities of two SecA proteins from the same microorganism have not been defined for any bacteria. Here, M. tuberculosis SecA1 and SecA2 are shown to bind ATP with high affinity, though the affinity of SecA1 for ATP is weaker than that of SecA2 or Escherichia coli SecA. Amino acid substitution of arginine or alanine for the conserved lysine in the Walker A motif of SecA2 eliminated ATP binding. We used the SecA2(K115R) variant to show that ATP binding was necessary for the SecA2 function of promoting intracellular growth of M. tuberculosis in macrophages. These results are the first to show the importance of ATPase activity in the function of accessory SecA2 proteins.

Tracking a serial killer: Integrating phylogenetic relationships, epidemiology, and geography for two invasive meningococcal disease outbreaks.

BACKGROUND: While overall rates of meningococcal disease have been declining in the United States for the past several decades, New York City (NYC) has experienced two serogroup C meningococcal disease outbreaks in 2005-2006 and in 2010-2013. The outbreaks were centered within drug use and sexual networks, were difficult to control, and required vaccine campaigns. METHODS: Whole Genome Sequencing (WGS) was used to analyze preserved meningococcal isolates collected before and during the two outbreaks. We integrated and analyzed epidemiologic, geographic, and genomic data to better understand transmission networks among patients. Betweenness centrality was used as a metric to understand the most important geographic nodes in the transmission networks. Comparative genomics was used to identify genes associated with the outbreaks. RESULTS: Neisseria meningitidis serogroup C (ST11/ET-37) was responsible for both outbreaks with each outbreak having distinct phylogenetic clusters. WGS did identify some misclassifications of isolates that were more distant from the outbreak strains, as well as those that should have been included based on high genomic similarity. Genomes for the second outbreak were more similar than the first and no polymorphism was found to either be unique or specific to either outbreak lineage. Betweenness centrality as applied to transmission networks based on phylogenetic analysis demonstrated that the outbreaks were transmitted within focal communities in NYC with few transmission events to other locations. CONCLUSIONS: Neisseria meningitidis is an ever changing pathogen and comparative genomic analyses can help elucidate how it spreads geographically to facilitate targeted interventions to interrupt transmission.

Draft Genome Sequences of Three ß-Lactam-Catabolizing Soil Proteobacteria.

Most antibiotics are derived from the soil, but their catabolism there, which is necessary to close the antibiotic carbon cycle, remains uncharacterized. We report the first draft genome sequences of soil Proteobacteria identified for subsisting solely on ß-lactams as their carbon sources. The genomes encode multiple ß-lactamases, although their antibiotic catabolic pathways remain enigmatic.

Genomic sequences of six botulinum neurotoxin-producing strains representing three clostridial species illustrate the mobility and diversity of botulinum neurotoxin genes.

The whole genomes for six botulinum neurotoxin-producing clostridial strains were sequenced to provide references for under-represented toxin types, bivalent strains or unusual toxin complexes associated with a bont gene. The strains include three Clostridium botulinum Group I strains (CDC 297, CDC 1436, and Prevot 594), a Group II C. botulinum strain (Eklund 202F), a Group IV Clostridium argentinense strain (CDC 2741), and a Group V Clostridium baratii strain (Sullivan). Comparisons of the Group I genomic sequences revealed close relationships and conservation of toxin gene locations with previously published Group I C. botulinum genomes. The bont/F6 gene of strain Eklund 202F was determined to be a chimeric toxin gene composed of bont/F1 and bont/F2. The serotype G strain CDC 2741 remained unfinished in 20 contigs with the bont/G located within a 1.15Mb contig, indicating a possible chromosomal location for this toxin gene. Within the genome of C. baratii Sullivan strain, direct repeats of IS1182 insertion sequence (IS) elements were identified flanking the bont/F7 toxin complex that may be the mechanism of bont insertion into C. baratii. Highlights of the six strains are described and release of their genomic sequences will allow further study of unusual neurotoxin-producing clostridial strains.

Pangenome Analysis of Burkholderia pseudomallei: Genome Evolution Preserves Gene Order despite High Recombination Rates.

The pangenomic diversity in Burkholderia pseudomallei is high, with approximately 5.8% of the genome consisting of genomic islands. Genomic islands are known hotspots for recombination driven primarily by site-specific recombination associated with tRNAs. However, recombination rates in other portions of the genome are also high, a feature we expected to disrupt gene order. We analyzed the pangenome of 37 isolates of B. pseudomallei and demonstrate that the pangenome is 'open', with approximately 136 new genes identified with each new genome sequenced, and that the global core genome consists of 4568±16 homologs. Genes associated with metabolism were statistically overrepresented in the core genome, and genes associated with mobile elements, disease, and motility were primarily associated with accessory portions of the pangenome. The frequency distribution of genes present in between 1 and 37 of the genomes analyzed matches well with a model of genome evolution in which 96% of the genome has very low recombination rates but 4% of the genome recombines readily. Using homologous genes among pairs of genomes, we found that gene order was highly conserved among strains, despite the high recombination rates previously observed. High rates of gene transfer and recombination are incompatible with retaining gene order unless these processes are either highly localized to specific sites within the genome, or are characterized by symmetrical gene gain and loss. Our results demonstrate that both processes occur: localized recombination introduces many new genes at relatively few sites, and recombination throughout the genome generates the novel multi-locus sequence types previously observed while preserving gene order.