Obesity is a growing health problem in Australian adults(1). According to the most recent national data in 2007–08, 25 % of individuals aged ≥ 18 years are obese and 37 % are overweight(2). Overweight and obesity are associated with a significant burden of comorbidity from chronic disease(3). Diet is a major modifiable determinant of obesity. Research concerning the effect of diet in obesity development and associated chronic disease has tended to focus on the role of individual nutrients and foods; however, development of methods to assess the overall quality of diet has recently received significant attention(Reference Newby and Tucker4, Reference Waijers, Feskens and Ocke5). Dietary indices are based on optimal dietary patterns and provide summary measures of overall diet quality using scoring systems determined by a priori dietary recommendations.
Few studies have investigated the longitudinal association between overall quality of diet and obesity and such information is lacking in Australia. We applied a recently developed food-based dietary index, the dietary guideline index (DGI score)(Reference McNaughton, Ball and Crawford6), to assess diet quality. The DGI score is based on the Dietary Guidelines for Australian Adults(7) and the Australia Guide to Healthy Eating(Reference Kellett, Smith and Schmerlaib8) and reflects current guidelines for optimal eating patterns. The aim of the present study was to investigate the associations between diet quality and change in BMI and waist circumference (WC) in participants of the Nambour Skin Cancer Study over 15 years of follow-up.
Study design and methods
Study population
The Nambour Skin Cancer Study involved a random, community-based sample of Australian men and women during their transition from young adulthood to older age. The design, study population, baseline data collection and follow-up methods have previously been described in detail(Reference Green, Williams and Neale9, Reference van der Pols, Heinen and Hughes10). In brief, participants were 1621 residents of the Nambour township, Queensland, who were originally randomly selected from the electoral roll (voting is compulsory in Australia) and had participated in the Nambour Skin Cancer Prevention Trial (1992–6) and were followed up until 2007. Almost all study participants (99·7 %) were of Caucasian descent.
BMI was measured at a study clinic by trained staff in 1992, 1996 and 2007; and WC was measured in 1992 and 2007. Participants with available dietary data were included in the present analyses if they had BMI data available for at least one of the three time points and WC data available for at least one of the two time points. Of the total 1621 participants (56 % women) who enrolled in the Nambour Skin Cancer Prevention Trial in 1992, the analytical cohort for this study consisted of 1231 individuals for BMI analyses and 1168 individuals for WC analyses who contributed to a total of 2767 and 1686 observations, respectively, over the 15-year follow-up period. For BMI analyses, the study population included 192 (16 %) participants with one observation, 542 (40 %) with two observations and the remaining 497 (44 %) participants had observations at all three time points. For WC, the study population included 650 (56 %) with one observation and the remaining 518 (44 %) participants had observations at two time points. Participants' age ranged from 25 to 75 years at baseline (age: men, 49·8 (sd 13·1); women, 48·4 (sd 12·0) years).
The study was approved by the Ethics Committee of the Queensland Institute of Medical Research (Brisbane, QLD, Australia) and all participants provided informed written consent.
Anthropometric examinations
Weight was measured in the clinic using calibrated electronic weighing scales. Calibration was performed each day using standardised weights. Participants, wearing light clothing without shoes or socks, were asked to stand straight and relaxed on the weighing scale until a stable reading was displayed. Weight was then recorded to the nearest 0·1 kg. Standing height was measured to the nearest 0·1 cm using a stationary stadiometer. Participants were measured twice, with the average of the two measurements calculated and used for analyses. BMI was calculated as weight in kg divided by height in m2 (kg/m2). WC was measured midway between the lower rib margin and the iliac crest in the horizontal plane using a flexible non-stretch fibreglass tape and was read to the nearest mm.
Dietary intake
Dietary intake was assessed using a self-administered semi-quantitative FFQ in 1992, 1996 and 2007. The FFQ was originally developed for the US Nurses' Health Study(Reference Willett, Sampson and Stampfer11, Reference Willett12) and adapted for the Australian setting and validated in the Nambour Study population(Reference McNaughton, Marks and Gaffney13–Reference Marks, Hughes and van der Pols15). Estimates of energy and alcohol showed reasonable to good correlation when compared with weighed food records (Spearman correlation coefficient 0·45 and 0·71, respectively)(Reference Marks, Hughes and van der Pols15). The FFQ used in 1992 and 1996 included 129 food items. Because newer food items have become available in recent years, the FFQ was revised for the 2007 data collection, resulting in an expanded FFQ containing 151 food items. Participants were asked to report their usual intake over the previous 6 months, of each item, with nine frequency response options ranging from ‘never’ to ‘four or more times per d’. Frequencies were converted to daily equivalents for statistical analyses. The amounts of foods were in household or common measures such as one slice, one tablespoon (15 ml) or one cup (250 ml), representing one standard serve for each food. For seasonal fruits and vegetables, the participants were asked to indicate how often these foods were eaten in a season, and a seasonality adjustment for these foods was applied in the calculation of food intakes. Nutrient intakes were calculated using Australian food composition data(16, 17).
Participants who did not indicate consumption frequencies for ≥ 10 % of the FFQ food items and those with reported energy intakes (EI) outside the recommended normal ranges (2100–14 700 kJ/d for women, 3360–16 800 kJ/d for men)(Reference Willett12) were a priori excluded from the analyses. Specifically, a total of seventy-three participants with extremely high EI were excluded (1992: 24 (2·2 %); 1996: 24 (2·1 %); 2007: 25 (3·5 %), and one participant (2007) with EI lower than the pre-defined allowable range.
Diet quality
We measured diet quality using DGI developed by McNaughton et al. (Reference McNaughton, Ball and Crawford6). The DGI was developed to reflect Dietary Guidelines for Australian Adults(7) and consisted of fifteen food items, including dietary indicators of vegetables and legumes, fruit, total cereals, meat and alternatives, total dairy, fluids, salt, saturated fat, alcoholic beverages, added sugars and ‘extra foods’. The indicators used were based on the dietary guidelines, cut-points, and food groupings by the Australian Guide to Healthy Eating, which provides age- and sex-specific recommendations for the consumption of food items in five core food groups (vegetables, fruits, cereals, meat and alternatives, and dairy) and ‘extra foods’(Reference Kellett, Smith and Schmerlaib8). Extra foods are defined as foods that are not essential to provide nutrient requirements and contain too much fat, sugar and salt(Reference Kellett, Smith and Schmerlaib8). For the calculation of all DGI components, the original method by McNaughton et al. (Reference McNaughton, Ball and Crawford6) was exactly replicated with two exceptions. Two components, salt use and fluid intake, were excluded from the calculation of the DGI, and the intake of fruits and vegetables was based on FFQ data rather than the short dietary questions, since data on these items were lacking in the Nambour FFQ in almost all the examination years.
Each component of the DGI was scored from 0 to 10, where 10 indicated that a participant was meeting the recommendation or had an optimal intake. For example, two servings per d of fruit (recommended amount) scored 10 points, one serving per d scored 5 points and no fruit consumption scored 0 points. The total DGI score was the sum of thirteen items so that the diet score had a possible range of 0–130, with a higher score reflecting increased compliance with the dietary guidelines.
Covariates
Details of sociodemographic and lifestyle characteristics were collected through self-completed questionnaires in 1992, 1996 and 2007. Participants were considered to have a medical condition if they answered ‘yes’ to the question: ‘Have you ever been told by a doctor/nurse that you have: glaucoma, gallstones, high cholesterol, high triglycerides, diabetes/high blood sugar, high blood pressure/hypertension, angina, heart attack, stroke, cancer?’. Physical activity was categorised based on self-reported engagement in walking (if no, categorised as sedentary; if yes, categorised as low physical activity level) or moderate or vigorous exercise in the past 2 weeks. Smoking status was ascertained based on the timing and number of pack-years smoked, calculated from the frequency and duration of cigarettes smoked, for each period preceding the assessment year of anthropometric measures and categorised as lifelong non-smoker, ex-smoker and current smoker (1–7 pack-years, >7 pack-years).
Alcohol consumption (g/d) was calculated from the FFQ data, and categorised based on the national alcohol guidelines(18): none, moderate ( ≤ 40 g/d men, ≤ 20 g/d women) and heavy (>40 g/d men, >20 g/d women). We also considered frequency of alcohol consumption based on the number of alcoholic drinks per week. Estimated EI (kJ/d) was used as a continuous variable.
Age (continuous), education (grade 12 or less, technical/college diploma, trade/apprenticeship, bachelor or higher, other) and occupation categorised as professionals (e.g. manager or teacher), para-professionals (e.g. technicians) and non-professionals (e.g. trade persons or clerks)(19) were assessed using the data collected at baseline in 1992. For women, self-reported current use of hormone replacement therapy (yes/no) and parity (nulliparous, one child, 2/3/4/5–7 children) were considered in the analyses.
Statistical analysis
The longitudinal associations between the quartiles of DGI score and change in BMI or WC were assessed using linear regression applying generalised estimating equations (GEE) approach(Reference Zeger, Liang and Albert20). BMI and WC were used as continuous outcome variables. Change in each of the anthropometric measures (BMI and WC) per year in quartiles of the DGI score was calculated by including an interaction term between the variable and time (year of observation as a continuous variable). A P value for the overall association between explanatory variable and change in anthropometric measures was derived from the likelihood ratio test for the interaction of DGI score by time. A P value for subgroup comparisons within quartiles of the DGI score was derived from a Wald-test based on parameter estimates and standard error from the GEE model.
To assess the longitudinal likelihood of change from normal weight to overweight in quartiles of the DGI score, logistic regression applying GEE approach was carried out in a subsample of the study population with at least two body measures. In these analyses, dichotomous outcomes indicating change from normal weight to overweight v. no change were considered. Overweight was defined as BMI ≥ 25 (kg/m2), or WC ≥ 94 (cm) in men and WC ≥ 80 (cm) in women.
The analyses were stratified by sex because the pattern of associations was expected to be different for men and women. The covariates were included as time-dependent (changing over time) or time-independent (constant over time) variables. For all time-dependent variables (medical condition, physical activity, alcohol consumption, smoking status, EI, and hormone replacement therapy use in women), data from 1992, 1996 and 2007 were used for BMI analyses and data from 1992 and 2007 for WC analyses. Age, education, occupation, parity in 2004 and weight loss attempts in the past 10 years were considered as time-independent variables.
Multivariable regression models were adjusted for confounders which were found to be significantly related to the selected anthropometric measures. To recognise the confounders, the covariates of interest were examined in a subsample of participants with complete data on selected covariates to make sure that change in outcomes is due to the covariates and not due to missing data in outcomes of interest. Covariates were considered as confounders if there was a ± 10 % change in the outcome variable after inclusion of the covariate in the multivariable model (P < 0·1). Estimated EI did not change the outcome variable and was not significantly correlated with the DGI score and therefore omitted from the multivariable analysis. Age (centred on the mean) and its squared value were used for adjustment in the multivariable model. BMI and WC at baseline (1992) were included in the multivariable models where appropriate and were used as categorical variables according to WHO criteria for overweight and obesity(21). We did not include weight loss attempts as confounder in the multivariable model due to the limited number of observations in the final examination year (2007), although in a subsample this variable was a confounder for the association between the DGI score and change in WC in both sexes and change in BMI in men.
In sensitivity analyses, we further adjusted the multivariable models for misreporting of EI, using the computed ratio of total EI to predicted BMR (EI/BMR) as proposed by Goldberg et al. (Reference Goldberg, Black and Jebb22) and Black(Reference Black23). Underreporting was estimated by calculating the confidence limits for the agreement between EI:BMR ratio and individual physical activity level according to recommendations for evaluating ‘habitual’ intake in individuals(Reference Black23). The 95 % CI were calculated as follows: exp ( ± 2 (S/100)), where , CVwB is the CV of the within-person variation for predictions of BMR using the Schofield equation(24) and CVtP is the coefficient of variation of the total variation in physical activity level. We considered 15 % of between-subject variation for physical activity level and 8·5 % of within-subject variation for estimated BMR which are suggested as an average substitute value by Black(Reference Black23). Individuals were classified as low energy reporters (EI:BMR < the lower 95 % CI) and non-low energy reporters (EI:BMR ≥ the upper 95 % CI). To address misreporting of EI, multivariable models were repeated after exclusion of low energy reporters or adjustment for EI:BMR ratio. However, given that the results did not materially change, we decided to present the results without this adjustment.
To test the representativeness of the study sample, the characteristics of the participants included in the analyses were compared with those excluded due to missing data on anthropometric measures or dietary intake using multiple logistic regressions applying the GEE approach. Values of P < 0·05 were considered statistically significant. All analyses were carried out using SAS statistical package version 9.1 (SAS Institute, Cary, NC, USA).
Results
The final cohort for this study consisted of 1231 individuals for BMI analyses and 1168 individuals for WC analyses who contributed to a total of 2767 and 1686 observations, respectively, over the 15-year follow-up period. Study participants who were excluded from the analysis due to missing data on body size or diet (n 390) were more likely to be physically active than their counterparts included in the analysis (P =0·02) (results not shown).
Baseline characteristics of the study participants according to the DGI score are displayed in Tables 1 and 2. Men and women who achieved the highest dietary quality tended to be older, consumed less alcohol and were less likely to have a medical condition. Men in the highest quartile of the DGI score were more likely to be non-smokers and moderately physically active. Women in the highest DGI score quartile were more likely to have a higher BMI at baseline. Education, occupation, EI and WC were not statistically different across the DGI quartile in both sexes at baseline.
* Continuous variables: linear trends were identified using regression analysis. Categorical variables: significant differences were assessed using the χ2 test.
† Current smokers with ≥ 1 pack-year; ex-smokers and non-smokers not shown.
‡ Moderate: ≤ 40 g/d men, heavy: >40 g/d men.
* Continuous variables: linear trends were identified using regression analysis. Categorical variables: significant differences were assessed using the χ2 test.
† Current smokers with ≥ 1 pack-year; ex-smokers and non-smokers not shown.
‡ Moderate: ≤ 20 g/d, heavy:>20 g/d women.
At study entry (1992), 40 % of the participants were overweight (BMI 25–30) and 16 % were obese (BMI ≥ 30). In terms of WC, 28 % were at an increased risk of metabolic complications (WC ≥ 94 cm in men, ≥ 80 cm in women) and another 28 % were at a substantially increased risk (WC ≥ 102 cm in men, ≥ 88 cm in women). Men gained on average 0·09 kg/m2 per year or 1·05 kg/m2 in total BMI, while women gained 0·16 kg/m2 per year or 2·4 kg/m2 in total BMI during the follow-up. In the same time period, there was an increase of 0·12 cm/year or 1·8 cm in total WC among men, and an increase of 0·25 cm/year or 3·8 cm in total WC among women. Compliance with Australian dietary guidelines, measured by DGI, was generally poor at each of the three examination years, although diet quality slightly improved over the assessment period (1992–7). Specifically, out of a possible total DGI score of 130, in men the mean DGI score increased from 71·9 (sd 16·9) in 1992 to 73·3 (sd 17·2) in 1996 and to 76·2 (sd 15·9) in 2007. In women, the corresponding figures were 80·6 (sd 16·4), 80·8 (sd 15·8) and 83·8 (sd 15·3) (both P-trend from GEE < 0·05).
In the present study, 3 % were classified as low energy reporters (EI/BMR ≤ 0·71) (data not shown). EI/BMR was inversely correlated with BMI (r − 0·24) and WC (r − 0·23) (both P = 0·0001), indicating higher prevalence of underreporting in heavier individuals.
The associations between the DGI score based on available data in each examination year and change in BMI and WC are shown in Tables 3 (men) and 4 (women). In men, the DGI score was inversely associated with change in BMI and WC in unadjusted models, indicating that men with a higher diet quality score gained less BMI and WC over time. The association between the DGI score and change in BMI or WC was attenuated in multivariable adjusted models, but sub-group differences were present. Compared to men in the reference category, men in the highest quartile had the lowest gain in BMI (0·05 v. 0·11 kg/m2 per year, P =0·01), and those in the third quartile had the smallest increase in WC during follow-up (0·04 v. 0·26 cm/year, P =0·04). In women, the DGI score was not associated with change in any of the anthropometric measures.
Ref., reference; WC, waist circumference.
* P value from likelihood ratio test for interaction of covariate by time.
† Values derived from the interaction between covariates and time using generalised estimating equation; data collected in 1992, 1996 and 2007 for BMI, and in 1992 and 2007 for WC.
‡ Comparison between categories: P values from Wald test based on parameter estimate and standard error from generalised estimating equations model.
§ Adjusted for baseline age, education, smoking status and frequency of alcohol consumption.
‖ Adjusted for baseline age, education, smoking status, frequency of alcohol consumption and baseline BMI (categorical).
¶ In 1996, 469 observations were omitted since WC was not measured.
** Adjusted for baseline age, education, smoking status, frequency of alcohol consumption, physical activity and baseline WC (categorical).
Ref., reference; WC, waist circumference.
* P value from likelihood ratio test for interaction of covariate by time.
† Values derived from the interaction between covariates and time using generalised estimating equation; data collected in 1992, 1996 and 2007 for BMI, and in 1992 and 2007 for WC.
‡ Comparison between categories: P values from Wald test based on parameter estimate and standard error from generalised estimating equations model.
§ Adjusted for baseline age, occupation, physical activity, smoking status, frequency of alcohol consumption, hormone replacement therapy (HRT) use and parity.
‖ Adjusted for baseline age, occupation, physical activity, smoking status, frequency of alcohol consumption, HRT use, parity and baseline BMI (categorical).
¶ In 1996, 612 observations were omitted since WC was not measured.
** Adjusted for baseline age, occupation, physical activity, smoking status, frequency of alcohol consumption, HRT use, parity, education and baseline WC (categorical).
In addition, we examined the likelihood of change from normal weight to overweight by quartile of the DGI score (Tables 5 and 6). Compared to the reference group, men categorised in higher quartiles had a lower risk of becoming overweight (BMI ≥ 25) or being at increased risk of metabolic complications (WC ≥ 94 cm) during the follow-up, but none of the OR reached statistical significance. We repeated this analysis, comparing the reference group to men with a higher diet quality score categorised as one group (quartiles 2–4 combined). Men with better diet quality were 37 % less likely to change from normal weight to overweight (BMI ≥ 25) and 40 % less likely to change from no risk to increased risk of metabolic complications (WC ≥ 94 cm) during the 15-year time period (data not shown). In women, there was no difference in likelihood to change from normal to overweight based on BMI or WC by ranking of the DGI score. Furthermore, in sensitivity analyses, we additionally adjusted for EI to account for residual confounding (data not shown). This adjustment did not change the results presented in Tables 3–6.
WC, waist circumference.
* Study participants with at least two measures.
† Adjusted for baseline age, education, smoking status and frequency of alcohol consumption.
‡ In 1996, 441 observations were omitted since WC was not measured.
§ Adjusted for baseline age, education, smoking status, frequency of alcohol consumption and physical activity.
WC, waist circumference.
* Study participants with at least two measures.
† Adjusted for baseline age, occupation, physical activity, smoking status, frequency of alcohol consumption, hormone replacement therapy (HRT) use and parity.
‡ In 1996, 564 observations were omitted since WC was not measured.
§ Adjusted for baseline age, occupation, physical activity, smoking status, frequency of alcohol consumption, HRT use, parity and education.
Discussion
To our knowledge, this is the first study to examine the association between diet quality and change in body mass and fat distribution over an extended period of time in a community-based sample of Australian men and women. On average, men and women gained 2·3 kg/m2 BMI and 3·0 cm WC over the 15-year follow-up period, with women having a substantially larger gain in both BMI and WC than men. Diet quality, measured by DGI, was generally poor at three examination years. This is similar to findings from national studies showing relatively poor compliance with dietary guidelines in the Australian population(Reference McNaughton, Ball and Crawford6, Reference Ball, Mishra and Thane25).
Dietary quality was shown to modify weight gain in men only. Specifically, men with a higher diet quality were less likely to gain BMI and WC over time and change from normal to overweight during follow-up. In women, diet quality was neither associated with change in BMI or WC nor associated with the likelihood of becoming overweight. As the DGI is an indicator for the overall quality of diet, we did not attempt to tease out the individual contributions of its components to change in anthropometric measures.
Comparison with other studies is not straightforward due to different approaches in the statistical analysis techniques and construction of diet quality indices, such as choice of components, assignment of foods to food groups, choice of cut-off values, adjustment for EI and decisions on the relative contribution of the individual components to the total score. Despite these differences in methodology, there are similarities between indices such as focus on moderating consumption of dietary fats while encouraging more complex carbohydrate in the diet.
Previous evidence from longitudinal studies on the association between diet quality indices and obesity indicators is limited and cross-sectional data do not provide consistent evidence of an association(Reference Togo, Osler and Sorensen26). Other studies we compared our results with were not adjusted for many of the confounders we examined in our analyses, except for a few(Reference Quatromoni, Pencina and Cobain27, Reference Gao, Beresford and Frank28).
In men, the inverse association between DGI quartiles and obesity outcomes was supported by findings from the Framingham Offspring Study in the USA which indicated that greater diet quality, measured by diet quality index, was associated with lower weight gain in men over a time period of 8 years(Reference Quatromoni, Pencina and Cobain27). Our findings in men also concurred with those from another prospective US study, showing an inverse association between diet quality measured by the healthy eating index and obesity outcomes(Reference Gao, Beresford and Frank28). This analysis, however, was not stratified by sex. Also, a cross-sectional analysis of the Health Professionals' Follow-up Study showed an inverse association between diet quality, assessed by healthy eating index, and BMI(Reference McCullough, Feskanich and Rimm29).
In women, the lack of association between diet quality and change in anthropometric measures is in line with a recent national cross-sectional study, using the DGI score(Reference McNaughton, Dunstan and Ball30), but it contrasts earlier cross-sectional(Reference McCullough, Feskanich and Stampfer31) and longitudinal reports(Reference Quatromoni, Pencina and Cobain27) indicating inverse associations between diet quality and obesity measures. Other previous cross-sectional studies investigating diet quality(Reference Toft, Kristoffersen and Lau32, Reference Rafferty, Anderson and McGee33) or dietary patterns(Reference Togo, Osler and Sorensen26) in relation to obesity have shown mixed results.
The lack of association in women may be related to energy underreporting in those with a higher BMI(Reference Livingstone and Black34). In sensitivity analyses, we adjusted for energy misreporting, but the results remained unchanged. However, residual confounding may occur, because the DGI reflects a range of dietary behaviours, not EI only(Reference McNaughton, Dunstan and Ball30). It may also result from reverse causality, that is, overweight individuals adopting a healthier diet to manage their weight(Reference Rennie, Johnson and Jebb35). Dieting is an issue mostly observed in female participants(Reference Saarni, Silventoinen and Rissanen36, Reference Timperio, Cameron-Smith and Burns37) and this was confirmed in the present study. Overall, 22 % of women reported that they had frequently tried to lose weight in the previous 10 years, while only 10 % of men reported to do so (data not shown).
There is also the possibility of the presence of other health-related behaviours and lifestyle variables not measured in this study which may be related to both changes in dietary intake and diet quality or obesity measures, and thus may have caused uncontrolled confounding.
Evidence from cross-sectional(Reference Togo, Osler and Sorensen26) and longitudinal(Reference Newby, Muller and Hallfrisch38–Reference Togo, Osler and Sorensen43) studies on the association between empirically derived dietary patterns and obesity outcomes has also showed mixed results. In longitudinal studies, some findings suggest that a dietary pattern characterised by reduced-fat dairy products, cereals and fruit is inversely associated with change in BMI among women and with WC in both men and women(Reference Newby, Muller and Hallfrisch38–Reference McNaughton, Mishra and Stephen40). Similarly, dietary patterns characterised by whole-grain cereals, fruits and vegetables have been reported to be predictive of less gain in BMI and weight(Reference Hu, Rimm and Stampfer41, Reference Schulz, Nothlings and Hoffmann42). In contrast, in a Danish study, no clear association was found between the food scores derived by pattern analysis and change in BMI or WC over a 5-year or 11-year follow-up(Reference Togo, Osler and Sorensen43).
The strengths of this study include direct measurement of BMI and WC at three consecutive time points, thus eliminating systematic bias due to self-reported body size. Further, a number of potential confounders were accounted for and sensitivity analyses performed to estimate the impact of underreporting. This study is one of few to prospectively examine association between diet quality and change in BMI and WC in Australian adults considering a range of confounding variables over time. This study was embedded in a longitudinal study of skin cancer; therefore, anthropometric and dietary assessment was not a primary reason for participation in the study, which may have reduced some participation bias. The dietary quality index selected for this study has been previously validated in a national sample in Australia(Reference McNaughton, Ball and Crawford6) and was reported to be inversely associated with some health indicators related to type-2 diabetes in men and women(Reference McNaughton, Dunstan and Ball30). Overall, the use of dietary quality indices has several advantages over data-driven measures of dietary patterns in assessing associations between diet and health outcomes(Reference Waijers, Feskens and Ocke5). They are based on the existing knowledge of optimal dietary intakes, provide summary measures of overall dietary quality and are easy for the public to understand and interpret(Reference Michels and Schulze44).
Possible limitations of our study warrant consideration. In longitudinal analyses, loss to follow-up is a matter of concern. However, a main advantage of the applied statistical model GEE is that all available data over the follow-up period could be used. Moreover, results from sensitivity analyses indicated that participants with missing data on WC or BMI in 2007 did not significantly differ in age- and sex-adjusted mean DGI score in 1992 from those with available anthropometric data in 2007 (data not shown). Therefore, we believe that missing data are not a major limitation in our study.
As a further measure to prevent participation bias, we have adjusted the multivariable models for factors which were different between subjects included in and excluded from the study(Reference Savitz45). Another limiting factor may be reliance on self-reported lifestyle characteristics, such as physical activity, which may be subject to recall bias(Reference Wareham46).
In conclusion, despite the increasing prevalence of obesity worldwide, little is known about the association between diet quality and change in obesity measures over extended periods of time. Our findings add longitudinal evidence to the knowledge that weight gain is a problem which is widespread in the adult population, affecting both men and women. In the present study, men with a higher DGI score were less likely to gain BMI and WC and to change from normal to overweight over time. Men with the lowest diet quality score gained more than two-fold in BMI and WC (non-significantly) as compared to their counterparts with highest ranking of the score. Our findings of no association between diet quality and change in obesity measures in women indicate that the current obesity epidemic does not discriminate in its effects across groups of adult women with different diet quality levels and this remains to be further elucidated.
Acknowledgements
We thank the Nambour Skin Cancer Study participants for their continued support. We also thank Maria Celia Hughes and Torukiri Ibiebele for their contributions to collection and preparation of the anthropometric and dietary data, and particularly Maria Celia Hughes for assistance in statistical analysis. This study was supported by National Health and Medical Research Council of Australia (data collection). S. A. was supported by a University of Queensland Graduate School Research Travel Grant, a Top-Up Scholarship from the Queensland Institute of Medical Research and a World Bank Scholarship. S. A. contributed to the statistical analysis and writing of the manuscript; J. C. v. d. P. conceived the analysis plan, technical assistance in statistical analysis and assistance in writing the manuscript; G. M. W. provided technical assistance in statistical analysis; G. C. M. contributed to technical assistance and assistance in writing the manuscript; and P. H. L. provided significant advice, technical assistance in interpretation and assistance in writing the manuscript. None of the authors had a conflict of interest.