Characteristics of Indonesian Stingless Bee Propolis and Study of Metabolomic Properties Based on Region and Species

The chemical compounds found in propolis vary according to plant sources, species, and geographical regions. To date, Indonesian propolis has not yet become standardized in terms of its chemical constituents. Thus, this study aimed to identify the presence of marker compounds and determine whether different classes of Indonesian propolis exist. In this study, yields, total polyphenol content (TPC), total flavonoid content (TFC), and antioxidants were measured. Identification of chemical compounds was carried out with Fourier-transform infrared (FTIR) spectroscopy and liquid chromatography-tandem mass spectrometry (LC-MS/MS). Metaboanalyst 6.0 was employed in conducting principal component analysis (PCA) and partial least squares discriminant analysis (PLS-DA) using the results of the FTIR and LC-MS/MS. The propolis with the highest TFC, TPC, and antioxidant activity was Geniotrigona thoracica from North Sumatra. The results of propolis compound mapping based on region with discriminant analysis revealed that types of propolis from Java have similar characteristics. Then, based on species, the types of propolis from Tetragonula laeviceps and Heterotrigona itama have special characteristics; the samples from these species can be grouped according to similar characteristics. In conclusion, 10 potential marker compounds were identified in Indonesian propolis, enabling regional and species-specific varieties of Indonesian propolis to be classified based on chemical composition mapping.


Introduction
Propolis, a natural resin produced by bees, is used in traditional medicine and in food, beverages, and supplements to improve health and prevent infection, inflammation, heart disease, and cancer [1,2].Propolis also has a range of applications across various industries due to its unique properties.Here are some notable uses of propolis in different sectors: cosmetics and personal care, pharmaceuticals, agriculture, and veterinary medicine [3][4][5][6].
The array of chemical compounds contained in propolis depends on plant sources, geographical location, environmental conditions, and the species of bee [7][8][9][10].These factors result in a wide diversity in the chemical compounds found in propolis across regions, the most common being polyphenols and flavonoids [4][5][6].Indonesia's high biodiversity is due to its position between two continental shelves, with the Wallace Line separating Asian and Australian faunal regions through the Lombok Strait between Bali and Lombok, and between Borneo and Sulawesi [11].
The metabolomic study investigated the metabolites present in propolis with a comprehensive analysis.The study identified and categorized a diverse range of propolis metabolites [12,13].Secondly, the study investigated the potential health benefits of the identified metabolites, such as antibacterial [14][15][16], antioxidant [15,[17][18][19], anti-inflammatory [20,21], and anti-cancer activities [22,23].Last, it examined how the metabolomic profile of propolis varies with its region [24][25][26] and botanical sources [27].The propolis extract in metabolomic research was identified using different analytical techniques, such as liquid chromatography-mass spectrometry (LC-MS) and gas chromatography-mass spectrometry (GC-MS), to identify and measure the amounts of metabolites [24,26,[28][29][30][31]. Nuclear Magnetic Resonance (NMR) spectroscopy is also used to determine the structure of metabolites [12,25,[32][33][34][35]. Fourier-transform infrared (FTIR) can be used to evaluate the functional group of the propolis compound [17].The data analysis of the metabolomic study employs bioinformatics techniques to scrutinize the data, discern metabolites, and elucidate the findings.The bioactivity assays involve the performance of in vitro and in vivo studies to evaluate the biological activities of the discovered metabolites.
According to the current literature, a key subject of interest is the profile of the chemical compounds in propolis.So, the point of this study was to find out if there are different types of Indonesian propolis based on region and species, as well as to see if there are different classes of Indonesian propolis.The establishment of stingless bee propolis will facilitate testing the quality, support marketing, and aid in increasing the production of Indonesian propolis.

Bee Species
The bee species was identified to confirm the origin of propolis.The morphology of stingless bees H. itama, G. thoracica, T. laeviceps, and T. biroi is presented in Figures 1-4.
Molecules 2024, 29, x FOR PEER REVIEW 2 of 22 medicine [3][4][5][6].The array of chemical compounds contained in propolis depends on plant sources, geographical location, environmental conditions, and the species of bee [7][8][9][10].These factors result in a wide diversity in the chemical compounds found in propolis across regions, the most common being polyphenols and flavonoids [4][5][6].Indonesia's high biodiversity is due to its position between two continental shelves, with the Wallace Line separating Asian and Australian faunal regions through the Lombok Strait between Bali and Lombok, and between Borneo and Sulawesi [11].
The metabolomic study investigated the metabolites present in propolis with a comprehensive analysis.The study identified and categorized a diverse range of propolis metabolites [12,13].Secondly, the study investigated the potential health benefits of the identified metabolites, such as antibacterial [14][15][16], antioxidant [15,[17][18][19], antiinflammatory [20,21], and anti-cancer activities [22,23].Last, it examined how the metabolomic profile of propolis varies with its region [24][25][26] and botanical sources [27].The propolis extract in metabolomic research was identified using different analytical techniques, such as liquid chromatography-mass spectrometry (LC-MS) and gas chromatography-mass spectrometry (GC-MS), to identify and measure the amounts of metabolites [24,26,[28][29][30][31]. Nuclear Magnetic Resonance (NMR) spectroscopy is also used to determine the structure of metabolites [12,25,[32][33][34][35]. Fourier-transform infrared (FTIR) can be used to evaluate the functional group of the propolis compound [17].The data analysis of the metabolomic study employs bioinformatics techniques to scrutinize the data, discern metabolites, and elucidate the findings.The bioactivity assays involve the performance of in vitro and in vivo studies to evaluate the biological activities of the discovered metabolites.
According to the current literature, a key subject of interest is the profile of the chemical compounds in propolis.So, the point of this study was to find out if there are different types of Indonesian propolis based on region and species, as well as to see if there are different classes of Indonesian propolis.The establishment of stingless bee propolis will facilitate testing the quality, support marketing, and aid in increasing the production of Indonesian propolis.

Bee Species
The bee species was identified to confirm the origin of propolis.The morphology of stingless bees H. itama, G. thoracica, T. laeviceps, and T. biroi is presented in Figures 1-4.

Total Phenolic Content (TPC) Results
The calculation of the gallic acid calibration curve resulted in the following equation: y = 0.0074x + 0.0376 (R 2 = 0.999).After the sample absorbance data were obtained, the original concentration of the sample in ppm could be calculated using the calibration curve equation.Subsequently, the conversion was carried out for the TPC of extracts of stingless bee propolis (EEP) from various Indonesian regions (Table 1).The TPC ranged from 32.99 ± 0.66 to 374.20 ± 0.51 mg of gallic acid equivalents (GAE)/g.Propolis from G. thoracica from Pematang Siantar, Sumatra (sample code Sum_GT4), had the highest concentration of TPC, whereas the lowest total value of polyphenol content was found in propolis from Barito Kuala, Kalimantan (sample code Kal_TL5).

Total Phenolic Content (TPC) Results
The calculation of the gallic acid calibration curve resulted in the following equation: y = 0.0074x + 0.0376 (R 2 = 0.999).After the sample absorbance data were obtained, the original concentration of the sample in ppm could be calculated using the calibration curve equation.Subsequently, the conversion was carried out for the TPC of extracts of stingless bee propolis (EEP) from various Indonesian regions (Table 1).The TPC ranged from 32.99 ± 0.66 to 374.20 ± 0.51 mg of gallic acid equivalents (GAE)/g.Propolis from G. thoracica from Pematang Siantar, Sumatra (sample code Sum_GT4), had the highest concentration of TPC, whereas the lowest total value of polyphenol content was found in propolis from Barito Kuala, Kalimantan (sample code Kal_TL5).

Total Phenolic Content (TPC) Results
The calculation of the gallic acid calibration curve resulted in the following equation: y = 0.0074x + 0.0376 (R 2 = 0.999).After the sample absorbance data were obtained, the original concentration of the sample in ppm could be calculated using the calibration curve equation.Subsequently, the conversion was carried out for the TPC of extracts of stingless bee propolis (EEP) from various Indonesian regions (Table 1).The TPC ranged from 32.99 ± 0.66 to 374.20 ± 0.51 mg of gallic acid equivalents (GAE)/g.Propolis from G. thoracica from Pematang Siantar, Sumatra (sample code Sum_GT4), had the highest concentration of TPC, whereas the lowest total value of polyphenol content was found in propolis from Barito Kuala, Kalimantan (sample code Kal_TL5).

Total Flavonoid Content (TFC) Results
The calculation of the quercetin calibration curve resulted in the following equation: y = 0.0069x − 0.0687 (R 2 = 0.9974).After the absorbance data were obtained from the sample, the original concentration of the sample in ppm could be calculated using the calibration curve equation.According to the conversion that was subsequently carried out, the TFC values of EEP for various islands are shown in Table 1.The TFC of Indonesian stingless bee propolis ranged from 20.66 ± 4.08 to 194.18 ± 0.98 mg QE/g.Propolis from Pematang Siantar, Sumatra (sample code Sum_GT4), had the highest TFC.Meanwhile, the lowest overall value of flavonoid content was found in propolis from Barito Kuala, Kalimantan (sample code Kal_TL5).

Antioxidant Activity Results
The results for the antioxidant activity of the EEP using the DPPH method are shown in Table 2, showing the lowest IC 50 for T. laeviceps propolis from Kendal, Central Java, at 118.1 ± 2362 ppm (moderate).The highest IC 50 activity was 11.06 ± 0.55 ppm (very strong) in G. thoracica propolis samples from Pematang Siantar, North Sumatra.The antioxidant activity of EEP determined using the ABTS method showed that IC 50 ranged from 17.29 ± 0.34 to 117.21 ± 10.76 ppm.The EEP G. thoracica from Pematang Siantar, North Sumatra, had the strongest activity.The samples with lowest activity were propolis samples T. biroi from Batang, Central Java.The classification of antioxidant strength in Table 2 is based on previous research, which classifies antioxidant strength as very strong (<50 ppm), strong (50-100 ppm), moderate (101-250 ppm), weak (250-500 ppm), and inactive (>500 ppm) [36].

Chemometric Analysis of FTIR Spectroscopies
The results of the chemometric analysis of the principal component analysis (PCA), partial least squares discriminant analysis (PLS-DA), and sparse partial least squares discriminant analysis (sPLS-DA) of 23 samples of propolis using Metaboanalyst 6.0 are shown in Figure 5.The score plots of PCA, PLS-DA, and SPLS-DA show consistent clustering and discriminant results.From the score plot of chemometrics analysis, the red area of Group 1 consists of all samples from the same species, H. itama.Its regions were Kendal (Java), Magelang (Java), Banjar Baru (Kalimantan), Kutai Kartanegara (Kalimantan), and Kampar (Sumatra).Group 2 with the purple color is G. thoracica from Kutai Kartanegara (Kalimantan) and Lebong (Sumatra).The blue area, Group 3, with its blue color, consists of eleven samples from species T. biroi, T. laeviceps, T. drescheri, and G. thoracica.This group included all T. laeviceps propolis samples.The region of samples in this group was Tanah Laut (Kalimantan), Batang (Java), Pandeglang (Java), Magelang (Java), Hulu Sungai Tengah (Kalimantan), and Barito Kuala (Kalimantan).The T. biroi propolis in Group 3 was from Tanah Laut (Kalimantan), Batang (Java), Bogor (Java), and Magelang (Java).The last species in this group comprised propolis (T.drescheri) from Batang (Java) and G. thoracica from Hulu Sungai Tengah (Kalimantan).The green area (Group 4) consists of T. biroi, T. clypearis (Nusa Tenggara), and T. biroi (Kalimantan).The last area, with a pink color (Group 5), consists of T. clypearis (Nusa Tenggara) and G. thoracica (Pematang Siantar).Table 3 describes each member in every group.The propolis spectra from FTIR per group are shown in Figures 6-10.They demonstrate that stingless bee propolis from the various species and regions represented in the study revealed a spectrum with different vibrations, characterized by transmittance in wavenumber.We can attribute this observation to the similar composition of certain species and inter-island propolis compounds.Table 4 displays the combined readings of the propolis vibrational bands and compounds identified from the FTIR results.FTIR analysis can be used in metabolomics studies and is considered an efficient and economical analysis technique for classifying bioactive metabolic profiles, but it has the limitation of not being able to identify marker molecules from the sample.We conducted LC-MS/MS analysis of propolis samples and considered the most suitable analytical technique for characterizing propolis bioactive metabolic profiles from those species, especially H. itama and T. laeviceps [37].

LC-MS/MS Analysis Result
The study analyzed the bioactive compounds of 11 samples of propolis (all H. itama and T. laeviceps) using LC-MS.The spectra of propolis G. thoracica from Pematang Siantar (sample code Sum_GT4), which has the highest phenolic and flavonoid content and also the best activity as an antioxidant, can be seen in Figure 11 and Table 5.

Chemometric Analysis of LC-MS/MS
This study employed LC-MS/MS to identify the chemical compounds of propolis samples.The chemometric analysis of the FTIR sample shown in Figure 5 shows that H. itama from Java (Kendal, Magelang), Kalimantan (Banjar Baru, Kutai Kartanegara), and Sumatra (Kampar) have a near score plot, so they cluster in Group 1. Propolis samples from T. laeviceps from Java (Pandeglang, Batang, and Magelang) and Kalimantan (Hulu Sungai Tengah and Barito Kuala) also have a near score plot, so they cluster in Group 3.
The chemometric analysis with MetaboAnalyst 6.0 software showed good discriminant analysis and metabolomic compounds.The results of OPLS-DA are shown in Figure 12.The identification of similar compounds that occurred in those species groups was carried out as a means of searching for marker compounds.We identified 15 distinct compounds in Indonesian propolis overall.This observation suggests that these

Chemometric Analysis of LC-MS/MS
This study employed LC-MS/MS to identify the chemical compounds of propolis samples.The chemometric analysis of the FTIR sample shown in Figure 5 shows that H. itama from Java (Kendal, Magelang), Kalimantan (Banjar Baru, Kutai Kartanegara), and Sumatra (Kampar) have a near score plot, so they cluster in Group 1. Propolis samples from T. laeviceps from Java (Pandeglang, Batang, and Magelang) and Kalimantan (Hulu Sungai Tengah and Barito Kuala) also have a near score plot, so they cluster in Group 3.
The chemometric analysis with MetaboAnalyst 6.0 software showed good discriminant analysis and metabolomic compounds.The results of OPLS-DA are shown in Figure 12.The identification of similar compounds that occurred in those species groups was carried out as a means of searching for marker compounds.We identified 15 distinct compounds in Indonesian propolis overall.This observation suggests that these compounds have the potential to act as marker compounds.Table 6 presents a list of these compounds.compounds have the potential to act as marker compounds.Table 6 presents a list of these compounds.

TPC and TFC
The values of TPC and TFC of propolis samples from G. thoracica, H. itama, T. biroi, T. laeviceps, and another stingless bee were varied.According to the study results, the TPC and TFC of this study were lower than previous research on bees of the genus Tetragonula spp.from Luwu, South Sulawesi, with a range of 269.57± 20.37-426.91 ± 61.08 mg GAE/g and 211.08 ± 10.14 to 324.43 ± 11.84 mg QE/g [38].However, the results of this study are in accordance with the TPC results of propolis samples from Anatolian Turkey, with varying values of 16.37-125.83mg GAE/g [39].This study range is equivalent to Polish propolis, with a TFC between 35.64 and 62.04 QE/g [40].In terms of TPC value, this study is slightly higher than the reported samples of Malaysian propolis species H. fimbriata, T. apicalis, and T. binghami, with TPC values of 13.21 ± 0.26, 7.60 ± 0.13, and 10.11 ± 0.19 mg RE/g [41].In

TPC and TFC
The values of TPC and TFC of propolis samples from G. thoracica, H. itama, T. biroi, T. laeviceps, and another stingless bee were varied.According to the study results, the TPC and TFC of this study were lower than previous research on bees of the genus Tetragonula spp.from Luwu, South Sulawesi, with a range of 269.57± 20.37-426.91 ± 61.08 mg GAE/g and 211.08 ± 10.14 to 324.43 ± 11.84 mg QE/g [38].However, the results of this study are in accordance with the TPC results of propolis samples from Anatolian Turkey, with varying values of 16.37-125.83mg GAE/g [39].This study range is equivalent to Polish propolis, with a TFC between 35.64 and 62.04 QE/g [40].In terms of TPC value, this study is slightly higher than the reported samples of Malaysian propolis species H. fimbriata, T. apicalis, and T. binghami, with TPC values of 13.21 ± 0.26, 7.60 ± 0.13, and 10.11 ± 0.19 mg RE/g [41].In this study, the TFC value obtained was higher than the results reported by Fikri et al. in 2019, using Banten, South Sulawesi, and South Kalimantan propolis; the flavonoid range was obtained as 0.76 ± 0.26-3.39± 1.08 mg QE/g [42].
Variations in TPC values in this study show that different results can occur in the same bee species, and differences in the region of origin of propolis also show variations in TPC content.The TPC and TFC vary greatly, which can occur due to the type of vegetation around the colony, bee species, origin, solvent, and/or extraction method used [42].This difference in variation can be used as a standardization reference for the TPC and TFC content of stingless bee propolis in Indonesia, because it has been taken from several regional points and samples of different species.The challenge is to obtain globally agreedupon standards for use in the TPC testing of stingless bee propolis internationally [43].Phenolic compounds are components related to biological activity whose levels can be used as standards in propolis extract content, for example, phenolic compounds (chlorogenic acid, caffeic acid, isochlorogenic acid, caffeic acid phenethyl ester, and artepillin C) and flavonoids (myrisetin, quercetin, kaempferol, apigenin, pinocembrin, and galangin), which are used as standards in Brazilian and Chinese propolis [44].
These values are also associated with extraction techniques, solvents, and extraction times [45].Previous research has shown that obtaining a high TPC requires the propolis sample to be extracted using a robust extraction method, such as ultrasound-assisted extraction (UAE), with an extraction time of less than 0.5 h [45].In contrast, achieving a high TFC requires that samples should be extracted using conventional methods of extraction, such as maceration, which typically involves an extraction time of more than 0.5 h [45].The optimal solvent for obtaining propolis extractions with high TFC and TPC is aqueous ethanol [45].
The antioxidant activity with ABTS methods was in line with previous research on stingless bee propolis from Malaysia, which showed the highest antioxidant activity found in the G. thoracica species, with a value of 64.98 ppm [49].
Research on propolis extracts from Malaysia, the species Tetrigona apicalis, H. itama, and G. thoracica, shows a strong correlation between total phenolic and flavonoid levels and antioxidant activity [49].The biological activity of propolis is chiefly due to the presence of phenolic and flavonoid compounds, as well as variations in the quantities of these components [50].The polyphenol and flavonoid content of propolis has been found to have a linear correlation with its antioxidant activity [38,49], with higher levels of polyphenols and flavonoids associated with increased antioxidant activity.The study shows a linear relationship between antioxidant activity and TPC and TFC.
The antioxidant activity of phenolic compounds is associated with their ability to eliminate free radicals by donating hydrogen atoms, electrons, or metal cations; this interaction capacity with free radicals is due to their structure (mainly due to the number and position of hydroxyl groups and the nature of substitution in the aromatic ring) but is also based on the binding of the compound with organic acids and sugars [51].
In this case, a discriminant analysis method using Metaboanalyst 6.0 was applied to determine the classification based on the characteristics of propolis across different species and geographical regions.The discriminant analysis method in this study was carried out by observing the distribution area of the entire spectrum, which spanned 4000-500 cm −1 .Mapping from the FTIR test using the discriminant analysis method showed similar characteristics for propolis compounds divided into five groups.According to the results, the propolis groups from the H. itama sample displayed similar characteristics to those groups.The H. itama bee collects nectar, pollen, and resin to produce honey, bee bread, and propolis from various types of plants near the colony.Previous research has found that this species of bee prefers to visit certain types of vegetation [62].H. itama prefers to forage in areas closer to the nesting site, where diverse food sources are found.The marked bees of H. itama prefer to forage on various resources available within a 500 m radius of the house yard.The resin value of 34.73% was material that was brought back to the hive by foragers [63].
This also happened in the propolis group from T. laeviceps, which has a near score plot.The prevalence of resin from Pinus sp. in the nests of T. laeviceps was found to be a primary source in the production of propolis [64].This preference for Pinus sp.resin might explain why T. laeviceps frequently visit places that offer both resin and pollen as their food source.
On the other hand, the third group of FTIR maps consists of mixed propolis samples.They are not grouped by the same species.This may be caused by the source of vegetation or geographical origin.Propolis T. biroi, T. laeviceps, and T. drescheri from Batang have a near score plot of PLSDA.This phenomenon also occurred in propolis samples belonging to T. biroi and T. laeviceps from Magelang, as well as T. laeviceps and G. thoracica from Hulu Sungai Tengah.This phenomenon can be evaluated by the fact that T. biroi, G. thoracica, T. clypearis, and T. drescheri exhibit preferences in foraging behavior around their nests.These preferences are influenced by factors such as the location of food sources and the types of plants that can provide resin for producing propolis.

Propolis Chemical Compounds and Potential Markers
The chemical compounds in propolis include a rich variety of phenolic molecules.Often referred to as polyphenols, these compounds possess at least one aromatic ring and one or more hydroxyl functional groups.Flavonoids, which represent the most abundant group of phenolic compounds, have structures based on a C6-C3-C6 skeleton and are subdivided into several classes that differ in the oxidation state of the central heterocyclic ring.The classes are comprised of flavones, flavonols, flavanones, chalcones, isoflavonoids, flavanols (catechins and tannins), and anthocyanidins.Non-flavonoids include phenolic acids, simple phenols, coumarins, xanthones, lignins, lignans, and stilbenes.Phenolic acids can be divided into benzoic acid derivatives (based on a C6-C1 skeleton) and cinnamic acid derivatives (based on a C6-C3 skeleton).The variability of chemical composition in propolis can be attributed to large numbers of phenolics from different classes, such as glycoside phenolic compounds, and highlights the challenges associated with the analysis of propolis samples [43].
In stingless bee propolis from sources across the globe, investigations of chemical composition to date have identified many polyphenolic compounds.Based on their molecular structure, these compounds can be divided into flavonoids and nonflavonoids.The flavonoid compounds include flavone (chrysin), flavonol (galangin), and flavanone (pinocembrin).The nonflavonoid compounds are phenolic acids, which can be further subdivided into benzoic acid derivatives, such as gallic acid and protocatechuic acid, and cinnamic acid derivatives, such as caffeic acid, p-coumaric acid, and ferulic acid [65].
Based on analysis using LC-MS/MS and MetaboAnalyst software, 10 propolis compounds were identified as potential markers for stingless bees from H. itama and T. laeviceps propolis.In addition to the identification of those species classes based on the results of compound mapping from LC-MS/MS using PLSDA, they contain similar chemical compounds with similar characteristics.
The 10 molecular marker compounds of H. itama propolis identified in this study are petunidin, gingerone C, L-β-aspartyl-L-leucine, flazine, prostaglandine F2a, mannitol 1phosphate, cucurbic acid, and galacticol.On the other hand, pinusolide and hydroxyprolyl-Isoleucine were identified as high in T. laeviceps propolis.Petunidin, a blue antochianidine flavonoid, was found in the bee pollen of Echium plantagineum [66] and Fuchsia excorticata [67].Gingerenone C is a polyphenol found in Zingiber officinalle [68].Flazine is an indole alkalloid compound from Brucea javanica [69].Cucurbic acid is a natural compound found in Solanum tuberosum, Solanum lycopersicum, Peponapis limitaris, Cucurbita moschata, and other plants in the Cucurbitaceae [70].Mannitol and gallacticol are sugar alcohols that can be found in bee products [71].L-β-aspartyl-L-leucine and hydroxyprolyl-Isoleucine are amino acid derivatives that are needed for the development of bee colonies.Pinusolide is a diterpene lactone, a metabolite of Agathis macrophylla, Pinus armandii, and other Pinus species.This compound is high in T. laeviceps propolis, because it supplies resin and pollen as a bee food source [72].
Different kinds of plants have different resins with different chemical compositions, and even different individuals within the same species might have different resins.Even more widespread is resin utilization in stingless bees; several species of stingless bees gather large quantities of resin, which they then utilize to bolster various parts of colony function [24][25][26].

Samples
We studied a total of 23 raw propolis samples collected from 17 cities in Indonesia (Figure 13).The species were H. itama (HI), T. laeviceps (TL), T. biroi (TB), G. thoracica (GT), and others stingless bees, T. clypearis and T. drescheri, with a code of SB.The curator of the Zoologicum Museum, University of Indonesia (UIMZ), Biota Collection Room (RKBUI), Department of Biology, FMIPA, Universitas Indonesia, carried out the identification of bee species.We examined the morphology of bee species using the digital microscope DeltaPix DPX M12000 (DeltaPix, Smoerum, Denmark).Table 7 provides the list of samples.

Sample Preparation
The raw materials for propolis were obtained after harvesting honey by squeezing or sucking it from the beehives.After harvesting the honey, the bee bread, bee pollen, and dead bees were separated from the beehive.Then, propolis was frozen at −18-20 • C for a day until it reached a hard texture or was no longer soft.Frozen propolis was then crushed using a chopper to obtain fine flakes.Propolis in the form of fine granules could be used for extraction.

Propolis Extraction
A total of 50 g of fine propolis was macerated with 250 mL of 96% ethanol.The resultant mixture was stirred at 355 rpm for 8 h.The extract was subsequently left overnight and then filtered using a Buchner funnel (Deschem, Changshu, China) lined with filter paper with a pore diameter of 10 µm.The filtrate was designated as EEP 96% (ethanol extract propolis 96%), while the cake was designated as resin.After filtering, EEP 96% was diluted with distilled water until a concentration of 70% ethanol was reached.The diluted EEP was incubated for 30 min at 50 • C in a water bath.The solution was then stored in the refrigerator overnight at 0-5 • C. The wax precipitate was filtered using a Buchner funnel to obtain a filtrate, which was then concentrated using a rotary evaporator to obtain a thick ethanol extract of propolis.For extracts that were still not thick, the extract was dried in a food dehydrator (FDH) at a temperature of 50 • C for 24 h.The ethanol extract of propolis (EEP) obtained was stored in the refrigerator until used for testing [38].

Total Phenolic Content (TPC) Measurements
Quantitative tests of TPC were carried out using the 96-well microplate method specified in previous research [38].A standard of gallic acid (50 mg) was dissolved in methanol to obtain a concentration of 1000 ppm, which was then diluted with water to obtain series concentrations of 12.5-250 ppm.
A 25 µL sample of EEP or gallic acid standard, added with 100 µL of Folin-Ciocalteu reagent (25%), and then incubated for 4 min.Then, we added 75 µL of 1 M Na 2 CO 3 solution and shook for 60 s, before incubating for 2 h at room temperature.Absorbance was measured at λ 765 nm using a 96-well microplate reader spectrophotometer (Thermo Scientific Multiskan GO, Waltham, MA, USA).The absorbance standard solution of gallic acid was measured at each concentration, and then the calibration curve equation was calculated to obtain y = a + bx and the R 2 value, where y is the absorbance of gallic acid or sample and x is phenolic compound levels expressed in GAE (gallic acid equivalent).The measurements were carried out in triplicate.The total phenolic value was expressed in mg GAE/g extract.

Total Flavonoid Content (TFC) Measurements
Quantitative tests of TFC were carried out using the AlCl 3 method specified in previous research [38].A standard solution of quercetin (10 mg) was dissolved in methanol to obtain a concentration of 1000 ppm, which was then diluted with water to obtain series concentrations of 12.5-250 ppm.A 20 µL sample of EEP or quercetin standard at different concentrations was pipetted and then mixed with 10%AlCl 3 (20 µL), 1M CH 3 COOK (20 µL), and water (140 mL) and shaken for 60 s.The mixture was left to incubate for 30 min at room temperature.The absorbance of the propolis sample was subsequently measured at a wavelength of 415 nm using a 96-well microplate reader spectrophotometer (Thermo Scientific Multiskan GO, USA).The measurements were carried out in triplicate.
The absorbance standard solution of quercetin was measured at each concentration, and then the calibration curve equation was calculated to obtain y = a + bx and the R 2 value, where y is the absorbance of quercetin or sample and x is flavonoid compound levels expressed in mg QE/g extract.

Antioxidant Activity with DPPH Assay
Antioxidant activity was conducted using the previous method, as specified in Table 8 and a 96-well microplate reader spectrophotometer (Thermo Scientific Multiskan GO, USA).The percentage of DPPH inhibition was calculated using the formula in Equation ( 1).Antioxidant activity is expressed by the IC 50 obtained from the calibration curve.

Antioxidant Activity with ABTS Assay
ABTS•+ were generated according to the previous method of Wołosiak et al. [73], with slight modification by mixing equal volumes of substrate solution (2,2 ′ -azino-bis(3ethylbenzothiazoline-6-sulfonic acid, 7 mM) and oxidant (potassium persulfate, 2.45 mM); the reaction was carried out in the dark for 16 h.After this time, the radical solution obtained was diluted with ethanol.A total of 50 µL of EEP solutions (with a series concentration of 50-625 ppm) was collected in test tubes, and then 50 µL of the radical solution was added (Table 9).The reaction was performed for a standard time (6 min), measured at 734 nm.

Fourier-Transform Infrared Analysis
Unprocessed raw propolis samples were analyzed with a Fourier-transform infrared (FTIR) alpha spectrophotometer (Bruker Optik, Ettlingen, Germany).Stingless bee propolis samples were placed on the ATR crystal.All spectra were acquired between 4000-400 cm −1 , with a resolution of 4 cm −1 and a signal speed of 32 waves/minute.Each sample was tested in triplicate.Data point tables (DPTs) were used to store IR spectra in ASCII format, which were converted to Microsoft Excel tables in the form of (.csv) files.The spectrum data used was total transmittance information [74].

LC-MS/MS Analysis
The LC-MS/MS analyses were performed using UHPLC Vanquish Tandem Q Exactive Plus Orbitrap HRMS (Thermo Scientific, Waltham, MA, USA) in the Advanced Research Laboratory, IPB University.The signal read on the LC-MS/MS device was analyzed using Thermo Scientific Xcalibur 4.2 software to obtain a chromatogram.Data generated by the XCalibur software was further processed using Compound Discoverer 3.2 software.The column was Accucore C18, 100 × 2.1 mm, 1.5 µm (Thermo Scientific Waltham, MA, USA), the temperature was set at 30 • C, and the flow rate was 0.2 mL/min.The analyses were performed using a linear gradient solvent system consisting of A:B (0.1% formic acid in H 2 O:0.1% formic acid in acetonitrile) as follows: t = 0-1 min 5% B; t = 1-25 min 5-95% B; t = 25-28 min 95% B; t = 28-30 min 5% B. EEP samples were filtered using a 0.2 µm PTFE membrane filter (Merck KGaA, Darmstadt, Germany).The injection volume was 5.0 µL, and MSE positive ionization ESI electrospray had an acquisition range of 100-1500 m/z [38].The spectra of each sample were saved as (.raw) or (.MzML).The data acquisition and interpretation were processed using mzcloud, chemspider, and the Human Metabolome Database (HMDB).

Chemometric Analysis
Chemometric analysis with Metaboanalyst 6.0 software was conducted by employing principal component analysis (PCA), partial least squares discriminant analysis (PLS-DA), and sparse partial least squares discriminant analysis (sPLS-DA) to map the chemical compounds in stingless bee Indonesian propolis.The FTIR data on the (.csv) and LC-MS/MS data (.mzML) files were statistically analyzed.The data from 23 samples were preprocessed before chemometric analysis.After obtaining the normalized data, PCA, PLS-DA, and sPLS-DA were conducted for integrative data analysis [75].

Conclusions
The conclusions drawn in this study are based on the analysis of propolis compounds using FTIR and LC-MS/MS.Our findings are in line with the reports of previous scholars from different sources.The metabolite compounds were correlated with bee species and phytogeographic sources.The results of propolis compound mapping from the FTIR test using discriminant analysis demonstrated that the types of propolis from Java have similar characteristics.Then, based on species, the types of propolis from T. laeviceps and H. itama have special characteristics; the samples from each species can be grouped according to similar characteristics.Using LC-MS/MS analysis, we identified 10 propolis compounds as potential markers for Indonesian propolis from H. itama and T. laeviceps.
In conclusion, a metabolomic study was conducted on Indonesian stingless bee propolis,

Table 2 .
Results of antioxidant activity using DPPH and ABTS methods.

Table 3 .
Grouping of stingless bee propolis from Indonesia.

Table 4 .
Propolis compound results with FTIR.

Table 4 .
Propolis compound results with FTIR.

Table 4 .
Propolis compound results with FTIR.

Table 6 .
List of potential marker compounds of propolis from H. itama and T. laeviceps.

Table 6 .
List of potential marker compounds of propolis from itama and T. laeviceps.

Table 7 .
List of propolis samples.

Table 8 .
Procedure for testing antioxidant activity using the DPPH method.Shaken for 60 s and incubated for 40 min at dark room.Absorbance was measured at λ 516 nm.

Table 9 .
Procedure for testing antioxidant activity using the ABTS method.Shaken for 60 s and incubated for 6 min.Absorbance was measured at λ 734 nm.