# Lebanese Population Exposure to Trace Elements via White Bread Consumption

Nada Lebbos, Claude Daou, Rosette R. Ouaini, Hanna H. Chebib, Michel Afram, Pierre Curmi, Laurence Dujourdy, Elias Bou-Maroun, Marie-Christine Chagnon

► To cite this version:

Nada Lebbos, Claude Daou, Rosette R. Ouaini, Hanna H. Chebib, Michel Afram, et al., Lebanese Population Exposure to Trace Elements via White Bread Consumption. Foods, 2019, 8 (11), pp.574. 10.3390/foods8110574. hal-02380322

HAL Id: hal-02380322

https://institut-agro-dijon.hal.science/hal-02380322

Submitted on 30 Apr 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Lebanese Population Exposure to Trace Elements via White Bread Consumption

Nada Lebbos<fnref n="1" />, Claude Daou<fnref n="2" />, Rosette Ouaini<fnref n="2" />, Hanna Chebib<fnref n="2" />, Michel Afram<fnref n="1" />, Pierre Curmi<fnref n="3" />, Laurence Dujourdy<fnref n="4" /> ID, Elias Bou-Maroun<fnref n="5" />,\* ID and Marie-Christine Chagnon<fnref n="6" />

Received: 17 October 2019; Accepted: 9 November 2019; Published: 14 November 2019

<figure><img src="image_1.png" /></figure>

check for
updates

**Abstract:** The objective of this study was to assess Lebanese population exposure to trace elements (TEs) via white pita consumption. A survey of white pita consumption was achieved among one thousand Lebanese individuals, grouped into adults (above 15 years old, men, and women) and young people (6–9 and 10–14 years old). The most consumed pita brands, labeled B1, B2, and B3, were selected. Levels of TEs (i.e., As, Cd, Co, Cr, Hg, Ni, and Pb) in B1, B2, B3 pitas were measured. The highest contents of TEs in pitas were: Ni (1292 µg/kg) and Co (91 µg/kg) in B1; As (400 µg/kg) and Cd (< 15 µg/kg) in B2; Cr (363 µg/kg), Pb (260 µg/kg), and Hg (0.89 µg/kg) in B3. The pita brand B3 was the source of the highest TEs exposure, except for Ni for which it was B1. Daily exposures to TEs due to the fact of pita consumption were compared to safety levels. There were no safety concerns for Hg, Cd, Cr or Co (except the 95th percentile of 6–9 years old). An excess of the Ni tolerable daily intake was observed for the most exposed populations. The very low margins of exposure for As and Pb suggest a worrying risk for the Lebanese population.

**Keywords:** survey; trace elements analysis; Lebanese pita; human exposure; risk characterization

## 1. Introduction

Bread is a source of energy and minerals, some of them being essential for the body. Essential minerals found in Lebanese bread are Cu, Fe, Mn, Zn, and Se which account for 36%, 25%, 30%, 22%, and 32% of respective recommended dietary allowance (RDA) values given for these elements [1].

Besides the quality of bread as nutriment, its safety must be assured and monitored in terms of putative contaminants like trace elements (TEs).

Indeed, trace elements represent the main food chemical contamination and are globally recognized as a public health hazard [2]. They occur either from natural process (lithogenic and pedogenic routes) and/or anthropogenic activities. The latter include industrial, agricultural (pesticides, fertilizers), and untreated wastewater [3]. This leads to the contamination of the different environmental matrices, including water, soil, air, and plants [4,5]. In general, plants can uptake TEs mainly by absorbing

<footnote n="1">¹ Lebanese Agricultural Research Institute LARI, Heavy metals Department, Beirut PO Box 90-1965, Lebanon; lebbosnada@hotmail.com (N.L.); lari@lari.gov.lb (M.A.)</footnote>

<footnote n="2">² Lebanese University, Faculty of Science II, Laboratoire d'Analyse Chimique (LAC), Fanar PO Box 90656, Lebanon; claudedaou@hotmail.com (C.D.); rouaini@ul.edu.lb (R.O.); hchebib@hotmail.com (H.C.)</footnote>

<footnote n="3">³ AgroSup Dijon, Univ. Bourgogne Franche-Comté, Biogéosciences UMR 6282, F-21000 Dijon, France; pierre.curmi@agrosupdijon.fr</footnote>

<footnote n="4">⁴ Service d'Appui à la recherche, AgroSup Dijon, F-21000 Dijon, France; laurence.dujourdy@agrosupdijon.fr</footnote>

<footnote n="5">⁵ AgroSup Dijon, Univ. Bourgogne Franche-Comté, PAM UMR A 02.102, Procédés Alimentaires et Microbiologiques, F-21000 Dijon, France</footnote>

<footnote n="6">⁶ Univ. Bourgogne Franche-Comté, INSERM U1231, NUTOX, Derttech "Packtox", AgroSup Dijon, F-21000 Dijon, France; marie-christine.chagnon@u-bourgogne.fr</footnote>

<footnote>\* Correspondence: elias.bou-maroun@agrosupdijon.fr; Tel.: +33-3-80-77-40-80</footnote>
Article

# Lebanese Population Exposure to Trace Elements via White Bread Consumption

Nada Lebbos<fnref n="1" />, Claude Daou<fnref n="2" />, Rosette Ouaini<fnref n="2" />, Hanna Chebib<fnref n="2" />, Michel Afram<fnref n="1" />, Pierre Curmi<fnref n="3" />, Laurence Dujourdy <fnref n="4" />[D], Elias Bou-Maroun <fnref n="5" />,*[D] and Marie-Christine Chagnon<fnref n="6" />

Received: 17 October 2019; Accepted: 9 November 2019; Published: 14 November 2019

<figure><img src="image_2.png" /></figure>

check for
updates

**Abstract:** The objective of this study was to assess Lebanese population exposure to trace elements (TEs) via white pita consumption. A survey of white pita consumption was achieved among one thousand Lebanese individuals, grouped into adults (above 15 years old, men, and women) and young people (6–9 and 10–14 years old). The most consumed pita brands, labeled B1, B2, and B3, were selected. Levels of TEs (i.e., As, Cd, Co, Cr, Hg, Ni, and Pb) in B1, B2, B3 pitas were measured. The highest contents of TEs in pitas were: Ni (1292 µg/kg) and Co (91 µg/kg) in B1; As (400 µg/kg) and Cd (< 15 µg/kg) in B2; Cr (363 µg/kg), Pb (260 µg/kg), and Hg (0.89 µg/kg) in B3. The pita brand B3 was the source of the highest TEs exposure, except for Ni for which it was B1. Daily exposures to TEs due to the fact of pita consumption were compared to safety levels. There were no safety concerns for Hg, Cd, Cr or Co (except the 95th percentile of 6–9 years old). An excess of the Ni tolerable daily intake was observed for the most exposed populations. The very low margins of exposure for As and Pb suggest a worrying risk for the Lebanese population.

**Keywords:** survey; trace elements analysis; Lebanese pita; human exposure; risk characterization

## 1. Introduction

Bread is a source of energy and minerals, some of them being essential for the body. Essential minerals found in Lebanese bread are Cu, Fe, Mn, Zn, and Se which account for 36%, 25%, 30%, 22%, and 32% of respective recommended dietary allowance (RDA) values given for these elements [1].

Besides the quality of bread as nutriment, its safety must be assured and monitored in terms of putative contaminants like trace elements (TEs).

Indeed, trace elements represent the main food chemical contamination and are globally recognized as a public health hazard [2]. They occur either from natural process (lithogenic and pedogenic routes) and/or anthropogenic activities. The latter include industrial, agricultural (pesticides, fertilizers), and untreated wastewater [3]. This leads to the contamination of the different environmental matrices, including water, soil, air, and plants [4,5]. In general, plants can uptake TEs mainly by absorbing

<footnote n="1">¹ Lebanese Agricultural Research Institute LARI, Heavy metals Department, Beirut PO Box 90-1965, Lebanon; lebbosnada@hotmail.com (N.L.); lari@lari.gov.lb (M.A.)</footnote>

<footnote n="2">² Lebanese University, Faculty of Science II, Laboratoire d'Analyse Chimique (LAC), Fanar PO Box 90656, Lebanon; claudedaou@hotmail.com (C.D.); rouaini@ul.edu.lb (R.O.); hchebib@hotmail.com (H.C.)</footnote>

<footnote n="3">³ AgroSup Dijon, Univ. Bourgogne Franche-Comté, Biogéosciences UMR 6282, F-21000 Dijon, France; pierre.curmi@agrosupdijon.fr</footnote>

<footnote n="4">⁴ Service d'Appui à la recherche, AgroSup Dijon, F-21000 Dijon, France; laurence.dujourdy@agrosupdijon.fr</footnote>

<footnote n="5">⁵ AgroSup Dijon, Univ. Bourgogne Franche-Comté, PAM UMR A 02.102, Procédés Alimentaires et Microbiologiques, F-21000 Dijon, France</footnote>

<footnote n="6">⁶ Univ. Bourgogne Franche-Comté, INSERM U1231, NUTOX, Derttech "Packtox", AgroSup Dijon, F-21000 Dijon, France; marie-christine.chagnon@u-bourgogne.fr</footnote>

<footnote>* Correspondence: elias.bou-maroun@agrosupdijon.fr; Tel.: +33-3-80-77-40-80</footnote>
them from polluted soils, and transporting them further to its different parts (roots, stem, leaves, and grains) [6]. The bioaccumulation of TEs in plants, especially those that enter the food chain, can be dangerous to human health due to the fact of their eventual toxicity [5].

Traces elements such as lead, arsenic, cadmium, and mercury are relevant toxic elements [7]. According to the existing literature, it is mandatory to quantify the contents of trace elements in foodstuffs in order to estimate dietary exposure, an essential step for risk assessment related to food consumption [2,8].

Currently, many studies have been undertaken in different countries to assess exposure to TEs occurring in their relative foodstuffs such as France [9], Spain [2], Italy [10], China [11], Saudi Arabia [12], and Lebanon [13]. In addition, the European Food Safety Authority (EFSA) panel on contaminants in the food chain has given its scientific opinion on the risks to public health related to the presence of Ni, Hg, and As in food [14–16]. Such studies should be representative of specific food as consumed by the population. In France, total diet studies (TDS) are national surveys routinely realized to assess public health risks associated with substances such as contaminants in food. They underline the necessity to reduce exposure to elements such as Cd, Al, Ni, As, and Co by a diversified diet (food items and origins) [17].

Nasreddine et al. [13] studied, using the TDS approach, the dietary exposure of the adult Lebanese urban population to six essential micro-nutriments (i.e., Co, Cu, Fe, Mn, Ni, and Zn) and two toxic heavy metals (i.e., Cd and Pb). This study constitutes a first estimate of the consumer exposure to trace elements through the diet in Lebanon. These authors determined that the average urban adult daily bread consumption was 136.8 g/day and that mean exposure was found to be satisfactory for most trace elements. However, the survey included only adult people living in the district of Beirut, thus preventing the extrapolation of the results to the whole country; indeed, agricultural areas are sections of high bread consumption.

Previous studies in Egypt, Iran, South Africa, Nigeria, and Portugal determined TEs contents in bread [18–22]. In Lebanon, Bou Khouzam and co-authors [1] determined the contents of 20 elements in three varieties of Lebanese bread (i.e., white, brown pita, and Saj bread) throughout the country during the wet and dry seasons. Toxic elements such as Cd and Pb were below maximum levels set by the Lebanese regulatory agency, (The Lebanese Standards Institution—Libnor) [23].

Bread is an integral part of the daily diet, but the types of bread and amounts ingested vary according to the age, sex, and socio-economic status of the individual. As a result, and as consumer behavior evolves over time, there is a need for up-to-date information on daily ETs intake via Lebanese bread consumption throughout Lebanon, including rural areas.

The Lebanese standards [23] provide only Pb, Cu, and Cd maximum levels in bread as food item contamination vectors. However, other TEs could be present [1,13]. As exposition represents the combination of consumption of bread potentially contaminated, it is of first interest to compare exposure with safety levels in order to monitor any safety preoccupation concerns for the Lebanese population. For this reason, this study was directed by a public institute (Lebanese Agricultural Research Institute—LARI) to assess exposure to trace elements via Lebanese pita consumption.

As, Cd, Co, Cr, Hg, Ni, and Pb were selected in this study based on previous studies of trace element composition of Lebanese bread [1,13] and based on the Lebanese bread standard Libnor, NL 240 [23]. In contrast to heavy metals known to be toxic, such as Pb, Cd or Hg, some elements, such as Cr and Co, are essential elements giving benefit to the body, the liver being a valuable source of elements [24]. Co is a transition metal with two oxidation states Co (II) and Co (III) and is essential as part of vitamin B12 (cobalamin) involved in folate and fatty acid and the metabolism of proteins and nucleic acids [25]. In regards to Cr, the Cr (III) form has been postulated to be necessary for the efficacy of the insulin-regulating metabolism of carbohydrates, lipids, and proteins [26]. However, EFSA [27] in 2014 concluded there was no evidence of a beneficial effect associated with Cr intake in healthy, normoglycemic subjects.
them from polluted soils, and transporting them further to its different parts (roots, stem, leaves, and grains) [6]. The bioaccumulation of TEs in plants, especially those that enter the food chain, can be dangerous to human health due to the fact of their eventual toxicity [5].

Traces elements such as lead, arsenic, cadmium, and mercury are relevant toxic elements [7]. According to the existing literature, it is mandatory to quantify the contents of trace elements in foodstuffs in order to estimate dietary exposure, an essential step for risk assessment related to food consumption [2,8].

Currently, many studies have been undertaken in different countries to assess exposure to TEs occurring in their relative foodstuffs such as France [9], Spain [2], Italy [10], China [11], Saudi Arabia [12], and Lebanon [13]. In addition, the European Food Safety Authority (EFSA) panel on contaminants in the food chain has given its scientific opinion on the risks to public health related to the presence of Ni, Hg, and As in food [14–16]. Such studies should be representative of specific food as consumed by the population. In France, total diet studies (TDS) are national surveys routinely realized to assess public health risks associated with substances such as contaminants in food. They underline the necessity to reduce exposure to elements such as Cd, Al, Ni, As, and Co by a diversified diet (food items and origins) [17].

Nasreddine et al. [13] studied, using the TDS approach, the dietary exposure of the adult Lebanese urban population to six essential micro-nutriments (i.e., Co, Cu, Fe, Mn, Ni, and Zn) and two toxic heavy metals (i.e., Cd and Pb). This study constitutes a first estimate of the consumer exposure to trace elements through the diet in Lebanon. These authors determined that the average urban adult daily bread consumption was 136.8 g/day and that mean exposure was found to be satisfactory for most trace elements. However, the survey included only adult people living in the district of Beirut, thus preventing the extrapolation of the results to the whole country; indeed, agricultural areas are sections of high bread consumption.

Previous studies in Egypt, Iran, South Africa, Nigeria, and Portugal determined TEs contents in bread [18–22]. In Lebanon, Bou Khouzam and co-authors [1] determined the contents of 20 elements in three varieties of Lebanese bread (i.e., white, brown pita, and Saj bread) throughout the country during the wet and dry seasons. Toxic elements such as Cd and Pb were below maximum levels set by the Lebanese regulatory agency, (The Lebanese Standards Institution—Libnor) [23].

Bread is an integral part of the daily diet, but the types of bread and amounts ingested vary according to the age, sex, and socio-economic status of the individual. As a result, and as consumer behavior evolves over time, there is a need for up-to-date information on daily ETs intake via Lebanese bread consumption throughout Lebanon, including rural areas.

The Lebanese standards [23] provide only Pb, Cu, and Cd maximum levels in bread as food item contamination vectors. However, other TEs could be present [1,13]. As exposition represents the combination of consumption of bread potentially contaminated, it is of first interest to compare exposure with safety levels in order to monitor any safety preoccupation concerns for the Lebanese population. For this reason, this study was directed by a public institute (Lebanese Agricultural Research Institute—LARI) to assess exposure to trace elements via Lebanese pita consumption.

As, Cd, Co, Cr, Hg, Ni, and Pb were selected in this study based on previous studies of trace element composition of Lebanese bread [1,13] and based on the Lebanese bread standard Libnor, NL 240 [23]. In contrast to heavy metals known to be toxic, such as Pb, Cd or Hg, some elements, such as Cr and Co, are essential elements giving benefit to the body, the liver being a valuable source of elements [24]. Co is a transition metal with two oxidation states Co (II) and Co (III) and is essential as part of vitamin B12 (cobalamin) involved in folate and fatty acid and the metabolism of proteins and nucleic acids [25]. In regards to Cr, the Cr (III) form has been postulated to be necessary for the efficacy of the insulin-regulating metabolism of carbohydrates, lipids, and proteins [26]. However, EFSA [27] in 2014 concluded there was no evidence of a beneficial effect associated with Cr intake in healthy, normoglycemic subjects.
Two risks may occur with trace elements: one is insufficient intake and one is excess intake. In this study, we focused on risk characterization due to the fact of excess intake via pita consumption.

Therefore, the objective of this study was to assess exposure to the selected trace elements (TEs) of the different Lebanese population categories via pita consumption. First, a country survey of pita consumption was achieved among one thousand individuals, grouped into adults (above 15 years old, men and women) and young people (6–9 and 10–14 years old). Second, levels of the seven selected trace elements in pitas from the three most consumed brands were measured. Third, using the survey data, population categories of exposures to TEs via white pita consumption were calculated and compared to the international health-based guidance values: tolerable daily intake (TDI), tolerable weekly intake (TWI) or toxicological reference points, such as bench mark dose limit (BMDL), to evaluate safety concerns.

## 2. Materials and Methods

### 2.1. Survey

A questionnaire survey about the daily consumption of pita bread was conducted on a target sample of 1000 people. They were selected randomly in 50 supermarkets over the five administrative regions of Lebanon, taking into consideration the population distribution [28]: 41% in Mount Lebanon, 20% in North Lebanon including Akkar, 17% in South Lebanon including Nabatieh, 13% in the Bekaa, and 9% in Beirut (Table 1 and Figure 1). The individuals were represented by 50% men, 50% women, and grouped as follows: 800 adults (> 15 years old) and 200 young people (between 6 and 15 years old).

<table><caption>**Table 1.** Sampling strategy according to the population distribution.</caption><thead><tr><th rowspan="2">Regions</th><th colspan="2">Population Distribution *</th><th colspan="2">Estimated Sample</th><th colspan="2">Surveyed Supermarkets</th><th colspan="2">Realized Questionnaires</th></tr><tr><th>n</th><th>%</th><th>n</th><th>%</th><th>n</th><th>%</th><th>n</th><th>%</th></tr></thead><tbody><tr><td>Beirut</td><td>361,366</td><td>9.6</td><td>90</td><td>9.0</td><td>5</td><td></td><td>82</td><td>8.3</td></tr><tr><td>Mount Lebanon</td><td>1,484,474</td><td>39.5</td><td>410</td><td>41.0</td><td>20</td><td></td><td>410</td><td>41.3</td></tr><tr><td>North Lebanon</td><td>763,712</td><td>20.3</td><td>200</td><td>20.0</td><td>10</td><td></td><td>200</td><td>20.2</td></tr><tr><td>South Lebanon</td><td>659,718</td><td>17.6</td><td>170</td><td>17.0</td><td>8</td><td></td><td>170</td><td>17.1</td></tr><tr><td>Bekaa</td><td>489,865</td><td>13.0</td><td>130</td><td>13.0</td><td>7</td><td></td><td>130</td><td>13.1</td></tr><tr><td>Total</td><td>3,759,136</td><td>100</td><td>1000</td><td>100</td><td>50</td><td></td><td>992</td><td>100</td></tr></tbody></table>

<figure><img src="image_2.png" /><figcaption>**Figure 1.** Location of the fifty surveyed supermarkets in Lebanon within the five administrative regions according to population density.</figcaption></figure>

<footnote>* Lebanese Population, Central Administration for Statistics (CAS) Living Conditions Survey 2007 [28].</footnote>
Two risks may occur with trace elements: one is insufficient intake and one is excess intake. In this study, we focused on risk characterization due to the fact of excess intake via pita consumption.

Therefore, the objective of this study was to assess exposure to the selected trace elements (TEs) of the different Lebanese population categories via pita consumption. First, a country survey of pita consumption was achieved among one thousand individuals, grouped into adults (above 15 years old, men and women) and young people (6–9 and 10–14 years old). Second, levels of the seven selected trace elements in pitas from the three most consumed brands were measured. Third, using the survey data, population categories of exposures to TEs via white pita consumption were calculated and compared to the international health-based guidance values: tolerable daily intake (TDI), tolerable weekly intake (TWI) or toxicological reference points, such as bench mark dose limit (BMDL), to evaluate safety concerns.

## 2. Materials and Methods

### 2.1. Survey

A questionnaire survey about the daily consumption of pita bread was conducted on a target sample of 1000 people. They were selected randomly in 50 supermarkets over the five administrative regions of Lebanon, taking into consideration the population distribution [28]: 41% in Mount Lebanon, 20% in North Lebanon including Akkar, 17% in South Lebanon including Nabatieh, 13% in the Bekaa, and 9% in Beirut (Table 1 and Figure 1). The individuals were represented by 50% men, 50% women, and grouped as follows: 800 adults (> 15 years old) and 200 young people (between 6 and 15 years old).

<table><caption>**Table 1.** Sampling strategy according to the population distribution.</caption><thead><tr><th rowspan="2">Regions</th><th colspan="2">Population Distribution *</th><th colspan="2">Estimated Sample</th><th colspan="2">Surveyed Supermarkets</th><th colspan="2">Realized Questionnaires</th></tr><tr><th>n</th><th>%</th><th>n</th><th>%</th><th>n</th><th>%</th><th>n</th><th>%</th></tr></thead><tbody><tr><td>Beirut</td><td>361,366</td><td>9.6</td><td>90</td><td>9.0</td><td>5</td><td></td><td>82</td><td>8.3</td></tr><tr><td>Mount Lebanon</td><td>1,484,474</td><td>39.5</td><td>410</td><td>41.0</td><td>20</td><td></td><td>410</td><td>41.3</td></tr><tr><td>North Lebanon</td><td>763,712</td><td>20.3</td><td>200</td><td>20.0</td><td>10</td><td></td><td>200</td><td>20.2</td></tr><tr><td>South Lebanon</td><td>659,718</td><td>17.6</td><td>170</td><td>17.0</td><td>8</td><td></td><td>170</td><td>17.1</td></tr><tr><td>Bekaa</td><td>489,865</td><td>13.0</td><td>130</td><td>13.0</td><td>7</td><td></td><td>130</td><td>13.1</td></tr><tr><td>Total</td><td>3,759,136</td><td>100</td><td>1000</td><td>100</td><td>50</td><td></td><td>992</td><td>100</td></tr></tbody></table>

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Location of the fifty surveyed supermarkets in Lebanon within the five administrative regions according to population density.</figcaption></figure>

<footnote>* Lebanese Population, Central Administration for Statistics (CAS) Living Conditions Survey 2007 [28].</footnote>
Data were collected over a two-month period, from 14 April 2017 to 5 June 2017. The team of interviewers was composed of 8 data collectors trained prior to the fieldwork by a supervisor. Participants were eligible for the study if they were Lebanese, had been living in Lebanon for at least fifteen years before the date of the survey, and did not have a chronic or serious illness.

In each supermarket, eligible participants were interviewed face-to-face by a data collector who instantly filled an offline questionnaire on an electronic tablet, using multiple-choice questions and figures (e.g., size of the pita). In the case of children, parents aided us in collecting survey data.

In addition to information on the daily consumption of pita bread, age, gender, body weight, and socio-economic status of the interviewed individuals, data were also collected on the level of education, address (urban or rural), job description, and health status (e.g., smoking habit, alcohol consumption and physical activity of the individual). Each questionnaire consisted of 44 questions (Supplementary Materials, Survey) which required approximately 7 min to complete. The collected answers were transmitted in the "iSurvey" application for data harmonization.

## 2.2. Chemicals

Standard solutions (1 g/L) of Hg, As, Co, Cr, Cd, Ni, and Pb were purchased from Fluka. In addition, HNO₃ 65%, H₂O₂ 30%, and the BCR (Community Bureau of Reference) certified reference material for bread trace elements (BCR-191) were purchased from Sigma-Aldrich. Deionized water, obtained from a BOECOpure UV/UF water purification system, was used in all experiments.

## 2.3. Sampling for Trace Elements Analysis

Based on the survey data, three Lebanese white pita brands were determined as the most consumed by the Lebanese population. These brands were tagged as B1, B2, and B3.

Three randomly sampled bags of each of the selected pita brands were collected. Each bread bag had seven pitas where three were randomly picked. Each pita was cut into 12 circular sectors with a plastic knife on a polyethylene plate to prevent trace element contamination. Two pieces were randomly drawn for analysis, creating 18 (3 bags × 3 pitas × 2 pieces) circular pita pieces for each brand. Before sample digestion, pita pieces were ground into powder using an agate mortar.

## 2.4. Sample Digestion for Trace Elements Analysis

Digestion of samples was done using an Anton Paar microwave (multiwave 3000—Rotor 8SXF 100, Graz, Austria). Samples of 0.5 g each were directly weighed after grinding in 50 mL Teflon tube. Then, 7 mL of concentrated HNO₃ and 1 mL of H₂O₂ were added to the preparation. The Teflon bombs were closed and placed in the digester at 200 °C for 15 min (1000 W). Once the Teflon tubes had cooled for at least 45 min, 3 mL of the digested samples were transferred to 25 mL polyethylene tubes and diluted with HNO₃ 1% before the trace elements analysis. Each sample was digested in triplicate.

## 2.5. Analysis of Trace Elements

The determination of trace elements (As, Co, Cd, Cr, Pb, and Ni) in pita samples was performed using an Atomic Absorption Spectrometer (AAS, Stafford House, UK) Zeeman Graphite Furnace GF95Z Thermo Electron corporation M series. Each pita sample was digested in triplicate and each digested sample was analyzed in triplicate. Mercury was analyzed by direct mercury analyzer DMA 80 Milestone (Sorisole, Italy) using 0.2 g of pita sample.

To guarantee the method reliability, a certified reference material (CRM; BCR-191) underwent the same digestion protocol and was analyzed at the same time as the pita samples. The recovery percentages obtained for the reference material were 93.2% for Hg, 110.8% for Pb, 86.3% for Ni, and 110.1% for Cd. For the other elements, As, Co, and Cr, for which reference materials were not available, bread samples were spiked by a known concentration (0.1 mg/L) of each of the above listed elements before digestion. The spiked bread was treated and digested in the same way as the bread
Data were collected over a two-month period, from 14 April 2017 to 5 June 2017. The team of interviewers was composed of 8 data collectors trained prior to the fieldwork by a supervisor. Participants were eligible for the study if they were Lebanese, had been living in Lebanon for at least fifteen years before the date of the survey, and did not have a chronic or serious illness.

In each supermarket, eligible participants were interviewed face-to-face by a data collector who instantly filled an offline questionnaire on an electronic tablet, using multiple-choice questions and figures (e.g., size of the pita). In the case of children, parents aided us in collecting survey data.

In addition to information on the daily consumption of pita bread, age, gender, body weight, and socio-economic status of the interviewed individuals, data were also collected on the level of education, address (urban or rural), job description, and health status (e.g., smoking habit, alcohol consumption and physical activity of the individual). Each questionnaire consisted of 44 questions (Supplementary Materials, Survey) which required approximately 7 min to complete. The collected answers were transmitted in the "iSurvey" application for data harmonization.

## 2.2. Chemicals

Standard solutions (1 g/L) of Hg, As, Co, Cr, Cd, Ni, and Pb were purchased from Fluka. In addition, HNO₃ 65%, H₂O₂ 30%, and the BCR (Community Bureau of Reference) certified reference material for bread trace elements (BCR-191) were purchased from Sigma-Aldrich. Deionized water, obtained from a BOECOpure UV/UF water purification system, was used in all experiments.

## 2.3. *Sampling for Trace Elements Analysis*

Based on the survey data, three Lebanese white pita brands were determined as the most consumed by the Lebanese population. These brands were tagged as B1, B2, and B3.

Three randomly sampled bags of each of the selected pita brands were collected. Each bread bag had seven pitas where three were randomly picked. Each pita was cut into 12 circular sectors with a plastic knife on a polyethylene plate to prevent trace element contamination. Two pieces were randomly drawn for analysis, creating 18 (3 bags × 3 pitas × 2 pieces) circular pita pieces for each brand. Before sample digestion, pita pieces were ground into powder using an agate mortar.

## 2.4. *Sample Digestion for Trace Elements Analysis*

Digestion of samples was done using an Anton Paar microwave (multiwave 3000—Rotor 8SXF 100, Graz, Austria). Samples of 0.5 g each were directly weighed after grinding in 50 mL Teflon tube. Then, 7 mL of concentrated HNO₃ and 1 mL of H₂O₂ were added to the preparation. The Teflon bombs were closed and placed in the digester at 200 °C for 15 min (1000 W). Once the Teflon tubes had cooled for at least 45 min, 3 mL of the digested samples were transferred to 25 mL polyethylene tubes and diluted with HNO₃ 1% before the trace elements analysis. Each sample was digested in triplicate.

## 2.5. *Analysis of Trace Elements*

The determination of trace elements (As, Co, Cd, Cr, Pb, and Ni) in pita samples was performed using an Atomic Absorption Spectrometer (AAS, Stafford House, UK) Zeeman Graphite Furnace GF95Z Thermo Electron corporation M series. Each pita sample was digested in triplicate and each digested sample was analyzed in triplicate. Mercury was analyzed by direct mercury analyzer DMA 80 Milestone (Sorisole, Italy) using 0.2 g of pita sample.

To guarantee the method reliability, a certified reference material (CRM; BCR-191) underwent the same digestion protocol and was analyzed at the same time as the pita samples. The recovery percentages obtained for the reference material were 93.2% for Hg, 110.8% for Pb, 86.3% for Ni, and 110.1% for Cd. For the other elements, As, Co, and Cr, for which reference materials were not available, bread samples were spiked by a known concentration (0.1 mg/L) of each of the above listed elements before digestion. The spiked bread was treated and digested in the same way as the bread
samples and analyzed using an AAS. The recovery percentages were 103.08% for As, 101.29% for Co, and 96.37% for Cr.

## 2.6. *Exposure Determination*

An exposure calculation was conducted for each respondent. They took into account the quantity of bread consumed by person, body weight, and the metal content of the consumed bread according to its brand. Results were expressed in terms of the median and 95th percentile using Ms-Excel for different population categories.

## 2.7. *Statistical Analysis*

In order to avoid a lack of representativeness of the selected samples of the surveyed population, a weighting adjustment technique was applied using auxiliary variables such as gender, age, and region by comparing the observed frequency distribution of each variable with its population distribution [29].

Data processing and statistical adjustments were analyzed using Sphinx IQ2 (Le Sphinx®, Chavanod, France) for survey data and R Studio version 1.0.153 with R version 3.5.0 [30,31].

Kruskal-Wallis non-parametric tests [32] completed by a multi-comparison Fisher's test ($\alpha$ = 0.05) were conducted to test for significant differences in consumption by region. Analysis of variances (ANOVA) completed by a multi-comparison Tukey's test ($\alpha$ = 0.05) were also executed to test differences in trace element contents among brands. They were performed with R using the "agricolae" package [33].

Median and percentile non-parametric tests ($\alpha$ = 0.05) [32] were realized to check significant differences in the level of trace element exposures between age and gender. They were performed using R and the "rcompanion" package [34].

Principal component analysis (PCA) was applied to examine the relationships among the trace element contents of the breads using XLSTAT 2018. The results are presented by loading and score plots.

Multivariate correspondence analysis was done to explore relationships among consumption behaviors and socioeconomic data [35].

# 3. Results

## 3.1. *Survey Results*

Based on a target sample of 1000 people, the surveyed sample consisted of 992 individuals with the following demographic characteristics: a sex ratio of 0.992, a proportion of 19.4% young people (8.9% children aged from 6 to 9 years old and 10.5% teenagers aged from 10 to 14 years old) and 80.6% adults, and 70.6% living in an urban area. Compared to the overall Lebanese population (where the sex ratio is 0.962, 23% are young people, 77% adults, and 87% live in urban areas [36]) there was a significant difference risk of 5% in terms of young people/adults and rural/urban ratios between the surveyed sample and the population distribution of Lebanon, which could lead to a bias. Consequently, a weighing adjustment in respect to these variables was performed to make a possible inference at the population level. The weighted numbers of respondents corresponded to 228 young people (i.e., 95 children and 133 teenagers) and 764 adults.

White pita was the type of bread most consumed according to the surveyed sample: 77% reported white pita consumption and 23% brown or other. We limited our study to the 762 individuals consuming mainly white pita (W). This population also consumed brown pita (B) and other types of pita (O) in smaller proportions. The relative proportions of these different consumption modes were as follows: W only > W + B > W + O > W + B + O with 88.8%, 8.9%, 1.8%, and 0.4%, respectively (Supplementary Materials, Table S1).

The survey makes it possible to estimate the consumption of white pita for the different categories of the population (Table 2). Adult men were the highest consumers followed by teenagers (10–14 years old) and women and children (6–9 years old): daily white pita consumption medians were 282, 143, 126,
samples and analyzed using an AAS. The recovery percentages were 103.08% for As, 101.29% for Co, and 96.37% for Cr.

## 2.6. *Exposure Determination*

An exposure calculation was conducted for each respondent. They took into account the quantity of bread consumed by person, body weight, and the metal content of the consumed bread according to its brand. Results were expressed in terms of the median and 95th percentile using Ms-Excel for different population categories.

## 2.7. Statistical Analysis

In order to avoid a lack of representativeness of the selected samples of the surveyed population, a weighting adjustment technique was applied using auxiliary variables such as gender, age, and region by comparing the observed frequency distribution of each variable with its population distribution [29].

Data processing and statistical adjustments were analyzed using Sphinx IQ2 (Le Sphinx®, Chavanod, France) for survey data and R Studio version 1.0.153 with R version 3.5.0 [30,31].

Kruskal-Wallis non-parametric tests [32] completed by a multi-comparison Fisher's test ($\alpha$ = 0.05) were conducted to test for significant differences in consumption by region. Analysis of variances (ANOVA) completed by a multi-comparison Tukey's test ($\alpha$ = 0.05) were also executed to test differences in trace element contents among brands. They were performed with R using the "agricolae" package [33].

Median and percentile non-parametric tests ($\alpha$ = 0.05) [32] were realized to check significant differences in the level of trace element exposures between age and gender. They were performed using R and the "rcompanion" package [34].

Principal component analysis (PCA) was applied to examine the relationships among the trace element contents of the breads using XLSTAT 2018. The results are presented by loading and score plots.

Multivariate correspondence analysis was done to explore relationships among consumption behaviors and socioeconomic data [35].

# 3. Results

## 3.1. Survey Results

Based on a target sample of 1000 people, the surveyed sample consisted of 992 individuals with the following demographic characteristics: a sex ratio of 0.992, a proportion of 19.4% young people (8.9% children aged from 6 to 9 years old and 10.5% teenagers aged from 10 to 14 years old) and 80.6% adults, and 70.6% living in an urban area. Compared to the overall Lebanese population (where the sex ratio is 0.962, 23% are young people, 77% adults, and 87% live in urban areas [36]) there was a significant difference risk of 5% in terms of young people/adults and rural/urban ratios between the surveyed sample and the population distribution of Lebanon, which could lead to a bias. Consequently, a weighing adjustment in respect to these variables was performed to make a possible inference at the population level. The weighted numbers of respondents corresponded to 228 young people (i.e., 95 children and 133 teenagers) and 764 adults.

White pita was the type of bread most consumed according to the surveyed sample: 77% reported white pita consumption and 23% brown or other. We limited our study to the 762 individuals consuming mainly white pita (W). This population also consumed brown pita (B) and other types of pita (O) in smaller proportions. The relative proportions of these different consumption modes were as follows: W only > W + B > W + O > W + B + O with 88.8%, 8.9%, 1.8%, and 0.4%, respectively (Supplementary Materials, Table S1).

The survey makes it possible to estimate the consumption of white pita for the different categories of the population (Table 2). Adult men were the highest consumers followed by teenagers (10–14 years old) and women and children (6–9 years old): daily white pita consumption medians were 282, 143, 126,
and 71 g/day, respectively. The 95th percentiles were 750 g/day for men and between 357 and 375 g/day for other categories. Maximum consumption was 1000 g/day for men and teenagers and 750 g/day for women and children. White pita brand consumption was significantly dependent on the age of the Lebanese population. Children consumed mostly the B2 or B3 brands, teenagers consumed mostly the B1 brands, and adults consumed mostly the B1 or B2 brands. The whole sample consumption by brand, including young people and adults, was B2 > B1 > B3 representing, respectively, 30.1%, 28.9%, and 17.1% of the total sample and the 46 other brands representing the remaining 24.9%.

<table><caption>**Table 2.** Daily white pita consumption according to population categories and repartition by brands.</caption><thead><tr><th rowspan="2">Population Categories</th><th colspan="2">Individuals</th><th colspan="6">Daily White Pita Consumption</th><th colspan="2">Brands Repartition</th></tr><tr><th>n</th><th>min</th><th>M</th><th>95th P</th><th>max</th><th>m</th><th>CI 95%</th><th>B1</th><th>B2</th><th>B3</th><th>Others *</th><th>Total</th></tr></thead><tbody><tr><td>Children</td><td>91</td><td>63</td><td>71</td><td>357</td><td>750</td><td>113</td><td>94-131</td><td>24</td><td>31</td><td>30</td><td>14</td><td>99</td></tr><tr><td>Teenagers</td><td>113</td><td>31</td><td>143</td><td>375</td><td>1000</td><td>171</td><td>147-196</td><td>40</td><td>29</td><td>12</td><td>19</td><td>100</td></tr><tr><td>Women</td><td>258</td><td>31</td><td>126</td><td>375</td><td>750</td><td>167</td><td>153-181</td><td>29</td><td>31</td><td>16</td><td>24</td><td>100</td></tr><tr><td>Men</td><td>300</td><td>31</td><td>282</td><td>750</td><td>1000</td><td>282</td><td>260-304</td><td>29</td><td>30</td><td>13</td><td>27</td><td>99</td></tr><tr><td>Young people & adults</td><td>762</td><td>31</td><td>143</td><td>500</td><td>1000</td><td>206</td><td>195-218</td><td>28.9</td><td>30.1</td><td>17.1</td><td>24.9</td><td>101</td></tr></tbody></table>

M = median; 95th P = 95th percentile; min = minimum; max = maximum; m = mean; CI 95% = 95% confidence interval. * "Others" correspond to an array of the 46 other pita brands.

White pita consumption varied according to region (Figure 2), and the highest consumption corresponded to South Lebanon and Nabatieh (355 g/day). There were no significant differences between the North, Mount Lebanon, and Bekaa regions, which were two-fold lower than South Lebanon and Nabatieh. The lowest consumption was observed in Beirut (135 g/day). This last value is in agreement with the adult population of Beirut studied by Nasreddine et al. [13].

<figure><img src="image_3.png" /><figcaption>**Figure 2.** Lebanese white pita consumption according to region. Letters represent homogeneous groups determined by a Kruskal-Wallis test and post-hoc Fisher's test for least significant difference.</figcaption></figure>

Regardless of the region, large standard deviations were observed in the distribution of bread consumption; the reason could be that adults over 65 years old eat more bread, especially in rural areas (in South Lebanon, more than 600 g per day can be consumed, see Figure 3).

In order to better explore the relationships between adult consumption and regions, a data-mining tool was used, the purpose of which was to represent in only one table or figure the most significant results from several cross analyses (based on multivariate correspondence analysis [35]). A pivot variable was first selected—the region—then other variables were crossed with it. Figure 3 shows the results of crossing region and consumption, gender, socio-professional category, smoking habits,
and 71 g/day, respectively. The 95th percentiles were 750 g/day for men and between 357 and 375 g/day for other categories. Maximum consumption was 1000 g/day for men and teenagers and 750 g/day for women and children. White pita brand consumption was significantly dependent on the age of the Lebanese population. Children consumed mostly the B2 or B3 brands, teenagers consumed mostly the B1 brands, and adults consumed mostly the B1 or B2 brands. The whole sample consumption by brand, including young people and adults, was B2 > B1 > B3 representing, respectively, 30.1%, 28.9%, and 17.1% of the total sample and the 46 other brands representing the remaining 24.9%.

<table><caption>**Table 2.** Daily white pita consumption according to population categories and repartition by brands.</caption><thead><tr><th rowspan="2">Population Categories</th><th colspan="2">Individuals</th><th colspan="6">Daily White Pita Consumption</th><th colspan="2">Brands Repartition</th></tr><tr><th>n</th><th>min</th><th>M</th><th>95th P</th><th>max</th><th>m</th><th>CI 95%</th><th>B1</th><th>B2</th><th>B3</th><th>Others *</th><th>Total</th></tr></thead><tbody><tr><td>Children</td><td>91</td><td>63</td><td>71</td><td>357</td><td>750</td><td>113</td><td>94-131</td><td>24</td><td>31</td><td>30</td><td>14</td><td>99</td></tr><tr><td>Teenagers</td><td>113</td><td>31</td><td>143</td><td>375</td><td>1000</td><td>171</td><td>147-196</td><td>40</td><td>29</td><td>12</td><td>19</td><td>100</td></tr><tr><td>Women</td><td>258</td><td>31</td><td>126</td><td>375</td><td>750</td><td>167</td><td>153-181</td><td>29</td><td>31</td><td>16</td><td>24</td><td>100</td></tr><tr><td>Men</td><td>300</td><td>31</td><td>282</td><td>750</td><td>1000</td><td>282</td><td>260-304</td><td>29</td><td>30</td><td>13</td><td>27</td><td>99</td></tr><tr><td>Young people & adults</td><td>762</td><td>31</td><td>143</td><td>500</td><td>1000</td><td>206</td><td>195-218</td><td>28.9</td><td>30.1</td><td>17.1</td><td>24.9</td><td>101</td></tr></tbody></table>

M = median; 95th P = 95th percentile; min = minimum; max = maximum; m = mean; CI 95% = 95% confidence interval. * "Others" correspond to an array of the 46 other pita brands.

White pita consumption varied according to region (Figure 2), and the highest consumption corresponded to South Lebanon and Nabatieh (355 g/day). There were no significant differences between the North, Mount Lebanon, and Bekaa regions, which were two-fold lower than South Lebanon and Nabatieh. The lowest consumption was observed in Beirut (135 g/day). This last value is in agreement with the adult population of Beirut studied by Nasreddine et al. [13].

<figure><img src="image_4.png" /><figcaption>**Figure 2.** Lebanese white pita consumption according to region. Letters represent homogeneous groups determined by a Kruskal-Wallis test and post-hoc Fisher's test for least significant difference.</figcaption></figure>

Regardless of the region, large standard deviations were observed in the distribution of bread consumption; the reason could be that adults over 65 years old eat more bread, especially in rural areas (in South Lebanon, more than 600 g per day can be consumed, see Figure 3).

In order to better explore the relationships between adult consumption and regions, a data-mining tool was used, the purpose of which was to represent in only one table or figure the most significant results from several cross analyses (based on multivariate correspondence analysis [35]). A pivot variable was first selected—the region—then other variables were crossed with it. Figure 3 shows the results of crossing region and consumption, gender, socio-professional category, smoking habits,
physical activity (4 categories: "No", "Yes, strenuous", "Yes, low", and "Yes, medium"), and area (urban or rural), using the most specific modalities as key views.

<figure><img src="image_4.png" /><figcaption>**Figure 3.** Results of crossings between regions, consumption, smoker or not, type of physical activity, socio-professional category, and area (urban or rural), using the most specific modalities as key views. Consumption is represented here by categories: less than 200 g/day (−200); from 200 to 400 g/day; from 400 to 600 g/day; and more than 600 g/day (600+). Sizes of circles are proportional to the number of observations by region.</figcaption></figure>

In South Lebanon and Nabatieh, where consumption was the highest, it correlated to rural areas with more agriculture and elderly populations. The other parts were characterized by urban areas with high concentrations of workers and retired people.

In addition, 41% of total consumers of Lebanese white pita bread were smokers (at least one cigarette per day, with a distribution of 68% men and 32% women) and 35% were alcohol drinkers (at least one glass per day, with a distribution of 67% men and 33% women, data not shown). Those behaviors may affect health and the detoxification process of alcohol and smoking are indeed well-known inducers of several xenobiotic metabolism enzymes.

Among the adult population, 85% reported engaging in physical activity which was qualified as low to medium by women and medium to high by men (Supplementary Materials Figure S1).

The correlation between consumption and health was also evaluated (data not shown). The pivot variable was still pita consumption per day and the other variables were gender, region, professional activity (student, retired, etc.), and health (smoking, physical activity).

From this study, two typologies of adult consumers can be drawn: i) Women, quite young (under 20 years old), non-smokers, and who felt very healthy according to the survey. They were
physical activity (4 categories: "No", "Yes, strenuous", "Yes, low", and "Yes, medium"), and area (urban or rural), using the most specific modalities as key views.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Results of crossings between regions, consumption, smoker or not, type of physical activity, socio-professional category, and area (urban or rural), using the most specific modalities as key views. Consumption is represented here by categories: less than 200 g/day (−200); from 200 to 400 g/day; from 400 to 600 g/day; and more than 600 g/day (600+). Sizes of circles are proportional to the number of observations by region.</figcaption></figure>

In South Lebanon and Nabatieh, where consumption was the highest, it correlated to rural areas with more agriculture and elderly populations. The other parts were characterized by urban areas with high concentrations of workers and retired people.

In addition, 41% of total consumers of Lebanese white pita bread were smokers (at least one cigarette per day, with a distribution of 68% men and 32% women) and 35% were alcohol drinkers (at least one glass per day, with a distribution of 67% men and 33% women, data not shown). Those behaviors may affect health and the detoxification process of alcohol and smoking are indeed well-known inducers of several xenobiotic metabolism enzymes.

Among the adult population, 85% reported engaging in physical activity which was qualified as low to medium by women and medium to high by men (Supplementary Materials Figure S1).

The correlation between consumption and health was also evaluated (data not shown). The pivot variable was still pita consumption per day and the other variables were gender, region, professional activity (student, retired, etc.), and health (smoking, physical activity).

From this study, two typologies of adult consumers can be drawn: i) Women, quite young (under 20 years old), non-smokers, and who felt very healthy according to the survey. They were
living in urban areas, Lebanon North-Mount Lebanon, with some level of pre-secondary education. The consumption of those aged 15-64 years old per day was approximately 163 g of pita (corresponding to 1 or 2 pitas). ii) Men between 40 and 59 years old who felt healthy. They were living in rural zones, with a technical secondary level education. Their consumption was approximately 267 g of pita per day, corresponding to 3 to 6 pitas for those 15-64 years old.

## 3.2. Trace Element Contents in Pita

The levels of trace elements in the most consumed white Lebanese pitas (B1, B2, and B3) are shown in Table 3. B1 contained the highest level of Ni (1292.1 ± 0.2 µg/kg) and Co (91 ± 3 µg/kg). B2 contained the highest level of As (400 ± 7 µg/kg), while B3 contained the highest level of Hg (0.89 ± 0.06 µg/kg), Cr (363 ± 10 µg/kg), and Pb (260 ± 81 µg/kg). Levels of Cd were lower than the maximum levels defined by the Lebanese standard (200 µg/kg of bread). In contrast, Pb exceeded the maximum level (200 µg/kg) in B2 and particularly B3 (260 ± 81 µg/kg).

<table><caption>**Table 3.** Trace element contents in white pita bread for the three major brands on the Lebanese market.</caption><thead><tr><th rowspan="2">Trace Elements</th><th colspan="3">Major Pita Brands</th><th rowspan="2">Limit of Detection (µg/kg)</th></tr><tr><th>B1</th><th>B2<br>(µg/kg Dry Weight)</th><th>B3</th></tr></thead><tbody><tr><td>Cd</td><td>&lt;LOQ</td><td>&lt;LOQ</td><td>&lt;LOQ</td><td>5</td></tr><tr><td>Hg</td><td>0.7 ± 0.1<sup>b</sup></td><td>0.70 ± 0.15<sup>b</sup></td><td>0.89 ± 0.06<sup>a</sup></td><td>0.04</td></tr><tr><td>Cr</td><td>&lt;LOQ</td><td>&lt;LOQ</td><td>363 ± 10</td><td>16</td></tr><tr><td>Co</td><td>91 ± 3<sup>a</sup></td><td>87 ± 5<sup>a</sup></td><td>84 ± 2<sup>a</sup></td><td>5</td></tr><tr><td>Ni</td><td>1292.1 ± 0.2<sup>a</sup></td><td>364.82 ± 0.02<sup>b</sup></td><td>&lt;LOQ</td><td>16</td></tr><tr><td>Pb</td><td>74 ± 5<sup>c</sup></td><td>203 ± 30<sup>b</sup></td><td>260 ± 81<sup>a</sup></td><td>21</td></tr><tr><td>As</td><td>235 ± 22<sup>c</sup></td><td>400 ± 7<sup>a</sup></td><td>321 ± 20<sup>b</sup></td><td>16</td></tr></tbody></table>

LOD = Limit of detection. Limit of quantification: LOQ = 3 × LOD. Values are the mean ± SD (n = 9/group) trace elements level. Statistical comparisons were done by Kruskal-Wallis with Tukey's post-hoc analysis. Significantly different (p < 0.05) values among brands per element are indicated by different letters.

In order to reveal the relationships among trace elements, a PCA was performed (Figure 4).

<figure><img src="image_5.png" /><figcaption>**Figure 4.** Principal component analysis showing trace element signatures of the B1, B2, and B3 pita brands.</figcaption></figure>

Principal component analysis is a useful technique for exploratory data analysis, allowing to better visualize the variation present in a dataset with many interrelated variables. In our case, the dataset
living in urban areas, Lebanon North-Mount Lebanon, with some level of pre-secondary education. The consumption of those aged 15–64 years old per day was approximately 163 g of pita (corresponding to 1 or 2 pitas). ii) Men between 40 and 59 years old who felt healthy. They were living in rural zones, with a technical secondary level education. Their consumption was approximately 267 g of pita per day, corresponding to 3 to 6 pitas for those 15–64 years old.

## 3.2. Trace Element Contents in Pita

The levels of trace elements in the most consumed white Lebanese pitas (B1, B2, and B3) are shown in Table 3. B1 contained the highest level of Ni (1292.1 ± 0.2 µg/kg) and Co (91 ± 3 µg/kg). B2 contained the highest level of As (400 ± 7 µg/kg), while B3 contained the highest level of Hg (0.89 ± 0.06 µg/kg), Cr (363 ± 10 µg/kg), and Pb (260 ± 81 µg/kg). Levels of Cd were lower than the maximum levels defined by the Lebanese standard (200 µg/kg of bread). In contrast, Pb exceeded the maximum level (200 µg/kg) in B2 and particularly B3 (260 ± 81 µg/kg).

<table><caption>**Table 3.** Trace element contents in white pita bread for the three major brands on the Lebanese market.</caption><thead><tr><th rowspan="2">Trace Elements</th><th colspan="3">Major Pita Brands</th><th rowspan="2">Limit of Detection (µg/kg)</th></tr><tr><th>B1</th><th>B2<br>(µg/kg Dry Weight)</th><th>B3</th></tr></thead><tbody><tr><td>Cd</td><td>&lt;LOQ</td><td>&lt;LOQ</td><td>&lt;LOQ</td><td>5</td></tr><tr><td>Hg</td><td>0.7 ± 0.1<sup>b</sup></td><td>0.70 ± 0.15<sup>b</sup></td><td>0.89 ± 0.06<sup>a</sup></td><td>0.04</td></tr><tr><td>Cr</td><td>&lt;LOQ</td><td>&lt;LOQ</td><td>363 ± 10</td><td>16</td></tr><tr><td>Co</td><td>91 ± 3<sup>a</sup></td><td>87 ± 5<sup>a</sup></td><td>84 ± 2<sup>a</sup></td><td>5</td></tr><tr><td>Ni</td><td>1292.1 ± 0.2<sup>a</sup></td><td>364.82 ± 0.02<sup>b</sup></td><td>&lt;LOQ</td><td>16</td></tr><tr><td>Pb</td><td>74 ± 5<sup>c</sup></td><td>203 ± 30<sup>b</sup></td><td>260 ± 81<sup>a</sup></td><td>21</td></tr><tr><td>As</td><td>235 ± 22<sup>c</sup></td><td>400 ± 7<sup>a</sup></td><td>321 ± 20<sup>b</sup></td><td>16</td></tr></tbody></table>

LOD = Limit of detection. Limit of quantification: LOQ = 3 × LOD. Values are the mean ± SD (n = 9/group) trace elements level. Statistical comparisons were done by Kruskal-Wallis with Tukey's post-hoc analysis. Significantly different (p < 0.05) values among brands per element are indicated by different letters.

In order to reveal the relationships among trace elements, a PCA was performed (Figure 4).

<figure><img src="image_6.png" /><figcaption>**Figure 4.** Principal component analysis showing trace element signatures of the B1, B2, and B3 pita brands.</figcaption></figure>

Principal component analysis is a useful technique for exploratory data analysis, allowing to better visualize the variation present in a dataset with many interrelated variables. In our case, the dataset
contained 27 observations (9 samples per brand) described by 7 variables (the trace elements content as mentioned in Table 3). Principal component analysis allows to see the overall “shape” of the data, identifying which samples are similar to one another and which are very different. It makes possible the identification of potential groups of samples that are similar and to work out which variables make one group different from another. When many variables correlate with one another, they will all contribute strongly to the same principal component [37].

Figure 4 shows typical plots from the PCA with a loading plot (on the left) and a score plot (on the right). The loading plot is used for interpreting relationships among variables and the score plot is used for interpreting relationships among observations.

All variables were well represented in the F1–F2 plot (sum of squared cosine above 0.75). The PCA shows that the following elements were correlated: Ni and Co, As and Pb as well as Cr and Hg. Lead and Co were anti-correlated, as well as Ni and Pb. There was no correlation between Cd and Co, Cd and Pb, Cr and As as well as Hg and As. The Pearson correlation matrix is provided in Supplementary Materials Table S2. We observed three separate groups: B1, B2, and B3. The dispersion in B2 and B3 was larger than in B1; this dispersion can be explained by a more homogeneous sampling for B1.

### 3.3. Trace Element Exposures

The population categories' exposure to Cd, Hg, Cr, Co, Ni, Pb, and As via white pita consumption according the most consumed brands is presented in Supplementary Materials Table S3. Using statistics, we observed that, for most of the cases, the B3 brand was the origin of the highest trace elements exposure. In regards to nickel, B1 was the main source of the highest exposure. We further decided to keep B3 TE data to calculate exposures (Table 4) using a worst-case scenario for all element levels, except for nickel for which the B1 data were retained.

<table><caption>**Table 4.** Population categories' exposure to trace elements via white pita consumption: Cd and Hg are expressed in µg/kg body weight/week while Cr, Co, Ni, Pb, and As are expressed in µg/kg body weight/day</caption><thead><tr><th rowspan="3">TEs</th><th colspan="8">Population Categories' Exposure to TEs via White Pita Consumption</th></tr><tr><th colspan="2">Children</th><th colspan="2">Teenagers</th><th colspan="2">Women</th><th colspan="2">Men</th></tr><tr><th>M <sup>1</sup></th><th>95th P <sup>1</sup></th><th>M <sup>1</sup></th><th>95th P <sup>1</sup></th><th>M <sup>1</sup></th><th>95th P <sup>1</sup></th><th>M <sup>1</sup></th><th>95th P <sup>1</sup></th></tr></thead><tbody><tr><td>Cd</td><td>0.19</td><td>1.00</td><td>0.20</td><td>0.56</td><td>0.13 <sup>a</sup></td><td>0.28 <sup>*</sup></td><td>0.19 <sup>b</sup></td><td>0.47 <sup>**</sup></td></tr><tr><td>Hg</td><td>0.03</td><td>0.15</td><td>0.03</td><td>0.09</td><td>0.02 <sup>a</sup></td><td>0.04 <sup>*</sup></td><td>0.03 <sup>b</sup></td><td>0.07 <sup>**</sup></td></tr><tr><td>Cr</td><td>1.74</td><td>9.00</td><td>1.79</td><td>5.06</td><td>2.47 <sup>a</sup></td><td>2.47 <sup>*</sup></td><td>1.72 <sup>b</sup></td><td>4.26 <sup>**</sup></td></tr><tr><td>Co</td><td>0.41</td><td>2.09</td><td>0.42</td><td>1.18</td><td>0.28 <sup>a</sup></td><td>0.57 <sup>*</sup></td><td>0.40 <sup>b</sup></td><td>0.99 <sup>**</sup></td></tr><tr><td>Ni</td><td>2.80</td><td>10.26</td><td>4.06</td><td>12.26</td><td>2.77 <sup>a</sup></td><td>10.26 <sup>*</sup></td><td>4.34 <sup>a</sup></td><td>12.26 <sup>*</sup></td></tr><tr><td>Pb</td><td>0.69</td><td>3.59</td><td>0.71</td><td>2.01</td><td>0.48 <sup>a</sup></td><td>1.77 <sup>*</sup></td><td>1.23 <sup>b</sup></td><td>3.04 <sup>**</sup></td></tr><tr><td>As</td><td>1.54</td><td>7.97</td><td>1.58</td><td>4.49</td><td>1.07 <sup>a</sup></td><td>2.19 <sup>*</sup></td><td>1.53 <sup>b</sup></td><td>3.77 <sup>**</sup></td></tr></tbody></table>

The B3 brand, the origin of the highest trace elements exposure, was chosen for comparison with toxicological reference values except for nickel for which the B1 brand was retained. ¹M = median, 95th P = 95th percentile. Statistical tests were performed separately for young people and adults for the median and the 95th percentile. Values with different letters (a,b) in each row indicate a significant difference among medians at *p* < 0.05. Values with different stars (*,**) in each row show significant difference among 95th percentiles at *p* < 0.05.

For young people, there was no significant difference among both populations in terms of exposure to any element. In contrast, we found significant differences between men and women for all elements except nickel. Men, at the 95th percentile, were always more exposed.

### 3.4. Risk Characterization (Intake Excess of TEs via Bread Consumption)

In regards to essential elements with nutritional value, risk assessment can be calculated by insufficient intake or intake excess. This work focused only on intake excess using toxicological reference points or safety levels for each element case by case in order to determine if there is a safety
contained 27 observations (9 samples per brand) described by 7 variables (the trace elements content as mentioned in Table 3). Principal component analysis allows to see the overall “shape” of the data, identifying which samples are similar to one another and which are very different. It makes possible the identification of potential groups of samples that are similar and to work out which variables make one group different from another. When many variables correlate with one another, they will all contribute strongly to the same principal component [37].

Figure 4 shows typical plots from the PCA with a loading plot (on the left) and a score plot (on the right). The loading plot is used for interpreting relationships among variables and the score plot is used for interpreting relationships among observations.

All variables were well represented in the F1–F2 plot (sum of squared cosine above 0.75). The PCA shows that the following elements were correlated: Ni and Co, As and Pb as well as Cr and Hg. Lead and Co were anti-correlated, as well as Ni and Pb. There was no correlation between Cd and Co, Cd and Pb, Cr and As as well as Hg and As. The Pearson correlation matrix is provided in Supplementary Materials Table S2. We observed three separate groups: B1, B2, and B3. The dispersion in B2 and B3 was larger than in B1; this dispersion can be explained by a more homogeneous sampling for B1.

### 3.3. Trace Element Exposures

The population categories' exposure to Cd, Hg, Cr, Co, Ni, Pb, and As via white pita consumption according the most consumed brands is presented in Supplementary Materials Table S3. Using statistics, we observed that, for most of the cases, the B3 brand was the origin of the highest trace elements exposure. In regards to nickel, B1 was the main source of the highest exposure. We further decided to keep B3 TE data to calculate exposures (Table 4) using a worst-case scenario for all element levels, except for nickel for which the B1 data were retained.

<table><caption>**Table 4.** Population categories' exposure to trace elements via white pita consumption: Cd and Hg are expressed in µg/kg body weight/week while Cr, Co, Ni, Pb, and As are expressed in µg/kg body weight/day</caption><thead><tr><th rowspan="3">TEs</th><th colspan="8">Population Categories' Exposure to TEs via White Pita Consumption</th></tr><tr><th colspan="2">Children</th><th colspan="2">Teenagers</th><th colspan="2">Women</th><th colspan="2">Men</th></tr><tr><th>M <sup>1</sup></th><th>95th P <sup>1</sup></th><th>M <sup>1</sup></th><th>95th P <sup>1</sup></th><th>M <sup>1</sup></th><th>95th P <sup>1</sup></th><th>M <sup>1</sup></th><th>95th P <sup>1</sup></th></tr></thead><tbody><tr><td>Cd</td><td>0.19</td><td>1.00</td><td>0.20</td><td>0.56</td><td>0.13 <sup>a</sup></td><td>0.28 <sup>*</sup></td><td>0.19 <sup>b</sup></td><td>0.47 <sup>**</sup></td></tr><tr><td>Hg</td><td>0.03</td><td>0.15</td><td>0.03</td><td>0.09</td><td>0.02 <sup>a</sup></td><td>0.04 <sup>*</sup></td><td>0.03 <sup>b</sup></td><td>0.07 <sup>**</sup></td></tr><tr><td>Cr</td><td>1.74</td><td>9.00</td><td>1.79</td><td>5.06</td><td>2.47 <sup>a</sup></td><td>2.47 <sup>*</sup></td><td>1.72 <sup>b</sup></td><td>4.26 <sup>**</sup></td></tr><tr><td>Co</td><td>0.41</td><td>2.09</td><td>0.42</td><td>1.18</td><td>0.28 <sup>a</sup></td><td>0.57 <sup>*</sup></td><td>0.40 <sup>b</sup></td><td>0.99 <sup>**</sup></td></tr><tr><td>Ni</td><td>2.80</td><td>10.26</td><td>4.06</td><td>12.26</td><td>2.77 <sup>a</sup></td><td>10.26 <sup>*</sup></td><td>4.34 <sup>a</sup></td><td>12.26 <sup>*</sup></td></tr><tr><td>Pb</td><td>0.69</td><td>3.59</td><td>0.71</td><td>2.01</td><td>0.48 <sup>a</sup></td><td>1.77 <sup>*</sup></td><td>1.23 <sup>b</sup></td><td>3.04 <sup>**</sup></td></tr><tr><td>As</td><td>1.54</td><td>7.97</td><td>1.58</td><td>4.49</td><td>1.07 <sup>a</sup></td><td>2.19 <sup>*</sup></td><td>1.53 <sup>b</sup></td><td>3.77 <sup>**</sup></td></tr></tbody></table>

The B3 brand, the origin of the highest trace elements exposure, was chosen for comparison with toxicological reference values except for nickel for which the B1 brand was retained. ¹M = median, 95th P = 95th percentile. Statistical tests were performed separately for young people and adults for the median and the 95th percentile. Values with different letters (a,b) in each row indicate a significant difference among medians at *p* < 0.05. Values with different stars (*,**) in each row show significant difference among 95th percentiles at *p* < 0.05.

For young people, there was no significant difference among both populations in terms of exposure to any element. In contrast, we found significant differences between men and women for all elements except nickel. Men, at the 95th percentile, were always more exposed.

### 3.4. Risk Characterization (Intake Excess of TEs via Bread Consumption)

In regards to essential elements with nutritional value, risk assessment can be calculated by insufficient intake or intake excess. This work focused only on intake excess using toxicological reference points or safety levels for each element case by case in order to determine if there is a safety
concern due to the fact of an intake excess related to bread consumption in Lebanon. Of course,
bread is not the only food consumed per day; however, in this study we show that the most exposed
population—men and teenagers—can consume up to 1000 g of bread per day.

Considering chromium levels in bread, represented mostly by Cr (III)—the major form in food [38], exposures to Cr ranged from 1.21 to 9 µg/kg bw/day which represents only 3% of the European TDI [27] for the most exposed population (children at the 95th percentile) (Table 5) in regards to pita consumption. Concerning cobalt, only the exposure related to children lead to an excess (131%) of the French TDI (Table 6) [39].

Regarding cadmium, the most exposed population were children at the 95th percentile followed by teenagers and then men and women. The most exposed population represented 40% of the European tolerable weekly intake (TWI) fixed by the EFSA [40] (Table 5).

Exposures to mercury were also very low (ranging from 0.02 to 0.15 µg/kg bw/week) even for the most exposed population (children at the 95th percentile) leading to 11.5% of the TWI fixed by the EFSA [15].

However, before concluding on safety in this survey in regards to Hg, Cd, Cr or Co, it is important
to take into account exposure from other food items in the context of the Lebanese population's
total diet.

Concerning nickel, the situation is already worrying when taking into account only pita
consumption. Teenagers were the most exposed population followed by children and then adults,
regardless of sex (Table 6). The European TDI was exceeded for all exposures for teenagers and men
and only at the 95th percentile (3 to 4 fold) for children and women.

Table 7 shows that the margins of exposure (MOE) in regard to lead, determined with a BMDL<sub>0.1</sub>
based on neurotoxic effect [41], were very low (between 0.08 and 0.58) whatever the population.
The most exposed populations in decreasing order were children, teenagers, men, and women.
The calculated MOEs suggest a safety level of preoccupation regardless of the population in terms of
lead exposure via bread consumption.

The same safety preoccupations can be observed with arsenic considering As is represented by
its inorganic form. The EFSA has concluded that 0.3 to 8 µg/kg bw/day should be used as a point of
reference [16]. Thus, these values were used for MOE calculation to assess risk characterization. Table 7
shows that the MOEs were very low regardless of the population in both scenarios. The most exposed
populations were children followed by teenagers and then men and women. The MOE data suggest a
worrying risk in regard to As exposure linked to pita consumption regardless of the population in
this survey.
concern due to the fact of an intake excess related to bread consumption in Lebanon. Of course, bread is not the only food consumed per day; however, in this study we show that the most exposed population—men and teenagers—can consume up to 1000 g of bread per day.

Considering chromium levels in bread, represented mostly by Cr (III)—the major form in food [38], exposures to Cr ranged from 1.21 to 9 µg/kg bw/day which represents only 3% of the European TDI [27] for the most exposed population (children at the 95th percentile) (Table 5) in regards to pita consumption. Concerning cobalt, only the exposure related to children lead to an excess (131%) of the French TDI (Table 6) [39].

Regarding cadmium, the most exposed population were children at the 95th percentile followed by teenagers and then men and women. The most exposed population represented 40% of the European tolerable weekly intake (TWI) fixed by the EFSA [40] (Table 5).

Exposures to mercury were also very low (ranging from 0.02 to 0.15 µg/kg bw/week) even for the most exposed population (children at the 95th percentile) leading to 11.5% of the TWI fixed by the EFSA [15].

However, before concluding on safety in this survey in regards to Hg, Cd, Cr or Co, it is important to take into account exposure from other food items in the context of the Lebanese population's total diet.

Concerning nickel, the situation is already worrying when taking into account only pita consumption. Teenagers were the most exposed population followed by children and then adults, regardless of sex (Table 6). The European TDI was exceeded for all exposures for teenagers and men and only at the 95th percentile (3 to 4 fold) for children and women.

Table 7 shows that the margins of exposure (MOE) in regard to lead, determined with a BMDL₀.₁ based on neurotoxic effect [41], were very low (between 0.08 and 0.58) whatever the population. The most exposed populations in decreasing order were children, teenagers, men, and women. The calculated MOEs suggest a safety level of preoccupation regardless of the population in terms of lead exposure via bread consumption.

The same safety preoccupations can be observed with arsenic considering As is represented by its inorganic form. The EFSA has concluded that 0.3 to 8 µg/kg bw/day should be used as a point of reference [16]. Thus, these values were used for MOE calculation to assess risk characterization. Table 7 shows that the MOEs were very low regardless of the population in both scenarios. The most exposed populations were children followed by teenagers and then men and women. The MOE data suggest a worrying risk in regard to As exposure linked to pita consumption regardless of the population in this survey.
<table><caption>**Table 5.** Population categories' exposure to Cd and Hg trace elements via B3 white pita consumption and comparison with the tolerable weekly intake (TWI). Exposure values are expressed in µg/kg body weight/week.</caption><thead><tr><th rowspan="3">Trace Elements</th><th colspan="12">Population Categories' Exposure to TEs via White Pita and Comparison with TWI</th></tr><tr><th colspan="4">Children</th><th colspan="4">Teenagers</th><th colspan="4">Women</th><th colspan="4">Men</th></tr><tr><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th></tr></thead><tbody><tr><td>Cd <sup>a</sup></td><td>0.19</td><td>7.6</td><td>1.0</td><td>40.0</td><td>0.20</td><td>8.0</td><td>0.56</td><td>22.4</td><td>0.13</td><td>5.2</td><td>0.28</td><td>11.2</td><td>0.19</td><td>7.6</td><td>0.47</td><td>18.8</td></tr><tr><td>Hg <sup>b</sup></td><td>0.03</td><td>2.3</td><td>0.15</td><td>11.5</td><td>0.03</td><td>2.4</td><td>0.09</td><td>6.7</td><td>0.02</td><td>1.5</td><td>0.04</td><td>3.1</td><td>0.03</td><td>2.3</td><td>0.07</td><td>5.6</td></tr></tbody></table>

The B3 brand, the origin of the highest trace elements exposure, was chosen for comparison with TWI values. M = median, 95thP = 95th percentile. <sup>a</sup> TWI for Cd = 2.5 µg/kg bw/week [40]. <sup>b</sup> TWI for methylmercury expressed as mercury, TWI for Hg = 1.3 µg/kg bw/week [15].

<table><caption>**Table 6.** Population categories exposure to Cr, Co, and Ni trace elements via white pita consumption and comparison with the tolerable daily intake (TDI). Exposure values are expressed in µg/kg body weight/day.</caption><thead><tr><th rowspan="3">Trace Elements</th><th colspan="12">Population Categories' Exposure to TEs via White Pita and Comparison with TDI</th></tr><tr><th colspan="4">Children</th><th colspan="4">Teenagers</th><th colspan="4">Women</th><th colspan="4">Men</th></tr><tr><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th></tr></thead><tbody><tr><td>Cr <sup>a</sup></td><td>1.74</td><td>0.58</td><td>9.0</td><td>3.0</td><td>1.79</td><td>0.6</td><td>5.06</td><td>1.69</td><td>1.21</td><td>0.40</td><td>2.47</td><td>0.82</td><td>1.72</td><td>0.57</td><td>4.26</td><td>1.42</td></tr><tr><td>Co <sup>b</sup></td><td>0.41</td><td>25.6</td><td>2.09</td><td>130.6</td><td>0.42</td><td>26.3</td><td>1.18</td><td>73.8</td><td>0.28</td><td>17.5</td><td>0.57</td><td>35.6</td><td>0.4</td><td>25.0</td><td>0.99</td><td>61.9</td></tr><tr><td>Ni <sup>c</sup></td><td>2.80</td><td>100</td><td>10.26</td><td>367</td><td>4.06</td><td>145</td><td>12.26</td><td>438</td><td>2.75</td><td>98</td><td>8.83</td><td>315</td><td>4.34</td><td>155</td><td>8.86</td><td>317</td></tr></tbody></table>

The B3 brand, the origin of the highest trace elements exposure, was chosen for comparison with TDI values, except for nickel for which the B1 brand was retained. M = median, 95thP = 95th percentile. <sup>a</sup> TDI for Cr (III) = 300 µg/kg bw/day [27]. <sup>b</sup> TDI for Co = 1.6 µg/kg bw/day [42]. <sup>c</sup> TDI for Ni = 2.8 µg/kg bw/day [14].

<table><caption>**Table 7.** Margin of exposure (MOE) to Pb and As trace elements via white pita consumption for the Lebanese population categories.</caption><thead><tr><th rowspan="3">TEs</th><th colspan="12">Population Categories' Exposure to TEs via White Pita and Comparison with the Margins of Exposure</th></tr><tr><th colspan="4">Children</th><th colspan="4">Teenagers</th><th colspan="4">Women</th><th colspan="4">Men</th></tr><tr><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th></tr></thead><tbody><tr><td>Pb <sup>a</sup></td><td>1.25</td><td>0.40</td><td>-</td><td>6.43</td><td>0.08</td><td>-</td><td>1.28</td><td>0.39</td><td>-</td><td>3.61</td><td>0.14</td><td>-</td><td>0.86</td><td>0.58</td><td>-</td><td>1.77</td><td>0.28</td><td>-</td></tr><tr><td>As <sup>b</sup></td><td>1.54</td><td>0.19</td><td>5.19</td><td>7.97</td><td>0.04</td><td>1.00</td><td>1.58</td><td>0.19</td><td>5.06</td><td>4.48</td><td>0.07</td><td>1.79</td><td>1.07</td><td>0.28</td><td>7.48</td><td>2.19</td><td>0.14</td><td>3.65</td></tr></tbody></table>

The B3 brand, the origin of the highest trace elements exposure, was chosen for the margin of exposure calculations. <sup>*</sup> Values of median and percentiles are expressed in µg/kg bw/day. <sup>a</sup> Pb margin of exposure calculated with BMDL<sub>0.1</sub> = 0.5 µg/kg bw/day based on neurotoxic effect [41]. <sup>b</sup> As margin of exposure calculated with BMDL<sub>0.1</sub> = 0.3 µg/kg bw/day for MOE 1 and BMDL<sub>0.1</sub> = 8 µg/kg bw/day for MOE 2 [16].
<table><caption>**Table 5.** Population categories' exposure to Cd and Hg trace elements via B3 white pita consumption and comparison with the tolerable weekly intake (TWI). Exposure values are expressed in µg/kg body weight/week.</caption><thead><tr><th rowspan="3">Trace Elements</th><th colspan="12">Population Categories' Exposure to TEs via White Pita and Comparison with TWI</th></tr><tr><th colspan="4">Children</th><th colspan="4">Teenagers</th><th colspan="4">Women</th><th colspan="4">Men</th></tr><tr><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th></tr></thead><tbody><tr><td>Cd <sup>a</sup></td><td>0.19</td><td>7.6</td><td>1.0</td><td>40.0</td><td>0.20</td><td>8.0</td><td>0.56</td><td>22.4</td><td>0.13</td><td>5.2</td><td>0.28</td><td>11.2</td><td>0.19</td><td>7.6</td><td>0.47</td><td>18.8</td></tr><tr><td>Hg <sup>b</sup></td><td>0.03</td><td>2.3</td><td>0.15</td><td>11.5</td><td>0.03</td><td>2.4</td><td>0.09</td><td>6.7</td><td>0.02</td><td>1.5</td><td>0.04</td><td>3.1</td><td>0.03</td><td>2.3</td><td>0.07</td><td>5.6</td></tr></tbody></table>

The B3 brand, the origin of the highest trace elements exposure, was chosen for comparison with TWI values. M = median, 95thP = 95th percentile. <sup>a</sup> TWI for Cd = 2.5 µg/kg bw/week [40]. <sup>b</sup> TWI for methylmercury expressed as mercury, TWI for Hg = 1.3 µg/kg bw/week [15].

<table><caption>**Table 6.** Population categories exposure to Cr, Co, and Ni trace elements via white pita consumption and comparison with the tolerable daily intake (TDI). Exposure values are expressed in µg/kg body weight/day.</caption><thead><tr><th rowspan="3">Trace Elements</th><th colspan="12">Population Categories' Exposure to TEs via White Pita and Comparison with TDI</th></tr><tr><th colspan="4">Children</th><th colspan="4">Teenagers</th><th colspan="4">Women</th><th colspan="4">Men</th></tr><tr><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th><th>M</th><th>M/TWI %</th><th>95thP</th><th>95thP/TWI %</th></tr></thead><tbody><tr><td>Cr <sup>a</sup></td><td>1.74</td><td>0.58</td><td>9.0</td><td>3.0</td><td>1.79</td><td>0.6</td><td>5.06</td><td>1.69</td><td>1.21</td><td>0.40</td><td>2.47</td><td>0.82</td><td>1.72</td><td>0.57</td><td>4.26</td><td>1.42</td></tr><tr><td>Co <sup>b</sup></td><td>0.41</td><td>25.6</td><td>2.09</td><td>130.6</td><td>0.42</td><td>26.3</td><td>1.18</td><td>73.8</td><td>0.28</td><td>17.5</td><td>0.57</td><td>35.6</td><td>0.4</td><td>25.0</td><td>0.99</td><td>61.9</td></tr><tr><td>Ni <sup>c</sup></td><td>2.80</td><td>100</td><td>10.26</td><td>367</td><td>4.06</td><td>145</td><td>12.26</td><td>438</td><td>2.75</td><td>98</td><td>8.83</td><td>315</td><td>4.34</td><td>155</td><td>8.86</td><td>317</td></tr></tbody></table>

The B3 brand, the origin of the highest trace elements exposure, was chosen for comparison with TDI values, except for nickel for which the B1 brand was retained. M = median, 95thP = 95th percentile. <sup>a</sup> TDI for Cr (III) = 300 µg/kg bw/day [27]. <sup>b</sup> TDI for Co = 1.6 µg/kg bw/day [42]. <sup>c</sup> TDI for Ni = 2.8 µg/kg bw/day [14].

<table><caption>**Table 7.** Margin of exposure (MOE) to Pb and As trace elements via white pita consumption for the Lebanese population categories.</caption><thead><tr><th rowspan="3">TEs</th><th colspan="12">Population Categories' Exposure to TEs via White Pita and Comparison with the Margins of Exposure</th></tr><tr><th colspan="4">Children</th><th colspan="4">Teenagers</th><th colspan="4">Women</th><th colspan="4">Men</th></tr><tr><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th><th>Value <sup>*</sup></th><th>M0E 1</th><th>M0E 2</th></tr></thead><tbody><tr><td>Pb <sup>a</sup></td><td>1.25</td><td>0.40</td><td>-</td><td>6.43</td><td>0.08</td><td>-</td><td>1.28</td><td>0.39</td><td>-</td><td>3.61</td><td>0.14</td><td>-</td><td>0.86</td><td>0.58</td><td>-</td><td>1.77</td><td>0.28</td><td>-</td></tr><tr><td>As <sup>b</sup></td><td>1.54</td><td>0.19</td><td>5.19</td><td>7.97</td><td>0.04</td><td>1.00</td><td>1.58</td><td>0.19</td><td>5.06</td><td>4.48</td><td>0.07</td><td>1.79</td><td>1.07</td><td>0.28</td><td>7.48</td><td>2.19</td><td>0.14</td><td>3.65</td></tr></tbody></table>

The B3 brand, the origin of the highest trace elements exposure, was chosen for the margin of exposure calculations. <sup>*</sup> Values of median and percentiles are expressed in µg/kg bw/day. <sup>a</sup> Pb margin of exposure calculated with BMDL<sub>0.1</sub> = 0.5 µg/kg bw/day based on neurotoxic effect [41]. <sup>b</sup> As margin of exposure calculated with BMDL<sub>0.1</sub> = 0.3 µg/kg bw/day for MOE 1 and BMDL<sub>0.1</sub> = 8 µg/kg bw/day for MOE 2 [16].
# 4. Discussion

## 4.1. Survey

The direct input of the answers to the questions asked by the data collectors on electronic tablets made it possible to calculate exposure values to TEs at the level of the individual insofar as we had all the necessary data at this scale: age, gender, body weight, daily intake, and brand of pita. The most recent works on this topic were based on a survey conducted in 2001 and for which the calculation of exposure from the intake was based on the average body weight of the whole surveyed sample [13,43].

This survey was based on a nationally representative survey, a real sample covering 992 Lebanese consumers. To our knowledge, this was the first survey assessing the pita diet in Lebanon in different Lebanese regions covering young people and adults. Previous studies focused on adults in the Beirut area only. There was a large standard deviation in the distribution of bread consumption among adults in all regions. This could be explained by the fact that people aged 65 and over eat a lot of bread.

The use of the questionnaire responses provided a good overview of pita consumption by region, age, and gender. It allowed for the drawing of typologies of people living in these regions too. The combination of estimates on consumption from the survey and the analytical results from pita experiments made it possible to do a risk evaluation.

The exposure study was done on a population consuming mainly white pita (762 individuals out of 992). The remaining population who did not consume white bread would be, in principle, more exposed to trace elements based on the article by Bou Khouzam et al. [1] which emphasizes that Lebanese brown bread is more contaminated than white.

## 4.2. Risk Characterization

In terms of risk characterization linked to the seven trace elements analyzed in the pita bread, a highly consumed foodstuff, two of the trace elements, Cd and Hg, appeared to not be a safety concern. The results for Cd are in accordance with the Nasreddine [13] study conducted on the total diet in an adult urban population. Bread may not be the main food involved in the highest exposure to Cd, as these authors showed that the main contributors to the dietary intake of Cd in a Beirut TDS were vegetables (46.8%) followed by breads and cereal-based products (30.9%). Same as for Cd, bread is not the main source of Hg, for which it is fish and seafood, which is a concern for high fish consumers [17] and, particularly, pregnant women [15]. We found between 0.70 and 0.89 µg/kg of Hg in Lebanese pita which is by far less than in the study of Mestek and co-workers [44] who found 13 µg/kg of Hg in bread in the Czech Republic.

In our study, total chromium was analyzed. The B3 pita contained the highest level of Cr (0.36 mg/kg of dry weight). It is known that the major form of Cr in food is Cr (III), and the major contributor of Cr is bread (8%) followed by bakery products. For children, it is milk (9%) followed by pasta (6%) [22]. In this study, the analytical data from the two brands (B1 and B2) were in accordance with Soares' [22] team who quantified the total chromium contents (47.3 ± 20.0 µg/kg) of dry weight white bread samples.

Compared to Cr (VI), Cr (III) is considered an essential element for humans and authorized in food supplements in Europe without any upper limit [45]. But, in 2010, the EFSA [46] recommended that total Cr does not exceed 250 µg/day, a value established by the WHO [26]. Furthermore, Cr (III) can complex with lactate and picolinate, causing concern as it has been demonstrated to be genotoxic [27]. In contrast, Cr (VI) is known to be highly toxic and carcinogenic (Group 1, [47]). It also induces skin dermatitis, then the interconversion will be of relevance in terms of risk assessment.

Even if data on the presence of Cr (VI) in food are missing, it is considered that food is a large reducing medium and that oxidation of Cr (III) to Cr (VI) would not occur. Additionally, most of the Cr (VI) is considered to be reduced in the stomach to Cr (III) [27].

Humans are exposed to Cr mainly via food, but it can be present also in water. Concerning bottled water in Lebanese markets [48], a study of the physicochemical and heavy metal parameters
# 4. Discussion

## 4.1. Survey

The direct input of the answers to the questions asked by the data collectors on electronic tablets made it possible to calculate exposure values to TEs at the level of the individual insofar as we had all the necessary data at this scale: age, gender, body weight, daily intake, and brand of pita. The most recent works on this topic were based on a survey conducted in 2001 and for which the calculation of exposure from the intake was based on the average body weight of the whole surveyed sample [13,43].

This survey was based on a nationally representative survey, a real sample covering 992 Lebanese consumers. To our knowledge, this was the first survey assessing the pita diet in Lebanon in different Lebanese regions covering young people and adults. Previous studies focused on adults in the Beirut area only. There was a large standard deviation in the distribution of bread consumption among adults in all regions. This could be explained by the fact that people aged 65 and over eat a lot of bread.

The use of the questionnaire responses provided a good overview of pita consumption by region, age, and gender. It allowed for the drawing of typologies of people living in these regions too. The combination of estimates on consumption from the survey and the analytical results from pita experiments made it possible to do a risk evaluation.

The exposure study was done on a population consuming mainly white pita (762 individuals out of 992). The remaining population who did not consume white bread would be, in principle, more exposed to trace elements based on the article by Bou Khouzam et al. [1] which emphasizes that Lebanese brown bread is more contaminated than white.

## 4.2. Risk Characterization

In terms of risk characterization linked to the seven trace elements analyzed in the pita bread, a highly consumed foodstuff, two of the trace elements, Cd and Hg, appeared to not be a safety concern. The results for Cd are in accordance with the Nasreddine [13] study conducted on the total diet in an adult urban population. Bread may not be the main food involved in the highest exposure to Cd, as these authors showed that the main contributors to the dietary intake of Cd in a Beirut TDS were vegetables (46.8%) followed by breads and cereal-based products (30.9%). Same as for Cd, bread is not the main source of Hg, for which it is fish and seafood, which is a concern for high fish consumers [17] and, particularly, pregnant women [15]. We found between 0.70 and 0.89 µg/kg of Hg in Lebanese pita which is by far less than in the study of Mestek and co-workers [44] who found 13 µg/kg of Hg in bread in the Czech Republic.

In our study, total chromium was analyzed. The B3 pita contained the highest level of Cr (0.36 mg/kg of dry weight). It is known that the major form of Cr in food is Cr (III), and the major contributor of Cr is bread (8%) followed by bakery products. For children, it is milk (9%) followed by pasta (6%) [22]. In this study, the analytical data from the two brands (B1 and B2) were in accordance with Soares' [22] team who quantified the total chromium contents (47.3 ± 20.0 µg/kg) of dry weight white bread samples.

Compared to Cr (VI), Cr (III) is considered an essential element for humans and authorized in food supplements in Europe without any upper limit [45]. But, in 2010, the EFSA [46] recommended that total Cr does not exceed 250 µg/day, a value established by the WHO [26]. Furthermore, Cr (III) can complex with lactate and picolinate, causing concern as it has been demonstrated to be genotoxic [27]. In contrast, Cr (VI) is known to be highly toxic and carcinogenic (Group 1, [47]). It also induces skin dermatitis, then the interconversion will be of relevance in terms of risk assessment.

Even if data on the presence of Cr (VI) in food are missing, it is considered that food is a large reducing medium and that oxidation of Cr (III) to Cr (VI) would not occur. Additionally, most of the Cr (VI) is considered to be reduced in the stomach to Cr (III) [27].

Humans are exposed to Cr mainly via food, but it can be present also in water. Concerning bottled water in Lebanese markets [48], a study of the physicochemical and heavy metal parameters
in 32 samples of bottled waters, 25 of these samples showed values below 1 µg/L for Cr, while the other samples showed values ranging between 1 and 4 µg/L in the drinking water. Although naturally occurring Cr (VI) is rarely found in the environment, its presence in ground and surface waters, as a consequence of oxidation of Cr (III)-bearing minerals by Mn (IV) oxides, has been reported recently [49]. Furthermore, in terms of the total quantity found in bread, we do not know its respective part coming from the wheat flour or water used during the bread-making process. Wheat are among species that demonstrate the ability to bio-accumulate chromium the most, just after tomato [5]. An average value of Cr below 0.05 mg/kg of dry weight was measured in wheat seeds in Swedish long-term soil-fertility experiments between 1967 and 2003 [50].

The recipient might also be a source of contamination if Cr is released. Indeed, migration of chromium from steel material occurred in the past mostly from Asia (China) and was above the upper legal limits for food contact material [51].

The EFSA panel recommended the generation of data using sensitive analytical methodologies which specifically measure the content of Cr (III) and Cr (VI) in food and drinking water in different EU Member States [27]; the issue is the same for Lebanon. However, Cr speciation is difficult to conduct, especially in complex matrices. Chromium speciation depends on several factors such as pH, concentration, oxidation state, and complexation with natural components present in the matrix such as peptide [49]. Soares and co-workers [22] mentioned the presence of Cr (VI) in bread after extraction of Cr by alkaline solutions. However, Novotnik et al. [38] demonstrated that the only Cr species in bread is Cr (III) by spiking with enriched isotopes of ⁵³Cr (III) and ⁵⁰Cr (VI), a more accurate quantification, and using HPLC-ICP-MS. Thus, we made a characterization in regards to Cr (III) and, in this way, there were no safety concerns in this study for Cr (III) in Lebanese pita bread.

For cobalt, the exposure of young people under 15 years old at the 95th percentile exceeded the TDI by a percentage of 131%. Bread and cereals are significant source of Co exposure (18.1%), with a daily intake of cobalt estimated to be around 11.4 µg/day according to the Nasreddine study [13]. Indeed, in their study they found 6.36 µg/kg of cobalt in bread and cereals. In this study, the level of Co was 91 µg/kg for the B1 pita brand, which is 7 fold more than in Australian bread and 13.5 fold higher than the data from the Nasreddine study [13] but 1.8 fold less than in Ethiopian bread [52].

In Lebanon, imported flour samples used in Lebanese white pita presented values for cobalt between 2 ± 2 and 33 ± 7 µg/kg which is far less than in the pita. Cobalt contamination can also occur during bread preparation via the water used. The mean level of Co in tap water is approximately 9 µg/L (IRAL, personal communication). In 2013, the Environmental Working Group's (EWG) mission detected Co in only one sample (6 analyzed) and found 1.01 µg/L [53]. Furthermore, cobalt can also be released from food contact materials (FCMs) as it is used in ceramics [54]; thus, we can wonder if cobalt migration can occur from recipients used for making bread in Lebanon.

Several forms of cobalt (II) (sulphate, dichloride) are classified as presumed human carcinogens via the inhalatory route [55]. The IARC classified the soluble form of cobalt II in class 2B but as metal, Co was classified in 2A [56]. Although Co (II) salts are able to induce genotoxicity in vitro and in vivo after oral or parenteral exposure; however, there is a lack of carcinogenesis studies in humans and animals following an oral route [57]. The EFSA concluded that it cannot be excluded that cobalt could have a non-threshold toxicological effect [58]. In France, ANSES made the recommendation to reduce food exposure to cobalt in infants under 3 years old [39]. Thus, for cobalt, the daily intake from all dietary sources may be quite pertinent to monitor in order to determine the total exposure of young people under 15 years old in Lebanon.

In regards to nickel, the B1 brand contained the highest value (1292 µg/kg), which was above values found in studies for Nigerian, Ethiopian, and Spanish bread [52,59,60] but under those obtained in Iranian bread [61]. However, for another brand the value (365 µg/kg) was in the same range.

The average daily intake of nickel was estimated at 126.27 µg/day bread and cereals contributing to a 54.73 µg/day Ni intake in the Lebanese adult urban population [13].
in 32 samples of bottled waters, 25 of these samples showed values below 1 µg/L for Cr, while the other samples showed values ranging between 1 and 4 µg/L in the drinking water. Although naturally occurring Cr (VI) is rarely found in the environment, its presence in ground and surface waters, as a consequence of oxidation of Cr (III)-bearing minerals by Mn (IV) oxides, has been reported recently [49]. Furthermore, in terms of the total quantity found in bread, we do not know its respective part coming from the wheat flour or water used during the bread-making process. Wheat are among species that demonstrate the ability to bio-accumulate chromium the most, just after tomato [5]. An average value of Cr below 0.05 mg/kg of dry weight was measured in wheat seeds in Swedish long-term soil-fertility experiments between 1967 and 2003 [50].

The recipient might also be a source of contamination if Cr is released. Indeed, migration of chromium from steel material occurred in the past mostly from Asia (China) and was above the upper legal limits for food contact material [51].

The EFSA panel recommended the generation of data using sensitive analytical methodologies which specifically measure the content of Cr (III) and Cr (VI) in food and drinking water in different EU Member States [27]; the issue is the same for Lebanon. However, Cr speciation is difficult to conduct, especially in complex matrices. Chromium speciation depends on several factors such as pH, concentration, oxidation state, and complexation with natural components present in the matrix such as peptide [49]. Soares and co-workers [22] mentioned the presence of Cr (VI) in bread after extraction of Cr by alkaline solutions. However, Novotnik et al. [38] demonstrated that the only Cr species in bread is Cr (III) by spiking with enriched isotopes of ⁵³Cr (III) and ⁵⁰Cr (VI), a more accurate quantification, and using HPLC-ICP-MS. Thus, we made a characterization in regards to Cr (III) and, in this way, there were no safety concerns in this study for Cr (III) in Lebanese pita bread.

For cobalt, the exposure of young people under 15 years old at the 95th percentile exceeded the TDI by a percentage of 131%. Bread and cereals are significant source of Co exposure (18.1%), with a daily intake of cobalt estimated to be around 11.4 µg/day according to the Nasreddine study [13]. Indeed, in their study they found 6.36 µg/kg of cobalt in bread and cereals. In this study, the level of Co was 91 µg/kg for the B1 pita brand, which is 7 fold more than in Australian bread and 13.5 fold higher than the data from the Nasreddine study [13] but 1.8 fold less than in Ethiopian bread [52].

In Lebanon, imported flour samples used in Lebanese white pita presented values for cobalt between 2 ± 2 and 33 ± 7 µg/kg which is far less than in the pita. Cobalt contamination can also occur during bread preparation via the water used. The mean level of Co in tap water is approximately 9 µg/L (IRAL, personal communication). In 2013, the Environmental Working Group's (EWG) mission detected Co in only one sample (6 analyzed) and found 1.01 µg/L [53]. Furthermore, cobalt can also be released from food contact materials (FCMs) as it is used in ceramics [54]; thus, we can wonder if cobalt migration can occur from recipients used for making bread in Lebanon.

Several forms of cobalt (II) (sulphate, dichloride) are classified as presumed human carcinogens via the inhalatory route [55]. The IARC classified the soluble form of cobalt II in class 2B but as metal, Co was classified in 2A [56]. Although Co (II) salts are able to induce genotoxicity in vitro and in vivo after oral or parenteral exposure; however, there is a lack of carcinogenesis studies in humans and animals following an oral route [57]. The EFSA concluded that it cannot be excluded that cobalt could have a non-threshold toxicological effect [58]. In France, ANSES made the recommendation to reduce food exposure to cobalt in infants under 3 years old [39]. Thus, for cobalt, the daily intake from all dietary sources may be quite pertinent to monitor in order to determine the total exposure of young people under 15 years old in Lebanon.

In regards to nickel, the B1 brand contained the highest value (1292 µg/kg), which was above values found in studies for Nigerian, Ethiopian, and Spanish bread [52,59,60] but under those obtained in Iranian bread [61]. However, for another brand the value (365 µg/kg) was in the same range.

The average daily intake of nickel was estimated at 126.27 µg/day bread and cereals contributing to a 54.73 µg/day Ni intake in the Lebanese adult urban population [13].
In a French TDS study regarding the infant population, the value determined in bread was approximately 69.5 µg/kg [39]. The contribution of drinking water to total exposure to nickel was very small across dietary surveys and age classes [14]. We do not have an explanation, but this high level could be due to the recipient used by bakers. Using the TDI 2.8 µg/kg bw/day [14] for oral chronic exposure, the data demonstrate that the percentage of Ni TDI was exceeded 3 to 4 fold for all populations with the 10–14 years old range being the most exposed. These data are in accordance with the TDS studies also showing a safety concern in France for very young people (<3 years old) [39].

Nickel compounds and metallic nickel are classified as human carcinogens, respectively, group 1 and group 2B after inhalation [62]. Available studies in animals suggest that exposure to nickel salts lead to renal effects and increases neonatal mortality with the kidney being the principal target. Reproductive and developmental toxicity were the critical effects after chronic oral exposure retained for risk characterization. However, oral exposure to nickel salts could be at the origin of contact dermatitis for sensitive populations. Indeed, oral nickel can also cause an oral and chronic effect in people dermally sensitized [14].

In humans, absorption of Ni varies according to the method of exposure (in water it is 40 fold more absorbed compared to food; 0.7% and the bioavailability varies with the food type [63]). There are no maximum levels in food for nickel; in contrast, a value of 20 µg/L is fixed in water intended for human consumption and in mineral waters in Europe [64,65].

In this study, arsenic exposure was more worrying, as we found a high level of arsenic for all brand analyzed (235–400 µg/kg) which is around 63 fold higher than those found in a bread study in the UK [66] and 32 fold higher than those mentioned in a French infant TDS study [39]. Bou Khouzam et al. [1] mentioned seasonal variation in the arsenic concentration in bread (33 µg/kg in dry season, which is 4 fold higher than in wet season). Arsenic was not studied by Nasreddine et al. [13] and calculated MOEs were very low for all bakeries and the population. These data are in accordance with the data and risk characterization in 19 countries in Europe [16], for example, in France where there is a safety concern with MOEs lower than 36 [39], suggesting that the situation in Lebanon is worrying.

The inorganic form of As, which is known to be toxic, may cause human malignancies. It was classified as carcinogenic (group 1) in Reference [67]. The toxicity, due to the inorganic form, occurs by an epigenetic mechanism, and there is increasing evidence that early life exposure affecting fetal and infant growth [68] may cause chronic disease later in life [69]. This highlights the urgent need to obtain data on the effect of the toxicity of As on development [39]. Arsenic passes through the placental barrier in humans [70], resulting in similar levels in both the fetus and mother. Both the inorganic form and its methylated metabolites are transferred [71]. The effect of As on the endocrine system has also been investigated for many years showing a disruption on the gonadal, adrenal, and thyroid endocrine systems [72]. Furthermore, it has been shown a toxic effect due to an interaction between Cd and As [73,74]. Drinking water, crops irrigated with contaminated water, and food prepared with contaminated water are sources of exposure [75].

The EFSA's recommendation is to reduce inorganic arsenic exposure in food and develop robust validated methods to determine the inorganic form [16]. We do not distinguish between organic and inorganic form of As in this study. It is known that arsenic intake is higher from solid foods than from liquids.

Furthermore, speciation plays a major role in determining the amount of arsenic absorbed after consumption of As contaminated food as well as bioavailability with significant interspecies differences for organic arsenic. Regards to inorganic form, absorption can be high as demonstrated with contaminated rice [76,77].

Representative data on speciation are indeed scarce. The EFSA panel assumed a proportion of the contaminant to inorganic form (the toxic one) to vary between 50 (best case) to 100% (worst case) of the total arsenic reported in food commodities in Europe, other than fish and seafood [16]. Furthermore, in bread, 92% of the inorganic form has been found [78]. Cereals and cereal-based products were identified as largely contributing to the daily exposure to inorganic arsenic, as it was quantified as 53 µg
In a French TDS study regarding the infant population, the value determined in bread was approximately 69.5 µg/kg [39]. The contribution of drinking water to total exposure to nickel was very small across dietary surveys and age classes [14]. We do not have an explanation, but this high level could be due to the recipient used by bakers. Using the TDI 2.8 µg/kg bw/day [14] for oral chronic exposure, the data demonstrate that the percentage of Ni TDI was exceeded 3 to 4 fold for all populations with the 10–14 years old range being the most exposed. These data are in accordance with the TDS studies also showing a safety concern in France for very young people (<3 years old) [39].

Nickel compounds and metallic nickel are classified as human carcinogens, respectively, group 1 and group 2B after inhalation [62]. Available studies in animals suggest that exposure to nickel salts lead to renal effects and increases neonatal mortality with the kidney being the principal target. Reproductive and developmental toxicity were the critical effects after chronic oral exposure retained for risk characterization. However, oral exposure to nickel salts could be at the origin of contact dermatitis for sensitive populations. Indeed, oral nickel can also cause an oral and chronic effect in people dermally sensitized [14].

In humans, absorption of Ni varies according to the method of exposure (in water it is 40 fold more absorbed compared to food; 0.7% and the bioavailability varies with the food type [63]). There are no maximum levels in food for nickel; in contrast, a value of 20 µg/L is fixed in water intended for human consumption and in mineral waters in Europe [64,65].

In this study, arsenic exposure was more worrying, as we found a high level of arsenic for all brand analyzed (235–400 µg/kg) which is around 63 fold higher than those found in a bread study in the UK [66] and 32 fold higher than those mentioned in a French infant TDS study [39]. Bou Khouzam et al. [1] mentioned seasonal variation in the arsenic concentration in bread (33 µg/kg in dry season, which is 4 fold higher than in wet season). Arsenic was not studied by Nasreddine et al. [13] and calculated MOEs were very low for all bakeries and the population. These data are in accordance with the data and risk characterization in 19 countries in Europe [16], for example, in France where there is a safety concern with MOEs lower than 36 [39], suggesting that the situation in Lebanon is worrying.

The inorganic form of As, which is known to be toxic, may cause human malignancies. It was classified as carcinogenic (group 1) in Reference [67]. The toxicity, due to the inorganic form, occurs by an epigenetic mechanism, and there is increasing evidence that early life exposure affecting fetal and infant growth [68] may cause chronic disease later in life [69]. This highlights the urgent need to obtain data on the effect of the toxicity of As on development [39]. Arsenic passes through the placental barrier in humans [70], resulting in similar levels in both the fetus and mother. Both the inorganic form and its methylated metabolites are transferred [71]. The effect of As on the endocrine system has also been investigated for many years showing a disruption on the gonadal, adrenal, and thyroid endocrine systems [72]. Furthermore, it has been shown a toxic effect due to an interaction between Cd and As [73,74]. Drinking water, crops irrigated with contaminated water, and food prepared with contaminated water are sources of exposure [75].

The EFSA's recommendation is to reduce inorganic arsenic exposure in food and develop robust validated methods to determine the inorganic form [16]. We do not distinguish between organic and inorganic form of As in this study. It is known that arsenic intake is higher from solid foods than from liquids.

Furthermore, speciation plays a major role in determining the amount of arsenic absorbed after consumption of As contaminated food as well as bioavailability with significant interspecies differences for organic arsenic. Regards to inorganic form, absorption can be high as demonstrated with contaminated rice [76,77].

Representative data on speciation are indeed scarce. The EFSA panel assumed a proportion of the contaminant to inorganic form (the toxic one) to vary between 50 (best case) to 100% (worst case) of the total arsenic reported in food commodities in Europe, other than fish and seafood [16]. Furthermore, in bread, 92% of the inorganic form has been found [78]. Cereals and cereal-based products were identified as largely contributing to the daily exposure to inorganic arsenic, as it was quantified as 53 µg
up to 72 µg of arsenic per kilo for all cereals and cereal product categories. In wheat samples, it was demonstrated that 91% to 95% of As was in the inorganic form [77,78]. In wheat plants cultivated in the Bekaa region, the element As presented values between 288.3 ± 14.7 and 743.7 ± 25.4 mg/kg. In the irrigation water used in the Bekaa sites cultivated for wheat plants, the As content varied between 0.8 ± 0.3 µg/L and 1.2 ± 0.2 µg/L (IRAL data). Indeed, the level of As in drinking water in Lebanon varied between 1.2 and 4.5 µg/L [79]. Water used for making bread can be the main contributor to the level of arsenic in the pita with the main form recovered in drinking water as inorganic. It is important to note that a change in the specification of As can occur during food preparation [80].

Concerning the level of lead in the different brands: B2 and B3 contents were more than 2.5 fold higher than B1 (74, 203, and 260 µg/kg, respectively). The order of magnitude of B1 is in accordance with Turkish, Iranian, Ethiopian, and Spanish studies on bread [52,59] and 7 fold more higher than in data on the Czech Republic [44]. The B3 pita lead content was 32 fold higher than Lebanese data [13] which focused on Beirut. The authors noticed that the lead exposure in Beirut was even two fold lower compared to a previous study [43]. Thus, our data were not in accordance, but it took into account all the Lebanese regions including rural ones. The lead origin may be from wheat. Several samples of Lebanese wheat cultivated in the Bekaa region were studied, showing values of lead varying between 4.5 and 23 µg/kg but far less than in the bread, noting that the lead values in Lebanese wheat samples were lower than the values showing in wheat samples in Saudi Arabian markets (2.81 ± 0.08 mg/kg dry weight) [12]. While, the lead levels, measured in five wheat flour samples in Lebanese markets varied between 0.0045 ± 0.002 and 0.023 ± 0.016 mg/kg dry weight, they were lower than values presented in wheat flour samples in the Canary Islands and used for bread production (varying between 0.037 ± 0.013 and 0.056 ± 0.045 mg/kg fresh weight) [81].

Furthermore, lead can also be a contaminant issued from recipient food contact materials or water used for bread making.

Human data show that the developing nervous system is a critical target in children. Various water soluble and insoluble lead can induce tumors in rodents, especially renal tumors as a carcinogenic/promoter; brain gliomas were also observed. Classified lead compounds are probably carcinogenic to humans on the basis of limited evidence of carcinogenicity in humans and sufficient evidence in animals [56].

It is important to note that due to the long lead half-life, chronic exposure is of most concern and children absorb lead after ingestion to a much greater extent (70% of the amount of lead ingested) than adults (20%) [82].

All MOEs calculated with exposure only via pita consumption in this study were very low, but in the same range than in France and Europe [17,39,41], suggesting an exposure safety concern and that it is necessary to decrease Pb exposure, for example, by varying the diet due to the difficulty in reducing Pb contamination.

For both lead and arsenic, this contamination occurred for all bread brands analyzed in this study and there is a safety concern (Supplementary Materials Table S3).

**5. Conclusions**

This was the first study conducted all over Lebanon showing the exposure of different categories of Lebanese population to trace elements via bread. A survey of daily consumption of pita bread was conducted on a sample of 992 people. Among the four categories (i.e., children, teenagers, women, and men), white bread was the most consumed. Men consumed the highest quantity (282 g/day). The survey revealed that three brands of white pita were the most consumed: B1, B2, and B3.

Seven trace elements (i.e., As, Cd, Co, Cr, Hg, Ni, and Pb) were quantified in pitas coming from B1, B2, and B3. B1 was characterized by high levels of Ni and Co. B2 was characterized by the presence of As and Cd, and B3 was characterized by the presence of Cr, Hg, and Pb. Lead concentrations exceeded the authorized limits of the Lebanese standard for bread in the two brands. Further analysis should be
up to 72 µg of arsenic per kilo for all cereals and cereal product categories. In wheat samples, it was demonstrated that 91% to 95% of As was in the inorganic form [77,78]. In wheat plants cultivated in the Bekaa region, the element As presented values between 288.3 ± 14.7 and 743.7 ± 25.4 mg/kg. In the irrigation water used in the Bekaa sites cultivated for wheat plants, the As content varied between 0.8 ± 0.3 µg/L and 1.2 ± 0.2 µg/L (IRAL data). Indeed, the level of As in drinking water in Lebanon varied between 1.2 and 4.5 µg/L [79]. Water used for making bread can be the main contributor to the level of arsenic in the pita with the main form recovered in drinking water as inorganic. It is important to note that a change in the specification of As can occur during food preparation [80].

Concerning the level of lead in the different brands: B2 and B3 contents were more than 2.5 fold higher than B1 (74, 203, and 260 µg/kg, respectively). The order of magnitude of B1 is in accordance with Turkish, Iranian, Ethiopian, and Spanish studies on bread [52,59] and 7 fold more higher than in data on the Czech Republic [44]. The B3 pita lead content was 32 fold higher than Lebanese data [13] which focused on Beirut. The authors noticed that the lead exposure in Beirut was even two fold lower compared to a previous study [43]. Thus, our data were not in accordance, but it took into account all the Lebanese regions including rural ones. The lead origin may be from wheat. Several samples of Lebanese wheat cultivated in the Bekaa region were studied, showing values of lead varying between 4.5 and 23 µg/kg but far less than in the bread, noting that the lead values in Lebanese wheat samples were lower than the values showing in wheat samples in Saudi Arabian markets (2.81 ± 0.08 mg/kg dry weight) [12]. While, the lead levels, measured in five wheat flour samples in Lebanese markets varied between 0.0045 ± 0.002 and 0.023 ± 0.016 mg/kg dry weight, they were lower than values presented in wheat flour samples in the Canary Islands and used for bread production (varying between 0.037 ± 0.013 and 0.056 ± 0.045 mg/kg fresh weight) [81].

Furthermore, lead can also be a contaminant issued from recipient food contact materials or water used for bread making.

Human data show that the developing nervous system is a critical target in children. Various water soluble and insoluble lead can induce tumors in rodents, especially renal tumors as a carcinogenic/promoter; brain gliomas were also observed. Classified lead compounds are probably carcinogenic to humans on the basis of limited evidence of carcinogenicity in humans and sufficient evidence in animals [56].

It is important to note that due to the long lead half-life, chronic exposure is of most concern and children absorb lead after ingestion to a much greater extent (70% of the amount of lead ingested) than adults (20%) [82].

All MOEs calculated with exposure only via pita consumption in this study were very low, but in the same range than in France and Europe [17,39,41], suggesting an exposure safety concern and that it is necessary to decrease Pb exposure, for example, by varying the diet due to the difficulty in reducing Pb contamination.

For both lead and arsenic, this contamination occurred for all bread brands analyzed in this study and there is a safety concern (Supplementary Materials Table S3).

## 5. Conclusions

This was the first study conducted all over Lebanon showing the exposure of different categories of Lebanese population to trace elements via bread. A survey of daily consumption of pita bread was conducted on a sample of 992 people. Among the four categories (i.e., children, teenagers, women, and men), white bread was the most consumed. Men consumed the highest quantity (282 g/day). The survey revealed that three brands of white pita were the most consumed: B1, B2, and B3.

Seven trace elements (i.e., As, Cd, Co, Cr, Hg, Ni, and Pb) were quantified in pitas coming from B1, B2, and B3. B1 was characterized by high levels of Ni and Co. B2 was characterized by the presence of As and Cd, and B3 was characterized by the presence of Cr, Hg, and Pb. Lead concentrations exceeded the authorized limits of the Lebanese standard for bread in the two brands. Further analysis should be
performed on soil, wheat, flour, and water in order to determine the source of the contamination and to reduce exposure.

The B3 brand was at the origin of the highest trace element exposure, except for nickel for which B1 was retained. In terms of risk characterization, there was no safety concern due to the fact of pita consumption for all studied categories with Hg, Cd, Cr, and for Co except the 95th percentile of the children category (6–9 years old). Nickel exposure showed an excess of the TDI up to 4 fold for the most exposed populations (95th percentile). There were safety concerns related to As and Pb, as the margins of exposure were very low. To our knowledge, this is the first time that exposure to As was studied in the Lebanese population. To confirm these observed data, it is urgent to analyze the inorganic arsenic form levels in bread. Experiments are currently being run to study any cocktail effects (synergistic) using pertinent bioassays due to the exposure to the sum of trace element identified in bread. The ultimate objectives are to protect the Lebanese population by reducing exposures to populations particularly exposed to TE presenting safety concern and to make recommendations, as bread is only one part of the daily diet.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2304-8158/8/11/574/s1, Survey and questionnaire related to Lebanese bread consumption. Figure S1: Distribution of physical activity intensity by gender in adult population. Table S1-a: Cross-tabulation of consumed pita type occurrences for the total 992 respondents in the number of individuals. Table S1-b: Types of consumed bread for the subpopulation consuming mainly white pita (762 individuals). Table S2: Pearson correlation matrix of the principal component analysis. Table S3: Population categories' exposure to Cd, Hg, Cr, Co, Ni, Pb, and As trace elements via white pita consumption according to the most consumed brands.

**Author Contributions:** Conceptualization, C.D., R.O., H.C., P.C., E.B.-M. and M.-C.C.; Formal analysis, N.L., L.D. and E.B.-M.; Funding acquisition, C.D., M.A. and P.C.; Investigation, N.L.; Methodology, N.L., P.C., L.D., E.B.-M. and M.-C.C.; Project administration, M.A., P.C. and M.-C.C.; Supervision, C.D., R.O., H.C., M.A., P.C. and E.B.-M.; Validation, N.L. and E.B.-M.; Writing—original draft, N.L.; Writing—review and editing, C.D., R.O., H.C., P.C., L.D., E.B.-M. and M.-C.C.

**Funding:** This research was funded by the project PHC-Cèdre (bilateral French–Lebanese funding program) and the Lebanese Agricultural Research Institute (LARI).

**Acknowledgments:** The authors warmly thank and express their gratitude to Jacqueline Saad Harfouche for her help in the survey implementation.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Bou Khouzam, R.; Pohl, P.; Al Ayoubi, B.; Jaber, F.; Lobinski, R. Concentrations of toxic and essential elements in Lebanese bread. *Pure Appl. Chem.* **2012**, *84*, 181–190. [CrossRef]

2. Marín, S.; Pardo, O.; Sánchez, A.; Sanchis, Y.; Vélez, D.; Devesa, V.; Font, G.; Yusà, V. Assessment of metal levels in foodstuffs from the Region of Valencia (Spain). *Toxicol. Rep.* **2018**, *5*, 654–670. [CrossRef] [PubMed]

3. Chaza, C.; Rayane, S.; Sopheak, N.; Moomen, B.; Baghdad, O. Distribution of Organochlorine Pesticides and Heavy Metals in Lebanese Agricultural Soil: Case Study—Plain of Akkar. *Int. J. Environ. Res.* **2018**, *12*, 631–649. [CrossRef]

4. Baize, D.; Bellanger, L.; Tomassone, R. Relationships between concentrations of trace metals in wheat grains and soil. *Agron. Sustain. Dev.* **2009**, *29*, 297–312. [CrossRef]

5. Gupta, N.; Yadav, K.K.; Kumar, V.; Kumar, S.; Chadd, R.P.; Kumar, A. Trace elements in soil-vegetables interface: Translocation, bioaccumulation, toxicity and amelioration—A review. *Sci. Total Environ.* **2019**, *651*, 2927–2942. [CrossRef]

6. Mitra, G. Essential Plant Nutrients and Recent Concepts about their Uptake. In *Essential Plant Nutrients: Uptake, Use Efficiency, and Management*; Naeem, M., Ansari, A.A., Gill, S.S., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 3–36. ISBN 978-3-319-58841-4.

7. Counter, S.A.; Buchanan, L.H. Mercury exposure in children: A review. *Toxicol. Appl. Pharmacol.* **2004**, *198*, 209–230. [CrossRef]
performed on soil, wheat, flour, and water in order to determine the source of the contamination and to reduce exposure.

The B3 brand was at the origin of the highest trace element exposure, except for nickel for which B1 was retained. In terms of risk characterization, there was no safety concern due to the fact of pita consumption for all studied categories with Hg, Cd, Cr, and for Co except the 95th percentile of the children category (6–9 years old). Nickel exposure showed an excess of the TDI up to 4 fold for the most exposed populations (95th percentile). There were safety concerns related to As and Pb, as the margins of exposure were very low. To our knowledge, this is the first time that exposure to As was studied in the Lebanese population. To confirm these observed data, it is urgent to analyze the inorganic arsenic form levels in bread. Experiments are currently being run to study any cocktail effects (synergistic) using pertinent bioassays due to the exposure to the sum of trace element identified in bread. The ultimate objectives are to protect the Lebanese population by reducing exposures to populations particularly exposed to TE presenting safety concern and to make recommendations, as bread is only one part of the daily diet.

**Supplementary Materials:** The following are available online at http://www.mdpi.com/2304-8158/8/11/574/s1, Survey and questionnaire related to Lebanese bread consumption. Figure S1: Distribution of physical activity intensity by gender in adult population. Table S1-a: Cross-tabulation of consumed pita type occurrences for the total 992 respondents in the number of individuals. Table S1-b: Types of consumed bread for the subpopulation consuming mainly white pita (762 individuals). Table S2: Pearson correlation matrix of the principal component analysis. Table S3: Population categories' exposure to Cd, Hg, Cr, Co, Ni, Pb, and As trace elements via white pita consumption according to the most consumed brands.

**Author Contributions:** Conceptualization, C.D., R.O., H.C., P.C., E.B.-M. and M.-C.C.; Formal analysis, N.L., L.D. and E.B.-M.; Funding acquisition, C.D., M.A. and P.C.; Investigation, N.L.; Methodology, N.L., P.C., L.D., E.B.-M. and M.-C.C.; Project administration, M.A., P.C. and M.-C.C.; Supervision, C.D., R.O., H.C., M.A., P.C. and E.B.-M.; Validation, N.L. and E.B.-M.; Writing—original draft, N.L.; Writing—review and editing, C.D., R.O., H.C., P.C., L.D., E.B.-M. and M.-C.C.

**Funding:** This research was funded by the project PHC-Cèdre (bilateral French–Lebanese funding program) and the Lebanese Agricultural Research Institute (LARI).

**Acknowledgments:** The authors warmly thank and express their gratitude to Jacqueline Saad Harfouche for her help in the survey implementation.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

## References

1. Bou Khouzam, R.; Pohl, P.; Al Ayoubi, B.; Jaber, F.; Lobinski, R. Concentrations of toxic and essential elements in Lebanese bread. *Pure Appl. Chem.* **2012**, *84*, 181–190. [CrossRef]

2. Marín, S.; Pardo, O.; Sánchez, A.; Sanchis, Y.; Vélez, D.; Devesa, V.; Font, G.; Yusà, V. Assessment of metal levels in foodstuffs from the Region of Valencia (Spain). *Toxicol. Rep.* **2018**, *5*, 654–670. [CrossRef] [PubMed]

3. Chaza, C.; Rayane, S.; Sopheak, N.; Moomen, B.; Baghdad, O. Distribution of Organochlorine Pesticides and Heavy Metals in Lebanese Agricultural Soil: Case Study—Plain of Akkar. *Int. J. Environ. Res.* **2018**, *12*, 631–649. [CrossRef]

4. Baize, D.; Bellanger, L.; Tomassone, R. Relationships between concentrations of trace metals in wheat grains and soil. *Agron. Sustain. Dev.* **2009**, *29*, 297–312. [CrossRef]

5. Gupta, N.; Yadav, K.K.; Kumar, V.; Kumar, S.; Chadd, R.P.; Kumar, A. Trace elements in soil-vegetables interface: Translocation, bioaccumulation, toxicity and amelioration—A review. *Sci. Total Environ.* **2019**, *651*, 2927–2942. [CrossRef]

6. Mitra, G. Essential Plant Nutrients and Recent Concepts about their Uptake. In *Essential Plant Nutrients: Uptake, Use Efficiency, and Management*; Naeem, M., Ansari, A.A., Gill, S.S., Eds.; Springer International Publishing: Cham, Switzerland, 2017; pp. 3–36. ISBN 978-3-319-58841-4.

7. Counter, S.A.; Buchanan, L.H. Mercury exposure in children: A review. *Toxicol. Appl. Pharmacol.* **2004**, *198*, 209–230. [CrossRef]
8. (Ine)Van der Fels-Klerx, H.J.; Edwards, S.G.; Kennedy, M.C.; O'Hagan, S.; O'Mahony, C.; Scholz, G.; Steinberg, P.; Chiodini, A. A framework to determine the effectiveness of dietary exposure mitigation to chemical contaminants. *Food Chem. Toxicol.* **2014**, *74*, 360–371. [CrossRef]

9. Arnich, N.; Sirot, V.; Rivière, G.; Jean, J.; Noël, L.; Guérin, T.; Leblanc, J.-C. Dietary exposure to trace elements and health risk assessment in the 2nd French Total Diet Study. *Food Chem. Toxicol.* **2012**, *50*, 2432–2449. [CrossRef]

10. Barone, G.; Dambrosio, A.; Storelli, A.; Busco, A.; Ioanna, F.; Quaglia, N.C.; Giacominelli-Stuffler, R.; Storelli, M.M. Traditional Italian cheeses: Trace element levels and estimation of dietary intake. *J. Food Compos. Anal.* **2018**, *66*, 205–211. [CrossRef]

11. Tang, J.; Huang, Z.; Pan, X.-D. Exposure assessment of heavy metals (Cd, Hg, and Pb) by the intake of local foods from Zhejiang, China. *Environ. Geochem. Health* **2014**, *36*, 765–771. [CrossRef]

12. Ali, M.H.H.; Al-Qahtani, K.M. Assessment of some heavy metals in vegetables, cereals and fruits in Saudi Arabian markets. *Egypt. J. Aquat. Res.* **2012**, *38*, 31–37. [CrossRef]

13. Nasreddine, L.; Nashalian, O.; Naja, F.; Itani, L.; Parent-Massin, D.; Nabhani-Zeidan, M.; Hwalla, N. Dietary exposure to essential and toxic trace elements from a Total diet study in an adult Lebanese urban population. *Food Chem. Toxicol.* **2010**, *48*, 1262–1269. [CrossRef] [PubMed]

14. EFSA. Scientific Opinion on the risks to public health related to the presence of nickel in food and drinking water. *EFSA J. 2015*, *13*, 4002.

15. EFSA. Scientific Opinion on the risk for public health related to the presence of mercury and methylmercury in food. *EFSA J. 2012*, *10*, 2985.

16. EFSA. Scientific Opinion on Arsenic in Food. *EFSA J. 2009*, *7*, 1351. [CrossRef]

17. ANSES. Second French Total Diet Study (TDS 2) Report 1. Inorganic Contaminants, Minerals, Persistent Organic Pollutants, Mycotoxins and Phytoestrogens. ANSES Opinion. Available online: https://www.anses.fr/en/system/files/PASER2006sa0361Ra1EN.pdf (accessed on 26 July 2019).

18. Ahmed, M.T.; Abdel Hadi, E.-S.; El Samahy, S.; Youssof, K. The influence of baking fuel on residues of polycyclic aromatic hydrocarbons and heavy metals in bread. *J. Hazard. Mater.* **2000**, *80*, 1–8. [CrossRef]

19. Ghoreishy, F.; Salehi, M.; Fallahzade, J. Cadmium and lead in rice grains and wheat breads in Isfahan (Iran) and human health risk assessment. *Hum. Ecol. Risk Assess. An Int. J.* **2019**, *25*, 924–934. [CrossRef]

20. Mathebula, M.W.; Mandiwana, K.; Panichev, N. Speciation of chromium in bread and breakfast cereals. *Food Chem.* **2017**, *217*, 655–659. [CrossRef]

21. Oyekunle, J.A.O.; Adekunle, A.S.; Ogunfowokan, A.O.; Olutona, G.O.; Omolere, O.B. Bromate and trace metal levels in bread loaves from outlets within Ile-Ife Metropolis, Southwestern Nigeria. *Toxicol. Rep.* **2014**, *1*, 224–230. [CrossRef]

22. Soares, M.E.; Vieira, E.; De Lourdes Bastos, M. Chromium speciation analysis in bread samples. *J. Agric. Food Chem.* **2010**, *58*, 1366–1370. [CrossRef]

23. Libnor. *NL 240 Arabic Lebanese Bread*; Libnor: Beirut, Lebanon, 2011.

24. Kicińska, A.; Glichowska, P.; Mamak, M. Micro- and macroelement contents in the liver of farm and wild animals and the health risks involved in liver consumption. *Environ. Monit. Assess.* **2019**, *191*, 132. [CrossRef] [PubMed]

25. EFSA. Dietary Reference Values for nutrients Summary report. *EFSA Sup. Publ.* **2017**, *14*, e15121E.

26. World Health Organization. *Trace Elements in Human Nutrition and Health*; WHO: Geneva, Switzerland, 1996.

27. EFSA. Scientific Opinion on the risks to public health related to the presence of chromium in food and drinking water. *EFSA J. 2014*, *12*, 3595.

28. CAS Central Administration of Statistics. Living Conditions Survey. Available online: http://www.cas.gov.lb/index.php/demographic-and-social-en (accessed on 16 July 2019).

29. Bethlehem, J. *Applied Survey Methods—A Statistical Perspective*; Wiley: Hoboken, NJ, USA, 2009.

30. R Core Team. *R: A language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2018.

31. RStudio Team. *RStudio: Integrated Development for R*; RStudio Inc.: Boston, MA, USA, 2016.

32. Sprent, P.; Smeeton, N.C. *Applied Nonparametric Statistical Methods*; Chapman & Hall/CRC: Boca Raton, FL, USA, 1989.

33. Mendiburu, F. *De Agricolae: Statistical Procedures for Agricultural Research*; Scientific Research Publishing: Wuhan, China, 2019.
8. (Ine)Van der Fels-Klerx, H.J.; Edwards, S.G.; Kennedy, M.C.; O'Hagan, S.; O'Mahony, C.; Scholz, G.; Steinberg, P.; Chiodini, A. A framework to determine the effectiveness of dietary exposure mitigation to chemical contaminants. *Food Chem. Toxicol.* **2014**, *74*, 360–371. [CrossRef]

9. Arnich, N.; Sirot, V.; Rivière, G.; Jean, J.; Noël, L.; Guérin, T.; Leblanc, J.-C. Dietary exposure to trace elements and health risk assessment in the 2nd French Total Diet Study. *Food Chem. Toxicol.* **2012**, *50*, 2432–2449. [CrossRef]

10. Barone, G.; Dambrosio, A.; Storelli, A.; Busco, A.; Ioanna, F.; Quaglia, N.C.; Giacominelli-Stuffler, R.; Storelli, M.M. Traditional Italian cheeses: Trace element levels and estimation of dietary intake. *J. Food Compos. Anal.* **2018**, *66*, 205–211. [CrossRef]

11. Tang, J.; Huang, Z.; Pan, X.-D. Exposure assessment of heavy metals (Cd, Hg, and Pb) by the intake of local foods from Zhejiang, China. *Environ. Geochem. Health* **2014**, *36*, 765–771. [CrossRef]

12. Ali, M.H.H.; Al-Qahtani, K.M. Assessment of some heavy metals in vegetables, cereals and fruits in Saudi Arabian markets. *Egypt. J. Aquat. Res.* **2012**, *38*, 31–37. [CrossRef]

13. Nasreddine, L.; Nashalian, O.; Naja, F.; Itani, L.; Parent-Massin, D.; Nabhani-Zeidan, M.; Hwalla, N. Dietary exposure to essential and toxic trace elements from a Total diet study in an adult Lebanese urban population. *Food Chem. Toxicol.* **2010**, *48*, 1262–1269. [CrossRef] [PubMed]

14. EFSA. Scientific Opinion on the risks to public health related to the presence of nickel in food and drinking water. *EFSA J. 2015*, *13*, 4002.

15. EFSA. Scientific Opinion on the risk for public health related to the presence of mercury and methylmercury in food. *EFSA J. 2012*, *10*, 2985.

16. EFSA. Scientific Opinion on Arsenic in Food. *EFSA J. 2009*, *7*, 1351. [CrossRef]

17. ANSES. Second French Total Diet Study (TDS 2) Report 1. Inorganic Contaminants, Minerals, Persistent Organic Pollutants, Mycotoxins and Phytoestrogens. ANSES Opinion. Available online: https://www.anses.fr/en/system/files/PASER2006sa0361Ra1EN.pdf (accessed on 26 July 2019).

18. Ahmed, M.T.; Abdel Hadi, E.-S.; El Samahy, S.; Youssof, K. The influence of baking fuel on residues of polycyclic aromatic hydrocarbons and heavy metals in bread. *J. Hazard. Mater.* **2000**, *80*, 1–8. [CrossRef]

19. Ghoreishy, F.; Salehi, M.; Fallahzade, J. Cadmium and lead in rice grains and wheat breads in Isfahan (Iran) and human health risk assessment. *Hum. Ecol. Risk Assess. An Int. J.* **2019**, *25*, 924–934. [CrossRef]

20. Mathebula, M.W.; Mandiwana, K.; Panichev, N. Speciation of chromium in bread and breakfast cereals. *Food Chem.* **2017**, *217*, 655–659. [CrossRef]

21. Oyekunle, J.A.O.; Adekunle, A.S.; Ogunfowokan, A.O.; Olutona, G.O.; Omolere, O.B. Bromate and trace metal levels in bread loaves from outlets within Ile-Ife Metropolis, Southwestern Nigeria. *Toxicol. Rep.* **2014**, *1*, 224–230. [CrossRef]

22. Soares, M.E.; Vieira, E.; De Lourdes Bastos, M. Chromium speciation analysis in bread samples. *J. Agric. Food Chem.* **2010**, *58*, 1366–1370. [CrossRef]

23. Libnor. *NL 240 Arabic Lebanese Bread*; Libnor: Beirut, Lebanon, 2011.

24. Kicińska, A.; Glichowska, P.; Mamak, M. Micro- and macroelement contents in the liver of farm and wild animals and the health risks involved in liver consumption. *Environ. Monit. Assess.* **2019**, *191*, 132. [CrossRef] [PubMed]

25. EFSA. Dietary Reference Values for nutrients Summary report. *EFSA Sup. Publ.* **2017**, *14*, e15121E.

26. World Health Organization. *Trace Elements in Human Nutrition and Health*; WHO: Geneva, Switzerland, 1996.

27. EFSA. Scientific Opinion on the risks to public health related to the presence of chromium in food and drinking water. *EFSA J. 2014*, *12*, 3595.

28. CAS Central Administration of Statistics. Living Conditions Survey. Available online: http://www.cas.gov.lb/index.php/demographic-and-social-en (accessed on 16 July 2019).

29. Bethlehem, J. *Applied Survey Methods—A Statistical Perspective*; Wiley: Hoboken, NJ, USA, 2009.

30. R Core Team. *R: A language and Environment for Statistical Computing*; R Foundation for Statistical Computing: Vienna, Austria, 2018.

31. RStudio Team. *RStudio: Integrated Development for R*; RStudio Inc.: Boston, MA, USA, 2016.

32. Sprent, P.; Smeeton, N.C. *Applied Nonparametric Statistical Methods*; Chapman & Hall/CRC: Boca Raton, FL, USA, 1989.

33. Mendiburu, F. *De Agricolae: Statistical Procedures for Agricultural Research*; Scientific Research Publishing: Wuhan, China, 2019.
34. Mangiafico, S. *Rcompanion: Functions to Support Extension Education Program Evaluation*; The Comprehensive R Archive Network: New Jersey, NJ, USA, 2019.

35. Ganassali, S. *Enquêtes et Analyse de Données Avec Sphinx*; Pearson Education France: Montreuil, France, 2014; ISBN 978-2-326-00073-5.

36. CAS Central Administration for Statistics. Population of Lebanon. Life expectancy of Lebanon. Density of population in Lebanon. Available online: http://en.worldstat.info/Asia/Lebanon/ (accessed on 16 July 2019).

37. Brereton, R.G. *Applied Chemometrics for Scientists*; John Wiley & Sons: Hoboken, NJ, USA, 2007; ISBN 978-0-470-05777-3.

38. Novotnik, B.; Zuliani, T.; Ščančar, J.; Milačič, R. Chromate in food samples: An artefact of wrongly applied analytical methodology? *J. Anal. At. Spectrom.* **2013**, *28*, 558–566. [CrossRef]

39. ANSES Infant Total Diet Study (iTDS). Available online: https://www.anses.fr/en/system/files/ERCA2010SA0317RaEN.pdf (accessed on 26 July 2019).

40. EFSA. Cadmium in food—Scientific opinion of the Panel on Contaminants in the Food Chain. *EFSA J. 2009*, *7*, 980.

41. EFSA. Scientific Opinion on Lead in Food. *EFSA J. 2010*, *8*, 1570. [CrossRef]

42. Available online: https://www.anses.fr/fr/system/files/PASER2006sa0361.pdf (accessed on 26 July 2019).

43. Nasreddine, L.; Hwalla, N.; Sibai, A.; Hamzé, M.; Parent-Massin, D. Food consumption patterns in an adult urban population in Beirut, Lebanon. *Public Health Nutr.* **2006**, *9*, 194–203. [CrossRef]

44. Mestek, O.; Komínková, J.; Šantrůček, J.; Kačer, P.; Mališová, K.; Koplík, R. Analyses of trace metals, peptide ligands of trace metals and mercury speciation in home prepared bread. *Chem. Speciat. Bioavailab.* **2012**, *24*, 79–88. [CrossRef]

45. European Parliament. Directive 2002/46/EC of the European Parliament and of the Council of 10 June 2002 on the approximation of the laws of the Member States relating to food supplements (Text with EEA relevance). *Off. J. Eur. Communities* **2002**, *183*, 51–57.

46. EFSA. Scientific Opinion on the safety of chromium picolinate as a source of chromium added for nutritional purposes to foodstuff for particular nutritional uses and to foods intended for the general population. *EFSA J. 2010*, *8*, 1883.

47. IARC. IARC Monographs Volume 100C. In *Chromium (VI) Compounds*; International Agency for Research on Cancer: Lyon, France, 2012.

48. Semerjian, L.A. Quality assessment of various bottled waters marketed in Lebanon. *Environ. Monit. Assess.* **2011**, *172*, 275–285. [CrossRef] [PubMed]

49. Ščančar, J.; Milačič, R. A critical overview of Cr speciation analysis based on high performance liquid chromatography and spectrometric techniques. *J. Anal. At. Spectrom.* **2014**, *29*, 427–443. [CrossRef]

50. Kirchmann, H.; Mattsson, L.; Eriksson, J. Trace element concentration in wheat grain: Results from the Swedish long-term soil fertility experiments and national monitoring program. *Environ. Geochem. Health* **2009**, *31*, 561–571. [CrossRef] [PubMed]

51. RASFF. RASFF Portal. Available online: https://webgate.ec.europa.eu/rasff-window/portal/?event=openRecallList&StartRow=44601 (accessed on 26 July 2019).

52. Ergetie, Z.; Berhane, A.; Bayew, T.; Hymete, A. Comparative analysis of Trace metals in bread from five bakeries in Addis Ababa, Ethiopia. *Glob. Adv. Res. J. Pharm. Pharmacol.* **2013**, *2*, 363–367.

53. EWG. EWG's Tap Water Database: Contaminants in Your Water. Available online: https://www.ewg.org/tapwater/system-contaminant.php?pws=MO5010458&contamcode=1081?pws=MO5010458&contamcode=1081 (accessed on 18 July 2019).

54. Available online: https://www.anses.fr/fr/system/files/MCDA2010sa0095.pdf (accessed on 26 July 2019).

55. Available online: https://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2008:353:0001:1355:en:PDF (accessed on 26 July 2019).

56. IARC. IARC Monographs on the Evaluation of Carcinogenic Risks to Humans. In *Cobalt in Hard Metals and Cobalt Sulfate, Gallium Arsenide, Indium Phosphide and Vanadium Pentoxide*; the International Agency for Research on Cancer: Lyon, France, 2006.

57. EFSA. Scientific Opinion on the use of cobalt compounds as additives in animal nutrition. *EFSA J. 2009*, *7*, 1383. [CrossRef]
34. Mangiafico, S. *Rcompanion: Functions to Support Extension Education Program Evaluation*; The Comprehensive R Archive Network: New Jersey, NJ, USA, 2019.

35. Ganassali, S. *Enquêtes et Analyse de Données Avec Sphinx*; Pearson Education France: Montreuil, France, 2014; ISBN 978-2-326-00073-5.

36. CAS Central Administration for Statistics. Population of Lebanon. Life expectancy of Lebanon. Density of population in Lebanon. Available online: http://en.worldstat.info/Asia/Lebanon/ (accessed on 16 July 2019).

37. Brereton, R.G. *Applied Chemometrics for Scientists*; John Wiley & Sons: Hoboken, NJ, USA, 2007; ISBN 978-0-470-05777-3.

38. Novotnik, B.; Zuliani, T.; Ščančar, J.; Milačič, R. Chromate in food samples: An artefact of wrongly applied analytical methodology? *J. Anal. At. Spectrom.* **2013**, *28*, 558–566. [CrossRef]

39. ANSES Infant Total Diet Study (iTDS). Available online: https://www.anses.fr/en/system/files/ERCA2010SA0317RaEN.pdf (accessed on 26 July 2019).

40. EFSA. Cadmium in food—Scientific opinion of the Panel on Contaminants in the Food Chain. *EFSA J. 2009*, *7*, 980.

41. EFSA. Scientific Opinion on Lead in Food. *EFSA J. 2010*, *8*, 1570. [CrossRef]

42. Available online: https://www.anses.fr/fr/system/files/PASER2006sa0361.pdf (accessed on 26 July 2019).

43. Nasreddine, L.; Hwalla, N.; Sibai, A.; Hamzé, M.; Parent-Massin, D. Food consumption patterns in an adult urban population in Beirut, Lebanon. *Public Health Nutr.* **2006**, *9*, 194–203. [CrossRef]

44. Mestek, O.; Komínková, J.; Šantrůček, J.; Kačer, P.; Mališová, K.; Koplík, R. Analyses of trace metals, peptide ligands of trace metals and mercury speciation in home prepared bread. *Chem. Speciat. Bioavailab.* **2012**, *24*, 79–88. [CrossRef]

45. European Parliament. Directive 2002/46/EC of the European Parliament and of the Council of 10 June 2002 on the approximation of the laws of the Member States relating to food supplements (Text with EEA relevance). *Off. J. Eur. Communities 2002*, *183*, 51–57.

46. EFSA. Scientific Opinion on the safety of chromium picolinate as a source of chromium added for nutritional purposes to foodstuff for particular nutritional uses and to foods intended for the general population. *EFSA J. 2010*, *8*, 1883.

47. IARC. IARC Monographs Volume 100C. In *Chromium (VI) Compounds*; International Agency for Research on Cancer: Lyon, France, 2012.

48. Semerjian, L.A. Quality assessment of various bottled waters marketed in Lebanon. *Environ. Monit. Assess.* **2011**, *172*, 275–285. [CrossRef] [PubMed]

49. Ščančar, J.; Milačič, R. A critical overview of Cr speciation analysis based on high performance liquid chromatography and spectrometric techniques. *J. Anal. At. Spectrom.* **2014**, *29*, 427–443. [CrossRef]

50. Kirchmann, H.; Mattsson, L.; Eriksson, J. Trace element concentration in wheat grain: Results from the Swedish long-term soil fertility experiments and national monitoring program. *Environ. Geochem. Health* **2009**, *31*, 561–571. [CrossRef] [PubMed]

51. RASFF. RASFF Portal. Available online: https://webgate.ec.europa.eu/rasff-window/portal/?event=openRecallList&StartRow=44601 (accessed on 26 July 2019).

52. Ergetie, Z.; Berhane, A.; Bayew, T.; Hymete, A. Comparative analysis of Trace metals in bread from five bakeries in Addis Ababa, Ethiopia. *Glob. Adv. Res. J. Pharm. Pharmacol.* **2013**, *2*, 363–367.

53. EWG. EWG's Tap Water Database: Contaminants in Your Water. Available online: https://www.ewg.org/tapwater/system-contaminant.php?pws=MO5010458&contamcode=1081?pws=MO5010458&contamcode=1081 (accessed on 18 July 2019).

54. Available online: https://www.anses.fr/fr/system/files/MCDA2010sa0095.pdf (accessed on 26 July 2019).

55. Available online: https://eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=OJ:L:2008:353:0001:1355:en:PDF (accessed on 26 July 2019).

56. IARC. IARC Monographs on the Evaluation of Carcinogenic Risks to Humans. In *Cobalt in Hard Metals and Cobalt Sulfate, Gallium Arsenide, Indium Phosphide and Vanadium Pentoxide*; the International Agency for Research on Cancer: Lyon, France, 2006.

57. EFSA. Scientific Opinion on the use of cobalt compounds as additives in animal nutrition. *EFSA J. 2009*, *7*, 1383. [CrossRef]
58. EFSA. Scientific Opinion on safety and efficacy of cobalt compounds (E3) as feed additives for all animal species: Cobaltous acetate tetrahydrate, basic cobaltous carbonate monohydrate and cobaltous sulphate heptahydrate, based on a dossier submitted by TREAC EEIG. *EFSA J. 2012*, **10**, 2791.

59. Cuadrado, C.; Kumpulainen, J.; Carbajal, A.; Moreiras, O. Cereals Contribution to the Total Dietary Intake of Heavy Metals in Madrid, Spain. *J. Food Compos. Anal.* **2000**, *13*, 495–503. [CrossRef]

60. Onianwa, P.C.; Lawal, J.A.; Ogunkeye, A.A.; Orejimi, B.M. Cadmium and Nickel Composition of Nigerian Foods. *J. Food Compos. Anal.* **2000**, *13*, 961–969. [CrossRef]

61. Khaniki, G.R.J.; Yunesian, M.; Mahvi, A.H.; Nazmara, S. Trace metal contaminants in Iranian flat breads. *J. Agric. Soc. Sci. (Pak.)* **2005**, *1*, 301–303.

62. IARC Monographs on the Evaluation of Carcinogenic Risks to Humans. *Arsenic, Metals, Fibres and Dusts*; IARC Publication: Lyon, France, 2012.

63. ATSDR. *Toxicological Profile for Nickel*; Agency for Toxic Substances and Disease Registry, US Department of Health and Human Services: Atlanta, GA, USA, 2005.

64. European Commission. Commission Directive 2003/40/EC Commission Directive 2003/40/EC of 16 May 2003 establishing the list, concentration limits and labelling requirements for the constituents of natural mineral waters and the conditions for using ozone-enriched air for the treatment of natural mineral waters and spring waters. *Off. J. Eur. Communities* **2003**, *126*, 34–39.

65. European. Council Directive 98/83/EC Council Directive 98/83/EC of 3 November 1998 on the quality of water intended for human consumption. *Off. J. Eur. Communities* **1998**, *330*, 32–54.

66. Sadee, B.; Foulkes, M.E.; Hill, S.J. Coupled techniques for arsenic speciation in food and drinking water: A review. *J. Anal. At. Spectrom.* **2014**, *30*, 102–118. [CrossRef]

67. IARC. *Overall Evaluations of Carcinogenicity: An Updating of IARC Monographs Volumes 1 to 42*; IARC Publication: Lyon, France, 1987.

68. Vahter, M. Health Effects of Early Life Exposure to Arsenic. *Basic Clin. Pharmacol. Toxicol.* **2008**, *102*, 204–211. [CrossRef] [PubMed]

69. Godfrey, K.M.; Barker, D.J. Fetal nutrition and adult disease. *Am. J. Clin. Nutr.* **2000**, *71*, 1344S–1352S. [CrossRef] [PubMed]

70. Hall, M.; Gamble, M.; Slavkovich, V.; Liu, X.; Levy, D.; Cheng, Z.; van Geen, A.; Yunus, M.; Rahman, M.; Pilsner, J.R.; et al. Determinants of Arsenic Metabolism: Blood Arsenic Metabolites, Plasma Folate, Cobalamin, and Homocysteine Concentrations in Maternal–Newborn Pairs. *Environ. Health Perspect.* **2007**, *115*, 1503–1509. [CrossRef]

71. Concha, G.; Vogler, G.; Lezcano, D.; Nermell, B.; Vahter, M. Exposure to inorganic arsenic metabolites during early human development. *Toxicol. Sci.* **1998**, *44*, 185–190. [CrossRef]

72. Sun, H.-J.; Xiang, P.; Luo, J.; Hong, H.; Lin, H.; Li, H.-B.; Ma, L.Q. Mechanisms of arsenic disruption on gonadal, adrenal and thyroid endocrine systems in humans: A review. *Environ. Int.* **2016**. [CrossRef]

73. Lin, X.; Gu, Y.; Zhou, Q.; Mao, G.; Zou, B.; Zhao, J. Combined toxicity of heavy metal mixtures in liver cells. *J. Appl. Toxicol.* **2016**, *36*, 1163–1172. [CrossRef]

74. Su, H.; Li, Z.; Fiati Kenston, S.S.; Shi, H.; Wang, Y.; Song, X.; Gu, Y.; Barber, T.; Aldinger, J.; Zou, B.; et al. Joint Toxicity of Different Heavy Metal Mixtures after a Short-Term Oral Repeated-Administration in Rats. *Int. J. Environ. Res. Public Health* **2017**. [CrossRef]

75. Kapaj, S.; Peterson, H.; Liber, K.; Bhattacharya, P. Human Health Effects From Chronic Arsenic Poisoning–A Review. *J. Environ. Sci. Health Part A* **2006**, *41*, 2399–2428. [CrossRef]

76. Juhasz, A.L.; Smith, E.; Weber, J.; Rees, M.; Rofe, A.; Kuchel, T.; Sansom, L.; Naidu, R. In Vivo Assessment of Arsenic Bioavailability in Rice and Its Significance for Human Health Risk Assessment. *Environ. Health Perspect.* **2006**, *114*, 1826–1831. [CrossRef]

77. Raber, G.; Stock, N.; Hanel, P.; Murko, M.; Navratilova, J.; Francesconi, K.A. An improved HPLC–ICPMS method for determining inorganic arsenic in food: Application to rice, wheat and tuna fish. *Food Chem.* **2012**, *134*, 524–532. [CrossRef]

78. D'Amato, M.; Aureli, F.; Ciardullo, S.; Raggi, A.; Cubadda, F. Arsenic speciation in wheat and wheat products using ultrasound- and microwave-assisted extraction and anion exchange chromatography-inductively coupled plasma mass spectrometry. *J. Anal. At. Spectrom.* **2011**, *26*, 207–213. [CrossRef]

79. Lebanon Utilities. Lebanon Water Quality Report. Available online: http://www.lebanon-utilities.com/wp-content/uploads/2017/04/Lebanon-Water-Quality-Report-2016.pdf (accessed on 26 July 2019).
58. EFSA. Scientific Opinion on safety and efficacy of cobalt compounds (E3) as feed additives for all animal species: Cobaltous acetate tetrahydrate, basic cobaltous carbonate monohydrate and cobaltous sulphate heptahydrate, based on a dossier submitted by TREAC EEIG. *EFSA J. 2012*, **10**, 2791.

59. Cuadrado, C.; Kumpulainen, J.; Carbajal, A.; Moreiras, O. Cereals Contribution to the Total Dietary Intake of Heavy Metals in Madrid, Spain. *J. Food Compos. Anal.* **2000**, *13*, 495–503. [CrossRef]

60. Onianwa, P.C.; Lawal, J.A.; Ogunkeye, A.A.; Orejimi, B.M. Cadmium and Nickel Composition of Nigerian Foods. *J. Food Compos. Anal.* **2000**, *13*, 961–969. [CrossRef]

61. Khaniki, G.R.J.; Yunesian, M.; Mahvi, A.H.; Nazmara, S. Trace metal contaminants in Iranian flat breads. *J. Agric. Soc. Sci. (Pak.)* **2005**, *1*, 301–303.

62. IARC Monographs on the Evaluation of Carcinogenic Risks to Humans. *Arsenic, Metals, Fibres and Dusts*; IARC Publication: Lyon, France, 2012.

63. ATSDR. *Toxicological Profile for Nickel*; Agency for Toxic Substances and Disease Registry, US Department of Health and Human Services: Atlanta, GA, USA, 2005.

64. European Commission. Commission Directive 2003/40/EC Commission Directive 2003/40/EC of 16 May 2003 establishing the list, concentration limits and labelling requirements for the constituents of natural mineral waters and the conditions for using ozone-enriched air for the treatment of natural mineral waters and spring waters. *Off. J. Eur. Communities* **2003**, *126*, 34–39.

65. European. Council Directive 98/83/EC Council Directive 98/83/EC of 3 November 1998 on the quality of water intended for human consumption. *Off. J. Eur. Communities* **1998**, *330*, 32–54.

66. Sadee, B.; Foulkes, M.E.; Hill, S.J. Coupled techniques for arsenic speciation in food and drinking water: A review. *J. Anal. At. Spectrom.* **2014**, *30*, 102–118. [CrossRef]

67. IARC. *Overall Evaluations of Carcinogenicity: An Updating of IARC Monographs Volumes 1 to 42*; IARC Publication: Lyon, France, 1987.

68. Vahter, M. Health Effects of Early Life Exposure to Arsenic. *Basic Clin. Pharmacol. Toxicol.* **2008**, *102*, 204–211. [CrossRef] [PubMed]

69. Godfrey, K.M.; Barker, D.J. Fetal nutrition and adult disease. *Am. J. Clin. Nutr.* **2000**, *71*, 1344S–1352S. [CrossRef] [PubMed]

70. Hall, M.; Gamble, M.; Slavkovich, V.; Liu, X.; Levy, D.; Cheng, Z.; van Geen, A.; Yunus, M.; Rahman, M.; Pilsner, J.R.; et al. Determinants of Arsenic Metabolism: Blood Arsenic Metabolites, Plasma Folate, Cobalamin, and Homocysteine Concentrations in Maternal–Newborn Pairs. *Environ. Health Perspect.* **2007**, *115*, 1503–1509. [CrossRef]

71. Concha, G.; Vogler, G.; Lezcano, D.; Nermell, B.; Vahter, M. Exposure to inorganic arsenic metabolites during early human development. *Toxicol. Sci.* **1998**, *44*, 185–190. [CrossRef]

72. Sun, H.-J.; Xiang, P.; Luo, J.; Hong, H.; Lin, H.; Li, H.-B.; Ma, L.Q. Mechanisms of arsenic disruption on gonadal, adrenal and thyroid endocrine systems in humans: A review. *Environ. Int.* **2016**. [CrossRef]

73. Lin, X.; Gu, Y.; Zhou, Q.; Mao, G.; Zou, B.; Zhao, J. Combined toxicity of heavy metal mixtures in liver cells. *J. Appl. Toxicol.* **2016**, *36*, 1163–1172. [CrossRef]

74. Su, H.; Li, Z.; Fiati Kenston, S.S.; Shi, H.; Wang, Y.; Song, X.; Gu, Y.; Barber, T.; Aldinger, J.; Zou, B.; et al. Joint Toxicity of Different Heavy Metal Mixtures after a Short-Term Oral Repeated-Administration in Rats. *Int. J. Environ. Res. Public Health* **2017**. [CrossRef]

75. Kapaj, S.; Peterson, H.; Liber, K.; Bhattacharya, P. Human Health Effects From Chronic Arsenic Poisoning–A Review. *J. Environ. Sci. Health Part A* **2006**, *41*, 2399–2428. [CrossRef]

76. Juhasz, A.L.; Smith, E.; Weber, J.; Rees, M.; Rofe, A.; Kuchel, T.; Sansom, L.; Naidu, R. In Vivo Assessment of Arsenic Bioavailability in Rice and Its Significance for Human Health Risk Assessment. *Environ. Health Perspect.* **2006**, *114*, 1826–1831. [CrossRef]

77. Raber, G.; Stock, N.; Hanel, P.; Murko, M.; Navratilova, J.; Francesconi, K.A. An improved HPLC-ICPMS method for determining inorganic arsenic in food: Application to rice, wheat and tuna fish. *Food Chem.* **2012**, *134*, 524–532. [CrossRef]

78. D'Amato, M.; Aureli, F.; Ciardullo, S.; Raggi, A.; Cubadda, F. Arsenic speciation in wheat and wheat products using ultrasound- and microwave-assisted extraction and anion exchange chromatography-inductively coupled plasma mass spectrometry. *J. Anal. At. Spectrom.* **2011**, *26*, 207–213. [CrossRef]

79. Lebanon Utilities. Lebanon Water Quality Report. Available online: http://www.lebanon-utilities.com/wp-content/uploads/2017/04/Lebanon-Water-Quality-Report-2016.pdf (accessed on 26 July 2019).
80. Devesa, V.; Vélez, D.; Montoro, R. Effect of thermal treatments on arsenic species contents in food. *Food Chem. Toxicol.* **2008**, *46*, 1-8. [CrossRef] [PubMed]

81. Tejera, R.L.; Luis, G.; González-Weller, D.; Caballero, J.M.; Gutiérrez, Á.J.; Rubio, C.; Hardisson, A. Metals in wheat flour; comparative study and safety control. *Nutr. Hosp.* **2013**. [CrossRef]

82. ATSDR. *Environmental Medicine (CSEM) Lead Toxicity*; Agency for Toxic Substances and Disease Registry, US Department of Health and Human Services: Atlanta, GA, USA, 2010.

<figure><img src="image_6.png" /></figure>

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).
80. Devesa, V.; Vélez, D.; Montoro, R. Effect of thermal treatments on arsenic species contents in food. *Food Chem. Toxicol.* **2008**, *46*, 1–8. [CrossRef] [PubMed]

81. Tejera, R.L.; Luis, G.; González-Weller, D.; Caballero, J.M.; Gutiérrez, Á.J.; Rubio, C.; Hardisson, A. Metals in wheat flour; comparative study and safety control. *Nutr. Hosp.* **2013**. [CrossRef]

82. ATSDR. *Environmental Medicine (CSEM) Lead Toxicity*; Agency for Toxic Substances and Disease Registry, US Department of Health and Human Services: Atlanta, GA, USA, 2010.

<figure><img src="image_7.png" /></figure>

© 2019 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (http://creativecommons.org/licenses/by/4.0/).