Previous Article | Next Article ![]()
Applied and Environmental Microbiology, March 2005, p. 1346-1355, Vol. 71, No. 3
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.3.1346-1355.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Departamento de Biotecnología y Bioquímica, CINVESTAV, Irapuato, Guanajuato, México
Received 6 September 2004/ Accepted 19 October 2004
|
|
|---|
|
|
|---|
B. thuringiensis, Bacillus cereus, Bacillus anthracis, and Bacillus mycoides, as well as the recently described Bacillus pseudomycoides (41) and Bacillus weihenstephanensis (33), constitute the so-called B. cereus group. Several authors (4, 8, 19) have suggested that these species should constitute only one species, due to their high genetic similarity. From these species, B. thuringiensis is the most diverse, and its strains have been classified in 84 serovars (serovarieties) (32), including the recently described serovar jordanica (H serotype 71 [H71]) (29). Serotyping is still the most widely accepted subspecific classification technique for varieties of B. thuringiensis, even if strains from the same serovar do not necessarily share the biochemical, genetic, or toxicological attributes (3).
While some serovars, such as serovar israelensis (H14), include strains with practically the same attributes (2), other serovars include strains with a wide diversity of features. This is the case of serovar morrisoni (H8a,8a), which includes some strains with toxicity toward mosquito larvae (44), others toward coleopteran larvae (22), and some others toward lepidopteran larvae (13). On the other hand, strains from different serovars may show high biochemical, genetic, and toxicological similarity, such as strains IMR 81-1 (serovar malaysiensis), 11S2-1 (serovar canadensis), B 175 (serovar thompsoni), K6 (self agglutinated), and B 51(self agglutinated), highly similar to serovar israelensis (50). Additionally, serotyping is useless for nonmotile strains as well as the so-called self-agglutinated strains, besides the agglutination found in some B. cereus strains with H antigens (32, 42).
Alternative typing methods for B. thuringiensis strains have been tested, mostly based on molecular techniques, such as arbitrary primer-PCR technology (7, 18), ribosomal DNA restriction fragment length polymorphism (RFLP) (1, 48), and amplified fragment length polymorphism (AFLP) (45), among others (39, 58), most of them using a limited number of strains. Diversity of rRNA intergenic spacer sequences of 31 strains proved insufficient to discriminate between isolates (97 to 99% similarity) (6). On the other hand, ribotyping (16S, 23S, and 5S rRNA gene RFLP) of 80 serovars of B. thuringiensis showed a great diversity of patterns (27, 28), similar to the diversity found with fluorescent AFLP, when 34 B. thuringiensis serovars were analyzed along with strains of B. cereus and B. anthracis (21).
Repetitive extragenic palindromic sequence-based PCR analysis (Rep-PCR) is a DNA fingerprinting technique originally based on the design of PCR primers from Rep sequences found in the Escherichia coli and Salmonella typhimurium genomes (56). Amplicons obtained from contiguous Rep sequences generate distinctive electrophoretic patterns among different strains. Similar approaches use other repetitive sequences, such as the so-called ERIC and BOX sequences, developed for E. coli and S. typhimurium (24) and for Streptococcus pneumoniae, respectively (37). Rep-PCR fingerprint analysis of strains has proved to be simple, fast, and reproducible in a great variety of organisms (14, 36). However, this technique has been applied to organisms with little (if any) relationship with enterobacteria, that is, organisms with no homology whatsoever with the Rep sequences of E. coli, including some eukaryotic organisms (15, 38), which may indicate that these Rep-PCR analyses are arbitrary primer-PCR analyses, in those cases. This is the case of the Rep-PCR analysis of 28 B. thuringiensis serovars, using primers from the E. coli Rep sequence (46). We know now that this sequence is not found in the B. cereus group genomes.
This report presents the B. cereus (Bc)-Rep-PCR analysis of 125 B. thuringiensis strains, including 83 serovars, two biovars, and 26 native isolates, with primers specifically designed from a 26-bp Rep sequence found in the B. cereus group genomes.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. B. thuringiensis type strains from the IEBC, Institut Pasteur, Paris, France subjected to Bc-Rep-PCR fingerprinting
|
|
View this table: [in a new window] |
TABLE 2. Three groups of native B. thuringiensis isolates from the CINVESTAV-Irapuato stock collection subjected to Bc-Rep-PCR fingerprinting
|
Search for REPs in the B. cereus genome and primer design.
Due to the availability of the first B. cereus genome in 2003 (http://ergo.integratedgenomics.com/B_cereus.html), REP sequences were searched in this genome to design specific REP primers for the B. cereus group. All the extragenic sequences in the genome were analyzed with scripts written in Perl (http://www.perl.org/). Short REP sequences were combined to obtain larger ones until a highly conserved 26-bp sequence was found, showing the highest repeatability, in terms of both the number of repeats within the genome and the homology between the repeats. Its presence within the recently reported B. cereus group genomes was corroborated by searching the sequence in the European Molecular Biology Laboratory-European Bioinformatics Institute (EMBL-EBI) (http://www.ebi.ac.uk/fasta33/) and the National Center for Biotechnology Information (NCBI) (http://www.ncbi.nlm.nih.gov/) data banks. Direct and reverse primers were designed according to this sequence to amplify inter-REP regions.
Rep-PCR amplification conditions.
PCR mixtures were prepared as follows: 100 ng of template DNA, 300 ng of each primer, 5 mM MgCl2, 200 µM deoxynucleoside triphosphate mixture, and 2.5 U of Taq DNA polymerase (Invitrogen) to a 25-µl final volume. PCR amplifications were performed under the following conditions: an initial denaturation of 5 min at 94°C, followed by 34 cycles each of denaturation at 94°C for 1 min, annealing at 42°C for 1 min, and polymerization at 72°C for 1.5 min. Amplification was finished with an extension step at 72°C for 7 min. All PCR amplifications were performed with a Perkin-Elmer GeneAmp PCR System 2400. Amplified samples were kept at 20°C until electrophoretic analysis was performed.
Electrophoretic analysis.
Bc-Rep-PCR patterns were visualized by agarose gel electrophoresis. Aliquots of 10 µl each of the amplification products were loaded onto 1.2% agarose slabs (11 by 14 cm) and run in TAE (40 mM Tris-acetate, 1 mM EDTA) buffer at 2 V/cm during 5 h. Slabs were stained with 0.4 µg of ethidium bromide/ml and documented with a Gel Doc 2000 gel system (Bio-Rad). Molecular weight analysis of patterns was performed with the Quantity One version 4.2.1 software (Bio-Rad), with the 1-kb DNA ladder (Invitrogen) as a molecular weight marker.
Analysis of Bc-Rep-PCR patterns.
Polymorphic bands from all the Rep-PCR patterns were individually identified by their specific migration rates in the electrophoretic analyses. Once bands were properly and distinctively identified, binary (0/1) matrices were constructed to compare the patterns. Jaccard's similarity coefficients were generated by the SIMQUAL subroutine from the NTSYS-pc 2.02j (Applied Biostatistics, Inc.) package. Cluster analyses along with their corresponding dendrograms were generated by the unweighted-pair group method using average linkages (UPGMA), with the SAHN and TREE subroutines from the NTSYS-pc package.
|
|
|---|
|
View this table: [in a new window] |
TABLE 3. Frequency of the Bc-REP and designed primer sequences in five reported genomes of the B. cereus groupa
|
![]() View larger version (99K): [in a new window] |
FIG. 1. Bc-Rep-PCR fingerprint patterns of the B. cereus group strains. Lane 1, B. thuringiensis serovar israelensis; lane 2, B. thuringiensis LBIT-13; lane 3, B. cereus DSM31; lane 4, B. cereus CER81; lane 5, B. cereus CER 183; lane 6, B. mycoides IP-M 001; lane 7, B. anthracis 7702; lane 8, B. subtilis IP-S 001; lane 9, negative control; lane M, molecular weight marker (1-kb DNA ladder; Invitrogen).
|
![]() View larger version (46K): [in a new window] |
FIG. 2. Bc-Rep-PCR fingerprinting (A) and schematic representation (B) of 11 B. thuringiensis type strains. Lane 1, serovar amagiensis; lane 2, serovar israelensis; lane 3, serovar jinghongiensis; lane 4, serovar sumoyoshiensis; lane 5, serovar coreanensis; lane 6, serovar pakistani; lane 7, serovar konkukian; lane 8, serovar guiyangiensis; lane 9, serovar kurstaki; lane 10, serovar vazensis; lane 11, serovar brasiliensis; lane PB, polymorphic banding; lane M, molecular weight marker (1-kb DNA ladder; Invitrogen).
|
![]() View larger version (64K): [in a new window] |
FIG. 3. Bc-Rep-PCR fingerprinting of 17 different type strains of B. thuringiensis. (A) Lane 1, serovar sumiyoshiensis; lane 2, serovar fukuokaensis; lane 3, serovar dakota; lane 4, serovar sotto; lane 5, biovar dendrolimus; lane 6, serovar leesis; lane 7, serovar amagiensis; lane 8, serovar seoulensis. (B) Lane 1, serovar pakistani; lane 2, serovar alesti; lane 3, serovar kyushuensis; lane 4, serovar galleriae; lane 5, serovar aizawai; lane 6, serovar kurstaki; lane 7, serovar entomocidus; lane 8, biovar subtoxicus. Lanes M, molecular weight marker (1-kb DNA ladder; Invitrogen).
|
![]() View larger version (26K): [in a new window] |
FIG. 4. Dendrogram estimated from the Bc-Rep-PCR patterns obtained from 83 type serovars, two biovars, three isolates of B. thuringiensis, and five B. cereus group strains, using the Jaccard coefficient and UPGMA.
|
The mosquitocidal strains showed a Bc-Rep-PCR pattern identical to that of the mosquitocidal IPS-82 standard, including not only the serovar israelensis (H14) strains, but also one serovar kenyae strain (H4a,4c; strain LBIT-52) and one serovar entomocidus strain (H6; strain LBIT-58), as well as two self-agglutinated strains (LBIT-396 and LBIT-388) (Fig. 5). Additionally, another group of mosquitocidal strains, highly related to the first one, was represented by the serovar morrisoni PG-14 strain, along with strains LBIT-93 (serovar morrisoni; H8a,8b), LBIT-94 (entomocidus; H6), and LBIT-426 (not serotyped) (Fig. 5). Accordingly, other nonmosquitocidal strains, tested and characterized earlier (50), showed Bc-Rep-PCR patterns identical to that of israelensis (Fig. 6).
![]() View larger version (45K): [in a new window] |
FIG. 5. Bc-Rep-PCR fingerprint patterns of native mosquitocidal B. thuringiensis isolates. (A) Lane 1, serovar israelensis IPS-82; lane 2, LBIT-52 (H4a,4b); lane 3 LBIT-58 (H6); lane 4, LBIT-153 (H14); lane 5, LBIT-163 (not serotyped); lane 6, LBIT-201 (not serotyped). (B) Lane 1, LBIT-62 (not serotyped); lane 2, LBIT-388 (self agglutinated); lane 3, LBIT-393 (H 14); lane 4, LBIT-396 (self agglutinated); lane 5, LBIT-432 (not serotyped); lane 6, morrisoni PG14 (H8a,8b); lane 7, LBIT-93 (H8a,8b); lane 8, LBIT-94 (H6), lane 9, LBIT-426 (not serotyped). Lanes M, molecular weight marker (1-kb DNA ladder; Invitrogen); Kb, band size; PB, polymorphic band.
|
![]() View larger version (69K): [in a new window] |
FIG. 6. Bc-Rep-PCR fingerprint patterns of mosquitocidal B. thuringiensis strains. Lane 1, serovar morrisoni PG14 (H8a,8b); lane 2, serovar israelensis (H14); lane 3, serovar malaysiensis (H 36); lane 4, serovar canadensis 11S2.1 (H5a,5c); lane 5, serovar thompsoni B175 (H12); lane 6, K6 (self agglutinated); lane 7, B51 (self agglutinated); lane 8, serovar medellin 163.131 (H30); lane M, molecular weight marker (1-kb DNA ladder; Invitrogen).
|
![]() View larger version (52K): [in a new window] |
FIG. 7. Bc-Rep-PCR fingerprint patterns of ß-exotoxin-producing (A) and coleopteran-active (B) B. thuringiensis isolates. (A) Lane 1, serovar thuringiensis (T01 001); lane 2, LBIT-63 (H1); lane 3, LBIT-279 (H1); lane 4, LBIT-299 (H1); lane 5, LBIT-301 (H); lane 6, LBIT-398 (self agglutinated). (B) Lane 1, strain tenebrionis (T08 017); lane 2, strain san diego (H8a,8b); lane 3 LBIT-18 (H8a,8b); lane 4, LBIT-24 (H a,8b); lane 5, LBIT-73 (H8a,8b); lane 6, LBIT-74 (self agglutinated); lane 7, LBIT-196 (nonserotypable); lane 8, LBIT-358 (H8a,8b); lane 9, LBIT-419 (H8a,8b). Lanes M, molecular weight marker (1-kb DNA ladder; Invitrogen); Kb, band size; PB, polymorphic band.
|
|
|
|---|
Rep-PCR has been widely used on a variety of bacterial (and nonbacterial) species (5, 9, 15, 17, 38, 54) to characterize and identify strains. It has also been used for strains within the genus Bacillus (including B. thuringiensis), but based on the use of REPs found in other unrelated bacteria, such as the streptococcal BOX (10), the enterobacterial REP (35, 46), and the enterobacterial ERIC (35). In fact, a quick search for all these sequences in the six B. cereus group genomes reported showed no significant matches, indicating that actual BOX, REP, or ERIC analyses of these strains may be uncertain and should be reviewed. However, partial homology of enterobacterial REPs was found in Bacillus sporothermodurans, which allowed a real Rep-PCR analysis of the strains (20).
The presence of REP sequences in prokaryotes is common (36) and has been used for the design of species- or group-specific primers. That is the case of a 26-bp REP found in Neisseria spp., which allowed the design of specific primers for the analysis of N. gonorrhoeae and N. meningitides strains (57). Based on this approach, we looked for and found a 26-base REP common in the six B. cereus group reported genomes, which also include B. anthracis and B. thuringiensis. This REP (Bc-REP) allowed the design of two specific primers for the B. cereus group and proved their applicability by amplifying discrete and reproducible patterns in B. cereus, B. thuringiensis, B. anthracis, and B. mycoides strains. Their specificity was corroborated when a B. subtilis strain showed only faint bands and an undefined pattern, which may be caused to a partial homology with the Bc-REP.
The strong relationship between these species has been corroborated before, either by DNA hybridization (53), ribotyping (48), AFLP (21, 49, 55), or BOX-PCR (10) analyses of a number of strains. In all cases, strains from all these species intermingle within the same dendrogram, with, in general, the B. anthracis strains being the most homogeneous, the B. cereus and B. thuringiensis strains being the most diverse, and the B. mycoides strains being the least related to the rest. In our study, only three B. cereus, one B. anthracis, and one B. mycoides strains were analyzed; although all five strains intermingled in the same dendrogram with the B. thuringiensis strains, more strains from the other species were required to corroborate the same trend when Bc-REP-PCR analysis was used.
The main purpose of this report is the Rep-PCR characterization of the B. thuringiensis type strains, using specific primers for the B. cereus group. All the type strains were included, except for the most recent serotype described last year (29). Bc-REP-PCR fingerprinting of the type strains showed that practically all the serotypes displayed a distinct pattern. It also shows the putative phylogenetic relationship between the 83 serotypes and two biovars included in the analysis. Only a few strains showed identical patterns, such as serovars entomocidus and its biovar subtoxicus, both isolated by Heimpel (26) in Canada. Similar results were obtained by Phucharoen et al. (46) and Brousseau et al. (7), but they differ from the results obtained by Priest et al. (48) and Joung and Côté (27). On the other hand, Bc-REP-PCR patterns from the serovar sotto and its biovar dendrolimus slightly differ from each other, with the serovar sotto pattern identical to that of serovar dakota and the serovar leesis pattern the closest to that of biovar dendrolimus. Interestingly, ribotyping of these strains (27) indicated that while biovar dendrolimus and serovars leesis and dakota were phylogenetically related, serovar sotto was located in a separate group. Other serovars that shared the same Bc-REP-PCR pattern were serovars amagiensis and seoulensis. This is in agreement with a previous random amplified polymorphic DNA analysis of these strains (16); ribotyping also connects both strains in the same group (27).
Other serovars highly related by their Bc-REP-PCR pattern were serovars galleriae, aizawai, and kurstaki, which agrees with previously reported DNA hybridization and RFLP analyses (40, 48). These results may indicate that these associated and highly common serovars form a tight phylogenetically related group, whose segregation should be reviewed. Likewise, such segregation should be reviewed for the highly related serovars sumiyoshiensis and fukuokaensis, indiana and thompsoni; amagiensis, seoulensis, and kyushuensis; thuringiensis and sooncheon; azorensis and vazensis; and monterrey and oswaldocruzi. All these serovars appear closely related among each other in the phylogenetic dendrogram generated by the Bc-REP-PCR fingerprinting and by previously reported ribotyping (27). Other serovars highly related by their Bc-REP-PCR pattern, such as serovars silo and ostriniae, palmanyolensis and darmstadiensis, japonensis and kenyae, and colmeri and mexicanensis show less of a relationship by ribotyping analysis (27).
Serovars israelensis and malaysiensis also share the same Bc-REP-PCR pattern, which is in agreement with their high relationship (mosquitocidal specificity, cry gene content, and crystal morphology) (50); however, ribotyping is unable to recognize such a relationship and locate both strains in separate groups (27). Interestingly, other mosquitocidal strains with attributes practically identical to those of serovar israelensis such as K6 (AAT028), B51 (AATO21), canadensis 11S2.1 (T05A030), and serovar thompsoni B175 (T12007) (50) also display Bc-REP-PCR patterns identical to those of serovars israelensis and malaysiensis. On the other hand, the mosquitocidal serovar medellin, which has attributes different from those of serovar israelensis (43, 50), also shows a very different Bc-REP-PCR pattern. These results may indicate that the genomic relationship between the B. thuringiensis strains is not necessarily defined only by their toxic specificity, but by a series of attributes, such as cry gene content, crystal morphology, and plasmid pattern. This was also corroborated when native mosquitocidal isolates, highly related to serovar israelensis, showed identical Bc-REP-PCR patterns, even when some of these isolates were self agglutinated or belonged to a serotype different from that of israelensis. These results not only corroborate the reported genomic homogeneity of the serovar israelensis strains (2) but also imply that the same homogeneity occurs in other non-israelensis strains, as long as they share other attributes.
Another mosquitocidal strain, serovar morrisoni PG14, is known to share only some of the serovar israelensis characteristics (toxic specificity, some cry gene content, and crystal morphology) (25, 50); accordingly, its Bc-REP-PCR pattern is similar but not identical to that of serovar israelensis. Interestingly, its pattern is significantly different from that of the serovar morrisoni type strain. This serovar also includes the coleopteran-active strains tenebrionis and san diego, previously reported to be the same strain (30); however, contrary to the PG14 strain, the identical Bc-Rep-PCR pattern of both strains is very similar to that of the serovar morrisoni type strain. Also, similar to the results obtained with the mosquitocidal strains, the native coleopteran-active strains share the same Bc-Rep-PCR pattern as the serovar tenebrionis strain, most of them serotyped as serovar morrisoni but also including two nonserotypable strains. Genomic homogeneity may also occur in this group, similar to that observed with the serovar israelensis group. Interestingly, the same homogeneity was found in the group of native strains that produce ß-exotoxin and belong to serovar thuringiensis (H1) (except for the self-agglutinated LBIT-398 strain). These results may indicate that Bc-Rep-PCR fingerprinting of B. thuringiensis strains is useful not only to differentiate between serovars, but also to properly identify the nonserotypable strains and, most of all, to recognize more accurately the evolutionary relationship between strains, to whichever serovar they belong.
B. thuringiensis constitutes a genetically diverse species; the great number of strains known today may form distinctive groups, according to their phenetic and genetic traits. Serotyping has been a useful tool to try to discriminate those groups since its establishment in 1962 (12); however, as strains mounted up, shortcomings started to appear in the technique. Molecular tools have been developed in recent years, trying to offer a new typing alternative for B. thuringiensis strains and to recognize the actual phylogenetic relationships between subspecific groups. Bc-Rep-PCR offers a new tool to identify these groups, based on the use of specific primers designed from a REP sequence found in the B. cereus group. The potential of this technique was tested in this work and proved to be sensitive, specific, reproducible, and fast; it may become a standardized characterization procedure. It may also help in the establishment of a new subspecies-level classification of B. thuringiensis.
This work was partially supported by grant 35320-B, CONACYT, Mexico.
|
|
|---|
-endotoxin electron microscopy. Gen. Mol. Biol. 25:225-229.
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»