Electrochemical Determination of the “Furanic Index” in Honey

5-(hydroxymethyl)furan-2-carbaldehyde, better known as hydroxymethylfurfural (HMF), is a well-known freshness parameter of honey: although mostly absent in fresh samples, its concentration tends to increase naturally with aging. However, high quantities of HMF are also found in fresh but adulterated samples or honey subjected to thermal or photochemical stresses. In addition, HMF deserves further consideration due to its potential toxic effects on human health. The processes at the origin of HMF formation in honey and in other foods, containing saccharides and proteins—mainly non-enzymatic browning reactions—can also produce other furanic compounds. Among others, 2-furaldehyde (2F) and 2-furoic acid (2FA) are the most abundant in honey, but also their isomers (i.e., 3-furaldehyde, 3F, and 3-furoic acid, 3FA) have been found in it, although in small quantities. A preliminary characterization of HMF, 2F, 2FA, 3F, and 3FA by cyclic voltammetry (CV) led to hypothesizing the possibility of a comprehensive quantitative determination of all these compounds using a simple and accurate square wave voltammetry (SWV) method. Therefore, a new parameter able to provide indications on quality of honey, named “Furanic Index” (FI), was proposed in this contribution, which is based on the simultaneous reduction of all analytes on an Hg electrode to ca. −1.50 V vs. Saturated Calomel Electrode (SCE). The proposed method, validated, and tested on 10 samples of honeys of different botanical origin and age, is fast and accurate, and, in the case of strawberry tree honey (Arbutus unedo), it highlighted the contribution to the FI of the homogentisic acid (HA), i.e., the chemical marker of the floral origin of this honey, which was quantitatively reduced in the working conditions. Excellent agreement between the SWV and Reverse-Phase High-Performance Liquid Chromatography (RP-HPLC) data was observed in all samples considered.


Introduction
Honey is the sweet natural product that honeybees obtain by conversion of the nectar gathered from flowers. Its composition depends on the geographical, floral, and entomological origins and includes water (15-20%), sugars (80-85% w/w), nitrogenous compounds, enzymes, phenolic and volatile compounds, organic and amino acids, minerals, and vitamins. On the other hand, it is well known that seasonal, environmental, storage, and processing time and conditions may affect honey's composition. Although the nutritional and healing features of this food are appreciated, the possible presence of heavy metals (usually in traces) [1], some alkaloids [2], and reactive organic compounds, such as aldehydes, might imply a possible threat for the health of consumers. The 5-(hydroxymethyl)furan-2-carbaldehyde, better known as hydroxymethylfurfural (HMF), is produced during the

Electrochemical Behavior of HMF
The electrochemical behavior of HMF was studied in cyclic voltammetry by varying the nature of the WE (Pt, GC, Au, Hg), the solvent system (solvents: water or methanol; supporting electrolytes: NaH2PO4, Na2HPO4, LiCl, CH3COONa, NaHCO3, LiClO4, [(C4H9)4N]ClO4) and the potential scan rate (between 0.05 and 2.0 V s −1 ). Based on the cathodic nature of the responses of the analyte, a preliminary evaluation has allowed to discard both Pt and Au working electrodes. Hence, only Hg and GC electrodes were further considered in the optimization of the experimental conditions, as summarized below in Table 1. HMF solution 10 cm 3 , 1 mmol dm −3 The CV response of an aqueous 1 mmol dm −3 solution of HMF showed a well-defined and reproducible irreversible cathodic peak at potentials between −1.10 and −1.55 V vs. SCE, as a function of the nature of both the WE and of the supporting electrolytes (Table 2 and Figure 2). The electrochemical behavior of HMF was studied in cyclic voltammetry by varying the nature of the WE (Pt, GC, Au, Hg), the solvent system (solvents: water or methanol; supporting electrolytes: NaH 2 PO 4 , Na 2 HPO 4 , LiCl, CH 3 COONa, NaHCO 3 , LiClO 4 , [(C 4 H 9 ) 4 N]ClO 4 ) and the potential scan rate (between 0.05 and 2.0 V s −1 ). Based on the cathodic nature of the responses of the analyte, a preliminary evaluation has allowed to discard both Pt and Au working electrodes. Hence, only Hg and GC electrodes were further considered in the optimization of the experimental conditions, as summarized below in Table 1. The CV response of an aqueous 1 mmol dm −3 solution of HMF showed a well-defined and reproducible irreversible cathodic peak at potentials between −1.10 and −1.55 V vs. SCE, as a function of the nature of both the WE and of the supporting electrolytes (Table 2 and Figure 2). Table 2. Cathodic peak potential (E p,c ) and current density at the cathodic peak (J p,c ) of the irreversible reduction process of HMF by varying the nature of the WE (Hg or GC) and the nature and the pH of the solvent system (a 0.1 mol dm −3 water solution of five different supporting electrolytes). HMF concentration: 1 mmol dm −3 , AE: Pt, RE: SCE, potential scan rate: 0.10 V s −1 . No anodic peak was observed in anodic direct scan (WE = GC) or after the inversion of the scan direction after crossing the cathodic peak ( Figure 2a). Voltammograms acquired using Hg as WE at scan rates between 0.10 and 1.5 V s −1 showed a cathodic shift of Ep,c (from −1.48 to −1.54 V vs SCE) and a linear increase of the peak's current (from −3·10 −5 to −1·10 −4 A) (Figure 2b). No backward anodic peak has been observed in this range of scan rates. The linear decreasing of the ip,c/√(logVscan) ratio as a function of the cathodic shift of the Ep,c peak suggests a kinetic control on the overall electrochemical process [35]. The replacement of water with methanol as solvent did not cause meaningful changes in the morphology of the main cathodic peak previously observed. Indeed, in these conditions, the CV responses recorded on GC, in cathodic direct scan, showed the same irreversible peak (peak 1 in Figure 3) at about −1.50 V vs. SCE. When the scan direction was inverted (from −2.0 to 2.0 V), three ill-resolved and consecutive irreversible anodic peaks became evident ( Figure 3, peaks 2, 3 and 4, respectively, at Ep,a between 0.90 and 1.70 V vs. SCE). A further inversion of the direction of the potential scan, performed just after peak 4, caused the appearance of a low and irreversible cathodic peak at Ep,c of −1.00 V (peak 5 in Figure 3), which was never observed in the cathodic direct scan ( Figure  3). No anodic peak was observed in anodic direct scan (WE = GC) or after the inversion of the scan direction after crossing the cathodic peak ( Figure 2a). Voltammograms acquired using Hg as WE at scan rates between 0.10 and 1.5 V s −1 showed a cathodic shift of E p,c (from −1.48 to −1.54 V vs SCE) and a linear increase of the peak's current (from −3 × 10 −5 to −1 × 10 −4 A) (Figure 2b). No backward anodic peak has been observed in this range of scan rates. The linear decreasing of the i p,c / √ (logV scan ) ratio as a function of the cathodic shift of the E p,c peak suggests a kinetic control on the overall electrochemical process [35].

WE
The replacement of water with methanol as solvent did not cause meaningful changes in the morphology of the main cathodic peak previously observed. Indeed, in these conditions, the CV responses recorded on GC, in cathodic direct scan, showed the same irreversible peak (peak 1 in Figure 3) at about −1.50 V vs. SCE. When the scan direction was inverted (from −2.0 to 2.0 V), three ill-resolved and consecutive irreversible anodic peaks became evident ( Figure 3, peaks 2, 3 and 4, respectively, at E p,a between 0.90 and 1.70 V vs. SCE). A further inversion of the direction of the potential scan, performed just after peak 4, caused the appearance of a low and irreversible cathodic peak at E p,c of −1.00 V (peak 5 in Figure 3), which was never observed in the cathodic direct scan ( Figure 3).
The voltammetric responses of the solutions of HMF were in agreement with those reported by other authors who studied this analyte's behavior in aqueous solutions [25] as well as the reduction processes that involve the carbonylic group in 2F [35,36]. All the authors attributed the irreversible cathodic peak at about −1.50 V to the two-electron reduction of the aldehydic group to the corresponding alcohol, that, in the case of HMF, produces the 2,5-dihydroxymethylfurane. In addition, the meaningful cathodic shift of this irreversible reduction peak, observed passing by a solvent system at pH = 4.1 (supporting electrolyte: NaH 2 PO 4 ) to that at pH = 9.8 (supporting electrolyte: Na 2 HPO 4 ), is consistent with the observations of Ganesan et al. [36], because low pH values assist the reduction process. On the other hand, the most anodic peak, observed at 1.68 V vs. SCE in the direct anodic scan on GC in the methanol/LiClO 4 solvent system, was previously attributed to the irreversible oxidation of HMF to 2,5-diformylfurane (1.63 V vs. SCE, [25]), which was likely reducible, in the reverse cathodic scan, to HMF (cathodic peak at −1.00 V vs. SCE). Since the reduction of carbonylic group of HMF is the main (often the unique) peak and the most reproducible one present in the different solutions, we decided to further investigate it to find the optimal conditions for its quantitative determination in honey samples.  The voltammetric responses of the solutions of HMF were in agreement with those reported by other authors who studied this analyte's behavior in aqueous solutions [25] as well as the reduction processes that involve the carbonylic group in 2F [35,36]. All the authors attributed the irreversible cathodic peak at about −1.50 V to the two-electron reduction of the aldehydic group to the corresponding alcohol, that, in the case of HMF, produces the 2,5-dihydroxymethylfurane. In addition, the meaningful cathodic shift of this irreversible reduction peak, observed passing by a solvent system at pH = 4.1 (supporting electrolyte: NaH2PO4) to that at pH = 9.8 (supporting electrolyte: Na2HPO4), is consistent with the observations of Ganesan et al. [36], because low pH values assist the reduction process. On the other hand, the most anodic peak, observed at 1.68 V vs. SCE in the direct anodic scan on GC in the methanol/LiClO4 solvent system, was previously attributed to the irreversible oxidation of HMF to 2,5-diformylfurane (1.63 V vs. SCE, [25]), which was likely reducible, in the reverse cathodic scan, to HMF (cathodic peak at −1.00 V vs. SCE). Since the reduction of carbonylic group of HMF is the main (often the unique) peak and the most reproducible one present in the different solutions, we decided to further investigate it to find the optimal conditions for its quantitative determination in honey samples.

Electrochemical Behavior of Furanic Aldehydes and Acids
The electrochemical behavior of 2F, 3F, 2FA, and 3FA was studied in CV in the same conditions used for HMF. Table 3 reports the potentials of the cathodic peak measured using Hg and GC as working electrodes, and LiCl, LiClO4, CH3COONa, and Na2HPO4 as supporting electrolytes in water. All the species showed an irreversible cathodic reduction peak at potential between −1.46 and −1.84 V vs. SCE, using both the working electrodes, when LiCl and LiClO4 were used as supporting electrolytes. These peaks can be tentatively ascribed to the reductions of 2F and 3F to give the corresponding alcohols (i.e., 2and 3-hydroxymethylfurane, respectively) and the reductions of 2FA and 3FA toward 2F and 3F, respectively. On the other hand, the cathodic shift of the reduction process caused

Electrochemical Behavior of Furanic Aldehydes and Acids
The electrochemical behavior of 2F, 3F, 2FA, and 3FA was studied in CV in the same conditions used for HMF. Table 3 reports the potentials of the cathodic peak measured using Hg and GC as working electrodes, and LiCl, LiClO 4 , CH 3 COONa, and Na 2 HPO 4 as supporting electrolytes in water. All the species showed an irreversible cathodic reduction peak at potential between −1.46 and −1.84 V vs. SCE, using both the working electrodes, when LiCl and LiClO 4 were used as supporting electrolytes. These peaks can be tentatively ascribed to the reductions of 2F and 3F to give the corresponding alcohols (i.e., 2-and 3-hydroxymethylfurane, respectively) and the reductions of 2FA and 3FA toward 2F and 3F, respectively. On the other hand, the cathodic shift of the reduction process caused by the use of alkaline supporting electrolytes brought the E p,c of 2FA and 3FA out of the available cathodic window. This happened when the GC working electrode was used in the presence of CH 3 COONa as supporting electrolyte and, for both Hg and GC electrodes, when Na 2 HPO 4 was used. Hence, the possibility to measure at the same time all these analytes in honey is achievable only using neutral supporting electrolytes such as LiCl or LiClO 4 . Preliminary CV experiments allowed to ascertain the feasibility in a contemporary determination of HMF and of the principal furanic aldehydes and acids in honey. Unfortunately, CV is an unsuitable electrochemical technique for an accurate and sensitive determination of trace compounds. On the other hand, the amount of the most abundant of these analytes in fresh honey, the HMF, is typically of few tens of mg kg −1 , i.e., well below the concentration used for the preliminary CV tests. Hence, the quantitative method to be assessed had to be able to reach LoQ levels at least of few mg kg −1 , and the best choice to quantify trace amounts of electroactive organic analytes was the square wave voltammetry (SWV). The method has been optimized as a function of the general parameters, such as the nature of the working electrode (Hg or GC) and of the supporting electrolyte (LiCl or LiClO 4 ), as well as of specific SWV parameters, such as the pulse height (values between 2.5 and 250 mV), the frequency (between 1.5 and 150 Hz), and the modulation amplitude (between 1 and 40 mV). The optimized SWV parameters are pulse height, 4 mV; frequency, 15 Hz; modulation amplitude, 25 mV. Table 4 shows the optimization of both WE and supporting electrolyte in the development of the SWV method. HMF was chosen as a molecular benchmark, since it is the most abundant furanic compound present in honey. Despite the higher sensitivity of the GC electrode, when LiClO 4 was the supporting electrolyte, no response was obtained using LiCl. On the other hand, the Hg electrode could detect the benchmark HMF content using both the supporting electrolytes. In these cases, the density of current measured in the water/LiCl solvent system was roughly double the amount measured in water/LiClO 4 . Based on the results here obtained, the SWV method was tested also in the presence of 2F, 3F, 2FA, and 3FA. Additionally in this case, two WEs and two supporting electrolytes were used in these experiments. Table 5 reports the electrochemical parameters obtained for all the analytes in this cycle of measurements.
Data reported in Table 5 confirmed the preliminary results obtained for HMF. The E p,c values for the five analytes measured with Hg electrode were more reproducible and closer to each other in LiCl rather than in LiClO 4 . In addition, it was impossible to detect the cathodic peak of 3FA with GC electrode, when LiClO 4 was used as supporting electrolyte. Furthermore, the current densities measured for all analytes, when Hg was the working electrode, were much higher than those measured by GC electrode, and the same happened switching from LiClO 4 to LiCl as the supporting electrolyte. Hence, Hg as the working electrode and LiCl as the supporting electrolyte were the best choices for the SWV determination of HMF, 2F, 3F, 2FA, and 3FA in honey. If the extreme closeness of all the E p,c values made it almost impossible to discriminate the amount of each analyte, individually, as clearly demonstrated by Figure 4, it was however possible-and easy-to quantify the amounts of all these analytes at once, expressing them in terms of an equivalent amount of HMF. This observation allowed defining the "Furanic Index" (FI) as a freshness-related parameter for honey that is able to evaluate the overall amount of HMF, 2F, 3F, 2FA, and 3FA, as expressed in mg of HMF per kg of honey. From a quantitative viewpoint, the expression of FI is: FI (mg dm −3 ) = C HMF + C 2F ·(HMF/2F) + C 3F ·(HMF/3F) + C 2FA ·(HMF/2FA) + C 3FA ·(HMF/3FA) where C analyte is the amount (in mg dm −3 ) of each furanic compound, and the HMF/analyte ratio refers to the relevant molecular weights. Table 5. Cathodic peak potentials E p,c (V vs. SCE) and cathodic peak current densities J p,c (µA mm −2 ) in the SWV determination of HMF, 2F, 3F, 2FA, and 3FA (concentration: 1 mmol dm −3 each) by varying the nature of WE (Hg or GC) and that of the supporting electrolyte (LiCl or LiClO 4 , both 0.1 mol dm −3 in water). Potential scan rate: 0.10 V s −1 ; cathodic window of potential: between 0 and −2.0 V vs. SCE; AE: Pt; RE: SCE; pulse height: 4 mV; frequency: 15 Hz; modulation amplitude: 25 mV.

Validation
Validation of the proposed method has been accomplished in terms of LoD, LoQ, linearity, precision, and trueness. Table 6 reports the features describing the performances of the method proposed.

LoD and LoQ
These parameters were calculated, according to IUPAC recommendation, through the ULA1 method [37]. Four different 0.1 mol dm −3 LiCl aqueous solutions containing HMF in concentrations of 3 mg dm −3 , 5 mg dm −3 , 7 mg dm −3 , and 10 mg dm −3 , respectively, were prepared and analyzed by means of SWV in triplicate. The LoD obtained in this way is 0.6 mg dm −3 . This amount is significantly higher than that (i.e., 0.0021 mg dm −3 ) recently

Validation
Validation of the proposed method has been accomplished in terms of LoD, LoQ, linearity, precision, and trueness. Table 6 reports the features describing the performances of the method proposed.

LoD and LoQ
These parameters were calculated, according to IUPAC recommendation, through the ULA1 method [37]. Four different 0.1 mol dm −3 LiCl aqueous solutions containing HMF in concentrations of 3 mg dm −3 , 5 mg dm −3 , 7 mg dm −3 , and 10 mg dm −3 , respectively, were prepared and analyzed by means of SWV in triplicate. The LoD obtained in this way is 0.6 mg dm −3 . This amount is significantly higher than that (i.e., 0.0021 mg dm −3 ) recently measured by Sahli et al. [25], but in this case, the approach used for the calculation was not described. On the other hand, it is well known that in the past, it has been possible to associate to the general expression "detection limit" values spanning over more than three orders of magnitude [38]. Hence, even without arguing on the overall reliability of very low LoDs reported in the literature, it is well known that the application of the ULA methods, albeit providing sometimes LoDs higher than those obtained by more "optimistic" approaches [33], is one of the best choices to avoid both type-1 and type-2 decision errors. For the ULA1 approach, LoQ is three times the LoD value, i.e., 1.8 mg dm −3 .

Linearity
Linearity was measured over a concentration range between LoQ and 200 mg dm −3 . This is the typical interval of amounts observed for HMF in honeys of different geographical origin and age. As shown in Table 6, very good correlation coefficients (R 2 ) were observed in this case. In addition, no "hidden" deviations from linearity have been evidenced by the output of the graphical analysis performed on the residuals of the regression line. The LoD value was calculated according to [37]. 1 LoD calculated on a 50% solution of honey in a 0.1 mol dm −3 LiCl/water solvent system; 2 LoD calculated on pure honey; 3 parameter evaluated by analyzing four different honey samples five times in the same analytical session; 4 parameter evaluated by analyzing five times in five different analytical sessions over five months a honey sample exhibiting an intermediate FI with respect the concentration range measured; 5 evaluated by means of comparison with an independent analytical method (i.e., the RP-HPLC method, [33]) made on five different honey samples.

Precision
This parameter was evaluated on real samples in terms of repeatability and intermediate precision, which are both expressed as the percent relative standard deviation (RSD%). Repeatability was obtained by analyzing five times, in the same analytical session, four honey samples with a concentration range between 4.8 and 306 mg kg −1 . On the other hand, intermediate precision was evaluated by analyzing five times, in five different analytical sessions over five months, a honey sample showing an intermediate value of FI. As a function of the FI in the honey samples, the RSD% of repeatability ranged between 0.72% (FI of 306 mg kg −1 ) and 9.9% (FI of 4.8 mg kg −1 ), whereas the RSD% of the intermediate precision was of 3.0%, which was measured on a honey showing an FI of 160 mg kg −1 . All precision parameters were acceptable according to the Horwitz's theory [39].

Trueness
Given the absence of any Certified Reference Materials, but in consideration that an independent analytical method is currently available [33], trueness was evaluated through the comparison of results obtained from five real samples by both the proposed method and the RP-HPLC method previously published by this research group, which is able to singularly quantify HMF, 2F, 3F, 2FA, and 3FA in the same chromatographic run. The comparison of the experimental t (between 0.55 and 2.19) with the tabulated t (2.57 for 5 degrees of freedom, p = 0.95 in a two-tail t-test) allowed concluding that the proposed method was bias-free.

SWV Determination of the Furanic Index in Honeys 2.4.1. Method Application
The proposed method was tested on ten honey samples of different botanical origin (three from blossom, one from eucalyptus, one from thistle, and five from strawberry tree) and of different age (between 2014 and 2019). Figure 5 reports a typical SW voltammogram for one real sample.  Table 7 reports the amounts of FI measured in these samples. For comparison purposes, the last column of Table 7 also reported the sum of all analytes (i.e., HMF, 2F, 3F, 2FA, and 3FA) measured in the same samples by means of a literature RP-HPLC method [33], which are all expressed in equivalent amounts of HMF. It is possible to observe that the FI measured on blossom, eucalyptus, and thistle honey is never statistically different (two tails t-test, p = 0.95) from the amount measured by the RP-HPLC method. This fact is relevant, as it shows that the quicker SWV determination of the total amount of the furanic compounds is effective in the assessment of the freshness of the honey samples. In particular, RP-HPLC data substantiated the expected large predominance of the amount of HMF on those of the other analytes. As a matter of fact, HMF constitutes on average 93% of the total furanic species in these honeys. 2F and 2FA were sometimes quantified, at concentrations never higher than 3.0 mg kg −1 and 4.0 mg kg −1 , respectively, while 3F and 3FA were always below the relevant LoQs (i.e., 0.3 mg kg −1 and 0.09 mg kg −1 , respectively). The furanic indexes measured in these samples were always congruent to a recently produced and well-preserved honey: indeed, the measured FI ranged between 4.8 ± 0.4 mg kg −1 (B3 honey) and 22.0 ± 0.5 mg kg −1 (E honey).  Table 7 reports the amounts of FI measured in these samples. For comparison purposes, the last column of Table 7 also reported the sum of all analytes (i.e., HMF, 2F, 3F, 2FA, and 3FA) measured in the same samples by means of a literature RP-HPLC method [33], which are all expressed in equivalent amounts of HMF. It is possible to observe that the FI measured on blossom, eucalyptus, and thistle honey is never statistically different (two tails t-test, p = 0.95) from the amount measured by the RP-HPLC method. This fact is relevant, as it shows that the quicker SWV determination of the total amount of the furanic compounds is effective in the assessment of the freshness of the honey samples. In particular, RP-HPLC data substantiated the expected large predominance of the amount of HMF on those of the other analytes. As a matter of fact, HMF constitutes on average 93% of the total furanic species in these honeys. 2F and 2FA were sometimes quantified, at concentrations never higher than 3.0 mg kg −1 and 4.0 mg kg −1 , respectively, while 3F and 3FA were always below the relevant LoQs (i.e., 0.3 mg kg −1 and 0.09 mg kg −1 , respectively). The furanic indexes measured in these samples were always congruent to a recently produced and well-preserved honey: indeed, the measured FI ranged between 4.8 ± 0.4 mg kg −1 (B3 honey) and 22.0 ± 0.5 mg kg −1 (E honey). Table 7. Furanic index (FI) in ten honeys of different floral origin and age. Comparison with RP-HPLC data.

Sample
Year On the other hand, a very different situation was exhibited by the five strawberry tree honeys: in comparison to the RP-HPLC data, the amount of the FI is much higher (from 191% to 350%) than the sum of furanic species measured by RP-HPLC. The fact that this result only occurred for the samples of this botanical origin-irrespectively of the age of the honey led to supposing that in this honey, a compound normally absent in the other honeys was present, which was reduced at the same potential of all the analytes considered. As a matter of fact, the strawberry tree honey is rich in 2,5-dihydroxyphenylacetic acid, homogentisic acid (HA), and L-tyrosine catabolite [40], and this compound has been found until now only in this honey [41,42]. The electrochemical behavior of HA was studied in the same conditions used for the characterization of all the furanic compounds, and the voltammetric evidences acquired both in CV and in SWV confirmed that it is irreversibly reduced in the cathodic direct scan, at −1.48 V vs. SCE. Figure 6 shows a SW voltammogram of a solution 0.1 mol dm −3 of HA.  On the other hand, a very different situation was exhibited by the five strawberry tree honeys: in comparison to the RP-HPLC data, the amount of the FI is much higher (from 191% to 350%) than the sum of furanic species measured by RP-HPLC. The fact that this result only occurred for the samples of this botanical origin-irrespectively of the age of the honey led to supposing that in this honey, a compound normally absent in the other honeys was present, which was reduced at the same potential of all the analytes considered. As a matter of fact, the strawberry tree honey is rich in 2,5-dihydroxyphenylacetic acid, homogentisic acid (HA), and L-tyrosine catabolite [40], and this compound has been found until now only in this honey [41,42]. The electrochemical behavior of HA was studied in the same conditions used for the characterization of all the furanic compounds, and the voltammetric evidences acquired both in CV and in SWV confirmed that it is irreversibly reduced in the cathodic direct scan, at −1.48 V vs. SCE. Figure 6 shows a SW voltammogram of a solution 0.1 mol dm −3 of HA.  To verify the contribution of this molecule to the FI measured in strawberry tree honeys, HA was quantified in all the samples of this botanical origin by means of another RP-HPLC literature method [42]. The amounts of HA ranged between 92.5 ± 0.2 mg kg −1 and 241.2 ± 0.5 mg kg −1 , as shown in Table 8. This table also reports the theoretical amount of reducible species present in the strawberry tree honey, as measured by RP-HPLC, and expressed in equivalent HMF concentration, column RP-HPLC ΣHMF+Fs+FAs+HA. It is interesting to note that-for honey samples S1 and S3-this amount is not statistically different from that measured by SWV (criteria: two tails t-test, p = 0.95), whereas for the remaining samples, these data exhibit a slight bias (between −3.16%, sample S5 and −6.43%, sample S4) in comparison to the RP-HPLC one. However, this bias is well within the acceptability range for this level of concentration, as suggested by the AOAC International in its Peer Verified Methods Programs [43]. All data are expressed as average ± standard deviation, n = 3. 1 Data obtained according [42]; 2 sum of the concentrations of HMF, 2F, 3F, 2FA, and 3FA, obtained according to [33] and expressed as mg of HMF on kg of honey; 3 sum of the concentrations of HMF, 2F, 3F, 2FA, 3FA, and HA, obtained according to [33,42] and expressed as mg of HMF on kg of honey.

Honey Samples
The study was carried out on 10 honey samples of different botanical origin: strawberry tree (Arbutus unedo), eucalyptus (Eucalyptus camaldulensis), thistle (Galactites tomentosus), and blossom. Some of these honeys were directly provided by beekeepers; the others were commercial samples. All the samples were produced in Italy, almost all in the Sardinia Island, between 2014 and 2019. The floral origin was established on the basis of the indications of the producers and verified by melissopalynological analysis [44]. All samples were stored in the dark at 4 • C until SWV analysis.

SWV Determination of FI in Honey Samples
Ca. 5 g (exactly weighted on an analytical balance) of homogenized honey was diluted to 10 cm 3 by a 0.1 mol dm −3 LiCl aqueous solution and transferred in a polarographic cell in glass equipped with a four-hole gas-tight cover. The solution was carefully deaerated with Ar for at least 15 min; hence, a gentle flow of Ar ensured an inert blanket over the solution during the tests. Quantification was performed by internal calibration, by means of the multiple additions method (i.e., with three consecutive spikings of HMF, respectively 50%, 100%, and 150% of the amount of the furanic index evaluated by measurements performed by a literature HPLC method [33]). Each analytical evaluation was replicated at least three times.

Equipment
All the electrochemical measurements were performed using a CHI-650 electrochemical workstation (CH Instruments, Austin, TX, USA) interfaced with a computer with the specific software CHI-650, in a conventional three-electrodes voltammetric cell in glass, equipped with a gas-tight four-hole cover, positioned inside a Faraday cage. Four different working electrodes (WE) were used in this study: a 2 mm diameter glassy carbon (GC) disk, a 2 mm diameter Pt disk, a ca. 1 mm diameter Au hemisphere electrode, and a ca. 1 mm diameter Hg/Au amalgam hemisphere electrode, respectively. A 5 mm Pt disk was always the auxiliary electrode (AE), and a saturated calomel electrode (SCE) was the reference electrode (RE); all potential values given were referred. Before use, all of the electrode surfaces of WE and AE were polished with alumina powder (0.3 µm of diameter) and rinsed with ultrapure water. All the experiments were performed at room temperature in an Ar saturated solution.
The HPLC equipment used for confirmatory purposes was a Series 200 apparatus (Perkin Elmer, Milan, Italy) formed by a binary pump, by an UV-vis variable wavelength detector, and by a sampling valve equipped with a 20 mm 3 sample loop. Separation was accomplished on an Alltima C18 column 250 mm × 4.6 mm, 5 µm particle size (Alltech, Sedriano, Italy) fitted with a guard cartridge packed with the same stationary phase. Data were processed by a Turbochrom Workstation Software (Perkin Elmer, Milan, Italy).
An Ultra-turrax mixer model T18 (IKA, Staufen, Germany) was used to homogenize honey samples before analysis.

Conclusions
Although HMF has been a well-recognized freshness index of honeys for years, it is possible that the aging and/or imperfect storage conditions allow the formation also of different furanic compounds, such as furaldehydes and furoic acids. In order to provide a thorough evaluation of any reactive furanic-based compounds formed in honeys by degradation of this matrix, a rapid and accurate SWV method has been developed, optimized, and validated to measure the amount of HMF as well as of 2-furaldehyde, 2F, 3-furaldehyde, 3F, 2-furoic acid, 2FA, and 3-furoic acid, 3FA. All these analytes underwent an irreversible reduction on Hg electrode at potentials close to −1.50 V vs. SCE. Hence, a new chemical index (i.e., the so-called "Furanic Index", FI) was proposed for a quick and accurate assessment of the freshness of honeys. FI expresses, as mg of HMF per kg, the total amount of HMF, 2F, 3F, 2FA, and 3FA found in honey. The method exhibits satisfactory LoD and LoQ values, excellent linearity inside the operative range of concentration (i.e., between 2 and 200 mg dm −3 ). The method was tested on ten real honey samples that were different in age and floral origin. An excellent accuracy has been demonstrated for blossom as well as eucalyptus and thistle honeys, while for strawberry tree samples, homogentisic acid, the chemical marker of this floral origin, also contributed to FI as it was quantitatively reduced in the analytical conditions. However, RP-HPLC determination of HA allows evaluating the freshness of these samples by FI. So, it is possible to believe that the proposed method could be advantageously used for screening purposes due its simplicity and its speed.

Institutional Review Board Statement: Not applicable.
Informed Consent Statement: Not applicable.

Data Availability Statement:
The data used to support the findings of this study are available from the corresponding author upon request.