Previous Article | Next Article ![]()
Applied and Environmental Microbiology, December 2005, p. 8107-8114, Vol. 71, No. 12
0099-2240/05/$08.00+0 doi:10.1128/AEM.71.12.8107-8114.2005
Copyright © 2005, American Society for Microbiology. All Rights Reserved.
Department of Biochemistry and Molecular Biology, Oswaldo Cruz Institute, FIOCRUZ, Rio de Janeiro, Brazil,1 Department of Biological Sciences, National School of Public Science, FIOCRUZ, Rio de Janeiro, Brazil,2 Department of Bacteriology, Oswaldo Cruz Institute, FIOCRUZ, Rio de Janeiro, Brazil,3 Department of Biology, Université de Moncton, Moncton, Canada,4 New Brunswick Research and Productivity Council, Fredericton, Canada5
Received 6 July 2005/ Accepted 2 September 2005
|
|
|---|
|
|
|---|
Some strains of B. cereus cause food poisoning and other infections (29). Two principal types of food poisoning, emetic and diarrheal, have been described (38). The emetic type is induced by a small, cyclic, heat-stable peptide, cereulide, which causes vomiting a few hours after ingestion (2, 31). Evidence for a clonal population structure of cereulide-producing emetic B. cereus has recently been presented (14). Diarrheal food poisoning has been attributed to the production of enterotoxins during vegetative growth of B. cereus within the small intestine (18, 31). Five enterotoxins have been reported to play a role in diarrheal food poisoning: three single gene products encoded by the sequences entFM, cytK, and bceT and two multicomponent protein complexes, namely, hemolysin BL (HBL) and nonhemolytic enterotoxin (NHE) (7, 18, 19).
It has been reported that both B. cereus and B. thuringiensis are capable of growth within the hemocoel of various insect species and may provoke septicemia, although the former species is not generally considered to be a true pathogen of insects (37). The entomopathogenic characteristics of B. thuringiensis result from the production of toxins (Cry, Cyt, and Vip3A) encoded by genes carried on plasmids (37, 40). If these plasmids are lost, B. thuringiensis can no longer be clearly distinguished from B. cereus. Expression of the Cry and Cyt crystal proteins, which are synthesized when the bacterium enters the stationary phase or the sporulation phase, has never been reported for any B. cereus isolate (37). On the other hand, little is known concerning the interspecies distribution of the plasmid carrying the vip3A sequence, which is expressed during vegetative growth (15). The limited data produced to date from PCR- and colony hybridization-based investigations have shown that the vip3A sequence is present in a low number (15 to 30%) of B. thuringiensis isolates (15, 20, 35).
Analysis of repetitive-element sequence polymorphism has been suggested as a means of differentiating between different members of the B. cereus complex. Indeed, PCR based analysis of a variable-number tandem repeat (VNTR) sequence in the vrrA gene was demonstrated to allow the differentiation of a limited number of B. cereus, B. thuringiensis, B. mycoides, and B. anthracis isolates simply by comparing the sizes of amplicons on agarose gels (4, 27).
The first objective of this study was to assess the distribution of genes encoding putative virulence factors, specifically, the Vip3a, BceT, PclA, Sph, NheB, and NheC components of the NHE complex and HblA of the HBL complex, in a panel of B. cereus and B. thuringiensis strains comprised mainly of Brazilian isolates. A second objective was to examine the variable region of the vrrA sequence for fragment length polymorphisms, which might prove useful for the specific identification of B. cereus and B. thuringiensis. The results of this study are discussed in the context of the ongoing debate over the relationship between these two taxa.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Distributions of virulence genes, VNTR types, and origins of the Bacillus species examined
|
|
View this table: [in a new window] |
TABLE 2. Primers and amplification conditions used in this study
|
|
|
|---|
The rationale for selecting the particular sequences examined was based first on the fact that the presence of these genes in a variety of strains of both taxa from other geographical locations had been established previously and second on our hope to extend upon previous data concerning the distribution of the vip3A sequence within these taxa. This issue required clarification, given that data from colony screening experiments with unidentified Bacillus isolates (15) suggested that the vip3A sequence could be present in some B. cereus strains.
As reported in previous PCR-based screening studies, the enterotoxin genes examined herein were widespread throughout both species and were present in strains isolated from food, clinical material, insects, and soil environments (Table 1). The presence of the HBL operon, examined via amplification of the hblA gene, encoding the B component of the tripartite enterotoxin, was detected at a level of 90% in both the total B. thuringiensis population and the Brazilian subset of strains (Table 1). The common occurrence of the HBL operon in B. thuringiensis had been demonstrated previously via PCR-based detection of the hblCD genes in 87% (65 of 74) of the isolates examined by Gaviria Rivera et al. (17) and in 92% (38 of 41) of the isolates investigated by Hansen and Hendriksen (21). Unfortunately, no information was provided as to the geographical origin of the strains used in either of those studies.
Amplicons of the expected size (883 bp) were detected for 86% of the total B. cereus population examined and for 90% of the Brazilian isolates (Table 1). These values were markedly higher than the levels recorded, using the same primers, by Prüb et al. (33) or Mäntynen and Lindström (30), who reported detection frequencies of 43% (10 of 23 isolates) and 52% (26 of 50 isolates), respectively. It is pertinent to note that in the aforementioned studies, the bulk of the strains examined were of European origin, indicating that geographic origin may influence the degree to which the hblA sequence is maintained among populations of B. cereus. It is of interest that in spite of previous studies showing a relatively limited distribution of sequences encoding the HBL complex protein among B. cereus isolates (21, 30, 31), this molecule has been classified as a hemolysin based on its ability to lyse mammalian erythrocytes. Insects do not have red blood cells, and the near universal carriage of these genes by the insect pathogen B. thuringiensis and by most of the B. cereus strains isolated from insects (18/19) in this study suggests that the actual biological function of this molecule is linked to the establishment of infections in insects, rather than to the lysis of red blood cells in potential mammalian hosts.
The sequence denoted bceT was reported to encode a single-component toxin (BceT) which exhibited Vero cell cytotoxicity and to which has also been attributed a role in the diarrheal syndrome (3). More recently, it has been suggested that the bceT sequence actually represents a cloning artifact, formed by the ligation of four individual restriction fragments derived from B. cereus genomic DNA (22). An additional study indicated that the BceT toxin was most likely an experimental artifact that probably could not contribute to food poisoning (10). Nevertheless, other workers have demonstrated that it is possible to detect at least part of the original published sequence in both B. cereus and B. thuringiensis, and for that reason, we decided to include that region of the sequence as a target in our study using previously described primers (17). In common with our data concerning the hblA sequence, the frequency of distribution of the bceT gene was found to be broadly similar throughout both taxa (Table 1). Specifically, we detected the expected 593-bp amplicon in 86% of the entire B. thuringiensis population, in 84% of the Brazilian isolates, in 78% of the total B. cereus population, and in 69% of the Brazilian B. cereus strains. Interestingly, a higher level of detection (95%) was observed among the non-Brazilian B. cereus strains, which had been isolated predominantly from insects (Table 1). Thus, in common with previous studies (17, 21), we would conclude that this sequence is broadly distributed throughout both taxa.
In contrast, the nheBC sequence was detected at a markedly lower frequency among the B. cereus strains analyzed (60% in the total test population and 50% in Brazilian strains) than among the B. thuringiensis samples (81% in both the total test and Brazilian populations) (Table 1). The level of carriage observed in our panel of B. thuringiensis isolates is consistent with the studies of Gaviria et al. (17) and Hansen and Hendriksen (21), who reported frequencies of 100% and 80%, respectively, within their smaller strain panels. Conversely, the observation that only 60% of our B. cereus isolates possessed these sequences seemed at odds with a previous investigation which demonstrated the production of the NheA protein in a substantial proportion (92%) of the 194 strains examined (36). Yet, the level seen in our study is similar to that reported by Hansen and Hendriksen (21), who detected the genes nheA, nheB, and nheC in only 59% of the 41 strains used in their study. A more profound analysis of the distribution of these sequences among our strain panel revealed that the apparently limited presence of these genes was strongly influenced by the extremely low level of detection (11%) among the Brazilian strains isolated from foodstuffs (Table 1). Moreover, the levels of detection noted among the environmental and clinical strains, both Brazilian and non-Brazilian, were in the range of 79% (Table 1), a value that falls between the levels reported in the previous studies. Taken together, these data suggested that despite having been defined as a toxin involved in diarrheal food poisoning in B. cereus, the nonhemolytic enterotoxin is probably not essential for virulence. On the other hand, the widespread distribution of these sequences among B. thuringiensis strains of diverse geographical and environmental origins indicated that this molecule fulfils an important biological function for this bacterium, as speculated for the BL hemolysin.
In comparison to the other sequences reported to occur in isolates of both taxa, the distribution among B. thuringiensis isolates of the plc gene, which encodes phosphatildylcholine-specific phospholipase C, has received relatively little attention. Data from the only other study to have investigated this topic showed the presence of this sequence in 100% of the 16 strains examined (20). Employing the same primer set and amplification conditions as used previously, we observed the presence of the plc gene in 64% of our B. thuringiensis strains. This sequence was detected more frequently in the non-Brazilian isolates (15/19; 78.9%) than in the Brazilian subset (21/38; 55.2%). Interestingly, the level in the Brazilian isolates is substantially lower than that reported previously, which questions the notion that this sequence has a universal distribution within this taxon. In the case of B. cereus, we detected the plc sequence in only two isolates, both of which were Brazilian. The first of these (LFB287), was isolated from contaminated infant milk formula, while the second strain (LFB1181) was recovered from an outbreak of avian septicemia at a zoo in São Paulo (Table 1). This result was somewhat unexpected, since this gene was detected in 100% (3/3) of the B. cereus isolates studied previously (20). It is possible that the PCR-negative strains lacked the plc sequence. Alternatively, the lack of amplification may have arisen from sequence variation in the region selected for primer annealing. Further study, using primers designed to amplify alternative regions of this sequence, will be necessary to clarify this issue.
The sph sequence, encoding sphingomyelinase, was present in all strains analyzed by PCR. This observation supported the detection of this sequence as a good indicator of the presence of B. cereus group strains, as previously proposed (25), and suggested that this molecule fulfils an essential biological function for members of both taxa. It is of considerable interest that the expression of a number of extracellular virulence factors identified in both B. thuringiensis and B. cereus, including sphingomyelinase, HBL, NHE, and phosphatildylcholine-specific phospholipase C, has been shown to be pleiotropically regulated by the PlcR protein (1). The sequences encoding these putative virulence factors are dispersed throughout the chromosome, rather than being concentrated within a pathogenicity island. This has led to the suggestion that these species are derived from a common ancestral bacterium containing the plcR regulon and associated regulated sequences (1), which further supports the notion that these bacteria represent variants of the same species.
Using PCR-based detection methods, it was demonstrated that the vip3A gene showed a limited distribution (23 to 37%) among B. thuringiensis isolates of undefined geographical origin (20, 35). However, a more recent study reported a significantly greater distribution level of 75% among a population of 40 Chinese strains (40). Using our predominantly Brazilian strain panel, we detected the vip3A sequence in only 7% of the total test population (Table 1). All four isolates were Brazilian and belonged to the Morrisoni (LFB679, -780, and -781) and Kurstaki (LFB263) serotypes, both of which have been reported to present toxicity towards lepidopteran insects (37). However, a correlation between serotype and the presence of this sequence was not established, given that two additional isolates from each serotype tested negative for the presence of the gene. In the case of B. cereus, the vip3A sequence was not detected in any of the 63 strains examined in our study. Given the absence of any other published studies examining the presence of this sequence in B. cereus, we are unable to state that our data represent conclusive evidence that this sequence is found only in B. thuringiensis. However, we do not consider it unreasonable to make such a suggestion, in the hope that it may stimulate additional research on this topic.
It has been reported recently that the vip3A sequence is carried on a 31.8-kb plasmid (40). Conversely, evidence indicating that the gene is chromosomally located has been presented (16), although that conclusion was based on the examination of a single isolate. The widely divergent values recorded for carriage of this gene among strains of B. thuringiensis may indicate that the acquisition of this gene, possibly via horizontal gene transfer from other soil microbes, was a relatively recent event in the evolution of B. thuringiensis. Furthermore, the absence of this sequence, as well as those encoding for insecticidal crystal proteins, in B. cereus indicates that the inability to effectively separate these taxa genotypically, using increasingly complex molecular techniques, may be a reflection of our inability to identify the genetic basis for both the acquisition and maintenance of the plasmids carrying such sequences.
The presence of sequence polymorphism in the vrrA gene, resulting from variability in the number of tandem repeats, had previously been reported to allow the separation of different members of the B. cereus complex and even allowed differentiation between individual isolates of B. anthracis and B. cereus (4, 27). In the present study, we amplified this region of the genome and observed a total of five differently sized amplicons (Fig. 1), designated VNTR types. The most commonly encountered VNTR type in both B. cereus and B. thuringiensis was type 3 (150 bp), which was detected in approximately 45% of strains from both taxa (Table 1). Type 2 (135 bp) was observed in 38% of the B. thuringiensis isolates and at a level of 22% among the B. cereus strains (Table 1). The three other VNTR types were detected at far lower levels among the strain panel, with type 1 being the least frequently detected (Table 1). Interestingly, an amplicon of the same size (90 bp) was also detected in the single B. mycoides isolate included in the study (Table 1). No correlation between VTNR type and the origin of any strain could be demonstrated, and the suitability of variation in this region for separating the two taxa was clearly not established. Indeed, our finding that the same two VNTR types were more or less equally distributed between the two taxa could be interpreted as further evidence in support of the hypothesis that these two taxa represent the same species.
![]() View larger version (101K): [in a new window] |
FIG. 1. Agarose gel electrophoresis of representative VNTR amplicons in DNA extracted from Bacillus strains. Lane M, 50-bp DNA ladder (Sigma) (molecular sizes are indicated on the left); lane A, VNTR profile 5 (LFB596); lane B, VNTR profile 1 (LFB770); lane C, VNTR profile 2 (LFB856); lane D, profile 3 (LFB677); lane E, VNTR profile 4 (LFB1126).
|
Second, our data extend the observations of various other research groups who have assessed the distribution of the putative virulence factors among these two closely related taxa. Numerous studies have attributed roles to these products from the viewpoint of their action upon animal cells or whole animals. Intriguingly, almost all studies performed to date have observed that the genes encoding the so-called mammalian virulence factors occur more frequently in B. thuringiensis, which is considered to be nonpathogenic towards mammals. Taken as a whole, these data suggest that the majority of B. cereus isolates could be considered to be nonpathogenic and that those strains which can cause infections do so because theyhave retained the majority of the sequences found in B. thuringiensis. However, the dearth of documented cases of mammalian infections caused by B. thuringiensis (26) clearly demonstrates that possession of these sequences does not necessarily correlate with virulence.
Finally, our examination of polymorphisms due to VNTR in the variable region of the vrrA sequence clearly demonstrated the close taxonomic proximity of B. cereus and B. thuringiensis, based on the theory that analysis of VNTR provides functional and evolutionary information concerning genetic relationships in microbial species (39). This observation will undoubtedly contribute to the continuing debate over whether these bacteria should actually be regarded as variants of the same species rather than as two distinct taxa.
We are grateful to Alberto Dávila (DBBM/IOC/FIOCRUZ) for hishelp in designing the primers for the amplification of the vip3A sequence.
|
|
|---|
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»