Skoči na glavni sadržaj

Izvorni znanstveni članak

From the Autochthonous Grape Varieties of the Kastav Region (Croatia) to the Belica Wine

Tomislav Pavlešić orcid id ; University of Rijeka, Faculty of Health Studies, Viktora Cara Emina 5, 51000 Rijeka, Croatia
Lara Saftić Martinović orcid id ; University of Rijeka, Department of Biotechnology, Radmile Matejčić 2, 51000 Rijeka, Croatia
Željka Peršurić orcid id ; University of Zagreb, Faculty of Chemical Engineering and Technology, Trg Marka Marulića 19, 10000 Zagreb, Croatia
Edi Maletić orcid id ; University of Zagreb, Faculty of Agriculture, Svetošimunska 25, 10000 Zagreb, Croatia
Maja Žulj Mihaljević orcid id ; University of Zagreb, Faculty of Agriculture, Svetošimunska 25, 10000 Zagreb, Croatia
Domagoj Stupić orcid id ; University of Zagreb, Faculty of Agriculture, Svetošimunska 25, 10000 Zagreb, Croatia
Željko Andabaka orcid id ; University of Zagreb, Faculty of Agriculture, Svetošimunska 25, 10000 Zagreb, Croatia
Zoran Grgić orcid id ; University of Zagreb, Faculty of Agriculture, Svetošimunska 25, 10000 Zagreb, Croatia
Sandra Kraljević Pavelić orcid id ; University of Rijeka, Faculty of Health Studies, Viktora Cara Emina 5, 51000 Rijeka, Croatia

Puni tekst: engleski pdf 860 Kb

str. 11-20

preuzimanja: 115


Preuzmi JATS datoteku

Prilozi: FTB-60-11_S1.pdf


Research background. Coastal region of Croatia is rich in autochthonous grape varieties. Many of them have been almost abandoned, such as the autochthonous varieties of Kastav (Croatia), used for the production of the Kastavska Belica wine. Therefore, the rationale of the presented study is to characterize autochthonous grape varieties Verdić, Mejsko belo, Jarbola, Divjaka and Brajkovac. In addition, we performed a molecular characterization of the corresponding Belica wines.
Experimental approach. Firstly, the genetic origin and ampelographic and economic characteristics of five autochthonous grape varieties were determined. Standard physicochemical profiles and phenolic components of 12 wines from different producers were determined by liquid chromatography coupled to triple quadrupole mass spectrometer (LC-QQQ-MS). Fourier-transform infrared spectroscopy (FTIR) was used for determination of standard physicochemical parameters.
Results and conclusions. Ampelographic analysis, which includes the data on producing characteristics and cluster and berry composition of the varieties, revealed significant differences between the analysed grape varieties. Results of the physicochemical analysis of the Belica wine showed that all wines met the requirements needed for the production of quality and top quality wines labelled with protected designation of origin (PDO) in Croatian coastal region. The LC-QQQ-MS analysis confirmed the presence of different phenolic components in the Belica wines, where the most prominent phenols were flavonoids from the flavan-3-ol group. Overall, these results showed that autochthonous grapes from the Kastav region can be used for production of wines with added market value due to a growing demand for autochthonous products on the global market.
Novelty and scientific contribution. The presented results give scientific insight and a basis for further determination of the optimal cultivation technology aimed to take advantage of the best characteristics of each variety for production of a wine with desirable features.

Ključne riječi

Belica wine, autochthonous wines, autochthonous grape varieties, polyphenols, FTIR

Hrčak ID:



Podaci na drugim jezicima: hrvatski

Posjeta: 711 *


A significant number of grape varieties deserve revitalization, due to their varietal characteristics that may also include resilience in the context of climate changes. In addition, the global market recognizes typical and autochthonous products, such as wines of indigenous grape varieties, often through high prices. In particular, EU has recognized the importance of traditional product sector, not only as a way to strengthen the local economy, but also as a way to generally develop a sustainability system (1).

Croatia is rich in indigenous grape varieties. Unfortunately, many of them are still neglected or scientifically uncharacterized. Successful example of a revitalized, nearly forgotten variety is Tribidrag (syn. Crljenak kaštelanski/Primitivo/Zinfandel/Kratošija), of which only 22 vines were found (2) near Kaštela. After revelation of its true identity, links to Plavac mali and eastern Adriatic origin (2), the renaissance and the resurgence of its production in Croatia began and increased demand for its planting material has been documented (3).

Malvazija istarska, Mejsko belo, Divjaka and Jarbola are also unique varieties, grown exclusively in Croatia, whilst synonyms of Verdić are spread over a wider area of the Northern Adriatic coast (4) under the names of Teran bijeli (5), Glera, Prosecco (6) and Beli Teran (7) in Slovenia and Prosecco tondo (8) in Italy. Apart from the vineyards, they are also conserved in the National Collection of Autochthonous Varieties at the University of Zagreb, Faculty of Agriculture, Croatia. However, that is not the case with Brajkovac, a variety mentioned back in 1853 (9), whose varietal status in terms of its uniqueness and possible synonyms/homonyms has not been evaluated before.

With these background data as a rationale for our study, we set the unique interdisciplinary approach for characterization of the autochthonous grape varieties and the corresponding wine. We tested our experimental approach to characterization of the autochthonous varieties of the Kastav region and corresponding Belica wine. Belica wine is a mixture of Mejsko belo, Verdić, Divjaka, Jarbola and Brajkovac grape varieties and belongs to the group of wine made from neglected and somewhat endangered varieties. Some of the varieties have indeed been on the verge of extinction. The Brajkovac variety occurs sporadically in some older vineyards and is used in a small percentage in Belica wine. The Jarbola variety is also present in small amounts. Varieties Divjaka and Mejsko belo are very important for the production of Belica wine and can be found exclusively in the vineyards of the Kastav region.

In this interdisciplinary approach, the ampelographic and economic analysis of varieties that can usually be found in the Belica wine, genetic background of the tested varieties and molecular components of the corresponding Belica wine will be determined. Standard wine chemical parameters are evaluated as well. The Fourier-transform infrared spectroscopy (FTIR) analysis was done as a confirmation of the results along with major phenolic compound evaluation by liquid chromatography coupled to triple-quadrupole mass spectrometry (LC-QQQ-MS).



The ampelographic and genetic research was conducted in the autumn of 2017 at the time of harvest. For wine analysis, total of 12 Belica wine samples (year 2017) were obtained from the local producers of the Kastav area (Croatia) (Fig. S1). The total number of vines ranged from 108–2050 pieces per vineyard of an individual producer. All producers use selected yeast (EC 1118, Saccharomyces cerevisiae) for the Belica wine production. The total area of the vineyards and the number of vines on which these twelve producers produce grapes for Belica wine is 13 442 m2 and 7285 vine pieces, respectively. The total annual production of Belica wine from these twelve producers is 6470 L. The predominant grape growing system (cultivation form) is single-legged or double-legged Guyot and the substrates are Kober 5 BB and SO4. The maximum yield per hectare for quality wines is 12 000 kg or 8400 L of wine per hectare, or for premium wines 11 000 kg or 6600 L of wine per hectare. Since each producer has different amounts of varieties in the vineyard and no standard guideline has been agreed on so far, it was not possible to determine the exact amount of each variety in the Belica wine.

Ampelographic and economic analyses of the varieties

Ampelographic analysis includes many parameters for grape variety characterization (morphology, phenology, production characteristic, etc.). This study included only parameters important for the production characteristics of the varieties. During the harvest, five healthy and vigorous vines were selected. Total yield and number of clusters were measured by picking, numbering and weighing the clusters from each vine. From the total mass of grapes, ten clusters were randomly sampled for further ampelographic analysis of cluster and berry composition. Cluster composition analysis encompassed measuring dimensions and mass of cluster according to Maletić et al. (10) and basic chemical composition of must. Berry composition parameters were chosen according to Rustioni et al. (11). A total of thirteen samples of berries were included for each variety. Berry composition analysis encompassed measurement of dimensions and mass of berry, and skin and seed mass. Skin and seeds of each sample were crushed, placed on paper, dried for two weeks at room temperature and then weighed.

The dimensions of clusters and berries were measured on graph paper. The mass of clusters and berries was measured with a precision laboratory balance (PS 4500.R2.M; Radwag, Radom, Poland). The basic chemical analysis of must comprised the analysis of sugar content (Brix scale, 2352 MASTER-53T refractometer; Atago, Tokyo, Japan) according to OIV-MA-AS2-02 method (12), total acid concentration (g/L expressed as tartaric acid equivalents) according to OIV-MA-AS313-01 method (13) and pH value (Lab 850 pH meter; Schott Instruments, Mainz, Germany) according to OIV-MA-AS313-15 method (14). These three parameters represent the most important quality parameters in wine production. Data obtained by weighing berries and their parts were used to calculate the mass of flesh and mass fraction of skin and seeds in berry composition. Parameter of yield per vine was used to calculate the economic value of grape production. First, grape production was expressed in yield per hectare and multiplied with average price for one kilogram of grapes. Economic analysis was expressed in total revenue, which included variable costs and gross margin.

Genetic analysis

For variety identification and confirmation, eight vines were sampled. DNA was extracted using the peqGOLD Plant DNA mini kit (PEQLAB Biotechnologie GmbH, Erlangen, Germany) according to the manufacturer's instructions. Nine microsatellite (simple sequence repeats, SSR) primers recommended for routine variety distinction of grapevine were used (15). Polymerase chain reaction (PCR) amplifications were carried out in an Applied Biosystems Veriti™ thermal cycler (Thermo Fisher Scientific, Foster City, CA, USA). The list and information about the used primers, as well as detailed information on multiplex PCR reactions performed, are described by Žulj Mihaljević et al. (4). Amplified products were separated using an Applied Biosystems 3130 genetic analyzer (Thermo Fisher Scientific) with GeneScan™ 500 LIZ® size standard. Sizing of the fragments was performed using GeneMapper v. 4.0 software (15). The obtained SSR profiles were compared to internal microsatellite database comprising profiles on 9 common loci from European Vitis database (16) as well as published SSR profiles from other research (4). Data were standardized and compared as described previously (4).

Reagents and materials used in wine analysis

Sodium potassium l(+)-tartrate tetrahydrate was obtained from VWR Chemicals (Vienna, Austria). Potassium iodide was obtained from BDH Prolabo Chemicals (Leuven, Belgium). Iodine, sodium hydroxide (1 and 0.1 M), sodium thiosulfate (0.1 M) and sodium hydroxide pellets were obtained from Gram-mol (Zagreb, Croatia). Sulphuric acid (96%), starch (p.a., soluble) and phenolphthalein were obtained from Kemika (Zagreb, Croatia). Bromothymol blue was obtained from Merck (Darmstadt, Germany). (+)-Catechin, (-)-epicatechin, 2,5-dihydroxybenzoic acid (2,5-DHBA), 3,4-dihydroxybenzoic acid (3,4-DHBA), 3-hydroxytyrosol, caffeic acid, ellagic acid, quercetin, naringenin, luteolin-7-O-glucoside, pinobanksin, p-coumaric acid and syringic acid were obtained from Sigma-Aldrich, Merck (St. Louis, MO, USA). Gallic acid was obtained from Alfa Aesar (Thermo Fischer Scientific, Tewksbury, MA, USA). Ferulic acid and resveratrol were obtained from Extrasynthese (Genay, France). Honeywell research chemicals (Charlotte, NC, USA) supplied ethanol (HPLC grade) and acetonitrile (LC-MS grade). Sigma-Aldrich supplied formic acid (LC-MS grade), ultrapure water (LC-MS grade) and ethanol (96%).

Standard wine analysis

Standard chemical parameters determined for the Belica wine samples were as follows: alcohol, reducing sugar and ash content, pH, total titratable and volatile acidity, and free and total sulphur dioxide content. Total alcohol content was determined by use of the electric ebulliometer (Exacta+Optech Labcenter S.p.A., San Prospero, Modena, Italy) (17). For the determination of pH according to OIV-MA-AS313-15 method (14), a pH meter Lab 860 (SI Analytics GmbH, Mainz, Germany) was used. Total acidity in the samples was determined according to modified OIV-MA-AS313-01 method (13). The only modification in the method was the use of the ultrasonic bath for the elimination of carbon dioxide from the wine instead of a vacuum flask and a water pump. The results were expressed as tartaric acid equivalents. The determination of reducing sugar content in wine samples was carried out according to the method developed by Rebelein (18). This is shortened iodometric method based on the titration of iodine formed in the reaction of potassium iodide and unused copper cation (left after the reaction of reducing sugars with alkaline copper sulphate) with sodium thiosulphate. The volatile acid in the samples (expressed as acetic acid equivalents) was determined according to the modified OIV-MA-AS313-02 method (19). The only modification in the method was the use of the ultrasonic bath for the elimination of carbon dioxide from the wine instead of a vacuum flask and a water pump. The ash in the wine samples was determined according to OIV-MA-AS2-04 method (20). Free and total sulfur dioxide were determined by the titration with a standard solution of iodine. Methods were developed according to OIV-MA-AS323-04B method (21) and the rapid method by Ripper (22). For free SO2 determination, 5 mL of diluted H2SO4 (1:3) with 2 mL of 1% starch were added to 50 mL of the sample and titrated with 0.01 M solution of J2 until a blue colour appeared. The consumption of J2 was multiplied by a factor of 12.8 and the results were expressed in mg/L of free SO2 in the sample. For total SO2 determination, 25 mL of 1 M NaOH were added to 50 mL sample and left to stand for 15 min. Afterwards, 10 mL of diluted H2SO4 (1:3) and 2 mL of 1% starch were added and titrated with 0.01 M J2 solution until a blue colour appeared. The titration consumption was multiplied by a factor of 12.8 and the results were expressed in mg/L of total SO2 in the sample.

Fourier-transform infrared spectroscopy analysis

Infrared spectra were recorded with WineScanTM FTIR spectrometer (FOSS, Hillerød, Denmark) within mid-IR (1000-5000 cm-1) range. The used samples were directly collected from the bottle without any pretreatment. Calibrations that are part of the WineScanTM FTIR allowed simultaneous analysis of major wine quality parameters such as alcohol, total acidity, volatile acidity and reducing sugars.

LC-QQQ-MS analysis

Wine samples were diluted twice with 10% ethanol solution, filtered through Chromafil cellulose acetate microfilters (0.45 µm, 25 mm; Macherey-Nagel, Düren, Germany) and analysed. LC-QQQ-MS analysis was performed with Agilent 1260 series HPLC chromatograph equipped with a degasser, binary pump, auto-sampler and column oven coupled to an Agilent 6460 triple quadrupole mass spectrometer equipped with jet stream electrospray (AJS ESI) source (Agilent Technologies, Palo Alto, CA, USA). For chromatographic separation, Zorbax SB-C18, rapid resolution HT, 6·107 Pa column (2.1 mm×50 mm i.d, 1.8 µm; Agilent Technologies) was used. The mobile phases were 0.1% formic acid in LC-MS grade water (A), and 0.1% formic acid in acetonitrile (B). Details of the method used for quantification of flavonoids and phenolic acids are described in our previous publication (23). Parameters for calibration curves of the analysed phenolic compounds (linearity, limit of detection (LOD), limit of quantification (LOQ) and coefficient of determination (R2) used for quantification of phenolic compounds are given in Table S1.

Statistical analysis

The obtained data were statistically processed using the SAS software, v. 9.3. (24). Statistical analysis included descriptive statistics (average, minimum and maximum value), analysis of variance (one-way ANOVA) and comparison of mean values (Duncan’s multiple-range test). Principal component analysis (PCA) was constructed using Python library Scikit-learn v. 0.20.3 (25) was used for both classifiers.


Ampelographic characteristics of varieties

All analysed varieties are white skin grapes with specific morphological characteristics. The results of one-way ANOVA of production characteristics showed significant differences in most of the parameters except for cluster mass, skin mass and mass fraction of skin and seeds (Table 1). Verdić had the largest cluster by dimensions (average length (188±23) mm, average width (116±15) mm), while Jarbola had the smallest clusters by dimension (average length (126±22) mm, average width (79±14) mm) (Table 2). Even though cluster mass was not significant, the ANOVA comparison of the mean values showed that a difference between varieties exists. Mejsko belo had the greatest and Jarbola the smallest cluster mass. Our analysis of the berry composition showed that Mejsko belo had the largest berry ((16.6±1.6) mm average length and (14.9±2.0) mm average width), while Jarbola had the smallest one ((13.9±1.3) mm average length and (12.5±1.2) mm average width) (Table 3). Variety Verdić also had the highest values for majority of other parameters: berry mass, flesh, skin and seed mass, and mass fraction of skin. All average values of the analysed varieties are very close to the average values of 22 383 data from the sample analysis in the study of Rustioni et al. (11). For example, the values of berry length/width ratio for Verdić and Brajkovac are the same as the average in the mentioned study. Varieties Mejsko belo and Brajkovac have the highest value of this parameter, which confirms the characteristic morphological ovoid shape of the berry.

Table 1 Results of one-way ANOVA for cluster and berry composition parameters
Table 2 Results of descriptive statistics (mean values with standard deviation) and comparison of mean values (Duncan’s multiple-range test) for cluster parameters of five autochthonous grape varieties of Kastav region (Croatia)*
Grape varietyl(cluster)/mmb(cluster)/mmm(cluster)/g
Mejsko belo(164±32)b(113±31)a(280±118)a

*Sample of 10 clusters. Different letters show a statistically significant difference ​​between varieties at p<0.05 (Duncan’s multiple range test)

Table 3 Results of descriptive statistics (mean values with standard deviation) and comparison of mean values (Duncan’s multiple-range test) for berry dimension parameters of five autochthonous grape varieties of Kastav region (Croatia)*
Grape varietyl(berry)/mmb(berry)/mml/bm(berry)/gm(flesh)/gm(skin)/gm(seed)/gw(skin)/%w(seed)/%
Mejsko belo(16.61±1.6)a(14.94±1.9)a(1.12±0.1)a(62.09±2.3)a(58.88±1.8)a(2.14±0.9)a(1.06±0.1)b(3.42±1.4)a(1.72±0.2)ab

*Sample of thirteen berry. Different letters show a statistically significant difference ​​between varieties at p<0.05 (Duncan’s multiple range test)

Considering the results of the analysed production characteristics of varieties, these also varied (Table 4). According to the cluster and berry analyses, Verdić showed to be the variety of large cluster, due to lower load (ten bunches per vine) and yield per vine (2.66 kg), and had the best quality of must. Additionally, both Verdić and Mejsko belo varieties had the highest sugar content (17.4 °Brix), and thus can be confirmed as the varieties with the best qualitative potential. Total acid concentration is very important quality parameter because it affects biochemical stability and organoleptic character of wine. Low acid concentration in must (3 to 5 g/L, like in Verdić and Mejsko belo) is not favourable from the technological aspect as it indicates artificial acidification of wine.

Table 4 Production and economical characteristics of five autochthonous grape varieties of the Kastav region (Croatia)
Grape variety(bunch)Total sugar/°Brixγ(total acids)/(g/L)pH/kgAverage price/(€/kg)Revenue/(€/ha)Variable costs/(€/ha)Gross margin/(€/ha)
Divjaka15. 889.605606.3618 283.24
Jarbola12. 610.325150.848003.96
Mejsko belo12.217.43.513.223.140.6612 278.975493.766672.61
Verdić10. 505.175248.085898.81
Brajkovac16.015.25.863.153.561.2025 311.605708.7219 705.24

Varieties Brajkovac and Divjaka in the observed year had the highest number of clusters and the highest yield per vine, over 3 kg. Such a high yield also affected the quality of the must, so both varieties had the lowest sugar content and the highest concentration of total acids (Table 4). On the other hand, high concentration of total acids, especially in Divjaka had a great variety potential for natural correction of acid concentration in Belica wine as a variety blend wine. In order to test the potential grouping of the samples and get deep insight into the differences between Belica wines, we performed PCA analysis (Fig. 1). Figs. 1a and 1b show the PCA projection of all analysed varieties and their potential grouping in the space of major components. The first two components (PC1 and PC2) describe 97.4% of the total variability. These results confirm a large difference between grape varieties in terms of standard chemical parameters.

Fig. 1 Distribution of elements (variables and samples) in the space of principal component 1 (PC1) and principal component 2 (PC2) when used as variables: a and b) standard physicochemical parameters of grape varieties, c and d) standard physicochemical parameters of Belica wines, e and f) phenolic compounds in Belica wines

Overall, from a technological point of view, in the future more attention should be paid to the reduction of yields of all varieties in order to achieve a higher concentration of sugar in the must but keep optimal acid concentration. Yield reduction in all varieties should be achieved by a combination of stronger pruning to maturity (leaving a smaller number of buds) and subsequent thinning of the clusters (after flowering or before the beginning of the ripening of clusters) and early defoliation. The results showed better insight into the agrobiological properties of these varieties, unexplored so far. However, we need to emphasize that this is a one-year study, and that results can differ depending on the year of production. Therefore, for the final confirmation of the obtained results, a multi-year research should be carried out under the above environmental conditions. The analysed varieties differ in their production characteristics (yields per vine and hectare), but also in economic characteristics (sale price, revenue and gross margin). The relationship between yield and selling price shows the highest profitability of growing varieties Brajkovac and Divjaka. Despite the relatively high yield of the variety Mejsko belo, it achieves lower revenues and coverage contribution due to the lowest selling price.

It can be assumed that the market design (branding) of Kastav Belica wine will contribute to the increase of the total income of wines of a mixture of autochthonous varieties, as well as individual varietiy wines.

Results of genetic analysis

Belica is a blend (cuvée or mixture) of several grape varieties. The following grape varieties can be found in Belica, most of which are found only in the Kastav region: Mejsko belo, Verdić, Divjaka, Jarbola, Malvazija istarska and Brajkovac. Among them, four (Mejsko belo, Verdić, Malvazija istarska and Divjaka) are the most represented. Malvazija istarska was not analysed in this paper, as it is not autochthonous variety of Kastav area, but of Istrian peninsula. All eight accessions were successfully amplified (Table 5). As expected, true-to-type nature of Mejsko belo, Verdić, Divjaka and Jarbola was confirmed after comparison with internal, previously harmonized SSR database (4) containing more than 2000 nonredundant grape genotypes. Four vines were sampled and assumed under the name Brajkovac; however, three different genetic profiles were obtained. Accessions labelled BRAJ_ORIG1 and BRAJ_ORIG2 showed to be identical to Duranija and Mejsko belo, respectively. Nevertheless, two accessions (BRAJ_A and BRAJ_B) showed identical profiles on all nine analysed SSR markers. The obtained genetic profile did not match any other previously known variety/genotype from the SSR database nor has this name appeared in foreign literature, thus confirming the unique status of this variety. These two accessions were further considered as true Brajkovac and were subject of further analysis. This result enables next conservation steps needed for the preservation of this variety, like its inclusion in National collection of autochthonous grapevine varieties settled at the Faculty of Agriculture in Zagreb, Croatia. Also, financial support of local authorities for propagation of the planting material will be beneficial for producers that will accordingly be able to plant those varieties in future.

Table 5 Microsatellite profiles and their genetic match on 9 SSR loci for eight analysed accessions. Alleles are presented as base pairs
BRAJ_ORIG2*141143237245175177187193234256222224239253247256272272Mejsko belo
Mejsko belo141143237245175177187193234256222224239253247256272272Mejsko belo

*The autochthonous grape variety Brajkovac samples

Belica wine analysis

According to physicochemical parameters, all analysed Belica wines meet the requirements for quality wine produced under the label protected designation of origin (PDO) on the Croatian coast (Table 6). FTIR analysis confirmed these results (Table S2). Particularly, the alcohol volume fractions in Belica wine samples were between 11.4 and 13.5% (Table 6). As for the sugar concentration, the results show that all the analysed samples are within the limits of dry wine production (1.05-3.89 g/L). Ash concentration was in the range 0.8 to 1.9 g/L. According to the regulation established for PDO Croatian coastal region (26), the minimum concentration of ash for white wines in the category of quality wine with controlled geographical origin is 1.4 g/L. From the obtained results, it is evident that all except one (sample B10) analysed wine samples meet the given criteria. Results of the analysis of total acids in wine showed that all the analysed samples are above the minimum legal limit of wine production (3.5 g/L), whereas the lowest measured concentration of total acids in Belica wine samples was 4.50 g/L. Furthermore, the upper limit for the volatile acid concentration was 1.1 g/L. It is evident that all analysed wines meet the given criteria. Results of the determination of the free sulfur in wine showed that in samples B1 to B4, B6 and B12, the measured concentration was too low, and the wines were in a state of oxidation. After the analysis of total sulfur in wine, all relevant samples of Belica wine complied with the regulation on wine, whereas the upper limit for white and rosé wines is 200 mg/L of free SO2 (Table 6). The PCA projection of the standard physicochemical parameters (Figs. 1c and 1d) shows moderate variability among Belica wine samples. The first two components (PC1 and PC2) describe 61.6% of the total variability. These results show that a more uniform production of Belica wine should be pursued.

Table 6 Results of standard analyses of Belica wine samples from the Kastav region (Croatia)
Wine sampleφ(alcohol)/%γ(sugar)/(g/L)γ(ash)/(g/L)pHγ(acidity)total/(g/L)γ(acidity)volatile/(g/L)γ(SO2)free/(mg/L)γ(SO2)total/(mg/L)

During the winemaking process, a number of chemical modifications that occur significantly affect the final phenolic profile of the wine, for example, grape ripeness, processing methods and environmental factors. Therefore, a systematic quantitative analysis of phenolic components in wine can provide reliable data on their quantity and type. During the handling and ripening of the grapes, the composition of the polyphenols in the wine may change. Phenolic analysis is therefore crucial to draw conclusions about the winemaking process, as well as about the final wine quality.

Phenolic content in Belica wines was analysed by LC-QQQ-MS. The results showed that the most common groups of polyphenols are hydroxycinnamic acids and flavan-3-ols, which is in line with previous research on white wines (2729). Among all hydroxycinnamic acids, caffeic acid was present at highest concentrations, ranging from (0.8±0.4) to (9.1±0.1) mg/L (Table 7). Lukić et al. (30) conducted qualitative research of phenols in different white wines and reported different trends. Namely, in their research, the phenolic acid with the highest concentration was the gallic acid, with concentrations up to (16.68±15.30) mg/L in Muscat Blanc wine. Our results report significantly lower concentrations of this phenolic acid in the analysed samples ((0.39±0.00) to (1.9±1.1) mg/L). However, the results of the caffeic and ferulic acid concentrations were similar to the results of Lukić et al. (30). Rochetti et al. (31) reported lower concentrations of caffeic acid in Chardonnay wines (up to (0.26±0.15) mg/L), which is more in line with our results. Along with ferulic acid, their results for syringic and ferulic acid concentrations were also in line with our data.

Table 7 Mass concentration of specific phenolic acids and flavonoids in Belica wine samples obtained by LC-QQQ-MS method
2,5-DHBA3,4-DHBAcaffeic acidellagic acidferulic acidgallic acidp-coumaric acidsyringic acid(+)-catechin(-)-epicatechin3-hydroxy

Results are expressed as a mean value±standard deviation (S.D.), N=2

As for flavan-3-ols, catechins were found in the highest concentrations, with (+)catechin being dominant ((4.4±0.1) to (17.0±2.2) mg/L) (Table 7). The catechins are responsible for the bitterness in wine (32-34) but also for wine health properties (35). In addition to sensory properties, they are also important as antioxidants, i.e. factors that protect the wine from oxidation during maturation (36, 37). These results are in line with the work by Lukić et al. (30), who reported concentrations of epicatechin and catechin in white wines in the range from (3.10±2.60) to (17.92±13.10) and from (1.51±1.82) to (3.54±2.99), respectively. Rochetti et al. (31) reported lower concentrations of catechin in Chardonnay wines, ranging from (1.19±0.49) to (6.81±2.53) mg/L. In addition, we also reported higher hydroxytyrosol content in Belica wines, up to (2.3±0.3) mg/L, while the highest concentration in Chardonnay wines was (0.9±0.2) mg/L.

Figs. 1e and 1f show the PCA projection of LC-QQQ-MS quantitative phenolic analysis where all analysed phenols were used as variables. The first two components describe 47.3% of the total variability. Although most wine samples are grouped centrally, samples B5, B6, B10 and B12 contribute to a larger variability of the system. Large differences may be observed in single polyphenol concentrations among samples, showing a need for the development of a more uniform method of Belica wine production.


The tested interdisciplinary approach for characterization of selected autochthonous grape varieties and corresponding wine samples proved useful in the assessment of important parameters for branding and quality assessment. Analyses of ampelographic characteristics of Verdić, Mejsko belo, Jarbola, Divjaka and Brajkovac varieties used for the production of Belica wine in the region of Kastav (Croatia) showed high genetic variability among the confirmed grape varieties. Confirmed variability certainly impacts the productivity and economic aspect of the Belica wine production. Results of this study are crucial in the determination of an optimal cultivation technology as it is required by modern trends in grape and wine production. Microsatellite genetic profile and uniqueness of the new variety Brajkovac was confirmed. The knowledge of standard wine characteristics coupled with molecular analyses may be used to evaluate the best characteristics of each variety and establish the production of a wine with desirable characteristics. Recording and monitoring of typical molecular composition of wine through different years, often called molecular profiling may be an important tool for standardization and/or monitoring of the technological process of the production, protection and branding of this autochthonous Belica wine. Current global market trends, indeed, emphasize local, specific and autochthonous products, increasing their demand. Currently, there are large differences in the production characteristics of tested varieties, which is determined by differences in yield and selling prices. Further work on branding the Kastavska Belica wine might increase this specific wine quality and provide benefits to the producers.


Supplementary materials are available at


We appreciate the access to the equipment in possession of University of Rijeka within the project ’Research Infrastructure for Campus-based Laboratories at University of Rijeka’, financed by the European Regional Development Fund (ERDF). We also greatly acknowledge the Association Belica and wine producers of the Town Kastav who supported this research and kindly provided samples for experimental analyses. In addition, we thank Ivana Lučin, University of Rijeka, Faculty of Engineering, for assistance in Phyton-based principal component analysis.


[1] Financial disclosure FUNDING

Funding sources are University of Rijeka research support uniri-biomed-18-133 as well as the project ‘Mass spectrometry characterisation of Belica wine’ financed by the town of Kastav to S.K.P.

[2] Conflicts of interest CONFLICT OF INTEREST

The authors have no conflict of interest to declare.



A sustainable bioeconomy for Europe - Strengthening the connection between economy, society and the environment: Updated bioeconomy strategy. Luxembourg: European Commission, Directorate-General for Research and Innovation, Publications Office; 2018.


Maletić E, Pejić I, Karoglan Kontić J, Piljac J, Dangl GS, Vokurka A, et al. Zinfandel, Dobričić, and Plavac Mali: The genetic relationship among three cultivars of the Dalmatian coast of Croatia. Am J Enol Vitic. 2004;55(2):174–80.


Andabaka Ž, Stupić D, Marković Z, Preiner D. New trends in plant material production of autochthonous grapevine cultivars. Glas Zast Bilja. 2011;34(1):46–56. [in Croatian]


Žulj Mihaljević M, Maletić E, Preiner D, Zdunić G, Bubola M, Zyprian E, et al. Genetic diversity, population structure, and parentage analysis of Croatian grapevine germplasm. Genes (Basel). 2020;11(7):737. PubMed:


Maletić E, Sefc MK, Steinkellner H, Karoglan Kontić J, Pejić I. Genetic characterization of Croatian grape cultivars and detection of synonymous cultivars in neighboring regions. Vitis. 1999;38(2):79–83.


Štajner N, Korošec-Koruza Z, Rusjan D, Javornik B. Microsatellite genotyping of old Slovenian grapevine varieties (Vitis vinifera L.) of the Primorje (coastal) winegrowing region. Vitis. 2008;47(4):201–4.


Štajner N, Rusjan D, Korošec-Koruza Z, Javornik B. Genetic characterization of old Slovenian grapevine varieties of Vitis vinifera L. by microsatellite genotyping. Am J Enol Vitic. 2011;62:250–5.


De Lorenzis G, Imazio S, Biagini B, Failla O, Scienza A. Pedigree reconstruction of the Italian grapevine aglianico (Vitis vinifera L.) from Campania. Mol Biotechnol. 2013;54:634–42. PubMed:


Sabljar M. Names of grape vine varieties around Rijeka (Fiume) II. Gospodarske novine. 1853;4:15. (in Croatian).


Maletić E, Karoglan Kontić J, Preiner D, Šimon S, Staver M, Pejić I. Ampelographic and genetic studies into ‘Teran’/‘Refošk’ grapes in Istria (Croatia) – One or two varieties? Mitt Klosterneubg. 2014;64(2):54–62.


Rustioni L, Cola D, Maghradze D, Abashidze E, Argirou A, Aroutiounian R, et al. Description of the Vitis vinifera L. phenotypic variability in eno-carpological traits by a Euro-Asiatic collaborative network among ampelographic collections. Vitis. 2019;58(1):37–46.


OIV-MA-AS2-02. Evaluation by refractometry of the sugar concentration in grape musts, concentrated grape musts and rectified concentrated grape musts (Oeno 466/2012). Compendium of international methods of wine and must analysis, vol. 1. Paris, France: International Organisation of Vine and Wine (OIV); 2021. pp. 61-82.


OIV-MA-AS313-01. Total acidity (Oeno 551/2015). Compendium of international methods of wine and must analysis, vol. 1. Paris, France: International Organisation of Vine and Wine (OIV); 2021. pp. 433-5.


OIV-MA-AS313-15. pH (A31 revised by Oeno 438-2011). Compendium of international methods of wine and must analysis, vol. 1. Paris, France: International Organisation of Vine and Wine (OIV); 2021. pp. 491-3.


GeneMapper® software, v. 4.0, reference and troubleshooting guide, Applied Biosystems, Foster City, CA, USA; 2005. Available from:


Maul E, Sudharma KN, Kecke S, Marx G, Müller C, Audeguin L, et al. The European Vitis database ( - A technical innovation through an online uploading and interactive modification. Vitis. 2012;51(2):79–85.


Sotlář R, Lampíř L. Effect of yeasts on aromatic profiles of wines from ‘Cabernet’ grapevine cultivars. Kvasny Prum. 2017;63(3):139–47.


Rebelein H. Rapid method for the determination of the alcohol sugar and the total sulfur dioxide contents by distillation in wine and fruit juices and also for determining blood alcohol. Chem Mikrobiol Technol Lebensm. 1973;2:112–21.


OIV-MA-AS313-02. Volatile acidity. Compendium of international methods of wine and must analysis, vol. 1. Paris, France: International Organisation of Vine and Wine (OIV); 2021. pp. 436-40.


OIV-MA-AS2-04. Ash. Compendium of international methods of wine and must analysis, vol. 1. Paris, France: International Organisation of Vine and Wine (OIV); 2021. pp. 93-4.


OIV-MA-AS323-04B. Sulfur dioxide (iodometry). Compendium of international methods of wine and must analysis, vol. 2. Paris, France: International Organisation of Vine and Wine (OIV); 2021. pp. 424-6.


Ripper M. The sulphurous acid in wine and its determination (Communication from Schmitt’s laboratory, [examination office, chemical test station and hygienic institute] in Wiesbaden.). J Prakt Chemie. 1892;46(1):428-73 (in German).


Saftić L, Peršurić Ž, Fornal E, Pavlešić T, Kraljević Pavelić S. Targeted and untargeted LC-MS polyphenolic profiling and chemometric analysis of propolis from different regions of Croatia. J Pharm Biomed Anal. 2019;165:162–72. PubMed:


SAS. v. 9.3, SAS Institute Inc, Cary, NC, USA; 2012. Available from:


Pedregosa F, Varoquaux G, Gramfort A, Michel V, Thirion B, Grisel O, et al. Scikit-learn: Machine learning in Python. J Mach Learn Res. 2011;12:2825–30.


Regulation (EU) No 1308/2013 of the European Parliament and of the Council of 17 December 2013 establishing a common organisation of the markets in agricultural products and repealing Council Regulations (EEC) No 922/72, (EEC) No 234/79, (EC) No 1037/2001 and (EC) No 1234/2007. Off J EU. 2013;L347:671-854. Available from:


Alvarez-Casas M, Pajaro M, Lores M, Garcia-Jares C. Polyphenolic composition and antioxidant activity of Galician monovarietal wines from native and experimental non-native white grape varieties. Int J Food Prop. 2016;19(10):2307–21.


Somers TC, Vérette E. Phenolic composition of natural wine types. In: Linskens HF, Jackson JF, editors. Wine analysis. Modern methods of plant analysis, vol 6. Berlin, Germany: Springer; 1988. pp. 219–57.


Hernanz D, Recamales ÁF, González-Miret ML, Gómez-Míguez MJ, Vicario IM, Heredia FJ. Phenolic composition of white wines with a prefermentative maceration at experimental and industrial scale. J Food Eng. 2007;80(1):327–35.


Lukić I, Radeka S, Budić-Leto I, Bubola M, Vrhovsek U. Targeted UPLC-QqQ-MS/MS profiling of phenolic compounds for differentiation of monovarietal wines and corroboration of particular varietal typicity concepts. Food Chem. 2019;300:125251. PubMed:


Rocchetti G, Gatti M, Bavaresco L, Lucini L. Untargeted metabolomics to investigate the phenolic composition of Chardonnay wines from different origins. J Food Compos Anal. 2018;71:87–93.


Fernandes I, Pérez-Gregorio R, Soares S, Mateus N, de Freitas V. Wine flavonoids in health and disease prevention. Molecules. 2017;22(2):292. PubMed:


Oliveira CM, Silva Ferreira AC, de Freitas V, Silva AMS. Oxidation mechanisms occurring in wines. Food Res Int. 2011;44(5):1115–26.


Kallithraka S, Bakker J, Clifford MN. Evaluation of bitterness and astringency of (+)-catechin and (-)-epicatechin in red wine and in model solution. J Sens Stud. 1997;12(1):25–37.


Plavša T, Jurinjak N, Antunović D, Peršurić Đ, Kovačević Ganić K. The influence of skin maceration time on the phenolic composition and antioxidant activity of red wine Teran (Vitis vinifera L.). Food Technol Biotechnol. 2012;50(2):152–8.


Tarko T, Duda-Chodak A, Sroka P, Siuta M. The impact of oxygen at various stages of vinification on the chemical composition and the antioxidant and sensory properties of white and red wines. Int J Food Sci. 2020;2020:7902974. PubMed:


Nikolantonaki M, Julien P, Coelho C, Roullier-Gall C, Ballester J, Schmitt-Kopplin P, et al. Impact of glutathione on wines oxidative stability: A combined sensory and metabolomic study. Front Chem. 2018;6:182. PubMed:

Floating objects

Table S1 LC-QQQ-MS parameters for phenolic acid and flavonoid analysis of Belica wine samples
Compoundt/minPolaritym/z precursor ionm/z product ion*Collision energy/VY=ax+bLinearity range/(µg/mL)R2LOD/(µg/mL)LOQ/(µg/mL)
Slope (a)Intercept (b)
Phenolic acid

*Quantifier ions are underlined. LOD=limit of detection, LOQ=limit of quantification, DHBA=dihydroxybenzoic acid

Table S2 Results of the Fourier-transform infrared spectroscopy (FTIR) analysis of Belica wine samples from Kastav region (Croatia)
φ(alcohol)/%γ(glucose and
γ(total acidity)/
γ(malic acid)/
γ(lactic acid)/
γ(volatile acidity)/(g/L)pHSpecific
Fig. S1 Map of the Kastav area (Croatia). Red circles indicate the locations of the vineyards from which Belica grapes were collected

This display is generated from NISO JATS XML with jats-html.xsl. The XSLT engine is libxslt.