Precise mapping of quantitative trait loci for resistance to southern leaf blight, caused by Cochliobolus heterostrophus race O, and flowering time using advanced intercross maize lines

The intermated B73 × Mo17 (IBM) population, an advanced intercross recombinant inbred line population derived from a cross between the maize lines B73 (susceptible) and Mo17 (resistant), was evaluated in four environments for resistance to southern leaf blight (SLB) disease caused by Cochliobolus heterostrophus race O. Two environments were artificially inoculated, while two were not inoculated and consequently had substantially lower disease pressure. Four common SLB resistance quantitative trait loci (QTL) were identified in all environments, two in bin 3.04 and one each in bins 1.10 and 8.02/3. There was no significant correlation between disease resistance and days to anthesis. A direct comparison was made between SLB QTL detected in two populations, independently derived from the same parental cross: the IBM advanced intercross population and a conventional recombinant inbred line population. Several QTL for SLB resistance were detected in both populations, with the IBM providing between 5 and, in one case, 50 times greater mapping resolution.     

Global insight into the distribution of velvet-like B protein in <Emphasis Type="Italic">Cochliobolus</Emphasis> species and implication in pathogenicity and fungicide resistance

The Cochliobolus genus consist of over 55 species among which the 5 most devastating are Cochliobolus carbonum, Cochliobolus heterostrophus, Cochliobolus miyabeanus, Crocus sativus and Cochliobolus lunatus causing damages in sorghum, wheat, rice, maize, cassava and soybean estimated at over 10 billion USD per annum worldwide. The dynamic pathogenicity of Cochliobolus species and the plethora of infected hosts is determined by the evolution of virulence determinants such as the velvet-like B protein (VelB). Nonetheless, the knowledge on the distribution of Cochliobolus VelB and its implication in pathogenicity and fungicide resistance are often lacking. By scanning through the annotated genomes of C. lunatus, C. heterostrophus, C. carbonum, C. victoriae, C. sativus and C. miyabeanus, it is revealed that the numbers of ortholog VelB and cognates vary. By using the phylogenetic approach, it is established that the diversification rates among velvet-domain-containing proteins for phytopathogenic Cochliobolus species could impact differently on their oxidant and fungicide resistance potentials, ability to form appressoria-like structures and infection pegs during infection. This study provides new insights into the pathogenicity evolution of Cochliobolus species at the VelB locus which is relevant for designing effective strategies for durable management of Cochliobolus diseases.     

Detection of Cochliobolus heterostrophus races in South China

Cochliobolus heterostrophus is the causal pathogen of the southern corn leaf blight (SCLB). There are three known races: race O, race C and race T. To determine which C. heterostrophus races comprise the field population in southern China and to assess diversity of these strains in terms of virulence, 200 isolates from diseased plants were collected in nine provinces/municipalities. All were race O, that is, no race T or race C isolates were found. Sixty race O isolates that sporulated well were chosen for further analysis. Virulence was measured using the integral optical density (IOD) of leaf lesions on four maize inbred lines. UPGMA cluster analysis of AFLP markers was applied to the 60 race O isolates plus control race O, T and C strains. Phylogenetic distribution, geographic location and virulence were not correlated. These results can provide valuable information for guidance in early warning and disease control.     

Locally Confined Clonal Complexes of Mycobacterium ulcerans in Two Buruli Ulcer Endemic Regions of Cameroon

Background Mycobacterium ulcerans is the causative agent of the necrotizing skin disease Buruli ulcer (BU), which has been reported from over 30 countries worldwide. The majority of notified patients come from West African countries, such as Côte d’Ivoire, Ghana, Benin and Cameroon. All clinical isolates of M. ulcerans from these countries are closely related and their genomes differ only in a limited number of single nucleotide polymorphisms (SNPs).ConclusionsOur comparative genomic analysis revealed that M. ulcerans clones diversify locally by the accumulation of SNPs. Case isolates coming from more recently emerging BU endemic areas, such as the Mapé river basin, may be less diverse than populations from longer standing disease foci, such as the Nyong river basin. Exchange of strains between distinct endemic areas seems to be rare and local clonal complexes can be easily distinguished by whole genome sequencing.     

Conservation of XYN11A and XYN11B xylanase genes in Bipolaris sorghicola,Cochliobolus sativus,Cochliobolus heterostrophus,and Cochliobolus spicifer

Two types of xylanase gene, XYN11A (XYL1) and XYN11B (XYL2), were amplified by PCR and partially sequenced in four phytopathogenic species of the ascomycete fungal genus Cochliobolus (anamorph genus Bipolaris). Three of the species, C. heterostrophus (B. maydis), C. sativus (B. sorokiniana), and Bipolaris sorghicola (no teleomorph known), are interrelated; the fourth, C. spicifer (B. spicifera), was found, through analysis of the 5.8S RNA and internal transcribed spacer (ITS) sequences of its ribosomal DNA, to be more distantly related to the other three. Isolates from all four species contain orthologous XYN11A and XYN11B genes, but a set of laboratory strains of C. heterostrophus gave no product corresponding to the XYN11B gene. The patterns of evolution of the two xylanase genes and ribosomal DNA sequences are mutually consistent; the results indicate that the two genes were present in the common ancestor of all Cochliobolus species and are evolving independently of each other. Received: 12 July 2001 / Accepted: 6 March 2002     

Characterization and potential evolutionary impact of transposable elements in the genome of Cochliobolus heterostrophus

Background

Cochliobolus heterostrophus is a dothideomycete that causes Southern Corn Leaf Blight disease. There are two races, race O and race T that differ by the absence (race O) and presence (race T) of ~ 1.2-Mb of DNA encoding genes responsible for the production of T-toxin, which makes race T much more virulent than race O. The presence of repetitive elements in fungal genomes is considered to be an important source of genetic variability between different species.

Results

A detailed analysis of class I and II TEs identified in the near complete genome sequence of race O was performed. In total in race O, 12 new families of transposons were identified. In silico evidence of recent activity was found for many of the transposons and analyses of expressed sequence tags (ESTs) demonstrated that these elements were actively transcribed. Various potentially active TEs were found near coding regions and may modify the expression and structure of these genes by acting as ectopic recombination sites. Transposons were found on scaffolds carrying polyketide synthase encoding genes, responsible for production of T-toxin in race T. Strong evidence of ectopic recombination was found, demonstrating that TEs can play an important role in the modulation of genome architecture of this species. The Repeat Induced Point mutation (RIP) silencing mechanism was shown to have high specificity in C. heterostrophus, acting only on transposons near coding regions.

Conclusions

New families of transposons were identified. In C. heterostrophus, the RIP silencing mechanism is efficient and selective. The co-localization of effector genes and TEs, therefore, exposes those genes to high rates of point mutations. This may accelerate the rate of evolution of these genes, providing a potential advantage for the host. Additionally, it was shown that ectopic recombination promoted by TEs appears to be the major event in the genome reorganization of this species and that a large number of elements are still potentially active. So, this study provides information about the potential impact of TEs on the evolution of C. heterostrophus.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-536) contains supplementary material, which is available to authorized users.     

Identification and Analysis of Genome-Wide SNPs Provide Insight into Signatures of Selection and Domestication in Channel Catfish (Ictalurus punctatus)

Domestication and selection for important performance traits can impact the genome, which is most often reflected by reduced heterozygosity in and surrounding genes related to traits affected by selection. In this study, analysis of the genomic impact caused by domestication and artificial selection was conducted by investigating the signatures of selection using single nucleotide polymorphisms (SNPs) in channel catfish (Ictalurus punctatus). A total of 8.4 million candidate SNPs were identified by using next generation sequencing. On average, the channel catfish genome harbors one SNP per 116 bp. Approximately 6.6 million, 5.3 million, 4.9 million, 7.1 million and 6.7 million SNPs were detected in the Marion, Thompson, USDA103, Hatchery strain, and wild population, respectively. The allele frequencies of 407,861 SNPs differed significantly between the domestic and wild populations. With these SNPs, 23 genomic regions with putative selective sweeps were identified that included 11 genes. Although the function for the majority of the genes remain unknown in catfish, several genes with known function related to aquaculture performance traits were included in the regions with selective sweeps. These included hypoxia-inducible factor 1β· HIFιβ ¨ and the transporter gene ATP-binding cassette sub-family B member 5 (ABCB5). HIF1β· is important for response to hypoxia and tolerance to low oxygen levels is a critical aquaculture trait. The large numbers of SNPs identified from this study are valuable for the development of high-density SNP arrays for genetic and genomic studies of performance traits in catfish.     

Complete Genomic Sequence and Comparative Analysis of the Genome Segments of Sweet Potato Chlorotic Stunt Virus in China

Background

Sweet potato chlorotic stunt virus (family Closteroviridae, genus Crinivirus) features a large bipartite, single-stranded, positive-sense RNA genome. To date, only three complete genomic sequences of SPCSV can be accessed through GenBank. SPCSV was first detected from China in 2011, only partial genomic sequences have been determined in the country. No report on the complete genomic sequence and genome structure of Chinese SPCSV isolates or the genetic relation between isolates from China and other countries is available.

Methodology/Principal Findings

The complete genomic sequences of five isolates from different areas in China were characterized. This study is the first to report the complete genome sequences of SPCSV from whitefly vectors. Genome structure analysis showed that isolates of WA and EA strains from China have the same coding protein as isolates Can181-9 and m2-47, respectively. Twenty cp genes and four RNA1 partial segments were sequenced and analyzed, and the nucleotide identities of complete genomic, cp, and RNA1 partial sequences were determined. Results indicated high conservation among strains and significant differences between WA and EA strains. Genetic analysis demonstrated that, except for isolates from Guangdong Province, SPCSVs from other areas belong to the WA strain. Genome organization analysis showed that the isolates in this study lack the p22 gene.

Conclusions/Significance

We presented the complete genome sequences of SPCSV in China. Comparison of nucleotide identities and genome structures between these isolates and previously reported isolates showed slight differences. The nucleotide identities of different SPCSV isolates showed high conservation among strains and significant differences between strains. All nine isolates in this study lacked p22 gene. WA strains were more extensively distributed than EA strains in China. These data provide important insights into the molecular variation and genomic structure of SPCSV in China as well as genetic relationships among isolates from China and other countries.     

Novel Conserved Genotypes Correspond to Antibiotic Resistance Phenotypes of E. coli Clinical Isolates

Current efforts to understand antibiotic resistance on the whole genome scale tend to focus on known genes even as high throughput sequencing strategies uncover novel mechanisms. To identify genomic variations associated with antibiotic resistance, we employed a modified genome-wide association study; we sequenced genomic DNA from pools of E. coli clinical isolates with similar antibiotic resistance phenotypes using SOLiD technology to uncover single nucleotide polymorphisms (SNPs) unanimously conserved in each pool. The multidrug-resistant pools were genotypically similar to SMS-3-5, a previously sequenced multidrug-resistant isolate from a polluted environment. The similarity was evenly spread across the entire genome and not limited to plasmid or pathogenicity island loci. Among the pools of clinical isolates, genomic variation was concentrated adjacent to previously reported inversion and duplication differences between the SMS-3-5 isolate and the drug-susceptible laboratory strain, DH10B. SNPs that result in non-synonymous changes in gyrA (encoding the well-known S83L allele associated with fluoroquinolone resistance), mutM, ligB, and recG were unanimously conserved in every fluoroquinolone-resistant pool. Alleles of the latter three genes are tightly linked among most sequenced E. coli genomes, and had not been implicated in antibiotic resistance previously. The changes in these genes map to amino acid positions in alpha helices that are involved in DNA binding. Plasmid-encoded complementation of null strains with either allelic variant of mutM or ligB resulted in variable responses to ultraviolet light or hydrogen peroxide treatment as markers of induced DNA damage, indicating their importance in DNA metabolism and revealing a potential mechanism for fluoroquinolone resistance. Our approach uncovered evidence that additional DNA binding enzymes may contribute to fluoroquinolone resistance and further implicate environmental bacteria as a reservoir for antibiotic resistance.     

Single mating type-specific genes and their 3′ UTRs control mating and fertility in Cochliobolus heterostrophus

To determine the number of proteins required for mating type (MAT) locus-regulated control of mating in Cochliobolus heterostrophus, MAT fragments of various sizes were expressed in MAT deletion strains. As little as 1.5?kb of MAT sequence, encoding a single unique protein in each mating type (MAT-1 and MAT-2), conferred mating ability, although an additional 160?bp of 3^′ UTR was needed for production of ascospores. No other mating type-specific genes involved in mating identity or fertility were found. Thus, although homologs of the C. heterostrophus MAT-1 and MAT-2 genes exist in the filamentous ascomycetes Neurospora crassa and Podospora anserina, C. heterostrophus does not appear to have mating type-specific homologs of two additional genes required by both N. crassa and P.?anserina for successful sexual reproduction. Three genes were identified in the common DNA flanking the MAT locus: a gene encoding a GTPase-activating protein and an ORF of unknown function lie 5^′ while a β-glucosidase encoding gene lies found 3^′. None of these genes appears to be involving in the mating process.     

Genome-Wide Genetic Diversity and Differentially Selected Regions among Suffolk,Rambouillet, Columbia,Polypay, and Targhee Sheep

Sheep are among the major economically important livestock species worldwide because the animals produce milk, wool, skin, and meat. In the present study, the Illumina OvineSNP50 BeadChip was used to investigate genetic diversity and genome selection among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds from the United States. After quality-control filtering of SNPs (single nucleotide polymorphisms), we used 48,026 SNPs, including 46,850 SNPs on autosomes that were in Hardy-Weinberg equilibrium and 1,176 SNPs on chromosome × for analysis. Phylogenetic analysis based on all 46,850 SNPs clearly separated Suffolk from Rambouillet, Columbia, Polypay, and Targhee, which was not surprising as Rambouillet contributed to the synthesis of the later three breeds. Based on pair-wise estimates of F _ST, significant genetic differentiation appeared between Suffolk and Rambouillet (F _ST = 0.1621), while Rambouillet and Targhee had the closest relationship (F _ST = 0.0681). A scan of the genome revealed 45 and 41 differentially selected regions (DSRs) between Suffolk and Rambouillet and among Rambouillet-related breed populations, respectively. Our data indicated that regions 13 and 24 between Suffolk and Rambouillet might be good candidates for evaluating breed differences. Furthermore, ovine genome v3.1 assembly was used as reference to link functionally known homologous genes to economically important traits covered by these differentially selected regions. In brief, our present study provides a comprehensive genome-wide view on within- and between-breed genetic differentiation, biodiversity, and evolution among Suffolk, Rambouillet, Columbia, Polypay, and Targhee sheep breeds. These results may provide new guidance for the synthesis of new breeds with different breeding objectives.     

Comparative Genomic Analysis of Escherichia coli O157:H7 Isolated from Super-Shedder and Low-Shedder Cattle

Cattle are the primary reservoir of the foodborne pathogen Escherichia coli O157:H7, with the concentration and frequency of E. coli O157:H7 shedding varying substantially among individual hosts. The term ‘‘super-shedder” has been applied to cattle that shed ≥10⁴ cfu E. coli O157:H7/g of feces. Super-shedders have been reported to be responsible for the majority of E. coli O157:H7 shed into the environment. The objective of this study was to determine if there are phenotypic and/or genotypic differences between E. coli O157:H7 isolates obtained from super-shedder compared to low-shedder cattle. From a total of 784 isolates, four were selected from low-shedder steers and six isolates from super-shedder steers (4.01–8.45 log cfu/g feces) for whole genome sequencing. Isolates were phage and clade typed, screened for substrate utilization, pH sensitivity, virulence gene profiles and Stx bacteriophage insertion (SBI) sites. A range of 89–2473 total single nucleotide polymorphisms (SNPs) were identified when sequenced strains were compared to E. coli O157:H7 strain Sakai. More non-synonymous SNP mutations were observed in low-shedder isolates. Pan-genomic and SNPs comparisons did not identify genetic segregation between super-shedder or low-shedder isolates. All super-shedder isolates and 3 of 4 of low-shedder isolates were typed as phage type 14a, SBI cluster 3 and SNP clade 2. Super-shedder isolates displayed increased utilization of galactitol, thymidine and 3-O-β-D-galactopyranosyl-D-arabinose when compared to low-shedder isolates, but no differences in SNPs were observed in genes encoding for proteins involved in the metabolism of these substrates. While genetic traits specific to super-shedder isolates were not identified in this study, differences in the level of gene expression or genes of unknown function may still contribute to some strains of E. coli O157:H7 reaching high densities within bovine feces.     

Genomic Characteristics of Chinese Borrelia burgdorferi Isolates

In China, B. burgdorferi, B.garinii, B. afzelii and B. yangtze sp. nov have been reported; B.garinii and B. afzelii are the main pathogenic genotypes. But until now only one Chinese strain was reported with whole genome sequence. In order to further understand the genomic characteristics and diversity of Chinese Borrelia strains, 5 isolates from China were sequenced and compared with the whole genome sequences of strains in other areas. The results showed a high degree of conservation within the linear chromosome of Chinese strains, whereas plasmid showed a much larger diversity according to the majority genomic information of plasmids. The genome sequences of the five Chinese strains were compared with the corresponding reference strains, respectively, according to the genospecies. Pairwise analysis demonstrates that there are only 70 SNPs between the genomes of CS4 and B31. However, there are many more SNPs between the genomes of QX-S13 and VS116, PD91 and PBi, FP1 and PKo, R9 and Pko, respectively. Gene comparison showed some important different genes. OspA was one of the important different genes. Comparative genomic studies have found that OspA gene sequences of PD91 and R9 had great differences compared with the sequence of B31. OspA gene sequence of R9 had a 96bp deletion; OspA gene of PD91 had two deletions: 9bp and 10 bp. To conclude, we showed the genomic characteristics of four genotype Chinese B. burgdorferi strains. The genomic sequence of B. yangtze sp. nov and differences from B. valaisiana were first reported. Comparative analysis of Chinese strains with the different Borrelia species from other areas will help us to understand evolution and pathogenesis of Chinese Borrelia burgdorferi strains.     

Comparative genomic analysis of Mycobacterium tuberculosis clinical isolates

Background

Due to excessive antibiotic use, drug-resistant Mycobacterium tuberculosis has become a serious public health threat and a major obstacle to disease control in many countries. To better understand the evolution of drug-resistant M. tuberculosis strains, we performed whole genome sequencing for 7 M. tuberculosis clinical isolates with different antibiotic resistance profiles and conducted comparative genomic analysis of gene variations among them.

Results

We observed that all 7 M. tuberculosis clinical isolates with different levels of drug resistance harbored similar numbers of SNPs, ranging from 1409–1464. The numbers of insertion/deletions (Indels) identified in the 7 isolates were also similar, ranging from 56 to 101. A total of 39 types of mutations were identified in drug resistance-associated loci, including 14 previously reported ones and 25 newly identified ones. Sixteen of the identified large Indels spanned PE-PPE-PGRS genes, which represents a major source of antigenic variability. Aside from SNPs and Indels, a CRISPR locus with varied spacers was observed in all 7 clinical isolates, suggesting that they might play an important role in plasticity of the M. tuberculosis genome. The nucleotide diversity (Л value) and selection intensity (dN/dS value) of the whole genome sequences of the 7 isolates were similar. The dN/dS values were less than 1 for all 7 isolates (range from 0.608885 to 0.637365), supporting the notion that M. tuberculosis genomes undergo purifying selection. The Л values and dN/dS values were comparable between drug-susceptible and drug-resistant strains.

Conclusions

In this study, we show that clinical M. tuberculosis isolates exhibit distinct variations in terms of the distribution of SNP, Indels, CRISPR-cas locus, as well as the nucleotide diversity and selection intensity, but there are no generalizable differences between drug-susceptible and drug-resistant isolates on the genomic scale. Our study provides evidence strengthening the notion that the evolution of drug resistance among clinical M. tuberculosis isolates is clearly a complex and diversified process.

Electronic supplementary material

The online version of this article (doi: 10.1186/1471-2164-15-469) contains supplementary material, which is available to authorized users.     

Genomic and Phenotypic Diversity of Coastal Vibrio cholerae Strains Is Linked to Environmental Factors

Studies of Vibrio cholerae diversity have focused primarily on pathogenic isolates of the O1 and O139 serotypes. However, autochthonous environmental isolates of this species routinely display more extensive genetic diversity than the primarily clonal pathogenic strains. In this study, genomic and metabolic profiles of 41 non-O1/O139 environmental isolates from central California coastal waters and four clinical strains are used to characterize the core genome and metabolome of V. cholerae. Comparative genome hybridization using microarrays constructed from the fully sequenced V. cholerae O1 El Tor N16961 genome identified 2,787 core genes that approximated the projected species core genome within 1.6%. Core genes are almost universally present in strains with widely different niches, suggesting that these genes are essential for persistence in diverse aquatic environments. In contrast, the dispensable genes and phenotypic traits identified in this study should provide increased fitness for certain niche environments. Environmental parameters, measured in situ during sample collection, are correlated to the presence of specific dispensable genes and metabolic capabilities, including utilization of mannose, sialic acid, citrate, and chitosan oligosaccharides. These results identify gene content and metabolic pathways that are likely selected for in certain coastal environments and may influence V. cholerae population structure in aquatic environments.     

Variation,Sex, and Social Cooperation: Molecular Population Genetics of the Social Amoeba Dictyostelium discoideum

Dictyostelium discoideum is a eukaryotic microbial model system for multicellular development, cell–cell signaling, and social behavior. Key models of social evolution require an understanding of genetic relationships between individuals across the genome or possibly at specific genes, but the nature of variation within D. discoideum is largely unknown. We re-sequenced 137 gene fragments in wild North American strains of D. discoideum and examined the levels and patterns of nucleotide variation in this social microbial species. We observe surprisingly low levels of nucleotide variation in D. discoideum across these strains, with a mean nucleotide diversity (π) of 0.08%, and no strong population stratification among North American strains. We also do not find any clear relationship between nucleotide divergence between strains and levels of social dominance and kin discrimination. Kin discrimination experiments, however, show that strains collected from the same location show greater ability to distinguish self from non-self than do strains from different geographic areas. This suggests that a greater ability to recognize self versus non-self may arise among strains that are more likely to encounter each other in nature, which would lead to preferential formation of fruiting bodies with clonemates and may prevent the evolution of cheating behaviors within D. discoideum populations. Finally, despite the fact that sex has rarely been observed in this species, we document a rapid decay of linkage disequilibrium between SNPs, the presence of recombinant genotypes among natural strains, and high estimates of the population recombination parameter ρ. The SNP data indicate that recombination is widespread within D. discoideum and that sex as a form of social interaction is likely to be an important aspect of the life cycle.     

Comparative Genome Analysis Provides Insights into the Pathogenicity of Flavobacterium psychrophilum

Flavobacterium psychrophilum is a fish pathogen in salmonid aquaculture worldwide that causes cold water disease (CWD) and rainbow trout fry syndrome (RTFS). Comparative genome analyses of 11 F. psychrophilum isolates representing temporally and geographically distant populations were used to describe the F. psychrophilum pan-genome and to examine virulence factors, prophages, CRISPR arrays, and genomic islands present in the genomes. Analysis of the genomic DNA sequences were complemented with selected phenotypic characteristics of the strains. The pan genome analysis showed that F. psychrophilum could hold at least 3373 genes, while the core genome contained 1743 genes. On average, 67 new genes were detected for every new genome added to the analysis, indicating that F. psychrophilum possesses an open pan genome. The putative virulence factors were equally distributed among isolates, independent of geographic location, year of isolation and source of isolates. Only one prophage-related sequence was found which corresponded to the previously described prophage 6H, and appeared in 5 out of 11 isolates. CRISPR array analysis revealed two different loci with dissimilar spacer content, which only matched one sequence in the database, the temperate bacteriophage 6H. Genomic Islands (GIs) were identified in F. psychrophilum isolates 950106-1/1 and CSF 259–93, associated with toxins and antibiotic resistance. Finally, phenotypic characterization revealed a high degree of similarity among the strains with respect to biofilm formation and secretion of extracellular enzymes. Global scale dispersion of virulence factors in the genomes and the abilities for biofilm formation, hemolytic activity and secretion of extracellular enzymes among the strains suggested that F. psychrophilum isolates have a similar mode of action on adhesion, colonization and destruction of fish tissues across large spatial and temporal scales of occurrence. Overall, the genomic characterization and phenotypic properties may provide new insights to the mechanisms of pathogenicity in F. psychrophilum.     

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.