Next Article in Journal
Microcystin-Bound Protein Patterns in Different Cultures of Microcystis aeruginosa and Field Samples
Next Article in Special Issue
Effect of Various Compounds Blocking the Colony Pigmentation on the Aflatoxin B1 Production by Aspergillus flavus
Previous Article in Journal
The Possible Mechanisms Involved in Degradation of Patulin by Pichia caribbica
Previous Article in Special Issue
Aflatoxin B1 and M1 Degradation by Lac2 from Pleurotus pulmonarius and Redox Mediators
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Effects of Milk Yield, Feed Composition, and Feed Contamination with Aflatoxin B1 on the Aflatoxin M1 Concentration in Dairy Cows’ Milk Investigated Using Monte Carlo Simulation Modelling

RIKILT Wageningen University & Research, P.O. Box 230, Wageningen 6700 AE, The Netherlands
*
Author to whom correspondence should be addressed.
Toxins 2016, 8(10), 290; https://doi.org/10.3390/toxins8100290
Submission received: 5 September 2016 / Revised: 26 September 2016 / Accepted: 30 September 2016 / Published: 9 October 2016

Abstract

:
This study investigated the presence of aflatoxin M1 (AfM1) in dairy cows’ milk, given predefined scenarios for milk production, compound feed (CF) contamination with aflatoxin B1 (AfB1), and inclusion rates of ingredients, using Monte Carlo simulation modelling. The model simulated a typical dairy farm in the Netherlands. Six different scenarios were considered, based on two lactation and three CF composition scenarios. AfB1 contamination of the CF was based on results from the Dutch national monitoring programme for AfB1 in feed materials from 2000 until 2010. Monitoring data from feed materials used in CF production for dairy cattle in the Netherlands were used. Additionally, AfB1 contamination data from an incident in maize in 2013 were used. In each scenario, five different transfer equations of AfB1 from feed to AfM1 in the milk were used, and 1000 iterations were run for each scenario. The results showed that under these six scenarios, the weekly farm concentration of AfM1 in milk was above the EC threshold in less than 1% of the iterations, with all five transfer equations considered. However, this increased substantially in weeks when concentrations from the contaminated maize batch were included, and up to 28.5% of the iterations exceeded the EC threshold. It was also observed that an increase in the milk production had a minimal effect on the exceedance of the AfM1 threshold due to an apparent dilution effect. Feeding regimes, including the composition of CF and feeding roughages of dairy cows, should be carefully considered based on the potential AfM1 contamination of the farm’s milk.

1. Introduction

Aflatoxin B1 (AfB1) is a genotoxic and carcinogenic mycotoxin that is produced by fungi, in particular Aspergillus flavus and Aspergillus parasiticus. AfB1 can be metabolized to aflatoxin M1 (AfM1) by cows, sheep, and goats. Like AfB1, AfM1 is also considered to be genotoxic and carcinogenic to animals and humans [1]. Elevated concentrations of AfB1 in feed result in elevated levels of AfM1 in milk and milk products. Due to their toxic effects on human and animal health, the presence of AfB1 and AfM1 in foodstuffs is strictly regulated within the European Union (EU). Regulation (EC) No 1881/2006 sets the maximum levels for aflatoxin M1 in dairy products, i.e., 0.05 µg/kg for raw milk, heat-treated milk, and milk for the manufacture of milk-based products. In addition, a maximum level for AfB1 in all feed materials was set at 0.02 mg/kg, as well as for compound feed for cattle, sheep, and goats, with the exception of dairy cattle, dairy sheep, and dairy goats at 0.005 mg/kg (Directive 2002/32/EC, Consolidated version 27 February 2015), all relative to a feed with a moisture content of 12%.
Historically, aflatoxins have been mainly found in products originating from countries with tropical weather conditions favourable for the growth of Aspergillus spp., like India, Brazil, and Colombia. However, during the last decade, serious contaminations of maize with aflatoxins have been reported in southern Europe. These include maize grown in 2003 in Italy [2] and maize from the 2012 harvest from the Balkan region [3]. Since climate change is increasingly affecting the formation of mycotoxins in Europe [4], these recent incidents are most likely not the last cases of aflatoxins in maize for European farmers.
Maize is a commonly used ingredient for feeding dairy cows. Feed manufacturers produce feed by mixing and grinding maize and other feed ingredients. In recent decades, the maize consumption of dairy cows has increased due to the prices of raw materials, compound feed composition, and the increase in the amount of concentrated feeds in diets. For example, in the Netherlands, the inclusion rates of maize in compound feed for dairy cows has increased from approximately 1% to 18% in the period 2010–2013 [5]. Additionally, the milk production of dairy cows is steadily increasing, which will coincide with an increased intake of feed, particularly compound feed, by cows.
The abovementioned developments have possibly increased the total exposure of dairy cows to aflatoxins, which might in turn lead to a higher probability of dairy cows’ milk to be contaminated by AfM1.
The aim of this study was to estimate the AfM1 contamination in dairy cows’ milk, using transfer modelling under multiple scenarios of compound feed composition, feed contamination with AfB1, feed consumption, and milk yield. The inclusion rate of maize in compound feed and the contamination of maize by AfB1 were modelled by Monte Carlo simulations, and hence the AfM1 concentration in dairy cows’ milk with the intake distribution could be investigated. Monte Carlo simulation is a computerized mathematical technique that is often used for quantitative analyses and decision-making. The strength of this simulation technique is that it provides the decision-maker with a range of possible outcomes together with their probabilities.

2. Results

2.1. Transfer of AfB1 in Feed to AfM1 in Milk

The transfer of AfB1 in feed to AfM1 in milk was modelled for a typical Dutch dairy farm using transfer modelling combined with Monte Carlo modelling (1000 iterations) under six different scenarios. Namely, three compound feed (CF) composition scenarios, under two different milk yield scenarios (normal and extreme lactation). Five different transfer equations (Materials & Methods Table 4) obtained from literature were used to model the transfer from AfB1 in feed to AfM1 in milk. AfB1 contamination in feed ingredients was modelled in line with the results from the Dutch monitoring data, except for weeks 25 and 26, during which data from a contaminated maize batch were used [3].
Model output is weekly resolved, and for each week, the percentage of simulations (from the 1000 iterations) which resulted in an exceedance of the EC limit of 0.05 µg/kg for AfM1 in milk was calculated. For each scenario, the maximum of these weekly percentage exceedances are shown in Table 1. The week with the maximum percentage exceedance coincided with the use of highly contaminated maize in compound feed. In Table 1, the results of the week with the highest percentage exceedance rate without using the contaminated maize batch are also shown (in italics). Less than 1% of all the weekly simulations were above the EC limit for AfM1 in milk, when only monitoring data were used.
The transfer equation from Veldman et al. [7] resulted in the highest percentage of simulations above the EC limit, with an exceedance in 28.3% of the weekly simulations. The use of a low-protein compound feed (CF composition Scenario 2) resulted in the highest percentage of simulations above the EC limit for all transfer equations. This is in line with the fact that the low-protein compound feed has a high maize inclusion rate. With most of the transfer models, however, no clear differences could be observed in the amount of simulations above the EC threshold between the two lactation scenarios. In the extreme lactation scenario, all cows start lactating at the same time, as opposed to having different cows starting their lactation cycle on different weeks (normal lactation). Under the extreme lactation scenario, even when the cows consume highly contaminated feed in the same weeks (weeks 25 & 26) during their lactation peak, a higher transfer rate of AfM1 has not resulted in an increased exceedance rate due to the high volume of milk produced in the farm, and even a decrease when using the model of Van Eijkeren et al. [9]. In our model, the farm weekly milk production in week 25 amounts to 12,000 kg under the normal lactation milk yield scenario. The extreme lactation scenario resulted in 25% higher milk yield in the same week. Hence, the reason for the similar number of simulations above the EC threshold for both lactation scenarios is probably due to dilution. This dilution effect is most clear for the scenarios using the transfer equation provided in the EFSA opinion [10] from Pettersson [11], where the maximum weekly percentage of simulations above the EC limit is lower under the extreme lactation scenarios, than under the normal lactation scenarios. In the equation from Pettersson [11], the concentration of AfM1 in milk is only dependent on the total AfB1 intake (Materials & Methods Table 4). In the remaining scenarios, the transfer is dependent both on the total AfB1 intake and on the milk yield, and hence this dilution effect is less obvious due to the interaction between the intake and milk yield. Consequently, when comparing the percent of simulations above the EC threshold of a scenario under normal lactation and the corresponding scenario under extreme lactation, the difference is minimal. This implies that a higher milk production has an overall minimal effect on the concentration of AfM1 in the milk of the farm. When using the model from Van Eijkeren et al. [9], the dilution effect is such that with the extreme lactation the probability of AfM1 exceedance is even lower than with the normal lactation.
For each weekly model output, the mean concentration of AfM1 in the milk produced at the farm over 1000 iterations was also calculated. The maximum values for all these weekly mean concentrations are presented in Table 2. For all the scenarios modelled, the maximum weekly mean falls within the EC limit for AfM1 in milk. In some of the considered scenarios, this maximum of the weekly mean concentrations can reach up to 0.04 μg/kg (Table 2). Concentrations of AfM1 in milk as high as 0.32 µg/kg were modelled, which is 6.4 times the EC limit. Notably, this high concentration was calculated under CF Scenario 3, using the transfer equation from Veldman et al. [7]; it also coincided with the weeks when contaminated maize was used.

2.2. Effect of Milk Yield and Feed Intake

Additional scenarios other than the standard scenarios were modelled in order to investigate the interaction between high/low yielding cows and high/low feed intake. The yearly milk yield per cow as modelled under the standard scenarios (Section 2.1) was equivalent to 9111 kg, with high and low yielding cows producing 11,845 and 6378 kg per year, respectively. Similar to the standard scenarios, each week, the percentage of simulations above the EC limit of 0.05 µg/kg for AfM1 in milk was calculated. The maximum values of these weekly percentages for each additional scenario are shown in Table 3. In comparison with Table 2 (specifically CF composition 1 under normal lactation), the high feed (HF) scenarios resulted in a higher weekly percent of simulations exceeding the EC limit. This is true for both the high milk yield (HY) and the low milk yield (LY) cows; however, the effect is larger for the high milk yield scenarios. The low feed (LF) scenarios resulted in a lower percent of simulations exceeding the EC limit, even for simulations with a high milk yield (HY). Additionally, the maximum weekly percent of simulations above the EC limit also coincided with the weeks when contaminated maize was used. If these weeks are excluded, the maximum percent of simulations above the limit is drastically reduced (percent shown in italics).
The mean AfM1 concentration as modelled for the whole farm was calculated every week as the mean of the 1000 iterations. The maximum of these weekly mean AfM1 concentrations in milk followed a similar trend when compared to the standard feed/standard milk yield scenario. The maximum of the weekly mean concentrations in all four scenarios did not exceed the EC limit; however, the highest farm AfM1 concentration modelled was 0.22 µg/kg under the high feed and low yield scenario (HF_LY) when using the equation by Veldman et al. [7] during the weeks with contaminated maize.
The final scenario modelled included corn silage as an additional source of AfB1 in the daily diet. The AfB1 contamination in corn silage was set at 1 µg/kg, being the limit of quantification. When compared to the corresponding scenario without the inclusion of contaminated silage (CF composition 1 under normal lactation), the percentage of simulations above the EC limit for AfM1 in milk is approximately doubled. Nevertheless, excluding weeks 25 and 26, the percent of simulations above the AfM1 limit is still less than 1%. This shows that this additional contamination at a low level does not contribute significantly to raising the daily AfB1 intake above the threshold required to exceed the AfM1 threshold in milk.

2.3. Effect of the Variation in the Transfer Equations

Under the same milk yield scenario, and the same compound composition, the equation used for the transfer of AfB1 in feed to AfM1 in milk played an important role. The exceedance percentage was up to six times higher, depending on which transfer equation is used. The equations established by Masoero et al. [6], Veldman et al. [7], Britzi et al. [8], and Van Eijkeren et al. [9] relate the concentration of AfM1 in milk to the daily milk yield and to the daily intake of AfB1. The equation from Pettersson in 1998 [11]—presented in the EFSA opinion in 2004 [10] on AfB1 in animal feed—related the AfM1 in milk to the daily intake of AfB1, irrespective of the daily milk yield. The transfer rate in the equations from Masoero et al. [6] and Van Eijkeren et al. [9] depends similarly on milk yield, and hence similar results were obtained when using these two transfer equations. On the other hand, the equation from Veldman et al. [7] results in an overall higher transfer rate. Hence, the modelled concentration of AfM1 in milk when using the equation from Veldman et al. [7] will always be the highest. The equation from Britzi et al. [8] has the highest dependence of transfer rate on milk yield. Hence, the modelled concentration of AfM1 in milk will vary greatly depending on the phase in the lactation cycle. On the contrary, in the model described by Van Eijkeren et al. [9] the dependence is very small.

3. Discussion

This study aimed at investigating the potential exceedance of the maximum limit for AfM1 concentrations in milk, under different scenarios for AfB1 contamination of feed ingredients, inclusion rates of ingredients in compound feed, milk production, and transfer rates. Monte Carlo simulation modelling was used to investigate the AfM1 contamination of milk for the wide range of situations that realistically occur in practice. The results show that milk production of the farm and milk yield appeared to have minimal effects on the outcomes; a higher milk production did not result in increased concentrations of AfM1 in the farm milk, due to a potential dilution effect. While, according to transfer equations, a higher milk yield will result in a higher transfer rate, the concentration in the milk produced on the farm changes only slightly due to an increased milk production.
The Monte Carlo approach was applied to simulate the AfB1 contamination in each individual feed ingredient, using monitored mean and standard deviation AfB1 concentrations in the respective feed ingredient monitored over a 10-year period. Thus, the monitoring results were used to fit a distribution for AfB1 in each feed ingredient, and this distribution was in turn used as an input in the Monte Carlo simulation. The same approach was used to simulate the compound feed composition for Scenario 3, provided a known minimum and maximum percentage inclusion for each ingredient. This approach allowed us to simulate 1000 different contamination/composition pairs for the daily feed intake for AfB1, changing every two weeks. With 10-year monitoring data of AfB1 contamination in feed ingredients and the guidelines available for feed composition, less than 0.6% of all the simulations for individual two-week periods were above the EC limit for AfM1 concentration in milk, with a maximum concentration of 0.08 µg/kg (over all five transfer equations). This resulted from either a high inclusion rate of maize with an AfB1 concentration below the EC limit in compound feed, or from a low inclusion rate of maize above the EC limit for AfB1 in the compound feed. Nevertheless, it can be concluded that given the current practices, the probability of exceeding the EC limit of AfM1 in dairy milk is very low. In our model we also tested the effect of corn silage contaminated with AfB1 on the AfM1 content in milk. The contamination of corn silage was set to the limit of quantification of AfB1. With an inclusion rate of corn silage in the dairy cow diet of 27% and a low AfB1 contamination level, a minimal change in the AfM1 content in milk was observed.
Our model concludes that under current inclusion rates of ingredients in compound feed and current AfB1 levels in feed ingredients, the probability of exceeding this EC limit is very low, and only a highly contaminated batch would be cause for concern. When considering the same lactation scenario and the same compound feed composition scenario, an increase in the milk yield and the daily feed intake by a factor of 1.3 resulted in a maximum increase of 0.2% in the probability of exceeding the EC limit in milk. Hence, using the current guidelines for compound feed composition, and the current limits on AfB1 in feed materials, an increase in milk yield does not appear to increase the probability of exceeding the limit of AfM1 in milk.
In 2013, a contaminated shipment of maize intended for feed materials was imported into the Netherlands [3]. The batch had a mean AfB1 concentration of 50.2 (±36.1) µg/kg, which is much higher than the EC legal limit for using maize as an ingredient (being 0.02 mg/kg at a moisture content of 12%). This batch was, however, not found to be contaminated during regular monitoring and was used for the production of a compound feed for dairy cattle. Such a batch was included in our transfer model, also using Monte Carlo simulations for the simulation of the AfB1 contamination of the feed ingredients and for the inclusion rate of each ingredient for a two-week period (weeks 25 and 26). Under normal lactation, 5%–28% of the simulations exceeded the EC limit for AfM1 in milk in the respective weeks. However, the probability, when including this contaminated feed in our model, was considerably dependent on the transfer equation used.
The transfer rate used to set the EC legislative limit in feed for dairy cows was 1%–2%. However, several studies found higher transfer rates of aflatoxins in cows with a higher milk yield [6,7,12] and in early/mid-lactating cows [7,8,13,14,15]. This was incorporated in our model through a 45-week lactation cycle, with milk yield varying through the cycle. The transfer rate varied, possibly due to differences between the metabolisms of the cows, the milk yield, and the source of contamination. In fact, the source of the contamination, milk yield, and cow breed in the studies available varied. Veldman et al. [7] used contaminated groundnut meal, Britzi et al. [8] and Masoero et al. [6] used contaminated corn meal, and the model by Van Eijkeren et al. [9] was fitted to data by Frobish et al. [16] using contaminated cottonseed. Concerning breed, Britzi et al. [8] carried out their study on Israeli Holstein cows (high-yield cows with an average yield of 11,400 kg milk/cow), Masoero et al. [6] carried out their study on Holstein cows, and Veldman et al. [7] did not specify which cows were used. The maximum transfer rate from Veldman et al. [7] and Britzi et al. [8] was about 6%; however, in the model set up by Van Eijkeren et al. [9], the maximum transfer rate was 3.2%. When the model set up by Van Eijkeren et al. [9] was applied to the results from Masoero et al. [6] and Veldman et al. [7], the model did not agree with the data. The rate of exceedance of the EC threshold varied considerably depending on the equation used. However, it is unclear which of these equations is most suitable within our model, and hence all should be considered.

4. Conclusions

Given the ranges in the available transfer equations, only when a highly contaminated batch was included in the model was there a high possibility of exceeding the EC threshold of AfM1 in milk. Additionally, this depended on the inclusion rate of maize in the compound feed. Regarding the increased use of maize in feed for dairy cows, and the increasing milk yield, we hypothesised higher concentrations of AfM1 in dairy cows’ milk. This study showed that under current practise and for the current limits on feed ingredients, an increased milk yield alone will not affect the contamination rate of milk with AfM1. However, to some extent, an increased use of maize in compound feed, combined with a higher contamination of maize with AfB1, will indeed increase the probability of exceedance of the EC limit for AfM1 in milk. Therefore, composition of compound feeds for dairy cows should be carefully performed, and should include information on the potential of AfB1 contamination of the ingredients used, as was shown by Van der Fels-Klerx and Bouzembrak [17], so as to comply with the EC limit for AfB1 in the final feeds.

5. Materials and Methods

5.1. Model

A simulation model was developed in order to estimate the distribution of AfM1 concentrations in dairy milk in a typical farm in the Netherlands. The model was developed in MATLAB R2015b. The model outline, given below, is based on a Monte Carlo simulation (1000 iterations) of the daily intake of aflatoxin B1 from compound feed. Monte Carlo simulation was used to assess the wide variety of possibilities that realistically occur in each of these steps, i.e., in compound feed composition, contamination of the feed ingredients with AfB1, milk production, and transfer of AFB1 to AFM1 in milk.
Compound feed composition scenarios & [AfB1] in compound feed ingredients   i Distribution of daily intake of AfB1 + Different scenarios for milk yield & Different transfer equations Distribution of [AfM1] in milk
In the model, it is assumed that all the cows in the farm are housed indoors throughout the year. The modelled farm was set up based on the survey by Driehuis et al. [18] among 24 dairy farms across the Netherlands. An average herd size of 69 cows was assumed, with an average total daily feed intake (and standard deviation) of 18.7 (1.3) kg DM/cow, of which 4.3 (0.2) kg DM/cow is compound feed.
The milk yield per cow was modelled according to an incomplete gamma model developed by Wood [19] and discussed in Olori et al. [20]. The lactation period per cow was set to 45 weeks, followed by a four-week dry period. The milk yield for the whole farm was modelled in two different milk yield scenarios, namely (1) an individual cow started a new lactation cycle each week (normal lactation scenario); and (2) as an extreme scenario, all cows started the lactation cycle on the same day (extreme lactation scenario).
The concentration of AfM1 in milk ([AfM1]milk) can be modelled through several equations; Table 4 presents the equations published in the scientific literature. Masoero et al. [6],Veldman et al. [7], and Britzi et al. [8] experimentally related the concentration of AfM1 in milk to the total intake of AfB1 (total intakeAfB1) and the daily milk yield (Ymilk). Van Eijkeren et al. [9] estimated the concentration of AfM1 at steady state, also from the total intake of AfB1 and the daily milk yield. The EFSA opinion from 2004 [10] on AfB1 in animal feed uses the equation from Pettersson [11], which is based on a collection of different experimental studies.

5.2. Input Data

5.2.1. AfB1 Contamination in Compound Feed

Data on aflatoxin concentrations in maize and other feed ingredients used in compound feed production for dairy cows, over the years 2000–2010, were extracted from the national monitoring database named KAP (Quality of Agricultural Products) in the Netherlands. These included 3427 records in total. All measurements below the limit of detection (LOD) of 1 µg/kg are recorded as zero in the database. A summary table of the data is presented in Table 5. In addition to the information from the KAP database, data from a recent incident in 2013 [3] of maize from the Balkan region, being contaminated above the maximum level of 20 µg/kg, were also considered. These data were from a shipment of maize dated from March 2013, intended for feed production in the Netherlands, as described by De Rijk et al. [3].

5.2.2. Compound Feed (CF) Composition

Data available for compound feed composition in the Netherlands for lactating dairy cows, including the recipe and the rate to which individual ingredients are used in the compound feed, are presented in Table 6 [21]. These include the mean rate of use of each ingredient in 2014, for each of high- and low-protein compound feed for dairy cattle, as well as the range (minimum–maximum) to which each ingredient is used. Based on these data, three different scenarios for compound feed composition were defined. In the high- and low-protein scenarios, the fixed percentage of each ingredient was used in each simulation, whereas in the third scenario, the minimum and maximum ranges were used to simulate 1000 different compound feed compositions, with composition percentage within the suggested range.
The data from Table 4 and Table 5 were combined to calculate the total intake of AfB1 from compound feed (CF) composed of ingredients i, using the equation below. It is assumed that a compound feed consists of 85% dry matter (DM) [22].
Total daily intake AfB 1 = daily feed intake · %   CF in feed 0.85 · ( n = 1 i [ AfB 1 ] i · %   i   in CF ) .
Additionally, it is assumed that a fresh batch of compound feed is used every two weeks and that all cows eat the same amount of feed, meaning that over a two-week period the total daily intake of AfB1 is constant for all cows. Figure 1 shows an example of the total daily intake for the 1000 iterations, changing every two weeks.

5.3. Scenarios

The Monte Carlo simulation model was run for each of the six different basis scenarios, with two lactation scenarios by three compound feed composition scenarios, using each of the five different transfer equations (Table 4). In each scenario, 1000 iterations were used. AfB1 contamination data from KAP were used for all weeks, except for weeks 25 and 26, when data from the contaminated maize batch was used. Under the extreme lactation scenario, all cows have maximum milk yield during weeks 25 and 26, hence coinciding with the fortnight of contaminated maize.
In order to investigate the interaction between high/low yielding cows and high/low feed intake, additional scenarios were also modelled. High and low yielding cows were modelled as 1.3 and 0.7 times, respectively, of the milk yield used in the scenarios above. In accordance with the average total feed intake (and standard deviation) presented by Driehuis et al. [18], total daily feed intake for high and low feeding cows is modelled as 23 (high feed, HF) and 14.5 kg·DM/cow (low feed, LF), respectively. These four feed/milk yield pairs were modelled with compound feed Scenario 1 (high-protein diet), and the normal lactation scenario. The final scenario modelled, also the worst case scenario, assumed contaminated corn silage at the limit of quantification (LOQ) of AfB1. Information from the KAP database shows that the contamination of corn silage in the Netherlands from 2001 until 2010 is below the LOQ, and hence setting the contamination at the LOQ represents the worst case scenario under current conditions. The percentage of corn silage in the total daily feed intake was set at 27% [18], with a contamination of AfB1 of 1 µg/kg and an assumed dry matter content of 30% [22].

Acknowledgments

The authors acknowledge financial contribution of the Netherlands Ministry of Economic Affairs through the Topsector project AF-14225, and the contribution of the private partners in this project. Leo van Raamsdonk and Ron Hoogenboom, RIKILT Wageningen University and Research, are thanked for their suggestions made to this manuscript.

Author Contributions

H.J.v.d.F.-K. initiated the study; H.J.v.d.F.-K. and L.C. designed the study; L.C. did the programing and running of the models; H.J.v.d.F.-K. and L.C. evaluated the results, and wrote the paper.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. International Agency for Research on Cancer. IARC Monographs on the Evaluation on Carcinogenic Risks to Humans: Aflatoxin M1; International Agency for Research on Cancer: Lyon, France, 1993; Volume 56. [Google Scholar]
  2. Piva, G.; Battilani, P.; Pietri, A. Emerging issues in Southern Europe: Aflatoxins in Italy. In The Mycotoxin Factbook; Barug, D., Bhatnagar, D., van Egmond, H.P., van der Kamp, J.W., van Osenbruggen, W.A., Visconti, A., Eds.; Wageningen Academic Publishers: Wageningen, The Netherlands, 2006; pp. 139–153. [Google Scholar]
  3. De Rijk, T.; van Egmond, H.; van der Fels-Klerx, H.J.; Herbes, R.; de Nijs, M.; Samson, R.; Slate, A.; van der Spiegel, M. A study of the 2013 Western European issue of aflatoxin contamination of maize from the Balkan area. World Mycotoxin J. 2015, 8, 641–651. [Google Scholar] [CrossRef]
  4. Battilani, P.; Toscano, P.; van der Fels-Klerx, H.J.; Moretti, A.; Camardo Leggieri, M.; Brera, C.; Rortais, A.; Goumperis, T.; Robinson, T. Aflatoxin B1 contamination in maize in Europe increases due to climate change. Sci. Rep. 2016, 6. [Google Scholar] [CrossRef] [PubMed]
  5. Remmelink, G.; Wageningen UR Livestock Research, Lelystad, The Netherlands. Use of maize in feeds for cows in The Netherlands in the period 2010–2013. Personal communication, 2013. [Google Scholar]
  6. Masoero, F.; Gallo, A.; Moschini, M.; Piva, G.; Diaz, D. Carry-over of aflatoxin from feed to milk in dairy cows with low or high somatic cell counts. Animal 2007, 1, 1344–1350. [Google Scholar] [CrossRef] [PubMed]
  7. Veldman, A.; Meijs, J.A.C.; Borggreve, G.J.; Heeres-van der Tol, J.J. Carry-over of aflatoxin from cows’ food to milk. Anim. Sci. 1992, 55, 163–168. [Google Scholar] [CrossRef]
  8. Britzi, M.; Friedman, S.; Miron, J.; Solomon, R.; Cuneah, O.; Shimshoni, J.; Soback, S.; Ashkenazi, R.; Armer, S.; Shlosberg, A. Carry-over of aflatoxin B1 to aflatoxin M1 in high yielding Israeli cows in mid- and late-lactation. Toxins 2013, 5, 173–183. [Google Scholar] [CrossRef] [PubMed]
  9. Van Eijkeren, J.C.H.; Bakker, M.I.; Zeilmaker, M.J. A simple steady-state model for carry-over of aflatoxins from feed to cow’s milk. Food Addit. Contam. 2006, 23, 833–838. [Google Scholar] [CrossRef] [PubMed]
  10. European Food Safety Authority. Opinion of the scientific panel on contaminants in the food chain on a request from the Commission related to aflatoxin B1 as undesirable substance in animal feed. EFSA J. 2004, 39, 1–27. [Google Scholar]
  11. Pettersson, H. Concerning Swedish derogation on aflatoxin. Complement to the Memo of 97-03-03 on “Carry-over of aflatoxin from feedingstuffs to milk”. Department of Animal Nutrition and Management, Swedish University of Agricultural Sciences: Uppsala, Sweden, 1998. [Google Scholar]
  12. Japanese Food Safety Commission. Mycotoxin Evaluation Report Aflatoxin M1 in Milk and Aflatoxin B1 in Feed; Food Safety Commission: Tokyo, Japan, 2013.
  13. Diaz, D.; Hagler, W., Jr.; Blackwelder, J.; Eve, J.; Hopkins, B.; Anderson, K.; Jones, F.; Whitlow, L. Aflatoxin Binders II: Reduction of aflatoxin M1 in milk by sequestering agents of cows consuming aflatoxin in feed. Mycopathologia 2004, 157, 233–241. [Google Scholar] [CrossRef] [PubMed]
  14. Sumantri, I.; Murti, T.W.; van der Poel, A.F.B.; Boehm, J.; Agus, A. Carry-over of aflatoxin B1 feed into aflatoxin M1 milk in dairy cows treated with natural sources of aflatoxin and bentonite. J. Indones. Trop. Anim. Agric. 2012, 37, 271. [Google Scholar] [CrossRef]
  15. Kutz, R.E.; Sampson, J.D.; Pompeu, L.B.; Ledoux, D.R.; Spain, J.N.; Vázquez-Añón, M.; Rottinghaus, G.E. Efficacy of Solis, NovasilPlus, and MTB-100 to reduce aflatoxin M1 levels in milk of early to mid lactation dairy cows fed aflatoxin B1. J. Dairy Sci. 2009, 92, 3959–3963. [Google Scholar] [CrossRef] [PubMed]
  16. Frobish, R.A.; Bradley, B.D.; Wagner, D.D.; Long-Bradley, P.E.; Hairston, H. Aflatoxin residues in milk of dairy cows after ingestion of naturally contaminated grain. J. Food Prot. 1986, 49, 781–785. [Google Scholar]
  17. Van der Fels-Klerx, H.J.; Bouzembrak, Y. Modelling approach to limit aflatoxin B1 contamination in dairy cattle compound feed. World Mycotoxin J. 2016, 9, 455–464. [Google Scholar] [CrossRef]
  18. Driehuis, F.; Spanjer, M.C.; Scholten, J.M.; te Giffel, M.C. Occurrence of mycotoxins in feedstuffs of dairy cows and estimation of total dietary intakes. J. Dairy Sci. 2008, 91, 4261–4271. [Google Scholar] [CrossRef] [PubMed]
  19. Wood, P.D.P. Algebraic model of the lactation curve in cattle. Nature 1967, 216, 164–165. [Google Scholar] [CrossRef]
  20. Olori, V.E.; Brotherstone, S.; Hill, W.G.; McGuirk, B.J. Fit of standard models of the lactation curve to weekly records of milk production of cows in a single herd. Livest. Prod. Sci. 1999, 58, 55–63. [Google Scholar] [CrossRef]
  21. Van der Fels-Klerx, H.J.; Adamse, P.; de Nijs, M.; de Jong, J.; Bikker, P. A model for risk-based monitoring of contaminants in feed ingredients. Food Control 2016, in press. [Google Scholar] [CrossRef]
  22. The Professionnal Nutrient Managament Group (The Industry). Tried and Tested: Feed Planning for Cattle and Sheep. Available online: http://www.nutrientmanagement.org/assets/12028 (accessed on 21 September 2015).
Figure 1. Daily intake of AfB1 from dairy cows’ compound feed over for 1000 simulations over the whole lactation period.
Figure 1. Daily intake of AfB1 from dairy cows’ compound feed over for 1000 simulations over the whole lactation period.
Toxins 08 00290 g001
Table 1. Maximum weekly percentage * of simulations above the threshold of AfM1 in milk from the whole farm. Numbers in italics represent the maximum percent of simulations in all weeks excluding weeks 25 and 26 (when contaminated maize was used).
Table 1. Maximum weekly percentage * of simulations above the threshold of AfM1 in milk from the whole farm. Numbers in italics represent the maximum percent of simulations in all weeks excluding weeks 25 and 26 (when contaminated maize was used).
CF Composition ScenarioMilk Yield ScenarioTransfer Model
Masoero et al. [6]Veldman et al. [7]Britzi et al. [8]Van Eijkeren et al. [9]Pettersson from EFSA Opinion [10]
1normal4.90.016.50.37.00.16.00.012.50.2
extreme4.80.016.30.38.90.14.70.112.30.1
2normal11.20.028.30.515.00.113.70.123.00.2
extreme11.90.128.50.317.30.111.20.122.80.2
3normal8.60.320.90.611.60.310.30.316.60.5
extreme7.50.118.90.311.20.27.20.214.80.2
* calculated as (the number of simulations in which the farm milk concentration is above 0.05 µg/kg) /1000 simulations × 100.
Table 2. Maximum of weekly mean AfM1 concentrations (µg/kg) in milk from the whole farm (over all iterations) *.
Table 2. Maximum of weekly mean AfM1 concentrations (µg/kg) in milk from the whole farm (over all iterations) *.
CF Composition ScenarioMilk Yield ScenarioTransfer Model
Masoero et al. [6]Veldman et al. [7]Britzi et al. [8]Van Eijkeren et al. [9]Pettersson [11] from EFSA Opinion [10]
1normal0.0150.0280.0180.0170.029
extreme0.0150.0280.0200.0150.029
2normal0.0220.0400.0260.0250.037
extreme0.0220.0410.0290.0220.037
3normal0.0180.0330.0210.0200.032
extreme0.0170.0310.0220.0170.031
* Including the contaminated batch in weeks 25 & 26. Please note that the maximum of the highest of the weekly mean concentrations is always seen with the contaminated batch.
Table 3. Maximum weekly percentage * of simulations above the AfM1 threshold in milk from the whole farm (HF = High Feed; LY = Low Yield). Numbers in italics represent the maximum percent of simulations in all weeks excluding weeks 25 and 26 (when contaminated maize was used).
Table 3. Maximum weekly percentage * of simulations above the AfM1 threshold in milk from the whole farm (HF = High Feed; LY = Low Yield). Numbers in italics represent the maximum percent of simulations in all weeks excluding weeks 25 and 26 (when contaminated maize was used).
Feed and Yield ScenarioMasoero et al. [6]Veldman et al. [7]Britzi et al. [8]Van Eijkeren et al. [9]Pettersson [11] from EFSA Opinion [10]
HF_HY8.70.123.40.514.50.36.90.016.70.3
HF_LY6.10.122.20.410.00.115.00.317.20.3
LF_HY1.80.011.30.16.10.01.00.07.30.1
LF_LY0.60.010.60.12.70.05.40.08.00.1
contaminated silage8.60.133.00.613.60.312.00.223.80.4
* calculated as (the number of simulations in which the farm milk concentration is above 0.05 µg/kg) /1000 simulations × 100.
Table 4. Equations used for modelling the transfer of AfB1 in feed to AfM1 in dairy milk.
Table 4. Equations used for modelling the transfer of AfB1 in feed to AfM1 in dairy milk.
EquationSource
[ AfM 1 ] milk ( μ g AfM 1 / kg milk ) = Total   intake AfB 1 · Y milk 1 · [ Y milk · 0.13 0.26 ] / 100 [6]
[ AfM 1 ] milk ( μ g AfM 1 / kg milk ) = Total   intake AfB 1 · Y milk 1 · [ Y milk · 0.077 0.326 ] / 100 [7]
[ AfM 1 ] milk ( μ g AfM 1 / kg milk ) = Total   intake AfB 1 · Y milk 1 · [ 0.5154 · e Y milk · 0.0521 ] / 100 [8]
[ AfM 1 ] milk   ( μ g AfM 1 / kg milk ) = Total   intake AfB 1 · [ Y milk · 0.032 ] · ( 17 + Y milk ) 1 [9]
[ AfM 1 ] milk ( ng AfM 1 / kg milk ) = Total   intake AfB 1 · 0.787 + 10.95 [10]
Table 5. Summary data for AfB1 concentration (µg/kg) in individual feed ingredients.
Table 5. Summary data for AfB1 concentration (µg/kg) in individual feed ingredients.
Ingredient# of RecordsMin. [AfB1]Max. [AfB1]Mean [AfB1]Std. Deviation
Wheat 134601.40.0190.148
Barley 115501.50.0160.144
Corn 176801150.6535.52
Triticale 12401.90.0790.388
Rye 1230000
Soybean meal 175105.000.03950.269
Sunflower scrap 113607.51.091.72
Palm kernel 1484026.00.5153.08
Rapeseed scrap 1710000
Corn gluten feed 1517042.00.4612.03
Flour 120000
Citrus pulp 111401.50.0220.168
Dried beet pulp 1340000
Molasses 120000
Contaminated maize 2726.216850.236.1
1 Data from the KAP database, 2 Data from De Rijk et al. [3].
Table 6. Feed composition (%) under a high-protein (Scenario 1) and a low-protein diet (Scenario 2). The general guidelines for feed composition are provided in Scenario 3 as the minimum and maximum percentages for each feed material.
Table 6. Feed composition (%) under a high-protein (Scenario 1) and a low-protein diet (Scenario 2). The general guidelines for feed composition are provided in Scenario 3 as the minimum and maximum percentages for each feed material.
Feed IngredientsHigh-ProteinLow-ProteinMinimumMaximum
CF composition scenario123
Total grains 2060
Wheat56035
Barley02.3900
Corn10.2415.06035
Triticale1.21.9015
Rye01015
Soybean meal14.960.23030
Sunflower seed meal4.53.83025
Palm kernel15.0115020
Rapeseed meal7.945.54030
Corn gluten feed3.671030
Flour00.04020
Dried beet pulp0.087.86040
Citrus pulp03.37025
Molasses1.51.55010

Share and Cite

MDPI and ACS Style

Van der Fels-Klerx, H.J.; Camenzuli, L. Effects of Milk Yield, Feed Composition, and Feed Contamination with Aflatoxin B1 on the Aflatoxin M1 Concentration in Dairy Cows’ Milk Investigated Using Monte Carlo Simulation Modelling. Toxins 2016, 8, 290. https://doi.org/10.3390/toxins8100290

AMA Style

Van der Fels-Klerx HJ, Camenzuli L. Effects of Milk Yield, Feed Composition, and Feed Contamination with Aflatoxin B1 on the Aflatoxin M1 Concentration in Dairy Cows’ Milk Investigated Using Monte Carlo Simulation Modelling. Toxins. 2016; 8(10):290. https://doi.org/10.3390/toxins8100290

Chicago/Turabian Style

Van der Fels-Klerx, H. J., and Louise Camenzuli. 2016. "Effects of Milk Yield, Feed Composition, and Feed Contamination with Aflatoxin B1 on the Aflatoxin M1 Concentration in Dairy Cows’ Milk Investigated Using Monte Carlo Simulation Modelling" Toxins 8, no. 10: 290. https://doi.org/10.3390/toxins8100290

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