Hostname: page-component-78c5997874-ndw9j Total loading time: 0 Render date: 2024-11-05T02:59:13.375Z Has data issue: false hasContentIssue false

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

Published online by Cambridge University Press:  12 September 2022

Torsten Krause
Affiliation:
Department of Safety and Quality of Milk and Fish, Max Rubner-Institut (MRI), Hermann-Weigmann-Straße 1, 24103 Kiel, Germany
Jan-Louis Moenning
Affiliation:
Department Safety in the Food Chain, German Federal Institute for Risk Assessment (BfR), Max-Dohrn-Straße 8-10, 10589 Berlin, Germany
Julika Lamp
Affiliation:
Department of Safety and Quality of Milk and Fish, Max Rubner-Institut (MRI), Hermann-Weigmann-Straße 1, 24103 Kiel, Germany
Ronald Maul
Affiliation:
Department of Safety and Quality of Milk and Fish, Max Rubner-Institut (MRI), Hermann-Weigmann-Straße 1, 24103 Kiel, Germany
Hans Schenkel
Affiliation:
Department of Animal Nutrition, University of Hohenheim, Emil-Wolff-Str. 10, 70599 Stuttgart, Germany
Peter Fürst
Affiliation:
Chemical and Veterinary Analytical Institute Münsterland-Emscher-Lippe (CVUA-MEL), Joseph-König-Straße 40, 48147 Münster, Germany
Robert Pieper
Affiliation:
Department Safety in the Food Chain, German Federal Institute for Risk Assessment (BfR), Max-Dohrn-Straße 8-10, 10589 Berlin, Germany
Jorge Numata*
Affiliation:
Department Safety in the Food Chain, German Federal Institute for Risk Assessment (BfR), Max-Dohrn-Straße 8-10, 10589 Berlin, Germany
*
*Corresponding author: Jorge Numata, email [email protected]
Rights & Permissions [Opens in a new window]

Abstract

Polychlorinated dibenzo-para-dioxins (PCDDs) and dibenzofurans (PCDFs) (collectively and colloquially referred to as ‘dioxins’) as well as polychlorinated biphenyls (PCBs) are persistent and ubiquitous environmental contaminants that may unintentionally enter and accumulate along the food chain. Owing to their chronic toxic effects in humans and bioaccumulative properties, their presence in feed and food requires particular attention. One important exposure pathway for consumers is consumption of milk and dairy products. Their transfer from feed to milk has been studied for the past 50 years to quantify the uptake and elimination kinetics. We extracted transfer parameters (transfer rate, transfer factor, biotransfer factor and elimination half-lives) in a machine-readable format from seventy-six primary and twenty-nine secondary literature items. Kinetic data for some toxicologically relevant dioxin congeners and the elimination half-lives of dioxin-like PCBs are still not available. A well-defined selection of transfer parameters from literature was statistically analysed and shown to display high variability. To understand this variability, we discuss the data with an emphasis on influencing factors, such as experimental conditions, cow performance parameters and metabolic state. While no universal interpretation could be derived, a tendency for increased transfer into milk is apparently connected to an increase in milk yield and milk fat yield as well as during times of body fat mobilisation, for example during the negative energy balance after calving. Over the past decades, milk yield has increased to over 40 kg/d during high lactation, so more research is needed on how this impacts feed to food transfer for PCDD/Fs and PCBs.

Type
Review Article
Creative Commons
Creative Common License - CCCreative Common License - BYCreative Common License - NCCreative Common License - SA
This is an Open Access article, distributed under the terms of the Creative Commons Attribution-NonCommercial-ShareAlike licence (https://creativecommons.org/licenses/by-nc-sa/4.0/), which permits non-commercial re-use, distribution, and reproduction in any medium, provided the same Creative Commons licence is included and the original work is properly cited. The written permission of Cambridge University Press must be obtained for commercial re-use.
Copyright
© The Author(s), 2022. Published by Cambridge University Press on behalf of The Nutrition Society

Introduction

The term ‘dioxins’ is colloquially used for a group of highly toxic organochlorine compounds consisting of seventy-five polychlorinated dibenzo-para-dioxins (PCDDs) and 135 polychlorinated dibenzofurans (PCDFs) differing in number and position of benzylic chlorine atoms. They arise as unintentional industrial by-products of, for example, waste incineration, cement works or metal industry and are formed to a lesser degree during forest fires and volcanic eruptions(Reference Safe1Reference Fiedler3).

Another type of organochlorine compound is represented by the group of polychlorinated biphenyls (PCBs) consisting of 209 congeners, each differing in position and number of chlorine atoms. Because of their useful properties, such as flame retardancy and electrical insulating capacity, they were industrially produced beginning in the 1930s(Reference Breivik, Sweetman and Pacyna4). Popular brands of PCB mixtures in the past included Aroclor, Chlophen and Kanechlor, which were commonly produced until their production was banned. The total worldwide production of PCBs is estimated to have been 1.3–1.4 million metric tonnes(Reference Breivik, Sweetman and Pacyna5). Despite the ban on intentional production, unintentionally formed PCBs were recently found in paint colourants(Reference Jahnke and Hornbuckle6) and polymers such as silicones(Reference Herkert, Jahnke and Hornbuckle7). Furthermore, PCDD/Fs and PCBs are unintentionally formed in many combustion processes(Reference Huang and Buekens8).

Among PCBs, twelve congeners are referred to as dioxin-like PCBs (dl-PCBs) owing to their molecular conformational and toxicological similarity to dioxins. The remaining PCBs exhibit different toxicological properties and are therefore labelled non-dioxin-like PCBs (ndl-PCBs). From an analytical perspective, detailed full-congener ndl-PCB measurements are time consuming and expensive. Therefore, six of these ndl-PCBs (Ballschmiter congener identification numbers 28, 52, 101, 138, 153 and 180(Reference Ballschmiter and Zell9)) are used as indicator ndl-PCBs for risk assessment purposes(Reference Weber, Herold and Hollert10). These six indicator ndl-PCBs represented approximately 50% of the total ndl-PCB content in food samples, and at least one of the six indicator ndl-PCBs was quantified in 68.4% of feed and 82.6% of food samples(11).

Both PCDD/Fs and PCBs are highly lipophilic and extremely stable compounds, with environmental degradation half-life periods up to several decades(Reference Sinkkonen and Paasivirta12). Once released, they persist in the environment, undergo transport and may end up in foods of animal origin mainly via contaminated feed and ingested soil(Reference Hoogenboom, Traag and Fernandes13). PCDD/Fs and PCBs are still considered one of the main contaminants in the food chain(Reference Hoogenboom, Malisch and van Leeuwen14). Owing to their lipophilic properties, organochlorine chemicals undergo bioaccumulation, especially in fatty tissues and other fat-containing foods(Reference MacLachlan and Bhula15). Owing to their ubiquitous occurrence, human exposure to PCDD/Fs and PCBs is inevitable, particularly through foods of animal origin, including ‘milk and milk products’(Reference Knutsen, Alexander and Barregård16). Understanding the dynamics of feed to food transfer(17) (historically also called ‘carry-over’) of these contaminants is thus crucial to perform appropriate health risk assessment.

Chronic and high exposure to PCDD/Fs and dl-PCBs can induce endocrine and reproductive interferences(Reference Bruner-Tran, Duleba and Taylor18), impair immunological reactions and cause developmental disruption as well as cancer(Reference Schecter and Gasiewicz19). Nevertheless, the most sensitive effect (i.e. the effect occurring at the lowest dose in the most sensitive population) is related to a reduction in sperm concentration for people exposed in early years(Reference Mínguez-Alarcón, Sergeyev and Burns20). The effects that PCDD/Fs and dl-PCB have on human health result mainly from the toxic activation of the aryl hydrocarbon receptor (AhR)(Reference Okey21). AhR influences numerous metabolic pathways(Reference Abel and Haarmann-Stemmann22). Of all possible congeners, only seven PCDDs, ten PCDFs and twelve dl-PCBs can bind to the AhR(Reference Birnbaum23), making these compounds very relevant for research and risk assessment. The toxicological potency of these twenty-nine congeners is based on 2,3,7,8-TCDD, also known as Seveso-Dioxin, which is the lead compound with the highest toxicity of this group. According to the congener-specific toxicity relative to 2,3,7,8-TCDD, toxic equivalency factors (TEF) were derived to ease toxicological risk assessment. Internationally, TEFs for PCDD/Fs were defined by NATO/CCMS in 1988 as I-TEF and for PCDD/Fs and dl-PCBs by the World Health Organization (WHO) in 1997(Reference Van den Berg, Birnbaum and Bosveld24,Reference Kutz, Barnes and Bottimore25) . The WHO TEFs were again revised in 2005(Reference Van den Berg, Birnbaum and Denison26). Using the TEF concept, the concentration of the different congeners can be summarised as toxicity-weighted masses or toxic equivalent (TEQ), which is used for risk assessment in feed and food worldwide as well as in management measures (i.e. the setting of maximum levels/content in food and feed). Owing to their toxicological properties, WHO International Agency for Research on Cancer (IARC) has listed 2,3,7,8-TCDD, 2,3,4,7,8-PeCDD and PCB-126 in Group 1 as ‘carcinogenic to humans’(27). Furthermore, the international Stockholm Convention on Persistent Organic Pollutants (2001) lists intentionally and unintentionally produced PCDD/Fs and PCBs as part of the so-called dirty dozen, a compilation of persistent organic pollutants known to cause adverse effects in humans and the environment(Reference Lallas28).

The toxicity of non-dioxin-like PCBs (ndl-PCBs) cannot be summarised through the interaction with a single receptor. Adverse effects on the thyroid, liver and brain biochemistry, immunotoxicity, oestrogenicity, reproductive and neurodevelopmental effects were reported in laboratory animals after exposure to ndl-PCBs(29). Since ndl-PCBs are frequently accompanied by PCDD/Fs and dl-PCBs, as well as other toxic compounds, establishing an unequivocal congener-specific toxicological profile is not straightforward(30).

In addition to the inevitable environmental background levels, recurrent contamination incidents with PCDD/Fs and PCBs have affected the safety of foods of animal origin during the past decades(Reference Malisch31). The first reported food chain incident occurred in the 1950s in the United States, where thousands of chickens died of ‘chicken-oedema-disease’, caused by diets containing PCDD/F-contaminated fat(Reference Sanger, Scott and Hamdy32,Reference Higginbotham, Huang and Firestone33) . Also in recent times, feed and food incidents have occurred frequently with different sources of the PCDD/Fs and PCBs. For example: contaminated fat intended as feed material in Belgium(Reference Bernard, Hermans and Broeckaert34), contaminated citrus pulp from Brazil causing elevated PCDD/F levels in European milk products(Reference Malisch35), bakery waste used as animal feed in Germany contaminated during the drying process(Reference Hoogenboom, Bovee and Portier36), German beet pulp contaminated during the drying process leading to elevated milk levels in the Netherlands(Reference Hoogenboom, Traag and Fernandes13), choline chloride premix in animal feed contaminated with PCDD/Fs and PCBs by pine saw dust(Reference Llerena, Abad and Caixach37), flooding of the river Elbe and subsequent contamination of flooded grazing areas in 2003(Reference Stachel, Christoph and Gotz38), recycled feeding fat from gelatine production contaminated by a broken filter(Reference Hoogenboom, Eijkeren and Zeilmaker39) and the Irish dioxin case in 2008 caused by PCB containing fuel used to dry animal feed, which led to the culling of 5707 cattle and 170 605 pigs(Reference Tlustos, Sheridan and O’Sullivan40). As suggested above, many of these feed and food incidents affected cattle(Reference Brambilla, Fochi and Falce41Reference Piskorska-Pliszczynska, Maszewski and Mikolajczyk43), leading to elevated PCDD/F and dl-PCB levels in the milk. This is particularly concerning, as a bulk of the PCDD/F and PCB contamination in humans stems from the consumption of milk and milk products(Reference Knutsen, Alexander and Barregård16,Reference Fürst, Beck and Theelen44) . Each source of PCDD/Fs and PCBs is associated with characteristic congener patterns that can be used to forensically relate contamination incidents and infer their possible sources(Reference Hoogenboom, Malisch and van Leeuwen14). However, it should be noted that passage through the animal alters the original source congener pattern in a way that can be predicted using congener-specific toxicokinetic models.

To protect the consumers, organisations such as the European Commission have defined action and maximum levels for the PCDD/F and dl-PCB content as well as maximum levels for ndl-PCB content in different animal feed and foods(4547). Following the risk-assessment-derived tolerable weekly intake of 2 pg WHO-PCDD/F-PCB-TEQ/(kg body weight × week)(Reference Knutsen, Alexander and Barregård16), maximum levels might be revised (i.e. lowered) in the near future. Although emissions of PCDD/Fs have declined by 45% between 1990 and 2012 in the European Union(48), our understanding of their toxicology has led to lower TWIs, meaning that they remain an important group of contaminants for research and regulation. Exceedance of maximum levels leads to non-marketable food and feed products. In the case of contamination of cattle feed and subsequent contamination of milk, the length of the depuration period depends on the initial intake and on the daily output via milk. In a case described by Malisch et al.(Reference Malisch35), this depuration period may take 1 year before reaching levels below maximum levels. Consequently, the resulting milk may not be marketable for a long period of time, resulting in considerable economic damages(Reference Lascano Alcoser, Velthuis and Hoogenboom49). To improve the risk management of recurrent PCDD/F and PCB incidents, reliable data and models on the transfer behaviour of all toxicologically relevant congeners are needed. In the case of contaminated animal feed, knowledge about the transfer of these congeners from feed to milk is necessary to estimate the extent of milk contamination. In the case of a contamination event, information about the elimination kinetics including the half-life of the congeners is indispensable to predict the length of the depuration period needed for the cow’s milk concentrations to fall below legally binding maximum levels.

Since the 1970s, many studies have analysed the transfer of organochlorines, specifically PCDD/Fs and PCBs from feed to milk, which have been partly reviewed elsewhere(Reference Akkan, Oltmanns and Kalberlah50Reference Amutova, Delannoy and Baubekova55). Since the early 1990s, attempts have been made to develop toxicokinetic models to quantitatively predict the transfer of different congeners into milk, which is discussed in detail in part II of this review(Reference Moenning, Krause, Lamp, Maul, Schenkel, Fürst, Pieper and Numata56). However, no in-depth comparison of data and knowledge in terms of their applicability to risk analysis has yet been made. In this review, seventy-six studies and secondary literature on this topic were evaluated with a focus on the quality and usability of the results for risk assessment. The main goals of the review are:

  • evaluate the availability of transfer parameters (such as transfer rates and elimination half-lives) for all toxicologically relevant congeners (seven PCDDs, ten PCDFs, twelve dl-PCBs) as well as the indicator ndl-PCBs in terms of their applicability for risk assessment;

  • identify future research focus areas regarding missing data with respect to the factors that influence transfer, for example metabolic status and breed.

Variability of experiments and reports on the transfer of PCDD/Fs and PCBs into milk

For this review, seventy-six experimental and theoretical studies from the past 50 years, as well as twenty-nine secondary literature items on the transfer of PCDD/Fs and PCBs into milk of ruminants, were consulted. The complete dataset is compiled in machine-readable form as Supplementary Materials Part B. In these studies, various experimental approaches were followed to derive kinetic parameters to describe the feed to milk transfer. We present an overview of the variability between the studies in terms of exposure scenarios, animal breeds (Table 1) and data integrity and how these differences influence the comparability and our choice of studies.

Table 1. Cattle breeds identified in the various transfer studies

Exposure scenarios in dairy cattle

Regarding the basic experimental design, the studies can be classified into field, incident, controlled feeding and mass balance studies. In field studies, in a stable or on pasture, the daily intake of contaminants could be assumed as being steady, but especially for grazing animals, this assumption is uncertain. Incidents and case reports refer to short-term contamination events where the origin and extent of the exposure may be uncertain.

In controlled feeding studies, the daily amount of contaminant uptake is known; contaminants are either spiked into the daily ration or they are already present in known amounts, with the total feed intake recorded. Alternatively, they are applied as a bolus (e.g. gelatine capsule) into the rumen. The time of exposure in feeding studies varied from 1 d(Reference Slob, Olling and Derks62,Reference Gardner, Righter and Roach75,Reference Jones, Safe and Morcom79,Reference Jones, Safe and Morcom80,Reference Heeschen and Blüthgen87,Reference Ewers88,Reference Platonow, Funnell and Bullock95,Reference Platonow and Meads97,Reference Arnott, Bullock and Platonow98,Reference Olling, Derks and Berende102,Reference Olling, Berende and Derks103,Reference Safe, Platonow and Hutzinger120,Reference Sips, van Eijkeren and Bode124Reference Ruoff and Fouzy126) to several weeks (<28 d)(Reference Steimer, Blüthgen and Heeschen60,Reference Jensen and Hummel77,Reference Willett and Liu78,Reference Arnott, Bullock and Platonow98,Reference McLachlan and Richter113) (Table 2). In some feeding studies, the animals were exposed to PCDD/Fs and PCBs during their dry period(Reference Roos, Berende and Traag61,Reference Ewers88,Reference Ruoff91,Reference Schulz93,Reference Schulz, Wiesmüller and Appuhn94,Reference Klein, Drochner and Forschner100,Reference Driesen, Lerch and Siegenthaler114,Reference Perry, Everson and Hendrix115,Reference Perry, Everson and Hendrix119) . Since elimination via lactation is not possible in this case, PCDD/Fs and PCBs accumulate in the gravid cows’ body fat – and likely in the unborn foetus(Reference Hirako127). Mass balance studies are a subcategory of feeding studies in which feed of a known contaminant concentration (background or higher) is given over a prolonged time of several weeks, with the intention of approaching a steady state between contaminant intake and elimination through milk and other excretion pathways(Reference Thomas, Sweetman and Jones65Reference Winters, Fries and Lorber68,Reference Fries, Paustenbach and Mather72,Reference Fries, Paustenbach and Luksemburg73,Reference McLachlan, Thoma and Reissinger109Reference McLachlan111,Reference McLachlan and Richter113,Reference Driesen, Lerch and Siegenthaler114,Reference Fries128) . Furthermore, the transfer of contaminants into milk has been studied after intravenous(Reference Safe, Platonow and Hutzinger120) and intramuscular(Reference Vrecl, Ursic and Pogacnik125) applications. The dose of PCDD/Fs and PCBs administered in feeding studies has decreased since the 1970s, likely driven by improvements in analytical detection capabilities (e.g. limit of quantification).

Table 2. Studies used in assessment of transfer rates (TR), transfer factors (TF) and biotransfer factors (BTF) from feed to cow milk in chronological order

* Converted from weeks to days by multiplying by 7.

Finally, at least four studies were available(Reference Travis and Hattemer-Frey136Reference Webster and Connett139) where a theoretical approach was used to assess the exposure and transfer of PCDD/Fs and PCBs in dairy cattle. Although different exposure scenarios in the available studies provide limited comparability of data, all these experimental approaches have their specific advantages to gain insight into the fate of PCDD/F and PCB congeners in dairy cattle and to develop predictive models.

Cattle breeds and other ruminants

Studies that presented feed to milk transfer parameters were based on various lactating dairy cows and beef cattle. Specific cattle breeds were named in some studies, but not in all (Table 1). On the basis of their specific genetic potential, the milk yield differs between cattle breeds. Owing to breeding measures over the past decades, these performance differences have become more evident for specialised dairy cow breeds like Holstein cattle(Reference VandeHaar and St-Pierre140). The potential influences of several breed-associated performance and metabolic parameters on the transfer of xenobiotics in milk of ruminants have already been shown in a modelling approach(Reference MacLachlan141). The authors of this study found that the elimination half-lives of xenobiotics can decrease significantly with increasing milk yield and increasing milk fat content, which was particularly evident for substances with slow liver metabolic rates. But, to date, no systematic study has been carried out to evaluate a potential influence of dairy cattle breed on PCDD/F and PCB transfer.

The extrapolation of the data from previous studies with lower milk-yield breeds to high-performance dairy cows needs to be addressed in the future and should be adopted in toxicokinetic models. In addition, the transferability of data might be further limited by other factors such as metabolic state, body weight and body fat content as well as milk and milk fat production, and will be discussed below.

In addition to cattle, studies have also been carried out with other ruminants, including dairy sheep(Reference Vrecl, Ursic and Pogacnik125), dairy goats(Reference Safe, Platonow and Hutzinger120,Reference Sips, van Eijkeren and Bode124,Reference Ruoff and Fouzy126,Reference Arstila, Reggiani and Sorvari142Reference Fournier, Rychen and Marchand148) and dairy buffaloes(Reference Diletti, Ceci and Ippoliti149). While the milk fat content of goats is very similar to that of cattle, it is more than twice as high in the milk of sheep and buffalo(Reference Park and Haenlein150). Regarding the lipophilic properties of PCDD/Fs and PCBs, these different milk fat contents could influence the transfer rates from feed to milk. Even though there are currently no experimental data available on this aspect, the composition of milk should be considered comparing transfer data in different ruminant species. Furthermore, Rychen et al. (2008) hypothesised goats to be a valuable ‘lactating animal model’ with a transfer behaviour of PCDD/Fs into milk similar to cows(Reference Rychen, Jurjanz and Toussaint52). Therefore, the aforementioned studies on goats will be referenced in the following if appropriate.

Data integrity

A crucial point in evaluating and comparing transfer studies is how and what data were published, as well as how to handle ambiguous or missing data. This is especially important when essential data for transfer modelling, such as milk and milk fat yields, are not reported. Furthermore, data expression in a variety of units and dimensions makes harmonisation challenging. Most studies expressed the amount of PCDD/F and PCB on the basis of milk fat corresponding to regulatory maximum residue levels or daily excretion via milk. However, in some studies the amount was expressed on a whole milk basis(Reference Jensen and Hummel77,Reference Platonow, Funnell and Bullock95,Reference Arnott, Bullock and Platonow98,Reference Perry, Everson and Hendrix115,Reference Arstila, Reggiani and Sorvari142,Reference Willett151) , not necessarily providing information about the respective milk fat content.

Another important aspect of data integrity is the length of the exposure period. During the initial first weeks of exposure, contaminant levels in milk fat increase rapidly, while after weeks of exposure, the increase flattens to slowly approach a steady state. Transfer parameters calculated at the end of a short exposure period will lead to an underestimation compared with steady-state or near-steady-state values. With a constant intake of PCDD/Fs and PCBs, a near steady state was postulated in lactating cows after 17 and 21 d for PCDDF and PCB TEQs(Reference Huwe and Smith82,Reference Huwe and Smith135) , 40–60 d for Aroclor 1254(Reference Fries70) and for specific congeners after 14–42 d(Reference Lorenzi, Angelone and Ferretti86), 21–28 d(Reference Ewers and Heeschen89,Reference Ewers, Heeschen and Reichmuth131) , 28 d(Reference Fries, Paustenbach and Mather72,Reference Fries, Paustenbach and Luksemburg73) , 33 d(Reference Hoogenboom, Traag and Fernandes13) and 40–70 d(Reference Blüthgen, Heeschen and Ruoff134). Predictive modelling suggested it would take approximately 150–200 d to reach a real steady state(Reference Hoogenboom, Traag and Fernandes13,Reference Hoogenboom, Zeilmaker and van Eijkeren42) . However, this might not be reachable owing to metabolic changes during lactation(Reference MacLachlan and Bhula15,Reference Brambilla, Fochi and Falce41) . Therefore, an exposure period of at least 28 d is suggested for transfer studies of lipophilic persistent organic pollutants(Reference MacLachlan and Bhula15,Reference MacLachlan141) . The time to reach steady state for PCDD/Fs and PCBs was also calculated for dairy buffaloes(Reference Diletti, Ceci and Ippoliti149) and goats(Reference Costera, Feidt and Marchand144,Reference Ounnas, Feidt and Toussaint146,Reference Fournier, Rychen and Marchand148) .

Owing to limited analytical performance and limited toxicological assessments, early studies did not differentiate between single congeners. Instead, commercially available PCB-mixtures, such as Aroclor 1254, were often used as a reference for the quantification of PCBs, and results were expressed as Aroclor 1254 equivalents(Reference Fries, Marrow and Gordon71,Reference Gardner, Righter and Roach75,Reference Willett and Liu78,Reference Willett, Liu and Durst81,Reference Platonow, Funnell and Bullock95,Reference Saschenbrecker, Funnell and Platonow96,Reference Arnott, Bullock and Platonow98,Reference Perry, Everson and Hendrix115,Reference Perry, Everson and Hendrix119,Reference Willett151) . Furthermore, some of the aforementioned studies with Aroclor 1254 were later re-evaluated to adjust for different approaches in quantification(Reference Willett, Liu and Fries121). In earlier studies, another analytical challenge could have been the inefficient chromatographic performance of the available gas chromatography columns to fully separate isomers(Reference McLachlan, Thoma and Reissinger109,Reference McLachlan110) .

In literature, data for single congeners are sometimes omitted in favour of the sum of the parameters. The WHO2005 TEQ(Reference Van den Berg, Birnbaum and Denison26) describes our current understanding of the relative toxicity of the congeners and has its justification in risk analysis, being currently used for action and maximum levels in feed and food. Nevertheless, any future changes in the TEF for a single congener will practically invalidate entire study outcomes for which the congener-specific data are not available. Furthermore, as transfer parameters are congener specific with a wide range of variation, transfer data derived for mixtures like Aroclor 1254, or based on sum-TEQ values, are restricted to the respective congener profiles present in the contamination source. Additionally, owing to the congener-specific transfer parameters, the congener profile in milk fat is not identical to that of the source. For these reasons, future publications are advised to report congener-specific data.

Another critical point in assessing previous studies is that transfer parameters are sometimes depicted graphically but not necessarily in text or tables(Reference Fries, Paustenbach and Mather72,Reference McLachlan and Richter113,Reference Fürst, Krause and Hein122,Reference Costera, Feidt and Marchand144) . However, graphical depictions of milk fat concentrations of PCDD/Fs and PCBs during exposure and depuration could be used retrospectively to derive transfer parameters if these were not given in papers, for example Refs.(Reference Slob, Olling and Derks62,Reference Fries70,Reference Fries, Marrow and Gordon71,Reference Furr, Mertens and Gutenmann74,Reference Gardner, Righter and Roach75,Reference Saschenbrecker, Funnell and Platonow96,Reference Arnott, Bullock and Platonow98,Reference Traag, Mengelers and Kan105,Reference Arstila, Reggiani and Sorvari142) . For this review, we extracted data from plots using tools like the WebPlotDigitizer(Reference Rohatgi152), which adds another level of uncertainty. While it is certainly possible to calculate transfer parameters from the tabulated data given in some studies(Reference Piskorska-Pliszczynska, Maszewski and Mikolajczyk43,Reference Jones, Safe and Morcom79,Reference Ewers88,Reference Ruoff91,Reference Platonow, Funnell and Bullock95,Reference Derks, Berende and Olling153) , we saw those as derivative calculations and thus beyond the scope of this review. Aside from milk fat concentrations, some studies include PCDD/F and PCB concentrations in blood, faeces, body fat compartments or organs, enabling a better understanding of the kinetics, which can also be exploited with toxicokinetic modelling. On the basis of these limitations, it is clear that comparison of data from different studies is not trivial.

Kinetic parameters to characterise the feed-to-milk transfer behaviour

Studies with regard to feed-to-milk transfer of PCDD/Fs and PCBs were selected on the basis of the following criteria on study design and data quality:

Of the initial 104 primary and secondary references (of which seventy-six are primary), nineteen remained for further consideration to evaluate congener-specific feed-to-milk transfer parameters after using the above-mentioned criteria. The majority were feeding studies performed under controlled conditions (Table 2). Two further studies analysed the transfer in herds of twenty-six and 1604 animals, respectively. The first study dealt with the background contamination via fresh grass(Reference Kerst, Waller and Reifenhäuser106), while the second was a case study of a contamination incident via food supplements(Reference Brambilla, Fochi and Falce41). While daily contaminant input and feed consumption could not be as well quantified in these studies as in feeding studies, both provide average transfer parameters based on large numbers of individuals.

Exact knowledge of the contaminant exposure amount and its duration is not necessary to derive half-lives during the depuration phase. Because data on milk elimination half-lives are very limited, data from other ruminants like goats (data points marked in Figs. 4 and 5 with a red dot) were included if available (Table 3), as these are also seen as ‘lactating animal models’ elsewhere(Reference Rychen, Jurjanz and Toussaint52). Whenever possible, data from individual animals instead of mean values from a group were used for plotting and discussing the transfer parameters: transfer rate, transfer factor, biotransfer factor and elimination half-lives. However, data reported as greater or smaller than a certain value were neglected. Whenever a range was reported, the arithmetic mean of the range was calculated.

Table 3. Studies used in assessment of α and β elimination half-lives in milk of ruminants in chronological order

Transfer rates

Transfer rates (TR), historically called carry-over rates or carry-over ratio, describe the fraction of congener intake with the diet (mass or mole) that is excreted with the milk. While TRs can be calculated for any given time period during an experiment or an incident, they reach a maximum when a steady state between constant intake and output is reached. The TR is used to assess the balances of the mass flow.

(1) $${\it{TR}}\left[ {\rm{\% }} \right] = {{{\it{DailyExcretionViaMilk}}\left[ {{{{\it{ng}}} \over {\it{d}}}} \right]} \over {{\it{DailyIntakeViaFeed}}\left[ {{{{\it{ng}}} \over {\it{d}}}} \right]}} \times 100{\rm{\% }}$$

Re-arranging this formula enables prediction of the daily excretion via milk for a given TR and a known daily intake via feed, that is,

(2) $$DailyExcretionViaMilk\left[ {{{ng} \over d}} \right] = TR\left[ \% \right]/100\% \times DailyIntakeViaFeed\left[ {{{ng} \over d}} \right]$$

This prediction from the re-arrangement of equation (1) can be seen as the simplest kind of model(Reference Numata, Schafft and Spolders154) and is only valid if the steady state has been reached and the conditions of the prediction are the same as in the determination of TR. TR is a ratio of amounts of contaminant, and is thus insensitive to the basis (e.g. dry matter for feed or milk fat for milk). Most studies expressed the transfer of contaminants into milk as TRs, and thus for all PCDD/F and PCB congeners TRs are available, as seen in the boxplot of Fig. 1. The statistical values behind this boxplot can be found in Supplementary Materials Part A, Table S1. The data used for this analysis are derived from the transfer studies listed in Table 2. Owing to differences in the reported transfer rates for the same feeding study, only the more recent parameters for PCDD/Fs and PCBs by McLachlan(Reference McLachlan, Thoma and Reissinger109Reference McLachlan111) were used for data analysis. We have also included transfer rates with a negative sign, such as those reported by Brambilla et al. (2008)(Reference Brambilla, Fochi and Falce41) and Lorenzi et al. (2020)(Reference Lorenzi, Angelone and Ferretti86). Although physically meaningless, a negative transfer rate can result from a background exposure correction containing measurement and other errors. Thomas et al. (1999) similarly reported a transfer rate above 100% for PCB-118, which will be discussed at a later point in this review(Reference Thomas, Sweetman and Jones66).

Fig. 1. Boxplots of transfer rates (TRs) for cows. Scatter points represent available data points from selected literature (refer to Table 2). Boxes (N > 5) are defined as the interquartile range (IQR) between 25th percentile (Q1) and 75th percentile (Q3) of the data according to the standard method. The black line in the box represents the median. Whiskers include data within 1.5 times of IQR below Q1 and above Q3. Plot generated with Python 3.10 using the Seaborn, Matplotlib, Numpy and Pandas libraries.

According to this analysis, TRs of highly chlorinated ndl-PCBs and dl-PCBs are, with some exceptions, higher compared with PCDD/Fs. Altogether, congener-specific TRs lie in the range between <1% and 100%. Despite the consultation of data from comparable studies (Table 2), coefficients of variation (CV) were between 29% and 204% (Table S1), with a trend of decreasing CV at higher TRs. The wide range of calculated TRs (Fig. 1) can be due to several reasons. For instance, the data of certain studies were of low accuracy, owing to low analytical performance at lower TRs; important experimental parameters were not considered, such as contaminant source-dependent bioavailability; or the inter-individual variance between animals was very high for unknown reasons. Moreover, for several congeners, only a very limited number of studies investigating a small number of animals exist, which are thus error prone. In any case, for the purpose of modelling, the large variability in the data negatively affects the predictive power.

Transfer factors

For regulatory purposes and risk assessment, the transfer factor (TF), also known as bioconcentration factor (BCF), carry-over factor (CoF) or accumulation factor (AF), has played a key role until now. This dimensionless factor describes the quotient of the congener concentration in milk (fat) and its concentration in the feed.

(3) $$TF = {{{\it{ConcentratonInMilkFat}}\left[ {{{{\it{ng}}} \over {{\it{kg}}}}} \right]} \over {{\it{ConcentrationInFeed}}\left[ {{{{\it{ng}}} \over {{\it{kg}}}}} \right]}}{\rm{\;}}$$

TF is a ratio of concentrations of a contaminant, and is thus dependent on the basis used. The standard basis from animal nutrition is compound feed with 88% dry matter for feed and milk fat for milk, as used in regulatory standards. The numerical value of TF will vary if other bases are used (e.g. wet weight or just mineral components of compound feed). TF has the disadvantage of having a strong dependence on, for example, milk fat yield. Furthermore, the concept of TF can be applied to describe the transfer and storage of contaminants into fat tissue. This formula can also be re-arranged so that it can be used to predict the concentration in milk fat given concentration in feed and the TF, that is,

(4) $$ConcentrationInMilkfat\left[ {{{ng} \over g}} \right] = TF\; \times \;ConcentrationInFeed\left[ {{{ng} \over g}} \right]$$

TR and TF are related through the feed efficiency (FE), or the ratio of milk fat produced per kg of feed (88% dry matter), such that:

(5) \begin{align}{\it{TR}}\left[ {\rm{\% }} \right] &= {\it{TF}} \times {\it{FE}}\\ &= {{{\it{ConcentrationInMilkfat}}\left[ {{{{\it{ng}}} \over {{\it{kg}}}}} \right]} \over {{\it{ConcentrationInFeed}}\left[ {{{{\it{ng}}} \over {{\it{kg}}}}} \right]}} \times {{{\it{milkfat}}\left[ {{{{\it{kg}}} \over {\it{d}}}} \right]} \over {{\it{feed}}\left[ {{{{\it{kg}}} \over {\it{d}}}} \right]}}\\ &\quad \times 100\% \end{align}

While the conversion between TR and TF is given by equation (5), the data in each study do not always allow for it to be effectively used. In this review, we have opted for only reproducing published values. TFs are less well documented than TRs. For example, there are no data given for some dl-PCBs, and for the lower chlorinated ndl-PCBs there is only one publication by Ewers (1987)(Reference Ewers88) that calculated TFs on the basis of a study by Tuinstra et al. (1981)(Reference Tuinstra, Vreman and Roos57). While Tuinstra et al. (1981) published enough data to calculate TFs, no data were given for PCB-28. It is therefore unknown how Ewers (1987) calculated the TF value, and thus this value was omitted. Furthermore, Connett and Webster (1987)(Reference Connett and Webster58), McLachlan et al. (1990)(Reference McLachlan, Thoma and Reissinger109) and Huwe and Smith (2003)(Reference Huwe and Smith82) calculated TFs on the basis of wet weight feed intake, which is not comparable to dry weight intake. Hence, these TFs were not used for further data analysis.

TFs based on studies listed in Table 2 are depicted in Supplementary Materials Part A, Table S2 and plotted in Fig. 2. The coefficient of variation of congener-specific TFs lies between 15% and 154%, again with a trend of lower CV for congeners with higher TFs. The reason for the wide range of calculated TFs (Fig. 2) may be similar to those presented above for TRs.

Fig. 2. Boxplots of transfer factors (TFs) in 88% dry matter basis for the feed and milk fat basis for the milk of cows. Scatter points represent available data points from selected literature (refer to Table 2). Boxes (N > 5) are defined as the interquartile range (IQR) between 25th percentile (Q1) and 75th percentile (Q3) of the data according to the standard method. The black line in the box represents the median. Whiskers include data within 1.5 times of IQR below Q1 and above Q3. Plot generated with Python 3.10 using the Seaborn, Matplotlib, Numpy and Pandas libraries.

Biotransfer factors

Another frequently reported parameter is the biotransfer factor (BTF). Contrary to TF, and similarly to TR, it is not restricted to an exposure from a single source (e.g. feed) but can also account for contamination through multiple pathways such as air, wall paint chippings, soil, etc. BTF is a ratio of concentration in whole milk and the daily intake (dose) of a contaminant. Usually, the BTF is calculated on the basis of whole milk instead of milk fat, deviating from the standard for TF. Moreover, the BTF is not dimensionless and has units of time/mass, such as d/kg. In non-physiological approaches (see part II of this review(Reference Moenning, Krause, Lamp, Maul, Schenkel, Fürst, Pieper and Numata56)) feed-to-milk BTFs are estimated on the basis of substance-specific octanol-water partition coefficients. However, BTF is more commonly applied in, for example, aquatic ecotoxicology (where the contaminant sources are multiple and diffuse) and is calculated as

(6) $${\it{BTF}}\left[ {{{\it{d}} \over {{\it{kg}}}}} \right] = {{{\it{ConcentrationInMilk}}\left[ {{{{\it{ng}}} \over {{\it{kg}}}}} \right]} \over {{\it{DailyIntake}}\left[ {{{{\it{ng}}} \over {\it{d}}}} \right]}}.$$

Again, this formula can be re-arranged, so that we can predict the concentration in milk given the daily intake and the BTF, that is,

(7) $$ConcentrationInMilk\left[ {{{ng} \over {kg}}} \right] = BTF\left[ {{d \over {kg}}} \right] \times DailyIntake\left[ {{{ng} \over d}} \right].$$

When feed is the only contamination source, the BTF in whole milk basis is related to TF in milk fat basis and TR by

(8) $${\it{TF}} = {\it{BTF}}\left[ {{{\it{d}} \over {{\it{kg}}}}} \right]/{\it{FatfractionInMilk}} \times {\it{feed}}\left[ {{{{\it{kg}}} \over {\it{d}}}} \right]$$
(9) $$\;{\it{TR}}\left[ {\it{\% }} \right] = {\it{BTF}}\left[ {{{\it{d}} \over {{\it{kg}}}}} \right] \times {\it{milk}}\left[ {{{{\it{kg}}} \over {\it{d}}}} \right] \times 100\% $$

In this review, we have opted not to transform BTF into TF or TR because of incomplete reported data (since this would have required assumptions). Data on BTFs are even more sparse than for TFs, and as shown in Fig. 3, BTFs are missing for most dl-PCBs. Rosenbaum et al. (2009)(Reference Rosenbaum, McKone and Jolliet130) retroactively calculated BTFs on the basis of various studies(Reference Tuinstra, Vreman and Roos57,Reference Thomas, Sweetman and Jones65,Reference Thomas, Sweetman and Jones66,Reference Firestone, Clower and Borsetti76,Reference Heeschen and Blüthgen87,Reference McLachlan, Thoma and Reissinger109Reference McLachlan111,Reference McLachlan and Richter113) , which were also used in the depiction of BTFs in Fig. 3. Average congener-specific BTFs given in a report by the California Office of Environmental Health Hazard Assessment(155) were omitted, as these values were based on several experimental studies(Reference Slob, Olling and Derks62,Reference Schuler, Schmid and Schlatter63,Reference Thomas, Sweetman and Jones66,Reference Lorber, Fries and Winters67,Reference Thomas, Sweetman and Lohmann85,Reference Kerst, Waller and Reifenhäuser106,Reference McLachlan, Thoma and Reissinger109,Reference McLachlan and Richter113,Reference Huwe and Smith135) . The data behind Fig. 3 can be found in Supplementary Materials Part A, Table S3. The coefficients of variation varied between 22% and 140%, which is in line with the previous observations for TFs and TRs. The reasons for this variability, also apparent in Fig. 3, may be similar to those presented above for TRs.

Fig. 3. Boxplots of biotransfer factors (BTFs) in whole milk basis for the milk of cows. Scatter points represent available data points from selected literature (refer to Table 2). Boxes (N > 5) are defined as the interquartile range (IQR) between 25th percentile (Q1) and 75th percentile (Q3) of the data according to the standard method. The black line in the box represents the median. Whiskers include data within 1.5 times of IQR below Q1 and above Q3. Plot generated with Python 3.10 using the Seaborn, Matplotlib, Numpy and Pandas libraries.

Fig. 4. Boxplots of α elimination half-lives. Scatter points represent available data points from selected literature (refer to Table 3) mostly for cows, with red dots indicating data from goats. Boxes (N > 5) are defined as the interquartile range (IQR) between 25th percentile (Q1) and 75th percentile (Q3) of the data according to the standard method. The black line in the box represents the median. Whiskers include data within 1.5 times of IQR below Q1 and above Q3. Plot generated with Python 3.10 using the Seaborn, Matplotlib, Numpy and Pandas libraries.

Fig. 5. Boxplots of β elimination half-lives. Scatter points represent available data points from selected literature (refer to Table 3), mostly for cows with red dots indicating data from goats. Boxes (N > 5) are defined as the interquartile range (IQR) between 25th percentile (Q1) and 75th percentile (Q3) of the data according to the standard method. The black line in the box represents the median. Whiskers include data within 1.5 times of IQR below Q1 and above Q3. Plot generated with Python 3.10 using the Seaborn, Matplotlib, Numpy and Pandas libraries.

Lastly, an appeal is presented to the next generation of risk assessors and risk managers: resist the temptation to use TF, TR and BTF as simple multiplicative factors to translate between maximum levels in food and maximum levels of contaminants in contaminated source materials (e.g. feed). Opt instead for full toxicokinetic models that capture the dynamics, physiological and metabolic variables explicitly (see part II of this review(Reference Moenning, Krause, Lamp, Maul, Schenkel, Fürst, Pieper and Numata56)). For all their disadvantages, these transfer parameters nevertheless remain a good way of comparing studies, as we do in Figs. 13.

Elimination half-lives

Knowledge of the elimination half-lives or rate constants is of utmost importance for risk analysis. Elimination half-lives are necessary to estimate depuration time needed to reach concentrations in milk which are compliant with regulatory requirements after a contaminant exposure incident in a herd. At this point, an explanation is provided as to why multiple elimination half-lives are observed when the dairy cows are allowed a depuration phase after a period of exposure to lipophilic contaminants. The faster α elimination half-life dominates initially because the contaminants are readily available from the blood to be transferred into milk. This is sometimes called ‘elimination phase’ or ‘initial phase’. Any further flow of lipophilic contaminants into milk has to come from contaminants remobilised from fat, which is a slow compartment and makes the terminal β elimination half-life slower (slower = long in time units). This is referred to as the ‘distribution phase’ or ‘terminal phase’. The term ‘distribution phase’ is not recommended for the β phase, since it coincides with the name given in the pharmacokinetics of medicinal compounds to the α phase. To better understand the physical origin of the α and β elimination half-lives, please refer to part II of this review on Toxicokinetic Predictive Models for Risk Assessment(Reference Moenning, Krause, Lamp, Maul, Schenkel, Fürst, Pieper and Numata56).

Unfortunately, most studies neglect inclusion of a depuration period, as animals were slaughtered directly after exposure to assess the distribution of contaminants into various body compartments. Even if a depuration phase was included in an experiment, derived kinetics were often limited to the later and slower β phase, as the time resolution was not sufficient to observe the initial fast elimination phase. Hence, kinetic parameters for fast α elimination phase are not available for most congeners.

α half-lives derived from the initial fast elimination phase are plotted in Fig. 4, with the statistical data listed in Supplementary Materials Part A, Table S4. However, the information on α half-lives is sparse at best. For the majority of the PCDD/F congeners, only two studies contained values for the α elimination phase(Reference Driesen, Lerch and Siegenthaler114,Reference Huwe and Smith135) , whereas kinetics for 1,2,3,7,8-PeCDF, 1,2,3,4,7,8,9-HpCDF and several PCBs is still unknown. Furthermore, the only available data for 2,3,7,8-TCDF α half-lives, as well as an indication for PCB-101 (<4 d) in dairy animals, are derived from a study with four alpine goats(Reference Fournier, Rychen and Marchand148). Three congeners were studied in more detail, namely 2,3,7,8-TCDD, 1,2,3,7,8-PeCDD and 2,3,4,7,8-PeCDF(Reference Roos, Berende and Traag61,Reference Jones, Safe and Morcom79,Reference Ruoff91,Reference Driesen, Lerch and Siegenthaler114,Reference Tuinstra, Roos and Berende116,Reference Huwe and Smith135,Reference Fournier, Rychen and Marchand148) . According to these data, half-lives range from less than 1 d for OCDF and up to 2 weeks for 1,2,3,4,6,7,8-HpCDD. Furthermore, the data show high variation, especially for 2,3,7,8-TCDD.

β half-lives are distinctly longer than α half-lives, ranging from days to several months, and are plotted in Fig. 5 on the basis of statistical data depicted in Supplementary Materials, Table S5. Again, data for several dl-PCBs, ndl-PCB-28 and 1,2,3,7,8-PeCDF are unavailable. Furthermore, data for OCDF, PCB-169, PCB-52 and PCB-101 are available from only one feeding(Reference Huwe and Smith135) and one field study(Reference Tremolada, Guazzoni and Parolini108), respectively. Additional data for cows were derived from feeding studies(Reference Firestone, Clower and Borsetti76,Reference Huwe and Smith82,Reference Ewers88,Reference Ruoff91,Reference Driesen, Lerch and Siegenthaler114,Reference Huwe and Smith135) , single exposure studies(Reference Jensen and Hummel77,Reference Olling, Derks and Berende102) , exposure during the dry phase(Reference Roos, Berende and Traag61,Reference Ruoff91) , field studies(Reference Derks, Berende and Olling104,Reference Tremolada, Guazzoni and Parolini108) and incident studies(Reference Malisch35,Reference Brambilla, Fochi and Falce41) and for alpine goats(Reference Fournier, Rychen and Marchand148) as listed in Table 3. Half-lives from the study by Roos et al. (1991)(Reference Roos, Berende and Traag61) were re-evaluated by Tuinstra et al. (1992)(Reference Tuinstra, Roos and Berende116), which were used for the present data analysis. Half-lives expressed in weeks(Reference Malisch35) or months(Reference Klein, Drochner and Forschner100) were converted to days by multiplying the value by 7 or 30, respectively.

The heterogeneity of the experimental conditions and methodical approaches under which the data have been derived could be one reason for the variability evident in Fig. 5. Furthermore, the high variance for PCB-180 could be explained with the inclusion of β half-lives from goats, which are shorter compared with half-lives derived from cows(Reference Klein, Drochner and Forschner100,Reference Tremolada, Guazzoni and Parolini108) . Fournier et al. (2013) explained this with the small body fat compartment and the relatively high milk fat yield seen in goats compared with cows(Reference Fournier, Rychen and Marchand148).

Aside from congener-specific data, α half-lives for sum parameters such as Aroclor 1254(Reference Fries70,Reference Fries, Marrow and Gordon71,Reference Willett and Liu78,Reference Saschenbrecker, Funnell and Platonow96) , ∑PCDDF/s, ∑dl-PCBs and ∑ndl-PCBs(Reference Fournier, Rychen and Marchand148), I-TEQ(Reference Roos, Berende and Traag61), WHO1998-PCDD/F-TEQ and WHO1998-dl-PCB-TEQ(Reference Huwe and Smith82,Reference Huwe and Smith135) as well as WHO2005-PCDD/F-TEQ(Reference Hoogenboom, Zeilmaker and van Eijkeren42) are available. Furthermore, β half-lives are available for Aroclor 1254(Reference Fries70,Reference Fries, Marrow and Gordon71,Reference Willett and Liu78,Reference Perry, Everson and Hendrix119,Reference Hendriks, Smítková and Huijbregts156) , ∑PCDDF/s and ∑dl-PCBs(Reference Fournier, Rychen and Marchand148) ∑ndl-PCBs(Reference Ewers88,Reference Fournier, Rychen and Marchand148) , I-TEQ(Reference Malisch35,Reference Roos, Berende and Traag61,Reference Olling, Derks and Berende102,Reference Traag, Mengelers and Kan105) , WHO1998-PCDD/F-TEQ and WHO1998-dl-PCB-TEQ(Reference Huwe and Smith82,Reference Huwe and Smith135) as well as WHO2005-PCDD/F-TEQ(Reference Brambilla, Fochi and Falce41,Reference Hoogenboom, Zeilmaker and van Eijkeren42) .

Proposing transfer parameters for the sum of congeners may be useful for practical applications, but as sum parameters (such as TEF) change, these models lose their relevance. Furthermore, congener profiles in foods of animal origin and feed can be used as forensic tools to identify contamination sources(Reference Hoogenboom, Malisch and van Leeuwen14). Congener-specific reporting is mandatory for risk analysis.

A number of transfer parameters for this substance class can be summarised, but they are rarely all reported in a single study. As a consequence, for the parameters mentioned, not many reliable data points are available. In addition, the parameters display high variability. Possible factors influencing transfer parameters in experimental studies are discussed in depth according to the cow’s metabolic state and the contaminant properties in the following chapters.

Factors influencing transfer parameters in experimental studies

Factors that enhance or restrict the transfer from feed to milk can be classified into two categories. The first category is the cow with its individual metabolic state, health and performance parameters(Reference MacLachlan141). The second category is determined by the contaminants, their physico-chemical properties and their interaction with the animal metabolism(Reference Bonate157).

Factors stemming from the cow’s metabolic state

Especially in high-performance cows, the energy requirements cannot be completely covered by the feed intake during early lactation. The resulting metabolic state of negative energy balance (NEB) needs to be offset by mobilisation of body reserves (mainly body fat). The tides are turned approximately 80–120 d after calving into a metabolic state with positive energy balance (PEB) thanks to increased feed intake capacity, leading to growing body fat reserves during later lactation. Thus, the contribution of mobilised body fat to total milk fat production is proportionally higher during the NEB in early lactation. During the later lactational stages in PEB, mainly de novo fat synthesis from ruminally produced short chain fatty acids and (in part) dietary fatty acids are the main source of milk fat(Reference Thomas, Sweetman and Jones65,Reference Stoop, Bovenhuis and Heck158) . Body fat mobilisation in NEB and weight gain in PEB are assumed to be of major importance when discussing the release and deposition of lipophilic contaminants in lactating ruminants(Reference Rychen, Jurjanz and Toussaint52,Reference Rychen, Jurjanz and Fournier54,Reference Rossi, Bertuzzi and Vitali123,Reference McLachlan159) . On the one hand, during NEB lipophilic contaminants might be released from body fat depots and excreted via milk(Reference Thomas, Sweetman and Jones66,Reference Klein99) . On the other hand, in PEB initial concentrations of contaminants in the body fat would be ‘diluted’ in the fat compartment, and the reservoir for storage of these compounds increases(Reference Thomas, Sweetman and Jones66). Furthermore, distribution kinetics between blood and body fat depend on the blood/fat surface area, which increases with growth of body fat. A greater blood/fat surface area supposedly increases the transfer coefficient for blood/fat diffusion of lipophilic contaminants and affects their elimination through milk(Reference McLachlan159). As a consequence, differentiation of lactational stages into ‘early’, ‘mid’ and ‘late’ should be made to understand possible differences in transfer rate and α or β elimination half-lives (Table 4).

Table 4. Overview of the lactational status defined by days in milk (DIM) of cows in selected studies at beginning to end of the exposure phase to discuss the influence on transfer parameters

* Calculated from month p.p. multiplied by 30 d.

However, the influence of the metabolic state on the transfer of PCDD/Fs and PCBs into milk is still ambiguous in view of several field observations and feeding studies. For example, no increase of PCDD/Fs in milk fat was observed during nutritionally induced body fat mobilisation(Reference Schulz, Wiesmüller and Appuhn94). Three studies(Reference Thomas, Sweetman and Jones66,Reference Willett and Liu78,Reference Ewers88) analysed the influence of NEB and PEB on the transfer of PCBs to milk during lactation. Thomas et al. (1999) conducted a long-term mass balance study in which the cows transitioned metabolic phases and described higher TRs for several PCBs during NEB than in PEB phase(Reference Thomas, Sweetman and Jones66). Similarly, Willet and Liu (1982) observed that excretion of PCBs was elevated for cows early in lactation compared with cows late in lactation(Reference Willett and Liu78). In contrast, Ewers (1987) described divergent effects of the lactational stage on the transfer of ndl-PCBs at two different dosage levels(Reference Ewers88). At the lower dosage, the TRs of ndl-PCBs were lower in early lactation, while no effect from the lactation period was evident with a 67% higher dose of the same ndl-PCB mixture on TRs, TFs and elimination half-lives. On the basis of the current data, no conclusive statement can be made on the influence of the lactational stage.

Regarding milk elimination half-lives, Brambilla et al. (2008) observed very short β half-lives compared with other studies in a herd of dairy cows after PCDD/F contamination. The authors attributed these findings to the physiological status of the herd, with 25% of the cows likely being in a negative energetic balance(Reference Brambilla, Fochi and Falce41). Fries et al. (1972) observed a tendency towards lower mean first-order rate constants for Aroclor 1254 depuration (half-lives were greater) for cows in early lactation compared with mid- and late lactation (Table 4)(Reference Fries70). However, no statistics were analysed to underline this statement. Importantly, a follow-up feeding study with nine cows in different stages of lactation saw no consistent relation between the first-order rate constants and the rate of body weight change(Reference Fries, Marrow and Gordon71).

Depuration kinetics of PCBs were further studied under the influence of a thyroprotein supplement, which induced body fat mobilisation and increased milk and milk fat production. During the thyroprotein treatment, elimination of PCBs in milk fat increased compared with a control group(Reference Willett and Liu78). Furthermore, Farries (1990) studied the excretion of ndl-PCBs in fifty dairy cows fed according to their energetic needs. During the depuration phase, half of the cows were purposely undersupplied by 30%, resulting in weight loss and a decrease in milk yield. The resulting fat mobilisation also caused a slight initial increase of PCB content in milk fat, whereas the PCB milk fat content in the energetically balanced cows further decreased(Reference Farries59).

During the dry period before lactation, the most important elimination pathway for the excretion of PCDD/Fs and PCBs through milk is not available, leading to increased bioaccumulation in body fat of the pregnant cow and the foetus. Several feeding studies described the exposure of animals during the dry period and observed transfer and elimination kinetics in milk postpartum with different results. Several studies observed a rapid decline of PCBs and PCDD/Fs in colostrum(Reference Ewers88,Reference Schulz93,Reference Schulz, Wiesmüller and Appuhn94,Reference Klein, Drochner and Forschner100,Reference Tuinstra, Roos and Berende116) , whereas Hirako (2008) described a decline for most PCDD/Fs congeners while the dl-PCB content in colostrum and milk lipids was nearly constant(Reference Hirako127). Klein et al. (1992) stated that high milk yields promoted fast depuration postpartum(Reference Klein, Drochner and Forschner100), and Tuinstra et al. (1992) postulated that the depuration kinetics for PCDD/Fs postpartum could be different later in lactation(Reference Tuinstra, Roos and Berende116). Furthermore, Ewers et al. (1987) estimated TRs on the basis of the accumulative input and output, and stated that these TRs were higher in early lactational state compared with TRs derived from cows later in lactation(Reference Ewers88). Differences in body fat dynamics, PCDD/F and dl-/ndl-PCB absorption or metabolic clearance rates may also be attributed to parity(Reference Driesen, Lerch and Siegenthaler114).

Body weight and body fat content

Aside from dynamic changes in body fat due to metabolism, the size of the fat compartments is an important factor for distribution and storage of lipophilic substances like PCDD/Fs and PCBs. However, quantification of fat compartments in a living animal is challenging. Hence, proxies, such as body weight or verbal description of the body fat content, are used to explain the influence of fat compartment sizes on transfer kinetics.

Fries et al. (1973) concluded that there was no connection between body weight of the cows and elimination half-lives in milk(Reference Fries, Marrow and Gordon71). In a follow-up study in 1999, the same authors observed no connection between transfer rates and body weight(Reference Fries, Paustenbach and Mather72), whereas by contrast, McLachlan and Richter (1998) suggested that TR could be influenced by body weight, albeit not strongly(Reference McLachlan and Richter113).

Ewers et al.(Reference Ewers88Reference Ewers, Reichmuth and Wetzel90,Reference Ewers, Heeschen and Reichmuth131) noted that TR increased with decreasing body fat percentage, though milk yields and dosage levels also need to be considered. On the basis of data from the feeding study published by Ewers et al.(Reference Ewers88,Reference Ewers and Heeschen89,Reference Ewers, Heeschen and Reichmuth131) , Heeschen et al. (1993) stated that the biological half-lives of PCBs in cows with an elevated body fat percentage are twice as high compared with animals with an estimated low body fat percentage, regardless of milk fat yield(Reference Heeschen, Ewers and Blüthgen132).

In another study, shorter half-lives for a cow with low body fat compared with three other cows were calculated. The authors postulated that with a lower body fat weight the volume for distribution of lipophilic substances would also be small, and thus the transfer from body fat to milk would be elevated(Reference Roos, Berende and Traag61,Reference Tuinstra, Roos and Berende116) . A similar observation regarding the decline of PCB concentrations in milk of two cows with different body weights was described by Klein (1991)(Reference Klein99).

Milk and milk fat production

The lipid content of milk depends on the lactation stage, individual animal, breed and feed composition. Furthermore, the daily milk fat production is related to the daily milk yield and the metabolic state of the cow. Both of these performance parameters have been steadily increased through breeding measures for dairy cows(160), a trend that also reflects in data of the cited transfer studies (Table 5).

Table 5. Milk and milk fat production of selected studies in chronological order

* Calculated from L/d by multiplying with a density of 1.04 kg/l and rounded to one decimal place.

Owing to the lipophilic nature of PCDD/Fs and PCBs, these compounds have a high affinity to milk lipids. It is therefore reasonable to assume that the daily milk fat yield has an effect on the transfer parameters. McLachlan et al. (1998) stated that the interaction of PCDD/F TRs and lactation rate, among other factors, is not very pronounced(Reference McLachlan and Richter113). Other studies mentioned a positive correlation between milk yield(Reference Fries, Paustenbach and Mather72,Reference Huwe and Smith135) or milk fat yield(Reference Fries, Paustenbach and Luksemburg73) and the TRs of PCDD/Fs into milk. TFs of PCDDF/s were found to be less dependent on milk production, but could be related to milk fat content(Reference Fries, Paustenbach and Mather72).

In a feeding study, total recovery of ndl-PCBs (PCB-138, PCB-153 and PCB-180) in milk was higher in a cow with higher milk yield but also lower body fat percentage compared with another cow with less milk yield and higher body fat percentage(Reference Heeschen and Blüthgen87). In a follow-up experiment with three cows, TRs increased with higher milk yields and lower body fat percentages(Reference Ewers88,Reference Ewers and Heeschen89,Reference Ewers, Heeschen and Reichmuth131,Reference Heeschen, Ewers and Blüthgen132) . These observations coincide with a recent feeding study with primiparous (lower milk fat yield, lower body weight) and multiparous (higher milk fat yield, higher body weight) cows, suggesting that lower milk fat yields could be one reason for lower PCDD/F and dl-/ndl-PCB TRs(Reference Driesen, Lerch and Siegenthaler114).

The influence of milk fat and milk yield on the depuration phase was discussed in several studies. Klein et al. (1992) stated that differences in milk fat production at comparable milk yields did not clearly correlate with PCB excretion rates(Reference Klein, Drochner and Forschner100). Heeschen et al. (1993) assumed that the mean milk fat production had no effect(Reference Heeschen, Ewers and Blüthgen132), and Rossi et al. (2010) concluded that averaged daily milk yields per cow did not affect the rate of PCB excretion(Reference Rossi, Bertuzzi and Vitali123). Furthermore, Fries et al. (1973) found no significant correlation between milk fat production per cow and PCB excretion rates(Reference Fries, Marrow and Gordon71).

However, Klein et al. (1992) also stated that elimination half-lives from cows, with comparable body weight and PCB milk fat content before depuration, were shorter with higher milk yields(Reference Klein, Drochner and Forschner100). Recently, Driesen et al. (2022) derived overall shorter PCB elimination half-lives for multiparous cows with higher milk fat yields compared with primiparous cows. However, a dependency between elimination half-lives calculated for PCDD/Fs and milk fat yield was less conclusive, with half-lives for PCDD/F TEQs being even shorter for primiparous cows(Reference Driesen, Lerch and Siegenthaler114).

Overall, variation of milk fat production within singular studies could be too small to observe significant effects on transfer parameters. Therefore, the influence of this factor cannot be assessed conclusively on the basis of the current data situation.

A final potentially important factor is the milk fat origin in terms of its proportion from de novo synthesised milk fat to remobilised body fat, as discussed above in the section on metabolic state.

Animal health

Little is known on whether animal health status has an impact on the transfer of PCDD/Fs and PCBs into milk. Absorption of PCBs in the gastrointestinal tract might be restricted by diarrhoea as the passage of feed is accelerated. Hence, the transfer into milk and tissue fat is reduced compared with healthy cows(Reference Furr, Mertens and Gutenmann74). Mastitis, on the other hand, might promote the transfer of PCDD/Fs into the milk, as the mammary gland blood–milk barrier becomes more permeable(Reference Fries, Paustenbach and Mather72). Field studies and incident investigations on PCDD/F and PCB contamination rarely mention the animal health status. Thus, a possible influence on the results of the studies included in this review can neither be excluded nor confirmed.

In turn, PCDD/F and dl-PCB exposure may have effects on animal health, such as oxidative damage combined with reduced levels of antioxidants in the blood(Reference Spagnuolo, Sarubbi and Rossetti161,Reference Spagnuolo, Cigliano and Nebbia162) and chromosome fragility(Reference Genualdo, Perucatti and Iannuzzi163). Exposure to elevated levels of TCDD was even suspected to be the cause of abortions and stillbirths in cattle(Reference Davies, Desouza and Waghela164). However, the European Food Safety Authority (EFSA) concluded that all the above-mentioned studies are not fit for risk assessment owing to missing data or exposure to other potential hazardous substances(Reference Knutsen, Alexander and Barregård16).

No adverse effects on the conversion of feed to milk or body weight were observed in dairy cattle exposed to high amounts of 1 g/d (1.7 mg/kg body weight/d(Reference Willett, Liu and Fries121)) Aroclor 1254(Reference Willett, Liu and Durst81) over several days. Furthermore, during in vitro fermentation experiments on micro-organism activity from the rumen of Holstein steers, no effect in terms of dry matter disappearance was observed in presence of different PCB concentrations. The toxic effects of ndl-PCBs on laboratory animals are generally difficult to identify, since they rarely appear in a pure form. Even less information is available on specific effects of ndl-PCBs on the health status of dairy cows. It can be expected that very high doses of pure ndl-PCB 180 might have neurotoxic, endocrine and behavioural effects(Reference Viluksela, Heikkinen and van der Ven165), but such doses are not expected in practice for farm animals.

Factors stemming from the contaminants

Dosage

Owing to substantial progress in analytical techniques, the limits of quantification for PCDD/Fs and PCBs in milk fat have steadily decreased during the past decades. Therefore, smaller daily doses could be applied during feeding studies while ensuring the quality of the derived transfer kinetics. This development is also in accordance with regulatory efforts that helped to reduce the overall exposure of ruminants to PCDD/Fs and PCBs through feed in the past decades. Furthermore, several feeding studies in the 1970s and 1980s applied varying daily concentrations to evaluate the correlation between dosage level and transfer into milk and or milk fat.

Platonow et al. (1971) noticed an approximately ten-fold increase of Aroclor 1254 in whole milk, and a more than ten-fold increase in cream, with a ten-fold increase of the daily dosage(Reference Platonow, Funnell and Bullock95). However, Arnott et al. (1977) stated, on the basis of single- and multiple-dosage studies with Aroclor 1254, that proportionally less PCBs are excreted with the milk at higher doses (higher dose, lower TR)(Reference Arnott, Bullock and Platonow98). The tentative mechanism proposed was a dose-dependent decrease in absorption.

Willet et al. (1990) re-evaluated the quantification of Aroclor 1254 in milk fat from earlier studies and found a significant linear correlation between dosage and PCB concentrations in milk fat(Reference Willett, Liu and Fries121). Contrastingly, an increase of TRs for the three indicator PCBs 138, 153 and 180 was observed at higher exposure and decreasing body fat percentage(Reference Ewers88,Reference Ewers and Heeschen89) .

For PCDD/Fs, only one study stated an increase in 2,3,7,8-TCDD residues in milk and cream in proportion to the daily dose(Reference Jensen and Hummel77). Dose-dependent absorption(Reference Savvateeva, Numata and Pieper166) or self-induced metabolism of the contaminants(Reference Inui, Itoh and Yamamoto167) are two possible mechanisms that could be responsible for the lack of linearity between the oral dose rate and the milk fat correlation response. However, these mechanisms have been proposed for other species; therefore, the conclusions regarding the influence of the dosage are still not clear for dairy cattle.

Contaminant source-dependent bioavailability

Ruminants may be exposed to PCDD/Fs and PCBs via various sources, such as roughage, concentrates and mineral supplements, but also wall paints and wood preservatives, soil, sewage sludge and fly ash. The bioavailability of the contaminants from different sources (matrices, Table 6), that is, primarily the absorption of PCDD/Fs and PCBs by the ruminants, needs to be considered when discussing transfer parameters. Furthermore, PCBs have also been administered intravenously(Reference Safe, Platonow and Hutzinger120) as well as intramuscularly(Reference Vrecl, Ursic and Pogacnik125) to sheep.

Table 6. Transfer studies grouped according to the contamination source

On the basis of two successive mass balance studies, McLachlan et al. (1998) concluded that digestive tract absorption of PCDD/Fs from sewage sludge, as the contamination source for grass silage, was comparable to grass silage contaminated by atmospheric deposition(Reference McLachlan and Richter113). The bioavailability of PCDD/Fs in fly ash collected from a filter of a municipal solid waste incinerator and administered intra-ruminally was lower than the value estimated for grass from a pasture near a large municipal solid waste incinerator(Reference Slob, Olling and Derks62). TFs of PCDD/Fs derived from a feeding study with pentachlorophenol-treated wood(Reference Fries, Paustenbach and Mather72) compared with a study with contaminated magnesium mineral supplement suggested a reduced bioavailability of PCDD/Fs from pentachlorophenol-treated wood(Reference Huwe and Smith82).

Fürst et al. (1993) stated that soil intake of grazing cows is either less than assumed or that the bioavailability of PCDD/Fs from soil is lower than expected(Reference Fürst, Krause and Hein122). A lower bioavailability of several PCDD/Fs and PCBs from soil was also hypothesised by Driesen et al. (2022), who calculated tendentially higher TRs for a control group fed grass-silage compared with another group fed grass-silage mixed with contaminated soil(Reference Driesen, Lerch and Siegenthaler114). However, according to several studies with goats and cows, bioavailability of PCBs from soil is comparable to maize silage(Reference Fournier, Rychen and Marchand148), hay(Reference Ounnas, Feidt and Toussaint146) or grain(Reference Jones, Safe and Morcom80). Compared with a spiked oil matrix, the bioavailability of PCBs from a sandy soil contaminated by a fire is slightly reduced(Reference Feidt, Ounnas and Julien-David147).

The absorption of PCDD/Fs from contaminated lime, when used as a neutraliser for acidic citrus pulp, is comparable to that of grass. During this process, the lime is completely dissolved and the PCDD/Fs are absorbed in the citrus pulp(Reference Malisch35). Slob et al. (1995) hypothesised increased lipid content and smaller particle sizes as factors that could increase bioavailability(Reference Slob, Olling and Derks62). An overview of various TFs according to different matrices as carriers for the contamination is given by Huwe et al. (2005)(Reference Huwe and Smith135), in which the substantially lower TFs for fly ash and grass contaminated with municipal solid waste(Reference Slob, Olling and Derks62) are attributed to particulate binding(Reference Huwe and Smith135). Furthermore, bioavailability of PCDD/Fs and PCBs also depends on the physico-chemical properties of the congeners, such as degree of chlorination.

Degree of chlorination and partition coefficients

The octanol-water partition coefficient ${K_{{\rm{ow}}}}$ , as well as its logarithm ${\rm{log}}\;\left( {{K_{{\rm{ow}}}}} \right)$ , are widely used measures of the lipophilicity of a compound. ${K_{{\rm{ow}}}}$ correlates with the in vivo thermodynamic equilibrium partition (or preference) of the congeners among water-rich and fat-rich tissues. In turn, the number of chlorine atoms in PCDD/F and PCB congeners correlates with ${K_{{\rm{ow}}}}$ . With more chlorine atoms in a congener, its ${K_{{\rm{ow}}}}$ increases, and its solubility preference for the aqueous phase decreases(Reference Shiu, Doucette and Gobas168). All PCDD/Fs and PCBs have a ${K_{{\rm{ow}}}}$ larger than 1, which implies that these substances are lipophilic and hydrophobic, hence their distribution into body and milk fat. ${K_{{\rm{ow}}}}\;$ values are used to predict toxicokinetic model parameters, since they correlate with the preference of contaminants for different compartments, such as blood, tissues and udder fat in the cow. The TRs(Reference Welsch-Pausch and McLachlan64,Reference Lorenzi, Angelone and Ferretti86,Reference McLachlan, Thoma and Reissinger109,Reference McLachlan and Richter113,Reference Costera, Feidt and Marchand144,Reference McLachlan169) , BTFs(Reference Fries, Paustenbach and Mather72) and bioavailability(Reference Olling, Derks and Berende102) of 2,3,7,8-substituted PCDD/Fs, with the exception of 2,3,7,8-TCDF and 1,2,3,7,8-PCDF, decrease with the number of chlorine substituents. This behaviour can be attributed to a decrease in absorption over the gastrointestinal wall(Reference Tremolada, Guazzoni and Parolini108,Reference Jilg, Mueller and Hagenmaier112) , due to the increased hydrophobicity of highly chlorinated PCDD/Fs(Reference Schuler, Schmid and Schlatter63).

Furthermore, lower chlorinated PCDD/Fs migrate principally into the lipoprotein fraction of blood, while higher chlorinated PCDD/Fs bind more readily to other blood plasma proteins with a prolonged dwell time in systemic circulation(Reference Hirako127,Reference Patterson, Fürst and Henderson170) . In addition, the extent of deposition of PCDD/Fs into subcutaneous fat was found to be inversely correlated to the degree of chlorination. This could explain the faster elimination rates of heptachlorinated PCDD/Fs compared with lower-chlorinated congeners(Reference Olling, Derks and Berende102). A faster distribution of lower-chlorinated PCDD/Fs from liver to fat tissues was also observed by Richter and McLachlan (2001), who found that higher-chlorinated PCDD/Fs were accumulated in the liver(Reference Richter and McLachlan171), likely owing to the low affinity of these compounds to lipoproteins(Reference Hirako127,Reference McLachlan169) . Other studies found no correlations between elimination half-lives and the number of chlorine substituents(Reference Tuinstra, Roos and Berende116,Reference Huwe and Smith135) .

Fluxes of PCBs in faeces of cows showed that absorption in the gastrointestinal tract also decreased with increasing degree of chlorination(Reference McLachlan111). Furthermore, the residence times of PCBs in the cow, calculated as body burden divided by amount absorbed, increased with the degree of chlorination(Reference Thomas, Sweetman and Jones65). TRs derived in studies with lactating alpine goats also showed no clear correlation(Reference Ounnas, Feidt and Toussaint146) or indications for an inverse correlation with degree of chlorination(Reference Costera, Feidt and Marchand144). However, aside from the degree of chlorination and ${K_{{\rm{ow}}}}$ , the substitution pattern, and thus metabolic stability, also needs to be considered when discussing effects on transfer parameters.

Metabolism of PCDD/Fs and PCBs

PCDD/Fs and PCBs are known to be persistent in the environment, with half-lives of up to several decades(Reference Sinkkonen and Paasivirta12). However, some congeners are more readily enzymatically metabolised by the cow, with the velocity and extent being largely dependent on the chlorination pattern. One metabolic pathway is the typical phase-1 metabolism, that is, hydroxylation of PCBs in the animal and subsequent excretion via urine and, to a lesser extent, via milk(Reference Gardner, Righter and Roach75,Reference Safe, Platonow and Hutzinger120) .

However, metabolism of PCBs is congener specific. One hypothesis is that PCBs with a 4,4′-substitution pattern and, to a lesser degree, with a 2,3,5-substitution pattern, are less likely to be metabolised(Reference McLachlan111). Another hypothesis is that PCBs with adjacent meta- and para-hydrogen atoms are more chemically stable, and PCBs with adjacent ortho- and meta-hydrogen atoms are only partially metabolised(Reference Thomas, Sweetman and Jones65). However, both hypotheses have their shortcomings, as there are PCBs that do not behave as predicted by these rules in terms of substitution patterns(Reference Thomas, Sweetman and Jones65).

PCDD/F congeners with a 2,3,7,8-chlorine substitution pattern are largely metabolically stable and are transferred into the milk in unchanged form, while other congeners can be metabolised in the animal(Reference Welsch-Pausch and McLachlan64,Reference Firestone, Clower and Borsetti76,Reference McLachlan, Thoma and Reissinger109) . Additionally, 2,3,7,8-substituted furan congeners without chlorine at 4 and 6 position(Reference Fries, Paustenbach and Mather72,Reference Fries, Paustenbach and Luksemburg73) , such as 2,3,7,8-TCDF(Reference Olling, Derks and Berende102) and 1,2,3,7,8-PeCDF(Reference McLachlan and Richter113), are susceptible to metabolism. In vitro fermentation experiments showed no significant metabolic degradation of PCBs(Reference Willett, Loerch and Willett172) or PCDD/Fs(Reference Fries, Paustenbach and Luksemburg73) by rumen microorganisms.

The transfer parameters for dioxins and PCBs in experimental studies can be influenced by several factors, some attributable to the cow physiology and others related to the chemical contaminant. There are some contradictions in the literature regarding the assessment of the influence of the factors responsible for the variability in transfer parameters, often aggravated by incomplete information. Comparability is not trivial, since there are many different possible study designs, due to many variables, such as animal species and breed, milk yield, type of congener administered and duration of supplementation. On the basis of existing gaps in the data and the multitude of different studies, the literature data are not always applicable and comparable to produce predictive models. However, several toxicokinetic predictive models have been proposed and successfully parametrised in the literature, and shown to agree well at least with the specific experimental data they are based on.

Conclusions

The present review provides a comprehensive overview of the transfer of PCDD/Fs and PCBs from feed into the milk of dairy cows, and reveals several key aspects that may influence their transfer kinetics. We have identified data gaps for dairy cows regarding congener-specific BTFs, especially for PCBs, and many congener-specific elimination half-lives. Even when transfer parameters are available, there is a high degree of variability between studies, sometimes even among studies with similar experimental setups. Nevertheless, our compilation and graphical depiction of transfer parameters testify to the strong differences in kinetics for the various PCDD/F and PCB congeners. Thus, there is a necessity to separately quantify the transfer of as many congeners as possible to perform risk assessment and risk management.

It seems likely that a considerable amount of variability in transfer parameters (i.e. TR, TF, BTF and half-lives) for the PCDD/F and PCB congeners summarised here is due to the considerable heterogeneity in study design and data quality. However, a classification or weighting of influence factors is not possible, as most studies were not explicitly designed to evaluate the influence of, for example, lactational stage, breed, source or route of exposure on the transfer kinetics. The literature contains valuable indications on the dependency of transfer parameters on the aforementioned factors, but conclusions on such dependencies are sometimes based on a low sample size, so that the bias of individual animals may be confused for systematic variation. Future studies should try to address the variability by using an appropriate number of animals with clearly defined physiological conditions (e.g. breed, milk and milk fat yield, lactational stage, NEB or PEB phase, body weight/body fat weight and animal health status). Furthermore, the influence of metabolic status on the transfer of PCDD/Fs and PCBs into milk is of special interest and yet difficult to pinpoint. Milk yield and body fat percentage are intertwined with the NEB and PEB phases of lactation but are also important factors on their own. High variability between individuals further complicates analyses. To study the effect of the metabolic status on the transfer of PCDD/Fs and PCBs, it is advisable to expose a large enough group of animals postpartum during a NEB phase and the same group again in high lactation during the PEB phase, giving enough time for depuration between both exposure phases.

Over the past decades, milk yield has increased from less than 30 kg/d to over 40 kg/d during high lactation(Reference VandeHaar and St-Pierre140). The question arises of how existing data and models represent transfer under the ranges of the current milk yields. At the moment, existing data do not allow for a conclusive answer to this question, as experimental data for modern highly productive cows are scarce. Advanced breeding efforts will lead to further increases in milk yield, but also impact other areas of dairy cow genetics and metabolism. At the same time, ethical requirements will demand a reduction in animal numbers in forthcoming feeding trials. Thus, the question arises in how far (i) data from previous studies can be utilised for predicting the transfer of contaminants in future breeds and (ii) valid conclusions can be drawn from smaller animal studies. Both issues are closely related to the advancement of prediction tools and toxicokinetic modelling approaches. These topics are the subject of the part II of the present review article(Reference Moenning, Krause, Lamp, Maul, Schenkel, Fürst, Pieper and Numata56).

In this review, we have highlighted where the gaps in experimental data lie and hope to aid future researchers to design the experiments of the future (and the models they engender) to support PCDD/F and PCB risk assessment and risk management of milk products.

Acknowledgements

We would like to thank Florian Ober for R scripts and for the preparation of the whisker boxplots as well as Ashley Gedon for proofreading.

Financial Support

H.S. and P.F. acknowledge the support from the Maria Sibylla Merian Fellowship; the Fellowship had no role in the design, analysis or writing of this article. All authors received no specific grant from any funding agency, commercial or not-for-profit sectors.

Conflict of interest

None.

Authorship

Torsten Krause: literature search, data compilation and analysis, writing (original draft), figures

Jan-Louis Moenning: literature search, modelling aspects, writing (review and editing), figures

Julika Lamp: literature search, data gathering, introduction, animal health and physiology aspects, figures

Ronald Maul: conceptualisation, writing (review and editing)

Hans Schenkel: conceptualisation, literature search

Peter Fürst: conceptualisation, literature search

Robert Pieper: animal health and physiology aspects, conceptualisation, writing (review and editing)

Jorge Numata: supervision, conceptualisation, modelling aspects, writing (original draft, review and editing), figures

Supplementary material

To view supplementary material for this article, please visit https://doi.org/10.1017/S0954422422000178.

References

Safe, S (1990) Polychlorinated biphenyls (PCBs), dibenzo-p-dioxins (PCDDs), dibenzofurans (PCDFs), and related compounds: environmental and mechanistic considerations which support the development of toxic equivalency factors (TEFs). Crit Rev Toxicol 21, 5188.CrossRefGoogle ScholarPubMed
Kulkarni, PS, Crespo, JG & Afonso, CAM (2008) Dioxins sources and current remediation technologies – a review. Environ Int 34, 139153.CrossRefGoogle ScholarPubMed
Fiedler, H (1996) Sources of PCDD/PCDF and impact on the environment. Chemosphere 32, 5564.CrossRefGoogle ScholarPubMed
Breivik, K, Sweetman, A, Pacyna, JM, et al. (2002) Towards a global historical emission inventory for selected PCB congeners – a mass balance approach: 2. Emissions. Sci Total Environ 290, 199224.CrossRefGoogle ScholarPubMed
Breivik, K, Sweetman, A, Pacyna, JM, et al. (2007) Towards a global historical emission inventory for selected PCB congeners – a mass balance approach: 3. An update. Sci Total Environ 377, 296307.CrossRefGoogle ScholarPubMed
Jahnke, JC & Hornbuckle, KC (2019) PCB emissions from paint colorants. Environ Sci Technol 53, 51875194.CrossRefGoogle ScholarPubMed
Herkert, NJ, Jahnke, JC & Hornbuckle, KC (2018) Emissions of tetrachlorobiphenyls (PCBs 47, 51, and 68) from polymer resin on kitchen cabinets as a non-Aroclor source to residential air. Environ Sci Technol 52, 51545160.CrossRefGoogle ScholarPubMed
Huang, H & Buekens, A (1995) On the mechanisms of dioxin formation in combustion processes. Chemosphere 31, 40994117.CrossRefGoogle Scholar
Ballschmiter, K & Zell, M (1980) Analysis of polychlorinated-biphenyls (PCB) by glass-capillary gas-chromatography – composition of technical Aroclor-PCB and Clophen-PCB mixtures. Fresen Z Anal Chem 302, 2031.CrossRefGoogle Scholar
Weber, R, Herold, C, Hollert, H, et al. (2018) Life cycle of PCBs and contamination of the environment and of food products from animal origin. Environ Sci Pollut Res 25, 1632516343.CrossRefGoogle ScholarPubMed
European Food Safety Authority (2012) Update of the monitoring of levels of dioxins and PCBs in food and feed. EFSA J 10, 2832.Google Scholar
Sinkkonen, S & Paasivirta, J (2000) Degradation half-life times of PCDDs, PCDFs and PCBs for environmental fate modeling. Chemosphere 40, 943949.CrossRefGoogle ScholarPubMed
Hoogenboom, R, Traag, W, Fernandes, A, et al. (2015) European developments following incidents with dioxins and PCBs in the food and feed chain. Food Control 50, 670683.CrossRefGoogle Scholar
Hoogenboom, RLAP, Malisch, R, van Leeuwen, SPJ, et al. (2020) 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 746, 141098.CrossRefGoogle ScholarPubMed
MacLachlan, DJ & Bhula, R (2008) Estimating the residue transfer of pesticides in animal feedstuffs to livestock tissues, milk and eggs: a review. Aust J Exp Agric 48, 589598.CrossRefGoogle Scholar
Knutsen, HK, Alexander, J, Barregård, L, et al. (2018) Risk for animal and human health related to the presence of dioxins and dioxin-like PCBs in feed and food. EFSA J 16, e05333.Google Scholar
FAO/WHO Codex Alimentarius (2013) Guidelines on the Application of Risk Assessment for Feed. http://www.fao.org/input/download/standards/13322/CXG_080e.pdf (accessed 11 March 2021).Google Scholar
Bruner-Tran, KL, Duleba, AJ, Taylor, HS, et al. (2016) Developmental toxicant exposure is associated with transgenerational adenomyosis in a murine model. Biol Reprod 95, 110.CrossRefGoogle ScholarPubMed
Schecter, A & Gasiewicz, TA (2003) Dioxins and Health. Hoboken: John Wiley & Sons, Inc. CrossRefGoogle Scholar
Mínguez-Alarcón, L, Sergeyev, O, Burns, JS, et al. (2017) A longitudinal study of peripubertal serum organochlorine concentrations and semen parameters in young men: the Russian children’s study. Environ Health Perspect 125, 460466.CrossRefGoogle ScholarPubMed
Okey, AB (2007) An aryl hydrocarbon receptor odyssey to the shores of toxicology: the Deichmann lecture, International Congress of Toxicology-XI. Toxicol Sci 98, 538.CrossRefGoogle Scholar
Abel, J & Haarmann-Stemmann, T (2010) An introduction to the molecular basics of aryl hydrocarbon receptor biology. Biol Chem 391, 12351248.CrossRefGoogle Scholar
Birnbaum, LS (1994) The mechanism of dioxin toxicity: relationship to risk assessment. Environ Health Perspect 102, Suppl 9, 157167.CrossRefGoogle ScholarPubMed
Van den Berg, M, Birnbaum, L, Bosveld, AT, et al. (1998) Toxic equivalency factors (TEFs) for PCBs, PCDDs, PCDFs for humans and wildlife. Environ Health Perspect 106, 775792.CrossRefGoogle ScholarPubMed
Kutz, FW, Barnes, DG, Bottimore, DP, et al. (1990) The international toxicity equivalency factor (I-TEF) method of risk assessment for complex mixtures of dioxins and related compounds. Chemosphere 20, 751757.CrossRefGoogle Scholar
Van den Berg, M, Birnbaum, LS, Denison, M, et al. (2006) The 2005 World Health Organization reevaluation of human and mammalian toxic equivalency factors for dioxins and dioxin-like compounds. Toxicol Sci 93, 223241.CrossRefGoogle ScholarPubMed
International Agency for Research on Cancer (2012) Chemical Agents and Related Occupations – IARC Monographs on the Evaluation of Carcinogenic Risks to Humans. vol. 100F. Lyon: World Health Organization, International Agency for Research on Cancer.Google Scholar
Lallas, PL (2001) The Stockholm convention on persistent organic pollutants. Am J Int Law 95, 692708.CrossRefGoogle Scholar
EFSA Panel on Contaminants in the Food Chain (2005) Opinion of the Scientific Panel on contaminants in the food chain [CONTAM] related to the presence of non dioxin-like polychlorinated biphenyls (PCB) in feed and food. EFSA J 3, 284.CrossRefGoogle Scholar
World Health Organization (2016) Safety Evaluation of Certain Food Additives and Contaminants, Supplement 1: Non-dioxin-like Polychlorinated Biphenyls/Prepared by the Eightieth Meeting of the Joint FAO/WHO Expert Committee on Food Additives (JECFA): World Health Organization, Geneva.Google Scholar
Malisch, R (2017) Incidents with dioxins and PCBs in food and feed-Investigative work, risk management and economic consequences. J Environ Protect 8, 744785.CrossRefGoogle Scholar
Sanger, VL, Scott, L, Hamdy, A, et al. (1958) Alimentary toxemia in chickens. J Am Vet Med Assoc 133, 172176.Google ScholarPubMed
Higginbotham, GR, Huang, A, Firestone, D, et al. (1968) Chemical and toxicological evaluations of isolated and synthetic chloro derivatives of dibenzo-p-dioxin. Nature 220, 702703.CrossRefGoogle ScholarPubMed
Bernard, A, Hermans, C, Broeckaert, F, et al. (1999) Food contamination by PCBs and dioxins. Nature 401, 231232.CrossRefGoogle ScholarPubMed
Malisch, R (2000) Increase of the PCDD/F-contamination of milk, butter and meat samples by use of contaminated citrus pulp. Chemosphere 40, 10411053.CrossRefGoogle ScholarPubMed
Hoogenboom, R, Bovee, T, Portier, L, et al. (2004) The German bakery waste incident; use of a combined approach of screening and confirmation for dioxins in feed and food. Talanta 63, 12491253.CrossRefGoogle ScholarPubMed
Llerena, JJ, Abad, E, Caixach, J, et al. (2003) An episode of dioxin contamination in feedingstuff: the choline chloride case. Chemosphere 53, 679683.CrossRefGoogle Scholar
Stachel, B, Christoph, EH, Gotz, R, et al. (2006) Contamination of the alluvial plain, feeding-stuffs and foodstuffs with polychlorinated dibenzo-p-dioxins, polychlorinated dibenzofurans (PCDD/Fs), dioxin-like polychlorinated biphenyls (DL-PCBs) and mercury from the River Elbe in the light of the flood event in August 2002. Sci Total Environ 364, 96112.CrossRefGoogle ScholarPubMed
Hoogenboom, LAP, Eijkeren, JCHV, Zeilmaker, MJ, et al. (2007) A novel source for dioxins present in recycled fat from gelatin production. Chemosphere 68, 814823.CrossRefGoogle ScholarPubMed
Tlustos, C, Sheridan, M, O’Sullivan, D, et al. (2012) The dioxin contamination incident in Ireland, 2008: analytical results and congener patterns. Food Addit Contam: Part A 29, 128138.CrossRefGoogle ScholarPubMed
Brambilla, G, Fochi, I, Falce, M, et al. (2008) PCDD and PCDF depletion in milk from dairy cows according to the herd metabolic scenario. Chemosphere 73, S216S219.CrossRefGoogle Scholar
Hoogenboom, R, Zeilmaker, M, van Eijkeren, J, et al. (2010) Kaolinic clay derived PCDD/Fs in the feed chain from a sorting process for potatoes. Chemosphere 78, 99105.CrossRefGoogle ScholarPubMed
Piskorska-Pliszczynska, J, Maszewski, S, Mikolajczyk, S, et al. (2017) Elimination of dioxins in milk by dairy cows after the long-term intake of contaminated sugar beet pellets. Food Addit Contam: Part A 34, 842852.Google ScholarPubMed
Fürst, P, Beck, H & Theelen, R (1992) Assessment of human intake of PCDDs and PCDFs from different environmental sources. Toxic Subst J 12, 133150.Google Scholar
European Commission (2011) Commission Regulation (EU) No 1259/2011 of 2 December 2011 amending Regulation (EC) No 1881/2006 as regards maximum levels for dioxins, dioxin-like PCBs and non dioxin-like PCBs in foodstuffs Text with EEA relevance.Google Scholar
European Commission (2012) Commission Regulation (EU) No 277/2012 of 28 March 2012 amending Annexes I and II to Directive 2002/32/EC of the European Parliament and of the Council as regards maximum levels and action thresholds for dioxins and polychlorinated biphenyls Text with EEA relevance.Google Scholar
European Commission (2011) 2011/516/EU: Commission Recommendation of 23 August 2011 on the reduction of the presence of dioxins, furans and PCBs in feed and food Text with EEA relevance.Google Scholar
European Commission (2021) On the Review and Update of the Second European Union Implementation Plan in Accordance with Article 8 (4) of Regulation No 850/2004 on Persistent Organic Pollutants. no. COM(2018) 848 final. Brussels.Google Scholar
Lascano Alcoser, VH, Velthuis, AG, Hoogenboom, LA, et al. (2011) Financial impact of a dioxin incident in the Dutch dairy chain. J Food Prot 74, 967979.CrossRefGoogle ScholarPubMed
Akkan, Z, Oltmanns, J & Kalberlah, F (2004) Literaturstudie zur Ermittlung des Depositionswertes von Dioxinen, Furanen und dioxinähnlichen PCB. Karlsruhe: Verlag-Auslieferung der LfU bei JVA.Google Scholar
Hoogenboom, LAP (2005) Behavior of Polyhalogenated and Polycyclic Aromatic Hydrocarbons in Food-Producing Animals. In: Reviews in Food and Nutrition Toxicity, pp. 269299 [VR Preedy and RR Watson, editors]. London: Taylor and Francis.Google Scholar
Rychen, G, Jurjanz, S, Toussaint, H, et al. (2008) Dairy ruminant exposure to persistent organic pollutants and excretion to milk. Animal 2, 312323.Google ScholarPubMed
Hennecke, D, Rolf, A & Becker, L (2011) Expositionsbetrachtung und Beurteilung des Transfers von Dioxinen, dioxinähnlichen PCB und PCB. vol. 57, Texte: Umweltbundesamt.Google Scholar
Rychen, G, Jurjanz, S, Fournier, A, et al. (2014) Exposure of ruminants to persistent organic pollutants and potential of decontamination. Environ Sci Pollut Res 21, 64406447.CrossRefGoogle ScholarPubMed
Amutova, F, Delannoy, M, Baubekova, A, et al. (2021) Transfer of persistent organic pollutants in food of animal origin – meta-analysis of published data. Chemosphere 262, 128351.CrossRefGoogle ScholarPubMed
Moenning, JL, Krause, T, Lamp, J, Maul, R, Schenkel, H, Fürst, P, Pieper, R & Numata, J (in press) Transfer of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs) and polychlorinated biphenyls (PCBs) from Oral exposure into cow’s milk – Part II: Toxicokinetic predictive models for risk assessment. Nutrition Research Reviews.Google Scholar
Tuinstra, LGMT, Vreman, K, Roos, AH, et al. (1981) Excretion of certain chlorobiphenyls into the milk fat after oral administration. Neth Milk Dairy J (Neth) 35, 147157.Google Scholar
Connett, P & Webster, T (1987) An estimation of the relative human exposure to 2,3,7,8-TCDD emissions via inhalation and ingestion of cow’s milk. Chemosphere 16, 20792084.Google Scholar
Farries, E (1990) PCB in Milch und Fleisch: eine Belastung mit Polychlorierten Biphenylen ist nur schwer wieder zu beseitigen. Tierzüchter 42, 539541.Google Scholar
Steimer, J, Blüthgen, A, Heeschen, W, et al. (1990) Untersuchungen zur Beeinflussung der Ausscheidung von Aflatoxin M1 durch polychlorierte Biphenyle beim laktierenden Rind. Kiel Milchwirtsch Forschungsber 42, 543552.Google Scholar
Roos, AH, Berende, PLM, Traag, WA, et al. (1991) Uitscheiding van polygechloreerde dibenzo-p-dioxinen en -furanen in melk van koeien, die tijdens de droogstand zijn gevoerd met dioxinen besmet voer: RIKILT - Institute of Food Safety, Wageningen, Netherlands.Google Scholar
Slob, W, Olling, M, Derks, HJGM, et al. (1995) Congener-specific bioavailability of PCDD/Fs and coplanar PCBs in cows: laboratory and field measurements. Chemosphere 31, 38273838.CrossRefGoogle ScholarPubMed
Schuler, F, Schmid, P & Schlatter, C (1997) Transfer of airborne polychlorinated dibenzo-p-dioxins and dibenzofurans into dairy milk. J Agric Food Chem 45, 41624167.CrossRefGoogle Scholar
Welsch-Pausch, K & McLachlan, MS (1998) Fate of airborne polychlorinated dibenzo-p-dioxins and dibenzofurans in an agricultural ecosystem. Environ Pollut 102, 129137.CrossRefGoogle Scholar
Thomas, GO, Sweetman, AJ & Jones, KC (1999) Metabolism and body-burden of PCBs in lactating dairy cows. Chemosphere 39, 15331544.CrossRefGoogle ScholarPubMed
Thomas, GO, Sweetman, AJ & Jones, KC (1999) Input–output balance of polychlorinated biphenyls in a long-term study of lactating dairy cows. Environ Sci Technol 33, 104112.Google Scholar
Lorber, M, Fries, G, Winters, D, et al. (2000) A study of the mass balance of dioxins and furans in lactating cows in background conditions. Part 2: mass balance and bioconcentration factors. Organohalogen Compd 46, 326329.Google Scholar
Winters, D, Fries, GF, Lorber, M, et al. (2000) A study of the mass balance of dioxins and furans in lactating cows in background conditions. Part 1: study design and analysis of feed. Organohalogen Compd 46, 534537.Google Scholar
Tato, L, Tremolada, P, Ballabio, C, et al. (2011) Seasonal and spatial variability of polychlorinated biphenyls (PCBs) in vegetation and cow milk from a high altitude pasture in the Italian Alps. Environ Pollut 159, 26562664.CrossRefGoogle ScholarPubMed
Fries, GF (1972) Polychlorinated biphenyl residues in milk of environmentally and experimentally contaminated cows. Environ Health Perspect 1, 5559.CrossRefGoogle ScholarPubMed
Fries, GF, Marrow, GS & Gordon, CH (1973) Long-term studies of residue retention and excretion by cows fed a polychlorinated biphenyl (Aroclor 1254). J Agric Food Chem 21, 117121.CrossRefGoogle ScholarPubMed
Fries, GF, Paustenbach, DJ, Mather, DB, et al. (1999) A congener specific evaluation of transfer of chlorinated dibenzo-p-dioxins and dibenzofurans to milk of cows following ingestion of pentachlorophenol-treated wood. Environ Sci Technol 33, 11651170.CrossRefGoogle Scholar
Fries, GF, Paustenbach, DJ & Luksemburg, WJ (2002) Complete mass balance of dietary polychlorinated dibenzo-p-dioxins and dibenzofurans in dairy cattle and characterization of the apparent synthesis of hepta- and octachlorodioxins. J Agric Food Chem 50, 42264231.CrossRefGoogle ScholarPubMed
Furr, AK, Mertens, DR, Gutenmann, WH, et al. (1974) Fate of polychlorinated biphenyls, metals, and other elements in papers fed to lactating cows. J Agric Food Chem 22, 954959.CrossRefGoogle ScholarPubMed
Gardner, AM, Righter, HF & Roach, JAG (1976) Excretion of hydroxylated polychlorinated biphenyl metabolites in cows’ milk. J Assoc Off Anal Chem 59, 273277.Google ScholarPubMed
Firestone, D, Clower, M, Borsetti, AP, et al. (1979) Polychlorodibenzo-p-dioxin and pentachlorophenol residues in milk and blood of cows fed technical pentachlorophenol. J Agric Food Chem 27, 11711177.Google ScholarPubMed
Jensen, DJ & Hummel, RA (1982) Secretion of TCDD in milk and cream following the feeding of TCDD to lactating dairy cows. Bull Environ Contam Toxicol 29, 440446.CrossRefGoogle ScholarPubMed
Willett, LB & Liu, TTY (1982) Effects of thyroprotein on excretion of polychlorinated biphenyls by lactating cows. J Dairy Sci 65, 7280.CrossRefGoogle ScholarPubMed
Jones, D, Safe, S, Morcom, E, et al. (1987) Bioavailability of tritiated 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) administered to Holstein dairy cows. Chemosphere 16, 17431748.CrossRefGoogle Scholar
Jones, D, Safe, S, Morcom, E, et al. (1989) Bioavailability of grain and soil-borne tritiated 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) administered to lactating Holstein cows. Chemosphere 18, 12571263.CrossRefGoogle Scholar
Willett, LB, Liu, TTY, Durst, HI, et al. (1987) Health and productivity of dairy cows fed polychlorinated biphenyls. Fundam Appl Toxicol 9, 6068.Google ScholarPubMed
Huwe, JK & Smith, DJ (2003) Uptake and excretion of PCDD/Fs from a magnesium mineral supplement fed to dairy cows. Organohalogen Compd 64, 419422.Google Scholar
Hoogenboom, RLAP, Klop, A, Herbes, R, et al. (2015) Carry-over of polychlorinated dibenzo-p-dioxins and dibenzofurans (PCDD/Fs) and polychlorinated biphenyls (PCBs) in dairy cows fed smoke contaminated maize silage or sugar beet pulp. Chemosphere 137, 214220.CrossRefGoogle ScholarPubMed
Bertocchi, L, Ghidini, S, Fedrizzi, G, et al. (2015) Case-study and risk management of dioxins and PCBs bovine milk contaminations in a high industrialized area in Northern Italy. Environ Sci Pollut Res 22, 111.CrossRefGoogle Scholar
Thomas, GO, Sweetman, AJ, Lohmann, R, et al. (1998) Derivation and field testing of air-milk and feed-milk transfer factors for PCBs. Environ Sci Technol 32, 35223528.Google Scholar
Lorenzi, V, Angelone, B, Ferretti, E, et al. (2020) PCDD/Fs, dl-PCBs, and ndl-PCBs in dairy cows: carryover in milk from a controlled feeding study. J Agric Food Chem 68, 22012213.Google ScholarPubMed
Heeschen, W & Blüthgen, A (1986) Aktuelle Probleme und Ergebnisse der Carry over-Forschung: Hexachlorcyclohexan (HCH) und polychlorierte Biphenyle (PCB). Probleme des Übergangs unerwünschter Stoffe in und auf Futtermitteln in die Nahrungskette: Informationsveranstaltung der Arbeitsgruppe “Carry over” am 7. und 8. Oktober 1986 in der Bundesforschungsanstalt für Landwirtschaft Braunschweig-Völkenrode (FAL), 193211.Google Scholar
Ewers, C (1987) Untersuchungen über Aufnahme, Verbleib und Ausscheidung (Carry over) Polychlorierter Biphenyle (PCB) beim laktierenden Rind. Kiel: Agrarwissenschaftliche Fakultät der Christian-Albrechts-Universität zu Kiel.Google Scholar
Ewers, C & Heeschen, W (1987) Zum Carry over polychlorierter Biphenyle aus Futtermitteln in die Milch. Milk Sci Int – Milchwissenschaft 42, 757761.Google Scholar
Ewers, C, Reichmuth, J, Wetzel, S, et al. (1989) Untersuchungen zum Verbleib lipophiler persistenter Organochlorverbindungen im Körper laktierender Rinder nach oraler Aufnahme am Beispiel des PCB-Kongeners Nr. 153. Kiel Milchwirtsch Forschungsber 41, 7595.Google Scholar
Ruoff, U (1995) Untersuchungen zum Übergang ausgewählter polychlorierter Dibenzo-para-dioxine und -furane nach oraler Supplementierung in die Milch laktierender Kühe. Kiel: Agrarwissenschaftliche Fakultät der Christian-Albrechts-Universität zu Kiel.Google Scholar
Blüthgen, A (2006) Übergang der Indikator-PCB Kongenere nach oraler Zufuhr in das Milchfett laktierender Kühe: Jahresbericht der Bundesforschungsanstalt für Ernährung und Lebensmittel. Karlsruhe: Bundesforschungsanstalt für Ernährung und Lebensmittel.Google Scholar
Schulz, AJ (2005) Auswirkungen originär Dioxin-belasteten Grundfutters auf die Dioxingehalte in Milch und Schlachtkörpern von Rindern und Schafen. Hannover: Institut für Tierernährung der Tierärztlichen Hochschule Hannover.Google Scholar
Schulz, AJ, Wiesmüller, T, Appuhn, H, et al. (2005) Dioxin concentration in milk and tissues of cows and sheep related to feed and soil contamination. J Anim Physiol Anim Nutr 89, 7278.CrossRefGoogle ScholarPubMed
Platonow, NS, Funnell, HS, Bullock, DH, et al. (1971) Fate of polychlorinated biphenyls in dairy products processed from the milk of exposed cows. J Dairy Sci 54, 13051308.CrossRefGoogle ScholarPubMed
Saschenbrecker, PW, Funnell, HS & Platonow, NS (1972) Persistence of polychlorinated biphenyls in the milk of exposed cows. Vet Rec 90, 100102.CrossRefGoogle ScholarPubMed
Platonow, NS & Meads, EB (1975) Distribution and excretion of two chlorinated biphenyl isomers: 4-chlorobiphenyl and decachlorobiphenyl in lactating bovine. Can J Comp Med 39, 104.Google ScholarPubMed
Arnott, DR, Bullock, DH & Platonow, NS (1977) Kinetics of polychlorinated biphenyls (Aroclor 1254) in lactating bovines and their distribution in dairy products. J Food Protect 40, 296299.CrossRefGoogle ScholarPubMed
Klein, U (1991) Verlauf der Konzentration einiger polychlorierter Biphenyle in Milch, Blut und Geweben belasteter Milchkühe bei Variation von Leistung und Fütterung. Hannover: Tierärztliche Hochschule Hannover.Google Scholar
Klein, U, Drochner, W, Forschner, E, et al. (1992) PCB-Konzentrationen in Haaren, Blut, Geweben und Ausscheidungen chronisch belasteter Milchkühe, Mastbullen und Kälber. 1. Mitteilung: Verlauf und Beziehung der Konzentration einiger hochchlorierter Biphenyle in Haaren, Milch und Gewebenbei belasteten Ammenkühen und konventionell gehaltenen Milchrindern. Dtsch Tierärztl Wochenschr 99, 242248.Google Scholar
Olling, M, Derks, HJGM, Berende, PLM, et al. (1990) De toxicokinetiek van polychloor-dibenzo-pdloxlnen en -furanen in de lacterende koe (I). Bilthoven: Rijksinstituut voor Volksgezondheid en Milieuhygiene.Google Scholar
Olling, M, Derks, HJGM, Berende, PLM, et al. (1991) Toxicokinetics of eight 13C-labelled polychlorinated dibenzo-p-dioxins and – furans in lactating cows. Chemosphere 23, 13771385.CrossRefGoogle Scholar
Olling, M, Berende, PLM, Derks, HJGM, et al. (1992) Biobeschikbaarheid dioxinen en furanen uit vliegstof bij de koe. Rijksinstituut voor Volksgezondheid en Milieuhygiene, Bilthoven.Google Scholar
Derks, HJGM, Berende, PLM, Olling, M, et al. (1991) De toxicokinetiek van PCDD’s en PCDF’s in de lacterende koe (II). no. 328904002. Bilthoven: Rijksinstituut voor Volksgezondheid en Milieuhygiene.Google Scholar
Traag, W, Mengelers, MJB, Kan, CA, et al. (1999) Studies on the uptake and carry over of polychlorinated dibenzodioxins and dibenzofurans from contaminated citrus pulp pellets to cows milk. Organohalogen Compd 201, 201203.Google Scholar
Kerst, M, Waller, U, Reifenhäuser, W, et al. (2004) Carry-over rates of dioxin-like PCB from grass to cows’ milk. Organohalogen Compd 66, 24402444.Google Scholar
Körner, W, Kerst, M, Waller, U, et al. (2007) Untersuchung und Bewertung von Proben aus verschiedenen Umweltkompartimenten auf PCDD/PCDF sowie PCB unter Berücksichtigung der neuen WHO-Toxizitätsäquivalenzfaktoren. Augsburg: Bayerisches Landesamt für Umwelt.Google Scholar
Tremolada, P, Guazzoni, N, Parolini, M, et al. (2014) Predicting PCB concentrations in cow milk: validation of a fugacity model in high-mountain pasture conditions. Sci Total Environ 487, 471480.CrossRefGoogle ScholarPubMed
McLachlan, MS, Thoma, H, Reissinger, M, et al. (1990) PCDD/F in an agricultural food chain Part 1: PCDD/F mass balance of a lactating cow. Chemosphere 20, 10131020.CrossRefGoogle Scholar
McLachlan, MS (1992) Das Verhalten hydrophober chlororganischer Verbindungen in laktierenden Rindern, Fakultät für Biologie. Bayreuth: Chemie und Geowissenschaften der Universität Bayreuth.Google Scholar
McLachlan, MS (1993) Mass balance of polychlorinated biphenyls and other organochlorine compounds in a lactating cow. J Agric Food Chem 41, 474480.Google Scholar
Jilg, T, Mueller, W, Hagenmaier, H, et al. (1992) Carry-over of polychlorinated dibenzodioxins (PCDD) and dibenzofurans (PCDF) from feed into the milk of the dairy cow. Agribiol Res 45, 303310.Google Scholar
McLachlan, M & Richter, W (1998) Uptake and transfer of PCDD/Fs by cattle fed naturally contaminated feedstuffs and feed contaminated as a result of sewage sludge application. 1. Lactating cows. J Agric Food Chem 46, 11661172.Google Scholar
Driesen, C, Lerch, S, Siegenthaler, R, et al. (2022) Accumulation and decontamination kinetics of PCBs and PCDD/Fs from grass silage and soil in a transgenerational cow-calf setting. Chemosphere 296, 133951.CrossRefGoogle Scholar
Perry, TW, Everson, RJ, Hendrix, KS, et al. (1984) In utero exposure of bovine fetuses to polychlorinated biphenyls. J Dairy Sci 67, 224228.CrossRefGoogle ScholarPubMed
Tuinstra, LGMT, Roos, AH, Berende, PLM, et al. (1992) Excretion of polychlorinated dibenzo-p-dioxins and -furans in milk of cows fed on dioxins in the dry period. J Agric Food Chem 40, 17721776.CrossRefGoogle Scholar
Gill, IJ, Roberts, GS & Galvin, JW (1992) Management of land and livestock contaminated with polychlorinated biphenyls. Aust Vet J 69, 155158.CrossRefGoogle ScholarPubMed
Ungemach, LC (2013) Die Dioxinbelastung von Schlachtrindern aus einer Färsenvornutzung auf exponierten Grünland nach einer Ausmast mit unbelastetem Futter. Hannover: Tierärztliche Hochschule Hannover.Google Scholar
Perry, TW, Everson, RJ, Hendrix, KS, et al. (1981) Dietary Aroclor 1254 in the milk fat of lactating beef cattle. J Dairy Sci 64, 22622265.CrossRefGoogle ScholarPubMed
Safe, S, Platonow, N & Hutzinger, O (1975) Metabolism of chlorobiphenyls in the goat and cow. J Agric Food Chem 23, 259261.CrossRefGoogle ScholarPubMed
Willett, LB, Liu, TTY & Fries, GF (1990) Reevaluation of polychlorinated biphenyl concentrations in milk and body fat of lactating cows. J Dairy Sci 73, 21362142.CrossRefGoogle ScholarPubMed
Fürst, P, Krause, GHM, Hein, D, et al. (1993) PCDD/PCDF in cow’s milk in relation to their levels in grass and soil. Chemosphere 27, 13491357.Google Scholar
Rossi, F, Bertuzzi, T, Vitali, A, et al. (2010) Monitoring of the declining trend of Polychlorobifenyls concentration in milk of contaminated dairy cows. Ital J Anim Sci 9, 8892.CrossRefGoogle Scholar
Sips, AJAM, van Eijkeren, JCH, Bode, W, et al. (1999) Toxicokinetics of lindane (γ-HCH) and Several PCB Congeners in the Goat after Oral Dosing. Bilthoven: Rijksinstituut voor Volksgezondheid en Milieu (National Institute of Public Health and the Environment).Google Scholar
Vrecl, M, Ursic, M, Pogacnik, A, et al. (2005) Excretion pattern of co-planar and non-planar tetra- and hexa-chlorobiphenyls in ovine milk and faeces. Toxicol Appl Pharmacol 204, 170174.CrossRefGoogle ScholarPubMed
Ruoff, U & Fouzy, SM (2006) Distribution of PCDDs/PCDFs into milk and organs of egyptian Baladi goats after oral supplementation of dioxins. Kieler milchwirtschaftliche Forschungsberichte 58, 516.Google Scholar
Hirako, M (2008) Transfer and accumulation of persistent organochlorine compounds from bovine dams to newborn and suckling calves. J Agric Food Chem 56, 67686774.CrossRefGoogle ScholarPubMed
Fries, G (1996) Transport of chlorinated dioxin and furan contaminants in pentachlorophenol treated wood to milk and adipose tissue of dairy cattle. Organohalogen Compd 29, 447451.Google Scholar
Fries, GF (1987) Assessment of potential residues in foods derived from animals exposed to TCDD-contaminated soil. Chemosphere 16, 21232128.CrossRefGoogle Scholar
Rosenbaum, RK, McKone, TE & Jolliet, O (2009) CKow: a dynamic model for chemical transfer to meat and milk. Environ Sci Technol 43, 81918198.CrossRefGoogle ScholarPubMed
Ewers, C, Heeschen, W, Reichmuth, J, et al. (1987) Untersuchungen zum Carry over persistenter polychlorierter Biphenyle (PCB) aus Futtermitteln in die Milch. Kieler milchwirtschaftliche Forschungsberichte 39, 123138.Google Scholar
Heeschen, W, Ewers, C & Blüthgen, A (1993) Zum Carry over von Polychlorierten Biphenylen (PCB). Schriftenreihe des Bundesministers für Ernährung, Landwirtschaft und Forsten. Münster: Landeswirtschaftsverlag GmbH, pp. 123132.Google Scholar
Blüthgen, A, Heeschen, W & Ruoff, U (1994) Zum Carry-over toxikologisch relevanter Polychlordibenzodioxin- und -furankogenere in die Milch laktierender Kühe. Kieler milchwirtschaftliche Forschungsberichte 46, 139150.Google Scholar
Blüthgen, A, Heeschen, W & Ruoff, U (1995) Untersuchungen zum Übergang ausgewählter polychlorierter Dibenzo-para-dioxine und -furane (PCDD/F) nach oraler Supplementierung in die Milch laktierender Kühe. Milk Science International – Milchwissenschaft 50, 508513.Google Scholar
Huwe, JK & Smith, DJ (2005) Laboratory and on-farm studies on the bioaccumulation and elimination of dioxins from a contaminated mineral supplement fed to dairy cows. J Agric Food Chem 53, 23622370.CrossRefGoogle ScholarPubMed
Travis, CC & Hattemer-Frey, HA (1987) Human exposure to 2,3,7,8-TCDD. Chemosphere 16, 23312342.CrossRefGoogle Scholar
Stevens, JB & Gerbec, EN (1988) Dioxin in the agricultural food chain. Risk Anal 8, 329335.CrossRefGoogle ScholarPubMed
McKone, TE & Ryan, PB (1989) Human exposures to chemicals through food chains: an uncertainty analysis. Environ Sci Technol 23, 11541163.CrossRefGoogle Scholar
Webster, T & Connett, P (1990) The use of bioconcentration factors in estimating the 2,3,7,8-TCDD content of cow’s milk. Chemosphere 20, 779786.CrossRefGoogle Scholar
VandeHaar, MJ & St-Pierre, N (2006) Major advances in nutrition: relevance to the sustainability of the dairy industry. J Dairy Sci 89, 12801291.CrossRefGoogle Scholar
MacLachlan, DJ (2009) Influence of physiological status on residues of lipophilic xenobiotics in livestock. Food Addit Contam: Part A 26, 692712.CrossRefGoogle ScholarPubMed
Arstila, AU, Reggiani, G, Sorvari, TE, et al. (1981) Elimination of 2,3,7,8-tetrachlorodibenzo-p-dioxine in goat milk. Toxicol Lett 9, 215219.CrossRefGoogle ScholarPubMed
Grova, N, Feidt, C, Laurent, C, et al. (2002) [14C] Milk, urine and faeces excretion kinetics in lactating goats after an oral administration of [14C]polycyclic aromatic hydrocarbons. Int Dairy J 12, 10251031.CrossRefGoogle Scholar
Costera, A, Feidt, C, Marchand, P, et al. (2006) PCDD/F and PCB transfer to milk in goats exposed to a long-term intake of contaminated hay. Chemosphere 64, 650657.CrossRefGoogle ScholarPubMed
Costera Pastor, A (2007) Transfert des polluants organiques persistants (POP) du fourrage vers le lait chez le ruminant, Sciences agricoles. Nancy: Institut National Polytechnique de Lorraine.Google Scholar
Ounnas, F, Feidt, C, Toussaint, H, et al. (2010) Polychlorinated biphenyl and low polybrominated diphenyl ether transfer to milk in lactating goats chronically exposed to contaminated soil. Environ Sci Technol 44, 26822688.CrossRefGoogle ScholarPubMed
Feidt, C, Ounnas, F, Julien-David, D, et al. (2013) Relative bioavailability of soil-bound polychlorinated biphenyls in lactating goats. J Dairy Sci 96, 39163923.CrossRefGoogle ScholarPubMed
Fournier, A, Rychen, G, Marchand, P, et al. (2013) Polychlorinated biphenyl (PCB) decontamination kinetics in lactating goats (Capra hircus) following a contaminated corn silage exposure. J Agric Food Chem 61, 71567164.CrossRefGoogle ScholarPubMed
Diletti, G, Ceci, R, Ippoliti, C, et al. (2014) PCDD/F and dl-PCB transfer to milk in buffaloes exposed to contaminated feed. Organohalogen Compd 76, 15571560.Google Scholar
Park, YW & Haenlein, GFW (2006) Overview of milk of non-Bovine mammals. In Handbook of Milk of Non-Bovine Mammals, pp. 39.CrossRefGoogle Scholar
Willett, LB (1975) The excretory behavior of polychlorinated biphenyls in lactating cows fed normal and thyroprotein containing rations. J Dairy Sci 58, 765.Google Scholar
Rohatgi, A (2020) Webplotdigitizer: Version 4.4.Google Scholar
Derks, HJGM, Berende, PLM, Olling, M, et al. (1994) Pharmacokinetic modeling of polychlorinated dibenzo-p-dioxins (PCDDs) and furans (PCDFs) in cows. Chemosphere 28, 711715.CrossRefGoogle Scholar
Numata, J, Schafft, H, Spolders, M, et al. (2016) Toxikokinetische Modellierung als Erweiterung des klassischen Transfer-Versuchs. 128 VDLUFA-Kongress, Kongressband 2016 Rostock: Vorträge zum Generalthema: Anforderungen an die Verwertung von Reststoffen in der Landwirtschaft 73, 387392.Google Scholar
California Office of Environmental Health Hazard Assessment (OEHHA) (2012) Technical Support Document for Exposure Assessment and Stochastic Analysis. Appendix K – meat, milk, and egg transfer coefficients. K-1–K-32   .Google Scholar
Hendriks, AJ, Smítková, H & Huijbregts, MAJ (2007) A new twist on an old regression: transfer of chemicals to beef and milk in human and ecological risk assessment. Chemosphere 70, 4656.CrossRefGoogle Scholar
Bonate, PL (2011) Pharmacokinetic-Pharmacodynamic Modeling and Simulation, 2 ed.. Boston, MA: Springer.CrossRefGoogle Scholar
Stoop, WM, Bovenhuis, H, Heck, JML, et al. (2009) Effect of lactation stage and energy status on milk fat composition of Holstein-Friesian cows. J Dairy Sci 92, 14691478.CrossRefGoogle ScholarPubMed
McLachlan, MS (1994) Model of the fate of hydrophobic contaminants in cows. Environ Sci Technol 28, 24072414.Google ScholarPubMed
FAOSTAT Statistical Database (2020) Food and Agriculture Organization of the United Nations. Rome. https://www.fao.org/faostat/en/#home Google Scholar
Spagnuolo, MS, Sarubbi, F, Rossetti, C, et al. (2011) Effect of dioxin exposure on several indices of blood redox status in lactating buffalo cows. J Dairy Res 78, 154159.CrossRefGoogle ScholarPubMed
Spagnuolo, MS, Cigliano, L, Nebbia, C, et al. (2012) Analysis of plasma indices of redox homeostasis in dairy cows reared in polluted areas of Piedmont (northern Italy). Sci Total Environ 433, 450455.CrossRefGoogle ScholarPubMed
Genualdo, V, Perucatti, A, Iannuzzi, A, et al. (2012) Chromosome fragility in river buffalo cows exposed to dioxins. J Appl Genet 53, 221226.CrossRefGoogle ScholarPubMed
Davies, FG, Desouza, C, Waghela, S, et al. (1985) Possible dioxin poisoning in cattle. Vet Rec 117, 207.CrossRefGoogle ScholarPubMed
Viluksela, M, Heikkinen, P, van der Ven, LT, et al. (2014) Toxicological profile of ultrapure 2,2′,3,4,4′,5,5′-heptachlorbiphenyl (PCB 180) in adult rats. PLoS One 9, e104639.CrossRefGoogle ScholarPubMed
Savvateeva, D, Numata, J, Pieper, R, et al. (2020) Physiologically based toxicokinetic models and in silico predicted partition coefficients to estimate tetrachlorodibenzo-p-dioxin transfer from feed into growing pigs. Arch Toxicol 94, 187196.CrossRefGoogle ScholarPubMed
Inui, H, Itoh, T, Yamamoto, K, et al. (2014) Mammalian cytochrome P450-dependent metabolism of polychlorinated dibenzo-p-dioxins and coplanar polychlorinated biphenyls. Int J Mol Sci 15, 1404414057.CrossRefGoogle ScholarPubMed
Shiu, WY, Doucette, W, Gobas, FA, et al. (1988) Physical–chemical properties of chlorinated dibenzo-p-dioxins. Environ Sci Technol 22, 651658.CrossRefGoogle Scholar
McLachlan, MS (1997) A simple model to predict accumulation of PCDD/Fs in an agricultural food chain. Chemosphere 34, 12631276.CrossRefGoogle Scholar
Patterson, DG, Fürst, P, Henderson, LO, et al. (1989) Partitioning of in vivo bound PCDDs/PCDFs among various compartments in whole blood. Chemosphere 19, 135142.CrossRefGoogle Scholar
Richter, W & McLachlan, MS (2001) Uptake and transfer of PCDD/Fs by cattle fed naturally contaminated feedstuffs and feed contaminated as a result of sewage sludge application. 2. Nonlactating cows. J Agric Food Chem 49, 58575865.CrossRefGoogle ScholarPubMed
Willett, KL, Loerch, SC & Willett, LB (1989) Effects of halogenated hydrocarbons on rumen microorganisms. J Vet Diagn Invest 1, 120123.CrossRefGoogle ScholarPubMed
Figure 0

Table 1. Cattle breeds identified in the various transfer studies

Figure 1

Table 2. Studies used in assessment of transfer rates (TR), transfer factors (TF) and biotransfer factors (BTF) from feed to cow milk in chronological order

Figure 2

Table 3. Studies used in assessment of α and β elimination half-lives in milk of ruminants in chronological order

Figure 3

Fig. 1. Boxplots of transfer rates (TRs) for cows. Scatter points represent available data points from selected literature (refer to Table 2). Boxes (N > 5) are defined as the interquartile range (IQR) between 25th percentile (Q1) and 75th percentile (Q3) of the data according to the standard method. The black line in the box represents the median. Whiskers include data within 1.5 times of IQR below Q1 and above Q3. Plot generated with Python 3.10 using the Seaborn, Matplotlib, Numpy and Pandas libraries.

Figure 4

Fig. 2. Boxplots of transfer factors (TFs) in 88% dry matter basis for the feed and milk fat basis for the milk of cows. Scatter points represent available data points from selected literature (refer to Table 2). Boxes (N > 5) are defined as the interquartile range (IQR) between 25th percentile (Q1) and 75th percentile (Q3) of the data according to the standard method. The black line in the box represents the median. Whiskers include data within 1.5 times of IQR below Q1 and above Q3. Plot generated with Python 3.10 using the Seaborn, Matplotlib, Numpy and Pandas libraries.

Figure 5

Fig. 3. Boxplots of biotransfer factors (BTFs) in whole milk basis for the milk of cows. Scatter points represent available data points from selected literature (refer to Table 2). Boxes (N > 5) are defined as the interquartile range (IQR) between 25th percentile (Q1) and 75th percentile (Q3) of the data according to the standard method. The black line in the box represents the median. Whiskers include data within 1.5 times of IQR below Q1 and above Q3. Plot generated with Python 3.10 using the Seaborn, Matplotlib, Numpy and Pandas libraries.

Figure 6

Fig. 4. Boxplots of α elimination half-lives. Scatter points represent available data points from selected literature (refer to Table 3) mostly for cows, with red dots indicating data from goats. Boxes (N > 5) are defined as the interquartile range (IQR) between 25th percentile (Q1) and 75th percentile (Q3) of the data according to the standard method. The black line in the box represents the median. Whiskers include data within 1.5 times of IQR below Q1 and above Q3. Plot generated with Python 3.10 using the Seaborn, Matplotlib, Numpy and Pandas libraries.

Figure 7

Fig. 5. Boxplots of β elimination half-lives. Scatter points represent available data points from selected literature (refer to Table 3), mostly for cows with red dots indicating data from goats. Boxes (N > 5) are defined as the interquartile range (IQR) between 25th percentile (Q1) and 75th percentile (Q3) of the data according to the standard method. The black line in the box represents the median. Whiskers include data within 1.5 times of IQR below Q1 and above Q3. Plot generated with Python 3.10 using the Seaborn, Matplotlib, Numpy and Pandas libraries.

Figure 8

Table 4. Overview of the lactational status defined by days in milk (DIM) of cows in selected studies at beginning to end of the exposure phase to discuss the influence on transfer parameters

Figure 9

Table 5. Milk and milk fat production of selected studies in chronological order

Figure 10

Table 6. Transfer studies grouped according to the contamination source

Supplementary material: File

Krause et al. supplementary material

Krause et al. supplementary material 1

Download Krause et al. supplementary material(File)
File 905.1 KB
Supplementary material: File

Krause et al. supplementary material

Krause et al. supplementary material 2

Download Krause et al. supplementary material(File)
File 52.6 KB