Phylogeographic reconstruction of a bacterial species with high levels of lateral gene transfer

Whole genome sequence typing to investigate the Apophysomyces outbreak following a tornado in Joplin, Missouri, 2011

Case reports of Apophysomyces spp. in immunocompetent hosts have been a result of traumatic deep implantation of Apophysomyces spp. spore-contaminated soil or debris. On May 22, 2011 a tornado occurred in Joplin, MO, leaving 13 tornado victims with Apophysomyces trapeziformis infections as a result of lacerations from airborne material. We used whole genome sequence typing (WGST) for high-resolution phylogenetic SNP analysis of 17 outbreak Apophysomyces isolates and five additional temporally and spatially diverse Apophysomyces control isolates (three A. trapeziformis and two A. variabilis isolates). Whole genome SNP phylogenetic analysis revealed three clusters of genotypically related or identical A. trapeziformis isolates and multiple distinct isolates among the Joplin group; this indicated multiple genotypes from a single or multiple sources. Though no linkage between genotype and location of exposure was observed, WGST analysis determined that the Joplin isolates were more closely related to each other than to the control isolates, suggesting local population structure. Additionally, species delineation based on WGST demonstrated the need to reassess currently accepted taxonomic classifications of phylogenetic species within the genus Apophysomyces.

Reproductive clonality of pathogens: a perspective on pathogenic viruses, bacteria, fungi, and parasitic protozoa

We propose that clonal evolution in micropathogens be defined as restrained recombination on an evolutionary scale, with genetic exchange scarce enough to not break the prevalent pattern of clonal population structure, a definition already widely used for all kinds of pathogens, although not clearly formulated by many scientists and rejected by others. The two main manifestations of clonal evolution are strong linkage disequilibrium (LD) and widespread genetic clustering ("near-clading"). We hypothesize that this pattern is not mainly due to natural selection, but originates chiefly from in-built genetic properties of pathogens, which could be ancestral and could function as alternative allelic systems to recombination genes ("clonality/sexuality machinery") to escape recombinational load. The clonal framework of species of pathogens should be ascertained before any analysis of biomedical phenotypes (phylogenetic character mapping). In our opinion, this model provides a conceptual framework for the population genetics of any micropathogen.

Developing insights into the mechanisms of evolution of bacterial pathogens from whole-genome sequences

Evolution of bacterial pathogen populations has been detected in a variety of ways including phenotypic tests, such as metabolic activity, reaction to antisera and drug resistance and genotypic tests that measure variation in chromosome structure, repetitive loci and individual gene sequences. While informative, these methods only capture a small subset of the total variation and, therefore, have limited resolution. Advances in sequencing technologies have made it feasible to capture whole-genome sequence variation for each sample under study, providing the potential to detect all changes at all positions in the genome from single nucleotide changes to large-scale insertions and deletions. In this review, we focus on recent work that has applied this powerful new approach and summarize some of the advances that this has brought in our understanding of the details of how bacterial pathogens evolve.

Grappling with Proteus: population level approaches to understanding microbial diversity

The emerging fields of microbial population genetics and genomics provide an avenue to study the ecological rules that govern how communities form, function, and evolve. Our struggle to understand the causes and consequences of microbial diversity stems from our inability to define ecologically and evolutionarily meaningful units of diversity. The 16S rRNA-based tools that have been so useful in charting microbial diversity may lack sufficient sensitivity to answer many questions about the ecology and evolution of microbes. Examining genetic diversity with increased resolution is vital to understanding the forces shaping community structure. Population genetic analyses enabled by whole genome sequencing, multilocus sequence analyses, or single-nucleotide polymorphism analyses permit the testing of hypotheses pertaining to the geographic distribution, migration, and habitat preference of specific microbial lineages. Furthermore, these approaches can reveal patterns of gene exchange within and between populations and communities. Tools from microbial population genetics and population genomics can be used to increase the resolution with which we measure microbial diversity, enabling a focus on the scale of genetic diversity at which ecological processes impact evolutionary events. This tighter focus promises to improve our understanding of the causes and consequences of microbial community structure.

Out of the ground: aerial and exotic habitats of the melioidosis bacterium Burkholderia pseudomallei in grasses in Australia

Melioidosis is an emerging infectious disease of humans and animals in the tropics caused by the soil bacterium Burkholderia pseudomallei. Despite high fatality rates, the ecology of B.pseudomallei remains unclear. We used a combination of field and laboratory studies to investigate B.pseudomallei colonization of native and exotic grasses in northern Australia. Multivariable and spatial analyses were performed to determine significant predictors for B.pseudomallei occurrence in plants and soil collected longitudinally from field sites. In plant inoculation experiments, the impact of B.pseudomallei upon these grasses was studied and the bacterial load semi-quantified. Fluorescence in situ hybridization and confocal laser scanning microscopy were performed to localize the bacteria in plants. Burkholderia pseudomallei was found to inhabit not only the rhizosphere and roots but also aerial parts of specific grasses. This raises questions about the potential spread of B.pseudomallei by grazing animals whose droppings were found to be positive for these bacteria. In particular, B.pseudomallei readily colonized exotic grasses introduced to Australia for pasture. The ongoing spread of these introduced grasses creates new habitats suitable for B.pseudomallei survival and may be an important factor in the evolving epidemiology of melioidosis seen both in northern Australia and elsewhere globally.

High confidence prediction of essential genes in Burkholderia cenocepacia

Background

Essential genes are absolutely required for the survival of an organism. The identification of essential genes, besides being one of the most fundamental questions in biology, is also of interest for the emerging science of synthetic biology and for the development of novel antimicrobials. New antimicrobial therapies are desperately needed to treat multidrug-resistant pathogens, such as members of the Burkholderia cepacia complex.

Methodology/Principal Findings

We hypothesize that essential genes may be highly conserved within a group of evolutionary closely related organisms. Using a bioinformatics approach we determined that the core genome of the order Burkholderiales consists of 649 genes. All but two of these identified genes were located on chromosome 1 of Burkholderia cenocepacia. Although many of the 649 core genes of Burkholderiales have been shown to be essential in other bacteria, we were also able to identify a number of novel essential genes present mainly, or exclusively, within this order. The essentiality of some of the core genes, including the known essential genes infB, gyrB, ubiB, and valS, as well as the so far uncharacterized genes BCAL1882, BCAL2769, BCAL3142 and BCAL3369 has been confirmed experimentally in B. cenocepacia.

Conclusions/Significance

We report on the identification of essential genes using a novel bioinformatics strategy and provide bioinformatics and experimental evidence that the large majority of the identified genes are indeed essential. The essential genes identified here may represent valuable targets for the development of novel antimicrobials and their detailed study may shed new light on the functions required to support life.

Evaluation of recombinant proteins of Burkholderia mallei for serodiagnosis of glanders

Glanders is a contagious disease caused by the Gram-negative bacillus Burkholderia mallei. The number of equine glanders outbreaks has increased steadily during the last decade. The disease must be reported to the Office International des Epizooties, Paris, France. Glanders serodiagnosis is hampered by the considerable number of false positives and negatives of the internationally prescribed tests. The major problem leading to the low sensitivity and specificity of the complement fixation test (CFT) and enzyme-linked immunosorbent assay (ELISA) has been linked to the test antigens currently used, i.e., crude preparations of whole cells. False-positive results obtained from other diagnostic tests utilizing crude antigens lead to financial losses to animal owners, and false-negative results can turn a risk into a possible threat. In this study, we report on the identification of diagnostic targets using bioinformatics tools for serodiagnosis of glanders. The identified gene sequences were cloned and expressed as recombinant proteins. The purified recombinant proteins of B. mallei were used in an indirect ELISA format for serodiagnosis of glanders. Two recombinant proteins, 0375H and 0375TH, exhibited 100% sensitivity and specificity for glanders diagnosis. The proteins also did not cross-react with sera from patients with the closely related disease melioidosis. The results of this investigation highlight the potential of recombinant 0375H and 0375TH proteins in specific and sensitive diagnosis of glanders.

Development and validation of Burkholderia pseudomallei-specific real-time PCR assays for clinical, environmental or forensic detection applications

The bacterium Burkholderia pseudomallei causes melioidosis, a rare but serious illness that can be fatal if untreated or misdiagnosed. Species-specific PCR assays provide a technically simple method for differentiating B. pseudomallei from near-neighbor species. However, substantial genetic diversity and high levels of recombination within this species reduce the likelihood that molecular signatures will differentiate all B. pseudomallei from other Burkholderiaceae. Currently available molecular assays for B. pseudomallei detection lack rigorous validation across large in silico datasets and isolate collections to test for specificity, and none have been subjected to stringent quality control criteria (accuracy, precision, selectivity, limit of quantitation (LoQ), limit of detection (LoD), linearity, ruggedness and robustness) to determine their suitability for environmental, clinical or forensic investigations. In this study, we developed two novel B. pseudomallei specific assays, 122018 and 266152, using a dual-probe approach to differentiate B. pseudomallei from B. thailandensis, B. oklahomensis and B. thailandensis-like species; other species failed to amplify. Species specificity was validated across a large DNA panel (>2,300 samples) comprising Burkholderia spp. and non-Burkholderia bacterial and fungal species of clinical and environmental relevance. Comparison of assay specificity to two previously published B. pseudomallei-specific assays, BurkDiff and TTS1, demonstrated comparable performance of all assays, providing between 99.7 and 100% specificity against our isolate panel. Last, we subjected 122018 and 266152 to rigorous quality control analyses, thus providing quantitative limits of assay performance. Using B. pseudomallei as a model, our study provides a framework for comprehensive quantitative validation of molecular assays and provides additional, highly validated B. pseudomallei assays for the scientific research community.

Genomic and SNP analyses demonstrate a distant separation of the hospital and community-associated clades of Enterococcus faecium

Recent studies have pointed to the existence of two subpopulations of Enterococcus faecium, one containing primarily commensal/community-associated (CA) strains and one that contains most clinical or hospital-associated (HA) strains, including those classified by multi-locus sequence typing (MLST) as belonging to the CC17 group. The HA subpopulation more frequently has IS16, pathogenicity island(s), and plasmids or genes associated with antibiotic resistance, colonization, and/or virulence. Supporting the two clades concept, we previously found a 3–10% difference between four genes from HA-clade strains vs. CA-clade strains, including 5% difference between pbp5-R of ampicillin-resistant, HA strains and pbp5-S of ampicillin-sensitive, CA strains. To further investigate the core genome of these subpopulations, we studied 100 genes from 21 E. faecium genome sequences; our analyses of concatenated sequences, SNPs, and individual genes all identified two distinct groups. With the concatenated sequence, HA-clade strains differed by 0–1% from one another while CA clade strains differed from each other by 0–1.1%, with 3.5–4.2% difference between the two clades. While many strains had a few genes that grouped in one clade with most of their genes in the other clade, one strain had 28% of its genes in the CA clade and 72% in the HA clade, consistent with the predicted role of recombination in the evolution of E. faecium. Using estimates for Escherichia coli, molecular clock calculations using sSNP analysis indicate that these two clades may have diverged ≥1 million years ago or, using the higher mutation rate for Bacillus anthracis, ∼300,000 years ago. These data confirm the existence of two clades of E. faecium and show that the differences between the HA and CA clades occur at the core genomic level and long preceded the modern antibiotic era.

Burkholderia pseudomallei in unchlorinated domestic bore water, Tropical Northern Australia

To determine whether unchlorinated bore water in northern Australia contained Burkholderia pseudomallei organisms, we sampled 55 bores; 18 (33%) were culture positive. Multilocus sequence typing identified 15 sequence types. The B. pseudomallei sequence type from 1 water sample matched a clinical isolate from a resident with melioidosis on the same property.

Miocene dispersal drives island radiations in the palm tribe Trachycarpeae (Arecaceae)

The study of three island groups of the palm tribe Trachycarpeae (Arecaceae/Palmae) permits both the analysis of each independent radiation and comparisons across the tribe to address general processes that drive island diversification. Phylogenetic relationships of Trachycarpeae were inferred from three plastid and three low-copy nuclear genes. The incongruent topological position of Brahea in CISP5 was hypothesized to be caused by a gene duplication event and was addressed using uninode coding. The resulting phylogenetic trees were well-resolved and the genera were all highly supported except for Johannesteijsmannia and Serenoa. Divergence time analysis estimated the stem of the tribe to be approximately 86 Ma and the crown to be 38 Ma, indicating that significant extinction may have occurred along this branch. Historical biogeographic analysis suggested that Trachycarpeae are of southern North American, Central American, or Caribbean origin and supports previous hypotheses of a Laurasian origin. The biogeography and disjunctions within the tribe were interpreted with respect to divergence times, the fossil record, and geological factors such as the formation of the Greater Antilles--Aves Ridge, the Bering and the North Atlantic land bridges, tectonic movement in Southeast Asia, climatic shifts between the Eocene and Pliocene, and volcanism in the Pacific basin. In considering the three major island radiations within Trachycarpeae, Miocene dispersal appears to have been the driving force in allopatric speciation and is highlighted here as an emerging pattern across the tree of life.

Diversity of 16S-23S rDNA internal transcribed spacer (ITS) reveals phylogenetic relationships in Burkholderia pseudomallei and its near-neighbors

Length polymorphisms within the 16S-23S ribosomal DNA internal transcribed spacer (ITS) have been described as stable genetic markers for studying bacterial phylogenetics. In this study, we used these genetic markers to investigate phylogenetic relationships in Burkholderia pseudomallei and its near-relative species. B. pseudomallei is known as one of the most genetically recombined bacterial species. In silico analysis of multiple B. pseudomallei genomes revealed approximately four homologous rRNA operons and ITS length polymorphisms therein. We characterized ITS distribution using PCR and analyzed via a high-throughput capillary electrophoresis in 1,191 B. pseudomallei strains. Three major ITS types were identified, two of which were commonly found in most B. pseudomallei strains from the endemic areas, whereas the third one was significantly correlated with worldwide sporadic strains. Interestingly, mixtures of the two common ITS types were observed within the same strains, and at a greater incidence in Thailand than Australia suggesting that genetic recombination causes the ITS variation within species, with greater recombination frequency in Thailand. In addition, the B. mallei ITS type was common to B. pseudomallei, providing further support that B. mallei is a clone of B. pseudomallei. Other B. pseudomallei near-neighbors possessed unique and monomorphic ITS types. Our data shed light on evolutionary patterns of B. pseudomallei and its near relative species.

Epidemiological tracking and population assignment of the non-clonal bacterium, Burkholderia pseudomallei

Rapid assignment of bacterial pathogens into predefined populations is an important first step for epidemiological tracking. For clonal species, a single allele can theoretically define a population. For non-clonal species such as Burkholderia pseudomallei, however, shared allelic states between distantly related isolates make it more difficult to identify population defining characteristics. Two distinct B. pseudomallei populations have been previously identified using multilocus sequence typing (MLST). These populations correlate with the major foci of endemicity (Australia and Southeast Asia). Here, we use multiple Bayesian approaches to evaluate the compositional robustness of these populations, and provide assignment results for MLST sequence types (STs). Our goal was to provide a reference for assigning STs to an established population without the need for further computational analyses. We also provide allele frequency results for each population to enable estimation of population assignment even when novel STs are discovered. The ability for humans and potentially contaminated goods to move rapidly across the globe complicates the task of identifying the source of an infection or outbreak. Population genetic dynamics of B. pseudomallei are particularly complicated relative to other bacterial pathogens, but the work here provides the ability for broad scale population assignment. As there is currently no independent empirical measure of successful population assignment, we provide comprehensive analytical details of our comparisons to enable the reader to evaluate the robustness of population designations and assignments as they pertain to individual research questions. Finer scale subdivision and verification of current population compositions will likely be possible with genotyping data that more comprehensively samples the genome. The approach used here may be valuable for other non-clonal pathogens that lack simple group-defining genetic characteristics and provides a rapid reference for epidemiologists wishing to track the origin of infection without the need to compile population data and learn population assignment algorithms.

Burkholderia pseudomallei is a soil-dwelling bacterium that can infect a large range of hosts. In humans, B. pseudomallei causes melioidosis, and typical routes of entry include open wounds, inhalation, or ingestion. Clinical features are diverse, although pneumonia and abscess formation are common. High rates of recombination within the genome of this bacterium have confounded attempts to match clinical samples to geographically defined populations. Here we provide a reference that simplifies source attribution issues. We applied population assignment software to previously generated sequence data from seven B. pseudomallei genes to define the major geographic populations within this species. We evaluated the robustness of our results by comparison with two additional population assignment programs. We present the likelihood that each variant is assigned to a particular geographic population. This information can be used to assign novel B. pseudomallei isolates to a geographic population without needing to learn and run cumbersome population assignment applications. This method can also be used for other bacteria that are difficult to source-attribute due to high levels of genomic variation and recombination.

Whole genome sequence analysis of Cryptococcus gattii from the Pacific Northwest reveals unexpected diversity

A recent emergence of Cryptococcus gattii in the Pacific Northwest involves strains that fall into three primarily clonal molecular subtypes: VGIIa, VGIIb and VGIIc. Multilocus sequence typing (MLST) and variable number tandem repeat analysis appear to identify little diversity within these molecular subtypes. Given the apparent expansion of these subtypes into new geographic areas and their ability to cause disease in immunocompetent individuals, differentiation of isolates belonging to these subtypes could be very important from a public health perspective. We used whole genome sequence typing (WGST) to perform fine-scale phylogenetic analysis on 20 C. gattii isolates, 18 of which are from the VGII molecular type largely responsible for the Pacific Northwest emergence. Analysis both including and excluding (289,586 SNPs and 56,845 SNPs, respectively) molecular types VGI and VGIII isolates resulted in phylogenetic reconstructions consistent, for the most part, with MLST analysis but with far greater resolution among isolates. The WGST analysis presented here resulted in identification of over 100 SNPs among eight VGIIc isolates as well as unique genotypes for each of the VGIIa, VGIIb and VGIIc isolates. Similar levels of genetic diversity were found within each of the molecular subtype isolates, despite the fact that the VGIIb clade is thought to have emerged much earlier. The analysis presented here is the first multi-genome WGST study to focus on the C. gattii molecular subtypes involved in the Pacific Northwest emergence and describes the tools that will further our understanding of this emerging pathogen.

Molecular investigations of a locally acquired case of melioidosis in Southern AZ, USA

Melioidosis is caused by Burkholderia pseudomallei, a Gram-negative bacillus, primarily found in soils in Southeast Asia and northern Australia. A recent case of melioidosis in non-endemic Arizona was determined to be the result of locally acquired infection, as the patient had no travel history to endemic regions and no previous history of disease. Diagnosis of the case was confirmed through multiple microbiologic and molecular techniques. To enhance the epidemiological analysis, we conducted several molecular genotyping procedures, including multi-locus sequence typing, SNP-profiling, and whole genome sequence typing. Each technique has different molecular epidemiologic advantages, all of which provided evidence that the infecting strain was most similar to those found in Southeast Asia, possibly originating in, or around, Malaysia. Advancements in new typing technologies provide genotyping resolution not previously available to public health investigators, allowing for more accurate source identification.

Melioidosis is a bacterial disease caused by percutaneous inoculation, aspiration or ingestion of the soil bacteria Burkholderia pseudomallei. Melioidosis is primarily found in Southeast Asia and Northern Australia, and, to a lesser degree, nearby regions. A recent case of melioidosis in Southwestern United States (Southern Arizona) prompted a detailed epidemiological and molecular investigation to discover the source of infection. The authors describe the use of multiple genomic analysis tools to aid in this investigation. The results of these analyses uniformly identified Southeast Asia as the source of the strain that infected the patient, however the epidemiological investigation had determined the patient had no international travel or known exposures to Southeast Asian products. New cutting edge technologies, such as next generation sequencing, are quickly being adapted into epidemiologic investigations, particularly for cases and outbreaks of unknown origin, although older, mature technologies with larger existing databases will still be needed for appropriate comparative analyses.

Molecular characterization of clinical Burkholderia pseudomallei isolates from India

Multilocus sequence typing of seven isolates of Burkholderia pseudomallei from India showed considerable diversity, with six different sequence types. Possible dissemination of melioidosis by historical trading routes is supported by links to strains from Southeast Asia, China, and Africa and the presence of the Burkholderia mallei allele of the bimA gene.

SNIT: SNP identification for strain typing

With ever-increasing numbers of microbial genomes being sequenced, efficient tools are needed to perform strain-level identification of any newly sequenced genome. Here, we present the SNP identification for strain typing (SNIT) pipeline, a fast and accurate software system that compares a newly sequenced bacterial genome with other genomes of the same species to identify single nucleotide polymorphisms (SNPs) and small insertions/deletions (indels). Based on this information, the pipeline analyzes the polymorphic loci present in all input genomes to identify the genome that has the fewest differences with the newly sequenced genome. Similarly, for each of the other genomes, SNIT identifies the input genome with the fewest differences. Results from five bacterial species show that the SNIT pipeline identifies the correct closest neighbor with 75% to 100% accuracy. The SNIT pipeline is available for download at http://www.bhsai.org/snit.html.

The aftermath of the Western Australian melioidosis outbreak

Melioidosis became a notifiable disease in Western Australia (WA) 2 years after the West Kimberley melioidosis outbreak. Two cases of melioidosis caused by the outbreak genotype of Burkholderia pseudomallei (National Collection of Type Cultures [NCTC] 13177) occurred in 1998 and 1999 in persons who visited the outbreak location at the time. No other infections caused by the outbreak strain have been recorded in WA since that time, despite an average of four culture-positive cases per year. Sporadic cases of melioidosis often follow tropical storms and cyclones during summer, and they have been detected outside the endemic area when cyclones travel far inland. In 2007, environmental isolates resembling NCTC 13177 were found 500 km east of the outbreak location after unusually severe weather. Recent whole-genome analysis places NCTC 13177 genetically close to other Australian isolates. Additional biogeographic and ecological studies are needed to establish the relative importance of environmental cofactors in disease pathogenesis.

A very early-branching Staphylococcus aureus lineage lacking the carotenoid pigment staphyloxanthin

Here we discuss the evolution of the northern Australian Staphylococcus aureus isolate MSHR1132 genome. MSHR1132 belongs to the divergent clonal complex 75 lineage. The average nucleotide divergence between orthologous genes in MSHR1132 and typical S. aureus is approximately sevenfold greater than the maximum divergence observed in this species to date. MSHR1132 has a small accessory genome, which includes the well-characterized genomic islands, νSAα and νSaβ, suggesting that these elements were acquired well before the expansion of the typical S. aureus population. Other mobile elements show mosaic structure (the prophage φSa3) or evidence of recent acquisition from a typical S. aureus lineage (SCCmec, ICE6013 and plasmid pMSHR1132). There are two differences in gene repertoire compared with typical S. aureus that may be significant clues as to the genetic basis underlying the successful emergence of S. aureus as a pathogen. First, MSHR1132 lacks the genes for production of staphyloxanthin, the carotenoid pigment that confers upon S. aureus its characteristic golden color and protects against oxidative stress. The lack of pigment was demonstrated in 126 of 126 CC75 isolates. Second, a mobile clustered regularly interspaced short palindromic repeat (CRISPR) element is inserted into orfX of MSHR1132. Although common in other staphylococcal species, these elements are very rare within S. aureus and may impact accessory genome acquisition. The CRISPR spacer sequences reveal a history of attempted invasion by known S. aureus mobile elements. There is a case for the creation of a new taxon to accommodate this and related isolates.

Dissection of the Burkholderia intracellular life cycle using a photothermal nanoblade

Burkholderia pseudomallei and Burkholderia thailandensis are related pathogens that invade a variety of cell types, replicate in the cytoplasm, and spread to nearby cells. We have investigated temporal and spatial requirements for virulence determinants in the intracellular life cycle, using genetic dissection and photothermal nanoblade delivery, which allows efficient placement of bacterium-sized cargo into the cytoplasm of mammalian cells. The conserved Bsa type III secretion system (T3SS(Bsa)) is dispensable for invasion, but is essential for escape from primary endosomes. By nanoblade delivery of B. thailandensis we demonstrate that all subsequent events in intercellular spread occur independently of T3SS(Bsa) activity. Although intracellular movement was essential for cell-cell spread by B. pseudomallei and B. thailandensis, neither BimA-mediated actin polymerization nor the formation of membrane protrusions containing bacteria was required for B. thailandensis. Surprisingly, the cryptic (fla2) flagellar system encoded on chromosome 2 of B. thailandensis supported rapid intracellular motility and efficient cell-cell spread. Plaque formation by both pathogens was dependent on the activity of a type VI secretion system (T6SS-1) that functions downstream from T3SS(Bsa)-mediated endosome escape. A remarkable feature of Burkholderia is their ability to induce the formation of multinucleate giant cells (MNGCs) in multiple cell types. By infection and nanoblade delivery, we observed complete correspondence between mutant phenotypes in assays for cell fusion and plaque formation, and time-course studies showed that plaque formation represents MNGC death. Our data suggest that the primary means for intercellular spread involves cell fusion, as opposed to pseudopod engulfment and bacterial escape from double-membrane vacuoles.

The utility of high-resolution melting analysis of SNP nucleated PCR amplicons--an MLST based Staphylococcus aureus typing scheme

High resolution melting (HRM) analysis is gaining prominence as a method for discriminating DNA sequence variants. Its advantage is that it is performed in a real-time PCR device, and the PCR amplification and HRM analysis are closed tube, and effectively single step. We have developed an HRM-based method for Staphylococcus aureus genotyping. Eight single nucleotide polymorphisms (SNPs) were derived from the S. aureus multi-locus sequence typing (MLST) database on the basis of maximized Simpson's Index of Diversity. Only G↔A, G↔T, C↔A, C↔T SNPs were considered for inclusion, to facilitate allele discrimination by HRM. In silico experiments revealed that DNA fragments incorporating the SNPs give much higher resolving power than randomly selected fragments. It was shown that the predicted optimum fragment size for HRM analysis was 200 bp, and that other SNPs within the fragments contribute to the resolving power. Six DNA fragments ranging from 83 bp to 219 bp, incorporating the resolution optimized SNPs were designed. HRM analysis of these fragments using 94 diverse S. aureus isolates of known sequence type or clonal complex (CC) revealed that sequence variants are resolved largely in accordance with G+C content. A combination of experimental results and in silico prediction indicates that HRM analysis resolves S. aureus into 268 “melt types” (MelTs), and provides a Simpson's Index of Diversity of 0.978 with respect to MLST. There is a high concordance between HRM analysis and the MLST defined CCs. We have generated a Microsoft Excel key which facilitates data interpretation and translation between MelT and MLST data. The potential of this approach for genotyping other bacterial pathogens was investigated using a computerized approach to estimate the densities of SNPs with unlinked allelic states. The MLST databases for all species tested contained abundant unlinked SNPs, thus suggesting that high resolving power is not dependent upon large numbers of SNPs.

Next-generation sequencing of Coccidioides immitis isolated during cluster investigation

Next-generation sequencing enables use of whole-genome sequence typing (WGST) as a viable and discriminatory tool for genotyping and molecular epidemiologic analysis. We used WGST to confirm the linkage of a cluster of Coccidioides immitis isolates from 3 patients who received organ transplants from a single donor who later had positive test results for coccidioidomycosis. Isolates from the 3 patients were nearly genetically identical (a total of 3 single-nucleotide polymorphisms identified among them), thereby demonstrating direct descent of the 3 isolates from an original isolate. We used WGST to demonstrate the genotypic relatedness of C. immitis isolates that were also epidemiologically linked. Thus, WGST offers unique benefits to public health for investigation of clusters considered to be linked to a single source.

Population genomics of Sinorhizobium medicae based on low-coverage sequencing of sympatric isolates

We investigated the genomic diversity of a local population of the symbiotic bacterium Sinorhizobium medicae, isolated from the roots of wild Medicago lupulina plants, in order to assess genomic diversity, to identify genomic regions influenced by duplication, deletion or strong selection, and to explore the composition of the pan-genome. Partial genome sequences of 12 isolates were obtained by Roche 454 shotgun sequencing (average 5.3 Mb per isolate) and compared with the published sequence of S. medicae WSM 419. Homologous recombination appears to have less impact on the polymorphism patterns of the chromosome than on the chromid pSMED01 and megaplasmid pSMED02. Moreover, pSMED02 is a hot spot of insertions and deletions. The whole chromosome is characterized by low sequence polymorphism, consistent with the high density of housekeeping genes. Similarly, the level of polymorphism of symbiosis genes (low) and of genes involved in polysaccharide synthesis (high) may reflect different selection. Finally, some isolates carry genes that may confer adaptations that S. medicae WSM 419 lacks, including homologues of genes encoding rhizobitoxine synthesis, iron uptake, response to autoinducer-2, and synthesis of distinct polysaccharides. The presence or absence of these genes was confirmed by PCR in each of these 12 isolates and a further 27 isolates from the same population. All isolates had rhizobitoxine genes, while the other genes were co-distributed, suggesting that they may be on the same mobile element. These results are discussed in relation to the ecology of Medicago symbionts and in the perspective of population genomics studies.

Lack of phylogeographic structure in the freshwater cyanobacterium Microcystis aeruginosa suggests global dispersal

BACKGROUND

Free-living microorganisms have long been assumed to have ubiquitous distributions with little biogeographic signature because they typically exhibit high dispersal potential and large population sizes. However, molecular data provide contrasting results and it is far from clear to what extent dispersal limitation determines geographic structuring of microbial populations. We aimed to determine biogeographical patterns of the bloom-forming freshwater cyanobacterium Microcystis aeruginosa. Being widely distributed on a global scale but patchily on a regional scale, this prokaryote is an ideal model organism to study microbial dispersal and biogeography.

METHODOLOGY/PRINCIPAL FINDINGS

The phylogeography of M. aeruginosa was studied based on a dataset of 311 rDNA internal transcribed spacer (ITS) sequences sampled from six continents. Richness of ITS sequences was high (239 ITS types were detected). Genetic divergence among ITS types averaged 4% (maximum pairwise divergence was 13%). Preliminary analyses revealed nearly completely unresolved phylogenetic relationships and a lack of genetic structure among all sequences due to extensive homoplasy at multiple hypervariable sites. After correcting for this, still no clear phylogeographic structure was detected, and no pattern of isolation by distance was found on a global scale. Concomitantly, genetic differentiation among continents was marginal, whereas variation within continents was high and was mostly shared with all other continents. Similarly, no genetic structure across climate zones was detected.

CONCLUSIONS/SIGNIFICANCE

The high overall diversity and wide global distribution of common ITS types in combination with the lack of phylogeographic structure suggest that intercontinental dispersal of M. aeruginosa ITS types is not rare, and that this species might have a truly cosmopolitan distribution.

Molecular phylogeny of Burkholderia pseudomallei from a remote region of Papua New Guinea

Background

The island of New Guinea is located midway between the world's two major melioidosis endemic regions of Australia and Southeast Asia. Previous studies in Papua New Guinea have demonstrated autochthonous melioidosis in Balimo, Western province. In contrast to other regions of endemicity, isolates recovered from both environmental and clinical sources demonstrate narrow genetic diversity over large spatial and temporal scales.

Methodology/Principal Findings

We employed molecular typing techniques to determine the phylogenetic relationships of these isolates to each other and to others worldwide to aid in understanding the origins of the Papua New Guinean isolates. Multi-locus sequence typing of the 39 isolates resolved three unique sequence types. Phylogenetic reconstruction and Structure analysis determined that all isolates were genetically closer to those from Australia than those from Southeast Asia. Gene cluster analysis however, identified a Yersinia-like fimbrial gene cluster predominantly found among Burkholderia pseudomallei derived from Southeast Asia. Higher resolution VNTR typing and phylogenetic reconstruction of the Balimo isolates resolved 24 genotypes with long branch lengths. These findings are congruent with long term persistence in the region and a high level of environmental stability.

Conclusions/Significance

Given that anthropogenic influence has been hypothesized as a mechanism for the dispersal of B. pseudomallei, these findings correlate with limited movement of the indigenous people in the region. The palaeogeographical and anthropogenic history of Australasia and the results from this study indicate that New Guinea is an important region for the further study of B. pseudomallei origins and dissemination.

The epidemiology and clinical spectrum of melioidosis: 540 cases from the 20 year Darwin prospective study

Background

Over 20 years, from October 1989, the Darwin prospective melioidosis study has documented 540 cases from tropical Australia, providing new insights into epidemiology and the clinical spectrum.

Principal Findings

The principal presentation was pneumonia in 278 (51%), genitourinary infection in 76 (14%), skin infection in 68 (13%), bacteremia without evident focus in 59 (11%), septic arthritis/osteomyelitis in 20 (4%) and neurological melioidosis in 14 (3%). 298 (55%) were bacteremic and 116 (21%) developed septic shock (58 fatal). Internal organ abscesses and secondary foci in lungs and/or joints were common. Prostatic abscesses occurred in 76 (20% of 372 males). 96 (18%) had occupational exposure to Burkholderia pseudomallei. 118 (22%) had a specific recreational or occupational incident considered the likely infecting event. 436 (81%) presented during the monsoonal wet season. The higher proportion with pneumonia in December to February supports the hypothesis of infection by inhalation during severe weather events. Recurrent melioidosis occurred in 29, mostly attributed to poor adherence to therapy. Mortality decreased from 30% in the first 5 years to 9% in the last five years (p<0.001). Risk factors for melioidosis included diabetes (39%), hazardous alcohol use (39%), chronic lung disease (26%) and chronic renal disease (12%). There was no identifiable risk factor in 20%. Of the 77 fatal cases (14%), 75 had at least one risk factor; the other 2 were elderly. On multivariate analysis of risk factors, age, location and season, the only independent predictors of mortality were the presence of at least one risk factor (OR 9.4; 95% CI 2.3–39) and age ≥50 years (OR 2.0; 95% CI 1.2–2.3).

Conclusions

Melioidosis should be seen as an opportunistic infection that is unlikely to kill a healthy person, provided infection is diagnosed early and resources are available to provide appropriate antibiotics and critical care.

Melioidosis is an occupationally and recreationally acquired infection important in Southeast Asia and northern Australia. Recently cases have been reported from more diverse locations globally. The responsible bacterium, Burkholderia pseudomallei, is considered a potential biothreat agent. Risk factors predisposing to melioidosis are well recognised, most notably diabetes. The Darwin prospective melioidosis study has identified 540 cases of melioidosis over 20 years and analysis of the epidemiology and clinical findings provides important new insights into this disease. Risk factors identified in addition to diabetes, hazardous alcohol use and chronic renal disease include chronic lung disease, malignancies, rheumatic heart disease, cardiac failure and age ≥50 years. Half of patients presented with pneumonia and septic shock was common (21%). The decrease in mortality from 30% in the first 5 years of the study to 9% in the last five years is attributed to earlier diagnosis and improvements in intensive care management. Of the 77 fatal cases (14%), all had known risk factors for melioidosis. This supports the most important conclusion of the study, which is that melioidosis is very unlikely to kill a healthy person, provided the infection is diagnosed early and resources are available to provide appropriate antibiotics and critical care where required.

BurkDiff: a real-time PCR allelic discrimination assay for Burkholderia pseudomallei and B. mallei

A real-time PCR assay, BurkDiff, was designed to target a unique conserved region in the B. pseudomallei and B. mallei genomes containing a SNP that differentiates the two species. Sensitivity and specificity were assessed by screening BurkDiff across 469 isolates of B. pseudomallei, 49 isolates of B. mallei, and 390 isolates of clinically relevant non-target species. Concordance of results with traditional speciation methods and no cross-reactivity to non-target species show BurkDiff is a robust, highly validated assay for the detection and differentiation of B. pseudomallei and B. mallei.

High-redundancy draft sequencing of 15 clinical and environmental Burkholderia strains

The Gram-negative Burkholderia genus includes several species of intracellular bacterial pathogens that pose substantial risk to humans. In this study, we have generated draft genome sequences of 15 strains of B. oklahomensis, B. pseudomallei, B. thailandensis, and B. ubonensis to an average sequence read coverage of 25- to 40-fold.

A genomic survey of positive selection in Burkholderia pseudomallei provides insights into the evolution of accidental virulence

Certain environmental microorganisms can cause severe human infections, even in the absence of an obvious requirement for transition through an animal host for replication ("accidental virulence"). To understand this process, we compared eleven isolate genomes of Burkholderia pseudomallei (Bp), a tropical soil microbe and causative agent of the human and animal disease melioidosis. We found evidence for the existence of several new genes in the Bp reference genome, identifying 282 novel genes supported by at least two independent lines of supporting evidence (mRNA transcripts, database homologs, and presence of ribosomal binding sites) and 81 novel genes supported by all three lines. Within the Bp core genome, 211 genes exhibited significant levels of positive selection (4.5%), distributed across many cellular pathways including carbohydrate and secondary metabolism. Functional experiments revealed that certain positively selected genes might enhance mammalian virulence by interacting with host cellular pathways or utilizing host nutrients. Evolutionary modifications improving Bp environmental fitness may thus have indirectly facilitated the ability of Bp to colonize and survive in mammalian hosts. These findings improve our understanding of the pathogenesis of melioidosis, and establish Bp as a model system for studying the genetics of accidental virulence.

Within-host evolution of Burkholderia pseudomallei in four cases of acute melioidosis

Little is currently known about bacterial pathogen evolution and adaptation within the host during acute infection. Previous studies of Burkholderia pseudomallei, the etiologic agent of melioidosis, have shown that this opportunistic pathogen mutates rapidly both in vitro and in vivo at tandemly repeated loci, making this organism a relevant model for studying short-term evolution. In the current study, B. pseudomallei isolates cultured from multiple body sites from four Thai patients with disseminated melioidosis were subjected to fine-scale genotyping using multilocus variable-number tandem repeat analysis (MLVA). In order to understand and model the in vivo variable-number tandem repeat (VNTR) mutational process, we characterized the patterns and rates of mutations in vitro through parallel serial passage experiments of B. pseudomallei. Despite the short period of infection, substantial divergence from the putative founder genotype was observed in all four melioidosis cases. This study presents a paradigm for examining bacterial evolution over the short timescale of an acute infection. Further studies are required to determine whether the mutational process leads to phenotypic alterations that impact upon bacterial fitness in vivo. Our findings have important implications for future sampling strategies, since colonies in a single clinical sample may be genetically heterogeneous, and organisms in a culture taken late in the infective process may have undergone considerable genetic change compared with the founder inoculum.