Adherence to Mediterranean Diet in Croatia: Lessons Learned Today for a Brighter Tomorrow

Non-communicable diseases (NCD) and lifestyle, particularly diet, have a close relationship. Based on the recent statistics, Croatian men and women lead in European overweight lists, which implies pessimistic prognosis in terms of incidence and prevalence of NCDs in the future. One of the possible solutions to overcome weight problems is turn to traditional balanced and sustainable diets, such as the Mediterranean diet. In this study, we assessed adherence towards Mediterranean diet using a validated questionnaire in an online survey and associated adherence scores with several demographic and anthropometric data. Based on the results of a validated Mediterranean Diet Adherence Screener (N = 3326), we assessed the adherence score to be 7.6 ± 2.5. The score tended to depend on sex, residence, age, education, income, and body mass index (BMI); indeed, women, residents of a coastal part of the country, older volunteers, those possessing a higher education degree, those with higher income, and those with lower BMI were associated with higher scores. As income was one of the significant findings related to higher adherence scores, we developed a dietary plan complying with Mediterranean diet principles that, on average, costed less than the average traditional balanced diet menu. Taken together, this study brought new findings regarding target groups who need to be encouraged to make lifestyle changes, and highlighted the first steps on how to make them.


Introduction
Non-communicable diseases (NCDs) are a leading cause of death globally, accounting for roughly 75% of deaths in 2019. This trend is more pronounced for high-income and upper-middle-income countries, where nine of the top ten estimated leading causes of deaths are attributed to NCDs [1]. The latest statistics indicate that, in 2017, 56 million people died, of which 48% were cardiovascular (17.8 mil) and cancer (9.6 mil) patients [2]. Regarding the incidence and outcome of those diseases, lifestyle has been identified as one of the key drivers. Smoking, excess body mass, lower physical activity, and poor nutrition are the most important constituents of lifestyle contributing to overall morbidity [3][4][5]. During the past 2 years, the COVID-19 pandemic emerged, accounting for more than 5 million deaths so far, and excess mortality of more than 80 deaths per 100,000 compared to the 5-year average in several countries [6]. Although COVID-19 is a communicable disease, a similar association was found, where worse lifestyle and obesity as a comorbidity led to worse disease outcomes [7,8]. Based on these data, we can suggest that targeting lifestyle is key to better disease prevention.
According to EUROSTAT, in 2019, Croatian males and females led the negative lists of prevalence of overweight individuals in Europe [9]. Moreover, recent reports suggest that Croatia has one of the highest prevalence in European countries for childhood overweight

Ethics
The Ethics committee of the Institute for Medical Research and Occupational Health approved the study, assuring the high principles of ethics in scientific research, privacy of results, and anonymity of participants.

Mediterranean Diet Adherence Questionnaire, Its Refinement, and Adaptation for Online Use
For the assessment of MedD adherence, we used a validated 14-point Mediterranean Diet Adherence Screener developed by Schröder et al. [29]. The questionnaire covers intake of olive oil, vegetables, fruit, red meat, fat, beverages, wine, legumes, seafood, pastry, nuts, and dishes based on tomato and olive oil, as well as preferences of olive oil, and white and red meat intake. Each question was scored with 0 or 1 based on participants' answers. The range of possible MedD score was 0 to 14.
The whole questionnaire was refined and adapted for online assessment using the Surveyplanet app using Pro features. The questionnaire refinement included food photographs to help the participants estimate portion size. Most of the questions were accompanied with a few photographs of representative food items (Figure 1). Representative food items and portion size were defined according to the national food consumption data [30]. The food photos were taken by a professional photographer using a digital camera (Canon) mounted on a tripod with constant lightning and angle according to guidelines by Ocké et al. [31].

Study Population
To recruit participants, the invitation link to participate in a study was sent using several social media platforms, non-governmental organizations' mailing lists, and posters and flyers with a QR code put up in public places, initiating non-probability snowball sampling. Before taking part in a study, participants read the informed consent and agreed to participate. The study was limited to one response per device, disabling multiple answers from the same person, including those who declined to participate in a study (N = 16). In addition to questions to assess MedD score, we collected basic anthropometric and socio-demographic data, as well as personal perception of adherence to MedD. The survey was active for a whole year (March 2020-March 2021).
To assure the reliability of the results, we approached participants who gave their consent and performed a repeatability study with the same set of questions. The repeatability study was performed at least 3 weeks from taking part in the original survey. Figure 1. Photographs of representative food items based on EFSA's national food consumption data and mass/volume for questions in the Mediterranean diet adherence questionnaire: (A) olive oil, (B) vegetables: chard, cucumber, pepper, carrot, broccoli, and lettuce, (C) fruit: blueberry and grapes, (D) meat: grilled minced meat "ćevapčići", Kranj sausage, and Frankfurt sausage, (E) butter, (F) sweetened beverage, (G) red wine, (H) legumes: beans, broad beans, lentils, and chickpeas, (I) seeds and nuts: sunflower seeds, pumpkin seeds, flax seeds, walnuts, hazelnuts, and almonds, (J) fish and seafood: gilt-head bream, hake, European pilchard, salmon fillet, and mixed seafood.

Study Population
To recruit participants, the invitation link to participate in a study was sent using several social media platforms, non-governmental organizations' mailing lists, and posters and flyers with a QR code put up in public places, initiating non-probability snowball sampling. Before taking part in a study, participants read the informed consent and agreed to participate. The study was limited to one response per device, disabling multiple answers from the same person, including those who declined to participate in a study (N = 16). In addition to questions to assess MedD score, we collected basic anthropometric and socio-demographic data, as well as personal perception of adherence to MedD. The survey was active for a whole year (March 2020-March 2021).
To assure the reliability of the results, we approached participants who gave their consent and performed a repeatability study with the same set of questions. The repeatability study was performed at least 3 weeks from taking part in the original survey.

Price Assessment for MedD
To determine average daily MedD cost, we used two different approaches. The first approach was to compare the prices of a MedD dietary pattern and a general healthy dietary pattern. Criteria for both dietary patterns were taken from the Dietary Guidelines for Americans 2020-2025 [32]. Diet plans for both dietary patterns were set to 2000 kcal per day. The following criteria were met for MedD: 184 g (6.5 oz) of protein foods, 27 g of oils, 170 g (6 oz) of grains, 2.5 cups of fruit and vegetables, and 2 cups of dairy per day. The healthy diet plan differed from MedD in the amount of fruit (2 cups vs. 2.5 cups), and dairy (3 cups vs. 2 cups). Additionally, the healthy diet plan differed in the lower total amount of protein food (5.5 oz vs. 6.5 oz in MedD), but also in the lower amount of protein food subgroup-seafood (8 oz vs. 15 oz in MedD per week). Average price per food group was estimated using data on the recommended intake according to the Dietary Guidelines (weight and frequency) and prices of selected locally consumed food. Locally consumed foods were selected from the Mediterranean food list for each food category [33]. For each food item, price per gram was multiplied by the grams that represent the portion (cups or ounce) of that food item. Price per portion of all food items in a certain food group was calculated, and consequently an average price for that food group was determined. Finally, the average price of a food group was multiplied by the recommended number of portions. We used prices from a nation-wide retailer providing the same food items and equal prices.
The second approach was to create 7-day MedD menus using the above-mentioned criteria, as well as the recommendation of weekly amounts for vegetable subgroups (darkgreen vegetables; red and orange vegetables; beans, peas and lentil; starchy vegetables and others) and protein foods subgroups (meat, poultry, eggs; seafood and nuts, seeds, soy products). The price for a 7-day MedD was averaged by the price for each day.

Statistical Analysis
Statistical evaluation was done using R statistical programming language version 4.1.0. Missing values from the data were imputed using the R package mice [34]. Basic statistical parameters were obtained using descriptive statistics. Basic statistical parameters were calculated by applying the basic statistic method and frequency tables. The correlation of various features was done by calculating Spearman's correlation. MedD scores were divided into quantiles and labelled Low (0-4), Intermediate (5-9), and High (10)(11)(12)(13)(14). The normality of continuous variables was checked using the Shapiro-Wilks W test. Afterwards, the difference between features between two quantiles was assessed by a Wilcoxon rank sum test, and among multiple groups by the Kruskal-Wallis test. A two-proportion Z-test was used to assess the difference in proportions within different features in each MedD score quantile. The influence of features on MedD scores was tested by the ordinal logistic regression model. The Wilcoxon signed-rank test was used to test the differences in MedD score before and after the initial study. Benjamini-Hochberg correction was applied to p-values, and the significance level in all of the tests was set to p < 0.05.

Results
The total study population consisted of 3342 volunteers, of which 3326 agreed to participate in the study (Croatian population in 2021 is 3.88 mil); 77.3% of volunteers were female, 22.7% male, while the average age was 37.4 ± 12 years. In total, 66.7% of volunteers possessed a university or higher degree, 33.3% had other educational levels. Half of the volunteers were never smokers, 23.2% were ex-smokers, and 26.8% were current smokers. The majority of participants were from the continental area, around 73%, while the rest were from the coastal area. The average MedD score for all participants was 7.6 ± 2.5.
The average MedD was significantly higher in participants from the coastal area, at 8.2 ± 2.3, compared to those from the continental area, at 7.4 ± 2.5, (Wilcoxon rank-sum test, p < 0.001). Higher proportions of participants from the coastal area were detected in higher MedD scores in comparison to the proportion of participants from the continental area ( Figure 2A), well above the calculated average MedD score. Women also had higher MedD scores when compared to men ( Figure 2B). The same trend was observed for participants above the median of 36 years of age ( Figure 2C). Higher proportions of participants whose income was in the upper two income quartiles, based on the National Bureau of Statistics data, had a MedD score above the median, while those whose income was in the lower two income quartiles tended to have a MedD score below the median ( Figure 2D).
Nutrients 2022, 14, x FOR PEER REVIEW 5 of 15 area (Figure 2A), well above the calculated average MedD score. Women also had higher MedD scores when compared to men ( Figure 2B). The same trend was observed for participants above the median of 36 years of age ( Figure 2C). Higher proportions of participants whose income was in the upper two income quartiles, based on the National Bureau of Statistics data, had a MedD score above the median, while those whose income was in the lower two income quartiles tended to have a MedD score below the median ( Figure  2D). Division of MedD scores into low (0-4), intermediate (5)(6)(7)(8)(9), and high (10-14) groups further validated these findings. There is a significant difference according to the different areas of residence, sex, and income in the high and low MedD score groups (two-proportion Z-test, p < 0.05) ( Figure 3A-C). On the other hand, for the intermediate MedD score groups, there was no significant difference according to any of the analyzed features. Age and body mass index (BMI) values were significantly different between three groups of MedD scores (Kruskal-Wallis test, p < 0.001) ( Figure 3D). Participants from the low MedD group were significantly younger compared to the intermediate and high MedD score group (one-sided Wilcoxon rank-sum test, p < 0.05). On the contrary, the BMI values were the highest in the low MedD group compared to the other MedD score groups (one-sided Wilcoxon rank-sum test, p < 0.05).
Spearman's correlation analysis showed a significant correlation between MedD score and age, BMI, level of education, income, smoking habits, and self-score in the complete cohort study (Figure 4). When analyzing the cohort separately based on the residency area, only education and smoking habit did not significantly correlate with the MedD score in participants from coastal parts of Croatia. The strongest significant and the most positive correlation is between self-score and MedD score (Spearman's correlation = 0.6), followed by correlation of income and education (Spearman's correlation = 0.44) in the complete cohorts. A similar trend was observed when separating the study group according to the residential area. Division of MedD scores into low (0-4), intermediate (5)(6)(7)(8)(9), and high (10-14) groups further validated these findings. There is a significant difference according to the different areas of residence, sex, and income in the high and low MedD score groups (two-proportion Z-test, p < 0.05) ( Figure 3A-C). On the other hand, for the intermediate MedD score groups, there was no significant difference according to any of the analyzed features. Age and body mass index (BMI) values were significantly different between three groups of MedD scores (Kruskal-Wallis test, p < 0.001) ( Figure 3D). Participants from the low MedD group were significantly younger compared to the intermediate and high MedD score group (one-sided Wilcoxon rank-sum test, p < 0.05). On the contrary, the BMI values were the highest in the low MedD group compared to the other MedD score groups (one-sided Wilcoxon rank-sum test, p < 0.05).
Spearman's correlation analysis showed a significant correlation between MedD score and age, BMI, level of education, income, smoking habits, and self-score in the complete cohort study (Figure 4). When analyzing the cohort separately based on the residency area, only education and smoking habit did not significantly correlate with the MedD score in participants from coastal parts of Croatia. The strongest significant and the most positive correlation is between self-score and MedD score (Spearman's correlation = 0.6), followed by correlation of income and education (Spearman's correlation = 0.44) in the complete cohorts. A similar trend was observed when separating the study group according to the residential area.      Due to the high correlation of these features, self-score and education were excluded from ordinal logistic regression for analysing the influence of different features on the Due to the high correlation of these features, self-score and education were excluded from ordinal logistic regression for analysing the influence of different features on the MedD score. The most significant features with the highest odds ratio detected with ordinal logistic regression are coastal areas (1.63, 95% CI [1.423-1.867]) and being female (1.343, 95% CI [1.156-1.559]) ( Figure 5). Only the smoking status did not have a significant influence on the MedD score. Participants who gave consent and their contact details (605 out of 3326) were asked to answer the same set of questions to assess the reliability of the results. We received 335 replies; of those, 29 people changed their dietary habits so they were excluded from the repeatability study, leaving us with a 50.6% response rate. No significant difference was detected between the MedD score obtained before and after the initial analysis (Wilcoxon signed-rank test, p = 0.73) ( Figure 6). Participants who gave consent and their contact details (605 out of 3326) were asked to answer the same set of questions to assess the reliability of the results. We received 335 replies; of those, 29 people changed their dietary habits so they were excluded from the repeatability study, leaving us with a 50.6% response rate. No significant difference was detected between the MedD score obtained before and after the initial analysis (Wilcoxon signed-rank test, p = 0.73) ( Figure 6). Due to the high correlation of these features, self-score and education were excluded from ordinal logistic regression for analysing the influence of different features on the MedD score. The most significant features with the highest odds ratio detected with ordinal logistic regression are coastal areas (1.63, 95% CI [1.423-1.867]) and being female (1.343, 95% CI [1.156-1.559]) ( Figure 5). Only the smoking status did not have a significant influence on the MedD score. Participants who gave consent and their contact details (605 out of 3326) were asked to answer the same set of questions to assess the reliability of the results. We received 335 replies; of those, 29 people changed their dietary habits so they were excluded from the repeatability study, leaving us with a 50.6% response rate. No significant difference was detected between the MedD score obtained before and after the initial analysis (Wilcoxon signed-rank test, p = 0.73) ( Figure 6).  In order to provide a weekly MedD dietary plan based on the calculation described in Section 2.3, and to target the low-income group as seen in Figures 2D, 3C and 5, we found that an average daily MedD plan per person is on average cheaper (EUR 9.47) compared to a traditional balanced plan (EUR 11.25). Moreover, in Figure 7, we managed to tailor a weekly MedD plan that would require an average of EUR 6.98 per day by choosing more affordable food items, thus potentiating its affordability. Figure 6. The results of the repeatability study with box plots and the distribution of Mediterranean diet score before and after the initial survey. No statistical differences were observed (Wilcoxon signed-rank test, p = 0.73).
In order to provide a weekly MedD dietary plan based on the calculation described in Section 2.3, and to target the low-income group as seen in Figures 2D, 3C and 5, we found that an average daily MedD plan per person is on average cheaper (EUR 9.47) compared to a traditional balanced plan (EUR 11.25). Moreover, in Figure 7, we managed to tailor a weekly MedD plan that would require an average of EUR 6.98 per day by choosing more affordable food items, thus potentiating its affordability.

Discussion
Weight control is an obvious global epidemiological problem for many countries, including Croatia [9, 35,36]. Energy-dense and nutrient-poor food items are characteristic of the Western diet, which is a key factor, along with a sedentary lifestyle, for obesity development. Studies found that obese individuals, among others, suffer from induced metabolic stress, which alters the gut microbiota, while a pro-inflammatory environment stimulates senescence. Therefore, being overweight or obese indicates a significant risk for development of many other non-communicable diseases, such as CVD, cancer, and diabetes [37][38][39]. Mediterranean diet, over the last few decades, has come into the spotlight due to having multiple beneficial properties to humans, including improvements in gut microbiota, vascular ageing, inflammatory and other health-related biomarkers, preventing certain cancer types, reduction in overall mortality, and combating CVD, thus potentially leading to substantial annual savings [19,21,[40][41][42][43][44][45][46][47][48][49][50]. In addition to many health and nutrition benefits, the Mediterranean diet, as a mainly plant-based diet, contributes to lower greenhouse emissions and a smaller water footprint, while promoting biodiversity. From 2010, UNESCO adopted MedD as an Intangible Cultural Heritage of Humanity, thus promoting the Mediterranean lifestyle as an invaluable socio-cultural achievement characterized by socialization, moderation, physical activity, and adequate rest. Finally, MedD contributes to positive local economic returns by protecting traditional production procedures and reducing the dependence on food import of local communities [20,25].
To enable rapid assessment of MedD within the population, we adapted a Mediterranean diet adherence questionnaire for online reporting. Since one of the major limitations in food consumption self-reporting is the over/underestimation of food portions, we

Discussion
Weight control is an obvious global epidemiological problem for many countries, including Croatia [9, 35,36]. Energy-dense and nutrient-poor food items are characteristic of the Western diet, which is a key factor, along with a sedentary lifestyle, for obesity development. Studies found that obese individuals, among others, suffer from induced metabolic stress, which alters the gut microbiota, while a pro-inflammatory environment stimulates senescence. Therefore, being overweight or obese indicates a significant risk for development of many other non-communicable diseases, such as CVD, cancer, and diabetes [37][38][39]. Mediterranean diet, over the last few decades, has come into the spotlight due to having multiple beneficial properties to humans, including improvements in gut microbiota, vascular ageing, inflammatory and other health-related biomarkers, preventing certain cancer types, reduction in overall mortality, and combating CVD, thus potentially leading to substantial annual savings [19,21,[40][41][42][43][44][45][46][47][48][49][50]. In addition to many health and nutrition benefits, the Mediterranean diet, as a mainly plant-based diet, contributes to lower greenhouse emissions and a smaller water footprint, while promoting biodiversity. From 2010, UNESCO adopted MedD as an Intangible Cultural Heritage of Humanity, thus promoting the Mediterranean lifestyle as an invaluable socio-cultural achievement characterized by socialization, moderation, physical activity, and adequate rest. Finally, MedD contributes to positive local economic returns by protecting traditional production procedures and reducing the dependence on food import of local communities [20,25].
To enable rapid assessment of MedD within the population, we adapted a Mediterranean diet adherence questionnaire for online reporting. Since one of the major limitations in food consumption self-reporting is the over/underestimation of food portions, we accompanied most of the questions with photographs of representative food items and their portion sizes. Online surveys have the ability to reach many people in a short period of time, with flexibility, and also low administration costs [51]; they proved to be quite useful, particularly in the period of heavy restrictions of mobility and accessibility due to the COVID-19 pandemic. Moreover, based on our results, the online survey proved to be reliable in assessing MedD score. On the other hand, females, people with a higher education degree, and those with high-income are more prone to participate in such a survey [52], which was also the case in our study.
After analysis of the results, we found that adherence to MedD is higher in people residing in coastal parts of Croatia, in women, in people older than 36, and in people with higher income. The score significantly correlated with education level and was inversely correlated with BMI. Such trends are in accordance with similar studies [28,47,[53][54][55][56][57][58][59][60][61][62][63][64][65][66][67]. Overall, the score of 7.6 ± 2.5 out of 14 could be considered as a medium MedD adherence level, which was in accordance with other Mediterranean [53][54][55][56][57] and non-Mediterranean countries [58,68,69]; however, certain specific populations might show poor to average adherence [70,71]. Still, we have to interpret results with the limitation that this study group was not representative for Croatia, since it comprised of > 75% women and > 50% of people with a university degree. Taking these factors into account, together with the results of the study by Pfeifer et al. [55] where it was concluded that Croatians tended to increase MedD adherence during pandemics, we could speculate that the actual MedD adherence score for Croatia is in fact lower than 7.6.
Concerning trends indicating abandoning MedD patterns have shown a significant decrease in Mediterranean adequacy index (MAI) over the last 50 years. On top of that, of top ten countries with biggest MAI decrease, nine were Mediterranean countries and Japan [27]. It is interesting to point out Japan, where the Okinawan diet is also characterized to be one of the drivers for longevity; however, the trends of westernization led to dietary changes, followed by a sharp increase in obesity and metabolic syndrome [72][73][74]. On the other hand, more optimistic trends are seen in the last decade, indicating the stabilization of MAI values which in some countries rose, while MedD adherence remained more or less the same [27,57,63].
When inspecting certain food items, we observed that 85.5% of volunteers did not meet the required red wine consumption, 69.9% of volunteers did not meet the required fish and seafood consumption, and 59.5% of volunteers did not meet the required olive oil intake. On the other hand, 83% of volunteers met the criteria for sweetened beverages intake, and 76.4% met the criteria for the intake of fat. These data are in line with the most recent paper investigating dietary preferences in the Croatian population [55]. In terms of health benefits, red wine is abundant in different polyphenols with potentially beneficial properties, while its consumption is associated with the "French paradox" and the reduction of coronary heart disease [75]. Fish and seafood is regarded as a source of high-quality protein and fat, particularly omega-3 fatty acids that were shown to reduce serum triglycerides and the risk of coronary heart disease [76]. Olive oil is highly valued because of its fatty acid and polyphenolic composition, which are key in risk reductions of CVD [77,78]. Based on the results of self-awareness towards complying with MedD patterns, the volunteers were quite accurate in their perception (2.96 ± 0.98 out of 5 (59%) vs. 7.6 ± 2.5 out of 14 (54%)), so we could conclude that they showed good general knowledge on how a MedD menu should be made.
If we go back to the turn of the 19th to 20th century, mainland Croatia's diet was based mostly on stews, soups, pork and its processed meat, poultry, dairy products, and cakes. Coastal Croatia's diet was based mostly on stews, a lot of vegetables, olive oil, fish, dairy products, and occasionally meat [79]. These traditional dietary patterns became less pronounced in younger generations [59,64,71,[80][81][82][83][84], which was also observed in our study. Interestingly, parents' diet also has a high impact on children's diet quality [85], which is key when advancing from planned kindergarten menus to primary schools where not all children follow the nutritionists'-tailored menus [65]. Therefore, particular groups, such as adolescents, students, and pregnant women, as well as men, should be considered when planning further awareness raising and educational campaigns, which should also be focused on the financial aspects of diet planning.
Based on our logistic regression analysis, we found that sex, residence, and income were major predictors for MedD score. Assuming that a person would like to keep her/his area of residence, and that sex is rather unchangeable, we wanted to target the lowincome population for refinement of their dietary plan based on the MedD principles. Our calculations actually showed that the MedD plan is on average slightly less expensive compared to the healthy omnivore plan; however, we did not compare it to unbalanced dietary plans that are often the reason for weight gain. Moreover, using the dietary plan tailoring, it is possible to eat according to MedD principles on a budget. Our calculations have the limitation of choosing prices from a nation-wide retailer, and some costs might be lower or higher in some other shops or during special offers. A major obstacle in dietary planning is the seasonality of fruit and vegetables, and although they are now available the whole year long, their price varies considerably in the winter period. Still, with a slight change in meat protein intake (pork instead of veal or mussels instead of prawns) and choosing different fruit (pears, oranges, kiwi, etc.) and vegetables (cabbage, broccoli, etc.), it is possible to make corrections within the plan, thus not altering the budget. Previous metaanalysis found that healthier food options costed slightly more (USD 0.29 per serving), and that the top quantile compared to the low quantile costed USD 1.54 per 2000 kcal more [86]. As a possible solution on how to tackle this issue, application of higher tax rates for "unhealthy" food items would contribute to price balancing, and potentially lead to substantial savings in healthcare system costs by preventing many NCDs [87]. Although nutrient-dense food is more expensive per calorie, recent findings [88][89][90] suggest that, with proper planning, it is possible to combine food items that are nutrient rich, affordable, and culturally acceptable. It is therefore necessary to initiate research that would tackle the assumptions that a "healthy" diet means higher expenses, and data presented in Figure 7 goes in line with such recommendations.

Conclusions
To conclude, in this study we showed that the Croatian population moderately practices adherence to Mediterranean diet. We also identified that more focus in promoting healthy aspects of MedD should be made towards men, people residing in the continental part of Croatia, adolescents and young adults, as well as toward people with incomes below the national median. Based on the available literature, the situation in Croatia reflects trends in other Mediterranean countries, as well as "Western countries", so similar approaches to improve nutrition quality could be implemented. By turning to MedD and improving lifestyle, it is possible to save money, prevent diseases, and even navigate ageing [91]; are not these enough to take action? Author Contributions: M.G., conceptualization, funding acquisition, investigation, methodology, writing-original draft; K.M., methodology, formal analysis, writing-review and editing, visualization; G.G., investigation, writing-review and editing; I.R., methodology, writing-review and editing, supervision; P.Š., methodology, formal analysis, writing-review and editing, visualization; R.K., methodology, formal analysis, writing-review and editing, supervision; M.B., investigation, methodology, writing-original draft, writing-review and editing. All authors have read and agreed to the published version of the manuscript. Informed Consent Statement: Online informed consent was obtained from all subjects involved in the study.

Data Availability Statement:
The original contributions generated for this study are included in the article. Further inquiries can be directed to the corresponding author upon reasonable request.