Introduction
Eating disorders are complex and multi-determined illnesses, which are influenced by genetic as well as environmental factors (Striegel-Moore & Bulik, Reference Striegel-Moore and Bulik2007). Twin studies, which partition the variance of a given trait or disorder into additive genetic (A), common environmental (i.e. factors shared among members of a twin pair, or C), and unique environmental (i.e. factors not shared by co-twins, or E) influences, have provided significant insight into the influence of genes and environment on anorexia nervosa (AN) and bulimia nervosa (BN; see Bulik et al. Reference Bulik, Sullivan, Wade and Kendler2000, for a review). However, only one published study (Javaras et al. Reference Javaras, Laird, Reichborn-Kjennerud, Bulik, Pope and Hudson2008) has investigated the heritability of binge eating disorder (BED) specifically. These authors found that, in two separate samples (one from the USA, and one from Norway), additive genetic factors significantly contributed to BED liability. In both samples, the best-fitting model included only additive genetic and unique environmental components; common environment was not significantly associated with BED liability. Variance estimates in the US sample were a2=0.57 (95% CI 0.30–0.77) and e2=0.43. In the Norwegian sample, results were comparable (a2=0.39 and e2=0.61).
Although the Javaras et al. (Reference Javaras, Laird, Reichborn-Kjennerud, Bulik, Pope and Hudson2008) study was the first to investigate the genetic epidemiology of the DSM-IV diagnosis of BED (APA, 2000), a small number of prior investigations evaluated the heritability of specific aspects of the disorder. Reichborn-Kjennerud et al. (Reference Reichborn-Kjennerud, Bulik, Kendler, Roysamb, Maes, Tambs and Harris2003) investigated binge eating with loss of control and found that among females, 48% of the variance in binge eating was accounted for by additive genetic factors (95% CI 0.17–0.60); 3% was associated with common environment (95% CI 0.00–0.29), and 49% was associated with unique environmental factors (95% CI 0.40–0.59). In a subsequent study, Reichborn-Kjennerud and colleagues (Reference Reichborn-Kjennerud, Bulik, Tambs and Harris2004) examined binge eating occurring with loss of control and in the absence of compensatory behaviors. They found that 41% of the variance in this phenotype, which more closely approximated BED, was due to additive genetic influences (95% CI 0.31–0.50), and 59% was associated with unshared environment (95% CI 0.50–0.69).
The current study attempts to extend this previous work on the genetic epidemiology of BED first by replicating Javaras and colleagues' (Reference Javaras, Laird, Reichborn-Kjennerud, Bulik, Pope and Hudson2008) study and evaluating A, C, and E influences on the latent liability to BED, as assessed by DSM-IV criteria (APA, 2000). A second aim of the current study is to investigate genetic and environmental influences on the specific symptoms of BED, which include binge eating, loss of control, frequency and duration of binges, binge characteristics (e.g. eating until feeling uncomfortably full), distress regarding binge eating, and the absence of compensatory behaviors. Previous item-level studies of AN (Mazzeo et al. Reference Mazzeo, Mitchell, Bulik, Reichborn-Kjennerud, Kendler and Neale2009) and BN (Mazzeo et al. Reference Bulik, Sullivan, Wade and Kendlerin press) found that heritability point estimates for individual symptoms of these disorders varied in magnitude.
Item-level biometric analyses can provide valuable information about the overall disorder. For example, if the criteria for a specific diagnosis vary in their respective genetic and environmental contributions to liability, use of a sum score can provide misleading results about the overall heritability of a given diagnosis. Use of sum scores can result in a loss of information, as a composite score is rarely a ‘pure’ indicator of a latent trait (Neale et al. Reference Neale, Lubke, Aggen and Dolan2005; Van den Berg et al. Reference Van den Berg, Glas and Boomsma2007). Further, data regarding environmental and genetic influences on specific eating-disorder symptoms can be quite useful regarding assessment and diagnosis of such complex phenotypes.
Symptom-level analysis might be particularly important in the assessment of BED, as the criteria for this diagnosis remain debated. For example, as Bulik and colleagues (Reference Bulik, Brownley and Shapiro2007) and others (Latner & Clyne, Reference Latner and Clyne2008; Striegel-Moore & Franko, Reference Striegel-Moore and Franko2008) have noted, terms such as ‘binge’ and ‘loss of control’ are difficult for clinicians, patients, and researchers to operationalize. Further, as is the case for BN, few empirical data are available to support the validity of the current binge frequency and duration standards (Bulik et al. Reference Bulik, Brownley and Shapiro2007; Latner & Clyne, Reference Latner and Clyne2008; Striegel-Moore & Franko, Reference Striegel-Moore and Franko2008). Thus, the current study uses an item-factor model to investigate BED (see Fig. 1). This model estimates overall A, C, and E influences on the latent continuous BED construct as well as A, C, and E contributions to each item assessing BED.
Fig. 1. Only four items per twin are shown. Λ1.1 denotes the factor loading for the first item, for twin 1. A1.1 denotes the residual variance due to A for the first item, for twin 1. r a.1=1.0 or 0.5 and denotes the twin correlation for the residual components due to A for the first item, etc. r c.1=1.0 and denotes the twin correlation for the residual components due to C for the first item, etc.
Method
Participants
Twins in this study are from the population-based Virginia Twin Registry (VTR; Kendler & Prescott, Reference Kendler and Prescott1999), which now constitutes part of the Mid-Atlantic Twin Registry (MATR). A complete description of this sample is provided elsewhere (Mazzeo et al. Reference Mazzeo, Landt-Slof, Jones, Mitchell, Kendler, Neale, Aggen and Bulik2006, Reference Bulik, Sullivan, Wade and Kendlerin press). The current study focuses on monozygotic (MZ, n=614) and (dizygotic, DZ, n=410) same-sex female twins. Participants' mean age was 40.44 years (s.d.=8.34).
Measures
BED symptoms
Items assessing BED, based on DSM-IV criteria (APA, 2000) and adapted from the Structured Clinical Interview for DSM-IV (SCID-I; First et al. Reference First, Spitzer, Gibbon and Williams1997), were administered in a self-report questionnaire (see Table 1 for an item summary). If participants reported no lifetime experiences of binge eating, they skipped all subsequent items related to this behavior. A total of 207 women endorsed binge-eating behavior. This item was scored dichotomously; all other items were assessed using ordinal scales. The item evaluating binge frequency per month used a free-response format. These responses were trichotomized (0–1, 2–7, and ⩾8 binges per month) for model-fitting purposes. Seven binary items assessing compensatory behaviors following a binge were summed to create one item, in order to assess the criterion of absence of these behaviors.
Table 1. BED item numbers, corresponding interview questions, and scoring
* If participants indicated that they had never binged, they skipped to the next section of the questionnaire.
Statistical analyses
Biometric modeling, which utilizes twin data to decompose the variance of a given trait or traits into additive genetic (A), common environmental (C), and unique environmental (E) influences, was conducted using Mx (Neale et al. Reference Neale, Boker, Xie and Maes2003). Additive genetic effects are assumed to contribute twice as much to the covariance between MZ twins as DZ twins because, generally, MZ twins share all of their genes, while DZ twins share approximately half. Of note, this assumes an absence of genetic dominance, epistasis, or assortative mating, all of which can influence heritability estimates obtained from twin models (Plomin et al. Reference Plomin, DeFries, McClearn and McGuffin2008). For example, assortative mating may lead to an overestimation of the influence due to C, and to an underestimation of variance due to A, in the classic twin model. Further, the presence of epistasis may lead to either an over- or under-estimation of A (see Keller et al. Reference Keller, Medland, Duncan, Hatemi, Neale, Maes and Eaves2008, for a review). Shared environmental influences are assumed to be equivalent for MZ and DZ twins, while specific environmental influences are assumed to be uncorrelated for both MZ and DZ twin pairs.
An item-factor analytic modeling approach (Neale et al. Reference Neale, Aggen, Maes, Kubarych and Schmitt2006a) was used to investigate BED symptoms. This method, as described in detail elsewhere, has been applied to the analysis of AN (Mazzeo et al. Reference Mazzeo, Mitchell, Bulik, Reichborn-Kjennerud, Kendler and Neale2009) and BN symptoms (Mazzeo et al. Reference Bulik, Sullivan, Wade and Kendlerin press). In this model, each item's observed variance is decomposed into two types: that shared with the latent factor and a residual component (item-specific variance plus measurement error). Residual variance for each item was partitioned into A, C, and E influences. Confidence intervals (CIs) were estimated to examine the statistical precision of the A, C, and E parameters.
Due to the design of the survey, (which involved skip patterns related to binge eating, as noted above), not all participants provided data for each item. Selection effects resulting from the use of ‘stem’ items (in this case, binge eating) can impact covariance estimates, which are essential for fitting the factor model (Neale et al. Reference Neale, Aggen, Maes, Kubarych and Schmitt2006a). However, maximum-likelihood joint analysis of data from stem and follow-up questions from twin pairs overcomes this problem, because the covariance between the stem item and the co-twin's follow-up items is available (Neale et al. Reference Neale, Harvey, Maes, Sullivan and Kendler2006b).
Total heritability for each individual item (i) was computed as the product of the item's squared factor loading (λ) and a2 for the latent trait, added to the product of one minus the item's squared factor loading and the amount of the item's residual variance due to A. This equation, where λi is the factor loading for the ith item is as follows:
Similarly, total shared and unique environmental influences on each item were computed using this equation, respectively substituting c2 or e2 and residual variance due to C or E.
Results
BED
A total of 32 women met criteria for BED. An ACE model, using MML estimation, was fitted to the BED data. The MZ correlation (r=0.58) was greater than the DZ correlation (r=0.35), providing evidence for additive genetic effects. The largest proportion of variance was due to A (a2=0.45, 95% CI 0.00–0.92), followed by E (e2=0.42, 95% CI 0.00–1.0). Last, C contributed little to the variance in liability to BED (c2=0.13, 95% CI 0.00–0.96).
The CIs for A and C included 0, indicating that A and C individually were non-significant. However, wide CIs can be a result of sparse data. Given the current study's sample size, parameters from the full ACE model are more likely to represent the true model than either submodel (Sullivan & Eaves, Reference Sullivan and Eaves2002). The fit statistics for this model were as follows: Akaike's Information Criterion=−51.68, −2 log-likelihood (2225 df)=4398.32.
This model also provides information about the variance of, and covariance among, BED criteria, through the estimation of factor loadings and residual variances, which are also partitioned into A, C, and E. All factor loadings were high in magnitude, ranging from 0.62 to 0.93 (see Table 2).
Table 2. BED item factor loadings, residual variances, and heritability estimates (95% confidence intervals)
The majority of the variance in all items was due to E, except for item 5, which had a slightly higher proportion of variance due to A. Two items (3 and 5), assessing loss of control during binges and distress due to bingeing, had heritability estimates of 0.43. Proportions of variance due to A ranged from 0.37–0.38 for binge eating (item 1), frequency of binges (item 2), and binge characteristics (item 6). The remaining two items (4 and 7), assessing duration of binges and the absence of compensatory behaviors, had heritability point estimates of 0.33 and 0.29, respectively. In addition, item 7 had the largest proportion of variance due to C (19%), relative to the other items; however, as can be seen in Table 2, no items had a large proportion of variance accounted for by common environment.
Discussion
This study was the first to investigate genetic and environmental influences on liability to BED using an item-factor approach. Use of this methodology enables examination of this disorder in three distinct ways. First, the factor loadings provide data regarding how strongly each symptom is related to a single underlying construct (or diagnosis, in this case). Second, the biometric analysis of the liability to BED enables comparison of current results with those derived from the other extant study of BED at this level (Javaras et al. Reference Javaras, Laird, Reichborn-Kjennerud, Bulik, Pope and Hudson2008). Third, the item-level analysis of the additive genetic, common environment and specific environment influences on specific BED symptoms indicates whether these individual components of BED are differentially influenced by these etiological factors.
The factor analysis results indicated that each of the DSM-IV criteria for BED was strongly related to the underlying latent liability to BED. It is particularly noteworthy that these factor loadings were higher than those previously obtained in symptom-level analyses of AN and BN (Mazzeo et al. Reference Mazzeo, Mitchell, Bulik, Reichborn-Kjennerud, Kendler and Neale2009, Reference Bulik, Sullivan, Wade and Kendlerin press). This suggests that each of these items is strongly related to a unidimensional latent construct. Thus, although specific controversies about the BED diagnosis remain, such as difficulties with operationally defining binge eating and determining appropriate cut-offs for binge frequency and duration, each of these symptoms does appear to capture a significant degree of variance in liability to this diagnosis.
Biometric analysis indicated that overall liability to BED was moderately influenced by additive genetic and unique environmental factors. Our CIs for BED variance components were quite wide. Indeed, they were much wider for the overall disorder than for the symptoms, probably because each item was more likely to be endorsed individually than altogether, reflecting the low prevalence of BED in this sample. However, the magnitude of the variance components estimates obtained in the current study are consistent with those obtained in the only other previous study of the heritability of BED (Javaras et al. Reference Javaras, Laird, Reichborn-Kjennerud, Bulik, Pope and Hudson2008), and provide additional evidence of the role of genetic factors in the etiology of this disorder. Current findings are also consistent with a previous investigation of broadly defined BN (Bulik et al. Reference Bulik, Sullivan and Kendler1998) which found that an AE model provided the best fit to the data. However, it should be noted that although these additive genetic effects on liability to BED were significant, their absolute value was moderate. Thus, it is important that researchers continue to investigate potential gene×environment interactions and correlations that might be particularly relevant to individuals with high genetic risk for BED.
Item-level analyses suggest that the specific/unique aspects of the BED symptoms all appeared to be best accounted for by a mixture of moderately strong additive genetic factors and unique environmental influences. The results differ somewhat from the previous item-level analyses conducted on the component symptoms of BN (Mazzeo et al. Reference Bulik, Sullivan, Wade and Kendlerin press) and AN (Mazzeo et al. Reference Mazzeo, Mitchell, Bulik, Reichborn-Kjennerud, Kendler and Neale2009), which found greater variability in estimates of A, C and E among specific diagnostic criteria. For example, in the item-factor analysis of BN, it was found that vomiting was more strongly influenced by additive genetic factors than the other BN symptoms. In contrast, the undue influence of weight on self-evaluation was more strongly influenced by common environmental factors than the other BN symptoms.
Overall, current results highlighting the role of the unique environment in liability to BED symptoms are consistent with those of Wade et al. (Reference Wade, Bergin, Martin, Gillespie and Fairburn2006), which found that this component of variance was significantly associated with liability to the number of lifetime eating-disordered behaviors endorsed in the Australian Twin Registry. Subsequent studies should investigate specific experiences which have been proposed as triggers for the development of eating disorders, particularly among individuals with known genetic risk for AN, BN, or BED, including trauma and encouragement to diet (Striegel-Moore & Kearney-Cooke, Reference Striegel-Moore and Kearney-Cooke1994; Pike et al. Reference Pike, Wilfley, Hilbert, Fairburn, Dohm and Striegel-Moore2006).
Importantly, however, the unique environment variance component includes measurement error. This should be considered in the interpretation of the results, and might also relate to some of the small differences identified among specific symptoms. For example, the item assessing the various compensatory behaviors (item 7) includes a range of behaviors (participants had to report not using these behaviors to meet BED criteria). It is interesting that the lowest a2 (0.29) and the second largest e2 (0.52) values were obtained for this item. The unique environment component of variance was also relatively high for items assessing binge characteristics (item 6, e2=0.56), the duration of binge-eating behaviors (item 4, e2=0.57), and their frequency (item 2, e2=0.52). Several authors have noted that appropriate thresholds for specific BED criteria, including binge frequency and duration, are difficult to define empirically (e.g. Bulik et al. Reference Bulik, Brownley and Shapiro2007; Latner & Clyne, Reference Latner and Clyne2008; Striegel-Moore & Franko, Reference Striegel-Moore and Franko2008). This difficulty in clearly defining criteria may increase measurement error for these items, thereby contributing to the size of the unique environmental component found for several items in the current study.
Several limitations of this study should be noted. First, because the sample included only European-American female twins, results might not generalize to men, non-twins, or members of other cultural groups. Second, like most twin studies, our sample includes a greater proportion of MZ twins, and individuals with relatively high levels of education (Lykken et al. Reference Lykken, McGue and Tellegen1987). Third, like many large-scale population-based twin studies, this investigation relied on self-report questionnaire data. Despite the limitations of this method, extant research has demonstrated acceptable levels of agreement between self-report, interview, and daily monitoring data in the assessment of eating disorders (e.g. Grilo et al. Reference Grilo, Masheb and Wilson2001). Last, single-item measures and ‘stem’ questions were used to evaluate BED symptoms, and symptoms were assessed over the lifespan. These measurement issues can contribute to attenuated reliability, decrease power, and increased recall bias (e.g. Crocker & Algina, Reference Crocker and Algina1986), which are potential confounds. However, the survey was designed in this manner to maximize reliability and validity while also limiting participant burden in this large population-based study.
Despite these possible limitations, this study has several strengths, including the use of a large, population-based sample, and the use of methodology which enabled examination of the factor structure as well as symptom and diagnostic-level modeling. Diagnostic criteria for psychiatric diagnoses in general, and eating disorders in particular, continue to be debated. These methods are useful for informing assessment and diagnostic techniques, particularly in the case of complex phenotypes. In the current study, BED symptoms were strongly related to latent liability to the diagnosis, providing support for the prevailing conceptualization of BED as a distinct diagnosis.
Acknowledgements
This research was supported by the National Institutes of Health Grants MH-068520 (Mazzeo), MH-20030 (Mitchell, Neale PI), MH-001553 (Bulik PI), MH-65322 (Aggen & Neale; Neale PI), and MH-40828 (Kendler). We are very grateful to the twins for their participation.
Declaration of Interest
None.
