# Food avoidance in anorexia nervosa: associated and predicting factors

L. Di Lodovico, C. Vansteene, D. Poupon, P. Gorwood, P. Duriez, Nathalie Godart, Sébastien Guillaume, Sylvain Lambert, F. Chevallier-Latreuille, Brigitte Remy, et al.

► To cite this version:

L. Di Lodovico, C. Vansteene, D. Poupon, P. Gorwood, P. Duriez, et al. Food avoidance in anorexia nervosa: associated and predicting factors. Eating and Weight Disorders - Studies on Anorexia, Bulimia and Obesity, 2023, 28 (1), pp.24. 10.1007/s40519-023-01545-4. hal-04087489

HAL Id: hal-04087489

https://u-paris.hal.science/hal-04087489

Submitted on 11 May 2023

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
# Food avoidance in anorexia nervosa: associated and predicting factors

L. Di Lodovico<fnref n="1" />,<fnref n="2" /> · C. Vansteene<fnref n="1" />,<fnref n="2" /> · D. Poupon<fnref n="1" /> · P. Gorwood<fnref n="1" />,<fnref n="2" /> · P. Duriez<fnref n="1" />,<fnref n="2" />,<fnref n="3" /> on behalf of the FFAB

Received: 12 July 2022 / Accepted: 25 January 2023
© The Author(s) 2023

## Abstract

**Purpose** Qualitative food avoidance is a significant issue in patients with anorexia nervosa (AN) and restoring diet diversity is an important part of the treatment process. We aimed to identify clinical factors which drive food avoidance and predict its maintenance in patients with AN.

**Methods** In this multicentre longitudinal study, 130 female outpatients with AN were assessed before and after 4 months of care in clinical centres specialized in AN. We assessed levels of avoidance of 16 food items, as well as body mass index (BMI), eating disorder severity, symptoms of depression and anxiety, emotional state, daily-life functioning, and body image perception.

**Results** We found that qualitative food avoidance was associated with the clinical severity of AN, anxiety and mood dimensions, and BMI- and body image-related factors. A younger age at onset predicted the maintenance of food avoidance after 4 months of treatment. Additional exploratory analyses suggested that anxiety and negative affect caused food avoidance more than the opposite.

**Conclusion** Qualitative food avoidance can be an indicator of illness severity. During treatment, focusing on reducing anxiety and negative affect may be a way to indirectly reduce food avoidance and restore diet diversity.

**Level of evidence** Level III: Evidence obtained from cohort or case-control analytic studies.

**Keywords** Anorexia nervosa · Food avoidance · Eating disorders · Depression · Anxiety

## Introduction

Anorexia nervosa (AN) is the eating disorder (ED) with the highest mortality rate [1]. It is a complex psychiatric disorder characterized by altered eating-related behaviours and body image distortions [2].

Patients with AN typically limit or avoid food intake [3]. High-calorie foods are their main target: patients rate their desire to eat high-calorie foods significantly lower than low-calorie foods, while no difference is observed in healthy controls [4]. Patients with AN, either restrictive (AN-R) or binge-eating/purging (AN-BP) type, tend to choose low-fat

foods over high-fat foods, and to undervalue the tastiness of high-fat foods [5]. Fat intake correlates with self-reported preference for high-fat food, and both are lower in patients with AN than in healthy controls [6]. Underlying processes are involved in food avoidance, altering responses to food. For example, an eye-tracking study highlighted that patients with AN avoid maintaining attention on food cues, which potentially facilitates restrictive eating [7]. These underlying processes seem to persist in weight-restored patients. Their *explicit*, self-reported, desire to eat high-calorie foods was indeed higher than in currently underweight patients and similar to that of healthy controls, but their *implicit* desire to eat high-calorie foods (assessed through reaction times) was as low as in currently underweight patients [8].

Nutritional rehabilitation is a key element in the treatment of AN [9]. As a result of their resistance to eating a variety of foods, nutrient needs are not met in patients with AN [10, 11]. This is why, while increasing food intake is important to restore weight, increasing diversity in food selections is also essential to restore nutritional status [9]. This is especially important in the long term as diet diversity is predictive of

<footnote>P. Duriez
p.duriez@ghu-paris.fr</footnote>

<footnote n="1">¹ GHU Paris Psychiatry and Neurosciences, Clinic of Mental Illnesses and Brain Disorders, 75014 Paris, France</footnote>

<footnote n="2">² Institute of Psychiatry and Neuroscience of Paris (IPNP), INSERM U1266, University of Paris, 75014 Paris, France</footnote>

<footnote n="3">³ CMME, GHU Paris Psychiatry and Neurosciences, 100 Rue de La Santé, 75674 Paris Cedex 14, France</footnote>
weight maintenance [12]. However, the literature on this topic is scarce [9]. While food avoidance is a main issue in AN and nutritional rehabilitation a main challenge during treatment, there has been limited research on factors driving food avoidance and predicting its variability in patients being treated for AN.

In the present prospective longitudinal exploratory study, we examined clinical characteristics and self-rated avoidance of 16 food items before and after 4 months of treatment. Our hypothesis was that food avoidance correlates with clinical severity of AN and other factors. As the avoided high-calorie foods may vary from one patient to another, we used a principal component analysis (PCA) to homogenize the results, expecting a single component to represent avoidance of high-calorie foods. Then, we tested our hypothesis through three aspects: (1) direct association at baseline between food avoidance, clinical severity of AN and potential associated factors; (2) colinear evolution over time; (3) capacity of clinical severity score and other factors to predict food avoidance maintenance versus successful food reintroduction. Identifying associated factors could help tackle food avoidance more efficiently during treatment.

# Methods

## Participants

Female outpatients with AN were screened for inclusion in 13 centres specialized in ED throughout France, as described in more detail elsewhere [13, 14]. Recruitment took place from February 2015 to July 2016. All patients were assessed during a face-to-face interview with a psychiatrist (who had at least 5 years of experience in ED) and were included when fulfilling the DSM5 criteria for AN [2]. Exclusion criteria were: not being affiliated to a social security system, not being fluent in French, being illiterate, not knowing how to use a computer, or presenting with dementia or delirium. Initially, 221 outpatients were included. Twenty-one patients were excluded because mandatory clinical data were missing, and 70 were lost to follow-up (35%). A total of 130 outpatients were therefore included in the present analyses.

Patients were assessed at admission (T1) and approximately 4 months later (T2). The average time period between first and second evaluations was 132 days (SD = 97.9). To address this variability, the delay between visits was included in the analyses.

Patient care can vary from one centre to another, but it consistently includes a multidisciplinary approach involving both a psychiatrist and/or a psychologist and a nutritionist or a dietician. All patients are offered at least one evidence-based psychotherapy for ED (cognitive-behavioural therapy, interpersonal therapy, family therapy, multifamily therapy),

and psychotropic drugs are prescribed when needed (primarily antidepressants).

Participants who did not attend the follow-up visit had a centre effect ($\chi^2 = 29.257, df = 12, p = 0.004$), and were characterized by a higher initial (16.128, SD = 2.966; $F = 5.116, p = 0.025$), minimum (13.919, SD = 2.159; $F = 8.397, p = 0.004$), and maximum (21.970, SD = 5.666; $F = 4.441, p = 0.036$) BMI, and lower positive (26.05, SD = 7.986; $F = 6.534, p = 0.011$) and negative (26.04, SD = 6.641; $F = 55.982, p < 0.001$) affect. Other variables did not differ [13, 14].

## Instruments

Clinical assessment included questions about current, subjective ideal, and minimum and maximum lifetime (since puberty, if present) body mass index (BMI), age at onset of AN, educational level, working activity, and the presence of a familial history of this disorder. For educational level, working activity and familial history of ED, to simplify comparisons, we divided patients into groups, i.e. university graduates versus below, working full or half-time versus not, and having at least one relative at the first or second degree diagnosed with anorexia nervosa or bulimia nervosa versus none.

Specific questionnaires and tests were provided to every patient with an established diagnosis after obtaining their consent during the first visit, at admission. All tests were repeated during the second visit, about 4 months later.

ED symptomatology was assessed using a French version of the Eating Attitudes Test-26 (EAT) [15, 16]. Three subdivisions distinguish “dieting” (13 items), “bulimia” (6 items), and “oral control” (7 items). From six-point Likert scales (from “never” to “always”), items are scored from 0 to 3 (three out of the six possible answers are rated 0). EAT total score ranges from 0 to 78. A score above 20 indicates problematic eating behaviours and a high level of concern about dieting and body weight [15].

Food avoidance was assessed through an ad hoc questionnaire constructed with a psychologist specialized in ED patients’ eating behaviours. Patients were asked to rate their levels of avoidance of 16 food items: butter, starchy foods, fries, cheese, pastries, cold meats, ham, red meat, white meat, white fish, 0% fat dairy produce, green vegetables, tomatoes, fresh fruits (except bananas), dried fruits, and whole wheat bread. Avoidance was rated on a three-point scale from 0 to 2 (0 = “I never avoid it”, 1 = “I sometimes avoid it”, 2 = “I always avoid it”).

Depression and anxiety scores were measured with the Hospital Anxiety and Depressive Scale (HADS), a self-report instrument with seven questions devoted to depression and a further seven to anxiety [17]. This instrument provides quantitative and qualitative data as, for both depression

<figure><img src="image_2.png" /></figure>
and anxiety, a score above 8 has been validated for current depressive or anxiety disorder [18]. Because the presence of a depressive or anxious disorder was evaluated both at the beginning and at the end of the protocol, we computed the number of patients in remission from these conditions during the second visit (patients with a score above 8 at T1 and a score below 8 at T2).

Patients' emotional state was assessed using the Positive and Negative Affect Schedule (PANAS), a 10-item self-report questionnaire [19]. Each item is rated on a 5-point scale, from 1 indicating that the word does "not at all" characterize the patient, to 5 meaning it "very much" does. Both scores range from 10 to 50, with higher scores indicating higher levels of positive and negative affect, respectively.

The Work and Social Adjustment Scale (WSAS) [20] assesses the level of impairment in the ability to work, home management, to engage in social and private leisure activities, and maintain close relationships. The maximum possible score is 40, with lower scores representing better functionality.

The body image perception test was based on a diagram representing the progression of ten female silhouettes, each corresponding to a specific BMI [21]. Patients were instructed to choose the silhouette that best represented their current body. Higher scores indicate higher perceived BMI and, within the scope of this study, stronger body distortion.

## Data analysis

Statistical analyses were performed using Jamovi 1.6.23 for Windows [22] and R version 4.2.1. Significance threshold was *p* < 0.05. Normal distribution was initially checked using the Kolmogorov-Smirnov test. When variables did not have a normal distribution (*p* > 0.05), we used non-parametric tests. To limit the risk of type I errors due to multiple comparisons, we controlled the false discovery rate (FDR) using the Benjamini-Hochberg adjustment with a FDR of 5% [23, 24]. Multiple regressions were then used with variables found significant.

A principal component analysis (PCA) with varimax rotation was initially performed for dimensionality reduction. Input variables were the 16 rates of avoidance of the 16 different food items at T1. Both scree plot [25] and parallel analysis [26, 27] suggested to retain three components. Factor scores for each of the three components were computed for each patient.

To test associations between factor scores of the three components and clinical characteristics, we performed Pearson's correlations for continuous variables and Mann-Whitney *U* tests for categorical variables.

To assess the evolution of factor scores and clinical characteristics between T1 and T2, we used Wilcoxon tests for

continuous variables and Chi-squared tests for categorical variables.

We computed the categorical variable "successful food reintroduction", defined as the presence of avoidance (avoidance ≥ 1) at T1 and the absence of avoidance (avoidance = 0) at T2 for at least one high-fat food (butter, fries, cheese, pastries, cold meats). We then compared profiles of patients with or without successful food reintroduction with Mann-Whitney *U* tests for continuous variables and with Chi-squared tests for categorical variables.

Finally, we explored the causal relationship between food avoidance and other variables with cross-lagged panel models using the R package "Lavaan".

## Results

The final sample included 130 patients aged between 11 and 52 years old (mean age = 25.1, SD = 10.9), with a mean age at onset of 17.2 years old (SD = 4.9), an average illness duration of 7.9 years (SD = 9.5), and an average BMI of 15.3 (SD = 1.9). The sample included 74 patients (56.9%) diagnosed with AN-R type, and 56 (43.1%) with AN-BP type.

## Food avoidance at baseline

The three principal components obtained from the PCA explained, respectively, 20.1%, 15.3% and 14.5% of the variance. Component loadings indicate that the first principal component mostly reflects avoidance of high-calorie foods, the second one avoidance of animal-based foods, and the third one avoidance of fruits and vegetables (Table S1). For clarity, we refer to the first, second and third principal components as PC-caloric, PC-animal, and PC-vegetal, respectively. Higher levels of food avoidance translate into higher factor scores.

PC-caloric avoidance factor score at T1 positively correlated with minimum lifetime BMI (*r* = 0.308, *p* < 0.001), perceived BMI (*r* = 0.226, *p* = 0.010) and perceptual distortion (*r* = 0.232, *p* = 0.008), and with EAT total (*r* = 0.520, *p* < 0.001), HADS anxiety (*r* = 0.230, *p* = 0.008) and PANAS negative (*r* = 0.211, *p* = 0.016) scores. It negatively correlated with the difference between current and minimum lifetime BMI (*r* = -0.286, *p* = 0.002; Table 1). A linear regression indicated a significant effect of EAT total score (Wald χ² = 20.412, *p* < 0.001) and minimum lifetime BMI (Wald χ² = 4.070, *p* = 0.046; Fig. 1).

PC-animal avoidance factor score at T1 negatively correlated with subjective ideal BMI (*r* = -0.253, *p* = 0.006). It positively correlated with minimum lifetime BMI (*r* = 0.251, *p* = 0.007), perceived BMI (*r* = 0.263, *p* = 0.003), perceptual distortion (*r* = 0.275, *p* = 0.002), EAT total (*r* = 0.408, *p* < 0.001), HADS anxiety
<table><caption>**Table 1** Factors correlated to, or associated with, factor scores of food avoidance in 130 patients with anorexia nervosa at baseline</caption><thead><tr><th rowspan="2">Patients' characteristics</th><th colspan="3">PC-caloric avoidance factor score at T2</th><th colspan="3">PC-animal avoidance factor score at T2</th><th colspan="3">PC-vegetal avoidance factor score at T2</th></tr><tr><th>r</th><th>U</th><th>p</th><th>r</th><th>U</th><th>p</th><th>r</th><th>U</th><th>p</th></tr></thead><tbody><tr><td>Age</td><td>-0.106</td><td>0.229</td><td></td><td>-0.124</td><td>0.159</td><td></td><td>0.089</td><td>0.316</td><td></td></tr><tr><td>Education</td><td></td><td>1789</td><td>0.308</td><td></td><td>1633</td><td>0.077</td><td></td><td>1861</td><td>0.500</td></tr><tr><td>Working</td><td></td><td>1003</td><td>0.837</td><td></td><td>1016</td><td>0.911</td><td></td><td>961</td><td>0.610</td></tr><tr><td>Familial history of ED</td><td></td><td>1432</td><td>0.187</td><td></td><td>1536</td><td>0.440</td><td></td><td>1680</td><td>0.987</td></tr><tr><td>Subtype</td><td></td><td>2010</td><td>0.772</td><td></td><td>2064</td><td>0.972</td><td></td><td>1516</td><td>0.009</td></tr><tr><td>Age at onset</td><td>-0.105</td><td></td><td>0.234</td><td>-0.173</td><td></td><td>0.049</td><td>-0.008</td><td></td><td>0.927</td></tr><tr><td>Illness duration</td><td>-0.068</td><td></td><td>0.444</td><td>-0.053</td><td></td><td>0.546</td><td>0.105</td><td></td><td>0.233</td></tr><tr><td>Current BMI</td><td>0.008</td><td></td><td>0.925</td><td>-0.007</td><td></td><td>0.939</td><td>0.092</td><td></td><td>0.295</td></tr><tr><td>Minimum lifetime BMI</td><td>0.308</td><td></td><td><strong>&lt;0.001</strong></td><td>0.251</td><td></td><td><strong>0.007</strong></td><td>0.195</td><td></td><td>0.037</td></tr><tr><td>Maximum lifetime BMI</td><td>0.120</td><td></td><td>0.201</td><td>0.022</td><td></td><td>0.818</td><td>0.077</td><td></td><td>0.410</td></tr><tr><td>BMI max-current</td><td>0.104</td><td></td><td>0.270</td><td>0.002</td><td></td><td>0.983</td><td>0.042</td><td></td><td>0.652</td></tr><tr><td>BMI current-min</td><td>-0.286</td><td></td><td><strong>0.002</strong></td><td>-0.231</td><td></td><td><strong>0.013</strong></td><td>-0.131</td><td></td><td>0.162</td></tr><tr><td>BMI max-min</td><td>-0.037</td><td></td><td>0.696</td><td>-0.108</td><td></td><td>0.252</td><td>-0.022</td><td></td><td>0.816</td></tr><tr><td>Perceived BMI</td><td>0.226</td><td></td><td><strong>0.010</strong></td><td>0.263</td><td></td><td><strong>0.003</strong></td><td>0.028</td><td></td><td>0.753</td></tr><tr><td>Perceptual distortion</td><td>0.232</td><td></td><td><strong>0.008</strong></td><td>0.275</td><td></td><td><strong>0.002</strong></td><td>-0.047</td><td></td><td>0.594</td></tr><tr><td>Subjective ideal BMI</td><td>-0.092</td><td></td><td>0.321</td><td>-0.253</td><td></td><td><strong>0.006</strong></td><td>-0.053</td><td></td><td>0.565</td></tr><tr><td>EAT Total</td><td>0.520</td><td></td><td><strong>&lt;0.001</strong></td><td>0.408</td><td></td><td><strong>&lt;0.001</strong></td><td>0.263</td><td></td><td>0.003</td></tr><tr><td>EAT Dieting</td><td>0.557</td><td></td><td><strong>&lt;0.001</strong></td><td>0.441</td><td></td><td><strong>&lt;0.001</strong></td><td>0.230</td><td></td><td>0.008</td></tr><tr><td>EAT Bulimia</td><td>0.266</td><td></td><td><strong>0.002</strong></td><td>0.191</td><td></td><td>0.029</td><td>0.190</td><td></td><td>0.031</td></tr><tr><td>EAT Oral</td><td>0.325</td><td></td><td><strong>&lt;0.001</strong></td><td>0.260</td><td></td><td><strong>0.003</strong></td><td>0.221</td><td></td><td>0.011</td></tr><tr><td>HADS anxiety score</td><td>0.230</td><td></td><td><strong>0.008</strong></td><td>0.285</td><td></td><td><strong>0.001</strong></td><td>0.213</td><td></td><td>0.015</td></tr><tr><td>HADS anxiety syndrome</td><td></td><td>653</td><td>0.128</td><td></td><td>692</td><td>0.215</td><td></td><td>827</td><td>0.799</td></tr><tr><td>HADS depression score</td><td>0.100</td><td></td><td>0.258</td><td>0.172</td><td></td><td>0.051</td><td>0.089</td><td></td><td>0.314</td></tr><tr><td>HADS depression syndrome</td><td></td><td>1536</td><td>0.037</td><td></td><td>1579</td><td>0.061</td><td></td><td>1776</td><td>0.356</td></tr><tr><td>PANAS positive</td><td>0.005</td><td></td><td>0.953</td><td>-0.128</td><td></td><td>0.148</td><td>-0.030</td><td></td><td>0.734</td></tr><tr><td>PANAS negative</td><td>0.211</td><td></td><td><strong>0.016</strong></td><td>0.245</td><td></td><td><strong>0.005</strong></td><td>0.197</td><td></td><td>0.025</td></tr><tr><td>WSAS</td><td>0.165</td><td></td><td>0.061</td><td>0.081</td><td></td><td>0.359</td><td>0.137</td><td></td><td>0.120</td></tr></tbody></table>

BMI body mass index, EAT eating attitudes test-26, ED eating disorder, HADS Hospital Anxiety and Depressive Scale, *p* *p*-value, PANAS positive and negative affect schedule, PC principal component, *r* Pearson's *r*, U Mann-Whitney *U*, WSAS Work and Social Adjustment Scale. Bold numbers indicate significant *p*-values after implementation of the Benjamini-Hochberg procedure

(r = 0.285, *p* = 0.001) and PANAS negative (r = 0.245, *p* = 0.005) scores (Table 1). Effects of minimum lifetime BMI (Wald χ² = 11.526, *p* < 0.001) and perceptual distortion (Wald χ² = 7.037, *p* = 0.009) remained significant in the linear regression, while a tendency was observed for subjective ideal BMI (Wald χ² = 3.553, *p* = 0.062).

None of the correlations with PC-vegetal avoidance factor score at T1 was significant (Table 1).

## Evolution between T1 and T2

Factor scores, current BMI, perceived BMI, and EAT, HADS, PANAS, and WSAS scores all significantly evolved between T1 and T2 (Table S2).

## Predicting factors of food avoidance

We identified predicting factors of food avoidance through Pearson's correlations between factor scores at T2 and patients' clinical characteristics at T1 (Table 2).

Avoidance of high-calorie foods at T2 was predicted by higher EAT (r = 0.434, *p* < 0.001), HADS anxiety (r = 0.281, *p* = 0.001), PANAS negative (r = 0.276, *p* = 0.001) and WSAS (r = 0.212, *p* = 0.016) scores at T1. It was also predicted by a younger age at onset (r = -0.263, *p* = 0.003), a higher perceived BMI (r = 0.257, *p* = 0.003) and more perceptual distortion (r = 0.241, *p* = 0.006; Table 2). Results from the linear regression indicated a significant effect of EAT total score (Wald χ² = 14.458, *p* < 0.001) and age at onset (Wald χ² = 8.662, *p* = 0.004; Fig. 1).

<figure><img src="image_3.png" /></figure>
<figure><img src="image_4.png" /></figure>

**Fig. 1** Factors associated with the description (T1) and evolution (T2) of the component “high-calorie food avoidance” according to multiple regressions. High-calorie food avoidance at baseline (T1) and after 4 months of treatment (T2) correspond to PC-caloric avoidance factor scores at T1 and T2, respectively. Successful high-fat food reintroduction was defined as the presence of avoidance at T1 and the absence of avoidance at T2 for at least one high-fat food (butter, fries, cheese, pastries, or cold meats). Factors were assessed at T1. The left set of arrows therefore depicts associations between

Avoidance of animal-based foods at T2 was predicted by higher EAT ($r$=0.263, $p$=0.002), and HADS anxiety ($r$=0.211, $p$=0.016) scores, and by a lower PANAS positive score ($r$=−0.217, $p$=0.013). It was also predicted by a younger age at onset ($r$=−0.308, $p$<0.001), a lower subjective ideal BMI ($r$=−0.243, $p$=0.008), a higher perceived BMI ($r$=0.209, $p$=0.017) and more perceptual distortion ($r$=0.220, $p$=0.012; Table 2). Age at onset (Wald $\chi^2$=4.642, $p$=0.033) and EAT total score (Wald $\chi^2$=4.213, $p$=0.042) retained a significant effect in the linear regression.

No factor significantly predicted avoidance of fruits and vegetables at T2 (Table 2).

Improvement of avoidance

Factor scores covaried with different variables. Greater reduction of PC-caloric avoidance factor score was strongly associated with greater BMI increase ($r$=−0.227, $p$=0.009) and greater reductions of perceptual distortion ($r$=0.252, $p$=0.004) and EAT score ($r$=0.383, $p$<0.001). A greater reduction of PC-caloric avoidance factor score was observed in patients in remission of depression ($U$=1150, $p$=0.011). Only the effect of EAT score remained significant in the linear regression (Wald $\chi^2$=10.980, $p$=0.001). Greater reduction of PC-animal

high-calorie food avoidance and factors at T1, while the right set of arrows indicates factors at T1 which predict high-calorie food avoidance or successful high-fat food reintroduction at T2. Asterisks indicate factors which covary with high-calorie food avoidance between T1 and T2. Arrow thickness reflects *p*-values obtained in the multiple regressions. AN anorexia nervosa, BMI body mass index, EAT Eating Attitudes Test Total Score, HADS Hospital Anxiety and Depressive Scale, *p* *p*-value, PANAS positive and negative affect schedule, WSAS Work and Social Adjustment Scale

avoidance factor score strongly correlated with a greater reduction of EAT score ($r$=0.324, $p$<0.001). PC-vegetal avoidance factor score did not covary with any factors (Fig. 1; Table 3).

Patients who successfully reintroduced at least one high-fat food had, at T1, lower HADS anxiety ($U$=1243, $p$=0.007) and EAT bulimia ($U$=1222, $p$=0.005) scores. EAT oral score also improved more in patients who successfully reintroduced at least one high-fat food than in those who did not ($U$=1215, $p$=0.004; Table 4). EAT oral score (Wald $\chi^2$=7.941, $p$=0.005) and HADS anxiety score (Wald $\chi^2$=3.883, $p$=0.049) retained a significant effect in the logistic regression (Fig. 1).

Cross-lagged panel models

Because avoidance factor scores correlate with other variables at both T1 and T2, we performed exploratory analyses to further investigate the causal relationship between high-calorie food avoidance and variables of interest, namely anxiety, depression and negative affect (Fig. 2). The cross-lagged paths suggest that anxiety and negative affect cause food avoidance, and not the opposite. For depression, the relationship was not significant, suggesting that food avoidance and depression do not cause each other.
<table><caption>**Table 2** Predicting factors of food avoidance</caption><thead><tr><th rowspan="2">Patients' characteristics at T1</th><th colspan="3">PC-caloric avoidance factor score at T2</th><th colspan="3">PC-animal avoidance factor score at T2</th><th colspan="3">PC-vegetal avoidance factor score at T2</th></tr><tr><th>r</th><th>U</th><th>p</th><th>r</th><th>U</th><th>p</th><th>r</th><th>U</th><th>p</th></tr></thead><tbody><tr><td>Age</td><td>-0.088</td><td>0.321</td><td></td><td>-0.097</td><td>0.273</td><td></td><td>0.049</td><td>0.581</td><td></td></tr><tr><td>Education</td><td></td><td>1957</td><td>0.829</td><td></td><td>1924</td><td>0.708</td><td></td><td>1659</td><td>0.100</td></tr><tr><td>Working</td><td></td><td>1030</td><td>0.991</td><td></td><td>1011</td><td>0.882</td><td></td><td>911</td><td>0.383</td></tr><tr><td>Familial history of ED</td><td></td><td>1622</td><td>0.748</td><td></td><td>1527</td><td>0.412</td><td></td><td>1387</td><td>0.120</td></tr><tr><td>Subtype</td><td></td><td>1773</td><td>0.160</td><td></td><td>1740</td><td>0.119</td><td></td><td>1648</td><td>0.046</td></tr><tr><td>Age at onset</td><td>-0.263</td><td><strong>0.003</strong></td><td></td><td>-0.308</td><td><strong>&lt;0.001</strong></td><td></td><td>-0.070</td><td>0.426</td><td></td></tr><tr><td>Illness duration</td><td>0.034</td><td>0.700</td><td>0.047</td><td></td><td>0.598</td><td>0.092</td><td></td><td>0.299</td><td></td></tr><tr><td>Current BMI</td><td>0.014</td><td>0.879</td><td>-0.008</td><td></td><td>0.932</td><td>0.082</td><td></td><td>0.352</td><td></td></tr><tr><td>Minimum lifetime BMI</td><td>0.144</td><td>0.126</td><td>0.124</td><td></td><td>0.186</td><td>0.162</td><td></td><td>0.083</td><td></td></tr><tr><td>Maximum lifetime BMI</td><td>0.159</td><td>0.089</td><td>-0.050</td><td></td><td>0.596</td><td>-0.014</td><td></td><td>0.879</td><td></td></tr><tr><td>BMI max-current</td><td>0.151</td><td>0.107</td><td>-0.070</td><td></td><td>0.456</td><td>-0.052</td><td></td><td>0.579</td><td></td></tr><tr><td>BMI current-min</td><td>-0.118</td><td>0.209</td><td>-0.101</td><td></td><td>0.281</td><td>-0.103</td><td></td><td>0.274</td><td></td></tr><tr><td>BMI max-min</td><td>0.088</td><td>0.351</td><td>-0.115</td><td></td><td>0.222</td><td>-0.098</td><td></td><td>0.295</td><td></td></tr><tr><td>Perceived BMI</td><td>0.257</td><td><strong>0.003</strong></td><td>0.209</td><td></td><td><strong>0.017</strong></td><td>-0.009</td><td></td><td>0.915</td><td></td></tr><tr><td>Perceptual distortion</td><td>0.241</td><td><strong>0.006</strong></td><td>0.220</td><td></td><td><strong>0.012</strong></td><td>-0.065</td><td></td><td>0.462</td><td></td></tr><tr><td>Subjective ideal BMI</td><td>-0.193</td><td>0.036</td><td>-0.243</td><td></td><td><strong>0.008</strong></td><td>-0.005</td><td></td><td>0.955</td><td></td></tr><tr><td>EAT Total</td><td>0.434</td><td><strong>&lt;0.001</strong></td><td>0.263</td><td></td><td><strong>0.002</strong></td><td>0.091</td><td></td><td>0.303</td><td></td></tr><tr><td>EAT Dieting</td><td>0.465</td><td><strong>&lt;0.001</strong></td><td>0.282</td><td></td><td><strong>0.001</strong></td><td>0.084</td><td></td><td>0.342</td><td></td></tr><tr><td>EAT Bulimia</td><td>0.242</td><td><strong>0.006</strong></td><td>0.172</td><td></td><td>0.050</td><td>0.096</td><td></td><td>0.278</td><td></td></tr><tr><td>EAT Oral</td><td>0.250</td><td><strong>0.004</strong></td><td>0.130</td><td></td><td>0.142</td><td>0.040</td><td></td><td>0.652</td><td></td></tr><tr><td>HADS anxiety score</td><td>0.281</td><td><strong>0.001</strong></td><td>0.211</td><td></td><td><strong>0.016</strong></td><td>-0.013</td><td></td><td>0.884</td><td></td></tr><tr><td>HADS anxiety syndrome</td><td></td><td>622</td><td>0.080</td><td></td><td>856</td><td>0.965</td><td></td><td>827</td><td>0.799</td></tr><tr><td>HADS depression score</td><td>0.098</td><td>0.266</td><td>0.141</td><td></td><td>0.109</td><td>0.119</td><td></td><td>0.179</td><td></td></tr><tr><td>HADS depression syndrome</td><td></td><td>1605</td><td>0.080</td><td></td><td>1693</td><td>0.185</td><td></td><td>1594</td><td>0.072</td></tr><tr><td>PANAS positive</td><td>-0.075</td><td>0.394</td><td>-0.217</td><td></td><td><strong>0.013</strong></td><td>-0.146</td><td></td><td>0.098</td><td></td></tr><tr><td>PANAS negative</td><td>0.276</td><td><strong>0.001</strong></td><td>0.140</td><td></td><td>0.112</td><td>-0.119</td><td></td><td>0.178</td><td></td></tr><tr><td>WSAS</td><td>0.212</td><td><strong>0.016</strong></td><td>0.095</td><td></td><td>0.282</td><td>0.042</td><td></td><td>0.635</td><td></td></tr></tbody></table>

BMI body mass index, EAT eating attitudes test-26, ED eating disorder, HADS Hospital Anxiety and Depressive Scale, *p* *p*-value, PANAS positive and negative affect schedule, PC principal component, *r* Pearson's *r*, T1 at baseline, T2 after 4 months of treatment, U Mann-Whitney *U*, WSAS Work and Social Adjustment Scale. Bold numbers indicate significant *p*-values after implementation of the Benjamini-Hochberg procedure

Discussion

After distributing food avoidance in three domains (high caloric, animal-based and vegetal), we found that food avoidance, especially for high-calorie foods, was associated with the clinical severity of AN (as indicated by EAT scores) and with mood and anxiety dimensions. A younger age at onset was associated with the maintenance of food avoidance. We also observed associations between food avoidance and BMI- and body image-related factors such as minimum lifetime BMI, subjective ideal BMI and perceptual distortion.

Correlations between food avoidance and EAT scores confirm that food avoidance is associated with AN severity and validate our ad hoc questionnaire. These correlations were expected as food restriction typically results from the concerns about weight which are characteristic of AN [28].

Food restriction can be quantitative or qualitative: while patients can exclude certain foods from their diet (qualitative restriction), they can also restrict the amount of food they eat without reducing diet diversity (quantitative restriction) [29]. Qualitative restriction is the focus of the present study, as we assess the avoidance of specific food items. We observed that, in patients who reintroduced at least one high-fat food into their diet, EAT scores were lower at baseline and improved more during treatment, once again confirming that food avoidance is related to illness severity.

We found that a younger age at onset was associated with more maintenance of food avoidance. This observation concurs with the literature, as a younger age at onset is known to be associated with more severe symptomatology [30, 31]. The impact of age at onset on illness severity may stem from the aetiology of AN, which involves a complex

<figure><img src="image_5.png" /></figure>
<table><caption>**Table 3** Clinical characteristics covarying with factor scores of food avoidance between two visits of 130 patients with anorexia nervosa</caption><thead><tr><th rowspan="2">Patients' characteristics (T2-T1)</th><th colspan="3">PC-caloric avoidance factor score difference (T2-T1)</th><th colspan="3">PC-animal avoidance factor score difference (T2-T1)</th><th colspan="3">PC-vegetal avoidance factor score difference (T2-T1)</th></tr><tr><th>r</th><th>U</th><th>p</th><th>r</th><th>U</th><th>p</th><th>r</th><th>U</th><th>p</th></tr></thead><tbody><tr><td>Age</td><td>0.028</td><td></td><td>0.754</td><td>0.031</td><td></td><td>0.727</td><td>-0.034</td><td></td><td>0.700</td></tr><tr><td>Current BMI</td><td>-0.227</td><td></td><td><strong>0.009</strong></td><td>-0.071</td><td></td><td>0.419</td><td>0.172</td><td></td><td>0.050</td></tr><tr><td>Perceived BMI</td><td>0.075</td><td></td><td>0.396</td><td>0.019</td><td></td><td>0.827</td><td>0.035</td><td></td><td>0.691</td></tr><tr><td>Perceptual distortion</td><td>0.252</td><td></td><td><strong>0.004</strong></td><td>0.065</td><td></td><td>0.466</td><td>-0.100</td><td></td><td>0.257</td></tr><tr><td>EAT Total</td><td>0.383</td><td></td><td><strong>&lt;0.001</strong></td><td>0.324</td><td></td><td><strong>&lt;0.001</strong></td><td>0.125</td><td></td><td>0.157</td></tr><tr><td>EAT Dieting</td><td>0.418</td><td></td><td><strong>&lt;0.001</strong></td><td>0.351</td><td></td><td><strong>&lt;0.001</strong></td><td>0.088</td><td></td><td>0.320</td></tr><tr><td>EAT Bulimia</td><td>0.255</td><td></td><td><strong>0.003</strong></td><td>0.230</td><td></td><td><strong>0.008</strong></td><td>0.119</td><td></td><td>0.177</td></tr><tr><td>EAT Oral</td><td>0.198</td><td></td><td>0.024</td><td>0.161</td><td></td><td>0.068</td><td>0.118</td><td></td><td>0.181</td></tr><tr><td>HADS anxiety score</td><td>0.161</td><td></td><td>0.067</td><td>0.155</td><td></td><td>0.079</td><td>0.043</td><td></td><td>0.629</td></tr><tr><td>Anxiety remission</td><td></td><td>826</td><td>0.078</td><td></td><td>910</td><td>0.221</td><td></td><td>1084</td><td>0.920</td></tr><tr><td>HADS depression score</td><td>0.182</td><td></td><td>0.039</td><td>0.103</td><td></td><td>0.245</td><td>-0.044</td><td></td><td>0.618</td></tr><tr><td>Depression remission</td><td></td><td>1150</td><td><strong>0.011</strong></td><td></td><td>1334</td><td>0.115</td><td></td><td>1588</td><td>0.818</td></tr><tr><td>PANAS positive</td><td>-0.042</td><td></td><td>0.633</td><td>-0.073</td><td></td><td>0.410</td><td>-0.078</td><td></td><td>0.381</td></tr><tr><td>PANAS negative</td><td>0.131</td><td></td><td>0.137</td><td>0.200</td><td></td><td>0.022</td><td>0.146</td><td></td><td>0.098</td></tr><tr><td>WSAS</td><td>0.071</td><td></td><td>0.424</td><td>0.119</td><td></td><td>0.176</td><td>0.099</td><td></td><td>0.262</td></tr><tr><td>Delay between visits</td><td>0.024</td><td></td><td>0.784</td><td>0.072</td><td></td><td>0.417</td><td>-0.059</td><td></td><td>0.509</td></tr></tbody></table>

BMI body mass index, EAT eating attitudes test-26, HADS Hospital Anxiety and Depressive Scale, *p* *p*-value, PANAS positive and negative affect schedule, PC principal component, *r* Pearson's *r*, T1 at baseline, T2 after 4 months of treatment, *U* Mann-Whitney *U*, WSAS Work and Social Adjustment Scale. Bold numbers indicate significant *p*-values after implementation of the Benjamini-Hochberg procedure

interplay between genetics and the environment [32] and differs depending on age at onset; indeed, early- and typical-onset AN show distinct genetic correlation patterns with risk factors for the disease [33].

While the multiple regressions did not highlight these results, we initially observed an association between food avoidance and anxiety, depression and negative affect. Anxiety and mood disorders are common comorbidities in AN [34]. The causal link between food avoidance and mood disorder has not been clearly established: depression may promote food avoidance while food avoidance, e.g. depriving oneself of high-fat foods, may worsen depressive symptoms. Likewise, it is not clear whether anxiety triggers food avoidance or the opposite. On one hand, fastidiously screening the caloric content of foods and avoiding the consumption of calorie-rich foods is a strategy of patients with AN to inefficiently alleviate their anxiety [3]. On the other hand, food avoidance may contribute to anxiety and depression: because of reduced dietary intake, micronutrient status is often altered in patients with AN [11, 35]. Some of the most frequent deficiencies are vitamin B9 and selenium deficits. Both of these elements are essential for neuronal function, and their deficiencies have been linked to depression and anxiety [11, 36]. Such correlations between AN severity, nutritional status, and anxiety and depression, are at the origin of the conceptualization of AN as a metabo-psychiatric disorder [37]. Our additional exploratory analyses using

cross-lagged panel models suggested that anxiety and negative affect caused high-calorie food avoidance, and not the opposite. Depression correlated with high-calorie food avoidance at T2 but the cross-lagged paths were not significant, suggesting that depression is not an actor in the maintenance of food avoidance. This implies that, during treatment, focusing on reducing anxiety and negative affect (but not depression) may be a way to indirectly reduce food avoidance.

Our results also indicated that food avoidance was associated with BMI- and body image-related factors. Correlations with perceived BMI and perceptual distortion make sense since body image distortion is one of the core characteristics of AN [2]. The association of lower minimum lifetime BMI with lower levels of food avoidance was more surprising, as a lower minimum lifetime BMI indicates more severe AN. In an attempt to explain our finding, we examined the difference between current and minimum lifetime BMI, because patients were enrolled in our study at different stages of illness and the difference between current and minimum lifetime BMI can reflect the benefit of care, i.e. weight gain. We observed that a bigger difference between current and minimum lifetime BMI was indeed associated with less food avoidance. Taken together, this suggests that food avoidance could depend not only on illness severity but also on recovery status. In other words, levels of food avoidance were not lower in patients whose minimum BMI
<table><caption>**Table 4** Characteristics of 130 patients with anorexia nervosa who reintroduced at least one high-fat food (versus did not) after 4 months of treatment</caption><thead><tr><th rowspan="2"></th><th colspan="3">Successful food reintroduction (N=39)</th><th colspan="3">No food reintroduction (N=91)</th><th colspan="3">Statistics</th></tr><tr><th>Mean</th><th>SD</th><th>%</th><th>Mean</th><th>SD</th><th>%</th><th>χ<sup>2</sup></th><th>U</th><th>p</th></tr></thead><tbody><tr><td>Age</td><td>24.18</td><td>9.80</td><td></td><td>25.47</td><td>11.37</td><td></td><td></td><td>1604</td><td>0.386</td></tr><tr><td>Education (high)</td><td></td><td>30.8</td><td></td><td></td><td>44.4</td><td>2.11</td><td></td><td></td><td>0.146</td></tr><tr><td>Working (presently)</td><td></td><td>22.6</td><td></td><td></td><td>21.5</td><td>0.01</td><td></td><td></td><td>0.903</td></tr><tr><td>Familial history of ED (yes)</td><td></td><td>23.1</td><td></td><td></td><td>31.5</td><td>0.93</td><td></td><td></td><td>0.336</td></tr><tr><td>Subtype (restrictive)</td><td></td><td>69.2</td><td></td><td></td><td>52.8</td><td>3.00</td><td></td><td></td><td>0.083</td></tr><tr><td>Age at onset</td><td>18.21</td><td>6.72</td><td></td><td>16.78</td><td>3.78</td><td></td><td></td><td>1726</td><td>0.807</td></tr><tr><td>Illness duration</td><td>5.97</td><td>6.54</td><td></td><td>8.69</td><td>10.50</td><td></td><td></td><td>1510</td><td>0.178</td></tr><tr><td>BMI</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>15.48</td><td>1.93</td><td></td><td>15.23</td><td>1.88</td><td></td><td></td><td>1609</td><td>0.402</td></tr><tr><td>Difference</td><td>1.15</td><td>1.67</td><td></td><td>1.07</td><td>1.68</td><td></td><td></td><td>1707</td><td>0.734</td></tr><tr><td>Minimum lifetime BMI</td><td>13.28</td><td>1.58</td><td></td><td>13.40</td><td>1.82</td><td></td><td></td><td>1325</td><td>0.559</td></tr><tr><td>Maximum lifetime BMI</td><td>21.11</td><td>3.62</td><td></td><td>21.29</td><td>3.36</td><td></td><td></td><td>1345</td><td>0.642</td></tr><tr><td>BMI max-current</td><td>5.56</td><td>3.41</td><td></td><td>5.95</td><td>3.12</td><td></td><td></td><td>1274</td><td>0.372</td></tr><tr><td>BMI current-min</td><td>2.28</td><td>1.79</td><td></td><td>1.94</td><td>1.50</td><td></td><td></td><td>1322</td><td>0.548</td></tr><tr><td>BMI max-min</td><td>7.84</td><td>3.62</td><td></td><td>7.89</td><td>3.27</td><td></td><td></td><td>1308</td><td>0.494</td></tr><tr><td>Perceived BMI</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>19.46</td><td>2.82</td><td></td><td>19.71</td><td>3.20</td><td></td><td></td><td>1734</td><td>0.838</td></tr><tr><td>Difference</td><td>1.28</td><td>2.30</td><td></td><td>0.79</td><td>2.51</td><td></td><td></td><td>1571</td><td>0.298</td></tr><tr><td>Perceptual distortion</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>1.27</td><td>0.20</td><td></td><td>1.30</td><td>0.21</td><td></td><td></td><td>1556</td><td>0.268</td></tr><tr><td>Difference</td><td>-0.01</td><td>0.15</td><td></td><td>-0.03</td><td>0.17</td><td></td><td></td><td>1474</td><td>0.127</td></tr><tr><td>Subjective ideal BMI</td><td>17.69</td><td>1.77</td><td></td><td>17.28</td><td>2.06</td><td></td><td></td><td>1341</td><td>0.377</td></tr><tr><td>EAT total</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>30.54</td><td>17.54</td><td></td><td>37.00</td><td>15.73</td><td></td><td></td><td>1383</td><td>0.047</td></tr><tr><td>Difference</td><td>-12.90</td><td>12.58</td><td></td><td>-6.54</td><td>12.85</td><td></td><td></td><td>1288</td><td>0.013</td></tr><tr><td>EAT dieting</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>16.00</td><td>10.85</td><td></td><td>19.75</td><td>10.05</td><td></td><td></td><td>1413</td><td>0.066</td></tr><tr><td>Difference</td><td>-6.54</td><td>7.94</td><td></td><td>-3.41</td><td>7.35</td><td></td><td></td><td>1370</td><td>0.040</td></tr><tr><td>EAT bulimia</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>6.23</td><td>4.39</td><td></td><td>8.51</td><td>4.42</td><td></td><td></td><td>1222</td><td><strong>0.005</strong></td></tr><tr><td>Difference</td><td>-2.46</td><td>3.75</td><td></td><td>-1.36</td><td>3.49</td><td></td><td></td><td>1497</td><td>0.155</td></tr><tr><td>EAT oral</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>8.31</td><td>4.70</td><td></td><td>8.75</td><td>5.00</td><td></td><td></td><td>1693</td><td>0.678</td></tr><tr><td>Difference</td><td>-3.90</td><td>3.62</td><td></td><td>-1.77</td><td>4.71</td><td></td><td></td><td>1215</td><td><strong>0.004</strong></td></tr><tr><td>HADS anxiety</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>11.74</td><td>4.25</td><td></td><td>13.96</td><td>3.92</td><td></td><td></td><td>1243</td><td><strong>0.007</strong></td></tr><tr><td>Remission</td><td></td><td></td><td>25.6</td><td></td><td></td><td>11.0</td><td>4.50</td><td></td><td>0.034</td></tr><tr><td>HADS depression</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>8.49</td><td>4.13</td><td></td><td>9.29</td><td>3.55</td><td></td><td></td><td>1548</td><td>0.249</td></tr><tr><td>Remission</td><td></td><td></td><td>28.2</td><td></td><td></td><td>25.3</td><td>0.12</td><td></td><td>0.728</td></tr><tr><td>PANAS positive</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>29.23</td><td>6.84</td><td></td><td>29.26</td><td>6.87</td><td></td><td></td><td>1709</td><td>0.741</td></tr><tr><td>Difference</td><td>2.87</td><td>5.93</td><td></td><td>0.71</td><td>6.08</td><td></td><td></td><td>1489</td><td>0.147</td></tr><tr><td>PANAS negative</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>33.64</td><td>8.73</td><td></td><td>37.00</td><td>7.34</td><td></td><td></td><td>1360</td><td>0.035</td></tr><tr><td>Difference</td><td>-4.79</td><td>7.42</td><td></td><td>-3.03</td><td>8.15</td><td></td><td></td><td>1542</td><td>0.238</td></tr><tr><td>WSAS</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>At T1</td><td>21.56</td><td>8.26</td><td></td><td>24.23</td><td>8.06</td><td></td><td></td><td>1408</td><td>0.062</td></tr><tr><td>Difference</td><td>-4.85</td><td>9.31</td><td></td><td>-4.01</td><td>10.33</td><td></td><td></td><td>1711</td><td>0.747</td></tr><tr><td>Delay between visits</td><td>126.13</td><td>89.08</td><td></td><td>136.21</td><td>101.77</td><td></td><td></td><td>1719</td><td>0.778</td></tr></tbody></table>

BMI body mass index, EAT eating attitudes test-26, ED eating disorder, HADS Hospital Anxiety and Depressive Scale, *p* *p*-value, PANAS positive and negative affect schedule, SD standard deviation, T1 at baseline, U Mann–Whitney U, WSAS Work and Social Adjustment Scale, χ² Chi-squared test. Bold numbers indicate significant *p*-values after implementation of the Benjamini–Hochberg procedure

<figure><img src="image_6.png" /></figure>
<figure><img src="image_7.png" /><figcaption>**Fig. 2** Cross-lagged panel models depicting the causal relationship between high-calorie food avoidance and A anxiety, B depression, and C negative affect. These models depict the synchronous correlations between the two variables at T1 and at T2, the autoregressive paths of each variable between T1 and T2, and the cross-lagged paths. Numbers indicate regression coefficients. ***$p<0.001$, **$p<0.01$, *$p<0.05$</figcaption></figure>

were less severe, but in those who regained some weight. To confirm this theory, it would have been interesting to have more accurate information about the time trend of patients' illness and weight history.

## Strengths and limits

We hereby present a longitudinal and multicentric study conducted in a sample including both teenagers and adults. This study focuses on a topic which, although central in AN, is understudied, namely food avoidance and more specifically its qualitative dimension. It is understudied to such an extent that no consensual questionnaire assesses food avoidance, hence the need to build one. Our ad hoc questionnaire allowed to detect expected time changes and correlated with illness severity, suggesting that it is a pertinent tool. A potential limitation of this tool is its

subjectivity, especially as insight is impaired in AN [13]. To mitigate this potential bias, it could be interesting to combine this questionnaire to physiological measures such as skin conductance response or pupil size. Another aspect to consider is that our sample was characterized by a wide range of age and illness duration, as it included both adolescents, young adults and older adults (up to 52 years old), while length of the disorder ranged between a couple of months and 34 years. While this variability is a strength (better representativeness), it can also be a limitation. However, in our analyses, we did not find any significant effects of age or illness duration. We performed additional analyses (Table S3) in which we ran our main analyses again, this time in adolescents and adults separately. Associated factors were similar in both groups. Overall, fewer factors were significant in each group than in the whole sample, but this may be due to the reduced statistical power caused by the smaller sample sizes. Other limitations can be considered. Firstly, principal components are not always easy to interpret. Factor scores were computed from the 16 ratings of food avoidance and therefore relied on indirect measures. Nevertheless, component loadings clearly indicated that our three principal components reflected avoidance of high-calorie foods, animal-based foods, and fruits and vegetables, respectively, which corresponds to a relatively simple and intuitive classification. Secondly, the variables included in our regression analyses were not independent. This multicollinearity, although limited (variance inflation factor VIF < 2), may weaken the statistical power of our regression models. Thirdly, a control group could have helped identify factors and associations that are characteristic of AN. Fourthly, some additional measures could have brought some interesting information. Food restriction can be qualitative or quantitative, but our measures only assess self-reported qualitative food restriction, and we do not have data about the nutritional status of patients and the reality of food avoidance. Also, our body image perception test uses drawn silhouettes that are rather minimalistic and may not be as ecologically valid as tests using patients' own silhouettes like in other studies [38], although the present test is less constraining. Finally, it is interesting to note that, as food choices are impacted by the sociocultural context (e.g. cultural values, lifestyles, food movements) [39], our sample exclusively made of French female outpatients may not be representative of all patients with AN.

## What is already known on this subject?

Food avoidance in AN can consist not only in reducing the amount of ingested food, but also in limiting diet diversity. This disordered eating behaviour is responsible for unmet
nutrient needs and is a challenge in the treatment of AN, but research on the topic is scarce.

## What this study adds

* This prospective longitudinal and multicentre study assesses food avoidance in teenagers and adults with AN before and after 4 months of care.

* Even though qualitative food avoidance is not a diagnostic criterion for AN, the present study shows that it could be an informative indicator of AN severity.

* It also suggests that improving anxiety or negative affect may be a leverage to reduce food avoidance.

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1007/s40519-023-01545-4.

**Acknowledgements** The protocol was organized through the FFAB (French Federation Anorexia Bulimia), which is the national network of health care providers working in eating disorders in France. The FFAB network is the national network of health care providers working in eating disorders in France and includes the following as co-authors in the context of this research: Pr Nathalie Godart, Paris, Pr Sébastien Guillaume, Montpellier, Dr. Sylvain Lambert, Nantes, Dr F. Chevallier-Latreuille, Rennes, Dr Brigitte Remy, Paris, Dr Q. Barrois, Dijon, Dr. M. Delorme, Bordeaux, Pr Catherine Massoubre, Saint-Etienne, Pr Vincent Dodin, Lilles, Dr Guillaume Lavoisy, Paris, Dr Sophie Criquillion, Paris, Dr Sylvan Iceta, Lyon, Dr C Fayollet, Paris, Pr Philippe Nubukpo, Limoges, Dr Florat Bat, Marseille.

**Author contributions** Conceptualization: LDL, CV, PG, PD; formal analysis: DP, PG, PD; investigation: LDL, CV, PG, PD; methodology: PG, PD; project administration: PG, PD; writing—original draft: DP, PG, PD; writing—review and editing: LDL, CV, DP, PG, PD. All authors read and approved the final manuscript.

**Funding** There was no funding source for this project.

**Data availability** The data generated and analysed during the current study are available from the corresponding author on request.

## Declarations

**Conflict of interest** Philip Gorwood received, during the last 5 years, fees for presentations at congresses or participation in scientific boards from Alcediag-Alcen, Angelini, GSK, Janssen, Lundbeck, Otsuka, SAGE and Servier. Other authors declare that they have no competing interests.

**Ethical approval** The study protocol was approved by Comité de Protection des Personnes Ile de France III (EUDRACT No: 2008-A008 17-48; CPP NoA: m5355-2-2592). The study was conducted according to ethics recommendations from the Declaration of Helsinki (World Medical Association, 2013).

**Informed consent** All patients gave written informed consent prior to participation. All data were recorded anonymously.

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

## References

1. Arcelus J, Mitchell AJ, Wales J, Nielsen S (2011) Mortality rates in patients with anorexia nervosa and other eating disorders. A meta-analysis of 36 studies. Arch Gen Psychiatry 68:724–731. https://doi.org/10.1001/archgenpsychiatry.2011.74

2. American Psychiatric Association (2013) Diagnostic and Statistical Manual of Mental Disorders, Fifth Edition, DSM-5

3. Murray SB, Strigo IA (2018) Anorexia nervosa, neuroimaging research, and the contextual salience of food cues: the food approach-avoidance conundrum. Int J Eat Disord 51:822–825. https://doi.org/10.1002/eat.22883

4. Stoner SA, Fedoroff IC, Andersen AE, Rolls BJ (1996) Food preferences and desire to eat in anorexia and bulimia nervosa. Int J Eat Disord 19:13–22. https://doi.org/10.1002/(SICI)1098-108X(199601)19:1%3c13::AID-EAT3%3e3.0.CO;2-Z

5. Uniacke B, Slattery R, Walsh BT et al (2020) A comparison of food-based decision-making between restricting and binge-eating/ purging subtypes of anorexia nervosa. Int J Eat Disord 53:1751–1756. https://doi.org/10.1002/eat.23359

6. Schebendach JE, Uniacke B, Walsh BT et al (2019) Fat preference and fat intake in individuals with and without anorexia nervosa. Appetite 139:35–41. https://doi.org/10.1016/j.appet.2019.04.008

7. Werthmann J, Simic M, Konstantellou A et al (2019) Same, same but different: Attention bias for food cues in adults and adolescents with anorexia nervosa. Int J Eat Disord 52:681–690. https://doi.org/10.1002/eat.23064

8. Cowdrey FA, Finlayson G, Park RJ (2013) Liking compared with wanting for high- and low-calorie foods in anorexia nervosa: aberrant food reward even after weight restoration. Am J Clin Nutr 97:463–470. https://doi.org/10.3945/ajcn.112.046011

9. Marzola E, Nasser JA, Hashim SA et al (2013) Nutritional rehabilitation in anorexia nervosa: review of the literature and implications for treatment. BMC Psychiatry 13:290. https://doi.org/10.1186/1471-244X-13-290

10. Hadigan CM, Anderson EJ, Miller KK et al (2000) Assessment of macronutrient and micronutrient intake in women with anorexia nervosa. Int J Eat Disord 28:284–292. https://doi.org/10.1002/1098-108x(200001)28:3%3c284::aid-eat5%3e3.0.co;2-g

11. Achamrah N, Coëffier M, Rimbert A et al (2017) Micronutrient status in 153 patients with anorexia nervosa. Nutrients 9:E225. https://doi.org/10.3390/nu9030225

12. Schebendach JE, Mayer LE, Devlin MJ et al (2008) Dietary energy density and diet variety as predictors of outcome in anorexia nervosa. Am J Clin Nutr 87:810–816. https://doi.org/10.1093/ajcn/87.4.810

13. Gorwood P, Duriez P, Lengvenyte A et al (2019) Clinical insight in anorexia nervosa: associated and predictive factors. Psychiatry Res 281:112561. https://doi.org/10.1016/j.psychres.2019.112561
14. Duriez P, Kaya Lefèvre H, Di Lodovico L et al (2021) Increased cognitive flexibility mediates the improvement of eating disorders symptoms, depressive symptoms and level of daily life functioning in patients with anorexia nervosa treated in specialised centres. *Eur Eat Disord Rev* 29:600–610. https://doi.org/10.1002/erv.2829

15. Garner DM, Garfinkel PE (1979) The Eating Attitudes Test: an index of the symptoms of anorexia nervosa. *Psychol Med* 9:273–279. https://doi.org/10.1017/s0033291700030762

16. Leichner P, Steiger H, Puentes-Neuman G et al (1994) Validation of an eating attitude scale in a French-speaking Quebec population. *Can J Psychiatry* 39:49–54. https://doi.org/10.1177/070674379403900110

17. Zigmond AS, Snaith RP (1983) The hospital anxiety and depression scale. *Acta Psychiatr Scand* 67:361–370. https://doi.org/10.1111/j.1600-0447.1983.tb09716.x

18. Bjelland I, Dahl AA, Haug TT, Neckelmann D (2002) The validity of the Hospital Anxiety and Depression Scale. An updated literature review. *J Psychosom Res* 52:69–77. https://doi.org/10.1016/s0022-3999(01)00296-3

19. Crawford JR, Henry JD (2004) The positive and negative affect schedule (PANAS): construct validity, measurement properties and normative data in a large non-clinical sample. *Br J Clin Psychol* 43:245–265. https://doi.org/10.1348/0144665031752934

20. Mundt JC, Marks IM, Shear MK, Greist JH (2002) The Work and Social Adjustment Scale: a simple measure of impairment in functioning. *Br J Psychiatry* 180:461–464. https://doi.org/10.1192/bjp.180.5.461

21. Williamson DA, Cubic BA, Gleaves DH (1993) Equivalence of body image disturbances in anorexia and bulimia nervosa. *J Abnorm Psychol* 102:177–180

22. (2021) The Jamovi Project. Jamovi

23. Benjamini Y, Hochberg Y (1995) Controlling the false discovery rate: a practical and powerful approach to multiple testing. *J Roy Stat Soc: Ser B (Methodol)* 57:289–300. https://doi.org/10.1111/j.2517-6161.1995.tb02031.x

24. Lee S, Lee DK (2018) What is the proper way to apply the multiple comparison test? *Korean J Anesthesiol* 71:353–360. https://doi.org/10.4097/kja.d.18.00242

25. Cattell RB (1966) The scree test for the number of factors. Multivariate Behav Res 1:245–276. https://doi.org/10.1207/s15327906mbr0102_10

26. Horn JL (1965) A rationale and test for the number of factors in factor analysis. *Psychometrika* 30:179–185. https://doi.org/10.1007/BF02289447

27. Franklin SB, Gibson DJ, Robertson PA et al (1995) Parallel analysis: a method for determining significant principal components. *J Veg Sci* 6:99–106. https://doi.org/10.2307/3236261

28. Treasure J, Duarte TA, Schmidt U (2020) Eating disorders. *Lancet* 395:899–911. https://doi.org/10.1016/S0140-6736(20)30059-3

29. Coniglio KA, Becker KR, Franko DL et al (2017) Won't stop or can't stop? Food restriction as a habitual behavior among individuals with anorexia nervosa or atypical anorexia nervosa. *Eur Eat Disord Rev* 26:144–147. https://doi.org/10.1016/j.eatbeh.2017.03.005

30. Grilo CM, Udo T (2021) Examining the significance of age of onset in persons with lifetime anorexia nervosa: comparing child, adolescent, and emerging adult onsets in nationally representative U.S. study. *Int J Eat Disord* 54:1632–1640. https://doi.org/10.1002/eat.23580

31. Jenkins ZM, Chait LM, Cistullo L, Castle DJ (2020) A comparison of eating disorder symptomatology, psychological distress and psychosocial function between early, typical and later onset anorexia nervosa. *J Eat Disord* 8:56. https://doi.org/10.1186/s40337-020-00337-w

32. Zipfel S, Giel KE, Bulik CM et al (2015) Anorexia nervosa: aetiology, assessment, and treatment. *Lancet Psychiatry* 2:1099–1111. https://doi.org/10.1016/S2215-0366(15)00356-9

33. Watson HJ, Thornton LM, Yilmaz Z et al (2021) Common genetic variation and age at onset of anorexia nervosa. *Biol Psychiatry* Global Open Sci. https://doi.org/10.1016/j.bpsgos.2021.09.001

34. Marucci S, Ragione LD, De Iaco G et al (2018) Anorexia nervosa and comorbid psychopathology. *Endocr Metab Immune Disord Drug Targets* 18:316–324. https://doi.org/10.2174/1871530318666180213111637

35. Hanachi M, Dicembre M, Rives-Lange C et al (2019) Micronutrients deficiencies in 374 severely malnourished anorexia nervosa inpatients. *Nutrients* 11:E792. https://doi.org/10.3390/nu11040792

36. Mikkelsen K, Stojanovska L, Apostolopoulos V (2016) The effects of vitamin B in depression. *Curr Med Chem* 23:4317–4337. https://doi.org/10.2174/0929867323666160920110810

37. Watson HJ, Yilmaz Z, Thornton LM et al (2019) Genome-wide association study identifies eight risk loci and implicates metabopychiatric origins for anorexia nervosa. *Nat Genet* 51:1207–1214. https://doi.org/10.1038/s41588-019-0439-2

38. Couton C, Gorwood P, Pham-Scottez A et al (2022) Pupil psychosensory reflex in response to own and standardised silhouettes in patients with anorexia nervosa. *Eur Eat Disord Rev.* https://doi.org/10.1002/erv.2881

39. Monterrosa EC, Frongillo EA, Drewnowski A et al (2020) Socio-cultural influences on food choices and implications for sustainable healthy diets. *Food Nutr Bull* 41:59S-73S. https://doi.org/10.1177/0379572120975874

**Publisher's Note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.