Binge eating behaviour (BE) is characterized by: (i) eating, in a discrete period of time, an amount of food that is definitely larger than what most people would eat in a similar period of time under similar circumstances; and (ii) a sense of lack of control over eating( 1 ). BE is the primary diagnostic criterion of several disorders, i.e. binge eating disorder, bulimia nervosa, some subtypes of anorexia nervosa, and eating disorder not otherwise specified( 1 ).
BE is associated with excess energy intake, contributing to obesity and its complications( Reference Raymond, Bartholome and Lee 2 ). Diet therapy is the first-line treatment for obesity but about three-quarters of obese people under dietary treatment regain weight within one year( Reference Wing and Phelan 3 ). BE is associated with weight regain( Reference Elfhag and Rossner 4 – Reference Keranen, Savolainen and Reponen 6 ), failed weight loss( Reference de Zwaan, Nutzinger and Schoenbeck 7 – Reference Moroshko, Brennan and O’Brien 9 ) and anxiety and depression( Reference Bulik, Sullivan and Kendler 10 , Reference Pagoto, Bodenlos and Kantor 11 ). Early identification of BE is thus important in persons seeking to lose or maintain weight( Reference Andrade, Coutinho and Silva 12 ).
The prevalence of BE among individuals participating in weight-loss programmes has been reported to vary from 23 to 55 %( Reference Linde, Jeffery and Levy 13 ), but recent studies report lower estimates (17–19 %)( Reference Pacanowski, Senso and Oriogun 14 , Reference de Man Lapidoth, Ghaderi and Norring 15 ). The high variability of these estimates is partly explained by the heterogeneity of the study populations, the generally low sample sizes and the use of different methods to detect BE. Moreover, few and conflicting data are available on the association of BE with sex, age and nutritional status as most of the available studies have been performed in young obese women( Reference Presnell, Pells and Stout 5 , Reference de Zwaan, Nutzinger and Schoenbeck 7 , Reference de Man Lapidoth, Ghaderi and Norring 15 – Reference Loh and Chan 17 ).
The present study, performed in a large sample of participants to a weight loss or maintenance programme, aimed to: (i) estimate the prevalence of BE; and (ii) test whether sex, age, BMI, educational level, smoking, physical activity and marital status are risk factors for BE.
Methods
Study design
We performed a cross-sectional study on 7524 consecutive adults who self-referred to the International Center for the Assessment of Nutritional Status (ICANS, University of Milan) from July 2008 to April 2014, in order to participate to a structured weight-loss or weight-maintenance programme. On the same day, they underwent a clinical examination, an anthropometric assessment and a structured interview by a trained dietitian, in order to obtain information about marital status, educational level, smoking and structured physical activity. The latter was investigated by asking subjects the following questions: ‘Do you practise any structured physical activity?’ and ‘How many hours per week do you spend on this activity?’ Individuals who devoted ≥2 h/week to any structured physical activity were considered as active. Moreover, all subjects filled in four questionnaires to evaluate the presence of BE and anxious and depressive symptoms. In order to be eligible for the study, the subjects had to be at least 18 years old and free of major organ disease (e.g. heart failure, cancer, neurological disease).
Anthropometric assessment
Anthropometric measurements were performed following international guidelines( Reference Lohman, Roche and Reynaldo 18 ). Body weight was measured to the nearest 100 g using a SECA 700 scale and height was measured to the nearest 0·1 cm using a SECA 217 vertical stadiometer (SECA, Hamburg, Germany). BMI was calculated as [weight (kg)]/[height (m)]2 and classified according to the WHO cut-offs (BMI=18·5–24·9 kg/m2, normal weight; BMI=25·0–29·9 kg/m2, overweight; BMI=30·0–34·9 kg/m2, obesity class 1; BMI=35·0–39·9 kg/m2, obesity class 2; and BMI≥40· kg/m2, obesity class 3)( 19 ).
Psychological assessment
Eating behaviour was evaluated using the Italian version of the Binge Eating Scale (BES)( Reference Gormally, Black and Daston 20 , Reference Di Bernardo, Barciulli and Ricca 21 ). BE was defined as a BES score ≥18( Reference Marcus, Wing and Hopkins 8 ). Anxious symptoms were evaluated using the Italian version of Spielberg et al.’s State-Trait Anxiety Inventory (STAI)( Reference Spielberg, Gorsuch and Lushene 22 ). Individuals with a STAI-X1 (state anxiety) score ≥99th percentile (raw score ≥65 for men and ≥71 for women) were excluded from the study. Anxiety was defined as a STAI-X2 (trait anxiety) score ≥95th percentile (raw score ≥56 for men and ≥62 for women)( Reference Sanavio 23 ). Depression was evaluated using the Italian Depression Questionnaire (DQ) and defined as a DQ score ≥95th percentile (raw score ≥10 for men and ≥15 for women)( Reference Sanavio 23 , Reference Balsamo and Saggino 24 ). Questionnaires were considered invalid when more of 10 % of items were missing( Reference Sanavio 23 ).
Statistical analysis
Most continuous variables had non-Gaussian distributions and all are reported as 25th, 50th and 75th percentiles. Discrete variables are reported as counts and percentages. A Poisson working regression model (PWRM) with robust 95 % confidence interval was used to estimate prevalence and prevalence ratios (PR)( Reference Lumley, Kronmal and Ma 25 , Reference Barros and Hirakata 26 ). A PWRM was used because a binomial regression model (BRM) failed to converge for some of the regressions of interest( Reference Lumley, Kronmal and Ma 25 , Reference Barros and Hirakata 26 ). Expectedly, the PR estimated by the PWRM and by the BRM were similar in all cases where both could computed( Reference Lumley, Kronmal and Ma 25 , Reference Barros and Hirakata 26 ). Uni- and multivariable PWRM were used to evaluate the associations of BE with sex, age, BMI, educational level, smoking status, physical activity and marital status. Univariable PWRM were used to model the association of anxiety and depression with sex, age and BMI. The outcome variables of all models (BE, anxiety, depression) were discrete (0=no; 1=yes). The covariates were coded as follows: (i) sex (discrete, 0=male; 1=female); (ii) age (continuous, years); (iii) BMI (continuous, kg/m2); (iv) educational level (discrete, 0=<high school; 1=≥high school); (v) smoking status (discrete, 0=no; 1=yes); (vi) physical activity (discrete, 0=no; 1=yes); and (vii) marital status (discrete, 0=unmarried; 1=married). Prevalence ratios and marginal probabilities were calculated from PWRM( Reference Williams 27 , Reference Muller and MacLehose 28 ). Univariable and multivariable fractional polynomials were used to test whether the relationships of continuous predictors with the outcomes were non-linear( Reference Royston and Sauerbrei 29 ). As there was a statistically significant but practically irrelevant improvement in model fit following an inverse transformation of BMI in the multivariable PWRM having BE as outcome, BMI was kept linear with the benefit of making the relationships more understandable to a clinical audience( Reference Hosmer, Lemeshow and Sturdivant 30 ).
Results
Study population
Out of 7524 eligible subjects, seventy-one had a STAI-XI score ≥99th percentile and were excluded from analysis as specified by the study protocol. Moreover, 523 subjects had at least one invalid questionnaire. In particular, 396 did not complete the BES questionnaire, seventy-six had missing answers in the STAI-X2 questionnaire and fifty-one had missing answers in the DQ. Thus data from 6930 subjects were available for analysis and their measurements are given in Table 1. The participants were mostly women (72 %) and aged 18–81 years (median age=46 years).
n, number of subjects; P50, 50th percentile; P25, 25th percentile; P75, 75th percentile; STAI-X1, State-Trait Anxiety Inventory (state anxiety) score; SDS, standard deviation score; STAI-X2, State-Trait Anxiety Inventory (trait anxiety) score; DQ, Italian Depression Questionnaire score; BES, Binge Eating Scale score.
Table 2 reports the anthropometric status, marital status, educational level, smoking status, physical activity, and the frequency of BE, anxiety and depression in all subjects as well as stratified by sex.
n, number of subjects; %, percentage of subjects; BES, Binge Eating Scale score; STAI-X1, State-Trait Anxiety Inventory (state anxiety) score; STAI-X2, State-Trait Anxiety Inventory (trait anxiety) score; DQ, Italian Depression Questionnaire score.
* Subjects with STAI-X1≥99th percentile were excluded from the study.
Prevalence of and risk factors for binge eating
The prevalence of BE in the pooled sample was 17 % (95 % CI 16 %, 18 %). Such prevalence was higher in women compared with men (PR=2·49, 95 % CI 2·12, 2·92, P<0·001; corresponding to probabilities of 20 % v. 8 %), decreased with increasing age (PR=0·99, 95 % CI 0·98, 0·99 for each 1-year increase, P<0·001) and increased with increasing BMI (PR=1·03, 95 % CI 1·03, 1·04 for each 1-kg/m2 increase, P<0·001; univariable PWRM).
Table 3 reports the PR obtained from the multivariable PWRM. As educational level was not associated with BE prevalence (model 1 of Table 3), it was removed from the final model (model 2 of Table 3). The PR of the remaining covariates changed very little from model 2 to model 1. In the final model (model 2 of Table 3), age (PR=0·99 for each 1-year increase), physical activity (PR=0·89) and being married (PR=0·88) were inversely associated with BE prevalence. On the contrary, being a woman (PR=2·70) and smoking (PR=1·15) were directly associated with the prevalence of BE. It should be noted that the 95 % CI of the PR associated with smoking (1·02, 1·29), physical activity (0·80, 0·99) and marital status (0·79, 0·99) were wide, suggesting that larger samples of subjects are needed to estimate these effects with acceptable precision.
Values are prevalence ratios (PR) with robust 95 % confidence intervals obtained from a multivariable Poisson working model; n, number of subjects.
*P<0·05, **P<0·001.
Figure 1 plots the prevalence of BE estimated from model 2 of Table 3 as a function of age and BMI in married and non-smoking men and women. It can be seen that the prevalence of BE is higher in females across all ages and BMI levels.
Prevalence of anxiety
The prevalence of anxiety in the pooled sample was 4 % (95 % CI 3 %, 4 %). Such prevalence was similar for women v. men (PR=1·04, 95 % CI 0·80, 1·36, P=0·75) but increased with increasing age (PR=1·03, 95 % CI 1·01, 1·04, P<0·001 for 1-year increase) and BMI (PR=1·05, 95 % CI 1·03, 1·07, P<0·001 for 1-kg/m2 increase; univariable PWRM). Expectedly, the prevalence of anxiety was higher in subjects with than in those without BE (PR=4·20, 95 % CI 3·40, 5·40, P<0·001; corresponding to probabilities of 11 % v. 3 %).
Prevalence of depression
The prevalence of depression in the pooled sample was 3 % (95 % CI 3 %, 4 %). It was lower in women compared with men (PR=0·57, 95 % CI 0·43, 0·75, P<0·001), did not change with age (PR=0·99, 95 % CI 0·99, 1·00, P=0·24 for 1-year increase) and increased with BMI (PR=1·06, 95 % CI 1·04, 1·08, P<0·001 for 1-kg/m2 increase; univariable PWRM). The prevalence of depression was higher subjects with than in those without BE (PR=5·00, 95 % CI 3·90, 6·50, P<0·001; corresponding to probabilities of 9 % v. 2 %).
Lastly, the prevalence of both anxious and depressive symptoms was higher in BE than in non-BE subjects (5 %, 95 % CI 3 %, 7 % v. 1 %, 95 % CI 0 %, 1 %, P<0·001; univariable PWRM).
Discussion
To our knowledge, the present study is the largest performed so far to estimate the prevalence of and the risk factors for BE in men and women enrolling into a weight loss or maintenance programme and the first study to systematically evaluate potential risk factors for BE besides sex, age and BMI.
The study adds to the available literature in that it was performed on a large sample of men and women with wide variation in both age (18 to 81 years) and BMI (from normal weight to class 3 obesity). However, it has several limitations. First, as it was performed in subjects enrolled at a nutritional research centre, its findings are unlikely to extend to the general population. Second, because of its cross-sectional nature, it cannot test whether changes in BE are associated with changes in body weight and a cohort study is needed to test this hypothesis. Third, we did not separate moderate from severe BE, as it is sometimes done using BES, as in our view this is more coherent with the screening psychometric properties of BES.
Our estimate of 17 % for the prevalence of BE is in agreement with recent studies( Reference Pacanowski, Senso and Oriogun 14 , Reference de Man Lapidoth, Ghaderi and Norring 15 ) but lower than that reported by other studies( Reference Linde, Jeffery and Levy 13 , Reference de Zwaan and Mitchell 31 ). The most likely explanation for this difference is that we did not restrict our analysis to young obese women as done by most studies( Reference de Zwaan, Nutzinger and Schoenbeck 7 , Reference de Zwaan and Mitchell 31 , Reference Skinner, Haines and Austin 32 ). Including males, elders and normal-weight subjects, as we deliberately did here, is expected to decrease the overall prevalence of BE. Up to now, few studies have evaluated men( Reference Presnell, Pells and Stout 5 , Reference Pacanowski, Senso and Oriogun 14 – Reference Loh and Chan 17 ) and normal-weight subjects( Reference Pacanowski, Senso and Oriogun 14 , Reference Loh and Chan 17 ). In order to allow a comparison with the available studies, performed mostly on obese women, the prevalence of BE in women in the present study was 14 % (95 % CI 12 %, 16 %) for normal-weight, 18 % (95 % CI 17 %, 20 %) for overweight, 24 % (95 % CI 22 %, 27 %) for class 1 obesity, 31 % (95 % CI 27, 35 %) for class 2 obesity and 32 % (95 % CI 26, 39 %) for class 3 obesity. On the other hand, in men we found a prevalence of BE of 1 % (95 % CI 1 %, 2 %), 5 % (95 % CI 3 %, 7 %), 10 % (95 % CI 8 %, 13 %), 12 % (95 % CI 8 %, 16 %), 17 % (95 % CI 9 %, 24 %) in normal-weight, overweight and 1, 2, 3 obesity classes, respectively. Our estimate of 17 % for the prevalence of BE is also lower than what we reported in a smaller sample of the same population (n 1472)( Reference Bertoli, Spadafranca and Bes-Rastrollo 33 ). Much of this difference stems from the fact that we defined BE as a BES score ≥18 instead of a BES score ≥17, as we did in our previous analysis. The reason for this choice is that, even though the BES≥17 cut-off point is used by most Italian studies, the BES≥18 cut-off point allows international comparisons( Reference Marcus, Wing and Hopkins 8 , Reference Dalle Grave, Calugi and Petroni 34 ). It is noteworthy how a small change in the cut-off point used to diagnose BE can make its prevalence to vary so widely.
In agreement with other studies( Reference Hudson, Hiripi and Pope 35 , Reference Preti, Girolamo and Vilagut 36 ), the prevalence of BE was higher in women than in men, decreased with increasing age and increased with increasing BMI. Although some studies did not detect between-sex differences in the prevalence of BE( Reference de Man Lapidoth, Ghaderi and Norring 15 , Reference Kelly-Weeder, Phillips and Leonard 37 ), other studies reported a greater prevalence of BE among women( Reference Presnell, Pells and Stout 5 , Reference Bertoli, Spadafranca and Bes-Rastrollo 33 , Reference Dalle Grave, Calugi and Petroni 34 ). Interestingly, there was virtually no change in the association of sex with BE after correction for confounders. In keeping with other studies( Reference de Zwaan, Nutzinger and Schoenbeck 7 , Reference de Zwaan and Mitchell 31 ), we found that the prevalence of BE decreased with increasing age. The effect of age on BE prevalence was similar at univariable and multivariable analysis, with an increase of 1 year in age associated with decrease of 1 % in BE prevalence. As expected( Reference Presnell, Pells and Stout 5 , Reference Pacanowski, Senso and Oriogun 14 ), the prevalence of BE was directly associated with BMI. After correction for multiple confounders, an increase of 1 kg/m2 in BMI was associated with an increase of 5 % in BE prevalence. Very few data are presently available on the prevalence of BE in normal-weight subjects( Reference Goldschmidt, Le Grange and Powers 38 , Reference Dingemans and van Furth 39 ). We found that BE was present also among normal-weight subjects, especially among women, although to a much lesser extent than among overweight and obese subjects.
Contrarily to other researchers( Reference Pacanowski, Senso and Oriogun 14 , Reference de Man Lapidoth, Ghaderi and Norring 15 , Reference Reagan and Hersch 40 ), we found that being married was a protective factor for BE. Other variables being equal, a married subject in our study had in fact a decrease of 12 % in BE prevalence. A direct association of BE with smoking was reported by other researchers( Reference Kelly-Weeder, Phillips and Leonard 37 , Reference Anzengruber, Klump and Thornton 41 ). In the present study, smoking was associated with a 15 % increase in the prevalence of BE. Interestingly, and contrarily to other researchers( Reference Burmeister and Carels 42 ), we found an inverse association between BE and physical activity. Being active was in fact associated with an 11 % lower prevalence of BE. In agreement with most( Reference Linde, Jeffery and Levy 13 , Reference Pacanowski, Senso and Oriogun 14 , Reference Loh and Chan 17 , Reference de Zwaan and Mitchell 31 ) but not all( Reference Presnell, Pells and Stout 5 , Reference de Zwaan, Nutzinger and Schoenbeck 7 ) studies, we found that anxiety and depression were more prevalent in subjects with than in those without BE.
In conclusion, BE is common in individuals seeking to lose or maintain weight. The prevalence of BE is higher in young obese women. However, BE is present also in men, elders and normal-weight subjects. These findings reinforce the need to evaluate the eating behaviour of individuals starting weight loss or maintenance programmes.
Acknowledgements
Acknowledgements: The authors wish to thank the ICANS research staff, and especially Maria Cazzulani, Diana Osio, Lidia Lewandowski, Valentina Giustizieri, Silvia Filippini, Anna Pirovano and Giovanna Croce, for their help during this study. Financial support: This work was supported by ICANS internal grants. Conflict of interest: None. Authorship: S.B., M.G.S. and A.B. conceived the study; S.B., V.P., V.B. and A.L. carried out the study; G.B. and A.L. performed statistical analysis; S.B., A.L. and V.P. drafted the manuscript; all authors had final approval of the submitted and published versions. Ethics of human subject participation: The study was carried out according to the Declaration of Helsinki and each participant gave written informed consent to participate.