The complete genome sequence of Bacillus thuringiensis Al Hakam

Bt GS57 Interaction With Gut Microbiota Accelerates Spodoptera exigua Mortality

The Beet armyworm Spodoptera exigua (Lepidoptera: Noctuidae, Spodoptera) is an important global polyphagous pest. Pathogen infection could destroy the intestinal microbial homeostasis of insects, leading to the death of the host. However, the effect of the host intestinal microbial community on the insecticidal effect of Bacillus thuringiensis is rarely studied. In this study, the genome characteristics of Bt GS57 and the diversity and functions of the gut bacteria in S. exigua are investigated using crystal morphology, biological activity, and Illumina HiSeq high-throughput sequencing. The total size of the Bt GS57 genome is 6.17 Mbp with an average G + C content of 35.66%. Furthermore, the Bt GS57 genome contains six cry genes: cry1Ca, cry1Da, cry2Ab, cry9Ea, cry1Ia, and cry1Aa, and a vegetative insecticidal protein gene vip3Aa. The Bt GS57 strain can produce biconical crystals, mainly expressing 70 kDa and 130 kDa crystal proteins. The LC₅₀ value of the Bt GS57 strain against the S. exigua larvae was 0.339 mg mL^–1. Physiological and biochemical reactions showed that Bt GS57 belongs to B.t. var. thuringiensis. In addition, we found that B. thuringiensis can cause a dynamic change in the gut microbiota of S. exigua, with a significant reduction in bacterial diversity and a substantial increase in bacterial load. In turn, loss of gut microbiota significantly decreased the B. thuringiensis susceptibility of S. exigua larvae. Our findings reveal the vital contribution of the gut microbiota in B. thuringiensis-killing activity, providing new insights into the mechanisms of B. thuringiensis pathogenesis in insects.

Putative plasmid prophages of Bacillus cereus sensu lato may hold the key to undiscovered phage diversity

Bacteriophages are bacterial viruses and the most abundant biological entities on Earth. Temperate bacteriophages can form prophages stably maintained in the host population: they either integrate into the host genome or replicate as plasmids in the host cytoplasm. As shown, tailed temperate bacteriophages may form circular plasmid prophages in many bacterial species of the taxa Firmicutes, Gammaproteobacteria and Spirochaetes. The actual number of such prophages is thought to be underestimated for two main reasons: first, in bacterial whole genome-sequencing assemblies, they are difficult to distinguish from actual plasmids; second, there is an absence of experimental studies which are vital to confirm their existence. In Firmicutes, such prophages appear to be especially numerous. In the present study, we identified 23 genomes from species of the Bacillus cereus group that were deposited in GenBank as plasmids and may belong to plasmid prophages with little or no homology to known viruses. We consider these putative prophages worth experimental assays since it will broaden our knowledge of phage diversity and suggest that more attention be paid to such molecules in all bacterial sequencing projects as this will help in identifying previously unknown phages.

A putative exosporium lipoprotein GBAA0190 of Bacillus anthracis as a potential anthrax vaccine candidate

BACKGROUND

Bacillus ancthracis causes cutaneous, pulmonary, or gastrointestinal forms of anthrax. B. anthracis is a pathogenic bacterium that is potentially to be used in bioterrorism because it can be produced in the form of spores. Currently, protective antigen (PA)-based vaccines are being used for the prevention of anthrax, but it is necessary to develop more safe and effective vaccines due to their prolonged immunization schedules and adverse reactions.

METHODS

We selected the lipoprotein GBAA0190, a potent inducer of host immune response, present in anthrax spores as a novel potential vaccine candidate. Then, we evaluated its immune-stimulating activity in the bone marrow-derived macrophages (BMDMs) using enzyme-linked immunosorbent assay (ELISA) and Western blot analysis. Protective efficacy of GBAA0190 was evaluated in the guinea pig (GP) model.

RESULTS

The recombinant GBAA0190 (r0190) protein induced the expression of various inflammatory cytokines including tumor necrosis factor-α (TNF-α), interleukin-6 (IL-6), monocyte chemoattractant protein-1 (MCP-1), and macrophage inflammatory protein-1α (MIP-1α) in the BMDMs. These immune responses were mediated through toll-like receptor 1/2 via activation of mitogen-activated protein (MAP) kinase and Nuclear factor-κB (NF-κB) pathways. We demonstrated that not only immunization of r0190 alone, but also combined immunization with r0190 and recombinant PA showed significant protective efficacy against B. anthracis spore challenges in the GP model.

CONCLUSIONS

Our results suggest that r0190 may be a potential target for anthrax vaccine.

Bacterial pathogens: threat or treat (a review on bioactive natural products from bacterial pathogens)

Covering: up to the second quarter of 2020 Threat or treat? While pathogenic bacteria pose significant threats, they also represent a huge reservoir of potential pharmaceuticals to treat various diseases. The alarming antimicrobial resistance crisis and the dwindling clinical pipeline urgently call for the discovery and development of new antibiotics. Pathogenic bacteria have an enormous potential for natural products drug discovery, yet they remained untapped and understudied. Herein, we review the specialised metabolites isolated from entomopathogenic, phytopathogenic, and human pathogenic bacteria with antibacterial and antifungal activities, highlighting those currently in pre-clinical trials or with potential for drug development. Selected unusual biosynthetic pathways, the key roles they play (where known) in various ecological niches are described. We also provide an overview of the mode of action (molecular target), activity, and minimum inhibitory concentration (MIC) towards bacteria and fungi. The exploitation of pathogenic bacteria as a rich source of antimicrobials, combined with the recent advances in genomics and natural products research methodology, could pave the way for a new golden age of antibiotic discovery. This review should serve as a compendium to communities of medicinal chemists, organic chemists, natural product chemists, biochemists, clinical researchers, and many others interested in the subject.

Tracking Antimicrobial Resistance Determinants in Diarrheal Pathogens: A Cross-Institutional Pilot Study

Infectious diarrhea affects over four billion individuals annually and causes over a million deaths each year. Though not typically prescribed for treatment of uncomplicated diarrheal disease, antimicrobials serve as a critical part of the armamentarium used to treat severe or persistent cases. Due to widespread over- and misuse of antimicrobials, there has been an alarming increase in global resistance, for which a standardized methodology for geographic surveillance would be highly beneficial. To demonstrate that a standardized methodology could be used to provide molecular surveillance of antimicrobial resistance (AMR) genes, we initiated a pilot study to test 130 diarrheal pathogens (Campylobacter spp., Escherichia coli, Salmonella, and Shigella spp.) from the USA, Peru, Egypt, Cambodia, and Kenya for the presence/absence of over 200 AMR determinants. We detected a total of 55 different determinants conferring resistance to ten different categories of antimicrobials: genes detected in ≥ 25 samples included bla_TEM, tet(A), tet(B), mac(A), mac(B), aadA1/A2, strA, strB, sul1, sul2, qacEΔ1, cmr, and dfrA1. The number of determinants per strain ranged from none (several Campylobacter spp. strains) to sixteen, with isolates from Egypt harboring a wider variety and greater number of genes per isolate than other sites. Two samples harbored carbapenemase genes, bla_OXA-48 or bla_NDM. Genes conferring resistance to azithromycin (ere(A), mph(A)/mph(K), erm(B)), a first-line therapeutic for severe diarrhea, were detected in over 10% of all Enterobacteriaceae tested: these included >25% of the Enterobacteriaceae from Egypt and Kenya. Forty-six percent of the Egyptian Enterobacteriaceae harbored genes encoding CTX-M-1 or CTX-M-9 families of extended-spectrum β-lactamases. Overall, the data provide cross-comparable resistome information to establish regional trends in support of international surveillance activities and potentially guide geospatially informed medical care.

A Survey of Antimicrobial Resistance Determinants in Category A Select Agents, Exempt Strains, and Near-Neighbor Species

A dramatic increase in global antimicrobial resistance (AMR) has been well documented. Of particular concern is the dearth of information regarding the spectrum and prevalence of AMR within Category A Select Agents. Here, we performed a survey of horizontally and vertically transferred AMR determinants among Category A agents and their near neighbors. Microarrays provided broad spectrum screening of 127 Francisella spp., Yersinia spp., and Bacillus spp. strains for the presence/absence of 500+ AMR genes (or families of genes). Detecting a broad variety of AMR genes in each genus, microarray analysis also picked up the presence of an engineered plasmid in a Y. pestis strain. High resolution melt analysis (HRMA) was also used to assess the presence of quinolone resistance-associated mutations in 100 of these strains. Though HRMA was able to detect resistance-causing point mutations in B. anthracis strains, it was not capable of discriminating these point mutations from other nucleotide substitutions (e.g., arising from sequence differences in near neighbors). Though these technologies are well-established, to our knowledge, this is the largest survey of Category A agents and their near-neighbor species for genes covering multiple mechanisms of AMR.

Formate-nitrite transporters carrying nonprotonatable amide amino acids instead of a central histidine maintain pH-dependent transport

Microbial formate-nitrite transporter-type proteins (FNT) exhibit dual transport functionality. At neutral pH, electrogenic anion currents are detectable, whereas upon acidification transport of the neutral, protonated monoacid predominates. Physiologically, FNT-mediated proton co-transport is vital when monocarboxylic acid products of the energy metabolism, such as l-lactate, are released from the cell. Accordingly, Plasmodium falciparum malaria parasites can be killed by small-molecule inhibitors of PfFNT. Two opposing hypotheses on the site of substrate protonation are plausible. The proton relay mechanism postulates proton transfer from a highly conserved histidine centrally positioned in the transport path. The dielectric slide mechanism assumes decreasing acidity of substrates entering the lipophilic vestibules and protonation via the bulk water. Here, we defined the transport mechanism of the FNT from the amoebiasis parasite Entamoeba histolytica, EhFNT, and also show that BtFdhC from Bacillus thuringiensis is a functional formate transporter. Both FNTs carry a nonprotonatable amide amino acid, asparagine or glutamine, respectively, at the central histidine position. Despite having a nonprotonatable residue, EhFNT displayed the same substrate selectivity for larger monocarboxylates including l-lactate, a low substrate affinity as is typical for FNTs, and, strikingly, proton motive force-dependent transport as observed for PfFNT harboring a central histidine. These results argue against a proton relay mechanism, indicating that substrate protonation must occur outside of the central histidine region, most likely in the vestibules. Furthermore, EhFNT is the sole annotated FNT in the Entamoeba genome suggesting that it could be a putative new drug target with similar utility as that of the malarial PfFNT.

Genes under positive selection in the core genome of pathogenic Bacillus cereus group members

In this comparative genomics study our aim was to unravel genes under positive selection in the core genome of the Bacillus cereus group. Indeed, the members of this group share close genetic relationships but display a rather large phenotypic and ecological diversity, providing a unique opportunity for studying how genomic changes reflect ecological adaptation during the divergence of a bacterial group. For this purpose, we screened ten completely sequenced genomes of four pathogenic Bacillus species, finding that 254 out of 3093 genes have codon sites with d_N/d_S (ω) values above one. These results remained unchanged after having disentangled the confounding effects of recombination and selection signature in a Bayesian framework. The presumably adaptive nucleotide polymorphisms are distributed over a wide range of biological functions, such as antibiotic resistance, DNA repair, nutrient uptake, metabolism, cell wall assembly and spore structure. Our results indicate that adaptation to animal hosts, whether as pathogens, saprophytes or symbionts, is the major driving force in the evolution of the Bacillus cereus group. Future work should seek to understand the evolutionary dynamics of both core and accessory genes in an integrative framework to ultimately unravel the key networks involved in host adaptation.

Bacillus thuringiensis: a successful insecticide with new environmental features and tidings

Bacillus thuringiensis (Bt) is known as the most successful microbial insecticide against different orders of insect pests in agriculture and medicine. Moreover, Bt toxin genes also have been efficiently used to enhance resistance to insect pests in genetically modified crops. In light of the scientific advantages of new molecular biology technologies, recently, some other new potentials of Bt have been explored. These new environmental features include the toxicity against nematodes, mites, and ticks, antagonistic effects against plant and animal pathogenic bacteria and fungi, plant growth-promoting activities (PGPR), bioremediation of different heavy metals and other pollutants, biosynthesis of metal nanoparticles, production of polyhydroxyalkanoate biopolymer, and anticancer activities (due to parasporins). This review comprehensively describes recent advances in the Bt whole-genome studies, the last updated known Bt toxins and their functions, and application of cry genes in plant genetic engineering. Moreover, the review thoroughly describes the new features of Bt which make it a suitable cell factory that might be used for production of different novel valuable bioproducts.

Hot, humid air decontamination of a C-130 aircraft contaminated with spores of two acrystalliferous Bacillus thuringiensis strains, surrogates for Bacillus anthracis

AIM

To develop test methods and evaluate survival of Bacillus thuringiensis kurstaki cry(-) HD-1 and B. thuringiensis Al Hakam spores after exposure to hot, humid air inside of a C-130 aircraft.

METHODS AND RESULTS

Bacillus thuringiensis spores were either pre-inoculated on 1 × 2 or 2 × 2 cm substrates or aerosolized inside the cargo hold of a C-130 and allowed to dry. Dirty, complex surfaces (10 × 10 cm) swabbed after spore dispersal showed a deposition of 8-10 log10 m(-2) through the entire cargo hold. After hot, humid air decontamination at 75-80°C, 70-90% relative humidity for 7 days, 87 of 98 test swabs covering 0·98 m(2) , showed complete spore inactivation. There was a total of 1·67 log10 live CFU detected in 11 of the test swabs. Spore inactivation in the 98 test swabs was measured at 7·06 log10 m(-2) .

CONCLUSIONS

Laboratory test methods for hot, humid air decontamination were scaled for a large-scale aircraft field test. The C-130 field test demonstrated that hot, humid air can be successfully used to decontaminate an aircraft.

SIGNIFICANCE AND IMPACT OF THE STUDY

Transition of a new technology from research and development to acquisition at a Technology Readiness Level 7 is unprecedented.

Genomic and transcriptomic insights into the efficient entomopathogenicity of Bacillus thuringiensis

Bacillus thuringiensis has been globally used as a microbial pesticide for over 70 years. However, information regarding its various adaptions and virulence factors and their roles in the entomopathogenic process remains limited. In this work, we present the complete genomes of two industrially patented Bacillus thuringiensis strains (HD-1 and YBT-1520). A comparative genomic analysis showed a larger and more complicated genome constitution that included novel insecticidal toxicity-related genes (ITRGs). All of the putative ITRGs were summarized according to the steps of infection. A comparative genomic analysis showed that highly toxic strains contained significantly more ITRGs, thereby providing additional strategies for infection, immune evasion, and cadaver utilization. Furthermore, a comparative transcriptomic analysis suggested that a high expression of these ITRGs was a key factor in efficient entomopathogenicity. We identified an active extra urease synthesis system in the highly toxic strains that may aid B. thuringiensis survival in insects (similar to previous results with well-known pathogens). Taken together, these results explain the efficient entomopathogenicity of B. thuringiensis. It provides novel insights into the strategies used by B. thuringiensis to resist and overcome host immune defenses and helps identify novel toxicity factors.

Test methods and response surface models for hot, humid air decontamination of materials contaminated with dirty spores of Bacillus anthracis ∆Sterne and Bacillus thuringiensis Al Hakam

AIMS

To develop test methods and evaluate survival of Bacillus anthracis ∆Sterne or Bacillus thuringiensis Al Hakam on materials contaminated with dirty spore preparations after exposure to hot, humid air using response surface modelling.

METHODS AND RESULTS

Spores (>7 log10 ) were mixed with humic acid + spent sporulation medium (organic debris) or kaolin (dirt debris). Spore samples were then dried on five different test materials (wiring insulation, aircraft performance coating, anti-skid, polypropylene, and nylon). Inoculated materials were tested with 19 test combinations of temperature (55, 65, 75°C), relative humidity (70, 80, 90%) and time (1, 2, 3 days). The slowest spore inactivation kinetics was on nylon webbing and/or after addition of organic debris.

CONCLUSIONS

Hot, humid air effectively decontaminates materials contaminated with dirty Bacillus spore preparations; debris and material interactions create complex decontamination kinetic patterns; and B. thuringiensis Al Hakam is a realistic surrogate for B. anthracis.

SIGNIFICANCE AND IMPACT OF THE STUDY

Response surface models of hot, humid air decontamination were developed which may be used to select decontamination parameters for contamination scenarios including aircraft.

Whole Genome Phylogeny of Bacillus by Feature Frequency Profiles (FFP)

Fifty complete Bacillus genome sequences and associated plasmids were compared using the “feature frequency profile” (FFP) method. The resulting whole-genome phylogeny supports the placement of three Bacillus species (B. thuringiensis, B. anthracis and B. cereus) as a single clade. The monophyletic status of B. anthracis was strongly supported by the analysis. FFP proved to be more effective in inferring the phylogeny of Bacillus than methods based on single gene sequences [16s rRNA gene, GryB (gyrase subunit B) and AroE (shikimate-5-dehydrogenase)] analyses. The findings of FFP analysis were verified using kSNP v2 (alignment-free sequence analysis method) and Harvest suite (core genome sequence alignment method).

Genome Sequence of Bacillus endophyticus and Analysis of Its Companion Mechanism in the Ketogulonigenium vulgare-Bacillus Strain Consortium

Bacillus strains have been widely used as the companion strain of Ketogulonigenium vulgare in the process of vitamin C fermentation. Different Bacillus strains generate different effects on the growth of K. vulgare and ultimately influence the productivity. First, we identified that Bacillus endophyticus Hbe603 was an appropriate strain to cooperate with K. vulgare and the product conversion rate exceeded 90% in industrial vitamin C fermentation. Here, we report the genome sequencing of the B. endophyticus Hbe603 industrial companion strain and speculate its possible advantage in the consortium. The circular chromosome of B. endophyticus Hbe603 has a size of 4.87 Mb with GC content of 36.64% and has the highest similarity with that of Bacillus megaterium among all the bacteria with complete genomes. By comparing the distribution of COGs with that of Bacillus thuringiensis, Bacillus cereus and B. megaterium, B. endophyticus has less genes related to cell envelope biogenesis and signal transduction mechanisms, and more genes related to carbohydrate transport and metabolism, energy production and conversion, as well as lipid transport and metabolism. Genome-based functional studies revealed the specific capability of B. endophyticus in sporulation, transcription regulation, environmental resistance, membrane transportation, extracellular proteins and nutrients synthesis, which would be beneficial for K. vulgare. In particular, B. endophyticus lacks the Rap-Phr signal cascade system and, in part, spore coat related proteins. In addition, it has specific pathways for vitamin B12 synthesis and sorbitol metabolism. The genome analysis of the industrial B. endophyticus will help us understand its cooperative mechanism in the K. vulgare-Bacillus strain consortium to improve the fermentation of vitamin C.

Comparative analysis of genomics and proteomics in Bacillus thuringiensis 4.0718

Bacillus thuringiensis is a widely used biopesticide that produced various insecticidal active substances during its life cycle. Separation and purification of numerous insecticide active substances have been difficult because of the relatively short half-life of such substances. On the other hand, substances can be synthetized at different times during development, so samples at different stages have to be studied, further complicating the analysis. A dual genomic and proteomic approach would enhance our ability to identify such substances, and particularily using mass spectrometry-based proteomic methods. The comparative analysis for genomic and proteomic data have showed that not all of the products deduced from the annotated genome could be identified among the proteomic data. For instance, genome annotation results showed that 39 coding sequences in the whole genome were related to insect pathogenicity, including five cry genes. However, Cry2Ab, Cry1Ia, Cytotoxin K, Bacteriocin, Exoenzyme C3 and Alveolysin could not be detected in the proteomic data obtained. The sporulation-related proteins were also compared analysis, results showed that the great majority sporulation-related proteins can be detected by mass spectrometry. This analysis revealed Spo0A~P, SigF, SigE(+), SigK(+) and SigG(+), all known to play an important role in the process of spore formation regulatory network, also were displayed in the proteomic data. Through the comparison of the two data sets, it was possible to infer that some genes were silenced or were expressed at very low levels. For instance, found that cry2Ab seems to lack a functional promoter while cry1Ia may not be expressed due to the presence of transposons. With this comparative study a relatively complete database can be constructed and used to transform hereditary material, thereby prompting the high expression of toxic proteins. A theoretical basis is provided for constructing highly virulent engineered bacteria and for promoting the application of proteogenomics in the life sciences.

High-resolution melt analysis for rapid comparison of bacterial community compositions

In the study of bacterial community composition, 16S rRNA gene amplicon sequencing is today among the preferred methods of analysis. The cost of nucleotide sequence analysis, including requisite computational and bioinformatic steps, however, takes up a large part of many research budgets. High-resolution melt (HRM) analysis is the study of the melt behavior of specific PCR products. Here we describe a novel high-throughput approach in which we used HRM analysis targeting the 16S rRNA gene to rapidly screen multiple complex samples for differences in bacterial community composition. We hypothesized that HRM analysis of amplified 16S rRNA genes from a soil ecosystem could be used as a screening tool to identify changes in bacterial community structure. This hypothesis was tested using a soil microcosm setup exposed to a total of six treatments representing different combinations of pesticide and fertilization treatments. The HRM analysis identified a shift in the bacterial community composition in two of the treatments, both including the soil fumigant Basamid GR. These results were confirmed with both denaturing gradient gel electrophoresis (DGGE) analysis and 454-based 16S rRNA gene amplicon sequencing. HRM analysis was shown to be a fast, high-throughput technique that can serve as an effective alternative to gel-based screening methods to monitor microbial community composition.

The impact of inducing germination of Bacillus anthracis and Bacillus thuringiensis spores on potential secondary decontamination strategies

AIMS

Decontamination and remediation of a site contaminated by the accidental or intentional release of fully virulent Bacillus anthracis spores are difficult, costly and potentially damaging to the environment. Development of novel decontamination strategies that have minimal environmental impacts remains a high priority. Although ungerminated spores are amongst the most resilient organisms known, once exposed to germinants, the germinating spores, in some cases, become susceptible to antimicrobial environments. We evaluated the concept that once germinated, B. anthracis spores would be less hazardous and significantly easier to remediate than ungerminated dormant spores.

METHODS AND RESULTS

Through in vitro germination and sensitivity assays, we demonstrated that upon germination, B. anthracis Ames spores and Bacillus thuringiensis Al Hakam spores (serving as a surrogate for B. anthracis) become susceptible to environmental stressors. The majority of these germinated B. anthracis and B. thuringiensis spores were nonviable after exposure to a defined minimal germination-inducing solution for prolonged periods of time. Additionally, we examined the impact of potential secondary disinfectant strategies including bleach, hydrogen peroxide, formaldehyde and artificial UV-A, UV-B and UV-C radiation, employed after a 60-min germination-induction step. Each secondary disinfectant employs a unique mechanism of killing; as a result, germination-induction strategies are better suited for some secondary disinfectants than others.

CONCLUSIONS

These results provide evidence that the deployment of an optimal combination strategy of germination-induction/secondary disinfection may be a promising aspect of wide-area decontamination following a B. anthracis contamination event.

SIGNIFICANCE AND IMPACT OF THE STUDY

By inducing spores to germinate, our data confirm that the resulting cells exhibit sensitivities that can be leveraged when paired with certain decontamination measures. This increased susceptibility could be exploited to devise more efficient and safe decontamination measures and may obviate the need for more stringent methods that are currently in place.

Bacillus thuringiensis HD-1 Cry- : development of a safe, non-insecticidal simulant for Bacillus anthracis

AIMS

A representative simulant for spores of Bacillus anthracis is needed for field testing. Bacillus thuringiensis is gaining recognition as a suitable organism. A strain that does not form the insecticidal, parasporal crystals that are characteristic of this species is a more accurate physical representative of B. anthracis spores. We developed noninsecticidal derivatives of two isolates of B. thuringiensis HD-1.

METHODS AND RESULTS

Two plasmid-cured derivatives of B. thuringiensis HD-1, unable to make crystal toxins ('Cry(-) '), were isolated. These isolates and the existing Cry(-) strain, B. thuringiensis Al Hakam, were probed with PCR assays against the known insecticidal genes cry, vip and cyt. Their genomic DNA was sequenced to demonstrate a lack of insecticidal genes. This was confirmed by bioassays against a number of invertebrate species. Real-time PCR assays were developed to identify the B. thuringiensis HD-1 Cry(-) derivatives and an effective differential and selective medium was assessed.

CONCLUSIONS

All three Cry(-) isolates are devoid of known insecticidal determinants. The B. thuringiensis HD-1 Cry(-) derivatives can easily be recovered from soil and identified by PCR with some selectivity.

SIGNIFICANCE AND IMPACT OF THE STUDY

The B. thuringiensis HD-1 Cry(-) derivatives represent accurate, nongenetically manipulated simulants for B. anthracis with excellent human and environmental safety records.

Evolution and dynamics of megaplasmids with genome sizes larger than 100 kb in the Bacillus cereus group

BACKGROUND

Plasmids play a crucial role in the evolution of bacterial genomes by mediating horizontal gene transfer. However, the origin and evolution of most plasmids remains unclear, especially for megaplasmids. Strains of the Bacillus cereus group contain up to 13 plasmids with genome sizes ranging from 2 kb to 600 kb, and thus can be used to study plasmid dynamics and evolution.

RESULTS

This work studied the origin and evolution of 31 B. cereus group megaplasmids (>100 kb) focusing on the most conserved regions on plasmids, minireplicons. Sixty-five putative minireplicons were identified and classified to six types on the basis of proteins that are essential for replication. Twenty-nine of the 31 megaplasmids contained two or more minireplicons. Phylogenetic analysis of the protein sequences showed that different minireplicons on the same megaplasmid have different evolutionary histories. Therefore, we speculated that these megaplasmids are the results of fusion of smaller plasmids. All plasmids of a bacterial strain must be compatible. In megaplasmids of the B. cereus group, individual minireplicons of different megaplasmids in the same strain belong to different types or subtypes. Thus, the subtypes of each minireplicon they contain may determine the incompatibilities of megaplasmids. A broader analysis of all 1285 bacterial plasmids with putative known minireplicons whose complete genome sequences were available from GenBank revealed that 34% (443 plasmids) of the plasmids have two or more minireplicons. This indicates that plasmid fusion events are general among bacterial plasmids.

CONCLUSIONS

Megaplasmids of B. cereus group are fusion of smaller plasmids, and the fusion of plasmids likely occurs frequently in the B. cereus group and in other bacterial taxa. Plasmid fusion may be one of the major mechanisms for formation of novel megaplasmids in the evolution of bacteria.

Comparative genome analysis of Enterobacter cloacae

The Enterobacter cloacae species includes an extremely diverse group of bacteria that are associated with plants, soil and humans. Publication of the complete genome sequence of the plant growth-promoting endophytic E. cloacae subsp. cloacae ENHKU01 provided an opportunity to perform the first comparative genome analysis between strains of this dynamic species. Examination of the pan-genome of E. cloacae showed that the conserved core genome retains the general physiological and survival genes of the species, while genomic factors in plasmids and variable regions determine the virulence of the human pathogenic E. cloacae strain; additionally, the diversity of fimbriae contributes to variation in colonization and host determination of different E. cloacae strains. Comparative genome analysis further illustrated that E. cloacae strains possess multiple mechanisms for antagonistic action against other microorganisms, which involve the production of siderophores and various antimicrobial compounds, such as bacteriocins, chitinases and antibiotic resistance proteins. The presence of Type VI secretion systems is expected to provide further fitness advantages for E. cloacae in microbial competition, thus allowing it to survive in different environments. Competition assays were performed to support our observations in genomic analysis, where E. cloacae subsp. cloacae ENHKU01 demonstrated antagonistic activities against a wide range of plant pathogenic fungal and bacterial species.

D-cycloserine or similar physiochemical compounds may be uniquely suited for use in Bacillus anthracis spore decontamination strategies

AIMS

As observed in the aftermath of the anthrax attacks of 2001, decontamination and remediation of a site contaminated by the accidental or intentional release of Bacillus anthracis spores is difficult, costly and potentially damaging to the environment. The identification of novel strategies that neutralize the threat of spores while minimizing environmental damage remains a high priority. We investigated the efficacy of d-cycloserine (DCS), an antibiotic and inhibitor of the spore-associated enzyme (alanine racemase) responsible for converting l-alanine to d-alanine, as a spore germination enhancer and antimicrobial agent.

METHODS AND RESULTS

We characterized the impact of DCS exposure on both germinating spores and vegetative cells of fully virulent B. anthracis by evaluating spore germination kinetics, determining the minimum inhibitory concentrations (MICs) required to affect growth of the bacteria and performing macrophage viability assays. DCS enhanced germination induced by l-alanine and also efficiently killed the newly germinated spores. Furthermore, DCS proved nontoxic to macrophages at concentrations that provided protection from the killing effects of spores. Similar tests were conducted with Bacillus thuringiensis (subspecies kurstaki and Al Hakam) to determine its potential as a possible surrogate for B. anthracis field trials. Bacillus thuringiensis spores responded in a similar manner to B. anthracis spores when exposed to DCS.

CONCLUSIONS

These results further support that DCS augments the germination response of spores in the presence of l-alanine but also reveal that DCS is bactericidal towards germinating spores.

SIGNIFICANCE AND IMPACT OF THE STUDY

DCS (or similar compounds) may be uniquely suited for use as part of decontamination strategies by augmenting the induction of spore germination and then rendering the germinated spores nonviable.

Complete Genome Sequence of Bacillus thuringiensis Strain 407 Cry-

Bacillus thuringiensis is an insect pathogen that has been used widely as a biopesticide. Here, we report the genome sequence of strain 407 Cry-, which is used to study the genetic determinants of pathogenicity. The genome consists of a 5.5-Mb chromosome and nine plasmids, including a novel 502-kb megaplasmid.

Genetic barcodes for improved environmental tracking of an anthrax simulant

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

Thousands of missed genes found in bacterial genomes and their analysis with COMBREX

Background

The dramatic reduction in the cost of sequencing has allowed many researchers to join in the effort of sequencing and annotating prokaryotic genomes. Annotation methods vary considerably and may fail to identify some genes. Here we draw attention to a large number of likely genes missing from annotations using common tools such as Glimmer and BLAST.

Results

By analyzing 1,474 prokaryotic genome annotations in GenBank, we identify 13,602 likely missed genes that are homologs to non-hypothetical proteins, and 11,792 likely missed genes that are homologs only to hypothetical proteins, yet have supporting evidence of their protein-coding nature from COMBREX, a newly created gene function database. We also estimate the likelihood that each potential missing gene found is a genuine protein-coding gene using COMBREX.

Conclusions

Our analysis of the causes of missed genes suggests that larger annotation centers tend to produce annotations with fewer missed genes than smaller centers, and many of the missed genes are short genes <300 bp. Over 1,000 of the likely missed genes could be associated with phenotype information available in COMBREX. 359 of these genes, found in pathogenic organisms, may be potential targets for pharmaceutical research. The newly identified genes are available on COMBREX’s website.

Reviewers

This article was reviewed by Daniel Haft, Arcady Mushegian, and M. Pilar Francino (nominated by David Ardell).

Investigation of the Cry4B-prohibitin interaction in Aedes aegypti cells

Bacillus thuringiensis (Bt) produces insecticidal toxins active against insects. Cry4B, one of the major insecticidal toxins produced by Bt subsp. israelensis, is highly toxic to mosquitoes in the genus Aedes: the major vectors of dengue, yellow fever, and chikungunya. Previous work has shown that Cry4B binds to several mid-gut membrane proteins in Aedes aegypti larvae including prohibitin, a protein recently identified as a receptor that also mediates entry of dengue virus into Aedes cells. This study confirms the interaction between Cry4B and prohibitin by co-immunoprecipitation analysis and demonstrates colocalization of prohibitin and Cry4B by confocal microscopy. While activated Cry4B toxin showed high larvicidal activity, it was not cytotoxic to two Aedes cell lines, allowing determination of its effect on dengue virus infectivity in the absence of Cry4B-induced cell lysis. Pre-exposure of Aedes cells to Cry4B resulted in a significant reduction in the number of infected cells compared to untreated cells.

Genomic characterization of the Bacillus cereus sensu lato species: backdrop to the evolution of Bacillus anthracis

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.

Mining new crystal protein genes from Bacillus thuringiensis on the basis of mixed plasmid-enriched genome sequencing and a computational pipeline

We have designed a high-throughput system for the identification of novel crystal protein genes (cry) from Bacillus thuringiensis strains. The system was developed with two goals: (i) to acquire the mixed plasmid-enriched genomic sequence of B. thuringiensis using next-generation sequencing biotechnology, and (ii) to identify cry genes with a computational pipeline (using BtToxin_scanner). In our pipeline method, we employed three different kinds of well-developed prediction methods, BLAST, hidden Markov model (HMM), and support vector machine (SVM), to predict the presence of Cry toxin genes. The pipeline proved to be fast (average speed, 1.02 Mb/min for proteins and open reading frames [ORFs] and 1.80 Mb/min for nucleotide sequences), sensitive (it detected 40% more protein toxin genes than a keyword extraction method using genomic sequences downloaded from GenBank), and highly specific. Twenty-one strains from our laboratory's collection were selected based on their plasmid pattern and/or crystal morphology. The plasmid-enriched genomic DNA was extracted from these strains and mixed for Illumina sequencing. The sequencing data were de novo assembled, and a total of 113 candidate cry sequences were identified using the computational pipeline. Twenty-seven candidate sequences were selected on the basis of their low level of sequence identity to known cry genes, and eight full-length genes were obtained with PCR. Finally, three new cry-type genes (primary ranks) and five cry holotypes, which were designated cry8Ac1, cry7Ha1, cry21Ca1, cry32Fa1, and cry21Da1 by the B. thuringiensis Toxin Nomenclature Committee, were identified. The system described here is both efficient and cost-effective and can greatly accelerate the discovery of novel cry genes.

Complete genome sequence of the highly hemolytic strain Bacillus cereus F837/76

Highly hemolytic strain Bacillus cereus F837/76 was isolated in 1976 from a contaminated prostate wound. The complete nucleotide sequence of this strain reported here counts nearly 36,500 single-nucleotide differences from the closest sequenced strain, Bacillus thuringiensis Al Hakam. F827/76 also contains a 10-kb plasmid that was not detected in the Al Hakam strain.

Selectivity, directionality, and promiscuity in peptide processing from a Bacillus sp. Al Hakam cyclodehydratase

The thiazole/oxazole-modified microcins (TOMMs) represent a burgeoning class of ribosomal natural products decorated with thiazoles and (methyl)oxazoles originating from cysteines, serines, and threonines. The ribosomal nature of TOMMs allows for the generation of derivative products from mutations in the amino acid sequence of the precursor peptide, which ultimately manifest in differing structures and, sometimes, biological functions. Employing a TOMM system for the purpose of creating new structures and functions via combinatorial biosynthesis requires processing machinery that can tolerate highly variable substrates. In this study, TOMM enzymatic promiscuity was assessed using a currently uncharacterized cluster in Bacillus sp. Al Hakam. As determined by Fourier transform tandem mass spectrometry (FT-MS/MS), azole rings were formed in both a regio- and chemoselective fashion. Cognate and noncognate precursor peptides were modified in an overall C- to N-terminal directionality, which to date is unique among characterized ribosomal natural products. Studies focused on the inherent promiscuity of the biosynthetic machinery elucidated a modest bias for glycine at the preceding (-1) position and a remarkable flexibility in the following (+1) position, even allowing for the incorporation of charged amino acids and bisheterocyclization. Two unnatural substrates were utilized as the conclusive test of substrate flexibility, of which both were processed in a predictable fashion. A greater understanding of substrate processing and enzymatic tolerance toward unnatural substrates will prove beneficial when designing combinatorial libraries to screen for artificial TOMMs that exhibit desired activities.

Sequence Analysis of Inducible Prophage phIS3501 Integrated into the Haemolysin II Gene of Bacillus thuringiensis var israelensis ATCC35646

Diarrheic food poisoning by bacteria of the Bacillus cereus group is mostly due to several toxins encoded in the genomes. One of them, cytotoxin K, was recently identified as responsible for severe necrotic syndromes. Cytotoxin K is similar to a class of proteins encoded by genes usually annotated as haemolysin II (hlyII) in the majority of genomes of the B. cereus group. The partially sequenced genome of Bacillus thuringiensis var israelensis ATCC35646 contains several potentially induced prophages, one of them integrated into the hlyII gene. We determined the complete sequence and established the genomic organization of this prophage-designated phIS3501. During induction of excision of this prophage with mitomycin C, intact hlyII gene is formed, thus providing to cells a genetic ability to synthesize the active toxin. Therefore, this prophage, upon its excision, can be implicated in the regulation of synthesis of the active toxin and thus in the virulence of bacterial host. A generality of selection for such systems in bacterial pathogens is indicated by the similarity of this genetic arrangement to that of Staphylococcus aureus  β-haemolysin.

Diversity and natural functions of antibiotics produced by beneficial and plant pathogenic bacteria

Soil- and plant-associated environments harbor numerous bacteria that produce antibiotic metabolites with specific or broad-spectrum activities against coexisting microorganisms. The function and ecological importance of antibiotics have long been assumed to yield a survival advantage to the producing bacteria in the highly competitive but resource-limited soil environments through direct suppression. Although specific antibiotics may enhance producer persistence when challenged by competitors or predators in soil habitats, at subinhibitory concentrations antibiotics exhibit a diversity of other roles in the life history of the producing bacteria. Many processes modulated by antibiotics may be inherently critical to the producing bacterium, such as the acquisition of substrates or initiation of developmental changes that will ensure survival under stressful conditions. Antibiotics may also have roles in more complex interactions, including in virulence on host plants or in shaping the outcomes of multitrophic interactions. The innate functions of antibiotics to producing bacteria in their native ecosystem are just beginning to emerge, but current knowledge already reveals a breadth of activities well beyond the historical perspective of antibiotics as weaponry in microbial conflicts.

Response of Bacillus thuringiensis Al Hakam Endospores to Gas Dynamic Heating in a Shock Tube

Experiments were conducted in a gas-driven shock tube to investigate shock wave-induced damage to Bacillus thuringiensis Al Hakam endospores over a wide range of post-shock temperatures in non-oxidative gas environments. The results were compared with previous studies on B. atrophaeus and B. subtilis and demonstrate that B. thuringiensis Al Hakam exhibited a qualitatively similar response to rapid shock heating, even though this strain has a significantly different endospore structure. B. thuringiensis Al Hakam endospores were nebulized into an aqueous aerosol, which was loaded into the Stanford aerosol shock tube, and subjected to shock waves of controlled strength. Endospores experienced uniform test temperatures between 500 and 1000 K and pressures ranging from 2 atm to 7 atm for approximately 2.5 ms. During this process the bio-aerosol was monitored using in situ time-resolved laser absorption and scattering diagnostics. Additionally, shock-treated bio-aerosol samples were extracted for ex situ analysis including viability plating, flow cytometry and scanning electron microscopy (SEM) imaging. B. thuringiensis Al Hakam endospores lost the ability to form colonies at post-shock temperatures above 500 K while significant breakdown in morphology was observed only for post-shock temperatures above 700 K. While viability loss and endospore morphological deterioration adhere to a similar framework across all endospore species studied, phenomena unique to B. thuringiensis Al Hakam were noted in the SEM images and optical extinction data. This initial characterization of the response of B. thuringiensis Al Hakam spores treated with shock/blast waves shows that these methods have potential for spore inactivation and detection.

Developing an integrated proteo-genomic approach for the characterisation of biomarkers for the identification of Bacillus anthracis

Bacillus anthracis is the causative agent of anthrax, an acute and often fatal disease in humans. Due to the high genomic relatedness within the Bacillus cereus group of species it is a challenge to identify B. anthracis consistently. Alternative strategies such as proteomics coupled with mass spectrometry (MS) provide a powerful approach for biomarker discovery. However, validating and evaluating these markers, particularly for genetically homogeneous species such as B. anthracis are challenging. The objective of this study is to develop a robust biomarker discovery and validation pipeline, using proteomic methodology combined with in silico and molecular approaches, to determine a biomarker list, using B. anthracis as a model. In this exploratory study we profiled the proteome of B. anthracis and genetically related species using GeLC-Liquid Chromatography MS/MS (GeLC-LC MS/MS), identifying peptides that could be used to detect B. anthracis. Peptides were filtered to remove low quality identifications. Using comparative bioinformatic approaches, matching and searching against genomic sequence data a shortlist of peptide biomarkers was determined and validated using DNA sequencing, against a panel of closely related strains, to determine marker specificity. Further validation was performed using MS quantitation methods to assess sensitivity and specificity. A biomarker discovery pipeline was successfully developed in this study, comprising four distinct stages: proteome profiling, comparative bioinformatic validation, DNA sequencing and MS validation. Using the pipeline, 5379 peptides specific for Bacillus species and 36 peptides specific for B. anthracis were identified and validated. The 36 peptides, representing 30 proteins were derived from over 15 different clusters of orthologous group categories, including proteins involved in transcription, energy production/conservation as well as multifunctional proteins. We demonstrated that the peptide biomarkers identified in this study could be detected in a complex background, in which 0.1 μg of protein extract from B. anthracis was spiked into 9.90 μg of B. cereus protein extracts. The integration of both stable non-redundant peptides with molecular methodology for marker discovery and validation, improves the robustness of identifying and characterising candidate biomarkers for the identification of bacteria such as B. anthracis.

A pangenomic study of Bacillus thuringiensis

Bacillus thuringiensis (B. thuringiensis) is a soil-dwelling Gram-positive bacterium and its plasmid-encoded toxins (Cry) are commonly used as biological alternatives to pesticides. In a pangenomic study, we sequenced seven B. thuringiensis isolates in both high coverage and base-quality using the next-generation sequencing platform. The B. thuringiensis pangenome was extrapolated to have 4196 core genes and an asymptotic value of 558 unique genes when a new genome is added. Compared to the pangenomes of its closely related species of the same genus, B. thuringiensis pangenome shows an open characteristic, similar to B. cereus but not to B. anthracis; the latter has a closed pangenome. We also found extensive divergence among the seven B. thuringiensis genome assemblies, which harbor ample repeats and single nucleotide polymorphisms (SNPs). The identities among orthologous genes are greater than 84.5% and the hotspots for the genome variations were discovered in genomic regions of 2.3-2.8Mb and 5.0-5.6Mb. We concluded that high-coverage sequence assemblies from multiple strains, before all the gaps are closed, are very useful for pangenomic studies.

Validation of the intact zwittermicin A biosynthetic gene cluster and discovery of a complementary resistance mechanism in Bacillus thuringiensis

Zwittermicin A (ZmA) is a hybrid polyketide-nonribosomal peptide produced by certain Bacillus cereus group strains. It displays broad-spectrum antimicrobial activity. Its biosynthetic pathway in B. cereus has been proposed through analysis of the nonribosomal peptide synthetase (NRPS) and polyketide synthase (PKS) modules involved in ZmA biosynthesis. In this study, we constructed a bacterial artificial chromosome (BAC) library from Bacillus thuringiensis subsp. kurstaki strain YBT-1520 genomic DNA. The presence of known genes involved in the biosynthesis of ZmA in this BAC library was investigated by PCR techniques. Nine positive clones were identified, two of which (covering an approximately 60-kb region) could confer ZmA biosynthesis ability upon B. thuringiensis BMB171 after simultaneous transfer into this host by two compatible shuttle BAC vectors. Another previously unidentified gene cluster, named zmaWXY, was found to improve the yield of ZmA and was experimentally defined to function as a ZmA resistance transporter which expels ZmA from the cells. Putative transposase genes were detected on the flanking regions of the two gene clusters (the ZmA synthetic cluster and zmaWXY), which suggests a mobile nature of these two gene clusters. The intact ZmA gene cluster was validated, and a resistance mechanism complementary to that for zmaR (the previously identified ZmA self-resistance gene) was revealed. This study also provided a straightforward strategy to isolate and identify a huge gene cluster from Bacillus.

Rapidly progressive, fatal, inhalation anthrax-like infection in a human: case report, pathogen genome sequencing, pathology, and coordinated response

CONTEXT

Ten years ago a bioterrorism event involving Bacillus anthracis spores captured the nation's interest, stimulated extensive new research on this pathogen, and heightened concern about illegitimate release of infectious agents. Sporadic reports have described rare, fulminant, and sometimes fatal cases of pneumonia in humans and nonhuman primates caused by strains of Bacillus cereus , a species closely related to Bacillus anthracis.

OBJECTIVES

To describe and investigate a case of rapidly progressive, fatal, anthrax-like pneumonia and the overwhelming infection caused by a Bacillus species of uncertain provenance in a patient residing in rural Texas.

DESIGN

We characterized the genome of the causative strain within days of its recovery from antemortem cultures using next-generation sequencing and performed immunohistochemistry on tissues obtained at autopsy with antibodies directed against virulence proteins of B anthracis and B cereus.

RESULTS

We discovered that the infection was caused by a previously unknown strain of B cereus that was closely related to, but genetically distinct from, B anthracis . The strain contains a plasmid similar to pXO1, a genetic element encoding anthrax toxin and other known virulence factors. Immunohistochemistry demonstrated that several homologs of B anthracis virulence proteins were made in infected tissues, likely contributing to the patient's death.

CONCLUSIONS

Rapid genome sequence analysis permitted us to genetically define this strain, rule out the likelihood of bioterrorism, and contribute effectively to the institutional response to this event. Our experience strongly reinforced the critical value of deploying a well-integrated, anatomic, clinical, and genomic strategy to respond rapidly to a potential emerging, infectious threat to public health.

Transcription of gene in an acrystalliferous strain of Bacillus thuringiensis XBU001 positively regulated by the metalloprotease camelysin gene at the onset of stationary phase

The calY gene, encoding metalloprotease camelysin in the Bacillus thuringiensis acrystalliferous strain XBU001, was amplified and sequenced. The camelysin from the calY sequence was 199 amino acids in size (c. 22 000 Da). The temperature-sensitive plasmid pKESX was used to construct a metalloprotease camelysin-deficient strain of B. thuringiensis. The calY gene was replaced by an erythromycin-resistant gene in KCTF. Sodium dodecyl sulfate polyacrylamide gel electrophoresis and MS analysis showed that the metalloprotease InhA was not expressed after knocking out the gene calY. The temperature-sensitive plasmid pKPC was used to construct a metalloprotease camelysin complementation strain KCTFC. The InhA protein was found in KCTFC. Analysis of the expression of InhA in the wild-type strain KCTF12, camelysin-deficient and complementation strains indicated that inhA expression depended on camelysin. Although camelysin did not directly regulate the expression of the InhA through binding to the promoter of the inhA, the results suggest that camelysin can positively regulate the expression of the InhA protein.

Complete genome sequence of Bacillus thuringiensis serovar finitimus strain YBT-020

Bacillus thuringiensis is a gram-positive, spore-forming bacterium that forms parasporal crystals at the onset of the sporulation phase of its growth. Here, we report the complete genome sequence of B. thuringiensis serovar finitimus strain YBT-020, whose parasporal crystals consist of Cry26Aa and Cry28Aa crystal proteins and are located between the exosporium and the spore coat and remain adhering to the spore after sporulation.

Characterization of a novel temperate phage originating from a cereulide-producing Bacillus cereus strain

A novel temperate bacteriophage was isolated from a Bacillus cereus cereulide-producing strain and named vB_BceS-IEBH. vB_BceS-IEBH belongs to the Siphoviridae family. The complete genome sequence (53 kb) was determined and annotated. Eighty-seven ORFs were detected and for 28, a putative function was assigned using the ACLAME database. vB_BceS-IEBH replicates as a plasmid in the prophage state. Accordingly, a 9-kb plasmid-like region composed of 13 ORFs was identified. A fragment of around 2000 bp comprising an ORF encoding a putative plasmid replication protein was shown to be self-replicating in Bacillus thuringiensis. Mass spectrometry analysis of the purified vB_BceS-IEBH particle identified 8 structural proteins and enabled assignment of a supplementary ORF as being part of the morphogenesis module. Genome analysis further illustrates the diversity of mobile genetic elements and their plasticity within the B. cereus group.

Inference of homologous recombination in bacteria using whole-genome sequences

Bacteria and archaea reproduce clonally, but sporadically import DNA into their chromosomes from other organisms. In many of these events, the imported DNA replaces an homologous segment in the recipient genome. Here we present a new method to reconstruct the history of recombination events that affected a given sample of bacterial genomes. We introduce a mathematical model that represents both the donor and the recipient of each DNA import as an ancestor of the genomes in the sample. The model represents a simplification of the previously described coalescent with gene conversion. We implement a Monte Carlo Markov chain algorithm to perform inference under this model from sequence data alignments and show that inference is feasible for whole-genome alignments through parallelization. Using simulated data, we demonstrate accurate and reliable identification of individual recombination events and global recombination rate parameters. We applied our approach to an alignment of 13 whole genomes from the Bacillus cereus group. We find, as expected from laboratory experiments, that the recombination rate is higher between closely related organisms and also that the genome contains several broad regions of elevated levels of recombination. Application of the method to the genomic data sets that are becoming available should reveal the evolutionary history and private lives of populations of bacteria and archaea. The methods described in this article have been implemented in a computer software package, ClonalOrigin, which is freely available from http://code.google.com/p/clonalorigin/.

The genome of a Bacillus isolate causing anthrax in chimpanzees combines chromosomal properties of B. cereus with B. anthracis virulence plasmids

Anthrax is a fatal disease caused by strains of Bacillus anthracis. Members of this monophyletic species are non motile and are all characterized by the presence of four prophages and a nonsense mutation in the plcR regulator gene. Here we report the complete genome sequence of a Bacillus strain isolated from a chimpanzee that had died with clinical symptoms of anthrax. Unlike classic B. anthracis, this strain was motile and lacked the four prohages and the nonsense mutation. Four replicons were identified, a chromosome and three plasmids. Comparative genome analysis revealed that the chromosome resembles those of non-B. anthracis members of the Bacillus cereus group, whereas two plasmids were identical to the anthrax virulence plasmids pXO1 and pXO2. The function of the newly discovered third plasmid with a length of 14 kbp is unknown. A detailed comparison of genomic loci encoding key features confirmed a higher similarity to B. thuringiensis serovar konkukian strain 97-27 and B. cereus E33L than to B. anthracis strains. For the first time we describe the sequence of an anthrax causing bacterium possessing both anthrax plasmids that apparently does not belong to the monophyletic group of all so far known B. anthracis strains and that differs in important diagnostic features. The data suggest that this bacterium has evolved from a B. cereus strain independently from the classic B. anthracis strains and established a B. anthracis lifestyle. Therefore we suggest to designate this isolate as "B. cereus variety (var.) anthracis".

Genome-wide single nucleotide polymorphism typing method for identification of Bacillus anthracis species and strains among B. cereus group species

As an issue of biosecurity, species-specific genetic markers have been well characterized. However, Bacillus anthracis strain-specific information is currently not sufficient for traceability to identify the origin of the strain. By using genome-wide screening using short read mapping, we identified strain-specific single nucleotide polymorphisms (SNPs) among B. anthracis strains including Japanese isolates, and we further developed a simplified 80-tag SNP typing method for the primary investigation of traceability. These 80-tag SNPs were selected from 2,965 SNPs on the chromosome and the pXO1 and pXO2 plasmids from a total of 19 B. anthracis strains, including the available genome sequences of 17 strains in the GenBank database and 2 Japanese isolates that were sequenced in this study. Phylogenetic analysis based on 80-tag SNP typing showed a higher resolution power to discriminate 12 Japanese isolates rather than the 25 loci identified by multiple-locus variable-number tandem-repeat analysis (MLVA). In addition, the 80-tag PCR testing enabled the discrimination of B. anthracis from other B. cereus group species, helping to identify whether a suspected sample originates from the intentional release of a bioterrorism agent or environmental contamination with a virulent agent. In conclusion, 80-tag SNP typing can be a rapid and sufficient test for the primary investigation of strain origin. Subsequent whole-genome sequencing will reveal apparent strain-specific genetic markers for traceability of strains following an anthrax outbreak.

Understanding the evolutionary relationships and major traits of Bacillus through comparative genomics

BACKGROUND

The presence of Bacillus in very diverse environments reflects the versatile metabolic capabilities of a widely distributed genus. Traditional phylogenetic analysis based on limited gene sampling is not adequate for resolving the genus evolutionary relationships. By distinguishing between core and pan-genome, we determined the evolutionary and functional relationships of known Bacillus.

RESULTS

Our analysis is based upon twenty complete and draft Bacillus genomes, including a newly sequenced Bacillus isolate from an aquatic environment that we report for the first time here. Using a core genome, we were able to determine the phylogeny of known Bacilli, including aquatic strains whose position in the phylogenetic tree could not be unambiguously determined in the past. Using the pan-genome from the sequenced Bacillus, we identified functional differences, such as carbohydrate utilization and genes involved in signal transduction, which distinguished the taxonomic groups. We also assessed the genetic architecture of the defining traits of Bacillus, such as sporulation and competence, and showed that less than one third of the B. subtilis genes are conserved across other Bacilli. Most variation was shown to occur in genes that are needed to respond to environmental cues, suggesting that Bacilli have genetically specialized to allow for the occupation of diverse habitats and niches.

CONCLUSIONS

The aquatic Bacilli are defined here for the first time as a group through the phylogenetic analysis of 814 genes that comprise the core genome. Our data distinguished between genomic components, especially core vs. pan-genome to provide insight into phylogeny and function that would otherwise be difficult to achieve. A phylogeny may mask the diversity of functions, which we tried to uncover in our approach. The diversity of sporulation and competence genes across the Bacilli was unexpected based on previous studies of the B. subtilis model alone. The challenge of uncovering the novelties and variations among genes of the non-subtilis groups still remains. This task will be best accomplished by directing efforts toward understanding phylogenetic groups with similar ecological niches.

Prevalence of Bacillus anthracis-like organisms and bacteriophages in the intestinal tract of the earthworm Eisenia fetida

Stable infection of Bacillus anthracis laboratory strains with environmental bacteriophages confers survival phenotypes in soil and earthworm intestinal niches (R. Schuch and V. A. Fischetti, PLoS One 4:e6532, 2009). Here, the natural occurrence of two such B. anthracis-infective bacteriophages, Wip1 and Wip4, was examined in the intestines of Eisenia fetida earthworms as part of a 6-year longitudinal study at a Pennsylvania forest site. The Wip1 tectivirus was initially dominant before being supplanted by the Wip4 siphovirus, which was then dominant for the next 3 years. In a host range analysis of a wide-ranging group of Bacillus species and related organisms, Wip1 and Wip4 were both infective only toward B. anthracis and certain B. cereus strains. The natural host of Wip4 remained constant for 3 years and was a B. cereus strain that expressed a B. anthracis-like surface polysaccharide at septal positions on the cell surface. Next, a novel metagenomic approach was used to determine the extent to which such B. cereus- and B. anthracis-like strains are found in worms from two geographical locations. Three different enrichment strategies were used for metagenomic DNA isolation, based either on the ability of B. cereus sensu lato to form heat-resistant spores, the sensitivity of B. anthracis to the PlyG lysin, or the selective amplification of environmental phages cocultured with B. anthracis. Findings from this work indicate that B. cereus sensu lato and its phages are common inhabitants of earthworm intestines.

Cloning of circular DNAs from microorganisms using a novel plasmid capture system

Plasmid capture system (PCS) facilitates cloning and manipulation of circular double-stranded DNA. We recently developed an improved PCS (PCS-LZ) to clone relatively large DNA molecules of 30-150 kb. The PCS-LZ donor consists of a mini-F replicon and a kanamycin resistance marker between Tn7 left and Tn7 right ends. Both the replicon and marker gene of the PCS-LZ donor are transferred into target plasmid DNAs by in vitro transposition, followed by replication in E. coli. Colonies are tested for lacZ expression by blue/white screening. Circular DNAs were obtained from plasmids of Bacillus thuringiensis, genome segments of Cotesia glomerata bracovirus and polymorphic genomes of Autographa californica nucleopolyhedrovirus. PCS-LZ is a powerful tool for use in genomic analysis and mutagenesis in microorganisms including invertebrate pathogens.

The genome and variation of Bacillus anthracis

The Bacillus anthracis genome reflects its close genetic ties to Bacillus cereus and Bacillus thuringiensis but has been shaped by its own unique biology and evolutionary forces. The genome is comprised of a chromosome and two large virulence plasmids, pXO1 and pXO2. The chromosome is mostly co-linear among B. anthracis strains and even with the closest near neighbor strains. An exception to this pattern has been observed in a large inversion in an attenuated strain suggesting that chromosome co-linearity is important to the natural biology of this pathogen. In general, there are few polymorphic nucleotides among B. anthracis strains reflecting the short evolutionary time since its derivation from a B. cereus-like ancestor. The exceptions to this lack of diversity are the variable number tandem repeat (VNTR) loci that exist in genic and non genic regions of the chromosome and both plasmids. Their variation is associated with high mutability that is driven by rapid insertion and deletion of the repeats within an array. A notable example is found in the vrrC locus which is homologous to known DNA translocase genes from other bacteria.