Obesity indexes and total mortality among elderly subjects at high cardiovascular risk: The PREDIMED study

Obesity Indexes and Total Mortality among Elderly Subjects at High Cardiovascular Risk: The PREDIMED Study Miguel A. Martı ´nez-Gonza ´lez 1,2,3 *, Ana Garcı ´a-Arellano 1,2,3 , Estefanı ´a Toledo 1,2,3 , Maira Bes-Rastrollo 1,2,3 ,Mo ´nica Bullo ´ 1,2,4 , Dolores Corella 1,5 , Montserrat Fito 1,6 , Emilio Ros 1,7 , Rosa Maria Lamuela-Ravento ´s 1,2,8 , Javier Rekondo 1,2,9 , Enrique Go ´mez-Gracia 1,2,10 , Miquel Fiol 1,11 , Jose Manuel Santos-Lozano 1,12 , Lluis Serra-Majem 1,2,13 , J. Alfredo Martı ´nez 1,2,14 , Sonia Eguaras 1,3 , Guillermo Sa ´ez-Tormo 1,2,15 , Xavier Pinto ´ 1,2,16 , Ramon Estruch 1,2,17 1Centro de Investigacio ´n Biome ´dica en Red de Fisiopatologı ´a de la Obesidad y Nutricio ´n, Instituto de Salud Carlos II, Madrid, Spain, 2The PREDIMED (Prevencio ´ncon Dieta Mediterra ´nea) Network (RD 06/0045) Instituto de Salud Carlos III, Madrid, Spain, 3Department of Preventive Medicine and Public Health, University of NavarraOsasunbidea, Servicio Navarro de Salud, Pamplona, Spain, 4Human Nutrition Department, Hospital Universitari Sant Joan, Institut d’Investigacio ´Sanitaria Pere Virgili, Universitat Rovira i Virgili, Reus, Spain, 5Department of Preventive Medicine, University of Valencia, Valencia, Spain, 6Cardiovascular and Nutrition Research Group, Institut de Recerca Hospital del Mar, Barcelona, Spain, 7Lipid Clinic, Department of Endocrinology and Nutrition, Institut d’Investigacions Biome `diques August Pi i Sunyer, Hospital Clinic, University of Barcelona, Barcelona, Spain, 8Department of Nutrition and Food Science, School of Pharmacy, Xarxa de Refere `ncia en Tecnologia dels Aliments, Instituto de Investigacio ´n en Nutricio ´n y Seguridad Alimentaria, University of Barcelona, Barcelona, Spain, 9Department of Cardiology, University Hospital of Alava, Vitoria, Spain, 10 Department of Preventive Medicine, University of Malaga, Malaga, Spain, 11 Institute of Health Sciences, University of Balearic Islands, and Hospital Son Espases, Palma de Mallorca, Spain, 12 Department of Family Medicine, Primary Care Division of Seville, San Pablo Health Center, Seville, Spain, 13 Department of Clinical Sciences, University of Las Palmas de Gran Canaria, Las Palmas, Spain, 14 Department of Nutrition and Food Sciences, Physiology and Toxicology, University of Navarra, Pamplona, Spain, 15 Department of Biochemistry and Molecular Biology-Service of Clinical AnalysisUniversity of Valencia, Valencia, Spain, 16 Lipids and Vascular Risk Unit, Internal Medicine, Hospital Universitario de Bellvitge, Hospitalet de Llobregat, Barcelona, Spain, 17 Department of Internal Medicine Institut d’Investigacions Biome `diques August Pi i Sunyer, Hospital Clinic, University of Barcelona, Barcelona, Spain Abstract Background: Different indexes of regional adiposity have been proposed for identifying persons at higher risk of death. Studies specifically assessing these indexes in large cohorts are scarce. It would also be interesting to know whether a dietary intervention may counterbalance the adverse effects of adiposity on mortality. Methods: We assessed the association of four different anthropometric indexes (waist-to-height ratio (WHtR), waist circumference (WC), body mass index (BMI) and height) with all-cause mortality in 7447 participants at high cardiovascular risk from the PREDIMED trial. Forty three percent of them were men (55 to 80 years) and 57% were women (60 to 80 years). All of them were initially free of cardiovascular disease. The recruitment took place in 11 recruiting centers between 2003 and 2009. Results: After adjusting for age, sex, smoking, diabetes, hypertension, intervention group, family history of coronary heart disease, and leisure-time physical activity, WC and WHtR were found to be directly associated with a higher mortality after 4.8 years median follow-up. The multivariable-adjusted HRs for mortality of WHtR (cut-off points: 0.60, 0.65, 0.70) were 1.02 (0.78–1.34), 1.30 (0.97–1.75) and 1.55 (1.06–2.26). When we used WC (cut-off points: 100, 105 and 110 cm), the multivariable adjusted Hazard Ratios (HRs) for mortality were 1.18 (0.88–1.59), 1.02 (0.74–1.41) and 1.57 (1.19–2.08). In all analyses, BMI exhibited weaker associations with mortality than WC or WHtR. The direct association between WHtR and overall mortality was consistent within each of the three intervention arms of the trial. Conclusions: Our study adds further support to a stronger association of abdominal obesity than BMI with total mortality among elderly subjects at high risk of cardiovascular disease. We did not find evidence to support that the PREDIMED intervention was able to counterbalance the harmful effects of increased adiposity on total mortality. Trial Registration: Controlled-Trials.com ISRCTN35739639 Citation: Martı ´nez-Gonza ´lez MA, Garcı ´a-Arellano A, Toledo E, Bes-Rastrollo M, Bullo ´M, et al. (2014) Obesity Indexes and Total Mortality among Elderly Subjects at High Cardiovascular Risk: The PREDIMED Study. PLoS ONE 9(7): e103246. doi:10.1371/journal.pone.0103246 Editor: Jose Vina, University of Valencia, Spain Received February 5, 2014; Accepted June 27, 2014; Published July 29, 2014 Copyright: ß2014 Martı ´nez-Gonza ´lez et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited. PLOS ONE | www.plosone.org 1 July 2014 | Volume 9 | Issue 7 | e103246 Funding: Supported by the official funding agency for biomedical research of the Spanish government, Instituto de Salud Carlos III (ISCIII), through grants provided to research networks specifically developed for the trial (RTIC G03/140, to Dr. Estruch; RTIC RD 06/0045, to Dr. Martı ´nez-Gonza ´lez and through Centro de Investigacio ´n Biome ´dica en Red de Fisiopatologı ´a de la Obesidad y Nutricio ´n [CIBERobn]), and by grants from Centro Nacional de Investigaciones Cardiovasculares (CNIC 06/2007), Fondo de Investigacio ´n Sanitaria–Fondo Europeo de Desarrollo Regional (PI04-2239, PI 05/2584, CP06/00100, PI07/0240, PI07/1138, PI07/0954, PI 07/0473, PI10/01407, PI10/02658, PI11/01647, and P11/02505), Ministerio de Ciencia e Innovacio ´n (AGL-2009-13906-C02 and AGL2010-22319-C03), Fundacio ´n Mapfre 2010, Consejerı ´ade Salud de la Junta de Andalucı ´a (PI0105/2007), Public Health Division of the Department of Health of the Autonomous Government of Catalonia, Generalitat Valenciana (ACOMP06109, GVACOMP2010-181, GVACOMP2011-151, CS2010-AP-111, and CS2011-AP-042), and Regional Government of Navarra (P27/2011). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript. Competing Interests: The authors have read the journal’s policy and have the following conflicts: Dr. SalasSalvado ´reports serving on the board of and receiving grant support through his institution from the International Nut and Dried Fruit Council; receiving consulting fees from Danone; and receiving grant support through his institution from Eroski and Nestle ´. Dr. Ros reports serving on the board of and receiving travel support, as well as grant support through his institution, from the California Walnut Commission; serving on the board of the Flora Foundation (Unilever); receiving lecture fees from Danone; receiving payment for the development of educational presentations from Recordati; and receiving grant support through his institution from Takeda, Daiichi Sankyo, Nutrexpa, Feiraco, Unilever, and Karo Bio. Dr. Lamuela-Raventos reports serving on the board of and receiving lecture fees from FIVIN; receiving lecture fees from Cerveceros de Espan ˜a; and receiving lecture fees and travel support from PepsiCo. Dr. Serra-Majem reports serving on the boards of the Mediterranean Diet Foundation and the Beer and Health Foundation. Dr. Pinto ´reports serving on the board of and receiving grant support through his institution from the Residual Risk Reduction Initiative (R3i) Foundation; serving on the board of Omegafort; serving on the board of and receiving payment for the development of educational presentations, as well as grant support through his institution, from Ferrer; receiving lecture fees from Danone; receiving payment for the development of educational presentations from Menarini; and receiving grant support through his institution from Kowa, Unilever, and Karo Bio. Dr. Estruch reports serving on the board of and receiving lecture fees from the Research Foundation on Wine and Nutrition (FIVIN); serving on the boards of the Beer and Health Foundation and the European Foundation for Alcohol Research (ERAB); and receiving lecture fees from Cerveceros de Espan ˜a. The other authors declare no conflict of interest relevant to the content of this manuscript. This does not alter the authors’ adherence to PLOS ONE policies on sharing data and materials. * Email: mam[email protected] Introduction The dose-response relationship between increased levels of overweight/obesity and all-cause mortality is controversial. A recent meta-analysis reported that overweight (defined as abody mass index [BMI] of 25 to 30 kg/m 2 ) was associated with significantly lower all-cause mortality in comparison with the normal weight category and that grade I obesity (BMI between 30 and 35) was associated with a non-significantly lower risk of mortality [1]. Other previous meta-analyses also suggested inverse or null associations, specially among elderly subjects [2,3]. In contrast with these reports, there is a high biological plausibility to expect a direct association between excess body weight and allcause mortality even at moderate levels of overweight. Some potential biases may have attenuated the association. They include insufficient adjustment for confounding by smoking, reverse causation bias due to disease-related weight loss [4], effect modification by age [5,6], with an attenuation of the detrimental effects of overweight or even a reversion to beneficial effects in the elderly [7]. In addition, studies using BMI as the single relevant measurement of adiposity and not other aspects of body composition such as visceral fat or fat distribution, may miss the true dose-response curve between adiposity and all-cause mortality. Alternative measures of abdominal obesity are known to be superior to BMI in showing a positive association with the risk of cardiovascular disease or diabetes [8,9,10,11]. In addition, a recent study has clearly shown significant associations between waist circumference or other indexes of abdominal obesity and allcause mortality in white subjects [12]. Furthermore, recent studies have suggested that the BMI should be refined by measuring additional indexes of fat distribution namely waist circumference (WC), waist-to-hip ratio (WHR) or waist-to height ratio (WHtR) to better identify higher-risk subjects [13,14]. WC has been often used to appraise abdominal obesity, but WC does not take differences in height into account, and subjects with a given WC will have more abdominal fat than taller subjects with the same WC [15]. The WHtR represents a further step because it also takes height into account. A systematic review found mean areas under receiving operator characteristic curves of 0.704, 0.693 and 0.671 for WHtR, WC and BMI, respectively as predictors of diabetes or cardiovascular disease [10]. We assessed the association between several anthopometric indexes (BMI, WC, WHtR, Height) and all-cause mortality in a Mediterranean cohort of elderly subjects at high cardiovascular risk included in the PREDIMED trial. The primary objective of the PREDIMED trial was to test the effect of a Mediterraneanstyle diet in primary cardiovascula prevention. Methods The protocol for the PREDIMED trial is available as supporting information at www.predimed.es. We have previously reported the design, objectives and methods of the PREDIMED trial in a specific publication [16]. The PREDIMED study is a multicenter, randomized, primary cardiovascular prevention trial conducted in Spain (www.predimed.es). The final results were reported in 2013 [17]. We randomly assigned participants to one of three diets: a Mediterranean diet (MeDiet) supplemented with extra-virgin olive oil, a MeDiet supplemented with mixed nuts, or a control group (allocated to receive advice to reduce dietary fat). All analyses were stratified by intervention group and adjusted for potential confounding. During the analysis (with the exception of spline models, see below) we used the lowest group of each anthopometric index as the reference category. The protocol of this study was approved by the Institutional Review Board of the University of Navarra and the Institutional Review Board of the Hospital Clinic (Barcelona). The Data Protection Agency is the Spanish National Agency (Agencia Espan˜ola de Proteccio´n de Datos), name of the file: PREDIMED, the responsible person for the file is Miguel A. Martı ´nez-Gonza´lez. The trial is registered at http://www.controlled-trials.com/ ISRCTN 35739639, and all clinical investigation have been conducted according to the principles expressed in the Declaration of Helsinki. Subjects who met entry requirements agreed to participate and provided written informed consent. Subjects We enrolled 7,447 participants. Forty three percent of them were men (55 to 80 years) and 57% were women (60 to 80 years). All of them were initially free of cardiovascular disease, but at high cardiovascular risk because they had at least three major cardiovascular risk factors out of six candidate risk factors or, alternatively, they were type 2 diabetics. The 6 candidate risk factors considered were: overweight/obesity, hypertension, elevated low-density lipoprotein Obesity Indexes and Mortality among Elderly Subjects at PREDIMED Study PLOS ONE | www.plosone.org 2 July 2014 | Volume 9 | Issue 7 | e103246 cholesterol, low high-density lipoprotein cholesterol, current smoking, or family history of premature coronary heart disease. The specific cut-off points for these factors and the exclusion criteria have been previously described [16]. Participants can be assumed to be on stable weights at the time of recruitment for the trial. Energy restriction was not part of the PREDIMED nutritional intervention. The recruitment took place in 11 recruiting centers between 2003 and 2009. Eighty-nine percent of candidate subjects who met entry requirements agreed to participate and provided written informed consent. Measurements Registered nurses who had been previously trained and certified to implement the PREDIMED protocol directly measured weight, height and WC of participants as previously described [16,18,19]. Height (m) and weight (kg) were measured with light clothing and no shoes with calibrated scales and a wall-mounted stadiometer, respectively; BMI was calculated as the weight in kilograms divided by the square of the height in meters; WC was measured midway between the lowest rib and the iliac crest using an anthropometric tape; in other previous studies, the protocol for WC measurement specified that WC should be measured at the minimal waist (33%), midpoint (26%) and umbilicus (27%). We selected the midpoint. However, the available evidence from a 2008 meta-analysis suggests that WC measurement protocol has no substantial influence on the association between WC, all-cause and CVD mortality, CVD and diabetes [20]; the WHtR was calculated as WC divided by height, both in centimeters. Blood pressure was measured in triplicate using a validated semiautomatic oscillometer with a 5-minute interval between each measurement and the subject in a sitting position (Omron HEM-705CP, Hoofddorp, The Netherlands). Hypertension was defined as a systolic blood pressure $140 mm Hg, a diastolic blood pressure $90 mm Hg, or the use of antihypertensive therapy. Confounders assessment Participants underwent a baseline interview that included the evaluation of cardiovascular risk factors and physician diagnoses of hypertension, diabetes and hipercolesterolemia. At the same time we gathered information about medical, socio-demographic, anthropometric, and lifestyle variables. We used the Minnesota validated physical activity questionnaire to assess leisure-time physical activity [21,22], Time spent in several activities in minutes per day was multiplied by its typical energy expenditure, expressed in metabolic equivalent tasks (METs), then summed over all activities to yield a METs-min/d score for each participant. Taking into account that the relationship between leisure-time physical activity (METs-min/d) and mortality was not linear in our data, we used a polynomial model to adjust for METs-min/d, adding a quadratic term to the multivariable model. The age range in the validation studies was 18–60 years. Dietary habits were collected through a semi-quantitative 137-item Food Frequency Questionnaire previously validated in Spain [23]. Statistical analysis We examined baseline characteristics of participants in each variable of interest according to quartiles or to predefined categories of anthropometric indexes. To compare means or percentages of each variable across quartiles (or pre-defined categories) of anthropometric indexes we used one-way ANOVA and chi-squared tests, respectively. We used Cox regression models to assess the Hazard Ratios (HR) and their 95% confidence intervals for total mortality according to quartiles or to pre-defined categories of each anthropometric index. For the multiple-adjusted model, the following potential confounders (all of them measured at baseline) were considered: age, sex, smoking, diabetes status, hypertensive status, intervention group and family history of CHD. To assess the dose-response shape between adiposity and mortality we used restricted cubic splines models with 4 degrees of freedom, and considered the point associated with the lowest mortality as the reference value for each anthropometric index. Finally, we also obtained the hazard ratios and their 95% confidence intervals for total mortality according to defined cut-off points of each anthropometric index. A p value ,0.05 was considered statistically significant. Analyses were performed using STATA version 12.1 (StataCorp, College Station, TX, USA). Results Table 1 shows baseline characteristics of the 7447 participants of PREDIMED according to quartiles (or pre-defined categories) of anthropometric indexes. We studied baseline characteristics by four variables of interest: WHtR, WC, BMI and height. Baseline risk factors (hypertension, overweight, type-2-diabetes) increased across increasing quartiles, as expected, with the exceptions of dyslipidemia and family history of premature CHD, that decreased with increasing quartiles. We also observed that adherence to the MeDiet tended to decrease across increasing quartiles of anthropometric indexes. Leisure-time physical activity followed the same inverse trend across adiposity indexes, but not across height quartiles, where a direct association was apparent. There were other differences in variables such as smoking or hypertension but they did not follow a consistent pattern across successive quartiles. The relationship between quartiles of anthropometric indexes and the risk of all-cause mortality is reflected in Table 2. We observed that the hazard ratio tended to be highest in the top quartiles for all anthropometric indexes. For WHtR, the hazard ratios slightly and non-significantly decreased in the second quartile and then increased in the upper quartiles. Thus, the ageand sex-adjusted HRs (95% confidence intervals) for categories of low, moderate, and high/very high WHtR compared to the lowest (reference) category were 0.88 (0.64–1.20), 1.01 (0.75–1.37) and 1.39 (1.02–1.88), respectively with a statistically significant linear trend (p = 0.027). In multivariable-adjusted models, the HRs were 0.98 (0.72–1.35), 1.01 (0.74– 1.38), and 1.44 (1.05–1.97) for categories of low, moderate, and high/very high WHtR, respectively, with a significant linear trend (p = 0.026). Although the p for interaction with sex was not statistically significant (p for interaction WHtR 6sex = 0.34), when we separated men and women, no significant association was observed for men. But among women, a significant linear trend remained apparent (p = 0.045 and p = 0.026 for age-adjusted, and multivariable-adjusted models, respectively). We observed that the hazard ratios for WC also increased across succesive quartiles. The ageand sex-adjusted HRs (95% confidence intervals) for quartiles of low, moderate, and high/ very high WC compared to the lowest quartile (ref.) were 1.01 (0.72–1.42), 1.17 (0.85–1.61) and 1.42 (1.03–1.96) respectively. The linear trend was statistically significant (p = 0.020). After multivariable adjustment, the HRs were 1.03 (0.73–1.46), 1.18 (0.86–1.63) and 1.36 (0.98–1.90) with a significant linear trend (p = 0.046). Although the p for interaction with sex was not statistically significant (p for interaction WC 6sex = 0.53), when Obesity Indexes and Mortality among Elderly Subjects at PREDIMED Study PLOS ONE | www.plosone.org 3 July 2014 | Volume 9 | Issue 7 | e103246 Table 1. Baseline Characteristics of Participants according to categories of waist-to-height ratio, waist circumference, body mass index and height. Waist to height Ratio (WHtR) Characteristic Lowest quartile N = 1865 Quartiles 2–3 N = 3734 Highest quartile N = 1848 Waist-to-height ratio (mean 6SD) 0.5560.03 0.6360.02 0.7160.04 Waist circumference – cm (mean 6SD) 8967.3 10166.2 11168.2 Body mass index (mean 6SD) { 26.462.5 29.962.8 33.763.3 Female sex – (%) { 51.3 53.8 71.2 Age-, – yr (mean 6SD) { 66.066.1 67.066.2 67.966.1 Smoking – (%) { Never 56.8 58.5 71.3 Former smoker 25.6 26.8 19.4 Current 17.6 14.7 9.3 Overweight (BMI $25) – (%) { 75.4 97.4 99.9 Obesity (BMI $30) – (%) { 7.2 46.0 89.2 Hypertension – (%) { 77.4 83.3 87.1 Type-2 diabetes – (%) { 45.3 48.3 52.4 Dyslipidemia – (%) 73.4 71.8 72.2 Family history of premature CHD – (%) { 26.1 21.3 20.8 Leisure-time physical activity (METS min/d) (mean 6SD) { 263.36252.6 242.26248.2 173.56193.4 Intervention group – (%) { MeDiet+EVOO 33.6 34.7 33.6 MeDiet+nuts 36.0 32.9 30.0 Control 30.4 32.4 36.4 MeDiet Adherence score (mean 6SD) { 9.062.0 8.662.0 8.362.0 Waist Circumference (WC) Characteristic Lowest quartileN = 2029 Quartiles 2–3N = 3661 Highest quartileN = 1757 Waist circumference – cm (mean 6SD) 8865.6 10163.6 11466.5 Waist-to-height ratio (mean 6SD) 0.5660.44 0.6360.38 0.7060.54 Body mass index (mean 6SD) { 27.062.9 29.963.0 33.463.5 Female sex – (%) { 75.3 53.5 45.4 Age-, – yr (mean 6SD) 67.166.0 67.066.3 66.766.3 Smoking – (%) { Never 71.8 59.3 53.3 Former smoker 16.2 26.4 30.8 Current 12.1 14.3 15.9 Overweight (BMI $25) – (%) { 78.2 97.0 99.7 Obesity (BMI $30) – (%) { 14.0 46.4 86.3 Hypertension – (%) { 79.8 82.9 85.8 Type-2 diabetes – (%) { 43.7 49.3 52.4 Dyslipidemia – (%) { 75.7 72.0 68.9 Family history of premature CHD – (%) { 25.9 21.9 19.4 Leisure-time physical activity (METS min/d) (mean 6SD) (mean 6SD) { 230.36223.8 244.16249.9 201.46231.5 Intervention group – (%) MeDiet+EVOO 34.8 34.2 33.4 MeDiet+nuts 33.9 32.9 31.9 Control 31.4 32.9 34.7 MeDiet Adherence score (mean 6SD) { 8.962.0 8.662.1 8.462.0 Body Mass Index (BMI) Characteristics , 25 kg/m 2 N = 563 . 25 & , 30 kg/m 2 N = 3402 . 30 kg/m 2 N = 3482 Obesity Indexes and Mortality among Elderly Subjects at PREDIMED Study PLOS ONE | www.plosone.org 4 July 2014 | Volume 9 | Issue 7 | e103246 Table 1. Cont. Body Mass Index (BMI) Characteristics , 25 kg/m 2 N = 563 . 25 & , 30 kg/m 2 N = 3402 . 30 kg/m 2 N = 3482 Body mass index (mean 6SD) 23.561.2 27.761.4 33.262.7 Waist-to-height ratio (mean 6SD) { 0.5460.05 0.6060.05 0.6760.06 Waist circumference – cm (mean 6SD) { 8868.1 9668.2 10669.1 Female sex – (%) { 55.1 51.8 63.4 Age-, – yr (mean 6SD) | 67.466.4 67.066.3 66.966.1 Smoking – (%) { Never 56.0 57.3 66.1 Former smoker 24.3 26.8 22.6 Current 19.7 16.0 11.3 Hypertension – (%) { 72.3 80.5 86.7 Type-2 diabetes – (%) { 64.5 46.2 48.2 Dyslipidemia – (%) 73.8 72.8 71.7 Family history of premature CHD – (%) { 30.9 21.4 22.0 Leisure-time physical activity (METS min/d) (mean 6SD) (mean 6SD) { 258.56265.5 260.66256.2 196.36211.7 Intervention group – (%) | MeDiet+EVOO 34.6 33.9 34.3 MeDiet+nuts 36.2 34.2 31.2 Control 29.1 31.9 34.5 MeDiet Adherence score (mean 6SD) { 8.862.0 8.862.0 8.462.0 Height Characteristics according to height. Lowest quartile N = 1960 Quartiles 2–3 N = 3859 Highest quartile N = 1628 Height – m (mean 6SD) 1.4960.03 1.6060.04 1.7360.04 Waist-to-height ratio (mean 6SD) { 0.6560.07 0.6360.06 0.6160.06 Waist circumference – cm (mean 6SD) { 97610.3 100610.5 10569.7 Body mass index (mean 6SD) { 30.864.1 29.963.8 29.263.4 Female sex – (%) { 97.1 59.5 5.0 Age-, – yr (mean 6SD) { 66.965.9 67.066.0 64.966.5 Smoking – (%) { Never 87.9 62.1 27.3 Former smoker 7.5 24.1 46.8 Current 4.6 13.8 25.9 Overweight (BMI $25) – (%) { 95.1 91.9 90.9 Obesity (BMI $30) – (%) { 55.2 46.2 39.1 Hypertension – (%) { 86.2 82.6 78.9 Type-2 diabetes – (%) { 46.5 48.6 50.8 Dyslipidemia – (%) | 76.4 72.3 67.2 Family history of premature CHD – (%) | 24.4 22.4 19.9 Leisure-time physical activity (METS min/d) (mean 6SD) (mean 6SD) { 176.96171.6 231.76233.9 291.66299.4 Intervention group – (%) | MeDiet+EVOO 33.8 34.3 34.3 MeDiet+nuts 30.3 33.3 32.4 Control 34.3 35.3 30.4 MeDiet Adherence score (mean 6SD) { 8.462.0 8.662.0 8.862.1 { :p,0.001. | :p,0.05. doi:10.1371/journal.pone.0103246.t001 Obesity Indexes and Mortality among Elderly Subjects at PREDIMED Study PLOS ONE | www.plosone.org 5 July 2014 | Volume 9 | Issue 7 | e103246 Table 2. Hazard Ratios (95% confidence intervals) for total mortality according to quartiles of the waist-to-height ratio, waist circumference, height and categories of body mass index. Quartiles of waist-to-height ratio HR according waist-to-height-ratio 1 (lowest) 2 3 4 (highest) Limits 0.30 to 0.59 0.59 to 0.63 0.63 to 0.67 0.67 to 1.00 Number of deaths 83 75 89 101 Person-years 8188 8132 8059 7642 P for trend Age-, sex-adjusted HR 1 (ref.) 0.88 (0.64–1.20) 1.01 (0.75–1.37) 1.39 (1.02–1.88) 0.027 Multivariable adjusted 1 1 (ref.) 0.98 (0.72–1.35) 1.01 (0.74–1.38) 1.44 (1.05–1.97) 0.026 Men Number of deaths 60 51 54 49 Person-years 3861 4022 3498 2185 Age-adjusted HR 1 (ref.) 0.80 (0.55–1.16) 0.85 (0.58–1.22) 1.28 (0.87–1.88) 0.344 Multivariable adjusted 1 1 (ref.) 0.85 (0.58–1.24) 0.76 (0.52–1.12) 1.31 (0.87–1.97) 0.430 Women Number of deaths 23 24 35 52 Person-years 4327 4110 4561 5457 Age-adjusted HR 1 (ref.) 1.07 (0.59–1.94) 1.42 (0.81–2.51) 1.65 (0.94–2.88) 0.045 Multivariable adjusted 1 1 (ref.) 1.14 (0.63–2.05) 1.38 (0.78–2.45) 1.78 (1.02–3.11) 0.026 Quartiles of waist circumference HR according waist circumference 1 (lowest) 2 3 4 (highest) Limits 50 to 94 94.3 to 100 100.6 to 107 107.5 to 177 Number of deaths 75 70 97 106 Person-years 9107 7476 8168 7271 P for trend Age-, sex-adjusted HR 1 (ref.) 1.01 (0.72–1.42) 1.17 (0.85–1.61) 1.42 (1.03–1.96) 0.020 Multivariable adjusted 1 1 (ref.) 1.03 (0.73–1.46) 1.18 (0.86–1.63) 1.36 (0.98–1.90) 0.046 Men Number of deaths 36 38 66 74 Person-years 2185 3158 4179 4045 Age-adjusted HR 1 (ref.) 0.73 (0.46–1.15) 0.92 (0.61–1.38) 1.07 (0.72–1.61) 0.323 Multivariable adjusted 1 1 (ref.) 0.74 (0.47–1.19) 0.88 (0.58–1.34) 1.04 (0.67–1.61) 0.461 Women Number of deaths 39 32 31 32 Person-years 6921 4318 3989 3226 Age-adjusted HR 1 (ref.) 1.42 (0.87–2.30) 1.43 (0.86–2.40) 1.88 (1.13–3.15) 0.018 Multivariable adjusted 1 1 (ref.) 1.42 (0.87–2.30) 1.40 (0.84–2.34) 2.02 (1.21–3.38) 0.012 Quartiles of Body mass index HR according body mass index 1 (lowest) 2 3 4 (highest) Limits 17.78 to 27.22 27.23 to 29.73 29.74 to 32.45 32.46 to 52.97 Number of deaths 98 78 82 90 Person-years 9033 8002 8083 7902 P for trend Age-, sex-adjusted HR 1 (ref.) 0.80 (0.60–1.08) 0.92 (0.68–1.24) 1.22 (0.91–1.65) 0.159 Multivariable adjusted 1 1 (ref.) 0.80 (0.59–1.08) 0.84 (0.63–1.14) 1.14 (0.85–1.53) 0.371 Men Number of deaths 67 55 48 44 Person-years 3756 3950 3567 2293 Age-, sex-adjusted HR 1 (ref.) 0.80 (0.56–1.14) 0.83 (0.57–1.20) 1.24 (0.85–1.83) 0.451 Multivariable adjusted 1 1 (ref.) 0.72 (0.50–1.04) 0.67 (0.46–0.99) 1.05 (0.70–1.56) 0.858 Women Number of deaths 31 23 34 46 Obesity Indexes and Mortality among Elderly Subjects at PREDIMED Study PLOS ONE | www.plosone.org 6 July 2014 | Volume 9 | Issue 7 | e103246 we separated men and women, the age-adjusted HRs were not significant among men. Among women the age -adjusted HRs were 1.42 (0.87–2.30), 1.43 (0.86–2.40) and 1.88 (1.13–3.15), respectively. After multivariable-adjustment, these HRs were 1.42 (0.87–2.30), 1.40 (0.84–2.34) and 2.02 (1.21–3.38). The linear trend among women was significant in both age-adjusted and multivariable-adjusted models (p = 0.018 and 0.012 respectively). For BMI the HRs increased across successive quartiles although the results were not significant. The ageand sex-adjusted HRs (95% confidence intervals) for quartiles of low, moderate, and high/very high BMI compared to the lowest quartile (ref.) were 0.80 (0.60–1.08), 0.92 (0.68–1.24) and 1.22 (0.91– 1.65). After multivariable adjustment, these HRs were 0.80 (0.59–1.08), 0.84 (0.63–1.14) and 1.14 (0.85–1.53). The linear trend was not significant in any of both cases. The interaction with sex was not significant (p = 0.77). When we separated men and women, no significant association was observed for either category although among women the association was stronger. When we classified the sample according to quartiles of height we also observed increased hazard ratios across successive quartiles. The p for interaction with sex was not significant (p = 0.35). When we studied each sex separately we observed that the HRs increased in both cases across successive quartiles. However among men this increase was higher, although in both cases, after multivariable adjustment, we observed no significant association between height and total mortality. We conducted similar analyses using a priori defined cut-off points instead of quartiles. We observed in Table 3 that hazard ratios tended to be highest in the upper categories of the anthropometric indexes, in the same way as we observed in Table 2. We observed also that the hazard ratio decreased in the second category and then it increased in the upper categories. When we analyzed separately each anthropometric index, we found this same trend. In the case of WHtR (cut-off points: 0.60, 0.65, 0.70) we observed that the hazard ratio monotonically increased. The age and sex-adjusted Hazard ratios (95% confidence intervals) for each category were 1.01 (0.77–1.32), 1.28 (0.96–1.71), and 1.58 (1.11–2.26), for WHtR categories of 0.60–0.65, 0.65–0.70 and .0.70 with respect to ,0.65, with a significant direct linear trend (p = 0.009). After multivariable adjustment, these HRs were 1.02 (0.78–1.34), 1.30 (0.97–1.75) and 1.55 (1.06–2.26), respectively with a significant linear trend (p = 0.013). Although the p for interaction with sex was not significant (p = 0.45), when we separated men and women, we observed no significant associations among men after multivariable adjustment. However among women, the multivariableadjusted HRs for successive categories were 1.34 (0.81–2.22), 1.82 (1.09–3.06) and 1.84 (1.03–3.29). The linear trend among women was significant in both age-adjusted and multivariable-adjusted Table 2. Cont. Quartiles of Body mass index HR according body mass index 1 (lowest) 2 3 4 (highest) Person-years 4277 4052 4516 5609 Age-adjusted HR 1 (ref.) 0.80 (0.47–1.36) 1.06 (0.64–1.73) 1.22 (0.77–1.95) 0.245 Multivariable adjusted 1 1 (ref.) 0.79 (0.45–1.40) 1.04 (0.64–1.69) 1.22 (0.74–2.00) 0.264 Quartiles of height HR according height 1 (lowest) 2 3 4 (highest) Limits (cm) 135 to 153 154 to 159 160 to 167 168 to 191 Number of deaths 70 64 109 105 Person-years 8518 7746 8947 6810 P for trend Age-, sex-adjusted HR 1 (ref.) 0.99 (0.70–1.41) 1.12 (0.75–1.65) 1.32 (0.86–2.04) 0.096 Multivariable adjusted 1 1 (ref.) 0.99 (0.70–1.41) 1.05 (0.71–1.54) 1.18 (0.76–1.83) 0.380 Men Limits (cm) 142 to 163 164 to 167 168 to 172 173 to 191 Number of deaths 61 49 56 48 Person-years 3883 3185 3711 2787 Age-adjusted HR 1 (ref.) 1.14 (0.79–1.66) 1.16 (0.81–1.66) 1.52 (1.04–2.22) 0.046 Multivariable adjusted 1 1 (ref.) 1.06 (0.73–1.53) 1.03 (0.72–1.48) 1.35 (0.91–2.02) 0.213 Women Limits (cm) 135 to 150 151 to 154 155 to 159 160 to 178 Number of deaths 39 33 38 24 Person-years 4767 4797 5213 3678 Age-adjusted HR 1 (ref.) 0.95 (0.60–1.52) 1.13 (0.71–1.79) 1.20 (0.71–2.06) 0.430 Multivariable adjusted 1 1 (ref.) 0.83 (0.52–1.34) 1.02 (0.64–1.65) 1.03 (0.59–1.81) 0.790 The PREDIMED study 2003–2010. 1 Adjusted for age, (sex, when pertinent) smoking, diabetes status, hypertensive status, intervention group and family history of CHD. All estimates are stratified for study center. The interaction term between age and each antropometric index was not statistically significant in all analyses: p = 0.34 (for waist-toheight ratio); p = 0.53 (for waist circumference); p = 0.35 (for height) and p = 0.77 (for body mass index). doi:10.1371/journal.pone.0103246.t002 Obesity Indexes and Mortality among Elderly Subjects at PREDIMED Study PLOS ONE | www.plosone.org 7 July 2014 | Volume 9 | Issue 7 | e103246 Table 3. Hazard Ratios (95% confidence intervals) for total mortality according to categories (a priori defined cut-off points) of the waist-to-height ratio and waist circumference. Waist-to-height ratio Limits ,= 0.60 .0.60 to 0.65 .0.65 to 0.70 .0.70 Number of deaths 104 106 87 51 Person-years 10841 10485 6921 3774 P for trend Age-, sex-adjusted HR 1 (ref.) 1.01 (0.77–1.32) 1.28 (0.96–1.71) 1.58 (1.11–2.26) 0.009 Multivariable adjusted 1 1 (ref.) 1.02 (0.78–1.34) 1.30 (0.97–1.75) 1.55 (1.06–2.26) 0.013 Men Number of deaths 75 69 46 24 Person-years 5230 4887 2613 837 Age-adjusted HR 1 (ref.) 0.91 (0.66–1.25) 1.08 (0.74–1.56) 1.62 (1.01–2.62) 0.148 Multivariable adjusted 1 1 (ref.) 0.82 (0.59–1.13) 1.01 (0.69–1.49) 1.59 (0.94–2.67) 0.238 Women Number of deaths 29 37 41 27 Person-years 5611 5598 4308 2937 Age-adjusted HR 1 (ref.) 1.29 (0.78–2.13) 1.70 (1.02–2.83) 1.67 (0.94–2.97) 0.036 Multivariable adjusted 1 1 (ref.) 1.34 (0.81–2.22) 1.82 (1.09–3.06) 1.84 (1.03–3.29) 0.016 Waist circumference Limits (cm) , = 100 . 100 to 105 . 105 to 110 . 110 Number of deaths 119 72 57 100 Person-years 13000 6632 5651 6737 P for trend Age-adjusted HR 1(ref.) 1.14 (0.85–1.53) 1.04 (0.76–1.43) 1.69 (1.29–2.23) 0.002 Multivariable adjusted 1 1(ref.) 1.18 (0.88–1.59) 1.02 (0.74–1.41) 1.57 (1.19–2.08) 0.008 Men Limits , = 100 . 100 to 105 . 105 to 110 . 110 Number of deaths 68 45 36 65 Person-years 4721 3108 2642 3096 Age-adjusted HR 1 (ref.) 1.01 (0.69–1.46) 0.89 (0.59–1.33) 1.43 (1.01–2.02) 0.093 Multivariable adjusted 1 1 (ref.) 0.98 (0.67–1.42) 0.87 (0.58–1.31) 1.32 (0.92–1.90) 0.204 Women Limits , =95 . 95 to 100 . 100 to 105 . 105 Number of deaths 40 28 18 48 Person-years 6982 3558 3285 4630 Age-adjusted HR 1 (ref.) 1.45 (0.88–2.39) 0.99 (0.55–1.80) 1.93 (1.21–3.07) 0.012 Multivariable adjusted 1 1 (ref.) 1.40 (0.85–2.31) 1.08 (0.60–1.94) 1.93 (1.21–3.07) 0.012 Height Limits (cm) , = 165 . 165 to 170 . 170 to 175 . 175 Number of deaths 175 90 52 31 Person-years 16303 8849 4854 2014 P for trend Age-, sex-adjusted HR 1 (ref.) 0.99 (0.77–1.28) 1.10 (0.80–1.52) 1.57 (1.07–2.31) 0.038 Multivariable adjusted 1 1 (ref.) 0.95 (0.73–1.23) 1.05 (0.77–1.45) 1.42 (0.97–2.08) 0.247 Men Limits (cm) , = 165 . 165 to 170 . 170 to1.75 . 175 Number of deaths 91 56 39 28 Person-years 5510 4012 2762 1282 Age-adjusted HR 1 (ref.) 0.96 (0.69–1.34) 1.11 (0.76–1.62) 1.73 (1.13–2.63) 0.046 Multivariable adjusted 1 1 (ref.) 0.89 (0.64–1.25) 1.03 (0.70–1.52) 1.57 (0.99–2.51) 0.213 Women Limits (cm) . 155 . 155 to 160 . 160 to 1.65 . 165 Number of deaths 84 34 13 3 Obesity Indexes and Mortality among Elderly Subjects at PREDIMED Study PLOS ONE | www.plosone.org 8 July 2014 | Volume 9 | Issue 7 | e103246 models (p = 0.036 and 0.016, respectively). This trend was also apparent in the spline analysis (Figure 1). In the case of WC (cut-off points: 100, 105 and 110 cm) we observed that the hazard ratios increased across successively increasing categories with the exception of one of the intermediate groups (105 to 110 cm). The age-, and sex-adjusted HRs (95% confidence intervals) were 1.14 (0.85–1.53), 1.04 (0.76–1.43) and 1.69 (1.29–2.23) with a significant linear trend (p = 0.002). After multivariable adjustment, these HRs were 1.18 (0.88–1.59), 1.02 (0.74–1.41) and 1.57 (1.19–2.08) with a significant linear trend (p = 0.008). The p for interaction with sex was not significant (p = 0.73). When we divided the sample into men and women, no significant association was observed after multivariable adjustment among men. But the association remained apparent among women (cut-off points: 95,100, 105), with multivariable-adjusted HRs of 1.40 (0.85–2.31), 1.08 (0.60–1.94) and 1.93 (1.21–3.07), with a significant linear trend (p = 0.012). The risk of total mortality across the different categories of BMI (,25, 25–30, 30–35, .35) is shown in the Table 4. We observed that the Hazard ratios decreased in the second group and then they increased in the upper categories. This trend was also apparent in the spline analysis (Figure 2). The p for interaction with sex was not significant (p = 0.43). The adjusted HRs obtained were not significant in men nor women when we stratified the sample by sex. We repeated these analyses after excluding deaths in the first follow-up year and also after excluding ever smokers. In Table 3. Cont. Height Person-years 10793 4837 2092 732 Age-adjusted HR 1 (ref.) 1.07 (0.71–1.60) 1.11 (0.61–2.01) 0.97 (0.29–3.22) 0.430 Multivariable adjusted 1 1 (ref.) 1.06 (0.69–1.61) 0.97 (0.51–1.87) 1.04 (0.33–3.26) 0.790 The PREDIMED study 2003–2010. 1 Adjusted for age, smoking, diabetes status, hypertensive status, intervention group and family history of CHD. All estimates are stratified for study center. The interaction term was not statistically significant in case of waist-to-height ratio (p = 0.45). When the analyses were restricted to deaths occurring after .= 2 yr follow-up, the only noticeable change was for the category .0.65 to 0.70 which exhibited a higher and significant elevated risk of death in the overall sample (multivariable-adjusted HR = 1.43, 95% CI: 1.03–1.97 in the total sample; 1.14 (0.77–2.42) in men, and 1.88 (1.09–3.24) in women). doi:10.1371/journal.pone.0103246.t003 Figure 1. Relative risk of death according to Waist-to-height ratio. Restricted cubic spline model adjusted for age, smoking, diabetes status, hypertensive status, intervention group and family history of CHD. doi:10.1371/journal.pone.0103246.g001 Obesity Indexes and Mortality among Elderly Subjects at PREDIMED Study PLOS ONE | www.plosone.org 9 July 2014 | Volume 9 | Issue 7 | e103246