HPLC Analysis of Phenolic Compounds and Flavonoids with Overlapping Peaks

Abstract

The identification and quantification of phenolic compounds and flavonoids in various natural food products is typically conducted using HPLC analysis. Their analysis is particularly complex since most natural food products contain a large number of different phenolic compounds, many of which have similar chemical characteristics such as polarity, which makes complete separation of all eluents extremely difficult. In this work we present and validate a method for the quantitative determination of the concentration of two compounds with similar retention times, i.e. they show overlapping peaks in a mixed solution. Two pairs of phenolic compounds were investigated: caffeic and vanillic acids and ferulic and p-coumaric acids. This technique takes advantage of the different absorbances of the two phenolic compounds in the eluent at various wavelengths and can be used for the quantitative determination of the concentration of these compounds even if they are not separated in the HPLC column. The presented method could be used to interpret the results of HPLC analysis of food products which possess a vast spectrum of phenolic compounds and flavonoids.

Headnote

SUMMARY

The identification and quantification of phenolic compounds and flavonoids in various natural food products is typically conducted using HPLC analysis. Their analysis is particularly complex since most natural food products contain a large number of different phenolic compounds, many of which have similar chemical characteristics such as polarity, which makes complete separation of all eluents extremely difficult. In this work we present and validate a method for the quantitative determination of the concentration of two compounds with similar retention times, i.e. they show overlapping peaks in a mixed solution. Two pairs of phenolic compounds were investigated: caffeic and vanillic acids and ferulic and p-coumaric acids. This technique takes advantage of the different absorbances of the two phenolic compounds in the eluent at various wavelengths and can be used for the quantitative determination of the concentration of these compounds even if they are not separated in the HPLC column. The presented method could be used to interpret the results of HPLC analysis of food products which possess a vast spectrum of phenolic compounds and flavonoids.

Key words: HPLC analysis, UV-Vis absorbance, quantification of phenolic compounds and flavonoids, overlapping peaks

(ProQuest: ... denotes formulae omitted.)

INTRODUCTION

Phenolic compounds and flavonoids are a class of natural products commonly found in food products (1-4) including herbs (5), wine (6-8), beer (9), olive oil (10-13), fruits (14,15) and honey (16-20). Despite being present in relatively small concentrations, these compounds are known to impart beneficial properties to these food products such as antimicrobial, food preservation and antioxidant properties (8,17,19,21-27). The amount and type of these compounds depends primarily on the product type and location, and in the case of honey, floral sources, so they can also sometimes serve as chemical fingerprints to trace the geographic and botanical origins of the food products.

The identification and quantification of phenolic compounds and flavonoids in food products is typically conducted using HPLC analysis with a UV-Vis diode array detector (DAD) (7,25,28-37). The regular modus operandi involves the isolation and extraction of phenolic compounds from the food product, followed by an HPLC run using a gradient mobile phase consisting of two or more reagents, which are typically a polar organic solvent such as methanol or acetonitrile and a weak acid such as phosphoric or acetic acid (30,31,38). The analytes are then identified and quantified by comparison against standard solutions. While this method is perfectly valid and accurate for certain food products, it may however prove to be insufficient for the analysis of products such as olive oil, wine and honey, which contain a considerably large assortment of natural products, most of which are chemically related and have similar polarity. This can make separation of peaks problematic, resulting in some cases in amalgamated peaks, which makes it difficult to determine the exact concentration of certain compounds, or indeed, in some situations, even to simply ascertain their presence in food products, particularly if most of the peaks in the spectrum are unidentified. In such scenarios, it is extremely unlikely that an analysis based solely on a single HPLC spectrum is

sufficient to obtain a completely accurate and reliable characterization and quantification of these compounds.

In view of this, the objective of this work is to propose a method that can be used to identify and quantify with a high degree of certainty fifteen phenolic compounds commonly found in a variety of natural food products ranging from honey and olive oil to fruit juices. The specific aim is the determination of the concentration of phenolic compounds that have overlapping peaks by taking advantage of their diverse absorbances at different wavelengths. Accurate determination of the individual concentrations of phenolic compounds having peaks with identical retention times in a mixture is the ultimate objective.

MATERIALS AND METHODS

Preparation of standards and phenolic mixtures

Standard solutions were prepared for the fifteen investigated phenolic compounds and flavonoids, namely: kaempferol, luteolin, phenylacetic acid (all Alfa Aesar, Ward Hill, MA, USA), apigenin, chrysin, quercetin, p-coumaric acid, naringenin (all Sigma-Aldrich, Merck, Darmstadt, Germany), ferulic, syringic, vanillic, caffeic, ellagic, gallic and benzoic acids (Acros Organics, Geel, Belgium). The solutions were prepared by dissolving the standards in HPLC grade methanol (ultragradient grade; Carlo Erba, Milan, Italy) to produce stock solutions of 100 mg/L, which were then used to prepare 50, 40, 30, 20 and 10 mg/L solutions for the standard plots. In addition, a mixture containing 30 mg/L of each phenolic compound in methanol was also prepared. Two mixtures of p-coumaric and ferulic acid were prepared, one with equal concentrations of 50 mg/L each and the other with 30 mg/L p-coumaric and 70 mg/L ferulic acid. Another similar set of mixtures was also prepared using vanillic and caffeic acids. We measured the absorbance of the samples of 100 mg/L of each phenolic compound and flavonoid with a UV-Vis spectrophotometer (UV-2600; Shimadzu, Tokyo, Japan) at wavelengths between 180 and 480 nm to find the optimum wavelength for the HPLC-DAD measurements.

HPLC analysis

The HPLC analysis of the phenolic compounds and flavonoids was conducted using a Waters 2695 Alliance HPLC system (Waters Inc., Milford, CT, USA), equipped with a UV-Vis DAD. The separation was conducted using a Waters Sunfire™ C18 reverse-phase chromatography column, 250 mm length, 4.6 mm width, and particle size 5 pm. The phenolic standard solutions and mixtures were injected into the system using an autoinjector. Different isocratic and gradient mobile phases were tested at different flow rates and column temperatures in order to find a suitable separation method for the standards.

The gradient method that was eventually chosen following a series of preliminary studies uses a mixture of acetonitrile (mobile phase A, HPLC grade >99.9 %; Honeywell Seelze, Germany) and phosphoric acid (mobile phase B), which was prepared by dropwise addition of 85 % orthophosphoric acid (Sigma-Aldrich, Merck, Darmstadt, Germany) to HPLC grade water (Carlo Erba) until pH=2 was reached. The total runtime of the method was 60 min and the concentration gradient was varied as follows: a) initially 5 % A and 95 % B, b) 15 min 35 % A and 65 % B, c) 20 min 35 % A and 65 % B, d) 30 min 40 % A and 60 % B, e) 35 min 40 % A and 60 % B, f) 40 min 50 % A and 50 % B, g) 52 min 70 % A and 30 % B and h) 60 min 5 % A and 95 % B. A constant flow rate of 0.5 mL/min and a temperature of 5 °C were used. Following the analysis of the UV-Vis spectra of the individual phenolic standards, three wavelengths (210, 280 and 360 nm) were chosen for analysis in this investigation using the HPLC-DAD.

RESULTS AND DISCUSSION

Fig. 1 shows the chromatograms of the solution containing all 15 phenolic compounds obtained at wavelengths of 210, 280 and 360 nm. Table 1 and Table 2 show the retention times and calibration constants based on the area and height of peaks for each phenolic compound for every wavelength, respectively.

As one can observe from the chromatograms in Fig. 1, the gradient method used here separates most phenolic compounds reasonably well with most of them showing distinct and sharp individual peaks. Moreover, while all phenolic compounds show peaks at 210 and 280 nm (except for phenylacetic acid at 280 nm), luteolin, ferulic acid, caffeic acid, p-coumaric acid, ellagic acid, apigenin, kaempferol and quercetin also show peaks at 360 nm. These results are in accordance with those obtained from the initial tests conducted using a UV-Vis spectrophotometer to determine the choice of wavelengths.

However, as one may observe in Fig. 1, there are also two pairs of phenolic compounds which have identical retention times and, hence overlapping peaks: vanillic and caffeic acids at 21.0 min and ferulic and p-coumaric acids at 24.5 min. Yet this drawback does not necessarily mean that it is impossible to determine the individual concentrations of these phenolic compounds. As evident from the values in Table 1 and Table 2, each phenolic compound has a different absorption profile. It is possible to take advantage of this property to determine the concentration of each phenolic compound in the mixture by using the standardization constants of the individual phenolic compounds and the total absorbance of the phenolic mixture at different wavelengths.

The method proposed here operates under the assumption that the total area of the peak at a given wavelength, Ærotai, is equal to the sum of the individual areas of the phenolic compounds, /4ph, making up the peak, Ph¡ and Phj, at the same wavelength, A¡ :

... /1/

This relationship is valid for all wavelengths and thus Eq. 1 can be used to generate the following simultaneous equations for the peaks obtained at two different wavelengths:

.. /2/

and

/3/

These equations can be expanded to incorporate the terms defining the concentrations (^Ph) of the phenolic compounds and the standardization gradient (k^ ) and y-intercept constants (bļ ), which are related to the area (Api ), through the following equation:

... /4/

thus:

... /5/

and

.. /6/

As one may observe from Eqs. 5 and 6, the terms ľ Ph1 and ľ Ph2 are common for both equations and thus, since all the other terms are known, one may obtain the values for these concentrations by solving the two simultaneous equations. The final values for ľPh1 and ľPh2 may be expressed as follows:

... /7/

and

... /8/

These equations may be used to calculate the concentrations of p-coumaric and ferulic acids since these two phenolics have very similar retention times and absorb to different extents at all of the three wavelengths used here. In the case of vanillic and caffeic acids, the problem is simpler since while the latter absorbs at all three wavelengths, the former absorbs only at ļ=210 and 280 nm. Therefore, Eqs. 7 and 8 may be simplified as follows to calculate the concentrations of these phenolics when considering a wavelength of 210 or 280 nm in conjunction with the 360 nm wavelength, since in the latter case Ai^is equal to Ah2:

... /9/

and

... /10/

where Ph1 is the phenolic compound that absorbs at both evaluated wavelengths, in this case caffeic acid, Ph2 is the other, i.e. vanillic acid, A1 is the wavelength at which both phenolics absorb, in this case 210 or 280 nm, and X2 is the wavelength at which only one phenolic compound absorbs, in this case 360 nm.

In order to validate the effectiveness of this method, Eqs. 7-10 were applied to determine the concentrations of two mixtures of vanillic and caffeic acids (mixtures 1 and 2) and ferulic and p-coumaric acid (mixtures 3 and 4) with known concentrations. The concentrations of these mixtures were calculated using the peak areas from three data sets of wavelengths: 210-280, 210-360 and 280-360 nm and the results are presented in Table 3. A comparison between the real and the calculated concentrations of the mixtures is also shown in Fig. 2.

It is evident from the data in Table 3 and the plot in Fig. 2 that the values obtained through the equations are extremely similar to the actual concentrations of the individual phenolic acids making up each of the four mixed solutions. In fact, in the case of the mixtures of p-coumaric and ferulic acids, the calculated values were all within ±0.5 mg/L of the actual values, indicating a high degree of accuracy. On the other hand, in the cases of vanillic and caffeic acid mixtures, there are slightly bigger discrepancies between the points, although overall the average predictions of each combination of wavelengths are still very close to the actual values.

These results confirm the validity of Eqs. 7-10 for calculating the concentrations of HPLC analytes with overlapping peaks based on their varying absorbances at different wavelengths. In theory, such a technique should also be applicable to peak height and peak area; however, this is only the case if the phenolic compounds in question possess exactly the same retention times. In the cases presented here the two pairs of phenolic compounds have extremely similar but not exact retention times. This means that while a single large peak is obtained for the mixture, it is wider as well as higher than the individual peaks and thus while the cumulative peak areas of the individual phenolic compounds conform to the assumption presented in Eq. 1, the same cannot be said for the cumulative peak heights:

... /11/

where H represents the peak height. In fact, this is evident from the results presented in Table 4 and Fig. 3, where calculations corresponding to Eqs. 7-10 but based on peak height are presented. As one may observe, the calculated values obtained with this method consistently underestimate the phenolic concentration by a large extent, hence confirming the inadmissibility of this method when applied to peak height data.

At this point it is important to mention the advantages of using the method described here to analyze UV-Vis HPLC-DAD results. Although typically it is important to ensure that all the analytes separate completely, this is not always so easily achieved, particularly in the case of natural food products such as honey (also evident from previous works (39-41)), which are known to contain over fifty different types of phenolic compounds and flavonoids. In such cases, finding a gradient method which is capable of achieving complete separation of all constituents is almost impossible, especially since many of these phenolic compounds have extremely similar chemical composition and polarities. By using the method presented here one may possibly circumvent this problem, particularly if like in the case described here, the gradient method is capable of completely separating the majority of phenolic compounds, and therefore there is no need to develop another method solely to separate a couple of peaks. Moreover, the equations described in this methodology can also be used to conduct a qualitative analysis in order to determine if any unknown compounds have overlapping peaks with the target compounds under analysis. If using the equations to calculate the concentrations of two phenolic compounds over multiple pairs of wavelength combinations results in different calculated values, then this is indicative of the presence of possibly a third, unknown eluent contributing to the peak area. On the other hand, if all combinations of wavelengths return the same concentrations, then this confirms that only the two phenolics in question are present at this retention time. Currently, the standard method used to counteract this problem is to either use multiple UV-Vis absorpti6n-based HPLC protocols with different gradient methods and/or mobile phases such as that employed by Gupta et al. (35), or else to validate the initial HPLC results using additional detectors such as a mass spectrometer (39,42-44). The m20hod proposed in this work eliminates the need of using such t20hniques as a lalidation meth60 for a UV-Vis absorption-based HPLC analysi s. This would facilitate the analysis of complex solutions since all the results required for this analysis may be obtained from a single HPLC run. However, it should be emphasized that the technique proposed here would replace these techniques for validation and quantification purposes only, and that the use of additional methods such as MS-HPLC is still required for the eventual characterization and identification of any unknown compounds in natural products. Furthermore, this technique could also be potentially employed as a quality control method for the analysis of synthetic products containing phenolic compounds and flavonoids. In such cases where the constituents are already known, a partial HPLC separation coupled with the method applied here could be sufficient to quantify the individual phenolic compound content.

CONCLUSIONS

In this work, we presented and validated an HPLC analysis method that can be used to find the concentrations of eluents with similar retention times in a mixture. The analysis was conducted on a mixture of fifteen phenolic compounds, with two pairs of phenolic compounds having peaks with nearly identical retention times, using UV-Vis absorbance measurements from an HPLC-DAD. The results obtained from the equations used to calculate the concentrations based on the peak area standardization constants of the individual phenolic compounds showed excellent agreement with the known concentrations of the mixtures and indicated that this technique could be a viable method to quantitatively analyze the concentrations of such eluents. It is envisaged that this technique could be applied for HPLC analysis of food products such as olive oil, fruit juices and honey, which have a vast spectrum of phenolic compounds and flavonoids with similar chemical characteristics and thus yield complex chromatograms that are extremely difficult to interpret accurately.

ACKNOWLEDGEMENTS AND FUNDING

This work was partly supported by the COST CA15118 Mathematical and Computer Science Methods for Food Science and Industry (FoodMC).

CONFLICT OF INTEREST

The authors declare no conflict of interest.

ORCID ID

L. Mizzi © https://ordd.org/0000-0002-7650-1173 C. Chatzitzika © https://ordd.org/0000-0002-2701-108X R. Gatt https://ordd.org/0000-0002-1950-743X V. Valdramidis https://ordd.org/0000-0001-6486-3890

Sidebar

Received: 2 June 2019

Accepted: 17 February 2020

Footnote

* Corresponding author:

E-mail: [email protected], [email protected]

References

REFERENCES

1. Bravo L. Polyphenols: Chemistry, dietary sources, metabolism, and nutritional significance. Nutr Rev. 1998;56(11):317-33. https://doi.org/10.1111/j.1 753-4887.1998.tb01670.x

2. Proestos C, Boziaris IS, Nychas GJE, Komaitis M. Analysis of flavonoids and phenolic acids in Greek aromatic plants: Investigation of their antioxidant capacity and antimicrobial activity. Food Chem. 2006;95(4):664-71. https://doi.org/10.1016/j.foodchem.2005.01.049

3. Shahidi F, Yeo JD. Insoluble-bound phenolics in food. Molecules. 2016;21(9):1216. https://doi.org/10.3390/molecules21091216

4. Pérez-Gregorio RM, García-Falcón MS, Simal-Gándara J, Rodrigues AS, Almeida DPF. Identification and quantification of flavonoids in traditional cultivars of red and white onions at harvest. J Food Compos Anal. 2010;23(6):592-8. https://doi.org/10.1016/j.jfca.2009.08.013

5. Wojdyło A, Oszmiański J, Czemerys R. Antioxidant activity and phenolic compounds in 32 selected herbs. Food Chem. 2007;105(3):940-9. https://doi.org/10.1016/j.foodchem.2007.04.038

6. Lasekan O, Buettner A, Christlbauer M. Investigation of important odorants of palm wine (Elaeisguineensis). Food Chem. 2007;105(1):15-23. https://doi.org/10.1016/j.foodchem.2006.12.052

7. Mattila P, Kumpulainen J. Determination of free and total phenolic acids in plant-derived foods by HPLC with diode-array detection. J Agric Food Chem. 2002;50(13):3660-7. https://doi.org/10.1021/jf020028p

8. Rodríguez Vaquero MJ, Alberto MR, Manca de Nadra MC. Antibacterial effect of phenolic compounds from different wines. Food Control. 2007;18(2):93-101. https://doi.org/10.1016/j.foodcont.2005.08.010

9. García AA, Cancho Grande B, Simal Gándara J. Development of a rapid method based on solid-phase extraction and liquid chromatography with ultraviolet absorbance detection for the determination of polyphenols in alcohol-free beers. J Chromatogr A. 2004;1054(1-2):175-80. https://doi.org/10.1016/j.chroma.2004.07.092

10.Paradiso VM, Clemente A, Summo C, Pasqualone A, Caponio F. Towards green analysis of virgin olive oil phenolic compounds: Extraction by a natural deep eutectic solvent and direct spectrophotometric detection. Food Chem. 2016;212:43-7. https://doi.org/10.1016/j.foodchem.2016.05.082

11. Tuck KL, Hayball PJ. Major phenolic compounds in olive oil: Metabolism and health effects. J Nutr Biochem. 2002;13(11):636-44. https://doi.org/10.1016/S0955-2863(02)00229-2

12. Bianco A, Uccella N. Biophenolic components of olives. Food Res Int. 2000;33(6):475-85. https://doi.org/10.1016/S0963-9969(00)00072-7

13. Reboredo-Rodríguez P, Valli E, Bendini A, Di Lecce G, Simal-Gándara J, Toschi TG. A widely used spectrophotometric assay to quantify olive oil biophenols according to the health claim (EU Reg. 432/2012). Eur J Lipid Sci Technol. 2016;118(10):1593-9. https://doi.org/10.1002/ejlt.201500313

14. Abbasi AM, Guo X, Fu X, Zhou L, Chen Y, Zhu Y, et al. Comparative assessment of phenolic content and in vitro antioxidant capacity in the pulp and peel of mango cultivars. Int J Mol Sci. 2015;16(6):13507-27. https://doi.org/10.3390/ijms160613507

15. Häkkinen SH, Kärenlampi SO, Mykkänen HM, Heinonen IM, Törrönen AR. Ellagic acid content in berries: Influence of domestic processing and storage. Eur Food Res Technol. 2000;212:75-80. https://doi.org/10.1007/s002170000184

16. Pyrzynska K, Biesaga M. Analysis of phenolic acids and flavonoids in honey. Trends Anal Chem. 2009;28(7):893-902. https://doi.org/10.1016/j.trac.2009.03.015

17. Stephens JM, Schlothauer RC, Morris BD, Yang D, Fearnley L, Greenwood DR, Loomes KM. Phenolic compounds and methylglyoxal in some New Zealand manuka and kanuka honeys. Food Chem. 2010;120(1):78-86. https://doi.org/10.1016/j.foodchem.2009.09.074

18. Silici S, Sagdic O, Ekici L. Total phenolic content, antiradical, antioxidant and antimicrobial activities of Rhododendron honeys. Food Chem. 2010;121 (1):238-43. https://doi.org/10.1016/j.foodchem.2009.11.078

19. da Silva PM, Gauche C, Gonzaga LV, Costa ACO, Fett R. Honey: Chemical composition, stability and authenticity. Food Chem. 2016;196:309-23. https://doi.org/10.1016/j.foodchem.2015.09.051

20. Escriche I, Kadar M, Juan-Borrás M, Domenech E. Suitability of antioxidant capacity, flavonoids and phenolic acids for floral authentication of honey. Impact of industrial thermal treatment. Food Chem. 2014;142:135-43. https://doi.org/10.1016/j.foodchem.2013.07.033

21. Lou Z, Wang H, Rao S, Sun J, Ma C, Li J. p-Coumaric acid kills bacteria through dual damage mechanisms. Food Control. 2012;25(2):550-4. https://doi.org/10.1016/j.foodcont.2011.11.022

22. Alves MJ, Ferreira ICFR, Froufe HJC, Abreu RMV, Martins A, Pintado M. Antimicrobial activity of phenolic compounds identified in wild mushrooms, SAR analysis and docking studies. J Appl Microbiol. 2013;115(2):346-57. https://doi.org/10.1111/jam .12196

23. Alvarez-Suarez JM, Giampieri F, González-Paramás AM, Damiani E, Astolfi P, Martinez-Sanchez G, et al. Phenolics from monofloral honeys protect human erythrocyte membranes against oxidative damage. Food Chem Toxicol. 2012;50(5):1508-16. https://doi.org/10.1016/j.fct.2012.01.042

24. Challacombe CA, Abdel-Aal ESM, Seetharaman K, Duizer LM. Influence of phenolic acid content on sensory perception of bread and crackers made from red or white wheat. J Cereal Sci. 2012;56(2):181-8. https://doi.org/10.1016/j.jcs.2012.03.006

25. Aljadi AM, Yusoff KM. Isolation and identification of phenolic acids in Malaysian honey with antibacterial properties. Turk J Med Sci. 2003;33(4):229-36.

26. Estevinho L, Pereira AP, Moreira L, Dias LG, Pereira E. Antioxidant and antimicrobial effects of phenolic compounds extracts of Northeast Portugal honey. Food Chem Toxicol. 2008;46(12):3774-9. https://doi.org/10.1016/j.fct.2008.09.062

27. Wahdan HAL. Causes of the antimicrobial activity of honey. Infection. 1998;26(1):26-31. https://doi.org/10.1007/BF02768748

28. Peñarrieta JM, Alvarado JA, Åkesson B, Bergenståhl B. Separation of phenolic compounds from foods by reversed-phase high performance liquid chromatography. Rev Boliv Química. 2007;24(1):1-4.

29. Hadjmohammadi MR, Nazari S, Kamel K. Determination of flavonoid markers in honey with SPE and LC using experimental design. Chroma. 2009;69(11 -12):1291 -7. https://doi.org/10.1365/s10337-009-1073-4

30. Khalil MI, Alam N, Moniruzzaman M, Sulaiman SA, Gan SH. Phenolic acid composition and antioxidant properties of Malaysian honeys. J Food Sci. 2011;76(6):C921-8. https://doi.org/10.1111/j.1750-3841.2011.02282.x

31. Moniruzzaman M, Yung An C, Rao PV, Hawlader MNI, Azlan SABM, Sulaiman SA, Gan SH. Identification of phenolic acids and flavonoids in monofloral honey from Bangladesh by high performance liquid chromatography: Determination of antioxidant capacity. Biomed Res Int. 2014;2014:Article ID 737490. https://doi.org/10.1155/2014/737490

32. Moniruzzaman M, Khalil MI, Sulaiman SA, Gan SH. Physicochemical and antioxidant properties of Malaysian honeys produced by Apis cerana, Apis dorsata and Apis mellifera. BMC Complement Altern Med. 2013;13:43. https://doi.org/10.1186/1472-6882-13-43

33. Ferreres F, Tomás-Barberán FA, Gil MI, Tomás-Lorente F. An HPLC technique for flavonoid analysis in honey. J Sci Food Agric. 1991;56(1):49-56. https://doi.org/10.1002/jsfa.2740560106

34. Andrade P, Ferreres F, Amaral MT. Analysis of honey phenolic acids by HPLC, its application to honey botanical characterization. J Liq Chromatogr Relat Technol. 1997;20(14):2281 -8. https://doi.org/10.1080/10826079708006563

35. Gupta M, Sasmal S, Majumdar S, Mukherjee A. HPLC profiles of standard phenolic compounds present in medicinal plants. IJPPR. 2012;4(3):162-7.

36. Nour V, Trandafir I, Cosmulescu S. HPLC determination of phenolic acids, flavonoids and juglone in walnut leaves. J Chromatogr Sci. 2013;51(9):883-90. https://doi.org/10.1093/chromsci/bms180

37. Tumbas VT, Mandić AI, Ćetković GS, D strok signilas SM, Čanadanović-Brunet JM. HPLC analysis of phenolic acids in mountain germander (Teucrium montanum L.) extracts. Acta Period Technol. 2004;35:265-73. https://doi.org/10.2298/APT0435265T

38. Michalkiewicz A, Biesaga M, Pyrzynska K. Solid-phase extraction procedure for determination of phenolic acids and some flavonols in honey. J Chromatogr A. 2008;1187(1 -2):18-24. https://doi.org/10.1016/j.chroma.2008.02.001

39. Kečkeš S, Gašić U, Ćirković Veličković T, Milojković-Opsenica D, Natić M, Tešić Ž. The determination of phenolic profiles of Serbian unifloral honeys using ultra-high-performance liquid chromatography/high resolution accurate mass spectrometry. Food Chem. 2013;138(1):32-40. https://doi.org/10.1016/j.foodchem.2012.10.025

40. Tomás-Barberán FA, Martos I, Ferreres F, Radovic BS, Anklam E. HPLC flavonoid profiles as markers for the botanical origin of European unifloral honeys. J Sci FoodAgric. 2001;81(5):485-96. https://doi.org/10.1002/jsfa.836

41. Yao L, Datta N, Tomás-Barberán FA, Ferreres F, Martos I, Singanusong R. Flavonoids, phenolic acids and abscisic acid in Australian and New Zealand Leptospermum honeys. Food Chem.2003;81 (2):159-68. https://doi.org/10.1016/S0308-8146(02)00388-6

42. Pulcini P, Allegrini F, Festuccia N. Fast SPE extraction and LC-ESI-MS-MS analysis of flavonoids and phenolic acids in honey. Apiacta. 2006;41:21-7.

43. Trautvetter S, Koelling-Speer I, Speer K. Confirmation of phenolic acids and flavonoids in honeys by UPLC-MS. Apidologie. 2009;40(2):140-50. https://doi.org/10.1051/apido/2008072

44. Figueiredo-González M, Regueiro J, Cancho-Grande B, Simal-Gándara J. Garnacha Tintorera-based sweet wines: Detailed phenolic composition by HPLC/DAD-ESI/MS analysis. Food Chem. 2014;143;282-92. https://doi.org/10.1016/j.foodchem.2013.07.120