Next Article in Journal
Microencapsulates of Blue Maize Polyphenolics as a Promising Ingredient in the Food and Pharmaceutical Industry: Characterization, Antioxidant Properties, and In Vitro-Simulated Digestion
Previous Article in Journal
Effect of a Dairy Cow’s Feeding System on the Flavor of Raw Milk: Indoor Feeding or Grazing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Safety Evaluation and Probabilistic Health Risk Assessment of Cow Milk Produced in Northern Italy According to Dioxins and PCBs Contamination Levels

1
Department of Food and Drug, University of Parma, 43124 Parma, Italy
2
Istituto Zooprofilattico Sperimentale della Lombardia e dell’Emilia Romagna “Bruno Ubertini”, 25124 Brescia, Italy
3
Chemical Department, Istituto Zooprofilattico Sperimentale della Lombardia e dell’Emilia-Romagna “Bruno Ubertini”, 40127 Bologna, Italy
4
Unità Organizzativa Veterinaria, Direzione Generale Welfare Regione Lombardia, 20124 Milano, Italy
5
Area Sanità Veterinaria e Igiene degli Alimenti, Settore Prevenzione Collettiva e Sanità Pubblica, Direzione Generale Cura della Persona, Salute e Welfare, Regione Emilia Romagna, 40127 Bologna, Italy
*
Author to whom correspondence should be addressed.
Foods 2023, 12(9), 1869; https://doi.org/10.3390/foods12091869
Submission received: 3 April 2023 / Revised: 27 April 2023 / Accepted: 29 April 2023 / Published: 30 April 2023
(This article belongs to the Section Food Quality and Safety)

Abstract

:
Contamination levels of dioxins and polychlorinated biphenyls (PCBs) were monitored over 2018–2021 in 214 bovine milk samples from farms located in two regions in northern Italy (Lombardy and Emilia-Romagna). The average concentrations of the sum of dioxins and dioxin-like PCBs (0.78 ± 0.55 pg TEQ/g fat) and six non-dioxin-like PCBs (6.55 ± 2.24 ng/g fat) were largely below the maximum, and action limits established at European level, confirming a decreasing trend observed both locally and across Europe in recent years. The impact of contamination levels on chronic dietary exposure of the Italian population to dioxins and PCBs was found to be highly variable based on the type of cow milk (skimmed, semi-skimmed, or whole-fat milk) and the population age group considered. Indeed, a first-tier screening of the potential exposure via determinist methods allowed for the identification of the youngest population as the group with the worst risk profile. The refinement of exposure assessment via Monte Carlo probabilistic methods suggested that, at the less pessimistic middle-bound simulation scenario, infants, toddlers, and children consuming whole cow milk may be exposed to dioxins and PCBs levels above the toxicological reference values with a probability of 76, 56, and 22%, respectively.

1. Introduction

Polychlorinated dibenzo-p-dioxins and furans (PCDD/Fs, also known as dioxins) and polychlorinated biphenyls (PCBs) are persistent organic pollutants widely distributed in the environment. Even though they include a wide range of chemical compounds, they all have low degradation potentials and high lipophilic character. Furthermore, PCDD/Fs and PCBs have similar chronic toxic effects on animals and humans, potentially causing developmental and reproductive disorders, endocrine and immune system impairment, teratogenesis, and tumorigenesis [1,2,3,4].
The most toxic congeners are 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) and 1,2,3,7,8-pentachlorodibenzo-p-dioxin (PCDD) [5], with an estimated half-life in the human body of approximately 10 years after the initial exposure [2,6].
Due to their high resistance to metabolic degradation and lipophilic behavior, these compounds easily accumulate in food and feed chains [7]. Milk and derived products were reported as the first dietary contributors to the total body burden of dioxins and PCBs for infants, toddlers, and children (28–50%), but also impactful for older European population groups (7–25%) [3].
The presence of dioxins and PCBs in the fat fraction of cow milk in considerable amounts is the consequence of the active excretion of these compounds by the animals, which use this mechanism as an elimination route [8]. Sources of dioxins and PCBs for cattle are pasture grass for free-range animals and mainly silage feedstuffs and feed additives for intensively farmed ones [9,10,11]. More specifically, certain types of feed such as citrus pulp, clay minerals, breadcrumbs, potato peelings, corn, sugar beet pulp, and certain fats and oils have been found to contain high levels of PCDD/Fs and PCBs [12,13,14,15,16]. Dried materials such as corn, sugar beet pulp, grass, alfalfa, and breadcrumbs are particularly susceptible to contamination due to the possible formation of PCDD/Fs if they are not dried properly [17]. In turn, the source of feed contamination may be ascribed to emissions from both industrial and non-industrial sources (e.g., coal burning in furnaces, domestic waste burning, and car exhausts) [18]. In this setting, it has been widely demonstrated that concentrations in animal-derived products are frequently higher than in the feed material carrying them [18,19,20]. Many factors related to the physiology, biology, and management of animals, including age, grazing habits, breeding characteristics, lactation stage, body weight, body fat content, and milkfat production, may contribute to the presence of dioxins and PCBs in animal tissues and milk [21,22,23,24]. However, there is not enough conclusive evidence to determine the specific influence of these factors since most studies do not focus on investigating them [24]. Conversely, greater emphasis has been placed on how the contamination of milk with dioxins and PCBs can vary depending on geographic location and season. Even small, periodic changes in contamination levels have been observed, but these were mostly attributed to fluctuations in the concentrations of contaminants in the feed consumed by the cattle, as well as differences in pollutant exposure between indoor and outdoor environments during different seasons [25,26,27].
The milk and dairy sectors are an important part of the Italian economy, with an annual production of thirteen million tons of bovine milk, mainly intended for fresh consumption or cheese making [28]. Specifically, approximately 20,000 cow farms are spread across the northern Italian regions of Lombardy and Emilia-Romagna, which alone account for more than 60% of national milk production [29]. The Italian population, therefore, has a high consumption of milk and dairy products, and exposure scenarios to dioxins and PCBs through the consumption of these food items may be of particular concern compared to other European citizens.
To reduce the concentrations of dioxins and PCBs along the whole food chain, risk management measures have been implemented in the European Union. First, maximum limits (MLs) and action limits (ALs) were established to prevent contaminated foodstuffs and feedstuffs from being placed on the market and to improve the identification and removal of pollution sources. Although, in the last few years, these measures have led to a general reduction in human assimilation of dioxins through the diet [3], it should be noted that MLs have been set according to the so-called ‘strict but feasible’ principle, i.e., they have been established at the 90th percentile of their frequency distribution in each specific food item without following toxicologically driven criteria [3,30]. Moreover, despite the reduction in dietary exposure to dioxins and PCBs, the European Food Safety Authority (EFSA) has recently estimated that the tolerable weekly intake (TWI) of these contaminants continues to be exceeded for most of the European population and, therefore, the Authority has decided to drastically reduce the TWI of the sum of PCDD/Fs and dioxin-like (dl)-PCB to 2 pg toxic equivalent (TEQ)/kg bw/week [31]. As a result, the European Commission lowered the MLs for several food commodities in 2022. For milk and dairy, MLs of 2 and 4 pg TEQ/g fat of PCDD/Fs and PCDD/Fs + dl-PCB, respectively, and of 40 ng/g fat for non-dioxin-like (ndl)-PCBs are currently in force [32].
In light of new TWI and MLs limits, as well as changes in the dietary habits of the Italian population, it is therefore important to pursue risk management strategies through the ongoing monitoring of dioxin and PCB contamination levels in food and feed to measure the public health impact of dietary exposure to these contaminants and to identify those segments of the population that are at increased risk of exposure and potentially long-term toxic effects. In this context, probabilistic risk assessment might provide more confident and realistic representations of the actual risk levels to which consumers are exposed, accounting for the variability and uncertainty of exposure estimates.
In this work, the results of the 2018–2021 monitoring plan on dioxin and PCB residues carried out in the northern Italian regions of Lombardy and Emilia Romagna are presented, with the aim of evaluating the occurrence, spatial distributions, and temporal variations of these contaminants in bovine milk produced in this territory. A second-level objective was to perform an exposure assessment of dioxins and PCBs intake through the consumption of milk with different fat contents by different population age groups in Italy and verify to what extent the consumption of milk contributes to the toxicological reference values using both deterministic and probabilistic (Monte Carlo method) approaches, providing a realistic scenario of the possible risks.

2. Materials and Methods

2.1. Milk Samples

Milk samples were collected from 2018 to 2021 under the framework of the Italian residue monitoring plan (RMP). In those years, the Italian RMP was structured according to Council Directive 96/23/EC [33] and Commission Decision 97/747/EC [34], as well as Regulation (EU) 625/2017/EU [35]. A total of 214 bovine milk samples were obtained from 187 dairy farms in the Lombardy and Emilia-Romagna regions of Northern Italy (Figure 1). Milk samples were collected from bulk farm tanks by the official competent authorities according to Regulation (EU) 644/2017/EU [36]. Specifically, bulk milk was thoroughly mixed in order to ensure and homogeneous distribution of the contaminants, and an aggregate sample of 1 L was obtained from each farm. Once collected, the samples were stored at +4 °C and delivered within 24 h to the official local laboratory of the Istituto Zooprofilattico Sperimentale della Lombardia e dell’Emilia-Romagna, certified under ISO/IEC 17025:2017 and accredited for quantitative determination of PCDD/Fs (2,3,7,8-tetracholodibezodioxin (TCDD), 1,2,3,7,8-pentachlorodibenzodioxin (PeCDD), 1,2,3,4,7,8-hexachlorodibenzodioxin (HxCDD), 1,2,3,6,7,8-HxCDD, 1,2,3,7,8,9-HxCDD, 1,2,3,4,6,7,8-heptaclorodibenzodioxin (HpCDD), 1,2,3,4,6,7,8,9-octochlorodibenzodioxin (OCDD), 2,3,7,8-tetracholodibezofuran (TCDF), 2,3,4,7,8-pentacholodibezofuran (PeCDF), 1,2,3,7,8-PeCDF, 2,3,4,6,7,8-hexacholodibezofuran HxCDF, 1,2,3,7,8,9-HxCDF, 1,2,3,6,7,8-HxCDF, 1,2,3,4,7,8-HxCDF, 1,2,3,4,7,8,9-heptacholodibezofuran (HpCDF), 1,2,3,4,6,7,8-HpCDF, 1,2,3,4,6,7,8,9-octacholodibezofuran (OCDF), dl-PCBs (PCB 77, 81, 105, 114, 118, 123, 126, 156, 157, 167, 169, 189) and ndl-PCBs (PCB 28, 52, 101, 138, 153, 180), using high-resolution gas chromatography coupled with high-resolution mass spectrometry (HRGC-HRMS).

2.2. Reagents and Standards

Ethyl-acetate, toluene, and n-hexane were purchased from Carlo Erba Reagents (Milan, Italy); nonane from Promochem (LGC Standards, Teddington, UK) and dichloromethane from Romil Ltd. (Cambridge, UK). All solvents were picograde.
Pre-packed multi-layer silica, alumina, and carbon columns were produced by FMS (Fluid Management System, Billerica, MA, USA). The 13C-labeled recovery, clean-up, and standard injection solutions were provided by CIL (Cambridge Isotope Laboratories, Andover, MA, USA). For PCDD/Fs, EDF-9999 Method 1613 calibration solutions (CS1-CS5) were used. For PCB calibration, an in-house curve was prepared using PCB MIX-75 (Dr. Ehrenstorfer, Augsburg, Germany) and 13C-labeled solutions EC-4995 and EC-4978 (Cambridge Isotope Laboratories, Andover, MA, USA).
All solvents and reagents used for the analyses were tested to ensure the absence of contaminants at the levels of interest (i.e., below one-fifth of the limit of quantification (LOQ) for PCDD/Fs and below one-tenth of the LOQ for dl-PCB and ndl-PCB).

2.3. Analytical Method

The methods used for the determination of PCDD/Fs and PCBs in milk samples were previously described by Lorenzi et al. [20,37]. Briefly, milk samples were homogenized and freeze-dried (Freeze Dryer Martin Christ Gefriertrocknungsanlagen, Osterode am Harz, Germany). Following further homogenization, a portion of each sample was tested with the Soxhlet method to determine the lipid content, while 8–10 g of the samples were mixed with diatomaceous earth and spiked with the 13C-labeled internal standard EDF-8999 (Cambridge Isotope Laboratories, Andover, MA, USA) and the 13C-labeled internal standard EC-4995 (Cambridge Isotope Laboratories, Andover, MA, USA). An accelerated solvent extractor (ASE® 300, Sunnyvale, Dionex, CA, USA) was used for fat extraction with toluene (2 cycles of 20 min at +135 °C and 1500 psi). The solvent was filtered using anhydrous sodium sulfate and evaporated at +50 °C with a rotatory evaporator. The lipid extracts were dried overnight at 70 °C ± 5 °C. Fat content was determined gravimetrically and compared to the value obtained from the Soxhlet method to confirm ASE extraction efficiency.
Lipid extracts were solubilized in 5 mL of hexane/dichloromethane solution (1:1, v/v), spiked with the standard clean-up solution EC-4978 (Cambridge Isotope Laboratories, Andover, MA, USA), containing three 13C-labeled PCB congeners, and diluted with 20 mL of hexane. The diluted extracts were purified using pre-packed silica columns, acidified with sulfuric acid, and eluted with n-hexane. The purification fractions were concentrated to 0.5 mL in the TurboVap evaporator (Zymark, Hopkinton, MA, USA) and loaded into the Power-Prep system (Fluid Management System, Watertown, MA, USA) equipped with silica, alumina, and carbon columns. Toluene (50 mL) was used to elute PCDD/Fs from the carbon column, while n-hexane and a mixture of hexane/dichloromethane solution (9:1, v/v) were used for PCB elution from the alumina column. The final extracts were dried using the TurboVap evaporator and a vacuum concentrator (Genevac, Ipswich, UK).
The PCDD/F fraction was dissolved in 10 μL of 1:100 ED-2521 injection solution (Cambridge Isotope Laboratories, Andover, MA, USA), and the PCB fraction was dissolved in 20 μL of 1:50 EC-4979 injection solution (Cambridge Isotope Laboratories, Andover, MA, USA). HRGC-HRMS analysis was carried out using two TRACE GC ULTRA gas chromatography coupled with a DFS (double focusing system) high-resolution magnetic scan system (Thermo Fisher Scientific, Waltham, MA, USA) or with AutoSpec high-resolution mass spectrometer (Micromass/Waters, Manchester, UK) [24]. A DB5 MS capillary column (60 m × 0.25 mm, 0.25 µm, Thermo Fisher Scientific, Waltham, MA, USA) was used for PCDD/F separation, and a TR-PCB 8MS capillary column (50 m × 0.25 mm, 0.25 µm, Thermo Fisher Scientific, Waltham, MA, USA) for PCB separation. The mass spectrometer operated in selected ion monitoring mode with a mass resolution of over 10,000.
Results were expressed in pg/g fat for PCDD/Fs and dl-PCBs and ng/g fat for ndl-PCBs. Toxic equivalent values (TEQ) were calculated using the 2005 World Health Organization Toxic Equivalency Factors (WHO-TEFs) [5].

2.4. Quality Control

The laboratory participated in different international proficiency tests and inter-laboratory studies, obtaining satisfactory results.
To ensure the absence of cross-contamination and interferences, blank samples were processed with each batch of five samples. For every ten samples, a fortified sample was processed to ensure analyte recovery. The stability of response factors for all congeners was checked daily using calibration verification standard solutions for PCDD/Fs and PCBs and assessing the acceptability of deviations from the calibration curves [38,39]. The recovery values of the internal standard congeners were determined for each sample.
Duplicate analyses were performed on each sample, and the compliance of the absolute difference between the two measurements with the repeatability limits, obtained during the validation study, was used to estimate the overall repeatability of the method.
LOQ and recovery values are reported in Supplementary Table S1.

2.5. Summary and Statistics

Contamination levels of dioxins and dl-PCBs were reported both as actual mass fraction values (pg/g fat) and toxic equivalency values (pg TEQ/g fat), while those of ndl-PCBs as actual mass fraction values (ng/g fat).
Non-detected values, i.e., values below the LOQ, were substituted with the LOQ (upper bound concentrations, UB), ½ LOQ (middle bound concentrations, MB), or 0 (lower bound concentrations, LB). Univariate descriptive statistics were applied to the UB-data matrix to identify potential differences among the tested contaminants. The Student’s T test was employed to evaluate differences between the two Italian regions from which samples were collected (Lombardy and Emilia-Romagna). The Analysis of Variance (ANOVA) followed by Tukey’s post hoc test was used to identify the provinces of each Italian region characterized by the highest concentrations of contaminants and to explore significant differences over the 2018–2021 sampling years. A Spearman’s correlation analysis was performed to investigate potential links between the fat content and dioxin concentrations in milk samples. PCDD/F + dl-PCB concentrations (pg/g fat) of each sample were multiplied by the relative fat content (g fat/100 g) to calculate the overall PCDD/Fs + dl-PCBs in 100 g of milk, and then these results were correlated with the respective lipid content. Correlation coefficients (r) higher or lower than 0.6 were considered indicative of the existence of a positive or negative correlation, respectively. The statistical significance was set at p ≤ 0.05 for all tests.
Statistical analysis was performed using OriginPro 2021 software (v. 9.8.0.200, OriginLab, Northampton, MA, USA).

2.6. Exposure Assessment and Risk Characterisation

2.6.1. Screening via Deterministic Approach

The estimated weekly intakes (EWIs) of PCDD/Fs + dl-PCBs through milk consumption were calculated using the measured PCDD/Fs + dl-PCBs mean (UB) concentrations, milkfat concentrations, and mean or 95th percentile chronic dietary consumption data of cow milk by the Italian population.
Consumption data of cow milk (consumers only) were retrieved from a more recent Italian survey on food consumption SCAI IV included within the EFSA Comprehensive Food Consumption Database [40,41,42]. Infants (up to 12 months old), toddlers (13–36 months old), children (37 months–9 years old), adolescents (10–17 years old), adults (18–65 years old), and the elderly (>65 years old) were the six Italian population age groups considered.
The cow milk collected in this study was supposed to be destined for the production and commercialization of milk with different fat percentages. Therefore, dietary consumption data of skimmed, semi-skimmed, and whole cow milk were retrieved from Level 6 of the FoodEx2 hierarchical classification system included in the EFSA Database and used to simulate three different chronic exposure scenarios per population group. For the calculation of EWIs, skimmed, semi-skimmed, and whole cow milk were supposed to contain 0.3, 1.8, and 3.5 g lipids/100 g, respectively, by using Equation (1):
EWI (pg TEQ/kg bw/week) = C × M × F ÷ 100,
where C is the mean concentration of PCDD/Fs + dl-PCBs (pg TEQ/g fat); M is the mean or 95th percentile chronic weekly consumption of skimmed, semi-skimmed, or whole cow milk (g/kg bw/week); F is the fat content of skimmed, semi-skimmed, or whole cow milk samples (g fat/100 g).
Finally, the risk was characterized by comparing the calculated EWI with the TWI of PCDD/Fs + dl-PCBs (equal to 2 pg/kg bw/week) by using Equation (2):
TWI (%) = (EWI ÷TWI) × 100.

2.6.2. Refinement via Probabilistic Methods (Monte Carlo Simulations)

The Monte Carlo (MC) probabilistic method was applied to simulate the distribution probabilities of EWI values and the related %TWI only for those deterministic scenarios indicating potential causes of concern, as suggested by EFSA [43]. Rather than using single point estimates, the distribution frequency of both the PCDD/Fs + dl-PCB and chronic consumption data of cow milk (i.e., assumptions) were considered. In this way, information concerning uncertainty, variability, and probability in health risk evaluation was obtained.
First, the best-fit probability distributions for the PCDD/Fs + dl-PCB contamination levels were tested using the modified Kolmogorov–Smirnov goodness-of-fit test (p > 0.05). In this case, the LB, MB, and UB concentrations were individually used to account for uncertainty embedded within the data. Based on the results, the LB, MB, and UB concentrations were fitted to the lognormal distribution. No distribution frequencies were set for milkfat, and fixed values at 1.8 and 3.5 g/100 g for semi-skimmed and whole cow milk were used. The overall distribution frequency of dietary consumption data for milk (expressed as g/kg bw/week) was supposed to be lognormal based on the literature data [44,45].
After that, MCSs were run randomly and repeatedly using 50,000 iterations for each population subgroup using Equations (1) and (2) so as to forecast a set of theoretically stable distributions of PCDD/F + dl-PCB EWIs and relative contributions to TWI (%), as well as their probabilities at different risk levels (10–90th percentiles).
All the analyses were performed using the Crystal Ball software (v. 11.1.2.4, Oracle©, Austin, TX, USA).

3. Results and Discussion

3.1. Spatial and Temporal Distribution of Dioxins and PCBs in Milk

The average concentrations of dioxins and PCBs measured in the 214 cow milk samples are reported in Table 1, where they were statistically summarized by the Italian region (i.e., Lombardy vs. Emilia-Romagna) and the year of sampling (i.e., 2018 vs. 2019 vs. 2020 vs. 2021). Individual concentrations of each of the 35 measured PCDD/F, dl-PCB, and ndl-PCB congeners are provided in Supplementary Tables S2 and S3. Regardless of both the region and the year, the average concentrations of PCDD/Fs, dl-PCBs, PCDD/Fs + dl-PCBs, and ndl-PCBs were found to be 0.21 pg TEQ/g fat, 0.54 pg TEQ/g fat, 0.75 pg TEQ/g fat, and 0.64 ng/g fat, respectively (Table 1).
No statistically significant differences (p > 0.05) concerning concentrations in all the tested contaminants were found according to the four years of sampling. On the contrary, milk samples from Emilia-Romagna were found to be significantly less contaminated (p ≤ 0.05) by both dioxins and PCBs than those collected from the Lombardy region. Indeed, as can be observed in Table 1, milk samples from Emilia-Romagna presented a 20% lower concentration of PCDD/Fs (0.19 vs. 0.23 pg TEQ/g fat), a 40% lower concentration of dl-PCBs (0.41 vs. 0.66 pg TEQ/g fat), 30% lower concentrations of PCDD/Fs + dl-PCBs (0.60 vs. 0.88 pg TEQ/g fat), and 7% lower concentrations of ndl-PCBs (6.25 vs. 6.71 ng/g fat).
A deeper analysis of data was conducted to investigate whether the increased contamination of milk with PCDD/Fs + dl-PCBs was proportional to the lipid content of milk. A Spearman correlation analysis between the sample contamination levels and the lipid content was performed. Interestingly, the analysis showed a positive correlation between contamination levels and fat content only in milk samples from Emilia-Romagna (r = 0.63, p ≤ 0.05), while no significant correlations were found in samples from Lombardy (r = 0.11, p > 0.05). These findings suggest that Emilia-Romagna may have a relatively consistent and stable level of background contamination to the extent that the fat content of milk from this region could be somewhat indicative of the degree of dioxin contamination. In contrast, milk samples from Lombardy may be affected by background and point contamination, which could explain why some samples with low-fat content showed high contamination levels while some samples with high-fat content showed low contamination levels. This situation may be attributed to the higher level of industrialization in Lombardy, the existence of contaminated sites that produced PCBs in the past, and the proximity of industrial facilities to dairy farms [46].
All milk samples were compliant with the European MLs of PCDD/Fs of 2 pg TEQ/g fat and ALs of 1.75 pg TEQ/g fat, as well as with the European MLs of ndl-PCBs of 40 ng/g fat. Nonetheless, one sample from Lombardy exceeded the MLs of 4 pg TEQ/g fat established for PCDD/Fs + dl-PCBs since contaminated by 6.28 pg TEQ/g fat, while the other four samples from Lombardy were above the ALs of 2 pg TEQ/g fat established for dl-PCBs (Figure 2).
When comparing cow milk contamination levels found in the present work with those reported elsewhere in the world, a quite variable scenario emerged. For example, in eight Brazilian states, cow milk samples were found to have a PCDD/Fs contamination level of 1.36 pg TEQ/g fat, which is approximately seven times higher than the concentration found in this study [10]. The authors attributed the higher contamination levels to the increased proximity of farms to urban centers and the degree of industrialization in the area [10]. In contrast, cow milk samples collected from seven different regions of Chile showed concentrations of PCDD/Fs + dl-PCBs up to 0.45 pg TEQ/g fat, two times lower than those found in this study [47]. The authors confirmed that higher contamination levels were observed in milk samples collected from more populated and industrialized areas [47].
When comparing data to those reported by Asian countries, it was observed that the PCDD/F levels of milk samples from traditional markets around Taiwan tended to be three-four times higher (0.59–0.89 pg TEQ/g fat) than those found in this study [48,49] but no significant differences in terms of dl-PCBs were identified [49]. Significantly lower PCDD/F + dl-PCB concentrations were recorded in milk samples from Malaysia (0.16 pg TEQ/g fat) [50], while concentrations of 0.13 ± 0.10 were reported in Chinese milk samples by Zhang et al. in 2015. However, since that value was expressed on a fresh weight rather than a lipid basis, it should be considered higher than those found in the present work [51].
In Europe, the EFSA reported a mean concentration of PCDD/Fs + dl-PCBs of 1.91 pg TEQ/g fat in milk and dairy products in 2012, which is higher compared to the contamination levels found in this study but includes contributions from milk-derived products (that usually have higher contamination levels than raw milk due to a larger amount of fat) [3]. In the same period, Marin et al. reported an average concentration of PCDD/Fs + dl-PCBs in milk and dairy products from Spain ranging from 0.25 pg TEQ/g fat to 2.13 pg TEQ/g fat [52].
Very high PCDD/F + dl-PCB contamination levels (up to 2–3 pg TEQ/g fat) were also reported in cow milk samples from Italy [7,53], as well as contamination levels up to 5.36 pg TEQ/g fat in buffalo milk from a contaminated site in southern Italy [54]. However, it is important to note that these data were calculated using the old toxic equivalency factors (TEFs) established in 1998 rather than those established in 2005 and currently employed [5]. Therefore, comparisons should be conducted with caution, especially since these surveys were conducted approximately 15 years ago and may not reflect the reduction trend observed in Europe in recent years.
Dioxin and PCB contamination levels in cow milk from Lombardy and Emilia-Romagna collected in this work (2018–2021) and the former monitoring plan (2012–2014) [37] were further compared. As can be observed from Table 2, concentrations of all the analyzed contaminants were significantly lower in 2018–2021 than in the 2012–2014 milk samples, with the highest reduction in dl-PCBs and PCDD/Fs + dl-PCBs actual mass values (pg/g fat or ng/g fat) recorded in milk samples from the Lombardy region (–70% in Lombardy vs. –53% in Emilia-Romagna). Similarly, compared to the 2012–2014 concentrations, the 2018–2021 ndl-PCB values fell to 6.71 ng/g fat (–37%) in Lombardy and to 6.25 ng/g fat (–19%) in Emilia Romagna. The highest PCDD/Fs decrease per year was instead observed in samples from Emilia Romagna (–44%) (Table 2).
The total TEQ values (pg TEQ/g fat) of both PCDD/Fs and dl-PCBs in the 2018–2021 milk samples from the two regions under study declined in direct proportion but less strongly than the actual mass values. Indeed, compared to the 2012–2014 outcomes, Lombardy and Emilia-Romagna saw PCB levels drop by 42% and 37%, respectively (Table 2), suggesting that this drop was mostly due to less toxic PCDD/F and dl-PCB congeners (i.e., congeners having lower TEF values and, hence, a lower impact on total TEQ values).

Variation of Dioxins and PCB Congener Patterns

Since each dioxin and PCB congener bears the imprint of the contamination source (geographical location and industrial activities) and is characterized by a different toxicological effect, monitoring dioxin and PCB congener patterns in food is crucial to understand better their presence in the environment, contamination sources, and the potential risks to human health. Numerous studies have suggested that the chemical species of dioxins and PCBs may vary significantly depending on the foods being considered [3]. Regarding cow milk, dioxin, and PCB chemical species might differ based on a number of variables, including the feeding regimen, age, and breed of the animal, as well as the ambient conditions in which the cattle are grown, the sampling, and the analytical techniques employed for the quantification [55].
To provide complementary information, dioxin and PCB congener profiles in cow milk from the two Italian areas were examined for their variation based on the sample sites and investigated by plotting both the relative percentage contribution of each congener to the total PCDD/Fs + dl-PCBs TEQ values (Figure 3) and the absolute mass fraction concentrations (pg/g fat) (Figure 4).
As seen in Figure 3, PCB 126 was found to be responsible for 66% and 57% of the overall PCDD/Fs + dl-PCBs TEQ values in milk samples from Lombardy and Emilia-Romagna, respectively. These results contrast with those of Barone et al., according to whom cow milk showed a predominance of PCDD/Fs [56], but support earlier findings from our team, who found that PCB 126 was the predominant dl-PCB congener in milk samples in northern Italy [37,57]. Additionally, previous studies from other world regions identified PCB 126 and the other non-ortho PCBs as the congeners having the most impact on the total TEQ values in milk [47,58,59,60,61].
The 2,3,4,7,8-PeCDF was found to be the chemical species with the second largest contribution (13% in samples from Lombardy and 9% in samples from Emilia-Romagna). On the other hand, two mono-ortho species, PCB 189 and PCB 105, showed a greater contribution to total TEQ values in milk samples from Emilia-Romagna (Figure 3).
When analyzing PCDD/F patterns expressed on mass fraction concentrations (pg/g fat), particularly evident differences of 2,3,4,7,8-PeCDF, 1,2,3,4,7,8-HxCDF, 1,2,3,4,6,7,8-HpCDD, and 1,2,3,4,6,7,8,9-OCDD were noted, whose levels were greater in milk from Lombardy (Figure 4A). By plotting contributions to TEQ values, the differences of the latter species were not evident since they were obscured by other congeners dominating the pattern.
PCDDs have typically been linked to agricultural chemicals such as chlorophenols, whereas more mixed PCDD/F profiles have been associated with both waste incineration and the usage of chlorinated agricultural pesticides [55]. Nonetheless, it can be challenging to accurately connect the PCDD and PCDF congener profiles of a specific food or feed to its sources because pollution sources might be numerous and varied.
As highlighted in Figure 3, milk samples from Lombardy also had a more pronounced pattern of two dl-PCB congeners, corresponding to PCB 105 and PCB 118 (Figure 4B). In addition to particular PCB-containing technical formulations used in the past, pollution emissions from the steel sector, cement furnaces, and incineration facilities have been linked to the prevalence of these chemical species in food and feed [62]. Northern Italy hosted chemical industries producing chlorine compounds and actually houses a huge number of industrial facilities potentially responsible for the emission of different environmental pollutants [63]. Therefore, as already reported, the risk of contamination of local feed and food with PCDD/Fs and PCBs may be particularly high [37,57,64].
Differences in ndl-PCB patterns among the two regions were particularly evident for three congeners, namely PCB 138, PCB 153, and PCB 180 (Figure 4C). Regardless of the sampling areas, these chemical species were the most abundant among the six ndl-PCB indicators, possibly because of their greater chlorination degree and longer persistence in the tissues of food-producing animals [55].

3.2. Dietary Exposure to Dioxins and PCBs through Milk Consumption

3.2.1. Screening-level Assessment via Deterministic Calculations

The UB exposure levels to PCDD/Fs + dl-PCBs due to the consumption of whole, semi-skimmed, and skimmed cow’s milk in all Italian age groups are presented in Table 3. The EWIs ranged from 0.02 to 15.23 pg TEQ/kg bw/week, with the lowest values calculated for mean elderly consumers ingesting skimmed milk and the highest values calculated for high infant consumers (95th percentile) ingesting whole milk. Globally, the EWIs of dioxins and PCBs were found to be higher across all population age groups due to the consumption of whole milk, while, regardless of the type of milk consumed, infants, toddlers, and children were found to be the most exposed population groups (Table 3).
A more thorough investigation revealed that for adolescents, adults, and the elderly, the EWIs from whole and semi-skimmed milk consumption were roughly the same (due to a compensation effect between the higher ingestions of semi-skimmed milk and the higher fat content of whole milk). On the other hand, despite the intermediate dietary consumption rates of skimmed milk, the EWIs associated with the consumption of this product were an order of magnitude lower (Table 3).
The youngest segment of the European population was found to be much more likely to be exposed to dioxins and PCBs than the other groups by other authors. For example, Austrian children were reported to ingest 0.77 pg TEQ/kg bw/day of dioxins and PCBs through the whole diet, with milk and dairy contributing to 65% of the overall intake (EWI of 3.50 pg TEQ/kg bw/week) [65]. For French children, EWIs of 0.11–0.63 and 0.21–0.69 (0.77–4.41 and 1.47–4.83 pg TEQ/kg bw/week) were reported in relation to the consumption of cow milk and milk plus dairy products, respectively [66].
As for Italian children, the whole diet was estimated to be responsible for a weekly intake of 1.98–4.98 pg TEQ/kg bw/day, with milk alone contributing to 10% (EWI of 1.38–3.48 pg TEQ/kg bw/week) and, when combined with dairy products, to 34% (EWI of 4.71–11.85 pg TEQ/kg bw/week) of the weekly ingestions of these contaminants [67]. Similarly, in another study, 1.50–2.64 pg/TEQ/kg bw/week of dioxins and PCBs were estimated to be taken by Italian children through the consumption of milk and dairy products [56]. Hence, based on the calculated EWIs and the conventional deterministic risk characterization analysis summarized in Figure 5, a potential health risk concern related to the consumption of whole and semi-skimmed milk emerged, with infants, toddlers, and children being the populations at higher risk. Indeed, average intakes significantly exceeded the toxicological reference values, being 125–320% (95th percentile: 390–762%) and 85–175% (95th percentile: 192–409%) of the TWI for infants and toddlers, respectively (Figure 5). As for children, the contribution of whole and semi-skimmed milk consumption to the TWI was lower (41–82%; 95th percentile: 95–213%) but still concerning, given that many foods, besides milk, provide significant amounts of dioxins and PCBs (Figure 5).
Our risk characterization provides more alarming results compared to those previously reported in Italy, according to which cow milk contributed roughly 20–30% of the TWI [37,53]. Nevertheless, it should be noted that a direct comparison can be misleading since these studies used old dioxin TEF values or compared the exposure level to the Tolerable Daily Intake (TDI) of 2 pg WHO-TEQ/kg bw, which is seven times lower than the actual TWI. Additionally, the conservative nature and simplifications of the deterministic exposure assessment, which means lacking information about the variability and uncertainty in both food consumption statistics and contamination values, typically lead to incomplete or inaccurate (i.e., overestimated) exposure levels, giving information on the possibility but not the probability of risk.

3.2.2. Refined Assessment via Probabilistic Modelling

Monte Carlo simulations were run to estimate in a more accurate and realistic way the distribution of EWIs and the contribution of milk consumption to the TWI and, therefore, to quantify the probability of the health risk occurring. Only those exposure scenarios of particular concern resulting from the deterministic assessment were simulated. Moreover, since the imputation method for the treatment of data below the LOQ might have a large effect on the forecasted outcomes, simulations were run three separate times using UB, MB, or LB concentrations one at a time.
Figure 6 presents the variability in the form of distribution frequencies of EWIs at the MB scenario for the most vulnerable population consuming semi-skimmed and whole milk. Average EWIs were estimated to range from 5.74 ± 7.77 pg/TEQ/kg bw/week for infants consuming whole milk to 0.74 ± 0.97 pg/TEQ/kg bw/week in children consuming semi-skewed milk. The median (P50) MB values were significantly lower when compared to the mean values (Figure 6) because of the right-skewness of the distribution frequency curve, which resulted as such because of the right-skewness of both consumption data and contamination levels used as input values for simulations. This condition led to greater robustness of EWIs at the left tail of the curve (lower percentiles). On the contrary, a greater degree of uncertainty is associated with the interpretation of EWIs beyond the 90th percentile since they are not supported by an adequate number of observations. In addition, from the comparison of EWIs at different 10% increment percentiles (Table 4), it was observed that there was a significant amount of uncertainty due to the method for the imputation of data below the LOQ. Indeed, the MB estimates were roughly 15% higher or lower than the LB or UB estimates, respectively. All these findings suggest a possible overestimation of the higher percentile intakes.
Similar considerations also apply to the outcomes resulting from the Monte Carlo-simulated risk characterization associated with chronic exposure to dioxins and PCBs. The MB 95th percentile contribution percentages of milk ingestion to the TWI ranged, on average, from a minimum of 22% for children ingesting semi-skimmed milk to a maximum of 169% for infants ingesting whole milk (Table 5), being significantly lower than those resulting from the deterministically calculated worst case scenario (Figure 5). Nevertheless, as it can be observed from Table 5, the consumption of whole milk by infants was already above the TWI at the MB 40th percentile, indicating that this population group is at particular risk due to the intake of dioxins and PCBs. At the 80th percentile, only the consumption of semi-skimmed milk by children was below the TWI.
To quantify the uncertainty around when each exposure scenario can occur, the cumulative probability distributions of milk intake percentage contributions to the TWI for infants, toddlers, and children were plotted in Figure 7. Ultimately, the probability for infants to exceed the TWI was 76% for the ingestion of whole milk and 35% for the ingestion of semi-skimmed milk, which means that three out of four and one out of three infants are at a very high risk of chronic toxicity. Moreover, 56% and 26% of toddlers, as well as 22% and 8% of children consuming whole and semi-skimmed milk, respectively, can be exposed to dioxin and PCB levels above the toxicological reference values (Figure 7). Such high percentages were already reported for infants and toddlers through the whole diet [68], but our findings might have more severe implications from a health perspective, given that they are linked to the sole consumption of cow milk, while the diet includes many other sources of these contaminants such as meat, fish, and eggs [3].
Despite these results and the high confidence in probabilistic predictions, there are a few limitations related to the applied methodology that need to be mentioned. In particular, regardless of the substitution method applied, the precision of the measured contaminants may have been biased by the left censorship of these data. Dietary intakes of milk by the general Italian population may not accurately reflect those of the inhabitants of northern Italy due to potential regional variations in dietary habits. Moreover, it was assumed that the entire amount of cow milk consumed originated from local farms in Lombardy and Emilia-Romagna, while commercial products can usually be a mixture of milk from different European (or non-European) countries. The unavailability of data concerning the distribution of the population’s body weight might have been another important limiting factor, which, during probabilistic calculations, has not made it possible to evaluate its impact on the output of the exposure assessment. Finally, the potential variation in the bioavailability of dioxins and PCBs due to milk processing and/or consumer age was not taken into consideration.

4. Conclusions

Despite confirming a significant downward trend in dioxin and PCB concentrations within European food and feed chains, this study raises important questions about the potential health risks associated with dietary long-term exposure to these contaminants for the youngest population. By applying probabilistic techniques, a more accurate and realistic picture of exposure levels to dioxins and PCBs and associated risks were achieved, which allowed estimating that the sole consumption of whole cow milk (albeit broadly compliant with maximum EU limits) may lead to exceeding the tolerable weekly intake in approximately 76, 56, and 22% of infants, toddlers, and children, respectively. This output, which is mainly the consequence of the recent seven-fold reduction in the tolerable weekly intake of dioxins and PCB, clearly worsens the overall risk profile of the young population and underscores the need for a multi-pronged approach involving regulatory measures and public education to reduce dietary exposure further. In this setting, monitoring and surveillance systems represent an essential tool to identify contamination sources and track human exposure. On the other hand, effective communication of both risks and benefits associated with the consumption of milk by vulnerable consumers is necessary. This may include promoting a varied diet and alternating with infant formulae, follow-up formulae, and other baby foods for which stricter limits of dioxins and PCBs are in force.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/foods12091869/s1, Table S1: Limit of quantification (LOQ) and recovery values of the 35 measured PCDD/F and PCB congeners. Table S2: Mean concentrations of the 35 measured PCDD/F and PCB congeners in milk samples from the two Italian regions. Table S3: Mean concentrations of the 35 measured PCDD/F and PCB congeners in milk samples collected during the 4-year sampling plan.

Author Contributions

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

Funding

This research did not receive any specific funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The dataset generated during the current study will be made available upon reasonable request.

Acknowledgments

The authors would like to thank the official veterinarians who carried out the monitoring plan in Lombardy and Emilia Romagna regions.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. IARC Polychlorinated biphenyls and polybrominated biphenyls. In IARC Monograph on the Evaluation of Carcinogenic Risks to Human; International Agency for Research on Cancer: Lyon, France, 2016; Volume 107, pp. 1–501.
  2. IARC Polychlorinated dibenzo-para-dioxins and polychlorinated dibenzofurans. In IARC Monograph on the Evaluation of Carcinogenic Risks to Human; International Agency for Research on Cancer: Lyon, France, 1997; Volume 69, pp. 1–639.
  3. EFSA Update of the monitoring of levels of dioxins and PCBs in food and feed. EFSA J. 2012, 10, 2832. [CrossRef]
  4. WHO. Preventing Disease through Healty Environments. Exposure to Dioxins and Dioxin-like Substances: A Major Public Health Concern; WHO: Geneva, Switzerland, 2019; Available online: https://www.who.int/publications/i/item/WHO-CED-PHE-EPE-19.4.4 (accessed on 11 February 2023).
  5. Van den Berg, M.; Birnbaum, L.S.; Denison, M.; De Vito, M.; Farland, W.; Feeley, M.; Fiedler, H.; Hakansson, H.; Hanberg, A.; Haws, L.; et al. The 2005 World Health Organization reevaluation of human and mammalian toxic equivalency factors for dioxins and dioxin-like compounds. Toxicol. Sci. 2006, 93, 223–241. [Google Scholar] [CrossRef] [PubMed]
  6. Milbrath, M.O.G.; Wenger, Y.; Chang, C.W.C.W.; Emond, C.; Garabrant, D.; Gillespie, B.W.; Jolliet, O. Apparent half-lives of dioxins, furans, and polychlorinated biphenyls as a function of age, body fat, smoking status, and breast-feeding. Environ. Health Perspect. 2009, 117, 417–425. [Google Scholar] [CrossRef] [PubMed]
  7. Esposito, M.; Cavallo, S.; Serpe, F.P.; D’Ambrosio, R.; Gallo, P.; Colarusso, G.; Pellicanò, R.; Baldi, L.; Guarino, A.; Serpe, L. Levels and congener profiles of polychlorinated dibenzo-p-dioxins, polychlorinated dibenzofurans and dioxin-like polychlorinated biphenyls in cow’s milk collected in Campania, Italy. Chemosphere 2009, 77, 1212–1216. [Google Scholar] [CrossRef]
  8. Adekunte, A.O.; Tiwari, B.K.; O’Donnell, C.P. Exposure assessment of dioxins and dioxin-like PCBs in pasteurised bovine milk using probabilistic modelling. Chemosphere 2010, 81, 509–516. [Google Scholar] [CrossRef] [PubMed]
  9. Malisch, R.; Kotz, A. Dioxins and PCBs in feed and food—Review from European perspective. Sci. Total Environ. 2014, 491, 2–10. [Google Scholar] [CrossRef]
  10. Rocha, D.A.M.; Torres, J.P.M.; Reichel, K.; Novotny, E.H.; Estrella, L.F.; Medeiros, R.O.; Netto, A.D.P. Determination of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs) in Brazilian cow milk. Sci. Total Environ. 2016, 572, 177–184. [Google Scholar] [CrossRef]
  11. Weber, R.; Herold, C.; Hollert, H.; Kamphues, J.; Blepp, M.; Ballschmiter, K. Reviewing the relevance of dioxin and PCB sources for food from animal origin and the need for their inventory, control and management. Environ. Sci. Eur. 2018, 30, 1–42. [Google Scholar] [CrossRef]
  12. Schulz, A.J.; Wiesmüller, T.; Appuhn, H.; Stehr, D.; Severin, K.; Landmann, D.; Kamphues, J. Dioxin concentration in milk and tissues of cows and sheep related to feed and soil contamination. J. Anim. Physiol. Anim. Nutr. 2005, 89, 72–78. [Google Scholar] [CrossRef]
  13. Hoogenboom, R.L.A.P.; Van Eijkeren, J.C.H.; Zeilmaker, M.J.; Mengelers, M.J.B.; Herbes, R.; Immerzeel, J.; Traag, W.A. A novel source for dioxins present in recycled fat from gelatin production. Chemosphere 2007, 68, 814–823. [Google Scholar] [CrossRef]
  14. Hoogenboom, R.L.A.P. Animal feed contamination by dioxins, polychlorinated biphenyls (PCBs) and brominated flame retardant. In Animal Feed Contamination—Effects on Livestock and Food Safety; Fink-Gremmels, J., Ed.; Woodhead Publishing Series in Food Science, Technology and Nutrition; Woodhead Publishing Limited: Sawston, UK, 2012; pp. 131–183. [Google Scholar] [CrossRef]
  15. Hoogenboom, R.; Zeilmaker, M.; van Eijkeren, J.; Kan, K.; Mengelers, M.; Luykx, D.; Traag, W. Kaolinic clay derived PCDD/Fs in the feed chain from a sorting process for potatoes. Chemosphere 2010, 78, 99–105. [Google Scholar] [CrossRef] [PubMed]
  16. Eljarrat, E.; Caixach, J.; Rivera, J. Determination of PCDDs and PCDFs in different animal feed ingredients. Chemosphere 2002, 46, 1403–1407. [Google Scholar] [CrossRef] [PubMed]
  17. Adamse, P.; Van der Fels-Klerx, H.J. (Ine.; Schoss, S.; de Jong, J.; Hoogenboom, R.L.A.P. Concentrations of dioxins and dioxin-like PCBs in feed materials in the Netherlands, 2001–2011. Food Addit. Contam. Part A Chem. Anal. Control. Expo. Risk Assess. 2015, 32, 1301–1311. [Google Scholar] [CrossRef] [PubMed]
  18. Pajurek, M.; Warenik-Bany, M.; Mikolajczyk, S. Dioxin transfer simulation from feed to animal tissues and risk assessment. Chemosphere 2023, 313, 137379. [Google Scholar] [CrossRef]
  19. Amutova, F.; Delannoy, M.; Baubekova, A.; Konuspayeva, G.; Jurjanz, S. Transfer of persistent organic pollutants in food of animal origin—Meta-analysis of published data. Chemosphere 2021, 262, 128351. [Google Scholar] [CrossRef]
  20. Lorenzi, V.; Angelone, B.; Ferretti, E.; Galli, A.; Tonoli, M.; Donati, M.; Fusi, F.; Zanardi, G.; Ghidini, S.; Bertocchi, L. PCDD/Fs, DL-PCBs, and NDL-PCBs in Dairy Cows: Carryover in Milk from a Controlled Feeding Study. J. Agric. Food Chem. 2020, 68, 2201–2213. [Google Scholar] [CrossRef]
  21. Kim, M.; Kim, D.G.; Bong, Y.H.; Jang, J.H.; Son, S.W. Concentrations of PCDD/Fs, dioxin-like PCBs, PBDEs, and hexachlorobenzene in fat samples from cattle of different ages and gender in Korea. Food Chem. 2013, 138, 1786–1791. [Google Scholar] [CrossRef]
  22. Girolami, F.; Spalenza, V.; Benedetto, A.; Manzini, L.; Badino, P.; Abete, M.C.; Nebbia, C. Comparative liver accumulation of dioxin-like compounds in sheep and cattle: Possible role of AhR-mediated xenobiotic metabolizing enzymes. Sci. Total Environ. 2016, 571, 1222–1229. [Google Scholar] [CrossRef]
  23. Benedetto, A.; Brizio, P.; Guaraldo, P.; Stella, C.; Cappa, C.; Baioni, E.; Spalenza, V.; Nebbia, C.; Abete, M.C. Dioxins, DL-PCB and NDL-PCB accumulation profiles in livers from sheep and cattle reared in North-western Italy. Chemosphere 2016, 152, 92–98. [Google Scholar] [CrossRef]
  24. Krause, T.; Moenning, J.-L.; Lamp, J.; Maul, R.; Schenkel, H.; Fürst, P.; Pieper, R.; Numata, J. Transfer of polychlorinated dibenzo- p -dioxins and dibenzofurans (PCDD/Fs) and polychlorinated biphenyls (PCBs) from oral exposure into cow’s milk—Part I: State of knowledge and uncertainties. Nutr. Res. Rev. 2022, 4, 1–23. [Google Scholar] [CrossRef]
  25. Lake, I.R.; Foxall, C.D.; Fernandes, A.; Lewis, M.; Rose, M.; White, O.; Dowding, A. Seasonal variations in the levels of PCDD/Fs, PCBs and PBDEs in cows’ milk. Chemosphere 2013, 90, 72–79. [Google Scholar] [CrossRef] [PubMed]
  26. Giannico, O.V.; Fraganelli, G.R.; Baldacci, S.; Desiante, F.; Pellegrino, A.; Basile, F.C.; Franco, E.; Diletti, G.; Conversano, M. Dioxins and PCBs contamination in milk and dairy products from Province of Taranto (Puglia Region, Southern Italy): A six years spatio-temporal monitoring study. Ann. Ist. Super. Sanità 2021, 57, 233–238. [Google Scholar] [CrossRef] [PubMed]
  27. O’Donovan, J.V.; O’Farrell, K.J.; O’Mahony, P.; Buckley, J.F. Temporal trends in dioxin, furan and polychlorinated biphenyl concentrations in bovine milk from farms adjacent to industrial and chemical installations over a 15 year period. Vet. J. 2011, 190, e117–e121. [Google Scholar] [CrossRef] [PubMed]
  28. Italian National Institute of Statistics 2022 Yearly Survey on Milk and Dairy Products. Available online: http://dati.istat.it/Index.aspx?QueryId=25270&lang=en (accessed on 11 February 2023).
  29. Italian National Zootechnical Registry Consistency of Cattle and Buffalo Farms. Available online: https://www.vetinfo.it/j6_statistiche/#/report-pbi/1 (accessed on 11 February 2023).
  30. Hoogenboom, R.L.A.P.; Traag, W.; Fernandes, A.; Rose, M. European developments following incidents with dioxins and PCBs in the food and feed chain. Food Control 2015, 50, 670–683. [Google Scholar] [CrossRef]
  31. EFSA Risk for animal and human health related to the presence of dioxins and dioxin-like PCBs in feed and food. EFSA J. 2018, 16, e05333. [CrossRef]
  32. European Commission Commission Regulation (EU) 2022/2002 of 21 October 2022 amending Regulation (EC) No 1881/2006 as regards maximum levels of dioxins and dioxin-like PCBs in certain foodstuffs. Off. J. Eur. Union 2022, L 274/64. Available online: http://data.europa.eu/eli/reg/2022/2002/oj (accessed on 30 March 2023).
  33. Council of the European Union Council Directive 96/23/EC of 29 April 1996 on measures to monitor certain substances and residues thereof in live animals and animal products and repealing Directives 85/358/EEC and 86/469/EEC and Decisions 89/187/EEC and91/664/EEC. Off. J. Eur. Communities 1996, L 125/10. Available online: http://data.europa.eu/eli/dir/1996/23/oj (accessed on 30 March 2023).
  34. Commission of the European Communities Commission Decision 97/747/EC of 27 October 1997 fixing the levels and frequencies of sampling provided for by Council Directive 96/23/EC for the monitoring of certain substances and residues thereof in certain animal products. Off. J. Eur. Communities 1997, L 303/12. Available online: https://eur-lex.europa.eu/eli/dec/1997/747/oj (accessed on 30 March 2023).
  35. European Commission Regulation (EU) 2017/625 on official controls and other official activities performed to ensure the application of food and feed law, rules on animal health and welfare, plant health and plant protection products. Off. J. Eur. Union 2017, L 95/1. Available online: http://data.europa.eu/eli/reg/2017/625/oj (accessed on 30 March 2023).
  36. European Commission Commission Regulation (EU) 2017/644 of 5 April 2017 laying down methods of sampling and analysis for the control of levels of dioxins, dioxin-like PCBs and non-dioxin-like PCBs in certain foodstuffs and repealing Regulation (EU) No 589/2014. Off. J. Eur. Union 2017, L 92/9. Available online: http://data.europa.eu/eli/reg/2017/644/oj (accessed on 30 March 2023).
  37. Lorenzi, V.; Ghidini, S.; Angelone, B.; Ferretti, E.; Menotta, S.; Fedrizzi, G.; Varisco, G.; Foschini, S.; Diegoli, G.; Bertocchi, L. Three years of monitoring of PCDD/F, DL-PCB and NDL-PCB residues in bovine milk from Lombardy and Emilia Romagna regions (Italy): Contamination levels and human exposure assessment. Food Control 2016, 68, 45–54. [Google Scholar] [CrossRef]
  38. United States Environmental Protection Agency Office of Water Engineering and Analysis Division US EPA. Method 1613/B Tetra-through octa-chlorinated dioxins and furans by isotope dilution HRGC/HRMS: Method 1613 Revision B October 1994; U.S. Environmental Protection Agency Office of Water Engineering and Analysis Division: Washington, DC, USA, 1994; pp. 1–89. Available online: https://nepis.epa.gov/Exe/ZyPDF.cgi/20002GR6.PDF?Dockey=20002GR6.PDF (accessed on 11 February 2023).
  39. United States Environmental Protection Agency Office of Water Engineering and Analysis Division US EPA. Method 1668/C Chlorinated Biphenyl Congeners in Water, Soil, Sediment, Biosolids, and Tissue by HRGC/HRMS: Method 1668 Revision C April 2010; U.S. Environmental Protection Agency Office of Water Engineering and Analysis Division: Washington, DC, USA, 2010; pp. 1–118. Available online: https://www.epa.gov/sites/default/files/2015-09/documents/method_1668c_2010.pdf (accessed on 11 February 2023).
  40. Turrini, A.; Sette, S.; Le Donne, C.; Piccinelli, R.; D’Addezio, L.; Mistura, L.; Ferrari, M.; Catasta, G.; Martone, D. Italian national dietary survey on children population from three months up to nine years old—IV SCAI CHILD. EFSA Support. Publ. 2021, 18, 7087E. [Google Scholar] [CrossRef]
  41. Turrini, A.; Le Donne, C.; Piccinelli, R.; D’Addezio, L.; Mistura, L.; Sette, S.; Martone, D.; Comendador Azcarraga, F.J.; Ferrari, M.; Catasta, G. Italian national dietary survey on adult population from 10 up to 74 years old—IV SCAI ADULT. EFSA Support. Publ. 2022, 19, 7559E. [Google Scholar] [CrossRef]
  42. EFSA The EFSA Comprehensive Food Consumption Database. Available online: https://www.efsa.europa.eu/en/data-report/food-consumption-data (accessed on 12 February 2023).
  43. EFSA Guidance on the use of probabilistic methodology for modelling dietary exposure to pesticide residues. EFSA J. 2012, 10, 2839. [CrossRef]
  44. Chang, W.H.; Chen, H.L.; Lee, C.C. Dietary exposure assessment to perchlorate in the Taiwanese population: A risk assessment based on the probabilistic approach. Environ. Pollut. 2020, 267, 115486. [Google Scholar] [CrossRef]
  45. Shin, E.S.; Park, M.K.; Kim, G.; Barghi, M.; Choi, S.D.; Yang, J.; Chang, Y.S. Dietary exposure and potential human health risk of dioxins in South Korea: Application of deterministic and probabilistic methods. Chemosphere 2022, 291, 133018. [Google Scholar] [CrossRef]
  46. Italian Institute for Environmental Protection and Research Contaminated Sites of National Interest (SIN). Available online: https://www.isprambiente.gov.it/en/activities/soil-and-territory/copy_of_contaminated-sites/contaminated-sites-of-national-interest-sin?set_language=en (accessed on 21 April 2023).
  47. Pizarro-Aránguiz, N.; Galbán-Malagón, C.J.; Ruiz-Rudolph, P.; Araya-Jordan, C.; Maddaleno, A.; San Martin, B. Chemosphere Occurrence, variability and human exposure to Polychlorinated Dibenzo-p-dioxins (PCDDs ), Polychlorinated Dibenzofurans (PCDFs) and Dioxin-Like Polychlorinated Biphenyls (DL-PCBs) in dairy products from Chile during the 2011–2013 survey. Chemosphere 2015, 126, 78–87. [Google Scholar] [CrossRef]
  48. Wang, I.C.; Wu, Y.L.; Lin, L.F.; Chang-Chien, G.P. Human dietary exposure to polychlorinated dibenzo-p-dioxins and polychlorinated dibenzofurans in Taiwan. J. Hazard. Mater. 2009, 164, 621–626. [Google Scholar] [CrossRef]
  49. Chang, J.W.; Liao, P.C.; Lee, C.C. Dietary intake of PCDD/Fs and dioxin-like PCBs from fresh foods around Taiwan. J. Food Drug Anal. 2012, 20, 12. [Google Scholar] [CrossRef]
  50. Leong, Y.H.; Chiang, P.N.; Jaafar, H.J.; Gan, C.Y.; Majid, M.I.A. Contamination of food samples from Malaysia with polychlorinated dibenzo-p-dioxins and dibenzofurans and estimation of human intake. Food Addit. Contam. Part A 2014, 31, 711–718. [Google Scholar] [CrossRef]
  51. Zhang, L.; Yin, S.; Wang, X.; Li, J.; Zhao, Y.; Li, X.; Shen, H.; Wu, Y. Assessment of dietary intake of polychlorinated dibenzo-p-dioxins and dibenzofurans and dioxin-like polychlorinated biphenyls from the Chinese Total Diet Study in 2011. Chemosphere 2015, 137, 178–184. [Google Scholar] [CrossRef] [PubMed]
  52. Marin, S.; Villalba, P.; Diaz-Ferrero, J.; Font, G.; Yusà, V. Congener profile, occurrence and estimated dietary intake of dioxins and dioxin-like PCBs in foods marketed in the Region of Valencia (Spain). Chemosphere 2011, 82, 1253–1261. [Google Scholar] [CrossRef] [PubMed]
  53. Fattore, E.; Fanelli, R.; Turrini, A.; Di Domenico, A. Current dietary exposure to polychlorodibenzo-p-dioxins, polychlorodibenzofurans, and dioxin-like polychlorobiphenyls in Italy. Mol. Nutr. Food Res. 2006, 50, 915–921. [Google Scholar] [CrossRef] [PubMed]
  54. Esposito, M.; Serpe, F.P.; Neugebauer, F.; Cavallo, S.; Gallo, P.; Colarusso, G.; Baldi, L.; Iovane, G.; Serpe, L. Contamination levels and congener distribution of PCDDs, PCDFs and dioxin-like PCBs in buffalo’s milk from Caserta province (Italy). Chemosphere 2010, 79, 341–348. [Google Scholar] [CrossRef] [PubMed]
  55. Hoogenboom, R.L.A.P.; Malisch, R.; van Leeuwen, S.P.J.; Vanderperren, H.; Hove, H.; Fernandes, A.; Schächtele, A.; Rose, M. Congener patterns of polychlorinated dibenzo-p-dioxins, dibenzofurans and biphenyls as a useful aid to source identification during a contamination incident in the food chain. Sci. Total Environ. 2020, 746, 141098. [Google Scholar] [CrossRef] [PubMed]
  56. Barone, G.; Storelli, A.; Busco, A.; Mallamaci, R.; Storelli, M.M. Polychlorinated dioxins, furans (PCDD/Fs) and dioxin-like polychlorinated biphenyls (dl-PCBs) in food from Italy: Estimates of dietaryintake and assessment. J. Food Sci. 2021, 86, 4741–4753. [Google Scholar] [CrossRef] [PubMed]
  57. Bertocchi, L.; Ghidini, S.; Fedrizzi, G.; Lorenzi, V. Case-study and risk management of dioxins and PCBs bovine milk contaminations in a high industrialized area in Northern Italy. Environ. Sci. Pollut. Res. 2015, 22, 9775–9785. [Google Scholar] [CrossRef]
  58. Konuspayeva, G.; Faye, B.; De Pauw, E.; Focant, J.F. Levels and trends of PCDD/Fs and PCBs in camel milk (Camelus bactrianus and Camelus dromedarius) from Kazakhstan. Chemosphere 2011, 85, 351–360. [Google Scholar] [CrossRef]
  59. Chirollo, C.; Ceruso, M.; Pepe, T.; Vassallo, A.; Marrone, R.; Severino, L.; Anastasio, A. Levels and congeners distribution of dioxins, furans and dioxin-like PCBs in buffaloes adipose tissues sampled in vivo and milk. CYTA J. Food 2018, 16, 1109–1114. [Google Scholar] [CrossRef]
  60. Ahmadkhaniha, R.; Nodehi, R.N.; Rastkari, N.; Aghamirloo, H.M. Polychlorinated biphenyls (PCBs) residues in commercial pasteurized cows’ milk in Tehran, Iran. J. Environ. Health Sci. Eng. 2017, 15, 1–6. [Google Scholar] [CrossRef]
  61. Turrio-Baldassarri, L.; Alivernini, S.; Carasi, S.; Casella, M.; Fuselli, S.; Iacovella, N.; Iamiceli, A.L.; La Rocca, C.; Scarcella, C.; Battistelli, C.L. PCB, PCDD and PCDF contamination of food of animal origin as the effect of soil pollution and the cause of human exposure in Brescia. Chemosphere 2009, 76, 278–285. [Google Scholar] [CrossRef] [PubMed]
  62. Pajurek, M.; Warenik-Bany, M.; Mikolajczyk, S. Feed as a source of dioxins and PCBs. Chemosphere 2022, 308, 136243. [Google Scholar] [CrossRef] [PubMed]
  63. European Environment Agency Industrial Reporting Database. Available online: https://www.eea.europa.eu/data-and-maps/data/industrial-reporting-under-the-industrial-6 (accessed on 24 February 2023).
  64. Ghidini, S.; Varrà, M.O.; Bertocchi, L.; Fusi, F.; Angelone, B.; Ferretti, E.; Foschini, S.; Giacometti, B.; Fedrizzi, G.; Menotta, S.; et al. The influence of different production systems on dioxin and PCB levels in chicken eggs from Emilia-Romagna and Lombardy regions (Italy) over 2017–2019 and consequent dietary exposure assessment. Food Addit. Contam. Part A Chem. Anal. Control. Expo. Risk Assess. 2022, 39, 130–148. [Google Scholar] [CrossRef] [PubMed]
  65. Rauscher-Gabernig, E.; Mischek, D.; Moche, W.; Prean, M. Dietary intake of dioxins, furans and dioxin-like PCBs in Austria. Food Addit. Contam. Part A Chem. Anal. Control. Expo. Risk Assess. 2013, 30, 1770–1779. [Google Scholar] [CrossRef]
  66. Sirot, V.; Tard, A.; Venisseau, A.; Brosseaud, A.; Marchand, P.; Le Bizec, B.; Leblanc, J.C. Dietary exposure to polychlorinated dibenzo-p-dioxins, polychlorinated dibenzofurans and polychlorinated biphenyls of the French population: Results of the second French Total Diet Study. Chemosphere 2012, 88, 492–500. [Google Scholar] [CrossRef]
  67. Diletti, G.; Scortichini, G.; Abete, M.C.; Binato, G.; Candeloro, L.; Ceci, R.; Chessa, G.; Conte, A.; Di Sandro, A.; Esposito, M.; et al. Intake estimates of dioxins and dioxin-like polychlorobiphenyls in the Italian general population from the 2013-2016 results of official monitoring plans in food. Sci. Total Environ. 2018, 627, 11–19. [Google Scholar] [CrossRef]
  68. Hulin, M.; Sirot, V.; Vasseur, P.; Mahe, A.; Leblanc, J.C.; Jean, J.; Marchand, P.; Venisseau, A.; Le Bizec, B.; Rivière, G. Health risk assessment to dioxins, furans and PCBs in young children: The first French evaluation. Food Chem. Toxicol. 2020, 139, 111292. [Google Scholar] [CrossRef]
Figure 1. Geographical distribution of bovine farms of the northern Italian regions (Lombardy, red dots; Emilia-Romagna, yellow dots) from which milk samples (N = 214) were collected within 2018–2021.
Figure 1. Geographical distribution of bovine farms of the northern Italian regions (Lombardy, red dots; Emilia-Romagna, yellow dots) from which milk samples (N = 214) were collected within 2018–2021.
Foods 12 01869 g001
Figure 2. Frequency distribution of milk samples (N = 214) by contamination levels of PCDD/Fs (pg TEQ/g fat), dl-PCBs (pg TEQ/g fat), PCDD/Fs + dl-PCBs (pg TEQ/g fat), and ndl-PCBs (ng/g fat). Vertical dashed yellow and red lines mark the exceedance of action levels and maximum levels according to European Union requirements.
Figure 2. Frequency distribution of milk samples (N = 214) by contamination levels of PCDD/Fs (pg TEQ/g fat), dl-PCBs (pg TEQ/g fat), PCDD/Fs + dl-PCBs (pg TEQ/g fat), and ndl-PCBs (ng/g fat). Vertical dashed yellow and red lines mark the exceedance of action levels and maximum levels according to European Union requirements.
Foods 12 01869 g002
Figure 3. Relative percentage contribution of the 29 PCDD/F and dl-PCB congeners to the total PCDD/F + dl-PCB TEQ values (pg TEQ/g fat). Lombardy in red, Emilia Romagna in yellow.
Figure 3. Relative percentage contribution of the 29 PCDD/F and dl-PCB congeners to the total PCDD/F + dl-PCB TEQ values (pg TEQ/g fat). Lombardy in red, Emilia Romagna in yellow.
Foods 12 01869 g003
Figure 4. Radar chart showing congener profiles of PCDD/Fs (pg/g fat, (A)), dl-PCBs (pg/g fat, (B)), and ndl-PCBs (ng/g fat, (C)) in milk samples. Lombardy in red, Emilia Romagna in yellow.
Figure 4. Radar chart showing congener profiles of PCDD/Fs (pg/g fat, (A)), dl-PCBs (pg/g fat, (B)), and ndl-PCBs (ng/g fat, (C)) in milk samples. Lombardy in red, Emilia Romagna in yellow.
Foods 12 01869 g004
Figure 5. Mean and 95th percentile contribution percentages of whole, semi-skimmed, and skimmed milk consumption by different Italian population age groups to the TWI of PCDD/Fs + dl-PCBs calculated via deterministic methods.
Figure 5. Mean and 95th percentile contribution percentages of whole, semi-skimmed, and skimmed milk consumption by different Italian population age groups to the TWI of PCDD/Fs + dl-PCBs calculated via deterministic methods.
Foods 12 01869 g005
Figure 6. Probability distribution functions calculated via probabilistic Monte Carlo simulations of EWIs of PCDD/Fs + dl-PCBs for infants, toddlers, and children consuming semi-skimmed and whole milk.
Figure 6. Probability distribution functions calculated via probabilistic Monte Carlo simulations of EWIs of PCDD/Fs + dl-PCBs for infants, toddlers, and children consuming semi-skimmed and whole milk.
Foods 12 01869 g006
Figure 7. Cumulative probability distribution chart of risk (contribution percentage to the TWI) for infants, toddlers, and children consuming semi-skimmed and whole milk resulting from probabilistic Monte Carlo simulations.
Figure 7. Cumulative probability distribution chart of risk (contribution percentage to the TWI) for infants, toddlers, and children consuming semi-skimmed and whole milk resulting from probabilistic Monte Carlo simulations.
Foods 12 01869 g007
Table 1. Summary of PCDD/F, dl-PCB, and ndl-PCB mean concentrations in cow milk samples collected in Lombardy and Emilia-Romagna regions (Italy) over the 2018–2021 sampling years. Concentrations of PCDD/Fs, dl-PCBs, and PCDD/Fs + dl-PCBs are reported both in pg/g fat and in pg TEQ/g fat (in square brackets). Concentrations of ndl-PCBs are reported in ng/g fat.
Table 1. Summary of PCDD/F, dl-PCB, and ndl-PCB mean concentrations in cow milk samples collected in Lombardy and Emilia-Romagna regions (Italy) over the 2018–2021 sampling years. Concentrations of PCDD/Fs, dl-PCBs, and PCDD/Fs + dl-PCBs are reported both in pg/g fat and in pg TEQ/g fat (in square brackets). Concentrations of ndl-PCBs are reported in ng/g fat.
∑PCDD/Fs∑dl-PCBs∑PCDD/Fs + dl-PCBs∑ndl-PCBs
MeanSDMedianMeanSDMedianMeanSDMedianMeanSDMedian
Region
Lombardy2.07 [0.23]0.5 [0.06]1.92 [0.21]764 [0.66]1238 [0.61]546 [0.53]766 [0.88]1238 [0.63]547 [0.75]6.712.656.00
Emilia-Romagna1.77 [0.19]0.39 [0.04]1.63 [0.18]740 [0.41]534.54 [0.28]597.94 [0.34]741.76 [0.60]534.55 [0.30]599.71 [0.52]6.251.076.00
Year
20181.9 [0.22]0.39 [0.06]1.78 [0.20]1015 [0.59]1258 [0.47]706 [0.46]1017 [0.81]1258 [0.50]707 [0.69]6.661.706.00
20191.94 [0.21]0.47 [0.05]1.82 [0.20]1079 [0.63]995 [0.74]837 [0.49]1081 [0.82]995 [0.77]838 [0.70]6.793.606.00
20202.07 [0.23]0.48 [0.07]1.95 [0.20]883 [0.56]489 [0.39]700 [0.44]885 [0.79]490 [0.42]703 [0.69]6.280.826.00
20212.04 [0.20]0.44 [0.03]1.96 [0.19]783 [0.43]437 [0.25]628 [0.35]785 [0.63]437 [0.27]630 [0.55]6.120.326.00
Table 2. Comparison of PCDD/F, dl-PCB, and ndl-PCB mean concentrations in cow milk samples collected in Lombardy and Emilia-Romagna regions (Italy) over the 2018–2021 sampling years with the previous 2012–2014 survey [37]. Concentrations of PCDD/Fs, dl-PCBs, and PCDD/Fs + dl-PCBs are reported both in pg/g fat and in pg TEQ/g fat (in square brackets). Concentrations of ndl-PCBs are reported in ng/g fat.
Table 2. Comparison of PCDD/F, dl-PCB, and ndl-PCB mean concentrations in cow milk samples collected in Lombardy and Emilia-Romagna regions (Italy) over the 2018–2021 sampling years with the previous 2012–2014 survey [37]. Concentrations of PCDD/Fs, dl-PCBs, and PCDD/Fs + dl-PCBs are reported both in pg/g fat and in pg TEQ/g fat (in square brackets). Concentrations of ndl-PCBs are reported in ng/g fat.
Sampling Years2012–20142018–2021Variation (%)
RegionLombardyEmilia-RomagnaLombardyEmilia-RomagnaLombardyEmilia-Romagna
∑PCDD/Fs3.21 [0.35]3.31 [0.34]2.07 [0.23]1.77 [0.19]−36 [−34]−47 [−44]
∑DL-PCBs2513 [1.14]1591 [0.65]764 [0.66]740 [0.41]−70 [−42]−53 [−37]
∑PCDD/Fs + dl-PCBs2517 [1.49]1595 [0.98]766 [0.88]742 [0.60]−70 [−41]−53 [−39]
∑6 ndl-PCBs10.597.736.716.25−37−19
Table 3. Estimated weekly intake (EWI, pg TEQ/kg bw/week) to PCDD/Fs + dl-PCBs of different Italian population age groups due to the chronic mean and high (95th percentile, in square brackets) consumption of whole (3.5 g fat/100 g), semi-skimmed (1.8 g fat/100 g) and skimmed (0.3 g fat/100 g) cow milk containing average concentrations of 0.78 pg TEQ/g fat of PCDD/Fs + dl-PCBs.
Table 3. Estimated weekly intake (EWI, pg TEQ/kg bw/week) to PCDD/Fs + dl-PCBs of different Italian population age groups due to the chronic mean and high (95th percentile, in square brackets) consumption of whole (3.5 g fat/100 g), semi-skimmed (1.8 g fat/100 g) and skimmed (0.3 g fat/100 g) cow milk containing average concentrations of 0.78 pg TEQ/g fat of PCDD/Fs + dl-PCBs.
Population GroupType of MilkWeekly ConsumptionEWI
(g/kg bw/week)(pg TEQ/kg bw week)
InfantsSemi-skimmed178.5 [555.87]2.51 [7.80]
Whole234.15 [557.89]6.39 [15.23]
ToddlersSemi-skimmed120.68 [272.86]1.69 [3.83]
Whole128.1 [299.32]3.50 [8.17]
ChildrenSkimmed53.76 [141.19]0.06 [0.09]
Semi-skimmed58.08 [135.1]0.83 [1.90]
Whole59.08 [156.1]1.61 [4.26]
AdolescentsSkimmed20.65 [58.66]0.05 [0.14]
Semi-skimmed22.75 [47.39]0.32 [0.67]
Whole17.08 [45.71]0.47 [1.25]
AdultsSkimmed11.13 [25.41]0.03 [0.06]
Semi-skimmed14.42 [29.33]0.20 [0.41]
Whole8.12 [23.45]0.22 [0.64]
ElderlySkimmed 10.50 [25.41]0.02 [0.06]
Semi-skimmed16.87 [30.66]0.24 [0.43]
Whole9.24 [27.23]0.25 [0.74]
Table 4. Percentiles from Monte Carlo simulations of estimated weekly intake (EWI, pg TEQ/kg bw/week) of PCDD/Fs + dl-PCBs through milk consumption.
Table 4. Percentiles from Monte Carlo simulations of estimated weekly intake (EWI, pg TEQ/kg bw/week) of PCDD/Fs + dl-PCBs through milk consumption.
Population GroupSimulation ScenarioMilk TypePercentiles
10th20th30th40th50th60th70th80th90th
Infants LBSemi-skimmed0.170.30.450.640.881.221.732.584.48
Whole0.681.11.562.12.783.654.926.9711.25
MBSemi-skimmed0.220.380.570.791.071.462.032.975.08
Whole0.91.421.972.593.384.385.767.9812.63
UBSemi-skimmed0.280.460.670.921.241.682.313.355.66
Whole1.171.772.43.113.965.056.528.8813.71
ToddlersLBSemi-skimmed0.350.520.690.881.111.391.782.393.61
Whole0.741.081.441.842.332.893.684.927.39
MBSemi-skimmed0.280.420.580.750.961.221.592.163.31
Whole0.580.881.211.561.992.533.284.446.75
UBSemi-skimmed0.360.530.720.881.111.391.782.383.59
Whole0.761.111.461.862.322.93.694.927.3
ChildrenLBSemi-skimmed0.090.150.210.270.360.470.630.881.42
Whole0.270.410.570.740.961.241.632.263.55
MBSemi-skimmed0.120.190.260.350.450.580.751.031.61
Whole0.210.340.470.630.821.081.442.013.24
UBSemi-skimmed0.270.420.570.750.971.241.642.253.51
Whole0.280.460.670.921.241.682.313.355.66
Table 5. Percentiles from Monte Carlo simulations of the risk characterization of PCDD/Fs + dl-PCBs expressed as percentage contribution to the Tolerable weekly intake (% TWI).
Table 5. Percentiles from Monte Carlo simulations of the risk characterization of PCDD/Fs + dl-PCBs expressed as percentage contribution to the Tolerable weekly intake (% TWI).
Population GroupSimulation ScenarioMilk TypePercentiles
10th20th30th40th50th60th70th80th90th
InfantsLBSemi-skimmed8.5415.0522.6231.8144.1160.9186.47128.79224.04
Whole33.8354.9877.88104.98138.96182.39246.23348.6562.56
MBSemi-skimmed11.2319.1628.3639.4253.6873.14101.52148.56254.19
Whole45.1571.1698.31129.71169218.78287.98398.29631.44
UBSemi-skimmed13.7523.1533.5746.2461.9783.81115.72167.39283.14
Whole58.7188.56119.85155.49197.89252.29326.13443.77685.62
ToddlersLBSemi-skimmed17.7125.9634.4343.9355.2669.4289.06119.5180.31
Whole37.0954.2572.0291.85115144.59184.04245.99369.72
MBSemi-skimmed13.9521.2128.8637.4647.8961.1879.51108.23165.22
Whole28.9444.2459.9177.9799.42126.37164.16221.88337.53
UBSemi-skimmed17.9226.3634.8144.255.4669.7189.12119.23179.64
Whole38.0955.6773.0392.82116.04144.94184.53244.97364.84
ChildrenLBSemi-skimmed4.617.3610.2813.7418.0323.6231.5944.1470.96
Whole13.4220.6328.2537.1647.9961.9181.47112.98177.36
MBSemi-skimmed6.089.5913.1417.2622.2828.7937.5751.7480.64
Whole10.5916.9123.7331.4941.1653.9971.85100.38160.07
UBSemi-skimmed7.9611.8915.9120.4425.8832.8642.4057.3887.00
Whole13.7321.0928.837.6548.5362.5782.11112.62176.55
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

Varrà, M.O.; Lorenzi, V.; Zanardi, E.; Menotta, S.; Fedrizzi, G.; Angelone, B.; Gasparini, M.; Fusi, F.; Foschini, S.; Padovani, A.; et al. Safety Evaluation and Probabilistic Health Risk Assessment of Cow Milk Produced in Northern Italy According to Dioxins and PCBs Contamination Levels. Foods 2023, 12, 1869. https://doi.org/10.3390/foods12091869

AMA Style

Varrà MO, Lorenzi V, Zanardi E, Menotta S, Fedrizzi G, Angelone B, Gasparini M, Fusi F, Foschini S, Padovani A, et al. Safety Evaluation and Probabilistic Health Risk Assessment of Cow Milk Produced in Northern Italy According to Dioxins and PCBs Contamination Levels. Foods. 2023; 12(9):1869. https://doi.org/10.3390/foods12091869

Chicago/Turabian Style

Varrà, Maria Olga, Valentina Lorenzi, Emanuela Zanardi, Simonetta Menotta, Giorgio Fedrizzi, Barbara Angelone, Mara Gasparini, Francesca Fusi, Stefano Foschini, Anna Padovani, and et al. 2023. "Safety Evaluation and Probabilistic Health Risk Assessment of Cow Milk Produced in Northern Italy According to Dioxins and PCBs Contamination Levels" Foods 12, no. 9: 1869. https://doi.org/10.3390/foods12091869

APA Style

Varrà, M. O., Lorenzi, V., Zanardi, E., Menotta, S., Fedrizzi, G., Angelone, B., Gasparini, M., Fusi, F., Foschini, S., Padovani, A., & Ghidini, S. (2023). Safety Evaluation and Probabilistic Health Risk Assessment of Cow Milk Produced in Northern Italy According to Dioxins and PCBs Contamination Levels. Foods, 12(9), 1869. https://doi.org/10.3390/foods12091869

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