Previous Article | Next Article ![]()
Applied and Environmental Microbiology, September 2002, p. 4216-4224, Vol. 68, No. 9
0099-2240/02/$04.00+0 DOI: 10.1128/AEM.68.9.4216-4224.2002
Copyright © 2002, American Society for Microbiology. All Rights Reserved.
Nestlé Research Center, Lausanne, Switzerland,1 Centre for Agricultural ResearchGhent, Department of Animal Product Quality, B-9090 Melle, Belgium2
Received 6 December 2001/ Accepted 30 May 2002
|
|
|---|
|
|
|---|
According to its original description (20), an important characteristic of B. sporothermodurans is its ability to produce highly heat-resistant spores (HRS) that may survive sterilization (115 to 120°C for 15 to 20 min) or UHT treatment (135 to 142°C for a few seconds). Huemer et al. (16) found D values at 140°C that varied between 3.4 and 7.9 s and z values that varied between 13.1 and 14.2°C for spore preparations from the original stock culture, but they also observed a significant decrease in the heat resistance after multiple laboratory culture passages. Surviving spores of B. sporothermodurans can germinate and multiply in products to a maximal concentration of ca. 105 cells/ml. Although the vegetative cells are not pathogenic (9, 10) and do not cause significant visual or taste deviations, their presence in sterilized and UHT-treated products is considered undesirable, as such products do not meet the legal requirements established by the European Union (1).
For detection and identification of B. sporothermodurans in raw and consumer milk, a PCR-based method, now called HRS-PCR (22), was developed by Herman et al. (13). Several molecular methods have been used successfully to differentiate and characterize B. sporothermodurans strains, including PCR methods like random amplified polymorphic DNA analysis, repetitive extragenic palindromic (REP)-PCR, and 16S ribosomal DNA (rDNA) sequence analysis (12, 18, 20). All these techniques showed that the different B. sporothermodurans strains isolated so far from European UHT-treated and sterilized milk are phylogenetically very closely related, forming the so-called HRS clone, named after the initial description of this highly heat-resistant spore-forming organism (8). Recently, B. sporothermodurans strains have also been isolated from raw milk and from animal feed (4, 22, 24). The majority of these farm isolates reacted negatively in the HRS-PCR of Herman et al. (13) but could be assigned to B. sporothermodurans by a polyphasic approach and/or a new 16S rDNA-based specific PCR identification test (22).
The purpose of the present study was to evaluate the genetic diversity of B. sporothermodurans isolates by using two molecular typing techniques (ribotyping and REP-PCR) targeted at different genomic sites as recommended by different authors (14, 25). REP-PCR has been recognized to be highly discriminatory for the differentiation of many bacterial species, including Escherichia coli and B. sporothermodurans (5, 11, 12, 18). In addition, ribotyping has proved to be useful for differentiation at the subspecies level, as demonstrated for species like Pseudomonas aeruginosa, Vibrio cholerae, Listeria monocytogenes, and Clostridium botulinum (3, 7, 19, 23). Both DNA fingerprinting techniques were applied to B. sporothermodurans strains isolated from a wide range of geographic areas, including Asia, South America, and Europe, and from different sources, as outlined above. The resulting genetic relationships found among the strains should allow workers to determine the relative contributions of different possible contamination sources.
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. B. sporothermodurans and B. oleronius strains used in this study and their originsa
|
Automated ribotyping and data analysis.
All B. sporothermodurans isolates were processed at least two times by using the automated microbial characterization system (RiboPrinter; Qualicon, Inc., Wilmington, Del.) and restriction enzymes EcoRI and PvuII (Qualicon) as described elsewhere (2). All ribopatterns with similarity coefficients higher than 0.93 were considered identical by the RiboPrinter software and were grouped together based on the position and intensity of the bands to form a ribogroup (a set of isolates with indistinguishable ribotypes). Further refinement of the automated ribogrouping was performed by visual evaluation of closely related ribopatterns, which resulted in merger or separation of ribogroups.
REP-PCR.
Total genomic DNA from purified B. sporothermodurans strains was isolated by using the method of Pitcher et al. (21), as slightly modified by Heyndrickx et al. (15), and the DNA concentration was determined with a spectrophotometer. Total genomic DNA (25 ng) was used as the template in REP-PCR performed with primers REP1R-I (5' IIIICGICGICATCIGGC3') and REP2-I (5'ICGICTTATCIGGCCTAC3') (Isogen Bioscience bv), and the REP fragments were separated by denaturing polyacrylamide gel electrophoresis and silver stained as described by Herman et al. (12). All strains were analyzed in the same REP-PCR experiment and on the same gel to minimize possible variations in patterns caused by experimentation (12). A REP-PCR pattern could not be obtained for strain MB 359.
Banding pattern data analyses.
Digitized images of the gel obtained with the RiboPrinter were converted with the GelConvert program (Qualicon) and were analyzed by the unweighted pair group arithmetic (UPGMA) clustering algorithm computed by the GelCompar software (version 4.1; Applied Maths, Sint-Martens-Latem, Belgium) by using the Pearson product moment correlation coefficient with optimization of 1%. REP-PCR patterns on the silver-stained gels were scanned with a flat-bed scanner (Agfa SnapScan1236S; Agfa-Gevaert N.V., Mortsel, Belgium), and images were analyzed by UPGMA clustering computed by BioNumerics software (version 2.0; Applied Maths) by using the Pearson product moment correlation coefficient with optimization of 1%. Normalized REP-patterns were also visually classified in REP groups.
Combined clustering of the two ribopatterns and the REP-PCR profiles was performed by using the BioNumerics 2.0 software. To do this, the converted ribopatterns obtained with EcoRI or PvuII were introduced into the BioNumerics 2.0 program, and each pattern was linked with the strain database. Combination of the three experiment types was implemented in such a way that the typing methods used (REP-PCR, EcoRI riboprinting, and PvuII riboprinting) were weighted 2:1:1. For UPGMA clustering for this combination of experiment types, the same cluster analysis settings that were used for each experiment type separately were used (i.e., Pearson product moment correlation coefficient). Finally, multidimensional scaling of the combined cluster analysis allowed visually interpretable grouping of the strains in a three-dimensional plot.
|
|
|---|
Ribotyping.
Molecular characterization of B. sporothermodurans strains and the B. oleronius reference strain was performed by automated ribotyping by using the restriction enzymes PvuII and EcoRI. Sixteen distinct major ribogroups were determined with PvuII (ribogroups A1 to A16) (Table 1 and Fig. 1). The three largest ribogroups, ribogroups A1, A5, and A13, contained 21, 2, and 2 strains, respectively (Tables 1 and 2), whereas the other 13 ribogroups were each represented by a single strain. Based on the cluster analysis performed with the ribopatterns of 39 strains (including B. sporothermodurans type strain MB 581 and the reference strain of B. oleronius), a dendrogram was constructed (Fig. 1). Two main clusters, clusters P1 and P2, were discerned visually and by cluster analysis at similarity levels of 91 and 71%, respectively. These two clusters exhibited 58% similarity to each other. The majority of the UHT isolates (21 of 25 isolates) were found in cluster P1, which included one ribogroup. These UHT isolates appeared to be closely related to each other, as reflected by the high similarity values (91 to 99%). All these isolates produced a typical pattern consisting of at least seven conserved bands at ca. 7.2, 8.3, 9.3, 11.0, 13.0, 25.0, and 45.0 kbp. In cluster P2, which was defined at a lower similarity level (71%), 14 different ribogroups were found; there were typical bands at ca. 5.4, 7.1, 7.5, 8.5, 9.0, 11.0, and 25.0 kbp in the patterns of most of the strains, while 13.0- and 45.0-kbp fragments were missing compared to the patterns of cluster P1 strains. Remarkably, UHT strains MB 372, MB 373, and MB 374, all of which originated from Germany, had clearly distinct patterns and were members of cluster P2 along with all 13 farm isolates. Type strain MB 581 was located separately between clusters P1 and P2.
![]() View larger version (60K): [in a new window] |
FIG. 1. Dendrogram of 38 strains of B. sporothermodurans and one strain of B. oleronius obtained after restriction with PvuII. The strain designations correspond to those shown in Table 1. The scale bar indicates the percentage of similarity.
|
|
View this table: [in a new window] |
TABLE 2. Relationships between the different most important PvuII and EcoRI ribogroupsa
|
![]() View larger version (64K): [in a new window] |
FIG. 2. Dendrogram of 38 strains of B. sporothermodurans and one strain of B. oleronius obtained after restriction with EcoRI. The strain designations correspond to those shown in Table 1. The scale bar indicates the percentage of similarity.
|
B. oleronius MB 397 was ribotyped as an outlier. The PvuII- or EcoRI-generated ribopatterns of B. oleronius MB 397 were clearly different from those obtained with all the B. sporothermodurans strains (similarity levels, 40 and 57%, respectively).
REP-PCR.
Molecular characterization of B. sporothermodurans strains and the B. oleronius reference strain was performed by REP-PCR by using high-resolution separation of the bands by polyacrylamide gel electrophoresis and silver staining. Fourteen distinct major REP groups were visually determined (REP groups C1 to C14) (Table 1 and Fig. 3). The two largest REP groups, REP groups C1 and C2, contained 3 and 22 strains, respectively (Table 1), whereas the other 12 REP groups were each represented by a single strain. In comparison with B. oleronius, all B. sporothermodurans strains were characterized by a conserved major band at ca. 1,020 bp. Furthermore, strains belonging to REP groups C1 and C2, which contained all UHT isolates and only one farm isolate (MB 1505), were characterized by major conserved bands at ca. 875, 730, and 600 bp. Remarkably, the German UHT strains, MB 372 to MB 374, belonged to the same REP group, REP group C1, which could be differentiated from REP group C2, which contained the other UHT strains from different countries, by the absence of a major conserved band at ca. 850 bp in REP group C2. The largest REP group (REP group C2) could be further subdivided into four subgroups on the basis of minor polymorphisms (presence or absence of minor bands at various molecular weights). These four subgroups, subgroups C2a to C2d (Table 1 and Fig. 3), contained 8, 2, 2, and 10 strains, respectively.
![]() View larger version (93K): [in a new window] |
FIG. 3. Dendrogram of 37 strains of B. sporothermodurans and one strain of B. oleronius obtained by REP-PCR. The strain and REP group designations correspond to those shown in Table 1. The scale bar indicates the percentage of similarity as determined with the Pearson coefficient. The dark zones under the metric scale represent parts of the REP patterns omitted in the numerical analysis.
|
The REP pattern of B. oleronius reference strain MB 397 was clearly different from the main B. sporothermodurans pattern.
Combined analysis of ribotyping and REP-PCR patterns.
The overall or consensus genetic relatedness among the B. sporothermodurans strains was inferred from a combined numerical analysis of the ribopatterns obtained with PvuII and EcoRI and the REP-PCR patterns by performing a UPGMA cluster analysis and by three-dimensional scaling. In the cluster analysis (data not shown), 21 of the 24 UHT strains included clustered together at a minimal similarity level of 81%. In this cluster, 19 UHT strains clustered together at 86% similarity, with strain MB 582 and type strain MB 581 joining at 84 and 81%, respectively. Conversely, 10 of the 13 farm strains and the three German UHT strains (MB 372 to MB 374) produced a cluster at a lower level of similarity (70%), exhibiting only 67% similarity to the major cluster of UHT strains. The remaining three farm strains (MB 1316, MB 1495, and MB 1499) clustered together at 76% similarity and exhibited a lower level of similarity (64%) to all other B. sporothermodurans strains.
The results of the combined analysis of the data by three-dimensional scaling are shown in Fig. 4. In this nonhierarchical presentation of the relationships among the strains, essentially the same groups were obtained as in the hierarchical cluster analysis explained above. All UHT strains except the three German UHT strains formed a compact group, while all farm strains formed a very diffuse group clearly separated from the group of UHT strains. The only remarkable difference from the UPGMA cluster analysis was the grouping of the three German UHT isolates at the outer edge of the diffuse group of farm strains, showing a somewhat closer relationship to the other UHT isolates.
![]() ![]() View larger version (54K): [in a new window] |
FIG. 4. Visual three-dimensional representation of the combined clustering of the two ribotyping patterns (EcoRI and PvuII) and REP-PCR profiles of 37 strains of B. sporothermodurans and one strain of B. oleronius, obtained by multidimensional scaling of the cluster as explained in text. The strain designations correspond to those shown in Table1.
|
|
|
|---|
In a polyphasic typing approach, separate clustering analyses of PvuII and EcoRI ribopatterns and of REP-PCR patterns were largely consistent with each other and revealed the existence of two main clusters; one homogeneous group contained all (REP-PCR) or most (ribotyping) of the UHT isolates, and the second, more diverse group comprised the farm isolates (Fig. 1 to 3). The high level of genetic homology of most of the UHT isolates was further shown by a combined analysis of all molecular typing data in this study, both by cluster analysis and by three-dimensional scaling, which revealed a very compact cluster or group of isolates. The close genetic relationship of the UHT isolates suggests a clonal origin (HRS clone), which is particularly remarkable since B. sporothermodurans strains were isolated from UHT-treated and sterilized milk samples produced on three different continents.
The three German isolates were the only UHT strains whose genetic characteristics were quite different from those of the majority of the UHT strains. In a combined cluster analysis of all molecular typing data obtained in this study, the three German UHT strains clustered with the farm strains. A three-dimensional scaling analysis of all molecular typing data showed that these strains were at the border of the diffuse group of farm strains and directed to the compact group of UHT isolates. The latter observation and the REP-PCR and EcoRI ribotyping cluster analysis data suggest that the German UHT isolates have a remote genetic relationship with the HRS clone. This could also suggest that the extreme resistance of spores to sterilization temperatures is restricted to particular clones that have a possible common ancestor.
In contrast to the homogeneity found for the majority of the UHT isolates, the ribopatterns and REP patterns of the B. sporothermodurans strains isolated from animal feed (feed concentrate, silage, soy) and raw milk were much more diverse. Most of the ribogroups and all 12 REP types for the farm isolates were represented by a single strain (Table 1). Also, the two HRS-PCR-positive farm strains, which originated from feed concentrate and silage, produced patterns that were different from each other and from the patterns of the main group of UHT isolates (HRS clone) in the ribotyping analysis. Overall, it seems that there is no 100% concordance between a positive result in the HRS-PCR analysis and an HRS clone pattern determined by molecular typing.
Milk powder has been suggested to be a possible source of contamination of heat-treated dairy products (8). Since in some plants, UHT-treated or sterilized milk is prepared from imported milk powder, this practice could explain the spread of the same B. sporothermodurans HRS clone over different continents. Within one country, one can envisage a contamination route via raw milk that has been contaminated through animal feed at the farm level. Based on this assumption, one would expect to find B. sporothermodurans strains having similar ribopatterns or REP-PCR patterns in animal feed or raw milk, as well as in UHT-treated or sterilized milk. However, the combined analysis of all typing data definitely showed that none of the farm isolates genetically resembled any of the UHT isolates. Since all farm isolates in this study originated from Belgium, they probably represent only a limited part of the natural genetic diversity of B. sporothermodurans. The data presented here confirm the hypothesis that the regular occurrence of contaminated UHT-treated and sterilized milk in some European dairy plants in the mid-1990s, as well as the present sporadic occurrence of contamination, can also be caused by circulation of the HRS clone within and between UHT production units. Occasionally, contamination of UHT-treated milk by a new genetic type (e.g., a type originating at the farm level) occurs, as exemplified here by the German UHT isolates. At present, the data obtained in this study do not favor or eliminate any of the potential contamination routes mentioned above.
In conclusion, this molecular typing study showed that a few clones of B. sporothermodurans, including the so-called HRS clone, have been and are still responsible for the contamination of UHT-treated and sterilized milk and milk products due to the production of highly heat-resistant spores. In particular, the HRS clone has spread over several European countries and even between continents. The strains isolated from UHT-treated and sterilized milk show a close genetic relationship, suggesting a common ancestry for the production of highly heat-resistant spores. An intriguing question which emerges, is whether the capacity to produce highly heat-resistant spores that allow survival after certain heat treatments is restricted to the subgroup of UHT isolates or whether it is a property more widespread in B sporothermodurnans, including the farm isolates. Although heat resistance is not an absolute spore property and is influenced by several factors, such as repeated laboratory cultivation (16), preliminary studies indicate that some farm isolates produce spores with remarkably high heat resistance. However, it remains to be determined whether these spores can also survive UHT treatment (O. Guillaume-Gentil, P. Keijzer, P. Scheldeman, and M. Heyndrickx, unpublished data). The molecular typing techniques used in this study demonstrated the great genetic heterogeneity of B. sporothermodurans isolates from dairy farms, even though they had been isolated in only one country. Because of this observed heterogeneity, the original taxonomic description of B. sporothermodurans, which was based on only a few genetically homogeneous UHT isolates (20), may no longer be adequate.
We thank E. Engels for performing REP-PCR, P. De Vos (University of Ghent, Ghent, Belgium) for the use of the BioNumerics software, and E. Bidlas for determination of ribopatterns.
|
|
|---|
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»