Search

Staphylococcus aureus is a Gram-positive pathogen responsible for tremendous morbidity and mortality. As with most bacteria, S. aureus requires iron to cause disease, and it can acquire iron from host hemoglobin. The current model for staphylococcal hemoglobin-iron acquisition proposes that S. aureus binds hemoglobin through the surface-exposed hemoglobin receptor IsdB. IsdB removes heme from bound hemoglobin and transfers this cofactor to other proteins of the Isd system, which import and degrade heme to release iron in the cytoplasm. Here we demonstrate that the individual components of the Isd system are required for growth on low nanomolar concentrations of hemoglobin as a sole source of iron. An in-depth study of hemoglobin binding by IsdB revealed key residues that are required for hemoglobin binding. Further, we show that these residues are necessary for heme extraction from hemoglobin and growth on hemoglobin as a sole iron source. These processes are found to contribute to the pathogenicity of S. aureus in a murine model of infection. Together these results build on the model for Isd-mediated hemoglobin binding and heme-iron acquisition during the pathogenesis of S. aureus infection.

Microbial virulence is a complex and often multifactorial phenotype, intricately linked to a pathogen's evolutionary trajectory. Toxicity, the ability to destroy host cell membranes, and adhesion, the ability to adhere to human tissues, are the major virulence factors of many bacterial pathogens, including Staphylococcus aureus. Here, we assayed the toxicity and adhesiveness of 90 MRSA (methicillin resistant S. aureus) isolates and found that while there was remarkably little variation in adhesion, toxicity varied by over an order of magnitude between isolates, suggesting different evolutionary selection pressures acting on these two traits. We performed a genome-wide association study (GWAS) and identified a large number of loci, as well as a putative network of epistatically interacting loci, that significantly associated with toxicity. Despite this apparent complexity in toxicity regulation, a predictive model based on a set of significant single nucleotide polymorphisms (SNPs) and insertion and deletions events (indels) showed a high degree of accuracy in predicting an isolate's toxicity solely from the genetic signature at these sites. Our results thus highlight the potential of using sequence data to determine clinically relevant parameters and have further implications for understanding the microbial virulence of this opportunistic pathogen.

Vancomycin-intermediate Staphylococcus aureus (VISA) is currently defined as having minimal inhibitory concentration (MIC) of 4-8 µg/ml. VISA evolves through changes in multiple genetic loci with at least 16 candidate genes identified in clinical and in vitro-selected VISA strains. We report a whole-genome comparative analysis of 49 vancomycin-sensitive S. aureus and 26 VISA strains. Resistance to vancomycin was determined by broth microdilution, Etest, and population analysis profile-area under the curve (PAP-AUC). Genome-wide association studies (GWAS) of 55,977 single-nucleotide polymorphisms identified in one or more strains found one highly significant association (P = 8.78 E-08) between a nonsynonymous mutation at codon 481 (H481) of the rpoB gene and increased vancomycin MIC. Additionally, we used a database of public S. aureus genome sequences to identify rare mutations in candidate genes associated with VISA. On the basis of these data, we proposed a preliminary model called ECM+RMCG for the VISA phenotype as a benchmark for future efforts. The model predicted VISA based on the presence of a rare mutation in a set of candidate genes (walKR, vraSR, graSR, and agrA) and/or three previously experimentally verified mutations (including the rpoB H481 locus) with an accuracy of 81% and a sensitivity of 73%. Further, the level of resistance measured by both Etest and PAP-AUC regressed positively with the number of mutations present in a strain. This study demonstrated 1) the power of GWAS for identifying common genetic variants associated with antibiotic resistance in bacteria and 2) that rare mutations in candidate gene, identified using large genomic data sets, can also be associated with resistance phenotypes., (© The Author(s) 2014. Published by Oxford University Press on behalf of the Society for Molecular Biology and Evolution.)

The whale shark (Rhincodon typus) is the largest extant species of fish, belonging to the order Orectolobiformes. It is listed as a "vulnerable" species on the International Union for Conservation of Nature (IUCN)'s Red List of Threatened Species, which makes it an important species for conservation efforts. We report here the first complete sequence of the mitochondrial genome (mitogenome) of the whale shark obtained by next-generation sequencing methods. The assembled mitogenome is a 16,875 bp circle, comprising of 13 protein-coding genes, two rRNA genes, 22 tRNA genes and a control region. We also performed comparative analysis of the whale shark mitogenome to the available mitogenome sequences of 17 other shark species, four from the order Orectolobiformes, five from Lamniformes and eight from Carcharhiniformes. The nucleotide composition, number and arrangement of the genes in whale shark mitogenome are the same as found in the mitogenomes of the other members of the order Orectolobiformes and its closest orders Lamniformes and Carcharhiniformes, although the whale shark mitogenome had a slightly longer control region. The availability of mitogenome sequence of whale shark will aid studies of molecular systematics, biogeography, genetic differentiation, and conservation genetics in this species., (Copyright © 2014 Elsevier B.V. All rights reserved.)

In order to better characterize the Bacillus anthracis typing phage AP50c, we designed a genetic screen to identify its bacterial receptor. Insertions of the transposon mariner or targeted deletions of the structural gene for the S-layer protein Sap and the sporulation genes spo0A, spo0B, and spo0F in B. anthracis Sterne resulted in phage resistance with concomitant defects in phage adsorption and infectivity. Electron microscopy of bacteria incubated with AP50c revealed phage particles associated with the surface of bacilli of the Sterne strain but not with the surfaces of Δsap, Δspo0A, Δspo0B, or Δspo0F mutants. The amount of Sap in the S layer of each of the spo0 mutant strains was substantially reduced compared to that of the parent strain, and incubation of AP50c with purified recombinant Sap led to a substantial reduction in phage activity. Phylogenetic analysis based on whole-genome sequences of B. cereus sensu lato strains revealed several closely related B. cereus and B. thuringiensis strains that carry sap genes with very high similarities to the sap gene of B. anthracis. Complementation of the Δsap mutant in trans with the wild-type B. anthracis sap or the sap gene from either of two different B. cereus strains that are sensitive to AP50c infection restored phage sensitivity, and electron microscopy confirmed attachment of phage particles to the surface of each of the complemented strains. Based on these data, we postulate that Sap is involved in AP50c infectivity, most likely acting as the phage receptor, and that the spo0 genes may regulate synthesis of Sap and/or formation of the S layer.

Rickettsia prowazekii is a notable intracellular pathogen, the agent of epidemic typhus, and a potential biothreat agent. We present here whole-genome sequence data for four strains of R. prowazekii, including one from a flying squirrel.

Unlabelled: Aeromonas hydrophila has increasingly been implicated as a virulent and antibiotic-resistant etiologic agent in various human diseases. In a previously published case report, we described a subject with a polymicrobial wound infection that included a persistent and aggressive strain of A. hydrophila (E1), as well as a more antibiotic-resistant strain of A. hydrophila (E2). To better understand the differences between pathogenic and environmental strains of A. hydrophila, we conducted comparative genomic and functional analyses of virulence-associated genes of these two wound isolates (E1 and E2), the environmental type strain A. hydrophila ATCC 7966(T), and four other isolates belonging to A. aquariorum, A. veronii, A. salmonicida, and A. caviae. Full-genome sequencing of strains E1 and E2 revealed extensive differences between the two and strain ATCC 7966(T). The more persistent wound infection strain, E1, harbored coding sequences for a cytotoxic enterotoxin (Act), a type 3 secretion system (T3SS), flagella, hemolysins, and a homolog of exotoxin A found in Pseudomonas aeruginosa. Corresponding phenotypic analyses with A. hydrophila ATCC 7966(T) and SSU as reference strains demonstrated the functionality of these virulence genes, with strain E1 displaying enhanced swimming and swarming motility, lateral flagella on electron microscopy, the presence of T3SS effector AexU, and enhanced lethality in a mouse model of Aeromonas infection. By combining sequence-based analysis and functional assays, we characterized an A. hydrophila pathotype, exemplified by strain E1, that exhibited increased virulence in a mouse model of infection, likely because of encapsulation, enhanced motility, toxin secretion, and cellular toxicity., Importance: Aeromonas hydrophila is a common aquatic bacterium that has increasingly been implicated in serious human infections. While many determinants of virulence have been identified in Aeromonas, rapid identification of pathogenic versus nonpathogenic strains remains a challenge for this genus, as it is for other opportunistic pathogens. This paper demonstrates, by using whole-genome sequencing of clinical Aeromonas strains, followed by corresponding virulence assays, that comparative genomics can be used to identify a virulent subtype of A. hydrophila that is aggressive during human infection and more lethal in a mouse model of infection. This aggressive pathotype contained genes for toxin production, toxin secretion, and bacterial motility that likely enabled its pathogenicity. Our results highlight the potential of whole-genome sequencing to transform microbial diagnostics; with further advances in rapid sequencing and annotation, genomic analysis will be able to provide timely information on the identities and virulence potential of clinically isolated microorganisms.

Chlamydia psittaci is an obligate intracellular bacterium. Interest in Chlamydia stems from its high degree of virulence as an intestinal and pulmonary pathogen across a broad range of animals, including humans. C. psittaci human pulmonary infections, referred to as psittacosis, can be life-threatening, which is why the organism was developed as a bioweapon in the 20th century and is listed as a CDC biothreat agent. One remarkable recent result from comparative genomics is the finding of frequent homologous recombination across the genome of the sexually transmitted and trachoma pathogen Chlamydia trachomatis. We sought to determine if similar evolutionary dynamics occurred in C. psittaci. We analyzed 20 C. psittaci genomes from diverse strains representing the nine known serotypes of the organism as well as infections in a range of birds and mammals, including humans. Genome annotation revealed a core genome in all strains of 911 genes. Our analyses showed that C. psittaci has a history of frequently switching hosts and undergoing recombination more often than C. trachomatis. Evolutionary history reconstructions showed genome-wide homologous recombination and evidence of whole-plasmid exchange. Tracking the origins of recombinant segments revealed that some strains have imported DNA from as-yet-unsampled or -unsequenced C. psittaci lineages or other Chlamydiaceae species. Three ancestral populations of C. psittaci were predicted, explaining the current population structure. Molecular clock analysis found that certain strains are part of a clonal epidemic expansion likely introduced into North America by South American bird traders, suggesting that psittacosis is a recently emerged disease originating in New World parrots.

The large number of sexually transmitted diseases and ocular trachoma cases that are caused globally each year by Chlamydia trachomatis has made this organism a World Health Organization priority for vaccine development. However, there is no gene transfer system for Chlamydia to help identify potential vaccine targets. To accelerate discoveries toward this goal, here we analyzed the broadest diversity of C. trachomatis genomes to date, including 25 geographically dispersed clinical and seven reference strains representing 14 of the 19 known serotypes. Strikingly, all 32 genomes were found to have evidence of DNA acquisition by homologous recombination in their history. Four distinct clades were identified, which correspond to all C. trachomatis disease phenotypes: lymphogranuloma venereum (LGV; Clade 1); noninvasive urogenital infections (Clade 2); ocular trachoma (Clade 3); and protocolitis (Clade 4; also includes some noninvasive urogenital infections). Although the ancestral relationship between clades varied, most strains acted as donor and recipient of recombination with no evidence for barriers to genetic exchange. The niche-specific LGV and trachoma clades have undergone less recombination, although the opportunity for mixing with strains from other clades that infect the rectal and ocular mucosa, respectively, is evident. Furthermore, there are numerous occasions for gene conversion events through sequential infections at the same anatomic sites. The size of recombinant segments is relatively small (~357 bp) compared with in vitro experiments of various C. trachomatis strains but is consistent with in vitro estimates for other bacterial species including Escherichia coli and Helicobacter pylori. Selection has also played a crucial role during the diversification of the organism. Clade 2 had the lowest nonsynonymous to synonymous ratio (dN/dS) but the highest effect of recombination, which is consistent with the widespread occurrence of synonymous substitutions in recombined genomic segments. The trachoma Clade 3 had the highest dN/dS estimates, which may be caused by an increased effect of genetic drift from niche specialization and a reduced effective population size. The degree of drift, selection, and recombination in C. trachomatis suggests that the challenge will remain to identify genomic regions that are stable and cross protective for the development of an efficacious vaccine.

Background: Spontaneous Bacillus anthracis mutants resistant to infection by phage AP50c (AP50R) exhibit a mucoid colony phenotype and secrete an extracellular matrix., Methods: Here we utilized a Roche/454-based whole genome sequencing approach to identify mutations that are candidates for conferring AP50c phage resistance, followed by genetic deletion and complementation studies to validate the whole genome sequence data and demonstrate that the implicated gene is necessary for AP50c phage infection., Results: Using whole genome sequence data, we mapped the relevant mutations in six AP50R strains to csaB. Eleven additional spontaneous mutants, isolated in two different genetic backgrounds, were screened by PCR followed by Sanger sequencing of the csaB gene. In each spontaneous mutant, we found either a non-synonymous substitution, a nonsense mutation, or a frame-shift mutation caused by single nucleotide polymorphisms or a 5 base pair insertion in csaB. All together, 5 and 12 of the 17 spontaneous mutations are predicted to yield altered full length and truncated CsaB proteins respectively. As expected from these results, a targeted deletion or frame-shift mutations introduced into csaB in a different genetic background, in a strain not exposed to AP50c, resulted in a phage resistant phenotype. Also, substitution of a highly conserved histidine residue with an alanine residue (H270A) in CsaB resulted in phage resistance, suggesting that a functional CsaB is necessary for phage sensitivity. Conversely, introduction of the wild type allele of csaB in cis into the csaB deletion mutant by homologous recombination or supplying the wild type CsaB protein in trans from a plasmid restored phage sensitivity. The csaB mutants accumulated cell wall material and appeared to have a defective S-layer, whereas these phenotypes were reverted in the complemented strains., Conclusions: Taken together, these data suggest an essential role for csaB in AP50c phage infection, most likely in phage adsorption. (The whole genome sequences generated from this study have been submitted to GenBank under SRA project ID: SRA023659.1 and sample IDs: AP50 R1: SRS113675.1, AP50 R2: SRS113676.1, AP50 R3: SRS113728.1, AP50 R4: SRS113733.1, AP50 R6: SRS113734.1, JB220 Parent: SRS150209.1, JB220 Mutant: SRS150211.1).

Methicillin-resistant Staphylococcus aureus (MRSA) strains with reduced susceptibility to vancomycin (MIC of 4 to 8 μg/ml) are referred to as vancomycin-intermediate S. aureus (VISA). In this study, we characterized two isogenic USA300 S. aureus isolates collected sequentially from a single patient with endocarditis where the S. aureus isolate changed from being susceptible to vancomycin (VSSA) (1 μg/ml) to VISA (8 μg/ml). In addition, the VISA isolate lost beta-lactamase activity and showed increased resistance to daptomycin and linezolid. The two strains did not differ in growth rate, but the VISA isolate had a thickened cell wall and was less autolytic. Transcriptome sequencing (RNA-seq) analysis comparing the two isolates grown to late exponential phase showed significant differences in transcription of cell surface protein genes (spa, SBI [second immunoglobulin-binding protein of S. aureus], and fibrinogen-binding proteins), regulatory genes (agrBCA, RNAIII, sarT, and saeRS), and others. Using whole-genome shotgun resequencing, we identified 6 insertion/deletion mutations between the VSSA and VISA isolates. A protein phosphatase 2C (PP2C) family phosphatase had a 6-bp (nonframeshift) insertion mutation in a highly conserved metal binding domain. Complementation of the clinical VISA isolate with a wild-type copy of the PP2C gene reduced the vancomycin and daptomycin MICs and increased autolytic activity, suggesting that this gene contributed to the reduced vancomycin susceptibility phenotype acquired in vivo. Creation of de novo mutants from the VSSA strain resulted in different mutations, demonstrating that reduced susceptibility to vancomycin in USA300 strains can occur via multiple routes, highlighting the complex nature of the VISA phenotype.

The key genes required for Bacillus anthracis to cause anthrax have been acquired recently by horizontal gene transfer. To understand the genetic background for the evolution of B. anthracis virulence, we obtained high-redundancy genome sequences of 45 strains of the Bacillus cereus sensu lato (s.l.) species that were chosen for their genetic diversity within the species based on the existing multilocus sequence typing scheme. From the resulting data, we called more than 324,000 new genes representing more than 12,333 new gene families for this group. The core genome size for the B. cereus s.l. group was ∼1750 genes, with another 2150 genes found in almost every genome constituting the extended core. There was a paucity of genes specific and conserved in any clade. We found no evidence of recent large-scale gene loss in B. anthracis or for unusual accumulation of nonsynonymous DNA substitutions in the chromosome; however, several B. cereus genomes isolated from soil and not previously associated with human disease were degraded to various degrees. Although B. anthracis has undergone an ecological shift within the species, its chromosome does not appear to be exceptional on a macroscopic scale compared with close relatives.

Background: Recent years have shown a marked increase in the use of next-generation sequencing technologies for quantification of gene expression (RNA sequencing, RNA-Seq). The expression level of a gene is a function of both its rate of transcription and RNA decay, and the influence of mRNA decay rates on gene expression in genome-wide studies of Gram-positive bacteria is under-investigated., Results: In this work, we employed RNA-Seq in a genome-wide determination of mRNA half-lives in the Gram-positive bacterium Bacillus cereus. By utilizing a newly developed normalization protocol, RNA-Seq was used successfully to determine global mRNA decay rates at the single nucleotide level. The analysis revealed positional degradation patterns, with mRNAs being degraded from both ends of the molecule, indicating that both 5' to 3' and 3' to 5' directions of RNA decay are present in B. cereus. Other operons showed segmental degradation patterns where specific ORFs within polycistrons were degraded at variable rates, underlining the importance of RNA processing in gene regulation. We determined the half-lives for more than 2,700 ORFs in B. cereus ATCC 10987, ranging from less than one minute to more than fifteen minutes, and showed that mRNA decay rate correlates globally with mRNA expression level, GC content, and functional class of the ORF., Conclusions: To our knowledge, this study presents the first global analysis of mRNA decay in a bacterium at single nucleotide resolution. We provide a proof of principle for using RNA-Seq in bacterial mRNA decay analysis, revealing RNA processing patterns at the single nucleotide level.

Background: Although genome-wide transcriptional analysis has been used for many years to study bacterial gene expression, many aspects of the bacterial transcriptome remain undefined. One example is antisense transcription, which has been observed in a number of bacteria, though the function of antisense transcripts, and their distribution across the bacterial genome, is still unclear., Methodology/principal Findings: Single-stranded RNA-seq results revealed a widespread and non-random pattern of antisense transcription covering more than two thirds of the B. anthracis genome. Our analysis revealed a variety of antisense structural patterns, suggesting multiple mechanisms of antisense transcription. The data revealed several instances of sense and antisense expression changes in different growth conditions, suggesting that antisense transcription may play a role in the ways in which B. anthracis responds to its environment. Significantly, genome-wide antisense expression occurred at consistently higher levels on the lagging strand, while the leading strand showed very little antisense activity. Intrasample gene expression comparisons revealed a gene dosage effect in all growth conditions, where genes farthest from the origin showed the lowest overall range of expression for both sense and antisense directed transcription. Additionally, transcription from both strands was verified using a novel strand-specific assay. The variety of structural patterns we observed in antisense transcription suggests multiple mechanisms for this phenomenon, suggesting that some antisense transcription may play a role in regulating the expression of key genes, while some may be due to chromosome replication dynamics and transcriptional noise., Conclusions/significance: Although the variety of structural patterns we observed in antisense transcription suggest multiple mechanisms for antisense expression, our data also clearly indicate that antisense transcription may play a genome-wide role in regulating the expression of key genes in Bacillus species. This study illustrates the surprising complexity of prokaryotic RNA abundance for both strands of a bacterial chromosome.

Background: Chlamydia trachomatis (Ct) is the most common cause of bacterial sexually transmitted diseases (STD) worldwide. While commercial nucleic acid amplification tests (NAAT) are available for Ct, none are rapid or inexpensive enough to be used at the point-of-care (POC). Towards the first Ct POC NAAT, we developed a microfluidic assay that simultaneously interrogates nine Ct loci in 20 minutes., Methodology and Principal Findings: Endocervical samples were selected from 263 women at high risk for Ct STDs (∼35% prevalence). A head-to-head comparison was performed with the Roche-Amplicor NAAT. 129 (49.0%) and 88 (33.5%) samples were positive by multiplex and Amplicor assays, respectively. Sequencing resolved 71 discrepant samples, confirming 53 of 53 positive multiplex samples and 12 of 18 positive Amplicor samples. The sensitivity and specificity were 91.5% and 100%, and 62.4% and 95.9%, respectively, for multiplex and Amplicor assays. Positive and negative predictive values were 100% and 91%, and 94.1% and 68.6%, respectively., Conclusions: This is the first rapid multiplex approach to Ct detection, and the assay was also found to be superior to a commercial NAAT. In effect, nine simultaneous reactions significantly increased sensitivity and specificity. Our assay can potentially increase Ct detection in globally diverse clinical settings at the POC.

Chlamydiae are evolutionarily well-separated bacteria that live exclusively within eukaryotic host cells. They include important human pathogens such as Chlamydia trachomatis as well as symbionts of protozoa. As these bacteria are experimentally challenging and genetically intractable, our knowledge about them is still limited. In this study, we obtained the genome sequences of Simkania negevensis Z, Waddlia chondrophila 2032/99, and Parachlamydia acanthamoebae UV-7. This enabled us to perform the first comprehensive comparative and phylogenomic analysis of representative members of four major families of the Chlamydiae, including the Chlamydiaceae. We identified a surprisingly large core gene set present in all genomes and a high number of diverse accessory genes in those Chlamydiae that do not primarily infect humans or animals, including a chemosensory system in P. acanthamoebae and a type IV secretion system. In S. negevensis, the type IV secretion system is encoded on a large conjugative plasmid (pSn, 132 kb). Phylogenetic analyses suggested that a plasmid similar to the S. negevensis plasmid was originally acquired by the last common ancestor of all four families and that it was subsequently reduced, integrated into the chromosome, or lost during diversification, ultimately giving rise to the extant virulence-associated plasmid of pathogenic chlamydiae. Other virulence factors, including a type III secretion system, are conserved among the Chlamydiae to variable degrees and together with differences in the composition of the cell wall reflect adaptation to different host cells including convergent evolution among the four chlamydial families. Phylogenomic analysis focusing on chlamydial proteins with homology to plant proteins provided evidence for the acquisition of 53 chlamydial genes by a plant progenitor, lending further support for the hypothesis of an early interaction between a chlamydial ancestor and the primary photosynthetic eukaryote.

Chlamydia psittaci is a highly prevalent avian pathogen and the cause of a potentially lethal zoonosis, causing life-threatening pneumonia in humans. We report the genome sequences of C. psittaci 6BC, the prototype strain of the species, and C. psittaci Cal10, a widely used laboratory strain.

The following series of concise summaries addresses the evolution of infectious agents in relation to sex in animals and humans from the perspective of three specific questions: (1) what have we learned about the likely origin and phylogeny, up to the establishment of the infectious agent in the genital econiche, including the relative frequency of its sexual transmission; (2) what further research is needed to provide additional knowledge on some of these evolutionary aspects; and (3) what evolutionary considerations might aid in providing novel approaches to the more practical clinical and public health issues facing us currently and in the future?, (© 2011 New York Academy of Sciences.)

The endospore-forming Gram-positive pathogen Bacillus anthracis is responsible for the usually fatal disease, inhalational anthrax. The success of this pathogen is dependent on its ability to subvert elements of the innate immune system of its animal hosts. B. anthracis spores, which are the main infective agent, are engulfed and germinate in patrolling alveolar macrophages. In order for the infection to progress, the resulting vegetative cells must resist the antimicrobial oxidative burst mounted by the host NADPH oxidase complex. The response of B. anthracis to this and other macrophage-related stresses is therefore of major importance to the success of this pathogen, and consequently we have analysed the superoxide and peroxide stress stimulons of B. anthracis strain UM23C1-2 by means of a combined transcriptomics and proteomics approach. The results show distinct patterns of expression in response to paraquat (endogenous superoxide) and hydrogen peroxide stress. While the main response to paraquat is the induction of iron uptake pathways, the response to peroxide predominantly involves the induction of protection and repair mechanisms. Comparisons between the responses of B. anthracis and related soil bacterium, B. subtilis, reveal differences that are likely to be relevant to their respective habitats., (Copyright © 2011 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.)

Chlamydia pecorum is an obligate intracellular bacterial pathogen that causes diverse disease in a wide variety of economically important mammals. We report the finished complete genome sequence of C. pecorum E58, the type strain for the species.

Background: Chlamydia trachomatis is an obligate intracellular bacterial parasite, which causes several severe and debilitating diseases in humans. This study uses comparative genomic analyses of 12 complete published C. trachomatis genomes to assess the contribution of recombination and selection in this pathogen and to understand the major evolutionary forces acting on the genome of this bacterium., Results: The conserved core genes of C. trachomatis are a large proportion of the pan-genome: we identified 836 core genes in C. trachomatis out of a range of 874-927 total genes in each genome. The ratio of recombination events compared to mutation (ρ/θ) was 0.07 based on ancestral reconstructions using the ClonalFrame tool, but recombination had a significant effect on genetic diversification (r/m=0.71). The distance-dependent decay of linkage disequilibrium also indicated that C. trachomatis populations behaved intermediately between sexual and clonal extremes. Fifty-five genes were identified as having a history of recombination and 92 were under positive selection based on statistical tests. Twenty-three genes showed evidence of being under both positive selection and recombination, which included genes with a known role in virulence and pathogencity (e.g., ompA, pmps, tarp). Analysis of inter-clade recombination flux indicated non-uniform currents of recombination between clades, which suggests the possibility of spatial population structure in C. trachomatis infections., Conclusions: C. trachomatis is the archetype of a bacterial species where recombination is relatively frequent yet gene gains by horizontal gene transfer (HGT) and losses (by deletion) are rare. Gene conversion occurs at sites across the whole C. trachomatis genome but may be more often fixed in genes that are under diversifying selection. Furthermore, genome sequencing will reveal patterns of serotype specific gene exchange and selection that will generate important research questions for understanding C. trachomatis pathogenesis., Reviewers: This article was reviewed by Dr. Jeremy Selengut, Dr. Lee S. Katz (nominated by Dr. I. King Jordan) and Dr. Arcady Mushegian., (© 2011 Joseph et al; licensee BioMed Central Ltd.)

Unlabelled: Chlamydia trachomatis is an obligate intracellular bacterium that causes a diversity of severe and debilitating diseases worldwide. Sporadic and ongoing outbreaks of lymphogranuloma venereum (LGV) strains among men who have sex with men (MSM) support the need for research on virulence factors associated with these organisms. Previous analyses have been limited to single genes or genomes of laboratory-adapted reference strain L(2)/434 and outbreak strain L(2)b/UCH-1/proctitis. We characterized an unusual LGV strain, termed L(2)c, isolated from an MSM with severe hemorrhagic proctitis. L(2)c developed nonfusing, grape-like inclusions and a cytotoxic phenotype in culture, unlike the LGV strains described to date. Deep genome sequencing revealed that L(2)c was a recombinant of L(2) and D strains with conserved clustered regions of genetic exchange, including a 78-kb region and a partial, yet functional, toxin gene that was lost with prolonged culture. Indels (insertions/deletions) were discovered in an ftsK gene promoter and in the tarp and hctB genes, which encode key proteins involved in replication, inclusion formation, and histone H1-like protein activity, respectively. Analyses suggest that these indels affect gene and/or protein function, supporting the in vitro and disease phenotypes. While recombination has been known to occur for C. trachomatis based on gene sequence analyses, we provide the first whole-genome evidence for recombination between a virulent, invasive LGV strain and a noninvasive common urogenital strain. Given the lack of a genetic system for producing stable C. trachomatis mutants, identifying naturally occurring recombinants can clarify gene function and provide opportunities for discovering avenues for genomic manipulation., Importance: Lymphogranuloma venereum (LGV) is a prevalent and debilitating sexually transmitted disease in developing countries, although there are significant ongoing outbreaks in Australia, Europe, and the United States among men who have sex with men (MSM). Relatively little is known about LGV virulence factors, and only two LGV genomes have been sequenced to date. We isolated an LGV strain from an MSM with severe hemorrhagic proctitis that was morphologically unique in tissue culture compared with other LGV strains. Bioinformatic and statistical analyses identified the strain as a recombinant of L(2) and D strains with highly conserved clustered regions of genetic exchange. The unique culture morphology and, more importantly, disease phenotype could be traced to the genes involved in recombination. The findings have implications for bacterial species evolution and, in the case of ongoing LGV outbreaks, suggest that recombination is a mechanism for strain emergence that results in significant disease pathology.

Background: OmniLog™ phenotype microarrays (PMs) have the capability to measure and compare the growth responses of biological samples upon exposure to hundreds of growth conditions such as different metabolites and antibiotics over a time course of hours to days. In order to manage the large amount of data produced from the OmniLog™ instrument, PheMaDB (Phenotype Microarray DataBase), a web-based relational database, was designed. PheMaDB enables efficient storage, retrieval and rapid analysis of the OmniLog™ PM data., Description: PheMaDB allows the user to quickly identify records of interest for data analysis by filtering with a hierarchical ordering of Project, Strain, Phenotype, Replicate, and Temperature. PheMaDB then provides various statistical analysis options to identify specific growth pattern characteristics of the experimental strains, such as: outlier analysis, negative controls analysis (signal/background calibration), bar plots, pearson's correlation matrix, growth curve profile search, k-means clustering, and a heat map plot. This web-based database management system allows for both easy data sharing among multiple users and robust tools to phenotype organisms of interest., Conclusions: PheMaDB is an open source system standardized for OmniLog™ PM data. PheMaDB could facilitate the banking and sharing of phenotype data. The source code is available for download at http://phemadb.sourceforge.net.

Background: Despite the decades-long use of Bacillus atrophaeus var. globigii (BG) as a simulant for biological warfare (BW) agents, knowledge of its genome composition is limited. Furthermore, the ability to differentiate signatures of deliberate adaptation and selection from natural variation is lacking for most bacterial agents. We characterized a lineage of BGwith a long history of use as a simulant for BW operations, focusing on classical bacteriological markers, metabolic profiling and whole-genome shotgun sequencing (WGS)., Results: Archival strains and two "present day" type strains were compared to simulant strains on different laboratory media. Several of the samples produced multiple colony morphotypes that differed from that of an archival isolate. To trace the microevolutionary history of these isolates, we obtained WGS data for several archival and present-day strains and morphotypes. Bacillus-wide phylogenetic analysis identified B. subtilis as the nearest neighbor to B. atrophaeus. The genome of B. atrophaeus is, on average, 86% identical to B. subtilis on the nucleotide level. WGS of variants revealed that several strains were mixed but highly related populations and uncovered a progressive accumulation of mutations among the "military" isolates. Metabolic profiling and microscopic examination of bacterial cultures revealed enhanced growth of "military" isolates on lactate-containing media, and showed that the "military" strains exhibited a hypersporulating phenotype., Conclusions: Our analysis revealed the genomic and phenotypic signatures of strain adaptation and deliberate selection for traits that were desirable in a simulant organism. Together, these results demonstrate the power of whole-genome and modern systems-level approaches to characterize microbial lineages to develop and validate forensic markers for strain discrimination and reveal signatures of deliberate adaptation.

Before the anthrax letter attacks of 2001, the developing field of microbial forensics relied on microbial genotyping schemes based on a small portion of a genome sequence. Amerithrax, the investigation into the anthrax letter attacks, applied high-resolution whole-genome sequencing and comparative genomics to identify key genetic features of the letters' Bacillus anthracis Ames strain. During systematic microbiological analysis of the spore material from the letters, we identified a number of morphological variants based on phenotypic characteristics and the ability to sporulate. The genomes of these morphological variants were sequenced and compared with that of the B. anthracis Ames ancestor, the progenitor of all B. anthracis Ames strains. Through comparative genomics, we identified four distinct loci with verifiable genetic mutations. Three of the four mutations could be directly linked to sporulation pathways in B. anthracis and more specifically to the regulation of the phosphorylation state of Spo0F, a key regulatory protein in the initiation of the sporulation cascade, thus linking phenotype to genotype. None of these variant genotypes were identified in single-colony environmental B. anthracis Ames isolates associated with the investigation. These genotypes were identified only in B. anthracis morphotypes isolated from the letters, indicating that the variants were not prevalent in the environment, not even the environments associated with the investigation. This study demonstrates the forensic value of systematic microbiological analysis combined with whole-genome sequencing and comparative genomics.

We performed whole-genome amplification followed by hybridization of custom-designed resequencing arrays to resequence 303 kb of genomic sequence from a worldwide panel of 39 Bacillus anthracis strains. We used an efficient algorithm contained within a custom software program, UniqueMER, to identify and mask repetitive sequences on the resequencing array to reduce false-positive identification of genetic variation, which can arise from cross-hybridization. We discovered a total of 240 single nucleotide variants (SNVs) and showed that B. anthracis strains have an average of 2.25 differences per 10,000 bases in the region we resequenced. Common SNVs in this region are found to be in complete linkage disequilibrium. These patterns of variation suggest there has been little if any historical recombination among B. anthracis strains since the origin of the pathogen. This pattern of common genetic variation suggests a framework for recognizing new or genetically engineered 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.

New sequencing technologies have made the production of bacterial genome sequences increasingly easy, and it can be confidently forecasted that vast genomic databases will be generated in the next few years. Here, we detail how collections of bacterial genomes from a particular species (population genomics libraries) have already been used to improve the design of several diagnostic assays for bacterial pathogens. Genome sequencing itself is also becoming more commonly used for epidemiological, forensic and clinical investigations. There is an opportunity for the further development of bioinformatic tools to bring even further value to bacterial diagnostic genomics., (Copyright © 2010 Elsevier Ltd. All rights reserved.)

Background: Despite the global threat caused by arthropod-borne viruses, there is not an efficient method for screening vector populations to detect novel viral sequences. Current viral detection and surveillance methods based on culture can be costly and time consuming and are predicated on prior knowledge of the etiologic agent, as they rely on specific oligonucleotide primers or antibodies. Therefore, these techniques may be unsuitable for situations when the causative agent of an outbreak is unknown., Methodology/principal Findings: In this study we explored the use of high-throughput pyrosequencing for surveillance of arthropod-borne RNA viruses. Dengue virus, a member of the positive strand RNA Flavivirus family that is transmitted by several members of the Aedes genus of mosquitoes, was used as a model. Aedes aegypti mosquitoes experimentally infected with dengue virus type 1 (DENV-1) were pooled with noninfected mosquitoes to simulate samples derived from ongoing arbovirus surveillance programs. Using random-primed methods, total RNA was reverse-transcribed and resulting cDNA subjected to 454 pyrosequencing., Conclusions/significance: In two types of samples, one with 5 adult mosquitoes infected with DENV-1- and the other with 1 DENV-1 infected mosquito and 4 noninfected mosquitoes, we identified DENV-1 DNA sequences. DENV-1 sequences were not detected in an uninfected control pool of 5 adult mosquitoes. We calculated the proportion of the Ae. aegypti metagenome contributed by each infecting Dengue virus genome (p(IP)), which ranged from 2.75×10(-8) to 1.08×10(-7). DENV-1 RNA was sufficiently concentrated in the mosquito that its detection was feasible using current high-throughput sequencing instrumentation. We also identified some of the components of the mosquito microflora on the basis of the sequence of expressed RNA. This included members of the bacterial genera Pirellula and Asaia, various fungi, and a potentially uncharacterized mycovirus.

Background: The anthrax letter attacks of 2001 highlighted the need for rapid identification of biothreat agents not only for epidemiological surveillance of the intentional outbreak but also for implementing appropriate countermeasures, such as antibiotic treatment, in a timely manner to prevent further casualties. It is clear from the 2001 cases that survival may be markedly improved by administration of antimicrobial therapy during the early symptomatic phase of the illness; i.e., within 3 days of appearance of symptoms. Microbiological detection methods are feasible only for organisms that can be cultured in vitro and cannot detect all genetic modifications with the exception of antibiotic resistance. Currently available immuno or nucleic acid-based rapid detection assays utilize known, organism-specific proteins or genomic DNA signatures respectively. Hence, these assays lack the ability to detect novel natural variations or intentional genetic modifications that circumvent the targets of the detection assays or in the case of a biological attack using an antibiotic resistant or virulence enhanced Bacillus anthracis, to advise on therapeutic treatments., Methodology/principal Findings: We show here that the Roche 454-based pyrosequencing can generate whole genome draft sequences of deep and broad enough coverage of a bacterial genome in less than 24 hours. Furthermore, using the unfinished draft sequences, we demonstrate that unbiased identification of known as well as heretofore-unreported genetic modifications that include indels and single nucleotide polymorphisms conferring antibiotic and phage resistances is feasible within the next 12 hours., Conclusions/significance: Second generation sequencing technologies have paved the way for sequence-based rapid identification of both known and previously undocumented genetic modifications in cultured, conventional and newly emerging biothreat agents. Our findings have significant implications in the context of whole genome sequencing-based routine clinical diagnostics as well as epidemiological surveillance of natural disease outbreaks caused by bacterial and viral agents.

Background: Multiple locus sequence typing (MLST) has become a central genotyping strategy for analysis of bacterial populations. The scheme involves de novo sequencing of 6-8 housekeeping loci to assign unique sequence types. In this work we adapted MLST to a rapid microfluidics platform in order to enhance speed and reduce laboratory labor time., Methodology/principal Findings: Using two integrated microfluidic devices, DNA was purified from 100 Bacillus cereus soil isolates, used as a template for multiplex amplification of 7 loci and sequenced on forward and reverse strands. The time on instrument from loading genomic DNA to generation of electropherograms was only 1.5 hours. We obtained full-length sequence of all seven MLST alleles from 84 representing 46 different Sequence Types. At least one allele could be sequenced from a further 15 strains. The nucleotide diversity of B. cereus isolated in this study from one location in Rockville, Maryland (0.04 substitutions per site) was found to be as great as the global collection of isolates., Conclusions/significance: Biogeographical investigation of pathogens is only one of a panoply of possible applications of microfluidics based MLST; others include microbiologic forensics, biothreat identification, and rapid characterization of human clinical samples.

While new sequencing technologies have ushered in an era where microbial genomes can be easily sequenced, the goal of routinely producing high-quality draft and finished genomes in a cost-effective fashion has still remained elusive. Due to shorter read lengths and limitations in library construction protocols, shotgun sequencing and assembly based on these technologies often results in fragmented assemblies. Correspondingly, while draft assemblies can be obtained in days, finishing can take many months and hence the time and effort can only be justified for high-priority genomes and in large sequencing centers. In this work, we revisit this issue in light of our own experience in producing finished and nearly-finished genomes for a range of microbial species in a small-lab setting. These genomes were finished with surprisingly little investments in terms of time, computational effort and lab work, suggesting that the increased access to sequencing might also eventually lead to a greater proportion of finished genomes from small labs and genomics cores.

Background: New DNA sequencing technologies have enabled detailed comparative genomic analyses of entire genera of bacterial pathogens. Prior to this study, three species of the enterobacterial genus Yersinia that cause invasive human diseases (Yersinia pestis, Yersinia pseudotuberculosis, and Yersinia enterocolitica) had been sequenced. However, there were no genomic data on the Yersinia species with more limited virulence potential, frequently found in soil and water environments., Results: We used high-throughput sequencing-by-synthesis instruments to obtain 25- to 42-fold average redundancy, whole-genome shotgun data from the type strains of eight species: Y. aldovae, Y. bercovieri, Y. frederiksenii, Y. kristensenii, Y. intermedia, Y. mollaretii, Y. rohdei, and Y. ruckeri. The deepest branching species in the genus, Y. ruckeri, causative agent of red mouth disease in fish, has the smallest genome (3.7 Mb), although it shares the same core set of approximately 2,500 genes as the other members of the species, whose genomes range in size from 4.3 to 4.8 Mb. Yersinia genomes had a similar global partition of protein functions, as measured by the distribution of Cluster of Orthologous Groups families. Genome to genome variation in islands with genes encoding functions such as ureases, hydrogenases and B-12 cofactor metabolite reactions may reflect adaptations to colonizing specific host habitats., Conclusions: Rapid high-quality draft sequencing was used successfully to compare pathogenic and non-pathogenic members of the Yersinia genus. This work underscores the importance of the acquisition of horizontally transferred genes in the evolution of Y. pestis and points to virulence determinants that have been gained and lost on multiple occasions in the history of the genus.

We characterize CdSe/ZnS quantum dot (QD) binding to genetically modified bacteriophage as a model for bacterial detection. Interactions among QDs, lambda (lambda) phage, and Escherichia coli are examined by several cross-validated methods. Flow and image-based cytometry clarify fluorescent labeling of bacteria, with image-based cytometry additionally reporting the number of decorated phage bound to cells. Transmission electron microscopy, image-based cytometry, and electrospray differential mobility analysis allow quantization of QDs attached to each phage (4-17 QDs) and show that lambda phage used in this study exhibits enhanced QD binding to the capsid by nearly a factor of four compared to bacteriophage T7. Additionally, the characterization methodology presented can be applied to the quantitative characterization of other fluorescent nanocrystal-biological conjugates., (2009 Wiley Periodicals, Inc.)

Zoonotic infections are a growing threat to global health. Chlamydia pneumoniae is a major human pathogen that is widespread in human populations, causing acute respiratory disease, and has been associated with chronic disease. C. pneumoniae was first identified solely in human populations; however, its host range now includes other mammals, marsupials, amphibians, and reptiles. Australian koalas (Phascolarctos cinereus) are widely infected with two species of Chlamydia, C. pecorum and C. pneumoniae. Transmission of C. pneumoniae between animals and humans has not been reported; however, two other chlamydial species, C. psittaci and C. abortus, are known zoonotic pathogens. We have sequenced the 1,241,024-bp chromosome and a 7.5-kb cryptic chlamydial plasmid of the koala strain of C. pneumoniae (LPCoLN) using the whole-genome shotgun method. Comparative genomic analysis, including pseudogene and single-nucleotide polymorphism (SNP) distribution, and phylogenetic analysis of conserved genes and SNPs against the human isolates of C. pneumoniae show that the LPCoLN isolate is basal to human isolates. Thus, we propose based on compelling genomic and phylogenetic evidence that humans were originally infected zoonotically by an animal isolate(s) of C. pneumoniae which adapted to humans primarily through the processes of gene decay and plasmid loss, to the point where the animal reservoir is no longer required for transmission.

Chlamydia trachomatis is a global cause of blinding trachoma and sexually transmitted infections (STIs). We used comparative genomics of the family Chlamydiaceae to select conserved housekeeping genes for C. trachomatis multilocus sequencing, characterizing 19 reference and 68 clinical isolates from 6 continental/subcontinental regions. There were 44 sequence types (ST). Identical STs for STI isolates were recovered from different regions, whereas STs for trachoma isolates were restricted by continent. Twenty-nine of 52 alleles had nonuniform distributions of frequencies across regions (p<0.001). Phylogenetic analysis showed 3 disease clusters: invasive lymphogranuloma venereum strains, globally prevalent noninvasive STI strains (ompA genotypes D/Da, E, and F), and nonprevalent STI strains with a trachoma subcluster. Recombinant strains were observed among STI clusters. Single nucleotide polymorphisms (SNPs) were predictive of disease specificity. Multilocus and SNP typing can now be used to detect diverse and emerging C. trachomatis strains for epidemiologic and evolutionary studies of trachoma and STI populations worldwide.

We sought to identify proteins in the Bacillus anthracis spore, conserved in other strains of the closely related Bacillus cereus group, that elicit an immune response in mammals. Two high throughput approaches were used. First, an in silico screening identified 200 conserved putative B. anthracis spore components. A total of 192 of those candidate genes were expressed and purified in vitro, 75 of which reacted with the rabbit immune sera generated against B. anthracis spores. The second approach was to screen for cross-reacting antigens in the spore proteome of 10 diverse B. cereus group strains. Two-dimensional electrophoresis resolved more than 200 protein spots in each spore preparation. About 72% of the protein spots were found in all the strains. 18 of these conserved proteins reacted against anti-B. anthracis spore rabbit immune sera, two of which (alanine racemase, Dal-1 and the methionine transporter, MetN) overlapped the set of proteins identified using the in silico screen. A conserved repeat domain protein (Crd) was the most immunoreactive protein found broadly across B. cereus sensu lato strains. We have established an approach for finding conserved targets across a species using population genomics and proteomics. The results of these screens suggest the possibility of a multiepitope antigen for broad host range diagnostics or therapeutics against Bacillus spore infection.

Unlabelled: DIYA (Do-It-Yourself Annotator) is a modular and configurable open source pipeline software, written in Perl, used for the rapid annotation of bacterial genome sequences. The software is currently used to take DNA contigs as input, either in the form of complete genomes or the result of shotgun sequencing, and produce an annotated sequence in Genbank file format as output., Availability: Distribution and source code are available at (https://sourceforge.net/projects/diyg/)., Supplementary Information: Supplementary data are available at Bioinformatics online.

The pathogenic bacterium Bacillus anthracis has become the subject of intense study as a result of its use in a bioterrorism attack in the United States in September and October 2001. Previous studies suggested that B. anthracis Ames Ancestor, the original Ames fully virulent plasmid-containing isolate, was the ideal reference. This study describes the complete genome sequence of that original isolate, derived from a sample kept in cold storage since 1981.

The development of vaccines against Campylobacter jejuni would be facilitated by the ability to perform phase II challenge studies. However, molecular mimicry of the lipooligosaccharide (LOS) of most C. jejuni strains with human gangliosides presents safety concerns about the development of Guillain-Barré syndrome. Clinical isolates of C. jejuni that appeared to lack genes for the synthesis of ganglioside mimics were identified by DNA probe analyses. Two clinical isolates from Southeast Asia (strains BH-01-0142 and CG8421) were determined to express the LOS type containing N-acetyl quinovosamine. No ganglioside structures were observed to be present in the LOSs of these strains, and pyrosequence analyses of the genomes of both strains confirmed the absence of genes involved in ganglioside mimicry. The capsule polysaccharide (CPS) of BH-01-0142 was determined to be composed of galactose (Gal), 6-deoxy-ido-heptose, and, in smaller amounts, D-glycero-D-ido-heptose, and the CPS of CG8421 was observed to contain Gal, 6-deoxy-altro-heptose, N-acetyl-glucosamine, and minor amounts of 6-deoxy-3-O-Me-altro-heptose. Both CPSs were shown to carry O-methyl-phosphoramidate. The two genomes contained strain-specific zones, some of which could be traced to a plasmid origin, and both contained a large chromosomal insertion related to the CJEI3 element of C. jejuni RM1221. The genomes of both strains shared a high degree of similarity to each other and, with the exception of the capsule locus of CG8421, to the type strain of the HS3 serotype, TGH9011.

The genome sequence of a Bacillus anthracis-specific clear plaque mutant phage, AP50c, contains 31 open reading frames spanning 14,398 bp, has two mutations compared to wild-type AP50t, and has a colinear genome architecture highly similar to that of gram-positive Tectiviridae phages. Spontaneous AP50c-resistant B. anthracis mutants exhibit a mucoid colony phenotype.

Bacillus thuringiensis serovar Monterrey strain BGSC 4AJ1 produced a microscopically visible capsule that reacted with a fluorescent antibody specific for the poly-gamma-d-glutamic acid (PGA) capsule of Bacillus anthracis. PGA capsule biosynthesis genes with 75%, 81%, 72%, 65% and 63% similarity, respectively, to those of the B. anthracis capBCADE cluster were present on a plasmid (pAJ1-1). Strain BGSC 4AJ1, together with five strains of Bacillus cereus that hybridized to a PGA cap gene probe, were analyzed phylogenetically using six housekeeping genes of a B. cereus multilocus sequence typing scheme. Bacillus thuringiensis BGSC 4AJ1 shared four identical alleles with B. anthracis and was the second most closely related to this bacterium of the 674 isolates in the multilocus sequence typing database. The other cap+ strains were distributed among various lineages of Clade 1 of the B. cereus group.

The ability to distinguish microbial pathogens from closely related but nonpathogenic strains is key to understanding the population biology of these organisms. In this regard, Bacillus anthracis, the bacterium that causes inhalational anthrax, is of interest because it is closely related and often difficult to distinguish from other members of the B. cereus group that can cause diverse diseases. We employed custom-designed resequencing arrays (RAs) based on the genome sequence of Bacillus anthracis to generate 422 kb of genomic sequence from a panel of 41 Bacillus cereus sensu lato strains. Here we show that RAs represent a "one reaction" genotyping technology with the ability to discriminate between highly similar B. anthracis isolates and more divergent strains of the B. cereus s.l. Clade 1. Our data show that RAs can be an efficient genotyping technology for pre-screening the genetic diversity of large strain collections to selected the best candidates for whole genome sequencing.

Motivation: New, high-throughput sequencing technologies have made it feasible to cheaply generate vast amounts of sequence information from a genome of interest. The computational reconstruction of the complete sequence of a genome is complicated by specific features of these new sequencing technologies, such as the short length of the sequencing reads and absence of mate-pair information. In this article we propose methods to overcome such limitations by incorporating information from optical restriction maps., Results: We demonstrate the robustness of our methods to sequencing and assembly errors using extensive experiments on simulated datasets. We then present the results obtained by applying our algorithms to data generated from two bacterial genomes Yersinia aldovae and Yersinia kristensenii. The resulting assemblies contain a single scaffold covering a large fraction of the respective genomes, suggesting that the careful use of optical maps can provide a cost-effective framework for the assembly of genomes., Availability: The tools described here are available as an open-source package at ftp://ftp.cbcb.umd.edu/pub/software/soma

In Chlamydiaceae, the nucleotide sequence between the 5S rRNA gene and the gene for subunit F of the Na(+)-translocating NADH-quinone reductase (nqrF or dmpP) has varied lengths and gene contents. We analyzed this site in 45 Chlamydiaceae strains having diverse geographical and pathological origins and including members of all nine species.

The genetic relatedness of the Bacillus anthracis typing phages Gamma and Cherry was determined by nucleotide sequencing and comparative analysis. The genomes of these two phages were identical except at three variable loci, which showed heterogeneity within individual lysates and among Cherry, Wbeta, Fah, and four Gamma bacteriophage sequences.

Background: Bacillus anthracis is considered to be a recently emerged clone within the Bacillus cereus sensu lato group. The B. anthracis genome sequence contains four putative lambdoid prophages. We undertook this study in order to understand whether the four prophages are unique to B. anthracis and whether they produce active phages., Results: More than 300 geographically and temporally divergent isolates of B. anthracis and its near neighbors were screened by PCR for the presence of specific DNA sequences from each prophage region. Every isolate of B. anthracis screened by PCR was found to produce all four phage-specific amplicons whereas none of the non-B. anthracis isolates, produced more than one phage-specific amplicon. Excision of prophages could be detected by a PCR based assay for attP sites on extra-chromosomal phage circles and for attB sites on phage-excised chromosomes. SYBR-green real-time PCR assays indicated that prophage excision occurs at very low frequencies (2 x 10(-5) - 8 x 10(-8)/cell). Induction with mitomycin C increased the frequency of excision of one of the prophages by approximately 250 fold. All four prophages appear to be defective since, mitomycin C induced culture did not release any viable phage particle or lyse the cells or reveal any phage particle under electron microscopic examination., Conclusion: The retention of all four putative prophage regions across all tested strains of B. anthracis is further evidence of the very recent emergence of this lineage and the prophage regions may be useful for differentiating the B. anthracis chromosome from that of its neighbors. All four prophages can excise at low frequencies, but are apparently defective in phage production.

Horizontal gene transfer (HGT) has a far more significant role than gene duplication in bacterial evolution. This has recently been illustrated by work demonstrating the importance of HGT in the emergence of bacterial metabolic networks, with horizontally acquired genes being placed in peripheral pathways at the outer branches of the networks.