CVD are the number one cause of death globally, representing 31% of all deaths. By 2030, almost 23·6 million people are predicted to die from CVD annually( 1 ). CVD is reaching epidemic proportions consistent with the worldwide sharp increase in diabetes and obesity( Reference Zimmet, Alberti and Shaw 2 ). Obesity relates to the total fat mass of an individual and is usually approximated by surrogate measurements such as BMI, waist circumference or waist-to-hip ratio (WHR). The latter two measures are more correlated to visceral fat content (central adiposity) than is BMI and are considered strong risk factors for diabetes and cardiometabolic risk( Reference Palacios, Perez and Guzman 3 , Reference Ashwell, Gunn and Gibson 4 ).
Risk factors associated with CVD include age, sex, ethnicity, family history, unhealthy diet, physical inactivity, smoking, alcohol consumption, obesity, hypertension, dyslipidaemia, impaired glucose and insulin metabolism, pro-thrombotic and pro-inflammatory state, and psychological and socio-economic factors( Reference Kahn, Buse and Ferrannini 5 , Reference Grundy 6 ). Interestingly, most of these are associated with lifestyle and can be modified. A recent meta-analysis provided strong evidence that long-term regular lifestyle modification significantly improves many of the intermediate CVD risk factors( Reference Yamaoka and Tango 7 ). In addition to environmental effects, a strong genetic component has also been shown to contribute to cardiometabolic disease risk( Reference Visscher, Brown and McCarthy 8 , Reference Speliotes, Willer and Berndt 9 ).
The Mediterranean diet represents a combination of many different foods in various consumption frequencies, combined with high levels of physical activity and moderate alcohol consumption, which suggests a more holistic approach to diet and lifestyle that can be incorporated in the term ‘pattern’. The ‘dietary pattern’ approach overcomes both the conceptual and methodological limitations of the ‘single nutrients’ approach, taking into account the antagonistic, interactive or synergistic relationships between nutrients, as people eat meals consisting of a variety of foods and not isolated nutrients( 10 ).
The Seven Countries Study was the first major study to investigate diet and lifestyle along with other risk factors for CVD( Reference Keys, Menotti and Karvonen 11 ). One of the regions studied within Greece was Crete, and the study included male residents of rural villages located inland from the northern coast, east of the capital of Heraklion. Cretan farmers were known for their longevity and held the title of the ‘gold standard’ of health status globally, with a very low prevalence of both CVD and cancer.
The inhabitants of mountainous Mylopotamos villages on Mount Idi, Crete, live in a geographically isolated high-altitude area and have preserved their customs, traditions and (generally) their anthropological type over the centuries. Anecdotally, this population is known to enjoy good health in old age despite a diet high in animal fat. In the present study, we established a well-phenotyped, population-based cohort in order to formally characterize the dietary and lifestyle patterns of this population and link those to traits of cardiometabolic relevance. In particular, we assessed the prevalence of diabetes and cardiovascular risk factors such as hypertension, hyperlipidaemia and obesity in this isolated population, analysed its particular dietary habits, and investigated the interactions among clinical, dietary and biochemical factors that may influence health.
Methods
Study population
The Minoan Isolates (MANOLIS) target population was adults originating from the mountainous region of Mylopotamos, on Mount Idi in Crete, which includes Anogia Municipality (2379 inhabitants), Zoniana Municipal Area (1117 inhabitants) and Livadia Local Community (1739 inhabitants). We contacted the local health centre of Anogia and with the cooperation of doctors and nurses from the regional clinics in each village, a total of 1553 volunteers were recruited between June 2010 and November 2012 (Fig. 1). According to the 2011 census, 29·7 % of the total population of the region participated in the MANOLIS study.
In regional medical centres, blood samples were collected from the antecubital vein by trained health professionals following blood pressure measurement. All blood draws were conducted in the morning after an overnight fast and compliance was high.
Serum and plasma were isolated following centrifugation at 3000 rpm for 10 min and stored at −80°C for future analyses. All participants were interviewed by trained personnel using a standard questionnaire. The questionnaire included clinical assessment and anthropometric measurement sections, as well as questions about demographic and lifestyle characteristics, detailed medical history, physical activity and dietary habits. The medical history and corresponding medication were cross-checked and confirmed by trained clinicians.
Anthropometric and biochemical measurements
The anthropometric measurements, including weight, height, waist and hip circumferences, were conducted by trained dietitians using standardized techniques and equipment. Body weight was measured using calibrated analogue scales. For height and other anthropometric measurements a stable tape was used.
BMI was calculated as weight divided by height2 (kg/m2). BMI was used to classify participants as underweight (BMI<18·5 kg/m2), normal weight (18·5≤BMI≤24·9 kg/m2), overweight (25·0≤BMI≤29·9 kg/m2) or obese (BMI≥30·0 kg/m2). WHR was also calculated.
Biochemical factors were assessed using enzymatic colorimetric assays and included glucose (hexokinase method), total cholesterol (cholesterol oxidase–phenol aminophenazone method), HDL cholesterol (HDL-C) and TAG (glycerol-3-phosphate oxidase–phenol aminophenazone method). Insulin was measured via chemiluminescence and LDL cholesterol (LDL-C) levels were calculated according to the Friedewald equation( Reference Friedewald, Levy and Fredrickson 12 ).
Demographic and lifestyle information
Information about the marital status (married, separated/divorced, single or widowed) of the participants was collected. Their educational level was defined as no education, primary, secondary or tertiary education, and the total number of years of education was recorded. Smoking habits were also collected. Participants were classified as current smokers, ex-smokers or never smoked. Pack-years was used as a unit for measuring the amount a person had smoked over a long period of time and was calculated by multiplying the number of packs of cigarettes (twenty cigarettes per pack) smoked per day by the number of years the person had smoked.
Physical activity assessment
Physical activity was evaluated using a brief self-reported validated questionnaire: the Harokopio Physical Activity Questionnaire (HAPAQ)( Reference Yannakoulia, Melistas and Solomou 13 ). The questionnaire is based on the metabolic equivalents of all activities of the preceding week, including activities at work, leisure time and rest or sleep, thus allowing the prediction of mean physical activity level and mean daily energy expenditure. The compendium of physical activities was used to quantify the intensity and consequently to calculate energy expenditure from the activities identified in the questionnaire using the method of metabolic equivalents( Reference Ainsworth, Haskell and Herrmann 14 ).
Based on the HAPAQ and its corresponding tool for estimating the necessary variables from the data obtained through the questionnaire, a script, the Physical Activity Program, was developed for time and memory saving purposes in order to derive these estimations directly and save them in one file (http://www.well.ox.ac.uk/~wrayner/tools/). This script was written in Perl and can be globally used to transform physical activity data into physical activity level, total energy expenditure and any other intermediate calculations desired.
Dietary assessment
Dietary information was collected via a validated, semi-quantitative, seventy-six-item FFQ( Reference Bountziouka, Bathrellou and Giotopoulou 15 ). All participants reported their daily, weekly or monthly average intake of several foods during the last year. Then the frequency of consumption was quantified on the basis of servings per week according to the dietary guidelines for adults in Greece( 16 ). Alcohol consumption was also assessed and calculated as drinks per day.
The above-mentioned transformations (foods and food groups in servings) were calculated programmatically.
Energy intake (EI) was extracted from the FFQ. The sex, weight and age data of participants were used to estimate their BMR with the Schofield equations and EI:BMR was calculated in order to identify and exclude energy under-reporters with the Goldberg et al. method( Reference Goldberg, Black and Jebb 17 , 18 ). Thus, a physiologically plausible EI:BMR value for our study is 1·15. Individuals whose EI:BMR was below this cut-off point were defined as under-reporters and excluded from subsequent dietary analyses (n final 1183).
Statistical analysis
Principal component analysis (PCA) was conducted to identify underlying dietary patterns( Reference Mardia, Kent and Bibby 19 ). In order for PCA to be effective in assessing food patterns, strong correlations between the predictive variables should exist. The correlation matrix of the food variables used in the present analysis showed that there were several correlation coefficients (r) >│0·4│, indicating that the variables were highly correlated. Moreover, the Kaiser–Meyer–Olkin test of sphericity and Barlett’s criterion was 0·76, implying high interrelationships between food variables and suitability of the data set for PCA. The orthogonal rotation (varimax option) was used to derive optimal non-correlated components (dietary patterns). From the entire database fifteen foods or food groups were used. To decide the number of components to retain the Kaiser criterion was used, according to which the number of components that should be retained from PCA is equal to the number of eigenvalues that are greater than 1. In our analysis, a four-component (food pattern) solution was selected. Based on the fact that factor loadings/correlation coefficients represent the correlation of each predicting variable with the dietary pattern score, higher absolute values indicate that the variable contributes more to the construction of this particular pattern. The dietary patterns were named according to scores of the predicting variables that correlated most with the component/pattern (>│0·4│). The median was used to divide pattern scoring into two groups: low and high. The PCA was performed using the statistical software package IBM SPSS Statistics for Windows, Version 21.0. The data handling, basic processing and descriptive characteristics analysis were carried out in R. Outliers were defined as values that exceed 3 sd above or below the mean and were removed from subsequent analyses.
Hypothesis testing was performed using the χ 2 test for categorical variables and Student’s t test or the Mann–Whitney for continuous variables in two different groups.
Multiple linear regression modelling was used to examine the associations between cardiometabolic indices and dietary patterns. In particular, three different models were applied: Model 1 was unadjusted; Model 2 included adjustment for age, sex, BMI (with the exception of BMI and WHR, in which analyses BMI was not a covariate), corresponding medication (for hypertension, hyperlipidaemia and diabetes), smoking in pack-years, physical activity level and EI; and Model 3, as per Model 2 with additional adjustment for the rest of the dietary patterns as these patterns are statistically independent, so it is possible for an individual to have high or low scores on more than one pattern at the same time and therefore all the patterns might act in concert( Reference Hu 20 ). The results from the linear regression models are presented as beta coefficients (β) and 95 % confidence intervals. The linear regressions were performed in IBM SPSS Statistics for Windows, Version 21.0. In all statistical analyses, the level of nominal significance was set at P=0·05.
Results
The anthropometric and clinical characteristics of the population are presented in Table 1. In total, 1553 individuals were recruited (55·8 % women and 44·2 % men) with a mean age of 61·6 years. Men had higher waist circumference (103·6 v. 97·0 cm) and WHR (0·98 v. 0·89) compared with women. The mean value for systolic blood pressure (SBP) was 138 mmHg, with men presenting significantly a higher level than women (140·1 v. 136·7 mmHg). Only 17·3 % of the participants were categorized as normal weight while the vast majority (82·7 %) were categorized as overweight or obese. Vegetables, refined cereals, fruits, meat and dairy were the food groups of the highest weekly consumption. Obesity, hypertension and hyperlipidaemia were the most prevalent diseases in this population (Table 2), with women showing higher morbidity than men.
WC, waist circumference; HC, hip circumference; WHR, waist-to-hip ratio; SBP, systolic blood pressure; DBP, diastolic blood pressure.
P values <0·05 are indicated in bold font; Student’s t test for continuous variables, χ 2 test for categorical variables.
P values <0·05 are indicated in bold font; χ 2 test.
The sociodemographic and lifestyle characteristics of the population are presented in Table 3. The majority (70·6 %) of participants were married, while a larger proportion of the women were widowed (27·8 %) compared with the men (3·5 %). The majority of participants reported primary education (63·0 %) and 10·6 % were illiterate. Men had significantly more total years of education than women (6·7 v. 5·5 years). A sex-specific difference was found for both current and former smokers, with the proportion of men being four and ten times higher, respectively, than that of women. Moreover, men were found to smoke for significantly longer, presented earlier starting and later quitting age points, and smoked more packs per year than women. Sixty per cent of the population had never smoked, the vast majority of whom were women (81·6 %).
P values <0·05 are indicated in bold font; χ 2 test for categorical variables, Mann–Whitney test for continuous variables.
From the initial fifty foods and food groups, fifteen were included in the PCA because of their high intercorrelation level. The Kaiser–Meyer–Olkin criterion was equal to 0·76 and the four different diet components explained 51·5 % of the total variance of the examined variables. Given that higher absolute values of the score coefficients derived from PCA (see online supplementary material, Supplemental Table 1) indicate that the food contributes more to the development of the component, the components were defined as follows: (i) the ‘local’ diet (component 1). which includes homemade pies, legumes, red meat, full-fat cheese and refined cereals; (ii) ‘high fat and sugar’ diet (component 2) that consists of chocolate, milky and starchy sweets, simple sugars, ready-made pies and cheese-flavoured snacks/chips; (iii) the ‘Greek café/tavern’ pattern (component 3) that is characterized by coffee and alcohol consumption, olives as a side for alcohol and fries; and (iv) the ‘olive oil, fruits and vegetables’ pattern (component 4).
The data from biochemical analyses showed that women had significantly higher mean serum levels of HDL-C than men, while men had significantly higher mean glucose, LDL-C and TAG levels compared with women. Notably, mean total cholesterol levels were above the normal range (>200 mg/dl)( 21 ) for this population (Table 4). The median was used to divide pattern scoring into two groups: low and high. The group with higher scores of the local pattern had higher mean glucose, while the groups with higher scores of the remaining patterns had lower glucose levels. Mean HDL-C was lower in the group having high scores of the Greek café/tavern pattern, while LDL-C and total cholesterol levels were found to be higher in the group with high scores of both the high fat and sugar and Greek café/tavern patterns. The group with higher scoring on the Greek café/tavern pattern presented higher TAG levels compared with the group with lower scoring on the same pattern and the group with higher scoring on the oil, fruits and vegetables pattern had higher insulin levels.
HDL-C, HDL cholesterol; LDL-C, LDL cholesterol; TC, total cholesterol.
P values <0·05 are indicated in bold font; Student’s t test for continuous variables.
* Low v. high compliance as defined by the median of the corresponding pattern score.
The results of the multi-adjusted linear regression analysis (Table 5) showed that component 1, representing the local typical diet, was initially associated with BMI, WHR, SBP and diastolic blood pressure (DBP), but the associations were attenuated after adjusting for confounding factors. The local typical diet was found to be consistently associated with serum glucose levels (β=4·861, P<0·001; β=4·618, P<0·001 after adjusting for confounders; β=4·026, P<0·001 after adjusting for confounders and the rest of the dietary patterns). Component 2 with high fat and sugar content was found to be negatively associated with BMI, WHR, HDL-C and SBP, but all these associations became non-significant after adjusting for confounding factors. In addition, the high fat and sugar dietary pattern was found to be associated with reduced glucose levels even after adjustment for confounders but not when the remaining diet was taken into account.
WHR, waist-to-hip ratio; HDL-C, HDL cholesterol; SBP, systolic blood pressure; DBP, diastolic blood pressure; TC, total cholesterol; LDL-C, LDL cholesterol. P values <0·05 are indicated in bold font.
* Model 1: unadjusted.
† Model 2: adjusted for age, sex, BMI (where applicable), medication (where applicable), physical activity level, smoking and energy intake.
‡ Model 3: adjusted for age, sex, BMI (where applicable), medication (where applicable), physical activity level, smoking, energy intake and the remaining food patterns scores.
Component 3, representing the Greek café/tavern pattern, was associated with reduced HDL-C but the effect was non-significant after adjusting for confounding factors. A consistent association of the Greek café/tavern pattern was found with increased WHR (β=0·020, P<0·001; β=0·012, P<0·001 after adjusting for confounders; β=0·012, P<0·001 after adjusting for confounders and the rest of the dietary patterns), DBP (β=0·903, P=0·002; β=0·846, P=0·014 after adjusting for confounders; β=1·015, P=0·005 after adjusting for confounders and the rest of the dietary patterns), total cholesterol (β=3·078, P=0·010; β=4·883, P=0·001 after adjusting for confounders; β=5·398, P<0·001 after adjusting for confounders and the rest of the dietary patterns) and LDL-C (β=3·977, P<0·001; β=3·933, P=0·001 after adjusting for confounders; β=4·401, P=0·001 after adjusting for confounders and the rest of the dietary patterns).
Component 4, capturing olive oil, fruit and vegetable consumption, was significantly associated with higher BMI (β=0·625, P<0·001; β=0·350, P=0·009 after adjusting for confounders; β=0·362, P=0·010 after adjusting for confounders and the rest of the dietary patterns), WHR (β=0·013, P<0·001; β=0·011, P<0·001 after adjusting for confounders; β=0·011, P=<0·001 after adjusting for confounders and the rest of the dietary patterns) and insulin levels, but not after adjustment for confounding factors for the latter association. Olive oil, fruit and vegetable consumption was found to be nominally associated with decreased SBP only after adjustment for confounders (β=−1·252, P 0·010) and taking into account both confounders and the rest of the diet (β=−1·146, P=0·025).
Discussion
In the present work we describe the collection of a population-based cohort from the isolated Mylopotamos villages on Crete, as part of the Hellenic Isolated Cohorts (HELIC) MANOLIS study (www.helic.org). The mean age of this genetically homogeneous( Reference Panoutsopoulou, Hatzikotoulas and Xifara 22 ) population from mountainous Mylopotamos was 61·6 years old, denoting an elderly population.
The MANOLIS population cohort was found to be overweight to marginally obese (mean BMI=29·5 kg/m2). This is in keeping with other Greek( Reference Tyrovolas, Zeimbekis and Bountziouka 23 , Reference Kanoni and Dedoussis 24 ) and European( Reference Dedoussis, Kanoni and Mariani 25 ) elderly populations. The prevalence of obesity was 44·0 % and women showed a higher tendency to be overweight than men. Generally, women were found to present higher morbidity than men, which could be attributed to the fact that the mean age for women was significantly higher. Of the population, 82·7 % was overweight or obese, higher than the corresponding percentage (57 %) from the same area (Anogia) in 1990( Reference Kafatos, Diacatou and Labadarios 26 ). One reason for this difference may be the study design as the previous study included only 167 elderly volunteers from Anogia, while our study assessed 1553 participants from different villages including Anogia, Zoniana, Livadia, Gonies, Axos, Astiraki, Kalivos and many others from the Mylopotamos area. Moreover, an increase in physical inactivity due to the changing modes of transportation over the years may also explain the observed shift, as the fundamental cause of obesity and overweight is the energy imbalance between energy consumed and energy expended. Hypertension affected almost half of the MANOLIS population (49·5 %), which is 20–25 % lower when compared with other Greek studies in the elderly( Reference Tyrovolas, Zeimbekis and Bountziouka 23 , Reference Kanoni and Dedoussis 24 ) and another isolated population( Reference Missoni 27 ), but 8 % higher than the prevalence reported by a study in one of the participating villages (Anogia), 20 years ago( Reference Kafatos, Diacatou and Labadarios 26 ). The Greek Health Randomized Aging Study (GHRAS) is a study of 782 elderly Greeks, randomly recruited in Athens, with the aim of investigating the interactions among socio-economic, lifestyle, biochemical and psychological factors determining the health status of elderly Greeks( Reference Kanoni and Dedoussis 24 ). The discrepancies for hyperlipidaemia (37·7 %) were in the same direction but of extended magnitude compared with GHRAS and the Mediterranean Islands Study (MEDIS), which found a higher prevalence of the disease (87 and 54 %, respectively). The prevalence of type 2 diabetes (12·3 %) in the MANOLIS population cohort is comparable with that in the general Greek population( Reference Athyros, Bouloukos and Pehlivanidis 28 ) and in other isolated cohorts like the one on the island of Vis in Croatia( Reference Missoni 27 ). However, the MANOLIS population does not suffer from associated complications of the disease such as neuropathy, nephropathy, eye disease, etc. Their most prevalent diseases can be grouped into two categories: one related to CVD risk (obesity, hypertension, hyperlipidaemia, type 2 diabetes) and the other related to degenerative diseases which have age-related pathologies (arthritis, presbyopia, osteoporosis, weak hearing, cataracts).
Meat intake of the MANOLIS population was three to eight times higher than that of the Attica general population( Reference Panagiotakos, Pitsavos and Chrysohoou 29 ) and other Greek elderly populations( Reference Tyrovolas, Zeimbekis and Bountziouka 23 , Reference Kanoni and Dedoussis 24 ). This is consistent with the geography of this mountainous isolated area and the tradition that has been preserved for centuries, with the local primary occupation being livestock farming. The increased alcohol consumption mainly by men reveals the social and cultural aspect of drinking. Men frequenting traditional cafés/bars drink and treat with their homemade alcoholic drink (raki) all visitors as a sign of hospitality.
The vast majority of the population was married (70·6 %), while less than 1 % were separated/divorced, demonstrating the traditional character of the population for which family is the core of social life, a characteristic that is weaker in populations living in more urban areas of Greece( Reference Kanoni and Dedoussis 24 , Reference Panagiotakos, Zeimbekis and Boutziouka 30 ). Most of the participants had attained primary education, with men more educated than women, a characteristic similar to an urban elderly population (GHRAS) but not to a younger urban Greek population( Reference Athyros, Bouloukos and Pehlivanidis 28 ) which has almost double the number of total years of education. This is a characteristic typical of rural areas, where agriculture is the main occupation. Educational attainment in the MANOLIS population is further curtailed by participant age and the effect of the onset of World War II.
Men were heavy smokers while the majority of women had never smoked. The percentage of male current smokers in MANOLIS is three times that in GHRAS (4·7 %), while the proportion of females who are current and ex-smokers is much lower than in GHRAS (6 v. 16 %)( Reference Kanoni and Dedoussis 24 ). Again this demonstrates the traditional nature of the particular society that preserves etiquettes and standards of the past. However, smoking in men was considerably lower (38 %) compared with 57 % in 1960, when a Cretan population was studied for the first time as part of the Seven Countries Study( Reference Moschandreas, Kafatos and Aravanis 31 ).
Serum glucose levels were the same as in GHRAS( Reference Kanoni and Dedoussis 24 ), while HDL-C levels in MANOLIS were found to be slightly (7–9 mg/dl) lower than those of other elderly Greek and European populations( Reference Kanoni and Dedoussis 24 , Reference Dedoussis, Kanoni and Mariani 25 ) although the levels appear consistent with previous studies in the area over the past 30 years( Reference Manios, Antonopoulou and Kaliora 32 – Reference Sandker, Kromhout and Aravanis 35 ). Our population from Mylopotamos had LDL-C levels that were 25–45 mg/dl lower compared with the mean level of both studies among the elderly in Attica( Reference Kanoni and Dedoussis 24 ) and Anogia( Reference Manios, Antonopoulou and Kaliora 32 ). Total cholesterol levels were found to be 20–40 mg/dl lower than those of other studies from Greece and the local area( Reference Kanoni and Dedoussis 24 , Reference Manios, Antonopoulou and Kaliora 32 , Reference Kafatos, Schlienger and Deslypere 34 , Reference Sandker, Kromhout and Aravanis 35 ), and resembled the levels found from the Seven Countries Study( Reference Moschandreas, Kafatos and Aravanis 31 ) in the 1960s. This highlights to some extent the Cretan paradigm of life expectancy and health.
The Mediterranean diet, one of the most studied dietary patterns, has been linked to increased survival, reduced risk of cancers and CVD and increased longevity( Reference Sofi, Abbate and Gensini 36 – Reference Knoops, de Groot and Kromhout 38 ), and has historically been derived from Crete, the origin of our study population. However, taking into account the geographical diversity of this island that affects the availability of food and living conditions, the dietary and lifestyle habits of a mountainous population may differ substantially from those of a valley population. In order to characterize the unique diet of this particular mountainous population from Mylopotamos, PCA, an a posteriori technique, was used to investigate patterns in the dietary data. A posteriori patterns are generated by empirically obtained data at hand( Reference Trichopoulos and Lagiou 39 ) and thus they do not necessarily identify ‘healthy’ patterns, but patterns that depict the commonly recognized eating habits of the particular population. The extracted components are based on subjective decisions. The derived patterns found were: the ‘local’ diet pattern, which represents a traditional diet of the native population; the ‘high fat and sugar’ diet pattern, which in addition to high sweet and fat consumption includes the preference for processed and pre-packed food; the ‘Greek café/tavern’ pattern that is characterized by coffee and alcohol consumption, fries and olives as side dishes for alcohol; and the ‘olive oil, fruits and vegetables’ pattern with olive oil, fresh fruit and vegetable consumption.
The associations of these patterns with several cardiovascular risk factors like BMI, WHR, total cholesterol, HDL-C, LDL-C, SBP, DBP, glucose and insulin were evaluated, after removing energy under-reporters. Under-reporting is more frequent among overweight individuals( Reference McCrory, Howarth and Roberts 40 ), like our study population.
Data analysis revealed no association of any of the derived dietary patterns with HDL-C and insulin. Others like Fung et al.( Reference Fung, Rimm and Spiegelman 41 ) observed a negative association between insulin and the prudent dietary pattern that included fruits, vegetables, poultry, fish, whole grains and legumes, and a positive association with the Western pattern diet characterized by higher consumption of red and processed meat, French fries, eggs, high-fat dairy products, sweets and refined grains.
The olive oil, fruits and vegetables diet pattern was negatively associated with SBP, suggesting that the antioxidant compounds( Reference Su, Rowley and Itsiopoulos 42 ) of this specific pattern might play a protective role. Olive oil has a dominant presence in the Mediterranean diet, which has been inversely related to SBP and DBP( Reference Psaltopoulou, Naska and Orfanos 43 , Reference Papamichael, Karatzi and Papaioannou 44 ). Similar findings have been published before( Reference Alonso, Ruiz-Gutierrez and Martinez-Gonzalez 45 , Reference Ferrara, Raimondi and d’Episcopo 46 ). However, the olive oil, fruits and vegetables diet pattern was also found to be associated with high BMI and WHR, an effect that could be attributed to the oil, which is the most energy-dense component of the pattern. When the compliance to this pattern (scores) was divided into tertiles, it was the medium tertile, not the higher tertile, that was significantly associated with increased risk for overweight/obesity. The latter might be due to the lack of discriminating capacity of this pattern considering the broadly high consumption of olive oil, fruits and vegetables in the area. In addition, it can be attributed to the fact that the majority of the specific population is overweight or obese, as well as to the previously mentioned protective role against other cardiometabolic risk factors such as blood pressure. In the Di@bet.es study( Reference Gutierrez-Repiso, Soriguer and Rojo-Martinez 47 ), a negative association was found between the prevalence of metabolic complications and olive oil intake, which explains the so-called ‘Spanish paradox’ where the morbidity and mortality rates for CVD are lower than expected from the high prevalence of metabolic complications, whether or not related to obesity( Reference Serra-Majem, Ribas and Tresserras 48 ).
The local diet pattern was associated only with higher glucose levels, in line with other studies where dietary patterns characterized by a higher intake of animal fat, mainly red meat and meat products, saturated fat, cholesterol and full-fat dairy products have been associated with increased CVD risk and death rate( Reference Panagiotakos, Pitsavos and Chrysohoou 29 , Reference Walker, Rhubart-Berg and McKenzie 49 ). This pattern was not found to be significantly associated with any other cardiometabolic risk factors, which seems surprising considering the unhealthy nature of the diet. This could be explained by the counterpoising effect of low carbohydrate and high fat content in increasing HDL-C levels( Reference Katan 50 ), an independent cardioprotective factor. Furthermore, the absence from the local pattern of any pre-packed or processed products reduces the load of trans-fatty acids. These compounds have been shown to have well-established adverse cardiovascular effects( Reference Michas, Micha and Zampelas 51 ) and decrease HDL-C( Reference Mozaffarian, Katan and Ascherio 52 ). In addition, the HELIC-MANOLIS population was found to have a low-frequency cardioprotective genetic variant associated with increased HDL-C levels, at a frequency forty times higher than the general Greek population( Reference Tachmazidou, Dedoussis and Southam 53 ). This shepherd population used to spend most of their time high in the mountains with their herds, so goat meat and milk/cheese were in abundance and in fact the only foods available to them to eat at times.
The Greek café/tavern pattern, a unique and very characteristic pattern for the culture of this population, was found to be associated with both BMI and WHR, demonstrating the relationship of this pattern with adiposity, as well as total cholesterol and DBP, indices that contribute more to the cardiometabolic risk. This finding is in absolute accordance with the established knowledge about cardiovascular risk factors( 1 ) taking into account the alcohol consumption as a component of this pattern, being in the same direction as observed in the ATTICA study where a diet pattern indicated by alcohol consumption was associated with CVD risk( Reference Panagiotakos, Pitsavos and Chrysohoou 29 ).
Our study had an observational, cross-sectional design and therefore shares the limitations of this approach. It is not appropriate to draw causal effect implications or to generalize based on the results from this particular Mylopotamos population. Residual confounding may also exist because of unused/unmeasured variables. Another limitation of the study is the self-reported medical history and medication data. In addition, the FFQ was not validated in this specific population, although it has been validated in the general Greek population( Reference Bountziouka, Bathrellou and Giotopoulou 15 ). The subjective nature of the decision making in components extraction is a limitation of PCA application. Lastly, we have not corrected our results for the number of tests carried out. We have tested nine outcomes of interest and hence the stringent Bonferroni-corrected P value is 0·0056 (=0·05/9). Using this significance threshold, the associations of the local dietary pattern with glucose levels, the Greek café/tavern pattern with WHR and cholesterol levels, and the olive oil, fruits and vegetables diet pattern with WHR, remain significant. The associations of the Greek café/tavern pattern with DBP and LDL-C levels, and of the olive oil, fruits and vegetables diet pattern with BMI and SBP, become suggestive as they no longer reach nominal significance.
Conclusion
The Cretan MANOLIS population is an isolated, elderly population with hypertension, obesity and hyperlipidaemia as the most prevalent diseases, the percentages of which resemble those of the general population. It can be described as a traditional society with family at the core of social life. The very characteristic dietary patterns, local and Greek café/tavern, identified in our study, are associated with overweight/obesity and other cardiovascular risk factors such as raised glucose, blood pressure and cholesterol levels, predicating their characterization as unhealthy. Consequently, the guidelines for decreasing cardiovascular risk should be applied to this population, through providing further information, raising awareness and policy making in these isolated areas. The burden of dietary compounds like high fat and red meat consumption in this population is not depicted in the intermediate cardiovascular manifestations to the expected extent, indicating compensatory procedures that require further investigation through gene–environment interaction analyses.
Acknowledgements
Acknowledgements: The authors thank the residents of the Mylopotamos villages for taking part. The MANOLIS study is dedicated to the memory of Manolis Giannakakis, 1978–2010. The MANOLIS study has been supported by many individuals who have contributed to sample collection (including Olina Balafouti, Christina Batzaki, Georgios Daskalakis, Eleni Emmanouil, Chrisoula Giannakaki, Margarita Giannakopoulou, Anastasia Kaparou, Vasiliki Kariakli, Stella Koinaki, Dimitra Kokori, Maria Konidari, Hara Koundouraki, Dimitris Koutoukidis, Eirini Mamalaki, Eirini Mpamiaki, George Pananikolaou, Maria Tsoukana, Dimitra Tzakou, Katerina Vosdogianni and Niovi Xenaki), data entry (Eleni Zengini, Thanos Antonos and Dimitra Papagrigoriou) and research administration (Anja Kolb-Kokocinski, Carol Smee and Danielle Walker). Financial support: This work was funded by the Wellcome Trust (grant number 098051) and the European Research Council (grant number ERC-2011-StG 280559-SEPI). The Wellcome Trust had no role in the design, analysis or writing of this article. The European Research Council had no role in the design, analysis or writing of this article. This work has taken place under the auspices of the Alexander S. Onassis Public Benefit Foundation Scholarship Program for Hellenes (grant number GZI 049/2012–2013). The Alexander S. Onassis Public Benefit Foundation had no role in the design, analysis or writing of this article. Conflict of interest: None. Authorship: E.Z. and G.D. contributed equally to this work. G.D. formulated the research question(s). E.Z. and G.D designed the study. A.-E.F., P.S., V.K., C.K. and E.T. carried out the study. A.-E.F., N.W.R., A.M. and A.G. analysed the data. A.-E.F. wrote the article. Ethics of human subject participation: This study was conducted according to the guidelines laid down in the Declaration of Helsinki and all procedures involving human subjects were approved by the Bioethics Committee of Harokopio University of Athens. Written informed consent was obtained from all participants.
Supplementary material
To view supplementary material for this article, please visit https://doi.org/10.1017/S1368980016003207