Discrimination of Bacillus cereus and Bacillus thuringiensis with 16S rRNA and gyrB gene based PCR primers and sequencing of their annealing sites

AIMS: To evaluate the possibility for discrimination of Bacillus cereus and B. thuringiensis using 16S rRNA and gyrB gene based PCR methods, and to obtain the sequences of the primer annealing sites so that the PCR results may be explained. METHODS AND RESULTS: Based on the sequence difference in the variable region (V1) of 16S rRNA and in the gyrB gene between B. cereus and B. thuringiensis, PCR primers specific to these Bacillus spp. were designed. When these primers were used to discriminate B. cereus and B. thuringiensis, six of 82 B. cereus strains were identified as B. thuringiensis while 67 of 73 B. thuringiensis strains were identified as B. cereus. Sequence analysis of the primer annealing sites showed that there is no clear-cut difference in the V1 region of 16S rRNA, and in the gyrB gene, between B. cereus and B. thuringiensis strains. CONCLUSIONS: Although 16S rDNA based probes and gyrB gene based PCR primers have been suggested for the discrimination of B. cereus and B. thuringiensis strains, when a large number of Bacillus strains was tested, results showed that discrimination between B. cereus and B. thuringiensis is difficult. Therefore, to distinguish B. thuringiensis from B. cereus, a single feature, such as the presence of a parasporal crystal protein or cry gene, may sometimes be reliable. SIGNIFICANCE AND IMPACT OF THE STUDY: Discrimination between B. cereus and B. thuringiensis is a challenging debate to which this paper makes a contribution.     

Cloning and nucleotide sequence analysis of gyrB of Bacillus cereus, B. thuringiensis, B. mycoides, and B. anthracis and their application to the detection of B. cereus in rice

As 16S rRNA sequence analysis has proven inadequate for the differentiation of Bacillus cereus from closely related species, we employed the gyrase B gene (gyrB) as a molecular diagnostic marker. The gyrB genes of B. cereus JCM 2152(T), Bacillus thuringiensis IAM 12077(T), Bacillus mycoides ATCC 6462(T), and Bacillus anthracis Pasteur #2H were cloned and sequenced. Oligonucleotide PCR primer sets were designed from within gyrB sequences of the respective bacteria for the specific amplification and differentiation of B. cereus, B. thuringiensis, and B. anthracis. The results from the amplification of gyrB sequences correlated well with results obtained with the 16S rDNA-based hybridization study but not with the results of their phenotypic characterization. Some of the reference strains of both B. cereus (three serovars) and B. thuringiensis (two serovars) were not positive in PCR amplification assays with gyrB primers. However, complete sequencing of 1.2-kb gyrB fragments of these reference strains showed that these serovars had, in fact, lower homology than their originally designated species. We developed and tested a procedure for the specific detection of the target organism in boiled rice that entailed 15 h of preenrichment followed by PCR amplification of the B. cereus-specific fragment. This method enabled us to detect an initial inoculum of 0.24 CFU of B. cereus cells per g of boiled rice food homogenate without extracting DNA. However, a simple two-step filtration step is required to remove PCR inhibitory substances.     

Strategy for amplification and sequencing of insecticidal <Emphasis Type="Italic">cry1A</Emphasis> genes from <Emphasis Type="Italic">Bacillus thuringiensis</Emphasis>

A feasible and fully described strategy, with a detailed list of primers, for amplifying, cloning and sequencing known and potentially novel cry1A genes harboured by a Bacillus thuringiensis strain was successfully established. Based on the analysis of conserved regions of the cry1A genes, the 1AF and 1UR oligonucleotide primers were designed to amplify the whole open reading frame of these genes. The PCR products obtained revealed the successful amplification of cry1A genes from 13 B. thuringiensis strains. These bacteria were previously known to harbour at least one cry1A gene. An Argentinean B. thuringiensis isolate INTA Mo1-12 was randomly chosen for cloning and sequencing of cry1A genes by using a primer set developed in this study. Both nucleotide and amino acid sequences similarity analysis revealed that cry1Aa and cry1Ac from B. thuringiensis INTA Mo1-12 are new natural variants, showing several differences with the other known cry1A subclasses. These genes were named by the B. thuringiensis Pesticidal Crystal Protein Nomenclature Committee as cry1Aa15 and cry1Ac21 respectively.     

Detection of chromosomally located and plasmid-borne genes on 20 kb DNA fragments in parasporal crystals from <Emphasis Type="Italic">Bacillus thuringiensis</Emphasis>

The association of 20 kb heterologous DNA fragments with the parasporal crystals from native and recombinant Bacillus thuringiensis strains was analyzed, respectively. The cry2Aa10 gene cloned in plasmid pHC39 was transformed into B. thuringiensis subsp. kurstaki strains CryˉB and HD73, producing recombinant strains CryˉB(pHC39) and HD73(pHC39). SDS-PAGE and scanning electron microscopy analyses demonstrated that the recombinant CryˉB(pHC39) produced cuboidal crystals of Cry2Aa10 protoxin, while recombinant HD73(pHC39) produced both bipyramidal crystals of Cry1Ac1 protoxin and cuboidal crystals of Cry2Aa10 protoxin. Bioassay results proved that recombinant HD73(pHC39) showed higher insecticidal activity to Helicoverpa armigera than CryˉB(pHC39). It was found that 20 kb DNA fragments were present in bipyramidal and cuboidal crystals from both native and recombinant strains, and the 20 kb heterologous DNAs contained chromosome-specific and resident large plasmid-borne DNA fragments, suggesting the 20 kb heterologous DNA fragment embodied in crystals came randomly from the bacterial chromosomal and plasmid genome. This was the first investigation devoted exclusively on the origin of 20 kb DNA fragments in the parasporal crystals of B. thuringiensis. The data provides a basis for further investigation of the origin of 20 kb DNAs in the crystals and the interaction of DNA and protoxins.     

Correlation between serovars of Bacillus thuringiensis and type I beta-exotoxin production

beta-Exotoxin is a thermostable metabolite produced by some strains of Bacillus thuringiensis. Because of vertebrate toxicity, most commercial preparations of B. thuringiensis are prepared from isolates that do not produce beta-exotoxin. The aim of the present study was to find out the possible relationship between serovars of B. thuringiensis and beta-exotoxin production. A specific HPLC assay for type I beta-exotoxin has been used to detect this exotoxin in supernatants from final whole cultures of 100 strains belonging to four serovars of B. thuringiensis: thuringiensis, kurstaki, aizawai, and morrisoni. For each serovar, 25 strains randomly chosen from two Spanish collections were analyzed. Frequency of beta-exotoxin production was higher in B. thuringiensis serovar thuringiensis, whereas only two strains from serovar kurstaki showed beta-exotoxin production. None of the 25 strains belonging to serovars aizawai and morrisoni was found to produce this compound. Along with data from other studies, serovars can be classified as "common," "seldom," or "rare" beta-exotoxin producers. The serovar-dependent beta-exotoxin production is discussed in relation to the evolutionary process of serovar differentiation, the plasmid compatibility and limited plasmid exchange between serovars, and with the serovar-dependent regulation of plasmid-encoded genes.     

Detection and phylogenic analysis of one anthrax virulence plasmid pXO1 conservative open reading frame ubiquitous presented within Bacillus cereus group strains

The presence of one of the anthrax virulence plasmid pXO1 conserved fragments was analyzed in 24 Bacillus cereus and B. thuringiensis strains, including 6 B. thuringiensis subspecies, by polymerase chain reactions. Twelve out of 24 strains showed PCR-positive for an ORF101 homologous sequence. Two pXO1-ORF101-like fragments from a B. cereus B-4ac and a commercial B. thuringiensis kurstaki HD1 were cloned, sequenced and expressed in Escherichia coli. Toxicity assays revealed that the product encoded by the pXO1-ORF101-like fragment had no impact on either Vero cells or Chinese Hamster Ovary cells, suggesting that this fragment probably not contribute to enterotoxic activity. Sequence alignment of the pXO1-ORF101 from three Bacillus anthracis and ORF101-like fragments from other 12 B. cereus group isolates indicated high identity (more than 90%) and the presence of subgroup- and strain-specific SNPs among these fragments.     

dnaJ and gyrB gene sequence relationship among species and strains of genus Streptococcus

The dnaJ and gyrB nucleotide sequences were determined for members of the genus Streptococcus. The average similarity between the species tested was 76.4% (69.7-100%) for dnaJ and 75.9 (70.1-98.7%) for gyrB. These data indicated that the dnaJ and gyrB genes are more divergent and more discriminatory than the 16S rDNA gene. Furthermore, the variation in the dnaJ nucleotide sequences among the mitis group was greater than that of the gyrB nucleotide sequences, especially between Streptococcus pneumoniae and Streptococcus mitis. Subsequently, the high discrimination power of dnaJ within the mitis group was confirmed. Thus, we conclude that the dnaJ and gyrB genes are efficient alternative targets for the classification of the genus Streptococcus, and that dnaJ is suitable for phylogenetic analysis of closely related Streptococcus strains.     

PCR amplification and direct sequencing of gyrB genes with universal primers and their application to the detection and taxonomic analysis of Pseudomonas putida strains.

Degenerate PCR primers, UP-1 and UP-2r, for the amplification of DNA gyrase subunit B genes (gyrB) were designed by using consensus amino acid sequences of gyrases from Escherichia coli, Pseudomonas putida, and Bacillus subtilis. In addition to the degenerate sequences, these primers have sequences at the 5' end which allow direct sequencing of amplified PCR products. With these primers, DNA segments of the predicted size were amplified from a variety of gram-negative and gram-positive genera. The nucleotide sequences of the amplified gyrB DNA from three P. putida strains were determined directly from the amplified fragments. The base substitution frequency of gyrB between the strains of P. putida was much higher than that of the 16S rRNA gene. With a specific set of PCR primers, it was possible to amplify gyrB fragments selectively from P. putida or its subgroups. The direct sequencing method of gyrB developed in this study provides a rapid and convenient system for bacterial identification, taxonomic analysis, and monitoring of bacteria in the natural environment.     

Phylogenetic analysis of Bacillus subtilis and related taxa based on partial gyrA gene sequences

Partial gyrA sequences were determined for twelve strains belonging to Bacillus amyloliquefaciens, B. atrophaeus, B. licheniformis, B. mojavensis,B. subtilis subsp. subtilis, B. subtilissubsp. spizizenii and B. vallismortis. The average nucleotide and translated amino acid similarities for the seven type strains were 83.7 and 95.1%, respectively, whereas the corresponding value for the 16S rRNA sequences was 99.1%. All of the type strains were sharply separated; the closest relationship was found between B. atrophaeus and B. mojavensis which shared a nucleotide similarity of 95.8%. Phylogenetic trees were inferred from gyrA nucleotide sequences using the neighbor-joining, Fitch–Margoliash and maximum parsimony algorithms. The test strains were divided into four groups, which generally reflected results previously reported in restriction digest and DNA-DNA hybridization studies. It is concluded from the comparative sequence analysis that the gyrA sequences provide a firm framework for the rapid and accurate classification and identification of Bacillus subtilis and related taxa.     

Presence of methylated adenine in GATC sequences in chromosomal DNAs from Campylobacter species.

We digested chromosomal DNAs from 12 Campylobacter strains (C. jejuni, 4 strains; C. coli, 2 strains; C. fetus subsp. fetus, 2 strains; C. hyointestinalis, 2 strains; and C. upsaliensis, 2 strains) and from 4 Helicobacter strains (H. pylori, 2 strains; and H. mustelae, 2 strains) with HindIII, SstI, BamHI, DpnI, MboI, and Sau3AI. Restriction fragments were then separated by electrophoresis in 1% agarose or 10% polyacrylamide gels. Only DNAs from three Campylobacter species (C. jejuni, C. coli, and C. upsaliensis) were digested with DpnI (an enzyme that recognizes only methylated adenine in GATC sequences). We used MboI and Sau3AI to confirm these findings.     

Molecular characterization of Bacillus thuringiensis strains from Argentina

Bacillus thuringiensis INTA 7-3, INTA 51-3, INTA Mo9-5 and INTA Mo14-4 strains were obtained from Argentina and characterized by determination of serotype, toxicity, plasmid composition, insecticidal gene content ( cry and vip ) and the cloning of the single- vip3A gene of the INTA Mo9-5 strain. The serotype analysis identified the serovars tohokuensis and darmstadiensis for the INTA 51-3 and INTA Mo14-4 strains, respectively, whereas the INTA Mo9-5 strain was classified as "autoagglutinated". In contrast to the plasmid patterns of INTA 7-3, INTA 51-3 and INTA Mo9-5 (which were similar to B. thuringiensis HD-1 strain), strain INTA Mo14-4 showed a unique plasmid array. PCR analysis of the four strains revealed the presence of cry genes and vip3A genes. Interestingly, it was found that B. thuringiensis 4Q7 strain, which is a plasmid cured strain, contained vip3A genes indicating the presence of these insecticidal genes in the chromosome. Bioassays towards various lepidopteran species revealed that B. thuringiensis INTA Mo9-5 and INTA 7-3 strains were highly active. In particular, the mean LC(50) obtained against A. gemmatalis larvae with the INTA Mo9-5 and INTA 7-3 strains were 7 (5.7-8.6) and 6.7 (5.6-8.0) ppm, respectively. The INTA Mo14-4 strain was non-toxic and strain INTA 51-3 showed only a weak larvicidal activity.     

土壤中可编码乌头酸异构酶的芽胞杆菌菌株筛选及鉴定

【目的】以芽胞杆菌(Bacillus)为筛选对象,分离土壤中可编码乌头酸异构酶(aconitate isomerase,AI)的革兰氏阳性(Gram positive,G+)菌株,以丰富对AI分布的科学认识,为其生物学功能研究奠定理论与材料基础。【方法】采用土样高温预处理法、含反式乌头酸(trans-aconitic acid,TAA)唯一碳源的ACO固体平板培养法,结合16S rDNA基因序列同源性分析,筛选能够编码AI的芽胞杆菌目的菌株。【结果】共分离得到22株能够利用TAA碳源的细菌菌株,成功鉴定了其中的16株,分别为巨大芽胞杆菌(Bacillus megaterium) 2株,阿氏芽胞杆菌(Bacillus aryabhattai) 7株,短小芽胞杆菌(Bacillus pumilus) 1株,未鉴定到种的芽胞杆菌(Bacillus sp.) 6株;且它们所含AI编码基因与已知AI基因在序列上存在差异。【结论】首次证明可编码AI的芽胞杆菌细菌种类具有多样性,暗示G+细菌广泛编码AI的可能性,更新了AI几乎只在G–细菌中分布的观点,为后续深入挖掘AI基因及其生物学功能研究提供更多可用微生物资源。     

Comparative analysis of the 16S to 23S ribosomal intergenic spacer sequences of Bacillus thuringiensis strains and subspecies and of closely related species.

Bacillus thuringiensis spacer regions between the 16S and 23S rRNAs were amplified with conserved primers, designated 19-mer and 23-mer primers. A spacer region of 144 bp was determined for all of 6 B. thuringiensis strains, 7 B. thuringiensis subspecies, and 11 B. thuringiensis field isolates, as well as for the closely related species Bacillus cereus and Bacillus anthracis. Computer analysis and alignment of nucleotide sequences identified three mutations and one deletion in the intergenic spacer region (ISR) of B. thuringiensis subsp. kurstaki HD-1 when compared with ISR sequences from other subspecies. The same differences were identified between the ISR of B. thuringiensis strains and the ISR of B. cereus and B. anthracis. These minor differences do not seem to be sufficient to allow the design of a species-specific oligonucleotide probe.     

Extended genetic analysis of Brazilian isolates of Bacillus cereus and Bacillus thuringiensis

Multiple locus sequence typing (MLST) was undertaken to extend the genetic characterization of 29 isolates of Bacillus cereus and Bacillus thuringiensis previously characterized in terms of presence/absence of sequences encoding virulence factors and via variable number tandem repeat (VNTR). Additional analysis involved polymerase chain reaction for the presence of sequences (be, cytK, inA, pag, lef, cya and cap), encoding putative virulence factors, not investigated in the earlier study. MLST analysis ascribed novel and unique sequence types to each of the isolates. A phylogenetic tree was constructed from a single sequence of 2,838 bp of concatenated loci sequences. The strains were not monophyletic by analysis of any specific housekeeping gene or virulence characteristic. No clear association in relation to source of isolation or to genotypic profile based on the presence or absence of putative virulence genes could be identified. Comparison of VNTR profiling with MLST data suggested a correlation between these two methods of genetic analysis. In common with the majority of previous studies, MLST was unable to provide clarification of the basis for pathogenicity among members of the B. cereus complex. Nevertheless, our application of MLST served to reinforce the notion that B. cereus and B. thuringiensis should be considered as the same species.     

A facile analytical method for the identification of protease gene profiles from Bacillus thuringiensis strains

Five pairs of degenerate universal primers have been designed to identify the general protease gene profiles from some distinct Bacillus thuringiensis strains. Based on the PCR amplification patterns and DNA sequences of the cloned fragments, it was noted that the protease gene profiles of the three distinct strains of B. thuringiensis subsp. kurstaki HD73, tenebrionis and israelensis T14001 are varied. Seven protease genes, neutral protease B (nprB), intracellular serine protease A (ispA), extracellular serine protease (vpr), envelope-associated protease (prtH), neutral protease F (nprF), thermostable alkaline serine protease and alkaline serine protease (aprS), with known functions were identified from three distinct B. thuringiensis strains. In addition, five DNA sequences with unknown functions were also identified by this facile analytical method. However, based on the alignment of the derived protein sequences with the protein domain database, it suggested that at least one of these unknown genes, yunA, might be highly protease-related. Thus, the proposed PCR-mediated amplification design could be a facile method for identifying the protease gene profiles as well as for detecting novel protease genes of the B. thuringiensis strains.     

Kinetics of plasmid transfer among <Emphasis Type="Italic">Bacillus cereus </Emphasis>group strains within lepidopteran larvae

The cry toxin encoding plasmid pHT73 was transferred from Bacillus thuringiensis subspecies kurstaki KT0 to six B. cereus group strains in three lepidopteran (Spodoptera exigua, Plutella xyllostella and Helicoverpa armigera) larvae by conjugation. The conjugation kinetics of the plasmid was precisely studied during the larval infection using a new protocol. The infections were performed with both vegetative and sporulated strains. However, larval death only occurred when infections were made with spore and toxin preparations. Likewise, spore germinations of both donor and recipient strains were only observed in killed larvae, 44–56 h post-infection. Accordingly, kinetics showed that gene transfer between B. thuringiensis strain KT0 and other B. cereus strains only took place in dead larvae among vegetatively growing bacteria. The conjugational transfer ratios varied among different strain combinations and different larvae. The highest transfer ratio reached 5.83 × 10⁻⁶ CFU/donor between the KT0 and the AW05R recipient in Helicoverpa armigera, and all transconjugants gained the ability to produce the insecticidal crystal. These results indicated that horizontal gene transfer among B. cereus group strains might play a key role for the acquisition of extra plasmids and evolution of these strains in toxin susceptible insect larvae.     

Use of oligonucleotide probes to study the relatedness of delta-endotoxin genes among Bacillus thuringiensis subspecies and strains

Fifteen Bacillus thuringiensis strains representing 13 serotypes were screened with five oligodeoxyribonucleotide probes specific for certain regions of two published sequences and one unpublished sequence of B. thuringiensis delta-endotoxin genes. Of the 15 cultures, 14 hybridized with at least one probe; the B. thuringiensis subsp. thompsoni strain alone did not hybridize. Two B. thuringiensis subsp. kurstaki strains of commercial interest, HD-1 and NRD-12, were found to be so closely related as to be indistinguishable with this technique; the same situation was found with strains from B. thuringiensis subspp. dendrolimus and sotto. Five strains were identified as probably containing only one endotoxin gene. A probe specific for the gene from the B. thuringiensis subsp. kurstaki HD-73 strain hybridized to only 3 of the 15 cultures tested. The hybridization data suggest that the DNA sequences coding for the C-terminal region of the endotoxin protein are as well conserved as those coding for the N-terminal toxic portion.     

A comparative study of enzyme variation in Bacillus cereus and Bacillus thuringiensis

Thirty-two strains of Bacillus spp. were examined in a multilocus enzyme study by agarose gel electrophoresis. The organisms were Bacillus thuringiensis (21 strains, B. cereus (8), including two of var. mycoides, and B. megaterium (3). Strains having similar enzyme variants were grouped into zymovars. A total of 10 of 11 enzyme loci studied were polymorphic and 27 zymovars were distinguished among the 32 strains. The results were subjected to numerical analysis, phenetic affinities and genetic distances between the strains were calculated. The numerical analysis was unable to differentiate between B. thuringiensis and B. cereus. Our results indicated that based on this multilocus enzyme study these zymovars should be considered as belonging to the same species. A mycoides variant of B. cereus was the most distinctive strain studied and clearly belonged to a separate species, B. mycoides. The technique also allowed for identification of contamination and mislabelling of strains.     

A comparative study of enzyme variation in Bacillus cereus and Bacillus thuringiensis

Thirty-two strains of Bacillus spp. were examined in a multilocus enzyme study by agarose gel electrophoresis. The organisms were Bacillus thuringiensis (21 strains), B. cereus (8), including two of var. mycoides , and B. megaterium (3). Strains having similar enzyme variants were grouped into zymovars. A total of 10 of 11 enzyme loci studied were polymorphic and 27 zymovars were distinguished among the 32 strains. The results were subjected to numerical analysis, phenetic affinities and genetic distances between the strains were calculated. The numerical analysis was unable to differentiate between B. thuringiensis and B. cereus . Our results indicated that based on this multilocus enzyme study these zymovars should be considered as belonging to the same species. A mycoides variant of B. cereus was the most distinctive strain studied and clearly belonged to a separate species, B. mycoides. The technique also allowed for identification of contamination and mislabelling of strains.