Hostname: page-component-78c5997874-g7gxr Total loading time: 0 Render date: 2024-11-10T06:42:21.319Z Has data issue: false hasContentIssue false

The mountainous Cretan dietary patterns and their relationship with cardiovascular risk factors: the Hellenic Isolated Cohorts MANOLIS study

Published online by Cambridge University Press:  19 December 2016

Aliki-Eleni Farmaki*
Affiliation:
Department of Nutrition and Dietetics, School of Health Science and Education, Harokopio University, 70 El Venizelou Avenue, 17671 Athens, Greece
Nigel William Rayner
Affiliation:
Wellcome Trust Sanger Institute, Cambridge, UK Wellcome Trust Centre for Human Genetics, University of Oxford, Oxford, UK Oxford Centre for Diabetes, Endocrinology and Metabolism, Churchill Hospital, Headington, UK
Angela Matchan
Affiliation:
Wellcome Trust Sanger Institute, Cambridge, UK
Panagiota Spiliopoulou
Affiliation:
Department of Nutrition and Dietetics, School of Health Science and Education, Harokopio University, 70 El Venizelou Avenue, 17671 Athens, Greece
Arthur Gilly
Affiliation:
Wellcome Trust Sanger Institute, Cambridge, UK
Vasiliki Kariakli
Affiliation:
Department of Nutrition and Dietetics, School of Health Science and Education, Harokopio University, 70 El Venizelou Avenue, 17671 Athens, Greece
Chrysoula Kiagiadaki
Affiliation:
Anogia Medical Centre, Anogia, Crete, Greece
Emmanouil Tsafantakis
Affiliation:
Anogia Medical Centre, Anogia, Crete, Greece
Eleftheria Zeggini
Affiliation:
Wellcome Trust Sanger Institute, Cambridge, UK
George Dedoussis
Affiliation:
Department of Nutrition and Dietetics, School of Health Science and Education, Harokopio University, 70 El Venizelou Avenue, 17671 Athens, Greece
*
*Corresponding author: Email alelfarmaki@gmail.com
Rights & Permissions [Opens in a new window]

Abstract

Objective

We carried out de novo recruitment of a population-based cohort (MANOLIS study) and describe the specific population, which displays interesting characteristics in terms of diet and health in old age, through deep phenotyping.

Design

Cross-sectional study where anthropometric, biochemical and clinical measurements were taken in addition to interview-based completion of an extensive questionnaire on health and lifestyle parameters. Dietary patterns were derived through principal component analysis based on a validated FFQ.

Setting

Geographically isolated Mylopotamos villages on Mount Idi, Crete, Greece.

Subjects

Adults (n 1553).

Results

Mean age of the participants was 61·6 years and 55·8 % were women. Of the population, 82·7 % were overweight or obese with a significantly different prevalence between overweight men and women (43·4 v. 34·7 %, P=0·002). The majority (70·6 %) of participants were married, while a larger proportion of women were widowed than men (27·8 v. 3·5 %, P<0·001). Smoking was more prevalent in men (38·7 v. 8·2 %, P<0·001), as 88·8% of women had never smoked. Four dietary patterns emerged as characteristic of the population; these were termed ‘local’, ‘high fat and sugar, ‘Greek café/tavern’ and ‘olive oil, fruits and vegetables’. Individuals more adherent to the local dietary pattern presented higher blood glucose (β=4·026, P<0·001). Similarly, individuals with higher compliance with the Greek café/tavern pattern had higher waist-to-hip ratio (β=0·012, P<0·001), blood pressure (β=1·015, P=0·005) and cholesterol (β=5·398, P<0·001).

Conclusions

Profiling of the MANOLIS elderly population identifies unique unhealthy dietary patterns that are associated with cardiometabolic indices.

Type
Research Papers
Copyright
Copyright © The Authors 2016 

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.

Fig. 1 Recruitment flowchart

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.

Table 1 Anthropometric, clinical and dietary characteristics of the study population by sex: adults from isolated Mylopotamos villages, Crete, Greece, the Hellenic Isolated Cohorts MANOLIS study

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.

Table 2 Disease prevalence in the study population by sex: adults from isolated Mylopotamos villages, Crete, Greece, the Hellenic Isolated Cohorts MANOLIS study

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 %).

Table 3 Sociodemographic and lifestyle characteristics of the study population by sex: adults from isolated Mylopotamos villages, Crete, Greece, the Hellenic Isolated Cohorts MANOLIS study

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.

Table 4 Biochemical analyses according to low and high adherence to each dietary pattern in the study population: adults from isolated Mylopotamos villages, Crete, Greece, the Hellenic Isolated Cohorts MANOLIS study

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.

Table 5 Results from multiple linear regression analysis of the association between dietary patterns and cardiovascular risk factors among the study population (n 1183): adults from isolated Mylopotamos villages, Crete, Greece, the Hellenic Isolated Cohorts MANOLIS study

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

Footnotes

These authors contributed equally to this work.

References

1. World Health Organization (2015) Cardiovascular Diseases (CVDs). Fact Sheet no. 317. Geneva: WHO.Google Scholar
2. Zimmet, P, Alberti, KG & Shaw, J (2001) Global and societal implications of the diabetes epidemic. Nature 414, 782787.Google Scholar
3. Palacios, C, Perez, CM, Guzman, M et al. (2011) Association between adiposity indices and cardiometabolic risk factors among adults living in Puerto Rico. Public Health Nutr 14, 17141723.Google Scholar
4. Ashwell, M, Gunn, P & Gibson, S (2012) Waist-to-height ratio is a better screening tool than waist circumference and BMI for adult cardiometabolic risk factors: systematic review and meta-analysis. Obes Rev 13, 275286.Google Scholar
5. Kahn, R, Buse, J, Ferrannini, E et al. (2005) The metabolic syndrome: time for a critical appraisal: joint statement from the American Diabetes Association and the European Association for the Study of Diabetes. Diabetes Care 28, 22892304.Google Scholar
6. Grundy, SM (2007) Metabolic syndrome: a multiplex cardiovascular risk factor. J Clin Endocrinol Metab 92, 399404.Google Scholar
7. Yamaoka, K & Tango, T (2012) Effects of lifestyle modification on metabolic syndrome: a systematic review and meta-analysis. BMC Med 10, 138.Google Scholar
8. Visscher, PM, Brown, MA, McCarthy, MI et al. (2012) Five years of GWAS discovery. Am J Hum Genet 90, 724.Google Scholar
9. Speliotes, EK, Willer, CJ, Berndt, SI et al. (2010) Association analyses of 249,796 individuals reveal 18 new loci associated with body mass index. Nat Genet 42, 937948.Google Scholar
10. National Research Council, National Committee on Diet and Health (1989) Diet and Health: Implications for Reducing Chronic Disease Risk. Washington, DC: National Academies Press.Google Scholar
11. Keys, A, Menotti, A, Karvonen, MJ et al. (1986) The diet and 15-year death rate in the seven countries study. Am J Epidemiol 124, 903915.Google Scholar
12. Friedewald, WT, Levy, RI & Fredrickson, DS (1972) Estimation of the concentration of low-density lipoprotein cholesterol in plasma, without use of the preparative ultracentrifuge. Clin Chem 18, 499502.Google Scholar
13. Yannakoulia, M, Melistas, L, Solomou, E et al. (2007) Association of eating frequency with body fatness in pre- and postmenopausal women. Obesity (Silver Spring) 15, 100106.Google Scholar
14. Ainsworth, BE, Haskell, WL, Herrmann, SD et al. (2011) 2011 Compendium of Physical Activities: a second update of codes and MET values. Med Sci Sports Exerc 43, 15751581.Google Scholar
15. Bountziouka, V, Bathrellou, E, Giotopoulou, A et al. (2012) Development, repeatability and validity regarding energy and macronutrient intake of a semi-quantitative food frequency questionnaire: methodological considerations. Nutr Metab Cardiovasc Dis 22, 659667.Google Scholar
16. Ministry of Health and Welfare, Supreme Scientific Health Council (1999) Dietary guidelines for adults in Greece. Arch Hellenic Med 16, 516524.Google Scholar
17. Goldberg, GR, Black, AE, Jebb, SA et al. (1991) Critical evaluation of energy intake data using fundamental principles of energy physiology: 1. Derivation of cut-off limits to identify under-recording. Eur J Clin Nutr 45, 569581.Google Scholar
18. World Health Organization (1985) Energy and Protein Requirements. Report of a Joint FAO/WHO/UNU Expert Consultation. WHO Technical Report Series no. 724. Geneva: WHO.Google Scholar
19. Mardia, KV, Kent, JT & Bibby, JM (1979) Multivariate Analysis. New York: Academic Press.Google Scholar
20. Hu, FB (2002) Dietary pattern analysis: a new direction in nutritional epidemiology. Curr Opin Lipidol 13, 39.Google Scholar
21. Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol In Adults (2001) Executive Summary of The Third Report of The National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, And Treatment of High Blood Cholesterol In Adults (Adult Treatment Panel III). JAMA 285, 24862497.Google Scholar
22. Panoutsopoulou, K, Hatzikotoulas, K, Xifara, DK et al. (2014) Genetic characterization of Greek population isolates reveals strong genetic drift at missense and trait-associated variants. Nat Commun 5, 5345.Google Scholar
23. Tyrovolas, S, Zeimbekis, A, Bountziouka, V et al. (2009) Factors associated with the prevalence of diabetes mellitus among elderly men and women living in Mediterranean Islands: the MEDIS Study. Rev Diabet Stud 6, 5463.CrossRefGoogle ScholarPubMed
24. Kanoni, S & Dedoussis, GV (2008) Design and descriptive characteristics of the GHRAS: the Greek Health Randomized Aging Study. Med Sci Monit 14, CR204CR212.Google Scholar
25. Dedoussis, GV, Kanoni, S, Mariani, E et al. (2008) Mediterranean diet and plasma concentration of inflammatory markers in old and very old subjects in the ZINCAGE population study. Clin Chem Lab Med 46, 990996.Google Scholar
26. Kafatos, A, Diacatou, A, Labadarios, D et al. (1993) Nutrition status of the elderly in Anogia, Crete, Greece. J Am Coll Nutr 12, 685692.Google Scholar
27. Missoni, S (2009) Metabolic syndrome among the inhabitants of the island of Vis. Coll Antropol 33, 12811287.Google Scholar
28. Athyros, VG, Bouloukos, VI, Pehlivanidis, AN et al. (2005) The prevalence of the metabolic syndrome in Greece: the MetS-Greece Multicentre Study. Diabetes Obes Metab 7, 397405.Google Scholar
29. Panagiotakos, D, Pitsavos, C, Chrysohoou, C et al. (2009) Dietary patterns and 5-year incidence of cardiovascular disease: a multivariate analysis of the ATTICA study. Nutr Metab Cardiovasc Dis 19, 253263.Google Scholar
30. Panagiotakos, DB, Zeimbekis, A, Boutziouka, V et al. (2007) Long-term fish intake is associated with better lipid profile, arterial blood pressure, and blood glucose levels in elderly people from Mediterranean islands (MEDIS epidemiological study). Med Sci Monit 13, CR307CR312.Google Scholar
31. Moschandreas, J, Kafatos, A, Aravanis, C et al. (2005) Long-term predictors of survival for the Seven Countries Study cohort from Crete: from 1960 to 2000. Int J Cardiol 100, 8591.Google Scholar
32. Manios, Y, Antonopoulou, S, Kaliora, AC et al. (2005) Dietary intake and biochemical risk factors for cardiovascular disease in two rural regions of Crete. J Physiol Pharmacol 56, Suppl. 1, 171181.Google Scholar
33. Aravanis, C, Mensink, RP, Corcondilas, A et al. (1988) Risk factors for coronary heart disease in middle-aged men in Crete in 1982. Int J Epidemiol 17, 779783.Google Scholar
34. Kafatos, A, Schlienger, JL, Deslypere, JP et al. (1991) Nutritional status: serum lipids. Euronut SENECA investigators. Eur J Clin Nutr 45, Suppl. 3, 5361.Google Scholar
35. Sandker, GW, Kromhout, D, Aravanis, C et al. (1993) Serum cholesteryl ester fatty acids and their relation with serum lipids in elderly men in Crete and The Netherlands. Eur J Clin Nutr 47, 201208.Google ScholarPubMed
36. Sofi, F, Abbate, R, Gensini, GF et al. (2010) Accruing evidence on benefits of adherence to the Mediterranean diet on health: an updated systematic review and meta-analysis. Am J Clin Nutr 92, 11891196.Google Scholar
37. Trichopoulou, A, Costacou, T, Bamia, C et al. (2003) Adherence to a Mediterranean diet and survival in a Greek population. N Engl J Med 348, 25992608.CrossRefGoogle Scholar
38. Knoops, KT, de Groot, LC, Kromhout, D et al. (2004) Mediterranean diet, lifestyle factors, and 10-year mortality in elderly European men and women: the HALE project. JAMA 292, 14331439.Google Scholar
39. Trichopoulos, D & Lagiou, P (2001) Dietary patterns and mortality. Br J Nutr 85, 133134.Google Scholar
40. McCrory, MA, Howarth, NC, Roberts, SB et al. (2011) Eating frequency and energy regulation in free-living adults consuming self-selected diets. J Nutr 141, 148153.Google Scholar
41. Fung, TT, Rimm, EB, Spiegelman, D et al. (2001) Association between dietary patterns and plasma biomarkers of obesity and cardiovascular disease risk. Am J Clin Nutr 73, 6167.Google Scholar
42. Su, Q, Rowley, KG, Itsiopoulos, C et al. (2002) Identification and quantitation of major carotenoids in selected components of the Mediterranean diet: green leafy vegetables, figs and olive oil. Eur J Clin Nutr 56, 11491154.Google Scholar
43. Psaltopoulou, T, Naska, A, Orfanos, P et al. (2004) Olive oil, the Mediterranean diet, and arterial blood pressure: the Greek European Prospective Investigation into Cancer and Nutrition (EPIC) study. Am J Clin Nutr 80, 10121018.Google Scholar
44. Papamichael, CM, Karatzi, KN, Papaioannou, TG et al. (2008) Acute combined effects of olive oil and wine on pressure wave reflections: another beneficial influence of the Mediterranean diet antioxidants? J Hypertens 26, 223229.Google Scholar
45. Alonso, A, Ruiz-Gutierrez, V & Martinez-Gonzalez, MA (2006) Monounsaturated fatty acids, olive oil and blood pressure: epidemiological, clinical and experimental evidence. Public Health Nutr 9, 251257.Google Scholar
46. Ferrara, LA, Raimondi, AS, d’Episcopo, L et al. (2000) Olive oil and reduced need for antihypertensive medications. Arch Intern Med 160, 837842.Google Scholar
47. Gutierrez-Repiso, C, Soriguer, F, Rojo-Martinez, G et al. (2014) Variable patterns of obesity and cardiometabolic phenotypes and their association with lifestyle factors in the Di@bet.es study. Nutr Metab Cardiovasc Dis 24, 947955.CrossRefGoogle Scholar
48. Serra-Majem, L, Ribas, L, Tresserras, R et al. (1995) How could changes in diet explain changes in coronary heart disease mortality in Spain? The Spanish paradox. Am J Clin Nutr 61, 6 Suppl., 1351S1359S.Google Scholar
49. Walker, P, Rhubart-Berg, P, McKenzie, S et al. (2005) Public health implications of meat production and consumption. Public Health Nutr 8, 348356.Google Scholar
50. Katan, MB (1998) Effect of low-fat diets on plasma high-density lipoprotein concentrations. Am J Clin Nutr 67, 3 Suppl., 573S576S.Google Scholar
51. Michas, G, Micha, R & Zampelas, A (2014) Dietary fats and cardiovascular disease: putting together the pieces of a complicated puzzle. Atherosclerosis 234, 320328.Google Scholar
52. Mozaffarian, D, Katan, MB, Ascherio, A et al. (2006) Trans fatty acids and cardiovascular disease. N Engl J Med 354, 16011613.Google Scholar
53. Tachmazidou, I, Dedoussis, G, Southam, L et al. (2013) A rare functional cardioprotective APOC3 variant has risen in frequency in distinct population isolates. Nat Commun 4, 2872.Google Scholar
Figure 0

Fig. 1 Recruitment flowchart

Figure 1

Table 1 Anthropometric, clinical and dietary characteristics of the study population by sex: adults from isolated Mylopotamos villages, Crete, Greece, the Hellenic Isolated Cohorts MANOLIS study

Figure 2

Table 2 Disease prevalence in the study population by sex: adults from isolated Mylopotamos villages, Crete, Greece, the Hellenic Isolated Cohorts MANOLIS study

Figure 3

Table 3 Sociodemographic and lifestyle characteristics of the study population by sex: adults from isolated Mylopotamos villages, Crete, Greece, the Hellenic Isolated Cohorts MANOLIS study

Figure 4

Table 4 Biochemical analyses according to low and high adherence to each dietary pattern in the study population: adults from isolated Mylopotamos villages, Crete, Greece, the Hellenic Isolated Cohorts MANOLIS study

Figure 5

Table 5 Results from multiple linear regression analysis of the association between dietary patterns and cardiovascular risk factors among the study population (n 1183): adults from isolated Mylopotamos villages, Crete, Greece, the Hellenic Isolated Cohorts MANOLIS study

Supplementary material: File

Farmaki supplementary material

Supplementary Table

Download Farmaki supplementary material(File)
File 14.1 KB