Next Article in Journal
Consumption of Plant-Derived Phenolic Acids Modulates Hunger and Satiety Responses Due to Chemical Interactions with Enteroendocrine Mediators
Previous Article in Journal
Enhancement of Nutrient, Trace Element, and Organic Selenium Contents of Ratooning Rice Grains and Straw Through Foliar Application of Selenite
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Spore-Forming Clostridia in Raw Cow Milk from Northern Italy: A Trend Analysis over the Past 20 Years

1
National Reference Center for Bovine Milk Quality, Via A. Bianchi 9, 25124 Brescia, Italy
2
“Produzione Primaria” Department, Istituto Zooprofilattico Sperimentale della Lombardia e dell’Emilia Romagna “Bruno Ubertini” (IZSLER), Via A. Bianchi 9, 25124 Brescia, Italy
*
Author to whom correspondence should be addressed.
Foods 2024, 13(22), 3638; https://doi.org/10.3390/foods13223638
Submission received: 17 October 2024 / Revised: 7 November 2024 / Accepted: 12 November 2024 / Published: 14 November 2024
(This article belongs to the Section Dairy)

Abstract

:
Clostridium species are known for their impact on animal and human health, but also for the spoilage of foodstuffs. Their spores contaminate milk and result in germination and gas production, the latter being particularly evident in the cheeses that suffer severe depreciation. To address this issue, the Primary Production Department of the IZSLER institute in Brescia, Italy conducts the Most Probable Number (MPN) method on bovine milk samples collected from Northern Italian dairies between 2004 and 2023. This approach leverages two semi-quantitative protocols, S2 and S3, to detect Clostridium species spore forms upon customer request. Here, we would like to present an a-posteriori analysis on the results of the S2 and S3 protocols. The goal of this study is to highlight the differences between these two methods and provide evidence of the actual decrease in Clostridium species in raw cow milk over a 20-year period. Our analysis shows that client demand for S2 has progressively decreased, while S3’s has remained constant, and both protocols reveal a significant reduction in positives; furthermore, S3’s greater sensitivity made it more responsive to environmental changes. This highlights the necessity of choosing the appropriate testing protocol that accounts for both regulatory standards and environmental factors. Overall, our findings underscore the importance of continued monitoring to manage Clostridium species contamination and ensure milk quality.

1. Introduction

The genus of Clostridium consists of gram-positive spore-forming anaerobic bacteria responsible for quality defects in food and feed, especially in hard and semi-hard cheese such as Grana Padano PDO (Protected Designation of Origin), Gouda, Emmental, and Gruyere [1]. The main species involved in cheese defects are Clostridium tyrobutyricum, Clostridium butyricum, Clostridium sporogenes, and Clostridium beijerinckii, all of which are included in the Butyric Acid Bacteria group (BABs), as they produce butyric acid. In particular, C. tyrobutyricum and C. beijerinckii are the main gas producers; hence, they are primarily responsible for late defects [2].
BABs convert lactate into gases (H2 and CO2), acetic acid, and butyric acid causing Late Blowing Defects (LBDs), cracks, cavities called “eyes”, shredding, and openings in the central part of the form, changes in flavor, and a spongy consistency that decrease the commercial value of cheese [3,4]. All these defects can occur at a very low concentration of spores, approximately 100 to 1000 clostridial spores per liter of milk [5,6]. This is why it is critical to pay close attention to clostridial milk contamination from a variety of sources, including crops, bedding, bulk tanks, feces, poor handling, contaminated equipment, and soils [7,8].
Clostridium spores are a major problem, as they are resistant to various cheese-manufacturing procedures. For example, pasteurization kills most vegetative bacterial cells, but not clostridial spores, which remain quiescent until germination [9,10]. A few weeks into ripening, the physical and chemical conditions of cheese become optimal for reactivation of the clostridial spores, leading to LBDs. Reducing the spore concentration in milk is possible through outcropping—a natural stage in some cheese production—or other methods such as bactofugation, microfiltration, or the addition of bacterial strains that compete with Clostridium for lactic acid metabolism [11]. The use of preservatives or adjuvant technology like nitrates and lysozyme (500 U/mL cheese milk) is also effective, especially if the clostridial spore concentration is less than 300 per liter [12,13]. However, these substances are permitted only under specific productions (e.g., up to 2.5 g of lysozyme per 100 kg of milk in Grana Padano PDO production [14]) or are completely forbidden (e.g., in Parmigiano Reggiano PDO production [15]).
In 2015, 15–35% of Grana Padano PDO production displayed non-compliant forms that presented defects partially caused by Clostridium [16]. In 2022, these defects decreased to 2% [17,18,19], but still led to an economic loss of EUR ~36 million [20]. This highlights the importance of detecting clostridial spores in milk samples for a prediction of the ripening outcomes. In 2016, Brändle at al. proposed a classification of the available laboratory methods to quantify the spore levels in milk [21]. Most of these approaches are culture-dependent and rely on plate counting or a Multiplex-PCR to distinguish the different Clostridium species present in the samples and ensure a precise quantification [4,11]. However, the high costs and small number of samples that can be analyzed in a day pose significant challenges [16,22]. The rapid, automated methods for spore detection that are commercially available (e.g., AMP6000) have a higher selectivity medium, but often require expensive equipment that most laboratories cannot afford [23,24].
The Most Probable Number (MPN) method on a liquid medium is currently the most widely used method for assessing the contamination of cow’s milk by sporulating forms of Clostridium spp. The process involves seeding the milk sample into tubes, pasteurizing either before or after the seeding, and adding a paraffin plug. If the spores in the sample germinate and produce CO2, the plug will detach from the surface, indicating contamination [21]. Notably, the MPN method cannot differentiate between germinated Clostridium species or determine their relevance to food spoilage.
In the Primary Production Department of the IZSLER institute in Brescia, Italy, the MPN method, as reported in the International Organization for Standardization (ISO) 7218:2013 [25], is employed to semi-quantitatively assess clostridial contamination in milk. This method is a commercial parameter within the Milk Quality Payment framework [26,27], determining milk premiums or penalties based on quality. It also adheres to Regulation (EC) No 853/2004 of the European Parliament and Council, which sets specific hygiene rules for food of animal origin, including dairy products [28].
The purpose of this study is to analyze the Milk Quality Payment data collected by IZSLER’s Primary Production Department and the National Reference Centre for Cow Milk Quality for the years 2004–2023 applying two MPN protocols, S2 and S3. Although conceptually similar, these two protocols differ in sensitivity according to the number of dilutions performed on the bovine samples: two for S2 and three for S3. In addition, the S2 protocol is used as a general screening test, whereas the S3 protocol is mainly used on milk intended for cheesemaking, which therefore needs to follow more accurate controls regarding sporal contamination. Here, we provide an overview of the trend of Clostridium species in the field samples submitted to the institute by breeders and dairies in Northern Italy. Since more than 40% of lactating cows with dirty udders can increase the average contamination of the spores in milk by 15% [29,30,31], evaluating these trends helps with understanding how the evolution of farm practices impacts the presence of sporigenic bacteria in milk. Even though BABs bacteria causing late defects in dairy products do not pose a risk to consumers [21,32], the Clostridium genus remains potentially toxic, neuro-toxigenic in addition to food spoiling [33]. Our 20-year review of the S2 and S3 analyses for the detection of Clostridium species shows that the constant application of the MPN method has helped control the spread of this bacterial genus [34].

2. Material and Methods

2.1. Sample Collection

As the alterations induced by LBDs affect the quality of long-ripened cheeses made from raw cow milk, we investigated the hygienic and commercial parameters of the raw milk samples destined for cheese-making. All the raw milk samples analyzed and reported in this study were taken by samplers, delivered to IZSLER by private companies in a refrigerated state, processed on the day of arrival, and then discarded, as their shelf-life is 3 days.
The samples, or matrices, included in the analysis were cow milk, boiler milk, raw mass milk, and unpasteurized bovine whole milk. In contrast, goat milk, goat mass milk, sheep milk, buffalo milk, mare milk, bovine whey, pasteurized bovine milk, Ultra-High Temperature (UHT) bovine milk, and cream were excluded (Table 1).
Over 95% of the raw milk samples originated from Northern Italian provinces, specifically Lombardy and Emilia Romagna, particularly in the Grana Padano PDO area (Table 2). Consequently, our analysis provides insights on the contamination levels and regional prevalence of Clostridium species within this specific geographical region.

2.2. MPN Protocols

The MPN method presented here was performed according to two analogous protocols, S2 and S3, based on customer request. Both protocols analyze raw bovine milk samples diluted at 1:10 and 1:100 ratios; to these two dilutions in the S3 protocol a third analytical step is added where undiluted samples are analyzed.
The raw milk is diluted with a BUtyricum and TIrobutyricum (BUTI) medium, which consists of 0.1% casein tryptic digest, 0.1% meat extract, 0.03% yeast extract, 0.08% sodium acetate, 0.05% sodium chloride, 0.01% soluble starch, 25 mL of 50% sodium lactate, and 975 mL of demineralized water. In addition, in the S3 application of the MPN method, 0.3 mL of 50% sodium lactate, which promotes the germination of Clostridium spp. spores, must be added. Completely dissolve the reagents in demineralized water preheated to 50 °C; adjust the pH to 6.1 (37% hydrochloric acid is often needed) and sterilize the solution by autoclaving at 121 °C for 15 min. Store it at 4 °C for up to 30 days, away from direct light to prevent photo-oxidation.
A total of 10 mL of raw milk (undiluted, 1:10, 1:100) are seeded in sterile 16 mL × 160 mm glass tubes containing 2 mL of a paraffin and Vaseline oil mixture (3:1 ratio). For the S3 undiluted samples, 0.3 mL in 50% solution sodium lactate (C3H5NaO3) is directly added to the glass tube.
After seeding the medium, the tubes are pasteurized at 85 °C for 20 min and cooled at 15 °C for 10 min. This process causes the paraffin and Vaseline oil mixture to melt, then solidify, forming a non-permeable gas plug that seals the sample in a tight anaerobic environment. Following incubation at 37 °C for 7 days, the tubes are checked individually to assess whether the oil plug had detached from the underlying medium due to CO2 and H2 production, indicating the presence of germinated spores in the sample.
The MPN method explained here [25,27,35,36], can be performed manually or with automatic sample preparers such as the ones employed by the company Skalar; in our case, the analyses were carried out using automatic procedures [37]. Each protocol is performed in triplicate. For every 60 samples analyzed, a control is analyzed to verify that the machine is not affected by entrainment and that the medium in use is not contaminated. The standard consists of sterile distilled water, added in the same dilution as the analyzed milk sample.

2.3. MPN Validation

The MPN values were calculated with the Thomas formula [38]:
M P N = P n × T
where P is the number of positives in the set of seeded tubes, n is the volume of the inoculated sample (L), and T is the total volume of the reaction tube (L). The MPN method presented here was validated for accuracy, precision, specificity, and the Limit of Detection (LOD) as required by ISO 9001:2000, in force at the time of validation and recently updated to 17025:2018 [39]. The uncertainty of the estimated MPN value measured was expressed through the fiducial limits with a 95% confidence level (P = 95%, K ≈ 2) in the S2 protocol, as specified in ISO 7218:2007 Amd 1:2013 [25]. In contrast, for S3, the standards do not specify explicit calculation methods. Therefore, we determined the uncertainty as the sum of the relative uncertainties associated with dilution and the MPN estimation, as described in the bibliography [38]. To assess the robustness of the method, several factors were taken into account for their influence on the analytical result such as the incubation time.

2.4. Data Analysis

The sensitivity of the MPN method was set at 308 spores/L for the S2 protocol and 31 spores/L for the S3 protocol. Both methods have an upper detection limit of 11,000 spores/L, as required by ISO 7218:2007 Amd 1:2013 [25]. S3 has higher sensitivity due to the inclusion of an undiluted sample, enabling the detection of lower spore concentrations.
To make the classification of the positive samples more interpretable, the results were divided into weak (Pos) and strong (Pos+) positives based on the spore count ranges. For the S2 protocol, the samples within the 308–966 spores/L range were classified as Pos, while those with more than 967 spores/L were classified as Pos+. For the S3 protocol, the samples with 31–200 spores/L were considered Pos, and those with 201 spores/L or more were Pos+. This approach allowed us to capture meaningful variations without introducing unnecessary complexity.
The data were analyzed using one-way ANOVA (α = 0.01) on monthly frequencies, with post-hoc comparisons made using the Tukey–Kramer test (α = 0.01) to compare each year’s relative frequency with the 2004 average. The Shapiro–Wilk test (α = 0.01) was applied to check for normality within each year and category. The outliers were excluded using the Z statistic. All the analyses were performed using the Microsoft Excel statistical package.

3. Results

3.1. The Demand for S2 Analysis Declines as S3 Requests Increase

We applied the MPN method, specifically the S2 and S3 protocols, to assess spore presence and germination in the raw milk samples from the Northern Italian dairies delivered to IZSLER. Between 2004 and 2023, a total of 1,017,606 raw milk samples were analyzed with the S2 protocol, and 57,575 with the S3 protocol, averaging approximately 50,000 and 3000 samples per year, respectively (Figure 1A,B). The customer demand for an S2 analysis began declining in 2007 and had decreased by 31.4% by 2023. Meanwhile, the S3 requests remained steady, except for a notable increase between 2017 and 2019.
Grouping the data by season, no significant variations were observed in the samples delivered for an S2 analysis (Figure 1C). However, the demand for the S3 protocol showed wider fluctuations across the different seasons (Figure 1D).

3.2. Samples Positive for Clostridium Contamination Decreased in the Past 20 Years

By monitoring gas production and paraffin–Vaseline cap detachment in the tubes, we recorded the number of positive results at each dilution, along with the corresponding MPN values and the associated 95% confidence limits (Table 3).
Following this, we compared the average percentages of the total positive samples (Pos tot) each year with those from 2004 to identify any significant differences (Figure 2A). This revealed that when using the S2 protocol, the Pos tot decreased by ~16% starting in 2017, while with the S3 protocol the reduction was ~17% starting in 2015. Then, we separated the weak positive samples (Pos) from the strong positives (Pos+) for both S2 and S3, and carried out the same analysis (Figure 2B,C).
Since 2004, the S2 protocol has shown a decrease of ~8% in the number of Pos and ~4% in the number of Pos+, with the first significant peaks observed in 2017 and 2011, respectively. In contrast, the S3 protocol revealed a decrease of about 21% in the number of Pos+, with an increase in the number of Pos in 2008 and 2011. Further information can be found in the Supplementary Materials.
Finally, we examined the seasonal variations in the Pos tot, Pos, and Pos+ results. The S2 protocol showed no significant differences across the seasons (Figure 3A). However, the more sensitive S3 protocol revealed a significant increase in the Pos tot from July to September, except for in 2021, where the highest values occurred in October–December (Figure 3B). When breaking down the Pos tot into Pos and Pos+, similar patterns emerged. S2 detected no significant seasonal variation (Figure 3C,D). In contrast, S3 showed an increase in Pos from July to September, though this trend was not sustained by the Pos+, with the exception of October–December of 2019 (Figure 3E,F).

4. Discussion

The Clostridium genus consists of gram-positive spore-forming anaerobic bacteria that are responsible for quality defects in hard and semi-hard cheese [1,3,4,7]. Thus, the detection of the presence and concentration of Clostridium contamination in milk is crucial to prevent economic losses and safeguard public health [33]. In the Quality Milk Payment framework [26] in Italy, the commercial milk limit for clostridial spores is set at 100 spores/L to ensure product quality and safety [5,40]. Since a spore count cannot be measured directly, the institute identifies butyric Clostridium in milk by applying either the S2 or S3 protocols of the Most Probable Number (MPN) method [4,16,41] to quantify the spores per liter of milk samples.

4.1. S2 and S3 Demand Is Tied to Regulatory Requirements, Dairy Practices, and External Events

In this study, we reviewed 20 years of data on raw milk samples delivered to IZSLER and analyzed them using the S2 and S3 protocols. We examined the sample sizes over the years as well as the customer preferences for these protocols, and we identified the trends in Clostridium contamination. This analysis provides insight into how clostridial spore levels have evolved over time as well as how customer demand for testing protocols has shifted, and highlights the patterns related to annual and seasonal variations.
The first trend identified over the past 20 years lies in the customer demand for the S2 and S3 protocols. Our data show a significantly higher number of requests for the S2 protocol compared to S3 (~50,000 vs. 3000 samples per year, Figure 1A,B). This is largely explained by the fact that S2 is a mandatory screening parameter of the Quality Milk Payment Regulation [28], with all milk producers being required to take and analyze at least two samples per month. In contrast, S3 is not compulsory, but it is more sensitive and is thus preferred by the dairies producing medium-long ripening cheese and those using silage-fed cows.
Secondly, we observed a steady decline in S2 demand, while the requests for the S3 protocol remained constant, except for a brief dip between 2017 and 2019 (Figure 1A,B). These trends align with the broader shifts in dairy farming practices. Since 2010, there has been a gradual reduction in dairy cow herds, with smaller farms being replaced by larger ones that produce more milk, often destined for cheese-making [42,43]. This change in farming practices has influenced both the number of milk samples delivered for analysis and the type of tests requested, leading to a decrease in S2 demand and an increase in the demand for the more sensitive S3 protocol.
Comparing the seasonal demands for the protocols, the S2 requests remained constant across the seasons, while the S3 requests exhibited much greater variability (Figure 1C,D). This difference is mainly due to the distinct nature of the two methods. The S2 protocol, being mandatory under the Quality Milk Payment Regulation, is consistently requested and routinely performed. In contrast, S3 is an optional protocol used by diaries at their discretion, often depending on the specific needs related to cheese-making or other factors.
Interestingly, our analysis revealed an anomaly in the number of milk samples analyzed with the S2 protocol during 2020, which can be attributed to the disruptions caused by the COVID-19 pandemic. Specifically, in March and April of 2020, the regions of Lombardy and Emilia Romagna granted temporary exemptions for the analysis of milk samples as part of the Quality Milk Payment regulations. This exemption led to a sharp decrease in the number of S2 protocol samples analyzed during this period (Figure 1A,C). In contrast, the S3 protocol, which was not subject to regulatory requirements, did not show any fluctuations in sample numbers during this time (Figure 1B,D). This indicates that while the pandemic and subsequent regulatory changes impacted the S2 protocol, the demand for the S3 protocol remained constant.

4.2. Regulatory Changes and Seasonal Variation Impact on Clostridium-Positive Milk Samples

Our data show a reduction in the total positive samples (Pos tot) observed in both the S2 and S3 protocols (Figure 2A), likely reflecting the impact of the stricter animal welfare and hygiene regulations introduced from 2017 to 2018 [44,45,46]. These measures, such as mandatory udder cleaning, have been proven to reduce clostridia spore transfer in raw milk [47,48], despite the fact that many factors may affect together or independently the ripening process of hard and semi-hard cheese, such as dairy technology, autochthonous starters, and additives. Despite the similar behavior of the two protocols, the results diverged when comparing the total (Pos tot), weak (Pos), and strong (Pos+) positives. S2 detected fewer contaminated samples for both the Pos and Pos+ (~8% and ~4%, Figure 2B), while S3 recorded a more pronounced drop in total positives, mainly due to the reduction in the number of Pos+ (21%, Figure 2C), and an almost unchanged trend in low positives. This highlights the greater sensitivity of the S3 protocol in detecting lower spore concentrations.
This conclusion is further emphasized when comparing the positivity differentials of each year to 2004 (Figure 4). All the positive samples identified through S2 follow relatively consistent trends (Figure 4A). In contrast, S3 reveals more pronounced differences across the Pos tot, Pos, and Pos+, with the Pos+ having a dominant influence on the Pos tot (Figure 4B). Notably, between 2011 and 2015, an increase in the Pos seemed to offset the drop in Pos+. It is likely that, during this period, dairy companies began implementing strategies to combat Clostridium contamination. While these actions appeared to target the strong positives (Pos+), the rise in weaker positives (Pos) suggests that these initial countermeasures were not fully effective in eradicating contamination but may have reduced its severity.
The trend towards a significant decrease in spore levels could provide valuable information for dairies, especially those increasingly inclined to replace silage with more environmentally friendly feeds, such as hay from stable meadows or other types of fodder. This shift could eventually reduce the need for lysozyme and promote the introduction of lactic flora, which can naturally inhibit the germination of the few remaining spores [49].
Seasonal patterns further highlight the differences between the two protocols.
While S2 showed no significant seasonal changes, S3 revealed a noticeable rise in the Pos tot, Pos, and Pos+ during the July–September period, with occasional annual peaks in October–December (Figure 3B,D,F). This suggests that S3 is more responsive to the environmental changes that affect spore contamination.
For the milk samples processed through the S2 protocol, there was no consistent seasonal trend (Figure 3A,C,D), but deviations from the mean occurred in specific years, likely due to climatic factors such as humidity and elevated temperatures. Interestingly, the raw cow milk tested with the S3 protocol during the summer showed a significant increase in Clostridium spores, probably due to the elevated temperatures and the cows’ physiological changes. This is supported by the studies demonstrating that environmental factors can significantly influence microbial levels [16]. Research showing a correlation between increased Clostridium spores and warmer climatic conditions supports the idea that temperature, rather than fixed seasonal variability, may play a major role in fluctuating contamination levels [50].
Furthermore, the changes observed in positivity rates during the warmer seasons may be attributed not only to the improved survival and adaptation of Clostridium spores in higher temperatures but also to the metabolic and physical changes in cows during the summer [51]. Additionally, it is essential to consider how evolving climatic conditions—both on farms and in dairies—can impact contamination levels through factors such as dust, humidity, and the type of bedding used.

5. Conclusions

Butyric acid bacteria are the main culprits of various types of defects in cheeses. Our a-posteriori analysis aims to describe how their prevalence in raw cow milk has changed over a 20-year period. The S2 and S3 protocols of the MPN method, applied to the bovine samples delivered to IZSLER from its area of competence, have shown reductions in Clostridium contamination in raw milk from 2004 to 2023. This reflects the positive impact of hygiene and animal welfare practices. While both methods demonstrate overall improvements, the S3 protocol, with its higher sensitivity, is more effective at detecting low spore concentrations, making it especially suitable for the targeted monitoring of milk destined for cheese-making. Finally, climatic factors, such as humidity and environmental dust, significantly influence contamination, underscoring the importance of regular monitoring.
In conclusion, our results highlight the need to choose appropriate testing protocols that consider both regulatory standards and environmental factors. This ensures the effective monitoring of barn hygiene conditions and quality control, providing clues to ensuring the welfare of the sampled animals and the quality of the resulting dairy products.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods13223638/s1, Table S1: Positivity differential between 2004 and the comparison years for total positives (Pos tot), weak positives (Pos) and strong positives (Pos+). Asterisks (*) in black and red indicate a decrease and increase in positives, respectively, Figure S1: 20-year annual relative frequency of positive samples of the S2 (A) and S3 (B) analyses. Negative, weak positive (Pos) and strong positive (Pos+) samples are shown as green line, yellow bars and red bars respectively.

Author Contributions

Conceptualization and project administration, P.D.; methodology, formal analysis, and writing—original draft preparation, review, and editing, A.G.; data curation, L.G.; investigation and resources, P.B.; validation, G.B.; visualization, G.Z.; supervision, F.P.; writing—review and editing, P.D., G.B. and G.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The original contributions presented in the study are included in the article and Supplementary Material, further inquiries can be directed to the corresponding author.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Bassi, D.; Fontana, C.; Gazzola, S.; Pietta, E.; Puglisi, E.; Cappa, F.; Cocconcelli, P.S. Draft Genome Sequence of Clostridium Tyrobutyricum Strain UC7086, Isolated from Grana Padano Cheese with Late-Blowing Defect. Genome Announc. 2013, 1, e00614-13. [Google Scholar] [CrossRef] [PubMed]
  2. Arias, C.; Oliete, B.; Seseña, S.; Jimenez, L.; Pérez-Guzmán, M.D.; Arias, R. Importance of On-Farm Management Practices on Lactate-Fermenting Clostridium Spp. Spore Contamination of Manchega Ewe Milk: Determination of Risk Factors and Characterization of Clostridium Population. Small Rumin. Res. 2013, 111, 120–128. [Google Scholar] [CrossRef]
  3. Cocolin, L.; Innocente, N.; Biasutti, M.; Comi, G. The Late Blowing in Cheese: A New Molecular Approach Based on PCR and DGGE to Study the Microbial Ecology of the Alteration Process. Int. J. Food Microbiol. 2004, 90, 83–91. [Google Scholar] [CrossRef]
  4. Cremonesi, P.; Vanoni, L.; Silvetti, T.; Morandi, S.; Brasca, M. Identification of Clostridium Beijerinckii, Cl. Butyricum, Cl. Sporogenes, Cl. Tyrobutyricum Isolated from Silage, Raw Milk and Hard Cheese by a Multiplex PCR Assay. J. Dairy Res. 2012, 79, 318–323. [Google Scholar] [CrossRef] [PubMed]
  5. Brändle, J.; Heinzle, L.; Fraberger, V.; Berta, J.; Zitz, U.; Schinkinger, M.; Stocker, W.; Kneifel, W.; Domig, K.J. Novel Approach to Enumerate Clostridial Endospores in Milk. Food Control 2018, 85, 318–326. [Google Scholar] [CrossRef]
  6. Reindl, A.; Dzieciol, M.; Hein, I.; Wagner, M.; Zangerl, P. Enumeration of Clostridia in Goat Milk Using an Optimized Membrane Filtration Technique. J. Dairy Sci. 2014, 97, 6036–6045. [Google Scholar] [CrossRef]
  7. Zucali, M.; Bava, L.; Colombini, S.; Brasca, M.; Decimo, M.; Morandi, S.; Tamburini, A.; Crovetto, G.M. Management Practices and Forage Quality Affecting the Contamination of Milk with Anaerobic Spore-Forming Bacteria. J. Sci. Food Agric. 2015, 95, 1294–1302. [Google Scholar] [CrossRef]
  8. Vissers, M.M.M.; Driehuis, F.; Te Giffel, M.C.; De Jong, P.; Lankveld, J.M.G. Concentrations of Butyric Acid Bacteria Spores in Silage and Relationships with Aerobic Deterioration. J. Dairy Sci. 2007, 90, 928–936. [Google Scholar] [CrossRef]
  9. Bassi, D.; Fontana, C.; Zucchelli, S.; Gazzola, S.; Cocconcelli, P.S. TaqMan Real Time-Quantitative PCR Targeting the Phosphotransacetylase Gene for Clostridium Tyrobutyricum Quantification in Animal Feed, Faeces, Milk and Cheese. Int. Dairy J. 2013, 33, 75–82. [Google Scholar] [CrossRef]
  10. McHugh, A.J.; Feehily, C.; Hill, C.; Cotter, P.D. Detection and Enumeration of Spore-Forming Bacteria in Powdered Dairy Products. Front. Microbiol. 2017, 8, 109. [Google Scholar] [CrossRef]
  11. Cecere, P.; Gatto, F.; Cortimiglia, C.; Bassi, D.; Lucchini, F.; Cocconcelli, P.S.; Pompa, P.P. Colorimetric Point-of-Care Detection of Clostridium Tyrobutyricum Spores in Milk Samples. Biosensors 2021, 11, 293. [Google Scholar] [CrossRef] [PubMed]
  12. Stadhouders, J.; van den Berg, G. Use of Lysozyme for the Prevention of Butyric Acid Fermentation in Gouda Cheese. Limited Effect of the Enzyme. In MILK the Vital Force, Presented at the XXII International Dairy Congress, The Hague, The Netherlands, 29 September–3 October 1986; Springer: Dordrecht, The Netherlands, 1986; p. 77. ISBN 978-94-009-3733-8. [Google Scholar]
  13. Publications Office of the European Union. CELEX1, Commission Regulation (EU) No 1129/2011 of 11 November 2011 Amending Annex II to Regulation (EC) No 1333/2008 of the European Parliament and of the Council by Establishing a Union List of Food Additives (Text with EEA Relevance). Available online: http://op.europa.eu/en/publication-detail/-/publication/28cb4a37-b40e-11e3-86f9-01aa75ed71a1 (accessed on 27 March 2023).
  14. Disciplinare Grana Padano, D.O.P. Available online: https://www.granapadano.it/it-it/disciplinare-grana-padano-dop.aspx (accessed on 8 February 2023).
  15. Disciplinare e Normative del Consorzio—Parmigiano Reggiano. Available online: https://www.parmigianoreggiano.com/it/consorzio-disciplinare-normative/ (accessed on 8 February 2023).
  16. Arnaboldi, S.; Benevenia, R.; Bertasi, B.; Galuppini, E.; Mangeri, L.; Tilola, M.; Bassi, D.; Cocconcelli, P.S.; Stroppa, A.; Varisco, G. Validation of a Real-Time PCR Method on Pta Gene for Clostridium Tyrobutyricum Quantification in Milk. Food Control 2021, 130, 108250. [Google Scholar] [CrossRef]
  17. Borreani, G.; Tabacco, E. Low Permeability to Oxygen of a New Barrier Film Prevents Butyric Acid Bacteria Spore Formation in Farm Corn Silage. J. Dairy Sci. 2008, 91, 4272–4281. [Google Scholar] [CrossRef]
  18. Colombari, G.; Borreani, G.; Crovetto, G.M. Effect of Ensiling Alfalfa at Low and High Dry Matter on Production of Milk Used to Make Grana Cheese. J. Dairy Sci. 2001, 84, 2494–2502. [Google Scholar] [CrossRef]
  19. Morandi, S.; Silvetti, T.; Brasca, M. Content and Spatial Distribution of Dairy-Related Clostridium Spores in Grana Padano Cheese during the Ripening Period. LWT 2022, 167, 113850. [Google Scholar] [CrossRef]
  20. Giraffa, G. The Microbiota of Grana Padano Cheese. A Review. Foods 2021, 10, 2632. [Google Scholar] [CrossRef]
  21. Brändle, J.; Domig, K.J.; Kneifel, W. Relevance and Analysis of Butyric Acid Producing Clostridia in Milk and Cheese. Food Control 2016, 67, 96–113. [Google Scholar] [CrossRef]
  22. Floris, I.; Martucci, F.; Romano, A.; Marello, G.; Ligotti, C.; Bianchi, D.M. Multiplex-PCR Detection of Clostridium Tyrobutyricum, Clostridium Butyricum, and Clostridium Sporogenes in Raw Milk for Cheesemaking. Life 2024, 14, 1093. [Google Scholar] [CrossRef]
  23. Burtscher, J.; Rudavsky, T.; Zitz, U.; Domig, K.J. Specificity of the AMP-6000 Method for Enumerating Clostridium Endospores in Milk. Foods 2024, 13, 1192. [Google Scholar] [CrossRef]
  24. SY-LAB—AMP-6000. Available online: https://microbiology.sylab.com/products/p/show/Product/product/amp-6000R.html (accessed on 26 April 2024).
  25. UNE EN ISO 7218:2008/A1:2013; Microbiology of Food and Animal Feeding Stuffs—General Requirements and Guidance for Microbiological Examinations—Amendment 1 (ISO 7218:2007/Amd 1:2013). ISO: Geneva, Switzerland, 2013. Available online: https://www.en-standard.eu/une-en-iso-7218-2008-a1-2013-microbiology-of-food-and-animal-feeding-stuffs-general-requirements-and-guidance-for-microbiological-examinations-amendment-1-iso-7218-2007-amd-1-2013/ (accessed on 24 April 2024).
  26. CLAL—Pagamento Latte-Qualita’. Available online: https://www.clal.it/?section=plq (accessed on 6 February 2023).
  27. IDF. Factsheet 22/2022: Enumeration of Butyric Acid Forming (Cheese Spoiling) Clostridia—Methodical Considerations. Available online: https://shop.fil-idf.org/products/idf-factsheet-22-2022-enumeration-of-butyric-acid-forming-cheese-spoiling-clostridia-methodical-considerations (accessed on 6 December 2022).
  28. Regulation (EC) No 853/2004 of the European Parliament and of the Council of 29 April 2004 Laying Down Specific Hygiene Rules for Food of Animal Origin; European Union: Brussels, Belgium, 2004; Volume 139.
  29. Bava, L.; Colombini, S.; Zucali, M.; Decimo, M.; Morandi, S.; Silvetti, T.; Brasca, M.; Tamburini, A.; Crovetto, G.M.; Sandrucci, A. Efficient Milking Hygiene Reduces Bacterial Spore Contamination in Milk. J. Dairy Res. 2017, 84, 322–328. [Google Scholar] [CrossRef] [PubMed]
  30. Burtscher, J.; Rudavsky, T.; Zitz, U.; Neubauer, V.; Domig, K.J. Importance of Pre-Milking Udder Hygiene to Reduce Transfer of Clostridial Spores from Teat Skin to Raw Milk. Microorganisms 2023, 11, 1337. [Google Scholar] [CrossRef] [PubMed]
  31. Lavilla, M.; De Luis, R.; Pérez, M.D.; Calvo, M.; Sánchez, L. Selection of High Affine Peptide Ligands for Detection of Clostridium Tyrobutyricum Spores. J. Microbiol. Methods 2009, 79, 214–219. [Google Scholar] [CrossRef] [PubMed]
  32. Doyle, C.J.; Gleeson, D.; Jordan, K.; Beresford, T.P.; Ross, R.P.; Fitzgerald, G.F.; Cotter, P.D. Anaerobic Sporeformers and Their Significance with Respect to Milk and Dairy Products. Int. J. Food Microbiol. 2015, 197, 77–87. [Google Scholar] [CrossRef]
  33. Nishihara, M.; Takahashi, H.; Sudo, T.; Kyoi, D.; Kawahara, T.; Ikeuchi, Y.; Fujita, T.; Kuda, T.; Kimura, B.; Yanahira, S. Multilocus Variable-Number of Tandem Repeat Analysis (MLVA) for Clostridium tyrobutyricum Strains Isolated from Cheese Production Environment. Int. J. Food Microbiol. 2014, 190, 61–65. [Google Scholar] [CrossRef]
  34. Prost, H.; Hugues, B. Dénombrement Des Micro-Organismes par la Technique du Nombre le Plus Probable (Indice NPP): Emploi Statistique de Cet Indice. In Annales de Falsifications de L’Expertise Chimique; La Société: Les Alluets-le-Roi, France, 1982; Volume 75, pp. 185–207. [Google Scholar]
  35. McBride, G.B.; McWhirter, J.L.; Dalgety, M.H. Uncertainty in Most Probable Number Calculations for Microbiological Assays. J. AOAC Int. 2003, 86, 1084–1088. [Google Scholar]
  36. Weenk, G.H. Chapter 1 Microbiological Assessment of Culture Media: Comparison and Statistical Evaluation of Methods. Prog. Ind. Microbiol. 2003, 37, 1–23. [Google Scholar]
  37. Skalar|Laboratory Automation. Available online: https://www.skalar.com/ (accessed on 20 December 2023).
  38. Niemelä, S.I. Uncertainty of Quantitative Determinations Derived by Cultivation of Microorganisms; MIKES Publication; Centre of Metrology and Accreditation (MIKES): Helsinki, Finland, 2003; ISBN 978-952-5209-76-1. [Google Scholar]
  39. «UNI CEI EN ISO/IEC 17025:2018 - UNI Ente Italiano di Normazione». Available online: https://store.uni.com/uni-cei-en-iso-iec-17025-2018 (accessed on 5 November 2024).
  40. Heyndrickx, M. The Importance of Endospore-Forming Bacteria Originating from Soil for Contamination of Industrial Food Processing. Appl. Environ. Soil Sci. 2011, 2011, e561975. [Google Scholar] [CrossRef]
  41. Esame Routine Sporigeni—Istituto Zooprofilattico Sperimentale Della Lombardia e Dell’Emilia Romagna. Available online: https://www.izsler.it/carta_servizi/esame-routine-sporigeni-2/ (accessed on 8 March 2024).
  42. CLAL TESEO—UE-27: Analisi Delle Aziende Agricole Europee. Available online: https://teseo.clal.it/?section=aziende-agricole (accessed on 4 October 2024).
  43. I Risultati: Produzione e Consumo. Grana Padano. Available online: https://www.granapadano.it/it-it/il-consorzio-di-tutela/i-risultati-produzione-e-consumo/ (accessed on 7 October 2024).
  44. Regolamento (UE) 2017/625 del Parlamento Europeo e del Consiglio, del 15 Marzo 2017, Relativo ai Controlli Ufficiali e Alle Altre Attività Ufficiali Effettuati per Garantire L’Applicazione Della Legislazione Sugli Alimenti e sui Mangimi, Delle Norme Sulla Salute e Sul Benessere Degli Animali, Sulla Sanità Delle Piante Nonché sui Prodotti Fitosanitari, Recante Modifica dei Regolamenti (CE) n. 999/2001, (CE) n. 396/2005, (CE) n. 1069/2009, (CE) n. 1107/2009, (UE) n. 1151/2012, (UE) n. 652/2014, (UE) 2016/429 e (UE) 2016/2031 del Parlamento Europeo e del Consiglio, dei Regolamenti (CE) n. 1/2005 e (CE) n. 1099/2009 del Consiglio e Delle Direttive 98/58/CE, 1999/74/CE, 2007/43/CE, 2008/119/CE e 2008/120/CE del Consiglio, e che Abroga i Regolamenti (CE) n. 854/2004 e (CE) n. 882/2004 del Parlamento Europeo e del Consiglio, le Direttive 89/608/CEE, 89/662/CEE, 90/425/CEE, 91/496/CEE, 96/23/CE, 96/93/CE e 97/78/CE del Consiglio e la Decisione 92/438/CEE del Consiglio (Regolamento sui Controlli Ufficiali) (Testo Rilevante ai fini del SEE); European Union: Brussels, Belgium, 2017; Volume 095.
  45. Neviani, E.; Bottari, B.; Lazzi, C.; Gatti, M. New Developments in the Study of the Microbiota of Raw-Milk, Long-Ripened Cheeses by Molecular Methods: The Case of Grana Padano and Parmigiano Reggiano. Front. Microbiol. 2013, 4, 36. [Google Scholar] [CrossRef]
  46. Bittante, G.; Cologna, N.; Cecchinato, A.; De Marchi, M.; Penasa, M.; Tiezzi, F.; Endrizzi, I.; Gasperi, F. Monitoring of Sensory Attributes Used in the Quality Payment System of Trentingrana Cheese. J. Dairy Sci. 2011, 94, 5699–5709. [Google Scholar] [CrossRef]
  47. von Keyserlingk, M.A.G.; Rushen, J.; de Passillé, A.M.; Weary, D.M. Invited Review: The Welfare of Dairy Cattle—Key Concepts and the Role of Science. J. Dairy Sci. 2009, 92, 4101–4111. [Google Scholar] [CrossRef]
  48. Council Directive 98/58/EC of 20 July 1998 Concerning the Protection of Animals Kept for Farming Purposes. 2019. Available online: https://www.legislation.gov.uk/eudr/1998/58 (accessed on 5 November 2024).
  49. Bermúdez, J.; González, M.J.; Olivera, J.A.; Burgueño, J.A.; Juliano, P.; Fox, E.M.; Reginensi, S.M. Seasonal Occurrence and Molecular Diversity of Clostridia Species Spores along Cheesemaking Streams of 5 Commercial Dairy Plants. J. Dairy Sci. 2016, 99, 3358–3366. [Google Scholar] [CrossRef] [PubMed]
  50. Calamari, L.; Morera, P.; Bani, P.; Minuti, A.; Basiricò, L.; Vitali, A.; Bernabucci, U. Effect of Hot Season on Blood Parameters, Fecal Fermentative Parameters, and Occurrence of Clostridium Tyrobutyricum Spores in Feces of Lactating Dairy Cows. J. Dairy Sci. 2018, 101, 4437–4447. [Google Scholar] [CrossRef] [PubMed]
  51. Bertocchi, L.; Vitali, A.; Lacetera, N.; Nardone, A.; Varisco, G.; Bernabucci, U. Seasonal Variations in the Composition of Holstein Cow’s Milk and Temperature–Humidity Index Relationship. Animal 2014, 8, 667–674. [Google Scholar] [CrossRef]
Figure 1. The 20-year trends in the number of samples delivered to and analyzed at IZSLER. Panels (A,B) display the annual demand for the S2 and S3 protocols, respectively. Panels (C,D) illustrate the seasonal demand for the same protocols.
Figure 1. The 20-year trends in the number of samples delivered to and analyzed at IZSLER. Panels (A,B) display the annual demand for the S2 and S3 protocols, respectively. Panels (C,D) illustrate the seasonal demand for the same protocols.
Foods 13 03638 g001
Figure 2. The 20-year annual trend in positive samples from S2 and S3 analyses. Panel (A) compares total positive samples (Pos tot) for both protocols. Panels (B,C) show trends in weak (Pos) and strong (Pos+) positive samples identified with the S2 and S3, respectively. Asterisks (*) denote years with significant differences from 2004 in the relative frequency of positive samples.
Figure 2. The 20-year annual trend in positive samples from S2 and S3 analyses. Panel (A) compares total positive samples (Pos tot) for both protocols. Panels (B,C) show trends in weak (Pos) and strong (Pos+) positive samples identified with the S2 and S3, respectively. Asterisks (*) denote years with significant differences from 2004 in the relative frequency of positive samples.
Foods 13 03638 g002
Figure 3. Seasonal distribution of positive samples using the S2 and S3 protocols. Trends for total (Pos tot), weak (Pos), and strong (Pos+) positive samples across seasons are displayed in the panels (A,C,D) for the S2 protocol and (B,E,F) for the S3 protocol, respectively.
Figure 3. Seasonal distribution of positive samples using the S2 and S3 protocols. Trends for total (Pos tot), weak (Pos), and strong (Pos+) positive samples across seasons are displayed in the panels (A,C,D) for the S2 protocol and (B,E,F) for the S3 protocol, respectively.
Foods 13 03638 g003
Figure 4. The 20-year annual trend in positive samples from S2 (A) and S3 (B) analyses. Total (Pos tot), weak (Pos), and strong (Pos+) positive samples are reported as blue bars, yellow lines, and orange bars, respectively.
Figure 4. The 20-year annual trend in positive samples from S2 (A) and S3 (B) analyses. Total (Pos tot), weak (Pos), and strong (Pos+) positive samples are reported as blue bars, yellow lines, and orange bars, respectively.
Foods 13 03638 g004
Table 1. Bovine milk samples submitted for Clostridium spore analysis (2004–2023). The table summarizes all bovine milk samples delivered to IZSLER by private owners between 2004 and 2023 to assess the presence or absence of Clostridium spores through S2 (A) and S3 (B) methods. Samples not relevant to this study, such as milk from non-bovine species or milk that had undergone processes affecting spore survival or germination were excluded.
Table 1. Bovine milk samples submitted for Clostridium spore analysis (2004–2023). The table summarizes all bovine milk samples delivered to IZSLER by private owners between 2004 and 2023 to assess the presence or absence of Clostridium spores through S2 (A) and S3 (B) methods. Samples not relevant to this study, such as milk from non-bovine species or milk that had undergone processes affecting spore survival or germination were excluded.
S2_2004–2023 (A)
Yearn° Samples Receivedn° Samples Excludedn° Samples Included
202342,8803942,841
202245,0235744,966
202147,0034046,963
202041,6343041,604
201950,9372750,910
201848,9824048,942
201749,3385249,286
201648,8722848,844
201549,693649,687
201451,1703251,138
201352,2505052,200
201253,3257353,252
201152,17814652,032
201049,76810949,659
200949,9967249,924
200849,7976549,732
200752,99320852,785
200661,26425261,012
200559,56218159,381
200462,56311562,448
S3_2004–2023 (B)
Yearn° Samples Receivedn° Samples Excludedn° Samples Included
20233020972923
2022232602326
2021289212891
20203035203015
201946713554316
201851313804751
201758244105414
201630212262795
20152203232180
20142116201996
20131744131731
2012201912018
20112031232008
20102244102234
2009222832223
2008199851993
2007298862982
20063599193580
2005322623224
2004287512874
Table 2. Territorial distribution of raw bovine milk samples and analysis with the S2 (A) and S3 (B) protocols. The provinces included are BS (Brescia), BG (Bergamo), CR (Cremona), LO (Lodi), MI (Milano), MN (Mantua), VR (Verona), SO (Sondrio), PC (Piacenza), LC (Lecco), and PV (Pavia) with the right-hand column displaying the prevalence percentage that the provinces considered cover.
Table 2. Territorial distribution of raw bovine milk samples and analysis with the S2 (A) and S3 (B) protocols. The provinces included are BS (Brescia), BG (Bergamo), CR (Cremona), LO (Lodi), MI (Milano), MN (Mantua), VR (Verona), SO (Sondrio), PC (Piacenza), LC (Lecco), and PV (Pavia) with the right-hand column displaying the prevalence percentage that the provinces considered cover.
S2_2004–2023 (A)
YearSampling Provinces% Samples by OriginYear
2023BG, BS, CR, LO, MI, MN, VR97.90%2023
2022BG, BS, CR, LO, MI, MN, VR97.87%2022
2021BG, BS, CR, LO, MI, MN, VR97.63%2021
2020BG, BS, CR, LO, MI, MN, VR97.76%2020
2019BG, BS, CR, LO, MI, MN, VR98.21%2019
2018BG, BS, CR, LO, MI, MN, VR97.27%2018
2017BG, BS, CR, LO, MI, MN, VR97.35%2017
2016BG, BS, CR, LO, MI, MN, VR97.46%2016
2015BG, BS, CR, MN, VR97.55%2015
2014BG, BS, CR, MN, VR97.56%2014
2013BG, BS, CR, MN, VR97.86%2013
2012BG, BS, CR, MN, VR97.75%2012
2011BG, BS, CR, MN, VR97.77%2011
2010BG, BS, CR, MN, VR97.52%2010
2009BG, BS, CR, MN, VR97.54%2009
2008BG, BS, CR, MN, VR96.96%2008
2007BG, BS, CR, MN, VR96.12%2007
2006BG, BS, CR, LO, MI, MN, VR97.10%2006
2005BG, BS, CR, MN, VR98.16%2005
2004BG, BS, CR, MN, SO, VR98.16%2004
S3_2004–2023 (B)
YearSampling Provinces% Samples by Origin
2023BG, BS, PC97.43%
2022BG, BS,96.82%
2021BG, BS, MN97.23%
2020BG, BS, CR96.62%
2019BG, BS, LC, MI, MN97.38%
2018BG, BS, LC, MI97.20%
2017BG, BS, LC, MI97.8%
2016BG, BS, LC, MI, MN97.46%
2015BG, BS96.93%
2014BG, BS, CR96.40%
2013BG, BS, CR, MN95.96%
2012BG, BS, MN95.14%
2011BG, BS, MN, SO97.86%
2010BG, BS, MN, SO96.96%
2009BG, BS, SO96.85%
2008BG, BS, CR, SO98.60%
2007BG, BS, CR, SO98.56%
2006BG, BS, CR, SO98.16%
2005BG, BS, CR, SO97.08%
2004BG, BS, CR, PV96.48%
Table 3. Clostridium-positive samples following S2 (A) and S3 (B) analyses. The number of positives (Ps) for each dilution of the sample, MPN values (spore/L), and 95% confidence intervals are reported in A and B for samples processed with the S2 and S3 protocols, respectively. The S2 protocol includes two sample dilutions (1:10 and 1:100), while the S3 protocol also considers an undiluted sample.
Table 3. Clostridium-positive samples following S2 (A) and S3 (B) analyses. The number of positives (Ps) for each dilution of the sample, MPN values (spore/L), and 95% confidence intervals are reported in A and B for samples processed with the S2 and S3 protocols, respectively. The S2 protocol includes two sample dilutions (1:10 and 1:100), while the S3 protocol also considers an undiluted sample.
S2
PMPN
(Spore/L)
Confidence Intervals (Lower-Upper)
1 mL0.1 mL
00<308n.s.
0130866–1438
1036378–1697
11742159–3470
121140244–5327
20966207–4514
211508323–7047
222099449–9814
303015645–14,094
3149241053–23,015
3287041862–40,685
33>11,000n.s.
S3
PMPN
(Spore/L)
Confidence Intervals (Lower-Upper)
10 mL1 mL0.1 mL
000<310–110
010314–230
100365–270
1017217–300
1107418–310
12011035–370
2009221–400
20114042–490
21015043–500
21120070–600
22021071–620
30023055–970
30139093–1600
31043095–1900
311750190–3000
3121200360–3700
320930220–4000
3211500440–5100
3222200720–6400
3302400560–10,000
3314600960–22,000
33211,0002500–49,000
333>11,0003600–n.s.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Guaita, A.; Gambi, L.; Baresi, P.; Paterlini, F.; Bolzoni, G.; Zanardi, G.; Daminelli, P. Spore-Forming Clostridia in Raw Cow Milk from Northern Italy: A Trend Analysis over the Past 20 Years. Foods 2024, 13, 3638. https://doi.org/10.3390/foods13223638

AMA Style

Guaita A, Gambi L, Baresi P, Paterlini F, Bolzoni G, Zanardi G, Daminelli P. Spore-Forming Clostridia in Raw Cow Milk from Northern Italy: A Trend Analysis over the Past 20 Years. Foods. 2024; 13(22):3638. https://doi.org/10.3390/foods13223638

Chicago/Turabian Style

Guaita, Arianna, Lorenzo Gambi, Pierluigi Baresi, Franco Paterlini, Giuseppe Bolzoni, Giorgio Zanardi, and Paolo Daminelli. 2024. "Spore-Forming Clostridia in Raw Cow Milk from Northern Italy: A Trend Analysis over the Past 20 Years" Foods 13, no. 22: 3638. https://doi.org/10.3390/foods13223638

APA Style

Guaita, A., Gambi, L., Baresi, P., Paterlini, F., Bolzoni, G., Zanardi, G., & Daminelli, P. (2024). Spore-Forming Clostridia in Raw Cow Milk from Northern Italy: A Trend Analysis over the Past 20 Years. Foods, 13(22), 3638. https://doi.org/10.3390/foods13223638

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop