Parental sociodemographic data. The 55 fathers had a mean age of 47.41 years (SD = 4.67) [37;60] and the 80 mothers had a mean age of 45.81 years (SD = 5.49) [35;61]. Professional and socio-economic status [38], 40.40% of fathers were in the upper to high-middle class category, and 32.10% belonged to the middle-class category. Regarding mothers, 48.80% belonged to the middle-class category and 26.20% were in the lower-middle class category.
Parental levels of emotional distress, grief and burden (in Supplementary Table 1). Mothers showed significantly higher scores of emotional distress than fathers for anxiety and depression, while fathers indicated significantly higher scores of alexithymia than mothers (p < 0.01). There were no significant differences between maternal and paternal scores of grief (MIV-Grief) and burden (ECI- burden).
Adolescents’ sociodemographic and clinical data. The majority of them were girls (96.4%, 81/84) with a mean age of 15.98 years (SD = 2.32). 25% were enrolled in primary education, 35.7% in secondary education (30/84), 7.2% in higher education (missing data for 31%). A minority was adopted (3) and 28.6% had parents who were separated. 46.40% had one sibling, 31% had two, 9.5% had three and 13.1% were only children.
The sample had a mean age of AN onset of 14.61 years (SD = 2.38) [8;23], a mean illness duration of 1.84 years (SD = 1,67) [0,2;7,9], a BMI of 14.3 upon admission (SD = 1,39) [11,1;18,9] and a minimum lifetime BMI of 13.61 (SD = 1,55) [9,9;18,5]. 53.6% met criteria for restrictive type AN (AN-R). The MR-GOAS mean score was 4.90 (SD = 1.41) [2,04; 7,89], suggesting poor clinical state.
Burden (Supplementary Table 2): For both mothers and fathers, their level of emotional distress and their child’s clinical state were both linked with their level of burden. For mothers, higher levels of anxiety (HADS-anxiety score; p = 0.03) and of depression (BDI score; p = 0.01), a longer duration of their child’s illness (p = 0.04) and a worse child clinical state (MR-GOAS score, p<0.0001) were all significantly associated with a higher maternal burden score (ECI – Burden score). The maternal level of alexithymia (BVAQ score) tended to be associated with a higher level of maternal burden (p = 0.07).
For fathers, a higher level of anxiety (HADS-anxiety score; p = 0.001) and a worse child clinical state (MR-GOAS score; tendency p = 0.06) were significantly related to a higher level of paternal burden (ECI score).
Grief (Supplementary Table 3): For both mothers and fathers, their level of emotional distress and their child’s clinical state were related to their level of grief. For mothers, higher levels of anxiety (HADS- anxiety score; p = 0.005), depression (BDI score; p<0.0001) and alexithymia (BVAQ score; p = 0.001), and a worse clinical state for their child (MR-GOAS score; p = 0.004)) were all significantly associated with a higher current grief score (lower MIV-Grief score). For fathers, higher levels of anxiety (HADS- anxiety score; p = 0.002) and depression (BDI score; p<0.003), and a worse clinical state for their child (MR-GOAS score; p = 0.03) were significantly associated with a higher current grief score (lower MIV-Grief score).
For mothers’ Model 1: on the basis of the univariate analysis, the following elements were introduced as explanatory variables: anxiety, depression, alexithymia, duration of illness and child’s clinical state. Higher maternal burden was related to a worse child clinical state (MR-GOAS score, β = -0.35, p = 0.00). However, it was not associated with mothers’ anxiety (HADS-anxiety score, β = 0.00, p = 0.96), depression (BDI score, β = 0.22, p = 0.21) or alexithymia (BVAQ score, β = 0.97, p = 0.40), nor with child AN duration (β = 0.15, p = 0.15). This model explained 18% of the variance (R² = 0.18).
For fathers’ Model 1: We also performed a paternal model with the following explanatory variables: anxiety, depression and child clinical state. A higher paternal level of burden was significantly associated with a higher level of fathers’ anxiety (HADS-anxiety score, β = 0.48, p = 0.00) and worse child clinical state (MR-GOAS score, β = -0.33, p = 0.00). Fathers’ level of burden was not associated with their level of depression (BDI score, β = 0.04, p = 0.71). This model explained 27% of the variance (R² = 0.27).
For mothers’ Model 2: the following elements were introduced as explanatory variables, on the basis of the univariate analysis: anxiety, depression, alexithymia, and worse child’s clinical state. Higher grief was explained by higher levels of mothers’ depression (BDI score, β = -0.55, p = 0.00) and alexithymia (BVAQ score, β = -0.24, p = 0.01), and by children’s worse clinical state (MR-GOAS score, β = 0.36, p = 0.00), but not by mothers’ anxiety levels (HADS-anxiety score, β = 0.24, p = 0.11). This model explained 35% of the variance (R² = 0.35).
For fathers’ Model 2: a similar paternal model was performed for grief, and the following explanatory variables were introduced: anxiety, depression, and worse child clinical state.
A higher paternal level of grief was explained by higher anxiety (HADS-anxiety score, β = -0.36, p = 0.001) and depression (BDI score, β = -0.25, p = 0.04) levels, and children’s worse clinical state (MR-GOAS score, β = 0.33, p = 0.001). This model explained 31% of the variance (R² = 0.31).
A third model was performed in order to test the relationship between burden and grief for mothers and then fathers, taking into account potential contributions identified in models 1 and 2.
For mothers’ Model 3: higher mothers’ burden levels were explained by higher grief levels (ECI score, β = -0.31, p = 0.01) and worse child clinical state (MR-GOAS score, β = -0.24, p = 0.03).
The other variables were not significant: mothers’ anxiety (HADS-anxiety score, β = 0.10, p = 0.56), depression (BDI score, β = 0.03, p = 0.85) and alexithymia levels (BVAQ score, β = 0.02, p = 0.85), as well as children’s illness duration (β = 0.15, p = 0.14) were not linked with mothers’ levels of burden. This model explained 24% of the variance (R² = 0.24).
For Fathers’ Model 3: Fathers’ higher burden levels were explained by higher grief levels (ECI score, β = -0.35, p = 0.01) and higher levels of anxiety (HADS-anxiety, β = 0.34, p = 0.01) but not by the other variables which were not significant: depression levels (BDI score, β = -0.04, p = 0.71) and child clinical state (MR-GOAS score, β = -0.21, p = 0.01). This model explained 34% of the variance (R² = 0.34).