|
|
||||||||
Milk Quality Improvement Program, Department of Food Science, Cornell University, Ithaca, NY 14853
1 Corresponding author: kjb4{at}cornell.edu
| ABSTRACT |
|---|
|
|
|---|
Key Words: milk Bacillus Paenibacillus bacterial subtyping
| INTRODUCTION |
|---|
|
|
|---|
Bacillus and related genera represent a heterogeneous taxonomic group of gram-positive, rod-shaped, endospore-forming bacteria (Ash et al., 1991) characterized by a wide range of genetic and physiological traits, with variation in cell wall structure, DNA base composition (% mol G + C) and ability to metabolize carbohydrates (Sneath, 1986; Claus and Fritze, 1989). A number of bacteria initially classified as Bacillus spp. recently have been reclassified into new genera (Pettipher et al., 1997; Eiroa et al., 1999; Guinebretiere et al., 2001; Fromm and Boor, 2004; Huck et al., 2007), including a number that are recognized as contributing to food and beverage spoilage [e.g., Paenibacillus (Ash et al., 1993) and Alicyclobacillus (Wisotzkey et al., 1992)]. Although conventional phenotypic characterization methods (e.g., cellular and spore morphology, gram staining, and biochemical tests) are useful for initial identification and characterization of Bacillus and related genera (Logan and Berkeley, 1984), these methods may not provide the sensitive differentiation needed to facilitate tracking of contaminants through processing systems, particularly for the 16S ribosomal RNA (rRNA) group 3 bacilli (i.e., Paenibacillus spp.; Pires and Seldin, 1997; Odumeru et al., 1999; Guinebretiere et al., 2001; Teng et al., 2003; Durak et al., 2006).
Molecular subtyping methods, which allow the differentiation of subtypes within a bacterial species, can provide sensitive tools for tracking contamination sources (Kabuki et al., 2004), including endospore-forming bacteria (Durak et al., 2006; Huck et al., 2007). A variety of molecular subtyping methods [e.g., ribotyping, pulsed-field gel electrophoresis (PFGE), and different PCR-based methods] have been used for subtyping bacterial pathogens and spoilage microorganisms (Wiedmann, 2002). Deoxyribonucleic acid sequencing-based subtyping approaches are emerging as rapid, cost-effective, and reproducible methods (Aires-de-Sousa et al., 2006). Specifically, a newly developed subtyping method for Bacillus and related genera, which is based on DNA sequencing of the rpoB gene, provides a tool for reproducible and economic subtype characterization that can be used to identify the sources and transmission of these spoilage microorganisms from raw materials to finished products (Durak et al., 2006; Huck et al., 2007).
The objective of this study was to determine whether heat-resistant spoilage bacteria present in pasteurized fluid milk products were also present in individual tank truck loads of raw milk used to manufacture those products. To that end, we used the rpoB sequence-based subtyping method (Durak et al., 2006) to characterize the diversity of psychrotolerant endospore-forming bacteria in raw milk tank truck loads and in pasteurized products from 2 fluid milk plants. In addition, molecular subtyping data for isolates collected in this study were compared with subtyping data from isolates collected in a previous study, which traced psychrotolerant endospore-forming bacteria from commingled raw milk silos to pasteurized packaged products (Huck et al., 2007).
| MATERIALS AND METHODS |
|---|
|
|
|---|
300 mL) were collected by each plants raw milk receiving personnel from each compartment of each tank truck during 1 d (site A; Figure 1
300 mL) was collected by plant personnel from each of 2 commingled raw milk silos by using the silos sampling ports (site B; Figure 1
|
Bacterial Isolation
For each 200-mL milk sample, bacterial colonies present on SPC and PBC plated at d 14 or 16 postprocessing were visually examined, and colonies with distinct morphologies (ranging from 1 to 3 colonies per sample) were selected and streaked for purity on brain heart infusion agar. Purified isolates were characterized for gram reaction by using a 3-step Gram Stain Kit (Becton, Dickinson and Co., Sparks, MD) and subsequently frozen at –80°C in 15% glycerol.
rpoB Sequencing
A total of 101 isolates were further characterized by molecular subtyping (39 and 62 isolates from plants C and M, respectively). All isolates from plant C were obtained from samples plated at d 14. For plant M, isolates from heat-treated raw milk samples were obtained from samples plated for SPC at d 14 postheating, whereas isolates from pasteurized milk samples were obtained from samples plated on d 16 postprocessing.
Subtyping was performed by determining the DNA sequence for a 632-nucleotide (nt) fragment of the rpoB gene, which encodes the beta subunit of RNA polymerase, by using the rpoB PCR primers described by Drancourt et al. (2004) as detailed by Durak et al. (2006). This method was selected because it allows for phylogenetic characterization of isolates in addition to subtype identification, and because it is more economical than most banding pattern-based methods such as ribotyping or PFGE (Kabuki et al., 2004). Briefly, purified DNA for PCR was prepared for each isolate from overnight brain heart infusion cultures by using protocol D of the QIAmp DNA Mini Kit (Qiagen Inc., Valencia, CA). The rpoB PCR products were purified by using the QIAquick PCR Purification Kit (Qiagen Inc.) and quantified by using a NanoDrop ND-1000 spectrophotometer (Nano-drop Technologies, Wilmington, DE). Bidirectional sequencing with PCR primers and Big Dye Terminator chemistry was performed at Cornell Universitys BioResource Center (Ithaca, NY). Deoxyribonucleic acid sequences were assembled and proofread in SeqMan (Lasergene, DNAstar, Madison, WI) and only high-quality, double-stranded sequence data were used for further analyses. Sequences were aligned in MegAlign (Lasergene, DNAstar) and 632-nt rpoB fragments (corresponding to nt 2,455 to 3,086 of the 3,534-nt rpoB open reading frame of Bacillus cereus ATCC 10987; GenBank AE017194, locus tag BCE_0102) were used for subsequent analyses. The software DnaSP, version 4.0 (Rozas and Rozas, 1999) was used for rpoB allele assignment, and different allelic types (AT) were assigned to gene sequences that differed from each other by one or more nt. Isolates with 2 different AT were considered to represent 2 different molecular subtypes.
Cluster Analyses of rpoB Sequencing Data
Identification of related rpoB AT and confirmation of genus and species for the isolates subtyped in this study were accomplished by aligning the 31 rpoB sequences representing the 31 rpoB AT identified in this study along with rpoB sequences representing the 102 rpoB AT previously characterized by Durak et al. (2006) and Huck et al. (2007) (GenBank Accession Numbers EF156897 to EF156925 and EF156999 to EF15702). A phylogenetic tree was constructed by using the neighbor-joining method with 2,000 bootstrap replicates; Streptococcus pyogenes (NC_004070) was used as an outgroup. Phylogenetic analysis was performed by using PAUP, version 4.0 beta 10 (Sinauer Associates, Inc., Sutherland, MA). The 102 previously identified rpoB AT used to construct phylogenetic trees in this study had previously been used to assign isolates to genus (Bacillus and Paenibacillus) and species, where possible, by using both rpoB and 16S rDNA sequence data (Huck et al., 2007). The resulting phylogenetic trees support the genus and species identification of the isolates characterized here.
16S rDNA Sequencing
Although the majority of isolates collected in the current study could be characterized to genus and species through phylogenetic analyses with previous rpoB sequence data (Durak et al., 2006; Huck et al., 2007) as described above, 16S rDNA sequencing was used to confirm the genus and species identification of rpoB AT that had not been identified in any of our previous studies. Specifically, one isolate representing each newly identified Bacillus or Paenibacillus spp. rpoB AT was characterized by sequencing the 3' end of the 16S rDNA. In addition, isolates for which the rpoB primers (i.e., those described by Drancourt et al., 2004) did not successfully amplify the targeted rpoB fragments were similarly classified to genus or species with 16S rDNA sequencing. The 16S rDNA PCR was performed with primers PEU7 (Rothman et al., 2002) and DG74 (Griesen et al., 1994) as previously described (Fromm and Boor, 2004; Durak et al., 2006). Subsequent DNA sequencing of PCR products was performed with primers PEU7 and P3SH (Ralyea et al., 1998). Purification of genomic DNA, PCR product purification, quantification, sequencing, and sequence assembly were performed as described above for rpoB sequencing. To account for the presence of multiple copies of rRNA operons with different sequences in a given isolate (Klappenbach et al., 2001), 16S rDNA sequences were assembled by using nt ambiguity codes as described by the Nomenclature Committee of the International Union of Biochemistry and Molecular Biology. Final partial 3' 16S rDNA sequences were used for similarity searches against the National Center for Biotechnology Information nt sequence database with the Blast Local Alignment Search Tool (BLAST; McGinnis and Madden, 2004). Genus or species assignments, or both, for specific 16S rDNA sequences were based on the top matches returned by the BLAST search.
Data Curation
Isolate characteristics and sequence data, along with relevant isolation and sample information, can be accessed through the Pathogen Tracker database at www.pathogentracker.net (Fugett, 2006). GenBank accession numbers for 16S rDNA and rpoB sequences of isolates representing each rpoB AT identified in this study are listed in Table 1
.
|
| RESULTS |
|---|
|
|
|---|
Microbiological Analysis of Heat-Treated Raw and Commercially Pasteurized Milk Samples
Total bacterial counts, as measured by SPC and PBC, of heat-treated raw and commercially pasteurized milk samples that had been stored at 6°C were plotted to assess bacterial numbers over time in samples from plants C and M (Figures 2
and 3
). Within each plant, patterns of bacterial counts, as measured by SPC and PBC, were similar for heat-treated raw milk samples and commercially pasteurized milk samples, confirming previous observations (Huck et al., 2007) that the 80°C heat treatment used to treat all raw milk samples prior to refrigerated storage at 6°C leaves a residual bacterial population in the raw milk with growth characteristics similar to those in commercially pasteurized milk. Further, because the laboratory-heated milk was filled aseptically into sterile bottles, the fact that bacterial growth patterns in the paired commercial and laboratory-heated samples were essentially identical also suggests that minimal additional post-HTST pasteurization contamination with organisms that are able to grow at 6°C occurred in the commercially processed milk products.
|
|
Psychrotrophic bacterial counts measured over the shelf life of the heat-treated samples of raw milk collected from tank trucks differed significantly (P < 0.05; ANOVA) between plants (Figure 3
). For plant C, PBC for the raw milk samples that had been heat treated and then held at 6°C increased from 1.23 average log cfu/mL on d 1 (<1.00 to 1.60 log cfu/mL; 8 samples) to 4.20 average log cfu/mL on d 14 (2.51 to 5.89 log cfu/ mL; 8 samples; Figure 3
). For plant M, PBC for these heat-treated raw milk samples increased from <1.00 log cfu/mL on d 1 for all 13 samples to 5.76 average log cfu/mL on d 14 (3.96 to 6.73 log cfu/mL; 13 samples). Psychrotrophic bacterial counts for heat-treated samples of raw milk collected from 2 silos sampled at plant M increased from <1.00 log cfu/mL on d 1 (2 samples) to 6.51 average log cfu/mL on d 14 (6.43 and 6.58 log cfu/mL for samples 1 and 2).
Standard plate counts for milk from packaged fin-ished products at plant C increased from 2.17 average log cfu/mL on d 1 (2.08 to 2.30 log cfu/mL; 5 samples) to 4.60 average log cfu/mL on d 14 (3.43 to 5.82 log cfu/ mL; 5 samples). Standard plate counts for milk from packaged finished products at plant M increased from 2.72 average log cfu/mL on d 3 (2.54 to 2.89 log cfu/mL; 6 samples) to 5.90 average log cfu/mL on d 16 (4.79 to 6.68 log cfu/mL; 6 samples; Figure 2
).
Psychrotrophic bacterial counts for milk from packaged finished products at plant C increased from 1.45 average log cfu/mL on d 1 (1.00 to 2.18 log cfu/mL; 5 samples) to 4.64 average log cfu/mL on d 14 (3.40 to 6.78 log cfu/mL; 5 samples). Psychrotrophic bacterial counts for milk from packaged finished products at plant M increased from <1.00 log cfu/mL on d 3 for all 6 samples to 6.42 average log cfu/mL on d 16 (5.68 to 6.78 log cfu/mL; 6 samples; Figure 3
).
rpoB Sequencing
A total of 101 isolates were selected from d 14 and 16 SPC plates for rpoB sequence subtyping (representing 39 and 62 isolates from plants C and M, respectively). Among these 101 isolates, 88 yielded rpoB products that could be successfully sequenced. For 3 isolates, sufficient growth could not be obtained from cultures that had been frozen; therefore, subtyping was not performed. The remaining 10 isolates, for which the rpoB primers (i.e., those described by Drancourt et al., 2004) did not successfully amplify the targeted rpoB fragments, were subsequently characterized by sequencing the 3' end of the 16S rRNA gene.
Overall, 31 rpoB AT (i.e., subtypes based on different rpoB sequences) were differentiated among the 88 isolates (32 and 56 isolates from plants C and M, respectively). Eight AT were isolated from both plants, whereas 10 and 14 AT were isolated either at plant C or at plant M, respectively. A total of 17 rpoB AT represented only a single isolate, including 7 from plant C and 10 from plant M. Ten rpoB AT (representing 29.5% of all isolates subtyped) represented
2 and
5 isolates. Four rpoB AT (AT 2, 15, 27, and 32) represented more than 5 isolates; these 4 AT represented 51.1% of all isolates identified.
Cluster Analysis of rpoB Allelic Types
The rpoB sequences were used in phylogenetic analyses to identify groups of related subtypes (i.e., clusters) and to confirm genus and species identification of isolates. Specifically, 31 rpoB sequences, representing the 31 rpoB AT identified in this study, were aligned with rpoB sequences for the 102 rpoB AT previously characterized by our group (Durak et al., 2006; Huck et al., 2007); this alignment of 133 sequences was then used to construct a phylogenetic tree (Figure 4
). All rpoB AT identified in this study grouped into 2 distinct and well-supported clusters representing the genera Bacillus and Paenibacillus, consistent with our previous findings (Huck et al., 2007). Ten of the rpoB AT identified among the 88 isolates characterized here had not been identified in our previous studies (Durak et al., 2006; Huck et al., 2007). These 10 AT were classified as Bacillus (3 AT) and Paenibacillus spp. (7 AT) AT; all of these 10 "new" rpoB AT grouped closely with Bacillus and Paenibacillus spp. clusters previously identified by Huck et al. (2007; Table 1
) and could thus be assigned to genus, species, and phylogenetic cluster.
|
All of the 10 isolates, for which the rpoB primers (i.e., those described by Drancourt et al., 2004) did not successfully amplify the targeted rpoB fragments, were successfully amplified by 16S rDNA PCR. Six isolates shown by staining to be gram-positive rods were identified as 2 Bacillus spp. and 4 Paenibacillus spp. The remaining 4 isolates were not identified as gram-positive rods and were classified by 16S rDNA sequencing as 3 Pseudomonas spp. and 1 Enterococcus sp.
Raw Milk Contamination Patterns
A total of 23 different rpoB AT were identified among the 63 isolates that were obtained from 23 heat-treated raw milk samples; raw milk samples represented 8 samples collected from tanker trucks at plant C and 13 samples collected from tanker trucks at plant M as well as 2 samples from raw milk silos at plant M. For plant C, 12 different rpoB AT were identified among the 22 isolates subtyped (Table 2
); only 2 rpoB AT were found in more than one truck load at this plant, including AT 15 (2 loads) and AT 32 (2 loads; Table 2
). For plant M, 15 different rpoB AT were identified among the 35 subtyped isolates obtained from the raw milk samples collected from tank truck loads. Four rpoB AT were identified in more than one truck load, including AT 2 (3 loads), AT 15 (8 loads), AT 27 (2 loads), and AT 32 (2 loads; Table 3
). Four rpoB AT were identified among the 6 isolates collected from heat-treated samples of raw milk obtained from the raw milk silos at plant M. Although one of these rpoB AT (AT 40) was isolated only from the raw milk silo samples at this plant, 3 AT (AT 15, 21, and 27) had also been isolated from the tank truck samples at plant M (AT 15 was identified in both silos, whereas AT 21 and 27 were identified only in silos 2 and 3, respectively).
|
|
Contamination Patterns from Raw Milk to Pasteurized Packaged Products
A total of 14 rpoB AT were identified among the 25 isolates obtained from 11 packaged product samples (representing 5 and 6 samples from plants C and M, respectively). Nine and 10 rpoB AT were identified among isolates from plant C (10 isolates) and plant M (15 isolates), respectively. Five of these rpoB AT (AT 2, 15, 19, 27, and 35), representing 64% of the packaged product isolates, were identified at both plants.
At plant M, 5 rpoB AT (AT 2, 15, 21, 27, and 35) were found in both heat-treated raw milk and finished product samples. Similarly, at plant C, 3 rpoB AT (AT 15, 27, and 46) were found in both heat-treated raw milk and finished product samples. Interestingly, at plant M, isolates with AT 15 and 27 were also found in raw milk samples from tank trucks, in raw milk samples from silos, and in pasteurized packaged products, suggesting transmission of these AT from incoming tank truck raw milk loads to finished products in both plants. Allelic type 15 was also found in raw and pasteurized milk samples in plant M in a previous study that included 3 samplings in 3 subsequent weeks (Huck et al., 2007). Furthermore, AT 2, which was found in plant M only in pasteurized milk samples in each of 3 wk in a previous study (Huck et al., 2007), was found in 4 of the 6 finished product samples and 3 raw milk samples (Table 3
). Importantly, all 3 of these AT (AT 2, 15, and 27) represent a specific Paenibacillus cluster (cluster XXII; Table 1
). Although this cluster has not yet been classified to the species level, our data indicate that this Paenibacillus spp. is a major cause of fluid milk spoilage, particularly because cluster XXII represents 77.3% of the 88 isolates characterized here as well as 44.7% of 385 previously characterized isolates from heat-treated raw and pasteurized fluid milk samples (Huck et al., 2007).
| DISCUSSION |
|---|
|
|
|---|
rpoB Sequencing Allows for Rapid and Economical Characterization and Subtyping of Endospore-Forming Dairy Spoilage Organisms
Application of molecular subtyping methods, including PFGE, ribotyping, and other banding pattern-based methods, has been critical for gaining a better understanding of the transmission of foodborne pathogens and for reducing transmission of these pathogens (Swaminathan et al., 2001). Some endospore-forming food spoilage microorganisms have been characterized by banding pattern-based subtyping methods (Schraft et al., 1996; Svensson et al., 1999), but these subtyping methods have been applied less frequently to the exploration of spoilage microorganism transmission, at least partly because of the cost of many modern molecular subtyping methods (Olive and Bean, 1999). Recent advances in the development of DNA sequencing-based subtyping methods, including the development of an rpoB sequencing-based subtyping method for Bacillus and Paenibacillus spp., now offer improved opportunities to probe the transmission and ecology of these organisms with an economical subtyping method. In addition, the development and availability of an rpoB sequence database that includes a diverse collection of >500 Bacillus and Paenibacillus spp. isolates (this study; Durak et al., 2006; Huck et al., 2007) comprising >110 rpoB AT that have also been characterized by 16S rDNA sequencing allows for characterization of many Bacillus and Paenibacillus isolates to the genus and species level. Even if species-level identification is not possible, which is currently typical for many Paenibacillus isolates, rpoB sequence data still allow for identification of distinct Paenibacillus clusters, which are likely to represent different species. Furthermore, these sequencing-based subtyping strategies allow simple Web-based data curation as well as reliable and robust comparison of subtype data generated by different research groups (Aanensen and Spratt, 2005).
Paenibacillus spp. Are Major Contributors to Fluid Milk Spoilage
Overall, our data indicate that a considerable subtype diversity of endospore-forming bacteria, particularly those representing the genera Bacillus and Paenibacillus, are present in raw and processed milk. In conjunction with previous studies (Huck et al., 2007), our data specifically indicate that a large proportion of endospore-forming bacteria isolated from raw and pasteurized fluid milk represent Paenibacillus spp. In particular, most of the Paenibacillus isolates from milk characterized here and previously (Huck et al., 2007) represent a specific Paenibacillus subgroup (cluster XXII). Bacillus spp., as well as Paenibacillus spp., are known for their ability to form heat-resistant endospores that can survive heat-treatment processing conditions (Collins, 1981; Huck et al., 2007). Although a number of studies have identified Bacillus spp. in milk from the farm (Crielly et al., 1994; Sutherland and Murdoch, 1994; Lukasova et al., 2001) and the processing plant, including in processing plant raw milk supplies (Lin et al., 1998; Huck et al., 2007) and in pasteurized packaged products (Griffiths and Phillips, 1990; Lin et al., 1998; Douglas et al., 2000; Fromm and Boor, 2004; Durak et al., 2006; Huck et al., 2007), Paenibacillus has only recently been identified as an important contributor to fluid milk spoilage that can be found in both raw and pasteurized fluid milk (Huck et al., 2007). The present study provides direct subtype-based evidence that Paenibacillus spp. present in raw milk (particularly those representing cluster XXII) can be transmitted to finished pasteurized products. Similar to Bacillus spp., Paenibacillus spp. that are able to grow under refrigeration temperatures (typically <6°C) are of particular concern as spoilage microorganisms in fluid milk products. Bacterial growth patterns in the heat-treated raw milk samples reported here and previously (Huck et al., 2007) indicate that Paenibacillus spp. can grow in milk stored at 6°C. This growth behavior is consistent with reports that a number of Bacillus spp., including B. cereus (Dufrenne et al., 1995; te Giffel et al., 1995), Bacillus weihenstephanensis (Pacova et al., 2003), and Bacillus circulans (Shehata and Collins, 1971), as well as Paenibacillus spp. (Huck et al., 2007) are psychro-tolerant.
The data reported here indicate that Paenibacillus and Bacillus spp. are commonly present in the raw milk entering processing plants. These microorganisms may contaminate raw milk at various points during production (e.g., milking, raw milk storage and handling on the farm) as well as possibly during loading of tank trucks and transport. Specifically, aerobic spore-formers of the genus Bacillus are ubiquitous and have been isolated from dairy farm environments, including silage (te Giffel et al., 2002), pasture (Slaghuis et al., 1997), soil (Christiansson et al., 1999), and fecal material (Labots et al., 1965). Environmental contamination of udders (Lukasova et al., 2001) and milking equipment (McKinnon and Pettipher, 1983) may lead to the presence of Bacillus spp., including B. licheniformis and B. cereus, in bulk tanks on dairy farms (Griffiths and Phillips, 1990; Crielly et al., 1994). In addition, a previous study by our group (Huck et al., 2007) has found preliminary evidence for the persistence of selected Paenibacillus and Bacillus subtypes in processing plants. Thus, contamination of fluid milk with endospore-forming spoilage microorganisms may also occur at the processing plant, either pre- or postpasteurization.
| CONCLUSIONS |
|---|
|
|
|---|
| ACKNOWLEDGEMENTS |
|---|
|
|
|---|
Received for publication March 13, 2007. Accepted for publication June 13, 2007.
| REFERENCES |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
M. L. Ranieri, J. R. Huck, M. Sonnen, D. M. Barbano, and K. J. Boor High temperature, short time pasteurization temperatures inversely affect bacterial numbers during refrigerated storage of pasteurized fluid milk J Dairy Sci, October 1, 2009; 92(10): 4823 - 4832. [Abstract] [Full Text] [PDF] |
||||
![]() |
M. L. Ranieri and K. J. Boor Short communication: Bacterial ecology of high-temperature, short-time pasteurized milk processed in the United States J Dairy Sci, October 1, 2009; 92(10): 4833 - 4840. [Abstract] [Full Text] [PDF] |
||||
![]() |
J. R. Huck, M. Sonnen, and K. J. Boor Tracking Heat-Resistant, Cold-Thriving Fluid Milk Spoilage Bacteria from Farm to Packaged Product J Dairy Sci, March 1, 2008; 91(3): 1218 - 1228. [Abstract] [Full Text] [PDF] |
||||
| ||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
| HOME | HELP | FEEDBACK | SUBSCRIPTIONS | ARCHIVE | SEARCH | TABLE OF CONTENTS |