Determination of the Mineral Composition and Toxic Element Contents of Propolis by Near Infrared Spectroscopy

The potential of near infrared spectroscopy (NIR) with remote reflectance fiber-optic probes for determining the mineral composition of propolis was evaluated. This technology allows direct measurements without prior sample treatment. Ninety one samples of propolis were collected in Chile (Bio-Bio region) and Spain (Castilla-León and Galicia regions). The minerals measured were aluminum, calcium, iron, potassium, magnesium, phosphorus, and some potentially toxic trace elements such as zinc, chromium, nickel, copper and lead. The modified partial least squares (MPLS) regression method was used to develop the NIR calibration model. The determination coefficient (R2) and root mean square error of prediction (RMSEP) obtained for aluminum (0.79, 53), calcium (0.83, 94), iron (0.69, 134) potassium (0.95, 117), magnesium (0.70, 99), phosphorus (0.94, 24) zinc (0.87, 10) chromium (0.48, 0.6) nickel (0.52, 0.7) copper (0.64, 0.9) and lead (0.70, 2) in ppm. The results demonstrated that the capacity for prediction can be considered good for wide ranges of potassium, phosphorus and zinc concentrations, and acceptable for aluminum, calcium, magnesium, iron and lead. This indicated that the NIR method is comparable to chemical methods. The method is of interest in the rapid prediction of potentially toxic elements in propolis before consumption.

of lead content, using the Graphite Furnace Atomic Adsorption Spectrometry (GF AAS) method [24]. Thus, using flame atomic absorption spectrometry Formicki et al. [12] determined the levels of cadmium, iron, magnesium, nickel, lead and zinc in various bee products collected in Poland; Pierini et al. [13] used a electroanalytical method for quantification of lead in Argentine propolis; Serra-Bonvehí and Orantes-Bermejo [25] determined, among others, the levels of arsenic, cadmium, mercury and lead in samples of propolis collected in southern Spain by ICP-atomic emission spectrometry and Finger et al. [26], studied the content of cadmium, chromium and lead in Brazilian propolis from Paraná by electrothermal atomization and flame atomic absorption spectrometry after calcination in a muffle furnace. Moreover, Gong et al. [11] studied the relationships between the geographical origin and the content of calcium, aluminum, magnesium, potassium, iron, sodium, zinc, manganese, strontium, copper, chromium, nickel and toxic elements like arsenic, cadmium and lead determined by inductively coupled plasma atomic emission spectrometry after microwave digestion of Chinese propolis. However, there are few studies evaluating the potential of near infrared spectroscopy (NIR) for quantitative analysis of propolis. Visible/near infrared spectroscopy (Vis/NIRS) has been applied for the analysis of chrysin and galangin in Chinese propolis [27] and to explore its applicability for the determination of antioxidant properties [28].The detection of propolis falsification by the addition of flavonoid glycosides and tree latex has been carried out by Fourier Transform-NIR [29]. NIR has been also used to confirm the identity of isolated beeswax propolis [30]. Regarding the mineral characterization and to ensure the safety and quality of this product for marketing in local and international markets, this paper proposes a fast method for determining the mineral composition of propolis by using near infrared spectroscopy (NIR), because this analytical technique is rapid, non-destructive, and requires little or no sample preparation [31]. Therefore, in this paper, propolis samples from Chile and Galicia and Castilla-León (Spain) were analyzed with the aim of develop a rapid method of analysis to quantify some mineral and trace elements, using near-infrared spectroscopy with a remote reflectance fibre-optic probe applied directly to the sample.

Propolis Samples
Propolis samples (N = 91) have been directly collected by beekeepers in Chile (52 samples: Bio-Bio region) and Spain (39 samples: Galicia and Castilla-León regions), from both organic and conventional farms (the samples are all from different farms, and are not duplicated). Samples were collected mostly with mesh and employing the scraping technique, keeping them frozen until used in the laboratory. For the NIR analysis, of all the 91 samples, 71 were employed for the denominated calibration set, and the other 20 were used for the external validation.

Chemical Analysis of the Mineral Composition
The mineral composition of propolis was determined using Inductively Coupled Plasma Optical Emission Spectrometry (ICP-OES) for aluminum, calcium, iron, potassium, magnesium, sodium and potassium, and Inductively Coupled Plasma Mass Spectrometry spectrometry (ICP-MS) for zinc, chromium, nickel, copper and lead. Propolis samples were crushed in the laboratory before their analysis with a grinder (Knifetec 1095 Sample Mill, Foss Tecator, Hohnaa, Sweden). Prior to analysis of mineral elements, the mineralization of the samples (0.2 g) in a microwave system (Ethos Sel Milestone, Ontario, ON, Canada) was performed and subsequently introduced into a high pressure capsule. In a first phase 5 mL of HNO3was added and a power of 1000 watts was applied for 5 min. Once the sample was cool, another 5 mL of HNO3 and 1 mL of 30% H2O2 were added, applying a power of 1000 W for 10 min. The sample was cooled to room temperature, made up to 100 mL with distilled water and stored at 4 °C until analysis.
The ICP-OES determinations were carried out using Ultima 2 equipment (JobinYvon, Paris, NJ, USA), performing the calibration with certified standard solutions, ranging from 0.5 to 10 mg/kg. Detection limits were 0.1 mg/kg in solution. The ICP-MS determinations were performed on an Elan 6000 instrument (Perkin-Elmer, Wellesley, Massachusett, MA, USA). Do do this an internal standard (Sc, Y, Ho Ge 20 and 100 µg/kg) was added to an aliquot of previously prepared sample. The calibration was performed with certified standard solutions, also adding to the same internal standard for the calibration samples and adapted to a range of 10 to 200 µg/kg. Detection limits were 0.1 µg/kg in solution.

NIR Spectroscopy
A Foss NIRSystem 5000 (Hillerod, Denmark) with a standard 1.5 m 210/210 bundle fibre-optic probe (Ref. n° R6539-A) was used. Figure 1 shows a diagram of the NIRS system, where the main components (optical system, sample module, fiber optic probe) can be observed. The reflectance detectors receive radiation from diffuse scattering of the sample. These detectors (four PbS elements) are positioned at 45° from the surface of the sample (to minimize the specular reflectance). The used spectral range is 1100-2000 nm, since above 2000 nm significant signal attenuation occurs because of the strong absorption of the -OH groups that may be present in the optical fiber. The probe employs a remote reflectance system and uses a ceramic plate as reference. The window is of quartz, with a 5 cm × 5 cm surface area. The measurement of the spectra was carried out using NIRS technology and a remote reflectance fibre-optic probe that was applied directly to samples of crushed propolis. The spectra were recorded at intervals of 2 nm, performing 32 scans for both the reference and samples. To minimise sampling errors, all the samples were analysed in triplicate.

Principal Component Analysis
Principal component analysis (PCA) was carried with the NIR spectral data from samples of the calibration set using WinISI II version 1.50 (Intrasoft International, LLC, Silver Spring, Maryland, MD, USA). This analysis transforms the original variables (wavelengths) into new axes called principal components, which are orthogonal, so that the data set presented on these axes are uncorrelated with each other. Spectra were pretreated by different techniques, including multiplicative scatter correction (MSC), standard normal variate (SNV), DT (Detrend) or SNV-DT [32]. These treatments allow minimization of the scattering effect, since mainly the shift of the maximum and the width changes of the spectra were considered.

Modified Partial Least Squares Regression
The modified partial least squares (MPLS) regression method was used to obtain the NIR equations for the minerals quantified in propolis. The principal aim was generate models that allow the prediction of these components in the propolis matrix. Partial least squares (PLS) regression is similar to principal component regression (PCR), but uses both reference data and spectral information to form the factors useful for fitting purposes. To optimize the multivariate regression equations, the spectral scattering effects were taken into account with several mathematical treatments (MSC; SNV; DT or SNV-DT) [32]. MPLS is often more stable and accurate than the standard partial least squares algorithm (PLS). Calibrations were performed by modified partial least squares regression (MPLS) for each component, after removing the samples for spectral (criterion H) or chemical reasons (criterion T). The criterion H (Mahalanobis distance) explains the difference of the spectrum of an unknown sample from the mean spectrum of the set of samples. Samples with an H-value higher than 3 are considered as a different population and are eliminated. Furthermore, the risk of there being mistakes in the equations under practical conditions is very low or almost null when the standardised H statistic (Mahalanobis distance) is used during routine analysis of unknown samples. Using the T > 2.5 criterion, samples that were different from the population owing to the chemical criterion were removed from the set. In MPLS, the NIR residuals at each wavelength, obtained after each factor has been calculated, are standardized (dividing by the standard deviations of the residuals at each wavelength) before calculating the next factor [33]. In order to select the optimal number of factors and to avoid over fitting, cross-validation is recommended [34]. The calibration set is divided into several groups for the cross-validation. Each group is then validated using a calibration developed on the other samples. Validation errors generated are combined into a root mean square error of cross-validation RMSECV [35]. It has been reported that the RMSECV is the best single estimate for the prediction capability of the equation and that this statistic is similar to the average root mean square error of prediction (RMSEP) from 10 randomly-chosen prediction sets. In all cases, cross-validation was performed by splitting the population into six groups. Squared correlation coefficient was considered for the realization of models. The squared correlation coefficient for predicted versus measured compositions in cross-validation and the ratio of standard deviation (SD) to RMSECV of data set have been used. This ratio (of the SD to the RMSECV) is called the ratio of performance to deviation (RPD). This ratio is desired to be larger than 2 for a good calibration [36]. An RPD ratio less than 1.5 indicates poor predictions and the model cannot be used for further prediction. The statistics used to select the best equations were R 2 , determination coefficient; RMSEC, root mean square error of calibration; RMSECV, root mean square error of cross-validation [34].

Chemical Composition
The propolis samples had wide ranges in the mineral composition and high standard deviation values for the studied elements. Table 1 shows the mean concentration and the range of values of the mineral composition (aluminum, calcium, iron, potassium, magnesium, sodium, phosphorus, zinc, lead, chromium, nickel and copper) of all the propolis samples according to the geographical origin obtained by the chemical reference methods (ICP-OES and ICP-MS). Also, a comparative study of the mineral composition of propolis from other geographical origins is shown in Table 2.
The possibility of correlation between all elements was investigated. Potassium, the most abundant mineral compound, was significantly (p < 0.01) and negatively correlated with aluminum (r = −0.29) and iron (r = −0.32), and positively correlated with phosphorus (r = 0.28), lead (r = 0.33) and copper (p < 0.05, r = 0.26). Iron was positively correlated (p < 0.01) with aluminum (r = 0.57), as observed for potassium, and also with calcium (r = 0.36), magnesium (r = 0.35) and zinc (r = 0.31). Magnesium showed a strong correlation with calcium (r = 0.97) and a significant correlation (p < 0.01) with phosphorus (r = 0.31). Finally, regarding the potentially toxic elements beside the significant correlation between lead and potassium, a strong correlation was found between chromium and nickel (r = 0.90).This result differs from previously reported by Finger et al. [26] who found a significant correlation between calcium and potassium and from those reported by Formicki et al. [12] who found a significant but negative correlation between iron and magnesium. These two works also found a positive correlation between cadmium and lead that they attributed to a common polluting source. However in this work the cadmium was not analyzed but a strong correlation between chromium and nickel was found, metals that were not analyzed in those works.

NIR Calibration Equations
For the calibration equations, a set of 71 samples of propolis from different sources (Galicia, Castilla-León, Chile) was used. Initially, the principal component analysis (PCA) was performed. In all cases, the spectral variability explained above was 99%. Using both criteria (criterion H and T) six samples were deleted for the calibration of aluminum, 11 for the calcium, six for iron, 13 for potassium, seven for magnesium, five for phosphorus, 10 for chromium, 13 for copper, four for nickel, 12 for lead and seven for zinc. Calibrations were performed by modified partial least squares regression (MPLS), using the spectral data and chemical data matrix (obtained with ICP-OES, ICP-MS). The best of the different mathematical treatments, the concentration range, standard deviation and the calibration parameters for all elements are shown in Table 3. The quantification of the chemical elements using NIR technology is possible thanks to the distinct associations of these elements with the organic and inorganic material and with the water molecules, because of relative strong absorption of the overtones and combination modes of OH, sulphates and carbonate groups. Since 1981, there have been several reports on mineral elements in plants determined using NIR [41][42][43][44]; if NIRS can be used for determining mineral concentrations this is due to the association between minerals and organic functional groups or the organic matrix [45]. When NIR radiation is absorbed by molecules the energy is converted to molecular vibration energy. Molecules which are infrared (IR)-active are those which undergo a change in the dipole moment during transition, this means that bonds commonly found in biological systems such as C-H, O-H and N-H bonds are IR active. The prediction of trace elements by NIRS in agricultural products has been reported even less frequently and has always been used in the context of plants, only a few reports were found related to the use of NIRS for macro and trace minerals in both grasses and hay samples [46], in botanical fractions of semi-arid grassland [47] or legumes [48]. In the case of the propolis, the most likely association of minerals with organic matter is through flavonoids. Propolis presents a general composition of 40%-50% resins and balsams (that contain, in turn, 50% flavonoids and phenolic acids), 30%-40% wax, 10% volatile oils, 5% pollen and 5% minerals and other organic compounds. The chemical structures of flavonoids favour the formation of very stable complexes with heavy metals, making propolis a hyperaccumulator of these kinds of elements [18][19][20]. The correlation between the concentration and that measured at different wavelengths is given by the expression: y = β0 + β1 X λ1 + β2 X λ2 + β3 X λ3 +…+ βn X λv , where β are the coefficients and X λ1 , X λ2 , X λ3 ,… X λn , are the wavelengths at which the correlation of the concentration of the components is maximum (in + or − value). The higher values of those β coefficients for each of the parameters studied (Al, Ca, Fe, K, Mg, P, Cr, Cu, Ni, Pb, Zn) together with the wavelength where these coefficients showed their maxima and The results showed that it is possible to determine the composition of aluminum, calcium, iron, potassium, magnesium, phosphorus, chromium, copper, nickel, lead and zinc in the ranges indicated in samples of propolis of different origins (from Chile and two Spanish regions) by direct application of a NIR fiber-optic probe on crushed propolis samples without prior treatment or manipulation.

Internal Validation (Prediction)
Models' evaluations were performed by cross-validation. In this method, the set of calibration samples is divided into a series of subsets, in our case six. Of these, five were taken for the calibration set and one for the prediction set. The process is repeated as many times as there are sets, so that all pass through the calibration set and the prediction set. Using this process, we validated the models used and checked their prediction capacities. Figure 2 shows the correlation of the values obtained in the laboratory with respect to those predicted by NIR with remote reflectance fibre-optic probe for aluminum, calcium, iron, potassium, magnesium, phosphorus, zinc, chromium, nickel, copper and lead in propolis. The prediction capacity of the model obtained was evaluated with the ratio performance deviation (RPD) [36].This parameter is defined as the relationship between the standard deviation of the chemical method (SD ref) and RMSECV, root mean square error of cross-validation encountered in the NIRS model. Table 3 shows the values obtained for RPD parameter that were comprehended between 3.1 for calcium and 1.0 for nickel. Therefore, NIRS technology presents a capacity for prediction that is interesting for the determination of mineral composition in samples of unknown propolis.   The proposed NIR method was verified by applying the developed chemometric model to 20 new samples of different compositions (called external validation set samples). The recording of the spectra in triplicate and the average spectra was considered. Then, calibration equations obtained in this work were applied and predicted values were compared with reference data for aluminum, calcium, chromium, copper, iron, potassium, magnesium, nickel, phosphorus, lead and zinc determined by OES-ICP and ICP-MS spectrometry. Table 4 shows the results of external validation corresponding with the prediction of mineral composition of 20 independent samples.

Conclusions
NIR methodology with a remote reflectance fibre-optic probe is presented as an effective method of analysis for determining some mineral and toxic trace elements in crushed propolis. This method is applicable to samples with a wide range of contents of aluminum, calcium, iron, potassium, magnesium, phosphorus, zinc, chromium, nickel, copper and lead. This is the first work on the quantification of mineral elements in a propolis matrix with NIR technology. The method is particularly interesting for the prediction of potentially toxic elements such as zinc, copper and lead, since it allows a previous detection in a short time (3 or 4 min). This ensures the safety and quality of propolis for commercialization in national and international markets.