Next Article in Journal
Single-Cell Transcriptomics Reveals Conserved Regulatory Networks in Human and Mouse Interneuron Development
Next Article in Special Issue
Current Levels of Environmental Exposure to Cadmium in Industrialized Countries as a Risk Factor for Kidney Damage in the General Population: A Comprehensive Review of Available Data
Previous Article in Journal
Small-Cell Lung Cancer—An Update on Targeted and Immunotherapies
Previous Article in Special Issue
CoNiCrMo Particles, but Not TiAlV Particles, Activate the NLRP3 Inflammasome in Periprosthetic Cells
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Reference Values on Children’s Hair for 28 Elements (Heavy Metals and Essential Elements) Based on a Pilot Study in a Representative Non-Contaminated Local Area

Institute of Bioengineering, Miguel Hernandez University of Elche, 03202 Elche, Spain
*
Author to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(9), 8127; https://doi.org/10.3390/ijms24098127
Submission received: 13 March 2023 / Revised: 19 April 2023 / Accepted: 25 April 2023 / Published: 1 May 2023
(This article belongs to the Special Issue Heavy Metal Toxicity in Humans 2.0)

Abstract

:
Studies have been published, and laboratories offer services of measuring elements in hair as biomarkers of environmental exposure and/or control of essential elements (trace or macro). These reported values can have only sense if compared with adopted reference values. In this work, we propose provisional reference values based on a pilot child population. The concentrations of 28 elements were measured in children’s hair samples. An observational, descriptive, cross-sectional study was conducted in a typical child population in the Mediterranean region void of excessive pollution problems to analyze 419 hair samples of children aged 3–12 years. Children were selected by a simple random method from eight primary education schools in different municipal districts, which included urban, rural and industrial areas. Samples of around 100 mg were washed and acid digested by an optimized procedure. All measures were performed using ICP-MS with Sc, Y and Re as internal standards. The statistical analysis was performed by two approaches: (a) considering all the data and (b) without outliers (second-order atypical data) to compare them with other published studies. The distribution curves in all the elements studied were asymmetric and did not fit the theoretical normality distributions. Therefore, the analysis based on percentiles was more appropriate. In most elements, only slight differences were observed with sex or age, which did not justify proposing separate reference ranges. From the results of this study, provisional reference values are proposed following two criteria: (a) simple application of the table of percentiles built by removing outlier values and (b) values after a detailed analysis case-by-case, considering other data as the distribution profile and other published data of each element. Although the pilot sample was from a limited area, it was carefully selected to be representative of a general non-contaminated population. With this limitation, the proposed reference values might be useful for researchers and physicians until a wider geographical study is available for a large number of elements.

1. Introduction

Harm to health by pollution caused by elements of concern, such as Pb, Hg, Cd, Tl and As, among others, is well-known, and the toxicity of such elements has been extensively studied in recent decades. Interest has been shown in assessing exposure to these elements and their effects on children and their development. The route for exposure might be from ingestion in food [1], drinks [2] or contaminated particles in the air [3] and others. The nomenclature of “trace elements” is frequently used for elements presented in biological samples at low concentrations, either for toxic or biologically essential elements. Essential minerals may be elements encountered in small amounts (Fe, Mn, Cu, I, Co, Zn, F and Se) and so-called macro-minerals, which are required in large amounts (Ca, P, Mg, Na, K and Cl) [4].
Some trace elements are considered probably or possibly essential as some beneficial effects have been suggested (i.e., Si and B), and some elements are considered essential by some food agencies but not others, such as Cr. The concentrations of essential minerals in the organism must be balanced, and they may harm health if lacking or found in excess. Some examples are night blindness related to lack of Zn in infant diets, especially in developing countries [5]; lack of Fe, also known as iron deficiency anemia; and problems related with hepatic damage caused by an excess of Cu [6]. They are all frequent problems but rarely harm the health severely. Given this delicate balance, it is important to monitor exposure levels in specific populations and vulnerable age groups, such as children.
Therefore, biomonitoring and assessing incidence on health and on infants’ development is a world concern and a priority for toxic and essential elements. Hair samples have been used as both a non-invasive method and a biomarker of exposure to analyze these elements. Routine methods require between 100 and 150 mg of hair to undergo chemical methods [7]. These samples are ideal for mega-analyses with healthy populations to conduct population studies of interest for public health issues.
Special efforts have been made to evaluate and biomonitor mercury in infants, and relevant studies have been conducted with Pb and Cd. There is particular concern for areas whose gastronomy is rich in sea products since bioaccumulation capacity may play a role. By focusing on Hg, the series of studies performed with children who live on the Seychelles and Faroe Islands is likely the most widely studied background available [8,9,10,11,12,13,14,15,16,17,18]. Some more limited studies exist, which have been conducted in other geographical areas, such as the Mediterranean, where a great deal of fish and seafood is eaten [19,20].
Efforts have been made for larger cohort studies in relationship to neurodevelopment and Hg exposure in Central Europe and Mediterranean countries [19,20]. Other studies have been performed in Japan with a general population, which essentially focused on Hg and other elements, such as Pb or Cd [21,22,23,24] or specifically on a child population [22,23,25,26]. A biomonitoring study in adults in a mining area in Colombia was reported in relationship with renal impairment [27]. Studies focused on children have been reported in Malaysia [28], China [29,30], Brazil [31,32], Poland [33], Czech [34] and Slovenia [35] as well as related to Au mining in Perú [36,37]. Studies relating Hg exposure with autism and other neurological problems are a focus of current attention [38,39,40].
For individual acute-phase studies, known tests that use blood and urine samples can be appropriate. The use of other biological sample, such as autopsy tissues, have been considered, especially for biomonitoring chronic exposure. However, an evaluation in a case study in a population potentially exposed to metals near a hazardous waste incinerator concluded that using autopsy tissues did not show advantages in comparison to other appropriate samples for biomonitoring, such as blood [41].
Yet when studying the problem in large populations, it is necessary to rely on more efficient user-friendly non-invasive biomarkers as far as sampling is concerned. Hair covers all the needed requirements: (1) It is a tissue of an accumulation type that grows 1 cm every month on average. It enables studies to evaluate chronic exposure and allows a temporal relationship to be established in accordance with the hair segment under study. (2) Obtaining hair samples does not involve aggressive invasive methods but merely requires simple techniques with minimum personnel training. (3) No special means or specific care are required to handle, conserve or transport samples. (4) It is a low-cost practice.
Nonetheless, the main problem of human biomonitoring in hair is that reference values are lacking for most elements in the general population, which is stressed even more for children.
Relevant studies of biomonitoring in hair have been performed with Hg and some with other elements of concern, such as Pb and Cd. However, very few studies have been conducted in multi-element biomonitoring in children’s hair [42,43]. Studies with a limited number of elements have been reported, mainly monitoring elements, such as Pb, Cd, MN, Cr, Fe, Ni, Zn, Al and As in addition to Hg [29,44,45].
In some cases, studies are in relationship to smoking, dietary or industrial activity in the area; gender influence; or military attacks [46,47,48,49,50]. More recently, attention has been paid to assess heavy metals in relation to autism [48,49] and other developmental impairment in children and adolescents [51,52,53] as well as the relationship with secondhand cigarette smoke exposure [54], but most studies deal with only a few elements, mainly Hg, Cd, Al, Co, Cu, Ni and Cr. One study in Ethiopia evaluated up to 41 elements in the hair of children. Although this was limited to 81 subjects [55], it is a useful recent study.
For this work, a pilot child population in the Mediterranean was selected. The study area was the municipality of Elche, in the province of Alicante (east Spain), which encompasses urban, rural, farming and industrial areas. The industrial areas neither are heavy industrial areas, nor are they heavily polluted, and this province has no agricultural overproduction areas. Therefore, the population of this area can be considered to be representative of a typical general population of the Mediterranean region without special heavy contamination, and, in this way, we were able to apply a viable logistic to obtain a sufficient statistical number of subjects and an appropriate distribution inside the population.
An observational, descriptive, cross-sectional study was conducted. The total population, aged under 14 years, was 34,109 individuals of both sexes, according to the census (Table 1), from which a representative sample was obtained of 419 children in the range of 3–12 years old. The concentrations of 28 mineral elements were determined in hair samples using ICP-MS technology. The distribution curves for each element were described, and the descriptive statistics and sex and age distributions were analyzed. From the results of this study, provisional reference values are proposed, which may be used in clinical diagnoses and epidemiological studies by public health professionals using hair element concentrations as biomarkers of environmental exposure.

2. Results

2.1. General Statistics and Distribution

Hair samples from 419 children were selected following the statistic and inclusion/exclusion criteria described in Material and Methods for the design of the study. Samples were washed, dried, digested and analyzed using ICP-MS for measuring the concentrations of 28 elements based on calibration curves with data from certified standards (Table S1).
The general statistics of the obtained data are presented in Table 2. The measuring unit is micrograms of element per gram of dry hair samples.
All the analyzed elements presented asymmetric distribution curves (ϒ1 > 0.5) that did not fit the theoretical normality distributions (K-S < 0.05), which makes the mean and standard deviation (SD) weak if used as central tendency and dispersion parameters. All the distribution curves were leptokurtic (ϒ2 > 0.5), and none were mesokurtic. Zn and Al were the elements whose kurtosis coefficients came closest to zero: ϒ2 = 2, for both. By way of example, the heavy metals Hg, Pb, Cd and As are shown in Figure 1. As examples, the most symmetric curves found (Zn, Fe and Se) are presented in Figure 2.
The proportion of the SD values vs. the X ¯ A values of the elements was high, with %SD that exceeded 50%. The variability was very high. The elements with the lowest %SD were Zn and Fe (Zn: 9.48 ± 7.34 µg/g; and Fe: 483.4 ± 291.8 µg/g), with 60% and 77%, respectively. In certain cases (i.e., As) the SD tripled the mean value. This broadness in the ranges of concentrations occurred with most elements, except for U, Be, Ag and Tl, whose scarcity in the samples did not enable suitable observations to be made. Thus, it was considered more suitable to base the data analysis on calculating the percentiles. Such percentile-based analyses considering all data are shown in Table 3 and might be considered for proposing reference values based on the 5–95th percentiles.
However, due to the high asymmetry of the distribution and the existence of very extreme atypical values, approaches considering eliminating such values were considered as follows.

2.2. Analysis by Eliminating Extreme Values

The quartiles technique was used to determine criteria to select the limits of the extreme values (see the criteria in the Section 4). Both degrees of atypicality were calculated by excess: first and second degrees (Table S2). As the calculated representative population sample had to contain 380 individuals, to maintain N above this figure, we selected the second-degree value (excess) as the limit of atypicality. All the values over this limit were considered atypical.
Table S3 shows the general statistics after eliminating the atypical second-degree (excess) values and the values below the LOD. The distribution curves of essential elements, such as Zn and Se, became more symmetric with asymmetry indices close to zero (Zn: ϒ1 = 0.4; and Se: ϒ1 = −0.3. When considering the example of Fe, despite being asymmetric (ϒ1 = 0.9), the contrast was evident when comparing the A panels with the B panels in Figure 2. The asymmetric distribution remained in the heavy metals, even after eliminating the atypical values (ϒ1 > 0.5).
This makes the mean and the SD weak as central tendency and dispersion measures. Therefore, we calculated a table of the values of the percentiles after eliminating the atypical values and those below the LOD (Table S4).

2.3. Crossing of Age and Sex Variables

The distribution of the elements in the study population did not meet the normality assumption. Therefore, for the crossing of variables, non-parametric tests were applied: Kruskal–Wallis for age and the Mann–Whitney U test for sex.

2.3.1. Age and Study Elements

The ages of the children in years were grouped into three mutually exclusive 3-year class intervals: 3–6 years (n = 156); 7–9 years (n = 162); and 10–12 years (n = 101). Differences were found among these groups in the concentrations of the following elements for age: Be, Mg, Al, K, Ca, Mn, Co, Cu, Zn, As, Se, Sr, Ba, Bi, Tl, Ag, Au and U (p < 0.01) (Table 4).
The age variable was processed as ordinal categorical and numerical. As a numerical variable, the degree of the correlation found with the study elements was measured using Spearman’s Rho coefficient, and correlation coefficients were found within a range from −0.266 (As was the most negative) to 0.311 (Sr was the most positive).
In all the elements for which a correlation with age was found, with both a falling and rising tendency, the correlation coefficients were <0.5, that is, the correlations were weak, scarce or null, which means that, although they were observed, they were not strong enough to explain the upward/downward tendency of an element according to age. Moreover, this tendency was influenced by other variables.

2.3.2. Sex and Study Elements

Differences in the concentrations were observed as higher in boys for B, Na, K, Zn and As but were higher in girls for Mg, Ca, Co, Ni, Cu, Sr and Ba (p < 0.01) with no differences for Al, Cr, Mn, Fe, Se, Mo, Cd, Hg, Pb and Bi (Table 5).
From the results of this study, we considered that the differences by age and sex were not sufficient to justify separate analysis for proposing reference values inside the child population.

3. Discussion

The main purpose of this paper is to propose reference values for concentrations of elements in the hair of children. We based this work on the statistical data obtained from the hair samples of a representative population in the Mediterranean region that were analyzed to obtain concentrations of heavy metals and essential elements by ICP-MS technology.
The distribution curves per element were described, and the descriptive statistics and the sex and age distributions were analyzed to propose the references values, which are presented in this discussion. Attempts to establish reference values of a wide number of multiple elements in hair have been published regarding university students in Poland [56] and schoolboys in Rome [40], and an interesting review was published [57].

3.1. Statistical Considerations

In many of the published studies, especially on Hg, authors have eliminated extreme values in order to reduce the asymmetry of the distribution as much as possible, with the goal of better applying the parametric statistical calculations.
We considered it appropriate to include all data obtained from the sample of the population and, therefore, to use non-parametric statistics, which is recommended for the analysis of data not adjusted to a theoretical normal distribution (see the general statistics in Table 2 and percentiles in Table 3). Nevertheless, to compare our data with those published by other authors, in this work, we also show descriptive statistical analyses after eliminating extreme (excess) and below the LOD (zero) values (Table S3) and the corresponding table of percentiles (Table S4).
Even doing this last, the asymmetry is maintained, so non-parametric statistics continue to be the technique of choice, and the mean and standard deviation lose strength. Therefore, we proposed reference values based on the percentiles.

3.2. Sex and Age Variables

While it is true that differences were found in the concentrations of several elements in relation to the age and sex variables (Table 4 and Table 5) and that these differences fit with those found by other authors in related works [42,43] we consider that they are not sufficient to justify proposing tables of references values of for each age and sex group.
The margins of reference that we propose are in a sufficiently wide range to include the entire study population in the same table. For practical purposes, we propose a single table of references, without different margins in relation to age and sex.

3.3. Criteria to Propose Reference Values for Elements in Hair of Children: Criteria

The main objective of this work was to propose a provisional tentative table of reference values for the concentrations of heavy metals and essential elements using hair as the matrix analysis of the biomarkers of exposure. Selecting these reference margins was based on statistical analysis without considering clinical criteria. To this end, a representative population sample was selected: school children. We studied two approaches to obtain reference values:
(1)
First proposal: Based on the statistics of the data presented herein, mainly using 5–95th percentiles after eliminating outliers (Table S4) or in some cases from the percentile table considering the total population of study (Table 3). The distribution curves of all the studied elements were completely asymmetric, and none fitted theoretical normality. Hence, the means and standard deviations had a low statistical value, and the percentiles were considered more suitable. Although the 5–95% range of percentiles might be the recommendable criteria, professionals can use other ranges that they consider to be more suitable in each case.
(2)
Second proposal: Setting values for the specific analysis of each element by considering, additionally to the percentile table, other data as follows: (i) the characteristics of the statistical distribution shown with and without eliminating extreme outliers and (ii) comparing our results in published works with a similar population in a case-by-case approach. The intention here was for them to be of general potential use. A descriptive analysis of the criteria applied for each element is shown in Supplementary Table S5. The main bibliographic references used to apply such criteria are shown inside that table.
The proposed values are summarized in Table 6 for both criteria.

4. Materials and Methods

4.1. Design of Study

4.1.1. Type of Study, Place and Time

An observational, descriptive, cross-sectional study was conducted in the municipality of Elche, located in the province of Alicante, east Spain. Selection of study subjects: a simple random method; time: between September and December 2006.

4.1.2. The Total Study Population and Population Sample

The population data, obtained from the Spanish National Institute of Statistics (INE), were tabulated by age groups in mutually exclusive 4-year intervals. The total population, aged under 14 years, constituted the total study population in the municipality of Elche.
There were 34,109 individuals of both sexes, according to the census updated in September 2006, with this being the field sampling starting time. The data are shown in Table 1. The population sample was selected from elementary school children under the age of 12 years old. Conformity of Ethical Committee of University Miguel Hernandez (reference: 2017.167.E.OEP, Ref.: IB.EVG.02.17).
Sample calculation: The sample was calculated using the statistical formula recommended for finite populations: n = N*Z2α*p*q/d2*(N − 1) + Z2α*p*, where N = 34,109 was the total study population of children in the area. A confidence level of 95% (which corresponds with a preset value of Zα = 1.96) was considered. The prevalence of the problem was unknown at the time; therefore, the ratio value was also unknown. Hence, we used p = 0.5 (50%) to maximize the sample size; q = 1–9 = 0.5, and the standard precision or acceptable limit error of the sample was 5% (d = 0.05). On this basis, a value of n = 380 individual was estimated. Increasing the sample size by 5% is usually recommended for dealing with any unforeseen absence of responses, records, processing errors, samples that are not useful, etc. In our study, the final increase in samples for contingencies was 10%.
The guest sample was made of 1200 individuals preselected to be invited to participate. From them, the collected sample were 422 individuals attempting to be representative of the different subareas and age–sex distribution. Finally, we had 419 useful samples producing data.
Inclusion criteria: children aged 3–12 years; both sexes; residents in the municipality of Elche; any ethnic group; not undergoing any hair treatments (dyes, treatments to avoid hair falling out, revitalizing products or similar) in the 3 months prior to sample collection; with sufficient hair length to obtain a useful sample; and their parents had signed the informed consent.
Sampling method: probabilistic, randomly simple. Schools were selected within each school district. The cooperation of the local Education Authorities was fundamental in this phase of the work. The project was presented to the Students–Parents Association, school directors and teachers. The distribution of the data collection instrument (survey and informed consent) for the guest sample was performed randomly. The selection of subjects was conducted on site, randomly from among those who attended on the day of sample collection and from among those whose parents had provided authorized written permission to participate in the study.
Sample collection: 419 useful samples were collected in eight schools, one school per district. The weight and height were measured, and the age and sex were recorded. A hair sample was obtained by cutting a small section of hair from the occipital area, close to the scalp. Each sample was kept in a self-closing plastic bag and was coded.

4.1.3. Study Variables

Sex: male and female; Age groups: 3–6; 7–9; and 10–12 years; Elements (µg/g in dry hair): Be, B, Na, Mg, Al, K, Ca, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, As, Se, Sr, Mo, Cd, Ba, Hg, Tl, Pb, Bi, Au, Ag and U.

4.2. Analytical Procedures

4.2.1. Hair Sample Processing

Hair samples (0.1–0.3 g each) were immersed in 200 mL Triton X 100 at 2% for 1 min in an ultrasound bath. They were rinsed three times with ultrapure water, each time for 1 min. Hair was dried for 24 h at 60·°C. Hair samples were cut with surgical steel scissors into ≤1 mm pieces to facilitate digestion. Scissors were previously treated with 0.6% nitric acid for 10 min. Approximately 0.1 g of sample was placed into pre-weighted Falcon tubes, and the exact weights were noted. Then, they were centrifuged at 5000 rpm for 5 min for moving hair pieces to the bottom of the tube. Next, 1 mL of HNO3 65% Suprapur® was added, and the samples were left in an acid digestion for 48 h.
After digestion, they were transferred to clean tubes, and 9 mL of ultrapure water was added to adjust to 10 mL to give a nitric acid concentration of up to 6.5%. Finally, 5 mL of each sample was placed in previously prepared tubes with 0.5 mL of the internal standard mix to obtain a final volume of 5.5 mL in each tube. The remaining 5 mL was stored as backup.
Several washing procedures were previously checked, and a decision was made to apply an intensive (ultrasound and detergent) but short (1 min) treatment to eliminate superficial materials and to avoid removing the very soluble ions, such as Na and Mg, from inside the hair as these are not tightly fixed to the hair tissue.

4.2.2. ICP MS Measures

The equipment used was an Inductively Coupled Plasma Mass Spectrometer (ICP-MS), Agilent® Model 7500A (Agilent Technologies Spain S.L., Madrid, Spain). The sampling pump was programmed to 1 mL/min, 30 s for washing tubing, 30 s for stabilizing detection and then 3 s by measuring per mass in duplicate (total time of around 2 min per sample). Approximately 160 samples can be measured in 8 h. Therefore, the 419 hair samples and calibration standard solutions were measured in three batches. For each batch, the procedure was repeated: first, the ICP-MS and TUNE standards were measured, followed by hair samples and blank controls. A standard was introduced every 20 samples.

4.2.3. ICP-MS Certified Standards and Calibration

The certified standard solution used for calibrating the equipment and for the semiquantitative analysis was Agilent Technologies Tuning solution 10 µg/L Li, Y, Ce, Tl, Co (Agilent Technologies Spain S.L., Madrid, Spain)). Sc (m/z 45), Y (m/z 89) and Re (m/z 185) ICP standards 1000 mg/L, CertiPUR® (Merck, Darmstadt, Germany)), were used as internal standards (ISTD). For this purpose, a mix solution was prepared that contained the three elements at (100 µg/L in 5% HNO3). A 0.5 mL volume of the mix was added to 5 mL of each sample or to a calibration standard.
Calibration standard for the full quantitative data analysis. Twelve standard solutions were prepared by serial dilution (ratio 1:2) from CertiPUR® Merck ICP multi-element standard solution X (23 elements in nitric acid) and the mercury ICP/MS standard by dilution with nitric acid 65% Suprapur® (Merck). They were used for the calibration curves of each element under study Table S1. For Cu and Zn, three other standards at higher concentrations were prepared using zinc sulfate 7-hydrate PANREAC® and copper sulfate pentahydrate PANREAC® (Panreac Química S.L.U., Barcelona, Spain) because the concentrations of these elements in many hair samples exceeded the ICP multi-element standard solution concentration. To 5 mL of each calibration standard solution, a 0.5 mL volume of the mix of the internal standards was added formed by a solution containing 100 µg/L of Sc, Y and Re (CertiPUR®, Merck) in 5% HNO3).
Calibration curves. Calibration curves were obtained by the mathematical function equation that best conformed to the calibration curve points within the range of concentrations within which all the elements were found in the hair samples.
Laboratory blank control and standard control. Three blank controls were made by each set of analyzed samples and were subjected to the same process as the hair samples. For checking the equipment variability over time, the Tuning solution and a calibration solution were measured as a control every 20 samples.
Validation Reference Materials. Validation of the methodology was performed testing with several reference materials: IAEA-086 (un-spiked), obtained from International Atomic Energy, Analytical Quality Control Service, Vienna, Austria; NIES No. 13 obtained from National Institute for Environmental Studies as Certified Reference Material No. 13 “Human Hair”; and ERM-DB001 obtained from European Reference material, Join Research Center, Institute for Reference Materials and Measurements).
Data analysis. The counts obtained for elements Be, B, Na, Mg, Al, K, Ca, V, Cr, Mn, Fe, Co, Ni, Cu, Zn, As, Se, Sr, Mo, Cd, Ba, Hg, Tl, Pb and Bi were analyzed by the full quantitative method included in the equipment’s software. For this purpose, the signals, corrected by the internal standard, were compared with the calibration curves deduced in the standard solutions per element. However, for Ag, Au and U, the semiquantitative analysis option of the equipment’s software was used based on the calibrations performed with the aforementioned tune solutions. The data obtained with the ICP-MS computer software were copied into an an electronic sheet to transform the units mathematically from µg/L into µg/g of dry hair according to the weight of each sample. Data were copied to the SSPS software (IBM-SPSS Statistics 18-22) to be statistically analyzed.

4.3. Statistical Approaches

Asymmetry indices and normality test. The Fisher asymmetry coefficient (ϒ1) and kurtosis coefficient (ϒ2) for each element were calculated to observe the behavior of the distribution curves. The K-S normality test was applied to check if the data distribution in the study population met the normality assumption.
Estimation of central tendency and dispersion. Several central tendency and dispersion parameters were calculated as follows: arithmetic mean ( X ¯ A) with its respective 95% confidence intervals; standard deviation (SD); median (Me); range of concentrations (minimum–maximum); and percentile values 5%, 10%, 15%, 25%, 50%, 65%, 75%, 85%, 90% and 95%.
To calculate the geometric mean ( X ¯ G) and to compare our data with those published by other authors, the values that equaled zero (below the limit of detection, LOD) were removed and excluded from the calculations.
Data tabulation: Data were grouped into tables per categories, calculated using the formula of Sturges: class interval = range/K, where the K value = 1 + 3.322 (Log 10N), where K is the number of categories and N is the total study population. Histograms were constructed according to these tables.
Outlier values: The elements distribution in the study population did not meet the normality assumption (see Section 2.1). For this reason, the standard deviation lost strength as a statistical parameter of dispersion in our data. Therefore, the quartiles technique was used to identify the upper boundary outliers, both first-degree and second-degree. For that, the value of the interquartile range (IQR) was calculated by the difference between the third quartile (Q3) and the first one (Q1): IQR = Q3 − Q1. The limits of atypical values, were obtained as follows:
First-degree limit of atypical value:
Excess: >Q3 + (IQR × 1.5).
Second-degree limit of atypical value:
Excess: >Q3 + (IQR × 3).

5. Conclusions

Several studies have been reported that have used measures of biomarkers in hair, especially with Hg and as well as studies with heavy metals, such as Pb and Cd. Very few studies have been published with multi-elemental values in children’s hair. Special consideration is needed for Hg. The Mediterranean population is a high fish and shellfish consumer, which explains the higher level compared with other populations. Therefore, the value based on the 95th percentile for Hg is likely only representative for similar high-fish-consuming populations.
For all elements, the proposed value suggested with the second criteria takes this factor into consideration for suggesting a value for a more general potential application. The data and arguments supporting this are discussed in Table S5, case by case, for each element. The proposed reference values for the concentrations of elements in the hair of children must be cautiously considered as provisional until obtaining values obtained with more universal geographical distribution studies and multi-center-based data. Nevertheless, the area of study and the sample were carefully selected to be representative of the general child population, including individuals in in different rural, urban and industrial areas and with, in each subarea (school), a careful distribution of age, sex and other random characteristics, including children of different social, economic, cultural and ethnic conditions.
Therefore, we consider that the sample should be very similar to a general child population of the whole country and likely also to other similar Mediterranean countries. Our approach allowed us to apply viable logistics to obtain an appropriate sample. Moreover, the values proposed with the second criteria consider other data that allow us to consider it reasonable to be applied provisionally for other wider populations.
In any case, the reference proposed in Table 6 for concentrations of elements in hair of children should be useful as a reference for biomarkers of exposure by environmental health professionals and epidemiologists at least for preliminary references for comparison of data of individuals or groups.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24098127/s1. References [58,59,60,61,62,63,64,65,66,67,68,69,70,71,72,73,74,75,76,77,78,79,80,81,82] are cited in the Supplementary Materials

Author Contributions

Conceptualization, E.V. and R.R.; methodology, R.R. and C.A.; software, R.R. and C.E.; validation, J.E. and M.A.S.; formal analysis, investigation, resources, data curation and writing—original draft preparation, E.V. and R.R. writing—review and editing, M.A.S. and C.E.; visualization, supervision, project administration and funding acquisition, E.V. All authors have read and agreed to the published version of the manuscript.

Funding

This work was partly supported by a grant from the Regional Government (Generalitat Valenciana GV 2008-FP7. Cap4. Sección 10 Centro gestor 61 DG Salud Pública. línea: T58940), by institutional annual budgets of the University Miguel Hernández of Elche and partially support by “Asociación Biotox”, Alicante Spain.

Institutional Review Board Statement

The procedure were approved by the accredited Ethical Committee of the Miguel Hernandez University, under the registration code 2017.167.E.OEP, and approval code IB.EVG.02.17, to the projects: “Desarrollo y aplicación de biomarcadores en menores, mediante métodos no cruentos, para establecer la exposición ambiental a metales pesados para evaluar los riegos en el desarrollo y salud infantile”.

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study. WE had their parents’ approval, and the nice collaborations of the Students–Parents Association, school directors and teachers.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Otim, E.O.; Kasozi, K.I.; Zirintunda, G.; Matama, K. Investigating vegetable contamination in indigent communities by heavy metals: A case of food safety in Bushenyi, Uganda. Environ. Bioavailab. 2021, 33, 292–300. [Google Scholar] [CrossRef]
  2. Otim, E.O.; Chen, R.; Otim, O. Applying multivariate analysis to characterize waragi spirits from Acoli, Uganda, by their metal contents. Heliyon 2019, 5, e01417. [Google Scholar] [CrossRef] [PubMed]
  3. Kumari, S.; Jain, M.K.; Elumalai, S.P. Assessment of Pollution and Health Risks of Heavy Metals in Particulate Matter and Road Dust Along the Road Network of Dhanbad, India. J. Health Pollut. 2021, 11, 210305. [Google Scholar] [CrossRef] [PubMed]
  4. Zoroddu, M.A.; Aashet, J.; Crisponi, G.; Medici, S.; Peana, M.; Nurchi, V.M. The essential metals for humans: A brief overview. J. Inorg. Biochem. 2019, 195, 120–129. [Google Scholar] [CrossRef]
  5. Kazi, T.G.; Afridi, H.; Kazi, G.H.; Jamali, M.K.; Arain, M.B.; Jalbani, N. Evaluation of essential and toxic metals by ultrasound-assisted acid leaching from scalp hair samples of children with macular degeneration patients. Clin. Chim. Acta 2006, 369, 52–60. [Google Scholar] [CrossRef]
  6. Scheinberg, I.H.; Sternlieb, I. Wilson disease and idiopathic copper toxicosis. Am. J. Clin. Nutr. 1996, 63, 842S–845S. [Google Scholar] [CrossRef]
  7. Legrand, M.; Sousa Passos, C.J.; Mergler, D.; Chan, H.M. Biomonitoring of Mercury Expositure with Single Hair. Environ. Sci. Technol. 2005, 39, 4594–4598. [Google Scholar] [CrossRef]
  8. Cernichiari, E.; Myers, G.J.; Ballatori, N.; Zareba, G.; Vyas, J.; Clarkson, T. The biological monitoring of prenatal exposure to methylmercury. Neurotoxicology 2007, 28, 1015–1022. [Google Scholar] [CrossRef]
  9. Cernichiari, E.; Toribara, T.Y.; Liang, L.; Marsh, D.O.; Berlin, M.W.; Myers, G.J.; Cox, C.; Shamlaye, C.F.; Choisy, O.; Davidson, P. The biological monitoring of mercury in the Seychelles study. Neurotoxicology 1995, 16, 613–628. [Google Scholar]
  10. Davidson, P.W.; Palumbo, D.; Myers, G.J.; Cox, C.; Shamlaye, C.F.; Sloane-Reeves, J.; Cernichiari, E.; Wilding, G.E.; Clarkson, T. Neurovelopmental Outcomes of Seychellois Children from the Pilot Cohort al 108 Month Following Prenatal Exposure to Methylmercury from a Maternal Fish Diet. Environ. Res. Sect. A 2000, 84, 1–11. [Google Scholar] [CrossRef]
  11. Debes, F.; Bultz-Jorgensen, E.; Weihe, P.; White, F.R.; Grandjean, P. Impact of prenatal methylmercury esposure on neuribehavioral funtion at age 14 year. Neurotoxicol. Teratol. 2006, 28, 363–375. [Google Scholar] [CrossRef] [PubMed]
  12. Debes, F.; Weihe, P.; Grandjean, P. Cognitive deficits at age 22 years associated with prenatal exposure to methylmercury. Cortex 2016, 74, 358–369. [Google Scholar] [CrossRef] [PubMed]
  13. Grandjean, P.; Weihe, P.; White, R.; Bebes, F.; Araki, S.; Yokoyama, K.; Murata, K.; Sorensen, N.; Dahl, R.; Jorgensen, P. Cognitive déficit in 7 year-old children with prenatal exposure to methylmercury. Neurotoxicol. Teratol. 1997, 19, 417–428. [Google Scholar] [CrossRef] [PubMed]
  14. Myers, G.J.; Davidson, P.W.; Palumbo, D.; Shamlaye, C.; Cox, C.; Cernichiari, E.; Clarkson, T.W. Secondary Analisis from the Seycheles Child Developmente Study: The Chile Behavior Ckecklist. Environ. Res. 2000, 84, 12–19. [Google Scholar] [CrossRef]
  15. Myers, G.J.; Marsh, D.O.; Davidson, P.W.; Cox, C.; Shamlaye, C.F.; Tanner, M.; Choi, A.; Cernichiari, E.; Choisy, O.; Clarkson, T.W. Main neurodevelopmental study of Seychellois children following in utero exposure to methylmercury from a maternal fish diet: Outcome at six months. Neurotoxicology 1995, 16, 653–664. [Google Scholar]
  16. Myers, G.J.; Thurston, S.W.; Pearson, A.T.; Davidson, P.W.; Cox, C.; Shamlaye, C.F.; Cernichiari, E.; Clarkson, T.W. Postnatal exposure to methyl mercury from fish consumption: A review and new data from the Seychelles Child Development Study. Neurotoxicology 2009, 30, 338–349. [Google Scholar] [CrossRef]
  17. Weihe, P.; Grandjean, P.; Debes, F.; White, R. Health implications for Faroe Islander of heavy metal and PCBs from pilot whales. Sci. Total Environ. 1996, 186, 141–148. [Google Scholar] [CrossRef]
  18. Weihe, P.; Grandjean, P.; Jørgensen, P.J. Application of hair-mercury analysis to determine the impact of a seafood advisory. Environ. Res. 2005, 97, 201–208. [Google Scholar] [CrossRef]
  19. Nadal, M.; Bocio, A.; Schuhmacher, M.; Domingo, J.L. Monitoring metals in the population living in the vicinity of a hazardous waste incinerator: Levels in hair of school children. Biol. Trace Elem. Res. 2005, 104, 203–213. [Google Scholar] [CrossRef]
  20. Torrente, M.; Colomina, M.-T.; Domingo, J. Metal Concentrations in Hair and Cognitive Assessment in an Adolescent Population. Biol. Trace Elem. Res. 2005, 104, 215. [Google Scholar] [CrossRef]
  21. Schwedler, G.; Seiwert, M.; Fiddicke, U.; Ißleb, S.; Hölzer, J.; Nendza, J.; Wilhelm, M.; Wittsiepe, J.; Koch, H.M.; Schindler, B.K.; et al. Human biomonitoring pilot study DEMOCOPHES in Germany: Contribution to a harmonized European approach. Int. J. Hyg. Environ. Health 2017, 220, 686–696. [Google Scholar] [CrossRef] [PubMed]
  22. Barbone, F.; Rosolen, V.; Mariuz, M.; Parpinel, M.; Casetta, A.; Sammartano, F.; Ronfani, L.; Vecchi Brumatti, L.; Bin, M.; Castriotta, L.; et al. Prenatal mercury exposure and child neurodevelopment outcomes at 18 months: Results from the Mediterranean PHIME cohort. Int. J. Hyg. Environ. Health 2018. [Google Scholar] [CrossRef] [PubMed]
  23. Matsumoto, H.; Koya, G.; Takeuchi, T. Fetal Minamata disease: A neuropathological study of two cases of intrauterine intoxication by a methyl mercury compound. J. Neuropathol. Exp. Neurol. 1965, 24, 563–574. [Google Scholar] [CrossRef] [PubMed]
  24. Ohba, T.; Kurokawa, N.; Nakai, K.; Shimada, M.; Suzuki, K.; Sugawara, N.; Kameo, S.; Satoh, C.; Satoh, H. Permanent waving does not change mercury concentration in the proximal segment of hair close to scalp. Tohoku J. Exp. Med. 2008, 214, 69–78. [Google Scholar] [CrossRef]
  25. Kusanagi, E.; Takamura, H.; Chen, S.J.; Adachi, M.; Hoshi, N. Children’s Hair Mercury Concentrations and Seafood Consumption in Five Regions of Japan. Arch. Environ. Contam. Toxicol. 2018, 74, 259–272. [Google Scholar] [CrossRef]
  26. Tatsuta, N.; Nakai, K.; Iwai-Shimada, M.; Suzuki, T.; Satoh, H.; Murata, K. Total mercury levels in hair of children aged 7years before and after the Great East Japan Earthquake. Sci. Total Environ. 2017, 596, 207–211. [Google Scholar] [CrossRef]
  27. Rodriguez-Villamizar, L.A.; Medina, O.M.; Flórez-Vargas, O.; Vilanova, E.; Idrovo, A.J.; Araque-Rodriguez, S.A.; Henao, J.A.; Sánchez-Rodríguez, L.H. Chemical Element Mixtures and Kidney Function in Mining and Non-Mining Settings in Northern Colombia. Int. J. Environ. Res. Public Health 2023, 20, 2321. [Google Scholar] [CrossRef]
  28. Abdul Samad, N.I.; Md Isa, Z.; Hod, R. Mercury Hair Concentration among Primary School Children in Malaysia. Children 2017, 4, 109. [Google Scholar] [CrossRef]
  29. Li, Z.; Su, H.; Wang, L.; Hu, D.; Zhang, L.; Fang, J.; Jin, M.; Fiati Kenston, S.S.; Song, X.; Shi, H.; et al. Epidemiological Study on Metal Pollution of Ningbo in China. Int. J. Environ. Res. Public Health 2018, 15, 424. [Google Scholar] [CrossRef]
  30. Jia, Q.; Zhu, X.; Hao, Y.; Yang, Z.; Wang, Q.; Fu, H.; Yu, H. Mercury in soil, vegetable and human hair in a typical mining area in China: Implication for human exposure. J. Environ. Sci. 2018, 68, 73–82. [Google Scholar] [CrossRef]
  31. Dos Santos Freitas, J.; da Costa Brito Lacerda, E.M.; da Silva Martins, I.C.V.; Rodrigues, D., Jr.; Bonci, D.M.O.; Cortes, M.I.T.; Corvelo, T.C.O.; Ventura, D.F.; de Lima Silveira, L.C.; da Conceição Nascimento Pinheiro, M.; et al. Cross-sectional study to assess the association of color vision with mercury hair concentration in children from Brazilian Amazonian riverine communities. Neurotoxicology 2018, 65, 60–67. [Google Scholar] [CrossRef] [PubMed]
  32. Barbosa, A.C.; Silva, S.R.; Dórea, J.G. Concentration of mercury in hair of indigenous mothers and infants from the Amazon basin. Arch. Environ. Contamitation Toxicol. 1998, 34, 100–105. [Google Scholar] [CrossRef] [PubMed]
  33. Marcinek-Jacel, M.; Albińska, J.; Pawlaczyk, A.; Szynkowska, M.I. The impact of demographic factors, behaviors and environmental exposure to mercury content in the hair of the population living in the region of Lodz (central Poland). Environ. Toxicol. Pharm. 2017, 55, 196–201. [Google Scholar] [CrossRef]
  34. Forysová, K.; Pinkr-Grafnetterová, A.; Malý, M.; Krsková, A.; Mráz, J.; Kašparová, L.; Čejchanová, M.; Sochorová, L.; Rödlová, S.; Černá, M. Urinary Cadmium and Cotinine Levels and Hair Mercury Levels in Czech Children and Their Mothers Within the Framework of the COPHES/DEMOCOPHES Projects. Arch. Environ. Contam. Toxicol. 2017, 73, 421–430. [Google Scholar] [CrossRef]
  35. Kobal, A.B.; Snoj Tratnik, J.; Mazej, D.; Fajon, V.; Gibičar, D.; Miklavčič, A.; Kocman, D.; Kotnik, J.; Sešek Briški, A.; Osredkar, J.; et al. Exposure to mercury in susceptible population groups living in the former mercury mining town of Idrija, Slovenia. Environ. Res. 2017, 152, 434–445. [Google Scholar] [CrossRef] [PubMed]
  36. Barenys, M.; Boix, N.; Farran-Codina, A.; Palma-Linares, I.; Montserrat, R.; Curto, A.; Gomez-Catalan, J.; Ortiz, P.; Deza, N.; Llobet, J.M. Heavy metal and metalloids intake risk assessment in the diet of a rural population living near a gold mine in the Peruvian Andes (Cajamarca). Food Chem. Toxicol. 2014, 71, 254–263. [Google Scholar] [CrossRef] [PubMed]
  37. Weinhouse, C.; Ortiz, E.J.; Berky, A.J.; Bullins, P.; Hare-Grogg, J.; Rogers, L.; Morales, A.M.; Hsu-Kim, H.; Pan, W.K. Hair Mercury Level is Associated with Anemia and Micronutrient Status in Children Living Near Artisanal and Small-Scale Gold Mining in the Peruvian Amazon. Am. J. Trop. Med. Hyg. 2017, 97, 1886–1897. [Google Scholar] [CrossRef]
  38. Skalny, A.V.; Simashkova, N.V.; Skalnaya, A.A.; Klyushnik, T.P.; Bjørklund, G.; Skalnaya, M.G.; Tinkov, A.A. Assessment of gender and age effects on serum and hair trace element levels in children with autism spectrum disorder. Metab. Brain Dis. 2017, 32, 1675–1684. [Google Scholar] [CrossRef]
  39. Gustin, K.; Tofail, F.; Mehrin, F.; Levi, M.; Vahter, M.; Kippler, M. Methylmercury exposure and cognitive abilities and behavior at 10years of age. Environ. Int. 2017, 102, 97–105. [Google Scholar] [CrossRef]
  40. Wahlberg, K.; Love, T.M.; Pineda, D.; Engström, K.; Watson, G.E.; Thurston, S.W.; Yeates, A.J.; Mulhern, M.S.; McSorley, E.M.; Strain, J.J.; et al. Maternal polymorphisms in glutathione-related genes are associated with maternal mercury concentrations and early child neurodevelopment in a population with a fish-rich diet. Environ. Int. 2018, 115, 142–149. [Google Scholar] [CrossRef]
  41. Domingo, J.L.; García, F.; Nadal, M.; Schuhmacher, M. Autopsy tissues as biological monitors of human exposure to environmental pollutants. A case study: Concentrations of metals and PCDD/Fs in subjects living near a hazardous waste incinerator. Environ. Res. 2017, 154, 269–274. [Google Scholar] [CrossRef] [PubMed]
  42. Park, H.-S.; Shin, K.-O.; Kim, J.-S. Assessment of Reference Values for Hair Minerals of Korean Preschool Children. Biol. Trace Elem. Res. 2007, 116, 119–130. [Google Scholar] [CrossRef]
  43. Senafonte, O.; Violante, N. Caroli S Assessement of refereence values of elements in human hair of urban schoolboys. J. Trace Elem. Med. Biol. 2000, 14, 6–13. [Google Scholar] [CrossRef] [PubMed]
  44. Rodrigues, J.L.G.; Araújo, C.F.S.; Dos Santos, N.R.; Bandeira, M.J.; Anjos, A.L.S.; Carvalho, C.F.; Lima, C.S.; Abreu, J.N.S.; Mergler, D.; Menezes-Filho, J.A. Airborne manganese exposure and neurobehavior in school-aged children living near a ferro-manganese alloy plant. Environ. Res. 2018, 167, 66–77. [Google Scholar] [CrossRef]
  45. Jursa, T.; Stein, C.R.; Smith, D.R. Determinants of Hair Manganese, Lead, Cadmium and Arsenic Levels in Environmentally Exposed Children. Toxics 2017, 22, 19. [Google Scholar] [CrossRef]
  46. Zhu, Y.; Li, Z.; Pang, Y.; Huo, W.; Li, N.; Li, Z.; Zhang, J.; Ye, R.; Wang, B. Association Between Chronic Exposure to Tobacco Smoke and Accumulation of Toxic Metals in Hair Among Pregnant Women. Biol. Trace Elem. Res. 2018, 185, 302–310. [Google Scholar] [CrossRef] [PubMed]
  47. Zhou, J.; Liang, J.; Hu, Y.; Zhang, W.; Liu, H.; You, L.; Zhang, W.; Gao, M.; Zhou, J. Exposure risk of local residents to copper near the largest flash copper smelter in China. Sci. Total Environ. 2018, 630, 453–461. [Google Scholar] [CrossRef] [PubMed]
  48. Xie, W.; Peng, C.; Wang, H.; Chen, W. Health Risk Assessment of Trace Metals in Various Environmental Media, Crops and Human Hair from a Mining Affected Area. Int. J. Environ. Res. Public Health 2017, 14, 1595. [Google Scholar] [CrossRef]
  49. Vigeh, M.; Yokoyama, K.; Matsukawa, T.; Shinohara, A.; Shariat, M.; Ohtani, K. Effects of Hair Metals on Body Weight in Iranian Children Aged 20 to 36 Months. Iran. J. Public Health 2017, 46, 1018–1027. [Google Scholar]
  50. Manduca, P.; Diab, S.Y.; Qouta, S.R.; Albarqouni, N.M.; Punamaki, R.L. A cross sectional study of the relationship between the exposure of pregnant women to military attacks in 2014 in Gaza and the load of heavy metal contaminants in the hair of mothers and newborns. BMJ Open 2017, 7, e014035. [Google Scholar] [CrossRef]
  51. Eqani, S.A.M.A.S.; Khuram, F.; Alamdar, A.; Tahir, A.; Shah, S.T.A.; Nasir, A.; Javed, S.; Bibi, N.; Hussain, A.; Rasheed, H.; et al. Environmental exposure pathway analysis of trace elements and autism risk in Pakistani children population. Sci. Total Environ. 2020, 712, 136471. [Google Scholar] [CrossRef]
  52. Amadi, C.N.; Orish, C.N.; Frazzoli, C.; Orisakwe, O.E. Association of autism with toxic metals: A systematic review of case-control studies. Pharm. Biochem. Behav. 2022, 212, 173313. [Google Scholar] [CrossRef] [PubMed]
  53. Yasuda, H.; Tsutsui, T. Metallomics analysis for early assessment and individualized intervention of neurodevelopmental disorders. Metallomics 2022, 14, mfac067. [Google Scholar] [CrossRef] [PubMed]
  54. Ballbè, M.; Fu, M.; Masana, G.; Pérez-Ortuño, R.; Gual, A.; Gil, F.; Olmedo, P.; García-Algar, Ó.; Pascual, J.A.; Fernández, E. Passive exposure to electronic cigarette aerosol in pregnancy: A case study of a family. Environ. Res. 2023, 216 Pt 1, 114490. [Google Scholar] [CrossRef]
  55. Astolfi, M.L.; Pietris, G.; Mazzei, C.; Marconi, E.; Canepari, S. Element Levels and Predictors of Exposure in the Hair of Ethiopian Children. Int. J. Environ. Res. Public Health. 2020, 17, 8652. [Google Scholar] [CrossRef]
  56. Chojnacka, K.; Zielińska, A.; Gorecka, H.; Dobrzański, Z.; Gorecki, H. Reference values for hair minerals of Polish students. Environ. Toxicol. Pharmacol. 2010, 29, 314–319. [Google Scholar] [CrossRef]
  57. Mikulewicz, M.; Chojnacka, K.; Gedrange, T.; Górecki, H. Reference values of elements in human hair: A systematic review. Environ. Toxicol. Pharmacol. 2013, 36, 1077–1086. [Google Scholar] [CrossRef]
  58. Fido, A.; Al-Saad, S. Toxic trace elements in the hair of children with autism. Autism 2005, 9, 290–298. [Google Scholar] [CrossRef]
  59. Blaurock-Busch, E.; Amin, O.R.; Rabah, T. Heavy metals and trace elements in hair and urine of a sample of arab children with autistic spectrum disorder. Maedica 2011, 6, 247–257. [Google Scholar]
  60. Elenge, M.M.; Aubry, J.C.; Jacob, L.; De Brouwer, C. Heavy metal in hair samples of 109 non-industrial (miners) population in Katanga. Sante 2011, 21, 41–46. [Google Scholar] [CrossRef]
  61. Ferré-Huguet, N.; Nadal, M.; Schuhmacher, M.; Domingo, J.L. Monitoring Metals in Blood and Hair of the Population Living Near a Hazardous Waste Incinerator: Temporal Trend. Biol. Trace Elem. Res. 2009, 128, 191–199. [Google Scholar] [CrossRef] [PubMed]
  62. Gil, F.; Hernández, A.F.; Márquez, C.; Femia, P.; Olmedo, P.; López-Guarnido, O.; Pla, A. Biomonitorization of cadmium, chromium, manganese, nickel and lead in whole blood, urine, axillary hair and saliva in an occupationally exposed population. Sci. Total Environ. 2011, 409, 1172–1180. [Google Scholar] [CrossRef] [PubMed]
  63. Carneiro, M.F.; Moresco, M.B.; Chagas, G.R.; de Oliveira Souza, V.C.; Rhoden, C.R.; Barbosa, F., Jr. Assessment of Trace Elements in Scalp Hair of a Young Urban Population in Brazil. Biol. Trace Elem. Res. 2011, 143, 815–824. [Google Scholar] [CrossRef] [PubMed]
  64. Blaurock-Busch, E.; Amin, O.R.; Dessoki, H.H.; Rabah, T. Toxic Metals and Essential Elements in Hair and Severity of Symptoms among Children with Autism. Maedica 2012, 7, 38–48. [Google Scholar] [PubMed]
  65. Priya, L.; Geetha, A. Level of Trace Elements (Copper, Zinc, Magnesium and Selenium) and Toxic Elements (Lead and Mercury) in the Hair and Nail of Children with Autism. Biol. Trace Elem. Res. 2011, 142, 148–158. [Google Scholar] [CrossRef] [PubMed]
  66. Doğan-Sağlamtimur, N.; Kumbur, H. Metals (Hg, Pb, Cu, and Zn) Bioaccumulation in Sediment, Fish, and Human Scalp Hair: A Case Study from the City of Mersin Along the Southern Coast of Turkey. Biol. Trace Elem. Res. 2010, 136, 55–70. [Google Scholar] [CrossRef]
  67. Garí, M.; Grimalt, J.O.; Torrent, M.; Sunyer, J. Influence of socio-demographic and diet determinants on the levels of mercury in preschool children from a Mediterranean island. Environ. Pollut. 2013, 182, 291–298. [Google Scholar] [CrossRef]
  68. Díez, S.; Delgado, S.; Aguilera, I.; Astray, J.; Pérez-Gómez, B.; Torrent, M.; Sunyer, J.; Bayona, J. Prenatal and Early Childhood Exposure to Mercury and Methylmercury in Spain, a High-Fish-Consumer Country. Arch. Environ. Contam. Toxicol. 2009, 56, 615–622. [Google Scholar] [CrossRef]
  69. Freire, C.; Ramos, R.; Lopez-Espinosa, M.J.; Díez, S.; Vioque, J.; Ballester, F.; Fernández, M.F. Hair mercury levels, fish consumption, and cognitive development in preschool children from Granada, Spain. Environ. Res. 2010, 110, 96–104. [Google Scholar] [CrossRef]
  70. Batista, J.; Schuhmacher, M.; Domingo, J.L.; Corbella, J. Mercury in hair for a child population from Tarragona Province, Spain. Sci. Total Environ. 1996, 193, 143–148. [Google Scholar] [CrossRef]
  71. Croes, K.; De Coster, S.; De Galan, S.; Morrens, B.; Loots, I.; Van de Mieroop, E.; Nelen, V.; Isabelle Sioen, I.; Bruckers, L.; Nawrot, T.; et al. Health effects in the Flemish population in relation to low levels ofmercury exposure: From organ to transcriptome level. Int. J. Hyg. Environ. Health 2014, 217, 239–247. [Google Scholar] [CrossRef] [PubMed]
  72. Pesch, A.; Wilhelm, M.; Rostek, U.; Schmitz, N.; Weishoff-Houben, M.; Ranft, U.; Idel, H. Mercury concentrations in urine, scalp hair, and saliva in children from Germany. J. Expo. Anal. Environ. Epidemiol. 2002, 12, 252–258. [Google Scholar] [CrossRef]
  73. McDowell, M.; Dillon, C.; Osterloh, J.; Bolger, P.M.; Pellizzari, E.; Fernando, R.; Montes de Oca, R.; Schober, S.; Sinks, T.; Jones, R.; et al. Hair Mercury Levels in U.S. Children and Women of Childbearing Age: Reference Range Data from NHANES 1999–2000. Environ. Health Perspect. 2004, 112, 1165–1171. [Google Scholar] [CrossRef] [PubMed]
  74. Tian, W.; Egeland, G.; Sobol, I.; Chan, H. Mercury hair concentrations and dietary exposure among Inuit preschool children in Nunavut, Canada. Environ. Int. 2011, 37, 42–48. [Google Scholar] [CrossRef]
  75. Murata, K.; Budtz-Jørgensen, E.; Grandjean, P. Benchmark Dose Calculations for Methylmercury- Associated Delays on Evoked Potential Latencies in Two Cohorts of Children. Risk Anal. 2002, 22, 465–474. [Google Scholar] [CrossRef] [PubMed]
  76. Cordier, S.; Garel, M.; Mandereau, L.; Morcel, H.; Doineau, P.; Gosme-Seguret, S.; Josse, D.; White, R.; and Amiel-Tison, C. Neurodevelopmental Investigations among Methylmercury-Exposed Children in French Guiana. Environ. Res. Sect. A 2002, 89, 1–11. [Google Scholar] [CrossRef]
  77. Drasch, G.; Böse-O’Reilly, S.; Beinhoff, C.; Roider, G.; Maydl, S. The Mt. Diwata study on the Philippines 1999 - assessing mercury intoxication of the population by small scale gold mining. Sci. Total Environ. 2001, 267, 151–168. [Google Scholar] [CrossRef]
  78. 78. Agah, H.; Leermakers, M.; Gao, Y.; Fatemi, S.M.R.; Katal, M.M.; Baeyens, W.; Elskens, M. Mercury accumulation in fish species from the Persian Gulf and in human hair from fishermen. Environ. Monit. Assess. 2010, 169, 203–216. [Google Scholar] [CrossRef]
  79. Barbosa, A.C.; Jardim, W.; Dórea, J.G.; Fosberg, B.; Souza, J. Hair mercury speciation as a function of gender, age, and body mass index in inhabitants of the Negro River basin, Amazon, Brazil. Arch. Environ. Contam. Toxicol. 2001, 40, 439–444. [Google Scholar]
  80. (*) UNE (2019) Global Mercury Assessment Technical Background Report 2018. Published 2019. UN Environment. Programme (Geneva, Switzerland). 2023. Available online: https://www.unep.org/resources/publication/global-mercury-assessment-technical-background-report (accessed on 28 April 2023).
  81. Sanna, E.; De Micco, A.; Vallascas, E. Evaluation of Association between Biomarkers of Lead Exposure in Sardinian Children (Italy). Biol. Trace Elem. Res. 2011, 143, 1383–1392. [Google Scholar] [CrossRef]
  82. Seifert, B.; Becker, K.; Helm, D.; Krause, C.; Schulz, C.; Seiwert, M. The German Environmental Survey 1990/1992 (GerES II): Reference concentrations of selected environmental pollutants in blood, urine, hair, house dust, drinking water and indoor air. J. Expo. Anal. Environ. Epidemiol. 2000, 10, 552–565. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Examples of the distribution curves of several toxic elements. The geometric mean ( X ¯ G), median and arithmetic mean are indicated by long dashed, dotted and continuous lines. B Panels show distribution splitting in more range categories than A panels, and do not consider extreme values.
Figure 1. Examples of the distribution curves of several toxic elements. The geometric mean ( X ¯ G), median and arithmetic mean are indicated by long dashed, dotted and continuous lines. B Panels show distribution splitting in more range categories than A panels, and do not consider extreme values.
Ijms 24 08127 g001
Figure 2. Example of distribution curves of several oligoelements. The geometric mean ( X ¯ G), median and arithmetic mean are indicated by long dashed, dotted and continuous lines. B Panels show distribution splitting in more range categories than A Panels and do not consider extreme values.
Figure 2. Example of distribution curves of several oligoelements. The geometric mean ( X ¯ G), median and arithmetic mean are indicated by long dashed, dotted and continuous lines. B Panels show distribution splitting in more range categories than A Panels and do not consider extreme values.
Ijms 24 08127 g002
Table 1. Child population according to sex and age (years old).
Table 1. Child population according to sex and age (years old).
0–405–0910–14Total
Male61755463584817,486
Female58595247551716,623
34,109
Source: Spanish National Institute of Statistics (INE). The population is indicated in 5-year groups in the 2006 review of the municipal census of Elche.
Table 2. General statistics (n = 419). Values are expressed as µg/g of hair.
Table 2. General statistics (n = 419). Values are expressed as µg/g of hair.
Elem.m/zLOD Min–Max X ¯ A ± SDCI 95% Meϒ1ϒ2K-S (p)
Be90.00030.000–0.0130.001 ± 0.0020.001–0.0010.000418<0.01
B110.080.00–15.802.51 ± 2.032.31–2.702.3328<0.01
Na232.80.0–1077142.7 ± 148.7128–15797.728<0.01
Mg240.130.70–355.228.7 ± 29.825.9–31.618.7434<0.01
Al270.390.00–43.397.68 ± 7.806.93–8.435.9812<0.01
K390.560.00–315.134.5 ± 40.930.6–38.52139<0.01
Ca430.80.000–4366406 ± 446363–449278.9425<0.01
V510.00040.000–6.0540.13 ± 0.3460.097–0.1630.07714219<0.01
Cr530.0050.000–17.560.215 ± 0.9780.121–0.3090.13515248<0.01
Mn550.0040.000–1.1390.159 ± 0.1450.15–0.170.11928<0.01
Fe561.220.000–69.1149.48 ± 7.348.7–10.28.78314<0.01
Co590.0010.000–0.7540.016 ± 0.0540.011–0.0220.00610118<0.01
Ni600.0050.000–3.40.278 ± 0.3450.245–0.3110.216531<0.01
Cu630.050.700–863.6105.7 ± 126.594–11867.239<0.01
Zn660.1310.400–1882483.4 ± 291.8455–511485.512<0.01
As750.0030.000–2.1420.034 ± 0.1180.022–0.0450.0215254<0.01
Se820.050.020–19.510.55 ± 1.160.43–0.660.4813189<0.01
Sr880.020.000–50.1984.68 ± 4.734.23–5.143.14320<0.01
Mo950.0010.000–1.1740.051 ± 0.090.043–0.0600.03910111<0.01
Cd1070.0020.000–1.0640.045 ± 0.0750.038–0.0520.025890<0.01
Ba1110.0020.000–3.2570.355 ± 0.3280.32–0.390.289317<0.01
Hg1370.0030.014–48.1223.58 ± 4.773.13–4.042.09427<0.01
Tl2050.00010.000–0.0060.0003 ± 0.0010.0002–0.00030.000425<0.01
Pb2020.030.000–55.3132.19 ± 41.81–2.581.17887<0.01
Bi2050.00010.000–1.9940.032 ± 0.1520.017–0.0460.0059101<0.01
Ag2080.0030.000–0.0380.001 ± 0.0030.001–0.0010876<0.01
Au2090.00020.000–1.270.029 ± 0.1180.018–0.0400758<0.01
U238---------------------
LOD: limit of detection; X ¯ A: arithmetic mean; X ¯ G: geometric mean; SD: standard deviations; Me: median. ϒ1: asymmetry; ϒ2: kurtosis; and KS: Kolmogorov–Smirnov test.
Table 3. Percentiles in the entire study population (n = 419). Values as µg/g of hair.
Table 3. Percentiles in the entire study population (n = 419). Values as µg/g of hair.
5101525506575859095
Be0.0000.0000.0000.0000.0000.0000.0010.0020.0020.004
B0.000.080.191.172.332.863.594.314.896.11
Na6.012.122.344.297.7146.3189.3249.1334.7439.6
Mg2.54.25.98.618.731.239.956.367.881.9
Al0.000.000.220.95.988.8111.9414.518.8324.78
K0.21.83.67.32133.245.36682.3114.4
Ca23.842.371.3127.4278.9422567723.1863.61133.8
V0.0070.0140.020.0320.0770.1180.1510.1960.2370.314
Cr0.000.000.000.0240.1350.1770.1990.2420.2790.357
Mn0.0060.0170.0320.0670.1190.1630.2110.2920.3450.45
Fe0.421.081.975.48.7810.8512.7215.116.921.27
Co0.000.000.000.000.0060.0120.0160.0230.0280.041
Ni0.0000.0150.0230.0620.2160.2860.3650.4680.5590.714
Cu3.45.47.616.267.298.8142.2201.7266371.1
Zn39.577.5131.9285.9485.5568.5634.6759.9851.1978.3
As0.0020.0040.0070.0110.020.0260.0330.0430.0520.068
Se0.040.080.140.360.480.530.560.610.650.79
Sr0.30.50.691.243.145.3179.1310.9913.37
Mo0.0020.0050.0070.020.0390.0520.0630.0770.090.111
Cd0.0020.0040.0070.0130.0250.0380.0490.0780.10.146
Ba0.0210.0360.0560.120.2890.4070.4990.6450.7130.901
Hg0.080.190.320.742.093.414.676.768.5810.98
Tl0.0000.0000.0000.0000.0000.0000.0000.0010.0010.001
Pb0.080.150.210.501.171.922.543.604.337.39
Bi0.0000.0010.0010.0020.0050.0080.0120.020.0380.088
Ag0.0000.0000.0000.0000.0000.0000.0010.0010.0020.003
Au0.0000.0000.0000.0000.0000.0000.0000.0060.0660.171
U0.0000.0000.0000.0000.0000.0000.0000.0000.0000.000
Table 4. Age group, element Kruskal–Wallis test and bivariant correlation (Spearman’s Rho).
Table 4. Age group, element Kruskal–Wallis test and bivariant correlation (Spearman’s Rho).
Rank Average (Median/Arithmetic Mean in µg/g)Bivariant Correlation
Elem.Age 3–6 (n=156)Age 7–9 (n=162)Age 10–12 (n=101)p Value *Rho Spearmanp *
Be194 (0.00/0.00)213 (0.00/0.001)231 (0.00/0.001)0.0140.1230.012
B197 (2.27/2.32)223 (2.42/2.75)210 (2.25/2.42)0.1500.0620.207
Na209 (92/148)204 (95/136)221 (102/146)0.5490.0320.519
Mg178 (14/21)213 (21/28)255 (32/41)0.0000.2840.000
Al234 (8/10)189 (5/6)207 (5/7)0.004−0.1430.003
K231 (29/45)199 (20/30)195 (19/26)0.020−0.1660.001
Ca183 (227/291)213 (305/430)247 (409/546)0.0000.2470.000
V217 (0.08/0.11)200 (0.07/0.16)215 (0.09/0.11)0.391−0.0090.849
Cr211 (0.15/0.24)205 (0.13/0.24)216 (0.14/0.14)0.774−0.0160.737
Mn245 (0.15/0.21)183 (0.11/0.13)199 (0.11/0.14)0.000−0.2090.000
Fe226 (9.7/10)195 (8.2/9)210 (8.2/10)0.083−0.0890.069
Co215 (0.008/0.019)193 (0.005/0.016)231 (0.008/0.014)0.0360.0180.720
Ni196 (0.19/0.25)210 (0.22/0.30)232 (0.25/0.28)0.0630.1090.026
Cu180 (52/74)210 (65/114)256 (126/142)0.0000.2390.000
Zn168 (389/391)227 (529/524)248 (567/561)0.0000.2750.000
As246 (0.024/0.035)200 (0.019/0.042)171 (0.015/0.018)0.000−0.2660.000
Se184 (0.45/0.52)221 (0.49/0.60)233 (0.51/0.50)0.0020.1730.000
Sr175 (2.32/3.41)214 (3.82/4.72)257 (5.44/6.59)0.0000.3110.000
Mo225 (0.04/0.06)204 (0.04/0.05)196 (0.04/0.04)0.128−0.0860.079
Cd218 (0.03/0.04)203 (0.02/0.04)210 (0.03/0.06)0.542−0.0480.330
Ba190 (0.24/0.31)209 (0.29/035)243 (0.36/0.44)0.0030.1970.000
Hg212 (2.10/3.70)199 (1.79/2.98)225 (2.46/4.36)0.2190.0340.489
Pb223 (1.36/2.54)194 (1.08/1.56)215 (1.21/2.67)0.101−0.0610.215
Bi244 (0.007/0.037)193 (0.004/0.035)185 (0.003/0.018)0.000−0.2450.000
Tl203 (0.0002/0.000)218 (0.0003/0.000)207 (0.0002/0.000)0.300.0240.629
Ag225 (0.001/0.000)202 (0.001/0.000)200 (0.001/0.000)0.04−0.1350.006
Au201 (0.0160.000)227 (0.00330.000)196 (0.00420.000)0.010.0830.090
U210 (0.000/0.000)210 (0.000/0.000)209 (0.000/0.000)0.73−0.0380.442
* p value, level of significance < 0.05.
Table 5. Elements according to sex: Mann–Whitney U test. Rank average (median/arithmetic mean in µg/g).
Table 5. Elements according to sex: Mann–Whitney U test. Rank average (median/arithmetic mean in µg/g).
Elem.Male (n = 202)Female (n = 217)p *
Be202 (0.001/0.000)214 (0.001/0.000)0.440
B230 (2.5/2.8)191 (2.0/2.3)0.001
Na230 (110/165)191 (82/122)0.001
Mg186 (16/24)232 (26/33)0.000
Al199 (5.5/7.1)221 (6.4/8.2)0.062
K237 (28/41)185 (17/29)0.000
Ca186 (243/326)232 (345/481)0.000
V181 (0.05/0.09)237 (0.11/0.17)0.000
Cr212 (0.14/0.29)208 (0.13/0.14)0.727
Mn218 (0.13/0.17)203 (0.12/0.15)0.204
Fe214 (8.8/9.7)206 (8.8/9.9)0.534
Co185 (0.004/0.014)234 (0.010/0.018)0.000
Ni196 (0.18/0.24)223 (0.24/0.31)0.026
Cu192 (57/87)227 (89/123)0.003
Zn227 (528/527)194 (462/443)0.005
As240 (0.024/0.030)182 (0.015/0.037)0.000
Se213 (0.47/0.48)207 (0.48/0.60)0.595
Sr184 (2.4/3.8)234 (4.4/5.5)0.000
Mo199 (0.04/0.05)220 (0.04/0.06)0.068
Cd209 (0.025/0.046)211 (0.026/0.044)0.823
Ba194 (0.24/0.32)225 (0.33/0.39)0.007
Hg215 (2.1/3.8)205 (2.1/3.4)0.392
Pb210 (1.15/2.42)210 (1.17/1.98)0.959
Bi213 (0.005/0.024)207 (0.005/0.038)0.609
Tl220 (0.000/0.000)201 (0.000/0.000)0.021
Ag193 (0.001/0.000)226 (0.001/0.000)0.000
Au211 (0.031/2.094)209 (0.027/0.000)0.857
U209 (0.000/0.000)211 (0.000/0.000)0.172
* p value, level of significance < 0.05.
Table 6. Proposal of reference values. The proposed values are indicated as a limit of maximum value or a range (minimum–maximum) for the essential elements. Values are expressed as µg element/g hair, equivalent to mg/kg as indicated by some authors. Data in the first column are values based on the simple criteria of 5–95th percentiles, eliminating the outliers’ values. The proposed values in the right column (in bold) are modulated considering, case by case, the distribution in our statistics and the other published data as discussed in Supplementary Table S5.
Table 6. Proposal of reference values. The proposed values are indicated as a limit of maximum value or a range (minimum–maximum) for the essential elements. Values are expressed as µg element/g hair, equivalent to mg/kg as indicated by some authors. Data in the first column are values based on the simple criteria of 5–95th percentiles, eliminating the outliers’ values. The proposed values in the right column (in bold) are modulated considering, case by case, the distribution in our statistics and the other published data as discussed in Supplementary Table S5.
Based on 5–95 PC (µg/g)Proposal Considering Other Data (µg/g)
Be 0–0.004<0.01
B0.2–6.1<6
Na8–4403–600
Mg0.4–264–100
Al0.6–26<20
K2–922–200
Ca35–110020–1000
V0.01–0.28<0.4
Cr0.01–0.35<0.4
Mn0.01–0.40.03–1.1
Fe2–212–20
Co0.001–0.04<0.04
Ni0.02–0.7<0.70
Cu3–3305–100
Zn39–97040–850
As0.01–0.06<0.1
Se0.1–0.80.1–1
Sr0.3–13<5
Mo0.01–0.1<0.1
Cd0.003–0.11<0.1
Ba0.02–0.9<0.9
Hg0.08–9<3
Tl0–0.001<0.001
Pb0.09–5.2<2
Bi0.001–0.02<0.03
Ag0–0.01<0.01
Au0–0.2<0.2
U<0.001<0.001
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ruiz, R.; Estevan, C.; Estévez, J.; Alcaide, C.; Sogorb, M.A.; Vilanova, E. Reference Values on Children’s Hair for 28 Elements (Heavy Metals and Essential Elements) Based on a Pilot Study in a Representative Non-Contaminated Local Area. Int. J. Mol. Sci. 2023, 24, 8127. https://doi.org/10.3390/ijms24098127

AMA Style

Ruiz R, Estevan C, Estévez J, Alcaide C, Sogorb MA, Vilanova E. Reference Values on Children’s Hair for 28 Elements (Heavy Metals and Essential Elements) Based on a Pilot Study in a Representative Non-Contaminated Local Area. International Journal of Molecular Sciences. 2023; 24(9):8127. https://doi.org/10.3390/ijms24098127

Chicago/Turabian Style

Ruiz, Roberto, Carmen Estevan, Jorge Estévez, Carolina Alcaide, Miguel A. Sogorb, and Eugenio Vilanova. 2023. "Reference Values on Children’s Hair for 28 Elements (Heavy Metals and Essential Elements) Based on a Pilot Study in a Representative Non-Contaminated Local Area" International Journal of Molecular Sciences 24, no. 9: 8127. https://doi.org/10.3390/ijms24098127

APA Style

Ruiz, R., Estevan, C., Estévez, J., Alcaide, C., Sogorb, M. A., & Vilanova, E. (2023). Reference Values on Children’s Hair for 28 Elements (Heavy Metals and Essential Elements) Based on a Pilot Study in a Representative Non-Contaminated Local Area. International Journal of Molecular Sciences, 24(9), 8127. https://doi.org/10.3390/ijms24098127

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop