- Research article
- Open Access
- Open Peer Review
Comparative oesophageal cancer risk assessment of hot beverage consumption (coffee, mate and tea): the margin of exposure of PAH vs very hot temperatures
BMC Cancer volume 18, Article number: 236 (2018)
Consumption of very hot (> 65 °C) beverages is probably associated with increased risk of oesophageal cancer. First associations were reported for yerba mate and it was initially believed that high content of polycyclic aromatic hydrocarbons (PAH) might explain the risk. Later research on other beverage groups such as tea and coffee, which are also consumed very hot, found associations with increased risk of oesophageal cancer as well. The risk may therefore not be inherent in any compound contained in mate, but due to temperature. The aim of this study was to quantitatively assess the risk of PAH in comparison with the risk of the temperature effect using the margin of exposure (MOE) methodology.
The human dietary benzo[a]pyrene (BaP) and PAH4 (sum of benzo[a]pyrene, benzo[a]anthracene, chrysene, and benzo[b]fluoranthene) exposure through consumption of coffee, mate, and tea was estimated. The oesophageal cancer risk assessment for both PAH and temperature was conducted using the MOE approach.
Considering differences in the transfer of the PAH from the leaves of mate and tea or from the ground coffee to the infusion, and considering the different preparation methods, exposures may vary considerably. The average individual exposure in μg/kg bw/day arising from consumption of 1 cup (0.2 L) of infusion was highest for mate (2.85E-04 BaP and 7.22E-04 PAH4). The average per capita exposure in μg/kg bw/day was as follows: coffee (4.21E-04 BaP, 4.15E-03 PAH4), mate (4.26E-03 BaP, 2.45E-02 PAH4), and tea (8.03E-04 BaP, 4.98E-03 PAH4). For all individual and population-based exposure scenarios, the average MOE for BaP and PAH4 was > 100,000 independent of beverage type. MOE values in this magnitude are considered as a very low risk. On the contrary, the MOE for the temperature effect was estimated as < 1 for very hot drinking temperatures, corroborating epidemiological observations about a probable oesophageal cancer risk caused by this behaviour.
The temperature effect but not PAH exposure may pose an oesophageal cancer risk. Consumer education on risks associated with consumption of ‘very hot’ beverages and policy measures to threshold serving temperatures should be discussed.
Consumption of very hot beverages (> 65 °C) is probably associated with an increased risk of oesophageal cancer [1, 2]. The first associations were found for the beverage mate (yerba mate). Mate is a herbal tisane-like beverage widely consumed in some South American countries, where the incidence of oesophageal squamous cell carcinoma (ESSC) is high, and notably higher than in other Latin American countries [3,4,5]. Currently, mate consumption is also gaining popularity worldwide. The infusion is made from dried leaves of Ilex paraguariensis A.St.-Hil. . Many mechanisms for carcinogenicity have been put forward. These include thermal injury, hyperthermia induced endogenous formation of nitrosamines, or impairment of mucosal barrier to entry of potential carcinogenic substances such as polycyclic aromatic hydrocarbons (PAH) . Later research on other beverage groups such as Camelia sinensis tea and coffee, which are also consumed very hot, resulted in reports of positive associations with increased risk of oesophageal cancer as well [8,9,10,11]. Considering that there is little evidence that these beverages contain high amounts of PAH , it has become the prevailing hypothesis that the risk might not be inherent to natural compounds contained in mate, or to compounds produced during processing, but rather due to a temperature effect. The first observations of associations between cancer and mate consumption may be explained by the fact that it is typically consumed extremely hot, with temperatures above 70–85 °C in some countries and regions, and in these areas, consumption levels are very high [1, 12]. For this reason, the International Agency for Research on Cancer (IARC) modified its evaluation from “hot mate drinking”, which was classified as probably carcinogenic (group 2A) since 1991 , into “very hot beverages at above 65 °C drinking” (group 2A) during its recent monograph (Vol. 116) meeting in 2016 . Nevertheless, an additional role for the PAH contents in mate on oesophageal cancer risk cannot be completely discarded yet.
PAH belong to a large group of over 100 different types of polyarenes that arise from incomplete combustion of organic matter [14,15,16]. Besides tobacco smoke, diet is the main source of PAH exposure in humans . PAH contamination in food results from processing practices such as smoking, roasting, grilling, or drying during preparation and also from environmental uptake . BaP is the traditional marker for PAH contamination but the European Food Safety Authority (EFSA) in 2008 recommended the use of both BaP and PAH4 (the sum of benzo[a]anthracene, chrysene, benzo[b]fluoranthene and BaP) as markers for the occurrence and toxicity of PAH in foods . The IARC classifies the compounds in the PAH4 subgroup as follows: BaP is carcinogenic to humans (group 1), while benzo[a]anthracene, chrysene and benzo[b]fluoranthene are “possibly carcinogenic to humans” (group 2B) . PAH4 is useful for evaluation of products that do not contain BaP. Following this, the EU established a regulation for the maximum levels of PAH4 in some food matrices . However, limits for coffee, mate and tea were not included in the regulation . According to epidemiological findings, the consumption of PAH contaminated food has been linked to the occurrence of malignant tumours of various sites including the oesophagus [3, 22,23,24].
To disentangle the role of PAH and high temperature, comparative risk assessment of both is needed. The EFSA has developed and recommends an approach known as the margin of exposure (MOE), for which doses of substances that have been observed to cause low but measurable harmful responses are compared with relevant substance specific dietary intake estimates in humans.
In this study, we used the MOE methodology to quantitatively assess the risk of PAH from mate, tea, or coffee, in comparison with the risk of high temperatures. We used occurrence data of BaP and PAH4 compounds in coffee, mate and tea (beans and leaves) and infusions to estimate exposure arising from moderate and heavy consumption.
Research on toxicity and occurrence data
Toxicity data on PAH were obtained by a computer-assisted literature search by researchers with qualifications in food science, chemistry, toxicology, epidemiology and cancer risk assessment. Searches were carried out in the following databases: PubMed, Toxnet and ChemIDplus (U.S. National Library of Medicine, Bethesda, MD), Web of Science (Thomson Scientific, Philadelphia, PA), and IPCS/INCHEM (International Programme on Chemical Safety/Chemical Safety Information from Intergovernmental Organizations, WHO, Geneva, Switzerland). Searches on the occurrence of PAH were carried out in August 2016 in the following databases: PubMed (US National Library of Medicine, Bethseda, MD), Web of Science (Thomson Scientific, Philadelphia, PA), Google scholar (Google, NY, USA) and SciELO - Scientific Electronic Library Online (FAPESP - BIREME, São Paulo SP - Brazil) using the key words ‘polycyclic aromatic hydrocarbons’ or ‘PAH AND coffee’ or ‘PAH AND tea’ or ‘PAH AND mate’ or ‘BaP AND coffee’ or ‘BaP AND tea’ or ‘BaP AND mate’ or ‘PAH4 AND coffee’ or ‘PAH4 AND tea’ or ‘PAH4 AND mate’. Additionally, own analytical data of 28 tea, 26 coffee and 9 mate samples were included. The samples were submitted between 2013 and 2016 as official samples in the German Federal State Baden-Württemberg and were analysed using the German reference procedure .
Efforts were made to include all available studies; this was accomplished by a hand search of the reference lists of all articles for any relevant studies not included in the databases. The references, including abstracts, were imported into Mendeley (Mendeley Inc., NY, USA) and the relevant articles were manually identified and obtained in full text. We did not identify any article, which was available as abstract only or which we were not able to obtain in full text. No unpublished study was identified.
To provide a meta analysis over all PAH studies, the analytical data of all single studies were combined and descriptive statistics were calculated. No weighting between studies was conducted. The values of non-detectable results were set at the limit of detection to obtain the final distributions. All distributions were separately calculated for the analytical data based on beans/leaves as well as for data based on infusion analyses (i.e. the beverages were analysed after brewing in the final form to be consumed). The infusion data for tea were insufficient to conduct such an analysis, however.
For data on temperature effect, the references cited in IARC monograph Vol. 116  were used and no additional literature searches were conducted.
Approach for risk assessment
The risk assessment was conducted according to the harmonised approach of the EFSA for the risk assessment of substances that are genotoxic and carcinogenic . The methodology for quantitative risk assessment was based on previous studies for comparative risk assessment of alcohol and drugs [27, 28]. Computation of exposure scenarios considered the amount of leaves/beans used for preparation of infusions in real life [6, 12, 29,30,31].
The MOE approach was used for the risk assessment [26, 32]. The benchmark dose (BMD), derived from animal cancer data by mathematical modelling within the observed range of experimental data, is recommended as a standardized reference point. To obtain the MOE, the Benchmark Dose Lower Confidence Limit (BMDL) of 10% was selected. The BMDL is an estimate of the lowest dose that is 95% certain to cause no more than a 10% cancer incidence in rodents. In general, BMDLs are used as the statistical lower confidence limits of benchmark doses to derive “safe” exposure levels .
The BMD doses for PAH are available in the literature, and the following BMDL values were used as points of departure for the cancer risk assessment: 70 μg/kg bw/day for BaP and 340 μg/kg bw/day for PAH4 based on the most recent EFSA report .
For temperature, no BMD modelling results were available in the literature so that own modelling was conducted based on literature studies. No adequate human data for dose-response modelling was identified. Two co-carcinogenicity studies in mice and rats were available. The study in mice reported the incidences of oesophageal lesions (hyperplasia, dysplasia, or oesophageal tumours) induced by N-nitrosodiethylamine (NDEA) in the presence of hot water (at 70 °C but not at 60 °C). The mouse study, however, did not provide raw data results suitable for BMD modelling. The study in rats reported data on N-nitrosomethylbenzylamine (NMBA)-induced oesophageal tumours, which was significantly increased at 65 °C versus 55 °C and the control group .
MOEs were calculated by dividing the reference point, i.e. the BMDL, by the estimated human intakes, with low MOE implying larger risks for humans. A threshold of 10,000 is often used to define public health risks, which considers species differences, human variability and additional uncertainties for substances that are genotoxic and carcinogenic .
The BMD(L) values were either taken from the literature search, or additionally BMD and BMDL values were calculated using the US EPA’s BMDS 220.127.116.11 software (available at the US EPA website: http://www.epa.gov/ncea/bmds/index.html). The exposures and MOE were then calculated using the software package @Risk for Excel Version 7.5.0 (Palisade Corporation, Ithaca, NY, USA). Monte Carlo simulations were performed with 10,000 iterations using Latin Hypercube sampling and Mersenne Twister random number generator. The distribution functions and detailed calculation methodology is specified in Additional file 1: Tables S1-S2 online.
Population-based dietary intake assessment and exposure scenarios
The EFSA harmonized approach has also been used for the dietary intake assessment analysis . The BaP and PAH4 content of coffee, tea and mate were estimated based on data from our literature review. Special exposure scenarios were developed for light, moderate and heavy drinkers. In order to assess individual exposure scenarios, average daily consumption at 1 cup (0.2 L) for light drinking, 2 cups (0.4 L) for moderate drinking, and one and three litres for heavy and very heavy consumers were considered. Figures of per capita consumption for population-based assessments were obtained from Ref. . The PAH exposure due to beverage consumption was then calculated considering these different individual as well as population-based exposure scenarios.
Occurrence of BaP and PAH4 in coffee, mate and tea
The occurrence data regarding BaP and PAH4 in coffee, tea and mate from 54 studies in total are summarized in Table 1. The final distributions over all data are provided in Table 2. The average contents of BaP in various matrices were as follows: mate leaves 142.9 μg/kg, tea leaves 8.71 μg/kg, coffee beans 4.12 μg/kg, mate infusion 0.072 μg/L and coffee infusion 0.004 μg/L; while the average content of PAH4 was as follows: mate leaves 872.6 μg/kg, tea leaves 150.53 μg/kg, coffee beans 22.29 μg/kg, mate infusion 0.214 μg/L and coffee infusion 0.041 μg/L (Table 2). Among the matrices of beans and leaves, mate leaves contained the highest amount of BaP and PAH4 while coffee beans had the lowest. The low number of studies on infusions must be noted, which may introduce bias especially as a majority of the infusion studies reported were on tea (n = 4, with most samples below the limit of detection making it impossible to make a distribution estimate) compared to only single studies on mate (n = 1) and coffee (n = 1).
The exposure for different individual and population-based scenarios was separately calculated for the analytical data based on beans/leaves considering the extraction percentage (Table 3) and the data based on direct analyses of infusions. The exposure scenarios are summarized in Table 4. The full distributions are provided in Additional file 1: Table S3.
Based on results for consumption of 1 cup (0.2 L) using the data for beans and leaves analyses, mate had the highest BaP exposures at 1.94E-03 μg/kg bw/day followed by coffee at 1.88E-04 μg/kg bw/day. Similarly, PAH4 exposure was highest for mate, 1.11E-02 μg/kg bw/day, and lowest for tea, 1.02E-03 μg/kg bw/day. For infusion analyses, the BaP exposure was highest for mate, 2.85E-04 μg/kg bw/day, and lowest for coffee while PAH4 exposure was lowest for mate, 7.22E-04 μg/kg bw/day as shown in Table 4. The per capita consumption data are reported as amount of beans/leaves only, so that for this approach only the analytical data from beans/leaves were used. The per capita exposures were as follows: coffee (4.21E-04 μg/kg bw/day BaP and 4.15E-03 μg/kg bw/day PAH4), mate (4.26E-03 μg/kg bw/day BaP and 2.45E-02 μg/kg bw/day PAH4), tea (8.03E-04 μg/kg bw/day BaP and 4.98E-03 μg/kg bw/day PAH4) as shown in Table 4.
Benchmark dose (BMD) modelling for temperature
The Additional file 2: data appendix provided as supplementary material shows the raw results of benchmark dose-response modelling of a study in rats. The results are summarized in Table 5. The BMD for NMBA-induced oesophageal tumours was 64 °C and the BMDL was 56 °C (Table 5).
Cancer risk assessment using the margin of exposure (MOE) approach
Based on results for consumption of 1 cup (0.2 L) in the matrices of beans and leaves, tea had the average highest MOE from BaP at 4.50E+06 followed by coffee at 4.48E+06 while MOE from PAH4 was highest for tea, 3.35E+06 and lowest for mate, 1.13E+06. For infusions, the BaP MOE was highest for coffee 4.61E+06 and lowest for mate 2.45E+06 while MOE from PAH4 was highest for mate, 6.85E+06 as shown in Table 4.
The per capita MOE was in the order of 1.97E+06, 1.39E+06, 4.91E+05 for coffee, tea and mate, respectively, and at the 95th percentile the order of MOE does not change (Table 4).
The margins of exposure for BaP (Fig. 1) and PAH4 (Fig. 2) for the different individual drinking scenarios are compared to the MOE threshold of 10,000. Only the extreme consumption of 3 L per day (meaning that the complete daily liquid requirement is fulfilled with these hot beverages) in worst-case scenarios (95th percentile contamination) would lead to MOEs below this threshold for all three beverages. Due to the high variance in the distributions, no significant beverage-specific influence was detected. Figure 3 shows the margin of exposure estimated for temperature. For 25 °C and 85 °C, the MOE would be 2.3 and 0.7, respectively. Hence, the MOE is below 1 at very hot temperatures.
The competing role of PAH and hot temperatures in hot beverages of widespread consumption has been hard to disentangle. This study reports a comprehensive worldwide systematic overview of the cancer risk assessment for PAH in comparison with very hot (> 65 °C) coffee, tea and mate. To our knowledge, this is the first study to quantitatively estimate cancer risk of PAH in mate compared to coffee and tea, and also to provide a comparison with the risk of the temperature effect using the margin of exposure methodology.
Mate is traditionally consumed either as hot mate or cold mate. The hot mate infusion consumed in South America is made by placing 20–50 g of mate in a vessel where very hot water (70–85 °C) is slowly poured over the material and after each pouring, the water is sucked through a special drinking straw fitted with a filter on the end immersed in the mate infusion [6, 12, 31]. Cold mate may be consumed in the same way or in a glass like a regular drink, at 4–8 °C . PAH in mate are produced by fire or high temperature exposure during traditional drying using direct fire and/or during roasting of the mate leaves [36,37,38,39,40]. The commercial processing of mate may involve two successive drying stages namely an initial rapid drying process at 400–750 °C by use of direct flames followed by final drying at 90–350 °C in rotating cylinders heated by burning wood before grinding .
PAH in the widely consumed beverage tea (Camellia sinensis) may result from air pollution in tea farms and during processing which may involve drying over burning wood, coal or oil . Coffee is also another non-alcoholic beverage consumed worldwide, and another possible source of PAH that can form through roasting, which is a crucial step in the processing of coffee for the development of aroma, flavour and colour [35, 41,42,43]. In green coffee, however, both absence and presence of PAH have been reported [43,44,45] and the eventual presence has been associated with air pollution rather than the high temperatures used during drying of the seeds.
Our finding of MOE values of more than 100,000 for PAH and less than 1 for very hot temperatures provides quantitative confirmation that “very hot” drinking temperatures are probably causally related to oesophageal cancer whereas low temperature mate is not. The very high MOE for PAH, which suggests that the human dosage is considerably below the lowest effective dose in animals, also makes the assumption implausible that PAH and temperature may act synergistically. Moreover, the very similar PAH exposure levels from mate and tea consumption offer indirect further support for IARC’s evaluation. In countries, such as the UK, where tea is consumed typically at lower temperature than the one of mate consumption in South America, increased rates of oesophageal cancer have not been observed. Also, it might be a matter of temperature dosage, mate drinkers might be drinking many more drinks per day than tea drinkers that possibly drink 1–2 cups per ‘session’, whereas a hot mate drinker can easily drink 1 l worth of water in a ‘session’, thus exposing the oesophagus many more times to very hot temperatures. Our assessment excludes assumptions of a beverage-specific effect, and confirms the limited epidemiological evidence that suggested a risk for very hot beverages independent of type. The mechanism of carcinogenicity of very hot beverages deserves further studies.
A limitation of our study is the fact that for some jurisdictions only a sparsity of data was available for PAH. Apart from Germany where data for all the beverage categories was available, most studies in other countries focused on only one type of beverage probably due to high prevalence of consumption in those countries. There were only extremely limited studies on PAH in infusions available, for mate and coffee only one single study. Due to the lack of data, even such single studies were included in our meta-analysis, specifically because we judge such direct beverage analyses more meaningful than bean/leaves analyses. However, we note that the same limitations generally apply to the assessments of the international agencies EFSA and Food and Agriculture Organization/WHO Expert Committee on Food Additives (JECFA), whose approaches we specifically followed.
In some samples, the occurrence of BaP was reported to be zero or below the limit of detection and this necessitated the use of the LOD value in order to calculate descriptive statistics. The sensitivity of the analytical methods used could explain the different LOD values reported. Analysis of PAH in foods such as coffee, mate and tea is complex due to the many compounds present in such food matrices, requiring sample preparative procedures such as liquid-liquid extraction (LLE), solid phase extraction (SPE), liquid phase microextraction (LPME), membrane assisted solvent microextraction, solid phase microextraction (SPME) and stir bar sorptive extraction (SBSE) . The extraction processes not only make the analysis expensive but may also compromise on recovery of PAH . Coupling the preparation step with analytical equipment is feasible such as SPME-HPLC or head space SPME-HPLC . Separation and quantification of the PAH in infusions is by liquid chromatography with fluorescent detection (LC-FLD) [12, 41, 48] or by GC-MS , while in ground coffee LC-FLD [44, 50], isotope dilution GC-MS , GC-MS/MS , or LC-UV , and in tea leaves, LC-UV , or GC-MS  have been utilized.
The substitution of “non detectable” data with LOD is a conservative approach taken and therefore the risk would have been overestimated (in relation to the approach of setting such values as 0). Some studies had limited sample sizes as low as 1 and this could skew meta-analysis results. Nevertheless, we think our distributions are robust. For example, no different judgments would arise with the calculation method of setting non detectable values as 0 (data not shown).
Most studies report the occurrence of PAH in mate and tea leaves and coffee beans (n = 48; 89%). However, since the raw leaves and coffee are usually not directly ingested, estimation of PAH concentrations in the beverages is necessary for exposure estimation. The differences in the manner of processing of the leaves for mate and tea and coffee beans and preparation of the beverages influences the levels of PAH in the beverages. The contribution of individual PAH to the total PAH concentration is dependent on its solubility. While the physicochemical data suggest that BaP and PAH are poorly soluble in water, the literature reports very variable extraction percentages from zero to above 50%. Furthermore, it must be considered that hot mate is typically consumed by pouring hot water on the same leaves repeated times, and very often replacing some of the leaves by fresh ones to keep the mate ‘strong’, therefore, the actual amount extracted from the leaves might be higher for mate than for the other beverages that are extracted only once. The actual preparation processes also influence the levels of PAH found in the drink, for example, low levels were found in cold mate compared to hot mate infusions . Literature reports extraction of BaP and PAH4 in the range 0.6–50% and 0–123% respectively (Table 3). The high transfer in tea reported has been linked to its high contents of essential oils, which may act as co-solvents that alter the physical and chemical properties of the water used for tea making and increase the affinity of PAH for the aqueous infusion [54, 55].
An investigation by Lin et al.  to identify the factors affecting the transfer of PAH singled out tea variety, tea/water ratio, and brewing time to affect the levels. The transfer average of all studies was 27% (which is very similar to the value of 26% used by EFSA for coffee ), but due to the probabilistic calculation method, the whole range of literature data was taken into account in our case. Another problem is the variety of brewing methods and amounts of tea leaves or coffee beans used per L of water. This may considerably vary depending on brewing method and individual preferences (e.g. see overview of brewing methods in ). For this reason, we have defined reasonable minimum and maximum values and used “no-knowledge” distributions in the probabilistic calculations to randomly consider the whole range of brewing methods. The appropriateness of this approach has been indirectly validated because the resulting exposure data are in reasonable agreement between both calculations methods (based on direct infusion analyses and recalculated from beans/leaves), so that we decided to average the values from both methods for the final data presentation (Figs. 1 and 2).
The higher average content of PAH and consequently the exposures due to mate and tea could be due to the similarities in processing of leaves that involve roasting over wood and that the leaves have larger surface areas that accumulate PAH from air . The age of the raw leaves collected from the farms also influences the level of PAH with the older cultivated mate having higher PAH as a result of longer exposure to environmental contamination .
Low amounts of low toxicity PAH have been identified in green coffee, mostly due to exposure to air pollution during sun drying of the beans. The roasting methods can also affect PAH formation in coffee. Drum roasters over direct fire enables uneven and higher temperatures in some parts of the beans while in fluidized bed roasters, frequently used currently, roasting temperatures do not exceed 240–260 °C, never reaching the temperature used for scorching mate. In general, the amount of PAH in coffee beans tend to increase in very dark roasts, only, when the beans are almost carbonized . Due to their low solubility, PAH concentrations in the beverage should be low. Paper filter should retain lipophilic compounds as opposed to extraction methods that do not use such filters such as espresso, Italian (Mocha), Turkish and Scandinavian boiled coffee . Higher amounts of PAH, however, have been identified in torrefacto coffee (roasted with sugar) . There was lack of data for PAH4 in some studies besides BaP; and in one study, occurrence of chrysene in tea leaves was reported as a sum of chrysene and triphenylene due to lack of resolution between the two PAH . However, this did not lead to higher values of PAH compared to other studies. Nevertheless, it could have contributed to risk overestimation in tea.
Compared to the large number of studies on PAH in hot beverages, and also to the toxicity data on PAH available, there is a sparsity of data on the temperature effect.
The weak data basis for thresholds for temperature are because temperature effect studies rely on self-reporting and there is subjectivity in the notion of different temperatures . The human epidemiological data provide evidence of a probable association between very hot beverages and oesophageal cancer, suggesting ~ 2-fold increase with various hot drinks [1, 2]. Importantly, a large pooled analysis was conducted that reported estimates of association between mate drinking and oesophageal cancer stratifying by temperature and showed little evidence of association when mate was consumed warm, regardless of the amount consumed, thus supporting a role of very hot temperatures and not mate per se .
Only two animal experiments were judged by IARC as adequate for evaluation, but they are only co-carcinogenicity studies. However, the animal experiments provide excellent corroborating evidence into a threshold, and the BMD of the animal experiment is in excellent agreement with the judgment of the IARC experts, that “very hot beverages” leading to a concern for oesophageal cancer, may start above 65 °C. However, there exists the problem of risk extrapolation from animal to human data especially in this case. The animals in the Li et al.  study received very hot water with an amount of 1 mL/kg for 5 times per week. Extrapolated to a human with a bodyweight of 60 kg this would amount to 60 mL of very hot beverage, which is not an unreasonable assumption, even less than one cup size. For this reason, we think that the temperature threshold of the animal study may be directly transferrable to humans. The normal practise in regulatory toxicology would be to add safety factors of 10 for both inter-species extrapolation and for intra-species differences to suggest acceptable daily intakes (ADI). In this case, we do not believe that it would make any sense to add safety factors for example to the BMD of 64 °C to derive such an ADI, because the ADI would clearly lie in a temperature zone (room temperature and below) that is perceived as being completely without risk. Till better studies become available, we think it is currently sensible to use the threshold value of IARC of 65 °C as “acceptable temperature”, which is also in excellent agreement to the BMD of 64 °C calculated in this study.
PAH are ubiquitous contaminants present in air, soil, water, foods; therefore, their role in carcinogenesis cannot be overlooked. The PAH exposure arising from consumption of coffee, mate and tea to cumulative exposure is rather small, according to EFSA, less than 10% of the combined exposure from all food groups. Nevertheless, contamination regulations usually demand that PAH contamination should be as low as technologically achievable; therefore, any amount of contamination in beverages should be avoided. Considering the limitations of the toxicological data but in light of the epidemiological evidence, the authors believe that risk management should be prioritized regarding the temperature problem rather than the PAH problem. Potential measures may include consumer education to drink at lower temperatures. However, similarly to alcohol consumption, changing behaviours about beverage temperatures might be difficult due to cultural traditions and preferences, since consumers may be more willing to tolerate higher risks for intentional behaviours (such as high temperature drinking) than for unintentional risks (such as PAH contamination) . Considering the increasing use of commercial hot beverage preparation machines, policy measures regarding serving temperatures (e.g. maximum temperatures for commercial hot beverage preparation machines) could also be effective.
Acceptable daily intake
Benchmark dose for a benchmark response of 10%
Benchmark Dose Lower Confidence Limit for a benchmark response of 10% (The BMDL is an estimate of the lowest dose that is 95% certain to cause no more than a 10% effect)
European Food and Safety Authority
Esophageal squamous cell carcinoma
Gas chromatography-mass spectrometry
Gas chromatography-tandem mass spectrometry
International Agency for Research on Cancer
International Programme on Chemical Safety/Chemical Safety Information
Joint FAO/WHO Expert Committee on Food Additives
Liquid phase microextraction
Margin of exposure
Polycyclic aromatic hydrocarbons
sum of benzo[a]pyrene, chrysene, benzo[b]fluroanthracene and benzo[a]anthracene
Stir bar sorptive extraction
Solid phase extraction
Solid phase microextraction
Toxic Equivalent Factors
- US EPA:
Unites States Environmental Protection Agency
IARC. IARC Monographs on the Evaluation of Carcinogenic risks to Humans, Vol. 116, Drinking coffee, maté, and very hot beverages. Lyon, France; in press.
Loomis D, Guyton KZ, Grosse Y, Lauby-Secretan B, El Ghissassi F, Bouvard V, et al. Carcinogenicity of drinking coffee, mate, and very hot beverages. Lancet Oncol. 2016;17:877–8.
Lubin JH, De Stefani E, Abnet CC, Acosta G, Boffetta P, Victora C, et al. Mate drinking and esophageal squamous cell carcinoma in South America: pooled results from two large multicenter case-control studies. Cancer Epidemiol Biomark Prev. 2014;23:107–16.
Fagundes RB, Abnet CC, Strickland PT, Kamangar F, Roth MJ, Taylor PR, et al. Higher urine 1-hydroxy pyrene glucuronide (1-OHPG) is associated with tobacco smoke exposure and drinking maté in healthy subjects from Rio Grande do Sul, Brazil. BMC Cancer. 2006;6:139.
Szymańska K, Matos E, Hung RJ, Wünsch-Filho V, Eluf-Neto J, Menezes A, et al. Drinking of maté and the risk of cancers of the upper aerodigestive tract in Latin America: a case-control study. Cancer Causes Control. 2010;21:1799–806.
Heck CI, De Mejia EG. Yerba mate tea (Ilex Paraguariensis): a comprehensive review on chemistry, health implications, and technological considerations. J Food Sci. 2007;72:R138–51.
Chen Y, Tong Y, Yang C, Gan Y, Sun H, Bi H, et al. Consumption of hot beverages and foods and the risk of esophageal cancer: a meta-analysis of observational studies. BMC Cancer. 2015;15:449.
Nubia M, Victoria CG, Crespi M, Saul C, Braga NM, Correa P. Hot mate drinking and precancerous lesions of the oesophagus: an endoscopic survey in southern Brazil. Int J Cancer. 1987;39:708–9.
Rolón PA, Castellsagué X, Benz M, Hot MN. Cold mate drinking and esophageal cancer in Paraguay. Cancer Epidemiol Biomark Prev. 1995;4:595–605.
Castellsagué X, Muñoz N, De Stefani E, Victora CG, Castelletto R, Rolón PA. Influence of mate drinking, hot beverages and diet on esophageal cancer risk in South America. Int J Cancer. 2000;88:658–64.
Victora CG, Munoz N, Day NE, Barcelos LB, Peccin DA, Braga NM. Hot beverages and esophageal cancer in southern Brazil: a case- control study. Int J Cancer. 1987;39(6):710.
Thea AE, Ferreira D, Brumovsky LA, Polycyclic SME. Aromatic hydrocarbons (PAHs) in yerba mate ( Ilex Paraguariensis St. Hil ) traditional infusions (mate and terere). Food Control. 2016;60:215–20.
IARC. IARC Monographs on the evaluation of Carcinogenic risks to Humans Vol. 51, coffee, tea, maté, methylxanthines and methylglyoxal. Lyon, France. 1991. Available online: http://monographs.iarc.fr/ENG/Monographs/vol51/mono51.pdf
Wenzl T, Simon R, Anklam E, Kleiner J. Analytical methods for polycyclic aromatic hydrocarbons (PAHs) in food and the environment needed for new food legislation in the European Union. TrAC - Trends Anal Chem. 2006;25:716–25.
Howard JW, Fazio T. Review of polycyclic aromatic hydrocarbons in foods. J Agric Food Chem. 1969;17:527–31.
Mumtaz MM, George JD, Gold KW, Cibulas WDC. ATSDR evaluation of health effects of chemicals. IV. Polycyclic aromatic hydrocarbons (PAHs): understanding a complex problem. Toxicol Ind Health. 1996;12(12):742–971.
Xia Z, Duan X, Qiu W, Liu D, Wang B, Tao S, et al. Health risk assessment on dietary exposure to polycyclic aromatic hydrocarbons (PAHs) in Taiyuan, China. Sci Total Environ. 2010;408:5331–7.
Bansal V, Review KK. Of PAH contamination in food products and their health hazards. Environ Int. 2015;84:26–38.
Authority EFS. Polycyclic aromatic hydrocarbons in food. EFSA J. 2008;724. https://doi.org/10.2903/j.efsa.2008.724.
IARC. List of classifications, Volumes 1–119. http://monographs.iarc.fr/ENG/Classification/latest_classif.php.
European Commission. Commission regulation (EU) no 835/2011 of 19 august 2011 amending regulation (EC) no 1881/2006 as regards maximum levels for polycyclic aromatic hydrocarbons in foodstuffs. Off J Eur Union. 2011;215:4–8.
Denissenko MF, Pao A, Tang M, Pfeifer GP. Preferential formation of benzo[a]pyrene adducts at lung cancer mutational hotspots in P53. Science. 1996;274:430–2.
Shen J, Terry MB, Gammon MD, Gaudet MM, Teitelbaum SL, Eng SM, et al. IGHMBP2 Thr671Ala polymorphism might be a modifier for the effects of cigarette smoking and PAH - DNA adducts to breast cancer risk. Breast Cancer Research Treat. 2006;99:1–7.
Olsson AC, Fevotte J, Fletcher T, Cassidy A, Mannetje A, Zaridze D, et al. Occupational exposure to polycyclic aromatic hydrocarbons and lung cancer risk : a multicenter study in Europe. Occup Env Med. 2010;67:98–103.
Anon. Amtliche Sammlung von Untersuchungsverfahren. BVL L 07.00–40. Untersuchung von Lebensmitteln - Bestimmung von Benzo(a)pyren in geräucherten und mit Raucharomen hergestellten Fleischerzeugnissen. Beuth-Verlag, Berlin, Germany. 2004.
EFSA. Opinion Of the scientific committee on a request from EFSA related to a harmonised approach for risk assessment of substances which are both genotoxic and carcinogenic. EFSA J. 2005;282:1–31. https://doi.org/10.2903/j.efsa.2005.282.
Lachenmeier DW, Rehm J. Comparative risk assessment of alcohol, tobacco, cannabis and other illicit drugs using the margin of exposure approach. Sci Rep. 2015;5:8126.
Pflaum T, Hausler T, Baumung C, Ackermann S, Kuballa T, Rehm J, et al. Carcinogenic compounds in alcoholic beverages: an update. Arch Toxicol. 2016;90:2349–67.
ISO 3103:1980. Tea - Preparation of liquor for use in sensory tests. Geneva: International Organization for Standardization. https://www.iso.org/standard/8250.html.
ISO 6668:2008. Green coffee - Preparation of samples for use in sensory analysis. Geneva: International Organization for Standardization. https://www.iso.org/standard/44609.html.
da Silveira TFF, Meinhart AD, Ballus CA, Godoy HT. The effect of the duration of infusion, temperature, and water volume on the rutin content in the preparation of mate tea beverages: an optimization study. Food Res Int. 2014;60:241–5.
EPA. The use of the benchmark dose approach in health risk assessment. Washington DC: EPA/630/R-94/007. EPA/630/R-94/007; 1995.
Kodell RL. Replace the NOAEL and LOAEL with the BMDL01 and BMDL10. Environ Ecol Stat. 2009;16:3–12.
Li ZG, Shimada Y, Sato F, Maeda M, Itami A, Kaganoi J, et al. Promotion effects of hot water on N-nitrosomethylbenzylamine-induced esophageal tumorigenesis in F344 rats. Oncol Rep. 2003;10:421–6. https://doi.org/10.3892/or.10.2.421.
Houessou JK, Maloug S, Leveque AS, Delteil C, Heyd B, Camel V. Effect of roasting conditions on the polycyclic aromatic hydrocarbon content in ground Arabica coffee and coffee brew. J Agric Food Chem. 2007;55:9719–26.
Kamangar F, Schantz MM, Abnet CC, Fagundes RB, Dawsey SM. High levels of carcinogenic polycyclic aromatic hydrocarbons in mate drinks. Cancer Epidemiol Biomark Prev. 2008;17:1262–8.
Vieira MA, Maraschin M, Rovaris ÂA, De RD, Amboni MC, Pagliosa CM, et al. Occurrence of polycyclic aromatic hydrocarbons throughout the processing stages of erva-mate (Ilex Paraguariensis). Food Addit Contam. 2010;27:776–82.
Kowalski BJ, Rigdon A, Cochran J. Analytical method for polycyclic aromatic hydrocarbons (PAHs) in yerba mate tea using modified QuEChERS, solid phase extraction and GC-TOFMS and GC-MS/MS. 2015. https://www.restek.com/pdfs/FFAN2086-UNV.pdf.
Golozar A, Fagundes RB, Etemadi A, Schantz MM, Kamangar F, Abnet CC, et al. Significant variation in the concentration of carcinogenic polycyclic aromatic hydrocarbons in yerba mate samples by brand, batch, and processing method. Environ Sci Technol. 2012;46:13488–93.
Winkler G, Damaty M, El SC, Lachenmeier DW. Mate - A “ new ” caffeinated raw material for the food and beverage industry. Ernahrungs Umschau. 2015;61:160–3. https://doi.org/10.4455/eu.2014.027.
Tfouni SAV, Serrate CS, Leme FM, Camargo MCR, Teles CRA, Cipolli KMVAB, et al. Polycyclic aromatic hydrocarbons in coffee brew: influence of roasting and brewing procedures in two Coffea cultivars. LWT Food Sci Technol. 2013;50:526–30. https://doi.org/10.1016/j.lwt.2012.08.015.
Sadowska-Rociek A, Surma M, Cieślik E. Determination of polycyclic aromatic hydrocarbons in coffee and coffee substitutes using dispersive SPE and gas chromatography-mass spectrometry. Food Anal Methods. 2015;8:109–21.
Perrone D, Rajoy RRS, Farah A. Effect of coffee roasting on the formation of polycyclic aromatic hydrocarbons. Dublin, Ireland: Proc. 2010 EFFoST Annu. Meet; 2010.
Badolato ESG, Martins MS, Aued-Pimentel S, Alaburda J, Kumagai EE, Baptista GG, et al. Sistematic study of benzo[a]pyrene in coffee samples. J Braz Chem Soc. 2006;17:989–93.
Moreira DP. Efeito da torrefação sobre aspectos benéficos e prejudiciais do café à saúde Humana: uma abordagem holística. DSc Thesis. Universidade Federal do Rio de Janeiro, Chemistry Institute, Rio de Janeiro, Brazil. 2009. p. 190.
Zuin VG, Montero L, Bauer C, Popp P. Stir bar sorptive extraction and high-performance liquid chromatography-fluorescence detection for the determination of polycyclic aromatic hydrocarbons in mate teas. J Chromatogr A. 2005;1091:2–10.
Viñas P, Campillo N, Aguinaga N, Pérez-Cánovas E, Hernández-Córdoba M. Use of headspace solid-phase microextraction coupled to liquid chromatography for the analysis of polycyclic aromatic hydrocarbons in tea infusions. J Chromatogr A. 2007;1164:10–7.
Adisa A, Jimenez A, Woodham C, Anthony K, Saleh MA. Determination of polycyclic aromatic hydrocarbons in dry tea. J Environ Sci Heal. 2016;50:552–9.
Duedahl-Olesen L, Navaratnam MA, Jewula J, Jensen AH. PAH in some Brands of tea and Coffee. Polycycl Aromat Compd. 2014;35:74–90.
Ishizaki A, Saito K, Hanioka N, Narimatsu S, Kataoka H. Determination of polycyclic aromatic hydrocarbons in food samples by automated on-line in-tube solid-phase microextraction coupled with high-performance liquid chromatography-fluorescence detection. J Chromatogr A. 2010;1217:5555–63.
Pissinatti R, Nunes CM, de Souza AG, Junqueira RG, de Souza SVC. Simultaneous analysis of 10 polycyclic aromatic hydrocarbons in roasted coffee by isotope dilution gas chromatography-mass spectrometry: optimization, in-house method validation and application to an exploratory study. Food Control. 2015;51:140–8.
Sadowska-Rociek A, Surma M. Comparison of different modifications on QuEChERS sample preparation method for PAHs determination in black, green, red and white tea. Env Sci Pollut Res. 2014;21:1326–38.
Bishnoi NR, Mehta U, Umashanker Sain GGP. Quantification of polycyclic aromatic hydrocarbons in tea and coffee samples of Mumbai city (India) by high performance liquid chromatography. Environ Monit Assess. 2005;107:399–406.
Schlemitz S, Pfannhauser W. Analysis of nitro-PAHs in food matrices by on-line reduction and high performance liquid chromatography. Food Addit Contam. 1996;13:969–77.
Pincemaille J, Schummer C, Heinen E, Determination MG. Of polycyclic aromatic hydrocarbons in smoked and non-smoked black teas and tea infusions. Food Chem. 2013;145:807–13.
Lin D, Zhu L. Polycyclic aromatic hydrocarbons: pollution and source analysis of a black tea. J Agric Food Chem. 2004;52:8268–71.
Farah A. Coffee constituents. In: Chu Y, editor Coffee: Emerging Health Effects and Disease Prevention IFT Press/Wiley-Blackwell; 2012. p. 21. https://doi.org/10.1002/9781119949893.ch2.
Kayali-Sayadi MN, Rubio-Barroso S, Cuesta-Jimenez MP, Polo-Díez LM. A new method for the determination of selected PAHs in coffee brew samples by HPLC with fluorimetric detection and solid-phase extraction. J Liq Chromatogr Relat Technol. 1999;22:615–7.
Ziegenhals K, Jira W, Speer K. Polycyclic aromatic hydrocarbons (PAH) in various types of tea. Eur Food Res Technol. 2008;228:83–91.
Rehm J, Lachenmeier DW, Room R. Why does society accept a higher risk for alcohol than for other voluntary or involuntary risks? BMC Med. 2014;12:189.
Houessou JK, Benac C, Delteil C, Camel V. Determination of polycyclic aromatic hydrocarbons in coffee brew using solid-phase extraction. J Agric Food Chem. 2005;53:871–9.
Grover IS, Sharma R, Singh S, Pal B. Polycyclic aromatic hydrocarbons in some grounded coffee brands. Environ Monit Assess. 2013;185:6459–63.
Iwegbue CMA, Agadaga H, Bassey FI, Overah LC, Tesi GO, Nwajei GE. Concentrations and profiles of polycyclic aromatic hydrocarbons in some commercial Brands of tea-, coffee-, and cocoa-based food drinks in Nigeria. Int J Food Prop. 2015;18:2124–33.
Stanciu G, Birghila S, Dobrinas S. Residues of polycyclic aromatic hydrocarbons in different types of coffee. Sci Study Res. 2008;9:325–9.
Hpas A, Dieta NA, População DA, Cristiane M, De Camargo R, Cecília M, et al. Chá-Mate e Café Como Fontes de Hidrocarbonetos Policíclicos. Ciênc Tecnol Aliment. 2002;22:49–53.
Lee K, Shin H. Determination of polycyclic aromatic hydrocarbons in commercial roasted coffee beans. Food Sci Biotechnol. 2010;19:1435–40.
Jimenez A, Adisa A, Woodham C, Saleh M. Determination of polycyclic aromatic hydrocarbons in roasted coffee. J Environ Sci Health B. 2014;49:828–35. https://doi.org/10.1080/03601234.2014.938552.
Schulz CM, Fritz H, Ruthenschrör A. Occurrence of 15 + 1 EU priority polycyclic aromatic hydrocarbons (PAH) in various types of tea (Camellia Sinensis) and herbal infusions. Food Addit Contam. 2014;31:1723–35.
Lin D, Tu Y, Zhu L. Concentrations and health risk of polycyclic aromatic hydrocarbons in tea. Food Chem Toxicol. 2005;43:41–8.
Singh S, Vashishth A. PAHs in some brands of tea. Environ Monit Assess. 2011;177:35–8.
Martena MJ, Grutters MMP, De Groot HNKE, Monitoring RI. Of polycyclic aromatic hydrocarbons (PAH) in food supplements containing botanicals and other ingredients on the Dutch market. Food Addit Contam. 2011;28:925–42.
Shi Y, Wu H, Wang C, Guo X, Du J, Du L. Determination of polycyclic aromatic hydrocarbons in coffee and tea samples by magnetic solid-phase extraction coupled with HPLC-FLD. Food Chem. 2016;199:75–80.
Fiedler H. PCDD/PCDF, chlorinated pesticides and PAH in Chinese teas. Chemosphere. 2002;46:1429–33.
Khiadani M, Amin MM, Beik FM, Ebrahimi A, Farhadkhani M, Mohammadi-Moghadam F. Determination of polycyclic aromatic hydrocarbons concentration in eight brands of black tea which are used more in Iran. Int J Environ Health Eng. 2013;2:40. https://doi.org/10.4103/2277-9183.122427.
Garcia Londoño VA, Reynoso CM, Resnik SL. Polycyclic aromatic hydrocarbons (PAHs) survey on tea (Camellia Sinensis) commercialized in Argentina. Food Control. 2015;50:31–7.
Drabova L, Pulkrabova J, Kalachova K, Tomaniova M, Kocourek V, Hajslova J. Rapid determination of polycyclic aromatic hydrocarbons (PAHs ) in tea using two-dimensional gas chromatography coupled with time of flight mass spectrometry. Talanta. 2012;100:207–16.
Li X, Li N, Luo H, Lin L, Zou Z, Jia Y, et al. A novel synchronous fluorescence spectroscopic approach for the rapid determination of three polycyclic aromatic hydrocarbons in tea with simple microwave-assisted pretreatment of sample. J Agric Food Chem. 2011;59:5899–905.
Kayali-Sayadi MN, Rubio-Barroso S, Cuesta-Jimenez MP, Polo-Díez LM. Rapid determination of polycyclic aromatic hydrocarbons in tea infusion samples by high-performance liquid chromatography and fluorimetric detection based on solid-phase extraction. Analyst. 1998;123:2145–8.
De KN, Ton S, Stegen GHD. Rapid determination of benzoaIpyrene in roasted coffee and coffee brew by high-performance liquid chromatography with fluorescence detection. J Agric Food Chem. 1987;35:545–9.
Mohammadi-Moghadam F, Amin M, Beik F, Ebrahimi A, Farhadkhani M, Khiadani (Hajian) M. Determination of polycyclic aromatic hydrocarbons concentration in eight brands of black tea which are used more in Iran. Int J Environ Health Eng. 2013;2:40.
Maier HG. Carcinogenic compound content in coffee beans. Café Cacao Thé. 1991;35:133–42.
No specific funding was received for this study. The DAAD (German Academic Exchange Service) is thanked for granting a PhD scholarship to AOO, funding no. 57250960.
Availability of data and materials
All data generated or analysed during this study are included in this published article [and its supplementary information files].
Ethics approval and consent to participate
Not applicable (no experiments on live vertebrates (or higher invertebrates), humans or human samples).
Consent for publication
The authors declare that they have no competing interests
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary Tables S1-S3. Supplementary Table S1. Distribution functions as input for probabilistic analysis; Supplementary Table S2. Detailed calculation methodology for probabilistic risk assessment of PAH due to beverage consumption; Supplementary Table S3. Raw results of probabilistic estimation of BaP and PAH4 exposure and margin of exposure (MOE) using 10,000 iterations. (DOC 250 kb)
Data appendix with raw result for benchmark dose-response modelling. Supplementary Data S1. BMD modelling for mean number of tumours per rat from Li et al. (2003); Supplementary Data S2. BMD modelling for mean volume of tumours (mm3) from Li et al. (2003). (DOC 81 kb)